Integrated analysis of lncRNA, miRNA and mRNA profiles reveals potential lncRNA functions during early HIV infection
Journal of Translational Medicine volume 19, Article number: 135 (2021)
Long noncoding RNAs (lncRNAs) can regulate gene expression in a cis-regulatory fashion or as “microRNA sponges”. However, the expression and functions of lncRNAs during early human immunodeficiency virus (HIV) infection (EHI) remain unclear.
3 HAART-naive EHI patients and 3 healthy controls (HCs) were recruited in this study to perform RNA sequencing and microRNA (miRNA) sequencing. The expression profiles of lncRNAs, mRNAs and miRNAs were obtained, and the potential roles of lncRNAs were analysed based on discovering lncRNA cis-regulatory target mRNAs and constructing lncRNA–miRNA–mRNA competing endogenous RNA (ceRNA) networks. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on 175 lncRNA-associated differentially expressed (DE) mRNAs to investigate the potential functions of DE lncRNAs in ceRNA networks.
A total of 242 lncRNAs, 1240 mRNAs and 21 mature known miRNAs were determined as differentially expressed genes in HAART-naive EHI patients compared to HCs. Among DE lncRNAs, 44 lncRNAs were predicted to overlap with 41 target mRNAs, and 107 lncRNAs might regulate their nearby DE mRNAs. Two DE lncRNAs might regulate their cis-regulatory target mRNAs BTLA and ZAP70, respectively, which were associated with immune activation. In addition, the ceRNA networks comprised 160 DE lncRNAs, 21 DE miRNAs and 175 DE mRNAs. Seventeen DE lncRNAs were predicted to regulate HIF1A and TCF7L2, which are involved in the process of HIV-1 replication. Twenty DE lncRNAs might share miRNA response elements (MREs) with FOS, FOSB and JUN, which are associated with both immune activation and HIV-1 replication.
This study revealed that lncRNAs might play a critical role in HIV-1 replication and immune activation during EHI. These novel findings are helpful for understanding of the pathogenesis of HIV infection and provide new insights into antiviral therapy.
Early human immunodeficiency virus infection (EHI) reflects the period following viral entry during which viremia bursts occur until decreasing to a stable viral load level approximately 6 months post infection [1,2,3,4,5]. During EHI, viremia increases exponentially, while proinflammatory cytokines are produced by innate immune cells and coincide with mounting adaptive immune responses. Meanwhile, the viral reservoir is formed early in infection. It has been reported that the events occurring during EHI are critical in determining the transmission rates, the course of disease progression, and HIV-related morbidity and mortality . To date, the pathogenesis of EHI remains unclear. Therefore, exploring the mechanism of HIV pathogenesis during EHI is helpful for the design of therapeutic strategies and vaccine development.
Long noncoding RNAs (lncRNAs), which are longer than 200 nt and do not encode proteins, have attracted much attention in recent years [7, 8]. LncRNAs can act as guides, scaffolds, decoys, signalling molecules or via genomic targeting, acting as cis-or trans-regulatory elements or through antisense interference to play vital roles in a number of biological processes [9,10,11,12]. Moreover, lncRNAs bind to specific miRNAs and regulate their functions by acting as “miRNA sponges”, which are known as competing endogenous RNAs (ceRNAs) [13, 14]. For example, lncRNA MT1JP regulates the expression of FBXW7 by competitively binding to miR-92a-3p, which in turn affects the progression of gastric cancer . Therefore, it will be of interest to explore the functional roles of lncRNAs in different research fields.
Several HIV-related transcriptome analyses have focused on the potential functions of lncRNAs, such as the differential expression of lncRNAs during HIV-1 and HIV-2 infection , the impact of lncRNAs on HIV replication , the roles of lncRNAs in the establishment and reversal of latent HIV infection [18, 19] and the relationships between lncRNAs and disease progression . For example, lnc-RNF125-1, lnc-TRAF5-1 and lnc-TRDMT1-1 might control HIV replication by regulating the expression and function of nearby mRNAs in a cis-regulatory fashion . Furthermore, 3 lncRNAs, PVT1, RP11347C18.3, and RP11-539 L10.2, are dysregulated genes in HIV latency and might be targets for HIV latency reversal . However, the potential roles of lncRNAs during EHI have not been systematically investigated.
In this study, we analysed the lncRNA, miRNA, and mRNA expression profiles from 3 HAART-naive EHI patients and 3 healthy controls (HCs) using RNA sequencing (RNA-Seq) and microRNA (miRNA) sequencing (miRNA-Seq). The potential roles of lncRNAs were analysed based on the prediction of lncRNA cis-regulatory target mRNAs and by constructing lncRNA–miRNA–mRNA ceRNA networks. We found several DE lncRNAs might involve in HIV-1 replication and immune activation during EHI. Our novel findings will benefit further exploration of the pathogenesis of HIV and propose new potential targets for therapeutic strategies.
Materials and methods
This study was approved by the Medical Research Ethics Committee of the First Affiliated Hospital of China Medical University (Shenyang, China). All participants provided written informed consent for participation in this study.
The HAART-naive EHI patients in our study were recruited from a large-scale, prospective high-risk men who had sex with men (MSM) cohort in Liaoning, China. In this cohort, HIV-1-negative MSM were followed up every 8 weeks. A fourth-generation enzyme-linked immunosorbent assay (ELISA) was used to detect HIV-1 infection status, and western blotting was applied for further validation. HIV-1 RNA amplification testing was performed on ELISA-negative samples. HIV infection time was estimated as the previous 14 days from the date on which the sample testing result was RNA-positive and ELISA-negative or the midpoint between the last negative and first positive results of the ELISA screening tests . The stage within 180 days post HIV infection was defined as EHI.
Three HAART-naive EHI patients and 3 HCs were enrolled for RNA-seq and miRNA-seq. For the 3 EHI patients, the level of CD4+ T cell counts and viral loads (VLs) were similar. The 3 EHI patients and 3 HCs were recruited from the same high-risk MSM cohort mentioned above and all participants were males with similar ages and none of them were infected by herpes simplex virus (HSV), syphilis, hepatitis B virus (HBV) or hepatitis C virus (HCV). Also, none of them were in use of alcohol, recreational drugs or nicotine products. The clinical characteristics of subjects were listed in Additional file 1: Table S1. In addition, we performed RT-qPCR validation on 10 EHI patients and 7 HCs, whom were recruited from the same cohort as mentioned above. The age, lifestyle and co-morbidities status were no statistically different between EHI patients and HCs (Additional file 1: Table S1).
Peripheral blood mononuclear cells (PBMCs) from EHI patients and HCs were isolated by density gradient centrifugation with Ficoll–Hypaque (GE Healthcare), cryopreserved in foetal calf serum (Gibco) supplemented with 10% DMSO (Sigma) and stored in liquid nitrogen.
Anticoagulated whole blood of HIV-1 infected patients was collected at each follow-up visit after they were diagnosed as HIV-1 infection. The CD4+T cell count was determined using FACS Calibur™ flow cytometer (Becton–Dickinson) and plasma viral load (VL) was detected by the COBAS AmpliPrep/COBAS TaqMan HIV-1 test (Roche). A Syphilis screening serological test was performed with the Rapid Plasma Reagin (RPR) test (Shanghai Kehua Bio-engineering Co). The Treponema Pallidum particle (TPPA) assay (Serodia TPPA) was used to confirm the positive RPR results. Subjects with both the TPPA and RPR plasma positive were considered as current infection. Herpes simplex virus type 2 (HSV-2)-specific immunoglobulin G (IgG) antibody testing was performed to detect HSV-2 infection using ELISA methods (HerpeSelect-2 ELISA IgG; Focus Diagnostics). Wondfo HIV-HCV-TP-HBsAg Multi-Test Kits (Colloidal Gold) were used to detect HBV and HCV infection.
RNA isolation and quality control
Total RNA was isolated from PBMCs of 3 HAART-naive EHI patients and 3 HCs using TRIzol reagent (Life Technologies). TurboDnase (Life Technologies) was used to eliminate DNA. RNA concentrations were measured with a NanoDrop ND-2000 instrument (Thermo Fisher Scientific), and RNA quality was estimated by denaturing agarose gel electrophoresis for subsequent sequencing. A NanoDrop ND-1000 instrument (Thermo Fisher Scientific) was to ensure quality control of miRNA concentrations.
Ribo-Zero rRNA removal kits (Illumina) were used to exclude ribosomal RNA molecules (rRNA). RNA sequencing libraries were constructed with the TruSeq Stranded Total RNA Library Prep Kit (Illumina) following the manufacturer’s instructions, and the quality was evaluated by a Bioanalyser 2100 system (Agilent Technologies). Then, single-stranded DNA molecules were clustered and sequenced for 150 cycles on an Illumina HiSeq 2500 sequencing system (Illumina). Finally, paired-end reads were acquired from the Illumina HiSeq sequencer, and quality control was performed with Q30. Removal of 3′ adaptor‐trimming and low‐quality reads was conducted on Cutadapt software (v1.9.3). Alignment of the high-quality clean reads with the human reference genome (UCSC hg19) was performed on hisat2 (v2.0.4) software (http://ccb.jhu.edu/software/hisat2/index.shtml).
Additionally, cuffdiff software (part of cufflinks) was used to determine the FPKM (fragments per kilobase of exon per million fragments mapped) as the expression profiles of mRNAs and lncRNAs with the guidance of the gtf gene annotation file.
NEB Next Multiplex Small RNA Library Prep Set for Illumina (New England Biolabs) was performed to construct miRNA sequencing libraries. The quality of the libraries was evaluated by a Bioanalyser 2100 system (Agilent Technologies). Then, the single-stranded DNA molecules from denatured libraries were captured on Illumina flow cells and amplified in situ as clusters. A HiSeq 4000 sequencing system (Illumina) was used for sequencing for 50 cycles according to the manufacturer’s instructions. Cutadapt software (v1.9.3) was used to trim the adaptor sequences from the sequencing data and retain the adaptor-trimmed reads (> = 15 nt). The trimmed reads were merged and used to predict novel miRNAs with miRDeep2 software (v184.108.40.206) . The alignment of the trimmed reads was used both miRBase (http://www.mirbase.org) to obtain known pre-miRNAs and Novoalign software (v3.02.12) (http://www.novocraft.com/main/index.php) for newly predicted pre-miRNAs. The primal expression levels of the miRNA were determined via the numbers of tags on each mature miRNA. The TPM (tag counts per million aligned miRNAs) method was used for the standardization of read counts.
Analysis of DE lncRNA cis-regulatory target mRNAs
The LncRNA Disease database (http://www.cuilab.cn/lncrnadisease)  was used to identify overlapping mRNAs of non-intergenic DE lncRNAs, which including bidirectional, exon sense-overlapping, intron sense-overlapping, intronic antisense and natural antisense category DE lncRNAs. Cis-regulatory target mRNAs that were involved in HIV infection were determined by the HIV Interaction Database (https://www.ncbi.nlm.nih.gov/genome/viruses/retroviruses/hiv-1/interactions/),  LncRNA Disease database and Genecard database (https://www.genecards.org/). Moreover, the neighbouring DE mRNAs were screened within a genomic distance of 300 kb upstream or downstream of each DE lncRNA transcription start and stop site.
Construction of ceRNA networks
LncLocator  and iLoc-LncRNA  were used to screen cytoplasmic lncRNAs because only cytoplasmic lncRNAs can function as ceRNAs. Therefore, lncRNAs that were predicted to be located in the nucleus by these two software programs were excluded from analysis. In addition, DE miRNAs that had records in miRbase were screened to construct the ceRNA networks. DE miRNAs that were potential targets of DE lncRNAs were identified in DIANA tools-LncBasev.2 . Notably, the target miRNAs of some DE lncRNAs that were not recorded in DIANA tools-LncBasev.2 were predicted in miRDB databases . Subsequently, DE mRNAs that were potential targets of DE miRNAs were predicted using TargetScan  and miRDB databases. CeRNA networks were constructed according to the above analysis and visualized using Cytoscape software (v3.6.0) . Furthermore, some experimentally validated miRNA-mRNAs were queried in the miRTarBase database .
Gene Ontology (GO) annotations and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment
GO and KEGG pathway enrichment analyses were performed on the target mRNAs involved in the ceRNA networks using the Database for Annotation, Visualization and Integrated Discovery (DAVID, Version 6.8 Beta) (http://david.abcc.ncifcrf.gov/)  online functional annotation tool. GO analysis included three categories: biological process (BP), molecular function (MF), and cellular component (CC). The p-value < 0.05 indicated statistically significant enrichment.
Real‑time quantitative PCR validation
Total RNA was isolated with the miRNeasy Mini Kit (Qiagen). The RNAs was reverse transcribed with the PrimpScript ®RT reagent Kit (TAKARA) for validation of lncRNAs and mRNAs and with the Mir-X™ miRNA First-Strand Synthesis Kit (TAKARA) for validation of miRNAs. TB Green® Premix Ex Taq™ II (TAKARA) was used for RT-qPCR. GAPDH and U6 were employed as endogenous controls for lncRNA/mRNA RT-qPCR and miRNA RT-qPCR, respectively. The relative expression of lncRNAs, mRNAs and miRNAs was calculated by the 2−△△Ct method.
A nonparametric Mann–Whitney U test was performed to investigate differentially expressed genes. Data were analysed, and graphs were created using GraphPad Prism v8.0. Two-tailed p-value of less than 0.05 was considered statistically significant.
Identification of differentially expressed lncRNAs, mRNAs and miRNAs between HAART‑naive EHI patients and HCs
A total of 17,235 lncRNAs were detected (Additional file 2: Table S2). Overall, 2619 lncRNAs were only identified in the HC group, 1379 lncRNAs were only identified in the EHI patients, and 13,237 lncRNAs were identified in both groups (Fig. 1a). Thirty-six upregulated lncRNAs and 206 downregulated lncRNAs were differentially expressed (p < 0.05, fold change ≥ 2) (Fig. 1b and c). For mRNAs, a total of 16,410 mRNAs were detected (Additional file 3: Table S3). 752 mRNAs were identified only in the HC group, 406 mRNAs were identified only in the EHI patients, and 15,252 mRNAs were identified in both groups (Fig. 1d). A total of 344 upregulated mRNAs and 896 downregulated mRNAs were differentially expressed (p < 0.05, fold change ≥ 2) (Fig. 1e and f). Finally, a total of 1090 mature miRNAs were obtained via miRNA-Seq (Additional file 4: Table S4). Among these miRNAs, 269 miRNAs were identified only in the HC group, 170 miRNAs were identified only in the EHI patients, and 651 miRNAs were detected in both groups (Fig. 1g). Twenty-one DE miRNAs, including 7 significantly upregulated miRNAs and 14 significantly downregulated miRNAs (p < 0.05, fold change ≥ 1.5), were identified from 691 known miRNAs recorded in miRbase (Fig. 1h and i).
Features of DE lncRNAs expression profiling
We further characterized the expression profiling of 242 DE lncRNAs. First, DE lncRNAs were categorized into six groups based on the chromosomal position relationship of the associated coding genes, namely, bidirectional, exon sense-overlapping, intergenic, intron sense-overlapping, intronic antisense, and natural antisense, with constituent ratios of 8.68% (21/242), 6.61% (16/242), 73.55% (178/242), 2.48% (6/242), 4.13% (10/242) and 4.55% (11/242), respectively (Fig. 2a). Second, the length range of lncRNAs was from 131 nt to over 6000 nt. The majority of the DE lncRNAs were distributed in 3 intervals: 501–1000, 1001–2000 and 2001–3000 (Fig. 2b). Third, DE lncRNAs were located on 22 autosomes and the X chromosome (Fig. 2c). A higher number of DE lncRNAs were located on chromosomes 1, 8, 11, and 16. Furthermore, expression levels of most DE lncRNAs were low.
Exploration of the potential roles of lncRNAs through cis-regulatory mechanisms
LncRNAs regulate the expression of nearby protein-coding genes in a cis-regulatory fashion. First, we investigated non-intergenic lncRNAs that overlapped with mRNAs. Forty-four lncRNAs were predicted to overlap with 41 cis-regulatory target mRNAs (Additional file 5: Table S5). Among these target mRNAs, 6 mRNAs were differentially expressed and might be targets of 6 DE lncRNAs. These 6 lncRNA-mRNA pairs were all consistently regulated, including 3 pairs up and 3 pairs down (Additional file 5: Table S5). Furthermore, 4 DE mRNAs that were possibly involved in HIV infection exhibited the same regulatory direction as 4 lncRNAs. Specifically, FCAR and HIST1H2BJ overlapped with 2 exon sense-overlapping lncRNAs, and CYB5R3 and AP3M2 overlapped with 2 bidirectional lncRNAs (Table 1).
Second, we predicted the potential target mRNAs nearby DE lncRNAs. We found 21 lncRNA-mRNA pairs from 36 upregulated lncRNAs and 86 lncRNA-mRNA pairs from 206 downregulated lncRNAs (Additional file 5: Table S5). Among these lncRNA-mRNA pairs, 18 pairs were inversely regulated. Three upregulated intergenic lncRNAs (ENST00000364880, ENST00000571722 and ENST00000577988) were predicted to target the downregulated mRNA GRAP (Table 1). Four downregulated lncRNAs (ENST00000456129, ENST00000566814, ENST00000423278 and ENST00000430080) were predicted to target the downregulated mRNA SEZ6L (Table 1). In addition, the downregulated intergenic lncRNAs ENST00000492960 and TCONS_00006930 were predicted to target the downregulated mRNAs ZAP70 and BTLA, respectively, which were associated with T-cell activation and HIV-1 replication in HIV infection (Table 1).
Exploration of the potential roles of lncRNAs through ceRNA networks
To further elucidate the potential interactions of DE lncRNAs, miRNAs and mRNAs involved in EHI patients, lncRNA-miRNA-mRNA networks were constructed. A total of 160 DE lncRNAs, 21 DE miRNAs and 175 DE mRNAs were selected to construct 2 ceRNA networks: one contained 136 downregulated lncRNAs, 7 upregulated miRNAs and 97 downregulated mRNAs, and the other contained 24 upregulated lncRNAs, 14 downregulated miRNAs and 78 upregulated mRNAs (Fig. 3). We found that most lncRNAs shared miRNA response elements (MREs) of multiple miRNAs, and most miRNAs regulated multiple mRNAs. For example, 88 downregulated lncRNAs might function as ceRNAs of the upregulated hsa-miR-484 targeting 14 downregulated mRNAs (Fig. 3a). Seventeen upregulated lncRNAs might share MREs of the downregulated hsa-miR-548ah-3p with 37 upregulated mRNAs (Fig. 3b). Three lncRNAs, ENST00000504409, TCONS_00006930 and NR_024368, were predicted to have MRE binding sites for all 7 upregulated miRNAs (Fig. 3a). Furthermore, lncRNA ENST00000564287 was predicted to share MREs of 13 downregulated miRNAs (Fig. 3b). Moreover, 31 experimentally validated miRNA-mRNA interactions were found in the ceRNA networks and might be related to 134 DE lncRNAs (Additional file 6: Table S6).
GO and KEGG enrichment analysis of lncRNA-associated mRNAs in ceRNA networks
To explore the potential roles of lncRNAs in lncRNA-associated ceRNA networks, 175 DE mRNAs involved in the ceRNA networks were selected for GO and KEGG enrichment analysis (Additional file 7: Table S7). The top 5 GO terms were transcription factor activity, RNA polymerase II core promoter proximal region sequence-specific binding (GO: 0000982), transcription factor activity, sequence-specific DNA binding (GO: 0003700), sequence-specific DNA binding (GO: 0043565), transcription factor binding (GO: 0008134), and transcription factor complex (GO: 0005667) in ascending order of p-values (Fig. 4a). Five upregulated mRNAs, HIF1A, TCF7L2, FOS, FOSB and JUN occurred in several GO terms. Seventeen DE lncRNAs were predicted to share the MRE with hsa-miR-548ah-3p and in turn regulate expression of HIF1A and TCF7L2, which were reported to impact HIV-1 replication (Fig. 5).
In addition, 5 KEGG pathways were identified among the DE mRNAs involved in the lncRNA-associated ceRNA networks (Fig. 4b). We observed that 5 genes, JUN, PLCG1, ICOS, FOS and CD28, were enriched in the T cell receptor signalling pathway (hsa04660). Among these mRNAs, the proteins encoded by the FOS, JUN and FOSB genes are the main subunits of activating protein 1 (AP-1), a vital transcription factor in HIV replication and the immune response. Four DE miRNAs, miR-101-3p, miR-139-5p, miR-548ah-3p and miR-27b-3p, were found to target the above 3 mRNAs. In addition, a total of 20 DE lncRNAs were predicted to share MREs with FOS, JUN and FOSB and might regulate AP-1 as ceRNAs (Fig. 5). Note that hsa-miR-101-3p/FOS and hsa-miR-139-5p/FOS/JUN in ceRNA networks were experimentally validated (Additional file 6: Table S6). Therefore, these results provide reliable evidence for further exploring the functions of lncRNAs in EHI.
Real‑time quantitative PCR validation
The lncRNAs, miRNAs and mRNAs mentioned above, which were predicted to play important roles during early HIV infection, were validated by RT-qPCR. A total of 15 DE RNA transcripts were validated, including 6 mRNAs, 5 lncRNAs, and 4 miRNAs. The results of 12 (80%, 12/15) RNA transcriptions were consistent with the sequencing data (p < 0.05) (Fig. 6), except for lncRNA ENST00000364880, ENST00000492960 and miRNA miR-101-3p (data not shown).
In this study, we analysed lncRNA, miRNA and mRNA profiles and explored the potential functions of lncRNAs during EHI. A total of 1240 mRNAs, 242 lncRNAs and 21 known miRNAs were identified as differentially expressed during EHI. We found that several lncRNAs might regulate ZAP70, BTLA and GRAP in a cis-regulatory fashion. Some lncRNAs were predicted to act as ceRNAs to regulate expression of HIF1A, TCF7L2, FOS, FOSB and JUN. Some target mRNAs here found regulated directly or indirectly by lncRNAs have been previously reported to impact HIV replication and immune activation. Therefore, our findings may provide new ideas for further study to explore the roles of lncRNAs during EHI and help to understand the interactions between the host and HIV virus.
LncRNAs are considered to be lower expressed compared to protein coding mRNAs, and the expression of lncRNAs is spatiotemporal, and tissue-and cell-specific [33,34,35,36]. Dysregulated lncRNAs have been found to play important roles in various diseases [37, 38], including HIV infection. Several lncRNAs have been found to impact HIV-1 viral replication, the establishment of latency and the HIV-associated immune response . For instance, lncRNA MALAT1 (metastasis-associated lung adenocarcinoma transcript 1) , uc002yug.2  and HEAL (HIV-1-enhanced lncRNA)  were reported to activate HIV transcription, while NRON (noncoding repressor of nuclear factor of activated T cells) [43, 44], NEAT1 (nuclear paraspeckle assembly transcript 1) [45, 46] and AK130181  were found to repress HIV transcription. In addition, LINC00173 was shown to regulate the cytokine secretion of T cells and further affect HIV-associated immune responses . However, the functions of lncRNAs during EHI have not been reported thus far. In this study, we predicted that some lncRNAs might play important roles in HIV-1 replication and immune activation during EHI.
In this study, we speculated that some lncRNAs might regulate genes in a cis-regulatory fashion during EHI. For example, lncRNA TCONS_00006930, which was identified as a cell cycle‑associated lncRNA in endometriosis  was downregulated in EHI patients and may regulate the expression of BTLA. BTLA and ZAP70 have been found to regulate T cell activation during HIV infection [50,51,52]. In addition, ZAP-70 is required for efficient cell to cell HIV transmission and for the formation of virological synapses that promote HIV replication . In our study, we predicted that lncRNA ENST00000492960 might regulate the expression of ZAP-70 and in turn control HIV-1 replication and impact T-cell activation. Furthermore, 3 lncRNAs were predicted to regulate the expression of GRAP. GRAP plays a negative regulatory role in lymphocyte proliferation induced by T-cell receptor, interleukin-2 production and c-fos induction . These 3 lncRNAs might belong to the small nucleolar RNA (snoRNA)-ended lncRNAs (sno-lncRNAs) family. It has been reported that sno-lncRNAs may play important regulatory roles in RNA splicing and ribosomal RNA (rRNA) biogenesis [55,56,57]. However, the potential regulatory relationships between the 3 lncRNAs and GRAP need further investigation.
LncRNAs can act as ceRNAs or “miRNA sponges”, for example, lncRNA GAS5 acts on miR-873, resulting in the suppression of HIV replication . In our study, we found that 20 lncRNAs might regulate the expression of FOS, JUN, FOSB, HIF1A and TCF7L2 mRNAs as ceRNAs and then affect HIV replication and immune activation during EHI. FOS, JUN and FOSB code the common protein components (c-Fos, c-Jun and FosB) of activating protein 1 (AP-1) and they are strong activators of AP-1 . It has been reported that HIV gp120 induces endogenous c-fos and c-jun expression to enhance the function of AP-1 . AP-1 performs vital functions in T-cell activation, Th differentiation, T-cell anergy and exhaustion [61,62,63], affecting the establishment of HIV viral latency . In addition, several studies found that c-FOS was conducive to HIV-1 replication , while FosB reduced the promoter activity of the HIV long terminal repeat (LTR) . Interestingly, HIF1A also can promotes HIV replication in CD4+ T cells  and the HIV-1 accessory protein Vpr was shown to activate the LTR in an HIF1A-dependent manner in human microglial cells [68, 69]. However, TCF7L2 (TCF-4) can repress HIV replication in multiple cell types [70,71,72,73]. In astrocytes, TCF-4 was reported to interact with β-catenin and SMAR1 and then bind to LTR to decrease basal HIV transcription , and TCF-4 also decreases Sp-1-mediated transcription of the HIV promotor . Taken together, our findings revealed that lncRNAs may act as ceRNAs to impact HIV-1 replication and the host viral immune response during EHI.
In our study, 19 lncRNA/miR-101-3p/miR-139-5p/miR-548ah-3p/FOS, 6 lncRNA/miR-139-5p/JUN, 8 lncRNA/miR-27b-3p/FOSB and 17 lncRNA/miR-548ah-3p/HIF1A/ TCF7L2 were listed in the ceRNA networks. Among these 20 DE lncRNAs, lncRNA ENST00000602461 and ENST00000602507 were different transcriptions of lncRNA MIR222HG (Lnc-Ang362) which is the host transcript for miR-221 and miR-222. LncRNA MIR222HG were reported to mediate the proliferation of vascular smooth muscle cells  and facilitate the development of castration-resistant prostate cancer . In addition, lncRNA ENST00000602478, ENST00000488545 and ENST00000587762 were differentially expressed in coxsackievirus A16 (CVA16) infection, human gastric cancer and astrocytoma, respectively [76,77,78]. These lncRNAs may be new targets for the diagnosis and treatment. However, the functions of these 20 lncRNAs in HIV-1 infection were not reported previously and the mechanisms need to be further studied.
However, our study has several limitations. First, RNA-seq and miRNA-seq were performed on total PBMCs that contain multiple cell types; thus, these results represent the most prominent or common biological processes during EHI. Further studies are needed to identify the roles of lncRNAs in specialized cell types. Second, we predicted the potential roles of lncRNAs during EHI by bioinformatics analysis. Although some experimentally validated gene regulatory networks appeared in our analysis, it is still necessary to validate these results with related experiments. Third, the sample size is small because RNA samples from HAART-naive EHI patients are difficult to collect. A larger sample size is needed to verify our results.
For the first time, our study identified a lncRNA profile and predicted several lncRNAs that might play a role in HIV-1 replication and immune activation during EHI. These novel findings contribute to our understanding of the pathogenesis of HIV infection during EHI and provide new insights into antiviral therapy and vaccine development.
Availability of data and materials
All relevant data and materials during this study are included in this published article and its supplementary information files.
Human immunodeficiency virus 1
Acquired immunodeficiency syndrome
Early HIV infection
Long noncoding RNA
Competing endogenous RNA
Peripheral blood mononuclear cell
Kyoto Encyclopedia of Genes and Genomes
MiR response element
B and T lymphocyte attenuator
Zeta chain of T cell receptor-associated protein kinase 70
GRB2-related adaptor protein
Transcription factor 7-like 2
Activating protein 1
Real-time quantitative PCR
Alter G, Teigen N, Ahern R, Streeck H, Meier A, Rosenberg ES, et al. Evolution of innate and adaptive effector cell functions during acute HIV-1 infection. J Infect Dis. 2007;195(10):1452–60.
Cohen MS, Shaw GM, McMichael AJ, Haynes BF. Acute HIV-1 Infection. N Engl J Med. 2011;364(20):1943–54.
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(1):11–23.
Altfeld M, Gale M Jr. Innate immunity against HIV-1 infection. Nat Immunol. 2015;16(6):554–62.
Salido J, Ruiz MJ, Trifone C, Figueroa MI, Caruso MP, Gherardi MM, et al. Phenotype, polyfunctionality, and antiviral activity of in vitro stimulated CD8(+) T-cells from HIV(+) subjects who initiated cART at different time-points after acute infection. Front Immunol. 2018;9:2443.
Wood E, Kerr T, Rowell G, Montaner JS, Phillips P, Korthuis PT, et al. Does this adult patient have early HIV infection?: The rational clinical examination systematic review. JAMA. 2014;312(3):278–85.
Derrien T, Johnson R, Bussotti G, Tanzer A, Djebali S, Tilgner H, et al. The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 2012;22(9):1775–89.
Zhang XD, Huang GW, Xie YH, He JZ, Guo JC, Xu XE, et al. The interaction of lncRNA EZR-AS1 with SMYD3 maintains overexpression of EZR in ESCC cells. Nucleic Acids Res. 2018;46(4):1793–809.
Quinn JJ, Chang HY. Unique features of long non-coding RNA biogenesis and function. Nat Rev Genet. 2016;17(1):47–62.
Wang KC, Chang HY. Molecular mechanisms of long noncoding RNAs. Mol Cell. 2011;43(6):904–14.
Guenzl PM, Barlow DP. Macro lncRNAs: a new layer of cis-regulatory information in the mammalian genome. RNA Biol. 2012;9(6):731–41.
Marchese FP, Raimondi I, Huarte M. The multidimensional mechanisms of long noncoding RNA function. Genome Biol. 2017;18(1):206.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505(7483):344–52.
Fatica A, Bozzoni I. Long non-coding RNAs: new players in cell differentiation and development. Nat Rev Genet. 2014;15(1):7–21.
Zhang G, Li S, Lu J, Ge Y, Wang Q, Ma G, et al. LncRNA MT1JP functions as a ceRNA in regulating FBXW7 through competitively binding to miR-92a-3p in gastric cancer. Mol Cancer. 2018;17(1):87.
Biswas S, Haleyurgirisetty M, Ragupathy V, Wang X, Lee S, Hewlett I, et al. Differentially expressed host long intergenic noncoding RNA and mRNA in HIV-1 and HIV-2 infection. Sci Rep. 2018;8(1):2546.
Trypsteen W, Mohammadi P, Van Hecke C, Mestdagh P, Lefever S, Saeys Y, et al. Differential expression of lncRNAs during the HIV replication cycle: an underestimated layer in the HIV-host interplay. Sci Rep. 2016;6(1):36111.
Nair M, Sagar V, Pilakka-Kanthikeel S. Gene-expression reversal of lncRNAs and associated mRNAs expression in active vs latent HIV infection. Sci Rep. 2016;6:34862.
Trypsteen W, White CH, Mukim A, Spina CA, De Spiegelaere W, Lefever S, et al. Long non-coding RNAs and latent HIV - A search for novel targets for latency reversal. PLoS ONE. 2019;14(11):e0224879.
Chen C, Lu X, Wu N. RNA sequencing of CD4 T-cells reveals the relationships between lncRNA-mRNA co-expression in elite controller vs HIV-positive infected patients. PeerJ. 2020;8:e8911.
Fiebig EW, Wright DJ, Rawal BD, Garrett PE, Schumacher RT, Peddada L, et al. Dynamics of HIV viremia and antibody seroconversion in plasma donors: implications for diagnosis and staging of primary HIV infection. AIDS. 2003;17(13):1871–9.
Friedlander MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40(1):37–52.
Chen G, Wang Z, Wang D, Qiu C, Liu M, Chen X, et al. LncRNADisease: a database for long-non-coding RNA-associated diseases. Nucleic Acids Res. 2013;41:D983-6.
Ako-Adjei D, Fu W, Wallin C, Katz KS, Song G, Darji D, et al. HIV-1, human interaction database: current status and new features. Nucleic Acids Res. 2015;43:D566-70.
Cao Z, Pan X, Yang Y, Huang Y, Shen HB. The lncLocator: a subcellular localization predictor for long non-coding RNAs based on a stacked ensemble classifier. Bioinformatics. 2018;34(13):2185–94.
Su ZD, Huang Y, Zhang ZY, Zhao YW, Wang D, Chen W, et al. iLoc-lncRNA: predict the subcellular location of lncRNAs by incorporating octamer composition into general PseKNC. Bioinformatics. 2018;34(24):4196–204.
Paraskevopoulou MD, Vlachos IS, Karagkouni D, Georgakilas G, Kanellos I, Vergoulis T, et al. DIANA-LncBase v2: indexing microRNA targets on non-coding transcripts. Nucleic Acids Res. 2016;44(D1):D231–8.
Chen Y, Wang X. miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020;48(D1):D127–31.
Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4:e05005.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Chou CH, Shrestha S, Yang CD, Chang NW, Lin YL, Liao KW, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2018;46(D1):D296–302.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Cabili MN, Trapnell C, Goff L, Koziol M, Tazon-Vega B, Regev A, et al. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 2011;25(18):1915–27.
Li F, Xiao Y, Huang F, Deng W, Zhao H, Shi X, et al. Spatiotemporal-specific lncRNAs in the brain, colon, liver and lung of macaque during development. Mol Biosyst. 2015;11(12):3253–63.
Fan Z, Zhao M, Joshi PD, Li P, Zhang Y, Guo W, et al. A class of circadian long non-coding RNAs mark enhancers modulating long-range circadian gene regulation. Nucleic Acids Res. 2017;45(10):5720–38.
Gloss BS, Dinger ME. The specificity of long noncoding RNA expression. Biochim Biophys Acta. 2016;1859(1):16–22.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs insights into functions. Nat Rev Genet. 2009;10(3):155–9.
Zhao Z, Zhang M, Ying J, Hu X, Zhang J, Zhou Y, et al. Significance of genetic polymorphisms in long non-coding RNA AC079767.4 in tuberculosis susceptibility and clinical phenotype in Western Chinese Han population. Sci Rep. 2017;7(1):965.
Shen L, Wu C, Zhang J, Xu H, Liu X, Wu X, et al. Roles and potential applications of lncRNAs in HIV infection. Int J Infect Dis. 2020;92:97–104.
Qu D, Sun W-W, Li L, Ma L, Sun L, Jin X, et al. Long noncoding RNA MALAT1 releases epigenetic silencing of HIV-1 replication by displacing the polycomb repressive complex 2 from binding to the LTR promoter. Nucleic Acids Res. 2019;47(6):3013–27.
Huan C, Li Z, Ning S, Wang H, Yu X, Zhang W. Long noncoding RNA uc002yug.2 activates HIV-1 latency through regulation of mRNA levels of various RUNX1 isoforms and increased Tat expression. J Virol. 2018;92(9):e01844-17.
Chao TC, Zhang Q, Li Z, Tiwari SK, Qin Y, Yau E, et al. The long noncoding RNA HEAL regulates HIV-1 replication through epigenetic regulation of the HIV-1 promoter. mBio. 2019;10(5):e02016.
Imam H, Bano AS, Patel P, Holla P, Jameel S. The lncRNA NRON modulates HIV-1 replication in a NFAT-dependent manner and is differentially regulated by early and late viral proteins. Sci Rep. 2015;5:8639.
Li J, Chen C, Ma X, Geng G, Liu B, Zhang Y, et al. Long noncoding RNA NRON contributes to HIV-1 latency by specifically inducing tat protein degradation. Nat Commun. 2016;7:11730.
Jin C, Peng X, Xie T, Lu X, Liu F, Wu H, et al. Detection of the long noncoding RNAs nuclear-enriched autosomal transcript 1 (NEAT1) and metastasis associated lung adenocarcinoma transcript 1 in the peripheral blood of HIV-1-infected patients. HIV Med. 2016;17(1):68–72.
Liu H, Hu PW, Couturier J, Lewis DE, Rice AP. HIV-1 replication in CD4+ T cells exploits the down-regulation of antiviral NEAT1 long non-coding RNAs following T cell activation. Virology. 2018;522:193–8.
Li H, Chi X, Li R, Ouyang J, Chen Y. A novel lncRNA, AK130181, contributes to HIV-1 latency by regulating viral promoter-driven gene expression in primary CD4(+) T cells. Mol Ther Nucleic Acids. 2020;20:754–63.
Postler TS, Pantry SN, Desrosiers RC, Ghosh S. Identification and characterization of a long non-coding RNA up-regulated during HIV-1 infection. Virology. 2017;511:30–9.
Liu J, Wang Q, Zhang R, Zhang C, Lin J, Huang X. Identification of LINC01279 as a cell cycle associated long noncoding RNA in endometriosis with GBA analysis. Mol Med Rep. 2018;18(4):3850–8.
Hung CH, Thomas L, Ruby CE, Atkins KM, Morris NP, Knight ZA, et al. HIV-1 Nef assembles a Src family kinase-ZAP-70/Syk-PI3K cascade to downregulate cell-surface MHC-I. Cell Host Microbe. 2007;1(2):121–33.
Dumont C, Blanchard N, Di Bartolo V, Lezot N, Dufour E, Jauliac S, et al. TCR/CD3 down-modulation and zeta degradation are regulated by ZAP-70. J Immunol. 2002;169(4):1705–12.
Zhang Z, Xu X, Lu J, Zhang S, Gu L, Fu J, et al. B and T lymphocyte attenuator down-regulation by HIV-1 depends on type I interferon and contributes to T-cell hyperactivation. J Infect Dis. 2011;203(11):1668–78.
Sol-Foulon N, Sourisseau M, Porrot F, Thoulouze MI, Trouillet C, Nobile C, et al. ZAP-70 kinase regulates HIV cell-to-cell spread and virological synapse formation. EMBO J. 2007;26(2):516–26.
Shen R, Ouyang YB, Qu CK, Alonso A, Sperzel L, Mustelin T, et al. Grap negatively regulates T-cell receptor-elicited lymphocyte proliferation and interleukin-2 induction. Mol Cell Biol. 2002;22(10):3230–6.
Yin QF, Yang L, Zhang Y, Xiang JF, Wu YW, Carmichael GG, et al. Long Noncoding RNAs with snoRNA Ends. Mol Cell. 2012;48(2):219–30.
Xing YH, Yao RW, Zhang Y, Guo CJ, Jiang S, Xu G, et al. SLERT regulates DDX21 rings associated with Pol I. Cell. 2017;269(4):664–78.
Xing YH, Ling-Ling C. Processing and roles of snoRNA-ended longnoncoding RNAs. Crit Rev Biochem Mol Biol. 2018;53(6):596–606.
Chen L, Chen L, Zuo L, Gao Z, Shi Y, Yuan P, et al. Short communication: long noncoding RNA GAS5 inhibits HIV-1 replication through interaction with miR-873. AIDS Res Hum Retroviruses. 2018;34(6):544–9.
Hess J, Angel P, Schorpp-Kistner M. AP-1 subunits: quarrel and harmony among siblings. J Cell Sci. 2004;117(Pt 25):5965–73.
Gibellini D, Re MC, Bassini A, Guidotti L, Catani L, La Placa M, et al. HIV-1 gp120 induces the activation of both c-fos and c-jun immediate-early genes in HEL megakaryocytic cells. Br J Haematol. 1999;104(1):81–6.
Foletta VC, Segal DH, Cohen DR. Transcriptional regulation in the immune systemall roads lead to AP-1. J Leukoc Biol. 1998;63(2):139–52.
Yukawa M, Jagannathan S, Vallabh S, Kartashov AV, Chen X, Weirauch MT, et al. AP-1 activity induced by co-stimulation is required for chromatin opening during T cell activation. J Exp Med. 2020;217(1):e20182009.
Atsaves V, Leventaki V, Rassidakis GZ, Claret FX. AP-1 transcription factors as regulators of immune responses in cancer. Cancers. 2019;11(7):1037.
Duverger A, Wolschendorf F, Zhang M, Wagner F, Hatcher B, Jones J, et al. An AP-1 binding site in the enhancer/core element of the HIV-1 promoter controls the ability of HIV-1 to establish latent infection. J Virol. 2013;87(4):2264–77.
Roebuck KA, Gu DS, Kagnoff MF. Activating protein-1 cooperates with phorbol ester activation signals to increase HIV-1 expression. AIDS. 1996;10(8):819–26.
Rizzi C, Crippa MP, Jeeninga RE, Berkhout B, Blasi F, Poli G, et al. Pertussis toxin B-oligomer suppresses IL-6 induced HIV-1 and chemokine expression in chronically infected U1 cells via inhibition of activator protein 1. J Immunol. 2006;176(2):999–1006.
Duette G, Pereyra Gerber P, Rubione J, Perez PS, Landay AL, Crowe SM, et al. Induction of HIF-1α by HIV-1 infection in CD4 + T cells promotes viral replication and drives extracellular vesicle-mediated inflammation. mBio. 2018;9(5):e00757-18.
Deshmane SL, Mukerjee R, Fan S, Del Valle L, Michiels C, Sweet T, et al. Activation of the oxidative stress pathway by HIV-1 Vpr leads to induction of hypoxia-inducible factor 1alpha expression. J Biol Chem. 2009;284(17):11364–73.
Deshmane SL, Amini S, Sen S, Khalili K, Sawaya BE. Regulation of the HIV-1 promoter by HIF-1alpha and Vpr proteins. Virol J. 2011;8:477.
Kumar A, Zloza A, Moon RT, Watts J, Tenorio AR, Al-Harthi L. Active beta-catenin signaling is an inhibitory pathway for human immunodeficiency virus replication in peripheral blood mononuclear cells. J Virol. 2008;82(6):2813–20.
Narasipura SD, Henderson LJ, Fu SW, Chen L, Kashanchi F, Al-Harthi L. Role of beta-catenin and TCF/LEF family members in transcriptional activity of HIV in astrocytes. J Virol. 2012;86(4):1911–21.
Henderson LJ, Narasipura SD, Adarichev V, Kashanchi F, Al-Harthi L. Identification of novel T cell factor 4 (TCF-4) binding sites on the HIV long terminal repeat which associate with TCF-4, beta-catenin, and SMAR1 to repress HIV transcription. J Virol. 2012;86(17):9495–503.
Rossi A, Mukerjee R, Ferrante P, Khalili K, Amini S, Sawaya BE. Human immunodeficiency virus type 1 Tat prevents dephosphorylation of Sp1 by TCF-4 in astrocytes. J Gen Virol. 2006;87(Pt 6):1613–23.
Leung A, Trac C, Jin W, Lanting L, Akbany A, Saetrom P, et al. Novel long noncoding RNAs are regulated by angiotensin II in vascular smooth muscle cells. Circ Res. 2013;113(3):266–78.
Sun T, Du SY, Armenia J, Qu F, Fan J, Wang X, et al. Expression of lncRNA MIR222HG co-transcribed from the miR-221/222 gene promoter facilitates the development of castration-resistant prostate cancer. Oncogenesis. 2018;7(3):30.
Shi Y, Tu H, Chen X, Zhang Y, Chen L, Liu Z, et al. The long non-coding RNA expression profile of Coxsackievirus A16 infected RD cells identified by RNA-seq. Virol Sin. 2016;31(2):131–41.
Ding J, Li D, Gong M, Wang J, Huang X, Wu T, et al. Expression and clinical significance of the long non-coding RNA PVT1 in human gastric cancer. Onco Targets Ther. 2014;7:1625–30.
Zhi F, Wang Q, Xue L, Shao N, Wang R, Deng D, et al. The use of three long non-coding RNAs as potential prognostic indicators of astrocytoma. PLoS ONE. 2015;10(8):e0135242.
We thank all the participants in this study and the native English-speaking medical editors from the Nature Research Editing Service (http://authorservices.springerna ture.com/language-editing/).
This work was supported by the Mega-projects of National Science Research for the 13th Five-Year Plan (2017ZX10201101), National Natural Science Foundation of China (81871637) and Chinese Academy of Medical Sciences Research Unit (No. 2019RU017), China Medical University.
Ethics approval and consent to participate
This study was reviewed and approved by the Medical Research Ethics Committee of the First Affiliated Hospital of China Medical University. All participants provided written informed consent prior to research participation.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Clinical characteristics of participants
LncRNA expression profiling.
mRNA expression profiling.
miRNA expression profiling.
Predictive cis-target mRNAs of differentially expressed lncRNAs.
Experimentally validated miRNA-mRNA interactions from miRTarBase in ceRNA network.
Enrichment analyses of target mRNAs in ceRNA neworks.
About this article
Cite this article
Ma, L., Zhang, H., Zhang, Y. et al. Integrated analysis of lncRNA, miRNA and mRNA profiles reveals potential lncRNA functions during early HIV infection. J Transl Med 19, 135 (2021). https://doi.org/10.1186/s12967-021-02802-9