# Resequencing Microarray Technology for Genotyping Human Papillomavirus in Cervical Smears

Nicolas Berthet, Michael Falguières, Claudia Filippone, Chloé Bertolus, Christine Bole-Feysot, Sylvain Brisse, Antoine Gessain, Isabelle Heard, Michel Favre

► To cite this version:

Nicolas Berthet, Michael Falguières, Claudia Filippone, Chloé Bertolus, Christine Bole-Feysot, et al.. Resequencing Microarray Technology for Genotyping Human Papillomavirus in Cervical Smears. PLoS ONE, 2014, 9 (11), pp.e109301. 10.1371/journal.pone.0109301. hal-01344730

HAL Id: hal-01344730

https://hal.sorbonne-universite.fr/hal-01344730v1

Submitted on 12 Jul 2016

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Resequencing Microarray Technology for Genotyping Human Papillomavirus in Cervical Smears

Nicolas Berthet¹,²*¤, Michael Falguières³, Claudia Filippone¹,², Chloé Bertolus⁴,⁵, Christine Bole-Feysot⁶, Sylvain Brisse⁷, Antoine Gessain¹,², Isabelle Heard³,⁸, Michel Favre³,⁹

¹Institut Pasteur, Epidémiologie et Physiopathologie des Virus Oncogènes, 75724 Paris Cedex 15, France, ²Centre National de la Recherche Scientifique, UMR 3569, 75724 Paris Cedex 15, France, ³Institut Pasteur, Centre National de Référence des Papillomavirus, 75724 Paris Cedex 15, France, ⁴AP-HP, Hôpital Pitié-Salpêtrière, Service de Chirurgie Maxillo-Faciale et Stomatologie, 75013 Paris, France, ⁵UPMC, Université Paris 06, CIMI-Paris, UMR5 CR7, INSERM U1135, 75005, Paris, France, ⁶Imagine, Institut des maladies génétiques - Plateforme Génomique, Hôpital Necker - Enfants Malades, 75743 Paris cedex 15, France, ⁷Institut Pasteur, Plate-forme Génotypage des Pathogènes et Santé Publique, 75724 Paris Cedex 15, France, ⁸UPMC Université Paris 06, Groupe hospitalier Pitié-Salpêtrière, Paris Cedex 13, France, ⁹Institut Pasteur, Unité de Génétique, Papillomavirus et Cancer humain, 75724 Paris Cedex 15, France

<figure><img src="image_2.png" /></figure>

## Abstract

There are more than 40 human papillomaviruses (HPVs) belonging to the alpha genus that cause sexually transmitted infections; these infections are among the most frequent and can lead to condylomas and anogenital intra-epithelial neoplasia. At least 18 of these viruses are causative agents of anogenital carcinomas. We evaluated the performance of a resequencing microarray for the detection and genotyping of alpha HPV of clinical significance using cloned HPV DNA. To reduce the number of HPV genotypes tiled on microarray, we used reconstructed ancestral sequences (RASs) as they are more closely related to the various genotypes than the current genotypes are among themselves. The performance of this approach was tested by genotyping with a set of 40 cervical smears already genotyped using the commercial PapilloCheck kit. The results of the two tests were concordant for 70% (28/40) of the samples and compatible for 30% (12/40). Our findings indicate that RASs were able to detect and identify one or several HPV in clinical samples. Associating RASs with homonym sequences improved the genotyping of HPV present in cases of multiple infection. In conclusion, we demonstrate the diagnostic potential of resequencing technology for genotyping of HPV, and illustrate its value both for epidemiological studies and for monitoring the distribution of HPV in the post-vaccination era.

**Citation:** Berthet N, Falguières M, Filippone C, Bertolus C, Bole-Feysot C, et al. (2014) Resequencing Microarray Technology for Genotyping Human Papillomavirus in Cervical Smears. PLoS ONE 9(11): e109301. doi:10.1371/journal.pone.0109301

**Editor:** Richard C. Willson, University of Houston, United States of America

**Received** November 21, 2013; **Accepted** September 10, 2014; **Published** November 10, 2014

**Copyright:** © 2014 Berthet et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This program was supported by the Programme Transversal de Recherche DEVA No. 246 financed by the Institut Pasteur (Paris, France), sponsored by the Fondation Total-Institut Pasteur and the Elisabeth Taub award from the Académie Nationale de Médécine. This study was funded by the French Institut National du Cancer (Canceropole 2011-1VADS 03-IP1). The funders had no role in study design, data analysis or preparation of the manuscript.

**Competing Interests:** The authors have declared that no competing interests exist.

* Email: nicolas.berthet@ird.fr

¤ Current address: CIRMF, Département Zoonose et Maladies Emergentes, BP769 Franceville, Gabon

## Introduction

More than 170 human papillomavirus (HPVs) have been characterized and they are classified into the alpha, beta, gamma, mu, and nu genera [1]. Some alpha papillomaviruses presenting a mucous tropism cause some of the most frequent sexually transmitted infections: around 300 million women worldwide carry a HPV infection of the uterine cervix. About 40 HPV genotypes able to infect genital mucous membranes have been identified. Low-risk mucosal HPVs (LR HPVs) are responsible for condylomas found in 1–2% of the sexually active population. High-risk HPVs (HR HPVs) are causative agents of more than 80% of high-grade cervical intra-epithelial neoplasia [2]. Persistent lesions associated with HR HPVs may evolve towards invasive cervical carcinomas. HPV DNA sequences are detected in the vast majority of adenocarcinomas, adenosquamous carcinomas and squamous cell carcinomas of the cervix, all of which are preceded by premalignant lesions [3]. HR HPVs also contributes to the pathogenesis of other anogenital and aero-digestive tract cancers [4].

The need for effective screening for cervical cancer is now recognized worldwide and testing for HPVs is a potentially powerful approach to the detection of cervical or anogenital lesions. Other than the hybrid capture II and the Cervista detection tests [5], most of the available tests for HPVs are based on PCR methods using degenerate or consensus primers. All these methods allow the detection, and some allow genotyping, of a wide spectrum of HPVs. In addition to their clinical value, robust methods of detection are required for epidemiological studies of HPV infections. In particular, they could be used to study the prevalence and distribution of circulating HPV genotypes in vaccinated and non-vaccinated populations.

Several low or high-density microarrays have been developed and successfully used for the detection of several viral families such as *Papillomaviridae*. Most of them are based on long oligonucleotide probes and virus identification is based only on hybridization patterns. This approach leads to a lack of information for pathogen identification at the single-pair resolution level [6,7]. However, this limitation has since been overcome with the use of the resequencing microarrays (RMAs) method. Resequencing
microarray (RMA) technology involves generating a DNA/RNA sequence corresponding to the pathogen(s) present in a sample in a single assay. It allows detection of a sequence with a tolerance of divergence of up to 10–15% from those tiled on the microarray [8,9,10]. We have previously used RMAs for the detection of a large panel of viral and bacterial pathogens by resequencing the hybridised DNA. This approach has been successfully applied to diverse clinical situations such as the establishment of a rapid laboratory diagnostic test for pandemic influenza viruses [11,12], the characterization of different strains of the monkeypox virus [13,14], the genotyping of several members of the Rhabdoviridae family [8] and the identification and characterization of arboviruses and hemorrhagic fever viruses in biological samples [15,16].

The aim of this study was to develop a multiplex bioassay method, based on RMA technology, for the detection of mucosal HPV genotypes of public health importance. Because of the wide genetic diversity of HPVs and to reduce the number of genotypes tiled on the microarray, we inferred reconstructed ancestral sequences (RASs); these sequences are more closely related to different genotypes than such genotypes are to each other. The advantages of using reconstructed ancestral sequences for RMAs has been described and illustrated using pathogen identification with *Enterobacteriaceae rpoB* sequences used as models [17]. We went on to apply the RMA method developed to detection of HPV DNA sequences in clinical samples.

<figure><img src="image_3.png" /><figcaption>**Figure 1. Phylogenetic tree of 38 alpha HPV isolates, constructed to assess the resequencing microarray approach.** The tree was constructed from the 416-base sequences of part of the L1 ORF in the different HPVs. The tree was generated by the neighbor-joining method using BioNumerics v5.10 software (Applied-Maths, Belgium). HPV sequences tiled on RMA (red squares), location of RASs on the tree (blue circles) and cloned HPVs tested with RMAs (green triangles) are indicated.
doi:10.1371/journal.pone.0109301.g001</figcaption></figure>
# Materials and Methods

## Principle and content of resequencing microarray

A third generation resequencing microarray, named VirID v3.0, containing 724 viral sequences distributed into 27 families and 85 genera (Table S1), was used in this study. Among the many viral sequences tiled, a set of 38 different HPV sequences was used. The size of the tiled region ranged from 455 to 482 nucleotides corresponding to part of the L1 open reading frame. This set includes sequences in the following genotypes: HPV6, 11, 16, 18, 26, 30, 31, 33, 34, 35, 40, 42, 43, 44, 45, 51, 53, 54, 56, 58, 59, 65, 66, 67, 68, 69, 70, 73, 74, 82, 85, 88, 91, 95, 97, 101, 103 and 108. In addition, ten reconstructed ancestral sequences (RASs) were tiled on the microarray. These sequences are described in Tables S1 and S2. For each sequence tiled, the principle of RMA is to interrogate each single base of the unknown sequence to be detected with a set of eight appropriate 25-mer probes. Two probes among the eight (four for each sense of the selected sequence) correspond to perfect matches at the central (13th) position of the probe, and the other six probes represent all possible mismatches at the same position.

## Ancestral sequences reconstruction

RASs were designed based on phylogenetic analysis of nucleotide sequences of a portion of 416 bases in the L1 ORF of 38 alpha-HPV; there were neither insertions nor deletions in this stretch. A neighbor-joining tree was obtained using software BioNumerics v5.10 (Applied-Maths, Belgium). Ancestral sequences were reconstructed by maximum likelihood using the PAML v4 software [18]. The nucleotide substitution model used was K80 with gamma (number of categories of distinct substitution rates) and kappa (transition/transversion ratio) parameters estimated.

## Extraction and amplification of HPV DNA

Viral DNA was obtained from cloned HPV DNA and from cervical smears provided by the National Reference Laboratory for Human Papillomavirus at Institut Pasteur (Paris). The cloned HPV DNA used for validation assays was amplified and purified as previously described [19]. Total DNA was extracted from cervical smears using the NucleoSpin Tissue extraction kit according to the

recommendations of the supplier (Macherey Nagel, Hoerdt, France). Following extraction, samples were digested with Plasmid-Safe DNase for 12 hours at 37°C according to the manufacturer's instructions to enrich for circular viral DNA (Epicentre). All DNAs were amplified using the Repli-g Mini Kit according to the manufacturer's instructions (Qiagen). The concentration of amplified DNAs was determined with the Quant-it kit (Life technologies) and aliquots of 10 ug were used for hybridization in RMA assays.

## Ethical considerations

Liquid-based cytology (LBC) samples were collected from women attending organized cervical screening in 16 pilot sites in France. After completion of cytology, residual LBC samples were anonymized and sent to the French HPV Reference Laboratory for genotyping. The residual material would otherwise have been discarded. According to French regulations of biomedical research, a written or verbal informed consent and an ethical approval are not required for such studies. However, the study was given favorable ethical opinion by the Comité de Recherche Clinique at the Institut Pasteur (Avis n° 200-51). The collection of human samples was declared to the French Research Ministry (n°DC-2010-1197, collection n°1) in accordance with French regulations.

## Hybridization on RMA and data analysis

Amplified DNA was fragmented and labeled using the GeneChip Resequencing Assay Kit (Affymetrix Inc., Santa Clara, CA). After overnight hybridization at 45°C, the RMA was washed, stained and scanned according to Affymetrix' instructions. The raw image file (.DAT), obtained after the scan, was converted to a fluorescence intensity file (.CEL). Bases were called by the GeneChip Sequence Analysis Software (GSEQ v 4.1) which uses a derivative of the ABACUS base-calling algorithm [20]. This algorithm consists of an automated statistical method that analyses raw RMA hybridization data and optimizes the base-calling process. For each base-call, a quality score is calculated from the difference between the best fitting model and the second best one. If this score is below a chosen threshold, an undetermined base "N" is assigned to that position. Sequences were output in FASTA

<table><caption>**Table 1.** Resequencing microarray assays with cloned HPV DNA.</caption><thead><tr><td>HPV species</td><td>HPV genotype</td><td>Call rate<sup>a</sup> (%)</td><td>Accuracy<sup>b</sup> (%)</td><td>BLAST score<sup>c</sup></td></tr></thead><tbody><tr><td>α1</td><td>42</td><td>96.5</td><td>99.5</td><td>780</td></tr><tr><td>α5</td><td>51</td><td>94.4</td><td>99.3</td><td>749</td></tr><tr><td>α6</td><td>53</td><td>93.6</td><td>98.6</td><td>717</td></tr><tr><td></td><td>56</td><td>95.7</td><td>99.5</td><td>751</td></tr><tr><td></td><td>66</td><td>95.4</td><td>100</td><td>751</td></tr><tr><td>α7</td><td>18</td><td>94.5</td><td>99.5</td><td>765</td></tr><tr><td></td><td>45</td><td>94.2</td><td>99.8</td><td>764</td></tr><tr><td>α9</td><td>16</td><td>94.9</td><td>100</td><td>767</td></tr><tr><td>α10</td><td>6</td><td>94.7</td><td>100</td><td>758</td></tr><tr><td></td><td>11</td><td>94.4</td><td>100</td><td>738</td></tr><tr><td>α13</td><td>54</td><td>92.1</td><td>100</td><td>704</td></tr><tr><td>γ</td><td>108</td><td>97.7</td><td>100</td><td>809</td></tr></tbody></table>

<footnote><sup>a</sup>Ratio between the number of bases determined and the sequence length.</footnote>

<footnote><sup>b</sup>Ratio between the number of correctly determined bases and the total number of determined bases.</footnote>

<footnote><sup>c</sup>The score of the highest scoring high-scoring segment pair (HSP) from that database sequence.</footnote>

<footnote>doi:10.1371/journal.pone.0109301.t001</footnote>
<table><caption>**Table 2.** Identification of HPV16 and other HPVs with homonymous and non homonymous sequences.</caption><thead><tr><th rowspan="3">Species</th><th colspan="10">Cloned HPV tested</th></tr><tr><th colspan="2">16</th><th colspan="2">66</th><th colspan="2">16+66</th><th colspan="2">16+66+51+108</th><th rowspan="2">HPV identified</th></tr><tr><th>Tiled HPV</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>HPV identified</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>HPV identified</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>Call rate value<sup>a</sup></th></tr></thead><tbody><tr><td>α5</td><td>51</td><td>15.2</td><td>Neg</td><td>N.A.</td><td>12.8</td><td>Neg</td><td>N.A.</td><td>22.5</td><td>Neg</td><td>N.A.</td><td>94.4</td><td>749</td><td>51</td></tr><tr><td>α6</td><td>82</td><td>18.7</td><td>Neg</td><td>N.A.</td><td>13.0</td><td>53.6</td><td>66</td><td>21.7</td><td>Neg</td><td>N.A.</td><td>45.3</td><td>147</td><td>51</td></tr><tr><td>α6</td><td>56</td><td>24.1</td><td>Neg</td><td>N.A.</td><td>41.4</td><td>134</td><td>66</td><td>46.6</td><td>75.2</td><td>66</td><td>46.0</td><td>118</td><td>66</td></tr><tr><td>α6</td><td>53</td><td>17.1</td><td>Neg</td><td>N.A.</td><td>31.9</td><td>59</td><td>66</td><td>34.7</td><td>Neg</td><td>N.A.</td><td>38.5</td><td>69.8</td><td>66</td></tr><tr><td>α6</td><td>66</td><td>28.4</td><td>Neg</td><td>N.A.</td><td>95.4</td><td>751</td><td>66</td><td>91.3</td><td>677</td><td>66</td><td>95.4</td><td>751</td><td>66</td></tr><tr><td>α7</td><td>30</td><td>26.9</td><td>51.8</td><td>16</td><td>29.9</td><td>101</td><td>66</td><td>41.9</td><td>145</td><td>66</td><td>36.9</td><td>205</td><td>66</td></tr><tr><td>α7</td><td>45</td><td>21.6</td><td>53.6</td><td>16</td><td>12.9</td><td>Neg</td><td>N.A.</td><td>25.3</td><td>59</td><td>Und.</td><td>22.4</td><td>60.8</td><td>Und.</td></tr><tr><td>α9</td><td>59</td><td>23.9</td><td>44.6</td><td>16</td><td>10.3</td><td>Neg</td><td>N.A.</td><td>25.6</td><td>Neg</td><td>N.A.</td><td>26.4</td><td>41</td><td>16</td></tr><tr><td>α9</td><td>16</td><td>94.9</td><td>767</td><td>16</td><td>23.1</td><td>77</td><td>66</td><td>90.3</td><td>688</td><td>16</td><td>94.9</td><td>767</td><td>16</td></tr><tr><td>α10</td><td>33</td><td>25.2</td><td>50</td><td>16</td><td>20.7</td><td>44.6</td><td>66</td><td>32.4</td><td>44.6</td><td>16</td><td>32.6</td><td>94.5</td><td>16</td></tr><tr><td>α10</td><td>35</td><td>29.6</td><td>Neg</td><td>N.A.</td><td>14.8</td><td>59</td><td>66</td><td>33.8</td><td>101</td><td>16</td><td>30.8</td><td>64.4</td><td>16</td></tr><tr><td>α13</td><td>74</td><td>22.8</td><td>46.4</td><td>16</td><td>9.8</td><td>Neg</td><td>N.A.</td><td>22.9</td><td>44.6</td><td>16</td><td>22.6</td><td>48.2</td><td>16</td></tr><tr><td>α13</td><td>54</td><td>23.3</td><td>44.6</td><td>16</td><td>12.7</td><td>Neg</td><td>N.A.</td><td>24.4</td><td>46.4</td><td>Und.</td><td>24.6</td><td>50</td><td>Und.</td></tr><tr><td>α16</td><td>108</td><td>18.1</td><td>Neg</td><td>N.A.</td><td>14.1</td><td>Neg</td><td>N.A.</td><td>18.9</td><td>Neg</td><td>N.A.</td><td>97.7</td><td>809</td><td>108</td></tr></tbody></table>

<footnote><sup>a</sup>Ratio between the number of bases determined and the sequence length.</footnote>

<footnote>Values in bold correspond to the call rate value, NBlast score and result of BLAST of HPVs tested with the homonymous sequence. Neg: Negative; N.A: Not applicable; Und.: Undetermined.
doi:10.1371/journal.pone.0109301.t002</footnote>
<table><caption>**Table 3.** Identification of HPV18 and other HPVs with homonymous and non homonymous sequences.</caption><thead><tr><th rowspan="3">Species</th><th colspan="8">Cloned HPV</th><th colspan="2">18+53</th><th colspan="2">18+53+6+42</th></tr><tr><th colspan="2">18</th><th colspan="2">18+45</th><th colspan="2">18+53</th><th colspan="2">18+53</th><th colspan="2">18+53+6+42</th></tr><tr><th>Tiled HPV</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>HPV identified</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>HPV identified</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>HPV identified</th><th>Call rate value<sup>a</sup></th><th>NBlast score</th><th>HPV identified</th></tr></thead><tbody><tr><td>α1</td><td>42</td><td>11.2</td><td>Neg</td><td>N.A.</td><td>12.5</td><td>Neg</td><td>N.A.</td><td>23.5</td><td>Neg</td><td>N.A.</td><td>96.6</td><td>780</td><td>42</td></tr><tr><td>α6</td><td>53</td><td>4</td><td>Neg</td><td>N.A.</td><td>9.5</td><td>Neg</td><td>N.A.</td><td><strong>92.7</strong></td><td><strong>702</strong></td><td><strong>53</strong></td><td><strong>93.6</strong></td><td><strong>717</strong></td><td><strong>53</strong></td></tr><tr><td>α7</td><td>56</td><td>10.6</td><td>Neg</td><td>N.A.</td><td>18.0</td><td>Neg</td><td>N.A.</td><td>28.0</td><td>46.4</td><td>53</td><td>29.5</td><td>41</td><td>53</td></tr><tr><td></td><td>18</td><td><strong>96.6</strong></td><td><strong>801</strong></td><td><strong>18</strong></td><td><strong>93.7</strong></td><td><strong>751</strong></td><td><strong>18</strong></td><td><strong>89.1</strong></td><td><strong>666</strong></td><td><strong>18</strong></td><td><strong>94.5</strong></td><td><strong>765</strong></td><td><strong>18</strong></td></tr><tr><td></td><td>45</td><td>47.9</td><td>168</td><td>18</td><td><strong>94.2</strong></td><td><strong>764</strong></td><td><strong>45</strong></td><td>31.3</td><td>136</td><td>18</td><td>52.9</td><td>204</td><td>18</td></tr><tr><td></td><td>59</td><td>20.2</td><td>Neg</td><td>N.A.</td><td>27.5</td><td>46.4</td><td>Und.</td><td>17.0</td><td>Neg</td><td>N.A.</td><td>26.9</td><td>Neg</td><td>N.A.</td></tr><tr><td></td><td>68</td><td>20.8</td><td>75.2</td><td>18</td><td>28.3</td><td>123</td><td>18</td><td>22.7</td><td>57.2</td><td>18</td><td>27.8</td><td>128</td><td>18</td></tr><tr><td></td><td>70</td><td>21.4</td><td>Neg</td><td>N.A.</td><td>31.9</td><td>62.8</td><td>45</td><td>23.1</td><td>51.8</td><td>53</td><td>31.5</td><td>Neg</td><td>N.A.</td></tr><tr><td>α8</td><td>40</td><td>17.2</td><td>Neg</td><td>N.A.</td><td>26.0</td><td>Neg</td><td>N.A.</td><td>14.3</td><td>Neg</td><td>N.A.</td><td>20.0</td><td>Neg</td><td>N.A.</td></tr><tr><td>α9</td><td>43</td><td>11.5</td><td>Neg</td><td>N.A.</td><td>27.7</td><td>55.4</td><td>Und.</td><td>13.2</td><td>Neg</td><td>N.A.</td><td>23.0</td><td>42.8</td><td>Und.</td></tr><tr><td></td><td>16</td><td>17.1</td><td>46.4</td><td>18</td><td>21.8</td><td>51.8</td><td>18</td><td>16.9</td><td>Neg</td><td>N.A.</td><td>32.6</td><td>77</td><td>6</td></tr><tr><td></td><td>31</td><td>15.0</td><td>69.8</td><td>18</td><td>23.3</td><td>64.4</td><td>Und</td><td>19.9</td><td>57.2</td><td>53</td><td>28.8</td><td>122</td><td>Und.</td></tr><tr><td></td><td>33</td><td>14.5</td><td>Neg</td><td>N.A.</td><td>21.7</td><td>50</td><td>18</td><td>19.4</td><td>42.8</td><td>53</td><td>26.4</td><td>66.2</td><td>6</td></tr><tr><td></td><td>35</td><td>16.5</td><td>Neg</td><td>N.A.</td><td>25.3</td><td>46.4</td><td>18</td><td>17.3</td><td>Neg</td><td>N.A.</td><td>38.3</td><td>48.2</td><td>42</td></tr><tr><td></td><td>58</td><td>10.6</td><td>Neg</td><td>N.A.</td><td>16.6</td><td>Neg</td><td>N.A.</td><td>13.2</td><td>Neg</td><td>N.A.</td><td>20.2</td><td>Neg</td><td>N.A.</td></tr><tr><td></td><td>67</td><td>23.0</td><td>62.6</td><td>18</td><td>26.4</td><td>112</td><td>18</td><td>18.9</td><td>Neg</td><td>N.A.</td><td>32.1</td><td>98.2</td><td>18</td></tr><tr><td>α10</td><td>6</td><td>17.4</td><td>54.8</td><td>18</td><td>21.3</td><td>60.8</td><td>Und.</td><td>23.0</td><td>Neg</td><td>N.A.</td><td><strong>94.7</strong></td><td><strong>758</strong></td><td><strong>6</strong></td></tr><tr><td>α11</td><td>73</td><td>15.9</td><td>Neg</td><td>N.A.</td><td>19.2</td><td>57.2</td><td>18</td><td>13.4</td><td>Neg</td><td>N.A.</td><td>29.2</td><td>42.8</td><td>18</td></tr><tr><td>α13</td><td>54</td><td>21.1</td><td>Neg</td><td>N.A.</td><td>28.2</td><td>103</td><td>45</td><td>16.8</td><td>Neg</td><td>N.A.</td><td>34.1</td><td>Neg</td><td>N.A.</td></tr></tbody></table>

<footnote><sup>a</sup>Ratio between the number of bases determined and the sequence length.</footnote>

<footnote>Values in bold correspond to the call rate value, NBlast score and result of BLAST of HPVs tested with the homonymous sequence. Neg: Negative; N.A: Not applicable; Und.: Undetermined.</footnote>

<footnote>doi:10.1371/journal.pone.0109301.t003</footnote>
<table><caption>**Table 4.** Identification of HPVs belonging to the α6 family based on homonymous and reconstructed ancestral sequences.</caption><thead><tr><th rowspan="2">HPV sequence tiled on RMA</th><th rowspan="2"></th><th colspan="3">HPV tested</th></tr><tr><th>53</th><th>56</th><th>66</th></tr></thead><tbody><tr><td rowspan="2">HPV53</td><td>Call rate<sup>a</sup></td><td><strong>93.6<sup>a</sup></strong></td><td>35.6</td><td>38.5</td></tr><tr><td>BLAST score</td><td>717</td><td>111</td><td>69.8</td></tr><tr><td rowspan="2">HPV56</td><td>Call rate<sup>a</sup></td><td>29.5</td><td><strong>95.7</strong></td><td>46</td></tr><tr><td>BLAST score</td><td>41</td><td>751</td><td>118</td></tr><tr><td rowspan="2">HPV66</td><td>Call rate<sup>a</sup></td><td>34.9</td><td>43.9</td><td><strong>95.4</strong></td></tr><tr><td>BLAST score</td><td>42.8</td><td>145</td><td>751</td></tr><tr><td rowspan="2">RAS_811</td><td>Call rate<sup>a</sup></td><td>59.1</td><td>71.2</td><td>75.7</td></tr><tr><td>BLAST score</td><td>107</td><td>197</td><td>306</td></tr></tbody></table>

format; for each HPV sequence obtained, the call rate value was calculated as the ratio between the number of bases determined (A or T or C or G) and the entire sequence length. The accuracy of the RMA process was determined as the ratio between the number of correctly determined bases and the total number of determined bases, by comparison with the known sequence of the strains tested.

## Real-time PCR for HPV16 and HPV18

Quantitative PCR tests for HPV16 and HPV18 were performed using DNA obtained from cervical smears (5 ul) or amplified products (1/1000 dilution). The DNA was mixed with 100 pmol of each primer (HPV16-for 5' GAACCGAAACCGGTTAGTATAA 3'; HPV16-rev 5' ATGTATAGTTGTTTGCAGCTCTGT 3'; HPV18-for 5' GGACCGAAAACGGTGTATATAA 3'; HPV18-rev 5' CAGTGAAGTGTTTCAGTTCGGT 3') and 15 pmol of a probe (HPV16-probe 5' Fam-CATTTTATGCACCAAAAGAGAACTGCAATGTTTC-BHQ13'; HPV18-probe 5' Tamra ATGTGAGAAACACACCACAATACTATGGCGCG BHQ2 3') in a final reaction volume of 20 µl, using a LC Taqman Master kit (Roche Diagnostic, Inc.). The LightCycler qPCR assay was performed as described by Schmitz et al. [21].

## PapilloCheck HPV genotyping test

Clinical samples were tested using the PapilloCheck HPV genotyping kit (Greiner BioOne, Frickenhausen, Germany). This test involves PCR amplification with fluorescent primers (CY5 fluorophore) specific for the HPV E1 gene and the cellular housekeeping gene ADAT1. The amplification products are then hybridized to a DNA chip containing sequences for 18 high-risk genotypes (HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 70, 73, 82) and six low-risk genotypes (HPV 6, 11, 40, 42, 43 and 54/55). The PapilloCheck test was used according to the recommendations of the supplier.

## Results

### 1. Validation of RMA using cloned HPV DNA

The third generation of the “VirID” RMA (v3.0) was used for the detection and identification of a broad panel of HPV sequences. Thirty-eight HPV sequences were selected and tiled on the microarray, including 32 from the α genus (HPV6, 11, 16, 18, 26, 30, 31, 33, 34, 35, 40, 42, 43, 44, 45, 51, 53, 54, 56, 58, 59, 66, 67, 68, 69, 70, 73, 74, 82, 85, 91, and 97) and covering all the nine families found in cervical lesions and condylomas. The other

six sequences belonged to the γ genus (HPV65, 88, 95, 101, 103, and 108). Reconstructed Ancestral Sequences (RAS) were designed through phylogenetic analysis of selected α genus sequences. Sequences at nodes at different phylogenetic depths on the tree were chosen as RASs (Fig. 1); we checked that all the RASs selected branched directly or very closely to their respective node, with near-zero branch lengths.

Of the 38 HPV sequences deposited on the chip, 11 α-HPV sequences and one γ-HPV sequence were tested using cloned HPV genomes (Table 1): all the 12 sequences were detected with call rate values between 92.1 and 97.7%, with a resequencing accuracy of 98.6 to 100%. Confirmation of the HPV genotype was obtained following BLASTN analysis, with scores between 704 (for HPV54) to 809 (for HPV108) (Table 1).

### 2. Genotyping of HPVs with homonymous sequences

Because several sequences were tiled on the chip, we tested whether each HPV genotype could be detected following hybridization either to its homonymous sequence, or to related sequences belonging to the same or different species. Various cloned HPVs were tested, alone or in mixtures with HPVs belonging to the same or different species (Tables 2 and 3). HPV16 and HPV66 were perfectly identified by their homonymous sequences, whether they were on their own or in mixtures; they were both identified by non-homonymous sequences, including HPV30 and HPV33. HPV66 was also detected by the sequences of HPV53, HPV56 and HPV82, whereas HPV16 was detected by the sequences of HPV45 and HPV59 (Table 2). With a mixture of both HPV16 and HPV66, only the most closely related HPV was detected by a non-homonymous sequence (Table 2). HPV18 was tested alone and in mixtures with HPV45 or HPV53: the three viruses were perfectly detected by their corresponding homonymous sequences. Moreover, they were simultaneously identified by several other non-homonymous sequences (Table 3): for example, when two HPVs belonging to the same species (such as HPV18 and HPV45) were tested together, no competition was observed for their identification either by homonymous or related sequences. Mixtures containing HPV16, 51, 66, and 108 (Table 2) or HPV6, 18, 42, and 53 (Table 3) were similarly studied, and in all cases the various HPVs were distinguished and identified by both homonymous and non-homonymous sequences.

<footnote><sup>a</sup>The call rate value for HPVs tested with its homonymous sequence is indicated in boldface.
doi:10.1371/journal.pone.0109301.t004</footnote>
<table><caption>**Table 5.** Identification of HPVs by reconstructed ancestral sequences.</caption><thead><tr><th rowspan="2">Mixture</th><th colspan="10">Reconstructed ancestral sequences tiled on the microarray<sup>a</sup></th></tr><tr><th>HPV Genotype (species)</th><th>805 (α 8)</th><th>806 (α 1)</th><th>808 (α 10)</th><th>809 (α 7)</th><th>810 (α 7)</th><th>811 (α 6)</th><th>812 (α 5)</th><th>813 (α 5)</th><th>814 (α 9)</th><th>815 (α 9)</th></tr></thead><tbody><tr><td rowspan="4">A</td><td>16 (α9)</td><td rowspan="4">43.3<sup>b</sup></td><td rowspan="4"></td><td rowspan="4"></td><td rowspan="4"></td><td rowspan="4"></td><td>59.9<sup>b</sup></td><td rowspan="4">88.9<sup>b</sup></td><td>57<sup>b</sup></td><td>39.9<sup>b</sup></td></tr><tr><td>51 (α 5)</td><td rowspan="3">75.7<sup>b</sup></td></tr><tr><td>66 (α 6)</td></tr><tr><td>108 (γ 6)</td></tr><tr><td rowspan="5">B</td><td>6 (α 10)</td><td rowspan="5">54.3<sup>b</sup></td><td rowspan="5"></td><td rowspan="5"></td><td rowspan="5"></td><td rowspan="5"></td><td rowspan="5">54.1<sup>b</sup></td><td rowspan="5">79.6<sup>b</sup></td><td>48.6<sup>b</sup></td><td>45.7<sup>b</sup></td></tr><tr><td>18 (α 7)</td></tr><tr><td>42 (α 1)</td></tr><tr><td>53 (α 6)</td></tr><tr><td>51.7<sup>b</sup></td><td>99.1<sup>b</sup></td></tr></tbody></table>

### 3. HPV genotyping with reconstructed ancestral sequences

We wanted to determine whether all HPVs belonging to a single species could be detected with a RAS specific for this species. Within any family, the RAS is more closely related to the sequences of the various HPVs than the sequences of the HPVs are to each other (data not shown). We tested HPV53, 56 and 66, all members of species alpha 6, and call rates of 59.1 to 75.7 were obtained with RAS-811 (Table 4). These values were much higher than those observed with non-homonymous sequences: 29.5 to 46.0. BLAST analyses of all the sequences obtained confirmed the identification of the genotypes in all cases.

We next investigated whether HPVs that belong to different species could be identified with RASs. We tested a mixture of HPV16, 51, 66 and 108 (mixture A) and one of HPV6, 18, 42 and 53 (mixture B) with reconstructed RASs tiled on the microarray (Table 5). HPV16 belonging to the alpha 9 species was detected by RASs of the same alpha 9 species (RAS 814 and 815) and by RASs designed for the alpha 1 (RAS 806) or alpha 5 (RAS 812) species. By contrast, the alpha 5 HPV51 and the alpha 6 HPV66 were only detected with RASs corresponding to the same species (RAS 813 and RAS 811, respectively). Similarly, HPV18 that belongs to the alpha 7 species was detected with two RASs of the same species (RAS809 and 810). HPV6 that belongs to the alpha 10 species was identified by RASs of the same (RAS 810) and other (RAS 814 and 815) species. HPV42 and 53 were detected by RASs of families alpha 1 (RAS 806) and alpha 6 (RAS 811), respectively. HPV18 and 45 were tested together, and each HPV was detected by its corresponding RAS (RAS809 for HPV45 and RAS810 for HPV18) without any ambiguity (data not shown). These findings indicate that one single RAS can detect and identify several HPVs; however, the association of RASs with the homonym sequences improved the genotyping of HPVs on their own or in mixtures.

### 4. Detection and genotyping of HPVs in clinical samples by RMAs

We assessed whether this RMA method could be used to identify HPVs in clinical samples. A set of 40 DNAs extracted from cervical smears with normal (n = 3) and abnormal (Atypical Squamous Cells of Unknown Significance (ASC-US), Atypical squamous cells cannot exclude HSIL (ASC-H), Low-grade Squamous Intraepithelial Lesion (LSIL) and High-grade Squamous Intraepithelial Lesion (HSIL)) cytology were provided by the French HPV Reference Laboratory (Table 6). For each sample, the HPV status was first determined with the PapilloCheck genotyping kit. Most of the samples contained HPV16 or HPV18 on their own or with other LR or HR HPVs. The copy number per cell was between 0.001 and 71 for HPV16 and 2 and 362 for HPV18, as determined by quantitative real-time PCR (Table 6). The DNA preparations were treated with PlasmidSafe DNase and amplified with phi29 polymerase, and then subjected to genotyping by RMAs with homonymous sequences and RASs. All assays with the Papillo-Check genotyping kit and RMAs were performed only one time except with sample n°2012-312. This DNA was independently amplified and hybridized on RMAs three times. Call rate and standard deviation (SD) values for HPV16 were determined with homonymous, non-homonymous and reconstructed ancestral sequences. SD values ranged from 0.9 to 1.9% according to the kind of considered sequences and were similar to values obtained with a triplicate assay using cloned HPV16 DNA (table 7).

The PapilloCheck genotyping kit can detect and identify 18 HR HPV (HPV 16, 18, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 70, 73, 82) and six LR HPV (HPV 6, 11, 40, 42, 43 and 54/

<footnote><sup>a</sup>The species corresponding to the RAS is given in parenthesis.</footnote>

<footnote><sup>b</sup>The call rate value for each RAS is indicated when a sequence was identified.
doi:10.1371/journal.pone.0109301.t005</footnote>
<table><caption>**Table 6.** Genotyping of HPV from cervical smears using the PapilloCheck kit and resequecing microarray.</caption><thead><tr><th rowspan="3">Samples</th><th rowspan="3">Cytology<sup>a</sup></th><th colspan="2">copy/cell<sup>b</sup></th><th colspan="3">HPV identified<sup>c</sup></th></tr><tr><th colspan="2">Papillocheck</th><th colspan="2">Resequencing Microarray (RMA)</th></tr><tr><th>HPV16</th><th>HPV18</th><th>Based on homonyme sequences</th><th>Based on reconstructed ancestral sequences</th></tr></thead><tbody><tr><td>2012-1557</td><td>N</td><td></td><td>Neg</td><td>Neg</td><td>Neg</td><td>Neg</td></tr><tr><td>2012-1558</td><td>N</td><td></td><td>Neg</td><td>Neg</td><td>Neg</td><td>Neg</td></tr><tr><td>2012-1559</td><td>N</td><td></td><td>Neg</td><td>Neg</td><td>Neg</td><td>Neg</td></tr><tr><td>2012-1930</td><td>ASCUS</td><td></td><td>Neg</td><td>Neg</td><td>Neg</td><td>Neg</td></tr><tr><td>2012-1931</td><td>ASCUS</td><td></td><td>Neg</td><td><strong>103</strong></td><td>Neg</td><td>Neg</td></tr><tr><td>2012-620</td><td>ASCUS</td><td>1.6</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-277</td><td>ASCUS</td><td>1.7</td><td>16, 42, 56</td><td>16, 42, 56</td><td>16, 42</td><td>16, 42</td></tr><tr><td>2012-1660</td><td>ASCUS</td><td></td><td>35, 45, 53, 56</td><td><strong>30, 35, 45, 53, 56, 67, 73</strong></td><td>35, 45, 53</td><td>35, 45, 53</td></tr><tr><td>2012-1649</td><td>ASCUS</td><td></td><td>6, 39, 42, 44/55, 45</td><td>6, 42, 45, 55, <strong>101</strong></td><td>6, 42, 45, 55</td><td>6, 42, 45, 55</td></tr><tr><td>2012-1620</td><td>ASCUS</td><td></td><td>40, 43, 44/55, 52, 56</td><td>40, 43, 52, 55, 56</td><td>43, 52, 55, 56</td><td>43, 52, 55, 56</td></tr><tr><td>2012-1642</td><td>ASCUS</td><td></td><td>40, 42, 51, 53, 59, 73</td><td><strong>35, 40, 42, 51, 53, 59, 67, 73</strong></td><td>40, 42, 51, 53, 59, <strong>67, 73</strong></td><td>40, 42, 51, 53, 59, <strong>67, 73</strong></td></tr><tr><td>2012-249</td><td>ASC-H</td><td>0.9</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-328</td><td>LSIL</td><td>0.5</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-146</td><td>LSIL</td><td></td><td>361.5</td><td>18</td><td><strong>18, 54, 74</strong></td><td><strong>18, 54, 74</strong></td></tr><tr><td>2012-133</td><td>LSIL</td><td></td><td>11.3</td><td>18</td><td>18</td><td>18</td></tr><tr><td>2012-352</td><td>LSIL</td><td>4.2</td><td>16, 35</td><td>16, 35, <strong>67</strong></td><td>16, 43, 67</td><td>16, 43, 67</td></tr><tr><td>2012-354</td><td>LSIL</td><td>84</td><td>16, 39</td><td>16, <strong>34, 39</strong></td><td>16, 39</td><td>16, 39</td></tr><tr><td>2012-601</td><td>LSIL</td><td></td><td>2.9</td><td>18, 33, 66</td><td>18, 33, 66</td><td>18, 33, 66</td></tr><tr><td>2012-351</td><td>LSIL</td><td>3.02</td><td></td><td>16, <strong>33, 52, 66</strong></td><td>16, <strong>34, 66</strong></td><td>16, <strong>34, 66</strong></td></tr><tr><td>2012-1940</td><td>LSIL</td><td></td><td></td><td>40, 42, 52, 66</td><td><strong>30, 40, 42, 52, 54, 66, 82, 91</strong></td><td>40, 52, <strong>82, 66</strong></td></tr><tr><td>2012-1607</td><td>LSIL</td><td></td><td></td><td>42, 51, 53, 66</td><td>42, 51, 53, 66</td><td>42, 51, 66</td></tr><tr><td>2012-445</td><td>LSIL</td><td>21</td><td></td><td>16, 39, 42, 59, 66</td><td>16, <strong>32, 39, 42, 59, 66</strong></td><td>16, 39, 42, 59, 66</td></tr><tr><td>2012-1673</td><td>LSIL</td><td></td><td></td><td><strong>33, 42, 51, 52, 66</strong></td><td>42, 51, 52, <strong>54, 66</strong></td><td>51, 52, 66</td></tr><tr><td>2012-451</td><td>HSIL</td><td>0.81</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-495</td><td>HSIL</td><td>0.47</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-622</td><td>HSIL</td><td>3.7</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-234</td><td>HSIL</td><td>3.6</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-312</td><td>HSIL</td><td>2.6</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-499</td><td>HSIL</td><td>5.0</td><td>16</td><td>16</td><td>16</td><td>16</td></tr><tr><td>2012-248</td><td>HSIL</td><td></td><td>2.1</td><td>18</td><td>18</td><td>18</td></tr><tr><td>2012-314</td><td>HSIL</td><td>0.8</td><td></td><td>16, <strong>42</strong></td><td>16</td><td>16</td></tr><tr><td>2012-243</td><td>HSIL</td><td>0.05</td><td></td><td>16, 56</td><td>16, 56</td><td>16, 56</td></tr><tr><td>2012-266</td><td>HSIL</td><td>0.001</td><td></td><td>16, 56</td><td>16, 56</td><td>16, 56</td></tr><tr><td>2012-610</td><td>HSIL</td><td>70.6</td><td></td><td>16, <strong>56</strong></td><td>16</td><td>16</td></tr><tr><td>2012-613</td><td>HSIL</td><td>4.02</td><td></td><td>16, 59</td><td>16, <strong>54, 59</strong></td><td>16, 59</td></tr><tr><td>2012-292</td><td>HSIL</td><td>3.06</td><td></td><td>16, <strong>40, 42</strong></td><td>16, 42</td><td>16</td></tr><tr><td>2012-1945</td><td>HSIL</td><td></td><td></td><td>39, 56, 66</td><td><strong>54, 56, 66</strong></td><td>39, 66</td></tr><tr><td>2012-453</td><td>HSIL</td><td>0.52</td><td></td><td>16, 33, <strong>42, 59</strong></td><td>16, 33, 59</td><td>16, 33</td></tr><tr><td>2012-504</td><td>HSIL</td><td></td><td>58.8</td><td>18, 39, 45, <strong>52, 56, 59, 70</strong></td><td>18, 39, 45, 56, 59</td><td>18, 52, 45, 56, 45</td></tr></tbody></table>

<footnote><sup>a</sup>ASCUS: Atypical Squamous Cells of Unknown Significance; ASCH: Atypical squamous cells cannot exclude HSIL; LSIL: Low-grade Squamous Intraepithelial Lesion; HSIL: High-grade Squamous Intraepithelial Lesion;</footnote>

<footnote><sup>b</sup>As determined by quantitative real-time PCR.</footnote>

<footnote><sup>c</sup>Neg: Negative; HPV genotypes detected by RMA and unidentifiable by the PapilloCheck kit are indicated in bold; HPV genotypes not detected by at least one of the two methods are underlined.</footnote>

<footnote>doi:10.1371/journal.pone.0109301.t006</footnote>
<table><caption>**Table 7.** Repeatability of RMA assay with cloned HPV16 DNA and biological sample with HPV16 monoinfection.</caption><thead><tr><th rowspan="3">HPV sequence tiled on RMA</th><th colspan="4">HPV16</th></tr><tr><th colspan="2">Cloned HPV DNA</th><th colspan="2">Sample n°312</th></tr><tr><th>Call rate<sup>a</sup></th><th>Standard deviation<sup>b</sup></th><th>Call rate<sup>a</sup></th><th>Standard deviation<sup>b</sup></th></tr></thead><tbody><tr><td>Homonymous</td><td>HPV16</td><td>89.6</td><td>2.2</td><td>91.4</td><td>1.6</td></tr><tr><td rowspan="2">Non-homonymous</td><td>HPV30</td><td>29.2</td><td>0.45</td><td>26.1</td><td>0.9</td></tr><tr><td>HPV33</td><td>28.0</td><td>1.4</td><td>18.7</td><td>1.6</td></tr><tr><td>Reconstructed Ancestral Sequence</td><td>RAS814</td><td>53</td><td>2.8</td><td>35.2</td><td>1.9</td></tr></tbody></table>

55). The two methods PapilloCheck genotyping and RMA genotyping with homonymous sequences gave concordant results for 28/40 (70%) of the samples, and compatible results (one or two HPV missed by one or other method) for 12/40 (30%) of the samples. A value of k = 0.36 was obtained with the kappa test in interrogation of agreement between these two tests. PapilloCheck did not detect some HPVs found by RMAs (HPV 35 in sample 2012-1642; HPV73 in sample 2012-1660; HPV82 in sample 2012-1940); the RMAs did not detect some strains identified by PapilloCheck (HPV39 in sample 2012-1649, HPV33 and HPV52 in sample 2012-351, HPV33 in sample 2012-1673, HPV42 in samples 2012-314 and 2012-453, HPV56 in sample 2012-610, HPV40 in sample 2012-292, and HPV52 and 70 in sample 2012-504). RMAs allowed two closely related HPVs, HPV44 and 55, to be distinguished and genotyped, although they are indistinguishable by PapilloCheck (samples n° 2012-1620 and 2012-1649). No discordant result was observed (Table 6).

The RMA method with homonymous sequences appeared to be highly sensitive: HPV16 was detected in a sample in which there was only 1 copy per 1000 cells (sample 2012-266). Eighteen HPVs not identifiable with the PapilloCheck genotyping kit were detected by RMAs in 14 samples (indicated in boldface in Table 6). In particular, two gamma HPVs (HPV 101 and 103) were identified in samples 2012-1931 and 2012-1649; these HPVs were first detected in cervico-vaginal smears [22].

HPV genotyping following RMAs by homonymous sequences and RASs gave compatible results, although fewer HPVs were identified with RASs, particularly when the number of HPVs detected by homonymous sequences was high. However, in one case, HPV39 first identified by PapilloCheck, was detected only by RASs (sample 2012-1945) (Table 6).

Discussion

Multiple broad-spectrum PCR methods have been developed for the detection of the alpha-HPV genus and most are based on HPV sequences in the L1 open reading frame. Here, we show that resequencing microarray (RMA) technology was feasible and effective for genotyping human papillomaviruses, and is thus a potentially valuable tool for epidemiology. We have demonstrated that the use of homonymous HPV sequences and reconstructed ancestral sequences (RASs), either separately or combined, allows precise molecular identification of different HPV genotypes. We validated the method first with cloned HPV DNA and subsequently for detection of HPVs in clinical samples. The RMA method used involved treatment of cellular DNA with Plasmid-Safe DNase and the universal Phi 29-based amplification; this methodology was highly sensitive for the genotyping of HPVs

in clinical samples, with a detection limit of one copy of HPV16 per 1000 cells.

This RMA method allows the detection and characterization of HPV genotypes other than those tiled on the microarray. Indeed, the call rate values were a function of the percentage of divergence between the tiled sequence and the sequence present in the sample. We found that there was no substantial loss of detection signal when the tiled sequence diverged by up to 10–15% from the sequence of the pathogen in the sample. Moreover, the RASs were designed to minimizing the divergence between the tiled and the viral sequences potentially present in clinical samples. Indeed, the use of RASs allowed an improvement in the detection of several HPVs. Nevertheless, the nucleotide diversity of HPV sequences is substantial and not a single RAS was able to detect all the different HPVs present in cases of multiple infections.

We compared this RMA method to a commercial HPV genotyping kit. It must be stressed that clinical samples with more than one HPV were intentionally selected to compare the two methods. Some HPV types identified in clinical samples by the RMA method were not identified by the PapilloCheck test. This discrepancy is in part explained by the number of HPV genotypes, limited to 24, which can be identified by the PapilloCheck test. The HPVs detected by RMAs but not by PapilloCheck tests include both LR HPVs (HPV 32, 34, 54, 74, and 91) and HR HPVs (HPV30, 34, and 67). Also, HR HPV35, 73, and 82, found in some clinical samples by RMAs, were not detected by the PapilloCheck test. Conversely, HPV33, 39, 52, 56, and 70 identifiable in some samples by the PapilloCheck kit were not found by the RMA method. This presumably reflects different sensitivities of the two tests for some HPV genotypes.

Our study has some limitations. First, due to the substantial genetic diversity of HPVs with more than 170 HPV genotypes already identified [1], not all alpha or gamma HPVs were studied, and HPVs belonging to beta, mu and nu genera were not studied. Nevertheless we were able to identify HPV101 and 103 that belong to the gamma 6 species, and that were initially detected in cervico-vaginal smears [20]. To increase the coverage of the known HPV diversity, and to facilitate detection of novel genotypes, the number of RASs should be increased. Second, in the current approach, only a part of the L1 gene for each HPV genotype was targeted and tiled on the microarray. The L1 gene may be deleted from the virus associated with some genital cancers, and this would lead to false negative results. The E6 and E7 genes have been described to be highly expressed in high-grade cervical lesions and cancer. Therefore, new generation RMA tests could include sequences corresponding to the E6 and E7 genes, and this would allow detection of both viral DNA and transcripts.

<footnote><sup>a</sup>Ratio between the number of bases determined and the sequence length.</footnote>

<footnote><sup>b</sup>Microarray assays were done in triplicate.</footnote>

<footnote>doi:10.1371/journal.pone.0109301.t007</footnote>
The cost of a high-density microarray could appear as a limit because this technology required the synthesis of a physical mask. However, the final cost can decrease if the demand increases, as may be the case in the context of a clinical diagnostic use. Moreover, the use of a RMA of the smallest size only dedicated to detect and characterize HPVs will help to reduce the global cost of assay. Adding of supplementary homonyms and reconstructed ancestral sequences will allow developing countries to obtain a better tool for the screening of HPVs. Based on the sequences obtained, RMAs may be useful for molecular epidemiological studies of HPV infection, particularly in geographical areas where the distribution of circulating genotypes has not yet been investigated. It can be used to identify variants of HPV genotypes that differ from the prototype until 10–15%. RMA-based tests may thus be informative about the prevalence and distribution of diverse HPV genotypes in the post-vaccination era.

## Supporting Information

<caption>**Table S1 Content description of VirID V3.0 RMA.**
(DOCX)</caption>

### References

1. de Villiers EM (2013) Cross-roads in the classification of papillomaviruses. Virology 445: 2–10.

2. Heard I, Tondeur L, Arowas L, Falguieres M, Demazoin MC, et al. (2013) Human papillomavirus types distribution in organised cervical cancer screening in France. PLoS One 8: e79372.

3. Bosch FX, Lorincz A, Munoz N, Meijer CJ, Shah KV (2002) The causal relation between human papillomavirus and cervical cancer. J Clin Pathol 55: 244–265.

4. Mighty KK, Laimins LA (2014) The role of human papillomaviruses in oncogenesis. Recent Results Cancer Res 193: 135–148.

5. Saslow D, Solomon D, Lawson HW, Killackey M, Kulasingam SL, et al. (2012) American Cancer Society, American Society for Colposcopy and Cervical Pathology, and American Society for Clinical Pathology screening guidelines for the prevention and early detection of cervical cancer. Am J Clin Pathol 137: 516–542.

6. Albrecht V, Chevallier A, Magnone V, Barbry P, Vandenbos F, et al. (2006) Easy and fast detection and genotyping of high-risk human papillomavirus by dedicated DNA microarrays. J Virol Methods 137: 236–244.

7. Lin ZH, Shen XH, Jin Z, Kim Y, Lee E, et al. (2005) Human papillomavirus genotyping by oligonucleotide microarray and p16 expression in uterine cervical intraepithelial neoplasm and in invasive carcinoma in Korean women. Pathol Int 55: 491–496.

8. Dacheux I, Berthet N, Dissard G, Holmes EC, Delmas O, et al. (2010) Application of broad-spectrum resequencing microarray for genotyping rhabdoviruses. J Virol 84: 9557–9574.

9. Wang Z, Daum LT, Vora GJ, Metzgar D, Walter EA, et al. (2006) Identifying influenza viruses with resequencing microarrays. Emerg Infect Dis 12: 638–646.

10. Wang Z, Malanoski AP, Lin B, Kidd C, Long NC, et al. (2008) Resequencing microarray probe design for typing genetically diverse viruses: human rhinoviruses and enteroviruses. BMC Genomics 9: 577.

11. Berthet N, Leclercq I, Dublineau A, Shigematsu S, Burguiere AM, et al. (2010) High-density resequencing DNA microarrays in public health emergencies. Nat Biotechnol 28: 25–27.

<caption>**Table S2 RASs tiled on the resequencing microarray.**
(DOCX)</caption>

### Acknowledgments

This program was supported by the *Programme Transversal de Recherche* DEVA No. 246 financed by the Institut Pasteur (Paris, France), sponsored by the Fondation Total-Institut Pasteur and the Elisabeth Taub award from the *Académie Nationale de Médécine*. We thank the group “Genotyping of Pathogens” for the use of the Affymetrix station and F. Fichenick for microarray experiments. This study was funded by the French Institut National du Cancer (Canceropole 2011-1VADS 03-IP1). The funders had no role in study design, data analysis or preparation of the manuscript.

### Author Contributions

Conceived and designed the experiments: NB AG M. Favre. Performed the experiments: NB M. Falguieres SB. Analyzed the data: NB CF CB CBF AG IH M. Favre. Contributed reagents/materials/analysis tools: SB. Wrote the paper: NB AG M. Favre.

12. Leclercq I, Berthet N, Batejat C, Rousseaux C, Dickinson P, et al. (2010) Use of consensus sequences for the design of high density resequencing microarrays: the influenza virus paradigm. BMC Genomics 11: 586.

13. Berthet N, Dickinson P, Filliol I, Reinhardt AK, Batejat C, et al. (2008) Massively parallel pathogen identification using high-density microarrays. Microb Biotechnol 1: 79–86.

14. Berthet N, Nakoune E, Whist E, Selekon B, Burguiere AM, et al. (2011) Maculopapular lesions in the Central African Republic. Lancet 378: 1354.

15. Berthet N, Paulous S, Coffey LL, Frenkiel MP, Moltini I, et al. (2012) Resequencing microarray method for molecular diagnosis of human arboviral diseases. J Clin Virol.

16. Filippone C, Marianneau P, Murri S, Mollard N, Zupanc TA, et al. (2012) Molecular diagnostic and genetic characterization of highly pathogenic viruses: application during Crimean Congo hemorrhagic fever virus outbreaks in Eastern Europe and Middle-East. Clinical Microbiology and Infection.

17. Berthet N, Deletoile A, Passet V, Kennedy GC, Manuguerra JC, et al. (2010) Reconstructed ancestral sequences improve pathogen identification using resequencing DNA microarrays. PLoS One 5: e15243.

18. Yang Z (2007) PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol 24: 1586–1591.

19. Mendoza JA, Jacob Y, Cassonnet P, Favre M (2006) Human papillomavirus type 5 E6 oncoprotein represses the transforming growth factor beta signaling pathway by binding to SMAD3. J Virol 80: 12420–12424.

20. Cutler DJ, Zwick ME, Carrasquillo MM, Yohn CT, Tobin KP, et al. (2001) High-throughput variation detection and genotyping using microarrays. Genome Res 11: 1913–1925.

21. Schmitz M, Scheungraber C, Herrmann J, Teller K, Gajda M, et al. (2009) Quantitative multiplex PCR assay for the detection of the seven clinically most relevant high-risk HPV types. J Clin Virol 44: 302–307.

22. Chen Z, Schiffman M, Herrero R, Desalle R, Burk RD (2007) Human papillomavirus (HPV) types 101 and 103 isolated from cervicovaginal cells lack an E6 open reading frame (ORF) and are related to gamma-papillomaviruses. Virology 360: 447–453.