Mapping of HPV transcripts in four human cervical lesions using RNAseq suggests quantitative rearrangements during carcinogenic progression

Two classes of Human papillomaviruses (HPV) infect the anogenital track: high risk viruses that are associated with risk of cervical cancer and low risk types that drive development of benign lesions, such as condylomas. In the present study, we established quantitative transcriptional maps of the viral genome in clinical lesions associated with high risk HPV16 or low risk HPV6b. Marked qualitative and quantitative changes in the HPV16 transcriptome were associated with progression from low to high grade lesions. Specific transcripts encoding essential regulatory proteins such as E7, E2, E1^E4 and E5 were identified. We also identified intrinsic differences between the HPV6b-associated condyloma transcript map and that of the HPV16-associated low grade CIN specifically regarding promoter usage. Characterization and quantification of HPV transcripts in patient samples thus establish the impact of viral transcriptional regulation on the status of HPV-associated lesions and may therefore help in defining new biologically-relevant prognosis markers.