Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A genome-wide association study identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we use mouse enhancer assays, three mouse enhancer knockouts, and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterize a sequence, PEC7, which overlaps the AIS-associated variant, and find it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, or deletion of both sequences lead to a kinky tail phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicates several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold 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 the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Impaired glycine neurotransmission causes adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is the most common form of spinal deformity, affecting millions of adolescents worldwide, but it lacks a defined theory of etiopathogenesis. Because of this, treatment of AIS is limited to bracing and/or invasive surgery after onset. Preonset diagnosis or preventive treatment remains unavailable. Here, we performed a genetic analysis of a large multicenter AIS cohort and identified disease-causing and predisposing variants of SLC6A9 in multigeneration families, trios, and sporadic patients. Variants of SLC6A9, which encodes glycine transporter 1 (GLYT1), reduced glycine-uptake activity in cells, leading to increased extracellular glycine levels and aberrant glycinergic neurotransmission. Slc6a9 mutant zebrafish exhibited discoordination of spinal neural activities and pronounced lateral spinal curvature, a phenotype resembling human patients. The penetrance and severity of curvature were sensitive to the dosage of functional glyt1. Administration of a glycine receptor antagonist or a clinically used glycine neutralizer (sodium benzoate) partially rescued the phenotype. Our results indicate a neuropathic origin for "idiopathic" scoliosis, involving the dysfunction of synaptic neurotransmission and central pattern generators (CPGs), potentially a common cause of AIS. Our work further suggests avenues for early diagnosis and intervention of AIS in preadolescents.

RAB1A haploinsufficiency phenocopies the 2p14-p15 microdeletion and is associated with impaired neuronal differentiation

Hereditary spastic parapareses (HSPs) are clinically heterogeneous motor neuron diseases with variable age of onset and severity. Although variants in dozens of genes are implicated in HSPs, much of the genetic basis for pediatric-onset HSP remains unexplained. Here, we re-analyzed clinical exome-sequencing data from siblings with HSP of unknown genetic etiology and identified an inherited nonsense mutation (c.523C>T [p.Arg175Ter]) in the highly conserved RAB1A. The mutation is predicted to produce a truncated protein with an intact RAB GTPase domain but without two C-terminal cysteine residues required for proper subcellular protein localization. Additional RAB1A mutations, including two frameshift mutations and a mosaic missense mutation (c.83T>C [p.Leu28Pro]), were identified in three individuals with similar neurodevelopmental presentations. In rescue experiments, production of the full-length, but not the truncated, RAB1a rescued Golgi structure and cell proliferation in Rab1-depleted cells. In contrast, the missense-variant RAB1a disrupted Golgi structure despite intact Rab1 expression, suggesting a dominant-negative function of the mosaic missense mutation. Knock-down of RAB1A in cultured human embryonic stem cell-derived neurons resulted in impaired neuronal arborization. Finally, RAB1A is located within the 2p14-p15 microdeletion syndrome locus. The similar clinical presentations of individuals with RAB1A loss-of-function mutations and the 2p14-p15 microdeletion syndrome implicate loss of RAB1A in the pathogenesis of neurodevelopmental manifestations of this microdeletion syndrome. Our study identifies a RAB1A-related neurocognitive disorder with speech and motor delay, demonstrates an essential role for RAB1a in neuronal differentiation, and implicates RAB1A in the etiology of the neurodevelopmental sequelae associated with the 2p14-p15 microdeletion syndrome.

Interrogating Causal Effects of Body Composition and Puberty-Related Risk Factors on Adolescent Idiopathic Scoliosis: A Two-Sample Mendelian Randomization Study

Adolescent idiopathic scoliosis (AIS) is the most common form of pediatric musculoskeletal disorder. Observational studies have pointed to several risk factors for AIS, but almost no evidence exists to support their causal association with AIS. Here, we applied Mendelian randomization (MR), known to limit bias from confounding and reverse causation, to investigate causal associations between body composition and puberty-related exposures and AIS risk in Europeans and Asians. For our two-sample MR studies, we used single nucleotide polymorphisms (SNPs) associated with body mass index (BMI), waist-hip ratio, lean mass, childhood obesity, bone mineral density (BMD), 25-hydroxyvitamin D (25OHD), age at menarche, and pubertal growth in large European genome-wide association studies (GWAS), and with adult osteoporosis risk and age of menarche in Biobank Japan. We extracted estimates of the aforementioned SNPs on AIS risk from the European or Asian subsets of the largest multiancestry AIS GWAS (N = 7956 cases/88,459 controls). The results of our inverse variance-weighted (IVW) MR estimates suggest no causal association between the aforementioned risk factors and risk of AIS. Pleiotropy-sensitive MR methods yielded similar results. However, restricting our analysis to European females with AIS, we observed a causal association between estimated BMD and the risk of AIS (IVW odds ratio for AIS = 0.1, 95% confidence interval 0.01 to 0.7, p = 0.02 per SD increase in estimated BMD), but this association was no longer significant after adjusting for BMI, body fat mass, and 25OHD and remained significant after adjusting for age at menarche in multivariable MR. In conclusion, we demonstrated a protective causal effect of BMD on AIS risk in females of European ancestry, but this effect was modified by BMI, body fat mass, and 25OHD levels. Future MR studies using larger AIS GWAS are needed to investigate small effects of the aforementioned exposures on AIS. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis

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 the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout 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 in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

Do Children With Medicaid Insurance Have Increased Revision Rates 5 Years After Posterior Spinal Fusions?

BACKGROUND

Socioeconomic disparities in musculoskeletal care are increasingly recognized, however, no studies to date have investigated the role of the insurance carrier on outcomes after posterior spinal fusion (PSF) with segmental spinal instrumentation for adolescent idiopathic scoliosis (AIS).

METHODS

A US insurance dataset was queried using the PearlDiver Mariner software for all patients aged 10 to 18 undergoing PSF for a primary diagnosis of AIS between 2010 and 2020. Age, sex, geographic region, number of levels fused, and baseline medical comorbidities were queried. Complications occurring within 90 days of the index surgery were queried using the International Classification of Diseases, Ninth Revision (ICD-9) and International Classification of Diseases, 10th Revision (ICD-10) codes. Revision surgery was also queried up to 5 years after the index PSF. Categorical variables were compared using the Fisher χ 2 tests and continuous variables were compared using independent t tests. All-cause revision within 5 years was compared using the Kaplan-Meier analysis and a log-rank test. Significance was set at P -value <0.05.

RESULTS

A total of 10,794 patients were identified with 9006 (83.4%) patients with private insurance and 1788 (16.6%) patients insured by Medicaid. The mean follow-up in the database was 5.36±3 years for patients with private insurance and 4.78±2.9 years for patients with Medicaid insurance ( P <0.001). Children with AIS and Medicaid insurance had a significantly higher prevalence of asthma, hypertension, and obesity. A larger percentage of children with Medicaid insurance (41.3%) underwent a ≥13-level PSF compared with privately insured children (34.5%) ( P <0.001). Medicaid patients did not experience higher odds of postoperative complications; in addition, revision surgeries occurred in 1.1% and 1.8% of patients with private insurance and Medicaid insurance, respectively at 5 years postoperatively ( P =0.223).

CONCLUSION

Despite worse baseline comorbidities and longer fusion constructs, AIS patients insured with Medicaid did not have higher rates of complications or revisions at 5-year follow-up versus privately insured patients.

LEVEL OF EVIDENCE

Level III-retrospective cohort study.

Rare variant enrichment analysis supports GREB1L as a contributory driver gene in the etiology of Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is characterized by aplasia of the female reproductive tract; the syndrome can include renal anomalies, absence or dysgenesis, and skeletal anomalies. While functional models have elucidated several candidate genes, only WNT4 (MIM: 603490) variants have been definitively associated with a subtype of MRKH with hyperandrogenism (MIM: 158330). DNA from 148 clinically diagnosed MRKH probands across 144 unrelated families and available family members from North America, Europe, and South America were exome sequenced (ES) and by family-based genomics analyzed for rare likely deleterious variants. A replication cohort consisting of 442 Han Chinese individuals with MRKH was used to further reproduce GREB1L findings in diverse genetic backgrounds. Proband and OMIM phenotypes annotated using the Human Phenotype Ontology were analyzed to quantitatively delineate the phenotypic spectrum associated with GREB1L variant alleles found in our MRKH cohort and those previously published. This study reports 18 novel GREB1L variant alleles, 16 within a multiethnic MRKH cohort and two within a congenital scoliosis cohort. Cohort-wide analyses for a burden of rare variants within a single gene identified likely damaging variants in GREB1L (MIM: 617782), a known disease gene for renal hypoplasia and uterine abnormalities (MIM: 617805), in 16 of 590 MRKH probands. GREB1L variant alleles, including a CNV null allele, were found in 8 MRKH type 1 probands and 8 MRKH type II probands. This study used quantitative phenotypic analyses in a worldwide multiethnic cohort to identify and strengthen the association of GREB1L to isolated uterine agenesis (MRKH type I) and syndromic MRKH type II.

Evidence of causality of low body mass index on risk of adolescent idiopathic scoliosis: a Mendelian randomization study

Introduction

Adolescent idiopathic scoliosis (AIS) is a disorder with a three-dimensional spinal deformity and is a common disease affecting 1-5% of adolescents. AIS is also known as a complex disease involved in environmental and genetic factors. A relation between AIS and body mass index (BMI) has been epidemiologically and genetically suggested. However, the causal relationship between AIS and BMI remains to be elucidated.

Material and methods

Mendelian randomization (MR) analysis was performed using summary statistics from genome-wide association studies (GWASs) of AIS (Japanese cohort, 5,327 cases, 73,884 controls; US cohort: 1,468 cases, 20,158 controls) and BMI (Biobank Japan: 173430 individual; meta-analysis of genetic investigation of anthropometric traits and UK Biobank: 806334 individuals; European Children cohort: 39620 individuals; Population Architecture using Genomics and Epidemiology: 49335 individuals). In MR analyses evaluating the effect of BMI on AIS, the association between BMI and AIS summary statistics was evaluated using the inverse-variance weighted (IVW) method, weighted median method, and Egger regression (MR-Egger) methods in Japanese.

Results

Significant causality of genetically decreased BMI on risk of AIS was estimated: IVW method (Estimate (beta) [SE] = -0.56 [0.16], p = 1.8 × 10-3), weighted median method (beta = -0.56 [0.18], p = 8.5 × 10-3) and MR-Egger method (beta = -1.50 [0.43], p = 4.7 × 10-3), respectively. Consistent results were also observed when using the US AIS summary statistic in three MR methods; however, no significant causality was observed when evaluating the effect of AIS on BMI.

Conclusions

Our Mendelian randomization analysis using large studies of AIS and GWAS for BMI summary statistics revealed that genetic variants contributing to low BMI have a causal effect on the onset of AIS. This result was consistent with those of epidemiological studies and would contribute to the early detection of AIS.

Deletion of Pax1 scoliosis-associated regulatory elements leads to a female-biased tail abnormality

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is sexually dimorphic, with increased incidence in females. A GWAS identified a female-specific AIS susceptibility locus near the PAX1 gene. Here, we used mouse enhancer assays, three mouse enhancer knockouts and subsequent phenotypic analyses to characterize this region. Using mouse enhancer assays, we characterized a sequence, PEC7, that overlaps the AIS-associated variant, and found it to be active in the tail tip and intervertebral disc. Removal of PEC7 or Xe1, a known sclerotome enhancer nearby, and deletion of both sequences led to a kinky phenotype only in the Xe1 and combined (Xe1+PEC7) knockouts, with only the latter showing a female sex dimorphic phenotype. Extensive phenotypic characterization of these mouse lines implicated several differentially expressed genes and estrogen signaling in the sex dimorphic bias. In summary, our work functionally characterizes an AIS-associated locus and dissects the mechanism for its sexual dimorphism.

Molecular Dissection of Somatic Skeletal Disease in Neurofibromatosis Type 1

Neurofibromatosis type 1 (NF1) is a tumor predisposition syndrome caused by heterozygous NF1 gene mutations. Patients with NF1 present with pleiotropic somatic secondary manifestations, including development of bone pseudarthrosis after fracture. Somatic NF1 gene mutations were reproducibly identified in patient-derived pseudarthrosis specimens, suggesting a local mosaic cell population including somatic pathologic cells. The somatic cellular pathogenesis of NF1 pseudarthroses remains unclear, though defects in osteogenesis have been posited. Here, we applied time-series single-cell RNA-sequencing (scRNA-seq) to patient-matched control and pseudarthrosis-derived primary bone stromal cells (BSCs). We show that osteogenic specification to an osteoblast progenitor cell population was evident for control bone-derived cells and haploinsufficient pseudarthrosis-derived cells. Similar results were observed for somatic patient fracture-derived NF1^-/- cells; however, expression of genetic pathways associated with skeletal mineralization were significantly reduced in NF1^-/- cells compared with fracture-derived NF1^+/- cells. In mice, we show that Nf1 expressed in bone marrow osteoprogenitors is required for the maintenance of the adult skeleton. Results from our study implicate impaired Clec11a-Itga11-Wnt signaling in the pathogenesis of NF1-associated skeletal disease. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Germline Saturation Mutagenesis Induces Skeletal Phenotypes in Mice

Proper embryonic and postnatal skeletal development require coordination of myriad complex molecular mechanisms. Disruption of these processes, through genetic mutation, contributes to variation in skeletal development. We developed a high-throughput N-ethyl-N-nitrosourea (ENU)-induced saturation mutagenesis skeletal screening approach in mice to identify genes required for proper skeletal development. Here, we report initial results from live-animal X-ray and dual-energy X-ray absorptiometry (DXA) imaging of 27,607 G3 mice from 806 pedigrees, testing the effects of 32,198 coding/splicing mutations in 13,020 genes. A total of 39.7% of all autosomal genes were severely damaged or destroyed by mutations tested twice or more in the homozygous state. Results from our study demonstrate the feasibility of in vivo mutagenesis to identify mouse models of skeletal disease. Furthermore, our study demonstrates how ENU mutagenesis provides opportunities to create and characterize putative hypomorphic mutations in developmentally essential genes. Finally, we present a viable mouse model and case report of recessive skeletal disease caused by mutations in FAM20B. Results from this study, including engineered mouse models, are made publicly available via the online Mutagenetix database. © 2021 American Society for Bone and Mineral Research (ASBMR).

Thousands of induced germline mutations affecting immune cells identified by automated meiotic mapping coupled with machine learning

Forward genetic studies use meiotic mapping to adduce evidence that a particular mutation, normally induced by a germline mutagen, is causative of a particular phenotype. Particularly in small pedigrees, cosegregation of multiple mutations, occasional unawareness of mutations, and paucity of homozygotes may lead to erroneous declarations of cause and effect. We sought to improve the identification of mutations causing immune phenotypes in mice by creating Candidate Explorer (CE), a machine-learning software program that integrates 67 features of genetic mapping data into a single numeric score, mathematically convertible to the probability of verification of any putative mutation-phenotype association. At this time, CE has evaluated putative mutation-phenotype associations arising from screening damaging mutations in ∼55% of mouse genes for effects on flow cytometry measurements of immune cells in the blood. CE has therefore identified more than half of genes within which mutations can be causative of flow cytometric phenovariation in Mus musculus The majority of these genes were not previously known to support immune function or homeostasis. Mouse geneticists will find CE data informative in identifying causative mutations within quantitative trait loci, while clinical geneticists may use CE to help connect causative variants with rare heritable diseases of immunity, even in the absence of linkage information. CE displays integrated mutation, phenotype, and linkage data, and is freely available for query online.

Saturation mutagenesis defines novel mouse models of severe spine deformity

Embryonic formation and patterning of the vertebrate spinal column requires coordination of many molecular cues. After birth, the integrity of the spine is impacted by developmental abnormalities of the skeletal, muscular and nervous systems, which may result in deformities, such as kyphosis and scoliosis. We sought to identify novel genetic mouse models of severe spine deformity by implementing in vivo skeletal radiography as part of a high-throughput saturation mutagenesis screen. We report selected examples of genetic mouse models following radiographic screening of 54,497 mice from 1275 pedigrees. An estimated 30.44% of autosomal genes harbored predicted damaging alleles examined twice or more in the homozygous state. Of the 1275 pedigrees screened, 7.4% presented with severe spine deformity developing in multiple mice, and of these, meiotic mapping implicated N-ethyl-N-nitrosourea alleles in 21% of pedigrees. Our study provides proof of concept that saturation mutagenesis is capable of discovering novel mouse models of human disease, including conditions with skeletal, neural and neuromuscular pathologies. Furthermore, we report a mouse model of skeletal disease, including severe spine deformity, caused by recessive mutation in Scube3. By integrating results with a human clinical exome database, we identified a patient with undiagnosed skeletal disease who harbored recessive mutations in SCUBE3, and we demonstrated that disease-associated mutations are associated with reduced transactivation of Smad signaling in vitro. All radiographic results and mouse models are made publicly available through the Mutagenetix online database with the goal of advancing understanding of spine development and discovering novel mouse models of human disease.

Summary: We report selected mouse models of spine deformity following mutagenesis across 30% of autosomal genes, results of which are made publicly available to advance understanding of spine development and disease.

Genomic characterization of the adolescent idiopathic scoliosis-associated transcriptome and regulome

Adolescent idiopathic scoliosis (AIS), a sideways curvature of the spine, is the most common pediatric musculoskeletal disorder, affecting ~3% of the population worldwide. However, its genetic bases and tissues of origin remain largely unknown. Several genome-wide association studies (GWAS) have implicated nucleotide variants in non-coding sequences that control genes with important roles in cartilage, muscle, bone, connective tissue and intervertebral disks (IVDs) as drivers of AIS susceptibility. Here, we set out to define the expression of AIS-associated genes and active regulatory elements by performing RNA-seq and chromatin immunoprecipitation-sequencing against H3 lysine 27 acetylation in these tissues in mouse and human. Our study highlights genetic pathways involving AIS-associated loci that regulate chondrogenesis, IVD development and connective tissue maintenance and homeostasis. In addition, we identify thousands of putative AIS-associated regulatory elements which may orchestrate tissue-specific expression in musculoskeletal tissues of the spine. Quantification of enhancer activity of several candidate regulatory elements from our study identifies three functional enhancers carrying AIS-associated GWAS SNPs at the ADGRG6 and BNC2 loci. Our findings provide a novel genome-wide catalog of AIS-relevant genes and regulatory elements and aid in the identification of novel targets for AIS causality and treatment.

Novel homozygous variant in WISP3 in a family with unrecognized progressive pseudorheumatoid dysplasia

We present the use of whole-genome sequencing to correctly diagnose progressive pseudorheumatoid dysplasia in patients with atypical clinical and radiologic findings and prior diagnosis of juvenile idiopathic arthritis.

Diagnostic yield and clinical impact of exome sequencing in early-onset scoliosis (EOS)

BACKGROUND

Early-onset scoliosis (EOS), defined by an onset age of scoliosis less than 10 years, conveys significant health risk to affected children. Identification of the molecular aetiology underlying patients with EOS could provide valuable information for both clinical management and prenatal screening.

METHODS

In this study, we consecutively recruited a cohort of 447 Chinese patients with operative EOS. We performed exome sequencing (ES) screening on these individuals and their available family members (totaling 670 subjects). Another cohort of 13 patients with idiopathic early-onset scoliosis (IEOS) from the USA who underwent ES was also recruited.

RESULTS

After ES data processing and variant interpretation, we detected molecular diagnostic variants in 92 out of 447 (20.6%) Chinese patients with EOS, including 8 patients with molecular confirmation of their clinical diagnosis and 84 patients with molecular diagnoses of previously unrecognised diseases underlying scoliosis. One out of 13 patients with IEOS from the US cohort was molecularly diagnosed. The age at presentation, the number of organ systems involved and the Cobb angle were the three top features predictive of a molecular diagnosis.

CONCLUSION

ES enabled the molecular diagnosis/classification of patients with EOS. Specific clinical features/feature pairs are able to indicate the likelihood of gaining a molecular diagnosis through ES.

The cartilage matrisome in adolescent idiopathic scoliosis

The human spinal column is a dynamic, segmented, bony, and cartilaginous structure that protects the neurologic system and simultaneously provides balance and flexibility. Children with developmental disorders that affect the patterning or shape of the spine can be at risk of neurologic and other physiologic dysfunctions. The most common developmental disorder of the spine is scoliosis, a lateral deformity in the shape of the spinal column. Scoliosis may be part of the clinical spectrum that is observed in many developmental disorders, but typically presents as an isolated symptom in otherwise healthy adolescent children. Adolescent idiopathic scoliosis (AIS) has defied understanding in part due to its genetic complexity. Breakthroughs have come from recent genome-wide association studies (GWAS) and next generation sequencing (NGS) of human AIS cohorts, as well as investigations of animal models. These studies have identified genetic associations with determinants of cartilage biogenesis and development of the intervertebral disc (IVD). Current evidence suggests that a fraction of AIS cases may arise from variation in factors involved in the structural integrity and homeostasis of the cartilaginous extracellular matrix (ECM). Here, we review the development of the spine and spinal cartilages, the composition of the cartilage ECM, the so-called "matrisome" and its functions, and the players involved in the genetic architecture of AIS. We also propose a molecular model by which the cartilage matrisome of the IVD contributes to AIS susceptibility.

TBX6 missense variants expand the mutational spectrum in a non-Mendelian inheritance disease

Congenital scoliosis (CS) is a birth defect with variable clinical and anatomical manifestations due to spinal malformation. The genetic etiology underlying about 10% of CS cases in the Chinese population is compound inheritance by which the gene dosage is reduced below that of haploinsufficiency. In this genetic model, the trait manifests as a result of the combined effect of a rare variant and common pathogenic variant allele at a locus. From exome sequencing (ES) data of 523 patients in Asia and two patients in Texas, we identified six TBX6 gene-disruptive variants from 11 unrelated CS patients via ES and in vitro functional testing. The in trans mild hypomorphic allele was identified in 10 of the 11 subjects; as anticipated these 10 shared a similar spinal deformity of hemivertebrae. The remaining case has a homozygous variant in TBX6 (c.418C>T) and presents a more severe spinal deformity phenotype. We found decreased transcriptional activity and abnormal cellular localization as the molecular mechanisms for TBX6 missense loss-of-function alleles. Expanding the mutational spectrum of TBX6 pathogenic alleles enabled an increased molecular diagnostic detection rate, provided further evidence for the gene dosage-dependent genetic model underlying CS, and refined clinical classification.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Distal chromosome 16p11.2 duplications containing SH2B1 in patients with scoliosis

INTRODUCTION

Adolescent 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.

METHODS

Exome 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.

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 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).

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.

A multi-ethnic meta-analysis confirms the association of rs6570507 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is the most common type of spinal deformity and has a significant genetic background. Genome-wide association studies (GWASs) identified several susceptibility loci associated with AIS. Among them is a locus on chromosome 6q24.1 that we identified by a GWAS in a Japanese cohort. The locus is represented by rs6570507 located within GPR126. To ensure the association of rs6570507 with AIS, we conducted a meta-analysis using eight cohorts from East Asia, Northern Europe and USA. The analysis included a total of 6,873 cases and 38,916 controls and yielded significant association (combined P = 2.95 × 10^-20; odds ratio = 1.22), providing convincing evidence of the worldwide association between rs6570507 and AIS susceptibility. In silico analyses strongly suggested that GPR126 is a susceptibility gene at this locus.

An international meta-analysis confirms the association of BNC2 with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a common spinal deformity with the prevalence of approximately 3%. We previously conducted a genome-wide association study (GWAS) using a Japanese cohort and identified a novel locus on chromosome 9p22.2. However, a replication study using multi-population cohorts has not been conducted. To confirm the association of 9p22.2 locus with AIS in multi-ethnic populations, we conducted international meta-analysis using eight cohorts. In total, we analyzed 8,756 cases and 27,822 controls. The analysis showed a convincing evidence of association between rs3904778 and AIS. Seven out of eight cohorts had significant P value, and remaining one cohort also had the same trend as the seven. The combined P was 3.28 × 10⁻¹⁸ (odds ratio = 1.19, 95% confidence interval = 1.14–1.24). In silico analyses suggested that BNC2 is the AIS susceptibility gene in this locus.

Summary of the first inaugural joint meeting of the International Consortium for scoliosis genetics and the International Consortium for vertebral anomalies and scoliosis, March 16-18, 2017, Dallas, Texas

Scoliosis represents the most common musculoskeletal disorder in children and affects approximately 3% of the world population. Scoliosis is separated into two major phenotypic classifications: congenital and idiopathic. Idiopathic scoliosis is defined as a curvature of the spine of 10° or greater visualized on plane radiograph and does not have associated vertebral malformations (VM). "Congenital" scoliosis (CS) due to malformations in vertebrae is frequently associated with other birth defects. Recently, significant advances have been made in understanding the genetic basis of both conditions. There is evidence that both conditions are etiologically related. A 2-day conference entitled "Genomic Approaches to Understanding and Treating Scoliosis" was held at Scottish Rite Hospital for Children in Dallas, Texas, to synergize research in this field. This first combined, multidisciplinary conference featured international scoliosis researchers in basic and clinical sciences. A major outcome of the conference advancing scoliosis research was the proposal and subsequent vote in favor of merging the International Consortium for Vertebral Anomalies and Scoliosis (ICVAS) and International Consortium for Scoliosis Genetics (ICSG) into a single entity called International Consortium for Spinal Genetics, Development, and Disease (ICSGDD). The ICSGDD is proposed to meet annually as a forum to synergize multidisciplinary spine deformity research.

Genome-Wide Association Meta-Analysis Reveals Novel Juvenile Idiopathic Arthritis Susceptibility Loci

OBJECTIVE

Juvenile idiopathic arthritis (JIA) is the most common childhood rheumatic disease and has a strong genomic component. To date, JIA genetic association studies have had limited sample sizes, used heterogeneous patient populations, or included only candidate regions. The aim of this study was to identify new associations between JIA patients with oligoarticular disease and those with IgM rheumatoid factor (RF)-negative polyarticular disease, which are clinically similar and the most prevalent JIA disease subtypes.

METHODS

Three cohorts comprising 2,751 patients with oligoarticular or RF-negative polyarticular JIA were genotyped using the Affymetrix Genome-Wide SNP Array 6.0 or the Illumina HumanCoreExome-12+ Array. Overall, 15,886 local and out-of-study controls, typed on these platforms or the Illumina HumanOmni2.5, were used for association analyses. High-quality single-nucleotide polymorphisms (SNPs) were used for imputation to 1000 Genomes prior to SNP association analysis.

RESULTS

Meta-analysis showed evidence of association (P < 1 × 10^-6 ) at 9 regions: PRR9_LOR (P = 5.12 × 10^-8 ), ILDR1_CD86 (P = 6.73 × 10^-8 ), WDFY4 (P = 1.79 × 10^-7 ), PTH1R (P = 1.87 × 10^-7 ), RNF215 (P = 3.09 × 10^-7 ), AHI1_LINC00271 (P = 3.48 × 10^-7 ), JAK1 (P = 4.18 × 10^-7 ), LINC00951 (P = 5.80 × 10^-7 ), and HBP1 (P = 7.29 × 10^-7 ). Of these, PRR9_LOR, ILDR1_CD86, RNF215, LINC00951, and HBP1 were shown, for the first time, to be autoimmune disease susceptibility loci. Furthermore, associated SNPs included cis expression quantitative trait loci for WDFY4, CCDC12, MTP18, SF3A1, AHI1, COG5, HBP1, and GPR22.

CONCLUSION

This study provides evidence of both unique JIA risk loci and risk loci overlapping between JIA and other autoimmune diseases. These newly associated SNPs are shown to influence gene expression, and their bounding regions tie into molecular pathways of immunologic relevance. Thus, they likely represent regions that contribute to the pathology of oligoarticular JIA and RF-negative polyarticular JIA.

Fine-mapping the MHC locus in juvenile idiopathic arthritis (JIA) reveals genetic heterogeneity corresponding to distinct adult inflammatory arthritic diseases

Objectives

Juvenile idiopathic arthritis (JIA) is a heterogeneous group of diseases, comprising seven categories. Genetic data could potentially be used to help redefine JIA categories and improve the current classification system. The human leucocyte antigen (HLA) region is strongly associated with JIA. Fine-mapping of the region was performed to look for similarities and differences in HLA associations between the JIA categories and define correspondences with adult inflammatory arthritides.

Methods

Dense genotype data from the HLA region, from the Immunochip array for 5043 JIA cases and 14 390 controls, were used to impute single-nucleotide polymorphisms, HLA classical alleles and amino acids. Bivariate analysis was performed to investigate genetic correlation between the JIA categories. Conditional analysis was used to identify additional effects within the region. Comparison of the findings with those in adult inflammatory arthritic diseases was performed.

Results

We identified category-specific associations and have demonstrated for the first time that rheumatoid factor (RF)-negative polyarticular JIA and oligoarticular JIA are genetically similar in their HLA associations. We also observe that each JIA category potentially has an adult counterpart. The RF-positive polyarthritis association at HLA-DRB1 amino acid at position 13 mirrors the association in adult seropositive rheumatoid arthritis (RA). Interestingly, the combined oligoarthritis and RF-negative polyarthritis dataset shares the same association with adult seronegative RA.

Conclusions

The findings suggest the value of using genetic data in helping to classify the categories of this heterogeneous disease. Mapping JIA categories to adult counterparts could enable shared knowledge of disease pathogenesis and aetiology and facilitate transition from paediatric to adult services.

Mutations in genes encoding condensin complex proteins cause microcephaly through decatenation failure at mitosis

Martin et al. report that biallelic mutations in NCAPD2, NCAPH, or NCAPD3, encoding subunits of condensin complexes, cause microcephaly. Frequent anaphase chromatin bridge formation observed in apical neural progenitors during neurogenesis are the consequence of failed sister chromatid disentanglement during chromosome compaction.

Compaction of chromosomes is essential for accurate segregation of the genome during mitosis. In vertebrates, two condensin complexes ensure timely chromosome condensation, sister chromatid disentanglement, and maintenance of mitotic chromosome structure. Here, we report that biallelic mutations in NCAPD2, NCAPH, or NCAPD3, encoding subunits of these complexes, cause microcephaly. In addition, hypomorphic Ncaph2 mice have significantly reduced brain size, with frequent anaphase chromatin bridge formation observed in apical neural progenitors during neurogenesis. Such DNA bridges also arise in condensin-deficient patient cells, where they are the consequence of failed sister chromatid disentanglement during chromosome compaction. This results in chromosome segregation errors, leading to micronucleus formation and increased aneuploidy in daughter cells. These findings establish “condensinopathies” as microcephalic disorders, with decatenation failure as an additional disease mechanism for microcephaly, implicating mitotic chromosome condensation as a key process ensuring mammalian cerebral cortex size.

Regulatory polymorphisms modulate the expression of HLA class II molecules and promote autoimmunity

Targeted sequencing of sixteen SLE risk loci among 1349 Caucasian cases and controls produced a comprehensive dataset of the variations causing susceptibility to systemic lupus erythematosus (SLE). Two independent disease association signals in the HLA-D region identified two regulatory regions containing 3562 polymorphisms that modified thirty-seven transcription factor binding sites. These extensive functional variations are a new and potent facet of HLA polymorphism. Variations modifying the consensus binding motifs of IRF4 and CTCF in the XL9 regulatory complex modified the transcription of HLA-DRB1, HLA-DQA1 and HLA-DQB1 in a chromosome-specific manner, resulting in a 2.5-fold increase in the surface expression of HLA-DR and DQ molecules on dendritic cells with SLE risk genotypes, which increases to over 4-fold after stimulation. Similar analyses of fifteen other SLE risk loci identified 1206 functional variants tightly linked with disease-associated SNPs and demonstrated that common disease alleles contain multiple causal variants modulating multiple immune system genes.

Meta-analysis of shared genetic architecture across ten pediatric autoimmune diseases

Genome-wide association studies (GWASs) have identified hundreds of susceptibility genes, including shared associations across clinically distinct autoimmune diseases. We performed an inverse χ(2) meta-analysis across ten pediatric-age-of-onset autoimmune diseases (pAIDs) in a case-control study including more than 6,035 cases and 10,718 shared population-based controls. We identified 27 genome-wide significant loci associated with one or more pAIDs, mapping to in silico-replicated autoimmune-associated genes (including IL2RA) and new candidate loci with established immunoregulatory functions such as ADGRL2, TENM3, ANKRD30A, ADCY7 and CD40LG. The pAID-associated single-nucleotide polymorphisms (SNPs) were functionally enriched for deoxyribonuclease (DNase)-hypersensitivity sites, expression quantitative trait loci (eQTLs), microRNA (miRNA)-binding sites and coding variants. We also identified biologically correlated, pAID-associated candidate gene sets on the basis of immune cell expression profiling and found evidence of genetic sharing. Network and protein-interaction analyses demonstrated converging roles for the signaling pathways of type 1, 2 and 17 helper T cells (TH1, TH2 and TH17), JAK-STAT, interferon and interleukin in multiple autoimmune diseases.

A PAX1 enhancer locus is associated with susceptibility to idiopathic scoliosis in females

Idiopathic scoliosis (IS) is a common paediatric musculoskeletal disease that displays a strong female bias. By performing a genome-wide association study (GWAS) of 3,102 individuals, we identify significant associations with 20p11.22 SNPs for females (P=6.89 × 10⁻⁹) but not males (P=0.71). This association with IS is also found in independent female cohorts from the United States of America and Japan (overall P=2.15 × 10⁻¹⁰, OR=1.30 (rs6137473)). Unexpectedly, the 20p11.22 IS risk alleles were previously associated with protection from early-onset alopecia, another sexually dimorphic condition. The 174-kb associated locus is distal to PAX1, which encodes paired box 1, a transcription factor involved in spine development. We identify a sequence in the associated locus with enhancer activity in zebrafish somitic muscle and spinal cord, an activity that is abolished by IS-associated SNPs. We thus identify a sexually dimorphic IS susceptibility locus, and propose the first functionally defined candidate mutations in an enhancer that may regulate expression in specific spinal cells.

Girls are tenfold more likely than boys to require surgical treatment for idiopathic scoliosis, a common paediatric skeletal disorder. Here, Sharma et al. identify the first sexually dimorphic idiopathic scoliosis risk locus, and demonstrate that it may play a role in the regulation of spinal cells.

Neurofibromin deficiency-associated transcriptional dysregulation suggests a novel therapy for tibial pseudoarthrosis in NF1

Neurofibromatosis type 1 (NF1) is an autosomal dominant disease caused by mutations in NF1. Among the earliest manifestations is tibial pseudoarthrosis and persistent nonunion after fracture. To further understand the pathogenesis of pseudoarthrosis and the underlying bone remodeling defect, pseudoarthrosis tissue and cells cultured from surgically resected pseudoarthrosis tissue from NF1 individuals were analyzed using whole-exome and whole-transcriptome sequencing as well as genomewide microarray analysis. Genomewide analysis identified multiple genetic mechanisms resulting in somatic biallelic NF1 inactivation; no other genes with recurring somatic mutations were identified. Gene expression profiling identified dysregulated pathways associated with neurofibromin deficiency, including phosphoinositide 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling pathways. Unlike aggressive NF1-associated malignancies, tibial pseudoarthrosis tissue does not harbor a high frequency of somatic mutations in oncogenes or other tumor-suppressor genes, such as p53. However, gene expression profiling indicates that pseudoarthrosis tissue has a tumor-promoting transcriptional pattern, despite lacking tumorigenic somatic mutations. Significant overexpression of specific cancer-associated genes in pseudoarthrosis highlights a potential for receptor tyrosine kinase inhibitors to target neurofibromin-deficient pseudoarthrosis and promote proper bone remodeling and fracture healing.

Genomic Analyses of Patients With Unexplained Early-Onset Scoliosis

STUDY DESIGN

To test for rare genetic mutations, a cohort of patients with unexplained early-onset scoliosis (EOS) was screened using high-density microarray genotyping. A cohort of patients with adolescent idiopathic scoliosis (AIS) was similarly screened and the results were compared.

SUMMARY OF BACKGROUND DATA

Patients with scoliosis in infancy or early childhood (EOS) are at high risk for progressive deformity and associated problems including respiratory compromise. Early-onset scoliosis is frequently associated with genetic disorders but many patients present with nonspecific clinical features and without an associated diagnosis. The authors hypothesized that EOS in these patients may be caused by rare genetic mutations detectable by next-generation genomic methods.

METHODS

The researchers identified 24 patients with unexplained EOS from pediatric orthopedic clinics. They genotyped them, along with 39 connecting family members, using the Illumina OmniExpress-12, version 1.0 beadchip. Resulting genotypes were analyzed for chromosomal changes, specifically copy number variation and absence of heterozygosity. They screened 482 adolescent idiopathic scoliosis (AIS) patients and 744 healthy controls, who were similarly genotyped with the same beadchip, for chromosomal changes identified in the EOS cohort.

RESULTS

Copy number variation and absence of heterozygosity analyses revealed a genetic diagnosis of chromosome 15q24 microdeletion syndrome in 1 patient and maternal uniparental disomy of chromosome 14 in a second one. Prior genetic testing and clinical evaluations had been negative in both cases. A large novel chromosome 10 deletion was likely causal in a third EOS patient. These mutations identified in the EOS patients were absent in AIS patients and controls, and thus were not associated with AIS or found in asymptomatic individuals.

CONCLUSIONS

These data underscore the usefulness of updated genetic evaluations including high-density microarray-based genotyping and other next-generation methods in patients with unexplained EOS, even when prior genetic studies were negative. These data also suggest the intriguing possibility that other mutations detectable by whole genome sequencing, as well as epigenetic effects, await discovery in the EOS population.

Extreme growth failure is a common presentation of ligase IV deficiency

Ligase IV syndrome is a rare differential diagnosis for Nijmegen breakage syndrome owing to a shared predisposition to lympho-reticular malignancies, significant microcephaly, and radiation hypersensitivity. Only 16 cases with mutations in LIG4 have been described to date with phenotypes varying from malignancy in developmentally normal individuals, to severe combined immunodeficiency and early mortality. Here, we report the identification of biallelic truncating LIG4 mutations in 11 patients with microcephalic primordial dwarfism presenting with restricted prenatal growth and extreme postnatal global growth failure (average OFC -10.1 s.d., height -5.1 s.d.). Subsequently, most patients developed thrombocytopenia and leucopenia later in childhood and many were found to have previously unrecognized immunodeficiency following molecular diagnosis. None have yet developed malignancy, though all patients tested had cellular radiosensitivity. A genotype-phenotype correlation was also noted with position of truncating mutations corresponding to disease severity. This work extends the phenotypic spectrum associated with LIG4 mutations, establishing that extreme growth retardation with microcephaly is a common presentation of bilallelic truncating mutations. Such growth failure is therefore sufficient to consider a diagnosis of LIG4 deficiency and early recognition of such cases is important as bone marrow failure, immunodeficiency, and sometimes malignancy are long term sequelae of this disorder.

Susceptibility to childhood-onset rheumatoid arthritis: investigation of a weighted genetic risk score that integrates cumulative effects of variants at five genetic loci

OBJECTIVE

Children with childhood-onset rheumatoid arthritis (RA) include those with rheumatoid factor or anti-citrullinated protein antibody-positive juvenile idiopathic arthritis. To test the hypothesis that adult-onset RA-associated variants are also associated with childhood-onset RA, we investigated RA-associated variants at 5 loci in a cohort of patients with childhood-onset RA. We also assessed the cumulative association of these variants in susceptibility to childhood-onset RA using a weighted genetic risk score (wGRS).

METHODS

A total of 155 children with childhood-onset RA and 684 healthy controls were genotyped for 5 variants in the PTPN22, TRAF1/C5, STAT4, and TNFAIP3 loci. High-resolution HLA-DRB1 genotypes were available for 149 cases and 373 controls. We tested each locus for association with childhood-onset RA via logistic regression. We also computed a wGRS for each subject, with weights based on the natural log of the published odds ratios (ORs) for the alleles investigated, and used logistic regression to test the wGRS for association with childhood-onset RA.

RESULTS

Childhood-onset RA was associated with TNFAIP3 rs10499194 (OR 0.60 [95% confidence interval 0.44-0.83]), PTPN22 rs2476601 (OR 1.61 [95% confidence interval 1.11-2.31]), and STAT4 rs7574865 (OR 1.41 [95% confidence interval 1.06-1.87]) variants. The wGRS was significantly different between cases and controls (P < 2 × 10(-16) ). Individuals in the third to fifth quintiles of wGRS had a significantly increased disease risk compared to baseline (individuals in the first quintile). Higher wGRS was associated with increased risk of childhood-onset RA, especially among males.

CONCLUSION

The magnitude and direction of the association between TNFAIP3, STAT4, and PTPN22 variants and childhood-onset RA are similar to those observed in RA, suggesting that adult-onset RA and childhood-onset RA share common genetic risk factors. Using a wGRS, we have demonstrated the cumulative association of RA-associated variants with susceptibility to childhood-onset RA.

Genetic variants in GPR126 are associated with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is the most common pediatric skeletal disease. We previously reported a locus on chromosome 10q24.31 associated with AIS susceptibility in Japanese using a genome-wide association study (GWAS) consisting of 1,033 cases and 1,473 controls. To identify additional AIS-associated loci, we expanded the study by adding X-chromosome SNPs in the GWAS and increasing the size of the replication cohorts. Through a stepwise association study including 1,819 cases and 25,939 controls, we identified a new susceptibility locus on chromosome 6q24.1 in Japanese (P = 2.25 × 10(-10); odds ratio (OR) = 1.28). The most significantly associated SNP, rs6570507, was in GPR126 (encoding G protein-coupled receptor 126). Its association was replicated in Han Chinese and European-ancestry populations (combined P = 1.27 × 10(-14); OR = 1.27). GPR126 was highly expressed in cartilage, and the knockdown of gpr126 in zebrafish caused delayed ossification of the developing spine. Our results should provide insights into the etiology and pathogenesis of AIS.

Identification of a homozygous PSTPIP1 mutation in a patient with a PAPA-like syndrome responding to canakinumab treatment

BACKGROUND

Pyogenic sterile arthritis, pyoderma gangrenosum, and acne (PAPA) syndrome (OMIM 604416) is a rare autosomal dominant inherited autoinflammatory syndrome characterized by pyogenic sterile arthritis and less frequently accompanied by pyoderma gangrenosum and acne. It is associated with dominant missense mutations in the proline-serine-threonine phosphatase-interacting protein 1 gene (PSTPIP1) located on chromosome 15. The patient was diagnosed as having features of a PAPA-like syndrome in which cutaneous manifestations, such as pyoderma gangrenosum and acne fulminans, predominated.

OBSERVATIONS

Sequencing of the PSTPIP1 gene was performed in the patient and his extended family. The patient's DNA analysis revealed a homozygous nucleotide exchange c.773G>C in the PSTPIP1 gene, leading to the substitution of glycine 258 by alanine (p.Gly258Ala), a previously reported heterozygous polymorphism. Heterozygous changes were identified in both of the patient's parents and in 7 other family members, all of whom were asymptomatic. The patient was treated with canakinumab, a human anti-interleukin 1β monoclonal antibody, which led to rapid remission of the symptoms.

CONCLUSIONS

To our knowledge, this is the first reported case of the resolution of dermatological symptoms associated with a PAPA-like syndrome using canakinumab treatment. Further study of the p.Gly258Ala variant is warranted to determine whether this mutation has a role in causing an apparently recessive cutaneous syndrome resembling PAPA syndrome.

A novel method for analyzing genetic association with longitudinal phenotypes

Knowledge of genes influencing longitudinal patterns may offer information about predicting disease progression. We developed a systematic procedure for testing association between SNP genotypes and longitudinal phenotypes. We evaluated false positive rates and statistical power to localize genes for disease progression. We used genome-wide SNP data from the Framingham Heart Study. With longitudinal data from two real studies unrelated to Framingham, we estimated three trajectory curves from each study. We performed simulations by randomly selecting 500 individuals. In each simulation replicate, we assigned each individual to one of the three trajectory groups based on the underlying hypothesis (null or alternative), and generated corresponding longitudinal data. Individual Bayesian posterior probabilities (BPPs) for belonging to a specific trajectory curve were estimated. These BPPs were treated as a quantitative trait and tested (using the Wald test) for genome-wide association. Empirical false positive rates and power were calculated. Our method maintained the expected false positive rate for all simulation models. Also, our method achieved high empirical power for most simulations. Our work presents a method for disease progression gene mapping. This method is potentially clinically significant as it may allow doctors to predict disease progression based on genotype and determine treatment accordingly.

Somatic gain-of-function mutations in PIK3CA in patients with macrodactyly

Macrodactyly is a discrete congenital anomaly consisting of enlargement of all tissues localized to the terminal portions of a limb, typically within a ‘nerve territory’. The classic terminology for this condition is ‘lipofibromatous hamartoma of nerve’ or Type I macrodactyly. The peripheral nerve, itself, is enlarged both in circumference and in length. It is not related to neurofibromatosis (NF1), nor is it associated with vascular malformations, such as in the recently reported CLOVES syndrome. The specific nerve pathophysiology in this form of macrodactyly has not been well described and a genetic etiology for this specific form of enlargement is unknown. To identify the genetic cause of macrodactyly, we used whole-exome sequencing to identify somatic mutations present in the affected nerve of a single patient. We confirmed a novel mutation in PIK3CA (R115P) present in the patient's affected nerve tissue but not in blood DNA. Sequencing PIK3CA exons identified gain-of-function mutations (E542K, H1047L or H1047R) in the affected tissue of five additional unrelated patients; mutations were absent in blood DNA available from three patients. Immunocytochemistry confirmed AKT activation in cultured cells from the nerve of a macrodactyly patient. Additionally, we found that the most abnormal structure within the involved nerve in a macrodactylous digit is the perineurium, with additional secondary effects on the axon number and size. Thus, isolated congenital macrodactyly is caused by somatic activation of the PI3K/AKT cell-signaling pathway and is genetically and biochemically related to other overgrowth syndromes.

Meier-Gorlin syndrome: growth and secondary sexual development of a microcephalic primordial dwarfism disorder

Meier-Gorlin syndrome (MGS) is a rare autosomal recessive disorder characterized by primordial dwarfism, microtia, and patellar aplasia/hypoplasia. Recently, mutations in the ORC1, ORC4, ORC6, CDT1, and CDC6 genes, encoding components of the pre-replication complex, have been identified. This complex is essential for DNA replication and therefore mutations are expected to impair cell proliferation and consequently could globally reduce growth. However, detailed growth characteristics of MGS patients have not been reported, and so this is addressed here through study of 45 MGS patients, the largest cohort worldwide. Here, we report that growth velocity (length) is impaired in MGS during pregnancy and first year of life, but, thereafter, height increases in paralleled normal reference centiles, resulting in a mean adult height of -4.5 standard deviations (SD). Height is dependent on ethnic background and underlying molecular cause, with ORC1 and ORC4 mutations causing more severe short stature and microcephaly. Growth hormone therapy (n = 9) was generally ineffective, though in two patients with significantly reduced IGF1 levels, growth was substantially improved by GH treatment, with 2SD and 3.8 SD improvement in height. Growth parameters for monitoring growth in future MGS patients are provided and as well we highlight that growth is disproportionately affected in certain structures, with growth related minor genital abnormalities (42%) and mammary hypoplasia (100%) frequently present, in addition to established effects on ears and patellar growth.

Genome-wide association analysis of juvenile idiopathic arthritis identifies a new susceptibility locus at chromosomal region 3q13

OBJECTIVE

In a genome-wide association study of Caucasian patients with juvenile idiopathic arthritis (JIA), we have previously described findings limited to autoimmunity loci shared by JIA and other diseases. The present study was undertaken to identify novel JIA-predisposing loci using genome-wide approaches.

METHODS

The discovery cohort consisted of Caucasian JIA cases (n = 814) and local controls (n = 658) genotyped on the Affymetrix Genome-Wide SNP 6.0 Array, along with 2,400 out-of-study controls. In a replication study, we genotyped 10 single-nucleotide polymorphisms (SNPs) in 1,744 cases and 7,010 controls from the US and Europe.

RESULTS

Analysis within the discovery cohort provided evidence of associations at 3q13 within C3orf1 and near CD80 (rs4688011) (odds ratio [OR] 1.37, P = 1.88 × 10(-6) ) and at 10q21 near JMJD1C (rs647989 [OR 1.59, P = 6.1 × 10(-8) ], rs12411988 [OR 1.57, P = 1.16 × 10(-7) ], and rs10995450 [OR 1.31, P = 6.74 × 10(-5) ]). Meta-analysis provided further evidence of association for these 4 SNPs (P = 3.6 × 10(-7) for rs4688011, P = 4.33 × 10(-5) for rs6479891, P = 2.71 × 10(-5) for rs12411988, and P = 5.39 × 10(-5) for rs10995450). Gene expression data on 68 JIA cases and 23 local controls showed cis expression quantitative trait locus associations for C3orf1 SNP rs4688011 (P = 0.024 or P = 0.034, depending on the probe set) and JMJD1C SNPs rs6479891 and rs12411988 (P = 0.01 or P = 0.04, depending on the probe set and P = 0.008, respectively). Using a variance component liability model, it was estimated that common SNP variation accounts for approximately one-third of JIA susceptibility.

CONCLUSION

Genetic association results and correlated gene expression findings provide evidence of JIA association at 3q13 and suggest novel genes as plausible candidates in disease pathology.

Rare and common variants in CARD14, encoding an epidermal regulator of NF-kappaB, in psoriasis

Psoriasis is a common inflammatory disorder of the skin and other organs. We have determined that mutations in CARD14, encoding a nuclear factor of kappa light chain enhancer in B cells (NF-kB) activator within skin epidermis, account for PSORS2. Here, we describe fifteen additional rare missense variants in CARD14, their distribution in seven psoriasis cohorts (>6,000 cases and >4,000 controls), and their effects on NF-kB activation and the transcriptome of keratinocytes. There were more CARD14 rare variants in cases than in controls (burden test p value = 0.0015). Some variants were only seen in a single case, and these included putative pathogenic mutations (c.424G>A [p.Glu142Lys] and c.425A>G [p.Glu142Gly]) and the generalized-pustular-psoriasis mutation, c.413A>C (p.Glu138Ala); these three mutations lie within the coiled-coil domain of CARD14. The c.349G>A (p.Gly117Ser) familial-psoriasis mutation was present at a frequency of 0.0005 in cases of European ancestry. CARD14 variants led to a range of NF-kB activities; in particular, putative pathogenic variants led to levels >2.5× higher than did wild-type CARD14. Two variants (c.511C>A [p.His171Asn] and c.536G>A [p.Arg179His]) required stimulation with tumor necrosis factor alpha (TNF-α) to achieve significant increases in NF-kB levels. Transcriptome profiling of wild-type and variant CARD14 transfectants in keratinocytes differentiated probably pathogenic mutations from neutral variants such as polymorphisms. Over 20 CARD14 polymorphisms were also genotyped, and meta-analysis revealed an association between psoriasis and rs11652075 (c.2458C>T [p.Arg820Trp]; p value = 2.1 × 10(-6)). In the two largest psoriasis cohorts, evidence for association increased when rs11652075 was conditioned on HLA-Cw*0602 (PSORS1). These studies contribute to our understanding of the genetic basis of psoriasis and illustrate the challenges faced in identifying pathogenic variants in common disease.

Hierarchy of risk of childhood-onset rheumatoid arthritis conferred by HLA-DRB1 alleles encoding the shared epitope

OBJECTIVE

Associations between shared epitope (SE)-encoding HLA-DRB1 alleles and rheumatoid arthritis (RA) are well established. However, only a limited number of studies have investigated these alleles in patients with childhood-onset RA, which is defined as rheumatoid factor- and/or anti-citrullinated protein antibody-positive juvenile idiopathic arthritis. The aims of this study were to investigate the largest cohort of patients with childhood-onset RA for association with SE alleles and to determine whether there is a hierarchy of risk based on the amino acid sequence of the SE.

METHODS

High-resolution HLA-DRB1 genotypes were obtained for 204 patients with childhood-onset RA and 373 healthy control subjects. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated for different SE-encoding HLA-DRB1 alleles. In addition, genotype ORs were calculated for combinations of SE alleles classified into S(2) , S(3P) , or L alleles, based on amino acid sequences in position 70-74 of the DRβ1 chain, as proposed by Tezenas du Montcel et al.

RESULTS

We confirmed associations between HLA-DRB1 SE alleles and childhood-onset RA (76% of patients carried 1 or 2 SE alleles compared with 46% of control subjects; OR 3.81, 95% CI 2.4-6.0, P < 1 × 10(-7) ). We also observed associations between individual SE alleles (HLA-DRB1*0101, *0401, *0404, *0405, *0408, and *1001) and childhood-onset RA. Genotype-specific risk estimates suggested a hierarchy of risk, with the highest risk among individuals heterozygous for S(2) /S(3P) (OR 22.3, 95% CI 9.9-50.5, P < 0.0001).

CONCLUSION

We confirm the association between SE-encoding HLA-DRB1 alleles and susceptibility to childhood-onset RA. The excess risk conferred by carriage of the combination of S(2) and S(3P) risk alleles suggests that children with DRβ1 chains containing the KRAA and QRRAA or RRRAA sequences are especially susceptible to RA.

TDT-HET: a new transmission disequilibrium test that incorporates locus heterogeneity into the analysis of family-based association data

BACKGROUND

Locus heterogeneity is one of the most documented phenomena in genetics. To date, relatively little work had been done on the development of methods to address locus heterogeneity in genetic association analysis. Motivated by Zhou and Pan's work, we present a mixture model of linked and unlinked trios and develop a statistical method to estimate the probability that a heterozygous parent transmits the disease allele at a di-allelic locus, and the probability that any trio is in the linked group. The purpose here is the development of a test that extends the classic transmission disequilibrium test (TDT) to one that accounts for locus heterogeneity.

RESULTS

Our simulations suggest that, for sufficiently large sample size (1000 trios) our method has good power to detect association even the proportion of unlinked trios is high (75%). While the median difference (TDT-HET empirical power - TDT empirical power) is approximately 0 for all MOI, there are parameter settings for which the power difference can be substantial. Our multi-locus simulations suggest that our method has good power to detect association as long as the markers are reasonably well-correlated and the genotype relative risk are larger. Results of both single-locus and multi-locus simulations suggest our method maintains the correct type I error rate.Finally, the TDT-HET statistic shows highly significant p-values for most of the idiopathic scoliosis candidate loci, and for some loci, the estimated proportion of unlinked trios approaches or exceeds 50%, suggesting the presence of locus heterogeneity.

CONCLUSIONS

We have developed an extension of the TDT statistic (TDT-HET) that allows for locus heterogeneity among coded trios. Benefits of our method include: estimates of parameters in the presence of heterogeneity, and reasonable power even when the proportion of linked trios is small. Also, we have extended multi-locus methods to TDT-HET and have demonstrated that the empirical power may be high to detect linkage. Last, given that we obtain PPBs, we conjecture that the TDT-HET may be a useful method for correctly identifying linked trios. We anticipate that researchers will find this property increasingly useful as they apply next-generation sequencing data in family based studies.

Common polymorphisms in human lysyl oxidase genes are not associated with the adolescent idiopathic scoliosis phenotype

Background

Although adolescent idiopathic scoliosis affects approximately 3% of adolescents, the genetic contributions have proven difficult to identify. Work in model organisms, including zebrafish, chickens, and mice, has implicated the lysyl oxidase family of enzymes in the development of scoliosis. We hypothesized that common polymorphisms in the five human lysyl oxidase genes (LOX, LOXL1, LOXL2, LOXL3, and LOXL4) may be associated with the phenotype of adolescent idiopathic scoliosis.

Methods

This was a case-control genetic association study. A total of 112 coding and tag SNPs in LOX, LOXL1, LOXL2, LOXL3, and LOXL4 were genotyped in a discovery cohort of 138 cases and 411 controls. Genotypes were tested for association with adolescent idiopathic scoliosis by logistic regression with a two degree of freedom genotypic model and gender as a covariate. Fourteen SNPs with p < 0.1 in the discovery phase were genotyped in an independent replication cohort of 400 cases and 506 controls.

Results

No evidence for significant association was found between coding or tag SNPs in LOX, LOXL1, LOXL2, LOXL3, and LOXL4 and the phenotype of adolescent idiopathic scoliosis.

Conclusions

Despite suggestive evidence in model organisms, common variants and known coding SNPs in the five human lysyl oxidase genes do not confer increased genotypic risk for adolescent idiopathic scoliosis. The above methodology does not address rare variants or individually private mutations in these genes, and future research may focus on this area.

Association of TALS developmental disorder with defect in minor splicing component U4atac snRNA

The spliceosome, a ribonucleoprotein complex that includes proteins and small nuclear RNAs (snRNAs), catalyzes RNA splicing through intron excision and exon ligation to produce mature messenger RNAs, which, in turn serve as templates for protein translation. We identified four point mutations in the U4atac snRNA component of the minor spliceosome in patients with brain and bone malformations and unexplained postnatal death [microcephalic osteodysplastic primordial dwarfism type 1 (MOPD 1) or Taybi-Linder syndrome (TALS); Mendelian Inheritance in Man ID no. 210710]. Expression of a subgroup of genes, possibly linked to the disease phenotype, and minor intron splicing were affected in cell lines derived from TALS patients. Our findings demonstrate a crucial role of the minor spliceosome component U4atac snRNA in early human development and postnatal survival.