LAUSR

Genomic experiments analyzing human papillomaviruses (HPVs) require a carefully selected list of sequences as a reference database to map millions of reads. The available sources, such as the Papillomavirus Episteme (PaVE), are organized based on variations in the L1 gene rather than the whole HPV sequence. Moreover, the PaVE process uses complex multiple sequence alignments containing hundreds or thousands of sequences. These issues complicate the generation of a reference database for genomics, leading to the generation of per-analysis-defined databases. Here, we propose a de novo strategy considering all HPV sequences reported in the NCBI database to define a subset of highly representative HPV sequences. The strategy is based on oligonucleotide frequency profiling of the whole sequence followed by hierarchical clustering. Using data from HPV capture experiments, we demonstrate that this strategy selects suitable sequences as a reference database to map most mappable reads unambiguously. We provide some recommendations to improve HPV mapping. The generated .fasta files can be accessed at https://github.com/vtrevino/HPV-Ref-Genomes .