Rewiring of genetic networks in response to DNA damage

Oxidative stress-induced YAP1 expression is regulated by NCE102, CDA2, and BCS1

Maintaining cellular homeostasis in the face of stress conditions is vital for the overall well-being of an organism. Reactive oxygen species (ROS) are among the most potent cellular stressors and can disrupt the internal redox balance, giving rise to oxidative stress. Elevated levels of ROS can severely affect biomolecules and have been associated with a range of pathophysiological conditions. In response to oxidative stress, yeast activator protein-1 (Yap1p) undergoes post-translation modification that results in its nuclear accumulation. YAP1 has a key role in oxidative detoxification by promoting transcription of numerous antioxidant genes. In this study, we identified previously undescribed functions for NCE102, CDA2, and BCS1 in YAP1 expression in response to oxidative stress induced by hydrogen peroxide (H2O2). Deletion mutant strains for these candidates demonstrated increased sensitivity to H2O2. Our follow-up investigation linked the activity of these genes to YAP1 expression at the level of translation. Under oxidative stress, global cap-dependent translation is inhibited, prompting stress-responsive genes like YAP1 to employ alternative modes of translation. We provide evidence that NCE102, CDA2, and BCS1 contribute to cap-independent translation of YAP1 under oxidative stress.

Histone variant H2A.Z is needed for efficient transcription-coupled NER and genome integrity in UV challenged yeast cells

The genome of living cells is constantly challenged by DNA lesions that interfere with cellular processes such as transcription and replication. A manifold of mechanisms act in concert to ensure adequate DNA repair, gene expression, and genome stability. Bulky DNA lesions, such as those induced by UV light or the DNA-damaging agent 4-nitroquinoline oxide, act as transcriptional and replicational roadblocks and thus represent a major threat to cell metabolism. When located on the transcribed strand of active genes, these lesions are handled by transcription-coupled nucleotide excision repair (TC-NER), a yet incompletely understood NER sub-pathway. Here, using a genetic screen in the yeast Saccharomyces cerevisiae, we identified histone variant H2A.Z as an important component to safeguard transcription and DNA integrity following UV irradiation. In the absence of H2A.Z, repair by TC-NER is severely impaired and RNA polymerase II clearance reduced, leading to an increase in double-strand breaks. Thus, H2A.Z is needed for proficient TC-NER and plays a major role in the maintenance of genome stability upon UV irradiation.

A method for mining condition-specific co-expressed genes in Camellia sinensis based on k-means clustering

BACKGROUND

As one of the world's most important beverage crops, tea plants (Camellia sinensis) are renowned for their unique flavors and numerous beneficial secondary metabolites, attracting researchers to investigate the formation of tea quality. With the increasing availability of transcriptome data on tea plants in public databases, conducting large-scale co-expression analyses has become feasible to meet the demand for functional characterization of tea plant genes. However, as the multidimensional noise increases, larger-scale co-expression analyses are not always effective. Analyzing a subset of samples generated by effectively downsampling and reorganizing the global sample set often leads to more accurate results in co-expression analysis. Meanwhile, global-based co-expression analyses are more likely to overlook condition-specific gene interactions, which may be more important and worthy of exploration and research.

RESULTS

Here, we employed the k-means clustering method to organize and classify the global samples of tea plants, resulting in clustered samples. Metadata annotations were then performed on these clustered samples to determine the "conditions" represented by each cluster. Subsequently, we conducted gene co-expression network analysis (WGCNA) separately on the global samples and the clustered samples, resulting in global modules and cluster-specific modules. Comparative analyses of global modules and cluster-specific modules have demonstrated that cluster-specific modules exhibit higher accuracy in co-expression analysis. To measure the degree of condition specificity of genes within condition-specific clusters, we introduced the correlation difference value (CDV). By incorporating the CDV into co-expression analyses, we can assess the condition specificity of genes. This approach proved instrumental in identifying a series of high CDV transcription factor encoding genes upregulated during sustained cold treatment in Camellia sinensis leaves and buds, and pinpointing a pair of genes that participate in the antioxidant defense system of tea plants under sustained cold stress.

CONCLUSIONS

To summarize, downsampling and reorganizing the sample set improved the accuracy of co-expression analysis. Cluster-specific modules were more accurate in capturing condition-specific gene interactions. The introduction of CDV allowed for the assessment of condition specificity in gene co-expression analyses. Using this approach, we identified a series of high CDV transcription factor encoding genes related to sustained cold stress in Camellia sinensis. This study highlights the importance of considering condition specificity in co-expression analysis and provides insights into the regulation of the cold stress in Camellia sinensis.

Mouse models to explore the biological and organismic role of DNA polymerase beta

Gene knock-out (KO) mouse models for DNA polymerase beta (Polβ) revealed that loss of Polβ leads to neonatal lethality, highlighting the critical organismic role for this DNA polymerase. While biochemical analysis and gene KO cell lines have confirmed its biochemical role in base excision repair and in TET-mediated demethylation, more long-lived mouse models continue to be developed to further define its organismic role. The Polb-KO mouse was the first of the Cre-mediated tissue-specific KO mouse models. This technology was exploited to investigate roles for Polβ in V(D)J recombination (variable-diversity-joining rearrangement), DNA demethylation, gene complementation, SPO11-induced DNA double-strand break repair, germ cell genome stability, as well as neuronal differentiation, susceptibility to genotoxin-induced DNA damage, and cancer onset. The revolution in knock-in (KI) mouse models was made possible by CRISPR/cas9-mediated gene editing directly in C57BL/6 zygotes. This technology has helped identify phenotypes associated with germline or somatic mutants of Polβ. Such KI mouse models have helped uncover the importance of key Polβ active site residues or specific Polβ enzyme activities, such as the PolbY265C mouse that develops lupus symptoms. More recently, we have used this KI technology to mutate the Polb gene with two codon changes, yielding the PolbL301R/V303R mouse. In this KI mouse model, the expressed Polβ protein cannot bind to its obligate heterodimer partner, Xrcc1. Although the expressed mutant Polβ protein is proteolytically unstable and defective in recruitment to sites of DNA damage, the homozygous PolbL301R/V303R mouse is viable and fertile, yet small in stature. We expect that this and additional targeted mouse models under development are poised to reveal new biological and organismic roles for Polβ.

Gene set correlation enrichment analysis for interpreting and annotating gene expression profiles

Pathway analysis, including nontopology-based (non-TB) and topology-based (TB) methods, is widely used to interpret the biological phenomena underlying differences in expression data between two phenotypes. By considering dependencies and interactions between genes, TB methods usually perform better than non-TB methods in identifying pathways that include closely relevant or directly causative genes for a given phenotype. However, most TB methods may be limited by incomplete pathway data used as the reference network or by difficulties in selecting appropriate reference networks for different research topics. Here, we propose a gene set correlation enrichment analysis method, Gscore, based on an expression dataset-derived coexpression network to examine whether a differentially expressed gene (DEG) list (or each of its DEGs) is associated with a known gene set. Gscore is better able to identify target pathways in 89 human disease expression datasets than eight other state-of-the-art methods and offers insight into how disease-wide and pathway-wide associations reflect clinical outcomes. When applied to RNA-seq data from COVID-19-related cells and patient samples, Gscore provided a means for studying how DEGs are implicated in COVID-19-related pathways. In summary, Gscore offers a powerful analytical approach for annotating individual DEGs, DEG lists, and genome-wide expression profiles based on existing biological knowledge.

Genetic Interactions in Various Environmental Conditions in Caenorhabditis elegans

Although it is well known that epistasis plays an important role in many evolutionary processes (e.g., speciation, evolution of sex), our knowledge on the frequency and prevalent sign of epistatic interactions is mainly limited to unicellular organisms or cell cultures of multicellular organisms. This is even more pronounced in regard to how the environment can influence genetic interactions. To broaden our knowledge in that respect we studied gene-gene interactions in a whole multicellular organism, Caenorhabditis elegans. We screened over one thousand gene interactions, each one in standard laboratory conditions, and under three different stressors: heat shock, oxidative stress, and genotoxic stress. Depending on the condition, between 7% and 22% of gene pairs showed significant genetic interactions and an overall sign of epistasis changed depending on the condition. Sign epistasis was quite common, but reciprocal sign epistasis was extremally rare. One interaction was common to all conditions, whereas 78% of interactions were specific to only one environment. Although epistatic interactions are quite common, their impact on evolutionary processes will strongly depend on environmental factors.

Role of genetic architecture in phenotypic plasticity

Phenotypic plasticity, the ability of an organism to display different phenotypes across environments, is widespread in nature. Plasticity aids survival in novel environments. Herein, we review studies from yeast that allow us to start uncovering the genetic architecture of phenotypic plasticity. Genetic variants and their interactions impact the phenotype in different environments, and distinct environments modulate the impact of genetic variants and their interactions on the phenotype. Because of this, certain hidden genetic variation is expressed in specific genetic and environmental backgrounds. A better understanding of the genetic mechanisms of phenotypic plasticity will help to determine short- and long-term responses to selection and how wide variation in disease manifestation occurs in human populations.

Mapping the Genetic Interaction Network of PARP inhibitor Response

Genetic interactions have long informed our understanding of the coordinated proteins and pathways that respond to DNA damage in mammalian cells, but systematic interrogation of the genetic network underlying that system has yet to be achieved. Towards this goal, we measured 147,153 pairwise interactions among genes implicated in PARP inhibitor (PARPi) response. Evaluating genetic interactions at this scale, with and without exposure to PARPi, revealed hierarchical organization of the pathways and complexes that maintain genome stability during normal growth and defined changes that occur upon accumulation of DNA lesions due to cytotoxic doses of PARPi. We uncovered unexpected relationships among DNA repair genes, including context-specific buffering interactions between the minimally characterized AUNIP and BRCA1-A complex genes. Our work thus establishes a foundation for mapping differential genetic interactions in mammalian cells and provides a comprehensive resource for future studies of DNA repair and PARP inhibitors.

Environmental challenge rewires functional connections among human genes

A fundamental question in biology is how a limited number of genes combinatorially govern cellular responses to environmental changes. While the prevailing hypothesis is that relationships between genes, processes, and ontologies could be plastic to achieve this adaptability, quantitatively comparing human gene functional connections between specific environmental conditions at scale is very challenging. Therefore, it remains unclear whether and how human genetic interaction networks are rewired in response to changing environmental conditions. Here, we developed a framework for mapping context-specific genetic interactions, enabling us to measure the plasticity of human genetic architecture upon environmental challenge for ~250,000 interactions, using cell cycle interruption, genotoxic perturbation, and nutrient deprivation as archetypes. We discover large-scale rewiring of human gene relationships across conditions, highlighted by dramatic shifts in the functional connections of epigenetic regulators (TIP60), cell cycle regulators (PP2A), and glycolysis metabolism. Our study demonstrates that upon environmental perturbation, intra-complex genetic rewiring is rare while inter-complex rewiring is common, suggesting a modular and flexible evolutionary genetic strategy that allows a limited number of human genes to enable adaptation to a large number of environmental conditions.

MOBILE pipeline enables identification of context-specific networks and regulatory mechanisms

Robust identification of context-specific network features that control cellular phenotypes remains a challenge. We here introduce MOBILE (Multi-Omics Binary Integration via Lasso Ensembles) to nominate molecular features associated with cellular phenotypes and pathways. First, we use MOBILE to nominate mechanisms of interferon-γ (IFNγ) regulated PD-L1 expression. Our analyses suggest that IFNγ-controlled PD-L1 expression involves BST2, CLIC2, FAM83D, ACSL5, and HIST2H2AA3 genes, which were supported by prior literature. We also compare networks activated by related family members transforming growth factor-beta 1 (TGFβ1) and bone morphogenetic protein 2 (BMP2) and find that differences in ligand-induced changes in cell size and clustering properties are related to differences in laminin/collagen pathway activity. Finally, we demonstrate the broad applicability and adaptability of MOBILE by analyzing publicly available molecular datasets to investigate breast cancer subtype specific networks. Given the ever-growing availability of multi-omics datasets, we envision that MOBILE will be broadly useful for identification of context-specific molecular features and pathways.

A multi-scale map of protein assemblies in the DNA damage response

The DNA damage response (DDR) ensures error-free DNA replication and transcription and is disrupted in numerous diseases. An ongoing challenge is to determine the proteins orchestrating DDR and their organization into complexes, including constitutive interactions and those responding to genomic insult. Here, we use multi-conditional network analysis to systematically map DDR assemblies at multiple scales. Affinity purifications of 21 DDR proteins, with/without genotoxin exposure, are combined with multi-omics data to reveal a hierarchical organization of 605 proteins into 109 assemblies. The map captures canonical repair mechanisms and proposes new DDR-associated proteins extending to stress, transport, and chromatin functions. We find that protein assemblies closely align with genetic dependencies in processing specific genotoxins and that proteins in multiple assemblies typically act in multiple genotoxin responses. Follow-up by DDR functional readouts newly implicates 12 assembly members in double-strand-break repair. The DNA damage response assemblies map is available for interactive visualization and query (ccmi.org/ddram/).

Bayesian Machine Learning Enables Identification of Transcriptional Network Disruptions Associated with Drug-Resistant Prostate Cancer

Survival rates of patients with metastatic castration-resistant prostate cancer (mCRPC) are low due to lack of response or acquired resistance to available therapies, such as abiraterone (Abi). A better understanding of the underlying molecular mechanisms is needed to identify effective targets to overcome resistance. Given the complexity of the transcriptional dynamics in cells, differential gene expression analysis of bulk transcriptomics data cannot provide sufficient detailed insights into resistance mechanisms. Incorporating network structures could overcome this limitation to provide a global and functional perspective of Abi resistance in mCRPC. Here, we developed TraRe, a computational method using sparse Bayesian models to examine phenotypically driven transcriptional mechanistic differences at three distinct levels: transcriptional networks, specific regulons, and individual transcription factors (TF). TraRe was applied to transcriptomic data from 46 patients with mCRPC with Abi-response clinical data and uncovered abrogated immune response transcriptional modules that showed strong differential regulation in Abi-responsive compared with Abi-resistant patients. These modules were replicated in an independent mCRPC study. Furthermore, key rewiring predictions and their associated TFs were experimentally validated in two prostate cancer cell lines with different Abi-resistance features. Among them, ELK3, MXD1, and MYB played a differential role in cell survival in Abi-sensitive and Abi-resistant cells. Moreover, ELK3 regulated cell migration capacity, which could have a direct impact on mCRPC. Collectively, these findings shed light on the underlying transcriptional mechanisms driving Abi response, demonstrating that TraRe is a promising tool for generating novel hypotheses based on identified transcriptional network disruptions.

SIGNIFICANCE

The computational method TraRe built on Bayesian machine learning models for investigating transcriptional network structures shows that disruption of ELK3, MXD1, and MYB signaling cascades impacts abiraterone resistance in prostate cancer.

Complexity and self-organization in the evolution of cell polarization

Cellular life exhibits order and complexity, which typically increase over the course of evolution. Cell polarization is a well-studied example of an ordering process that breaks the internal symmetry of a cell by establishing a preferential axis. Like many cellular processes, polarization is driven by self-organization, meaning that the macroscopic pattern emerges as a consequence of microscopic molecular interactions at the biophysical level. However, the role of self-organization in the evolution of complex protein networks remains obscure. In this Review, we provide an overview of the evolution of polarization as a self-organizing process, focusing on the model species Saccharomyces cerevisiae and its fungal relatives. Moreover, we use this model system to discuss how self-organization might relate to evolutionary change, offering a shift in perspective on evolution at the microscopic scale.

Systematic comparison of differential expression networks in MTB mono-, HIV mono- and MTB/HIV co-infections for drug repurposing

The synergy between human immunodeficiency virus (HIV) and Mycobacterium tuberculosis (MTB) could accelerate the deterioration of immunological functions. Previous studies have explored the pathogenic mechanisms of HIV mono-infection (HMI), MTB mono-infection (MMI) and MTB/HIV co-infection (MHCI), but their similarities and specificities remain to be profoundly investigated. We thus designed a computational framework named IDEN to identify gene pairs related to these states, which were then compared from different perspectives. MMI-related genes showed the highest enrichment level on a greater number of chromosomes. Genes shared by more states tended to be more evolutionarily conserved, posttranslationally modified and topologically important. At the expression level, HMI-specific gene pairs yielded higher correlations, while the overlapping pairs involved in MHCI had significantly lower correlations. The correlation changes of common gene pairs showed that MHCI shared more similarities with MMI. Moreover, MMI- and MHCI-related genes were enriched in more identical pathways and biological processes, further illustrating that MTB may play a dominant role in co-infection. Hub genes specific to each state could promote pathogen infections, while those shared by two states could enhance immune responses. Finally, we improved the network proximity measure for drug repurposing by considering the importance of gene pairs, and approximately ten drug candidates were identified for each disease state.

Applying differential network analysis to longitudinal gene expression in response to perturbations

Differential Network (DN) analysis is a method that has long been used to interpret changes in gene expression data and provide biological insights. The method identifies the rewiring of gene networks in response to external perturbations. Our study applies the DN method to the analysis of RNA-sequencing (RNA-seq) time series datasets. We focus on expression changes: (i) in saliva of a human subject after pneumococcal vaccination (PPSV23) and (ii) in primary B cells treated ex vivo with a monoclonal antibody drug (Rituximab). The DN method enabled us to identify the activation of biological pathways consistent with the mechanisms of action of the PPSV23 vaccine and target pathways of Rituximab. The community detection algorithm on the DN revealed clusters of genes characterized by collective temporal behavior. All saliva and some B cell DN communities showed characteristic time signatures, outlining a chronological order in pathway activation in response to the perturbation. Moreover, we identified early and delayed responses within network modules in the saliva dataset and three temporal patterns in the B cell data.

Gene Differential Co-Expression Networks Based on RNA-Seq: Construction and Its Applications

Gene co-expression network (GCN) becomes an increasingly important tool in omics data analysis. A great challenge for GCN construction is that the sample size is far lower than the number of genes. Traditional methods rely on considerable samples. Moreover, association signals are likely weak, nonlinear and stochastic, which are difficult to be identified among thousands of candidates. In this paper, the gray correlation coefficient (GCC) is introduced, and a novel method to construct gene differential co-expression networks (GDCNs) is proposed. Based on the GDCNs, three measures are proposed to explore informative genes. The proposed method can make full use of the information provided by a handful of samples and overcome the shortages of GCNs, which can evaluate the changes of co-expression relationships that are possibly triggered by treatments. Based on RNA-seq data of Brassica napus, GDCNs under multiple experimental conditions are constructed and investigated. It is found that the GCC-based method is very robust to data processing. The GDCNs facilitate the inference of gene functions and the identification of informative genes that are responsible for stress responsiveness. The GDCN-based approaches integrate the 'guilt by association' and the 'guilt by rewiring' rules, which provide alternative tools for omics data analysis.

Exploring conditional gene essentiality through systems genetics approaches in yeast

An essential gene encodes for a cellular function that is required for viability. Although viability is a straightforward phenotype to analyze in yeast, defining a gene as essential is not always trivial. Gene essentiality has generally been studied in specific laboratory strains and under standard growth conditions, however, essentiality can vary across species, strains, and environments. Recent systematic studies of gene essentiality revealed that two sets of essential genes exist: core essential genes that are always required for viability and conditional essential genes that vary in essentiality in different genetic and environmental contexts. Here, we review recent advances made in the systematic analysis of gene essentiality in yeast and discuss the properties that distinguish core from context-dependent essential genes.

Dynamic rewiring of biological activity across genotype and lineage revealed by context-dependent functional interactions

BACKGROUND

Coessentiality networks derived from CRISPR screens in cell lines provide a powerful framework for identifying functional modules in the cell and for inferring the roles of uncharacterized genes. However, these networks integrate signal across all underlying data and can mask strong interactions that occur in only a subset of the cell lines analyzed.

RESULTS

Here, we decipher dynamic functional interactions by identifying significant cellular contexts, primarily by oncogenic mutation, lineage, and tumor type, and discovering coessentiality relationships that depend on these contexts. We recapitulate well-known gene-context interactions such as oncogene-mutation, paralog buffering, and tissue-specific essential genes, show how mutation rewires known signal transduction pathways, including RAS/RAF and IGF1R-PIK3CA, and illustrate the implications for drug targeting. We further demonstrate how context-dependent functional interactions can elucidate lineage-specific gene function, as illustrated by the maturation of proreceptors IGF1R and MET by proteases FURIN and CPD.

CONCLUSIONS

This approach advances our understanding of context-dependent interactions and how they can be gleaned from these data. We provide an online resource to explore these context-dependent interactions at diffnet.hart-lab.org.

Application of modular response analysis to medium- to large-size biological systems

The development of high-throughput genomic technologies associated with recent genetic perturbation techniques such as short hairpin RNA (shRNA), gene trapping, or gene editing (CRISPR/Cas9) has made it possible to obtain large perturbation data sets. These data sets are invaluable sources of information regarding the function of genes, and they offer unique opportunities to reverse engineer gene regulatory networks in specific cell types. Modular response analysis (MRA) is a well-accepted mathematical modeling method that is precisely aimed at such network inference tasks, but its use has been limited to rather small biological systems so far. In this study, we show that MRA can be employed on large systems with almost 1,000 network components. In particular, we show that MRA performance surpasses general-purpose mutual information-based algorithms. Part of these competitive results was obtained by the application of a novel heuristic that pruned MRA-inferred interactions a posteriori. We also exploited a block structure in MRA linear algebra to parallelize large system resolutions.

RNAi Screening Uncovers a Synthetic Sick Interaction between CtIP and the BARD1 Tumor Suppressor

Human CtIP is best known for its role in DNA end resection to initiate DNA double-strand break repair by homologous recombination. Recently, CtIP has also been shown to protect reversed replication forks from nucleolytic degradation upon DNA replication stress. However, still little is known about the DNA damage response (DDR) networks that preserve genome integrity and sustain cell survival in the context of CtIP insufficiency. Here, to reveal such potential buffering relationships, we screened a DDR siRNA library in CtIP-deficient cells to identify candidate genes that induce synthetic sickness/lethality (SSL). Our analyses unveil a negative genetic interaction between CtIP and BARD1, the heterodimeric binding partner of BRCA1. We found that simultaneous disruption of CtIP and BARD1 triggers enhanced apoptosis due to persistent replication stress-induced DNA lesions giving rise to chromosomal abnormalities. Moreover, we observed that the genetic interaction between CtIP and BARD1 occurs independently of the BRCA1-BARD1 complex formation and might be, therefore, therapeutical relevant for the treatment of BRCA-defective tumors.

Identifying biomarkers for breast cancer by gene regulatory network rewiring

Background

Mining gene regulatory network (GRN) is an important avenue for addressing cancer mechanism. Mutations in cancer genome perturb GRN and cause a rewiring in an orchestrated network. Hence, the exploration of gene regulatory network rewiring is significant to discover potential biomarkers and indicators for discriminating cancer phenotypes.

Results

Here, we propose a new bioinformatics method of identifying biomarkers based on network rewiring in different states. It firstly reconstructs GRN in different phenotypic conditions from gene expression data with a priori background network. We employ the algorithm based on path consistency algorithm and conditional mutual information to delete false-positive regulatory interactions between independent nodes/genes or not closely related gene pairs. And then a differential gene regulatory network (D-GRN) is constructed from the rewiring parts in the two phenotype-specific GRNs. Community detection technique is then applied for D-GRN to detect functional modules. Finally, we apply logistic regression classifier with recursive feature elimination to select biomarker genes in each module individually. The extracted feature genes result in a gene set of biomarkers with impressing ability to distinguish normal samples from controls. We verify the identified biomarkers in external independent validation datasets. For a proof-of-concept study, we apply the framework to identify diagnostic biomarkers of breast cancer. The identified biomarkers obtain a maximum AUC of 0.985 in the internal sample classification experiments. And these biomarkers achieve a maximum AUC of 0.989 in the external validations.

Conclusion

In conclusion, network rewiring reveals significant differences between different phenotypes, which indicating cancer dysfunctional mechanisms. With the development of sequencing technology, the amount and quality of gene expression data become available. Condition-specific gene regulatory networks that are close to the real regulations in different states will be established. Revealing the network rewiring will greatly benefit the discovery of biomarkers or signatures for phenotypes. D-GRN is a general method to meet this demand of deciphering the high-throughput data for biomarker discovery. It is also easy to be extended for identifying biomarkers of other complex diseases beyond breast cancer.

Identifying Gene Network Rewiring Based on Partial Correlation

It is an important task to learn how gene regulatory networks change under different conditions. Several Gaussian graphical model-based methods have been proposed to deal with this task by inferring differential networks from gene expression data. However, most existing methods define the differential networks as the difference of precision matrices, which may include false differential edges caused by the change of conditional variances. In addition, prior information about the condition-specific networks and the differential networks can be obtained from other domains. It is useful to incorporate prior information into differential network analysis. In this study, we propose a new differential network analysis method to address the above challenges. Instead of using the precision matrices, we define the differential networks as the difference of partial correlations, which can exclude the spurious differential edges due to the variants of conditional variances. Furthermore, prior information from multiple hypothesis testing is incorporated using a weighted fused penalty. Simulation studies show that our method outperforms the competing methods. We also apply our method to identify the differential network between luminal A and basal-like subtypes of breast cancers and the differential network between acute myeloid leukemia tumors and normal samples. The hub genes in the differential networks identified by our method carry out important biological functions.

Gene utility recapitulates chromosomal aberrancies in advanced stage neuroblastoma

Neuroblastoma (NB) is the most common extracranial solid tumor in children. Although only a few recurrent somatic mutations have been identified, chromosomal abnormalities, including the loss of heterozygosity (LOH) at the chromosome 1p and gains of chromosome 17q, are often seen in the high-risk cases. The biological basis and evolutionary forces that drive such genetic abnormalities remain enigmatic. Here, we conceptualize the Gene Utility Model (GUM) that seeks to identify genes driving biological signaling via their collective gene utilities and apply it to understand the impact of those differentially utilized genes on constraining the evolution of NB karyotypes. By employing a computational process-guided flow algorithm to model gene utility in protein–protein networks that built based on transcriptomic data, we conducted several pairwise comparative analyses to uncover genes with differential utilities in stage 4 NBs with distinct classification. We then constructed a utility karyotype by mapping these differentially utilized genes to their respective chromosomal loci. Intriguingly, hotspots of the utility karyotype, to certain extent, can consistently recapitulate the major chromosomal abnormalities of NBs and also provides clues to yet identified predisposition sites. Hence, our study not only provides a new look, from a gene utility perspective, into the known chromosomal abnormalities detected by integrative genomic sequencing efforts, but also offers new insights into the etiology of NB and provides a framework to facilitate the identification of novel therapeutic targets for this devastating childhood cancer.

The CWI Pathway: A Versatile Toolbox to Arrest Cell-Cycle Progression

Cell-signaling pathways are essential for cells to respond and adapt to changes in their environmental conditions. The cell-wall integrity (CWI) pathway of Saccharomyces cerevisiae is activated by environmental stresses, compounds, and morphogenetic processes that compromise the cell wall, orchestrating the appropriate cellular response to cope with these adverse conditions. During cell-cycle progression, the CWI pathway is activated in periods of polarized growth, such as budding or cytokinesis, regulating cell-wall biosynthesis and the actin cytoskeleton. Importantly, accumulated evidence has indicated a reciprocal regulation of the cell-cycle regulatory system by the CWI pathway. In this paper, we describe how the CWI pathway regulates the main cell-cycle transitions in response to cell-surface perturbance to delay cell-cycle progression. In particular, it affects the Start transcriptional program and the initiation of DNA replication at the G1/S transition, and entry and progression through mitosis. We also describe the involvement of the CWI pathway in the response to genotoxic stress and its connection with the DNA integrity checkpoint, the mechanism that ensures the correct transmission of genetic material and cell survival. Thus, the CWI pathway emerges as a master brake that stops cell-cycle progression when cells are coping with distinct unfavorable conditions.

Discovery of putative tumor suppressors from CRISPR screens reveals rewired lipid metabolism in acute myeloid leukemia cells

CRISPR knockout fitness screens in cancer cell lines reveal many genes whose loss of function causes cell death or loss of fitness or, more rarely, the opposite phenotype of faster proliferation. Here we demonstrate a systematic approach to identify these proliferation suppressors, which are highly enriched for tumor suppressor genes, and define a network of 145 such genes in 22 modules. One module contains several elements of the glycerolipid biosynthesis pathway and operates exclusively in a subset of acute myeloid leukemia cell lines. The proliferation suppressor activity of genes involved in the synthesis of saturated fatty acids, coupled with a more severe loss of fitness phenotype for genes in the desaturation pathway, suggests that these cells operate at the limit of their carrying capacity for saturated fatty acids, which we confirm biochemically. Overexpression of this module is associated with a survival advantage in juvenile leukemias, suggesting a clinically relevant subtype.

CRISPR-based knockout screens in cancer cells have suggested the existence of proliferation suppressor genes (PSG). Here, the authors develop an approach to systematically identify them, and reveal a PSG module involved in fatty acid synthesis and tumour suppression in acute myeloid leukemia cell lines.

Importance-Penalized Joint Graphical Lasso (IPJGL): differential network inference via GGMs

MOTIVATION

Differential network inference is a fundamental and challenging problem to reveal gene interactions and regulation relationships under different conditions. Many algorithms have been developed for this problem; however, they do not consider the differences between the importance of genes, which may not fit the real-world situation. Different genes have different mutation probabilities, and the vital genes associated with basic life activities have less fault tolerance to mutation. Equally treating all genes may bias the results of differential network inference. Thus, it is necessary to consider the importance of genes in the models of differential network inference.

RESULTS

Based on the Gaussian graphical model with adaptive gene importance regularization, we develop a novel Importance-Penalized Joint Graphical Lasso method (IPJGL) for differential network inference. The presented method is validated by the simulation experiments as well as the real datasets. Furthermore, to precisely evaluate the results of differential network inference, we propose a new metric named APC2 for the differential levels of gene pairs. We apply IPJGL to analyze the TCGA colorectal and breast cancer datasets and find some candidate cancer genes with significant survival analysis results, including SOST for colorectal cancer and RBBP8 for breast cancer. We also conduct further analysis based on the interactions in the Reactome database and confirm the utility of our method.

AVAILABILITY AND IMPLEMENTATION

R source code of Importance-Penalized Joint Graphical Lasso is freely available at https://github.com/Wu-Lab/IPJGL.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

CorDiffViz: an R package for visualizing multi-omics differential correlation networks

BACKGROUND

Differential correlation networks are increasingly used to delineate changes in interactions among biomolecules. They characterize differences between omics networks under two different conditions, and can be used to delineate mechanisms of disease initiation and progression.

RESULTS

We present a new R package, CorDiffViz, that facilitates the estimation and visualization of differential correlation networks using multiple correlation measures and inference methods. The software is implemented in R, HTML and Javascript, and is available at https://github.com/sqyu/CorDiffViz . Visualization has been tested for the Chrome and Firefox web browsers. A demo is available at https://diffcornet.github.io/CorDiffViz/demo.html .

CONCLUSIONS

Our software offers considerable flexibility by allowing the user to interact with the visualization and choose from different estimation methods and visualizations. It also allows the user to easily toggle between correlation networks for samples under one condition and differential correlations between samples under two conditions. Moreover, the software facilitates integrative analysis of cross-correlation networks between two omics data sets.

Bayesian Structure Learning in Multi-layered Genomic Networks

Integrative network modeling of data arising from multiple genomic platforms provides insight into the holistic picture of the interactive system, as well as the flow of information across many disease domains including cancer. The basic data structure consists of a sequence of hierarchically ordered datasets for each individual subject, which facilitates integration of diverse inputs, such as genomic, transcriptomic, and proteomic data. A primary analytical task in such contexts is to model the layered architecture of networks where the vertices can be naturally partitioned into ordered layers, dictated by multiple platforms, and exhibit both undirected and directed relationships. We propose a multi-layered Gaussian graphical model (mlGGM) to investigate conditional independence structures in such multi-level genomic networks in human cancers. We implement a Bayesian node-wise selection (BANS) approach based on variable selection techniques that coherently accounts for the multiple types of dependencies in mlGGM; this flexible strategy exploits edge-specific prior knowledge and selects sparse and interpretable models. Through simulated data generated under various scenarios, we demonstrate that BANS outperforms other existing multivariate regression-based methodologies. Our integrative genomic network analysis for key signaling pathways across multiple cancer types highlights commonalities and differences of p53 integrative networks and epigenetic effects of BRCA2 on p53 and its interaction with T68 phosphorylated CHK2, that may have translational utilities of finding biomarkers and therapeutic targets.

The conserved Tpk1 regulates non-homologous end joining double-strand break repair by phosphorylation of Nej1, a homolog of the human XLF

The yeast cyclic AMP-dependent protein kinase A (PKA) is a ubiquitous serine-threonine kinase, encompassing three catalytic (Tpk1-3) and one regulatory (Bcy1) subunits. Evidence suggests PKA involvement in DNA damage checkpoint response, but how DNA repair pathways are regulated by PKA subunits remains inconclusive. Here, we report that deleting the tpk1 catalytic subunit reduces non-homologous end joining (NHEJ) efficiency, whereas tpk2-3 and bcy1 deletion does not. Epistatic analyses revealed that tpk1, as well as the DNA damage checkpoint kinase (dun1) and NHEJ factor (nej1), co-function in the same pathway, and parallel to the NHEJ factor yku80. Chromatin immunoprecipitation and resection data suggest that tpk1 deletion influences repair protein recruitments and DNA resection. Further, we show that Tpk1 phosphorylation of Nej1 at S298 (a Dun1 phosphosite) is indispensable for NHEJ repair and nuclear targeting of Nej1 and its binding partner Lif1. In mammalian cells, loss of PRKACB (human homolog of Tpk1) also reduced NHEJ efficiency, and similarly, PRKACB was found to phosphorylate XLF (a Nej1 human homolog) at S263, a corresponding residue of the yeast Nej1 S298. Together, our results uncover a new and conserved mechanism for Tpk1 and PRKACB in phosphorylating Nej1 (or XLF), which is critically required for NHEJ repair.

A computational strategy for metabolic network construction based on the overlapping ratio: Study of patients' metabolic responses to different dialysis patterns

BACKGROUND

Uremia is a worldwide epidemic disease and poses a serious threat to human health. Both maintenance hemodialysis (HD) and maintenance high flux hemodialysis (HFD) are common treatments for uremia and are generally used in clinical applications. In-depth exploration of patients' metabolic responses to different dialysis patterns can facilitate the understanding of pathological alterations associated with uremia and the effects of different dialysis methods on uremia, which may be used for future personalized therapy. However, due to variations of multiple factors (i.e., genetic, epigenetic and environment) in the process of disease treatments, identification of the similarities and differences in plasma metabolite changes in uremic patients in response to HD and HFD remains challenging.

METHODS

In this study, a computational strategy for metabolic network construction based on the overlapping ratio (MNC-OR) was proposed for disease treatment effect research. In MNC-OR, the overlapping ratio was introduced to measure metabolic reactions and to construct metabolic networks for analysis of different treatment options. Then, MNC-OR was employed to analyze HD-pattern-dependent changes in plasma metabolites to explore the pathological alterations associated with uremia and the effectiveness of different dialysis patterns (i.e., HD and HFD) on uremia. Based on the networks constructed by MNC-OR, two network analysis techniques, namely, similarity analysis and difference analysis of network topology, were used to find the similarity and differences in metabolic signals in patients under treatment with either HD or HFD, which can facilitate the understanding of pathological alterations associated with uremia and provide the guidance for personalized dialysis therapy.

RESULTS

Similarity analysis of network topology suggested that abnormal energy metabolism, gut metabolism and pyrimidine metabolism might occur in uremic patients, and maintenance of both HFD and HD therapies have beneficial effects on uremia. Then, difference analysis of network topology was employed to extract the crucial information related to HD-pattern-dependent changes in plasma metabolites. Experimental results indicated that the amino acid metabolism was closer to the normal status in HFD-treated patients; however, in HD-treated patients, the ability of antioxidation showed greater reduction, and the protein O-GlcNAcylation level was higher. Our findings demonstrate the potential of MNC-OR for explaining the metabolic similarities and differences of patients in response to different dialysis methods, thereby contributing to the guidance of personalized dialysis therapy.

Statistical mechanics meets single-cell biology

Single-cell omics is transforming our understanding of cell biology and disease, yet the systems-level analysis and interpretation of single-cell data faces many challenges. In this Perspective, we describe the impact that fundamental concepts from statistical mechanics, notably entropy, stochastic processes and critical phenomena, are having on single-cell data analysis. We further advocate the need for more bottom-up modelling of single-cell data and to embrace a statistical mechanics analysis paradigm to help attain a deeper understanding of single-cell systems biology.

Mapping the multiscale structure of biological systems

Biological systems are by nature multiscale, consisting of subsystems that factor into progressively smaller units in a deeply hierarchical structure. At any level of the hierarchy, an ever-increasing diversity of technologies can be applied to characterize the corresponding biological units and their relations, resulting in large networks of physical or functional proximities-e.g., proximities of amino acids within a protein, of proteins within a complex, or of cell types within a tissue. Here, we review general concepts and progress in using network proximity measures as a basis for creation of multiscale hierarchical maps of biological systems. We discuss the functionalization of these maps to create predictive models, including those useful in translation of genotype to phenotype, along with strategies for model visualization and challenges faced by multiscale modeling in the near future. Collectively, these approaches enable a unified hierarchical approach to biological data, with application from the molecular to the macroscopic.

Environmental robustness of the global yeast genetic interaction network

Phenotypes associated with genetic variants can be altered by interactions with other genetic variants (GxG), with the environment (GxE), or both (GxGxE). Yeast genetic interactions have been mapped on a global scale, but the environmental influence on the plasticity of genetic networks has not been examined systematically. To assess environmental rewiring of genetic networks, we examined 14 diverse conditions and scored 30,000 functionally representative yeast gene pairs for dynamic, differential interactions. Different conditions revealed novel differential interactions, which often uncovered functional connections between distantly related gene pairs. However, the majority of observed genetic interactions remained unchanged in different conditions, suggesting that the global yeast genetic interaction network is robust to environmental perturbation and captures the fundamental functional architecture of a eukaryotic cell.

Genetic interaction analysis in microbial pathogens: unravelling networks of pathogenesis, antimicrobial susceptibility and host interactions

Genetic interaction (GI) analysis is a powerful genetic strategy that analyzes the fitness and phenotypes of single- and double-gene mutant cells in order to dissect the epistatic interactions between genes, categorize genes into biological pathways, and characterize genes of unknown function. GI analysis has been extensively employed in model organisms for foundational, systems-level assessment of the epistatic interactions between genes. More recently, GI analysis has been applied to microbial pathogens and has been instrumental for the study of clinically important infectious organisms. Here, we review recent advances in systems-level GI analysis of diverse microbial pathogens, including bacterial and fungal species. We focus on important applications of GI analysis across pathogens, including GI analysis as a means to decipher complex genetic networks regulating microbial virulence, antimicrobial drug resistance and host-pathogen dynamics, and GI analysis as an approach to uncover novel targets for combination antimicrobial therapeutics. Together, this review bridges our understanding of GI analysis and complex genetic networks, with applications to diverse microbial pathogens, to further our understanding of virulence, the use of antimicrobial therapeutics and host-pathogen interactions. .

IMCC: A Novel Quantitative Approach Revealing Variation of Global Modular Map and Local Inter-Module Coordination Among Differential Drug's Targeted Cerebral Ischemic Networks

Stroke is a common disease characterized by multiple genetic dysfunctions. In this complex disease, detecting the strength of inter-module coordination (genetic community interaction) and subsequent modular rewiring is essential to characterize the reactive biosystematic variation (biosystematic perturbation) brought by multiple-target drugs, whose effects are achieved by hitting on a series of targets (target profile) jointly. Here, a quantitative approach for inter-module coordination and its transition, named as IMCC, was developed. Applying IMCC to mouse cerebral ischemia–related gene microarray, we investigated a holistic view of modular map and its rewiring from ischemic stroke to drugs (baicalin, BA; ursodeoxycholic acid, UA; and jasminoidin, JA) perturbation states and locally identified the cooperative pathological module pair and its dissection. Our result suggested the global modular map in cerebral ischemia exhibited a characteristic “core–periphery” architecture, and this architecture was rewired by the effective drugs heterogeneously: BA and UA converged modules into an intensively connected integrity, whereas JA diverged partial modules and widened the remaining inter-module paths. Locally, the PMP dissociation brought by drugs contributed to the reversion of the pathological condition: the focus of the cellular function shift from survival after nervous system injury into development and repair, including neurotrophin regulation, hormone releasing, and chemokine signaling activation. The core targets and mechanisms were validated by in vivo experiments. Overall, our result highlights the holistic inter-module coordination rearrangement rather than a target or a single module that brings phenotype alteration. This strategy may lead to systematically explore detailed variation of inter-module pharmacological action mode of multiple-target drugs, which is the principal problem of module pharmacology for network-based drug discovery.

FIREWORKS: a bottom-up approach to integrative coessentiality network analysis

Genetic coessentiality analysis, a computational approach which identifies genes sharing a common effect on cell fitness across large-scale screening datasets, has emerged as a powerful tool to identify functional relationships between human genes. However, widespread implementation of coessentiality to study individual genes and pathways is limited by systematic biases in existing coessentiality approaches and accessibility barriers for investigators without computational expertise. We created FIREWORKS, a method and interactive tool for the construction and statistical analysis of coessentiality networks centered around gene(s) provided by the user. FIREWORKS incorporates a novel bias reduction approach to reduce false discoveries, enables restriction of coessentiality analyses to custom subsets of cell lines, and integrates multiomic and drug-gene interaction datasets to investigate and target contextual gene essentiality. We demonstrate the broad utility of FIREWORKS through case vignettes investigating gene function and specialization, indirect therapeutic targeting of "undruggable" proteins, and context-specific rewiring of genetic networks.

Differential Co-Expression Analyses Allow the Identification of Critical Signalling Pathways Altered during Tumour Transformation and Progression

Biological systems respond to perturbations through the rewiring of molecular interactions, organised in gene regulatory networks (GRNs). Among these, the increasingly high availability of transcriptomic data makes gene co-expression networks the most exploited ones. Differential co-expression networks are useful tools to identify changes in response to an external perturbation, such as mutations predisposing to cancer development, and leading to changes in the activity of gene expression regulators or signalling. They can help explain the robustness of cancer cells to perturbations and identify promising candidates for targeted therapy, moreover providing higher specificity with respect to standard co-expression methods. Here, we comprehensively review the literature about the methods developed to assess differential co-expression and their applications to cancer biology. Via the comparison of normal and diseased conditions and of different tumour stages, studies based on these methods led to the definition of pathways involved in gene network reorganisation upon oncogenes' mutations and tumour progression, often converging on immune system signalling. A relevant implementation still lagging behind is the integration of different data types, which would greatly improve network interpretability. Most importantly, performance and predictivity evaluation of the large variety of mathematical models proposed would urgently require experimental validations and systematic comparisons. We believe that future work on differential gene co-expression networks, complemented with additional omics data and experimentally tested, will considerably improve our insights into the biology of tumours.

Computational identification of biomarker genes for lung cancer considering treatment and non-treatment studies

BACKGROUND

Lung cancer is the number one cancer killer in the world with more than 142,670 deaths estimated in the United States alone in the year 2019. Consequently, there is an overreaching need to identify the key biomarkers for lung cancer. The aim of this study is to computationally identify biomarker genes for lung cancer that can aid in its diagnosis and treatment. The gene expression profiles of two different types of studies, namely non-treatment and treatment, are considered for discovering biomarker genes. In non-treatment studies healthy samples are control and cancer samples are cases. Whereas, in treatment studies, controls are cancer cell lines without treatment and cases are cancer cell lines with treatment.

RESULTS

The Differentially Expressed Genes (DEGs) for lung cancer were isolated from Gene Expression Omnibus (GEO) database using R software tool GEO2R. A total of 407 DEGs (254 upregulated and 153 downregulated) from non-treatment studies and 547 DEGs (133 upregulated and 414 downregulated) from treatment studies were isolated. Two Cytoscape apps, namely, CytoHubba and MCODE, were used for identifying biomarker genes from functional networks developed using DEG genes. This study discovered two distinct sets of biomarker genes - one from non-treatment studies and the other from treatment studies, each set containing 16 genes. Survival analysis results show that most non-treatment biomarker genes have prognostic capability by indicating low-expression groups have higher chance of survival compare to high-expression groups. Whereas, most treatment biomarkers have prognostic capability by indicating high-expression groups have higher chance of survival compare to low-expression groups.

CONCLUSION

A computational framework is developed to identify biomarker genes for lung cancer using gene expression profiles. Two different types of studies - non-treatment and treatment - are considered for experiment. Most of the biomarker genes from non-treatment studies are part of mitosis and play vital role in DNA repair and cell-cycle regulation. Whereas, most of the biomarker genes from treatment studies are associated to ubiquitination and cellular response to stress. This study discovered a list of biomarkers, which would help experimental scientists to design a lab experiment for further exploration of detail dynamics of lung cancer development.

Whole transcriptomic network analysis using Co-expression Differential Network Analysis (CoDiNA)

Biological and medical sciences are increasingly acknowledging the significance of gene co-expression-networks for investigating complex-systems, phenotypes or diseases. Typically, complex phenotypes are investigated under varying conditions. While approaches for comparing nodes and links in two networks exist, almost no methods for the comparison of multiple networks are available and—to best of our knowledge—no comparative method allows for whole transcriptomic network analysis. However, it is the aim of many studies to compare networks of different conditions, for example, tissues, diseases, treatments, time points, or species. Here we present a method for the systematic comparison of an unlimited number of networks, with unlimited number of transcripts: Co-expression Differential Network Analysis (CoDiNA). In particular, CoDiNA detects links and nodes that are common, specific or different among the networks. We developed a statistical framework to normalize between these different categories of common or changed network links and nodes, resulting in a comprehensive network analysis method, more sophisticated than simply comparing the presence or absence of network nodes. Applying CoDiNA to a neurogenesis study we identified candidate genes involved in neuronal differentiation. We experimentally validated one candidate, demonstrating that its overexpression resulted in a significant disturbance in the underlying gene regulatory network of neurogenesis. Using clinical studies, we compared whole transcriptome co-expression networks from individuals with or without HIV and active tuberculosis (TB) and detected signature genes specific to HIV. Furthermore, analyzing multiple cancer transcription factor (TF) networks, we identified common and distinct features for particular cancer types. These CoDiNA applications demonstrate the successful detection of genes associated with specific phenotypes. Moreover, CoDiNA can also be used for comparing other types of undirected networks, for example, metabolic, protein-protein interaction, ecological and psychometric networks. CoDiNA is publicly available as an R package in CRAN (https://CRAN.R-project.org/package=CoDiNA).

Computational solutions for modeling and controlling plant response to abiotic stresses: a review with focus on iron deficiency

Computational solutions enable plant scientists to model protein-mediated stress responses and characterize novel gene functions that coordinate responses to a variety of abiotic stress conditions. Recently, density functional theory was used to study proteins active sites and elucidate enzyme conversion mechanisms involved in iron deficiency responsive signaling pathways. Computational approaches for protein homology modeling and the kinetic modeling of signaling pathways have also resolved the identity and function in proteins involved in iron deficiency signaling pathways. Significant changes in gene relationships under other stress conditions, such as heat or drought stress, have been recently identified using differential network analysis, suggesting that stress tolerance is achieved through asynchronous control. Moreover, the increasing development and use of statistical modeling and systematic modeling of transcriptomic data have provided significant insight into the gene regulatory mechanisms associated with abiotic stress responses. These types of in silico approaches have facilitated the plant science community's future goals of developing multi-scale models of responses to iron deficiency stress and other abiotic stress conditions.

Genome-Wide Dynamic Evaluation of the UV-Induced DNA Damage Response

Genetic screens in Saccharomyces cerevisiae have allowed for the identification of many genes as sensors or effectors of DNA damage, typically by comparing the fitness of genetic mutants in the presence or absence of DNA-damaging treatments. However, these static screens overlook the dynamic nature of DNA damage response pathways, missing time-dependent or transient effects. Here, we examine gene dependencies in the dynamic response to ultraviolet radiation-induced DNA damage by integrating ultra-high-density arrays of 6144 diploid gene deletion mutants with high-frequency time-lapse imaging. We identify 494 ultraviolet radiation response genes which, in addition to recovering molecular pathways and protein complexes previously annotated to DNA damage repair, include components of the CCR4-NOT complex, tRNA wobble modification, autophagy, and, most unexpectedly, 153 nuclear-encoded mitochondrial genes. Notably, mitochondria-deficient strains present time-dependent insensitivity to ultraviolet radiation, posing impaired mitochondrial function as a protective factor in the ultraviolet radiation response.

Regulation of alternative polyadenylation in the yeast Saccharomyces cerevisiae by histone H3K4 and H3K36 methyltransferases

Adjusting DNA structure via epigenetic modifications, and altering polyadenylation (pA) sites at which precursor mRNA is cleaved and polyadenylated, allows cells to quickly respond to environmental stress. Since polyadenylation occurs co-transcriptionally, and specific patterns of nucleosome positioning and chromatin modifications correlate with pA site usage, epigenetic factors potentially affect alternative polyadenylation (APA). We report that the histone H3K4 methyltransferase Set1, and the histone H3K36 methyltransferase Set2, control choice of pA site in Saccharomyces cerevisiae, a powerful model for studying evolutionarily conserved eukaryotic processes. Deletion of SET1 or SET2 causes an increase in serine-2 phosphorylation within the C-terminal domain of RNA polymerase II (RNAP II) and in the recruitment of the cleavage/polyadenylation complex, both of which could cause the observed switch in pA site usage. Chemical inhibition of TOR signaling, which causes nutritional stress, results in Set1- and Set2-dependent APA. In addition, Set1 and Set2 decrease efficiency of using single pA sites, and control nucleosome occupancy around pA sites. Overall, our study suggests that the methyltransferases Set1 and Set2 regulate APA induced by nutritional stress, affect the RNAP II C-terminal domain phosphorylation at Ser2, and control recruitment of the 3′ end processing machinery to the vicinity of pA sites.

An integrative framework of heterogeneous genomic data for cancer dynamic modules based on matrix decomposition

Cancer progression is dynamic, and tracking dynamic modules is promising for cancer diagnosis and therapy. Accumulated genomic data provide us an opportunity to investigate the underlying mechanisms of cancers. However, as far as we know, no algorithm has been designed for dynamic modules by integrating heterogeneous omics data. To address this issue, we propose an integrative framework for dynamic module detection based on regularized nonnegative matrix factorization method (DrNMF) by integrating the gene expression and protein interaction network. To remove the heterogeneity of genomic data, we divide the samples of expression profiles into groups to construct gene co-expression networks. To characterize the dynamics of modules, the temporal smoothness framework is adopted, in which the gene co-expression network at the previous stage and protein interaction network are incorporated into the objective function of DrNMF via regularization. The experimental results demonstrate that DrNMF is superior to state-of-the-art methods in terms of accuracy. For breast cancer data, the obtained dynamic modules are more enriched by the known pathways, and can be used to predict the stages of cancers and survival time of patients. The proposed model and algorithm provide an effective integrative analysis of heterogeneous genomic data for cancer progression.

Dynamic rewiring of the human interactome by interferon signaling

BACKGROUND

The type I interferon (IFN) response is an ancient pathway that protects cells against viral pathogens by inducing the transcription of hundreds of IFN-stimulated genes. Comprehensive catalogs of IFN-stimulated genes have been established across species and cell types by transcriptomic and biochemical approaches, but their antiviral mechanisms remain incompletely characterized. Here, we apply a combination of quantitative proteomic approaches to describe the effects of IFN signaling on the human proteome, and apply protein correlation profiling to map IFN-induced rearrangements in the human protein-protein interaction network.

RESULTS

We identify > 26,000 protein interactions in IFN-stimulated and unstimulated cells, many of which involve proteins associated with human disease and are observed exclusively within the IFN-stimulated network. Differential network analysis reveals interaction rewiring across a surprisingly broad spectrum of cellular pathways in the antiviral response. We identify IFN-dependent protein-protein interactions mediating novel regulatory mechanisms at the transcriptional and translational levels, with one such interaction modulating the transcriptional activity of STAT1. Moreover, we reveal IFN-dependent changes in ribosomal composition that act to buffer IFN-stimulated gene protein synthesis.

CONCLUSIONS

Our map of the IFN interactome provides a global view of the complex cellular networks activated during the antiviral response, placing IFN-stimulated genes in a functional context, and serves as a framework to understand how these networks are dysregulated in autoimmune or inflammatory disease.

The kinase Isr1 negatively regulates hexosamine biosynthesis in S. cerevisiae

The S. cerevisiae ISR1 gene encodes a putative kinase with no ascribed function. Here, we show that Isr1 acts as a negative regulator of the highly-conserved hexosamine biosynthesis pathway (HBP), which converts glucose into uridine diphosphate N-acetylglucosamine (UDP-GlcNAc), the carbohydrate precursor to protein glycosylation, GPI-anchor formation, and chitin biosynthesis. Overexpression of ISR1 is lethal and, at lower levels, causes sensitivity to tunicamycin and resistance to calcofluor white, implying impaired protein glycosylation and reduced chitin deposition. Gfa1 is the first enzyme in the HBP and is conserved from bacteria and yeast to humans. The lethality caused by ISR1 overexpression is rescued by co-overexpression of GFA1 or exogenous glucosamine, which bypasses GFA1's essential function. Gfa1 is phosphorylated in an Isr1-dependent fashion and mutation of Isr1-dependent sites ameliorates the lethality associated with ISR1 overexpression. Isr1 contains a phosphodegron that is phosphorylated by Pho85 and subsequently ubiquitinated by the SCF-Cdc4 complex, largely confining Isr1 protein levels to the time of bud emergence. Mutation of this phosphodegron stabilizes Isr1 and recapitulates the overexpression phenotypes. As Pho85 is a cell cycle and nutrient responsive kinase, this tight regulation of Isr1 may serve to dynamically regulate flux through the HBP and modulate how the cell's energy resources are converted into structural carbohydrates in response to changing cellular needs.

Integration of multiple biological contexts reveals principles of synthetic lethality that affect reproducibility

Synthetic lethal screens have the potential to identify new vulnerabilities incurred by specific cancer mutations but have been hindered by lack of agreement between studies. In the case of KRAS, we identify that published synthetic lethal screen hits significantly overlap at the pathway rather than gene level. Analysis of pathways encoded as protein networks could identify synthetic lethal candidates that are more reproducible than those previously reported. Lack of overlap likely stems from biological rather than technical limitations as most synthetic lethal phenotypes are strongly modulated by changes in cellular conditions or genetic context, the latter determined using a pairwise genetic interaction map that identifies numerous interactions that suppress synthetic lethal effects. Accounting for pathway, cellular and genetic context nominates a DNA repair dependency in KRAS-mutant cells, mediated by a network containing BRCA1. We provide evidence for why most reported synthetic lethals are not reproducible which is addressable using a multi-faceted testing framework.

Dynamic genetic architecture of yeast response to environmental perturbation shed light on origin of cryptic genetic variation

Cryptic genetic variation could arise from, for example, Gene-by-Gene (G-by-G) or Gene-by-Environment (G-by-E) interactions. The underlying molecular mechanisms and how they influence allelic effects and the genetic variance of complex traits is largely unclear. Here, we empirically explored the role of environmentally influenced epistasis on the suppression and release of cryptic variation by reanalysing a dataset of 4,390 haploid yeast segregants phenotyped on 20 different media. The focus was on 130 epistatic loci, each contributing to segregant growth in at least one environment and that together explained most (69–100%) of the narrow sense heritability of growth in the individual environments. We revealed that the epistatic growth network reorganised upon environmental changes to alter the estimated marginal (additive) effects of the individual loci, how multi-locus interactions contributed to individual segregant growth and the level of expressed genetic variance in growth. The estimated additive effects varied most across environments for loci that were highly interactive network hubs in some environments but had few or no interactors in other environments, resulting in changes in total genetic variance across environments. This environmentally dependent epistasis was thus an important mechanism for the suppression and release of cryptic variation in this population. Our findings increase the understanding of the complex genetic mechanisms leading to cryptic variation in populations, providing a basis for future studies on the genetic maintenance of trait robustness and development of genetic models for studying and predicting selection responses for quantitative traits in breeding and evolution.

Many biological traits are polygenic, with complex interplay between underlying genes and the surrounding environment. As a result, individuals with the same allele might have distinctive phenotypes due to differences in the polygenic background and/or the environment. Such differences often create additional genetic variation that is highly relevant to quantitative and evolutionary genetics by limiting our ability to accurately predict the phenotypes in medical or agricultural applications and providing opportunities for long term evolution. Previously, yeast growth regulating genes were found to be organised in large interacting networks. Here, we found that these networks were reorganised upon environmental changes, and that this resulted in altered effect sizes of individual genes, and how the whole network contributed to growth and the level of total genetic variance, providing a basis for future studies on the genetic maintenance of trait robustness and development of genetic models for studying and predicting selection responses for quantitative traits.

Genetic interactions derived from high-throughput phenotyping of 6589 yeast cell cycle mutants

Over the last 30 years, computational biologists have developed increasingly realistic mathematical models of the regulatory networks controlling the division of eukaryotic cells. These models capture data resulting from two complementary experimental approaches: low-throughput experiments aimed at extensively characterizing the functions of small numbers of genes, and large-scale genetic interaction screens that provide a systems-level perspective on the cell division process. The former is insufficient to capture the interconnectivity of the genetic control network, while the latter is fraught with irreproducibility issues. Here, we describe a hybrid approach in which the 630 genetic interactions between 36 cell-cycle genes are quantitatively estimated by high-throughput phenotyping with an unprecedented number of biological replicates. Using this approach, we identify a subset of high-confidence genetic interactions, which we use to refine a previously published mathematical model of the cell cycle. We also present a quantitative dataset of the growth rate of these mutants under six different media conditions in order to inform future cell cycle models.

Personalized Network Modeling of the Pan-Cancer Patient and Cell Line Interactome

PURPOSE

Personalized network inference on diverse clinical and in vitro model systems across cancer types can be used to delineate specific regulatory mechanisms, uncover drug targets and pathways, and develop individualized predictive models in cancer.

METHODS

We developed TransPRECISE (personalized cancer-specific integrated network estimation model), a multiscale Bayesian network modeling framework, to analyze the pan-cancer patient and cell line interactome to identify differential and conserved intrapathway activities, to globally assess cell lines as representative models for patients, and to develop drug sensitivity prediction models. We assessed pan-cancer pathway activities for a large cohort of patient samples (> 7,700) from the Cancer Proteome Atlas across ≥ 30 tumor types, a set of 640 cancer cell lines from the MD Anderson Cell Lines Project spanning 16 lineages, and ≥ 250 cell lines' response to > 400 drugs.

RESULTS

TransPRECISE captured differential and conserved proteomic network topologies and pathway circuitry between multiple patient and cell line lineages: ovarian and kidney cancers shared high levels of connectivity in the hormone receptor and receptor tyrosine kinase pathways, respectively, between the two model systems. Our tumor stratification approach found distinct clinical subtypes of the patients represented by different sets of cell lines: patients with head and neck tumors were classified into two different subtypes that are represented by head and neck and esophagus cell lines and had different prognostic patterns (456 v 654 days of median overall survival; P = .02). High predictive accuracy was observed for drug sensitivities in cell lines across multiple drugs (median area under the receiver operating characteristic curve > 0.8) using Bayesian additive regression tree models with TransPRECISE pathway scores.

CONCLUSION

Our study provides a generalizable analytic framework to assess the translational potential of preclinical model systems and to guide pathway-based personalized medical decision making, integrating genomic and molecular data across model systems.

Genetic interaction profiles of regulatory kinases differ between environmental conditions and cellular states

Cell growth and quiescence in eukaryotic cells is controlled by an evolutionarily conserved network of signaling pathways. Signal transduction networks operate to modulate a wide range of cellular processes and physiological properties when cells exit proliferative growth and initiate a quiescent state. How signaling networks function to respond to diverse signals that result in cell cycle exit and establishment of a quiescent state is poorly understood. Here, we studied the function of signaling pathways in quiescent cells using global genetic interaction mapping in the model eukaryotic cell, Saccharomyces cerevisiae (budding yeast). We performed pooled analysis of genotypes using molecular barcode sequencing (Bar-seq) to test the role of ~4,000 gene deletion mutants and ~12,000 pairwise interactions between all non-essential genes and the protein kinase genes TOR1, RIM15, and PHO85 in three different nutrient-restricted conditions in both proliferative and quiescent cells. We detect up to 10-fold more genetic interactions in quiescent cells than proliferative cells. We find that both individual gene effects and genetic interaction profiles vary depending on the specific pro-quiescence signal. The master regulator of quiescence, RIM15, shows distinct genetic interaction profiles in response to different starvation signals. However, vacuole-related functions show consistent genetic interactions with RIM15 in response to different starvation signals, suggesting that RIM15 integrates diverse signals to maintain protein homeostasis in quiescent cells. Our study expands genome-wide genetic interaction profiling to additional conditions, and phenotypes, and highlights the conditional dependence of epistasis.