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.

poly(UG)-tailed RNAs in genome protection and epigenetic inheritance.

Mobile genetic elements threaten genome integrity in all organisms. RDE-3 (also known as MUT-2) is a ribonucleotidyltransferase that is required for transposon silencing and RNA interference in Caenorhabditis elegans1-4. When tethered to RNAs in heterologous expression systems, RDE-3 can add long stretches of alternating non-templated uridine (U) and guanosine (G) ribonucleotides to the 3' termini of these RNAs (designated poly(UG) or pUG tails)5. Here we show that, in its natural context in C. elegans, RDE-3 adds pUG tails to targets of RNA interference, as well as to transposon RNAs. RNA fragments attached to pUG tails with more than 16 perfectly alternating 3' U and G nucleotides become gene-silencing agents. pUG tails promote gene silencing by recruiting RNA-dependent RNA polymerases, which use pUG-tailed RNAs (pUG RNAs) as templates to synthesize small interfering RNAs (siRNAs). Our results show that cycles of pUG RNA-templated siRNA synthesis and siRNA-directed pUG RNA biogenesis underlie double-stranded-RNA-directed transgenerational epigenetic inheritance in the C. elegans germline. We speculate that this pUG RNA-siRNA silencing loop enables parents to inoculate progeny against the expression of unwanted or parasitic genetic elements.

Cells of the adult human heart.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

Neither cardiac mitochondrial DNA variation nor copy number contribute to congenital heart disease risk.

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.

An ancient founder mutation located between ROBO1 and ROBO2 is responsible for increased microtia risk in Amerindigenous populations.

Microtia is a congenital malformation that encompasses mild hypoplasia to complete loss of the external ear, or pinna. Although the contribution of genetic variation and environmental factors to microtia remains elusive, Amerindigenous populations have the highest reported incidence. Here, using both transmission disequilibrium tests and association studies in microtia trios (parents and affected child) and microtia cohorts enrolled in Latin America, we map an ∼10-kb microtia locus (odds ratio = 4.7; P = 6.78e-18) to the intergenic region between Roundabout 1 (ROBO1) and Roundabout 2 (ROBO2) (chr3: 78546526 to 78555137). While alleles at the microtia locus significantly increase the risk of microtia, their penetrance is low (<1%). We demonstrate that the microtia locus contains a polymorphic complex repeat element that is expanded in affected individuals. The locus is located near a chromatin loop region that regulates ROBO1 and ROBO2 expression in induced pluripotent stem cell­derived neural crest cells. Furthermore, we use single nuclear RNA sequencing to demonstrate ROBO1 and ROBO2 expression in both fibroblasts and chondrocytes of the mature human pinna. Because the microtia allele is enriched in Amerindigenous populations and is shared by some East Asian subjects with craniofacial malformations, we propose that both populations share a mutation that arose in a common ancestor prior to the ancient migration of Eurasian populations into the Americas and that the high incidence of microtia among Amerindigenous populations reflects the population bottleneck that occurred during the migration out of Eurasia.

Ablation of lysophosphatidic acid receptor 1 attenuates hypertrophic cardiomyopathy in a mouse model.

Myocardial fibrosis is a key pathologic feature of hypertrophic cardiomyopathy (HCM). However, the fibrotic pathways activated by HCM-causing sarcomere protein gene mutations are poorly defined. Because lysophosphatidic acid is a mediator of fibrosis in multiple organs and diseases, we tested the role of the lysophosphatidic acid pathway in HCM. Lysphosphatidic acid receptor 1 (LPAR1), a cell surface receptor, is required for lysophosphatidic acid mediation of fibrosis. We bred HCM mice carrying a pathogenic myosin heavy-chain variant (403+/-) with Lpar1-ablated mice to create mice carrying both genetic changes (403+/- LPAR1 -/-) and assessed development of cardiac hypertrophy and fibrosis. Compared with 403+/- LPAR1WT, 403+/- LPAR1 -/- mice developed significantly less hypertrophy and fibrosis. Single-nucleus RNA sequencing of left ventricular tissue demonstrated that Lpar1 was predominantly expressed by lymphatic endothelial cells (LECs) and cardiac fibroblasts. Lpar1 ablation reduced the population of LECs, confirmed by immunofluorescence staining of the LEC markers Lyve1 and Ccl21a and, by in situ hybridization, for Reln and Ccl21a. Lpar1 ablation also altered the distribution of fibroblast cell states. FB1 and FB2 fibroblasts decreased while FB0 and FB3 fibroblasts increased. Our findings indicate that Lpar1 is expressed predominantly by LECs and fibroblasts in the heart and is required for development of hypertrophy and fibrosis in an HCM mouse model. LPAR1 antagonism, including agents in clinical trials for other fibrotic diseases, may be beneficial for HCM.

Discordant clinical features of identical hypertrophic cardiomyopathy twins.

Hypertrophic cardiomyopathy (HCM) is a disease of heart muscle, which affects ∼1 in 500 individuals and is characterized by increased left ventricular wall thickness. While HCM is caused by pathogenic variants in any one of eight sarcomere protein genes, clinical expression varies considerably, even among patients with the same pathogenic variant. To determine whether background genetic variation or environmental factors drive these differences, we studied disease progression in 11 pairs of monozygotic HCM twins. The twin pairs were followed for 5 to 14 y, and left ventricular wall thickness, left atrial diameter, and left ventricular ejection fraction were collected from echocardiograms at various time points. All nine twin pairs with sarcomere protein gene variants and two with unknown disease etiologies had discordant morphologic features of the heart, demonstrating the influence of nonhereditable factors on clinical expression of HCM. Whole genome sequencing analysis of the six monozygotic twins with discordant HCM phenotypes did not reveal notable somatic genetic variants that might explain their clinical differences. Discordant cardiac morphology of identical twins highlights a significant role for epigenetics and environment in HCM disease progression.

Pathogenesis of Cardiomyopathy Caused by Variants in ALPK3, an Essential Pseudokinase in the Cardiomyocyte Nucleus and Sarcomere.

BACKGROUND: ALPK3 encodes α-kinase 3, a muscle-specific protein of unknown function. ALPK3 loss-of-function variants cause cardiomyopathy with distinctive clinical manifestations in both children and adults, but the molecular functions of ALPK3 remain poorly understood. METHODS: We explored the putative kinase activity of ALPK3 and the consequences of damaging variants using isogenic human induced pluripotent stem cell-derived cardiomyocytes, mice, and human patient tissues. RESULTS: Multiple sequence alignment of all human α-kinase domains and their orthologs revealed 4 conserved residues that were variant only in ALPK3, demonstrating evolutionary divergence of the ALPK3 α-kinase domain sequence. Phosphoproteomic evaluation of both ALPK3 kinase domain inhibition and overexpression failed to detect significant changes in catalytic activity, establishing ALPK3 as a pseudokinase. Investigations into alternative functions revealed that ALPK3 colocalized with myomesin proteins (MYOM1, MYOM2) at both the nuclear envelope and the sarcomere M-band. ALPK3 loss-of-function variants caused myomesin proteins to mislocalize and also dysregulated several additional M-band proteins involved in sarcomere protein turnover, which ultimately impaired cardiomyocyte structure and function. CONCLUSIONS: ALPK3 is an essential cardiac pseudokinase that inserts in the nuclear envelope and the sarcomere M-band. Loss of ALPK3 causes mislocalization of myomesins, critical force-buffering proteins in cardiomyocytes, and also dysregulates M-band proteins necessary for sarcomere protein turnover. We conclude that ALPK3 cardiomyopathy induces ventricular dilatation caused by insufficient myomesin-mediated force buffering and hypertrophy by impairment of sarcomere proteostasis.

Filamin C Cardiomyopathy Variants Cause Protein and Lysosome Accumulation.

[Figure: see text].

Mechanisms of Congenital Heart Disease Caused by NAA15 Haploinsufficiency.

[Figure: see text].

Genetic and Phenotypic Landscape of Peripartum Cardiomyopathy.

BACKGROUND: Peripartum cardiomyopathy (PPCM) occurs in ≈1:2000 deliveries in the United States and worldwide. The genetic underpinnings of PPCM remain poorly defined. Approximately 10% of women with PPCM harbor truncating variants in TTN (TTNtvs). Whether mutations in other genes can predispose to PPCM is not known. It is also not known if the presence of TTNtvs predicts clinical presentation or outcomes. Nor is it known if the prevalence of TTNtvs differs in women with PPCM and preeclampsia, the strongest risk factor for PPCM. METHODS: Women with PPCM were retrospectively identified from several US and international academic centers, and clinical information and DNA samples were acquired. Next-generation sequencing was performed on 67 genes, including TTN, and evaluated for burden of truncating and missense variants. The impact of TTNtvs on the severity of clinical presentation, and on clinical outcomes, was evaluated. RESULTS: Four hundred sixty-nine women met inclusion criteria. Of the women with PPCM, 10.4% bore TTNtvs (odds ratio=9.4 compared with 1.2% in the reference population; Bonferroni-corrected P [P*]=1.2×10-46). We additionally identified overrepresentation of truncating variants in FLNC (odds ratio=24.8, P*=7.0×10-8), DSP (odds ratio=14.9, P*=1.0×10-8), and BAG3 (odds ratio=53.1, P*=0.02), genes not previously associated with PPCM. This profile is highly similar to that found in nonischemic dilated cardiomyopathy. Women with TTNtvs had lower left ventricular ejection fraction on presentation than did women without TTNtvs (23.5% versus 29%, P=2.5×10-4), but did not differ significantly in timing of presentation after delivery, in prevalence of preeclampsia, or in rates of clinical recovery. CONCLUSIONS: This study provides the first extensive genetic and phenotypic landscape of PPCM and demonstrates that predisposition to heart failure is an important risk factor for PPCM. The work reveals a degree of genetic similarity between PPCM and dilated cardiomyopathy, suggesting that gene-specific therapeutic approaches being developed for dilated cardiomyopathy may also apply to PPCM, and that approaches to genetic testing in PPCM should mirror those taken in dilated cardiomyopathy. Last, the clarification of genotype/phenotype associations has important implications for genetic counseling.

Abnormal Left-Hemispheric Sulcal Patterns Correlate with Neurodevelopmental Outcomes in Subjects with Single Ventricular Congenital Heart Disease.

Neurodevelopmental abnormalities are the most common noncardiac complications in patients with congenital heart disease (CHD). Prenatal brain abnormalities may be due to reduced oxygenation, genetic factors, or less commonly, teratogens. Understanding the contribution of these factors is essential to improve outcomes. Because primary sulcal patterns are prenatally determined and under strong genetic control, we hypothesized that they are influenced by genetic variants in CHD. In this study, we reveal significant alterations in sulcal patterns among subjects with single ventricle CHD (n = 115, 14.7 ± 2.9 years [mean ± standard deviation]) compared with controls (n = 45, 15.5 ± 2.4 years) using a graph-based pattern-analysis technique. Among patients with CHD, the left hemisphere demonstrated decreased sulcal pattern similarity to controls in the left temporal and parietal lobes, as well as the bilateral frontal lobes. Temporal and parietal lobes demonstrated an abnormally asymmetric left-right pattern of sulcal basin area in CHD subjects. Sulcal pattern similarity to control was positively correlated with working memory, processing speed, and executive function. Exome analysis identified damaging de novo variants only in CHD subjects with more atypical sulcal patterns. Together, these findings suggest that sulcal pattern analysis may be useful in characterizing genetically influenced, atypical early brain development and neurodevelopmental risk in subjects with CHD.

ORE identifies extreme expression effects enriched for rare variants.

MOTIVATION: Non-coding rare variants (RVs) may contribute to Mendelian disorders but have been challenging to study due to small sample sizes, genetic heterogeneity and uncertainty about relevant non-coding features. Previous studies identified RVs associated with expression outliers, but varying outlier definitions were employed and no comprehensive open-source software was developed. RESULTS: We developed Outlier-RV Enrichment (ORE) to identify biologically-meaningful non-coding RVs. We implemented ORE combining whole-genome sequencing and cardiac RNAseq from congenital heart defect patients from the Pediatric Cardiac Genomics Consortium and deceased adults from Genotype-Tissue Expression. Use of rank-based outliers maximized sensitivity while a most extreme outlier approach maximized specificity. Rarer variants had stronger associations, suggesting they are under negative selective pressure and providing a basis for investigating their contribution to Mendelian disorders. AVAILABILITY AND IMPLEMENTATION: ORE, source code, and documentation are available at https://pypi.python.org/pypi/ore under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cardiovascular homeostasis dependence on MICU2, a regulatory subunit of the mitochondrial calcium uniporter.

Comparative analyses of transcriptional profiles from humans and mice with cardiovascular pathologies revealed consistently elevated expression of MICU2, a regulatory subunit of the mitochondrial calcium uniporter complex. To determine if MICU2 expression was cardioprotective, we produced and characterized Micu2-/- mice. Mutant mice had left atrial enlargement and Micu2-/- cardiomyocytes had delayed sarcomere relaxation and cytosolic calcium reuptake kinetics, indicating diastolic dysfunction. RNA sequencing (RNA-seq) of Micu2-/- ventricular tissues revealed markedly reduced transcripts encoding the apelin receptor (Micu2-/- vs. wild type, P = 7.8 × 10-40), which suppresses angiotensin II receptor signaling via allosteric transinhibition. We found that Micu2-/- and wild-type mice had comparable basal blood pressures and elevated responses to angiotensin II infusion, but that Micu2-/- mice exhibited systolic dysfunction and 30% lethality from abdominal aortic rupture. Aneurysms and rupture did not occur with norepinephrine-induced hypertension. Aortic tissue from Micu2-/- mice had increased expression of extracellular matrix remodeling genes, while single-cell RNA-seq analyses showed increased expression of genes related to reactive oxygen species, inflammation, and proliferation in fibroblast and smooth muscle cells. We concluded that Micu2-/- mice recapitulate features of diastolic heart disease and define previously unappreciated roles for Micu2 in regulating angiotensin II-mediated hypertensive responses that are critical in protecting the abdominal aorta from injury.

Identification of pathogenic gene mutations in LMNA and MYBPC3 that alter RNA splicing.

Genetic variants that cause haploinsufficiency account for many autosomal dominant (AD) disorders. Gene-based diagnosis classifies variants that alter canonical splice signals as pathogenic, but due to imperfect understanding of RNA splice signals other variants that may create or eliminate splice sites are often clinically classified as variants of unknown significance (VUS). To improve recognition of pathogenic splice-altering variants in AD disorders, we used computational tools to prioritize VUS and developed a cell-based minigene splicing assay to confirm aberrant splicing. Using this two-step procedure we evaluated all rare variants in two AD cardiomyopathy genes, lamin A/C (LMNA) and myosin binding protein C (MYBPC3). We demonstrate that 13 LMNA and 35 MYBPC3 variants identified in cardiomyopathy patients alter RNA splicing, representing a 50% increase in the numbers of established damaging splice variants in these genes. Over half of these variants are annotated as VUS by clinical diagnostic laboratories. Familial analyses of one variant, a synonymous LMNA VUS, demonstrated segregation with cardiomyopathy affection status and altered cardiac LMNA splicing. Application of this strategy should improve diagnostic accuracy and variant classification in other haploinsufficient AD disorders.

Genome sequencing as a first-line genetic test in familial dilated cardiomyopathy.

PURPOSE: We evaluated genome sequencing (GS) as an alternative to multigene panel sequencing (PS) for genetic testing in dilated cardiomyopathy (DCM). METHODS: Forty-two patients with familial DCM underwent PS and GS, and detection rates of rare single-nucleotide variants and small insertions/deletions in panel genes were compared. Loss-of-function variants in 406 cardiac-enriched genes were evaluated, and an assessment of structural variation was performed. RESULTS: GS provided broader and more uniform coverage than PS, with high concordance for rare variant detection in panel genes. GS identified all PS-identified pathogenic or likely pathogenic variants as well as two additional likely pathogenic variants: one was missed by PS due to low coverage, the other was a known disease-causing variant in a gene not included on the panel. No loss-of-function variants in the extended gene set met clinical criteria for pathogenicity. One BAG3 structural variant was classified as pathogenic. CONCLUSION: Our data support the use of GS for genetic testing in DCM, with high variant detection accuracy and a capacity to identify structural variants. GS provides an opportunity to go beyond suites of established disease genes, but the incremental yield of clinically actionable variants is limited by a paucity of genetic and functional evidence for DCM association.

A gene-centric strategy for identifying disease-causing rare variants in dilated cardiomyopathy.

PURPOSE: We evaluated strategies for identifying disease-causing variants in genetic testing for dilated cardiomyopathy (DCM). METHODS: Cardiomyopathy gene panel testing was performed in 532 DCM patients and 527 healthy control subjects. Rare variants in 41 genes were stratified using variant-level and gene-level characteristics. RESULTS: A majority of DCM cases and controls carried rare protein-altering cardiomyopathy gene variants. Variant-level characteristics alone had limited discriminative value. Differentiation between groups was substantially improved by addition of gene-level information that incorporated ranking of genes based on literature evidence for disease association. The odds of DCM were increased to nearly 9-fold for truncating variants or high-impact missense variants in the subset of 14 genes that had the strongest biological links to DCM (P <0.0001). For some of these genes, DCM-associated variants appeared to be clustered in key protein functional domains. Multiple rare variants were present in many family probands, however, there was generally only one "driver" pathogenic variant that cosegregated with disease. CONCLUSION: Rare variants in cardiomyopathy genes can be effectively stratified by combining variant-level and gene-level information. Prioritization of genes based on their a priori likelihood of disease causation is a key factor in identifying clinically actionable variants in cardiac genetic testing.

Robust identification of mosaic variants in congenital heart disease.

Mosaicism due to somatic mutations can cause multiple diseases including cancer, developmental and overgrowth syndromes, neurodevelopmental disorders, autoinflammatory diseases, and atrial fibrillation. With the increased use of next generation sequencing technology, multiple tools have been developed to identify low-frequency variants, specifically from matched tumor-normal tissues in cancer studies. To investigate whether mosaic variants are implicated in congenital heart disease (CHD), we developed a pipeline using the cancer somatic variant caller MuTect to identify mosaic variants in whole-exome sequencing (WES) data from a cohort of parent/affected child trios (n = 715) and a cohort of healthy individuals (n = 416). This is a novel application of the somatic variant caller designed for cancer to WES trio data. We identified two cases with mosaic KMT2D mutations that are likely pathogenic for CHD, but conclude that, overall, mosaicism detectable in peripheral blood or saliva does not account for a significant portion of CHD etiology.

Cardiac myosin binding protein C regulates postnatal myocyte cytokinesis.

Homozygous cardiac myosin binding protein C-deficient (Mybpc(t/t)) mice develop dramatic cardiac dilation shortly after birth; heart size increases almost twofold. We have investigated the mechanism of cardiac enlargement in these hearts. Throughout embryogenesis myocytes undergo cell division while maintaining the capacity to pump blood by rapidly disassembling and reforming myofibrillar components of the sarcomere throughout cell cycle progression. Shortly after birth, myocyte cell division ceases. Cardiac MYBPC is a thick filament protein that regulates sarcomere organization and rigidity. We demonstrate that many Mybpc(t/t) myocytes undergo an additional round of cell division within 10 d postbirth compared with their wild-type counterparts, leading to increased numbers of mononuclear myocytes. Short-hairpin RNA knockdown of Mybpc3 mRNA in wild-type mice similarly extended the postnatal window of myocyte proliferation. However, adult Mybpc(t/t) myocytes are unable to fully regenerate the myocardium after injury. MYBPC has unexpected inhibitory functions during postnatal myocyte cytokinesis and cell cycle progression. We suggest that human patients with homozygous MYBPC3-null mutations develop dilated cardiomyopathy, coupled with myocyte hyperplasia (increased cell number), as observed in Mybpc(t/t) mice. Human patients, with heterozygous truncating MYBPC3 mutations, like mice with similar mutations, have hypertrophic cardiomyopathy. However, the mechanism leading to hypertrophic cardiomyopathy in heterozygous MYBPC3(+/-) individuals is myocyte hypertrophy (increased cell size), whereas the mechanism leading to cardiac dilation in homozygous Mybpc3(-/-) mice is primarily myocyte hyperplasia.