Modulation of microRNA-mRNA target pairs by human papillomavirus 16 oncoproteins

The E6 and E7 proteins are the major oncogenic drivers encoded by high-risk human papillomaviruses (HPVs). While many aspects of the transforming activities of these proteins have been extensively studied, there are fewer studies that have investigated how HPV E6/E7 expression affects the expression of cellular noncoding RNAs. The goal of our study was to investigate HPV16 E6/E7 modulation of cellular microRNA (miR) levels and to determine the potential consequences for cellular gene expression. We performed deep sequencing of small and large cellular RNAs in primary undifferentiated cultures of human foreskin keratinocytes (HFKs) with stable expression of HPV16 E6/E7 or a control vector. After integration of the two data sets, we identified 51 differentially expressed cellular miRs associated with the modulation of 1,456 potential target mRNAs in HPV16 E6/E7-expressing HFKs. We discovered that the degree of differential miR expression in HFKs expressing HPV16 E6/E7 was not necessarily predictive of the number of corresponding mRNA targets or the potential impact on gene expression. Additional analyses of the identified miR-mRNA pairs suggest modulation of specific biological activities and biochemical pathways. Overall, our study supports the model that perturbation of cellular miR expression by HPV16 E6/E7 importantly contributes to the rewiring of cellular regulatory circuits by the high-risk HPV E6 and E7 proteins that contribute to oncogenic transformation. IMPORTANCE High-risk human papillomaviruses (HPVs) are the causative agents of almost all cervical cancers and many other cancers, including anal, vaginal, vulvar, penile, and oropharyngeal cancers. Despite the availability of efficacious HPV vaccines, it is critical to determine how HPVs cause cancer, as many people remain unvaccinated and the vaccine does not prevent cancer development in individuals who are already infected. Two HPV proteins, E6 and E7, are the major drivers of cancer development, and much remains to be learned about how the expression of these viral proteins reprograms infected cells, ultimately resulting in cancer development. Small, noncoding human RNAs, termed microRNAs (miRs), regulate gene expression and have been implicated in almost all human cancers, including HPVassociated cancers. Our study provides a comprehensive analysis of how E6 and E7 alter the expression of human miRs and how this potentially impacts cellular gene expression, which may contribute to cancer development.