Rewiring of genetic networks in response to DNA damage

The Chromone Alkaloid, Rohitukine, Affords Anti-Cancer Activity via Modulating Apoptosis Pathways in A549 Cell Line and Yeast Mitogen Activated Protein Kinase (MAPK) Pathway

The field of cancer research and treatment has made significant progress, yet we are far from having completely safe, efficient and specific therapies that target cancer cells and spare the healthy tissues. Natural compounds may reduce the problems related to cancer treatment. Currently, many plant products are being used to treat cancer. In this study, Rohitukine, a natural occurring chromone alkaloid extracted from Dysoxylum binectariferum, was investigated for cytotoxic properties against budding yeast as well as against lung cancer (A549) cells. We endeavored to specifically study Rohitukine in S. cerevisiae in the context of MAPK pathways as yeast probably represents the experimental model where the organization and regulation of MAPK pathways are best understood. MAPK are evolutionarily conserved protein kinases that transfer extracellular signals to the machinery controlling essential cellular processes like growth, migration, differentiation, cell division and apoptosis. We aimed at carrying out hypothesis driven studies towards targeting the important network of cellular communication, a critical process that gets awry in cancer. Employing mutant strains of genetic model system Saccharomyces cerevisiae. S. cerevisiae encodes five MAPKs involved in control of distinct cellular responses such as growth, differentiation, migration and apoptosis. Our study involves gene knockouts of Slt2 and Hog1 which are functional homologs of human ERK5 and mammalian p38 MAPK, respectively. We performed cytotoxicity assay to evaluate the effect of Rohitukine on cell viability and also determined the effects of drug on generation of reactive oxygen species, induction of apoptosis and expression of Slt2 and Hog1 gene at mRNA level in the presence of drug. The results of this study show a differential effect in the activity of drug between the WT, Slt2 and Hog1 gene deletion strain indicating involvement of MAPK pathway. Further, we investigated Rohitukine induced cytotoxic effects in lung cancer cells and stimulated the productions of ROS after exposure for 24 hrs. Results from western blotting suggest that Rohitukine triggered apoptosis in A549 cell line through upregulation of p53, caspase9 and down regulation of Bcl-2 protein. The scope of this study is to understand the mechanism of anticancer activity of Rohitukine to increase the repertoire of anticancer drugs, so that problem created by emergence of resistance towards standard anticancer compounds can be alleviated.

Identification of Biomarker and Co-Regulatory Motifs in Lung Adenocarcinoma Based on Differential Interactions

Changes in intermolecular interactions (differential interactions) may influence the progression of cancer. Specific genes and their regulatory networks may be more closely associated with cancer when taking their transcriptional and post-transcriptional levels and dynamic and static interactions into account simultaneously. In this paper, a differential interaction analysis was performed to detect lung adenocarcinoma-related genes. Furthermore, a miRNA-TF (transcription factor) synergistic regulation network was constructed to identify three kinds of co-regulated motifs, namely, triplet, crosstalk and joint. Not only were the known cancer-related miRNAs and TFs (let-7, miR-15a, miR-17, TP53, ETS1, and so on) were detected in the motifs, but also the miR-15, let-7 and miR-17 families showed a tendency to regulate the triplet, crosstalk and joint motifs, respectively. Moreover, several biological functions (i.e., cell cycle, signaling pathways and hemopoiesis) associated with the three motifs were found to be frequently targeted by the drugs for lung adenocarcinoma. Specifically, the two 4-node motifs (crosstalk and joint) based on co-expression and interaction had a closer relationship to lung adenocarcinoma, and so further research was performed on them. A 10-gene biomarker (UBC, SRC, SP1, MYC, STAT3, JUN, NR3C1, RB1, GRB2 and MAPK1) was selected from the joint motif, and a survival analysis indicated its significant association with survival. Among the ten genes, JUN, NR3C1 and GRB2 are our newly detected candidate lung adenocarcinoma-related genes. The genes, regulators and regulatory motifs detected in this work will provide potential drug targets and new strategies for individual therapy.

Functional Analysis and Characterization of Differential Coexpression Networks

Differential coexpression analysis is emerging as a complement to conventional differential gene expression analysis. The identified differential coexpression links can be assembled into a differential coexpression network (DCEN) in response to environmental stresses or genetic changes. Differential coexpression analyses have been successfully used to identify condition-specific modules; however, the structural properties and biological significance of general DCENs have not been well investigated. Here, we analyzed two independent Saccharomyces cerevisiae DCENs constructed from large-scale time-course gene expression profiles in response to different situations. Topological analyses show that DCENs are tree-like networks possessing scale-free characteristics, but not small-world. Functional analyses indicate that differentially coexpressed gene pairs in DCEN tend to link different biological processes, achieving complementary or synergistic effects. Furthermore, the gene pairs lacking common transcription factors are sensitive to perturbation and hence lead to differential coexpression. Based on these observations, we integrated transcriptional regulatory information into DCEN and identified transcription factors that might cause differential coexpression by gain or loss of activation in response to different situations. Collectively, our results not only uncover the unique structural characteristics of DCEN but also provide new insights into interpretation of DCEN to reveal its biological significance and infer the underlying gene regulatory dynamics.

Algorithms for modeling global and context-specific functional relationship networks

Functional genomics has enormous potential to facilitate our understanding of normal and disease-specific physiology. In the past decade, intensive research efforts have been focused on modeling functional relationship networks, which summarize the probability of gene co-functionality relationships. Such modeling can be based on either expression data only or heterogeneous data integration. Numerous methods have been deployed to infer the functional relationship networks, while most of them target the global (non-context-specific) functional relationship networks. However, it is expected that functional relationships consistently reprogram under different tissues or biological processes. Thus, advanced methods have been developed targeting tissue-specific or developmental stage-specific networks. This article brings together the state-of-the-art functional relationship network modeling methods, emphasizes the need for heterogeneous genomic data integration and context-specific network modeling and outlines future directions for functional relationship networks.

The cancer cell map initiative: defining the hallmark networks of cancer

Progress in DNA sequencing has revealed the startling complexity of cancer genomes, which typically carry thousands of somatic mutations. However, it remains unclear which are the key driver mutations or dependencies in a given cancer and how these influence pathogenesis and response to therapy. Although tumors of similar types and clinical outcomes can have patterns of mutations that are strikingly different, it is becoming apparent that these mutations recurrently hijack the same hallmark molecular pathways and networks. For this reason, it is likely that successful interpretation of cancer genomes will require comprehensive knowledge of the molecular networks under selective pressure in oncogenesis. Here we announce the creation of a new effort, The Cancer Cell Map Initiative (CCMI), aimed at systematically detailing these complex interactions among cancer genes and how they differ between diseased and healthy states. We discuss recent progress that enables creation of these cancer cell maps across a range of tumor types and how they can be used to target networks disrupted in individual patients, significantly accelerating the development of precision medicine.

Cancer type-dependent genetic interactions between cancer driver alterations indicate plasticity of epistasis across cell types

Cancers, like many diseases, are normally caused by combinations of genetic alterations rather than by changes affecting single genes. It is well established that the genetic alterations that drive cancer often interact epistatically, having greater or weaker consequences in combination than expected from their individual effects. In a stringent statistical analysis of data from > 3,000 tumors, we find that the co-occurrence and mutual exclusivity relationships between cancer driver alterations change quite extensively in different types of cancer. This cannot be accounted for by variation in tumor heterogeneity or unrecognized cancer subtypes. Rather, it suggests that how genomic alterations interact cooperatively or partially redundantly to driver cancer changes in different types of cancers. This re-wiring of epistasis across cell types is likely to be a basic feature of genetic architecture, with important implications for understanding the evolution of multicellularity and human genetic diseases. In addition, if this plasticity of epistasis across cell types is also true for synthetic lethal interactions, a synthetic lethal strategy to kill cancer cells may frequently work in one type of cancer but prove ineffective in another.

DINGO: differential network analysis in genomics

MOTIVATION

Cancer progression and development are initiated by aberrations in various molecular networks through coordinated changes across multiple genes and pathways. It is important to understand how these networks change under different stress conditions and/or patient-specific groups to infer differential patterns of activation and inhibition. Existing methods are limited to correlation networks that are independently estimated from separate group-specific data and without due consideration of relationships that are conserved across multiple groups.

METHOD

We propose a pathway-based differential network analysis in genomics (DINGO) model for estimating group-specific networks and making inference on the differential networks. DINGO jointly estimates the group-specific conditional dependencies by decomposing them into global and group-specific components. The delineation of these components allows for a more refined picture of the major driver and passenger events in the elucidation of cancer progression and development.

RESULTS

Simulation studies demonstrate that DINGO provides more accurate group-specific conditional dependencies than achieved by using separate estimation approaches. We apply DINGO to key signaling pathways in glioblastoma to build differential networks for long-term survivors and short-term survivors in The Cancer Genome Atlas. The hub genes found by mRNA expression, DNA copy number, methylation and microRNA expression reveal several important roles in glioblastoma progression.

AVAILABILITY AND IMPLEMENTATION

R Package at: odin.mdacc.tmc.edu/∼vbaladan.

CONTACT

veera@mdanderson.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Conserved Circular Network of Coregulated Lipids Modulates Innate Immune Responses

Lipid composition affects the biophysical properties of membranes that provide a platform for receptor-mediated cellular signaling. To study the regulatory role of membrane lipid composition, we combined genetic perturbations of sphingolipid metabolism with the quantification of diverse steps in Toll-like receptor (TLR) signaling and mass spectrometry-based lipidomics. Membrane lipid composition was broadly affected by these perturbations, revealing a circular network of coregulated sphingolipids and glycerophospholipids. This evolutionarily conserved network architecture simultaneously reflected membrane lipid metabolism, subcellular localization, and adaptation mechanisms. Integration of the diverse TLR-induced inflammatory phenotypes with changes in lipid abundance assigned distinct functional roles to individual lipid species organized across the network. This functional annotation accurately predicted the inflammatory response of cells derived from patients suffering from lipid storage disorders, based solely on their altered membrane lipid composition. The analytical strategy described here empowers the understanding of higher-level organization of membrane lipid function in diverse biological systems.

Revealing Pathway Dynamics in Heart Diseases by Analyzing Multiple Differential Networks

Development of heart diseases is driven by dynamic changes in both the activity and connectivity of gene pathways. Understanding these dynamic events is critical for understanding pathogenic mechanisms and development of effective treatment. Currently, there is a lack of computational methods that enable analysis of multiple gene networks, each of which exhibits differential activity compared to the network of the baseline/healthy condition. We describe the iMDM algorithm to identify both unique and shared gene modules across multiple differential co-expression networks, termed M-DMs (multiple differential modules). We applied iMDM to a time-course RNA-Seq dataset generated using a murine heart failure model generated on two genotypes. We showed that iMDM achieves higher accuracy in inferring gene modules compared to using single or multiple co-expression networks. We found that condition-specific M-DMs exhibit differential activities, mediate different biological processes, and are enriched for genes with known cardiovascular phenotypes. By analyzing M-DMs that are present in multiple conditions, we revealed dynamic changes in pathway activity and connectivity across heart failure conditions. We further showed that module dynamics were correlated with the dynamics of disease phenotypes during the development of heart failure. Thus, pathway dynamics is a powerful measure for understanding pathogenesis. iMDM provides a principled way to dissect the dynamics of gene pathways and its relationship to the dynamics of disease phenotype. With the exponential growth of omics data, our method can aid in generating systems-level insights into disease progression.

A sensitised RNAi screen reveals a ch-TOG genetic interaction network required for spindle assembly

How multiple spindle assembly pathways are integrated to drive bipolar spindle assembly is poorly understood. We performed an image-based double RNAi screen to identify genes encoding Microtubule-Associated Proteins (MAPs) that interact with the highly conserved ch-TOG gene to regulate bipolar spindle assembly in human cells. We identified a ch-TOG centred network of genetic interactions which promotes ensures centrosome-mediated microtubule polymerisation, leading to the incorporation of microtubules polymerised by all pathways into a bipolar structure. Our genetic screen also reveals that ch-TOG maintains a dynamic microtubule population, in part, through modulating HSET activity. ch-TOG ensures that spindle assembly is robust to perturbation but sufficiently dynamic such that spindles can explore a diverse shape space in search of structures that can align chromosomes.

Integrative network-based approach identifies key genetic elements in breast invasive carcinoma

BACKGROUND

Breast cancer is a genetically heterogeneous type of cancer that belongs to the most prevalent types with a high mortality rate. Treatment and prognosis of breast cancer would profit largely from a correct classification and identification of genetic key drivers and major determinants driving the tumorigenesis process. In the light of the availability of tumor genomic and epigenomic data from different sources and experiments, new integrative approaches are needed to boost the probability of identifying such genetic key drivers. We present here an integrative network-based approach that is able to associate regulatory network interactions with the development of breast carcinoma by integrating information from gene expression, DNA methylation, miRNA expression, and somatic mutation datasets.

RESULTS

Our results showed strong association between regulatory elements from different data sources in terms of the mutual regulatory influence and genomic proximity. By analyzing different types of regulatory interactions, TF-gene, miRNA-mRNA, and proximity analysis of somatic variants, we identified 106 genes, 68 miRNAs, and 9 mutations that are candidate drivers of oncogenic processes in breast cancer. Moreover, we unraveled regulatory interactions among these key drivers and the other elements in the breast cancer network. Intriguingly, about one third of the identified driver genes are targeted by known anti-cancer drugs and the majority of the identified key miRNAs are implicated in cancerogenesis of multiple organs. Also, the identified driver mutations likely cause damaging effects on protein functions. The constructed gene network and the identified key drivers were compared to well-established network-based methods.

CONCLUSION

The integrated molecular analysis enabled by the presented network-based approach substantially expands our knowledge base of prospective genomic drivers of genes, miRNAs, and mutations. For a good part of the identified key drivers there exists solid evidence for involvement in the development of breast carcinomas. Our approach also unraveled the complex regulatory interactions comprising the identified key drivers. These genomic drivers could be further investigated in the wet lab as potential candidates for new drug targets. This integrative approach can be applied in a similar fashion to other cancer types, complex diseases, or for studying cellular differentiation processes.

Stress-induced DNA damage: a case study in diffuse large B-cell lymphoma

DNA damage is one of the mechanisms of mutagenesis. Sequence integrity may be affected by the action of thermal changes, chemical agents, both endogenous and exogenous, and other environmental issues. Abnormally high mutation rates are referred to as genomic instability: a phenomenon closely related to the onset of cancer. Mutant genotypes may be able to confer some kind of selective advantage on subclonal cell populations, leading them to multiply until dominance in a localized tissue environment that later becomes the tumour. Cellular stress, especially that of oxidative and ionic nature, is a recognized trigger for DNA-damaging processes. A physico-chemical model has shown that high hysteresis rates in DNA denaturation curves may be indicative of dissipative processes inducing DNA damage, thus potentially leading to uncontrolled mutagenesis and genome instability. We here study selectively to what extent this phenomenon may occur by analysing the sequence length and composition effects on the thermodynamic behaviour and the presence of hysteresis in pressure-driven DNA denaturation; pronounced hysteresis in the denaturation/renaturation curves may indicate thermal susceptibility to DNA damage. In particular, we consider highly mutated regions of the genome characterized in diffuse large B-cell lymphoma on a recent whole exome next-generation sequencing effort.

Increased signaling entropy in cancer requires the scale-free property of protein interaction networks

One of the key characteristics of cancer cells is an increased phenotypic plasticity, driven by underlying genetic and epigenetic perturbations. However, at a systems-level it is unclear how these perturbations give rise to the observed increased plasticity. Elucidating such systems-level principles is key for an improved understanding of cancer. Recently, it has been shown that signaling entropy, an overall measure of signaling pathway promiscuity, and computable from integrating a sample's gene expression profile with a protein interaction network, correlates with phenotypic plasticity and is increased in cancer compared to normal tissue. Here we develop a computational framework for studying the effects of network perturbations on signaling entropy. We demonstrate that the increased signaling entropy of cancer is driven by two factors: (i) the scale-free (or near scale-free) topology of the interaction network, and (ii) a subtle positive correlation between differential gene expression and node connectivity. Indeed, we show that if protein interaction networks were random graphs, described by Poisson degree distributions, that cancer would generally not exhibit an increased signaling entropy. In summary, this work exposes a deep connection between cancer, signaling entropy and interaction network topology.

Differential genetic interactions of yeast stress response MAPK pathways

Genetic interaction screens have been applied with great success in several organisms to study gene function and the genetic architecture of the cell. However, most studies have been performed under optimal growth conditions even though many functional interactions are known to occur under specific cellular conditions. In this study, we have performed a large-scale genetic interaction analysis in Saccharomyces cerevisiae involving approximately 49 × 1,200 double mutants in the presence of five different stress conditions, including osmotic, oxidative and cell wall-altering stresses. This resulted in the generation of a differential E-MAP (or dE-MAP) comprising over 250,000 measurements of conditional interactions. We found an extensive number of conditional genetic interactions that recapitulate known stress-specific functional associations. Furthermore, we have also uncovered previously unrecognized roles involving the phosphatase regulator Bud14, the histone methylation complex COMPASS and membrane trafficking complexes in modulating the cell wall integrity pathway. Finally, the osmotic stress differential genetic interactions showed enrichment for genes coding for proteins with conditional changes in phosphorylation but not for genes with conditional changes in gene expression. This suggests that conditional genetic interactions are a powerful tool to dissect the functional importance of the different response mechanisms of the cell.

Identification of ncRNAs as potential therapeutic targets in multiple sclerosis through differential ncRNA - mRNA network analysis

Background

Several studies have revealed a potential role for both small nucleolar RNAs (snoRNAs) and microRNAs (miRNAs) in the physiopathology of relapsing-remitting multiple sclerosis (RRMS). This potential implication has been mainly described through differential expression studies. However, it has been suggested that, in order to extract additional information from large-scale expression experiments, differential expression studies must be complemented with differential network studies. Thus, the present work is aimed at the identification of potential therapeutic ncRNA targets for RRMS through differential network analysis of ncRNA – mRNA coexpression networks. ncRNA – mRNA coexpression networks have been constructed from both selected ncRNA (specifically miRNAs, snoRNAs and sdRNAs) and mRNA large-scale expression data obtained from 22 patients in relapse, the same 22 patients in remission and 22 healthy controls. Condition-specific (relapse, remission and healthy) networks have been built and compared to identify the parts of the system most affected by perturbation and aid the identification of potential therapeutic targets among the ncRNAs.

Results

All the coexpression networks we built present a scale-free topology and many snoRNAs are shown to have a prominent role in their architecture. The differential network analysis (relapse vs. remission vs. controls’ networks) has revealed that, although both network topology and the majority of the genes are maintained, few ncRNA – mRNA links appear in more than one network. We have selected as potential therapeutic targets the ncRNAs that appear in the disease-specific network and were found to be differentially expressed in a previous study.

Conclusions

Our results suggest that the diseased state of RRMS has a strong impact on the ncRNA – mRNA network of peripheral blood leukocytes, as a massive rewiring of the network happens between conditions. Our findings also indicate that the role snoRNAs have in targeted gene silencing is a widespread phenomenon. Finally, among the potential therapeutic target ncRNAs, SNORA40 seems to be the most promising candidate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1396-5) contains supplementary material, which is available to authorized users.

Extracting high confidence protein interactions from affinity purification data: at the crossroads

Deriving protein-protein interactions from data generated by affinity-purification and mass spectrometry (AP-MS) techniques requires application of scoring methods to measure the reliability of detected putative interactions. Choosing the appropriate scoring method has become a major challenge. Here we apply six popular scoring methods to the same AP-MS dataset and compare their performance. The comparison was carried out for six distinct datasets from human, fly and yeast, which focus on different biological processes and differ in their coverage of the proteome. Results show that the performance of a given scoring method may vary substantially depending on the dataset. Disturbingly, we find that the high confidence (HC) PPI networks built by applying the six scoring methods to the same raw AP-MS dataset display very poor overlap, with only 1.7-4.1% of the HC interactions present in all the networks built, respectively, from the proteome-wide human, fly or yeast datasets. Various properties of the shared versus unique interactions in each network, including biases in protein abundance, suggest that current scoring methods are able to eliminate only the most obvious contaminants, but still fail to reliably single out specific interactions from the large body of spurious associations detected in the AP-MS experiments.

BIOLOGICAL SIGNIFICANCE

The fast progress in AP-MS techniques has prompted the development of a multitude of scoring methods, which are relied upon to remove contaminants and non-specific binders. Choosing the appropriate scoring scheme for a given AP-MS dataset has become a major challenge. The comparative analysis of 6 of the most popular scoring methods, presented here, reveals that overall these methods do not perform as expected. Evidence is provided that this is due to 3 closely related issues: the high 'noise' levels of the raw AP-MS data, the limited capacity of current scoring methods to deal with such high noise levels, and the biases introduced using Gold Standard datasets to benchmark the scoring functions and threshold the networks. For the field to move forward, all three issues will have to be addressed. This article is part of a Special Issue entitled: Protein dynamics in health and disease. Guest Editors: Pierre Thibault and Anne-Claude Gingras.

Fitness profiling links topoisomerase II regulation of centromeric integrity to doxorubicin resistance in fission yeast

Doxorubicin, a chemotherapeutic agent, inhibits the religation step of topoisomerase II (Top2). However, the downstream ramifications of this action are unknown. Here we performed epistasis analyses of top2 with 63 genes representing doxorubicin resistance (DXR) genes in fission yeast and revealed a subset that synergistically collaborate with Top2 to confer DXR. Our findings show that the chromatin-regulating RSC and SAGA complexes act with Top2 in a cluster that is functionally distinct from the Ino80 complex. In various DXR mutants, doxorubicin hypersensitivity was unexpectedly suppressed by a concomitant top2 mutation. Several DXR proteins showed centromeric localization, and their disruption resulted in centromeric defects and chromosome missegregation. An additional top2 mutation could restore centromeric chromatin integrity, suggesting a counterbalance between Top2 and these DXR factors in conferring doxorubicin resistance. Overall, this molecular basis for mitotic catastrophe associated with doxorubicin treatment will help to facilitate drug combinatorial usage in doxorubicin-related chemotherapeutic regimens.

How to train your microbe: methods for dynamically characterizing gene networks

Gene networks regulate biological processes dynamically. However, researchers have largely relied upon static perturbations, such as growth media variations and gene knockouts, to elucidate gene network structure and function. Thus, much of the regulation on the path from DNA to phenotype remains poorly understood. Recent studies have utilized improved genetic tools, hardware, and computational control strategies to generate precise temporal perturbations outside and inside of live cells. These experiments have, in turn, provided new insights into the organizing principles of biology. Here, we introduce the major classes of dynamical perturbations that can be used to study gene networks, and discuss technologies available for creating them in a wide range of microbial pathways.

Differential co-expression network centrality and machine learning feature selection for identifying susceptibility hubs in networks with scale-free structure

Background

Biological insights into group differences, such as disease status, have been achieved through differential co-expression analysis of microarray data. Additional understanding of group differences may be achieved by integrating the connectivity structure of the differential co-expression network and per-gene differential expression between phenotypic groups. Such a global differential co-expression network strategy may increase sensitivity to detect gene-gene interactions (or expression epistasis) that may act as candidates for rewiring susceptibility co-expression networks.

Methods

We test two methods for inferring Genetic Association Interaction Networks (GAIN) incorporating both differential co-expression effects and differential expression effects: a generalized linear model (GLM) regression method with interaction effects (reGAIN) and a Fisher test method for correlation differences (dcGAIN). We rank the importance of each gene with complete interaction network centrality (CINC), which integrates each gene’s differential co-expression effects in the GAIN model along with each gene’s individual differential expression measure. We compare these methods with statistical learning methods Relief-F, Random Forests and Lasso. We also develop a mixture model and permutation approach for determining significant importance score thresholds for network centralities, Relief-F and Random Forest. We introduce a novel simulation strategy that generates microarray case–control data with embedded differential co-expression networks and underlying correlation structure based on scale-free or Erdos-Renyi (ER) random networks.

Results

Using the network simulation strategy, we find that Relief-F and reGAIN provide the best balance between detecting interactions and main effects, plus reGAIN has the ability to adjust for covariates and model quantitative traits. The dcGAIN approach performs best at finding differential co-expression effects by design but worst for main effects, and it does not adjust for covariates and is limited to dichotomous outcomes. When the underlying network is scale free instead of ER, all interaction network methods have greater power to find differential co-expression effects. We apply these methods to a public microarray study of the differential immune response to influenza vaccine, and we identify effects that suggest a role in influenza vaccine immune response for genes from the PI3K family, which includes genes with known immunodeficiency function, and KLRG1, which is a known marker of senescence.

Electronic supplementary material

The online version of this article (doi:10.1186/s13040-015-0040-x) contains supplementary material, which is available to authorized users.

Linking tumor mutations to drug responses via a quantitative chemical-genetic interaction map

There is an urgent need in oncology to link molecular aberrations in tumors with therapeutics that can be administered in a personalized fashion. One approach identifies synthetic-lethal genetic interactions or dependencies that cancer cells acquire in the presence of specific mutations. Using engineered isogenic cells, we generated a systematic and quantitative chemical-genetic interaction map that charts the influence of 51 aberrant cancer genes on 90 drug responses. The dataset strongly predicts drug responses found in cancer cell line collections, indicating that isogenic cells can model complex cellular contexts. Applying this dataset to triple-negative breast cancer, we report clinically actionable interactions with the MYC oncogene, including resistance to AKT-PI3K pathway inhibitors and an unexpected sensitivity to dasatinib through LYN inhibition in a synthetic lethal manner, providing new drug and biomarker pairs for clinical investigation. This scalable approach enables the prediction of drug responses from patient data and can accelerate the development of new genotype-directed therapies.

SIGNIFICANCE

Determining how the plethora of genomic abnormalities that exist within a given tumor cell affects drug responses remains a major challenge in oncology. Here, we develop a new mapping approach to connect cancer genotypes to drug responses using engineered isogenic cell lines and demonstrate how the resulting dataset can guide clinical interrogation.

Dynamic epistasis under varying environmental perturbations

Epistasis describes the phenomenon that mutations at different loci do not have independent effects with regard to certain phenotypes. Understanding the global epistatic landscape is vital for many genetic and evolutionary theories. Current knowledge for epistatic dynamics under multiple conditions is limited by the technological difficulties in experimentally screening epistatic relations among genes. We explored this issue by applying flux balance analysis to simulate epistatic landscapes under various environmental perturbations. Specifically, we looked at gene-gene epistatic interactions, where the mutations were assumed to occur in different genes. We predicted that epistasis tends to become more positive from glucose-abundant to nutrient-limiting conditions, indicating that selection might be less effective in removing deleterious mutations in the latter. We also observed a stable core of epistatic interactions in all tested conditions, as well as many epistatic interactions unique to each condition. Interestingly, genes in the stable epistatic interaction network are directly linked to most other genes whereas genes with condition-specific epistasis form a scale-free network. Furthermore, genes with stable epistasis tend to have similar evolutionary rates, whereas this co-evolving relationship does not hold for genes with condition-specific epistasis. Our findings provide a novel genome-wide picture about epistatic dynamics under environmental perturbations.

Yeast PP4 interacts with ATR homolog Ddc2-Mec1 and regulates checkpoint signaling

Mec1-Ddc2 (ATR-ATRIP) controls the DNA damage checkpoint and shows differential cell-cycle regulation in yeast. To find regulators of Mec1-Ddc2, we exploited a mec1 mutant that retains catalytic activity in G2 and recruitment to stalled replication forks, but which is compromised for the intra-S phase checkpoint. Two screens, one for spontaneous survivors and an E-MAP screen for synthetic growth effects, identified loss of PP4 phosphatase, pph3Δ and psy2Δ, as the strongest suppressors of mec1-100 lethality on HU. Restored Rad53 phosphorylation accounts for part, but not all, of the pph3Δ-mediated survival. Phosphoproteomic analysis confirmed that 94% of the mec1-100-compromised targets on HU are PP4 regulated, including a phosphoacceptor site within Mec1 itself, mutation of which confers damage sensitivity. Physical interaction between Pph3 and Mec1, mediated by cofactors Psy2 and Ddc2, is shown biochemically and through FRET in subnuclear repair foci. This establishes a physical and functional Mec1-PP4 unit for regulating the checkpoint response.

Genetic Interaction Scoring Procedure for Bacterial Species

A genetic interaction occurs when the phenotype of an organism carrying two mutant genes differs from what should have been observed given their independent influence. Such unexpected outcome indicates a mechanistic connection between the perturbed genes, providing a key source of functional information about the cell. Large-scale screening for genetic interactions involves measuring phenotypes of single and double mutants, which for microorganisms is usually done by automated analysis of images of ordered colonies. Obtaining accurate colony sizes, and using them to identify genetic interactions from such screens remains a challenging and time-consuming task. Here, we outline steps to compute genetic interaction scores in E. coli by measuring colony sizes from plate images, performing normalisation, and quantifying the strength of the effect.

Single-cell and multivariate approaches in genetic perturbation screens

Large-scale genetic perturbation screens are a classical approach in biology and have been crucial for many discoveries. New technologies can now provide unbiased quantification of multiple molecular and phenotypic changes across tens of thousands of individual cells from large numbers of perturbed cell populations simultaneously. In this Review, we describe how these developments have enabled the discovery of new principles of intracellular and intercellular organization, novel interpretations of genetic perturbation effects and the inference of novel functional genetic interactions. These advances now allow more accurate and comprehensive analyses of gene function in cells using genetic perturbation screens.

Predictive medicine: outcomes, challenges and opportunities in the Synergy-COPD project

BACKGROUND

Chronic Obstructive Pulmonary Disease (COPD) is a major challenge for healthcare. Heterogeneities in clinical manifestations and in disease progression are relevant traits in COPD with impact on patient management and prognosis. It is hypothesized that COPD heterogeneity results from the interplay of mechanisms governing three conceptually different phenomena: 1) pulmonary disease, 2) systemic effects of COPD and 3) co-morbidity clustering.

OBJECTIVES

To assess the potential of systems medicine to better understand non-pulmonary determinants of COPD heterogeneity. To transfer acquired knowledge to healthcare enhancing subject-specific health risk assessment and stratification to improve management of chronic patients.

METHOD

Underlying mechanisms of skeletal muscle dysfunction and of co-morbidity clustering in COPD patients were explored with strategies combining deterministic modelling and network medicine analyses using the Biobridge dataset. An independent data driven analysis of co-morbidity clustering examining associated genes and pathways was done (ICD9-CM data from Medicare, 13 million people). A targeted network analysis using the two studies: skeletal muscle dysfunction and co-morbidity clustering explored shared pathways between them.

RESULTS

(1) Evidence of abnormal regulation of pivotal skeletal muscle biological pathways and increased risk for co-morbidity clustering was observed in COPD; (2) shared abnormal pathway regulation between skeletal muscle dysfunction and co-morbidity clustering; and, (3) technological achievements of the projects were: (i) COPD Knowledge Base; (ii) novel modelling approaches; (iii) Simulation Environment; and, (iv) three layers of Clinical Decision Support Systems.

CONCLUSIONS

The project demonstrated the high potential of a systems medicine approach to address COPD heterogeneity. Limiting factors for the project development were identified. They were relevant to shape strategies fostering 4P Medicine for chronic patients. The concept of Digital Health Framework and the proposed roadmap for its deployment constituted relevant project outcomes.

Identifying functional gene regulatory network phenotypes underlying single cell transcriptional variability

Recent analysis of single-cell transcriptomic data has revealed a surprising organization of the transcriptional variability pervasive across individual neurons. In response to distinct combinations of synaptic input-type, a new organization of neuronal subtypes emerged based on transcriptional states that were aligned along a gradient of correlated gene expression. Individual neurons traverse across these transcriptional states in response to cellular inputs. However, the regulatory network interactions driving these changes remain unclear. Here we present a novel fuzzy logic-based approach to infer quantitative gene regulatory network models from highly variable single-cell gene expression data. Our approach involves developing an a priori regulatory network that is then trained against in vivo single-cell gene expression data in order to identify causal gene interactions and corresponding quantitative model parameters. Simulations of the inferred gene regulatory network response to experimentally observed stimuli levels mirrored the pattern and quantitative range of gene expression across individual neurons remarkably well. In addition, the network identification results revealed that distinct regulatory interactions, coupled with differences in the regulatory network stimuli, drive the variable gene expression patterns observed across the neuronal subtypes. We also identified a key difference between the neuronal subtype-specific networks with respect to negative feedback regulation, with the catecholaminergic subtype network lacking such interactions. Furthermore, by varying regulatory network stimuli over a wide range, we identified several cases in which divergent neuronal subtypes could be driven towards similar transcriptional states by distinct stimuli operating on subtype-specific regulatory networks. Based on these results, we conclude that heterogenous single-cell gene expression profiles should be interpreted through a regulatory network modeling perspective in order to separate the contributions of network interactions from those of cellular inputs.

Rewiring of genetic networks in response to modification of genetic background

Genome-scale genetic interaction networks are progressively contributing to map the molecular circuitry that determines cellular behavior. To what extent this mapping changes in response to different environmental or genetic conditions is, however, largely unknown. Here, we assembled a genetic network using an in silico model of metabolism in yeast to explicitly ask how separate genetic backgrounds alter network structure. Backgrounds defined by single deletions of metabolically active enzymes induce strong rewiring when the deletion corresponds to a catabolic gene, evidencing a broad redistribution of fluxes to alternative pathways. We also show how change is more pronounced in interactions linking genes in distinct functional modules and in those connections that present weak epistasis. These patterns reflect overall the distributed robustness of catabolism. In a second class of genetic backgrounds, in which a number of neutral mutations accumulate, we dominantly observe modifications in the negative interactions that together with an increase in the number of essential genes indicate a global reduction in buffering. Notably, neutral trajectories that originate considerable changes in the wild-type network comprise mutations that diminished the environmental plasticity of the corresponding metabolism, what emphasizes a mechanistic integration of genetic and environmental buffering. More generally, our work demonstrates how the specific mechanistic causes of robustness influence the architecture of multiconditional genetic interaction maps.

Widespread genetic epistasis among cancer genes

Quantitative genetic epistasis has been hypothesized to be an important factor in the development and progression of complex diseases. Cancers in particular are driven by the accumulation of mutations that may act epistatically during the course of the disease. However, as cancer mutations are uncovered at an unprecedented rate, determining which combinations of genetic alterations interact to produce cancer phenotypes remains a challenge. Here we show that by using combinatorial RNAi screening in cell culture, dense and often previously undetermined interactions among cancer genes were revealed by assessing gene pairs that are frequently co-altered in primary breast cancers. These interacting gene pairs are significantly associated with survival time when co-altered in patients, indicating that genetic interaction mapping may be leveraged to improve risk assessment. As many of these interacting gene pairs involve known drug targets, personalized treatment regimens may be improved by overlaying genetic interactions with mutational profiling.

Integration of Network Biology and Imaging to Study Cancer Phenotypes and Responses

Ever growing "omics" data and continuously accumulated biological knowledge provide an unprecedented opportunity to identify molecular biomarkers and their interactions that are responsible for cancer phenotypes that can be accurately defined by clinical measurements such as in vivo imaging. Since signaling or regulatory networks are dynamic and context-specific, systematic efforts to characterize such structural alterations must effectively distinguish significant network rewiring from random background fluctuations. Here we introduced a novel integration of network biology and imaging to study cancer phenotypes and responses to treatments at the molecular systems level. Specifically, Differential Dependence Network (DDN) analysis was used to detect statistically significant topological rewiring in molecular networks between two phenotypic conditions, and in vivo Magnetic Resonance Imaging (MRI) was used to more accurately define phenotypic sample groups for such differential analysis. We applied DDN to analyze two distinct phenotypic groups of breast cancer and study how genomic instability affects the molecular network topologies in high-grade ovarian cancer. Further, FDA-approved arsenic trioxide (ATO) and the ND2-SmoA1 mouse model of Medulloblastoma (MB) were used to extend our analyses of combined MRI and Reverse Phase Protein Microarray (RPMA) data to assess tumor responses to ATO and to uncover the complexity of therapeutic molecular biology.

Gene coexpression measures in large heterogeneous samples using count statistics

With the advent of high-throughput technologies making large-scale gene expression data readily available, developing appropriate computational tools to process these data and distill insights into systems biology has been an important part of the "big data" challenge. Gene coexpression is one of the earliest techniques developed that is still widely in use for functional annotation, pathway analysis, and, most importantly, the reconstruction of gene regulatory networks, based on gene expression data. However, most coexpression measures do not specifically account for local features in expression profiles. For example, it is very likely that the patterns of gene association may change or only exist in a subset of the samples, especially when the samples are pooled from a range of experiments. We propose two new gene coexpression statistics based on counting local patterns of gene expression ranks to take into account the potentially diverse nature of gene interactions. In particular, one of our statistics is designed for time-course data with local dependence structures, such as time series coupled over a subregion of the time domain. We provide asymptotic analysis of their distributions and power, and evaluate their performance against a wide range of existing coexpression measures on simulated and real data. Our new statistics are fast to compute, robust against outliers, and show comparable and often better general performance.

DiffNet: automatic differential functional summarization of dE-MAP networks

The study of genetic interaction networks that respond to changing conditions is an emerging research problem. Recently, Bandyopadhyay et al. (2010) proposed a technique to construct a differential network (dE-MAPnetwork) from two static gene interaction networks in order to map the interaction differences between them under environment or condition change (e.g., DNA-damaging agent). This differential network is then manually analyzed to conclude that DNA repair is differentially effected by the condition change. Unfortunately, manual construction of differential functional summary from a dE-MAP network that summarizes all pertinent functional responses is time-consuming, laborious and error-prone, impeding large-scale analysis on it. To this end, we propose DiffNet, a novel data-driven algorithm that leverages Gene Ontology (go) annotations to automatically summarize a dE-MAP network to obtain a high-level map of functional responses due to condition change. We tested DiffNet on the dynamic interaction networks following MMS treatment and demonstrated the superiority of our approach in generating differential functional summaries compared to state-of-the-art graph clustering methods. We studied the effects of parameters in DiffNet in controlling the quality of the summary. We also performed a case study that illustrates its utility.

Evolutionarily conserved genetic interactions with budding and fission yeast MutS identify orthologous relationships in mismatch repair-deficient cancer cells

Background

The evolutionarily conserved DNA mismatch repair (MMR) system corrects base-substitution and insertion-deletion mutations generated during erroneous replication. The mutation or inactivation of many MMR factors strongly predisposes to cancer, where the resulting tumors often display resistance to standard chemotherapeutics. A new direction to develop targeted therapies is the harnessing of synthetic genetic interactions, where the simultaneous loss of two otherwise non-essential factors leads to reduced cell fitness or death. High-throughput screening in human cells to directly identify such interactors for disease-relevant genes is now widespread, but often requires extensive case-by-case optimization. Here we asked if conserved genetic interactors (CGIs) with MMR genes from two evolutionary distant yeast species (Saccharomyces cerevisiae and Schizosaccharomyzes pombe) can predict orthologous genetic relationships in higher eukaryotes.

Methods

High-throughput screening was used to identify genetic interaction profiles for the MutSα and MutSβ heterodimer subunits (msh2Δ, msh3Δ, msh6Δ) of fission yeast. Selected negative interactors with MutSβ (msh2Δ/msh3Δ) were directly analyzed in budding yeast, and the CGI with SUMO-protease Ulp2 further examined after RNA interference/drug treatment in MSH2-deficient and -proficient human cells.

Results

This study identified distinct genetic profiles for MutSα and MutSβ, and supports a role for the latter in recombinatorial DNA repair. Approximately 28% of orthologous genetic interactions with msh2Δ/msh3Δ are conserved in both yeasts, a degree consistent with global trends across these species. Further, the CGI between budding/fission yeast msh2 and SUMO-protease Ulp2 is maintained in human cells (MSH2/SENP6), and enhanced by Olaparib, a PARP inhibitor that induces the accumulation of single-strand DNA breaks. This identifies SENP6 as a promising new target for the treatment of MMR-deficient cancers.

Conclusion

Our findings demonstrate the utility of employing evolutionary distance in tractable lower eukaryotes to predict orthologous genetic relationships in higher eukaryotes. Moreover, we provide novel insights into the genome maintenance functions of a critical DNA repair complex and propose a promising targeted treatment for MMR deficient tumors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-014-0068-4) contains supplementary material, which is available to authorized users.

Measuring genetic interactions in human cells by RNAi and imaging

Observation of how genetic interactions modulate phenotypes is a powerful method for dissecting their underlying molecular and functional networks. Whereas in model organisms genetic interaction studies are well established, systematic analysis of genetic interactions in human cells has remained challenging. Here we provide a detailed protocol for large-scale mapping of genetic interactions in human cells using a high-throughput phenotyping approach. Pairwise gene product depletion is induced by siRNA-mediated knockdown, and the resulting phenotypes are quantified by automated imaging and computational analysis to provide the basis for detecting genetic interactions between all pairs of genes tested. The whole workflow, depending on the size of the experiment, takes 3 or more weeks to complete.

esyN: network building, sharing and publishing

The construction and analysis of networks is increasingly widespread in biological research. We have developed esyN ("easy networks") as a free and open source tool to facilitate the exchange of biological network models between researchers. esyN acts as a searchable database of user-created networks from any field. We have developed a simple companion web tool that enables users to view and edit networks using data from publicly available databases. Both normal interaction networks (graphs) and Petri nets can be created. In addition to its basic tools, esyN contains a number of logical templates that can be used to create models more easily. The ability to use previously published models as building blocks makes esyN a powerful tool for the construction of models and network graphs. Users are able to save their own projects online and share them either publicly or with a list of collaborators. The latter can be given the ability to edit the network themselves, allowing online collaboration on network construction. esyN is designed to facilitate unrestricted exchange of this increasingly important type of biological information. Ultimately, the aim of esyN is to bring the advantages of Open Source software development to the construction of biological networks.

SANTA: quantifying the functional content of molecular networks

Linking networks of molecular interactions to cellular functions and phenotypes is a key goal in systems biology. Here, we adapt concepts of spatial statistics to assess the functional content of molecular networks. Based on the guilt-by-association principle, our approach (called SANTA) quantifies the strength of association between a gene set and a network, and functionally annotates molecular networks like other enrichment methods annotate lists of genes. As a general association measure, SANTA can (i) functionally annotate experimentally derived networks using a collection of curated gene sets and (ii) annotate experimentally derived gene sets using a collection of curated networks, as well as (iii) prioritize genes for follow-up analyses. We exemplify the efficacy of SANTA in several case studies using the S. cerevisiae genetic interaction network and genome-wide RNAi screens in cancer cell lines. Our theory, simulations, and applications show that SANTA provides a principled statistical way to quantify the association between molecular networks and cellular functions and phenotypes. SANTA is available from http://bioconductor.org/packages/release/bioc/html/SANTA.html.

The genomic landscape of compensatory evolution

The Genomic Landscape of Compensatory Evolution Laboratory selection experiment explains how organisms compensate for the loss of genes during evolution, and reveals the deleterious side-effects of this process when adapting to novel environments.

Adaptive evolution is generally assumed to progress through the accumulation of beneficial mutations. However, as deleterious mutations are common in natural populations, they generate a strong selection pressure to mitigate their detrimental effects through compensatory genetic changes. This process can potentially influence directions of adaptive evolution by enabling evolutionary routes that are otherwise inaccessible. Therefore, the extent to which compensatory mutations shape genomic evolution is of central importance. Here, we studied the capacity of the baker's yeast genome to compensate the complete loss of genes during evolution, and explored the long-term consequences of this process. We initiated laboratory evolutionary experiments with over 180 haploid baker's yeast genotypes, all of which initially displayed slow growth owing to the deletion of a single gene. Compensatory evolution following gene loss was rapid and pervasive: 68% of the genotypes reached near wild-type fitness through accumulation of adaptive mutations elsewhere in the genome. As compensatory mutations have associated fitness costs, genotypes with especially low fitnesses were more likely to be subjects of compensatory evolution. Genomic analysis revealed that as compensatory mutations were generally specific to the functional defect incurred, convergent evolution at the molecular level was extremely rare. Moreover, the majority of the gene expression changes due to gene deletion remained unrestored. Accordingly, compensatory evolution promoted genomic divergence of parallel evolving populations. However, these different evolutionary outcomes are not phenotypically equivalent, as they generated diverse growth phenotypes across environments. Taken together, these results indicate that gene loss initiates adaptive genomic changes that rapidly restores fitness, but this process has substantial pleiotropic effects on cellular physiology and evolvability upon environmental change. Our work also implies that gene content variation across species could be partly due to the action of compensatory evolution rather than the passive loss of genes.

While core cellular processes are generally conserved during evolution, the constituent genes differ somewhat between related species with similar lifestyles. Why should this be so? In this work, we propose that gene loss may initially be deleterious, but organisms can recover fitness by the accumulation of compensatory mutations elsewhere in the genome. To investigate this process in the laboratory, we investigated 180 haploid yeast strains, each of which initially displayed slow growth owing to the deletion of a single gene. Laboratory evolutionary experiments revealed that defects in a broad range of molecular processes can readily be compensated during evolution. Genomic analyses and functional assays demonstrated that compensatory evolution generates hidden genetic and physiological variation across parallel evolving lines, which can be revealed when the environment changes. Strikingly, despite nearly full recovery of fitness, the wild-type genomic expression pattern is generally not restored. Based on these results, we argue that genomes undergo major changes not simply to adapt to external conditions but also to compensate for previously accumulated deleterious mutations.

Prediction of dynamical drug sensitivity and resistance by module network rewiring-analysis based on transcriptional profiling

Revealing functional reorganization or module rewiring between modules at network levels during drug treatment is important to systematically understand therapies and drug responses. The present article proposed a novel model of module network rewiring to characterize functional reorganization of a complex biological system, and described a new framework named as module network rewiring-analysis (MNR) for systematically studying dynamical drug sensitivity and resistance during drug treatment. MNR was used to investigate functional reorganization or rewiring on the module network, rather than molecular network or individual molecules. Our experiments on expression data of patients with Hepatitis C virus infection receiving Interferon therapy demonstrated that consistent module genes derived by MNR could be directly used to reveal new genotypes relevant to drug sensitivity, unlike the other differential analyses of gene expressions. Our results showed that functional connections and reconnections among consistent modules bridged by biological paths were necessary for achieving effective responses of a drug. The hierarchical structures of the temporal module network can be considered as spatio-temporal biomarkers to monitor the efficacy, efficiency, toxicity, and resistance of the therapy. Our study indicates that MNR is a useful tool to identify module biomarkers and further predict dynamical drug sensitivity and resistance, characterize complex dynamic processes for therapy response, and provide biologically systematic clues for pharmacogenomic applications.

Modeling dynamic functional relationship networks and application to ex vivo human erythroid differentiation

MOTIVATION

Functional relationship networks, which summarize the probability of co-functionality between any two genes in the genome, could complement the reductionist focus of modern biology for understanding diverse biological processes in an organism. One major limitation of the current networks is that they are static, while one might expect functional relationships to consistently reprogram during the differentiation of a cell lineage. To address this potential limitation, we developed a novel algorithm that leverages both differentiation stage-specific expression data and large-scale heterogeneous functional genomic data to model such dynamic changes. We then applied this algorithm to the time-course RNA-Seq data we collected for ex vivo human erythroid cell differentiation.

RESULTS

Through computational cross-validation and literature validation, we show that the resulting networks correctly predict the (de)-activated functional connections between genes during erythropoiesis. We identified known critical genes, such as HBD and GATA1, and functional connections during erythropoiesis using these dynamic networks, while the traditional static network was not able to provide such information. Furthermore, by comparing the static and the dynamic networks, we identified novel genes (such as OSBP2 and PDZK1IP1) that are potential drivers of erythroid cell differentiation. This novel method of modeling dynamic networks is applicable to other differentiation processes where time-course genome-scale expression data are available, and should assist in generating greater understanding of the functional dynamics at play across the genome during development.

AVAILABILITY AND IMPLEMENTATION

The network described in this article is available at http://guanlab.ccmb.med.umich.edu/stageSpecificNetwork.

iRegulon: from a gene list to a gene regulatory network using large motif and track collections

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

Gene regulatory networks control developmental, homeostatic, and disease processes by governing precise levels and spatio-temporal patterns of gene expression. Determining their topology can provide mechanistic insight into these processes. Gene regulatory networks consist of interactions between transcription factors and their direct target genes. Each regulatory interaction represents the binding of the transcription factor to a specific DNA binding site near its target gene. Here we present a computational method, called iRegulon, to identify master regulators and direct target genes in a human gene signature, i.e. a set of co-expressed genes. iRegulon relies on the analysis of the regulatory sequences around each gene in the gene set to detect enriched TF motifs or ChIP-seq peaks, using databases of nearly 10.000 TF motifs and 1000 ChIP-seq data sets or “tracks”. Next, it associates enriched motifs and tracks with candidate transcription factors and determines the optimal subset of direct target genes. We validate iRegulon on ENCODE data, and use it in combination with RNA-seq and ChIP-seq data to map a p53 downstream network with new predicted co-factors and targets. iRegulon is available as a Cytoscape plugin, supporting human, mouse, and Drosophila genes, and provides access to hundreds of cancer-related TF-target subnetworks or “regulons”.

Functional genomics platform for pooled screening and generation of mammalian genetic interaction maps

Systematic genetic interaction maps in microorganisms are powerful tools for identifying functional relationships between genes and for defining the function of uncharacterized genes. We have recently implemented this strategy in mammalian cells as a two-stage approach. First, genes of interest are robustly identified in a pooled genome-wide screen using complex shRNA libraries. Second, phenotypes for all pairwise combinations of 'hit' genes are measured in a double-shRNA screen and used to construct a genetic interaction map. Our protocol allows for rapid pooled screening under various conditions without a requirement for robotics, in contrast to arrayed approaches. Each round of screening can be implemented in ∼2 weeks, with additional time for analysis and generation of reagents. We discuss considerations for screen design, and we present complete experimental procedures, as well as a full computational analysis suite for the identification of hits in pooled screens and generation of genetic interaction maps. Although the protocol outlined here was developed for our original shRNA-based approach, it can be applied more generally, including to CRISPR-based approaches.

Endocrinology: advances through omics and related technologies

The rapid development of new omics technologies to measure changes at genetic, transcriptomic, proteomic, and metabolomics levels together with the evolution of methods to analyze and integrate the data at a systems level are revolutionizing the study of biological processes. Here we discuss how new approaches using omics technologies have expanded our knowledge especially in nontraditional models. Our increasing knowledge of these interactions and evolutionary pathway conservation facilitates the use of nontraditional species, both invertebrate and vertebrate, as new model species for biological and endocrinology research. The increasing availability of technology to create organisms overexpressing key genes in endocrine function allows manipulation of complex regulatory networks such as growth hormone (GH) in transgenic fish where disregulation of GH production to produce larger fish has also permitted exploration of the role that GH plays in testis development, suggesting that it does so through interactions with insulin-like growth factors. The availability of omics tools to monitor changes at nearly any level in any organism, manipulate gene expression and behavior, and integrate data across biological levels, provides novel opportunities to explore endocrine function across many species and understand the complex roles that key genes play in different aspects of the endocrine function.

An integrative network algorithm identifies age-associated differential methylation interactome hotspots targeting stem-cell differentiation pathways

Epigenetic changes have been associated with ageing and cancer. Identifying and interpreting epigenetic changes associated with such phenotypes may benefit from integration with protein interactome models. We here develop and validate a novel integrative epigenome-interactome approach to identify differential methylation interactome hotspots associated with a phenotype of interest. We apply the algorithm to cancer and ageing, demonstrating the existence of hotspots associated with these phenotypes. Importantly, we discover tissue independent age-associated hotspots targeting stem-cell differentiation pathways, which we validate in independent DNA methylation data sets, encompassing over 1000 samples from different tissue types. We further show that these pathways would not have been discovered had we used a non-network based approach and that the use of the protein interaction network improves the overall robustness of the inference procedure. The proposed algorithm will be useful to any study seeking to identify interactome hotspots associated with common phenotypes.

A role for nuclear envelope-bridging complexes in homology-directed repair

Persistent double-strand DNA breaks (DSBs) are recruited to the nuclear periphery, where they induce formation of associated nuclear envelope–spanning LINC complexes made up of the SUN protein Sad1 and the KASH protein Kms1. The LINC complex couples DSBs within the nucleus to cytoplasmic microtubules, which alters DSB repair pathway choice.

Unless efficiently and faithfully repaired, DNA double-strand breaks (DSBs) cause genome instability. We implicate a Schizosaccharomyces pombe nuclear envelope–spanning linker of nucleoskeleton and cytoskeleton (LINC) complex, composed of the Sad1/Unc84 protein Sad1 and Klarsicht/Anc1/SYNE1 homology protein Kms1, in the repair of DSBs. An induced DSB associates with Sad1 and Kms1 in S/G2 phases of the cell cycle, connecting the DSB to cytoplasmic microtubules. DSB resection to generate single-stranded DNA and the ATR kinase drive the formation of Sad1 foci in response to DNA damage. Depolymerization of microtubules or loss of Kms1 leads to an increase in the number and size of DSB-induced Sad1 foci. Further, Kms1 and the cytoplasmic microtubule regulator Mto1 promote the repair of an induced DSB by gene conversion, a type of homology-directed repair. kms1 genetically interacts with a number of genes involved in homology-directed repair; these same gene products appear to attenuate the formation or promote resolution of DSB-induced Sad1 foci. We suggest that the connection of DSBs with the cytoskeleton through the LINC complex may serve as an input to repair mechanism choice and efficiency.

An RNA polymerase II-coupled function for histone H3K36 methylation in checkpoint activation and DSB repair

Histone modifications are major determinants of DNA double-strand break (DSB) response and repair. Here we elucidate a DSB repair function for transcription-coupled Set2 methylation at H3 lysine 36 (H3K36me). Cells devoid of Set2/H3K36me are hypersensitive to DNA-damaging agents and site-specific DSBs, fail to properly activate the DNA-damage checkpoint, and show genetic interactions with DSB-sensing and repair machinery. Set2/H3K36me3 is enriched at DSBs, and loss of Set2 results in altered chromatin architecture and inappropriate resection during G1 near break sites. Surprisingly, Set2 and RNA polymerase II are programmed for destruction after DSBs in a temporal manner – resulting in H3K36me3 to H3K36me2 transition that may be linked to DSB repair. Finally, we show a requirement of Set2 in DSB repair in transcription units – thus underscoring the importance of transcription-dependent H3K36me in DSB repair.

H2A.Z-dependent regulation of cohesin dynamics on chromosome arms

Structural maintenance of chromosomes (SMC) complexes and DNA topoisomerases are major determinants of chromosome structure and dynamics. The cohesin complex embraces sister chromatids throughout interphase, but during mitosis most cohesin is stripped from chromosome arms by early prophase, while the remaining cohesin at kinetochores is cleaved at anaphase. This two-step removal of cohesin is required for sister chromatids to separate. The cohesin-related Smc5/6 complex has been studied mostly as a determinant of DNA repair via homologous recombination. However, chromosome segregation fails in Smc5/6 null mutants or cells treated with small interfering RNAs. This also occurs in Smc5/6 hypomorphs in the fission yeast Schizosaccharomyces pombe following genotoxic and replication stress, or topoisomerase II dysfunction, and these mitotic defects are due to the postanaphase retention of cohesin on chromosome arms. Here we show that mitotic and repair roles for Smc5/6 are genetically separable in S. pombe. Further, we identified the histone variant H2A.Z as a critical factor to modulate cohesin dynamics, and cells lacking H2A.Z suppress the mitotic defects conferred by Smc5/6 dysfunction. Together, H2A.Z and the SMC complexes ensure genome integrity through accurate chromosome segregation.

Eaf5/7/3 form a functionally independent NuA4 submodule linked to RNA polymerase II-coupled nucleosome recycling

The NuA4 histone acetyltransferase complex is required for gene regulation, cell cycle progression, and DNA repair. Dissection of the 13-subunit complex reveals that the Eaf7 subunit bridges Eaf5 with Eaf3, a H3K36me3-binding chromodomain protein, and this Eaf5/7/3 trimer is anchored to NuA4 through Eaf5. This trimeric subcomplex represents a functional module, and a large portion exists in a native form outside the NuA4 complex. Gene-specific and genome-wide location analyses indicate that Eaf5/7/3 correlates with transcription activity and is enriched over the coding region. In agreement with a role in transcription elongation, the Eaf5/7/3 trimer interacts with phosphorylated RNA polymerase II and helps its progression. Loss of Eaf5/7/3 partially suppresses intragenic cryptic transcription arising in set2 mutants, supporting a role in nucleosome destabilization. On the other hand, loss of the trimer leads to an increase of replication-independent histone exchange over the coding region of transcribed genes. Taken together, these results lead to a model where Eaf5/7/3 associates with elongating polymerase to promote the disruption of nucleosomes in its path, but also their refolding in its wake.

Parallel profiling of the transcriptome, cistrome, and epigenome in the cellular response to ionizing radiation

The DNA damage response (DDR) is a vast signaling network that is robustly activated by DNA double-strand breaks, the critical lesion induced by ionizing radiation (IR). Although much of this response operates at the protein level, a critical component of the network sustains many DDR branches by modulating the cellular transcriptome. Using deep sequencing, we delineated three layers in the transcriptional response to IR in human breast cancer cells: changes in the expression of genes encoding proteins or long noncoding RNAs, alterations in genomic binding by key transcription factors, and dynamics of epigenetic markers of active promoters and enhancers. We identified protein-coding and previously unidentified noncoding genes that were responsive to IR, and demonstrated that IR-induced transcriptional dynamics was mediated largely by the transcription factors p53 and nuclear factor κB (NF-κB) and was primarily dependent on the kinase ataxia-telangiectasia mutated (ATM). The resultant data set provides a rich resource for understanding a basic, underlying component of a critical cellular stress response.

Direct estimation of differential networks

It is often of interest to understand how the structure of a genetic network differs between two conditions. In this paper, each condition-specific network is modeled using the precision matrix of a multivariate normal random vector, and a method is proposed to directly estimate the difference of the precision matrices. In contrast to other approaches, such as separate or joint estimation of the individual matrices, direct estimation does not require those matrices to be sparse, and thus can allow the individual networks to contain hub nodes. Under the assumption that the true differential network is sparse, the direct estimator is shown to be consistent in support recovery and estimation. It is also shown to outperform existing methods in simulations, and its properties are illustrated on gene expression data from late-stage ovarian cancer patients.