Progress and perspective of TBX6 gene in congenital vertebral malformations

Zebrafish as a model to investigate a biallelic gain-of-function variant in MSGN1, associated with a novel skeletal dysplasia syndrome

BACKGROUND/OBJECTIVES

Rare genetic disorders causing specific congenital developmental abnormalities often manifest in single families. Investigation of disease-causing molecular features are most times lacking, although these investigations may open novel therapeutic options for patients. In this study, we aimed to identify the genetic cause in an Iranian patient with severe skeletal dysplasia and to model its molecular function in zebrafish embryos.

RESULTS

The proband displays short stature and multiple skeletal abnormalities, including mesomelic dysplasia of the arms with complete humero-radio-ulna synostosis, arched clavicles, pelvic dysplasia, short and thin fibulae, proportionally short vertebrae, hyperlordosis and mild kyphosis. Exome sequencing of the patient revealed a novel homozygous c.374G > T, p.(Arg125Leu) missense variant in MSGN1 (NM_001105569). MSGN1, a basic-Helix-Loop-Helix transcription factor, plays a crucial role in formation of presomitic mesoderm progenitor cells/mesodermal stem cells during early developmental processes in vertebrates. Initial in vitro experiments show protein stability and correct intracellular localization of the novel variant in the nucleus and imply retained transcription factor function. To test the pathogenicity of the detected variant, we overexpressed wild-type and mutant msgn1 mRNA in zebrafish embryos and analyzed tbxta (T/brachyury/ntl). Overexpression of wild-type or mutant msgn1 mRNA significantly reduces tbxta expression in the tailbud compared to control embryos. Mutant msgn1 mRNA injected embryos depict a more severe effect, implying a gain-of-function mechanism. In vivo analysis on embryonic development was performed by clonal msgn1 overexpression in zebrafish embryos further demonstrated altered cell compartments in the presomitic mesoderm, notochord and pectoral fin buds. Detection of ectopic tbx6 and bmp2 expression in these embryos hint to affected downstream signals due to Msgn1 gain-of-function.

CONCLUSION

In contrast to loss-of-function effects described in animal knockdown models, gain-of-function of MSGN1 explains the only mildly affected axial skeleton of the proband and rather normal vertebrae. In this context we observed notochord bending and potentially disruption of pectoral fin buds/upper extremity after overexpression of msgn1 in zebrafish embryos. The latter might result from Msgn1 function on mesenchymal stem cells or on chondrogenesis in these regions. In addition, we detected ectopic tbx6 and bmp2a expression after gain of Msgn1 function in zebrafish, which are interconnected to short stature, congenital scoliosis, limb shortening and prominent skeletal malformations in patients. Our findings highlight a rare, so far undescribed skeletal dysplasia syndrome associated with a gain-of-function mutation in MSGN1 and hint to its molecular downstream effectors.

A de novo variant in RERE causes autistic behavior by disrupting related genes and signaling pathway

Autism spectrum disorder (ASD) is a highly variable neurodevelopmental disorder that typically manifests childhood, characterized by a triad of symptoms: impaired social interaction, communication difficulties, and restricted interests with repetitive behaviors. De novo variants in related genes can cause ASD. We present the case of a 6-year-old Chinese boy with autistic behavior, including language communication impairments, intellectual disabilities, stunted development, and irritability in social interactions. Using Sanger sequencing, we confirmed a pathogenic in the RERE gene (NM_012102.4) (c.3732delC, p.Tyr1245Thrfs*12; EX21; Het). Subsequently, we generated an RERE point mutation cell line (ReMut) using CRISPR/Cas9 Targeted Genome Editing. Immunofluorescence was conducted to determine the location of the mutant RERE. RNA-sequencing and mass spectrometry analyses were performed to elucidate the ASD-related genes and signaling pathways disrupted by this variant in RERE. We identified 3790 differentially expressed genes and 684 differentially expressed proteins. The SHH signaling pathway was found to be downregulated, and the Hippo pathway was upregulated in ReMut. Genes implicated in autism, such as CNTNAP2, STX1A, FARP2, and GPC1, were significantly downregulated. Simultaneously, we noted alterations in HDAC1 and HDAC2, which are members of the WHHERE complex, suggesting their role in the pathogenesis of this patient. In conclusion, we report a de novo variant in RERE associated with autistic behavior. The finding that ASD is associated with RERE variants underscore the role of genetic factors in ASD and provides insights regarding the mechanisms underlying RERE variants in disease onset.

Pudgy mouse rib deformities emanate from abnormal paravertebral longitudinal cartilage/bone accumulations

The pudgy (pu/pu) mouse, caused by a recessive mutation in the Notch family Delta like-3 gene (Dll3), has severe rib, vertebral body and intervertebral disc abnormalities. Using whole-mount preparations and serial histologic sections we demonstrate: 1) localized paravertebral longitudinal cartilage/bone accumulations (PVLC/BAs) invariably associated with branched, fused and asymmetrically spaced ribs that emanate from it laterally; 2) abnormal rib formation immediately adjacent to abnormal vertebral body and intervertebral disc formation in asymmetric right/left fashion; and 3) patterns of rib deformation that differ in each mouse. Normal BALB/c embryo and age-matched non-affected pu/+ mice assessments allow for pu/pu comparisons. The Dll3 Notch family gene is involved in normal somitogenesis via the segmentation clock mechanism. Although pathogenesis of rib deformation is initially triggered by the Dll3 gene mutation, these findings of abnormal asymmetric costo-vertebral region structure imply that differing patterns cannot be attributed to this single gene mutation alone. All findings implicate a dual mechanism of malformation: the Dll3 gene mutation leading to subtle timing differences in traveling oscillation waves of the segmentation clock and further subsequent misdirection of tissue formation by altered chemical reaction-diffusion and epigenetic landscape responses. PVLC/BAs appear as primary supramolecular structures underlying severe rib malformation associated both with time-sensitive segmentation clock mutations and subsequent reactions.

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.

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.

Spondylocostal dysplasia and brachydactyly associated with TBX6 and IHH variants: A case report

We report a preterm male neonate presenting with a short trunk, short neck, low hairline, deformed ears, preauricular skin tag, penoscrotal transposition (PT), palmar crease, short and broad fingers and toes (brachydactyly), hypoplastic and deep‐set nails, metatarsal abductus, and cross‐fused, small echogenic kidneys. Radiologic findings and genetic studies are consistent with spondylocostal dysostosis (SCD) and autosomal dominant brachydactyly. This is the first case report of spondylocostal dysostosis and brachydactyly associated with TBX6 and IHH variants. We reviewed the literature and compared our patient's phenotype with previously reported cases of SCD.

Daw1 regulates the timely onset of cilia motility during development

Motile cilia generate cell propulsion and extracellular fluid flows that are crucial for airway clearance, fertility and left-right patterning. Motility is powered by dynein arm complexes that are assembled in the cytoplasm then imported into the cilium. Studies in Chlamydomonas reinhardtii showed that ODA16 is a cofactor which promotes dynein arm import. Here, we demonstrate that the zebrafish homolog of ODA16, Daw1, facilitates the onset of robust cilia motility during development. Without Daw1, cilia showed markedly reduced motility during early development; however, motility subsequently increased to attain close to wild-type levels. Delayed motility onset led to differential effects on early and late cilia-dependent processes. Remarkably, abnormal body axis curves, which formed during the first day of development due to reduced cilia motility, self-corrected when motility later reached wild-type levels. Zebrafish larva therefore possess the ability to survey and correct body shape abnormalities. This work defines Daw1 as a factor which promotes the onset of timely cilia motility and can explain why human patients harboring DAW1 mutations exhibit significant laterality perturbations but mild airway and fertility complications.

Summary: Daw1 promotes the onset of timely cilia motility for robust axial straightening during zebrafish development.

Craniocervical abnormalities in osteogenesis imperfecta type V

Craniocervical abnormalities in osteogenesis imperfecta (OI) such as basilar invagination or cervical kyphosis can cause severe neurological morbidity. These abnormalities may be more frequent in OI type V compared with other OI subtypes of similar disease severity, underlining the importance of screening in this group.

INTRODUCTION

Craniocervical abnormalities in osteogenesis imperfecta (OI) can cause severe neurological morbidity. Although radiological cranial base abnormalities in OI have been well described in the literature, there are limited data on these abnormalities in OI type V and their association with clinical sequelae.

METHODS

A retrospective case series on patients with craniocervical abnormalities in OI type V at our institution.

RESULTS

Craniocervical abnormalities were present in 7 of 37 patients with OI type V (19%). For 5 patients (age at last follow-up: 5 to 26 years; 2 females), sufficient information was available for inclusion in the case series. All had genetically confirmed OI type V. Age range at diagnosis of the craniocervical abnormality was 1 day to 18 years. Basilar invagination was present in 3 patients; 2 had cervical kyphosis. Dysplasia of upper cervical vertebrae or base of skull was seen in 3 patients. The severity of the craniocervical abnormality did not clearly correlate with the severity of the OI phenotype. Three patients required surgical intervention (ages 7, 11, and 26 years) due to compression of the spinal cord or brainstem. Craniocervical abnormalities were detected incidentally or on screening in 3 patients, and only 2 had significant positive findings on neurological examination.

CONCLUSION

A variety of craniocervical abnormalities are seen in OI type V including dysplasia of the cervical vertebrae. These cases highlight the importance of screening patients with OI type V with lateral skull and cervical spine x-rays throughout childhood and after skeletal maturity.

Integration of genetic, transcriptomic, and clinical data provides insight into 16p11.2 and 22q11.2 CNV genes

BACKGROUND

Deletions and duplications of the multigenic 16p11.2 and 22q11.2 copy number variant (CNV) regions are associated with brain-related disorders including schizophrenia, intellectual disability, obesity, bipolar disorder, and autism spectrum disorder (ASD). The contribution of individual CNV genes to each of these identified phenotypes is unknown, as well as the contribution of these CNV genes to other potentially subtler health implications for carriers. Hypothesizing that DNA copy number exerts most effects via impacts on RNA expression, we attempted a novel in silico fine-mapping approach in non-CNV carriers using both GWAS and biobank data.

METHODS

We first asked whether gene expression level in any individual gene in the CNV region alters risk for a known CNV-associated behavioral phenotype(s). Using transcriptomic imputation, we performed association testing for CNV genes within large genotyped cohorts for schizophrenia, IQ, BMI, bipolar disorder, and ASD. Second, we used a biobank containing electronic health data to compare the medical phenome of CNV carriers to controls within 700,000 individuals in order to investigate the full spectrum of health effects of the CNVs. Third, we used genotypes for over 48,000 individuals within the biobank to perform phenome-wide association studies between imputed expressions of individual 16p11.2 and 22q11.2 genes and over 1500 health traits.

RESULTS

Using large genotyped cohorts, we found individual genes within 16p11.2 associated with schizophrenia (TMEM219, INO80E, YPEL3), BMI (TMEM219, SPN, TAOK2, INO80E), and IQ (SPN), using conditional analysis to identify upregulation of INO80E as the driver of schizophrenia, and downregulation of SPN and INO80E as increasing BMI. We identified both novel and previously observed over-represented traits within the electronic health records of 16p11.2 and 22q11.2 CNV carriers. In the phenome-wide association study, we found seventeen significant gene-trait pairs, including psychosis (NPIPB11, SLX1B) and mood disorders (SCARF2), and overall enrichment of mental traits.

CONCLUSIONS

Our results demonstrate how integration of genetic and clinical data aids in understanding CNV gene function and implicates pleiotropy and multigenicity in CNV biology.

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.

Identification of novel FBN1 variations implicated in congenital scoliosis

Congenital scoliosis (CS) is a form of scoliosis caused by congenital vertebral malformations. Genetic predisposition has been demonstrated in CS. We previously reported that TBX6 loss-of-function causes CS in a compound heterozygous model; however, this model can explain only 10% of CS. Many monogenic and polygenic CS genes remain to be elucidated. In this study, we analyzed exome sequencing (ES) data of 615 Chinese CS from the Deciphering Disorders Involving Scoliosis and COmorbidities (DISCO) project. Cosegregation studies for 103 familial CS identified a novel heterozygous nonsense variant, c.2649G>A (p.Trp883Ter) in FBN1. The association between FBN1 and CS was then analyzed by extracting FBN1 variants from ES data of 574 sporadic CS and 828 controls; 30 novel variants were identified and prioritized for further analyses. A mutational burden test showed that the deleterious FBN1 variants were significantly enriched in CS subjects (OR = 3.9, P = 0.03 by Fisher’s exact test). One missense variant, c.2613A>C (p.Leu871Phe) was recurrent in two unrelated CS subjects, and in vitro functional experiments for the variant suggest that FBN1 may contribute to CS by upregulating the transforming growth factor beta (TGF-β) signaling. Our study expanded the phenotypic spectrum of FBN1, and provided nove insights into the genetic etiology of CS.

Comparative analysis of serum proteome in congenital scoliosis patients with TBX6 haploinsufficiency - a first report pointing to lipid metabolism

Congenital scoliosis (CS) is a three‐dimensional deformity of the spine affecting quality of life. We have demonstrated TBX6 haploinsufficiency is the most important contributor to CS. However, the pathophysiology at the protein level remains unclear. Therefore, this study was to explore the differential proteome in serum of CS patients with TBX6 haploinsufficiency. Sera from nine CS patients with TBX6 haploinsufficiency and nine age‐ and gender‐matched healthy controls were collected and analysed by isobaric tagged relative and absolute quantification (iTRAQ) labelling coupled with mass spectrometry (MS). In total, 277 proteins were detected and 20 proteins were designated as differentially expressed proteins, which were submitted to subsequent bioinformatics analysis. Gene Ontology classification analysis showed the biological process was primarily related to ‘cellular process’, molecular function ‘structural molecule activity’ and cellular component ‘extracellular region’. IPA analysis revealed ‘LXR/RXR activation’ was the top pathway, which is a crucial pathway in lipid metabolism. Hierarchical clustering analysis generated two clusters. In summary, this study is the first proteomic research to delineate the total and differential serum proteins in TBX6 haploinsufficiency‐caused CS. The proteins discovered in this experiment may serve as potential biomarkers for CS, and lipid metabolism might play important roles in the pathogenesis of CS.