Dynamic shifts of functional diversity through climate-Resilient strategies and farmland restoration in a mountain protected area.

Land-use land-cover (LULC) change contributes to major ecological impacts, particularly in areas undergoing land abandonment, inducing modifications on habitat structure and species distributions. Alternative land-use policies are potential solutions to alleviate the negative impacts of contemporary tendencies of LULC change on biodiversity. This work analyzes these tendencies in the Montesinho Natural Park (Portugal), an area representative of European abandoned mountain rural areas. We built ecological niche models for 226 species of vertebrates (amphibians, reptiles, birds, and mammals) and vascular plants, using a consensus modelling approach available in the R package 'biomod2'. We projected the models to contemporary (2018) and future (2050) LULC scenarios, under four scenarios aiming to secure relevant ecosystem services and biodiversity conservation for 2050: an afforestation and a rewilding scenario, focused on climate-smart management strategies, and a farmland and an agroforestry recovery scenario, based on re-establishing human traditional activities. We quantified the influences of these scenarios on biodiversity through species habitat suitability changes for 2018-2050. We analyzed how these management strategies could influence indices of functional diversity (functional richness, functional evenness and functional dispersion) within the park. Habitat suitability changes revealed complementary patterns among scenarios. Afforestation and rewilding scenarios benefited more species adapted to habitats with low human influence, such as forests and open woodlands. The highest functional richness and dispersion was predicted for rewilding scenarios, which could improve landscape restoration and provide opportunities for the expansion and recolonization of forest areas by native species. The recovery of traditional farming and agroforestry activities results in the lowest values of functional richness, but these strategies contribute to complex landscape matrices with diversified habitats and resources. Moreover, this strategy could offer opportunities for fire suppression and increase landscape fire resistance. An integrative approach reconciling rewilding initiatives with the recovery of extensive agricultural and agroforestry activities is potentially an harmonious strategy for supporting the provision of ecosystem services while securing biodiversity conservation and functional diversity within the natural park.

Forecasting habitat suitability and niche shifts of two global maize pests: Ostrinia furnacalis and Ostrinia nubilalis (Lepidoptera: Crambidae).

BACKGROUND: Ostrinia furnacalis (ACB) and Ostrinia nubilalis (ECB) are devastating pests of the agricultural crop maize worldwide. However, little is known about their potential distribution and niche shifts during their global invasion. Since long-term selection to past climate variability has shaped their historical niche breadth, such niche shifts may provide an alternative basis for understanding their responses to present and future climate change. By integrating the niche unfilling, stability, and expansion situations into a single framework, our study quantifies the patterns of niche shift in the spatial distribution of these two pests during the different periods. RESULTS: Our results show that the overall suitable habitats of ACB and ECB in the future decrease but highly and extremely suitable habitat will become more widespread, suggesting these two insects may occur more frequently in specific regions. Compared with Southeast Asia and Australia, the ACB niche in China exhibited expansion rather than unfilling. For ECB, initial niches have a tendency to be retained in Eurasia despite there also being potential for expansion in North America. The niche equivalency and similarity test results further indicate that niche shifts were common for both ACB and ECB in different survival regions during their colonization of new habitat and their suitable habitat changes during the paleoclimate were associated with climatic changes. CONCLUSIONS: These findings improve our understanding of the ecological characteristics of ACB and ECB worldwide, and will be useful in the development of prevention and control strategies for two insect pests worldwide. © 2024 Society of Chemical Industry.

Can upwelling regions be potential thermal refugia for marine fishes during climate warming?

Species are expected to migrate to higher latitudes as warming intensifies due to anthropogenic climate change since physiological mechanisms have been adapted to maximize fitness under specific temperatures. However, literature suggests that upwellings could act as thermal refugia under climate warming protecting marine ecosystem diversity. This research aimed to predict the effects of climate warming on commercial and non-commercial fish species reported in official Mexican documents (>200 species) based on their thermal niche to observe if upwellings can act as potential thermal refugia. Present (2000-2014) and Representative Concentration Pathway (6.0 and 8.5) scenarios (2040-2050 and 2090-2100) have been considered for this work. Current and future suitability patterns, species distribution, richness, and turnover were calculated using the minimum volume ellipsoids as algorithm. The results in this study highlight that beyond migration to higher latitudes, upwelling regions could protect marine fishes, although the mechanism differed between the innate characteristics of upwellings. Most modeled species (primarily tropical fishes) found refuge in the tropical upwelling in Northern Yucatan. However, the highest warming scenario overwhelmed this region. In contrast, the Baja California region lies within the Eastern Boundary Upwelling Systems. While the area experiences an increase in suitability, the northern regions have a higher upwelling intensity acting as environmental barriers for many tropical species. Conversely, in the southern regions where upwelling is weaker, species tend to congregate and persist even during elevated warming, according to the turnover analysis. These findings suggest that tropicalization in higher latitudes may not be as straightforward as previously assumed. Nevertheless, climate change affects numerous ecosystem features, such as trophic relationships, phenology, and other environmental variables not considered here. In addition, uncertainty still exists about the assumption of increasing intensity of upwelling systems.

Leaf functional traits and ecological niche of Fagus grandifolia and Oreomunnea mexicana in natural forests and plantings as a proxy of climate change.

PREMISE: Functional traits reflect species' responses to environmental variation and the breadth of their ecological niches. Fagus grandifolia and Oreomunnea mexicana have restricted distribution in upper montane cloud forests (1700-2000 m a.s.l.) in Mexico. These species were introduced into plantings at lower elevations (1200-1600 m a.s.l.) that have climates predicted for montane forests in 2050 and 2070. The aim was to relate morphological leaf traits to the ecological niche structure of each species. METHODS: Leaf functional traits (leaf area, specific leaf area [SLA], thickness, and toughness) were analyzed in forests and plantings. Atmospheric circulation models and representative concentration pathways (RCPs: 2.6, 4.5, 8.5) were used to assess future climate conditions. Trait-niche relationships were analyzed by measuring the Mahalanobis distance (MD) from the forests and the plantings to the ecological niche centroid (ENC). RESULTS: For both species, leaf area and SLA were higher and toughness lower in plantings at lower elevation relative to those in higher-elevation forests, and thickness was similar. Leaf traits varied with distance from sites to the ENC. Forests and plantings have different environmental locations regarding the ENC, but forests are closer (MD 0.34-0.58) than plantings (MD 0.50-0.70) for both species. CONCLUSIONS: Elevation as a proxy for expected future climate conditions influenced the functional traits of both species, and trait patterns related to the structure of their ecological niches were consistent. The use of distances to the ENC is a promising approach to explore variability in species' functional traits and phenotypic responses in optimal versus marginal environmental conditions.

Identifying Areas of Invasion Risk and Changes in the Ecological Niche Occupied by the Coffee Leaf Miner Leucoptera coffeella (Lepidoptera: Lyonetiidae).

Insects of economic importance such as Leucoptera coffeella can cause high defoliation in plants and reduce crop yields. We aimed to identify changes in the ecological niche and potential zones of the invasion. Occurrence records were obtained from databases and bibliography. WorldClim V2.0 bioclimatic layers were used. For the modeling of the potential distribution, the kuenm R package was used by executing the Maxent algorithm. The potential distribution models suggested greatest environmental suitability extends from Europe, South Asia, and Central and South Africa, showing the "tropical and subtropical moist broadleaf forests" as the ecoregion that presents the greatest probability of the presence of L. coffeella. The potential distribution model projected in the invaded area agrees with the known distribution in the region (America), although the results show that it is occupying environmental spaces not present in the area of origin. This species presented a large proportion of the invaded niche that overlaps the native niche and is colonizing new environmental conditions in the invaded area relative to its native distribution (Africa). This information could be used in the planning of coffee crops on the American continent.

The cumulative niche approach: A framework to assess the performance of ecological niche model projections.

Ecological Niche Models (ENMs) are often used to project species distributions within alien ranges and in future climatic scenarios. However, ENMs depend on species-environment equilibrium, which may be absent for actively expanding species. We present a novel framework to estimate whether species have reached environmental equilibrium in their native and alien ranges. The method is based on the estimation of niche breadth with the accumulation of species occurrences. An asymptote will indicate exhaustive knowledge of the realised niches. We demonstrate the CNA framework for 26 species of mammals, amphibians, and birds. Possible outcomes of the framework include: (1) There is enough data to quantify the native and alien realised niches, allowing us to calculate niche expansion between the native and alien ranges, also indicating that ENMs can be reliably projected to new environmental conditions. (2) The data in the native range is not adequate but an asymptote is reached in the alien realised niche, indicating low confidence in our ability to evaluate niche expansion in the alien range but high confidence in model projections to new environmental conditions within the alien range. (3) There is enough data to quantify the native realised niche, but not enough knowledge about the alien realised niche, hindering the reliability of projections beyond sampled conditions. (4) Both the native and alien ranges do not reach an asymptote, and thus few robust conclusions about the species' niche or future projections can be made. Our framework can be used to detect species' environmental equilibrium in both the native and alien ranges, to quantify changes in the realised niche during the invasion processes, and to estimate the likely accuracy of model projections to new environmental conditions.

Climate change may drive the distribution of tribe Zyginelline pests in China and the Indo-China Peninsula to shift towards higher latitude river-mountain systems.

BACKGROUND: Tribe Zyginelline leafhoppers can transmit plant viruses and are important pests that affect agriculture, forestry, and animal husbandry, causing serious economic losses. The potential distribution patterns of Zyginellini will change under climate change. Therefore, the best-performing random forest and maximum entropy models among 12 commonly used ecological niche models, alongside an ensemble model, were selected to predict the changes in habitat suitability distribution of Zyginellini under current and future climate scenarios [represented by two shared socio-economic pathways (SSPs), namely SSP126 and SSP585, for three periods (2050s, 2070s, and 2090s)] in China and the Indo-China Peninsula for the first time. RESULTS: The results revealed that the distribution of Zyginellini was mainly dominated by minimum temperature of coldest month. Under current and future climate scenarios, Zyginellini was mostly distributed southeast of the 400 mm equivalent precipitation line in China, and Vietnam. Under the future SSP126 scenario, the alert areas will mainly be concentrated in Hunan, Jiangxi, Zhejiang, Anhui, and Hebei in China, alongside Myanmar and Thailand in the Indo-China Peninsula. Meanwhile, in the SSP585 scenario, the alert areas in China will increase, whereas there will be little change in the Indo-China Peninsula. Interestingly, from the current to the future, the cores of Zyginelline distribution occurred around rivers and mountains, and shifted from Guizhou along the Yuanjiang River system to higher latitudes in Hunan. CONCLUSION: Zyginellini prefers higher latitude river-mountain systems under climate change. Our results will contribute to effective pest control strategies and biogeographical research for Zyginellini alongside other Cicadellidae insects. © 2023 Society of Chemical Industry.

Predicting current and future potential distributions of the greater bandicoot rat (Bandicota indica) under climate change conditions.

BACKGROUND: Rodent infestation is a global problem. Rodents cause huge harm to agriculture, forestry, and animal husbandry around the world and spread various zoonoses. In this study, we simulated the potentially suitable habitats of Bandicota indica and predicted the impact of future climate change on its distribution under different socio-economic pathway scenarios of CMIP6 using a parameter-optimized maximum entropy (MaxEnt) model. RESULTS: The average area under the receiver operating characteristic curve (AUC) value (0.958 ± 0.006) after ten repetitions proved the high accuracy of the MaxEnt model. Model results show that the annual mean temperature (≥ 15.93 °C), isothermality (28.52-80.49%), annual precipitation (780.13-3863.13 mm), precipitation of the warmest quarter (≥ 204.37 mm), and nighttime light (≥ 3.38) were important limiting environmental variables for the distribution of B. indica. Under current climate conditions, the projected potential suitable habitats for B. indica were mainly in India, China, Myanmar, Thailand, and Vietnam, which cover a total area of 301.70 × 104 km2 . The potentially suitable areas of B. indica in the world will expand under different future climate change scenarios by 1.61-17.65%. CONCLUSIONS: These results validate the potential influence of climate change on the distribution of B. indica and aid in understanding the linkages between B. indica niches and the relevant environment, thereby identifying urgent management areas where interventions may be necessary to develop feasible early warning and prevention strategies to protect against this rodent's spread. © 2023 Society of Chemical Industry.

The risk of vector transmission of Trypanosoma cruzi remains high in the State of Paraná

BACKGROUND Monitoring and analysing the infection rates of the vector of Trypanosoma cruzi, that causes Chagas disease, helps assess the risk of transmission. OBJECTIVES A study was carried out on triatomine in the State of Paraná, Brazil, between 2012 and 2021 and a comparison was made with a previous study. This was done to assess the risk of disease transmission. METHODS Ecological niche models based on climate and landscape variables were developed to predict habitat suitability for the vectors as a proxy for risk of occurrence. FINDINGS A total of 1,750 specimens of triatomines were recorded, of which six species were identified. The overall infection rate was 22.7%. The areas with the highest risk transmission of T. cruzi are consistent with previous predictions in municipalities. New data shows that climate models are more accurate than landscape models. This is likely because climate suitability was higher in the previous period. MAIN CONCLUSION Regardless of uneven sampling and potential biases, risk remains high due to the wide presence of infected vectors and high environmental suitability for vector species throughout the state and, therefore, improvements in public policies aimed at wide dissemination of knowledge about the disease are recommended to ensure the State remains free of Chagas disease.

Microclimate reveals the true thermal niche of forest plant species.

Species distributions are conventionally modelled using coarse-grained macroclimate data measured in open areas, potentially leading to biased predictions since most terrestrial species reside in the shade of trees. For forest plant species across Europe, we compared conventional macroclimate-based species distribution models (SDMs) with models corrected for forest microclimate buffering. We show that microclimate-based SDMs at high spatial resolution outperformed models using macroclimate and microclimate data at coarser resolution. Additionally, macroclimate-based models introduced a systematic bias in modelled species response curves, which could result in erroneous range shift predictions. Critically important for conservation science, these models were unable to identify warm and cold refugia at the range edges of species distributions. Our study emphasizes the crucial role of microclimate data when SDMs are used to gain insights into biodiversity conservation in the face of climate change, particularly given the growing policy and management focus on the conservation of refugia worldwide.

Climate change should drive mammal defaunation in tropical dry forests.

Human-induced climate change has intensified negative impacts on socioeconomic factors, the environment, and biodiversity, including changes in rainfall patterns and an increase in global average temperatures. Drylands are particularly at risk, with projections suggesting they will become hotter, drier, and less suitable for a significant portion of their species, potentially leading to mammal defaunation. We use ecological niche modelling and community ecology biodiversity metrics to examine potential geographical range shifts of non-volant mammal species in the largest Neotropical dryland, the Caatinga, and evaluate impacts of climate change on mammal assemblages. According to projections, 85% of the mammal species will lose suitable habitats, with one quarter of species projected to completely lose suitable habitats by 2060. This will result in a decrease in species richness for more than 90% of assemblages and an increase in compositional similarity to nearby assemblages (i.e., reduction in spatial beta diversity) for 70% of the assemblages. Small-sized mammals will be the most impacted and lose most of their suitable habitats, especially in highlands. The scenario is even worse in the eastern half of Caatinga where habitat destruction already prevails, compounding the threats faced by species there. While species-specific responses can vary with respect to dispersal, behavior, and energy requirements, our findings indicate that climate change can drive mammal assemblages to biotic homogenization and species loss, with drastic changes in assemblage trophic structure. For successful long-term socioenvironmental policy and conservation planning, it is critical that findings from biodiversity forecasts are considered.

Phylogeography and Genetic Structure of the Bush Cricket Decma fissa (Orthoptera, Tettigoniidae) in Southern China.

Decma fissa is the most widely distributed species of the genus Decma occuring in southern China. This study presents the first phylogeographic work of D. fissa based on COI, Cytb and ITS sequence. We examined genetic diversity with ITS and mitochondrial sequence respectively, and phylogenetic work was based on the mitochondrial data. A high-level genetic diversity was revealed based on mitochondrial data but a low-level diversity was shown with ITS sequence. For the mitochondrial data, divergence time analysis displayed five lineages. Based on the Mantel test, geographic and genetic distances among D. fissa populations revealed a significant positive correlation. Bayesian skyline plot (BSP) analyses implied that none of three major lineages of D. fissa was seemingly affected by the last glacial maximum (LGM, 0.015-0.025 Mya). Ecological niche modeling was used to predict the distribution of D. fissa in four periods (LGM, Mid-Holocene, current and 2070) in China. Analysis of the ancestral area reconstruction indicated that D. fissa occurred in the South China area.

Different facets of the same niche: Integrating citizen science and scientific survey data to predict biological invasion risk under multiple global change drivers.

Citizen science initiatives have been increasingly used by researchers as a source of occurrence data to model the distribution of alien species. Since citizen science presence-only data suffer from some fundamental issues, efforts have been made to combine these data with those provided by scientifically structured surveys. Surprisingly, only a few studies proposing data integration evaluated the contribution of this process to the effective sampling of species' environmental niches and, consequently, its effect on model predictions on new time intervals. We relied on niche overlap analyses, machine learning classification algorithms and ecological niche models to compare the ability of data from citizen science and scientific surveys, along with their integration, in capturing the realized niche of 13 invasive alien species in Italy. Moreover, we assessed differences in current and future invasion risk predicted by each data set under multiple global change scenarios. We showed that data from citizen science and scientific surveys captured similar species niches though highlighting exclusive portions associated with clearly identifiable environmental conditions. In terrestrial species, citizen science data granted the highest gain in environmental space to the pooled niches, determining an increased future biological invasion risk. A few aquatic species modelled at the regional scale reported a net loss in the pooled niches compared to their scientific survey niches, suggesting that citizen science data may also lead to contraction in pooled niches. For these species, models predicted a lower future biological invasion risk. These findings indicate that citizen science data may represent a valuable contribution to predicting future spread of invasive alien species, especially within national-scale programmes. At the same time, citizen science data collected on species poorly known to citizen scientists, or in strictly local contexts, may strongly affect the niche quantification of these taxa and the prediction of their future biological invasion risk.

Corrigendum: Relocation of bioclimatic suitability of Portuguese grapevine varieties under climate change scenarios.

[This corrects the article DOI: 10.3389/fpls.2023.974020.].

Space-for-time substitutions in climate change ecology and evolution.

In an epoch of rapid environmental change, understanding and predicting how biodiversity will respond to a changing climate is an urgent challenge. Since we seldom have sufficient long-term biological data to use the past to anticipate the future, spatial climate-biotic relationships are often used as a proxy for predicting biotic responses to climate change over time. These 'space-for-time substitutions' (SFTS) have become near ubiquitous in global change biology, but with different subfields largely developing methods in isolation. We review how climate-focussed SFTS are used in four subfields of ecology and evolution, each focussed on a different type of biotic variable - population phenotypes, population genotypes, species' distributions, and ecological communities. We then examine the similarities and differences between subfields in terms of methods, limitations and opportunities. While SFTS are used for a wide range of applications, two main approaches are applied across the four subfields: spatial in situ gradient methods and transplant experiments. We find that SFTS methods share common limitations relating to (i) the causality of identified spatial climate-biotic relationships and (ii) the transferability of these relationships, i.e. whether climate-biotic relationships observed over space are equivalent to those occurring over time. Moreover, despite widespread application of SFTS in climate change research, key assumptions remain largely untested. We highlight opportunities to enhance the robustness of SFTS by addressing key assumptions and limitations, with a particular emphasis on where approaches could be shared between the four subfields.

Unveiling the Potential Distribution of the Highly Threatened Madeira Pipistrelle (Pipistrellus maderensis): Do Different Evolutionary Significant Units Exist?

The isolation of islands has played a significant role in shaping the unique evolutionary histories of many species of flora and fauna, including bats. One notable example is the Madeira pipistrelle (Pipistrellus maderensis), which inhabits the Macaronesian archipelagos of the Azores, Madeira, and the Canary Islands. Despite the high biogeographic and conservation importance of this species, there is limited information on its ecology and evolutionary history across different archipelagos. In our study, we employed species distribution models (SDMs) to identify suitable habitats for the Madeira pipistrelle and determine the environmental factors influencing its distribution. Additionally, we conducted molecular comparisons using mitochondrial DNA data from various Macaronesian islands. Molecular analyses provided compelling evidence for the presence of distinct Evolutionary Significant Units on the different archipelagos. We identified distinct haplotypes in the populations of Madeira and the Canary Islands, with a genetic distance ranging from a minimum of 2.4% to a maximum of 3.3% between samples from different archipelagos. In support of this, SDMs highlighted relevant dissimilarities between the environmental requirements of the populations of the three archipelagos, particularly the climatic niche. Our research demonstrates that deeper investigations that combine ecological, morphological, and genetic areas are necessary to implement tailored conservation strategies.

Phylogeographic and Paleoclimatic Modelling Tools Improve Our Understanding of the Biogeographic History of Hierophis viridiflavus (Colubridae).

Phylogeographic and paleoclimatic modelling studies have been combined to infer the role of Pleistocene climatic oscillations as drivers of the genetic structure and distribution of Mediterranean taxa. For the European whip snake, Hierophis viridiflavus, previous studies based on paleoclimatic modelling have depicted a low reliability in the pattern of past climatic suitability across the central Mediterranean Basin, which barely fits the species' genetic structure. In this study, we combined phylogeographic and paleoclimatic modelling tools to improve our understanding of the biogeographic history of H. viridiflavus, particularly extending the sampling and phylogeographic inferences to previously under-sampled regions. Phylogeographic analyses recovered two major clades that diverged at the beginning of the Pleistocene and had diversified in different ways by the late Pleistocene: the east clade (composed of three subclades) and the west clade (with no further structure). Paleoclimatic models highlighted the temperate character of H. viridiflavus, indicating range contractions during both the last inter-glacial and last glacial maximum periods. Range expansions from southern-located climatic refugia likely occurred in the Bølling-Allerød and Middle Holocene periods, which are supported by signals of demographic growth in the west clade and South-East-North subclade. Overall, this work improves our understanding of the historical biogeography of H. viridiflavus, providing further insights into the evolutionary processes that occurred in the Mediterranean Basin hotspot.

Projected northward shifts in eastern red-backed salamanders due to changing climate.

Many species' distributions are being impacted by the acceleration of climate change. Amphibians in particular serve numerous ecosystem functions and are useful indicators of environmental change. Understanding how their distributions have been impacted by climate change and will continue to be impacted is thus important to overall ecosystem health. Plethodon cinereus (Eastern Red-Backed Salamander) is a widespread species of lungless salamander (Plethodontidae) that ranges across northeastern North America. To better understand future potential lungless salamander range shifts, we quantify environmental favorability, the likelihood of membership in a set of sites where environmental conditions are favorable for a species, for P. cinereus in multiple time periods, and examine shifts in the species' distribution. First, utilizing a large data set of georeferenced records, we assessed which bioclimatic variables were associated with environmental favorability in P. cinereus. We then used species distribution modeling for two time periods (1961-1980 and 2001-2020) to determine whether there was a regional shift in environmental favorability in the past 60 years. Models were then used to project future distributions under eight climate change scenarios to quantify potential range shifts. Shifts were assessed using fuzzy logic, avoiding thresholds that oversimplify model predictions into artificial binary outputs. We found that P. cinereus presence is strongly associated with environmental stability. There has been a substantial northward shift in environmental favorability for P. cinereus between 1961-1980 and 2001-2020. This shift is predicted to continue by 2070, with larger shifts under higher greenhouse gas emission scenarios. As climate change accelerates, it is differentially impacting species but has especially strong impacts on dispersal-limited species. Our results show substantial northward shifts in climatic favorability in the last 60 years for P. cinereus, which are likely to be exacerbated by ongoing climate change. Since P. cinereus is dispersal-limited, these models may imply local extirpations along the southern modern range with limited northward dispersal. Continued monitoring of amphibians in the field will reveal microclimatic effects associated with climate change and the accuracy of the model predictions presented here.

The past, present, and future of ecogeographic isolation between closely related Aquilegia plants.

Quantifying the strength of the ecogeographic barrier is an important aspect of plant speciation research, and serves as a practical step to understanding the evolutionary trajectory of plants under climate change. Here, we quantified the extent of ecogeographic isolation in four closely related Aquilegia species that radiated in the Mountains of SW China and adjacent regions, often lacking intrinsic barriers. We used environmental niche models to predict past, present, and future species potential distributions and compared them to determine the degree of overlap and ecogeographic isolation. Our investigation found significant ecological differentiation in all studied species pairs except A. kansuensis and A. ecalacarata. The current strengths of ecogeographic isolation are above 0.5 in most cases. Compared with current climates, most species had an expanding range in the Last Glacial Maximum, the Mid Holocene, and under four future climate scenarios. Our results suggested that ecogeographic isolation contributes to the diversification and maintenance of Aquilegia species in the Mountains of northern and SW China and would act as an essential reproductive barrier in the future.

Cryptosporidiosis threat under climate change in China: prediction and validation of habitat suitability and outbreak risk for human-derived Cryptosporidium based on ecological niche models.

BACKGROUND: Cryptosporidiosis is a zoonotic intestinal infectious disease caused by Cryptosporidium spp., and its transmission is highly influenced by climate factors. In the present study, the potential spatial distribution of Cryptosporidium in China was predicted based on ecological niche models for cryptosporidiosis epidemic risk warning and prevention and control. METHODS: The applicability of existing Cryptosporidium presence points in ENM analysis was investigated based on data from monitoring sites in 2011-2019. Cryptosporidium occurrence data for China and neighboring countries were extracted and used to construct the ENMs, namely Maxent, Bioclim, Domain, and Garp. Models were evaluated based on Receiver Operating Characteristic curve, Kappa, and True Skill Statistic coefficients. The best model was constructed using Cryptosporidium data and climate variables during 1986‒2010, and used to analyze the effects of climate factors on Cryptosporidium distribution. The climate variables for the period 2011‒2100 were projected to the simulation results to predict the ecological adaptability and potential distribution of Cryptosporidium in future in China. RESULTS: The Maxent model (AUC = 0.95, maximum Kappa = 0.91, maximum TSS = 1.00) fit better than the other three models and was thus considered the best ENM for predicting Cryptosporidium habitat suitability. The major suitable habitats for human-derived Cryptosporidium in China were located in some high-population density areas, especially in the middle and lower reaches of the Yangtze River, the lower reaches of the Yellow River, and the Huai and the Pearl River Basins (cloglog value of habitat suitability > 0.9). Under future climate change, non-suitable habitats for Cryptosporidium will shrink, while highly suitable habitats will expand significantly (χ2 = 76.641, P < 0.01; χ2 = 86.836, P < 0.01), and the main changes will likely be concentrated in the northeastern, southwestern, and northwestern regions. CONCLUSIONS: The Maxent model is applicable in prediction of Cryptosporidium habitat suitability and can achieve excellent simulation results. These results suggest a current high risk of transmission and significant pressure for cryptosporidiosis prevention and control in China. Against a future climate change background, Cryptosporidium may gain more suitable habitats within China. Constructing a national surveillance network could facilitate further elucidation of the epidemiological trends and transmission patterns of cryptosporidiosis, and mitigate the associated epidemic and outbreak risks.