A mid‐Pleistocene rainforest corridor enabled synchronous invasions of the Atlantic Forest by Amazonian anole lizards

Shifts in the geographic distribution of habitats over time can promote dispersal and vicariance, thereby influencing large‐scale biogeographic patterns and ecological processes. An example is that of transient corridors of suitable habitat across disjunct but ecologically similar regions, which have been associated with climate change over time. Such connections likely played a role in the assembly of tropical communities, especially within the highly diverse Amazonian and Atlantic rainforests of South America. Although these forests are presently separated by open and dry ecosystems, paleoclimatic and phylogenetic evidence suggest that they have been transiently connected in the past. However, little is known about the timing, magnitude and the distribution of former forest connections. We employ sequence data at multiple loci from three codistributed arboreal lizards (Anolis punctatus, Anolis ortonii and Polychrus marmoratus) to infer the phylogenetic relationships among Amazonian and Atlantic Forest populations and to test alternative historical demographic scenarios of colonization and vicariance using coalescent simulations and approximate Bayesian computation (ABC). Data from the better‐sampled Anolis species support colonization of the Atlantic Forest from eastern Amazonia. Hierarchical ABC indicates that the three species colonized the Atlantic Forest synchronously during the mid‐Pleistocene. We find support of population bottlenecks associated with founder events in the two Anolis, but not in P. marmoratus, consistently with their distinct ecological tolerances. Our findings support that climatic fluctuations provided key opportunities for dispersal and forest colonization in eastern South America through the cessation of environmental barriers. Evidence of species‐specific histories strengthens assertions that biological attributes play a role in responses to shared environmental change.     

Postglacial climate changes and rise of three ecotypes of harbour porpoises,Phocoena phocoena,in western Palearctic waters

Despite no obvious barriers to gene flow in the marine realm, environmental variation and ecological specializations can lead to genetic differentiation in highly mobile predators. Here, we investigated the genetic structure of the harbour porpoise over the entire species distribution range in western Palearctic waters. Combined analyses of 10 microsatellite loci and a 5085 base‐pair portion of the mitochondrial genome revealed the existence of three ecotypes, equally divergent at the mitochondrial genome, distributed in the Black Sea (BS), the European continental shelf waters, and a previously overlooked ecotype in the upwelling zones of Iberia and Mauritania. Historical demographic inferences using approximate Bayesian computation (ABC) suggest that these ecotypes diverged during the last glacial maximum (c. 23–19 kilo‐years ago, kyrbp ). ABC supports the hypothesis that the BS and upwelling ecotypes share a more recent common ancestor (c. 14 kyrbp ) than either does with the European continental shelf ecotype (c. 28 kyrbp ), suggesting they probably descended from the extinct populations that once inhabited the Mediterranean during the glacial and post‐glacial period. We showed that the two Atlantic ecotypes established a narrow admixture zone in the Bay of Biscay during the last millennium, with highly asymmetric gene flow. This study highlights the impacts that climate change may have on the distribution and speciation process in pelagic predators and shows that allopatric divergence can occur in these highly mobile species and be a source of genetic diversity.     

Mitochondrial introgressive hybridization following a demographic expansion in the tomato frogs of Madagascar,genus Dyscophus

Madagascar is a biodiversity hotspot with a unique fauna and flora largely endemic at the species level and highly threatened by habitat destruction. The processes underlying population‐level differentiation in Madagascar's biota are poorly understood and have been proposed to be related to Pleistocene climatic cycles, yet the levels of genetic divergence observed are often suggestive of ancient events. We combined molecular markers of different variability to assess the phylogeography of Madagascar's emblematic tomato frogs (Dyscophus guineti and D. antongilii) and interpret the observed pattern as resulting from ancient and recent processes. Our results suggest that the initial divergence between these taxa is probably old as reflected by protein‐coding nuclear genes and by a strong mitochondrial differentiation of the southernmost population. Dramatic changes in their demography appear to have been triggered by the end of the last glacial period and possibly by the short return of glacial conditions known as the 8K event. This dramatic change resulted in an approximately 50‐fold reduction of the effective population size in various populations of both species. We hypothesize these species' current mitochondrial DNA diversity distribution reflects a swamping of the mitochondrial genetic diversity of D. guineti by that of D. antongilii previous to the populations' bottlenecks during the Holocene, and probably as a consequence of D. antongilii demographic expansion approximately 1 million years ago. Our data support the continued recognition of D. antongilii and D. guineti as separate species and flag D. guineti as the more vulnerable species to past and probably also future environmental changes.     

Postglacial range shift and demographic expansion of the marine intertidal snail Batillaria attramentaria

To address the impacts of past climate changes, particularly since the last glacial period, on the history of the distribution and demography of marine species, we investigated the evolutionary and demographic responses of the intertidal batillariid gastropod, Batillaria attramentaria, to these changes, using the snail as a model species in the northwest Pacific. We applied phylogeographic and divergence population genetic approaches to mitochondrial COI sequences from B. attramentaria. To cover much of its distributional range, 197 individuals collected throughout Korea and 507 publically available sequences (mostly from Japan) were used. Finally, a Bayesian skyline plot (BSP) method was applied to reconstruct the demographic history of this species. We found four differentiated geographic groups around Korea, confirming the presence of two distinct, geographically subdivided haplogroups on the Japanese coastlines along the bifurcated routes of the warm Tsushima and Kuroshio Currents. These two haplogroups were estimated to have begun to split approximately 400,000 years ago. Population divergence analysis supported the hypothesis that the Yellow Sea was populated by a northward range expansion of a small fraction of founders that split from a southern ancestral population since the last glacial maximum (LGM: 26,000–19,000 years ago), when the southern area became re‐submerged. BSP analyses on six geographically and genetically defined groups in Korea and Japan consistently demonstrated that each group has exponentially increased approximately since the LGM. This study resolved the phylogeography of B. attramentaria as a series of events connected over space and time; while paleoceanographic conditions determining the connectivity of neighboring seas in East Asia are responsible for the vicariance of this species, the postglacial sea‐level rise and warming temperatures have played a crucial role in rapid range shifts and broad demographic expansions of its populations.     

Stability and acceleration of phenotypic evolution in spiny rats (Trinomys,Echimyidae) across different environments

Trinomys is a very diverse mammalian genus of eastern Brazil, comprising ten species. Eight of them have remained in their original biome, the Atlantic Forest, while two others occupy the semi‐arid Caatinga. One of the species from the Atlantic Forest, T. eliasi, also inhabits adjacent restingas, an open vegetation ecosystem on sandy sediments of marine origin. Here we assess the impact of colonization of very different environments on the rates of phenotypic change in Trinomys. Multivariate analyses showed that species that had remained in the Atlantic Forest maintained a strong cohesion in skull shape and size. By contrast, colonization of the Caatinga was associated with remarkable differentiation and a strong increase in the rates of phenotypic change, 7–17 times higher than in the Atlantic Forest. Trinomys eliasi from restinga differed from forest samples and other species by a recent change in morphology. Our findings showed that occupation of remarkably different ecosystems had a significant impact in the phenotypic evolution of Trinomys.     

Phylogeography of Liquidambar styraciflua (Altingiaceae) in Mesoamerica: survivors of a Neogene widespread temperate forest (or cloud forest) in North America?

We investigate the genetic variation between populations of the American sweetgum (Liquidambar styraciflua), a tree species with a disjunct distribution between northeastern Texas and Mexico, by analyzing sequences of two chloroplast DNA plastid regions in Mesoamerica. Our results revealed phylogeographical structure, with private haplotypes distributed in unique environmental space at either side of the Trans‐Mexican Volcanic Belt, and a split in the absence of gene flow dating back ca. 4.2–1.4 million years ago (MYA). Species distribution modeling results fit a model of refugia along the Gulf and Atlantic coasts but the present ranges of US and Mesoamerican populations persisted disjunct during glacial/interglacial cycles. Divergence between the US and Mesoamerican (ca. 8.4–2.8 MYA) populations of L. styraciflua and asymmetrical gene flow patterns support the hypothesis of a long‐distance dispersal during the Pliocene, with fragmentation since the most recent glacial advance (120,000 years BP) according to coalescent simulations and high effective migration rates from Mesoamerica to the USA and close to zero in the opposite direction. Our findings implicate the Trans‐Mexican Volcanic Belt as a porous barrier driving genetic divergence of L. styraciflua, corresponding with environmental niche differences, during the Pliocene to Quaternary volcanic arc episode 3.6 MYA, and a Mesoamerican origin of populations in the USA.     

Genetic diversity of two mitochondrial DNA genes in Spirometra erinaceieuropaei (Cestoda: Diphyllobothridae) from Poland

The tapeworm species Spirometra erinaceieuropaei was documented mainly in Asia and Europe. In recent years, plerocercoid larvae (spargana) of this parasite have been found in different hosts in north‐eastern Poland. The evolutionary history and way of S. erinaceieuropaei spreading across Eurasia have been not described yet. However, this phenomenon could be closely related to the evolutionary history and migration routes of studied tapeworm host species. We investigated the genetic variability and divergence pattern among S. erinaceieuropaei populations in intermediate and paratenic hosts from north‐eastern Poland based on complete mitochondrial sequences of cytochrome b (cytb) and cytochrome c oxidase subunit I (cox1) genes. Analysis of 319 consolidated sequences of these two genes showed no genetic structure across study area. Comparison of sequences from Poland and China showed distinct separation of S. erinaceieuropaei populations from these two regions. They split from their common ancestor approximately 28.6 million years ago. Demographic expansion of Polish population of S. erinaceieuropaei started from glacial refugia approximately 12.5 thousand years ago, and recent population expansion has been observed in the tapeworm population from north‐eastern Poland.     

Ecological regime shift drives declining growth rates of sea turtles throughout the West Atlantic

Somatic growth is an integrated, individual‐based response to environmental conditions, especially in ectotherms. Growth dynamics of large, mobile animals are particularly useful as bio‐indicators of environmental change at regional scales. We assembled growth rate data from throughout the West Atlantic for green turtles, Chelonia mydas, which are long‐lived, highly migratory, primarily herbivorous mega‐consumers that may migrate over hundreds to thousands of kilometers. Our dataset, the largest ever compiled for sea turtles, has 9690 growth increments from 30 sites from Bermuda to Uruguay from 1973 to 2015. Using generalized additive mixed models, we evaluated covariates that could affect growth rates; body size, diet, and year have significant effects on growth. Growth increases in early years until 1999, then declines by 26% to 2015. The temporal (year) effect is of particular interest because two carnivorous species of sea turtles—hawksbills, Eretmochelys imbricata, and loggerheads, Caretta caretta—exhibited similar significant declines in growth rates starting in 1997 in the West Atlantic, based on previous studies. These synchronous declines in productivity among three sea turtle species across a trophic spectrum provide strong evidence that an ecological regime shift (ERS) in the Atlantic is driving growth dynamics. The ERS resulted from a synergy of the 1997/1998 El Niño Southern Oscillation (ENSO)—the strongest on record—combined with an unprecedented warming rate over the last two to three decades. Further support is provided by the strong correlations between annualized mean growth rates of green turtles and both sea surface temperatures (SST) in the West Atlantic for years of declining growth rates (r = −.94) and the Multivariate ENSO Index (MEI) for all years (r = .74). Granger‐causality analysis also supports the latter finding. We discuss multiple stressors that could reinforce and prolong the effect of the ERS. This study demonstrates the importance of region‐wide collaborations.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.     

Predicting current and future global distributions of whale sharks

The Vulnerable (IUCN) whale shark spans warm and temperate waters around the globe. However, their present‐day and possible future global distribution has never been predicted. Using 30 years (1980–2010) of whale shark observations recorded by tuna purse‐seiners fishing in the Atlantic, Indian and Pacific Oceans, we applied generalized linear mixed‐effects models to test the hypothesis that similar environmental covariates predict whale shark occurrence in all major ocean basins. We derived global predictors from satellite images for chlorophyll a and sea surface temperature, and bathymetric charts for depth, bottom slope and distance to shore. We randomly generated pseudo‐absences within the area covered by the fisheries, and included fishing effort as an offset to account for potential sampling bias. We predicted sea surface temperatures for 2070 using an ensemble of five global circulation models under a no climate‐policy reference scenario, and used these to predict changes in distribution. The full model (excluding standard deviation of sea surface temperature) had the highest relative statistical support (wAIC_c = 0.99) and explained ca. 60% of the deviance. Habitat suitability was mainly driven by spatial variation in bathymetry and sea surface temperature among oceans, although these effects differed slightly among oceans. Predicted changes in sea surface temperature resulted in a slight shift of suitable habitat towards the poles in both the Atlantic and Indian Oceans (ca. 5°N and 3–8°S, respectively) accompanied by an overall range contraction (2.5–7.4% and 1.1–6.3%, respectively). Predicted changes in the Pacific Ocean were small. Assuming that whale shark environmental requirements and human disturbances (i.e. no stabilization of greenhouse gas emissions) remain similar, we show that warming sea surface temperatures might promote a net retreat from current aggregation areas and an overall redistribution of the species.     

Late Quaternary climatic vegetational shifts in an ecological transition zone of northern Madagascar: insights from genetic analyses of two endemic rodent species

The Loky‐Manambato region, located in northern Madagascar, is a biotically rich contact zone between different forest biomes. Local current forest cover is composed of both humid and dry formations, which show elevational stratification. A recent phylogeographical study of a regional dry forest rodent, Eliurus carletoni (subfamily Nesomyinae), found genetic evidence of forest contractions between 18 750 and 7500 years BP, which based on extrapolation of the pollen subfossil record, was thought to be associated with an expansion of local humid forests. Herein, we conduct a genetic test of this hypothesis and focused on populations on two neighbouring massifs of forest‐dependent rodent species, one associated with low‐elevation dry forests (E. carletoni) and the other with higher elevation humid forests (Eliurus tanala). Using mitochondrial markers and a combination of traditional and coalescent‐based phylogeographical, historical demographic and population genetic methods, we found evidence of historical connections between populations of E. tanala. Adjacent populations of E. carletoni and E. tanala exhibit opposite historical demographic patterns, and for both, evidence suggests that historical demographic events occurred within the last 25 000 years BP. These findings strongly support the proposed late Quaternary shifts in the floristic composition of the Loky‐Manambato region.     

Collections from the mesophytic zone off Bermuda reveal three species of Kallymeniaceae (Gigartinales,Rhodophyta) in genera with transoceanic distributions

A molecular survey of red algae collected by technical divers and submersibles from 90 m in the mesophotic zone off the coast of Bermuda revealed three species assignable to the Kallymeniaceae. Two of the species are representative of recently described genera centered in the western Pacific in Australia and New Zealand, Austrokallymenia and Psaromenia and the third from the Mediterranean Sea and the eastern Atlantic, Nothokallymenia. A phylogenetic analysis of concatenated mitochondrial (COI‐5P) and chloroplast (rbcL) genes, as well as morphological characteristics, revealed that two are shown to be new species with distant closest relatives (N. erosa and Psaromenia septentrionalis), while the third represents a deep water western Atlantic species now moved to an Australasian genus (A. westii).     

Contrasting responses to Pleistocene climate changes: a case study of two sister species Allium cyathophorum and A. spicata (Amaryllidaceae) distributed in the eastern and western Qinghai–Tibet Plateau

It has been hypothesized that species occurring in the eastern and the western Qinghai–Tibet Plateau (QTP) responded differently to climate changes during the Pleistocene. Here, we test this hypothesis by phylogeographic analysis of two sister species, Allium cyathophorum and A. spicata. We sequenced two chloroplast DNA (cpDNA) fragments (accD‐psaI and the rpl16 intron) of 150 individuals, and the nuclear (ITS) region of 114 individuals, from 19 populations throughout the distributional ranges of these species. The divergence between the two species was dated at 779 ‐ 714 thousand years before the present and was likely initiated by the most major glaciation in the QTP. Analysis of chlorotype diversity showed that A. spicata, the species occurring in the western QTP, contains much lower genetic diversity (0.25) than A. cyathophorum (0.93), which is distributed in the eastern QTP. Moreover, multiple independent tests suggested that the A. spicata population had expanded recently, while no such expansion was detected in A. cyathophorum, indicating a contrasting pattern of responses to Pleistocene climate changes. These findings highlight the importance of geographical topography in determining how species responded to the climate changes that took place in the QTP during the Pleistocene.     

Multispecies genetic structure and hybridization in the Betula genus across Eurasia

Boreal and cool temperate forests are the major land cover of northern Eurasia, and information about continental‐scale genetic structure and past demographic history of forest species is important from an evolutionary perspective and has conservation implications. However, although many population genetic studies of forest tree species have been conducted in Europe or Eastern Asia, continental‐scale genetic structure and past demographic history remain poorly known. Here, we focus on the birch genus Betula, which is commonly distributed in boreal and cool temperate forests, and examine 129 populations of two tetraploid and four diploid species collected from Iceland to Japan. All individuals were genotyped at seven to 18 nuclear simple sequence repeats (nSSRs). Pairwise among the six species ranged from 0.285 to 0.903, and genetic differentiation among them was clear. structure analysis suggested that Betula pubescens is an allotetraploid and one of the parental species was Betula pendula. In both species pairs of B. pendula and B. plathyphylla, and B. pubescens and B. ermanii, genetic diversity was highest in central Siberia. A hybrid zone was detected around Lake Baikal for eastern and western species pairs regardless of ploidy level. Approximate Bayesian computation suggested that the divergence of B. pendula and B. platyphylla occurred around the beginning of the last ice age (36 300 years BP, 95% CI: 15 330–92 700) and hybridization between them was inferred to have occurred after the last glacial maximum (1614 years BP, 95% CI: 561–4710), with B. pendula providing a higher contribution to hybrids.     

Uncommon southwest swells trigger sea urchin disease outbreaks in Eastern Atlantic archipelagos

Recurrent sea urchin mass mortality has recently affected eastern Atlantic populations of the barren‐forming sea urchin Diadema africanum. This new episode of die‐off affords the opportunity to determine common meteorological and oceanographic conditions that may promote disease outbreaks. The population dynamics of this sea urchin species are well known—urchin barrens have persisted for many decades along most of the coastlines off the archipelagos of Madeira, Selvages, and the Canary Islands, where they limit macroalgae biomass growth. However, this new and explosive mortality event decimated the sea urchin population by 93% on Tenerife and La Palma Islands. Two severe episodes of southwestern rough sea that led to winter storms, in February 2010 (Xynthia) and February 2018 (Emma), preceded both mass mortality events. The autumn and winter months of those years were anomalous and characterized by swells with an average wave height above 2 m that hit the south and southwest sides of the islands. The amoeba Paramoeba brachiphila was the only pathogen isolated this time from the moribund and dead sea urchins, suggesting that the amoeba was the primary cause of the mortality. This new sea urchin die‐off event supports the “killer‐storm” hypothesis that has been already described for western Atlantic coasts. These anomalous southwest storms during winters generate pronounced underwater sediment movement and large‐scale vertical mixing, detected in local tide gauge, which may promote paramoebiasis. This study presents valuable insights about climate‐mediated changes in disease frequency and its impacts on the future of coastal marine ecosystems in the Atlantic.     

Genomic replacement of native Cobitis lutheri with introduced C. tetralineata through a hybrid swarm following the artificial connection of river systems

River connections via artificial canals will bring about secondary contacts between previously isolated fish species. Here, we present a genetic consequence of such a secondary contact between Cobitis fish species, C. lutheri in the Dongjin River, and C. tetralineata in the Seomjin River in Korea. The construction of water canals about 80 years ago has unidirectionally introduced C. tetralineata into the native habitat of C. lutheri, and then these species have hybridized in the main stream section of the Dongjin River. According to the divergence population genetic analyses of DNA sequence data, the two species diverged about 3.3 million years ago, which is interestingly coincident with the unprecedented paleoceanographic change that caused isolations of the paleo‐river systems in northeast Asia due to sea‐level changes around the late Pliocene. Multilocus genotypic data of nine microsatellites and three nuclear loci revealed an extensively admixed structure in the hybrid zone with a high proportion of various post‐F1 hybrids. Surprisingly, pure native C. lutheri was absent in the hybrid zone in contrast to the 7% of pure C. tetralineata. Such a biased proportion must have resulted from the dominant influence of continually introducing C. tetralineata on the native C. lutheri which has no supply of natives from other tributaries to the hybrid zone due to numerous low‐head dams. In addition, mating experiments indicated that there is no discernible reproductive isolation between them. All the results suggest that the gene pool of native C. lutheri is being rapidly replaced by that of continually introducing C. tetralineata through a hybrid swarm for the last 80 years, which will ultimately lead to the genomic extinction of natives in this hybrid zone.     

Life history tactics shape amphibians’ demographic responses to the North Atlantic Oscillation

Over the last three decades, climate abnormalities have been reported to be involved in biodiversity decline by affecting population dynamics. A growing number of studies have shown that the North Atlantic Oscillation (NAO) influences the demographic parameters of a wide range of plant and animal taxa in different ways. Life history theory could help to understand these different demographic responses to the NAO. Indeed, theory states that the impact of weather variation on a species’ demographic traits should depend on its position along the fast–slow continuum. In particular, it is expected that NAO would have a higher impact on recruitment than on adult survival in slow species, while the opposite pattern is expected occur in fast species. To test these predictions, we used long‐term capture–recapture datasets (more than 15,000 individuals marked from 1965 to 2015) on different surveyed populations of three amphibian species in Western Europe: Triturus cristatus, Bombina variegata, and Salamandra salamandra. Despite substantial intraspecific variation, our study revealed that these three species differ in their position on a slow–fast gradient of pace of life. Our results also suggest that the differences in life history tactics influence amphibian responses to NAO fluctuations: Adult survival was most affected by the NAO in the species with the fastest pace of life (T. cristatus), whereas recruitment was most impacted in species with a slower pace of life (B. variegata and S. salamandra). In the context of climate change, our findings suggest that the capacity of organisms to deal with future changes in NAO values could be closely linked to their position on the fast–slow continuum.     

Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species

Prioritizing and making efficient conservation plans for threatened populations requires information at both evolutionary and ecological timescales. Nevertheless, few studies integrate multidisciplinary approaches, mainly because of the difficulty for conservationists to assess simultaneously the evolutionary and ecological status of populations. Here, we sought to demonstrate how combining genetic and demographic analyses allows prioritizing and initiating conservation plans. To do so, we combined snapshot microsatellite data and a 30‐year‐long demographic survey on a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale. Our results revealed low levels of genetic diversity and weak effective population sizes (<63 individuals) in all populations. We further detected severe bottlenecks dating back to the last centuries (200–800 years ago), which may explain the differentiation of certain populations. The demographic survey revealed a general decrease in the spatial distribution and abundance of P. toxostoma over the last three decades. We conclude that demo‐genetic approaches are essential for (1) identifying populations for which both evolutionary and ecological extinction risks are high; and (2) proposing conservation plans targeted toward these at risk populations, and accounting for the evolutionary history of populations. We suggest that demo‐genetic approaches should be the norm in conservation practices.