Climatic stability and genetic divergence in the tropical insular lizard Anolis krugi, the Puerto Rican 'Lagartijo Jardinero de la Montaña'

History of Diversification and Adaptation from North to South Revealed by Genomic Data: Guanacos from the Desert to Sub-Antarctica

The increased availability of quality genomic data has greatly improved the scope and resolution of our understanding of the recent evolutionary history of wild species adapted to extreme environments and their susceptibility to anthropogenic impacts. The guanaco (Lama guanicoe), the largest wild ungulate in South America, is a good example. The guanaco is well adapted to a wide range of habitats, including the Sechura Desert, the high Andes Mountains to the north, and the extreme temperatures and conditions of Navarino Island to the south. Guanacos also have a long history of overexploitation by humans. To assess the evolutionary impact of these challenging habitats on the genomic diversity, we analyzed 38 genomes (∼10 to 16×) throughout their extensive latitudinal distribution from the Sechura and Atacama Desert to southward into Tierra del Fuego Island. These included analyses of patterns of unique differentiation in the north and geographic region further south with admixture among L. g. cacsilensis and L. g. guanicoe. Our findings provide new insights on the divergence of the subspecies ∼800,000 yr BP and document two divergent demographic trajectories and to the initial expansion of guanaco into the more southern portions of the Atacama Desert. Patagonian guanacos have experienced contemporary reductions in effective population sizes, likely the consequence of anthropogenic impacts. The lowest levels of genetic diversity corresponded to their northern and western limits of distribution and some varying degrees of genetic differentiation. Adaptive genomic diversity was strongly linked with environmental variables and was linked with colonization toward the south followed by adaptation.

Climate warming will increase chances of hybridization and introgression between two Takydromus lizards (Lacertidae)

Coexisting species may experience population and range changes alone or jointly in response to environmental change. Here, we used six climate variables and ten modeling algorithms to predict the distribution of two Takydromus species (T. septentrionalis and T. sexlineatus) in China. We identified the sympatric and allopatric areas by comparing projections between the two species based on habitat suitability under present and future climate scenarios. We constructed the hypervolumes of six climate variables for the two species and then evaluated overlaps between hypervolumes. From this study, we know the following. First, minimum temperature of coldest month contributes the most to the prediction of habitat suitability. Second, habitats suitable for the two species will shift northward in response to climate warming. Third, the range of T. sexlineatus will expand across the four future time intervals before 2,100, namely the 2021-2040, 2041-2060, 2061-2080, and 2081-2100 intervals, under both Shared socioeconomic pathway (SSP) 245 and SSP585 scenarios, and the range of T. septentrionalis will also expand in the future except at the 2081-2100 interval under the SSP585 scenario. Fourth, the sympatric areas will contract or expand under the SSP245 scenario and expand across the four future time intervals before 2,100 under the SSP585 scenario. Fifth, the niche hypervolumes of the two species partially overlapped, and the differences in niche centroid show some degree of niche differentiation between the two species. These results allow to conclude that climate warming will not only drive the northward drift of sympatric areas but also increase the size of these areas if nothing is done to limit the emission of greenhouse gases. Given the existence of hybridization and introgression between T. septentrionalis and T. sexlineatus in the field where they coexist, we also conclude that climate warming will increase chances of hybridization and introgression between the two species.

Climate as a driver of tropical insular diversity: comparative phylogeography of two ecologically distinctive frogs in Puerto Rico

The effects of late Quaternary climate on distributions and evolutionary dynamics of insular species are poorly understood in most tropical archipelagoes. We used ecological niche models under past and current climate to derive hypotheses regarding how stable climatic conditions shaped genetic diversity in two ecologically distinctive frogs in Puerto Rico. Whereas the Mountain Coquí, Eleutherodactylus portoricensis, is restricted to montane forest in the Cayey and Luquillo Mountains, the Red-eyed Coquí, E. antillensis, is a habitat generalist distributed across the entire Puerto Rican Bank (Puerto Rico and the Virgin Islands, excluding St. Croix). To test our hypotheses, we conducted phylogeographic and population genetic analyses based on mitochondrial and nuclear loci of each species across their range in Puerto Rico. Patterns of population differentiation in E. portoricensis, but not in E. antillensis, supported our hypotheses. For E. portoricensis, these patterns include: individuals isolated by long-term unsuitable climate in the Río Grande de Loíza Basin in eastern Puerto Rico belong to different genetic clusters; past and current climate strongly predicted genetic differentiation; and Cayey and Luquillo Mountains populations split prior to the last interglacial. For E. antillensis, these patterns include: genetic clusters did not fully correspond to predicted long-term unsuitable climate; and past and current climate weakly predicted patterns of genetic differentiation. Genetic signatures in E. antillensis are consistent with a recent range expansion into western Puerto Rico, possibly resulting from climate change and anthropogenic influences. As predicted, regions with a large area of long-term suitable climate were associated with higher genetic diversity in both species, suggesting larger and more stable populations. Finally, we discussed the implications of our findings for developing evidence-based management decisions for E. portoricensis, a taxon of special concern. Our findings illustrate the role of persistent suitable climatic conditions in promoting the persistence and diversification of tropical island organisms.

Biogeography and evolution of a widespread Central American lizard species complex: Norops humilis, (Squamata: Dactyloidae)

Background

Caribbean anole lizards (Dactyloidae) have frequently been used as models to study questions regarding biogeography and adaptive radiations, but the evolutionary history of Central American anoles (particularly those of the genus Norops) has not been well studied. Previous work has hypothesized a north-to-south dispersal pattern of Central American Norops, but no studies have examined dispersal within any Norops lineages. Here we test two major hypotheses for the dispersal of the N. humilis/quaggulus complex (defined herein, forming a subset within Savage and Guyer’s N. humilis group).

Results

Specimens of the N. humilis group were collected in Central America, from eastern Mexico to the Canal Zone of Panama. Major nodes were dated for comparison to the geologic history of Central America, and ancestral ranges were estimated for the N. humilis/quaggulus complex to test hypothesized dispersal patterns. These lineages displayed a northward dispersal pattern. We also demonstrate that the N. humilis/quaggulus complex consists of a series of highly differentiated mitochondrial lineages, with more conserved nuclear evolution. The paraphyly of the N. humilis species group is confirmed. A spatial analysis of molecular variance suggests that current populations are genetically distinct from one another, with limited mitochondrial gene flow occurring among sites.

Conclusions

The observed south-to-north colonization route within the Norops humilis/quaggulus complex represents the first evidence of a Norops lineage colonizing in a south-to-north pattern, (opposite to the previously held hypothesis for mainland Norops). One previously described taxon (N. quaggulus) was nested within N. humilis, demonstrating the paraphyly of this species; while our analyses also reject the monophyly of the Norops humilis species group (sensu Savage and Guyer), with N. tropidonotus, N. uniformis, and N. marsupialis being distantly related to/highly divergent from the N. humilis/quaggulus complex. Our work sheds light on mainland anole biogeography and past dispersal events, providing a pattern to test against other groups of mainland anoles.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0391-4) contains supplementary material, which is available to authorized users.

Relationship between wild greylag and European domestic geese based on mitochondrial DNA

The origins of the European domestic goose are uncertain. The available information comes from archaeological findings and historical literature, but genetic evidence has hitherto been scarce. The domestic goose in Europe is derived from the greylag goose (Anser anser), but it is not known where the initial domestication took place and which of the two subspecies of greylag goose was ancestral. We aimed to determine the amount and geographical distribution of genetic diversity in modern populations of greylag geese as well as in different breeds of the domestic goose to make inferences about goose domestication. We studied DNA sequence variation in the mitochondrial control region of greylag geese from multiple populations across Europe and western Asia as well as specimens of domestic geese representing 18 modern breeds and individuals not belonging to any recognised breed. Our results show notable differences in genetic diversity between different greylag goose populations and the presence of six mitochondrial haplogroups which show a degree of geographical partitioning. The genetic diversity of the domestic goose is low, with 84% of sampled individuals having one of two major closely related haplotypes, suggesting that modern European domestic geese may derive from a narrow genetic base. The site of domestication remains unresolved, but domestic geese in Turkey were unusually diverse, indicating the importance of further sampling in the vicinity of the eastern Mediterranean and the Near East. There appears to be past or ongoing hybridisation between greylags and domestic geese in particular areas, consistent with field observations.

Niche conservatism and the future potential range of Epipactis helleborine (Orchidaceae)

The aim of the present study was to evaluate the current distribution of suitable niches for the invasive orchid species, Epipactis helleborine, and to estimate the possibility of its further expansion. Moreover, niche modeling tools were used to explain its rapid expansion in North America and to test the niche conservatism of the species. The maximum entropy method was used to create models of the suitable niche distribution. A database of E. helleborine localities was prepared based on the examination of herbarium specimens, information from electronic databases as well as data gathered during field works. The differences between the niches occupied by native and invasive populations were evaluated using the niche overlap and niche identity test indexes. Moreover, the coverage of the most suitable habitats for the species was measured for three future scenarios as well as for the present time model. Populations of E. helleborine occupy North American west coast habitats very similar to those preferred by native, Eurasian populations, while the expansion in the east coast is related to the niche shift. The created models of suitable niche distribution indicate that the species does not realize its potential niche in the native range. The total surface of the habitats potentially available for E. helleborine will decrease in all climate change scenarios created for 2080.

Preliminary genetic analysis supports cave populations as targets for conservation in the endemic endangered Puerto Rican boa (Boidae: Epicrates inornatus)

The endemic Puerto Rican boa (Epicrates inornatus) has spent 42 years on the Endangered Species List with little evidence for recovery. One significant impediment to effective conservation planning has been a lack of knowledge of the distribution of genetic variability in the species. It has previously been suggested that boas might best be protected around caves that harbor large populations of bats. Prior study has found Puerto Rican boas at relatively high densities in and around bat caves, which they use both to feed and seek shelter. However, it is unknown whether these behaviorally distinctive populations represent a distinct evolutionary lineage, or (conversely) whether caves harbor representative genetic diversity for the species across the island. We provide the first genetic study of the Puerto Rican boa, and we examine and compare genetic diversity and divergence among two cave populations and two surface populations of boas. We find three haplogroups and an apparent lack of phylogeographic structure across the island. In addition, we find that the two cave populations appear no less diverse than the two surface populations, and harbor multiple mtDNA lineages. We discuss the conservation implications of these findings, including a call for the immediate protection of the remaining cave-associated populations of boas.

Genetic drift or natural selection? Hybridization and asymmetric mitochondrial introgression in two Caribbean lizards (Anolis pulchellus and Anolis krugi)

Hybridization and gene introgression can occur frequently between closely related taxa, but appear to be rare phenomena among members of the species-rich West Indian radiation of Anolis lizards. We investigated the pattern and possible mechanism of introgression between two sister species from Puerto Rico, Anolis pulchellus and Anolis krugi, using mitochondrial (ND2) and nuclear (DNAH3, NKTR) DNA sequences. Our findings demonstrated extensive introgression of A. krugi mtDNA (k-mtDNA) into the genome of A. pulchellus in western Puerto Rico, to the extent that k-mtDNA has mostly or completely replaced the native mtDNA of A. pulchellus on this part of the island. We proposed two not mutually exclusive scenarios to account for the interspecific matings between A. pulchellus and A. krugi. We inferred that hybridization events occurred independently in several populations, and determined that k-mtDNA haplotypes harboured in individuals of A. pulchellus can be assigned to four of the five major mtDNA clades of A. krugi. Further, the spatial distribution of k-mtDNA clades in the two species is largely congruent. Based on this evidence, we concluded that natural selection was the probable driving mechanism for the extensive k-mtDNA introgression into A. pulchellus. Our two nuclear data sets yielded different results. DNAH3 showed reciprocal monophyly of A. pulchellus and A. krugi, indicating no effect of hybridization on this marker. In contrast, the two species shared nine NKTR alleles, probably due to incomplete lineage sorting. Our study system will provide an excellent opportunity to experimentally assess the behavioural and ecological mechanisms that can lead to hybridization in closely related taxa.

Assessing the genetic landscape of a contact zone: the case of European hare in northeastern Greece

The European hare populations of the Balkan Peninsula comprise two divergent phylogenetic lineages with discrete geographical distribution slightly overlapping in the area of northeastern Greece and Bulgaria. Here we elucidate their contact zone, by defining the spatial distributional pattern of the two highly divergent groups, detecting individuals of hybrid origin, and identifying genetic barriers present in the area of their co-existence. Specimens from northeastern Greece were assayed for lineage assignment and population genetic inference based on a 511 bp fragment of mitochondrial DNA control region and allelic data from 10 microsatellite loci. Bayesian analyses on original and simulated genotypes were performed allowing for the contact zone delineation. Our results indicate high genetic diversity in both nuclear and mitochondrial DNA, strong population structure and non random spatial distribution of the differentiated gene pools. The information provided by the two types of molecular markers yielded consistent results. This study comprises a fine scale analysis of the contact zone between the two evolutionary lineages of European brown hares in northeastern Greece. Specific questions on the spatial patterns where addressed for the first time. Furthermore, hypotheses regarding the presence of hybrids were also tested. As a result, interpretive power to the diversity patterns observed today in the Balkans was added and previously overlooked aspects of the species biology were highlighted.

Sea level, topography and island diversity: phylogeography of the Puerto Rican Red-eyed Coquí, Eleutherodactylus antillensis

Quaternary climatic oscillations caused changes in sea level that altered the size, number and degree of isolation of islands, particularly in land-bridge archipelagoes. Elucidating the demographic effects of these oscillations increases our understanding of the role of climate change in shaping evolutionary processes in archipelagoes. The Puerto Rican Bank (PRB) (Puerto Rico and the Eastern Islands, which comprise Vieques, Culebra, the Virgin Islands and associated islets) in the eastern Caribbean Sea periodically coalesced during glaciations and fragmented during interglacial periods of the quaternary. To explore population-level consequences of sea level changes, we studied the phylogeography of the frog Eleutherodactylus antillensis across the archipelago. We tested hypotheses encompassing vicariance and dispersal narratives by sequencing mtDNA (c. 552 bp) of 285 individuals from 58 localities, and four nuDNA introns (totalling c. 1633 bp) from 173 of these individuals. We found low support for a hypothesis of divergence of the Eastern Islands populations prior to the start of the penultimate interglacial c. 250 kya, and higher support for a hypothesis of colonization of the Eastern Islands from sources in eastern Puerto Rico during the penultimate and last glacial period, when a land bridge united the PRB. The Río Grande de Loíza Basin in eastern Puerto Rico delineates a phylogeographic break. Haplotypes shared between the PRB and St. Croix (an island c. 105 km south-east of this archipelago) likely represent human-mediated introductions. Our findings illustrate how varying degrees of connectivity and isolation influence the evolution of tropical island organisms.

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.

Synergy between Allopatry and Ecology in Population Differentiation and Speciation

The general diversity pattern of the Caribbean anole radiation has been described in detail; however, the actual mechanisms at the origin of their diversification remain controversial. In particular, the role of ecological speciation, and the relative importance of divergence in allopatry and in parapatry, is debated. We describe the genetic structure of anole populations across lineage contact zones and ecotones to investigate the effect of allopatric divergence, natural selection, and the combination of both factors on population differentiation. Allopatric divergence had no significant impact on differentiation across the lineage boundary, while a clear bimodality in genetic and morphological characters was observed across an ecotone within a single lineage. Critically, the strongest differentiation was observed when allopatry and ecology act together, leading to a sharp reduction in gene flow between two lineages inhabiting different habitats. We suggest that, for Caribbean anoles to reach full speciation, a synergistic combination of several historical and ecological factors may be requisite.