Genetic diversity from pre-bottleneck to recovery in two sympatric pinniped species in the Northwest Atlantic

Conservation successes of the past several decades provide natural settings to study post-bottleneck evolutionary processes in species undergoing recovery. Here, we study the impact of demographic change on genetic diversity in parallel natural experiments of historical decline and subsequent recovery in two sympatric pinniped species in the Northwest Atlantic, the gray seal (Halichoerus grypus atlantica) and harbor seal (Phoca vitulina concolor). We compare genetic diversity at the mitochondrial control region today to diversity in archaeological specimens, which represent the populations prior to the regional bounties of the late 1800s to mid-1900s that drastically reduced population sizes and led to local extirpations. We further assess genetic diversity throughout recovery, using biological collections from ongoing long-term studies of both species. Overall, the genetic data are consistent with the historical presence of large, genetically diverse populations of pinnipeds prior to human exploitation, and suggest that gray seals were more dramatically impacted by historical bottlenecks than harbor seals in the Northwest Atlantic. Current mitochondrial diversity in both species is relatively high, and we observe little change over the past several decades during a period of roughly parallel rapid population increases. However, there remain large differences in haplotype composition between pinniped populations of pre-exploitation and today, a lasting genetic signature of historical exploitation that is likely to persist into the future.     

A hitchhiker’s guide to the Maritimes: anthropogenic transport facilitates long‐distance dispersal of an invasive marine crab to Newfoundland

Aim To determine timing, source and vector for the recent introduction of the European green crab, Carcinus maenas (Linnaeus, 1758), to Newfoundland using multiple lines of evidence. Location Founding populations in Placentia Bay, Newfoundland, Canada and potential source populations in the north‐west Atlantic (NWA) and Europe. Methods We analysed mitochondrial and microsatellite genetic data from European and NWA populations sampled during 1999–2002 to determine probable source locations and vectors for the Placentia Bay introduction discovered in 2007. We also analysed Placentia Bay demographic data and shipping records to look for congruent patterns with genetic analyses. Results Demographic data and surveys suggested that C. maenas populations are established and were in Placentia Bay for several years (c. 2002) prior to discovery. Genetic data corroboratively suggested central/western Scotian Shelf populations (e.g., Halifax) as the likely source area for the anthropogenic introduction. These Scotian Shelf populations were within an admixture zone made up of genotypes from both the earlier (early 1800s) and later (late 1900s) introductions of the crab to the NWA from Europe. Placentia Bay also exhibited this mixed ancestry. Probable introduction vectors included vessel traffic and shipping, especially vessels carrying ballast water. Main conclusions Carcinus maenas overcame considerable natural barriers (i.e., coastal and ocean currents) via anthropogenic transport to become established and abundant in Newfoundland. Our study thus demonstrates how non‐native populations can be important secondary sources of introduction especially when aided by human transport. Inference of source populations was possible owing to the existence of an admixture zone in central/western Nova Scotia made up of southern and northern genotypes corresponding with the crab’s two historical introductions. Coastal vessel traffic was found to be a likely vector for the crab’s spread to Newfoundland. Our study demonstrates that there is considerable risk for continued introduction or reintroduction of C. maenas throughout the NWA.     

Life-stage-dependent supergene haplotype frequencies and metapopulation neutral genetic patterns of Atlantic cod,Gadus morhua,from Canada's Northern cod stock region and adjacent areas

Among highly migratory fish species, nursery areas occupied by juveniles often differ from adult habitats. To better understand the spatial dynamics of Canada's Northern cod stock, juveniles caught off the east coast of Newfoundland and Labrador were compared to adults from the same region as well as individuals from other areas in Atlantic Canada using double-digest restriction site-associated DNA sequencing–derived single nucleotide polymorphisms. A reduced proportion of homozygotes with a chromosomal inversion located in linkage group 1 (LG1) was detected between juvenile and adult samples in the Northern cod stock region, potentially indicating age-dependent habitat use or ontogenetic selection for attributes associated with the many genes located in LG1. No selectively neutral genetic differences were found between samples from the Northern cod stock; nevertheless, significant differences were found between some of these samples and cod collected from St. Pierre Bank, Bay of Fundy, Browns Bank and the southern Scotian Shelf. Clustering analysis of variants at neutral loci provided evidence for three major genetic units: (a) the Newfoundland Atlantic Coast, (b) eastern and southern Gulf of St. Lawrence and Burgeo Bank and (c) the Bay of Fundy, Browns Bank and southern Scotian Shelf. Both adaptive and neutral population structure within the Northern cod stock should be considered by managers to promote demographic rebuilding of the stock.     

Microsatellite variation in cod Gadus morhua throughout its geographic range

Previous genetic studies using neutral markers such as allozymes, mtDNA and minisatellite loci have demonstrated varying amounts of population structure in cod Gadus morhua throughout the Atlantic. Microsatellite loci, which are potentially the most informative of presently available neutral genetic markers, have been applied extensively within western and eastern Atlantic areas but not on a range-wide basis. In the present study, six microsatellite DNA loci were used to screen cod samples from nine locations throughout the geographic range from the Scotian Shelf in the West Atlantic to the Barents and Baltic Seas in the east. Overall F _ST value was 0·03 ( P = < 0·001) across all samples. Statistically significant population differences over all loci combined were evident between more geographically distant samples, using either heterogeneity tests or F _ST analysis, with at least one locus showing significant differences between all samples (prior to Bonferroni correction). A significant correlation was observed between genetic and geographical distance, suggesting a higher level of historical and contemporary gene flow between adjacent populations than more distant populations. Samples from either end of the geographic range (Scotian Shelf and Baltic Sea) were particularly distinct when analysed using the STRUCTURE programme and also showed a high level of self-assignment when individuals of either the Scotian Shelf or Baltic Sea were tested against the entire data set. The present microsatellite study demonstrates a high level of geographic population structure between the western Atlantic, middle and eastern Atlantic and Baltic Sea, and thus, the findings should be useful in devising overall management and conservation strategies for the species.     

Recent Demographic History and Present Fine-Scale Structure in the Northwest Atlantic Leatherback (Dermochelys coriacea) Turtle Population

The leatherback turtle Dermochelys coriacea is the most widely distributed sea turtle species in the world. It exhibits complex life traits: female homing and migration, migrations of juveniles and males that remain poorly known, and a strong climatic influence on resources, breeding success and sex-ratio. It is consequently challenging to understand population dynamics. Leatherbacks are critically endangered, yet the group from the Northwest Atlantic is currently considered to be under lower risk than other populations while hosting some of the largest rookeries. Here, we investigated the genetic diversity and the demographic history of contrasted rookeries from this group, namely two large nesting populations in French Guiana, and a smaller one in the French West Indies. We used 10 microsatellite loci, of which four are newly isolated, and mitochondrial DNA sequences of the control region and cytochrome b. Both mitochondrial and nuclear markers revealed that the Northwest Atlantic stock of leatherbacks derives from a single ancestral origin, but show current genetic structuration at the scale of nesting sites, with the maintenance of migrants amongst rookeries. Low nuclear genetic diversities are related to founder effects that followed consequent bottlenecks during the late Pleistocene/Holocene. Most probably in response to climatic oscillations, with a possible influence of early human hunting, female effective population sizes collapsed from 2 million to 200. Evidence of founder effects and high numbers of migrants make it possible to reconsider the population dynamics of the species, formerly considered as a metapopulation model: we propose a more relaxed island model, which we expect to be a key element in the currently observed recovering of populations. Although these Northwest Atlantic rookeries should be considered as a single evolutionary unit, we stress that local conservation efforts remain necessary since each nesting site hosts part of the genetic diversity and species history.     

Genomic signatures of population bottleneck and recovery in Northwest Atlantic pinnipeds

Population increases over the past several decades provide natural settings in which to study the evolutionary processes that occur during bottleneck, growth, and spatial expansion. We used parallel natural experiments of historical decline and subsequent recovery in two sympatric pinniped species in the Northwest Atlantic, the gray seal (Halichoerus grypus atlantica) and harbor seal (Phoca vitulina vitulina), to study the impact of recent demographic change in genomic diversity. Using restriction site‐associated DNA sequencing, we assessed genomic diversity at over 8,700 polymorphic gray seal loci and 3,700 polymorphic harbor seal loci in samples from multiple cohorts collected throughout recovery over the past half‐century. Despite significant differences in the degree of genetic diversity assessed in the two species, we found signatures of historical bottlenecks in the contemporary genomes of both gray and harbor seals. We evaluated temporal trends in diversity across cohorts, as well as compared samples from sites at both the center and edge of a recent gray seal range expansion, but found no significant change in genomewide diversity following recovery. We did, however, find that the variance and degree of allele frequency change measured over the past several decades were significantly different from neutral expectations of drift under population growth. These two cases of well‐described demographic history provide opportunities for critical evaluation of current approaches to simulating and understanding the genetic effects of historical demographic change in natural populations.     

Genetic Structure of Capelin (Mallotus villosus) in the Northwest Atlantic Ocean

Capelin (Mallotus villosus) is a commercially exploited, key forage-fish species found in the boreal waters of the North Pacific and North Atlantic Oceans. We examined the population structure of capelin throughout their range in the Canadian northwest Atlantic Ocean using genetic-based methods. Capelin collected at ten beach and five demersal spawning locations over the period 2002 through 2008 (N = 3,433 fish) were genotyped using six polymorphic microsatellite loci. Temporally distinct samples were identified at three beach spawning locations: Chance Cove, Little Lawn and Straitsview, Newfoundland. Four capelin stocks are assumed for fisheries management in the northwest Atlantic Ocean based on meristics, morphometrics, tag returns, and seasonal distribution patterns. Our results suggested groupings that were somewhat different than the assumed structure, and indicate at least seven genetically defined populations arising from two ancestral populations. The spatial mosaic of capelin from each of the two basal cluster groups explains much of the observed geographic variability amongst neighbouring samples. The genetic-defined populations were resolved at Jost’s D _est ≥ 0.01 and were composed of fish collected 1) in the Gulf of St. Lawrence, 2) along the south and east coasts of Newfoundland, 3) along coastal northern Newfoundland and southern Labrador, 4) along coastal northern Labrador, 5) near the Saguenay River, and at two nearshore demersal spawning sites, 6) one at Grebes Nest off Bellevue Beach on the east coast of Newfoundland, and 7) one off the coast of Labrador at Domino Run. Moreover, the offshore demersal spawners on the Scotian Shelf and Southeast Shoal appeared to be related to the inshore demersal spawners at Grebes Nest and in Domino Run and to beach spawners from the Gulf of St. Lawrence.     

A nuclear DNA basis for shelf- and bank-scale population structure in northwest Atlantic cod (Gadus morhua): Labrador to Georges Bank

Variation at five microsatellite DNA loci scored in ≈ 1300 individuals provided evidence of genetic structure among 14 cod ( Gadus morhua ) populations spanning the range of the species in the northwest (NW) Atlantic. Using D _A and D _SW measures of genetic distance, as well as F _ST and R _ST measures of population structure, differences are revealed among populations at continental shelf scales (NE Newfoundland Shelf, Grand Banks, Flemish Cap, Scotian Shelf, Georges Bank) where regions are separated by submarine saddles, channels and trenches. However, we also provide evidence of genetic structure at spawning-bank scales consistent with variation in oceanographic features and in the spatiotemporal distribution of spawning, each of which may represent barriers to gene flow among geographically contiguous populations inhabiting a highly advective environment. The differences described are consistent with postdispersal spawning fidelity to natal areas, a behaviour that may be facilitated by topographically induced gyre-like circulations that can act as retention mechanisms. Significant degrees of substructure among neighbouring and contiguous cod populations may be most easily explained by the associated oceanographic features and processes that conceivably form the template for the evolution of the structure. We suggest that bathymetric and hydrodynamic structure represents a rational starting point for developing hypotheses to examine the processes that lead to the genetic structuring of marine fish species.     

Bay-scale population structure in coastal Atlantic cod in Labrador and Newfoundland, Canada

Polymorphisms at five microsatellite DNA loci provide evidence that Atlantic cod Gadus morhua inhabiting Gilbert Bay, Labrador are genetically distinguishable from offshore cod on the north-east Newfoundland shelf and from inshore cod in Trinity Bay, Newfoundland. Antifreeze activity in the blood suggests that Gilbert Bay cod overwinter within the Bay. Gilbert Bay cod are also smaller (weight and length) for their age and consequently less fecund for their age, than cod elsewhere within the northern cod complex. The productivity and recruitment potential of coastal cod off Labrador may thus be much lower than that of offshore northern cod or of inshore cod farther south, implying that a more conservative management strategy may be required for cod from coastal Labrador than traditionally practised for northern cod inhabiting less harsh environments. Relatively high F _ST and R _ST measures of population structure suggest that important barriers to gene flow exist among five components that include two inshore (Gilbert and Trinity Bay) and three offshore cod aggregations on the north-east Newfoundland Shelf and the Grand Bank. D _A and D _SW estimates of genetic distance that involve Gilbert Bay cod are approximately three- and 10–fold larger, respectively, than estimates not involving Gilbert Bay cod. The differences between inshore cod from Gilbert Bay and Trinity Bay raise the possibility that other genetically distinguishable coastal populations may exist, or may have existed prior to the northern cod fishery collapse. Harvesting strategies for northern cod should recognize the existence of genetic diversity between inshore and offshore components as well as among coastal components.     

The Potential for Spatial Distribution Indices to Signal Thresholds in Marine Fish Biomass

The frequently observed positive relationship between fish population abundance and spatial distribution suggests that changes in distribution can be indicative of trends in abundance. If contractions in spatial distribution precede declines in spawning stock biomass (SSB), spatial distribution reference points could complement the SSB reference points that are commonly used in marine conservation biology and fisheries management. When relevant spatial distribution information is integrated into fisheries management and recovery plans, risks and uncertainties associated with a plan based solely on the SSB criterion would be reduced. To assess the added value of spatial distribution data, we examine the relationship between SSB and four metrics of spatial distribution intended to reflect changes in population range, concentration, and density for 10 demersal populations (9 species) inhabiting the Scotian Shelf, Northwest Atlantic. Our primary purpose is to assess their potential to serve as indices of SSB, using fisheries independent survey data. We find that metrics of density offer the best correlate of spawner biomass. A decline in the frequency of encountering high density areas is associated with, and in a few cases preceded by, rapid declines in SSB in 6 of 10 populations. Density-based indices have considerable potential to serve both as an indicator of SSB and as spatially based reference points in fisheries management.     

Population genetic structure of the endangered crayfish Austropotamobius pallipes in France based on microsatellite variation: biogeographical inferences and conservation implications

1. One important goal in conservation biology is to characterise evolutionary lineages within endangered species before management decisions are taken. Here, we assess population differentiation in the freshwater crayfish Austropotamobius pallipes, an endangered species endemic to western Europe and provide valuable information for the conservation of French populations. 2. Analysis of five microsatellite loci in 44 populations revealed very different within population levels of genetic diversity (0.000 < H₀ < 0.564). Two groups, corresponding to northern and southern French populations, showed a high degree of genetic differentiation in both allele frequencies and allele sizes. Comparison of these results with previous studies of A. pallipes strongly suggests that the divergence between northern and southern populations could have occurred during the last glaciation period of the Pleistocene from one Atlantic and one Mediterranean refuge. 3. Evidence for genetic admixture between these two lineages was revealed by correspondence analyses in southern populations, probably as the result of artificial translocations. 4. French populations appeared significantly differentiated among the different river drainages and were highly structured within rivers. The impact of population size, population bottlenecks and founder events on the population genetic differentiation are discussed. 5. Based on these results, we propose the designation of two evolutionarily significant units for A. pallipes in France. Our data also support the maintenance of separate demographic management strategies for crayfish inhabiting different river systems. However, genetic analyses will have to be combined with demographic and ecological data for sustainable conservation programmes.     

Chronicle of an extinction foretold: genetic properties of an extremely small population of Iberolacerta monticola

The fragmentation and destruction of natural habitats by human intervention is producing a continuous and inexorable reduction of the size of populations in multitude of species all over the world. Small and isolated populations face higher extinction risks, due to demographic and environmental stochasticity, and also because of several genetic threats, among which inbreeding is considered the most important one. For many of these species, the extinction of a population is an irreversible event, so that determining the immediate importance of these risk factors and understanding their interactions is crucial for conservation plans. Iberolacerta monticola is a small lacertid endemic to the northwestern Iberian Peninsula, distributed mainly across moderate/high altitude mountainous regions. Some populations are found nearly at sea-level, though, in fluvial valleys with relict Atlantic forests, in the severely fragmented western part of its range. One of them has been dramatically reduced over the last 30 years, and presently is on the brink of extinction. Using microsatellite nuclear markers, we obtained different measures of genetic variation at this site, together with demographic and breeding data. Both the level of heterozygosity and the number of alleles per locus indicate that the level of variation in this population is comparatively high, and the average inbreeding coefficient is very low. Individuals appear healthy and long-lived, and are related by a few different lines of descent. These findings are discussed in the context of current theories and experimental evidence of associative overdominance and purging of the genetic load of populations, with special emphasis on the evolutionary potential of recovery of small evolutionary units.     

Genome‐scale sampling suggests cryptic epigenetic structuring and insular divergence in Canada lynx

Determining the molecular signatures of adaptive differentiation is a fundamental component of evolutionary biology. A key challenge is to identify such signatures in wild organisms, particularly between populations of highly mobile species that undergo substantial gene flow. The Canada lynx (Lynx canadensis) is one species where mainland populations appear largely undifferentiated at traditional genetic markers, despite inhabiting diverse environments and displaying phenotypic variation. Here, we used high‐throughput sequencing to investigate both neutral genetic structure and epigenetic differentiation across the distributional range of Canada lynx. Newfoundland lynx were identified as the most differentiated population at neutral genetic markers, with demographic modelling suggesting that divergence from the mainland occurred at the end of the last glaciation (20–33 KYA). In contrast, epigenetic structure revealed hidden levels of differentiation across the range coincident with environmental determinants including winter conditions, particularly in the peripheral Newfoundland and Alaskan populations. Several biological pathways related to morphology were differentially methylated between populations, suggesting that epigenetic modifications might explain morphological differences seen between geographically peripheral populations. Our results indicate that epigenetic modifications, specifically DNA methylation, are powerful markers to investigate population differentiation in wild and non‐model systems.     

Mitochondrial phylogeography and demographic history of the vicuña: implications for conservation

The vicu?a (Vicugna vicugna; Miller, 1924) is a conservation success story, having recovered from near extinction in the 1960s to current population levels estimated at 275,000. However, lack of information about its demographic history and genetic diversity has limited both our understanding of its recovery and the development of science-based conservation measures. To examine the evolution and recent demographic history of the vicu?a across its current range and to assess its genetic variation and population structure, we sequenced mitochondrial DNA from the control region (CR) for 261 individuals from 29 populations across Peru, Chile and Argentina. Our results suggest that populations currently designated as Vicugna vicugna vicugna and Vicugna vicugna mensalis comprise separate mitochondrial lineages. The current population distribution appears to be the result of a recent demographic expansion associated with the last major glacial event of the Pleistocene in the northern (18 to 22 degrees S) dry Andes 14-12,000 years ago and the establishment of an extremely arid belt known as the 'Dry Diagonal' to 29 degrees S. Within the Dry Diagonal, small populations of V. v. vicugna appear to have survived showing the genetic signature of demographic isolation, whereas to the north V. v. mensalis populations underwent a rapid demographic expansion before recent anthropogenic impacts.     

Selection maintains MHC diversity through a natural population bottleneck

A perceived consequence of a population bottleneck is the erosion of genetic diversity and concomitant reduction in individual fitness and evolutionary potential. Although reduced genetic variation associated with demographic perturbation has been amply demonstrated for neutral molecular markers, the effective management of genetic resources in natural populations is hindered by a lack of understanding of how adaptive genetic variation will respond to population fluctuations, given these are affected by selection as well as drift. Here, we demonstrate that selection counters drift to maintain polymorphism at a major histocompatibility complex (MHC) locus through a population bottleneck in an inbred island population of water voles. Before and after the bottleneck, MHC allele frequencies were close to balancing selection equilibrium but became skewed by drift when the population size was critically low. MHC heterozygosity generally conformed to Hardy-Weinberg expectations except in one generation during the population recovery where there was a significant excess of heterozygous genotypes, which simulations ascribed to strong differential MHC-dependent survival. Low allelic diversity and highly skewed frequency distributions at microsatellite loci indicated potent genetic drift due to a strong founder affect and/or previous population bottlenecks. This study is a real-time examination of the predictions of fundamental evolutionary theory in low genetic diversity situations. The findings highlight that conservation efforts to maintain the genetic health and evolutionary potential of natural populations should consider the genetic basis for fitness-related traits, and how such adaptive genetic diversity will vary in response to both the demographic fluctuations and the effects of selection.     

Heterogeneity within the native range: population genetic analyses of sympatric invasive and noninvasive clades of the freshwater invading copepod Eurytemora affinis

Invasive species are often composed of highly differentiated populations or sibling species distributed across their native ranges. This study analysed patterns of distribution and the evolutionary and demographic histories of populations within the native range of the copepod species complex Eurytemora affinis. Genetic structure was analysed for samples from 17 locations from both the invaded and native ranges in the St Lawrence River drainage basin, using 652 base pairs of the mitochondrial cytochrome oxidase subunit I gene. This study revealed a high degree of heterogeneity in genetic structure and habitat type in the native range, as well as a bias in the sources of invasive populations. Two genetically distinct clades showed a pattern of niche partitioning within the St Lawrence basin. The noninvasive North Atlantic clade primarily occupied the central portion of the St Lawrence Middle Estuary, whereas the invasive Atlantic clade was more prevalent along the margins, in the upstream reaches of the estuary and downstream salt marshes. Habitat partitioning and genetic subdivision was also present within the Atlantic clade. The freshwater populations were genetically more proximate to the Atlantic clade populations in the estuary than to those in the salt marsh, suggesting the estuary as the source of the invasive populations. The freshwater invading populations showed evidence of a modest population bottleneck. Populations from both clades showed genetic signatures of demographic population expansions that preceded the timing of the last glacial maximum, supporting the St Lawrence as a secondary contact zone between the two clades. Additional analyses on physiological and evolutionary properties of populations in the native range, along with analysis of the selection regime within native habitats, might yield insights into the evolutionary potential to invade.     

Panmictic population structure in the hooded seal (Cystophora cristata)

Two putative populations of hooded seals (Cystophora cristata) occur in the North Atlantic. The Greenland Sea population pup and breed on the pack ice near Jan Mayen ('West Ice') while the Northwest Atlantic population is thought to pup in the Davis Strait, in the Gulf of St. Lawrence (the 'Gulf'), and off southern Labrador or northeast Newfoundland (the 'Front'). We used microsatellite profiling of 300 individuals at 13 loci and mitochondrial DNA sequencing of the control region of 123 individuals to test for genetic differentiation between these four breeding herds. We found no significant genetic differences between breeding areas, nor evidence for cryptic nor higher level genetic structure in this species. The Greenland Sea breeding herd was genetically most distant from the Northwest Atlantic breeding areas; however, the differences were statistically nonsignificant. Our data therefore suggest that the world's hooded seals comprise a single panmictic genetic population.     

Northern areas as refugia for temperate species under current climate warming: Atlantic salmon (Salmo salar L.) as a model in northern Europe

In this work, patterns of geographical genetic diversity in Atlantic salmon Salmo salar were studied across the whole Atlantic Arc; whether these patterns (and thus genetic population structure) were affected by water temperatures was also evaluated. Salmo salar populations were characterized using microsatellite loci and then analysed with reference to ocean surface temperature data from across the region. Analysis showed the presence of a latitudinal cline of genetic variability (higher in northern areas) and water temperatures (sea surface temperatures) determining genetic population structure (the latter in combination with genetic drift in southern populations). Under the current global change scenario, northern areas of Europe would constitute refugia for diversity in the future. This is effectively the inverse of what appears to have happened in glacial refugia during the last glacial maximum. From this perspective, the still abundant and large northern populations S. salar should be considered as precious as the small almost relict southern ones and given appropriate protection. Careful management of the species, coordinated across countries and latitudes, is needed in order to avoid its extinction in Europe.     

Phylogeography and genetic differentiation along the distributional range of the orchid Epidendrum fulgens: a Neotropical coastal species not restricted to glacial refugia

Aim Phylogeographical studies in the Brazilian Atlantic Forest (BAF) have mostly included species associated with forest habitats, whereas taxa associated with grassland and sand‐dune plant communities have so far been largely overlooked. This study examines the phylogeography of the orchid Epidendrum fulgens, which occurs on coastal sand dunes and granitic outcrops, in order to identify major genetic divergences or disjunctions across the range of the species and to investigate the genetic signatures of past range contractions and expansions. Location Southern and south‐eastern seashore vegetation along the BAF biome, and granitic and arenitic outcrops that occur in the subtropical grassland plant communities located south of the BAF. Methods Nine nuclear and four plastid microsatellite loci were used to genotype 424 individuals from 16 populations across the distributional range of E. fulgens. For both sets of markers, we estimated genetic diversity and population differentiation, testing for a north–south gradient of genetic diversity. The plastid haplotype network and a Bayesian assignment analysis of nuclear markers were used to infer population structure. Past demographic changes were investigated using a coalescence approach. Results A deep disjunction was found between northern populations within the BAF and southern populations outside the BAF that occur on granitic and arenitic outcrops. Recent demographic reductions were detected in northern populations on coastal sands. Such demographic changes were not expected for those populations, as previous studies with forest species had found evidence of population expansion in the same areas. Higher genetic diversity was found in southern populations on granite, in contrast to patterns observed in previous studies of forest species. Main conclusions The results are consistent with the long‐term persistence of E. fulgens. Bottlenecks were detected in populations from areas where population expansion events have been detected in other plant (and animal) species, suggesting that forest expansion after the Last Glacial Maximum played a role in the population fragmentation and decrease in genetic diversity in E. fulgens. A substantial genetic division in E. fulgens corresponds to the ‘Portal de Torres’, a region that demarcates the northern limits of subtropical grassland plant communities and the southern limits of the BAF.     

Impacts of the Toba eruption and montane forest expansion on diversification in Sumatran parachuting frogs (Rhacophorus)

Catastrophic events, such as volcanic eruptions, can have profound impacts on the demographic histories of resident taxa. Due to its presumed effect on biodiversity, the Pleistocene eruption of super‐volcano Toba has received abundant attention. We test the effects of the Toba eruption on the diversification, genetic diversity, and demography of three co‐distributed species of parachuting frogs (Genus Rhacophorus) on Sumatra. We generate target‐capture data (~950 loci and ~440,000 bp) for three species of parachuting frogs and use these data paired with previously generated double digest restriction‐site associated DNA (ddRADseq) data to estimate population structure and genetic diversity, to test for population size changes using demographic modelling, and to estimate the temporal clustering of size change events using a full‐likelihood Bayesian method. We find that populations around Toba exhibit reduced genetic diversity compared with southern populations, and that northern populations exhibit a shift in effective population size around the time of the eruption (~80 kya). However, we infer a stronger signal of expansion in southern populations around ~400 kya, and at least two of the northern populations may have also expanded at this time. Taken together, these findings suggest that the Toba eruption precipitated population declines in northern populations, but that the demographic history of these three species was also strongly impacted by mid‐Pleistocene forest expansion during glacial periods. We propose local rather than regional effects of the Toba eruption, and emphasize the dynamic nature of diversification on the Sunda Shelf.