RNA-seq analysis in equine papillomavirus type 2-positive carcinomas identifies affected pathways and potential cancer markers as well as viral gene expression and splicing events
Equine papillomavirus type 2 (EcPV2) was discovered only recently, but it is found consistently in the context of genital squamous cell carcinomas (SCCs). Since neither cell cultures nor animal models exist, the characterization of this potential disease agent relies on the analysis of patient materials. To analyse the host and viral transcriptome in EcPV2-affected horses, genital tissue samples were collected from horses with EcPV2-positive lesions as well as from healthy EcPV2-negative horses. It was determined by RNA-seq analysis that there were 1957 differentially expressed (DE) host genes between the SCC and control samples. These genes were most abundantly related to DNA replication, cell cycle, extracellular matrix (ECM)–receptor interaction and focal adhesion. By comparison to other cancer studies, MMP1 and IL8 appeared to be potential marker genes for the development of SCCs. Analysis of the viral reads revealed the transcriptional activity of EcPV2 in all SCC samples. While few reads mapped to the structural viral genes, the majority of reads mapped to the non-structural early (E) genes, in particular to E6, E7 and E2/E4. Within these reads a distinct pattern of splicing events, which are essential for the expression of different genes in PV infections, was observed. Additionally, in one sample the integration of EcPV2 DNA into the host genome was detected by DNA-seq and confirmed by PCR. In conclusion, while host MMP1 and IL8 expression and the presence of EcPV2 may be useful markers in genital SCCs, further research on EcPV2-related pathomechanisms may focus on cell cycle-related genes, the viral genes E6, E7 and E2/E4, and integration events.
ScaseT, BrandtS, KainzbauerC, SykoraS, BijmholtS et al. Equus caballus papillomavirus-2 (EcPV-2): an infectious cause for equine genital cancer?. Equine Vet J2010; 42:738–745 [View Article]
LangeCE, ToblerK, LehnerA, GrestP, WelleMM et al. EcPV2 DNA in equine papillomas and in situ and invasive squamous cell carcinomas supports papillomavirus etiology. Vet Pathol2013; 50:686–692 [View Article]
BogaertL, WillemsenA, VanderstraetenE, BrachoMA, De BaereC et al. EcPV2 DNA in equine genital squamous cell carcinomas and normal genital mucosa. Vet Microbiol2012; 158:33–41 [View Article]
KnightCG, MundayJS, RosaBV, KiupelM. Persistent, widespread papilloma formation on the penis of a horse: a novel presentation of equine papillomavirus type 2 infection. Vet Dermatol2011; 22:570–574 [View Article]
SykoraS, SamekL, SchönthalerK, PalmF, BorzacchielloG et al. EcPV-2 is transcriptionally active in equine SCC but only rarely detectable in swabs and semen from healthy horses. Vet Microbiol Elsevier B.V 2012; 158:194–198 [View Article]
LangeCE, ToblerK, AckermannM, FavrotC. Identification of two novel equine papillomavirus sequences suggests three genera in one cluster. Vet Microbiol Elsevier B.V 2011; 149:85–90 [View Article]
KnightCG, DunowskaM, MundayJS, Peters-KennedyJ, RosaBV. Comparison of the levels of Equus caballus papillomavirus type 2 (EcPV-2) DNA in equine squamous cell carcinomas and non-cancerous tissues using quantitative PCR. Vet Microbiol Elsevier B.V 2013; 166:257–262 [View Article]
NewkirkKM, HendrixDVH, AnisEA, RohrbachBW, EhrhartEJ et al. Detection of papillomavirus in equine periocular and penile squamous cell carcinoma. J Vet Diagn Invest2014; 26:131–135 [View Article]
van den TopJGB, HarkemaL, LangeC, EnsinkJM, van de LestCHA et al. Expression of p53, Ki67, EcPV2- and EcPV3 DNA, and viral genes in relation to metastasis and outcome in equine penile and preputial squamous cell carcinoma. Equine Vet J2015; 47:188–195 [View Article]
Van den TopJGB, de HeerN, KleinWR, EnsinkJM. Penile and preputial squamous cell carcinoma in the horse: a retrospective study of treatment of 77 affected horses. Equine Vet J Blackwell Publishing Ltd 2008; 40:533–537 [View Article]
MairTS, WalmsleyJP, PhillipsTJ. Surgical treatment of 45 horses affected by squamous cell carcinoma of the penis and prepuce. Equine Vet J Blackwell Publishing Ltd 2000; 32:406–410 [View Article]
SykoraS, JindraC, HoferM, SteinbornR, BrandtS. Equine papillomavirus type 2: an equine equivalent to human papillomavirus 16?. Vet J Elsevier Ltd 2017; 225:3–8 [View Article]
McBrideAA, WarburtonA. The role of integration in oncogenic progression of HPV-associated cancers. PLoS Pathog Public Library of Science 2017; 13:e1006211 [View Article]
NakajimaK, KanekoT, AizuT, NakanoH, MatsuzakiY et al. Signet-ring cutaneous squamous cell carcinoma arising on the back of the finger. Case Rep Dermatol Karger Publishers 2013; 5:215–218 [View Article]
Espinosa de los MonterosA, Aguirre-SanceledonioM, RamírezGA, CastroP, RodríguezF. Signet-ring squamous cell carcinoma in a dog. Vet Rec2003; 153:90–92 [View Article]
DemellawyDE, OnumaK, AlowamiS. Signet ring squamous cell carcinoma-the forgotten variant: case report and review of the literature. J Cutan Pathol John Wiley & Sons, Ltd 2011; 38:306–308 [View Article]
YuanZQ, NicolsonL, MarchettiB, GaultEA, CampoMS et al. Transcriptional changes induced by bovine papillomavirus type 1 in equine fibroblasts. J Virol American Society for Microbiology 2008; 82:6481–6491 [View Article]
SantegoetsLAM, van BaarsR, TerlouA, Heijmans-AntonissenC, SwagemakersSMA et al. Different DNA damage and cell cycle checkpoint control in low- and high-risk human papillomavirus infections of the vulva. Int J Cancer John Wiley & Sons 2012; 130:2874–2885 [View Article]
ClarkeB, ChettyR. Cell cycle aberrations in the pathogenesis of squamous cell carcinoma of the uterine cervix. Gynecol Oncol2001; 82:238–246 [View Article]
MoodyCA, LaiminsLA. Human papillomavirus oncoproteins: pathways to transformation. Nat Rev Cancer Nature Publishing Group 2010; 10:550–560 [View Article]
YimEK, ParkJS. The role of HPV E6 and E7 oncoproteins in HPV-associated cervical carcinogenesis. Cancer Res Treat Korean Cancer Association 2005; 37:319–324 [View Article]
BrennaSMF, SyrjänenKJ. Regulation of Cell Cycles Is of Key Importance in Human Papillomavirus (Hpv)-Associated Cervical Carcinogenesis. Sao Paulo Medical Journal121, 5th ed. Associação Paulista de Medicina; 2003 pp 128–132
GriffinH, SonejiY, Van BaarsR, AroraR, JenkinsD et al. Stratification of HPV-induced cervical pathology using the virally encoded molecular marker E4 in combination with p16 or MCM. Mod Pathol2015; 28:977–993 [View Article]
von Knebel DoeberitzM, ReuschenbachM, SchmidtD, BergeronC. Biomarkers for cervical cancer screening: the role of p16INK4a to highlight transforming HPV infections. Expert Rev Proteomics Taylor & Francis 2012; 9:149–163 [View Article]
LuP, WeaverVM, WerbZ. The extracellular matrix: a dynamic niche in cancer progression. J Cell Biol Rockefeller University Press 2012; 196:395–406 [View Article]
GérardC, GoldbeterA. The balance between cell cycle arrest and cell proliferation: control by the extracellular matrix and by contact inhibition. Interface Focus2014; 4:20130075–5 [View Article]
Ben KhalifaY, TeissierS, TanMKM, PhanQT, DaynacM et al. The human papillomavirus E6 oncogene represses a cell adhesion pathway and disrupts focal adhesion through degradation of TAp63β upon transformation. Lambert P, editor PLoS Pathog Public Library of Science 2011; 7:e1002256 [View Article]
TongX, HowleyPM. The bovine papillomavirus E6 oncoprotein interacts with paxillin and disrupts the actin cytoskeleton. Proc Natl Acad Sci USA National Academy of Sciences 1997; 94:4412–4417 [View Article]
HollowayA, StoreyA. A conserved C-terminal sequence of high-risk cutaneous beta-human papillomavirus E6 proteins alters localization and signalling of β1-integrin to promote cell migration. J Gen Virol2014; 95:123–134 [View Article]
KanehisaM, GotoS, SatoY, KawashimaM, FurumichiM et al. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res2014; 42:D199–D205 [View Article]
AkgülB, García-EscuderoR, GhaliL, PfisterHJ, FuchsPG et al. The E7 protein of cutaneous human papillomavirus type 8 causes invasion of human keratinocytes into the dermis in organotypic cultures of skin. Cancer Res American Association for Cancer Research 2005; 65:2216–2223 [View Article]
HsiaoYP, YangJH, WuWJ, LinMH, SheuGT. E6 and E7 of human papillomavirus type 18 and UVB irradiation corporately regulate interleukin-6 and interleukin-8 expressions in basal cell carcinoma. Exp Dermatol2013; 22:672–674 [View Article]
RajkumarT, SabithaK, VijayalakshmiN, ShirleyS, BoseMV et al. Identification and validation of genes involved in cervical tumourigenesis. BMC Cancer BioMed Central 2011; 11:80 [View Article]
RegeziJA, DekkerNP, RamosDM, LiX, Macabeo-OngM et al. Proliferation and invasion factors in HIV-associated dysplastic and nondysplastic oral warts and in oral squamous cell carcinoma: an immunohistochemical and RT-PCR evaluation. Oral Surg Oral Med Oral Pathol Oral Radiol Endod2002; 94:724–731 [View Article]
RyzhakovaOS, SolovyevaNI, metalloproteinasesM. Matrix metalloproteinases (MMP)—MMP-1,-2,-9 and endogenous regulators of their activity in cervical carcinoma cell lines transformed by the E7 oncogene HPV16 and HPV18. Biochemistry (Moscow) Supplement Series B: Biomedical Chemistry2012; 6:247–253 [View Article]
WalkerJ, SmileyLC, IngramD, RomanA. Expression of human papillomavirus type 16 E7 is sufficient to significantly increase expression of angiogenic factors but is not sufficient to induce endothelial cell migration. Virology2011; 410:283–290 [View Article]
YuanZ, GaultEA, CampoMS, NasirL. Upregulation of equine matrix metalloproteinase 1 by bovine papillomavirus type 1 is through the transcription factor activator protein-1. J Gen Virol2011; 92:2608–2619 [View Article]
PruittKD, BrownGR, HiattSM, Thibaud-NissenF, AstashynA et al. RefSeq: an update on mammalian reference sequences. Nucleic Acids Res Oxford University Press 2014; 42:D756–D763 [View Article]
ChenJ, XueY, PoidingerM, LimT, ChewSH et al. Mapping of HPV transcripts in four human cervical lesions using RNAseq suggests quantitative rearrangements during carcinogenic progression. Virology Elsevier 2014; 462-463:14–24 [View Article]
LiY, WangX, NiT, WangF, LuW et al. Human papillomavirus type 58 genome variations and RNA expression in cervical lesions. J Virol2013; 87:9313–9322 [View Article]
XueXY, MajerciakV, UberoiA, KimBH, GotteD et al. The full transcription map of mouse papillomavirus type 1 (MmuPV1) in mouse wart tissues. Hearing P, editor PLoS Pathog2017; 13:e1006715–1006736 [View Article]
JohanssonC, SchwartzS. Regulation of human papillomavirus gene expression by splicing and polyadenylation. Nat Rev Microbiol Nature Publishing Group 2013; 11:239–251 [View Article]
Van DoorslaerK, TanQ, XirasagarS, BandaruS, GopalanV et al. The papillomavirus Episteme: a central resource for papillomavirus sequence data and analysis. Nucleic Acids Res Oxford University Press 2013; 41:D571–D578 [View Article]
ThorlandEC, MyersSL, GostoutBS, SmithDI. Common fragile sites are preferential targets for HPV16 integrations in cervical tumors. Oncogene Nature Publishing Group 2003; 22:1225–1237 [View Article]
XuF, CaoM, ShiQ, ChenH, WangY et al. Integration of the full-length HPV16 genome in cervical cancer and CaSKi and SiHa cell lines and the possible ways of HPV integration. Virus Genes Springer US 2015; 50:210–220 [View Article]
HoCM, LeeBH, ChangSF, ChienTY, HuangSH et al. Integration of human papillomavirus correlates with high levels of viral oncogene transcripts in cervical carcinogenesis. Virus Res Elsevier B.V 2011; 161:124–130 [View Article]
LiB, DeweyCN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics BioMed Central 2011; 12:323 [View Article]
RobinsonMD, McCarthyDJ, SmythGK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics2009; 26:139–140 [View Article]
SubramanianA, TamayoP, MoothaVK, MukherjeeS, EbertBL et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA National Academy of Sciences 2005; 102:15545–15550 [View Article]
MoothaVK, LindgrenCM, ErikssonKF, SubramanianA, SihagS et al. PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet2003; 34:267–273 [View Article]
HuangDW, ShermanBT, LempickiRA. Systematic and integrative analysis of large gene lists using David bioinformatics resources. Nature Protocols Nature Publishing Group 2009; 4:44–57 [View Article]
KimD, PerteaG, TrapnellC, PimentelH, KelleyR et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol BioMed Central 2013; 14:R36 [View Article]
ParisDBBP, KuijkEW, RoelenBAJ, StoutTAE. Establishing reference genes for use in real-time quantitative PCR analysis of early equine embryos. Reprod Fertil Dev CSIRO PUBLISHING 2011; 23:353–363 [View Article]
CappelliK, FelicettiM, CapomaccioS, SpinsantiG, SilvestrelliM et al. Exercise induced stress in horses: selection of the most stable reference genes for quantitative RT-PCR normalization. BMC Mol Biol BioMed Central 2008; 9:49 [View Article]
VandesompeleJ, De PreterK, PattynF, PoppeB, Van RoyN et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol BioMed Central 2002; 3:research0034.1 [View Article]
RNA-seq analysis in equine papillomavirus type 2-positive carcinomas identifies affected pathways and potential cancer markers as well as viral gene expression and splicing events