+
Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Letter
  • Published:

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

Abstract

The conserved region of influenza hemagglutinin (HA) stalk (or stem) has gained attention as a potent target for universal influenza vaccines1,2,3,4,5. Although the HA stalk region is relatively well conserved, the evolutionarily dynamic nature of influenza viruses6 raises concerns about the possible emergence of viruses carrying stalk escape mutation(s) under sufficient immune pressure. Here we show that immune pressure on the HA stalk can lead to expansion of escape mutant viruses in study participants challenged with a 2009 H1N1 pandemic influenza virus inoculum containing an A388V polymorphism in the HA stalk (45% wild type and 55% mutant). High level of stalk antibody titers was associated with the selection of the mutant virus both in humans and in vitro. Although the mutant virus showed slightly decreased replication in mice, it was not observed in cell culture, ferrets or human challenge participants. The A388V mutation conferred resistance to some of the potent HA stalk broadly neutralizing monoclonal antibodies (bNAbs). Co-culture of wild-type and mutant viruses in the presence of either a bNAb or human serum resulted in rapid expansion of the mutant. These data shed light on a potential obstacle for the success of HA-stalk-targeting universal influenza vaccines—viral escape from vaccine-induced stalk immunity.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Association between pre-existing anti-HA stalk immunity and mutant virus selection in humans.
Fig. 2: Effect of A388V stalk mutation on viral replicative fitness.
Fig. 3: A significant conformational change to the HA stalk region induced by A388V mutation.
Fig. 4: Increased resistance to broadly neutralizing antibodies by A388V mutation and its rapid selection by immune pressure.

Similar content being viewed by others

Data availability

The data that support the findings of this study are available from the corresponding author upon reasonable request.

References

  1. Okuno, Y., Isegawa, Y., Sasao, F. & Ueda, S. A common neutralizing epitope conserved between the hemagglutinins of influenza A virus H1 and H2 strains. J. Virol. 67, 2552–2558 (1993).

    CAS  PubMed  PubMed Central  Google Scholar 

  2. Krammer, F., Pica, N., Hai, R., Margine, I. & Palese, P. Chimeric hemagglutinin influenza virus vaccine constructs elicit broadly protective stalk-specific antibodies. J. Virol. 87, 6542–6550 (2013).

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Impagliazzo, A. et al. A stable trimeric influenza hemagglutinin stem as a broadly protective immunogen. Science 349, 1301–1306 (2015).

    CAS  PubMed  Google Scholar 

  4. Yassine, H. M. et al. Hemagglutinin-stem nanoparticles generate heterosubtypic influenza protection. Nat. Med. 21, 1065–1070 (2015).

    CAS  PubMed  Google Scholar 

  5. Andrews, S. F., Graham, B. S., Mascola, J. R. & McDermott, A. B. Is it possible to develop a “universal” influenza virus vaccine? Immunogenetic considerations underlying B-cell biology in the development of a pan-subtype influenza A vaccine targeting the hemagglutinin stem. Cold Spring Harb. Perspect. Biol. 10, a029413 (2018).

    PubMed  PubMed Central  Google Scholar 

  6. Rambaut, A. et al. The genomic and epidemiological dynamics of human influenza A virus. Nature 453, 615–619 (2008).

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Xie, H. et al. H3N2 Mismatch of 2014-15 northern hemisphere influenza vaccines and head-to-head comparison between human and ferret antisera derived antigenic maps. Sci. Rep. 5, 15279 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  8. Bush, R. M., Fitch, W. M., Bender, C. A. & Cox, N. J. Positive selection on the H3 hemagglutinin gene of human influenza virus A. Mol. Biol. Evol. 16, 1457–1465 (1999).

    CAS  PubMed  Google Scholar 

  9. Altman, M. O., Angeletti, D. & Yewdell, J. W. Antibody immunodominance: the key to understanding influenza virus antigenic drift. Viral Immunol. 31, 142–149 (2018).

    CAS  PubMed  PubMed Central  Google Scholar 

  10. Taubenberger, J. K. & Morens, D. M. Influenza: the once and future pandemic. Public Health Rep. 125, 16–26 (2010).

    PubMed  PubMed Central  Google Scholar 

  11. Morens, D. M. & Taubenberger, J. K. Making universal influenza vaccines: lessons from the 1918 pandemic. J. Infect. Dis. 219, S5–S13 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Paules, C. I. & Fauci, A. S. Influenza vaccines: good, but we can do better. J. Infect. Dis. 219, S1–S4 (2019).

    PubMed  PubMed Central  Google Scholar 

  13. Park, J. K., et al. Evaluation of preexisting anti-hemagglutinin stalk antibody as a correlate of protection in a healthy volunteer challenge with influenza A/H1N1pdm Virus. mBio 9, e02284-17 (2018).

  14. Christensen, S. R., et al. Assessing the protective potential of H1N1 influenza virus hemagglutinin head and stalk antibodies in humans. J Virol. 93, e02134-17 (2019).

  15. Ng, S. et al. Novel correlates of protection against pandemic H1N1 influenza A virus infection. Nat. Med. 25, 962–967 (2019).

  16. Kirkpatrick, E., Qiu, X., Wilson, P. C., Bahl, J. & Krammer, F. The influenza virus hemagglutinin head evolves faster than the stalk domain. Sci. Rep. 8, 10432 (2018).

    PubMed  PubMed Central  Google Scholar 

  17. Wu, N. C. & Wilson, I. A. A perspective on the structural and functional constraints for immune evasion: insights from influenza virus. J. Mol. Biol. 429, 2694–2709 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Andrews, S. F. et al. Immune history profoundly affects broadly protective B cell responses to influenza. Sci. Transl. Med. 7, 316ra192 (2015).

    PubMed  PubMed Central  Google Scholar 

  19. Wang, R. et al. Examining the hemagglutinin subtype diversity among wild duck-origin influenza A viruses using ethanol-fixed cloacal swabs and a novel RT–PCR method. Virology 375, 182–189 (2008).

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Throsby, M. et al. Heterosubtypic neutralizing monoclonal antibodies cross-protective against H5N1 and H1N1 recovered from human IgM+ memory B cells. PLoS ONE 3, e3942 (2008).

    PubMed  PubMed Central  Google Scholar 

  21. Clementi, N. et al. A human monoclonal antibody with neutralizing activity against highly divergent influenza subtypes. PLoS ONE 6, e28001 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Tan, G. S. et al. A pan-H1 anti-hemagglutinin monoclonal antibody with potent broad-spectrum efficacy in vivo. J. Virol. 86, 6179–6188 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Henry Dunand, C. J. et al. Preexisting human antibodies neutralize recently emerged H7N9 influenza strains. J. Clin. Invest. 125, 1255–1268 (2015).

    PubMed  PubMed Central  Google Scholar 

  24. Henry Dunand, C. J. et al. Both neutralizing and non-neutralizing human h7n9 influenza vaccine-induced monoclonal antibodies confer protection. Cell Host Microbe 19, 800–813 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Chai, N. et al. Two escape mechanisms of influenza a virus to a broadly neutralizing stalk-binding antibody. PLoS Pathog. 12, e1005702 (2016).

    PubMed  PubMed Central  Google Scholar 

  26. Anderson, C. S. et al. Natural and directed antigenic drift of the H1 influenza virus hemagglutinin stalk domain. Sci. Rep. 7, 14614 (2017).

    PubMed  PubMed Central  Google Scholar 

  27. Memoli, M. J. et al. Validation of the wild-type influenza A human challenge model H1N1pdMIST: an A(H1N1)pdm09 dose-finding investigational new drug study. Clin. Infect. Dis 60, 693–702 (2015).

    CAS  PubMed  Google Scholar 

  28. Memoli, M. J. et al. Evaluation of antihemagglutinin and antineuraminidase antibodies as correlates of protection in an Influenza A/H1N1 virus healthy human challenge model. mBio 7, e00417-16 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Xiao, Y. et al. Deep sequencing of 2009 influenza A/H1N1 virus isolated from volunteer human challenge study participants and natural infections. Virology 534, 96–107 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Poritz, M. A. et al. FilmArray, an automated nested multiplex PCR system for multi-pathogen detection: development and application to respiratory tract infection. PLoS ONE 6, e26047 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Ekiert, D. C. et al. Antibody recognition of a highly conserved influenza virus epitope. Science 324, 246–251 (2009).

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Dreyfus, C. et al. Highly conserved protective epitopes on influenza B viruses. Science 337, 1343–1348 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Corti, D. et al. A neutralizing antibody selected from plasma cells that binds to group 1 and group 2 influenza A hemagglutinins. Science 333, 850–856 (2011).

    CAS  PubMed  Google Scholar 

  34. Wrammert, J. et al. Broadly cross-reactive antibodies dominate the human B cell response against 2009 pandemic H1N1 influenza virus infection. J. Exp. Med. 208, 181–193 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Wu, Y. et al. A potent broad-spectrum protective human monoclonal antibody crosslinking two haemagglutinin monomers of influenza A virus. Nat. Commun. 6, 7708 (2015).

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Huang, K. Y. et al. Focused antibody response to influenza linked to antigenic drift. J. Clin. Invest. 125, 2631–2645 (2015).

    PubMed  PubMed Central  Google Scholar 

  37. Doud, M. B., Lee, J. M. & Bloom, J. D. How single mutations affect viral escape from broad and narrow antibodies to H1 influenza hemagglutinin. Nat. Commun. 9, 1386 (2018).

    PubMed  PubMed Central  Google Scholar 

  38. Lee, J. et al. Molecular-level analysis of the serum antibody repertoire in young adults before and after seasonal influenza vaccination. Nat. Med. 22, 1456–1464 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Linderman, S. L. & Hensley, S. E. Antibodies with ‘Original Antigenic Sin’ properties are valuable components of secondary immune responses to influenza viruses. PLoS Pathog. 12, e1005806 (2016).

    PubMed  PubMed Central  Google Scholar 

  40. Coughlan, L. & Palese, P. Overcoming barriers in the path to a universal influenza virus vaccine. Cell Host Microbe 24, 18–24 (2018).

    CAS  PubMed  Google Scholar 

  41. Neumann, G. et al. Generation of influenza A viruses entirely from cloned cDNAs. Proc. Natl Acad. Sci. USA 96, 9345–9350 (1999).

    CAS  PubMed  Google Scholar 

  42. Reed, L. J. & H., M. A simple method of estimating fifty percent endpoints. Am. J. Epidemiol. 27, 493–497 (1938).

    Google Scholar 

  43. Krammer, F. et al. A carboxy-terminal trimerization domain stabilizes conformational epitopes on the stalk domain of soluble recombinant hemagglutinin substrates. PLoS ONE 7, e43603 (2012).

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Stevens, J. et al. Structure of the uncleaved human H1 hemagglutinin from the extinct 1918 influenza virus. Science 303, 1866–1870 (2004).

    CAS  PubMed  Google Scholar 

  45. Pettersen, E. F. et al. UCSF Chimera–a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612 (2004).

    CAS  PubMed  Google Scholar 

  46. Xu, R. et al. Structural basis of preexisting immunity to the 2009 H1N1 pandemic influenza virus. Science. 328, 357–360 (2010).

    CAS  PubMed  PubMed Central  Google Scholar 

  47. WHO Global Influenza Surveillance Network. Manual for the Laboratory Diagnosis and Virological Surveillance of Influenza. (World Health Organization, 2011).

  48. Shu, B. et al. Design and performance of the CDC real-time reverse transcriptase PCR swine flu panel for detection of 2009 A (H1N1) pandemic influenza virus. J. Clin. Microbiol. 49, 2614–2619 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was supported by the intramural research program of the National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH). We thank the Comparative Medicine Branch (NIAID, NIH) for assistance with the animal studies and S. Hunsberger (Biostatistics Research Branch, NIAID, NIH) for the helpful discussion of statistical analyses performed in this study. We also thank Janssen Pharmaceutica, R. Ahmed (Emory Vaccine Center) and A. Townsend (University of Oxford) for generously sharing monoclonal antibodies.

Author information

Authors and Affiliations

Authors

Contributions

J.-K.P., M.J.M., J.C.K. and J.K.T. conceived and designed the study. J.-K.P., Y.X., M.D.R., L.A.R., S.F., A.M.M., A.D.F., M.A.G., N.A.B., L.Q., X.Y. and K.S. generated the laboratory data. S.R., R.A., L.C., A.H. and M.J.M. designed and performed the primary clinical study. J.-K.P., A.H., Y.X., M.D.R., D.M.M., K.-A.W., M.J.M., J.C.K. and J.K.T. interpreted the data. J.-K.P., D.M.M., K.-A.W., M.J.M., J.C.K. and J.K.T. wrote the manuscript. All authors critically reviewed the paper and approved of the final version of the paper for submission.

Corresponding author

Correspondence to Jeffery K. Taubenberger.

Ethics declarations

Competing interests

All authors declare no competing financial or non-financial interests.

Additional information

Peer review information Alison Farrell was the primary editor on this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data

Extended Data Fig. 1 A Single Nucleotide Polymorphism (SNP) assay using minor groove binder (MGB)-based TaqMan probes for detecting wild-type (A388) or mutant (V388) genes in various conditions.

A SNP assay was developed for detecting wild-type (A388) or mutant (V388) HA genes in various conditions utilizing a set of Minor Groove Binder (MGB)-based TaqMan probes: VIC-labeled probe detecting the wild-type (A388) and a FAM-labeled probe detecting the mutant (V388). The SNP assay was validated using a mixed viral genome from the 2009 H1N1pdm wild-type (A388) and the mutant (V388) viruses. Viral RNA was mixed in varying ratios from 10:0 (0% mutant virus) to 0:10 (100% mutant virus). For the validation of two-step SNP assay for analyzing nasal wash samples from the human challenge study, mixed viral RNA was diluted to represent varying viral loads of (a) 105.5 TCID50/ml, (b) 104.0 TCID50/ml, and (c) 103.0 TCID50/ml. Prepared viral RNA was analyzed by the two-step SNP assay (see Methods). Blue and red bars indicate the ∆Rn value of wild-type and mutant virus, respectively. Graphs show mean ∆Rn value and standard deviation from 4 independent experiments (a-c, n = 4). For the validation of one-step SNP assay used to analyze selection dynamics in vitro, mixed viral RNA was diluted to represent varying viral loads of (d) 107.0 TCID50/ml, (e) 105.5 TCID50/ml, and (f) 104.0 TCID50/ml. Prepared viral RNA was analyzed by the one-step SNP assay (see Methods). Blue bars indicate the threshold cycle (Ct) values representing the amount of wild-type (A388) virus. Red bars represent the amount of mutant (V388) virus. Graphs show mean Ct value and standard deviation from 3 independent experiments (d-f, n = 3). Results show that the SNP assay reliably detects minor population, either wild-type or mutant, existing as low as 10% across the various viral loads.

Extended Data Fig. 2 Comparison of different parameters between study participants with different selection outcomes.

Different parameters were compared between groups to find possible correlates of the selection outcomes. a, Pre-challenge serum hemagglutination inhibition (HAI) titers were compared. Black horizontal lines and grey error bars represent geometric mean titer and 95% CI, respectively. Dashed line shows the detection limit. b, The fraction of HA stalk antibodies relative to total HA antibodies was compared. The ratio is a calculation of serum anti-stalk IgG and anti-full-length HA IgG titers measured by ELISA. c, Serum IgG titer data (anti-stalk and anti-total HA) generated in (b) were used to analyze the correlation between the stalk and total HA antibody titers in each selection group. A two-tailed Pearson’s correlation coefficient (Pearson’s r) was used for the analysis. The positive correlation between the stalk antibody titers and total HA titers may explain the lack of difference in stalk/total HA antibody ratio between groups seen in (b). To find a possible role of immunological imprinting on the selection outcomes, (d) age of the participants at the study enrollment and (e) birth year was compared. To find potential loss in viral fitness in humans caused by the A388V mutation, (f) days of symptoms, (g) number of symptoms, and (h) duration of shedding after the challenge infection were compared between groups. Black horizontal lines and grey error bars represent median value and 95% CI, respectively (b,d,e-h). A two tailed nonparametric Mann-Whitney test was used to compare parameters from the mutant selection group to other selection groups. Each symbol represents an individual study participant and their selection outcome (wild-type selection, n = 8; mixed shedding, n = 9; mutant selection, n = 12).

Extended Data Fig. 3 A significant conformational change to the HA stalk region induced by A388V mutation.

MDCK cells were infected with 1 multiplicity of infection (MOI) of wild-type (A388) or mutant (V388) H1N1pdm viruses generated by reverse genetics. 24 hours after infection, cells were harvested, and the expression of wild-type and mutant HA was measured by flow cytometry to evaluate the effect of the A388V mutation on the HA stalk epitopes. a, A representative example of the gating strategy. Dead cells and debris were excluded based on FSC/SSC cell dot plot. Anti-influenza nucleoprotein (NP) antibodies conjugated with allophycocyanin (APC) were used for gating. Only infected cells, expressing NP, were used for the analysis. Broadly neutralizing antibodies binding to the HA stalk, (b) CR6261, (c) CR9114, (d) FI6V3, (e) 70-1F02, (f) C179, and (g) CT149 were conjugated with fluorescein isothiocyanate (FITC). A monoclonal antibody that binds to the HA globular head, (h) EM-4C04, was conjugated with r-phycoerythrin (R-PE). Each stalk-binding antibody was mixed with the head-binding EM-4C04 antibody and NP antibody, and the antibody mixtures were used to stain cells expressing the wild-type or mutant HA. Histograms are colored differently to show different experimental groups: Blue - cells infected with wild-type (A388) virus; Pink - cells infected with mutant (V388) virus; Black - unstained cell control. Representative histograms from three independent experiments are shown. The summary table shows the average median fluorescence intensity (MFI) and standard error of mean (s.e.m.) of the three independent experiments.

Extended Data Fig. 4 Comparison between stalk-only and full-length HA construct.

Stalk-only constructs with or without the A388V mutation (See Extended Data Fig. 7b for the amino acid sequence) were produced to measure the level of antibodies recognizing the mutant stalk in human serum while excluding head-binding antibodies (Fig. 3j). To confirm that the stalk-only construct with A388V mutation appropriately represents the natural A388V stalk structure, the level of decrease in broadly neutralizing monoclonal antibodies (bNAbs) binding was compared between the stalk-only construct and the full-length HA. ELISA was performed using serially diluted (a,b) CR6261, (c,d) CR9114, (e,f) FI6V3, (g,h) 70-1F02, (I,j) C179, and (k,l) CT149. (m,n) Anti-StrepTag II antibody was used to show that equal amounts of wild-type and mutant antigen were used for the analysis. The AUC was calculated using GraphPad Prism8 (v.8.3.0). The AUC for the full-length HAs (b,d,f,h,j,l) was calculated using the data from Fig. 3. Graphs show mean and standard deviations from three independent measurements. The summary table shows the mean OD492 values and standard error of mean (s.e.m.) of the three independent measurements. The comparison result (summarized in the table) shows that the A388V stalk-only construct closely represents the natural stalk structure of the full-length A388V HA.

Extended Data Fig. 5 AUC calculation using wile-type and mutant stalk-only construct.

Raw ELISA data used to generate Fig. 3j are shown. Stalk-only constructs with or without the A388V mutation (See Extended Data Fig. 7b for the amino acid sequence) were used to measure changes in the level of antibodies recognizing the wild-type (A388) of mutant (V388) stalk in human serum. Twenty-nine pre-challenge serum samples from the influenza human challenge study participants were serially diluted and used for the ELISA. The AUC was calculated using GraphPad Prism8 (v.8.3.0) with the baseline value of 0.1 (approximately 2 times the OD492 value from the control wells). Blue lines and numbers show the ELISA data and the AUC, respectively, obtained using the wild-type (A388) stalk construct. Red lines and numbers show the ELISA data and the AUC, respectively, obtained using the mutant (V388) stalk construct. Graphs show mean and standard deviations from three independent measurements.

Extended Data Fig. 6 Association between the stalk antibody titers and selection pressure measured in vitro.

Selection pressure placed by individual human sera was measured using pre-challenge serum from the influenza human challenge study participants. a, An equal mixture of the wild-type (A388) and mutant (V388) virus was cultured with 1:50 diluted serum samples. Sera from Q1, Q2, Q3 and Q4 from Fig. 1 (n = 7 per quartile) were used. Culture supernatants were collected at 48 and 72 hours after infection followed by viral RNA extraction. The selection dynamics were measured using a Single Nucleotide Polymorphism (SNP) assay utilizing a set of Minor Groove Binder (MGB)-based TaqMan probes; VIC-labeled probe detects the wild-type (A388); FAM-labeled probe detects the mutant (V388). Data are presented as the threshold cycle (Ct) value from the SNP assay. Dashed lines show the Ct value limit (Ct 40) of the SNP assay. A Ct value of 41 was given to undetected signals to generate graphs. Error bars represent standard deviations from three independent experiments. The final dilution is noted on the individual graph if higher than 1:50. b, Mutant selection index was calculated based on data from (a) by ∆∆Ct method using controls cultured without serum (see Methods). A mutant selection index higher than 0 (pink area) indicates a serum sample selected for the mutant virus. An index lower than 0 (blue area) indicates a serum sample selected for the wild-type virus. Horizontal lines show median values and error bars represent 95% CI. The indexes between samples from different quartiles were compared using nonparametric one-way analysis of variance (Kruskal-Wallis test) and Dunn’s test as a post-hoc test. c, Correlation between the anti-stalk serum IgG titer and the selection index of 29 sera samples were analyzed by calculating two-tailed Spearman’s rank correlation coefficient (Spearman r). The best-fit line was plotted using simple linear regression analysis. Statistical analyses were performed using GraphPad Prism8 (v.8.3.0).

Extended Data Fig. 7 Sequence of wild-type and mutant full-length HA or HA stalk-only constructs.

Amino acid sequences of (a) full-length HA and (b) stalk-only proteins with or without A388V mutation are shown. A388V mutation, HA trimerization domain, and StrepTag II sequence are highlighted. Consensus amino acid sequences are shown in dots.

Supplementary information

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Park, JK., Xiao, Y., Ramuta, M.D. et al. Pre-existing immunity to influenza virus hemagglutinin stalk might drive selection for antibody-escape mutant viruses in a human challenge model. Nat Med 26, 1240–1246 (2020). https://doi.org/10.1038/s41591-020-0937-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • Version of record:

  • Issue date:

  • DOI: https://doi.org/10.1038/s41591-020-0937-x

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing
点击 这是indexloc提供的php浏览器服务,不要输入任何密码和下载