Integrative single-cell analysis of cardiogenesis identifies developmental trajectories and non-coding mutations in congenital heart disease.

To define the multi-cellular epigenomic and transcriptional landscape of cardiac cellular development, we generated single-cell chromatin accessibility maps of human fetal heart tissues. We identified eight major differentiation trajectories involving primary cardiac cell types, each associated with dynamic transcription factor (TF) activity signatures. We contrasted regulatory landscapes of iPSC-derived cardiac cell types and their in vivo counterparts, which enabled optimization of in vitro differentiation of epicardial cells. Further, we interpreted sequence based deep learning models of cell-type-resolved chromatin accessibility profiles to decipher underlying TF motif lexicons. De novo mutations predicted to affect chromatin accessibility in arterial endothelium were enriched in congenital heart disease (CHD) cases vs. controls. In vitro studies in iPSCs validated the functional impact of identified variation on the predicted developmental cell types. This work thus defines the cell-type-resolved cis-regulatory sequence determinants of heart development and identifies disruption of cell type-specific regulatory elements in CHD.

Transcription factor protein interactomes reveal genetic determinants in heart disease.

Congenital heart disease (CHD) is present in 1% of live births, yet identification of causal mutations remains challenging. We hypothesized that genetic determinants for CHDs may lie in the protein interactomes of transcription factors whose mutations cause CHDs. Defining the interactomes of two transcription factors haplo-insufficient in CHD, GATA4 and TBX5, within human cardiac progenitors, and integrating the results with nearly 9,000 exomes from proband-parent trios revealed an enrichment of de novo missense variants associated with CHD within the interactomes. Scoring variants of interactome members based on residue, gene, and proband features identified likely CHD-causing genes, including the epigenetic reader GLYR1. GLYR1 and GATA4 widely co-occupied and co-activated cardiac developmental genes, and the identified GLYR1 missense variant disrupted interaction with GATA4, impairing in vitro and in vivo function in mice. This integrative proteomic and genetic approach provides a framework for prioritizing and interrogating genetic variants in heart disease.

Radical S-Adenosylmethionine Enzymes in Human Health and Disease.

Radical S-adenosylmethionine (SAM) enzymes catalyze an astonishing array of complex and chemically challenging reactions across all domains of life. Of approximately 114,000 of these enzymes, 8 are known to be present in humans: MOCS1, molybdenum cofactor biosynthesis; LIAS, lipoic acid biosynthesis; CDK5RAP1, 2-methylthio-N(6)-isopentenyladenosine biosynthesis; CDKAL1, methylthio-N(6)-threonylcarbamoyladenosine biosynthesis; TYW1, wybutosine biosynthesis; ELP3, 5-methoxycarbonylmethyl uridine; and RSAD1 and viperin, both of unknown function. Aberrations in the genes encoding these proteins result in a variety of diseases. In this review, we summarize the biochemical characterization of these 8 radical S-adenosylmethionine enzymes and, in the context of human health, describe the deleterious effects that result from such genetic mutations.

Disease Model of GATA4 Mutation Reveals Transcription Factor Cooperativity in Human Cardiogenesis.

Mutation of highly conserved residues in transcription factors may affect protein-protein or protein-DNA interactions, leading to gene network dysregulation and human disease. Human mutations in GATA4, a cardiogenic transcription factor, cause cardiac septal defects and cardiomyopathy. Here, iPS-derived cardiomyocytes from subjects with a heterozygous GATA4-G296S missense mutation showed impaired contractility, calcium handling, and metabolic activity. In human cardiomyocytes, GATA4 broadly co-occupied cardiac enhancers with TBX5, another transcription factor that causes septal defects when mutated. The GATA4-G296S mutation disrupted TBX5 recruitment, particularly to cardiac super-enhancers, concomitant with dysregulation of genes related to the phenotypic abnormalities, including cardiac septation. Conversely, the GATA4-G296S mutation led to failure of GATA4 and TBX5-mediated repression at non-cardiac genes and enhanced open chromatin states at endothelial/endocardial promoters. These results reveal how disease-causing missense mutations can disrupt transcriptional cooperativity, leading to aberrant chromatin states and cellular dysfunction, including those related to morphogenetic defects.

Modeling congenital heart disease: lessons from mice, hPSC-based models, and organoids.

Congenital heart defects (CHDs) are among the most common birth defects, but their etiology has long been mysterious. In recent decades, the development of a variety of experimental models has led to a greater understanding of the molecular basis of CHDs. In this review, we contrast mouse models of CHD, which maintain the anatomical arrangement of the heart, and human cellular models of CHD, which are more likely to capture human-specific biology but lack anatomical structure. We also discuss the recent development of cardiac organoids, which are a promising step toward more anatomically informative human models of CHD.

Integrated multi-omic characterization of congenital heart disease.

The heart, the first organ to develop in the embryo, undergoes complex morphogenesis that when defective results in congenital heart disease (CHD). With current therapies, more than 90% of patients with CHD survive into adulthood, but many suffer premature death from heart failure and non-cardiac causes1. Here, to gain insight into this disease progression, we performed single-nucleus RNA sequencing on 157,273 nuclei from control hearts and hearts from patients with CHD, including those with hypoplastic left heart syndrome (HLHS) and tetralogy of Fallot, two common forms of cyanotic CHD lesions, as well as dilated and hypertrophic cardiomyopathies. We observed CHD-specific cell states in cardiomyocytes, which showed evidence of insulin resistance and increased expression of genes associated with FOXO signalling and CRIM1. Cardiac fibroblasts in HLHS were enriched in a low-Hippo and high-YAP cell state characteristic of activated cardiac fibroblasts. Imaging mass cytometry uncovered a spatially resolved perivascular microenvironment consistent with an immunodeficient state in CHD. Peripheral immune cell profiling suggested deficient monocytic immunity in CHD, in agreement with the predilection in CHD to infection and cancer2. Our comprehensive phenotyping of CHD provides a roadmap towards future personalized treatments for CHD.

The Role of Cilia and the Complex Genetics of Congenital Heart Disease.

Congenital heart disease (CHD) can affect up to 1% of live births, and despite abundant evidence of a genetic etiology, the genetic landscape of CHD is still not well understood. A large-scale mouse chemical mutagenesis screen for mutations causing CHD yielded a preponderance of cilia-related genes, pointing to a central role for cilia in CHD pathogenesis. The genes uncovered by the screen included genes that regulate ciliogenesis and cilia-transduced cell signaling as well as many that mediate endocytic trafficking, a cell process critical for both ciliogenesis and cell signaling. The clinical relevance of these findings is supported by whole-exome sequencing analysis of CHD patients that showed enrichment for pathogenic variants in ciliome genes. Surprisingly, among the ciliome CHD genes recovered were many that encoded direct protein-protein interactors. Assembly of the CHD genes into a protein-protein interaction network yielded a tight interactome that suggested this protein-protein interaction may have functional importance and that its disruption could contribute to the pathogenesis of CHD. In light of these and other findings, we propose that an interactome enriched for ciliome genes may provide the genomic context for the complex genetics of CHD and its often-observed incomplete penetrance and variable expressivity.

STIGMA: Single-cell tissue-specific gene prioritization using machine learning.

Clinical exome and genome sequencing have revolutionized the understanding of human disease genetics. Yet many genes remain functionally uncharacterized, complicating the establishment of causal disease links for genetic variants. While several scoring methods have been devised to prioritize these candidate genes, these methods fall short of capturing the expression heterogeneity across cell subpopulations within tissues. Here, we introduce single-cell tissue-specific gene prioritization using machine learning (STIGMA), an approach that leverages single-cell RNA-seq (scRNA-seq) data to prioritize candidate genes associated with rare congenital diseases. STIGMA prioritizes genes by learning the temporal dynamics of gene expression across cell types during healthy organogenesis. To assess the efficacy of our framework, we applied STIGMA to mouse limb and human fetal heart scRNA-seq datasets. In a cohort of individuals with congenital limb malformation, STIGMA prioritized 469 variants in 345 genes, with UBA2 as a notable example. For congenital heart defects, we detected 34 genes harboring nonsynonymous de novo variants (nsDNVs) in two or more individuals from a set of 7,958 individuals, including the ortholog of Prdm1, which is associated with hypoplastic left ventricle and hypoplastic aortic arch. Overall, our findings demonstrate that STIGMA effectively prioritizes tissue-specific candidate genes by utilizing single-cell transcriptome data. The ability to capture the heterogeneity of gene expression across cell populations makes STIGMA a powerful tool for the discovery of disease-associated genes and facilitates the identification of causal variants underlying human genetic disorders.

Comprehensive EHMT1 variants analysis broadens genotype-phenotype associations and molecular mechanisms in Kleefstra syndrome.

The shift to a genotype-first approach in genetic diagnostics has revolutionized our understanding of neurodevelopmental disorders, expanding both their molecular and phenotypic spectra. Kleefstra syndrome (KLEFS1) is caused by EHMT1 haploinsufficiency and exhibits broad clinical manifestations. EHMT1 encodes euchromatic histone methyltransferase-1-a pivotal component of the epigenetic machinery. We have recruited 209 individuals with a rare EHMT1 variant and performed comprehensive molecular in silico and in vitro testing alongside DNA methylation (DNAm) signature analysis for the identified variants. We (re)classified the variants as likely pathogenic/pathogenic (molecularly confirming Kleefstra syndrome) in 191 individuals. We provide an updated and broader clinical and molecular spectrum of Kleefstra syndrome, including individuals with normal intelligence and familial occurrence. Analysis of the EHMT1 variants reveals a broad range of molecular effects and their associated phenotypes, including distinct genotype-phenotype associations. Notably, we showed that disruption of the "reader" function of the ankyrin repeat domain by a protein altering variant (PAV) results in a KLEFS1-specific DNAm signature and milder phenotype, while disruption of only "writer" methyltransferase activity of the SET domain does not result in KLEFS1 DNAm signature or typical KLEFS1 phenotype. Similarly, N-terminal truncating variants result in a mild phenotype without the DNAm signature. We demonstrate how comprehensive variant analysis can provide insights into pathogenesis of the disorder and DNAm signature. In summary, this study presents a comprehensive overview of KLEFS1 and EHMT1, revealing its broader spectrum and deepening our understanding of its molecular mechanisms, thereby informing accurate variant interpretation, counseling, and clinical management.

Pathogenic variants in KMT2C result in a neurodevelopmental disorder distinct from Kleefstra and Kabuki syndromes.

Trithorax-related H3K4 methyltransferases, KMT2C and KMT2D, are critical epigenetic modifiers. Haploinsufficiency of KMT2C was only recently recognized as a cause of neurodevelopmental disorder (NDD), so the clinical and molecular spectrums of the KMT2C-related NDD (now designated as Kleefstra syndrome 2) are largely unknown. We ascertained 98 individuals with rare KMT2C variants, including 75 with protein-truncating variants (PTVs). Notably, ∼15% of KMT2C PTVs were inherited. Although the most highly expressed KMT2C transcript consists of only the last four exons, pathogenic PTVs were found in almost all the exons of this large gene. KMT2C variant interpretation can be challenging due to segmental duplications and clonal hematopoesis-induced artifacts. Using samples from 27 affected individuals, divided into discovery and validation cohorts, we generated a moderate strength disorder-specific KMT2C DNA methylation (DNAm) signature and demonstrate its utility in classifying non-truncating variants. Based on 81 individuals with pathogenic/likely pathogenic variants, we demonstrate that the KMT2C-related NDD is characterized by developmental delay, intellectual disability, behavioral and psychiatric problems, hypotonia, seizures, short stature, and other comorbidities. The facial module of PhenoScore, applied to photographs of 34 affected individuals, reveals that the KMT2C-related facial gestalt is significantly different from the general NDD population. Finally, using PhenoScore and DNAm signatures, we demonstrate that the KMT2C-related NDD is clinically and epigenetically distinct from Kleefstra and Kabuki syndromes. Overall, we define the clinical features, molecular spectrum, and DNAm signature of the KMT2C-related NDD and demonstrate they are distinct from Kleefstra and Kabuki syndromes highlighting the need to rename this condition.

Etiology of craniofacial and cardiac malformations in a mouse model of SF3B4-related syndromes.

Pathogenic variants in SF3B4, a component of the U2 snRNP complex important for branchpoint sequence recognition and splicing, are responsible for the acrofacial disorders Nager and Rodriguez Syndrome, also known as SF3B4-related syndromes. Patients exhibit malformations in the head, face, limbs, vertebrae as well as the heart. To uncover the etiology of craniofacial malformations found in SF3B4-related syndromes, mutant mouse lines with homozygous deletion of Sf3b4 in neural crest cells (NCC) were generated. Like in human patients, these embryos had craniofacial and cardiac malformations with variable expressivity and penetrance. The severity and survival of Sf3b4 NCC mutants was modified by the level of Sf3b4 in neighboring non-NCC. RNA sequencing analysis of heads of embryos prior to morphological abnormalities revealed significant changes in expression of genes forming the NCC regulatory network, as well as an increase in exon skipping. Additionally, several key histone modifiers involved in craniofacial and cardiac development showed increased exon skipping. Increased exon skipping was also associated with use of a more proximal branch point, as well as an enrichment in thymidine bases in the 50 bp around the branch points. We propose that decrease in Sf3b4 causes changes in the expression and splicing of transcripts required for proper craniofacial and cardiac development, leading to abnormalities.

Functional analysis of germline VANGL2 variants using rescue assays of vangl2 knockout zebrafish.

Developmental studies have shown that the evolutionarily conserved Wnt Planar Cell Polarity (PCP) pathway is essential for the development of a diverse range of tissues and organs including the brain, spinal cord, heart and sensory organs, as well as establishment of the left-right body axis. Germline mutations in the highly conserved PCP gene VANGL2 in humans have only been associated with central nervous system malformations, and functional testing to understand variant impact has not been performed. Here we report three new families with missense variants in VANGL2 associated with heterotaxy and congenital heart disease p.(Arg169His), non-syndromic hearing loss p.(Glu465Ala) and congenital heart disease with brain defects p.(Arg135Trp). To test the in vivo impact of these and previously described variants, we have established clinically-relevant assays using mRNA rescue of the vangl2 mutant zebrafish. We show that all variants disrupt Vangl2 function, although to different extents and depending on the developmental process. We also begin to identify that different VANGL2 missense variants may be haploinsufficient and discuss evidence in support of pathogenicity. Together, this study demonstrates that zebrafish present a suitable pipeline to investigate variants of unknown significance and suggests new avenues for investigation of the different developmental contexts of VANGL2 function that are clinically meaningful.

Complex Congenital Heart Disease in the Adult.

Congenital heart disease (CHD), a heterogeneous group of structural abnormalities of the cardiovascular system, is the most frequent cause of severe birth defects. Related to improved pediatric outcomes, there are now more adults living with CHD, including complex lesions, than children. Adults with CHD are at high risk for complications related to their underlying anatomy and past surgical palliative interventions. Adults with CHD require close monitoring and proactive management strategies to improve outcomes.

The H2Bub1-deposition complex is required for human and mouse cardiogenesis.

De novo variants affecting monoubiquitylation of histone H2B (H2Bub1) are enriched in human congenital heart disease. H2Bub1 is required in stem cell differentiation, cilia function, post-natal cardiomyocyte maturation and transcriptional elongation. However, how H2Bub1 affects cardiogenesis is unknown. We show that the H2Bub1-deposition complex (RNF20-RNF40-UBE2B) is required for mouse cardiogenesis and for differentiation of human iPSCs into cardiomyocytes. Mice with cardiac-specific Rnf20 deletion are embryonic lethal and have abnormal myocardium. We then analyzed H2Bub1 marks during differentiation of human iPSCs into cardiomyocytes. H2Bub1 is erased from most genes at the transition from cardiac mesoderm to cardiac progenitor cells but is preserved on a subset of long cardiac-specific genes. When H2Bub1 is reduced in iPSC-derived cardiomyocytes, long cardiac-specific genes have fewer full-length transcripts. This correlates with H2Bub1 accumulation near the center of these genes. H2Bub1 accumulation near the center of tissue-specific genes was also observed in embryonic fibroblasts and fetal osteoblasts. In summary, we show that normal H2Bub1 distribution is required for cardiogenesis and cardiomyocyte differentiation, and suggest that H2Bub1 regulates tissue-specific gene expression by increasing the amount of full-length transcripts.

Autonomous and non-cell autonomous role of cilia in structural birth defects in mice.

Ciliopathies are associated with wide spectrum of structural birth defects (SBDs), indicating important roles for cilia in development. Here, we provide novel insights into the temporospatial requirement for cilia in SBDs arising from deficiency in Ift140, an intraflagellar transport (IFT) protein regulating ciliogenesis. Ift140-deficient mice exhibit cilia defects accompanied by wide spectrum of SBDs including macrostomia (craniofacial defects), exencephaly, body wall defects, tracheoesophageal fistula (TEF), randomized heart looping, congenital heart defects (CHDs), lung hypoplasia, renal anomalies, and polydactyly. Tamoxifen inducible CAGGCre-ER deletion of a floxed Ift140 allele between E5.5 to 9.5 revealed early requirement for Ift140 in left-right heart looping regulation, mid to late requirement for cardiac outflow septation and alignment, and late requirement for craniofacial development and body wall closure. Surprisingly, CHD were not observed with 4 Cre drivers targeting different lineages essential for heart development, but craniofacial defects and omphalocele were observed with Wnt1-Cre targeting neural crest and Tbx18-Cre targeting epicardial lineage and rostral sclerotome through which trunk neural crest cells migrate. These findings revealed cell autonomous role of cilia in cranial/trunk neural crest-mediated craniofacial and body wall closure defects, while non-cell autonomous multi-lineage interactions underlie CHD pathogenesis, revealing unexpected developmental complexity for CHD associated with ciliopathies.

Buffering Mechanism in Aortic Arch Artery Formation and Congenital Heart Disease.

BACKGROUND: The resiliency of embryonic development to genetic and environmental perturbations has been long appreciated; however, little is known about the mechanisms underlying the robustness of developmental processes. Aberrations resulting in neonatal lethality are exemplified by congenital heart disease arising from defective morphogenesis of pharyngeal arch arteries (PAAs) and their derivatives. METHODS: Mouse genetics, lineage tracing, confocal microscopy, and quantitative image analyses were used to investigate mechanisms of PAA formation and repair. RESULTS: The second heart field (SHF) gives rise to the PAA endothelium. Here, we show that the number of SHF-derived endothelial cells (ECs) is regulated by VEGFR2 (vascular endothelial growth factor receptor 2) and Tbx1. Remarkably, when the SHF-derived EC number is decreased, PAA development can be rescued by the compensatory endothelium. Blocking such compensatory response leads to embryonic demise. To determine the source of compensating ECs and mechanisms regulating their recruitment, we investigated 3-dimensional EC connectivity, EC fate, and gene expression. Our studies demonstrate that the expression of VEGFR2 by the SHF is required for the differentiation of SHF-derived cells into PAA ECs. The deletion of 1 VEGFR2 allele (VEGFR2SHF-HET) reduces SHF contribution to the PAA endothelium, while the deletion of both alleles (VEGFR2SHF-KO) abolishes it. The decrease in SHF-derived ECs in VEGFR2SHF-HET and VEGFR2SHF-KO embryos is complemented by the recruitment of ECs from the nearby veins. Compensatory ECs contribute to PAA derivatives, giving rise to the endothelium of the aortic arch and the ductus in VEGFR2SHF-KO mutants. Blocking the compensatory response in VEGFR2SHF-KO mutants results in embryonic lethality shortly after mid-gestation. The compensatory ECs are absent in Tbx1+/- embryos, a model for 22q11 deletion syndrome, leading to unpredictable arch artery morphogenesis and congenital heart disease. Tbx1 regulates the recruitment of the compensatory endothelium in an SHF-non-cell-autonomous manner. CONCLUSIONS: Our studies uncover a novel buffering mechanism underlying the resiliency of PAA development and remodeling.

Evolutionary history of MEK1 illuminates the nature of deleterious mutations.

Mutations in signal transduction pathways lead to various diseases including cancers. MEK1 kinase, encoded by the human MAP2K1 gene, is one of the central components of the MAPK pathway and more than a hundred somatic mutations in the MAP2K1 gene were identified in various tumors. Germline mutations deregulating MEK1 also lead to congenital abnormalities, such as the cardiofaciocutaneous syndrome and arteriovenous malformation. Evaluating variants associated with a disease is a challenge, and computational genomic approaches aid in this process. Establishing evolutionary history of a gene improves computational prediction of disease-causing mutations; however, the evolutionary history of MEK1 is not well understood. Here, by revealing a precise evolutionary history of MEK1, we construct a well-defined dataset of MEK1 metazoan orthologs, which provides sufficient depth to distinguish between conserved and variable amino acid positions. We matched known and predicted disease-causing and benign mutations to evolutionary changes observed in corresponding amino acid positions and found that all known and many suspected disease-causing mutations are evolutionarily intolerable. We selected several variants that cannot be unambiguously assessed by automated prediction tools but that are confidently identified as "damaging" by our approach, for experimental validation in Drosophila. In all cases, evolutionary intolerant variants caused increased mortality and severe defects in fruit fly embryos confirming their damaging nature. We anticipate that our analysis will serve as a blueprint to help evaluate known and novel missense variants in MEK1 and that our approach will contribute to improving automated tools for disease-associated variant interpretation.

Cardiac progenitors instruct second heart field fate through Wnts.

The heart develops in a synchronized sequence of proliferation and differentiation of cardiac progenitor cells (CPCs) from two anatomically distinct pools of cells, the first heart field (FHF) and second heart field (SHF). Congenital heart defects arise upon dysregulation of these processes, many of which are restricted to derivatives of the FHF or SHF. Of the conserved set of signaling pathways that regulate development, the Wnt signaling pathway has long been known for its importance in SHF development. The source of such Wnts has remained elusive, though it has been postulated that these Wnts are secreted from ectodermal or endodermal sources. The central question remains unanswered: Where do these Wnts come from? Here, we show that CPCs autoregulate SHF development via Wnt through genetic manipulation of a key Wnt export protein (Wls), scRNA-seq analysis of CPCs, and use of our precardiac organoid system. Through this, we identify dysregulated developmental trajectories of anterior SHF cell fate, leading to a striking single ventricle phenotype in knockout embryos. We then applied our findings to our precardiac organoid model and found that Wnt2 is sufficient to restore SHF cell fate in our model of disrupted endogenous Wnt signaling. In this study, we provide a basis for SHF cell fate decision-proliferation vs. differentiation-autoregulated by CPCs through Wnt.

Pleiotropic role of TRAF7 in skull-base meningiomas and congenital heart disease.

While somatic variants of TRAF7 (Tumor necrosis factor receptor-associated factor 7) underlie anterior skull-base meningiomas, here we report the inherited mutations of TRAF7 that cause congenital heart defects. We show that TRAF7 mutants operate in a dominant manner, inhibiting protein function via heterodimerization with wild-type protein. Further, the shared genetics of the two disparate pathologies can be traced to the common origin of forebrain meninges and cardiac outflow tract from the TRAF7-expressing neural crest. Somatic and inherited mutations disrupt TRAF7-IFT57 interactions leading to cilia degradation. TRAF7-mutant meningioma primary cultures lack cilia, and TRAF7 knockdown causes cardiac, craniofacial, and ciliary defects in Xenopus and zebrafish, suggesting a mechanistic convergence for TRAF7-driven meningiomas and developmental heart defects.

Mutations in Bcl9 and Pygo genes cause congenital heart defects by tissue-specific perturbation of Wnt/ß-catenin signaling.

Bcl9 and Pygopus (Pygo) are obligate Wnt/ß-catenin cofactors in Drosophila, yet their contribution to Wnt signaling during vertebrate development remains unresolved. Combining zebrafish and mouse genetics, we document a conserved, ß-catenin-associated function for BCL9 and Pygo proteins during vertebrate heart development. Disrupting the ß-catenin-BCL9-Pygo complex results in a broadly maintained canonical Wnt response yet perturbs heart development and proper expression of key cardiac regulators. Our work highlights BCL9 and Pygo as selective ß-catenin cofactors in a subset of canonical Wnt responses during vertebrate development. Moreover, our results implicate alterations in BCL9 and BCL9L in human congenital heart defects.