Compound Heterozygosity for Null Mutations and a Common Hypomorphic Risk Haplotype in TBX6 Causes Congenital Scoliosis

Butterfly Vertebra: A Retrospective Study of 30 Patients

BACKGROUND

Butterfly vertebra (BV) is a rare congenital spinal anomaly for which there is a paucity of large-scale retrospective studies and established guidelines for treatment. The objective of this study was to elucidate the clinical characteristics, imaging findings, and therapeutic approaches for BV.

METHODS

We conducted a retrospective analysis of 30 patients diagnosed with BV at our hospital from 2009 to 2023, examining clinical data, imaging findings, and clinical interventions.

RESULTS

The analysis comprised a cohort of 30 patients, consisting of 15 males and 15 females, with a mean age of 27.63 ± 19.84 years. Imaging studies indicated that the majority of vertebral bodies affected by BV were single-segmented (63.3%, 19/30) and less commonly multi-segmented (36.7%, 11/30). These findings frequently coexisted with other medical conditions, most notably spinal scoliosis (76.7%, 23/30). Furthermore, the study identified a range of spinal abnormalities among patients, including hemivertebral deformity (30.0%, 9/30), spinal cleft (10.0%, 3/30), lumbar disc protrusion or herniation (10.0%, 3/30), vertebral slippage (10.0%, 3/30), thoracic kyphosis deformity (6.67%, 2/30), vertebral fusion deformity (6.67%, 2/30), compressive fractures (3.3%, 1/30), and vertebral developmental anomalies (3.3%, 1/30). Clinical intervention resulted in symptom relief for 23 nonsurgical patients through lifestyle modifications, analgesic use, and physical therapy. Seven surgical patients underwent appropriate surgical procedures, leading to satisfaction and adherence to regular postoperative follow-up appointments.

CONCLUSIONS

BV is a rare vertebral anomaly that can be easily misdiagnosed due to its similarity to other diseases. Consequently, it is imperative to enhance vigilance in the differential diagnosis process in order to promptly recognize BV. Furthermore, in cases where patients present with additional associated radiographic findings, a thorough evaluation is typically warranted and timely measures should be taken for treatment.

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.

Molecular landscape of congenital vertebral malformations: recent discoveries and future directions

Vertebral malformations (VMs) pose a significant global health problem, causing chronic pain and disability. Vertebral defects occur as isolated conditions or within the spectrum of various congenital disorders, such as Klippel-Feil syndrome, congenital scoliosis, spondylocostal dysostosis, sacral agenesis, and neural tube defects. Although both genetic abnormalities and environmental factors can contribute to abnormal vertebral development, our knowledge on molecular mechanisms of numerous VMs is still limited. Furthermore, there is a lack of resource that consolidates the current knowledge in this field. In this pioneering review, we provide a comprehensive analysis of the latest research on the molecular basis of VMs and the association of the VMs-related causative genes with bone developmental signaling pathways. Our study identifies 118 genes linked to VMs, with 98 genes involved in biological pathways crucial for the formation of the vertebral column. Overall, the review summarizes the current knowledge on VM genetics, and provides new insights into potential involvement of biological pathways in VM pathogenesis. We also present an overview of available data regarding the role of epigenetic and environmental factors in VMs. We identify areas where knowledge is lacking, such as precise molecular mechanisms in which specific genes contribute to the development of VMs. Finally, we propose future research avenues that could address knowledge gaps.

The role of Notch signaling pathway in metabolic bone diseases

Metabolic bone diseases is the third most common endocrine diseases after diabetes and thyroid diseases. More than 500 million people worldwide suffer from metabolic bone diseases. The generation and development of bone metabolic diseases is a complex process regulated by multiple signaling pathways, among which the Notch signaling pathway is one of the most important pathways. The Notch signaling pathway regulates the differentiation and function of osteoblasts and osteoclasts, and affects the process of cartilage formation, bone formation and bone resorption. Genetic mutations in upstream and downstream of Notch signaling genes can lead to a series of metabolic bone diseases, such as Alagille syndrome, Adams-Oliver syndrome and spondylocostal dysostosis. In this review, we analyzed the mechanisms of Notch ligands, Notch receptors and signaling molecules in the process of signal transduction, and summarized the progress on the pathogenesis and clinical manifestations of bone metabolic diseases caused by Notch gene mutation. We hope to draw attention to the role of the Notch signaling pathway in metabolic bone diseases and provide new ideas and approaches for the diagnosis and treatment of metabolic bone diseases.

Prenatally diagnosed 16p11.2 copy number variations by SNP Array: A retrospective case series

OBJECTIVE

The 16p11.2 copy number variations (CNVs) are increasingly recognized as one of the most frequent genomic disorders, with a broad spectrum of phenotypes. The fetal phenotype associated with 16p11.2 CNVs is poorly described. The current study presents prenatal series of 16p11.2 CNVs and provides a better understanding of this submicroscopic imbalance in prenatal diagnosis.

METHOD

Retrospective case series were extracted from a single tertiary referral center performing prenatal single nucleotide polymorphism (SNP) array from April 2017 to December 2021. The maternal demographics, indication for amniocentesis, ultrasound findings, SNP array results, inheritance of the CNVs, and pregnancy outcomes were studied.

RESULTS

We indentified 30 fetuses carrying 16p11.2 CNVs, representing 0.35% (30/8578) of prenatal SNP array results. The series included 17 fetuses with a proximal deletion, 7 with a distal deletion, 4 with a proximal duplication, and 2 with a distal duplication. Prenatal ultrasound anomalies were reported in 80% of these cases. The most common presentation was vertebralanomalies (9/30). Other features noted in more than one fetus were increased nuchal translucency/nuchal fold (NT/NF) (5/30), absent/hypoplastic nasal bone (3/30), polyhydramnios (3/30), ventricular septal defect (VSD) (2/30), unilateral mild ventriculomegaly (2/30), fetal growth restriction (FGR) (2/30), right aortic arch (2/30). All the 9 vertebralanomalies were present in fetuses harboring proximal deletion (9/17). Familial transmission was confirmed in 44% of cases (11/25) and termination of pregnancy was requested in 62.1% (18/29) of cases.

CONCLUSION

The 16p11.2 CNVs can have variable prenatal phenotypes and these CNVs are frequently inherited from parents with a milder or normal phenotype. Our results underline that vertebral deformities were frequent in cases of 16p11.2 proximal deletion, and further demonstrate the incomplete penetrance of the CNVs.

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.

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.

Association of LBX1 Gene Methylation Level with Disease Severity in Patients with Idiopathic Scoliosis: Study on Deep Paravertebral Muscles

Idiopathic scoliosis (IS) is a multifactorial disease with a genetic background. The association of Ladybird Homeobox 1 (LBX1) polymorphisms with IS has been proven in multiple studies. However, the epigenetic mechanisms have not been evaluated. This study aimed to evaluate the LBX1 methylation level in deep paravertebral muscles in order to analyze its association with IS occurrence and/or IS severity. Fifty-seven IS patients and twenty non-IS patients were examined for the paravertebral muscles’ methylation level of the LBX1 promoter region. There was no significant difference in methylation level within paravertebral muscles between patients vs. controls, except for one CpG site. The comparison of the paravertebral muscles’ LBX1 promoter region methylation level between patients with a major curve angle of ≤70° vs. >70° revealed significantly higher methylation levels in 17 of 23 analyzed CpG sequences at the convex side of the curvature in patients with a major curve angle of >70° for the reverse strand promoter region. The association between LBX1 promoter methylation and IS severity was demonstrated. In patients with severe IS, the deep paravertebral muscles show an asymmetric LBX1 promoter region methylation level, higher at the convex scoliosis side, which reveals the role of locally acting factors in IS progression.

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.

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.

Identification of Copy Number Variants in a Southern Chinese Cohort of Patients with Congenital Scoliosis

Congenital scoliosis (CS) is a lateral curvature of the spine resulting from congenital vertebral malformations (CVMs) and affects 0.5–1/1000 live births. The copy number variant (CNV) at chromosome 16p11.2 has been implicated in CVMs and recent studies identified a compound heterozygosity of 16p11.2 microdeletion and TBX6 variant/haplotype causing CS in multiple cohorts, which explains about 5–10% of the affected cases. Here, we studied the genetic etiology of CS by analyzing CNVs in a cohort of 67 patients with congenital hemivertebrae and 125 family controls. We employed both candidate gene and family-based approaches to filter CNVs called from whole exome sequencing data. This identified 12 CNVs in four scoliosis-associated genes (TBX6, NOTCH2, DSCAM, and SNTG1) as well as eight recessive and 64 novel rare CNVs in 15 additional genes. Some candidates, such as DHX40, NBPF20, RASA2, and MYSM1, have been found to be associated with syndromes with scoliosis or implicated in bone/spine development. In particular, the MYSM1 mutant mouse showed spinal deformities. Our findings suggest that, in addition to the 16p11.2 microdeletion, other CNVs are potentially important in predisposing to CS.

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.

Exome Sequencing Reveals De Novo Variants in Congenital Scoliosis

Congenital scoliosis (CS) is a lateral curvature of the spine characterized by the presence of vertebral anomalies. Pathogenic genetic variants in the TBX6 gene are one of the causes of CS. However, since many clinically diagnosed cases of CS are without known TBX6 gene variations, this study aims to uncover new genes related to disease susceptibility of CS by exome sequencing (ES). This study employed ES in a cohort of 5 Japanese patients with CS and their healthy parents or a sister for a total of 16 samples among 5 families. Variant interpretation was performed using SIFT, PolyPhen-2, Mutation Taster, and CADD. Four de novo variants were identified by ES and confirmed by Sanger sequencing: 1 frameshift variant (SHISA3) and 3 missense variants (AGBL5, HDAC4, and PDE2A). ES also uncovered 1 homozygous variant in the MOCOS gene. All of these variants were predicted to be deleterious by SIFT, PolyPhen-2, Mutation Taster, and/or CADD. The number of de novo variants identified in this study was exactly what would be expected by chance. Additional functional studies or gathering matched patients using Gene Matcher are needed.

A Survey of Compound Heterozygous Variants in Pediatric Cancers and Structural Birth Defects

Compound heterozygous (CH) variants occur when two recessive alleles are inherited and the variants are located at different loci within the same gene in a given individual. CH variants are important contributors to many different types of recessively inherited diseases. However, many studies overlook CH variants because identification of this type of variant requires knowing the parent of origin for each nucleotide. Using computational methods, haplotypes can be inferred using a process called "phasing," which estimates the chromosomal origin of most nucleotides. In this paper, we used germline, phased, whole-genome sequencing (WGS) data to identify CH variants across seven pediatric diseases (adolescent idiopathic scoliosis: n = 16, congenital heart defects: n = 709, disorders of sex development: n = 79, ewing sarcoma: n = 287, neuroblastoma: n = 259, orofacial cleft: n = 107, and syndromic cranial dysinnervation: n = 172), available as parent-child trios in the Gabriella Miller Kids First Data Resource Center. Relatively little is understood about the genetic underpinnings of these diseases. We classified CH variants as "potentially damaging" based on minor allele frequencies (MAF), Combined Annotation Dependent Depletion scores, variant impact on transcription or translation, and gene-level frequencies in the disease group compared to a healthy population. For comparison, we also identified homozygous alternate (HA) variants, which affect both gene copies at a single locus; HA variants represent an alternative mechanism of recessive disease development and do not require phasing. Across all diseases, 2.6% of the samples had a potentially damaging CH variant and 16.2% had a potentially damaging HA variant. Of these samples with potentially damaging variants, the average number of genes per sample was 1 with a CH variant and 1.25 with a HA variant. Across all samples, 5.1 genes per disease had a CH variant, while 35.6 genes per disease had a HA variant; on average, only 4.3% of these variants affected common genes. Therefore, when seeking to identify potentially damaging variants of a putatively recessive disease, CH variants should be considered as potential contributors to disease development. If CH variants are excluded from analysis, important candidate genes may be overlooked.

Altered Cogs of the Clock: Insights into the Embryonic Etiology of Spondylocostal Dysostosis

Spondylocostal dysostosis (SCDO) is a rare heritable congenital condition, characterized by multiple severe malformations of the vertebrae and ribs. Great advances were made in the last decades at the clinical level, by identifying the genetic mutations underlying the different forms of the disease. These were matched by extraordinary findings in the Developmental Biology field, which elucidated the cellular and molecular mechanisms involved in embryo body segmentation into the precursors of the axial skeleton. Of particular relevance was the discovery of the somitogenesis molecular clock that controls the progression of somite boundary formation over time. An overview of these concepts is presented, including the evidence obtained from animal models on the embryonic origins of the mutant-dependent disease. Evidence of an environmental contribution to the severity of the disease is discussed. Finally, a brief reference is made to emerging in vitro models of human somitogenesis which are being employed to model the molecular and cellular events occurring in SCDO. These represent great promise for understanding this and other human diseases and for the development of more efficient therapeutic approaches.

Increased TBX6 gene dosages induce congenital cervical vertebral malformations in humans and mice

BACKGROUND

Congenital vertebral malformations (CVMs) manifest with abnormal vertebral morphology. Genetic factors have been implicated in CVM pathogenesis, but the underlying pathogenic mechanisms remain unclear in most subjects. We previously reported that the human 16p11.2 BP4-BP5 deletion and its associated TBX6 dosage reduction caused CVMs. We aim to investigate the reciprocal 16p11.2 BP4-BP5 duplication and its potential genetic contributions to CVMs.

METHODS AND RESULTS

Patients who were found to carry the 16p11.2 BP4-BP5 duplication by chromosomal microarray analysis were retrospectively analysed for their vertebral phenotypes. The spinal assessments in seven duplication carriers showed that four (57%) presented characteristics of CVMs, supporting the contention that increased TBX6 dosage could induce CVMs. For further in vivo functional investigation in a model organism, we conducted genome editing of the upstream regulatory region of mouse Tbx6 using CRISPR-Cas9 and obtained three mouse mutant alleles (Tbx6up1 to Tbx6up3 ) with elevated expression levels of Tbx6. Luciferase reporter assays showed that the Tbx6up3 allele presented with the 160% expression level of that observed in the reference (+) allele. Therefore, the homozygous Tbx6up3/up3 mice could functionally mimic the TBX6 dosage of heterozygous carriers of 16p11.2 BP4-BP5 duplication (approximately 150%, ie, 3/2 gene dosage of the normal level). Remarkably, 60% of the Tbx6up3/up3 mice manifested with CVMs. Consistent with our observations in humans, the CVMs induced by increased Tbx6 dosage in mice mainly affected the cervical vertebrae.

CONCLUSION

Our findings in humans and mice consistently support that an increased TBX6 dosage contributes to the risk of developing cervical CVMs.

Human and mouse studies establish TBX6 in Mendelian CAKUT and as a potential driver of kidney defects associated with the 16p11.2 microdeletion syndrome

Congenital anomalies of the kidney and urinary tract (CAKUTs) are the most common cause of chronic kidney disease in children. Human 16p11.2 deletions have been associated with CAKUT, but the responsible molecular mechanism remains to be illuminated. To explore this, we investigated 102 carriers of 16p11.2 deletion from multi-center cohorts, among which we retrospectively ascertained kidney morphologic and functional data from 37 individuals (12 Chinese and 25 Caucasian/Hispanic). Significantly higher CAKUT rates were observed in 16p11.2 deletion carriers (about 25% in Chinese and 16% in Caucasian/Hispanic) than those found in the non-clinically ascertained general populations (about 1/1000 found at autopsy). Furthermore, we identified seven additional individuals with heterozygous loss-of-function variants in TBX6, a gene that maps to the 16p11.2 region. Four of these seven cases showed obvious CAKUT. To further investigate the role of TBX6 in kidney development, we engineered mice with mutated Tbx6 alleles. The Tbx6 heterozygous null (i.e., loss-of-function) mutant (Tbx6^+/‒) resulted in 13% solitary kidneys. Remarkably, this incidence increased to 29% in a compound heterozygous model (Tbx6^mh/‒) that reduced Tbx6 gene dosage to below haploinsufficiency, by combining the null allele with a novel mild hypomorphic allele (mh). Renal hypoplasia was also frequently observed in these Tbx6-mutated mouse models. Thus, our findings in patients and mice establish TBX6 as a novel gene involved in CAKUT and its gene dosage insufficiency as a potential driver for kidney defects observed in the 16p11.2 microdeletion syndrome.

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.

Retrospective analysis of fetal vertebral defects: Associated anomalies, etiologies, and outcome

Our objectives were to describe fetal cases of vertebral defects (VD), assess the diagnostic yield of fetal chromosomal analysis for VD and determine which investigations should be performed when evaluating fetal VD. We performed a retrospective chart review for fetuses with VD seen between 2006 and 2015. Cases were identified from CHU Sainte-Justine's prenatal clinic visits, postmortem fetal skeletal surveys, and medical records. Cases with neural tube defects were excluded. Sixty-six fetuses with VD were identified at a mean gestational age of 20 weeks. Forty-seven (71.2%) had associated antenatal anomalies, most commonly genitourinary, skeletal/limb, and cardiac anomalies. Thirteen mothers (19.7%) had pregestational diabetes (95% CI [10.1%-29.3%]). Fifty-three cases had chromosomal analysis. Three had abnormal results (5.6%): trisomy 13, trisomy 22, and 9q33.1q34.11 deletion. Thirty-four (51.5%) pregnancies were terminated, one led to intrauterine fetal demise and 31 (46.9%) continued to term. Of 27 children who survived the neonatal period, 21 had congenital scoliosis and 3 had spondylocostal dysostosis. Seven had developmental delay. In conclusion, prenatal evaluation of fetuses with VD should include detailed morphological assessment (including fetal echocardiogram), maternal diabetes screening, and chromosomal microarray if non-isolated. Our findings provide guidance about management and counseling after a diagnosis of fetal VD.

Whole-Exome Sequencing Identified a TBX6 Loss of Function Mutation in a Patient with Distal Vaginal Atresia

STUDY OBJECTIVE

The purpose of this study was to determine if there are any genetic changes with whole-exome sequencing associated with distal vaginal atresia.

DESIGN

This was a retrospective genetics study of 5 patients who presented with distal vaginal atresia who were recruited between 2017 and 2018. Whole-exome sequencing was performed in each subject with distal vaginal atresia. Sanger sequencing was used to confirm the potential causative genetic mutation.

SETTING

Beijing Obstetrics and Gynecology Hospital, Capital Medical University, Beijing, China.

PARTICIPANTS AND MAIN OUTCOME MEASURES

The main outcome measure was the rare mutations potentially associated with distal vaginal atresia in 5 patients.

RESULTS

A truncating mutation c.266delC (p.P89Rfs*5) in the T-box transcription factor 6 (TBX6) gene, which is highly expressed in the human vagina, was identified in 1 patient using whole-exome sequencing. The deletion of the 16p11.2 region containing the TBX6 locus has also been reported previously to have the clinical feature of Müllerian agenesis. This mutation was paternally inherited by the patient. This truncating mutation was absent from all of the databases we checked, suggesting that the variant is rare and pathogenic.

CONCLUSION

We showed, to our knowledge, for the first time, that the mutation in TBX6 might be associated with human distal vaginal atresia.

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

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

Key Regulators of Cardiovascular Differentiation and Regeneration: Harnessing the Potential of Direct Reprogramming to Treat Heart Failure

Cardiovascular diseases remain a leading cause of death worldwide, with the number of patients with heart failure increasing rapidly in aging societies. As adult cardiomyocytes are terminally differentiated cells and opportunities for heart transplantation are very limited, regenerative medicine may become a game changer in heart failure treatment. To develop strategies for generating cardiomyocytes, vascular cells, and other supporting cells, it is necessary to understand the mechanism of cardiovascular differentiation during development and from pluripotent stem cells. Master regulators for cardiovascular differentiation can generate new cardiomyocytes and vascular cells directly from other differentiated cells such as fibroblasts. Fibroblasts can be directly reprogrammed into cardiomyocytes by overexpressing a combination of 3 cardiac-specific transcription factors (Gata4, Mef2c, Tbx5) both in vitro and in vivo, which restores cardiac function after myocardial infarction in mice. Moreover, a direct reprogramming-based approach can be used to identify new key regulators of the cardiovascular mesoderm, which can differentiate into all 3 types of cardiovascular cells including cardiomyocytes, endothelial cells, and smooth muscle cells. This review provides a perspective on how key regulators for cardiovascular differentiation and regeneration can be identified and used to develop new treatments for heart failure.

TBX6 compound inheritance leads to congenital vertebral malformations in humans and mice

Congenital vertebral malformations (CVMs) are associated with human TBX6 compound inheritance that combines a rare null allele and a common hypomorphic allele at the TBX6 locus. Our previous in vitro evidence suggested that this compound inheritance resulted in a TBX6 gene dosage of less than haploinsufficiency (i.e. <50%) as a potential mechanism of TBX6-associated CVMs. To further investigate this pathogenetic model, we ascertained and collected 108 Chinese CVM cases and found that 10 (9.3%) of them carried TBX6 null mutations in combination with common hypomorphic variants at the second TBX6 allele. For in vivo functional verification and genetic analysis of TBX6 compound inheritance, we generated both null and hypomorphic mutations in mouse Tbx6 using the CRISPR-Cas9 method. These Tbx6 mutants are not identical to the patient variants at the DNA sequence level, but instead functionally mimic disease-associated TBX6 variants. Intriguingly, as anticipated by the compound inheritance model, a high penetrance of CVM phenotype was only observed in the mice with combined null and hypomorphic alleles of Tbx6. These findings are consistent with our experimental observations in humans and supported the dosage effect of TBX6 in CVM etiology. In conclusion, our findings in the newly collected human CVM subjects and Tbx6 mouse models consistently support the contention that TBX6 compound inheritance causes CVMs, potentially via a gene dosage-dependent mechanism. Furthermore, mouse Tbx6 mutants mimicking human CVM-associated variants will be useful models for further mechanistic investigations of CVM pathogenesis in the cases associated with TBX6.

Bi-allelic loss of function variants of TBX6 causes a spectrum of malformation of spine and rib including congenital scoliosis and spondylocostal dysostosis

Background

Congenital scoliosis (CS) is a common vertebral malformation. Spondylocostal dysostosis (SCD) is a rare skeletal dysplasia characterised by multiple vertebral malformations and rib anomalies. In a previous study, a compound heterozygosity for a null mutation and a risk haplotype composed by three single-nucleotide polymorphisms in TBX6 have been reported as a disease-causing model of CS. Another study identified bi-allelic missense variants in a SCD patient. The purpose of our study is to identify TBX6 variants in CS and SCD and examine their pathogenicity.

Methods

We recruited 200 patients with CS or SCD and investigated TBX6 variants. We evaluated the pathogenicity of the variants by in silico prediction and in vitro experiments.

Results

We identified five 16p11.2 deletions, one splice-site variant and five missense variants in 10 patients. In vitro functional assays for missense variants identified in the previous and present studies demonstrated that most of the variants caused abnormal localisation of TBX6 proteins. We confirmed mislocalisation of TBX6 proteins in presomitic mesoderm cells induced from SCD patient-derived iPS cells. In induced cells, we found decreased mRNA expressions of TBX6 and its downstream genes were involved in somite formation. All CS patients with missense variants had the risk haplotype in the opposite allele, while a SCD patient with bi-allelic missense variants did not have the haplotype.

Conclusions

Our study suggests that bi-allelic loss of function variants of TBX6 cause a spectrum of phenotypes including CS and SCD, depending on the severity of the loss of TBX6 function.

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

INTRODUCTION

Adolescent idiopathic scoliosis (AIS) is a common musculoskeletal disorder with strong evidence for a genetic contribution. CNVs play an important role in congenital scoliosis, but their role in idiopathic scoliosis has been largely unexplored.

METHODS

Exome sequence data from 1197 AIS cases and 1664 in-house controls was analysed using coverage data to identify rare CNVs. CNV calls were filtered to include only highly confident CNVs with >10 average reads per region and mean log-ratio of coverage consistent with single-copy duplication or deletion. The frequency of 55 common recurrent CNVs was determined and correlated with clinical characteristics.

RESULTS

Distal chromosome 16p11.2 microduplications containing the gene SH2B1 were found in 0.7% of AIS cases (8/1197). We replicated this finding in two additional AIS cohorts (8/1097 and 2/433), resulting in 0.7% (18/2727) of all AIS cases harbouring a chromosome 16p11.2 microduplication, compared with 0.06% of local controls (1/1664) and 0.04% of published controls (8/19584) (p=2.28×10^-11, OR=16.15). Furthermore, examination of electronic health records of 92 455 patients from the Geisinger health system showed scoliosis in 30% (20/66) patients with chromosome 16p11.2 microduplications containing SH2B1 compared with 7.6% (10/132) of controls (p=5.6×10^-4, OR=3.9).

CONCLUSIONS

Recurrent distal chromosome 16p11.2 duplications explain nearly 1% of AIS. Distal chromosome 16p11.2 duplications may contribute to scoliosis pathogenesis by directly impairing growth or by altering expression of nearby genes, such as TBX6. Individuals with distal chromosome 16p11.2 microduplications should be screened for scoliosis to facilitate early treatment.

TBX6-associated congenital scoliosis (TACS) as a clinically distinguishable subtype of congenital scoliosis: further evidence supporting the compound inheritance and TBX6 gene dosage model

PURPOSE

To characterize clinically measurable endophenotypes, implicating the TBX6 compound inheritance model.

METHODS

Patients with congenital scoliosis (CS) from China(N = 345, cohort 1), Japan (N = 142, cohort 2), and the United States (N = 10, cohort 3) were studied. Clinically measurable endophenotypes were compared according to the TBX6 genotypes. A mouse model for Tbx6 compound inheritance (N = 52) was investigated by micro computed tomography (micro-CT). A clinical diagnostic algorithm (TACScore) was developed to assist in clinical recognition of TBX6-associated CS (TACS).

RESULTS

In cohort 1, TACS patients (N = 33) were significantly younger at onset than the remaining CS patients (P = 0.02), presented with one or more hemivertebrae/butterfly vertebrae (P = 4.9 × 10^‒8), and exhibited vertebral malformations involving the lower part of the spine (T8-S5, P = 4.4 × 10^‒3); observations were confirmed in two replication cohorts. Simple rib anomalies were prevalent in TACS patients (P = 3.1 × 10^‒7), while intraspinal anomalies were uncommon (P = 7.0 × 10^‒7). A clinically usable TACScore was developed with an area under the curve (AUC) of 0.9 (P = 1.6 × 10^‒15). A Tbx6^{-/mh (mild-hypomorphic)} mouse model supported that a gene dosage effect underlies the TACS phenotype.

CONCLUSION

TACS is a clinically distinguishable entity with consistent clinically measurable endophenotypes. The type and distribution of vertebral column abnormalities in TBX6/Tbx6 compound inheritance implicate subtle perturbations in gene dosage as a cause of spine developmental birth defects responsible for about 10% of CS.

Noncoding rare variants of TBX6 in congenital anomalies of the kidney and urinary tract

Congenital anomalies of the kidney and urinary tract (CAKUT) are a wide range of congenital structural renal defects. CAKUT is the leading cause of chronic renal failure and end-stage renal disease in children. Studies in humans and animal models have confirmed the large genetic contribution to CAKUT. The previous evidence suggested that human TBX6 coding mutations might cause CAKUT via gene-dosage insufficiency. However, the potential involvement of TBX6 noncoding mutations in CAKUT remains to be elucidated. Here, we described DNA sequencing and copy-number analysis of TBX6 in 269 Chinese subjects with CAKUT. Interestingly, we identified two heterozygous noncoding variants of TBX6 in sporadic subjects with CAKUT: one is c.769-7delT, from a subject with duplex renal and collecting system, and the other is a 3' untranslated region (3'-UTR) variant (c.1392C>T) from a subject with unilateral renal hypoplasia. These two TBX6 noncoding variants are novel and extremely rare, respectively, in human populations archived in the ExAC database. The mini-gene splicing assay showed that the TBX6 c.769-7delT variant significantly reduced the splicing efficiency of TBX6 intron 5 when compared to the wild-type control. In this work, we identified a novel splicing variant of TBX6 in human CAKUT. Our experimental observations suggested that the TBX6 noncoding variant can affect gene expression and may potentially be involved in human CAKUT.

An N-Ethyl-N-Nitrosourea (ENU) Mutagenized Mouse Model for Autosomal Dominant Nonsyndromic Kyphoscoliosis Due to Vertebral Fusion

Kyphosis and scoliosis are common spinal disorders that occur as part of complex syndromes or as nonsyndromic, idiopathic diseases. Familial and twin studies implicate genetic involvement, although the causative genes for idiopathic kyphoscoliosis remain to be identified. To facilitate these studies, we investigated progeny of mice treated with the chemical mutagen N-ethyl-N-nitrosourea (ENU) and assessed them for morphological and radiographic abnormalities. This identified a mouse with kyphoscoliosis due to fused lumbar vertebrae, which was inherited as an autosomal dominant trait; the phenotype was designated as hereditary vertebral fusion (HVF) and the locus as Hvf. Micro-computed tomography (μCT) analysis confirmed the occurrence of nonsyndromic kyphoscoliosis due to fusion of lumbar vertebrae in HVF mice, consistent with a pattern of blocked vertebrae due to failure of segmentation. μCT scans also showed the lumbar vertebral column of HVF mice to have generalized disc narrowing, displacement with compression of the neural spine, and distorted transverse processes. Histology of lumbar vertebrae revealed HVF mice to have irregularly shaped vertebral bodies and displacement of intervertebral discs and ossification centers. Genetic mapping using a panel of single nucleotide polymorphic (SNP) loci arranged in chromosome sets and DNA samples from 23 HVF (eight males and 15 females) mice, localized Hvf to chromosome 4A3 and within a 5-megabase (Mb) region containing nine protein coding genes, two processed transcripts, three microRNAs, five small nuclear RNAs, three large intergenic noncoding RNAs, and 24 pseudogenes. However, genome sequence analysis in this interval did not identify any abnormalities in the coding exons, or exon-intron boundaries of any of these genes. Thus, our studies have established a mouse model for a monogenic form of nonsyndromic kyphoscoliosis due to fusion of lumbar vertebrae, and further identification of the underlying genetic defect will help elucidate the molecular mechanisms involved in kyphoscoliosis. © 2018 The Authors. JBMR Plus is published by Wiley Periodicals, Inc. on behalf of the American Society for Bone and Mineral Research.

Screening of known disease genes in congenital scoliosis

Background

Congenital scoliosis (CS) is defined as a lateral curvature of the spine due to the vertebral malformations and has an incidence of 0.5–1/1,000 births. We previously examined TBX6 in Japanese CS patients and revealed that approximately 10% of CS was caused by TBX6 mutations. However, the genetic cause of remaining CS is unknown.

Methods

We recruited 78 CS patients without TBX6 mutations and major comorbidities, and investigated the genes previously reported to be associated with CS and congenital vertebral malformations by whole‐exome sequencing.

Results

We identified the compound heterozygous missense variants in LFNG in one patient. No likely disease‐causing variants were identified in other patients, however. LFNG encodes a GlcNAc‐transferase. The LFNG variants showed loss of their enzyme function.

Conclusions

A LFNG mutation is reported in a case of spondylocostal dysostosis (SCD), a skeletal dysplasia with severe malformations of vertebra and rib. The CS patient with LFNG mutations had multiple vertebral malformations including hemivertebrae, butterfly vertebrae, and block vertebrae, and rib malformations. LFNG mutations may cause a spectrum of phenotypes including CS and SCD. The current list of known disease genes could explain only a small fraction of genetic cause of CS.

Tbx6 Induces Nascent Mesoderm from Pluripotent Stem Cells and Temporally Controls Cardiac versus Somite Lineage Diversification

The mesoderm arises from pluripotent epiblasts and differentiates into multiple lineages; however, the underlying molecular mechanisms are unclear. Tbx6 is enriched in the paraxial mesoderm and is implicated in somite formation, but its function in other mesoderms remains elusive. Here, using direct reprogramming-based screening, single-cell RNA-seq in mouse embryos, and directed cardiac differentiation in pluripotent stem cells (PSCs), we demonstrated that Tbx6 induces nascent mesoderm from PSCs and determines cardiovascular and somite lineage specification via its temporal expression. Tbx6 knockout in mouse PSCs using CRISPR/Cas9 technology inhibited mesoderm and cardiovascular differentiation, whereas transient Tbx6 expression induced mesoderm and cardiovascular specification from mouse and human PSCs via direct upregulation of Mesp1, repression of Sox2, and activation of BMP/Nodal/Wnt signaling. Notably, prolonged Tbx6 expression suppressed cardiac differentiation and induced somite lineages, including skeletal muscle and chondrocytes. Thus, Tbx6 is critical for mesoderm induction and subsequent lineage diversification.

The coexistence of copy number variations (CNVs) and single nucleotide polymorphisms (SNPs) at a locus can result in distorted calculations of the significance in associating SNPs to disease

With the recent advance in genome-wide association studies (GWAS), disease-associated single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) have been extensively reported. Accordingly, the issue of incorrect identification of recombination events that can induce the distortion of multi-allelic or hemizygous variants has received more attention. However, the potential distorted calculation bias or significance of a detected association in a GWAS due to the coexistence of CNVs and SNPs in the same genomic region may remain under-recognized. Here we performed the association study within a congenital scoliosis (CS) cohort whose genetic etiology was recently elucidated as a compound inheritance model, including mostly one rare variant deletion CNV null allele and one common variant non-coding hypomorphic haplotype of the TBX6 gene. We demonstrated that the existence of a deletion in TBX6 led to an overestimation of the contribution of the SNPs on the hypomorphic allele. Furthermore, we generalized a model to explain the calculation bias, or distorted significance calculation for an association study, that can be 'induced' by CNVs at a locus. Meanwhile, overlapping between the disease-associated SNPs from published GWAS and common CNVs (overlap 10%) and pathogenic/likely pathogenic CNVs (overlap 99.69%) was significantly higher than the random distribution (p < 1 × 10^-6 and p = 0.034, respectively), indicating that such co-existence of CNV and SNV alleles might generally influence data interpretation and potential outcomes of a GWAS. We also verified and assessed the influence of colocalizing CNVs to the detection sensitivity of disease-associated SNP variant alleles in another adolescent idiopathic scoliosis (AIS) genome-wide association study. We proposed that detecting co-existent CNVs when evaluating the association signals between SNPs and disease traits could improve genetic model analyses and better integrate GWAS with robust Mendelian principles.

Genome-wide compound heterozygosity analysis highlighted 4 novel susceptibility loci for congenital heart disease in Chinese population

Genome-wide association studies (GWASs) have achieved great success in deciphering the genetic cause of congenital heart disease (CHD). However, the heritability of CHD remains to be clarified, and numerous genetic factors responsible for occurrence of CHD are yet unclear. In this study, we performed a genome-wide search for relaxed forms of compound heterozygosity (CH) in association with CHD using our existing GWAS data including 2265 individuals (957 CHD cases and 1308 controls). CollapsABEL was used to iteratively test the association between the CH genotype and the CHD phenotype in a sliding window manner. We highlighted 17 genetic loci showing suggestive CH-like associations with CHD (P < 5 × 10^-8 ), among which 4 genetic loci had expression quantitative trait loci (eQTL) effects in blood (P_eQTL  < 0.01). After conditional association analysis, each loci had only 1 independently effective signal reaching the significance threshold (rs2071477/rs3129299 at 6p21.32, P = 2.47 × 10^-10 ; rs10773097/rs2880921 at 12q24.31, P = 3.30 × 10^-8 ; rs73032040/rs7259476 at 19q13.11, P = 1.14 × 10^-8 ; rs10416386/rs4239517 at 19q13.31, P = 1.15 × 10^-9 ), together explained 7.83% of the CHD variance. Among these 4 associated loci, outstanding candidates for CHD-associated genes included UBC, CFM2, ZNF302, LYPD3 and CADM4. Although replication studies with larger sample size are warranted, the first CH GWAS of CHD may extend our current knowledge of the genetic contributions to CHD in the Han Chinese population.

Intrauterine phenotypic features associated with 16p11.2 recurrent microdeletions

OBJECTIVE

To investigate the detection rate of 16p11.2 recurrent microdeletions in fetuses with abnormal ultrasound findings and determine the common abnormal ultrasound findings in fetuses carrying the deletion.

METHODS

This study reviewed 2262 consecutive fetuses with abnormal ultrasound findings who underwent prenatal chromosomal microarray analysis between October 2014 and December 2016. Cases carrying the 16p11.2 recurrent microdeletion were further genetically analyzed, and their clinical features were reviewed.

RESULTS

The 16p11.2 recurrent microdeletion was identified in 12 fetuses, who had skeletal malformations (5/12), cardiovascular malformations (4/12), or isolated ultrasound markers (3/12). Approximately 0.5% (12/2262) of the fetuses with abnormal ultrasound findings harbored the deletion. The 5 fetuses with skeletal malformations displayed vertebral defects, particularly in the hemivertebra and butterfly vertebra. The detection rate of the 16p11.2 recurrent microdeletion was statistically significant (P < .05) among fetuses with skeletal malformations (3.6%, 5/140), fetuses with cardiovascular malformations (1.1%, 4/367), and fetuses with isolated ultrasound markers (0.4%, 3/702).

CONCLUSION

The most frequent ultrasound findings in fetuses with 16p11.2 recurrent microdeletions are skeletal malformations (particularly vertebral malformations), followed by cardiovascular malformations, and isolated ultrasound markers.

Response to Lefebvre et al

Congenital scoliosis (CS) is a common vertebral malformation with incidence of up to 1 of 1000 births worldwide. Recently, TBX6 has been reported as the first disease gene for CS: about 10% of CS patients are compound heterozygotes of rare null mutations and a common haplotype composed by 3 SNPs in TBX6. Lefebvre et al in this journal reported that 2 patients with spondylocostal dysostosis (SCD), a rare skeletal dysplasia affecting spine and ribs also have TBX6 mutations: 1 carried the microdeletion and a rare missense variant, and another 2 rare missense variants. We investigated the pathogenicity of the 3 missense variants in SCD by a luciferase assay. The results were negative for the proposal of Lefebvre et al. We consider these 2 SCD patients are more probably compound heterozygotes of null mutations and a common risk haplotype just as CS patients with TBX6 mutations.

Notch Signaling in Development, Tissue Homeostasis, and Disease

Notch signaling is an evolutionarily highly conserved signaling mechanism, but in contrast to signaling pathways such as Wnt, Sonic Hedgehog, and BMP/TGF-β, Notch signaling occurs via cell-cell communication, where transmembrane ligands on one cell activate transmembrane receptors on a juxtaposed cell. Originally discovered through mutations in Drosophila more than 100 yr ago, and with the first Notch gene cloned more than 30 yr ago, we are still gaining new insights into the broad effects of Notch signaling in organisms across the metazoan spectrum and its requirement for normal development of most organs in the body. In this review, we provide an overview of the Notch signaling mechanism at the molecular level and discuss how the pathway, which is architecturally quite simple, is able to engage in the control of cell fates in a broad variety of cell types. We discuss the current understanding of how Notch signaling can become derailed, either by direct mutations or by aberrant regulation, and the expanding spectrum of diseases and cancers that is a consequence of Notch dysregulation. Finally, we explore the emerging field of Notch in the control of tissue homeostasis, with examples from skin, liver, lung, intestine, and the vasculature.