• All
  • Article Title
  • Abstract
  • Keywords
  • Author
  • Institution
  • Open Access
    Original Article

    Prevalence of HPV, EBV, MMTV and SV40 co-infection in women with invasive breast cancer: as a potential new viral etiology

    Kaoutar Anouar Tadlaoui 1
    Amal Souiri 1,2
    Soukayna Alaoui Sosse 1
    Kawtar Nabil 1
    Mustapha Benhessou 1,3
    Moulay Mustapha Ennaji 1*

    Explor Med. 2025;6:1001304 DOI: https://doi.org/10.37349/emed.2025.1001304

    Received: October 13, 2024 Accepted: March 10, 2025 Published: April 16, 2025

    Academic Editor: Andrea Nicolini, University of Pisa, Italy

    This article belongs to the special issue Breast Cancer: Basic and Clinical Advances

    Abstract

    Aim:

    Breast cancer (BC) is the leading cause of female cancer-related death worldwide. The high incidence of BC has sparked interest in the viral agents role in its development. Identifying co-infection involving potential oncogenic viruses, such as human papillomavirus (HPV), Epstein-Barr virus (EBV), mouse mammary tumor virus (MMTV), and Simian virus 40 (SV40), may improve early detection and treatment strategies of BC. However, the involvement of these viral co-infections in invasive breast cancer (IBC) has not been elucidated.

    Methods:

    To investigate this eventual co-infection, screening of viral DNAs isolated from 120 fresh IBC tissue biopsies was performed using a polymerase chain reaction. Statistical analysis were conducted to assess the correlation between viral prevalence and IBC clinicopathological features.

    Results:

    Our findings revealed the prevalence of EBV (67.1%), HPV (45.7%), MMTV (72.9%), and SV40 (22.9%) in IBC samples. Co-infection rates were as follows: EBV/HPV (17.14%), EBV/MMTV (22.86%), EBV/SV40 (4.28%), and HPV/MMTV (11.43%). Triple infection with EBV, HPV, and MMTV was observed in 5.71% of cases. Statistically significant associations were identified between: EBV and histological grade, tumour size and stage, and progesterone receptor (PR) and estrogen receptor (ER) status; HPV and histological grade, tumour stage, and PR status; MMTV and histological type, as well as PR and ER status; and SV40 and histological grade and PR status.

    Conclusions:

    Although the prevalence of HPV, EBV, MMTV, and SV40 coinfection in Moroccan IBC patients is low, their potential synergistic role in breast carcinogenesis needs to be further investigated, in order to identify new etiological factors for BC.

    Keywords

    Invasive breast cancer, HPV, EBV, MMTV, SV40, coinfection, viral etiology, PCR

    Introduction

    Breast cancer (BC) is the most common cancer in women worldwide with alarming incidence and mortality rates. The etiology of BC is multifactorial, involving host genetic factors, viral infections, environmental factors, and hormonal factors related to estrogen and progesterone. Namely, prolonged use of oral contraceptives, early age of menarche, late age of menopause, and postmenopausal obesity [1]. However, there are other co-factors coupled to hormone-sensitive viruses, which induce malignant transformation of virally infected breast cells. Among the main oncoviral agents responsible for the development of BC are human papillomavirus (HPV), Epstein-Barr virus (EBV), mouse mammary tumor virus (MMTV), and Simian virus 40 (SV40), but there is no conclusive evidence for a causal role of co-infection with these viruses in BC [24].

    HPV and breast cancer

    HPVs are small non-enveloped particles of viruses containing a double-stranded DNA genome of approximately 8 kb [5]. They encompass 200 genotypes, some of which infect mucosal and skin epithelial cells, and about 50 genotypes cause genitourinary infections [6], leading to the development of specific cancers. HPV is the leading common cause of sexually transmitted genital infections worldwide and the main causative factor of cervical cancer, one of the most common cancers in women [7]. In Morocco, it is the second most common cancer in women after BC [8].

    HPV genotypes are classified according to their oncogenic potential: high-risk oncogenic HPV and low-risk oncogenic HPV. Numerous studies have shown that approximately 99.7% of cervical cancers are associated with high-risk oncogenic HPVs infections, particularly HPV 16 and 18 [9]. These two are frequently observed in high-grade intraepithelial lesions and neoplasia [10]. Thus, cell transformation and maintenance of the malignant phenotype are linked in particular to the expression of the two viral oncoproteins E6 and E7 of HPV as well as to the persistence of HPV infection, which can potentially progress to precancerous lesions and then to invasive malignancies, leading to carcinogenesis [7].

    The oncogenic HPV infections are also associated with other types of cancer, such as cancer of the penis, anus, vulva, vagina, neck, head, oropharynx [11, 12]. Moreover, several studies have highlighted the involvement of HPV in BC tissues and cell lines. Shortly after the discovery of the association between HPV and BC, evidence of the frequent involvement of HPV 16 in multiple invasive and metastatic BC was reported [13, 14]. HPV genotypes 16, 18, and 33 have also been identified in BC from very different populations [15]. Thus, it has been proposed that infection with high-risk oncogenic HPVs may drive cell invasion and metastasis in BC via Id-1, a family of helix-loop-helix transcription factors [14, 16].

    EBV and breast cancer

    EBV is a ubiquitous human virus, infecting 90% of the population, and is linked to several human malignancies, including BC [17]. EBV is a γ-herpesvirus whose genome, 184 kb in length of DNA, encodes about 100 genes [18].

    The Horiuchi team identified EBV genetic material for the first time in Japanese women with BC in 1994 [19]. To date, more than 50 studies on EBV and BC have followed. EBV has been reported to be involved in the development of invasive breast cancer (IBC) [20]. Fina and coworkers described a predominance of the EBV1 subtype in BC samples [21]. The link between BC and EBV provides information not only on the viral etiology of BC, but also on the early diagnosis, therapy, and prevention of this disease.

    Moreover, a current study indicates that EBV and HPV coexist in 32% of the same BC samples in women [2]. In the same study, viral DNA for EBV, HPV, and MMTV was detected in more than 50% of Australian BC cases [2].

    MMTV and breast cancer

    MMTV is a betaretrovirus, known as a milk-borne agent that causes mouse mammary cancer [22]. 70% of the complete MMTV-like virus genome found, in human BC tissues was sequenced and showed 91–99% homology with MMTV from murine mammary tumors [23, 24]. Although the MMTV gene sequence has been detected in human BC, its role in humans is still disputable [22]. However, evidence from Stewart and Chen (2022) [25] suggests that MMTV zoonosis may contribute to the geographic variation in BC incidence. A recent study reported the identification of MMTV, EBV, and HPV in the DNA of 50 fresh IBC samples, and the presence of more than one viral type in the same IBC sample [2].

    SV40 and breast cancer

    The SV40 virus has been added to the growing list of viruses suspected of being associated with BC. It is a monkey-origin polyomavirus, involved in several human tumors. SV40 carcinogenesis is based on the expression of the large viral antigen T (Tag) oncoprotein, which induces progression and apoptosis of infected cells by promoting the IGF-I (insulin-like growth factor-1) receptor, c-MET, and MAPK (mitogen-activated protein kinase) signaling pathways, as well as by inhibiting the two tumor suppressors p53 and pRB (retinoblastoma protein) [26].

    In 2002, Wong et al. [3] described novel associations between HPV, EBV, and SV40 in some human cancers, including the association of EBV and BC, but no relationship between SV40 and BC was reported. However, testing for SV40 in 109 BC specimens revealed the presence of SV40 Tag sequences in 22% of the cases [27].

    As well, Hachana et al. [28] identified SV40-positive DNA in BC tissue with a low rate. The prevalence of SV40 antibodies in serum samples from women with the combination is elaborate and substantial but cannot be considered conclusive.

    Therefore, to investigate the obvious presence or absence of viral coinfection in IBC, 120 IBC cases were screened for HPV, EBV, MMTV, and SV40 by molecular detection of DNA from each virus by polymerase chain reaction (PCR). In addition, statistical correlation analysis of the clinicopathologic characteristics of the patients and the candidate viruses was performed.

    Materials and methods

    Breast cancer patients and specimens

    Our study was performed on 120 fresh biopsies, 70 tumor tissue samples and 50 matched normal tissue samples, obtained from patients with histopathologically confirmed IBC after surgery in the Onco-Gynecology Department of the Mohammed IV Oncology Center, Casablanca, Morocco, over two years from January 2019 to December 2021. After collection, biopsies were immediately stored at –80°C until use. Patients who underwent chemotherapy and/or radiotherapy were excluded.

    The Ethics Committee for Biomedical Research of the Faculty of Medicine and Pharmacy of Casablanca, Morocco (Reference 3/2018 on 30.04.2018) approved the study protocol. In accordance with the ethical rules, free and informed consent was acquired from all patients who were included.

    Clinicopathological parameters of the recruited cases were collected according to STROCSS guidelines to enhance the quality and relevance of the information acquired [29].

    IBC patients clinicopathological characteristics are described in Table 1. The mean age of the patients was 45 (± 10.99) years (range: 29–78 years). The most common age range was ≥ 40 years, accounting for 74.3% of cases. Regarding histological findings, 38 patients (54.3%) had grade I or II, and 32 (45.7%) had grade III. 65 patients (92.9%) had invasive ductal carcinoma, and only five (7.1%) had invasive lobular carcinoma. In 52 cases (74.3%), the tumor size was superior to 5 cm. The majority of patients (57; 81.4%) were in stages T1 and T2. Regarding hormonal factors, 39 patients (55.7%) showed a positive progesterone receptor (PR) and 29 (41.4%) showed a positive estrogen receptor (ER).

    Clinicopathological features of invasive breast cancer patients (n = 70)

    FeaturesnPercentage (%)
    Age at diagnosis (years): mean age = 45 (± 10.99)
    < 401825.7
    ≥ 405274.3
    Histological grade
    I–II3854.3
    III3245.7
    Histological type
    Invasive ductal6592.9
    Invasive lobular57.1
    Tumor size (in cm)
    ≤ 51825.7
    > 55274.3
    Tumor stage
    T1–T25781.4
    T3–T41318.6
    Progesterone receptor (PR)
    Positive3955.7
    Negative3144.3
    Estrogen receptor (ER)
    Positive2941.4
    Negative4158.6
    Display full size

    n: number of patients

    DNA extraction

    DNA was extracted from frozen tissue subsamples of ≤ 25 mg taken from our IBC biopsies using the Purelink Invitrogen® Genomic DNA Minikit (Thermo Fisher Scientific, USA) according to the manufacturer’s instructions. DNA extracts were analyzed directly or stored at –20°C until use.

    Qualitative and quantitative DNA analysis

    The NanoDrop 2000 spectrophotometer (Technologies, Wilmington, DE, USA) was used to evaluate the concentration and purity of the extracted DNA. To assess the quality of the extraction, the integrity of the DNA extracts and the absence of PCR inhibitors, the β-globin gene used as a reference housekeeping gene was amplified by PCR using the primers PCO4 and GH20 (Table 2) as described previously [30].

    Primers sequences used for PCR amplification

    GenePrimers sequences (5–3)Tm (°C)PCR product size (bp)
    β-globinGH20-F 5'-GAAGAGCCAAGGACAGGTAC-3'54268
    PC04-R 5'-CAACTTCATCCACGTTCACC-3'
    HPVMY09-F 5'-CGTCCMARRGGAWACTGATC-3'55450
    MY11-R 5'-GCMCAGGGWCATAAYAATGG-3'
    GP5+-F 5'-TTTGTTACTGTGGTAGATACTAC-3'42150
    GP6+-R 5'-GAAAAATAAACTGTAAATCATATTC-3'
    EBVEBNA-2 OP-F 5'-GCGGGTGGAGGGAAAGG-3'58168
    EBNA-2 OP-R 5'-GTCAGCCAAGGGACGCG-3'
    EBNA-2 IP-F 5'-AGGCTGCCCACCCTGAGGAT-3'66168
    EBNA-2 IP-R 5'-GCCACCTGGCAGCCCTAAAG-3'
    SV40SVTAGP1-F 5'-TTAGCAATTCTGAAGGAAAGTCCTTG-3'52126
    SVTAGP3-F 5'-ACCTGTTTTGCTCAGAAG-3'
    MMTVMMTV1-F 5'-TGCGCCTTCCCTGACCAAGGG-3'54356
    MMTV2-R 5'-GTAACACAGGCAGATGTAGG-3'
    Display full size

    F: forward primer; R: reverse primer; PCR: polymerase chain reaction; Tm: melting temperature; HPV: human papillomavirus; EBV: Epstein-Barr virus; MMTV: mouse mammary tumor virus; SV40: Simian virus 40. M = A + C, R = A + G, W= A + T, Y = C + T

    Viral DNA detection

    DNA extracts from IBC samples were screened for viral DNA of each HPV, EBV, MMTV, and SV40. This analysis was performed by PCR amplification, using specific primer pairs (Table 2). Gene sequences of HPV and EBV were amplified by nested PCR using appropriate consensus primer sets, as previously described [31, 32]. However, a single-PCR round was performed to amplify MMTV-like and SV40 viral DNA using specific primer sets, as previously described [33, 34].

    Briefly, the PCR reaction mixture with a total volume of 25 µL consisted of 6.5 µL of ultrapure water, 2 µL of forward and reverse primers, and 12.5 µL of Taq PCR master mix (Vazyme Biotech, China) and, alternatively, 2 µL of concentrated DNA (100 ng/µL) for the first PCR, and 2 μL of the PCR product of the first reaction for nested-PCR. PCR amplification was performed on a Perkin Elmer 2400® thermal cycler (Scientific Support, Inc, Hayward, CA) under the following cycling conditions: initial denaturation at 94°C/10 min, 35 cycles of denaturation at 94°C/1 min, annealing of the primer at the corresponding melting temperature (Tm) for 1 min (Table 2), elongation at 72°C/1 min, followed by a final incubation at 72°C/10 min. A negative control containing the amplification reaction mixture without DNA was performed for each reaction.

    To evaluate the outcome of each PCR, all PCR products were analyzed by 2% agarose gel electrophoresis at 70 V for 1.5 h, stained with ethidium bromide, then visualized under UV light (SERVA Electrophoresis GmbH, Heidelberg, Germany).

    Statistical analysis

    Statistical analysis of the data was generated using Jamovi software (The Jamovi project, 2021). The correlation between DNA-positive viruses and IBC cases and the clinico-pathologic parameters was carried out by Chi-square test or Fisher’s exact test (when one of the theoretical numbers is less than 5). The difference was considered statistically significant when the P < 0.05.

    Results

    Oncoviral DNA prevalence

    In the present study, the β-globin gene was amplified in all 120 specimens (70 IBCs and 50 matched normal tissues), and the entire set of DNA samples was therefore suitable for further analysis. Our findings show that the prevalence of viral DNA sequences was significantly higher in IBC than in matched normal samples. Viral infection was detected in 61/70 (87.1%) of IBC tissues and 10/50 (20%) of matched normal tissues. Single and nested-PCR amplification of viral DNA sequences of each virus, namely EBV, HPV, MMTV, and SV40 revealed the detection of EBV DNA in (47/70; 67.1%) of IBC cases, HPV DNA in (32/70; 45.7%), MMTV DNA in (51/70; 72.9%), and SV40 DNA in (16/70; 22.9%). However, in the matched normal tissue samples, EBV, HPV, MMTV, and SV40 viral DNA was weakly detected in (6/50; 12%), (4/50; 8%), (8/50; 16%), and (0/50; 0%) of cases, respectively (Figure 1).

    Statistical analysis of the PCR results of the cases showed a significant difference between the viral prevalence of all viruses enrolled in IBC tissues and matched normal tissues (P < 0.001) (Figure 1).

    Prevalence of (A) EBV, (B) HPV, (C) MMTV, and (D) SV40 sequences in invasive breast cancer tumor tissues (n = 70) vs in matched normal tissues (n = 50). *P-value < 0.05 is considered significant. HPV: human papillomavirus; EBV: Epstein-Barr virus; MMTV: mouse mammary tumor virus; SV40: Simian virus 40

    Association between EBV, HPV, MMTV, and SV40 infections and clinico-pathological features of IBC patients

    Statistical analysis of the association between EBV, HPV, MMTV, and SV40 infection and clinico-pathological characteristics was evaluated and presented in Table 3. Regarding age at diagnosis, there is obviously no significant difference between the age and the IBC samples positive for all viruses studied. Statistical analysis clearly showed a significant difference between the presence of EBV viral DNA and histological grade (P = 0.022), tumor size and tumor stage (P < 0.001), positive PR (P = 0.003), and ER (P < 0.001). Thus, HPV infection is significantly associated with histological grade (P < 0.001), tumor stage (P = 0.004), and positive PR (P < 0.001). As well, MMTV infection is significantly associated with the type of invasive ductal carcinoma (P < 0.001), and the presence of PR and ER (P < 0.001). Finally, statistical analysis revealed that SV40 was also significantly associated with histological grade (P < 0.001), and positive PR (P = 0.023).

    Distribution of EBV, HPV, MMTV, and SV40 positive cases according to clinicopathological features of IBC patients (n = 70)

    FeaturesTotal
    n (%)
    EBV positiveHPV positiveMMTV positiveSV40 positive
    n (%)P-valuen (%)P-valuen (%)P-valuen (%)P-value
    Age at diagnosis (years)
    < 4018 (25.71)12 (66.7)0.9605 (27.8)0.07616 (88.9)0.1236 (33.3)0.219
    ≥ 4052 (74.29)35 (67.3)27 (51.9)35 (67.3)10 (19.2)
    Histological grade
    I–II38 (54.29)30 (78.9)0.022*26 (68.4)< 0.001*30 (78.9)0.21216 (42.1)< 0.001*
    III32 (45.71)17 (53.1)6 (18.8)21 (65.6)0 (0)
    Histological type
    Invasive ductal65 (92.86)45 (69.2)0.32230 (46.2)> 0.9950 (76.9)< 0.001*16 (24.6)0.207
    Invasive lobular5 (7.14)2 (40.0)2 (40.0)1 (20.0)0 (0)
    Tumor size (in cm)
    ≤ 518 (25.71)4 (22.2)< 0.001*5 (27.8)0.07615 (83.3)0.0712 (11.1)0.209
    > 552 (74.29)43 (82.7)27 (51.9)36 (69.2)14 (26.9)
    Tumor stage
    T1–T257 (81.43)45 (78.9)< 0.001*21 (36.8)0.004*41 (71.9)> 0.9913 (22.8)> 0.99
    T3–T413 (18.57)2 (15.4)11 (84.6)10 (76.9)3 (23.1)
    Progesterone receptor (PR)
    Positive39 (55.71)32 (82.1)0.003*25 (64.1)< 0.001*35 (89.7)< 0.001*13 (33.3)0.023*
    Negative31 (44.29)15 (48.4)7 (22.6)16 (51.6)3 (9.7)
    Estrogen receptor (ER)
    Positive29 (41.43)26 (89.7)< 0.001*13 (44.8)0.90113 (44.8)< 0.001*7 (24.1)0.830
    Negative41 (58.57)21 (51.2)19 (46.3)38 (92.7)9 (22.0)
    Display full size

    *P < 0.05 is considered significant. HPV: human papillomavirus; EBV: Epstein-Barr virus; MMTV: mouse mammary tumor virus; SV40: Simian virus 40; IBC: invasive breast cancer

    Viral co-infection of EBV, HPV, MMTV and SV40 in IBC

    Our findings demonstrated that EBV, HPV, MMTV, and SV40 viral sequences may all be present as individual or multiple viruses in IBC and in some normal breast tissues, as shown in Figure 2. The results showed that mono-infection with EBV, HPV, MMTV, or SV40 was detected in 21.43%, 11.43%, 32.86%, and 18.57% of cases respectively. However, co-infection of EBV and HPV was detected in 12 (17.14%), EBV and MMTV co-infection in 16 (22.86%), co-infection of EBV and SV40 in 3 (4.28%), and co-infection of HPV and MMTV in 8 (11.43%) of the IBC specimens (Figure 2). Triple infection with EBV, HPV, and MMTV was detected in only 4 (5.71%) of the samples. There was no viral coinfection of all viral sequences in any of the IBC samples. No significant statistical correlation was disclosed between the prevalence of co-infection and IBC disease findings.

    Diagram of the prevalence of EBV, HPV, MMTV, and SV40 co-infection in IBC tumor tissues. HPV: human papillomavirus; EBV: Epstein-Barr virus; MMTV: mouse mammary tumor virus; SV40: Simian virus 40; IBC: invasive breast cancer

    Discussion

    To the best of our knowledge, this is the first study to investigate the presence/co-presence of HPV, EBV, MMTV-like, and SV40 viruses in IBC tumor tissue in Morocco. Our study revealed that 87.1% of samples were positive for these oncoviruses.

    Of particular interest, the prevalence of HPV in BC cases differed between countries worldwide. In Morocco, previous studies showed a prevalence of HPV DNA in around 25% of BC tumors [35], whereas our data showed a higher percentage of 45.7%. In fact, a recent assessment of HPV-related cancer incidence in Morocco, published by HPV center information, showed that the high level was attributed to BC at 63% compared to other types of non-cervical cancers [36]. This finding is consistent with the increasing frequency of HPV DNA in BC patients worldwide.

    Various research on the investigation of HPV in BC revealed that 29% of BC samples carried high-risk HPV genotypes [37]. According to a recent study, in 855 cases of BC from the TCGA database, there were 30 (3.5%) low-risk and 20 (2.3%) high-risk HPV subtypes. The high-risk genotypes included HPV18 (48%), HPV113 (24%), HPV16 (10%), and HPV52 (10%) [38]. Most HPV types detected in BC studies belong to the high-risk HPV genotypes. However, a low-risk type, HPV6, was presented in co-infection with another high-risk type in a study of 700 samples from Thai women with BC [39]. Khan et al. [40] also detected HPV6 in 46% of BC samples of Japanese patients. de Villiers’ study reported that the most prevalent type in IBC was HPV11, followed by HPV6, while other types were detected at a low percentage, including HPV16, 23, 27, 57, 3, 15, 24, and 87 [41]. Moreover, several studies indicate that high-risk HPV types 16, 18, and 33 have been identified in BC patients from a diverse range of demographic groups. According to a PCR cohort research, HPV18 was found in 55% of IBC specimens, followed by HPV16 in 13% [38]. According to Hennig et al. [42], 46% of women with BC tested positive for HPV-16. Damin et al. [43] showed that Brazilian women with BC were mainly carriers of HPV16 (56%) and HPV18 (40%). 72 patients from China and Japan showed HPV33-positive BC [44]. Kan et al. [45] showed that 48% of BC tissues tested positive for HPV18. According to research by Fernandes et al. [46], HPV51 was found in 30.77% of BC samples, with HPV genotypes 18 and 33 following closely behind with 23.08%.

    In fact, according to research on HPV16 DNA in BC samples, HPV16 is more frequently found in invasive and metastatic BC than in in-situ BC [14]. Id-1, a family of helix-loop-helix transcription factors, has been proposed as a potential mechanism through which high-risk HPV infection can cause cell invasion and metastasis in BC [14, 16].

    On the other hand, a single EBV infection was detected in 67.1% of BC tissues and in 12% of normal tissues. Various EBV genes have been identified in BC in a wide range of countries. Richardson et al. [47] have summarized the results of these studies in a meta-analysis. The percentage of EBV positives detected by molecular methods in BC patients ranged from 24% to 100%, while in normal controls it ranged from 0% to 33%. However, the percentage is higher when tested by serology [48]. In Morocco, the high incidence of EBV infection in BC cases, as reported by Gihbid et al. [49] in a recent study, supports the idea that EBV infection may play a role in the promotion or progression of BC as an etiologic factor or cofactor in the oncogenic mechanism that raises the risk of various BC subtypes. The disparate outcomes might be brought on by sample source restrictions or geographic variations in the prevalence of EBV infection and biopsies. As a result, the detection rate varied according to the patient’s ethnicity, the study’s approach, and the sample’s demographics [50].

    MMTV-like involvement in BC has also been evaluated. In our current study, viral sequences were detected in 72.9% of cases. Two previous Moroccan studies showed relatively independent frequencies of MMTV-like in BC (57.14% and 56%, respectively) [51, 52]. Indeed, MMTV-like infection has been shown to be strongly associated with the inflammatory phenotype of BC. Unfortunately, the inflammatory status of our sample was not investigated. Overall, the prevalence of MMTV-like virus in BC is variable worldwide, ranging from 0 to 74% of cases [53].

    Further, SV40 DNA sequences were detected in 22.9% of tumors and were not found in the corresponding non-tumor tissues. The positive rate was in line with the previous study conducted by Hachana et al. [28]. However, Martini et al. [54] concluded that the low prevalence of SV40 in samples from BC women and healthy women suggests no association between SV40 infection and BC. SV40 is likely a passenger not involved in the development of BC, in addition to possible endogenous contamination by SV40 sequences. Despite this, numerous epidemiological studies do not rule out a link between SV40 and cancer prevalence, since the virus can transform human cells.

    Otherwise, the prevalence of viral genomes in tumors was analyzed for possible association with clinicopathological parameters. In our overall population, none of the four viruses was correlated with age. There was, however, a positive association between almost all viruses and histological grade, except for MMTV. Furthermore, our results demonstrated the importance of tumor stage for HPV and EBV. On the other hand, tumor size was only associated with EBV.

    Interestingly, our results showed a trend toward a significant association between MMTV and histological type, as well as PRs and ERs. Indeed, MMTV was detected in 86.2% of invasive ductal carcinoma. Consistent with our study, Mazzanti et al. [55] reported the presence of an MMTV-like env sequence in 82% of ductal carcinomas in-situ. In contrast to our results, a study conducted in Morocco by Slaoui et al. [51], revealed no significant correlation between the presence of MMTV-like sequences and histological subtype and hormone receptors. In another study, Taneja et al. [22], similarly found that the presence of PR correlated significantly with MMTV-like positive BC, concluding that steroid hormones and MMTV may play a causal role in the development of BC.

    Furthermore, our results showed a significant correlation between EBV and histological grade, tumor size, as well as PR and ER expression. In agreement with our results, Zhang et al. [50] found that EBV infection correlated significantly with tumor size and ER and PR (all P < 0.05). In contrast, they found no correlation with histological grade [50]. ER/PR receptors play a well-established role in the management and prognosis of BC [56].

    Our findings analysis corroborates those of Kroupis et al. [57] and Kouloura et al. [58], who revealed a statistically significant relationship between HPV-positive BCs and PR-positive expression (P = 0.041) and histological grade (P = 0.005), respectively. However, they disagree with the outcomes of the research conducted by Elagali et al. [59].

    Moreover, only 4 (5.71%) of the IBC samples in our investigation had coexistence of EBV, HPV, and MMTV infection, and there was no significant correlation between these findings and the clinicopathological characteristics. This result differs slightly from that of research conducted by Gupta et al. [60], who found triple co-infection in just 3% of samples of human triple-negative BC. The findings of Naushad et al. [4] also showed that 2.4% of samples were positive for all three viruses in Pakistani BC patients. Consistently, another study by Glenn et al. [2] in Australia, showed that the rate of co-infection of these three viruses in BC appears to be extremely low.

    Moreover, available evidence clearly indicates that certain viruses play a role in the development of cancer. Above all, data on the maintenance of malignancy after loss of viral “presence” are enough to support the theory of viral cancer emergence. The hit-and-run theory will be put to the test by observing the results of long-term vaccination. According to this theory, viruses may play a triggering role at an early stage of carcinogenesis, such as cancer initiation (the “Hit”), and then the viral genome can be completely lost after the host cell has acquired many mutations and cancer has reached a developed stage (the “Run”). Alternately, viral infection at an advanced stage of cancer development could promote oncogenic features including growth, invasiveness, angiogenesis, and metastasis [61]. This theory explains the multi-stage manner in which virus-induced BC develops and progresses.

    Furthermore, it has been demonstrated that HPV-cancer cells can only proliferate when the E6 and E7 oncogenes are expressed in cases of cervical, head, and neck cancer [62]. Other cancers such as BC and prostate cancer may be experiencing “hit-and-run” occurrences [63, 64]. BC development is initiated by E6 and E7 oncogenes, but as mutations are generated during carcinogenesis, the viral genes are eliminated since they are no longer needed [64]. In addition, Arbach et al. [61] reviewed several studies showing that many EBV-positive tumor cells do not contain EBV genomes, that the carcinogenic phenotype is unaffected by the loss of EBV episomes, and that BC exhibit genomic heterogeneity.

    If the “hit-and-run” theory proves true, viruses could be responsible for more cancers than previously thought, and vaccination programs could therefore help prevent many cancers.

    Although the prevalence of HPV, EBV, and MMTV co-infection is low, the synergistic carcinogenic impact of the different viruses in BC pathogenesis remains possible. We propose that the initiation and advancement of various malignancies, including BC, can be caused by interactions between the oncoproteins of HPVs (E6 and E7), and those of EBV (LMP1 and/or EBNA1) via EMT pathway, as well as with MMTV-like virus [60, 65, 66]. By promoting IFN signaling and APOBEC-mediated mutagenesis in BC, the MMTV-like virus contributed to the growth of tumors and the initiation and progression of BC [60].

    Furthermore, it has been reported that a chronic inflammatory state brought on by the co-infection of another virus may enhance EBV, HPV, or MMTV-induced carcinogenesis. According to evidence on coinfection, virus-infected cells may provide a favorable environment for subsequent viral infection by producing inflammatory cytokines. These conclusions contribute to a better knowledge of the immune system’s limitations in order to create more effective vaccinations to prevent viral infections responsible for BC [4].

    The findings of these investigations called for additional research using a large sample size and appropriate cell biologic tools, in order to clarify the likely causative roles of these viruses in BC as well as those of other viruses. Moreover, a potential limitation of our study is that the selected MMTV primers may not entirely capture all exogenous MMTV variants, which could lead to under-detection in some samples, and may require alternative primer strategies for future investigations.

    The data obtained on single and multiple HPV, EBV, MMTV, and SV40 infections in IBC tissues consolidate the idea of a viral etiology of IBC, but do not imply a causal link between co-infection of the four viruses and the development of IBC. Further studies on the roles and the action mechanisms of these viruses are required to support the perspective of a viral etiology of BC.

    Abbreviations

    BC:

    breast cancer

    EBV:

    Epstein-Barr virus

    ER:

    estrogen receptor

    HPV:

    human papillomavirus

    IBC:

    invasive breast cancer

    MMTV:

    mouse mammary tumor virus

    PCR:

    polymerase chain reaction

    PR:

    progesterone receptor

    SV40:

    Simian virus 40

    Tag:

    antigen T

    Declarations

    Acknowledgments

    Special acknowledgment for support of Faculty of Sciences and Techniques, Mohammedia, University Hassan II of Casablanca and Ministry of National Education, Vocational Training, Higher Education and Scientific Research in Morocco.

    Author contributions

    KAT: Conceptualization, Data curation, Methodology, Formal analysis, Software, Investigation, Project administration, Writing—original draft. AS: Writing—review & editing, Project administration. SAS: Visualization, Validation. KN: Investigation. MB: Resources. MME: Project administration, Supervision. All authors read and approved the submitted version.

    Conflicts of interest

    The authors declare that they have no conflicts of interest.

    Ethical approval

    The study protocol was approved by the Ethics Committee for Biomedical Research of the Faculty of Medicine and Pharmacy of Casablanca, Morocco (Reference 3/2018 on 30.04.2018).

    Consent to participate

    Free and informed consent was obtained from all patients involved in this research and the confidentiality of their personal information was well-respected according to ethical rules.

    Consent to publication

    Not applicable.

    Availability of data and materials

    The datasets generated and/or analyzed in the course of the present study are available from the corresponding author upon reasonable request.

    Funding

    Not applicable.

    Copyright

    © The Author(s) 2025.

    Publisher’s note

    Open Exploration maintains a neutral stance on jurisdictional claims in published institutional affiliations and maps. All opinions expressed in this article are the personal views of the author(s) and do not represent the stance of the editorial team or the publisher.

    References

    Horn J, Vatten LJ. Reproductive and hormonal risk factors of breast cancer: a historical perspective. Int J Womens Health. 2017;9:26572. [DOI] [PubMed] [PMC]
    Glenn WK, Heng B, Delprado W, Iacopetta B, Whitaker NJ, Lawson JS. Epstein-Barr virus, human papillomavirus and mouse mammary tumour virus as multiple viruses in breast cancer. PLoS One. 2012;7:e48788. [DOI] [PubMed] [PMC]
    Wong M, Pagano JS, Schiller JT, Tevethia SS, Raab-Traub N, Gruber J. New associations of human papillomavirus, Simian virus 40, and Epstein-Barr virus with human cancer. J Natl Cancer Inst. 2002;94:18326. [DOI] [PubMed]
    Naushad W, Surriya O, Sadia H. Prevalence of EBV, HPV and MMTV in Pakistani breast cancer patients: A possible etiological role of viruses in breast cancer. Infect Genet Evol. 2017;54:2307. [DOI] [PubMed]
    Doorslaer KV, Chen Z, Bernard H, Chan PKS, DeSalle R, Dillner J, et al. ICTV Virus Taxonomy Profile: Papillomaviridae. J Gen Virol. 2018;99:98990. [DOI] [PubMed] [PMC]
    Graham SV. The human papillomavirus replication cycle, and its links to cancer progression: a comprehensive review. Clin Sci (Lond). 2017;131:220121. [DOI] [PubMed]
    Tadlaoui KA, Hassou N, Bennani B, Ennaji MM. Emergence of oncogenic high-risk human papillomavirus types and cervical cancer. In: Ennaji MM, editor. Emerging and Reemerging Viral Pathogens. Elsevier; 2020. pp.539–70. [DOI]
    Tadlaoui KA, Ennaji MM. The molecular mechanism of novel oncogenes dysregulating signaling pathways associated with cervical carcinoma. In: Ennaji MM, editor. Immunological Implications and Molecular Diagnostics of Genitourinary Cancer. Elsevier; 2020. pp. 19–31. [DOI]
    Belglaiaa E. Génotypage moléculaire des papillomavirus humains chez des femmes à risque de cancer du col de l'utérus: implication pour le dépistage et la prévention [dissertation]. Agadir: Université Ibn Zohr; 2015. French.
    Omland T, Lie KA, Akre H, Sandlie LE, Jebsen P, Sandvik L, et al. Recurrent respiratory papillomatosis: HPV genotypes and risk of high-grade laryngeal neoplasia. PLoS One. 2014;9:e99114. [DOI] [PubMed] [PMC]
    de Martel C, Plummer M, Vignat J, Franceschi S. Worldwide burden of cancer attributable to HPV by site, country and HPV type. Int J Cancer. 2017;141:66470. [DOI] [PubMed] [PMC]
    Serrano B, Brotons M, Bosch FX, Bruni L. Epidemiology and burden of HPV-related disease. Best Pract Res Clin Obstet Gynaecol. 2018;47:1426. [DOI] [PubMed]
    Lonardo AD, Venuti A, Marcante ML. Human papillomavirus in breast cancer. Breast Cancer Res Treat. 1992;21:95100. [DOI] [PubMed]
    Yasmeen A, Bismar TA, Kandouz M, Foulkes WD, Desprez P, Moustafa AA. E6/E7 of HPV type 16 promotes cell invasion and metastasis of human breast cancer cells. Cell Cycle. 2007;6:203842. [DOI] [PubMed]
    Amarante MK, Watanabe MAE. The possible involvement of virus in breast cancer. J Cancer Res Clin Oncol. 2009;135:32937. [DOI] [PubMed]
    Akil N, Yasmeen A, Kassab A, Ghabreau L, Darnel AD, Moustafa AA. High-risk human papillomavirus infections in breast cancer in Syrian women and their association with Id-1 expression: a tissue microarray study. Br J Cancer. 2008;99:4047. [DOI] [PubMed] [PMC]
    Lawson JS, Günzburg WH, Whitaker NJ. Viruses and human breast cancer. Future Microbiol. 2006;1:3351. [DOI] [PubMed]
    Muhammad U, Yasir H, Mukhtiar A. Does Epstein–Barr Virus Participate in the Development of Breast Cancer? A Brief and Critical Review with Molecular Evidences. Biomedical and Biotechnology Research Journal (BBRJ). 2020;4:28592. [DOI]
    Horiuchi K, Mishima K, Ohsawa M, Aozasa K. Carcinoma of stomach and breast with lymphoid stroma: localisation of Epstein-Barr virus. J Clin Pathol. 1994;47:53840. [DOI] [PubMed] [PMC]
    Labrecque LG, Barnes DM, Fentiman IS, Griffin BE. Epstein-Barr virus in epithelial cell tumors: a breast cancer study. Cancer Res. 1995;55:3945. [PubMed]
    Fina F, Romain S, Ouafik L, Palmari J, Ayed FB, Benharkat S, et al. Frequency and genome load of Epstein-Barr virus in 509 breast cancers from different geographical areas. Br J Cancer. 2001;84:78390. [DOI] [PubMed] [PMC]
    Taneja P, Frazier DP, Kendig RD, Maglic D, Sugiyama T, Kai F, et al. MMTV mouse models and the diagnostic values of MMTV-like sequences in human breast cancer. Expert Rev Mol Diagn. 2009;9:42340. [DOI] [PubMed] [PMC]
    Melana SM, Nepomnaschy I, Sakalian M, Abbott A, Hasa J, Holland JF, et al. Characterization of viral particles isolated from primary cultures of human breast cancer cells. Cancer Res. 2007;67:89605. [DOI] [PubMed]
    Lawson JS, Glenn WK, Salmons B, Ye Y, Heng B, Moody P, et al. Mouse mammary tumor virus-like sequences in human breast cancer. Cancer Res. 2010;70:357685. [DOI] [PubMed]
    Stewart AFR, Chen H. Revisiting the MMTV Zoonotic Hypothesis to Account for Geographic Variation in Breast Cancer Incidence. Viruses. 2022;14:559. [DOI] [PubMed] [PMC]
    Shang D, Zhou T, Zhuang X, Wu Y, Liu H, Tu Z. Molecular dissection on inhibition of Ras-induced cellular senescence by small t antigen of SV40. Cell Mol Life Sci. 2022;79:242. [DOI] [PubMed] [PMC]
    Alibek K, Kakpenova A, Mussabekova A, Sypabekova M, Karatayeva N. Role of viruses in the development of breast cancer. Infect Agent Cancer. 2013;8:32. [DOI] [PubMed] [PMC]
    Hachana M, Trimeche M, Ziadi S, Amara K, Korbi S. Evidence for a role of the Simian Virus 40 in human breast carcinomas. Breast Cancer Res Treat. 2009;113:4358. [DOI] [PubMed]
    Agha R, Abdall-Razak A, Crossley E, Dowlut N, Iosifidis C, Mathew G, et al.; STROCSS Group. STROCSS 2019 Guideline: Strengthening the reporting of cohort studies in surgery. Int J Surg. 2019;72:15665. [DOI] [PubMed]
    Saiki RK, Bugawan TL, Horn GT, Mullis KB, Erlich HA. Analysis of enzymatically amplified beta-globin and HLA-DQ alpha DNA with allele-specific oligonucleotide probes. Nature. 1986;324:1636. [DOI] [PubMed]
    Lee SH, Vigliotti VS, Vigliotti JS, Pappu S. Validation of human papillomavirus genotyping by signature DNA sequence analysis. BMC Clin Pathol. 2009;9:3. [DOI] [PubMed] [PMC]
    van Baarle D, Hovenkamp E, Kersten MJ, Klein MR, Miedema F, van Oers MH. Direct Epstein-Barr virus (EBV) typing on peripheral blood mononuclear cells: no association between EBV type 2 infection or superinfection and the development of acquired immunodeficiency syndrome-related non-Hodgkin’s lymphoma. Blood. 1999;93:394955. [PubMed]
    Pogo BG, Melana SM, Holland JF, Mandeli JF, Pilotti S, Casalini P, et al. Sequences homologous to the mouse mammary tumor virus env gene in human breast carcinoma correlate with overexpression of laminin receptor. Clin Cancer Res. 1999;5:210811. [PubMed]
    Li R, Mannon RB, Kleiner D, Tsokos M, Bynum M, Kirk AD, et al. BK virus and SV40 co-infection in polyomavirus nephropathy. Transplantation. 2002;74:1497504. [DOI] [PubMed]
    ElAmrani A, Gheit T, Benhessou M, McKay-Chopin S, Attaleb M, Sahraoui S, et al. Prevalence of mucosal and cutaneous human papillomavirus in Moroccan breast cancer. Papillomavirus Res. 2018;5:1505. [DOI] [PubMed] [PMC]
    Bruni L, Albero G, Serrano B, Mena M, Collado JJ, Gómez D, et al. Human Papillomavirus and Related Diseases Report [Internet]. ICO/IARC Information Centre on HPV and Cancer; © 2025 [cited 2024 Oct 13]. Available from: https://hpvcentre.net/statistics/reports/MAR.pdf
    Wang T, Chang P, Wang L, Yao Q, Guo W, Chen J, et al. The role of human papillomavirus infection in breast cancer. Med Oncol. 2012;29:4855. [DOI] [PubMed]
    Lawson JS, Glenn WK, Salyakina D, Delprado W, Clay R, Antonsson A, et al. Human Papilloma Viruses and Breast Cancer. Front Oncol. 2015;5:277. [DOI] [PubMed] [PMC]
    Ngamkham J, Karalak A, Chaiwerawattana A, Sornprom A, Thanasutthichai S, Sukarayodhin S, et al. Prevalence of Human Papillomavirus Infection in Breast Cancer Cells from Thai Women. Asian Pac J Cancer Prev. 2017;18:183945. [DOI] [PubMed] [PMC]
    Khan NA, Castillo A, Koriyama C, Kijima Y, Umekita Y, Ohi Y, et al. Human papillomavirus detected in female breast carcinomas in Japan. Br J Cancer. 2008;99:40814. [DOI] [PubMed] [PMC]
    de Villiers EM, Sandstrom RE, Hausen Hz, Buck CE. Presence of papillomavirus sequences in condylomatous lesions of the mamillae and in invasive carcinoma of the breast. Breast Cancer Res. 2005;7:R111. [DOI] [PubMed] [PMC]
    Hennig EM, Suo Z, Thoresen S, Holm R, Kvinnsland S, Nesland JM. Human papillomavirus 16 in breast cancer of women treated for high grade cervical intraepithelial neoplasia (CIN III). Breast Cancer Res Treat. 1999;53:12135. [DOI] [PubMed]
    Damin APS, Karam R, Zettler CG, Caleffi M, Alexandre COP. Evidence for an association of human papillomavirus and breast carcinomas. Breast Cancer Res Treat. 2004;84:1317. [DOI] [PubMed]
    Yu Y, Morimoto T, Sasa M, Okazaki K, Harada Y, Fujiwara T, et al. HPV33 DNA in premalignant and malignant breast lesions in Chinese and Japanese populations. Anticancer Res. 1999;19:505761. [PubMed]
    Kan C, Iacopetta BJ, Lawson JS, Whitaker NJ. Identification of human papillomavirus DNA gene sequences in human breast cancer. Br J Cancer. 2005;93:9468. [DOI] [PubMed] [PMC]
    Fernandes A, Bianchi G, Feltri AP, Pérez M, Correnti M. Presence of human papillomavirus in breast cancer and its association with prognostic factors. Ecancermedicalscience. 2015;9:548. [DOI] [PubMed] [PMC]
    Richardson AK, Currie MJ, Robinson BA, Morrin H, Phung Y, Pearson JF, et al. Cytomegalovirus and Epstein-Barr virus in breast cancer. PLoS One. 2015;10:e0118989. [DOI] [PubMed] [PMC]
    Lawson JS, Salmons B, Glenn WK. Oncogenic Viruses and Breast Cancer: Mouse Mammary Tumor Virus (MMTV), Bovine Leukemia Virus (BLV), Human Papilloma Virus (HPV), and Epstein-Barr Virus (EBV). Front Oncol. 2018;8:1. [DOI] [PubMed] [PMC]
    Gihbid A, Amrani AE, Mouh FZ, Gheit T, Benhessou M, Amrani M, et al. Prevalence of Polyomaviruses and Herpesviruses in Moroccan Breast Cancer. Pathogens. 2023;12:640. [DOI] [PubMed] [PMC]
    Zhang W, Zhu G, Xu W, Zhang W, Wang H, Wang Y, et al. Association of PD-1/PD-L1 expression and Epstein--Barr virus infection in patients with invasive breast cancer. Diagn Pathol. 2022;17:61. [DOI] [PubMed] [PMC]
    Slaoui M, Mzibri ME, Razine R, Qmichou Z, Attaleb M, Amrani M. Detection of MMTV-Like sequences in Moroccan breast cancer cases. Infect Agent Cancer. 2014;9:37. [DOI] [PubMed] [PMC]
    Saif I, Bouziyane A, Mzibri ME, Kerroumi ME, Benhessou M, Ennaji MM. Prevalence of Human Papillomavirus and Mouse Mammary Tumor Virus Like DNAs in Tumors from Moroccan Breast Cancer Patients. Indian J Gynecol Oncolog. 2022;20:38. [DOI]
    Amarante MK, de Sousa Pereira N, Vitiello GAF, Watanabe MAE. Involvement of a mouse mammary tumor virus (MMTV) homologue in human breast cancer: Evidence for, against and possible causes of controversies. Microb Pathog. 2019;130:28394. [DOI] [PubMed]
    Martini F, Mazzoni E, Corallini A, Taronna A, Querzoli P, Magri E, et al. Breast cancer and simian virus 40 infection. Epidemiology. 2013;24:4645. [DOI] [PubMed]
    Mazzanti CM, Hamad MA, Fanelli G, Scatena C, Zammarchi F, Zavaglia K, et al. A mouse mammary tumor virus env-like exogenous sequence is strictly related to progression of human sporadic breast carcinoma. Am J Pathol. 2011;179:208390. [DOI] [PubMed] [PMC]
    Hu R, Dawood S, Holmes MD, Collins LC, Schnitt SJ, Cole K, et al. Androgen receptor expression and breast cancer survival in postmenopausal women. Clin Cancer Res. 2011;17:186774. [DOI] [PubMed] [PMC]
    Kroupis C, Markou A, Vourlidis N, Dionyssiou-Asteriou A, Lianidou ES. Presence of high-risk human papillomavirus sequences in breast cancer tissues and association with histopathological characteristics. Clin Biochem. 2006;39:72731. [DOI] [PubMed]
    Kouloura A, Nicolaidou E, Misitzis I, Panotopoulou E, Kassiani T, Smyrniotis V, et al. HPV infection and breast cancer. Results of a microarray approach. Breast. 2018;40:1659. Erratum in: Breast. 2018;41:136. [DOI] [PubMed]
    Elagali AM, Suliman AA, Altayeb M, Dannoun AI, Parine NR, Sakr HI, et al. Human papillomavirus, gene mutation and estrogen and progesterone receptors in breast cancer: a cross-sectional study. Pan Afr Med J. 2021;38:43. [DOI] [PubMed] [PMC]
    Gupta I, Ulamec M, Peric-Balja M, Ramic S, Moustafa AA, Vranic S, et al. Presence of high-risk HPVs, EBV, and MMTV in human triple-negative breast cancer. Hum Vaccin Immunother. 2021;17:445766. [DOI] [PubMed] [PMC]
    Arbach H, Viglasky V, Lefeu F, Guinebretière J, Ramirez V, Bride N, et al. Epstein-Barr virus (EBV) genome and expression in breast cancer tissue: effect of EBV infection of breast cancer cells on resistance to paclitaxel (Taxol). J Virol. 2006;80:84553. [DOI] [PubMed] [PMC]
    Jubair L, Fallaha S, McMillan NAJ. Systemic Delivery of CRISPR/Cas9 Targeting HPV Oncogenes Is Effective at Eliminating Established Tumors. Mol Ther. 2019;27:20919. [DOI] [PubMed] [PMC]
    Chen AC, Waterboer T, Keleher A, Morrison B, Jindal S, McMillan D, et al. Human papillomavirus in benign prostatic hyperplasia and prostatic adenocarcinoma patients. Pathol Oncol Res. 2011;17:6137. [DOI] [PubMed]
    Antonsson A, Spurr TP, Chen AC, Francis GD, McMillan NAJ, Saunders NA, et al. High prevalence of human papillomaviruses in fresh frozen breast cancer samples. J Med Virol. 2011;83:215763. [DOI] [PubMed]
    Malki MI, Gupta I, Fernandes Q, Aboulkassim T, Yasmeen A, Vranic S, et al. Co-presence of Epstein-Barr virus and high-risk human papillomaviruses in Syrian colorectal cancer samples. Hum Vaccin Immunother. 2020;16:24037. [DOI] [PubMed] [PMC]
    Cyprian FS, Al-Farsi HF, Vranic S, Akhtar S, Moustafa AA. Epstein-Barr Virus and Human Papillomaviruses Interactions and Their Roles in the Initiation of Epithelial-Mesenchymal Transition and Cancer Progression. Front Oncol. 2018;8:111. [DOI] [PubMed] [PMC]