The chromatin accessibility landscape of primary human cancers

The Promise for Histone Methyltransferase Inhibitors for Epigenetic Therapy in Clinical Oncology: A Narrative Review

Epigenetic processes are essential for normal development and the maintenance of tissue-specific gene expression in mammals. Changes in gene expression and malignant cellular transformation can result from disruption of epigenetic mechanisms, and global disruption in the epigenetic landscape is a key feature of cancer. The study of epigenetics in cancer has revealed that human cancer cells harbor both genetic alterations and epigenetic abnormalities that interplay at all stages of cancer development. Unlike genetic mutations, epigenetic aberrations are potentially reversible through epigenetic therapy, providing a therapeutically relevant treatment option. Histone methyltransferase inhibitors are emerging as an epigenetic therapy approach with great promise in the field of clinical oncology. The recent accelerated approval of the enhancer of zeste homolog 2 (EZH2; also known as histone-lysine N-methyltransferase EZH2) inhibitor tazemetostat for metastatic or locally advanced epithelioid sarcoma marks the first approval of such a compound for the treatment of cancer. Many other histone methyltransferase inhibitors are currently in development, some of which are being tested in clinical studies. This review focuses on histone methyltransferase inhibitors, highlighting their potential in the treatment of cancer. We also discuss the role for such epigenetic drugs in overcoming epigenetically driven drug resistance mechanisms, and their value in combination with other therapeutic approaches such as immunotherapy.

Navigating the path to distant metastasis

How ‘difficult’ is it for somatic evolution to produce a cell that is capable of leaving the primary tumor and growing in a distant organ? In this issue, Reiter et al. assess genetic diversity across metastatic lesions and identify a tight selective bottleneck preceding distant metastasis.

Allele-specific DNA methylation is increased in cancers and its dense mapping in normal plus neoplastic cells increases the yield of disease-associated regulatory SNPs

BACKGROUND

Mapping of allele-specific DNA methylation (ASM) can be a post-GWAS strategy for localizing regulatory sequence polymorphisms (rSNPs). The advantages of this approach, and the mechanisms underlying ASM in normal and neoplastic cells, remain to be clarified.

RESULTS

We perform whole genome methyl-seq on diverse normal cells and tissues and three cancer types. After excluding imprinting, the data pinpoint 15,112 high-confidence ASM differentially methylated regions, of which 1838 contain SNPs in strong linkage disequilibrium or coinciding with GWAS peaks. ASM frequencies are increased in cancers versus matched normal tissues, due to widespread allele-specific hypomethylation and focal allele-specific hypermethylation in poised chromatin. Cancer cells show increased allele switching at ASM loci, but disruptive SNPs in specific classes of CTCF and transcription factor binding motifs are similarly correlated with ASM in cancer and non-cancer. Rare somatic mutations affecting these same motif classes track with de novo ASM. Allele-specific transcription factor binding from ChIP-seq is enriched among ASM loci, but most ASM differentially methylated regions lack such annotations, and some are found in otherwise uninformative "chromatin deserts."

CONCLUSIONS

ASM is increased in cancers but occurs by a shared mechanism involving disruptive SNPs in CTCF and transcription factor binding sites in both normal and neoplastic cells. Dense ASM mapping in normal plus cancer samples reveals candidate rSNPs that are difficult to find by other approaches. Together with GWAS data, these rSNPs can nominate specific transcriptional pathways in susceptibility to autoimmune, cardiometabolic, neuropsychiatric, and neoplastic diseases.

3Scover: Identifying Safeguard TF from Cell Type-TF Specificity Network by an Extended Minimum Set Cover Model

Transcription factors (TFs) define cellular identity either by activating target cell program or by silencing donor program as demonstrated by intensive cell reprogramming studies. Here, we propose an extended minimum set cover model with stable selection (3Scover) to systematically identify silencing TFs, named safeguard TFs, from omics data. First, a cell type-TF specificity network is constructed to systematically link cell types with their specifically expressed TFs. Then we search the minimum TF set to cover this network with “many but one specificity” characteristic and integrate many subsampling models for a stable solution. 3Scover identified 30 safeguard TFs in human and mouse. These safeguard TFs are significantly enriched in the experimentally discovered reprogramming panel with their protein-protein interactors. In addition, they tend to interact closely with chromatin regulators, negatively regulate transcription, and function earlier in development. Collectively, 3Scover allows us to probe master TFs and combinatorial regulation in controlling cell identity.

•

Cell type-TF specificity networks reveal the relationships among TF and cell identity

•

3SCover extracts safeguard TFs by “many but one specificity” and parsimony principle

•

Safeguard TFs are enriched in reprogramming panel and interact closely with CR

•

Safeguard TFs are conserved in mouse and human

Super-enhancer acquisition drives oncogene expression in triple negative breast cancer

Triple Negative Breast Cancer (TNBC) is a heterogeneous disease lacking known molecular drivers and effective targeted therapies. Cytotoxic chemotherapy remains the mainstay of treatment for TNBCs, which have significantly poorer survival rates compared to other breast cancer subtypes. In addition to changes within the coding genome, aberrant enhancer activity is a well-established contributor to tumorigenesis. Here we use H3K27Ac chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map the active cis-regulatory landscape in TNBC. We identify distinct disease subtypes associated with specific enhancer activity, and over 2,500 unique superenhancers acquired by tumor cells but absent from normal breast tissue. To identify potential actionable disease drivers, we probed the dependency on genes that associate with tumor-specific enhancers by CRISPR screening. In this way we identify a number of tumor-specific dependencies, including a previously uncharacterized dependency on the TGFβ pseudo-receptor BAMBI.

Joint reconstruction of cis-regulatory interaction networks across multiple tissues using single-cell chromatin accessibility data

The rapid accumulation of single-cell chromatin accessibility data offers a unique opportunity to investigate common and specific regulatory mechanisms across different cell types. However, existing methods for cis-regulatory network reconstruction using single-cell chromatin accessibility data were only designed for cells belonging to one cell type, and resulting networks may be incomparable directly due to diverse cell numbers of different cell types. Here, we adopt a computational method to jointly reconstruct cis-regulatory interaction maps (JRIM) of multiple cell populations based on patterns of co-accessibility in single-cell data. We applied JRIM to explore common and specific regulatory interactions across multiple tissues from single-cell ATAC-seq dataset containing ~80 000 cells across 13 mouse tissues. Reconstructed common interactions among 13 tissues indeed relate to basic biological functions, and individual cis-regulatory networks show strong tissue specificity and functional relevance. More importantly, tissue-specific regulatory interactions are mediated by coordination of histone modifications and tissue-related TFs, and many of them may reveal novel regulatory mechanisms.

Analytical Approaches for ATAC-seq Data Analysis

ATAC-seq, the assay for transposase-accessible chromatin using sequencing, is a quick and efficient approach to investigating the chromatin accessibility landscape. Investigating chromatin accessibility has broad utility for answering many biological questions, such as mapping nucleosomes, identifying transcription factor binding sites, and measuring differential activity of DNA regulatory elements. Because the ATAC-seq protocol is both simple and relatively inexpensive, there has been a rapid increase in the availability of chromatin accessibility data. Furthermore, advances in ATAC-seq protocols are rapidly extending its breadth to additional experimental conditions, cell types, and species. Accompanying the increase in data, there has also been an explosion of new tools and analytical approaches for analyzing it. Here, we explain the fundamentals of ATAC-seq data processing, summarize common analysis approaches, and review computational tools to provide recommendations for different research questions. This primer provides a starting point and a reference for analysis of ATAC-seq data. © 2020 Wiley Periodicals LLC.

MEDEA: analysis of transcription factor binding motifs in accessible chromatin

Deciphering the interplay between chromatin accessibility and transcription factor (TF) binding is fundamental to understanding transcriptional regulation, control of cellular states, and the establishment of new phenotypes. Recent genome-wide chromatin accessibility profiling studies have provided catalogs of putative open regions, where TFs can recognize their motifs and regulate gene expression programs. Here, we present motif enrichment in differential elements of accessibility (MEDEA), a computational tool that analyzes high-throughput chromatin accessibility genomic data to identify cell-type-specific accessible regions and lineage-specific motifs associated with TF binding therein. To benchmark MEDEA, we used a panel of reference cell lines profiled by ENCODE and curated by the ENCODE Project Consortium for the ENCODE-DREAM Challenge. By comparing results with RNA-seq data, ChIP-seq peaks, and DNase-seq footprints, we show that MEDEA improves the detection of motifs associated with known lineage specifiers. We then applied MEDEA to 610 ENCODE DNase-seq data sets, where it revealed significant motifs even when absolute enrichment was low and where it identified novel regulators, such as NRF1 in kidney development. Finally, we show that MEDEA performs well on both bulk and single-cell ATAC-seq data. MEDEA is publicly available as part of our Glossary-GENRE suite for motif enrichment analysis.

B-cell lymphoma 2 family genes show a molecular pattern of spatiotemporal heterogeneity in gynaecologic and breast cancer

Objectives

BCL2 family proteins have been widely studied over the past decade due to their essential roles in apoptosis, oncogenesis and anti‐cancer therapy. However, the similarities and differences in the spatial pattern of the BCL2 gene family within the context of chromatin have not been well characterized. We sought to fill this knowledge gap by assessing correlations between gene alteration, gene expression, chromatin accessibility, and clinical outcomes in gynaecologic and breast cancer.

Materials and methods

In this study, the molecular characteristics of the BCL2 gene family in gynaecologic cancer were systematically analysed by integrating multi‐omics datasets, including transcriptomics, chromatin accessibility, copy number variation, methylomics and clinical outcome.

Results

We evaluated spatiotemporal associations between long‐range regulation peaks and tumour heterogeneity. Differential expression of the BCL2 family was coupled with widespread chromatin accessibility changes in gynaecologic cancer, accompanied by highly heterogeneous distal non‐coding accessibility surrounding the BCL2L1 gene loci. A relationship was also identified between gene expression, gene amplification, enhancer signatures, DNA methylation and overall patient survival. Prognostic analysis implied clinical correlations with BAD, BIK and BAK1. A shared protein regulatory network was established in which the co‐mutation signature of TP53 and PIK3CA was linked to the BCL2L1 gene.

Conclusions

Our results provide the first systematic identification of the molecular features of the BCL2 family under the spatial pattern of chromatin in gynaecologic and breast cancer. These findings broaden the therapeutic scope of the BCL2 family to the non‐coding region by including a significantly conserved distal region overlaying an enhancer.

Unique and assay specific features of NOMe-, ATAC- and DNase I-seq data

Chromatin accessibility maps are important for the functional interpretation of the genome. Here, we systematically analysed assay specific differences between DNase I-seq, ATAC-seq and NOMe-seq in a side by side experimental and bioinformatic setup. We observe that most prominent nucleosome depleted regions (NDRs, e.g. in promoters) are roboustly called by all three or at least two assays. However, we also find a high proportion of assay specific NDRs that are often ‘called’ by only one of the assays. We show evidence that these assay specific NDRs are indeed genuine open chromatin sites and contribute important information for accurate gene expression prediction. While technically ATAC-seq and DNase I-seq provide a superb high NDR calling rate for relatively low sequencing costs in comparison to NOMe-seq, NOMe-seq singles out for its genome-wide coverage allowing to not only detect NDRs but also endogenous DNA methylation and as we show here genome wide segmentation into heterochromatic B domains and local phasing of nucleosomes outside of NDRs. In summary, our comparisons strongly suggest to consider assay specific differences for the experimental design and for generalized and comparative functional interpretations.

Clinical implications of chromatin accessibility in human cancers

Assay for transposase-accessible chromatin using sequencing (ATAC-seq) has not yet been widely used in cancer research. Clinical implications of chromatin accessibility assessed by ATAC-seq profiling in human cancers especially in a large patient cohort is largely unknown. In this study, we analyzed ATAC-seq data in 404 cancer patients from the Cancer Genome Atlas, representing the largest cancer patient cohort with ATAC-seq data, and correlated chromatin accessibility with patient demographics, tumor histology, molecular subtypes, and survival. Our results showed that chromatin accessibility varies from chromosome to chromosome, and is different in different genomic regions along the same chromosome. Chromatin accessibility especially on the X chromosome is strongly dependent on patient sex, but not much on patient age or tumor stage. Striking difference in chromatin accessibility is observed between lung adenocarcinoma and lung squamous cell carcinoma, the two most common histological subgroups in lung cancer. Furthermore, chromatin accessibility was different between basal and non-basal breast cancer. Finally, we identified prognostic peaks in the promoter regions that were significantly correlated with survival. In particular, we identified six peaks in the ESR1 gene promoter region in the ATAC-seq profiling and found that the peak about 247 bp away from the transcription start site was significantly associated with better survival. In conclusion, our study provides an alternative mechanism underlying tumor prognosis.

Epigenetic driver mutations in ARID1A shape cancer immune phenotype and immunotherapy

Whether mutations in cancer driver genes directly affect cancer immune phenotype and T cell immunity remains a standing question. ARID1A is a core member of the polymorphic BRG/BRM-associated factor chromatin remodeling complex. ARID1A mutations occur in human cancers and drive cancer development. Here, we studied the molecular, cellular, and clinical impact of ARID1A aberrations on cancer immunity. We demonstrated that ARID1A aberrations resulted in limited chromatin accessibility to IFN-responsive genes, impaired IFN gene expression, anemic T cell tumor infiltration, poor tumor immunity, and shortened host survival in many human cancer histologies and in murine cancer models. Impaired IFN signaling was associated with poor immunotherapy response. Mechanistically, ARID1A interacted with EZH2 via its carboxyl terminal and antagonized EZH2-mediated IFN responsiveness. Thus, the interaction between ARID1A and EZH2 defines cancer IFN responsiveness and immune evasion. Our work indicates that cancer epigenetic driver mutations can shape cancer immune phenotype and immunotherapy.

ATAC-seq identifies thousands of extrachromosomal circular DNA in cancer and cell lines

Extrachromosomal circular DNAs (eccDNAs) are somatically mosaic and contribute to intercellular heterogeneity in normal and tumor cells. Because short eccDNAs are poorly chromatinized, we hypothesized that they are sequenced by tagmentation in ATAC-seq experiments without any enrichment of circular DNA. Indeed, ATAC-seq identified thousands of eccDNAs in cell lines that were validated by inverse PCR and by metaphase FISH. ATAC-seq in gliomas and glioblastomas identify hundreds of eccDNAs, including one containing the well-known EGFR gene amplicon from chr7. More than 18,000 eccDNAs, many carrying known cancer driver genes, are identified in a pan-cancer analysis of ATAC-seq libraries from 23 tumor types. Somatically mosaic eccDNAs are identified by ATAC-seq even before amplification is recognized by genome-wide copy number variation measurements. Thus, ATAC-seq is a sensitive method to detect eccDNA present in a tumor at the pre-amplification stage and can be used to predict resistance to therapy.

Are Parallel Proliferation Pathways Redundant?

Are the receptor tyrosine kinase (RTK) and JAK-STAT-driven proliferation pathways 'parallel' or 'redundant'? And what about those of K-Ras4B versus N-Ras? 'Parallel' proliferation pathways accomplish a similar drug resistance outcome. Thus, are they 'redundant'? In this paper, it is argued that there is a fundamental distinction between 'parallel' and 'redundant'. Cellular proliferation pathways are influenced by the genome sequence, 3D organization and chromatin accessibility, and determined by protein availability prior to cancer emergence. In the opinion presented, if they operate the same downstream protein families, they are redundant; if evolutionary-independent, they are parallel. Thus, RTK and JAK-STAT-driven proliferation pathways are parallel; those of Ras isoforms are redundant. Our Precision Medicine Call to map cancer proliferation pathways is vastly important since it can expedite effective therapeutics.

Lineage-Specific Epigenomic and Genomic Activation of Oncogene HNF4A Promotes Gastrointestinal Adenocarcinomas

Gastrointestinal adenocarcinomas (GIAC) of the tubular gastrointestinal (GI) tract including esophagus, stomach, colon, and rectum comprise most GI cancers and share a spectrum of genomic features. However, the unified epigenomic changes specific to GIAC are poorly characterized. Using 907 GIAC samples from The Cancer Genome Atlas, we applied mathematical algorithms to large-scale DNA methylome and transcriptome profiles to reconstruct transcription factor (TF) networks and identify a list of functionally hyperactive master regulator (MR) TF shared across different GIAC. The top candidate HNF4A exhibited prominent genomic and epigenomic activation in a GIAC-specific manner. A complex interplay between the HNF4A promoter and three distal enhancer elements was coordinated by GIAC-specific MRTF including ELF3, GATA4, GATA6, and KLF5. HNF4A also self-regulated its own promoter and enhancers. Functionally, HNF4A promoted cancer proliferation and survival by transcriptional activation of many downstream targets, including HNF1A and factors of interleukin signaling, in a lineage-specific manner. Overall, our study provides new insights into the GIAC-specific gene regulatory networks and identifies potential therapeutic strategies against these common cancers. SIGNIFICANCE: These findings show that GIAC-specific master regulatory transcription factors control HNF4A via three distal enhancers to promote GIAC cell proliferation and survival. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/80/13/2722/F1.large.jpg.

ATAC-seq normalization method can significantly affect differential accessibility analysis and interpretation

BACKGROUND

Chromatin dysregulation is associated with developmental disorders and cancer. Numerous methods for measuring genome-wide chromatin accessibility have been developed in the genomic era to interrogate the function of chromatin regulators. A recent technique which has gained widespread use due to speed and low input requirements with native chromatin is the Assay for Transposase-Accessible Chromatin, or ATAC-seq. Biologists have since used this method to compare chromatin accessibility between two cellular conditions. However, approaches for calculating differential accessibility can yield conflicting results, and little emphasis is placed on choice of normalization method during differential ATAC-seq analysis, especially when global chromatin alterations might be expected.

RESULTS

Using an in vivo ATAC-seq data set generated in our recent report, we observed differences in chromatin accessibility patterns depending on the data normalization method used to calculate differential accessibility. This observation was further verified on published ATAC-seq data from yeast. We propose a generalized workflow for differential accessibility analysis using ATAC-seq data. We further show this workflow identifies sites of differential chromatin accessibility that correlate with gene expression and is sensitive to differential analysis using negative controls.

CONCLUSIONS

We argue that researchers should systematically compare multiple normalization methods before continuing with differential accessibility analysis. ATAC-seq users should be aware of the interpretations of potential bias within experimental data and the assumptions of the normalization method implemented.

Leveraging mouse chromatin data for heritability enrichment informs common disease architecture and reveals cortical layer contributions to schizophrenia

Genome-wide association studies have implicated thousands of noncoding variants across common human phenotypes. However, they cannot directly inform the cellular context in which disease-associated variants act. Here, we use open chromatin profiles from discrete mouse cell populations to address this challenge. We applied stratified linkage disequilibrium score regression and evaluated heritability enrichment in 64 genome-wide association studies, emphasizing schizophrenia. We provide evidence that mouse-derived human open chromatin profiles can serve as powerful proxies for difficult to obtain human cell populations, facilitating the illumination of common disease heritability enrichment across an array of human phenotypes. We demonstrate that signatures from discrete subpopulations of cortical excitatory and inhibitory neurons are significantly enriched for schizophrenia heritability with maximal enrichment in cortical layer V excitatory neurons. We also show that differences between schizophrenia and bipolar disorder are concentrated in excitatory neurons in cortical layers II-III, IV, and V, as well as the dentate gyrus. Finally, we leverage these data to fine-map variants in 177 schizophrenia loci nominating variants in 104/177. We integrate these data with transcription factor binding site, chromatin interaction, and validated enhancer data, placing variants in the cellular context where they may modulate risk.

Hotspots of Aberrant Enhancer Activity in Fibrolamellar Carcinoma Reveal Candidate Oncogenic Pathways and Therapeutic Vulnerabilities

Fibrolamellar carcinoma (FLC) is a rare, therapeutically intractable liver cancer that disproportionately affects youth. Although FLC tumors exhibit a distinct gene expression profile, the chromatin regulatory landscape and the genes most critical for tumor cell survival remain unclear. Here, we use chromatin run-on sequencing to discover ∼7,000 enhancers and 141 enhancer hotspots activated in FLC relative to nonmalignant liver. Bioinformatic analyses reveal aberrant ERK/MEK signaling and candidate master transcriptional regulators. We also define the genes most strongly associated with hotspots of FLC enhancer activity, including CA12 and SLC16A14. Treatment of FLC cell models with inhibitors of CA12 or SLC16A14 independently reduce cell viability and/or significantly enhance the effect of the MEK inhibitor cobimetinib. These findings highlight molecular targets for drug development, as well as drug combination approaches.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

Characterization of the chromatin accessibility in an Alzheimer's disease (AD) mouse model

Background

The pathological hallmarks of Alzheimer’s disease (AD) involve alterations in the expression of numerous genes associated with transcriptional levels, which are determined by chromatin accessibility. Here, the landscape of chromatin accessibility was studied to understand the outline of the transcription and expression of AD-associated metabolism genes in an AD mouse model.

Methods

The assay for transposase-accessible chromatin by sequencing (ATAC-seq) was used to investigate the AD-associated chromatin reshaping in the APPswe/PS1dE9 (APP/PS1) mouse model. ATAC-seq data in the hippocampus of 8-month-old APP/PS1 mice were generated, and the relationship between chromatin accessibility and gene expression was analyzed in combination with RNA sequencing. Gene ontology (GO) analysis was applied to elucidate biological processes and signaling pathways altered in APP/PS1 mice. Critical transcription factors were identified; alterations in chromatin accessibility were further confirmed using chromatin immunoprecipitation assays.

Results

We identified 1690 increased AD-associated chromatin-accessible regions in the hippocampal tissues of APP/PS1 mice. These regions were enriched in genes related to diverse signaling pathways, including the PI3K-Akt, Hippo, TGF-β, and Jak-Stat signaling pathways, which play essential roles in regulating cell proliferation, apoptosis, and inflammatory responses. A total of 1003 decreased chromatin-accessible regions were considered to be related with declined AD-associated biological processes including cellular response to hyperoxia and insulin stimulus, synaptic transmission, and positive regulation of autophagy. In the APP/PS1 hippocampus, 1090 genes were found to be upregulated and 1081 downregulated. Interestingly, enhanced ATAC-seq signal was found in approximately 740 genes, with 43 exhibiting upregulated mRNA levels. Several genes involved in AD development were found to have a significantly increased expression in APP/PS1 mice compared to controls, including Sele, Clec7a, Cst7, and Ccr6. The signatures of numerous transcription factors, including Olig2, NeuroD1, TCF4, and NeuroG2, were found enriched in the AD-associated accessible chromatin regions. The transcription-activating marks of H3K4me3 and H3K27ac were also found increased in the promoters of these genes. These results indicate that the mechanism for the upregulation of genes could be attributed to the enrichment of open chromatin regions with transcription factors motifs and the histone marks H3K4me3 and H3K27ac.

Conclusion

Our study reveals that alterations in chromatin accessibility may be an initial mechanism in AD pathogenesis.

Supplementary information

Supplementary information accompanies this paper at 10.1186/s13195-020-00598-2.

Chromatin structure regulates cancer-specific alternative splicing events in primary HPV-related oropharyngeal squamous cell carcinoma

Human papillomavirus-related oropharyngeal squamous cell carcinoma (HPV+ OPSCC) represents a unique disease entity within head and neck cancer with rising incidence. Previous work has shown that alternative splicing events (ASEs) are prevalent in HPV+ OPSCC, but further validation is needed to understand the regulation of this process and its role in these tumours. In this study, eleven ASEs (GIT2, CTNNB1, MKNK2, MRPL33, SIPA1L3, SNHG6, SYCP2, TPRG1, ZHX2, ZNF331, and ELOVL1) were selected for validation from 109 previously published candidate ASEs to elucidate the post-transcriptional mechanisms of oncogenesis in HPV+ disease. In vitro qRT-PCR confirmed differential expression of 9 of 11 ASE candidates, and in silico analysis within the TCGA cohort confirmed 8 of 11 candidates. Six ASEs (MRPL33, SIPA1L3, SNHG6, TPRG1, ZHX2, and ELOVL1) showed significant differential expression across both methods. Further evaluation of chromatin modification revealed that ASEs strongly correlated with cancer-specific distribution of acetylated lysine 27 of histone 3 (H3K27ac). Subsequent epigenetic treatment of HPV+ HNSCC cell lines (UM-SCC-047 and UPCI-SCC-090) with JQ1 not only induced downregulation of cancer-specific ASE isoforms, but also growth inhibition in both cell lines. The UPCI-SCC-090 cell line, with greater ASE expression, also showed more significant growth inhibition after JQ1 treatment. This study confirms several novel cancer-specific ASEs in HPV+OPSCC and provides evidence for the role of chromatin modifications in regulation of alternative splicing in HPV+OPSCC. This highlights the role of epigenetic changes in the oncogenesis of HPV+OPSCC, which represents a unique, unexplored target for therapeutics that can alter the global post-transcriptional landscape.

The tubal epigenome - An emerging target for ovarian cancer

Ovarian cancer is the most lethal gynecologic malignancy in the United States. The mortality of this disease is primarily attributed to challenges in early detection and therapeutic resistance. Recent studies indicate that the majority of high-grade serous ovarian carcinomas (HGSCs) originate from aberrant fallopian tube epithelial (FTE) cells. This shift in thinking about ovarian cancer pathogenesis has been met with an effort to identify the early genetic and epigenetic changes that underlie the transformation of normal FTE cells and prompt them to migrate and colonize the ovary, ultimately giving rise to aggressive HGSC. While identification of these early changes is important for biomarker discovery, the emergence of epigenetic alterations in FTE chromatin may also provide new opportunities for early detection, prevention, and therapeutic intervention. Here we provide a comprehensive overview of the current knowledge regarding early epigenetic reprogramming that precedes HGSC tumor development, the way that these alterations affect both intrinsic and extrinsic tumor properties, and how the epigenome may be targeted to thwart HGSC tumorigenesis.

3D ATAC-PALM: super-resolution imaging of the accessible genome

To image the accessible genome at nanometer scale in situ, we developed 3D ATAC-PALM which integrates Assay for Transposase-Accessible Chromatin with visualization, PALM super-resolution imaging and lattice light-sheet microscopy. Multiplexed with Oligopaint DNA-FISH, RNA-FISH and protein fluorescence, 3D ATAC-PALM connected microscopy and genomic data, revealing spatially-segregated accessible chromatin domains (ACDs) that enclose active chromatin and transcribed genes. Using these methods to analyze genetically perturbed cells, we demonstrated that genome architectural protein CTCF prevents excessive clustering of accessible chromatin and decompacts ACDs. These results highlight 3D ATAC-PALM as a useful tool to probe the structure and organizing mechanism of the genome.

TERT promoter mutations and telomeres during tumorigenesis

Telomerase regulation and telomere shortening act as a strong tumor suppressor mechanism in human somatic cells. Point mutations in the promoter of telomerase reverse transcriptase (TERT) are the most frequent non-coding mutation in cancer. These TERT promoter mutations (TPMs) create de novo ETS factor binding sites upstream of the start codon of the gene, which can be bound by different ETS factors. TPMs can occur early during tumorigenesis and are thought to be among the first mutations in melanoma, glioblastoma and hepatocellular carcinoma. Despite their association with increased TERT levels, TPMs do not prohibit telomere shortening and TPM-harboring cancers present with short telomeres. Their short telomere length combined with their high prevalence and specificity for cancer makes TPMs an attractive target for future therapeutic exploitation of telomerase inhibition and telomere deprotection-induced cell death.

Specific chromatin landscapes and transcription factors couple breast cancer subtype with metastatic relapse to lung or brain

Background

Few somatic mutations have been linked to breast cancer metastasis, whereas transcriptomic differences among primary tumors correlate with incidence of metastasis, especially to the lungs and brain. However, the epigenomic alterations and transcription factors (TFs) which underlie these alterations remain unclear.

Methods

To identify these, we performed RNA-seq, Chromatin Immunoprecipitation and sequencing (ChIP-seq) and Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) of the MDA-MB-231 cell line and its brain (BrM2) and lung (LM2) metastatic sub-populations. We incorporated ATAC-seq data from TCGA to assess metastatic open chromatin signatures, and gene expression data from human metastatic datasets to nominate transcription factor biomarkers.

Results

Our integrated epigenomic analyses found that lung and brain metastatic cells exhibit both shared and distinctive signatures of active chromatin. Notably, metastatic sub-populations exhibit increased activation of both promoters and enhancers. We also integrated these data with chromosome conformation capture coupled with ChIP-seq (HiChIP) derived enhancer-promoter interactions to predict enhancer-controlled pathway alterations. We found that enhancer changes are associated with endothelial cell migration in LM2, and negative regulation of epithelial cell proliferation in BrM2. Promoter changes are associated with vasculature development in LM2 and homophilic cell adhesion in BrM2. Using ATAC-seq, we identified a metastasis open-chromatin signature that is elevated in basal-like and HER2-enriched breast cancer subtypes and associates with worse prognosis in human samples. We further uncovered TFs associated with the open chromatin landscapes of metastatic cells and whose expression correlates with risk for metastasis. While some of these TFs are associated with primary breast tumor subtypes, others more specifically correlate with lung or brain metastasis.

Conclusions

We identify distinctive epigenomic properties of breast cancer cells that metastasize to the lung and brain. We also demonstrate that signatures of active chromatin sites are partially linked to human breast cancer subtypes with poor prognosis, and that specific TFs can independently distinguish lung and brain relapse.

Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations

Identifying cancer-relevant mutations in noncoding regions is challenging due to the large numbers of such mutations, their low levels of recurrence, and difficulties in interpreting their functional impact. To uncover genes that are dysregulated due to somatic mutations in cis, we build upon the concept of differential allele-specific expression (ASE) and introduce methods to identify genes within an individual's cancer whose ASE differs from what is found in matched normal tissue. When applied to breast cancer tumor samples, our methods detect the known allele-specific effects of copy number variation and nonsense-mediated decay. Further, genes that are found to recurrently exhibit differential ASE across samples are cancer relevant. Genes with cis mutations are enriched for differential ASE, and we find 147 potentially functional noncoding mutations cis to genes that exhibit significant differential ASE. We conclude that differential ASE is a promising means for discovering gene dysregulation due to cis noncoding mutations.

Comprehensive pharmacogenomic characterization of gastric cancer

BACKGROUND

Gastric cancer is among the most lethal human malignancies. Previous studies have identified molecular aberrations that constitute dynamic biological networks and genomic complexities of gastric tumors. However, the clinical translation of molecular-guided targeted therapy is hampered by challenges. Notably, solid tumors often harbor multiple genetic alterations, complicating the development of effective treatments.

METHODS

To address such challenges, we established a comprehensive dataset of molecularly annotated patient derivatives coupled with pharmacological profiles for 60 targeted agents to explore dynamic pharmacogenomic interactions in gastric cancers.

RESULTS

We identified lineage-specific drug sensitivities based on histopathological and molecular subclassification, including substantial sensitivities toward VEGFR and EGFR inhibition therapies in diffuse- and signet ring-type gastric tumors, respectively. We identified potential therapeutic opportunities for WNT pathway inhibitors in ALK-mutant tumors, a significant association between PIK3CA-E542K mutation and AZD5363 response, and transcriptome expression of RNF11 as a potential predictor of response to gefitinib.

CONCLUSIONS

Collectively, our results demonstrate the feasibility of drug screening combined with tumor molecular characterization to facilitate personalized therapeutic regimens for gastric tumors.

epiCOLOC: Integrating Large-Scale and Context-Dependent Epigenomics Features for Comprehensive Colocalization Analysis

High-throughput genome-wide epigenomic assays, such as ChIP-seq, DNase-seq and ATAC-seq, have profiled a huge number of functional elements across numerous human tissues/cell types, which provide an unprecedented opportunity to interpret human genome and disease in context-dependent manner. Colocalization analysis determines whether genomic features are functionally related to a given search and will facilitate identifying the underlying biological functions characterizing intricate relationships with queries for genomic regions. Existing colocalization methods leveraged diverse assumptions and background models to assess the significance of enrichment, however, they only provided limited and predefined sets of epigenomic features. Here, we comprehensively collected and integrated over 44,385 bulk or single-cell epigenomic assays across 53 human tissues/cell types, such as transcription factor binding, histone modification, open chromatin and transcriptional event. By classifying these profiles into hierarchy of tissue/cell type, we developed a web portal, epiCOLOC (http://mulinlab.org/epicoloc or http://mulinlab.tmu.edu.cn/epicoloc), for users to perform context-dependent colocalization analysis in a convenient way.

TRlnc: a comprehensive database for human transcriptional regulatory information of lncRNAs

Long noncoding RNAs (lncRNAs) have been proven to play important roles in transcriptional processes and biological functions. With the increasing study of human diseases and biological processes, information in human H3K27ac ChIP-seq, ATAC-seq and DNase-seq datasets is accumulating rapidly, resulting in an urgent need to collect and process data to identify transcriptional regulatory regions of lncRNAs. We therefore developed a comprehensive database for human regulatory information of lncRNAs (TRlnc, http://bio.licpathway.net/TRlnc), which aimed to collect available resources of transcriptional regulatory regions of lncRNAs and to annotate and illustrate their potential roles in the regulation of lncRNAs in a cell type-specific manner. The current version of TRlnc contains 8 683 028 typical enhancers/super-enhancers and 32 348 244 chromatin accessibility regions associated with 91 906 human lncRNAs. These regions are identified from over 900 human H3K27ac ChIP-seq, ATAC-seq and DNase-seq samples. Furthermore, TRlnc provides the detailed genetic and epigenetic annotation information within transcriptional regulatory regions (promoter, enhancer/super-enhancer and chromatin accessibility regions) of lncRNAs, including common SNPs, risk SNPs, eQTLs, linkage disequilibrium SNPs, transcription factors, methylation sites, histone modifications and 3D chromatin interactions. It is anticipated that the use of TRlnc will help users to gain in-depth and useful insights into the transcriptional regulatory mechanisms of lncRNAs.

From reads to insight: a hitchhiker's guide to ATAC-seq data analysis

Assay of Transposase Accessible Chromatin sequencing (ATAC-seq) is widely used in studying chromatin biology, but a comprehensive review of the analysis tools has not been completed yet. Here, we discuss the major steps in ATAC-seq data analysis, including pre-analysis (quality check and alignment), core analysis (peak calling), and advanced analysis (peak differential analysis and annotation, motif enrichment, footprinting, and nucleosome position analysis). We also review the reconstruction of transcriptional regulatory networks with multiomics data and highlight the current challenges of each step. Finally, we describe the potential of single-cell ATAC-seq and highlight the necessity of developing ATAC-seq specific analysis tools to obtain biologically meaningful insights.

CRISPR Screens Identify Essential Cell Growth Mediators in BRAF Inhibitor-resistant Melanoma

BRAF is a serine/threonine kinase that harbors activating mutations in ∼7% of human malignancies and ∼60% of melanomas. Despite initial clinical responses to BRAF inhibitors, patients frequently develop drug resistance. To identify candidate therapeutic targets for BRAF inhibitor resistant melanoma, we conduct CRISPR screens in melanoma cells harboring an activating BRAF mutation that had also acquired resistance to BRAF inhibitors. To investigate the mechanisms and pathways enabling resistance to BRAF inhibitors in melanomas, we integrate expression, ATAC-seq, and CRISPR screen data. We identify the JUN family transcription factors and the ETS family transcription factor ETV5 as key regulators of CDK6, which together enable resistance to BRAF inhibitors in melanoma cells. Our findings reveal genes contributing to resistance to a selective BRAF inhibitor PLX4720, providing new insights into gene regulation in BRAF inhibitor resistant melanoma cells.

Chromatin accessibility analysis reveals regulatory dynamics of developing human retina and hiPSC-derived retinal organoids

Retinal organoids (ROs) derived from human induced pluripotent stem cells (hiPSCs) provide potential opportunities for studying human retinal development and disorders; however, to what extent ROs recapitulate the epigenetic features of human retinal development is unknown. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics over long-term human retinal and RO development. Our results showed that ROs recapitulated the human retinogenesis to a great extent, but divergent chromatin features were also discovered. We further reconstructed the transcriptional regulatory network governing human and RO retinogenesis in vivo. Notably, NFIB and THRA were identified as regulators in human retinal development. The chromatin modifications between developing human and mouse retina were also cross-analyzed. Notably, we revealed an enriched bivalent modification of H3K4me3 and H3K27me3 in human but not in murine retinogenesis, suggesting a more dedicated epigenetic regulation on human genome.

Genomic basis for RNA alterations in cancer

Transcript alterations often result from somatic changes in cancer genomes1. Various forms of RNA alterations have been described in cancer, including overexpression2, altered splicing3 and gene fusions4; however, it is difficult to attribute these to underlying genomic changes owing to heterogeneity among patients and tumour types, and the relatively small cohorts of patients for whom samples have been analysed by both transcriptome and whole-genome sequencing. Here we present, to our knowledge, the most comprehensive catalogue of cancer-associated gene alterations to date, obtained by characterizing tumour transcriptomes from 1,188 donors of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)5. Using matched whole-genome sequencing data, we associated several categories of RNA alterations with germline and somatic DNA alterations, and identified probable genetic mechanisms. Somatic copy-number alterations were the major drivers of variations in total gene and allele-specific expression. We identified 649 associations of somatic single-nucleotide variants with gene expression in cis, of which 68.4% involved associations with flanking non-coding regions of the gene. We found 1,900 splicing alterations associated with somatic mutations, including the formation of exons within introns in proximity to Alu elements. In addition, 82% of gene fusions were associated with structural variants, including 75 of a new class, termed 'bridged' fusions, in which a third genomic location bridges two genes. We observed transcriptomic alteration signatures that differ between cancer types and have associations with variations in DNA mutational signatures. This compendium of RNA alterations in the genomic context provides a rich resource for identifying genes and mechanisms that are functionally implicated in cancer.

The Length and Distribution of Plasma Cell-Free DNA Fragments in Stroke Patients

A number of studies have shown that plasma cell-free DNA is closely related to the risk of stroke, but the fragmentation status of plasma cell-free DNA and its clinical application value in ischemic stroke are still unclear. In this study, 48 patients with new ischemic stroke and 20 healthy subjects were enrolled. The second-generation high-throughput sequencing technique was used to study the plasma cell-free fragment length and regional distribution of the subjects. As noted in our results, the ratio of plasma cell-free DNA fragments in the disease group was significantly greater than that of the healthy group in the 300–400 bp range; conversely for fragments at the 75–250 bp range, the ratio of plasma cell-free DNA fragments in the patient group was apparently lower than that of the healthy group. In-depth analysis of the proportion of fragments distributed on each component of the genome was carried out. Our results recorded that the plasma cell-free DNA fragments in the disease group were inclined to the EXON, CpG islands, and ALU regions in contrast to that of the healthy group. In particular, fragments within the 300–400 bp range of the disease group were enrichment in the regions of EXON, INTRON, INTERGENIC, LINE, Fragile, ALU, and CpG islands. In summary, our findings suggested that the intracellular DNA degradation profiles could be applied to distinguish the stroke group and the healthy group, which provided a theoretical basis for the clinical diagnosis and prognosis of stroke by profiling the characteristic of plasma cell-free DNA fragments.

Chromatin accessibility dynamics in a model of human forebrain development

Forebrain development is characterized by highly synchronized cellular processes, which, if perturbed, can cause disease. To chart the regulatory activity underlying these events, we generated a map of accessible chromatin in human three-dimensional forebrain organoids. To capture corticogenesis, we sampled glial and neuronal lineages from dorsal or ventral forebrain organoids over 20 months in vitro. Active chromatin regions identified in human primary brain tissue were observed in organoids at different developmental stages. We used this resource to map genetic risk for disease and to explore evolutionary conservation. Moreover, we integrated chromatin accessibility with transcriptomics to identify putative enhancer-gene linkages and transcription factors that regulate human corticogenesis. Overall, this platform brings insights into gene-regulatory dynamics at previously inaccessible stages of human forebrain development, including signatures of neuropsychiatric disorders.

Enhancer RNAs in cancer: regulation, mechanisms and therapeutic potential

Enhancers are distal genomic elements critical for gene regulation and cell identify control during development and diseases. Many human cancers were found to associate with enhancer malfunction, due to genetic and epigenetic alterations, which in some cases directly drive tumour growth. Conventionally, enhancers are known to provide DNA binding motifs to recruit transcription factors (TFs) and to control target genes. However, recent progress found that most, if not all, active enhancers pervasively transcribe noncoding RNAs that are referred to as enhancer RNAs (eRNAs). Increasing evidence points to functional roles of at least a subset of eRNAs in gene regulation in both normal and cancer cells, adding new insights into the action mechanisms of enhancers. eRNA expression was observed to be widespread but also specific to tumour types and individual patients, serving as opportunities to exploit them as potential diagnosis markers or therapeutic targets. In this review, we discuss the brief history of eRNA research, their functional mechanisms and importance in cancer gene regulation, as well as their therapeutic and diagnostic values in cancer. We propose that further studies of eRNAs in cancer will offer a promising 'eRNA targeted therapy' for human cancer intervention.

Candidate Cancer Driver Mutations in Distal Regulatory Elements and Long-Range Chromatin Interaction Networks

A comprehensive catalog of cancer driver mutations is essential for understanding tumorigenesis and developing therapies. Exome-sequencing studies have mapped many protein-coding drivers, yet few non-coding drivers are known because genome-wide discovery is challenging. We developed a driver discovery method, ActiveDriverWGS, and analyzed 120,788 cis-regulatory modules (CRMs) across 1,844 whole tumor genomes from the ICGC-TCGA PCAWG project. We found 30 CRMs with enriched SNVs and indels (FDR < 0.05). These frequently mutated regulatory elements (FMREs) were ubiquitously active in human tissues, showed long-range chromatin interactions and mRNA abundance associations with target genes, and were enriched in motif-rewiring mutations and structural variants. Genomic deletion of one FMRE in human cells caused proliferative deficiencies and transcriptional deregulation of cancer genes CCNB1IP1, CDH1, and CDKN2B, validating observations in FMRE-mutated tumors. Pathway analysis revealed further sub-significant FMREs at cancer genes and processes, indicating an unexplored landscape of infrequent driver mutations in the non-coding genome.

CATA: a comprehensive chromatin accessibility database for cancer

Accessible chromatin refers to the active regions of a chromosome that are bound by many transcription factors (TFs). Changes in chromatin accessibility play a critical role in tumorigenesis. With the emergence of novel methods like Assay for Transposase-accessible Chromatin Sequencing, a sequencing method that maps chromatin-accessible regions (CARs) and enables the computational analysis of TF binding at chromatin-accessible sites, the regulatory landscape in cancer can be dissected. Herein, we developed a comprehensive cancer chromatin accessibility database named CATA, which aims to provide available resources of cancer CARs and to annotate their potential roles in the regulation of genes in a cancer type-specific manner. In this version, CATA stores 2 991 163 CARs from 23 cancer types, binding information of 1398 TFs within the CARs, and provides multiple annotations about these regions, including common single nucleotide polymorphisms (SNPs), risk SNPs, copy number variation, somatic mutations, motif changes, expression quantitative trait loci, methylation and CRISPR/Cas9 target loci. Moreover, CATA supports cancer survival analysis of the CAR-associated genes and provides detailed clinical information of the tumor samples. Database URL: CATA is available at http://www.xiejjlab.bio/cata/.

A genome-wide view of the de-differentiation of central nervous system endothelial cells in culture

Vascular endothelial cells (ECs) derived from the central nervous system (CNS) variably lose their unique barrier properties during in vitro culture, hindering the development of robust assays for blood-brain barrier (BBB) function, including drug permeability and extrusion assays. In previous work (Sabbagh et al., 2018) we characterized transcriptional and accessible chromatin landscapes of acutely isolated mouse CNS ECs. In this report, we compare transcriptional and accessible chromatin landscapes of acutely isolated mouse CNS ECs versus mouse CNS ECs in short-term in vitro culture. We observe that standard culture conditions are associated with a rapid and selective loss of BBB transcripts and chromatin features, as well as a greatly reduced level of beta-catenin signaling. Interestingly, forced expression of a stabilized derivative of beta-catenin, which in vivo leads to a partial conversion of non-BBB CNS ECs to a BBB-like state, has little or no effect on gene expression or chromatin accessibility in vitro.

Interpreting pathways to discover cancer driver genes with Moonlight

Cancer driver gene alterations influence cancer development, occurring in oncogenes, tumor suppressors, and dual role genes. Discovering dual role cancer genes is difficult because of their elusive context-dependent behavior. We define oncogenic mediators as genes controlling biological processes. With them, we classify cancer driver genes, unveiling their roles in cancer mechanisms. To this end, we present Moonlight, a tool that incorporates multiple -omics data to identify critical cancer driver genes. With Moonlight, we analyze 8000+ tumor samples from 18 cancer types, discovering 3310 oncogenic mediators, 151 having dual roles. By incorporating additional data (amplification, mutation, DNA methylation, chromatin accessibility), we reveal 1000+ cancer driver genes, corroborating known molecular mechanisms. Additionally, we confirm critical cancer driver genes by analysing cell-line datasets. We discover inactivation of tumor suppressors in intron regions and that tissue type and subtype indicate dual role status. These findings help explain tumor heterogeneity and could guide therapeutic decisions.

Molecular mechanisms of epigenetic inheritance: Possible evolutionary implications

Recently interest in multi-generational epigenetic phenomena have been fuelled by highly reproducible intergenerational and transgenerational inheritance paradigms in several model organisms. Such paradigms are essential in order to begin to use genetics to unpick the mechanistic bases of how epigenetic information may be transmitted between generations; indeed great strides have been made towards understanding these mechanisms. Far less well understood is the relationship between epigenetic inheritance, ecology and evolution. In this review I focus on potential connections between laboratory studies of transgenerational epigenetic inheritance phenomena and evolutionary processes that occur in natural populations. In the first section, I consider whether transgenerational epigenetic inheritance might provide an advantage to organisms over the short term in adapting to their environment. Second, I consider whether epigenetic changes can contribute to the evolution of species by contributing to stable phenotypic variation within a population. Finally I discuss whether epigenetic changes could influence evolution by either directly or indirectly promoting DNA sequence changes that could impact phenotypic divergence. Additionally, I will discuss how epigenetic changes could influence the evolution of human cancer and thus be directly relevant for the development of this disease.

The Promises and Challenges of Toxico-Epigenomics: Environmental Chemicals and Their Impacts on the Epigenome

BACKGROUND

It has been estimated that a substantial portion of chronic and noncommunicable diseases can be caused or exacerbated by exposure to environmental chemicals. Multiple lines of evidence indicate that early life exposure to environmental chemicals at relatively low concentrations could have lasting effects on individual and population health. Although the potential adverse effects of environmental chemicals are known to the scientific community, regulatory agencies, and the public, little is known about the mechanistic basis by which these chemicals can induce long-term or transgenerational effects. To address this question, epigenetic mechanisms have emerged as the potential link between genetic and environmental factors of health and disease.

OBJECTIVES

We present an overview of epigenetic regulation and a summary of reported evidence of environmental toxicants as epigenetic disruptors. We also discuss the advantages and challenges of using epigenetic biomarkers as an indicator of toxicant exposure, using measures that can be taken to improve risk assessment, and our perspectives on the future role of epigenetics in toxicology.

DISCUSSION

Until recently, efforts to apply epigenomic data in toxicology and risk assessment were restricted by an incomplete understanding of epigenomic variability across tissue types and populations. This is poised to change with the development of new tools and concerted efforts by researchers across disciplines that have led to a better understanding of epigenetic mechanisms and comprehensive maps of epigenomic variation. With the foundations now in place, we foresee that unprecedented advancements will take place in the field in the coming years. https://doi.org/10.1289/EHP6104.

Epigenetics Analysis and Integrated Analysis of Multiomics Data, Including Epigenetic Data, Using Artificial Intelligence in the Era of Precision Medicine

To clarify the mechanisms of diseases, such as cancer, studies analyzing genetic mutations have been actively conducted for a long time, and a large number of achievements have already been reported. Indeed, genomic medicine is considered the core discipline of precision medicine, and currently, the clinical application of cutting-edge genomic medicine aimed at improving the prevention, diagnosis and treatment of a wide range of diseases is promoted. However, although the Human Genome Project was completed in 2003 and large-scale genetic analyses have since been accomplished worldwide with the development of next-generation sequencing (NGS), explaining the mechanism of disease onset only using genetic variation has been recognized as difficult. Meanwhile, the importance of epigenetics, which describes inheritance by mechanisms other than the genomic DNA sequence, has recently attracted attention, and, in particular, many studies have reported the involvement of epigenetic deregulation in human cancer. So far, given that genetic and epigenetic studies tend to be accomplished independently, physiological relationships between genetics and epigenetics in diseases remain almost unknown. Since this situation may be a disadvantage to developing precision medicine, the integrated understanding of genetic variation and epigenetic deregulation appears to be now critical. Importantly, the current progress of artificial intelligence (AI) technologies, such as machine learning and deep learning, is remarkable and enables multimodal analyses of big omics data. In this regard, it is important to develop a platform that can conduct multimodal analysis of medical big data using AI as this may accelerate the realization of precision medicine. In this review, we discuss the importance of genome-wide epigenetic and multiomics analyses using AI in the era of precision medicine.

Characterization of cancer omics and drug perturbations in panels of lung cancer cells

To better understand the disruptions of transcriptional regulations and gene expression in lung cancers, we constructed a multi-omics catalogue of the responses of lung cancer cells to a series of chemical compounds. We generated and analyzed 3,240 RNA-seq and 3,393 ATAC-seq libraries obtained from 23 cell lines treated with 95 well-annotated compounds. To demonstrate the power of the created multi-omics resource, we attempted to identify drugs that could induce the designated changes alone or in combination. The basal multi-omics information was first integrated into co-expression modules. Among these modules, we identified a stress response module that may be a promising drug intervention target, as new combinations of compounds that could be used to regulate this module and the consequent phenotypic appearance of cancer cells have been identified. We believe that the multi-omics profiles generated in this study and the strategy used to stratify them will lead to more rational and efficient development of anticancer drugs.

Which came first, the transcriptional regulator or its target genes? An evolutionary perspective into the construction of eukaryotic regulons

Eukaryotic regulons are regulatory units formed by a set of genes under the control of the same transcription factor (TF). Despite the functional plasticity, TFs are highly conserved and recognize the same DNA sequences in different organisms. One of the main factors that confer regulatory specificity is the distribution of the binding sites of the TFs along the genome, allowing the configuration of different transcriptional regulatory networks (TRNs) from the same regulator. A similar scenario occurs between tissues of the same organism, where a TRN can be rewired by epigenetic factors, modulating the accessibility of the TF to its binding sites. In this article we discuss concepts that can help to formulate testable hypotheses about the construction of regulons, exploring the presence and absence of the elements that form a TRN throughout the evolution of an ancestral lineage. This article is part of a Special Issue entitled: Transcriptional Profiles and Regulatory Gene Networks edited by Dr. Federico Manuel Giorgi and Dr. Shaun Mahony.

Single-nucleus chromatin accessibility reveals intratumoral epigenetic heterogeneity in IDH1 mutant gliomas

The presence of genome-wide DNA hypermethylation is a hallmark of lower grade gliomas (LGG) with isocitrate dehydrogenase (IDH) mutations. Further molecular classification of IDH mutant gliomas is defined by the presence (IDHmut-codel) or absence (IDHmut-noncodel) of hemizygous codeletion of chromosome arms 1p and 19q. Despite the DNA hypermethylation seen in bulk tumors, intra-tumoral heterogeneity at the epigenetic level has not been thoroughly analyzed. To address this question, we performed the first epigenetic profiling of single cells in a cohort of 5 gliomas with IDH1 mutation using single nucleus Assay for Transposase-Accessible Chromatin with high-throughput sequencing (snATAC-seq). Using the Fluidigm HT IFC microfluidics platform, we generated chromatin accessibility maps from 336 individual nuclei, and identified variable promoter accessibility of non-coding RNAs in LGGs. Interestingly, local chromatin structures of several non-coding RNAs are significant factors that contribute to heterogeneity, and show increased promoter accessibility in IDHmut-noncodel samples. As an example for clinical significance of this result, we identify CYTOR as a poor prognosis factor in gliomas with IDH mutation. Open chromatin assay points to differential accessibility of non-coding RNAs as an important source of epigenetic heterogeneity within individual tumors and between molecular subgroups. Rare populations of nuclei that resemble either IDH mutant molecular group co-exist within IDHmut-noncodel and IDHmut-codel groups, and along with non-coding RNAs may be an important issue to consider for future studies, as they may help guide predict treatment response and relapse.

A web-based explorer for the data is available at shiny.turcanlab.org.