Prediction of the potential geographic distribution of the ectomycorrhizal mushroom Tricholoma matsutake under multiple climate change scenarios

Predicting the Impact of Climate Change on the Future Distribution of Paederus fuscipes Curtis, 1826, in China Based on the MaxEnt Model

Paederus fuscipes Curtis, 1826, belongs to the Coleoptera order, Staphylinidae family, and Paederus genus (Fabricius, 1775). It has a wide distribution and strong invasive and environmental adaptation capabilities. As a predatory natural enemy of agricultural and forestry pests, understanding its suitable habitat is crucial for the control of other pests. This study, for the first time, uses the MaxEnt model and ArcGIS software, combining known distribution information of P. fuscipes and climate environmental factors to predict the current and future suitable habitat distribution of this insect. The key environmental variables affecting the distribution of P. fuscipes have been identified as mean diurnal range (mean of monthly (max temp-min temp)) (bio2), isothermality (Bio2/Bio7) (*100) (bio3), minimum temperature of the coldest month (bio6), temperature annual range (bio5-bio6) (bio7), mean temperature of the driest quarter (bio9), mean temperature of the coldest quarter (bio11), precipitation of the wettest month (bio13), precipitation of the driest month (bio14), and precipitation seasonality (coefficient of variation) (bio15). The highly suitable areas for P. fuscipes in China are mainly distributed in the hilly regions of Shandong, the North China Plain, and the middle and lower reaches of the Yangtze River Plain, with a total suitable area of 118.96 × 104 km2, accounting for 12.35% of China's total area. According to future climate change scenarios, it is predicted that the area of highly and lowly suitable regions will significantly decrease, while moderately suitable regions will increase (except for the 2090s, SSP2-4.5 scenario). These research findings provide important theoretical support for pest control and ecological conservation applications.

Prediction of potential distributions of Morina kokonorica and Morina chinensis in China

Changes in the habitats of species can provide insights into the impact of climate change on their habitats. Species in the genus Morina (Morinoideae) are perennial herbaceous plants that are mainly distributed in the South Asian Mountains and Eastern Mediterranean. In China, there are four species and two varieties of this genus distributed across the Yunnan, Sichuan, Qinghai, and Gansu provinces. This study used the optimal MaxEnt model to simulate past, current, and future potentially suitable habitats of Morina kokonorica and Morina chinensis. Seventy data of M. kokonorica occurrences and 3 of M. chinensis were used in the model to predict potentially suitable habitats. The model prediction results indicated that both M. kokonorica and M. chinensis exhibited trends of northward migration to higher latitudes and westward migration along the Himalayas to higher elevations, suggesting that the northern valleys of Hengduan Mountains and northern and eastern parts of the Himalayas were potential refugia for M. kokonorica, and the potential refugia for M. chinensis was located in the eastern part of Qinghai-Tibet Plateau. The results of this niche analysis showed that the two species had higher levels of interspecific competition and that the environmental adaptability of M. chinensis was stronger. This research could help further understand the response pattern of Morina to environmental change, to understand the adaptability of species to the environment, and promote the protection of species.

Comparative Assessment of Habitat Suitability and Niche Overlap of Three Cytospora Species in China

The plant pathogenic fungus Cytospora is notoriously known for causing woody plant canker diseases, resulting in substantial economic losses to biological forests and fruit trees worldwide. Despite their strong negative ecological impact, the existing and prospective distribution patterns of these plant pathogens in China, according to climate change, have received little attention. In this study, we chose three widely dispersed and seriously damaging species, namely, Cytospora chrysosperma, Cytospora mali, and Cytospora nivea, which are the most common species that damage the Juglans regia, Malus domestica, Eucalyptus, Pyrus sinkiangensis, Populus spp., and Salix spp. in China. We utilized ecological niche modeling to forecast their regional distribution in China under four climate change scenarios (present, SSP 126, SSP 370, and SSP 585). The results show that temperature-related climate factors limit the current distribution ranges of the three species. Currently, the three studied species are highly suitable for northeast, northwest, north, and southwest China. Under future climate scenarios, the distribution ranges of the three species are projected to increase, and the centers of the adequate distribution areas of the three species are expected to shift to high-latitude regions. The three species coexist in China, primarily in the northwest and north regions. The ecological niches of C. chrysosperma and C. nivea are more similar. The distribution range of C. mali can reach the warmer and wetter eastern region, whereas C. chrysosperma and C. nivea are primarily found in drought-prone areas with little rainfall. Our findings can help farmers and planners develop methods to avoid the spread of Cytospora spp. and calculate the costs of applying pesticides to reduce contamination and boost yields.

The ecological suitability area of Cirsium lineare (Thunb.) Sch.-Bip. under future climate change in China based on MaxEnt modeling

Many kinds of medicinal ingredients occur in Cirsium lineare that have good clinical efficacy, conferring on this species its high medicinal development value. However, with a rapidly changing global climate, it is increasingly imperative to study the factors affecting the habitat distribution and survival of species. We predicted the current and future distribution areas of suitable habitats for C. lineare, analyzed the importance of environmental variables in influencing habitat shifts, and described the alterations to suitable habitats of C. lineare in different periods (modern, 2050s, and 2070s) and scenarios (RCP2.6, RCP4.5, and RCP8.5). The results show that, under the current climate, the total suitable area of C. lineare is about 2,220,900 km2, of which the highly suitable portion amounts to ca. 292,600 km2. The minimum temperature of the coldest month, annual precipitation, and mean daily temperature range are the chief environmental variables affecting the distribution of habitat for C. lineare. In the same period, with rising greenhouse gas emission concentrations, the total suitable area will increase. In general, under future climate change, the suitable habitat for C. lineare will gradually migrate to the west and north, and its total suitable area will also expand. The results of this experiment can be used for the conservation and management of the wild resources of C. lineare. We can choose suitable growth areas to protect the medicinal resources of C. lineare through in situ conservation and artificial breeding.

Predictive Modeling Analysis for the Quality Indicators of Matsutake Mushrooms in Different Transport Environments

Matsutake mushrooms, known for their high value, present challenges due to their seasonal availability, difficulties in harvesting, and short shelf life, making it crucial to extend their post-harvest preservation period. In this study, we developed three quality predictive models of Matsutake mushrooms using three different methods. The quality changes of Matsutake mushrooms were experimentally analyzed under two cases (case A: Temperature control and sealing measures; case B: Alteration of gas composition) with various parameters including the hardness, color, odor, pH, soluble solids content (SSC), and moisture content (MC) collected as indicators of quality changes throughout the storage period. Prediction models for Matsutake mushroom quality were developed using three different methods based on the collected data: multiple linear regression (MLR), support vector regression (SVR), and an artificial neural network (ANN). The comparative results reveal that the ANN outperforms MLR and SVR as the optimal model for predicting Matsutake mushroom quality indicators. To further enhance the ANN model's performance, optimization techniques such as the Levenberg-Marquardt, Bayesian regularization, and scaled conjugate gradient backpropagation algorithm techniques were employed. The optimized ANN model achieved impressive results, with an R-Square value of 0.988 and an MSE of 0.099 under case A, and an R-Square of 0.981 and an MSE of 0.164 under case B. These findings provide valuable insights for the development of new preservation methods, contributing to the assurance of a high-quality supply of Matsutake mushrooms in the market.

A brighter shade of future climate on Himalayan musk deer Moschus leucogaster

Himalayan musk deer (Moschus leucogaster) is classified as an endangered species by IUCN with a historically misunderstood distribution due to misidentification with other species of musk deer, Moschus spp. Taking advantage of recent genetic analyses confirming the species of various populations in Nepal and China, we produced an accurate estimate of the species' current and future distribution under multiple climate change scenarios. We collected high-quality occurrence data using systematic surveys of various protected areas of Nepal to train species distribution models. The most influential determinants of the distribution of Himalayan musk deer were precipitation of the driest quarter, temperature seasonality, and annual mean temperature. These variables, and precipitation in particular, determine the vegetation type and structure in the Himalaya, which is strongly correlated with the distribution of Himalayan musk deer. We predicted suitable habitats between the Annapurna and Kanchenjunga region of Nepal Himalaya as well as the adjacent Himalaya in China. Under multiple climate change scenarios, the vast majority (85-89%) of current suitable sites are likely to remain suitable and many new areas of suitable habitat may emerge to the west and north of the current species range in Nepal and China. Two-thirds of current and one-third of future suitable habitats are protected by the extensive network of protected areas in Nepal. The projected large gains in suitable sites may lead to population expansion and conservation gains, only when the threat of overexploitation and population decline is under control.

Environmental and climatic drivers of phenotypic evolution and distribution changes in a widely distributed subfamily of subterranean mammals

How environmental factors shape species morphology and distributions is a key issue in ecology, especially in similar environments. Species of Myospalacinae exhibit widespread distribution spanning the eastern Eurasian steppe and the extreme adaptation to the subterranean environment, providing an excellent opportunity for investigating species responses to environmental changes. At the national scale, we here use geometric morphometric and distributional data to assess the environmental and climatic drivers of morphological evolution and distribution of Myospalacinae species in China. Based on phylogenetic relationships of Myospalacinae species constructed using genomic data in China, we integrate geometric morphometrics and ecological niche models to reveal the interspecific variation of skull morphology, trace the ancestral state, and assess factors influencing interspecific variation. Our approach further allows us to project future distributions of Myospalacinae species throughout China. We found that the interspecific morphology variations were mainly concentrated in the temporal ridge, premaxillary-frontal suture, premaxillary-maxillary suture, and molars, and the skull morphology of the two current species in Myospalacinae followed the ancestral state; temperature and precipitation were important environmental variables influencing skull morphology. Elevation, temperature annual range, and precipitation of warmest quarter were identified as dominant factors affecting the distribution of Myospalacinae species in China, and their suitable habitat area will decrease in the future. Collectively, environmental and climate changes have an effect on skull phenotypes of subterranean mammals, highlighting the contribution of phenotypic differentiation in similar environments in the formation of species phenotypes. Climate change will further shrink their habitats under future climate assumptions in the short-term. Our findings provide new insights into effects of environmental and climate change on the morphological evolution and distribution of species as well as a reference for biodiversity conservation and species management.

Invasive Trends of Spartina alterniflora in the Southeastern Coast of China and Potential Distributional Impacts on Mangrove Forests

Mangrove forests are one of the most productive and seriously threatened ecosystems in the world. The widespread invasion of Spartina alterniflora has seriously imperiled the security of mangroves as well as coastal mudflat ecosystems. Based on a model evaluation index, we selected RF, GBM, and GLM as a predictive model for building a high-precision ensemble model. We used the species occurrence records combined with bioclimate, sea-land topography, and marine environmental factors to predict the potentially suitable habitats of mangrove forests and the potentially suitable invasive habitats of S. alterniflora in the southeastern coast of China. We then applied the invasion risk index (IRI) to assess the risk that S. alterniflora would invade mangrove forests. The results show that the suitable habitats for mangrove forests are mainly distributed along the coastal provinces of Guangdong, Hainan, and the eastern coast of Guangxi. The suitable invasive habitats for S. alterniflora are mainly distributed along the coast of Zhejiang, Fujian, and relatively less in the southern provinces. The high-risk areas for S. alterniflora invasion of mangrove forests are concentrated in Zhejiang and Fujian. Bioclimate variables are the most important variables affecting the survival and distribution of mangrove forests and S. alterniflora. Among them, temperature is the most important environmental variable determining the large-scale distribution of mangrove forests. Meanwhile, S. alterniflora is more sensitive to precipitation than temperature. Our results can provide scientific insights and references for mangrove forest conservation and control of S. alterniflora.

Comparison between optimized MaxEnt and random forest modeling in predicting potential distribution: A case study with Quasipaa boulengeri in China

Random forest (RF) and MaxEnt models are shallow machine learning approaches that perform well in predicting species' potential distributions. RF models can produce robust results with the default automatic configuration in most cases, but it is necessary for MaxEnt to optimize the model settings to improve the performance, and the predictive performance difference between optimized MaxEnt and RF is uncertain. To explore this issue, the potential distribution of the endangered amphibian Quasipaa boulengeri in China was predicted using optimized MaxEnt and RF models. A total of 408 occurrence data were selected, 1000 locations were generated as pseudo-absence data by the geographic distance method, and 10,000 sites were selected as background data by creating a bias file. Partial ROC at different thresholds and success rate curves were used to compare the predictive performances between optimized MaxEnt and RF. Our results showed that the RF and optimized MaxEnt models both had good performance in predicting the potential distribution of Q. boulengeri, with the RF model performing slightly better whether based on partial ROC or success rate curves. Furthermore, the core suitable habitat regions of Q. boulengeri identified by RF and MaxEnt were similar and were all located in the Sichuan, Chongqing, Hubei, Hunan, and Guizhou provinces. However, the RF model produced a habitat suitability map with higher discrimination and greater heterogeneity. Temperature annual range, mean temperature of the driest quarter, and annual precipitation were the vital environmental variables limiting the distribution of Q. boulengeri. The RF model is the stronger machine learner. We believe it may be more applicable in predicting the native potential distributions of species with sufficient occurrence data, given the additional predictive detail, the simplicity of use, the computational time involved, and the operational complexity.

Response of distribution patterns of two closely related species in Taxus genus to climate change since last inter-glacial

Climate change affects the species spatio-temporal distribution deeply. However, how climate affects the spatio-temporal distribution pattern of related species on the large scale remains largely unclear. Here, we selected two closely related species in Taxus genus Taxus chinensis and Taxus mairei to explore their distribution pattern. Four environmental variables were employed to simulate the distribution patterns using the optimized Maxent model. The results showed that the highly suitable area of T. chinensis and T. mairei in current period was 1.616 × 10⁵ km² and 3.093 × 10⁵ km², respectively. The distribution area of T. chinensis was smaller than that of T. mairei in different periods. Comparison of different periods shown that the distribution area of the two species was almost in stasis from LIG to the future periods. Temperature and precipitation were the main climate factors that determined the potential distribution of the two species. The centroids of T. chinensis and T. mairei were in Sichuan and Hunan provinces in current period, respectively. In the future, the centroid migration direction of the two species would shift towards northeast. Our results revealed that the average elevation distribution of T. chinensis was higher than that of T. mairei. This study sheds new insights into the habitat preference and limiting environment factors of the two related species and provides a valuable reference for the conservation of these two threatened species.

Potential distribution of Blumea balsamifera in China using MaxEnt and the ex situ conservation based on its effective components and fresh leaf yield

Blumea balsamifera is a famous Chinese Minority Medicine, which has a long history in Miao, Li, Zhuang, and other minority areas. In recent years, due to the influence of natural and human factors, the distribution area of B. balsamifera resources has a decreasing trend. Therefore, it is very important to analyze the suitability of B. balsamifera in China. Following three climate change scenarios (SSP1-2.6, SSP2-4.5, and SSP5-8.5) under 2050s and 2070s, geographic information technology (GIS) and maximum entropy model (MaxEnt) were used to simulate the ecological suitability of B. balsamifera. The contents of L-borneol and total flavonoids of B. balsamifera in different populations were determined by gas chromatography (GC) and ultraviolet spectrophotometry (UV). The results showed that the key environmental variables affecting the distribution of B. balsamifera were mean temperature of coldest quarter (6.18-26.57 ℃), precipitation of driest quarter (22.46-169.7 mm), annual precipitation (518.36-1845.29 mm), and temperature seasonality (291.31-878.87). Under current climate situation, the highly suitable habitat was mainly located western Guangxi, southern Yunnan, most of Hainan, southwestern Guizhou, southwestern Guangdong, southeastern Fujian, and western Taiwan, with a total area of 24.1 × 10⁴ km². The areas of the moderately and poorly suitable habitats were 27.57 × 10⁴ km² and 42.43 × 10⁴ km², respectively. Under the future climate change scenarios, the areas of the highly, moderately, and poorly suitable habitats of B. balsamifera showed a significant increasing trend, the geometric center of the total suitable habitats of B. balsamifera would move to the northeast. In recent years, the planting area of B. balsamifera has been reduced on a large scale in Guizhou, and its ex situ protection is imperative. By comparison, the content of L-borneol, total flavonoids and fresh leaf yield had no significant difference between Guizhou and Hainan (P > 0.05), which indicated that Hainan is one of the best choice for ex situ protection of B. balsamifera.

Effects of climate change on the distribution of wild Akebia trifoliata

Understanding the impacts and constraints of climate change on the geographical distribution of wild Akebia trifoliata is crucial for its sustainable management and economic development as a medicinal material or fruit. In this study, according to the first‐hand information obtained from field investigation, the distribution and response to climate change of A. trifoliata were studied by the MaxEnt model and ArcGIS. The genetic diversity and population structure of 21 natural populations of A. trifoliata were studied by simple sequence repeat (SSR) markers. The results showed that the most important bioclimatic variable limiting the distribution of A. trifoliata was the Mean Temperature of Coldest Quarter (bio11). Under the scenarios SSP1‐2.6 and SSP2‐4.5, the suitable area of A. trifoliata in the world will remain stable, and the suitable area will increase significantly under the scenarios of SSP3‐7.0 and SSP5‐8.5. Under the current climate scenario, the suitable growth regions of A. trifoliata in China were 79.9–122.7°E and 21.5–37.5°N. Under the four emission scenarios in the future, the geometric center of the suitable distribution regions of Akebia trifoliata in China will move to the north. The clustering results of 21 populations of A. trifoliata analyzed by SSR markers showed that they had a trend of evolution from south to north.

Predicting the potential impacts of climate change on the endangered Caesalpinia bonduc (L.) Roxb in Benin (West Africa)

Caesalpinia bonduc (L.) Roxb is a medicinal plant with high therapeutic values but declared extinct in the wild in Benin. This study explored the potential distribution and climatic suitability of the species under the present-day and future conditions in Benin, based on two Representative Concentration Pathways (RCP4.5 and RCP8.5) at the 2055-time horizon. The occurrence data were recorded in the distribution area of the species in Benin. These data were supplemented with those from the Global Biodiversity Information Facility (GBIF, www.gbif.org) website and the literature. A total of 23 environmental variables (15 bioclimatic data and 8 biophysical data) were used. The Bioclimatic variables for temperature and humidity were downloaded from Africlim site at 1 km resolution. The biophysical variables concern population, elevation, slope, landcover, wetland, distance to river, soil and distance to dwellings data that are downloaded respectively from DIVA-GIS, ISRIC and SEDAC website at different resolution. A correlation test has been applied to eliminate the highly correlated variables (r ≥ 0.9) using Pearson correlation coefficient. Species distribution modelling data were processed using five algorithms namely Random Forest (RF), Boosted Regression Trees (BRT), Maximum entropy (MAXENT), Generalized Linear Models (GLM) and Generalized Additive Models (GAM). The results showed that all models performed well with the area under the curve (AUC) values greater than 0.9. The RF, GLM, and GAM models predicted an increase in the suitable areas for the cultivation of the species. BRT and MaxEnt showed a substantial decrease in the suitable areas based on the two scenarios but this reduction is more observed with the MaxEnt model. These results show that climate change and human pressures will have significant effects on the distribution of C. bonduc throughout Benin. Sustainable management measures are necessary for C. bonduc and should be integrated in development policies to preserve the population of the species from total extinction in Benin.

Climate change threatens native potential agroforestry plant species in Brazil

Climate change is one of the main drivers of species extinction in the twentyfirst-century. Here, we  (1) quantify potential changes in species' bioclimatic area of habitat (BAH) of 135 native potential agroforestry species from the Brazilian flora, using two different climate change scenarios (SSP2-4.5 and SSP5-8.5) and dispersal scenarios, where species have no ability to disperse and reach new areas (non-dispersal) and where species can migrate within the estimated BAH (full dispersal) for 2041–2060 and 2061–2080. We then (2) assess the preliminary conservation status of each species based on IUCN criteria. Current and future potential habitats for species were predicted using MaxEnt, a machine-learning algorithm used to estimate species' probability distribution. Future climate is predicted to trigger a mean decline in BAH between 38.5–56.3% under the non-dispersal scenario and between 22.3–41.9% under the full dispersal scenario for 135 native potential agroforestry species. Additionally, we found that only 4.3% of the studied species could be threatened under the IUCN Red List criteria B1 and B2. However, when considering the predicted quantitative habitat loss due to climate change (A3c criterion) the percentages increased between 68.8–84.4% under the non-dispersal scenario and between 40.7–64.4% under the full dispersal scenario. To lessen such threats, we argue that encouraging the use of these species in rural and peri-urban agroecosystems are promising, complementary strategies for their long-term conservation.

Determination of Diversity, Distribution and Host Specificity of Korean Laccaria Using Four Approaches

The genus Laccaria (Hydnangiaceae, Agaricales) plays an important role in forest ecosystems as an ectomycorrhizal fungus, contributing to nutrient cycles through symbiosis with many types of trees. Though understanding Laccaria diversity and distribution patterns, as well as its association with host plants, is fundamental to constructing a balanced plant diversity and conducting effective forest management, previous studies have not been effective in accurately investigating, as they relied heavily on specimen collection alone. To investigate the true diversity and distribution pattern of Laccaria species and determine their host types, we used four different approaches: specimen-based analysis, open database search (ODS), NGS analysis, and species-specific PCR (SSP). As a result, 14 Laccaria species have been confirmed in Korea. Results regarding the species distribution pattern were different between specimen-based analysis and SSP. However, when both were integrated, the exact distribution pattern of each Laccaria species was determined. In addition, the SSP revealed that many Laccaria species have a wide range of host types. This study shows that using these four different approaches is useful in determining the diversity, distribution, and host of ECM fungi. Furthermore, results obtained for Laccaria will serve as a baseline to help understand the role of ECM fungi in forest management in response to climate change.

Metagenomic Analysis of Bacterial and Fungal Communities Inhabiting Shiro Dominant Soils of Two Production Regions of Tricholoma Matsutake S. Ito & S. Imai in Korea

Tricholoma matsutake is an ectomycorrhizal fungus that has obligate symbiotic relationships with Pinus densiflora. Its fruiting body has a distinctive flavor and is traded at a high price. Thus, it has been a significant source of income for rural communities in Korea. We hypothesized that biotic factors considerably influence the formation of the T. matsutake mushroom, and the soils producing T. matsutake share similar microbial characteristics. Therefore, the present study aimed to detect the specific fungal and bacterial groups in T. matsutake production soils (shiro+) and nonproduction soils (shiro−) of the Bonghwa and Yanyang regions via next-generation sequencing. In a total of 15 phyla, 36 classes, 234 genera of bacteria, six phyla, 29 classes, and 164 genera of fungi were detected from four samples at both sites. The species diversity of shiro+ soils was lower than the shiro− samples in both the fungal and bacterial groups. In addition, we did not find high similarities in the microbial communities between the shiro+ soils of the two regions. However, in the resulting differences between the fungal communities categorized by their trophic assembly, we found a distinguishable compositional pattern in the fungal communities from the shiro+ soils and the shiro− soils of the two sites. Thus, the similarity among the microbial communities in the forest soils may be due to the fact that the microbial communities in the T. matsutake dominant soils are closely associated with biotic factors and abiotic factors such as soil properties.

Spatiotemporal dynamics of habitat suitability for the Ethiopian staple crop, Eragrostis tef (teff), under changing climate

Teff (Eragrostis tef (Zucc.) Trotter) is a staple, ancient food crop in Ethiopia. Its growth is affected by climate change, so it is essential to understand climatic effects on its habitat suitability in order to design countermeasures to ensure food security. Based on the four Representative Concentration Pathway emission scenarios (i.e., RCP2.6, RCP4.5, RCP6.0 and RCP8.5) set by the Intergovernmental Panel on Climate Change (IPCC), we predicted the potential distribution of teff under current and future scenarios using a maximum entropy model (Maxent). Eleven variables were selected out of 19, according to correlation analysis combined with their contribution rates to the distribution. Simulated accuracy results validated by the area under the curve (AUC) had strong predictability with values of 0.83–0.85 for current and RCP scenarios. Our results demonstrated that mean temperature in the coldest season, precipitation seasonality, precipitation in the cold season and slope are the dominant factors driving potential teff distribution. Proportions of suitable teff area, relative to the total study area were 58% in current climate condition, 58.8% in RCP2.6, 57.6% in RCP4.5, 59.2% in RCP6.0, and 57.4% in RCP8.5, respectively. We found that warmer conditions are correlated with decreased land suitability. As expected, bioclimatic variables related to temperature and precipitation were the best predictors for teff suitability. Additionally, there were geographic shifts in land suitability, which need to be accounted for when assessing overall susceptibility to climate change. The ability to adapt to climate change will be critical for Ethiopia’s agricultural strategy and food security. A robust climate model is necessary for developing primary adaptive strategies and policy to minimize the harmful impact of climate change on teff.

Quo vadis? Historical distribution and impact of climate change on the worldwide distribution of the Australasian fungus Clathrus archeri (Phallales, Basidiomycota)

Clathrus archeri is a fungus native to Australia and New Zealand that has started to expand into Europe, and it is considered a potentially invasive species. In this study, we examine the historical occurrence, current geographical range and potential future changes in the distribution of C. archeri using worldwide distribution data. Ecological modelling was used to assess the locations of the potential climatic niches of C. archeri within both its native and introduced ranges in the past, present and future. Our study clearly shows that the coverage of suitable habitats of this fungus has decreased since the last glacial maximum, and anthropogenic climate changes are accelerating the process of niche loss. The highest rate of C. archeri range contraction is expected in Australia, where the fungus should be considered a threatened species in the future. Highly valuable habitats will be available in Tasmania and New Zealand. However, a significant expansion rate of C. archeri will still probably be observed in Europe, where the climatic conditions preferred by the fungus will allow its quick expansion northeast into the continent.

Impact of Past and Future Climate Change on the Potential Distribution of an Endangered Montane Shrub Lonicera oblata and Its Conservation Implications

Climate change is an important driver of biodiversity patterns and species distributions, understanding how organisms respond to climate change will shed light on the conservation of endangered species. In this study, we modeled the distributional dynamics of a critically endangered montane shrub Lonicera oblata in response to climate change under different periods by building a comprehensive habitat suitability model considering the effects of soil and vegetation conditions. Our results indicated that the current suitable habitats for L. oblata are located scarcely in North China. Historical modeling indicated that L. oblata achieved its maximum potential distribution in the last interglacial period which covered southwest China, while its distribution area decreased for almost 50% during the last glacial maximum. It further contracted during the middle Holocene to a distribution resembling the current pattern. Future modeling showed that the suitable habitats of L. oblata contracted dramatically, and populations were fragmentedly distributed in these areas. As a whole, the distribution of L. oblata showed significant migration northward in latitude but no altitudinal shift. Several mountains in North China may provide future stable climatic areas for L. oblata, particularly, the intersections between the Taihang and Yan mountains. Our study strongly suggested that the endangered montane shrub L. oblata are sensitive to climate change, and the results provide new insights into the conservation of it and other endangered species.

Predicted climate change will increase the truffle cultivation potential in central Europe

Climate change affects the distribution of many species, including Burgundy and Périgord truffles in central and southern Europe, respectively. The cultivation potential of these high-prized cash crops under future warming, however, remains highly uncertain. Here we perform a literature review to define the ecological requirements for the growth of both truffle species. This information is used to develop niche models, and to estimate their cultivation potential in the Czech Republic under current (2020) and future (2050) climate conditions. The Burgundy truffle is already highly suitable for cultivation on ~ 14% of agricultural land in the Czech Republic (8486 km²), whereas only ~ 8% of the warmest part of southern Moravia are currently characterised by a low suitability for Périgord truffles (6418 km²). Though rising temperatures under RCP8.5 will reduce the highly suitable cultivation areas by 7%, the 250 km² (3%) expansion under low-emission scenarios will stimulate Burgundy truffles to benefit from future warming. Doubling the moderate and expanding the highly suitable land by 352 km² in 2050, the overall cultivation potential for Périgord truffles will rise substantially. Our findings suggest that Burgundy and Périgord truffles could become important high-value crops for many regions in central Europe with alkaline soils. Although associated with uncertainty, long-term investments in truffle cultivation could generate a wide range of ecological and economic benefits.

Past, current, and potential future distributions of unique genetic diversity in a cold-adapted mountain butterfly

AIM

Climatic changes throughout the Pleistocene have strongly modified species distributions. We examine how these range shifts have affected the genetic diversity of a montane butterfly species and whether the genetic diversity in the extant populations is threatened by future climate change.

LOCATION

Europe.

TAXON

Erebia epiphron Lepidoptera: Nymphalidae.

METHODS

We analyzed mtDNA to map current genetic diversity and differentiation of E. epiphron across Europe to identify population refugia and postglacial range shifts. We used species distribution modeling (SDM) to hindcast distributions over the last 21,000 years to identify source locations of extant populations and to project distributions into the future (2070) to predict potential losses in genetic diversity.

RESULTS

We found substantial genetic diversity unique to specific regions within Europe (total number of haplotypes = 31, number of unique haplotypes = 27, H d = 0.9). Genetic data and SDM hindcasting suggest long-term separation and survival of discrete populations. Particularly, high rates of unique diversity in postglacially colonized sites in England (H d = 0.64) suggest this population was colonized from a now extinct cryptic refugium. Under future climate change, SDMs predict loss of climate suitability for E. epiphron, particularly at lower elevations (<1,000 meters above sea level) equating to 1 to 12 unique haplotypes being at risk under climate scenarios projecting 1°C and 2-3°C increases respectfully in global temperature by 2070.

MAIN CONCLUSIONS

Our results suggest that historical range expansion and retraction processes by a cold-adapted mountain species caused diversification between populations, resulting in unique genetic diversity which may be at risk if distributions of cold-adapted species shrink in future. Assisted colonizations of individuals from at-risk populations into climatically suitable unoccupied habitat might help conserve unique genetic diversity, and translocations into remaining populations might increase their genetic diversity and hence their ability to adapt to future climate change.

The Global Potential Distribution of Invasive Plants: Anredera cordifolia under Climate Change and Human Activity Based on Random Forest Models

The potential distribution of the invasive plant Anredera cordifolia (Tenore) Steenis was predicted by Random Forest models under current and future climate-change pathways (i.e., RCP4.5 and RCP8.5 of 2050s and the 2070s). Pearson correlations were used to select variables; the prediction accuracy of the models was evaluated by using AUC, Kappa, and TSS. The results show that suitable future distribution areas are mainly in Southeast Asia, Eastern Oceania, a few parts of Eastern Africa, Southern North America, and Eastern South America. Temperature is the key climatic factor affecting the distribution of A. cordifolia. Important metrics include mean temperature of the coldest quarter (0.3 °C ≤ Bio11 ≤ 22.9 °C), max temperature of the warmest month (17.1 °C ≤ Bio5 ≤ 35.5 °C), temperature annual range (10.7 °C ≤ Bio7 ≤ 33 °C), annual mean air temperature (6.8 °C ≤ Bio1 ≤ 24.4 °C), and min temperature of coldest month (−2.8 °C ≤ Bio6 ≤ 17.2 °C). Only one precipitation index (Bio19) was important, precipitation of coldest quarter (7 mm ≤ Bio19 ≤ 631 mm). In addition, areas with strong human activities are most prone to invasion. This species is native to Brazil, but has been introduced in Asia, where it is widely planted and has escaped from cultivation. Under the future climate scenarios, suitable habitat areas of A. cordifolia will expand to higher latitudes. This study can provide a reference for the rational management and control of A. cordifolia.

Projected distribution and climate refugia of endangered Kashmir musk deer Moschus cupreus in greater Himalaya, South Asia

Kashmir musk deer Moschus cupreus (KMD) are the least studied species of musk deer. We compiled genetically validated occurrence records of KMD to construct species distribution models using Maximum Entropy. We show that the distribution of KMD is limited between central Nepal on the east and north-east Afghanistan on the west and is primarily determined by precipitation of driest quarter, annual mean temperature, water vapor, and precipitation during the coldest quarter. Precipitation being the most influential determinant of distribution suggests the importance of pre-monsoon moisture for growth of the dominant vegetation, Himalayan birch Betula utilis and Himalayan fir Abies spectabilis, in KMD's preferred forests. All four Representative Concentration Pathway Scenarios result an expansion of suitable habitat in Uttarakhand, India, west Nepal and their associated areas in China in 2050s and 2070s but a dramatic loss of suitable habitat elsewhere (Kashmir region and Pakistan-Afghanistan border). About 1/4th of the current habitat will remain as climate refugia in future. Since the existing network of protected areas will only include a tiny fraction (4%) of the climatic refugia of KMD, the fate of the species will be determined by the interplay of more urgent short-term forces of poaching and habitat degradation and long-term forces of climate change.

Non-Pessimistic Predictions of the Distributions and Suitability of Metasequoia glyptostroboides under Climate Change Using a Random Forest Model

Metasequoia glyptostroboides Hu & W. C. Cheng, which is a remarkable rare relict plant, has gradually been reduced to its current narrow range due to climate change. Understanding the comprehensive distribution of M. glyptostroboides under climate change on a large spatio-temporal scale is of great significance for determining its forest adaptation. In this study, based on 394 occurrence data and 10 bioclimatic variables, the global potential distribution of M. glyptostroboides under eight different climate scenarios (i.e., the past three, the current one, and the next four) from the Quaternary glacial to the future was simulated by a random forest model built with the biomod2 package. The key bioclimatic variables affecting the distribution of M. glyptostroboides are BIO2 (mean diurnal range), BIO1 (annual mean temperature), BIO9 (mean temperature of driest quarter), BIO6 (min temperature of coldest month), and BIO18 (precipitation of warmest quarter). The result indicates that the temperature affects the potential distribution of M. glyptostroboides more than the precipitation. A visualization of the results revealed that the current relatively suitable habitats of M. glyptostroboides are mainly distributed in East Asia and Western Europe, with a total area of approximately 6.857 × 106 km2. With the intensification of global warming in the future, the potential distribution and the suitability of M. glyptostroboides have a relatively non-pessimistic trend. Whether under the mild (RCP4.5) and higher (RCP8.5) emission scenarios, the total area of suitable habitats will be wider than it is now by the 2070s, and the habitat suitability will increase to varying degrees within a wide spatial range. After speculating on the potential distribution of M. glyptostroboides in the past, the glacial refugia of M. glyptostroboides were inferred, and projections regarding the future conditions of these places are expected to be optimistic. In order to better protect the species, the locations of its priority protected areas and key protected areas, mainly in Western Europe and East Asia, were further identified. Our results will provide theoretical reference for the long-term management of M. glyptostroboides, and can be used as background information for the restoration of other endangered species in the future.

Spatiotemporal evolution and impacts of climate change on bamboo distribution in China

Understanding the impact and restriction of climate change on potential distribution of bamboo forest is crucial for sustainable management of bamboo forest and bamboo-based economic development. In this study, climatic variables and maximum entropy model were used to simulate the potential distribution of bamboo forest in China under the future climate scenarios. Seven climatic variables, such as Spring precipitation, Summer precipitation, Autumn precipitation, average annual relative humidity, Autumn average temperature, average annual temperature range and annual total radiation, were selected as input variables of maximum entropy model based on the relative importance of those climate variables for predicting bamboo forest presence. The suitable ranges of the seven climatic variables for potential distribution of bamboo forest were 337-794 mm, 496-705 mm, 213-929 mm, 74.3%-83.4%, 16.6-23.8 °C, 2.3-10.1 °C and 3.2 × 10⁴-4.3 × 10⁴ W m^-2, respectively. Under RCP4.5 and RCP8.5 climate scenarios, the suitable area of bamboo forest growth first increased and then decreased, and showed range contractions towards the interior and expansions towards southwest in China. The results of the present study can serve as a useful reference to dynamic monitoring of the spatial distribution and sustainable utilization of bamboo forest in the future under climate change.

Assessing Habitat Suitability of Parasitic Plant Cistanche deserticola in Northwest China under Future Climate Scenarios

: Cistanche deserticola Ma, a perennial parasitic herb of family Orobanchaceae, is mainly parasitic on the roots of the Haloxylon ammodendron Bunge. In view of this special parasitic relationship, we applied random forest (RF) model to forecast potential geographic distribution, and developed a comprehensive habitat suitability model by integrating bioclimatic and soil factors to assess the suitable distribution of C. deserticola and H. ammodendron across China in 2050s and 2070s under RCP2.6, RCP4.5, and RCP8.5, respectively. We modeled the core potential geographic distribution of C. deserticola by overlaying the distribution of these two species, and analyzed the spatial distribution pattern and migration trend of C. deserticola by using the standard deviational ellipse. In addition, we evaluated the accuracy of RF model through three evaluation indexes, and analyzed the dominant climate factors. The results showed that the core potential distribution areas of C. deserticola are distributed in the Xinjiang Uygur Autonomous Region, the junction of Shaanxi–Gansu–Ningxia provinces, and the Inner Mongolia Autonomous Region. The spatial dispersion would intensify with the increasing of emission scenarios, and the geographical habitat is moving towards higher latitude. Among the three evaluation indexes, the area under the ROC curve (AUC) and True Skill Statistic (TSS) have better assessment results. The main bioclimatic factors affecting the distribution are min temperature of coldest month (Bio6), annual precipitation (Bio12), precipitation of wettest month (Bio13), precipitation of wettest quarter (Bio16), and precipitation of warmest quarter (Bio18), among which the importance of precipitation factors is greater than temperature factors. More importantly, the results of this study could provide some guidance for the improvement of desert forest system, the protection of endangered species and the further improvement of the ecological environment.

Future climate change likely to reduce the Australian plague locust (Chortoicetes terminifera) seasonal outbreaks

Climate is a major limiting factor for insect distributions and it is expected that a changing climate will likely alter spatial patterns of pest outbreaks. The Australian plague locust (APL) Chortoicetes terminifera, is the most economically important locust species in Australia. Invasions cause large scale economic damage to agricultural crops and pastures. Understanding the regional-scale and long-term dynamics is a prerequisite to develop effective control and preventive management strategies. In this study, we used a 32-year locust survey database to uncover the relationship between historical bioclimatic variables and spatial seasonal outbreaks by developing two machine learning species distribution models (SDMs), random forest and boosted regression trees. The explanatory variables were ranked by contribution to the generated models. The bio-climate models were then projected into a future climate change scenario (RCP8.5) using downscaled 34 global climate models (GCMs) to assess how climate change may alter APL seasonal distribution patterns in eastern Australia. Our results show that the model for the distribution of spring outbreaks performed better than those for summer and autumn, based on statistical evaluation criteria. The spatial models of seasonal outbreaks indicate that the areas subject to APL outbreaks were likely to decrease in all seasons. Multi-GCM ensemble means show the largest decrease in area was for spring outbreaks, reduced by 93-94% by 2071-2090, while the area of summer outbreaks decreased by 78-90%, and 67-74% for autumn outbreaks. The bioclimatic variables could explain 78-98% outbreak areas change. This study represents an important step toward the assessment of the effects of the changing climate on locust outbreaks and can help inform future priorities for regional mitigation efforts in the context of global climate change in eastern Australia.

Predicting the impacts of climate change, soils and vegetation types on the geographic distribution of Polyporus umbellatus in China

Polyporus umbellatus is a fungus that has been used medically as a diuretic for thousands of years in China. To evaluate the impacts of climatic change on the distribution of P. umbellatus, we selected the annual mean air temperature, isothermality, minimum temperature of the coldest month, annual temperature range, annual precipitation and precipitation seasonality and used observations from the 2000s and simulated values from two future periods (2041 to 2060 and 2061 to 2080) to build an ensemble model (EM); then, we developed a comprehensive habitat suitability model by integrating soil and vegetation conditions into the EM to assess the distribution of suitable P. umbellatus habitats across China in the 2000s and the two future periods. Our results show that annual precipitation and annual mean air temperature together largely determine the distribution of P. umbellatus and those suitable P. umbellatus habitats generally occur in areas with an optimal annual precipitation of approximately 1000 mm and an optimal annual mean air temperature of approximately 13 °C. In other words, P. umbellatus requires a humid and cool environment for growth. In addition, brown soils with a granular structure and low acidity are more suitable for P. umbellatus. Furthermore, we have observed that the distribution of P. umbellatus is usually associated with the presence of coniferous, mixed coniferous, and broad-leaved forests, suggesting that these vegetation types are suitable habitats for P. umbellatus. In the future, annual precipitation and annual mean air temperature will continue to increase, consequently increasing the availability of habitats suitable for P. umbellatus in northeastern and southwestern China but likely leading to a degradation of suitable P. umbellatus habitats in central China.

Modeling the Effect of Climate Change on the Potential Distribution of Qinghai Spruce (Picea crassifolia Kom.) in Qilian Mountains

Qinghai spruce forests play a key role in water conservation in the dry region of northwest China. So, it is necessary to understand the impacts of climate change on the species to implement adaptation strategies. Based on the four-emission scenario (i.e., RCP2.6 (Representative Concentration Pathway), RCP4.5, RCP6.0 and RCP8.5) set by the Intergovernmental Panel on Climate Change (IPCC) fifth assessment report, in the study, we predicted the potential distribution of Qinghai spruce (Picea crassifolia Kom.) under current and future scenarios using a maximum entropy (Maxent) model. Seven variables, selected from 22 variables according to correlation analysis combining with their contribution rates to the distribution, are used to simulate the potential distribution of the species under current and future scenarios. Simulated results are validated by area under the operating characteristic curve (AUC). Results demonstrate that elevation, mean temperature of wettest quarter, annual mean temperature, and mean diurnal range are more important in dominating the potential distribution of Qinghai spruce. Ratios of the suitable area to the total study area are 34.3% in current climate condition, 34% in RCP2.6, 33.9% in RCP4.5, 33.8% in RCP6.0, and 30.5% in RCP8.5, respectively. The warmer the climate condition is, the more area of higher suitable classification is changed to that of lower suitable classification. The ratios of real distribution area in simulated unsuitable class to the real distribution area change from 4.3% (60.7 km2) in the current climate to 13% (185 km2) in RCP8.5, suggesting that the real distribution area may decrease in the future. We conclude that there is a negative effect of climate change on the distribution of Qinghai spruce forest. The result can help decision-makers to draft adaptation countermeasures based on climate change.

Global ecological regionalization of 15 Illicium species: nature sources of shikimic acid

BACKGROUND

Illicium plants are relevant officinal and ornamental species that are native in Eastern Asia, and they are the main sources of shikimic acid. Shikimic acid is an important component of Tamiflu, which is recognized for its ability to resist avian influenza by the World Health Organization. To determine areas where 15 Illicium species can be grown and to understand the importance of species diversity, we should enhance the prediction of suitable areas.

METHODS

In this study, the global potential distribution of 15 Illicium species was predicted using a geographic information system for global medicinal plants.

RESULTS

Results showed that the possible suitable areas for these plants in China covered 1357.68 × 104 km2 (56%), and the second-largest area spanning 527.42 × 104 km2 was found in the United States. Illicium verum Hook, an edible species with the highest shikimic acid content among them, grew in areas of 59.92 × 104 (48%), 64.04 × 104 (19%), and 60.53 × 104 km2(18%) in China, the United States, and Brazil, respectively. Illicium.difengpi B. N. Chamg, an endangered species, was distributed in an area of 19.03 × 104 km2 or 95% of the total area in China.

CONCLUSIONS

This research provided a guarantee for the demand of Tamiflu, presented strategies that helped protect endangered species, and provided a reference for species cultivation and introduction.

Root-associated bacteria influencing mycelial growth of Tricholoma matsutake (pine mushroom)

Tricholoma matsutake is an ectomycorrhizal fungus usually associated with Pinus densiflora in South Korea. Fruiting bodies (mushrooms) of T. matsutake are economically important due to their attractive aroma; yet, T. matsutake is uncultivatable and its habitat is rapidly being eradicated due to global climate change. Root-associated bacteria can influence the growth of ectomycorrhizal fungi that co-exist in the host rhizosphere and distinctive bacterial communities are associated with T. matsutake. In this study, we investigated how these bacterial communities affect T. matsutake growth by isolating bacteria from the roots of P. densiflora colonized by ectomycorrhizae of T. matsutake and co-culturing rootassociated bacteria with T. matsutake isolates. Thirteen species of bacteria (27 isolates) were found in pine roots, all belonging to the orders Bacillales or Burkholderiales. Two species in the genus Paenibacillus promoted the growth of T. matsutake in glucose poor conditions, likely using soluble metabolites. In contrast, other bacteria suppressed the growth of T. matsutake using both soluble and volatile metabolites. Antifungal activity was more frequent in glucose poor conditions. In general, pine rhizospheres harbored many bacteria that had a negative impact on T. matsutake growth and the few Paenibacillus species that promoted T. matsutake growth. Paenibacillus species, therefore, may represent a promising resource toward successful cultivation of T. matsutake.

Predicting Shifts in the Suitable Climatic Distribution of Walnut (Juglans regia L.) in China: Maximum Entropy Model Paves the Way to Forest Management

Cultivation of woody oil plants in environmentally suitable habitats is a successful ecological solution for oil development and forest management. In this study, we predicted the influences of future climate change on the potentially suitable climatic distribution of an important woody oil plant species (walnut; Juglans regia L.) in China based on given climate change scenarios and the maximum entropy (MaxEnt) model. The MaxEnt model showed that the minimum temperature of the coldest month and annual precipitation were the most important determinant variables limiting the geographical distribution of J. regia. We have found that the current suitable environmental habitat of J. regia is mainly distributed in central and southwestern China. Results of the MaxEnt model showed that global warming in the coming half-century may lead to an increase in the area size of environmentally suitable habitats for J. regia in China, indicating more lands available for artificial cultivation and oil production. However, those suitable habitat gains may be practically inaccessible due to over-harvest and urban development, and effective management strategies are urgently needed to establish those forests. This research will provide theoretical suggestions for the protection, cultivation management, and sustainable utilization of J. regia resources to face the challenge of global climate change.

Effect of fruiting body bacteria on the growth of Tricholoma matsutake and its related molds

Tricholoma matsutake (pine mushroom, PM) is a prized mushroom in Asia due to its unique flavor and pine aroma. The fruiting body of PM forms only in its natural habitat (pine forest), and little is known regarding the natural conditions required for successful generation of the fruiting bodies in this species. Recent studies suggest that microbial interactions may be associated with the growth of PM; however, there have been few studies of the bacterial effects on PM growth. In this study, we surveyed which bacteria can directly and indirectly promote the growth of PM by using co-cultures with PM and molds associated with the fruiting body. Among 16 bacterial species isolated from the fruiting body, some species significantly influenced the mycelial growth of PM and molds. Most bacteria negatively affected PM growth and exhibited various enzyme activities, which suggests that they use the fruiting body as nutrient source. However, growth-promoting bacteria belonging to the Dietzia, Ewingella, Pseudomonas, Paenibacillus, and Rodococcus were also found. In addition, many bacteria suppressed molds, which suggests an indirect positive effect on PM as a biocontrol agent. Our results provide important insights toward a better understanding of the microbial interactions in the fruiting body of PM, and indicate that growth-promoting bacteria may be an important component in successful cultivation of PM.

Ectomycorrhization of Tricholoma matsutake with Quercus aquifolioides affects the endophytic microbial community of host plant

Tricholoma matsutake (S. Ito et Imai) is an ectomycorrhizal basidiomycete associated with Pinaceae and Fagaceae trees in the Northern Hemisphere. It is still unknown whether the symbiotic relationship with this ectomycorrhiza could affect the host plant's endophytic microbial community. In this study, we used high throughput sequencing to analyze the endophytic microbial communities of different Quercus aquifolioides tissues with or without T. matsutake partner. About 35,000 clean reads were obtained per sample, representing 34 bacterial phyla and 7 fungal phyla. We observed 3980 operational taxonomic units (OTUs) of bacteria and 457 OTUs of fungi at a 97% similarity level. Three bacterial phyla, Proteobacteria, Cyanobacteria, and Bacteroidetes, and the fungal phylum Ascomycota were dominant in all tissues. The relative abundance of these taxa differed significantly between Q. aquifolioides tissues with and without T. matsutake partner (p < 0.05). The bacterial genus Pseudomonas and the fungal genus Cryptosporiopsis were more abundant in mycorrhized roots than in control roots. This study showed that the community structure and dominant species of endophytic microbial communities in Q. aquifolioides tissues might be altered by colonization with T. matsutake. This work provides a new insight into the interactions between ectomycorrhizal fungus and host plant.

Land Cover and Climate Change May Limit Invasiveness of Rhododendron ponticum in Wales

Invasive plant species represent a serious threat to biodiversity precipitating a sustained global effort to eradicate or at least control the spread of this phenomenon. Current distribution ranges of many invasive species are likely to be modified in the future by land cover and climate change. Thus, invasion management can be made more effective by forecasting the potential spread of invasive species. Rhododendron ponticum (L.) is an aggressive invasive species which appears well suited to western areas of the UK. We made use of MAXENT modeling environment to develop a current distribution model and to assess the likely effects of land cover and climatic conditions (LCCs) on the future distribution of this species in the Snowdonia National park in Wales. Six global circulation models (GCMs) and two representative concentration pathways (RCPs), together with a land cover simulation for 2050 were used to investigate species' response to future environmental conditions. Having considered a range of environmental variables as predictors and carried out the AICc-based model selection, we find that under all LCCs considered in this study, the range of R. ponticum in Wales is likely to contract in the future. Land cover and topographic variables were found to be the most important predictors of the distribution of R. ponticum. This information, together with maps indicating future distribution trends will aid the development of mitigation practices to control R. ponticum.

Habitat suitability and movement corridors of grey wolf (Canis lupus) in Northern Pakistan

Habitat suitability models are useful to understand species distribution and to guide management and conservation strategies. The grey wolf (Canis lupus) has been extirpated from most of its historic range in Pakistan primarily due to its impact on livestock and livelihoods. We used non-invasive survey data from camera traps and genetic sampling to develop a habitat suitability model for C. lupus in northern Pakistan and to explore the extent of connectivity among populations. We detected suitable habitat of grey wolf using a maximum entropy approach (Maxent ver. 3.4.0) and identified suitable movement corridors using the Circuitscape 4.0 tool. Our model showed high levels of predictive performances, as seen from the values of area under curve (0.971±0.002) and true skill statistics (0.886±0.021). The main predictors for habitat suitability for C. lupus were distances to road, mean temperature of the wettest quarter and distance to river. The model predicted ca. 23,129 km2 of suitable areas for wolf in Pakistan, with much of suitable habitat in remote and inaccessible areas that appeared to be well connected through vulnerable movement corridors. These movement corridors suggest that potentially the wolf range can expand in Pakistan's Northern Areas. However, managing protected areas with stringent restrictions is challenging in northern Pakistan, in part due to heavy dependence of people on natural resources. The habitat suitability map provided by this study can inform future management strategies by helping authorities to identify key conservation areas.