Study of association and molecular analysis of human papillomavirus in breast cancer of Indian patients: Clinical and prognostic implication

Objectives Human papillomavirus (HPV) causes tumors primarily Cervical cancer. Recently, inconsistent reports came up in Breast cancer (BC) too. In India, despite treatment 70,218 BC patients die each year. So, we explored the association of HPV, if any, with BC prognosis in Indian pre-therapeutic (PT) and Neo-adjuvant chemotherapy (NACT) patients with subsequent analysis of HPV profile. Methods HPV prevalence was checked and analysis of physical status, copy number, genome variation, promoter methylation and expression (mRNA and protein) of the prevalent subtype was done. Results High prevalence of HPV was observed in both PT (64.0%) and NACT (71.0%) cases with significant association with younger (20-45 yrs) PT patients. Interestingly, HPV infection was significantly increased from adjacent normal breast (9.5%, 2/21), fibro adenomas (30%, 3/10) to tumors (64.8%, 203/313) samples. In both PT and NACT cases, HPV16 was the most prevalent subtype (69.0%) followed by HPV18 and HPV33. Survival analysis illustrated hrHPV infected PT patients had worst prognosis. So, detailed analysis of HPV16 profile was done which showed Europian-G350 as the most frequent HPV16 variant along with high rate of integration. Moreover, low copy number and hyper-methylation of P97 early promoter were concordant with low HPV16 E6 and E7 mRNA and protein expression. Notably, four novel variations (KT020838, KT020840, KT020841 and KT020839) in the LCR region and two (KT020836 and KT020837) in the E6 region were identified for the first time along with two novel E6^E7*I (KU199314) and E6^E7*II (KU199315) fusion transcript variants. Conclusion Thus, significant association of hrHPV with prognosis of Indian BC patients led to additional investigation of HPV16 profile. Outcomes indicated a plausible role of HPV in Indian BC patients.