Genetic signatures for lineage/sublineage classification of HPV16, 18, 52 and 58 variants

Comparative Analysis of HPV16 Variants in the Untranslated Regulatory Region, L1, and E6 Genes among Vaccinated and Unvaccinated Young Women: Assessing Vaccine Efficacy and Viral Diversity

HPV16 is occasionally detected in vaccinated women who received the bivalent HPV16/18 vaccine, usually at low viral loads. This study explored potential differences in HPV16 variants between vaccinated and unvaccinated women. HPV16-postive viral loads were detected in 1.9% (17/875) and 13% (162/760) of vaccinated and unvaccinated women, respectively, showcasing the vaccine's high efficacy. The L1, E6, and URR regions of HPV16 were sequenced from genital swabs from 16 vaccinated and 25 unvaccinated women in the HAVANA (HPV Among Vaccinated And Non-vaccinated Adolescents) study. The majority of HPV16 variants from vaccinated and unvaccinated women clustered similarly with sub-lineages A1 and A2. Additionally, a separate cluster within lineage A was found, with the variants sharing the L1-located SNP A753G (synonymous) and the URR-located SNP T340C, which did not occur in the other variants. Furthermore, four variants from vaccinated women had relatively long branches, but were not characterized by specific SNPs. The frequency of G712A in the URR was the only SNP observed to be marginally higher among vaccinated women than unvaccinated women. Non-synonymous SNPs T266A in the FG-loop of L1 and L83V in E6 were common among variants from vaccinated and unvaccinated women, but present in similar frequencies. In conclusion, the detection of HPV16 in vaccinated (and unvaccinated) women seemed to be the result of random circulation within this study population.

Genetic diversity and bioinformatic analysis in the L1 gene of HPV genotypes 31, 33, and 58 circulating in women with normal cervical cytology

BACKGROUND

HPV-31, -33, and -58, along with HPV-45 and -52, account for almost 11% of HPV-associated cancers. Our previous studies showed that after HPV-16 and -51, HPV-58 was common and HPV-31 was as frequent as HPV-18 among Iranian women with normal cytology. Hence, in this study, we aimed to investigate the intra-type variations in L1 genes of HPV-58, -31, and -33 to find the predominant lineages circulating in women with normal cytology.

METHODS

Complete coding sequencing of the L1 gene was amplified and nucleotide and amino acid sequences were compared to those of the references. The selective pressure on L1 protein and whether the variations of the L1 genes embed in L1 loops, or N-glycosylated sites were also investigated.

RESULTS

B1, A, and A1 (sub)lineages were common in the HPV-58, -33, and -31 samples, respectively. Ninety nucleotide mutations were observed. Twenty nine nucleotide changes corresponded to nonsynonymous substitutions in which seventeen mutations were located in L1 loops. Only one codon position in HPV-58 sequences was found as the positive selection. No difference was observed in N-glycosylation sites between reference and understudied amino acid sequences.

CONCLUSION

In the current study, we reported, for the first time, the (sub) lineages, amino acid, and genetic diversity in the L1 gene of circulating HPV-58, -33, and -31, in women with normal cytology, in Iran. Such studies can not only have epidemiological values, but also aid to set vaccination programs.

Malignancy and viral infections in Sub-Saharan Africa: A review

The burden of malignancy related to viral infection is increasing in Sub-Saharan Africa (SSA). In 2018, approximately 2 million new cancer cases worldwide were attributable to infection. Prevention or treatment of these infections could reduce cancer cases by 23% in less developed regions and about 7% in developed regions. Contemporaneous increases in longevity and changes in lifestyle have contributed to the cancer burden in SSA. African hospitals are reporting more cases of cancer related to infection (e.g., cervical cancer in women and stomach and liver cancer in men). SSA populations also have elevated underlying prevalence of viral infections compared to other regions. Of 10 infectious agents identified as carcinogenic by the International Agency for Research on Cancer, six are viruses: hepatitis B and C viruses (HBV and HCV, respectively), Epstein-Barr virus (EBV), high-risk types of human papillomavirus (HPV), Human T-cell lymphotropic virus type 1 (HTLV-1), and Kaposi's sarcoma herpesvirus (KSHV, also known as human herpesvirus type 8, HHV-8). Human immunodeficiency virus type 1 (HIV) also facilitates oncogenesis. EBV is associated with lymphomas and nasopharyngeal carcinoma; HBV and HCV are associated with hepatocellular carcinoma; KSHV causes Kaposi's sarcoma; HTLV-1 causes T-cell leukemia and lymphoma; HPV causes carcinoma of the oropharynx and anogenital squamous cell cancer. HIV-1, for which SSA has the greatest global burden, has been linked to increasing risk of malignancy through immunologic dysregulation and clonal hematopoiesis. Public health approaches to prevent infection, such as vaccination, safer injection techniques, screening of blood products, antimicrobial treatments and safer sexual practices could reduce the burden of cancer in Africa. In SSA, inequalities in access to cancer screening and treatment are exacerbated by the perception of cancer as taboo. National level cancer registries, new screening strategies for detection of viral infection and public health messaging should be prioritized in SSA's battle against malignancy. In this review, we discuss the impact of carcinogenic viruses in SSA with a focus on regional epidemiology.

Analysis of the genomic diversity of human papillomavirus type 31 in cervical samples reveals the presence of novel sublineages in clade C

Human papillomavirus 31 (HPV31) is the fourth most frequent high-risk HPV (HR-HPV) genotype identified in cervical cancer (CC) worldwide and in Mexico. It has been recently classified into three lineages (A, B, and C) and eight sublineages (A1, A2, B1, B2, and C1 - C4). Here, we report the complete genomic sequences of 14 HPV31 isolates from cervical samples, and these were compared with viral genome sequences from the GenBank database for phylogenetic and genetic distance analysis. The formation of two novel clades within the C lineage (proposed as C5 and C6) was observed, with a well-defined variant-specific mutational pattern. The smallest average pairwise distance was 0.71% for lineages A and B, 0.94% for lineages A and C, and 1.01% for lineages B and C, and between sublineages, these values were 0.21% for clade A, 0.29% for clade B, and 0.24% for clade C. The isolates were grouped into the sublineages A1, B2, C1-C3, and C6. This is the first report on the whole-genome diversity of HPV31 in Mexico.

Concordance between the BD Onclarity and Roche cobas assays for detection of HPV DNA in a Chinese population

As cervical cancer screening shifts from cytology to human papillomavirus (HPV) testing, a major issue involves validating more HPV tests. In recent years, some HPV tests are used for clinical performance verification in China. The purpose of this study was to explore whether the BD Onclarity (Becton, Dickinson and Company)HPV assay differs from the Roche cobas (Roche Molecular Systems)HPV assay, as determined using 944 cervical samples, including 588 with sequencing results. In the nucleic acid assay accuracy verification, the assays showed excellent concordance for detection of HPV16 (κ = 0.93, 95% confidence interval [CI]: 0.89-0.97) and HPV18 (κ = 0.90, 95% CI: 0.83-0.97), and very good concordance for the 12 other high-risk types (HPV31/33/35/39/45/51/52/56/58/59/66/68, κ = 0.79, 95% CI: 0.75-0.83). The overall agreement for HPV DNA detection between Onclarity and cobas was very good (κ = 0.7755). No difference for ≥CIN2 sensitivity was observed between Onclarity and cobas (both 96.5%), whereas the ≥CIN2 specificity for detection of Onclarity (16.6%, 95% CI: 13.7-19.9) was higher than that of cobas (11.5%, 95% CI: 9.1-14.5). Onclarity exhibited comparable screening performance and triage efficiency compared to cobas in the detection of cervical disease in Chinese women.

Machine Learning Prediction of Adenovirus D8 Conjunctivitis Complications from Viral Whole-Genome Sequence

Objective

To obtain complete DNA sequences of adenoviral (AdV) D8 genome from patients with conjunctivitis and determine the relation of sequence variation to clinical outcomes.

Design

This study is a post hoc analysis of banked conjunctival swab samples from the BAYnovation Study, a previously conducted, randomized controlled clinical trial for AdV conjunctivitis.

Participants

Ninety-six patients with AdV D8-positive conjunctivitis who received placebo treatment in the BAYnovation Study were included in the study.

Methods

DNA from conjunctival swabs was purified and subjected to whole-genome viral DNA sequencing. Adenovirus D8 variants were identified and correlated with clinical outcomes, including 2 machine learning methods.

Main Outcome Measures

Viral DNA sequence and development of subepithelial infiltrates (SEIs) were the main outcome measures.

Results

From initial sequencing of 80 AdV D8-positive samples, full adenoviral genome reconstructions were obtained for 71. A total of 630 single-nucleotide variants were identified, including 156 missense mutations. Sequence clustering revealed 3 previously unappreciated viral clades within the AdV D8 type. The likelihood of SEI development differed significantly between clades, ranging from 83% for Clade 1 to 46% for Clade 3. Genome-wide analysis of viral single-nucleotide polymorphisms failed to identify single-gene determinants of outcome. Two machine learning models were independently trained to predict clinical outcome using polymorphic sequences. Both machine learning models correctly predicted development of SEI outcomes in a newly sequenced validation set of 16 cases (P = 1.5 × 10-5). Prediction was dependent on ensemble groups of polymorphisms across multiple genes.

Conclusions

Adenovirus D8 has ≥ 3 prevalent molecular substrains, which differ in propensity to result in SEIs. Development of SEIs can be accurately predicted from knowledge of full viral sequence. These results suggest that development of SEIs in AdV D8 conjunctivitis is largely attributable to pathologic viral sequence variants within the D8 type and establishes machine learning paradigms as a powerful technique for understanding viral pathogenicity.

A Straightforward HPV16 Lineage Classification Based on Machine Learning

Human Papillomavirus (HPV) is the causal agent of 5% of cancers worldwide and the main cause of cervical cancer and it is also associated with a significant percentage of oropharyngeal and anogenital cancers. More than 60% of cervical cancers are caused by HPV16 genotype, which has been classified into lineages (A, B, C, and D). Lineages are related to the progression of cervical cancer and the current method to assess lineages is by building a Maximum Likelihood Tree (MLT); which is slow, it cannot assess poor sequenced samples, and annotation is done manually. In this study, we have developed a new model to assess HPV16 lineage using machine learning tools. A total of 645 HPV16 genomes were analyzed using Genome-Wide Association Study (GWAS), which identified 56 lineage-specific Single Nucleotide Polymorphisms (SNPs). From the SNPs found, training-test models were constructed using different algorithms such as Random Forest (RF), Support Vector Machine (SVM), and K-nearest neighbor (KNN). A distinct set of HPV16 sequences (n = 1,028), whose lineage was previously determined by MLT, was used for validation. The RF-based model allowed a precise assignment of HPV16 lineage, showing an accuracy of 99.5% in the known lineage samples. Moreover, the RF model could assess lineage to 273 samples that MLT could not determine. In terms of computer consuming time, the RF-based model was almost 40 times faster than MLT. Having a fast and efficient method for assigning HPV16 lineages, could facilitate the implementation of lineage classification as a triage or prognostic marker in the clinical setting.

Epstein-Barr Virus in Burkitt Lymphoma in Africa Reveals a Limited Set of Whole Genome and LMP-1 Sequence Patterns: Analysis of Archival Datasets and Field Samples From Uganda, Tanzania, and Kenya

Epstein-Barr virus (EBV) is associated with endemic Burkitt lymphoma (eBL), but the contribution of EBV variants is ill-defined. Studies of EBV whole genome sequences (WGS) have identified phylogroups that appear to be distinct for Asian versus non-Asian EBV, but samples from BL or Africa, where EBV was first discovered, are under-represented. We conducted a phylogenetic analysis of EBV WGS and LMP-1 sequences obtained primarily from BL patients in Africa and representative non-African EBV from other conditions or regions using data from GenBank, Sequence Read Archive, or Genomic Data Commons for the Burkitt Lymphoma Genome Sequencing Project (BLGSP) to generate data to support the use of a simpler biomarker of geographic or phenotypic associations. We also investigated LMP-1 patterns in 414 eBL cases and 414 geographically matched controls in the Epidemiology of Burkitt Lymphoma in East African children and minors (EMBLEM) study using LMP-1 PCR and Sanger sequencing. Phylogenetic analysis revealed distinct genetic patterns of African versus Asian EBV sequences. We identified 281 single nucleotide variations (SNVs) in LMP-1 promoter and coding region, which formed 12 unique patterns (A to L). Nine patterns (A, AB, C, D, F, I, J, K and L) predominated in African EBV, of which four were found in 92% of BL samples (A, AB, D, and H). Predominant patterns were B and G in Asia and H in Europe. EBV positivity in peripheral blood was detected in 95.6% of EMBLEM eBL cases versus 79.2% of the healthy controls (odds ratio [OR] =3.83; 95% confidence interval 2.06-7.14). LMP-1 was successfully sequenced in 66.7% of the EBV DNA positive cases but in 29.6% of the controls (ORs ranging 5-11 for different patterns). Four LMP-1 patterns (A, AB, D, and K) were detected in 63.1% of the cases versus 27.1% controls (ORs ranges: 5.58-11.4). Dual strain EBV infections were identified in WGS and PCR-Sanger data. In conclusion, EBV from Africa is phylogenetically separate from EBV in Asia. Genetic diversity in LMP-1 formed 12 patterns, which showed promising geographic and phenotypic associations. Presence of multiple strain infection should be considered in efforts to refine or improve EBV markers of ancestry or phenotype.

Lay Summary

Epstein-Barr virus (EBV) infection, a ubiquitous infection, contributes to the etiology of both Burkitt Lymphoma (BL) and nasopharyngeal carcinoma, yet their global distributions vary geographically with no overlap. Genomic variation in EBV is suspected to play a role in the geographical patterns of these EBV-associated cancers, but relatively few EBV samples from BL have been comprehensively studied. We sought to compare phylogenetic patterns of EBV genomes obtained from BL samples in Africa and from tumor and non-tumor samples from elsewhere. We concluded that EBV obtained from BL in Africa is genetically separate from EBV in Asia. Through comprehensive analysis of nucleotide variations in EBV's LMP-1 gene, we describe 12 LMP-1 patterns, two of which (B and G) were found mostly in Asia. Four LMP-1 patterns (A, AB, D, and F) accounted for 92% of EBVs sequenced from BL in Africa. Our results identified extensive diversity of EBV, but BL in Africa was associated with a limited number of variants identified, which were different from those identified in Asia. Further research is needed to optimize the use of PCR and sequencing to study LMP-1 diversity for classification of EBV variants and for use in epidemiologic studies to characterize geographic and/or phenotypic associations of EBV variants with EBV-associated malignancies, including eBL.

Genetic Diversity of HPV 16 and HPV 18 Based on Partial Long Control Region in Iranian Women

Background

Human papillomavirus (HPV) 16 and HPV 18 account for 75% of all cervical cancers. The L1 gene, encoding the major surface protein (MSP), is used to classify HPV types (lineages and sublineages), genotypes, and intratypic variants. Therefore, this study aimed to investigate the lineages, sublineages, genetic variabilities, and mutation effects on transcription factor binding sites by using partial sequences of the HPV 16 and HPV 18 long control regions (LCRs) in these samples.

Materials and Methods

After DNA isolation from 56 positive samples, the LCR of HPV 16 and HPV 18 were amplified using specific primers, and phylogenetic trees were drawn through MEGA X. Compared to the reference sequences, single nucleotide polymorphisms (SNPs) were identified. The transcription binding sites were also evaluated using the online PROMO database.

Results

The LCRs of 52 samples were successfully sequenced. Overall, 81.58% of all HPV 16 variants belonged to the D1 sublineage, followed by A4 (13.16%), A1 (2.63%), and C1 (2.63%) sublineages. All HPV 18 isolates belonged to A sublineage, 92.85% to A3 sublineage, and 7.15% to A4 sublineage. Out of 27 SNPs in the HPV 16 LCR, A7382T, T7384G, C7387T, C7393G, A7431G, T7448C, and C7783A were HPV 16-specific. Also, among 14 SNPs in the HPV 18 LCR, C7577A and A7943T were not previously reported. An insertion (C) between 7432 and 7433 positions was identified in all studied HPV 16 variants. Besides, most of the HPV 16 mutations were embedded in the YY1, TFIID, Oct-2, and NF-1 binding sites, while c-Fos and MBF1, as the most common binding sites, were affected by HPV 18 LCR mutations.

Conclusion

The present results showed that D1 and A3 were the dominant sublineages of HPV 16 and HPV 18, respectively. Therefore, women infected with these variants need to be examined in further longitudinal studies to obtain more information about the oncogenic potential of these dominant variants in Iran. Besides, YY1, TFIID, Oct-2, NF-1, c-Fos, and MBF1 were the most frequent binding sites, which were influenced by the mutations.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.