Viral population analysis and minority-variant detection using short read next-generation sequencing

RNA viruses within infected individuals exist as a population of evolutionary-related variants. Owing to evolutionary change affecting the constitution of this population, the frequency and/or occurrence of individual viral variants can show marked or subtle fluctuations. Since the development of massively parallel sequencing platforms, such viral populations can now be investigated to unprecedented resolution. A critical problem with such analyses is the presence of sequencing-related errors that obscure the identification of true biological variants present at low frequency. Here, we report the development and assessment of the Quality Assessment of Short Read (QUASR) Pipeline (http://sourceforge.net/projects/quasr) specific for virus genome short read analysis that minimizes sequencing errors from multiple deep-sequencing platforms, and enables post-mapping analysis of the minority variants within the viral population. QUASR significantly reduces the error-related noise in deep-sequencing datasets, resulting in increased mapping accuracy and reduction of erroneous mutations. Using QUASR, we have determined influenza virus genome dynamics in sequential samples from an in vitro evolution of 2009 pandemic H1N1 (A/H1N1/09) influenza from samples sequenced on both the Roche 454 GSFLX and Illumina GAIIx platforms. Importantly, concordance between the 454 and Illumina sequencing allowed unambiguous minority-variant detection and accurate determination of virus population turnover in vitro

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

Pandemic H1N1 virus transmission and shedding dynamics in index case households of a prospective Vietnamese cohort

Objectives: Influenza household transmission studies are required to guide prevention strategies but most passively recruit index cases that seek healthcare. We investigated A(H1N1)pdm09 transmission in a household-based cohort during 2009. Methods: Health-workers visited 270 households weekly, and collected swabs from influenza-like-illness cases. If A(H1N1)pdm09 was RT-PCR-confirmed, all household members had symptoms assessed and swabs collected daily for 10-15 days. Viral RNA was quantified and sequenced and serology performed on pre-pandemic sera. Results: Index cases were detected in 20 households containing 81 people. 98.5% lacked A(H1N1)pdm09 neutralizing antibodies in pre-pandemic sera. Eleven (18.6%, 95% CI 10.7-30.4%) of 59 contacts were infected. Virus genetic diversity within households was negligible and less than between households. Index and secondary cases were distributed between mothers, daughters and sons, and had similar virus-RNA shedding and symptom dynamics. Fathers were rarely infected. Five secondary cases (45%) had no apparent symptoms and three shed virus before symptoms. Secondary infection was associated with index case wet cough (OR 1.56, 95% CI 1.22-1.99). Conclusions: In this cohort of A(H1N1)pdm09 susceptible persons, virus sequencing was capable of discriminating household from community transmission. Household transmission involved mothers and children but rarely fathers. Asymptomatic or pre-symptomatic shedding was common. © 2014 The Authors

Pandemic H1N1 virus transmission and shedding dynamics in index case households of a prospective Vietnamese cohort

OBJECTIVES: Influenza household transmission studies are required to guide prevention strategies but most passively recruit index cases that seek healthcare. We investigated A(H1N1)pdm09 transmission in a household-based cohort during 2009. METHODS: Health-workers visited 270 households weekly, and collected swabs from influenza-like-illness cases. If A(H1N1)pdm09 was RT-PCR-confirmed, all household members had symptoms assessed and swabs collected daily for 10-15 days. Viral RNA was quantified and sequenced and serology performed on pre-pandemic sera. RESULTS: Index cases were detected in 20 households containing 81 people. 98.5% lacked A(H1N1)pdm09 neutralizing antibodies in pre-pandemic sera. Eleven (18.6%, 95% CI 10.7-30.4%) of 59 contacts were infected. Virus genetic diversity within households was negligible and less than between households. Index and secondary cases were distributed between mothers, daughters and sons, and had similar virus-RNA shedding and symptom dynamics. Fathers were rarely infected. Five secondary cases (45%) had no apparent symptoms and three shed virus before symptoms. Secondary infection was associated with index case wet cough (OR 1.56, 95% CI 1.22-1.99). CONCLUSIONS: In this cohort of A(H1N1)pdm09 susceptible persons, virus sequencing was capable of discriminating household from community transmission. Household transmission involved mothers and children but rarely fathers. Asymptomatic or pre-symptomatic shedding was common