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

Autor(es)

Castelli, Erick C. ^{[1, 2]} ; Paz, Michelle A. ^[1] ; Souza, Andreia S. ^[1] ; Ramalho, Jaqueline ^[1] ; Mendes-Junior, Celso Teixeira ^[3]

Afiliação do(s) autor(es)

^[1] Sao Paulo State Univ UNESP, Mol Genet & Bioinformat Lab, Expt Res Unit UNIPEX, Sch Med, Botucatu, SP - Brazil

^[2] Sao Paulo State Univ UNESP, Pathol Dept, Sch Med, BR-18618687 Botucatu, SP - Brazil

^[3] Univ Sao Paulo, Dept Quim, Fac Filosofia Ciencias & Letras Ribeirao Preto, BR-14040901 Ribeirao Preto, SP - Brazil

• Tipo de documento: Artigo Científico
• Fonte: HUMAN IMMUNOLOGY; v. 79, n. 9, p. 678-684, SEP 2018.
• Citações Web of Science: 4

Resumo

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. Ma-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. (AU)

• Processo FAPESP: 13/17084-2 - Estudo da variabilidade, controle da expressão e história evolutiva dos genes de Classe I do MHC humano
• Beneficiário: Erick da Cruz Castelli
• Modalidade de apoio: Auxílio à Pesquisa - Regular