Analysis of the mechanisms mediating tumor-specific changes in gene expression in human liver tumors

DNA methylation of transposons pattern aging differences across a diverse cohort of dogs from the Dog Aging Project

Within a species, larger individuals often have shorter lives and higher rates of age-related disease. Despite this well-known link, we still know little about underlying age-related epigenetic differences, which could help us better understand inter-individual variation in aging and the etiology, onset, and progression of age-associated disease. Dogs exhibit this negative correlation between size, health, and longevity and thus represent an excellent system in which to test the underlying mechanisms. Here, we quantified genome-wide DNA methylation in a cohort of 864 dogs in the Dog Aging Project. Age strongly patterned the dog epigenome, with the majority (66% of age-associated loci) of regions associating age-related loss of methylation. These age effects were non-randomly distributed in the genome and differed depending on genomic context. We found the LINE1 (long interspersed elements) class of TEs (transposable elements) were the most frequently hypomethylated with age (FDR < 0.05, 40% of all LINE1 regions). This LINE1 pattern differed in magnitude across breeds of different sizes- the largest dogs lost 0.26% more LINE1 methylation per year than the smallest dogs. This suggests that epigenetic regulation of TEs, particularly LINE1s, may contribute to accelerated age and disease phenotypes within a species. Since our study focused on the methylome of immune cells, we looked at LINE1 methylation changes in golden retrievers, a breed highly susceptible to hematopoietic cancers, and found they have accelerated age-related LINE1 hypomethylation compared to other breeds. We also found many of the LINE1s hypomethylated with age are located on the X chromosome and are, when considering X chromosome inactivation, counter-intuitively more methylated in males. These results have revealed the demethylation of LINE1 transposons as a potential driver of inter-species, demographic-dependent aging variation.

H3K27me3-modulated Hofbauer cell BMP2 signalling enhancement compensates for shallow trophoblast invasion in preeclampsia

BACKGROUND

Preeclampsia (PE) is a common hypertensive pregnancy disorder associated with shallow trophoblast invasion. Although bone morphogenetic protein 2 (BMP2) has been shown to promote trophoblast invasion in vitro, its cellular origin and molecular regulation in placenta, as well as its potential role in PE, has yet to be established. Additionally, whether BMP2 and/or its downstream molecules could serve as potential diagnostic or therapeutic targets for PE has not been explored.

METHODS

Placentas and sera from PE and healthy pregnant women were subjected to multi-omics analyses, immunoblots, qPCR, and ELISA assays. Immortalized trophoblast cells, primary cultures of human trophoblasts, and first-trimester villous explants were used for in vitro experiments. Adenovirus expressing sFlt-1 (Ad Flt1)-induced PE rat model was used for in vivo studies.

FINDINGS

We find globally decreased H3K27me3 modifications and increased BMP2 signalling in preeclamptic placentas, which is negatively correlated with clinical manifestations. BMP2 is derived from Hofbauer cells and epigenetically regulated by H3K27me3 modification. BMP2 promotes trophoblast invasion and vascular mimicry by upregulating BMP6 via BMPR1A-SMAD2/3-SMAD4 signalling. BMP2 supplementation alleviates high blood pressure and fetal growth restriction phenotypes in Ad Flt1-induced rat PE model.

INTERPRETATION

Our findings demonstrate that epigenetically regulated Hofbauer cell-derived BMP2 signalling enhancement in late gestation could serve as a compensatory response for shallow trophoblast invasion in PE, suggesting opportunities for diagnostic marker and therapeutic target applications in PE clinical management.

FUNDING

National Key Research and Development Program of China (2022YFC2702400), National Natural Science Foundation of China (82101784, 82171648, 31988101), and Natural Science Foundation of Shandong Province (ZR2020QH051, ZR2020MH039).

SETD2 regulates chromatin accessibility and transcription to suppress lung tumorigenesis

SETD2, a H3K36 trimethyltransferase, is the most frequently mutated epigenetic modifier in lung adenocarcinoma, with a mutation frequency of approximately 9%. However, how SETD2 loss of function promotes tumorigenesis remains unclear. Using conditional Setd2-KO mice, we demonstrated that Setd2 deficiency accelerated the initiation of KrasG12D-driven lung tumorigenesis, increased tumor burden, and significantly reduced mouse survival. An integrated chromatin accessibility and transcriptome analysis revealed a potentially novel tumor suppressor model of SETD2 in which SETD2 loss activates intronic enhancers to drive oncogenic transcriptional output, including the KRAS transcriptional signature and PRC2-repressed targets, through regulation of chromatin accessibility and histone chaperone recruitment. Importantly, SETD2 loss sensitized KRAS-mutant lung cancer to inhibition of histone chaperones, the FACT complex, or transcriptional elongation both in vitro and in vivo. Overall, our studies not only provide insight into how SETD2 loss shapes the epigenetic and transcriptional landscape to promote tumorigenesis, but they also identify potential therapeutic strategies for SETD2 mutant cancers.

PI3Kδ activity controls plasticity and discriminates between EMT and stemness based on distinct TGFβ signaling

The stem cells involved in formation of the complex human body are epithelial cells that undergo apicobasal polarization and form a hollow lumen. Epithelial plasticity manifests as epithelial to mesenchymal transition (EMT), a process by which epithelial cells switch their polarity and epithelial features to adopt a mesenchymal phenotype. The connection between the EMT program and acquisition of stemness is now supported by a substantial number of reports, although what discriminates these two processes remains largely elusive. In this study, based on 3D organoid culture of hepatocellular carcinoma (HCC)-derived cell lines and AAV8-based protein overexpression in the mouse liver, we show that activity modulation of isoform δ of phosphoinositide 3-kinase (PI3Kδ) controls differentiation and discriminates between stemness and EMT by regulating the transforming growth factor β (TGFβ) signaling. This study provides an important tool to control epithelial cell fate and represents a step forward in understanding the development of aggressive carcinoma.

Direct Conversion of Human Endothelial Cells Into Liver Cancer-Forming Cells Using Nonintegrative Episomal Vectors

Liver cancer is an aggressive cancer associated with a poor prognosis. Development of therapeutic strategies for liver cancer requires fundamental research using suitable experimental models. Recent progress in direct reprogramming technology has enabled the generation of many types of cells that are difficult to obtain and provide a cellular resource in experimental models of human diseases. In this study, we aimed to establish a simple one-step method for inducing cells that can form malignant human liver tumors directly from healthy endothelial cells using nonintegrating episomal vectors. To screen for factors capable of inducing liver cancer-forming cells (LCCs), we selected nine genes and one short hairpin RNA that suppresses tumor protein p53 (TP53) expression and introduced them into human umbilical vein endothelial cells (HUVECs), using episomal vectors. To identify the essential factors, we examined the effect of changing the amounts and withdrawing individual factors. We then analyzed the proliferation, gene and protein expression, morphologic and chromosomal abnormality, transcriptome, and tumor formation ability of the induced cells. We found that a set of six factors, forkhead box A3 (FOXA3), hepatocyte nuclear factor homeobox 1A (HNF1A), HNF1B, lin-28 homolog B (LIN28B), MYCL proto-oncogene, bHLH transcription factor (L-MYC), and Kruppel-like factor 5 (KLF5), induced direct conversion of HUVECs into LCCs. The gene expression profile of these induced LCCs (iLCCs) was similar to that of human liver cancer cells, and these cells effectively formed tumors that resembled human combined hepatocellular-cholangiocarcinoma following transplantation into immunodeficient mice. Conclusion: We succeeded in the direct induction of iLCCs from HUVECs by using nonintegrating episomal vectors. iLCCs generated from patients with cancer and healthy volunteers will be useful for further advancements in cancer research and for developing methods for the diagnosis, treatment, and prognosis of liver cancer.

Transcriptional down-regulation of metabolic genes by Gdown1 ablation induces quiescent cell re-entry into the cell cycle

Here, Jishage et al. show that hepatocyte-specific ablation of RNA polymerase II (Pol II)-associated Gdown1 leads to down-regulation of highly expressed genes involved in plasma protein synthesis and metabolism, a concomitant cell cycle re-entry associated with induction of cell cycle-related genes (including cyclin D1). Their findings establish an important physiological function for a Pol II regulatory factor (Gdown1) in the maintenance of normal liver cell transcription through constraints on cell cycle re-entry of quiescent hepatocytes.

Liver regeneration and metabolism are highly interconnected. Here, we show that hepatocyte-specific ablation of RNA polymerase II (Pol II)-associated Gdown1 leads to down-regulation of highly expressed genes involved in plasma protein synthesis and metabolism, a concomitant cell cycle re-entry associated with induction of cell cycle-related genes (including cyclin D1), and up-regulation of p21 through activation of p53 signaling. In the absence of p53, Gdown1-deficient hepatocytes show a severe dysregulation of cell cycle progression, with incomplete mitoses, and a premalignant-like transformation. Mechanistically, Gdown1 is associated with elongating Pol II on the highly expressed genes and its ablation leads to reduced Pol II recruitment to these genes, suggesting that Pol II redistribution may facilitate hepatocyte re-entry into the cell cycle. These results establish an important physiological function for a Pol II regulatory factor (Gdown1) in the maintenance of normal liver cell transcription through constraints on cell cycle re-entry of quiescent hepatocytes.

Transcriptomic profiling identifies novel mechanisms of transcriptional regulation of the cytochrome P450 (Cyp)3a11 gene

Cytochrome P450 (CYP)3A is the most abundant CYP enzyme in the human liver, and a functional impairment of this enzyme leads to unanticipated adverse reactions and therapeutic failures; these reactions result in the early termination of drug development or the withdrawal of drugs from the market. The transcriptional regulation mechanism of the Cyp3a gene is not fully understood and requires a thorough investigation. We mapped the transcriptome of the Cyp3a gene in a mouse model. The Cyp3a gene was induced using the mPXR activator pregnenolone-16alpha-carbonitrile (PCN) and was subsequently downregulated using lipopolysaccharide (LPS). Our objective was to identify the transcription factors (TFs), epigenetic modulators and molecular pathways that are enriched or repressed by PCN and LPS based on a gene set enrichment analysis. Our analysis shows that 113 genes were significantly upregulated (by at least 1.5-fold) with PCN treatment, and that 834 genes were significantly downregulated (by at least 1.5-fold) with LPS treatment. Additionally, the targets of the 536 transcription factors were enriched by a combined treatment of PCN and LPS, and among these, 285 were found to have binding sites on Cyp3a11. Moreover, the repressed targets of the epigenetic markers HDAC1, HDAC3 and EZH2 were further suppressed by LPS treatment and were enhanced by PCN treatment. By identifying and contrasting the transcriptional regulators that are altered by PCN and LPS, our study provides novel insights into the transcriptional regulation of CYP3A in the liver.

Transcriptomic analysis and mutational status of IDH1 in paired primary-recurrent intrahepatic cholangiocarcinoma

Background

Effective target therapies for intrahepatic cholangiocarcinoma (ICC) have not been identified so far. One of the reasons may be the genetic evolution from primary (PR) to recurrent (REC) tumors. We aim to identify peculiar characteristics and to select potential targets specific for recurrent tumors.

Eighteen ICC paired PR and REC tumors were collected from 5 Italian Centers. Eleven pairs were analyzed for gene expression profiling and 16 for mutational status of IDH1. For one pair, deep mutational analysis by Next Generation Sequencing was also carried out. An independent cohort of patients was used for validation.

Results

Two class-paired comparison yielded 315 differentially expressed genes between REC and PR tumors. Up-regulated genes in RECs are involved in RNA/DNA processing, cell cycle, epithelial to mesenchymal transition (EMT), resistance to apoptosis, and cytoskeleton remodeling. Down-regulated genes participate to epithelial cell differentiation, proteolysis, apoptotic, immune response, and inflammatory processes. A 24 gene signature is able to discriminate RECs from PRs in an independent cohort; FANCG is statistically associated with survival in the chol-TCGA dataset. IDH1 was mutated in the RECs of five patients; 4 of them displayed the mutation only in RECs. Deep sequencing performed in one patient confirmed the IDH1 mutation in REC.

Conclusions

RECs are enriched for genes involved in EMT, resistance to apoptosis, and cytoskeleton remodeling. Key players of these pathways might be considered druggable targets in RECs. IDH1 is mutated in 30% of RECs, becoming both a marker of progression and a target for therapy.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4829-0) contains supplementary material, which is available to authorized users.

DNMT1 is a required genomic regulator for murine liver histogenesis and regeneration

DNA methyltransferase 1 (DNMT1) is an essential regulator maintaining both epigenetic reprogramming during DNA replication and genome stability. We investigated the role of DNMT1 in the regulation of postnatal liver histogenesis under homeostasis and stress conditions. We generated Dnmt1 conditional knockout mice (Dnmt1(Δalb) ) by crossing Dnmt1(fl/fl) with albumin-cyclization recombination transgenic mice. Serum, liver tissues, and primary hepatocytes were collected from 1-week-old to 20-week old mice. The Dnmt1(Δalb) phenotype was assessed by histology, confocal and electron microscopy, biochemistry, as well as transcriptome and methylation profiling. Regenerative growth was induced by partial hepatectomy and exposure to carbon tetrachloride. The impact of Dnmt1 knockdown was also analyzed in hepatic progenitor cell lines; proliferation, apoptosis, DNA damage, and sphere formation were assessed. Dnmt1 loss in postnatal hepatocytes caused global hypomethylation, enhanced DNA damage response, and initiated a senescence state causing a progressive inability to maintain tissue homeostasis and proliferate in response to injury. The liver regenerated through activation and repopulation from progenitors due to lineage-dependent differences in albumin-cyclization recombination expression, providing a basis for selection of less mature and therefore less damaged hepatic progenitor cell progeny. Consistently, efficient knockdown of Dnmt1 in cultured hepatic progenitor cells caused severe DNA damage, cell cycle arrest, senescence, and cell death. Mx1-cyclization recombination-driven deletion of Dnmt1 in adult quiescent hepatocytes did not affect liver homeostasis.

CONCLUSION

These results establish the indispensable role of DNMT1-mediated epigenetic regulation in postnatal liver growth and regeneration; Dnmt1(Δalb) mice provide a unique experimental model to study the role of senescence and the contribution of progenitor cells to physiological and regenerative liver growth. (Hepatology 2016;64:582-598).

Locally disordered methylation forms the basis of intratumor methylome variation in chronic lymphocytic leukemia

Intratumoral heterogeneity plays a critical role in tumor evolution. To define the contribution of DNA methylation to heterogeneity within tumors, we performed genome-scale bisulfite sequencing of 104 primary chronic lymphocytic leukemias (CLLs). Compared with 26 normal B cell samples, CLLs consistently displayed higher intrasample variability of DNA methylation patterns across the genome, which appears to arise from stochastically disordered methylation in malignant cells. Transcriptome analysis of bulk and single CLL cells revealed that methylation disorder was linked to low-level expression. Disordered methylation was further associated with adverse clinical outcome. We therefore propose that disordered methylation plays a similar role to that of genetic instability, enhancing the ability of cancer cells to search for superior evolutionary trajectories.

The R2TP chaperone complex: its involvement in snoRNP assembly and tumorigenesis

R2TP was originally identified in yeast Saccharomyces cerevisiae as Hsp90 interacting complex, and is composed of four different proteins: Rvb1, Rvb2, Tah1, and Pih1. This complex is well-conserved in eukaryotes, and is involved in many cellular processes such as snoRNP biogenesis, RNA polymerase assembly, PIKK signaling, and apoptosis. An increasing number of research related to R2TP suggests a linkage of its function with tumorigenesis. In this review, we provide an overview of several recent studies on R2TP that are related to cell proliferation and carcinogenesis, and propose a possible role of R2TP in tumorigenesis through regulating snoRNA/snoRNP biogenesis.

Differences among brain tumor stem cell types and fetal neural stem cells in focal regions of histone modifications and DNA methylation, broad regions of modifications, and bivalent promoters

Background

Aberrational epigenetic marks are believed to play a major role in establishing the abnormal features of cancer cells. Rational use and development of drugs aimed at epigenetic processes requires an understanding of the range, extent, and roles of epigenetic reprogramming in cancer cells. Using ChIP-chip and MeDIP-chip approaches, we localized well-established and prevalent epigenetic marks (H3K27me3, H3K4me3, H3K9me3, DNA methylation) on a genome scale in several lines of putative glioma stem cells (brain tumor stem cells, BTSCs) and, for comparison, normal human fetal neural stem cells (fNSCs).

Results

We determined a substantial “core” set of promoters possessing each mark in every surveyed BTSC cell type, which largely overlapped the corresponding fNSC sets. However, there was substantial diversity among cell types in mark localization. We observed large differences among cell types in total number of H3K9me3+ positive promoters and peaks and in broad modifications (defined as >50 kb peak length) for H3K27me3 and, to a lesser extent, H3K9me3. We verified that a change in a broad modification affected gene expression of CACNG7. We detected large numbers of bivalent promoters, but most bivalent promoters did not display direct overlap of contrasting epigenetic marks, but rather occupied nearby regions of the proximal promoter. There were significant differences in the sets of promoters bearing bivalent marks in the different cell types and few consistent differences between fNSCs and BTSCs.

Conclusions

Overall, our “core set” data establishes sets of potential therapeutic targets, but the diversity in sets of sites and broad modifications among cell types underscores the need to carefully consider BTSC subtype variation in epigenetic therapy. Our results point toward substantial differences among cell types in the activity of the production/maintenance systems for H3K9me3 and for broad regions of modification (H3K27me3 or H3K9me3). Finally, the unexpected diversity in bivalent promoter sets among these multipotent cells indicates that bivalent promoters may play complex roles in the overall biology of these cells. These results provide key information for forming the basis for future rational drug therapy aimed at epigenetic processes in these cells.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-724) contains supplementary material, which is available to authorized users.

Genomic profiling of hepatocellular adenomas reveals recurrent FRK-activating mutations and the mechanisms of malignant transformation

Hepatocellular adenomas (HCA) are benign liver tumors predominantly developed in women using oral contraceptives. Here, exome sequencing identified recurrent somatic FRK mutations that induce constitutive kinase activity, STAT3 activation, and cell proliferation sensitive to Src inhibitors. We also found uncommon recurrent mutations activating JAK1, gp130, or β-catenin. Chromosome copy number and methylation profiling revealed patterns that correlated with specific gene mutations and tumor phenotypes. Finally, integrative analysis of HCAs transformed to hepatocellular carcinoma revealed β-catenin mutation as an early alteration and TERT promoter mutations as associated with the last step of the adenoma-carcinoma transition. In conclusion, we identified the genomic diversity in benign hepatocyte proliferation, several therapeutic targets, and the key genomic determinants of the adenoma-carcinoma transformation sequence.

Development of biomarkers for screening hepatocellular carcinoma using global data mining and multiple reaction monitoring

Hepatocellular carcinoma (HCC) is one of the most common and aggressive cancers and is associated with a poor survival rate. Clinically, the level of alpha-fetoprotein (AFP) has been used as a biomarker for the diagnosis of HCC. The discovery of useful biomarkers for HCC, focused solely on the proteome, has been difficult; thus, wide-ranging global data mining of genomic and proteomic databases from previous reports would be valuable in screening biomarker candidates. Further, multiple reaction monitoring (MRM), based on triple quadrupole mass spectrometry, has been effective with regard to high-throughput verification, complementing antibody-based verification pipelines. In this study, global data mining was performed using 5 types of HCC data to screen for candidate biomarker proteins: cDNA microarray, copy number variation, somatic mutation, epigenetic, and quantitative proteomics data. Next, we applied MRM to verify HCC candidate biomarkers in individual serum samples from 3 groups: a healthy control group, patients who have been diagnosed with HCC (Before HCC treatment group), and HCC patients who underwent locoregional therapy (After HCC treatment group). After determining the relative quantities of the candidate proteins by MRM, we compared their expression levels between the 3 groups, identifying 4 potential biomarkers: the actin-binding protein anillin (ANLN), filamin-B (FLNB), complementary C4-A (C4A), and AFP. The combination of 2 markers (ANLN, FLNB) improved the discrimination of the before HCC treatment group from the healthy control group compared with AFP. We conclude that the combination of global data mining and MRM verification enhances the screening and verification of potential HCC biomarkers. This efficacious integrative strategy is applicable to the development of markers for cancer and other diseases.

Identifying master regulators of cancer and their downstream targets by integrating genomic and epigenomic features

Vast amounts of molecular data characterizing the genome, epigenome and transcriptome are becoming available for a variety of cancers. The current challenge is to integrate these diverse layers of molecular biology information to create a more comprehensive view of key biological processes underlying cancer. We developed a biocomputational algorithm that integrates copy number, DNA methylation, and gene expression data to study master regulators of cancer and identify their targets. Our algorithm starts by generating a list of candidate driver genes based on the rationale that genes that are driven by multiple genomic events in a subset of samples are unlikely to be randomly deregulated. We then select the master regulators from the candidate driver and identify their targets by inferring the underlying regulatory network of gene expression. We applied our biocomputational algorithm to identify master regulators and their targets in glioblastoma multiforme (GBM) and serous ovarian cancer. Our results suggest that the expression of candidate drivers is more likely to be influenced by copy number variations than DNA methylation. Next, we selected the master regulators and identified their downstream targets using module networks analysis. As a proof-of-concept, we show that the GBM and ovarian cancer module networks recapitulate known processes in these cancers. In addition, we identify master regulators that have not been previously reported and suggest their likely role. In summary, focusing on genes whose expression can be explained by their genomic and epigenomic aberrations is a promising strategy to identify master regulators of cancer.

Integrative genomics: liver regeneration and hepatocellular carcinoma

Numerous genome wide profiles of gene expression changes in human hepatocellular carcinoma (HCC), compared to normal liver tissue, have been reported. Hierarchical clustering of these data reveal distinct patterns, which underscore conservation between human disease and mouse models of HCC, as well as suggest specific classification of subtypes within the heterogeneous disease of HCC. Global profiling of gene expression in mouse liver, challenged by partial hepatectomy to regenerate, reveals alterations in gene expression that occur in response to acute injury, inflammation, and re-entry into cell cycle. When we integrated datasets of gene expression changes in mouse models of HCC and those that are altered at specific times of liver regeneration, we saw shared, conserved alterations in gene expression within specific biological pathways, both up-regulated, for example, cell cycle, cell death, and cellular development, or down-regulated, for example, vitamin and mineral metabolism, lipid metabolism, and molecular transport. Additional molecular mechanisms shared by liver regeneration and HCC, as yet undiscovered, may have important implications in tumor development and recurrence. These comparisons may offer a way to judge how liver resection, in the treatment of HCC, introduces challenges to care of the disease. Further, uncovering the pathways conserved in inflammatory response, hypertrophy, proliferation, and architectural remodeling of the liver, which are shared in liver regeneration and HCC, versus those specific to tumor development and progression in HCC, may reveal new biomarkers or potential therapeutic targets in HCC.

Pluripotent stem cells escape from senescence-associated DNA methylation changes

Pluripotent stem cells evade replicative senescence, whereas other primary cells lose their proliferation and differentiation potential after a limited number of cell divisions, and this is accompanied by specific senescence-associated DNA methylation (SA-DNAm) changes. Here, we investigate SA-DNAm changes in mesenchymal stromal cells (MSC) upon long-term culture, irradiation-induced senescence, immortalization, and reprogramming into induced pluripotent stem cells (iPSC) using high-density HumanMethylation450 BeadChips. SA-DNAm changes are highly reproducible and they are enriched in intergenic and nonpromoter regions of developmental genes. Furthermore, SA-hypomethylation in particular appears to be associated with H3K9me3, H3K27me3, and Polycomb-group 2 target genes. We demonstrate that ionizing irradiation, although associated with a senescence phenotype, does not affect SA-DNAm. Furthermore, overexpression of the catalytic subunit of the human telomerase (TERT) or conditional immortalization with a doxycycline-inducible system (TERT and SV40-TAg) result in telomere extension, but do not prevent SA-DNAm. In contrast, we demonstrate that reprogramming into iPSC prevents almost the entire set of SA-DNAm changes. Our results indicate that long-term culture is associated with an epigenetically controlled process that stalls cells in a particular functional state, whereas irradiation-induced senescence and immortalization are not causally related to this process. Absence of SA-DNAm in pluripotent cells may play a central role for their escape from cellular senescence.

A transgenic mouse marking live replicating cells reveals in vivo transcriptional program of proliferation

Most adult mammalian tissues are quiescent, with rare cell divisions serving to maintain homeostasis. At present, the isolation and study of replicating cells from their in vivo niche typically involves immunostaining for intracellular markers of proliferation, causing the loss of sensitive biological material. We describe a transgenic mouse strain, expressing a CyclinB1-GFP fusion reporter, that marks replicating cells in the S/G2/M phases of the cell cycle. Using flow cytometry, we isolate live replicating cells from the liver and compare their transcriptome to that of quiescent cells to reveal gene expression programs associated with cell proliferation in vivo. We find that replicating hepatocytes have reduced expression of genes characteristic of liver differentiation. This reporter system provides a powerful platform for gene expression and metabolic and functional studies of replicating cells in their in vivo niche.

The histone demethylase Jhdm1a regulates hepatic gluconeogenesis

Hepatic gluconeogenesis is required for maintaining blood glucose homeostasis; yet, in diabetes mellitus, this process is unrestrained and is a major contributor to fasting hyperglycemia. To date, the impacts of chromatin modifying enzymes and chromatin landscape on gluconeogenesis are poorly understood. Through catalyzing the removal of methyl groups from specific lysine residues in the histone tail, histone demethylases modulate chromatin structure and, hence, gene expression. Here we perform an RNA interference screen against the known histone demethylases and identify a histone H3 lysine 36 (H3K36) demethylase, Jhdm1a, as a key negative regulator of gluconeogenic gene expression. In vivo, silencing of Jhdm1a promotes liver glucose synthesis, while its exogenous expression reduces blood glucose level. Importantly, the regulation of gluconeogenesis by Jhdm1a requires its demethylation activity. Mechanistically, we find that Jhdm1a regulates the expression of a major gluconeogenic regulator, C/EBPα. This is achieved, at least in part, by its USF1-dependent association with the C/EBPα promoter and its subsequent demethylation of dimethylated H3K36 on the C/EBPα locus. Our work provides compelling evidence that links histone demethylation to transcriptional regulation of gluconeogenesis and has important implications for the treatment of diabetes.

Histones are small proteins that are essential for packaging and ordering genetic information (DNA) into high-order chromatin structures. Methylation of specific lysine residues of histones alters chromatin structure, serving as an important epigenetic mechanism for regulation of gene expression. The dynamic nature of histone methylation is controlled by a balance of methyltransferases and demethylases. We have discovered here that the demethylase Jhdm1a negatively regulates gluconeogenesis (de novo glucose synthesis) through suppressing the expression of two rate-limiting gluconeogenic enzymes. Gluconeogenesis is required for maintaining blood glucose homeostasis; yet, in diabetes mellitus, this process is unrestrained and is a major contributor to hyperglycemia. Indeed, we have found that manipulation of Jhdm1a level in liver affects glucose production in normal mice and hyperglycemia in diabetic mice. Mechanistically, Jhdm1a actively removes dimethyl groups from histone H3K36 along the locus of a key gluconeogenic regulator, C/EBPα, which in turn results in decreased C/EBPα expression. Our findings thus identify histone demethylation as a novel regulatory mechanism for gluconeogenesis and have important implications for the treatment of diabetes.

Workflows for microarray data processing in the Kepler environment

BACKGROUND

Microarray data analysis has been the subject of extensive and ongoing pipeline development due to its complexity, the availability of several options at each analysis step, and the development of new analysis demands, including integration with new data sources. Bioinformatics pipelines are usually custom built for different applications, making them typically difficult to modify, extend and repurpose. Scientific workflow systems are intended to address these issues by providing general-purpose frameworks in which to develop and execute such pipelines. The Kepler workflow environment is a well-established system under continual development that is employed in several areas of scientific research. Kepler provides a flexible graphical interface, featuring clear display of parameter values, for design and modification of workflows. It has capabilities for developing novel computational components in the R, Python, and Java programming languages, all of which are widely used for bioinformatics algorithm development, along with capabilities for invoking external applications and using web services.

RESULTS

We developed a series of fully functional bioinformatics pipelines addressing common tasks in microarray processing in the Kepler workflow environment. These pipelines consist of a set of tools for GFF file processing of NimbleGen chromatin immunoprecipitation on microarray (ChIP-chip) datasets and more comprehensive workflows for Affymetrix gene expression microarray bioinformatics and basic primer design for PCR experiments, which are often used to validate microarray results. Although functional in themselves, these workflows can be easily customized, extended, or repurposed to match the needs of specific projects and are designed to be a toolkit and starting point for specific applications. These workflows illustrate a workflow programming paradigm focusing on local resources (programs and data) and therefore are close to traditional shell scripting or R/BioConductor scripting approaches to pipeline design. Finally, we suggest that microarray data processing task workflows may provide a basis for future example-based comparison of different workflow systems.

CONCLUSIONS

We provide a set of tools and complete workflows for microarray data analysis in the Kepler environment, which has the advantages of offering graphical, clear display of conceptual steps and parameters and the ability to easily integrate other resources such as remote data and web services.

Novel DNA binding domain-based assays for detection of methylated and nonmethylated DNA

Methylation of cytosine in the context of CpG dinucleotides is an epigenetic phenomenon in eukaryotes that plays important roles in genome function and transcription regulation. Aberrant changes in DNA methylation is an important feature of several human diseases such as cancer and neurological disorders. These discoveries have opened a new field of new therapies and diagnostics. During recent years, there has been a revolution in DNA methylation analysis technologies. This article focuses on methods with which to study DNA methylation that employ protein domains that specifically recognize either 5-methyl-cytosine in the CpG context or nonmethylated DNA, and methods developed for the detection of 5-hydroxymethylcytosine, the recently described epigenetic mark known as the sixth base of the epigenome.

Functional heterologous expression and purification of a mammalian methyl-CpG binding domain in suitable yield for DNA methylation profiling assays

DNA methylation is a major epigenetic modification in mammalian cells, and patterns involving methylation of cytosine bases, known as CpG methylation, have been implicated in the development of many types of cancer. Methyl binding domains (MBDs) excised from larger mammalian methyl-CpG-binding proteins specifically recognize methyl-cytosine bases of CpG dinucleotides in duplex DNA. Previous molecular diagnostic studies involving MBDs have employed Escherichia coli for protein expression with either low soluble yields or the use of time-consuming denaturation-renaturation purification procedures to improve yields. Efficient MBD-based diagnostics require expression and purification methods that maximize protein yield and minimize time and resource expenditure. This study is a systematic optimization analysis of MBD expression using both SDS-PAGE and microscopy and it provides a comparison of protein yield from published procedures to that from the conditions found to be optimal in these experiments. Protein binding activity and specificity were verified using a DNA electrophoretic mobility shift assay, and final protein yield was improved from the starting conditions by a factor of 65 with a simple, single-step purification.

Cancer associated epigenetic transitions identified by genome-wide histone methylation binding profiles in human colorectal cancer samples and paired normal mucosa

BACKGROUND

Despite their well-established functional roles, histone modifications have received less attention than DNA methylation in the cancer field. In order to evaluate their importance in colorectal cancer (CRC), we generated the first genome-wide histone modification profiles in paired normal colon mucosa and tumor samples.

METHODS

Chromatin immunoprecipitation and microarray hybridization (ChIP-chip) was used to identify promoters enriched for histone H3 trimethylated on lysine 4 (H3K4me3) and lysine 27 (H3K27me3) in paired normal colon mucosa and tumor samples from two CRC patients and for the CRC cell line HT29.

RESULTS

By comparing histone modification patterns in normal mucosa and tumors, we found that alterations predicted to have major functional consequences were quite rare. Furthermore, when normal or tumor tissue samples were compared to HT29, high similarities were observed for H3K4me3. However, the differences found for H3K27me3, which is important in determining cellular identity, indicates that cell lines do not represent optimal tissue models. Finally, using public expression data, we uncovered previously unknown changes in CRC expression patterns. Genes positive for H3K4me3 in normal and/or tumor samples, which are typically already active in normal mucosa, became hyperactivated in tumors, while genes with H3K27me3 in normal and/or tumor samples and which are expressed at low levels in normal mucosa, became hypersilenced in tumors.

CONCLUSIONS

Genome wide histone modification profiles can be used to find epigenetic aberrations in genes associated with cancer. This strategy gives further insights into the epigenetic contribution to the oncogenic process and may identify new biomarkers.

Q-PCR in combination with ChIP assays to detect changes in chromatin acetylation

Quantitative polymerase chain reaction (Q-PCR) allows for the accurate and reproducible determination of the amount of target DNA in a sample through the measurement of PCR product accumulation in "real time." This method determines starting target DNA quantity over a large assay dynamic range and requires no post-PCR sample manipulation. When used in combination with the method of chromatin immunoprecipitation (ChIP), the amount of protein binding to a specific region of DNA can be accurately and rapidly determined. A method for quantifying the presence of acetylated histones H3 and H4 on different regions of a target locus using Q-PCR after ChIP is described.

Integrative epigenomic and genomic analysis of malignant pheochromocytoma

Epigenomic and genomic changes affect gene expression and contribute to tumor development. The histone modifications trimethylated histone H3 lysine 4 (H3K4me3) and lysine 27 (H3K27me3) are epigenetic regulators associated to active and silenced genes, respectively and alterations of these modifications have been observed in cancer. Furthermore, genomic aberrations such as DNA copy number changes are common events in tumors. Pheochromocytoma is a rare endocrine tumor of the adrenal gland that mostly occurs sporadic with unknown epigenetic/genetic cause. The majority of cases are benign. Here we aimed to combine the genome-wide profiling of H3K4me3 and H3K27me3, obtained by the ChIP-chip methodology, and DNA copy number data with global gene expression examination in a malignant pheochromocytoma sample. The integrated analysis of the tumor expression levels, in relation to normal adrenal medulla, indicated that either histone modifications or chromosomal alterations, or both, have great impact on the expression of a substantial fraction of the genes in the investigated sample. Candidate tumor suppressor genes identified with decreased expression, a H3K27me3 mark and/or in regions of deletion were for instance TGIF1, DSC3, TNFRSF10B, RASSF2, HOXA9, PTPRE and CDH11. More genes were found with increased expression, a H3K4me3 mark, and/or in regions of gain. Potential oncogenes detected among those were GNAS, INSM1, DOK5, ETV1, RET, NTRK1, IGF2, and the H3K27 trimethylase gene EZH2. Our approach to associate histone methylations and DNA copy number changes to gene expression revealed apparent impact on global gene transcription, and enabled the identification of candidate tumor genes for further exploration.

Analysis of post-transcriptional regulation using the FunREG method

An increasing number of arguments, including altered microRNA expression, support the idea that post-transcriptional deregulation participates in gene disturbances found in diseased tissues. To evaluate this hypothesis, we developed a method which facilitates post-transcriptional investigations in a wide range of human cells and experimental conditions. This method, called FunREG (functional, integrated and quantitative method to measure post-transcriptional regulation), connects lentiviral transduction with a fluorescent reporter system and quantitative PCR. Using FunREG, we efficiently measured post-transcriptional regulation mediated either by selected RNA sequences or regulatory factors (microRNAs), and then evaluated the contribution of mRNA decay and translation efficiency in the observed regulation. We demonstrated the existence of gene-specific post-transcriptional deregulation in liver tumour cells, and also reported a molecular link between a transcript variant abrogating HDAC6 (histone deacetylase 6) regulation by miR-433 and a rare familial genetic disease. Because FunREG is sensitive, quantitative and easy to use, many applications can be envisioned in fundamental and pathophysiological research.

BeadDataPackR: A Tool to Facilitate the Sharing of Raw Data from Illumina BeadArray Studies

Microarray technologies have been an increasingly important tool in cancer research in the last decade, and a number of initiatives have sought to stress the importance of the provision and sharing of raw microarray data. Illumina BeadArrays provide a particular problem in this regard, as their random construction simultaneously adds value to analysis of the raw data and obstructs the sharing of those data.We present a compression scheme for raw Illumina BeadArray data, designed to ease the burdens of sharing and storing such data, that is implemented in the BeadDataPackR BioConductor package (http://bioconductor.org/packages/release/bioc/html/BeadDataPackR.html). It offers two key advantages over off-the-peg compression tools. First it uses knowledge of the data formats to achieve greater compression than other approaches, and second it does not need to be decompressed for analysis, but rather the values held within can be directly accessed.

Impairment of pre-mRNA splicing in liver disease: Mechanisms and consequences

Pre-mRNA splicing is an essential step in the process of gene expression in eukaryotes and consists of the removal of introns and the linking of exons to generate mature mRNAs. This is a highly regulated mechanism that allows the alternative usage of exons, the retention of intronic sequences and the generation of exonic sequences of variable length. Most human genes undergo splicing events, and disruptions of this process have been associated with a variety of diseases, including cancer. Hepatocellular carcinoma (HCC) is a molecularly heterogeneous type of tumor that usually develops in a cirrhotic liver. Alterations in pre-mRNA splicing of some genes have been observed in liver cancer, and although still scarce, the available data suggest that splicing defects may have a role in hepatocarcinogenesis. Here we briefly review the general mechanisms that regulate pre-mRNA splicing, and discuss some examples that illustrate how this process is impaired in liver tumorigenesis, and may contribute to HCC development. We believe that a more thorough examination of pre-mRNA splicing is still needed to accurately draw the molecular portrait of liver cancer. This will surely contribute to a better understanding of the disease and to the development of new effective therapies.

Histone H3 lysine 27 methylation asymmetry on developmentally-regulated promoters distinguish the first two lineages in mouse preimplantation embryos

First lineage specification in the mammalian embryo leads to formation of the inner cell mass (ICM) and trophectoderm (TE), which respectively give rise to embryonic and extraembryonic tissues. We show here that this first differentiation event is accompanied by asymmetric distribution of trimethylated histone H3 lysine 27 (H3K27me3) on promoters of signaling and developmentally-regulated genes in the mouse ICM and TE. A genome-wide survey of promoter occupancy by H3K4me3 and H3K27me3 indicates that both compartments harbor promoters enriched in either modification, and promoters co-enriched in trimethylated H3K4 and H3K27 linked to developmental and signaling functions. The majority of H3K4/K27me3 co-enriched promoters are distinct between the two lineages, primarily due to differences in the distribution of H3K27me3. Derivation of embryonic stem cells leads to significant losses and gains of H3K4/K27me3 co-enriched promoters relative to the ICM, with distinct contributions of (de)methylation events on K4 and K27. Our results show histone trimethylation asymmetry on promoters in the first two developmental lineages, and highlight an epigenetic skewing associated with embryonic stem cell derivation.

BeadDatapackR: A Tool to Facilitate the Sharing of Raw Data from Illumina BeadArray Studies

Microarray technologies have been an increasingly important tool in cancer research in the last decade, and a number of initiatives have sought to stress the importance of the provision and sharing of raw microarray data. Illumina BeadArrays provide a particular problem in this regard, as their random construction simultaneously adds value to analysis of the raw data and obstructs the sharing of those data.

We present a compression scheme for raw Illumina BeadArray data, designed to ease the burdens of sharing and storing such data, that is implemented in the BeadDataPackR BioConductor package ( http://bioconductor.org/packages/release/bioc/html/BeadDataPackR.html ). It offers two key advantages over off-the-peg compression tools. First it uses knowledge of the data formats to achieve greater compression than other approaches, and second it does not need to be decompressed for analysis, but rather the values held within can be directly accessed.

New tools for functional genomic analysis

For the past decade, the development of genomic technology has revolutionized modern biological research and drug discovery. Functional genomic analyses enable biologists to perform analysis of genetic events on a global scale and they have been widely used in gene discovery, biomarker determination, disease classification, and drug target identification. In this article, we provide an overview of the current and emerging tools involved in genomic studies, including expression arrays, microRNA arrays, array CGH, ChIP-on-chip, methylation arrays, mutation analysis, genome-wide association studies, proteomic analysis, integrated functional genomic analysis and related bioinformatic and biostatistical analyses. Using human liver cancer as an example, we provide further information of how these genomic approaches can be applied in cancer research.

Fast genomic muChIP-chip from 1,000 cells

A new method for rapid genome-wide μChIP-chip from as few as 1,000 cells.

Genome-wide location analysis of histone modifications and transcription factor binding relies on chromatin immunoprecipitation (ChIP) assays. These assays are, however, time-consuming and require large numbers of cells, hindering their application to the analysis of many interesting cell types. We report here a fast microChIP (μChIP) assay for 1,000 cells in combination with microarrays to produce genome-scale surveys of histone modifications. μChIP-chip reliably reproduces data obtained by large-scale assays: H3K9ac and H3K9m3 enrichment profiles are conserved and nucleosome-free regions are revealed.

E2F in vivo binding specificity: comparison of consensus versus nonconsensus binding sites

We have previously shown that most sites bound by E2F family members in vivo do not contain E2F consensus motifs. However, differences between in vivo target sites that contain or lack a consensus E2F motif have not been explored. To understand how E2F binding specificity is achieved in vivo, we have addressed how E2F family members are recruited to core promoter regions that lack a consensus motif and are excluded from other regions that contain a consensus motif. Using chromatin immunoprecipitation coupled with DNA microarray analysis (ChIP-chip) assays, we have shown that the predominant factors specifying whether E2F is recruited to an in vivo binding site are (1) the site must be in a core promoter and (2) the region must be utilized as a promoter in that cell type. We have tested three models for recruitment of E2F to core promoters lacking a consensus site, including (1) indirect recruitment, (2) looping to the core promoter mediated by an E2F bound to a distal motif, and (3) assisted binding of E2F to a site that weakly resembles an E2F motif. To test these models, we developed a new in vivo assay, termed eChIP, which allows analysis of transcription factor binding to isolated fragments. Our findings suggest that in vivo (1) a consensus motif is not sufficient to recruit E2Fs, (2) E2Fs can bind to isolated regions that lack a consensus motif, and (3) binding can require regions other than the best match to the E2F motif.

Using ChIP-chip technology to reveal common principles of transcriptional repression in normal and cancer cells

We compared 12 different cell populations, including embryonic stem cells before and during differentiation into embryoid bodies as well as various types of normal and tumor cells to determine if pluripotent versus differentiated cell types use different mechanisms to establish their transcriptome. We first identified genes that were not expressed in the 12 different cell populations and then determined which of them were regulated by histone methylation, DNA methylation, at the step of productive elongation, or by the inability to establish a preinitiation complex. For these experiments, we performed chromatin immunoprecipitation using antibodies to H3me3K27, H3me3K9, 5-methyl-cytosine, and POLR2A. We found that (1) the percentage of low expressed genes bound by POLR2A, H3me3K27, H3me3K9, or 5-methyl-cytosine is similar in all 12 cell types, regardless of differentiation or neoplastic state; (2) a gene is generally repressed by only one mechanism; and (3) distinct classes of genes are repressed by certain mechanisms. We further characterized two transitioning cell populations, 3T3 cells progressing from G0/G1 into S phase and mES cells differentiating into embryoid bodies. We found that the transient regulation through the cell cycle was achieved predominantly by changes in the recruitment of the general transcriptional machinery or by post-POLR2A recruitment mechanisms. In contrast, changes in chromatin silencing were critical for the permanent changes in gene expression in cells undergoing differentiation.