Strategies for Integrated Analysis of Genetic, Epigenetic, and Gene Expression Variation in Cancer: Addressing the Challenges

Inter- and transgenerational heritability of preconception chronic stress or alcohol exposure: Translational outcomes in brain and behavior

Chronic stress and alcohol (ethanol) use are highly interrelated and can change an individual's behavior through molecular adaptations that do not change the DNA sequence, but instead change gene expression. A recent wealth of research has found that these nongenomic changes can be transmitted across generations, which could partially account for the "missing heritability" observed in genome-wide association studies of alcohol use disorder and other stress-related neuropsychiatric disorders. In this review, we summarize the molecular and behavioral outcomes of nongenomic inheritance of chronic stress and ethanol exposure and the germline mechanisms that could give rise to this heritability. In doing so, we outline the need for further research to: (1) Investigate individual germline mechanisms of paternal, maternal, and biparental nongenomic chronic stress- and ethanol-related inheritance; (2) Synthesize and dissect cross-generational chronic stress and ethanol exposure; (3) Determine cross-generational molecular outcomes of preconception ethanol exposure that contribute to alcohol-related disease risk, using cancer as an example. A detailed understanding of the cross-generational nongenomic effects of stress and/or ethanol will yield novel insight into the impact of ancestral perturbations on disease risk across generations and uncover actionable targets to improve human health.

Assessment of the Feasibility of a Future Integrated Larger-Scale Epidemiological Study to Evaluate Health Risks of Air Pollution Episodes in Children

Air pollution exposure can lead to exacerbation of respiratory disorders in children. Using sensitive biomarkers helps to assess the impact of air pollution on children’s respiratory health and combining protein, genetic and epigenetic biomarkers gives insights on their interrelatedness. Most studies do not contain such an integrated approach and investigate these biomarkers individually in blood, although its collection in children is challenging. Our study aimed at assessing the feasibility of conducting future integrated larger-scale studies evaluating respiratory health risks of air pollution episodes in children, based on a qualitative analysis of the technical and logistic aspects of a small-scale field study involving 42 children. This included the preparation, collection and storage of non-invasive samples (urine, saliva), the measurement of general and respiratory health parameters and the measurement of specific biomarkers (genetic, protein, epigenetic) of respiratory health and air pollution exposure. Bottlenecks were identified and modifications were proposed to expand this integrated study to a higher number of children, time points and locations. This would allow for non-invasive assessment of the impact of air pollution exposure on the respiratory health of children in future larger-scale studies, which is critical for the development of policies or measures at the population level.

Bayesian tensor factorization-drive breast cancer subtyping by integrating multi-omics data

Breast cancer is a highly heterogeneous disease. Subtyping the disease and identifying the genomic features driving these subtypes are critical for precision oncology for breast cancer. This study focuses on developing a new computational approach for breast cancer subtyping. We proposed to use Bayesian tensor factorization (BTF) to integrate multi-omics data of breast cancer, which include expression profiles of RNA-sequencing, copy number variation, and DNA methylation measured on 762 breast cancer patients from The Cancer Genome Atlas. We applied a consensus clustering approach to identify breast cancer subtypes using the factorized latent features by BTF. Subtype-specific survival patterns of the breast cancer patients were evaluated using Kaplan-Meier (KM) estimators. The proposed approach was compared with other state-of-the-art approaches for cancer subtyping. The BTF-subtyping analysis identified 17 optimized latent components, which were used to reveal six major breast cancer subtypes. Out of all different approaches, only the proposed approach showed distinct survival patterns (p < 0.05). Statistical tests also showed that the identified clusters have statistically significant distributions. Our results showed that the proposed approach is a promising strategy to efficiently use publicly available multi-omics data to identify breast cancer subtypes.

Methylation-eQTL Analysis in Cancer Research

MOTIVATION

DNA methylation is a key epigenetic factor regulating gene expression. While promoter methylation has been well studied, recent publications have revealed that functionally important methylation also occurs in intergenic and distal regions, and varies across genes and tissue types. Given the growing importance of inter-platform integrative genomic analyses, there is an urgent need to develop methods to discover and characterize gene-level relationships between methylation and expression.

RESULTS

We introduce a novel sequential penalized regression approach to identify methylation-expression quantitative trait loci (methyl-eQTLs), a term that we have coined to represent, for each gene and tissue type, a sparse set of CpG loci best explaining gene expression and accompanying weights indicating direction and strength of association. Using TCGA and MD Anderson colorectal cohorts to build and validate our models, we demonstrate our strategy better explains expression variability than current commonly used gene-level methylation summaries. The methyl-eQTLs identified by our approach can be used to construct gene-level methylation summaries that are maximally correlated with gene expression for use in integrative models, and produce a tissue-specific summary of which genes appear to be strongly regulated by methylation. Our results introduce an important resource to the biomedical community for integrative genomics analyses involving DNA methylation.

AVAILABILITY AND IMPLEMENTATION

We produce an R Shiny app (https://rstudio-prd-c1.pmacs.upenn.edu/methyl-eQTL/) that interactively presents methyl-eQTL results for colorectal, breast, and pancreatic cancer. The source R code for this work is provided in the supplement.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of relevant genetic alterations in cancer using topological data analysis

Large-scale cancer genomic studies enable the systematic identification of mutations that lead to the genesis and progression of tumors, uncovering the underlying molecular mechanisms and potential therapies. While some such mutations are recurrently found in many tumors, many others exist solely within a few samples, precluding detection by conventional recurrence-based statistical approaches. Integrated analysis of somatic mutations and RNA expression data across 12 tumor types reveals that mutations of cancer genes are usually accompanied by substantial changes in expression. We use topological data analysis to leverage this observation and uncover 38 elusive candidate cancer-associated genes, including inactivating mutations of the metalloproteinase ADAMTS12 in lung adenocarcinoma. We show that ADAMTS12-/- mice have a five-fold increase in the susceptibility to develop lung tumors, confirming the role of ADAMTS12 as a tumor suppressor gene. Our results demonstrate that data integration through topological techniques can increase our ability to identify previously unreported cancer-related alterations.

Harmonization of quality metrics and power calculation in multi-omic studies

Multi-omic studies combine measurements at different molecular levels to build comprehensive models of cellular systems. The success of a multi-omic data analysis strategy depends largely on the adoption of adequate experimental designs, and on the quality of the measurements provided by the different omic platforms. However, the field lacks a comparative description of performance parameters across omic technologies and a formulation for experimental design in multi-omic data scenarios. Here, we propose a set of harmonized Figures of Merit (FoM) as quality descriptors applicable to different omic data types. Employing this information, we formulate the MultiPower method to estimate and assess the optimal sample size in a multi-omics experiment. MultiPower supports different experimental settings, data types and sample sizes, and includes graphical for experimental design decision-making. MultiPower is complemented with MultiML, an algorithm to estimate sample size for machine learning classification problems based on multi-omic data.

Tightly integrated genomic and epigenomic data mining using tensor decomposition

Motivation

Complex diseases such as cancers often involve multiple types of genomic and/or epigenomic abnormalities. Rapid accumulation of multiple types of omics data demands methods for integrating the multidimensional data in order to elucidate complex relationships among different types of genomic and epigenomic abnormalities.

Results

In the present study, we propose a tightly integrated approach based on tensor decomposition. Multiple types of data, including mRNA, methylation, copy number variations and somatic mutations, are merged into a high-order tensor which is used to develop predictive models for overall survival. The weight tensors of the models are constrained using CANDECOMP/PARAFAC (CP) tensor decomposition and learned using support tensor machine regression (STR) and ridge tensor regression (RTR). The results demonstrate that the tensor decomposition based approaches can achieve better performance than the models based individual data type and the concatenation approach.

Supplementary information

Supplementary data are available at Bioinformatics online.

Dissecting the sources of gene expression variation in a pan-cancer analysis identifies novel regulatory mutations

Although the catalog of cancer-associated mutations in protein-coding regions is nearly complete for all major cancer types, an assessment of regulatory changes in cancer genomes and their clinical significance remain largely preliminary. Adopting bottom-up approach, we quantify the effects of different sources of gene expression variation in a cohort of 3899 samples from 10 cancer types. We find that copy number alterations, epigenetic changes, transcription factors and microRNAs collectively explain, on average, only 31–38% and 18–26% expression variation for cancer-associated and other genes, respectively, and that among these factors copy number alteration has the highest effect. We show that the genes with systematic, large expression variation that could not be attributed to these factors are enriched for pathways related to cancer hallmarks. Integrating whole genome sequencing data and focusing on genes with systematic expression variation we identify novel, recurrent regulatory mutations affecting known cancer genes such as NKX2-1 and GRIN2D in multiple cancer types. Nonetheless, at a genome-wide scale proportions of gene expression variation attributed to recurrent point mutations appear to be modest so far, especially when compared to that attributed to copy number changes – a pattern different from that observed for other complex diseases and traits. We suspect that, owing to plasticity and redundancy in biological pathways, regulatory alterations show complex combinatorial patterns, modulating gene expression in cancer genomes at a finer scale.

Genome-wide functional association networks: background, data & state-of-the-art resources

The vast amount of experimental data from recent advances in the field of high-throughput biology begs for integration into more complex data structures such as genome-wide functional association networks. Such networks have been used for elucidation of the interplay of intra-cellular molecules to make advances ranging from the basic science understanding of evolutionary processes to the more translational field of precision medicine. The allure of the field has resulted in rapid growth of the number of available network resources, each with unique attributes exploitable to answer different biological questions. Unfortunately, the high volume of network resources makes it impossible for the intended user to select an appropriate tool for their particular research question. The aim of this paper is to provide an overview of the underlying data and representative network resources as well as to mention methods of integration, allowing a customized approach to resource selection. Additionally, this report will provide a primer for researchers venturing into the field of network integration.

Bioinformatic analysis of bacteria and host cell dual RNA-sequencing experiments

Bacterial pathogens subvert host cells by manipulating cellular pathways for survival and replication; in turn, host cells respond to the invading pathogen through cascading changes in gene expression. Deciphering these complex temporal and spatial dynamics to identify novel bacterial virulence factors or host response pathways is crucial for improved diagnostics and therapeutics. Dual RNA sequencing (dRNA-Seq) has recently been developed to simultaneously capture host and bacterial transcriptomes from an infected cell. This approach builds on the high sensitivity and resolution of RNA sequencing technology and is applicable to any bacteria that interact with eukaryotic cells, encompassing parasitic, commensal or mutualistic lifestyles. Several laboratory protocols have been presented that outline the collection, extraction and sequencing of total RNA for dRNA-Seq experiments, but there is relatively little guidance available for the detailed bioinformatic analyses required. This protocol outlines a typical dRNA-Seq experiment, based on a Chlamydia trachomatis-infected host cell, with a detailed description of the necessary bioinformatic analyses with currently available software tools.

CAFE MOCHA: An Integrated Platform for Discovering Clinically Relevant Molecular Changes in Cancer-An Example of Distant Metastasis- and Recurrence-Linked Classifiers in Head and Neck Squamous Cell Carcinoma

PURPOSE

With large amounts of multidimensional molecular data on cancers generated and deposited into public repositories such as The Cancer Genome Atlas and International Cancer Genome Consortium, a cancer type agnostic and integrative platform will help to identify signatures with clinical relevance. We devised such a platform and showcase it by identifying a molecular signature for patients with metastatic and recurrent (MR) head and neck squamous cell carcinoma (HNSCC).

METHODS

We devised a statistical framework accompanied by a graphical user interface-driven application, Clinical Association of Functionally Established MOlecular CHAnges ( CAFE MOCHA; https://github.com/binaypanda/CAFEMOCHA), to discover molecular signatures linked to a specific clinical attribute in a cancer type. The platform integrates mutations and indels, gene expression, DNA methylation, and copy number variations to discover a classifier first and then to predict an incoming tumor for the same by pulling defined class variables into a single framework that incorporates a coordinate geometry-based algorithm called complete specificity margin-based clustering, which ensures maximum specificity. CAFE MOCHA classifies an incoming tumor sample using either its matched normal or a built-in database of normal tissues. The application is packed and deployed using the install4j multiplatform installer. We tested CAFE MOCHA in HNSCC tumors (n = 513) followed by validation in tumors from an independent cohort (n = 18) for discovering a signature linked to distant MR.

RESULTS

CAFE MOCHA identified an integrated signature, MR44, associated with distant MR HNSCC, with 80% sensitivity and 100% specificity in the discovery stage and 100% sensitivity and 100% specificity in the validation stage.

CONCLUSION

CAFE MOCHA is a cancer type and clinical attribute agnostic statistical framework to discover integrated molecular signatures.

Stepwise approach to SNP-set analysis illustrated with the Metabochip and colorectal cancer in Japanese Americans of the Multiethnic Cohort

BACKGROUND

Common variants have explained less than the amount of heritability expected for complex diseases, which has led to interest in less-common variants and more powerful approaches to the analysis of whole-genome scans. Because of low frequency (low statistical power), less-common variants are best analyzed using SNP-set methods such as gene-set or pathway-based analyses. However, there is as yet no clear consensus regarding how to focus in on potential risk variants following set-based analyses. We used a stepwise, telescoping approach to analyze common- and rare-variant data from the Illumina Metabochip array to assess genomic association with colorectal cancer (CRC) in the Japanese sub-population of the Multiethnic Cohort (676 cases, 7180 controls). We started with pathway analysis of SNPs that are in genes and pathways having known mechanistic roles in colorectal cancer, then focused on genes within the pathways that evidenced association with CRC, and finally assessed individual SNPs within the genes that evidenced association. Pathway SNPs downloaded from the dbSNP database were cross-matched with Metabochip SNPs and analyzed using the logistic kernel machine regression approach (logistic SNP-set kernel-machine association test, or sequence kernel association test; SKAT) and related methods.

RESULTS

The TGF-β and WNT pathways were associated with all CRC, and the WNT pathway was associated with colon cancer. Individual genes demonstrating the strongest associations were TGFBR2 in the TGF-β pathway and SMAD7 (which is involved in both the TGF-β and WNT pathways). As partial validation of our approach, a known CRC risk variant in SMAD7 (in both the TGF-β and WNT pathways: rs11874392) was associated with CRC risk in our data. We also detected two novel candidate CRC risk variants (rs13075948 and rs17025857) in TGFBR2, a gene known to be associated with CRC risk.

CONCLUSIONS

A stepwise, telescoping approach identified some potentially novel risk variants associated with colorectal cancer, so it may be a useful method for following up on results of set-based SNP analyses. Further work is required to assess the statistical characteristics of the approach, and additional applications should aid in better clarifying its utility.

Multilevel omic data clustering reveals variable contribution of methylator phenotype to integrative cancer subtypes

AIM

We aimed to assess to what extent CpG island methylator phenotype (CIMP) contributes to cancer subtypes obtained by multilevel omic data analysis.

MATERIALS & METHODS

16 The Cancer Genome Atlas datasets encompassing three data layers in 4688 tumor samples were analyzed. We identified cancer integrative subtypes (ISs) by the use of similarity network fusion and consensus clustering. CIMP high (CIMP-H) associated ISs were profiled by gene sets and transcriptional regulators enrichment analysis.

RESULTS & CONCLUSION

In nine out of 16 cancer datasets CIMP-H clusters significantly overlaped with unique ISs. The contribution of CIMP-H on integrative molecular profiling is variable; therefore, only in a subset of cancer types does CIMP-H contribute to homogenous integrative subtype. CIMP-H associated ISs are heterogenous groups with regard to deregulated pathways and transcriptional regulators.

Cross-generational effects of alcohol dependence in humans on HRAS and TP53 methylation in offspring

AIM

We hypothesized that cross-generational effects of alcohol exposure could alter DNA methylation and expression of the HRAS oncogene and TP53 tumor suppressor gene that drive cancer development.

METHODS

DNA methylation of the HRAS and TP53 genes was tested in samples from young participants (Mean age of 13.4 years).

RESULTS

Controlling for both personal use and maternal use of substances during pregnancy, familial alcohol dependence was associated with hypomethylation of CpG sites in the HRAS promoter region and hypermethylation of the TP53 gene.

CONCLUSION

The results suggest that ancestral exposure to alcohol can have enduring effects that impact epigenetic processes such as DNA methylation that controls expression of genes that drive cancer development such as HRAS and TP53.

iGC-an integrated analysis package of gene expression and copy number alteration

BACKGROUND

With the advancement in high-throughput technologies, researchers can simultaneously investigate gene expression and copy number alteration (CNA) data from individual patients at a lower cost. Traditional analysis methods analyze each type of data individually and integrate their results using Venn diagrams. Challenges arise, however, when the results are irreproducible and inconsistent across multiple platforms. To address these issues, one possible approach is to concurrently analyze both gene expression profiling and CNAs in the same individual.

RESULTS

We have developed an open-source R/Bioconductor package (iGC). Multiple input formats are supported and users can define their own criteria for identifying differentially expressed genes driven by CNAs. The analysis of two real microarray datasets demonstrated that the CNA-driven genes identified by the iGC package showed significantly higher Pearson correlation coefficients with their gene expression levels and copy numbers than those genes located in a genomic region with CNA. Compared with the Venn diagram approach, the iGC package showed better performance.

CONCLUSION

The iGC package is effective and useful for identifying CNA-driven genes. By simultaneously considering both comparative genomic and transcriptomic data, it can provide better understanding of biological and medical questions. The iGC package's source code and manual are freely available at https://www.bioconductor.org/packages/release/bioc/html/iGC.html .

The well-accepted notion that gene amplification contributes to increased expression still remains, after all these years, a reasonable but unproven assumption

“Gene amplification causes overexpression” is a longstanding and well-accepted concept in cancer genetics. However, raking the whole literature, we find only statistical analyses showing a positive correlation between gene copy number and expression level, but do not find convincing experimental corroboration for this notion, for most of the amplified oncogenes in cancers. Since an association does not need to be an actual causal relation, in our opinion, this widespread notion still remains a reasonable but unproven assumption awaiting experimental verification.