Phylogeography and ecological niche modeling reveal evolutionary history of Leiolepis ocellata (Squamata,Leiolepidae)

Leiolepis ocellata is a lizard species distributing in topographically diverse habitats in northern Thailand. To explore its evolutionary history, 113 samples of L. ocellata were collected from 11 localities covering its distributional range in northern Thailand, and sequenced for mtDNA fragments (Cyt b and ND2). Pairwise comparisons across sampling localities yielded significant genetic differentiation (F _ST and Jost''s D) but no clear pattern of isolation by distance could be demonstrated based on the Mantel test. Phylogenetic and network analyses highlighted six haplogroups. Their divergence times were estimated to occur during the Pleistocene, much more recent than major orogenic events affecting northern Thailand. Instead, the results suggested that lineage divergences, of particularly eastern and western haplogroups of the region, coincided with the major rivers in the region (Yom river and Ping river, respectively), indicating vicariance in response to riverine barriers. Furthermore, ecological niche modeling suggested an expansion of suitable habitats of L. ocellata, when LGM‐liked conditions. This expansion potentially facilitated their dispersal among adjacent localities leading to lineage diversification and genetic admixture, after the riverine divergence.     

Integrating coalescent and ecological niche modeling in comparative phylogeography

Understanding the factors that contribute to the formation of population genetic structure is a central goal of phylogeographic research, but achieving this goal can be complicated by the stochastic variance inherent to genetic processes. Statistical approaches to testing phylogeographic hypotheses accommodate this stochasticity by evaluating competing models of putative historical population structure, often by simulating null distributions of the expected variance. The effectiveness of these tests depends on the biological realism of the models. Information from the fossil record can aid in reconstructing the historical distributions of some taxa. However, for the majority of taxa, which lack sufficient fossils, paleodistributional modeling can provide valuable spatial-geographic data concerning ancestral distributions. Paleodistributional models are generated by projecting ecological niche models, which predict the current distribution of each species, onto a model of past climatic conditions. Here, we generate paleodistributional models describing the suitable habitat during the last glacial maximum for lineages from the mesic forests of the Pacific Northwest of North America, and use these models to generate alternative phylogeographic hypotheses. Coalescent simulations are then used to test these hypotheses to improve our understanding of the historical events that promoted the formation of population genetic structure in this ecosystem. Results from Pacific Northwest mesic forest organisms demonstrate the utility of these combined approaches. Paleodistribution models and population genetic structure are congruent across three amphibian lineages, suggesting that they have responded in a concerted manner to environmental change. Two other species, a willow and a water vole, despite being currently codistributed and having similar population genetic structure, were predicted by the paleodistributional model to have had markedly different distributions during the last glacial maximum. This suggests that congruent phylogeographic patterns can arise from incongruent ancestral distributions. Paleodistributional models introduce a much-needed spatial-geographic perspective to statistical phylogeography. In conjunction with coalescent models of population genetic structure, they have the potential to improve our understanding of the factors that promote population divergence and ultimately produce regional patterns of biodiversity.     

Genetic structure and ecological niche segregation of Indian gray mongoose (Urva edwardsii) in Iran

Combining genetic data with ecological niche models is an effective approach for exploring climatic and nonclimatic environmental variables affecting spatial patterns of intraspecific genetic variation. Here, we adopted this combined approach to evaluate genetic structure and ecological niche of the Indian gray mongoose (Urva edwardsii) in Iran, as the most western part of the species range. Using mtDNA, we confirmed the presence of two highly differentiated clades. Then, we incorporated ensemble of small models (ESMs) using climatic and nonclimatic variables with genetic data to assess whether genetic differentiation among clades was coupled with their ecological niche. Climate niche divergence was also examined based on a principal component analysis on climatic factors only. The relative habitat suitability values predicted by the ESMs for both clades revealed their niche separation. Between‐clade climate only niche comparison revealed that climate space occupied by clades is similar to some extent, but the niches that they utilize differ between the distribution ranges of clades. We found that in the absence of evidence for recent genetic exchanges, distribution models suggest the species occurs in different niches and that there are apparent areas of disconnection across the species range. The estimated divergence time between the two Iranian clades (4.9 Mya) coincides with the uplifting of the Zagros Mountains during the Early Pliocene. The Zagros mountain‐building event seems to have prevented the distribution of U. edwardsii populations between the western and eastern parts of the mountains as a result of vicariance events. Our findings indicated that the two U. edwardsii genetic clades in Iran can be considered as two conservation units and can be utilized to develop habitat‐specific and climate change‐integrated management strategies.     

The importance (or lack thereof) of niche divergence to the maintenance of a northern species complex: the case of the long‐toed salamander (Ambystoma macrodactylum Baird)

The relative importance of ecological vs. nonecological factors for the origin and maintenance of species is an open question in evolutionary biology. Young lineages – such as the distinct genetic groups that make up the ranges of many northern species – represent an opportunity to study the importance of ecological divergence during the early stages of diversification. Yet, few studies have examined the extent of niche divergence between lineages in previously glaciated regions and the role of ecology in maintaining the contact zones between them. In this study, we used tests of niche overlap in combination with ecological niche models to explore the extent of niche divergence between lineages of the long‐toed salamander (Ambystoma macrodactylum Baird) species complex and to determine whether contact zones correspond to (divergent) niche limits. We found limited evidence for niche divergence between the different long‐toed salamander lineages, substantial overlap in the predicted distribution of suitable climatic space for all lineages and range limits that are independent of niche limits. These results raise questions as to the importance of ecological divergence to the development of this widespread species complex and highlight the potential for non‐ecological factors to play a more important role in the maintenance of northern taxa.     

Phylogenetic niche conservatism and the evolutionary basis of ecological speciation

Phylogenetic niche conservatism (PNC) typically refers to the tendency of closely related species to be more similar to each other in terms of niche than they are to more distant relatives. This has been implicated as a potential driving force in speciation and other species‐richness patterns, such as latitudinal gradients. However, PNC has not been very well defined in most previous studies. Is it a pattern or a process? What are the underlying endogenous (e.g. genetic) and exogenous (e.g. ecological) factors that cause niches to be conserved? What degree of similarity is necessary to qualify as PNC? Is it possible for the evolutionary processes causing niches to be conserved to also result in niche divergence in different habitats? Here, we revisit these questions, codifying a theoretical and operational definition of PNC as a mechanistic evolutionary process resulting from several factors. We frame this both from a macroevolutionary and population‐genetic perspective. We discuss how different axes of physical (e.g. geographic) and environmental (e.g. climatic) heterogeneity interact with the fundamental process of PNC to produce different outcomes of ecological speciation. We also review tests for PNC, and suggest ways that these could be improved or better utilized in future studies. Ultimately, PNC as a process has a well‐defined mechanistic basis in organisms, and future studies investigating ecological speciation would be well served to consider this, and frame hypothesis testing in terms of the processes and expected patterns described herein. The process of PNC may lead to patterns where niches are conserved (more similar than expected), constrained (divergent within a limited subset of available niches), or divergent (less similar than expected), based on degree of phylogenetic relatedness between species.     

Taxonomy and ecological niche modeling: Implications for the conservation of wood partridges (genus Dendrortyx)

Ecological niche models (ENMs) have a wide range of biological applications, particularly in conservation. To build these models, two sources of information are needed: occurrence records for the species of interest and environmental variables. However, taxonomic limits are often unclear, and the selection of occurrence data depends on the species concept being used. In this study we generated ENMs based on different taxonomic levels within the Dendrortyx group, which is comprised of three species and several subspecies; we analyzed the geographic and ecological distribution patterns and discuss the implications for the biogeography and conservation of this group. Our results suggest that the area with suitable climate depends on the taxonomic category used in the model, which in turn affects the interpretation of the importance of different biogeographic barriers and introduces variation into the potential differentiation of Dendrortyx. In terms of conservation, Dendrortyx macroura and Dendrortyx leucophrys are in a low risk category, that of “least concern,” although they may be amended to a higher category when their allopatric lineages are considered as the units for modeling. We suggest carrying out an a priori taxonomic analysis to facilitate the empirical identification of the units to be modeled in order to allow for a better ecological and biogeographic interpretation and more sound conservation policies.     

Phylogeography and ecological niche modelling implicate coastal refugia and trans‐alpine dispersal of a New Zealand fungus beetle

The Last Glacial Maximum (LGM) severely restricted forest ecosystems on New Zealand’s South Island, but the extent of LGM distribution for forest species is still poorly understood. We used mitochondrial DNA phylogeography (COI) and ecological niche modelling (ENM) to identify LGM refugia for the mycophagous beetle Agyrtodes labralis (Leiodidae), a forest edge species widely distributed in the South Island. Both the phylogenetic analyses and the ENM indicate that A. labralis refuged in Kaikoura, Nelson, and along much of the South Island’s west coast. Phylogeography of this species indicates that recolonization of the largely deforested east and southeast South Island occurred in a west–east direction, with populations moving through the Southern Alps, and that the northern refugia participated little in interglacial population expansion. This contradicts published studies of other New Zealand species, in which recolonization occurs in a north–south fashion from many of the same refugia.     

Spatial distribution of sand fly species (Psychodidae: Phlebtominae), ecological niche,and climatic regionalization in zoonotic foci of cutaneous leishmaniasis,southwest of Iran

Cutaneous leishmaniasis (CL) is a complex vector‐borne disease caused by Leishmania parasites that are transmitted by the bite of several species of infected female phlebotomine sand flies. Monthly factor analysis of climatic variables indicated fundamental variables. Principal component‐based regionalization was used for recognition of climatic zones using a clustering integrated method that identified five climatic zones based on factor analysis. To investigate spatial distribution of the sand fly species, the kriging method was used as an advanced geostatistical procedure in the ArcGIS modeling system that is beneficial to design measurement plans and to predict the transmission cycle in various regions of Khuzestan province, southwest of Iran. However, more than an 80% probability of P. papatasi was observed in rainy and temperate bio‐climatic zones with a high potential of CL transmission. Finding P. sergenti revealed the probability of transmission and distribution patterns of a non‐native vector of CL in related zones. These findings could be used as models indicating climatic zones and environmental variables connected to sand fly presence and vector distribution. Furthermore, this information is appropriate for future research efforts into the ecology of Phlebotomine sand flies and for the prevention of CL vector transmission as a public health priority.     

Ecotone formation through ecological niche construction: the role of biodiversity and species interactions

Rapid changes in species composition, also known as ecotones, can result from various causes including rapid changes in environmental conditions, or physiological thresholds. The possibility that ecotones arise from ecological niche construction by ecosystem engineers has received little attention. In this study, we investigate how the diversity of ecosystem engineers, and their interactions, can give rise to ecotones. We build a spatially explicit dynamical model that couples a multispecies community and its abiotic environment. We use numerical simulations and analytical techniques to determine the biotic and abiotic conditions under which ecotone emergence is expected to occur, and the role of biodiversity therein. We show that the diversity of ecosystem engineers can lead to indirect interactions through the modification of their shared environment. These interactions, which can be either competitive or mutualistic, can lead to the emergence of discrete communities in space, separated by sharp ecotones where a high species turnover is observed. Considering biodiversity is thus critical when studying the influence of species–environment interactions on the emergence of ecotones. This is especially true for the wide range of species that have small to moderate effects on their environment. Our work highlights new mechanisms by which biodiversity loss could cause significant changes in spatial community patterns in changing environments.     

High invasion potential of Hydrilla verticillata in the Americas predicted using ecological niche modeling combined with genetic data

Ecological niche modeling is an effective tool to characterize the spatial distribution of suitable areas for species, and it is especially useful for predicting the potential distribution of invasive species. The widespread submerged plant Hydrilla verticillata (hydrilla) has an obvious phylogeographical pattern: Four genetic lineages occupy distinct regions in native range, and only one lineage invades the Americas. Here, we aimed to evaluate climatic niche conservatism of hydrilla in North America at the intraspecific level and explore its invasion potential in the Americas by comparing climatic niches in a phylogenetic context. Niche shift was found in the invasion process of hydrilla in North America, which is probably mainly attributed to high levels of somatic mutation. Dramatic changes in range expansion in the Americas were predicted in the situation of all four genetic lineages invading the Americas or future climatic changes, especially in South America; this suggests that there is a high invasion potential of hydrilla in the Americas. Our findings provide useful information for the management of hydrilla in the Americas and give an example of exploring intraspecific climatic niche to better understand species invasion.     

Phylogeography and paleodistribution models of a widespread birch (Betula platyphylla Suk.) across East Asia: Multiple refugia,multidirectional expansion,and heterogeneous genetic pattern

Widespread tree species cover large geographical areas and play important roles in various vegetation types. Understanding how these species responded to historical climatic changes is important for understanding community assembly mechanisms with evolutionary and conservation implications. However, the location of refugial areas and postglacial history of widespread trees in East Asia remain poorly known. We combined microsatellite data (63 populations, 1756 individuals) and ecological niche modeling to examine the range‐wide population diversity, genetic structure, and historical demography of a pioneer tree species, Asian white birch (Betula platyphylla Suk.) across East Asia. We found a north‐to‐south trend of declining genetic diversity and five clusters, corresponding to geographical regions. Different clusters were inferred to have diverged through Pleistocene climatic oscillations and have different expansion routes, leading to genetic admixture in some populations. Ecological niche models indicated that the distribution of B. platyphylla during the last glacial maximum still had a large latitude span with slight shifts toward southeast, and northern populations had more variable distribution ranges than those in the south during later climatic oscillations. Our results reflect the relatively stable distribution through the last glacial–interglacial cycles and recent multidirectional expansion of B. platyphylla, providing new hypotheses for the response pattern of widespread tree species to climate change. The gradual genetic pattern from northeast to southwest and alternative distribution dynamics possibly resulted from environmental differences caused by latitude and topographic heterogeneity.     

Extending ecological niche models to the past 120 000 years corroborates the lack of strong phylogeographic structure in the Crested Drongo (Dicrurus forficatus forficatus) on Madagascar

We conduct a phylogeographic study of the Crested Drongo (Dicrurus forficatus forficatus), a broadly distributed bird species on Madagascar. We first determined the demographic and spatial pattern inferred from mitochondrial and nuclear data, and then compared these results with predictions from a present to 0.120‐Myr‐old reconstruction of the spatial dynamics of the range of D. f. forficatus on Madagascar, enabling putative areas of stability (lineage persistence) to be detected. Weak genetic structure along an east–west pattern and comparatively low genetic diversity were recovered, with strong evidence of population expansion found at all ten loci sampled. The palaeoclimatic distribution models over the past 0.120 Myr suggest the presence of extensive areas of suitable climate in the east and west for the species since its colonization of Madagascar, a result in strong concordance with the spatial and genetic signal derived from our multilocus data set. © 2013 The Linnean Society of London     

Phylogenomic,morphological, and niche differentiation analyses unveil species delimitation and evolutionary history of endangered maples in Acer series Campestria (Sapindaceae)

Accurate species delimitation is crucial for biodiversity conservation. The Acer series Campestria comprises four species, A. campestre L., A. miyabei Maxim., A. miaotaiense P. C. Tsoong, and A. yangjuechi Fang & P. L. Chiu. To clarify controversies over the taxonomic status of the latter three endangered species, we undertook phylogenomic, morphological, and niche differentiation analyses in series Campestria. Our coalescent species tree of 544 and 77 single-copy nuclear genes supported series Campestria as monophyletic, with A. yangjuechi having the closest relationship with A. miaotaiense. However, in the plastome-derived tree based on 64 protein coding sequences, the four species did not cluster together, and each of them grouped with some other sympatric Acer species. Given this nuclear-cytoplasmic conflict, we hypothesize that A. yangjuechi have been subject to nuclear gene introgression and plastid (pt) capture involving another sympatric maple, that is, A. amplum Rehder. Principal component analysis and machine learning based on morphological data could not separate A. yangjuechi and A. miaotaiense, but they both could be clearly distinguished from A. miyabei. Moreover, the niche overlap tests of the two more widespread species, A. miyabei and A. miaotaiense, showed they clearly occupy distinct niches. Overall, we conclude that A. miyabei and A. miaotaiense are distinct species, while A. yangjuechi (endemic to Mt. Tianmu/East China) should be treated as a subspecies of A. miaotaiense. Our study points out that multiple lines of phylogenomic, morphological, and ecological evidence prove highly useful in species delimitation. Additionally, our results should help to inform conservation measures for endangered species of the genus Acer/series Campestria in East Asia.     

Combined ecological niche modelling and molecular phylogeography revealed the evolutionary history of Hordeum marinum (Poaceae)--niche differentiation, loss of genetic diversity, and speciation in Mediterranean Quaternary refugia

The Hordeum marinum species group consists of two annual grasses of western Eurasian saline meadows or marshes. The two grasses split in the Quaternary about two million years ago. Hordeum marinum and the diploid of Hordeum gussoneanum (2x) co-occur throughout the Mediterranean basin, while the autotetraploid cytotype of H. gussoneanum (4x) overlaps with its diploid progenitor geographically only in the utmost Eastern Mediterranean, extending from there eastwards into Asia. Using chloroplast sequences of the trnL-F region, six newly developed chloroplast microsatellite loci, ecological predictive models based on climate data, and the present geographical distribution of the two species we analysed differentiation processes in the H. marinum group. The chloroplast data indicated clear differences in the history of both species. For H. marinum we found a subdivision between genetically variable populations from the Iberian Peninsula and the more uniform populations from the remaining Mediterranean. As an explanation, we assume Pleistocene fragmentation of an earlier widespread population and survival in an Iberian and a Central Mediterranean glacial refuge. Chloroplast variation was completely absent within the cytotypes of H. gussoneanum, indicating a severe and recent genetic bottleneck. Due to this lack of chloroplast variation only the combination of ecological habitat modelling with molecular data analyses allowed conclusions about the history of this taxon. The distribution areas of the two cytotypes of H. gussoneanum overlap today in parts of Turkey, indicating an area with similar climate conditions during polyploid formation. However, after its origin the polyploid cytotype underwent a pronounced ecological shift, compared to its diploid progenitor, allowing it to colonize mountainous inland habitats between the Mediterranean basin and Afghanistan. The extant sympatric occurrence of H. marinum and H. gussoneanum 2x in the Mediterranean region is interpreted as a result of secondary contact after fast Holocene range expansion out of different ice age refugia.     

Early stages of divergence: phylogeography, climate modeling, and morphological differentiation in the South American lizard Liolaemus petrophilus (Squamata: Liolaemidae)

This study examines the phylogeographic structure within the Patagonian lizard Liolaemus petrophilus and tests for patterns of between-clade morphological divergence and sexual dimorphism, as well as demographic and niche changes associated with Pleistocene climate changes. We inferred intraspecific relationships, tested hypotheses for historical patterns of population expansion, and incorporated ecological niche modeling (ENM) with standard morphological and geometric morphometric analyses to examine between-clade divergence as indirect evidence for adaptation to different niches. The two inferred haploclades diverged during the early Pleistocene with the Southern clade depicting the genetic signature of a recent population increase associated with expanding niche envelope, whereas the Northern clade shows stable populations in a shrinking niche envelope. The combination of molecular evidence for postisolation demographic change and ENM, suggest that the two haploclades have responded differently to Pleistocene climatic events.