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.

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.

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.

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.