Influence of Sex on Respiratory Syncytial Virus Genotype Infection Frequency and Nasopharyngeal Microbiome.

Respiratory syncytial virus (RSV) has a significant health burden in children, older adults, and the immunocompromised. However, limited effort has been made to identify emergence of new RSV genotypes' frequency of infection and how the combination of nasopharyngeal microbiome and viral genotypes impact RSV disease outcomes. In an observational cohort designed to capture the first infant RSV infection, we employed multi-omics approaches to sequence 349 RSV complete genomes and matched nasopharyngeal microbiomes, during which the 2012/2013 season was dominated by RSV-A, whereas 2013 and 2014 was dominated by RSV-B. We found non-G-72nt-duplicated RSV-A strains were more frequent in male infants (P = 0.02), whereas G-72nt-duplicated genotypes (which is ON1 lineage) were seen equally in both males and females. DESeq2 testing of the nasal microbiome showed Haemophilus was significantly more abundant in infants with RSV-A infection compared to infants with RSV-B infection (adjusted P = 0.002). In addition, the broad microbial clustering of the abundant genera was significantly associated with infant sex (P = 0.03). Overall, we show sex differences in infection by RSV genotype and host nasopharyngeal microbiome, suggesting an interaction between host genetics, virus genotype, and associated nasopharyngeal microbiome. IMPORTANCE Respiratory syncytial virus (RSV) is one of the leading causes of lower respiratory tract infections in young children and is responsible for high hospitalization rates and morbidity in infants and the elderly. To understand how the emergence of RSV viral genotypes and viral-respiratory microbiome interactions contribute to infection frequency and severity, we utilized an observational cohort designed to capture the first infant RSV infection we employed multi-omics approaches to sequence 349 RSV complete genomes and matched nasopharyngeal microbiomes. We found non-G-72nt-duplicated RSV-A genotypes were more frequent in male infants, whereas G-72nt-duplicated RSV-A strains (ON1 lineage) were seen equally in both males and females. Microbiome analysis show Haemophilus was significantly more abundant in infants with RSV-A compared to infants with RSV-B infection and the microbial clustering of the abundant genera was associated with infant sex. Overall, we show sex differences in RSV genotype-nasopharyngeal microbiome, suggesting an interaction host genetics-virus-microbiome interaction.

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.

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.].

Sequencing and analysis of globally obtained human parainfluenza viruses 1 and 3 genomes.

Human Parainfluenza viruses (HPIV) type 1 and 3 are important causes of respiratory tract infections in young children globally. HPIV infections do not confer complete protective immunity so reinfections occur throughout life. Since no effective vaccine is available for the two virus subtypes, comprehensive understanding of HPIV-1 and HPIV-3 genetic and epidemic features is important for diagnosis, prevention, and treatment of HPIV-1 and HPIV-3 infections. Relatively few whole genome sequences are available for both HPIV-1 and HPIV-3 viruses, so our study sought to provide whole genome sequences from multiple countries to further the understanding of the global diversity of HPIV at a whole-genome level. We collected HPIV-1 and HPIV-3 samples and isolates from Argentina, Australia, France, Mexico, South Africa, Switzerland, and USA from the years 2003-2011 and sequenced the genomes of 40 HPIV-1 and 75 HPIV-3 viruses with Sanger and next-generation sequencing with the Ion Torrent, Illumina, and 454 platforms. Phylogenetic analysis showed that the HPIV-1 genome is evolving at an estimated rate of 4.97 × 10-4 mutations/site/year (95% highest posterior density 4.55 × 10-4 to 5.38 × 10-4) and the HPIV-3 genome is evolving at a similar rate (3.59 × 10-4 mutations/site/year, 95% highest posterior density 3.26 × 10-4 to 3.94 × 10-4). There were multiple genetically distinct lineages of both HPIV-1 and 3 circulating on a global scale. Further surveillance and whole-genome sequencing are greatly needed to better understand the spatial dynamics of these important respiratory viruses in humans.

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.

Human Monoclonal Antibodies Potently Neutralize Zika Virus and Select for Escape Mutations on the Lateral Ridge of the Envelope Protein.

The mosquito-borne Zika virus (ZIKV) has been causing epidemic outbreaks on a global scale. Virus infection can result in severe disease in humans, including microcephaly in newborns and Guillain-Barré syndrome in adults. Here, we characterized monoclonal antibodies isolated from a patient with an active Zika virus infection that potently neutralized virus infection in Vero cells at the nanogram-per-milliliter range. In addition, these antibodies enhanced internalization of virions into human leukemia K562 cells in vitro, indicating their possible ability to cause antibody-dependent enhancement of disease. Escape variants of the ZIKV MR766 strain to a potently neutralizing antibody, AC10, exhibited an amino acid substitution at residue S368 in the lateral ridge region of the envelope protein. Analysis of publicly availably ZIKV sequences revealed the S368 site to be conserved among the vast majority (97.6%) of circulating strains. We validated the importance of this residue by engineering a recombinant virus with an S368R point mutation that was unable to be fully neutralized by AC10. Four out of the 12 monoclonal antibodies tested were also unable to neutralize the virus with the S368R mutation, suggesting this region to be an important immunogenic epitope during human infection. Last, a time-of-addition infection assay further validated the escape variant and showed that all monoclonal antibodies inhibited virus binding to the cell surface. Thus, the present study demonstrates that the lateral ridge region of the envelope protein is likely an immunodominant, neutralizing epitope.IMPORTANCE Zika virus (ZIKV) is a global health threat causing severe disease in humans, including microcephaly in newborns and Guillain-Barré syndrome in adults. Here, we analyzed the human monoclonal antibody response to acute ZIKV infection and found that neutralizing antibodies could not elicit Fc-mediated immune effector functions but could potentiate antibody-dependent enhancement of disease. We further identified critical epitopes involved with neutralization by generating and characterizing escape variants by whole-genome sequencing. We demonstrate that the lateral ridge region, particularly the S368 amino acid site, is critical for neutralization by domain III-specific antibodies.

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.

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.

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.

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.

Sequences of Zika Virus Genomes from a Pediatric Cohort in Nicaragua.

We report here the whole-genome sequence of 11 Zika virus (ZIKV) samples from six pediatric patients in Nicaragua. Serum samples were collected, and ZIKV was isolated in tissue culture. Both serum and virus isolates were sequenced. The consensus ZIKV genomes are greater than 99% identical to each other.

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.

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.

Large-Scale Complete-Genome Sequencing and Phylodynamic Analysis of Eastern Equine Encephalitis Virus Reveals Source-Sink Transmission Dynamics in the United States.

Eastern equine encephalitis virus (EEEV) has a high case-fatality rate in horses and humans, and Florida has been hypothesized to be the source of EEEV epidemics for the northeastern United States. To test this hypothesis, we sequenced complete genomes of 433 EEEV strains collected within the United States from 1934 to 2014. Phylogenetic analysis suggested EEEV evolves relatively slowly and that transmission is enzootic in Florida, characterized by higher genetic diversity and long-term local persistence. In contrast, EEEV strains in New York and Massachusetts were characterized by lower genetic diversity, multiple introductions, and shorter local persistence. Our phylogeographic analysis supported a source-sink model in which Florida is the major source of EEEV compared to the other localities sampled. In sum, this study revealed the complex epidemiological dynamics of EEEV in different geographic regions in the United States and provided general insights into the evolution and transmission of other avian mosquito-borne viruses in this region.IMPORTANCE Eastern equine encephalitis virus (EEEV) infections are severe in horses and humans on the east coast of the United States with a >90% mortality rate in horses, an ∼33% mortality rate in humans, and significant brain damage in most human survivors. However, little is known about the evolutionary characteristics of EEEV due to the lack of genome sequences. By generating large collection of publicly available complete genome sequences, this study comprehensively determined the evolution of the virus, described the epidemiological dynamics of EEEV in different states in the United States, and identified Florida as one of the major sources. These results may have important implications for the control and prevention of other mosquito-borne viruses in the Americas.

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.

Whole-Genome Sequences of Zika Virus FLR Strains after Passage in Vero or C6/36 Cells.

We report 26 complete genomes of Zika virus (ZIKV) isolated after passaging the Zika virus strain FLR in mosquito (C6/36) and mammalian (Vero) cell lines. The consensus ZIKV genomes we recovered show greater than 99% nucleotide identify with each other and with the FLR strain used as input.

Identification of Dezidougou Virus in a DAK AR 41524 Zika Virus Stock.

We report here the complete genome of a Dezidougou virus (DEZV) isolated from a passaged culture of the Zika virus strain DAK AR 41524. The consensus DEZV sequence we recovered shows 99% nucleotide similarity using BLASTN to a previously reported DEZV (accession no. JQ675604.1). The current sequence has additional repeat regions as well as a deleted repeat region, which we confirmed by Sanger sequencing, that were not present in the originally published sequence, JQ675604.1.

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.