key: cord-0741812-avnf6w3m
authors: Montano, M; Rarick, M; Sebastiani, P; Brinkmann, P; Russell, M; Navis, A; Wester, C; Thior, I; Essex, M
title: Gene-expression profiling of HIV-1 infection and perinatal transmission in Botswana
date: 2006-05-04
journal: Genes Immun
DOI: 10.1038/sj.gene.6364297
sha: 0d28bf4dd63a40ee27a56df7c573e556159cff56
doc_id: 741812
cord_uid: avnf6w3m

Perinatal transmission of human immunodeficiency virus (HIV)-1 represents a major problem in many regions of the world, especially Southern Africa. With the exception of viral and proviral load, the role for maternal cofactors in perinatal transmission outcome is largely unknown. In this study, an assessment was made of peripheral blood mononuclear cells (PBMC) gene-expression profiles to better understand transcriptional changes associated with HIV-1 infection and perinatal transmission among young adult mothers with infants in Botswana. Peripheral blood mononuclear cells specimens were used from 25 HIV+ drug naive and 20 HIV− healthy mothers, similar in age and location, collected in 1999–2000 and 2003, and processed with the exact same methods, as previously described. Expression profiling of 22 277 microarray gene probes implicated a broad initiation of innate response gene-sets, including toll-like receptor, interferon-stimulated and antiviral RNA response pathways in association with maternal HIV-1 infection. Maternal transmission status was further associated with host genes that influence RNA processing and splicing patterns. In addition to real-time polymerase chain reaction validation of specific genes, enriched category validation of PBMC profiles was conducted using two independent data sets for either HIV-1 infection or an unrelated RNA virus, severe acute respiratory virus infection. HIV-1 pathogen-specific host profiles should prove a useful tool in infection and transmission intervention efforts worldwide. SUPPLEMENTARY INFORMATION: The online version of this article (doi:10.1038/sj.gene.6364297) contains supplementary material, which is available to authorized users.

Various immunologic and virologic factors influence perinatal transmission; but their relative contributions are difficult to assess, suggesting that transmission is multifactorial. [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] Nevertheless, the use of antiretroviral drugs, based on early results from the PACTG 076 Study, 15 has substantially reduced the number of perinatally acquired AIDS cases in the US. A significant challenge, however, still remains within Africa and other resource-limited locales, wherein approximately 700 000 new infections are estimated to have occurred among children in 2003 (WHO). Southern Africa has the highest prevalence of human immunodeficiency virus (HIV)-1 infection in the world (http://www.unaids.org). The predominant HIV in this region is subtype C (HIV-1C). 16 A serosurveillance study conducted in 2003 among pregnant women in Botswana indicated that 37.4% were infected and that among women aged [25] [26] [27] [28] [29] 49 .7% were infected. 17 The role of maternal viral load as a strong predictor of perinatal transmission outcome has been well established, although we and others have observed that there is a substantial overlap in the detectable range of plasma and genital fluid associated viral load between those who transmit virus and those who do not. [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] In many cases, a clear threshold of virus has not been identified above, which perinatal transmission is 100% nor has a threshold been consistently identified below which transmission does not occur. 19, 27, 32, 34 Relatively few studies have examined the role of maternal genetic factors on transmission, and have generally been limited to coreceptor/ligand polymorphism. 11, 35, 36 Notably, previous studies have implicated a role for the chemokine/coreceptors SDF1/CXCR4 and RANTES/CCR5 in HIV-1 disease progression, and polymorphisms in SDF1 and CCR5 have been associated with perinatal transmission. [36] [37] [38] Also, cellular immune response factors have been previously hypothesized to influence perinatal transmission, specifically interleukin (IL)-10 and interferon g (IFN-g). 39 These genetic studies to date have tended to focus on a small number of candidate genes and have not taken advantage of the genomic survey approach that is possible with the use of gene-expression profiling. Such profiles have been used successfully to investigate HIV infection in US subjects. 40, 41 However, no such data exist among African subjects, which account for most infections worldwide. Nor is it known whether gene-expression profiles are associated with transmission outcome.

Human immunodeficiency virus infection has been previously associated with a complex set of host-virus interactions that variably contribute to influence overall host response. 42, 43 The considerable genetic variation present within the human population allows for the possibility of differential host effects on viral replication and immune response. [44] [45] [46] To date, many polymorphic genes predominantly within US-based populations have been described that influence HIV disease progression and infection (termed ARGs: AIDS restriction genes) and they include CCR5, [47] [48] [49] CCR2, 50,51 CCL5/RANTES, 52 CXCL12/SDF-1, 53 CXCR6, 54 CCL2/MCP-1, CCL7/ MCP-3, CCL11/Eotaxin, 55 IL-10, 56 IFN-g, 57 HLA 58-60 and KIR3DS1, 61 and more recently, the RNA editing gene APOBEC3G. 62 (see review by O'Brien and Nelson 63 ).

The present study was undertaken to survey a broad array of gene probes (22 277) , to identify genes that were differentially expressed during HIV-1C infection and to determine whether the expression differences could also be associated with perinatal transmission outcome. We examined the relation between gene-expression profiles and mother-to-infant transmission of HIV-1C among women and infants that were cross-sectionally identified within Botswana.

Sampling for microarray analysis was representative of initial population Table 1 compares the viral load statistics for the specimen groups used for microarray analysis with the previously described maternal HIV-1 infection and transmitter status larger data set. 18 The groups tested had similar viral loads, and were similar in age and nationality (all citizens of Botswana). The specimens consisted of HIV þ mothers (n ¼ 25), including a subset of transmitter (TRs) mothers, TRs (n ¼ 11), nontransmitter (NTRs) mothers, NTRs (n ¼ 14) and an HIVÀ control population of mothers (n ¼ 20). TRs, NTRs and seronegative controls did not differ significantly with regard to location, age, clinical status or condition, parity, Cesarean section experience, breast-feeding practice, or prevalence of sexually transmitted diseases (Montano et al 18 

A statistical method implemented in the program BADGE 64-66 was used to identify differentially expressed genes associated with infection or maternal transmitter status in four grouped comparisons: HIVÀ vs HIV þ (group 1); HIVÀ vs TR (group 2); HIVÀ vs NTR (group 3) and HIV NTR vs TR (group 4). Differentially expressed genes within these groups were then subjected to gene category over-representation analysis (termed 'enrichment') using a standalone version of EASE 67 that contains annotated gene-sets available in the Gene Onotology (GO) database (see Methods) and additional literature query based gene-sets, as described. 68 Enrichment of gene categories was further annotated into biological themes and plotted for each comparison group based on the enrichment significance score (Figure 1 , note that all categories shown were significant, that is, below P ¼ 0.05). See also Supplementary Figures 1 and 2 for the entire list of significant categories with associated dendrograms (http://idisk.mac.com/monty and alan-Public/HIV-neg-pos-analvsis/web-supplement/index. html).

Expression profiles for maternal infection with HIV-1C differed from the seronegative reference subjects (group 1 comparison) in enriched gene-sets specifically associated with IMMUNE RESPONSE (P ¼ 1.38E-09) and enriched gene-sets specifically associated with RNA (mRNA metabolism including processing and editing, P ¼ 3.46E-06). The immune response categories included many genes that overlapped with other significant categories including ANTIVIRAL (P ¼ 0.0001) and INTERFERON (P ¼ 0.0038). The profile of biologically enriched categories also differed significantly between TR and NTR mothers compared with HIVÀ controls. For example, the enriched RNA categories were predominant in HIV-1C TR mothers and markedly distinguished this group of subjects from the control subjects (group 2 comparison, P ¼ 4.70E-06). By contrast, the enriched IMMUNE RESPONSE categories were predominant in HIV-1C NTR subjects in comparison with controls (group 3 comparison, P ¼ 1.41E-12). The comparison of HIV-1C TRs with HIV-1C NTRs (group 4 comparison) resembled the comparison of HIV-1C TRs with HIVÀ, although the magnitude of significance for RNA categories was reduced, possibly due to sample size (compare P ¼ 0.0003 for group 4 with P ¼ 3.46E-06 for group 1 comparison). As described in the Methods section, estimates of significance were based on two steps: the first was to determine the probability for differential expression in each group comparison. The second step was then to evaluate the 'enrichment' for differentially expressed gene-sets (biological categories) using a modified Fisher's exact test. The results are shown graphically in Figure 1 , for categories with significance greater than P ¼ 0.05 in groups 1-4 comparisons.

Human immunodeficiency virus-1C infection increased expression of immune response genes associated with Toll-like receptor and interferon g pathways and RNA processing genes associated with interferon response Based on the significance score for biological categories identified in each comparison group shown in Figure 1 , and to visualize trends in gene expression, profiles for all differentially expressed gene-sets were systematically converted into heatmaps, and representative gene expression data from selected heatmaps were displayed as boxplots (for complete set of heatmaps for each comparison group, see Supplementary Transmitter mothers displayed a broad reduction in RNA processing genes, except for antiviral RNA associated genes The most prominent categories identified in the group 2 comparison (TR vs HIVÀ control subjects) were RNA associated gene-sets (e.g. mRNA metabolism, P ¼ 1.36E-07) with most genes downregulated. Similarly, most genes in this category were also downregulated in NTR subjects, compared with the HIVÀ subjects. A notable exception within this categorical comparison was a subset of genes upregulated in association with RNA binding activity (P ¼ 2.3E-09) and interferon induced antiviral response. Interferon induced RNA response genes included OAS1-3, OASL, ADAR and MX1. The OAS genes encode essential proteins involved in the innate immune response to viral infection. These molecules activate latent RNase L, which results in viral RNA degradation and the inhibition of viral replication. These specific genes tended to be upregulated in NTR subjects compared with HIVÀ controls (see Supplementary Figure 3 , heatmap for group 3, RNA binding category http://idisk.mac.com/monty and alan-Public/ HIV-neg-pos-analysis/web-supplement/index.html), but as a category did not reach significance in the group 3 comparison -due, in part, to the presence of fewer genes within each category and the low expression trend difference in comparison with HIVÀ control subjects.

Nontransmitter mothers displayed a more robust immune response profile than transmitter mothers, particularly in genes associated with antiviral and interferon activity The most prominent biological categories identified in the group 3 comparison (NTR vs HIVÀ subjects) were IMMUNE RESPONSE (P ¼ 2.78E-10), and ANTIVIRAL (P ¼ 5.28E-06), with most genes upregulated. There was a notable absence of significant RNA processing categories in the group 3 comparison, due to a lack of sufficient differential expression among RNA gene-sets in the seronegative controls compared with NTRs (this was in contrast with the group 2 (TR vs HIVÀ) comparison, see Figure 1 ). The ANTIVIRAL genes induced included CCL5, and several interferon induced antiviral RNA response genes including MX1, OAS1-3, IFI35, IFI44, PRKR. Representative group 3 upregulated genes in IMMUNE RESPONSE (P ¼ 2.78E-10) included innate response and RNA editing genes such as IRF7, CCR5, MYD88, ADAR and CCL4. The tendency for most genes to be upregulated in this category (group 3) accounted for the higher significance, in contrast with the group 2 comparison (TR vs HIVÀ) with IMMUNE RESPONSE significance lower at P ¼ 6.16E-06.

Nontransmitter mothers displayed altered expression of RNA processing and splicing associated genes and displayed two expression subclusters that differed by viral load In addition to evaluating NTR and TR profiles to a negative control population, they were compared directly to each other. Although the sample size was small (14 NTR vs 11 TR), significant categories were also identified in the TR vs NTR (group 4) comparison, implicating RNA associated gene-sets representing RNA processing and splicing. Although the gene-sets clustered together, we noted that the expression profile for the NTR subjects (when compared to TR subjects) contained two subsets of specimens. Figure 3a shows the NTR subsets, termed NTR-hi and NTR-lo, which were correlated with expression trends in a heatmap for the category RNA PROCESSING. The subsets differed in their viral load range, with one subset exhibiting a relatively lower viral load (NTR-lo, mean VL ¼ 3.88 logs) and the other subset exhibiting a significantly higher viral load (NTR-hi, mean VL ¼ 4.84 logs) (Figure 3b ). The NTR-hi expression profile resembled levels present in TRs (mean VL ¼ 4.70) in many (but not all) gene-sets evaluated in the group 4 comparisons. Therefore, gene expression levels for some categories within NTRs appeared to differ in association with viral burden and these categories tended to be related to RNA processing functions. In contrast with NTRs, the TR subset of specimens did not appear to contain gene expression subsets that differed significantly based on viral burden (data not shown).

The observation of NTR subsets associated with viral load prompted us to directly evaluate gene-sets that displayed differential expression in association with viral load. Significant enrichment for gene-sets was observed for both groups 1 and 3 comparisons. For the group 1 comparison (HIV-1 þ vs HIVÀ), representative categories associated with viral load were generally consistent with broad innate response (see Supplementary Figure 4 for all categories with P-value less than 0.05 and associated heatmaps http://idisk.mac.com/monty and alan-Public/HIV-neg-pos-analysis/web-supplement/index. html) and included gene-sets associated with interferon signaling (e.g. HEMATOPOIETIN Specific gene validation by real-time polymerase chain reaction To validate expression, unique genes within IMMUNE RESPONSE enriched categories were measured for quantitative RNA expression using real-time reverse transcriptase-polymerase chain reaction. To accomplish this, specimens were chosen based on available RNA and included 6 HIVÀ and 19 HIV þ specimens that partially overlapped with specimens used in the microarray assessment (four of six HIVÀ and nine of 19 HIV þ ).

Results were normalized to endogenous 18S RNA levels to control for RNA quantity and are shown in Figure 4 . Log-transformed differences between HIVÀ and HIV þ levels were evaluated by t-test and were highly significant: ADAR (P ¼ 0001), APOBEC3G (P ¼ 0078), MX1 (Po0.0001), IRF7(Po0.0001), MYD88 (P ¼ 0002), RANTES (Po0.0001) and STAT 1(P ¼ 0001).

Because we were unable to identify a second population of drug naive subjects in Botswana (based on existing ethical guidelines); we therefore chose to independently validate the host profiles observed in this study by asking whether the enriched biological categories representing the HIV signature among peripheral blood mononuclear cells (PBMCs) in this study were comparable to other PBMC-based profiles for infection. To this end, we identified two studies: one comparing HIV-1 infection in PBMCs with healthy controls in a US Army cohort 40 and a second study evaluating host response in PBMCs to an unrelated pathogen, that is, severe acute respiratory virus (SARS). 69 Although the HIV-1 retrovirus and the SARS coronavirus are both plus-strand RNA viruses, their mode of infection, viral life cycle and pathogenic sequelae are distinct (for a review see Montano and Williamson 70 and Ziebuhr ,71 ). All three data sets were analyzed using the exact same overlapping gene list (8793 gene probes), namely, our HIV infection data set (set 1), the US Army HIV infection data set (set 2) and the SARS infection data set (set 3). As shown in Figure 5 , IMMUNE RESPONSE categories for HIV infection (both in the Botswana and the US Army data sets) included the same top four categories (top 20 are shown) in contrast with the RNA associated categories that were predominant in the SARS infection data set. This supports the view that HIV infection in distinct data sets elicited many (but not all) of the same enriched category gene-sets in contrast to host response to infection with a distinct viral pathogen (i.e. SARS).

We identified patterns of expression for specific gene-sets that were related to infection status, transmission outcome and viral burden. In this study, HIV-1C infection appeared to be associated with the differential expression of multiple gene-sets representing a broad innate response, characterized by an activation of TLR, interferon and antiviral RNA response pathways. These response pathways are functionally related. 72, 73 Recent studies indicate that TLRs trigger interferon associated genes (e.g. IFN-a/b) to initiate adaptive immunity by providing a link between the innate and adaptive immune response to infection, 74, 75 with subsequent influence on the expression of co-stimulatory molecules (e.g. CD80/86, class II), 76 CTL/CD8 þ T cell effector activity [77] [78] [79] and antigen presentation. 80 Our data indicate that HIV-1C infection was associated with a differential expression of innate response genes in the TLR pathway, including IL-1A, MYD88, RIP2/RIPK2, IRAK3/IRAKM, TRIF/TICAM1, NFKB2 and IP-10/ CXCL10. Upregulation of the adaptor protein MYD88 and TRIF suggested that both MYD88 dependent and MYD88 independent (TRIF mediated) pathways are engaged in HIV-1C infection. Members of the TLR family of receptors mediate innate immune response to a broad range of microbial ligands via activation of members of the REL, IRF and STAT transcription factor families and their respective target genes. MYD88 dependent effectors include proinflammatory and chemotactic cytokines, whereas the MYD88 independent effectors are associated with type I/II interferons and stimulation of the JAK-STAT pathway. Studies in vitro have shown that HIV-1 RNA can activate TLR signaling 81 and that microbial TLR engagement can activate HIV-1 transcription 82,83 and proinflammatory chemokine release. 84 HIV-1C infection was also associated with the differential expression of interferon-stimulated genes (e.g. STAT1, TRIM22, MX1, ISGF3G, IRF2, IRF7, IFI27, CXCR3 and PRKR. Interferons are a family of proteins produced in response to viral infection (notably RNA viruses) and/ or microbial activation through TLRs and various other cytokine signaling pathways, including RNA degradation and editing responses. Interferon ligand-receptor interactions stimulate JAK-STAT signaling that induce various IRFs, which in turn upregulate host chemotatic effector genes (e.g. IP-10, MIG, I-TAC, MCP-1) and multiple antiviral RNA response effectors. 85 STATs are activated by multiple cytokines and interferons (e.g. IFNa/b, IFN-g). 86 Multiple STATs are activated by HIV-1 infection in vitro 87 and chronic HIV-1 infection in vivo. 88 Mechanistic studies in vitro also implicate a role for IRFs in HIV-1 expression through the HIV-1 long terminal repeat (LTR) [89] [90] [91] and the HIV-1 transactivator protein pTAT, 92 suggesting that host induction of interferon and antiviral RNA response may be beneficial to the virus by influencing replication. Interferon stimulated genes in the peripheral blood have also been detected in acute infection using an SIV/HIV-1 chimeric virus, SHIV89.6 93 and have been detected in lymph node biopsies from HIV-1 infected subjects. 41 Our data also indicated a differential expression of interferon associated antiviral RNA response genes (e.g. MX1, PKR, OAS, ADAR, AP0BEC3G). Many of these genes are associated with type I/II interferon response (for a review see Samuel 94 ) and may influence transmission. 95 Activation of these genes are often associated with interferon-induced response to RNA viral infection. However, the upregulation of the RNA editing gene APOBEC3G that we observe (see Figures 2 and 4) has not previously reported. However, an examination of the distribution of predicted regulatory elements within the promoter region for APOBEC3G suggests the presence of multiple IRF binding sites (data not shown). Initial in vitro studies of HIV-1 infected cell lines did not show activation of APOBEC3G expression, potentially suggesting differences between cell types or in vitro/in vivo differences. Interestingly, APOBEC3G has been associated with G-to-A hypermutations of coding sequence for viral and cellular genes and restricts viral replication, in the absence of the HIV-1 vif gene. 96 Hypermutation of transmitted HIV sequences implicating RNA editing activity has been noted among newborns in Tanzania. 97 The activation of APOBEC3G and other antiviral RNA response genes may in part represent an ancient innate response to invading viral RNA that is engaged in addition to adaptive immunity. 98 Also implicated in these profiles were a subset of genes coding for SR proteins associated with RNA splicing, notably in the TR subset compared with HIVÀ and in the two subsets of NTR subjects. Some of these proteins have been previously associated with HIV-1 RNA splice-site selection in vitro, 99 and evidence for activation of different SR proteins has also been described during HIV infection in vitro 100, 101 and in vivo. 102 The relative abundance of SR proteins may influence the pattern of viral and host gene expression to influence local viral production and potentially transmission likelihood. Interferon induced RNA editing and degradation has been linked to SR protein activity. 103 It is unclear to what extent the differential expression of genes associated with antiviral response represent a host limitation on viral replication or promote viral replication. The mechanism(s) that promote and regulate HIV-1 RNA processing and alternative splicing configurations need to be directly studied in relation to viral replication to better understand their potential role in transmission outcome. We speculate that during infection in vivo that there may be increased antiviral and RNA editing gene expression that, in turn, influence the ratio of various splicing factors (by influencing SR protein activity), thereby shifting the host and viral transcriptome in favor of viral replication. The outcome of this process may contribute to increased transmission likelihood. In this study, the role of viral load in explaining the differential expression of gene-sets was examined. Interestingly, most genes identified in association with HIV-1 infection and/or transmission were not associated with viral load, although a subset of genes associated with immune response/interferon were present in correlation with viral burden and NTRs also exhibited two subsets of RNA processing genes that differed in relation to viral load. Collectively, this may suggest that infection induced a broad innate immune response that was sensitive to infection but was predominantly insensitive to viral burden in the peripheral blood. Alternatively, innate response genes and viral burden are more closely related in specific cell subsets or in different biological compartments and are not apparent in PBMCs, which represent mixed cell subsets. Overall, the findings in this study document that HIV-1C infection and transmission status were associated with the expression of different functional groups of genes that form a bridge between the innate and adaptive immune response. Our findings point to specific gene-sets associated with innate immune response, RNA processing and splicing, antiviral RNA response. Furthermore, the presence of common features of host response induced during HIV-1 infection in different settings, despite ethnographic, gender and viral subtype differences, and in contradistinction with other infections (see Figure 5 ), seems promising. These data raise enthusiasm for the potential of utilizing gene profiling to detect and characterize pathogen-specific host responses. Direct evaluation of specific gene role(s) in local infection, and expression monitoring of these genes in at risk subjects, may help augment efforts to both understand pathogen-specific host response and intervene in viral-host mechanisms engaged during both HIV infection and perinatal transmission.

The study population consisted of a cross-sectional group of 20 HIVÀ Botswana mothers with a mean age ¼ 27 (ranging from 18 to 38) and 25 HIV-1 þ Botswana mothers with a mean age ¼ 25 (ranging from 17 to 44) living within four different study sites. We have previously described the viral burden of the HIV-1 þ subjects in relation to transmission outcome. 18 The HIV þ mothers who participated were not identified for the study until after their infants were 2-5 months old. The HIV-1 seronegative subject specimens were collected through co-enrollment in an ongoing substudy in 2003-2004 using the exact same specimen collection method (see below and Montano et al.) 18 TRs, NTRs and seronegative controls did not differ significantly with regard to age, clinical status or condition, parity, Cesarean section experience, breast-feeding practice or prevalence of sexually transmitted diseases (Montano et al. 18 and data not shown). All subjects were asymptomatic, but CD4 counts were unavailable at the time of specimen collection. Informed consent was obtained from all participating subjects. The Botswana Health Ministry, as well as the institutional review boards at the Harvard School of Public Health and the Boston University School of Medicine, approved this study.

Specimen collection and RNA processing All PBMCs were collected from mothers and infants, as described. 18 Approximately, 1 ml of blood was placed into a cryovial tube containing 4 ml of RNA/DNA stabilizing reagent (Roche), inverted and stored at À801C. Peripheral blood mononuclear cells were processed to obtain total nucleic acid using a modified protocol of the mRNA isolation kit for Blood (Roche), then processed for RNA isolation using the RNeasy extraction kit (Qiagen, Valencia, CA, USA). Trace DNA was removed from the RNA samples using the DNA Free kit (Ambion, Austin, TX, USA) and the RNA concentrations were determined using the NanoDrop-1000 (NanoDrop Technologies, Rockland, DE, USA). Five micrograms of each total RNA specimen was provided to the Boston University Microarray Facility for labeling, amplification and hybridization to a U133A 2.0 chip from Affymetrix (Santa Clara, CA, USA). Hybridization signals were read using an Affymetrix Genechip Scanner 3000 and processed with the statistical software GCOS v1.2.1, and raw intensity values were scaled to a target ¼ 500. The 22 777 gene probes were filtered based on at least 25% present calls resulting in data for 11705 gene probes. The specimens utilized for microarray analysis were a representative random sampling based solely on specimen availability and RNA quantity.

Human ADAR, APOBEC3G, MX1, MYD88, STAT1, RANTES, IRF7 and 18 s rRNA (endogenous control) were measured using Assays-on-Demand (AoD) from Applied Biosystems Inc. (Foster City, CA, USA). Fifty nanograms of each sample RNA was used per reaction in duplicate and were normalized to each sample's corresponding 18s fluorescence value. Normalized values were log-transformed (due to skewing in some values) and then P-values for significance were evaluated using a two-sample Student's t-test.

Two independent PBMC profile Gene Ontology (GEO) data sets were identified (GSE2171 and GSE1739) representing an HIV-1 infection series and an SARS infection series, respectively. Since both of these data sets used a smaller focused chip (8793 gene probes) than the chip used in this study (22 777 gene probes), we identified the overlapping gene probes (8793) present in our (this study) HIV þ /HIVÀ Botswana data set, the US Army HIV þ /HIVÀ data set (GSE2171) and the SARS þ / SARSÀ data set (GSE1739). BADGE and EASE analysis was conducted on each infection/control series for each of the three data sets, as described in the statistical analysis section below. The top 20 categories representing enriched gene sets are shown. The entire category list (B200 categories) with significance greater than P ¼ 0.05 (adjusted Fisher's exact) is available in the web supplement. The Affymetrix data sets can be accessed under the GEO accession number GSE4124.

The arrays data sets were analyzed for differential expression based on HIV status (seronegative, seropositive) and transmitter status (TR; NTR) using BADGE (Bayesian Analysis of Differential Gene Expression) version 1.0, a computer program implementing a Bayesian approach to identify differentially expressed genes across experimental conditions. [64] [65] [66] The statistical procedure is described in detail in Supplemental Methods. The algorithm is also described in Klings et al 64 and Sebastiani et al 66, 104 and is available from (http://people.bu.edu/sebas/software.htm). The differential expression of each gene in two conditions is estimated by the fold change and measured by the probability that the fold change exceeds a fixed threshold, conditional on the data. To compute this probability, BADGE uses model averaging to gain robustness over model misspecifications. The current implementation of BADGE uses two models for the gene expression data: log-normal and gamma distributions. By combining both models, BADGE gains robustness and reproducibility over simpler analyses. In the analysis we used an expected false positive rate of 1%, and choose those genes that changed expression by at least 1.5-fold.

Once differentially expressed genes were identified by BADGE, the biologically enriched categories were identified, as recently described, 105 by implementing a stand-alone version of the EASE statistical software. 67 This program computes a modified Fisher's exact probability score for observing the frequency of biological category associated with a variable (infection, transmission), compared with the likelihood of seeing that category by chance given the total number of gene probes in the data set. An adjusted score is then reported representing the upper bound of the distribution of Jackknife Fisher exact probabilities for observing an enriched biological category. For more detail, see Hosack et al. 67 

Risk for perinatal HIV-1 transmission according to maternal immunologic, virologic, placental factors

Maternal viral load, vertical transmission of HIV-1: an important factor but not the only one. The European Collaborative Study

Perinatal transmission of the human immunodeficiency virus type 1 to infants of seropositive women in Zaire

Vertical transmission of human immunodeficiency virus (HIV) infection. Reactivity of maternal sera with glycoprotein 120 and 41 peptides from HIV type 1

Lack of association between maternal antibodies to V3 loop peptides and maternal-infant HIV-1 transmission

Comparison of mother-to-child transmission rates in Ugandan women with subtype A versus D HIV-1 who received single-dose nevirapine prophylaxis: HIV Network For Prevention Trials 012

The vertical transmission of human immunodeficiency virus type 1: molecular and biological properties of the virus

Preferential in-utero transmission of HIV-1 subtype C as compared to HIV-1 subtype A or D

HIV type 1 chemokine receptor usage in mother-to-child transmission

Concordance between the CC chemokine receptor 5 genetic determinants that alter risks of transmission and disease progression in children exposed perinatally to human immunodeficiency virus

Co-receptor usage of HIV-1 primary isolates, viral burden, and CCR5 genotype in motherto-child HIV-1 transmission

Genotypic characterization of human immunodeficiency virus type 1 isolated from vertically infected children with antiretroviral therapy experience

HIV-1 genotypic zidovudine drug resistance and the risk of maternal-infant transmission in the women and infants transmission study. The Women and Infants Transmission Study Group

Evolution and biological characterization of human immunodeficiency virus type 1 subtype E gp120 V3 sequences following horizontal and vertical virus transmission in a single family

Reduction of maternal-infant transmission of human immunodeficiency virus type 1 with zidovudine treatment. Pediatric AIDS Clinical Trials Group Protocol 076 Study Group

Estimated global distribution and regional spread of HIV-1 genetic subtypes in the year 2000

Botswana 2003 second generation HIV/AIDS surveillance: a technical report

Comparative prediction of perinatal human immunodeficiency virus type 1 transmission, using multiple virus load markers

Maternal virus load during pregnancy and mother-to-child transmission of human immunodeficiency virus type 1: the French perinatal cohort studies. SEROGEST Cohort Group

Maternal HIV-1 viral load and vertical transmission of infection: the Ariel Project for the prevention of HIV transmission from mother to infant

Maternal viral load, zidovudine treatment, and the risk of transmission of human immunodeficiency virus type 1 from mother to infant. Pediatric AIDS Clinical Trials Group Protocol 076 Study Group

Short-course zidovudine for prevention of perinatal infection

Effect of maternal CD4+ cell count, acquired immunodeficiency syndrome, and viral load on disease progression in infants with perinatally acquired human immunodeficiency virus type 1 infection. New York City Perinatal HIV Transmission Collaborative Study Group

Influence of other maternal variables on the relationship between maternal virus load and mother-toinfant transmission of human immunodeficiency virus type 1

Maternal plasma human immunodeficiency virus type 1 RNA level: a determinant and projected threshold for mother-to-child transmission

Identification of levels of maternal HIV-1 RNA associated with risk of perinatal transmission. Effect of maternal zidovudine treatment on viral load

Maternal levels of plasma human immunodeficiency virus type 1 RNA and the risk of perinatal transmission. Women and Infants Transmission Study Group

Effect of pregnancy and zidovudine therapy on viral load in HIV-1-infected women

Maternal viral load and timing of mother-to-child HIV transmission, Bangkok, Thailand. Bangkok Collaborative Perinatal HIV Transmission Study Group

Maternal virus load and perinatal human immunodeficiency virus type 1 subtype E transmission, Thailand. Bangkok Collaborative Perinatal HIV Transmission Study Group

Maternal plasma viral RNA levels determine marked differences in mother-to-child transmission rates of HIV-1 and HIV-2 in The Gambia. MRC/Gambia Government/University College London Medical School working group on mother-child transmission of HIV

Serum level of maternal human immunodeficiency virus (HIV) RNA, infant mortality, and vertical transmission of HIV in Zimbabwe

Maternal cell-free viremia in the natural history of perinatal HIV-1 transmission: a meta-analysis

Short-course zidovudine for perinatal HIV-1 transmission in Bangkok, Thailand: a randomised controlled trial. Bangkok Collaborative Perinatal HIV Transmission Study Group

Genetic determinants of pediatric HIV-1 infection: vertical transmission and disease progression among children

Maternal SDF1 3 0 A polymorphism is associated with increased perinatal human immunodeficiency virus type 1 transmission

CCR5 promoter polymorphisms in a Kenyan perinatal human immunodeficiency virus type 1 cohort: association with increased 2-year maternal mortality

A polymorphism in the regulatory region of the CC-chemokine receptor 5 gene influences perinatal transmission of human immunodeficiency virus type 1 to African-American infants

Interferon-gamma and interleukin-10 production among HIV-1-infected and uninfected infants of HIV-1-infected mothers

Functional genomic relationships in HIV-1 disease revealed by gene-expression profiling of primary human peripheral blood mononuclear cells

Functional genomic analysis of the response of HIV-1-infected lymphatic tissue to antiretroviral therapy

Host determinants in HIV infection and disease. Part 1: Cellular and humoral immune responses

Host determinants in HIV infection and disease. Part 2: Genetic factors and implications for antiretroviral therapeutics

HIV and AIDS: 20 years of science

The effect of genetic variation in chemokines and their receptors on HIV transmission and progression to AIDS

The influence of HLA genotype on AIDS

Genetic restriction of HIV-1 infection and progression to AIDS by a deletion allele of the CKR5 structural gene. Hemophilia Growth and Development Study, Multicenter AIDS Cohort Study, Multicenter Hemophilia Cohort Study, San Francisco City Cohort, ALIVE Study

Genetic acceleration of AIDS progression by a promoter variant of CCR5

Reduced risk of AIDS lymphoma in individuals heterozygous for the CCR5-delta32 mutation

Contrasting genetic influence of CCR2 and CCR5 variants on HIV-1 infection and disease progression. Hemophilia Growth and Development Study (HGDS)

CCR2-64I allele and genotype association with delayed AIDS progression in African women. University of Nairobi Collaboration for HIV Research

Modulating influence on HIV/AIDS by interacting RANTES gene variants

Genetic restriction of AIDS pathogenesis by an SDF-1 chemokine gene variant. ALIVE Study, Hemophilia Growth and Development Study (HGDS)

Genetic influence of CXCR6 chemokine receptor alleles on PCP-mediated AIDS progression among African Americans

MCP-1-MCP-3-Eotaxin gene cluster influences HIV-1 transmission

Genetic restriction of HIV-1 pathogenesis to AIDS by promoter alleles of IL10

A tumor necrosis factor-alpha-inducible promoter variant of interferon-gamma accelerates CD4+ T cell depletion in human immunodeficiency virus-1-infected individuals

Effect of a single amino acid change in MHC class I molecules on the rate of progression to AIDS

HLA and HIV-1: heterozygote advantage and B*35-Cw*04 disadvantage [see comment

HLA class I homozygosity accelerates disease progression in human immunodeficiency virus type 1 infection

Epistatic interaction between KIR3DS1 and HLA-B delays the progression to AIDS

APOBEC3G genetic variants and their influence on the progression to AIDS

Human genes that limit AIDS

Differential gene expression in pulmonary artery endothelial cells exposed to sickle cell plasma

Design and Analysis of Screening Experiments with Microarrays

The Analysis of Gene Expression Data: Methods and Software

Identifying biological themes within lists of genes with EASE

HIV-1 burden influences host response to coinfection with Neisseria gonorrhoeae in vitro

Expression profile of immune response genes in patients with severe acute respiratory syndrome

The molecular virology of HIV-1

Molecular biology of severe acute respiratory syndrome coronavirus

TLR signalling and activation of IRFs: revisiting old friends from the NF-kappaB pathway

IRF3 mediates a TLR3/TLR4-specific antiviral gene program

Links between innate and adaptive immunity via type I interferon

Regulation of the type I IFN induction: a current view

Genetic analysis of innate immunity: identification and function of the TIR adapter proteins

Induction of bystander T cell proliferation by viruses and type I interferon in vivo

Cross-priming of CD8+ T cells stimulated by virusinduced type I interferon

Combined TLR and CD40 triggering induces potent CD8+ T cell expansion with variable dependence on type I IFN

Regulation of phagosome maturation by signals from toll-like receptors [see comment

Species-specific recognition of singlestranded RNA via toll-like receptor 7 and 8

The induction of Tolllike receptor tolerance enhances rather than suppresses HIV-1 gene expression in transgenic mice

Rac1 and Toll-IL-1 receptor domain-containing adapter protein mediate Toll-like receptor 4 induction of HIV-long terminal repeat

Stimulation of toll-like receptor 2 in mononuclear cells from HIV-infected patients induces chemokine responses: possible pathogenic consequences

IRF family of transcription factors as regulators of host defense

Jaks and STATs: biological implications

Human immunodeficiency virus type 1 (HIV-1) induces activation of multiple STATs in CD4+ cells of lymphocyte or monocyte/macrophage lineages

Constitutive activation of STATs upon in vivo human immunodeficiency virus infection

IRF regulation of HIV-1 long terminal repeat activity

On the role of interferon regulatory factors in HIV-1 replication

Modulation of human immunodeficiency virus 1 replication by interferon regulatory factors

HIV-1 Tat reprograms immature dendritic cells to express chemoattractants for activated T cells and macrophages

Gene expression profiling of host response in models of acute HIV infection

Antiviral actions of interferons

Production of interferons and beta-chemokines by placental trophoblasts of HIV-1-infected women

Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein [see comment]

Hypermutation of HIV type 1 genomes isolated from infants soon after vertical infection

RNA interference: the molecular immune system

Differential effects of the SR proteins 9G8, SC35, ASF/SF2, and SRp40 on the utilization of the A1 to A5 splicing sites of HIV-1 RNA

The expression of the essential nuclear splicing factor SC35 is altered by human immunodeficiency virus infection

Identification and characterization of differentially expressed mRNAs in HIV type 1-infected human T cells

Gene expression and viral prodution in latently infected, resting CD4+ T cells in viremic versus aviremic HIV-infected individuals

Involvement of SR proteins in mRNA surveillance

Genetic dissection and prognostic modeling of overt stroke in sickle cell anemia

Evidence for cross-regulated cytokine response in human PBMCs exposed to whole gonococcal bacteria, in vitro

Supplementary Information accompanies the paper on Genes and Immunity's website

Genomic profile of HIV-1 infection and transmission M Montano et al