A Maximum Entropy Model of the Distribution of Dengue Serotype in Mexico

Pathogen strain diversity is an important driver of the trajectory of epidemics. The role of bioclimatic factors on the spatial distribution of dengue virus (DENV) serotypes has, however, not been previously studied. Hence, we developed municipality-scale environmental suitability maps for the four dengue virus serotypes using maximum entropy modeling. We fit climatic variables to municipality presence records from 2012 to 2020 in Mexico. Bioclimatic variables were explored for their environmental suitability to different DENV serotypes, and the different distributions were visualized using three cutoff probabilities representing 90%, 95%, and 99% sensitivity. Municipality-level results were then mapped in ArcGIS. The overall accuracy for the predictive models was 0.69, 0.68, 0.75, and 0.72 for DENV-1, DENV-2, DENV-3, and DENV-4, respectively. Important predictors of all DENV serotypes were the growing degree days for December, January, and February, which are an indicator of higher temperatures and the precipitation of the wettest month. The minimum temperature of the coldest month between -5°C and 20°C was found to be suitable for DENV-1 and DENV-2 serotypes. Respectively, above 700-900 mm of rainfall, the suitability for DENV-1 and DENV-2 begins to decline, while higher humidity still favors DENV-3 and DENV-4. The sensitivity concerning the suitability map was developed for Mexico. DENV-1, DENV-2, DENV-3, and DENV-4 serotypes will be found more commonly in the municipalities classified as suitable based on their respective sensitivity of 91%, 90%, 89%, and 85% in Mexico. As the microclimates continue to change, specific bioclimatic indices may be used to monitor potential changes in DENV serotype distribution. The suitability for DENV-1 and DENV-2 is expected to increase in areas with lower minimum temperature ranges, while DENV-3 and DENV-4 will likely increase in areas that experience higher humidity. Ongoing surveillance of municipalities with predicted suitability of 89% and 85% should be expanded to account for the accurate DENV serotype prevalence and association between bioclimatic parameters.