Diversifying Selection Analysis Predicts Antigenic Evolution of 2009 Pandemic H1N1 Influenza A Virus in Humans

Replication in the presence of dengue convalescent serum impacts Zika virus neutralization sensitivity and fitness

Introduction

Flaviviruses like dengue virus (DENV) and Zika virus (ZIKV) are mosquito-borne viruses that cause febrile, hemorrhagic, and neurological diseases in humans, resulting in 400 million infections annually. Due to their co-circulation in many parts of the world, flaviviruses must replicate in the presence of pre-existing adaptive immune responses targeted at serologically closely related pathogens, which can provide protection or enhance disease. However, the impact of pre-existing cross-reactive immunity as a driver of flavivirus evolution, and subsequently the implications on the emergence of immune escape variants, is poorly understood. Therefore, we investigated how replication in the presence of convalescent dengue serum drives ZIKV evolution.

Methods

We used an in vitro directed evolution system, passaging ZIKV in the presence of serum from humans previously infected with DENV (anti-DENV) or serum from DENV-naïve patients (control serum). Following five passages in the presence of serum, we performed next-generation sequencing to identify mutations that arose during passaging. We studied two non-synonymous mutations found in the anti-DENV passaged population (E-V355I and NS1-T139A) by generating individual ZIKV mutants and assessing fitness in mammalian cells and live mosquitoes, as well as their sensitivity to antibody neutralization.

Results and discussion

Both viruses had increased fitness in Vero cells with and without the addition of anti-DENV serum and in human lung epithelial and monocyte cells. In Aedes aegypti mosquitoes-using blood meals with and without anti-DENV serum-the mutant viruses had significantly reduced fitness compared to wild-type ZIKV. These results align with the trade-off hypothesis of constrained mosquito-borne virus evolution. Notably, only the NS1-T139A mutation escaped neutralization, while E-V335I demonstrated enhanced neutralization sensitivity to neutralization by anti-DENV serum, indicating that neutralization escape is not necessary for viruses passaged under cross-reactive immune pressures. Future studies are needed to assess cross-reactive immune selection in humans and relevant animal models or with different flaviviruses.

Characterizing Complex Populations of Endogenous Adeno-Associated Viruses by Single-Genome Amplification

The isolation of adeno-associated virus (AAV) genomes from biomaterials at the molecular level has traditionally relied on polymerase chain reaction-based and cloning-based techniques. However, when applied to samples containing multiple species, traditional techniques for isolating viral genomes can amplify artificial recombinants and introduce polymerase misincorporation errors. In this study, we describe AAV single-genome amplification (AAV-SGA): a powerful technique to isolate, amplify, and sequence single AAV genomes from mammalian genomic DNA, which can then be used to construct vectors for gene therapy. We used AAV-SGA to precisely isolate 15 novel AAV genomes belonging to AAV clades A, D, and E and the Fringe outgroup. This technique also enables investigations of AAV population dynamics and recombination events to provide insights into virus-host interactions and virus biology. Using AAV-SGA, we identified regional heterogeneity within AAV populations from different lobes of the liver of a rhesus macaque and found evidence of frequent genomic recombination between AAV populations. This study highlights the strengths of AAV-SGA and demonstrates its capability to provide valuable insights into the biology and diversity of AAVs.

Mutation in Hemagglutinin Antigenic Sites in Influenza A pH1N1 Viruses from 2015-2019 in the United States Mountain West, Europe, and the Northern Hemisphere

H1N1 influenza A virus is a respiratory pathogen that undergoes antigenic shift and antigenic drift to improve viral fitness. Tracking the evolutionary trends of H1N1 aids with the current detection and the future response to new viral strains as they emerge. Here, we characterize antigenic drift events observed in the hemagglutinin (HA) sequence of the pandemic H1N1 lineage from 2015-2019. We observed the substitutions S200P, K147N, and P154S, together with other mutations in structural, functional, and/or epitope regions in 2015-2019 HA protein sequences from the Mountain West region of the United States, the larger United States, Europe, and other Northern Hemisphere countries. We reconstructed multiple phylogenetic trees to track the relationships and spread of these mutations and tested for evidence of selection pressure on HA. We found that the prevalence of amino acid substitutions at positions 147, 154, 159, 200, and 233 significantly changed throughout the studied geographical regions between 2015 and 2019. We also found evidence of coevolution among a subset of these amino acid substitutions. The results from this study could be relevant for future epidemiological tracking and vaccine prediction efforts. Similar analyses in the future could identify additional sequence changes that could affect the pathogenicity and/or infectivity of this virus in its human host.

Isolating Natural Adeno-Associated Viruses from Primate Tissues with a High-Fidelity Polymerase

Adeno-associated viruses (AAVs) are advantageous as gene-transfer vectors due to their favorable biological and safety characteristics, with discovering novel AAV variants being key to improving this treatment platform. To date, researchers have isolated over 200 AAVs from natural sources using PCR-based methods. We compared two modern DNA polymerases and their utility for isolating and amplifying the AAV genome. Compared to the HotStar polymerase, the higher-fidelity Q5 Hot Start High-Fidelity DNA Polymerase provided more precise and accurate amplification of the input AAV sequences. The lower-fidelity HotStar DNA polymerase introduced mutations during the isolation and amplification processes, thus generating multiple mutant capsids with variable bioactivity compared to the input AAV gene. The Q5 polymerase enabled the successful discovery of novel AAV capsid sequences from human and nonhuman primate tissue sources. Novel AAV sequences from these sources showed evidence of positive evolutionary selection. This study highlights the importance of using the highest fidelity DNA polymerases available to accurately isolate and characterize AAV genomes from natural sources to ultimately develop more effective gene therapy vectors.

Prediction model for the spread of the COVID-19 outbreak in the global environment

COVID-19 has long become a worldwide pandemic. It is responsible for the death of over two million people and posed an economic recession. This paper studies the spread pattern of COVID-19, aiming to establish a prediction model for this event. We harness Data Mining and Machine Learning methodologies to train regression models to predict the number of confirmed cases in a spatial-temporal space. We introduce an innovative concept ‒ the Center of Infection Mass (CoIM) ‒ adapted from the field of physics. We empirically evaluated our model on western European countries, based on the CoIM index and other features, and showed that a relatively high accurate prediction of the spread can be obtained. Our contribution is twofold: first, we introduced a prediction methodology and proved empirically that a prediction can be made even to the range of over a month; second, we showed promise in adopting the CoIM index to prediction models, when models that adopt the CoIM yield significantly better results than those that discard it. By applying our model, and better controlling the inherent tradeoff between life-saving and economy, we believe that decision-makers can take close to optimal measures. Thus, this methodology may contribute to public welfare.

Database and Analytical Resources for Viral Research Community

Many public databases and analytical resources are available to facilitate virology research. The Virus Pathogen Database and Analysis Resource (ViPR, see “Relevant Websites section”) and Influenza Research Database (IRD, see “Relevant Websites section”) are comprehensive and highly curated repositories of genome and protein sequence records and annotations, protein structures, immune epitopes, and epidemiological and surveillance data about human and related viral pathogens. These data are acquired from public repositories, direct submissions and in-house bioinformatics analyses. The resources offer seamless integration of data, analytics and visualization, and are freely available without cost or restriction to facilitate diagnostics and therapeutics development for priority pathogens.

Interim estimates of 2018/19 vaccine effectiveness against influenza A(H1N1)pdm09, Canada, January 2019

Using a test-negative design, the Canadian Sentinel Practitioner Surveillance Network assessed interim 2018/19 vaccine effectiveness (VE) against predominant influenza A(H1N1)pdm09 viruses. Adjusted VE was 72% (95% confidence interval: 60 to 81) against medically attended, laboratory-confirmed influenza A(H1N1)pdm09 illness. This substantial vaccine protection was observed in all age groups, notably young children who appeared to be disproportionately affected. Sequence analysis identified heterogeneity in emerging clade 6B.1 viruses but no dominant drift variant.

Clinical Performance of the AllplexTM Respiratory Panel 1 Test Compared to SimplexaTM Flu A/B and RSV for Detection of Influenza Virus and Respiratory Syncytial Virus Infection Including Their Subtyping

OBJECTIVE

TheAllplexTM Respiratory Panel 1 (ARP) is a new assay based on a real-time polymerase chain reaction (RT-PCR) for the detection of influenza A (Flu A), influenza B virus (Flu B), and respiratory syncytial virus (RSV), including subtyping by multiple detection temperature (MuDT) technology. We evaluated the performance of the Allplex Respiratory Panel compared to the SimplexaTM Flu A/B & RSV assay (SP) and other diagnostic tools.

MATERIALS AND METHODS

A total of 372 samples were collected from patients at the Korea University Guro Hospital in Seoul, Korea. All samples were tested for influenza virus and RSV by ARP, SP, and an in-house RT-PCR.

RESULTS

The sensitivity of ARP was 95.56, 100, and 95.24% for Flu A, Flu B, and RSV, respectively. The specificity of ARP was 100, 100, and 100% for Flu A, Flu B, and RSV, respectively. SP had sensitivities and specificities of 98.89 and 100% for Flu A, 100 and 100% for Flu B, and 100 and 100% for RSV.

CONCLUSION

The Allplex panelshowed high sensitivity, specificity, positive predictive, and negative predictive values for the detection of Flu A, Flu B, and RSV. This assay is fast and easy to perform because it takes only about 150 min and there is no need for post-PCR electrophoresis. The ARP can be used as a reliable and convenient assay in clinical laboratories.

The influenza virus hemagglutinin head evolves faster than the stalk domain

The limited ability of current influenza virus vaccines to protect from antigenically drifted or shifted viruses creates a public health problem that has led to the need to develop effective, broadly protective vaccines. While current influenza virus vaccines mostly induce an immune response against the immunodominant and variable head domain of the hemagglutinin, the major surface glycoprotein of the virus, the hemagglutinin stalk domain has been identified to harbor neutralizing B-cell epitopes that are conserved among and even between influenza A virus subtypes. A complete understanding of the differences in evolution between the main target of current vaccines and this more conserved stalk region are missing. Here, we performed an evolutionary analysis of the stalk domains of the hemagglutinin of pre-pandemic seasonal H1N1, pandemic H1N1, seasonal H3N2, and influenza B viruses and show quantitatively for the first time that the stalk domain is evolving at a rate that is significantly slower than that of the head domain. Additionally, we found that the cross-reactive epitopes in the stalk domain targeted by broadly neutralizing monoclonal antibodies are evolving at an even slower rate compared to the full head and stalk regions of the protein. Finally, a fixed-effects likelihood selection analysis was performed for these virus groups in both the head and stalk domains. While several positive selection sites were found in the head domain, only a single site in the stalk domain of pre-pandemic seasonal H1 hemagglutinin was identified at amino acid position 468 (H1 numbering from methionine). This site is not located in or close to the epitopes of cross-reactive anti-stalk monoclonal antibodies. Furthermore, we found that changes in this site do not significantly impact virus binding or neutralization by human anti-stalk antibodies, suggesting that some positive selection in the stalk domain is independent of immune pressures. We conclude that, while the stalk domain does evolve over time, this evolution is slow and, historically, is not directed to aid in evading neutralizing antibody responses.

Applying Unique Molecular Identifiers in Next Generation Sequencing Reveals a Constrained Viral Quasispecies Evolution under Cross-Reactive Antibody Pressure Targeting Long Alpha Helix of Hemagglutinin

To overcome yearly efforts and costs for the production of seasonal influenza vaccines, new approaches for the induction of broadly protective and long-lasting immune responses have been developed in the past decade. To warrant safety and efficacy of the emerging crossreactive vaccine candidates, it is critical to understand the evolution of influenza viruses in response to these new immune pressures. Here we applied unique molecular identifiers in next generation sequencing to analyze the evolution of influenza quasispecies under in vivo antibody pressure targeting the hemagglutinin (HA) long alpha helix (LAH). Our vaccine targeting LAH of hemagglutinin elicited significant seroconversion and protection against homologous and heterologous influenza virus strains in mice. The vaccine not only significantly reduced lung viral titers, but also induced a well-known bottleneck effect by decreasing virus diversity. In contrast to the classical bottleneck effect, here we showed a significant increase in the frequency of viruses with amino acid sequences identical to that of vaccine targeting LAH domain. No escape mutant emerged after vaccination. These results not only support the potential of a universal influenza vaccine targeting the conserved LAH domains, but also clearly demonstrate that the well-established bottleneck effect on viral quasispecies evolution does not necessarily generate escape mutants.

Antigenic cartography of H1N1 influenza viruses using sequence-based antigenic distance calculation

BACKGROUND

The ease at which influenza virus sequence data can be used to estimate antigenic relationships between strains and the existence of databases containing sequence data for hundreds of thousands influenza strains make sequence-based antigenic distance estimates an attractive approach to researchers. Antigenic mismatch between circulating strains and vaccine strains results in significantly decreased vaccine effectiveness. Furthermore, antigenic relatedness between the vaccine strain and the strains an individual was originally primed with can affect the cross-reactivity of the antibody response. Thus, understanding the antigenic relationships between influenza viruses that have circulated is important to both vaccinologists and immunologists.

RESULTS

Here we develop a method of mapping antigenic relationships between influenza virus stains using a sequence-based antigenic distance approach (SBM). We used a modified version of the p-all-epitope sequence-based antigenic distance calculation, which determines the antigenic relatedness between strains using influenza hemagglutinin (HA) genetic coding sequence data and provide experimental validation of the p-all-epitope calculation. We calculated the antigenic distance between 4838 H1N1 viruses isolated from infected humans between 1918 and 2016. We demonstrate, for the first time, that sequence-based antigenic distances of H1N1 Influenza viruses can be accurately represented in 2-dimenstional antigenic cartography using classic multidimensional scaling. Additionally, the model correctly predicted decreases in cross-reactive antibody levels with 87% accuracy and was highly reproducible with even when small numbers of sequences were used.

CONCLUSION

This work provides a highly accurate and precise bioinformatics tool that can be used to assess immune risk as well as design optimized vaccination strategies. SBM accurately estimated the antigenic relationship between strains using HA sequence data. Antigenic maps of H1N1 virus strains reveal that strains cluster antigenically similar to what has been reported for H3N2 viruses. Furthermore, we demonstrated that genetic variation differs across antigenic sites and discuss the implications.

Temporal distribution and genetic variants in influenza A(H1N1)pdm09 virus circulating in Mexico, seasons 2012 and 2013

The 2012 and 2013 annual influenza epidemics in Mexico were characterized by presenting different seasonal patterns. In 2012 the A(H1N1)pdm09 virus caused a high incidence of influenza infections after a two-year period of low circulation; whereas the 2013 epidemic presented circulation of the A(H1N1)pdm09 virus throughout the year. We have characterized the molecular composition of the Hemagglutinin (HA) and Neuraminidase (NA) genes of the A(H1N1)pdm09 virus from both epidemic seasons, emphasizing the genetic characteristics of viruses isolated from Yucatan in Southern Mexico. The molecular analysis of viruses from the 2012 revealed that all viruses from Mexico were predominantly grouped in clade 7. Strikingly, the molecular characterization of viruses from 2013 revealed that viruses circulating in Yucatan were genetically different to viruses from other regions of Mexico. In fact, we identified the occurrence of two genetic variants containing relevant mutations at both the HA and NA surface antigens. There was a difference on the temporal circulation of each genetic variant, viruses containing the mutations HA-A141T / NA-N341S were detected in May, June and July; whereas viruses containing the mutations HA-S162I / NA-L206S circulated in August and September. We discuss the significance of these novel genetic changes.

Antigenicity of the 2015-2016 seasonal H1N1 human influenza virus HA and NA proteins

Antigenic drift of the hemagglutinin (HA) and neuraminidase (NA) influenza virus proteins contributes to reduced vaccine efficacy. To analyze antigenic drift in human seasonal H1N1 viruses derived from the 2009 pandemic H1N1 virus (pH1N1-like viruses) accounts for the limited effectiveness (around 40%) of vaccination against pH1N1-like viruses during the 2015-2016 season, nasal washes/swabs collected from adult subjects in the Rochester, NY area, were used to sequence and isolate the circulating viruses. The HA and NA proteins from viruses circulating during the 2015-2016 season encoded eighteen and fourteen amino acid differences, respectively, when compared to A/California/04/2009, a strain circulating at the origin of the 2009 pandemic. The circulating strains belonged to subclade 6B.1, defined by HA amino acid substitutions S101N, S179N, and I233T. Hemagglutination-inhibiting (HAI) and HA-specific neutralizing serum antibody (Ab) titers from around 50% of pH1N1-like virus-infected subjects and immune ferrets were 2-4 fold lower for the 2015-2016 circulating strains compared to the vaccine strain. In addition, using a luminex-based mPlex HA assay, the binding of human sera from subjects infected with pH1N1-like viruses to the HA proteins from circulating and vaccine strains was not identical, strongly suggesting antigenic differences in the HA protein. Additionally, NA inhibition (NAI) Ab titers in human sera from pH1N1-like virus-infected subjects increased after the infection and there were measurable antigenic differences between the NA protein of circulating strains and the vaccine strain using both ferret and human antisera. Despite having been vaccinated, infected subjects exhibited low HAI Ab titers against the vaccine and circulating strains. This suggests that poor responses to the H1N1 component of the vaccine as well as antigenic differences in the HA and NA proteins of currently circulating pH1N1-like viruses could be contributing to risk of infection even after vaccination.

Hajj vaccinations-facts, challenges, and hope

Vaccination is an effective preventive measure that has been used in the unique Hajj pilgrimage setting to control the transmission of infectious diseases. The current vaccination policy applied during Hajj is reviewed herein, highlighting the effectiveness of the approaches applied and identifying research gaps that need to be filled in order to improve the development and dissemination of Hajj vaccination strategies.

Genome-Wide Analysis of Evolutionary Markers of Human Influenza A(H1N1)pdm09 and A(H3N2) Viruses May Guide Selection of Vaccine Strain Candidates

Here we analyzed whole-genome sequences of 3,969 influenza A(H1N1)pdm09 and 4,774 A(H3N2) strains that circulated during 2009–2015 in the world. The analysis revealed changes at 481 and 533 amino acid sites in proteins of influenza A(H1N1)pdm09 and A(H3N2) strains, respectively. Many of these changes were introduced as a result of random drift. However, there were 61 and 68 changes that were present in relatively large number of A(H1N1)pdm09 and A(H3N2) strains, respectively, that circulated during relatively long time. We named these amino acid substitutions evolutionary markers, as they seemed to contain valuable information regarding the viral evolution. Interestingly, influenza A(H1N1)pdm09 and A(H3N2) viruses acquired non-overlapping sets of evolutionary markers. We next analyzed these characteristic markers in vaccine strains recommended by the World Health Organization for the past five years. Our analysis revealed that vaccine strains carried only few evolutionary markers at antigenic sites of viral hemagglutinin (HA) and neuraminidase (NA). The absence of these markers at antigenic sites could affect the recognition of HA and NA by human antibodies generated in response to vaccinations. This could, in part, explain moderate efficacy of influenza vaccines during 2009–2014. Finally, we identified influenza A(H1N1)pdm09 and A(H3N2) strains, which contain all the evolutionary markers of influenza A strains circulated in 2015, and which could be used as vaccine candidates for the 2015/2016 season. Thus, genome-wide analysis of evolutionary markers of influenza A(H1N1)pdm09 and A(H3N2) viruses may guide selection of vaccine strain candidates.

Phylodynamics of H1N1/2009 influenza reveals the transition from host adaptation to immune-driven selection

Influenza A H1N1/2009 virus that emerged from swine rapidly replaced the previous seasonal H1N1 virus. Although the early emergence and diversification of H1N1/2009 is well characterized, the ongoing evolutionary and global transmission dynamics of the virus remain poorly investigated. To address this we analyse >3,000 H1N1/2009 genomes, including 214 full genomes generated from our surveillance in Singapore, in conjunction with antigenic data. Here we show that natural selection acting on H1N1/2009 directly after introduction into humans was driven by adaptation to the new host. Since then, selection has been driven by immunological escape, with these changes corresponding to restricted antigenic diversity in the virus population. We also show that H1N1/2009 viruses have been subject to regular seasonal bottlenecks and a global reduction in antigenic and genetic diversity in 2014.