Multi-task learning sparse group lasso: a method for quantifying antigenicity of influenza A(H1N1) virus using mutations and variations in glycosylation of Hemagglutinin

MAIVeSS: streamlined selection of antigenically matched, high-yield viruses for seasonal influenza vaccine production

Vaccines are the main pharmaceutical intervention used against the global public health threat posed by influenza viruses. Timely selection of optimal seed viruses with matched antigenicity between vaccine antigen and circulating viruses and with high yield underscore vaccine efficacy and supply, respectively. Current methods for selecting influenza seed vaccines are labor intensive and time-consuming. Here, we report the Machine-learning Assisted Influenza VaccinE Strain Selection framework, MAIVeSS, that enables streamlined selection of naturally circulating, antigenically matched, and high-yield influenza vaccine strains directly from clinical samples by using molecular signatures of antigenicity and yield to support optimal candidate vaccine virus selection. We apply our framework on publicly available sequences to select A(H1N1)pdm09 vaccine candidates and experimentally confirm that these candidates have optimal antigenicity and growth in cells and eggs. Our framework can potentially reduce the optimal vaccine candidate selection time from months to days and thus facilitate timely supply of seasonal vaccines.

A sequence-based machine learning model for predicting antigenic distance for H3N2 influenza virus

Introduction

Seasonal influenza A H3N2 viruses are constantly changing, reducing the effectiveness of existing vaccines. As a result, the World Health Organization (WHO) needs to frequently update the vaccine strains to match the antigenicity of emerged H3N2 variants. Traditional assessments of antigenicity rely on serological methods, which are both labor-intensive and time-consuming. Although numerous computational models aim to simplify antigenicity determination, they either lack a robust quantitative linkage between antigenicity and viral sequences or focus restrictively on selected features.

Methods

Here, we propose a novel computational method to predict antigenic distances using multiple features, including not only viral sequence attributes but also integrating four distinct categories of features that significantly affect viral antigenicity in sequences.

Results

This method exhibits low error in virus antigenicity prediction and achieves superior accuracy in discerning antigenic drift. Utilizing this method, we investigated the evolution process of the H3N2 influenza viruses and identified a total of 21 major antigenic clusters from 1968 to 2022.

Discussion

Interestingly, our predicted antigenic map aligns closely with the antigenic map generated with serological data. Thus, our method is a promising tool for detecting antigenic variants and guiding the selection of vaccine candidates.

Development of PREDAC-H1pdm to model the antigenic evolution of influenza A/(H1N1) pdm09 viruses

The Influenza A (H1N1) pdm09 virus caused a global pandemic in 2009 and has circulated seasonally ever since. As the continual genetic evolution of hemagglutinin in this virus leads to antigenic drift, rapid identification of antigenic variants and characterization of the antigenic evolution are needed. In this study, we developed PREDAC-H1pdm, a model to predict antigenic relationships between H1N1pdm viruses and identify antigenic clusters for post-2009 pandemic H1N1 strains. Our model performed well in predicting antigenic variants, which was helpful in influenza surveillance. By mapping the antigenic clusters for H1N1pdm, we found that substitutions on the Sa epitope were common for H1N1pdm, whereas for the former seasonal H1N1, substitutions on the Sb epitope were more common in antigenic evolution. Additionally, the localized epidemic pattern of H1N1pdm was more obvious than that of the former seasonal H1N1, which could make vaccine recommendation more sophisticated. Overall, the antigenic relationship prediction model we developed provides a rapid determination method for identifying antigenic variants, and the further analysis of evolutionary and epidemic characteristics can facilitate vaccine recommendations and influenza surveillance for H1N1pdm.

Influenza vaccination-induced H3 stalk-reactive memory B-cell clone expansion

Current vaccine formulations elicit a recall immune response against viruses by targeting epitopes on the globular head of hemagglutinin (HA), and stalk-reactive antibodies are rarely found. However, stalk-specific memory B-cell expansion after influenza vaccination is poorly understood. In this study, B cells were isolated from individuals immunized with seasonal tetravalent influenza vaccines at days 0 and 28 for H7N9 stimulation in vitro. Plasma and supernatants were collected for the analysis of anti-HA IgG using ELISA and a Luminex assay. Memory B cells were positively enriched, and total RNA was extracted for B cell receptor (BCR) H-CDR3 sequencing. All subjects displayed increased anti-H3 antibody secretion after vaccination, whereas no increase in cH5/3-reactive IgG levels was detected. The number of shared memory B-cell clones among individuals dropped dramatically from 593 to 37. Four out of 5 subjects displayed enhanced frequencies of the VH3-23 and VH3-30 genes, and one exhibited an increase in the frequency of VH1-18, which are associated with the stalk of HA. An increase in H3 stalk-specific antibodies produced by B cells stimulated with H7N9 viruses was detected after vaccination. These results demonstrated that H3 stalk-specific memory B cells can expand and secrete antibodies that bind to the stalk in vitro, although no increase in serum H3 stalk-reactive antibodies was found after vaccination, indicating potential for developing a universal vaccine strategy.

Selecting Representative Samples From Complex Biological Datasets Using K-Medoids Clustering

Rapid growth of single-cell sequencing techniques enables researchers to investigate almost millions of cells with diverse properties in a single experiment. Meanwhile, it also presents great challenges for selecting representative samples from massive single-cell populations for further experimental characterization, which requires a robust and compact sampling with balancing diverse properties of different priority levels. The conventional sampling methods fail to generate representative and generalizable subsets from a massive single-cell population or more complicated ensembles. Here, we present a toolkit called Cookie which can efficiently select out the most representative samples from a massive single-cell population with diverse properties. This method quantifies the relationships/similarities among samples using their Manhattan distances by vectorizing all given properties and then determines an appropriate sample size by evaluating the coverage of key properties from multiple candidate sizes, following by a k-medoids clustering to group samples into several clusters and selects centers from each cluster as the most representatives. Comparison of Cookie with conventional sampling methods using a single-cell atlas dataset, epidemiology surveillance data, and a simulated dataset shows the high efficacy, efficiency, and flexibly of Cookie. The Cookie toolkit is implemented in R and is freely available at https://wilsonimmunologylab.github.io/Cookie/.

Antigenic characterization of influenza and SARS-CoV-2 viruses

Antigenic characterization of emerging and re-emerging viruses is necessary for the prevention of and response to outbreaks, evaluation of infection mechanisms, understanding of virus evolution, and selection of strains for vaccine development. Primary analytic methods, including enzyme-linked immunosorbent/lectin assays, hemagglutination inhibition, neuraminidase inhibition, micro-neutralization assays, and antigenic cartography, have been widely used in the field of influenza research. These techniques have been improved upon over time for increased analytical capacity, and some have been mobilized for the rapid characterization of the SARS-CoV-2 virus as well as its variants, facilitating the development of highly effective vaccines within 1 year of the initially reported outbreak. While great strides have been made for evaluating the antigenic properties of these viruses, multiple challenges prevent efficient vaccine strain selection and accurate assessment. For influenza, these barriers include the requirement for a large virus quantity to perform the assays, more than what can typically be provided by the clinical samples alone, cell- or egg-adapted mutations that can cause antigenic mismatch between the vaccine strain and circulating viruses, and up to a 6-month duration of vaccine development after vaccine strain selection, which allows viruses to continue evolving with potential for antigenic drift and, thus, antigenic mismatch between the vaccine strain and the emerging epidemic strain. SARS-CoV-2 characterization has faced similar challenges with the additional barrier of the need for facilities with high biosafety levels due to its infectious nature. In this study, we review the primary analytic methods used for antigenic characterization of influenza and SARS-CoV-2 and discuss the barriers of these methods and current developments for addressing these challenges.

Librator: a platform for the optimized analysis, design, and expression of mutable influenza viral antigens

Artificial mutagenesis and protein engineering have laid the foundation for antigenic characterization and universal vaccine design for influenza viruses. However, many methods used in this process require manual sequence editing and protein expression, limiting their efficiency and utility in high-throughput applications. More streamlined in silico tools allowing researchers to properly analyze and visualize influenza viral protein sequences with accurate nomenclature are necessary to improve antigen design and productivity. To address this need, we developed Librator, a system for analyzing and designing custom protein sequences of influenza virus hemagglutinin (HA) and neuraminidase (NA) glycoproteins. Within Librator's graphical interface, users can easily interrogate viral sequences and phylogenies, visualize antigen structures and conservation, mutate target residues and design custom antigens. Librator also provides optimized fragment design for Gibson Assembly of HA and NA expression constructs based on peptide conservation of all historical HA and NA sequences, ensuring fragments are reusable and compatible across related subtypes, thereby promoting reagent savings. Finally, the program facilitates single-cell immune profiling, epitope mapping of monoclonal antibodies and mosaic protein design. Using Librator-based antigen construction, we demonstrate that antigenicity can be readily transferred between HA molecules of H3, but not H1, lineage viruses. Altogether, Librator is a valuable tool for analyzing influenza virus HA and NA proteins and provides an efficient resource for optimizing recombinant influenza antigen synthesis.