Association between genetic polymorphisms and risk of adolescent idiopathic scoliosis in case-control studies: a systematic review.

BACKGROUND: Adolescent idiopathic scoliosis (AIS) is a structural lateral spinal curvature of ≥10° with rotation. Approximately 2%-3% of children across populations are affected with AIS, and this condition is responsible for ~$3 billion in costs within the USA. Although AIS is believed to have a strong genetic contribution, clinical translation of identified genetic variants has stalled. METHODS: The databases MEDLINE (via PubMed), Embase, Google Scholar and Ovid MEDLINE were searched and limited to articles in English. Title and abstract, full-text and data extraction screening was conducted through Covidence, followed by data transfer to a custom REDCap database. Studies containing variant-level data using genome-wide methodology as well as validation studies of genome-wide methods were considered. Quality assessment was conducted using Q-Genie. RESULTS: 33 studies were included, including 9 genome-wide association studies, 4 whole exome sequencing and 20 validation studies. Combined, these studies included data from >35,000 cases and >67,000 controls, not including validation cohorts. Additionally, results from six meta-analyses containing novel cohorts were also reported. All included study cohorts were from populations of primarily East Asian or Caucasian descent. Quality assessment found that overall study quality was high and control group selection was moderate. The highest number of reported associations were in single nucleotide polymorphisms (SNPs) in or near LBX1, LBX1-AS1, GPR126/ADGRG6 or BNC2. CONCLUSION: AIS risk may be influenced by specific SNPs, particularly those in/near LBX1 and GPR126. Translatability of study findings is unknown due to an underrepresentation of most ethnic groups as well as few identified genome-wide studies. Further studies may benefit from increased cohort diversity and thorough evaluation of control cohort groups.

Mutations in KIF7 implicated in idiopathic scoliosis in humans and axial curvatures in zebrafish.

Idiopathic scoliosis (IS) is a spinal disorder affecting up to 3% of otherwise healthy children. IS has a strong familial genetic component and is believed to be genetically complex due to significant variability in phenotype and heritability. Previous studies identified putative loci and variants possibly contributing to IS susceptibility, including within extracellular matrix, cilia, and actin networks, but the genetic architecture and underlying mechanisms remain unresolved. Here, we used whole-exome sequencing from three affected individuals in a multigenerational family with IS and identified 19 uncommon variants (minor allele frequency < 0.05). Genotyping of additional family members identified a candidate heterozygous variant (H1115Q, G>C, rs142032413) within the ciliary gene KIF7, a regulator within the hedgehog (Hh) signaling pathway. Resequencing of the second cohort of unrelated IS individuals and controls identified several severe mutations in KIF7 in affected individuals only. Subsequently, we generated a mutant zebrafish model of kif7 using CRISPR-Cas9. kif7co63/co63 zebrafish displayed severe scoliosis, presenting in juveniles and progressing through adulthood. We observed no deformities in the brain, Reissner fiber, or central canal cilia in kif7co63/co63 embryos, although alterations were seen in Hh pathway gene expression. This study suggests defects in KIF7-dependent Hh signaling, which may drive pathogenesis in a subset of individuals with IS.

Cytoskeletal Keratins Are Overexpressed in a Zebrafish Model of Idiopathic Scoliosis.

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 kif7co63/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) kif7co63/co63 zebrafish with and without scoliosis. Additionally, we sequenced kif7co63/co63, kif7co63/+, 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 kif7co63/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 kif7co63/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.

Rare Coding Variants in Patients with Non-Syndromic Vestibular Dysfunction.

Vertigo due to vestibular dysfunction is rare in children. The elucidation of its etiology will improve clinical management and the quality of life of patients. Genes for vestibular dysfunction were previously identified in patients with both hearing loss and vertigo. This study aimed to identify rare, coding variants in children with peripheral vertigo but no hearing loss, and in patients with potentially overlapping phenotypes, namely, Meniere's disease or idiopathic scoliosis. Rare variants were selected from the exome sequence data of 5 American children with vertigo, 226 Spanish patients with Meniere's disease, and 38 European-American probands with scoliosis. In children with vertigo, 17 variants were found in 15 genes involved in migraine, musculoskeletal phenotypes, and vestibular development. Three genes, OTOP1, HMX3, and LAMA2, have knockout mouse models for vestibular dysfunction. Moreover, HMX3 and LAMA2 were expressed in human vestibular tissues. Rare variants within ECM1, OTOP1, and OTOP2 were each identified in three adult patients with Meniere's disease. Additionally, an OTOP1 variant was identified in 11 adolescents with lateral semicircular canal asymmetry, 10 of whom have scoliosis. We hypothesize that peripheral vestibular dysfunction in children may be due to multiple rare variants within genes that are involved in the inner ear structure, migraine, and musculoskeletal disease.

Genetic animal modeling for idiopathic scoliosis research: history and considerations.

PURPOSE: Idiopathic scoliosis (IS) is defined as a structural lateral spinal curvature ≥ 10° in otherwise healthy children and is the most common pediatric spinal deformity. IS is known to have a strong genetic component; however, the underlying etiology is still largely unknown. Animal models have been used historically to both understand and develop treatments for human disease, including within the context of IS. This intended audience for this review is clinicians in the fields of musculoskeletal surgery and research. METHODS: In this review article, we synthesize current literature of genetic animal models of IS and introduce considerations for researchers. RESULTS: Due to complex genetic and unique biomechanical factors (i.e., bipedalism) hypothesized to contribute to IS in humans, scoliosis is a difficult condition to replicate in model organisms. CONCLUSION: We advocate careful selection of animal models based on the scientific question and introduce gaps and limitations in the current literature. We advocate future research efforts to include animal models with multiple characterized genetic or environmental perturbations to reflect current understanding of the human condition.

Whole Exome Sequencing of 23 Multigeneration Idiopathic Scoliosis Families Reveals Enrichments in Cytoskeletal Variants, Suggests Highly Polygenic Disease.

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.

Severity of Idiopathic Scoliosis Is Associated with Differential Methylation: An Epigenome-Wide Association Study of Monozygotic Twins with Idiopathic Scoliosis.

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Idiopathic Scoliosis Families Highlight Actin-Based and Microtubule-Based Cellular Projections and Extracellular Matrix in Disease Etiology.

Idiopathic scoliosis (IS) is a structural lateral spinal curvature of ≥10° that affects up to 3% of otherwise healthy children and can lead to life-long problems in severe cases. It is well-established that IS is a genetic disorder. Previous studies have identified genes that may contribute to the IS phenotype, but the overall genetic etiology of IS is not well understood. We used exome sequencing to study five multigenerational families with IS. Bioinformatic analyses identified unique and low frequency variants (minor allele frequency ≤5%) that were present in all sequenced members of the family. Across the five families, we identified a total of 270 variants with predicted functional consequences in 246 genes, and found that eight genes were shared by two families. We performed GO term enrichment analyses, with the hypothesis that certain functional annotations or pathways would be enriched in the 246 genes identified in our IS families. Using three complementary programs to complete these analyses, we identified enriched categories that include stereocilia and other actin-based cellular projections, cilia and other microtubule-based cellular projections, and the extracellular matrix (ECM). Our results suggest that there are multiple paths to IS and provide a foundation for future studies of IS pathogenesis.

Sequencing of the TBX6 Gene in Families with Familial Idiopathic Scoliosis.

STUDY DESIGN: A hypothesis-driven study was conducted in a familial cohort to determine the potential association between variants within the TBX6 gene and Familial Idiopathic Scoliosis (FIS). OBJECTIVE: To determine if variants within exons of the TBX6 gene segregate with the FIS phenotype within a sample of families with FIS. SUMMARY OF BACKGROUND DATA: Idiopathic Scoliosis (IS) is a structural curvature of the spine whose underlying genetic etiology has not been established. IS has been reported to occur at a higher rate than expected in family members of individuals with congenital scoliosis (CS), suggesting that the two diseases might have a shared etiology. The TBX6 gene on chromosome 16p, essential to somite development, has been associated with CS in a Chinese population. Previous studies have identified linkage to this locus in families with FIS, and specifically with rs8060511, located in an intron of the TBX6 gene. METHODS: Parent-offspring trios from 11 families (13 trios, 42 individuals) with FIS were selected for Sanger sequencing of the TBX6 gene. Trios were selected from a large population of families with FIS in which a genome-wide scan had resulted in linkage to 16p. RESULTS: Sequencing analyses of the subset of families resulted in the identification of five coding variants. Three of the five variants were novel; the remaining two variants were previously characterized and account for 90% of the observed variants in these trios. In all cases, there was no correlation between transmission of the TBX6 variant allele and FIS phenotype. However, an analysis of regulatory markers in osteoblasts showed that rs8060511 is in a putative enhancer element. CONCLUSIONS: Although this study did not identify any TBX6 coding variants that segregate with FIS, we identified a variant that is located in a potential TBX6 enhancer element. Therefore, further investigation of the region is needed.

Exome sequencing identifies a rare HSPG2 variant associated with familial idiopathic scoliosis.

Idiopathic scoliosis occurs in 3% of individuals and has an unknown etiology. The objective of this study was to identify rare variants that contribute to the etiology of idiopathic scoliosis by using exome sequencing in a multigenerational family with idiopathic scoliosis. Exome sequencing was completed for three members of this multigenerational family with idiopathic scoliosis, resulting in the identification of a variant in the HSPG2 gene as a potential contributor to the phenotype. The HSPG2 gene was sequenced in a separate cohort of 100 unrelated individuals affected with idiopathic scoliosis and also was examined in an independent idiopathic scoliosis population. The exome sequencing and subsequent bioinformatics filtering resulted in 16 potentially damaging and rare coding variants. One of these variants, p.Asn786Ser, is located in the HSPG2 gene. The variant p.Asn786Ser also is overrepresented in a larger cohort of idiopathic scoliosis cases compared with a control population (P = 0.024). Furthermore, we identified additional rare HSPG2 variants that are predicted to be damaging in two independent cohorts of individuals with idiopathic scoliosis. The HSPG2 gene encodes for a ubiquitous multifunctional protein within the extracellular matrix in which loss of function mutation are known to result in a musculoskeletal phenotype in both mouse and humans. Based on these results, we conclude that rare variants in the HSPG2 gene potentially contribute to the idiopathic scoliosis phenotype in a subset of patients with idiopathic scoliosis. Further studies must be completed to confirm the effect of the HSPG2 gene on the idiopathic scoliosis phenotype.