Skip to main content

A comprehensive HPV-STI NGS assay for detection of 29 HPV types and 14 non-HPV sexually transmitted infections



Sexually transmitted infections (STIs) are prevalent throughout the world and impose a significant burden on individual health and public health systems. Missed diagnosis and late treatment of STIs can lead to serious complications such as infertility and cervical cancer. Although sexually transmitted co-infections are common, most commercial assays target one or a few STIs. The HPV-STI ChapterDx Next Generation Sequencing (NGS) assay detects and quantifies 29 HPVs and 14 other STIs in a single-tube and single-step PCR reaction and can be applied to tens to thousands of samples in a single sequencing run.


A cohort of 274 samples, previously analyzed by conventional cytology/histology and Roche cobas HPV Test, were analyzed by ChapterDx HPV-STI NGS assay for detection of 43 HPV and STI. A set of 43 synthetic control DNA fragments for 43 HPV and STI were developed to evaluate the limit of detection, specificity, and sensitivity of ChapterDx HPV-STI NGS assay.


The assay was evaluated in this study, and the limit of detection was 100% at 50 copies for all targets, and 100%, 96%, 88% at 20 copies for 34, 8, and 1 target, respectively. The performance of this assay has been compared to Roche cobas HPV test, showing an overall agreement of 97.5% for hr-HPV, and 98.5% for both, HPV16 and HPV18. The assay also detected all HPV-infected CIN2/3 with 100% agreement with Roche cobas HPV results. Moreover, several co-infections with non-HPV STIs, such as C. trachomatis, T. vaginalis, M. genitalium, and HSV2 were identified.


The ChapterDx HPV-STI NGS assay is a user-friendly, easy to automate and cost-efficient assay, which provides accurate and comprehensive results for a wide spectrum of HPVs and STIs.


Sexually transmitted infections (STIs) are prevalent globally and World Health Organization (WHO) estimates that annually there are 376 million new cases of four curable STIs: Chlamydia trachomatis, Neisseria gonorrhoeae, Treponema pallidum (syphilis), and Trichomonas vaginalis. Prevalence of several viral STIs is similarly high with an estimated 500 million people infected with herpes simplex type 2 (HSV2), and approximately 291 million women harboring human papillomaviruses (HPV) [1]. STI prevalence varies by region, gender. STI epidemics cause a profound impact on physical and psychological health worldwide [1, 2].

HPVs are a large group of viruses with more than 200 completely characterized types, and new HPV types being continuously found [3]. Among them, about 40 types known to infect human anogenital tract are grouped into high risk (hr)-HPV and low risk (lr)-HPV types based on their oncogenic potential [4, 5]. Hr-HPVs are the causative agents of cervical cancer and have been detected in 99.7% of cervical cancers [6, 7]. Fourteen hr-HPV types are considered to cause most cervical cancers worldwide. Nevertheless, other hr-types have also been reported to cause cervical cancer [4, 5]. HPV oncogenic potential is highly type-dependent as some types are more oncogenic. Many studies have shown that extended genotyping information holds important clinical value for patient management and triage of HPV positive women [8,9,10].

HPV viral load quantification has been shown to be a predictor of infection persistence or clearance [11]. It has also been reported that HPV viral load is associated with risk of persistent infection and precancer [12,13,14,15], indicating viral load as a candidate risk marker [16]. Currently, viral load is measured mainly by quantitative polymerase chain reaction (qPCR) on a limited number of hr-HPV types such as HPV16 and HPV18. An assay that can detect and quantify many HPVs in a single-tube PCR reaction is a needed tool for association studies of viral load and precancer.

Besides HPVs, several bacteria and viruses infect human anogenital area. Earlier studies indicated that non-HPV STIs could facilitate the entry of HPV virions into host cells, change in the immunological response pathways, and decrease the host’s ability to clear HPV infection [16,17,18]. Co-infections between hr-HPVs and non-HPV STIs such as Chlamydia trachomatis, Neisseria gonorrhoeae, Trichomonas vaginalis, Ureaplasma urealyticum, and HSV2 have been reported to be associated with HPV persistence, cervical dysplastic and neoplastic lesion [19,20,21,22,23,24,25,26]. Co-infections between HPVs and non-HPV STIs are widespread. Since most commercial kits detect one or two STIs, multiple tests are required for comprehensive STI identification.

This work presents the evaluation of the ChapterDx HPV-STI NGS assay (NGS assay in short) that amplifies 29 HPVs and 14 STIs (Table 1) with species/type-specific primers in a single-tube and single-step PCR reaction followed by NGS (Fig. 1). This assay can detect and quantify these 43 targets from tens to thousands of samples in a single Illumina sequencing run (depending on the sequencing platform and sequencing kit).

Table 1 The list of 43a HPV/STI in the NGS panel
Fig. 1
figure 1

A schematic illustration of ChapterDx HPV-STI NGS workflow

Materials and methods

Clinical HPV screening specimens

A cohort of 274 blind samples that had previously been analyzed by conventional cytology/histology and Roche cobas HPV Test (Roche, Pleasanton, CA, USA) was selected for this study. Within this cohort, 102, 101 and 57 samples were negative for intraepithelial lesions or malignancy (NILM), low-grade squamous intraepithelial lesion (LSIL) and high-grade squamous intraepithelial lesion (HSIL), respectively; and five were invasive carcinoma and nine samples had no cytology results. The available histological results for this cohort were: cervical intraepithelial neoplasia grade 1 (CIN1) (N = 4), CIN2 (N = 1) and CIN3 (N = 58), 3 squamous cell carcinoma (SCC) and 2 adenocarcinoma (AdC). HPV samples were shipped from Oncogenic Viruses Service, National and Regional HPV Reference Laboratory, National Institute of Infectious Diseases-ANLIS “Dr. Malbrán”, Buenos Aires, Argentina to Chapter Diagnostics in Menlo Park, California, USA, for comprehensive HPV-STI analysis.

Roche cobas HPV testing

The 274-sample cohort were analyzed by Roche cobas HPV automated platform [27, 28]. The amplification of target DNA is performed by real-time PCR technology delivering genotype results for HPV16 and HPV18 (individually), along with a simultaneous, pooled result for other 12 high-risk genotypes. The system uses an internal control (ß-globin cellular) to prevent false negative results.

ChapterDx HPV-STI NGS kit

ChapterDx HPV-STI NGS kit targets 19 hr-HPVs, 10 lr-HPVs, 14 non-HPV STIs, and GAPDH as internal control. Table 1 shows the HPV genotypes and non-HPV STIs covered by the kit.

Controls for analytical performance

A set of 43 synthetic control DNA fragments for all targets were synthesized by Integrated DNA Technologies (Coralville, IA). The controls were designed based on following references derived from PapillomaVirus Episteme (PaVE) database: HPV6 (KU298876.1), HPV11 (KU298879.1), HPV16 (LC193821.1), HPV18 (KU707825.1), HPV26 (KF444058.1), HPV31 (KU163582.1), HPV33 (KU050111.1), HPV35 (JX129488.1), HPV39 (KC470231.1), HPV40 (KU298895.1), HPV42 (JQ902118.1), HPV43 (HE962401.1), HPV44 (HE963128.1), HPV45 (KU049754.1), HPV51 (KU050115.1), HPV52 (KU298909.1), HPV53 (KU951265.1), HPV55 (KU298899.1), HPV56 (KU050122.1), HPV58 (KU550636.1), HPV59 (KC470266.1), HPV61 (KU298923.1), HPV66 (KU298928.1), HPV68a (KC470268.1), HPV68b (KU298932.1), HPV73 (KU298935.1), HPV81 (KU298939.1), HPV82 (AB027021.1), HPV83 (AF151983.1), Chlamydia trachomatis (serovars A to K) (CP018052.1), Chlamydia trachomatis serovar L1 (JN795429.1), Chlamydia trachomatis serovar L2 (JN795428.1), HSV1 (KX946970.1), HSV2 (KU886303.1), VZV (KX262866.1), Haemophilus ducreyi (CP015434.1), Mycoplasma genitalium (CP003773.1), Mycoplasma hominis (CP011538.1), Neisseria gonorrhoeae (LT592163.1), Treponema pallidum (CP016048.1), Trichomonas vaginalis (XM_001582357.1), Ureaplasma parvum (AP014584.1) and Ureaplasma urealyticum (CP041200.1). The DNA controls were pooled and diluted into 100, 50 and 20 copies per vial in a final volume of 5 µl. Each pool was spiked with 10 ng human genomic DNA. The pools were amplified and sequenced in 23 replicates for 100 copies and 24 replicates for 50 and 20 copies to assess the limit of detection (LoD) of the NGS assay. The acceptable consensus for an assay’s LoD at the lowest concentration is ≥ 95% for at which 19/20 replicates are positive [29].

DNA extraction and PCR

Genomic DNA extraction for 274-sample cohort were performed by Roche cobas HPV Test (cobas for short) platform. All the specimens were measured by Qubit 3 (ThermoFisher, CA, USA) to monitor the presence and concentration of DNA prior to PCR. One-step multiplex PCR was performed in 15 µl final volume in a 96-well plate on a Veriti thermocycler (ThermoFisher, CA, USA). The PCR reaction consisted of target-specific primers (29 HPV, 14 STI and internal control), barcoded universal primers, sample DNA, DNA polymerase, dNTPs and PCR buffer. The PCR conditions were according the ChapterDx HPV-STI NGS kit instructions.

Library preparation and next-generation sequencing

After PCR, the amplified products were pooled into a 1.5 ml tube (or a 15 ml tube depending on the sample size). A portion of the pooled amplicons was purified with SPRIbeads (Beckman Coulter, CA, USA) according to the manufacturer’s instructions. The purified pooled amplicons concentration was measured on a Qubit fluorometer, and the concentration was adjusted for sequencing according to Illumina library preparation instructions. The library was sequenced on an Illumina MiniSeq platform using an Illumina Mid-Output sequencing kit (Illumina, CA, USA). The workflow is depicted in Fig. 1.

Sequence data analysis software and interpretation

(1) Decoding: sequencing reads (FASTQ format) for forward and reverse reads and forward and reverse indexing reads are input in ChapterDx Analysis Software. The software assigns sequencing reads for a sample based on the sequence of both forward and reverse index reads with no mismatch; (2) Mapping: sequencing reads are mapped onto reference sequences using the Smith-Waterman algorithm with options as nucleotide match reward is 1, nucleic mismatch penalty is -3, cost to open a gap is 5, and cost to extend a gap is 2. Only the alignment of best match is kept for each sequencing read if the alignment score is beyond 60. Alignments of paired reads with identical reference sequences are kept. Reads are labeled to the reference based on their corresponding alignments; (3) Genotyping: a sample is called positive for a species/type when there are more reads than the preset cutoff reads labelled with the corresponding reference.

NGS detection threshold

For determining the cutoff detection threshold, samples from HPV LabNet international proficiency panel [30] as well as control HPV plasmids (HPV11, 16, 18, 31, 33, 35, 39, 45, 51, 58, 59 and 68) were used. The LabNet panel consists of 43 samples, containing different concentrations of HPV plamids of hr-HPV16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, 68a, 68b and lr-HPV6, 11. Several samples contain 2–4 HPV plasmids to mimic multiple co-infections. Based on ChapterDx HPV-STI NGS sequence data from LabNet international proficiency HPV panel and control HPV plasmids, the threshold for positive samples were set at 3 reads. Samples below the threshold were considered negative.

Statistical analysis

The HPV and STI species/types frequencies, cytology and histology were calculated by Microsoft excel. Overall agreement for 14 hr-HPV, HPV16, 18 and other 12 HPVs between ChapterDx HPV-STI NGS and Roche cobas HPV assays were calculated by Cohen’s Kappa coefficients with 95% confidence intervals (CIs).


Analytical performance of the assay

To evaluate the analytical performance and sensitivity of the assay, a set of 43 synthetic control DNA fragment pools of 100 copies (23 replicates), 50 copies (24 replicates), 20 copies (24 replicates) and 0 copies (blank) spiked with 10 ng human genomic DNA were analyzed to determine the LoD. 100 copies and 50 copies control DNA were detected in all replicates. For the 20 copies pool replicates, the detection rate was 100% (24/24) for 34 targets, 95.8% (23/24) for HPV6, 11, 39, 42, 43, 45, 56 and 83 (8 targets) and 87.5% (21/24) for Mycoplasma genitalium (Tables 2 and 3). The lowest concentration LoD is recommended at ≥ 95% (for 19/20 replicates) [29], therefore, the LoD was determined at 20 copies for all targets except for Mycoplasma genitalium. There was no sequencing reads for the blank sample. The 43 control DNA average sequencing depth for replicates of 100, 50 and 20 copies were 876X, 640X, and 309X (excluding GAPDH internal control).

Table 2 LoD results of 43 control fragment pools for 23 replicates of 100 copies and 24 replicates of 50 and 20 copies
Table 3 20 copies LoD for 9 control fragments below100% agreement

Agreement between ChapterDx HPV-STI NGS assay and Roche cobas HPV Test in cervical samples with cyto-histological diagnosis

To evaluate the agreement with the Roche cobas HPV test, 274 clinical samples were analyzed by ChapterDx HPV-STI NGS assay. Out of 274 samples, the results for 267 samples agreed with those by cobas HPV test for 14 hr-HPVs. Among those 267 samples, 241 were positive and 26 were negative for 14 hr-HPVs. The overall agreement for 14 hr-HPVs between both methods was 97.5% with kappa 0.867. The CIN1 (n = 4), CIN2/3 (n = 59) were positive for hr-HPV by both cobas and the NGS assays (Table 4). Of the 5 cancer samples, one squamous carcinoma sample and one adenocarcinoma sample were negative for 14 hr-HPV by both cobas and the NGS assay. The adenocarcinoma sample was positive for HPV44, 81 by the NGS assay, which was confirmed by reverse line blot assay [31].

Table 4 Overall concordance of 274 clinical samples between ChapterDx HPV-STI NGS assay and Roche cobas for 14 high-risk HPVs, stratified by histology and cytology

The overall agreement for HPV16 between cobas and the NGS assays was 98.5% (270/274) with a kappa value of 0.967 (Table 5). Two of the HPV16 samples were scored negative. However, one and two HPV16 reads were detected in each of these two samples, respectively. The overall agreement for HPV18 was 98.5% (70/274) with a kappa value of (0.901) and the overall agreement for other HPVs was 96% (263/274) with kappa value of 0.906. For the other HPVs, two HPV samples were also detected by the NGS assay as HPV52 and HPV31/52, respectively, but were below the NGS detection threshold and scored negative. Of the 29 HPV targets in the NGS panel, 28 were detected among the samples analyzed (except for HPV73). In this cohort, co-infections between 2 and 12 HPVs were detected in 146 samples (shown in Table 6).

Table 5 Concordance between ChapterDx HPV-STI NGS assay and Roche cobas assay for 14 high-risk HPVs
Table 6 Nr of 2 to 12 HPV co-infections in 274 samples

In the 274-sample cohort, ChapterDx HPV-STI NGS assay detected a total frequency of 205 non-HPV STIs, of which U. parvum (32.1%) and M. hominis (25.2%) comprised the majority of the non-HPV STIs (Table 7). 52.1% (143/274) of the samples were positive for 7 out of 14 non-HPV STIs, harboring a range of 1–4 STIs per sample (Table 8). HPV and non-HPV STIs co-infections were observed in 135/274 (49.3%) of the samples. Table 9 lists results by Roche cobas assay and the NGS assay for 9 samples with co-infections of multiple HPVs and non-HPV STIs, demonstrating detection of HPV and non-HPV STIs by the NGS assay.

Table 7 Other STIs detected in the 274 samples by ChapterDx HPV-STI NGS assay
Table 8 Nr of 2 to 4 non-HPV STI co-infections in 274 samples
Table 9 Comparison of ChapterDX HPV-STI NGS and Roche cobas for samples with combined HPV and STI co-infection

On average, 192 positive samples were sequenced per Illumina Mid-Output sequencing kit on Illumina MiniSeq platform including negative and positive controls. Additional file 1: Table S1 lists the estimated number of samples for different Illumina sequencing platforms and kits.


In the present study, the performance of ChapterDx HPV-STI NGS assay was compared to Roche cobas assay. Forty-three control DNA pools of 100 and 50 copies were detected in all replicates and for 20 copies replicates, M. genitalium had an LoD below 95%. The lower detection rate for M. genitalium could be due to synthetic error in control DNA, which affects amplification efficiency. In this study, the HPV52 control DNA contained several synthesizing sequence errors in the primer regions. The sequence errors were confirmed by Sanger sequencing. The HPV52 control DNA was re-synthesized, and the problem was resolved.

The 3 reads detection was consistent with Labnet HPV panel, control HPV plasmids and sample cohort analyzed by Roche cobas. The NGS assay did not generate background sequences or false positive results, showing that 3 reads is an acceptable detection threshold.

The NGS and cobas assays showed a high overall agreement, which indicates a good performance of the NGS assay on a set of well-studied samples; nevertheless, a larger study is required based on Meijer’s criteria [32] for further clinical evaluation. In this 274-sample cohort, none of the discrepant results involved HSIL, CIN1, CIN2/3 or cancer samples. The adenocarcinoma sample that was negative for 14 hr-HPVs by both methods, harbored HPV44, 81, which are lr-HPV types and not considered carcinogenic. Due to type-specific amplification, the NGS assay was able to detect a range of 2–12 HPV co-infections in 146/274 of samples allowing unbiased analysis of existing HPV types in each sample.

Earlier studies have shown evidence of association between hr-HPVs and non-HPV STIs (including commensal STIs) with cervical carcinogenesis [17,18,19,20,21,22,23,24,25,26, 33]. In this study, 52.1% samples in the 274-sample cohort harbored non-HPV STIs, of which 87.8% were commensal infections of three species: M. hominis, U. urealyticum and U. parvum. Moreover, almost half of the samples carried co-infections of HPV and non-HPV STIs. Thus, a comprehensive assay detecting multiple HPV and non-HPV STIs in a single test, is a useful tool to diagnose and study the correlation of hr-HPV and non-HPV co-infections. It is worth mentioning that the 14 non-HPV STIs in the NGS assay have different pathogenic characteristics and some of them are considered commensal species of the genital tract (such as M. hominis, U. urealyticum and U. parvum), existing in both healthy and symptomatic individuals [34]. Such commensal infections are interpreted in clinical context, which is out of the scope of this study to investigate the clinical impact and association of these commensal STI infections with disease.

HPV genotyping has been mostly focused on HPV16 and HPV18 due to their high prevalence in cervical cancer, although, other hr-types indicate high positive predictive value (PPV) [8, 10]. Earlier studies suggested that extended HPV genotyping provides a better risk stratification and identification of women at increased risk of cervical cancer by simply providing individual risk assessment of HPV positive women [35,36,37]. The results of this study show that co-infections of multiple HPVs and non-HPV STIs are widespread (Tables 7 and 8). Due to limitation of current widely used assays, many of those STIs can go undetected or more than one assay might be needed.

It has been reported that non-transient hr-HPV infections and high viral loads is a risk marker for dysplasia and carcinoma in situ in normal cytology [11,12,13,14,15,16]. Such studies have been typically focused on few high-risk types. The NGS assay uses type-specific HPV primers for amplification and generates quantification information for each type in the panel. The wide spectrum of STIs and their quantification information could be useful for studying the association between STIs and diseases. Investigating the association of broad-spectrum HPVs, their respective viral load, role of other STIs in relation to cancer would be the scope of another study.

There are several NGS-based HPV assays [38,39,40] that commonly use more than one round of PCR. ChapterDx HPV-STI NGS assay uses one-step amplification workflow, which significantly reduces time/labor, the risk of cross contamination and allows easier automation.


ChapterDx HPV-STI NGS Assay is cost-effective as it detects a much broader spectrum of HPV/STI than other commercial kits, at a similar or lower cost per sample, with proven analytical performance; therefore, it is a valuable tool for large epidemiology and surveillance studies. Once evaluated with larger, well-characterized cohorts for clinical validation, this assay could also be applied to cervical cancer screening programs.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. (WHO) WHO: Global health sector strategy on sexually transmitted infections, 2016–2021. 2018.

  2. Eisinger RW, Erbelding E, Fauci AS. Refocusing research on sexually transmitted infections. J Infect Dis. 2020;222(9):1432–4.

    Article  Google Scholar 

  3. Bernard HU, Burk RD, Chen Z, van Doorslaer K, zur Hausen H, de Villiers EM. Classification of papillomaviruses (PVs) based on 189 PV types and proposal of taxonomic amendments. Virology. 2010;401(1):70–9.

    Article  CAS  Google Scholar 

  4. Humans IWGotEoCRt: Biological agents. Volume 100 B. A review of human carcinogens. IARC Monogr Eval Carcinog Risks Hum. 2012;100(Pt B):1–441.

  5. Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, Snijders PJ, Meijer CJ. International Agency for Research on Cancer Multicenter Cervical Cancer Study G: epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348(6):518–27.

    Article  Google Scholar 

  6. Munoz N, Bosch FX, de Sanjose S, Tafur L, Izarzugaza I, Gili M, Viladiu P, Navarro C, Martos C, Ascunce N, et al. The causal link between human papillomavirus and invasive cervical cancer: a population-based case-control study in Colombia and Spain. Int J Cancer. 1992;52(5):743–9.

    Article  CAS  Google Scholar 

  7. Walboomers JM, Jacobs MV, Manos MM, Bosch FX, Kummer JA, Shah KV, Snijders PJ, Peto J, Meijer CJ, Munoz N. Human papillomavirus is a necessary cause of invasive cervical cancer worldwide. J Pathol. 1999;189(1):12–9.

    Article  CAS  Google Scholar 

  8. Cuzick J, Wheeler C. Need for expanded HPV genotyping for cervical screening. Papillomavirus Res. 2016;2:112–5.

    Article  Google Scholar 

  9. Joura EA, Ault KA, Bosch FX, Brown D, Cuzick J, Ferris D, Garland SM, Giuliano AR, Hernandez-Avila M, Huh W, et al. Attribution of 12 high-risk human papillomavirus genotypes to infection and cervical disease. Cancer Epidemiol Biomarkers Prev. 2014;23(10):1997–2008.

    Article  Google Scholar 

  10. Wright TC Jr, Parvu V, Stoler MH, Kodsi S, Eckert K, Yanson K, Cooper CK. HPV infections and cytologic abnormalities in vaccinated women 21–34 years of age: results from the baseline phase of the Onclarity trial. Gynecol Oncol. 2019;153(2):259–65.

    Article  Google Scholar 

  11. Trevisan A, Schlecht NF, Ramanakumar AV, Villa LL, Franco EL, Ludwig-McGill Study G. Human papillomavirus type 16 viral load measurement as a predictor of infection clearance. J Gen Virol. 2013;94(Pt 8):1850–7.

  12. Moberg M, Gustavsson I, Gyllensten U. Type-specific associations of human papillomavirus load with risk of developing cervical carcinoma in situ. Int J Cancer. 2004;112(5):854–9.

    Article  CAS  Google Scholar 

  13. Moberg M, Gustavsson I, Wilander E, Gyllensten U. High viral loads of human papillomavirus predict risk of invasive cervical carcinoma. Br J Cancer. 2005;92(5):891–4.

    Article  CAS  Google Scholar 

  14. Swan DC, Tucker RA, Tortolero-Luna G, Mitchell MF, Wideroff L, Unger ER, Nisenbaum RA, Reeves WC, Icenogle JP. Human papillomavirus (HPV) DNA copy number is dependent on grade of cervical disease and HPV type. J Clin Microbiol. 1999;37(4):1030–4.

    Article  CAS  Google Scholar 

  15. Xi LF, Hughes JP, Edelstein ZR, Kiviat NB, Koutsky LA, Mao C, Ho J, Schiffman M. Human Papillomavirus (HPV) type 16 and type 18 DNA loads at baseline and persistence of type-specific infection during a 2-year follow-up. J Infect Dis. 2009;200(11):1789–97.

    Article  CAS  Google Scholar 

  16. Xi LF, Hughes JP, Castle PE, Edelstein ZR, Wang C, Galloway DA, Koutsky LA, Kiviat NB, Schiffman M. Viral load in the natural history of human papillomavirus type 16 infection: a nested case-control study. J Infect Dis. 2011;203(10):1425–33.

    Article  Google Scholar 

  17. de Abreu AL, Malaguti N, Souza RP, Uchimura NS, Ferreira EC, Pereira MW, Carvalho MD, Pelloso SM, Bonini MG, Gimenes F, et al. Association of human papillomavirus, Neisseria gonorrhoeae and Chlamydia trachomatis co-infections on the risk of high-grade squamous intraepithelial cervical lesion. Am J Cancer Res. 2016;6(6):1371–83.

    PubMed  PubMed Central  Google Scholar 

  18. Verteramo R, Pierangeli A, Mancini E, Calzolari E, Bucci M, Osborn J, Nicosia R, Chiarini F, Antonelli G, Degener AM. Human Papillomaviruses and genital co-infections in gynaecological outpatients. BMC Infect Dis. 2009;9:16.

    Article  Google Scholar 

  19. Biernat-Sudolska M, Szostek S, Rojek-Zakrzewska D, Klimek M, Kosz-Vnenchak M. Concomitant infections with human papillomavirus and various mycoplasma and ureaplasma species in women with abnormal cervical cytology. Adv Med Sci. 2011;56(2):299–303.

    Article  CAS  Google Scholar 

  20. Discacciati MG, Gimenes F, Pennacchi PC, Faiao-Flores F, Zeferino LC, Derchain SM, Teixeira JC, Costa MC, Zonta M, Termini L, et al. MMP-9/RECK Imbalance: a mechanism associated with high-grade cervical lesions and genital infection by human papillomavirus and Chlamydia trachomatis. Cancer Epidemiol Biomarkers Prev. 2015;24(10):1539–47.

    Article  CAS  Google Scholar 

  21. Finan RR, Musharrafieh U, Almawi WY. Detection of Chlamydia trachomatis and herpes simplex virus type 1 or 2 in cervical samples in human papilloma virus (HPV)-positive and HPV-negative women. Clin Microbiol Infect. 2006;12(9):927–30.

    Article  CAS  Google Scholar 

  22. Paba P, Bonifacio D, Di Bonito L, Ombres D, Favalli C, Syrjanen K, Ciotti M. Co-expression of HSV2 and Chlamydia trachomatis in HPV-positive cervical cancer and cervical intraepithelial neoplasia lesions is associated with aberrations in key intracellular pathways. Intervirology. 2008;51(4):230–4.

    Article  CAS  Google Scholar 

  23. Samoff E, Koumans EH, Markowitz LE, Sternberg M, Sawyer MK, Swan D, Papp JR, Black CM, Unger ER. Association of Chlamydia trachomatis with persistence of high-risk types of human papillomavirus in a cohort of female adolescents. Am J Epidemiol. 2005;162(7):668–75.

    Article  Google Scholar 

  24. Smith JS, Bosetti C, Munoz N, Herrero R, Bosch FX, Eluf-Neto J, Meijer CJ, Van Den Brule AJ, Franceschi S, Peeling RW, et al. Chlamydia trachomatis and invasive cervical cancer: a pooled analysis of the IARC multicentric case-control study. Int J Cancer. 2004;111(3):431–9.

    Article  CAS  Google Scholar 

  25. Smith JS, Herrero R, Bosetti C, Munoz N, Bosch FX, Eluf-Neto J, Castellsague X, Meijer CJ, Van den Brule AJ, Franceschi S, et al. Herpes simplex virus-2 as a human papillomavirus cofactor in the etiology of invasive cervical cancer. J Natl Cancer Inst. 2002;94(21):1604–13.

    Article  CAS  Google Scholar 

  26. Zhang ZF, Graham S, Yu SZ, Marshall J, Zielezny M, Chen YX, Sun M, Tang SL, Liao CS, Xu JL, et al. Trichomonas vaginalis and cervical cancer. A prospective study in China. Ann Epidemiol. 1995;5(4):325–32.

    Article  CAS  Google Scholar 

  27. Castle PE, Sadorra M, Lau T, Aldrich C, Garcia FA, Kornegay J. Evaluation of a prototype real-time PCR assay for carcinogenic human papillomavirus (HPV) detection and simultaneous HPV genotype 16 (HPV16) and HPV18 genotyping. J Clin Microbiol. 2009;47(10):3344–7.

    Article  Google Scholar 

  28. Heideman DA, Hesselink AT, Berkhof J, van Kemenade F, Melchers WJ, Daalmeijer NF, Verkuijten M, Meijer CJ, Snijders PJ. Clinical validation of the cobas 4800 HPV test for cervical screening purposes. J Clin Microbiol. 2011;49(11):3983–5.

    Article  CAS  Google Scholar 

  29. FDA: Policy for Coronavirus Disease-2019 Tests During the Public Health Emergency (Revised), U.S. Department of Health and Human Services. Food and Drug Administration (FDA). 2020.

  30. Eklund C, Forslund O, Wallin KL, Dillner J. Global improvement in genotyping of human papillomavirus DNA: the 2011 HPV LabNet International Proficiency Study. J Clin Microbiol. 2014;52(2):449–59.

    Article  Google Scholar 

  31. van den Brule AJ, Pol R, Fransen-Daalmeijer N, Schouls LM, Meijer CJ, Snijders PJ. GP5+/6+ PCR followed by reverse line blot analysis enables rapid and high-throughput identification of human papillomavirus genotypes. J Clin Microbiol. 2002;40(3):779–87.

    Article  Google Scholar 

  32. Meijer CJ, Berkhof J, Castle PE, Hesselink AT, Franco EL, Ronco G, Arbyn M, Bosch FX, Cuzick J, Dillner J, et al. Guidelines for human papillomavirus DNA test requirements for primary cervical cancer screening in women 30 years and older. Int J Cancer. 2009;124(3):516–20.

    Article  CAS  Google Scholar 

  33. Wallin KL, Wiklund F, Luostarinen T, Angstrom T, Anttila T, Bergman F, Hallmans G, Ikaheimo I, Koskela P, Lehtinen M, et al. A population-based prospective study of Chlamydia trachomatis infection and cervical carcinoma. Int J Cancer. 2002;101(4):371–4.

    Article  CAS  Google Scholar 

  34. Horner P, Donders G, Cusini M, Gomberg M, Jensen JS, Unemo M. Should we be testing for urogenital Mycoplasma hominis, Ureaplasma parvum and Ureaplasma urealyticum in men and women? A position statement from the European STI Guidelines Editorial Board. J Eur Acad Dermatol Venereol. 2018;32(11):1845–51.

    Article  CAS  Google Scholar 

  35. Bottari F, Iacobone AD, Passerini R, Preti EP, Sandri MT, Cocuzza CE, Gary DS, Andrews JC. Human papillomavirus genotyping compared with a qualitative high-risk human papillomavirus test after treatment of high-grade cervical intraepithelial neoplasia: a systematic review. Obstet Gynecol. 2019;134(3):452–62.

    Article  CAS  Google Scholar 

  36. Demarco M, Hyun N, Carter-Pokras O, Raine-Bennett TR, Cheung L, Chen X, Hammer A, Campos N, Kinney W, Gage JC, et al. A study of type-specific HPV natural history and implications for contemporary cervical cancer screening programs. EClinicalMedicine. 2020;22:100293.

    Article  Google Scholar 

  37. Jentschke M, Soergel P, Hillemanns P. Importance of HPV genotyping for the screening, therapy and management of cervical neoplasias. Geburtshilfe Frauenheilkd. 2012;72(6):507–12.

    Article  CAS  Google Scholar 

  38. Ardhaoui M, Ennaifer E, De Matos Salim AC, Gomez FM, Laasili T, Boubaker MS, Guizani I. Nested PCR followed by NGS: Validation and application for HPV genotyping of Tunisian cervical samples. PLoS ONE. 2021;16(8):e0255914.

    Article  CAS  Google Scholar 

  39. Arroyo Muhr LS, Lagheden C, Lei J, Eklund C, Nordqvist Kleppe S, Sparen P, Sundstrom K, Dillner J. Deep sequencing detects human papillomavirus (HPV) in cervical cancers negative for HPV by PCR. Br J Cancer. 2020;123(12):1790–5.

    Article  CAS  Google Scholar 

  40. Wagner S, Roberson D, Boland J, Yeager M, Cullen M, Mirabello L, Dunn ST, Walker J, Zuna R, Schiffman M, et al. Development of the TypeSeq assay for detection of 51 human papillomavirus genotypes by next-generation sequencing. J Clin Microbiol. 2019;57(5):e01794-18.

    Article  Google Scholar 

Download references


Not applicable.

Author information

Authors and Affiliations



ZM, BG, and CW are involved in the overall study of this assay and manuscript writing. XC and ML participated in the LoD assay and manuscript writing. The INEI-ANLIS Malbrán investigators are responsible for the selection of clinical samples previously characterized, participate in the analysis, interpretation of the data and revision of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Chunlin Wang.

Ethics declarations

Ethics approval and consent to participate

This cohort is part of a screening study that is being carried out by the INEI-ANLIS Malbrán with the collaboration of the “Prof. A. Posadas” Hospital (Argentina); it was approved by the Research Ethics Board and the Ethics Committee of this hospital. Written informed consent was obtained before enrolling the participant women, and confidentiality and anonymity were guaranteed.

Consent for publication

Not applicable.

Competing interests

ZM, BG, and CW are shareholders and employees of Chapter Diagnostics. A patent has been issued for the technology. The INEI-ANLIS Malbrán investigators declare no competing interests; they did not receive any financial contribution from Chapter Diagnostics.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Comparison of Illumina sequencing platforms and sequencing kits for HPV-STI analysis.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ma, Z., Gharizadeh, B., Cai, X. et al. A comprehensive HPV-STI NGS assay for detection of 29 HPV types and 14 non-HPV sexually transmitted infections. Infect Agents Cancer 17, 9 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: