Search

Background: Heart 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., Method: We 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-Lep ob /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., Results: In contrast to the post-transverse aortic constriction C57Bl/6J lean mice, which developed pathological features of heart failure with reduced ejection fraction over time, the post-transverse 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., Conclusions: These 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., Competing Interests: D.J. Drucker has served as a consultant or speaker within the past 12 months to Altimmune, Amgen, Boehringer Ingelheim, Kallyope, Merck Research Laboratories, Novo Nordisk Inc, and Pfizer Inc. Neither D.J. Drucker or his family members hold issued stock directly or indirectly in any of these companies. He holds nonexercised options in Kallyope. The other authors report no conflicts.

Background: Recent advances in single cell sequencing have led to an increased focus on the role of cell-type composition in phenotypic presentation and disease progression. Cell-type composition research in the heart is challenging due to large, frequently multinucleated cardiomyocytes that preclude most single cell approaches from obtaining accurate measurements of cell composition. Our in silico studies reveal that ignoring cell type composition when calculating differentially expressed genes (DEGs) can have significant consequences. For example, a relatively small change in cell abundance of only 10% can result in over 25% of DEGs being false positives., Methods: We have implemented an algorithmic approach that uses snRNAseq datasets as a reference to accurately calculate cell type compositions from bulk RNAseq datasets through robust data cleaning, gene selection, and multi-sample cross-subject and cross-cell-type deconvolution. We applied our approach to cardiomyocyte-specific α1A adrenergic receptor (CM-α1A-AR) knockout mice. 8-12 week-old mice (either WT or CM-α1A-KO) were subjected to permanent left coronary artery (LCA) ligation or sham surgery (n=4 per group). Transcriptomes from the infarct border zones were collected 3 days later and analyzed using our algorithm to determine cell-type abundances, corrected differential expression calculations using DESeq2, and validated these findings using RNAscope., Results: Uncorrected DEGs for the CM-α1A-KO X LCA interaction term featured many cell-type specific genes such as Timp4 (fibroblasts) and Aplnr (cardiomyocytes) and overall GO enrichment for terms pertaining to cardiomyocyte differentiation (P=3.1E-4). Using our algorithm, we observe a striking loss of cardiomyocytes and gain in fibroblasts in the α1A-KO + LCA mice that was not recapitulated in WT + LCA animals, although we did observe a similar increase in macrophage abundance in both conditions. This recapitulates prior results that showed a much more severe heart failure phenotype in CM-α1A-KO + LCA mice. Following correction for cell-type, our DEGs now highlight a novel set of genes enriched for GO terms such as cardiac contraction (P=3.7E-5) and actin filament organization (P=6.3E-5)., Conclusions: Our algorithm identifies and corrects for cell-type abundance in bulk RNAseq datasets opening new avenues for research on novel genes and pathways as well as an improved understanding of the role of cardiac cell types in cardiovascular disease.

Tuning of genome structure and function is accomplished by chromatin-binding proteins, which determine the transcriptome and phenotype of the cell. Here we investigate how communication between extracellular stress and chromatin structure may regulate cellular mechanical behaviors. We demonstrate that histone H1.0, which compacts nucleosomes into higher-order chromatin fibers, controls genome organization and cellular stress response. We show that histone H1.0 has privileged expression in fibroblasts across tissue types and that its expression is necessary and sufficient to induce myofibroblast activation. Depletion of histone H1.0 prevents cytokine-induced fibroblast contraction, proliferation and migration via inhibition of a transcriptome comprising extracellular matrix, cytoskeletal and contractile genes, through a process that involves locus-specific H3K27 acetylation. Transient depletion of histone H1.0 in vivo prevents fibrosis in cardiac muscle. These findings identify an unexpected role of linker histones to orchestrate cellular mechanical behaviors, directly coupling force generation, nuclear organization and gene transcription., Competing Interests: Competing interestsThe authors declare no competing interests., (© The Author(s) 2024.)

Cardiometabolic disorders encompass a broad range of cardiovascular complications associated with metabolic dysfunction. These conditions have an increasing share in the health burden worldwide due to worsening endemic of hypertension, obesity, and diabetes. Previous studies have identified Tumor Protein p53-inducible Nuclear Protein 2 (Trp53inp2) as a molecular link between hyperglycemia and cardiac hypertrophy. However, its role in cardiac pathology has never been determined in vivo . In this study, we generated a cardiac specific knockout model of Trp53inp2 (Trp53inp2-cKO) and investigated the impact of Trp53inp2 inactivation on the pathogenesis of heart failure under mechanic or/and metabolic stresses. Based on echocardiography assessment, inactivation of Trp53inp2 in heart led to accelerated onset of HFrEF in response to pressure-overload, with significantly reduced ejection fraction and elevated heart failure marker genes comparing to the control mice. In contrast, inactivation of Trp53inp2 ameliorated cardiac dysfunction induced by combined stresses of high fat diet and moderate pressure overload (Cardiometabolic Disorder Model). Moreover, Trp53inp2 inactivation led to reduced expression of glucose metabolism genes in lean, pressure-overloaded hearts. However, the same set of genes were significantly induced in the Trp53inp2-cKO hearts under both mechanical and metabolic stresses. In summary, we have demonstrated for the first time that cardiomyocyte Trp53inp2 has diametrically differential roles in the pathogenesis of heart failure and glucose regulation under mechanical vs. mechanical plus metabolic stresses. This insight suggests that Trp53inp2 may exacerbate the cardiac dysfunction during pressure overload injury but have a protective effect in cardiac diastolic function in cardiometabolic disease., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (© 2023 Liu, Liu, Ren, Zhu, Bouso, Mamlouk, Rau, Wang and Gao.)

Competing Interests: Disclosures None.

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., (© 2023. The Author(s).)

Somatic polyploidization, an adaptation by which cells increase their DNA content to support growth, is observed in many cell types, including cardiomyocytes. Although polyploidization is believed to be beneficial, progression to a polyploid state is often accompanied by loss of proliferative capacity. Recent work suggests that genetics heavily influence cardiomyocyte ploidy. However, the developmental course by which cardiomyocytes reach their final ploidy state has only been investigated in select backgrounds. Here, we assessed cardiomyocyte number, cell cycle activity, and ploidy dynamics across two divergent mouse strains: C57BL/6J and A/J. Both strains are born and reach adulthood with comparable numbers of cardiomyocytes; however, the end composition of ploidy classes and developmental progression to reach the final state differ substantially. We expand on previous findings that identified Tnni3k as a mediator of cardiomyocyte ploidy and uncover a role for Runx1 in ploidy dynamics and cardiomyocyte cell division, in both developmental and injury contexts. These data provide novel insights into the developmental path to cardiomyocyte polyploidization and challenge the paradigm that hypertrophy is the sole mechanism for growth in the postnatal heart., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2023. Published by The Company of Biologists Ltd.)

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., (© 2022. The Author(s).)

Background: Human-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., Objective: To improve the electrophysiological maturation of hiPSC-CMs., Methods and Results: On 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 (I f ) and increased inward rectifier potassium currents (I K1 ), sodium currents (I Na ), and the rapidly and slowly activating delayed rectifier potassium currents (I Kr and I Ks , 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 I K1 by 100 μM BaCl 2 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 I K1 for conduction velocity and action potential duration., Conclusion: T3 + Dex treatment is an effective approach to rapidly enhance electrophysiological maturation of hiPSC-CMs., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

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.

Heart failure (HF) patients with deteriorating right ventricular (RV) structure and function have a nearly twofold increased risk of death compared with those without. Despite the well-established clinical risk, few studies have examined the molecular signature associated with this HF condition. The purpose of this study was to integrate morphological, molecular, and functional data with the transcriptome data set in the RV of a preclinical model of cardiometabolic HF. Ossabaw swine were fed either normal diet without surgery (lean control, n = 5) or Western diet and aortic-banding (WD-AB; n = 4). Postmortem RV weight was increased and positively correlated with lung weight in the WD-AB group compared with CON. Total RNA-seq was performed and gene expression profiles were compared and analyzed using principal component analysis, weighted gene co-expression network analysis, module enrichment analysis, and ingenuity pathway analysis. Gene networks specifically associated with RV hypertrophic remodeling identified a hub gene in MAPK8 (or JNK1) that was associated with the selective induction of the extracellular matrix (ECM) component fibronectin. JNK1 and fibronectin protein were increased in the right coronary artery (RCA) of WD-AB animals and associated with a decrease in matrix metalloproteinase 14 protein, which specifically degrades fibronectin. RCA fibronectin content was correlated with increased vascular stiffness evident as a decreased elastin elastic modulus in WD-AB animals. In conclusion, this study establishes a molecular and transcriptome signature in the RV using Ossabaw swine with cardiometabolic HF. This signature was associated with altered ECM regulation and increased vascular stiffness in the RCA, with selective dysregulation of fibronectin.

Background: The 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., Results: We 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., Conclusions: Unlike 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., Competing Interests: The authors have declared that no competing interests exist.

Matrix metalloproteinases (MMPs) and tissue inhibitors of metalloproteinases (TIMPs) are pivotal regulators of extracellular matrix (ECM) composition and could, due to their dynamic activity, function as prognostic tools for fibrosis and cardiac function in left ventricular diastolic dysfunction (LVDD) and heart failure with preserved ejection fraction (HFpEF). We conducted a systematic review on experimental animal models of LVDD and HFpEF published in MEDLINE or Embase. Twenty-three studies were included with a total of 36 comparisons that reported established LVDD, quantification of cardiac fibrosis and cardiac MMP or TIMP expression or activity. LVDD/HFpEF models were divided based on underlying pathology: hemodynamic overload (17 comparisons), metabolic alteration (16 comparisons) or ageing (3 comparisons). Meta-analysis showed that echocardiographic parameters were not consistently altered in LVDD/HFpEF with invasive hemodynamic measurements better representing LVDD. Increased myocardial fibrotic area indicated comparable characteristics between hemodynamic and metabolic models. Regarding MMPs and TIMPs; MMP2 and MMP9 activity and protein and TIMP1 protein levels were mainly enhanced in hemodynamic models. In most cases only mRNA was assessed and there were no correlations between cardiac tissue and plasma levels. Female gender, a known risk factor for LVDD and HFpEF, was underrepresented. Novel studies should detail relevant model characteristics and focus on MMP and TIMP protein expression and activity to identify predictive circulating markers in cardiac ECM remodeling.

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.

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.

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., Competing Interests: The authors declare no competing financial interests.

Isoproterenol is used widely for inducing heart failure in mice. Isoproterenol is a nonselective beta-adrenergic agonist. The acute model mimics stress-induced cardiomyopathy. The chronic model mimics advanced heart failure in humans. In this chapter, we describe a protocol that we used to induce heart failure in 100+ strains of inbred mice. Techniques on surgical pump implantation and echocardiography are described in detail. We also discuss the impact of drug dosage, duration, mortality, age, gender, and strain on cardiac remodeling responses. The success of model creation may be assessed by echocardiogram or molecular markers. This chapter may be relevant to those who are interested in using this heart failure model.

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., (Copyright © 2017 the American Physiological Society.)

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., (Copyright © 2017 The Author(s). Published by Elsevier Inc. All rights reserved.)

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.

A systems approach deconvolutes genes specific to and enriched in endothelium from whole-organ transcriptome data, with applications to other cell types and tissues., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Expression of a cohort of disease-associated genes, some of which are active in fetal myocardium, is considered a hallmark of transcriptional change in cardiac hypertrophy models. How this transcriptome remodeling is affected by the common genetic variation present in populations is unknown. We examined the role of genetics, as well as contributions of chromatin proteins, to regulate cardiac gene expression and heart failure susceptibility. We examined gene expression in 84 genetically distinct inbred strains of control and isoproterenol-treated mice, which exhibited varying degrees of disease. Unexpectedly, fetal gene expression was not correlated with hypertrophic phenotypes. Unbiased modeling identified 74 predictors of heart mass after isoproterenol-induced stress, but these predictors did not enrich for any cardiac pathways. However, expanded analysis of fetal genes and chromatin remodelers as groups correlated significantly with individual systemic phenotypes. Yet, cardiac transcription factors and genes shown by gain-/loss-of-function studies to contribute to hypertrophic signaling did not correlate with cardiac mass or function in disease. Because the relationship between gene expression and phenotype was strain specific, we examined genetic contribution to expression. Strikingly, strains with similar transcriptomes in the basal heart did not cluster together in the isoproterenol state, providing comprehensive evidence that there are different genetic contributors to physiological and pathological gene expression. Furthermore, the divergence in transcriptome similarity versus genetic similarity between strains is organ specific and genome-wide, suggesting chromatin is a critical buffer between genetics and gene expression., (Copyright © 2016 the American Physiological Society.)

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., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

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., (Copyright © 2016 by the American Society for Biochemistry and Molecular Biology, Inc.)

Background: Although 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 Results: Transcriptomic 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., Conclusions: BCAA 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., (© 2016 American Heart Association, Inc.)

Rationale: Only 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., Objective: We 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 Results: BALB/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., Conclusions: This 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., (© 2016 American Heart Association, Inc.)

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.

Human genome-wide association studies have identified thousands of loci associated with disease phenotypes. Genome-wide association studies also have become feasible using rodent models and these have some important advantages over human studies, including controlled environment, access to tissues for molecular profiling, reproducible genotypes, and a wide array of techniques for experimental validation. Association mapping with common mouse inbred strains generally requires 100 or more strains to achieve sufficient power and mapping resolution; in contrast, sample sizes for human studies typically are one or more orders of magnitude greater than this. To enable well-powered studies in mice, we have generated high-density genotypes for ∼175 inbred strains of mice using the Mouse Diversity Array. These new data increase marker density by 1.9-fold, have reduced missing data rates, and provide more accurate identification of heterozygous regions compared with previous genotype data. We report the discovery of new loci from previously reported association mapping studies using the new genotype data. The data are freely available for download, and Web-based tools provide easy access for association mapping and viewing of the underlying intensity data for individual loci., (Copyright © 2015 Rau et al.)

Purpose of Review: In contrast to many other human diseases, the use of genome-wide association studies (GWAS) to identify genes for heart failure (HF) has had limited success. We will discuss the underlying challenges as well as potential new approaches to understanding the genetics of common forms of HF., Recent Findings: Recent research using intermediate phenotypes, more detailed and quantitative stratification of HF symptoms, founder populations and novel animal models has begun to allow researchers to make headway toward explaining the genetics underlying HF using GWAS techniques., Summary: By expanding analyses of HF to improved clinical traits, additional HF classifications and innovative model systems, the intractability of human HF GWAS should be ameliorated significantly.

Insulin resistance (IR) is a complex trait with multiple genetic and environmental components. Confounded by large differences between the sexes, environment, and disease pathology, the genetic basis of IR has been difficult to dissect. Here we examine IR and related traits in a diverse population of more than 100 unique male and female inbred mouse strains after feeding a diet rich in fat and refined carbohydrates. Our results show dramatic variation in IR among strains of mice and widespread differences between sexes that are dependent on genotype. We uncover more than 15 genome-wide significant loci and validate a gene, Agpat5, associated with IR. We also integrate plasma metabolite levels and global gene expression from liver and adipose tissue to identify metabolite quantitative trait loci (mQTL) and expression QTL (eQTL), respectively. Our results provide a resource for analysis of interactions between diet, sex, and genetic background in IR., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Background: Chronic stress-induced cardiac pathology exhibits both a wide range in severity and a high degree of heterogeneity in clinical manifestation in human patients. This variability is contributed to by complex genetic and environmental etiologies within the human population. Genetic approaches to elucidate the genetics underlying the acquired forms of cardiomyopathies, including genome-wide association studies, have been largely unsuccessful, resulting in limited knowledge as to the contribution of genetic variations for this important disease., Methods and Results: Using the β-adrenergic agonist isoproterenol as a specific pathological stressor to circumvent the problem of etiologic heterogeneity, we performed a genome-wide association study for genes influencing cardiac hypertrophy and fibrosis in a large panel of inbred mice. Our analyses revealed 7 significant loci and 17 suggestive loci, containing an average of 14 genes, affecting cardiac hypertrophy, fibrosis, and surrogate traits relevant to heart failure. Several loci contained candidate genes which are known to contribute to Mendelian cardiomyopathies in humans or have established roles in cardiac pathology based on molecular or genetic studies in mouse models. In particular, we identify Abcc6 as a novel gene underlying a fibrosis locus by validating that an allele with a splice mutation of Abcc6 dramatically and rapidly promotes isoproterenol-induced cardiac fibrosis., Conclusions: Genetic variants significantly contribute to the phenotypic heterogeneity of stress-induced cardiomyopathy. Systems genetics is an effective approach to identify genes and pathways underlying the specific pathological features of cardiomyopathies. Abcc6 is a previously unrecognized player in the development of stress-induced cardiac fibrosis., (© 2014 American Heart Association, Inc.)

The protein phosphatase 1-like gene (PPM1l) was identified as causal gene for obesity and metabolic abnormalities in mice. However, the underlying mechanisms were unknown. In this report, we find PPM1l encodes an endoplasmic reticulum (ER) membrane targeted protein phosphatase (PP2Ce) and has specific activity to basal and ER stress induced auto-phosphorylation of Inositol-REquiring protein-1 (IRE1). PP2Ce inactivation resulted in elevated IRE1 phosphorylation and higher expression of XBP-1, CHOP, and BiP at basal. However, ER stress stimulated XBP-1 and BiP induction was blunted while CHOP induction was further enhanced in PP2Ce null cells. PP2Ce protein levels are significantly induced during adipogenesis in vitro and are necessary for normal adipocyte maturation. Finally, we provide evidence that common genetic variation of PPM11 gene is significantly associated with human lipid profile. Therefore, PPM1l mediated IRE1 regulation and downstream ER stress signaling is a plausible molecular basis for its role in metabolic regulation and disorder.

Genetic variations in blood cell parameters can impact clinical traits. We report here the mapping of blood cell traits in a panel of 100 inbred strains of mice of the Hybrid Mouse Diversity Panel (HMDP) using genome-wide association (GWA). We replicated a locus previously identified in using linkage analysis in several genetic crosses for mean corpuscular volume (MCV) and a number of other red blood cell traits on distal chromosome 7. Our peak for SNP association to MCV occurred in a linkage disequilibrium (LD) block spanning from 109.38 to 111.75 Mb that includes Hbb-b1, the likely causal gene. Altogether, we identified five loci controlling red blood cell traits (on chromosomes 1, 7, 11, 12, and 16), and four of these correspond to loci for red blood cell traits reported in a recent human GWA study. For white blood cells, including granulocytes, monocytes, and lymphocytes, a total of six significant loci were identified on chromosomes 1, 6, 8, 11, 12, and 15. An average of ten candidate genes were found at each locus and those were prioritized by examining functional variants in the HMDP such as missense and expression variants. These results provide intermediate phenotypes and candidate loci for genetic studies of atherosclerosis and cancer as well as inflammatory and immune disorders in mice.

Background: Network construction and analysis algorithms provide scientists with the ability to sift through high-throughput biological outputs, such as transcription microarrays, for small groups of genes (modules) that are relevant for further research. Most of these algorithms ignore the important role of non-linear interactions in the data, and the ability for genes to operate in multiple functional groups at once, despite clear evidence for both of these phenomena in observed biological systems., Results: We have created a novel co-expression network analysis algorithm that incorporates both of these principles by combining the information-theoretic association measure of the maximal information coefficient (MIC) with an Interaction Component Model. We evaluate the performance of this approach on two datasets collected from a large panel of mice, one from macrophages and the other from liver by comparing the two measures based on a measure of module entropy, Gene Ontology (GO) enrichment, and scale-free topology (SFT) fit. Our algorithm outperforms a widely used co-expression analysis method, weighted gene co-expression network analysis (WGCNA), in the macrophage data, while returning comparable results in the liver dataset when using these criteria. We demonstrate that the macrophage data has more non-linear interactions than the liver dataset, which may explain the increased performance of our method, termed Maximal Information Component Analysis (MICA) in that case., Conclusions: In making our network algorithm more accurately reflect known biological principles, we are able to generate modules with improved relevance, particularly in networks with confounding factors such as gene by environment interactions.

Obesity is a highly heritable disease driven by complex interactions between genetic and environmental factors. Human genome-wide association studies (GWAS) have identified a number of loci contributing to obesity; however, a major limitation of these studies is the inability to assess environmental interactions common to obesity. Using a systems genetics approach, we measured obesity traits, global gene expression, and gut microbiota composition in response to a high-fat/high-sucrose (HF/HS) diet of more than 100 inbred strains of mice. Here we show that HF/HS feeding promotes robust, strain-specific changes in obesity that are not accounted for by food intake and provide evidence for a genetically determined set point for obesity. GWAS analysis identified 11 genome-wide significant loci associated with obesity traits, several of which overlap with loci identified in human studies. We also show strong relationships between genotype and gut microbiota plasticity during HF/HS feeding and identify gut microbial phylotypes associated with obesity., (Copyright © 2013 Elsevier Inc. All rights reserved.)

We have developed an association-based approach using classical inbred strains of mice in which we correct for population structure, which is very extensive in mice, using an efficient mixed-model algorithm. Our approach includes inbred parental strains as well as recombinant inbred strains in order to capture loci with effect sizes typical of complex traits in mice (in the range of 5% of total trait variance). Over the last few years, we have typed the hybrid mouse diversity panel (HMDP) strains for a variety of clinical traits as well as intermediate phenotypes and have shown that the HMDP has sufficient power to map genes for highly complex traits with resolution that is in most cases less than a megabase. In this essay, we review our experience with the HMDP, describe various ongoing projects, and discuss how the HMDP may fit into the larger picture of common diseases and different approaches.

In this Emerging Science Review, we discuss a systems genetics strategy, which we call gene module association study (GMAS), as a novel approach complementing genome-wide association studies (GWAS), to understand complex diseases by focusing on how genes work together in groups rather than singly. The first step is to characterize phenotypic differences among a genetically diverse population. The second step is to use gene expression microarray (or other high-throughput) data from the population to construct gene coexpression networks. Coexpression analysis typically groups 20 000 genes into 20 to 30 modules containing tens to hundreds of genes, whose aggregate behavior can be represented by the module's "eigengene." The third step is to correlate expression patterns with phenotype, as in GWAS, only applied to eigengenes instead of single nucleotide polymorphisms. The goal of the GMAS approach is to identify groups of coregulated genes that explain complex traits from a systems perspective. From an evolutionary standpoint, we hypothesize that variability in eigengene patterns reflects the "good enough solution" concept, that biological systems are sufficiently complex so that many possible combinations of the same elements (in this case eigengenes) can produce an equivalent output, that is, a "good enough solution" to accomplish normal biological functions. However, when faced with environmental stresses, some "good enough solutions" adapt better than others, explaining individual variability to disease and drug susceptibility. If validated, GMAS may imply that common polygenic diseases are related as much to group interactions between normal genes, as to multiple gene mutations.