| Literature DB >> 33934119 |
Abstract
Immunogenetic studies in the past three decades have uncovered a broad range of human genetic factors that seem to influence heterosexual HIV-1 transmission in one way or another. In our own work that jointly evaluated both genetic and nongenetic factors in two African cohorts of cohabiting, HIV-1-discordant couples (donor and recipient pairs) at risk of transmission during quarterly follow-up intervals, relatively consistent findings have been seen with three loci (IL19, HLA-A, and HLA-B), although the effect size (i.e., odds ratio or hazards ratio) of each specific variant was quite modest. These studies offered two critical lessons that should benefit future research on sexually transmitted infections. First, in donor partners, immunogenetic factors (e.g., HLA-B*57 and HLA-A*36:01) that operate directly through HIV-1 viral load or indirectly through genital coinfections are equally important. Second, thousands of single-nucleotide polymorphisms previously recognized as "causal" factors for human autoimmune disorders did not appear to make much difference, which is somewhat puzzling as these variants are predicted or known to influence the expression of many immune response genes. Replicating these observations in additional cohorts is no longer feasible as the field has shifted its focus to early diagnosis, universal treatment, and active management of comorbidities.Entities:
Mesh:
Year: 2021 PMID: 33934119 PMCID: PMC8225584 DOI: 10.1038/s41435-021-00130-y
Source DB: PubMed Journal: Genes Immun ISSN: 1466-4879 Impact factor: 2.676
Figure 1.Several scenarios in the setting of heterosexual HIV-1 transmission among discordant couples.
When two main factors in the HIV-1-infected source partners (SPs) are assessed along with two well-documented factors in the recipient partners (also known as HIV-1-exposed, seronegative individuals or HESNs), couples can be divided into various composite risk groups. The group with extremely low risk (having all factors indicated by arrows) is expected to obscure the search for immunogenetic determinants, especially when duration of follow-up is short. GUI, genital ulcer/inflammation (clinical manifestations of genital co-infections); MMC, male medical circumcision (rates of uptake <5% in our cohorts).
Figure 2.Kaplan-Meier curves showing the impact of HIV-1 viral load (VL) in donor/index partners on heterosexual HIV-1 transmission among 420 HIV-1 discordant Zambian couples with longitudinal follow-up.
The three VL categories followed earlier strategies [83], and subjects with low VL (<10−4 RNA copies/ml of plasma) were treated as the reference (ref.) group. The estimates of relative hazards (RH) and 95% confidence interval (CI) were based on a Cox proportional hazards model. Graph is redrawn from data published in 2008 [17]. Couples that remained eligible for further assessment at various follow-up intervals are indicated below the curves.
Assembly of two study cohorts for immunogenetic studies.[a]
| Overall characteristics | Lusaka, Zambia | Kigali, Rwanda |
|---|---|---|
| Participants | HIV-1 discordant couples | HIV-1 discordant couples |
| Enrollment dates | 1995–2011 | 2002–2011 |
| Intervention | CVCT & cART | CVCT & cART |
| Follow-up visits | Quarterly | Quarterly |
| Genital co-infections[ | GUI | GUI |
| Total enrollment | 4,759 couples | 1,748 couples |
| Subset for immunogenetic studies[ | 994 pairs (21%) | 361 pairs (21%) |
Whenever applicable, data from other studies pertinent to HIV-1 transmission in African countries are also discussed in the text.
Genital ulcer/inflammation (GUI) was frequent enough for ancillary studies. Other abbreviations: CVCT, couple-based voluntary counseling and testing; cART, combination antiretroviral therapy.
Selection criteria: (i) ages between 18 at enrollment and 65 at end of follow-up, (ii) >12 months of follow-up before therapy or drop-out (while transmission-free), (iii) donor/source partner viral load >2,000 RNA copies/mL, and (iv) at least one known risk factor (e.g., unprotected sex or pregnancy) for transmission during study intervals.
Summary of consensus immunogenetic findings from two study cohorts.
| Locus | Variant (allele or aa residue) | Setting | Association | Reference(s) |
|---|---|---|---|---|
| rs12407485-A[ | HESNs | Resistance to HIV-1 acquisition | [ | |
| A*68:02[ | HESNs | Susceptibility to HIV-1 acquisition | [ | |
| B*57[ | SPs | HIV-1 transmission | [ | |
| Alleles carrying P2-Met[ | HESNs | Susceptibility to HIV-1 acquisition | [ |
Within a predicted enhancer region. According to a recent (March 2021) search in PhenoScanner (version 2), rs12407485 is a well-recognized eQTL, being associated with the expression of four genes (CD55, IL10, IL19 and IL24).
Confirmative results from the Rwandan cohort have not been published.
Mostly B*57:03 and often reflected by its persistent impact on viral load
Methionine at position 2 of a leader peptide that binds to HLA-E; also supported by in vitro assays [6].
Cohort-specific immunogenetic findings with various supporting evidence.
| Locus | Variant (allele or aa residue) | Setting | Association | Reference(s) |
|---|---|---|---|---|
| A*36:01[ | Zambian SPs | Promoting HIV-1 transmission | Ref. [ | |
| Alleles sharing[ | Both partners | Promoting HIV-1 transmission | Ref. [ | |
| Full-length allele *001[ | Zambian SPs | Promoting HIV-1 transmission | Ref. [ |
Observation was confirmed multiple times in expanded cohorts (e.g., Figure 2), with clear impact on viral load in the SP partners.
Supported by experimental evidence [8].
The only allele encoding a functional, full-length receptor, with implications for inflammation mediated through natural killer cells [7].
Figure 3.Kaplan-Meier curves showing the impact of HLA-A*36:01 in donor/index partners on heterosexual HIV-1 transmission and genital ulcer/inflammation (GUI) during longitudinal follow-up.
Prompted by our original report on 429 HIV-1 discordant Zambian couples [17], several analyses of expanded datasets (530 couples in panel a versus 639 couples in panel b) confirmed that HLA-A*36:01 in donor/index partners was persistently an unfavorable factor for HIV-1 transmission. This association was partially explained by elevated HIV-1 viral load setpoint in donor partners, as well as presence of GUI over time (panel c, 530 couples), the association with GUI was mostly seen in female-to-male transmission (panel d, with 273 couples) (combined data from two publications [17, 19], including recent updates). The estimates of relative hazards (RH) and 95% confidence interval (CI) were based on Cox proportional hazards models.
Multivariable models that jointly assess the relative impact of three risk factors (host and viral) on heterosexual HIV-1 transmission among Zambian couples.
| Factors in
model[ | All subjects (638 couples) | MTF subset (312 couples) | FTM subset (316 pairs) | |||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| RH | 95% CI | RH | 95% CI | RH | 95% CI | |||||||
| Donor VL[ | 634 | 1.8 | 1.6–2.2 | <0.0001 | 312 | 1.67 | 1.3–2.2 | <0.0001 | 316 | 1.9 | 1.5–2.5 | <0.0001 |
| A*36:01 (donor) | 63 | 1.6 | 1.2–2.2 | 0.004 | 28 | 1.44 | 0.9–2.3 | 0.130 | 35 | 1.8 | 1.2–2.9 | 0.009 |
| A*68:02 (recipient)[ | 99 | 1.5 | 1.1–2.0 | 0.006 | 58 | 1.51 | 1.0–2.2 | 0.029 | 41 | 1.5 | 0.9–2.4 | 0.099 |
Data from an expanded Zambian cohort (638 couples) (combined data from two publications [17, 19], including recent updates). Results were consistent with earlier findings based on analyses of paired recipient and index/donor partners (each couple was counted as one unit) [17], even when the cohort was split by the direction of transmission (male-to-female/MTF versus female-to-male/FTM).
HIV-1 viral load (VL) was treated as a continuous variable after log10-transformation. The summary statistics, including relative hazards (RH) and 95% confidence interval (CI), are based on Cox proportional hazards models (adjusted for factors in each model).
Data missing in several subjects. Of note, two earlier reports have portrayed the A2 supertype, including A*68:02, as protective against HIV-1 acquisition in cohorts from Kenya [81, 82].
Figure 4.A 3D graph showing the clustering of subjects from three countries.
In principal component analyses using genome-wide SNPs with low linkage disequilibrium (pairwise r2 <0.20), the first three dimensions (PC1-PC3) captured most of the genetic variations for subjects enrolled from Kigali (Rwanda), Lusaka (Zambia) and Birmingham (United States). The Birmingham cohort consists of African Americans and European Americans. Variance explained by each component is shown in parentheses (Tang et al., unpublished data).