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
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).
Anderson, R. M. & May, R. M. Coevolution of hosts and parasites. Parasitology 85, 411–426 (1982).
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).
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).
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).
Byrne, A. Q. et al. Host species is linked to pathogen genotype for the amphibian chytrid fungus (Batrachochytrium dendrobatidis). PLoS ONE 17, e0261047 (2022).
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).
Azat, C. Not just a pathogen: the importance of recognizing genetic variability to mitigate a wildlife pandemic. Mol. Ecol. Resour. 21, 1410–1412 (2021).
Drees, K. P. et al. Phylogenetics of a fungal invasion: origins and widespread dispersal of white-nose syndrome. mBio 8, e01941–17 (2017).
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).
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).
Laine, A.-L. Role of coevolution in generating biological diversity: spatially divergent selection trajectories. J. Exp. Bot. 60, 2957–2970 (2009).
Vander Wal, E. et al. Applying evolutionary concepts to wildlife disease ecology and management. Evolut. Appl. 7, 856–868 (2014).
Anderson, R. M. & May, R. M. Regulation and stability of host-parasite population interactions: I. regulatory processes. J. Anim. Ecol. 47, 219–247 (1978).
May, R. M. & Anderson, R. M. Regulation and stability of host-parasite population interactions: II. destabilizing processes. J. Anim. Ecol. 47, 249–267 (1978).
Dwyer, G., Dushoff, J. & Yee, S. H. The combined effects of pathogens and predators on insect outbreaks. Nature 430, 341–345 (2004).
Cassirer, E. F. et al. Spatio-temporal dynamics of pneumonia in bighorn sheep. J. Anim. Ecol. 82, 518–528 (2013).
Lips, K. R. Overview of chytrid emergence and impacts on amphibians. Philos. Trans. R. Soc. Lond. B Biol. Sci. 371, 20150465 (2016).
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).
Graham, A. L. et al. Fitness correlates of heritable variation in antibody responsiveness in a wild mammal. Science 330, 662–665 (2010).
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).
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).
Andreasen, V. & Dwyer, G. Seasonality and the coexistence of pathogen strains. Am. Nat. 201, 639–658 (2023).
Thompson, J. N. & Burdon, J. J. Gene-for-gene coevolution between plants and parasites. Nature 360, 121–125 (1992).
Brown, J. K. M. & Tellier, A. Plant-parasite coevolution: bridging the gap between genetics and ecology. Annu. Rev. Phytopathol. 49, 345–367 (2011).
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).
Vale, P. F. Killing them softly: managing pathogen polymorphism and virulence in spatially variable environments. Trends Parasitol. 29, 417–422 (2013).
Gavrilets, S. & Michalakis, Y. Effects of environmental heterogeneity on victim–exploiter coevolution. Evolution 62, 3100–3116 (2008).
Tellier, A. & Brown, J. K. M. Spatial heterogeneity, frequency-dependent selection and polymorphism in host-parasite interactions. BMC Evol. Biol. 11, 1–15 (2011).
Thompson, J. N. The Geographic Mosaic of Coevolution (University of Chicago Press, 2005).
Nuismer, S. L. Parasite local adaptation in a geographic mosaic. Evolution 60, 24–30 (2006).
Gandon, S. & Nuismer, S. Interactions between genetic drift, gene flow, and selection mosaics drive parasite local adaptation. Am. Nat. 173, 212–224 (2009).
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).
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).
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).
Mason, R. R. Dynamic behavior of Douglas-fir tussock moth populations in the Pacific northwest. Science 42, 182–191 (1996).
Hughes, K. M. Notes on the nuclear polyhedrosis viruses of tussock moths of the genus Orgyia (Lepidoptera). Can. Entomol. 108, 479–484 (1976).
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).
Hughes, K. M. & Addison, R. B. Two nuclear polyhedrosis viruses of Douglas-fir tussock moth. J. Invertebr. Pathol. 16, 196–204 (1970).
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).
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).
Cory, J. S. & Hoover, K. Plant-mediated effects in insect-pathogen interactions. Trends Ecol. Evol. 21, 278–286 (2006).
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).
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).
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).
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).
Brookes, M. H., Stark, R. W. & Campbell, R. W. The Douglas-fir tussock moth: a synthesis. 1585 (US Department of Agriculture, 1978).
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).
Watanabe, H. The host population. in (eds Fuxa, J. R. & Tanada, Y.) Epizootiology of Insect Diseases, 71–112 (Wiley-Interscience, 1987).
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).
Mitchell, R. G. Dispersal of early instars of the Douglas-fir tussock moth. Ann. Entomol. Soc. Am. 72, 291–297 (1979).
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).
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).
Duffy, M. A. & Sivars-Becker, L. Rapid evolution and ecological host-parasite dynamics. Ecol. Lett. 10, 44–53 (2007).
Gomes, M. G. M. et al. How host heterogeneity governs tuberculosis reinfection? Proc. R. Soc. B Biol. Sci. 279, 2473–2478 (2012).
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).
Kendall, B. E. et al. Why do populations cycle? a synthesis of statistical and mechanistic modeling approaches. Ecology 80, 1789–1805 (1999).
Abbott, K. C. & Dwyer, G. Food limitation and insect outbreaks: complex dynamics in plant-herbivore models. J. Anim. Ecol. 76, 1004–1014 (2007).
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).
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).
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).
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).
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).
Alizon, S. & van Baalen, M. Multiple infections, immune dynamics, and the evolution of virulence. Am. Nat. 172, E150–E168 (2008).
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).
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).
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).
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).
Berger, J. O., Liseo, B. & Wolpert, R. L. Integrated likelihood methods for eliminating nuisance parameters. Stat. Sci. 14, 1–28 (1999).
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).
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
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
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
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-025-67952-6
Source: Ecology - nature.com
