Disentangling the effects of geographic peripherality and habitat suitability on neutral and adaptive genetic variation in Swiss stone pine

It is generally accepted that the spatial distribution of neutral genetic diversity within a species’ native range mostly depends on effective population size, demographic history, and geographic position. However, it is unclear how genetic diversity at adaptive loci correlates with geographic peripherality or with habitat suitability within the ecological niche. Using exome‐wide genomic data and distribution maps of the Alpine range, we first tested whether geographic peripherality correlates with four measures of population genetic diversity at > 17,000 SNP loci in 24 Alpine populations (480 individuals) of Swiss stone pine (Pinus cembra) from Switzerland. To distinguish between neutral and adaptive SNP sets, we used four approaches (two gene diversity estimates, F_ST outlier test, and environmental association analysis) that search for signatures of selection. Second, we established ecological niche models for P. cembra in the study range and investigated how habitat suitability correlates with genetic diversity at neutral and adaptive loci. All estimates of neutral genetic diversity decreased with geographic peripherality, but were uncorrelated with habitat suitability. However, heterozygosity (H_e) at adaptive loci based on Tajima's D declined significantly with increasingly suitable conditions. No other diversity estimates at adaptive loci were correlated with habitat suitability. Our findings suggest that populations at the edge of a species' geographic distribution harbour limited neutral genetic diversity due to demographic properties. Moreover, we argue that populations from suitable habitats went through strong selection processes, are thus well adapted to local conditions, and therefore exhibit reduced genetic diversity at adaptive loci compared to populations at niche margins.     

Idiosyncratic responses to climate‐driven forest fragmentation and marine incursions in reed frogs from Central Africa and the Gulf of Guinea Islands

Organismal traits interact with environmental variation to mediate how species respond to shared landscapes. Thus, differences in traits related to dispersal ability or physiological tolerance may result in phylogeographic discordance among co‐distributed taxa, even when they are responding to common barriers. We quantified climatic suitability and stability, and phylogeographic divergence within three reed frog species complexes across the Guineo‐Congolian forests and Gulf of Guinea archipelago of Central Africa to investigate how they responded to a shared climatic and geological history. Our species‐specific estimates of climatic suitability through time are consistent with temporal and spatial heterogeneity in diversification among the species complexes, indicating that differences in ecological breadth may partly explain these idiosyncratic patterns. Likewise, we demonstrated that fluctuating sea levels periodically exposed a land bridge connecting Bioko Island with the mainland Guineo‐Congolian forest and that habitats across the exposed land bridge likely enabled dispersal in some species, but not in others. We did not find evidence that rivers are biogeographic barriers across any of the species complexes. Despite marked differences in the geographic extent of stable climates and temporal estimates of divergence among the species complexes, we recovered a shared pattern of intermittent climatic suitability with recent population connectivity and demographic expansion across the Congo Basin. This pattern supports the hypothesis that genetic exchange across the Congo Basin during humid periods, followed by vicariance during arid periods, has shaped regional diversity. Finally, we identified many distinct lineages among our focal taxa, some of which may reflect incipient or unrecognized species.     

Ecological niche models display nonlinear relationships with abundance and demographic performance across the latitudinal distribution of Astragalus utahensis (Fabaceae)

The potential for ecological niche models (ENMs) to accurately predict species' abundance and demographic performance throughout their geographic distributions remains a topic of substantial debate in ecology and biogeography. Few studies simultaneously examine the relationship between ENM predictions of environmental suitability and both a species' abundance and its demographic performance, particularly across its entire geographic distribution. Yet, studies of this type are essential for understanding the extent to which ENMs are a viable tool for identifying areas that may promote high abundance or performance of a species or how species might respond to future climate conditions. In this study, we used an ensemble ecological niche model to predict climatic suitability for the perennial forb Astragalus utahensis across its geographic distribution. We then examined relationships between projected climatic suitability and field‐based measures of abundance, demographic performance, and forecasted stochastic population growth (λ_s). Predicted climatic suitability showed a J‐shaped relationship with A. utahensis abundance, where low‐abundance populations were associated with low‐to‐intermediate suitability scores and abundance increased sharply in areas of high predicted climatic suitability. A similar relationship existed between climatic suitability and λ_s from the center to the northern edge of the latitudinal distribution. Patterns such as these, where density or demographic performance only increases appreciably beyond some threshold of climatic suitability, support the contention that ENM‐predicted climatic suitability does not necessarily represent a reliable predictor of abundance or performance across large geographic regions.     

Climatic niche divergence and habitat suitability of eight alien invasive weeds in China under climate change

Testing climatic niche divergence and modeling habitat suitability under conditions of climate change are important for developing strategies to limit the introduction and expansion of alien invasive weeds (AIWs) and providing important ecological and evolutionary insights. We assessed climatic niches in both native and invasive ranges as well as habitat suitability under climate change for eight representative Chinese AIWs from the American continent. We used climatic variables associated with occurrence records and developed ecological niche models with Maxent. Interestingly, the climatic niches of all eight AIWs diverged significantly between the native and invasive ranges (the American continent and China). Furthermore, the AIWs showed larger climatic niche breadths in the invasive ranges than in the native ranges. Our results suggest that climatic niche shifts between native and invasive ranges occurred. Thus, the occurrence records of both native and invasive regions must be considered when modeling and predicting the spatial distributions of AIWs under current and future climate scenarios. Owing to high habitat suitability, AIWs were more likely to expand into regions of low latitude, and future climate change was predicted to result in a shift in the AIWs in Qinghai and Tibet (regions of higher altitude) as well as Heilongjiang, Jilin, Liaoning, Inner Mongolia, and Gansu (regions of higher latitude). Our results suggest that we need measures to prevent and control AIW expansion at the country‐wide level.     

Adding more ecology into species delimitation: ecological niche models and phylogeography help define cryptic species in the black salamander (Aneides flavipunctatus)

Being able to efficiently and accurately delimit species is one of the most basic and important aspects of systematics because species are the fundamental unit of analysis in biogeography, ecology, and conservation. We present a rationale and approach for combining ecological niche modeling, spatially explicit analyses of environmental data, and phylogenetics in species delimitation, and we use our methodology in an empirical example focusing on Aneides flavipunctatus, the black salamander (Caudata: Plethodontidae), in California. We assess the relationships between genetic, environmental, and geographic distance among populations. We use 11 climatic variables and point locality data from public databases to create ecological niche models. The suitability of potential contact zones between parapatric lineages is also assessed using the data from ecological niche modeling. Phylogenetic analyses of portions of the mitochondrial genome reveal morphologically cryptic mitochondrial lineages in this species. In addition, we find that patterns of genetic divergence are strongly associated with divergence in the ecological niche. Our work demonstrates the ease and utility of using spatial analyses of environmental data and phylogenetics in species delimitation, especially for groups displaying fine-scaled endemism and cryptic species.     

Are species' responses to global change predicted by past niche evolution?

Predicting how and when adaptive evolution might rescue species from global change, and integrating this process into tools of biodiversity forecasting, has now become an urgent task. Here, we explored whether recent population trends of species can be explained by their past rate of niche evolution, which can be inferred from increasingly available phylogenetic and niche data. We examined the assemblage of 409 European bird species for which estimates of demographic trends between 1970 and 2000 are available, along with a species-level phylogeny and data on climatic, habitat and trophic niches. We found that species'' proneness to demographic decline is associated with slow evolution of the habitat niche in the past, in addition to certain current-day life-history and ecological traits. A similar result was found at a higher taxonomic level, where families prone to decline have had a history of slower evolution of climatic and habitat niches. Our results support the view that niche conservatism can prevent some species from coping with environmental change. Thus, linking patterns of past niche evolution and contemporary species dynamics for large species samples may provide insights into how niche evolution may rescue certain lineages in the face of global change.     

Historical demographic dynamics underlying local adaptation in the presence of gene flow

The range of a species is the result of the relative contribution of spatial tracking of environmental requirements and adaptation to ecological conditions outside the ancestral niche. The appearance of novel habitats caused by climatic oscillation can promote range expansion and accompanying demographic growth. The demographic dynamics of populations leave a signal in \ patterns. We modeled three competing scenarios pertaining to the circumstance of a range expansion by the Karoo Scrub‐Robin into newly available habitat resulting from the increasing aridification of southern Africa. Genetic variation was contrasted with the theoretical expectations of a spatial range expansion, and compared with data of a putative adaptive trait. We infer that this bird likely colonized the arid zone, as a consequence of adaptive evolution in a small peripheral population, followed by an expansion with recurrent exchange of migrants with the ancestral populations.     

How to understand species’ niches and range dynamics: a demographic research agenda for biogeography

Range dynamics causes mismatches between a species’ geographical distribution and the set of suitable environments in which population growth is positive (the Hutchinsonian niche). This is because source–sink population dynamics cause species to occupy unsuitable environments, and because environmental change creates non‐equilibrium situations in which species may be absent from suitable environments (due to migration limitation) or present in unsuitable environments that were previously suitable (due to time‐delayed extinction). Because correlative species distribution models do not account for these processes, they are likely to produce biased niche estimates and biased forecasts of future range dynamics. Recently developed dynamic range models (DRMs) overcome this problem: they statistically estimate both range dynamics and the underlying environmental response of demographic rates from species distribution data. This process‐based statistical approach qualitatively advances biogeographical analyses. Yet, the application of DRMs to a broad range of species and study systems requires substantial research efforts in statistical modelling, empirical data collection and ecological theory. Here we review current and potential contributions of these fields to a demographic understanding of niches and range dynamics. Our review serves to formulate a demographic research agenda that entails: (1) advances in incorporating process‐based models of demographic responses and range dynamics into a statistical framework, (2) systematic collection of data on temporal changes in distribution and abundance and on the response of demographic rates to environmental variation, and (3) improved theoretical understanding of the scaling of demographic rates and the dynamics of spatially coupled populations. This demographic research agenda is challenging but necessary for improved comprehension and quantification of niches and range dynamics. It also forms the basis for understanding how niches and range dynamics are shaped by evolutionary dynamics and biotic interactions. Ultimately, the demographic research agenda should lead to deeper integration of biogeography with empirical and theoretical ecology.     

Correlation between genetic diversity and environmental suitability: taking uncertainty from ecological niche models into account

The hindcast of shifts in the geographical ranges of species as estimated by ecological niche modelling (ENM) has been coupled with phylogeographical patterns, allowing the inference of past processes that drove population differentiation and genetic variability. However, more recently, some studies have suggested that maps of environmental suitability estimated by ENM may be correlated to species' abundance, raising the possibility of using environmental suitability to infer processes related to population demographic dynamics and genetic variability. In both cases, one of the main problems is that there is a wide variation in ENM development methods and climatic models. In this study, we analyse the relationship between heterozygosity (He) and environmental suitability from multiple ENMs for 25 population estimates for Dipteryx alata, a widely distributed, endemic tree species of the Cerrado region of central Brazil. We propose a new approach for generating a statistical distribution of correlations under randomly generated ENM. The confidence intervals from these distributions indicate how model selection with different properties affects the ability to detect a correlation of interest (e.g. the correlation between He and suitability). Additionally, our approach allows us to explore which particular ensemble of ENMs produces the better result for finding an association between environmental suitability and He. Caution is necessary when choosing a method or a climatic data set for modelling geographical distributions, but the new approach proposed here provides a conservative way to evaluate the ability of ensembles to detect patterns of interest.     

Pleistocene climatic changes drive diversification across a tropical savanna

Spatial responses of species to past climate change depend on both intrinsic traits (climatic niche breadth, dispersal rates) and the scale of climatic fluctuations across the landscape. New capabilities in generating and analysing population genomic data, along with spatial modelling, have unleashed our capacity to infer how past climate changes have shaped populations, and by extension, complex communities. Combining these approaches, we uncover lineage diversity across four codistributed lizards from the Australian Monsoonal Tropics and explore how varying climatic tolerances interact with regional climate history to generate common vs. disparate responses to late Pleistocene change. We find more divergent spatial structuring and temporal demographic responses in the drier Kimberley region compared to the more mesic and consistently suitable Top End. We hypothesize that, in general, the effects of species’ traits on sensitivity to climate fluctuation will be more evident in climatically marginal regions. If true, this points to the need in climatically marginal areas to craft more species‐(or trait)‐specific strategies for persistence under future climate change.     

Ecological Niche Transferability Using Invasive Species as a Case Study

Species distribution modeling is widely applied to predict invasive species distributions and species range shifts under climate change. Accurate predictions depend upon meeting the assumption that ecological niches are conserved, i.e., spatially or temporally transferable. Here we present a multi-taxon comparative analysis of niche conservatism using biological invasion events well documented in natural history museum collections. Our goal is to assess spatial transferability of the climatic niche of a range of noxious terrestrial invasive species using two complementary approaches. First we compare species’ native versus invasive ranges in environmental space using two distinct methods, Principal Components Analysis and Mahalanobis distance. Second we compare species’ native versus invaded ranges in geographic space as estimated using the species distribution modeling technique Maxent and the comparative index Hellinger’s I. We find that species exhibit a range of responses, from almost complete transferability, in which the invaded niches completely overlap with the native niches, to a complete dissociation between native and invaded ranges. Intermediate responses included expansion of dimension attributable to either temperature or precipitation derived variables, as well as niche expansion in multiple dimensions. We conclude that the ecological niche in the native range is generally a poor predictor of invaded range and, by analogy, the ecological niche may be a poor predictor of range shifts under climate change. We suggest that assessing dimensions of niche transferability prior to standard species distribution modeling may improve the understanding of species’ dynamics in the invaded range.     

Niche differentiation in a postglacial colonizer,the bank vole Clethrionomys glareolus

Species‐level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species‐wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid‐Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.     

Scale-sensitivity in the measurement and interpretation of environmental niches

Species environmental niches are central to ecology, evolution, and global change research, but their characterization and interpretation depend on the spatial scale (specifically, the spatial grain) of their measurement. We find that the spatial grain of niche measurement is usually uninformed by ecological processes and varies by orders of magnitude. We illustrate the consequences of this variation for the volume, position, and shape of niche estimates, and discuss how it interacts with geographic range size, habitat specialization, and environmental heterogeneity. Spatial grain significantly affects the study of niche breadth, environmental suitability, niche evolution, niche tracking, and climate change effects. These and other fields will benefit from a more mechanism-informed choice of spatial grain and cross-grain evaluations that integrate different data sources.     

Forecasting species ranges by statistical estimation of ecological niches and spatial population dynamics

Aim The study and prediction of species–environment relationships is currently mainly based on species distribution models. These purely correlative models neglect spatial population dynamics and assume that species distributions are in equilibrium with their environment. This causes biased estimates of species niches and handicaps forecasts of range dynamics under environmental change. Here we aim to develop an approach that statistically estimates process‐based models of range dynamics from data on species distributions and permits a more comprehensive quantification of forecast uncertainties. Innovation We present an approach for the statistical estimation of process‐based dynamic range models (DRMs) that integrate Hutchinson's niche concept with spatial population dynamics. In a hierarchical Bayesian framework the environmental response of demographic rates, local population dynamics and dispersal are estimated conditional upon each other while accounting for various sources of uncertainty. The method thus: (1) jointly infers species niches and spatiotemporal population dynamics from occurrence and abundance data, and (2) provides fully probabilistic forecasts of future range dynamics under environmental change. In a simulation study, we investigate the performance of DRMs for a variety of scenarios that differ in both ecological dynamics and the data used for model estimation. Main conclusions Our results demonstrate the importance of considering dynamic aspects in the collection and analysis of biodiversity data. In combination with informative data, the presented framework has the potential to markedly improve the quantification of ecological niches, the process‐based understanding of range dynamics and the forecasting of species responses to environmental change. It thereby strengthens links between biogeography, population biology and theoretical and applied ecology.     

Snakes on the Balearic Islands: An Invasion Tale with Implications for Native Biodiversity Conservation

Biological invasions are a major conservation threat for biodiversity worldwide. Islands are particularly vulnerable to invasive species, especially Mediterranean islands which have suffered human pressure since ancient times. In the Balearic archipelago, reptiles represent an outstanding case with more alien than native species. Moreover, in the last decade a new wave of alien snakes landed in the main islands of the archipelago, some of which were originally snake-free. The identification of the origin and colonization pathways of alien species, as well as the prediction of their expansion, is crucial to develop effective conservation strategies. In this study, we used molecular markers to assess the allochthonous status and the putative origin of the four introduced snake species (Hemorrhois hippocrepis, Malpolon monspessulanus, Macroprotodon mauritanicus and Rhinechis scalaris) as well as ecological niche models to infer their patterns of invasion and expansion based on current and future habitat suitability. For most species, DNA sequence data suggested the Iberian Peninsula as the potential origin of the allochthonous populations, although the shallow phylogeographic structure of these species prevented the identification of a restricted source-area. For all of them, the ecological niche models showed a current low habitat suitability in the Balearic, which is however predicted to increase significantly in the next few decades under climate change scenarios. Evidence from direct observations and spatial distribution of the first-occurrence records of alien snakes (but also lizards and worm lizards) suggest the nursery trade, and in particular olive tree importation from Iberian Peninsula, as the main pathway of introduction of alien reptiles in the Balearic islands. This trend has been reported also for recent invasions in NE Spain, thus showing that olive trees transplantation may be an effective vector for bioinvasion across the Mediterranean. The combination of molecular and ecological tools used in this study reveals a promising approach for the understanding of the complex invasion process, hence guiding conservation management actions.     

Altitude‐mediated soil properties,not geography or climatic distance,explain the distribution of a tropical endemic herb

Understanding the ecological processes that govern species'' range margins is a fundamental question in ecology with practical implications in conservation biology. The center‐periphery hypothesis predicts that organisms have higher abundance at the center of their geographic range. However, most tests of this hypothesis often used raster data, assuming that climatic conditions are consistent across one square km. This assumption is not always justified, particularly for mountainous species for which climatic conditions can vary widely across a small spatial scale. Previous studies rarely evenly sample occurrence data across the species'' distribution. In this study, we sampled an endemic perennial herb, Thunbergia atacorensis (Acanthanceae), throughout its range in West Africa using 54 plots and collected data on (a)biotic variables, the species density, leaf mass per area, and basal diameter. We built a structural equation model to test the direct and indirect effects of distance from geographic and climatic niche centers, and altitude on Thunbergia density as mediated by abiotic and biotic factors, population demographic structure, and individual size. Contrary to the prediction of the center‐periphery hypothesis, we found no significant effect of distance from geographic or climatic niche centers on plant density. This indicates that even the climatic center does not necessarily have optimal ecological conditions. In contrast, plant density varied with altitudinal gradient, but this was mediated by the effect of soil nitrogen and potassium which had positive effect on plant size. Surprisingly, we found no direct or mediating effect of interspecific competition on plant density. Altogether, our results highlight the role of geography, climatic, and ecological mismatch in predicting species distribution. Our study highlights that where altitudinal gradient is strong local‐scale heterogeneity in abiotic factors can play important role in shaping species range limits.     

The effects of phenotypic plasticity and local adaptation on forecasts of species range shifts under climate change

Species are the unit of analysis in many global change and conservation biology studies; however, species are not uniform entities but are composed of different, sometimes locally adapted, populations differing in plasticity. We examined how intraspecific variation in thermal niches and phenotypic plasticity will affect species distributions in a warming climate. We first developed a conceptual model linking plasticity and niche breadth, providing five alternative intraspecific scenarios that are consistent with existing literature. Secondly, we used ecological niche‐modeling techniques to quantify the impact of each intraspecific scenario on the distribution of a virtual species across a geographically realistic setting. Finally, we performed an analogous modeling exercise using real data on the climatic niches of different tree provenances. We show that when population differentiation is accounted for and dispersal is restricted, forecasts of species range shifts under climate change are even more pessimistic than those using the conventional assumption of homogeneously high plasticity across a species' range. Suitable population‐level data are not available for most species so identifying general patterns of population differentiation could fill this gap. However, the literature review revealed contrasting patterns among species, urging greater levels of integration among empirical, modeling and theoretical research on intraspecific phenotypic variation.     

On population abundance and niche structure

Recent published evidence indicates a negative correlation between density of populations and the distance of their environments to a suitably defined ‘niche centroid’. This empirical observation lacks theoretical grounds. We provide a theoretical underpinning for the empirical relationship between population density and position in niche space, and use this framework to understand the circumstances under which the relationship will fail. We propose a metapopulation model for the area of distribution, as a system of ordinary differential equations coupled with a dispersal kernel. We present an analytical approximation to the solution of the system as well as R code to solve the full model numerically. We use this tool to analyze various scenarios and assumptions. General and realistic demographic assumptions imply a good correlation between position in niche space and population abundance. Factors that modify this correlation are: transitory states, a heterogeneous spatial structure of suitability, and Allee effects. We also explain why the raw output of the niche modeling algorithm MaxEnt is not a good predictor of environmental suitability. Our results elucidate the empirical results for spatial patterns of population size in niche terms, and provide a theoretical basis for a structured theory of the niche.     

Environmental influences on spatial and temporal patterns of body-size variation in California ground squirrels (Spermophilus beecheyi)

Aim  In order to understand how ground squirrels ( Spermophilus beecheyi ) may respond to future environmental change, we investigated five biotic and environmental factors potentially responsible for explaining body-size variation in this species across California. We examined the concordance of spatial patterns with temporal body-size change since the last glacial maximum (LGM).
Location  California, western North America.
Methods  We quantified body size of modern populations of ground squirrels ( n  = 81) and used a model-selection approach to determine the best variables (sex, vegetation, number of congeners, temperature and/or precipitation) explaining geographical variation in body size among modern populations. We also quantified body size of one fossil population in northern California ( n  = 39) and compared temporal body-size change in S. beecheyi at this location since the LGM with model predictions.
Results  Body size of modern populations conformed to Bergmann's rule, with larger individuals in northern (wetter and cooler) portions of California. However, the models suggest that precipitation, rather than temperature or other variables, may best explain variation in body size across modern spatial gradients. Our conclusion is supported by the temporal data, demonstrating that the body size of S. beecheyi has increased in northern California since the LGM, concordant with precipitation but not temperature change in the region.
Main conclusions  Precipitation, rather than temperature, vegetation or number of congeneric species, was the main factor explaining both spatial and temporal patterns of body-size variation in S. beecheyi . The integration of space and time provides a powerful mechanism for predicting how local populations may respond to current and future climatic changes.     

Comparative and predictive phylogeography in the South American diagonal of open formations: Unravelling the biological and environmental influences on multitaxon demography

Phylogeography investigates historical drivers of the geographical distribution of intraspecific lineages. Particular attention has been given to ecological, climatic and geological processes in the diversification of the Neotropical biota. Several species sampled across the South American diagonal of open formations (DOF), comprising the Caatinga, Cerrado and Chaco biomes, experienced range shifts coincident with Quaternary climatic changes. However, comparative studies across different spatial, temporal and biological scales on DOF species are still meagre. Here, we combine phylogeographical model selection and machine learning predictive frameworks to investigate the influence of Pleistocene climatic changes on several plant and animal species from the DOF. We assembled mitochondrial/chloroplastic DNA sequences in public repositories and inferred the demographic responses of 44 species, comprising 70 intraspecific lineages of plants, lizards, frogs, spiders and insects. We then built a random forest model using biotic and abiotic information to identify the best predictors of demographic responses in the Pleistocene. Finally, we assessed the temporal synchrony of species demographic responses with hierarchical approximate Bayesian computation. Biotic variables related to population connectivity, gene flow and habitat preferences largely predicted how species responded to Pleistocene climatic changes, and demographic changes were synchronous primarily during the Middle Pleistocene. Although 22 (~31%) lineages underwent demographic expansion, presumably associated with the spread of aridity during the glacial Pleistocene periods, our findings suggest that nine lineages (~13%) exhibited the opposite response due to taxon-specific attributes.