A Novel Sample Selection Approach to Aid the Identification of Factors That Correlate With the Control of HIV-1 Infection

Individuals infected with HIV display varying rates of viral control and disease progression, with a small percentage of individuals being able to spontaneously control infection in the absence of treatment. In attempting to define the correlates associated with natural protection against HIV, extreme heterogeneity in the datasets generated from systems methodologies can be further complicated by the inherent variability encountered at the population, individual, cellular and molecular levels. Furthermore, such studies have been limited by the paucity of well-characterised samples and linked epidemiological data, including duration of infection and clinical outcomes. To address this, we selected 10 volunteers who rapidly and persistently controlled HIV, and 10 volunteers each, from two control groups who failed to control (based on set point viral loads) from an acute and early HIV prospective cohort from East and Southern Africa. A propensity score matching approach was applied to control for the influence of five factors (age, risk group, virus subtype, gender, and country) known to influence disease progression on causal observations. Fifty-two plasma proteins were assessed at two timepoints in the 1st year of infection. We independently confirmed factors known to influence disease progression such as the B^*57 HLA Class I allele, and infecting virus Subtype. We demonstrated associations between circulating levels of MIP-1α and IL-17C, and the ability to control infection. IL-17C has not been described previously within the context of HIV control, making it an interesting target for future studies to understand HIV infection and transmission. An in-depth systems analysis is now underway to fully characterise host, viral and immunological factors contributing to control.

Introduction

Individuals infected with HIV display varying rates of viral control and disease progression, with a small percentage being able to spontaneously control in vivo viral replication without the need for anti-retroviral treatment (ART) (1). Such exquisite control is likely to happen in the very early battle between host and virus in acute and early HIV infection (2). Our understanding of host-pathogen interactions and the mechanisms underpinning the immune response to HIV infection have been informed by studies of individuals who demonstrate an enhanced ability to control in vivo viral replication, and on non-pathogenic SIV infection in non-human primates (NHP) (3). However, many of the studies of HIV control are cross-sectional after set point viral load and control has been achieved. Many of these studies have been focussed in Clade A and C infection. A full understanding of the mechanisms governing such spontaneous control of infection has been hampered by the paucity of informative and linked samples coupled to technology with sufficient resolution to define this phenomenon.

Systems-based approaches have helped define novel factors driving disease progression and protection during infections such as tuberculosis (4, 5), yellow fever (6), malaria (7), and influenza (8). But their application to aid the definition of the drivers of spontaneous control in HIV has been limited. The gene signature analysis of early gut mucosal T cell responses to HIV-1 suggest that the absence of an inflammatory gene signature may define Long-term non-progressors (LTNPs) (9). But recent scRNA-Seq profiling during acute HIV infection in a limited number of treatment-naïve subjects from the Females Rising through Education, Support and Health (FRESH) (10) cohort described an interferon response gene signature before peak viraemia as well as the presence of gene modules associated with antiviral control (APOBEC3A, IFITM1, and IFITM3) in individuals able to naturally control infection (11). The post hoc integrated systems analysis to the RV144 trial samples also uncovered roles for Type I and II interfons, as well as IRF7 and mTORC1 in susceptibility to infection post-vaccination (12). The mammalian target of rapamycin metabolic pathway has also been shown to be key to enhanced CD8 activity in elite controllers (13). These studies highlight the potential to utilise systems methods to define the correlates associated with the control of HIV-1 infection.

Heterogeneity in the data generated using high throughput systems methodologies can be further complicated by the inherent variability encountered at the population, individual, cellular and molecular levels (14). Studies by Chowdhury et al. (13) and others (15, 16) have highlighted the diversity of transcriptional profiles that exist within a single subset of T lymphocytes that accounts in part for control of HIV infection. The control of HIV replication in vivo is multifactorial. Indeed viral control has been shown to be associated with age at infection, time post-infection, gender, HLA type, virus subtype and route of infection (17–22). Obtaining sufficient numbers of samples to allow for the control of all these confounders and the discovery of new correlates of disease trajectory poses a real challenge (23, 24).

We applied a unique approach to retrospectively classify HIV-infected individuals in order to aid the delineation of a profile associated with early and persistent in vivo control of HIV-1 replication in the absence of antiretroviral treatment. Using this approach, we defined three groups of HIV infected volunteers from Protocol C; a multisite early infection prospective cohort consisting of 613 participants recruited from nine clinical research centres in five African countries (25, 26) (Figure 1, also on https://dataspace.iavi.org/). These groups comprised volunteers with low (n = 10), medium (n = 10) and high (n = 10) set point viral load who were identified within days of their estimated date of HIV infection and followed over time for up to 7 years. Importantly, the low viral load volunteers showed rapid and persistent control of viral replication in the absence of treatment, and had sufficient samples available during the resolution of peak infection to enable the investigation of signatures associated with rapid and persistent HIV control. We present the profile for fifty-two soluble proteins in the acute phase of HIV infection across the three groups, demonstrating the potential to identify unique signatures associated with ART-naïve viral control using this selection approach.

FIGURE 1

Figure 1. Derivation by country of the 30 volunteers used in this study from the African incidence cohort [Protocol C (22)].

Methods

Ethics

This study was reviewed and approved by the following ethical review boards: the Kenya Medical Research Institute Ethical Review Committee, the Kenyatta National Hospital Ethical Review Committee of the University of Nairobi, the Rwanda National Ethics Committee, the Uganda Virus Research Institute Science and Ethics Committee (Currently the UVRI Research Ethics Committee) and the Uganda National Council of Science and Technology, the University of Cape Town Health Science Research and Ethics Committee, the Bio-Medical Research Ethics Committee at the University of KwaZulu Natal, the University of Zambia Research Ethics Committee, and the Emory University Institutional Review Board. Informed consent was obtained from all volunteers prior to the collection of study related resource. All methods were carried out in accordance with relevant guidelines and regulations.

Study Population and Selection Approach

Volunteers included in this study were selected from a historic acute and early HIV infection prospective cohort drawn from nine clinical research centres in South Africa, Zambia, Uganda, Kenya and Rwanda enrolled from 2006 to 2011 (Figure 1). Details of study characteristics, distributions, recruitment procedures, initial immunological methods and epidemiological profiling data of the Protocol C cohort are described elsewhere (18, 22, 25, 26).

Individuals from the study were ranked according to the magnitude of their mean viral load (Geometric mean) measurements taken between 9–36 months post-EDI (estimated day of infection), and divided into quartiles. Mean viral load was calculated for 362 of the 613 volunteers from the Protocol C cohort who did not receive antiretroviral treatment. A matching algorithm (27) based on the nearest neighbour was then applied to define the groups of volunteers for this study. Briefly, we selected 10 Low viral load volunteers (LVLVs) from the first quartile of the ranked dataset who had a visible period of dynamic control of infection demonstrated by the presence of a downslope in their viral load measurements, and who were able to control viral load to ≤2,000 copies/mL in the first 3 years of infection. The LVLVs were then matched with Intermediate viral load volunteers (IVLVs, n = 10) from the second and third quartiles, and High viral load volunteers (HVLVs, n = 10) from the fourth quartile of the ranked dataset. Volunteers were matched on age, clade, country, gender and risk group. Soluble proteins in plasma were assessed at two timepoints in the period following peak viraemia and within the 1st year of infection for each of the selected volunteers (Supplementary Table 1). All individuals assessed for this study were treatment naïve.

HLA Frequency Calculation

To determine the HLA I frequencies within the 362 ART-naïve Protocol C volunteers, two-digit allelic frequencies were calculated using the Los Alamos National Laboratory HLA frequency and Graphing tool (https://www.hiv.lanl.gov/content/immunology/hla). For each MHC Class I alleles with an allele frequency >5%, we compared the set point viral load of all positive volunteers with those of all negative volunteers. Statistical tests used are described in subsequent sections.

Plasma Analyte Quantification

Fifty-two soluble analytes were quantified in plasma using a combination of six Meso Scale Discovery (MSD) human V-PLEX panels including the Angiogenesis Panel 1 (VEGF-A, VEGF-C, VEGF-D, Tie-2, Flt-1, PIGF, bFGF), TH17 Panel 1 (IL-17A, IL-21, IL-31, IL-27, IL-23, IL-22, MIP-3α), Chemokine Panel 1 (Eotaxin, MIP-1β, Eotaxin-3, TARC, IP-10, MIP-1α, IL-8, MCP-1, MDC, MCP-4), Cytokine Panel 1 (GM-CSF, IL-1α, IL-5, IL-7, IL-12/IL-23p40, IL-15, IL-16, IL-17A, TNF-β, VEGF-A), Cytokine Panel 2 (IL-1RA, IL-3, IL-9, IL-17A/F, IL-17B, IL-17C, IL-17D, TSLP), Proinflammatory Panel 1 (IFN-γ, IL-1β, IL-2, IL-4, IL-6, IL-8, IL-10, IL-12p70, IL-13, TNF-α), and Vascular Injury Panel 2 (SAA, CRP, VCAM-1, ICAM-1). Cryopreserved plasma from the incidence study were thawed at room temperature and applied to the panels according to the manufacturer's protocol. Plates were read on the MSD plate reader model MESO QuickPlex SQ 120. All plasma samples for the study were thawed and run at the same time and grouped on plates in the order in which they were selected for the study to avoid intra-assay variability. Data was collected for two replicates per sample using the MSD software (Discovery Workbench Version 4.0). A five-parameter logistic regression formula was used to derive sample concentrations from the standard curves. Analytes below the lower limit of detection were assigned a concentration of half the lower limit of quantification (LLOQ).

Statistical Analysis

We analysed the MSD data using a non-parametric approach, because of the small sample size and the non-Gaussian distribution, as determined using the Shapiro-Wilk test.

Non-parametric analysis (Mann-Whitney test, comparing ranks) of the differences in VL measurements for individuals expressing HLA alleles was performed in Graphpad Prism 8 Software. P-values <0.05 were considered significant. Propensity score matching to define study populations was executed using the MatchIT package (27) in R. Non-parametric test for similarities in the age distribution between the study groups was performed in SPSS 24.

We computed descriptive summary statistics, including the median and inter–quartile range (IQR) and Spearman correlations and excluded MIP-3α from further analyses due to missing data (40%). The null hypothesis of the difference between the two time points was assessed using Wilcoxon Signed Ranks tests. We computed the (rank based) correlation matrix for each group (LVLVs, IVLVs, HVLVs) by averaging the concentrations over time and presented the correlation matrices as heat maps.

To investigate the association between the mean viral load of volunteers (which was used to define the study groups) and the concentration of proteins in peripheral blood we fitted a linear robust regression model where a function of the ranks of the residuals was used instead of the Euclidian distance in the least square estimation (28). SAA and CRP were excluded from the model due to singularity issues. Based on these univariate results, we selected the significant proteins and estimated a multivariate robust regression model after adjusting for multiple comparison (Bonferroni). We then fitted a multivariate regression model with the significant (p-value adjusted 0.05) proteins and removed those that eventually where highly correlated to avoid multicollinearity issues in the regression fitting.

R analysis was conducted using version 3.1.2, 2014, (available at https://www.R-project.org).

Results

The Outcome of Infection Is Linked to Gender, Viral Subtype, and the Expression of Immune Receptors on Lymphocytes

Set-point viral load represents a dynamic state of equilibrium between infecting virus and the immune response in the absence of complete elimination of the virus (29). It remains an important measure of disease progression. Set-point viral load was calculated for 362 of the 613 volunteers who did not receive antiretroviral treatment over at least 36 months of follow-up. Volunteers were then ranked and divided into equal quartiles to explore any associations with disease progression (Figure 2A). Median set point VL for the 362 volunteers was 26,061 copies/ml (IQR: 6,981–65,813 copies/mL). Mean CD4 counts calculated for the same period for all 362 volunteers was inversely correlated with mean viral load (r = −0.2892, p < 0.0001) (Figure 2B).

FIGURE 2

Figure 2. (A) Ranks of set point VL calculated for 362 volunteers between 9 and 36 months post infection showing distribution of Low viral load volunteers (LVLVs), Intermediate viral load volunteers (IVLVs) and High viral load volunteers (HVLVs). Ranked dataset was also divided into equal quartiles based on set point viral load. (B) Set point VL correlates inversely with mean CD4 counts calculated over the same period (Pearson's r = −0.2892; p < 0.0001). (C) Subtype and gender distribution of ART naïve volunteers by quartiles of the ranked dataset. Infecting subtype and gender both have a relationship with set point VL.

We examined the distribution of gender and viral subtype within our ranked dataset. In agreement with previous studies (18, 22) our analysis showed that that women had lower set point viral loads than men. There was also a higher representation of Subtype A in the lower quartiles, with the opposite being true for Subtype C infected subjects (Figure 2C).

To assess the impact of MHC on disease progression, we compared the influence of Class I alleles with an allele frequency > 5% on the set point viral load and found that individuals with B^*57 (p < 0.0001) and C^*04 (p = 0.0335) had lower and higher set point viral loads, respectively, compared with individuals lacking either HLA allele (Figures 3A,B).

FIGURE 3

Figure 3. Distribution of set point VL based on the presence of MHC Class I (A,B). Violin plots show individual data points as well as the 25th, 50th, and 75th percentiles. The null hypothesis was tested using a non-parametric unpaired test (Mann-Whitney p < 0.05 considered significant).

A Propensity-Based Approach to Sampling an HIV Incidence Cohort to Aid Systems Analysis

During untreated HIV infection, the rate of viral replication and set-point probably reflects the dynamic interaction between the virus and host responses (23). We defined Low Viral Load Volunteers (LVLVs) as those with a set point of <2,000 copies/mL in line with previous studies (1) and were able to identify 40 Protocol C volunteers within this category. Ten of the 40 LVLVs identified had sufficient samples and viral load measurements collected over the period following peak viraemia for analysis. There was a steady decrease in viral load in the LVLVs over the 12 months following peak viraemia in the absence of ART (Figure 4A, Supplementary Figure 1). We focussed initially on the period immediately following peak viraemia in an effort to describe correlates of the early control of viral replication.

FIGURE 4

Figure 4. Locally weighted scatterplot smoothing (LOWESS) curves showing the overall trend in VL (black circles) and CD4 counts (red squares) for Low viral load volunteers (LVLVs) (A), Intermediate viral load volunteers (IVLVs) (B) and High viral load volunteers (HVLVs) (C) over 36 months post infection. N = 10 for each group. The shaded region indicates the dynamic period of immune control observed in the LVLVs following peak viraemia.

The ability to control viral replication in vivo has been linked to factors such as age at infection, time post-infection, gender, HLA type, route of virus entry and HIV subtype (17–22). In an effort to control for the confounding effects of some of these factors, we utilised a propensity score matching approach (27), which allows the matching of persons in one group with persons in another group based on each case's propensity score. For all the volunteers in the first group (LVLVs), we selected an equal number of volunteers from quartiles 2 and 3 (Figure 2A) based on their propensity scores for age, time post infection, gender, route of virus entry and infecting virus subtype. We designated the group selected from quartiles 2 and 3 as Intermediate viral load volunteers (IVLVs). We applied the same matching approach to volunteers in quartile 4 to identify 10 High viral load volunteers (HVLVs), from the ranked dataset. It was impossible to match on HLA haplotypes due to the diversity of alleles represented in the cohort.

We successfully identified three distinct groups consisting of 10 individuals per group, from the ranked dataset that were matched on age, gender, risk group (route of infection), country and infecting subtype (Figures 4A–C, Table 1) and for whom samples were available at two timepoints within the initial phase of control of viral replication immediately after peak viraemia. Given the age of the cohort, sample availability within this period (obtained from dataspace.iavi.org) was a real challenge. Days post EDI was also considered during the selection of these early timepoints with matched timepoints no more than 6 months apart where possible (Supplementary Table 1).

TABLE 1

Table 1. Characteristics of groups selected from the HIV incidence cohort.

Concentrations of Soluble Markers in the Acute Phase Are Associated With Early and Sustained Control of in vivo Viral Load

We measured the levels of 52 soluble proteins in plasma at two timepoints following peak viraemia and report the median and IQR for the two timepoints for all volunteers (Table 2). For the most part, protein concentrations were not significantly different across the two early timepoints assessed with the exception of VCAM1 (p = 0.02) and IL-10 (p = 0.012) for the overall dataset, IL-6 (p = 0.037) and IL-17C (0.006) for LVLVs, SAA (p = 0.027) and CRP (p = 0.027) for IVLVs, and IFN-γ (p = 0.01) for HVLVs (group data also shown in Supplementary Figures 2, 3).

TABLE 2

Table 2. Descriptive statistics for all the proteins assessed.

To further investigate the association between the mean viral load of volunteers (which was used to define the study groups) and the concentration of proteins in peripheral blood we applied a univariate regression model where a function of the ranks of the residuals was used instead of the Euclidian distance in the least square estimation (28). The estimates for nine analytes that were significant are shown in Table 3 (GM-CSF, p < 0.001; MIP 1-α, p < 0.001; IL-8, p < 0.001; IFN-γ, p < 0.001; IL-2, p < 0.001; IL-13, p = 0.05; IL-17C, p < 0.001; IL-9 = 0.02; IL-31 = 0.02).

TABLE 3

Table 3. Model estimates from the univariate non-parametric model to investigate the association between the mean viral load of volunteers and the concentration of proteins in peripheral blood.

Based on the results of the univariate analyses, we generated a multivariate robust regression model using the proteins which were associated with mean viral load, after adjusting for multiple comparison (Bonferroni) and excluding highly correlated analytes to avoid issues of multicollinearity. The only protein that remained significantly associated with mean viral load was MIP1-α (p < 0.001) after p-value adjustment (Holmmel). IL-8 was excluded from the model because it was highly correlated with GM-CSF (Spearman correlation, 0.71).

Exploratory heatmaps based on the lower triangle Spearman correlation matrices and using the average plasma protein value between the two timepoints for each group suggests that differences exist in the relationships between different plasma proteins across the groups (Figure 5) with more frequent positively correlated proteins seen in IVLVs and HVLVs compared to LVLVs.

FIGURE 5

Figure 5. Correlograms of the correlations between 52 plasma protein concentrations for Low viral load volunteers (LVLVs), Intermediate viral load volunteers (IVLVs) and High viral load volunteers (HVLVs). Blue and red squares represent positive and negative correlations, respectively with darker colours indicating a greater magnitude of correlation.

Discussion

We present a unique approach to classifying individuals drawn from an acute and early HIV infection cohort that considers a range of factors known to have an impact on disease progression, to efficiently define the peripheral secretory profile associated with early and sustained control of in vivo viral replication. This selection approach enabled us to define the profile of 52 proteins deployed within the specific period of dynamic immunological control of viral replication for all volunteers in the absence of antiretroviral treatment. Expectedly our ranking approach show that measurements for CD4 cells, which are the first cells to become infected during transmission (30, 31) and continue to be a primary target for HIV-1 (32), were negatively correlated with set point viral load calculations.

HIV subtype has been shown to be associated with disease trajectory and outcome (21, 22, 33, 34) but limitations of study size and design have meant that such findings have been largely descriptive in nature. By matching the selected individuals from Protocol C based on other relevant confounding variables, we were able to address the independent contribution of viral subtype to the control of viral replication in vivo. We report similar observations to Price et al. (22) who performed a subtype-by-geographic-region covariate analysis on the whole Protocol C cohort and showed that Subtype A-infected volunteers were more likely to control viral load than Subtype-C infected volunteers. Also in this cohort, Amornkul et al. (18) show that subtype C is associated with faster progression to AIDS and CD4+ T cell decline compared to subtype A. Here we also show a higher representation of Subtype-A relative to Subtype-C infected subjects in the lower quartiles of the viral load-ranked ranked dataset of the same cohort.

Whilst the diversity of HLA types represented in the cohort did not permit complete matching of volunteers based on this factor, well-reported trends like the favourable influence of B^*57 on disease control were evident. The presence of the less studied C^*04 HLA Class I allele, which is reportedly associated with B^*35 on chromosome 6 (35, 36) appeared to have a less favourable influence on disease control in this study. This potentially deleterious effect of the C^*04 allele on disease progression has only been reported by a few studies that focused on either HIV-1 subtype B alone (35–37) or mostly C (38). Our data covering HIV-1 subtype A, C, D and recombinant viral subtypes suggests that this effect of C^*04 may apply regardless of the subtype of the virus. It has been suggested that this effect may be mitigated by the association of C^*04 with other deleterious HLA Class I alleles or with the killer-cell immunoglobulin-like receptor (KIR) KIR2DS4 (37). These observations provide some validation of the novel propensity matching method presented in this study.

Early HIV infection is characterised by a cytokine storm that is detectable at the levels of gene (11) and protein expression (39). We were able to focus on the dynamics of the plasma protein response following peak viraemia and during the 1st year of infection when the immune system is most actively involved in the control of viral replication. The fact that the circulating levels of most of the plasma proteins appeared to be largely stable in all volunteers with a few exceptions (IL-6, IL-17C, SAA, CRP, and IFN-γ) over the period assessed was somewhat surprising. For the most part the indicated analytes have been linked previously to HIV disease outcomes (40, 41), even if the strength of their associations remain poorly understood. A notable exception was the falling levels of IL-17C observed in LVLVs. In contrast to other members of the IL-17 family (IL-17A and IL-17F), IL-17C is predominantly produced by epithelial cells (42–44) and not leucocytes (45), with broad activity on epithelial cells, TH17 leucocytes (46) and monocytic lineage cells (43). Although not fully described within the context of HIV-1 infection, early release of IL-17C in other models of infection suggest a dual function in the regulation of both innate and adaptive immune responses (45). The decreasing levels of serum IL-17C seen in LVLVs may suggest a role in early recruitment and differentiation of innate and adaptive modulators in response to HIV infection-i.e., prior to peak viral load, and subsequent downregulation in those who eventually go on to control infection. The difference in its cellular source compared to other members of the IL-17 cytokine family makes it a potentially interesting target for further studies to examine a potential role in barrier immunity during and immediately after transmission at mucosal surfaces.

Our regression analysis suggests a relationship between the levels of nine plasma proteins including IL-17C in the period following peak viraemia and set point viral load, with MIP-1α being the most significantly associated with mean viral load in the multivariate analysis. MIP-1α is one of three well-characterised β-chemokines produced by immune cells including CD8 and CD4 T cells that have been implicated in the inhibition of HIV infection (47, 48). Whilst the evidence strongly suggests that T cell capacity to secrete β-chemokines like MIP-1α is strongly associated with in vivo and in vitro viral control (49–52), the relationship between the quantity of circulating soluble β chemokines and disease outcome appears to be less clear (53–55). Our results suggest that the concentration of MIP-1α in the period following peak viraemia is strongly associated with set point viral load.

Several groups have assessed the relationships between specific cytokines and the control of viral replication (56–58). Our exploratory analysis of the correlations between 52 proteins during the period following peak viraemia suggest that disease progression is underpinned by clear differences in the deployment of immune modulators. Whilst the functional implications of the negatively correlated proteins in LVLVs is not immediately clear, this and previous observations validate the unique selection approach presented here and the potential to apply it to support systems investigation of the correlates that define natural control of HIV-1 infection.

The impact of biological sex on the outcome of viral infections has been highlighted by other groups (59, 60). Our ranked data set showed that women had lower set point viral loads than men. In comparing the two groups (men and women), we identified differences in the levels of IL-7, a cytokine implicated in early T cell development, proliferation and differentiation (61) in the period studied. Levels of thymus and activation-regulated chemokine (TARC), a chemokine constitutively produced in the thymus and by keratinocytes and dendritic cells (62, 63) with powerful chemoattractant effect on T cells was also found to be higher in men in this study. Whilst recent studies examining these sex differences point to a possible role for increased levels of innate inflammatory cytokines (60) on disease progression during viral infections, our data may indicate differences in how T cells are activated/recruited and deployed in men and women. A recent study by El-Badry et al. highlights the impact of plasma levels of 17β-estradiol in women on T cell activation in the acute phase of HIV infection. It is worth noting the small sample size from which these observations were made; confirmation of our observations will be required.

Taken together with previous results, it is reasonable to state that whilst our results suggest an association between levels of soluble MIP-1α in the period of active immune suppression of viral replication and disease progression, they also support the notion that the mechanism of in vivo suppression of HIV is likely multifactorial (64). As such efforts aimed at reducing the potential for noise in datasets will go a long way to enable the definition of the correlates associated with immunological control of HIV.

We present our unique selection approach as a way to potentially counter some of the noise associated with extreme heterogeneity in datasets allowing for the application of high-resolution systems methodologies to define the correlates associated with natural control of HIV infection. Whilst the ranking and propensity-based selection methods presented may not directly predict correlates of natural protection against HIV-1, they enable the exclusion of any noise arising as a result of the factors that are controlled for in the study design. Given that HIV pathogenesis is multifactorial, the tendency for such noise to obscure valid observations is considered a real barrier to the application of high dimensional (or systems) analytical methods to aid the definition of the correlates of natural control (14). We demonstrate the utility of the ranking and matching approach by independently confirming known factors associated with disease progression, albeit in a small study population. We also identify novel soluble proteins such as IL-17C and illustrate differences in the pattern of deployment of peripheral cytokines that tally with disease progression. These findings demonstrate the utility of the unique ranking and selection approach for systems analysis is subsequent studies.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

This study was reviewed and approved by the following ethical review boards: the Kenya Medical Research Institute Ethical Review Committee, the Kenyatta National Hospital Ethical Review Committee of the University of Nairobi, the Rwanda National Ethics Committee, the Uganda Virus Research Institute Science and Ethics Committee (Currently the UVRI Research Ethics Committee) and the Uganda National Council of Science and Technology, the University of Cape Town Health Science Research and Ethics Committee, the Bio-Medical Research Ethics Committee at the University of KwaZulu Natal, the University of Zambia Research Ethics Committee, and the Emory University Institutional Review Board. Informed consent was obtained from all volunteers prior to the collection of study related resource. All methods were carried out in accordance with relevant guidelines and regulations. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JM was responsible for conceptualisation, methodology, formal analysis, data curation, sample application preparation, original draft preparation, review, and editing. EN was responsible for methodology development, sample application preparation, and original draft review. AF-S was responsible for statistical analysis. CS was responsible for running immunological assays. CK, JD, SB, GM, and PH were responsible for methodology development and original draft review. JH was responsible for assay review, processing and approval of sample application. DK, SJ, EM, and BA were responsible for methodology development, data curation, assay review and original draft review. EH, MP, ES, and JG were responsible for project conceptualisation, methodology development, data curation, granting sample access, review, and editing. The IAVI protocol C investigators were responsible for the initiation and successful completion of the Protocol C study. All authors contributed to the article and approved the submitted version.

List of the IAVI Protocol C Investigators

• Eduard J. Sanders, Centre for Geographic Medicine -Coast/KEMRI, Kenya; University of Oxford, UK.

• Omu Anzala, Kenya AIDS Vaccine Institute -Institute of Clinical Research, Kenya.

• Anatoli Kamali, Medical Research Council/Uganda Virus Research Institute, Uganda Research Unit on AIDS, Uganda.

• Etienne Karita, Project San Francisco, Rwanda.

• William Kilembe, Mubiana Inambao, Shabir Lakhi, Zambia Emory Research Project, Zambia.

• Susan Allen, Eric Hunter, Emory University, Georgia, USA.

• Vinodh Edward, The Aurum Institute, South Africa.

• Pat Fast, IAVI, New York, USA.

• Matt A. Price, IAVI, New York, USA; Department of Epidemiology and Biostatistics, University of California San Francisco, USA.

• Jill Gilmour, IAVI Human Immunology Laboratory, Imperial College, UK.

• Jianming Tang, Ryals Public Health Building, University of Alabama, USA.

• Fran Priddy, IAVI, New York, USA.

• Mary H. Latka, The Aurum Institute, South Africa.

• Linda-Gail Bekker, Desmond Tutu HIV Foundation, South Africa.

Funding

This work was made possible by IAVI, which was supported by funding from many donors, including USAID, the Bill and Melinda Gates Foundation, the Ministry of Foreign Affairs of Denmark, Irish Aid, the Ministry of Finance of Japan in partnership with The World Bank, the Ministry of Foreign Affairs of the Netherlands, the Norwegian Agency for Development Cooperation, the United Kingdom Department for International Development, and the US Agency for International Development (the full list of IAVI donors is available at: http://www.iavi.org).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.634832/full#supplementary-material

References

1. Gurdasani D, Iles L, Dillon DG, Young EH, Olson AD, Naranbhai V, et al. A systematic review of definitions of extreme phenotypes of HIV control and progression. AIDS. (2014) 28:149–62. doi: 10.1097/QAD.0000000000000049

CrossRef Full Text | Google Scholar

2. Yue L, Pfafferott KJ, Baalwa J, Conrod K, Dong CC, Chui C, et al. Transmitted virus fitness and host T cell responses collectively define divergent infection outcomes in two HIV-1 recipients. PLoS Pathog. (2015) 11:e1004565. doi: 10.1371/journal.ppat.1004565

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Sodora DL, Allan JS, Apetrei C, Brenchley JM, Douek DC, Else JG, et al. Toward an AIDS vaccine: lessons from natural simian immunodeficiency virus infections of African nonhuman primate hosts. Nat Med. (2009) 15:861–5. doi: 10.1038/nm.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Sweeney TE, Braviak L, Tato CM, Khatri P. Genome-wide expression for diagnosis of pulmonary tuberculosis: a multicohort analysis. Lancet Respir Med. (2016) 4:213–24. doi: 10.1016/S2213-2600(16)00048-5

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Zak DE, Penn-Nicholson A, Scriba TJ, Thompson E, Suliman S, Amon LM, et al. A blood RNA signature for tuberculosis disease risk: a prospective cohort study. Lancet. (2016) 387:2312–22. doi: 10.1016/S0140-6736(15)01316-1

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Querec TD, Akondy RS, Lee EK, Cao W, Nakaya HI, Teuwen D, et al. Systems biology approach predicts immunogenicity of the yellow fever vaccine in humans. Nat Immunol. (2009) 10:116–25. doi: 10.1038/ni.1688

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Kazmin D, Nakaya HI, Lee EK, Johnson MJ, van der Most R, van den Berg RA, et al. Systems analysis of protective immune responses to RTS,S malaria vaccination in humans. Proc Natl Acad Sci U S A. (2017) 114:2425–30. doi: 10.1073/pnas.1621489114

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Nakaya HI, Wrammert J, Lee EK, Racioppi L, Marie-Kunze S, Haining WN, et al. Systems biology of vaccination for seasonal influenza in humans. Nat Immunol. (2011) 12:786–95. doi: 10.1038/ni.2067

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sankaran S, Guadalupe M, Reay E, George MD, Flamm J, Prindiville T, et al. Gut mucosal T cell responses and gene expression correlate with protection against disease in long-term HIV-1-infected nonprogressors. Proc Natl Acad Sci U S A. (2005) 102:9860–65. doi: 10.1073/pnas.0503463102

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Dong KL, Moodley A, Kwon DS, Ghebremichael MS, Dong M, Ismail N, et al. Detection and treatment of Fiebig stage I HIV-1 infection in young at-risk women in South Africa: a prospective cohort study. Lancet HIV. (2018) 5:e35–44. doi: 10.1016/S2352-3018(17)30146-7

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Kazer SW, Aicher TP, Muema DM, Carroll SL, Ordovas-Montanes J, Miao VN, et al. Integrated single-cell analysis of multicellular immune dynamics during hyperacute HIV-1 infection. Nat Med. (2020) 26:511–8. doi: 10.1101/654814

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Fourati S, Ribeiro SP, Blasco Tavares Pereira Lopes F, Talla A, Lefebvre F, Cameron M, et al. Integrated systems approach defines the antiviral pathways conferring protection by the RV144 HIV vaccine. Nat Commun. (2019) 10:863. doi: 10.1038/s41467-019-08854-2

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chowdhury FZ, Ouyang Z, Buzon M, Walker BD, Lichterfeld M, Yu XG. Metabolic pathway activation distinguishes transcriptional signatures of CD8+ T cells from HIV-1 elite controllers. AIDS. (2018) 32:2669–77. doi: 10.1097/QAD.0000000000002007

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Sekaly R. Pulendran B. Systems biology in understanding HIV pathogenesis and guiding vaccine development. Curr Opin HIV AIDS. (2012) 7:1–3. doi: 10.1097/COH.0b013e32834e0667

CrossRef Full Text | Google Scholar

15. Pereyra F, Addo MM, Kaufmann DE, Liu Y, Miura T, Rathod A, et al. Genetic and immunologic heterogeneity among persons who control HIV infection in the absence of therapy. J Infect Dis. (2008) 197:563–71. doi: 10.1086/526786

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Saez-Cirion A, Sinet M, Shin SY, Urrutia A, Versmisse P, Lacabaratz C, et al. Heterogeneity in HIV suppression by CD8 T cells from HIV controllers: association with Gag-specific CD8 T cell responses. J Immunol. (2009) 182:7828–37. doi: 10.4049/jimmunol.0803928

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Rodriguez MA, Ding M, Ratner D, Chen Y, Tripathy SP, Kulkarni SS, et al. High replication fitness and transmission efficiency of HIV-1 subtype C from India: Implications for subtype C predominance. Virology. (2009) 385:416–24. doi: 10.1016/j.virol.2008.12.025

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Amornkul PN, Karita E, Kamali A, Rida WN, Sanders EJ, Lakhi S, et al. Disease progression by infecting HIV-1 subtype in a seroconverter cohort in sub-Saharan Africa. AIDS. (2013) 27:2775–86. doi: 10.1097/QAD.0000000000000012

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yue L, Prentice HA, Farmer P, Song W, He D, Lakhi S, et al. Cumulative impact of host and viral factors on HIV-1 viral-load control during early infection. J Virol. (2013) 87:708–15. doi: 10.1128/JVI.02118-12

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Carlson JM, Schaefer M, Monaco DC, Batorsky R, Claiborne DT, Prince J, et al. HIV transmission. Selection bias at the heterosexual HIV-1 transmission bottleneck. Science. (2014) 345:1254031. doi: 10.1126/science.1254031

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Venner CM, Nankya I, Kyeyune F, Demers K, Kwok C, Chen PL, et al. Infecting HIV-1 subtype predicts disease progression in women of Sub-Saharan Africa. EBioMedicine. (2016) 13:305–14. doi: 10.1016/j.ebiom.2016.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Price MA, Rida W, Kilembe W, Karita E, Inambao M, Ruzagira E, et al. HIV-1 viral control varies by viral subtype in a large cohort of African adults with incident HIV-1 infection. J Infect Dis. (2019) 220:432–41. doi: 10.1093/infdis/jiz127

PubMed Abstract | CrossRef Full Text | Google Scholar

23. McMichael AJ, Borrow P, Tomaras GD, Goonetilleke N, Haynes BF. The immune response during acute HIV-1 infection: clues for vaccine development. Nat Rev Immunol. (2010) 10:11–23. doi: 10.1038/nri2674

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Nakaya HI, Pulendran B. Systems vaccinology: its promise and challenge for HIV vaccine development. Curr Opin HIV AIDS. (2012) 7:24–31. doi: 10.1097/COH.0b013e32834dc37b

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kamali A, Price MA, Lakhi S, Karita E, Inambao M, Sanders EJ, et al. Creating an African HIV clinical research and prevention trials network: HIV prevalence, incidence and transmission. PLoS ONE. (2015) 10:e0116100. doi: 10.1371/journal.pone.0116100

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Price MA, Kilembe W, Ruzagira E, Karita E, Inambao M, Sanders EJ, et al. Cohort profile: IAVI's HIV epidemiology and early infection cohort studies in Africa to support vaccine discovery. Int J Epidemiol. (2020) 1–10. doi: 10.1093/ije/dyaa100

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Daniel Ho E KI, King G, Stuart EA. MatchIt: nonparametric preprocessing for parametric causal inference. J Stat Soft. (2011) 42:1–28. doi: 10.18637/jss.v042.i08

CrossRef Full Text | Google Scholar

28. Hettmansperger TP, McKean JW. Robust Nonparametric Statistical Methods. CRC Press (2010).

Google Scholar

29. Mei Y, Wang L, Holte SE. A comparison of methods for determining HIV viral set point. Stat Med. (2008) 27:121–39. doi: 10.1002/sim.3038

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Miller CJ, Li Q, Abel K, Kim EY, Ma ZM, Wietgrefe S, et al. Propagation and dissemination of infection after vaginal transmission of simian immunodeficiency virus. J Virol. (2005) 79:9217–27. doi: 10.1128/JVI.79.14.9217-9227.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Salazar-Gonzalez JF, Salazar MG, Keele BF, Learn GH, Giorgi EE, Li H, et al. Genetic identity, biological phenotype, and evolutionary pathways of transmitted/founder viruses in acute and early HIV-1 infection. J Exp Med. (2009) 206:1273–89. doi: 10.1084/jem.20090378

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Mattapallil JJ, Douek DC, Hill B, Nishimura Y, Martin M, Roederer M. Massive infection and loss of memory CD4+ T cells in multiple tissues during acute SIV infection. Nature. (2005) 434:1093–7. doi: 10.1038/nature03501

PubMed Abstract | CrossRef Full Text | Google Scholar

33. ElKassar N, Gress RE. An overview of IL-7 biology and its use in immunotherapy. J Immunotoxicol. (2010) 7:1–7. doi: 10.3109/15476910903453296

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Kaleebu P, French N, Mahe C, Yirrell D, Watera C, Lyagoba F, et al. Effect of human immunodeficiency virus (HIV) type 1 envelope subtypes A and D on disease progression in a large cohort of HIV-1-positive persons in Uganda. J Infect Dis. (2002) 185:1244–50. doi: 10.1086/340130

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Carrington M, Nelson GW, Martin MP, Kissner T, Vlahov D, Goedert JJ, et al. HLA and HIV-1: heterozygote advantage and B^*35-Cw^*04 disadvantage. Science. (1999) 283:1748–52. doi: 10.1126/science.283.5408.1748

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Gao X, Nelson GW, Karacki P, Martin MP, Phair J, Kaslow R, et al. Effect of a single amino acid change in MHC class I molecules on the rate of progression to AIDS. N Engl J Med. (2001) 344:1668–75. doi: 10.1056/NEJM200105313442203

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Olvera A, Perez-Alvarez S, Ibarrondo J, Ganoza C, Lama JR, Lucchetti A, et al. The HLA-C^*04: 01/KIR2DS4 gene combination and human leukocyte antigen alleles with high population frequency drive rate of HIV disease progression. AIDS. (2015) 29:507–17. doi: 10.1097/QAD.0000000000000574

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Kiepiela P, Leslie AJ, Honeyborne I, Ramduth D, Thobakgale C, Chetty S, et al. Dominant influence of HLA-B in mediating the potential co-evolution of HIV HLA. Nature. (2004) 432:769–75. doi: 10.1038/nature03113

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Keating SM, Jacobs ES, Norris PJ. Soluble mediators of inflammation in HIV and their implications for therapeutics and vaccine development. Cytokine Growth Factor Rev. (2012) 23:193–206. doi: 10.1016/j.cytogfr.2012.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Neaton JD, Neuhaus J, Emery S. Soluble biomarkers and morbidity and mortality among people infected with HIV: summary of published reports from 1997 to 2010. Curr Opin HIV AIDS. (2010) 5:480–90. doi: 10.1097/COH.0b013e32833ed75d

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Borges AH, O'Connor JL, Phillips AN, Neaton JD, Grund B, Neuhaus J, et al. Interleukin 6 is a stronger predictor of clinical events than high-sensitivity C-reactive protein or D-dimer during HIV infection. J Infect Dis. (2016) 214:408–16. doi: 10.1093/infdis/jiw173

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Johansen C, Riis JL, Gedebjerg A, Kragballe K, Iversen L. Tumor necrosis factor alpha-mediated induction of interleukin 17C in human keratinocytes is controlled by nuclear factor kappaB. J Biol Chem. (2011) 286:25487–94. doi: 10.1074/jbc.M111.240671

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ramirez-Carrozzi V, Sambandam A, Luis E, Lin Z, Jeet S, Lesch J, et al. IL-17C regulates the innate immune function of epithelial cells in an autocrine manner. Nat Immunol. (2011) 12:1159–66. doi: 10.1038/ni.2156

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Peng T, Chanthaphavong RS, Sun S, Trigilio JA, Phasouk K, Jin L, et al. Keratinocytes produce IL-17c to protect peripheral nervous systems during human HSV-2 reactivation. J Exp Med. (2017) 214:2315–29. doi: 10.1084/jem.20160581

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Nies JF, Panzer U. IL-17C/IL-17RE: emergence of a unique axis in TH17 biology. Front Immunol. (2020) 11:341. doi: 10.3389/fimmu.2020.00341

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Chang SH, Reynolds JM, Pappu BP, Chen G, Martinez GJ, Dong C. Interleukin-17C promotes Th17 cell responses autoimmune disease via interleukin-17 receptor E. Immunity. (2011) 35:611–21. doi: 10.1016/j.immuni.2011.09.010

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Cocchi F, DeVico AL, Garzino-Demo A, Arya SK, Gallo RC, Lusso P. Identification of RANTES, MIP-1 alpha, and MIP-1 beta as the major HIV-suppressive factors produced by CD8+ T cells. Science. (1995) 270:1811–5. doi: 10.1126/science.270.5243.1811

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Jennes W, Sawadogo S, Koblavi-Deme S, Vuylsteke B, Maurice C, Roels TH, et al. Positive association between beta-chemokine-producing T cells and HIV type 1 viral load in HIV-infected subjects in Abidjan, Cote d'Ivoire. AIDS Res Hum Retroviruses. (2002) 18:171–7. doi: 10.1089/08892220252781220

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Cocchi F, DeVico AL, Yarchoan R, Redfield R, Cleghorn F, Blattner WA, et al. Higher macrophage inflammatory protein (MIP)-1alpha and MIP-1beta levels from CD8+ T cells are associated with asymptomatic HIV-1 infection. Proc Natl Acad Sci U S A. (2000) 97:13812–7. doi: 10.1073/pnas.240469997

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Aung H, McKenna SM, Ketoff NR, Jones L, Wu M, Hejal R, et al. Dysregulation of beta-chemokines in the lungs of HIV-1-infected patients. J Acquir Immune Defic Syndr. (2001) 26:305–14. doi: 10.1097/00126334-200104010-00002

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Freel SA, Lamoreaux L, Chattopadhyay PK, Saunders K, Zarkowsky D, Overman RG, et al. Phenotypic and functional profile of HIV-inhibitory CD8 T cells elicited by natural infection and heterologous prime/boost vaccination. J Virol. (2010) 84:4998–5006. doi: 10.1128/JVI.00138-10

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Walker WE, Kurscheid S, Joshi S, Lopez CA, Goh G, Choi M, et al. Increased levels of macrophage inflammatory proteins result in resistance to R5-tropic HIV-1 in a subset of elite controllers. J Virol. (2015) 89:5502–14. doi: 10.1128/JVI.00118-15

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kakkanaiah VN, Ojo-Amaize EA, Peter JB. Concentrations of circulating beta-chemokines do not correlate with viral load in human immunodeficiency virus-infected individuals. Clin Diagn Lab Immunol. (1998) 5:499–502. doi: 10.1128/CDLI.5.4.499-502.1998

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Ye P, Kazanjian P, Kunkel SL, Kirschner DE. Lack of good correlation of serum CC-chemokine levels with human immunodeficiency virus-1 disease stage and response to treatment. J Lab Clin Med. (2004) 143:310–9. doi: 10.1016/j.lab.2004.01.012

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Ramalingam S, Kannangai R, Abraham OC, Subramanian S, Rupali P, Pulimood SA, et al. Chemokine profile among human immunodeficiency virus-1 (HIV-1) infected individuals from southern India. Indian J Med Res. (2008) 127:133–9.

PubMed Abstract | Google Scholar

56. Roberts L, Passmore JA, Williamson C, Little F, Bebell LM, Mlisana K, et al. Plasma cytokine levels during acute HIV-1 infection predict HIV disease progression. Aids. (2010) 24:819–31. doi: 10.1097/QAD.0b013e3283367836

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Huang X, Liu X, Meyers K, Liu L, Su B, Wang P, et al. Cytokine cascade and networks among MSM HIV seroconverters: implications for early immunotherapy. Sci Rep. (2016) 6:36234. doi: 10.1038/srep36234

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Jacobs ES, Keating SM, Abdel-Mohsen M, Gibb SL, Heitman JW, Inglis HC, et al. Cytokines elevated in HIV elite controllers reduce HIV Replication in vitro and modulate HIV restriction factor expression. J Virol. (2017) 91:e02051–16. doi: 10.1128/JVI.02051-16

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Fischer J, Jung N, Robinson N, Lehmann C. Sex differences in immune responses to infectious diseases. Infection. (2015) 43:399–403. doi: 10.1007/s15010-015-0791-9

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Takahashi T, Ellingson MK, Wong P, Israelow B, Lucas C, Klein J, et al. Sex differences in immune responses that underlie COVID-19 disease outcomes. Nature. (2020) 588:315–20. doi: 10.1038/s41586-020-2700-3

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Saeki H, Tamaki K. Thymus and activation regulated chemokine (TARC)/CCL17 and skin diseases. J Dermatol Sci. (2006) 43:75–84. doi: 10.1016/j.jdermsci.2006.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Vestergaard C, Deleuran M, Gesser B, Larsen CG. Thymus- and activation-regulated chemokine (TARC/CCL17) induces a Th2-dominated inflammatory reaction on intradermal injection in mice. Exp Dermatol. (2004) 13:265–71. doi: 10.1111/j.0906-6705.2004.00149.x

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Kanki PJ, Hamel DJ, Sankale JL, Hsieh C, Thior I, Barin F, et al. Human immunodeficiency virus type 1 subtypes differ in disease progression. J Infect Dis. (1999) 179:68–73. doi: 10.1086/314557

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Rubbert A, Weissman D, Combadiere C, Pettrone KA, Daucher JA, Murphy PM, et al. Multifactorial nature of noncytolytic CD8+ T cell-mediated suppression of HIV replication: beta-chemokine-dependent and -independent effects. AIDS Res Hum Retroviruses. (1997) 13:63–9. doi: 10.1089/aid.1997.13.63

PubMed Abstract | CrossRef Full Text | Google Scholar