Adolescent idiopathic scoliosis (AIS) is a lateral spinal curvature >10° with rotation that affects 2–3% of healthy children across populations. AIS is known to have a significant genetic component, and despite a handful of risk loci identified in unrelated individuals by GWAS and next-generation sequencing methods, the underlying etiology of the condition remains largely unknown. In this study, we performed exome sequencing of affected individuals within 23 multigenerational families, with the hypothesis that the occurrence of rare, low frequency, disease-causing variants will co-occur in distantly related, affected individuals. Bioinformatic filtering of uncommon, potentially damaging variants shared by all sequenced family members revealed 1448 variants in 1160 genes across the 23 families, with 132 genes shared by two or more families. Ten genes were shared by >4 families, and no genes were shared by all. Gene enrichment analysis showed an enrichment of variants in cytoskeletal and extracellular matrix related processes. These data support a model that AIS is a highly polygenic disease, with few variant-containing genes shared between affected individuals across different family lineages. This work presents a novel resource for further exploration in familial AIS genetic research.

Idiopathic scoliosis (IS) is a three-dimensional rotation of the spine >10 degrees with an unknown etiology. Our laboratory established a late-onset IS model in zebrafish (Danio rerio) containing a deletion in kif7. A total of 25% of kif7^co63/co63 zebrafish develop spinal curvatures and are otherwise developmentally normal, although the molecular mechanisms underlying the scoliosis are unknown. To define transcripts associated with scoliosis in this model, we performed bulk mRNA sequencing on 6 weeks past fertilization (wpf) kif7^co63/co63 zebrafish with and without scoliosis. Additionally, we sequenced kif7^co63/co63, kif7^co63/+, and AB zebrafish (n = 3 per genotype). Sequencing reads were aligned to the GRCz11 genome and FPKM values were calculated. Differences between groups were calculated for each transcript by the t-test. Principal component analysis showed that transcriptomes clustered by sample age and genotype. kif7 mRNA was mildly reduced in both homozygous and heterozygous zebrafish compared to AB. Sonic hedgehog target genes were upregulated in kif7^co63/co63 zebrafish over AB, but no difference was detected between scoliotic and non-scoliotic mutants. The top upregulated genes in scoliotic zebrafish were cytoskeletal keratins. Pankeratin staining of 6 wpf scoliotic and non-scoliotic kif7^co63/co63 zebrafish showed increased keratin levels within the zebrafish musculature and intervertebral disc (IVD). Keratins are major components of the embryonic notochord, and aberrant keratin expression has been associated with intervertebral disc degeneration (IVDD) in both zebrafish and humans. The role of increased keratin accumulation as a molecular mechanism associated with the onset of scoliosis warrants further study.

It is common practice for laboratories to discard clotted blood or freeze it for future DNA extraction after extracting serum from a serum-separating tube. If freezing for DNA extraction, the blood clot is not usually cryopreserved, which leads to cell membrane fragility. In this protocol, we describe steps to isolate high-quality nuclei from leukocytes derived from whole blood samples frozen without a cryoprotective medium. Nuclei isolated from this protocol were able to undergo ATAC (assay for transposase-accessible chromatin) sequencing to obtain chromatin accessibility data. We successfully characterized and isolated B cells and T cells from leukocytes isolated from previously frozen blood clot using Miltenyi’s gentleMACS Octo Dissociator coupled with flow sorting. Nuclei showed round, intact nuclear envelopes suitable for downstream applications, including bulk sequencing of nuclei or single-cell nuclei sequencing. We validated this protocol by performing bulk ATAC-seq.