Search

BackgroundHeart failure with preserved ejection fraction (HFpEF) is an emerging major unmet need and one of the most significant clinic challenges in cardiology. The pathogenesis of HFpEF is associated with multiple risk factors. Hypertension and metabolic disorders associated with obesity are the 2 most prominent comorbidities observed in patients with HFpEF. Although hypertension-induced mechanical overload has long been recognized as a potent contributor to heart failure with reduced ejection fraction, the synergistic interaction between mechanical overload and metabolic disorders in the pathogenesis of HFpEF remains poorly characterized.MethodWe investigated the functional outcome and the underlying mechanisms from concurrent mechanic and metabolic stresses in the heart by applying transverse aortic constriction in lean C57Bl/6J or obese/diabetic B6.Cg-Lepob/J (ob/ob) mice, followed by single-nuclei RNA-seq and targeted manipulation of a top-ranked signaling pathway differentially affected in the 2 experimental cohorts.ResultsIn contrast to the post-trans-aortic constriction C57Bl/6J lean mice, which developed pathological features of heart failure with reduced ejection fraction over time, the post-trans-aortic constriction ob/ob mice showed no significant changes in ejection fraction but developed characteristic pathological features of HFpEF, including diastolic dysfunction, worsened cardiac hypertrophy, and pathological remodeling, along with further deterioration of exercise intolerance. Single-nuclei RNA-seq analysis revealed significant transcriptome reprogramming in the cardiomyocytes stressed by both pressure overload and obesity/diabetes, markedly distinct from the cardiomyocytes singularly stressed by pressure overload or obesity/diabetes. Furthermore, glucagon signaling was identified as the top-ranked signaling pathway affected in the cardiomyocytes associated with HFpEF. Treatment with a glucagon receptor antagonist significantly ameliorated the progression of HFpEF-related pathological features in 2 independent preclinical models. Importantly, cardiomyocyte-specific genetic deletion of the glucagon receptor also significantly improved cardiac function in response to pressure overload and metabolic stress.ConclusionsThese findings identify glucagon receptor signaling in cardiomyocytes as a critical determinant of HFpEF progression and provide proof-of-concept support for glucagon receptor antagonism as a potential therapy for the disease.

Parallel to major changes in fatty acid and glucose metabolism, defect in branched-chain amino acid (BCAA) catabolism has also been recognized as a metabolic hallmark and potential therapeutic target for heart failure. However, BCAA catabolic enzymes are ubiquitously expressed in all cell types and a systemic BCAA catabolic defect is also manifested in metabolic disorder associated with obesity and diabetes. Therefore, it remains to be determined the cell-autonomous impact of BCAA catabolic defect in cardiomyocytes in intact hearts independent from its potential global effects. In this study, we developed two mouse models. One is cardiomyocyte and temporal-specific inactivation of the E1α subunit (BCKDHA-cKO) of the branched-chain α-ketoacid dehydrogenase (BCKDH) complex, which blocks BCAA catabolism. Another model is cardiomyocyte specific inactivation of the BCKDH kinase (BCKDK-cKO), which promotes BCAA catabolism by constitutively activating BCKDH activity in adult cardiomyocytes. Functional and molecular characterizations showed E1α inactivation in cardiomyocytes was sufficient to induce loss of cardiac function, systolic chamber dilation and pathological transcriptome reprogramming. On the other hand, inactivation of BCKDK in intact heart does not have an impact on baseline cardiac function or cardiac dysfunction under pressure overload. Our results for the first time established the cardiomyocyte cell autonomous role of BCAA catabolism in cardiac physiology. These mouse lines will serve as valuable model systems to investigate the underlying mechanisms of BCAA catabolic defect induced heart failure and to provide potential insights for BCAA targeted therapy.

Heart failure with preserved ejection fraction (HFpEF) exhibits a sex bias, being more common in women than men, and we hypothesize that mitochondrial sex differences might underlie this bias. As part of genetic studies of heart failure in mice, we observe that heart mitochondrial DNA levels and function tend to be reduced in females as compared to males. We also observe that expression of genes encoding mitochondrial proteins are higher in males than females in human cohorts. We test our hypothesis in a panel of genetically diverse inbred strains of mice, termed the Hybrid Mouse Diversity Panel (HMDP). Indeed, we find that mitochondrial gene expression is highly correlated with diastolic function, a key trait in HFpEF. Consistent with this, studies of a "two-hit" mouse model of HFpEF confirm that mitochondrial function differs between sexes and is strongly associated with a number of HFpEF traits. By integrating data from human heart failure and the mouse HMDP cohort, we identify the mitochondrial gene Acsl6 as a genetic determinant of diastolic function. We validate its role in HFpEF using adenoviral over-expression in the heart. We conclude that sex differences in mitochondrial function underlie, in part, the sex bias in diastolic function.

BackgroundHuman-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have emerged as a promising tool for disease modeling and drug development. However, hiPSC-CMs remain functionally immature, which hinders their utility as a model of human cardiomyocytes.ObjectiveTo improve the electrophysiological maturation of hiPSC-CMs.Methods and resultsOn day 16 of cardiac differentiation, hiPSC-CMs were treated with 100 nmol/L triiodothyronine (T3) and 1 μmol/L Dexamethasone (Dex) or vehicle for 14 days. On day 30, vehicle- and T3 + Dex-treated hiPSC-CMs were dissociated and replated either as cell sheets or single cells. Optical mapping and patch-clamp technique were used to examine the electrophysiological properties of vehicle- and T3 + Dex-treated hiPSC-CMs. Compared to vehicle, T3 + Dex-treated hiPSC-CMs had a slower spontaneous beating rate, more hyperpolarized resting membrane potential, faster maximal upstroke velocity, and shorter action potential duration. Changes in spontaneous activity and action potential were mediated by decreased hyperpolarization-activated current (If) and increased inward rectifier potassium currents (IK1), sodium currents (INa), and the rapidly and slowly activating delayed rectifier potassium currents (IKr and IKs, respectively). Furthermore, T3 + Dex-treated hiPSC-CM cell sheets (hiPSC-CCSs) exhibited a faster conduction velocity and shorter action potential duration than the vehicle. Inhibition of IK1 by 100 μM BaCl2 significantly slowed conduction velocity and prolonged action potential duration in T3 + Dex-treated hiPSC-CCSs but had no effect in the vehicle group, demonstrating the importance of IK1 for conduction velocity and action potential duration.ConclusionT3 + Dex treatment is an effective approach to rapidly enhance electrophysiological maturation of hiPSC-CMs.

Individuals with heart failure (HF) frequently present with comorbidities, including obesity, insulin resistance, hypertension, and dyslipidemia. Many patients with HF experience cardiogenic dementia, yet the pathophysiology of this disease remains poorly understood. Using a swine model of cardiometabolic HF (Western diet+aortic banding; WD-AB), we tested the hypothesis that WD-AB would promote a multidementia phenotype involving cerebrovascular dysfunction alongside evidence of Alzheimer's disease (AD) pathology. The results provide evidence of cerebrovascular insufficiency coupled with neuroinflammation and amyloidosis in swine with experimental cardiometabolic HF. Although cardiac ejection fraction was normal, indices of arterial compliance and cerebral blood flow were reduced, and cerebrovascular regulation was impaired in the WD-AB group. Cerebrovascular dysfunction occurred concomitantly with increased MAPK signaling and amyloidogenic processing (i.e., increased APP, BACE1, CTF, and Aβ40 in the prefrontal cortex and hippocampus) in the WD-AB group. Transcriptomic profiles of the stellate ganglia revealed the WD-AB group displayed an enrichment of gene networks associated with MAPK/ERK signaling, AD, frontotemporal dementia, and a number of behavioral phenotypes implicated in cognitive impairment. These provide potentially novel evidence from a swine model that cerebrovascular and neuronal pathologies likely both contribute to the dementia profile in a setting of cardiometabolic HF.

BackgroundThe majority of quantitative genetic models used to map complex traits assume that alleles have similar effects across all individuals. Significant evidence suggests, however, that epistatic interactions modulate the impact of many alleles. Nevertheless, identifying epistatic interactions remains computationally and statistically challenging. In this work, we address some of these challenges by developing a statistical test for polygenic epistasis that determines whether the effect of an allele is altered by the global genetic ancestry proportion from distinct progenitors.ResultsWe applied our method to data from mice and yeast. For the mice, we observed 49 significant genotype-by-ancestry interaction associations across 14 phenotypes as well as over 1,400 Bonferroni-corrected genotype-by-ancestry interaction associations for mouse gene expression data. For the yeast, we observed 92 significant genotype-by-ancestry interactions across 38 phenotypes. Given this evidence of epistasis, we test for and observe evidence of rapid selection pressure on ancestry specific polymorphisms within one of the cohorts, consistent with epistatic selection.ConclusionsUnlike our prior work in human populations, we observe widespread evidence of ancestry-modified SNP effects, perhaps reflecting the greater divergence present in crosses using mice and yeast.

In the mammalian heart, the left ventricle (LV) rapidly becomes more dominant in size and function over the right ventricle (RV) after birth. The molecular regulators responsible for this chamber-specific differential growth are largely unknown. We found that cardiomyocytes in the neonatal mouse RV had lower proliferation, more apoptosis, and a smaller average size compared with the LV. This chamber-specific growth pattern was associated with a selective activation of p38 mitogen-activated protein kinase (MAPK) activity in the RV and simultaneous inactivation in the LV. Cardiomyocyte-specific deletion of both the Mapk14 and Mapk11 genes in mice resulted in loss of p38 MAPK expression and activity in the neonatal heart. Inactivation of p38 activity led to a marked increase in cardiomyocyte proliferation and hypertrophy but diminished cardiomyocyte apoptosis, specifically in the RV. Consequently, the p38-inactivated hearts showed RV-specific enlargement postnatally, progressing to pulmonary hypertension and right heart failure at the adult stage. Chamber-specific p38 activity was associated with differential expression of dual-specific phosphatases (DUSPs) in neonatal hearts, including DUSP26. Unbiased transcriptome analysis revealed that IRE1α/XBP1-mediated gene regulation contributed to p38 MAPK-dependent regulation of neonatal cardiomyocyte proliferation and binucleation. These findings establish an obligatory role of DUSP/p38/IRE1α signaling in cardiomyocytes for chamber-specific growth in the postnatal heart.

Cardiovascular diseases are the leading cause of death worldwide. Complex diseases with highly heterogenous disease progression among patient populations, cardiovascular diseases feature multifactorial contributions from both genetic and environmental stressors. Despite significant effort utilizing multiple approaches from molecular biology to genome-wide association studies, the genetic landscape of cardiovascular diseases, particularly for the nonfamilial forms of heart failure, is still poorly understood. In the past decade, systems-level approaches based on omics technologies have become an important approach for the study of complex traits in large populations. These advances create opportunities to integrate genetic variation with other biological layers to identify and prioritize candidate genes, understand pathogenic pathways, and elucidate gene-gene and gene-environment interactions. In this review, we will highlight some of the recent progress made using systems genetics approaches to uncover novel mechanisms and molecular bases of cardiovascular pathophysiological manifestations. The key technology and data analysis platforms necessary to implement systems genetics will be described, and the current major challenges and future directions will also be discussed. For complex cardiovascular diseases, such as heart failure, systems genetics represents a powerful strategy to obtain mechanistic insights and to develop individualized diagnostic and therapeutic regiments, paving the way for precision cardiovascular medicine.

The development of new treatments for heart failure lack animal models that encompass the increasingly heterogeneous disease profile of this patient population. This report provides evidence supporting the hypothesis that Western Diet-fed, aortic-banded Ossabaw swine display an integrated physiological, morphological, and genetic phenotype evocative of cardio-metabolic heart failure. This new preclinical animal model displays a distinctive constellation of findings that are conceivably useful to extending the understanding of how pre-existing cardio-metabolic syndrome can contribute to developing HF.

Right ventricular dysfunction is highly prevalent across cardiopulmonary diseases and independently predicts death in both heart failure (HF) and pulmonary hypertension (PH). Progression towards right ventricular failure (RVF) can occur in spite of optimal medical treatment of HF or PH, highlighting current insufficient understanding of RVF molecular pathophysiology. To identify molecular mechanisms that may distinctly underlie RVF, we investigated the cardiac ventricular transcriptome of advanced HF patients, with and without RVF. Using an integrated systems genomic and functional biology approach, we identified an RVF-specific gene module, for which WIPI1 served as a hub and HSPB6 and MAP4 as drivers, and confirmed the ventricular specificity of Wipi1, Hspb6, and Map4 transcriptional changes in adult murine models of pressure overload induced RV- versus LV- failure. We uncovered a shift towards non-canonical autophagy in the failing RV that correlated with RV-specific Wipi1 upregulation. In vitro siRNA silencing of Wipi1 in neonatal rat ventricular myocytes limited non-canonical autophagy and blunted aldosterone-induced mitochondrial superoxide levels. Our findings suggest that Wipi1 regulates mitochondrial oxidative signaling and non-canonical autophagy in cardiac myocytes. Together with our human transcriptomic analysis and corroborating studies in an RVF mouse model, these data render Wipi1 a potential target for RV-directed HF therapy.

Crouzon syndrome is a congenital craniofacial disorder caused by mutations in the Fibroblast Growth Factor Receptor 2 (FGFR2). It is characterized by the premature fusion of cranial sutures, leading to a brachycephalic head shape, and midfacial hypoplasia. The aim of this study was to investigate the effect of the FGFR2 mutation on the microarchitecture of cranial bones at different stages of postnatal skull development, using the FGFR2C342Y mouse model. Apart from craniosynostosis, this model shows cranial bone abnormalities. High‐resolution synchrotron microtomography images of the frontal and parietal bone were acquired for both FGFR2C342Y/+ (Crouzon, heterozygous mutant) and FGFR2+/+ (control, wild‐type) mice at five ages (postnatal days 1, 3, 7, 14 and 21, n = 6 each). Morphometric measurements were determined for cortical bone porosity: osteocyte lacunae and canals. General linear model to assess the effect of age, anatomical location and genotype was carried out for each morphometric measurement. Histological analysis was performed to validate the findings. In both groups (Crouzon and wild‐type), statistical difference in bone volume fraction, average canal volume, lacunar number density, lacunar volume density and canal volume density was found at most age points, with the frontal bone generally showing higher porosity and fewer lacunae. Frontal bone showed differences between the Crouzon and wild‐type groups in terms of lacunar morphometry (average lacunar volume, lacunar number density and lacunar volume density) with larger, less dense lacunae around the postnatal age of P7–P14. Histological analysis of bone showed marked differences in frontal bone only. These findings provide a better understanding of the pathogenesis of Crouzon syndrome and will contribute to computational models that predict postoperative changes with the aim to improve surgical outcome. [ABSTRACT FROM AUTHOR]

Although the LiCoO2 (LCO) cathode material has been widely used in commercial lithium ion batteries (LIB) and shows high stability, LIB's improvements have several challenges that still need to be overcome. In this paper, we have studied the in-operando structural properties of LCO within battery cells using Bragg Coherent X-ray Diffraction Imaging to identify ways to optimise the LCO batteries' cycling. We have successfully reconstructed the X-ray scattering phase variation (a fingerprint of atomic displacement) within a ≈ (1.6 × 1.4 × 1.3) μm3 LCO nanocrystal across a charge/discharge cycle. Reconstructions indicate strained domains forming, expanding, and fragmenting near the surface of the nanocrystal during charging, with a determined maximum relative lattice displacements of 0.467 Å. While discharging, all domains replicate in reverse the effects observed from the charging states, but with a lower maximum relative lattice displacements of 0.226 Å. These findings show the inefficiency-increasing domain dynamics within LCO lattices during cycling. Although lithium cobalt(III) oxide (LCO) is a widely used and highly stable lithium-ion battery material, several challenges still need to be overcome. Here, the authors use in-operando Bragg Coherent X-ray Diffraction Imaging to visualize the structural properties of LCO in battery cells and identify inefficiency-increasing domain dynamics during cycling. [ABSTRACT FROM AUTHOR]

We describe a simple bioinformatics method for biomarker discovery that is based on the analysis of global transcript levels in a population of inbred mouse strains showing variation for disease-related traits. This method has advantages such as controlled environment and accessibility to heart and plasma tissue in the preclinical selection stage. We illustrate the approach by identifying candidate heart failure (HF) biomarkers by overlaying mouse transcriptome and clinical traits from 91 Hybrid Mouse Diversity Panel (HMDP) inbred strains and human HF transcriptome from the Myocardial Applied Genomics Network (MAGNet) consortium. We found that some of the top differentially expressed genes correlated with known human HF biomarkers, such as galectin-3 and tissue inhibitor of metalloproteinase 1. Using ELISA assays, we investigated one novel candidate, Glycoprotein NMB, in a mouse model of chronic β-adrenergic stimulation by isoproterenol (ISO) induced HF. We observed significantly lower GPNMB plasma levels in the ISO model compared to the control group (p-value = 0.007). In addition, we assessed GPNMB plasma levels among 389 HF cases and controls from the METabolic Syndrome In Men (METSIM) study. Lower levels of GPNMB were also observed in patients with HF from the METSIM study compared to non-HF controls (p-value < 0.0001). In summary, we have identified several candidate biomarkers for HF using the cardiac transcriptome data in a population of mice that may be directly relevant and applicable to human populations.

BackgroundGenetic diversity and the heterogeneous nature of cardiac fibroblasts (CFbs) have hindered characterization of the molecular mechanisms that regulate cardiac fibrosis. The Hybrid Mouse Diversity Panel offers a valuable tool to examine genetically diverse cardiac fibroblasts and their role in fibrosis.MethodsThree strains of mice (C57BL/6J, C3H/HeJ, and KK/HlJ) were selected from the Hybrid Mouse Diversity Panel and treated with either isoproterenol (ISO) or saline by an intraperitoneally implanted osmotic pump. After 21 days, cardiac function and levels of fibrosis were measured by echocardiography and trichrome staining, respectively. Activation and proliferation of CFbs were measured by in vitro and in vivo assays under normal and injury conditions. RNA sequencing was done on isolated CFbs from each strain. Results were analyzed by Ingenuity Pathway Analysis and validated by reverse transcription-qPCR, immunohistochemistry, and ELISA.ResultsISO treatment in C57BL/6J, C3H/HeJ, and KK/HlJ mice resulted in minimal, moderate, and extensive levels of fibrosis, respectively (n=7-8 hearts per condition). Isolated CFbs treated with ISO exhibited strain-specific increases in the levels of activation but showed comparable levels of proliferation. Similar results were found in vivo, with fibroblast activation, and not proliferation, correlating with the differential levels of cardiac fibrosis after ISO treatment. RNA sequencing revealed that CFbs from each strain exhibit unique gene expression changes in response to ISO. We identified Ltbp2 as a commonly upregulated gene after ISO treatment. Expression of LTBP2 was elevated and specifically localized in the fibrotic regions of the myocardium after injury in mice and in human heart failure patients.ConclusionsThis study highlights the importance of genetic variation in cardiac fibrosis by using multiple inbred mouse strains to characterize CFbs and their response to ISO treatment. Our data suggest that, although fibroblast activation is a response that parallels the extent of scar formation, proliferation may not necessarily correlate with levels of fibrosis. In addition, by comparing CFbs from multiple strains, we identified pathways as potential therapeutic targets and LTBP2 as a marker for fibrosis, with relevance to patients with underlying myocardial fibrosis.

A traditional approach to investigate the genetic basis of complex diseases is to identify genes with a global change in expression between diseased and healthy individuals. However, population heterogeneity may undermine the effort to uncover genes with significant but individual contribution to the spectrum of disease phenotypes within a population. Here we investigate individual changes of gene expression when inducing hypertrophy and heart failure in 100 + strains of genetically distinct mice from the Hybrid Mouse Diversity Panel (HMDP). We find that genes whose expression fold-change correlates in a statistically significant way with the severity of the disease are either up or down-regulated across strains, and therefore missed by a traditional population-wide analysis of differential gene expression. Furthermore, those "fold-change" genes are enriched in human cardiac disease genes and form a dense co-regulated module strongly interacting with the cardiac hypertrophic signaling network in the human interactome. We validate our approach by showing that the knockdown of Hes1, predicted as a strong candidate, induces a dramatic reduction of hypertrophy by 80-90% in neonatal rat ventricular myocytes. Our results demonstrate that individualized approaches are crucial to identify genes underlying complex diseases as well as to develop personalized therapies.

BackgroundThe implantation of mechanical circulatory support devices in heart failure patients is associated with a systemic inflammatory response, potentially leading to death from multiple organ dysfunction syndrome. Previous studies point to the involvement of many mechanisms, but an integrative hypothesis does not yet exist. Using time-dependent whole-genome mRNA expression in circulating leukocytes, we constructed a systems-model to improve mechanistic understanding and prediction of adverse outcomes.MethodsWe sampled peripheral blood mononuclear cells from 22 consecutive patients undergoing mechanical circulatory support device (MCS) surgery, at 5 timepoints: day -1 preoperative, and postoperative days 1, 3, 5, and 8. Clinical phenotyping was performed using 12 clinical parameters, 2 organ dysfunction scoring systems, and survival outcomes. We constructed a strictly phenotype-driven time-dependent non-supervised systems-representation using weighted gene co-expression network analysis, and annotated eigengenes using gene ontology, pathway, and transcription factor binding site enrichment analyses. Genes and eigengenes were mapped to the clinical phenotype using a linear mixed-effect model, with Cox models also fit at each timepoint to survival outcomes.ResultsWe inferred a 19-module network, in which most module eigengenes correlated with at least one aspect of the clinical phenotype. We observed a response of advanced heart failure patients to surgery orchestrated into stages: first, activation of the innate immune response, followed by anti-inflammation, and finally reparative processes such as mitosis, coagulation, and apoptosis. Eigengenes related to red blood cell production and extracellular matrix degradation became predictors of survival late in the timecourse corresponding to multiorgan dysfunction and disseminated intravascular coagulation.ConclusionsOur model provides an integrative representation of leukocyte biology during the systemic inflammatory response following MCS device implantation. It demonstrates consistency with previous hypotheses, identifying a number of known mechanisms. At the same time, it suggests novel hypotheses about time-specific targets.

Adult mammalian cardiomyocyte regeneration after injury is thought to be minimal. Mononuclear diploid cardiomyocytes (MNDCMs), a relatively small subpopulation in the adult heart, may account for the observed degree of regeneration, but this has not been tested. We surveyed 120 inbred mouse strains and found that the frequency of adult mononuclear cardiomyocytes was surprisingly variable (>7-fold). Cardiomyocyte proliferation and heart functional recovery after coronary artery ligation both correlated with pre-injury MNDCM content. Using genome-wide association, we identified Tnni3k as one gene that influences variation in this composition and demonstrated that Tnni3k knockout resulted in elevated MNDCM content and increased cardiomyocyte proliferation after injury. Reciprocally, overexpression of Tnni3k in zebrafish promoted cardiomyocyte polyploidization and compromised heart regeneration. Our results corroborate the relevance of MNDCMs in heart regeneration. Moreover, they imply that intrinsic heart regeneration is not limited nor uniform in all individuals, but rather is a variable trait influenced by multiple genes.

Cardiac failure has been widely associated with an increase in glucose utilization. The aim of our study was to identify factors that mechanistically bridge this link between hyperglycemia and heart failure. Here, we screened the Hybrid Mouse Diversity Panel (HMDP) for substrate-specific cardiomyocyte candidates based on heart transcriptional profile and circulating nutrients. Next, we utilized an in vitro model of rat cardiomyocytes to demonstrate that the gene expression changes were in direct response to substrate abundance. After overlaying candidates of interest with a separate HMDP study evaluating isoproterenol-induced heart failure, we chose to focus on the gene Trp53inp2 as a cardiomyocyte glucose utilization-specific factor. Trp53inp2 gene knockdown in rat cardiomyocytes reduced expression and protein abundance of key glycolytic enzymes. This resulted in reduction of both glucose uptake and glycogen content in cardiomyocytes stimulated with isoproterenol. Furthermore, this reduction effectively blunted the capacity of glucose and isoprotereonol to synergistically induce hypertrophic gene expression and cell size expansion. We conclude that Trp53inp2 serves as regulator of cardiomyocyte glycolytic activity and can consequently regulate hypertrophic response in the context of elevated glucose content.NEW & NOTEWORTHY Here, we apply a novel method for screening transcripts based on a substrate-specific expression pattern to identify Trp53inp2 as an induced cardiomyocyte glucose utilization factor. We further show that reducing expression of the gene could effectively blunt hypertrophic response in the context of elevated glucose content.

In this chapter we address the recent explosion in large multilevel population studies such as the METSIM study in humans as well as large panels of animal models such as the Hybrid Mouse Diversity Panel or the BXD set of recombinant inbred strains. These studies have harnessed the increasing affordability of large-scale high-throughput profiling to gather massive quantities of data. These datasets, spread across different -omics levels (genome, transcriptome, etc.), different tissues (e.g. heart, plasma, bone) and different environmental factors (e.g. diet, drugs) each individually have led to a number of novel findings relevant to a variety of complex diseases and other phenotypes. The analysis of these results, however, is often limited to individuals with a comprehensive understanding of database languages such as SQL. In this chapter, we describe the development of a GUI-based database analysis suite, using the Hybrid Mouse Diversity Panel as an example to lay out a series of methods for visualization and integration of large systems genetics datasets. The database is based on the Shiny suite of tools in R, and is transferrable to other SQL-based datasets.

BackgroundChronic myocardial infarction (MI) triggers pathological remodeling in the heart and cardiac nervous system. Abnormal function of the autonomic nervous system (ANS), including stellate ganglia (SG) and dorsal root ganglia (DRG) contribute to increased sympathoexcitation, cardiac dysfunction and arrythmogenesis. ANS modulation is a therapeutic target for arrhythmia associated with cardiac injury. However, the molecular mechanism involved in the pathological remodeling in ANS following cardiac injury remains to be established.Methods and resultsIn this study, we performed transcriptome analysis by RNA-sequencing in thoracic SG and (T1-T4) DRG obtained from Yorkshire pigs following either acute (3 to 5 hours) or chronic (8 weeks) myocardial infarction. By differential expression and weighted gene co-expression network analysis (WGCNA), we identified significant transcriptome changes and specific gene modules in the ANS tissues in response to myocardial infarction at either acute or chronic phases. Both differential expressed genes and the member genes of the WGCNA gene module associated with post-infarct condition were significantly enriched for inflammatory signaling and apoptotic cell death. Targeted validation analysis supported a significant induction of inflammatory and apoptotic signal in both SG and DRG following myocardial infarction, along with cellular evidence of apoptosis induction based on TUNEL analysis. Importantly, these molecular changes were observed specifically in the thoracic segments but not in their counterparts obtained from lumbar sections.ConclusionMyocardial injury leads to time-dependent global changes in gene expression in the innervating ANS. Induction of inflammatory gene expression and loss of neuron cell viability in SG and DRG are potential novel mechanisms contributing to abnormal ANS function which can promote cardiac arrhythmia and pathological remodeling in myocardium.

We previously reported a genetic analysis of heart failure traits in a population of inbred mouse strains treated with isoproterenol to mimic catecholamine-driven cardiac hypertrophy. Here, we apply a co-expression network algorithm, wMICA, to perform a systems-level analysis of left ventricular transcriptomes from these mice. We describe the features of the overall network but focus on a module identified in treated hearts that is strongly related to cardiac hypertrophy and pathological remodeling. Using the causal modeling algorithm NEO, we identified the gene Adamts2 as a putative regulator of this module and validated the predictive value of NEO using small interfering RNA-mediated knockdown in neonatal rat ventricular myocytes. Adamts2 silencing regulated the expression of the genes residing within the module and impaired isoproterenol-induced cellular hypertrophy. Our results provide a view of higher order interactions in heart failure with potential for diagnostic and therapeutic insights.

A systems approach deconvolutes genes specific to and enriched in endothelium from whole-organ transcriptome data, with applications to other cell types and tissues.

We aimed to understand the genetic control of cardiac remodeling using an isoproterenol-induced heart failure model in mice, which allowed control of confounding factors in an experimental setting. We characterized the changes in cardiac structure and function in response to chronic isoproterenol infusion using echocardiography in a panel of 104 inbred mouse strains. We showed that cardiac structure and function, whether under normal or stress conditions, has a strong genetic component, with heritability estimates of left ventricular mass between 61% and 81%. Association analyses of cardiac remodeling traits, corrected for population structure, body size and heart rate, revealed 17 genome-wide significant loci, including several loci containing previously implicated genes. Cardiac tissue gene expression profiling, expression quantitative trait loci, expression-phenotype correlation, and coding sequence variation analyses were performed to prioritize candidate genes and to generate hypotheses for downstream mechanistic studies. Using this approach, we have validated a novel gene, Myh14, as a negative regulator of ISO-induced left ventricular mass hypertrophy in an in vivo mouse model and demonstrated the up-regulation of immediate early gene Myc, fetal gene Nppb, and fibrosis gene Lgals3 in ISO-treated Myh14 deficient hearts compared to controls.

Transcriptome remodeling in heart disease occurs through the coordinated actions of transcription factors, histone modifications, and other chromatin features at pathology-associated genes. The extent to which genome-wide chromatin reorganization also contributes to the resultant changes in gene expression remains unknown. We examined the roles of two chromatin structural proteins, Ctcf (CCCTC-binding factor) and Hmgb2 (high mobility group protein B2), in regulating pathologic transcription and chromatin remodeling. Our data demonstrate a reciprocal relationship between Hmgb2 and Ctcf in controlling aspects of chromatin structure and gene expression. Both proteins regulate each others' expression as well as transcription in cardiac myocytes; however, only Hmgb2 does so in a manner that involves global reprogramming of chromatin accessibility. We demonstrate that the actions of Hmgb2 on local chromatin accessibility are conserved across genomic loci, whereas the effects on transcription are loci-dependent and emerge in concert with histone modification and other chromatin features. Finally, although both proteins share gene targets, Hmgb2 and Ctcf, neither binds these genes simultaneously nor do they physically colocalize in myocyte nuclei. Our study uncovers a previously unknown relationship between these two ubiquitous chromatin proteins and provides a mechanistic explanation for how Hmgb2 regulates gene expression and cellular phenotype. Furthermore, we provide direct evidence for structural remodeling of chromatin on a genome-wide scale in the setting of cardiac disease.

The Hybrid Mouse Diversity Panel (HMDP) is a collection of approximately 100 well-characterized inbred strains of mice that can be used to analyze the genetic and environmental factors underlying complex traits. While not nearly as powerful for mapping genetic loci contributing to the traits as human genome-wide association studies, it has some important advantages. First, environmental factors can be controlled. Second, relevant tissues are accessible for global molecular phenotyping. Finally, because inbred strains are renewable, results from separate studies can be integrated. Thus far, the HMDP has been studied for traits relevant to obesity, diabetes, atherosclerosis, osteoporosis, heart failure, immune regulation, fatty liver disease, and host-gut microbiota interactions. High-throughput technologies have been used to examine the genomes, epigenomes, transcriptomes, proteomes, metabolomes, and microbiomes of the mice under various environmental conditions. All of the published data are available and can be readily used to formulate hypotheses about genes, pathways and interactions.

BackgroundAlthough metabolic reprogramming is critical in the pathogenesis of heart failure, studies to date have focused principally on fatty acid and glucose metabolism. Contribution of amino acid metabolic regulation in the disease remains understudied.Methods and resultsTranscriptomic and metabolomic analyses were performed in mouse failing heart induced by pressure overload. Suppression of branched-chain amino acid (BCAA) catabolic gene expression along with concomitant tissue accumulation of branched-chain α-keto acids was identified as a significant signature of metabolic reprogramming in mouse failing hearts and validated to be shared in human cardiomyopathy hearts. Molecular and genetic evidence identified the transcription factor Krüppel-like factor 15 as a key upstream regulator of the BCAA catabolic regulation in the heart. Studies using a genetic mouse model revealed that BCAA catabolic defect promoted heart failure associated with induced oxidative stress and metabolic disturbance in response to mechanical overload. Mechanistically, elevated branched-chain α-keto acids directly suppressed respiration and induced superoxide production in isolated mitochondria. Finally, pharmacological enhancement of branched-chain α-keto acid dehydrogenase activity significantly blunted cardiac dysfunction after pressure overload.ConclusionsBCAA catabolic defect is a metabolic hallmark of failing heart resulting from Krüppel-like factor 15-mediated transcriptional reprogramming. BCAA catabolic defect imposes a previously unappreciated significant contribution to heart failure.

RationaleOnly a small portion of the known heritability of cardiovascular diseases, such as heart failure, can be explained based on single-gene mutations. Chromatin structure and regulation provide a substrate through which genetic differences in noncoding regions may affect cellular function and response to disease, but the mechanisms are unknown.ObjectiveWe conducted genome-wide measurements of DNA methylation in different strains of mice that are susceptible and resistant to isoproterenol-induced dysfunction to test the hypothesis that this epigenetic mark may play a causal role in the development of heart failure.Methods and resultsBALB/cJ and BUB/BnJ mice, determined to be susceptible and resistant to isoproterenol-induced heart failure, respectively, were administered the drug for 3 weeks via osmotic minipump. Reduced representational bisulfite sequencing was then used to compare the differences between the cardiac DNA methylomes in the basal state between strains and then after isoproterenol treatment. Single-base resolution DNA methylation measurements were obtained and revealed a bimodal distribution of methylation in the heart, enriched in lone intergenic CpGs and depleted from CpG islands around genes. Isoproterenol induced global decreases in methylation in both strains; however, the basal methylation pattern between strains shows striking differences that may be predictive of disease progression before environmental stress. The global correlation between promoter methylation and gene expression (as measured by microarray) was modest and revealed itself only with focused analyses of transcription start site and gene body regions (in contrast to when gene methylation was examined in toto). Modules of comethylated genes displayed correlation with other protein-based epigenetic marks, supporting the hypothesis that chromatin modifications act in a combinatorial manner to specify transcriptional phenotypes in the heart.ConclusionsThis study provides the first single-base resolution map of the mammalian cardiac DNA methylome and the first case-control analysis of the changes in DNA methylation with heart failure. The findings demonstrate marked genetic differences in DNA methylation that are associated with disease progression.

RNA splicing is a major contributor to total transcriptome complexity; however, the functional role and regulation of splicing in heart failure remain poorly understood. Here, we used a total transcriptome profiling and bioinformatic analysis approach and identified a muscle-specific isoform of an RNA splicing regulator, RBFox1 (also known as A2BP1), as a prominent regulator of alternative RNA splicing during heart failure. Evaluation of developing murine and zebrafish hearts revealed that RBFox1 is induced during postnatal cardiac maturation. However, we found that RBFox1 is markedly diminished in failing human and mouse hearts. In a mouse model, RBFox1 deficiency in the heart promoted pressure overload-induced heart failure. We determined that RBFox1 is a potent regulator of RNA splicing and is required for a conserved splicing process of transcription factor MEF2 family members that yields different MEF2 isoforms with differential effects on cardiac hypertrophic gene expression. Finally, induction of RBFox1 expression in murine pressure overload models substantially attenuated cardiac hypertrophy and pathological manifestations. Together, this study identifies regulation of RNA splicing by RBFox1 as an important player in transcriptome reprogramming during heart failure that influence pathogenesis of the disease.