Search

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near thePAX1gene. Here we sought to define the roles ofPAX1and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 9,161 individuals with AIS and 80,731 unaffected controls, significant association was identified with a variant inCOL11A1encoding collagen (α1) XI (rs3753841; NM_080629_c.4004C>T; p.(Pro1335Leu); P=7.07e−11, OR=1.118). Using CRISPR mutagenesis we generatedPax1knockout mice (Pax1−/−). In postnatal spines we found that Pax1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected inPax1−/−spines compared to wildtype. By genetic targeting we found that wildtypeCol11a1expression in growth plate cells (GPCs) suppresses expression ofPax1and ofMmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associatedCOL11A1P1335Lmutant. Further, we found that either knockdown of the estrogen receptor geneEsr2, or tamoxifen treatment, significantly alteredCol11a1andMmp3expression in GPCs. These studies support a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering aPax1-Col11a1-Mmp3signaling axis in the growth plate.

BACKGROUND Clubfoot, a congenital deformity that presents as a rigid, inward turning of the foot, affects approximately 1 in 1000 infants and occurs as an isolated birth defect in 80% of patients. Despite its high level of heritability, few causative genes have been identified, and mutations in known genes are only responsible for a small portion of clubfoot heritability. QUESTIONS/PURPOSES (1) Are any rare gene variants enriched (that is, shared) in unrelated patients with isolated clubfoot? (2) Are there other rare variants in the identified gene (Filamin B) in these patients with clubfoot? METHODS Whole-exome sequence data were generated from a discovery cohort of 183 unrelated probands with clubfoot and 2492 controls. Variants were filtered with minor allele frequency < 0.02 to identify rare variants as well as small insertions and deletions (indels) resulting in missense variants, nonsense or premature truncation, or in-frame deletions. A candidate deletion was then genotyped in another cohort of 974 unrelated patients with clubfoot (a replication cohort). Other rare variants in the candidate gene were also investigated. A segregation analysis was performed in multigenerational families of individuals with clubfoot to see if the genotypes segregate with phenotypes. Single-variant association analysis was performed using the Fisher two-tailed exact test (exact p values are presented to give an indication of the magnitude of the association). RESULTS There were no recurrent variants in the known genes causing clubfoot in this study. A three-base pair in-frame codon deletion of Filamin B (FLNB) (p.E1792del, rs1470699812) was identified in 1.6% (3 of 183) of probands with clubfoot in the discovery cohort compared with 0% of controls (0 of 2492) (odds ratio infinity (inf) [95% CI 5.64 to inf]; p = 3.18 x 10-5) and 0.0016% of gnomAD controls (2 of 125,709) (OR 1.01 x 103 [95% CI 117.42 to 1.64 x 104]; p = 3.13 x 10-8). By screening a replication cohort (n = 974 patients), we found two probands with the identical FLNB deletion. In total, the deletion was identified in 0.43% (5 of 1157) of probands with clubfoot compared with 0% of controls and 0.0016% of gnomAD controls (OR 268.5 [95% CI 43.68 to 2.88 x 103]; p = 1.43 x 10-9). The recurrent FLNB p.E1792del variant segregated with clubfoot, with incomplete penetrance in two families. Affected individuals were more likely to be male and have bilateral clubfoot. Although most patients had isolated clubfoot, features consistent with Larsen syndrome, including upper extremity abnormalities such as elbow and thumb hypermobility and wide, flat thumbs, were noted in affected members of one family. We identified 19 additional rare FLNB missense variants located throughout the gene in patients with clubfoot. One of these missense variants, FLNB p.G2397D, exhibited incomplete penetrance in one family. CONCLUSION A recurrent FLNB E1792 deletion was identified in 0.43% of 1157 isolated patients with clubfoot. Given the absence of any recurrent variants in our discovery phase (n = 183) for any of the known genes causing clubfoot, our findings support that novel and rare missense variants in FLNB in patients with clubfoot, although rare, may be among the most commonly known genetic causes of clubfoot. Patients with FLNB variants often have isolated clubfoot, but they and their family members may be at an increased risk of having additional clinical features consistent with Larsen syndrome. CLINICAL RELEVANCE Identification of FLNB variants may be useful for determining clubfoot recurrence risk and comorbidities.

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 ECM, cilia and actin networks, but the genetic architecture and underlying mechanisms remains unresolved. Here, we used whole exome sequencing from three affected individuals in a multigenerational family with IS and identified 19 uncommon variants (MAF C, rs142032413) within the ciliary gene KIF7, a regulator within the hedgehog (Hh) signaling pathway. Resequencing of a 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. kif7(co63/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 kif7(co63/co63) embryos, although alterations were seen in Hh pathway gene expression. This research suggests defects in KIF7-dependent Hh signaling may drive pathogenesis in a subset of individuals with IS.

IntroductionCongenital clubfoot is a common birth defect that affects at least 0.1% of all births. Nearly 25% cases are familial and the remaining are sporadic in inheritance. Copy number variants (CNVs) involving transcriptional regulators of limb development, including PITX1 and TBX4, have previously been shown to cause familial clubfoot, but much of the heritability remains unexplained.MethodsExome sequence data from 816 unrelated clubfoot cases and 2645 in-house controls were analysed using coverage data to identify rare CNVs. The precise size and location of duplications were then determined using high-density Affymetrix Cytoscan chromosomal microarray (CMA). Segregation in families and de novo status were determined using qantitative PCR.ResultsChromosome Xp22.33 duplications involving SHOX were identified in 1.1% of cases (9/816) compared with 0.07% of in-house controls (2/2645) (p=7.98×10−5, OR=14.57) and 0.27% (38/13592) of Atherosclerosis Risk in Communities/the Wellcome Trust Case Control Consortium 2 controls (p=0.001, OR=3.97). CMA validation confirmed an overlapping 180.28 kb duplicated region that included SHOX exons as well as downstream non-coding regions. In four of six sporadic cases where DNA was available for unaffected parents, the duplication was de novo. The probability of four de novo mutations in SHOX by chance in a cohort of 450 sporadic clubfoot cases is 5.4×10–10.ConclusionsMicroduplications of the pseudoautosomal chromosome Xp22.33 region (PAR1) containing SHOX and downstream enhancer elements occur in ~1% of patients with clubfoot. SHOX and regulatory regions have previously been implicated in skeletal dysplasia as well as idiopathic short stature, but have not yet been reported in clubfoot. SHOX duplications likely contribute to clubfoot pathogenesis by altering early limb development.

IntroductionAdolescent idiopathic scoliosis (AIS) is a common musculoskeletal disorder with strong evidence for a genetic contribution. CNVs play an important role in congenital scoliosis, but their role in idiopathic scoliosis has been largely unexplored.MethodsExome sequence data from 1197 AIS cases and 1664 in-house controls was analysed using coverage data to identify rare CNVs. CNV calls were filtered to include only highly confident CNVs with >10 average reads per region and mean log-ratio of coverage consistent with single-copy duplication or deletion. The frequency of 55 common recurrent CNVs was determined and correlated with clinical characteristics.ResultsDistal chromosome 16p11.2 microduplications containing the gene SH2B1 were found in 0.7% of AIS cases (8/1197). We replicated this finding in two additional AIS cohorts (8/1097 and 2/433), resulting in 0.7% (18/2727) of all AIS cases harbouring a chromosome 16p11.2 microduplication, compared with 0.06% of local controls (1/1664) and 0.04% of published controls (8/19584) (p=2.28×10−11, OR=16.15). Furthermore, examination of electronic health records of 92 455 patients from the Geisinger health system showed scoliosis in 30% (20/66) patients with chromosome 16p11.2 microduplications containing SH2B1 compared with 7.6% (10/132) of controls (p=5.6×10−4, OR=3.9).ConclusionsRecurrent distal chromosome 16p11.2 duplications explain nearly 1% of AIS. Distal chromosome 16p11.2 duplications may contribute to scoliosis pathogenesis by directly impairing growth or by altering expression of nearby genes, such as TBX6. Individuals with distal chromosome 16p11.2 microduplications should be screened for scoliosis to facilitate early treatment.

Summary Chiari I malformation (CM1), the displacement of the cerebellum through the foramen magnum into the spinal canal, is one of the most common pediatric neurological conditions. Individuals with CM1 can present with neurological symptoms, including severe headaches and sensory or motor deficits, often as a consequence of brainstem compression or syringomyelia (SM). We conducted whole-exome sequencing (WES) on 668 CM1 probands and 232 family members and performed gene-burden and de novo enrichment analyses. A significant enrichment of rare and de novo non-synonymous variants in chromodomain (CHD) genes was observed among individuals with CM1 (combined p = 2.4 × 10−10), including 3 de novo loss-of-function variants in CHD8 (LOF enrichment p = 1.9 × 10−10) and a significant burden of rare transmitted variants in CHD3 (p = 1.8 × 10−6). Overall, individuals with CM1 were found to have significantly increased head circumference (p = 2.6 × 10−9), with many harboring CHD rare variants having macrocephaly. Finally, haploinsufficiency for chd8 in zebrafish led to macrocephaly and posterior hindbrain displacement reminiscent of CM1. These results implicate chromodomain genes and excessive brain growth in CM1 pathogenesis.

Distal arthrogryposis (DA) is group of syndromes characterized by congenital joint contractures. Treatment development is hindered by the lack of vertebrate models. Here, we describe a zebrafish model in which a common MYH3 missense mutation (R672H) was introduced into the orthologous zebrafish gene smyhc1 (slow myosin heavy chain 1) (R673H). We simultaneously created a smyhc1 null allele (smyhc1 −), which allowed us to compare the effects of both mutant alleles on muscle and bone development, and model the closely related disorder, spondylocarpotarsal synostosis syndrome. Heterozygous smyhc1 R673H/+ embryos developed notochord kinks that progressed to scoliosis with vertebral fusions; motor deficits accompanied the disorganized and shortened slow‐twitch skeletal muscle myofibers. Increased dosage of the mutant allele in both homozygous smyhc1 R673H/R673H and transheterozygous smyhc1 R673H/− embryos exacerbated the notochord and muscle abnormalities, causing early lethality. Treatment of smyhc1 R673H/R673H embryos with the myosin ATPase inhibitor, para‐aminoblebbistatin, which decreases actin–myosin affinity, normalized the notochord phenotype. Our zebrafish model of MYH3‐associated DA2A provides insight into pathogenic mechanisms and suggests a beneficial therapeutic role for myosin inhibitors in treating disabling contractures., Introduction of the common distal arthrogryposis (DA)‐associated missense variant (MYH3R672H) into the analogous zebrafish gene (smyhc1R673H) largely reproduced the human DA phenotype. This model provides key mechanistic insight into the pathogenesis of developmental musculoskeletal disorders as well as exploration of novel therapeutics.

Introduction Congenital clubfoot is a common birth defect that affects at least 0.1% of all births. Nearly 25% cases are familial and the remaining are sporadic in inheritance. Copy number variants (CNVs) involving transcriptional regulators of limb development, including PITX1 and TBX4, have previously been shown to cause familial clubfoot, but much of the heritability remains unexplained. Methods Exome sequence data from 816 unrelated clubfoot cases and 2645 in-house controls were analysed using coverage data to identify rare CNVs. The precise size and location of duplications were then determined using high-density Affymetrix Cytoscan chromosomal microarray (CMA). Segregation in families and de novo status were determined using qantitative PCR. Results Chromosome Xp22.33 duplications involving SHOX were identified in 1.1% of cases (9/816) compared with 0.07% of in-house controls (2/2645) (p=7.98×10−5, OR=14.57) and 0.27% (38/13592) of Atherosclerosis Risk in Communities/the Wellcome Trust Case Control Consortium 2 controls (p=0.001, OR=3.97). CMA validation confirmed an overlapping 180.28 kb duplicated region that included SHOX exons as well as downstream non-coding regions. In four of six sporadic cases where DNA was available for unaffected parents, the duplication was de novo. The probability of four de novo mutations in SHOX by chance in a cohort of 450 sporadic clubfoot cases is 5.4×10–10. Conclusions Microduplications of the pseudoautosomal chromosome Xp22.33 region (PAR1) containing SHOX and downstream enhancer elements occur in ~1% of patients with clubfoot. SHOX and regulatory regions have previously been implicated in skeletal dysplasia as well as idiopathic short stature, but have not yet been reported in clubfoot. SHOX duplications likely contribute to clubfoot pathogenesis by altering early limb development.

Genetic factors predictive of severe adolescent idiopathic scoliosis (AIS) are largely unknown. To identify genetic variation associated with severe AIS, we performed an exome-wide association study of 457 severe AIS cases and 987 controls. We find a missense SNP in SLC39A8 (p.Ala391Thr, rs13107325) associated with severe AIS (P = 1.60 × 10−7, OR = 2.01, CI = 1.54–2.62). This pleiotropic SNP was previously associated with BMI, blood pressure, cholesterol, and blood manganese level. We replicate the association in a second cohort (841 cases and 1095 controls) resulting in a combined P = 7.02 × 10−14, OR = 1.94, CI = 1.63–2.34. Clinically, the minor allele of rs13107325 is associated with greater spinal curvature, decreased height, increased BMI and lower plasma manganese in our AIS cohort. Functional studies demonstrate reduced manganese influx mediated by the SLC39A8 p.Ala391Thr variant and vertebral abnormalities, impaired growth, and decreased motor activity in slc39a8 mutant zebrafish. Our results suggest the possibility that scoliosis may be amenable to dietary intervention., The majority of scoliosis is considered idiopathic with onset in adolescence (AIS) and has a genetic contribution. Here, the authors perform an exome wide association study of data from 457 severe AIS cases and 987 controls, and find a missense variant in SLC39A8 is associated with AIS.

INTRODUCTION: Adolescent idiopathic scoliosis (AIS) is a common musculoskeletal disorder with strong evidence for a genetic contribution. Copy number variants (CNVs) play an important role in congenital scoliosis, but their role in idiopathic scoliosis has been largely unexplored. METHODS: Exome sequence data from 1197 AIS cases and 1664 in-house controls was analyzed using coverage data to identify rare CNVs. CNV calls were filtered to include only highly confident CNVs with > 10 average reads per region and mean log-ratio of coverage consistent with single-copy duplication or deletion. The frequency of 55 common recurrent CNVs was determined and correlated with clinical characteristics. RESULTS: Distal chromosome 16p11.2 microduplications containing the gene SH2B1 were found in 0.7% of AIS cases (8/1197). We replicated this finding in two additional AIS cohorts (8/1097 and 2/433), resulting in 0.7% (18/2727) of all AIS cases harboring a chromosome 16p11.2 microduplication, compared to 0.06% of local controls (1/1664) and 0.04% of published controls (8/19584)(P=2.28×10(−11), OR=16.15). Furthermore, examination of electronic health records of 92,455 patients from the Geisinger health system showed scoliosis in 30% (20/66) patients with chromosome 16p11.2 microduplications containing SH2B1 compared to 7.6% (10/132) of controls (P=5.6×10(−4), OR=3.9). CONCLUSIONS: Recurrent distal chromosome 16p11.2 duplications explain nearly 1% of AIS. Distal chromosome 16p11.2 duplications may contribute to scoliosis pathogenesis by directly impairing growth or by altering expression of nearby genes, such as TBX6. Individuals with distal chromosome 16p11.2 microduplications should be screened for scoliosis to facilitate early treatment.

The 1000 Genomes Project set out to provide a comprehensive description of common human genetic variation by applying whole-genome sequencing to a diverse set of individuals from multiple populations. Here we report completion of the project, having reconstructed the genomes of 2,504 individuals from 26 populations using a combination of low-coverage whole-genome sequencing, deep exome sequencing, and dense microarray genotyping. We characterized a broad spectrum of genetic variation, in total over 88 million variants (84.7 million single nucleotide polymorphisms (SNPs), 3.6 million short insertions/deletions (indels), and 60,000 structural variants), all phased onto high-quality haplotypes. This resource includes >99% of SNP variants with a frequency of >1% for a variety of ancestries. We describe the distribution of genetic variation across the global sample, and discuss the implications for common disease studies., Wellcome Trust (London, England) (Core Award 090532/Z/09/Z), Wellcome Trust (London, England) (Senior Investigator Award 095552/Z/11/Z ), Wellcome Trust (London, England) (WT095908), Wellcome Trust (London, England) (WT109497), Wellcome Trust (London, England) (WT098051), Wellcome Trust (London, England) (WT086084/Z/08/Z), Wellcome Trust (London, England) (WT100956/Z/13/Z ), Wellcome Trust (London, England) (WT097307), Wellcome Trust (London, England) (WT0855322/Z/08/Z ), Wellcome Trust (London, England) (WT090770/Z/09/Z ), Wellcome Trust (London, England) (Major Overseas program in Vietnam grant 089276/Z.09/Z), Medical Research Council (Great Britain) (grant G0801823), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I02593X/1), Biotechnology and Biological Sciences Research Council (Great Britain) (grant BB/I021213/1), Zhongguo ke xue ji shu qing bao yan jiu suo. Office of 863 Programme of China (2012AA02A201), National Basic Research Program of China (2011CB809201), National Basic Research Program of China (2011CB809202), National Basic Research Program of China (2011CB809203), National Natural Science Foundation of China (31161130357), Shenzhen Municipal Government of China (grant ZYC201105170397A), Canadian Institutes of Health Research (grant 136855), Quebec Ministry of Economic Development, Innovation, and Exports (PSR-SIIRI-195), Germany. Bundesministerium für Bildung und Forschung (0315428A), Germany. Bundesministerium für Bildung und Forschung (01GS08201), Germany. Bundesministerium für Bildung und Forschung (BMBF-EPITREAT grant 0316190A), Deutsche Forschungsgemeinschaft (Emmy Noether Grant KO4037/1-1), Beatriu de Pinos Program (2006 BP-A 10144), Beatriu de Pinos Program (2009 BP-B 00274), Spanish National Institute for Health (grant PRB2 IPT13/0001-ISCIII-SGEFI/FEDER), Japan Society for the Promotion of Science (fellowship number PE13075), Marie Curie Actions Career Integration (grant 303772), Fonds National Suisse del la Recherche, SNSF, Scientifique (31003A_130342), National Center for Biotechnology Information (U.S.) (U54HG3067), National Center for Biotechnology Information (U.S.) (U54HG3273), National Center for Biotechnology Information (U.S.) (U01HG5211), National Center for Biotechnology Information (U.S.) (U54HG3079), National Center for Biotechnology Information (U.S.) (R01HG2898), National Center for Biotechnology Information (U.S.) (R01HG2385), National Center for Biotechnology Information (U.S.) (RC2HG5552), National Center for Biotechnology Information (U.S.) (U01HG6513), National Center for Biotechnology Information (U.S.) (U01HG5214), National Center for Biotechnology Information (U.S.) (U01HG5715), National Center for Biotechnology Information (U.S.) (U01HG5718), National Center for Biotechnology Information (U.S.) (U01HG5728), National Center for Biotechnology Information (U.S.) (U41HG7635), National Center for Biotechnology Information (U.S.) (U41HG7497), National Center for Biotechnology Information (U.S.) (R01HG4960), National Center for Biotechnology Information (U.S.) (R01HG5701), National Center for Biotechnology Information (U.S.) (R01HG5214), National Center for Biotechnology Information (U.S.) (R01HG6855), National Center for Biotechnology Information (U.S.) (R01HG7068), National Center for Biotechnology Information (U.S.) (R01HG7644), National Center for Biotechnology Information (U.S.) (DP2OD6514), National Center for Biotechnology Information (U.S.) (DP5OD9154), National Center for Biotechnology Information (U.S.) (R01CA166661), National Center for Biotechnology Information (U.S.) (R01CA172652), National Center for Biotechnology Information (U.S.) (P01GM99568), National Center for Biotechnology Information (U.S.) (R01GM59290), National Center for Biotechnology Information (U.S.) (R01GM104390), National Center for Biotechnology Information (U.S.) (T32GM7790), National Center for Biotechnology Information (U.S.) (R01HL87699), National Center for Biotechnology Information (U.S.) (R01HL104608), National Center for Biotechnology Information (U.S.) (T32HL94284), National Center for Biotechnology Information (U.S.) (HHSN268201100040C), National Center for Biotechnology Information (U.S.) (HHSN272201000025C), Lundbeck Foundation (grant R170-2014-1039, Simons Foundation (SFARI award SF51), National Science Foundation (U.S.) (Research Fellowship DGE-1147470)

Background Identifying insertion/deletion polymorphisms (INDELs) with high confidence has been intrinsically challenging in short-read sequencing data. Here we report our approach for improving INDEL calling accuracy by using a machine learning algorithm to combine call sets generated with three independent methods, and by leveraging the strengths of each individual pipeline. Utilizing this approach, we generated a consensus exome INDEL call set from a large dataset generated by the 1000 Genomes Project (1000G), maximizing both the sensitivity and the specificity of the calls. Results This consensus exome INDEL call set features 7,210 INDELs, from 1,128 individuals across 13 populations included in the 1000 Genomes Phase 1 dataset, with a false discovery rate (FDR) of about 7.0%. Conclusions In our study we further characterize the patterns and distributions of these exonic INDELs with respect to density, allele length, and site frequency spectrum, as well as the potential mutagenic mechanisms of coding INDELs in humans. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1333-7) contains supplementary material, which is available to authorized users.