Multicellular transcriptional analysis of mammalian heart regeneration

BACKGROUND: The inability of the adult mammalian heart to regenerate following injury represents a major barrier in cardiovascular medicine. In contrast, the neonatal mammalian heart retains a transient capacity for regeneration, which is lost shortly after birth. Defining the molecular mechanisms that govern regenerative capacity in the neonatal period remains a central goal in cardiac biology. Here, we assemble a transcriptomic framework of multiple cardiac cell populations during postnatal development and following injury, which enables comparative analyses of the regenerative (neonatal) versus nonregenerative (adult) state for the first time. METHODS: Cardiomyocytes, fibroblasts, leukocytes, and endothelial cells from infarcted and noninfarcted neonatal (P1) and adult (P56) mouse hearts were isolated by enzymatic dissociation and fluorescence-activated cell sorting at day 3 following surgery. RNA sequencing was performed on these cell populations to generate the transcriptome of the major cardiac cell populations during cardiac development, repair, and regeneration. To complement our transcriptomic data, we also surveyed the epigenetic landscape of cardiomyocytes during postnatal maturation by performing deep sequencing of accessible chromatin regions by using the Assay for Transposase-Accessible Chromatin from purified mouse cardiomyocyte nuclei (P1, P14, and P56). RESULTS: Profiling of cardiomyocyte and nonmyocyte transcriptional programs uncovered several injury-responsive genes across regenerative and nonregenerative time points. However, the majority of transcriptional changes in all cardiac cell types resulted from developmental maturation from neonatal stages to adulthood rather than activation of a distinct regeneration-specific gene program. Furthermore, adult leukocytes and fibroblasts were characterized by the expression of a proliferative gene expression network following infarction, which mirrored the neonatal state. In contrast, cardiomyocytes failed to reactivate the neonatal proliferative network following infarction, which was associated with loss of chromatin accessibility around cell cycle genes during postnatal maturation. CONCLUSIONS: This work provides a comprehensive framework and transcriptional resource of multiple cardiac cell populations during cardiac development, repair, and regeneration. Our findings define a regulatory program underpinning the neonatal regenerative state and identify alterations in the chromatin landscape that could limit reinduction of the regenerative program in adult cardiomyocytes

Vascular histone deacetylation by pharmacological HDAC inhibition

HDAC inhibitors can regulate gene expression by post-translational modification of histone as well as nonhistone proteins. Often studied at single loci, increased histone acetylation is the paradigmatic mechanism of action. However, little is known of the extent of genome-wide changes in cells stimulated by the hydroxamic acids, TSA and SAHA. In this article, we map vascular chromatin modifications including histone H3 acetylation of lysine 9 and 14 (H3K9/14ac) using chromatin immunoprecipitation (ChIP) coupled with massive parallel sequencing (ChIP-seq). Since acetylation-mediated gene expression is often associated with modification of other lysine residues, we also examined H3K4me3 and H3K9me3 as well as changes in CpG methylation (CpG-seq). RNA sequencing indicates the differential expression of ∼30% of genes, with almost equal numbers being up- and down-regulated. We observed broad deacetylation and gene expression changes conferred by TSA and SAHA mediated by the loss of EP300/CREBBP binding at multiple gene promoters. This study provides an important framework for HDAC inhibitor function in vascular biology and a comprehensive description of genome-wide deacetylation by pharmacological HDAC inhibition

RNA sequencing supports distinct reactive oxygen species-mediated pathways of apoptosis by high and low size mass fractions of Bay leaf (Lauris nobilis) in HT-29 cells

Anti-proliferative and pro-apoptotic effects of Bay leaf (Laurus nobilis) in mammalian cancer and HT-29 adenocarcinoma cells have been previously attributed to effects of polyphenolic and essential oil chemical species. Recently, we demonstrated differentiated growth-regulating effects of high (HFBL) versus low molecular mass (LFBL) aqueous fractions of bay leaf and now confirm by comparative effects on gene expression, that HFBL and LFBL suppress HT-29 growth by distinct mechanisms. Induction of intra-cellular lesions including DNA strand breakage by extra-cellular HFBL, invoked the hypothesis that iron-mediated reactive oxygen species with capacity to penetrate cell membrane, were responsible for HFBL-mediated effects, supported by equivalent effects of HFBL in combination with γ radiation. Activities of HFBL and LFBL were interpreted to reflect differentiated responses to iron-mediated reactive oxygen species (ROS), occurring either outside or inside cells. In the presence of LFBL, apoptotic death was relatively delayed compared with HFBL. ROS production by LFBL mediated p53-dependent apoptosis and recovery was suppressed by promoting G1/S phase arrest and failure of cellular tight junctions. In comparison, intra-cellular anti-oxidant protection exerted by LFBL was absent for extra-cellular HFBL (likely polysaccharide-rich), which potentiated more rapid apoptosis by producing DNA double strand breaks. Differentiated effects on expression of genes regulating ROS defense and chromatic condensation by LFBL versus HFBL, were observed. The results support ferrous iron in cell culture systems and potentially in vivo, can invoke different extra-cellular versus intra-cellular ROS-mediated chemistries, that may be regulated by exogenous, including dietary species

Bioinformatics-Based Identification of Expanded Repeats: A Non-reference Intronic Pentamer Expansion in RFC1 Causes CANVAS

Genomic technologies such as next-generation sequencing (NGS) are revolutionizing molecular diagnostics and clinical medicine. However, these approaches have proven inefficient at identifying pathogenic repeat expansions. Here, we apply a collection of bioinformatics tools that can be utilized to identify either known or novel expanded repeat sequences in NGS data. We performed genetic studies of a cohort of 35 individuals from 22 families with a clinical diagnosis of cerebellar ataxia with neuropathy and bilateral vestibular areflexia syndrome (CANVAS). Analysis of whole-genome sequence (WGS) data with five independent algorithms identified a recessively inherited intronic repeat expansion [(AAGGG)exp] in the gene encoding Replication Factor C1 (RFC1). This motif, not reported in the reference sequence, localized to an Alu element and replaced the reference (AAAAG)11 short tandem repeat. Genetic analyses confirmed the pathogenic expansion in 18 of 22 CANVAS-affected families and identified a core ancestral haplotype, estimated to have arisen in Europe more than twenty-five thousand years ago. WGS of the four RFC1-negative CANVAS-affected families identified plausible variants in three, with genomic re-diagnosis of SCA3, spastic ataxia of the Charlevoix-Saguenay type, and SCA45. This study identified the genetic basis of CANVAS and demonstrated that these improved bioinformatics tools increase the diagnostic utility of WGS to determine the genetic basis of a heterogeneous group of clinically overlapping neurogenetic disorders

Genome-wide changes conferred by trichostatin A, glucose and metformin in human vascular endothelial cells

© 2014 Dr. Haloom RafehiEndothelial dysfunction is a precursor of cardiovascular disease. It is characterised by the impaired production of the vasodilator nitric oxide and the development of a pro-thrombotic and pro-inflammatory endothelial state. Hyperglycaemia is a condition of elevated blood glucose levels as seen in diabetes and is a contributing factor to the development of endothelial dysfunction. The diabetes drug metformin reverses damage to the vascular endothelium, whereas trichostatin A (TSA), a prototypical HDAC inhibitor, has anti-inflammatory properties. While glucose, and TSA have been the focus of extensive research for decades, the genome-wide effects of these compounds in endothelial cells have not been adequately characterised. Recent advances in high throughput sequencing (HTS) technologies have fast-tracked advancements in the genome-wide study of transcriptional regulation. The ability to look beyond single genes and loci is now readily accessible and affordable. This study utilised HTS to study genome-wide gene changes induced by hyperglycaemia and TSA in human endothelial cells. This genome-wide approach led to the identification of novel mechanisms of action by both compounds, while qRT-PCR was used to show that metformin attenuates changes in expression in hyperglycaemic endothelial cells. HDAC inhibitors such as TSA, a diverse group of clinically relevant compounds, are thought to induce global histone hyperacetylation and therefore increased gene expression. In endothelial cells stimulated with TSA, the integration of gene expression changes (RNA-seq) with global histone acetylation data (ChIP-seq) identified that HDAC inhibition induces equal amounts of increases and decreases in gene expression and genome-wide histone acetylation. Histone deacetylation partially accounted for the anti-inflammatory effects of TSA. Deacetylation was dependent on the loss of histone acetyltransferases (HAT) binding at gene promoters. Indeed, the inhibition of structurally related HATs p300 and CBP prevented TSA-dependent reductions in expression at deacetylated genes such as the inflammatory cytokine IL6. A 50% reduction in CBP protein binding was also observed at the IL6 promoter. Transcription factors YY1, STAT2 and IRF3 were also associated with deacetylated promoters. The use of HTS also identified novel pathways in hyperglycaemic endothelial cells that may be involved in the development of endothelial dysfunction. Hyperglycaemic endothelial cells expressed interferon-response pathway genes such as MX1 and IFI27. Transcription factor analysis implicates the activation of STAT1 and IRF1. Co-treatment of hyperglycaemic cells with metformin prevented glucose-dependent changes in gene expression, including interferon-response genes. Indeed, the effects of metformin in endothelial cells were dependent on glucose levels. In normoglycaemic cells, metformin subtly regulated changes in gene expression. In contrast, metformin was strongly associated with the reversal of gene expression changes induced by hyperglycaemia. In this respect, metformin and TSA differ substantially in their anti-inflammatory properties. Metformin did not suppress pro-inflammatory gene expression in unstimulated endothelial cells, whereas TSA did. Interestingly, TSA and metformin were associated with the suppression of genes regulated by pro-inflammatory STAT and IRF family transcription factors. In conclusion, the use of HTS facilitated the identification of novel responses in endothelial cells that were previously under-appreciated. This has provided new insights into the mechanism of action of HDAC inhibitors. It also suggests a potential mechanism for the development of endothelial dysfunction associated with hyperglycaemia and its reversal by metformin. The use of HTS has clear applications in attaining a deeper understanding of disease mechanisms as well as drug action, which, coupled together, have the potential to improve rational drug design

HDAC Inhibition in Vascular Endothelial Cells Regulates the Expression of ncRNAs

While clinical and pre-clinical trials indicate efficacy of histone deacetylase (HDAC) inhibitors in disease mediated by dynamic lysine modification, the impact on the expression of non-coding RNAs (ncRNAs) remains poorly understood. In this study, we investigate high throughput RNA sequencing data derived from primary human endothelial cells stimulated with HDAC inhibitors suberanilohydroxamic acid (SAHA) and Trichostatin A (TSA). We observe robust regulation of ncRNA expression. Integration of gene expression data with histone 3 lysine 9 and 14 acetylation (H3K9/14ac) and histone 3 lysine 4 trimethylation (H3K4me3) datasets identified complex and class-specific expression of ncRNAs. We show that EP300 target genes are subject to histone deacetylation at their promoter following HDAC inhibition. This deacetylation drives suppression of protein-coding genes. However, long intergenic non-coding RNAs (lincRNAs) regulated by EP300 are activated following HDAC inhibition, despite histone deacetylation. This increased expression was driven by increased H3K4me3 at the gene promoter. For example, elevated promoter H3K4me3 increased lincRNA MALAT1 expression despite broad EP300-associated histone deacetylation. In conclusion, we show that HDAC inhibitors regulate the expression of ncRNA by complex and class-specific epigenetic mechanisms

HDAC inhibition in vascular endothelial cells regulates the expression of ncRNAs

While clinical and pre-clinical trials indicate efficacy of histone deacetylase (HDAC) inhibitors in disease mediated by dynamic lysine modification, the impact on the expression of non-coding RNAs (ncRNAs) remains poorly understood. In this study, we investigate high throughput RNA sequencing data derived from primary human endothelial cells stimulated with HDAC inhibitors suberanilohydroxamic acid (SAHA) and Trichostatin A (TSA). We observe robust regulation of ncRNA expression. Integration of gene expression data with histone 3 lysine 9 and 14 acetylation (H3K9/14ac) and histone 3 lysine 4 trimethylation (H3K4me3) datasets identified complex and class-specific expression of ncRNAs. We show that EP300 target genes are subject to histone deacetylation at their promoter following HDAC inhibition. This deacetylation drives suppression of protein-coding genes. However, long intergenic non-coding RNAs (lincRNAs) regulated by EP300 are activated following HDAC inhibition, despite histone deacetylation. This increased expression was driven by increased H3K4me3 at the gene promoter. For example, elevated promoter H3K4me3 increased lincRNA MALAT1 expression despite broad EP300-associated histone deacetylation. In conclusion, we show that HDAC inhibitors regulate the expression of ncRNA by complex and class-specific epigenetic mechanisms

MeCP2 interacts with chromosomal microRNAs in brain

Although methyl CpG binding domain protein-2 (MeCP2) is commonly understood to function as a silencing factor at methylated DNA sequences, recent studies also show that MeCP2 can bind unmethylated sequences and coordinate gene activation. MeCP2 displays broad binding patterns throughout the genome, with high expression levels similar to histone H1 in neurons. Despite its significant presence in the brain, only subtle gene expression changes occur in the absence of MeCP2. This may reflect a more complex regulatory mechanism of MeCP2 to complement chromatin binding. Using an RNA immunoprecipitation of native chromatin technique, we identify MeCP2 interacting microRNAs in mouse primary cortical neurons. In addition, comparison with mRNA sequencing data from Mecp2-null mice suggests that differentially expressed genes may indeed be targeted by MeCP2-interacting microRNAs. These findings highlight the MeCP2 interaction with microRNAs that may modulate its binding with chromatin and regulate gene expression

Systems approach to the pharmacological actions of HDAC inhibitors reveals EP300 activities and convergent mechanisms of regulation in diabetes

<p>Given the skyrocketing costs to develop new drugs, repositioning of approved drugs, such as histone deacetylase (HDAC) inhibitors, may be a promising strategy to develop novel therapies. However, a gap exists in the understanding and advancement of these agents to meaningful translation for which new indications may emerge. To address this, we performed systems-level analyses of 33 independent HDAC inhibitor microarray studies. Based on network analysis, we identified enrichment for pathways implicated in metabolic syndrome and diabetes (insulin receptor signaling, lipid metabolism, immunity and trafficking). Integration with ENCODE ChIP-seq datasets identified suppression of EP300 target genes implicated in diabetes. Experimental validation indicates reversal of diabetes-associated EP300 target genes in primary vascular endothelial cells derived from a diabetic individual following inhibition of HDACs (by SAHA), EP300, or <i>EP300</i> knockdown. Our computational systems biology approach provides an adaptable framework for the prediction of novel therapeutics for existing disease.</p