Limits of Bayesian skyline plot analysis of mtDNA sequences to infer historical demographies in Pacific herring (and other species)

Species-Specific Traits Shape Genetic Diversity During an Expansion-Contraction Cycle and Bias Demographic History Reconstruction

Species ranges are dynamic, experiencing expansions, contractions or shifts in response to habitat changes driven by extrinsic factors such as climate change or human activities. While existing research examines the genetic consequences of spatial processes, few studies integrate species-specific traits to analyse how habitat changes affect co-existing species. In this study, we address this gap by investigating how genetic diversity patterns vary among species with different traits (such as generation length, population density and dispersal) experiencing similar habitat changes. Using spatial simulations and a simpler panmictic population model, we investigate the temporal genetic diversity in refugium populations undergoing range expansion of their habitat, followed by stationary and contraction periods. By varying habitat contraction speed and species traits, we identified three distinct temporal dynamics of genetic diversity during contraction: (i) a decrease in genetic diversity, (ii) an initial increase followed by a decrease and (iii) a continuous increase throughout the contraction period. We show that genetic diversity trajectories during population decline can be predicted by comparing sampled population diversity to equilibrium values expected under expanded and contracted habitat ranges. Our study also challenges the belief that high genetic diversity in a refugium population is due to a recent and rapid habitat loss. Instead, we found contrasting effects of contraction speed on genetic diversity depending on the interaction between species-specific traits and the dynamics of habitat change. Finally, using simulated genetic data, we found that demographic histories inferred from effective population size estimates may vary across species, even when they experience similar habitat changes.

Genetic variation of endangered Jankowski’s Bunting (Emberiza jankowskii): High connectivity and a moderate history of demographic decline

Introduction

Continued discovery of “mismatch” patterns between population size and genetic diversity, involving wild species such as insects, amphibians, birds, mammals, and others, has raised issues about how population history, especially recent dynamics under human disturbance, affects currently standing genetic variation. Previous studies have revealed high genetic diversity in endangered Jankowski’s Bunting. However, it is unclear how the demographic history and recent habitat changes shape the genetic variation of Jankowski’s Bunting.

Methods

To explore the formation and maintenance of high genetic diversity in endangered Jankowski’s Bunting, we used a mitochondrial control region (partial mtDNA CR) and 15 nuclear microsatellite markers to explore the recent demographic history of Jankowski’s Bunting, and we compared the historical and contemporary gene flows between populations to reveal the impact of habitat change on population connectivity. Specifically, we aimed to test the following hypotheses: (1) Jankowski’s Bunting has a large historical Ne and a moderate demographic history; and (2) recent habitat change might have no significant impact on the species’ population connectivity.

Results

The results suggested that large historical effective population size, as well as severe but slow population decline, may partially explain the high observable genetic diversity. Comparison of historical (over the past 4Ne generations) and contemporary (1–3 generations) gene flow indicated that the connectivity between five local populations was only marginally affected by landscape changes.

Discussion

Our results suggest that high population connectivity and a moderate history of demographic decline are powerful explanations for the rich genetic variation in Jankowski’s Bunting. Although there is no evidence that the genetic health of Jankowski’s Bunting is threatened, the time-lag effects on the genetic response to recent environmental changes is a reminder to be cautious about the current genetic characteristics of this species. Where possible, factors influencing genetic variation should be integrated into a systematic framework for conducting robust population health assessments. Given the small contemporary population size, inbreeding, and ecological specialization, we recommend that habitat protection be maintained to maximize the genetic diversity and population connectivity of Jankowski’s Bunting.

Bayesian Skyline Plots disagree with range size changes based on Species Distribution Models for Holarctic birds

During the Quaternary, large climate oscillations impacted the distribution and demography of species globally. Two approaches have played a major role in reconstructing changes through time: Bayesian Skyline Plots (BSPs), which reconstruct population fluctuations based on genetic data, and Species Distribution Models (SDMs), which allow us to back-cast the range occupied by a species based on its climatic preferences. In this paper, we contrast these two approaches by applying them to a large data set of 102 Holarctic bird species, for which both mitochondrial DNA sequences and distribution maps are available, to reconstruct their dynamics since the Last Glacial Maximum (LGM). Most species experienced an increase in effective population size (N_e , as estimated by BSPs) as well as an increase in geographical range (as reconstructed by SDMs) since the LGM; however, we found no correlation between the magnitude of changes in N_e and range size. The only clear signal we could detect was a later and greater increase in N_e for wetland birds compared to species that live in other habitats, a probable consequence of a delayed and more extensive increase in the extent of this habitat type after the LGM. The lack of correlation between SDM and BSP reconstructions could not be reconciled even when range shifts were considered. We suggest that this pattern might be linked to changes in population densities, which can be independent of range changes, and caution that interpreting either SDMs or BSPs independently is problematic and potentially misleading.

Late Pleistocene Expansion of Small Murid Rodents across the Palearctic in Relation to the Past Environmental Changes

We investigated the evolutionary history of the striped field mouse to identify factors that initiated its past demographic changes and to shed light on the causes of its current genetic structure and trans-Eurasian distribution. We sequenced mitochondrial cyt b from 184 individuals, obtained from 35 sites in central Europe and eastern Mongolia. We compared genetic analyses with previously published historical distribution models and data on environmental and climatic changes. The past demographic changes displayed similar population trends in the case of recently expanded clades C1 and C3, with the glacial (MIS 3–4) expansion and postglacial bottleneck preceding the recent expansion initiated in the late Holocene and were related to environmental changes during the upper Pleistocene and Holocene. The past demographic trends of the eastern Asian clade C3 were correlated with changes in sea level and the formation of new land bridges formed by the exposed sea shelf during the glaciations. These data were supported by reconstructed historical distribution models. The results of our genetic analyses, supported by the reconstruction of the historical spatial distributions of the distinct clades, confirm that over time the local populations mixed as a consequence of environmental and climatic changes resulting from cyclical glaciation and the interglacial period during the Pleistocene.

mtDNAcombine: tools to combine sequences from multiple studies

BACKGROUND

Today an unprecedented amount of genetic sequence data is stored in publicly available repositories. For decades now, mitochondrial DNA (mtDNA) has been the workhorse of genetic studies, and as a result, there is a large volume of mtDNA data available in these repositories for a wide range of species. Indeed, whilst whole genome sequencing is an exciting prospect for the future, for most non-model organisms' classical markers such as mtDNA remain widely used. By compiling existing data from multiple original studies, it is possible to build powerful new datasets capable of exploring many questions in ecology, evolution and conservation biology. One key question that these data can help inform is what happened in a species' demographic past. However, compiling data in this manner is not trivial, there are many complexities associated with data extraction, data quality and data handling.

RESULTS

Here we present the mtDNAcombine package, a collection of tools developed to manage some of the major decisions associated with handling multi-study sequence data with a particular focus on preparing sequence data for Bayesian skyline plot demographic reconstructions.

CONCLUSIONS

There is now more genetic information available than ever before and large meta-data sets offer great opportunities to explore new and exciting avenues of research. However, compiling multi-study datasets still remains a technically challenging prospect. The mtDNAcombine package provides a pipeline to streamline the process of downloading, curating, and analysing sequence data, guiding the process of compiling data sets from the online database GenBank.

Chaotic genetic structure and past demographic expansion of the invasive gastropod Tritia neritea in its native range, the Mediterranean Sea

To better predict population evolution of invasive species in introduced areas it is critical to identify and understand the mechanisms driving genetic diversity and structure in their native range. Here, we combined analyses of the mitochondrial COI gene and 11 microsatellite markers to investigate both past demographic history and contemporaneous genetic structure in the native area of the gastropod Tritia neritea, using Bayesian skyline plots (BSP), multivariate analyses and Bayesian clustering. The BSP framework revealed population expansions, dated after the last glacial maximum. The haplotype network revealed a strong geographic clustering. Multivariate analyses and Bayesian clustering highlighted the strong genetic structure at all scales, between the Black Sea and the Adriatic Sea, but also within basins. Within basins, a random pattern of genetic patchiness was observed, suggesting a superimposition of processes involving natural biological effects (no larval phase and thus limited larval dispersal) and putative anthropogenic transport of specimens. Contrary to the introduced area, no isolation-by-distance patterns were recovered in the Mediterranean or the Black Seas, highlighting different mechanisms at play on both native and introduced areas, triggering unknown consequences for species’ evolutionary trajectories. These results of Tritia neritea populations on its native range highlight a mixture of ancient and recent processes, with the effects of paleoclimates and life history traits likely tangled with the effects of human-mediated dispersal.

Comparative phylogeography of two Agarophyton species in the New Zealand archipelago

Molecular studies have reported the coexistence of two species of Agarophyton in New Zealand: the newly described A.transtasmanicum with an apparently restricted distribution to some sites in the North Island, and the more widespread A.chilense. Here, we compared the distribution, genetic diversity, and structure of both Agarophyton species throughout the archipelago using sequences of the nuclear Internal Transcribed Spacer 2 (ITS2) marker. Agarophyton chilense's distribution was continuous and extensive along the North and South Islands, Stewart Island, and Chatham Island, and the genetic clusters were mostly concordant with boundaries between biogeographic regions. In contrast, specimens of A.transtasmanicum were collected in four sites broadly distributed in both the North and South Islands, with no clear spatial structure of the genetic diversity. Populations, where the species co-occurred, tended to display similar levels in genetic diversity for the two species. Demographic inferences supported a postglacial demographic expansion for two A.chilense genetic clusters, one present in the South Island and the eastern coast of the North Island, and the other present in northern South Island. A third genetic cluster located on the western coast of the North Island had a signature of long-term demographic stability. For A.transtasmanicum, the skyline plot also suggested a postglacial demographic expansion. Last, we developed a new molecular tool to quickly and easily distinguish between the two Agarophyton species, which could be used to ease future fine-scale population studies, especially in areas where the two species coexist.

Against all odds: a tale of marine range expansion with maintenance of extremely high genetic diversity

The displacement of species from equatorial latitudes to temperate locations following the increase in sea surface temperatures is among the significant reported consequences of climate change. Shifts in the distributional ranges of species result in fish communities tropicalisation, i.e., high latitude colonisations by typically low latitude distribution species. These movements create new interactions between species and new trophic assemblages. The Senegal seabream, Diplodus bellottii, may be used as a model to understand the population genetics of these invasions. In the last decades, this species has undergone an outstanding range expansion from its African area of origin to the Atlantic coast of the Iberian Peninsula, where now occurs abundantly. Mitochondrial and nuclear markers revealed a striking high haplotypic nucleotide and genetic diversity values, along with significant population differentiation throughout the present-day geographical range of the Senegal seabream. These results are not consistent with the central-marginal hypothesis, nor with the expectations of a leptokurtic distribution of individuals, as D. bellottii seems to be able to retain exceptional levels of diversity in marginal and recently colonised areas. We discuss possible causes for hyperdiversity and lack of geographical structure and subsequent implications for fisheries.

Reconstructing the population history of the sandy beach amphipod Haustorioides japonicus using the calibration of demographic transition (CDT) approach

Calibration of the molecular rate is one of the major challenges in marine population genetics. Although the use of an appropriate evolutionary rate is crucial in exploring population histories, calibration of the rate is always difficult because fossil records and geological events are rarely applicable for rate calibration. The acceleration of the evolutionary rate for recent coalescent events (or more simply, the time dependency of the molecular clock) is also a problem that can lead to overestimation of population parameters. Calibration of demographic transition (CDT) is a rate calibration technique that assumes a post-glacial demographic expansion, representing one of the most promising approaches for dealing with these potential problems in the rate calibration. Here, we demonstrate the importance of using an appropriate evolutionary rate, and the power of CDT, by using populations of the sandy beach amphipod Haustorioides japonicus along the Japanese coast of the northwestern Pacific Ocean. Analysis of mitochondrial sequences found that the most peripheral population in the Pacific coast of northeastern Honshu Island (Tohoku region) is genetically distinct from the other northwestern Pacific populations. By using the two-epoch demographic model and rate of temperature change, the evolutionary rate was modeled as a log-normal distribution with a median rate of 2.2%/My. The split-time of the Tohoku population was subsequently estimated to be during the previous interglacial period by using the rate distribution, which enables us to infer potential causes of the divergence between local populations along the continuous Pacific coast of Japan.

The Gambian epauletted fruit bat shows increased genetic divergence in the Ethiopian highlands and in an area of rapid urbanization

The Gambian epauletted fruit bat (Epomophorus gambianus) is an abundant species that roosts in both urban and rural settings. The possible role of E. gambianus as a reservoir host of zoonotic diseases underlines the need to better understand the species movement patterns. So far, neither observational nor phylogenetic studies have identified the dispersal range or behavior of this species. Comparative analyses of mitochondrial and nuclear markers from 20 localities across the known distribution of E. gambianus showed population panmixia, except for the populations in Ethiopia and southern Ghana (Accra and Ve-Golokwati). The Ethiopian population may be ancestral and is highly divergent to the species across the rest of its range, possibly reflecting isolation of an ancient colonization along an east-west axis. Mitochondrial haplotypes in the Accra population display a strong signature of a past bottleneck event; evidence of either an ancient or recent bottleneck using microsatellite data, however, was not detected. Demographic analyses identified population expansion in most of the colonies, except in the female line of descent in the Accra population. The molecular analyses of the colonies from Ethiopia and southern Ghana show gender dispersal bias, with the mitochondrial DNA fixation values over ten times those of the nuclear markers. These findings indicate free mixing of the species across great distances, which should inform future epidemiological studies.

Riverine barrier effects on population genetic structure of the Hanuman langur (Semnopithecus entellus) in the Nepal Himalaya

BACKGROUND

Past climatological events and contemporary geophysical barriers shape the distribution, population genetic structure, and evolutionary history of many organisms. The Himalayan region, frequently referred to as the third pole of the Earth, has experienced large-scale climatic oscillations in the past and bears unique geographic, topographic, and climatic areas. The influences of the Pleistocene climatic fluctuations and present-day geographical barriers such as rivers in shaping the demographic history and population genetic structure of organisms in the Nepal Himalaya have not yet been documented. Hence, we examined the effects of late-Quaternary glacial-interglacial cycles and riverine barriers on the genetic composition of Hanuman langurs (Semnopithecus entellus), a colobine primate with a wide range of altitudinal distribution across the Nepalese Himalaya, using the mitochondrial DNA control region (CR, 1090 bp) and cytochrome B (CYTB, 1140 bp) sequences combined with paleodistribution modeling.

RESULTS

DNA sequences were successfully retrieved from 67 non-invasively collected fecal samples belonging to 18 wild Hanuman langur troops covering the entire distribution range of the species in Nepal. We identified 37 haplotypes from the concatenated CR + CYTB (2230 bp) sequences, with haplotype and nucleotide diversities of 0.958 ± 0.015 and 0.0237 ± 0.0008, respectively. The troops were clustered into six major clades corresponding to their river-isolated spatial distribution, with the significantly high genetic variation among these clades confirming the barrier effects of the snow-fed Himalayan rivers on genetic structuring. Analysis of demographic history projected a decrease in population size with the onset of the last glacial maximum (LGM); and, in accordance with the molecular analyses, paleodistribution modeling revealed a range shift in its suitable habitat downward/southward during the LGM. The complex genetic structure among the populations of central Nepal, and the stable optimal habitat through the last interglacial period to the present suggest that the central mid-hills of Nepal served as glacial refugia for the Hanuman langur.

CONCLUSIONS

Hanuman langurs of the Nepal Himalaya region exhibit high genetic diversity, with their population genetic structure is strongly shaped by riverine barrier effects beyond isolation by distance; hence, this species demands detailed future phylogenetic study.

Natural history and molecular evolution of demersal Mediterranean sharks and skates inferred by comparative phylogeographic and demographic analyses

Background

The unique and complex paleoclimatic and paleogeographic events which affected the Mediterranean Sea since late Miocene deeply influenced the distribution and evolution of marine organisms and shaped their genetic structure. Following the Messinian salinity crisis and the sea-level fluctuations during the Pleistocene, several Mediterranean marine species developed deep genetic differentiation, and some underwent rapid radiation. Here, we consider two of the most prioritized groups for conservation in the light of their evolutionary history: sharks and rays (elasmobranchs). This paper deals with a comparative multispecies analysis of phylogeographic structure and historical demography in two pairs of sympatric, phylogenetically- and ecologically-related elasmobranchs, two scyliorhinid catsharks (Galeus melastomus, Scyliorhinus canicula) and two rajid skates (Raja clavata, Raja miraletus). Sampling and experimental analyses were designed to primarily test if the Sicilian Channel can be considered as effective eco-physiological barrier for Mediterranean demersal sympatric elasmobranchs.

Methods

The phylogeography and the historical demography of target species were inferred by analysing the nucleotide variation of three mitochondrial DNA markers (i.e., partial sequence of COI, NADH2 and CR) obtained from a total of 248 individuals sampled in the Western and Eastern Mediterranean Sea as well as in the adjacent northeastern Atlantic Ocean. Phylogeographic analysis was performed by haplotype networking and testing spatial genetic differentiation of samples (i.e., analysis of molecular variance and of principal components). Demographic history of Mediterranean populations was reconstructed using mismatch distribution and Bayesian Skyline Plot analyses.

Results

No spatial genetic differentiation was identified in either catshark species, while phylogeographic structure of lineages was identified in both skates, with R. miraletus more structured than R. clavata. However, such structuring of skate lineages was not consistent with the separation between Western and Eastern Mediterranean. Sudden demographic expansions occurred synchronously during the upper Pleistocene (40,000–60,000 years ago) in both skates and G. melastomus, likely related to optimal environmental conditions. In contrast, S. canicula experienced a slow and constant increase in population size over the last 350,000 years.

Discussion

The comparative analysis of phylogeographic and historical demographic patterns for the Mediterranean populations of these elasmobranchs reveals that historical phylogeographic breaks have not had a large impact on their microevolution. We hypothesize that interactions between environmental and ecological/physiological traits may have been the driving force in the microevolution of these demersal elasmobranch species in the Mediterranean rather than oceanographic barriers.

Population genetic structure and evolutionary history of Bale monkeys (Chlorocebus djamdjamensis) in the southern Ethiopian Highlands

Background

Species with a restricted geographic distribution, and highly specialized habitat and dietary requirements, are particularly vulnerable to extinction. The Bale monkey (Chlorocebus djamdjamensis) is a little-known arboreal, bamboo-specialist primate endemic to the southern Ethiopian Highlands. While most Bale monkeys inhabit montane forests dominated by bamboo, some occupy forest fragments where bamboo is much less abundant. We used mitochondrial DNA (mtDNA) sequences to analyse the genetic structure and evolutionary history of Bale monkeys covering the majority of their remaining distribution range. We analysed 119 faecal samples from their two main habitats, continuous forest (CF) and fragmented forests (FF), and sequenced 735 bp of the hypervariable region I (HVI) of the control region. We added 12 orthologous sequences from congeneric vervets (C. pygerythrus) and grivets (C. aethiops) as well as animals identified as hybrids, previously collected in southern Ethiopia.

Results

We found strong genetic differentiation (with no shared mtDNA haplotypes) between Bale monkey populations from CF and FF. Phylogenetic analyses revealed two distinct and highly diverged clades: a Bale monkey clade containing only Bale monkeys from CF and a green monkey clade where Bale monkeys from FF cluster with grivets and vervets. Analyses of demographic history revealed that Bale monkey populations (CF and FF) have had stable population sizes over an extended period, but have all recently experienced population declines.

Conclusions

The pronounced genetic structure and deep mtDNA divergence between Bale monkey populations inhabiting CF and FF are likely to be the results of hybridization and introgression of the FF population with parapatric Chlorocebus species, in contrast to the CF population, which was most likely not impacted by hybridization. Hybridization in the FF population was probably enhanced by an alteration of the bamboo forest habitat towards a more open woodland habitat, which enabled the parapatric Chlorocebus species to invade the Bale monkey's range and introgress the FF population. We therefore propose that the CF and FF Bale monkey populations should be managed as separate units when developing conservation strategies for this threatened species.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1217-y) contains supplementary material, which is available to authorized users.

High genetic diversity and low-population differentiation in the Patagonian sprat (Sprattus fuegensis) based on mitochondrial DNA

The Patagonian sprat, Sprattus fuegensis, is a small pelagic marine fish that inhabits the continental shelf along the coasts of Chilean Patagonian and Argentina, a distribution that was highly impacted during the Last Glacial Maximum (LGM). In order to identify how the LGM played a role on the current observed genetic diversity and population structure of S. fuegensis, we analyzed 1438 nucleotide positions from the control region of 335 individuals collected at 12 sites across its distribution. Genetic diversity and differentiation indices were calculated to identify population structure, and a Bayesian skyride plot (BSRP) reconstruction was carried out to infer the historic population dynamics. Extremely high genetic diversity was found at all locations analyzed, non-population structure was found across its distribution, and the BSRP showed two increases in effective population size over time. Our outcomes suggest that the current genetic diversity, population structure and population expansion may have occurred during the medium and late Pleistocene.

Population panmixia and demographic expansion of a highly piscivorous marine fish Scomberomorus niphonius

Population structure and demographic history of the Japanese Spanish mackerel Scomberomorus niphonius a highly piscivorous and migratory marine fish, were assessed using mitochondrial DNA control region sequences (n = 720) and microsatellite genotypes at five loci (n = 1331) for samples collected on Japanese coasts from 2001 to 2010. The population structure was panmictic and the haplotype and allele frequencies were temporally stable even during the recent recovery process. Demographic expansion was strongly supported throughout the Pleistocene, suggesting that the oscillating glacial and interglacial climate conditions in the Pleistocene had no substantial impact on the demographic history of S. niphonius.

Identification and analysis of mtDNA genomes attributed to Finns reveal long-stagnant demographic trends obscured in the total diversity

In Europe, modern mitochondrial diversity is relatively homogeneous and suggests an ubiquitous rapid population growth since the Neolithic revolution. Similar patterns also have been observed in mitochondrial control region data in Finland, which contrasts with the distinctive autosomal and Y-chromosomal diversity among Finns. A different picture emerges from the 843 whole mitochondrial genomes from modern Finns analyzed here. Up to one third of the subhaplogroups can be considered as Finn-characteristic, i.e. rather common in Finland but virtually absent or rare elsewhere in Europe. Bayesian phylogenetic analyses suggest that most of these attributed Finnish lineages date back to around 3,000–5,000 years, coinciding with the arrival of Corded Ware culture and agriculture into Finland. Bayesian estimation of past effective population sizes reveals two differing demographic histories: 1) the ‘local’ Finnish mtDNA haplotypes yielding small and dwindling size estimates for most of the past; and 2) the ‘immigrant’ haplotypes showing growth typical of most European populations. The results based on the local diversity are more in line with that known about Finns from other studies, e.g., Y-chromosome analyses and archaeology findings. The mitochondrial gene pool thus may contain signals of local population history that cannot be readily deduced from the total diversity.

Sea-level variations have influenced the demographic history of estuarine and freshwater fishes of the coastal plain of Paraná, Brazil

This study surveyed the mitochondrial haplotype diversity of nine freshwater fish species and two estuarine-marine species from the coastal basins and drainages of the highland plateaus of Paraná, Brazil. Portions of the cytochrome b gene or the control region were sequenced. The demographic history of each species was inferred using the Bayesian skyline method, mismatch distribution analysis and statistical neutrality tests. Demographic reconstruction analyses revealed a single pattern of variation in the effective population size (N_e ) among species. No dramatic changes in N_e were detected in upland species. By contrast, evidence of population expansion over the past 200 000 years was detected in all coastal plain and estuarine species. These findings correspond to periods of low sea-level (regressions) followed by a rapid increase in the sea-level by >100 m. The resulting reconnections and subsequent fragmentation and isolation between the estuarine and freshwater bodies were putatively relevant to the historical demography of the fish species in these areas.

Genetic diversity and population structure of Hyporhamphus sajori (Beloniformes: Hemiramphidae) inferred from mtDNA control region and msDNA markers

This paper presents preliminary data on the genetic diversity and population structure of Hyporhamphus sajori by analysing a 510 bp sequence in the mitochondrial DNA (mtDNA) control region and eight polymorphic microsatellite DNA loci. The H. sajori individuals from different locations were indistinguishable from one another based on mtDNA variation, as demonstrated with a neighbour-joining tree and minimum spanning network analysis. Low level of genetic diversity and the absence of population structure in H. sajori from the north-west Pacific Ocean, combined with negative indices for neutral evolution in these populations, suggest that H. sajori underwent a population expansion after a recent bottleneck. The Structure analysis, discriminant analysis of principal components (DAPC) and the pair-wise Φ_ST values after Bonferroni correction using eight microsatellite loci provided no clear inference on the genetic differentiation and thus no evidence of population structure of H. sajori. The genetic connectivity among locations might be due to fairly high gene flow via transport of eggs and larvae by the Kuroshio and Tsushima warm current. This study revealed low levels of genetic diversity and suggested high level of contemporary gene flow among populations of H. sajori in the East (Japan) Sea and the Pacific Ocean.

Phylogeography of the Koala, (Phascolarctos cinereus), and Harmonising Data to Inform Conservation

The Australian continent exhibits complex biogeographic patterns but studies of the impacts of Pleistocene climatic oscillation on the mesic environments of the Southern Hemisphere are limited. The koala (Phascolarctos cinereus), one of Australia’s most iconic species, was historically widely distributed throughout much of eastern Australia but currently represents a complex conservation challenge. To better understand the challenges to koala genetic health, we assessed the phylogeographic history of the koala. Variation in the maternally inherited mitochondrial DNA (mtDNA) Control Region (CR) was examined in 662 koalas sampled throughout their distribution. In addition, koala CR haplotypes accessioned to Genbank were evaluated and consolidated. A total of 53 unique CR haplotypes have been isolated from koalas to date (including 15 haplotypes novel to this study). The relationships among koala CR haplotypes were indicative of a single Evolutionary Significant Unit and do not support the recognition of subspecies, but were separated into four weakly differentiated lineages which correspond to three geographic clusters: a central lineage, a southern lineage and two northern lineages co-occurring north of Brisbane. The three geographic clusters were separated by known Pleistocene biogeographic barriers: the Brisbane River Valley and Clarence River Valley, although there was evidence of mixing amongst clusters. While there is evidence for historical connectivity, current koala populations exhibit greater structure, suggesting habitat fragmentation may have restricted female-mediated gene flow. Since mtDNA data informs conservation planning, we provide a summary of existing CR haplotypes, standardise nomenclature and make recommendations for future studies to harmonise existing datasets. This holistic approach is critical to ensuring management is effective and small scale local population studies can be integrated into a wider species context.

Panmixia in a Critically Endangered Fish: The Totoaba (Totoaba macdonaldi) in the Gulf of California

Conservation of the evolutionary legacy of endangered species is a key component for long-term persistence. Totoaba is a long-lived fish endemic to the Gulf of California and is considered critically endangered. There is currently a debate concerning its conservation status and whether it can be used as a fishery resource. Unfortunately, basic information on biological and genetic population structure of the species is lacking. We sampled 313 individuals and employed 16 microsatellite loci and 3 mitochondrial DNA markers (16S, 547 pb; COI, 619 pb; control region, 650 pb) to assess population structure and demography of totoaba in the Gulf of California, with samples from locations that encompass nearly all of its recognized geographic distribution. We could not reject a hypothesis of panmixia for totoaba, using nuclear or mitochondrial markers. Demographic analysis of mtDNA suggests a sudden population expansion model. The results have important implications for totoaba conservation because poaching is a significant conservation challenge and could have additive negative effects over the single population of totoaba in the Gulf of California.

Does asymmetric gene flow among matrilines maintain the evolutionary potential of the European eel?

Using evolutionary theory to predict the dynamics of populations is one of the aims of evolutionary conservation. In endangered species, with geographic range extending over continuous areas, the predictive capacity of evolutionary-based conservation measures greatly depends on the accurate identification of reproductive units. The endangered European eel (Anguilla anguilla) is a highly migratory fish species with declining population due to a steep recruitment collapse in the beginning of the 1980s. Despite punctual observations of genetic structure, the population is viewed as a single panmictic reproductive unit. To understand the possible origin of the detected structure in this species, we used a combination of mitochondrial and nuclear loci to indirectly evaluate the possible existence of cryptic demes. For that, 403 glass eels from three successive cohorts arriving at a single location were screened for phenotypic and genetic diversity, while controlling for possible geographic variation. Over the 3 years of sampling, we consistently identified three major matrilines which we hypothesized to represent demes. Interestingly, not only we found that population genetic models support the existence of those matriline-driven demes over a completely panmictic mode of reproduction, but also we found evidence for asymmetric gene flow amongst those demes. We uphold the suggestion that the detection of demes related to those matrilines reflect a fragmented spawning ground, a conceptually plausible consequence of the low abundance that the European eel has been experiencing for three decades. Furthermore, we suggest that this cryptic organization may contribute to the maintenance of the adaptive potential of the species.

Effects of geological changes and climatic fluctuations on the demographic histories and low genetic diversity of Squaliobarbus curriculus in Yellow River

The 104 samples of Squaliobarbus curriculus were collected from four localities in Yellow River and one region in Yangtze River. Analyses of the first hypervariable region of mitochondrial DNA control region of 555bp revealed only 15 polymorphism sites and defined 19 haplotypes. Low-to-moderate levels of haplotype diversity and low nucleotide diversity were observed in Yellow River populations (h=0.2529-0.7510, π=0.0712%-0.2197%). In contrast, Poyang Lake population showed high haplotype diversity and lower-middle nucleotide diversity (h=0.9636, π=0.5317%). Low genetic differentiation was estimated among Yellow River populations and significant level of genetic structure was detected between two rivers. Population genetic structure between two rivers was believed to be connected with geographical barriers and paleoclimatic events. The demographic history of S. curriculus in Yellow River examined by neutrality tests, mismatch distribution analysis, and Bayesian skyline analysis suggested a sudden and spatial population expansion dating to the Holocene. Climatic warming and changes of Yellow River course may have important effects on demographic facet of S. curriculus history. The same signal was also obtained on Poyang Lake population in late Pleistocene during the last interglacial period. During the period, the pronounced climatic change and the water system variation of PYL may have an important influence on the population.

Contemporary genetic structure and postglacial demographic history of the black scorpionfish, Scorpaena porcus, in the Mediterranean and the Black Seas

Understanding the distribution of genetic diversity in the light of past demographic events linked with climatic shifts will help to forecast evolutionary trajectories of ecosystems within the current context of climate change. In this study, mitochondrial sequences and microsatellite loci were analysed using traditional population genetic approaches together with Bayesian dating and the more recent approximate Bayesian computation scenario testing. The genetic structure and demographic history of a commercial fish, the black scorpionfish, Scorpaena porcus, was investigated throughout the Mediterranean and Black Seas. The results suggest that the species recently underwent population expansions, in both seas, likely concomitant with the warming period following the Last Glacial Maximum, 20 000 years ago. A weak contemporaneous genetic differentiation was identified between the Black Sea and the Mediterranean Sea. However, the genetic diversity was similar for populations of the two seas, suggesting a high number of colonizers entered the Black Sea during the interglacial period and/or the presence of a refugial population in the Black Sea during the glacial period. Finally, within seas, an east/west genetic differentiation in the Adriatic seems to prevail, whereas the Black Sea does not show any structured spatial genetic pattern of its population. Overall, these results suggest that the Black Sea is not that isolated from the Mediterranean, and both seas revealed similar evolutionary patterns related to climate change and changes in sea level.

Red River barrier and Pleistocene climatic fluctuations shaped the genetic structure of Microhyla fissipes complex (Anura: Microhylidae) in southern China and Indochina

South China and Indochina host striking species diversity and endemism. Complex tectonic and climatic evolutions appear to be the main drivers of the biogeographic patterns. In this study, based on the geologic history of this region, we test 2 hypotheses using the evolutionary history of Microhyla fissipes species complex. Using DNA sequence data from both mitochondrial and nuclear genes, we first test the hypothesis that the Red River is a barrier to gene flow and dispersal. Second, we test the hypothesis that Pleistocene climatic cycling affected the genetic structure and population history of these frogs. We detect 2 major genetic splits that associate with the Red River. Time estimation suggests that late Miocene tectonic movement associated with the Red River drove their diversification. Species distribution modeling (SDM) resolves significant ecological differences between sides of the Red River. Thus, ecological divergence also probably promoted and maintained the diversification. Genogeography, historical demography, and SDM associate patterns in southern China with climate changes of the last glacial maximum (LGM), but not Indochina. Differences in geography and climate between the 2 areas best explain the discovery. Responses to the Pleistocene glacial–interglacial cycling vary among species and regions.

Late Glacial Demographic Expansion Motivates a Clock Overhaul for Population Genetics

The molecular clock hypothesis is fundamental in evolutionary biology as by assuming constancy of the molecular rate it provides a timeframe for evolution. However, increasing evidence shows time dependence of inferred molecular rates with inflated values obtained using recent calibrations. As recent demographic calibrations are virtually non-existent in most species, older phylogenetic calibration points (>1 Ma) are commonly used, which overestimate demographic parameters. To obtain more reliable rates of molecular evolution for population studies, I propose the calibration of demographic transition (CDT) method, which uses the timing of climatic changes over the late glacial warming period to calibrate expansions in various species. Simulation approaches and empirical data sets from a diversity of species (from mollusk to humans) confirm that, when compared with other genealogy-based calibration methods, the CDT provides a robust and broadly applicable clock for population genetics. The resulting CDT rates of molecular evolution also confirm rate heterogeneity over time and among taxa. Comparisons of expansion dates with ecological evidence confirm the inaccuracy of phylogenetically derived divergence rates when dating population-level events. The CDT method opens opportunities for addressing issues such as demographic responses to past climate change and the origin of rate heterogeneity related to taxa, genes, time, and genetic information content.

Mitogenomic Meta-Analysis Identifies Two Phases of Migration in the History of Eastern Eurasian Sheep

Despite much attention, history of sheep (Ovis aries) evolution, including its dating, demographic trajectory and geographic spread, remains controversial. To address these questions, we generated 45 complete and 875 partial mitogenomic sequences, and performed a meta-analysis of these and published ovine mitochondrial DNA sequences (n = 3,229) across Eurasia. We inferred that O. orientalis and O. musimon share the most recent female ancestor with O. aries at approximately 0.790 Ma (95% CI: 0.637-0.934 Ma) during the Middle Pleistocene, substantially predating the domestication event (∼8-11 ka). By reconstructing historical variations in effective population size, we found evidence of a rapid population increase approximately 20-60 ka, immediately before the Last Glacial Maximum. Analyses of lineage expansions showed two sheep migratory waves at approximately 4.5-6.8 ka (lineages A and B: ∼6.4-6.8 ka; C: ∼4.5 ka) across eastern Eurasia, which could have been influenced by prehistoric West-East commercial trade and deliberate mating of domestic and wild sheep, respectively. A continent-scale examination of lineage diversity and approximate Bayesian computation analyses indicated that the Mongolian Plateau region was a secondary center of dispersal, acting as a "transportation hub" in eastern Eurasia: Sheep from the Middle Eastern domestication center were inferred to have migrated through the Caucasus and Central Asia, and arrived in North and Southwest China (lineages A, B, and C) and the Indian subcontinent (lineages B and C) through this region. Our results provide new insights into sheep domestication, particularly with respect to origins and migrations to and from eastern Eurasia.

Population subdivision of the surf clam Mactra chinensis in the East China Sea: Changjiang River outflow is not the sole driver

The northwestern Pacific, characterized by unique tectonic and hydrological settings, has greatly intrigued marine phylogeographers. However, current studies mostly focus on the influence of Pleistocene isolation of sea basins in population structure of species in the region, leaving the contribution of other factors (such as freshwater outflow and environmental gradients) largely unexploited. Here we shed light on the question by investigating phylogeography of the surf clam Mactra chinensis in the East China Sea (ECS). Genetic information was acquired from 501 specimens collected from its main distribution in the region, represented by mitochondrial cytochrome oxidase I (COI) and nine polymorphic microsatellite loci. A shallow and star-like phylogeny was revealed for all COI haplotypes, indicating the origin of populations from a single refugium. Although no divergent lineages existed, population subdivision was detected in both data sets. The most striking pattern was the significant differentiation between populations north and south of a biogeographic boundary—the Changjiang Estuary, suggesting a barrier effect of the freshwater outflow to gene flow. For the northern group, substructure was revealed by COI result as one southernmost population was significant different from other ones. Clear latitude gradations in allele frequencies were revealed by microsatellite analyses, likely influenced by environmental gradient factors such as temperature. Our results demonstrate that genetic subdivision can arise for populations within the ECS despite they have a single origin, and multiple mechanisms including Changjiang River outflow, environmental gradient factors and life-history traits may act in combination in the process.

Phylogeography of the sand dune ant Mycetophylax simplex along the Brazilian Atlantic Forest coast: remarkably low mtDNA diversity and shallow population structure

BACKGROUND

During past glacial periods, many species of forest-dwelling animals experienced range contractions. In contrast, species living outside such moist habitats appear to have reacted to Quaternary changes in different ways. The Atlantic Forest represents an excellent opportunity to test phylogeographic hypotheses, because it has a wide range of vegetation types, including unforested habitats covered predominantly by herbaceous and shrubby plants, which are strongly influenced by the harsh environment with strong wind and high insolation. Here, we investigated the distribution of genetic diversity in the endemic sand dune ant Mycetophylax simplex across its known range along the Brazilian coast, with the aim of contributing to the understanding of alternative phylogeographic patterns. We used partial sequences of the mitochondrial gene cytochrome oxidase I and nuclear gene wingless from 108 specimens and 51 specimens, respectively, to assess the phylogeography and demographic history of this species. To achieve this we performed different methods of phylogenetic and standard population genetic analyses.

RESULTS

The observed genetic diversity distribution and historical demographic profile suggests that the history of M. simplex does not match the scenario suggested for other Atlantic Forest species. Instead, it underwent demographic changes and range expansions during glacial periods. Our results show that M. simplex presents a shallow phylogeographic structure with isolation by distance among the studied populations, living in an almost panmictic population. Our coalescence approach indicates that the species maintained a stable population size until roughly 75,000 years ago, when it underwent a gradual demographic expansion that were coincident with the low sea-level during the Quaternary. Such demographic events were likely triggered by the expansion of the shorelines during the lowering of the sea level.

CONCLUSIONS

Our data suggest that over evolutionary time M. simplex did not undergo dramatic range fragmentation, but rather it likely persisted in largely interconnected populations. Furthermore, we add an important framework about how both glacial and interglacial events could positively affect the distribution and diversification of species. The growing number of contrasting phylogeographic patterns within and among species and regions have shown that Quaternary events influenced the distribution of species in more ways than first supposed.

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.

Phylogeography of amphi-boreal fish: tracing the history of the Pacific herring Clupea pallasii in North-East European seas

Background

The relationships between North Atlantic and North Pacific faunas through times have been controlled by the variation of hydrographic circumstances in the intervening Arctic Ocean and Bering Strait. We address the history of trans-Arctic connections in a clade of amphi-boreal pelagic fishes using genealogical information from mitochondrial DNA sequence data. The Pacific and Atlantic herrings (Clupea pallasii and C. harengus) have basically vicarious distributions in the two oceans since pre-Pleistocene times. However, remote populations of C. pallasii are also present in the border waters of the North-East Atlantic in Europe. These populations show considerable regional and life history differentiation and have been recognized in subspecies classification. The chronology of the inter-oceanic invasions and genetic basis of the phenotypic structuring however remain unclear.

Results

The Atlantic and Pacific herrings both feature high mtDNA diversities (large long-term population sizes) in their native basins, but an ocean-wide homogeneity of C. harengus is contrasted by deep east-west Pacific subdivision within Pacific C. pallasii. The outpost populations of C. pallasii in NE Europe are identified as members of the western Pacific C. pallasii clade, with some retained inter-oceanic haplotype sharing. They have lost diversity in colonization bottlenecks, but have also thereafter accumulated abundant new variation. The data delineate three phylogeographic groups within the European C. pallasii: herring from the inner White Sea; herring from the Mezen and Chesha Bays; and a strongly bottlenecked peripheral population in Balsfjord of the Norwegian Sea.

Conclusions

The NE European outposts of C. pallasii are judged to be early post-glacial colonists from the NW Pacific. A strong regional substructure has evolved since that time, in contrast to the apparent broad-scale uniformity maintained by herrings in their native basins. The structure only partly matches the previous biological concepts based on seasonal breeding stocks or geographical subspecies designations. The trans-Arctic herring phylogeography is notably similar to those of the amphi-boreal mollusk taxa Macoma and Mytilus, suggesting similar histories of inter-oceanic connections. We also considered the time dependency of molecular rates, critical for interpreting timing of relatively recent biogeographical events, by comparing the estimates from coding and non-coding mitochondrial regions of presumably different mutation dynamics.

What maintains the central North Pacific genetic discontinuity in Pacific herring?

Pacific herring show an abrupt genetic discontinuity in the central North Pacific that represents secondary contact between refuge populations previously isolated during Pleistocene glaciations. Paradoxically, high levels of gene flow produce genetic homogeneity among ocean-type populations within each group. Here, we surveyed variability in mtDNA control-region sequences (463 bp) and nine microsatellite loci in Pacific herring from sites across the North Pacific to further explore the nature of the genetic discontinuity around the Alaska Peninsula. Consistent with previous studies, little divergence (Φ(ST)  = 0.011) was detected between ocean-type populations of Pacific herring in the North West Pacific, except for a population in the Yellow Sea (Φ(ST)  = 0.065). A moderate reduction in genetic diversity for both mtDNA and microsatellites in the Yellow Sea likely reflects founder effects during the last colonization of this sea. Reciprocal monophyly between divergent mtDNA lineages (Φ(ST)  = 0.391) across the Alaska Peninsula defines the discontinuity across the North Pacific. However, microsatellites did not show a strong break, as eastern Bering Sea (EBS) herring were more closely related to NE Pacific than to NW Pacific herring. This discordance between mtDNA and microsatellites may be due to microsatellite allelic convergence or to sex-biased dispersal across the secondary contact zone. The sharp discontinuity between Pacific herring populations may be maintained by high-density blocking, competitive exclusion or hybrid inferiority.

High potential for using DNA from ancient herring bones to inform modern fisheries management and conservation

Pacific herring (Clupea pallasi) are an abundant and important component of the coastal ecosystems for the west coast of North America. Current Canadian federal herring management assumes five regional herring populations in British Columbia with a high degree of exchange between units, and few distinct local populations within them. Indigenous traditional knowledge and historic sources, however, suggest that locally adapted, distinct regional herring populations may have been more prevalent in the past. Within the last century, the combined effects of commercial fishing and other anthropogenic factors have resulted in severe declines of herring populations, with contemporary populations potentially reflecting only the remnants of a previously more abundant and genetically diverse metapopulation. Through the analysis of 85 archaeological herring bones, this study attempted to reconstruct the genetic diversity and population structure of ancient herring populations using three different marker systems (mitochondrial DNA (mtDNA), microsatellites and SNPs). A high success rate (91%) of DNA recovery was obtained from the extremely small herring bone samples (often <10 mg). The ancient herring mtDNA revealed high haplotype diversity comparable to modern populations, although population discrimination was not possible due to the limited power of the mtDNA marker. Ancient microsatellite diversity was also similar to modern samples, but the data quality was compromised by large allele drop-out and stuttering. In contrast, SNPs were found to have low error rates with no evidence for deviations from Hardy-Weinberg equilibrium, and simulations indicated high power to detect genetic differentiation if loci under selection are used. This study demonstrates that SNPs may be the most effective and feasible approach to survey genetic population structure in ancient remains, and further efforts should be made to screen for high differentiation markers.This study provides the much needed foundation for wider scale studies on temporal genetic variation in herring, with important implications for herring fisheries management, Aboriginal title rights and herring conservation.

Ancient divergence in the trans-oceanic deep-sea shark Centroscymnus crepidater

Unravelling the genetic structure and phylogeographic patterns of deep-sea sharks is particularly challenging given the inherent difficulty in obtaining samples. The deep-sea shark Centroscymnus crepidater is a medium-sized benthopelagic species that exhibits a circumglobal distribution occurring both in the Atlantic and Indo-Pacific Oceans. Contrary to the wealth of phylogeographic studies focused on coastal sharks, the genetic structure of bathyal species remains largely unexplored. We used a fragment of the mitochondrial DNA control region, and microsatellite data, to examine genetic structure in C. crepidater collected from the Atlantic Ocean, Tasman Sea, and southern Pacific Ocean (Chatham Rise). Two deeply divergent (3.1%) mtDNA clades were recovered, with one clade including both Atlantic and Pacific specimens, and the other composed of Atlantic samples with a single specimen from the Pacific (Chatham Rise). Bayesian analyses estimated this splitting in the Miocene at about 15 million years ago. The ancestral C. crepidater lineage was probably widely distributed in the Atlantic and Indo-Pacific Oceans. The oceanic cooling observed during the Miocene due to an Antarctic glaciation and the Tethys closure caused changes in environmental conditions that presumably restricted gene flow between basins. Fluctuations in food resources in the Southern Ocean might have promoted the dispersal of C. crepidater throughout the northern Atlantic where habitat conditions were more suitable during the Miocene. The significant genetic structure revealed by microsatellite data suggests the existence of present-day barriers to gene flow between the Atlantic and Pacific populations most likely due to the influence of the Agulhas Current retroflection on prey movements.