# Decoding the Genetic and Structural Features of HPV16 E5 Oncogene in Cervical Cancer Isolates from Pakistan: A Pilot Study

## TL;DR

This study is the first to report on HPV16 E5 gene mutations in cervical cancer patients from Pakistan, finding a variant with potential increased oncogenic effects.

## Contribution

The first report of HPV16 E5 variants in a Pakistani population, identifying a novel mutation with enhanced interaction with EGFR.

## Key findings

- Isolate 11 has mutations C3979A and G4042A leading to amino acid substitutions in the E5 protein.
- The mutations in isolate 11 increase interaction with EGFR, potentially enhancing downstream signaling.
- Most isolates clustered with the European-Asian lineage of HPV16.

## Abstract

Many anogenital cancers are caused by high-risk HPV. The most common subtype is HPV16, which is prevalent in the world, including Pakistan. Various amino acid residues in HPV16 E5 are associated with high cell cycle progression and proliferation. Lack of studies on HPV16E5 in Pakistan prompted the current study. This is the first report on the occurrence of pathogenic E5 variant of HPV16 in tissue sections obtained from invasive cervical cancerous patients in Pakistan.

A subset of 11 samples from HPV-positive biopsies were subjected to E5 gene amplification using PCR and analyzed using bioinformatics programs. The bioinformatics analysis detected mutations causing structural variations, which potentially contribute to the oncogenic properties of proteins.

The two-point mutations, C3979A and G4042A, observed in isolate 11 caused the substitution of isoleucine for leucine and valine at positions 44 and 65 in E5 protein. The rest of the isolates had Leu44Val65 amino acids. Intratypic variations and phylogenetic analysis revealed that the majority of the isolates were closely clustered with European-Asian lineage. Moreover, C3979A and G4042A contributed to higher degree of interactions with host receptors, i.e. EGFR.

This is the first study reporting HPV16 variants in a Pakistani population based on variations in the E5 region. Our findings indicate that isolate 11 has a strong interaction with the intracellular domain of EGFR, which may enhance the generation of downstream signals. Since this was a pilot study to explore E5 gene mutation, future studies with large samples are absolutely needed.