RNA sequencing analysis identifies novel spliced transcripts but does not indicate quantitative or qualitative changes of viral transcripts during progression of cottontail rabbit papillomavirus-induced tumours

Persistent infections with high-risk human papillomaviruses (HPVs) can result in the development of cancer of the cervix uteri and other malignancies. The underlying molecular mechanisms leading to the progression of HPV-induced lesions are, however, not well understood. Cottontail rabbit papillomavirus (CRPV) induces papillomas in domestic rabbits which progress at a very high rate to cancer. Using this model, we compared the transcriptional patterns of CRPV in papillomas and carcinomas by RNA sequencing (RNA-seq). The most abundant transcripts can encode E7, short E6 and E1^E4, followed by full-length E6, E2, E1 and E9^E2C. In addition, we identified two rare, novel splice junctions 7810/3714 and 1751/3065 in both papillomas and carcinomas which have been described for other papillomaviruses. Neither RNA-seq nor quantitative real-time PCR-based assays identified qualitative or quantitative changes of viral transcription between papillomas and carcinomas. In summary, our analyses confirmed that papillomaviruses have highly similar transcriptional patterns, but they do not suggest that changes in these patterns contribute to the progression of CRPV-induced tumours.