An optimized methodology for whole genome sequencing of RNA respiratory viruses from nasopharyngeal aspirates

VALINOTTO, LAURA E.; NABAES JODAR, MERCEDES S.; MARTI, MARCELO A.; GOYA, STEPHANIE; ROJO, GABRIEL L.; ZAIAT, JONATHAN J.; VIEGAS, MARIANA; TITTARELLI, ESTEFANIA; GRENINGER, ALEXANDER L.; MISTCHENKO, ALICIA S.; VALINOTTO, LAURA E.; NABAES JODAR, MERCEDES S.; MARTI, MARCELO A.; GOYA, STEPHANIE; ROJO, GABRIEL L.; ZAIAT, JONATHAN J.; VIEGAS, MARIANA; TITTARELLI, ESTEFANIA; GRENINGER, ALEXANDER L.; MISTCHENKO, ALICIA S.

PLOS ONE

PUBLIC LIBRARY SCIENCE

Año: 2018 vol. 13

1932-6203

Over the last decade, the number of viral genome sequences deposited in available databases has grown exponentially. However, sequencing methodology vary widely and many published works have relied on viral enrichment by viral culture or nucleic acid amplification with specific primers rather than through unbiased techniques such as metagenomics. The genome of RNA viruses is highly variable and these enrichment methodologies may be difficult to achieve or may bias the results. In order to obtain genomic sequences of human respiratory syncytial virus (HRSV) from positive nasopharyngeal aspirates diverse methodologies were evaluated and compared. A total of 29 nearly complete and complete viral genomes were obtained. The best performance was achieved with a DNase I treatment to the RNA directly extracted from the nasopharyngeal aspirate (NPA), sequence-independent single-primer amplification (SISPA) and library preparation performed with Nextera XT DNA Library Prep Kit with manual normalization. An average of 633,789 and 1,674,845 filtered reads per library were obtained with MiSeq and NextSeq 500 platforms, respectively. The higher output of NextSeq 500 was accompanied by the increasing of duplicated reads percentage generated during SISPA (from an average of 1.5% duplicated viral reads in MiSeq to an average of 74% in NextSeq 500). HRSV genome recovery was not affected by the presence or absence of duplicated reads but the computational demand during the analysis was increased. Considering that only samples with viral load E+06 copies/ml NPA were tested, no correlation between sample viral loads and number of total filtered reads was observed, nor with the mapped viral reads. The HRSV genomes showed a mean coverage of 98.46% with the best methodology. In addition, genomes of human metapneumovirus (HMPV), human rhinovirus (HRV) and human parainfluenza virus types 1?3 (HPIV1-3) were also obtained with the selected optimal methodology.