key: cord-0973611-cci7z0ab
authors: Guthmiller, Jenna J.; Han, Julianna; Utset, Henry A.; Li, Lei; Lan, Linda Yu-Ling; Henry, Carole; Stamper, Christopher T.; McMahon, Meagan; O’Dell, George; Fernández-Quintero, Monica L.; Freyn, Alec W.; Amanat, Fatima; Stovicek, Olivia; Gentles, Lauren; Richey, Sara T.; de la Peña, Alba Torrents; Rosado, Victoria; Dugan, Haley L.; Zheng, Nai-Ying; Tepora, Micah E.; Bitar, Dalia J.; Changrob, Siriruk; Strohmeier, Shirin; Huang, Min; García-Sastre, Adolfo; Liedl, Klaus R.; Bloom, Jesse D.; Nachbagauer, Raffael; Palese, Peter; Krammer, Florian; Coughlan, Lynda; Ward, Andrew B.; Wilson, Patrick C.
title: Broadly neutralizing antibodies target a haemagglutinin anchor epitope
date: 2021-12-23
journal: Nature
DOI: 10.1038/s41586-021-04356-8
sha: c4ab98c6e0e4d582c88f9747752f66a4fe6c4099
doc_id: 973611
cord_uid: cci7z0ab

Broadly neutralizing antibodies that target epitopes of haemagglutinin on the influenza virus have the potential to provide near universal protection against influenza virus infection(1). However, viral mutants that escape broadly neutralizing antibodies have been reported(2,3). The identification of broadly neutralizing antibody classes that can neutralize viral escape mutants is critical for universal influenza virus vaccine design. Here we report a distinct class of broadly neutralizing antibodies that target a discrete membrane-proximal anchor epitope of the haemagglutinin stalk domain. Anchor epitope-targeting antibodies are broadly neutralizing across H1 viruses and can cross-react with H2 and H5 viruses that are a pandemic threat. Antibodies that target this anchor epitope utilize a highly restricted repertoire, which encodes two public binding motifs that make extensive contacts with conserved residues in the fusion peptide. Moreover, anchor epitope-targeting B cells are common in the human memory B cell repertoire and were recalled in humans by an oil-in-water adjuvanted chimeric haemagglutinin vaccine(4,5), which is a potential universal influenza virus vaccine. To maximize protection against seasonal and pandemic influenza viruses, vaccines should aim to boost this previously untapped source of broadly neutralizing antibodies that are widespread in the human memory B cell pool.

Broadly neutralizing antibodies that target epitopes of haemagglutinin on the influenza virus have the potential to provide near universal protection against influenza virus infection 1 . However, viral mutants that escape broadly neutralizing antibodies have been reported 2, 3 . The identification of broadly neutralizing antibody classes that can neutralize viral escape mutants is critical for universal influenza virus vaccine design. Here we report a distinct class of broadly neutralizing antibodies that target a discrete membrane-proximal anchor epitope of the haemagglutinin stalk domain. Anchor epitope-targeting antibodies are broadly neutralizing across H1 viruses and can cross-react with H2 and H5 viruses that are a pandemic threat. Antibodies that target this anchor epitope utilize a highly restricted repertoire, which encodes two public binding motifs that make extensive contacts with conserved residues in the fusion peptide. Moreover, anchor epitope-targeting B cells are common in the human memory B cell repertoire and were recalled in humans by an oil-in-water adjuvanted chimeric haemagglutinin vaccine 4, 5 , which is a potential universal influenza virus vaccine. To maximize protection against seasonal and pandemic influenza viruses, vaccines should aim to boost this previously untapped source of broadly neutralizing antibodies that are widespread in the human memory B cell pool.

Antibodies to the major surface glycoprotein haemagglutinin (HA) are critical for providing protection against influenza virus infection 6, 7 . However, most HA-binding antibodies target variable epitopes of the HA head domain, which provide limited protection against antigenically similar influenza virus strains 3 . Vaccine formulations that preferentially induce antibodies to conserved epitopes of the HA head and stalk domains could provide broad and potent protection against a wide array of influenza viruses. Several leading universal influenza virus candidates are designed to induce antibodies specifically to the stalk domain, but the spectrum of distinct epitopes on the stalk targeted by the human B cell repertoire remains to be determined. By analysing the specificities of B cells targeting the H1 stalk through the generation of monoclonal antibodies (mAbs), our study reveals a new class of broadly neutralizing antibodies (bnAbs) to an underappreciated epitope where HA anchors itself into the viral membrane. Next-generation vaccine platforms should take advantage of this finding to elicit antibodies to the conserved anchor epitope.

To investigate the specificities of HA-specific antibodies, we generated 358 mAbs from plasmablasts and HA + memory B cells (MBCs) isolated from volunteers who were vaccinated against or naturally infected with seasonal influenza viruses or were participants in a phase I clinical trial of a chimeric HA (cHA) vaccine 4, 5 . Of all mAbs tested, nearly half targeted the HA stalk domain, 21% of which targeted the well-characterized central stalk (CS) epitope (Extended Data Fig. 1a) . Notably, stalk-binding mAbs were disproportionally isolated from the infected, 2009 pandemic H1N1 (pH1N1) monovalent inactivated influenza vaccine, and the cHA vaccine cohorts (Extended Data Fig. 1b) , as these exposure routes have previously been shown to induce antibody responses to the HA stalk 5, 8, 9 . To investigate which epitopes the remaining 79% of stalk-binding mAbs were targeting, we performed negative-stain electron microscopy with three non-CR9114 (ref. 10 ) competing stalk domain-binding mAbs: 047-09 4F04, 241 IgA 2F04 and 222-1C06. All three mAbs bound an epitope near the anchor of the HA stalk and were oriented at an upward angle towards the epitope (Fig. 1a , b, Extended Data Fig. 1c) , suggesting that this epitope may be partially obstructed by the lipid membrane and may only be exposed for antibody binding as the HA trimers flex on the membrane. FISW84, a recently identified stalk-binding mAb 11 , targets an epitope that overlaps with the three identified anchor-binding mAbs (Extended Data Fig. 1d) , suggesting that the anchor epitope is a common stalk epitope. Moreover, a proximal epitope was previously identified on group 2 viruses that is targeted by mAbs CR8020 (ref. 12 ) and CR8043 (ref. 13 ). Despite some overlap, the epitope targeted by the group 2 mAbs was considerably farther up and to the right on the HA stalk relative to the anchor epitope, and mAbs CR8020 and CR8043 targeted the stalk from above (Extended Data Fig. 1e, f) , at an angle and positioning more akin to antibodies targeting the CS epitope. mAbs binding the CS epitope (CR9114 (ref. 10 ) and FI6v3 (ref. 14 )) did not have overlapping footprints or compete for HA binding with the anchor-binding 047-09 4F04 mAb in an HA competition assay (Fig. 1c, d) . In total, we identified 50 distinct mAbs that competed for binding to the anchor epitope from a total of 21 individuals (Fig. 1d , Extended Data Tables 1, 2). Of these, 34 anchor-binding mAbs from 15 donors were isolated from the cHA vaccine trial (Extended Data Tables 1, 2).

Anchor-binding mAbs were broadly reactive and broadly neutralizing to pre-pandemic and post-pandemic H1N1 viruses and a swine-origin H1N2 virus (Fig. 1e , Extended Data Fig. 1g, h) . Notably, anchor-binding mAbs had similar neutralizing potency to pH1N1 relative to mAbs to the CS epitope (Fig. 1e , Extended Data Fig. 1i , j). Many stalk-targeting antibodies mediate protection via Fc-mediated functions, including antibody-dependent cellular cytotoxicity 15, 16 . Anchor epitope-binding mAbs largely did not possess antibody-dependent cellular cytotoxicity activity (15 out of 18; Extended Data Fig. 1k-l) , potentially due to the upward angle of approach of anchor-binding antibodies, which may position the Fc distantly from effector cells. Despite pan-H1 binding, anchor-targeting mAbs rarely cross-reacted with other HA subtypes tested, including H3, a group 2 subtype, other group 1 subtypes, including H2 and H5, and influenza B viruses (Fig. 1f, Extended Data Fig. 2a ). Despite this, the 222-1C06 mAb cross-reacted with H2 and H5 HA (Fig. 1g , Extended Data Fig. 2b ) and several anchor-binding antibodies could neutralize an H2N2 virus (Extended Data Fig. 2c ), suggesting that antibodies targeting the anchor epitope can cross-neutralize other group 1 influenza A viruses. We recently demonstrated that bnAbs to the HA stalk are often polyreactive 17 , which may limit the activation of B cells expressing bnAbs. Relative to mAbs that target the CS epitope, we identified that anchor-binding mAbs were proportionally less likely to be polyreactive and those that were polyreactive had weaker relative affinity for lipopolysaccharide (Extended Data Fig. 2d, e) . These data suggest that, although polyreactivity is selected for within the anti-anchor epitope B cell pool, it is not to the same level as B cells to the CS epitope.

H1N1 viruses have acquired several mutations within the HA stalk domain, probably due to antibody selective pressures or to increase stability. To understand whether these mutations affect antibody binding to the anchor epitope, we screened mAbs against naturally occurring and experimentally identified viral escape mutants of mAbs binding to the CS epitope (Extended Data Fig. 2f , g, Extended Data Table 3 ). Anchor epitope-binding mAbs were mostly unaffected by these mutants, whereas most of the CS-binding mAbs showed reduced binding to at least one mutant (Extended Data Fig. 2g) . Notably, most mAbs had reduced binding to A373V of HA2, which has recently been shown to preferentially grow in the presence of mAbs to the CS epitope 2 . While A373 is distant from the anchor epitope, the A373V mutation has been shown to affect the conformation of the HA stalk 2 , explaining the broad reduction of HA binding by antibodies targeting either stalk epitope. Anchor-binding mAbs only demonstrated a 10-30% reduction in binding (Extended Data Fig. 2g) , indicating that they are still likely to neutralize viruses carrying the A373V mutation. 

To test whether mAbs targeting the anchor epitope were protective in vivo, we administered a cocktail of five mAbs targeting the anchor epitope or the CS epitope prophylactically and therapeutically (48 h after infection) to mimic a polyclonal response against these epitopes and infected mice with a lethal dose of a mouse-adapted pH1N1 virus (A/Netherlands/602/2009; Supplementary Table 1 ). Mice that received a prophylactic or therapeutic anti-anchor cocktail at 5 mg/kg experienced 100% protection from weight loss and lethal infection (Extended Data Fig. 3a , b). No differences in lung viral titres were detected in mice that received the anti-anchor cocktail prophylactically relative to mice administered the negative control mAb (Extended Data Fig. 3c ). Anti-stalk antibodies do not provide sterilizing immunity but are known to neutralize subsequent rounds and limit disease progression 18 . As a result, lung viral titres do not necessarily reflect protection from morbidity and mortality. Finally, the anti-anchor cocktail given prophylactically provided 70% protection against lethal A/Fort Monmouth/1/1947 infection (Extended Data Fig. 3d) , a virus that circulated before the birth of any of the donors in our study (Extended Data Table 1 ). Since anchor-binding mAbs do not engage in antibody-dependent cellular cytotoxicity for the most part, antibodies that target the anchor epitope probably provide protection in vivo through direct neutralization of the virus. Together, these data indicate that antibodies to the anchor epitope are pan-H1 neutralizing and protective in vivo.

All anchor epitope-binding mAbs utilized one of four VH3 genes: VH3-23, VH3-30 or VH3-30-3, and VH3-48, in contrast to mAbs targeting the CS epitope, which commonly use VH1-69 (Fig. 2a , Extended Data Fig. 4a, b) . Anchor epitope-binding mAbs also utilized a highly restricted light chain repertoire, with all mAbs utilizing a combination of VK3-11 or VK3-15 combined with JK4 or JK5 (Fig. 2b , Extended Data Fig. 4c, d) . Furthermore, all but one light chain of the anchor-targeting mAbs were clonally related (Extended Data Table 2 ), indicating that the light chains were very similar across mAbs and study participants. We identified four distinct clonal expansions, with one public clone found across two donors (Extended Data Fig. 4e , f, Extended Data Table 2 ). Anchor epitope-targeting and CS-targeting mAbs exhibited similar levels of somatic hypermutations (Extended Data Fig. 4g ). The κ-chain complementarity-determining region 3 (K-CDR3) length of anchor epitope-binding mAbs was highly restricted, with all K-CDR3s being ten amino acids long (Extended Data Fig. 4h) .

To investigate the binding motif of anchor-targeting mAbs, we generated a high-resolution (3.38 Å) cryo-electron microscopy (cryo-EM) structure of 222-1C06 bound to A/California/7/2009 HA (Fig. 2c , Extended Data Fig. 4i-k) . The broad paratope of 222-1C06 made extensive contacts across the HA fusion peptide 19 (Fig. 2d, e) , which mediates viral membrane fusion with the host membrane. Binding of the fusion peptide was largely mediated by an NWP motif within the K-CDR3, a Y58 directly following the H-CDR2, and a W99 in the H-CDR3, with these HA-binding motifs acting independently and in combination via an aromatic pocket (Fig. 2f , g, Supplementary Table 2 ). Moreover, both the K-CDR3 NWP and the H-CDR2 Y58 were found in all the anchor-binding mAbs and were germline encoded (Fig. 2h, i) , which could have led to the selection of B cells utilizing these variable genes. Notably, FISW84 utilizes VH3-23 and VK3-15, and comparison of the paratopes showed that the NWP and Y58 motifs of FISW84 similarly form an aromatic pocket but orient towards the fusion peptide slightly differently than Fig. 4l ). Molecular dynamics simulations also showed that VH3-23-utilizing and VH3-30-utilizing mAbs from our study and FISW84 targeted the HA fusion peptide similarly to the cryo-EM structure of 222-1C06 via the aromatic pocket created by the K-CDR3 NWP and H-CDR2 Y58 motifs, albeit at different orientations (Extended Data Fig. 4m ). Crucially, molecular dynamic and cryo-EM analyses revealed numerous intra-Fab interactions of hydrophobic and aromatic amino acids, including p-stacking of K-CDR3-P95 with K-CDR3-W94 and H-CDR2-Y58 that rigidified the paratope (Extended Data Fig. 4n , o, Supplementary Table 3 ). Broad analysis of human, swine and avian H1 viruses revealed that the side-chain contacts of 222-1C06 were highly conserved (94-100% conserved; Fig. 2j ). In addition, the side-chain contacts were 100% conserved across 100 years of H1N1 virus evolution in humans (Extended Data Fig. 4p ). Deep mutational scanning of potential H1 viruses at these contact residues indicated substantial permissibility (Extended Data Fig. 4q ), although these mutations appear to not have been selected for in nature (Fig. 2j) . The five side-chain contacts of this broad epitope were also highly conserved across group 1 viruses, with the W343 contact being 100% conserved across all group 1 viruses (Fig. 2k) . Together, these data reveal that B cells targeting the anchor epitope utilized a highly restricted V(D)J gene repertoire, and the specific features within this repertoire made critical and extensive contacts with the conserved anchor epitope.

Owing to the restricted repertoire features of anchor-targeting mAbs, we next determined the relative proportion of B cells with these features by interrogating single-cell repertoire sequencing of HA-specific B cells isolated from 20 human participants following cHA vaccination 4, 5 . The cHA vaccine platform is intended to specifically induce antibodies to the stalk domain by retaining the stalk domain of H1 and replacing the head domain of HA with that of an avian subtype, either H8 (prime) or H5 (boost) for this trial 4, 5 . To investigate the proportion of B cells with anchor epitope-binding repertoire features, we selected B cells that used VH3-23/VH3-30/VH3-30-3/VH3-48, VK3-11/VK3-15, JK4/JK5, a K-CDR3 ten amino acids in length, and possessed an NWP motif within the K-CDR3. We also segregated B cells expressing VH1-69 and a κ-chain, which are commonly used by B cells that target the CS epitope. We identified that B cells with features of antibodies binding the anchor epitope were abundant within the human B cell repertoire, with 6% of all B cells identified fitting within this defined repertoire (Fig. 3a) . Moreover, all but one participant had at least one B cell with anchor-binding repertoire features (Fig. 3b ). In addition, 32 out of 33 mAbs generated from the selected anchor-targeting B cell list competed for HA binding with 047-09 4F04 (Extended Data Fig. 5a ), indicating that this population was greatly enriched for anchor-targeting B cells. Moreover, the anchor-binding B cells were highly mutated and were largely class-switched to IgG1 (Extended Data Fig. 5b , c), suggesting that these B cells were MBCs that had undergone affinity maturation and class-switch recombination. Together, these data indicate that the anchor epitope is a common target of the human MBC repertoire following cHA vaccination.

To confirm that anchor epitope-targeting mAbs were representative of the serum antibody response, we performed electron microscopy polyclonal epitope mapping (EMPEM) 20 with serum antibodies from participants 236 and 241 from the 2014 quadrivalent inactivated influenza vaccine cohort (Extended Data Table 1 ). Both participants had detectable antibodies targeting the anchor epitope at days 7 and 14 post-vaccination ( Fig. 3c , Extended Data Fig. 5d , e). Comparison of anchor epitope-binding polyclonal antibodies identified in participant 241 revealed that the 241 IgA 2F04 mAb strongly overlapped with the 241 polyclonal antibody (Fig. 3d) , suggesting that the polyclonal antibody derived from this clonal expansion. Together, these data indicate that humoral immunity against the anchor epitope is common within the MBC pool and polyclonal serum antibody response after vaccination.

To investigate whether participants enrolled in a phase I clinical trial of the cHA vaccine ( Fig. 4a ) mounted an antibody response to the anchor and the CS epitope, we adapted the competition ELISA to detect serum antibody responses that could compete for binding with 047-09 4F04 and CR9114, respectively. Three different vaccine formulations were used in this trial, with participants being primed with either a cH8/1 inactivated influenza vaccine (IIV) with an oil-in-water adjuvant (AS03) or a cH8/1 live-attenuated influenza vaccine followed by a boost with a cH5/1 IIV with or without AS03 (Fig. 4a) . Only participants who received the IIV + AS03 on either the prime or boost were capable of seroconverting against both the anchor and the CS epitopes (Fig. 4b , c, Extended Data Fig. 6a, b) . Participants who received the cH8/1 IIV + AS03 prime did not further increase serum antibodies to either stalk epitope after the cH5/1 + AS03 boost (Fig. 4c) , suggesting that these B cells were refractory to continued activation. Serum titres against the anchor and CS epitopes (Extended Data Fig. 6a 

neutralizing titres against a cH6/1N8 virus and an avian-swine H1N1 virus 5 , suggesting that the anchor-targeting and CS-targeting serum antibodies were responsible for neutralization. All but one participant in the IIV + AS03 groups seroconverted against the anchor epitope ( Fig. 4d ) and had higher titres against the anchor epitope than participants who received the 2009 monovalent inactivated influenza vaccine or seasonal influenza virus vaccines (Fig. 4e) . However, the precise role of the cHA immunogen versus the AS03 adjuvant in inducing anti-stalk antibody responses could not be resolved in our study. Notably, those individuals who received the IIV alone had weak plasmablast responses relative to those individuals who received the IIV + AS03 (ref. 4 ), suggesting that the oil-in-water adjuvant, not the cHA immunogen, was essential for robust activation of B cells and anti-anchor antibody responses. Moreover, considerable pre-existing antibodies may hinder recall of B cells to the stalk. Individuals first exposed to the 2009 pH1N1, a virus for which individuals had low pre-existing humoral immunity, had proportionally more isolated mAbs to the stalk and were more likely to have an anti-anchor mAb than individuals who had repeatedly been exposed to the pH1N1 virus in subsequent influenza seasons (Extended Data Fig. 6c, d) , therefore suggesting that pre-existing antibodies may limit antibody responses to the HA stalk. Despite robust recall of antibodies to the anchor and CS epitopes by the adjuvanted vaccines, titres decreased 1 year after vaccination (day 420; Extended Data Fig. 6e , f). Together, these data indicate that the adjuvanted inactivated cHA vaccine can robustly induce antibodies to multiple stalk epitopes, including the anchor.

Headless HA antigens, including mini-HA 21 , are attractive universal influenza virus vaccine antigens, as these antigens lack the immunodominant epitopes of the HA head 21,22 . However, only 1 out of 50 anchor-binding mAbs bound the mini-HA antigen 21 , whereas all anchor epitope-binding mAbs bound cH6/1 (Extended Data Fig. 6g) . Compared to full-length HA, the membrane-proximal region of the mini-HA splays by approximately 14.5 Å (ref. 21 ), which may disrupt the antigenicity of the anchor epitope. To understand whether anchor epitope-targeting antibodies could bind to the mini-HA in a more native setting, we generated a membrane-bound mini-HA and observed that mAbs targeting the anchor and CS epitopes readily bind to both the full-length membrane-bound A/California/7/2009 HA and the membrane-bound mini-HA (Extended Data Fig. 6h) , indicating that the mini-HA is antigenic when natively presented. Furthermore, we demonstrated that anchor epitope-targeting antibodies bound HA with a fibritin but not a GCN4 trimerization domain (Extended Data Fig. 6i ), highlighting selection of the trimerization domain as an important factor for vaccine design. Together, these data demonstrate that native-like HA antigens, such as the cHA vaccine, can recall MBCs that target the anchor epitope.

In this study, we identified a public class of bnAbs that target an epitope at the anchor of the HA stalk domain near the membrane. The anchor-targeting mAbs were public clonotypes across participants, with all antibodies possessing two conserved, germline-encoded binding motifs: a Y58 directly following the H-CDR2 and an NWP motif within the K-CDR3. The neutralizing activity of anchor epitope-targeting mAbs to pre-pH1N1 and post-pH1N1 viruses and a swine-origin H1-expressing virus indicates that the anchor epitope is an important target for pan-subtype neutralizing antibodies. As two of the last four influenza virus pandemics were caused by H1N1 viruses and a recent report has shown that antigenically novel H1-expressing viruses commonly spill over from swine into humans 23 , it is critically important to generate pan-H1 vaccines to prevent the next influenza pandemic. Moreover, the ability of anti-anchor antibodies to neutralize an H2 virus and the general conservation of contact residues suggests that anchor-targeting antibodies have the potential to acquire cross-neutralization against group 1 viruses. Our study highlights an additional broadly protective epitope of the HA stalk and provides guidance on how vaccines can be designed to drive broadly protective antibodies to multiple distinct epitopes, which can work cooperatively to provide optimal protection while avoiding the generation of antibody escape mutants. Notably, studies in the HIV-1 field have shown that bnAb monotherapy can lead to the development of antibody-resistant viral variants 24-26 , whereas combination bnAb therapy does not 27 . In addition, immune focusing towards a single epitope can lead to the generation and selection of viral escape mutants at these highly conserved epitopes 2,3,28,29 . Therefore, it is critical that future universal influenza virus vaccines elicit antibodies to multiple conserved epitopes to prevent the generation of bnAb escape viruses.

The angle of approach of anchor-binding mAbs and the proximity of the epitope to the viral membrane suggest that this epitope is typically obstructed, limiting antibody recognition and B cell activation. However, membrane-bound HA typically flexes between 0° and 25° and up to 52° on its threefold axis 11 , suggesting antibodies and B cells can easily access the anchor epitope during this flexing process. Moreover, H1 viruses possess a highly conserved glycosylation site on the HA stalk that lies above the anchor epitope 11, 30 , which may obstruct antibodies from targeting this epitope from above. Similarly, a neutralizing mAb to the Middle East respiratory syndrome coronavirus targets an epitope on S2 from an upward angle to avoid glycans 31 . Therefore, the upward angle of approach may be a common feature of antibodies that recognize epitopes below glycans.

Our study shows that humans have pre-existing immunity against the anchor epitope and influenza virus vaccination can recall MBCs to secrete antibodies to this epitope. However, vaccine HA antigens must have a native-like conformation near the transmembrane domain, as trimer splaying potentially due to the GCN4 trimerization domain ablates antibody binding to the anchor epitope. Moreover, our study reveals that the cHA vaccine strategy recalled MBCs to the anchor and CS epitopes, as these MBCs do not need to compete against MBCs that target the immunodominant variable HA head epitopes 3, 9, 17, 32 . Similarly, the mini-HA/headless HA vaccine strategy has the potential to also recall MBCs to multiple epitopes of the HA stalk domain, if displayed natively 22,33 . The addition of an oil-in-water AS03 adjuvant to the cHA was critical for recalling MBCs to the anchor epitope. Oil-in-water adjuvants largely function to emulsify antigen, which may prevent sequestration of antigen by circulating antibodies, increase delivery of free antigen to lymph nodes, and help to stimulate innate immune receptors 34, 35 . Notably, an oil-in-water adjuvanted inactivated H5N1 vaccine induced higher neutralizing antibody titres, antibodies to more HA epitopes, and higher avidity antibodies 36 . Therefore, the inclusion of oil-in-water adjuvants may have an important role in generating bnAbs and may improve vaccine effectiveness of seasonal influenza vaccines. Together, our study shows that influenza vaccination strategies, such as the cHA vaccine with the AS03 adjuvant, have the capability to robustly induce antibodies to the previously underappreciated anchor epitope and can provide broad protection against H1 viruses.

Any methods, additional references, Nature Research reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-021-04356-8. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Human peripheral blood mononuclear cells (PBMCs) and serum were obtained from multiple donors from multiple cohorts, which is outlined in Extended Data Table 1 

Human embryonic kidney HEK293T (female, #CRL-11268), Madin Darby canine kidney (MDCK; female, #CCL-34, NBL-2) and human A549 (#CCL-185) cells were purchased and authenticated by the American Type Culture Collection (ATCC). MDCK-SIAT1 cells were generated previously 37 . All cells were maintained in a humidified atmosphere of 5% CO 2 at 37 °C. HEK293T cells were maintained in advanced-DMEM supplemented with 2% ultra-low IgG fetal bovine serum (FBS; Invitrogen), 1% l-glutamine (Invitrogen) and 1% antibiotic-antimycotic (Invitrogen). MDCK, MDCK-SIAT1 and A549 cells were maintained in DMEM supplemented with 10% FBS (Invitrogen), 1% l-glutamine (Invitrogen) and 1% penicillin-streptomycin (Invitrogen). Jurkat cells expressing FcgRIIIa and FcgRI with a NFAT-driven luciferase reporter gene (#G7010) were acquired and validated by Promega and were directly used for antibody-dependent cellular cytotoxicity (ADCC) assays. Cell lines were not authenticated after receiving from suppliers and were not tested for mycoplasma.

Monoclonal antibodies were generated as previously described [38] [39] [40] . Peripheral blood was obtained from each donor approximately 7 days after vaccination or infection or obtained 28 + days post-vaccination. Lymphocytes were isolated and enriched for B cells using RosetteSep. Total plasmablasts (CD3 − CD19 + CD27 hi CD38 hi ; all cohorts except 2014 quadrivalent inactivated influenza vaccine (QIV)), IgG + plasmablasts (CD3 − CD19 + IgM − CD27 hi CD38 hi IgG + IgA − ; 2014 QIV), IgA + plasmablasts (CD3 − CD19 + IgM − CD27 hi CD38 hi IgG − IgA + ; 2014 QIV cohort), or HA + bait-sorted MBCs (CD3 − CD19 + CD27 + CD38 lo/+ HA + ; for 030-09M 1B06) were single-cell sorted into 96-well plates. Genes encoding immunoglobulin heavy and light chains were amplified by reverse transcriptase PCR (RT-PCR), sequenced, cloned into human IgG1, human κ-chain, or human λ expression vectors, and co-transfected into HEK293T cells. Secreted mAbs were purified from the supernatant using protein A agarose beads. For mAbs generated from the 2014 QIV cohort, mAb names include the original isotype of the sorted plasmablast, and all mAbs were expressed as human IgG1. cH5/1-binding B cells (CD19 + CD27 + cH5/1 + ) were sorted from donors 28 days after cH5/1 vaccination (NCT03300050). Cells were sorted with A/California/04/2009 HA probe (for 030-09M 1B06) or cH5/1 probe, each with a Y98F mutation to ablate non-specific binding to sialic acids on B cells. mAb heavy chain and light chain sequences were synthesized from single-cell RNA sequencing data of cH5/1-baited B cells (IDT), and cloned into the human IgG1, human κ-chain or human λ expression vectors. B cell clones were determined by aligning all the V(D)J sequences sharing identical progenitor sequences, as predicted by IgBLAST using our in-house software, VGenes. Consensus sequence analysis was performed using WebLogo 41 and sequence alignments were determined using Clustal Omega.

Influenza viruses used in all assays were grown in-house in specific pathogen free (SPF) eggs, harvested, purified and titred. Recombinant HA, cHA and mini-HA were obtained from BEI Resources or generated in-house. Recombinant HA mutant proteins used in Extended Data Fig. 2 were generated with identified mutations from the deep mutational scanning experiments (see below) or with known mutations that have arisen naturally or were identified in other studies 2,10,42-52 (Extended Data Table 3 ). All mutations were made on HA from A/California/7/2009 and were expressed in HEK293T cells and purified using Ni-NTA agarose beads (Qiagen).

High protein-binding microtitre plates (Costar) were coated with 8 haemagglutination units (HAU) of virus in carbonate buffer or with recombinant HA, including HA mutants described below, at 2 μg ml −1 in phosphate-buffered saline (PBS) overnight at 4 °C. Plates were washed the next morning with PBS 0.05% Tween and blocked with PBS containing 20% FBS for 1 h at 37 °C. Antibodies were then serially diluted 1:3 starting at 10 μg ml −1 and incubated for 1.5 h at 37 °C. Horseradish peroxidase (HRP)-conjugated goat anti-human IgG antibody diluted 1:1,000 ( Jackson Immuno Research) was used to detect binding of mAbs, and plates were subsequently developed with Super Aquablue ELISA substrate (eBiosciences). Absorbance was measured at 405 nm on a microplate spectrophotometer (Bio-Rad). To standardize the assays, control antibodies with known binding characteristics were included on each plate, and the plates were developed when the absorbance of the control reached 3.0 optical density (OD) units. All ELISAs were performed in duplicate twice. To define antibodies as targeting the H1 stalk, mAbs were tested for binding to cH5/1, which utilizes the head domain from H5-expressing viruses and the stalk domain from the pH1N1 virus 53 , and for haemagglutination inhibition (HAI) activity against pH1N1 (A/California/7/2009). mAbs that bound the cHA and lacked HAI activity were classified as those binding the HA stalk domain. To classify antigen specificity, mAbs that did not definitively bind to the HA head or stalk are listed as binding unknown HA + epitopes. Affinity measurements, as represented as dissociation constant (K d ) at a molar concentration (M), were calculated using Prism 9 (GraphPad) by performing a non-linear regression. All experiments were performed in duplicate and technically replicated twice.

Plates were coated with 50 μl of A/California/7/2009 HA at a concentration of 1 μg ml −1 and incubated overnight at 4 °C. To biotinylate the antibodies with known epitope specificities, CR9114 (CS epitope) and 047-09-4F04 (anchor epitope) were incubated at 4 °C with EZ-Link Sulfo-NHS-Biotin (Thermo Scientific) for 24 h or 48 h before use, respectively. After blocking the plates with PBS containing 20% FBS for 1 h at 37 °C, serum samples were incubated (starting dilution of 1:50 for human serum or 20 μg ml −1 for mAbs) in the coated wells for 2 h at room temperature. Either biotinylated CR9114 or 047-09 4F04 was then added at a concentration equal to twice its K d and incubated in the wells with the serum or mAbs for 2 h at room temperature. The biotinylated antibodies were desalted before addition to remove free biotin using Zeba spin desalting columns, 7 k MWCO (Thermo Scientific). After washing the plates, wells were incubated with HRP-conjugated streptavidin (Southern Biotech) at 37 °C for 1 h for detection of the biotinylated antibody. Super Aquablue ELISA substrate (eBiosciences) was then added, and absorbance was measured at 405 nm on a microplate spectrophotometer (Bio-Rad). To standardize the assays, biotinylated CR9114 or 047-09 4F04 was incubated in designated wells on each plate without any competing serum or mAb, and data were recorded when the absorbance of these wells reached an OD of 1-1.5 units. After subtracting background, percent competition by serum samples was then determined by dividing the observed OD of a sample by the OD reached by the positive control, subtracting this value from 1, and multiplying by 100. For the serum data, ODs were log transformed and analysed by non-linear regression to determine EC 50 values using Prism software (GraphPad). For Fig. 4 and Extended Data Fig. 5 , only donors with serum for all timepoints were included. All experiments were performed in duplicate and technically replicated twice.

High-protein binding microtitre plates (Costar) were coated with 10 μg ml −1 calf thymus double-stranded DNA (dsDNA; Thermo Fisher Scientific), 2 μg ml −1 Salmonella enterica serovar Typhimurium flagellin (Invitrogen), 5 μg ml −1 human insulin (Sigma-Aldrich), 10 μg ml −1 KLH (Invitrogen) and 10 μg ml −1 Escherichia coli LPS (Sigma-Aldrich) in PBS. Plates were coated with 10 μg ml −1 cardiolipin in 100% ethanol and allowed to dry overnight. Plates were washed with water and blocked with PBS, 0.05% Tween and 1 mM EDTA. mAbs were diluted 1 μg ml −1 in PBS and serially diluted fourfold and added to plates for 1.5 h. Plates were washed and goat anti-human IgG-HRP ( Jackson Immunoresearch) was diluted 1:2,000 in PBS, 0.05% Tween and 1 mM EDTA. Plates were washed with water and were blocked with PBS, 0.05% Tween and 1 mM EDTA for 5 min. Plates were washed again with water and were developed with Super Aquablue ELISA substrate (eBioscience) until the positive control mAb, 3H9 (ref. 54 ), reached an A 450 of 3. All experiments were performed in duplicate and technically replicated twice.

The mutant libraries used herein were previously described 55 . The libraries consist of all single amino acid mutations to A/WSN/1933 (H1N1). The experiments were performed by using biological triplicate libraries. The mutational antigenic profiling of 045-09 2B06, a CS epitope-binding mAb, was performed as previously outlined 56 . In brief, 10 6 tissue culture infectious dose 50 (TCID 50 ) of two of the virus library biological replicates was diluted in 1 ml in IGM (Opti-MEM supplemented with 0.01% FBS, 0.3% BSA and 100 mg ml −1 calcium chloride) and incubated with an equal volume of 045-09 2B06 antibody at a final concentration of 50 or 25 μg ml −1 for 1.5 h at 37 °C. MDCK-SIAT1 cells were infected with the virus antibody mixtures. Two hours post-infection, the media were removed, the cells washed with 1 ml PBS, and 2 ml of fresh IGM was added. Fifteen hours post-infection, viral RNA was extracted, reverse-transcribed using primers WSNHA-For (5′-AGCAAAAGCAGGGGAAAATAAAAACAAC-3′) and WSNHA-Rev (5′-AGTAGAAACAAGGGTGTTTTTCCTTATATTTCTG-3′), and PCR amplified according to the barcoded-subamplicon library preparation as previously described 55 . The overall fraction of virions that survive antibody neutralization was estimated using quantitative RT-PCR (qRT-PCR) targeting the viral nucleoprotein (NP) and cellular GAPDH as previously described 56 . Using tenfold serial dilutions of the virus libraries, we infected cells with no antibody selection to serve as a standard curve of infectivity. qPCR Ct values from the standard curve samples compared to the virus-antibody mix samples were determined for NP and GAPDH. We then generated a linear regression to fit the difference between the NP and GAPDH Ct values for the standard curve samples, and then used this curve to interpolate the fraction surviving for the antibody-virus selection samples. Across the three library replicates, the fraction of virus surviving antibody selection was 0.17, 0.1 and 0.14.

Illumina(R) deep sequencing data were analysed using dms_tools2 version 2.4.12 software package 57 , which can be found at https://github. com/jbloomlab/dms_tools2. The computer code used is at https:// github.com/jbloomlab/2B06_DMS, and the Jupyter notebook that performed most of the analysis is at https://github.com/jbloomlab/2B06_ DMS/blob/master/analysis_notebook.ipynb. The sequencing counts were processed to estimate the differential selection for each mutation, which is the log enrichment of that mutation in the antibody-selected condition versus the control 56 . The numerical measurements of the differential selection that 2B06 imposes on each mutation can be found at: https://github.com/jbloomlab/2B06_DMS/blob/master/results/ diffsel/tidy_diffsel.csv.

Amino acid preferences for the HA of A/WSN/1933 (H1N1) were previously determined 55 . In brief, deep mutational scanning was performed by passaging virus libraries at a low multiplicity of infection in MDCK-SIAT1 cells. Following deep sequencing of the resulting virus, amino acid preferences were determined using the Python package dms_tools (http://jbloomlab.github.io/dms_tools/), version 1.1.12. This program aligns subamplicon reads to the reference HA sequence, counts the number of mutations at each amino acid site, and determines amino acid preferences based on the mutation counts pre-selection and post-selection.

Microneutralization assays for mAb characterization were carried out as previously described 58, 59 . MDCK cells were maintained in DMEM supplemented with 10% FBS, 1% penicillin-streptomycin, and 1% l-glutamine at 37 °C with 5% CO 2 . The day before the experiment, 25,000 MDCK cells were added to each well of a 96-well plate. Serial twofold dilutions of mAb were mixed with an equal volume of 100 TCID 50 of virus for 1 h and added to MDCK cells for 1 h at 37 °C. The mixture was removed, and cells were cultured for 20 h at 37 °C with 1X MEM supplemented with 1 μg ml −1 tosyl phenylalanyl chloromethyl ketone (TPCK)-treated trypsin and appropriate mAb concentration. Cells were washed twice with PBS, fixed with 80% ice-cold acetone at −20 °C for at least 1 h, washed three times with PBS, blocked for 30 min with 3% BSA-PBS, and then treated for 30 min with 2% H 2 O 2 . Cells were incubated with a mouse anti-NP antibody (1:1,000; Millipore) in 3% BSA-PBS for 1 h at room temperature, followed by goat anti-mouse IgG HRP (1:1,000; Southern Biotech) in 3% BSA-PBS for 1 h at room temperature. The plates were developed with Super Aquablue ELISA substrate at 405 nm until virus-only controls reached an OD of 1. The signal from uninfected wells was averaged to represent 100% inhibition. The signal from infected wells without mAb was averaged to represent 0% inhibition. Duplication wells were used to calculate the mean and s.d. of neutralization, and the half-maximal inhibitory concentration (IC 50 ) was determined by a sigmoidal dose-response curve. The inhibition ratio (%) was calculated as: ((OD positive control -OD sample)/ (OD positive control -OD negative control)) × 100. The final IC 50 was determined using Prism software (GraphPad). All experiments were performed in duplicate and technically replicated twice.

Twenty thousand MDCK cells were seeded per well in a 96-well cell culture plate (Corning) and the cells were used the next morning to perform the neutralization assay. Antibody dilutions were prepared starting at 30 μg ml −1 with threefold subsequent dilutions in 1X MEM. Each respective dilution was mixed with 10,000 plaque-forming units (PFU) of cold-adapted A/Ann Arbor/6/1960 (H2N2) virus for 1 h at room temperature. After 1 h, cells were washed with PBS and 100 μl of antibody-virus mixture was added onto the cells for 1 h at 37 °C. Next, the antibody-virus mixture was removed and 60 μl of 1X MEM containing TPCK was added to each well. Of each corresponding antibody dilution, 60 μl was also added to each well and the cells were incubated at 33 °C for 3 days. On the third day, a haemagglutination assay was performed using turkey red blood cells to assess the HAU at each antibody dilution.

In vivo challenge infections mAb cocktails (Extended Data Fig. 2b) were passively transferred into 6-8-week-old female BALB/c mice ( Jackson Laboratories) by intraperitoneal injection of 0.2, 1 and 5 mg per kg mAb cocktail, which are further detailed in Supplementary Table 1 . Negative control mice received 5 mg per kg of the anthrax-specific mAb 003-15D03 as an isotype control. mAbs were administered 2 h before infection for prophylactic treatment and 48 h post-infection for therapeutic treatment. For prophylactic mAb studies with A/Netherlands/602/2009 (Extended Data Fig. 3a) , mice were anaesthetized with isoflurane and intranasally challenged with 10 lethal dose 50 (LD 50 ) of mouse-adapted A/Netherlands/602/2009 H1N1 virus, with 10 μl of virus administered into each nostril (20 μl total). For therapeutic treatment of A/Netherlands/602/2009 and prophylactic treatment of A/Fort Monmouth/1/1947, mice were anaesthetized with a ketamine-xylazinewater cocktail (0.15 mg ketamine per kg and 0.03 mg per kg xylazine; 100 ml intraperitoneally) and infected with 10 LD 50 of A/Netherlands/602/2009 or A/Fort Monmouth/1/1947. As a read out, survival and weight loss were monitored 1-2 times daily for 2 weeks. Mice were euthanized upon 25% weight loss or at the end of the experiment (14 days post-challenge). Five mice per condition per experiment with two biological replicates were utilized based on a previously performed power analysis. Data were pooled for analysis.

To determine differences in lung viral load, 5 mg per kg of the antibody cocktails was administered prophylactically as described above. Two hours after mAb administration, mice (n = 5 mice per group) were anaesthetized and intranasally challenged with 1 LD 50 of A/Netherlands/602/2009. Lungs were collected at day 3 and day 6 post-infection, homogenized and viral load was determined via plaque assay. All experiments were done in accordance with the University of Chicago and Icahn School of Medicine at Mount Sinai Institutional Animal Care and Use Committees. Animals studies were not blinded or randomized.

For determination of viral load in mouse lung tissues a standard plaque assay was performed. Confluent monolayers of MDCK cells were infected with serial dilutions of homogenized lung tissue ranging from 1:10 to 1:1,000,000 diluted in 1X MEM (1% penicillin-streptomycin antibiotics mix, 1% HEPES, 1% l-glutamine and 1% sodium-bicarbonate (Gibco)) for 1 h at 33 °C, with shaking every 15 min. Afterwards, an overlay containing 2% Oxoid agar (Thermo Fisher), H 2 O, 2X MEM, DEAE and TPCK-treated trypsin was added to the cells. The plates were incubated at 33 °C for 3 days and then fixed with 3.7% paraformaldehyde overnight at 4 °C. Plaques were visualized by immunostaining. Here, the agar overlay was removed and the plates blocked with 3% milk and PBS. The blocking solution was removed and primary antibody ((H1N1 guinea pig anti-sera (generated in house)) diluted 1:3,000 in 1% milk and PBS was added for 1 h. The plates were washed three times with PBS and secondary antibody (anti-mouse IgG H&L peroxidase-conjugated (Rockland) diluted 1:3,000 in 1% milk and PBS was added for 1 h. The plates were washed three times with PBS and developed by using KPL TrueBlue Peroxidase Substrate (SeraCare).

Antibody-dependent cellular cytotoxicity reporter assay A549 cells were maintained in DMEM supplemented with 10% FBS, 10 U ml −1 penicillin, and 10 mg ml −1 streptomycin) and were plated in 96-well, white-walled plates (Costar) at 2.5 × 10 5 cells per ml overnight at 37 °C with 5% CO 2 . The following day, cells were washed with PBS and infected with A/Netherlands/602/2009 at a multiplicity of infection of 5 in UltraMDCK media (Lonza) for 24 h in the absence of TPCK-treated trypsin. mAbs were serially diluted in assay buffer (RPMI 1640 supplemented with 4% ultra-low IgG FBS; Gibco), starting at 60 μg ml −1 and diluted threefold. Cell medium was aspirated and 25 μl of assay buffer and 25 μl of diluted antibody were added to each well. Jurkat cells expressing human FcgRIIIa with a NFAT-driven luciferase reporter gene (Promega) were diluted to 3 × 10 6 cells per ml, 25 μl of cells was added to each well and incubated at 37 °C with 5% CO 2 for 6 h. Plates were removed from the incubator and placed at room temperature for 15 min. Of the BioGlo luciferase substrate (Promega), 75 μl was added to each well and luminescence was read immediately using a Syngery H1 hybrid multimode microplate reader (Biotek). EC 50 values were determined using Prism 8 (GraphPad).

The footprints of three mAbs (FISW84 (PDB: 6HJQ), CR9114 (PDB: 4FQI) and FI6v3 (PDB: 3ZTN)) were mapped onto one HA protomer (A/California/4/2009, PDB: 4M4Y) using UCSF Chimera 60 and Adobe Photoshop. Negative-stain EM maps of HA-Fab complexes were aligned in UCSF Chimera and estimated footprints were mapped onto one HA protomer. Individual protomers of the HA trimer are indicated in different shades of grey.

Immune complexes were prepared by incubating Fab with HA (A/California/04/2009 with E376K or E376G stabilizing mutations) at greater than 3:1 molar ratio for 2 h at room temperature. Samples were deposited at approximately 10 μg ml −1 on glow-discharged, carbon-coated 400 mesh copper grids (Electron Microscopy Sciences) and stained with 2% w/v uranyl formate. Samples were imaged at ×52,000 magnification, 120 kV, on a Tecnai Spirit T12 microscope equipped with an Eagle CCD 4k camera (FEI) or ×62,000 magnification, 200 kV, on a Tecnai T20 microscope equipped with a CMOS 4k camera (TVIPS). Micrographs were collected with Leginon, single particles were processed with Appion, Relion and XQuartz, and footprints were mapped with UCSF Chimera, and figures were made with UCSF Chimera 60-63 .

222-1C06 and 045-09 2B05 Fabs were incubated at greater than 3:1 molar ratio with HA (A/California/7/2009, E376K) for 1 h at room temperature. 045-09 2B05 Fab, targeting the lateral patch, was added to the immune complex to induce particle tumbling and increase angular sampling 3 . Using a Thermo Fisher Vitrobot, the immune complex (0.5 mg ml −1 ) incubated with lauryl maltose neopentyl glycol (5 μM, Anatrace) was deposited onto glow-discharged Au 1.2/1.3 300 mesh grids (Electron Microscopy Sciences), blotted for 7 s, and plunge-frozen in liquid ethane. Samples were imaged at ×36,000 nominal magnification on a 200 kV Talos Arctica electron microscope (FEI) with a CETA 4k CMOS camera (FEI, total dose 49.92 e/Å 2 ) and Gatan K2 Summit detector in counting mode. 2,243 micrographs were collected, aligned and CTF-corrected using Leginon, MotionCor2 in Appion, and Patch-CTF in CryoSPARC2, respectively 61, 62, 64, 65 . In CryoSPARC2, particles were picked using apo HA templates, selected through reference-free 2D classification, and reconstructed through 3D classification and refinement. The final map resolved to a global 3.75 Å resolution with C3 symmetry and 44,224 particles. Figures were made in Prism 8 (GraphPad) and UCSF Chimera 60 .

A predicted model of 222-1C06 Fab was generated using abYsis (http:// www.abysis.org/abysis/) and docked into EM density along with an initial model of CA09 H1 HA + 045-09 2B05 (PDB: 7MEM). The initial model was iteratively refined using COOT and Rosetta 66, 67 . The final model was numbered using the H3 and Kabat numbering schemes. The final model and map were evaluated using MolProbity, EMRinger 68,69 , Phenix and the PDB validation server. After modelling the immune complex, we segmented the Fab density from HA in the cryo-EM map and mapped the footprint of the 222-1C06 model in the HA density. Cryo-EM data collection and refinement statistics are included in Extended Data Table 4 .

Human serum samples were heat-inactivated at 55 °C for 30 min before incubating on Capture Select IgG-Fc (ms) Affinity Matrix (Fisher) to Article bind IgG at 4 °C for 72 h on a rotator. Samples with IgG bound to resin were centrifuged at 4,000 rpm and supernatant was collected. IgG samples were washed three times with PBS followed by centrifugation to remove supernatant. Samples were buffer exchanged into buffer containing 100 mM Tris, 2 mM EDTA and 10 mM l-cysteine through centrifugation with Amicon filters, then incubated with papain for 4 h at 37 °C, shaking at 80 rpm. The digestion reactions were quenched with 50 mM iodoacetamide, buffer exchanged to TBS and separated by size-exclusion chromatography (SEC) with a Superdex 200 increase 10/300 column (GE Healthcare). Fab and undigested IgG were collected and concentrated, and 500 μg Fab was complexed with 10 μg HA for 18 h at room temperature. Reactions were purified by SEC and immune complexes were collected and concentrated. Negative-stain EM grids were prepared as described above.

HEK293T cells were plated into a six-well plate and transfected overnight with 0.2 μg of plasmid and 10 μg ml −1 PEI. After 12-16 h, media were replaced with PFHM-II (Gibco) and cells were rested for 3 days. Transfected cells were trypsinized, washed and aliquoted. Cells were stained with 10 μg ml −1 of individual mAbs for 30 min. Cells were washed and stained with anti-human IgG Fc-BV421 for 30 min. Cells were washed two times and run on a BD LSRFortessa and collected with BD FACSDiva software. Data were analysed using FlowJo v10.

Single-cell RNA sequencing and repertoire analysis cH5/1 + memory B cells (CD19 + CD27 + HA + ) were bulk sorted and partitioned into nanolitre-scale gel bead-in-emulsions (GEMs) to achieve single-cell resolution using the 10x Genomics Chromium Controller and according to the manufacturer's instruction (10x Genomics). The sorted single cells were processed according to 5′ gene expression and B cell immunoglobulin enrichment instruction to prepare the libraries for sequencing. Libraries were sequenced using an Illumina HiSeq 4000 at Northwestern University or an Illumina NextSeq 500 at the University of Chicago. Cellranger Single-Cell Software Suite (version 3.0) was used to perform sample demultiplexing, barcode processing, and single-cell 5′ and V(D)J counting, and Cellranger mkfastq was used to demultiplex raw base call (BCL) files into sample-specific fastq files. Subsequently, GRCh38-1.2.0 and cellranger-vdj-GRCh 38-alta-ensembl-2.0.0 were used as references for the transcriptome and V(D)J assembly, respectively. Cellranger counts and Cellranger vdj package were used to identify gene expression and assemble V(D) J pairs of antibodies.

Single-cell datasets were analysed using Seurat 3 toolkit (Version 3.2.0). We performed conventional pre-process steps for all 20 donors including cell quality control (QC), normalization, identification of highly variable features, data scaling and linear dimensional reduction. More specifically, we only kept cells with more than 200 and less then 2,500 detected genes for the QC step. We also filtered out cells with high mitochondrial gene expression using a 'softThreshold' function in the R package LinQ-View (version 0.99) 70 . We normalized the RNA data using conventional log normalization. We identified 2,000 highly variable genes for each dataset and performed principle component analysis (PCA) in linear dimensional reduction step. We then integrated all 20 single-cell datasets from vaccinated participants to remove batch effects using the Anchor method in Seurat 3. In this analysis, we filtered our dataset and only kept cells with both transcriptome and full length and paired heavy and light chain V(D)J sequences (n = 1,952). From these cells, we identified a group of 'VH1-69/κ' B cells that used the VH1-69 gene and κ-light chain, which is enriched for B cells targeting the BN stalk epitope. We also identified a group of 'anchor epitope'-specific B cells by the following rules: (1) VH locus: VH3-23, VH3-30, VH3-30-3 or VH-3-48; (2) VK locus: VK3-11 or VK3-15; (3) JK locus: JK4 or JK5; (4) K-CDR3 length equal to 10; and (5) a 'NWP' pattern in the K-CDR3 peptide.

Pan-H1 conservation models are based on consensus strains (listed in Supplementary Table 4 ) of distinct H1 clades isolated from humans, swine and avian sources, as described in Zhuang et al. 71 and inclusion of the Eurasian swine-like A/swine/Jiangsu/J004/2018 (ref. 23 ). To generate the group 1 HA conservation model, we selected one representative sequence for each group 1 HA subtype from FluDB (https://www. fludb.org/; Supplementary Table 5 ) according to a previous study 72 . A multiple sequence alignment from these HA protein sequences was generated using MUSCLE 73 and the conservation of each residue was quantified using an entropy model 41 . Seasonal H1 conservation models are based on consensus strains of H1N1 viruses (59 strains total) circulating between 1918-1957 and 1976-2019, which was previously described 3 . Amino acid alignments and H3 numbering were performed using Librator 74 and Burke and Smith HA numbering 72 .

To predict the structures of the investigated Fv fragments (222-1C06, FISW84, 241 IgA 2F04 and SFV009 3G01) with A/California/4/2009 E47G HA (PDB: 7MEM), we applied the program RosettaAntibody 67, 75, 76 . The Fvs were protonated using the Protonate 3D tool 77, 78 . Charge neutrality was ensured by utilizing the uniform background plasma approach in AMBER 79, 80 . Using the tleap tool of the AmberTools20 (ref. 81 ) package, the structure models were soaked in cubic water boxes of TIP3P water molecules with a minimum wall distance of 10 Å to the protein 82 . Parameters for all antibody models derive from the AMBER force field 14SB 83 . The Fvs were carefully equilibrated using a multistep equilibration protocol 84 .

To enhance the sampling of the conformational space, well-tempered bias-exchange metadynamics [85] [86] [87] simulations were performed in GROMACS 88, 89 with the PLUMED 2 implementation 90 . We chose metadynamics as it enhances sampling on predefined collective variables. The sampling is accelerated by a history-dependent bias potential, which is constructed in the space of the collective variables 85, 86, 91 . As collective variables, we used a well-established protocol, boosting a linear combination of sine and cosine of the ψ torsion angles of all six CDR loops calculated with functions MATHEVAL and COMBINE implemented in PLUMED 2 (ref. 90 ). As discussed previously, the ψ torsion angle captures conformational transitions comprehensively 92 . The underlying method presented in this paper has been validated in various studies against a large number of experimental results 93 . The simulations were performed at 300 K in an NpT ensemble using the GPU implementation of the pmemd module 94 to be as close to the experimental conditions as possible and to obtain the correct density distributions of both protein and water. We used a Gaussian height of 10.0 kJ mol −1 and a width of 0.3 rad. Gaussian deposition occurred every 1,000 steps and a biasfactor of 10 was used. 500 ns of bias-exchange metadynamics simulations were performed for the prepared Fv structures. The resulting trajectories were aligned to the whole Fv and clustered with the program cpptraj 80, 95 using the average linkage hierarchical clustering algorithm with a RMSD cut-off criterion of 1.2 Å resulting in a large number of clusters. The cluster representatives for the antibody fragments were equilibrated and simulated for 100 ns using the AMBER 20 (ref. 81 ) simulation package.

Molecular dynamics simulations were performed in an NpT ensemble using the pmemd.cuda module of AMBER 20 (ref. 80 ). Bonds involving hydrogen atoms were restrained with the SHAKE algorithm 96 , allowing a time step of 2.0 fs. Atmospheric pressure (1 bar) of the system was set by weak coupling to an external bath using the Berendsen algorithm 97 . The Langevin thermostat 98 was used to maintain the temperature during simulations at 300 K.

With the obtained trajectories, we performed a time-lagged independent component analysis (tICA) using the Python library PyEMMA 2, using a lag time of 10 ns. tICA was applied to identify the slowest movements of the investigated Fab fragments and consequently to obtain a kinetic discretization of the sampled conformational space 99 . On the basis of the tICA conformational spaces, thermodynamics and kinetics were calculated with a Markov-state model 100 of all six CDR loops by using PyEMMA 2. The resulting kinetically dominant ensemble in solution was further used to predict the interactions of H1 with the Fvs. To model the complex and to predict interactions in the binding interface, we used the crystal structure of the full-length influenza HA (PDB: 7MEM) as template structure. In addition, the obtained complex structure was further minimized and equilibrated.

All statistical analyses were performed using Prism software (GraphPad versions 8 and 9) or R. Sample sizes (n) for the number of mAbs tested are indicated in corresponding figures or in the centre of pie graphs. The number of biological repeats for experiments and specific tests for statistical significance used are indicated in the corresponding figure legends. P values less than or equal to 0.05 were considered significant: *P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001 and ****P < 0.0001.

Further information on research design is available in the Nature Research Reporting Summary linked to this paper.

Repertoire data generated from single-cell RNA sequencing data were deposited at Mendeley Data (https://data.mendeley.com/datasets/ jzsx489pmk/1). Accession numbers for all other anchor-targeting mAbs are included in Supplementary Table 6 . Electron microscopy maps were deposited to the Electron Microscopy Data Bank under accession IDs: EMD-25634-EMD-25646. The cryo-EM map and model of anchor and lateral patch Fabs binding H1 HA were deposited to the RCSB database with accession numbers EMD-25655/PDB 7T3D. All next-generation sequencing data for 045-09 2B06 deep mutational scanning and for the H1N1 mutational scanning can be found on the Sequence Read Archive under BioProject accession number PRJNA309339. The following Protein Data Bank accession numbers were downloaded and included in the paper: 6HJQ, 3SDY, 4NM8, 4M4Y, 4WE4, 4JTV, 4FQI, 3ZTN and 7MEM. Sera from the vaccine cohorts are unique to this study and are not publicly available. All source data are included with the paper. All other material is available on reasonable request to the corresponding authors. Source data are provided with this paper.

Corresponding author(s): Jenna Guthmiller, Andrew Ward, Patrick Wilson Last updated by author(s): Dec 3, 2021 Reporting Summary Nature Portfolio wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Portfolio policies, see our Editorial Policies and the Editorial Policy Checklist.

For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section.

The exact sample size (n) for each experimental group/condition, given as a discrete number and unit of measurement A statement on whether measurements were taken from distinct samples or whether the same sample was measured repeatedly

The statistical test(s) used AND whether they are one-or two-sided Only common tests should be described solely by name; describe more complex techniques in the Methods section.

A description of all covariates tested A description of any assumptions or corrections, such as tests of normality and adjustment for multiple comparisons A full description of the statistical parameters including central tendency (e.g. means) or other basic estimates (e.g. regression coefficient) AND variation (e.g. standard deviation) or associated estimates of uncertainty (e.g. confidence intervals)

For null hypothesis testing, the test statistic (e.g. F, t, r) with confidence intervals, effect sizes, degrees of freedom and P value noted For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Portfolio guidelines for submitting code & software for further information.

Policy information about availability of data All manuscripts must include a data availability statement. This statement should provide the following information, where applicable:

-Accession codes, unique identifiers, or web links for publicly available datasets -A description of any restrictions on data availability -For clinical datasets or third party data, please ensure that the statement adheres to our policy Repertoire data generated from single cell RNA-sequencing data is deposited at Mendeley Data (https://data.mendeley.com/datasets/jzsx489pmk/1). Accession numbers for all other anchor-targeting mAbs are included in Supplemental 

The pathway to a universal influenza vaccine

Pre-existing immunity to influenza virus hemagglutinin stalk might drive selection for antibody-escape mutant viruses in a human challenge model

First exposure to the pandemic H1N1 virus induced broadly neutralizing antibodies targeting hemagglutinin head epitopes

Immunogenicity of chimeric haemagglutinin-based, universal influenza virus vaccine candidates: interim results of a randomised, placebo-controlled, phase 1 clinical trial

A chimeric hemagglutinin-based universal influenza virus vaccine approach induces broad and long-lasting immunity in a randomized, placebo-controlled phase I trial

Novel correlates of protection against pandemic H1N1 influenza A virus infection

Pre-existing hemagglutinin stalk antibodies correlate with protection of lower respiratory symptoms in flu-infected transplant patients

Broadly cross-reactive antibodies dominate the human B cellresponse against 2009 pandemic H1N1 influenza virus infection

Immune history profoundly affects broadly protective B cell responses to influenza

Highly conserved protective epitopes on influenza B viruses

Influenza hemagglutinin membrane anchor

A highly conserved neutralizing epitope on group 2 influenza A viruses

A common solution to group 2 influenza virus neutralization

A neutralizing antibody selected from plasma cells that binds to group 1 and group 2 influenza A hemagglutinins

Broadly neutralizing anti-influenza antibodies require Fc receptor engagement for in vivo protection

Epitope specificity plays a critical role in regulating antibody-dependent cell-mediated cytotoxicity against influenza A virus

Polyreactive broadly neutralizing B cells are selected to provide defense against pandemic threat influenza viruses

In vitro neutralization is not predictive of prophylactic efficacy of broadly neutralizing monoclonal antibodies CR6261 and CR9114 against lethal H2 influenza virus challenge in mice

Structural transitions in influenza haemagglutinin at membrane fusion pH

Polyclonal epitope mapping reveals temporal dynamics and diversity of human antibody responses to H5N1 vaccination

Permissive secondary mutations enable the evolution of influenza oseltamivir resistance

An efficient method to generate monoclonal antibodies from human B cells

Rapid cloning of high-affinity human monoclonal antibodies against influenza virus

Rapid generation of fully human monoclonal antibodies specific to a vaccinating antigen

WebLogo: a sequence logo generator

How single mutations affect viral escape from broad and narrow antibodies to H1 influenza hemagglutinin

Antibody 27F3 broadly targets influenza A group 1 and 2 hemagglutinins through a further variation in VH1-69 antibody orientation on the HA stem

Functional balance of the hemagglutinin and neuraminidase activities accompanies the emergence of the 2009 H1N1 influenza pandemic

A broadly neutralizing anti-influenza antibody reveals ongoing capacity of haemagglutinin-specific memory B cells to evolve

Different genetic barriers for resistance to HA stem antibodies in influenza H3 and H1 viruses

A pan-H1 anti-hemagglutinin monoclonal antibody with potent broad-spectrum efficacy in vivo

Use of hemagglutinin stem probes demonstrate prevalence of broadly reactive group 1 influenza antibodies in human sera

Structural and genetic basis for development of broadly neutralizing influenza antibodies

A single amino acid in the stalk region of the H1N1pdm influenza virus HA protein affects viral fusion, stability and infectivity

Structural and functional bases for broad-spectrum neutralization of avian and human influenza A viruses

Novel inhibitors of influenza virus fusion: structure-activity relationship and interaction with the viral hemagglutinin

Influenza viruses expressing chimeric hemagglutinins: globular head and stalk domains derived from different subtypes

Structure and function of anti-DNA autoantibodies derived from a single autoimmune mouse

Accurate measurement of the effects of all amino-acid mutations on influenza hemagglutinin

Complete mapping of viral escape from neutralizing antibodies

Software for the analysis and visualization of deep mutational scanning data

Preexisting human antibodies neutralize recently emerged H7N9 influenza strains

Influenza infection in humans induces broadly cross-reactive and protective neuraminidase-reactive antibodies

UCSF Chimera-a visualization system for exploratory research and analysis

Automated molecular microscopy: the new Leginon system

Appion: an integrated, database-driven pipeline to facilitate EM image processing

RELION: implementation of a Bayesian approach to cryo-EM structure determination

MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy

cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination

Coot: model-building tools for molecular graphics

Modeling and docking of antibody structures with Rosetta

MolProbity: all-atom structure validation for macromolecular crystallography

EMRinger: side chain-directed model and map validation for 3D cryo-electron microscopy

Improved integration of single-cell transcriptome and surface protein expression by LinQ-View

Diversity and distribution of type A influenza viruses: an updated panorama analysis based on protein sequences

A recommended numbering scheme for influenza A HA subtypes

MUSCLE: multiple sequence alignment with high accuracy and high throughput

Librator, a platform for optimized sequence editing, design, and expression of influenza virus proteins

Blind prediction performance of RosettaAntibody 3.0: grafting, relaxation, kinematic loop modeling, and full CDR optimization

Improvements to robotics-inspired conformational sampling in Rosetta

Protonate3D: assignment of ionization states and hydrogen coordinates to macromolecular structures

Molecular Operating Environment (MOE)

Quantifying artifacts in Ewald simulations of inhomogeneous systems with a net charge

3rd PTRAJ and CPPTRAJ: software for processing and analysis of molecular dynamics trajectory data

Comparison of simple potential functions for simulating liquid water

ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB

A challenging system: free energy prediction for factor Xa

Well-tempered metadynamics: a smoothly converging and tunable free-energy method

Metadynamics enhanced Markov modeling of protein dynamics

Article

GROMACS: high performance molecular simulations through multi-level parallelism from laptops to supercomputers

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

PLUMED 2: new feathers for an old bird

Well-tempered metadynamics as a tool for characterizing multi-component, crystalline molecular machines

Stereochemistry of polypeptide chain configurations

Antibodies exhibit multiple paratope states influencing VH-VL domain orientations

Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh Ewald

Clustering molecular dynamics trajectories: 1. Characterizing the performance of different clustering algorithms

Settle: an analytical version of the SHAKE and RATTLE algorithm for rigid water models

Molecular-dynamics with coupling to an external bath

Generalized Langevin equation approach for atom/solidsurface scattering: general formulation for classical scattering off harmonic solids

PyEMMA 2: a software package for estimation, validation, and analysis of Markov models

Markov state models of biomolecular conformational dynamics

MIV -2009 pandemic H1N1 monovalent inactivated influenza vaccine

trivalent inactivated influenza vaccine

QIV -2014-2015 quadrivant inactivated influenza vaccine

LAIV -live-attenuated influenza vaccine; IIV -inactivated influenza vaccine; AS03 -Adjuvant System 03

All next generation sequencing data for 045-09 2B06 deep mutational scanning and for the H1N1 mutational scanning can be found on the Sequence Read Archive under BioProject accession number PRJNA309339. The following Protein Database accession numbers were downloaded and included in the manuscript -6HJQ

Replication All experiments were performed more than once, except the mouse viral titers. All replications were successful

Randomization Recipients for the cHA vaccine arms were randomized. Otherwise, no other experimental groups were used and therefore were not randomized

Blinding Serum competition ELISAs in Fig. 4 and Extended Data 6 were blinded. Otherwise, no other blinding was performed

Reporting for specific materials, systems and methods

We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, system or method listed is relevant to your study

Anti-human IgG Fc -BV421, clone M1310G05, Biolegend, Cat#410703 Mouse anti-influenza A nucleoprotein antibody, clone A3, biotin, Sigma/Millipore

Polyclonal Goat anti-mouse IgG HRP, Southern Biotech

Polyclonal Goat anti-mouse IgG H&L peroxidase-conjugated

Monoclonal antibodies from previous publications and synthesized in-house are listed below with reference

Monoclonal antibodies unique to this study and synthesized in-house are listed

Validation All commercial antibodies were validated by their manufacturers and were titrated in the lab to determine optimal concentration for experimentation. In-house generated monoclonal antibodies were validated in preliminary ELISAs to A/California/7/2009 H1N1 virus or recombinant H1 or cH5/1 proteins. MAb concentrations were standardized based on the assay and starting concentration is described in methods section

HEK293T cells -ATCC MDCK cells -ATCC A549 cells -ATCC Jurkat cells expressing FcgRIIIa with NFAT-drive luciferase reporter gene -Promega

No other authentication at the lab level was performed. 8 week old female Balb/c mice were used for challenge studies. Mice were housed in ABSL-2 conditions with 12-hour light/dark cycles and controlled temperature and humidity

Animal experiments were approved by the University of Chicago and Icahn School of Medicine at Mount Sinai IACUCs

Only birth date, prior influenza virus vaccination status, and sex were obtained or released to the authors for donors within these cohorts. For the cHA vaccine trial -In Groups 1, 2, and 4 about two-thirds of the subjects were female (70.0%, 66.7%, and 62.5%, respectively), compared to 2/5 (40.0%) in Group 3 and 5/10 (50.0%) in Group 5. Most were non-Hispanic or Latino (80% to 100% per group) and black or African American (66.7% to 87.5%). The median age at enrollment ranged from 26 to 31 years across groups, with minimum age ranging from 18 to 22 and maximum age from 29 to 38. The median weight ranged from 74.7 kg to 84.2 kg and the median height ranged from 165.5 cm to 174.0 cm. The most common pre-existing conditions across all treatment groups were immune system disorders (24/66, 36.4%), the majority being allergies (seasonal and food), followed by drug hypersensitivity. No subjects in the placebo Groups 3 and 5 reported taking prior medications and more subjects in Group 1 (5/19, 26.3%) reported taking prior medications

QIV, and 2009 pH1N1 infection had been recruited from prior studies from the local community. Participants for the cHA vaccine trial were recruited according to Bernstein et al. Lancet Infectious Disease 2020 and Nachbagauer et al. Nature Medicine 2020.The target population reflected the community at large at each of the participating study sites. Information regarding this trial was provided to potential subjects who have previously participated in vaccine trials conducted at the participating study sites. Ethics oversight The study protocol was approved by the IRB at the University of Chicago for all studies, the Emory University for the 2009 pH1N1 infection and 2009 MIV cohorts, and at Icahn School of Medicine at Mount Sinai

All manuscripts should comply with the ICMJE guidelines for publication of clinical research and a completed CONSORT checklist must be included with all submissions. Clinical trial registration NCT03300050 Study protocol Study Protocol can be accessed here

NC 27710 Studied period: Date of first enrollment

Outcomes Primary Objectives To assess the reactogenicity and safety through 28 days after each priming dose of cH8/1Nl LAIV (or placebo) and the booster dose of cH5/1Nl IIV +/-AS03A (or placebo) and through 28 days after each dose of IIV (cH8/1Nl IIV + AS03A and cH5/1Nl IIV + AS03A) (or placebo) in terms of rates of solicited local and general adverse events (AEs) through 7 days post-vaccination, unsolicited AEs through 28 days post vaccination, hematological and biochemical laboratory abnormalities up to Visit 13, and medically attended event (MAEs), laboratory-confirmed influenza-like illness (LC-ILi), potential immune-mediated disease (plMDs), and serious adverse events (SAEs) through Visit 13. Secondary Objectives For this study, serum antibody titers (EC50) for the anchor and central stalk epitopes was determined

Flow Cytometry Plots Confirm that: The axis labels state the marker and fluorochrome used (e.g. CD4-FITC)

Include numbers along axes only for bottom left plot of group (a 'group' is an analysis of identical markers)

All plots are contour plots with outliers or pseudocolor plots. A numerical value for number of cells or percentage

Methodology Sample preparation HEK293T cells were transfected with plasmids to express full length membrane-bound A/California/7/2009 H1 or membrane bound mini-HA. Cells were harvested 4 days later and stained with monoclonal antibodies of interest

Instrument BD LSRFortessa Software Data Collection -FACSDIva Anaylsis -FlowJo v10

Cell population abundance Cell population abundance is shown for representative panels and gates are identical across individual datasets. Gating strategy Cells were gated on using FSC and SSC. From this gate, mAb+ cells were gated on, based on the secondary Ab only stain. Gating strategy and examples are in Supplemental Fig. 1. For mAb+ gating

Tick this box to confirm that a figure exemplifying the gating strategy is provided in the Supplementary Information

We thank all of the volunteers who participated in this study; S. Andrews, R. Ahmed, J. Wrammert 

Extended Data Fig. 1 | Binding and neutralization features of anchor epitope-binding mAbs. Related to Fig. 1 . a, Proportion of HA + mAbs binding to distinct HA domains (left) and proportion of stalk-binding mAbs binding the CS domain (right). Number in the center of the pie graphs represent the number of mAbs tested. b, Proportion of mAbs per cohort that bind the HA stalk domain. c, Negative stain 2D class averages of 047-09 4F04, 241 IgA 2F04, and 222-1C06 binding to H1 (A/California/7/2009 E376K HA). Imaging of 047-09 4F04 was performed at ×52,000 normal magnification and of 222-1C06 and 241 IgA 2F04 at ×62,000 normal magnification. d, Overlay of 047-09 4F04, 241 IgA 2F04, 222-1C06, and FISW84 (PDB:6HJQ) Fabs binding the anchor epitope of A/California/4/2009 HA. e, Overlay of CR8020 (PDB:3SDY), CR8043 (PDB:4NM8), and FISW84 (PDB:6HJQ) modeled on A/California/7/2009 E376G (PDB:4M4Y). f, Footprints of anchor mAb 222-1C06 on H1 (top; PDB: 4M4Y) and CR8020 and CR8043 on H3 (bottom; PDB:4WE4). g, Heatmap of apparent affinity (K d ; M) of anchor-targeting mAbs binding to historical and recent H1N1 viruses. h, Neutralization potency of anchor-binding mAbs (n = 15) against H1expressing viruses. i, Representative microneutralization curves of anchor-(n = 42) and CS-binding (n = 29) mAbs against A/California/7/2009. j, IC 80 of anchor-and CS-binding mAbs against A/California/7/2009. k, ADCC activity of mAbs targeting the CS and anchor epitopes. Dashed line represents the limit of detection (L.O.D). l, ADCC potency of mAbs targeting the anchor (n = 18 mAbs) and CS (n = 8 mAbs) epitopes. Data in h, j, and l are represented as mean ± S.D. Data in j and l were analyzed using a two-tailed unpaired non-parametric Mann-Whitney test. Residues in blue are located on HA1 and residues in red are located on HA2. Outlines represent binding footprints of 047-09 4F04 (sky blue) and CR9114 (green). g, Heatmap of mAb binding to WT and mutant HAs shown as the proportion of signal relative to mAb binding to the WT HA. Data in a and d were analyzed using Fisher's Exact tests. Data in e were analyzed using a two-tailed unpaired non-parametric Mann-Whitney test. For a, b, and d, 10 mice per treatment group were used and data are pooled from two independent experiments. For c, 5 mice per treatment group and timepoint were used except for anchor cocktail group at dpi 3 only 4 mice were used. Data in a, b, and d are represented as mean ± S.D and data in c are represented by mean ± S.E.M. Kaplan-Meier curves in a, b, d were analyzed using a Log-rank Mantel-Cox test, and data in c were analyzed using multiple two-tailed unpaired non-parametric Kruskal-Wallis tests.

Extended Data Fig. 4 | See next page for caption. Fig. 4 | Additional repertoire and structural features of mAbs binding the anchor epitope. Related to Fig. 2 . a, VH locus usage by anchor-(n = 52 mAbs) and CS-binding mAbs (n = 37 mAbs). b, VH1 gene usage of mAbs targeting the CS epitope. c, VK locus usage by anchor-(n = 52 mAbs) and CS-binding mAbs (n = 37 mAbs). d, JK gene usage by anchor epitope-binding mAbs. e, Clonal expansions of anchor epitope-targeting mAbs. Numbers indicate heavy and light chain parings, which are described in Extended Data Table 2 . f, Heavy and light chain sequences of the public clone. g, h, Mutations (g) and CDR3 amino acid (AA) lengths (h) of heavy and light chains of mAbs binding the anchor (n = 52 mabs) or CS (n = 37 mAbs) epitopes. Data are mean ± S.D. Fig. 4. a, b , EC50s of serum antibodies competing for binding with 047-09 4F04 for binding to the anchor epitope (a) and CR9114 for binding to the CS epitope (b). a, b, Kinetics of serum antibody responses against the anchor (a) and CS (b) epitopes. Data are mean + S.D. c, d, Proportion of stalk + mAbs per donor (c) or proportion of donors with an isolated anchor mAb (d) upon first exposure to the pH1N1 virus (2009 MIV cohort) relative to donors who have repeatedly been exposed to pH1N1 (2010 TIV and 2014 QIV). Data in c are mean ± S.D. Data in c includes only donors with an isolated anti-stalk mAb, whereas d includes all donors. e, f, Antibody titers (EC 50 ) of serum antibodies collected on day 113 and day 420 against the anchor epitope (e) and the CS epitope (f). Lines connect titers from the same donor and each pair of symbols represents one donor. g, Proportion of anchor epitopebinding mAbs binding to cHA (cH6/1) or mini-HA (n = 50). h, Representative flow cytometry plots of mAbs binding to A/California/7/2009 Cal09 HA and mini-HA (left) and geometric mean fluorescence intensity (gMFI) of mAbs binding to Cal09 and mini-HA (right). Data represent the mean ± S.D. and each symbol represents an individual mAb. i, Proportion of anchor epitopetargeting mAbs binding to A/California/7/2009 recombinant HA with a GCN4 or fibritin trimerization domain (n = 50). For data in a, b, e, and f, Group 1 n = 10 participants, group 2 n = 7 participants, group 4 n = 7 participants, group 3&5 n = 6. For data in c, first exposure n = 7 donors and repeated exposure n = 4 donors. For data in d, first exposure n = 10 donors and repeated exposure n = 13 donors. Data in a, b were analyzed using a two-tailed two-way ANOVA testing for simple effects within rows, data in c and h were analyzed using a two-tailed unpaired non-parametric Mann-Whitney test, data in d, g, and i were analyzed using Fisher's Exact test, and data in c, d were analyzed using a two-tailed paired non-parametric Wilcoxon matched-pairs signed rank test. See also Supplementary Fig. 1 for gating strategy for panel h.

Extended Data Table 2 | Anchor epitope-binding mAb information. Related to Fig. 1 and 2 Extended Data Table 3 | Mutation information for Extended Data Fig. 2f, g Article Extended Data Table 4 | Cryo-EM data collection and refinement statistics for 222-1C06 binding H1. Related to Fig. 2