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

First copy number variant in trans with single nucleotide variant in CCN6 causing progressive pseudorheumatoid dysplasia revealed by genome sequencing and deep phenotyping in monozygotic twins

Biallelic pathogenic variants in CCN6 cause progressive pseudorheumatoid dysplasia (PPD), a rare skeletal dysplasia. The predominant features include noninflammatory progressive joint stiffness and enlargement, which are not unique to this condition. Nearly 100% of the reported variants are single nucleotide variants or small indels, and missing of a second variant has been reported. Genome sequencing (GS) covers various types of variants and deep phenotyping (DP) provides detailed and precise information facilitating genetic data interpretation. The combination of GS and DP improves diagnostic yield, especially in rare and undiagnosed diseases. We identified a novel compound heterozygote involving a disease-causing copy number variant (g.112057664_112064205del) in trans with a single nucleotide variant (c.624dup(p.Cys209MetfsTer21)) in CCN6 in a pair of monozygotic twins, through the methods of GS and DP. The twins had received three nondiagnostic results before. The g.112057664_112064205del variant was missed by all the tests, and the recorded phenotypes were inaccurate or even misleading. The twins were diagnosed with PPD, ending a 13-year diagnostic odyssey. There may be other patients with PPD experiencing underdiagnosis and misdiagnosis due to inadequate genetic testing or phenotyping methods. This case highlights the critical role of GS and DP in facilitating an accurate and timely diagnosis.

Rare variant association analyses reveal the significant contribution of carbohydrate metabolic disturbance in severe adolescent idiopathic scoliosis

BACKGROUND

Adolescent idiopathic scoliosis (AIS), the predominant genetic-influenced scoliosis, results in spinal deformities without vertebral malformations. However, the molecular aetiology of AIS remains unclear.

METHODS

Using genome/exome sequencing, we studied 368 patients with severe AIS (Cobb angle >40°) and 3794 controls from a Han Chinese cohort. We performed gene-based and pathway-based weighted rare variant association tests to assess the mutational burden of genes and established biological pathways. Differential expression analysis of muscle tissues from 14 patients with AIS and 15 controls was served for validation.

RESULTS

SLC16A8, a lactate transporter linked to retinal glucose metabolism, was identified as a novel severe AIS-associated gene (p=3.08E-06, false discovery rate=0.009). Most AIS cases with deleterious SLC16A8 variants demonstrated early onset high myopia preceding scoliosis. Pathway-based burden test also revealed a significant enrichment in multiple carbohydrate metabolism pathways, especially galactose metabolism. Patients with deleterious variants in these genes demonstrated a significantly larger spinal curve. Genes related to catabolic processes and nutrient response showed divergent expression between AIS cases and controls, reinforcing our genomic findings.

CONCLUSION

This study uncovers the pivotal role of genetic variants in carbohydrate metabolism in the development of AIS, unveiling new insights into its aetiology and potential treatment.

Deciphering the pathogenic role of rare RAF1 heterozygous missense mutation in the late-presenting DDH

Background

Developmental Dysplasia of the Hip (DDH) is a skeletal disorder where late-presenting forms often escape early diagnosis, leading to limb and pain in adults. The genetic basis of DDH is not fully understood despite known genetic predispositions.

Methods

We employed Whole Genome Sequencing (WGS) to explore the genetic factors in late-presenting DDH in two unrelated families, supported by phenotypic analyses and in vitro validation.

Results

In both cases, a novel de novo heterozygous missense mutation in RAF1 (c.193A>G [p.Lys65Glu]) was identified. This mutation impacted RAF1 protein structure and function, altering downstream signaling in the Ras/ERK pathway, as demonstrated by bioinformatics, molecular dynamics simulations, and in vitro validations.

Conclusion

This study contributes to our understanding of the genetic factors involved in DDH by identifying a novel mutation in RAF1. The identification of the RAF1 mutation suggests a possible involvement of the Ras/ERK pathway in the pathogenesis of late-presenting DDH, indicating its potential role in skeletal development.

Core planar cell polarity genes VANGL1 and VANGL2 in predisposition to congenital vertebral malformations

Congenital scoliosis (CS), affecting approximately 0.5 to 1 in 1,000 live births, is commonly caused by congenital vertebral malformations (CVMs) arising from aberrant somitogenesis or somite differentiation. While Wnt/ß-catenin signaling has been implicated in somite development, the function of Wnt/planar cell polarity (Wnt/PCP) signaling in this process remains unclear. Here, we investigated the role of Vangl1 and Vangl2 in vertebral development and found that their deletion causes vertebral anomalies resembling human CVMs. Analysis of exome sequencing data from multiethnic CS patients revealed a number of rare and deleterious variants in VANGL1 and VANGL2, many of which exhibited loss-of-function and dominant-negative effects. Zebrafish models confirmed the pathogenicity of these variants. Furthermore, we found that Vangl1 knock-in (p.R258H) mice exhibited vertebral malformations in a Vangl gene dose- and environment-dependent manner. Our findings highlight critical roles for PCP signaling in vertebral development and predisposition to CVMs in CS patients, providing insights into the molecular mechanisms underlying this disorder.

Somatic mutations in four novel genes contribute to homologous recombination deficiency in breast cancer: a real-world clinical tumor sequencing study

Breast cancers involving mutations in homologous recombination (HR) genes, most commonly BRCA1 and BRCA2 (BRCA1/2), respond well to PARP inhibitors and platinum-based chemotherapy. However, except for these specific HR genes, it is not clear which other mutations contribute to homologous recombination defects (HRD). Here, we performed next-generation sequencing of tumor tissues and matched blood samples from 119 breast cancer patients using the OncoScreen Plus panel. Genomic mutation characteristics and HRD scores were analyzed. In the HR genes, we found that BRCA1/2 and PLAB2 mutations were related to HRD. HRD was also detected in a subset of patients without germline or somatic mutations in BRCA1/2, PLAB2, or other HR-related genes. Notably, LRP1B, NOTCH3, GATA2, and CARD11 (abbreviated as LNGC) mutations were associated with high HRD scores in breast cancer patients. Furthermore, functional experiments demonstrated that silencing CARD11 and GATA2 impairs HR repair efficiency and enhances the sensitivity of tumor cells to olaparib treatment. In summary, in the absence of mutations in the HR genes, the sensitivity of tumor cells to PARP inhibitors and platinum-based chemotherapy may be enhanced in a subset of breast cancer patients with LNGC somatic mutations.

Unraveling the genetic architecture of congenital vertebral malformation with reference to the developing spine

Congenital vertebral malformation, affecting 0.13-0.50 per 1000 live births, has an immense locus heterogeneity and complex genetic architecture. In this study, we analyze exome/genome sequencing data from 873 probands with congenital vertebral malformation and 3794 control individuals. Clinical interpretation identifies Mendelian etiologies in 12.0% of the probands and reveals a muscle-related disease mechanism. Gene-based burden test of ultra-rare variants identifies risk genes with large effect sizes (ITPR2, TBX6, TPO, H6PD, and SEC24B). To further investigate the biological relevance of the genetic association signals, we perform single-nucleus RNAseq on human embryonic spines. The burden test signals are enriched in the notochord at early developmental stages and myoblast/myocytes at late stages, highlighting their critical roles in the developing spine. Our work provides insights into the developmental biology of the human spine and the pathogenesis of spine malformation.

A Study of Polish Family with Scoliosis and Limb Contractures Expands the MYH3 Disease Spectrum

A disease associated with malfunction of the MYH3 gene is characterised by scoliosis, contractures of the V fingers, knees and elbows, dysplasia of the calf muscles, foot deformity and limb length asymmetry. The aim of this study was to identify the cause of musculoskeletal deformities in a three-generation Polish family by exome sequencing. The segregation of the newly described c.866A>C variant of the MYH3 gene in the family indicates an autosomal dominant model of inheritance. The detected MYH3 variant segregates the disease within the family. The presented results expand the MYH3 disease spectrum and emphasize the clinical diagnostic challenge in syndromes harbouring congenital spine defects and joint contractures.

COL11A2 as a candidate gene for vertebral malformations and congenital scoliosis

Human vertebral malformations (VMs) have an estimated incidence of 1/2000 and are associated with significant health problems including congenital scoliosis (CS) and recurrent organ system malformation syndromes such as VACTERL (vertebral anomalies; anal abnormalities; cardiac abnormalities; tracheo-esophageal fistula; renal anomalies; limb anomalies). The genetic cause for the vast majority of VMs are unknown. In a CS/VM patient cohort, three COL11A2 variants (R130W, R1407L and R1413H) were identified in two patients with cervical VM. A third patient with a T9 hemivertebra and the R130W variant was identified from a separate study. These substitutions are predicted to be damaging to protein function, and R130 and R1407 residues are conserved in zebrafish Col11a2. To determine the role for COL11A2 in vertebral development, CRISPR/Cas9 was used to create a nonsense mutation (col11a2L642*) as well as a full gene locus deletion (col11a2del) in zebrafish. Both col11a2L642*/L642* and col11a2del/del mutant zebrafish exhibit vertebral fusions in the caudal spine, which form due to mineralization across intervertebral segments. To determine the functional consequence of VM-associated variants, we assayed their ability to suppress col11a2del VM phenotypes following transgenic expression within the developing spine. While wildtype col11a2 expression suppresses fusions in col11a2del/+ and col11a2del/del backgrounds, patient missense variant-bearing col11a2 failed to rescue the loss-of-function phenotype in these animals. These results highlight an essential role for COL11A2 in vertebral development and support a pathogenic role for two missense variants in CS.

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.

Variants in the SOX9 transactivation middle domain induce axial skeleton dysplasia and scoliosis

SOX9 is an essential transcriptional regulator of cartilage development and homeostasis. In humans, dysregulation of SOX9 is associated with a wide spectrum of skeletal disorders, including campomelic and acampomelic dysplasia, and scoliosis. The mechanism of how SOX9 variants contribute to the spectrum of axial skeletal disorders is not well understood. Here, we report four novel pathogenic variants of SOX9 identified in a large cohort of patients with congenital vertebral malformations. Three of these heterozygous variants are in the HMG and DIM domains, and for the first time, we report a pathogenic variant within the transactivation middle (TAM) domain of SOX9 . Probands with these variants exhibit variable skeletal dysplasia, ranging from isolated vertebral malformation to acampomelic dysplasia. We also generated a Sox9 hypomorphic mutant mouse model bearing a microdeletion within the TAM domain ( Sox9 ^Asp272del ). We demonstrated that disturbance of the TAM domain with missense mutation or microdeletion results in reduced protein stability but does not affect the transcriptional activity of SOX9. Homozygous Sox9 ^Asp272del mice exhibited axial skeletal dysplasia including kinked tails, ribcage anomalies, and scoliosis, recapitulating phenotypes observed in human, while heterozygous mutants display a milder phenotype. Analysis of primary chondrocytes and the intervertebral discs in Sox9 ^Asp272del mutant mice revealed dysregulation of a panel of genes with major contributions of the extracellular matrix, angiogenesis, and ossification-related processes. In summary, our work identified the first pathologic variant of SOX9 within the TAM domain and demonstrated that this variant is associated with reduced SOX9 protein stability. Our finding suggests that reduced SOX9 stability caused by variants in the TAM domain may be responsible for the milder forms of axial skeleton dysplasia in humans.

Retrospective Analysis of Associated Anomalies in 636 Patients with Operatively Treated Congenital Scoliosis

BACKGROUND

Congenital scoliosis is frequently associated with anomalies in multiple organ systems. However, the prevalence and distribution of associated anomalies remain unclear, and there is a large amount of variation in data among different studies.

METHODS

Six hundred and thirty-six Chinese patients who had undergone scoliosis correction surgery at Peking Union Medical College Hospital from January 2012 to July 2019 were recruited, as a part of the Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study. The medical data for each subject were collected and analyzed.

RESULTS

The mean age (and standard deviation) at the time of presentation for scoliosis was 6.4 ± 6.3 years, and the mean Cobb angle of the major curve was 60.8° ± 26.5°. Intraspinal abnormalities were found in 186 (30.3%) of 614 patients, with diastematomyelia being the most common anomaly (59.1%; 110 of 186). The prevalence of intraspinal abnormalities was remarkably higher in patients with failure of segmentation and mixed deformities than in patients with failure of formation (p < 0.001). Patients with intraspinal anomalies showed more severe deformities, including larger Cobb angles of the major curve (p < 0.001). We also demonstrated that cardiac anomalies were associated with remarkably worse pulmonary function, i.e., lower forced expiratory volume in the first second (FEV1), forced vital capacity (FVC), and peak expiratory flow (PEF). Additionally, we identified associations among different concomitant malformations. We found that patients with musculoskeletal anomalies of types other than intraspinal and maxillofacial were 9.2 times more likely to have additional maxillofacial anomalies.

CONCLUSIONS

In our cohort, comorbidities associated with congenital scoliosis occurred at a rate of 55%. To our knowledge, our study is the first to show that patients with congenital scoliosis and cardiac anomalies have reduced pulmonary function, as demonstrated by lower FEV1, FVC, and PEF. Moreover, the potential associations among concomitant anomalies revealed the importance of a comprehensive preoperative evaluation scheme.

LEVEL OF EVIDENCE

Diagnostic Level III. See Instructions for Authors for a complete description of levels of evidence.

Dissection of mendelian predisposition and complex genetic architecture of craniovertebral junction malformation

The craniovertebral junction (CVJ) is an anatomically complex region of the axial skeleton that provides protection of the brainstem and the upper cervical spinal cord. Structural malformation of the CVJ gives rise to life-threatening neurological deficits, such as quadriplegia and dyspnea. Unfortunately, genetic studies on human subjects with CVJ malformation are limited and the pathogenesis remains largely elusive. In this study, we recruited 93 individuals with CVJ malformation and performed exome sequencing. Manual interpretation of the data identified three pathogenic variants in genes associated with Mendelian diseases, including CSNK2A1, MSX2, and DDX3X. In addition, the contribution of copy number variations (CNVs) to CVJ malformation was investigated and three pathogenic CNVs were identified in three affected individuals. To further dissect the complex mutational architecture of CVJ malformation, we performed a gene-based rare variant association analysis utilizing 4371 in-house exomes as control. Rare variants in LGI4 (carrier rate = 3.26%, p = 3.3 × 10^-5) and BEST1 (carrier rate = 5.43%, p = 5.77 × 10^-6) were identified to be associated with CVJ malformation. Furthermore, gene set analyses revealed that extracellular matrix- and RHO GTPase-associated biological pathways were found to be involved in the etiology of CVJ malformation. Overall, we comprehensively dissected the genetic underpinnings of CVJ malformation and identified several novel disease-associated genes and biological pathways.

TBX6 as a cause of a combined skeletal-kidney dysplasia syndrome

TBX6 encodes transcription-factor box 6, a transcription factor critical to paraxial mesoderm segmentation and somitogenesis during embryonic development. TBX6 haploinsufficiency is believed to drive the skeletal and kidney phenotypes associated with the 16p11.2 deletion syndrome. Heterozygous and biallelic variants in TBX6 are associated with vertebral and rib malformations (TBX6-associated congenital scoliosis) and spondylocostal dysostosis, and heterozygous TBX6 variants are associated with increased risk of genitourinary tract malformations. Combined skeletal and kidney phenotypes in individuals harboring heterozygous or biallelic TBX6 variants are rare. Here, we present seven individuals with vertebral and rib malformations and structural kidney differences associated with heterozygous TBX6 gene deletion in trans with a hypomorphic TBX6 allele or biallelic TBX6 variants. Our case series highlights the association between TBX6 and both skeletal and kidney disease.

Phenotype expansion of variants affecting p38 MAPK signaling in hypospadias patients

Background

Hypospadias is a congenital anomaly of the male urogenital system. Genetics factors play an important role in its pathogenesis. To search for potential causal genes/variants for hypospadias, we performed exome sequencing in a pedigree with three patients across two generations and a cohort of 49 sporadic patients with hypospadias.

Results

A novel BRAF variant (NM_004333.6: c.362C > A) was found to co-segregate with the hypospadias phenotype in the disease pedigree. In cells overexpressing the BRAF mutant, the phosphorylation level of p38 MAPK was significantly increased as compared with the cells overexpressing the wild-type BRAF or RASopathy-related BRAF mutant. This variant further led to a reduced transcription level of the SRY gene, which is essential for the normal development of the male reproductive system. In the cohort of sporadic patients, we identified two additional variants in p38 MAPK signaling-related genes (TRIM67 and DAB2IP) potentially associated with hypospadias.

Conclusion

Our study expands the phenotypic spectrum of variants affecting p38 MAPK signaling toward the involvement of hypospadias.

Exploring the association between specific genes and the onset of idiopathic scoliosis: a systematic review

Background

Idiopathic Scoliosis (IS) is the most common spinal deformity in adolescents, accounting for 80% of all spinal deformities. However, the etiology remains uncertain in most cases, being identified as Adolescent Idiopathic Scoliosis (AIS). IS treatments range from observation and sport to bracing or surgery. Several risk factors including sex and familiarity, have been linked with IS. Although there are still many uncertainties regarding the cause of this pathology, several studies report a greater incidence of the defect in families in which at least one other first degree relative is affected. This study systematically reviews the available literature to identify the most significant genes or variants related to the development and onset of IS.

Methods

The research question was formulated using a PIOS approach on the following databases: Medline, Embase, Cinahl, Scopus, Web of Science and Google Scholar. The search was performed from July to August 2021, and articles from the inception of the database to August 2021 were searched.

Results

24 of the 919 initially identified studies were included in the present review. The 24 included studies observed a total of 16,316 cases and 81,567 controls. All the considered studies stated either the affected gene and/or specific SNPs. CHD7, SH2B1, ESR, CALM1, LBX1, MATN1, CHL1, FBN1 and FBN2 genes were associated with IS development.

Conclusions

Although association can be found in some candidate genes the field of research regarding genetic association with the onset of IS still requires more information.

The utility of hierarchical genetic testing in paediatric liver disease

BACKGROUND & AIMS

Genetic factors underlie a substantial proportion of paediatric liver diseases. Hereditary liver diseases have considerable genetic heterogeneity and variable clinical manifestations, which bring great challenges to clinical and molecular diagnoses. In this study, we investigated a group of paediatric patients with varying degrees of liver dysfunction using a hierarchical genetic testing strategy.

METHODS

We first applied a panel encompassing 166 known causal genes of liver disease. We then used exome sequencing (ES) in those patients whose cases remained undiagnosed to identify the genetic aetiology of their symptoms.

RESULTS

In total, we enrolled 131 unrelated paediatric patients with liver disease of Chinese Han ethnicity. We first applied targeted gene sequencing of 166 genes to all patients and yielded a diagnostic rate of 35.9% (47 of 131). Eighty-four patients who remained undiagnosed after target gene sequencing were subjected to ES. As a result, eight (8/84, 9.5%) of them obtained molecular diagnoses, including four patients suspected of abnormal bilirubin metabolism and four idiopathic cases. Non-typical genetic findings, including digenic inheritance and dual molecular diagnosis, were also identified. Through a comprehensive assessment of novel candidate variants of uncertain disease association, 11 patients of the remaining undiagnosed patients were able to obtain likely molecular diagnoses.

CONCLUSIONS

Our study presents evidence for the diagnostic utility of sequential genetic testing in a cohort of patients with paediatric liver disease. Our findings expand the understanding of the phenotypic and mutational spectrum underlying this heterogeneous group of diseases.

Understanding Musculoskeletal Disorders Through Next-Generation Sequencing

»

An insight into musculoskeletal disorders through advancements in next-generation sequencing (NGS) promises to maximize benefits and improve outcomes through improved genetic diagnosis.

»

The primary use of whole exome sequencing (WES) for musculoskeletal disorders is to identify functionally relevant variants.

»

The current evidence has shown the superiority of NGS over conventional genotyping for identifying novel and rare genetic variants in patients with musculoskeletal disorders, due to its high throughput and low cost.

»

Genes identified in patients with scoliosis, osteoporosis, osteoarthritis, and osteogenesis imperfecta using NGS technologies are listed for further reference.

Delineation of dual molecular diagnosis in patients with skeletal deformity

Background

Skeletal deformity is characterized by an abnormal anatomical structure of bone and cartilage. In our previous studies, we have found that a substantial proportion of patients with skeletal deformity could be explained by monogenic disorders. More recently, complex phenotypes caused by more than one genetic defect (i.e., dual molecular diagnosis) have also been reported in skeletal deformities and may complicate the diagnostic odyssey of patients. In this study, we report the molecular and phenotypic characteristics of patients with dual molecular diagnosis and variable skeletal deformities.

Results

From 1108 patients who underwent exome sequencing, we identified eight probands with dual molecular diagnosis and variable skeletal deformities. All eight patients had dual diagnosis consisting of two autosomal dominant diseases. A total of 16 variants in 12 genes were identified, 5 of which were of de novo origin. Patients with dual molecular diagnosis presented blended phenotypes of two genetic diseases. Mendelian disorders occurred more than once include Osteogenesis Imperfecta Type I (COL1A1, MIM:166200), Neurofibromatosis, Type I (NF1, MIM:162200) and Marfan Syndrome (FBN1, MIM:154700).

Conclusions

This study demonstrated the complicated skeletal phenotypes associated with dual molecular diagnosis. Exome sequencing represents a powerful tool to detect such complex conditions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02293-x.

Identification of Novel FBN2 Variants in a Cohort of Congenital Contractural Arachnodactyly

Congenital contractural arachnodactyly (CCA) is a rare autosomal dominant disorder of connective tissue characterized by crumpled ears, arachnodactyly, camptodactyly, large joint contracture, and kyphoscoliosis. The nature course of CCA has not been well-described. We aim to decipher the genetic and phenotypic spectrum of CCA. The cohort was enrolled in Beijing Jishuitan Hospital and Peking Union Medical College Hospital, Beijing, China, based on Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study (http://www.discostudy.org/). Exome sequencing was performed on patients’ blood DNA. A recent published CCA scoring system was validated in our cohort. Seven novel variants and three previously reported FBN2 variants were identified through exome sequencing. Two variants outside of the neonatal region of FBN2 gene were found. The phenotypes were comparable between patients in our cohort and previous literature, with arachnodactyly, camptodactyly and large joints contractures found in almost all patients. All patients eligible for analysis were successfully classified into likely CCA based on the CCA scoring system. Furthermore, we found a double disease-causing heterozygous variant of FBN2 and ANKRD11 in a patient with blended phenotypes consisting of CCA and KBG syndrome. The identification of seven novel variants broadens the mutational and phenotypic spectrum of CCA and may provide implications for genetic counseling and clinical management.

Exome-wide Analysis of De Novo and Rare Genetic Variants in Patients With Brain Arteriovenous Malformation

BACKGROUND AND OBJECTIVES

Brain arteriovenous malformation (bAVM) is a congenital disorder and a leading cause of hemorrhagic stroke. Germline genetic variants play an essential role in the pathogenesis of brain arteriovenous malformation. However, the biological relevance of the disease-associated genes identified in previous studies is elusive. In this study, we aim to systematically investigate the contribution of germline variants to bAVM and explore the critical molecular pathways underlying the pathogenesis of bAVM.

METHODS

Probands with sporadic bAVM were consecutively recruited into this study from November 2015 to November 2018 and underwent exome sequencing. The controls were aggregated from individuals who were not known to have vascular malformation and underwent exome sequencing for clinical or research purposes. The retained control dataset included 4609 individuals, including 251 individuals with parental samples sequenced. We firstly compared de novo variants in cases and controls and performed a pathway enrichment analysis. A gene-based rare variant association analysis was then performed to identify genes whose variants were significantly enriched in cases.

RESULTS

We collected an exome-sequenced bAVM cohort consisting of 152 trios and 40 singletons. By firstly focusing on de novo variants, we observed a significant mutational burden of de novo likely gene-disrupting variants in cases versus controls. By performing a pathway enrichment analysis of all nonsynonymous de novo variants identified in cases, we found the angiopoietin-like protein 8 (ANGPTL8) regulatory pathway to be significantly enriched in patients with bAVM. Through an exome-wide rare variant association analysis utilizing 4394 in-house exome data as controls, we identified SLC19A3 as a disease-associated gene for bAVM. In addition, we found that the SLC19A3 variants in cases are preferably located at the N' side of the SLC19A3 protein. These findings implicate a phenotypic extension of SLC19A3-related disorders with a domain-specific effect.

DISCUSSION

This study provides insights into the biological basis of bAVM by identifying novel molecular pathways and candidate genes.

DrABC: deep learning accurately predicts germline pathogenic mutation status in breast cancer patients based on phenotype data

Background

Identifying breast cancer patients with DNA repair pathway-related germline pathogenic variants (GPVs) is important for effectively employing systemic treatment strategies and risk-reducing interventions. However, current criteria and risk prediction models for prioritizing genetic testing among breast cancer patients do not meet the demands of clinical practice due to insufficient accuracy.

Methods

The study population comprised 3041 breast cancer patients enrolled from seven hospitals between October 2017 and 11 August 2019, who underwent germline genetic testing of 50 cancer predisposition genes (CPGs). Associations among GPVs in different CPGs and endophenotypes were evaluated using a case-control analysis. A phenotype-based GPV risk prediction model named DNA-repair Associated Breast Cancer (DrABC) was developed based on hierarchical neural network architecture and validated in an independent multicenter cohort. The predictive performance of DrABC was compared with currently used models including BRCAPRO, BOADICEA, Myriad, PENN II, and the NCCN criteria.

Results

In total, 332 (11.3%) patients harbored GPVs in CPGs, including 134 (4.6%) in BRCA2, 131 (4.5%) in BRCA1, 33 (1.1%) in PALB2, and 37 (1.3%) in other CPGs. GPVs in CPGs were associated with distinct endophenotypes including the age at diagnosis, cancer history, family cancer history, and pathological characteristics. We developed a DrABC model to predict the risk of GPV carrier status in BRCA1/2 and other important CPGs. In predicting GPVs in BRCA1/2, the performance of DrABC (AUC = 0.79 [95% CI, 0.74–0.85], sensitivity = 82.1%, specificity = 63.1% in the independent validation cohort) was better than that of previous models (AUC range = 0.57–0.70). In predicting GPVs in any CPG, DrABC (AUC = 0.74 [95% CI, 0.69–0.79], sensitivity = 83.8%, specificity = 51.3% in the independent validation cohort) was also superior to previous models in their current versions (AUC range = 0.55–0.65). After training these previous models with the Chinese-specific dataset, DrABC still outperformed all other methods except for BOADICEA, which was the only previous model with the inclusion of pathological features. The DrABC model also showed higher sensitivity and specificity than the NCCN criteria in the multi-center validation cohort (83.8% and 51.3% vs. 78.8% and 31.2%, respectively, in predicting GPVs in any CPG). The DrABC model implementation is available online at http://gifts.bio-data.cn/.

Conclusions

By considering the distinct endophenotypes associated with different CPGs in breast cancer patients, a phenotype-driven prediction model based on hierarchical neural network architecture was created for identification of hereditary breast cancer. The model achieved superior performance in identifying GPV carriers among Chinese breast cancer patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-022-01027-9.

Mutational spectrum of syndromic genes in sporadic brain arteriovenous malformation

Background

Brain arteriovenous malformations (BAVMs) are abnormal vessels that are apt to rupture, causing life-threatening intracranial hemorrhage (ICH). The estimated prevalence of BAVMs is 0.05% among otherwise healthy individuals. In this study, we aim to investigate the mutational spectrum of syndromic genes in sporadic BAVM.

Methods

We recruited a cohort of 150 patients with BAVM and performed whole-exome sequencing on their peripheral blood DNA. To explore the mutational spectrum of syndromic genes in sporadic brain arteriovenous malformation, we selected six genes according to the Online Mendelian Inheritance in Man (OMIM) and literature. All variants in the six candidate genes were extracted and underwent filtering for qualifying variants.

Results

There are a total of four patients with rare variants in hereditary hemorrhagic telangiectasia-related genes. In addition, we identified two patients have the variant of RASA1 gene in our database, which are also rare mutations that are absent from population databases. However, we did not find any patients with GNAQ mutations in our database.

Conclusions

In conclusion, we demonstrated that variants in syndromic vascular malformations play important roles in the etiology of sporadic BAVM.

Expanding the mutation and phenotype spectrum of MYH3-associated skeletal disorders

Pathogenic variants in MYH3 cause distal arthrogryposis type 2A and type 2B3 as well as contractures, pterygia and spondylocarpotarsal fusion syndromes types 1A and 1B. These disorders are ultra-rare and their natural course and phenotypic variability are not well described. In this study, we summarize the clinical features and genetic findings of 17 patients from 10 unrelated families with vertebral malformations caused by dominant or recessive pathogenic variants in MYH3. Twelve novel pathogenic variants in MYH3 (NM_002470.4) were identified: three of them were de novo or inherited in autosomal dominant way and nine were inherited in autosomal recessive way. The patients had vertebral segmentation anomalies accompanied with variable joint contractures, short stature and dysmorphic facial features. There was a significant phenotypic overlap between dominant and recessive MYH3-associated conditions regarding the degree of short stature as well as the number of vertebral fusions. All monoallelic variants caused significantly decreased SMAD3 phosphorylation, which is consistent with the previously proposed pathogenic mechanism of impaired canonical TGF-β signaling. Most of the biallelic variants were predicted to be protein-truncating, while one missense variant c.4244T>G,p.(Leu1415Arg), which was inherited in an autosomal recessive way, was found to alter the phosphorylation level of p38, suggesting an inhibition of the non-canonical pathway of TGF-β signaling. In conclusion, the identification of 12 novel pathogenic variants and overlapping phenotypes in 17 affected individuals from 10 unrelated families expands the mutation and phenotype spectrum of MYH3-associated skeletal disorders. We show that disturbances of canonical or non-canonical TGF-β signaling pathways are involved in pathogenesis of MYH3-associated skeletal fusion (MASF) syndrome.

PhenoApt leverages clinical expertise to prioritize candidate genes via machine learning

In recent years, exome sequencing (ES) has shown great utility in the diagnoses of Mendelian disorders. However, after rigorous filtering, a typical ES analysis still involves the interpretation of hundreds of variants, which greatly hinders the rapid identification of causative genes. Since the interpretations of ES data require comprehensive clinical analyses, taking clinical expertise into consideration can speed the molecular diagnoses of Mendelian disorders. To leverage clinical expertise to prioritize candidate genes, we developed PhenoApt, a phenotype-driven gene prioritization tool that allows users to assign a customized weight to each phenotype, via a machine-learning algorithm. Using the ability to rank causative genes in top-10 lists as an evaluation metric, baseline analysis demonstrated that PhenoApt outperformed previous phenotype-driven gene prioritization tools by a relative increase of 22.7%-140.0% in three independent, real-world, multi-center cohorts (cohort 1, n = 185; cohort 2, n = 784; and cohort 3, n = 208). Additional trials showed that, by adding weights to clinical indications, which should be explained by the causative gene, PhenoApt performance was improved by a relative increase of 37.3% in cohort 2 (n = 471) and 21.4% in cohort 3 (n = 208). Moreover, PhenoApt could assign an intrinsic weight to each phenotype based on the likelihood of its being a Mendelian trait using term frequency-inverse document frequency techniques. When clinical indications were assigned with intrinsic weights, PhenoApt performance was improved by a relative increase of 23.7% in cohort 2 and 15.5% in cohort 3. For the integration of PhenoApt into clinical practice, we developed a user-friendly website and a command-line tool.

A genotype-first analysis in a cohort of Mullerian anomaly

Müllerian anomaly (M.A.) is a group of congenital anatomic abnormalities caused by aberrations of the development process of the Müllerian duct. M.A. can either be isolated or be involved in Mendelian syndromes, such as Dandy-Walker syndrome, Holt-Oram syndrome and Bardet-Biedl syndrome, which are often associated with both uterus and kidney malformations. In this study, we applied a genotype-first approach to analyze the whole-exome sequencing data of 492 patients with M.A. Six potential pathogenic variants were found in five genes previously related to female urogenital deformities (PKD1, SON, SALL1, BMPR1B, ITGA8), which are partially overlapping with our patients' phenotypes. We further identified eight incidental findings in seven genes related to Mendelian syndromes without known association with reproductive anomalies (TEK, COL11A1, ANKRD11, LEMD3, DLG5, SPTB, BMP2), which represent potential phenotype expansions of these genes.

Aberrant interaction between mutated ADAMTSL2 and LTBP4 is associated with adolescent idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a complex spinal structure deformity with a prevalence of 1%-3%. Genetic and hereditary factors have been associated with the etiology of AIS. However, previous studies mainly focused on common single nucleotide polymorphisms which confer modest disease risk. Recently, rare variants in FBN1 and other extracellular matrix genes have been implicated in AIS, suggesting a potential overlapping disease etiology between AIS and hereditary connective tissue disorders (HCTD). In this study, we systematically analyzed rare variants in a set of HCTD-related genes in 302 AIS patients using exome sequencing. We firstly focused on pathogenic variants based on a monogenic inheritance and identified nine disease-associated variants in FBN1, COL11A1, COL11A2 and TGFBR2. We then explored the potential interactions between variants in different genes based on the case-control statistics. We identified three ADAMTSL2-LTBP4 variant pairs in three AIS patients and none in controls. Furthermore, we revealed that the variant pairs identified in these genes could affect the interaction between ADAMTSL2 and LTBP4 and upregulate TGF-β signaling pathway in human fibroblasts. Our findings supported that the aberrant interaction between mutated ADAMTSL2 and LTBP4 was associated with AIS.

The Mutational Landscape of PTK7 in Congenital Scoliosis and Adolescent Idiopathic Scoliosis

Depletion of ptk7 is associated with both congenital scoliosis (CS) and adolescent idiopathic scoliosis (AIS) in zebrafish models. However, only one human variant of PTK7 has been reported previously in a patient with AIS. In this study, we systemically investigated the variant landscape of PTK7 in 583 patients with CS and 302 patients with AIS from the Deciphering Disorders Involving Scoliosis and COmorbidities (DISCO) study. We identified a total of four rare variants in CS and four variants in AIS, including one protein truncating variant (c.464_465delAC) in a patient with CS. We then explored the effects of these variants on protein expression and sub-cellular location. We confirmed that the c.464_465delAC variant causes loss-of-function (LoF) of PTK7. In addition, the c.353C>T and c.2290G>A variants identified in two patients with AIS led to reduced protein expression of PTK7 as compared to that of the wild type. In conclusion, LoF and hypomorphic variants are associated with CS and AIS, respectively.

Disruptive NADSYN1 Variants Implicated in Congenital Vertebral Malformations

Genetic perturbations in nicotinamide adenine dinucleotide de novo (NAD) synthesis pathway predispose individuals to congenital birth defects. The NADSYN1 encodes the final enzyme in the de novo NAD synthesis pathway and, therefore, plays an important role in NAD metabolism and organ embryogenesis. Biallelic mutations in the NADSYN1 gene have been reported to be causative of congenital organ defects known as VCRL syndrome (Vertebral-Cardiac-Renal-Limb syndrome). Here, we analyzed the genetic variants in NADSYN1 in an exome-sequenced cohort consisting of patients with congenital vertebral malformations (CVMs). A total number of eight variants in NADSYN1, including two truncating variants and six missense variants, were identified in nine unrelated patients. All enrolled patients presented multiple organ defects, with the involvement of either the heart, kidney, limbs, or liver, as well as intraspinal deformities. An in vitro assay using COS-7 cells demonstrated either significantly reduced protein levels or disrupted enzymatic activity of the identified variants. Our findings demonstrated that functional variants in NADSYN1 were involved in the complex genetic etiology of CVMs and provided further evidence for the causative NADSYN1 variants in congenital NAD Deficiency Disorder.

Novel FGFR1 Variants Are Associated with Congenital Scoliosis

FGFR1 encodes a transmembrane cytokine receptor, which is involved in the early development of the human embryo and plays an important role in gastrulation, organ specification and patterning of various tissues. Pathogenic FGFR1 variants have been associated with Kallmann syndrome and hypogonadotropic hypogonadism. In our congenital scoliosis (CS) patient series of 424 sporadic CS patients under the framework of the Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study, we identified four unrelated patients harboring FGFR1 variants, including one frameshift and three missense variants. These variants were predicted to be deleterious by in silico prediction and conservation analysis. Signaling activities and expression levels of the mutated protein were evaluated in vitro and compared to that of the wild type (WT) FGFR1. As a result, the overall protein expressions of c.2334dupC, c.2339T>C and c.1261A>G were reduced to 43.9%, 63.4% and 77.4%, respectively. By the reporter gene assay, we observed significantly reduced activity for c.2334dupC, c.2339T>C and c.1261A>G, indicating the diminished FGFR1 signaling pathway. In conclusion, FGFR1 variants identified in our patients led to only mild disruption to protein function, caused milder skeletal and cardiac phenotypes than those reported previously.

The Distinct Performances of Ultrasound, Mammograms, and MRI in Detecting Breast Cancer in Patients With Germline Pathogenic Variants in Cancer Predisposition Genes

A proportion of up to 10% of breast cancer resulted from hereditary germline pathogenic variants (GPVs) in cancer predisposition genes (CPGs), which been demonstrated distinct clinical features and imaging manifestations. However, the performance of imaging modalities for breast cancer surveillance in CPG mutation-carriers is still unclear, especially in Asian women. A population of 3002 breast cancer patients who received germline genetic testing of CPGs was enrolled from three hospitals in China. In total, 343 (11.6%) patients were found to harbor GPVs in CPGs, including 137 (4.6%) in BRCA1 and 135 (4.6%) in BRCA2. We compared the performances of ultrasound, mammograms, MRI, and the combining strategies in CPG mutation carriers and non-carriers. As a result, the ultrasound showed a higher detection rate compared with mammograms regardless of the mutation status. However, its detection rate was lower in CPG mutation carriers than in non-carriers (93.2% vs 98.0%, P=2.1×10-4), especially in the BRCA1 mutation carriers (90.9% vs 98.0%, P=2.0×10-4). MRI presented the highest sensitivity (98.5%) and the lowest underestimation rate (14.5%) in CPG mutation carriers among ultrasound, mammograms, and their combination. Supplemental ultrasound or mammograms would add no significant value to MRI for detecting breast cancer (P>0.05). In multivariate logistic regression analysis, the family or personal cancer history could not replace the mutation status as the impact factor for the false-negative result and underestimation. In summary, clinicians and radiologists should be aware of the atypical imaging presentation of breast cancer in patients with GPVs in CPGs.

Deciphering the mutational signature of congenital limb malformations

Congenital limb malformations (CLMs) affect 1 in 500 live births. However, the value of exome sequencing (ES) for CLM is lacking. The purpose of this study was to decipher the mutational signature of CLM on an exome level. We enrolled a cohort of 66 unrelated probands (including 47 families) with CLM requiring surgical correction. ES was performed for all patients and available parental samples. A definite molecular diagnosis was achieved in 21 out of 66 (32%) patients. We identified 19 pathogenic or likely pathogenic single-nucleotide variants and three copy number variants, of which 11 variants were novel. We identified four variants of uncertain significance. Additionally, we identified RPL9 and UBA2 as novel candidate genes for CLM. By comparing the detailed phenotypic features, we expand the phenotypic spectrum of diastrophic dysplasia and chromosome 6q terminal deletion syndrome. We also found that the diagnostic rate was significantly higher in patients with a family history of CLM (p = 0.012) or more than one limb affected (p = 0.034). Our study expands our understanding of the mutational and phenotypic spectrum of CLM and provides novel insights into the genetic basis of these syndromes.

Heterozygous Recurrent Mutations Inducing Dysfunction of ROR2 Gene in Patients With Short Stature

PURPOSE

ROR2, a member of the ROR family, is essential for skeletal development as a receptor of Wnt5a. The present study aims to investigate the mutational spectrum of ROR2 in children with short stature and to identify the underlying molecular mechanisms.

METHODS

We retrospectively analyzed clinical phenotype and whole-exome sequencing (WES) data of 426 patients with short stature through mutation screening of ROR2. We subsequently examined the changes in protein expression and subcellular location in ROR2 caused by the mutations. The mRNA expression of downstream signaling molecules of the Wnt5a-ROR2 pathway was also examined.

RESULTS

We identified 12 mutations in ROR2 in 21 patients, including 10 missense, one nonsense, and one frameshift. Among all missense variants, four recurrent missense variants [c.1675G > A(p.Gly559Ser), c.2212C > T(p.Arg738Cys), c.1930G > A(p.Asp644Asn), c.2117G > A(p.Arg706Gln)] were analyzed by experiments in vitro. The c.1675G > A mutation significantly altered the expression and the cellular localization of the ROR2 protein. The c.1675G > A mutation also caused a significantly decreased expression of c-Jun. In contrast, other missense variants did not confer any disruptive effect on the biological functions of ROR2.

CONCLUSION

We expanded the mutational spectrum of ROR2 in patients with short stature. Functional experiments potentially revealed a novel molecular mechanism that the c.1675G > A mutation in ROR2 might affect the expression of downstream Wnt5a-ROR2 pathway gene by disturbing the subcellular localization and expression of the protein.

Exome sequencing reveals genetic architecture in patients with isolated or syndromic short stature

Short stature is among the most common endocrinological disease phenotypes of childhood and may occur as an isolated finding or in conjunction with other clinical manifestations. Although the diagnostic utility of clinical genetic testing in short stature has been implicated, the genetic architecture and the utility of genomic studies such as exome sequencing (ES) in a sizable cohort of patients with short stature have not been investigated systematically. In this study, we recruited 561 individuals with short stature from two centers in China during a 4-year period. We performed ES for all patients and available parents. All patients were retrospectively divided into two groups: an isolated short stature group (group I, n = 257) and an apparently syndromic short stature group (group II, n = 304). Causal variants were identified in 135 of 561 (24.1%) patients. In group I, 29 of 257 (11.3%) of the patients were solved by variants in 24 genes. In group II, 106 of 304 (34.9%) patients were solved by variants in 57 genes. Genes involved in fundamental cellular process played an important role in the genetic architecture of syndromic short stature. Distinct genetic architectures and pathophysiological processes underlie isolated and syndromic short stature.

Variants Affecting the C-Terminal of CSF1R Cause Congenital Vertebral Malformation Through a Gain-of-Function Mechanism

CSF1R encodes the colony-stimulating factor 1 receptor which regulates the proliferation, differentiation, and biological activity of monocyte/macrophage lineages. Pathogenic variants in CSF1R could lead to autosomal dominant adult-onset leukoencephalopathy with axonal spheroids and pigmented glia or autosomal recessive skeletal dysplasia. In this study, we identified three heterozygous deleterious rare variants in CSF1R from a congenital vertebral malformation (CVM) cohort. All of the three variants are located within the carboxy-terminal region of CSF1R protein and could lead to an increased stability of the protein. Therefore, we established a zebrafish model overexpressing CSF1R. The zebrafish model exhibits CVM phenotypes such as hemivertebral and vertebral fusion. Furthermore, overexpression of the mutated CSF1R mRNA depleted of the carboxy-terminus led to a higher proportion of zebrafish with vertebral malformations than wild-type CSF1R mRNA did (p = 0.03452), implicating a gain-of-function effect of the C-terminal variant. In conclusion, variants affecting the C-terminal of CSF1R could cause CVM though a potential gain-of-function mechanism.

Whole-genome methylation analysis reveals novel epigenetic perturbations of congenital scoliosis

Congenital scoliosis (CS) is a congenital disease caused by malformations of vertebrae. Recent studies demonstrated that DNA modification could contribute to the pathogenesis of disease. This study aims to identify epigenetic perturbations that may contribute to the pathogenesis of CS. Four CS patients with hemivertebra were enrolled and underwent spine correction operations. DNA was extracted from the hemivertebrae and spinal process collected from the specimen during the hemivertebra resection. Genome-wide DNA methylation profiling was examined at base-pair resolution using whole-genome bisulfite sequencing (WGBS). We identified 343 genes with hyper-differentially methylated regions (DMRs) and 222 genes with hypo-DMRs, respectively. These genes were enriched in the mitogen-activated protein kinase (MAPK) signaling pathway, calcium signaling pathway, and axon guidance in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and were enriched in positive regulation of cell morphogenesis involved in differentiation, regulation of cell morphogenesis involved in differentiation, and regulation of neuron projection development in Biological Process of Gene Ontology (GO-BP) terms. Hyper-DMR-related genes, including IGHG1, IGHM, IGHG3, RNF213, and GSE1, and hypo DMR-related genes, including SORCS2, COL5A1, GRID1, RGS3, and ROBO2, may contribute to the pathogenesis of hemivertebra. The aberrant DNA methylation may be associated with the formation of hemivertebra and congenital scoliosis.

Whole Exome Sequencing Uncovered the Genetic Architecture of Growth Hormone Deficiency Patients

PURPOSE

Congenital growth hormone deficiency (GHD) is a rare and etiologically heterogeneous disease. We aim to screen disease-causing mutations of GHD in a relatively sizable cohort and discover underlying mechanisms via a candidate gene-based mutational burden analysis.

METHODS

We retrospectively analyzed 109 short stature patients associated with hormone deficiency. All patients were classified into two groups: Group I (n=45) with definitive GHD and Group II (n=64) with possible GHD. We analyzed correlation consistency between clinical criteria and molecular findings by whole exome sequencing (WES) in two groups. The patients without a molecular diagnosis (n=90) were compared with 942 in-house controls for the mutational burden of rare mutations in 259 genes biologically related with the GH axis.

RESULTS

In 19 patients with molecular diagnosis, we found 5 possible GHD patients received known molecular diagnosis associated with GHD (NF1 [c.2329T>A, c.7131C>G], GHRHR [c.731G>A], STAT5B [c.1102delC], HRAS [c.187_207dup]). By mutational burden analysis of predicted deleterious variants in 90 patients without molecular diagnosis, we found that POLR3A (p = 0.005), SUFU (p = 0.006), LHX3 (p = 0.021) and CREB3L4 (p = 0.040) represented top genes enriched in GHD patients.

CONCLUSION

Our study revealed the discrepancies between the laboratory testing and molecular diagnosis of GHD. These differences should be considered when for an accurate diagnosis of GHD. We also identified four candidate genes that might be associated with GHD.

Cost-effectiveness analysis of using the TBX6-associated congenital scoliosis risk score (TACScore) in genetic diagnosis of congenital scoliosis

Background

We previously reported a novel clinically distinguishable subtype of congenital scoliosis (CS), namely, TBX6-associated congenital scoliosis (TACS). We further developed the TBX6-associated CS risk score (TACScore), a multivariate phenotype-based model to predict TACS according to the patient’s clinical manifestations. In this study, we aimed to evaluate whether using the TACScore as a screening method prior to performing whole-exome sequencing (WES) is more cost-effective than using WES as the first-line genetic test for CS.

Methods

We retrospectively collected the molecular data of 416 CS patients in the Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study. A decision tree was constructed to estimate the cost and the diagnostic time required for the two alternative strategies (TACScore versus WES). Bootstrapping simulations and sensitivity analyses were performed to examine the distributions and robustness of the estimates. The economic evaluation considered both the health care payer and the personal budget perspectives.

Results

From the health care payer perspective, the strategy of using the TACScore as the primary screening method resulted in an average cost of $1074.2 (95%CI: $1044.8 to $1103.5) and an average diagnostic duration of 38.7d (95%CI: 37.8d to 39.6d) to obtain a molecular diagnosis for each patient. In contrast, the corresponding values were $1169.6 (95%CI: $1166.9 to $1172.2) and 41.4d (95%CI: 41.1d to 41.7d) taking WES as the first-line test (P < 0.001). From the personal budget perspective, patients who were predicted to be positive by the TACScore received a result with an average cost of $715.1 (95%CI: $594.5 to $835.7) and an average diagnostic duration of 30.4d (95%CI: 26.3d to 34.6d). Comparatively, the strategy of WES as the first-line test was estimated to have significantly longer diagnostic time with an average of 44.0d (95%CI: 43.2d to 44.9d), and more expensive with an average of $1193.4 (95%CI: $1185.5 to $1201.3) (P < 0.001). In 100% of the bootstrapping simulations, the TACScore strategy was significantly less costly and more time-saving than WES. The sensitivity analyses revealed that the TACScore strategy remained cost-effective even when the cost per WES decreased to $8.8.

Conclusions

This retrospective study provides clinicians with economic evidence to integrate the TACScore into clinical practice. The TACScore can be considered a cost-effective tool when it serves as a screening test prior to performing WES.

Mutational burden and potential oligogenic model of TBX6-mediated genes in congenital scoliosis

BACKGROUND

Congenital scoliosis (CS) is a spinal deformity due to vertebral malformations. Although insufficiency of TBX6 dosage contributes to a substantial proportion of CS, the molecular etiology for the majority of CS remains largely unknown. TBX6-mediated genes involved in the process of somitogenesis represent promising candidates.

METHODS

Individuals affected with CS and without a positive genetic finding were referred to this study. Proband-only exome sequencing (ES) were performed on the recruited individuals, followed by analysis of TBX6-mediated candidate genes, namely MEOX1, MEOX2, MESP2, MYOD1, MYF5, RIPPLY1, and RIPPLY2.

RESULTS

A total of 584 patients with CS of unknown molecular etiology were recruited. After ES analysis, protein-truncating variants in RIPPLY1 and MYF5 were identified from two individuals, respectively. In addition, we identified five deleterious missense variants (MYOD1, n = 4; RIPPLY2, n = 1) in TBX6-mediated genes. We observed a significant mutational burden of MYOD1 in CS (p = 0.032) compared with the in-house controls (n = 1854). Moreover, a potential oligogenic disease-causing mode was proposed based on the observed mutational co-existence of MYOD1/MEOX1 and MYOD1/RIPPLY1.

CONCLUSION

Our study characterized the mutational spectrum of TBX6-mediated genes, prioritized core candidate genes/variants, and provided insight into a potential oligogenic disease-causing mode in CS.