TBX6 null variants and a common hypomorphic allele in congenital scoliosis

Prenatal Diagnosis of Chromosome 16p11.2 Microdeletion

(1) Objective: To investigate the prenatal diagnosis and genetic counseling for 16p11.2 microdeletion syndrome and to evaluate its pregnancy outcome. (2) Methods: This study included 4968 pregnant women who selected invasive prenatal diagnoses from 1 January 2017 to 1 August 2022. These 4698 pregnancies underwent chromosomal microarray analysis (CMA), data on 81 fetuses diagnosed with 16p11.2 microdeletion syndrome based on prenatal ultrasound features and genetic test results were recorded, and their pregnancy outcome was evaluated. (3) Results: 1.63% of fetuses (81/4968) were diagnosed with 16p11.2 microdeletion syndrome. Among these, there were skeletal malformations in 48.15% of the 81 fetuses, cardiovascular malformations in 30.86%, central nervous system malformations (CNS) in 11.11%, digestive system structural abnormalities in 6.17%, and isolated ultrasonography markers in 3.70%. (4) Conclusions: 16p11.2 microdeletion syndrome can display various systemic ultrasound abnormalities in the perinatal period but vertebral malformations are the most common. Our study is the first to report that TBX1 and CJA5 are associated with 16p11.2 microdeletion syndrome, expanding the disease spectrum of 16p11.2 microdeletion syndrome. In our study, the ventricular septal defect is the main feature of cardiac structural abnormalities caused by 16p11.2 microdeletion syndrome. In addition, our study highlights the use of CMA in 16p11.2 microdeletion syndrome, analyzed their genetic results, and evaluated the follow-up prognosis, which can be useful for prenatal diagnosis and genetic counseling.

Utility of trio-based prenatal exome sequencing incorporating splice-site and mitochondrial genome assessment in pregnancies with fetal ultrasound anomalies: prospective cohort study

OBJECTIVE

To evaluate the utility of trio-based prenatal exome sequencing (pES), incorporating splice-site and mitochondrial genome assessment, in the prenatal diagnosis of fetuses with ultrasound anomalies and normal copy-number variant sequencing (CNV-seq) results.

METHODS

This was a prospective study of 90 ongoing pregnancies with ultrasound anomalies that underwent trio-based pES after receiving normal CNV-seq results, from September 2020 to November 2021, in a single center in China. By using pES with a panel encompassing exome coding and splicing regions as well as mitochondrial genome for fetuses and parents, we identified the underlying genetic causes of fetal anomalies, incidental fetal findings and parental carrier status. Information on pregnancy outcome and the impact of pES findings on parental decision-making was collected.

RESULTS

Of the 90 pregnancies included, 28 (31.1%) received a diagnostic result that could explain the fetal ultrasound anomalies. The highest diagnostic yield was noted for brain abnormalities (3/6 (50.0%)), followed by hydrops (4/9 (44.4%)) and skeletal abnormalities (13/34 (38.2%)). Collectively, 34 variants of 20 genes were detected in the 28 diagnosed cases, with 55.9% (19/34) occurring de novo. Variants of uncertain significance (VUS) associated with fetal phenotypes were detected in six (6.7%) fetuses. Interestingly, fetal (n = 4) and parental (n = 3) incidental findings (IFs) were detected in seven (7.8%) cases. These included two fetuses carrying a de-novo likely pathogenic (LP) variant of the CIC and FBXO11 genes, respectively, associated with neurodevelopmental disorders, and one fetus with a LP variant in a mitochondrial gene. The remaining fetus presented with unilateral renal dysplasia and was incidentally found to carry a pathogenic PKD1 gene variant resulting in adult-onset polycystic kidney, which was later confirmed to be inherited from the mother. In addition, parental heterozygous variants associated with autosomal recessive diseases were detected in three families, including one with additional fetal diagnostic findings. Diagnostic results or fetal IFs contributed to parental decision-making about termination of the pregnancy in 26 families (26/72 (36.1%)), while negative pES results or identification of VUS encouraged 40 families (40/72 (55.6%)) to continue their pregnancy, which ended in a live birth in all cases.

CONCLUSION

Trio-based pES can provide additional genetic information for pregnancies with fetal ultrasound anomalies without a CNV-seq diagnosis. The incidental findings and parental carrier status reported by trio-based pES with splice-site and mitochondrial genome analysis extend its clinical application, but careful genetic counseling is warranted. © 2022 International Society of Ultrasound in Obstetrics and Gynecology.

Early-Onset and Severe Complex Hereditary Spastic Paraplegia Caused by De Novo Variants in SPAST

BACKGROUND

Familial hereditary spastic paraplegia (HSP)-SPAST (SPG4) typically presents with a pure HSP phenotype.

OBJECTIVE

The aim of this study was to delineate the genotypic and phenotypic spectrum of children with de novo HSP-SPAST.

METHODS

This study used a systematic cross-sectional analysis of clinical and molecular features.

RESULTS

We report the clinical and molecular spectrum of 40 patients with heterozygous pathogenic de novo variants in SPAST (age range: 2.2-27.7 years). We identified 19 unique variants (16/40 carried the same recurrent variant, p.Arg499His). Symptom onset was in early childhood (median: 11.0 months, interquartile range: 6.0 months) with significant motor and speech delay, followed by progressive ascending spasticity, dystonia, neurogenic bladder dysfunction, gastrointestinal dysmotility, and epilepsy. The mean Spastic Paraplegia Rating Scale score was 32.8 ± 9.7 (standard deviation).

CONCLUSIONS

These results confirm that de novo variants in SPAST lead to a severe and complex form of HSP that differs from classic familial pure HSP-SPAST. Clinicians should be aware of this syndrome in the differential diagnosis for cerebral palsy. © 2022 International Parkinson and Movement Disorder Society.

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.

Clinical genetics of spondylocostal dysostosis: A mini review

Spondylocostal dysostosis is a genetic defect associated with severe rib and vertebrae malformations. In recent years, extensive clinical and molecular diagnosis advancements enabled us to identify disease-causing variants in different genes for such severe conditions. The identification of novel candidate genes enabled us to understand the developmental biology and molecular and cellular mechanisms involved in the etiology of these rare diseases. Here, we discuss the clinical and molecular targets associated with spondylocostal dysostosis, including clinical evaluation, genes, and pathways involved. This review might help us understand the basics of such a severe disorder, which might help in proper clinical characterization and help in future therapeutic strategies.

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

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

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

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

Sequencing individual genomes with recurrent genomic disorder deletions: an approach to characterize genes for autosomal recessive rare disease traits

BACKGROUND

In medical genetics, discovery and characterization of disease trait contributory genes and alleles depends on genetic reasoning, study design, and patient ascertainment; we suggest a segmental haploid genetics approach to enhance gene discovery and molecular diagnostics.

METHODS

We constructed a genome-wide map for nonallelic homologous recombination (NAHR)-mediated recurrent genomic deletions and used this map to estimate population frequencies of NAHR deletions based on large-scale population cohorts and region-specific studies. We calculated recessive disease carrier burden using high-quality pathogenic or likely pathogenic variants from ClinVar and gnomAD. We developed a NIRD (NAHR deletion Impact to Recessive Disease) score for recessive disorders by quantifying the contribution of NAHR deletion to the overall allele load that enumerated all pairwise combinations of disease-causing alleles; we used a Punnett square approach based on an assumption of random mating. Literature mining was conducted to identify all reported patients with defects in a gene with a high NIRD score; meta-analysis was performed on these patients to estimate the representation of NAHR deletions in recessive traits from contemporary human genomics studies. Retrospective analyses of extant clinical exome sequencing (cES) were performed for novel rare recessive disease trait gene and allele discovery from individuals with NAHR deletions.

RESULTS

We present novel genomic insights regarding the genome-wide impact of NAHR recurrent segmental variants on recessive disease burden; we demonstrate the utility of NAHR recurrent deletions to enhance discovery in the challenging context of autosomal recessive (AR) traits and biallelic variation. Computational results demonstrate new mutations mediated by NAHR, involving recurrent deletions at 30 genomic regions, likely drive recessive disease burden for over 74% of loci within these segmental deletions or at least 2% of loci genome-wide. Meta-analyses on 170 literature-reported patients implicate that NAHR deletions are depleted from the ascertained pool of AR trait alleles. Exome reanalysis of personal genomes from subjects harboring recurrent deletions uncovered new disease-contributing variants in genes including COX10, ERCC6, PRRT2, and OTUD7A.

CONCLUSIONS

Our results demonstrate that genomic sequencing of personal genomes with NAHR deletions could dramatically improve allele and gene discovery and enhance clinical molecular diagnosis. Moreover, results suggest NAHR events could potentially enable human haploid genetic screens as an approach to experimental inquiry into disease biology.

De novo variants in EMC1 lead to neurodevelopmental delay and cerebellar degeneration and affect glial function in Drosophila

BACKGROUND

The endoplasmic reticulum (ER)-membrane protein complex (EMC) is a multi-protein transmembrane complex composed of 10 subunits that functions as a membrane-protein chaperone. Variants in EMC1 lead to neurodevelopmental delay and cerebellar degeneration. Multiple families with biallelic variants have been published, yet to date, only a single report of a monoallelic variant has been described, and functional evidence is sparse.

METHODS

Exome sequencing was used to investigate the genetic cause underlying severe developmental delay in three unrelated children. EMC1 variants were modeled in Drosophila, using loss-of-function (LoF) and overexpression studies. Glial-specific and neuronal-specific assays were used to determine whether the dysfunction was specific to one cell type.

RESULTS

Exome sequencing identified de novo variants in EMC1 in three individuals affected by global developmental delay, hypotonia, seizures, visual impairment and cerebellar atrophy. All variants were located at Pro582 or Pro584. Drosophila studies indicated that imbalance of EMC1-either overexpression or knockdown-results in pupal lethality and suggest that the tested homologous variants are LoF alleles. In addition, glia-specific gene dosage, overexpression or knockdown, of EMC1 led to lethality, whereas neuron-specific alterations were tolerated.

DISCUSSION

We establish de novo monoallelic EMC1 variants as causative of a neurological disease trait by providing functional evidence in a Drosophila model. The identified variants failed to rescue the lethality of a null allele. Variations in dosage of the wild-type EMC1, specifically in glia, lead to pupal lethality, which we hypothesize results from the altered stoichiometry of the multi-subunit protein complex EMC.

Prenatal Diagnosis and Outcomes in Fetuses with Hemivertebra

Background: There are few studies on the burden of chromosomal abnormalities and single gene disorders in fetal hemivertebra (HV). We aim to investigate the cytogenetic and monogenic risk and evaluate prenatal outcomes of fetal HV. Method: This study included fetuses diagnosed with HV divided into two groups: isolated HV and non-isolated HV. Data on other sonographic structural anomalies, chromosomal and sub-chromosomal abnormalities, monogenic variations detected by WES, and prenatal outcomes are recorded and reviewed. Results: Among 109 fetal HV cases, forty-seven (43.1%) non-isolated HV cases were associated with structural anomalies. Chromosomal test results were available in 58 cases, identifying six (10.3%) chromosomal aberrations involved in four isolated and two non-isolated HV. WES identified four (likely) pathogenic variants in three cases among 16 fetuses with HV, involving three novel variants, 1250G > T and c.1277G> inherited from parents, respectively, in DLL3 and c.7213C > A ** in the FLNB. The live birth rate (LB) was higher in the isolated fetal HV group than in the non-isolated group (67.7% (42/62) vs. 12.5% (12/47), p < 0.001). Conclusion: This study emphasizes the risk of cytogenetic abnormalities in isolated HV. WES yields a diagnostic rate of 18.3% in HV with normal CMA, probably aiding the prenatal counseling and management of fetal HV.

Network meta-analysis of the optimal time of applying enteral immunonutrition in esophageal cancer patients receiving esophagectomy

BACKGROUND

Enteral immunonutrition (EIN) has been extensively applied in cancer patients; however, its role in esophageal cancer (EC) patients receiving esophagectomy remains unclear. We performed this network meta-analysis to investigate the impact of EIN on patients undergoing surgery for EC and further determine the optimal time of applying EIN.

METHODS

We searched PubMed, EMBASE, Cochrane library, and China National Knowledgement Infrastructure (CNKI) to identify eligible studies. Categorical data was expressed as the odds ratio with 95% confidence interval (CI), and continuous data was expressed as mean difference (MD) with 95% CI. Meta-analysis with head-to-head approach and network meta-analysis was performed to evaluate the impact of EIN on clinical outcomes using RevMan 5.3 and ADDIS V.1.16.8 software. The surface under the cumulative ranking curve (SUCRA) was calculated to rank all nutritional regimes.

RESULTS

Total 14 studies involving 1071 patients were included. Meta-analysis with head-to-head approach indicated no difference between EIN regardless of the application time and standard EN (SEN); however, subgroup analyses found that postoperative EIN was associated with decreased incidence of total infectious complications (OR = 0.47; 95%CI = 0.26 to 0.84; p = 0.01) and pneumonia (OR = 0.47; 95%CI = 0.25 to 0.90; p = 0.02) and shortened the length of hospitalization (LOH) (MD =  - 1.01; 95%CI =  - 1.44 to - 0.57; p < 0.001) compared to SEN, which were all supported by network meta-analyses. Ranking probability analysis further indicated that postoperative EIN has the highest probability of being the optimal option in terms of these three outcomes.

CONCLUSIONS

Postoperative EIN should be preferentially utilized in EC patients undergoing esophagectomy because it has optimal potential of decreasing the risk of total infectious complications and pneumonia and shortening LOH.

OSF REGISTRATION NUMBER

10.17605/OSF.IO/KJ9UY.

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.

A novel CLCNKB variant in a Chinese family with classic Bartter syndrome and prenatal genetic diagnosis

Background

Type III Bartter syndrome (BS), often known as classic Bartter syndrome is caused by variants in CLCNKB gene, which encoding the basolateral chloride channel protein ClC‐Kb, and is characterized by renal salt wasting, hypokalemia, metabolic alkalosis, increased renin, and aldosterone levels.

Methods

A 2‐year‐old boy presented severe malnutrition, severe metabolic alkalosis and severe hypokalemia and was clinically diagnosed with BS. The trio exome sequencing (ES) was performed to discover the genetic cause of this patient, followed by validation using Sanger sequencing and quantitative polymerase chain reaction subsequently.

Results

The genetic analysis indicated that this patient with a compound heterozygous variants of CLCNKB gene including a novel nonsense variant c.876 T > A and a whole‐gene deletion. The two variants were inherited from his parents, respectively. Subsequently, target sequencing of CLCNKB gene was performed for next pregnancy, and prenatal genetic diagnosis was provided for the family.

Conclusions

The results of current study identified the compound heterozygous variants in a patient with classic BS. The novel variant expands the spectrum of CLCNKB variants in BS. Our study also indicates that ES is an alternative tool to simultaneously detect single‐nucleotide variants and copy‐number variants.

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.

Possible association of 16p11.2 copy number variation with altered lymphocyte and neutrophil counts

Recurrent copy-number variations (CNVs) at chromosome 16p11.2 are associated with neurodevelopmental diseases, skeletal system abnormalities, anemia, and genitourinary defects. Among the 40 protein-coding genes encompassed within the rearrangement, some have roles in leukocyte biology and immunodeficiency, like SPN and CORO1A. We therefore investigated leukocyte differential counts and disease in 16p11.2 CNV carriers. In our clinically-recruited cohort, we identified three deletion carriers from two families (out of 32 families assessed) with neutropenia and lymphopenia. They had no deleterious single-nucleotide or indel variant in known cytopenia genes, suggesting a possible causative role of the deletion. Noticeably, all three individuals had the lowest copy number of the human-specific BOLA2 duplicon (copy-number range: 3–8). Consistent with the lymphopenia and in contrast with the neutropenia associations, adult deletion carriers from UK biobank (n = 74) showed lower lymphocyte (Padj = 0.04) and increased neutrophil (Padj = 8.31e-05) counts. Mendelian randomization studies pinpointed to reduced CORO1A, KIF22, and BOLA2-SMG1P6 expressions being causative for the lower lymphocyte counts. In conclusion, our data suggest that 16p11.2 deletion, and possibly also the lowest dosage of the BOLA2 duplicon, are associated with low lymphocyte counts. There is a trend between 16p11.2 deletion with lower copy-number of the BOLA2 duplicon and higher susceptibility to moderate neutropenia. Higher numbers of cases are warranted to confirm the association with neutropenia and to resolve the involvement of the deletion coupled with deleterious variants in other genes and/or with the structure and copy number of segments in the CNV breakpoint regions.

Towards clinical applications of in vitro-derived axial progenitors

The production of the tissues that make up the mammalian embryonic trunk takes place in a head-tail direction, via the differentiation of posteriorly-located axial progenitor populations. These include bipotent neuromesodermal progenitors (NMPs), which generate both spinal cord neurectoderm and presomitic mesoderm, the precursor of the musculoskeleton. Over the past few years, a number of studies have described the derivation of NMP-like cells from mouse and human pluripotent stem cells (PSCs). In turn, these have greatly facilitated the establishment of PSC differentiation protocols aiming to give rise efficiently to posterior mesodermal and neural cell types, which have been particularly challenging to produce using previous approaches. Moreover, the advent of 3-dimensional-based culture systems incorporating distinct axial progenitor-derived cell lineages has opened new avenues toward the functional dissection of early patterning events and cell vs non-cell autonomous effects. Here, we provide a brief overview of the applications of these cell types in disease modelling and cell therapy and speculate on their potential uses in the future.

Biology in balance: human diploid genome integrity, gene dosage, and genomic medicine

The path to completion of the functional annotation of the haploid human genome reference build, exploration of the clan genomics hypothesis, understanding human gene and genome functional biology, and gene genome and organismal evolution, is in reach.

Clinical Implications of Chromosome 16 Copy Number Variation

Chromosome 16 is one of the gene-rich chromosomes; however, approximately 10% of the chromosome 16 sequence is composed of segmental copies, which renders this chromosome instable and predisposes it to rearrangements via frequent nonallelic homologous recombination. Microarray technologies have enabled the analysis of copy number variations (CNV), which may be associated with the risk of developing complex diseases. Through comparative genomic hybridisation in 1,298 patients, we detected 18 cases with chromosome 16 CNV. We identified 2recurrent CNV regions, including 1 at 16p13.11 in 4 patients and another at 16p11.2 in 7 patients. We also detected atypical chromosome 16 rearrangements in 7 patients. Furthermore, we noted an increased frequency of co-occurring genomic changes, supporting the two-hit hypothesis to explain the phenotypic variability in the clinical presentation of CNV syndromes. Our findings can contribute to the creation of a chromosome 16 disease map based on regions that may be associated with disease development.

Delineation of dual molecular diagnosis in patients with skeletal deformity

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Developments in Congenital Scoliosis and Related Research from 1992 to 2021: A Thirty-Year Bibliometric Analysis

BACKGROUND

Many studies, mainly original articles and reviews, have been reported on congenital scoliosis (CS), but there is a lack of bibliometric analyses. This study aimed to systematically analyze the developments and focuses in CS and related research fields.

METHODS

Data were retrieved from the Web of Science Core Collection database, and the top 100 most-cited studies were analyzed emphatically. The Web of Science Results Analysis and Citation Report was used to analyze different aspects of the literature. CiteSpace was used to analyze the cooperation network, reference co-citation, burst keywords, and burst citations.

RESULTS

The final analysis included 749 studies. CS and related research has been rapidly expanding. Several journals have published relevant studies and most-cited studies on this topic. Cooperation was noted among authors, institutions, and countries/regions in multiple instances. Surgical techniques (hemivertebra resection/posterior vertebral column resection/vertical expandable prosthetic titanium rib/double approach/grade 4 osteotomy) was one of the most common research focuses. In addition, research on genetics and molecular biology related to CS has become an emerging trend as a result of advances in basic science.

CONCLUSIONS

Over time, research on CS and in related fields has gained greater attention and has been expanding continuously, showing a trend toward globalization. We recommend that researchers focus on the progress of surgical techniques, advances in molecular biology and genetics, and characteristics of CS. The top clusters, most-cited articles, and references with the strongest burst citations should be studied further.

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

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

StrVCTVRE: A supervised learning method to predict the pathogenicity of human genome structural variants

Whole-genome sequencing resolves many clinical cases where standard diagnostic methods have failed. However, at least half of these cases remain unresolved after whole-genome sequencing. Structural variants (SVs; genomic variants larger than 50 base pairs) of uncertain significance are the genetic cause of a portion of these unresolved cases. As sequencing methods using long or linked reads become more accessible and SV detection algorithms improve, clinicians and researchers are gaining access to thousands of reliable SVs of unknown disease relevance. Methods to predict the pathogenicity of these SVs are required to realize the full diagnostic potential of long-read sequencing. To address this emerging need, we developed StrVCTVRE to distinguish pathogenic SVs from benign SVs that overlap exons. In a random forest classifier, we integrated features that capture gene importance, coding region, conservation, expression, and exon structure. We found that features such as expression and conservation are important but are absent from SV classification guidelines. We leveraged multiple resources to construct a size-matched training set of rare, putatively benign and pathogenic SVs. StrVCTVRE performs accurately across a wide SV size range on independent test sets, which will allow clinicians and researchers to eliminate about half of SVs from consideration while retaining a 90% sensitivity. We anticipate clinicians and researchers will use StrVCTVRE to prioritize SVs in probands where no SV is immediately compelling, empowering deeper investigation into novel SVs to resolve cases and understand new mechanisms of disease. StrVCTVRE runs rapidly and is publicly available.

Exacerbation of mild lung disorders to lethal pulmonary hypoplasia by a noncoding hypomorphic SNV in a lung-specific enhancer in trans to the frameshifting TBX4 variant

Variants involving TBX4 are associated with a wide variety of disorders, including pulmonary arterial hypertension, ischiocoxopodopatellar syndrome (ICPPS)/small patella syndrome (SPS), lethal lung developmental disorders (LLDDs) in neonates, heart defects, and prenatally lethal posterior amelia with pelvic and pulmonary hypoplasia syndrome. The objective of our study was to elucidate the wide variable phenotypic expressivity and incomplete penetrance in a three-generation family with a truncating variant in TBX4. In addition to exome and genome sequencing analyses, a candidate noncoding regulatory single nucleotide variant (SNV) within the lung-specific TBX4 enhancer was functionally tested using an in vitro luciferase reporter assay. A heterozygous frameshift variant c.1112dup (p.Pro372Serfs*14) in TBX4 was identified in patients with mild interstitial lung disease (1), bronchiolitis obliterans (1), recurrent pneumothorax (1), ICPPS/SPS (1), LLDD (2), and in unaffected individuals (4). In two deceased neonates with LLDD, we identified a noncoding SNV rs62069651-C located in trans to the mutated TBX4 allele that reduced the TBX4 promoter activity by 63% in the reporter assay. Our findings provide a functional evidence for the recently reported model of complex compound inheritance in which both TBX4 coding and in trans noncoding hypomorphic variants in the lung-specific enhancer of TBX4 contribute to LLDD.

Advances in single-cell sequencing and its application to musculoskeletal system research

In recent years, single-cell sequencing (SCS) technologies have continued to advance with improved operating procedures and reduced cost, leading to increasing practical adoption among researchers. These emerging technologies have superior abilities to analyse cell heterogeneity at a single-cell level, which have elevated multi-omics research to a higher level. In some fields of research, application of SCS has enabled many valuable discoveries, and musculoskeletal system offers typical examples. This article reviews some major scientific issues and recent advances in musculoskeletal system. In addition, combined with SCS technologies, the research of cell or tissue heterogeneity in limb development and various musculoskeletal system clinical diseases also provides new possibilities for treatment strategies. Finally, this article discusses the challenges and future development potential of SCS and recommends the direction of future applications of SCS to musculoskeletal medicine.

Diagnosis and treatment of congenital scoliosis

Background: Congenital scoliosis is caused by anomalies of the vertebra, such as hemivertebra or unsegmented bar, which result in asymmetric growth of the spine. The disruption of vertebra development during embryogenesis may be accompanied by other congenital multi-organ anomalies. The progression of the scoliotic curve may also hinder the development of other organs.Current Concepts: Hemivertebra excision and short spinal fusion have demonstrated favorable outcomes. However, the need for spinal growth and lung development has led to new treatment modalities. Growth-friendly surgeries, such as with a growing rod or vertical expandable rib-based distraction device, have demonstrated good results with curve correction while maintaining spinal growth. Although the outcome of conservative treatment for congenital scoliosis is questionable, casting may be effective as a “time-buying strategy” to delay the need for surgery.Discussion and Conclusion: It is essential to decide on a treatment plan considering the progression of the curve and growth of the spine and lungs through an individualized approach.

Clan genomics: From OMIM phenotypic traits to genes and biology

Clinical characterization of a patient phenotype has been the quintessential approach for elucidating a differential diagnosis and a hypothesis to explore a potential clinical diagnosis. This has resulted in a language of medicine and a semantic ontology, with both specialty- and subspecialty-specific lexicons, that can be challenging to translate and interpret. There is no 'Rosetta Stone' of clinical medicine such as the genetic code that can assist translation and interpretation of the language of genetics. Nevertheless, the information content embodied within a clinical diagnosis can guide management, therapeutic intervention, and potentially prognostic outlook of disease enabling anticipatory guidance for patients and families. Clinical genomics is now established firmly in medical practice. The granularity and informative content of a personal genome is immense. Yet, we are limited in our utility of much of that personal genome information by the lack of functional characterization of the overwhelming majority of computationally annotated genes in the haploid human reference genome sequence. Whereas DNA and the genetic code have provided a 'Rosetta Stone' to translate genetic variant information, clinical medicine, and clinical genomics provide the context to understand human biology and disease. A path forward will integrate deep phenotyping, such as available in a clinical synopsis in the Online Mendelian Inheritance in Man (OMIM) entries, with personal genome analyses.

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.

Disruptive NADSYN1 Variants Implicated in Congenital Vertebral Malformations

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

Hypomorphic alleles pose challenges in rare disease genomic variant interpretation

Exon skipping associated with an ATP7B intronic variant in a patient with Wilson's disease. (A) Sashimi plot visualization of aligned RNA sequencing data from proband liver tissue at ATP7B exons 14-13-12. The red track shows traditional RNA-seq data; the blue track shows RNA-seq enriched with exon capture (cDNA-cap) which achieves higher depth of protein-coding transcripts. The histogram indicates overall sequencing depth while arcs tabulate the number of junction-spanning reads supporting exon pairs. (B) The domain structure (top) and exon structure (bottom) of ATP7B. Loss of exon 13 (dashed box) would remove a transmembrane domain and disrupt the first phosphorylation domain.

Heterozygous loss of WBP11 function causes multiple congenital defects in humans and mice

The genetic causes of multiple congenital anomalies are incompletely understood. Here, we report novel heterozygous predicted loss-of-function (LoF) and predicted damaging missense variants in the WW domain binding protein 11 (WBP11) gene in seven unrelated families with a variety of overlapping congenital malformations, including cardiac, vertebral, tracheo-esophageal, renal and limb defects. WBP11 encodes a component of the spliceosome with the ability to activate pre-messenger RNA splicing. We generated a Wbp11 null allele in mouse using CRISPR-Cas9 targeting. Wbp11 homozygous null embryos die prior to E8.5, indicating that Wbp11 is essential for development. Fewer Wbp11 heterozygous null mice are found than expected due to embryonic and postnatal death. Importantly, Wbp11 heterozygous null mice are small and exhibit defects in axial skeleton, kidneys and esophagus, similar to the affected individuals, supporting the role of WBP11 haploinsufficiency in the development of congenital malformations in humans. LoF WBP11 variants should be considered as a possible cause of VACTERL association as well as isolated Klippel-Feil syndrome, renal agenesis or esophageal atresia.

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.

Saturation mutagenesis defines novel mouse models of severe spine deformity

Embryonic formation and patterning of the vertebrate spinal column requires coordination of many molecular cues. After birth, the integrity of the spine is impacted by developmental abnormalities of the skeletal, muscular and nervous systems, which may result in deformities, such as kyphosis and scoliosis. We sought to identify novel genetic mouse models of severe spine deformity by implementing in vivo skeletal radiography as part of a high-throughput saturation mutagenesis screen. We report selected examples of genetic mouse models following radiographic screening of 54,497 mice from 1275 pedigrees. An estimated 30.44% of autosomal genes harbored predicted damaging alleles examined twice or more in the homozygous state. Of the 1275 pedigrees screened, 7.4% presented with severe spine deformity developing in multiple mice, and of these, meiotic mapping implicated N-ethyl-N-nitrosourea alleles in 21% of pedigrees. Our study provides proof of concept that saturation mutagenesis is capable of discovering novel mouse models of human disease, including conditions with skeletal, neural and neuromuscular pathologies. Furthermore, we report a mouse model of skeletal disease, including severe spine deformity, caused by recessive mutation in Scube3. By integrating results with a human clinical exome database, we identified a patient with undiagnosed skeletal disease who harbored recessive mutations in SCUBE3, and we demonstrated that disease-associated mutations are associated with reduced transactivation of Smad signaling in vitro. All radiographic results and mouse models are made publicly available through the Mutagenetix online database with the goal of advancing understanding of spine development and discovering novel mouse models of human disease.

Summary: We report selected mouse models of spine deformity following mutagenesis across 30% of autosomal genes, results of which are made publicly available to advance understanding of spine development and disease.

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Copy Number Variations Analysis Identifies QPRT as a Candidate Gene Associated With Susceptibility for Solitary Functioning Kidney

Background

The lack of understanding of molecular pathologies of the solitary functioning kidney makes improving and strengthening the continuity of care between pediatric and adult nephrological patients difficult. Copy number variations (CNVs) account for a molecular cause of solitary functioning kidney, but characterization of the pathogenic genes remains challenging.

Methods

In our prospective cohort study, 99 fetuses clinically diagnosed with a solitary functioning kidney were enrolled and evaluated using chromosomal microarray analysis (CMA). The genetic drivers for the pathogenic CNVs were analyzed. We characterized QPRT localization in fetal kidneys using immunohistochemistry and its expression in adult kidneys using quantitative RT-PCR. Further, QPRT was knocked down using siRNA in human embryonic kidney (HEK293T) cells, and the cell cycle and proliferation were tested.

Results

Besides one Triple X syndrome and one Down syndrome, we identified a total of 45 CNVs out of 34 subjects. Among the 14 pathogenic CNVs, CNV 16p11.2 reached the highest number of records with the phenotype of kidney anomalies in the Decipher database. Among the 26 genes within the 16p11.2 region, as a key enzyme for nicotinamide adenine dinucleotide (NAD+) biosynthesis, QPRT was distinctly localized in renal tubules but was barely observed in renal interstitial and glomeruli in fetal kidneys. The loss of QPRT prevented cells’ efficient transition into S phase, affected cell-cycle progression, and abrogated proliferation of human embryonic kidney cells.

Conclusion

Our data suggest that QPRT is a candidate gene associated with susceptibility for solitary functioning kidney. The CNVs discovered in our study exhibit great potential for future applications in genetic counseling and pregnancy management.

Rare functional genetic variants in COL7A1, COL6A5, COL1A2 and COL5A2 frequently occur in Chiari Malformation Type 1

Chiari Malformation Type 1 (CM-1) is characterized by herniation of the cerebellar tonsils below the foramen magnum and the presence of headaches and other neurologic symptoms. Cranial bone constriction is suspected to be the most common biologic mechanism leading to CM-1. However, other mechanisms may also contribute, particularly in the presence of connective tissue disorders (CTDs), such as Ehlers Danlos Syndrome (EDS). Accumulating data suggest CM-1 with connective tissue disorders (CTD+) may have a different patho-mechanism and different genetic risk factors than CM-1 without CTDs (CTD-). To identify CM-1 genetic risk variants, we performed whole exome sequencing on a single large, multiplex family from Spain and targeted sequencing on a cohort of 186 unrelated adult, Caucasian females with CM-1. Targeted sequencing captured the coding regions of 21 CM-1 and EDS candidate genes, including two genes identified in the Spanish family. Using gene burden analysis, we compared the frequency of rare, functional variants detected in CM-1 cases versus publically available ethnically-matched controls from gnomAD. A secondary analysis compared the presence of rare variants in these genes between CTD+ and CTD- CM-1 cases. In the Spanish family, rare variants co-segregated with CM-1 in COL6A5, ADGRB3 and DST. A variant in COL7A1 was present in affected and unaffected family members. In the targeted sequencing analysis, rare variants in six genes (COL7A1, COL5A2, COL6A5, COL1A2, VEGFB, FLT1) were significantly more frequent in CM-1 cases compared to public controls. In total, 47% of CM-1 cases presented with rare variants in at least one of the four significant collagen genes and 10% of cases harbored variants in multiple significant collagen genes. Moreover, 26% of CM-1 cases presented with rare variants in the COL6A5 gene. We also identified two genes (COL7A1, COL3A1) for which the burden of rare variants differed significantly between CTD+ and CTD- CM-1 cases. A higher percentage of CTD+ patients had variants in COL7A1 compared to CTD+ patients, while CTD+ patients had fewer rare variants in COL3A1 than did CTD- patients. In summary, rare variants in several collagen genes are particularly frequent in CM-1 cases and those in COL6A5 co-segregated with CM-1 in a Spanish multiplex family. COL6A5 has been previously associated with musculoskeletal phenotypes, but this is the first association with CM-1. Our findings underscore the contribution of rare genetic variants in collagen genes to CM-1, and suggest that CM-1 in the presence and absence of CTD symptoms is driven by different genes.

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.

Altered Circulating Cell-free Mitochondrial DNA of Patients with Congenital Scoliosis

STUDY DESIGN

Case-control study.

OBJECTIVE

The aim of this study was to estimate the relationship between circulating cell-free DNA (ccf DNA) and clinical parameters of patients with congenital scoliosis (CS).

SUMMARY OF BACKGROUND DATA

CS is a complex spinal deformity characteristic of congenital vertebral malformations. Although numerous studies have centered on the etiology of CS, the cause of CS remains unclear. Previously, we reported that circulating cell-free DNA (ccf DNA) is altered in adolescent idiopathic scoliosis (AIS). However, the relationship between ccf DNA and the clinical parameters of patients with CS remains unclear.

METHODS

The plasma of peripheral blood from 35 patients with CS and 32 age-matched controls was collected for ccf DNA analysis. Quantitative PCR was used to detect ccf n-DNA and ccf mt-DNA levels, and correlation analyses between ccf n-DNA and ccf mt-DNA levels were conducted. Receiver-operating characteristic (ROC) curves were used to analyze the sensitivity and specificity of ccf n-DNA and ccf mt-DNA levels to different characteristics.

RESULTS

The plasma ccf mt-DNA levels of both ND1 and CYTC were significantly decreased in patients with CS compared with levels in controls both in total and by sex, whereas the plasma ccf n-DNA levels showed no significant difference. There is no difference in both ccf mt-DNA and ccf n-DNA between S-SDV and M-SDV according to The International Consortium for Vertebral Anomalies and Scoliosis (ICVAS) classification. The ROC curve analyses showed a reliable sensitivity and specificity of CS predicted by ccf mt-DNA levels in total but failed to distinguish different ICVAS types.

CONCLUSION

Significantly decreased plasma ccf mt-DNA levels were observed in patients with CS compared with those in controls. Although this finding has limited significance for clinical practice, it indicates that ccf mt-DNA may predict the onset or development of CS. Further studies should focus on the role of ccf mt-DNA in embryo development and whether ccf mt-DNAs could be considered as a marker for prenatal screening in development disorder like CS.Level of Evidence: 4.

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.

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.

Three case reports of patients indicating the diversity of molecular and clinical features of 16p11.2 microdeletion anomaly

BACKGROUND

16p11.2 microdeletion is a known chromosomal anomaly associated mainly with neurocognitive developmental delay, predisposition to obesity, and variable dysmorphism. Although this deletion is relatively rare among the general population, it is one of the serious known genetic aetiologies of obesity and autism spectrum disorder.

CASE PRESENTATION

This study presents three cases of deletions within the 16p11.2 region. Every child had mild variable craniofacial abnormalities, hand or foot anomalies and developmental and language delays. The first proband had obesity, epilepsy, moderate intellectual disability, aphasia, motor delay, hyperinsulinism, and café au lait spots. The second proband suffered from cardiac, pulmonary, and haematological problems. The third proband had motor and language delays, bronchial asthma, and umbilical hernia. Although each patient presented some features of the syndrome, the children differed in terms of their clinical pictures. Genetic diagnosis of 16p11.2 microdeletion syndrome was made in children at different ages based on multiplex ligation probe-dependent amplification analysis and/or microarray methods.

CONCLUSIONS

Our reports allow us to analyse and better understand the biology of 16p11.2 microdeletion throughout development. However, the variability of presented cases supports the alternate conclusion to this presented in available literature regarding 16p11.2 deletion, as we observed no direct cause-and-effect genotype/phenotype relationships. The reported cases indicate the key role of the interdisciplinary approach in 16p11.2 deletion diagnostics. The care of patients with this anomaly is based on regular health assessment and adjustment of nervous system development therapy.

Prevalence of congenital scoliosis in infants based on chest-abdomen X-ray films detected in the emergency department

PURPOSE

To investigate the prevalence of congenital scoliosis (CS) in infants based on chest-abdomen radiographs.

METHODS

A retrospective review was conducted on infants in the emergency department (ED) of a tertiary children's hospital between February 2008 and September 2019. Patients who had undergone chest-abdomen X-rays were included. All films from the enrolled patients were screened for CS. Their demographic characteristics, type, and location of the vertebral and rib anomalies, and concomitant defects of other systems were analyzed.

RESULTS

In total, 50,426 infants were enrolled; 89 (1.8‰) were diagnosed with CS, including 56 males and 33 females. There was no gender difference in CS prevalence. The visiting age of the CS patients (70 ± 98days) was significantly younger than that of the non-CS group (P < 0.05), with CS patients mainly visiting for digestive (53.9%) and respiratory symptoms (41.6%). Sixty-eight (76.4%) CS patients had main thoracic (T6-T11) vertebral malformations. Rib anomalies were documented in 27 (30.3%) patients, of which 14 had complex rib anomalies. Forty (44.9%) patients had concomitant defects of other organs, of which eight patients had two systemic abnormalities mixed. The most common extraspinal defects were imperforate anus (21, 23.6%) and congenital cardiac defects (17, 19.1%).

CONCLUSION

The prevalence of CS in infants based on chest-abdomen X-rays in the ED was 1.8‰. Both the vertebral and rib anomalies mainly affected the main thoracic region. The spine deformities in infants with concomitant defects of other organs could be identified earlier because of early-onset symptoms, which also bring out a selection bias in our analysis.

Clinical characterization of individuals with the distal 1q21.1 microdeletion

Distal 1q21.1 microdeletions have shown highly variable clinical expressivity and incomplete penetrance, with affected individuals manifesting a broad spectrum of nonspecific features. The goals of this study were to better describe the phenotypic spectrum of patients with distal 1q21.1 microdeletions and to compare the clinical features among affected individuals. We performed a retrospective chart review of 47 individuals with distal 1q21.1 microdeletions tested at a large clinical genetic testing laboratory, with most patients being clinically evaluated in the same children's hospital. Health information such as growth charts, results of imaging studies, developmental history, and progress notes were collected. Statistical analysis was performed using Fisher's exact test to compare clinical features among study subjects. Common features in our cohort include microcephaly (51.2%), seizures (29.8%), developmental delay (74.5%), failure to thrive (FTT) (68.1%), dysmorphic features (63.8%), and a variety of congenital anomalies such as cardiac abnormalities (23.4%) and genitourinary abnormalities (19.1%). Compared to prior literature, we found that seizures, brain anomalies, and FTT were more prevalent among our study cohort. Females were more likely than males to have microcephaly (p = 0.0199) and cardiac abnormalities (p = 0.0018). Based on existing genome-wide clinical testing results, at least a quarter of the cohort had additional genetic findings that may impact the phenotype of the individual. Our study represents the largest cohort of distal 1q21.1 microdeletion carriers available in the literature thus far, and it further illustrates the wide spectrum of clinical manifestations among symptomatic individuals. These results may allow for improved genetic counseling and management of affected individuals. Future studies may help to elucidate the underlying molecular mechanisms impacting the phenotypic variability observed with this microdeletion.

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.

Perturbations of genes essential for Müllerian duct and Wölffian duct development in Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser syndrome (MRKHS) is associated with congenital absence of the uterus, cervix, and the upper part of the vagina; it is a sex-limited trait. Disrupted development of the Müllerian ducts (MD)/Wölffian ducts (WD) through multifactorial mechanisms has been proposed to underlie MRKHS. In this study, exome sequencing (ES) was performed on a Chinese discovery cohort (442 affected subjects and 941 female control subjects) and a replication MRKHS cohort (150 affected subjects of mixed ethnicity from North America, South America, and Europe). Phenotypic follow-up of the female reproductive system was performed on an additional cohort of PAX8-associated congenital hypothyroidism (CH) (n = 5, Chinese). By analyzing 19 candidate genes essential for MD/WD development, we identified 12 likely gene-disrupting (LGD) variants in 7 genes: PAX8 (n = 4), BMP4 (n = 2), BMP7 (n = 2), TBX6 (n = 1), HOXA10 (n = 1), EMX2 (n = 1), and WNT9B (n = 1), while LGD variants in these genes were not detected in control samples (p = 1.27E-06). Interestingly, a sex-limited penetrance with paternal inheritance was observed in multiple families. One additional PAX8 LGD variant from the replication cohort and two missense variants from both cohorts were revealed to cause loss-of-function of the protein. From the PAX8-associated CH cohort, we identified one individual presenting a syndromic condition characterized by CH and MRKHS (CH-MRKHS). Our study demonstrates the comprehensive utilization of knowledge from developmental biology toward elucidating genetic perturbations, i.e., rare pathogenic alleles involving the same loci, contributing to human birth defects.