Search

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.

Abstract Cardiomyopathies have unresolved genotype–phenotype relationships and lack disease-specific treatments. Here we provide a framework to identify genotype-specific pathomechanisms and therapeutic targets to accelerate the development of precision medicine. We use human cardiac electromechanical in-silico modelling and simulation which we validate with experimental hiPSC-CM data and modelling in combination with clinical biomarkers. We select hypertrophic cardiomyopathy as a challenge for this approach and study genetic variations that mutate proteins of the thick (MYH7 R403Q/+) and thin filaments (TNNT2 R92Q/+, TNNI3 R21C/+) of the cardiac sarcomere. Using in-silico techniques we show that the destabilisation of myosin super relaxation observed in hiPSC-CMs drives disease in virtual cells and ventricles carrying the MYH7R403Q/+ variant, and that secondary effects on thin filament activation are necessary to precipitate slowed relaxation of the cell and diastolic insufficiency in the chamber. In-silico modelling shows that Mavacamten corrects the MYH7R403Q/+ phenotype in agreement with hiPSC-CM experiments. Our in-silico model predicts that the thin filament variants TNNT2R92Q/+ and TNNI3R21C/+ display altered calcium regulation as central pathomechanism, for which Mavacamten provides incomplete salvage, which we have corroborated in TNNT2R92Q/+ and TNNI3R21C/+ hiPSC-CMs. We define the ideal characteristics of a novel thin filament-targeting compound and show its efficacy in-silico. We demonstrate that hybrid human-based hiPSC-CM and in-silico studies accelerate pathomechanism discovery and classification testing, improving clinical interpretation of genetic variants, and directing rational therapeutic targeting and design.

IntroductionChagas cardiomyopathy, a disease caused by Trypanosoma cruzi (T. cruzi) infection, is a major contributor to heart failure in Latin America. There are significant gaps in our understanding of the mechanism for infection of human cardiomyocytes, the pathways activated during the acute phase of the disease, and the molecular changes that lead to the progression of cardiomyopathy.MethodsTo investigate the effects of T. cruzi on human cardiomyocytes during infection, we infected induced pluripotent stem cell-derived cardiomyocytes (iPSC-CM) with the parasite and analyzed cellular, molecular, and metabolic responses at 3 hours, 24 hours, and 48 hours post infection (hpi) using transcriptomics (RNAseq), proteomics (LC-MS), and metabolomics (GC-MS and Seahorse) analyses.ResultsAnalyses of multiomic data revealed that cardiomyocyte infection caused a rapid increase in genes and proteins related to activation innate and adaptive immune systems and pathways, including alpha and gamma interferons, HIF-1α signaling, and glycolysis. These responses resemble prototypic responses observed in pathogen-activated immune cells. Infection also caused an activation of glycolysis that was dependent on HIF-1α signaling. Using gene editing and pharmacological inhibitors, we found that T. cruzi uptake was mediated in part by the glucose-facilitated transporter GLUT4 and that the attenuation of glycolysis, HIF-1α activation, or GLUT4 expression decreased T. cruzi infection. In contrast, pre-activation of pro-inflammatory immune responses with LPS resulted in increased infection rates.ConclusionThese findings suggest that T. cruzi exploits a HIF-1α-dependent, cardiomyocyte-intrinsic stress-response activation of glycolysis to promote intracellular infection and replication. These chronic immuno-metabolic responses by cardiomyocytes promote dysfunction, cell death, and the emergence of cardiomyopathy.

Background Chronic Chagas Cardiomyopathy (CCC) usually develops between 10 and 20 years after the first parasitic infection and is one of the leading causes of end-stage heart failure in Latin America. Despite the great inter-individual variability in CCC susceptibility (only 30% of infected individuals ever present CCC), there are no known predictors for disease development in those chronically infected. Methodology/Principal findings We describe a new susceptibility locus for CCC through a GWAS analysis in the SaMi-Trop cohort, a population-based study conducted in a Chagas endemic region from Brazil. This locus was also associated with CCC in the REDS II Study. The newly identified locus (rs34238187, OR 0.73, p-value 2.03 x 10−9) spans a haplotype of approximately 30Kb on chromosome 18 (chr18: 5028302–5057621) and is also associated with 80 different traits, most of them blood protein traits significantly enriched for immune-related biological pathways. Hi-C data show that the newly associated locus is able to interact with chromatin sites as far as 10Mb on chromosome 18 in a number of different cell types and tissues. Finally, we were able to confirm, at the tissue transcriptional level, the immune-associated blood protein signature using a multi-tissue differential gene expression and enrichment analysis. Conclusions/Significance We suggest that the newly identified locus impacts CCC risk among T cruzi infected individuals through the modulation of a downstream transcriptional and protein signature associated with host-parasite immune response. Functional characterization of the novel risk locus is warranted. Author summary Chronic Chagas Cardiomyopathy (CCC) usually develops between 10 and 20 years after the first parasitic infection and is one of the leading causes of end-stage heart failure in Latin America. Despite the great inter-individual variability in CCC susceptibility, there are no known predictors for disease development in those chronically infected. We describe a new susceptibility locus for CCC through a GWAS analysis in the SaMi-Trop cohort, a population-based study conducted in a Chagas endemic region from Brazil. The newly identified locus on chromosome 18 is also associated with 80 different traits, most of them blood protein traits significantly enriched for immune-related biological pathways. Finally, we were able to confirm, at the tissue transcriptional level, the immune-associated blood protein signature using a multi-tissue differential gene expression and enrichment analysis. We suggest that the newly identified locus impacts CCC risk among T cruzi infected individuals through the modulation of a downstream transcriptional and protein signature associated with host-parasite immune response.

Chagas disease is a tropical zoonosis caused by Trypanosoma cruzi. After infection, the host present an acute phase, usually asymptomatic, in which an extensive parasite proliferation and intense innate immune activity occurs, followed by a chronic phase, characterized by low parasitemia and development of specific immunity. Most individuals in the chronic phase remain without symptoms or organ damage, a state called indeterminate IND form. However, 20 to 40% of individuals develop cardiac or gastrointestinal complications at any time in life. Cardiomyocytes have an important role in the development of Chronic Chagas Cardiomyopathy (CCC) due to transcriptional and metabolic alterations that are crucial for the parasite survival and replication. However, it still not clear why some infected individuals progress to a cardiomyopathy phase, while others remain asymptomatic. In this work, we used hiPSCs-derived cardiomyocytes (hiPSC-CM) to investigate patterns of infection, proliferation and transcriptional response in IND and CCC patients. Our data show that T. cruzi infection and proliferation efficiency do not differ significantly in PBMCs and hiPSC-CM from both groups. However, RNA-seq analysis in hiPSC-CM infected for 24 hours showed a significantly different transcriptional response to the parasite in cells from IND or CCC patients. Cardiomyocytes from IND showed significant differences in the expression of genes related to antigen processing and presentation, as well as, immune co-stimulatory molecules. Furthermore, the downregulation of collagen production genes and extracellular matrix components was significantly different in these cells. Cardiomyocytes from CCC, in turn, showed increased expression of mTORC1 pathway and unfolded protein response genes, both associated to increased intracellular ROS production. These data point to a differential pattern of response, determined by baseline genetic differences between groups, which may have an impact on the development of a chronic outcome with or without the presentation of cardiac symptoms.

Background Myocardial iron deficiency (MID) in heart failure (HF) remains largely unexplored. We aim to establish defining criterion for MID, evaluate its pathophysiological role, and evaluate the applicability of monitoring it non‐invasively in human explanted hearts. Methods and Results Biventricular tissue iron levels were measured in both failing (n=138) and non‐failing control (NFC, n=46) explanted human hearts. Clinical phenotyping was complemented with comprehensive assessment of myocardial remodeling and mitochondrial functional profiles, including metabolic and oxidative stress. Myocardial iron status was further investigated by cardiac magnetic resonance imaging. Myocardial iron content in the left ventricle was lower in HF versus NFC (121.4 [88.1–150.3] versus 137.4 [109.2–165.9] μg/g dry weight), which was absent in the right ventricle. With a priori cutoff of 86.1 μg/g d.w. in left ventricle, we identified 23% of HF patients with MID (HF‐MID) associated with higher NYHA class and worsened left ventricle function. Respiratory chain and Krebs cycle enzymatic activities were suppressed and strongly correlated with depleted iron stores in HF‐MID hearts. Defenses against oxidative stress were severely impaired in association with worsened adverse remodeling in iron‐deficient hearts. Mechanistically, iron uptake pathways were impeded in HF‐MID including decreased translocation to the sarcolemma, while transmembrane fraction of ferroportin positively correlated with MID. Cardiac magnetic resonance with T2* effectively captured myocardial iron levels in failing hearts. Conclusions MID is highly prevalent in advanced human HF and exacerbates pathological remodeling in HF driven primarily by dysfunctional mitochondria and increased oxidative stress in the left ventricle. Cardiac magnetic resonance demonstrates clinical potential to non‐invasively monitor MID.

Hypertrophic cardiomyopathy is one of the most common inherited cardiomyopathies and a leading cause of sudden cardiac death in young adults. Despite profound insights into the genetics, there is imperfect correlation between mutation and clinical prognosis, suggesting complex molecular cascades driving pathogenesis. To investigate this, we performed an integrated quantitative multi-omics (proteomic, phosphoproteomic, and metabolomic) analysis to illuminate the early and direct consequences of mutations in myosin heavy chain in engineered human induced pluripotent stem-cell-derived cardiomyocytes relative to late-stage disease using patient myectomies. We captured hundreds of differential features, which map to distinct molecular mechanisms modulating mitochondrial homeostasis at the earliest stages of pathobiology, as well as stage-specific metabolic and excitation-coupling maladaptation. Collectively, this study fills in gaps from previous studies by expanding knowledge of the initial responses to mutations that protect cells against the early stress prior to contractile dysfunction and overt disease.

Abstract The structural and functional maturation of human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) is essential for pharmaceutical testing, disease modeling, and ultimately therapeutic use. Multicellular 3D-tissue platforms have improved the functional maturation of hiPSC-CMs, but probing cardiac contractile properties in a 3D environment remains challenging, especially at depth and in live tissues. Using small-angle X-ray scattering (SAXS) imaging, we show that hiPSC-CMs matured and examined in a 3D environment exhibit a periodic spatial arrangement of the myofilament lattice, which has not been previously detected in hiPSC-CMs. The contractile force is found to correlate with both the scattering intensity (R 2 = 0.44) and lattice spacing (R 2 = 0.46). The scattering intensity also correlates with lattice spacing (R 2 = 0.81), suggestive of lower noise in our structural measurement than in the functional measurement. Notably, we observed decreased myofilament ordering in tissues with a myofilament mutation known to lead to hypertrophic cardiomyopathy (HCM). Our results highlight the progress of human cardiac tissue engineering and enable unprecedented study of structural maturation in hiPSC-CMs.

Abstract 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.

Summary: Childhood-onset myocardial hypertrophy and cardiomyopathic changes are associated with significant morbidity and mortality in early life, particularly in patients with Noonan syndrome, a multisystemic genetic disorder caused by autosomal dominant mutations in genes of the Ras-MAPK pathway. Although the cardiomyopathy associated with Noonan syndrome (NS-CM) shares certain cardiac features with the hypertrophic cardiomyopathy caused by mutations in sarcomeric proteins (HCM), such as pathological myocardial remodeling, ventricular dysfunction, and increased risk for malignant arrhythmias, the clinical course of NS-CM significantly differs from HCM. This suggests a distinct pathophysiology that remains to be elucidated. Here, through analysis of sarcomeric myosin conformational states, histopathology, and gene expression in left ventricular myocardial tissue from NS-CM, HCM, and normal hearts complemented with disease modeling in cardiomyocytes differentiated from patient-derived PTPN11N308S/+ induced pluripotent stem cells, we demonstrate distinct disease phenotypes between NS-CM and HCM and uncover cell cycle defects as a potential driver of NS-CM.

Background Human mutations in the X‐linked lysosome‐associated membrane protein‐2 (LAMP2) gene can cause a multisystem Danon disease or a primary cardiomyopathy characterized by massive hypertrophy, conduction system abnormalities, and malignant ventricular arrhythmias. We introduced an in‐frame LAMP2 gene exon 6 deletion mutation (denoted L2Δ6) causing human cardiomyopathy, into mouse LAMP2 gene, to elucidate its consequences on cardiomyocyte biology. This mutation results in in‐frame deletion of 41 amino acids, compatible with presence of some defective LAMP2 protein. Methods and Results Left ventricular tissues from L2Δ6 and wild‐type mice had equivalent amounts of LAMP2 RNA, but a significantly lower level of LAMP2 protein. By 20 weeks of age male mutant mice developed left ventricular hypertrophy which was followed by left ventricular dilatation and reduced systolic function. Cardiac electrophysiology and isolated cardiomyocyte studies demonstrated ventricular arrhythmia, conduction disturbances, abnormal calcium transients and increased sensitivity to catecholamines. Myocardial fibrosis was strikingly increased in 40‐week‐old L2Δ6 mice, recapitulating findings of human LAMP2 cardiomyopathy. Immunofluorescence and transmission electron microscopy identified mislocalization of lysosomes and accumulation of autophagosomes between sarcomeres, causing profound morphological changes disrupting the cellular ultrastructure. Transcription profile and protein expression analyses of L2Δ6 hearts showed significantly increased expression of genes encoding activators and protein components of autophagy, hypertrophy, and apoptosis. Conclusions We suggest that impaired autophagy results in cardiac hypertrophy and profound transcriptional reactions that impacted metabolism, calcium homeostasis, and cell survival. These responses define the molecular pathways that underlie the pathology and aberrant electrophysiology in cardiomyopathy of Danon disease.

LMNA is one of the leading causative genes of genetically inherited dilated cardiomyopathy (DCM). Unlike most DCM-causative genes, which encode sarcomeric or sarcomere-related proteins, LMNA encodes nuclear envelope proteins, lamin A and C, and does not directly associate with contractile function. However, a mutation in this gene could lead to the development of DCM. The molecular mechanism of how LMNA mutation contributes to DCM development remains largely unclear and yet to be elucidated. The objective of this study was to clarify the mechanism of developing DCM caused by LMNA mutation.Methods and Results: We assessed cardiomyocyte phenotypes and characteristics focusing on cell cycle activity in mice with Lmna mutation. Both cell number and cell size were reduced, cardiomyocytes were immature, and cell cycle activity was retarded in Lmna mutant mice at both 5 weeks and 2 years of age. RNA-sequencing and pathway analysis revealed “proliferation of cells” had the most substantial impact on Lmna mutant mice. Cdkn1a, which encodes the cell cycle regulating protein p21, was strongly upregulated in Lmna mutants, and upregulation of p21 was confirmed by Western blot and immunostaining. DNA damage, which is known to upregulate Cdkn1a, was more abundantly detected in Lmna mutant mice. To assess the proliferative capacity of cardiomyocytes, the apex of the neonate mouse heart was resected, and recovery from the insult was observed. A restricted cardiomyocyte proliferating capacity after resecting the apex of the heart was observed in Lmna mutant mice.Conclusions: Our results strongly suggest that loss of lamin function contributes to impaired cell proliferation through cell cycle defects. The inadequate inborn or responsive cell proliferation capacity plays an essential role in developing DCM with LMNA mutation.

Large whole-exome sequencing studies have suggested that the genetic architecture of syndromic congenital heart disease (CHD) is different from sporadic forms. Here, Watkins et al. estimate the relative contribution of damaging recessive and de novo genotypes to CHD in 2391 trios and find them to be associated with different gene functions.

Advanced paternal age associates with increased risk for psychiatric and developmental disorders in offspring. Here, Taylor et al. utilize parent-child trio exome sequencing data sets to estimate the contribution of paternal age-related de novo mutations to multiple disorders, including heart disease and schizophrenia.

ObjectiveThe myosin-binding protein C (MYBPC3) c.927-2A>G founder mutation accounts for >90% of sarcomeric hypertrophic cardiomyopathy (HCM) in Iceland. This cross-sectional observational study explored the penetrance and phenotypic burden among carriers of this single, prevalent founder mutation.MethodsWe studied 60 probands with HCM caused by MYBPC3 c.927-2A>G and 225 first-degree relatives. All participants underwent comprehensive clinical evaluation and relatives were genotyped.ResultsGenetic and clinical evaluation of relatives identified 49 genotype-positive (G+) relatives with left ventricular hypertrophy (G+/LVH+), 59 G+without LVH (G+/LVH−) and 117 genotype-negative relatives (unaffected). Compared with HCM probands, G+/LVH+ relatives were older at HCM diagnosis, had less LVH, a less prevalent diastolic dysfunction, fewer ECG abnormalities, lower serum N-terminal pro-B-type natriuretic peptide (NT-proBNP) and high-sensitivity cardiac troponin I levels, and fewer symptoms. The penetrance of HCM was influenced by age and sex; specifically, LVH was present in 39% of G+males but only 9% of G+females under age 40 years (p=0.015), versus 86% and 83%, respectively, after age 60 (p=0.89). G+/LVH− subjects had normal wall thicknesses, diastolic function and NT-proBNP levels, but subtle changes in LV geometry and more ECG abnormalities than their unaffected relatives.ConclusionsPhenotypic expression of the Icelandic MYBPC3 founder mutation varies by age, sex and proband status. Men are more likely to have LVH at a younger age, and disease manifestations were more prominent in probands than in relatives identified via family screening. G+/LVH− individuals had subtle clinical differences from unaffected relatives well into adulthood, indicating subclinical phenotypic expression of the pathogenic mutation.

The processes regulating cardiomyocyte (CM) maturation are unclear. Here, the authors show that serum response factor regulates CM maturation only in neonatal CMs through stage-specific chromatin occupancy that affects cell size, sarcomere and transverse-tubule organization, and mitochondria

Summary: Impulses generated by a multicellular, bioelectric signaling center termed the sinoatrial node (SAN) stimulate the rhythmic contraction of the heart. The SAN consists of a network of electrochemically oscillating pacemaker cells encased in a heterogeneous connective tissue microenvironment. Although the cellular composition of the SAN has been a point of interest for more than a century, the biological processes that drive the tissue-level assembly of the cells within the SAN are unknown. Here, we demonstrate that the SAN’s structural features result from a developmental process during which mesenchymal cells derived from a multipotent progenitor structure, the proepicardium, integrate with and surround pacemaker myocardium. This process actively remodels the forming SAN and is necessary for sustained electrogenic signal generation and propagation. Collectively, these findings provide experimental evidence for how the microenvironmental architecture of the SAN is patterned and demonstrate that proper cellular arrangement is critical for cardiac pacemaker biorhythmicity. : How the higher order, tissue-level architecture of the heart’s pacemaker region is patterned during development is poorly understood. Bressan et al. demonstrate that a cellular remodeling process, which depends on the integration of pacemaker muscle with proepicardial-derived mesenchymal cells, establishes the microenvironmental conditions required for rhythmic stimulation of the heart. Keywords: sinoatrial node, cardiac pacemaker cell, cardiac morphogenesis, proepicardium

Abstract Background Omentin-1, also known as Intelectin-1 (ITLN1), is an adipokine with plasma levels associated with diabetes, obesity, and coronary artery disease. Recent studies suggest that ITLN1 can mitigate myocardial ischemic injury but the expression of ITLN1 in the heart itself has not been well characterized. The purpose of this study is to discern the relationship between the expression pattern of ITLN1 RNA in the human heart and the level of circulating ITLN1 protein in plasma from the same patients following myocardial ischemia. Methods A large cohort of patients (n = 140) undergoing elective cardiac surgery for aortic valve replacement were enrolled in this study. Plasma and left ventricular biopsy samples were taken at the beginning of cardiopulmonary bypass and after an average of 82 min of ischemic cross clamp time. The localization of ITLN1 in epicardial adipose tissue (EAT) was also further characterized with immunoassays and cell fate transition studies. Results mRNA expression of ITLN1 decreases in left ventricular tissue after acute ischemia in human patients (mean difference 280.48, p = 0.001) whereas plasma protein levels of ITLN1 increase (mean difference 5.24, p

Summary: The degradation and release of cardiac myosin binding protein-C (cMyBP-C) upon cardiac damage may stimulate an inflammatory response and autoantibody (AAb) production. We determined whether the presence of cMyBP-C-AAbs associated with adverse cardiac function in cardiovascular disease patients. Importantly, cMyBP-C-AAbs were significantly detected in acute coronary syndrome patient sera upon arrival to the emergency department, particularly in ST-segment elevation myocardial infarction patients. Patients positive for cMyBP-C-AAbs had reduced left ventricular ejection fraction and elevated levels of clinical biomarkers of myocardial infarction. We conclude that cMyBP-C-AAbs may serve as early predictive indicators of deteriorating cardiac function and patient outcome in acute coronary syndrome patients prior to the infarction. Key Words: acute myocardial infarction, autoantibodies, cardiac myosin binding protein-c, cardiomyopathy

AMP-activated protein kinase (AMPK) is a metabolic enzyme that can be activated by nutrient stress or genetic mutations. Missense mutations in the regulatory subunit, PRKAG2, activate AMPK and cause left ventricular hypertrophy, glycogen accumulation, and ventricular pre-excitation. Using human iPS cell models combined with three-dimensional cardiac microtissues, we show that activating PRKAG2 mutations increase microtissue twitch force by enhancing myocyte survival. Integrating RNA sequencing with metabolomics, PRKAG2 mutations that activate AMPK remodeled global metabolism by regulating RNA transcripts to favor glycogen storage and oxidative metabolism instead of glycolysis. As in patients with PRKAG2 cardiomyopathy, iPS cell and mouse models are protected from cardiac fibrosis, and we define a crosstalk between AMPK and post-transcriptional regulation of TGFβ isoform signaling that has implications in fibrotic forms of cardiomyopathy. Our results establish critical connections among metabolic sensing, myocyte survival, and TGFβ signaling.

Abstract Background Individuals with hypertrophic cardiomyopathy (HCM), even when asymptomatic, are at‐risk for sudden cardiac death and stroke from arrhythmias, making it imperative to identify individuals affected by this familial disorder. Consensus guidelines recommend that first‐degree relatives (FDRs) of a person with HCM undergo serial cardiovascular evaluations. Methods We determined the uptake of family screening in patients with HCM and developed an online video intervention to facilitate family communication and screening. Family screening and genetic testing data were collected through a prospective quality improvement initiative, a standardized clinical assessment and management plan (SCAMP), utilized in an established cardiovascular genetics clinic. Patients were prescribed an online video if screening of their FDRs was incomplete and a pilot study on video utilization and family communication was conducted. Results Two‐hundred and sixteen probands with HCM were enrolled in SCAMP Phase I and 190 were enrolled in SCAMP Phase II. In both phases, probands reported that 51% of FDRs had been screened (382/749 in Phase I, 258/504 in Phase II). Twenty patients participated in a pilot study on video utilization and family communication. Nine participants reported watching the video and six participants reported sharing the video with relatives; however only one participant reported sharing the video with relatives who were not yet aware of the diagnosis of HCM in the family. Conclusion Despite care in a specialized cardiovascular genetics clinic, approximately one half of FDRs of patients with HCM remained unscreened. Online interventions and videos may serve as supplemental tools for patients communicating genetic risk information to relatives.

Congenital heart disease (CHD) is a disorder that occurs in ∼1% of live births. Here the authors describe a genome-wide allele-specific expression analyses in CHD patients, identifying five new genes involved in CHD and showing that paternally-expressed imprinted genes are monoallelic, while maternally-expressed imprinted genes are biallelic.

Heterozygous mice (αMHC403/+) expressing the human disease-causing mutation Arg403Gln exhibit cardinal features of hypertrophic cardiomyopathy (HCM) including hypertrophy, myocyte disarray, and increased myocardial fibrosis. Treatment of αMHC403/+mice with the L-type calcium channel (ICa-L) antagonist diltiazem has been shown to decrease left ventricular anterior wall thickness, cardiac myocyte hypertrophy, disarray, and fibrosis. However, the role of the ICa-L in the development of HCM is not known. In addition to maintaining cardiac excitation and contraction in myocytes, the ICa-L also regulates mitochondrial function through transmission of movement of ICa-L via cytoskeletal proteins to mitochondrial voltage-dependent anion channel. Here, the authors investigated the role of ICa-L in regulating mitochondrial function in αMHC403/+mice. Whole-cell patch clamp studies showed that ICa-L current inactivation kinetics were significantly increased in αMHC403/+cardiac myocytes, but that current density and channel expression were similar to wild-type cardiac myocytes. Activation of ICa-L caused a significantly greater increase in mitochondrial membrane potential and metabolic activity in αMHC403/+. These increases were attenuated with ICa-L antagonists and following F-actin or β-tubulin depolymerization. The authors observed increased levels of fibroblast growth factor-21 in αMHC403/+mice, and altered mitochondrial DNA copy number consistent with altered mitochondrial activity and the development of cardiomyopathy. These studies suggest that the Arg403Gln mutation leads to altered functional communication between ICa-L and mitochondria that is associated with increased metabolic activity, which may contribute to the development of cardiomyopathy. ICa-L antagonists may be effective in reducing the cardiomyopathy in HCM by altering metabolic activity.

How chromatin-remodeling complexes modulate gene networks to control organ-specific properties is not well understood. For example, Baf60c (Smarcd3) encodes a cardiac-enriched subunit of the SWI/SNF-like BAF chromatin complex, but its role in heart development is not fully understood. We found that constitutive loss of Baf60c leads to embryonic cardiac hypoplasia and pronounced cardiac dysfunction. Conditional deletion of Baf60c in cardiomyocytes resulted in postnatal dilated cardiomyopathy with impaired contractile function. Baf60c regulates a gene expression program that includes genes encoding contractile proteins, modulators of sarcomere function, and cardiac metabolic genes. Many of the genes deregulated in Baf60c null embryos are targets of the MEF2/SRF co-factor Myocardin (MYOCD). In a yeast two-hybrid screen, we identified MYOCD as a BAF60c interacting factor; we showed that BAF60c and MYOCD directly and functionally interact. We conclude that Baf60c is essential for coordinating a program of gene expression that regulates the fundamental functional properties of cardiomyocytes.

Despite an obvious role for consanguinity in congenital heart disease (CHD), most studies fail to document a monogenic model of inheritance except for few cases. We hereby describe a first-degree cousins consanguineous Lebanese family with 7 conceived children: 2 died in utero of unknown causes, 3 have CHD, and 4 have polydactyly. The aim of the study is to unveil the genetic variant(s) causing these phenotypes using next generation sequencing (NGS) technology. Targeted exome sequencing identified a heterozygous duplication in CSRP1 which leads to a potential frameshift mutation at position 154 of the protein. This mutation is inherited from the father, and segregates only with the CHD phenotype. The in vitro characterization demonstrates that the mutation dramatically abrogates its transcriptional activity over cardiac promoters like NPPA. In addition, it differentially inhibits the physical association of CSRP1 with SRF, GATA4, and with the newly described partner herein TBX5. Whole exome sequencing failed to show any potential variant linked to polydactyly, but revealed a novel missense mutation in TRPS1. This mutation is inherited from the healthy mother, and segregating only with the cardiac phenotype. Both TRPS1 and CSRP1 physically interact, and the mutations in each abrogate their partnership. Our findings add fundamental knowledge into the molecular basis of CHD, and propose the di-genic model of inheritance as responsible for such malformations.

Ryr1I4895T/WT (IT/+) mice express a knockin mutation corresponding to the human I4898T EC-uncoupling mutation in the type 1 ryanodine receptor/Ca2+ release channel (RyR1), which causes a severe form of central core disease ((CD). IT/+ mice exhibit a slowly progressive congenital myopathy, with neonatal respiratory stress, skeletal muscle weakness, impaired mobility, dorsal kyphosis, and hind limb paralysis. Lesions observed in myofibers from diseased mice undergo age-dependent transformation from minicores to cores and nemaline rods. Early ultrastructural abnormalities include sarcomeric misalignment, Z-line streaming, focal loss of cross-striations, and myofibrillar splitting and intermingling that may arise from defective myofibrillogenesis. However, manifestation of the disease phenotype is highly variable on a Sv129 genomic background. Quantitative RT-PCR shows an equimolar ratio of WT and mutant Ryrl transcripts within lT/+ myofibers and total RyR1 protein expression levels are normal. We propose a unifying theory in which the cause of core formation lies in functional heterogeneity among RyR1 tetramers. Random combinations of normal and either leaky or EC-uncoupled RyR subunits would lead to spatial differences in Ca2+ transients; the resulting heterogeneity of contraction among myofibrils would lead to focal, irreversible tearing and shearing, which would, over time, enlarge to form minicores, cores, and nemaline rods. The IT/+ mouse line is proposed to be a valid model of RyR1-related congenital myopathy, offering high potential for elucidation of the pathogenesis of skeletal muscle disorders arising from impaired EC coupling. [ABSTRACT FROM AUTHOR]

Although DNA methylation is the best characterized epigenetic mark, the mechanism by which it is targeted to specific regions in the genome remains unclear. Recent studies have revealed that local DNA methylation profiles might be dictated by cis-regulatory DNA sequences that mainly operate via DNA-binding factors. Consistent with this finding, we have recently shown that disruption of CTCF-binding sites by rare single nucleotide variants (SNVs) can underlie cis-linked DNA methylation changes in patients with congenital anomalies. These data raise the hypothesis that rare genetic variation at transcription factor binding sites (TFBSs) might contribute to local DNA methylation patterning. In this work, by combining blood genome-wide DNA methylation profiles, whole genome sequencing-derived SNVs from 247 unrelated individuals along with 133 predicted TFBS motifs derived from ENCODE ChIP-Seq data, we observed an association between the disruption of binding sites for multiple TFs by rare SNVs and extreme DNA methylation values at both local and, to a lesser extent, distant CpGs. While the majority of these changes affected only single CpGs, 24% were associated with multiple outlier CpGs within ±1kb of the disrupted TFBS. Interestingly, disruption of functionally constrained sites within TF motifs lead to larger DNA methylation changes at nearby CpG sites. Altogether, these findings suggest that rare SNVs at TFBS negatively influence TF-DNA binding, which can lead to an altered local DNA methylation profile. Furthermore, subsequent integration of DNA methylation and RNA-Seq profiles from cardiac tissues enabled us to observe an association between rare SNV-directed DNA methylation and outlier expression of nearby genes. In conclusion, our findings not only provide insights into the effect of rare genetic variation at TFBS on shaping local DNA methylation and its consequences on genome regulation, but also provide a rationale to incorporate DNA methylation data to interpret the functional role of rare variants.

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.

[This corrects the article DOI: 10.1371/journal.pone.0172995.].

Deep nucleotide sequencing enables the unbiased, broad-spectrum detection of viruses in clinical samples without requiring an a priori hypothesis for the source of infection. However, its use in clinical research applications is limited by low cost-effectiveness given that most of the sequencing information from clinical samples is related to the human genome, which renders the analysis of viral genomes challenging. To overcome this limitation we developed ViroFind, an in-solution target-enrichment platform for virus detection and discovery in clinical samples. ViroFind comprises 165,433 viral probes that cover the genomes of 535 selected DNA and RNA viruses that infect humans or could cause zoonosis. The ViroFind probes are used in a hybridization reaction to enrich viral sequences and therefore enhance the detection of viral genomes via deep sequencing. We used ViroFind to detect and analyze all viral populations in the brain of 5 patients with progressive multifocal leukoencephalopathy (PML) and of 18 control subjects with no known neurological disease. Compared to direct deep sequencing, by using ViroFind we enriched viral sequences present in the clinical samples up to 127-fold. We discovered highly complex polyoma virus JC populations in the PML brain samples with a remarkable degree of genetic divergence among the JC virus variants of each PML brain sample. Specifically for the viral capsid protein VP1 gene, we identified 24 single nucleotide substitutions, 12 of which were associated with amino acid changes. The most frequent (4 of 5 samples, 80%) amino acid change was D66H, which is associated with enhanced tissue tropism, and hence likely a viral fitness advantage, compared to other variants. Lastly, we also detected sparse JC virus sequences in 10 of 18 (55.5%) of control samples and sparse human herpes virus 6B (HHV6B) sequences in the brain of 11 of 18 (61.1%) control subjects. In sum, ViroFind enabled the in-depth analysis of all viral genomes in PML and control brain samples and allowed us to demonstrate a high degree of JC virus genetic divergence in vivo that has been previously underappreciated. ViroFind can be used to investigate the structure of the virome with unprecedented depth in health and disease state.

Cardiac β-myosin variants cause hypertrophic (HCM) or dilated (DCM) cardiomyopathy by disrupting sarcomere contraction and relaxation. The locations of variants on isolated myosin head structures predict contractility effects but not the prominent relaxation and energetic deficits that characterize HCM. During relaxation, pairs of myosins form interacting-heads motif (IHM) structures that with other sarcomere proteins establish an energy-saving, super-relaxed (SRX) state. Using a human β-cardiac myosin IHM quasi-atomic model, we defined interactions sites between adjacent myosin heads and associated protein partners, and then analyzed rare variants from 6112 HCM and 1315 DCM patients and 33,370 ExAC controls. HCM variants, 72% that changed electrostatic charges, disproportionately altered IHM interaction residues (expected 23%; HCM 54%, p=2.6×10−19; DCM 26%, p=0.66; controls 20%, p=0.23). HCM variant locations predict impaired IHM formation and stability, and attenuation of the SRX state - accounting for altered contractility, reduced diastolic relaxation, and increased energy consumption, that fully characterizes HCM pathogenesis.

AimsDilated cardiomyopathy (DCM) is an important cause of heart failure with a strong familial component. We performed an exome-wide array-based association study (EWAS) to assess the contribution of missense variants to sporadic DCM.Methods and results116,855 single nucleotide variants (SNVs) were analyzed in 2796 DCM patients and 6877 control subjects from 6 populations of European ancestry. We confirmed two previously identified associations with SNVs in BAG3 and ZBTB17 and discovered six novel DCM-associated loci (Q-valueConclusionWe identified eight loci independently associated with sporadic DCM. The functions of the best candidate genes at these loci suggest that proteostasis regulation might play a role in DCM pathophysiology.

The epicardium plays an important role in coronary vessel formation and Tgfbr3-/- mice exhibit failed coronary vessel development associated with decreased epicardial cell invasion. Immortalized Tgfbr3-/- epicardial cells display the same defects. Tgfbr3+/+ and Tgfbr3-/- cells incubated for 72 hours with VEH or ligands known to promote invasion via TGFβR3 (TGFβ1, TGFβ2, BMP2), for 72 hours were harvested for RNA-seq analysis. We selected for genes >2-fold differentially expressed between Tgfbr3+/+ and Tgfbr3-/- cells when incubated with VEH (604), TGFβ1 (515), TGFβ2 (553), or BMP2 (632). Gene Ontology (GO) analysis of these genes identified dysregulated biological processes consistent with the defects observed in Tgfbr3-/- cells, including those associated with extracellular matrix interaction. GO and Gene Regulatory Network (GRN) analysis identified distinct expression profiles between TGFβ1-TGFβ2 and VEH-BMP2 incubated cells, consistent with the differential response of epicardial cells to these ligands in vitro. Despite the differences observed between Tgfbr3+/+ and Tgfbr3-/- cells after TGFβ and BMP ligand addition, GRNs constructed from these gene lists identified NF-ĸB as a key nodal point for all ligands examined. Tgfbr3-/- cells exhibited decreased expression of genes known to be activated by NF-ĸB signaling. NF-ĸB activity was stimulated in Tgfbr3+/+ epicardial cells after TGFβ2 or BMP2 incubation, while Tgfbr3-/- cells failed to activate NF-ĸB in response to these ligands. Tgfbr3+/+ epicardial cells incubated with an inhibitor of NF-ĸB signaling no longer invaded into a collagen gel in response to TGFβ2 or BMP2. These data suggest that NF-ĸB signaling is dysregulated in Tgfbr3-/- epicardial cells and that NF-ĸB signaling is required for epicardial cell invasion in vitro. Our approach successfully identified a signaling pathway important in epicardial cell behavior downstream of TGFβR3. Overall, the genes and signaling pathways identified through our analysis yield the first comprehensive list of candidate genes whose expression is dependent on TGFβR3 signaling.

Abstract Aberrant organ development is associated with a wide spectrum of disorders, from schizophrenia to congenital heart disease, but systems‐level insight into the underlying processes is very limited. Using heart morphogenesis as general model for dissecting the functional architecture of organ development, we combined detailed phenotype information from deleterious mutations in 255 genes with high‐confidence experimental interactome data, and coupled the results to thorough experimental validation. Hereby, we made the first systematic analysis of spatio‐temporal protein networks driving many stages of a developing organ identifying several novel signaling modules. Our results show that organ development relies on surprisingly few, extensively recycled, protein modules that integrate into complex higher‐order networks. This design allows the formation of a complicated organ using simple building blocks, and suggests how mutations in the same genes can lead to diverse phenotypes. We observe a striking temporal correlation between organ complexity and the number of discrete functional modules coordinating morphogenesis. Our analysis elucidates the organization and composition of spatio‐temporal protein networks that drive the formation of organs, which in the future may lay the foundation of novel approaches in treatments, diagnostics, and regenerative medicine.

Postoperative ventricular dysfunction (VnD) occurs in 9-20% of coronary artery bypass graft (CABG) surgical patients and is associated with increased postoperative morbidity and mortality. Understanding genetic causes of postoperative VnD should enhance patient risk stratification and improve treatment and prevention strategies. We aimed to determine if genetic variants associate with occurrence of in-hospital VnD after CABG surgery.A genome-wide association study identified single nucleotide polymorphisms (SNPs) associated with postoperative VnD in male subjects of European ancestry undergoing isolated primary CABG surgery with cardiopulmonary bypass. VnD was defined as the need for ≥2 inotropes or mechanical ventricular support after CABG surgery. Validated SNPs were assessed further in two replication CABG cohorts and meta-analysis was performed.Over 100 SNPs were associated with VnD (P2.1) of developing in-hospital VnD after CABG surgery. However, three genetic loci identified by meta-analysis were more modestly associated with development of postoperative VnD. Studies of larger cohorts to assess these loci as well as to define other genetic mechanisms and related biology that link genetic variants to postoperative ventricular dysfunction are warranted.

Hypertrophic (HCM) and dilated (DCM) cardiomyopathies result from sarcomeric protein mutations, including cardiac troponin T (cTnT, TNNT2). We determined whether TNNT2 mutations cause cardiomyopathies by altering cTnT function or quantity; whether the severity of DCM is related to the ratio of mutant to wildtype cTnT; whether Ca(2+) desensitization occurs in DCM; and whether absence of cTnT impairs early embryonic cardiogenesis.We ablated Tnnt2 to produce heterozygous Tnnt2(+/-) mice, and crossbreeding produced homozygous null Tnnt2(-/-) embryos. We also generated transgenic mice overexpressing wildtype (TG(WT)) or DCM mutant (TG(K210Delta)) Tnnt2. Crossbreeding produced mice lacking one allele of Tnnt2, but carrying wildtype (Tnnt2(+/-)/TG(WT)) or mutant (Tnnt2(+/-)/TG(K210Delta)) transgenes. Tnnt2(+/-) mice relative to wildtype had significantly reduced transcript (0.82+/-0.06[SD] vs. 1.00+/-0.12 arbitrary units; p = 0.025), but not protein (1.01+/-0.20 vs. 1.00+/-0.13 arbitrary units; p = 0.44). Tnnt2(+/-) mice had normal hearts (histology, mass, left ventricular end diastolic diameter [LVEDD], fractional shortening [FS]). Moreover, whereas Tnnt2(+/-)/TG(K210Delta) mice had severe DCM, TG(K210Delta) mice had only mild DCM (FS 18+/-4 vs. 29+/-7%; p