A systems-level gene regulatory network model for Plasmodium falciparum.

Many of the gene regulatory processes of Plasmodium falciparum, the deadliest malaria parasite, remain poorly understood. To develop a comprehensive guide for exploring this organism's gene regulatory network, we generated a systems-level model of P. falciparum gene regulation using a well-validated, machine-learning approach for predicting interactions between transcription regulators and their targets. The resulting network accurately predicts expression levels of transcriptionally coherent gene regulatory programs in independent transcriptomic data sets from parasites collected by different research groups in diverse laboratory and field settings. Thus, our results indicate that our gene regulatory model has predictive power and utility as a hypothesis-generating tool for illuminating clinically relevant gene regulatory mechanisms within P. falciparum. Using the set of regulatory programs we identified, we also investigated correlates of artemisinin resistance based on gene expression coherence. We report that resistance is associated with incoherent expression across many regulatory programs, including those controlling genes associated with erythrocyte-host engagement. These results suggest that parasite populations with reduced artemisinin sensitivity are more transcriptionally heterogenous. This pattern is consistent with a model where the parasite utilizes bet-hedging strategies to diversify the population, rendering a subpopulation more able to navigate drug treatment.

Bone marrow microenvironments that contribute to patient outcomes in newly diagnosed multiple myeloma: A cohort study of patients in the Total Therapy clinical trials.

BACKGROUND: The tumor microenvironment (TME) is increasingly appreciated as an important determinant of cancer outcome, including in multiple myeloma (MM). However, most myeloma microenvironment studies have been based on bone marrow (BM) aspirates, which often do not fully reflect the cellular content of BM tissue itself. To address this limitation in myeloma research, we systematically characterized the whole bone marrow (WBM) microenvironment during premalignant, baseline, on treatment, and post-treatment phases. METHODS AND FINDINGS: Between 2004 and 2019, 998 BM samples were taken from 436 patients with newly diagnosed MM (NDMM) at the University of Arkansas for Medical Sciences in Little Rock, Arkansas, United States of America. These patients were 61% male and 39% female, 89% White, 8% Black, and 3% other/refused, with a mean age of 58 years. Using WBM and matched cluster of differentiation (CD)138-selected tumor gene expression to control for tumor burden, we identified a subgroup of patients with an adverse TME associated with 17 fewer months of progression-free survival (PFS) (95% confidence interval [CI] 5-29, 49-69 versus 70-82 months, χ2 p = 0.001) and 15 fewer months of overall survival (OS; 95% CI -1 to 31, 92-120 versus 113-129 months, χ2 p = 0.036). Using immunohistochemistry-validated computational tools that identify distinct cell types from bulk gene expression, we showed that the adverse outcome was correlated with elevated CD8+ T cell and reduced granulocytic cell proportions. This microenvironment develops during the progression of premalignant to malignant disease and becomes less prevalent after therapy, in which it is associated with improved outcomes. In patients with quantified International Staging System (ISS) stage and 70-gene Prognostic Risk Score (GEP-70) scores, taking the microenvironment into consideration would have identified an additional 40 out of 290 patients (14%, premutation p = 0.001) with significantly worse outcomes (PFS, 95% CI 6-36, 49-73 versus 74-90 months) who were not identified by existing clinical (ISS stage III) and tumor (GEP-70) criteria as high risk. The main limitations of this study are that it relies on computationally identified cell types and that patients were treated with thalidomide rather than current therapies. CONCLUSIONS: In this study, we observe that granulocyte signatures in the MM TME contribute to a more accurate prognosis. This implies that future researchers and clinicians treating patients should quantify TME components, in particular monocytes and granulocytes, which are often ignored in microenvironment studies.

Severe adult malaria is associated with specific PfEMP1 adhesion types and high parasite biomass.

The interplay between cellular and molecular determinants that lead to severe malaria in adults is unexplored. Here, we analyzed parasite virulence factors in an infected adult population in India and investigated whether severe malaria isolates impair endothelial protein C receptor (EPCR), a protein involved in coagulation and endothelial barrier permeability. Severe malaria isolates overexpressed specific members of the Plasmodium falciparum var gene/PfEMP1 (P. falciparum erythrocyte membrane protein 1) family that bind EPCR, including DC8 var genes that have previously been linked to severe pediatric malaria. Machine learning analysis revealed that DC6- and DC8-encoding var transcripts in combination with high parasite biomass were the strongest indicators of patient hospitalization and disease severity. We found that DC8 CIDRα1 domains from severe malaria isolates had substantial differences in EPCR binding affinity and blockade activity for its ligand activated protein C. Additionally, even a low level of inhibition exhibited by domains from two cerebral malaria isolates was sufficient to interfere with activated protein C-barrier protective activities in human brain endothelial cells. Our findings demonstrate an interplay between parasite biomass and specific PfEMP1 adhesion types in the development of adult severe malaria, and indicate that low impairment of EPCR function may contribute to parasite virulence.

Combining inferred regulatory and reconstructed metabolic networks enhances phenotype prediction in yeast.

Gene regulatory and metabolic network models have been used successfully in many organisms, but inherent differences between them make networks difficult to integrate. Probabilistic Regulation Of Metabolism (PROM) provides a partial solution, but it does not incorporate network inference and underperforms in eukaryotes. We present an Integrated Deduced And Metabolism (IDREAM) method that combines statistically inferred Environment and Gene Regulatory Influence Network (EGRIN) models with the PROM framework to create enhanced metabolic-regulatory network models. We used IDREAM to predict phenotypes and genetic interactions between transcription factors and genes encoding metabolic activities in the eukaryote, Saccharomyces cerevisiae. IDREAM models contain many fewer interactions than PROM and yet produce significantly more accurate growth predictions. IDREAM consistently outperformed PROM using any of three popular yeast metabolic models and across three experimental growth conditions. Importantly, IDREAM's enhanced accuracy makes it possible to identify subtle synthetic growth defects. With experimental validation, these novel genetic interactions involving the pyruvate dehydrogenase complex suggested a new role for fatty acid-responsive factor Oaf1 in regulating acetyl-CoA production in glucose grown cells.

Predicting antidisease immunity using proteome arrays and sera from children naturally exposed to malaria.

Malaria remains one of the most prevalent and lethal human infectious diseases worldwide. A comprehensive characterization of antibody responses to blood stage malaria is essential to support the development of future vaccines, sero-diagnostic tests, and sero-surveillance methods. We constructed a proteome array containing 4441 recombinant proteins expressed by the blood stages of the two most common human malaria parasites, P. falciparum (Pf) and P. vivax (Pv), and used this array to screen sera of Papua New Guinea children infected with Pf, Pv, or both (Pf/Pv) that were either symptomatic (febrile), or asymptomatic but had parasitemia detectable via microscopy or PCR. We hypothesized that asymptomatic children would develop antigen-specific antibody profiles associated with antidisease immunity, as compared with symptomatic children. The sera from these children recognized hundreds of the arrayed recombinant Pf and Pv proteins. In general, responses in asymptomatic children were highest in those with high parasitemia, suggesting that antibody levels are associated with parasite burden. In contrast, symptomatic children carried fewer antibodies than asymptomatic children with infections detectable by microscopy, particularly in Pv and Pf/Pv groups, suggesting that antibody production may be impaired during symptomatic infections. We used machine-learning algorithms to investigate the relationship between antibody responses and symptoms, and we identified antibody responses to sets of Plasmodium proteins that could predict clinical status of the donors. Several of these antibody responses were identified by multiple comparisons, including those against members of the serine enriched repeat antigen family and merozoite protein 4. Interestingly, both P. falciparum serine enriched repeat antigen-5 and merozoite protein 4 have been previously investigated for use in vaccines. This machine learning approach, never previously applied to proteome arrays, can be used to generate a list of potential seroprotective and/or diagnostic antigens candidates that can be further evaluated in longitudinal studies.

Molecular mechanisms of system responses to novel stimuli are predictable from public data.

Systems scale models provide the foundation for an effective iterative cycle between hypothesis generation, experiment and model refinement. Such models also enable predictions facilitating the understanding of biological complexity and the control of biological systems. Here, we demonstrate the reconstruction of a globally predictive gene regulatory model from public data: a model that can drive rational experiment design and reveal new regulatory mechanisms underlying responses to novel environments. Specifically, using ∼ 1500 publically available genome-wide transcriptome data sets from Saccharomyces cerevisiae, we have reconstructed an environment and gene regulatory influence network that accurately predicts regulatory mechanisms and gene expression changes on exposure of cells to completely novel environments. Focusing on transcriptional networks that induce peroxisomes biogenesis, the model-guided experiments allow us to expand a core regulatory network to include novel transcriptional influences and linkage across signaling and transcription. Thus, the approach and model provides a multi-scalar picture of gene dynamics and are powerful resources for exploiting extant data to rationally guide experimentation. The techniques outlined here are generally applicable to any biological system, which is especially important when experimental systems are challenging and samples are difficult and expensive to obtain-a common problem in laboratory animal and human studies.

PhosphoChain: a novel algorithm to predict kinase and phosphatase networks from high-throughput expression data.

MOTIVATION: Protein phosphorylation is critical for regulating cellular activities by controlling protein activities, localization and turnover, and by transmitting information within cells through signaling networks. However, predictions of protein phosphorylation and signaling networks remain a significant challenge, lagging behind predictions of transcriptional regulatory networks into which they often feed. RESULTS: We developed PhosphoChain to predict kinases, phosphatases and chains of phosphorylation events in signaling networks by combining mRNA expression levels of regulators and targets with a motif detection algorithm and optional prior information. PhosphoChain correctly reconstructed ∼78% of the yeast mitogen-activated protein kinase pathway from publicly available data. When tested on yeast phosphoproteomic data from large-scale mass spectrometry experiments, PhosphoChain correctly identified ∼27% more phosphorylation sites than existing motif detection tools (NetPhosYeast and GPS2.0), and predictions of kinase-phosphatase interactions overlapped with ∼59% of known interactions present in yeast databases. PhosphoChain provides a valuable framework for predicting condition-specific phosphorylation events from high-throughput data. AVAILABILITY: PhosphoChain is implemented in Java and available at http://virgo.csie.ncku.edu.tw/PhosphoChain/ or http://aitchisonlab.com/PhosphoChain

Transcriptomic classification of diffuse large B-cell lymphoma identifies a high-risk activated B-cell-like subpopulation with targetable MYC dysregulation.

Immunochemotherapy has been the mainstay of treatment for newly diagnosed diffuse large B-cell lymphoma (ndDLBCL) yet is inadequate for many patients. In this work, we perform unsupervised clustering on transcriptomic features from a large cohort of ndDLBCL patients and identify seven clusters, one called A7 with poor prognosis, and develop a classifier to identify these clusters in independent ndDLBCL cohorts. This high-risk cluster is enriched for activated B-cell cell-of-origin, low immune infiltration, high MYC expression, and copy number aberrations. We compare and contrast our methodology with recent DLBCL classifiers to contextualize our clusters and show improved prognostic utility. Finally, using pre-clinical models, we demonstrate a mechanistic rationale for IKZF1/3 degraders such as lenalidomide to overcome the low immune infiltration phenotype of A7 by inducing T-cell trafficking into tumors and upregulating MHC I and II on tumor cells, and demonstrate that TCF4 is an important regulator of MYC-related biology in A7.

The functional role of SUMO E3 ligase Mms21p in the maintenance of subtelomeric silencing in budding yeast.

In Saccharomyces cerevisiae, subtelomeric silencing is involved in the propagation of Silent Information Regulator (SIR) proteins toward euchromatin. Numerous mechanisms are involved in antagonizing the local spread of Sir-dependent silent chromatin into neighboring euchromatin. Here, we identified a novel role for sumoylation E3 ligase Mms21 in the maintenance of subtelomeric silencing. We found that disruption of E3 ligase activity of Mms21 results in the de-repression of subtelomeric silencing. Deletion of E3 ligase domain of Mms21 led to decreased binding of Sir2p, Sir3p and Sir4 at subtelomeric chromatins and increased H3K4 tri-methylation at telomere-distal euchromatin regions, correlating with increased gene expression in two subtelomeric reporter genes. In addition, a mms21Δsl mutant caused a severe growth defect in combination with htz1Δ deletion and showed an enhanced association of Htz1 with telomere proximal regions. Taken together, our findings suggest an important role of Mms21p; it contributes to subtelomeric silencing during the formation of a heterochromatin boundary.

All-codon scanning identifies p53 cancer rescue mutations.

In vitro scanning mutagenesis strategies are valuable tools to identify critical residues in proteins and to generate proteins with modified properties. We describe the fast and simple All-Codon Scanning (ACS) strategy that creates a defined gene library wherein each individual codon within a specific target region is changed into all possible codons with only a single codon change per mutagenesis product. ACS is based on a multiplexed overlapping mutagenesis primer design that saturates only the targeted gene region with single codon changes. We have used ACS to produce single amino-acid changes in small and large regions of the human tumor suppressor protein p53 to identify single amino-acid substitutions that can restore activity to inactive p53 found in human cancers. Single-tube reactions were used to saturate defined 30-nt regions with all possible codon changes. The same technique was used in 20 parallel reactions to scan the 600-bp fragment encoding the entire p53 core domain. Identification of several novel p53 cancer rescue mutations demonstrated the utility of the ACS approach. ACS is a fast, simple and versatile method, which is useful for protein structure-function analyses and protein design or evolution problems.