Widespread genetic epistasis among cancer genes

Detecting genetic effects on phenotype variability to capture gene-by-environment interactions: a systematic method comparison

Genetically associated phenotypic variability has been widely observed across organisms and traits, including in humans. Both gene-gene and gene-environment interactions can lead to an increase in genetically associated phenotypic variability. Therefore, detecting the underlying genetic variants, or variance Quantitative Trait Loci (vQTLs), can provide novel insights into complex traits. Established approaches to detect vQTLs apply different methodologies from variance-only approaches to mean-variance joint tests, but a comprehensive comparison of these methods is lacking. Here, we review available methods to detect vQTLs in humans, carry out a simulation study to assess their performance under different biological scenarios of gene-environment interactions, and apply the optimal approaches for vQTL identification to gene expression data. Overall, with a minor allele frequency (MAF) of less than 0.2, the squared residual value linear model (SVLM) and the deviation regression model (DRM) are optimal when the data follow normal and non-normal distributions, respectively. In addition, the Brown-Forsythe (BF) test is one of the optimal methods when the MAF is 0.2 or larger, irrespective of phenotype distribution. Additionally, a larger sample size and more balanced sample distribution in different exposure categories increase the power of BF, SVLM, and DRM. Our results highlight vQTL detection methods that perform optimally under realistic simulation settings and show that their relative performance depends on the phenotype distribution, allele frequency, sample size, and the type of exposure in the interaction model underlying the vQTL.

Synergistic epistasis among cancer drivers can rescue early tumors from the accumulation of deleterious passengers

Epistasis among driver mutations is pervasive and explains relevant features of cancer, such as differential therapy response and convergence towards well-characterized molecular subtypes. Furthermore, a growing body of evidence suggests that tumor development could be hampered by the accumulation of slightly deleterious passenger mutations. In this work, we combined empirical epistasis networks, computer simulations, and mathematical models to explore how synergistic interactions among driver mutations affect cancer progression under the burden of slightly deleterious passengers. We found that epistasis plays a crucial role in tumor development by promoting the transformation of precancerous clones into rapidly growing tumors through a process that is analogous to evolutionary rescue. The triggering of epistasis-driven rescue is strongly dependent on the intensity of epistasis and could be a key rate-limiting step in many tumors, contributing to their unpredictability. As a result, central genes in cancer epistasis networks appear as key intervention targets for cancer therapy.

SLKB: synthetic lethality knowledge base

Emerging CRISPR-Cas9 technology permits synthetic lethality (SL) screening of large number of gene pairs from gene combination double knockout (CDKO) experiments. However, the poor integration and annotation of CDKO SL data in current SL databases limit their utility, and diverse methods of calculating SL scores prohibit their comparison. To overcome these shortcomings, we have developed SL knowledge base (SLKB) that incorporates data of 11 CDKO experiments in 22 cell lines, 16,059 SL gene pairs and 264,424 non-SL gene pairs. Additionally, within SLKB, we have implemented five SL calculation methods: median score with and without background control normalization (Median-B/NB), sgRNA-derived score (sgRNA-B/NB), Horlbeck score, GEMINI score and MAGeCK score. The five scores have demonstrated a mere 1.21% overlap among their top 10% SL gene pairs, reflecting high diversity. Users can browse SL networks and assess the impact of scoring methods using Venn diagrams. The SL network generated from all data in SLKB shows a greater likelihood of SL gene pair connectivity with other SL gene pairs than non-SL pairs. Comparison of SL networks between two cell lines demonstrated greater likelihood to share SL hub genes than SL gene pairs. SLKB website and pipeline can be freely accessed at https://slkb.osubmi.org and https://slkb.docs.osubmi.org/, respectively.

ELISL: early-late integrated synthetic lethality prediction in cancer

MOTIVATION

Anti-cancer therapies based on synthetic lethality (SL) exploit tumour vulnerabilities for treatment with reduced side effects, by targeting a gene that is jointly essential with another whose function is lost. Computational prediction is key to expedite SL screening, yet existing methods are vulnerable to prevalent selection bias in SL data and reliant on cancer or tissue type-specific omics, which can be scarce. Notably, sequence similarity remains underexplored as a proxy for related gene function and joint essentiality.

RESULTS

We propose ELISL, Early-Late Integrated SL prediction with forest ensembles, using context-free protein sequence embeddings and context-specific omics from cell lines and tissue. Across eight cancer types, ELISL showed superior robustness to selection bias and recovery of known SL genes, as well as promising cross-cancer predictions. Co-occurring mutations in a BRCA gene and ELISL-predicted pairs from the HH, FGF, WNT, or NEIL gene families were associated with longer patient survival times, revealing therapeutic potential.

AVAILABILITY AND IMPLEMENTATION

Data: 10.6084/m9.figshare.23607558 & Code: github.com/joanagoncalveslab/ELISL.

Pairwise and higher-order epistatic effects among somatic cancer mutations across oncogenesis

Cancer occurs as a consequence of multiple somatic mutations that lead to uncontrolled cell growth. Mutual exclusivity and co-occurrence of mutations imply-but do not prove-that mutations exert synergistic or antagonistic epistatic effects on oncogenesis. Knowledge of these interactions, and the consequent trajectories of mutation and selection that lead to cancer has been a longstanding goal within the cancer research community. Recent research has revealed mutation rates and scaled selection coefficients for specific recurrent variants across many cancer types. However, there are no current methods to quantify the strength of selection incorporating pairwise and higher-order epistatic effects on selection within the trajectory of likely cancer genotoypes. Therefore, we have developed a continuous-time Markov chain model that enables the estimation of mutation origination and fixation (flux), dependent on somatic cancer genotype. Coupling this continuous-time Markov chain model with a deconvolution approach provides estimates of underlying mutation rates and selection across the trajectory of oncogenesis. We demonstrate computation of fluxes and selection coefficients in a somatic evolutionary model for the four most frequently variant driver genes (TP53, LRP1B, KRAS and STK11) from 565 cases of lung adenocarcinoma. Our analysis reveals multiple antagonistic epistatic effects that reduce the possible routes of oncogenesis, and inform cancer research regarding viable trajectories of somatic evolution whose progression could be forestalled by precision medicine. Synergistic epistatic effects are also identified, most notably in the somatic genotype TP53 LRP1B for mutations in the KRAS gene, and in somatic genotypes containing KRAS or TP53 mutations for mutations in the STK11 gene. Large positive fluxes of KRAS variants were driven by large selection coefficients, whereas the flux toward LRP1B mutations was substantially aided by a large mutation rate for this gene. The approach enables inference of the most likely routes of site-specific variant evolution and estimation of the strength of selection operating on each step along the route, a key component of what we need to know to develop and implement personalized cancer therapies.

High-Throughput Identification, Modeling, and Analysis of Cancer Driver Genes In Vivo

The vast number of genomic and molecular alterations in cancer pose a substantial challenge to uncovering the mechanisms of tumorigenesis and identifying therapeutic targets. High-throughput functional genomic methods in genetically engineered mouse models allow for rapid and systematic investigation of cancer driver genes. In this review, we discuss the basic concepts and tools for multiplexed investigation of functionally important cancer genes in vivo using autochthonous cancer models. Furthermore, we highlight emerging technical advances in the field, potential opportunities for future investigation, and outline a vision for integrating multiplexed genetic perturbations with detailed molecular analyses to advance our understanding of the genetic and molecular basis of cancer.

Synthetic lethal gene pairs: Experimental approaches and predictive models

Synthetic lethality (SL) refers to a genetic interaction in which the simultaneous perturbation of two genes leads to cell or organism death, whereas viability is maintained when only one of the pair is altered. The experimental exploration of these pairs and predictive modeling in computational biology contribute to our understanding of cancer biology and the development of cancer therapies. We extensively reviewed experimental technologies, public data sources, and predictive models in the study of synthetic lethal gene pairs and herein detail biological assumptions, experimental data, statistical models, and computational schemes of various predictive models, speculate regarding their influence on individual sample- and population-based synthetic lethal interactions, discuss the pros and cons of existing SL data and models, and highlight potential research directions in SL discovery.

Inferring Epistasis from Genetic Time-series Data

Epistasis refers to fitness or functional effects of mutations that depend on the sequence background in which these mutations arise. Epistasis is prevalent in nature, including populations of viruses, bacteria, and cancers, and can contribute to the evolution of drug resistance and immune escape. However, it is difficult to directly estimate epistatic effects from sampled observations of a population. At present, there are very few methods that can disentangle the effects of selection (including epistasis), mutation, recombination, genetic drift, and genetic linkage in evolving populations. Here we develop a method to infer epistasis, along with the fitness effects of individual mutations, from observed evolutionary histories. Simulations show that we can accurately infer pairwise epistatic interactions provided that there is sufficient genetic diversity in the data. Our method also allows us to identify which fitness parameters can be reliably inferred from a particular data set and which ones are unidentifiable. Our approach therefore allows for the inference of more complex models of selection from time-series genetic data, while also quantifying uncertainty in the inferred parameters.

Pervasive conditional selection of driver mutations and modular epistasis networks in cancer

Cancer driver mutations often display mutual exclusion or co-occurrence, underscoring the key role of epistasis in carcinogenesis. However, estimating the magnitude of epistasis and quantifying its effect on tumor evolution remains a challenge. We develop a method (Coselens) to quantify conditional selection on the excess of nonsynonymous substitutions in cancer genes. Coselens infers the number of drivers per gene in different partitions of a cancer genomics dataset using covariance-based mutation models and determines whether coding mutations in a gene affect selection for drivers in any other gene. Using Coselens, we identify 296 conditionally selected gene pairs across 16 cancer types in the TCGA dataset. Conditional selection affects 25%-50% of driver substitutions in tumors with >2 drivers. Conditionally co-selected genes form modular networks, whose structures challenge the traditional interpretation of within-pathway mutual exclusivity and across-pathway synergy, suggesting a more complex scenario where gene-specific across-pathway epistasis shapes differentiated cancer subtypes.

EPIMUTESTR: a nearest neighbor machine learning approach to predict cancer driver genes from the evolutionary action of coding variants

Discovering rare cancer driver genes is difficult because their mutational frequency is too low for statistical detection by computational methods. EPIMUTESTR is an integrative nearest-neighbor machine learning algorithm that identifies such marginal genes by modeling the fitness of their mutations with the phylogenetic Evolutionary Action (EA) score. Over cohorts of sequenced patients from The Cancer Genome Atlas representing 33 tumor types, EPIMUTESTR detected 214 previously inferred cancer driver genes and 137 new candidates never identified computationally before of which seven genes are supported in the COSMIC Cancer Gene Census. EPIMUTESTR achieved better robustness and specificity than existing methods in a number of benchmark methods and datasets.

Statistical tests for intra-tumour clonal co-occurrence and exclusivity

Tumour progression is an evolutionary process in which different clones evolve over time, leading to intra-tumour heterogeneity. Interactions between clones can affect tumour evolution and hence disease progression and treatment outcome. Intra-tumoural pairs of mutations that are overrepresented in a co-occurring or clonally exclusive fashion over a cohort of patient samples may be suggestive of a synergistic effect between the different clones carrying these mutations. We therefore developed a novel statistical testing framework, called GeneAccord, to identify such gene pairs that are altered in distinct subclones of the same tumour. We analysed our framework for calibration and power. By comparing its performance to baseline methods, we demonstrate that to control type I errors, it is essential to account for the evolutionary dependencies among clones. In applying GeneAccord to the single-cell sequencing of a cohort of 123 acute myeloid leukaemia patients, we find 1 clonally co-occurring and 8 clonally exclusive gene pairs. The clonally exclusive pairs mostly involve genes of the key signalling pathways.

Tumours typically display high levels of heterogeneity, not only between different tumours but also within a single one. Intra-tumour heterogeneity results from an evolutionary process, giving rise to different populations of cancer cells known as clones. How clones interact may affect tumour evolution, which in turn determines disease progression and treatment outcome. In practice, we may observe pairs of mutations that co-occur in clones or exclude each other more often than we would expect for a given cohort of patient samples. Exclusive pairs are suggestive that clones carrying one or the other mutation may cooperate in the evolutionary process. Targeting only one of them may then suffice to alter the tumour evolution. Therefore it is critical to have statistical methods which allow us to identify such pairs. GeneAccord is a novel statistical testing framework we developed especially to identify pairs of genes altered in distinct clones of the same tumour. Accounting for the evolutionary dependencies among clones emerged as critical to adequately control testing errors. In a cohort of 123 acute myeloid leukaemia patients, GeneAccord identified one clonally co-occurring and eight clonally exclusive gene pairs. The latter predominantly involved genes of key signalling pathways.

Model guided trait-specific co-expression network estimation as a new perspective for identifying molecular interactions and pathways

A wide variety of 1) parametric regression models and 2) co-expression networks have been developed for finding gene-by-gene interactions underlying complex traits from expression data. While both methodological schemes have their own well-known benefits, little is known about their synergistic potential. Our study introduces their methodological fusion that cross-exploits the strengths of individual approaches via a built-in information-sharing mechanism. This fusion is theoretically based on certain trait-conditioned dependency patterns between two genes depending on their role in the underlying parametric model. Resulting trait-specific co-expression network estimation method 1) serves to enhance the interpretation of biological networks in a parametric sense, and 2) exploits the underlying parametric model itself in the estimation process. To also account for the substantial amount of intrinsic noise and collinearities, often entailed by expression data, a tailored co-expression measure is introduced along with this framework to alleviate related computational problems. A remarkable advance over the reference methods in simulated scenarios substantiate the method's high-efficiency. As proof-of-concept, this synergistic approach is successfully applied in survival analysis, with acute myeloid leukemia data, further highlighting the framework's versatility and broad practical relevance.

Genome-wide Screens Implicate Loss of Cullin Ring Ligase 3 in Persistent Proliferation and Genome Instability in TP53-Deficient Cells

TP53 deficiency is the most common alteration in cancer; however, this alone is typically insufficient to drive tumorigenesis. To identify genes promoting tumorigenesis in combination with TP53 deficiency, we perform genome-wide CRISPR-Cas9 knockout screens coupled with proliferation and transformation assays in isogenic cell lines. Loss of several known tumor suppressors enhances cellular proliferation and transformation. Loss of neddylation pathway genes promotes uncontrolled proliferation exclusively in TP53-deficient cells. Combined loss of CUL3 and TP53 activates an oncogenic transcriptional program governed by the nuclear factor κB (NF-κB), AP-1, and transforming growth factor β (TGF-β) pathways. This program maintains persistent cellular proliferation, induces partial epithelial to mesenchymal transition, and increases DNA damage, genomic instability, and chromosomal rearrangements. Our findings reveal CUL3 loss as a key event stimulating persistent proliferation in TP53-deficient cells. These findings may be clinically relevant, since TP53-CUL3-deficient cells are highly sensitive to ataxia telangiectasia mutated (ATM) inhibition, exposing a vulnerability that could be exploited for cancer treatment.

•

Mixed-effect models with MEMcrispR applied to CRISPR screen analyses

•

Knockout of neddylation genes increases persistent proliferation in TP53^−/− cells

•

TP53^−/−,CUL3^−/− cells exhibit persistent proliferation and partial EMT phenotype

•

TP53^−/−,CUL3^−/− cells show increased DNA damage and display sensitivity to ATM inhibition

Genome-wide epistasis analysis for Alzheimer's disease and implications for genetic risk prediction

BACKGROUND

Single-nucleotide polymorphisms (SNPs) identified by genome-wide association studies only explain part of the heritability of Alzheimer's disease (AD). Epistasis has been considered as one of the main causes of "missing heritability" in AD.

METHODS

We performed genome-wide epistasis screening (N = 10,389) for the clinical diagnosis of AD using three popularly adopted methods. Subsequent analyses were performed to eliminate spurious associations caused by possible confounding factors. Then, candidate genetic interactions were examined for their co-expression in the brains of AD patients and analyzed for their association with intermediate AD phenotypes. Moreover, a new approach was developed to compile the epistasis risk factors into an epistasis risk score (ERS) based on multifactor dimensional reduction. Two independent datasets were used to evaluate the feasibility of ERSs in AD risk prediction.

RESULTS

We identified 2 candidate genetic interactions with PFDR <  0.05 (RAMP3-SEMA3A and NSMCE1-DGKE/C17orf67) and another 5 genetic interactions with PFDR <  0.1. Co-expression between the identified interactions supported the existence of possible biological interactions underlying the observed statistical significance. Further association of candidate interactions with intermediate phenotypes helps explain the mechanisms of neuropathological alterations involved in AD. Importantly, we found that ERSs can identify high-risk individuals showing earlier onset of AD. Combined risk scores of SNPs and SNP-SNP interactions showed slightly but steadily increased AUC in predicting the clinical status of AD.

CONCLUSIONS

In summary, we performed a genome-wide epistasis analysis to identify novel genetic interactions potentially implicated in AD. We found that ERS can serve as an indicator of the genetic risk of AD.

GWAS of three molecular traits highlights core genes and pathways alongside a highly polygenic background

Genome-wide association studies (GWAS) have been used to study the genetic basis of a wide variety of complex diseases and other traits. We describe UK Biobank GWAS results for three molecular traits-urate, IGF-1, and testosterone-with better-understood biology than most other complex traits. We find that many of the most significant hits are readily interpretable. We observe huge enrichment of associations near genes involved in the relevant biosynthesis, transport, or signaling pathways. We show how GWAS data illuminate the biology of each trait, including differences in testosterone regulation between females and males. At the same time, even these molecular traits are highly polygenic, with many thousands of variants spread across the genome contributing to trait variance. In summary, for these three molecular traits we identify strong enrichment of signal in putative core gene sets, even while most of the SNP-based heritability is driven by a massively polygenic background.

Positive epistasis between disease-causing missense mutations and silent polymorphism with effect on mRNA translation velocity

Epistasis refers to the dependence of a mutation on other mutation(s) and the genetic context in general. In the context of human disorders, epistasis complicates the spectrum of disease symptoms and has been proposed as a major contributor to variations in disease outcome. The nonadditive relationship between mutations and the lack of complete understanding of the underlying physiological effects limit our ability to predict phenotypic outcome. Here, we report positive epistasis between intragenic mutations in the cystic fibrosis transmembrane conductance regulator (CFTR)-the gene responsible for cystic fibrosis (CF) pathology. We identified a synonymous single-nucleotide polymorphism (sSNP) that is invariant for the CFTR amino acid sequence but inverts translation speed at the affected codon. This sSNP in cis exhibits positive epistatic effects on some CF disease-causing missense mutations. Individually, both mutations alter CFTR structure and function, yet when combined, they lead to enhanced protein expression and activity. The most robust effect was observed when the sSNP was present in combination with missense mutations that, along with the primary amino acid change, also alter the speed of translation at the affected codon. Functional studies revealed that synergistic alteration in ribosomal velocity is the underlying mechanism; alteration of translation speed likely increases the time window for establishing crucial domain-domain interactions that are otherwise perturbed by each individual mutation.

Leveraging phenotypic variability to identify genetic interactions in human phenotypes

Although thousands of loci have been associated with human phenotypes, the role of gene-environment (GxE) interactions in determining individual risk of human diseases remains unclear. This is partly because of the severe erosion of statistical power resulting from the massive number of statistical tests required to detect such interactions. Here, we focus on improving the power of GxE tests by developing a statistical framework for assessing quantitative trait loci (QTLs) associated with the trait means and/or trait variances. When applying this framework to body mass index (BMI), we find that GxE discovery and replication rates are significantly higher when prioritizing genetic variants associated with the variance of the phenotype (vQTLs) compared to when assessing all genetic variants. Moreover, we find that vQTLs are enriched for associations with other non-BMI phenotypes having strong environmental influences, such as diabetes or ulcerative colitis. We show that GxE effects first identified in quantitative traits such as BMI can be used for GxE discovery in disease phenotypes such as diabetes. A clear conclusion is that strong GxE interactions mediate the genetic contribution to body weight and diabetes risk.

Deep phylogeny of cancer drivers and compensatory mutations

Driver mutations (DM) are the genetic impetus for most cancers. The DM are assumed to be deleterious in species evolution, being eliminated by purifying selection unless compensated by other mutations. We present deep phylogenies for 84 cancer driver genes and investigate the prevalence of 434 DM across gene-species trees. The DM are rare in species evolution, and 181 are completely absent, validating their negative fitness effect. The DM are more common in unicellular than in multicellular eukaryotes, suggesting a link between these mutations and cell proliferation control. 18 DM appear as the ancestral state in one or more major clades, including 3 among mammals. We identify within-gene, compensatory mutations for 98 DM and infer likely interactions between the DM and compensatory sites in protein structures. These findings elucidate the evolutionary status of DM and are expected to advance the understanding of the functions and evolution of oncogenes and tumor suppressors.

Rochman et al. present deep phylogenies for 84 cancer driver genes and examine the prevalence of driver mutations across gene-species trees. Their results show that driver mutations are rare in species evolution and give insight into the evolution of driver mutations and oncogenes.

Genome-wide interaction analysis of pathological hallmarks in Alzheimer's disease

Genome-wide association studies have identified many loci associated with Alzheimer's dementia. However, these variants only explain part of the heritability of Alzheimer's disease (AD). As genetic epistasis can be a major contributor to the "missing heritability" of AD, we conducted genome-wide epistasis screening for AD pathologies in 2 independent cohorts. First, we performed a genome-wide epistasis study of AD-related brain pathologies (N_max = 1318) in ROS/MAP. Candidate interactions were validated using cerebrospinal fluid biomarkers of AD in ADNI (N_max = 1128). Further functional analysis tested the association of candidate interactions with neuroimaging phenotypes. For tau and amyloid-β pathology, we identified 2803 and 464 candidate SNP-SNP interactions, respectively. Associations of candidate SNP-SNP interactions with brain volume and white matter changes from neuroimages provides additional insights into their molecular functions. Transcriptional analysis supported possible gene-gene interactions identified by statistical screening through their co-expression in the brain. In summary, we outlined an exhaustive epistasis analysis to identify novel genetic interactions with potential roles in AD pathologies. We further delved into the functional relevance of candidate interactions by association with neuroimaging phenotypes and analysis of co-expression between corresponding gene pairs.

Clonal evolution driven by superdriver mutations

Background

Tumors are widely recognized to progress through clonal evolution by sequentially acquiring selectively advantageous genetic alterations that significantly contribute to tumorigenesis and thus are termned drivers. Some cancer drivers, such as TP53 point mutation or EGFR copy number gain, provide exceptional fitness gains, which, in time, can be sufficient to trigger the onset of cancer with little or no contribution from additional genetic alterations. These key alterations are called superdrivers.

Results

In this study, we employ a Wright-Fisher model to study the interplay between drivers and superdrivers in tumor progression. We demonstrate that the resulting evolutionary dynamics follow global clonal expansions of superdrivers with periodic clonal expansions of drivers. We find that the waiting time to the accumulation of a set of superdrivers and drivers in the tumor cell population can be approximated by the sum of the individual waiting times.

Conclusions

Our results suggest that superdriver dynamics dominate over driver dynamics in tumorigenesis. Furthermore, our model allows studying the interplay between superdriver and driver mutations both empirically and theoretically.

Link synthetic lethality to drug sensitivity of cancer cells

Synthetic lethal (SL) interactions occur when alterations in two genes lead to cell death but alteration in only one of them is not lethal. SL interactions provide a new strategy for molecular-targeted cancer therapy. Currently, there are few drugs targeting SL interactions that entered into clinical trials. Therefore, it is necessary to investigate the link between SL interactions and drug sensitivity of cancer cells systematically for drug development purpose. We identified SL interactions by integrating the high-throughput data from The Cancer Genome Atlas, small hairpin RNA data and genetic interactions of yeast. By integrating SL interactions from other studies, we tested whether the SL pairs that consist of drug target genes and the genes with genomic alterations are related with drug sensitivity of cancer cells. We found that only 6.26%∼34.61% of SL interactions showed the expected significant drug sensitivity using the pooled cancer cell line data from different tissues, but the proportion increased significantly to approximately 90% using the cancer cell line data for each specific tissue. From an independent pharmacogenomics data of 41 breast cancer cell lines, we found three SL interactions (ABL1-IFI16, ABL1-SLC50A1 and ABL1-SYT11) showed significantly better prognosis for the patients with both genes being altered than the patients with only one gene being altered, which partially supports the SL effect between the gene pairs. Our study not only provides a new way for unraveling the complex mechanisms of drug sensitivity but also suggests numerous potentially important drug targets for cancer therapy.

Causal Reasoning on Boolean Control Networks Based on Abduction: Theory and Application to Cancer Drug Discovery

Complex diseases such as Cancer or Alzheimer's are caused by multiple molecular perturbations leading to pathological cellular behavior. However, the identification of disease-induced molecular perturbations and subsequent development of efficient therapies are challenged by the complexity of the genotype-phenotype relationship. Accordingly, a key issue is to develop frameworks relating molecular perturbations and drug effects to their consequences on cellular phenotypes. Such framework would aim at identifying the sets of causal molecular factors leading to phenotypic reprogramming. In this article, we propose a theoretical framework, called Boolean Control Networks, where disease-induced molecular perturbations and drug actions are seen as topological perturbations/actions on molecular networks leading to cell phenotype reprogramming. We present a new method using abductive reasoning principles inferring the minimal causal topological actions leading to an expected behavior at stable state. Then, we compare different implementations of the algorithm and finally, show a proof-of-concept of the approach on a model of network regulating the proliferation/apoptosis switch in breast cancer by automatically discovering driver genes and their synthetic lethal drug target partner.

A Gaussian process-based definition reveals new and bona fide genetic interactions compared to a multiplicative model in the Gram-negative Escherichia coli

MOTIVATION

A digenic genetic interaction (GI) is observed when mutations in two genes within the same organism yield a phenotype that is different from the expected, given each mutation's individual effects. While multiplicative scoring is widely applied to define GIs, revealing underlying gene functions, it remains unclear if it is the most suitable choice for scoring GIs in Escherichia coli. Here, we assess many different definitions, including the multiplicative model, for mapping functional links between genes and pathways in E.coli.

RESULTS

Using our published E.coli GI datasets, we show computationally that a machine learning Gaussian process (GP)-based definition better identifies functional associations among genes than a multiplicative model, which we have experimentally confirmed on a set of gene pairs. Overall, the GP definition improves the detection of GIs, biological reasoning of epistatic connectivity, as well as the quality of GI maps in E.coli, and, potentially, other microbes.

AVAILABILITY AND IMPLEMENTATION

The source code and parameters used to generate the machine learning models in WEKA software were provided in the Supplementary information.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Animal models of well-differentiated/dedifferentiated liposarcoma: utility and limitations

Liposarcoma is a malignant neoplasm of fat tissue. Well-differentiated and dedifferentiated liposarcoma (WDL/DDL) represent the two most clinically observed histotypes occurring in middle-aged to older adults, particularly within the retroperitoneum or extremities. WDL/DDL are thought to represent the broad spectrum of one disease, as they are both associated with the amplification in the chromosomal 12q13-15 region that causes MDM2 and CDK4 overexpression, the most useful predictor for liposarcoma diagnosis. In comparison to WDL, DDL contains additional genetic abnormalities, principally coamplifications of 1p32 and 6q23, that increase recurrence and metastatic rate. In this review, we discuss the xenograft and transgenic animal models generated for studying progression of WDL/DDL, highlighting utilities and pitfalls in such approaches that can facilitate or impede the development of new therapies.

Genetic interaction analysis among oncogenesis-related genes revealed novel genes and networks in lung cancer development

The development of cancer is driven by the accumulation of many oncogenesis-related genetic alterations and tumorigenesis is triggered by complex networks of involved genes rather than independent actions. To explore the epistasis existing among oncogenesis-related genes in lung cancer development, we conducted pairwise genetic interaction analyses among 35,031 SNPs from 2027 oncogenesis-related genes. The genotypes from three independent genome-wide association studies including a total of 24,037 lung cancer patients and 20,401 healthy controls with Caucasian ancestry were analyzed in the study. Using a two-stage study design including discovery and replication studies, and stringent Bonferroni correction for multiple statistical analysis, we identified significant genetic interactions between SNPs in RGL1:RAD51B (OR=0.44, p value=3.27x10-11 in overall lung cancer and OR=0.41, p value=9.71x10-11 in non-small cell lung cancer), SYNE1:RNF43 (OR=0.73, p value=1.01x10-12 in adenocarcinoma) and FHIT:TSPAN8 (OR=1.82, p value=7.62x10-11 in squamous cell carcinoma) in our analysis. None of these genes have been identified from previous main effect association studies in lung cancer. Further eQTL gene expression analysis in lung tissues provided information supporting the functional role of the identified epistasis in lung tumorigenesis. Gene set enrichment analysis revealed potential pathways and gene networks underlying molecular mechanisms in overall lung cancer as well as histology subtypes development. Our results provide evidence that genetic interactions between oncogenesis-related genes play an important role in lung tumorigenesis and epistasis analysis, combined with functional annotation, provides a valuable tool for uncovering functional novel susceptibility genes that contribute to lung cancer development by interacting with other modifier genes.

The structure and evolution of eukaryotic chaperonin-containing TCP-1 and its mechanism that folds actin into a protein spring

Actin is folded to its native state in eukaryotic cytosol by the sequential allosteric mechanism of the chaperonin-containing TCP-1 (CCT). The CCT machine is a double-ring ATPase built from eight related subunits, CCT1–CCT8. Non-native actin interacts with specific subunits and is annealed slowly through sequential binding and hydrolysis of ATP around and across the ring system. CCT releases a folded but soft ATP-G-actin monomer which is trapped 80 kJ/mol uphill on the folding energy surface by its ATP-Mg2+/Ca2+ clasp. The energy landscape can be re-explored in the actin filament, F-actin, because ATP hydrolysis produces dehydrated and more compact ADP-actin monomers which, upon application of force and strain, are opened and closed like the elements of a spring. Actin-based myosin motor systems underpin a multitude of force generation processes in cells and muscles. We propose that the water surface of F-actin acts as a low-binding energy, directional waveguide which is recognized specifically by the myosin lever-arm domain before the system engages to form the tight-binding actomyosin complex. Such a water-mediated recognition process between actin and myosin would enable symmetry breaking through fast, low energy initial binding events. The origin of chaperonins and the subsequent emergence of the CCT–actin system in LECA (last eukaryotic common ancestor) point to the critical role of CCT in facilitating phagocytosis during early eukaryotic evolution and the transition from the bacterial world. The coupling of CCT-folding fluxes to the cell cycle, cell size control networks and cancer are discussed together with directions for further research.

Lipid Droplets in Cancer: Guardians of Fat in a Stressful World

Cancer cells possess remarkable abilities to adapt to adverse environmental conditions. Their survival during severe nutrient and oxidative stress depends on their capacity to acquire extracellular lipids and the plasticity of their mechanisms for intracellular lipid synthesis, mobilisation, and recycling. Lipid droplets, cytosolic fat storage organelles present in most cells from yeast to men, are emerging as major regulators of lipid metabolism, trafficking, and signalling in various cells and tissues exposed to stress. Their biogenesis is induced by nutrient and oxidative stress and they accumulate in various cancers. Lipid droplets act as switches that coordinate lipid trafficking and consumption for different purposes in the cell, such as energy production, protection against oxidative stress or membrane biogenesis during rapid cell growth. They sequester toxic lipids, such as fatty acids, cholesterol and ceramides, thereby preventing lipotoxic cell damage and engage in a complex relationship with autophagy. Here, we focus on the emerging mechanisms of stress-induced lipid droplet biogenesis; their roles during nutrient, lipotoxic, and oxidative stress; and the relationship between lipid droplets and autophagy. The recently discovered principles of lipid droplet biology can improve our understanding of the mechanisms that govern cancer cell adaptability and resilience to stress.

A practical guide to the handling and administration of personalized transcriptionally attenuated oncolytic adenoviruses (PTAVs)

The aim of this review is to provide practical information on the handling, storage, and administration procedures for personalized oncolytic adenoviruses (PTAVs), which have recently entered clinical trials. As described herein, personalized oncolytic viruses refer to transcriptionally attenuated (TA) type 5 adenoviruses that are engineered to carry one or more neoantigenic transgenes derived from patient tumors. Vials of personalized viruses should be stored at -60°C without refreezing after thawing to maintain infectivity. To prevent accidental exposure and transmission, full implementation of universal precautions for preparation, administration, and handling is required. Contaminated materials that come into contact with personalized viruses should be properly disposed of in accordance with local institutional procedures. Severely immunocompromised or pregnant healthcare workers should not prepare or administer personalized viruses or directly contact injection sites. Personalized viruses are administered subcutaneously and intratumorally; however, only subcutaneous injection will be considered in this review. The specific storage, handling, administration, and safety requirements for personalized viruses are easily managed in the context of a clinical trial following the directives from the study protocol.

Harnessing synthetic lethality to predict the response to cancer treatment

While synthetic lethality (SL) holds promise in developing effective cancer therapies, SL candidates found via experimental screens often have limited translational value. Here we present a data-driven approach, ISLE (identification of clinically relevant synthetic lethality), that mines TCGA cohort to identify the most likely clinically relevant SL interactions (cSLi) from a given candidate set of lab-screened SLi. We first validate ISLE via a benchmark of large-scale drug response screens and by predicting drug efficacy in mouse xenograft models. We then experimentally test a select set of predicted cSLi via new screening experiments, validating their predicted context-specific sensitivity in hypoxic vs normoxic conditions and demonstrating cSLi's utility in predicting synergistic drug combinations. We show that cSLi can successfully predict patients' drug treatment response and provide patient stratification signatures. ISLE thus complements existing actionable mutation-based methods for precision cancer therapy, offering an opportunity to expand its scope to the whole genome.

Fine mapping of epistatic genetic determinants of blood pressure on rat chromosome 5

BACKGROUND

Previously, using linkage analysis and substitution mapping, two closely-linked interactive blood pressure quantitative trait loci (QTLs), BP QTL1 and BP QTL2, were located within a 13.96 Mb region from 117894038 to 131853815 bp (RGSC 3.4 version) on rat chromosome 5 (RNO5). This was done by using a series of congenic strains consisting of genomic segments of the Dahl salt-sensitive (S) rat substituted with that of the normotensive Lewis (LEW) rat. The interactive nature of the two loci was further confirmed by the construction and characterization of a panel of S.LEW bicongenic strains and corresponding S.LEW monocongenic strains, which provided definitive evidence of epistasis (genetic interaction) between BP QTL1 (7.77 Mb) and BP QTL2 (4.18 Mb). The purpose of this work was to further map these interacting QTLs.

METHOD

A new panel of seven new S.LEW bicongenic strains was constructed and characterized for BP.

RESULTS

The data obtained from these new strains further resolved BP QTL1 from 7.77 to 2.93 Mb. Further, BP QTL2 was traceable as not being a single QTL, but a composite of at least three QTLs, LEW alleles at two of which located within 2.26 Mb and 175 kb lowered BP but the third one located within 1.31 Mb increased BP.

CONCLUSION

Lack of coding variation within any of the regions further mapped within the previous QTL2 suggests noncoding variation as likely responsible for the observed epistasis.

Looking beyond drivers and passengers in cancer genome sequencing data

Cancer arises as a result of acquired changes in the DNA sequence of the genome of somatic cells. A subset of the genetic changes, dubbed driver mutations, propels tumor growth, and the remaining changes are passengers, apparently inconsequential for neoplastic transformation. Massive genome sequencing of thousands of tumors from all major cancer types has enabled cataloging of the so-called driver and passenger mutations, and facilitated molecular classification of cancer, guiding precision medicine approach for the patients. Nonetheless, innovative analyses of cancer genomics data has led to novel, sometimes serendipitous findings that have aided to our understanding of other aspects of the biology of the disease and opened up new frontiers. For instance, emerging findings show that mutational patterns in cancer genomes can help detect signatures of known and novel DNA damage and repair processes, provide a likely chronological account of genomic changes in cancer genomes, and allow revisiting the models of cancer evolution. These findings have stimulated original approaches to identify disease etiology, stratify patients, target the disease, and monitor patient responses, complementing driver-mutation centric approaches. In this review, we discuss these emerging approaches and unexpected breakthroughs, and their implications for basic cancer research and clinical practices.

WeSME: uncovering mutual exclusivity of cancer drivers and beyond

Motivation

Mutual exclusivity is a widely recognized property of many cancer drivers. Knowledge about these relationships can provide important insights into cancer drivers, cancer-driving pathways and cancer subtypes. It can also be used to predict new functional interactions between cancer driving genes and uncover novel cancer drivers. Currently, most of mutual exclusivity analyses are preformed focusing on a limited set of genes in part due to the computational cost required to rigorously compute P -values.

Results

To reduce the computing cost and perform less restricted mutual exclusivity analysis, we developed an efficient method to estimate P -values while controlling the mutation rates of individual patients and genes similar to the permutation test. A comprehensive mutual exclusivity analysis allowed us to uncover mutually exclusive pairs, some of which may have relatively low mutation rates. These pairs often included likely cancer drivers that have been missed in previous analyses. More importantly, our results demonstrated that mutual exclusivity can also provide information that goes beyond the interactions between cancer drivers and can, for example, elucidate different mutagenic processes in different cancer groups. In particular, including frequently mutated, long genes such as TTN in our analysis allowed us to observe interesting patterns of APOBEC activity in breast cancer and identify a set of related driver genes that are highly predictive of patient survival. In addition, we utilized our mutual exclusivity analysis in support of a previously proposed model where APOBEC activity is the underlying process that causes TP53 mutations in a subset of breast cancer cases.

Availability and Implementation

http://www.ncbi.nlm.nih.gov/CBBresearch/Przytycka/index.cgi#wesme.

Contact

przytyck@ncbi.nlm.nih.gov.

Supplementary information

Supplementary data are available at Bioinformatics online.

Precision Oncology: The Road Ahead

Current efforts in precision oncology largely focus on the benefit of genomics-guided therapy. Yet, advances in sequencing techniques provide an unprecedented view of the complex genetic and nongenetic heterogeneity within individual tumors. Herein, we outline the benefits of integrating genomic and transcriptomic analyses for advanced precision oncology. We summarize relevant computational approaches to detect novel drivers and genetic vulnerabilities, suitable for therapeutic exploration. Clinically relevant platforms to functionally test predicted drugs/drug combinations for individual patients are reviewed. Finally, we highlight the technological advances in single cell analysis of tumor specimens. These may ultimately lead to the development of next-generation cancer drugs, capable of tackling the hurdles imposed by genetic and phenotypic heterogeneity on current anticancer therapies.

Epistatic interaction between the lipase-encoding genes Pnpla2 and Lipe causes liposarcoma in mice

Liposarcoma is an often fatal cancer of fat cells. Mechanisms of liposarcoma development are incompletely understood. The cleavage of fatty acids from acylglycerols (lipolysis) has been implicated in cancer. We generated mice with adipose tissue deficiency of two major enzymes of lipolysis, adipose triglyceride lipase (ATGL) and hormone-sensitive lipase (HSL), encoded respectively by Pnpla2 and Lipe. Adipocytes from double adipose knockout (DAKO) mice, deficient in both ATGL and HSL, showed near-complete deficiency of lipolysis. All DAKO mice developed liposarcoma between 11 and 14 months of age. No tumors occurred in single knockout or control mice. The transcriptome of DAKO adipose tissue showed marked differences from single knockout and normal controls as early as 3 months. Gpnmb and G0s2 were among the most highly dysregulated genes in premalignant and malignant DAKO adipose tissue, suggesting a potential utility as early markers of the disease. Similar changes of GPNMB and G0S2 expression were present in a human liposarcoma database. These results show that a previously-unknown, fully penetrant epistatic interaction between Pnpla2 and Lipe can cause liposarcoma in mice. DAKO mice provide a promising model for studying early premalignant changes that lead to late-onset malignant disease.

Liposarcoma is an often fatal adult-onset tumor of fat tissue. Lipolysis, the central pathway of fat tissue metabolism, has been implicated in cancer. We generated mice that were deficient in two key enzymes of lipolysis, adipose triglyceride lipase (ATGL) and hormone-sensitive lipase (HSL). Strikingly, all mice with combined ATGL and HSL deficiency developed liposarcoma by 11–14 months of age. No liposarcoma occurred in single knockout or normal controls. Transcriptome analysis revealed that a subset of genes is dysregulated by 3 months of age. Our study reveals a novel epistatic interaction in fat cells between these two lipase genes and that causes a unique form of liposarcoma in mice. The double knockout mice provide a novel tool to study the early stages of liposarcoma development, prognostic markers and preventive treatments.

Tissue-specific tumorigenesis: context matters

How can we treat cancer more effectively? Traditionally, tumours from the same anatomical site are treated as one tumour entity. This concept has been challenged by recent breakthroughs in cancer genomics and translational research that have enabled molecular tumour profiling. The identification and validation of cancer drivers that are shared between different tumour types, spurred the new paradigm to target driver pathways across anatomical sites by off-label drug use, or within so-called basket or umbrella trials which are designed to test whether molecular alterations in one tumour entity can be extrapolated to all others. However, recent clinical and preclinical studies suggest that there are tissue- and cell type-specific differences in tumorigenesis and the organization of oncogenic signalling pathways. In this Opinion article, we focus on the molecular, cellular, systemic and environmental determinants of organ-specific tumorigenesis and the mechanisms of context-specific oncogenic signalling outputs. Investigation, recognition and in-depth biological understanding of these differences will be vital for the design of next-generation clinical trials and the implementation of molecularly guided cancer therapies in the future.

Microscopy-Based High-Content Screening

Image-based screening is used to measure a variety of phenotypes in cells and whole organisms. Combined with perturbations such as RNA interference, small molecules, and mutations, such screens are a powerful method for gaining systematic insights into biological processes. Screens have been applied to study diverse processes, such as protein-localization changes, cancer cell vulnerabilities, and complex organismal phenotypes. Recently, advances in imaging and image-analysis methodologies have accelerated large-scale perturbation screens. Here, we describe the state of the art for image-based screening experiments and delineate experimental approaches and image-analysis approaches as well as discussing challenges and future directions, including leveraging CRISPR/Cas9-mediated genome engineering.

A genetic interaction map of cell cycle regulators

A combination of genome-scale RNA interference screening and genetic interaction analysis using process-directed phenotypes is used to assign components to specific pathways and complexes for modulators of mitosis and cytokinesis in Drosophila S2 cells.

Cell-based RNA interference (RNAi) is a powerful approach to screen for modulators of many cellular processes. However, resulting candidate gene lists from cell-based assays comprise diverse effectors, both direct and indirect, and further dissecting their functions can be challenging. Here we screened a genome-wide RNAi library for modulators of mitosis and cytokinesis in Drosophila S2 cells. The screen identified many previously known genes as well as modulators that have previously not been connected to cell cycle control. We then characterized ∼300 candidate modifiers further by genetic interaction analysis using double RNAi and a multiparametric, imaging-based assay. We found that analyzing cell cycle–relevant phenotypes increased the sensitivity for associating novel gene function. Genetic interaction maps based on mitotic index and nuclear size grouped candidates into known regulatory complexes of mitosis or cytokinesis, respectively, and predicted previously uncharacterized components of known processes. For example, we confirmed a role for the Drosophila CCR4 mRNA processing complex component l(2)NC136 during the mitotic exit. Our results show that the combination of genome-scale RNAi screening and genetic interaction analysis using process-directed phenotypes provides a powerful two-step approach to assigning components to specific pathways and complexes.

Gene network inference by fusing data from diverse distributions

Motivation: Markov networks are undirected graphical models that are widely used to infer relations between genes from experimental data. Their state-of-the-art inference procedures assume the data arise from a Gaussian distribution. High-throughput omics data, such as that from next generation sequencing, often violates this assumption. Furthermore, when collected data arise from multiple related but otherwise nonidentical distributions, their underlying networks are likely to have common features. New principled statistical approaches are needed that can deal with different data distributions and jointly consider collections of datasets.

Results: We present FuseNet, a Markov network formulation that infers networks from a collection of nonidentically distributed datasets. Our approach is computationally efficient and general: given any number of distributions from an exponential family, FuseNet represents model parameters through shared latent factors that define neighborhoods of network nodes. In a simulation study, we demonstrate good predictive performance of FuseNet in comparison to several popular graphical models. We show its effectiveness in an application to breast cancer RNA-sequencing and somatic mutation data, a novel application of graphical models. Fusion of datasets offers substantial gains relative to inference of separate networks for each dataset. Our results demonstrate that network inference methods for non-Gaussian data can help in accurate modeling of the data generated by emergent high-throughput technologies.

Availability and implementation: Source code is at https://github.com/marinkaz/fusenet.

Contact:blaz.zupan@fri.uni-lj.si

Supplementary information:Supplementary information is available at Bioinformatics online.

SNP-SNP interaction analysis of NF-κB signaling pathway on breast cancer survival

In breast cancer, constitutive activation of NF-κB has been reported, however, the impact of genetic variation of the pathway on patient prognosis has been little studied. Furthermore, a combination of genetic variants, rather than single polymorphisms, may affect disease prognosis. Here, in an extensive dataset (n = 30,431) from the Breast Cancer Association Consortium, we investigated the association of 917 SNPs in 75 genes in the NF-κB pathway with breast cancer prognosis. We explored SNP-SNP interactions on survival using the likelihood-ratio test comparing multivariate Cox' regression models of SNP pairs without and with an interaction term. We found two interacting pairs associating with prognosis: patients simultaneously homozygous for the rare alleles of rs5996080 and rs7973914 had worse survival (HRinteraction 6.98, 95% CI=3.3-14.4, P=1.42E-07), and patients carrying at least one rare allele for rs17243893 and rs57890595 had better survival (HRinteraction 0.51, 95% CI=0.3-0.6, P = 2.19E-05). Based on in silico functional analyses and literature, we speculate that the rs5996080 and rs7973914 loci may affect the BAFFR and TNFR1/TNFR3 receptors and breast cancer survival, possibly by disturbing both the canonical and non-canonical NF-κB pathways or their dynamics, whereas, rs17243893-rs57890595 interaction on survival may be mediated through TRAF2-TRAIL-R4 interplay. These results warrant further validation and functional analyses.

Corticosteroid receptor genes and childhood neglect influence susceptibility to crack/cocaine addiction and response to detoxification treatment

The aim of this study was to analyze hypotheses-driven gene-environment and gene-gene interactions in smoked (crack) cocaine addiction by evaluating childhood neglect and polymorphisms in mineralocorticoid and glucocorticoid receptor genes (NR3C2 and NR3C1, respectively). One hundred thirty-nine crack/cocaine-addicted women who completed 3 weeks of follow-up during early abstinence composed our sample. Childhood adversities were assessed using the Childhood Trauma Questionnaire (CTQ), and withdrawal symptoms were assessed using the Cocaine Selective Severity Assessment (CSSA) scale. Conditional logistic regression with counterfactuals and generalized estimating equation modeling were used to test gene-environment and gene-gene interactions. We found an interaction between the rs5522-Val allele and childhood physical neglect, which altered the risk of crack/cocaine addiction (Odds ratio = 4.0, P = 0.001). Moreover, a NR3C2-NR3C1 interaction (P = 0.002) was found modulating the severity of crack/cocaine withdrawal symptoms. In the post hoc analysis, concomitant carriers of the NR3C2 rs5522-Val and NR3C1 rs6198-G alleles showed lower overall severity scores when compared to other genotype groups (P-values ≤ 0.035). This gene-environment interaction is consistent with epidemiological and human experimental findings demonstrating a strong relationship between early life stress and the hypothalamic-pituitary-adrenal (HPA) axis dysregulation in cocaine addiction. Additionally, this study extended in crack/cocaine addiction the findings previously reported for tobacco smoking involving an interaction between NR3C2 and NR3C1 genes.

Concepts in solid tumor evolution

Evolutionary mechanisms in cancer progression give tumors their individuality. Cancer evolution is different from organismal evolution, however, and we discuss where concepts from evolutionary genetics are useful or limited in facilitating an understanding of cancer. Based on these concepts we construct and apply the simplest plausible model of tumor growth and progression. Simulations using this simple model illustrate the importance of stochastic events early in tumorigenesis, highlight the dominance of exponential growth over linear growth and differentiation, and explain the clonal substructure of tumors.