Dysregulation of miRNA expression and excitation in MEF2C autism patient hiPSC-neurons and cerebral organoids.

MEF2C is a critical transcription factor in neurodevelopment, whose loss-of-function mutation in humans results in MEF2C haploinsufficiency syndrome (MHS), a severe form of autism spectrum disorder (ASD)/intellectual disability (ID). Despite prior animal studies of MEF2C heterozygosity to mimic MHS, MHS-specific mutations have not been investigated previously, particularly in a human context as hiPSCs afford. Here, for the first time, we use patient hiPSC-derived cerebrocortical neurons and cerebral organoids to characterize MHS deficits. Unexpectedly, we found that decreased neurogenesis was accompanied by activation of a micro-(mi)RNA-mediated gliogenesis pathway. We also demonstrate network-level hyperexcitability in MHS neurons, as evidenced by excessive synaptic and extrasynaptic activity contributing to excitatory/inhibitory (E/I) imbalance. Notably, the predominantly extrasynaptic (e)NMDA receptor antagonist, NitroSynapsin, corrects this aberrant electrical activity associated with abnormal phenotypes. During neurodevelopment, MEF2C regulates many ASD-associated gene networks, suggesting that treatment of MHS deficits may possibly help other forms of ASD as well.

Epilepsy as a Novel Phenotype of BPTF-Related Disorders.

BACKGROUND: Neurodevelopmental disorder with dysmorphic facies and distal limb anomalies (NEDDFL) is associated to BPTF gene haploinsufficiency. Epilepsy was not included in the initial descriptions of NEDDFL, but emerging evidence indicates that epileptic seizures occur in some affected individuals. This study aims to investigate the electroclinical epilepsy features in individuals with NEDDFL. METHODS: We enrolled individuals with BPTF-related seizures or interictal epileptiform discharges (IEDs) on electroencephalography (EEG). Demographic, clinical, genetic, raw EEG, and neuroimaging data as well as response to antiseizure medication were assessed. RESULTS: We studied 11 individuals with a null variant in BPTF, including five previously unpublished ones. Median age at last observation was 9 years (range: 4 to 43 years). Eight individuals had epilepsy, one had a single unprovoked seizure, and two showed IEDs only. Key features included (1) early childhood epilepsy onset (median 4 years, range: 10 months to 7 years), (2) well-organized EEG background (all cases) and brief bursts of spikes and slow waves (50% of individuals), and (3) developmental delay preceding seizure onset. Spectrum of epilepsy severity varied from drug-resistant epilepsy (27%) to isolated IEDs without seizures (18%). Levetiracetam was widely used and reduced seizure frequency in 67% of the cases. CONCLUSIONS: Our study provides the first characterization of BPTF-related epilepsy. Early-childhood-onset epilepsy occurs in 19% of subjects, all presenting with a well-organized EEG background associated with generalized interictal epileptiform abnormalities in half of these cases. Drug resistance is rare.

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.

A Small De Novo CNV Deletion of the Paternal Copy of FOXF1, Leaving lncRNA FENDRR Intact, Provides Insight into Their Bidirectional Promoter Region.

Pathogenic single-nucleotide variants (SNVs) and copy-number variant (CNV) deletions involving the FOXF1 transcription factor gene or CNV deletions of its distant lung-specific enhancer are responsible for alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a rarely diagnosed lethal lung developmental disorder in neonates. In contrast to SNVs within FOXF1 and CNV deletions involving only the FOXF1 enhancer, larger-sized deletions involving FOXF1 and the adjacent, oppositely oriented lncRNA gene FENDRR have additionally been associated with hypoplastic left heart syndrome and single umbilical artery (SUA). Here, in an ACDMPV infant without any congenital heart defect or SUA, we identified a small 5 kb CNV deletion that removed the paternal allele of FOXF1 and its promoter, leaving FENDRR and its promoter intact. Reporter assay in the IMR-90 fetal cell line implied that the deletion may indeed not have significantly affected FENDRR expression. It also showed a polarization of the FOXF1-FENDRR inter-promoter region consisting of its ability to increase the transcription of FENDRR but not FOXF1. Interestingly, this transcription-stimulating activity was suppressed in the presence of the FOXF1 promoter. Our data shed more light on the interactions between neighboring promoters of FOXF1-FENDRR and possibly other divergently transcribed mRNA-lncRNA gene pairs.

Further refinement of the differentially methylated distant lung-specific FOXF1 enhancer in a neonate with alveolar capillary dysplasia.

Heterozygous SNVs or CNV deletions involving the FOXF1 gene, or its distant enhancer, are causative for 80-90% of cases of alveolar capillary dysplasia with misalignment of pulmonary veins. Recently, we proposed bimodal structure and parental functional dimorphism of the lung-specific FOXF1 enhancer, with Unit 1 having higher activity on the paternal chr16 and Unit 2 on the maternal chr16. Here, we describe a novel unusually sized pathogenic de novo copy-number variant deletion involving a portion of the FOXF1 enhancer on maternal chr16 that implies narrowing Unit 2 to an essential ~ 9-kb segment. Using a restrictase-based assay, we found that this enhancer segment is weakly methylated at ApT adenine, with about twice the frequency of methylation on the maternal versus paternal chr16. Our data provide further insight into the FOXF1 enhancer structure and function.

PhaseDancer: a novel targeted assembler of segmental duplications unravels the complexity of the human chromosome 2 fusion going from 48 to 46 chromosomes in hominin evolution.

Resolving complex genomic regions rich in segmental duplications (SDs) is challenging due to the high error rate of long-read sequencing. Here, we describe a targeted approach with a novel genome assembler PhaseDancer that extends SD-rich regions of interest iteratively. We validate its robustness and efficiency using a golden-standard set of human BAC clones and in silico-generated SDs with predefined evolutionary scenarios. PhaseDancer enables extension of the incomplete complex SD-rich subtelomeric regions of Great Ape chromosomes orthologous to the human chromosome 2 (HSA2) fusion site, informing a model of HSA2 formation and unravelling the evolution of human and Great Ape genomes.

Histologic features and decreased lung FOXF1 gene expression in severe bronchopulmonary dysplasia without a genetic diagnosis of alveolar capillary dysplasia.

We report the case of a preterm infant who died at 10 months of age with severe bronchopulmonary dysplasia (sBPD) with refractory pulmonary hypertension and respiratory failure who had striking histologic features compatible with the diagnosis of alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) but without genetic confirmation of the diagnosis. We further demonstrate dramatic reductions in lung FOXF1 and TMEM100 content in sBPD, suggesting common mechanistic links between ACDMPV and sBPD with impaired FOXF1 signaling.

Single Cell Multiomics Identifies Cells and Genetic Networks Underlying Alveolar Capillary Dysplasia.

Rationale: Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a lethal developmental disorder of lung morphogenesis caused by insufficiency of FOXF1 (forkhead box F1) transcription factor function. The cellular and transcriptional mechanisms by which FOXF1 deficiency disrupts human lung formation are unknown. Objectives: To identify cell types, gene networks, and cell-cell interactions underlying the pathogenesis of ACDMPV. Methods: We used single-nucleus RNA and assay for transposase-accessible chromatin sequencing, immunofluorescence confocal microscopy, and RNA in situ hybridization to identify cell types and molecular networks influenced by FOXF1 in ACDMPV lungs. Measurements and Main Results: Pathogenic single-nucleotide variants and copy-number variant deletions involving the FOXF1 gene locus in all subjects with ACDMPV (n = 6) were accompanied by marked changes in lung structure, including deficient alveolar development and a paucity of pulmonary microvasculature. Single-nucleus RNA and assay for transposase-accessible chromatin sequencing identified alterations in cell number and gene expression in endothelial cells (ECs), pericytes, fibroblasts, and epithelial cells in ACDMPV lungs. Distinct cell-autonomous roles for FOXF1 in capillary ECs and pericytes were identified. Pathogenic variants involving the FOXF1 gene locus disrupt gene expression in EC progenitors, inhibiting the differentiation or survival of capillary 2 ECs and cell-cell interactions necessary for both pulmonary vasculogenesis and alveolar type 1 cell differentiation. Loss of the pulmonary microvasculature was associated with increased VEGFA (vascular endothelial growth factor A) signaling and marked expansion of systemic bronchial ECs expressing COL15A1 (collagen type XV α 1 chain). Conclusions: Distinct FOXF1 gene regulatory networks were identified in subsets of pulmonary endothelial and fibroblast progenitors, providing both cellular and molecular targets for the development of therapies for ACDMPV and other diffuse lung diseases of infancy.

SNV/indel hypermutator phenotype in biallelic RAD51C variant: Fanconi anemia.

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C. Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5 and 37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the hypothesis of a rare phenomenon of hypermutation. Eight-hundred-thirty apparent de novo SNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component having an average VAF of 25% (range: 18.7-37.41%) (n = 446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range 9.55-18.69%) (n = 384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA. We propose that other FA groups should be investigated for genome-wide de novo variants.

SNV/indel hypermutator phenotype in biallelic RAD51C variant - Fanconi anemia.

We previously reported a fetus with Fanconi anemia (FA), complementation group O due to compound heterozygous variants involving RAD51C . Interestingly, the trio exome sequencing analysis also detected eight apparent de novo mosaic variants with variant allele fraction (VAF) ranging between 11.5%-37%. Here, using whole genome sequencing and a 'home-brew' variant filtering pipeline and DeepMosaic module, we investigated the number and signature of de novo heterozygous and mosaic variants and the rare phenomenon of hypermutation. Eight-hundred-thirty apparent dnSNVs and 21 de novo indels had VAFs below 37.41% and were considered postzygotic somatic mosaic variants. The VAFs showed a bimodal distribution, with one component with an average VAF of 25% (range: 18.7-37.41%) (n=446), representing potential postzygotic first mitotic events, and the other component with an average VAF of 12.5% (range: 9.55-18.69%) (n=384), describing potential second mitotic events. No increased rate of CNV formation was observed. The mutational pattern analysis for somatic single base substitution showed SBS40, SBS5, and SBS3 as the top recognized signatures. SBS3 is a known signature associated with homologous recombination-based DNA damage repair error. Our data demonstrate that biallelic RAD51C variants show evidence for defective genomic DNA damage repair and thereby result in a hypermutator phenotype with the accumulation of postzygotic de novo mutations, at least in the prenatal period. This 'genome hypermutator phenomenon' might contribute to the observed hematological manifestations and the predisposition to tumors in patients with FA, and pregnancy loss in general. We propose that other FA groups should be investigated for genome-wide de novo variants.

Prenatal Detection of a FOXF1 Deletion in a Fetus with ACDMPV and Hydronephrosis.

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a lethal lung developmental disorder caused by the arrest of fetal lung formation, resulting in neonatal death due to acute respiratory failure and pulmonary arterial hypertension. Heterozygous single-nucleotide variants or copy-number variant (CNV) deletions involving the FOXF1 gene and/or its lung-specific enhancer are found in the vast majority of ACDMPV patients. ACDMPV is often accompanied by extrapulmonary malformations, including the gastrointestinal, cardiac, or genitourinary systems. Thus far, most of the described ACDMPV patients have been diagnosed post mortem, based on histologic evaluation of the lung tissue and/or genetic testing. Here, we report a case of a prenatally detected de novo CNV deletion (~0.74 Mb) involving the FOXF1 gene in a fetus with ACDMPV and hydronephrosis. Since ACDMPV is challenging to detect by ultrasound examination, the more widespread implementation of prenatal genetic testing can facilitate early diagnosis, improve appropriate genetic counselling, and further management.

Cloud-native distributed genomic pileup operations.

MOTIVATION: Pileup analysis is a building block of many bioinformatics pipelines, including variant calling and genotyping. This step tends to become a bottleneck of the entire assay since the straightforward pileup implementations involve processing of all base calls from all alignments sequentially. On the other hand, a distributed version of the algorithm faces the intrinsic challenge of splitting reads-oriented file formats into self-contained partitions to avoid costly data exchange between computational nodes. RESULTS: Here, we present a scalable, distributed and efficient implementation of a pileup algorithm that is suitable for deploying in cloud computing environments. In particular, we implemented: (i) our custom data-partitioning algorithm optimized to work with the alignment reads, (ii) a novel and unique approach to process alignment events from sequencing reads using the MD tags, (iii) the source code micro-optimizations for recurrent operations, and (iv) a modular structure of the algorithm. We have proven that our novel approach consistently and significantly outperforms other state-of-the-art distributed tools in terms of execution time (up to 6.5× faster) and memory usage (up to 2× less), resulting in a substantial cloud cost reduction. SeQuiLa is a cloud-native solution that can be easily deployed using any managed Kubernetes and Hadoop services available in public clouds, like Microsoft Azure Cloud, Google Cloud Platform, or Amazon Web Services. Together with the already implemented distributed range join and coverage calculations, our package provides end-users with a unified SQL interface for convenient analyses of population-scale genomic data in an interactive way. AVAILABILITY AND IMPLEMENTATION: https://biodatageeks.github.io/sequila/.

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.

High-level gonosomal mosaicism for a pathogenic non-coding CNV deletion of the lung-specific FOXF1 enhancer in an unaffected mother of an infant with ACDMPV.

BACKGROUND: Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) results from haploinsufficiency of the mesenchymal transcription factor FOXF1 gene. To date, only one case of an ACDMPV-causative CNV deletion inherited from a very-low level somatic mosaic mother has been reported. METHODS: Clinical, histopathological, and molecular studies, including whole genome sequencing, chromosomal microarray analysis, qPCR, and Sanger sequencing, followed by in vitro fertilization (IVF) with preimplantation genetic testing (PGT) were used to study a family with a deceased neonate with ACDMPV. RESULTS: A pathogenic CNV deletion of the lung-specific FOXF1 enhancer in the proband was found to be inherited from an unaffected mother, 36% mosaic for this deletion in her peripheral blood cells. The qPCR analyses of saliva, buccal cells, urine, nail, and hair samples revealed 19%, 18%, 15%, 19%, and 27% variant allele fraction, respectively, indicating a high recurrence risk. Grandparental studies revealed that the deletion arose on the mother's paternal chromosome 16. PGT studies revealed 44% embryos with the deletion, reflecting high-level germline mosaicism. CONCLUSION: Our data further demonstrate the importance of parental testing in ACDMPV families and reproductive usefulness of IVF with PGT in families with high-level parental gonosomal mosaicism.

Revised time estimation of the ancestral human chromosome 2 fusion.

BACKGROUND: The reduction of the chromosome number from 48 in the Great Apes to 46 in modern humans is thought to result from the end-to-end fusion of two ancestral non-human primate chromosomes forming the human chromosome 2 (HSA2). Genomic signatures of this event are the presence of inverted telomeric repeats at the HSA2 fusion site and a block of degenerate satellite sequences that mark the remnants of the ancestral centromere. It has been estimated that this fusion arose up to 4.5 million years ago (Mya). RESULTS: We have developed an enhanced algorithm for the detection and efficient counting of the locally over-represented weak-to-strong (AT to GC) substitutions. By analyzing the enrichment of these substitutions around the fusion site of HSA2 we estimated its formation time at 0.9 Mya with a 95% confidence interval of 0.4-1.5 Mya. Additionally, based on the statistics derived from our algorithm, we have reconstructed the evolutionary distances among the Great Apes (Hominoidea). CONCLUSIONS: Our results shed light on the HSA2 fusion formation and provide a novel computational alternative for the estimation of the speciation chronology.

Variable expressivity in a four-generation ACDMPV family with a non-coding hypermorphic SNV in trans to the frameshifting FOXF1 variant.

Heterozygous single nucleotide variants (SNVs) or copy-number variant deletions involving FOXF1 or its distant lung-specific enhancer on chromosome 16q24.1 have been identified in 80-90% of patients with Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV), a lethal neonatal lung developmental disorder. We describe a four-generation family with a deceased ACDMPV neonate, her sibling from the electively terminated pregnancy, healthy mother with a history of pulmonary arterial hypertension (PAH), an unaffected aunt, an aunt deceased due to findings consistent with ACDMPV, and a reportedly unaffected grandmother, all with the frameshifting variant c.881_902dup (p.Gly302Profs*46) in FOXF1, and a deceased great-grandmother with a history of PAH. Genome sequencing analyses in the proband's unaffected mother revealed a non-coding putative regulatory SNV rs560517434-A within the lung-specific distant FOXF1 enhancer in trans to the FOXF1 frameshift mutation. Functional testing of this variant using an in vitro luciferase reporter assay showed that it increased FOXF1 promoter activity 10-fold. Our studies further demonstrate that non-coding SNVs in the FOXF1 enhancer region can rescue the lethal ACDMPV phenotype and support the compound inheritance gene dosage model.

Do paternal deletions involving the FOXF1 locus on chromosome 16q24.1 manifest with more severe non-lung anomalies?

Alveolar capillary dysplasia with misalignment of pulmonary veins (ACDMPV) is a rare lethal lung developmental disorder in neonates due to heterozygous loss-of-function of the mesenchymal transcription factor gene, FOXF1. Interestingly, unlike ACDMPV-causing point mutations in FOXF1 that can be inherited from the mother or father, causative copy-number variant (CNV) deletions arise de novo and almost exclusively on chromosome 16 inherited from the mother (n = 50 vs. n = 3). Here, we describe a fourth case of a de novo paternal CNV deletion encompassing FOXF1, its neighboring long non-coding RNA gene FENDRR, and their distant lung-specific enhancer, identified in a 21-week-old fetus with tetralogy of Fallot, gastrointestinal and genitourinary abnormalities, a single umbilical artery, and patchy histopathological findings of ACDMPV in autopsy lung. We also review the ACDMPV-causative CNV deletions detected prenatally and propose that the majority of paternal deletions manifest with more severe additional non-lung abnormalities.

Parental mosaicism for apparent de novo genetic variants: Scope, detection, and counseling challenges.

The disease burden of de novo mutations (DNMs) has been evidenced only recently when the common application of next-generation sequencing technologies enabled their reliable and affordable detection through family-based clinical exome or genome sequencing. Implementation of exome sequencing into prenatal diagnostics revealed that up to 63% of pathogenic or likely pathogenic variants associated with fetal structural anomalies are apparently de novo, primarily for autosomal dominant disorders. Apparent DNMs have been considered to primarily occur as germline or zygotic events, with consequently negligible recurrence risks. However, there is now evidence that a considerable proportion of them are in fact inherited from a parent mosaic for the variant. Here, we review the burden of DNMs in prenatal diagnostics and the influence of parental mosaicism on the interpretation of apparent DNMs and discuss the challenges with detecting and quantifying parental mosaicism and its effect on recurrence risk. We also describe new bioinformatic and technological tools developed to assess mosaicism and discuss how they improve the accuracy of reproductive risk counseling when parental mosaicism is detected.