Hla-mapper: An application to optimize the mapping of HLA sequences produced by massively parallel sequencing procedures

A challenging task when more than one HLA gene is evaluated together by second-generation sequencing is to achieve a reliable read mapping. The polymorphic and repetitive nature of HLA genes might bias the read mapping process, usually underestimating variability at very polymorphic segments, or ove...

Descripción completa

Detalles Bibliográficos
Autores: Castelli, Erick C. [UNESP], Paz, Michelle A. [UNESP], Souza, Andréia S. [UNESP], Ramalho, Jaqueline [UNESP], Mendes-Junior, Celso Teixeira
Tipo de recurso: artículo
Estado:Versión publicada
Fecha de publicación:2018
País:Brasil
Institución:Universidade Estadual Paulista (UNESP)
Repositorio:Repositório Institucional da UNESP
Idioma:inglés
OAI Identifier:oai:repositorio.unesp.br:11449/171178
Acceso en línea:http://dx.doi.org/10.1016/j.humimm.2018.06.010
http://hdl.handle.net/11449/171178
Access Level:acceso abierto
Palabra clave:Aligners
HLA
Mapping tool
MHC
Next Generation Sequencing (NGS)
Polymorphisms
Second Generation Sequencing
Typing
Variability
Descripción
Sumario:A challenging task when more than one HLA gene is evaluated together by second-generation sequencing is to achieve a reliable read mapping. The polymorphic and repetitive nature of HLA genes might bias the read mapping process, usually underestimating variability at very polymorphic segments, or overestimating variability at some segments. To overcome this issue we developed hla-mapper, which takes into account HLA sequences derived from the IPD-IMGT/HLA database and unpublished HLA sequences to apply a scoring system. This comprehends the evaluation of each read pair, addressing them to the most likely HLA gene they were derived from. Hla-mapper provides a reliable map of HLA sequences, allowing accurate downstream analysis such as variant calling, haplotype inference, and allele typing. Moreover, hla-mapper supports whole genome, exome, and targeted sequencing data. To assess the software performance in comparison with traditional mapping algorithms, we used three different simulated datasets to compare the results obtained with hla-mapper, BWA MEM, and Bowtie2. Overall, hla-mapper presented a superior performance, mainly for the classical HLA class I genes, minimizing wrong mapping and cross-mapping that are typically observed when using BWA MEM or Bowtie2 with a single reference genome.