Phylogeography of the flathead mullet Mugil cephalus in the north-west Pacific as inferred from the mtDNA control region

The population genetic structure and historical demography of the flathead mullet Mugil cephalus were investigated using the mtDNA control region (CR) sequences (909–1015 bp) of 126 individuals collected from seven locations in the north-west Pacific between 2005 and 2007. Haplotype diversity ( h = 0·9333–1·000) and nucleotide diversity (π = 0·0046–0·1467) varied greatly among the sampling locations. Phylogenetic analysis of the CR sequences indicated that M. cephalus in the north-west Pacific belongs to two highly divergent lineages (lineages 1 and 2), with the inferred population structure being closely associated with the distribution of both lineages. Two populations were identified, one from the East China Sea and the other from the South China Sea. The former samples were obtained from Taiwan and Qingdao of north China and associated with lineage 1 haplotypes. The latter samples were collected from the Philippines, Pearl River of South China and two samples from Japan, all of which were associated with lineage 2. Japanese samples from Okinawa and Yokosuka had different degrees of mixing between lineages 1 and 2. Historical demographic variables in both populations indicated that Pleistocene glaciations had a strong impact on M. cephalus in the north-west Pacific, resulting in a recent demographic decline of the East China Sea population but in demographic equilibrium for the South China Sea population. Japan appears to be a contact zone between lineages 1 and 2, but it may also be indicative of coexistence between resident and migratory populations. Further global studies are required to clarify the taxonomic status of this cosmopolitan species.     

East Asian allopatry and north Eurasian sympatry in Long‐tailed Tit lineages despite similar population dynamics during the late Pleistocene

Eurasia is a large continent characterized by heterogeneous environments. Glacial cycles during the late Pleistocene have had variable impacts on the avifauna across Eurasia. Bird populations from South‐East Asia show stability through the Last Glacial Maximum (LGM), while populations from Europe exhibit evidence of post‐LGM expansion. We investigated the phylogeography of the Long‐tailed Tit (Aegithalos caudatus), which spans the longitudinal breadth of Eurasia to test how climatic history and regional topographical complexity affected populations and diversification within the species complex. Our results show that two lineages from central and southern China (lineages C and D) segregate geographically, while lineages across northern Eurasia (lineage A and B) show substantial sympatry. Bayesian estimates for the timing of diversification suggest that the four lineages diverged during the middle Pleistocene, splitting in parallel and undergoing concurrent demographic histories since divergence. A. caudatus lineages experienced similar and synchronous population size dynamics during glacial cycles before the LGM. We conclude that the difference in geo‐topologic complexity may be an important factor that led to the variation in secondary admixture between northern Eurasian and eastern Asian lineages.     

Phylogeography and historical demography of the anadromous fish Leucopsarion petersii in relation to geological history and oceanography around the Japanese Archipelago

Phylogeographical patterns of marine and diadromous organisms are often influenced by dynamic ocean histories. For example, the marine realm around the Japanese Archipelago is an interesting area for phylogeographical research because of the wide variation in the environments driven by repeated shifts in sea level in the Quaternary. We analysed mitochondrial cyt b gene and nuclear myh6 gene sequences for individuals collected from throughout the range of the anadromous fish Leucopsarion petersii to assess the lineage divergence, phylogeographical pattern and historical demography in relation to geological history and oceanographic features around the archipelago. Leucopsarion petersii has two major lineages (the Japan Sea and Pacific Ocean lineages), which diverged during the late-early to middle Pleistocene. Geographical distributions of the two lineages were closely related to the pathways of the two warm currents, the Tsushima Current and the Kuroshio Current, that flow past the archipelago. Evidence of introgressive hybridization between these lineages was found at two secondary contact zones. Demographic tests suggested that the Japan Sea and Pacific Ocean lineages carried the genetic signal of different historical demographic processes, and these signals are probably associated with differences in habitat stability during recent glacial periods. The Japan Sea lineage has a larger body-size and more vertebrae, probably in relation to severe habitat conditions through Pleistocene climatic oscillations. Thus, the two lineages have long independent evolutionary histories, and the phylogeographical structure and demography of this species have been influenced both by historical events and the present-day oceanography around the Japanese Archipelago.     

Phylogeography study of Ammodytes personatus in Northwestern Pacific: Pleistocene isolation, temperature and current conducted secondary contact

To assess the role of historical process and contemporary factors in shaping population structures in Northwestern Pacific, mitochondrial control region sequences were analyzed to characterize the phylogeography and population structure of the Japanese sand lance Ammodytes personatus. A total of 429 individuals sampled from 17 populations through the species' range are sequenced. Two distinct lineages are detected, which might have been divergent in the Sea of Japan and Pacific costal waters of Japanese Island, during the low sea level. Significant genetic structure is revealed between the Kuroshio and Oyashio Currents. However, significant genetic structure is also detected in the Sea of Japan, contracting expected homogenization hypothesis in Tsushima Current. The haplotype frequency of lineages in both sides of Japanese Island and significant genetic structure between north and south groups revealed that the distribution of lineage B and north group were highly limited by the annual sea temperature. The lack of lineage B in Qingdao population with low sea temperature reflects the sea temperature barrier. Lack of genetic structure in the south group and north group populations indicated that ocean currents within groups facilitated the dispersal of A. personatus.     

Profound climatic effects on two East Asian black-throated tits (Ave: Aegithalidae), revealed by ecological niche models and phylogeographic analysis

Although a number of studies have assessed the effects of geological and climatic changes on species distributions in East Asian, we still have limited knowledge of how these changes have impacted avian species in south-western and southern China. Here, we aim to study paleo-climatic effects on an East Asian bird, two subspecies of black-throated tit (A. c. talifuensis-concinnus) with the combined analysis of phylogeography and Ecological Niche Models (ENMs). We sequenced three mitochondrial DNA markers from 32 populations (203 individuals) and used phylogenetic inferences to reconstruct the intra-specific relationships among haplotypes. Population genetic analyses were undertaken to gain insight into the demographic history of these populations. We used ENMs to predict the distribution of target species during three periods; last inter-glacial (LIG), last glacial maximum (LGM) and present. We found three highly supported, monophyletic MtDNA lineages and different historical demography among lineages in A. c. talifuensis-concinnus. These lineages formed a narrowly circumscribed intra-specific contact zone. The estimated times of lineage divergences were about 2.4 Ma and 0.32 Ma respectively. ENMs predictions were similar between present and LGM but substantially reduced during LIG. ENMs reconstructions and molecular dating suggest that Pleistocene climate changes had triggered and shaped the genetic structure of black-throated tit. Interestingly, in contrast to profound impacts of other glacial cycles, ENMs and phylogeographic analysis suggest that LGM had limited effect on these two subspecies. ENMs also suggest that Pleistocene climatic oscillations enabled the formation of the contact zone and thus support the refuge theory.     

Lineage admixture during postglacial range expansion is responsible for the increased gene diversity of Kalopanax septemlobus in a recently colonised territory

We aimed to reveal the effects of range expansion and subsequent lineage admixture from separated glacial refugia on genetic diversity of Kalopanax septemlobus in Japan, by combining nuclear microsatellite data and ecological niche modelling. Allelic richness and gene diversity were compared at the population and regional level. We also statistically examined these indices as a function of population accessibility to the last glacial maximum (LGM) palaeodistribution reconstructed by ecological niche modelling to test a simple range expansion scenario from glacial refugia. Genetic diversity was highest in the populations of southern Japan and gradually decreased towards the north. However, an additional centre of genetic diversity, when measured as gene diversity, was found in northern Honshu Island, where distinct lineages were shown to be in contact. Positive effects of population accessibility to the LGM range were detected in both diversity indices at different spatial scales. The combined data support independent postglacial range expansions towards the north from the edge populations on the exposed coastal shelf of Pacific and Sea of Japan in northern Honshu during the LGM, which subsequently resulted in markedly low genetic diversity in the northernmost extant range, Hokkaido. The regional increase in gene diversity in northern Honshu is likely to be the result of postglacial lineage admixture. Relative difference in the spatial scales best relating population genetic diversity with the LGM distribution can be explained by a higher rate of allelic richness diversity loss during range expansions and stronger effects of lineage admixture on gene diversity.     

Diversification and Demography of the Oriental Garden Lizard (Calotes versicolor) on Hainan Island and the Adjacent Mainland

The Oriental garden lizard (Calotes versicolor) is one of the few non-gekkonid lizards that are geographically widespread in the tropics. We investigated its population dynamics on Hainan Island and the adjacent mainland of China and Vietnam, focusing on the impact of cyclic upheaval and submergence of land bridges during the Pleistocene. Our Bayesian phylogenetic analysis reveals two mitochondrial lineages, A and B, which are estimated to have coalesced about 0.26 million years ago (95% credibility interval: 0.05–0.61 million years ago). Lineage A contains individuals mainly from central and southern Wuzhi Mountain on Hainan Island, whereas lineage B mainly comprises individuals from other sites on the island plus the adjacent mainland. The estimated coalescence times within lineages A (0.05 million years ago) and B (0.13 million years ago) fall within a period of cyclical land-bridge formation and disappearance in the Pleistocene. A spatial analysis of molecular variance identified two distinct population groupings: I, primarily containing lineage A, and II, mainly consisting of lineage B. However, haplotypes from lineages A and B occur sympatrically, suggesting that gene flow is ongoing. Neither Wuzhi Mountain nor Qiongzhou Strait and Gulf of Tonkin act as barriers to gene flow among C. versicolor populations. Analyses of the data using mismatch distributions and extended Bayesian skyline plots provide evidence of a relatively stable population size through time for Group I, and moderate population expansions and contractions during the end of the Pleistocene for Group II. We conclude that the phylogeographical patterns of C. versicolor are the combined product of Pleistocene sea-level oscillations and nonphysical barriers to gene flow.     

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri,Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.     

Phylogeography of the Northern Alligator Lizard (Squamata,Anguidae): Hidden diversity in a western endemic

Western North America includes the California Floristic Province and the Pacific Northwest, biologically diverse regions highlighted by a complex topography, geology, climate and history. A number of animals span these regions and show distinctive patterns of dispersal, vicariance and lineage diversification. Examining phylogeographic patterns in the fauna of this area aids in our understanding of the forces that have contributed to the generation and maintenance of regional biodiversity. Here, we investigate the biogeography and population structure of the Northern Alligator Lizard (Elgaria coerulea), a wide‐ranging anguid endemic to western North America. We sequenced two mtDNA fragments (ND2 and ND4) for 181 individuals across the range of the species and analysed these data with phylogenetic approaches to infer population and biogeographic history, and date major divergences within the taxon. We further used Bayesian clustering methods to assess major patterns of population structure and performed ecological niche modelling (ENM) to aid in our interpretation of geographic structure and diversification of E. coerulea lineages. Our phylogeographic examination of E. coerulea uncovered surprising diversity and structure, recovering 10 major lineages, each with substantial geographic substructure. While some divergences within the species are relatively old (Pliocene, 5.3–2.6 mya), most intraspecific variation appears to be of more recent origin (Pleistocene, 2.6 mya‐11,700 ya). Current diversity appears to have arisen in the Sierra Nevada Mountains and spread west and north since the Pliocene. Finally, our ENMs suggest that much of the Coast Ranges in California provided ideal habitat during the Last Glacial Maxima (LGM) that has since contracted dramatically and shifted northwards, whereas significant portions of the Sierra Nevada were unsuitable during the LGM and have since become more suitable. Interestingly, E. coerulea shares a number of genetic boundaries with other sympatric taxa, suggesting common historical events and geomorphological features have shaped the biota of this region.     

Phylogeographic Patterns of mtDNA Variation Revealed Multiple Glacial Refugia for the Frog Species Feirana taihangnica Endemic to the Qinling Mountains

Diversification patterns and demography of montane species are affected by Pleistocene climate fluctuations. Empirical cases from the Qinling Mountains (QM) region, which is a major biogeographic divider of East Asia, are few. We used DNA sequence data of the complete mitochondrial ND2 gene to detect effects of the Pleistocene glaciations on phylogeographic profiles of a frog species, Feirana taihangnica, which is endemic to the QM. Four distinct lineages consisting of seven sublineages were revealed. The strongest signal of biogeographical structure (F _ct = 0.971, P < 0.01) was found when populations were grouped according to these seven sublineages. One narrow secondary contact zone was detected in the middle QM between the lineage from middle QM and the lineage from eastern QM. Coalescent simulations indicated that this species colonized the QM region by a stepping-stone model. Divergences among lineages had likely been influenced by the uplift of the Tibetan Plateau during the late Miocene-to-late Pleistocene, as well as by the Pleistocene climatic cycles. Coalescent simulations also suggested that F. taihangnica populations have persisted through the Pleistocene glacial periods in multiple refugia across the QM region. Demographic analyses indicated that all lineages, except the lineage in the Funiu Mountains, have been experienced postglacial expansion of population size and distribution range. In conclusion, Pleistocene climate fluctuations and tectonic changes during the late Miocene-late Pleistocene have profoundly influenced the phylogeography and historical demography of F. taihangnica.     

Diversity and trans-arctic invasion history of mitochondrial lineages in the North Atlantic Macoma balthica complex (Bivalvia: Tellinidae)

The history of repeated inter- or transoceanic invasions in bivalve mollusks of the circumpolar Macoma balthica complex was assessed from mtDNA COIII sequences. The data suggest that four independent trans-Arctic invasions, from the Pacific, gave rise to the current lineage diversity in the North Atlantic. Unlike in many other prominent North Atlantic littoral taxa, no evidence for (postinvasion) trans-Atlantic connections was found in the M. balthica complex. The earliest branch of the mtDNA tree is represented by the temperate-boreal North American populations (=Macoma petalum), separated from the M. balthica complex proper in the Early Pliocene at latest. The ensuing trans-Arctic invasions established the North European M. b. rubra, which now prevails on the North Sea and northeast Atlantic coasts, about two million years ago, and the currently northwest Atlantic M. balthica lineage in the Canadian Maritimes, in the Middle Pleistocene. The final reinvasion(s) introduced a lineage that now prevails in a number of North European marginal seas and is still hardly distinguishable from North Pacific mtDNA (M. b. balthica). We used coalescence simulation analyses to assess the age of the latest invasion from the Pacific to the northeast Atlantic. The results refute the hypothesis of recent, human-mediated reintroductions between northeast Pacific and the North European marginal seas in historical times. Yet they also poorly fit the alternative hypotheses of an early postglacial trans-Arctic invasion (< 11 thousand years ago), or an invasion during the previous Eemian interglacial (120 thousand years ago). Divergence time estimates rather fall in the Middle Weichselian before the Last Glacial Maximum, in conflict with the conventional thinking of trans-Arctic biogeographical connections; an early Holocene reinvasion may still be regarded as the most plausible scenario. Today, the most recently invaded Pacific mtDNA lineage is found admixed with the earlier established European Atlantic "rubra" lineage in the Baltic Sea and in Barents Sea populations east of the Varanger peninsula, and it is practically exclusive in the White and Pechora seas. Yet mtDNA does not always constitute an unequivocal taxonomic marker at individual level; the marginal populations represent hybrid swarms of the Atlantic and Pacific lineages in their nuclear genes.     

Molecular data and ecological niche modelling reveal the Pleistocene history of a semi‐aquatic bug (Microvelia douglasi douglasi) in East Asia

This study investigated the Pleistocene history of a semi‐aquatic bug, Microvelia douglasi douglasi Scott, 1874 (Hemiptera: Veliidae) in East Asia. We used M. douglasi douglasi as a model species to explore the effects of historical climatic fluctuations on montane semi‐aquatic invertebrate species. Two hypotheses were developed using ecological niche models (ENMs). First, we hypothesized that M. douglasi douglasi persisted in suitable habitats in southern Guizhou, southern Yunnan, Hainan, Taiwan and southeast China during the LIG. After that, the populations expanded (Hypothesis 1). As the spatial prediction in the LGM was significantly larger than in the LIG, we then hypothesized that the population expanded during the LIG to LGM transition (Hypothesis 2). We tested these hypotheses using mitochondrial data (COI+COII) and nuclear data (ITS1 + 5.8S+ITS2). Young lineages, relatively deep splits, lineage differentiation among mountain ranges in central, south and southwest China and high genetic diversities were observed in these suitable habitats. Evidence of mismatch distributions and neutrality tests indicate that a population expansion occurred in the late Pleistocene. The Bayesian skyline plot (BSP) revealed an unusual population expansion that likely happened during the cooling transition between LIG and LGM. The results of genetic data were mostly consistent with the spatial predictions from ENM, a finding that can profoundly improve phylogeographic research. The ecological requirements of M. douglasi douglasi, together with the geographical heterogeneity and climatic fluctuations of Pleistocene in East Asia, could have shaped this unusual demographic history. Our study contributes to our knowledge of semi‐aquatic bug/invertebrate responses to Pleistocene climatic fluctuations in East Asia.     

Phylogeographic structure and deep lineage diversification of the red alga Chondrus ocellatus Holmes in the Northwest Pacific

A major goal of phylogeographic analysis using molecular markers is to understand the ecological and historical variables that influence genetic diversity within a species. Here, we used sequences of the mitochondrial Cox1 gene and nuclear internal transcribed spacer to reconstruct its phylogeography and demographic history of the intertidal red seaweed Chondrus ocellatus over most of its geographical range in the Northwest Pacific. We found three deeply separated lineages A, B and C, which diverged from one another in the early Pliocene–late Miocene (c. 4.5–7.7 Ma). The remarkably deep divergences, both within and between lineages, appear to have resulted from ancient isolations, accelerated by random drift and limited genetic exchange between regions. The disjunct distributions of lineages A and C along the coasts of Japan may reflect divergence during isolation in scattered refugia. The distribution of lineage B, from the South China Sea to the Korean Peninsula, appears to reflect postglacial recolonizations of coastal habitats. These three lineages do not coincide with the three documented morphological formae in C. ocellatus, suggesting that additional cryptic species may exist in this taxon. Our study illustrates the interaction of environmental variability and demographic processes in producing lineage diversification in an intertidal seaweed and highlights the importance of phylogeographic approaches for discovering cryptic marine biodiversity.     

In situ genetic differentiation in a Hispaniolan lizard (Ameiva chrysolaema): a multilocus perspective

A previous phylogeographic study of mitochondrial haplotypes for the Hispaniolan lizard Ameiva chrysolaema revealed deep genetic structure associated with seawater inundation during the late Pliocene/early Pleistocene and evidence of subsequent population expansion into formerly inundated areas. We revisit hypotheses generated by our previous study using increased geographic sampling of populations and analysis of three nuclear markers (alpha-enolase intron 8, alpha-cardiac-actin intron 4, and beta-actin intron 3) in addition to mitochondrial haplotypes (ND2). Large genetic discontinuities correspond spatially and temporally with historical barriers to gene flow (sea inundations). NCPA cross-validation analysis and Bayesian multilocus analyses of divergence times (IMa and MCMCcoal) reveal two separate episodes of fragmentation associated with Pliocene and Pleistocene sea inundations, separating the species into historically separate Northern, East-Central, West-Central, and Southern population lineages. Multilocus Bayesian analysis using IMa indicates asymmetrical migration from the East-Central to the West-Central populations following secondary contact, consistent with expectations from the more pervasive sea inundation in the western region. The West-Central lineage has a genetic signature of population growth consistent with the expectation of geographic expansion into formerly inundated areas. Within each lineage, significant spatial genetic structure indicates isolation by distance at comparable temporal scales. This study adds to the growing body of evidence that vicariant speciation may be the prevailing source of lineage accumulation on oceanic islands. Thus, prior theories of island biogeography generally underestimate the role and temporal scale of intra-island vicariant processes.     

Parallel responses of bees to Pleistocene climate change in three isolated archipelagos of the southwestern Pacific

The impacts of glacial cycles on the geographical distribution and size of populations have been explored for numerous terrestrial and marine taxa. However, most studies have focused on high latitudes, with only a few focused on the response of biota to the last glacial maximum (LGM) in equatorial regions. Here, we examine how population sizes of key bee fauna in the southwest Pacific archipelagos of Fiji, Vanuatu and Samoa have fluctuated over the Quaternary. We show that all three island faunas suffered massive population declines, roughly corresponding in time to the LGM, followed by rapid expansion post-LGM. Our data therefore suggest that Pleistocene climate change has had major impacts across a very broad tropical region. While other studies indicate widespread Holarctic effects of the LGM, our data suggest a much wider range of latitudes, extending to the tropics, where these climate change repercussions were important. As key pollinators, the inferred changes in these bee faunas may have been critical in the development of the diverse Pacific island flora. The magnitude of these responses indicates future climate change scenarios may have alarming consequences for Pacific island systems involving pollinator-dependent plant communities and agricultural crops.     

Distributional and demographic consequences of Pleistocene climate fluctuations for a marine demersal fish in the north-eastern Atlantic

Aim The Pleistocene glaciations were the most significant historical event during the evolutionary life span of most extant species. However, little is known about the consequences of these climate changes for the distribution and demography of marine animals of the north‐eastern Atlantic. The present study focuses on the phylogeographic and demographic patterns of the sand goby, Pomatoschistus minutus (Teleostei: Gobiidae), a small marine demersal fish. Location North‐eastern Atlantic, Mediterranean, Irish, North and Baltic seas. Methods Analysis was carried out by sequencing the mtDNA cytochrome b gene of sand gobies from 12 localities throughout the species’ range, and using this information in combination with published data of allozyme markers and mtDNA control region sequences. Several phylogenetic methods and a network analysis were used to explore the phylogeographic pattern. The historical demography of P. minutus was studied through a mismatch analysis and a Bayesian skyline plot. Results Reciprocal monophyly was found between a Mediterranean Sea (MS) clade and an Atlantic Ocean (AO) clade, both with a Middle Pleistocene origin. The AO Clade contains two evolutionary significant units (ESUs): the Iberian Peninsula (IB) Group and the North Atlantic (NA) Group. These two groups diverged during Middle Pleistocene glacial cycles. For the NA Group there is evidence for geographic sorting of the ancestral haplotypes with recent radiations in the Baltic Sea, Irish Sea, North Sea and Bay of Biscay. The demographic histories of the Mediterranean Clade and the two Atlantic ESUs were influenced mainly by expansions dated as occurring during the Middle Pleistocene glaciations and post‐Eem, respectively. Main conclusions The pre‐LGM (Last Glacial Maximum) subdivision signals were not erased for P. minutus during the LGM. Middle Pleistocene glaciations yielded isolated and differently evolving sets of populations. In contrast to the case for most other taxa, only the northern Atlantic group contributed to the post‐glacial recolonization. The historical demography of Mediterranean sand gobies was influenced mainly by Middle Pleistocene glaciations, in contrast to that of the Atlantic populations, which was shaped by Late Pleistocene expansions.     

Phylogeography of sugar kelp: Northern ice‐age refugia in the Gulf of Alaska

Many Northeast (NE) Pacific fishes and invertebrates survived Pleistocene glaciations in northern refugia, but the extent that kelps survived in northern areas is uncertain. Here, we test the hypothesis that populations of sugar kelp (Saccharina latissima) persisted in the Gulf of Alaska during ice‐age maxima when the western margin of the Cordilleran ice sheet covered coastal areas around the NE Pacific Ocean. We estimated genetic diversities within and phylogeographical relationships among 14 populations along 2,800 km in the NE Pacific and Bering Sea with partial sequences of mitochondrial DNA 5′‐cytochrome oxidase subunit I (COI, bp = 624, n = 543), chloroplast DNA ribulose‐1,5‐bisphosphate carboxylase large subunit‐3′ (rbcL, bp = 735, n = 514), and 11 microsatellite loci. Concatenated sequences of rbcL and COI showed moderate levels of within‐population genetic diversity (mean h = 0.200) but substantial differences among populations (Φ_ST = 0.834, p < .0001). Microsatellites showed moderate levels of heterozygosity within populations (mean H _E = 0.391). Kelps in the same organellar lineage tended to cluster together, regardless of geographic origins, as indicated in a principal coordinate analysis (PCoA) of microsatellite genotypes. The PCoA also showed evidence of nuclear hybridizations between co‐occurring organellar lineages. Individual admixture plots with population clusters of K = 2, 6, and 9 showed increasing complexity with considerable historical admixture between some clusters. A time‐calibrated phylogeny placed divergences between rbcL‐COI lineages at 1.4 million years at most. The time frames of mutation in the rbcL‐COI lineages and microsatellite population clusters differed among locations. The existence of ancient lineages in the Gulf of Alaska, moderate levels of genetic diversity, and the absence of departures from neutrality are consistent with northern refugia during multiple Croll‐Milankovitch climate cycles in the Pleistocene Epoch.     

Comparative Analysis of Japanese Three-Spined Stickleback Clades Reveals the Pacific Ocean Lineage Has Adapted to Freshwater Environments while the Japan Sea Has Not

Divergent selection and adaptive divergence can increase phenotypic diversification amongst populations and lineages. Yet adaptive divergence between different environments, habitats or niches does not occur in all lineages. For example, the colonization of freshwater environments by ancestral marine species has triggered adaptive radiation and phenotypic diversification in some taxa but not in others. Studying closely related lineages differing in their ability to diversify is an excellent means of understanding the factors promoting and constraining adaptive evolution. A well-known example of the evolution of increased phenotypic diversification following freshwater colonization is the three-spined stickleback. Two closely related stickleback lineages, the Pacific Ocean and the Japan Sea occur in Japan. However, Japanese freshwater stickleback populations are derived from the Pacific Ocean lineage only, suggesting the Japan Sea lineage is unable to colonize freshwater. Using stable isotope data and trophic morphology, we first show higher rates of phenotypic and ecological diversification between marine and freshwater populations within the Pacific Ocean lineage, confirming adaptive divergence has occurred between the two lineages and within the Pacific Ocean lineage but not in the Japan Sea lineage. We further identified consistent divergence in diet and foraging behaviour between marine forms from each lineage, confirming Pacific Ocean marine sticklebacks, from which all Japanese freshwater populations are derived, are better adapted to freshwater environments than Japan Sea sticklebacks. We suggest adaptive divergence between ancestral marine populations may have played a role in constraining phenotypic diversification and adaptive evolution in Japanese sticklebacks.     

Mitochondrial DNA phylogeography of the three-spined stickleback (Gasterosteus aculeatus) in Europe-evidence for multiple glacial refugia

An analysis of mitochondrial DNA sequence variation in 172 three-spined sticklebacks (Gasterosteus aculeatus) sampled across the European distribution range revealed three major evolutionary lineages occupying relatively large and separate geographic areas. The trans-Atlantic lineage comprised of populations spanning from the East Coast of USA to the continental Europe and was basal group to the other European lineages in the phylogeny. The European lineage included populations located in the Western and Eastern Europe, British Isles, Scandinavia as well as some parts of the Mediterranean region. The third lineage was specific to the Black Sea drainages. The within lineage structure was characterized by significant excess of low frequency haplotypes and star-like mtDNA genealogies, which suggest a recent population expansions to the formerly glaciated marine and freshwater environments. A coalescent-based method dated the splits between the major lineages to have occurred during the Saalian and Weichselian glaciations in the late Pleistocene, depending on the molecular clock calibration. The coalescent simulations further indicate high degree of genetic diversity within the lineages and a substantial increase in the genetic diversity in the European lineage relative to the ancestral level. In addition to the three major lineages, the freshwater populations in R. Neretva and L. Skadar in the Adriatic Sea coast region harboured unique and highly divergent haplotypes suggesting long independent histories of these populations. Evidence from mtDNA analyses suggests that these populations deserve a status of an evolutionary significant unit.     

DIFFERING LEVELS OF AMONG-POPULATION DIVERGENCE IN THE MITOCHONDRIAL DNA OF PERIODICAL CICADAS RELATED TO HISTORICAL BIOGEOGRAPHY

Mitochondrial DNA (mtDNA) haplotypes were determined for 118 individuals of 13-and 17-year periodical cicadas (genus Magicicada) collected from 16 localities throughout the Midwest and eastern United States. Two distinct mtDNA lineages, identified as A and B, differ by 2.5% based on analysis of fragment patterns and restriction maps. Observed levels of mtDNA diversity within each lineage are low compared to estimates for other taxa. The two lineages are regionally segregated, with the boundary line occurring at a latitude of approximately 33° North. The levels of mtDNA diversity and population genetic structure differ within the two lineages. There is a remarkably low level of mean mtDNA divergence and no genetic structure in lineage A, whereas lineage B exhibits an order of magnitude higher level of mtDNA diversity and significant genetic structure among sampled populations. The low level of mtDNA diversity in cicadas may be attributed to (1) a population bottleneck that most likely occurred during the Pleistocene, (2) recent colonization following the retreat of the glaciers and the expansion of deciduous forests, and/or (3) high among-family reproductive variance (as a consequence of large population size, high fecundity, aggregative behavior of adults, and clumping of eggs). The difference in mtDNA diversity and population genetic structure between the lineages suggests that they experienced different biogeographic histories; we relate this to Pleistocene changes.