Predicting West Nile virus transmission in North American bird communities using phylogenetic mixed effects models and eBird citizen science data

Background West Nile virus (WNV) is a mosquito-transmitted disease of birds that has caused bird population declines and can spill over into human populations. Previous research has identified bird species that infect a large fraction of the total pool of infected mosquitoes and correlate with human infection risk; however, these analyses cover small spatial regions and cannot be used to predict transmission in bird communities in which these species are rare or absent. Here we present a mechanistic model for WNV transmission that predicts WNV spread (R-0) in any bird community in North America by scaling up from the physiological responses of individual birds to transmission at the level of the community. We predict unmeasured bird species' responses to infection using phylogenetic imputation, based on these species' phylogenetic relationships with bird species with measured responses. Results We focused our analysis on Texas, USA, because it is among the states with the highest total incidence of WNV in humans and is well sampled by birders in the eBird database. Spatio-temporal patterns: WNV transmission is primarily driven by temperature variation across time and space, and secondarily by bird community composition. In Texas, we predicted WNV R-0 to be highest in the spring and fall when temperatures maximize the product of mosquito transmission and survival probabilities. In the most favorable months for WNV transmission (April, May, September and October), we predicted R-0 to be highest in the Piney Woods and Oak Woods & Prairies ecoregions of Texas, and lowest in the High Plains and South Texas Brush County ecoregions. Dilution effect: More abundant bird species are more competent hosts for WNV, and predicted WNV R-0 decreases with increasing species richness. Keystone species: We predicted that northern cardinals (Cardinalis cardinalis) are the most important hosts for amplifying WNV and that mourning doves (Zenaida macroura) are the most important sinks of infection across Texas. Conclusions Despite some data limitations, we demonstrate the power of phylogenetic imputation in predicting disease transmission in heterogeneous host communities. Our mechanistic modeling framework shows promise both for assisting future analyses on transmission and spillover in heterogeneous multispecies pathogen systems and for improving model transparency by clarifying assumptions, choices and shortcomings in complex ecological analyses.