REVIEW 4 major objections 6 minor 35 references
Could short-lasting non-specific immunity explain seasonal spacing of epidemic diseases?
T0 review · 4 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read Short-lived non-specific immunity can make epidemics alternate without any seasonal forcing.
desk verdict The core dynamical claim – that short-lived non-specific immunity can make SIRS systems oscillate and alternate epidemics – is real and worth referee time, but the PIV/seasonality part is a parameter-tuned illustration, not a test. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The key machinery is an extended SIRS model with a host state for general immunity, G: a host is either susceptible, infectious, or specifically immune for each disease, and additionally either protected or unprotected by a non-specific immunity that lasts on average TG time units. The non-specific immunity is switched on by any contact with any disease—even when specific immunity prevents infection—and blocks infection by all diseases while active. This state, inserted between exposure and infection, is what turns the model's endemic equilibrium into a Hopf bifurcation and sustains periodic epidemics; the duration TG acts as the interaction-strength control parameter.
What would settle it
A direct test would compare two groups of animals exposed to pathogen B either after a contact that did not produce a full infection with pathogen A or after a full A infection; if contact alone does not shorten or block B infection, the model's strongest-interaction regime is not realized. A second check is measuring the duration TG required for alternating epidemics in a two-pathogen experiment and comparing it with the model's predicted bifurcation threshold.
Extended reading notes
Core claim
The paper's central discovery is that adding a short non-specific immune state to the classical SIRS model destabilizes the steady endemic equilibrium and produces stable epidemic limit cycles. In the two-disease case, disease A's outbreak temporarily shields the host population from disease B; while B is suppressed, the pool of hosts susceptible to B grows, so when the non-specific protection fades B erupts in a large epidemic, and the roles reverse. This gives alternating recurrent epidemics for realistic parameters (long specific immunity, R0 around 2.5, non-specific immunity of a few weeks). Under a single seasonal modulation of transmissibility, the interaction spreads the epidemic peaks of otherwise identical diseases across the year, and a tuned two-disease version matches several qualitative features of US parainfluenza PIV-1/PIV-3 time series, including deeper PIV-3 valleys and delayed PIV-3 peaks in PIV-1 years.
Load-bearing premise
The load-bearing assumption is that every contact with a pathogen—even one that cannot cause infection because the host already has specific immunity—activates the short non-specific protection; if activation requires a productive infection, the disease interaction is weaker and the alternating-epidemic regime may shrink or disappear.
Editorial extensions
If this is right
- In a population with two similar respiratory pathogens, the model predicts spontaneous anti-phase epidemic cycles for a wide range of realistic transmission and immunity parameters.
- A single seasonal driver is sufficient to generate sequentially spaced epidemics of diseases that have identical dynamics; distinct environmental drivers for each virus are not required.
- The model reproduces the distinctive PIV-1/PIV-3 pattern: annual PIV-3 peaks, biennial PIV-1 peaks, deeper PIV-3 valleys, and delayed PIV-3 peaks in PIV-1 years.
- With more than a few diseases, sustained alternating epidemics require narrow epidemic spikes, so larger transmissibility and longer specific immunity increase the number of diseases that can be spaced in time.
- Without non-specific immunity, the same two-disease system with the same seasonal driver would produce annual outbreaks of both diseases, so the interaction is what explains the spacing.
Reading between the lines
- If the paper is right, observed seasonal 'niches' of respiratory viruses may be partly an emergent property of pathogen competition, and vaccination campaigns against one respiratory virus could temporarily suppress others via the same non-specific window.
- The mechanism implies that short cross-protective windows should be detectable in routine surveillance as negative correlations between unrelated respiratory virus activities at lags of weeks to months.
- An extension worth testing is whether adding a second, longer non-specific immune component can reproduce multi-year patterns such as RSV's biennial cycle and influenza's subtype alternations.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes an extension of the classical SIRS model to multiple diseases that interact through a short-lasting, non-specific immunity (G) state. The authors show, via direct simulation and local stability analysis, that sufficiently long non-specific immunity (TG) can produce recurrent alternating epidemic cycles even without seasonal forcing, and that with a single seasonal driver the model can produce sequential epidemic patterns with different phases. They apply the model to US parainfluenza data (PIV-1 and PIV-3), claiming to reproduce several qualitative features including biennial PIV-1 peaks, deeper PIV-3 valleys in PIV-1 years, delayed PIV-3 peaks, and asymmetric PIV-1 epidemic shape. The central claim is that a few-weeks-long non-specific immune response is sufficient to explain seasonal spacing of unrelated respiratory viruses.
Significance. If the central claim is correct, the paper offers a parsimonious mechanism for seasonally spaced epidemics without invoking a separate environmental driver for each virus, and it provides a concrete, testable model of disease interference via short-lived innate or club-cell-mediated immunity. The paper's strength is that the alternating-epidemic phenomenon is derived rather than fitted: the Hopf-bifurcation analysis and deterministic/stochastic simulations agree, and the model is simple enough to be reproduced from the text. The PIV-1/PIV-3 comparison is a genuine attempt to confront the model with data. However, the significance is qualified by the model's strong activation assumption and by the fact that the PIV reproduction involves substantial parameter tuning without uncertainty quantification.
major comments (4)
- [Sec. II, Multiple diseases] The model assumes that non-specific immunity is activated on every contact with any disease, even when host-specific immunity prevents infection, as stated in Sec. II: 'The non-specific immunity is activated (Goff→Gon) every time a host comes into contact with any of the diseases, also if infection is prevented by specific immunity.' This maximizes inter-disease coupling and is the key nonlinearity that generates the Hopf bifurcation behind alternating epidemics (Sec. III.1, Figs. 2–3). The only robustness check reported, 'twice as long, but not fully protective (50%)' with data not shown, varies the duration and strength of G but not the triggering condition. If activation requires a productive infection (S→I), hosts in R for one disease would not enter G upon contact with that disease, reducing the protected pool and weakening the postponement mechanism that builds large susceptible pools. The paper should test this alternative triggering rule and show whether the alternating-epidemic regime persists for biologically plausible parameters.
- [Sec. II, Model and Methods; Sec. III.1] The deterministic rate equations are never written out. The text describes the model verbally and gives parameter values, and Sec. III.1 states that a linear stability analysis was performed on a grid, but the equations for the multi-disease SIRS model with the G state are not displayed. Without the explicit equations, the Hopf-bifurcation calculation cannot be independently checked, and the reader cannot verify how the external reservoir, the seasonal forcing R0(t), and the symmetry-breaking initial conditions enter the deterministic system. Please provide the full set of ordinary differential equations, including the G-transition terms, and specify the numerical method used for the bifurcation analysis.
- [Sec. III.5, Fig. 6] The PIV-1/PIV-3 'reproduction' is a parameter-tuned qualitative match rather than an estimated fit. The parameters TR,PIV-1, TR,PIV-3, TG, the seasonal amplitude c, the phase offset (peak aligned to March 1), and the external reservoir size are all chosen to reproduce the observed pattern, yet no uncertainty quantification, sensitivity analysis, or systematic parameter search is provided. In particular, the statement that the same model 'without non-specific immunity' converges to annual epidemics for both diseases is a useful control, but it does not establish that the tuned parameter set is biologically plausible or that the fit is unique. To support the claim that the model 'reproduces multiple features,' the authors should show how sensitive the four listed features (biennial PIV-1, deep PIV-3 valleys, delayed PIV-3 peaks, asymmetric PIV-1 shape) are to variations in each tuned parameter, and ideally compare the model output to data quantitatively (e.g., peak timing and valley depth).
- [Sec. II, Simulation details; Fig. 2c] The stability diagram in Fig. 2c and the matching statement in Sec. III.1 ('The amplitudes measured in deterministic simulations match well with the calculated bifurcation lines') are not fully documented. The plotted quantity is log10(Imax/Imin), which is a spikiness measure rather than a true amplitude, and the white dashed line is described only as a bifurcation line. It would be helpful to show explicitly how the bifurcation line was computed from the linearized equations and to report the agreement between the predicted oscillation period/amplitude and the simulated values. Without this, the reader cannot separate the contribution of the Hopf calculation from the simulation-based color map.
minor comments (6)
- [Abstract] The abstract contains a typo: 'cloub cell' should be 'club cell'.
- [Fig. 6 caption] The caption states 'TR,P IV 1 = 170TI and TR,P IV 1 = 120TI'; the second should presumably be TR,P IV 3 = 120TI. Please correct the label.
- [Fig. 3 caption] The caption uses '∆T0' where the text and equations define ∆TR; this is inconsistent and should be fixed.
- [Sec. II, Model parameters] The notation 'R0(t) = 2 + c·sin(2πt/τ)' is clear, but the relationship between the base R0=2.5 used in the no-seasonality cases and the seasonal mean of 2 in Sec. III.4 and III.5 should be stated explicitly to avoid confusion.
- [Sec. IV, Discussion] The statement 'This seems to indicate, that the observed effect would not occur just because of behavioural changes caused by the sickness, such as staying home from work' is an interesting inference, but the model does not explicitly include behavior; please rephrase to avoid overclaiming.
- [Data, code and materials] The paper says the model is 'completely reproducible from the description in the main text,' but no code or data repository is given. Since the simulations and the extracted PIV time series are nontrivial, providing code and the digitized data would materially aid reproducibility.
Circularity Check
No significant circularity: the recurrent-epidemic result is an emergent model property, and the PIV comparison is explicitly a parameter-tuned reproduction rather than an out-of-sample prediction.
full rationale
The paper's central derivation chain is self-contained. The model is a multi-disease SIRS extension with a non-specific immunity compartment G; recurrent alternating epidemics are obtained from agent-based and deterministic simulations and from a linear-stability/Hopf-bifurcation analysis (Sec. III.1, Figs. 2-3), not by inserting the target pattern as an assumption. The seasonal-pattern results (Sec. III.4, Fig. 5) are generated by varying TG under fixed seasonal forcing and are not fit to specific empirical time series. The only data comparison, the PIV-1/PIV-3 application (Sec. III.5, Fig. 6), is described as reproducing observed features, not as predicting them from fixed constants; the caption states the chosen parameters (TG=4TI, TR,PIV1=170TI, TR,PIV3=120TI, R0(t)=2+0.2 sin, phase aligned to March 1, external reservoir 1/200,000). This is a post-hoc sufficiency demonstration with tuned parameters, which limits the strength of the empirical evidence but does not make the claim circular in the defined sense. There are no load-bearing self-citations: the paper cites prior related work by Rohani et al. [11,14,16] and explicitly acknowledges that similar interference phenomena were observed in those models, while the present contribution adds the SIRS/waning-immunity framework and bifurcation analysis. No uniqueness theorem is imported from the authors' own prior work, and no fitted parameter is renamed as a prediction. The core theoretical claim therefore stands as an emergent derivation from the model assumptions rather than a reduction to its inputs.
Assumptions & free parameters
free parameters (7)
- Basic reproduction number R0 (base) =
2.5 (nonseasonal case studies)
- R0 seasonal amplitude c =
0.1 or 0.2 in different figures
- Specific immunity duration TR =
100 TI (generic); 170 TI for PIV-1 and 120 TI for PIV-3
- Non-specific immunity duration TG =
4 TI in PIV case; scanned 0.1-15 TI elsewhere
- Seasonal phase offset =
Annual R0 maxima aligned to March 1
- External reservoir concentration =
1/200,000 of population (seasonal simulations); 1/500,000 otherwise
- Initial condition symmetry breaking =
1/1000 infectious with disease 1, 2/1000 with disease 2
assumptions (7)
- domain assumption Homogeneous, fully mixed host population with mass-action transmission (standard SIRS).
- ad hoc to paper Non-specific immunity is activated on every contact with any disease, even when infection is blocked by specific immunity.
- ad hoc to paper While active, non-specific immunity is 100% effective at blocking infection.
- domain assumption Specific immunity wanes with a mean duration TR, allowing reinfection (SIRS).
- domain assumption An external reservoir of constantly infectious hosts prevents stochastic extinction.
- domain assumption Seasonal forcing is a single sinusoidal modulation of R0 with period one year.
- domain assumption The extracted PIV-1/PIV-3 surveillance data from Fry et al. is sufficiently accurate for qualitative comparison.
Cite this review
Pith. "Pith review of Could short-lasting non-specific immunity explain seasonal spacing of epidemic diseases?." pith.science (2026). https://pith.science/paper/G45RJ6AB
@misc{pith2026190800925,
author = {Pith},
title = {Pith review of: Could short-lasting non-specific immunity explain seasonal spacing of epidemic diseases?},
year = {2026},
howpublished = {\url{https://pith.science/paper/G45RJ6AB}},
note = {Machine review of arXiv:1908.00925}
}
read the original abstract
Common respiratory viruses cause seasonal epidemics in sequential patterns. The underlying mechanisms for this pattern have been debated for some time. For influenza, contenders include temperature, humidity, and vitamin D levels. While such seasonal drivers may be sufficient to explain midwinter peaking it is unclear if other respiratory viruses peaking in spring, summer, or fall have individual environmental drivers. Here we present a dynamic model of interacting diseases, in an effort to explain observed seasonal patterns without requiring multiple seasonal drivers. Our model extends the classical SIRS-model to include multiple diseases that interact via a short-lasting, non-specific immunity component. This temporal protection is triggered by exposure to any epidemic disease and lasts only a couple of weeks. We show that the inhibiting disease-interaction allows recurrent epidemic behaviour without any seasonal driving. In the presence of a single seasonal driver our model parsimoniously predicts epidemic patterns such as sequentially occurring epidemic diseases. We also present a two-disease simulation reproducing multiple features observed in time series of parainfluenza strains PIV-3 and biennial PIV-1 epidemic patterns. Complex epidemic patterns of seasonal diseases may be explained by non-specific innate or cloub cell (T-cell) mediated immune responses that result in interaction between unrelated respiratory epidemic diseases.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
A new concept of the epidemic process of influenza a virus
RE Hope-Simpson and DB Golubev. A new concept of the epidemic process of influenza a virus. Epidemiology and infection, 99(01):5–54, 1987
work page 1987
-
[2]
Melatonin me- diates seasonal adjustments in immune function
Randy J Nelson and Deborah L Drazen. Melatonin me- diates seasonal adjustments in immune function. Repro- duction Nutrition Development , 39(3):383–398, 1999
work page 1999
-
[3]
Update: Influenza activity-united states
S Smith, L Blanton, K Kniss, D Mustaquim, C Steffens, C Reed, A Bramley, B Flannery, AM Fry, LA Grohskopf, et al. Update: Influenza activity-united states. MMWR. Morbidity and mortality weekly report , 64(48):1342, 2015
work page 2015
-
[4]
Seasonal trends of human parainfluenza viral infections: United states, 1990–2004
Alicia M Fry, Aaron T Curns, Kathryn Harbour, Lori Hutwagner, Robert C Holman, and Larry J Anderson. Seasonal trends of human parainfluenza viral infections: United states, 1990–2004. Clinical Infectious Diseases , 43(8):1016–1022, 2006
work page 1990
-
[5]
Henri Laurichesse, Daniel Dedman, John M Watson, and Maria C Zambon. Epidemiological features of parain- fluenza virus infections: laboratory surveillance in eng- land and wales, 1975–1997. European journal of epidemi- ology, 15(5):475–484, 1999
work page 1975
-
[6]
R0,TI andTR can be disease specific parameters
See individual figure for specific values. R0,TI andTR can be disease specific parameters. Unless stated otherwise, they are chosen to be identical for all diseases. The mean duration of specific immunity is rather short in most simulations (few years). This reflects the fact that immunity is often incomplete in viruses such as PIV[30] and RSV[31], which allow...
work page 1990
-
[7]
Molly E Counihan, David K Shay, Robert C Holman, Sara A Lowther, and Larry J Anderson. Human parain- fluenza virus-associated hospitalizations among children less than five years of age in the united states. The Pe- diatric infectious disease journal , 20(7):646–653, 2001
work page 2001
-
[8]
M Maes, W Stevens, S Scharpe, E Bosmans, F De Meyer, P D’Hondt, D Peeters, P Thompson, P Cosyns, L De Clerck, et al. Seasonal variation in peripheral blood leukocyte subsets and in serum interleukin-6, and solu- ble interleukin-2 and-6 receptor concentrations in normal volunteers. Experientia, 50(9):821–829, 1994
work page 1994
Show all 35 references
-
[9]
Seasonal dif- ferences in the rhythmicity of human male and female lymphocyte blastogenic responses
FN Boctor, RA Charmy, and EL Cooper. Seasonal dif- ferences in the rhythmicity of human male and female lymphocyte blastogenic responses. Immunological inves- tigations, 18(6):775–784, 1989
1989
-
[10]
Seasonal variation in host susceptibil- ity and cycles of certain infectious diseases
Scott F Dowell. Seasonal variation in host susceptibil- ity and cycles of certain infectious diseases. Emerging infectious diseases, 7(3):369, 2001
2001
-
[11]
Epidemic influenza and vitamin d
JJ Cannell, R Vieth, JC Umhau, MF Holick, WB Grant, S Madronich, CF Garland, and E Giovannucci. Epidemic influenza and vitamin d. Epidemiology and infection, 134 (06):1129–1140, 2006
2006
-
[12]
Population dynamic interference among childhood diseases
Pejman Rohani, David J Earn, B Finkenst¨ adt, and Bryan T Grenfell. Population dynamic interference among childhood diseases. Proceedings of the Royal So- ciety of London B: Biological Sciences , 265(1410):2033– 2041, 1998
1998
-
[13]
Cross- immunity between strains explains the dynamical pattern of paramyxoviruses
Samit Bhattacharyya, Per H Gesteland, Kent Korgen- ski, Ottar N Bjørnstad, and Frederick R Adler. Cross- immunity between strains explains the dynamical pattern of paramyxoviruses. Proceedings of the National Academy of Sciences, 112(43):13396–13400, 2015
2015
-
[14]
The effect of cross- immunity and seasonal forcing in a multi-strain epidemic model
Masashi Kamo and Akira Sasaki. The effect of cross- immunity and seasonal forcing in a multi-strain epidemic model. Physica D: Nonlinear Phenomena , 165(3):228– 241, 2002
2002
-
[15]
The dynamical im- plications of disease interference: correlations and co- existence
Yunxin Huang and Pejman Rohani. The dynamical im- plications of disease interference: correlations and co- existence. Theoretical population biology, 68(3):205–215, 2005
2005
-
[16]
Tracking the dynamics of pathogen interactions: model- ing ecological and immune-mediated processes in a two- pathogen single-host system
Daniel A Vasco, Helen J Wearing, and Pejman Rohani. Tracking the dynamics of pathogen interactions: model- ing ecological and immune-mediated processes in a two- pathogen single-host system. Journal of theoretical biol- ogy, 245(1):9–25, 2007. 8
2007
-
[17]
Ecological interference between fatal diseases
P Rohani, CJ Green, NB Mantilla-Beniers, and BT Gren- fell. Ecological interference between fatal diseases. Na- ture, 422(6934):885–888, 2003
2003
-
[18]
Dynamics of interacting diseases
Joaqu´ ın Sanz, Cheng-Yi Xia, Sandro Meloni, and Yamir Moreno. Dynamics of interacting diseases. Physical Re- view X , 4(4):041005, 2014
2014
-
[19]
Graham, Fabrice Carrat, Oliver Ratmann, and Bernard Cazelles
Anton Camacho, S´ ebastien Ballesteros, Andrea L. Graham, Fabrice Carrat, Oliver Ratmann, and Bernard Cazelles. Explaining rapid reinfections in multiple-wave influenza outbreaks: Tristan da cunha 1971 epidemic as a case study. Proceedings of the Royal Society of London B: Biol...
1971
-
[20]
Mod- eling within-host dynamics of influenza virus infection including immune responses
Kasia A Pawelek, Giao T Huynh, Michelle Quinlivan, Ann Cullinane, Libin Rong, and Alan S Perelson. Mod- eling within-host dynamics of influenza virus infection including immune responses. PLoS Comput Biol , 8(6): e1002588, 2012
2012
-
[21]
Dynamics of influenza virus infection and pathology
Roberto A Saenz, Michelle Quinlivan, Debra Elton, Shona MacRae, Anthony S Blunden, Jennifer A Mum- ford, Janet M Daly, Paul Digard, Ann Cullinane, Bryan T Grenfell, et al. Dynamics of influenza virus infection and pathology. Journal of virology , 84(8):3974– 3983, 2010
2010
-
[22]
Innate immunity and the inter-exposure interval de- termine the dynamics of secondary influenza virus infec- tion and explain observed viral hierarchies
Pengxing Cao, Ada WC Yan, Jane M Heffernan, Stephen Petrie, Robert G Moss, Louise A Carolan, Teagan A Guarnaccia, Anne Kelso, Ian G Barr, Jodie McVernon, et al. Innate immunity and the inter-exposure interval de- termine the dynamics of secondary influenza virus infec- tion and ...
2015
-
[23]
A model of non- specific immunity
Rustom Antia and Jacob C Koella. A model of non- specific immunity. Journal of theoretical biology , 168(2): 141–150, 1994
1994
-
[24]
The unspecific side of acquired immu- nity against infectious disease: Causes and consequences
Eric Muraille. The unspecific side of acquired immu- nity against infectious disease: Causes and consequences. Frontiers in Microbiology, 6, 2015
2015
-
[25]
Ecological and immunological determinants of influenza evolution
Neil M Ferguson, Alison P Galvani, and Robin M Bush. Ecological and immunological determinants of influenza evolution. Nature, 422(6930):428–433, 2003
2003
-
[26]
Training innate immunity: the changing concept of immunological memory in innate host defence
Mihai G Netea. Training innate immunity: the changing concept of immunological memory in innate host defence. European journal of clinical investigation , 43(8):881–884, 2013
2013
-
[27]
Trained immunity: a memory for innate host de- fense
Mihai G Netea, Jessica Quintin, and Jos WM van der Meer. Trained immunity: a memory for innate host de- fense. Cell host & microbe , 9(5):355–361, 2011
2011
-
[28]
Viral cross-reactivity and antigenic determinants recognized by human parain- fluenza virus type 1-specific cytotoxic t-cells
Vibhuti P Dave, Jane E Allan, Karen S Slobod, F Suzette Smith, Kevin W Ryan, Toru Takimoto, Ultan F Power, Allen Portner, and Julia L Hurwitz. Viral cross-reactivity and antigenic determinants recognized by human parain- fluenza virus type 1-specific cytotoxic t-cells. Virology,...
1994
-
[29]
Club cells surviving influenza a virus infection induce tempo- rary nonspecific antiviral immunity
Jennifer R Hamilton, David Sachs, Jean K Lim, Ryan A Langlois, Peter Palese, and Nicholas S Heaton. Club cells surviving influenza a virus infection induce tempo- rary nonspecific antiviral immunity. Proceedings of the National Academy of Sciences , 113(14):3861–3866, 2016
2016
-
[30]
Myxoviruses: Parainfluenza 1
RM Chanock, RH Parrott, KM Johnson, AZ Kapikian, and JA Am Bell. Myxoviruses: Parainfluenza 1. Amer- ican Review of Respiratory Disease , 88(3P2):152–166, 1963
1963
-
[31]
Hemagglutinin sequence clusters and the antigenic evo- lution of influenza a virus
Joshua B Plotkin, Jonathan Dushoff, and Simon A Levin. Hemagglutinin sequence clusters and the antigenic evo- lution of influenza a virus. Proceedings of the National Academy of Sciences , 99(9):6263–6268, 2002
2002
-
[32]
In- fectious diseases of humans: dynamics and control , vol- ume 28, chapter 2
Roy M Anderson, Robert M May, and B Anderson. In- fectious diseases of humans: dynamics and control , vol- ume 28, chapter 2. Wiley Online Library, 1992
1992
-
[33]
Immunity to and frequency of reinfection with respiratory syncytial virus
Caroline Breese Hall, Edward E Walsh, Christine E Long, and Kenneth C Schnabel. Immunity to and frequency of reinfection with respiratory syncytial virus. Journal of Infectious Diseases, 163(4):693–698, 1991
1991
-
[35]
Lyapunov functions and global stability for sir, sirs, and sis epi- demiological models
Andrei Korobeinikov and Graeme C Wake. Lyapunov functions and global stability for sir, sirs, and sis epi- demiological models. Applied Mathematics Letters , 15 (8):955–960, 2002
2002
-
[2011]
doi:10.1098/rspb.2011.0300
ISSN 0962-8452. doi:10.1098/rspb.2011.0300. URL http://rspb.royalsocietypublishing.org/content/ 278/1725/3635
2011
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.