- Research Article
- Open Access
Inconclusive role of human papillomavirus infection in breast cancer
Infectious Agents and Cancer volume 10, Article number: 36 (2015)
Epidemiological studies have examined the association between human papillomavirus (HPV) and breast cancer, but the findings are inconclusive. This study aimed to detect the prevalence of HPV in breast cancer tissue in patients from northeastern China and define the association between HPV and breast cancer using meta-analysis.
Polymerase chain reaction (PCR) was used to test cutaneous or mucosal HPV DNA sequence in 77 breast cancer samples and 77 corresponding adjacent normal tissues. The prevalence of HPV in breast cancer was estimated by pooling data from 38 studies. A meta-analysis of 16 case–control studies was conducted to investigate the association between HPV and breast cancer.
We did not find HPV DNA sequence in any of the 154 tissue specimens we tested. However, our meta-analysis revealed a HPV prevalence of 30.30 % (95 % confidence interval [CI] = 22.30–38.40 %) among breast cancer cases; most of these involved high-risk HPV types (35.50 %, 95 % CI = 25.00–46.10 %). HPV prevalence in breast cancer varied by geographic region, publication period, and PCR detection method. An increased risk of breast cancer was observed in association with exposure to HPV (odds ratio [OR] = 3.24, 95 % CI = 1.59–6.57), which was influenced by geographic region, HPV DNA source, PCR primer used, and publication period.
HPV, especially high-risk HPV types, may be associated with an increased risk of breast cancer, and this association varies dramatically among geographic regions.
Breast cancer is the leading cause of cancer death in females worldwide, and its incidence and mortality rates have been rising in many Asian countries . Although genetic susceptibility represents a risk factor for breast cancer, the etiologic factors that determine breast cancer risk have not been completely defined. Recently, increasing evidence has indicated that some viruses, especially human papillomavirus (HPV), may be involved in the pathogenesis of breast cancer [2–4].
It is well known that infection with specific types of HPV can cause cervical cancer. Additional evidence indicates a role for these high-risk HPV types in anal cancer, in up to 50 % of other anogenital cancers, and in 25–30 % of head and neck cancers . The association of HPV with breast cancer was indicated when Band et al. reported that HPV could immortalize normal human mammary epithelial cells and reduce their growth factor requirements . Later, HPV-16 DNA was detected in 29.4 % of breast carcinomas by polymerase chain reaction (PCR) . Up to now the presence of HPV DNA has been repeatedly detected in 0 to 86 % of breast cancer [2–4, 8–40]; this large range indicates an inconsistent association of HPV with breast cancer. Considering that multiple primers are used to detect HPV DNA in breast cancers, we hypothesized that this inconsistency may partially result from the heterogeneity of primers used in HPV detection. Although both type-specific primers and general primers have been used in previous studies, most of these were designed for mucosal HPV types [10, 11, 13, 16, 20, 21]. Thus, cutaneous HPV types might escape detection.
The aim of our study was to detect any cutaneous or mucosal HPV types in breast cancer using two set of primers: MY09/11 , which targets mucosal HPV types, and FAP 59/64 , which targets cutaneous HPV types. Given the controversy among reports regarding the relationship between HPV and breast cancer, we systematically reviewed published studies of HPV and breast cancer to investigate the global prevalence and type distribution of HPV in breast cancer, to determine potential factors related to HPV prevalence, and to examine the evidence regarding the association of HPV infection and breast cancer risk.
Of the 77 patients (age range: 39–74 years; mean: 49 years; SE = 9.43) selected for our molecular study, 61 had been diagnosed with invasive ductal carcinomas. The clinical characteristics of the breast cancer samples were shown in Table 1. None of the patients reported a family history of any other cancer. DNA isolated from all of these cases underwent detection of HPV infection by PCR. All samples were positive for β-globin gene amplification (Fig. 1a), indicating that DNA was available for molecular analysis. However, no HPV DNA was detected using either mucosal consensus primers (MY09/11) or cutaneous consensus primers (FAP59/64) in any of the 77 breast tumor tissues or 77 normal adjacent tissues (Fig. 1b and c). Weak bands between 400 and 500 bp appeared in two breast cancer samples, but the results of three extra PCR runs of the first-cycle products or DNA from these specimens were negative (data not shown).
The overall and type-specific prevalence of HPV in breast cancer and factors impacting the prevalence of HPV in breast cancer
A total of 326 published records were retrieved using the key words mentioned previously (Fig. 2). Among those publications, 38 studies [2–4, 7–40] were pooled to calculate the HPV prevalence in breast cancer. The selected studies comprised a total of 2569 breast cancer cases, and the majority came from Asia (51.96 %) and South America (20.32 %). The number of cases investigated in each study varied from 17 to 228, and the prevalence of HPV ranged from 0 to 86 %.
The pooled prevalence of HPV in breast cancer across the 38 studies was 30.30 % (95 % CI = 22.30–38.40 %). Stratified by region, Oceania yielded the highest HPV prevalence in breast cancer, with 44.30 % (95 % CI = 33.50–55.00 %), followed by Asia with 35.70 % (95 % CI = 20.60–50.80 %), Europe with 32.60 % (95 % CI = 6.50–58.70 %), South America with 14.60 % (95 % CI = 5.80–23.50 %), and North America with 10.70 % (95 % CI = 4.40–17.00 %). When the unconditional logit regression model was introduced to compare HPV prevalence in breast cancer among different geographic regions, the result remained consistent: Oceanian breast cancer patients exhibited the highest HPV prevalence, and the priority of this prevalence was statistically significant (P < 0.001) (Table 2). HPV prevalence was found to be higher when HPV DNA was extracted from fresh tissues (37.30 %, 95 % CI = 20.10–54.50 %) than when HPV DNA was extracted from paraffin-embedded tissues (27.20 %, 95 % CI = 17.90–36.40 %), but this priority of prevalence was statistically insignificant (p = 0.08) (Table 2).
A cumulative meta-analysis was conducted to investigate the impact of publication period on HPV prevalence in breast cancer, and we found that the HPV detection rate between 2007 and 2012 (22.60 %, 95 % CI = 14.40–30.70 %) was dramatically lower than that between 1992 and 2006 (42.20 %, 95 % CI = 25.60–58.80 %) (p < 0.001) (Table 2).
PCR was the dominant method used to detect HPV in breast cancer; it was used in 31.30 % (95 % CI = 22.60–39.90 %) of the studies we reviewed, while HPV prevalence was 5.30 % (95 % CI = 2.40–8.30 %) using other methods, such as in situ hybridization and Southern blotting (Table 3). Multiple primers, including broad-spectrum PCR primers, type-specific PCR primers, and a combination of these, were used in PCR-based methods of detecting HPV. The HPV prevalence rate was 30.00 % (95 % CI = 18.30–41.70 %) when broad-spectrum PCR primers were used; it was 31.90 % (95 % CI = 13.00–50.80 %) when type-specific PCR primers were used; and it was 35.30 % (95 % CI = 11.50 –59.00 %) when a combination of these two types of primers was used. Because most studies dealt with invasive ductal cancer as well as non-invasive cancer, the detection rates of HPV in invasive ductal cancer were pooled, yielding HPV prevalence of 32.90 % (95 % CI = 21.50–44.20 %) compared with non-invasive cancer, which showed a prevalence rate of 22.00 % (95 % CI = 12.40–31.50 %) (Table 3).
Although multiple HPV types were determined in breast cancer cases across the 38 studies in our meta-analysis, the five most common HPV types, in decreasing order of prevalence, were the following: HPV-16 (30.70 %; 95 % CI = 20.50–41.00 %), HPV-33 (25.60 %; 95 % CI = 7.10–44.00 %), HPV-11 (17.20 %; 95 % CI = 1.00–33.30 %), HPV-18 (9.90 %; 95 % CI = 6.20–13.70 %), and HPV-6 (6.30 %; 95 % CI = 2.00–10.50 %) (Table 4). Classified by HPV oncogenic features, the prevalence of high-risk HPV types was 35.50 % (95 % CI = 25.00–46.10 %), much higher than the prevalence of low-risk HPV types, which was 11.70 % (95 % CI = 5.80–17.70 %).
The association between HPV infection and breast cancer risk
A total of 16 case–control studies were included in our evaluation of the association between HPV infection and breast cancer risk (Fig. 1). Due to existing heterogeneity (p = 0.000 for heterogeneity test, I 2 = 63.9 %), a random model was chosen to calculate the combined OR and its 95 % CI across the 16 studies. When these 16 studies were pooled in the meta-analysis, results showed that HPV infection was associated with higher risk for breast cancer (OR =3.24, 95 % CI = 1.59–6.57, p < 0.001) (Fig. 3). Egger’s and Begg’s tests were performed to assess the publication bias; the results of both indicated no statistical significance (p = 0.16 and 0.45, respectively).
Considering that substantial heterogeneity was observed among the 16 studies for HPV infection and breast cancer, the prevalence of HPV infection in breast cancer was further evaluated by subgroup analysis. The analyzed subgroups were defined according to the main features of pooled studies: geographical region, DNA source, detection method (i.e., PCR primer), and publication period (Table 5). In the subgroup of geographical regions, between-study heterogeneity was absent (p = 0.611 for heterogeneity test, I 2 = 0.0 %) in Oceania and Europe, but pooled OR remained the same (OR = 3.16, 95 % CI = 1.27–7.90, p < 0.001). When ORs were grouped by DNA source, however, significantly higher risk for breast cancer was revealed among fresh-tissue methods (OR = 7.88, 95% CI = 3.99–15.60, p < 0.001) without between-study heterogeneity. Similarly, breast cancer risk was greater with HPV infection in studies published between 1992 and 2006 (OR = 4.12, 95% CI = 2.42–7.02, p < 0.001) compared with studies published between 2007 and 2012. HPV infection detected using broad-spectrum PCR primers was associated with higher risk for breast cancer (OR = 5.67, 95 % CI = 3.40–9.45, p < 0.001), while HPV infection detected using combined PCR primers showed the opposite result (OR = 0.68, 95 % CI = 0.32–1.45, p < 0.001).
The role of HPV in breast cancer has been controversial; thus, a PCR-based strategy as well as a meta-analysis were conducted to investigate the association between HPV and breast cancer. In present study, 77 breast cancer samples and 77 adjacent normal tissues all tested negative for HPV. However, in our meta-analysis of 38 studies from 20 countries published since 1992, we found the overall HPV prevalence in women with breast cancer to be 30.30 % (95 % CI = 22.30–38.40 %). Furthermore, the prevalence estimate varied by geographic region, publication period, and PCR detection method. Our meta-analysis of 16 pooled case–control studies revealed that HPV infection is associated with breast cancer risk.
In our molecular study, both cutaneous and mucosal primers were applied to detect HPV DNA in breast tumors, but no specimen was positive for HPV. Other studies [9, 37] from China in which HPV DNA was detected in breast tumors have reported similar findings, suggesting that the HPV infection rate in breast cancer in Chinese patients is low and that HPV infection rate varies by geographic region. To support this idea, 38 studies [2–4, 7–40] from different regions were pooled to calculate the prevalence of HPV infection in breast cancer and to explore the factors impacting it.
The prevalence of HPV infection in breast cancer ranged from 10.70 % in North America to 44.30 % in Oceania. The geographic differences of HPV infection rates may be related to the distinctly different cervical cancer burdens across different continents . In 1999, Hennig et al. reported that HPV-associated cervical neoplasia might be the original site of HPV infection from which the virus could be transported to the breast . We found that HPV prevalence in frozen specimens was higher than that in paraffin-embedded tissues, but the difference was statistically insignificant (P = 0.08). Interestingly, the HPV detection rate of studies published between 1992 and 2006 was remarkably higher than that of studies published between 2007 and 2012. This can be explained by the fact that the latter generally included more breast cancer cases and that the specificity of detection advanced over time, thus yielding a conservative result.
In present study, HPV prevalence in 77 breast cancer samples and 77 adjacent normal tissues was zero, which is lower than that estimated by the meta-analysis included in this report. However, this is consistent with the results showed by stratified analysis, i.e. , HPV infection rate varies by geographic region and publication period (shown in Table 2). The geographic differences of HPV infection rates may be related to the distinctly different cervical cancer burdens across different continents, since cervix is proposed to be the original site of HPV infection from which HPV could be transported to the breast. In addition, due to the limitations of statistical software, 11 studies (references 8–12, 19–21, 27, 29, 30), as well as the present PCR-based study, showing no HPV prevalence in breast cancer, were not included in the pool used to estimate overall HPV prevalence. Thus, our pooled results showed a relatively strong combined HPV prevalence.
PCR emerged as the most effective method of detecting HPV in breast cancer, and combining general primers and specific primers revealed the highest HPV prevalence, due to the diversity of HPV types that co-occur with breast cancer. In our meta-analysis, HPV frequency in invasive carcinomas was higher than that in non-invasive carcinomas, but the difference was not statistically significant, suggesting that HPV infection may not be directly associated with the progression of breast cancer. Although multiple HPV types have been detected in breast cancer, high-risk HPV types are more frequently detected than low-risk HPV types. These results imply that the use of primers targeting high-risk HPV will increase the detection rate of HPV in breast cancer tissue.
Although the pooled estimate from 16 randomized controlled trials revealed that HPV infection resulted in an increased risk of breast cancer, the association was not strong enough to prove that HPV plays a causal role in breast cancer development. To identify the definitive role of HPV in breast cancer development, further studies are required to investigate the temporal association between the virus and breast cancer, the integration in the host genome, the expression of viral oncoproteins, the traits of HPV-positive breast cancer, and the route of breast infection.
Our meta-analysis has several limitations. First, due to the limitations of statistical software, 11 studies [8–12, 19–21, 27, 29, 30], as well as the present PCR-based study, showing HPV prevalence in breast cancer was zero, were not included in the pool used to estimate overall HPV prevalence; thus, our pooled results showed a relatively strong combined HPV prevalence. Interestingly, of 885 cases from the 11 studies that detected no HPV, 306 [9, 12, 20] cases (34.5 %) were from Asia and in them PCR was the only method used to detect HPV DNA. These results reflect those of the current PCR-based study. Second, substantial heterogeneity was observed when the pooled HPV prevalence was calculated across different continents. But unconditional logit regression revealed variables, including geographic region and publication period, that can only explain part of the existing heterogeneity. Third, some evidence supports the idea that some clinical features, including age and estrogen receptors, are linked to the HPV detection rate in breast cancer , and we cannot rule out the possibility of such an effect. Detailed information about these clinical characteristics were unavailable in the studies used in our meta-analysis.
Although two meta-analyses [45, 46] have been published regarding the correlation between HPV infection and the risk of breast cancer, our meta-analysis included a greater number of publications as well as more detailed information on HPV prevalence in breast cancer. Thus, our study has yielded more valid results.
HPV infection, especially high-risk HPV types, may be associated with an increased risk of breast cancer, and the association varies dramatically among geographic regions.
A total of 77 breast cancer specimens and 77 corresponding paraneoplastic breast tissues were collected at the time of operation and were frozen immediately at −80 °C. These samples were obtained from patients diagnosed between February 2011 and May 2012 at the First Affiliated Hospital of Harbin Medical University (Harbin, China). All of the patients were women living in northeastern China, and none had a history of cervical carcinoma. Informed consent according to the criteria of Harbin Medical University was obtained from all patients, and this study was approved by the Institutional Research Board of Harbin Medical University (No. HMUIRB20120002).
DNA extraction and PCR amplification
DNA was extracted from fresh frozen breast tissues as previously described, using SDS lysis/proteinase-K digestion, phenol/chloroform extraction, and ethanol precipitation. The DNA quality was confirmed by the amplification of a 268-bp fragment of human β-globin gene using and primers G073 (5′-GAAGAGCCAAG GACAGGTAC-3′) and G074 (5′-CAACTTCATCCACGTTCACC-3′). To avoid contamination between samples, the cutting blade was changed after the sectioning of each sample.
DNA samples were then screened for the presence of HPV by PCR using two sets of primers: MY09/11 and FAP59/64. Consensus primers MY09 (5′-CGTCCMARRGGAWACTGATC-3′) and MY11 (5′-GCMCAGGGWCATAAY AATGG-3′) amplify a 450-bp fragment in the highly conserved region in L1 of mucosal HPV types (annealing temperature 50 °C). Degenerated primers FAP59 (5′-TAACWGTIGGICAYCCWTATT-3′) and FAP64 (5′-CCWATATCWVHCATIT CICCATC-3′) can detect a 478-bp fragment in the L1 region of both mucosal and cutaneous HPV types (annealing temperature 55 °C). PCR was performed in a volume of 20 μL, containing 100 ng of extracted DNA, 1.67 mM MgCl2, 200 μM dNTPs, 0.4 μM of each primer and 0.5 U Taq DNA polymerase (TaKaRa, Japan), and reaction buffer (10 × PCR buffer, TaKaRa, Japan), according to the following conditions: 5 min at 94 °C, then 40 cycles of 5 s at 98 °C, 30s at 50 or 55 °C, and 60 s at 72 °C, followed by a 5-min extension at 72 °C in an automated thermocycler (Eppendorf, Germany). The full genomes of HPV-16, −11, and −8 were cloned in the pBS322 plasmid (donated by Michelle A. Ozbun), and HPV-18-positive cervical cancer specimens were used as positive controls; sterile water was used as a negative control. The PCR products were separated on 1.5 % agarose gel (Blowest, Spain) and visualized with ethidium bromide staining after electrophoresis.
Search strategy and selection criteria for HPV and breast cancer
Eligible articles published from January 1989 to June 2013 were identified by an electronic PubMed search using the MeSH terms “human,” “papillomavirus,” and “breast carcinoma.” Additional articles from the reference lists of retrieved articles were also screened. The inclusion criteria of the studies were as follows: (1) detection of HPV in female breast cancer; (2) full text in English; (3) independent of other studies. The experimental results from the present molecular study were included in the meta-analysis.
Data extraction for HPV and breast cancer
Two authors (Jinyuan Li and Lanlan Wei) independently extracted data and reached a consensus on all issues. Disagreements were discussed and resolved according to the inclusion criteria and consensus. The following data were extracted from each report: the first author, year of publication, country of origin, ethnicity, sample size (cases and controls), HPV prevalence, DNA source, detection method, PCR primers used, histological subtypes tested, and HPV types detected.
An epidemiological review of the overall and type-specific HPV prevalence in breast cancer cases was conducted. Four high-risk HPV types were identified (HPV-16, −18, −31, −33), and two low-risk HPV types were identified (HPV-6 and −11). A meta-analysis was performed to explore the association between HPV infection and breast cancer risk in the form of a case (breast cancer tissues) versus control (normal breast cancer tissues, breast tumor adjacent tissues, or benign breast lesions) comparison. An unconditional logistic regression model was used to compare HPV prevalence according to the influential parameters: HPV DNA source, PCR primers used in detection, and publication period.
A fixed-effect model (Mantel-Haenszel method) and a random-effect model (DerSimonian and Laird method) were used to pool the case–control data. These two models provide similar results when between-study heterogeneity is absent; otherwise, the random-effect model is more appropriate. Between-study heterogeneity was tested using the χ2-based Q test, and heterogeneity was considered significant if p < 0.05. Subgroup analyses were further performed to explore the source of the existing heterogeneity. Publication bias was evaluated using the linear regression asymmetry tests designed by Egger et al. . Analyses were carried out using Stata version 11 software (StataCorp, College Station, TX).
Jemal A, Center MM, DeSantis C, Ward EM. Global patterns of cancer incidence and mortality rates and trends. Cancer Epidemiol Biomark Prev. 2010;19:1893–907.
Damin AP, Karam R, Zettler CG, Caleffi M, Alexandre CO. Evidence for an association of human papillomavirus and breast carcinomas. Breast Cancer Res Treat. 2004;84:131–7.
De Villiers EM, Sandstrom RE, zur Hausen H, Buck CE. Presence of papillomavirus sequences in condylomatous lesions of the mamillae and in invasive carcinoma of the breast. Breast Cancer Res. 2005;7:R1–11.
Heng B, Glenn WK, Ye Y, Tran B, Delprado W, Lutze-Mann L, et al. Human papilloma virus is associated with breast cancer. Br J Cancer. 2009;101:1345–50.
zur Hausen H. Papillomaviruses in the causation of human cancers - a brief historical account. Virology. 2009;384:260–5.
Vimla Band DZ, Kulesa V, Sager R. Human papilloma virus DNAs immortalize normal human mammary epithelial cells and reduce their growth factor requirements. Proc Natl Acad Sci U S A. 1990;87:463–7.
Di Lonardo A. AVaMLM: Human papillomavirus in breast cancer Breast Cancer Research and Treatment. 1992.
Czerwenka K, Heuss F, Hosmann JW, Manavi M, Lu Y, Jelincic D, et al. Human papilloma virus DNA: a factor in the pathogenesis of mammary Paget’s disease? Breast Cancer Res Treat. 1996;41:51–7.
Chang P, Wang T, Yao Q, Lv Y, Zhang J, Guo W, et al. Absence of human papillomavirus in patients with breast cancer in north-west China. Med Oncol. 2012;29:521–5.
Bratthauer GL, Tavassoli FA, O’Leary TJ. Etiology of breast carcinoma: no apparent role for papillomavirus types 6/11/16/18. Pathol Res Pract. 1992;188:384–6.
De Cremoux P, Thioux M, Lebigot I, Sigal-Zafrani B, Salmon R, Sastre-Garau X. No evidence of human papillomavirus DNA sequences in invasive breast carcinoma. Breast Cancer Res Treat. 2008;109:55–8.
Gopalkrishna V, Singh UR, Sodhani P, Sharma JK, Hedau ST, Mandal AK, et al. Absence of human papillomavirus DNA in breast cancer as revealed by polymerase chain reaction. Breast Cancer Res Treat. 1996;39:197–202.
Aguayo F, Khan N, Koriyama C, González C, Ampuero S, Padilla O, et al. Human papillomavirus and Epstein-Barr virus infections in breast cancer from chile. Infect Agent Cancer. 2011;6:7.
Akil N, Yasmeen A, Kassab A, Ghabreau L, Darnel AD, Al Moustafa AE. 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:404–7.
Antonsson A, Spurr TP, Chen AC, Francis GD, McMillan NA, Saunders NA, et al. High prevalence of human papillomaviruses in fresh frozen breast cancer samples. J Med Virol. 2011;83:2157–63.
Baltzell K, Buehring GC, Krishnamurthy S, Kuerer H, Shen HM, Sison JD. Limited evidence of human papillomavirus in [corrected] breast tissue using molecular in situ methods. Cancer. 2012;118:1212–20.
Choi YL, Cho EY, Kim JH, Nam SJ, Oh YL, Song SY, et al. Detection of human papillomavirus DNA by DNA chip in breast carcinomas of Korean women. Tumour Biol. 2007;28:327–32.
de Leon DC, Montiel DP, Nemcova J, Mykyskova I, Turcios E, Villavicencio V, et al. Human papillomavirus (HPV) in breast tumors: prevalence in a group of Mexican patients. BMC Cancer. 2009;9:26.
Herrera-Romano L, Fernandez-Tamayo N, Gomez-Conde E, Reyes-Cardoso JM, Ortiz-Gutierrez F, Ceballos G, et al. Absence of human papillomavirus sequences in epithelial breast cancer in a Mexican female population. Med Oncol. 2012;29:1515–7.
Hedau S, Kumar U, Hussain S, Shukla S, Pande S, Jain N, et al. Breast cancer and human papillomavirus infection: no evidence of HPV etiology of breast cancer in Indian women. BMC Cancer. 2011;11:27.
Hachana M, Ziadi S, Amara K, Toumi I, Korbi S, Trimeche M. No evidence of human papillomavirus DNA in breast carcinoma in Tunisian patients. Breast. 2010;19:541–4.
Duo D, Ghimenti C, Migliora P, Pavanelli MC, Mastracci L, Angeli G. Identification and characterization of human papillomavirus DNA sequences in Italian breast cancer patients by PCR and line probe assay reverse hybridization. Mol Med Report. 2008;1:673–7.
Frega A, Lorenzon L, Bononi M, De Cesare A, Ciardi A, Lombardi D, et al. Evaluation of E6 and E7 mRNA expression in HPV DNA positive breast cancer. Eur J Gynaecol Oncol. 2012;33:164–7.
Gumus M, Yumuk PF, Salepci T, Aliustaoglu M, Dane F, Ekenel M, et al. HPV DNA frequency and subset analysis in human breast cancer patients’ normal and tumoral tissue samples. J Exp Clin Cancer Res. 2006;25:515–21.
He Q, Zhang SQ, Chu YL, Jia XL, Wang XL. The correlations between HPV16 infection and expressions of c-erbB-2 and bcl-2 in breast carcinoma. Mol Biol Rep. 2009;36:807–12.
Herrera-Goepfert R, Khan NA, Koriyama C, Akiba S, Perez-Sanchez VM. High-risk human papillomavirus in mammary gland carcinomas and non-neoplastic tissues of Mexican women: no evidence supporting a cause and effect relationship. Breast. 2011;20:184–9.
Lindel K, Forster A, Altermatt HJ, Greiner R, Gruber G. Breast cancer and human papillomavirus (HPV) infection: no evidence of a viral etiology in a group of Swiss women. Breast. 2007;16:172–7.
Ong K, Koay ES, Putti TC. Detection of cutaneous HPV types 4 and 24 DNA sequences in breast carcinoma in Singaporean women of Asian ancestry. Pathology. 2009;41:436–42.
Silva Jr RG, da Silva BB. No evidence for an association of human papillomavirus and breast carcinoma. Breast Cancer Res Treat. 2011;125:261–4.
Wrede D, Luqmani YA, Coombes RC, Vousden KH. Absence of HPV 16 and 18 DNA in breast cancer. Br J Cancer. 1992;65:891–4.
Kan CY, Iacopetta BJ, Lawson JS, Whitaker NJ. Identification of human papillomavirus DNA gene sequences in human breast cancer. Br J Cancer. 2005;93:946–8.
Khan NA, Castillo A, Koriyama C, et al. Human papillomavirus detected in female breast carcinomas in Japan. Br J Cancer. 2008;99:408–14.
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:727–31.
Li T, Lu ZM, Guo M, Wu QJ, Chen KN, Xing HP, et al. p53 codon 72 polymorphism (C/G) and the risk of human papillomavirus-associated carcinomas in China. Cancer. 2002;95:2571–6.
Liu Y, Klimberg VS, Andrews NR, Hicks CR, Peng H, Chiriva-Internati M, et al. Human papillomavirus DNA is present in a subset of unselected breast cancers. J Hum Virol. 2001;4:329–34.
Mendizabal-Ruiz AP, Morales JA, Ramirez-Jirano LJ, Padilla-Rosas M, Moran-Moguel MC, Montoya-Fuentes H. Low frequency of human papillomavirus DNA in breast cancer tissue. Breast Cancer Res Treat. 2009;114:189–94.
Mou X, Chen L, Liu F, Shen Y, Wang H, Li Y, et al. Low prevalence of human papillomavirus (HPV) in Chinese patients with breast cancer. J Int Med Res. 2011;39:1636–44.
Sigaroodi A, Nadji SA, Naghshvar F, Nategh R, Emami H, Velayati AA. Human papillomavirus is associated with breast cancer in the north part of iran. ScientificWorldJournal. 2012;2012:837191.
Tsai JH, Tsai CH, Cheng MH, Lin SJ, Xu FL, Yang CC. Association of viral factors with non-familial breast cancer in Taiwan by comparison with non-cancerous, fibroadenoma, and thyroid tumor tissues. J Med Virol. 2005;75:276–81.
Widschwendter A, Brunhuber T, Wiedemair A, Mueller-Holzner E, Marth C. Detection of human papillomavirus DNA in breast cancer of patients with cervical cancer history. J Clin Virol. 2004;31:292–7.
Bauer HM, Ting Y, Greer CE, Chambers JC, Tashiro CJ, Chimera J, et al. Genital human papillomavirus infection in female university students as determined by a PCR-based method. JAMA. 1991;265:472–7.
Forslund O, Antonsson A, Nordin P, Stenquist B, Hansson BG. A broad range of human papillomavirus types detected with a general PCR method suitable for analysis of cutaneous tumours and normal skin. J Gen Virol. 1999;80(Pt 9):2437–43.
WHO/ICO Information Centre on HPV and Cervical Cancer (HPV Information Centre): . Human Papillomavirus and Related Cancers in China. Summary Report 2014. Available at www.hpvcentre.net/summaryreport.php
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:121–35.
Li N, Bi X, Zhang Y, Zhao P, Zheng T, Dai M. Human papillomavirus infection and sporadic breast carcinoma risk: a meta-analysis. Breast Cancer Res Treat. 2011;126:515–20.
Simoes PW, Medeiros LR, Simoes Pires PD, Edelweiss MI, Rosa DD, Silva FR, et al. Prevalence of human papillomavirus in breast cancer: a systematic review. Int J Gynecol Cancer. 2012;22:343–7.
Egger M, Davey Smith G, Schneider M, Minder C. Bias in meta-analysis detected by a simple, graphical test. BMJ. 1997;315:629–34.
This study was supported by grants from the National Natural Science Foundation of China (No. 30901706), the Natural Science Foundation of Heilongjiang Province (No. D201065 and D201051). The authors thank Michelle A. Ozbun (Molecular Genetics and Microbiology, University of New Mexico) for providing HPV-16 and HPV-11 clones in the pBS322 plasmid.
The authors declare that they have no competing interests.
LW conceived and designed the experiments; YZ, JL, YJ, MR, BP, and MC performed the experiments; all authors read and approve the manuscript.
Yi Zhou and Jinyuan Li are Co-first authors.