Activation of complex molecular programs in specific cell lineages governs mammalian heart development, from a primordial linear tube to a four-chamber organ. To characterize lineage-specific, spatiotemporal developmental programs, we performed single-cell RNA sequencing of >1,200 murine cells isolated at seven time points spanning embryonic day 9.5 (primordial heart tube) to postnatal day 21 (mature heart). Using unbiased transcriptional data, we classified cardiomyocytes, endothelial cells, and fibroblast-enriched cells, thus identifying markers for temporal and chamber-specific developmental programs. By harnessing these datasets, we defined developmental ages of human and mouse pluripotent stem-cell-derived cardiomyocytes and characterized lineage-specific maturation defects in hearts of mice with heterozygous mutations in Nkx2.5 that cause human heart malformations. This spatiotemporal transcriptome analysis of heart development reveals lineage-specific gene programs underlying normal cardiac development and congenital heart disease.

The transcriptome is subject to multiple changes during pathogenesis, including the use of alternate 5' start-sites that can affect transcription levels and output. Current RNA sequencing techniques can assess mRNA levels, but do not robustly detect changes in 5' start-site use. Here, we developed a transcriptome sequencing strategy that detects genome-wide changes in start-site usage (5'RNA-Seq) and applied this methodology to identify regulatory events that occur in hypertrophic cardiomyopathy (HCM). Compared with transcripts from WT mice, 92 genes had altered start-site usage in a mouse model of HCM, including four-and-a-half LIM domains protein 1 (Fhl1). HCM-induced altered transcriptional regulation of Fhl1 resulted in robust myocyte expression of a distinct protein isoform, a response that was conserved in humans with genetic or acquired cardiomyopathies. Genetic ablation of Fhl1 in HCM mice was deleterious, which suggests that Fhl1 transcriptional changes provide salutary effects on stressed myocytes in this disease. Because Fhl1 is a chromosome X-encoded gene, stress-induced changes in its transcription may contribute to gender differences in the clinical severity of HCM. Our findings indicate that 5'RNA-Seq has the potential to identify genome-wide changes in 5' start-site usage that are associated with pathogenic phenotypes.

Background

The contribution of somatic mosaicism, or genetic mutations arising after oocyte fertilization, to congenital heart disease (CHD) is not well understood. Further, the relationship between mosaicism in blood and cardiovascular tissue has not been determined.

Methods

We developed a new computational method, EM-mosaic (Expectation-Maximization-based detection of mosaicism), to analyze mosaicism in exome sequences derived primarily from blood DNA of 2530 CHD proband-parent trios. To optimize this method, we measured mosaic detection power as a function of sequencing depth. In parallel, we analyzed our cohort using MosaicHunter, a Bayesian genotyping algorithm-based mosaic detection tool, and compared the two methods. The accuracy of these mosaic variant detection algorithms was assessed using an independent resequencing method. We then applied both methods to detect mosaicism in cardiac tissue-derived exome sequences of 66 participants for which matched blood and heart tissue was available.

Results

EM-mosaic detected 326 mosaic mutations in blood and/or cardiac tissue DNA. Of the 309 detected in blood DNA, 85/97 (88%) tested were independently confirmed, while 7/17 (41%) candidates of 17 detected in cardiac tissue were confirmed. MosaicHunter detected an additional 64 mosaics, of which 23/46 (50%) among 58 candidates from blood and 4/6 (67%) of 6 candidates from cardiac tissue confirmed. Twenty-five mosaic variants altered CHD-risk genes, affecting 1% of our cohort. Of these 25, 22/22 candidates tested were confirmed. Variants predicted as damaging had higher variant allele fraction than benign variants, suggesting a role in CHD. The estimated true frequency of mosaic variants above 10% mosaicism was 0.14/person in blood and 0.21/person in cardiac tissue. Analysis of 66 individuals with matched cardiac tissue available revealed both tissue-specific and shared mosaicism, with shared mosaics generally having higher allele fraction.

Conclusions

We estimate that ~ 1% of CHD probands have a mosaic variant detectable in blood that could contribute to cardiac malformations, particularly those damaging variants with relatively higher allele fraction. Although blood is a readily available DNA source, cardiac tissues analyzed contributed ~ 5% of somatic mosaic variants identified, indicating the value of tissue mosaicism analyses.

The well-established manifestation of mitochondrial mutations in functional cardiac disease (e.g., mitochondrial cardiomyopathy) prompted the hypothesis that mitochondrial DNA (mtDNA) sequence and/or copy number (mtDNAcn) variation contribute to cardiac defects in congenital heart disease (CHD). MtDNAcns were calculated and rare, non-synonymous mtDNA mutations were identified in 1,837 CHD-affected proband-parent trios, 116 CHD-affected singletons, and 114 paired cardiovascular tissue/blood samples. The variant allele fraction (VAF) of heteroplasmic variants in mitochondrial RNA from 257 CHD cardiovascular tissue samples was also calculated. On average, mtDNA from blood had 0.14 rare variants and 52.9 mtDNA copies per nuclear genome per proband. No variation with parental age at proband birth or CHD-affected proband age was seen. mtDNAcns in valve/vessel tissue (320 ± 70) were lower than in atrial tissue (1,080 ± 320, p = 6.8E-21), which were lower than in ventricle tissue (1,340 ± 280, p = 1.4E-4). The frequency of rare variants in CHD-affected individual DNA was indistinguishable from the frequency in an unaffected cohort, and proband mtDNAcns did not vary from those of CHD cohort parents. In both the CHD and the comparison cohorts, mtDNAcns were significantly correlated between mother-child, father-child, and mother-father. mtDNAcns among people with European (mean = 52.0), African (53.0), and Asian haplogroups (53.5) were calculated and were significantly different for European and Asian haplogroups (p = 2.6E-3). Variant heteroplasmic fraction (HF) in blood correlated well with paired cardiovascular tissue HF (r = 0.975) and RNA VAF (r = 0.953), which suggests blood HF is a reasonable proxy for HF in heart tissue. We conclude that mtDNA mutations and mtDNAcns are unlikely to contribute significantly to CHD risk.

Recent human genetic studies have provided evidences that sporadic or inherited missense mutations in four-and-a-half LIM domain protein 1 (FHL1), resulting in alterations in FHL1 protein expression, are associated with rare congenital myopathies, including reducing body myopathy and Emery-Dreifuss muscular dystrophy. However, it remains to be clarified whether mutations in FHL1 cause skeletal muscle remodeling owing to gain- or loss of FHL1 function. In this study, we used FHL1-null mice lacking global FHL1 expression to evaluate loss-of-function effects on skeletal muscle homeostasis. Histological and functional analyses of soleus, tibialis anterior and sternohyoideus muscles demonstrated that FHL1-null mice develop an age-dependent myopathy associated with myofibrillar and intermyofibrillar (mitochondrial and sarcoplasmic reticulum) disorganization, impaired muscle oxidative capacity and increased autophagic activity. A longitudinal study established decreased survival rates in FHL1-null mice, associated with age-dependent impairment of muscle contractile function and a significantly lower exercise capacity. Analysis of primary myoblasts isolated from FHL1-null muscles demonstrated early muscle fiber differentiation and maturation defects, which could be rescued by re-expression of the FHL1A isoform, highlighting that FHL1A is necessary for proper muscle fiber differentiation and maturation in vitro. Overall, our data show that loss of FHL1 function leads to myopathy in vivo and suggest that loss of function of FHL1 may be one of the mechanisms underlying muscle dystrophy in patients with FHL1 mutations.

Damaging GATA6 variants cause cardiac outflow tract defects, sometimes with pancreatic and diaphragmic malformations. To define molecular mechanisms for these diverse developmental defects, we studied transcriptional and epigenetic responses to GATA6 loss of function (LoF) and missense variants during cardiomyocyte differentiation of isogenic human induced pluripotent stem cells. We show that GATA6 is a pioneer factor in cardiac development, regulating SMYD1 that activates HAND2, and KDR that with HAND2 orchestrates outflow tract formation. LoF variants perturbed cardiac genes and also endoderm lineage genes that direct PDX1 expression and pancreatic development. Remarkably, an exon 4 GATA6 missense variant, highly associated with extra-cardiac malformations, caused ectopic pioneer activities, profoundly diminishing GATA4, FOXA1/2, and PDX1 expression and increasing normal retinoic acid signaling that promotes diaphragm development. These aberrant epigenetic and transcriptional signatures illuminate the molecular mechanisms for cardiovascular malformations, pancreas and diaphragm dysgenesis that arise in patients with distinct GATA6 variants.

A genetic etiology is identified for one-third of patients with congenital heart disease (CHD), with 8% of cases attributable to coding de novo variants (DNVs). To assess the contribution of noncoding DNVs to CHD, we compared genome sequences from 749 CHD probands and their parents with those from 1,611 unaffected trios. Neural network prediction of noncoding DNV transcriptional impact identified a burden of DNVs in individuals with CHD (n = 2,238 DNVs) compared to controls (n = 4,177; P = 8.7 × 10^-4). Independent analyses of enhancers showed an excess of DNVs in associated genes (27 genes versus 3.7 expected, P = 1 × 10^-5). We observed significant overlap between these transcription-based approaches (odds ratio (OR) = 2.5, 95% confidence interval (CI) 1.1-5.0, P = 5.4 × 10^-3). CHD DNVs altered transcription levels in 5 of 31 enhancers assayed. Finally, we observed a DNV burden in RNA-binding-protein regulatory sites (OR = 1.13, 95% CI 1.1-1.2, P = 8.8 × 10^-5). Our findings demonstrate an enrichment of potentially disruptive regulatory noncoding DNVs in a fraction of CHD at least as high as that observed for damaging coding DNVs.