Identification of novel FBN1 variations implicated in congenital scoliosis

Mitochondrial genes modulate the phenotypic expression of congenital scoliosis syndrome caused by mutations in the TBXT gene

BACKGROUND

Congenital scoliosis (CS) is a spinal disorder caused by genetic-congenital vertebral malformations and may be associated with other congenital defects or may occur alone. It is genetically heterogeneous and numerous genes contributing to this disease have been identified. In addition, CS has a wide range of phenotypic and genotypic variability, which has been explained by the intervention of genetic factors like modifiers and environment genes. The aim of the present study was to determine the possible cause of CS in a Tunisian patient and to examine the association between mtDNA mutations and mtDNA content and CS.

METHODS

Here we performed Whole-Exome Sequencing (WES) in a patient presenting clinical features suggestive of severe congenital scoliosis syndrome. Direct sequencing of the whole mitochondrial DNA (mtDNA) was also performed in addition to copy number quantification in the blood of the indexed case. In silico prediction tools, 3D modeling and molecular docking approaches were used.

RESULTS

The WES revealed the homozygous missense mutation c.512A > G (p.H171R) in the TBXT gene. Bioinformatic analysis demonstrated that the p.H171R variant was highly deleterious and caused the TBXT structure instability. Molecular docking revealed that the p.H171R mutation disrupted the monomer stability which seemed to be crucial for maintaining the stability of the homodimer and consequently to the destabilization of the homodimer-DNA complex. On the other hand, we hypothesized that mtDNA can be a modifier factor, so, the screening of the whole mtDNA showed a novel heteroplasmic m.10150T > A (p.M31K) variation in the MT-ND3 gene. Further, qPCR analyses of the patient's blood excluded mtDNA depletion. Bioinformatic investigation revealed that the p.M31K mutation in the ND3 protein was highly deleterious and may cause the ND3 protein structure destabilization and could disturb the interaction between complex I subunits.

CONCLUSION

We described the possible role of mtDNA genetics on the pathogenesis of congenital scoliosis by hypothesizing that the presence of the homozygous variant in TBXT accounts for the CS phenotype in our patient and the MT-ND3 gene may act as a modifier gene.

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.

COL11A2 as a candidate gene for vertebral malformations and congenital scoliosis

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

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.

[Latest advances in the diagnosis and treatment of Marfan syndrome]

Marfan syndrome (MFS) is a multisystem connective tissue disease with autosomal dominant inheritance. It is mainly caused by FBN1 gene mutation and often has different clinical manifestations. Neonatal MFS is especially rare with severe conditions and a poor prognosis. At present, there is still no radical treatment method for MFS, but early identification, early diagnosis, and early treatment can effectively prolong the life span of patients. This article reviews the latest advances in the diagnosis and treatment of MFS.

Understanding Musculoskeletal Disorders Through Next-Generation Sequencing

»

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

»

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

»

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

»

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

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.

CRISPR/Cas9 in zebrafish: An attractive model for FBN1 genetic defects in humans

Background

Mutations in the fibrillin‐1 gene (FBN1) are associated with various heritable connective tissue disorders (HCTD). The most studied HCTD is Marfan syndrome. Ninety percent of Marfan syndrome is caused by mutations in the FBN1 gene. The zebrafish share high genetic similarity to humans, representing an ideal model for genetic research of human diseases. This study aimed to generate and characterize fbn1^+/− mutant zebrafish using the CRISPR/Cas9 gene‐editing technology.

Methods

CRISPR/Cas9 was applied to generate an fbn1 frameshift mutation (fbn1^+/−) in zebrafish. F1 fbn1^+/− heterozygotes were crossed with transgenic fluorescent zebrafish to obtain F2 fbn1^+/− zebrafish. Morphological abnormalities were assessed in F2 fbn1^+/− zebrafish by comparing with the Tuebingen (TU) wild‐type controls at different development stages.

Results

We successfully generated a transgenic line of fbn1^+/− zebrafish. Compared with TU wild‐type zebrafish, F2 fbn1^+/− zebrafish exhibited noticeably decreased pigmentation, increased lengths, slender body shape, and abnormal cardiac blood flow from atrium to ventricle.

Conclusion

We generated the first fbn1^+/− zebrafish model using CRISPR/Cas9 gene‐editing approach to mimic FBN1 genetic defects in humans, providing an attractive model of Marfan syndrome and a method to determine the pathogenicity of gene mutation sites.

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.

Proof-of-Concept: Antisense Oligonucleotide Mediated Skipping of Fibrillin-1 Exon 52

Marfan syndrome is one of the most common dominantly inherited connective tissue disorders, affecting 2–3 in 10,000 individuals, and is caused by one of over 2800 unique FBN1 mutations. Mutations in FBN1 result in reduced fibrillin-1 expression, or the production of two different fibrillin-1 monomers unable to interact to form functional microfibrils. Here, we describe in vitro evaluation of antisense oligonucleotides designed to mediate exclusion of FBN1 exon 52 during pre-mRNA splicing to restore monomer homology. Antisense oligonucleotide sequences were screened in healthy control fibroblasts. The most effective sequence was synthesised as a phosphorodiamidate morpholino oligomer, a chemistry shown to be safe and effective clinically. We show that exon 52 can be excluded in up to 100% of FBN1 transcripts in healthy control fibroblasts transfected with PMO52. Immunofluorescent staining revealed the loss of fibrillin 1 fibres with ~50% skipping and the subsequent re-appearance of fibres with >80% skipping. However, the effect of exon skipping on the function of the induced fibrillin-1 isoform remains to be explored. Therefore, these findings demonstrate proof-of-concept that exclusion of an exon from FBN1 pre-mRNA can result in internally truncated but identical monomers capable of forming fibres and lay a foundation for further investigation to determine the effect of exon skipping on fibrillin-1 function.

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

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

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

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

FBN1 promotes DLBCL cell migration by activating the Wnt/β-catenin signaling pathway and regulating TIMP1

The heterogeneity of diffuse large B-cell lymphoma (DLBCL) acts as a main barrier to identify the genetic basis of the disease and the choice of treatment. Differentially expressed genes (DEGs) from three mRNA expression profile datasets were screened using GEO2R, and bioinformatics analysis was performed on the DEGs. A total of six upregulated and 13 downregulated DEGs were identified. Among these, two hub genes with a high degree of correlation were selected. FBN1 and TIMP1 were identified via STRING analysis and validated by GEPIA. FBN1 and TIMP1 were highly expressed in DLBCL tissues. FBN1 expression was significantly higher in patients of the Ann Arbor stage group (III-IV), with higher IPI score (3-5), and in the non-GCB group. Patients with high TIMP1 expression were more frequently associated with B symptoms, Ann Arbor stage (III-IV), higher IPI score (3-5) and were in the non-GCB group. Furthermore, FBN1 siRNA decreased FBN1 and TIMP1 expression and downregulation of TIMP1 attenuated TIMP1 expression but not of FBN1. Migration of DLBCL cells reduced when treated with either FBN1 or TIMP1 siRNA. Moreover, FBN1 or TIMP1 siRNA decreased the expression of Wnt target genes. Simultaneous overexpression of TIMP1 resulted in an increase in these proteins. This confirmed that both FBN1 and TIMP1 were positively associated with DLBCL progression. Further analysis revealed that FBN1/TIMP1 interaction could improve DLBCL cell migration and regulate the Wnt signaling pathway. Although the underlying mechanisms regarding the interaction between FBN1 and TIMP1 requires further clarification, they might be potential therapeutic targets for DLBCL therapy.

Body Posture Defects and Body Composition in School-Age Children

The aim of the study was to assess the relationship between the shape of the anteriorposterior spinal curvature and body composition in schoolchildren. The study included 257 children, aged 11–12. Correct spinal curvature was established in 106 (41.08%) subjects. Other types included: decreased kyphosis and correct lordosis—40 participants (15.50%), correct kyphosis and decreased lordosis—24 individuals (9.30%), increased kyphosis and correct lordosis—17 subjects (6.59%), correct kyphosis and increased lordosis—22 children (8.53%), decreased kyphosis and decreased lordosis—32 people (12.40%), decreased kyphois and increased lordosis—four of the examined subjects (1.55%) increased kyphosis and lordosis—13 people (5.04%). In addition, 134 (51.94%) demonstrated scoliotic posture and eight (3.10%) scoliosis. There were significant relationships between the shape of the anteriorposterior curvatures and body composition in schoolchildren. Those with a strong body build (predominance of mesomorphs) were generally characterised by the correct formation of these curvatures. In contrast, lean subjects (with the predominance of ectomorphic factors) were more likely to experience abnormalities. No correlations with body composition were observed in the group with scoliotic posture or scoliosis. Both in the prevention and correction of postural defects, one should gradually move away from one-sided, usually one-system, therapeutic effects. An approach that takes into account both somatic and neurophysiological factors seems appropriate. With the correct body composition and structure, shaping the habit of correct posture is much easier.

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.

Genome-Wide Analysis of circular RNAs and validation of hsa_circ_0006719 as a potential novel diagnostic biomarker in congenital scoliosis patients

Congenital scoliosis (CS) is a form of spinal curvature resulting from anomalous development of vertebrae. Recent studies demonstrated that circRNAs could serve as potential biomarkers of disease diagnosis. Genome‐wide circRNAs expression in seven CS patients and three healthy controls was initially detected. Bioinformatics analysis was conducted to explore the potential pathological pathway of CS. Quantitative PCR (qPCR) was performed to validate the selected circRNAs in the replication cohort with 32 CS patients and 30 healthy controls. Logistic regression controlling for gender was conducted to compare the expression difference. Receiver operating characteristic (ROC) curve analysis was performed to evaluate the diagnostic value. Twenty‐two differentially expressed circRNAs were filtered from genome‐wide circRNA sequencing. Seven circRNAs were validated by qPCR. Only hsa_circ_0006719 was confirmed to have a higher expression level in the CS group than the healthy control group (P = 0.036). Receiver operating characteristic curve also suggested that hsa_circ_0006719 had significant diagnostic value for CS (AUC = 0.739, P = 0.001). We described the first study of circRNAs in CS and validated hsa_circ_0006719 as a potential novel diagnostic biomarker of CS.