TBX6 null variants and a common hypomorphic allele in congenital scoliosis

Comparative proteomics analysis for identifying the lipid metabolism related pathways in patients with Klippel-Feil syndrome

Background

Klippel-Feil syndrome (KFS) represents the rare and complex deformity characterized by congenital defects in the formation or segmentation of the cervical vertebrae. There is a wide gap in understanding the detailed mechanisms of KFS because of its rarity, heterogeneity, small pedigrees, and the broad spectrum of anomalies.

Methods

We recruited eight patients of Chinese Han ethnicity with KFS, five patients with congenital scoliosis (CS) who presented with congenital fusion of the thoracic or lumbar spine and without known syndrome or cervical deformity, and seven healthy controls. Proteomic analysis by data-independent acquisition (DIA) was performed to identify the differential proteome among the three matched groups and the data were analyzed by bioinformatics tools including Gene Ontology (GO) categories and Ingenuity Pathway Analysis (IPA) database, to explore differentially abundant proteins (DAPs) and canonical pathways involved in the pathogenesis of KFS.

Results

A total of 49 DAPs were detected between KFS patients and the controls, and moreover, 192 DAPs were identified between patients with KFS and patients with CS. Fifteen DAPs that were common in both comparisons were considered as candidate biomarkers for KFS, including membrane primary amine oxidase, noelin, galectin-3-binding protein, cadherin-5, glyceraldehyde-3-phosphate dehydrogenase, peroxiredoxin-1, CD109 antigen, and eight immunoglobulins. Furthermore, the same significant canonical pathways of LXR/RXR activation and FXR/RXR activation were observed in both comparisons. Seven of DAPs were apolipoproteins related to these pathways that are involved in lipid metabolism.

Conclusions

This study provides the first proteomic profile for understanding the pathogenesis and identifying predictive biomarkers of KFS. We detected 15 DAPs that were common in both comparisons as candidate predictive biomarkers of KFS. The lipid metabolism-related canonical pathways of LXR/RXR and FXR/RXR activation together with seven differentially abundant apolipoproteins may play significant roles in the etiology of KFS and provide possible pathogenesis correlation between KFS and CS.

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

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

Molecular Diagnosis of Craniosynostosis Using Targeted Next-Generation Sequencing

BACKGROUND

Genetic factors play an important role in the pathogenesis of craniosynostosis (CRS). However, the molecular diagnosis of CRS in clinical practice is limited because of its heterogeneous etiology.

OBJECTIVE

To investigate the genomic landscape of CRS in a Korean cohort and also to establish a practical diagnostic workflow by applying targeted panel sequencing.

METHODS

We designed a customized panel covering 34 CRS-related genes using in-solution hybrid capture method. We enrolled 110 unrelated Korean patients with CRS, including 40 syndromic and 70 nonsyndromic cases. A diagnostic pipeline was established by combining in-depth clinical reviews and multiple bioinformatics tools for analyzing single-nucleotide variants (SNV)s and copy number variants (CNV)s.

RESULTS

The diagnostic yield of the targeted panel was 30.0% (33/110). Twenty-five patients (22.7%) had causal genetic variations resulting from SNVs or indels in 9 target genes (TWIST1, FGFR3, TCF12, ERF, FGFR2, ALPL, EFNB1, FBN1, and SKI, in order of frequency). CNV analysis identified 8 (7.3%) additional patients with chromosomal abnormalities involving 1p32.3p31.3, 7p21.1, 10q26, 15q21.3, 16p11.2, and 17p13.3 regions; these cases mostly presented with syndromic clinical features.

CONCLUSION

The present study shows the wide genomic landscape of CRS, revealing various genetic factors for CRS pathogenesis. In addition, the results demonstrate that an efficient diagnostic workup using target panel sequencing provides great clinical utility in the molecular diagnosis of CRS.

Understanding the clinical manifestations of 16p11.2 deletion syndrome: a series of developmental case reports in children

BACKGROUND

Copy number variants (CNVs) are genetic rearrangements, such as deletions and duplications, which result in a deviation from the normal number of copies of a given gene segment. CNVs are implicated in many neuropsychiatric disorders. Deletions of the human chromosomal region 16p11.2 are one of the most common genetic linkages to autism spectrum disorders (ASD). However, ASD is not the only presenting feature, and many patients with 16p11.2 deletions present with a variable clinical spectrum.

METHODS

To better understand the nature and presentation of the syndrome throughout development, we present three different, unrelated clinical cases of children with 16p11.2 deletion and provide a detailed description of their clinical manifestations.

RESULTS

Cognitive and motor impairments were characteristic of all three patients with 16p11.2 deletion, despite the differences in the extent and clinical presentation of impairment. Two patients had a clinical diagnosis of ASD and one showed several ASD traits. In addition, two patients also had severe speech and language impairments, which is in line with previous reports on 16p11.2 phenotypes. Although epilepsy and obesity have been frequently associated with 16p11.2 deletion, only one patient had a diagnosis of epilepsy and none of the three cases were obese.

CONCLUSION

This variation in clinical phenotype renders correct clinical interpretation and diagnosis challenging. Therefore, it is critical to elucidate the variable clinical phenotypes of rare CNVs, including 16p11.2 deletions, to help guide clinical monitoring and counselling of patients and families.

Rare genetic causes of complex kidney and urological diseases

Although often considered a single-entity, chronic kidney disease (CKD) comprises many pathophysiologically distinct disorders that result in persistently abnormal kidney structure and/or function, and encompass both monogenic and polygenic aetiologies. Rare inherited forms of CKD frequently span diverse phenotypes, reflecting genetic phenomena including pleiotropy, incomplete penetrance and variable expressivity. Use of chromosomal microarray and massively parallel sequencing technologies has revealed that genomic disorders and monogenic aetiologies contribute meaningfully to seemingly complex forms of CKD across different clinically defined subgroups and are characterized by high genetic and phenotypic heterogeneity. Investigations of prevalent genomic disorders in CKD have integrated genetic, bioinformatic and functional studies to pinpoint the genetic drivers underlying their renal and extra-renal manifestations, revealing both monogenic and polygenic mechanisms. Similarly, massively parallel sequencing-based analyses have identified gene- and allele-level variation that contribute to the clinically diverse phenotypes observed for many monogenic forms of nephropathy. Genome-wide sequencing studies suggest that dual genetic diagnoses are found in at least 5% of patients in whom a genetic cause of disease is identified, highlighting the fact that complex phenotypes can also arise from multilocus variation. A multifaceted approach that incorporates genetic and phenotypic data from large, diverse cohorts will help to elucidate the complex relationships between genotype and phenotype for different forms of CKD, supporting personalized medicine for individuals with kidney disease.

Cooperation between ETS transcription factor ETV1 and histone demethylase JMJD1A in colorectal cancer

ETS variant 1 (ETV1) is an oncogenic transcription factor. However, its role in colorectal cancer has remained understudied. The present study demonstrated that ETV1 downregulation led to reduced HCT116 colorectal cancer cell growth and clonogenic activity. Furthermore, the ETV1 mRNA levels were enhanced in colorectal tumors and were associated with disease severity. In addition, ETV1 directly bound to Jumonji C domain-containing (JMJD) 1A, a histone demethylase known to promote colon cancer. ETV1 and JMJD1A, but not a catalytically inactive mutant thereof, cooperated in inducing the matrix metalloproteinase (MMP)1 gene promoter that was similar to the cooperation between ETV1 and another histone demethylase, JMJD2A. RNA-sequencing revealed multiple potential ETV1 target genes in HCT116 cells, including the FOXQ1 and TBX6 transcription factor genes. Moreover, JMJD1A co-regulated FOXQ1 and other ETV1 target genes, but not TBX6, whereas JMJD2A downregulation had no impact on FOXQ1 as well as TBX6 transcription. Accordingly, the FOXQ1 gene promoter was stimulated by ETV1 and JMJD1A in a cooperative manner, and both ETV1 and JMJD1A bound to the FOXQ1 promoter. Notably, the overexpression of FOXQ1 partially reversed the growth inhibitory effects of ETV1 ablation on HCT116 cells, whereas TBX6 impaired HCT116 cell growth and may thereby dampen the oncogenic activity of ETV1. The latter also revealed for the first time, to the best of our knowledge, a potential tumor suppressive function of TBX6. Taken together, the present study uncovered a ETV1/JMJD1A-FOXQ1 axis that may drive colorectal tumorigenesis.

Clinical genomics and contextualizing genome variation in the diagnostic laboratory

INTRODUCTION

The human genome contains the instructions for the development and biological homeostasis of the human organism and the genetic transmission of traits. Genome variation in human populations is the basis of evolution; individual or personal genomes vary tremendously, making each of us truly unique.

AREAS COVERED

Assaying this individual variation using genomic technologies has many applications in clinical medicine, from elucidating the biology of disease to designing strategies to ameliorate perturbations from homeostasis. Detecting pathogenic rare variation in a genome may provide a molecular diagnosis that can be informative for patient management and family healthcare.

EXPERT OPINION

Despite the increasing clinical use of unbiased genomic testing, including chromosome microarray analysis (CMA) with array comparative genomic hybridization (aCGH) or SNP arrays, clinical exome sequencing (cES), and whole-genome sequencing (WGS), to survey genome-wide for molecular aberrations, clinical acumen paired with an understanding of the limitations of each testing type will be needed to achieve molecular diagnoses. Potential opportunities for improving case solved rates, functionally annotating the majority of genes in the human genome, and further understanding genetic contributions to disease will empower clinical genomics and the precision medicine initiative.

New Insight of Circular RNAs in Human Musculoskeletal Diseases

Circular RNAs (circRNAs), a novel group of noncoding RNAs, are present in most eukaryotic cells. Different from messenger RNAs, circRNAs have a covalently closed single-stranded stable structure and often act in cell type and tissue-specific manners, indicating that they can be used as biomarkers. With the advance of high-throughput RNA sequencing technology and bioinformatics, a large number of circRNAs have been identified in association with musculoskeletal diseases, but the functions of most circRNAs have not been clarified. circRNAs regulate biological processes by adsorbing microRNA as "sponges," binding to proteins, acting as transcriptional regulators, and participating in translation of proteins. In this study, we discuss the latest understanding of biogenesis and gene regulatory mechanisms of circRNAs with special emphasis on new targets for musculoskeletal disease diagnosis and clinical treatment.

Muscle defects due to perturbed somite segmentation contribute to late adult scoliosis

Scoliosis is an abnormal bending of the body axis. Truncated vertebrae or a debilitated ability to control the musculature in the back can cause this condition, but in most cases the causative reason for scoliosis is unknown (idiopathic). Using mutants for somite clock genes with mild defects in the vertebral column, we here show that early defects in somitogenesis are not overcome during development and have long lasting and profound consequences for muscle fiber organization, structure and whole muscle volume. These mutants present only mild alterations in the vertebral column, and muscle shortcomings are uncoupled from skeletal defects. None of the mutants presents an overt musculoskeletal phenotype at larval or early adult stages, presumably due to compensatory growth mechanisms. Scoliosis becomes only apparent during aging. We conclude that adult degenerative scoliosis is due to disturbed crosstalk between vertebrae and muscles during early development, resulting in subsequent adult muscle weakness and bending of the body axis.

Applications and grants of National Natural Scientific Foundation of China's General Program in abnormalities and diseases of locomotor system: a ten-year review

Translational medicine is a branch of medical research that attempts to establish a more direct link between basic research and clinical practice. The improvement of clinical management cannot be separated from the progress and breakthrough made in basic biomedical research. As one of the main funding agencies for basic research in China, National Natural Scientific Foundation of China (NSFC) plays an active role in promotion of the development of basic medical research and translational medicine. General Program mainly supports scientists to carry out basic research on bottom-up based topics within the funding scope of NSFC to conduct innovative research and promote a balanced, coordinated and sustained development of all disciplines. Musculoskeletal injury and disease is a major threat to human health worldwide and has become the major cause of disability in China. In this review, we aim to further promote the clinical transformation by summarizing the trends and hotspots of basic research in this field based on the analysis of the data of General Program in abnormalities and diseases of locomotor system over the last decade, including analyses of the number of applications received and projects funded, the distribution of the keywords in applications and grants, distribution of applications and grants of General Program in the secondary application code of H06 and in the host institution.

Cost-effectiveness analysis of using the TBX6-associated congenital scoliosis risk score (TACScore) in genetic diagnosis of congenital scoliosis

Background

We previously reported a novel clinically distinguishable subtype of congenital scoliosis (CS), namely, TBX6-associated congenital scoliosis (TACS). We further developed the TBX6-associated CS risk score (TACScore), a multivariate phenotype-based model to predict TACS according to the patient’s clinical manifestations. In this study, we aimed to evaluate whether using the TACScore as a screening method prior to performing whole-exome sequencing (WES) is more cost-effective than using WES as the first-line genetic test for CS.

Methods

We retrospectively collected the molecular data of 416 CS patients in the Deciphering disorders Involving Scoliosis and COmorbidities (DISCO) study. A decision tree was constructed to estimate the cost and the diagnostic time required for the two alternative strategies (TACScore versus WES). Bootstrapping simulations and sensitivity analyses were performed to examine the distributions and robustness of the estimates. The economic evaluation considered both the health care payer and the personal budget perspectives.

Results

From the health care payer perspective, the strategy of using the TACScore as the primary screening method resulted in an average cost of $1074.2 (95%CI: $1044.8 to $1103.5) and an average diagnostic duration of 38.7d (95%CI: 37.8d to 39.6d) to obtain a molecular diagnosis for each patient. In contrast, the corresponding values were $1169.6 (95%CI: $1166.9 to $1172.2) and 41.4d (95%CI: 41.1d to 41.7d) taking WES as the first-line test (P < 0.001). From the personal budget perspective, patients who were predicted to be positive by the TACScore received a result with an average cost of $715.1 (95%CI: $594.5 to $835.7) and an average diagnostic duration of 30.4d (95%CI: 26.3d to 34.6d). Comparatively, the strategy of WES as the first-line test was estimated to have significantly longer diagnostic time with an average of 44.0d (95%CI: 43.2d to 44.9d), and more expensive with an average of $1193.4 (95%CI: $1185.5 to $1201.3) (P < 0.001). In 100% of the bootstrapping simulations, the TACScore strategy was significantly less costly and more time-saving than WES. The sensitivity analyses revealed that the TACScore strategy remained cost-effective even when the cost per WES decreased to $8.8.

Conclusions

This retrospective study provides clinicians with economic evidence to integrate the TACScore into clinical practice. The TACScore can be considered a cost-effective tool when it serves as a screening test prior to performing WES.

Neurodevelopmental trajectory and modifiers of 16p11.2 microdeletion: A follow-up study of four Chinese children carriers

Background

Neurodevelopmental disorders (NDDs) are a group of disorders with high genetic and phenotypic heterogeneities. The 16p11.2 microdeletion has been implicated as an important genetic risk factor for NDDs.

Methods

Multiple genetic tests were used to detect the 16p11.2 microdeletion from 918 Chinese children with NDDs. Targeted sequencing of genes in the 16p11.2 interval was performed in all carriers of the 16p11.2 microdeletion, and whole‐genome expression profiling analysis was performed for the patient carriers and normal carriers in their intra‐family.

Results

Three patients carrying the 16p11.2 microdeletion were screened out, indicating a frequency of 0.33% for the 16p11.2 microdeletion in this cohort. We reviewed the neurodevelopmental trajectories of the 16p11.2 microdeletion carriers from childhood to puberty and confirmed that this microdeletion was associated with abnormal neurodevelopment, with varied neurodevelopmental phenotypes. A differential PRRT2 genotype (rs10204, T>C) was identified between patients and normal carriers of the 16p11.2 microdeletion. Moreover, the determination of differential whole‐genome expression profiling demonstrated the destruction of the top‐ranked network in neurogenesis and accounted for observation of abnormal neurodevelopmental phenotypes in the 16p11.2 microdeletion carriers.

Conclusions

We have provided the frequency of the 16p11.2 microdeletion in a Chinese pediatric NDD cohort with a variable NDD phenotype from childhood to puberty, which is useful for Chinese geneticists/pediatricians to conduct the 16p11.2 microdeletion testing in children with NDDs.

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.

A Recurrent Rare SOX9 Variant (M469V) is Associated with Congenital Vertebral Malformations

OBJECTIVE

The genetic variations contributed to a substantial proportion of congenital vertebral malformations (CVM). SOX9 gene, a member of the SOX gene family, has been implicated in CVM. To study the SOX9 mutation in CVM patients is of great significance to explain the pathogenesis of scoliosis (the clinical manifestation of CVM) and to explore the pathogenesis of SOX9-related skeletal deformities.

METHODS

A total of 50 singleton patients with CVM were included in this study. Exome Sequencing (ES) was performed on all the patients. The recurrent candidate variant of SOX9 gene was validated by Sanger sequencing. Luciferase assay was performed to investigate the functional changes of this variant.

RESULTS

A recurrent rare heterozygous missense variant in SOX9 gene (NM_000346.3: c.1405A>G, p.M469V) which had not been reported previously was identified in three CVM patients who had the clinical findings of congenital scoliosis without deformities in other systems. This variant was absent from our in-house database and it was predicted to be deleterious (CADD = 24.5). The luciferase assay demonstrated that transactivation capacity of the mutated SOX9 protein was significantly lower than that of the wild-type for the two luciferase reporters (p = 0.0202, p = 0.0082, respectively).

CONCLUSION

This SOX9 mutation (p.M469V) may contribute to CVM without other systematic deformity, which provides important implications and better understanding of phenotypic variability in SOX9-related skeletal deformities.

Whole genome sequencing of familial isolated oesophagus atresia uncover shared structural variants

Background

Oesophageal atresia (OA) is a life-threatening developmental defect characterized by a lost continuity between the upper and lower oesophagus. The most common form is a distal connection between the trachea and the oesophagus, i.e. a tracheoesophageal fistula (TEF). The condition may be part of a syndrome or occurs as an isolated feature. The recurrence risk in affected families is increased compared to the population-based incidence suggesting contributing genetic factors.

Methods

To gain insight into gene variants and genes associated with isolated OA we conducted whole genome sequencing on samples from three families with recurrent cases affected by congenital and isolated TEF.

Results

We identified a combination of single nucleotide variants (SNVs), splice site variants (SSV) and structural variants (SV) annotated to altogether 100 coding genes in the six affected individuals.

Conclusion

This study highlights rare SVs among candidate gene variants in our individuals with OA and provides a gene framework for further investigations of genetic factors behind this malformation.

Notochordal Signals Establish Phylogenetic Identity of the Teleost Spine

The spine is a defining feature of the vertebrate body plan. However, broad differences in vertebral structures and morphogenetic strategies occur across vertebrate groups, clouding the homology between their developmental programs. Analysis of a zebrafish mutant, spondo, whose spine is dysmorphic, prompted us to reconstruct paleontological evidence, highlighting specific transitions during teleost spine evolution. Interestingly, the spondo mutant recapitulates characteristics present in basal fishes, not found in extant teleosts. Further analysis of the mutation implicated the teleost-specific notochord protein, Calymmin, as a key regulator of spine patterning in zebrafish. The mutation in cmn results in loss of notochord sheath segmentation, altering osteoblast migration to the developing spine, and increasing sensitivity to somitogenesis defects associated with congenital scoliosis in amniotes. These data suggest that signals from the notochord define the evolutionary identity of the spine and demonstrate how simple shifts in development can revert traits canalized for about 250 million years.

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

BACKGROUND

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

METHODS AND RESULTS

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

CONCLUSION

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

High methylation of lysine acetyltransferase 6B is associated with the Cobb angle in patients with congenital scoliosis

Background

The etiology of congenital scoliosis (CS) is complex and uncertain. Abnormal DNA methylation affects the growth and development of spinal development. In this study, we investigated the role of DNA methylation in CS.

Methods

The target region DNA methylation level in the peripheral blood of patients with CS was analyzed. Through in-depth analysis, genes closely related to the growth and development of the vertebra were identified. EdU staining was performed to verify the role of differentially expressed genes in chondrocyte proliferation.

Results

The hypermethylated KAT6B gene was observed in patients with CS, and was positively correlated with the Cobb angle. KAT6B was primarily expressed on chondrocytes. The promoter of KAT6B in CS patients was hypermethylated, and its expression was significantly reduced. Further mechanistic studies revealed that EZH2 mediated trimethylation of lysine 27 on histone H3 of the KAT6B promoter. Overexpression of KAT6B in CS-derived primary chondrocytes can significantly promote chondrocyte proliferation, which may be related to activation of the RUNX2/Wnt/β-catenin signaling pathway.

Conclusion

Epigenetic modification of KAT6B may be a cause of CS. If similar epigenetic modification abnormalities can be detected through maternal liquid biopsy screening, they may provide useful biomarkers for early screening and diagnosis of CS.

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

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

Progress of orthopaedic research in China over the last decade

Objective

To summarize the representative scientific achievements in the past decade, and discuss the future challenges and directions for orthopaedic research in China.

Methods

In this review, we used the data provided by National Natural Science Foundation of China (NSFC) for analysis.

Results

Over the last decade, NSFC has initiated various research programs with a total funding of over 1149 million RMB to support orthopaedic exploration. Under the strong support of NSFC, great progresses have been made in basic research, talent training, platform construction and the clinical translation in the field of orthopaedics in China.

Conclusion

In general, since the establishment of the Department of Health Sciences of NSFC 10 years ago, both the amount of funding and the scale of researchers in the field of orthopaedic research have increased substantially. Despite of several shortcomings in orthopaedic research, with continuous support from NSFC both in funding and in policy, we believe that the orthopaedic research in China will surely make steady and significant progress.

The translational potential of this article

This article summarizes the representative scientific achievements in the past decade and puts forward the future challenges and directions for orthopaedic research in China.

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.

Diagnostic yield and clinical impact of exome sequencing in early-onset scoliosis (EOS)

BACKGROUND

Early-onset scoliosis (EOS), defined by an onset age of scoliosis less than 10 years, conveys significant health risk to affected children. Identification of the molecular aetiology underlying patients with EOS could provide valuable information for both clinical management and prenatal screening.

METHODS

In this study, we consecutively recruited a cohort of 447 Chinese patients with operative EOS. We performed exome sequencing (ES) screening on these individuals and their available family members (totaling 670 subjects). Another cohort of 13 patients with idiopathic early-onset scoliosis (IEOS) from the USA who underwent ES was also recruited.

RESULTS

After ES data processing and variant interpretation, we detected molecular diagnostic variants in 92 out of 447 (20.6%) Chinese patients with EOS, including 8 patients with molecular confirmation of their clinical diagnosis and 84 patients with molecular diagnoses of previously unrecognised diseases underlying scoliosis. One out of 13 patients with IEOS from the US cohort was molecularly diagnosed. The age at presentation, the number of organ systems involved and the Cobb angle were the three top features predictive of a molecular diagnosis.

CONCLUSION

ES enabled the molecular diagnosis/classification of patients with EOS. Specific clinical features/feature pairs are able to indicate the likelihood of gaining a molecular diagnosis through ES.

Functional biology of the Steel syndrome founder allele and evidence for clan genomics derivation of COL27A1 pathogenic alleles worldwide

Previously we reported the identification of a homozygous COL27A1 (c.2089G>C; p.Gly697Arg) missense variant and proposed it as a founder allele in Puerto Rico segregating with Steel syndrome (STLS, MIM #615155); a rare osteochondrodysplasia characterized by short stature, congenital bilateral hip dysplasia, carpal coalitions, and scoliosis. We now report segregation of this variant in five probands from the initial clinical report defining the syndrome and an additional family of Puerto Rican descent with multiple affected adult individuals. We modeled the orthologous variant in murine Col27a1 and found it recapitulates some of the major Steel syndrome associated skeletal features including reduced body length, scoliosis, and a more rounded skull shape. Characterization of the in vivo murine model shows abnormal collagen deposition in the extracellular matrix and disorganization of the proliferative zone of the growth plate. We report additional COL27A1 pathogenic variant alleles identified in unrelated consanguineous Turkish kindreds suggesting Clan Genomics and identity-by-descent homozygosity contributing to disease in this population. The hypothesis that carrier states for this autosomal recessive osteochondrodysplasia may contribute to common complex traits is further explored in a large clinical population cohort. Our findings augment our understanding of COL27A1 biology and its role in skeletal development; and expand the functional allelic architecture in this gene underlying both rare and common disease phenotypes.

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

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

A de novo 2.2 Mb recurrent 17q23.1q23.2 deletion unmasks novel putative regulatory non-coding SNVs associated with lethal lung hypoplasia and pulmonary hypertension: a case report

BACKGROUND

Application of whole genome sequencing (WGS) enables identification of non-coding variants that play a phenotype-modifying role and are undetectable by exome sequencing. Recently, non-coding regulatory single nucleotide variants (SNVs) have been reported in patients with lethal lung developmental disorders (LLDDs) or congenital scoliosis with recurrent copy-number variant (CNV) deletions at 17q23.1q23.2 or 16p11.2, respectively.

CASE PRESENTATION

Here, we report a deceased newborn with pulmonary hypertension and pulmonary interstitial emphysema with features suggestive of pulmonary hypoplasia, resulting in respiratory failure and neonatal death soon after birth. Using the array comparative genomic hybridization and WGS, two heterozygous recurrent CNV deletions: ~ 2.2 Mb on 17q23.1q23.2, involving TBX4, and ~ 600 kb on 16p11.2, involving TBX6, that both arose de novo on maternal chromosomes were identified. In the predicted lung-specific enhancer upstream to TBX4, we have detected seven novel putative regulatory non-coding SNVs that were absent in 13 control individuals with the overlapping deletions but without any structural lung anomalies.

CONCLUSIONS

Our findings further support a recently reported model of complex compound inheritance of LLDD in which both non-coding and coding heterozygous TBX4 variants contribute to the lung phenotype. In addition, this is the first report of a patient with combined de novo heterozygous recurrent 17q23.1q23.2 and 16p11.2 CNV deletions.

The publication trends and hot spots of scoliosis research from 2009 to 2018: a 10-year bibliometric analysis

Background

This study aims to quantitatively and qualitatively investigate the trends in scoliosis research and evaluate research hotspots using bibliometric analysis.

Methods

All relevant publications on scoliosis from the period from 2009 to 2018 were extracted from the Web of Science and PubMed databases. Publication trends were analyzed using an Online analysis platform of literature metrology, Bibliographic Item Co-occurrence Matrix Builder (BICOMB), and CiteSpace software. Hotspots were analyzed and visualized using the gCLUTO software package.

Results

A total of 7,445 scoliosis research publications dated between 2009 and 2018 were found. The spine was the most popular journal in this field during this period. The United States maintained a top position in global scoliosis research throughout the 10 years and has had a pivotal influence, followed by China and Canada. Among all institutions, the University of California, San Francisco, was a leader in research collaboration. At the same time, Professors Yong Qiu and Lawrence G. Lenke made great achievements in scoliosis research. We analyzed the major Medical Subject Headings (MeSH) terms/MeSH subheadings and identified eight hotspots in scoliosis research.

Conclusions

We summarized the publication information of scoliosis-related literature in the 10 years from 2009 to 2018, including country and institution of origin, authors, and publication journal. We analyzed former research hotspots in the field of scoliosis and predicted future areas of interest. The development of various new orthopedic plants, artificial intelligence diagnosis, and genetic research will be future hotspots in scoliosis research.

Genetic Burden Contributing to Extremely Low or High Bone Mineral Density in a Senior Male Population From the Osteoporotic Fractures in Men Study (MrOS)

Worldwide, one in five men aged over 50 years will experience osteoporosis or a clinical bone fracture, with a greater fracture-related mortality rate than women. However, the genetic etiology of osteoporosis in men is still poorly understood. We aimed to identify the genetic variants and candidate genes associated with extremely low or high BMD for a better understanding of the biology underlying low bone density that may point to potential therapeutic targets for increasing bone mass. Subjects from the Osteoporotic Fractures in Men Study (MrOS) cohort were evaluated by age and BMI-adjusted total hip BMD. Those with BMD values 3 SDs away from the mean were selected and the remaining individuals whose adjusted BMD ranked at the highest or lowest 100 were included. Men with the lowest adjusted BMD (N = 98) and highest adjusted BMD (N = 110) were chosen for exome sequencing. Controls (N = 82) were men of Northern and Western European descent from the US Utah population of the 1000 Genomes Project. Fisher's exact test was performed to compare low- or high-BMD subjects with controls for single-gene associations. Additionally, sets of candidate genes causative of heritable disorders of connective tissue, including osteogenesis imperfecta (OI) and Ehlers-Danlos syndrome (EDS), were grouped for multigene and mutation burden analyses. No single-gene associations with rare variants were found for either the low BMD group (33 genes) or high BMD group (18 genes). In the group of OI genes, we detected a significant threefold increased accumulation of rare variants in low-BMD subjects compared with controls (p = 0.009). Additionally, genes associated with EDS had a twofold increased frequency in low-BMD subjects compared with controls (p = 0.03). These findings reveal a rare variant burden in OI and EDS disease genes at low BMD, which suggests a potential gene-panel approach to screen for multivariant associations in larger cohorts. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Dstyk mutation leads to congenital scoliosis-like vertebral malformations in zebrafish via dysregulated mTORC1/TFEB pathway

Congenital scoliosis (CS) is a complex genetic disorder characterized by vertebral malformations. The precise etiology of CS is not fully defined. Here, we identify that mutation in dual serine/threonine and tyrosine protein kinase (dstyk) lead to CS-like vertebral malformations in zebrafish. We demonstrate that the scoliosis in dstyk mutants is related to the wavy and malformed notochord sheath formation and abnormal axial skeleton segmentation due to dysregulated biogenesis of notochord vacuoles and notochord function. Further studies show that DSTYK is located in late endosomal/lysosomal compartments and is involved in the lysosome biogenesis in mammalian cells. Dstyk knockdown inhibits notochord vacuole and lysosome biogenesis through mTORC1-dependent repression of TFEB nuclear translocation. Inhibition of mTORC1 activity can rescue the defect in notochord vacuole biogenesis and scoliosis in dstyk mutants. Together, our findings reveal a key role of DSTYK in notochord vacuole biogenesis, notochord morphogenesis and spine development through mTORC1/TFEB pathway.

Congenital scoliosis is a complex genetic disorder characterized by vertebral malformation. Here, the authors demonstrate that loss of dstyk leads to scoliosis in zebrafish due to dysregulated biogenesis of notochord vacuoles and that DSTYK is required for lysosome biogenesis through mTORC1 regulation.

Development of coupling controlled polymerizations by adapter-ligation in mate-pair sequencing for detection of various genomic variants in one single assay

The diversity of disease presentations warrants one single assay for detection and delineation of various genomic disorders. Herein, we describe a gel-free and biotin-capture-free mate-pair method through coupling Controlled Polymerizations by Adapter-Ligation (CP-AL). We first demonstrated the feasibility and ease-of-use in monitoring DNA nick translation and primer extension by limiting the nucleotide input. By coupling these two controlled polymerizations by a reported non-conventional adapter-ligation reaction 3′ branch ligation, we evidenced that CP-AL significantly increased DNA circularization efficiency (by 4-fold) and was applicable for different sequencing methods but at a faction of current cost. Its advantages were further demonstrated by fully elimination of small-insert-contaminated (by 39.3-fold) with a ∼50% increment of physical coverage, and producing uniform genome/exome coverage and the lowest chimeric rate. It achieved single-nucleotide variants detection with sensitivity and specificity up to 97.3 and 99.7%, respectively, compared with data from small-insert libraries. In addition, this method can provide a comprehensive delineation of structural rearrangements, evidenced by a potential diagnosis in a patient with oligo-atheno-terato-spermia. Moreover, it enables accurate mutation identification by integration of genomic variants from different aberration types. Overall, it provides a potential single-integrated solution for detecting various genomic variants, facilitating a genetic diagnosis in human diseases.

Genetic and molecular mechanism for distinct clinical phenotypes conveyed by allelic truncating mutations implicated in FBN1

BACKGROUND

The molecular and genetic mechanisms by which different single nucleotide variant alleles in specific genes, or at the same genetic locus, cause distinct disease phenotypes often remain unclear. Allelic truncating mutations of FBN1 could cause either classical Marfan syndrome (MFS) or a more complicated phenotype associated with Marfanoid-progeroid-lipodystrophy syndrome (MPLS).

METHODS

We investigated a small cohort, encompassing two classical MFS and one MPLS subjects from China, whose clinical presentation included scoliosis potentially requiring surgical intervention. Targeted next generation sequencing was performed on all the participants. We analyzed the molecular diagnosis, clinical features, and the potential molecular mechanism involved in the MPLS subject in our cohort.

RESULTS

We report a novel de novo FBN1 mutation for the first Chinese subject with MPLS, a more complicated fibrillinopathy, and two subjects with more classical MFS. We further predict that the MPLS truncating mutation, and others previously reported, is prone to escape the nonsense-mediated decay (NMD), while MFS mutations are predicted to be subjected to NMD. Also, the MPLS mutation occurs within the glucogenic hormone asprosin domain of FBN1. In vitro experiments showed that the single MPLS mutation p.Glu2759Cysfs*9 appears to perturb proper FBN1 protein aggregation as compared with the classical MFS mutation p.Tyr2596Thrfs*86. Both mutations appear to upregulate SMAD2 phosphorylation in vitro.

CONCLUSION

We provide direct evidence that a dominant-negative interaction of FBN1 potentially explains the complex MPLS phenotypes through genetic and functional analysis. Our study expands the mutation spectrum of FBN1 and highlights the potential molecular mechanism for MPLS.

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

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

Mutations in ANKLE2, a ZIKA Virus Target, Disrupt an Asymmetric Cell Division Pathway in Drosophila Neuroblasts to Cause Microcephaly

The apical Par complex, which contains atypical protein kinase C (aPKC), Bazooka (Par-3), and Par-6, is required for establishing polarity during asymmetric division of neuroblasts in Drosophila, and its activity depends on L(2)gl. We show that loss of Ankle2, a protein associated with microcephaly in humans and known to interact with Zika protein NS4A, reduces brain volume in flies and impacts the function of the Par complex. Reducing Ankle2 levels disrupts endoplasmic reticulum (ER) and nuclear envelope morphology, releasing the kinase Ballchen-VRK1 into the cytosol. These defects are associated with reduced phosphorylation of aPKC, disruption of Par-complex localization, and spindle alignment defects. Importantly, removal of one copy of ballchen or l(2)gl suppresses Ankle2 mutant phenotypes and restores viability and brain size. Human mutational studies implicate the above-mentioned genes in microcephaly and motor neuron disease. We suggest that NS4A, ANKLE2, VRK1, and LLGL1 define a pathway impinging on asymmetric determinants of neural stem cell division.

Identification of novel FBN1 variations implicated in congenital scoliosis

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

Nablus syndrome: Easy to diagnose yet difficult to solve

Nablus syndrome was first described by the late Ahmad Teebi in 2000, and 13 individuals have been reported to date. Nablus syndrome can be clinically diagnosed based on striking facial features, including tight glistening skin with reduced facial expression, blepharophimosis, telecanthus, bulky nasal tip, abnormal external ear architecture, upswept frontal hairline, and sparse eyebrows. However, the precise genetic etiology for this rare condition remains elusive. Comparative microarray analyses of individuals with Nablus syndrome (including two mother-son pairs) reveal an overlapping 8q22.1 microdeletion, with a minimal critical region of 1.84 Mb (94.43-96.27 Mb). Whereas this deletion is present in all affected individuals, 13 individuals without Nablus syndrome (including two mother-child pairs) also have the 8q22.1 microdeletion that partially or fully overlaps the minimal critical region. Thus, the 8q22.1 microdeletion is necessary but not sufficient to cause the clinical features characteristic of Nablus syndrome. We discuss possible explanations for Nablus syndrome, including one-locus, two-locus, epigenetic, and environmental mechanisms. We performed exome sequencing for five individuals with Nablus syndrome. Although we failed to identify any deleterious rare coding variants in the critical region that were shared between individuals, we did identify one common SNP in an intronic region that was shared. Clearly, unraveling the genetic mechanism(s) of Nablus syndrome will require additional investigation, including genomic and RNA sequencing of a larger cohort of affected individuals. If successful, it will provide important insights into fundamental concepts such as variable expressivity, incomplete penetrance, and complex disease relevant to both Mendelian and non-Mendelian disorders.

A Genocentric Approach to Discovery of Mendelian Disorders

The advent of inexpensive, clinical exome sequencing (ES) has led to the accumulation of genetic data from thousands of samples from individuals affected with a wide range of diseases, but for whom the underlying genetic and molecular etiology of their clinical phenotype remains unknown. In many cases, detailed phenotypes are unavailable or poorly recorded and there is little family history to guide study. To accelerate discovery, we integrated ES data from 18,696 individuals referred for suspected Mendelian disease, together with relatives, in an Apache Hadoop data lake (Hadoop Architecture Lake of Exomes [HARLEE]) and implemented a genocentric analysis that rapidly identified 154 genes harboring variants suspected to cause Mendelian disorders. The approach did not rely on case-specific phenotypic classifications but was driven by optimization of gene- and variant-level filter parameters utilizing historical Mendelian disease-gene association discovery data. Variants in 19 of the 154 candidate genes were subsequently reported as causative of a Mendelian trait and additional data support the association of all other candidate genes with disease endpoints.

Association of rare non-coding SNVs in the lung-specific FOXF1 enhancer with a mitigation of the lethal ACDMPV phenotype

Haploinsufficiency of FOXF1 causes alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a lethal neonatal lung developmental disorder. We describe two similar heterozygous CNV deletions involving the FOXF1 enhancer and re-analyze FOXF1 missense mutation, all associated with an unexpectedly mitigated disease phenotype. In one case, the deletion of the maternal allele of the FOXF1 enhancer caused pulmonary hypertension and histopathologically diagnosed MPV without the typical ACD features. In the second case, the deletion of the paternal enhancer resulted in ACDMPV rather than the expected neonatal lethality. In both cases, FOXF1 expression in lung tissue was higher than usually seen or expected in patients with similar deletions, suggesting an increased activity of the remaining allele of the enhancer. Sequencing of these alleles revealed two rare SNVs, rs150502618-A and rs79301423-T, mapping to the partially overlapping binding sites for TFAP2s and CTCF in the core region of the enhancer. Moreover, in a family with three histopathologically-diagnosed ACDMPV siblings whose missense FOXF1 mutation was inherited from the healthy non-mosaic carrier mother, we have identified a rare SNV rs28571077-A within 2-kb of the above-mentioned non-coding SNVs in the FOXF1 enhancer in the mother, that was absent in the affected newborns and 13 unrelated ACDMPV patients with CNV deletions of this genomic region. Based on the low population frequencies of these three variants, their absence in ACDMPV patients, the results of reporter assay, RNAi and EMSA experiments, and in silico predictions, we propose that the described SNVs might have acted on FOXF1 enhancer as hypermorphs.

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

STUDY OBJECTIVE

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

DESIGN

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

SETTING

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

PARTICIPANTS AND MAIN OUTCOME MEASURES

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

RESULTS

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

CONCLUSION

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

Diagnostic interpretation of genetic studies in patients with primary immunodeficiency diseases: A working group report of the Primary Immunodeficiency Diseases Committee of the American Academy of Allergy, Asthma & Immunology

Genetic testing has become an integral component of the diagnostic evaluation of patients with suspected primary immunodeficiency diseases. Results of genetic testing can have a profound effect on clinical management decisions. Therefore clinical providers must demonstrate proficiency in interpreting genetic data. Because of the need for increased knowledge regarding this practice, the American Academy of Allergy, Asthma & Immunology Primary Immunodeficiency Diseases Committee established a work group that reviewed and summarized information concerning appropriate methods, tools, and resources for evaluating variants identified by genetic testing. Strengths and limitations of tests frequently ordered by clinicians were examined. Summary statements and tables were then developed to guide the interpretation process. Finally, the need for research and collaboration was emphasized. Greater understanding of these important concepts will improve the diagnosis and management of patients with suspected primary immunodeficiency diseases.

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

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

Germline 16p11.2 Microdeletion Predisposes to Neuroblastoma

Neuroblastoma is a cancer of the developing sympathetic nervous system. It is diagnosed in 600-700 children per year in the United States and accounts for 12% of pediatric cancer deaths. Despite recent advances in our understanding of this malignancy's complex genetic architecture, the contribution of rare germline variants remains undefined. Here, we conducted a genome-wide analysis of large (>500 kb), rare (<1%) germline copy number variants (CNVs) in two independent, multi-ethnic cohorts totaling 5,585 children with neuroblastoma and 23,505 cancer-free control children. We identified a 550-kb deletion on chromosome 16p11.2 significantly enriched in neuroblastoma cases (0.39% of cases and 0.03% of controls; p = 3.34 × 10^-9). Notably, this CNV corresponds to a known microdeletion syndrome that affects approximately one in 3,000 children and confers risk for diverse developmental phenotypes including autism spectrum disorder and other neurodevelopmental disorders. The CNV had a substantial impact on neuroblastoma risk, with an odds ratio of 13.9 (95% confidence interval = 5.8-33.4). The association remained significant when we restricted our analysis to individuals of European ancestry in order to mitigate potential confounding by population stratification (0.42% of cases and 0.03% of controls; p = 4.10 × 10^-8). We used whole-genome sequencing (WGS) to validate the deletion in paired germline and tumor DNA from 12 cases. Finally, WGS of four parent-child trios revealed that the deletion primarily arose de novo without maternal or paternal bias. This finding expands the clinical phenotypes associated with 16p11.2 microdeletion syndrome to include cancer, and it suggests that disruption of the 16p11.2 region may dysregulate neurodevelopmental pathways that influence both neurological phenotypes and neuroblastoma.

Mendelian Gene Discovery: Fast and Furious with No End in Sight

Gene discovery for Mendelian conditions (MCs) offers a direct path to understanding genome function. Approaches based on next-generation sequencing applied at scale have dramatically accelerated gene discovery and transformed genetic medicine. Finding the genetic basis of ∼6,000-13,000 MCs yet to be delineated will require both technical and computational innovation, but will rely to a larger extent on meaningful data sharing.

Prenatal Diagnosis of Fetuses With Increased Nuchal Translucency by Genome Sequencing Analysis

Background: Increased nuchal translucency (NT) is an important biomarker associated with increased risk of fetal structural anomalies. It is known to be contributed by a wide range of genetic etiologies from single-nucleotide variants to those affecting millions of base pairs. Currently, prenatal diagnosis is routinely performed by karyotyping and chromosomal microarray analysis (CMA); however, both of them have limited resolution. The diversity of the genetic etiologies warrants an integrated assay such as genome sequencing (GS) for comprehensive detection of genomic variants. Herein, we aim to evaluate the feasibility of applying GS in prenatal diagnosis for the fetuses with increased NT.

Methods: We retrospectively applied GS (> 30-fold) for fetuses with increased NT (≥3.5 mm) who underwent routine prenatal diagnosis. Detection of single-nucleotide variants, copy number variants, and structural rearrangements was performed simultaneously, and the results were integrated for interpretation in accordance with the guidelines of the American College of Medical Genetics and Genomics. Pathogenic or likely pathogenic (P/LP) variants were selected for validation and parental confirmation, when available.

Results: Overall, 50 fetuses were enrolled, including 34 cases with isolated increased NT and 16 cases with other fetal structural malformations. Routine CMA and karyotyping reported eight P/LP CNVs, yielding a diagnostic rate of 16.0% (8/50). In comparison, GS provided a twofold increase in diagnostic yield (32.0%, 16/50), including one mosaic turner syndrome, eight cases with microdeletions/microduplications, and seven cases with P/LP point mutations. Moreover, GS identified two cryptic insertions and two inversions. Follow-up study further demonstrated the potential pathogenicity of an apparently balanced insertion that disrupted an OMIM autosomal dominant disease-causing gene at the insertion site.

Conclusions: Our study demonstrates that applying GS in fetuses with increased NT can comprehensively detect and delineate the various genomic variants that are causative to the diseases. Importantly, prenatal diagnosis by GS doubled the diagnostic yield compared with routine protocols. Given a comparable turnaround time and less DNA required, our study provides strong evidence to facilitate GS in prenatal diagnosis, particularly in fetuses with increased NT.

Genome-wide association study identifies 14 previously unreported susceptibility loci for adolescent idiopathic scoliosis in Japanese

Adolescent idiopathic scoliosis (AIS) is the most common pediatric spinal deformity. Several AIS susceptibility loci have been identified; however, they could explain only a small proportion of AIS heritability. To identify additional AIS susceptibility loci, we conduct a meta-analysis of the three genome-wide association studies consisting of 79,211 Japanese individuals. We identify 20 loci significantly associated with AIS, including 14 previously not reported loci. These loci explain 4.6% of the phenotypic variance of AIS. We find 21 cis-expression quantitative trait loci-associated genes in seven of the fourteen loci. By a female meta-analysis, we identify additional three significant loci. We also find significant genetic correlations of AIS with body mass index and uric acid. The cell-type specificity analyses show the significant heritability enrichment for AIS in multiple cell-type groups, suggesting the heterogeneity of etiology and pathogenesis of AIS. Our findings provide insights into etiology and pathogenesis of AIS.

Adolescent idiopathic scoliosis (AIS) is a common pediatric disease leading to spinal deformities. Here, the authors report GWAS followed by genome-wide meta-analysis in up to 79,211 Japanese individuals, identifying 20 genetic loci for AIS, 14 of which were previously unreported, and perform in vitro validation for rs1978060.

Genome sequencing and implications for rare disorders

The practice of genomic medicine stands to revolutionize our approach to medical care, and to realize this goal will require discovery of the relationship between rare variation at each of the ~ 20,000 protein-coding genes and their consequent impact on individual health and expression of Mendelian disease. The step-wise evolution of broad-based, genome-wide cytogenetic and molecular genomic testing approaches (karyotyping, chromosomal microarray [CMA], exome sequencing [ES]) has driven much of the rare disease discovery to this point, with genome sequencing representing the newest member of this team. Each step has brought increased sensitivity to interrogate individual genomic variation in an unbiased method that does not require clinical prediction of the locus or loci involved. Notably, each step has also brought unique limitations in variant detection, for example, the low sensitivity of ES for detection of triploidy, and of CMA for detection of copy neutral structural variants. The utility of genome sequencing (GS) as a clinical molecular diagnostic test, and the increased sensitivity afforded by addition of long-read sequencing or other -omics technologies such as RNAseq or metabolomics, are not yet fully explored, though recent work supports improved sensitivity of variant detection, at least in a subset of cases. The utility of GS will also rely upon further elucidation of the complexities of genetic and allelic heterogeneity, multilocus rare variation, and the impact of rare and common variation at a locus, as well as advances in functional annotation of identified variants. Much discovery remains to be done before the potential utility of GS is fully appreciated.