Improved detection of artifactual viral minority variants in high-throughput sequencing data

High-throughput sequencing (HTS) of viral samples provides important information on the presence of viral minority variants. However, detection and accurate quantification is limited by the capacity to distinguish biological from artificial variation. In this study, errors related to the Illumina Hiseq2000 library generation and HTS process were investigated by determining minority variant frequencies in an influenza A/WSN/1933(H1N1) virus reversegenetics plasmid pool. Errors related to amplification and sequencing were determined using the same plasmid pool, by generation of infectious virus using reverse genetics followed by in duplo reverse-transcriptase PCR (RT-PCR) amplification and HTS in the same sequence run. Results showed that after 'best practice' quality control (QC), within the plasmid pool, 1 minority variant with a frequency >0.5% was identified, while 84 and 139 were identified in the RT-PCR amplified samples, indicating RT-PCR amplification artificially increased variation. Detailed analysis showed that artifactual minority variants could be identified by two major technical characteristics: their predominant presence in a single read orientation and uneven distribution of mismatches over the length of the reads. We demonstrate that by addi

Improved detection of artifactual viral minority variants in high-throughput sequencing data

High-throughput sequencing (HTS) of viral samples provides important information on the presence of viral minority variants. However, detection and accurate quantification is limited by the capacity to distinguish biological from artificial variation. In this study, errors related to the Illumina HiSeq2000 library generation and HTS process were investigated by determining minority variant frequencies in an influenza A/WSN/1933(H1N1) virus reverse-genetics plasmid pool. Errors related to amplification and sequencing were determined using the same plasmid pool, by generation of infectious virus using reverse genetics followed by in duplo reverse-transcriptase PCR (RT-PCR) amplification and HTS in the same sequence run. Results showed that after "best practice" quality control (QC), within the plasmid pool, one minority variant with a frequency >0.5% was identified, while 84 and 139 were identified in the RT-PCR amplified samples, indicating RT-PCR amplification artificially increased variation. Detailed analysis showed that artifactual minority variants could be identified by two major technical characteristics: their predominant presence in a single read orientation and uneven distribution of mismatches over the length of the reads. We demonstrate that by addition of two QC steps 95% of the artifactual minority variants could be identified. When our analysis approach was applied to three clinical samples 68% of the initially identified minority variants were identified as artifacts. Our study clearly demonstrated that, without additional QC steps, overestimation of viral minority variants is very likely to occur, m

Emergence of the virulence-associated PB2 E627K substitution in a fatal human case of highly pathogenic avian influenza virus A(H7N7) infection as determined by illumina ultra-deep sequencing

Avian influenza viruses are capable of crossing the species barrier and infecting humans. Although evidence of human-to-human transmission of avian influenza viruses to date is limited, evolution of variants toward more-efficient human-to-human transmission could result in a new influenza virus pandemic. In both the avian influenza A(H5N1) and the recently emerging avian influenza A(H7N9) viruses, the polymerase basic 2 protein (PB2) E627K mutation appears to be of key importance for human adaptation. During a large influenza A(H7N7) virus outbreak in the Netherlands in 2003, the A(H7N7) virus isolated from a fatal human case contained the PB2 E627K mutation as well as a hemagglutinin (HA) K416R mutation. In this study, we aimed to investigate whether these mutations occurred in the avian or the human host by Illumina Ultra-Deep sequencing of three previously uninvestigated clinical samples obtained from the fatal case. In addition, we investigated three chicken samples, two of which were obtained from the source farm. Results showed that the PB2 E627K mutation was not present in any of the chicken samples tested. Surprisingly, the avian samples were characterized by the presence of influenza virus defective RNA segments, suggestive for the synthesis of defective interfering viruses during infection in poultry. In the human samples, the PB2 E627K mutation was identified with increasing frequency during infection. Our results strongly suggest that human adaptation marker PB2 E627K has emerged during virus infection of a single human host, emphasizing the importance of reducing human exposure to avian influenza viruses to reduce the likelihood of viral adaptation to humans