CHD4 and SMYD1 repress common transcriptional programs in the developing heart.

Regulation of chromatin states is essential for proper temporal and spatial gene expression. Chromatin states are modulated by remodeling complexes composed of components that have enzymatic activities. CHD4 is the catalytic core of the nucleosome remodeling and deacetylase (NuRD) complex, which represses gene transcription. However, it remains to be determined how CHD4, a ubiquitous enzyme that remodels chromatin structure, functions in cardiomyocytes to maintain heart development. In particular, whether other proteins besides the NuRD components interact with CHD4 in the heart is controversial. Using quantitative proteomics, we identified that CHD4 interacts with SMYD1, a striated muscle-restricted histone methyltransferase that is essential for cardiomyocyte differentiation and cardiac morphogenesis. Comprehensive transcriptomic and chromatin accessibility studies of Smyd1 and Chd4 null embryonic mouse hearts revealed that SMYD1 and CHD4 repress a group of common genes and pathways involved in glycolysis, response to hypoxia, and angiogenesis. Our study reveals a mechanism by which CHD4 functions during heart development, and a previously uncharacterized mechanism regarding how SMYD1 represses cardiac transcription in the developing heart.

Molecular and Spatial Signatures of Mouse Embryonic Endothelial Cells at Single-Cell Resolution.

BACKGROUND: Mature endothelial cells (ECs) are heterogeneous, with subtypes defined by tissue origin and position within the vascular bed (ie, artery, capillary, vein, and lymphatic). How this heterogeneity is established during the development of the vascular system, especially arteriovenous specification of ECs, remains incompletely characterized. METHODS: We used droplet-based single-cell RNA sequencing and multiplexed error-robust fluorescence in situ hybridization to define EC and EC progenitor subtypes from E9.5, E12.5, and E15.5 mouse embryos. We used trajectory inference to analyze the specification of arterial ECs (aECs) and venous ECs (vECs) from EC progenitors. Network analysis identified candidate transcriptional regulators of arteriovenous differentiation, which we tested by CRISPR (clustered regularly interspaced short palindromic repeats) loss of function in human-induced pluripotent stem cells undergoing directed differentiation to aECs or vECs (human-induced pluripotent stem cell-aECs or human-induced pluripotent stem cell-vECs). RESULTS: From the single-cell transcriptomes of 7682 E9.5 to E15.5 ECs, we identified 19 EC subtypes, including Etv2+Bnip3+ EC progenitors. Spatial transcriptomic analysis of 15 448 ECs provided orthogonal validation of these EC subtypes and established their spatial distribution. Most embryonic ECs were grouped by their vascular-bed types, while ECs from the brain, heart, liver, and lung were grouped by their tissue origins. Arterial (Eln, Dkk2, Vegfc, and Egfl8), venous (Fam174b and Clec14a), and capillary (Kcne3) marker genes were identified. Compared with aECs, embryonic vECs and capillary ECs shared fewer markers than their adult counterparts. Early capillary ECs with venous characteristics functioned as a branch point for differentiation of aEC and vEC lineages. CONCLUSIONS: Our results provide a spatiotemporal map of embryonic EC heterogeneity at single-cell resolution and demonstrate that the diversity of ECs in the embryo arises from both tissue origin and vascular-bed position. Developing aECs and vECs share common venous-featured capillary precursors and are regulated by distinct transcriptional regulatory networks.

Identifying novel data-driven subgroups in congenital heart disease using multi-modal measures of brain structure.

Individuals with congenital heart disease (CHD) have an increased risk of neurodevelopmental impairments. Given the hypothesized complexity linking genomics, atypical brain structure, cardiac diagnoses and their management, and neurodevelopmental outcomes, unsupervised methods may provide unique insight into neurodevelopmental variability in CHD. Using data from the Pediatric Cardiac Genomics Consortium Brain and Genes study, we identified data-driven subgroups of individuals with CHD from measures of brain structure. Using structural magnetic resonance imaging (MRI; N = 93; cortical thickness, cortical volume, and subcortical volume), we identified subgroups that differed primarily on cardiac anatomic lesion and language ability. In contrast, using diffusion MRI (N = 88; white matter connectivity strength), we identified subgroups that were characterized by differences in associations with rare genetic variants and visual-motor function. This work provides insight into the differential impacts of cardiac lesions and genomic variation on brain growth and architecture in patients with CHD, with potentially distinct effects on neurodevelopmental outcomes.

Automated Identification of Germline de novo Mutations in Family Trios: A Consensus-Based Informatic Approach.

Accurate identification of germline de novo variants (DNVs) remains a challenging problem despite rapid advances in sequencing technologies as well as methods for the analysis of the data they generate, with putative solutions often involving ad hoc filters and visual inspection of identified variants. Here, we present a purely informatic method for the identification of DNVs by analyzing short-read genome sequencing data from proband-parent trios. Our method evaluates variant calls generated by three genome sequence analysis pipelines utilizing different algorithms-GATK HaplotypeCaller, DeepTrio and Velsera GRAF-exploring the assumption that a requirement of consensus can serve as an effective filter for high-quality DNVs. We assessed the efficacy of our method by testing DNVs identified using a previously established, highly accurate classification procedure that partially relied on manual inspection and used Sanger sequencing to validate a DNV subset comprising less confident calls. The results show that our method is highly precise and that applying a force-calling procedure to putative variants further removes false-positive calls, increasing precision of the workflow to 99.6%. Our method also identified novel DNVs, 87% of which were validated, indicating it offers a higher recall rate without compromising accuracy. We have implemented this method as an automated bioinformatics workflow suitable for large-scale analyses without need for manual intervention.

Direct reprogramming of non-limb fibroblasts to cells with properties of limb progenitors.

The early limb bud consists of mesenchymal limb progenitors derived from the lateral plate mesoderm (LPM). The LPM also gives rise to the mesodermal components of the flank and neck. However, the cells at these other levels cannot produce the variety of cell types found in the limb. Taking advantage of a direct reprogramming approach, we find a set of factors (Prdm16, Zbtb16, and Lin28a) normally expressed in the early limb bud and capable of imparting limb progenitor-like properties to mouse non-limb fibroblasts. The reprogrammed cells show similar gene expression profiles and can differentiate into similar cell types as endogenous limb progenitors. The further addition of Lin41 potentiates the proliferation of the reprogrammed cells. These results suggest that these same four factors may play pivotal roles in the specification of endogenous limb progenitors.

Functional dissection of human cardiac enhancers and noncoding de novo variants in congenital heart disease.

Rare coding mutations cause ∼45% of congenital heart disease (CHD). Noncoding mutations that perturb cis-regulatory elements (CREs) likely contribute to the remaining cases, but their identification has been problematic. Using a lentiviral massively parallel reporter assay (lentiMPRA) in human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs), we functionally evaluated 6,590 noncoding de novo variants (ncDNVs) prioritized from the whole-genome sequencing of 750 CHD trios. A total of 403 ncDNVs substantially affected cardiac CRE activity. A majority increased enhancer activity, often at regions with undetectable reference sequence activity. Of ten DNVs tested by introduction into their native genomic context, four altered the expression of neighboring genes and iPSC-CM transcriptional state. To prioritize future DNVs for functional testing, we used the MPRA data to develop a regression model, EpiCard. Analysis of an independent CHD cohort by EpiCard found enrichment of DNVs. Together, we developed a scalable system to measure the effect of ncDNVs on CRE activity and deployed it to systematically assess the contribution of ncDNVs to CHD.

The BMAL1/HIF2A heterodimer modulates circadian variations of myocardial injury.

Acute myocardial infarction stands as a prominent cause of morbidity and mortality worldwide1-6. Clinical studies have demonstrated that the severity of cardiac injury following myocardial infarction exhibits a circadian pattern, with larger infarct sizes and poorer outcomes in patients experiencing morning onset myocardial infarctions7-14. However, the molecular mechanisms that govern circadian variations of myocardial injury remain unclear. Here, we show that BMAL114-20, a core circadian transcription factor, orchestrates diurnal variability in myocardial injury. Unexpectedly, BMAL1 modulates circadian-dependent cardiac injury by forming a transcriptionally active heterodimer with a non-canonical partner, hypoxia-inducible factor 2 alpha (HIF2A)6,21-23, in a diurnal manner. Substantiating this finding, we determined the cryo-EM structure of the BMAL1/HIF2A/DNA complex, revealing a previously unknown capacity for structural rearrangement within BMAL1, which enables the crosstalk between circadian rhythms and hypoxia signaling. Furthermore, we identified amphiregulin (AREG) as a rhythmic transcriptional target of the BMAL1/HIF2A heterodimer, critical for regulating circadian variations of myocardial injury. Finally, pharmacologically targeting the BMAL1/HIF2A-AREG pathway provides effective cardioprotection, with maximum efficacy when aligned with the pathway's circadian trough. Our findings not only uncover a novel mechanism governing the circadian variations of myocardial injury but also pave the way for innovative circadian-based treatment strategies, potentially shifting current treatment paradigms for myocardial infarction.

An engineered human cardiac tissue model reveals contributions of systemic lupus erythematosus autoantibodies to myocardial injury.

Systemic lupus erythematosus (SLE) is a highly heterogenous autoimmune disease that affects multiple organs, including the heart. The mechanisms by which myocardial injury develops in SLE, however, remain poorly understood. Here we engineered human cardiac tissues and cultured them with IgG fractions containing autoantibodies from SLE patients with and without myocardial involvement. We observed unique binding patterns of IgG from two patient subgroups: (i) patients with severe myocardial inflammation exhibited enhanced binding to apoptotic cells within cardiac tissues subjected to stress, and (ii) patients with systolic dysfunction exhibited enhanced binding to the surfaces of viable cardiomyocytes. Functional assays and RNA sequencing (RNA-seq) revealed that IgGs from patients with systolic dysfunction exerted direct effects on engineered tissues in the absence of immune cells, altering tissue cellular composition, respiration and calcium handling. Autoantibody target characterization by phage immunoprecipitation sequencing (PhIP-seq) confirmed distinctive IgG profiles between patient subgroups. By coupling IgG profiling with cell surface protein analyses, we identified four pathogenic autoantibody candidates that may directly alter the function of cells within the myocardium. Taken together, these observations provide insights into the cellular processes of myocardial injury in SLE that have the potential to improve patient risk stratification and inform the development of novel therapeutic strategies.

Increased endothelial sclerostin caused by elevated DSCAM mediates multiple trisomy 21 phenotypes.

Trisomy 21 (T21), a recurrent aneuploidy occurring in 1:800 births, predisposes to congenital heart disease (CHD) and multiple extracardiac phenotypes. Despite a definitive genetic etiology, the mechanisms by which T21 perturbs development and homeostasis remain poorly understood. We compared the transcriptome of CHD tissues from 49 patients with T21 and 226 with euploid CHD (eCHD). We resolved cell lineages that misexpressed T21 transcripts by cardiac single-nucleus RNA sequencing and RNA in situ hybridization. Compared with eCHD samples, T21 samples had increased chr21 gene expression; 11-fold-greater levels (P = 1.2 × 10-8) of SOST (chr17), encoding the Wnt inhibitor sclerostin; and 1.4-fold-higher levels (P = 8.7 × 10-8) of the SOST transcriptional activator ZNF467 (chr7). Euploid and T21 cardiac endothelial cells coexpressed SOST and ZNF467; however, T21 endothelial cells expressed 6.9-fold more SOST than euploid endothelial cells (P = 2.7 × 10-27). Wnt pathway genes were downregulated in T21 endothelial cells. Expression of DSCAM, residing within the chr21 CHD critical region, correlated with SOST (P = 1.9 × 10-5) and ZNF467 (P = 2.9 × 10-4). Deletion of DSCAM from T21 endothelial cells derived from human induced pluripotent stem cells diminished sclerostin secretion. As Wnt signaling is critical for atrioventricular canal formation, bone health, and pulmonary vascular homeostasis, we concluded that T21-mediated increased sclerostin levels would inappropriately inhibit Wnt activities and promote Down syndrome phenotypes. These findings imply therapeutic potential for anti-sclerostin antibodies in T21.

Genome Sequencing is Critical for Forecasting Outcomes following Congenital Cardiac Surgery.

While genome sequencing has transformed medicine by elucidating the genetic underpinnings of both rare and common complex disorders, its utility to predict clinical outcomes remains understudied. Here, we used artificial intelligence (AI) technologies to explore the predictive value of genome sequencing in forecasting clinical outcomes following surgery for congenital heart defects (CHD). We report results for a cohort of 2,253 CHD patients from the Pediatric Cardiac Genomics Consortium with a broad range of complex heart defects, pre- and post-operative clinical variables and exome sequencing. Damaging genotypes in chromatin-modifying and cilia-related genes were associated with an elevated risk of adverse post-operative outcomes, including mortality, cardiac arrest and prolonged mechanical ventilation. The impact of damaging genotypes was further amplified in the context of specific CHD phenotypes, surgical complexity and extra-cardiac anomalies. The absence of a damaging genotype in chromatin-modifying and cilia-related genes was also informative, reducing the risk for adverse postoperative outcomes. Thus, genome sequencing enriches the ability to forecast outcomes following congenital cardiac surgery.

Tbx5 maintains atrial identity in post-natal cardiomyocytes by regulating an atrial-specific enhancer network.

Understanding how the atrial and ventricular heart chambers maintain distinct identities is a prerequisite for treating chamber-specific diseases. Here, we selectively knocked out (KO) the transcription factor Tbx5 in the atrial working myocardium to evaluate its requirement for atrial identity. Atrial Tbx5 inactivation downregulated atrial cardiomyocyte (aCM) selective gene expression. Using concurrent single nucleus transcriptome and open chromatin profiling, genomic accessibility differences were identified between control and Tbx5 KO aCMs, revealing that 69% of the control-enriched ATAC regions were bound by TBX5. Genes associated with these regions were downregulated in KO aCMs, suggesting they function as TBX5-dependent enhancers. Comparing enhancer chromatin looping using H3K27ac HiChIP identified 510 chromatin loops sensitive to TBX5 dosage, and 74.8% of control-enriched loops contained anchors in control-enriched ATAC regions. Together, these data demonstrate TBX5 maintains the atrial gene expression program by binding to and preserving the tissue-specific chromatin architecture of atrial enhancers.