# Non-invasive host transcriptome and HPV oncogene expression map the molecular landscape of HPV-driven cervical lesions

**Authors:** Mohamad Ammar Ayass, Naila Zaman, Natalya Griko, Victor Pashkov, Kevin Zhu, Ghulam Abbas, Melesse Ghelan, Ramya Ramankutty Nair, Tutku Okyay, Jin Zhang, Lina Abi-Mosleh

PMC · DOI: 10.1371/journal.pone.0339693 · PLOS One · 2026-03-19

## TL;DR

This study uses transcriptomic profiling to identify molecular changes in cervical lesions caused by HPV, offering new biomarkers for better risk assessment and treatment.

## Contribution

The study introduces stage-specific transcriptomic signatures and integrates HPV oncogene expression to improve cervical cancer screening.

## Key findings

- NILM showed upregulated growth signaling and immune suppression, indicating stealth infection.
- HSIL was marked by cell cycle hyperactivation and epithelial dedifferentiation, showing molecular discontinuity from earlier stages.
- Higher E6/E7 mRNA levels correlated with more severe lesions, suggesting their use as biomarkers.

## Abstract

Cervical cancer remains a significant global health burden. Current screening methods are not yet capable of detecting molecular alterations preceding cytological abnormalities. In this study, we performed integrative transcriptomic profiling of 132 HPV-positive cervical Pap Smear specimens (NILM, ASCUS, LSIL, HSIL), combining HPV genotyping, and E6/E7 mRNA quantification to map molecular progression. Our analysis revealed stage-specific signatures: NILM displayed a “stealth infection” profile marked by upregulated protein synthesis and growth signaling (EGFR/ERBB2) alongside immune suppression. ASCUS presented a critical tipping point with introduction of early oncogenic drivers (CCND1, SHH), while LSIL prioritized viral productivity with suppressed antimicrobial defenses (MPO, DEFA1). HSIL was distinct from earlier stages and defined by cell cycle hyperactivation (CDK1, PLK1), replication licensing (MCMs), and epithelial dedifferentiation. Pathway crosstalk analysis demonstrated minimal overlap between HSIL and earlier stages (OC < 0.07), highlighting molecular discontinuity during malignant transformation. Additionally, E6/E7 mRNA Ct levels were significantly associated with lesion severity (X2 = 24.407, df = 9, p = 0.003), indicating higher viral mRNA expression are associated with more severe cytological abnormalities. These findings highlight the transformative potential of transcriptomic profiling in cervical cancer prevention, offering stage-specific biomarkers to refine risk stratification. By integrating transcriptomics profiling with current clinical testing, clinicians can distinguish transient infections from high-risks lesions likely to progress. This combined approach addresses the critical limitations of morphology and DNA-based methods, enabling more precise therapeutic interventions and reducing unnecessary overtreatment, and the risk of undertreatment or dismissal of high-risk cases.

## Linked entities

- **Genes:** EGFR (epidermal growth factor receptor) [NCBI Gene 1956], ERBB2 (erb-b2 receptor tyrosine kinase 2) [NCBI Gene 2064], CCND1 (cyclin D1) [NCBI Gene 595], SHH (sonic hedgehog signaling molecule) [NCBI Gene 6469], MPO (myeloperoxidase) [NCBI Gene 4353], DEFA1 (defensin alpha 1) [NCBI Gene 1667], CDK1 (cyclin dependent kinase 1) [NCBI Gene 983], PLK1 (polo like kinase 1) [NCBI Gene 5347]
- **Diseases:** cervical cancer (MONDO:0002974)

## Full-text entities

- **Genes:** CXCR4 (C-X-C motif chemokine receptor 4) [NCBI Gene 7852] {aka CD184, D2S201E, FB22, HM89, HSY3RR, LCR1}, NEFH (neurofilament heavy chain) [NCBI Gene 4744] {aka CMT2CC, NFH}, DEFA4 (defensin alpha 4) [NCBI Gene 1669] {aka DEF4, HNP-4, HP-4, HP4}, DLG1 (discs large MAGUK scaffold protein 1) [NCBI Gene 1739] {aka DLGH1, SAP-97, SAP97, hdlg}, CD9 (CD9 molecule) [NCBI Gene 928] {aka BTCC-1, DRAP-27, MIC3, MRP-1, TSPAN-29, TSPAN29}, IL18 (interleukin 18) [NCBI Gene 3606] {aka IGIF, IL-18, IL-1g, IL1F4}, KIF2C (kinesin family member 2C) [NCBI Gene 11004] {aka CT139, KNSL6, MCAK}, H2BC5 (H2B clustered histone 5) [NCBI Gene 3017] {aka H2B.1B, H2B/a, H2B/b, H2B/g, H2B/h, H2B/k}, CMPK2 (cytidine/uridine monophosphate kinase 2) [NCBI Gene 129607] {aka IBGC10, NDK, TMPK2, TYKi, UMP-CMPK2}, RPS6 (ribosomal protein S6) [NCBI Gene 6194] {aka S6, eS6}, MRM2 (mitochondrial rRNA methyltransferase 2) [NCBI Gene 29960] {aka FJH1, FTSJ2, HEL97, MTDPS17, RRMJ2}, CXCL1 (C-X-C motif chemokine ligand 1) [NCBI Gene 2919] {aka FSP, GRO1, GROa, MGSA, MGSA-a, NAP-3}, IL17A (interleukin 17A) [NCBI Gene 3605] {aka CTLA-8, CTLA8, IL-17, IL-17A, IL17, ILA17}, CXCL8 (C-X-C motif chemokine ligand 8) [NCBI Gene 3576] {aka GCP-1, GCP1, IL8, LECT, LUCT, LYNAP}, IL22 (interleukin 22) [NCBI Gene 50616] {aka IL-21, IL-22, IL-D110, IL-TIF, ILTIF, TIFIL-23}, KRT1 (keratin 1) [NCBI Gene 3848] {aka AEI2, CK1, EHK, EHK1, EPPK, K1}, PTGS2 (prostaglandin-endoperoxide synthase 2) [NCBI Gene 5743] {aka COX-2, COX2, GRIPGHS, PGG/HS, PGHS-2, PHS-2}, RPL8 (ribosomal protein L8) [NCBI Gene 6132] {aka L8, uL2}, KDR (kinase insert domain receptor) [NCBI Gene 3791] {aka CD309, FLK1, VEGFR, VEGFR2}, PLK1 (polo like kinase 1) [NCBI Gene 5347] {aka PLK, STPK13}, FCER1G (Fc epsilon receptor Ig) [NCBI Gene 2207] {aka FCRG}, MTOR (mechanistic target of rapamycin kinase) [NCBI Gene 2475] {aka FRAP, FRAP1, FRAP2, RAFT1, RAPT1, SKS}, PCNA (proliferating cell nuclear antigen) [NCBI Gene 5111] {aka ATLD2}, IFI6 (interferon alpha inducible protein 6) [NCBI Gene 2537] {aka 6-16, FAM14C, G1P3, IFI-6-16, IFI616}, TYROBP (transmembrane immune signaling adaptor TYROBP) [NCBI Gene 7305] {aka DAP12, KARAP, PLOSL, PLOSL1}, POLR1B (RNA polymerase I subunit B) [NCBI Gene 84172] {aka A135, RPA135, RPA2, Rpo1-2, TCS4}, NT5E (5'-nucleotidase ecto) [NCBI Gene 4907] {aka CALJA, CD73, E5NT, NT, NT5, NTE}, SLC9A7 (solute carrier family 9 member A7) [NCBI Gene 84679] {aka MRX108, NHE-7, NHE7}, RFC4 (replication factor C subunit 4) [NCBI Gene 5984] {aka A1, MRMNS, RFC37}, SHH (sonic hedgehog signaling molecule) [NCBI Gene 6469] {aka HHG1, HLP3, HPE3, MCOPCB5, SMMCI, ShhNC}, CENPF (centromere protein F) [NCBI Gene 1063] {aka CENF, CILD31, PRO1779, STROMS, hcp-1}, DSG1 (desmoglein 1) [NCBI Gene 1828] {aka CDHF4, DG1, DSG, EPKHE, EPKHIA, PPKS1}, TXK (TXK tyrosine kinase) [NCBI Gene 7294] {aka BTKL, PSCTK5, PTK4, RLK, TKL}, CUL7 (cullin 7) [NCBI Gene 9820] {aka 3M1, CUL-7, KIAA0076, dJ20C7.5}, FCGR1A (Fc gamma receptor Ia) [NCBI Gene 2209] {aka CD64, CD64A, FCG1, FCGR1, FCRI, FcgammaRI}, ABHD5 (abhydrolase domain containing 5, lysophosphatidic acid acyltransferase) [NCBI Gene 51099] {aka CGI58, IECN2, NCIE2}, CXCL9 (C-X-C motif chemokine ligand 9) [NCBI Gene 4283] {aka CMK, Humig, MIG, SCYB9, crg-10}, KPNA2 (karyopherin subunit alpha 2) [NCBI Gene 3838] {aka IPOA1, PTAC58, QIP2, RCH1, SRP1-alpha, SRP1alpha}, WFDC12 (WAP four-disulfide core domain 12) [NCBI Gene 128488] {aka C20orf122, SWAM2, WAP2, dJ211D12.4}, FN1 (fibronectin 1) [NCBI Gene 2335] {aka CIG, ED-B, FINC, FN, FNZ, GFND}, NUCLEOLIN (nucleolin multifunctional protein) [NCBI Gene 4691] {aka C23, NCL, Nsr1}, CCRL2 (C-C motif chemokine receptor like 2) [NCBI Gene 9034] {aka ACKR5, CKRX, CRAM, CRAM-A, CRAM-B, HCR}, CTPS2 (CTP synthase 2) [NCBI Gene 56474] {aka GATD5B}, CTPS1 (CTP synthase 1) [NCBI Gene 1503] {aka CTPS, GATD5, GATD5A, IMD24}, SPINK5 (serine peptidase inhibitor Kazal type 5) [NCBI Gene 11005] {aka LEKTI, LETKI, NETS, NS, VAKTI}, LCE3A (late cornified envelope 3A) [NCBI Gene 353142] {aka LEP13}, LCE2D (late cornified envelope 2D) [NCBI Gene 353141] {aka LEP12, SPRL1A}, CRYGC (crystallin gamma C) [NCBI Gene 1420] {aka CCL, CRYG3, CTRCT2}, NFKB1 (nuclear factor kappa B subunit 1) [NCBI Gene 4790] {aka CVID12, EBP-1, KBF1, NF-kB, NF-kB1, NF-kappa-B1}, MYC (MYC proto-oncogene, bHLH transcription factor) [NCBI Gene 4609] {aka MRTL, MYCC, bHLHe39, c-Myc}, PALB2 (partner and localizer of BRCA2) [NCBI Gene 79728] {aka BROVCA5, FANCN, PNCA3}, KLK14 (kallikrein related peptidase 14) [NCBI Gene 43847] {aka KLK-L6}, RAD51 (RAD51 recombinase) [NCBI Gene 5888] {aka BRCC5, FANCR, HRAD51, HsRad51, HsT16930, MRMV2}, MGLL (monoglyceride lipase) [NCBI Gene 11343] {aka HU-K5, HUK5, MAGL, MGL}, GMPS (guanosine monophosphate synthase) [NCBI Gene 8833] {aka GATD7}, UBE2C (ubiquitin conjugating enzyme E2 C) [NCBI Gene 11065] {aka UBCH10, dJ447F3.2}, TNFAIP6 (TNF alpha induced protein 6) [NCBI Gene 7130] {aka TSG-6, TSG6}, GOSR2 (golgi SNAP receptor complex member 2) [NCBI Gene 9570] {aka Bos1, EPM6, GS27, MYOS}, TGM3 (transglutaminase 3) [NCBI Gene 7053] {aka TGE, UHS2}, ESR1 (estrogen receptor 1) [NCBI Gene 2099] {aka ER, ESR, ESRA, ESTRR, Era, NR3A1}
- **Diseases:** cervical lesion (MESH:D002575), HPV (MESH:D030361), infected (MESH:D007239), carcinogenesis (MESH:D063646), deaths (MESH:D003643), Cancer (MESH:D009369), prostate cancer (MESH:D011471), invasive cancer (MESH:D009362), viremia (MESH:D014766), Cervical cancer (MESH:D002583), HSIL (MESH:D000081483), dysplasia (MESH:D015792), bladder cancer (MESH:D001749), breast cancer (MESH:D001943), lymph node metastasis (MESH:D008207), proteoglycans (MESH:C536201), Pap (MESH:C535787), Cervical Intraepithelial Neoplasia (MESH:D002578), ASCUS (MESH:D065309), cytological abnormalities (MESH:D000014), inflammation (MESH:D007249), carcinogenic (MESH:D011230), Squamous (MESH:D002294)
- **Chemicals:** C34881 (-), reactive oxygen species (MESH:D017382), Pap (MESH:D010724), SYBR  Green (MESH:C098022), lipid (MESH:D008055), nucleotide (MESH:D009711), water (MESH:D014867), TRIzol (MESH:C411644)
- **Species:** Human papillomavirus (species) [taxon 10566], Halorubrum sp. PV6 (species) [taxon 634157], Human papillomavirus 16 (serotype) [taxon 333760], human papillomavirus 35 (serotype) [taxon 10587], human papillomavirus 11 (serotype) [taxon 10580], Homo sapiens (human, species) [taxon 9606]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC13001926/full.md

## Figures

5 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13001926/full.md

## References

72 references — full list in the complete paper: https://tomesphere.com/paper/PMC13001926/full.md

---
Source: https://tomesphere.com/paper/PMC13001926