Role of Epstein-Barr Virus and Human Papillomavirus Coinfection in Cervical Intraepithelial Neoplasia in Chinese Women Living With HIV

Given that only a small percentage of human papillomavirus (HPV)-positive women develop cancer, HPV is necessary but insufficient for carcinogenesis. Mucosally transmitted viral cofactors appear to contribute to HPV-related cervical cancer, such as Epstein-Barr virus (EBV), but previous studies have shown inconsistent outcomes. The exact role of EBV in cervical cancer remains unclear, and more studies are needed to determine its involvement. In this study, we describe the prevalence of EBV and HPV coinfection in HIV-positive women and explore how abnormal host immune status induced by viral coinfections modulates epithelial gene expression. We found a significant correlation between EBV-HPV coinfection and the incidence of high-grade cervical intraepithelial neoplasia (CIN2+). RNA sequencing indicated that CIN tissues coinfected with EBV and HPV led to significant changes in the gene expression of epithelial differentiation and development compared to normal tissues with HPV infection alone. In particular, several differentially expressed genes (DEGs) are closely associated with cancer, such as CACNG4, which was confirmed to be upregulated at both the mRNA and protein levels. Therefore, these findings provide some evidence that EBV may act as a cofactor or mediator in HPV-related cervical cancer. Specific genes or proteins, such as CACNG4, may serve as biomarkers that can risk stratify patients based on pathological changes in the cervix.

Introduction

Sexually transmitted infections (STIs) can cause cancer or have serious reproductive health consequences (WHO, 2018), especially in women. Human papillomavirus (HPV) infection is the most prevalent STI in the female population (Satterwhite et al., 2013), with an estimated prevalence of 11-12% worldwide (Forman et al., 2012). Cervical cancer is the main HPV-related cancer, and according to the WHO, there were an estimated 570,000 new cases worldwide in 2018 (WHO). Although the etiological role of HPV in cervical cancer is well established, more than 90% of HPV infections can be cleared within a few months to 2 years without any intervention (Moscicki et al., 2006). Only 10% of women with high-risk HPV will develop long-lasting infections that put them at risk for cancer (Centers for Disease Control and Prevention, 2020), indicating that HPV is necessary but may not be sufficient for carcinogenesis. In addition to HPV infection, several biological and environmental cofactors have been implicated in the development of cervical cancer, such as immune status and coinfection with other sexually transmissible pathogens (Castellsague et al., 2002; Bosch and de Sanjose, 2007). Detection of viral cofactors involved in malignancy progression can be used for HPV-related cervical cancer intervention.

Similar to HPV, Epstein-Barr virus (EBV) is a mucosally transmitted pathogen that is widespread in the human population (Higgins et al., 2007). This virus is able to persist in the human host through the establishment of a lifelong latent infection, which could contribute to host immune modulation (Djaoud et al., 2017). EBV can infect cervical epithelial cells (Sixbey et al., 1983), and viral DNA, RNA or proteins of EBV have been detected in exfoliated cells or tissues from cervical cancer or intraepithelial neoplasia (CIN) lesions (Khenchouche et al., 2013; Marinho-Dias et al., 2013). A correlation between EBV infection and abnormal cervical cytology has been shown; on one hand, the prevalence of EBV positivity increased with lesion severity (Silver et al., 2011; Aromseree et al., 2015); on the other hand, CIN or cervical carcinoma occurs more often among EBV positive women than women without EBV infection (Cameron et al., 2018; de Lima et al., 2018). Recently, one study (Cameron et al., 2020) from a cohort of HIV-seropositive women, reported an approximately fourfold increased risk of abnormal cervical cytology in women with cervical EBV and high-risk HPV compared to women with high-risk HPV alone, suggesting that cervical EBV appears to predict a greater risk of cervical dysplasia in HIV-infected women with a high-risk HPV infection. Although a few studies (Kahla et al., 2012; Blanco et al., 2020) have provided explanations for how EBV can act as a risk factor for HPV persistence and tumor development, the exact role of EBV in cervical cancer remains unclear. In addition, the opposite result has also been reported, showing a lack of association between EBV infection and HPV-related cervical cancer (Ammatuna et al., 2000). Therefore, more studies are needed to determine whether EBV has a role in cervical cancer or acts only as an opportunistic pathogen.

In contrast to people with a normal immune system, women living with HIV are at high risk of acquiring oncogenic HPV and are also susceptible to persistent and reactivated HPV infection and to developing cancer (Ahdieh et al., 2001; Camargo et al., 2014; Jin et al., 2017). In addition, EBV incidence is also higher in the HIV-infected population (Miller et al., 2006; Jin et al., 2017). Therefore, the increased comorbidity among these sexually transmitted viruses in HIV-positive individuals may provide a better model to explore the possible mechanism of coinfection-associated malignancies. In this context, we detected viral DNA in exfoliated cells from cervix to describe the prevalence of coinfection between EBV and HPV in HIV-positive women and to determine the association between viral coinfection and CIN using mRNA sequencing (mRNA-seq). If EBV is a risk factor for HPV-related cervical cancer, we would expect that certain alterations in genes or proteins under specific coinfection conditions can be found to be potential predictive biomarkers for tumorigenesis in cervical disease.

Materials and Methods

Study Population and Procedure

Participants were recruited from the AIDS Antiretroviral Therapy Department of the Third People’s Hospital of Kunming, Yunnan Province, China, between February and November 2019. The inclusion criteria were as follows: 1) HIV-positive women with ongoing antiretroviral therapy, 2) age of 18-64 years, 3) sexually active, 4) ability to both physically and mentally undergo cervical sampling and colposcopy examinations, and 5) intact cervix and no current pregnancy. The study protocol was approved by the Ethical Review Committees of National Cancer Center/Cancer Hospital Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China. All participants provided informed consent. A summary of the study procedure is shown in Figure 1.

FIGURE 1

Figure 1 A summary of the study procedure.

Sample Collection, HPV DNA Testing, and Cytology

Trained gynecologists inserted cytology brushes and collected exfoliated cervical cells after a pelvic examination, and cervical samples were collected into a prepared liquid medium (ThinPrep, Hologic, Marlborough, MA, USA). All collected samples were placed in a 4°C refrigerator and transferred to a third-party laboratory, KingMed Diagnostics, Kunming branch, for Cobas 4800 HPV testing (Roche, Basel, Switzerland) and Sansure PCR HPV test (Sansure Biotech, China). Sansure PCR HPV test (Zhao et al., 2019) is a PCR-based pioneered One-step Fast Release technology and has been approved by European Union Certificate. The test uses real-time fluorescent quantitative PCR to target and report 15 high-risk HPV types, including HPV-16, 18, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, and 68. Cycle threshold (Ct) values ≤ 39 is considered HPV positive and >39 is considered negative. Then the cytology slide was prepared and the cytological interpretation was made by experienced cytologists according to the Bethesda 2014 classification system (Diane Solomon, 2015).

Detection of EBV

Total DNA was extracted from exfoliated cervical cells from liquid medium (ThinPrep) using a QIAamp DNA Mini Kit (Qiagen, Germany) according to the protocol recommended by the manufacturer. To detect EBV, 5 µl of DNA was then used for subsequent RT-PCR (Premix Ex Taq, Takara, Janpa) with a CFX96 Real-Time PCR Detection System (Bio-Rad, USA). In addition, β-globin primers were used to verify the quality of the DNA and exclude false negative results, and a positive control was included in each RT-PCR assay to ensure the success of the reaction. RT-PCR was performed for 45 cycles of amplification, and samples were considered negative if the cycle threshold (Ct) value was greater than 40. The primers and probes used in this study are identical to those in previous studies (Lo et al., 1998; Wadowsky et al., 2003).

Colposcopy Examination and Pathology

Women positive for HPV on Cobas or Sansure PCR testing or with a cytological diagnosis of atypical squamous cells of undetermined significance (ASCUS) or higher (ASCUS+) were referred for a colposcopy examination. A directed biopsy was collected if any abnormality was identified under colposcopy. Endocervical curettage was performed if the squamocolumnar junction was invisible. Biopsy tissues were immediately immersed in 10% buffered formalin and transported to KingMed Diagnostics, Kunming branch for processing (fixed colposcopy-directed biopsy samples were dehydrated in graded ethanol, embedded in paraffin and sectioned, and hematoxylin and eosin (HE) staining was conducted according to standard protocols) and diagnosis by experienced pathologists blinded to other screening results. Pathology results were reported as negative, cervical intraepithelial neoplasia grade 1 (CIN1), grade 2 (CIN2), grade 3 (CIN3), microinvasive carcinoma (MIC), squamous cell carcinoma (SCC), adenocarcinoma in situ (AIS), or adenocarcinoma (ADC). Women with negative results by all tests (cytology <ASCUS and without HPV detection) were considered to be negative for the outcome of CIN2 or worse (CIN2+). Women with confirmed CIN2+ were recommended for treatment.

Transcriptional Profiling and Analysis

Sample Preparation

Total RNA was extracted from formalin-fixed paraffin embedded (FFPE) biopsy tissues using the RNeasy FFPE Kit (Qiagen, Germany) according to the manufacturer’s instructions. RNA purity was evaluated using a NanoPhotometer spectrophotometer (Implen, USA), and integrity was assessed using an RNA Nano 6000 Assay Kit with a Bioanalyzer 2100 system (Agilent Technologies, USA). Next, the RNA concentration was measured using a Qubit RNA Assay Kit with a Qubit 2.0 Fluorometer (Life Technologies, USA). A total amount of 1 µg of RNA per sample was used to generate sequencing libraries. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations at an elevated temperature in First Strand Synthesis Reaction Buffer. Then, a random hexamer primer and mRNA fragments as templates were used to synthesize the first-strand cDNA. Second-strand cDNA fragments were obtained using DNA Polymerase I and dNTP. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. The 3’end adenylated and adaptor-ligated cDNA were then purified with AMPure XP system (Beckman Coulter, USA) to select the fragments with the correct target size (370-420 bp). Finally, the cDNA libraries were enriched by PCR amplification with Phusion High-Fidelity DNA polymerase, universal PCR primers and Index (X) Primer. The library quality was assessed on an Agilent Bioanalyzer 2100 system.

Before sequencing, clustering of the index-coded samples was performed on a cBot Cluster Generation System using a TruSeq PE Cluster Kit v3-cBot-HS (Illumina) according to the manufacturer’s instructions. Then, the library preparations were sequenced on an Illumina HiSeq platform, and 150 bp paired-end reads were generated.

Quality Control and Read Mapping

Clean reads were obtained by removing reads containing adapters, reads containing poly-N and low-quality reads from the raw data. All subsequent analyses were based on filtered reads. Paired-end filtered reads were mapped to the homo sapiens sequence (ftp://ftp.ensembl.org/pub/release-98/fasta/homo_sapiens/) with HISAT2 v2.0.5.

mRNA Expression Analysis

FeatureCounts v1.5.0-p3 was used to count the read numbers mapped to each gene. To normalize the read count, fragments per kilobase of transcript per million mapped fragments (FPKM) were calculated for each gene. Then, the DESeq R package (1.16.1) was utilized to perform differential expression analysis of the two groups. The resulting P values were adjusted using Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted P value <0.05 or | log2-fold change | >0 were considered differentially expressed. Gene Ontology (GO) enrichment analysis of differentially expressed genes (DEGs) was implemented by the GOseq R package. GO terms with corrected P values less than 0.05 were considered significantly enriched by DEGs. Subsequently, we used the clusterProfiler R package to test the statistical enrichment of DEGs in KEGG pathways. The mRNA-seq data for this study can be found in the Sequence Read Archive (SRA) (the BioProject accession number is PRJNA701124, and the SRA records will be accessible with the following link upon publication: https://www.ncbi.nlm.nih.gov/sra/PRJNA701124).

Quantitative Real Time-PCR (qRT-PCR)

To validate the results from mRNA-seq, 14 DEGs were selected for qRT-PCR test. The involved genes and primers are shown in Table 1. RNA samples were the same as the ones used for mRNA-seq. qRT-PCR was performed using One Step SYBR PrimeScript PLUS RT-PCR kit (Takara, Japan) on a CFX96 Real-Time PCR Detection System (Bio-Rad, USA). The reaction was performed as follows: 42°C for 5 min and 95°C for 10 s, followed by 40 cycles of 95°C for 5 s and 60°C for 30 s. The dissociation curve analysis were used to confirm the specificity of amplification of each product and the absence of primer dimer. Each sample was performed in triplicates. The GAPDH gene was used as a reference gene. The relative expression level of each gene was measured by the method of 2^-ΔΔCt and was standardized by reference gene.

TABLE 1

Table 1 Sequences of RT-PCR primers used for validation.

Data Analysis in Public Database

We utilized GEPIA2 to access and analyze the public database on cancer tissues generated by The Cancer Genome Atlas (TCGA) program and normal tissues from the Genotype-Tissue Expression (GTEx) project.

Immunofluorescence (IF) Assays and Western Blotting (WB)

IF was performed on FFPE sections after heat-induced antigen retrieval. Briefly, sections were blocked for 1 h, incubated overnight at 4°C with a rabbit primary antibody against human CACNG4 (Immunoway, YM3404), and subsequently incubated with a secondary Alexa Fluor 568-conjugated goat anti-rabbit IgG antibody (Abcam, ab175696) for 1 h. Slides were mounted with mounting medium containing DAPI and imaged by a Leica SP8 laser scanning confocal microscope system. For WB, samples were pooled and extracted using the FFPE Total Protein Extraction Kit (Sangon Biotech, C500058) and then incubated with primary antibody against CACNG4 or β-actin (Cell Signaling, #3700). The proteins were visualized using HRP-conjugated secondary antibody (Cell Signaling, #7040 or #7076) and chemiluminescent HRP substrate (Bio-Rad).

Statistical Analysis

Descriptive statistical analysis was performed to summarize the prevalence of HPV and EBV infections and coinfections. Chi-square tests were used to compare the distribution of CIN2+ between the HPV infection and HPV-EBV coinfection groups. Student t-test was performed to evaluate the mRNA expression difference between two groups. p values less than 0.05 (two-sided) were considered statistically significant. Data analyses were performed using SPSS 20.0.

Results

Prevalence of HPV and EBV Infections and Coinfections in HIV-Positive Women

A total of 180 women were enrolled in this study; of these, four women with extremely high HIV viral load (≥10⁶ copies/ml) that reflected severe immunosuppression were excluded to avoid interference. Of 176 women, the median age was 40 (34–46) years. Among them, 61.9% have been married, 68.8% had education at junior high school or below, 45.4% were unemployed and 44.3% were farmer or rural migrant worker. The median age at sexual debut and the first delivery were 20 (19–22) and 24 (21–27) years, respectively. The median CD4 count was 540 (IQR 398–656) cells/µL, and the median HIV viral load was <50 (IQR <50–<50) copies/ml. Valid HPV test results were obtained from 176 women. The prevalence of HPV and EBV in HIV-positive women was 33.5% and 15.9%, respectively. In addition, the coinfection rate of EBV-HPV was 9.1%. Eleven participants were missing from the histopathological analyses due to loss to follow-up. Thus, a total of 165 women were involved in the distribution analysis of CIN2+ among different coinfection statuses, and Table 2 shows detailed information on the STI distribution among these women. Table 3 shows that EBV-HPV coinfected women were more likely to be diagnosed with CIN2+ (p=0.017) than those without coinfection.

TABLE 2

Table 2 Detailed information on HPV and EBV distribution among 165 HIV-positive women.

TABLE 3

Table 3 Distribution of CIN2+ in HPV- and HIV-positive women with or without EBV.

Comparison of General mRNA Expression Profiles Between Different Infection Statuses in HIV-Positive Women

To clarify the reasons behind the significant correlation of EBV and HPV coinfection and CIN2+, we used mRNA-seq to capture the transcriptome changes in pathological tissues from colposcopy-directed biopsy samples. Owing to the relatively small patient number for biopsy, 7 samples within two groups were included in this analysis. The individual demographical and clinical characteristics are presented in Supplementary Table 1. The CIN grade and histopathology of these samples are presented in Figure 2. Total mRNA from the coinfection group was compared to the HPV alone infection group, revealing that there were 187 significantly altered genes (Figure 3A). The resulting heatmap confirmed that patients with the same infection status could be grouped together based on similar gene expression profiles (Figure 3B). Next, a GO analysis was used to identify potential biological or molecular mechanisms that were significantly altered in the coinfection group. The most significantly enriched common GO terms were associated with epidermal cell and keratinocyte differentiation, skin and epidermis development, and peptide cross-linking (Figure 3C). Most genes involved in these significantly changed GO terms were downregulated in the coinfection condition. However, several upregulated DEGs were significantly enriched in biological process GO terms concerning the regulation of signaling receptor activity, protein maturation, positive regulation of signaling receptor activity, regulation of complement activation and protein activation cascade. These findings suggest that coinfection of EBV and HPV may cause increased virus-host interactions and have a greater effect on host immunity.

FIGURE 2

Figure 2 Pathology of cervical tissues involved in mRNA-seq.

FIGURE 3

Figure 3 Analysis of cervical tissue transcriptome data of different viral infection statuses. (A) Differences in gene expression in the EBV and HPV coinfection groups compared with the HPV infection alone group (red and green dots indicate DEGs); the numbers in plots indicate the numbers of DEGs. (B) Hierarchical cluster analysis of DEGs within different paired comparisons confirming that patients with each infection or coinfection status could be grouped together based on similar gene expression profiles (red, upregulated; blue, downregulated). (C) Numbers of DEGs belonging to numerous significantly altered GO enrichment terms (up, upregulated DEGs; down, downregulated DEGs; total, total DEGs).

Specific DEGs and Overexpression of CACNG4 in EBV-HPV-Associated CIN Lesions

By combining the findings from the identified DEGs and the GO analysis, we selected 14 DEGs that were included in the enriched GO terms related to EBV and HPV coinfection or were the most up- or downregulated DEGs. RT-PCR was performed to validate the expression profiles of these marker genes in the EBV-HPV group, and the results were consistent with the mRNA-seq data (Figures 4A, B). Notably, a significant increase in CACNG4 expression was observed among these genes. CACNG4, an L-type voltage-gated calcium channel gamma subunit, has been reported to be overexpressed in several cancers (Kuznetsova et al., 2007; Zhang et al., 2016; Slattery et al., 2018) and may promote calcium homeostasis to increase the survival and metastatic ability of tumor cells (Kanwar et al., 2020). This alteration in CACNG4 expression may be linked to CIN lesions. IF (Figure 4C) and WB (Figures 4D, E) analyses were further performed and confirmed the specific upregulation of CACNG4 in EBV-HPV-coinfected CIN lesions at the protein level. Next, we examined the levels of CACNG4 mRNA transcripts within cancer samples using the GEPIA2 website. A comparison of CACNG4 expression between tumor samples and normal tissue samples taken from TCGA and the GTEx databases showed that CACNG4 is upregulated in cervical cancer (Figure 4F). Taken together, these data suggest that CACNG4 may act as a potential biomarker for the early diagnosis of cervical cancer, as well as a target for modulating intracellular calcium homeostasis, which sustains tumor growth.

FIGURE 4

Figure 4 Selected mRNAs with significant changes and overexpression of CACNG4 in EBV-HPV-associated CIN lesions. Expression levels of the selected, significantly altered mRNAs of the mRNA-seq analysis (A) and RT-PCR validation (B). Each value was normalized to the GAPDH value and is presented as the mean ± SD, *p<0.05, compared with the HPV group. (C) A red fluorescence signal indicating the presence of CACNG4 was detected in cervical tissues. (D) Protein expression of CACNG4 relative to β-actin expression in cervical tissues. (E) Densitometric analysis of the immunoblots. The intensity of each band was normalized to the β-actin intensity. Due to the limited amount of clinical material, the EBV-HPV sample consisted of 4 pooled biological replicates and the HPV sample consisted of 3 pooled biological replicates in (D, E). (F) The expression of CACNG4 between cervical cancer samples and normal tissue samples taken from TCGA and the GTEx databases. TPM, transcripts per million.

Discussion

In this study, we showed the prevalence of HPV and EBV as single infections or coinfections among women living with HIV, and our results were similar to what has been previously reported (Camargo et al., 2014; Jin et al., 2017). Our findings also support that HIV-infected individuals have an increased risk for EBV infection compared with the general population (Petrara et al., 2012; Jin et al., 2017). Given that only a small percentage of HPV-positive women develop cancer, HPV is necessary but insufficient for carcinogenesis. Thus, mucosally transmitted viral cofactors appear to contribute to HPV-related cervical cancer. Human herpesvirus EBV shares similar routes and sites of infection with HPV. Genital coinfection of these viruses may disrupt the mucosal epithelial barrier (Hebner and Laimins, 2006) and the different viruses may influence each other in many ways, such as viral entry and clearance (Blanco et al., 2020). Consistent with previous reports (Santos et al., 2009; Kienka et al., 2019; Cameron et al., 2020; Joharinia et al., 2020), we revealed a significant association between EBV and HPV coinfection and CIN2+, suggesting that coinfection of EBV and HPV is a potential risk factor in carcinogenesis.

While there is a close relationship between EBV infection and HPV-related cervical cancers, the underlying mechanisms remain unclear. Several previous studies have examined gene expression changes in HPV-infected cell lines or tissues (Wong et al., 2006; Klymenko et al., 2017). However, no studies have analyzed the features of host gene expression in the presence of other viruses coinfected with HPV. In this study, exfoliated cells from the cervix were utilized to determine the viral infection status, followed by evaluating the global changes in the biopsy tissue using mRNA-seq. Due to the limited samples, we did not detect the viruses in these biopsy tissues. However, of particular interest is how abnormal host immune status induced by viral coinfections modulates epithelium gene expression. Based on our data, we revealed that obvious changes in the transcriptome in CIN tissues with EBV and HPV coinfection compared to normal tissues with HPV infection alone. Similar to previous transcriptomic studies (Klymenko et al., 2017) showing that HPV infection alone in epithelial cell lines abrogates differentiation and epithelial barrier formation, the current study also showed decreased epithelial and keratinocyte differentiation factors in the coinfection groups. For example, small proline-rich repeat (SPRR) protein family members, including SPRR2A, SPRR2B, SPRR2D, SPRR2E, SPRR2F, and SPRR2G, the late cornified envelope (LCE) gene cluster members LCE3D and LCE3E, and a number of keratin (KRT) protein-coding genes, such as KRT17, KRT23 and KRT24, were downregulated (Supplementary Table 2). These findings suggest that coinfection of EBV and HPV increased the effect of HPV on epithelial differentiation and development, but these changes may mainly show a profile expected of viral infection rather than tumor progression. Meanwhile, the presence of virus and viral gene expression is unknown in the mRNA-profiled specimens and the histopathology grade of the HPV+ and HPV+EBV+ samples are not matched, so the findings are compelling but further studies are needed to confirm that EBV is responsible for the gene expression changes observed.

Intriguingly, several DEGs involved in EBV and HPV coinfection attracted our attention. The GO analysis indicated that some DEGs, including CACNG4, C7, C1QA, C1QB, CD55, and SCGB3A1, were involved in the immune response, signaling receptor activity, and protein activation, suggesting that coinfection with EBV and HPV may cause interactions with the host. In particular, some of these DEGs are closely associated with cancer. For example, genes related to the GO term “components of complement,” such as C1QA, C1QB, and C7, may promote tumorigenesis in cells within the tumor microenvironment (Bulla et al., 2016). TFEC is a member of the microphthalmia (MiT) family, which regulates multiple physiological processes, including cell survival, differentiation, proliferation, invasion, metabolism, and DNA damage repair, thus contributing to the initiation and development of some cancer types (Haq and Fisher, 2011; Goding and Arnheiter, 2019). We confirmed these changes in mRNA levels (Figures 4A, B); unfortunately, due to the limited amount of clinical material, we could not evaluate all of their protein expression levels. We therefore selected CACNG4 as a candidate to evaluate its protein levels because of the highest upregulation of CACNG4 mRNA levels (mean fold change=8.55). Indeed, a consistent upregulation of CACNG4 was found in cervical tissues (Figures 4D, E). Of interest, CACNG4 encodes an L-type voltage-gated calcium channel (VGCC) γ4 subunit. One gene expression analysis revealed upregulation of the CACNG4 gene in human colon cancer (Slattery et al., 2018). Breast tumors with nodal metastases are also associated with higher protein expression of CACNG4 in clinical samples (Kuznetsova et al., 2007). Mechanistically, CACNG4 affects calcium influx by modulating VGCCs, which in turn regulate the homeostasis and metastasis of tumor cells (Kanwar et al., 2020). Since VGCCs of different types have been associated with several cancers (Monteith et al., 2012), the possible oncogenesis mechanisms are seemingly linked to calcium-dependent mitogenic signals of epidermal growth factor (Huang et al., 2004). As such, CACNG4 may serve as a potential predictor for tumorigenesis in cervical disease. Nevertheless, considering that the EBV-HPV coinfected cases involved in mRNA-seq all had CIN lesions in this study, we were unable to determine whether gene expression changes were simply associated with coinfection status but not cervical lesions. However, we noted that EBV-HPV coinfected women were more likely to be diagnosed with CIN2+ (p = 0.017) than women infected with HPV alone, and the precise mechanisms underlying the association between herpesvirus and HPV coinfection and carcinogenesis await further investigation in the future.

In conclusion, we reported transcriptional changes in pathological tissues from HIV-positive women with EBV and HPV coinfection. Coinfection is associated with precancerous lesions, which leads to changes in the gene expression of epithelial differentiation and development compared to normal tissues with HPV only. More important, several cancer-related DEGs were involved in CIN lesions with coinfection. These findings provide some evidence that EBV can act as a cofactor or mediator in HPV-related cervical cancer. Notably, specific genes or proteins, such as CACNG4, may serve as biomarkers that can risk stratify patients based on pathological changes in the cervix. However, the relatively small sample size is the limitation of our study. Further large-scale studies of the general population will allow us to validate and select specific biomarkers for early intervention in cervical cancer.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA701124.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethical Review Committees of National Cancer Center/Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

MF, RD, and YG performed the experiments and analyzed the data. HZ analyzed the data and assisted with writing the manuscript. YQ assisted with the project design and revised the manuscript. QL and FZ designed the project. MF, RD, QL, and FZ wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Science and Technology Major Project of Yunnan Province (2017ZF006), the CAMS Initiative for Innovative Medicine (2016-I2M-1-019), and the National Natural Sciences Foundation of China (81761128006 and 31670173).

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.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

For support in recruitment and patient sample collection, we thank all members of the AIDS Antiretroviral Therapy Department and Gynecological Department of the Third People’s Hospital of Kunming, China.

Supplementary Material

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

References

Ahdieh, L., Klein, R. S., Burk, R., Cu-Uvin, S., Schuman, P., Duerr, A., et al. (2001). Prevalence, Incidence, and Type-Specific Persistence of Human Papillomavirus in Human Immunodeficiency Virus (HIV)-Positive and HIV-Negative Women. J. Infect. Dis. 184, 682–690. doi: 10.1086/323081

PubMed Abstract | CrossRef Full Text | Google Scholar

Ammatuna, P., Giovannelli, L., Giambelluca, D., Mancuso, S., Rubino, E., Colletti, P., et al. (2000). Presence of Human Papillomavirus and Epstein-Barr Virus in the Cervix of Women Infected With the Human Immunodeficiency Virus. J. Med. Virol. 62, 410–415. doi: 10.1002/1096-9071(200012)62:4<410::AID-JMV3>3.0.CO;2-J

PubMed Abstract | CrossRef Full Text | Google Scholar

Aromseree, S., Pientong, C., Swangphon, P., Chaiwongkot, A., Patarapadungkit, N., Kleebkaow, P., et al. (2015). Possible Contributing Role of Epstein-Barr Virus (EBV) as a Cofactor in Human Papillomavirus (HPV)-Associated Cervical Carcinogenesis. J. Clin. Virol. 73, 70–76. doi: 10.1016/j.jcv.2015.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Blanco, R., Carrillo-Beltran, D., Osorio, J. C., Calaf, G. M., Aguayo, F. (2020). Role of Epstein-Barr Virus and Human Papillomavirus Coinfection in Cervical Cancer: Epidemiology, Mechanisms and Perspectives. Pathogens 9 (9), 685. doi: 10.3390/pathogens9090685

CrossRef Full Text | Google Scholar

Bosch, F. X., de Sanjose, S. (2007). The Epidemiology of Human Papillomavirus Infection and Cervical Cancer. Dis. Markers 23, 213–227. doi: 10.1155/2007/914823

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulla, R., Tripodo, C., Rami, D., Ling, G. S., Agostinis, C., Guarnotta, C., et al. (2016). C1q Acts in the Tumour Microenvironment as a Cancer-Promoting Factor Independently of Complement Activation. Nat. Commun. 7, 10346. doi: 10.1038/ncomms10346

PubMed Abstract | CrossRef Full Text | Google Scholar

Camargo, M., Soto-De Leon, S. C., Munoz, M., Sanchez, R., Pena-Herrera, D., Pineda-Pena, A. C., et al. (2014). Human Papillomavirus Detection in Women With and Without Human Immunodeficiency Virus Infection in Colombia. BMC Cancer 14, 451. doi: 10.1186/1471-2407-14-451

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameron, J. E., Dennis, D. C., Herrel, N. R., Chapple, A. G., Hagensee, M. E. (2020). Risk of Abnormal Cervical Cytology in HIV-Infected Women Testing Positive for Both Human Papillomavirus and Epstein-Barr Virus in Genital Tract Specimens. Cancer Causes Control 31, 365–375. doi: 10.1007/s10552-020-01287-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Cameron, J. E., Rositch, A. F., Vielot, N. A., Mugo, N. R., Kwatampora, J. K. L., Waweru, W., et al. (2018). Epstein-Barr Virus, High-Risk Human Papillomavirus and Abnormal Cervical Cytology in a Prospective Cohort of African Female Sex Workers. Sex Transm. Dis. 45, 666–672. doi: 10.1097/OLQ.0000000000000857

PubMed Abstract | CrossRef Full Text | Google Scholar

Castellsague, X., Bosch, F. X., Munoz, ,. N. (2002). Environmental Co-Factors in HPV Carcinogenesis. Virus Res. 89, 191–199. doi: 10.1016/S0168-1702(02)00188-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Centers for Disease Control and Prevention (2020) “Basic Information About HPV and Cancer”. Available at: https://www.cdc.gov/cancer/hpv/basic_info/index.htm.

Google Scholar

de Lima, M. A. P., Neto, P. J. N., Lima, L. P. M., Goncalves Junior, J., Teixeira Junior, A. G., Teodoro, I. P. P., et al. (2018). Association Between Epstein-Barr Virus (EBV) and Cervical Carcinoma: A Meta-Analysis. Gynecol. Oncol. 148, 317–328. doi: 10.1016/j.ygyno.2017.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Djaoud, Z., Guethlein, L. A., Horowitz, A., Azzi, T., Nemat-Gorgani, N., Olive, D., et al. (2017). Two Alternate Strategies for Innate Immunity to Epstein-Barr Virus: One Using NK Cells and the Other NK Cells and Gammadelta T Cells. J. Exp. Med. 214, 1827–1841. doi: 10.1084/jem.20161017

PubMed Abstract | CrossRef Full Text | Google Scholar

Forman, D., De Martel, C., Lacey, C. J., Soerjomataram, I., Lortet-Tieulent, J., Bruni, L., et al. (2012). Global Burden of Human Papillomavirus and Related Diseases. Vaccine 30 (Suppl 5), F12–F23. doi: 10.1016/j.vaccine.2012.07.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Goding, C. R., Arnheiter, H. (2019). MITF-The First 25 Years. Genes Dev. 33, 983–1007. doi: 10.1101/gad.324657.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Haq, R., Fisher, D. E. (2011). Biology and Clinical Relevance of the Micropthalmia Family of Transcription Factors in Human Cancer. J. Clin. Oncol. 29, 3474–3482. doi: 10.1200/JCO.2010.32.6223

PubMed Abstract | CrossRef Full Text | Google Scholar

Hebner, C. M., Laimins, L. A. (2006). Human Papillomaviruses: Basic Mechanisms of Pathogenesis and Oncogenicity. Rev. Med. Virol. 16, 83–97. doi: 10.1002/rmv.488

PubMed Abstract | CrossRef Full Text | Google Scholar

Higgins, C. D., Swerdlow, A. J., Macsween, K. F., Harrison, N., Williams, H., Mcaulay, K., et al. (2007). A Study of Risk Factors for Acquisition of Epstein-Barr Virus and its Subtypes. J. Infect. Dis. 195, 474–482. doi: 10.1086/510854

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J. B., Kindzelskii, A. L., Clark, A. J., Petty, H. R. (2004). Identification of Channels Promoting Calcium Spikes and Waves in HT1080 Tumor Cells: Their Apparent Roles in Cell Motility and Invasion. Cancer Res. 64, 2482–2489. doi: 10.1158/0008-5472.CAN-03-3501

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Z. Y., Liu, X., Ding, Y. Y., Zhang, Z. F., He, N. (2017). Cancer Risk Factors Among People Living With HIV/AIDS in China: A Systematic Review and Meta-Analysis. Sci. Rep. 7, 4890. doi: 10.1038/s41598-017-05138-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Joharinia, N., Faghihinejad, S., Seyedi, K., Farhadi, A., Hosseini, S. Y., Safaei, A., et al. (2020). Co-Existing of HSV1/2 or EBV Infection With the Presence of High-Risk HPV DNA in Cervical Lesions in the Southwest of Iran. Asian Pac. J. Cancer Prev. 21, 1459–1464. doi: 10.31557/APJCP.2020.21.5.1459

PubMed Abstract | CrossRef Full Text | Google Scholar

Kahla, S., Oueslati, S., Achour, M., Kochbati, L., Chanoufi, M. B., Maalej, M., et al. (2012). Correlation Between Ebv Co-Infection and HPV16 Genome Integrity in Tunisian Cervical Cancer Patients. Braz. J. Microbiol. 43, 744–753. doi: 10.1590/S1517-83822012000200039

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanwar, N., Carmine-Simmen, K., Nair, R., Wang, C., Moghadas-Jafari, S., Blaser, H., et al. (2020). Amplification of a Calcium Channel Subunit CACNG4 Increases Breast Cancer Metastasis. EBioMedicine 52, 102646. doi: 10.1016/j.ebiom.2020.102646

PubMed Abstract | CrossRef Full Text | Google Scholar

Khenchouche, A., Sadouki, N., Boudriche, A., Houali, K., Graba, A., Ooka, T., et al. (2013). Human Papillomavirus and Epstein-Barr Virus Co-Infection in Cervical Carcinoma in Algerian Women. Virol. J. 10, 340. doi: 10.1186/1743-422X-10-340

PubMed Abstract | CrossRef Full Text | Google Scholar

Kienka, T., Varga, M. G., Caves, J., Smith, J. S., Sivaraman, V. (2019). Epstein-Barr Virus, But Not Human Cytomegalovirus, is Associated With a High-Grade Human Papillomavirus-Associated Cervical Lesions Among Women in North Carolina. J. Med. Virol. 91, 450–456. doi: 10.1002/jmv.25336

PubMed Abstract | CrossRef Full Text | Google Scholar

Klymenko, T., Gu, Q., Herbert, I., Stevenson, A., Iliev, V., Watkins, G., et al. (2017). RNA-Seq Analysis of Differentiated Keratinocytes Reveals a Massive Response to Late Events During Human Papillomavirus 16 Infection, Including Loss of Epithelial Barrier Function. J. Virol. 91 (24), e01001–17. doi: 10.1128/JVI.01001-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuznetsova, E. B., Kekeeva, T. V., Larin, S. S., Zemliakova, V. V., Babenko, O. V., Nemtsova, M. V., et al. (2007). [Novel Methylation and Expression Markers Associated With Breast Cancer]. Mol. Biol. (Mosk) 41, 624–633. doi: 10.1134/S0026893307040061

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo, Y. M., Tein, M. S., Lau, T. K., Haines, C. J., Leung, T. N., Poon, P. M., et al. (1998). Quantitative Analysis of Fetal DNA in Maternal Plasma and Serum: Implications for Noninvasive Prenatal Diagnosis. Am. J. Hum. Genet. 62, 768–775. doi: 10.1086/301800

PubMed Abstract | CrossRef Full Text | Google Scholar

Marinho-Dias, J., Ribeiro, J., Monteiro, P., Loureiro, J., Baldaque, I., Medeiros, R., et al. (2013). Characterization of Cytomegalovirus and Epstein-Barr Virus Infection in Cervical Lesions in Portugal. J. Med. Virol. 85, 1409–1413. doi: 10.1002/jmv.23596

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, C. S., Berger, J. R., Mootoor, Y., Avdiushko, S. A., Zhu, H., Kryscio, R. J. (2006). High Prevalence of Multiple Human Herpesviruses in Saliva From Human Immunodeficiency Virus-Infected Persons in the Era of Highly Active Antiretroviral Therapy. J. Clin. Microbiol. 44, 2409–2415. doi: 10.1128/JCM.00256-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Monteith, G. R., Davis, F. M., Roberts-Thomson, S. J. (2012). Calcium Channels and Pumps in Cancer: Changes and Consequences. J. Biol. Chem. 287, 31666–31673. doi: 10.1074/jbc.R112.343061

PubMed Abstract | CrossRef Full Text | Google Scholar

Moscicki, A. B., Schiffman, M., Kjaer, S., Villa, L. L. (2006). Chapter 5: Updating the Natural History of HPV and Anogenital Cancer. Vaccine 24 (Suppl 3), S3/42–S3/51. doi: 10.1016/j.vaccine.2006.06.018

CrossRef Full Text | Google Scholar

Nayar, R., Wilbur, D. C. eds. (2015). The Bethesda System for Reporting Cervical Cytology. Definitions, Criteria, and Explanatory Notes. 3rd ed. Springer, Cham: Switzerland. (2015).

Google Scholar

Petrara, M. R., Cattelan, A. M., Zanchetta, M., Sasset, L., Freguja, R., Gianesin, K., et al. (2012). Epstein-Barr Virus Load and Immune Activation in Human Immunodeficiency Virus Type 1-Infected Patients. J. Clin. Virol. 53, 195–200. doi: 10.1016/j.jcv.2011.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos, N. B., Villanova, F. E., Andrade, P. M., Ribalta, J., Focchi, J., Otsuka, A. Y., et al. (2009). Epstein-Barr Virus Detection in Invasive and Pre-Invasive Lesions of the Uterine Cervix. Oncol. Rep. 21, 403–405. doi: 10.3892/or_00000236

PubMed Abstract | CrossRef Full Text | Google Scholar

Satterwhite, C. L., Torrone, E., Meites, E., Dunne, E. F., Mahajan, R., Ocfemia, M. C., et al. (2013). Sexually Transmitted Infections Among US Women and Men: Prevalence and Incidence Estimates 2008. Sex Transm. Dis. 40, 187–193. doi: 10.1097/OLQ.0b013e318286bb53

PubMed Abstract | CrossRef Full Text | Google Scholar

Silver, M. I., Paul, P., Sowjanya, P., Ramakrishna, G., Vedantham, H., Kalpana, B., et al. (2011). Shedding of Epstein-Barr Virus and Cytomegalovirus From the Genital Tract of Women in a Periurban Community in Andhra Pradesh, India. J. Clin. Microbiol. 49, 2435–2439. doi: 10.1128/JCM.02206-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Sixbey, J. W., Vesterinen, E. H., Nedrud, J. G., Raab-Traub, N., Walton, L. A., Pagano, J. S. (1983). Replication of Epstein-Barr Virus in Human Epithelial Cells Infected In Vitro. Nature 306, 480–483. doi: 10.1038/306480a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Slattery, M. L., Mullany, L. E., Sakoda, L. C., Wolff, R. K., Samowitz, W. S., Herrick, J. S. (2018). The MAPK-Signaling Pathway in Colorectal Cancer: Dysregulated Genes and Their Association With MicroRNAs. Cancer Inform 17, 1176935118766522. doi: 10.1177/1176935118766522

PubMed Abstract | CrossRef Full Text | Google Scholar

Wadowsky, R. M., Laus, S., Green, M., Webber, S. A., Rowe, D. (2003). Measurement of Epstein-Barr Virus DNA Loads in Whole Blood and Plasma by TaqMan PCR and in Peripheral Blood Lymphocytes by Competitive PCR. J. Clin. Microbiol. 41, 5245–5249. doi: 10.1128/JCM.41.11.5245-5249.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

WHO (2020). Global Strategy to Accelerate the Elimination of Cervical Cancer As A Public Health Problem. Available at: https://www.who.int/publications/i/item/9789240014107.

Google Scholar

WHO (2018) Report on Global Sexually Transmitted Infection Surveillance 2018. Available at: https://www.who.int/news-room/fact-sheets/detail/sexually-transmitted-infections-(stis).

Google Scholar

Wong, Y. F., Cheung, T. H., Tsao, G. S., Lo, K. W., Yim, S. F., Wang, V. W., et al. (2006). Genome-Wide Gene Expression Profiling of Cervical Cancer in Hong Kong Women by Oligonucleotide Microarray. Int. J. Cancer 118, 2461–2469. doi: 10.1002/ijc.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, X., Zhang, M., Hou, Y., Xu, L., Li, W., Zou, Z., et al. (2016). Single-Cell Analyses of Transcriptional Heterogeneity in Squamous Cell Carcinoma of Urinary Bladder. Oncotarget 7, 66069–66076. doi: 10.18632/oncotarget.11803

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, S., Zhao, X., Hu, S., Lu, J., Duan, X., Zhang, X., et al. (2019). Distribution of High-Risk Human Papillomavirus Genotype Prevalence and Attribution to Cervical Precancerous Lesions in Rural North China. Chin. J. Cancer Res. 31, 663–672. doi: 10.21147/j.issn.1000-9604.2019.04.10

PubMed Abstract | CrossRef Full Text | Google Scholar