The combination of niche modelling and landscape genetics (genomics) helps to disentangle processes that have shaped population structure in the evolutionary past and presence of species. Herein, we integrate a comprehensive genomic dataset with ecological parameters and niche modelling for the threatened Kaiser’s newt, a newt species adapted to mountain spring-ponds in Iran. Genomic analysis suggests the existence of two highly differentiated clades North and South of the Dez River. Genetic variation between the two clades (76.62%) was much greater than within clades (16.25%), suggesting that the Dez River prevented gene flow. River disconnectivity, followed by geographic distance, contributed mostly to genetic differentiation between populations. Environmental niche and landscape resistance had no significant influence. Though a significant difference between climatic niches occupied by each clade at the landscape-scale, habitat niches at the local-scale were equivalent. ‘Niche similarity analysis’ supported niche conservatism between the two clades despite the southward shift in the climatic niche of the Southern clade. Accordingly, populations of different clades may occupy different climatic niches within their ancestral niche. Our results indicate that the change of climatic conditions of geographically and genetically separated populations does not necessarily result in the shift of an ecological niche.
Background
The orogeny of the eastern Mediterranean region has substantially affected ecological speciation patterns, particularly of mountain-dwelling species. Mountain vipers of the genus Montivipera are among the paramount examples of Mediterranean neo-endemism, with restricted ranges in the mountains of Anatolia, the Levant, Caucasus, Alborz, and Zagros. Here we explore the phylogenetic and ecological diversification of Montivipera to reconstruct its ecological niche evolution and biogeographic history. Using 177 sequences of three mitochondrial genes, a dated molecular phylogeny of mountain vipers was reconstructed. Based on 320 occurrence points within the entire range of the genus and six climatic variables, ecological niches were modelled and used to infer ancestral niche occupancy. In addition, the biogeographic history and ancestral states of the species were reconstructed across climate gradients.
Results
Dated phylogenetic reconstruction revealed that the ancestor of mountain vipers split into two major clades at around 12.18 Mya followed by multiple vicariance events due to rapid orogeny. Montivipera colonised coastal regions from a mountain-dwelling ancestor. We detected a highly complex ecological niche evolution of mountain vipers to temperature seasonality, a variable that also showed a strong phylogenetic signal and high contribution in niche occupation.
Conclusion
Raising mountain belts in the Eastern Mediterranean region and subsequent remarkable changes in temperature seasonality have led to the formation of important centres of diversification and endemism in this biodiversity hotspot. High rates of niche conservatism, low genetic diversity, and segregation of ranges into the endemic distribution negatively influenced the adaptive capacity of mountain vipers. We suggest that these species should be considered as evolutionary significant units and priority species for conservation in Mediterranean mountain ecosystems.
The spatial scale of environmental layers is an important factor to consider in developing an understanding of ecological processes. This study employed Maxent modeling to investigate the geographic distribution of goitered gazelle, Gazella subgutturosa (Güldenstädt, 1780), in central Iran using uncorrelated variables at a spatial resolution of 250 m. We used spatial downscaling to downscale WorldClim data to 250-m resolution. We evaluated the sensitivity of the model to different grain and extent sizes from 250 m to 3 km. We compared the performance of the model at different scales using suitability indexes (AUC) and predicted habitat areas. Two models performed with AUC values higher than random (AUC un = 0.957, AUC pu = 0.953). The distribution of potential habitats at 250m grid size was strongly influenced by bioclimatic data, vegetation type and density, and elevation. There were few spatial divergences between uncorrelated and pruned models. The mean AUC across eight different spatial scales ranged from 0.936 to 0.959. There was a significant negative correlation between grain size and AUC (R 2 = 0.57). An increase in grain size increased the predicted habitat area. The extent size and AUC showed a positive correlation (R 2 = 0.18). Predicted suitability habitat also decreased as extent size increased (R 2 = 0.49). Spatial congruence AUC fluctuated within a small range and the maximum difference occurred between models of 1 × 1 and 2.5 × 2.5 km. These results showed that an increase in extent size is more accurate than an increase in grain size, and the maximum accuracy for predicting distribution of goitered gazelle in Iran was obtained if the grain size and extent size were 750 m.
Climate change is among the most important drivers of biodiversity decline through shift or shrinkage in suitable habitat of species. Mountain vipers of the genus
Montivipera
are under extreme risk from climate changes given their evolutionary history and geographic distribution. In this study, we divided all
Montivipera
species into three phylogenetic-geographic
Montivipera
clades (PGMC; Bornmuelleri, Raddei and Xanthina) and applied an ensemble ecological niche modelling (ENM) approach under different climatic scenarios to assess changes in projected suitable habitats of these species. Based on the predicted range losses, we assessed the projected extinction risk of the species relative to IUCN Red List Criteria. Our result revealed a strong decline in suitable habitats for all PGMCs (63.8%, 79.3% and 96.8% for Xanthina, Raddei and Bornmuelleri, respectively, by 2070 and under 8.5 RCP scenario) with patterns of altitudinal range shifts in response to projected climate change. We found that the mountains close to the Mediterranean Sea are exposed to the highest threats in the future (84.6 ± 9.1 percent range loss). We also revealed that disjunct populations of
Montivipera
will be additionally highly isolated and fragmented in the future. We argue that leveraging climate niche projections into the risk assessment provides the opportunity to implement IUCN criteria and better assess forthcoming extinction risks of species.
The Eurasian lynx (Lynx lynx) is one of the widespread felids in Eurasia; however, relatively little is known about the Asian subspecies, and especially the Iranian populations, which comprise the most southwestern part of its range. The current study aimed to assess the phylogenetic status of Iranian populations relative to other populations of Eurasia, by sequencing a 613 bp fragment of the mitochondrial control region. In total, 44 haplotypes were recorded from 83 sequences throughout Eurasia, two of which were found in Iran. The haplotype (H1) is dominant in all Iranian lynx populations and identical to specimens from SW Russia and central China. The second haplotype (H2) is unique and was recorded only from Ghazvin Province in the central Alborz Mountains. Both haplotypes occur in Ghazvin Province. The phylogenetic tree and a median-joining network identified four clades (i.e., East, West 1, West 2, and South). These results are congruent with previous studies and suggest that Eurasian lynx was restricted to the southern part of its range during the glacial maxima and expanded from there to East Asia and to Europe during several independent re-colonization events. The Caucasus region most like plays an important role as a refugium during glacial cycles.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.