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

ViroSpot microneutralization assay for antigenic characterization of human influenza viruses

The hemagglutination inhibition (HI) assay has been used for the antigenic characterization of influenza viruses for decades. However, the majority of recent seasonal influenza A viruses of the H3N2 subtype has lost the capacity to agglutinate erythrocytes of various species. The hemagglutination (HA) activity of other A(H3N2) strains is generally sensitive to the action of the neuraminidase inhibitor oseltamivir, which indicates that the neuraminidase and not the hemagglutinin is responsible for the HA activity. These findings complicate the antigenic characterization and selection of A(H3N2) vaccine strains, calling for alternative antigenic characterization assays. Here we describe the development and use of the ViroSpot microneutralization (MN) assay as a reliable and robust alternative for the HI assay. Serum neutralization of influenza A(H3N2) reference virus strains and epidemic isolates was determined by automated readout of immunostained cell monolayers, in a format designed to minimize the influence of infectious virus doses on serum neutralization titers. Neutralization of infection was largely independent from rates of viral replication and cell-to-cell transmission, facilitating the comparison of different virus isolates. Other advantages of the ViroSpot MN assay include its relative insensitivity to variation in test dose of infectious virus, automated capture and analyses of residual infection patterns, and compatibility with standardized large scale analyses. Using this assay, a number of epidemic influenza A(H3N2) strains that failed to agglutinate erythrocytes, were readily characterized antigenically

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