Homozygous Null TBX4 Mutations Lead to Posterior Amelia with Pelvic and Pulmonary Hypoplasia

Single cell RNA analysis uncovers the cell differentiation and functionalization for air breathing of frog lung

The evolution and development of vertebrate lungs have been widely studied due to their significance in terrestrial adaptation. Amphibians possess the most primitive lungs among tetrapods, underscoring their evolutionary importance in bridging the transition from aquatic to terrestrial life. However, the intricate process of cell differentiation during amphibian lung development remains poorly understood. Using single-cell RNA sequencing, we identify 13 cell types in the developing lungs of a land-dwelling frog (Microhyla fissipes). We elucidate the differentiation trajectories and mechanisms of mesenchymal cells, identifying five cell fates and their respective driver genes. Using temporal dynamics analyses, we reveal the gene expression switches of epithelial cells, which facilitate air breathing during metamorphosis. Furthermore, by integrating the published data from another amphibian and two terrestrial mammals, we illuminate both conserved and divergent cellular repertoires during the evolution of tetrapod lungs. These findings uncover the frog lung cell differentiation trajectories and functionalization for breathing in air and provide valuable insights into the cell-type evolution of vertebrate lungs.

Modelling human genetic disorders in Xenopus tropicalis

Recent progress in human disease genetics is leading to rapid advances in understanding pathobiological mechanisms. However, the sheer number of risk-conveying genetic variants being identified demands in vivo model systems that are amenable to functional analyses at scale. Here we provide a practical guide for using the diploid frog species Xenopus tropicalis to study many genes and variants to uncover conserved mechanisms of pathobiology relevant to human disease. We discuss key considerations in modelling human genetic disorders: genetic architecture, conservation, phenotyping strategy and rigour, as well as more complex topics, such as penetrance, expressivity, sex differences and current challenges in the field. As the patient-driven gene discovery field expands significantly, the cost-effective, rapid and higher throughput nature of Xenopus make it an essential member of the model organism armamentarium for understanding gene function in development and in relation to disease.

Diminished TMEM100 Expression in a Newborn With Acinar Dysplasia and a Novel TBX4 Variant: A Case Report

Acinar dysplasia (AcDys) of the lung is a rare lethal developmental disorder in neonates characterized by severe respiratory failure and pulmonary arterial hypertension refractory to treatment. Recently, abnormalities of TBX4-FGF10-FGFR2-TMEM100 signaling regulating lung development have been reported in patients with AcDys due to heterozygous single-nucleotide variants or copy-number variant deletions involving TBX4, FGF10, or FGFR2. Here, we describe a female neonate who died at 4 hours of life due to severe respiratory distress related to AcDys diagnosed by postmortem histopathologic evaluation. Genomic analyses revealed a novel deleterious heterozygous missense variant c.728A>C (p.Asn243Thr) in TBX4 that arose de novo on paternal chromosome 17. We also identified 6 candidate hypomorphic rare variants in the TBX4 enhancer in trans to TBX4 coding variant. Gene expression analyses of proband's lung tissue showed a significant reduction of TMEM100 expression with near absence of TMEM100 within the endothelium of arteries and capillaries by immunohistochemistry. These results support the pathogenicity of the detected TBX4 variant and provide further evidence that disrupted signaling between TBX4 and TMEM100 may contribute to severe lung phenotypes in humans, including AcDys.

Clinical application of next generation sequencing for Mendelian disease diagnosis in the Iranian population

Next-generation sequencing (NGS) has been proven to be one of the most powerful diagnostic tools for rare Mendelian disorders. Several studies on the clinical application of NGS in unselected cohorts of Middle Eastern patients have reported a high diagnostic yield of up to 48%, correlated with a high level of consanguinity in these populations. We evaluated the diagnostic utility of NGS-based testing across different clinical indications in 1436 patients from Iran, representing the first study of its kind in this highly consanguineous population. A total of 1075 exome sequencing and 361 targeted gene panel sequencing were performed over 8 years at a single clinical genetics laboratory, with the majority of cases tested as proband-only (91.6%). The overall diagnostic rate was 46.7%, ranging from 24% in patients with an abnormality of prenatal development to over 67% in patients with an abnormality of the skin. We identified 660 pathogenic or likely pathogenic variants, including 241 novel variants, associated with over 342 known genetic conditions. The highly consanguineous nature of this cohort led to the diagnosis of autosomal recessive disorders in the majority of patients (79.1%) and allowed us to determine the shared carrier status of couples for suspected recessive phenotypes in their deceased child(ren) when direct testing was not possible. We also highlight the observations of recessive inheritance of genes previously associated only with dominant disorders and provide an expanded genotype-phenotype spectrum for multiple less-characterized genes. We present the largest mutational spectrum of known Mendelian disease, including possible founder variants, throughout the Iranian population, which can serve as a unique resource for clinical genomic studies locally and beyond.

A novel non-recurrent CNV deletion involving TBX4 and leaving TBX2 intact causes congenital alveolar dysplasia

Congenital alveolar dysplasia (CAD) belongs to rare lethal lung developmental disorders (LLDDs) in neonates, manifesting with acute respiratory failure and pulmonary arterial hypertension refractory to treatment. The majority of CAD cases have been associated with copy-number variant (CNV) deletions at 17q23.1q23.2 or 5p12. Most CNV deletions at 17q23.1q23.2 were recurrent and encompassed two closely located genes, TBX4 and TBX2. In a few CAD cases, intragenic frameshifting deletions or single-nucleotide variants (SNVs) involved TBX4 but not TBX2. Here, we describe a male neonate who died at 27 days of life from acute respiratory failure caused by lung growth arrest along the spectrum of CAD confirmed by histopathological assessment. Trio-based genome sequencing revealed in the proband a novel non-recurrent ~1.07 Mb heterozygous CNV deletion at 17q23.2, encompassing TBX4 that arose de novo on the paternal chromosome. This is the first report of a larger-sized CNV deletion in a CAD patient involving TBX4 and leaving TBX2 intact. Our results, together with previous reports, indicate that perturbations of TBX4, rather than TBX2, cause severe lung phenotypes in humans.

Long-Term Effect of TBX4 Germline Mutation on Pulmonary Clinico-Histopathologic Phenotype

Tbx4 protein, expressed in mesenchyme of the developing lung, contributes to airway branching and distal lung growth. An association between pediatric onset of pulmonary arterial hypertension (PAH) and genetic variations coding for the T-box transcription factor 4 gene (TBX4) has been increasingly recognized. Tbx4-related PAH onset has a bimodal age distribution, including severe to lethal PAH in newborns and later onset PAH. We present an autopsy study of a 24-year-old male with a heterozygous TBX4 variant, who developed pulmonary arterial hypertension at age 12 years. This unique case highlights the complex pulmonary histopathology leading to lethal cardiopulmonary failure in the setting of TBX4 mutation.

Molecular Function and Contribution of TBX4 in Development and Disease

Over the past decade, recognition of the profound impact of the TBX4 (T-box 4) gene, which encodes a member of the evolutionarily conserved family of T-box-containing transcription factors, on respiratory diseases has emerged. The developmental importance of TBX4 is emphasized by the association of TBX4 variants with congenital disorders involving respiratory and skeletal structures; however, the exact role of TBX4 in human development remains incompletely understood. Here, we discuss the developmental, tissue-specific, and pathological TBX4 functions identified through human and animal studies and review the published TBX4 variants resulting in variable disease phenotypes. We also outline future research directions to fill the gaps in our understanding of TBX4 function and of how TBX4 disruption affects development.

Evolutionary genetics of flipper forelimb and hindlimb loss from limb development-related genes in cetaceans

Background

Cetacean hindlimbs were lost and their forelimb changed into flippers characterized by webbed digits and hyperphalangy, thus allowing them to adapt to a completely aquatic environment. However, the underlying molecular mechanism behind cetacean limb development remains poorly understood.

Results

In the present study, we explored the evolution of 16 limb-related genes and their cis-regulatory elements in cetaceans and compared them with that of other mammals. TBX5, a forelimb specific expression gene, was identified to have been under accelerated evolution in the ancestral branches of cetaceans. In addition, 32 cetacean-specific changes were examined in the SHH signaling network (SHH, PTCH1, TBX5, BMPs and SMO), within which mutations could yield webbed digits or an additional phalange. These findings thus suggest that the SHH signaling network regulates cetacean flipper formation. By contrast, the regulatory activity of the SHH gene enhancer—ZRS in cetaceans—was significantly lower than in mice, which is consistent with the cessation of SHH gene expression in the hindlimb bud during cetacean embryonic development. It was suggested that the decreased SHH activity regulated by enhancer ZRS might be one of the reasons for hindlimb degeneration in cetaceans. Interestingly, a parallel / convergent site (D42G) and a rapidly evolving CNE were identified in marine mammals in FGF10 and GREM1, respectively, and shown to be essential to restrict limb bud size; this is molecular evidence explaining the convergence of flipper-forelimb and shortening or degeneration of hindlimbs in marine mammals.

Conclusions

We did evolutionary analyses of 16 limb-related genes and their cis-regulatory elements in cetaceans and compared them with those of other mammals to provide novel insights into the molecular basis of flipper forelimb and hindlimb loss in cetaceans.

Identification and functional evaluation of GRIA1 missense and truncation variants in individuals with ID: An emerging neurodevelopmental syndrome

GRIA1 encodes the GluA1 subunit of α-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) receptors, which are ligand-gated ion channels that act as excitatory receptors for the neurotransmitter L-glutamate (Glu). AMPA receptors (AMPARs) are homo- or heteromeric protein complexes with four subunits, each encoded by different genes, GRIA1 to GRIA4. Although GluA1-containing AMPARs have a crucial role in brain function, the human phenotype associated with deleterious GRIA1 sequence variants has not been established. Subjects with de novo missense and nonsense GRIA1 variants were identified through international collaboration. Detailed phenotypic and genetic assessments of the subjects were carried out and the pathogenicity of the variants was evaluated in vitro to characterize changes in AMPAR function and expression. In addition, two Xenopus gria1 CRISPR-Cas9 F₀ models were established to characterize the in vivo consequences. Seven unrelated individuals with rare GRIA1 variants were identified. One individual carried a homozygous nonsense variant (p.Arg377Ter), and six had heterozygous missense variations (p.Arg345Gln, p.Ala636Thr, p.Ile627Thr, and p.Gly745Asp), of which the p.Ala636Thr variant was recurrent in three individuals. The cohort revealed subjects to have a recurrent neurodevelopmental disorder mostly affecting cognition and speech. Functional evaluation of major GluA1-containing AMPAR subtypes carrying the GRIA1 variant mutations showed that three of the four missense variants profoundly perturb receptor function. The homozygous stop-gain variant completely destroys the expression of GluA1-containing AMPARs. The Xenopus gria1 models show transient motor deficits, an intermittent seizure phenotype, and a significant impairment to working memory in mutants. These data support a developmental disorder caused by both heterozygous and homozygous variants in GRIA1 affecting AMPAR function.

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.

Tbx4 function during hindlimb development reveals a mechanism that explains the origins of proximal limb defects

We dissect genetically a gene regulatory network that involves the transcription factors Tbx4, Pitx1 and Isl1 acting cooperatively to establish the hindlimb bud, and identify key differences in the pathways that initiate formation of the hindlimb and forelimb. Using live image analysis of murine limb mesenchyme cells undergoing chondrogenesis in micromass culture, we distinguish a series of changes in cellular behaviours and cohesiveness that are required for chondrogenic precursors to undergo differentiation. Furthermore, we provide evidence that the proximal hindlimb defects observed in Tbx4 mutant mice result from a failure in the early differentiation step of chondroprogenitors into chondrocytes, providing an explanation for the origins of proximally biased limb defects.

Ttc30a affects tubulin modifications in a model for ciliary chondrodysplasia with polycystic kidney disease

Cilia are tubulin-based cellular appendages, and their dysfunction has been linked to a variety of genetic diseases. Ciliary chondrodysplasia is one such condition that can co-occur with cystic kidney disease and other organ manifestations. We modeled skeletal ciliopathies by mutating two established disease genes in Xenopus tropicalis frogs. Bioinformatic analysis identified ttc30a as a ciliopathy network component, and targeting it replicated skeletal malformations and renal cysts as seen in patients and the amphibian models. A loss of Ttc30a affected cilia by altering posttranslational tubulin modifications. Our findings identify TTC30A/B as a component of ciliary segmentation essential for cartilage differentiation and renal tubulogenesis. These findings may lead to novel therapeutic targets in treating ciliary skeletopathies and cystic kidney disease.

Skeletal ciliopathies (e.g., Jeune syndrome, short rib polydactyly syndrome, and Sensenbrenner syndrome) are frequently associated with nephronophthisis-like cystic kidney disease and other organ manifestations. Despite recent progress in genetic mapping of causative loci, a common molecular mechanism of cartilage defects and cystic kidneys has remained elusive. Targeting two ciliary chondrodysplasia loci (ift80 and ift172) by CRISPR/Cas9 mutagenesis, we established models for skeletal ciliopathies in Xenopus tropicalis. Froglets exhibited severe limb deformities, polydactyly, and cystic kidneys, closely matching the phenotype of affected patients. A data mining–based in silico screen found ttc30a to be related to known skeletal ciliopathy genes. CRISPR/Cas9 targeting replicated limb malformations and renal cysts identical to the models of established disease genes. Loss of Ttc30a impaired embryonic renal excretion and ciliogenesis because of altered posttranslational tubulin acetylation, glycylation, and defective axoneme compartmentalization. Ttc30a/b transcripts are enriched in chondrocytes and osteocytes of single-cell RNA-sequenced embryonic mouse limbs. We identify TTC30A/B as an essential node in the network of ciliary chondrodysplasia and nephronophthisis-like disease proteins and suggest that tubulin modifications and cilia segmentation contribute to skeletal and renal ciliopathy manifestations of ciliopathies in a cell type–specific manner. These findings have implications for potential therapeutic strategies.

Comparative genomic analysis of a naturally born serpentized pig reveals putative mutations related to limb and bone development

BACKGROUND

It is believed that natural selection acts on the phenotypical changes caused by mutations. Phenotypically, from fishes to amphibians to reptiles, the emergence of limbs greatly facilitates the landing of ancient vertebrates, but the causal mutations and evolutionary trajectory of this process remain unclear.

RESULTS

We serendipitously obtained a pig of limbless phenotype. Mutations specific to this handicapped pig were identified using genome re-sequencing and comparative genomic analysis. We narrowed down the causal mutations to particular chromosomes and even several candidate genes and sites, such like a mutation-containing codon in gene BMP7 (bone morphogenetic protein) which was conserved in mammals but variable in lower vertebrates.

CONCLUSIONS

We parsed the limbless-related mutations in the light of evolution. The limbless pig shows phenocopy of the clades before legs were evolved. Our findings might help deduce the emergence of limbs during vertebrate evolution and should be appealing to the broad community of human genetics and evolutionary biology.

Mutation profiling of a limbless pig reveals genome-wide regulation of RNA processing related to bone development

Mutation is the basis of phenotypic changes and serves as the source of natural selection. The development of limbs has been the milestone in vertebrate evolution. Several limb and bone-related genes were verified experimentally, but other indirect and regulatory factors of limb development remained untested, especially very few cases were observed in natural environment. We report a naturally born serpentized pig without hindlimbs. Whole genome sequencing followed by comparative genomic analysis revealed multiple interesting patterns on the handicapped pig-specific mutations. Although the bone-related genes are not directly subjected to mutations, other regulatory factors such as the RNA deaminase genes Adar are damaged in the handicapped pig, leading to the abolished A-to-I deamination in many functional, conserved genes as well as the bone-related genes. This is a precious case that the limbless phenotype is observed in naturally born non-model organisms. Our study broadened the generality of the limbless phenotype across mammals and extended the regulation of hindlimb development to other non-bone-related genes. Our knowledge of limb and bone-related mutations and regulation would also contribute to human genetics.

Biallelic variants in COPB1 cause a novel, severe intellectual disability syndrome with cataracts and variable microcephaly

BACKGROUND

Coat protein complex 1 (COPI) is integral in the sorting and retrograde trafficking of proteins and lipids from the Golgi apparatus to the endoplasmic reticulum (ER). In recent years, coat proteins have been implicated in human diseases known collectively as "coatopathies".

METHODS

Whole exome or genome sequencing of two families with a neuro-developmental syndrome, variable microcephaly and cataracts revealed biallelic variants in COPB1, which encodes the beta-subunit of COPI (β-COP). To investigate Family 1's splice donor site variant, we undertook patient blood RNA studies and CRISPR/Cas9 modelling of this variant in a homologous region of the Xenopus tropicalis genome. To investigate Family 2's missense variant, we studied cellular phenotypes of human retinal epithelium and embryonic kidney cell lines transfected with a COPB1 expression vector into which we had introduced Family 2's mutation.

RESULTS

We present a new recessive coatopathy typified by severe developmental delay and cataracts and variable microcephaly. A homozygous splice donor site variant in Family 1 results in two aberrant transcripts, one of which causes skipping of exon 8 in COPB1 pre-mRNA, and a 36 amino acid in-frame deletion, resulting in the loss of a motif at a small interaction interface between β-COP and β'-COP. Xenopus tropicalis animals with a homologous mutation, introduced by CRISPR/Cas9 genome editing, recapitulate features of the human syndrome including microcephaly and cataracts. In vitro modelling of the COPB1 c.1651T>G p.Phe551Val variant in Family 2 identifies defective Golgi to ER recycling of this mutant β-COP, with the mutant protein being retarded in the Golgi.

CONCLUSIONS

This adds to the growing body of evidence that COPI subunits are essential in brain development and human health and underlines the utility of exome and genome sequencing coupled with Xenopus tropicalis CRISPR/Cas modelling for the identification and characterisation of novel rare disease genes.

Maximizing CRISPR/Cas9 phenotype penetrance applying predictive modeling of editing outcomes in Xenopus and zebrafish embryos

CRISPR/Cas9 genome editing has revolutionized functional genomics in vertebrates. However, CRISPR/Cas9 edited F0 animals too often demonstrate variable phenotypic penetrance due to the mosaic nature of editing outcomes after double strand break (DSB) repair. Even with high efficiency levels of genome editing, phenotypes may be obscured by proportional presence of in-frame mutations that still produce functional protein. Recently, studies in cell culture systems have shown that the nature of CRISPR/Cas9-mediated mutations can be dependent on local sequence context and can be predicted by computational methods. Here, we demonstrate that similar approaches can be used to forecast CRISPR/Cas9 gene editing outcomes in Xenopus tropicalis, Xenopus laevis, and zebrafish. We show that a publicly available neural network previously trained in mouse embryonic stem cell cultures (InDelphi-mESC) is able to accurately predict CRISPR/Cas9 gene editing outcomes in early vertebrate embryos. Our observations can have direct implications for experiment design, allowing the selection of guide RNAs with predicted repair outcome signatures enriched towards frameshift mutations, allowing maximization of CRISPR/Cas9 phenotype penetrance in the F0 generation.

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.