Insights from ecological niche modeling on the taxonomic distinction and niche differentiation between the black‐spotted and red‐spotted tokay geckoes (Gekko gecko)

The black‐spotted tokay and the red‐spotted tokay are morphologically distinct and have largely allopatric distributions. The black‐spotted tokay is characterized by a small body size and dark skin with sundry spots, while the red‐spotted tokay has a relatively large body size and red spots. Based on morphological, karyotypic, genetic, and distribution differences, recent studies suggested their species status; however, their classifications remain controversial, and additional data such as ecological niches are necessary to establish firm hypotheses regarding their taxonomic status. We reconstructed their ecological niches models using climatic and geographic data. We then performed niche similarity tests (niche identity and background tests) and point‐based analyses to explore whether ecological differentiation has occurred, and whether such differences are sufficient to explain the maintenance of their separate segments of environmental ranges. We found that both niche models of the black‐ and the red‐spotted tokay had a good fit and a robust performance, as indicated by the high area under the curve (AUC) values (“black” = 0.982, SD = ± 0.002, “red” = 0.966 ± 0.02). Significant ecological differentiation across the entire geographic range was found, indicating that the involvement of ecological differentiation is important for species differentiation. Divergence along the environmental axes is highly associated with climatic conditions, with isothermality being important for the “black” form, while temperature seasonality, precipitation of warmest quarter, and annual temperature range together being important for the “red” form. These factors are likely important factors in niche differentiation between the two forms, which result in morphological replacement. Overall, beside morphological and genetic differentiation information, our results contribute to additional insights into taxonomic distinction and niche differentiation between the black‐ and the red‐spotted tokay.     

Patterns of differential introgression in a salamander hybrid zone: inferences from genetic data and ecological niche modelling

Hybrid zones have yielded considerable insight into many evolutionary processes, including speciation and the maintenance of species boundaries. Presented here are analyses from a hybrid zone that occurs among three salamanders –Plethodon jordani, Plethodon metcalfi and Plethodon teyahalee– from the southern Appalachian Mountains. Using a novel statistical approach for analysis of non‐clinal, multispecies hybrid zones, we examined spatial patterns of variation at four markers: single‐nucleotide polymorphisms (SNPs) located in the mtDNA ND2 gene and the nuclear DNA ILF3 gene, and the morphological markers of red cheek pigmentation and white flecks. Concordance of the ILF3 marker and both morphological markers across four transects is observed. In three of the four transects, however, the pattern of mtDNA is discordant from all other markers, with a higher representation of P. metcalfi mtDNA in the northern and lower elevation localities than is expected given the ILF3 marker and morphology. To explore whether climate plays a role in the position of the hybrid zone, we created ecological niche models for P. jordani and P. metcalfi. Modelling results suggest that hybrid zone position is not determined by steep gradients in climatic suitability for either species. Instead, the hybrid zone lies in a climatically homogenous region that is broadly suitable for both P. jordani and P. metcalfi. We discuss various selective (natural selection associated with climate) and behavioural processes (sex‐biased dispersal, asymmetric reproductive isolation) that might explain the discordance in the extent to which mtDNA and nuclear DNA and colour‐pattern traits have moved across this hybrid zone.     

Evidence of ecological niche shift in Rhododendron ponticum (L.) in Britain: Hybridization as a possible cause of rapid niche expansion

Biological invasions threaten global biodiversity and natural resources. Anticipating future invasions is central to strategies for combating the spread of invasive species. Ecological niche models are thus increasingly used to predict potential distribution of invasive species. In this study, we compare ecological niches of Rhododendron ponticum in its native (Iberian Peninsula) and invasive (Britain) ranges. Here, we test the conservation of ecological niche between invasive and native populations of R. ponticum using principal component analysis, niche dynamics analysis, and MaxEnt‐based reciprocal niche modeling. We show that niche overlap between native and invasive populations is very low, leading us to the conclusion that the two niches are not equivalent and are dissimilar. We conclude that R. ponticum occupies novel environmental conditions in Britain. However, the evidence of niche shift presented in this study should be treated with caution because of nonanalogue climatic conditions between native and invasive ranges and a small population size in the native range. We then frame our results in the context of contradicting genetic evidence on possible hybridization of this invasive species in Britain. We argue that the existing contradictory studies on whether hybridization caused niche shift in R. ponticum are not sufficient to prove or disprove this hypothesis. However, we present a series of theoretical arguments which indicate that hybridization is a likely cause of the observed niche expansion of R. ponticum in Britain.     

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.     

Development of diagnostic microsatellite markers from whole‐genome sequences of Ammodramus sparrows for assessing admixture in a hybrid zone

Studies of hybridization and introgression and, in particular, the identification of admixed individuals in natural populations benefit from the use of diagnostic genetic markers that reliably differentiate pure species from each other and their hybrid forms. Such diagnostic markers are often infrequent in the genomes of closely related species, and genomewide data facilitate their discovery. We used whole‐genome data from Illumina HiSeqS2000 sequencing of two recently diverged (600,000 years) and hybridizing, avian, sister species, the Saltmarsh (Ammodramus caudacutus) and Nelson's (A. nelsoni) Sparrow, to develop a suite of diagnostic markers for high‐resolution identification of pure and admixed individuals. We compared the microsatellite repeat regions identified in the genomes of the two species and selected a subset of 37 loci that differed between the species in repeat number. We screened these loci on 12 pure individuals of each species and report on the 34 that successfully amplified. From these, we developed a panel of the 12 most diagnostic loci, which we evaluated on 96 individuals, including individuals from both allopatric populations and sympatric individuals from the hybrid zone. Using simulations, we evaluated the power of the marker panel for accurate assignments of individuals to their appropriate pure species and hybrid genotypic classes (F1, F2, and backcrosses). The markers proved highly informative for species discrimination and had high accuracy for classifying admixed individuals into their genotypic classes. These markers will aid future investigations of introgressive hybridization in this system and aid conservation efforts aimed at monitoring and preserving pure species. Our approach is transferable to other study systems consisting of closely related and incipient species.     

From refugia to contact: Pine processionary moth hybrid zone in a complex biogeographic setting

Contact zones occur at the crossroad between specific dispersal routes and are facilitated by biogeographic discontinuities. Here, we focused on two Lepidoptera sister species that come in contact near the Turkish Straits System (TSS). We aimed to infer their phylogeographic histories in the Eastern Mediterranean and finely analyze their co‐occurrence and hybridization patterns in this biogeographic context. We used molecular mitochondrial and nuclear markers to study 224 individuals from 42 localities. We used discordances between markers and complementary assignment methods to identify and map hybrids and parental individuals. We confirmed the parapatric distribution of Thaumetopoea pityocampa (Lepidoptera: Notodontidae) in the west and Thaumetopoea wilkinsoni in the east and identified a narrow contact zone. We identified several glacial refugia of T. wilkinsoni in southern Turkey with a strong east–west differentiation in this species. Unexpectedly, T. pityocampa crossed the TSS and occur in northern Aegean Turkey and some eastern Greek islands. We found robust evidence of introgression between the two species in a restricted zone in northwestern Turkey, but we did not identify any F₁ individuals. The identified hybrid zone was mostly bimodal. The distributions and genetic patterns of the studied species were strongly influenced both by the Quaternary climatic oscillations and the complex geological history of the Aegean region. T. pityocampa and T. wilkinsoni survived the last glacial maximum in disjoint refugia and met in western Turkey at the edge of the recolonization routes. Expanding population of T. wilkinsoni constrained T. pityocampa to the western Turkish shore. Additionally, we found evidence of recurrent introgression by T. wilkinsoni males in several T. pityocampa populations. Our results suggest that some prezygotic isolation mechanisms, such as differences in timing of the adult emergences, might be a driver of the isolation between the sister species.     

Evidence for concerted movement of nuclear and mitochondrial clines in a lizard hybrid zone

Moving hybrid zones provide compelling examples of evolution in action, yet long‐term studies that test the assumptions of hybrid zone stability are rare. Using replicated transect samples collected over a 10‐year interval from 2002 to 2012, we find evidence for concerted movement of genetic clines in a plateau fence lizard hybrid zone (Sceloporus tristichus) in Arizona. Cline‐fitting analyses of SNP and mtDNA data both provide evidence that the hybrid zone shifted northward by approximately 2 km during the 10‐year interval. For each sampling period, the mtDNA cline centre is displaced from the SNP cline centre and maintaining an introgression distance of approximately 3 km. The northward expansion of juniper trees into the Little Colorado River Basin in the early 1900s provides a plausible mechanism for hybrid zone formation and movement, and a broadscale quantification of recent land cover change provides support for increased woody species encroachment at the southern end of the hybrid zone. However, population processes can also contribute to hybrid zone movement, and the current stability of the ecotone habitats in the centre of the hybrid zone suggests that movement could decelerate in the future.     

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.     

Predictive modeling for allopatric Strix (Strigiformes: Strigidae) owls in South America: determinants of their distributions and ecological niche‐based processes

Strix (Strigidae) is a worldwide genus of 17 owl species typical of forested habitats, including Rusty‐barred Owls (S. hylophila), Chaco Owls (S. chacoensis), and Rufous‐legged Owls (S. rufipes) in South America. These species are distributed allopatrically, but the ecological traits that determine their distributions remain largely unknown and their phylogenetic relationships are unclear. We used species distribution models (SDMs) to identify variables explaining their distribution patterns and test hypotheses about ecological divergence and conservatism based on niche overlap analysis. For Rusty‐barred Owls and Chaco Owls, climatic factors related to temperature played a major role, whereas a rainfall variable was more important for Rufous‐legged Owls. When niche overlaps were compared, accounting for regional similarities in the habitat available to each species, an ecological niche divergence process was supported for Chaco Owl‐Rusty‐barred Owl and Chaco Owl‐Rufous‐legged Owl, whereas a niche conservatism process was supported for Rusty‐barred Owl‐Rufous‐legged Owl. Different ecological requirements support current species delimitation, but they are in disagreement with the two main hypotheses currently envisaged about their phylogenetic relationships (Chaco Owls as the sister taxa of either Rufous‐legged Owls or Rusty‐barred Owls) and support a new phylogenetic hypothesis (Rufous‐legged Owls as sister taxa of Rusty‐barred Owls). Our findings suggest that speciation of Rusty‐barred Owls and Rufous‐legged Owls was a vicariant event resulting from Atlantic marine transgressions in southern South America in the Miocene, but their niche was conserved because habitat changed little in their respective ranges. In contrast, Chaco Owls diverged ecologically from the other two species as a result of their adaptations to the habitat they currently occupy. Ecological and historical approaches in biogeography can be embedded to explain distribution patterns, and results provided by SDMs can be used to infer historical and ecological processes in an integrative way.     

Clinal variation in a brown lemur (Eulemur spp.) hybrid zone: Combining morphological,genetic and climatic data to examine stability

Studies of hybrid zones can inform our understanding of reproductive isolation and speciation. Two species of brown lemur (Eulemur rufifrons and E. cinereiceps) form an apparently stable hybrid zone in the Andringitra region of south‐eastern Madagascar. The aim of this study was to identify factors that contribute to this stability. We sampled animals at 11 sites along a 90‐km transect through the hybrid zone and examined variation in 26 microsatellites, the D‐loop region of mitochondrial DNA, six pelage and nine morphological traits; we also included samples collected in more distant allopatric sites. Clines in these traits were noncoincident, and there was no increase in either inbreeding coefficients or linkage disequilibrium at the centre of the zone. These results could suggest that the hybrid zone is maintained by weak selection against hybrids, conforming to either the tension zone or geographical selection‐gradient model. However, a closer examination of clines in pelage and microsatellites indicates that these clines are not sigmoid or stepped in shape but instead plateau at their centre. Sites within the hybrid zone also occur in a distinct habitat, characterized by greater seasonality in precipitation and lower seasonality in temperature. Together, these findings suggest that the hybrid zone may follow the bounded superiority model, with exogenous selection favouring hybrids within the transitional zone. These findings are noteworthy, as examples supporting the bounded superiority model are rare and may indicate a process of ecologically driven speciation without geographical isolation.     

Spatial genetic features of eastern oysters (Crassostrea virginica Gmelin) in the Gulf of Mexico: northward movement of a secondary contact zone

The eastern oyster (Crassostrea virginica Gmelin) is an economically and ecologically valuable marine bivalve occurring in the Gulf of Mexico. This study builds upon previous research that identified two divergent populations of eastern oysters in the western Gulf of Mexico. Allelic and genotypic patterns from 11 microsatellite markers were used to assess genetic structure and migration between the previously described oyster populations in Texas. The main findings are as follows: (1) there are two distinct populations (F_ST = 0.392, P < 0.001) of oysters that overlap in the Corpus Christi/Aransas Bay estuarine complex in Texas, (2) the distribution of genotypes among individuals in the contact zone suggests limited hybridization between populations, (3) the variables of salinity, temperature, dissolved oxygen, turbidity and depth are not correlated with allele frequencies on reefs in the contact zone or when analyzed across Texas, and (4) there is little evidence of directional selection acting on the loci assayed here, although patterns at four markers suggested the influence of balancing selection based on outlier analyses. These results are consistent with long‐term historical isolation between populations, followed by secondary contact. Recent hydrological changes in the area of secondary contact may be promoting migration in areas that were previously inhospitable to eastern oysters, and observed differences in the timing of spawning may limit hybridization between populations. Comparison of these findings with the results of an earlier study of oysters in Texas suggests that the secondary contact zone has shifted approximately 27 km north, in as little as a 23‐year span.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Ecological niche modeling for conservation planning of an endemic snail in the verge of becoming a pest in cardamom plantations in the Western Ghats biodiversity hotspot

Conservation managers and policy makers are often confronted with a challenging dilemma of devising suitable strategies to maintain agricultural productivity while conserving endemic species that at the early stages of becoming pests of agricultural crops. Identification of environmental factors conducive to species range expansion for forecasting species distribution patterns will play a central role in devising management strategies to minimize the conflict between the agricultural productivity and biodiversity conservation. Here, we present results of a study that predicts the distribution of Indrella ampulla, a snail endemic to the Western Ghats biodiversity hotspot, which is becoming a pest in cardamom (Ellettaria cardamomum) plantations. We determined the distribution patterns and niche overlap between I. ampulla and Ellettaria cardamomum using maximum entropy (MaxEnt) niche modeling techniques under current and future (2020–2080) climatic scenarios. The results showed that climatic (precipitation of coldest quarter and isothermality) and soil (cation exchange capacity of soil [CEC]) parameters are major factors that determine the distribution of I. ampulla in Western Ghats. The model predicted cardamom cultivation areas in southern Western Ghats are highly sensitive to invasion of I. ampulla under both present and future climatic conditions. While the land area in the central Western Ghats is predicted to become unsuitable for I. ampulla and Ellettaria cardamomum in future, we found 71% of the Western Ghats land area is suitable for Ellettaria cardamomum cultivation and 45% suitable for I. ampulla, with an overlap of 35% between two species. The resulting distribution maps are invaluable for policy makers and conservation managers to design and implement management strategies minimizing the conflicts to sustain agricultural productivity while maintaining biodiversity in the region.     

Integrating Bayesian genomic cline analyses and association mapping of morphological and ecological traits to dissect reproductive isolation and introgression in a Louisiana Iris hybrid zone

Hybrid zones provide unique opportunities to examine reproductive isolation and introgression in nature. We utilized 45,384 single nucleotide polymorphism (SNP) loci to perform association mapping of 14 floral, vegetative and ecological traits that differ between Iris hexagona and Iris fulva, and to investigate, using a Bayesian genomic cline (BGC) framework, patterns of genomic introgression in a large and phenotypically diverse hybrid zone in southern Louisiana. Many loci of small effect size were consistently found to be associated with phenotypic variation across all traits, and several individual loci were revealed to influence phenotypic variation across multiple traits. Patterns of genomic introgression were quite heterogeneous throughout the Louisiana Iris genome, with I. hexagona alleles tending to be favoured over those of I. fulva. Loci that were found to have exceptional patterns of introgression were also found to be significantly associated with phenotypic variation in a small number of morphological traits. However, this was the exception rather than the rule, as most loci that were associated with morphological trait variation were not significantly associated with excess ancestry. These findings provide insights into the complexity of the genomic architecture of phenotypic differences and are a first step towards identifying loci that are associated with both trait variation and reproductive isolation in nature.     

Molecular data and ecological niche modelling reveal the Pleistocene history of a semi‐aquatic bug (Microvelia douglasi douglasi) in East Asia

This study investigated the Pleistocene history of a semi‐aquatic bug, Microvelia douglasi douglasi Scott, 1874 (Hemiptera: Veliidae) in East Asia. We used M. douglasi douglasi as a model species to explore the effects of historical climatic fluctuations on montane semi‐aquatic invertebrate species. Two hypotheses were developed using ecological niche models (ENMs). First, we hypothesized that M. douglasi douglasi persisted in suitable habitats in southern Guizhou, southern Yunnan, Hainan, Taiwan and southeast China during the LIG. After that, the populations expanded (Hypothesis 1). As the spatial prediction in the LGM was significantly larger than in the LIG, we then hypothesized that the population expanded during the LIG to LGM transition (Hypothesis 2). We tested these hypotheses using mitochondrial data (COI+COII) and nuclear data (ITS1 + 5.8S+ITS2). Young lineages, relatively deep splits, lineage differentiation among mountain ranges in central, south and southwest China and high genetic diversities were observed in these suitable habitats. Evidence of mismatch distributions and neutrality tests indicate that a population expansion occurred in the late Pleistocene. The Bayesian skyline plot (BSP) revealed an unusual population expansion that likely happened during the cooling transition between LIG and LGM. The results of genetic data were mostly consistent with the spatial predictions from ENM, a finding that can profoundly improve phylogeographic research. The ecological requirements of M. douglasi douglasi, together with the geographical heterogeneity and climatic fluctuations of Pleistocene in East Asia, could have shaped this unusual demographic history. Our study contributes to our knowledge of semi‐aquatic bug/invertebrate responses to Pleistocene climatic fluctuations in East Asia.     

Trophic niche ontogeny and palaeoecology of early Toarcian Stenopterygius (Reptilia: Ichthyosauria)

Reconstructing ecological niche shifts during ontogeny in extinct animals with no living analogues is difficult without exceptional fossil collections. Here we demonstrate how a previously identified ontogenetic shift in the size and shape of the dentition in the early Toarcian ichthyosaur Stenopterygius quadriscissus accurately predicts a particular dietary shift. The smallest S. quadriscissus fed on small, burst‐swimming fishes, with a steady shift towards faster moving fish and cephalopods with increasing body size. Larger adult specimens appear to have been completely reliant on cephalopods, with fish completely absent from gut contents shortly after onset of sexual maturity. This is consistent with a previously proposed ontogenetic niche shift based on tooth shape and body size, corroborating the idea that dental ontogeny may be a useful predictor of dietary shifts in marine reptiles. Applying the theoretical framework used here to other extinct species will improve the resolution of palaeoecological reconstructions, where appropriate sample sizes exist.     

Ecological niche comparison and molecular phylogeny segregate the invasive moss species Campylopus introflexus (Leucobryaceae,Bryophyta) from its closest relatives

The delimitation of the invasive moss species Campylopus introflexus from its closest relative, Campylopus pilifer, has been long debated based on morphology. Previous molecular phylogenetic reconstructions based on the nuclear ribosomal internal transcribed spacers (ITS) 1 and 2 showed that C. pilifer is split into an Old World and a New World lineage, but remained partly inconclusive concerning the relationships between these two clades and C. introflexus. Analyses of an extended ITS dataset displayed statistically supported incongruence between ITS1 and ITS2. ITS1 separates the New World clade of C. pilifer from a clade comprising C. introflexus and the Old World C. pilifer. Ancestral state reconstruction showed that this topology is morphologically supported by differences in the height of the dorsal costal lamellae in leaf cross‐section (despite some overlap). ITS2, in contrast, supports the current morphological species concept, i.e., separating C. introflexus from C. pilifer, which is morphologically supported by the orientation of the hyaline hair point at leaf apex as well as costal lamellae height. Re‐analysis of published and newly generated plastid atpB‐rbcL spacer sequences supported the three ITS lineages. Ecological niche modeling proved a useful approach and showed that all three molecular lineages occupy distinct environmental spaces that are similar, but undoubtedly not equivalent. In line with the ITS1 topology, the C. pilifer lineage from the New World occupies the most distinct environmental niche, whereas the niches of Old World C. pilifer and C. introflexus are very similar. Taking the inferences from ecological niche comparisons, phylogenetics, and morphology together, we conclude that all three molecular lineages represent different taxa that should be recognized as independent species, viz. C. introflexus, C. pilifer (Old World clade), and the reinstated C. lamellatus Mont. (New World clade).