1. Introduction
Human endogenous retroviruses (HERVs) are remnants of ancient viral infections integrated into the human genome throughout evolutionary history. These genetic elements are estimated to comprise nearly 8% of the human genome and are passed from generation to generation [
1]. Although most HERVs are now functionally inactive due to mutations and deletions, some have been implicated in disease processes, such as autoimmune and neurodegenerative (such as amyotrophic lateral sclerosis – ALS) disorders, several types of cancer and immune dysregulation and coagulopathy in critical COVID-19 [
2,
3,
4,
5,
6,
7,
8,
9].
HERV-K represents the latest classified family, encompassing the most recently integrated HERV groups within the human genome [
10]. HERV-Ks can be further divided into ten subfamilies (HML 1-10), with some retaining their structural integrity and ability to generate viral proteins [
11]. Importantly, their high sequence similarity with contemporary retroviruses, such as HIV-1, has led to investigations into potential roles in HIV infection and pathogenesis. HIV-1 can activate HERV-K (HML-2) mRNA and protein expression
in vitro [
12,
13,
14,
15]. Moreover, untreated PLWH display higher levels of HERV-K RNA in peripheral blood mononuclear cells (PBMCs) [
4,
13,
16] or plasma [
17,
18,
19] and higher titers of anti-HERV-K (HML-2) antibodies and HERV-specific CD8
+ T-cell responses [
19,
20,
21,
22] compared to antiretroviral therapy (ART)-treated individuals and/or HIV-uninfected healthy subjects. However, the exact nature of the relationship between HERV-K and HIV-1 expression and the pathogenic impact of HERV-K activation in the setting of HIV-1 infection remains controversial.
Some findings support that expression of HERV-K during HIV-1 infection could directly affect immunopathogenesis. One study showed that the anti-HERV-K (HML-2) transmembrane (TM) antibody response was positively correlated with HIV plasma viremia and inversely correlated with CD4
+ T-cell counts in viremic noncontrollers, suggesting that anti-TM antibodies could target HERV-K (HML-2) Env TM-expressing HIV-1-infected CD4
+ cells and accelerate CD4
+ T-cell depletion [
13,
17,
22]. Interestingly, the HERV-K Env TM protein contains an immunosuppressive domain that modulates the release of different cytokines by PBMCs
in vitro, including immunosuppressive IL-10, and may thus contribute to suppressing immune responses during HIV-1 infection [
15]. Other studies propose that HERV-K expression could also exacerbate HIV-1 infection by stimulating the innate immune system and the expression of type I interferon in macrophages and probably other immune cells and, therefore, enhancing proinflammatory signaling in a harmful way [
23].
Other findings, by contrast, support that increased expression of HERV-K and an enhanced immune response against HERV-K proteins may improve virologic control in HIV-1 infection. The induction of HERV-K antigen expression in HIV-1-infected cells could stimulate a HERV-K–specific cytotoxic CD8
+ T-cell response, resulting in the targeted elimination of reservoir HIV-infected cells [
15]. Some studies show an inverse correlation between HIV-1 plasma viral load and HERV-specific CD8
+ T-cell responses or antibodies against HERV-K Gag and further observed that HIV controllers (HICs), a group of rare individuals who maintain undetectable viral loads without ART, have more robust cellular and antibody responses against HERV-K than ART-suppressed, viremic noncontrollers, and uninfected subjects [
19,
20,
21,
24]. Moreover, the level of HERV-K RNA expression in PBMCs was negatively associated with the level of T-cell activation in untreated PLWH, suggesting that HERV-K expression does not contribute to pathologic immune activation in chronic HIV-1 infection [
16].
Previous studies revealed that LTR-directed transcription of both HERV-K (HML-2) and HIV-1 proviruses could be regulated by similar factors. One study demonstrated that HIV-1 Tat plays an important role in activating the expression of HERV-K (HML-2) in the setting of HIV-1 infection [
25]. Activation of HERV-K (HML-2) transcription by Tat occurs in a partially different way than that of HIV-1 because it involves Tat’s interaction with cellular factors NF-kB and NF-AT but not with Cyclin T1 or P-TEFb. Another study with sorted PBMCs showed that increased HERV-K (HML-2) RNA expression in the setting of HIV-1 infection resulted from multiple cell sources, including those not infected by HIV-1 [
4]. Thus, HIV-1 infection may also enhance HERV-K (HML-2) transcription in both infected and uninfected cells by indirect mechanisms that probably modulate the expression of key transcription factors, such as increased immune activation and inflammation [
23] and/or mechanisms influencing the cellular epigenetic environment [
26].
In a previous study conducted by our group, we discovered that the transcripts of two multifunctional cellular proteins, cyclin-dependent kinase (CDK) inhibitor 1A (CDKN1A/p21) and monocyte chemotactic protein-induced protein 1 (MCPIP1), were significantly elevated in HICs compared to ART-suppressed and HIV-1 negative control individuals [
27]. The p21 can block HIV-1 by suppressing the biosynthesis of deoxyribonucleoside triphosphates (dNTPs), which are essential for viral reverse transcription, and by inhibiting the activity of CDK9, which is crucial for HIV-1 viral mRNA transcription [
28,
29,
30,
31,
32]. MCPIP1 is critical in controlling inflammatory responses and maintaining immune homeostasis. It can block HIV-1 replication by promoting the degradation of several inflammatory cytokine mRNAs and viral mRNAs via its RNase activity [
33,
34]. Interestingly, p21 and MCPIP1 contribute to the limited activation of macrophages and T cells by downregulating NF-kB activity [
34,
35,
36,
37,
38,
39], which may also indirectly reduce the LTR-driven activation of HIV-1 and HERV-K proviruses.
We hypothesize that the upregulation of cellular restriction factors, such as p21 and MCPIP1, that negatively modulate HIV-1 transcription in HICs might also negatively affect the LTR-directed transcription of endogenous retroviruses. To test this hypothesis, we measured the ERVK-6 RNA expression levels (envelope protein fragment, also called HERV-K (HML-2) or HERV-K_7p22.1) in PBMCs of HIV controllers (n = 23), HIV-1 ART-suppressed subjects (n = 10), and HIV-1-negative individuals (n = 8) and correlated the transcript expression level of HERV-K with multiple HIV restriction factors.
3. Results
3.1. ERVK-6 mRNA levels are downregulated in HICs
We first evaluated the ERVK-6 transcript levels in the HIC and ART groups compared to the NEG group. We found a significant (P < 0.05) downregulation of ERVK-6 mRNA in the HIC (0.88 - mean fold-change), VC (0.84 - mean fold-change) and EC subgroups (0.86 - mean fold-change), while no significant differences in ERVK-6 expression levels were observed between the NEG and ART groups (0.94 - mean fold-change, P = 0.33) (
Figure 1). In comparison to the ART group, we found a significant downregulation of ERVK-6 mRNA in VC (0.89-mean fold-change, P = 0.02) and a lower, but not significant, expression in EC (0.94-mean fold-change, P = 0.70) and HIC (0.92-mean fold-change, P = 0.08) (
Figure 1). No significant differences in ERVK-6 expression levels were observed between the EC and VC subgroups (
P = 0.08) (
Figure 1).
3.2. ERVK-6 mRNA and cellular restriction factor levels are correlated
We next assessed whether reduced ERVK-6 transcripts in HICs could be associated with mRNA expression levels of 13 cellular restriction factors that act against HIV-1. Pearson's correlation of all groups combined showed that the mRNA expression of ERVK-6 was negatively correlated with both p21 (r = -0.48;
P = 0.0013) and MCPIP1 (r = -0.35;
P = 0.0248) (
Figure 2A-B), positively correlated with SERINC3 (r = 0.58;
P < 0.0001), SERINC5 (r = 0.60;
P < 0.0001), APOBEC3G (r = 0.43;
P = 0.005), IFITM2 (r = 0.43;
P = 0.004), SAMHD1 (r = 0.43;
P = 0.005), and SLFN11 (r = 0.40;
P = 0.009) (
Figure 2C-H) and was not significantly correlated (
P > 0.05) with the other five restriction factors (Mx1, Mx2, IFITM1, IFITM3, and Tetherin) analyzed (
Figure S1).
We also evaluated the correlation coefficient between restriction factors significantly correlated with ERVK-6 mRNA expression levels. We detected a significant negative correlation between p21 and SERINC3 (r = -0.53;
P = 0.0004) and APOBEC3G (r = - 0.40;
P = 0.009) (
Figure 3A). Significant positive correlations were observed between SERINC3 and APOBEC3G (r = 0.49;
P = 0.0013) and between SERINC5 and SLFN11 (r = 0.37;
P = 0.018), SAMHD1 (r = 0.45;
P = 0.003), and IFITM2 (r = 0.34;
P = 0.028). Finally, we also observed positive correlations between SLFN11 and APOBEC3G (r = 0.41; P = 0.008) (
Figure 3A).
Multivariate principal component analysis (PCA) comprising p21, MCPIP1, SERINC3, and SERINC5 variables is shown in
Figure 3B and 3C according to either their mRNA ERVK-6 expression levels divided into four different quartiles (
Figure 3C) or their group (
Figure 3D). Our results revealed that 80.6% of the total variance in response to the four restriction factors was expressed by two principal components. The first component (Dim1) represented 47.2%, while the second (Dim2) represented 33.4% of the total variance. The variables p21 and MCPIP1 contributed most to the first component (Dim1), while SERINC5 and SERINC3 contributed most to the second component (Dim2).
In general, we observed a strong association between the ERVK-6 expression quartiles and the HIV-1 control groups (X2 = 32.293, df = 9, P = 0.0001772). This association was even more evident between the NEG control group and the quartiles with the highest ERVK-6 expression (i.e., 7 in the 4th quartile, 2 in the 3rd quartile, 1 in the 2nd quartile, and 0 in the 1st quartile) (X2 = 16.056, df = 3, P = 0.001104). The opposite situation, also highly associated, was observed between the VC group and ERVK-6 expression (7 in the 1st quartile, 1 in the 2nd quartile, 2 in the 3rd quartile, and 0 in the 4th quartile) (X2 = 13.641, df = 3, P = 0.003437). This association was confirmed when examining the HIC group (EC + VC), i.e., 10 in the 1st quartile, 8 in the 2nd quartile, 4 in the 3rd quartile, and only 1 in the 4th quartile (X2 = 17.413, df = 3, P = 0.0005812).
The separation between the HIC (EC + VC) and control (NEG + ART) groups occurred almost perfectly along the first component (Dim1), with the first group to the right of the origin (1st and 4th quadrants) and the second to the left of the origin (2nd and 3rd quadrants), except for two EC and one VC individuals, in the 3rd and 2nd quadrants, respectively; two of them had low ERVK-6 expression (2nd quartile), and another had very low ERVK-6 expression (1st quartile). The same observation regarding the first component (Dim1) was made with the levels of ERVK-6 expression, with the individuals with the lowest expression (1st and 2nd quartile) being further to the right of the origin (1st and 4th quadrants) and those with the highest expression to the left of the origin (2nd and 3rd quadrants). Notably, the individuals with the highest ERVK-6 expression (4th quartile) were concentrated in the 2nd quadrant of the PCA, coinciding with the predominance of healthy control individuals (NEG), while the individuals with the lowest ERVK-6 expression (1st quartile) were concentrated in the 4th quadrant of the PCA, coinciding with the predominance of VC.
4. Discussion
Previous studies revealed that HERV-K (HML-2) RNA expression in PBMCs from untreated viremic-PLWH was higher than that in ART-treated PLWH and/or HIV-uninfected healthy subjects [
4,
13,
16] and that HICs have more robust cellular and antibody responses against HERV-K than ART-suppressed, viremic noncontrollers and uninfected subjects [
20,
21]. This led to the hypothesis that HERV-K expression in HICs could expose the immune system to HERV-K antigens that act as immutable targets for cytotoxic lysis, resulting in the virologic control of HIV-1. However, none of these studies analyzed HERV-K (HML-2) RNA expression in HICs. Our study reveals that HERV-K 6 (HML-2) RNA expression in PBMCs from HICs, particularly the VC subgroup, was significantly downregulated compared with both ART and NEG control groups.
We hypothesize that HICs may activate antiviral mechanisms that control the expression of both exogenous and endogenous retroviruses. Consistent with this notion, we detected that ERVK-6 RNA levels in PBMCs were negatively correlated with the expression levels of two cellular restriction factors, namely, p21 and MCPIP1. These proteins can block HIV-1 replication in macrophages and CD4
+ T cells [
29,
30,
33,
47,
48] and further limit aberrant immune activation [
34,
49,
50]. Interestingly, p21 and MCPIP1 limit the activation of macrophages and T cells by downregulating the activity of NF-κB [
34,
35,
36,
37,
38,
39], a transcription factor that stimulates the LTR-driven transcription of HIV-1 and HERV-K proviruses [
51,
52]. In a previous study, we showed that MCPIP1 and p21 mRNA and protein expression levels are upregulated in PBMCs from our HIC cohort [
27]. Therefore, we suggest that HICs can activate a homeostatic anti-inflammatory response that comprises antiviral factors, such as p21 and MCPIP1, to prevent excessive immune activation driven by residual HIV-1 replication. This negative homeostatic response inhibits the NF-κB pathway, which may reduce the efficiency of the LTR-driven transcription of HIV-1 and HERV-K proviruses.
Interestingly, we detected that ERVK-6 RNA expression in our cohort was positively correlated with two cellular restriction factors, SERINC3 and SERINC5, that may upregulate HERV-K (HML-2) RNA levels by activating the NF-κB pathway. SERINC3/5 were initially identified as cell restriction factors that can potently suppress HIV-1 infectivity by incorporating into budding viral particles and impairing subsequent virion fusion and infection of new target cells [
53,
54]. A recent study revealed that SERINC3/5 exhibited additional antiviral activities by forming a signaling complex with the mitochondrial antiviral signaling protein (MAVS) and TRAF6 at the mitochondria that increased the phosphorylation and activation of IRF3 and IκBα, thus cooperatively enhancing the NF-κB inflammatory pathway and type I IFN (IFN-I) production [
55]. Moreover, we detected a significant negative correlation between the expression of p21 and SERINC3, suggesting that these restriction factors may display antagonist transcriptional regulation and may be part of the mechanism that controls the repression and activation of HERVs.
Our findings also revealed that the expression of some ISGs analyzed here (APOBEC3G, SLFN11, and SAMDH1) was positively correlated with the expression of both ERVK-6 and SERINC3/5. Cytosolic HERV-K dsRNA/cDNA may activate different nucleic acid sensors of innate immunity, such as retinoic acid-inducible gene-I-like receptors (RIG-I and MDA5) [
56,
57,
58] and cGMPAMP synthase (cGAS) [
59]. Activation of those innate immune sensors leads to phosphorylation and activation of TANK-binding kinase 1 (TBK1) and IκB kinase-ϵ (IKKϵ) kinases, which in turn activate the transcription factors IRF3 and NF-κB, ultimately leading to upregulation of type I IFNs, ISGs, and other proinflammatory cytokines. Interestingly, SAMHD1 and SLFN11 belong to a subset of ISGs that may be directly induced by IRF3 [
60,
61], and APOBEC3G may be positively modulated by the NF-κB pathway [
62]. We propose that increased levels of SERINC3/5 and cytosolic HERV-K RNA/DNA may directly upregulate some ISGs by triggering the activation of IRF-3 and NF-κB, independent of the engagement of the type I IFN-JAK-STAT pathway.
The generation and maintenance of strong HERV-specific cellular and antibody responses in HICs seem inconsistent with the very low ERVK-6 mRNA expression detected here in PBMCs. Some findings, however, may explain these apparently conflicting results. Previous studies conducted in our and other cohorts of HICs showed that HIV-1 continues to replicate and evolve despite undetectable or extremely low levels of viremia [
41,
63,
64,
65,
66,
67,
68]. Moreover, HIV-1 replication in HICs seems to mainly occur in the lymph nodes [
68]. Thus, HERV-K expression in HICs may be mostly restricted to HIV-infected cells residing in lymph nodes, and expression of both HERV-K and HIV-1 antigens at those sites may stimulate the continuous generation of HERV-specific cellular and antibody responses. In contrast, most PBMCs and HIV-infected cells in the periphery display an antiviral transcriptional signature that limits the activation of both exogenous and endogenous retroviruses.
Persistent overexpression of IFN-α/β and ISGs is a hallmark of chronic HIV-1 infection and is associated with detrimental immune activation, bystander CD4
+ T-cell apoptosis, and progressive disease [
69,
70]. Some data suggest that augmented HERV expression may play an active role in exacerbating and perpetuating chronic inflammation in some autoimmune diseases via type I IFN and MAVS positive signaling feedback loops [
71]. On the other hand, HERV-K (HML-2) knockdown significantly downregulated genes containing interferon-stimulated response elements (ISREs) in their promoters in basal and IFN-γ-challenged macrophages and had implications for the paracrine activation of nearby cells following macrophage activation [
23]. Moreover, repression of endogenous retrovirus
in vivo alleviates tissue aging and, to some extent, organismal aging by attenuating innate immune responses [
59]. These findings suggest that repression of HERV-K expression in HICs may contribute to attenuating the activation of innate immune responses, preventing chronic overexpression of IFN-α/β/ISGs, and preserving immune homeostasis in those subjects.
Previous studies from our cohort and other groups showed that HICs displayed elevated levels of some key plasma inflammatory markers, such as sCD14, IP-10, IL-18, and D-dimer, and higher CD8
+ T-cell activation with respect to HIV-negative individuals [
40,
72,
73,
74,
75,
76,
77,
78]. Moreover, VCs displayed higher levels of immune activation and residual inflammation than ECs. The proinflammatory profile observed in HICs, however, is different from the aberrant inflammation observed in untreated noncontrollers [
79,
80]. Moreover, another study that performed gene expression analysis of PBMCs and purified CD4
+ and CD8
+ T cells from HICs showed that these individuals, compared to ART-treated individuals, have lower expression of several inflammatory genes, such as IL1A/B, IL-6, CXCL5, CXCL13, and CXCL1, but maintain upregulation of genes associated with cytotoxicity and the T-cell response [
81]. We support that the low levels of HERV-K 6 (HML-2) detected here in HICs may be linked to the upregulation of an anti-inflammatory response that limits aberrant inflammation while preserving efficient antiviral T-cell responses. Such homeostatic anti-inflammatory mechanisms would be more critical for VCs due to the higher levels of chronic HIV-1 antigenic stimulation with respect to ECs.