Late Glacial Demographic Expansion Motivates a Clock Overhaul for Population Genetics

Genetic diversity and demography of Bufo japonicus and B. torrenticola (Amphibia: Anura: Bufonidae) influenced by the Quaternary climate

The Quaternary climate affected the present species richness and geographic distribution patterns of amphibians by limiting their activities during the glacial period. The present study examined the phylogenetic relationships of Japanese toads (Bufo japonicus and B. torrenticola) and the demography of each lineage from the past to the present based on mitochondrial sequences and ecological niche models. Japanese toads are a monophyletic group with two main clades (clades A and B). Clade A represents B. j. formosus, including three clades (clades A1, A2, and A3). Clade B contains three clades, two of which corresponded to B. j. japonicus (clades B1 and B2) and the other to B. torrenticola. Clade B2 and B. torrenticola made a sister group, and, thus, B. j. japonicus is paraphyletic. Clades A and B diverged in the late Miocene 5.7 million years ago (Mya) during the period when the Japanese archipelago was constructed. The earliest divergence between the three clades of clade A was estimated at 1.8 Mya. Clades A1 and A2 may have diverged at 0.8 Mya, resulting from the isolation in the multiple different refugia; however, the effects of the glacial climate on the divergence events of clade A3 are unclear. Divergences within clade B occurred from the late Pliocene to the early Pleistocene (3.2-2.2 Mya). Niche similarity between the parapatric clade in clade B (clades B1 and B2) indicated their allopatric divergence. It was suggested that niche segregation between B. japonicus and B. torrenticola contributed to a rapid adaptation of B. torrenticola for lotic breeding. All clade of Japanese toads retreated to each refugium at a low elevation in the glacial period, and effective population sizes increased to construct the current populations after the Last Glacial Maximum. Furthermore, we highlight the areas of climate stability from the last glacial maximum to the present that have served as the refugia of Japanese toads and, thus, affected their present distribution patterns.

Coastal dunefields maintain pre-Holocene genetic structure in a rocky shore red alga

Most intertidal algae have limited dispersal potential, and areas that lack hard substratum suitable for attachment are thus expected to isolate regional populations from each other. Here, we used nuclear and mitochondrial genetic data to compare genetic structure in two co-distributed intertidal red algae with different dispersal potential along the South African coastline. Gelidium pristoides is divided into a south-eastern and a south-western evolutionary lineage separated by extensive, continuous sandy shoreline habitat adjacent to coastal dunefields. In contrast, Hypnea spicifera is genetically homogeneous throughout its range. In G. pristoides, the genetic breaks are associated with contemporary coastal dunefields. The age of the divergence event suggests that this may reflect the effect of older dispersal barriers, and that genetic structure was subsequently maintained by the formation of contemporary coastal dunefields.

Phylogeography and genetic diversity of the microbivalve Kidderia subquadrata, reveals new data from West Antarctic Peninsula

It is well established that Antarctic biodiversity has been strongly influenced by rapid climatic fluctuations during the Quaternary. Marine invertebrates from Antarctica constitute an interesting lens through which to study the impacts of the last glacial periods as glaciation impacted the distribution and intraspecific genetic variation of these animals. However, the impact on the spatial genetic distribution and historical demography of local processes in areas adjacent to the West Antarctic Peninsula (WAP) is less clear. Here we present new genetic information on the bivalve Kidderia subquadrata, a small mollusk that inhabits intertidal rocky island ecosystems throughout the WAP. Using a phylogeographical approach, we examined the spatial patterns of genetic diversity in this brooder species to test the hypothesis of strong genetic structure in incubating organisms and the hypothesis of glacial refugia in organisms with limited dispersion. We found evidence of strong genetic structure among populations of the WAP and a recent expansion in the South Shetland Islands. Our findings are concordant with the predictions that incubating organisms, abundant in Antarctica, present a strong genetic structure among their populations and also support the hypothesis of glacial refugia in organisms with limited dispersion. The effect of the coastal current pattern in the WAP is suggested as a driver to the local spatial dynamics of the genetic diversity distribution. Although genetic information about this microbivalve is still scarce, the knowledge reported here has increased our understanding of the evolutionary patterns of this organism that is endemic to the Southern Ocean.

Deep phylogeographic structure may indicate cryptic species within the Sparid genus Spondyliosoma

Two geographically nonoverlapping species are currently described within the sparid genus Spondyliosoma: Spondyliosoma cantharus (Black Seabream) occurring across Mediterranean and eastern Atlantic waters from NW Europe to Angola and S. emarginatum (Steentjie) considered endemic to southern Africa. To address prominent knowledge gaps this study investigated range-wide phylogeographic structure across both species. Mitochondrial DNA sequences revealed deep phylogeographic structuring with four regionally partitioned reciprocally monophyletic clades, a Mediterranean clade and three more closely related Atlantic clades [NE Atlantic, Angola and South Africa (corresponding to S. emarginatum)]. Divergence and distribution of the lineages reflects survival in, and expansion from, disjunct glacial refuge areas. Cytonuclear differentiation of S. emarginatum supports its validity as a distinct species endemic to South African waters. However, the results also indicate that S. cantharus may be a cryptic species complex wherein the various regional lineages represent established/incipient species. A robust multilocus genetic assessment combining morphological data and detailing interactions among lineages is needed to determine the full diversity within Spondyliosoma and the most adequate biological and taxonomic status.

Genetic analysis provides insights into species distribution and population structure in East Atlantic horse mackerel (Trachurus trachurus and T. capensis)

Two sister species of horse mackerel (Trachurus trachurus and T. capensis) are described that are intensively harvested in East Atlantic waters. To address long-standing uncertainties as to their respective geographical ranges, overlap and intraspecific population structure this study combined genetic (mitochondrial DNA and microsatellite) analysis and targeted sampling of the hitherto understudied West African coast. mtDNA revealed two reciprocally monophyletic clades corresponding to each species with interspecies nuclear differentiation supported by FST values. The T. trachurus clade was found across the north-east Atlantic down to Ghana but was absent from Angolan and South African samples. The T. capensis clade was found only in South Africa, Angola and a single Ghanaian individual. This pattern suggests that both species may overlap in the waters around Ghana. The potential for cryptic hybridization and/or indiscriminate harvesting of both species in the region is discussed. For T. capensis mtDNA supports high gene flow across the Benguela upwelling system, which fits with the species' ecology. The data add to evidence of a lack of significant genetic structure throughout the range of T. trachurus though the assumption of demographic panmixia is cautioned against. For both species, resolution of stock recruitment heterogeneity relevant to fishery management, as well as potential hybridization, will require more powerful genomic analyses.

Phylogeography of higher Diptera in glacial and postglacial grasslands in western North America

BACKGROUND

Pleistocene glaciations have had an important impact on the species distribution and community composition of the North American biota. Species survived these glacial cycles south of the ice sheets and/or in other refugia, such as Beringia. In this study, we assessed, using mitochondrial DNA from three Diptera species, whether flies currently found in Beringian grasslands (1) survived glaciation as disjunct populations in Beringia and in the southern refugium; (2) dispersed northward postglacially from the southern refugium; or (3) arose by a combination of the two. Samples were collected in grasslands in western Canada: Prairies in Alberta and Manitoba; the Peace River region (Alberta); and the southern Yukon Territory. We sequenced two gene regions (658 bp of cytochrome c oxidase subunit I, 510 bp of cytochrome b) from three species of higher Diptera: one with a continuous distribution across grassland regions, and two with disjunct populations between the regions. We used a Bayesian approach to determine population groupings without a priori assumptions and performed analysis of molecular variance (AMOVA) and exact tests of population differentiation (ETPD) to examine their validity. Molecular dating was used to establish divergence times.

RESULTS

Two geographically structured populations were found for all species: a southern Prairie and Peace River population, and a Yukon population. Although AMOVA did not show significant differentiation between populations, ETPD did. Divergence time between Yukon and southern populations predated the Holocene for two species; the species with an ambiguous divergence time had high haplotype diversity, which could suggest survival in a Beringian refugium.

CONCLUSIONS

Populations of Diptera in Yukon grasslands could have persisted in steppe habitats in Beringia through Pleistocene glaciations. Current populations in the region appear to be a mix of Beringian relict populations and, to a lesser extent, postglacial dispersal northward from southern prairie grasslands.

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.

Estimation of Evolutionary Rates of Mitochondrial DNA in Two Japanese Wood Mouse Species Based on Calibrations with Quaternary Environmental Changes

Reliable estimates of evolutionary rates of mitochondrial DNA might allow us to build realistic evolutionary scenarios covering broad time scales based on phylogenetic inferences. In the present study, we sought to obtain estimates of evolutionary rates in murine rodents using calibrations against historical biogeographic events. We first assumed that land-bridge-like structures that appeared intermittently at glacial maxima with 100,000-year intervals shaped the divergence patterns of cytochrome b (Cytb) sequences (1140 bp) of the larger Japanese wood mouse Apodemus speciosus. The comparison of sequences from peripheral remote islands that are separated from one another by deep straits allowed us to estimate mitochondrial DNA evolutionary rates (substitutions/site/million years) to be 0.027 to 0.036, with presumed calibrations from 140,000, 250,000, 350,000, and 440,000 years ago. Second, we addressed rapid expansion events inferred from analyses of the Cytb sequences of the lesser Japanese wood mouse A. argenteus. We detected five expansion signals in the dataset and established three categories based on the expansion parameter tau values: 3.9, 5.6-5.7, and 7.8-8.1. Considering that the climate became warmer 15,000, 53,000, and 115,000 years ago after preceding periods of rapid cooling, we calculated evolutionary rates to be 0.114, 0.047, and 0.031, respectively. This preliminary concept of the evolutionary rates on a time scale from 15,000 to 440,000 years ago for the wood mouse should be refined and tested in other species of murine rodents, including mice and rats.

Disentangling the genetic effects of refugial isolation and range expansion in a trans-continentally distributed species

In wide-ranging taxa with historically dynamic ranges, past allopatric isolation and range expansion can both influence the current structure of genetic diversity. Considering alternate historical scenarios involving expansion from either a single refugium or from multiple refugia can be useful in differentiating the effects of isolation and expansion. Here, we examined patterns of genetic variability in the trans-continentally distributed painted turtle (Chrysemys picta). We utilized an existing phylogeographic dataset for the mitochondrial control region and generated additional data from nine populations for the mitochondrial control region (n = 302) and for eleven nuclear microsatellite loci (n = 247). We created a present-day ecological niche model (ENM) for C. picta and hindcast this model to three reconstructions of historical climate to define three potential scenarios with one, two, or three refugia. Finally, we employed spatially-explicit coalescent simulations and an approximate Bayesian computation (ABC) framework to test which scenario best fit the observed genetic data. Simulations indicated that phylogeographic and multilocus population-level sampling both could differentiate among refugial scenarios, although inferences made using mitochondrial data were less accurate when a longer coalescence time was assumed. Furthermore, all empirical genetic datasets were most consistent with expansion from a single refugium based on ABC. Our results indicate a stronger role for post-glacial range expansion, rather than isolation in allopatric refugia followed by range expansion, in structuring diversity in this species. To distinguish among complex historical scenarios, we recommend explicitly modeling the effects of range expansion and evaluating alternate refugial scenarios for wide-ranging taxa.

Genetic structure and demographic inference of the regular sea urchin Sterechinus neumayeri (Meissner, 1900) in the Southern Ocean: The role of the last glaciation

One of the most relevant characteristics of the extant Southern Ocean fauna is its resiliency to survive glacial processes of the Quaternary. These climatic events produced catastrophic habitat reductions and forced some marine benthic species to move, adapt or go extinct. The marine benthic species inhabiting the Antarctic upper continental shelf faced the Quaternary glaciations with different strategies that drastically modified population sizes and thus affected the amount and distribution of intraspecific genetic variation. Here we present new genetic information for the most conspicuous regular sea urchin of the Antarctic continental shelf, Sterechinus neumayeri. We studied the patterns of genetic diversity and structure in this broadcast-spawner across three Antarctic regions: Antarctic Peninsula, the Weddell Sea and Adélie Land in East Antarctica. Genetic analyses based on mitochondrial and nuclear markers suggested that S. neumayeri is a single genetic unit around the Antarctic continent. The species is characterized by low levels of genetic diversity and exhibits a typical star-like haplotype genealogy that supports the hypothesis of a single in situ refugium. Based on two mutation rates standardized for this genus, the Bayesian Skyline plot analyses detected a rapid demographic expansion after the Last Glacial Maximum. We propose a scenario of rapid postglacial expansion and recolonization of Antarctic shallow areas from a less ice-impacted refugium where the species survived the LGM. Considering the patterns of genetic diversity and structure recorded in the species, this refugium was probably located in East Antarctica.

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.

Significant genetic differentiation of Gobiopterus lacustris, a newly recorded transparent goby in China

In our recent survey, the transparent small Lacustrine goby, Gobiopterus lacustris had reported as the endemic species of Luzon, Philippines, was identified as an abundant species in mangroves of Leizhou Peninsula, China. Here, high diversity and significant differentiation of five sites of samples representing the west and east populations were revealed by mitochondrial DNA sequences. Five haplotypes of 56 cytochrome oxidase subunit I (Cox1) with the lengths of 623 base pairs (bp) have the high pairwise identity (>98.8%). Moreover, a total of 31 haplotypes for 129 partial D-loop regions were clustered into two clades corresponding to the east and west sampling sites. The strong population structure was confirmed (ΦST = 0.43017, p < .0001) with high haplotype diversity (h = 0.880 ± 0.017) and low nucleotide diversity (p=.00484). Moreover, both the mismatch distribution analysis and neutral test of D-loop revealed that the west group might experience a recent demographic expansion. Lastly, the isolation-with-migration analysis supported the expansion and indicated that the east-west split happened at approximately 7.1 kyr ago. Given the distribution and diversity, G. lacustris could be a good model for the study of the sea-level fluctuations and coast evolution of the South China Sea.

Landscape genetics indicate recently increased habitat fragmentation in African forest-associated chafers

Today, indigenous forests cover less than 0.6% of South Africa's land surface and are highly fragmented. Most forest relicts are very small and typically occur in fire-protected gorges along the eastern Great Escarpment. Yet, they hold a unique and valuable fauna with high endemism and ancient phylogenetic lineages, fostered by long-term climatic stability and complex microclimates. Despite numerous studies on southern African vegetation cover, the current state of knowledge about the natural extension of indigenous forests is rather fragmentary. We use an integrated approach of population-level phylogeography and climatic niche modeling of forest-associated chafer species to assess connectivity and extent of forest habitats since the last glacial maximum. Current and past species distribution models ascertained potential fluctuations of forest distribution and supported a much wider potential current extension of forests based on climatic data. Considerable genetic admixture of mitochondrial and nuclear DNA among many populations and an increase in mean population mutation rate in Extended Bayesian Skyline Plots of all species indicated more extended or better connected forests in the recent past (<5 kya). Genetic isolation of certain populations, as revealed by population differentiation statistics (GST'), as well as landscape connectivity statistics and habitat succession scenarios suggests considerable loss of habitat connectivity. As major anthropogenic influence is likely, conservational actions need to be considered.

Secondary contact and asymmetrical gene flow in a cosmopolitan marine fish across the Benguela upwelling zone

The combination of oceanographic barriers and habitat heterogeneity are known to reduce connectivity and leave specific genetic signatures in the demographic history of marine species. However, barriers to gene flow in the marine environment are almost never impermeable which inevitably allows secondary contact to occur. In this study, eight sampling sites (five along the South African coastline, one each in Angola, Senegal and Portugal) were chosen to examine the population genetic structure and phylogeographic history of the cosmopolitan bluefish (Pomatomus saltatrix), distributed across a large South-east Atlantic upwelling zone. Molecular analyses were applied to mtDNA cytochrome b, intron AM2B1 and 15 microsatellite loci. We detected uncharacteristically high genetic differentiation (F_ST 0.15-0.20; P<0.001) between the fish sampled from South Africa and the other sites, strongly influenced by five outlier microsatellite loci located in conserved intergenic regions. In addition, differentiation among the remaining East Atlantic sites was detected, although mtDNA indicated past isolation with subsequent secondary contact between these East Atlantic populations. We further identified secondary contact, with unidirectional gene flow from South Africa to Angola. The directional contact is likely explained by a combination of the northward flowing offshore current and endogenous incompatibilities restricting integration of certain regions of the genome and limiting gene flow to the south. The results confirm that the dynamic system associated with the Benguela current upwelling zone influences species distributions and population processes in the South-east Atlantic.