EpiFactors: a comprehensive database of human epigenetic factors and complexes

Linking aberrant chromatin features in chronic lymphocytic leukemia to transcription factor networks

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

KMT1 family methyltransferases regulate heterochromatin-nuclear periphery tethering via histone and non-histone protein methylation

Euchromatic histone methyltransferases (EHMTs), members of the KMT1 family, methylate histone and non-histone proteins. Here, we uncover a novel role for EHMTs in regulating heterochromatin anchorage to the nuclear periphery (NP) via non-histone methylation. We show that EHMTs methylate and stabilize LaminB1 (LMNB1), which associates with the H3K9me2-marked peripheral heterochromatin. Loss of LMNB1 methylation or EHMTs abrogates heterochromatin anchorage at the NP We further demonstrate that the loss of EHMTs induces many hallmarks of aging including global reduction of H3K27methyl marks and altered nuclear morphology. Consistent with this, we observe a gradual depletion of EHMTs, which correlates with loss of methylated LMNB1 and peripheral heterochromatin in aging human fibroblasts. Restoration of EHMT expression reverts peripheral heterochromatin defects in aged cells. Collectively, our work elucidates a new mechanism by which EHMTs regulate heterochromatin domain organization and reveals their impact on fundamental changes associated with the intrinsic aging process.

Coexpression patterns define epigenetic regulators associated with neurological dysfunction

Coding variants in epigenetic regulators are emerging as causes of neurological dysfunction and cancer. However, a comprehensive effort to identify disease candidates within the human epigenetic machinery (EM) has not been performed; it is unclear whether features exist that distinguish between variation-intolerant and variation-tolerant EM genes, and between EM genes associated with neurological dysfunction versus cancer. Here, we rigorously define 295 genes with a direct role in epigenetic regulation (writers, erasers, remodelers, readers). Systematic exploration of these genes reveals that although individual enzymatic functions are always mutually exclusive, readers often also exhibit enzymatic activity (dual-function EM genes). We find that the majority of EM genes are very intolerant to loss-of-function variation, even when compared to the dosage sensitive transcription factors, and we identify 102 novel EM disease candidates. We show that this variation intolerance is driven by the protein domains encoding the epigenetic function, suggesting that disease is caused by a perturbed chromatin state. We then describe a large subset of EM genes that are coexpressed within multiple tissues. This subset is almost exclusively populated by extremely variation-intolerant genes and shows enrichment for dual-function EM genes. It is also highly enriched for genes associated with neurological dysfunction, even when accounting for dosage sensitivity, but not for cancer-associated EM genes. Finally, we show that regulatory regions near epigenetic regulators are genetically important for common neurological traits. These findings prioritize novel disease candidate EM genes and suggest that this coexpression plays a functional role in normal neurological homeostasis.

Prenatal arsenic exposure alters the placental expression of multiple epigenetic regulators in a sex-dependent manner

BACKGROUND

Prenatal exposure to arsenic has been linked to a range of adverse health conditions in later life. Such fetal origins of disease are frequently the result of environmental effects on the epigenome, leading to long-term alterations in gene expression. Several studies have demonstrated effects of prenatal arsenic exposure on DNA methylation; however the impact of arsenic on the generation and decoding of post-translational histone modifications (PTHMs) is less well characterized, and has not been studied in the context of prenatal human exposures.

METHODS

In the current study, we examined the effect of exposure to low-to-moderate levels of arsenic in a US birth cohort, on the expression of 138 genes encoding key epigenetic regulators in the fetal portion of the placenta. Our candidate genes included readers, writers and erasers of PTHMs, and chromatin remodelers.

RESULTS

Arsenic exposure was associated with the expression of 27 of the 138 epigenetic genes analyzed. When the cohort was stratified by fetal sex, arsenic exposure was associated with the expression of 40 genes in male fetal placenta, and only 3 non-overlapping genes in female fetal placenta. In particular, we identified an inverse relationship between arsenic exposure and expression of the gene encoding the histone methyltransferase, PRDM6 (p < 0.001). Mutation of PRDM6 has been linked to the congenital heart defect, patent ductus arteriosus.

CONCLUSIONS

Our findings suggest that prenatal arsenic exposure may have sex-specific effects on the fetal epigenome, which could plausibly contribute to its subsequent health impacts.

CpG traffic lights are markers of regulatory regions in human genome

Background

DNA methylation is involved in the regulation of gene expression. Although bisulfite-sequencing based methods profile DNA methylation at a single CpG resolution, methylation levels are usually averaged over genomic regions in the downstream bioinformatic analysis.

Results

We demonstrate that on the genome level a single CpG methylation can serve as a more accurate predictor of gene expression than an average promoter / gene body methylation. We define CpG traffic lights (CpG TL) as CpG dinucleotides with a significant correlation between methylation and expression of a gene nearby. CpG TL are enriched in all regulatory regions. Among all promoters, CpG TL are especially enriched in poised ones, suggesting involvement of DNA methylation in their regulation. Yet, binding of only a handful of transcription factors, such as NRF1, ETS, STAT and IRF-family members, could be regulated by direct methylation of transcription factor binding sites (TFBS) or its close proximity. For the majority of TF, an alternative scenario is more likely: methylation and inactivation of the whole regulatory element indirectly represses functional TF binding with a CpG TL being a reliable marker of such inactivation.

Conclusions

CpG TL provide a promising insight into mechanisms of enhancer activity and gene regulation linking methylation of single CpG to gene expression. CpG TL methylation can be used as reliable markers of enhancer activity and gene expression in applications, e.g. in clinic where measuring DNA methylation is easier compared to directly measuring gene expression due to more stable nature of DNA.

Electronic supplementary material

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

TBL1 is required for the mesenchymal phenotype of transformed breast cancer cells

The epithelial-to-mesenchymal transition (EMT) and its reversion (MET) are related to tumor cell dissemination and migration, tumor circulating cell generation, cancer stem cells, chemoresistance, and metastasis formation. To identify chromatin and epigenetic factors possibly involved in the process of EMT, we compare the levels of expression of epigenetic genes in a transformed human breast epithelial cell line (HMEC-RAS) versus a stable clone of the same cell line expressing the EMT master regulator ZEB1 (HMEC-RAS-ZEB1). One of the factors strongly induced in the HMEC-RAS-ZEB1 cells was Transducin beta-like 1 (TBL1), a component of the NCoR complex, which has both corepressor and coactivator activities. We show that TBL1 interacts with ZEB1 and that both factors cooperate to repress the promoter of the epithelial gene E-cadherin (CDH1) and to autoactivate the ZEB1 promoter. Consistent with its central role, TBL1 is required for mesenchymal phenotypes of transformed breast epithelial and breast cancer cell lines of the claudin-low subtype. Importantly, a high expression of the TBL1 gene correlates with poor prognosis and increased proportion of metastasis in breast cancer patients, indicating that the level of TBL1 expression can be used as a prognostic marker.

FusionGDB: fusion gene annotation DataBase

Gene fusion is one of the hallmarks of cancer genome via chromosomal rearrangement initiated by DNA double-strand breakage. To date, many fusion genes (FGs) have been established as important biomarkers and therapeutic targets in multiple cancer types. To better understand the function of FGs in cancer types and to promote the discovery of clinically relevant FGs, we built FusionGDB (Fusion Gene annotation DataBase) available at https://ccsm.uth.edu/FusionGDB. We collected 48 117 FGs across pan-cancer from three representative fusion gene resources: the improved database of chimeric transcripts and RNA-seq data (ChiTaRS 3.1), an integrative resource for cancer-associated transcript fusions (TumorFusions), and The Cancer Genome Atlas (TCGA) fusions by Gao et al. For these ∼48K FGs, we performed functional annotations including gene assessment across pan-cancer fusion genes, open reading frame (ORF) assignment, and retention search of 39 protein features based on gene structures of multiple isoforms with different breakpoints. We also provided the fusion transcript and amino acid sequences according to multiple breakpoints and transcript isoforms. Our analyses identified 331, 303 and 667 in-frame FGs with retaining kinase, DNA-binding, and epigenetic factor domains, respectively, as well as 976 FGs lost protein-protein interaction. FusionGDB provides six categories of annotations: FusionGeneSummary, FusionProtFeature, FusionGeneSequence, FusionGenePPI, RelatedDrug and RelatedDisease.

Transposable elements are regulated by context-specific patterns of chromatin marks in mouse embryonic stem cells

The majority of mammalian genomes are devoted to transposable elements (TEs). Whilst TEs are increasingly recognized for their important biological functions, they are a potential danger to genomic stability and are carefully regulated by the epigenetic system. However, the full complexity of this regulatory system is not understood. Here, using mouse embryonic stem cells, we show that TEs are suppressed by heterochromatic marks like H3K9me3, and are also labelled by all major types of chromatin modification in complex patterns, including bivalent activatory and repressive marks. We identified 29 epigenetic modifiers that significantly deregulated at least one type of TE. The loss of Setdb1, Ncor2, Rnf2, Kat5, Prmt5, Uhrf1, and Rrp8 caused widespread changes in TE expression and chromatin accessibility. These effects were context-specific, with different chromatin modifiers regulating the expression and chromatin accessibility of specific subsets of TEs. Our work reveals the complex patterns of epigenetic regulation of TEs.

Transposable elements (TEs) fulfill essential but poorly understood roles in genome organization and gene expression control. Here the authors show that the regulation of TEs occurs through overlapping epigenetic mechanisms that control the expression and chromatin signatures at TEs.

Epigenetic modifications acetylation and deacetylation play important roles in juvenile hormone action

BACKGROUND

Epigenetic modifications including DNA methylation and post-translational modifications of histones are known to regulate gene expression. Antagonistic activities of histone acetyltransferases (HATs) and histone deacetylases (HDACs) mediate transcriptional reprogramming during insect development as shown in Drosophila melanogaster and other insects. Juvenile hormones (JH) play vital roles in the regulation of growth, development, metamorphosis, reproduction and other physiological processes. However, our current understanding of epigenetic regulation of JH action is still limited. Hence, we studied the role of CREB binding protein (CBP, contains HAT domain) and Trichostatin A (TSA, HDAC inhibitor) on JH action.

RESULTS

Exposure of Tribolium castaneum cells (TcA cells) to JH or TSA caused an increase in expression of Kr-h1 (a known JH-response gene) and 31 or 698 other genes respectively. Knockdown of the gene coding for CBP caused a decrease in the expression of 456 genes including Kr-h1. Interestingly, the expression of several genes coding for transcription factors, nuclear receptors, P450 and fatty acid synthase family members that are known to mediate JH action were affected by CBP knockdown or TSA treatment.

CONCLUSIONS

These data suggest that acetylation and deacetylation mediated by HATs and HDACs play an important role in JH action.

The role of TCF3 as potential master regulator in blastemal Wilms tumors

Wilms tumors are the most common type of pediatric kidney tumors. While the overall prognosis for patients is favorable, especially tumors that exhibit a blastemal subtype after preoperative chemotherapy have a poor prognosis. For an improved risk assessment and therapy stratification, it is essential to identify the driving factors that are distinctive for this aggressive subtype. In our study, we compared gene expression profiles of 33 tumor biopsies (17 blastemal and 16 other tumors) after neoadjuvant chemotherapy. The analysis of this dataset using the Regulator Gene Association Enrichment algorithm successfully identified several biomarkers and associated molecular mechanisms that distinguish between blastemal and nonblastemal Wilms tumors. Specifically, regulators involved in embryonic development and epigenetic processes like chromatin remodeling and histone modification play an essential role in blastemal tumors. In this context, we especially identified TCF3 as the central regulatory element. Furthermore, the comparison of ChIP-Seq data of Wilms tumor cell cultures from a blastemal mouse xenograft and a stromal tumor provided further evidence that the chromatin states of blastemal cells share characteristics with embryonic stem cells that are not present in the stromal tumor cell line. These stem-cell like characteristics could potentially add to the increased malignancy and chemoresistance of the blastemal subtype. Along with TCF3, we detected several additional biomarkers that are distinctive for blastemal Wilms tumors after neoadjuvant chemotherapy and that may provide leads for new therapeutic regimens.

A pan-cancer analysis of driver gene mutations, DNA methylation and gene expressions reveals that chromatin remodeling is a major mechanism inducing global changes in cancer epigenomes

Background

Recent large-scale cancer sequencing studies have discovered many novel cancer driver genes (CDGs) in human cancers. Some studies also suggest that CDG mutations contribute to cancer-associated epigenomic and transcriptomic alterations across many cancer types. Here we aim to improve our understanding of the connections between CDG mutations and altered cancer cell epigenomes and transcriptomes on pan-cancer level and how these connections contribute to the known association between epigenome and transcriptome.

Method

Using multi-omics data including somatic mutation, DNA methylation, and gene expression data of 20 cancer types from The Cancer Genome Atlas (TCGA) project, we conducted a pan-cancer analysis to identify CDGs, when mutated, have strong associations with genome-wide methylation or expression changes across cancer types, which we refer as methylation driver genes (MDGs) or expression driver genes (EDGs), respectively.

Results

We identified 32 MDGs, among which, eight are known chromatin modification or remodeling genes. Many of the remaining 24 MDGs are connected to chromatin regulators through either regulating their transcription or physically interacting with them as potential co-factors. We identified 29 EDGs, 26 of which are also MDGs. Further investigation on target genes’ promoters methylation and expression alteration patterns of these 26 overlapping driver genes shows that hyper-methylation of target genes’ promoters are significantly associated with down-regulation of the same target genes and hypo-methylation of target genes’ promoters are significantly associated with up-regulation of the same target genes.

Conclusion

This finding suggests a pivotal role for genetically driven changes in chromatin remodeling in shaping DNA methylation and gene expression patterns during tumor development.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0425-z) contains supplementary material, which is available to authorized users.

Gene Signature of Proliferating Human Pancreatic α Cells

Pancreatic α cells proliferate at a low rate, and little is known about the control of this process. Here we report the characterization of human α cells by large-scale, single-cell RNA sequencing coupled with pseudotime ordering. We identified two large subpopulations and a smaller cluster of proliferating α cells with increased expression of genes involved in cell-cycle regulation. The proliferating α cells were differentiated, had normal levels of GCG expression, and showed no signs of cellular stress. Proliferating α cells were detected in both the G1S and G2M phases of the cell cycle. Human α cells proliferate at a fivefold higher rate than human β cells and express lower levels of the cell-cycle inhibitors CDKN1A and CDKN1C. Collectively, this study provides the gene signatures of human α cells and the genes involved in their cell division. The lower expression of two cell-cycle inhibitors in human α cells could account for their higher rate of proliferation compared with their insulin-producing counterparts.

Chromatin Accessibility Dynamics during iPSC Reprogramming

Cell-fate decisions remain poorly understood at the chromatin level. Here, we map chromatin remodeling dynamics during induction of pluripotent stem cells. ATAC-seq profiling of MEFs expressing Oct4-Sox2-Klf4 (OSK) reveals dynamic changes in chromatin states shifting from open to closed (OC) and closed to open (CO), with an initial burst of OC and an ending surge of CO. The OC loci are largely composed of genes associated with a somatic fate, while the CO loci are associated with pluripotency. Factors/conditions known to impede reprogramming prevent OSK-driven OC and skew OC-CO dynamics. While the CO loci are enriched for OSK motifs, the OC loci are not, suggesting alternative mechanisms for chromatin closing. Sap30, a Sin3A corepressor complex component, is required for the OC shift and facilitates reduced H3K27ac deposition at OC loci. These results reveal a chromatin accessibility logic during reprogramming that may apply to other cell-fate decisions.

The Epigenetic Factor Landscape of Developing Neocortex Is Regulated by Transcription Factors Pax6→ Tbr2→ Tbr1

Epigenetic factors (EFs) regulate multiple aspects of cerebral cortex development, including proliferation, differentiation, laminar fate, and regional identity. The same neurodevelopmental processes are also regulated by transcription factors (TFs), notably the Pax6→ Tbr2→ Tbr1 cascade expressed sequentially in radial glial progenitors (RGPs), intermediate progenitors, and postmitotic projection neurons, respectively. Here, we studied the EF landscape and its regulation in embryonic mouse neocortex. Microarray and in situ hybridization assays revealed that many EF genes are expressed in specific cortical cell types, such as intermediate progenitors, or in rostrocaudal gradients. Furthermore, many EF genes are directly bound and transcriptionally regulated by Pax6, Tbr2, or Tbr1, as determined by chromatin immunoprecipitation-sequencing and gene expression analysis of TF mutant cortices. Our analysis demonstrated that Pax6, Tbr2, and Tbr1 form a direct feedforward genetic cascade, with direct feedback repression. Results also revealed that each TF regulates multiple EF genes that control DNA methylation, histone marks, chromatin remodeling, and non-coding RNA. For example, Tbr1 activates Rybp and Auts2 to promote the formation of non-canonical Polycomb repressive complex 1 (PRC1). Also, Pax6, Tbr2, and Tbr1 collectively drive massive changes in the subunit isoform composition of BAF chromatin remodeling complexes during differentiation: for example, a novel switch from Bcl7c (Baf40c) to Bcl7a (Baf40a), the latter directly activated by Tbr2. Of 11 subunits predominantly in neuronal BAF, 7 were transcriptionally activated by Pax6, Tbr2, or Tbr1. Using EFs, Pax6→ Tbr2→ Tbr1 effect persistent changes of gene expression in cell lineages, to propagate features such as regional and laminar identity from progenitors to neurons.

Recurrent loss of heterozygosity correlates with clinical outcome in pancreatic neuroendocrine cancer

Pancreatic neuroendocrine tumors (pNETs) are uncommon cancers arising from pancreatic islet cells. Here we report the analysis of gene mutation, copy number, and RNA expression of 57 sporadic well-differentiated pNETs. pNET genomes are dominated by aneuploidy, leading to concordant changes in RNA expression at the level of whole chromosomes and chromosome segments. We observed two distinct patterns of somatic pNET aneuploidy that are associated with tumor pathology and patient prognosis. Approximately 26% of the patients in this series had pNETs with genomes characterized by recurrent loss of heterozygosity (LoH) of 10 specific chromosomes, accompanied by bi-allelic MEN1 inactivation and generally poor clinical outcome. Another ~40% of patients had pNETs that lacked this recurrent LoH pattern but had chromosome 11 LoH, bi-allelic MEN1 inactivation, and universally good clinical outcome. The somatic aneuploidy allowed pathogenic germline variants (e.g., ATM) to be expressed unopposed, with RNA expression patterns showing inactivation of downstream tumor suppressor pathways. No prognostic associations were found with tumor morphology, single gene mutation, or expression of RNAs reflecting the activity of immune, differentiation, proliferative or tumor suppressor pathways. In pNETs, single gene mutations appear to be less important than aneuploidy, with MEN1 the only statistically significant recurrently mutated driver gene. In addition, only one pNET in the series had clearly actionable single nucleotide variants (SNVs) (in PTEN and FLCN) confirmed by corroborating RNA expression changes. The two clinically relevant patterns of LoH described here define a novel oncogenic mechanism and a plausible route to genomic precision oncology for this tumor type.

Heterochromatin Protein 1α Mediates Development and Aggressiveness of Neuroendocrine Prostate Cancer

Neuroendocrine prostate cancer (NEPC) is a lethal subtype of prostate cancer arising mostly from adenocarcinoma via neuroendocrine transdifferentiation following androgen deprivation therapy. Mechanisms contributing to both NEPC development and its aggressiveness remain elusive. In light of the fact that hyperchromatic nuclei are a distinguishing histopathologic feature of NEPC, we utilized transcriptomic analyses of our patient-derived xenograft (PDX) models, multiple clinical cohorts, and genetically engineered mouse models to identify 36 heterochromatin-related genes that are significantly enriched in NEPC. Longitudinal analysis using our unique, first-in-field PDX model of adenocarcinoma-to-NEPC transdifferentiation revealed that, among those 36 heterochromatin-related genes, heterochromatin protein 1α (HP1α) expression increased early and steadily during NEPC development and remained elevated in the developed NEPC tumor. Its elevated expression was further confirmed in multiple PDX and clinical NEPC samples. HP1α knockdown in the NCI-H660 NEPC cell line inhibited proliferation, ablated colony formation, and induced apoptotic cell death, ultimately leading to tumor growth arrest. Its ectopic expression significantly promoted NE transdifferentiation in adenocarcinoma cells subjected to androgen deprivation treatment. Mechanistically, HP1α reduced expression of androgen receptor and RE1 silencing transcription factor and enriched the repressive trimethylated histone H3 at Lys9 mark on their respective gene promoters. These observations indicate a novel mechanism underlying NEPC development mediated by abnormally expressed heterochromatin genes, with HP1α as an early functional mediator and a potential therapeutic target for NEPC prevention and management.Significance: Heterochromatin proteins play a fundamental role in NEPC, illuminating new therapeutic targets for this aggressive disease. Cancer Res; 78(10); 2691-704. ©2018 AACR.

Expression profiling of cell-intrinsic regulators in the process of differentiation of human iPSCs into retinal lineages

BACKGROUND

Differentiation of human induced pluripotent stem cells (hiPSCs) into retinal lineages offers great potential for medical application. Therefore, it is of crucial importance to know the key intrinsic regulators of differentiation and the specific biomarker signatures of cell lineages.

METHODS

In this study, we used microarrays to analyze transcriptomes of terminally differentiated retinal ganglion cell (RGC) and retinal pigment epithelium (RPE) lineages, as well as intermediate retinal progenitor cells of optic vesicles (OVs) derived from hiPSCs. In our analysis, we specifically focused on the classes of transcripts that encode intrinsic regulators of gene expression: the transcription factors (TFs) and epigenetic chromatin state regulators. We applied two criteria for the selection of potentially important regulators and markers: firstly, the magnitude of fold-change of upregulation; secondly, the contrasted pattern of differential expression between OV, RGC and RPE lineages.

RESULTS

We found that among the most highly overexpressed TF-encoding genes in the OV/RGC lineage were three members of the Collier/Olfactory-1/Early B-cell family: EBF1, EBF2 and EBF3. Knockdown of EBF1 led to significant impairment of differentiation of hiPSCs into RGCs. EBF1 was shown to act upstream of ISL1 and BRN3A, the well-characterized regulators of RGC lineage specification. TF-encoding genes DLX1, DLX2 and INSM1 were the most highly overexpressed genes in the OVs, indicating their important role in the early stages of retinal differentiation. Along with MITF, the two paralogs, BHLHE41 and BHLHE40, were the most robust TF markers of RPE cells. The markedly contrasted expression of ACTL6B, encoding the component of chromatin remodeling complex SWI/SNF, discriminated hiPSC-derived OV/RGC and RPE lineages.

CONCLUSIONS

We identified novel, potentially important intrinsic regulators of RGC and RPE cell lineage specification in the process of differentiation from hiPSCs. We demonstrated the crucial role played by EBF1 in differentiation of RGCs. We identified intrinsic regulator biomarker signatures of these two retinal cell types that can be applied with high confidence to confirm the cell lineage identities.

Twin DNA Methylation Profiling Reveals Flare-Dependent Interferon Signature and B Cell Promoter Hypermethylation in Systemic Lupus Erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE) has limited monozygotic twin concordance, implying a role for pathogenic factors other than genetic variation, such as epigenetic changes. Using the disease-discordant twin model, we investigated genome-wide DNA methylation changes in sorted CD4+ T cells, monocytes, granulocytes, and B cells in twin pairs with at least 1 SLE-affected twin.

METHODS

Peripheral blood obtained from 15 SLE-affected twin pairs (6 monozygotic and 9 dizygotic) was processed using density-gradient centrifugation for the granulocyte fraction. CD4+ T cells, monocytes, and B cells were further isolated using magnetic beads. Genome-wide DNA methylation was analyzed using Infinium HumanMethylation450K BeadChips. When comparing probes from SLE-affected twins and co-twins, differential DNA methylation was considered statistically significant when the P value was less than 0.01 and biologically relevant when the median DNA methylation difference was >7%. Findings were validated by pyrosequencing and replicated in an independent case-control sample.

RESULTS

In paired analyses of twins discordant for SLE restricted to the gene promoter and start region, we identified 55, 327, 247, and 1,628 genes with differentially methylated CpGs in CD4+ T cells, monocytes, granulocytes, and B cells, respectively. All cell types displayed marked hypomethylation in interferon-regulated genes, such as IFI44L, PARP9, and IFITM1, which was more pronounced in twins who experienced a disease flare within the past 2 years. In contrast to what was observed in the other cell types, differentially methylated CpGs in B cells were predominantly hypermethylated, and the most important upstream regulators included TNF and EP300.

CONCLUSION

Hypomethylation of interferon-regulated genes occurs in all major cellular compartments in SLE-affected twins. The observed B cell promoter hypermethylation is a novel finding with potential significance in SLE pathogenesis.

Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes Meloidogyne spp., with a focus on the asexually reproducing species

Background

The renewed interest in epigenetics has led to the understanding that both the environment and individual lifestyle can directly interact with the epigenome to influence its dynamics. Epigenetic phenomena are mediated by DNA methylation, stable chromatin modifications and non-coding RNA-associated gene silencing involving specific proteins called epigenetic factors.

Multiple organisms, ranging from plants to yeast and mammals, have been used as model systems to study epigenetics. The interactions between parasites and their hosts are models of choice to study these mechanisms because the selective pressures are strong and the evolution is fast.

The asexually reproducing root-knot nematodes (RKN) offer different advantages to study the processes and mechanisms involved in epigenetic regulation. RKN genomes sequencing and annotation have identified numerous genes, however, which of those are involved in the adaption to an environment and potentially relevant to the evolution of plant-parasitism is yet to be discovered.

Results

Here, we used a functional comparative annotation strategy combining orthology data, mining of curated genomics as well as protein domain databases and phylogenetic reconstructions.

Overall, we show that (i) neither RKN, nor the model nematode Caenorhabditis elegans possess any DNA methyltransferases (DNMT) (ii) RKN do not possess the complete machinery for DNA methylation on the 6th position of adenine (6mA) (iii) histone (de)acetylation and (de)methylation pathways are conserved between C. elegans and RKN, and the corresponding genes are amplified in asexually reproducing RKN (iv) some specific non-coding RNA families found in plant-parasitic nematodes are dissimilar from those in C. elegans. In the asexually reproducing RKN Meloidogyne incognita, expression data from various developmental stages supported the putative role of these proteins in epigenetic regulations.

Conclusions

Our results refine previous predictions on the epigenetic machinery of model species and constitute the most comprehensive description of epigenetic factors relevant to the plant-parasitic lifestyle and/or asexual mode of reproduction of RKN. Providing an atlas of epigenetic factors in RKN is an informative resource that will enable researchers to explore their potential role in adaptation of these parasites to their environment.

Electronic supplementary material

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

Transcription profiling of peripheral B cells in antibody-positive primary Sjögren's syndrome reveals upregulated expression of CX3CR1 and a type I and type II interferon signature

B cells play a key role in the pathogenesis of primary Sjögren's syndrome (pSS). The aim of this study was to analyse the transcriptome of CD19+ B cells from patients with pSS and healthy controls to decipher the B cell-specific contribution to pSS. RNA from purified CD19+ B cells from 12 anti-SSA antibody-positive untreated female patients with pSS and 20 healthy blood donors was subjected to whole transcriptome sequencing. A false discovery rate corrected significance threshold of α < 0.05 was applied to define differential gene expression. As validation, gene expression in B cells from 17 patients with pSS and 16 healthy controls was analysed using a targeted gene panel. RNA-sequencing identified 4047 differentially expressed autosomal genes in pSS B cells. Upregulated expression of type I and type II interferon (IFN)-induced genes was observed, establishing an IFN signature in pSS B cells. Among the top upregulated and validated genes were CX3CR1, encoding the fractalkine receptor involved in regulation of B-cell malignancies, CCL5/RANTES and CCR1. Increased expression of several members of the TNF superfamily was also identified; TNFSF4/Ox40L, TNFSF10/TRAIL, TNFSF13B/BAFF, TNFRSF17/BCMA as well as S100A8 and -A9/calprotectin, TLR7, STAT1 and STAT2. Among genes with downregulated expression in pSS B cells were SOCS1 and SOCS3, CD70 and TNFAIP3/A20. We conclude that B cells from patients with anti-SSA antibody-positive pSS display immune activation with upregulated expression of chemokines, chemokine receptors and a prominent type I and type II IFN signature, while suppressors of cytokine signalling are downregulated. This adds insight into the autoimmune process and suggests potential targets for future functional studies.

Enhancer Remodeling and MicroRNA Alterations Are Associated with Acquired Resistance to ALK Inhibitors

Anaplastic lymphoma kinase (ALK) inhibitors are highly effective in patients with ALK fusion-positive lung cancer, but acquired resistance invariably emerges. Identification of secondary mutations has received considerable attention, but most cases cannot be explained by genetic causes alone, raising the possibility of epigenetic mechanisms in acquired drug resistance. Here, we investigated the dynamic changes in the transcriptome and enhancer landscape during development of acquired resistance to ALK inhibitors. Histone H3 lysine 27 acetylation (H3K27ac) was profoundly altered during acquisition of resistance, and enhancer remodeling induced expression changes in both miRNAs and mRNAs. Decreased H3K27ac levels and reduced miR-34a expression associated with the activation of target genes such as AXL. Panobinostat, a pan-histone deacetylase inhibitor, altered the H3K27ac profile and activated tumor-suppressor miRNAs such as miR-449, another member of the miR-34 family, and synergistically induced antiproliferative effects with ALK inhibitors on resistant cells, xenografts, and EML4-ALK transgenic mice. Paired analysis of patient samples before and after treatment with ALK inhibitors revealed that repression of miR-34a or miR-449a and activation of AXL were mutually exclusive of secondary mutations in ALK. Our findings indicate that enhancer remodeling and altered expression of miRNAs play key roles in cancer drug resistance and suggest that strategies targeting epigenetic pathways represent a potentially effective method for overcoming acquired resistance to cancer therapy.Significance: Epigenetic deregulation drives acquired resistance to ALK inhibitors in ALK-positive lung cancer. Cancer Res; 78(12); 3350-62. ©2018 AACR.

In silico cancer research towards 3R

BACKGROUND

Improving our understanding of cancer and other complex diseases requires integrating diverse data sets and algorithms. Intertwining in vivo and in vitro data and in silico models are paramount to overcome intrinsic difficulties given by data complexity. Importantly, this approach also helps to uncover underlying molecular mechanisms. Over the years, research has introduced multiple biochemical and computational methods to study the disease, many of which require animal experiments. However, modeling systems and the comparison of cellular processes in both eukaryotes and prokaryotes help to understand specific aspects of uncontrolled cell growth, eventually leading to improved planning of future experiments. According to the principles for humane techniques milestones in alternative animal testing involve in vitro methods such as cell-based models and microfluidic chips, as well as clinical tests of microdosing and imaging. Up-to-date, the range of alternative methods has expanded towards computational approaches, based on the use of information from past in vitro and in vivo experiments. In fact, in silico techniques are often underrated but can be vital to understanding fundamental processes in cancer. They can rival accuracy of biological assays, and they can provide essential focus and direction to reduce experimental cost.

MAIN BODY

We give an overview on in vivo, in vitro and in silico methods used in cancer research. Common models as cell-lines, xenografts, or genetically modified rodents reflect relevant pathological processes to a different degree, but can not replicate the full spectrum of human disease. There is an increasing importance of computational biology, advancing from the task of assisting biological analysis with network biology approaches as the basis for understanding a cell's functional organization up to model building for predictive systems.

CONCLUSION

Underlining and extending the in silico approach with respect to the 3Rs for replacement, reduction and refinement will lead cancer research towards efficient and effective precision medicine. Therefore, we suggest refined translational models and testing methods based on integrative analyses and the incorporation of computational biology within cancer research.

Neuron-specific alternative splicing of transcriptional machineries: Implications for neurodevelopmental disorders

The brain has long been known to display the most complex pattern of alternative splicing, thereby producing diverse protein isoforms compared to other tissues. Recent evidence indicates that many alternative exons are neuron-specific, evolutionarily conserved, and found in regulators of transcription including DNA-binding protein and histone modifying enzymes. This raises a possibility that neurons adopt unique mechanisms of transcription. Given that transcriptional machineries are frequently mutated in neurodevelopmental disorders with cognitive dysfunction, it is important to understand how neuron-specific alternative splicing contributes to proper transcriptional regulation in the brain. In this review, we summarize current knowledge regarding how neuron-specific splicing events alter the function of transcriptional regulators and shape unique gene expression patterns in the brain and the implications of neuronal splicing to the pathophysiology of neurodevelopmental disorders.

B-CAN: a resource sharing platform to improve the operation, visualization and integrated analysis of TCGA breast cancer data

Breast cancer is a high-risk heterogeneous disease with myriad subtypes and complicated biological features. The Cancer Genome Atlas (TCGA) breast cancer database provides researchers with the large-scale genome and clinical data via web portals and FTP services. Researchers are able to gain new insights into their related fields, and evaluate experimental discoveries with TCGA. However, it is difficult for researchers who have little experience with database and bioinformatics to access and operate on because of TCGA's complex data format and diverse files. For ease of use, we build the breast cancer (B-CAN) platform, which enables data customization, data visualization, and private data center. The B-CAN platform runs on Apache server and interacts with the backstage of MySQL database by PHP. Users can customize data based on their needs by combining tables from original TCGA database and selecting variables from each table. The private data center is applicable for private data and two types of customized data. A key feature of the B-CAN is that it provides single table display and multiple table display. Customized data with one barcode corresponding to many records and processed customized data are allowed in Multiple Tables Display. The B-CAN is an intuitive and high-efficient data-sharing platform.

Monitoring transcription initiation activities in rat and dog

The promoter landscape of several non-human model organisms is far from complete. As a part of FANTOM5 data collection, we generated 13 profiles of transcription initiation activities in dog and rat aortic smooth muscle cells, mesenchymal stem cells and hepatocytes by employing CAGE (Cap Analysis of Gene Expression) technology combined with single molecule sequencing. Our analyses show that the CAGE profiles recapitulate known transcription start sites (TSSs) consistently, in addition to uncover novel TSSs. Our dataset can be thus used with high confidence to support gene annotation in dog and rat species. We identified 28,497 and 23,147 CAGE peaks, or promoter regions, for rat and dog respectively, and associated them to known genes. This approach could be seen as a standard method for improvement of existing gene models, as well as discovery of novel genes. Given that the FANTOM5 data collection includes dog and rat matched cell types in human and mouse as well, this data would also be useful for cross-species studies.

P53 at the start of the 21st century: lessons from elephants

Crucial, natural protection against tumour onset in humans is orchestrated by the dynamic protein p53. The best-characterised functions of p53 relate to its cellular stress responses. In this review, we explore emerging insights into p53 activities and their functional consequences. We compare p53 in humans and elephants, in search of salient features of cancer protection.

Clinical significance of miRNA host gene promoter methylation in prostate cancer

Only a part of prostate cancer (PCa) patients has aggressive malignancy requiring adjuvant treatment after radical prostatectomy (RP). Biomarkers capable to predict biochemical PCa recurrence (BCR) after RP would significantly improve preoperative risk stratification and treatment decisions. MicroRNA (miRNA) deregulation has recently emerged as an important phenomenon in tumor development and progression, however, the mechanisms remain largely unstudied. In the present study, based on microarray profiling of DNA methylation in 9 pairs of PCa and noncancerous prostate tissues (NPT), host genes of miR-155-5p, miR-152-3p, miR-137, miR-31-5p, and miR-642a, -b were analyzed for promoter methylation in 129 PCa, 35 NPT, and 17 benign prostatic hyperplasia samples (BPH) and compared to the expression of mature miRNAs and their selected targets (DNMT1, KDM1A, and KDM5B). The Cancer Genome Atlas dataset was utilized for validation. Methylation of mir-155, mir-152, and mir-137 host genes was PCa-specific, and downregulation of miR-155-5p significantly correlated with promoter methylation. Higher KDM5B expression was observed in samples with methylated mir-155 or mir-137 promoters, whereas upregulation of KDM1A and DNMT1 was associated with mir-155 and mir-152 methylation status, respectively. Promoter methylation of mir-155, mir-152, and mir-31 was predictive of BCR-free survival in various Cox models and increased the prognostic value of clinicopathologic factors. In conclusion, methylated mir-155, mir-152, mir-137, and mir-31 host genes are promising diagnostic and/or prognostic biomarkers of PCa. Methylation status of particular miRNA host genes as independent variables or in combinations might assist physicians in identifying poor prognosis PCa patients preoperatively.

Insights into beta cell regeneration for diabetes via integration of molecular landscapes in human insulinomas

Although diabetes results in part from a deficiency of normal pancreatic beta cells, inducing human beta cells to regenerate is difficult. Reasoning that insulinomas hold the “genomic recipe” for beta cell expansion, we surveyed 38 human insulinomas to obtain insights into therapeutic pathways for beta cell regeneration. An integrative analysis of whole-exome and RNA-sequencing data was employed to extensively characterize the genomic and molecular landscape of insulinomas relative to normal beta cells. Here, we show at the pathway level that the majority of the insulinomas display mutations, copy number variants and/or dysregulation of epigenetic modifying genes, most prominently in the polycomb and trithorax families. Importantly, these processes are coupled to co-expression network modules associated with cell proliferation, revealing candidates for inducing beta cell regeneration. Validation of key computational predictions supports the concept that understanding the molecular complexity of insulinoma may be a valuable approach to diabetes drug discovery.

Diabetes results in part from a deficiency of functional pancreatic beta cells. Here, the authors study the genomic and epigenetic landscapes of human insulinomas to gain insight into possible pathways for therapeutic beta cell regeneration, highlighting epigenetic genes and pathways.

Epigenome alterations in aortic valve stenosis and its related left ventricular hypertrophy

Aortic valve stenosis is the most common cardiac valve disease, and with current trends in the population demographics, its prevalence is likely to rise, thus posing a major health and economic burden facing the worldwide societies. Over the past decade, it has become more than clear that our traditional genetic views do not sufficiently explain the well-known link between AS, proatherogenic risk factors, flow-induced mechanical forces, and disease-prone environmental influences. Recent breakthroughs in the field of epigenetics offer us a new perspective on gene regulation, which has broadened our perspective on etiology of aortic stenosis and other aortic valve diseases. Since all known epigenetic marks are potentially reversible this perspective is especially exciting given the potential for development of successful and non-invasive therapeutic intervention and reprogramming of cells at the epigenetic level even in the early stages of disease progression. This review will examine the known relationships between four major epigenetic mechanisms: DNA methylation, posttranslational histone modification, ATP-dependent chromatin remodeling, and non-coding regulatory RNAs, and initiation and progression of AS. Numerous profiling and functional studies indicate that they could contribute to endothelial dysfunctions, disease-prone activation of monocyte-macrophage and circulatory osteoprogenitor cells and activation and osteogenic transdifferentiation of aortic valve interstitial cells, thus leading to valvular inflammation, fibrosis, and calcification, and to pressure overload-induced maladaptive myocardial remodeling and left ventricular hypertrophy. This is especcialy the case for small non-coding microRNAs but was also, although in a smaller scale, convincingly demonstrated for other members of cellular epigenome landscape. Equally important, and clinically most relevant, the reported data indicate that epigenetic marks, particularly certain microRNA signatures, could represent useful non-invasive biomarkers that reflect the disease progression and patients prognosis for recovery after the valve replacement surgery.

SWI/SNF Infobase-An exclusive information portal for SWI/SNF remodeling complex subunits

Chromatin remodeling complexes facilitate the access of condensed genomic DNA during transcription, replication, and repair, by altering the histone-DNA contacts in the nucleosome structures. SWI/SNF (SWItch/Sucrose Non-Fermentable) family of ATP dependent chromatin remodeling complexes have been documented for their tumour suppressor function. Recent studies have reported the high frequency of cancer causing mutations in this protein family. There exist multiple subunits for this complex and can form context-dependent sub-complexes. The cataloguing of individual subunits of this complex is essential for understanding their specific functions and their mechanism of action during chromatin remodeling. This would also facilitate further studies to characterize cancer causing mutations in SWI/SNF subunits. In the current study, a database containing information on the subunits of SWI/SNF-α (BRG1/BRM-Associated Factors (BAF)) and SWI/SNF-β (Polybromo-Associated BAF (PBAF)) sub classes of SWI/SNF family has been curated and catalogued. The database hosts information on 27 distinct SWI/SNF subunits from 20 organisms spanning a wide evolutionary range of eukaryotes. A non-redundant set of 522 genes coding for SWI/SNF subunits have been documented in the database. A detailed annotation on each subunit, including basic protein/gene information, protein sequence, functional domains, homologs and missense mutations of human proteins have been provided with a user-friendly graphical interface. The SWI/SNF Infobase presented here, would be a first of its kind exclusive information portal on SWI/SNF complex subunits and would be a valuable resource for the research community working on chromatin remodeling. The database is available at http://scbt.sastra.edu/swisnfdb/index.php.

A computational method using the random walk with restart algorithm for identifying novel epigenetic factors

Epigenetic regulation has long been recognized as a significant factor in various biological processes, such as development, transcriptional regulation, spermatogenesis, and chromosome stabilization. Epigenetic alterations lead to many human diseases, including cancer, depression, autism, and immune system defects. Although efforts have been made to identify epigenetic regulators, it remains a challenge to systematically uncover all the components of the epigenetic regulation in the genome level using experimental approaches. The advances of constructing protein-protein interaction (PPI) networks provide an excellent opportunity to identify novel epigenetic factors computationally in the genome level. In this study, we identified potential epigenetic factors by using a computational method that applied the random walk with restart (RWR) algorithm on a protein-protein interaction (PPI) network using reported epigenetic factors as seed nodes. False positives were identified by their specific roles in the PPI network or by a low-confidence interaction and a weak functional relationship with epigenetic regulators. After filtering out the false positives, 26 candidate epigenetic factors were finally accessed. According to previous studies, 22 of these are thought to be involved in epigenetic regulation, suggesting the robustness of our method. Our study provides a novel computational approach which successfully identified 26 potential epigenetic factors, paving the way on deepening our understandings on the epigenetic mechanism.

Integrative CAGE and DNA Methylation Profiling Identify Epigenetically Regulated Genes in NSCLC

Lung cancer is the leading cause of cancer-related deaths worldwide. The majority of cancer driver mutations have been identified; however, relevant epigenetic regulation involved in tumorigenesis has only been fragmentarily analyzed. Epigenetically regulated genes have a great theranostic potential, especially in tumors with no apparent driver mutations. Here, epigenetically regulated genes were identified in lung cancer by an integrative analysis of promoter-level expression profiles from Cap Analysis of Gene Expression (CAGE) of 16 non-small cell lung cancer (NSCLC) cell lines and 16 normal lung primary cell specimens with DNA methylation data of 69 NSCLC cell lines and 6 normal lung epithelial cells. A core set of 49 coding genes and 10 long noncoding RNAs (lncRNA), which are upregulated in NSCLC cell lines due to promoter hypomethylation, was uncovered. Twenty-two epigenetically regulated genes were validated (upregulated genes with hypomethylated promoters) in the adenocarcinoma and squamous cell cancer subtypes of lung cancer using The Cancer Genome Atlas data. Furthermore, it was demonstrated that multiple copies of the REP522 DNA repeat family are prominently upregulated due to hypomethylation in NSCLC cell lines, which leads to cancer-specific expression of lncRNAs, such as RP1-90G24.10, AL022344.4, and PCAT7. Finally, Myeloma Overexpressed (MYEOV) was identified as the most promising candidate. Functional studies demonstrated that MYEOV promotes cell proliferation, survival, and invasion. Moreover, high MYEOV expression levels were associated with poor prognosis.Implications: This report identifies a robust list of 22 candidate driver genes that are epigenetically regulated in lung cancer; such genes may complement the known mutational drivers.Visual Overview: http://mcr.aacrjournals.org/content/molcanres/15/10/1354/F1.large.jpg Mol Cancer Res; 15(10); 1354-65. ©2017 AACR.

Transforming Cancer Epigenetics Using Nutritive Approaches and Noncoding RNAs

Cancer is considered one of the leading causes of death in the United States. Although preventive strategies, early detection, and improved treatment options have been developed, novel targets and therapeutics are still needed. Since concluding that cancer is mediated by genetic and epigenetic alterations of the cell, many research groups are now focusing on other means of prevention and therapy via nutrition, epigenetic mechanisms, and non-coding RNAs which have been shown to control gene expression and have many different functions at the cellular level. With the advent of high-throughput sequencing in human cancer, the potential to identify novel biomarkers and therapeutic targets of disease has increased tremendously and led to the identification of many non-coding RNAs that are dysregulated in various cancers. Gene expression and regulation is important in maintaining the homeostasis of normal tissues and cells. Not uncommonly, up- or down-regulation of particular genes are associated with cancer as a result of increased or decreased expression of transcriptional targets. This review focuses on the role of nutrition in cancer and the dysregulation of non-coding RNAs with particular emphasis on long non-coding RNAs and microRNAs in different cancer types.

Interactions Between Genetics, Lifestyle, and Environmental Factors for Healthcare

The occurrence and progression of diseases are strongly associated with a combination of genetic, lifestyle, and environmental factors. Understanding the interplay between genetic and nongenetic components provides deep insights into disease pathogenesis and promotes personalized strategies for people healthcare. Recently, the paradigm of systems medicine, which integrates biomedical data and knowledge at multidimensional levels, is considered to be an optimal way for disease management and clinical decision-making in the era of precision medicine. In this chapter, epigenetic-mediated genetics-lifestyle-environment interactions within specific diseases and different ethnic groups are systematically discussed, and data sources, computational models, and translational platforms for systems medicine research are sequentially presented. Moreover, feasible suggestions on precision healthcare and healthy longevity are kindly proposed based on the comprehensive review of current studies.

ARN: Analysis and Visualization System for Adipogenic Regulation Network Information

Adipogenesis is the process of cell differentiation through which preadipocytes become adipocytes. Lots of research is currently ongoing to identify genes, including their gene products and microRNAs, that correlate with fat cell development. However, information fragmentation hampers the identification of key regulatory genes and pathways. Here, we present a database of literature-curated adipogenesis-related regulatory interactions, designated the Adipogenesis Regulation Network (ARN, http://210.27.80.93/arn/), which currently contains 3101 nodes (genes and microRNAs), 1863 regulatory interactions, and 33,969 expression records associated with adipogenesis, based on 1619 papers. A sentence-based text-mining approach was employed for efficient manual curation of regulatory interactions from approximately 37,000 PubMed abstracts. Additionally, we further determined 13,103 possible node relationships by searching miRGate, BioGRID, PAZAR and TRRUST. ARN also has several useful features: i) regulatory map information; ii) tests to examine the impact of a query node on adipogenesis; iii) tests for the interactions and modes of a query node; iv) prediction of interactions of a query node; and v) analysis of experimental data or the construction of hypotheses related to adipogenesis. In summary, ARN can store, retrieve and analyze adipogenesis-related information as well as support ongoing adipogenesis research and contribute to the discovery of key regulatory genes and pathways.

Gene expression profiling of epigenetic chromatin modification enzymes and histone marks by cigarette smoke: implications for COPD and lung cancer

Chromatin-modifying enzymes mediate DNA methylation and histone modifications on recruitment to specific target gene loci in response to various stimuli. The key enzymes that regulate chromatin accessibility for maintenance of modifications in DNA and histones, and for modulation of gene expression patterns in response to cigarette smoke (CS), are not known. We hypothesize that CS exposure alters the gene expression patterns of chromatin-modifying enzymes, which then affects multiple downstream pathways involved in the response to CS. We have, therefore, analyzed chromatin-modifying enzyme profiles and validated by quantitative real-time PCR (qPCR). We also performed immunoblot analysis of targeted histone marks in C57BL/6J mice exposed to acute and subchronic CS, and of lungs from nonsmokers, smokers, and patients with chronic obstructive pulmonary disease (COPD). We found a significant increase in expression of several chromatin modification enzymes, including DNA methyltransferases, histone acetyltransferases, histone methyltransferases, and SET domain proteins, histone kinases, and ubiquitinases. Our qPCR validation data revealed a significant downregulation of Dnmt1, Dnmt3a, Dnmt3b, Hdac2, Hdac4, Hat1, Prmt1, and Aurkb We identified targeted chromatin histone marks (H3K56ac and H4K12ac), which are induced by CS. Thus CS-induced genotoxic stress differentially affects the expression of epigenetic modulators that regulate transcription of target genes via DNA methylation and site-specific histone modifications. This may have implications in devising epigenetic-based therapies for COPD and lung cancer.

Update of the FANTOM web resource: high resolution transcriptome of diverse cell types in mammals

Upon the first publication of the fifth iteration of the Functional Annotation of Mammalian Genomes collaborative project, FANTOM5, we gathered a series of primary data and database systems into the FANTOM web resource (http://fantom.gsc.riken.jp) to facilitate researchers to explore transcriptional regulation and cellular states. In the course of the collaboration, primary data and analysis results have been expanded, and functionalities of the database systems enhanced. We believe that our data and web systems are invaluable resources, and we think the scientific community will benefit for this recent update to deepen their understanding of mammalian cellular organization. We introduce the contents of FANTOM5 here, report recent updates in the web resource and provide future perspectives.

Preservation of methylated CpG dinucleotides in human CpG islands

BACKGROUND

CpG dinucleotides are extensively underrepresented in mammalian genomes. It is widely accepted that genome-wide CpG depletion is predominantly caused by an elevated CpG > TpG mutation rate due to frequent cytosine methylation in the CpG context. Meanwhile the CpG content in genomic regions called CpG islands (CGIs) is noticeably higher. This observation is usually explained by lower CpG > TpG substitution rates within CGIs due to reduced cytosine methylation levels.

RESULTS

By combining genome-wide data on substitutions and methylation levels in several human cell types we have shown that cytosine methylation in human sperm cells was strongly and consistently associated with increased CpG > TpG substitution rates. In contrast, this correlation was not observed for embryonic stem cells or fibroblasts. Surprisingly, the decreased sperm CpG methylation level was insufficient to explain the reduced CpG > TpG substitution rates in CGIs.

CONCLUSIONS

While cytosine methylation in human sperm cells is strongly associated with increased CpG > TpG substitution rates, substitution rates are significantly reduced within CGIs even after sperm CpG methylation levels and local GC content are controlled for. Our findings are consistent with strong negative selection preserving methylated CpGs within CGIs including intergenic ones.