Histone modification profiles are predictive for tissue/cell-type specific expression of both protein-coding and microRNA genes

Whole-genome landscape of histone H3K4me3 modification during sperm cell lineage development in tomato

BACKGROUND

During male gametogenesis of flowering plants, sperm cell lineage (microspores, generative cells, and sperm cells) differentiated from somatic cells and acquired different cell fates. Trimethylation of histone H3 on lysine 4 (H3K4me3) epigenetically contributes to this process, however, it remained unclear how H3K4me3 influences the gene expression in each cell type. Here, we conducted chromatin immunoprecipitation sequencing (ChIP-seq) to obtain a genome-wide landscape of H3K4me3 during sperm cell lineage development in tomato (Solanum lycopersicum).

RESULTS

We show that H3K4me3 peaks were mainly enriched in the promoter regions, and intergenic H3K4me3 peaks expanded as sperm cell lineage differentiated from somatic cells. H3K4me3 was generally positively associated with transcript abundance and served as a better indicator of gene expression in somatic and vegetative cells, compared to sperm cell lineage. H3K4me3 was mutually exclusive with DNA methylation at 3' proximal of the transcription start sites. The microspore maintained the H3K4me3 features of somatic cells, while generative cells and sperm cells shared an almost identical H3K4me3 pattern which differed from that of the vegetative cell. After microspore division, significant loss of H3K4me3 in genes related to brassinosteroid and cytokinin signaling was observed in generative cells and vegetative cells, respectively.

CONCLUSIONS

Our results suggest the asymmetric division of the microspore significantly reshapes the genome-wide distribution of H3K4me3. Selective loss of H3K4me3 in genes related to hormone signaling may contribute to functional differentiation of sperm cell lineage. This work provides new resource data for the epigenetic studies of gametogenesis in plants.

Studying phenotypic variation and DNA methylation across development, ecology and evolution in the clonal marbled crayfish: a paradigm for investigating epigenotype-phenotype relationships in macro-invertebrates

Animals can produce different phenotypes from the same genome during development, environmental adaptation and evolution, which is mediated by epigenetic mechanisms including DNA methylation. The obligatory parthenogenetic marbled crayfish, Procambarus virginalis, whose genome and methylome are fully established, proved very suitable to study this issue in detail. Comparison between developmental stages and DNA methylation revealed low expression of Dnmt methylation and Tet demethylation enzymes from the spawned oocyte to the 256 cell embryo and considerably increased expression thereafter. The global 5-methylcytosine level was 2.78% at mid-embryonic development and decreased slightly to 2.41% in 2-year-old adults. Genetically identical clutch-mates raised in the same uniform laboratory setting showed broad variation in morphological, behavioural and life history traits and differences in DNA methylation. The invasion of diverse habitats in tropical to cold-temperate biomes in the last 20 years by the marbled crayfish was associated with the expression of significantly different phenotypic traits and DNA methylation patterns, despite extremely low genetic variation on the whole genome scale, suggesting the establishment of epigenetic ecotypes. The evolution of marbled crayfish from its parent species Procambarus fallax by autotriploidy a few decades ago was accompanied by a significant increase in body size, fertility and life span, a 20% reduction of global DNA methylation and alteration of methylation in hundreds of genes, suggesting that epigenetic mechanisms were involved in speciation and fitness enhancement. The combined analysis of phenotypic traits and DNA methylation across multiple biological contexts in the laboratory and field in marbled crayfish may serve as a blueprint for uncovering the role of epigenetic mechanisms in shaping of phenotypes in macro-invertebrates.

Functional annotations of three domestic animal genomes provide vital resources for comparative and agricultural research

Gene regulatory elements are central drivers of phenotypic variation and thus of critical importance towards understanding the genetics of complex traits. The Functional Annotation of Animal Genomes consortium was formed to collaboratively annotate the functional elements in animal genomes, starting with domesticated animals. Here we present an expansive collection of datasets from eight diverse tissues in three important agricultural species: chicken (Gallus gallus), pig (Sus scrofa), and cattle (Bos taurus). Comparative analysis of these datasets and those from the human and mouse Encyclopedia of DNA Elements projects reveal that a core set of regulatory elements are functionally conserved independent of divergence between species, and that tissue-specific transcription factor occupancy at regulatory elements and their predicted target genes are also conserved. These datasets represent a unique opportunity for the emerging field of comparative epigenomics, as well as the agricultural research community, including species that are globally important food resources.

Binding of an X-Specific Condensin Correlates with a Reduction in Active Histone Modifications at Gene Regulatory Elements

Condensins are evolutionarily conserved protein complexes that are required for chromosome segregation during cell division and genome organization during interphase. In Caenorhabditis elegans, a specialized condensin, which forms the core of the dosage compensation complex (DCC), binds to and represses X chromosome transcription. Here, we analyzed DCC localization and the effect of DCC depletion on histone modifications, transcription factor binding, and gene expression using chromatin immunoprecipitation sequencing and mRNA sequencing. Across the X, the DCC accumulates at accessible gene regulatory sites in active chromatin and not heterochromatin. The DCC is required for reducing the levels of activating histone modifications, including H3K4me3 and H3K27ac, but not repressive modification H3K9me3. In X-to-autosome fusion chromosomes, DCC spreading into the autosomal sequences locally reduces gene expression, thus establishing a direct link between DCC binding and repression. Together, our results indicate that DCC-mediated transcription repression is associated with a reduction in the activity of X chromosomal gene regulatory elements.

Expression of Micro-RNA-492 (MiR-492) in Human Cervical Cancer Cell Lines is Upregulated by Transfection with Wild-Type P53, Irradiation, and 5-Fluorouracil Treatment In Vitro

Background

The status of p53 is critical to the chemoradiosensitivity of cervical cancer cells. Wild-type p53 is essential to orchestrate the cellular response to cytotoxic stimuli. Our previous data illustrated that cervical cancer patients whose specimens overexpressed microR-492 (miR-492) were highly sensitive to concurrent chemoradiation. Although p53 activation has been reported to upregulate miR-492 by a miRNA profiling assay in lung cancer cells, the transcriptional regulation of miR-492 in cervical cancer cells remains poorly understood. Therefore, we aimed to decipher the relationship between p53 and miR-492 in cervical cancer cells.

Material/Methods

The expression of p53 and miR-492 in cervical cancer cell lines was measured by western blot and real-time PCR. After cells were transfected with wild-type p53 plasmid or were treated by irradiation and 5-fluorouracil (5-FU), the expression changes of p53 as well as miR-492 were examined by western blot and real-time PCR. The putative p53 binding site of miR-492 was first analyzed by bioinformatics tools, then validated by chromatin immunoprecipitation and dual-luciferase reporter assays.

Results

We found that miR-492 was upregulated in cells with wild-type p53 compared to cells with mutant p53. Transfection of wild-type p53 plasmid or treatments with cytotoxic reagents including irradiation and 5-FU all induced miR-492 overexpression. Bioinformatics analysis and experimental validations further proved p53 interacted with miR-492 promoter directly.

Conclusions

In cervical cancer cells, p53 activated miR-492 expression transcriptionally.

HebbPlot: an intelligent tool for learning and visualizing chromatin mark signatures

BACKGROUND

Histone modifications play important roles in gene regulation, heredity, imprinting, and many human diseases. The histone code is complex and consists of more than 100 marks. Therefore, biologists need computational tools to characterize general signatures representing the distributions of tens of chromatin marks around thousands of regions.

RESULTS

To this end, we developed a software tool, HebbPlot, which utilizes a Hebbian neural network in learning a general chromatin signature from regions with a common function. Hebbian networks can learn the associations between tens of marks and thousands of regions. HebbPlot presents a signature as a digital image, which can be easily interpreted. Moreover, signatures produced by HebbPlot can be compared quantitatively. We validated HebbPlot in six case studies. The results of these case studies are novel or validating results already reported in the literature, indicating the accuracy of HebbPlot. Our results indicate that promoters have a directional chromatin signature; several marks tend to stretch downstream or upstream. H3K4me3 and H3K79me2 have clear directional distributions around active promoters. In addition, the signatures of high- and low-CpG promoters are different; H3K4me3, H3K9ac, and H3K27ac are the most different marks. When we studied the signatures of enhancers active in eight tissues, we observed that these signatures are similar, but not identical. Further, we identified some histone modifications - H3K36me3, H3K79me1, H3K79me2, and H4K8ac - that are associated with coding regions of active genes. Other marks - H4K12ac, H3K14ac, H3K27me3, and H2AK5ac - were found to be weakly associated with coding regions of inactive genes.

CONCLUSIONS

This study resulted in a novel software tool, HebbPlot, for learning and visualizing the chromatin signature of a genetic element. Using HebbPlot, we produced a visual catalog of the signatures of multiple genetic elements in 57 cell types available through the Roadmap Epigenomics Project. Furthermore, we made a progress toward a functional catalog consisting of 22 histone marks. In sum, HebbPlot is applicable to a wide array of studies, facilitating the deciphering of the histone code.

Molecular mechanism(s) involved in differential expression of vitamin C transporters along the intestinal tract

Mammalian cells utilize two transporters for the uptake of ascorbic acid (AA), Na⁺-dependent vitamin C transporter SVCT-1 and SVCT-2. In the intestine, these transporters are involved in AA absorption and are expressed at the apical and basolateral membrane domains of the polarized epithelia, respectively. Little is known about the differential expression of these two transporters along the anterior-posterior axis of the intestinal tract and the molecular mechanism(s) that dictate this pattern of expression. We used mouse and human intestinal cDNAs to address these issues. The results showed a significantly lower rate of carrier-mediated AA uptake by mouse colon than jejunum. This was associated with a significantly lower level of expression of SVCT-1 and SVCT-2 at the protein, mRNA, and heterogeneous nuclear RNA (hnRNA) levels in the colon than the jejunum, implying the involvement of transcriptional mechanism(s). Similarly, expression levels of SVCT-1 and SVCT-2 mRNA and hnRNA were significantly lower in human colon. We also examined the levels of expression of hepatocyte nuclear factor 1α and specificity protein 1, which drive transcription of the Slc23a1 and Slc23a2 promoters, respectively, and found them to be markedly lower in the colon. Furthermore, significantly lower levels of the activating markers for histone (H3) modifications [H3 trimethylation of lysine 4 (H3K4me3) and H3 triacetylation of lysine 9 (H3K9ac)] were observed in the Slc23a1 and Slc23a2 promoters in the colon. These findings show, for the first time, that SVCT-1 and SVCT-2 are differentially expressed along the intestinal tract and that this pattern of expression is, at least in part, mediated via transcriptional/epigenetic mechanisms.NEW & NOTEWORTHY Our findings show, for the first time, that transporters of the water-soluble vitamin ascorbic acid (i.e., the vitamin C transporters SVCT-1 and SVCT-2) are differentially expressed along the length of the intestinal tract and that the pattern of expression is mediated, at least in part, by transcriptional and epigenetic mechanism(s) affecting both Slc23a1 and Slc23a2 genes.

Striatal H3K27 Acetylation Linked to Glutamatergic Gene Dysregulation in Human Heroin Abusers Holds Promise as Therapeutic Target

BACKGROUND

Opiate abuse and overdose reached epidemic levels in the United States. However, despite significant advances in animal and in vitro models, little knowledge has been directly accrued regarding the neurobiology of the opiate-addicted human brain.

METHODS

We used postmortem human brain specimens from a homogeneous European Caucasian population of heroin users for transcriptional and epigenetic profiling, as well as direct assessment of chromatin accessibility in the striatum, a brain region central to reward and emotion. A rat heroin self-administration model was used to obtain translational molecular and behavioral insights.

RESULTS

Our transcriptome approach revealed marked impairments related to glutamatergic neurotransmission and chromatin remodeling in the human striatum. A series of biochemical experiments tracked the specific location of the epigenetic disturbances to hyperacetylation of lysine 27 of histone H3, showing dynamic correlations with heroin use history and acute opiate toxicology. Targeted investigation of GRIA1, a glutamatergic gene implicated in drug-seeking behavior, verified the increased enrichment of lysine-27 acetylated histone H3 at discrete loci, accompanied by enhanced chromatin accessibility at hyperacetylated regions in the gene body. Analogous epigenetic impairments were detected in the striatum of heroin self-administering rats. Using this translational model, we showed that bromodomain inhibitor JQ1, which blocks the functional readout of acetylated lysines, reduced heroin self-administration and cue-induced drug-seeking behavior.

CONCLUSIONS

Overall, our data suggest that heroin-related histone H3 hyperacetylation contributes to glutamatergic transcriptional changes that underlie addiction behavior and identify JQ1 as a promising candidate for targeted clinical interventions in heroin use disorder.

MicroRNA-Based Linkage between Aging and Cancer: from Epigenetics View Point

Ageing is a complex process and a broad spectrum of physical, psychological, and social changes over time. Accompanying diseases and disabilities, which can interfere with cancer treatment and recovery, occur in old ages. MicroRNAs (miRNAs) are a set of small non-coding RNAs, which have considerable roles in post-transcriptional regulation at gene expression level. In this review, we attempted to summarize the current knowledge of miRNAs functions in ageing, with mainly focuses on malignancies and all underlying genetic, molecular and epigenetics mechanisms. The evidences indicated the complex and dynamic nature of miRNA-based linkage of ageing and cancer at genomics and epigenomics levels which might be generally crucial for understanding the mechanisms of age-related cancer and ageing. Recently in the field of cancer and ageing, scientists claimed that uric acid can be used to regulate reactive oxygen species (ROS), leading to cancer and ageing prevention; these findings highlight the role of miRNA-based inhibition of the SLC2A9 antioxidant pathway in cancer, as a novel way to kill malignant cells, while a patient is fighting with cancer.

The upstreams and downstreams of H3K79 methylation by DOT1L

Histone modifications regulate key processes of eukaryotic genomes. Misregulation of the enzymes that place these modifications can lead to disease. An example of this is DOT1L, the enzyme that can mono-, di-, and trimethylate the nucleosome core on lysine 79 of histone H3 (H3K79). DOT1L plays a role in development and its misregulation has been implicated in several cancers, most notably leukemias caused by a rearrangement of the MLL gene. A DOT1L inhibitor is in clinical trials for these leukemias and shows promising results, yet we are only beginning to understand DOT1L's function and regulation in the cell. Here, we review what happens upstream and downstream of H3K79 methylation. H3K79 methylation levels are highest in transcribed genes, where H2B ubiquitination can promote DOT1L activity. In addition, DOT1L can be targeted to transcribed regions of the genome by several of its interaction partners. Although methylation levels strongly correlate with transcription, the mechanistic link between the two is unclear and probably context-dependent. Methylation of H3K79 may act through recruiting or repelling effector proteins, but we do not yet know which effectors mediate DOT1L's functions. Understanding DOT1L biology better will help us to understand the effects of DOT1L inhibitors and may allow the development of alternative strategies to target the DOT1L pathway.

Mesodermal iPSC-derived progenitor cells functionally regenerate cardiac and skeletal muscle

Conditions such as muscular dystrophies (MDs) that affect both cardiac and skeletal muscles would benefit from therapeutic strategies that enable regeneration of both of these striated muscle types. Protocols have been developed to promote induced pluripotent stem cells (iPSCs) to differentiate toward cardiac or skeletal muscle; however, there are currently no strategies to simultaneously target both muscle types. Tissues exhibit specific epigenetic alterations; therefore, source-related lineage biases have the potential to improve iPSC-driven multilineage differentiation. Here, we determined that differential myogenic propensity influences the commitment of isogenic iPSCs and a specifically isolated pool of mesodermal iPSC-derived progenitors (MiPs) toward the striated muscle lineages. Differential myogenic propensity did not influence pluripotency, but did selectively enhance chimerism of MiP-derived tissue in both fetal and adult skeletal muscle. When injected into dystrophic mice, MiPs engrafted and repaired both skeletal and cardiac muscle, reducing functional defects. Similarly, engraftment into dystrophic mice of canine MiPs from dystrophic dogs that had undergone TALEN-mediated correction of the MD-associated mutation also resulted in functional striatal muscle regeneration. Moreover, human MiPs exhibited the same capacity for the dual differentiation observed in murine and canine MiPs. The findings of this study suggest that MiPs should be further explored for combined therapy of cardiac and skeletal muscles.

Genome-wide analysis of DNA methylation and gene expression patterns in purified, uncultured human liver cells and activated hepatic stellate cells

BACKGROUND & AIMS

Liver fibrogenesis - scarring of the liver that can lead to cirrhosis and liver cancer - is characterized by hepatocyte impairment, capillarization of liver sinusoidal endothelial cells (LSECs) and hepatic stellate cell (HSC) activation. To date, the molecular determinants of a healthy human liver cell phenotype remain largely uncharacterized. Here, we assess the transcriptome and the genome-wide promoter methylome specific for purified, non-cultured human hepatocytes, LSECs and HSCs, and investigate the nature of epigenetic changes accompanying transcriptional changes associated with activation of HSCs.

MATERIAL AND METHODS

Gene expression profile and promoter methylome of purified, uncultured human liver cells and culture-activated HSCs were respectively determined using Affymetrix HG-U219 genechips and by methylated DNA immunoprecipitation coupled to promoter array hybridization. Histone modification patterns were assessed at the single-gene level by chromatin immunoprecipitation and quantitative PCR.

RESULTS

We unveil a DNA-methylation-based epigenetic relationship between hepatocytes, LSECs and HSCs despite their distinct ontogeny. We show that liver cell type-specific DNA methylation targets early developmental and differentiation-associated functions. Integrative analysis of promoter methylome and transcriptome reveals partial concordance between DNA methylation and transcriptional changes associated with human HSC activation. Further, we identify concordant histone methylation and acetylation changes in the promoter and putative novel enhancer elements of genes involved in liver fibrosis.

CONCLUSIONS

Our study provides the first epigenetic blueprint of three distinct freshly isolated, human hepatic cell types and of epigenetic changes elicited upon HSC activation.

Association analysis between the distributions of histone modifications and gene expression in the human embryonic stem cell

It is well known that histone modifications are associated with gene expression. In order to further study this relationship, 16 kinds of Chip-seq histone modification data and mRNA-seq data of the human embryonic stem cell H1 are chosen. The distributions of histone modifications in the regions flanking transcription start sites (TSSs) for highly expressed and lowly expressed genes are computed, respectively. And four types of distributions of histone modifications in regions flanking TSSs and the spatial patterning of the correlations between histone modifications and gene expression are detected. Our results suggest that the correlations between the regions overlapped by peaks are higher than the non-overlapped ones for each histone modification. In addition, to obtain the effect of the cooperative action of histone modification on gene expression, five histone modification clusters are found in highly expressed and lowly expressed genes, histone modification and gene expression interaction network is constructed. To further explore which region is the main target region for the specific histone modification, the human genes are divided into five functional regions. The results indicate that histone modifications are mostly located in the promoters of highly expressed genes versus the exons of lowly expressed genes, and exons have a smaller range of normalized tag counts than other gene elements in the two groups of genes. Finally, the type specificity and regional bias of histone modifications for 11 key transcription factor genes regulating the stem cell renewal are analyzed.

Histone methyltransferase Setd8 represses Gata2 expression and regulates erythroid maturation

Setd8 is the sole histone methyltransferase in mammals capable of monomethylating histone H4 lysine 20 (H4K20me1). Setd8 is expressed at significantly higher levels in erythroid cells than any other cell or tissue type, suggesting that Setd8 has an erythroid-cell-specific function. To test this hypothesis, stable Setd8 knockdown was established in extensively self-renewing erythroblasts (ESREs), a well-characterized, nontransformed model of erythroid maturation. Knockdown of Setd8 resulted in impaired erythroid maturation characterized by a delay in hemoglobin accumulation, larger mean cell area, persistent ckit expression, incomplete nuclear condensation, and lower rates of enucleation. Setd8 knockdown did not alter ESRE proliferation or viability or result in accumulation of DNA damage. Global gene expression analyses following Setd8 knockdown demonstrated that in erythroid cells, Setd8 functions primarily as a repressor. Most notably, Gata2 expression was significantly higher in knockdown cells than in control cells and Gata2 knockdown rescued some of the maturation impairments associated with Setd8 disruption. Setd8 occupies critical regulatory elements in the Gata2 locus, and knockdown of Setd8 resulted in loss of H4K20me1 and gain of H4 acetylation at the Gata2 1S promoter. These results suggest that Setd8 is an important regulator of erythroid maturation that works in part through repression of Gata2 expression.

A liver-enriched long non-coding RNA, lncLSTR, regulates systemic lipid metabolism in mice

Long non-coding RNAs (lncRNAs) constitute a significant portion of mammalian genome, yet the physiological importance of lncRNAs is largely unknown. Here, we identify a liver-enriched lncRNA in mouse that we term liver-specific triglyceride regulator (lncLSTR). Mice with a liver-specific depletion of lncLSTR exhibit a marked reduction in plasma triglyceride levels. We show that lncLSTR depletion enhances apoC2 expression, leading to robust lipoprotein lipase activation and increased plasma triglyceride clearance. We further demonstrate that the regulation of apoC2 expression occurs through an FXR-mediated pathway. LncLSTR forms a molecular complex with TDP-43 to regulate expression of Cyp8b1, a key enzyme in the bile acid synthesis pathway, and engenders an in vivo bile pool that induces apoC2 expression through FXR. Finally, we demonstrate that lncLSTR depletion can reduce triglyceride levels in a hyperlipidemia mouse model. Taken together, these data support a model in which lncLSTR regulates a TDP-43/FXR/apoC2-dependent pathway to maintain systemic lipid homeostasis.

Modelling epigenetic regulation of gene expression in 12 human cell types reveals combinatorial patterns of cell-type-specific genes

The maintenance of the diverse cell types in a multicellular organism is one of the fundamental mysteries of biology. Modelling the dynamic regulatory relationships between the histone modifications and the gene expression across the diverse cell types is essential for the authors to understand the mechanisms of the epigenetic regulation. Here, the authors thoroughly assessed the histone modification enrichment profiles at the promoters and constructed quantitative models between the histone modification abundances and the gene expression in 12 human cell types. The author's results showed that the histone modifications at the promoters exhibited remarkably cell-type-dependent variability in the cell-type-specific (CTS) genes. They demonstrated that the variable profiles of the modifications are highly predictive for the dynamic changes of the gene expression across all the cell types. Their findings revealed the close relationship between the combinatorial patterns of the histone modifications and the CTS gene expression. They anticipate that the findings and the methods they used in this study could provide useful information for the future studies of the regulatory roles of the histone modifications in the CTS genes.

Epigenetic modulation of intestinal cholesterol transporter Niemann-Pick C1-like 1 (NPC1L1) gene expression by DNA methylation

Intestinal NPC1L1 transporter is essential for cholesterol absorption and the maintenance of cholesterol homeostasis in the body. NPC1L1 is differentially expressed along the gastrointestinal tract with very low levels in the colon as compared with the small intestine. This study was undertaken to examine whether DNA methylation was responsible for segment-specific expression of NPC1L1. Treatment of mice with 5-azacytidine (i.p.) resulted in a significant dose-dependent increase in NPC1L1 mRNA expression in the colon. The lack of expression of NPC1L1 in the normal colon was associated with high levels of methylation in the area flanking the 3-kb fragment upstream of the initiation site of the mouse NPC1L1 gene in mouse colon as analyzed by EpiTYPER® MassARRAY®. The high level of methylation in the colon was observed in specific CpG dinucleotides and was significantly decreased in response to 5-azacytidine. Similar to mouse NPC1L1, 5-azacytidine treatment also increased the level of human NPC1L1 mRNA expression in the intestinal HuTu-80 cell line in a dose- and time-dependent manner. Silencing the expression of DNA methyltransferase DNMT1, -2, -3A, and -3B alone by siRNA did not affect NPC1L1 expression in HuTu-80 cells. However, the simultaneous attenuation of DNMT1 and -3B expression caused a significant increase in NPC1L1 mRNA expression as compared with control. Also, in vitro methylation of the human NPC1L1 promoter significantly decreased NPC1L1 promoter activity in human intestinal Caco2 cells. In conclusion, our data demonstrated for the first time that DNA methylation in the promoter region of the NPC1L1 gene appears to be a major mechanism underlying differential expression of NPC1L1 along the length of the gastrointestinal tract.

Regression hidden Markov modeling reveals heterogeneous gene expression regulation: a case study in mouse embryonic stem cells

BACKGROUND

Studies have shown the strong association between histone modification levels and gene expression levels. The detailed relationships between the two can vary substantially due to differential regulation, and hence a simple regression model may not be adequate. We apply a regression hidden Markov model (regHMM) to further investigate the potential multiple relationships between genes and histone methylation levels in mouse embryonic stem cells.

RESULTS

Seven histone methylation levels are used in the study. Averaged histone modifications over non-overlapping 200 bp windows on the range transcription starting site (TSS) ± 1 Kb are used as predictors, and in total 70 explanatory variables are generated. Based on regHMM results, genes segregated into two groups, referred to as State 1 and State 2, have distinct association strengths. Genes in State 1 are better explained by histone methylation levels with R(2)=.72 while those in State 2 have weaker association strength with R(2)=.38. The regression coefficients in the two states are not very different in magnitude except in the intercept,.25 and 1.15 for State 1 and State 2, respectively. We found specific GO categories that may be attributed to the different relationships. The GO categories more frequently observed in State 2 match those of housekeeping genes, such as cytoplasm, nucleus, and protein binding. In addition, the housekeeping gene expression levels are significantly less explained by histone methylation in mouse embryonic stem cells, which is consistent with the constitutive expression patterns that would be expected.

CONCLUSION

Gene expression levels are not universally affected by histone methylation levels, and the relationships between the two differ by the gene functions. The expression levels of the genes that perform the most common housekeeping genes' GO categories are less strongly associated with histone methylation levels. We suspect that additional biological factors may also be strongly associated with the gene expression levels in State 2. We discover that the effect of the presence of CpG island in TSS ± 1 Kb is larger in State 2.

ngs.plot: Quick mining and visualization of next-generation sequencing data by integrating genomic databases

BACKGROUND

Understanding the relationship between the millions of functional DNA elements and their protein regulators, and how they work in conjunction to manifest diverse phenotypes, is key to advancing our understanding of the mammalian genome. Next-generation sequencing technology is now used widely to probe these protein-DNA interactions and to profile gene expression at a genome-wide scale. As the cost of DNA sequencing continues to fall, the interpretation of the ever increasing amount of data generated represents a considerable challenge.

RESULTS

We have developed ngs.plot - a standalone program to visualize enrichment patterns of DNA-interacting proteins at functionally important regions based on next-generation sequencing data. We demonstrate that ngs.plot is not only efficient but also scalable. We use a few examples to demonstrate that ngs.plot is easy to use and yet very powerful to generate figures that are publication ready.

CONCLUSIONS

We conclude that ngs.plot is a useful tool to help fill the gap between massive datasets and genomic information in this era of big sequencing data.

Transcription factor binding sites prediction based on modified nucleosomes

In computational methods, position weight matrices (PWMs) are commonly applied for transcription factor binding site (TFBS) prediction. Although these matrices are more accurate than simple consensus sequences to predict actual binding sites, they usually produce a large number of false positive (FP) predictions and so are impoverished sources of information. Several studies have employed additional sources of information such as sequence conservation or the vicinity to transcription start sites to distinguish true binding regions from random ones. Recently, the spatial distribution of modified nucleosomes has been shown to be associated with different promoter architectures. These aligned patterns can facilitate DNA accessibility for transcription factors. We hypothesize that using data from these aligned and periodic patterns can improve the performance of binding region prediction. In this study, we propose two effective features, “modified nucleosomes neighboring” and “modified nucleosomes occupancy”, to decrease FP in binding site discovery. Based on these features, we designed a logistic regression classifier which estimates the probability of a region as a TFBS. Our model learned each feature based on Sp1 binding sites on Chromosome 1 and was tested on the other chromosomes in human CD4+T cells. In this work, we investigated 21 histone modifications and found that only 8 out of 21 marks are strongly correlated with transcription factor binding regions. To prove that these features are not specific to Sp1, we combined the logistic regression classifier with the PWM, and created a new model to search TFBSs on the genome. We tested the model using transcription factors MAZ, PU.1 and ELF1 and compared the results to those using only the PWM. The results show that our model can predict Transcription factor binding regions more successfully. The relative simplicity of the model and capability of integrating other features make it a superior method for TFBS prediction.

Epigenetic regulation of miR-21 in colorectal cancer: ITGB4 as a novel miR-21 target and a three-gene network (miR-21-ITGΒ4-PDCD4) as predictor of metastatic tumor potential

Previous studies have uncovered several transcription factors that determine biological alterations in tumor cells to execute the invasion-metastasis cascade, including the epithelial-mesenchymal transition (EMT). We sought to investigate the role of miR-21 in colorectal cancer regulation. For this purpose, miR-21 expression was quantified in a panel of colorectal cancer cell lines and clinical specimens. High expression was found in cell lines with EMT properties and in the vast majority of human tumor specimens. We demonstrate in a cell-specific manner the occupancy of MIR-21 gene promoter by AP-1 and ETS1 transcription factors and, for the first time, the pattern of histone posttranslational modifications necessary for miR-21 overexpression. We also show that Integrin-β4 (ITGβ4), exclusively expressed in polarized epithelial cells, is a novel miR-21 target gene and plays a role in the regulation of EMT, since it is remarkably de-repressed after transient miR-21 silencing and downregulated after miR-21 overexpression. miR-21-dependent change of ITGβ4 expression significantly affects cell migration properties of colon cancer cells. Finally, in a subgroup of tumor specimens, ROC curve analysis performed on quantitative PCR data sets for miR-21, ITGβ4, and PDCD4 shows that the combination of high miR-21 with low ITGβ4 and PDCD4 expression is able to predict presence of metastasis. In conclusion, miR-21 is a key player in oncogenic EMT, its overexpression is controlled by the cooperation of genetic and epigenetic alterations, and its levels, along with ITGβ4 and PDCD4 expression, could be exploited as a prognostic tool for CRC metastasis.

Finding associations among histone modifications using sparse partial correlation networks

Histone modifications are known to play an important role in the regulation of transcription. While individual modifications have received much attention in genome-wide analyses, little is known about their relationships. Some authors have built Bayesian networks of modifications, however most often they have used discretized data, and relied on unrealistic assumptions such as the absence of feedback mechanisms or hidden confounding factors. Here, we propose to infer undirected networks based on partial correlations between histone modifications. Within the partial correlation framework, correlations among two variables are controlled for associations induced by the other variables. Partial correlation networks thus focus on direct associations of histone modifications. We apply this methodology to data in CD4+ cells. The resulting network is well supported by common knowledge. When pairs of modifications show a large difference between their correlation and their partial correlation, a potential confounding factor is identified and provided as explanation. Data from different cell types (IMR90, H1) is also exploited in the analysis to assess the stability of the networks. The results are remarkably similar across cell types. Based on this observation, the networks from the three cell types are integrated into a consensus network to increase robustness. The data and the results discussed in the manuscript can be found, together with code, on http://spcn.molgen.mpg.de/index.html.

Predictive models of gene regulation from high-throughput epigenomics data

The epigenetic regulation of gene expression involves multiple factors. The synergistic or antagonistic action of these factors has suggested the existence of an epigenetic code for gene regulation. Highthroughput sequencing (HTS) provides an opportunity to explore this code and to build quantitative models of gene regulation based on epigenetic differences between specific cellular conditions. We describe a new computational framework that facilitates the systematic integration of HTS epigenetic data. Our method relates epigenetic signals to expression by comparing two conditions. We show its effectiveness by building a model that predicts with high accuracy significant expression differences between two cell lines, using epigenetic data from the ENCODE project. Our analyses provide evidence for a degenerate epigenetic code, which involves multiple genic regions. In particular, signal changes at the 1st exon, 1st intron, and downstream of the polyadenylation site are found to associate strongly with expression regulation. Our analyses also show a different epigenetic code for intron-less and intron-containing genes. Our work provides a general methodology to do integrative analysis of epigenetic differences between cellular conditions that can be applied to other studies, like cell differentiation or carcinogenesis.

Histone methylation analysis and pathway predictions in chickens after MDV infection

Marek's disease (MD) is a lymphoproliferative disease in chicken induced by Marek's disease virus (MDV). Although studies have focused on the genetic differences between the resistant and susceptible chicken, less is known about the role of epigenetic factors in MD. In this study, genome-wide histone modifications in the non-MHC-associated resistant and susceptible chicken lines were examined. We found that tri-methylation at histone H3 Lys4 (H3K4me3) enrichment is positively correlated with the expression of protein coding genes as well as microRNA (miRNA) genes, whereas tri-methylation at histone H3 Lys27 (H3K27me3) exhibits a negative correlation. By identifying line-specific histone modifications in MDV infection, we found unique H3K4me3 islands in the resistant chicken activated genes, which are related to immune response and cell adhesion. Interestingly, we also found some miRNAs from unique H3K27me3 patterns in the susceptible chickens that targeted genes involved in 5-hydroxytryptamine (5-HT)-receptor and adrenergic receptor pathways. In conclusion, dynamic line-specific histone modifications in response to MDV infection suggested that intrinsic epigenetic mechanisms may play a role in MD-resistance and -susceptibility.

Bivalent-like chromatin markers are predictive for transcription start site distribution in human

Deep sequencing of 5′ capped transcripts has revealed a variety of transcription initiation patterns, from narrow, focused promoters to wide, broad promoters. Attempts have already been made to model empirically classified patterns, but virtually no quantitative models for transcription initiation have been reported. Even though both genetic and epigenetic elements have been associated with such patterns, the organization of regulatory elements is largely unknown. Here, linear regression models were derived from a pool of regulatory elements, including genomic DNA features, nucleosome organization, and histone modifications, to predict the distribution of transcription start sites (TSS). Importantly, models including both active and repressive histone modification markers, e.g. H3K4me3 and H4K20me1, were consistently found to be much more predictive than models with only single-type histone modification markers, indicating the possibility of “bivalent-like” epigenetic control of transcription initiation. The nucleosome positions are proposed to be coded in the active component of such bivalent-like histone modification markers. Finally, we demonstrated that models trained on one cell type could successfully predict TSS distribution in other cell types, suggesting that these models may have a broader application range.

Computational inference of mRNA stability from histone modification and transcriptome profiles

Histone modifications play important roles in regulating eukaryotic gene expression and have been used to model expression levels. Here, we present a regression model to systematically infer mRNA stability by comparing transcriptome profiles with ChIP-seq of H3K4me3, H3K27me3 and H3K36me3. The results from multiple human and mouse cell lines show that the inferred unstable mRNAs have significantly longer 3′Untranslated Regions (UTRs) and more microRNA binding sites within 3′UTR than the inferred stable mRNAs. Regression residuals derived from RNA-seq, but not from GRO-seq, are highly correlated with the half-lives measured by pulse-labeling experiments, supporting the rationale of our inference. Whereas, the functions enriched in the inferred stable and unstable mRNAs are consistent with those from pulse-labeling experiments, we found the unstable mRNAs have higher cell-type specificity under functional constraint. We conclude that the systematical use of histone modifications can differentiate non-expressed mRNAs from unstable mRNAs, and distinguish stable mRNAs from highly expressed ones. In summary, we represent the first computational model of mRNA stability inference that compares transcriptome and epigenome profiles, and provides an alternative strategy for directing experimental measurements.