Heterogeneous DNA methylation contributes to tumorigenesis through inducing the loss of coexpression connectivity in colorectal cancer

Promoter Hypomethylation Upregulates ANXA2 Expression in Pancreatic Cancer and is Associated with Poor Prognosis

Pancreatic cancer (PC) is one of the world's most aggressive and deadly cancers, owing to non-specific early clinical symptoms, late-stage diagnosis, and poor survival. Therefore, it is critical to identify specific biomarkers for its early diagnosis. Annexin A2 (ANXA2) is a calcium-dependent phospholipid-binding protein that has been reported to be upregulated in several cancer types, making it an emerging biomarker and potential cancer therapeutic target. However, the mechanism underlying the regulation of ANXA2 overexpression is still unclear. It is well established that genetic and epigenetic alterations may lead to widespread dysregulation of gene expression. Hence, in this study, we focused on exploring the regulatory mechanism of ANXA2 by investigating the transcriptional profile, methylation pattern, somatic mutation, and prognostic value of ANXA2 in PC using several bioinformatics databases. Our results revealed that the expression levels of ANXA2 were remarkably increased in PC tissues comparing to normal tissues. Furthermore, the high expression of ANXA2 was significantly related to the poor prognosis of PC patients. More importantly, we demonstrated for the first time that the ANXA2 promoter is hypomethylated in PC tissues compared to normal tissues which may result in ANXA2 overexpression in PC. However, more experimental research is required to corroborate our findings.

Genome-Wide Correlation of DNA Methylation and Gene Expression in Postmortem Brain Tissues of Opioid Use Disorder Patients

BACKGROUND

Opioid use disorder (OUD) affects millions of people, causing nearly 50 000 deaths annually in the United States. While opioid exposure and OUD are known to cause widespread transcriptomic and epigenetic changes, few studies in human samples have been conducted. Understanding how OUD affects the brain at the molecular level could help decipher disease pathogenesis and shed light on OUD treatment.

METHODS

We generated genome-wide transcriptomic and DNA methylation profiles of 22 OUD subjects and 19 non-psychiatric controls. We applied weighted gene co-expression network analysis to identify genetic markers consistently associated with OUD at both transcriptomic and methylomic levels. We then performed functional enrichment for biological interpretation. We employed cross-omics analysis to uncover OUD-specific regulatory networks.

RESULTS

We found 6 OUD-associated co-expression gene modules and 6 co-methylation modules (false discovery rate <0.1). Genes in these modules are involved in astrocyte and glial cell differentiation, gliogenesis, response to organic substance, and response to cytokine (false discovery rate <0.05). Cross-omics analysis revealed immune-related transcription regulators, suggesting the role of transcription factor-targeted regulatory networks in OUD pathogenesis.

CONCLUSIONS

Our integrative analysis of multi-omics data in OUD postmortem brain samples suggested complex gene regulatory mechanisms involved in OUD-associated expression patterns. Candidate genes and their upstream regulators revealed in astrocyte, and glial cells could provide new insights into OUD treatment development.

Efficient weighted univariate clustering maps outstanding dysregulated genomic zones in human cancers

Motivation

Chromosomal patterning of gene expression in cancer can arise from aneuploidy, genome disorganization or abnormal DNA methylation. To map such patterns, we introduce a weighted univariate clustering algorithm to guarantee linear runtime, optimality and reproducibility.

Results

We present the chromosome clustering method, establish its optimality and runtime and evaluate its performance. It uses dynamic programming enhanced with an algorithm to reduce search-space in-place to decrease runtime overhead. Using the method, we delineated outstanding genomic zones in 17 human cancer types. We identified strong continuity in dysregulation polarity—dominance by either up- or downregulated genes in a zone—along chromosomes in all cancer types. Significantly polarized dysregulation zones specific to cancer types are found, offering potential diagnostic biomarkers. Unreported previously, a total of 109 loci with conserved dysregulation polarity across cancer types give insights into pan-cancer mechanisms. Efficient chromosomal clustering opens a window to characterize molecular patterns in cancer genome and beyond.

Availability and implementation

Weighted univariate clustering algorithms are implemented within the R package ‘Ckmeans.1d.dp’ (4.0.0 or above), freely available at https://cran.r-project.org/package=Ckmeans.1d.dp.

Supplementary information

Supplementary data are available at Bioinformatics online.

Lower connectivity of tumor coexpression networks is not specific to cancer

Global level network analysis of molecular links is necessary for systems level view of complex diseases like cancer. Using genome-wide expression datasets, we constructed and compared gene co-expression based specific networks of pre-cancerous tumors (adenoma) and cancerous tumors (carcinoma) with paired normal networks to assess for any possible changes in network connectivity. Previously, loss of connectivity was reported as a characteristic of cancer samples. Here, we observed that pre-cancerous conditions also had significantly less connections than paired normal samples. We observed a loss of connectivity trend for colorectal adenoma, aldosterone producing adenoma and uterine leiomyoma. We also showed that the loss of connectivity trend is not specific to positive or negative correlation based networks. Differential hub genes, which were the most highly differentially less connected genes in tumor, were mostly different between different datasets. No common gene list could be defined which underlies the lower connectivity of tumor specific networks. Connectivity of colorectal cancer methylation targets was different from other genes. Extracellular space related terms were enriched in negative correlation based differential hubs and common methylation targets of colorectal carcinoma. Our results indicate a systems level change of lower connectivity as cells transform to not only cancer but also pre-cancerous conditions. This systems level behavior could not be attributed to a group of genes.

Methylation Status of Transcriptional Modulatory Genes Associated with Colorectal Cancer in Northeast China

Background/Aims

Methylation status plays a causal role in carcinogenesis in targeted tissues. However, the relationship between the DNA methylation status of multiple genes in blood leukocytes and colorectal cancer (CRC) susceptibility as well as interactions between dietary factors and CRC risks are unclear.

Methods

We performed a case-control study with 466 CRC patients and 507 cancer-free controls to investigate the association among the methylation status of individual genes, multiple CpG site methylation (MCSM), multiple CpG site heterogeneous methylation and CRC susceptibility. Peripheral blood DNA methylation levels were detected by performing methylation-sensitive high-resolution melting.

Results

Total heterogeneous methylation of CA10 and WT1 conferred a significantly higher risk of CRC (adjusted odds ratio [OR_adjusted], 5.445; 95% confidence interval [CI], 3.075 to 9.643; OR_adjusted, 1.831; 95% CI, 1.100 to 3.047; respectively). Subjects with high-level MCSM (MCSM-H) status demonstrated a higher risk of CRC (OR_adjusted, 4.318; 95% CI, 1.529 to 12.197). Additionally, interactions between the high-level intake of fruit and CRH, WT1, and MCSM on CRC were statistically significant.

Conclusions

The gene methylation status of blood leukocytes may be associated with CRC risk. MCSM-H of blood leukocytes was associated with CRC, especially in younger people. Some dietary factors may affect hypermethylation status and influence susceptibility to CRC.

Genome-wide DNA methylation analysis reveals molecular subtypes of pancreatic cancer

Pancreatic cancer (PC) is the fourth leading cause of cancer deaths in the United States with a five-year patient survival rate of only 6%. Early detection and treatment of this disease is hampered due to lack of reliable diagnostic and prognostic markers. Recent studies have shown that dynamic changes in the global DNA methylation and gene expression patterns play key roles in the PC development; hence, provide valuable insights for better understanding the initiation and progression of PC. In the current study, we used DNA methylation, gene expression, copy number, mutational and clinical data from pancreatic patients. We independently investigated the DNA methylation and differential gene expression profiles between normal and tumor samples and correlated methylation levels with gene expression patterns. We observed a total of ~23-thousand differentially methylated CpG sites (Δβ≥0.1) between normal and tumor samples, where majority of the CpG sites are hypermethylated in PC, and this phenomenon is more prominent in the 5′UTRs and promoter regions compared to the gene bodies. Differential methylation is observed in genes associated with the homeobox domain, cell division and differentiation, cytoskeleton, epigenetic regulation and development, pancreatic development and pancreatic signaling and pancreatic cancer core signaling pathways. Correlation analysis suggests that methylation in the promoter region and 5′UTR has mostly negative correlations with gene expression while gene body and 3′UTR associated methylation has positive correlations. Regulatory element analysis suggests that HOX cluster and histone core proteins are upstream regulators of hypomethylation, while SMAD4, STAT4, STAT5B and zinc finger proteins (ZNF) are upstream regulators of hypermethylation. Non-negative matrix factorization (NMF) clustering of differentially methylated sites generated three clusters in PCs suggesting the existence of distinct molecular subtypes. Cluster 1 and cluster 2 showed samples enriched with clinical phenotypes like neoplasm histological grade and pathologic T-stage T3, respectively, while cluster 3 showed the enrichment of samples with neoplasm histological grade G1. To the best of our knowledge, this is the first genome-scale methylome analysis of PC data from TCGA. Our clustering analysis provides a strong basis for future work on the molecular subtyping of epigenetic regulation in pancreatic cancer.

Individualized network-based drug repositioning infrastructure for precision oncology in the panomics era

Advances in next-generation sequencing technologies have generated the data supporting a large volume of somatic alterations in several national and international cancer genome projects, such as The Cancer Genome Atlas and the International Cancer Genome Consortium. These cancer genomics data have facilitated the revolution of a novel oncology drug discovery paradigm from candidate target or gene studies toward targeting clinically relevant driver mutations or molecular features for precision cancer therapy. This focuses on identifying the most appropriately targeted therapy to an individual patient harboring a particularly genetic profile or molecular feature. However, traditional experimental approaches that are used to develop new chemical entities for targeting the clinically relevant driver mutations are costly and high-risk. Drug repositioning, also known as drug repurposing, re-tasking or re-profiling, has been demonstrated as a promising strategy for drug discovery and development. Recently, computational techniques and methods have been proposed for oncology drug repositioning and identifying pharmacogenomics biomarkers, but overall progress remains to be seen. In this review, we focus on introducing new developments and advances of the individualized network-based drug repositioning approaches by targeting the clinically relevant driver events or molecular features derived from cancer panomics data for the development of precision oncology drug therapies (e.g. one-person trials) to fully realize the promise of precision medicine. We discuss several potential challenges (e.g. tumor heterogeneity and cancer subclones) for precision oncology. Finally, we highlight several new directions for the precision oncology drug discovery via biotherapies (e.g. gene therapy and immunotherapy) that target the 'undruggable' cancer genome in the functional genomics era.

Towards integrated oncogenic marker recognition through mutual information-based statistically significant feature extraction: an association rule mining based study on cancer expression and methylation profiles

Background

Marker detection is an important task in complex disease studies. Here we provide an association rule mining (ARM) based approach for identifying integrated markers through mutual information (MI) based statistically significant feature extraction, and apply it to acute myeloid leukemia (AML) and prostate carcinoma (PC) gene expression and methylation profiles.

Methods

We first collect the genes having both expression and methylation values in AML as well as PC. Next, we run Jarque-Bera normality test on the expression/methylation data to divide the whole dataset into two parts: one that ollows normal distribution and the other that does not follow normal distribution. Thus, we have now four parts of the dataset: normally distributed expression data, normally distributed methylation data, non-normally distributed expression data, and non-normally distributed methylated data. A feature-extraction technique, "mRMR" is then utilized on each part. This results in a list of top-ranked genes. Next, we apply Welch t-test (parametric test) and Shrink t-test (non-parametric test) on the expression/methylation data for the top selected normally distributed genes and non-normally distributed genes, respectively. We then use a recent weighted ARM method, "RANWAR" to combine all/specific resultant genes to generate top oncogenic rules along with respective integrated markers. Finally, we perform literature search as well as KEGG pathway and Gene-Ontology (GO) analyses using Enrichr database for in silico validation of the prioritized oncogenes as the markers and labeling the markers as existing or novel.

Results

The novel markers of AML are {ABCB11↑∪KRT17↓} (i.e., ABCB11 as up-regulated, & KRT17 as down-regulated), and {AP1S1-∪KRT17↓∪NEIL2-∪DYDC1↓}) (i.e., AP1S1 and NEIL2 both as hypo-methylated, & KRT17 and DYDC1 both as down-regulated). The novel marker of PC is {UBIAD1¶∪APBA2‡∪C4orf31‡} (i.e., UBIAD1 as up-regulated and hypo-methylated, & APBA2 and C4orf31 both as down-regulated and hyper-methylated).

Conclusion

The identified novel markers might have critical roles in AML as well as PC. The approach can be applied to other complex disease.

Heterogeneous pattern of DNA methylation in developmentally important genes correlates with its chromatin conformation

Background

DNA methylation is a major epigenetic modification, playing a crucial role in the development and differentiation of higher organisms. DNA methylation is also known to regulate transcription by gene repression. Various developmental genes such as c-mos, HoxB5, Sox11, and Sry show tissue-specific gene expression that was shown to be regulated by promoter DNA methylation. The aim of the present study is to investigate the establishment of chromatin marks (active or repressive) in relation to heterogeneous methylation in the promoter regions of these developmentally important genes.

Results

Chromatin-immunoprecipitation (ChIP) assays were performed to immuno-precipitate chromatin by antibodies against both active (H3K4me3) and repressive (H3K9me3) chromatin regions. The analysis of ChIP results showed that both the percentage input and fold enrichment of activated chromatin was higher in tissues expressing the respective genes as compared to the tissues not expressing the same set of genes. This was true for all the genes selected for the study (c-mos, HoxB5, Sox11, and Sry). These findings illustrate that inconsistent DNA methylation patterns (sporadic, mosaic and heterogeneous) may also influence gene regulation, thereby resulting in the modulation of chromatin conformation.

Conclusions

These findings illustrate that various patterns of DNA methylation (asynchronous, mosaic and heterogeneous) correlates with chromatin modification, resulting in the gene regulation.

Linkage disequilibrium analysis of allelic heterogeneity in DNA methylation

Heterogeneity of DNA methylation status among alleles is observed in various cell types and is involved in epigenetic gene regulation and cancer biology. However, the individual methylation profile within each allele has not yet been examined at the whole-genome level. In the present study, we applied linkage disequilibrium analysis to the DNA methylation data obtained from whole-genome bisulfite sequencing studies in mouse germline and other types of cells. We found that the methylation status of 2 consecutive CpG sites showed deviation from equilibrium frequency toward concordant linkage (both methylated or both unmethylated) in germline cells. In the imprinting loci where methylation of constituent alleles is known, our analysis detected the deviation toward the concordant linkage as expected. In addition, we applied this analysis to the transitional zone between methylated and unmethylated regions and to the cells undergoing epigenetic reprogramming. In both cases, deviation to the concordant-linked alleles was conspicuous, indicating that the methylation pattern is not random but rather concordant within each allele. These results will provide the key to understanding the mechanism underlying allelic heterogeneity.

Transcriptome- and proteome-oriented identification of dysregulated eIF4G, STAT3, and Hippo pathways altered by PIK3CA H1047R in HER2/ER-positive breast cancer

PURPOSE

Phosphatidylinositol 3-kinase (PI3K)/AKT pathway aberrations are common in human breast cancer. Furthermore, PIK3CA mutations are commonly associated with resistance to anti-epidermal growth factor receptor 2 (HER2) or anti-estrogen receptor (ER) agents in HER2 or ER positive (HER2⁺/ER⁺) breast cancer. Hence, deciphering the underlying mechanisms of PIK3CA mutations in HER2⁺/ER⁺ breast cancer would provide novel insights into elucidating resistance to anti-HER2/ER therapies.

METHODS

In this study, we systematically investigated the biological consequences of PIK3CA ^H1047R in HER2⁺/ER⁺ breast cancer by uniquely incorporating mRNA transcriptomic data from The Cancer Genome Atlas and proteomic data from reverse-phase protein arrays.

RESULTS

Our integrative bioinformatics analyses revealed that several important pathways such as STAT3 and VEGF/hypoxia were selectively altered by PIK3CA ^H1047R in HER2⁺/ER⁺ breast cancer. Protein differential expression analysis indicated that an elevated eIF4G might promote tumor angiogenesis and growth via regulation of the hypoxia-activated switch in HER2⁺ PIK3CA ^H1047R breast cancer. We observed hypo-phosphorylation of EGFR in HER2⁺ PIK3CA ^H1047R breast cancer versus HER2⁺PIK3CA^wild-type (PIK3CA ^WT). In addition, ER and PIK3CA ^H1047R might cooperate to activate STAT3, MAPK, AKT, and Hippo pathways in ER⁺ PIK3CA ^H1047R breast cancer. A higher YAP_pS127 level was observed in ER⁺ PIK3CA ^H1047R patients than that in an ER⁺ PIK3CA ^WT subgroup. By examining breast cancer cell lines having both microarray gene expression and drug treatment data from the Genomics of Drug Sensitivity in Cancer and the Stand Up to Cancer datasets, we found that the elevated YAP1 mRNA expression was associated with the resistance of BCL-2 family inhibitors, but with the sensitivity to MEK/MAPK inhibitors in breast cancer cells.

CONCLUSIONS

In summary, these findings shed light on the functional consequences of PIK3CA ^H1047R-driven breast tumorigenesis and resistance to the existing therapeutic agents in HER2⁺/ER⁺ breast cancer.

Integration of genome scale data for identifying new players in colorectal cancer

Colorectal cancers (CRCs) display a wide variety of genomic aberrations that may be either causally linked to their development and progression, or might serve as biomarkers for their presence. Recent advances in rapid high-throughput genetic and genomic analysis have helped to identify a plethora of alterations that can potentially serve as new cancer biomarkers, and thus help to improve CRC diagnosis, prognosis, and treatment. Each distinct data type (copy number variations, gene and microRNAs expression, CpG island methylation) provides an investigator with a different, partially independent, and complementary view of the entire genome. However, elucidation of gene function will require more information than can be provided by analyzing a single type of data. The integration of knowledge obtained from different sources is becoming increasingly essential for obtaining an interdisciplinary view of large amounts of information, and also for cross-validating experimental results. The integration of numerous types of genetic and genomic data derived from public sources, and via the use of ad-hoc bioinformatics tools and statistical methods facilitates the discovery and validation of novel, informative biomarkers. This combinatory approach will also enable researchers to more accurately and comprehensively understand the associations between different biologic pathways, mechanisms, and phenomena, and gain new insights into the etiology of CRC.

Colorectal cancer DNA methylation patterns from patients in Manaus, Brazil

Background

DNA methylation is commonly linked with the silencing of the gene expression for many tumor suppressor genes. As such, determining DNA methylation patterns should aid, in times to come, in the diagnosis and personal treatment for various types of cancers. Here, we analyzed the methylation pattern from five colorectal cancer patients from the Amazon state in Brazil for four tumor suppressor genes, viz.: DAPK, CDH1, CDKN2A, and TIMP2 by employing a polymerase chain reaction (PCR) specific to methylation. Efforts in the study of colorectal cancer are fundamental as it is the third most of highest incidence in the world.

Results

Tumor biopsies were methylated in 1/5 (20 %), 2/5 (40 %), 4/5 (80 %), and 4/5 (80 %) for CDH1, CDKN2A, DAPK, and TIMP2 genes, respectively. The margin biopsies were methylated in 3/7 (43 %), 2/7 (28 %), 7/7 (100 %), and 6/7 (86 %) for CDH1, CDKN2A, DAPK, and TIMP2, respectively.

Conclusions

Our findings showed DAPK and TIMP2 to be methylated in most samples from both tumor tissues and adjacent non-neoplastic margins; thus presenting distinct methylation patterns. This emphasizes the importance of better understanding of the relation of these patterns with cancer in the context of different populations.

DREAMing: a simple and ultrasensitive method for assessing intratumor epigenetic heterogeneity directly from liquid biopsies

Many cancers comprise heterogeneous populations of cells at primary and metastatic sites throughout the body. The presence or emergence of distinct subclones with drug-resistant genetic and epigenetic phenotypes within these populations can greatly complicate therapeutic intervention. Liquid biopsies of peripheral blood from cancer patients have been suggested as an ideal means of sampling intratumor genetic and epigenetic heterogeneity for diagnostics, monitoring and therapeutic guidance. However, current molecular diagnostic and sequencing methods are not well suited to the routine assessment of epigenetic heterogeneity in difficult samples such as liquid biopsies that contain intrinsically low fractional concentrations of circulating tumor DNA (ctDNA) and rare epigenetic subclonal populations. Here we report an alternative approach, deemed DREAMing (Discrimination of Rare EpiAlleles by Melt), which uses semi-limiting dilution and precise melt curve analysis to distinguish and enumerate individual copies of epiallelic species at single-CpG-site resolution in fractions as low as 0.005%, providing facile and inexpensive ultrasensitive assessment of locus-specific epigenetic heterogeneity directly from liquid biopsies. The technique is demonstrated here for the evaluation of epigenetic heterogeneity at p14^ARF and BRCA1 gene-promoter loci in liquid biopsies obtained from patients in association with non-small cell lung cancer (NSCLC) and myelodysplastic/myeloproliferative neoplasms (MDS/MPN), respectively.