Oncogenic Signaling Pathways in The Cancer Genome Atlas.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFß signaling, p53 and ß-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

Large DNA Methylation Nadirs Anchor Chromatin Loops Maintaining Hematopoietic Stem Cell Identity.

Higher-order chromatin structure and DNA methylation are implicated in multiple developmental processes, but their relationship to cell state is unknown. Here, we find that large (>7.3 kb) DNA methylation nadirs (termed "grand canyons") can form long loops connecting anchor loci that may be dozens of megabases (Mb) apart, as well as inter-chromosomal links. The interacting loci cover a total of ∼3.5 Mb of the human genome. The strongest interactions are associated with repressive marks made by the Polycomb complex and are diminished upon EZH2 inhibitor treatment. The data are suggestive of the formation of these loops by interactions between repressive elements in the loci, forming a genomic subcompartment, rather than by cohesion/CTCF-mediated extrusion. Interestingly, unlike previously characterized subcompartments, these interactions are present only in particular cell types, such as stem and progenitor cells. Our work reveals that H3K27me3-marked large DNA methylation grand canyons represent a set of very-long-range loops associated with cellular identity.

The Immune Landscape of Cancer.

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

Low-input and single-cell methods for Infinium DNA methylation BeadChips.

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

BISCUIT: an efficient, standards-compliant tool suite for simultaneous genetic and epigenetic inference in bulk and single-cell studies.

Data from both bulk and single-cell whole-genome DNA methylation experiments are under-utilized in many ways. This is attributable to inefficient mapping of methylation sequencing reads, routinely discarded genetic information, and neglected read-level epigenetic and genetic linkage information. We introduce the BISulfite-seq Command line User Interface Toolkit (BISCUIT) and its companion R/Bioconductor package, biscuiteer, for simultaneous extraction of genetic and epigenetic information from bulk and single-cell DNA methylation sequencing. BISCUIT's performance, flexibility and standards-compliant output allow large, complex experimental designs to be characterized on clinical timescales. BISCUIT is particularly suited for processing data from single-cell DNA methylation assays, with its excellent scalability, efficiency, and ability to greatly enhance mappability, a key challenge for single-cell studies. We also introduce the epiBED format for single-molecule analysis of coupled epigenetic and genetic information, facilitating the study of cellular and tissue heterogeneity from DNA methylation sequencing.

Comparative epigenome analysis using Infinium DNA methylation BeadChips.

The arrival of the Infinium DNA methylation BeadChips for mice and other nonhuman mammalian species has outpaced the development of the informatics that supports their use for epigenetics study in model organisms. Here, we present informatics infrastructure and methods to allow easy DNA methylation analysis on multiple species, including domesticated animals and inbred laboratory mice (in SeSAMe version 1.16.0+). First, we developed a data-driven analysis pipeline covering species inference, genome-specific data preprocessing and regression modeling. We targeted genomes of 310 species and 37 inbred mouse strains and showed that genome-specific preprocessing prevents artifacts and yields more accurate measurements than generic pipelines. Second, we uncovered the dynamics of the epigenome evolution in different genomic territories and tissue types through comparative analysis. We identified a catalog of inbred mouse strain-specific methylation differences, some of which are linked to the strains' immune, metabolic and neurological phenotypes. By streamlining DNA methylation array analysis for undesigned genomes, our methods extend epigenome research to broad species contexts.

mLiftOver: harmonizing data across infinium DNA methylation platforms.

MOTIVATION: Infinium DNA methylation BeadChips are widely used for genome-wide DNA methylation profiling at the population scale. Recent updates to probe content and naming conventions in the EPIC version 2 (EPICv2) arrays have complicated integrating new data with previous Infinium array platforms, such as the MethylationEPIC (EPIC) and the HumanMethylation450 (HM450) BeadChip. RESULTS: We present mLiftOver, a user-friendly tool that harmonizes probe ID, methylation level, and signal intensity data across different Infinium platforms. It manages probe replicates, missing data imputation, and platform-specific bias for accurate data conversion. We validated the tool by applying HM450-based cancer classifiers to EPICv2 cancer data, achieving high accuracy. Additionally, we successfully integrated EPICv2 healthy tissue data with legacy HM450 data for tissue identity analysis and produced consistent copy number profiles in cancer cells. AVAILABILITY: mLiftOver is implemented R and available in the Bioconductor package SeSAMe (version 1.21.13+): https://bioconductor.org/packages/release/bioc/html/sesame.html. Analysis of EPIC and EPICv2 platform-specific bias and high-confidence mapping is available at https://github.com/zhou-lab/InfiniumAnnotationV1/raw/main/Anno/EPICv2/EPICv2ToEPIC_conversion.tsv.gz. The source code is available at https://github.com/zwdzwd/sesame/blob/devel/R/mLiftOver.R under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dupsifter: a lightweight duplicate marking tool for whole genome bisulfite sequencing.

SUMMARY: In whole genome sequencing data, polymerase chain reaction amplification results in duplicate DNA fragments coming from the same location in the genome. The process of preparing a whole genome bisulfite sequencing (WGBS) library, on the other hand, can create two DNA fragments from the same location that should not be considered duplicates. Currently, only one WGBS-aware duplicate marking tool exists. However, it only works with the output from a single tool, does not accept streaming input or output, and requires a substantial amount of memory relative to the input size. Dupsifter provides an aligner-agnostic duplicate marking tool that is lightweight, has streaming capabilities, and is memory efficient. AVAILABILITY AND IMPLEMENTATION: Source code and binaries are freely available at https://github.com/huishenlab/dupsifter under the MIT license. Dupsifter is implemented in C and is supported on macOS and Linux.

EHMT2 suppresses the variation of transcriptional switches in the mouse embryo.

EHMT2 is the main euchromatic H3K9 methyltransferase. Embryos with zygotic, or maternal mutation in the Ehmt2 gene exhibit variable developmental delay. To understand how EHMT2 prevents variable developmental delay we performed RNA sequencing of mutant and somite stage-matched normal embryos at 8.5-9.5 days of gestation. Using four-way comparisons between delayed and normal embryos we clarified what it takes to be normal and what it takes to develop. We identified differentially expressed genes, for example Hox genes that simply reflected the difference in developmental progression of wild type and the delayed mutant uterus-mate embryos. By comparing wild type and zygotic mutant embryos along the same developmental window we detected a role of EHMT2 in suppressing variation in the transcriptional switches. We identified transcription changes where precise switching during development occurred only in the normal but not in the mutant embryo. At the 6-somite stage, gastrulation-specific genes were not precisely switched off in the Ehmt2-/- zygotic mutant embryos, while genes involved in organ growth, connective tissue development, striated muscle development, muscle differentiation, and cartilage development were not precisely switched on. The Ehmt2mat-/+ maternal mutant embryos displayed high transcriptional variation consistent with their variable survival. Variable derepression of transcripts occurred dominantly in the maternally inherited allele. Transcription was normal in the parental haploinsufficient wild type embryos despite their delay, consistent with their good prospects. Global profiling of transposable elements revealed EHMT2 targeted DNA methylation and suppression at LTR repeats, mostly ERVKs. In Ehmt2-/- embryos, transcription over very long distances initiated from such misregulated 'driver' ERVK repeats, encompassing a multitude of misexpressed 'passenger' repeats. In summary, EHMT2 reduced transcriptional variation of developmental switch genes and developmentally switching repeat elements at the six-somite stage embryos. These findings establish EHMT2 as a suppressor of transcriptional and developmental variation at the transition between gastrulation and organ specification.

DNA methylation enables transposable element-driven genome expansion.

Multicellular eukaryotic genomes show enormous differences in size. A substantial part of this variation is due to the presence of transposable elements (TEs). They contribute significantly to a cell's mass of DNA and have the potential to become involved in host gene control. We argue that the suppression of their activities by methylation of the C-phosphate-G (CpG) dinucleotide in DNA is essential for their long-term accommodation in the host genome and, therefore, to its expansion. An inevitable consequence of cytosine methylation is an increase in C-to-T transition mutations via deamination, which causes CpG loss. Cytosine deamination is often needed for TEs to take on regulatory functions in the host genome. Our study of the whole-genome sequences of 53 organisms showed a positive correlation between the size of a genome and the percentage of TEs it contains, as well as a negative correlation between size and the CpG observed/expected (O/E) ratio in both TEs and the host DNA. TEs are seldom found at promoters and transcription start sites, but they are found more at enhancers, particularly after they have accumulated C-to-T and other mutations. Therefore, the methylation of TE DNA allows for genome expansion and also leads to new opportunities for gene control by TE-based regulatory sites.

Switching roles for DNA and histone methylation depend on evolutionary ages of human endogenous retroviruses.

We provide a comprehensive genomic and epigenomic map of the more than 500,000 endogenous retroviruses (ERVs) and fragments that populate the intergenic regions of the human genome. The repressive epigenetic marks associated with the ERVs, particularly long terminal repeats (LTRs), show a remarkable switch in silencing mechanisms, depending on the evolutionary age of the LTRs. Young LTRs tend to be CpG rich and are mainly suppressed by DNA methylation, whereas intermediate age LTRs are associated predominantly with histone modifications, particularly histone H3 lysine 9 (H3K9) methylation. Young LTRs can be reactivated by treatment with the DNA methylation inhibitor 5-aza-2'-deoxycytidine (5-aza-CdR) alone, but their level of expression is much increased by 5-aza-CdR treatment plus knockdown of one of several H3K9 methyltransferases or of the H3K27 methyltransferase EZH2. The removal of cytosine methylation led to rapid, widespread increases in H3K9me3 in the LTRs. Intermediate age LTRs had lower CpG densities and were not up-regulated by 5-aza-CdR treatment, but they were sensitive to knockdown of H3K9 methyltransferases. Unlike the situation in embryonic stem cells, the polycomb repressive complex (PRC2) has a minor role in LTR suppression by itself and is only a player after removal of cytosine methylation in the analyzed cancer cell line. Up-regulation of LTRs and induction of "viral mimicry" is rapidly becoming of interest for predicting cancer patient response to epigenetic therapies. Understanding the mechanism for LTR suppression is of major importance in order to improve patient treatment strategies.

Mother-child transmission of epigenetic information by tunable polymorphic imprinting.

Genomic imprinting mediated by DNA methylation restricts gene expression to a single allele determined by parental origin and is not generally considered to be under genetic or environmental influence. Here, we focused on a differentially methylated region (DMR) of approximately 1.9 kb that includes a 101-bp noncoding RNA gene (nc886/VTRNA2-1), which is maternally imprinted in ∼75% of humans. This is unlike other imprinted genes, which demonstrate monoallelic methylation in 100% of individuals. The DMR includes a CTCF binding site on the centromeric side defining the DMR boundary and is flanked by a CTCF binding site on the telomeric side. The centromeric CTCF binding site contains an A/C polymorphism (rs2346018); the C allele is associated with less imprinting. The frequency of imprinting of the nc886 DMR in infants was linked to at least two nongenetic factors, maternal age at delivery and season of conception. In a separate cohort, nc886 imprinting was associated with lower body mass index in children at 5 y of age. Thus, we propose that the imprinting status of the nc886 DMR is "tunable" in that it is associated with maternal haplotype and prenatal environment. This provides a potential mechanism for transmitting information, with phenotypic consequences, from mother to child.

novoBreak: local assembly for breakpoint detection in cancer genomes.

We present novoBreak, a genome-wide local assembly algorithm that discovers somatic and germline structural variation breakpoints in whole-genome sequencing data. novoBreak consistently outperformed existing algorithms on real cancer genome data and on synthetic tumors in the ICGC-TCGA DREAM 8.5 Somatic Mutation Calling Challenge primarily because it more effectively utilized reads spanning breakpoints. novoBreak also demonstrated great sensitivity in identifying short insertions and deletions.

SeSAMe: reducing artifactual detection of DNA methylation by Infinium BeadChips in genomic deletions.

We report a new class of artifacts in DNA methylation measurements from Illumina HumanMethylation450 and MethylationEPIC arrays. These artifacts reflect failed hybridization to target DNA, often due to germline or somatic deletions and manifest as incorrectly reported intermediate methylation. The artifacts often survive existing preprocessing pipelines, masquerade as epigenetic alterations and can confound discoveries in epigenome-wide association studies and studies of methylation-quantitative trait loci. We implement a solution, P-value with out-of-band (OOB) array hybridization (pOOBAH), in the R package SeSAMe. Our method effectively masks deleted and hyperpolymorphic regions, reducing or eliminating spurious reports of epigenetic silencing at oft-deleted tumor suppressor genes such as CDKN2A and RB1 in cases with somatic deletions. Furthermore, our method substantially decreases technical variation whilst retaining biological variation, both within and across HM450 and EPIC platform measurements. SeSAMe provides a light-weight, modular DNA methylation data analysis suite, with a performant implementation suitable for efficient analysis of thousands of samples.

Integrative Epigenetic Analysis Reveals Therapeutic Targets to the DNA Methyltransferase Inhibitor Guadecitabine (SGI-110) in Hepatocellular Carcinoma.

There is an urgent need to develop more effective therapies for hepatocellular carcinoma (HCC) because of its aggressiveness. Guadecitabine (SGI-110) is a second-generation DNA methyltransferase inhibitor (DNMTi), which is currently in clinical trials for HCC and shows greater stability and performance over first-generation DNMTis. In order to identify potential therapeutic targets of SGI-110 for clinical trials, HCC cell lines (SNU398, HepG2, and SNU475) were used to evaluate the effects of transient SGI-110 treatment by an integrative analysis of DNA methylation, nucleosome accessibility, gene expression profiles, and its clinical relevance by comparison to The Cancer Genome Atlas (TCGA) HCC clinical data. Each HCC cell line represents a different DNA methylation subtype of primary HCC tumors based on TCGA data. After SGI-110 treatment, all cell lines were sensitive to SGI-110 with prolonged antiproliferation effects. Expression of up-regulated genes, including tumor suppressors, was positively correlated with nucleosome accessibility and negatively correlated with gene promoter DNA methylation. Alternatively, expression of down-regulated genes, such as oncogenes, was negatively correlated with nucleosome accessibility and positively correlated with gene body DNA methylation. SGI-110 can also act as a dual inhibitor to down-regulate polycomb repressive complex 2 (PRC2) genes by demethylating their gene bodies, resulting in reactivation of PRC2 repressed genes without involvement of DNA methylation. Furthermore, it can up-regulate endogenous retroviruses to reactivate immune pathways. Finally, about 48% of frequently altered genes in primary HCC tumors can be reversed by SGI-110 treatment. CONCLUSION: Our integrative analysis has successfully linked the antitumor effects of SGI-110 to detailed epigenetic alterations in HCC cells, identified potential therapeutic targets, and provided a rationale for combination treatments of SGI-110 with immune checkpoint therapies.

Comprehensive characterization, annotation and innovative use of Infinium DNA methylation BeadChip probes.

Illumina Infinium DNA Methylation BeadChips represent the most widely used genome-scale DNA methylation assays. Existing strategies for masking Infinium probes overlapping repeats or single nucleotide polymorphisms (SNPs) are based largely on ad hoc assumptions and subjective criteria. In addition, the recently introduced MethylationEPIC (EPIC) array expands on the utility of this platform, but has not yet been well characterized. We present in this paper an extensive characterization of probes on the EPIC and HM450 microarrays, including mappability to the latest genome build, genomic copy number of the 3΄ nested subsequence and influence of polymorphisms including a previously unrecognized color channel switch for Type I probes. We show empirical evidence for exclusion criteria for underperforming probes, providing a sounder basis than current ad hoc criteria for exclusion. In addition, we describe novel probe uses, exemplified by the addition of a total of 1052 SNP probes to the existing 59 explicit SNP probes on the EPIC array and the use of these probes to predict ethnicity. Finally, we present an innovative out-of-band color channel application for the dual use of 62 371 probes as internal bisulfite conversion controls.

Vitamin C increases viral mimicry induced by 5-aza-2'-deoxycytidine.

Vitamin C deficiency is found in patients with cancer and might complicate various therapy paradigms. Here we show how this deficiency may influence the use of DNA methyltransferase inhibitors (DNMTis) for treatment of hematological neoplasias. In vitro, when vitamin C is added at physiological levels to low doses of the DNMTi 5-aza-2'-deoxycytidine (5-aza-CdR), there is a synergistic inhibition of cancer-cell proliferation and increased apoptosis. These effects are associated with enhanced immune signals including increased expression of bidirectionally transcribed endogenous retrovirus (ERV) transcripts, increased cytosolic dsRNA, and activation of an IFN-inducing cellular response. This synergistic effect is likely the result of both passive DNA demethylation by DNMTi and active conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) by ten-eleven translocation (TET) enzymes at LTR regions of ERVs, because vitamin C acts as a cofactor for TET proteins. In addition, TET2 knockout reduces the synergy between the two compounds. Furthermore, we show that many patients with hematological neoplasia are markedly vitamin C deficient. Thus, our data suggest that correction of vitamin C deficiency in patients with hematological and other cancers may improve responses to epigenetic therapy with DNMTis.

Comprehensive evaluation of extracellular small RNA isolation methods from serum in high throughput sequencing.

BACKGROUND: DNA and RNA fractions from whole blood, serum and plasma are increasingly popular analytes that are currently under investigation for their utility in the diagnosis and staging of disease. Small non-coding ribonucleic acids (sRNAs), specifically microRNAs (miRNAs) and their variant isoforms (isomiRs), and transfer RNA (tRNA)-derived small RNAs (tDRs) comprise a repertoire of molecules particularly promising in this regard. RESULTS: In this designed study, we compared the performance of various methods and kits for isolating circulating extracellular sRNAs (ex-sRNAs). ex-sRNAs from one healthy individual were isolated using five different isolation kits: Qiagen Circulating Nucleic Acid Kit, ThermoFisher Scientific Ambion TRIzol LS Reagent, Qiagen miRNEasy, QiaSymphony RNA extraction kit and the Exiqon MiRCURY RNA Isolation Kit. Each isolation method was repeated four times. A total of 20 small RNA sequencing (sRNAseq) libraries were constructed, sequenced and compared using a rigorous bioinformatics approach. The Circulating Nucleic Acid Kit had the greatest miRNA isolation variability, but had the lowest isolation variability for other RNA classes (isomiRs, tDRs, and other miscellaneous sRNAs (osRNA). However, the Circulating Nucleic Acid Kit consistently generated the fewest number of reads mapped to the genome, as compared to the best-performing method, Ambion TRIzol, which mapped 10% of the miRNAs, 7.2% of the tDRs and 23.1% of the osRNAs. The other methods performed intermediary, with QiaSymphony mapping 14% of the osRNAs, and miRNEasy mapping 4.6% of the tDRs and 2.9% of the miRNAs, achieving the second best kit performance rating overall. CONCLUSIONS: In summary, each isolation kit displayed different performance characteristics that could be construed as biased or advantageous, depending upon the downstream application and number of samples that require processing.

Hotspot mutations delineating diverse mutational signatures and biological utilities across cancer types.

BACKGROUND: An important step towards personalizing cancer treatment is to integrate heterogeneous evidences to catalog mutational hotspots that are biologically and therapeutically relevant and thus represent where targeted therapy would likely be beneficial. However, existing methods do not sufficiently delineate varying functionality of individual mutations within the same genes. RESULTS: We observed a large discordancy of mutation rates across different mutation subtypes and tumor types, and nominated 702 hotspot mutations in 549 genes in the Catalog of Somatic Mutations in Cancer (COSMIC) by considering context specific mutation characteristics such as genes, cancer types, mutation rates, mutation subtypes and sequence contexts. We observed that hotspot mutations were highly prevalent in Non CpG-island C/G transition and transversion sequence contexts in 10 tumor types, and specific insertion hotspot mutations were enriched in breast cancer and deletion hotspot mutations in colorectal cancer. We found that the hotspot mutations nominated by our approach were significantly more conserved than non-hotspot mutations in the corresponding cancer genes. We also examined the biological significance and pharmacogenomics properties of these hotspot mutations using data in the Cancer Genome Atlas (TCGA) and the Cancer Cell-Line Encyclopedia (CCLE), and found that 53 hotspot mutations are independently associated with diverse functional evidences in 1) mRNA and protein expression, 2) pathway activity, or 3) drug sensitivity and 82 were highly enriched in specific tumor types. We highlighted the distinct functional indications of hotspot mutations under different contexts and nominated novel hotspot mutations such as MAP3K4 A1199 deletion, NR1H2 Q175 insertion, and GATA3 P409 insertion as potential biomarkers or drug targets. CONCLUSION: We identified a set of hotspot mutations across 17 tumor types by considering the background mutation rate variations among genes, tumor subtypes, mutation subtypes, and sequence contexts. We illustrated the common and distinct mutational signatures of hotspot mutations among different tumor types and investigated their variable functional relevance under different contexts, which could potentially serve as a resource for explicitly selecting targets for diagnosis, drug development, and patient management.