Comparative Flavivirus-Host Protein Interaction Mapping Reveals Mechanisms of Dengue and Zika Virus Pathogenesis.

Mosquito-borne flaviviruses, including dengue virus (DENV) and Zika virus (ZIKV), are a growing public health concern. Systems-level analysis of how flaviviruses hijack cellular processes through virus-host protein-protein interactions (PPIs) provides information about their replication and pathogenic mechanisms. We used affinity purification-mass spectrometry (AP-MS) to compare flavivirus-host interactions for two viruses (DENV and ZIKV) in two hosts (human and mosquito). Conserved virus-host PPIs revealed that the flavivirus NS5 protein suppresses interferon stimulated genes by inhibiting recruitment of the transcription complex PAF1C and that chemical modulation of SEC61 inhibits DENV and ZIKV replication in human and mosquito cells. Finally, we identified a ZIKV-specific interaction between NS4A and ANKLE2, a gene linked to hereditary microcephaly, and showed that ZIKV NS4A causes microcephaly in Drosophila in an ANKLE2-dependent manner. Thus, comparative flavivirus-host PPI mapping provides biological insights and, when coupled with in vivo models, can be used to unravel pathogenic mechanisms.

From structure to systems: high-resolution, quantitative genetic analysis of RNA polymerase II.

RNA polymerase II (RNAPII) lies at the core of dynamic control of gene expression. Using 53 RNAPII point mutants, we generated a point mutant epistatic miniarray profile (pE-MAP) comprising ∼60,000 quantitative genetic interactions in Saccharomyces cerevisiae. This analysis enabled functional assignment of RNAPII subdomains and uncovered connections between individual regions and other protein complexes. Using splicing microarrays and mutants that alter elongation rates in vitro, we found an inverse relationship between RNAPII speed and in vivo splicing efficiency. Furthermore, the pE-MAP classified fast and slow mutants that favor upstream and downstream start site selection, respectively. The striking coordination of polymerization rate with transcription initiation and splicing suggests that transcription rate is tuned to regulate multiple gene expression steps. The pE-MAP approach provides a powerful strategy to understand other multifunctional machines at amino acid resolution.

A Quantitative Genetic Interaction Map of HIV Infection.

We have developed a platform for quantitative genetic interaction mapping using viral infectivity as a functional readout and constructed a viral host-dependency epistasis map (vE-MAP) of 356 human genes linked to HIV function, comprising >63,000 pairwise genetic perturbations. The vE-MAP provides an expansive view of the genetic dependencies underlying HIV infection and can be used to identify drug targets and study viral mutations. We found that the RNA deadenylase complex, CNOT, is a central player in the vE-MAP and show that knockout of CNOT1, 10, and 11 suppressed HIV infection in primary T cells by upregulating innate immunity pathways. This phenotype was rescued by deletion of IRF7, a transcription factor regulating interferon-stimulated genes, revealing a previously unrecognized host signaling pathway involved in HIV infection. The vE-MAP represents a generic platform that can be used to study the global effects of how different pathogens hijack and rewire the host during infection.

Phenotypic landscape of a bacterial cell.

The explosion of sequence information in bacteria makes developing high-throughput, cost-effective approaches to matching genes with phenotypes imperative. Using E. coli as proof of principle, we show that combining large-scale chemical genomics with quantitative fitness measurements provides a high-quality data set rich in discovery. Probing growth profiles of a mutant library in hundreds of conditions in parallel yielded > 10,000 phenotypes that allowed us to study gene essentiality, discover leads for gene function and drug action, and understand higher-order organization of the bacterial chromosome. We highlight new information derived from the study, including insights into a gene involved in multiple antibiotic resistance and the synergy between a broadly used combinatory antibiotic therapy, trimethoprim and sulfonamides. This data set, publicly available at http://ecoliwiki.net/tools/chemgen/, is a valuable resource for both the microbiological and bioinformatic communities, as it provides high-confidence associations between hundreds of annotated and uncharacterized genes as well as inferences about the mode of action of several poorly understood drugs.

An Mtb-Human Protein-Protein Interaction Map Identifies a Switch between Host Antiviral and Antibacterial Responses.

Although macrophages are armed with potent antibacterial functions, Mycobacterium tuberculosis (Mtb) replicates inside these innate immune cells. Determinants of macrophage intrinsic bacterial control, and the Mtb strategies to overcome them, are poorly understood. To further study these processes, we used an affinity tag purification mass spectrometry (AP-MS) approach to identify 187 Mtb-human protein-protein interactions (PPIs) involving 34 secreted Mtb proteins. This interaction map revealed two factors involved in Mtb pathogenesis-the secreted Mtb protein, LpqN, and its binding partner, the human ubiquitin ligase CBL. We discovered that an lpqN Mtb mutant is attenuated in macrophages, but growth is restored when CBL is removed. Conversely, Cbl-/- macrophages are resistant to viral infection, indicating that CBL regulates cell-intrinsic polarization between antibacterial and antiviral immunity. Collectively, these findings illustrate the utility of this Mtb-human PPI map for developing a deeper understanding of the intricate interactions between Mtb and its host.

The histone variant H2A.Z promotes splicing of weak introns.

Multiple lines of evidence implicate chromatin in the regulation of premessenger RNA (pre-mRNA) splicing. However, the influence of chromatin factors on cotranscriptional splice site usage remains unclear. Here we investigated the function of the highly conserved histone variant H2A.Z in pre-mRNA splicing using the intron-rich model yeast Schizosaccharomyces pombe Using epistatic miniarray profiles (EMAPs) to survey the genetic interaction landscape of the Swr1 nucleosome remodeling complex, which deposits H2A.Z, we uncovered evidence for functional interactions with components of the spliceosome. In support of these genetic connections, splicing-specific microarrays show that H2A.Z and the Swr1 ATPase are required during temperature stress for the efficient splicing of a subset of introns. Notably, affected introns are enriched for H2A.Z occupancy and more likely to contain nonconsensus splice sites. To test the significance of the latter correlation, we mutated the splice sites in an affected intron to consensus and found that this suppressed the requirement for H2A.Z in splicing of that intron. These data suggest that H2A.Z occupancy promotes cotranscriptional splicing of suboptimal introns that may otherwise be discarded via proofreading ATPases. Consistent with this model, we show that overexpression of splicing ATPase Prp16 suppresses both the growth and splicing defects seen in the absence of H2A.Z.

Functional organization of the S. cerevisiae phosphorylation network.

Reversible protein phosphorylation is a signaling mechanism involved in all cellular processes. To create a systems view of the signaling apparatus in budding yeast, we generated an epistatic miniarray profile (E-MAP) comprised of 100,000 pairwise, quantitative genetic interactions, including virtually all protein and small-molecule kinases and phosphatases as well as key cellular regulators. Quantitative genetic interaction mapping reveals factors working in compensatory pathways (negative genetic interactions) or those operating in linear pathways (positive genetic interactions). We found an enrichment of positive genetic interactions between kinases, phosphatases, and their substrates. In addition, we assembled a higher-order map from sets of three genes that display strong interactions with one another: triplets enriched for functional connectivity. The resulting network view provides insights into signaling pathway regulation and reveals a link between the cell-cycle kinase, Cak1, the Fus3 MAP kinase, and a pathway that regulates chromatin integrity during transcription by RNA polymerase II.

Global mapping of herpesvirus-host protein complexes reveals a transcription strategy for late genes.

Mapping host-pathogen interactions has proven instrumental for understanding how viruses manipulate host machinery and how numerous cellular processes are regulated. DNA viruses such as herpesviruses have relatively large coding capacity and thus can target an extensive network of cellular proteins. To identify the host proteins hijacked by this pathogen, we systematically affinity tagged and purified all 89 proteins of Kaposi's sarcoma-associated herpesvirus (KSHV) from human cells. Mass spectrometry of this material identified over 500 virus-host interactions. KSHV causes AIDS-associated cancers, and its interaction network is enriched for proteins linked to cancer and overlaps with proteins that are also targeted by HIV-1. We found that the conserved KSHV protein ORF24 binds to RNA polymerase II and brings it to viral late promoters by mimicking and replacing cellular TATA-box-binding protein (TBP). This is required for herpesviral late gene expression, a complex and poorly understood phase of the viral lifecycle.

A lipid E-MAP identifies Ubx2 as a critical regulator of lipid saturation and lipid bilayer stress.

Biological membranes are complex, and the mechanisms underlying their homeostasis are incompletely understood. Here, we present a quantitative genetic interaction map (E-MAP) focused on various aspects of lipid biology, including lipid metabolism, sorting, and trafficking. This E-MAP contains ∼250,000 negative and positive genetic interaction scores and identifies a molecular crosstalk of protein quality control pathways with lipid bilayer homeostasis. Ubx2p, a component of the endoplasmic-reticulum-associated degradation pathway, surfaces as a key upstream regulator of the essential fatty acid (FA) desaturase Ole1p. Loss of Ubx2p affects the transcriptional control of OLE1, resulting in impaired FA desaturation and a severe shift toward more saturated membrane lipids. Both the induction of the unfolded protein response and aberrant nuclear membrane morphologies observed in cells lacking UBX2 are suppressed by the supplementation of unsaturated FAs. Our results point toward the existence of dedicated bilayer stress responses for membrane homeostasis.

Genetic interaction mapping in mammalian cells using CRISPR interference.

We describe a combinatorial CRISPR interference (CRISPRi) screening platform for mapping genetic interactions in mammalian cells. We targeted 107 chromatin-regulation factors in human cells with pools of either single or double single guide RNAs (sgRNAs) to downregulate individual genes or gene pairs, respectively. Relative enrichment analysis of individual sgRNAs or sgRNA pairs allowed for quantitative characterization of genetic interactions, and comparison with protein-protein-interaction data revealed a functional map of chromatin regulation.

Hierarchical modularity and the evolution of genetic interactomes across species.

To date, cross-species comparisons of genetic interactomes have been restricted to small or functionally related gene sets, limiting our ability to infer evolutionary trends. To facilitate a more comprehensive analysis, we constructed a genome-scale epistasis map (E-MAP) for the fission yeast Schizosaccharomyces pombe, providing phenotypic signatures for ~60% of the nonessential genome. Using these signatures, we generated a catalog of 297 functional modules, and we assigned function to 144 previously uncharacterized genes, including mRNA splicing and DNA damage checkpoint factors. Comparison with an integrated genetic interactome from the budding yeast Saccharomyces cerevisiae revealed a hierarchical model for the evolution of genetic interactions, with conservation highest within protein complexes, lower within biological processes, and lowest between distinct biological processes. Despite the large evolutionary distance and extensive rewiring of individual interactions, both networks retain conserved features and display similar levels of functional crosstalk between biological processes, suggesting general design principles of genetic interactomes.

Genetic interaction mapping reveals a role for the SWI/SNF nucleosome remodeler in spliceosome activation in fission yeast.

Although numerous regulatory connections between pre-mRNA splicing and chromatin have been demonstrated, the precise mechanisms by which chromatin factors influence spliceosome assembly and/or catalysis remain unclear. To probe the genetic network of pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, we constructed an epistatic mini-array profile (E-MAP) and discovered many new connections between chromatin and splicing. Notably, the nucleosome remodeler SWI/SNF had strong genetic interactions with components of the U2 snRNP SF3 complex. Overexpression of SF3 components in ΔSWI/SNF cells led to inefficient splicing of many fission yeast introns, predominantly those with non-consensus splice sites. Deletion of SWI/SNF decreased recruitment of the splicing ATPase Prp2, suggesting that SWI/SNF promotes co-transcriptional spliceosome assembly prior to first step catalysis. Importantly, defects in SWI/SNF as well as SF3 overexpression each altered nucleosome occupancy along intron-containing genes, illustrating that the chromatin landscape both affects--and is affected by--co-transcriptional splicing.

Global landscape of HIV-human protein complexes.

Human immunodeficiency virus (HIV) has a small genome and therefore relies heavily on the host cellular machinery to replicate. Identifying which host proteins and complexes come into physical contact with the viral proteins is crucial for a comprehensive understanding of how HIV rewires the host's cellular machinery during the course of infection. Here we report the use of affinity tagging and purification mass spectrometry to determine systematically the physical interactions of all 18 HIV-1 proteins and polyproteins with host proteins in two different human cell lines (HEK293 and Jurkat). Using a quantitative scoring system that we call MiST, we identified with high confidence 497 HIV-human protein-protein interactions involving 435 individual human proteins, with ∼40% of the interactions being identified in both cell types. We found that the host proteins hijacked by HIV, especially those found interacting in both cell types, are highly conserved across primates. We uncovered a number of host complexes targeted by viral proteins, including the finding that HIV protease cleaves eIF3d, a subunit of eukaryotic translation initiation factor 3. This host protein is one of eleven identified in this analysis that act to inhibit HIV replication. This data set facilitates a more comprehensive and detailed understanding of how the host machinery is manipulated during the course of HIV infection.

Quantitative genetic-interaction mapping in mammalian cells.

Mapping genetic interactions (GIs) by simultaneously perturbing pairs of genes is a powerful tool for understanding complex biological phenomena. Here we describe an experimental platform for generating quantitative GI maps in mammalian cells using a combinatorial RNA interference strategy. We performed ∼11,000 pairwise knockdowns in mouse fibroblasts, focusing on 130 factors involved in chromatin regulation to create a GI map. Comparison of the GI and protein-protein interaction (PPI) data revealed that pairs of genes exhibiting positive GIs and/or similar genetic profiles were predictive of the corresponding proteins being physically associated. The mammalian GI map identified pathways and complexes but also resolved functionally distinct submodules within larger protein complexes. By integrating GI and PPI data, we created a functional map of chromatin complexes in mouse fibroblasts, revealing that the PAF complex is a central player in the mammalian chromatin landscape.

Differential genetic interactions of yeast stress response MAPK pathways.

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

A genetic interaction map of RNA-processing factors reveals links between Sem1/Dss1-containing complexes and mRNA export and splicing.

We used a quantitative, high-density genetic interaction map, or E-MAP (Epistatic MiniArray Profile), to interrogate the relationships within and between RNA-processing pathways. Due to their complexity and the essential roles of many of the components, these pathways have been difficult to functionally dissect. Here, we report the results for 107,155 individual interactions involving 552 mutations, 166 of which are hypomorphic alleles of essential genes. Our data enabled the discovery of links between components of the mRNA export and splicing machineries and Sem1/Dss1, a component of the 19S proteasome. In particular, we demonstrate that Sem1 has a proteasome-independent role in mRNA export as a functional component of the Sac3-Thp1 complex. Sem1 also interacts with Csn12, a component of the COP9 signalosome. Finally, we show that Csn12 plays a role in pre-mRNA splicing, which is independent of other signalosome components. Thus, Sem1 is involved in three separate and functionally distinct complexes.

Genetic interaction analysis of point mutations enables interrogation of gene function at a residue-level resolution: exploring the applications of high-resolution genetic interaction mapping of point mutations.

We have achieved a residue-level resolution of genetic interaction mapping - a technique that measures how the function of one gene is affected by the alteration of a second gene - by analyzing point mutations. Here, we describe how to interpret point mutant genetic interactions, and outline key applications for the approach, including interrogation of protein interaction interfaces and active sites, and examination of post-translational modifications. Genetic interaction analysis has proven effective for characterizing cellular processes; however, to date, systematic high-throughput genetic interaction screens have relied on gene deletions or knockdowns, which limits the resolution of gene function analysis and poses problems for multifunctional genes. Our point mutant approach addresses these issues, and further provides a tool for in vivo structure-function analysis that complements traditional biophysical methods. We also discuss the potential for genetic interaction mapping of point mutations in human cells and its application to personalized medicine.

Functional dissection of protein complexes involved in yeast chromosome biology using a genetic interaction map.

Defining the functional relationships between proteins is critical for understanding virtually all aspects of cell biology. Large-scale identification of protein complexes has provided one important step towards this goal; however, even knowledge of the stoichiometry, affinity and lifetime of every protein-protein interaction would not reveal the functional relationships between and within such complexes. Genetic interactions can provide functional information that is largely invisible to protein-protein interaction data sets. Here we present an epistatic miniarray profile (E-MAP) consisting of quantitative pairwise measurements of the genetic interactions between 743 Saccharomyces cerevisiae genes involved in various aspects of chromosome biology (including DNA replication/repair, chromatid segregation and transcriptional regulation). This E-MAP reveals that physical interactions fall into two well-represented classes distinguished by whether or not the individual proteins act coherently to carry out a common function. Thus, genetic interaction data make it possible to dissect functionally multi-protein complexes, including Mediator, and to organize distinct protein complexes into pathways. In one pathway defined here, we show that Rtt109 is the founding member of a novel class of histone acetyltransferases responsible for Asf1-dependent acetylation of histone H3 on lysine 56. This modification, in turn, enables a ubiquitin ligase complex containing the cullin Rtt101 to ensure genomic integrity during DNA replication.

Global landscape of protein complexes in the yeast Saccharomyces cerevisiae.

Identification of protein-protein interactions often provides insight into protein function, and many cellular processes are performed by stable protein complexes. We used tandem affinity purification to process 4,562 different tagged proteins of the yeast Saccharomyces cerevisiae. Each preparation was analysed by both matrix-assisted laser desorption/ionization-time of flight mass spectrometry and liquid chromatography tandem mass spectrometry to increase coverage and accuracy. Machine learning was used to integrate the mass spectrometry scores and assign probabilities to the protein-protein interactions. Among 4,087 different proteins identified with high confidence by mass spectrometry from 2,357 successful purifications, our core data set (median precision of 0.69) comprises 7,123 protein-protein interactions involving 2,708 proteins. A Markov clustering algorithm organized these interactions into 547 protein complexes averaging 4.9 subunits per complex, about half of them absent from the MIPS database, as well as 429 additional interactions between pairs of complexes. The data (all of which are available online) will help future studies on individual proteins as well as functional genomics and systems biology.

Global post-translational modification profiling of HIV-1-infected cells reveals mechanisms of host cellular pathway remodeling.

Viruses must effectively remodel host cellular pathways to replicate and evade immune defenses, and they must do so with limited genomic coding capacity. Targeting post-translational modification (PTM) pathways provides a mechanism by which viruses can broadly and rapidly transform a hostile host environment into a hospitable one. We use mass spectrometry-based proteomics to quantify changes in protein abundance and two PTM types-phosphorylation and ubiquitination-in response to HIV-1 infection with viruses harboring targeted deletions of a subset of HIV-1 genes. PTM analysis reveals a requirement for Aurora kinase activity in HIV-1 infection and identified putative substrates of a phosphatase that is degraded during infection. Finally, we demonstrate that the HIV-1 Vpr protein inhibits histone H1 ubiquitination, leading to defects in DNA repair.