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

Concurrent of compound heterozygous variant of a novel in-frame deletion and the common hypomorphic haplotype in TBX6 and inherited 17q12 microdeletion in a fetus

BACKGROUND

TBX6, a member of the T-box gene family, encodes the transcription factor box 6 that is critical for somite segmentation in vertebrates. It is known that the compound heterozygosity of disruptive variants in trans with a common hypomorphic risk haplotype (T-C-A) in the TBX6 gene contribute to 10% of congenital scoliosis (CS) cases. The deletion of chromosome 17q12 is a rare cytogenetic abnormality, which often leads to renal cysts and diabetes mellitus. However, the affected individuals often exhibit clinical heterogeneity and incomplete penetrance.

METHODS

We here present a Chinese fetus who was shown to have CS by ultrasound examination at 17 weeks of gestation. Trio whole-exome sequencing (WES) was performed to investigate the underlying genetic defects of the fetus. In vitro functional experiments, including western-blotting and luciferase transactivation assay, were performed to determine the pathogenicity of the novel variant of TBX6.

RESULTS

WES revealed the fetus harbored a compound heterozygous variant of c.338_340del (p.Ile113del) and the common hypomorphic risk haplotype of the TBX6 gene. In vitro functional study showed the p.Ile113del variant had no impact on TBX6 expression, but almost led to complete loss of its transcriptional activity. In addition, we identified a 1.85 Mb deletion on 17q12 region in the fetus and the mother. Though there is currently no clinical phenotype associated with this copy number variation in the fetus, it can explain multiple renal cysts in the pregnant woman.

CONCLUSIONS

This study is the first to report a Chinese fetus with a single amino acid deletion variant and a T-C-A haplotype of TBX6. The clinical heterogeneity of 17q12 microdeletion poses significant challenges for prenatal genetic counseling. Our results once again suggest the complexity of prenatal genetic diagnosis.

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.

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.

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.

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.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.

Genetics of non-isolated hemivertebra: A systematic review of fetal, neonatal, and infant cases

Hemivertebra is a congenital vertebral malformation caused by unilateral failure of formation during embryogenesis that may be associated with additional abnormalities. A systematic review was conducted to investigate genetic etiologies of non-isolated hemivertebra identified in the fetal, neonatal, and infant periods using PubMed, Cochrane database, Ovid Medline, and ClinicalTrials.gov from inception through May 2022 (PROSPERO ID CRD42021229576). The Human Phenotype Ontology database was accessed May 2022. Studies were deemed eligible for inclusion if they addressed non-isolated hemivertebra or genetic causes of non-isolated hemivertebra identified in the fetal, neonatal, or infant periods. Cases diagnosed clinically without molecular confirmation were included. Systematic review identified 23 cases of non-isolated hemivertebra with karyotypic abnormalities, 2 cases due to microdeletions, 59 cases attributed to single gene disorders, 18 syndromic cases without known genetic etiology, and 14 cases without a known syndromic association. The Human Phenotype Ontology search identified 49 genes associated with hemivertebra. Non-isolated hemivertebra is associated with a diverse spectrum of cytogenetic abnormalities and single gene disorders. Genetic syndromes were notably common. Frequently affected organ systems include musculoskeletal, cardiovascular, central nervous system, genitourinary, gastrointestinal, and facial dysmorphisms. When non-isolated hemivertebra is identified on prenatal ultrasound, the fetus must be assessed for associated anomalies and genetic counseling is recommended.

Molecular identification of T-box transcription factor 6 and prognostic assessment in patients with congenital scoliosis: A single-center study

Background

Congenital scoliosis (CS) is characterized by vertebral malformations. The precise etiology of CS is not fully defined. A compound inheritance of TBX6 was identified in 10% of patients with CS in Han Chinese and formed a distinguishable subtype named TBX6-associated congenital scoliosis (TACS).

Methods

To investigate the variants and risk haplotype of TBX6, we recruited 121 patients with CS at Beijing Children’s Hospital. We collected the clinical characteristics and surgical treatment options and followed their postoperative prognoses.

Results

Eight patients (6.6%) were molecularly diagnosed with TACS and carried the previously defined pathogenic TBX6 compound heterozygous variants. All the eight patients with TACS had the typical TACS clinical feature of hemivertebrae in the lower part of the spine. These patients received posterior hemivertebra resection combined with segmental fusion. Follow-ups revealed satisfactory correction without postoperative complications.

Conclusion

We observed a 6.6% prevalence of TACS in our CS cohort. Follow-ups further highlighted that surgical treatment of hemivertebra resection combined with segmental fusion performed well with prognosis for patients with TACS. This could provide valuable information for CS individuals with compound heterozygosity in TBX6.

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.

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.

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.

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.

Cross-sectional analysis and trend of vertebral and associated anomalies in Chinese congenital scoliosis population: a retrospective study of one thousand, two hundred and eighty nine surgical cases from 2010 to 2019

PURPOSE

The objective of this article is to report associated anomaly incidences of a large CS cohort and analyze interrelationships among vertebral anomaly types and associated abnormalities.

METHODS

We retrospectively searched and extracted medical records of 1289 CS inpatients surgically treated in our institute from January 2010 to December 2019. All patients have taken spine X-ray, CT, MRI, echocardiogram, urogenital ultrasound, and systemic physical examination. We analyzed information on demographics, CS types, and associated anomalies.

RESULTS

CS type was found to be 49.1% for failure of formation (FF), 19.5% for failure of segmentation (FS), and 31.4% for mixed defects (MD). Intraspinal defects were found in 29.4% patients (16.0% for FF, 45.4% for FS, 40.5% for MD), cardiac in 13.7% (12.3% for FF, 14.3% for FS, 15.6% for MD), genitourinary in 5.8% (4.1% for FF, 6.0% for FS, 8.4% for MD), gastrointestinal in 3.6% (4.7% for FF, 1.6% for FS, 3.0% for MD), and musculoskeletal in 16.4% (10.3% for FF, 19.9% for FS, 23.7% for MD). The intraspinal and musculoskeletal defect incidences were significantly higher in patients with failure of segmentation and mixed defects. We also observed a decreasing trend for intraspinal and musculoskeletal defect incidences as well as a tendency for more failure of formation and less failure of segmentation from 2010 to 2019.

CONCLUSIONS

The intraspinal and musculoskeletal defect incidences were higher in patients with failure of segmentation and mixed defects. Strong interrelationships were found between intraspinal and musculoskeletal defects and among cardiovascular, genitourinary, and gastrointestinal defects. From 2010 to 2019, the proportion of patients with failure of formation increased significantly, causing a decrease in the intraspinal and musculoskeletal defect incidences over time. Female sex, failure of segmentation, and mixed defects could be considered risk factors for more associated anomalies in CS individuals, which would help surgeons in medical management and prenatal consultation.

Genetic variants of TBX6 and TBXT identified in patients with congenital scoliosis in Southern China

Congenital scoliosis (CS) is a spinal deformity present at birth due to underlying congenital vertebral malformation (CVM) that occurs during embryonic development. Hemivertebrae is the most common anomaly that causes CS. Recently, compound heterozygosity in TBX6 has been identified in Northern Chinese, Japanese, and European CS patient cohorts, which explains about 7%-10% of the affected population. In this report, we recruited 67 CS patients characterized with hemivertebrae in the Southern Chinese population and investigated the TBX6 variant and risk haplotype. We found that two patients with hemivertebrae in the thoracic spine and one patient with hemivertebrae in the lumbar spine carry the previously defined pathogenic TBX6 compound heterozygous variants. In addition, whole exome sequencing of patients with CS and their family members identified a de novo missense mutation (c.G47T: p.R16L) in another member of the T-box family, TBXT. This rare mutation compromised the binding of TBXT to its target sequence, leading to reduced transcriptional activity, and exhibited dominant-negative effect on wild-type TBXT. Our findings further highlight the importance of T-box family genes in the development of congenital scoliosis.

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.

Low penetrance COL5A1 variants in a young patient with intracranial aneurysm and very mild signs of Ehlers-Danlos syndrome

Spontaneous cervical artery dissection (CeAD) is a major cause of ischemic stroke in young adults, whose genetic susceptibility factors are still largely unknown. Nevertheless, subtle ultrastructural connective tissue alterations (especially in the collagen fibril morphology) are recognized in a large proportion of CeAD patients, in which recent genetic investigations reported an enrichment of variants in genes associated with known connective tissue disorders. In this regard, COL5A1 variants have been reported in a small subset of CeAD patients, with or without classical Ehlers-Danlos syndrome (cEDS) features. We investigated a 22-year-old patient with intracranial aneurysm and mild connective tissue manifestations reminiscent of EDS. Whole-exome sequencing identified two COL5A1 missense variants in trans configuration: NM_000093.5:c.[1588G>A];[4135C>T], NP_000084.3:p.[(Gly530Ser)];[(Pro1379Ser)]. Functional assays demonstrated a significant decrease of collagen α1(V) chain expression in both heterozygous parents compared to control cells, and an additive effect of these two variants in the proband. Interestingly, both parents manifested very subtle EDS signs, such as atrophic scars, recurrent bone fractures, colonic diverticulosis, varicose veins, and osteoarthritis. Our findings emphasize the involvement of COL5A1 in the predisposition to vascular phenotypes and provide novel insights on the c.1588G>A variant, whose functional significance has not been definitely established. In fact, it was previously reported as both "disease modifying", and as a biallelic causative mutation (with heterozygous individuals showing subtle clinical signs of cEDS). We speculated that the c.1588G>A variant might lead to overt phenotype in combination with additional genetic "hits" lowering the collagen α1(V) chain expression below a hypothetical disease threshold.

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.

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.

KIAA1217: A novel candidate gene associated with isolated and syndromic vertebral malformations

Vertebral malformations (VMs) are caused by alterations in somitogenesis and may occur in association with other congenital anomalies. The genetic etiology of most VMs remains unknown and their identification may facilitate the development of novel therapeutic and prevention strategies. Exome sequencing was performed on both the discovery cohort of nine unrelated probands from the USA with VMs and the replication cohort from China (Deciphering Disorders Involving Scoliosis & COmorbidities study). The discovery cohort was analyzed using the PhenoDB analysis tool. Heterozygous and homozygous, rare and functional variants were selected and evaluated for their ClinVar, HGMD, OMIM, GWAS, mouse model phenotypes, and other annotations to identify the best candidates. Genes with candidate variants in three or more probands were selected. The replication cohort was analyzed by another in-house developed pipeline. We identified rare heterozygous variants in KIAA1217 in four out of nine probands in the discovery cohort and in five out of 35 probands in the replication cohort. Collectively, we identified 11 KIAA1217 rare variants in 10 probands, three of which have not been described in gnomAD and one of which is a nonsense variant. We propose that genetic variations of KIAA1217 may contribute to the etiology of VMs.