Integrated metabolic modelling reveals cell-type specific epigenetic control points of the macrophage metabolic network

in BMC Genomics (2015), 16(809),

Background: The reconstruction of context-specific metabolic models from easily and reliably measurable features such as transcriptomics data will be increasingly important in research and medicine ... [more ▼]

Background: The reconstruction of context-specific metabolic models from easily and reliably measurable features such as transcriptomics data will be increasingly important in research and medicine. Current reconstruction methods suffer from high computational effort and arbitrary threshold setting. Moreover, understanding the underlying epigenetic regulation might allow the identification of putative intervention points within metabolic networks. Genes under high regulatory load from multiple enhancers or super-enhancers are known key genes for disease and cell identity. However, their role in regulation of metabolism and their placement within the metabolic networks has not been studied. Methods: Here we present FASTCORMICS, a fast and robust workflow for the creation of high-quality metabolic models from transcriptomics data. FASTCORMICS is devoid of arbitrary parameter settings and due to its low computational demand allows cross-validation assays. Applying FASTCORMICS, we have generated models for 63 primary human cell types from microarray data, revealing significant differences in their metabolic networks. Results: To understand the cell type-specific regulation of the alternative metabolic pathways we built multiple models during differentiation of primary human monocytes to macrophages and performed ChIP-Seq experiments for histone H3 K27 acetylation (H3K27ac) to map the active enhancers in macrophages. Focusing on the metabolic genes under high regulatory load from multiple enhancers or super-enhancers, we found these genes to show the most cell type-restricted and abundant expression profiles within their respective pathways. Importantly, the high regulatory load genes are associated to reactions enriched for transport reactions and other pathway entry points, suggesting that they are critical regulatory control points for cell type-specific metabolism. Conclusions: By integrating metabolic modelling and epigenomic analysis we have identified high regulatory load as a common feature of metabolic genes at pathway entry points such as transporters within the macrophage metabolic network. Analysis of these control points through further integration of metabolic and gene regulatory networks in various contexts could be beneficial in multiple fields from identification of disease intervention strategies to cellular reprogramming. [less ▲]

Transcriptomics profiling of human SGBS adipogenesis

in Genomics Data (2014), 2

Obesity is an ever-growing epidemic where tissue homeostasis is influenced by the differentiation of adipocytes that function in lipid metabolism, endocrine and inflammatory processes. While this ... [more ▼]

Obesity is an ever-growing epidemic where tissue homeostasis is influenced by the differentiation of adipocytes that function in lipid metabolism, endocrine and inflammatory processes. While this differentiation process has been well-characterized in mice, limited data is available from human cells. Applying microarray expression profiling in the human SGBS pre-adipocyte cell line, we identified genes with differential expression during differentiation in combination with constraint-based modeling of metabolic pathway activity. Here we describe the experimental design and quality controls in detail for the gene expression and related results published by Galhardo et al. in Nucleic Acids Research 2014 associated with the data uploaded to NCBI Gene Expression Omnibus (). [less ▲]

Neurodegeneration by Activation of the Microglial Complement–Phagosome Pathway

in Journal of Neuroscience (2014), 34(25), 8546-8556

Systemic inflammatory reactions have been postulated to exacerbate neurodegenerative diseases via microglial activation. We now demonstrate in vivo that repeated systemic challenge of mice over four ... [more ▼]

Systemic inflammatory reactions have been postulated to exacerbate neurodegenerative diseases via microglial activation. We now demonstrate in vivo that repeated systemic challenge of mice over four consecutive days with bacterial LPS maintained an elevated microglial inflammatory phenotype and induced loss of dopaminergic neurons in the substantia nigra. The same total cumulative LPS dose given within a single application did not induce neurodegeneration. Whole-genome transcriptome analysis of the brain demonstrated that repeated systemic LPS application induced an activation pattern involving the classical complement system and its associated phagosome pathway. Loss of dopaminergic neurons induced by repeated systemic LPS application was rescued in complement C3-deficient mice, confirming the involvement of the complement system in neurodegeneration. Our data demonstrate that a phagosomal inflammatory response of microglia is leading to complement-mediated loss of dopaminergic neurons. [less ▲]

Integrated analysis of transcript-level regulation of metabolism reveals disease-relevant nodes of the human metabolic network

in Nucleic Acids Research (2013)

Metabolic diseases and comorbidities represent an ever-growing epidemic where multiple cell types impact tissue homeostasis. Here, the link between the metabolic and gene regulatory networks was studied ... [more ▼]

Metabolic diseases and comorbidities represent an ever-growing epidemic where multiple cell types impact tissue homeostasis. Here, the link between the metabolic and gene regulatory networks was studied through experimental and computational analysis. Integrating gene regulation data with a human metabolic network prompted the establishment of an open-sourced web portal, IDARE (Integrated Data Nodes of Regulation), for visualizing various gene-related data in context of metabolic pathways. Motivated by increasing availability of deep sequencing studies, we obtained ChIP-seq data from widely studied human umbilical vein endothelial cells. Interestingly, we found that association of metabolic genes with multiple transcription factors (TFs) enriched disease-associated genes. To demonstrate further extensions enabled by examining these networks together, constraintbased modeling was applied to data from human preadipocyte differentiation. In parallel, data on gene expression, genome-wide ChIP-seq profiles for peroxisome proliferator-activated receptor (PPAR) c, CCAAT/enhancer binding protein (CEBP) a, liver X receptor (LXR) and H3K4me3 and microRNA target identification for miR-27a, miR-29a and miR-222 were collected. Disease-relevant key nodes, including mitochondrial glycerol-phosphateacyltransferase (GPAM), were exposed from metabolic pathways predicted to change activity by focusing on association with multiple regulators. In both cell types, our analysis reveals the convergence of microRNAs and TFs within the branched chain amino acid (BCAA) metabolic pathway, possibly providing an explanation for its downregulation in obese and diabetic conditions. [less ▲]

Dataset integration identifies transcriptional regulation of microRNA genes by PPARgamma in differentiating mouse 3T3-L1 adipocytes

in Nucleic Acids Research (2012), 40(10), 4446-4460

Peroxisome proliferator-activated receptor gamma (PPARgamma) is a key transcription factor in mammalian adipogenesis. Genome-wide approaches have identified thousands of PPARgamma binding sites in mouse ... [more ▼]

Peroxisome proliferator-activated receptor gamma (PPARgamma) is a key transcription factor in mammalian adipogenesis. Genome-wide approaches have identified thousands of PPARgamma binding sites in mouse adipocytes and PPARgamma upregulates hundreds of protein-coding genes during adipogenesis. However, no microRNA (miRNA) genes have been identified as primary PPARgamma-targets. By integration of four separate datasets of genome-wide PPARgamma binding sites in 3T3-L1 adipocytes we identified 98 miRNA clusters with PPARgamma binding within 50 kb from miRNA transcription start sites. Nineteen mature miRNAs were upregulated >/=2-fold during adipogenesis and for six of these miRNA loci the PPARgamma binding sites were confirmed by at least three datasets. The upregulation of five miRNA genes miR-103-1 (host gene Pank3), miR-148b (Copz1), miR-182/96/183, miR-205 and miR-378 (Ppargc1b) followed that of Pparg. The PPARgamma-dependence of four of these miRNA loci was demonstrated by PPARgamma knock-down and the loci of miR-103-1 (Pank3), miR-205 and miR-378 (Ppargc1b) were also responsive to the PPARgamma ligand rosiglitazone. Finally, chromatin immunoprecipitation analysis validated in silico predicted PPARgamma binding sites at all three loci and H3K27 acetylation was analyzed to confirm the activity of these enhancers. In conclusion, we identified 22 putative PPARgamma target miRNA genes, showed the PPARgamma dependence of four of these genes and demonstrated three as direct PPARgamma target genes in mouse adipogenesis. [less ▲]

Dicer is associated with ribosomal DNA chromatin in mammalian cells

in PLoS ONE (2010), 5(8), 1-2

BACKGROUND: RNA silencing is a common term for pathways utilizing small RNAs as sequence-specific guides to repress gene expression. Components of the RNA silencing machinery are involved in different ... [more ▼]

BACKGROUND: RNA silencing is a common term for pathways utilizing small RNAs as sequence-specific guides to repress gene expression. Components of the RNA silencing machinery are involved in different aspects of chromatin function in numerous organisms. However, association of RNA silencing with chromatin in mammalian cells remains unclear. METHODOLOGY/PRINCIPAL FINDINGS: Immunostaining of mitotic chromosomes with antibodies visualizing either endogenous or ectopically expressed Dicer in mammalian cells revealed association of the protein with ribosomal DNA (rDNA) repeats. Chromatin immunoprecipitations and bisulfite sequencing experiments indicated that Dicer is associated with transcribed regions of both active and silenced genes in rDNA arrays of interphase chromosomes. Metabolic labeling of the mouse embryonic stem (ES) cells lacking Dicer did not reveal apparent defect in rRNA biogenesis though pre-rRNA synthesis in these cells was decreased, likely as a consequence of their slower growth caused by the loss of miRNAs. We analyzed in detail chromatin structure of rDNA but did not find any epigenetic changes at rDNA loci in Dicer(-/-) ES cells. Instead, we found that rDNA methylation is rather low in primary tissues, contrasting with rDNA methylation patterns in transformed cell lines. CONCLUSION/SIGNIFICANCE: We found that Dicer, a key component of RNA silencing pathways, can be detected in association with rDNA chromatin in mammalian cells. The role of this particular localization of Dicer is not readily apparent since the enzyme is associated with rDNA genes regardless of their transcriptional activity. However, localization of Dicer to the transcribed region suggests that transcription may contribute to the Dicer deposition at rDNA chromatin. We hypothesize that Dicer functions in maintaining integrity of rDNA arrays. [less ▲]

Spatio-temporal activation of chromatin on the human CYP24 gene promoter in the presence of 1alpha,25-Dihydroxyvitamin D3

in Journal of Molecular Biology (2005), 350(1), 65-77

The vitamin D3 24-hydroxylase gene (CYP24) is one of the most strongly induced genes known. Despite this, its induction by the hormone 1alpha,25-dihydroxyvitamin D3 (1alpha,25OH2D3) has been characterized ... [more ▼]

The vitamin D3 24-hydroxylase gene (CYP24) is one of the most strongly induced genes known. Despite this, its induction by the hormone 1alpha,25-dihydroxyvitamin D3 (1alpha,25OH2D3) has been characterized only partially. Therefore, we monitored the spatio-temporal, 1alpha,25OH2D3-dependent chromatin acetylation status of the human CYP24 promoter by performing chromatin immunoprecipitation (ChIP) assays with antibodies against acetylated histone 4. This was achieved by performing PCR on 25 contiguous genomic regions spanning the first 7.7 kb of the promoter. ChIP assays using antibodies against the 1alpha,25OH2D3 receptor (VDR) revealed that, in addition to the proximal promoter, three novel regions further upstream associated with VDR. Combined in silico/in vitro screening identified in three of the four promoter regions sequences resembling known VDREs and reporter gene assays confirmed the inducibility of these regions by 1alpha,25OH2D3)=. In contrast, the fourth VDR-associated promoter region did not contain any recognizable classical VDRE that could account for the presence of the protein on this region. However, re-ChIP assays monitored on all four promoter regions simultaneous association of VDR with retinoid X receptor, coactivator, mediator and RNA polymerase II proteins. These proteins showed a promoter region-specific association pattern demonstrating the complex choreography of the CYP24 gene promoter activation over 300 minutes. Thus, this study reveals new information concerning the regulation of the CYP24 gene by 1alpha,25OH2D3, and is a demonstration of the simultaneous participation of multiple, structurally diverse response elements in promoter activation in a living cell. [less ▲]