Mutational landscape and intra-host diversity of human papillomavirus type 16 long control region and E6 variants in cervical samples.

Human papillomavirus genotype 16 (HPV16) is the most frequent high-risk HPV (HR-HPV) identified in cervical precursor lesions and cervical cancer (CC) worldwide. The oncogenic potential of HPV16 is partly dependent on the lineage involved in the infection and the presence of clinically relevant mutations. In this report, we present the distribution of HR-HPV and the mutational profile and intra-host variability of HPV16 lineages, based on analysis of the long control region (LCR) and the E6 gene in samples with normal cytology (n = 39), squamous intraepithelial lesions (n = 25), and CC (n = 39). HR-HPV genotyping was performed using multiplex real-time PCR. HPV16 lineage assignments and mutation frequencies were determined by conventional PCR and Sanger DNA sequencing, and intra-patient viral populations were analyzed using next-generation sequencing (NGS). The most frequent HR-HPV type was HPV16, followed by HPV31 and HPV18. The frequency of HPV16 sublineages was A1/A2 > D2 > D3 and B1. Moreover, the most frequent mutations, both in samples from this study and in the available sequences from Mexican isolates in the GenBank database were LCR-G7518A, which is involved in carcinogenesis, and E6-T350G (producing L83V), associated with persistence of infection. Otherwise, deep sequencing revealed high conservation of viral lineages and mutations, independently of the stages studied. In conclusion, the high frequency and stability of these molecular markers, as well as the circulating viral lineages, could be related to the incidence of CC associated with HPV16. Hence, they deserve a broader analysis to determine the risk of specific populations for progression of the disease.