Antigenic cartography using variant-specific hamster sera reveals substantial antigenic variation among Omicron subvariants.

Severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) has developed substantial antigenic variability. As the majority of the population now has pre-existing immunity due to infection or vaccination, the use of experimentally generated animal immune sera can be valuable for measuring antigenic differences between virus variants. Here, we immunized Syrian hamsters by two successive infections with one of nine SARS-CoV-2 variants. Their sera were titrated against 16 SARS-CoV-2 variants, and the resulting titers were visualized using antigenic cartography. The antigenic map shows a condensed cluster containing all pre-Omicron variants (D614G, Alpha, Delta, Beta, Mu, and an engineered B.1+E484K variant) and considerably more diversity among a selected panel of Omicron subvariants (BA.1, BA.2, BA.4/BA.5, the BA.5 descendants BF.7 and BQ.1.18, the BA.2.75 descendant BN.1.3.1, the BA.2-derived recombinants XBB.2 and EG.5.1, and the BA.2.86 descendant JN.1). Some Omicron subvariants were as antigenically distinct from each other as the wildtype is from the Omicron BA.1 variant. Compared to titers measured in human sera, titers in hamster sera are of higher magnitude, show less fold change, and result in a more compact antigenic map topology. The results highlight the potential of sera from hamsters for the continued antigenic characterization of SARS-CoV-2.

Substitutions near the HA receptor binding site explain the origin and major antigenic change of the B/Victoria and B/Yamagata lineages.

Influenza B virus primarily infects humans, causing seasonal epidemics globally. Two antigenic variants-Victoria-like and Yamagata-like-were detected in the 1980s, of which the molecular basis of emergence is still incompletely understood. Here, the antigenic properties of a unique collection of historical virus isolates, sampled from 1962 to 2000 and passaged exclusively in mammalian cells to preserve antigenic properties, were determined with the hemagglutination inhibition assay and an antigenic map was built to quantify and visualize the divergence of the lineages. The antigenic map revealed only three distinct antigenic clusters-Early, Victoria, and Yamagata-with relatively little antigenic diversity in each cluster until 2000. Viruses with Victoria-like antigenic properties emerged around 1972 and diversified subsequently into two genetic lineages. Viruses with Yamagata-like antigenic properties evolved from one lineage and became clearly antigenically distinct from the Victoria-like viruses around 1988. Recombinant mutant viruses were tested to show that insertions and deletions (indels), as observed frequently in influenza B virus hemagglutinin, had little effect on antigenic properties. In contrast, amino-acid substitutions at positions 148, 149, 150, and 203, adjacent to the hemagglutinin receptor binding site, determined the main antigenic differences between the Early, Victoria-like, and Yamagata-like viruses. Surprisingly, substitutions at two of the four positions reverted in recent viruses of the Victoria lineage, resulting in antigenic properties similar to viruses circulating ∼50 y earlier. These data shed light on the antigenic diversification of influenza viruses and suggest there may be limits to the antigenic evolution of influenza B virus.

Durability of Cross-Neutralizing Antibodies 5.5 Months After Bivalent Coronavirus Disease 2019 Vaccine Booster.

We analyzed neutralizing antibodies in samples from ancestral + BA.1 and ancestral + BA.4/5 boosted individuals, collected around 5.5 months after booster. Titers of neutralizing antibodies generally decreased compared to a time point early after the bivalent booster immunization. This was more pronounced for individuals without infection history and for recently emerged Omicron variants.

Beneath the surface: Amino acid variation underlying two decades of dengue virus antigenic dynamics in Bangkok, Thailand.

Neutralizing antibodies are important correlates of protection against dengue. Yet, determinants of variation in neutralization across strains within the four dengue virus serotypes (DENV1-4) is imperfectly understood. Studies focus on structural DENV proteins, especially the envelope (E), the primary target of anti-DENV antibodies. Although changes in immune recognition (antigenicity) are often attributed to variation in epitope residues, viral processes influencing conformation and epitope accessibility also affect neutralizability, suggesting possible modulating roles of nonstructural proteins. We estimated effects of residue changes in all 10 DENV proteins on antigenic distances between 348 DENV collected from individuals living in Bangkok, Thailand (1994-2014). Antigenic distances were derived from response of each virus to a panel of twenty non-human primate antisera. Across 100 estimations, excluding 10% of virus pairs each time, 77 of 295 positions with residue variability in E consistently conferred antigenic effects; 52 were within ±3 sites of known binding sites of neutralizing human monoclonal antibodies, exceeding expectations from random assignments of effects to sites (p = 0.037). Effects were also identified for 16 sites on the stem/anchor of E which were only recently shown to become exposed under physiological conditions. For all proteins, except nonstructural protein 2A (NS2A), root-mean-squared-error (RMSE) in predicting distances between pairs held out in each estimation did not outperform sequences of equal length derived from all proteins or E, suggesting that antigenic signals present were likely through linkage with E. Adjusted for E, we identified 62/219 sites embedding the excess signals in NS2A. Concatenating these sites to E additionally explained 3.4% to 4.0% of observed variance in antigenic distances compared to E alone (50.5% to 50.8%); RMSE outperformed concatenating E with sites from any protein of the virus (ΔRMSE, 95%IQR: 0.01, 0.05). Our results support examining antigenic determinants beyond the DENV surface.

Vaccination with SARS-CoV-2 variants of concern protects mice from challenge with wild-type virus.

Vaccines against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) have been highly efficient in protecting against Coronavirus Disease 2019 (COVID-19). However, the emergence of viral variants that are more transmissible and, in some cases, escape from neutralizing antibody responses has raised concerns. Here, we evaluated recombinant protein spike antigens derived from wild-type SARS-CoV-2 and from variants B.1.1.7, B.1.351, and P.1 for their immunogenicity and protective effect in vivo against challenge with wild-type SARS-CoV-2 in the mouse model. All proteins induced high neutralizing antibodies against the respective viruses but also induced high cross-neutralizing antibody responses. The decline in neutralizing titers between variants was moderate, with B.1.1.7-vaccinated animals having a maximum fold reduction of 4.8 against B.1.351 virus. P.1 induced the most cross-reactive antibody responses but was also the least immunogenic in terms of homologous neutralization titers. However, all antigens protected from challenge with wild-type SARS-CoV-2 in a mouse model.

Author Correction: Ancient hepatitis B viruses from the Bronze Age to the Medieval period.

In Fig. 2 of this Letter, the 'E' and 'G' clade labels were inadvertently reversed, and in Table 2 the genotype of DA27 was 'D1' instead of 'D5'. These have been corrected online.

Ancient hepatitis B viruses from the Bronze Age to the Medieval period.

Hepatitis B virus (HBV) is a major cause of human hepatitis. There is considerable uncertainty about the timescale of its evolution and its association with humans. Here we present 12 full or partial ancient HBV genomes that are between approximately 0.8 and 4.5 thousand years old. The ancient sequences group either within or in a sister relationship with extant human or other ape HBV clades. Generally, the genome properties follow those of modern HBV. The root of the HBV tree is projected to between 8.6 and 20.9 thousand years ago, and we estimate a substitution rate of 8.04 × 10-6-1.51 × 10-5 nucleotide substitutions per site per year. In several cases, the geographical locations of the ancient genotypes do not match present-day distributions. Genotypes that today are typical of Africa and Asia, and a subgenotype from India, are shown to have an early Eurasian presence. The geographical and temporal patterns that we observe in ancient and modern HBV genotypes are compatible with well-documented human migrations during the Bronze and Iron Ages1,2. We provide evidence for the creation of HBV genotype A via recombination, and for a long-term association of modern HBV genotypes with humans, including the discovery of a human genotype that is now extinct. These data expose a complexity of HBV evolution that is not evident when considering modern sequences alone.

Immunogenicity of a 2-Dose Regimen of Moderna mRNA Beta/Omicron BA.1 Bivalent Variant Vaccine Boost in a Randomized Clinical Trial.

We compared the serologic responses of 1 dose versus 2 doses of a variant vaccine (Moderna mRNA-1273 Beta/Omicron BA.1 bivalent vaccine) in adults. A 2-dose boosting regimen with a variant vaccine did not increase the magnitude or the durability of the serological responses compared to a single variant vaccine boost.

Immunogenicity of the BA.1 and BA.4/BA.5 Severe Acute Respiratory Syndrome Coronavirus 2 Bivalent Boosts: Preliminary Results From the COVAIL Randomized Clinical Trial.

In a randomized clinical trial, we compare early neutralizing antibody responses after boosting with bivalent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) messenger RNA (mRNA) vaccines based on either BA.1 or BA.4/BA.5 Omicron spike protein combined with wild-type spike. Responses against SARS-CoV-2 variants exhibited the greatest reduction in titers against currently circulating Omicron subvariants for both bivalent vaccines.

US National Institutes of Health Prioritization of SARS-CoV-2 Variants.

Since late 2020, SARS-CoV-2 variants have regularly emerged with competitive and phenotypic differences from previously circulating strains, sometimes with the potential to escape from immunity produced by prior exposure and infection. The Early Detection group is one of the constituent groups of the US National Institutes of Health National Institute of Allergy and Infectious Diseases SARS-CoV-2 Assessment of Viral Evolution program. The group uses bioinformatic methods to monitor the emergence, spread, and potential phenotypic properties of emerging and circulating strains to identify the most relevant variants for experimental groups within the program to phenotypically characterize. Since April 2021, the group has prioritized variants monthly. Prioritization successes include rapidly identifying most major variants of SARS-CoV-2 and providing experimental groups within the National Institutes of Health program easy access to regularly updated information on the recent evolution and epidemiology of SARS-CoV-2 that can be used to guide phenotypic investigations.

Uptake of Phytoplankton-Derived Carbon and Cobalamins by Novel Acidobacteria Genera in Microcystis Blooms Inferred from Metagenomic and Metatranscriptomic Evidence.

Interactions between bacteria and phytoplankton can influence primary production, community composition, and algal bloom development. However, these interactions are poorly described for many consortia, particularly for freshwater bloom-forming cyanobacteria. Here, we assessed the gene content and expression of two uncultivated Acidobacteria from Lake Erie Microcystis blooms. These organisms were targeted because they were previously identified as important catalase producers in Microcystis blooms, suggesting that they protect Microcystis from H2O2. Metatranscriptomics revealed that both Acidobacteria transcribed genes for uptake of organic compounds that are known cyanobacterial products and exudates, including lactate, glycolate, amino acids, peptides, and cobalamins. Expressed genes for amino acid metabolism and peptide transport and degradation suggest that use of amino acids and peptides by Acidobacteria may regenerate nitrogen for cyanobacteria and other organisms. The Acidobacteria genomes lacked genes for biosynthesis of cobalamins but expressed genes for its transport and remodeling. This indicates that the Acidobacteria obtained cobalamins externally, potentially from Microcystis, which has a complete gene repertoire for pseudocobalamin biosynthesis; expressed them in field samples; and produced pseudocobalamin in axenic culture. Both Acidobacteria were detected in Microcystis blooms worldwide. Together, the data support the hypotheses that uncultured and previously unidentified Acidobacteria taxa exchange metabolites with phytoplankton during harmful cyanobacterial blooms and influence nitrogen available to phytoplankton. Thus, novel Acidobacteria may play a role in cyanobacterial physiology and bloom development. IMPORTANCE Interactions between heterotrophic bacteria and phytoplankton influence competition and successions between phytoplankton taxa, thereby influencing ecosystem-wide processes such as carbon cycling and algal bloom development. The cyanobacterium Microcystis forms harmful blooms in freshwaters worldwide and grows in buoyant colonies that harbor other bacteria in their phycospheres. Bacteria in the phycosphere and in the surrounding community likely influence Microcystis physiology and ecology and thus the development of freshwater harmful cyanobacterial blooms. However, the impacts and mechanisms of interaction between bacteria and Microcystis are not fully understood. This study explores the mechanisms of interaction between Microcystis and uncultured members of its phycosphere in situ with population genome resolution to investigate the cooccurrence of Microcystis and freshwater Acidobacteria in blooms worldwide.

Heterotrophic Bacteria Dominate Catalase Expression during Microcystis Blooms.

In the oligotrophic oceans, key autotrophs depend on "helper" bacteria to reduce oxidative stress from hydrogen peroxide (H2O2) in the extracellular environment. H2O2 is also a ubiquitous stressor in freshwaters, but the effects of H2O2 on autotrophs and their interactions with bacteria are less well understood in freshwaters. Naturally occurring H2O2 in freshwater systems is proposed to impact the proportion of microcystin-producing (toxic) and non-microcystin-producing (nontoxic) Microcystis in blooms, which influences toxin concentrations and human health impacts. However, how different strains of Microcystis respond to naturally occurring H2O2 concentrations and the microbes responsible for H2O2 decomposition in freshwater cyanobacterial blooms are unknown. To address these knowledge gaps, we used metagenomics and metatranscriptomics to track the presence and expression of genes for H2O2 decomposition by microbes during a cyanobacterial bloom in western Lake Erie in the summer of 2014. katG encodes the key enzyme for decomposing extracellular H2O2 but was absent in most Microcystis cells. katG transcript relative abundance was dominated by heterotrophic bacteria. In axenic Microcystis cultures, an H2O2 scavenger (pyruvate) significantly improved growth rates of one toxic strain while other toxic and nontoxic strains were unaffected. These results indicate that heterotrophic bacteria play a key role in H2O2 decomposition in Microcystis blooms and suggest that their activity may affect the fitness of some Microcystis strains and thus the strain composition of Microcystis blooms but not along a toxic versus nontoxic dichotomy. IMPORTANCE Cyanobacterial harmful algal blooms (CHABs) threaten freshwater ecosystems globally through the production of toxins. Toxin production by cyanobacterial species and strains during CHABs varies widely over time and space, but the ecological drivers of the succession of toxin-producing species remain unclear. Hydrogen peroxide (H2O2) is ubiquitous in natural waters, inhibits microbial growth, and may determine the relative proportions of Microcystis strains during blooms. However, the mechanisms and organismal interactions involved in H2O2 decomposition are unexplored in CHABs. This study shows that some strains of bloom-forming freshwater cyanobacteria benefit from detoxification of H2O2 by associated heterotrophic bacteria, which may impact bloom development.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed.

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

Individual Microcystis colonies harbour distinct bacterial communities that differ by Microcystis oligotype and with time.

Interactions between bacteria and phytoplankton in the phycosphere have impacts at the scale of whole ecosystems, including the development of harmful algal blooms. The cyanobacterium Microcystis causes toxic blooms that threaten freshwater ecosystems and human health globally. Microcystis grows in colonies that harbour dense assemblages of other bacteria, yet the taxonomic composition of these phycosphere communities and the nature of their interactions with Microcystis are not well characterized. To identify the taxa and compositional variance within Microcystis phycosphere communities, we performed 16S rRNA V4 region amplicon sequencing on individual Microcystis colonies collected biweekly via high-throughput droplet encapsulation during a western Lake Erie cyanobacterial bloom. The Microcystis phycosphere communities were distinct from microbial communities in whole water and bulk phytoplankton seston in western Lake Erie but lacked 'core' taxa found across all colonies. However, dissimilarity in phycosphere community composition correlated with sampling date and the Microcystis 16S rRNA oligotype. Several taxa in the phycosphere were specific to and conserved with Microcystis of a single oligotype or sampling date. Together, this suggests that physiological differences between Microcystis strains, temporal changes in strain phenotypes, and the composition of seeding communities may impact community composition of the Microcystis phycosphere.

The genetic and ecophysiological diversity of Microcystis.

Microcystis is a cyanobacterium that forms toxic blooms in freshwater ecosystems around the world. Biological variation among taxa within the genus is apparent through genetic and phenotypic differences between strains and via the spatial and temporal distribution of strains in the environment, and this fine-scale diversity exerts strong influence over bloom toxicity. Yet we do not know how varying traits of Microcystis strains govern their environmental distribution, the tradeoffs and links between these traits, or how they are encoded at the genomic level. Here we synthesize current knowledge on the importance of diversity within Microcystis and on the genes and traits that likely underpin ecological differentiation of taxa. We briefly review spatial and environmental patterns of Microcystis diversity in the field and genetic evidence for cohesive groups within Microcystis. We then compile data on strain-level diversity regarding growth responses to environmental conditions and explore evidence for variation of community interactions across Microcystis strains. Potential links and tradeoffs between traits are identified and discussed. The resulting picture, while incomplete, highlights key knowledge gaps that need to be filled to enable new models for predicting strain-level dynamics, which influence the development, toxicity and cosmopolitan nature of Microcystis blooms.

Global circulation patterns of seasonal influenza viruses vary with antigenic drift.

Understanding the spatiotemporal patterns of emergence and circulation of new human seasonal influenza virus variants is a key scientific and public health challenge. The global circulation patterns of influenza A/H3N2 viruses are well characterized, but the patterns of A/H1N1 and B viruses have remained largely unexplored. Here we show that the global circulation patterns of A/H1N1 (up to 2009), B/Victoria, and B/Yamagata viruses differ substantially from those of A/H3N2 viruses, on the basis of analyses of 9,604 haemagglutinin sequences of human seasonal influenza viruses from 2000 to 2012. Whereas genetic variants of A/H3N2 viruses did not persist locally between epidemics and were reseeded from East and Southeast Asia, genetic variants of A/H1N1 and B viruses persisted across several seasons and exhibited complex global dynamics with East and Southeast Asia playing a limited role in disseminating new variants. The less frequent global movement of influenza A/H1N1 and B viruses coincided with slower rates of antigenic evolution, lower ages of infection, and smaller, less frequent epidemics compared to A/H3N2 viruses. Detailed epidemic models support differences in age of infection, combined with the less frequent travel of children, as probable drivers of the differences in the patterns of global circulation, suggesting a complex interaction between virus evolution, epidemiology, and human behaviour.

Ancient human parvovirus B19 in Eurasia reveals its long-term association with humans.

Human parvovirus B19 (B19V) is a ubiquitous human pathogen associated with a number of conditions, such as fifth disease in children and arthritis and arthralgias in adults. B19V is thought to evolve exceptionally rapidly among DNA viruses, with substitution rates previously estimated to be closer to those typical of RNA viruses. On the basis of genetic sequences up to ∼70 years of age, the most recent common ancestor of all B19V has been dated to the early 1800s, and it has been suggested that genotype 1, the most common B19V genotype, only started circulating in the 1960s. Here we present 10 genomes (63.9-99.7% genome coverage) of B19V from dental and skeletal remains of individuals who lived in Eurasia and Greenland from ∼0.5 to ∼6.9 thousand years ago (kya). In a phylogenetic analysis, five of the ancient B19V sequences fall within or basal to the modern genotype 1, and five fall basal to genotype 2, showing a long-term association of B19V with humans. The most recent common ancestor of all B19V is placed ∼12.6 kya, and we find a substitution rate that is an order of magnitude lower than inferred previously. Further, we are able to date the recombination event between genotypes 1 and 3 that formed genotype 2 to ∼5.0-6.8 kya. This study emphasizes the importance of ancient viral sequences for our understanding of virus evolution and phylogenetics.

Epistatic interactions can moderate the antigenic effect of substitutions in haemagglutinin of influenza H3N2 virus.

We previously showed that single amino acid substitutions at seven positions in haemagglutinin determined major antigenic change of influenza H3N2 virus. Here, the impact of two such substitutions was tested in 11 representative H3 haemagglutinins to investigate context-dependence effects. The antigenic effect of substitutions introduced at haemagglutinin position 145 was fully independent of the amino acid context of the representative haemagglutinins. Antigenic change caused by substitutions introduced at haemagglutinin position 155 was variable and context-dependent. Our results suggest that epistatic interactions with contextual amino acids in the haemagglutinin can moderate the magnitude of antigenic change.

Limited airborne transmission of H7N9 influenza A virus between ferrets.

Wild waterfowl form the main reservoir of influenza A viruses, from which transmission occurs directly or indirectly to various secondary hosts, including humans. Direct avian-to-human transmission has been observed for viruses of subtypes A(H5N1), A(H7N2), A(H7N3), A(H7N7), A(H9N2) and A(H10N7) upon human exposure to poultry, but a lack of sustained human-to-human transmission has prevented these viruses from causing new pandemics. Recently, avian A(H7N9) viruses were transmitted to humans, causing severe respiratory disease and deaths in China. Because transmission via respiratory droplets and aerosols (hereafter referred to as airborne transmission) is the main route for efficient transmission between humans, it is important to gain an insight into airborne transmission of the A(H7N9) virus. Here we show that although the A/Anhui/1/2013 A(H7N9) virus harbours determinants associated with human adaptation and transmissibility between mammals, its airborne transmissibility in ferrets is limited, and it is intermediate between that of typical human and avian influenza viruses. Multiple A(H7N9) virus genetic variants were transmitted. Upon ferret passage, variants with higher avian receptor binding, higher pH of fusion, and lower thermostability were selected, potentially resulting in reduced transmissibility. This A(H7N9) virus outbreak highlights the need for increased understanding of the determinants of efficient airborne transmission of avian influenza viruses between mammals.

Genetic and antigenic characterisation of influenza A(H3N2) viruses isolated in Yokohama during the 2016/17 and 2017/18 influenza seasons.

BACKGROUND: Influenza A(H3N2) virus rapidly evolves to evade human immune responses, resulting in changes in the antigenicity of haemagglutinin (HA). Therefore, continuous genetic and antigenic analyses of A(H3N2) virus are necessary to detect antigenic mutants as quickly as possible. AIM: We attempted to phylogenetically and antigenically capture the epidemic trend of A(H3N2) virus infection in Yokohama, Japan during the 2016/17 and 2017/18 influenza seasons. METHODS: We determined the HA sequences of A(H3N2) viruses detected in Yokohama, Japan during the 2016/17 and 2017/18 influenza seasons to identify amino acid substitutions and the loss or gain of potential N-glycosylation sites in HA, both of which potentially affect the antigenicity of HA. We also examined the antigenicity of isolates using ferret antisera obtained from experimentally infected ferrets. RESULTS: Influenza A(H3N2) viruses belonging to six clades (clades 3C.2A1, 3C.2A1a, 3C.2A1b, 3C.2A2, 3C.2A3 and 3C.2A4) were detected during the 2016/17 influenza season, whereas viruses belonging to two clades (clades 3C.2A1b and 3C.2A2) dominated during the 2017/18 influenza season. The isolates in clades 3C.2A1a and 3C.2A3 lost one N-linked glycosylation site in HA relative to other clades. Antigenic analysis revealed antigenic differences among clades, especially clade 3C.2A2 and 3C.2A4 viruses, which showed distinct antigenic differences from each other and from other clades in the antigenic map. CONCLUSION: Multiple clades, some of which differed antigenically from others, co-circulated in Yokohama, Japan during the 2016/17 and 2017/18 influenza seasons.