Host Cell Factors in HIV Replication: Meta-Analysis of Genome-Wide Studies

doi:10.1371/journal.ppat.1000437

PLoS Pathog. 2009 May; 5(5): e1000437.

Published online 2009 May 29. doi: 10.1371/journal.ppat.1000437

PMCID: PMC2682202

Copyright Bushman et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Host Cell Factors in HIV Replication: Meta-Analysis of Genome-Wide Studies

Frederic D. Bushman,¹^* Nirav Malani,¹ Jason Fernandes,^2,³ Iván D'Orso,^2,³ Gerard Cagney,^3,^4,⁵ Tracy L. Diamond,⁶ Honglin Zhou,⁶ Daria J. Hazuda,⁶ Amy S. Espeseth,⁶ Renate König,⁷ Sourav Bandyopadhyay,⁸ Trey Ideker,⁸ Stephen P. Goff,⁹ Nevan J. Krogan,^3,⁴ Alan D. Frankel,^2,³ John A. T. Young,¹⁰ and Sumit K. Chanda⁷^*

¹Department of Microbiology, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania, United States of America

²Department of Biochemistry and Biophysics, University of California, San Francisco, California, United States of America

³HARC Center, San Francisco, California, United States of America

⁴Department of Cellular and Molecular Pharmacology, University of California, San Francisco, California, United States of America

⁵Conway Institute, University College Dublin, Ireland

⁶Department of Antiviral Research, Merck Research Laboratories, West Point, Pennsylvania, United States of America

⁷Infectious & Inflammatory Disease Center, Burnham Institute for Medical Research, La Jolla, California, United States of America

⁸Department of Bioengineering, University of California, San Diego, La Jolla, California, United States of America

⁹Howard Hughes Medical Institute, Department of Biochemistry and Molecular Biophysics, College of Physicians and Surgeons, Columbia University, New York, New York, United States of America

¹⁰Infectious Disease Laboratory, The Salk Institute for Biological Studies, La Jolla, California, United States of America

Glenn F. Rall, Editor

The Fox Chase Cancer Center, United States of America

* E-mail: bushman/at/mail.med.upenn.edu (FDB); Email: schanda/at/burnham.org (SKC)

This article has been cited by other articles in PMC.

Abstract

We have analyzed host cell genes linked to HIV replication that were identified in nine genome-wide studies, including three independent siRNA screens. Overlaps among the siRNA screens were very modest (<7% for any pairwise combination), and similarly, only modest overlaps were seen in pairwise comparisons with other types of genome-wide studies. Combining all genes from the genome-wide studies together with genes reported in the literature to affect HIV yields 2,410 protein-coding genes, or fully 9.5% of all human genes (though of course some of these are false positive calls). Here we report an “encyclopedia” of all overlaps between studies (available at http://www.hostpathogen.org), which yielded a more extensively corroborated set of host factors assisting HIV replication. We used these genes to calculate refined networks that specify cellular subsystems recruited by HIV to assist in replication, and present additional analysis specifying host cell genes that are attractive as potential therapeutic targets.

Introduction

Genome-wide screening technologies offer unprecedented opportunities for discovery [1], but each method is imperfect, so that correct calls will be mixed with false positives, and authentic functions will be missed at some frequency, yielding false negatives. For example, three small interfering RNA (siRNA) screens have been reported that interrogated most of the human genes for effects on HIV infection [2]–[5], but though these screens identified many cellular factors previously implicated in HIV replication, the overlap between any pair of screens was <7%. The siRNA method has many limitations [6],[7]. For a gene to be detected as important during HIV infection, it must meet the following criteria: 1) It must be possible to achieve biologically meaningful reduction in mRNA levels with the siRNA(s) used. 2) The protein must be sufficiently unstable to allow functionally significant reduction over the time course tested. 3) The knockdown must not be toxic. 4) The function targeted cannot be provided by multiple redundant factors. In addition, genes may be called mistakenly due to experimental errors during high throughput analysis or off-target activities of the siRNAs used. Furthermore, siRNAs that do pass all of the above hurdles and affect viral infection may target factors that act only indirectly. Other screening technologies are also fraught with experimental limitations. However, genes identified independently in multiple studies should have a greater chance of being correctly called. Here, we report a meta-analysis of nine genome-wide screens for cellular factors associated with HIV replication.

Genome-Wide Surveys Used in the Meta-Analysis

We analyzed human gene products identified as important for HIV infection in the nine screens presented in Table 1. Three screens (lists 1–3) used transfection of siRNAs to knock down >20,000 human genes, then assessed the efficiency of HIV infection. The König et al. study [3] (list 1) used 293T cells as targets and only examined the steps of uncoating through viral gene expression. This study used a relatively large number of siRNAs per gene and included extensive mapping of the effects of knockdown to steps in the HIV replication cycle. The Brass et al. [2] (list 2) and Zhou et al. [4] (list 3) studies used HeLa cells as targets, and had the advantage of examining all the steps of HIV replication, though with less redundant siRNA coverage. List 4 contains genes near human polymorphisms identified by Fellay et al. as associated with different HIV viral loads in patients [8]. List 5 is composed of genes encoding proteins found in HIV particles that budded out of monocyte-derived macrophages [9]. Lists 6–9 contain raw screening data on binding interactions between HIV proteins and cellular proteins identified using a pull-down mass spectrometry approach (lists 6–8, targeting Nef, Tat, and Rev) or yeast two-hybrid analysis (list 9, targeting IN) [10]. List 10 summarizes interactions between HIV and cellular proteins from the published literature [11]—the depth and quality of the literature is quite variable among these proposed cellular factors, but for this analysis all calls were treated equally. Lists 11 and 12 contain siRNA screen data for two additional viruses (influenza virus studied in fly cells [12], and West Nile Virus studied in human cells [13]), allowing comparison with HIV. Report S1 presents more detailed descriptions of the 12 data sets along with extensive overlap analysis (also available at http://www.hostpathogen.org).

Table 1

Gene Sets Analyzed.

Overview of Genes Proposed to Be Associated with HIV Infection

A total of 1,254 genes were called as important during HIV infection in at least one genome-wide survey (lists 1–9 above), representing about 5% of all human protein-coding genes (using the RefSeq total number of 25,157). One measure of the accuracy of the genome-wide methods is assessing the overlap with genes previously identified in published peer-reviewed studies of HIV (list 10). Comparing the genes called in the HIV interaction database (list 10, 1,434 genes) to those identified in the nine genome-wide surveys (lists 1–9) yielded an overlap of only 257 genes. The union of all genes called in the genome-wide studies and the National Center for Biotechnology Information (NCBI) interaction database (lists 1–10) contains a remarkable 2,393 human protein-coding genes associated with HIV infection, or 9.5% of all human genes.

Were the genes identified in the genome-wide screens (lists 1–9) even enriched at all for previously identified HIV interacting genes (list 10)? The significance of this overlap was assessed by comparison to a random distribution. For the NCBI list of HIV-interacting factors (list 10), 1,434 randomly selected genes were drawn 1,000 times with replacement from the background of all human genes, simulating the NCBI list, and 1,254 random genes were drawn 1,000 times from all human genes as well, simulating the genome-wide list. The overlaps of the 1,000 repetitions of the random draws were quantified and plotted (Figure 1A). The modal number of overlapping genes was 71, and no simulations showed an overlap of 257 or more genes, yielding a highly significant p-value (<0.001). Thus the overlap, though modest, is highly significant. Genes identified in lists 1–10 were analyzed in all pairwise combinations to identify genes in common between each pair (detailed data are in Report S1, pp. 5–70) and the significance tested by simulation or calculated as in [14]. An example is presented in Figure 1B, and the numbers of overlapping genes and their significance for all pairwise combinations of gene lists is shown in Table 2.

Figure 1

Overlap Analysis to Generate an Empirical p-Value.

Table 2

Statistical Analysis of Genes in Common between All Pairs of Genome-Wide Studies.

Analysis of Overlap among the Three siRNA Screens for Genes Affecting HIV Replication

The three siRNA screens (lists 1–3) together called 842 genes as diminishing HIV replication when knocked down, or 3.3% of all human protein-coding genes (Report S1, pp. 98–120). A total of 34 genes were called in at least two siRNA screens (Table 3). Three genes were called in all three screens (MED6, MED7, and RELA). The pairwise overlaps were statistically significant (p<0.024 for all pairs of screens), but the percentages of shared genes were quite modest, ranging from 3% to 6%. The Brass et al. and Zhou et al. screens (lists 2 and 3) both surveyed the entire HIV life cycle and studied infection in HeLa cells, and these two share the greatest overlap (6%). The three siRNA screens identified the NCBI genes as 13.3%–18.3% of the total, indicating highly significant enrichment (p<0.001), as reported previously.

Table 3

Genes Called in at Least Two siRNA Screens.

We then asked whether further enrichment relative to the NCBI HIV interaction database was achieved by examining human genes identified in at least two siRNA two screens. Of the 34 genes on two or more lists, 11 were previously reported in the HIV interaction database (NUP153, CCNT1, CTDP1, CHST1, CD4, CXCR4, TCEB3, JAK1, AKT1, DDX3X, and RELA), or 30% of the total, substantially higher than the 13%–18% identified in each single list alone. From this we infer that the newly identified genes called in two or more siRNA screens (Table 3) are more likely to be authentic new cellular cofactors for HIV infection. Twenty-nine out of the 34 genes were found to be expressed in cells or tissues expressing CD4 and coreceptor by transcriptional profiling analysis, and so competent for HIV entry. Of the remaining five, CCNT1 (cyclinT1) is known to be expressed in T cells and represents a false negative call in the expression data used. A comprehensive table of all genes identified in pairwise combinations of lists 1–12 is provided at the end of Report S1 together with the expression analysis (pp. 72–98).

Why did the three different siRNA screens yield such different gene lists? One possible explanation could be that the expression of host cell factors differed between the HeLa and 293T cells studied. However, analysis of transcriptional profiling data showed that >93% of the genes called as important for HIV infection in any one of the three studies were expressed in both cell types.

However, variation due to 1) experimental noise, 2) timing of sampling, and 3) different filtering criteria likely do explain some of the differences. Two replicates were available for analysis from the König et al. screen, allowing estimation of the variance. From this, the expected overlap for of the top 300 genes in replicate screens could be simulated. A test of two replicates or ten replicates per screen (Figure 1C) yielded 150 or 240 overlapping genes, illustrating how the high variance reduced the overlap, but replication improves it.

A second source of variation was differences between time points analyzed, which varied among the published siRNA screens. Although data were not available for multiple time points for the HIV screens, data were available for a screen of influenza virus infection at three time points (S. Chanda, unpublished data). Analysis demonstrated that variation between time points was of the same magnitude as variation within time points and partially independent.

A third source of variation is likely to be differences in the filtering thresholds used. We investigated the effects of different choices for the toxicity filter by reanalyzing the data of the König et al. screen using three different toxicity thresholds. In the first, no filter was applied (100% of genes were accepted for further analysis), in the second, only genes in the 50% least toxic group were considered, in the third, only the 20% least toxic genes were considered. For each set, the 300 genes with the strongest reduction in HIV infection after knockdown were extracted and the overlap among sets compared (Figure 1D). Fewer than 150 genes out of 300 overlapped between the 100% and 20% sets, and the maximum between any pair was 222 genes, indicating that the final gene set called is very sensitive to the toxicity threshold chosen.

Thus variations between replicates, between time points, and in filtering thresholds all likely contributed to the differences between siRNA screens. Further differences also likely arose from use of different siRNA libraries, cell types, and viral strains [5].

We next asked whether the three siRNA screens yielded host factors participating in similar cellular processes. We extracted overrepresented functional clusters for each screen (lists 1–3) using The Database for Annotation, Visualization, and Integrated Discovery (DAVID) tool [15]. Overrepresented clusters were then filtered for significance (p<0.06 based on a geometric mean of all terms in a group), redundancy, biological relevance, and specificity. This yielded 24 functional groups, most containing contributions from factors identified in two or more screens (Figure 2; genes are cataloged in Table S1). For example, all three screens were enriched in factors involved in “Nuclear Pore/Transport” (21–24 genes each), which likely facilitate the trafficking of HIV complexes between cellular compartments, including the nuclear import of the HIV pre-integration complex, export of viral RNAs, or possibly synthesis of other required factors. Other functional classes or complexes identified included “DNA Repair”, “Ubiquitin-associated”, “Mediator Complex/Transcription”, “RNA Binding”, “GTP Binding”, and “Helicase”. Thus, a functional analysis of the three screens revealed greater overlap in gene ontology (GO) categories than was seen for individual genes.

Figure 2

Gene Ontology Analysis of GO Groups Enriched among Genes Called in Two or More siRNA Screens.

Recovery of Well-Documented Host Genes Affecting HIV in Genome-Wide Studies

A variety of well-documented cellular cofactors for infection were identified in two out of three screens, including i) the binding and entry factors CD4 and CXCR4; ii) the NFkappaB subunit RELA; iii) the activating kinases AKT1 and JAK1; iv) the Vpr and Vif cofactor TCEB3/elonginB; and v) the Tat cofactor CCNT1/cyclinT1 [16] (which was also in the mass spectrometry study of Tat-associated proteins). The Rev cofactor DDX3X [17], an RNA helicase, was also identified in two out of three screens and in addition was found by mass spectrometry to bind to both Tat and Rev. The DNA repair factor MRE11 was also identified, which was previously implicated in HIV DNA circularization [18], though effects on HIV infection efficiency have not been reported previously.

A variety of further well-established factors were identified in one siRNA screen only. The well-studied viral budding factor TSG101 [19] was called in the Zhou et al. siRNA screen (list 3), and also identified as associated with HIV particles after release. The Rev cofactor XPO1/CRM1 [20] was called in the Zhou et al. siRNA screen but not the others.

Also instructive is analyzing the known HIV cofactors that were not identified. HLA-B57 and HLA-C have well-documented effects on viral set point and HIV disease progression [8],[21], but these were not detected in the siRNA screens, probably because the HLA proteins affect the immune response to HIV and not replication at the cellular level. The integration cofactor PSIP1/LEDGF/p75 was not identified, probably because only very complete knockdowns diminish HIV replication [22]–[30]. Several genes known to encode products important for HIV replication were identified in the initial screen of König et al., which yielded 4,019 candidates, but were not further validated in the filtered data set of 293 proteins. These included Sp1, a transcription factor known to bind the HIV LTR; the HIV Gag binding protein cyclophilin A (PPIA) [31]; and several integrin proteins, believed to assist in virus binding to cells (ITGB1, ITGB2, ITGB3). The ESCRT proteins are known to be important in HIV budding [32], but only VPS24 was identified. Another member of this complex, VPS53, was called in the Brass et al. study, and the initial unfiltered König et al. screen, but not in other studies. The RNA lariat debranching enzyme DBR1 was used as a positive control in the König et al. study, and is well known to affect reverse transcription [33], but DBR1 was not identified in any of the other studies. Thus, the recovery of already implicated host factors was generally good in the overlap analysis, providing confidence about the authenticity of the newly called genes. However, some well-documented factors were missed, indicating that other factors important for HIV replication were probably missed by the analysis.

Network Analysis of the HIV–Host Interactome

To identify the cellular subsystems recruited by HIV in more detail, we assembled a host–pathogen protein interaction network based on the gene products in lists 1–10 (Table 1). The network interaction map took advantage of binary protein binding relationships cataloged in curated literature-based protein–protein interaction databases (i.e., Bind, HPRD, MINT, Reactome, etc.). HIV–host interactions were predicted based on evidence compiled in the NCBI HIV interaction database. The resulting HIV–host network was comprised of 1,657 cellular proteins that formed interactions with other host cell factors or HIV-encoded proteins (p<10⁻⁵). Two hundred and ninety of these host proteins (“nodes”) were supported by experimental evidence from two or more datasets, reflecting a 35% enrichment of proteins that are called by multiple datasets in this analysis. We performed a further analysis to identify unusually dense network neighborhoods within this interactome map using a graph theoretic clustering algorithm (MCODE) [34]. This revealed 11 putative molecular clusters, ten of which could be associated with distinct biochemical or cellular functions.

Proteasome

A densely connected network of proteasome subunits was identified by the MCODE analysis (Figure 3A). The proteasome was prominent in the published siRNA screening data, and implicated in probable early steps of viral infection. In previous literature, the proteasome was shown to act negatively on HIV infection by destroying replication intermediates [35],[36]. The siRNA data indicate that the proteasome may also facilitate HIV infection. The mechanism is unclear and could be indirect–for example, reducing proteasome activity may alter cellular ubiquitin levels, and so affect HIV replication by altering the free ubiquitin pool.

Figure 3

Gene Clusters, Generated Using PPI and MCODE Analysis, Derived from the Full Set of Genes Implicated in HIV Infection (Lists 1–10).

Transcription/RNA Polymerase

Genes for subunits of RNA polymerase II and associated factors were identified in several different screens, yielding a densely connected network (Figure 3B). In some of the siRNA screens, the knockdown of Pol II subunits was mapped to the step of Tat transactivation. Because so many subunits were identified, the simplest interpretation is that reduced dosage of the full complex is responsible for the deficit in HIV replication.

Mediator Complex

Multiple subunits of the mediator complex were identified in two or more siRNA screens (Table 3 and Figure 3C). The mediator complex links transcriptional activator proteins to the RNA polymerase II basal transcription apparatus, thereby allowing transcriptional activation [37],[38]. The observation that so many subunits were identified suggests that activity of the complex as a whole is the target of siRNA modulation. Viral replication cycle mapping by Zhou et al. indicated that some of the mediator proteins were needed to support Tat-activated transcription, though studies in König et al. suggest a possible further role in reverse transcription. The data can be accommodated in a model where changes in dosage in the mediator complex are not toxic to cells, but where Tat-activated transcription is extremely sensitive to mediator dosage. Previously, mediator was shown to be important for Sp1-driven transcription, and Sp1 is required for transcription from the HIV LTR, suggesting possible involvement of Sp1 as well.

Tat Activation/Transcriptional Elongation

A dense network was formed containing the Tat cofactor cyclin T1 (CCNT1) [16], which was identified in two out of three siRNA screens and by mass spectrometry (Figure 3D). Together with its binding partner CDK9, which was identified as a Tat binding protein, cyclinT1 forms positive transcription elongation factor b (P-TEFb). The MCODE analysis links the P-TEFb complex and the elongin complex involved in transcriptional elongation. Another factor, the RNA Pol II carboxyl-terminal domain (CTD) phosphatase CTDP1, was also identified and was also previously associated with Tat activation. In addition, two STAT proteins, also involved in transcription and NFkappaB signaling and implicated in lentiviral infection [39], were identified (Figure 3E).

RNA Binding/Splicing

A large cluster of RNA binding and splicing proteins was identified in the MCODE analysis. Eleven of the cellular genes encode protein components of hnRNP complexes (HNR factors) that form on pre-mRNA and direct splicing and other activities. HNRNPU contains both an RNA binding domain and a DNA binding domain that mediates attachment to the nuclear scaffold, potentially linking sites of mRNA synthesis to specific sub-nuclear locations. Six further genes (SF3 factors) encode components of the splicing factor 3 a/b complex, which is involved in activating the U2 snRNP and promoting splicing. Three SNR proteins and two SF proteins were also identified and are implicated in RNA splicing and RNP formation. Several of these proteins were implicated in the literature to modulate Tat or Rev function (e.g., [40]–[42]), and seven direct binding interactions to these viral proteins were identified in the mass spectrometry data reported here.

Several observations also suggest possible connections of RNAP/splicing factors to the viral DNA integration step. Two components of the splicing factor SF3 bound integrase in the yeast two-hybrid data (list 9) [10]. The splicing protein SNW1/SKIIP1 was found by König et al. be selectively important at the integration step [3]. The integrase-interacting protein PSIP1/LEDGF/p75 appears to tether integrase to active transcription units [23],[29],[43],[44], and an alternatively spliced variant of this protein (p52) is involved in RNA metabolism [45]. Though indirect, these observations suggest a model in which splicing factors may help recruit integrase to active transcription units, which are favored for integration [46]–[48].

Another possible role of splicing factors is in maintaining the proper balance between spliced and incompletely spliced HIV RNAs. HIV replication requires multiply spliced messages (encoding Tat, Rev, and Nef), singly spliced messages (encoding Vif, Vpr, Env/Vpu, and a second form of Tat), and unspliced messages (encoding Gag and Gag-Pol). The unspliced RNA also serves as the genomic RNA. Alterations in dosage of splicing factors by siRNA knockdown may well diminish HIV replication by altering the ratios of the different HIV mRNA forms.

The BiP/GRP78/HSPA5 Chaperone

The BiP/GRP78/HSPA5 protein chaperone was identified in the network analysis (Figure 3G). BiP/GRP78/HSPA5 is a member of the HSP70 family that is involved in the folding and assembly of proteins in the endoplasmic reticulum. BiP has been implicated in interacting with newly synthesized HIV gp160 SU/TM precursor [49], and HSP70 family members have been proposed to interact with Gag, Tat, Vpr, and MA. The MCODE analysis connected BiP/GRP78/HSPA5 to a collection of nuclear proteins involved in splicing (PRPF8, SFPQ, and SNW1), nuclear matrix architecture (MATR3), and ubiquitylation (UBR5). Determining how these cellular proteins modulate the interactions of BiP/GRP78/HSPA5 with HIV proteins offers a potential route to better understanding protein folding and sorting during HIV replication.

The CCT Chaperone

The MCODE analysis identified subunits of the chaperone containing TCP1 (CCT) complex (Figure 3H). Subunits were identified in siRNA screens, in HIV particles after budding, and also as Tat binding proteins. This complex consists of two identical stacked rings of eight subunits. Unfolded proteins are thought to pass through the central cavity, and become folded in an ATP-dependent manner. The CCT chaperone has not previously been associated with HIV replication, and represents a new candidate for involvement in Tat activation and HIV budding.

Additional Densely Connected Clusters

Several further functions were identified, including proteins involved in t-RNA synthase function, transport, and one of unknown function (Figure 3 I–3K). The t-RNA synthase and transport complexes contained members that associated with Tat according to the mass spectrometry study, and the unknown complex contained a member binding to Rev, suggesting specific links to HIV replication.

Other Newly Identified Functions

Several sets of proteins were identified that were not called as densely connected networks but appear to be functionally related. The nuclear pore and associated factors were clustered in the initial MCODE network but were not sufficiently densely to emerge as a densely connected network. Proteins involved in nuclear import identified in two out of three siRNA screens included products of NUP153, RANBP2, TNPO3, and RGPD8. NUP153 and TNPO3 have been associated with the trafficking of HIV proteins previously [2], [3], [50]–[52]. RANBP2 is a giant gene encoding a product that accumulates at nuclear pores and binds to RAN, which is a small GTP-binding protein of the RAS superfamily. RANBP2 also contains FG repeats, a cyclophilin-related nucleoporin, and a domain that binds UBC9, the E2 for SUMO1 transfer. RGPD8 is named for “RANBP2-like and GRIP domain containing 8”. It too accumulates at the nuclear pore and is believed to assist in RNA and protein transport. The actions of NUP153 and TNPO3 have been mapped to nuclear import of the HIV preintegration complex in [2],[3],[51], and NUP153 has also been proposed to be involved in export of HIV Rev [53].

Three genes were identified that affect the microtubule system. MAP4 is a microtubule-associated protein that has not previously been studied in detail. MID1IP1 is a regulator of microtubule polymerization. CAV2 (caveolin 2) is involved in the formation of plasma membrane invaginations involved in a variety of cellular functions including signal transduction, cell growth, and apoptosis. Caveoli have also been implicated as interacting with the microtubule network [54]. Previous studies have suggested that HIV particles may traffic along microtubules to reach the nucleus [55]. Thus MAP4, MID1IP1, and possibly CAV2 are candidates for cofactors in this process. Other proteins were also called in two siRNA screens (ANAPC2, DMXL1, HMCN2, and IDH1) but are of unknown function (Table 3).

Identifying New Drug Targets

One of the main reasons for carrying out the screens for host factors is the hope of identifying new targets for HIV therapeutics. Several studies have indentified potentially “druggable” human proteins by cataloging families of InterPro domains where one member is the target of one or more small molecule inhibitors with drug-like properties. All members of the family are then proposed as potential drug targets [56],[57] (John Hogenesch, data available at http://www1.qiagen.com/Products/GeneSilencing/LibrarySiRna/SiRnaSets/HumanDruggableGenomesiRNASetV30.aspx?ShowInfo=1). In Report S1 (pp. 72–120), we annotate our overlap study for “druggable” targets by these criteria. Focusing on an updated version of the list from Hopkins and Groom [56], we found that eight of the 34 genes common to the two siRNA screens were called as potential drug targets (Table 3, column labeled “druggable”). Two of these eight are in fact known to be the targets of small molecules with activity against HIV. CXCR4 is the target of AMD3100 and related molecules [58],[59], and AKT1 is the target of miltefosine [60]. The fact that two out of eight genes called as druggable are known antiviral drug targets (at least in tissue culture) suggests that this analysis is yielding viable new targets. Of the additional genes encoding candidate drug targets, inhibitors have been reported for the kinases ADRBK1/GRK2 and JAK1. These can be tested for activity against HIV in cell culture. Annotation of the larger collection of genes that were found on two or more lists (lists 1–9) yielded a further 56 genes encoding potentially druggable cellular factors.

Summary

Analysis of genes called as important for HIV replication in multiple genome-wide screens yielded a list rich in well-known factors and also intriguing new candidates. Many important factors were surely missed by this approach, but at least some of the most promising new genes can be distilled from among the 9.5% of all human protein-coding genes now proposed to affect HIV infection. Many of the new genes can be linked into clusters, specifying cellular subsystems associated with HIV replication. Promising drug targets could be discerned among the best-documented new factors.

Methods

Overlap analysis and comparisons to random distributions were carried out using R [61]. The p-values for overlaps between lists were generated by comparison to results of random simulation and by calculation based on the hypergeometric distribution as in [14]. No correction was applied for multiple comparisons. Networks were generated using MCODE analysis on the binary interaction file and plotted in cytoscape. The global interaction network was judged to be statistically significant by comparison to random simulation (p<10⁻⁵). Subnetworks were selected that 1) contained at least two proteins from different studies and 2) showed high connectivity. The protein–protein interaction data for Nef, Tat, and Rev from the HARC Center were derived using previously described methods (LC MS-MS followed by database matching) [62]. A more thorough analysis of these (early stage) data will be reported elsewhere.

Updated versions of Report S1 documenting the overlap among genome-wide screens can be found at http://www.hostpathogen.org.

Supporting Information

Report S1

Gene List Comparison Report

(1.88 MB PDF)

Click here for additional data file.^{(1.7M, pdf)}

Table S1

Gene IDs from Functional Clusters in Figure 3

(0.09 MB XLS)

Click here for additional data file.^{(88K, xls)}

Acknowledgments

We particularly thank Yingyao Zhou for help with the statistical analysis of overlap among the siRNA screens. We thank Jenny Gu, Jiewei Xu, Marie Fahey, and Melanie Brewer for technical assistance with mass spectrometry and analyses.

Footnotes

The authors have declared that no competing interests exist.

This work was supported by NIH (http://www.nih.gov) grant AI52845 to FDB and the University of Pennsylvania Center for AIDS Research, by NIH grant AI72645 to JATY and SKC, and by NIH HARC Center grant P50 GM082250. ID is an amfAR (http://www.amfar.org) Mathilde Krim Fellow in Basic Biomedical Research (106988-43-RFNT). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Rines DR, Tu B, Miraglia L, Welch GL, Zhang J, et al. High-content screening of functional genomic libraries. Methods Enzymol. 2006;414:530–565. [PubMed]

2. Brass AL, Dykxhoorn DM, Benita Y, Yan N, Engelman A, et al. Identification of host proteins required for HIV infection through a functional genomic screen. Science. 2008;319(5865):921–926. doi:10.1126/science.1152725. [PubMed]

3. Konig R, Zhou Y, Elleder D, Diamond TL, Bonamy GM, et al. Global analysis of host-pathogen interactions that regulate early-stage HIV-1 replication. Cell. 2008;135(1):49–60. [PMC free article] [PubMed]

4. Zhou H, Xu M, Huang Q, Gates AT, Zhang XD, et al. Genome-scale RNAi screen for host factors required for HIV replication. Cell Host Microbe. 2008;4(5):495–504. doi:10.1016/j.chom.2008.10.004. [PubMed]

5. Goff SP. Knockdown screens to knockout HIV-1. Cell. 2008;135(3):417–420. [PubMed]

6. Hutvagner G, Zamore PD. RNAi: Nature abhors a double-strand. Curr Opin Genet Dev. 2002;12(2):225–232. [PubMed]

7. Konig R, Chiang CY, Tu BP, Yan SF, DeJesus PD, et al. A probability-based approach for the analysis of large-scale RNAi screens. Nat Methods. 2007;4(10):847–849. doi:10.1038/nmeth1089. [PubMed]

8. Fellay J, Shianna KV, Ge D, Colombo S, Ledergerber B, et al. A whole-genome association study of major determinants for host control of HIV-1. Science. 2007;317(5840):944–947. doi:10.1126/science.1143767. [PMC free article] [PubMed]

9. Chertova E, Chertov O, Coren LV, Roser JD, Trubey CM, et al. Proteomic and biochemical analysis of purified human immunodeficiency virus type 1 produced from infected monocyte-derived macrophages. J Virol. 2006;80(18):9039–9052. [PMC free article] [PubMed]

10. Studamire B, Goff SP. Host proteins interacting with the moloney murine leukemia virus integrase: Multiple transcriptional regulators and chromatin binding factors. Retrovirology. 2008;5:48. doi:10.1186/1742-4690-5-48. [PMC free article] [PubMed]

11. Fu W, Sanders-Beer BE, Katz KS, Maglott DR, Pruitt KD, et al. Human immunodeficiency virus type 1, human protein interaction database at NCBI. Nucleic Acids Res. 2009;37(Database issue):D417–D422. Epub 15 Oct 2008. doi:10.1093/nar/gkn708. [PMC free article] [PubMed]

12. Hao L, Sakurai A, Watanabe T, Sorensen E, Nidom CA, et al. Drosophila RNAi screen identifies host genes important for influenza virus replication. Nature. 2008;454(7206):890–893. doi:10.1038/nature07151. [PMC free article] [PubMed]

13. Krishnan MN, Ng A, Sukumaran B, Gilfoy FD, Uchil PD, et al. RNA interference screen for human genes associated with west nile virus infection. Nature. 2008;455(7210):242–245. [PMC free article] [PubMed]

14. Fury W, Batliwalla F, Gregersen PK, Li W. Overlapping probabilities of top ranking gene lists, hypergeometric distribution, and stringency of gene selection criterion. Conf Proc IEEE Eng Med Biol Soc. 2006;1:5531–5534. [PubMed]

15. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, et al. Database for annotation, visualization, and integrated discovery. Genome Bio. 2003;4:P3. [PubMed]

16. Wei P, Garber ME, Fang SM, Fischer WH, Jones KA. A novel CDK9-associated C-type cyclin interacts directly with HIV-1 tat and mediates its high-affinity, loop-specific binding to TAR RNA. Cell. 1998;92:451–462. [PubMed]

17. Yedavalli VS, Neuveut C, Chi YH, Kleiman L, Jeang KT. Requirement of DDX3 DEAD box RNA helicase for HIV-1 rev-RRE export function. Cell. 2004;119(3):381–392. [PubMed]

18. Kilzer JM, Stracker TH, Beitzel B, Meek K, Weitzman MD, et al. Roles of host cell factors in circularization of retroviral DNA. Virology. 2003;314:460–467. [PubMed]

19. Garrus JE, von Schwedler UK, Pornillos OW, Morham SG, Zavitz KH, et al. Tsg101 and the vacuolar protein sorting pathway are essential for hiv-1 budding. Cell. 2001;107(1):55–65. [PubMed]

20. Neville M, Stutz F, Lee L, Davis LI, Rosbash M. The importin-beta family member Crm1p bridges the interaction between rev and the nuclear pore complex during nuclear export. Curr Biol. 1997;7(10):767–775. [PubMed]

21. Matthews PC, Prendergast A, Leslie A, Crawford H, Payne R, et al. Central role of reverting mutations in HLA associations with human immunodeficiency virus set point. J Virol. 2008;82(17):8548–8559. doi:10.1128/JVI.00580-08. [PMC free article] [PubMed]

22. Llano M, Saenz DT, Meehan A, Wongthida P, Peretz M, et al. An essential role for LEDGF/p75 in HIV integration. Science. 2006;314(5798):461–464. [PubMed]

23. Marshall H, Ronen K, Berry C, Llano M, Sutherland H, et al. Role of PSIP1/LEDGF/p75 in lentiviral infectivity and integration targeting. PLoS ONE. 2007;2(12):e1340. doi:10.1371/journal.pone.0001340. [PMC free article] [PubMed]

24. Ciuffi A, Diamond T, Hwang Y, Marshall H, Bushman FD. Fusions of LEDGF/p75 to lambda repressor promote HIV DNA integration near lambda operators in vitro. Hum Gene Ther. 2006;17:960–967. [PubMed]

25. Llano M, Vanegas M, Fregoso O, Saenz D, Chung S, et al. LEDGF/p75 determines cellular trafficking of diverse lentiviral but not murine oncoretroviral integrase proteins and is a component of functional lentiviral preintegration complexes. J Virol. 2004;78(17):9524–9537. [PMC free article] [PubMed]

26. Cherepanov P. LEDGF/p75 interacts with divergent lentiviral integrases and modulates their enzymatic activity in vitro. Nucleic Acids Res. 2007;35(1):113–124. [PMC free article] [PubMed]

27. Cherepanov P, Maertens G, Proost P, Devreese B, Van Beeumen J, et al. HIV-1 integrase forms stable tetramers and associates with LEDGF/p75 protein in human cells. J Biol Chem. 2003;278:372–381. [PubMed]

28. Maertens G, Cherepanov P, Pluymers W, Busschots K, De Clercq E, et al. LEDGF/p75 is essential for nuclear and chromosomal targeting of HIV-1 integrase in human cells. J Biol Chem. 2003;278:33528–33539. [PubMed]

29. Ciuffi A, Llano M, Poeschla E, Hoffmann C, Leipzig J, et al. A role for LEDGF/p75 in targeting HIV DNA integration. Nat Med. 2005;11(12):1287–1289. [PubMed]

30. Turlure F, Devroe E, Silver PA, Engelman A. Human cell proteins and human immunodeficiency virus DNA integration. Front Biosci. 2004;9:3187–3208. [PubMed]

31. Gamble TR, Vajdos FF, Yoo S, Worthylake DK, Houseweart M, et al. Crystal structure of human cyclophillin A bound to the amino-terminal domain of HIV-1 capsid. Cell. 1996;87:1285–1294. [PubMed]

32. von Schwedler UK, Stuchell M, Muller B, Ward DM, Chung HY, et al. The protein network of HIV budding. Cell. 2003;114(6):701–713. [PubMed]

33. Ye Y, De Leon J, Yokoyama N, Naidu Y, Camerini D. DBR1 siRNA inhibition of HIV-1 replication. Retrovirology. 2005;2:63. doi:10.1186/1742-4690-2-63. [PMC free article] [PubMed]

34. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003;4:2. [PMC free article] [PubMed]

35. Schwartz O, V M, Friguet B, Arenzana-Seisdedos F, Heard J. Antiviral activity of the proteasome on incoming human immunodeficiency virus type 1. J Virol. 1998;72:3845–3850. [PMC free article] [PubMed]

36. Butler SL, Johnson EP, Bushman FD. HIV cDNA metabolism studied by fluorescence-monitored PCR: Notable stability of two-LTR circles. J Virol. 2002;76:3739–3747. [PMC free article] [PubMed]

37. Kuras L, Borggrefe T, Kornberg RD. Association of the mediator complex with enhancers of active genes. Proc Natl Acad Sci U S A. 2003;100(24):13887–13891. doi:10.1073/pnas.2036346100. [PubMed]

38. Kornberg RD. The molecular basis of eukaryotic transcription. Proc Natl Acad Sci U S A. 2007;104(32):12955–12961. doi:10.1073/pnas.0704138104. [PubMed]

39. Mohan M, Aye PP, Borda JT, Alvarez XJ, Lackner AA. Gastrointestinal disease in SIV-infected rhesus macaques is characterized by proinflammatory dysregulation of the IL-6-JAK-STAT3 pathway. Am J Pathol. 2007;171(6):1952–1965. [PubMed]

40. Pongoski J, Asai K, Cochrane A. Positive and negative modulation of human immunodeficiency virus type 1 rev function by cis and trans regulators of viral RNA splicing. J Virol. 2002;76(10):5108–5120. [PMC free article] [PubMed]

41. Roy BB, Hu J, Guo X, Russell RS, Guo F, et al. Association of RNA helicase a with human immunodeficiency virus type 1 particles. J Biol Chem. 2006;281(18):12625–12635. doi:10.1074/jbc.M510596200. [PubMed]

42. Li J, Tang H, Mullen TM, Westberg C, Reddy TR, et al. A role for RNA helicase A in post-transcriptional regulation of HIV type 1. Proc Natl Acad Sci U S A. 1999;96(2):709–714. [PubMed]

43. Shun MC, Raghavendra NK, Vandegraaff N, Daigle JE, Hughes S, et al. LEDGF/p75 functions downstream from preintegration complex formation to effect gene-specific HIV-1 integration. Genes Dev. 2007;21(14):1767–1778. [PubMed]

44. Ciuffi A, Diamond TL, Hwang Y, Marshall HM, Bushman FD. Modulating target site selection during human immunodeficiency virus DNA integration in vitro with an engineered tethering factor. Hum Gene Ther. 2006;17(9):960–967. [PubMed]

45. Ge H, Si Y, Roeder RG. Isolation of cDNAs encoding novel transcription coactivators p52 and p75 reveals an alternate regulatory mechanism of transcriptional activation. EMBO J. 1998;17:6723–6729. [PubMed]

46. Schroder AR, Shinn P, Chen H, Berry C, Ecker JR, et al. HIV-1 integration in the human genome favors active genes and local hotspots. Cell. 2002;110(4):521–529. [PubMed]

47. Mitchell RS, Beitzel BF, Schroder AR, Shinn P, Chen H, et al. Retroviral DNA integration: ASLV, HIV, and MLV show distinct target site preferences. PLoS Biol. 2004;2(8):e234. doi:10.1371/journal.pbio.0020234. [PMC free article] [PubMed]

48. Wu X, Li Y, Crise B, Burgess SM. Transcription start regions in the human genome are favored targets for MLV integration. Science. 2003;300(5626):1749–1751. [PubMed]

49. Earl PL, Moss B, Doms RW. Folding, interaction with GRP78-BiP, assembly, and transport of the human immunodeficiency virus type 1 envelope protein. J Virol. 1991;65(4):2047–2055. [PMC free article] [PubMed]

50. Luban J. HIV-1 infection: Going nuclear with TNPO3/Transportin-SR2 and integrase. Curr Biol. 2008;18(16):R710–R713. [PubMed]

51. Christ F, Thys W, De Rijck J, Gijsbers R, Albanese A, et al. Transportin-SR2 imports HIV into the nucleus. Curr Biol. 2008;18(16):1192–1202. [PubMed]

52. Varadarajan P, Mahalingam S, Liu P, Ng SB, Gandotra S, et al. The functionally conserved nucleoporins Nup124p from fission yeast and the human Nup153 mediate nuclear import and activity of the Tf1 retrotransposon and HIV-1 vpr. Mol Biol Cell. 2005;16(4):1823–1838. doi:10.1091/mbc.E04-07-0583. [PMC free article] [PubMed]

53. Zolotukhin AS, Felber BK. Nucleoporins nup98 and nup214 participate in nuclear export of human immunodeficiency virus type 1 rev. J Virol. 1999;73(1):120–127. [PMC free article] [PubMed]

54. Ito J, Kheirollah A, Nagayasu Y, Lu R, Kato K, et al. Apolipoprotein A-I increases association of cytosolic cholesterol and caveolin-1 with microtubule cytoskeletons in rat astrocytes. J Neurochem. 2006;97(4):1034–1043. doi:10.1111/j.1471-4159.2006.03805.x. [PubMed]

55. McDonald D, Vodicka MA, Lucero G, Svitkina TM, Borisy GG, et al. Visualization of the intracellular behavior of HIV in living cells. J Cell Biol. 2002;159(3):441–52. [PMC free article] [PubMed]

56. Hopkins AL, Groom CR. The druggable genome. Nat Rev Drug Discov. 2002;1(9):727–730. [PubMed]

57. Orth AP, Batalov S, Perrone M, Chanda SK. The promise of genomics to identify novel therapeutic targets. Expert Opin Ther Targets. 2004;8(6):587–596. [PubMed]

58. Schols D, Este JA, Henson G, De Clercq E. Bicyclams, a class of potent anti-HIV agents, are targeted at the HIV coreceptor fusin/CXCR-4. Antiviral Res. 1997;35(3):147–156. [PubMed]

59. Harrison JE, Lynch JB, Sierra LJ, Blackburn LA, Ray N, et al. Baseline resistance of primary HIV-1 strains to the CXCR4 inhibitor AMD3100. J Virol. 2008;82(23):11695–11704. doi:10.1128/JVI.01303-08. [PMC free article] [PubMed]

60. Chugh P, Bradel-Tretheway B, Monteiro-Filho CM, Planelles V, Maggirwar SB, et al. Akt inhibitors as an HIV-1 infected macrophage-specific anti-viral therapy. Retrovirology. 2008;5:11. doi:10.1186/1742-4690-5-11. [PMC free article] [PubMed]

61. R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing 2006

62. Krogan NJ, Cagney G, Yu H, Zhong G, Guo X, et al. Global landscape of protein complexes in the yeast saccharomyces cerevisiae. Nature. 2006;440(7084):637–643. doi:10.1038/nature04670. [PubMed]

Articles from PLoS Pathogens are provided here courtesy of
Public Library of Science