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.