Multiple Introductions and Antigenic Mismatch with Vaccines May Contribute to Increased Predominance of G12P[8] Rotaviruses in the United States.

Rotavirus is the leading global cause of diarrheal mortality for unvaccinated children under 5 years of age. The outer capsid of rotavirus virions consists of VP7 and VP4 proteins, which determine viral G and P types, respectively, and are primary targets of neutralizing antibodies. Successful vaccination depends upon generating broadly protective immune responses following exposure to rotaviruses presenting a limited number of G- and P-type antigens. Vaccine introduction resulted in decreased rotavirus disease burden but also coincided with the emergence of uncommon G and P genotypes, including G12. To gain insight into the recent predominance of G12P[8] rotaviruses in the United States, we evaluated 142 complete rotavirus genome sequences and metadata from 151 clinical specimens collected in Nashville, TN, from 2011 to 2013 through the New Vaccine Surveillance Network. Circulating G12P[8] strains were found to share many segments with other locally circulating strains but to have distinct constellations. Phylogenetic analyses of G12 sequences and their geographic sources provided evidence for multiple separate introductions of G12 segments into Nashville, TN. Antigenic epitopes of VP7 proteins of G12P[8] strains circulating in Nashville, TN, differ markedly from those of vaccine strains. Fully vaccinated children were found to be infected with G12P[8] strains more frequently than with other rotavirus genotypes. Multiple introductions and significant antigenic mismatch may in part explain the recent predominance of G12P[8] strains in the United States and emphasize the need for continued monitoring of rotavirus vaccine efficacy against emerging rotavirus genotypes.IMPORTANCE Rotavirus is an important cause of childhood diarrheal disease worldwide. Two immunodominant proteins of rotavirus, VP7 and VP4, determine G and P genotypes, respectively. Recently, G12P[8] rotaviruses have become increasingly predominant. By analyzing rotavirus genome sequences from stool specimens obtained in Nashville, TN, from 2011 to 2013 and globally circulating rotaviruses, we found evidence of multiple introductions of G12 genes into the area. Based on sequence polymorphisms, VP7 proteins of these viruses are predicted to present themselves to the immune system very differently than those of vaccine strains. Many of the sick children with G12P[8] rotavirus in their diarrheal stools also were fully vaccinated. Our findings emphasize the need for continued monitoring of circulating rotaviruses and the effectiveness of the vaccines against strains with emerging G and P genotypes.

Avian Influenza Viruses in Wild Birds: Virus Evolution in a Multihost Ecosystem.

Wild ducks and gulls are the major reservoirs for avian influenza A viruses (AIVs). The mechanisms that drive AIV evolution are complex at sites where various duck and gull species from multiple flyways breed, winter, or stage. The Republic of Georgia is located at the intersection of three migratory flyways: the Central Asian flyway, the East Africa/West Asia flyway, and the Black Sea/Mediterranean flyway. For six complete study years (2010 to 2016), we collected AIV samples from various duck and gull species that breed, migrate, and overwinter in Georgia. We found a substantial subtype diversity of viruses that varied in prevalence from year to year. Low-pathogenic AIV (LPAIV) subtypes included H1N1, H2N3, H2N5, H2N7, H3N8, H4N2, H6N2, H7N3, H7N7, H9N1, H9N3, H10N4, H10N7, H11N1, H13N2, H13N6, H13N8, and H16N3, and two highly pathogenic AIVs (HPAIVs) belonging to clade 2.3.4.4, H5N5 and H5N8, were found. Whole-genome phylogenetic trees showed significant host species lineage restriction for nearly all gene segments and significant differences in observed reassortment rates, as defined by quantification of phylogenetic incongruence, and in nucleotide sequence diversity for LPAIVs among different host species. Hemagglutinin clade 2.3.4.4 H5N8 viruses, which circulated in Eurasia during 2014 and 2015, did not reassort, but analysis after their subsequent dissemination during 2016 and 2017 revealed reassortment in all gene segments except NP and NS. Some virus lineages appeared to be unrelated to AIVs in wild bird populations in other regions, with maintenance of local AIVs in Georgia, whereas other lineages showed considerable genetic interrelationships with viruses circulating in other parts of Eurasia and Africa, despite relative undersampling in the area.IMPORTANCE Waterbirds (e.g., gulls and ducks) are natural reservoirs of avian influenza viruses (AIVs) and have been shown to mediate the dispersal of AIVs at intercontinental scales during seasonal migration. The segmented genome of influenza viruses enables viral RNA from different lineages to mix or reassort when two viruses infect the same host. Such reassortant viruses have been identified in most major human influenza pandemics and several poultry outbreaks. Despite their importance, we have only recently begun to understand AIV evolution and reassortment in their natural host reservoirs. This comprehensive study illustrates AIV evolutionary dynamics within a multihost ecosystem at a stopover site where three major migratory flyways intersect. Our analysis of this ecosystem over a 6-year period provides a snapshot of how these viruses are linked to global AIV populations. Understanding the evolution of AIVs in the natural host is imperative to mitigating both the risk of incursion into domestic poultry and the potential risk to mammalian hosts, including humans.

Zika Virus Antagonizes Type I Interferon Responses during Infection of Human Dendritic Cells.

Zika virus (ZIKV) is an emerging mosquito-borne flavivirus that is causally linked to severe neonatal birth defects, including microcephaly, and is associated with Guillain-Barre syndrome in adults. Dendritic cells (DCs) are an important cell type during infection by multiple mosquito-borne flaviviruses, including dengue virus, West Nile virus, Japanese encephalitis virus, and yellow fever virus. Despite this, the interplay between ZIKV and DCs remains poorly defined. Here, we found human DCs supported productive infection by a contemporary Puerto Rican isolate with considerable variability in viral replication, but not viral binding, between DCs from different donors. Historic isolates from Africa and Asia also infected DCs with distinct viral replication kinetics between strains. African lineage viruses displayed more rapid replication kinetics and infection magnitude as compared to Asian lineage viruses, and uniquely induced cell death. Infection of DCs with both contemporary and historic ZIKV isolates led to minimal up-regulation of T cell co-stimulatory and MHC molecules, along with limited secretion of inflammatory cytokines. Inhibition of type I interferon (IFN) protein translation was observed during ZIKV infection, despite strong induction at the RNA transcript level and up-regulation of other host antiviral proteins. Treatment of human DCs with RIG-I agonist potently restricted ZIKV replication, while type I IFN had only modest effects. Mechanistically, we found all strains of ZIKV antagonized type I IFN-mediated phosphorylation of STAT1 and STAT2. Combined, our findings show that ZIKV subverts DC immunogenicity during infection, in part through evasion of type I IFN responses, but that the RLR signaling pathway is still capable of inducing an antiviral state, and therefore may serve as an antiviral therapeutic target.

Reassortment of Influenza A Viruses in Wild Birds in Alaska before H5 Clade 2.3.4.4 Outbreaks.

Sampling of mallards in Alaska during September 2014-April 2015 identified low pathogenic avian influenza A virus (subtypes H5N2 and H1N1) that shared ancestry with highly pathogenic reassortant H5N2 and H5N1 viruses. Molecular dating indicated reassortment soon after interhemispheric movement of H5N8 clade 2.3.4.4, suggesting genetic exchange in Alaska or surrounds before outbreaks.

Influenza A virus evolution and spatio-temporal dynamics in Eurasian wild birds: a phylogenetic and phylogeographical study of whole-genome sequence data.

Low pathogenic avian influenza A viruses (IAVs) have a natural host reservoir in wild waterbirds and the potential to spread to other host species. Here, we investigated the evolutionary, spatial and temporal dynamics of avian IAVs in Eurasian wild birds. We used whole-genome sequences collected as part of an intensive long-term Eurasian wild bird surveillance study, and combined this genetic data with temporal and spatial information to explore the virus evolutionary dynamics. Frequent reassortment and co-circulating lineages were observed for all eight genomic RNA segments over time. There was no apparent species-specific effect on the diversity of the avian IAVs. There was a spatial and temporal relationship between the Eurasian sequences and significant viral migration of avian IAVs from West Eurasia towards Central Eurasia. The observed viral migration patterns differed between segments. Furthermore, we discuss the challenges faced when analysing these surveillance and sequence data, and the caveats to be borne in mind when drawing conclusions from the apparent results of such analyses.

Molecular characterization of respiratory syncytial viruses circulating in a paediatric cohort in Amman, Jordan.

Respiratory syncytial viruses (RSVs) are an important cause of mortality worldwide and a major cause of respiratory tract infections in children, driving development of vaccine candidates. However, there are large gaps in our knowledge of the local evolutionary and transmission dynamics of RSVs, particularly in understudied regions such as the Middle East. To address this gap, we sequenced the complete genomes of 58 RSVA and 27 RSVB samples collected in a paediatric cohort in Amman, Jordan, between 2010 and 2013. RSVA and RSVB co-circulated during each winter epidemic of RSV in Amman, and each epidemic comprised multiple independent viral introductions of RSVA and RSVB. However, RSVA and RSVB alternated in dominance across years, potential evidence of immunological interactions. Children infected with RSVA tended to be older than RSVB-infected children [30 months versus 22.4 months, respectively (P value = 0.02)], and tended to developed bronchopneumonia less frequently than those with RSVB, although the difference was not statistically significant (P value = 0.06). Differences in spatial patterns were investigated, and RSVA lineages were often identified in multiple regions in Amman, whereas RSVB introductions did not spread beyond a single region of the city, although these findings were based on small sample sizes. Multiple RSVA genotypes were identified in Amman, including GA2 viruses as well as three viruses from the ON1 sub-genotype that emerged in 2009 and are now the dominant genotype circulating worldwide. As vaccine development advances, further sequencing of RSV is needed to understand viral ecology and transmission, particularly in under-studied locations.

Evolutionary dynamics and molecular epidemiology of West Nile virus in New York State: 1999-2015.

Following its introduction into New York State (NYS) in 1999, West Nile virus (WNV; Flavivirus, Flaviviridae) underwent a rapid expansion throughout the USA and into Canada and Latin America. WNV has been characterized as being evolutionarily stable, with weak geographic structure, a dominance of purifying selection and limited adaptive change. We analyzed all available full-genome WNV sequences, focusing on the 543 available sequences from NYS, which included 495 newly sequenced 2000-15 isolates. In addition, we analyzed deep-sequencing data from 317 of these isolates. While our data are generally in agreement with the limited pace of evolutionary change and broad geographic and temporal mixing identified in other studies, we have identified some important exceptions. Most notably, there are 14 codons which demonstrated evidence of positive selection as determined by multiple models, including some positions with evidence of selection in NYS exclusively. Coincident with increased WNV activity, genotypes possessing one or more of these mutations, designated NY01, NY07, and NY10, have increased in prevalence in recent years and displaced historic strains. In addition, we have found a geographical bias with many of these mutations, which suggests selective pressures and adaptations could be regional. Lastly, our deep-sequencing data suggest both increased overall diversity in avian tissue isolates relative to mosquito isolates and multiple non-synonymous minority variants that are both host-specific and retained over time and space. Together, these data provide novel insight into the evolutionary pressures on WNV and the need for continued genetic surveillance and characterization of emergent strains.

Erratum: Evolutionary dynamics and molecular epidemiology of West Nile virus in New York State: 1999-2015.

[This corrects the article DOI: 10.1093/ve/vez020.].

Growth and adaptation of Zika virus in mammalian and mosquito cells.

The recent emergence of Zika virus (ZIKV) in the Americas coincident with increased caseloads of microcephalic infants and Guillain-Barre syndrome has prompted a flurry of research on ZIKV. Much of the research is difficult to compare or repeat because individual laboratories use different virus isolates, growth conditions, and quantitative assays. Here we obtained three readily available contemporary ZIKV isolates and the prototype Ugandan isolate. We generated stocks of each on Vero mammalian cells (ZIKVmam) and C6/36 mosquito cells (ZIKVmos), determined titers by different assays side-by-side, compared growth characteristics using one-step and multi-step growth curves on Vero and C6/36 cells, and examined plaque phenotype. ZIKV titers consistently peaked earlier on Vero cells than on C6/36 cells. Contemporary ZIKV isolates reached peak titer most quickly in a multi-step growth curve when the amplifying cell line was the same as the titering cell line (e.g., ZIKVmam titered on Vero cells). Growth of ZIKVmam on mosquito cells was particularly delayed. These data suggest that the ability to infect and/or replicate in insect cells is limited after growth in mammalian cells. In addition, ZIKVmos typically had smaller, more homogenous plaques than ZIKVmam in a standard plaque assay. We hypothesized that the plaque size difference represented early adaptation to growth in mammalian cells. We plaque purified representative-sized plaques from ZIKVmos and ZIKVmam. ZIKVmos isolates maintained the initial phenotype while plaques from ZIKVmam isolates became larger with passaging. Our results underscore the importance of the cells used to produce viral stocks and the potential for adaptation with minimal cell passages. In addition, these studies provide a foundation to compare current and emerging ZIKV isolates in vitro and in vivo.

Whole genome sequencing, variant analysis, phylogenetics, and deep sequencing of Zika virus strains.

The recent emergence of Zika virus (ZIKV) has been concentrated in the Caribbean, Southeastern United States, and South- and Central America; resulting in travel-based cases being reported around the globe. As multi-disciplinary collaborations are combatting the ZIKV outbreak, the need to validate the sequence of existing strains has become apparent. Here, we report high-quality sequence data for multiple ZIKV strains made publicly available through the National Institutes of Health- (NIH) funded biorepository, BEI Resources (www.beiresources.org). Next-generation sequencing, 3' rapid amplification of cDNA ends (RACE), and viral genome annotation pipelines generated GenBank sequence records for 16 BEI Resources strains. Minor variants, consensus mutations, and consensus insertions/deletions were identified within the viral stocks using next-generation sequencing (NGS) and consensus changes were confirmed with Sanger sequencing. Bioinformatics analyses of the sequencing results confirm that the virus stocks available to the scientific research community through BEI Resources adequately represent the viral population diversity of ZIKV.

Whole-genome sequencing and analyses identify high genetic heterogeneity, diversity and endemicity of rotavirus genotype P[6] strains circulating in Africa.

Rotavirus A (RVA) exhibits a wide genotype diversity globally. Little is known about the genetic composition of genotype P[6] from Africa. This study investigated possible evolutionary mechanisms leading to genetic diversity of genotype P[6] VP4 sequences. Phylogenetic analyses on 167 P[6] VP4 full-length sequences were conducted, which included six porcine-origin sequences. Of the 167 sequences, 57 were newly acquired through whole genome sequencing as part of this study. The other 110 sequences were all publicly-available global P[6] VP4 full-length sequences downloaded from GenBank. The strength of association between the phenotypic features and the phylogeny was also determined. A number of reassortment and mixed infections of RVA genotype P[6] strains were observed in this study. Phylogenetic analyses demostrated the extensive genetic diversity that exists among human P[6] strains, porcine-like strains, their concomitant clades/subclades and estimated that P[6] VP4 gene has a higher substitution rate with the mean of 1.05E-3 substitutions/site/year. Further, the phylogenetic analyses indicated that genotype P[6] strains were endemic in Africa, characterised by an extensive genetic diversity and long-time local evolution of the viruses. This was also supported by phylogeographic clustering and G-genotype clustering of the P[6] strains when Bayesian Tip-association Significance testing (BaTS) was applied, clearly supporting that the viruses evolved locally in Africa instead of spatial mixing among different regions. Overall, the results demonstrated that multiple mechanisms such as reassortment events, various mutations and possibly interspecies transmission account for the enormous diversity of genotype P[6] strains in Africa. These findings highlight the need for continued global surveillance of rotavirus diversity.

Differing epidemiological dynamics of Chikungunya virus in the Americas during the 2014-2015 epidemic.

Chikungunya virus (CHIKV) has been detected sporadically since the 1950s and includes three distinct co-circulating genotypes. In late 2013, the Asian genotype of CHIKV was responsible for the Caribbean outbreak (CO) that rapidly became an epidemic throughout the Americas. There is a limited understanding of the molecular evolution of CHIKV in the Americas during this epidemic. We sequenced 185 complete CHIKV genomes collected mainly from Nicaragua in Central America and Florida in the United States during the 2014-2015 Caribbean/Americas epidemic. Our comprehensive phylogenetic analyses estimated the epidemic history of the Asian genotype and the recent Caribbean outbreak (CO) clade, revealed considerable genetic diversity within the CO clade, and described different epidemiological dynamics of CHIKV in the Americas. Specifically, we identified multiple introductions in both Nicaragua and Florida, with rapid local spread of viruses in Nicaragua but limited autochthonous transmission in Florida in the US. Our phylogenetic analysis also showed phylogeographic clustering of the CO clade. In addition, we identified the significant amino acid substitutions that were observed across the entire Asian genotype during its evolution and examined amino acid changes that were specific to the CO clade. Deep sequencing analysis identified specific minor variants present in clinical specimens below-consensus levels. Finally, we investigated the association between viral phylogeny and geographic/clinical metadata in Nicaragua. To date, this study represents the largest single collection of CHIKV complete genomes during the Caribbean/Americas epidemic and significantly expands our understanding of the emergence and evolution of CHIKV CO clade in the Americas.

Author Correction: Genomic analysis of a large set of currently-and historically-important human adenovirus pathogens.

The original publication of this article [1] was missing the below author that made contributions to the research and the published article.

Genomic analysis of a large set of currently-and historically-important human adenovirus pathogens.

Human adenoviruses (HAdVs) are uniquely important "model organisms" as they have been used to elucidate fundamental biological processes, are recognized as complex pathogens, and are used as remedies for human health. As pathogens, HAdVs may effect asymptomatic or mild and severe symptomatic disease upon their infection of respiratory, ocular, gastrointestinal, and genitourinary systems. High-resolution genomic data have enhanced the understanding of HAdV epidemiology, with recombination recognized as an important and major pathway in the molecular evolution and genesis of emergent HAdV pathogens. To support this view and to actualize an algorithm for identifying, characterizing, and typing novel HAdVs, we determined the DNA sequence of 95 isolates from archives containing historically important pathogens and collections housing currently circulating strains to be sequenced. Of the 85 samples that were completely sequenced, 18 novel recombinants within species HAdV-B and D were identified. Two HAdV-D genomes were found to contain novel penton base and fiber genes with significant divergence from known molecular types. In this data set, we found additional isolates of HAdV-D53 and HAdV-D58, two novel genotypes recognized recently using genomics. This supports the thesis that novel HAdV genotypes are not limited to "one-time" appearances of the prototype but are of importance in HAdV epidemiology. These data underscore the significance of lateral genomic transfer in HAdV evolution and reinforce the potential public health impact of novel genotypes of HAdVs emerging in the population.

Isolation of type A influenza viruses from Red-necked Grebes (Podiceps grisegena).

Six type-A low pathogenic influenza viruses from 14 Red-necked Grebes (Podiceps grisegena) from Agassiz National Wildlife Refuge were sequenced. The grebe viruses were closely related to North American duck viruses. The genetic and temporal subtype consistency between the duck and grebe isolates suggest spillover events, potentially enhanced by feather eating.

The contrasting phylodynamics of human influenza B viruses.

A complex interplay of viral, host, and ecological factors shapes the spatio-temporal incidence and evolution of human influenza viruses. Although considerable attention has been paid to influenza A viruses, a lack of equivalent data means that an integrated evolutionary and epidemiological framework has until now not been available for influenza B viruses, despite their significant disease burden. Through the analysis of over 900 full genomes from an epidemiological collection of more than 26,000 strains from Australia and New Zealand, we reveal fundamental differences in the phylodynamics of the two co-circulating lineages of influenza B virus (Victoria and Yamagata), showing that their individual dynamics are determined by a complex relationship between virus transmission, age of infection, and receptor binding preference. In sum, this work identifies new factors that are important determinants of influenza B evolution and epidemiology.

Genomic and protein structural maps of adaptive evolution of human influenza A virus to increased virulence in the mouse.

Adaptive evolution is characterized by positive and parallel, or repeated selection of mutations. Mouse adaptation of influenza A virus (IAV) produces virulent mutants that demonstrate positive and parallel evolution of mutations in the hemagglutinin (HA) receptor and non-structural protein 1 (NS1) interferon antagonist genes. We now present a genomic analysis of all 11 genes of 39 mouse adapted IAV variants from 10 replicate adaptation experiments. Mutations were mapped on the primary and structural maps of each protein and specific mutations were validated with respect to virulence, replication, and RNA polymerase activity. Mouse adapted (MA) variants obtained after 12 or 20-21 serial infections acquired on average 5.8 and 7.9 nonsynonymous mutations per genome of 11 genes, respectively. Among a total of 115 nonsynonymous mutations, 51 demonstrated properties of natural selection including 27 parallel mutations. The greatest degree of parallel evolution occurred in the HA receptor and ribonucleocapsid components, polymerase subunits (PB1, PB2, PA) and NP. Mutations occurred in host nuclear trafficking factor binding sites as well as sites of virus-virus protein subunit interaction for NP, NS1, HA and NA proteins. Adaptive regions included cap binding and endonuclease domains in the PB2 and PA polymerase subunits. Four mutations in NS1 resulted in loss of binding to the host cleavage and polyadenylation specificity factor (CPSF30) suggesting that a reduction in inhibition of host gene expression was being selected. The most prevalent mutations in PB2 and NP were shown to increase virulence but differed in their ability to enhance replication and demonstrated epistatic effects. Several positively selected RNA polymerase mutations demonstrated increased virulence associated with >300% enhanced polymerase activity. Adaptive mutations that control host range and virulence were identified by their repeated selection to comprise a defined model for studying IAV evolution to increased virulence in the mouse.

Sequencing Bacillus anthracis typing phages gamma and cherry reveals a common ancestry.

The genetic relatedness of the Bacillus anthracis typing phages Gamma and Cherry was determined by nucleotide sequencing and comparative analysis. The genomes of these two phages were identical except at three variable loci, which showed heterogeneity within individual lysates and among Cherry, Wbeta, Fah, and four Gamma bacteriophage sequences.