in

Synthesizing selection mosaic theory and host-pathogen theory to explain large-scale pathogen coexistence


Abstract

Selection mosaic theory explains observations of polymorphism in host-pathogen interactions in terms of spatially variable natural selection but does not account for population dynamics. In contrast, classical host-pathogen theory easily explains observations of population cycles, but does not explain the persistence of pathogen polymorphism. Here, we synthesize these two frameworks to understand the effects of population cycles on pathogen polymorphism. We show that geographic variation in the frequency of two morphotypes of a baculovirus that infects the Douglas-fir tussock moth (Orgyia pseudotsugata) depends on the frequency of Douglas-fir (Pseudotsuga menziesii), an important tussock moth host tree. The morphotype frequency data are best explained by host-pathogen models that combine a selection mosaic with population cycles. In our model, population cycles intensify pathogen competition across a selection mosaic, leading to a strong effect of Douglas-fir frequency on morphotype frequency that matches the data. Models without host-pathogen cycles or a selection mosaic project only weak effects of varying Douglas-fir frequency. Our model further projects that a biopesticide made up of both viral morphotypes would be more effective than the current single-morphotype biopesticide, demonstrating that our synthesis of selection mosaic theory and host-pathogen theory provides useful insights into pest management.

Data availability

The raw data for the morphotype frequency dataset, field experiments, and line search results supporting the findings of this study are openly available in the GitHub repository at https://github.com/kpd19/Two_Pathogen_Evolution/, with a persistent identifier assigned to version 1.0.0 via Zenodo: https://doi.org/10.5281/zenodo.17574036. The Bayesian model outputs from Stan are very large and are available from the first author upon request. The data from the increased realizations from the line search results are also very large and are available from the first author upon request. The previously published data included as part of our morphotype frequency dataset can be found in Fig. 1 from Hughes37, Table 1 from Williams and Otvos65, and Table 1 from Williams et al.40. The National Forest Type Dataset for the continental United States was previously publicly available from the USDA Forest Service. The dataset is no longer hosted by the USDA and is available from the first author upon request. State and Province administrative boundaries for the United States and Canada used in the maps for this paper are publicly available for download from GADM v4.1 via https://geodata.ucdavis.edu/gadm/gadm4.1/shp/. Source data are provided with this paper.

Code availability

All code to perform the simulations and statistical analysis, as well as for plotting Figs. 1–6 and the supplementary data figures, is openly available in the GitHub repository at https://github.com/kpd19/Two_Pathogen_Evolution/, with a persistent identifier assigned to version 1.0.0 via Zenodo: https://doi.org/10.5281/zenodo.17574036.

References

  1. May, R. M. & Anderson, R. M. Epidemiology and genetics in the coevolution of parasites and hosts. Proc. R. Soc. Lond. B Biol. Sci. 219, 281–313 (1983).

    Google Scholar 

  2. Anderson, R. M. & May, R. M. Coevolution of hosts and parasites. Parasitology 85, 411–426 (1982).

    Google Scholar 

  3. Anis, E., Newell, T. K., Dyer, N. & Wilkes, R. P. Phylogenetic analysis of the wild-type strains of canine distemper virus circulating in the United States. Virol. J. 15, 1–10 (2018).

    Google Scholar 

  4. Troupin, C. et al. Large-scale phylogenomic analysis reveals the complex evolutionary history of rabies virus in multiple carnivore hosts. PLoS Pathog. 12, e1006041 (2016).

    Google Scholar 

  5. Karki, M., Rajak, K. K. & Singh, R. P. Canine morbillivirus (CDV): a review on current status, emergence and the diagnostics. Virusdisease 33, 309–321 (2022).

    Google Scholar 

  6. Byrne, A. Q. et al. Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis). PLoS ONE 17, e0261047 (2022).

    Google Scholar 

  7. Kaiser, K. & Pollinger, J. Batrachochytrium dendrobatidis shows high genetic diversity and ecological niche specificity among haplotypes in the maya mountains of belize. PLoS ONE 7, e32113 (2012).

    Google Scholar 

  8. Azat, C. Not just a pathogen: the importance of recognizing genetic variability to mitigate a wildlife pandemic. Mol. Ecol. Resour. 21, 1410–1412 (2021).

    Google Scholar 

  9. Drees, K. P. et al. Phylogenetics of a fungal invasion: origins and widespread dispersal of white-nose syndrome. mBio 8, e01941–17 (2017).

    Google Scholar 

  10. Ladner, J. T. et al. The population genetics of the causative agent of snake fungal disease indicate recent introductions to the USA. PLoS Biol. 20, e3001676 (2022).

    Google Scholar 

  11. Kerr, P. J. et al. Punctuated evolution of myxoma virus: rapid and disjunct evolution of a recent viral lineage in Australia. J. Virol. 93, e01994–18 (2019).

    Google Scholar 

  12. Laine, A.-L. Role of coevolution in generating biological diversity: spatially divergent selection trajectories. J. Exp. Bot. 60, 2957–2970 (2009).

    Google Scholar 

  13. Vander Wal, E. et al. Applying evolutionary concepts to wildlife disease ecology and management. Evolut. Appl. 7, 856–868 (2014).

    Google Scholar 

  14. Anderson, R. M. & May, R. M. Regulation and stability of host-parasite population interactions: I. regulatory processes. J. Anim. Ecol. 47, 219–247 (1978).

    Google Scholar 

  15. May, R. M. & Anderson, R. M. Regulation and stability of host-parasite population interactions: II. destabilizing processes. J. Anim. Ecol. 47, 249–267 (1978).

    Google Scholar 

  16. Dwyer, G., Dushoff, J. & Yee, S. H. The combined effects of pathogens and predators on insect outbreaks. Nature 430, 341–345 (2004).

    Google Scholar 

  17. Cassirer, E. F. et al. Spatio-temporal dynamics of pneumonia in bighorn sheep. J. Anim. Ecol. 82, 518–528 (2013).

    Google Scholar 

  18. Lips, K. R. Overview of chytrid emergence and impacts on amphibians. Philos. Trans. R. Soc. Lond. B Biol. Sci. 371, 20150465 (2016).

    Google Scholar 

  19. Cunningham, C. X. et al. Quantifying 25 years of disease-caused declines in tasmanian devil populations: host density drives spatial pathogen spread. Ecol. Lett. 24, 958–969 (2021).

    Google Scholar 

  20. Graham, A. L. et al. Fitness correlates of heritable variation in antibody responsiveness in a wild mammal. Science 330, 662–665 (2010).

    Google Scholar 

  21. New, L. F., Matthiopoulos, J., Redpath, S. & Buckland, S. T. Fitting models of multiple hypotheses to partial population data: investigating the causes of cycles in red grouse. Am. Nat. 174, 399–412 (2009).

    Google Scholar 

  22. Woods, S. A. & Elkinton, J. S. Bimodal patterns of mortality from nuclear polyhedrosis-virus in gypsy-moth (Lymantria dispar) populations. J. Invertebr. Pathol. 50, 151–157 (1987).

    Google Scholar 

  23. Andreasen, V. & Dwyer, G. Seasonality and the coexistence of pathogen strains. Am. Nat. 201, 639–658 (2023).

    Google Scholar 

  24. Thompson, J. N. & Burdon, J. J. Gene-for-gene coevolution between plants and parasites. Nature 360, 121–125 (1992).

    Google Scholar 

  25. Brown, J. K. M. & Tellier, A. Plant-parasite coevolution: bridging the gap between genetics and ecology. Annu. Rev. Phytopathol. 49, 345–367 (2011).

    Google Scholar 

  26. Tack, A. J. M., Thrall, P. H., Barrett, L. G., Burdon, J. J. & Laine, A.-L. Variation in infectivity and aggressiveness in space and time in wild host–pathogen systems: causes and consequences. J. Evol. Biol. 25, 1918–1936 (2012).

    Google Scholar 

  27. Vale, P. F. Killing them softly: managing pathogen polymorphism and virulence in spatially variable environments. Trends Parasitol. 29, 417–422 (2013).

    Google Scholar 

  28. Gavrilets, S. & Michalakis, Y. Effects of environmental heterogeneity on victim–exploiter coevolution. Evolution 62, 3100–3116 (2008).

    Google Scholar 

  29. Tellier, A. & Brown, J. K. M. Spatial heterogeneity, frequency-dependent selection and polymorphism in host-parasite interactions. BMC Evol. Biol. 11, 1–15 (2011).

    Google Scholar 

  30. Thompson, J. N. The Geographic Mosaic of Coevolution (University of Chicago Press, 2005).

  31. Nuismer, S. L. Parasite local adaptation in a geographic mosaic. Evolution 60, 24–30 (2006).

    Google Scholar 

  32. Gandon, S. & Nuismer, S. Interactions between genetic drift, gene flow, and selection mosaics drive parasite local adaptation. Am. Nat. 173, 212–224 (2009).

    Google Scholar 

  33. Papaïx, J., Burdon, J. J., Walker, E., Barrett, L. G. & Thrall, P. H. Metapopulation structure predicts population dynamics in the cakile maritima–alternaria brassicicola host-pathogen interaction. Am. Nat. 197, E55–E71 (2021).

    Google Scholar 

  34. Mihaljevic, J. R. et al. An empirical test of the role of small-scale transmission in large-scale disease dynamics. Am. Nat. 195, 616–635 (2020).

    Google Scholar 

  35. Dwyer, G., Dushoff, J., Elkinton, J. S. & Levin, S. A. Pathogen-driven outbreaks in forest defoliators revisited: building models from experimental data. Am. Nat. 156, 105–120 (2000).

    Google Scholar 

  36. Mason, R. R. Dynamic behavior of Douglas-fir tussock moth populations in the Pacific northwest. Science 42, 182–191 (1996).

    Google Scholar 

  37. Hughes, K. M. Notes on the nuclear polyhedrosis viruses of tussock moths of the genus Orgyia (Lepidoptera). Can. Entomol. 108, 479–484 (1976).

    Google Scholar 

  38. Jakubowska, A., van Oers, M. M., Otvos, I. S. & Vlak, J. M. Phylogenetic analysis of Orgyia pseudotsugata single-nucleocapsid nucleopolyhedrovirus. Virol. Sin. 22, 257–265 (2007).

    Google Scholar 

  39. Hughes, K. M. & Addison, R. B. Two nuclear polyhedrosis viruses of Douglas-fir tussock moth. J. Invertebr. Pathol. 16, 196–204 (1970).

    Google Scholar 

  40. Williams, H. L., Monge-Monge, K. S., Otvos, I. S., Reardon, R. & Ragenovich, I. Genotypic variation among Douglas-fir tussock moth nucleopolyhedrovirus (OpNPV) isolates in the western United States. J. Invertebr. Pathol. 108, 13–21 (2011).

    Google Scholar 

  41. Elderd, B. D. & Dwyer, G. Using insect baculoviruses to understand how population structure affects disease spread. in Wildlife disease ecology: linking theory to data and application, 225–261 (Cambridge University Press, 2019).

  42. Cory, J. S. & Hoover, K. Plant-mediated effects in insect-pathogen interactions. Trends Ecol. Evol. 21, 278–286 (2006).

    Google Scholar 

  43. Freedman, A. S., Huang, A. Y., Dixon, K. P., Polivka, C. & Dwyer, G. Effects of host-tree foliage on polymorphism in an insect pathogen. Am. Nat. Preprint at https://www.biorxiv.org/content/10.1101/2024.09.29.615698v1 (2025).

  44. Swor, K., Satyal, P., Poudel, A. & Setzer, W. N. Gymnosperms of Idaho: chemical compositions and enantiomeric distributions of essential oils of Abies lasiocarpa, Picea engelmannii, Pinus contorta, Pseudotsuga menziesii, and Thuja plicata. Molecules 28, 2477 (2023).

    Google Scholar 

  45. Blackwood, J. C., Streicker, D. G., Altizer, S. & Rohani, P. Resolving the roles of immunity, pathogenesis, and immigration for rabies persistence in vampire bats. Proc. Natl. Acad. Sci. USA 110, 20837–20842 (2013).

    Google Scholar 

  46. Martignoni, M. E. History of TM Biocontrol-1, the first registered virus-based produced for control of a forest insect. Am. Entomol. 45, 30–37 (1999).

    Google Scholar 

  47. Brookes, M. H., Stark, R. W. & Campbell, R. W. The Douglas-fir tussock moth: a synthesis. 1585 (US Department of Agriculture, 1978).

  48. Fleming-Davies, A. E., Dukic, V., Andreasen, V. & Dwyer, G. Effects of host heterogeneity on pathogen diversity and evolution. Ecol. Lett. 18, 1252–1261 (2015).

    Google Scholar 

  49. Watanabe, H. The host population. in (eds Fuxa, J. R. & Tanada, Y.) Epizootiology of Insect Diseases, 71–112 (Wiley-Interscience, 1987).

  50. Páez, D. J., Dukic, V., Dushoff, J., Fleming-Davies, A. & Dwyer, G. Eco-evolutionary theory and insect outbreaks. Am. Nat. 189, 616–629 (2017).

    Google Scholar 

  51. Mitchell, R. G. Dispersal of early instars of the Douglas-fir tussock moth. Ann. Entomol. Soc. Am. 72, 291–297 (1979).

    Google Scholar 

  52. Wetzel, W. C., Kharouba, H. M., Robinson, M., Holyoak, M. & Karban, R. Variability in plant nutrients reduces insect herbivore performance. Nature 539, 425–427 (2016).

    Google Scholar 

  53. Reilly, J. R. & Elderd, B. D. Effects of biological control on long-term population dynamics: identifying unexpected outcomes. J. Appl. Ecol. 51, 90–101 (2014).

    Google Scholar 

  54. Duffy, M. A. & Sivars-Becker, L. Rapid evolution and ecological host-parasite dynamics. Ecol. Lett. 10, 44–53 (2007).

    Google Scholar 

  55. Gomes, M. G. M. et al. How host heterogeneity governs tuberculosis reinfection? Proc. R. Soc. B Biol. Sci. 279, 2473–2478 (2012).

    Google Scholar 

  56. Kain, M. P., Cattadori, I. M. & Bolker, B. M. The evolutionary response of virulence to host heterogeneity: a general model with application to myxomatosis in rabbits co-infected with intestinal helminths. Evol. Ecol. Res. 19, 257–278 (2018).

    Google Scholar 

  57. Kendall, B. E. et al. Why do populations cycle? a synthesis of statistical and mechanistic modeling approaches. Ecology 80, 1789–1805 (1999).

    Google Scholar 

  58. Abbott, K. C. & Dwyer, G. Food limitation and insect outbreaks: complex dynamics in plant-herbivore models. J. Anim. Ecol. 76, 1004–1014 (2007).

    Google Scholar 

  59. Otvos, I. S., Cunningham, J. C. & Friskie, L. M. Aerial application of nuclear polyhedrosis virus against Douglas-fir tussock moth, Orgyia pseudostugata (Mcdunnough) (Lepidoptera: Lymantriidae). 1. Impact in the year of application. Can. Entomol. 119, 697–706 (1987).

    Google Scholar 

  60. Dixon, K. P. Interacting effects of host-pathogen ecology and evolution and climate change on outbreaks of a forest pest insect. Ph.D. thesis, University of Chicago (2024).

  61. Cotter, S. C., Kruuk, L. E. B. & Wilson, K. Costs of resistance: genetic correlations and potential trade-offs in an insect immune system. J. Evol. Biol. 17, 421–429 (2004).

    Google Scholar 

  62. Tabashnik, B. E. Resistance risk assessment: realized heritability of resistance to Bacillus thuringiensis in diamondback moth (Lepidoptera: Plutellidae), tobacco budworm (Lepidoptera: Noctuidae), and colorado potato beetle (Coleoptera: Chrysomelidae). J. Econ. Entomol. 85, 1551–1559 (1992).

    Google Scholar 

  63. Decanini, L. I., Collins, A. M. & Evans, J. D. Variation and heritability in immune gene expression by diseased honeybees. J. Hered. 98, 195–201 (2007).

    Google Scholar 

  64. Alizon, S. & van Baalen, M. Multiple infections, immune dynamics, and the evolution of virulence. Am. Nat. 172, E150–E168 (2008).

    Google Scholar 

  65. Williams, H. L. & Otvos, I. S. Genotypic variation and presence of rare genotypes among Douglas-fir tussock moth multicapsid nucleopolyhedrovirus (OpMNPV) isolates in British Columbia. J. Invertebr. Pathol. 88, 190–200 (2005).

    Google Scholar 

  66. Ruefenacht, B. et al. Conterminous us and alaska forest type mapping using forest inventory and analysis data. Photogramm. Eng. Remote Sens. 74, 1379–1388 (2008).

    Google Scholar 

  67. Beaudoin, A. et al. Mapping attributes of canada’s forests at moderate resolution through k nn and modis imagery. Can. J. Res. 44, 521–532 (2014).

    Google Scholar 

  68. Birch, C. P. D., Oom, S. P. & Beecham, J. A. Rectangular and hexagonal grids used for observation, experiment and simulation in ecology. Ecol. Model. 206, 347–359 (2007).

    Google Scholar 

  69. Berger, J. O., Liseo, B. & Wolpert, R. L. Integrated likelihood methods for eliminating nuisance parameters. Stat. Sci. 14, 1–28 (1999).

    Google Scholar 

  70. Kennedy, D., Dukic, V. & Dwyer, G. Pathogen growth in insect hosts: inferring the importance of different mechanisms using stochastic models and response-time data. Am. Nat. 184, 407–423 (2014).

    Google Scholar 

Download references

Acknowledgements

We are extremely grateful for the support of many dedicated and talented field technicians: Rhiannon Archerelle, Ari Freedman, Amy Gannon, Laurel Haavik, Sophia Horigan, Amy Huang, Alison Hunter, Jessica Johnson, August Kramer, Allie Kreitman, Kate-Lynne Logan, and Chelsea Miller. Special thanks to the field technicians who made the work possible in the Summer of 2020, as well as Mary Johnson and Luis Marmolejo, who managed shipping. We thank Roy Magelssen and Connie Mehmel at the Forestry Sciences Laboratory in Wenatchee, Washington, for providing important institutional and biological knowledge of the system, and Joe Mihaljevic for important guidance in doing transmission experiments. We would like to thank Cara Brook, Sarah Cobey, and Tim Wootton for providing important feedback on earlier drafts of this work. Computational Resources were provided by the Research Computing Center at the University of Chicago. Help in morphotyping isolates was provided by the Electron Microscopy Center at the University of Chicago. Our work was supported by EEID NSF grant DEB-2109774 to G.D. and V.D. K.P.D. was supported by the University of Chicago Data Science for Energy and Environmental Research (DSEER) training grant as part of an NSF Research Traineeship program (1735359) and the U.S. Department of Education Quantitative Ecology GAANN training grant (P200A150101). W.T.K. and K.P.D. received separate awards from the University of Chicago Hinds Fund for Student Research. Our work was further supported by a grant from the Western Wildlands Environmental Threat Assessment Center to C.M.P., by NIFA Biological Sciences grant 2019-67014-29919 to V.D., by an ARCS Foundation Fellowship to W.T.K., and by a Theodore Roosevelt Memorial Grant through the American Museum of Natural History to W.T.K.

Author information

Authors and Affiliations

Authors

Contributions

G.D. and K.P.D. planned and designed the research. W.T.K., C.M.P., G.B., G.D., and K.P.D. collected the data. K.P.D. analyzed the data. K.P.D., G.D., V.D., and C.M.P. contributed substantially to the discussion of the results and the writing of the manuscript.

Corresponding author

Correspondence to
Greg Dwyer.

Ethics declarations

Competing interests

The authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks Samuel Hulse and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Additional information

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

Supplementary information

Supplementary Information

Reporting Summary

Transparent Peer Review File

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Cite this article

Dixon, K.P., Koval, W.T., Polivka, C.M. et al. Synthesizing selection mosaic theory and host-pathogen theory to explain large-scale pathogen coexistence.
Nat Commun (2025). https://doi.org/10.1038/s41467-025-67952-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41467-025-67952-6


Source: Ecology - nature.com

Species richness and Spatial distribution of three Pieridae subfamilies across mainland China under past and future climates

A comprehensive approach to enhancing irrigation network management through the water accounting plus framework

Back to Top