Chloroplast DNA phylogeography of Fagus crenata (Fagaceae) in Japan

 CpDNA variation in Japanese beech, Fagus crenata (Fagaceae), was studied in 45 populations distributed throughout the species' range. Two cpDNA regions were sequenced: the non-coding region between the trnL (UAA) 5′exon and trnF (GAA), and the trnK region (including matK). Thirteen distinct cpDNA haplotypes were recognized and each haplotype was found to be geographically structured. Two major clades (I and II+III) were revealed in phylogenetic analyses among the haplotypes using F. sylvatica as an outgroup. The haplotypes of Clade I were distributed mainly along the Japan Sea side of the Japanese Archipelago, while those of Clade II+III occurred chiefly along the Pacific Ocean side. Consequently, the distribution of the two major cpDNA clades suggests that there were two migration routes in the history of F. crenata; one along the Japan Sea and the other along the Pacific Ocean side of the Japanese Islands. Received March 19, 2001 Accepted November 22, 2001     

Comparative phylogeography of unglaciated eastern North America

Regional phylogeographical studies involving co-distributed animal and plant species have been conducted for several areas, most notably for Europe and the Pacific Northwest of North America. Until recently, phylogeographical studies in unglaciated eastern North America have been largely limited to animals. As more studies emerge for diverse lineages (including plants), it seems timely to assess the phylogeography across this region: (i) comparing and contrasting the patterns seen in plants and animals; (ii) assessing the extent of pseudocongruence; and (iii) discussing the potential applications of regional phylogeography to issues in ecology, such as response to climatic change. Unglaciated eastern North America is a large, geologically and topographically complex area with the species examined having diverse distributions. Nonetheless, some recurrent patterns emerge: (i) maritime - Atlantic vs. Gulf Coast; (ii) Apalachicola River discontinuity; (iii) Tombigbee River discontinuity; (iv) the Appalachian Mountain discontinuity; (v) the Mississippi River discontinuity; and (vi) the Apalachicola River and Mississippi River discontinuities. Although initially documented in animals, most of these patterns are also apparent in plants, providing support for phylogeographical generalizations. These patterns may generally be attributable to isolation and differentiation during Pleistocene glaciation, but in some cases may be older (Pliocene). Molecular studies sometimes agree with longstanding hypotheses of glacial refugia, but also suggest additional possible refugia, such as the southern Appalachian Mountains and areas close to the Laurentide Ice Sheet. Many species exhibit distinct patterns that reflect the unique, rather than the shared, aspects of species' phylogeographical histories. Furthermore, similar modern phylogeographical patterns can result from different underlying causal factors operating at different times (i.e. pseudocongruence). One underemphasized component of pseudocongruence may result from the efforts of researchers to categorize patterns visually - similar patterns may, in fact, not fully coincide, and inferring agreement may obscure the actual patterns and lead to erroneous conclusions. Our modelling analyses indicate no clear spatial patterning and support the hypothesis that phylogeographical structure in diverse temperate taxa is complex and was not shaped by just a few barriers.     

Chloroplast DNA phylogeography and cytotype geography in autopolyploid Plantago media

In order to gain insight into the causes of parapatric diploid and tetraploid distributions in Plantago media chloroplast DNA (cpDNA) restriction site polymorphism was studied in 36 European populations. Parapatric distributions are often explained by adaptive differences between cytotypes to an underlying heterogeneity in environmental factors. Alternatively, such distribution patterns may be explained nonadaptively, through frequency-dependant production of hybrids with low fitness. However, nonadaptive explanations have been neglected in polyploid literature. In this study nine chloroplast haplotypes were found. Their phylogeny suggests that tetraploids arose at least three times from diploids. In general, related haplotypes were also geographically clustered, although there were some marked geographical discontinuities. In the Pyrenees, diploids and tetraploids carried diverged haplotypes throughout their parapatric ranges. At the contact zone the level of cpDNA introgression in a mixed diploid-tetraploid population was low. It is discussed that the cpDNA phylogeography supports the nonadaptive hypothesis that parapatric cytotype distributions may be explained by postPleistocene range expansions followed by mutual minority cytotype exclusion, due to hybrid unfitness.     

Chloroplast DNA phylogeography of the shrub Cistus ladanifer L. (Cistaceae) in the highly diverse Western Mediterranean region

This study investigated the phylogeographic structure of Cistus ladanifer, in order to locate its Quaternary refugia, reconstruct its recolonisation patterns and assess the role of geographical features (mountain ranges, rivers and the Strait of Gibraltar) as barriers to its seed flow and expansion through the Western Mediterranean. Thirty-eight populations were screened for length variation of polymorphic chloroplast simple sequence repeats (cpSSRs). Statistical analyses included estimation of haplotypic diversity, hierarchical analysis of molecular variation (amova) and fixation indices. Mantel tests, SAMOVA and BARRIER analyses were applied to evaluate the geographical partitioning of genetic diversity across the entire species range. Pollen data from bibliography were used to complement molecular inferences. Chlorotype diversity within populations was similar throughout the natural range of C. ladanifer (mean haplotypic diversity=0.32). High differentiation among populations was estimated (G(ST)=0.60). Our data suggest that the barriers of the Strait of Gibraltar and the Betic ranges may have favoured the divergence during glacial periods of four different lineages of populations inferred with SAMOVA. The main northward colonisation of in the Iberian Peninsula occurred from refugia in southwest Iberia. This process may have been influenced by human activities (forest clearance, livestock grazing and even commerce) in the Iberian Peninsula. In contrast, populations in the Betic area have conserved a specific haplotype.     

Cytoplasmic DNA variation and biogeography of Larix Mill. in Northeast Asia

Range‐wide variation in 54 populations of Dahurian larch (Larix gmelinii) and related taxa in Northeast Asia was assessed with four mitochondrial PCR‐RFLP and five chloroplast SSR markers. Eleven mitotypes and 115 chlorotypes were detected. The highest diversity was observed in the southern Russian Far East where hybrids of L. gmelinii, L. olgensis and L. kamtschatica are distributed. In contrast, only two mitotypes occurred in L. cajanderi and L. gmelinii. The Japanese larch (L. kaempferi) was found to be closely related to populations of L. kamtschatica inhabiting the Kuril Islands and South Sakhalin, populations from the northern part of Sakhalin being more closely related to continental species. In general, both mitochondrial (G_ST = 0.786; N_ST = 0.823) and chloroplast (G_ST = 0.144; R_ST = 0.432) markers showed a strong phylogeographical structure and evidence of isolation‐by‐distance. Yet both markers did not allow a clear delineation of species borders. In particular, and contrary to expectations, cpDNA was not significantly better than mtDNA at delineating species borders. This lack of concordance between morphological species and molecular markers could reflect extensive ancestral haplotype sharing and past and ongoing introgression. Finally the distribution of mtDNA and cpDNA variation suggests the presence of several refugia during Pleistocene glacial intervals. In particular, mtDNA and cpDNA reveal weak but visible differentiation between L. gmelinii and L. cajanderi, suggesting independent glacial histories of these species.     

Glacial vicariance in the Pacific Northwest: evidence from a lodgepole pine mitochondrial DNA minisatellite for multiple genetically distinct and widely separated refugia

The Canadian side of the Pacific Northwest was almost entirely covered by ice during the last glacial maximum, which has induced vicariance and genetic population structure for several plant and animal taxa. Lodgepole pine ( Pinus contorta Dougl. ex. Loud.) has a wide latitudinal and longitudinal distribution in the Pacific Northwest. Our main objective was to identify relictual signatures of glacial vicariance in the population structure of the species and search for evidence of distinct glacial refugia in the Pacific Northwest. A maternally inherited mitochondrial DNA minisatellite-like marker was used to decipher haplotype diversity in 91 populations of lodgepole pine located across the natural range. Overall population differentiation was sizeable ( G _ST = 0.365 and R _ST = 0.568). Four relatively homogeneous groups of populations, possibly representative of as many genetically distinct glacial populations, were identified for the two main subspecies, ssp. latifolia and ssp. contorta . For ssp. contorta , one glacial lineage is suggested to have been located at high latitudes and possibly off the coast of mainland British Columbia (BC), while the other is considered to have been located south of the ice sheet along the Pacific coast. For ssp. latifolia , two genetically distinct glacial populations probably occurred south of the ice sheet: in the area bounded by the Cascades and Rocky Mountains ranges, and on the eastern side of the Rockies. A possible fifth refugium located in the Yukon may have also been present for ssp. latifolia . Zones of contact between these ancestral lineages were also apparent in interior and northern BC. These results indicate the role of the Queen Charlotte Islands and the Alexander Archipelago as a refugial zone for some Pacific Northwest species and the vicariant role played by the Cascades and the American Rocky Mountains during glaciation.     

Of glaciers and refugia: a decade of study sheds new light on the phylogeography of northwestern North America

Glacial cycles have played a dominant role in shaping the genetic structure and distribution of biota in northwestern North America. The two major ice age refugia of Beringia and the Pacific Northwest were connected by major mountain chains and bordered by the Pacific Ocean. As a result, numerous refugial options were available for the regions taxa during glacial advances. We reviewed the importance of glaciations and refugia in shaping northwestern North America’s phylogeographic history. We also tested whether ecological variables were associated with refugial history. The recurrent phylogeographic patterns that emerged were the following: (i) additional complexity, i.e. refugia within refugia, in both Beringia and the Pacific Northwest; and (ii) strong evidence for cryptic refugia in the Alexander Archipelago and Haida Gwaii, the Canadian Arctic and within the ice‐sheets. Species with contemporary ranges that covered multiple refugia, or those with high dispersal ability, were significantly more likely to have resided in multiple refugia. Most of the shared phylogeographic patterns can be attributed to multiple refugial locales during the last glacial maximum or major physiographic barriers like rivers and glaciers. However, some of the observed patterns are much older and appear connected to the orogeny of the Cascade‐Sierra chain or allopatric differentiation during historic glacial advances. The emergent patterns from this review suggest we should refine the classic Beringian‐southern refugial paradigm for northwestern North American biota and highlight the ecological and evolutionary consequences of colonization from multiple refugia.     

Evidence for multiple sources of invasion and intraspecific hybridization in Brachypodium sylvaticum (Hudson) Beauv. in North America

We compared the levels and distribution of genetic diversity in Eurasian and North American populations of Brachypodium sylvaticum (Huds.) Beauv. (false brome), a newly invasive perennial bunchgrass in western North America. Our goals were to identify source regions for invasive populations, determine the number of independent invasion events, and assess the possibility that postinvasion bottlenecks and hybridization have affected patterns of genetic diversity in the invaded range. We tested the hypothesis that this Eurasian grass was accidentally introduced into two areas in Oregon and one site in California by examining nuclear microsatellites and chloroplast haplotype variation in 23 introduced and 25 native populations. In the invaded range, there was significantly lower allelic richness (R(S)), observed heterozygosity (H(O)) and within-population gene diversity (H(S)), although a formal test failed to detect a significant genetic bottleneck. Most of the genetic variation existed among populations in the native range but within populations in the invaded range. All of the allelic variation in the invaded range could be explained based on alleles found in western European populations. The distribution of identified genetic clusters in the North American populations and the unique alleles associated with them is consistent with two historical introductions in Oregon and a separate introduction to California. Further analyses of population structure indicate that intraspecific hybridization among genotypes from geographically distinct regions of western Europe occurred following colonization in Oregon. The California populations, however, are more likely to be derived from one or perhaps several genetically similar regions in the native range. The emergence and spread of novel recombinant genotypes may be facilitating the rapid spread of this invasive species in Oregon.     

Glacial refugia and the phylogeography of Steller's sea lion (Eumatopias jubatus) in the North Pacific

Mitochondrial DNA sequence data were used to examine the phylogeographic history of Steller's sea lions (Eumetopias jubatus) in relation to the presence of Plio-Pleistocene insular refugia. Cytochrome b and control region sequences from 336 Steller's sea lions reveal phylogenetic lineages associated with continental refugia south of the ice sheets in North America and Eurasia. Phylogenetic analysis suggests the genetic structure of E. jubatus is the result of Pleistocene glacial geology, which caused the elimination and subsequent reappearance of suitable rookery habitat during glacial and interglacial periods. The cyclic nature of geological change produced a series of independent population expansions, contractions and isolations that had analogous results on Steller's sea lions and other marine and terrestrial species. Our data show evidence of four glacial refugia in which populations of Steller's sea lions diverged. These events occurred from approximately 60,000 to 180,000 years BP and thus preceded the last glacial maximum.     

Comparative phylogeography of two North American 'glacial relict' crustaceans

The Pleistocene glaciations represent the most recent and dramatic series of habitat changes since the Cretaceous. The impact of these events was particularly acute for aquatic taxa with poor powers of dispersal, but few organisms have evolutionary histories more intimately entwined with the advance and retreat of ice than the 'glacial relicts'. In this study, we used a mitochondrial gene, cytochrome c oxidase subunit I (COI), to examine and compare the phylogeographical structure of two glacial relict crustaceans (Limnocalanus macrurus and members of the Mysis relicta species group) across North America. In both cases, we found a sharp phylogenetic division between populations from inland lakes formed during glacial retreat, and arctic lakes isolated from polar seas via isostatic rebound. However, the depth of this phylogenetic partition varied between taxa. In L. macrurus, nucleotide sequence divergence of 2.2% between these zones is consistent with its current status as a single morphologically variable species, but in Mysis the split occurred among recently described, morphologically conserved species, at a divergence of 8.2%. The disparity in the depth of divergence indicates a history of recurrent freshwater invasions from the arctic seas, in concordance with previous studies of Eurasian glacial relicts. However, we suggest further consideration of a largely overlooked explanation that could account for some of the discrepancies between molecular divergences and glaciation events. Many cladogenetic events could have occurred in arctic seas prior to the transition to inland waters, a possibility supported both by the complex physical and ionic history of arctic seas and by high marine and estuarine lineage diversity in the north.     

Phylogeography of ninespine sticklebacks (Pungitius pungitius) in North America: glacial refugia and the origins of adaptive traits

The current geographical distribution of the ninespine stickleback (Pungitius pungitius) was shaped in large part by the glaciation events of the Pleistocene epoch (2.6 Mya–10 Kya). Previous efforts to elucidate the phylogeographical history of the ninespine stickleback in North America have focused on a limited set of morphological traits, some of which are likely subject to widespread convergent evolution, thereby potentially obscuring relationships among populations. In this study, we used genetic information from both mitochondrial DNA (mtDNA) sequences and nuclear microsatellite markers to determine the phylogenetic relationships among ninespine stickleback populations. We found that ninespine sticklebacks in North America probably dispersed from at least three glacial refugia—the Mississippi, Bering, and Atlantic refugia—not two as previously thought. However, by applying a molecular clock to our mtDNA data, we found that these three groups diverged long before the most recent glacial period. Our new phylogeny serves as a critical framework for examining the evolution of derived traits in this species, including adaptive phenotypes that evolved multiple times in different lineages. In particular, we inferred that loss of the pelvic (hind fin) skeleton probably evolved independently in populations descended from each of the three putative North American refugia.     

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed‐dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3–4) and paternally inherited, pollen‐ and seed‐dispersed chloroplast DNA (cpDNA) (rpl16 and trnS‐trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G_ST = 0.738, N_ST = 0.771) indicating low levels of seed‐based gene flow and significant phylogeographical structure (N_ST > G_ST, P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north–central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north‐central part of the species’ range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north‐central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among‐population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.     

From glacial refugia to modern populations: new assemblages of organelle genomes generated by differential cytoplasmic gene flow in transcontinental black spruce

Assessing species' range-wide cytoplasmic diversity provides valuable insights regarding their dispersal and adaptive potential in a changing environment. Transcontinental chloroplast (cpDNA) and mitochondrial DNA (mtDNA) population structures were compared to identify putative ancestral and new cytoplasmic genome assemblages in black spruce (Picea mariana), a North American boreal conifer. Mean within-population diversity and allelic richness for cpSSR markers were 0.80 and 4.21, respectively, and diminished westward. Population differentiation based on G(ST) was lower for cpDNA than for mtDNA (G(ST) =0.104 and 0.645, respectively) but appeared comparable when estimated using Jost differentiation index (D=0.459 and 0.537, respectively). Further analyses resulted in the delineation of at least three genetically distinct cpDNA lineages partially congruent with those inferred from mtDNA data, which roughly corresponded to western, central and eastern Canada. Additionally, the patterns of variation in Alaska for both cpDNA and mtDNA markers suggested that black spruce survived the last glacial maximum in this northern region. The range-wide comparison of the geographic extent of cytoplasmic DNA lineages revealed that extensive pollen gene flow between ancestral lineages occurred preferentially from west to east during the postglacial expansion of the species, while seed-mediated gene flow remained geographically restricted. This differential gene flow promoted intraspecific cytoplasmic capture that generated new assemblages of cpDNA and mtDNA genomes during the Holocene. Hence, black spruce postglacial colonization unexpectedly resulted in an increase in genetic diversity with possible adaptive consequences.     

The evolutionary history of a mammal species with a highly fragmented range: the phylogeography of the European snow vole

The European snow vole Chionomys nivalis has a patchy distribution restricted to rocky habitats across southern Europe and the Near and Middle East. We carried out a phylogeographic study to provide a biogeographic scenario, based on molecular data, outlining the major processes that determined the current distribution of the species. The samples include 26 snow voles from 14 different populations across the entire species range from Spain to Anatolia and Israel. Nearly complete sequences (1037 bp) of the mitochondrial gene for cytochrome b were sequenced. Relationships among haplotypes were inferred with neighbour-joining, maximum likelihood, maximum parsimony analyses and minimum spanning network. An analysis of mismatch distribution was used to cast light on past demographic expansion. We found 22 different haplotypes that fall into six distinct lineages, all but one is supported by high bootstrap values with all methods. Four lineages are allopatric (Tatra Mts., Iberia, Balkans and Middle East) while divergent haplotypes from two lineages show sympatry in the Alps and the Apennines. The basal relationships of these lineages could not be established by any tree. The mean pairwise genetic distance between lineages ranges from 2.4 to 4.2%. The shape of the mismatch distribution indicated a past expansion event dating back to between 158 000 and 84 000 years ago. These data can be interpreted with the existence of southern glacial refugia (Iberia, Balkans, Middle East and Italy) and one additional northern glacial refugium. The lack of phylogenetic resolution among lineages and the shape of mismatch distribution are indicative of a simultaneous and rapid splitting due to a relatively fast initial expansion of populations. Moreover, the analysis supports the hypothesis of the European origin of C. nivalis and its subsequent eastward dispersion during the Middle Pleistocene.     

Chloroplast DNA phylogeography suggests a West African centre of origin for the baobab, Adansonia digitata L. (Bombacoideae, Malvaceae)

The African baobab ( Adansonia digitata L.) is an emblematic, culturally important, and physically huge tropical tree species whose natural geographical distribution comprises most of tropical Africa, but also small patches of southern Arabia, and several Atlantic and Indian Ocean islands surrounding the African continent, notably including Madagascar. We analysed the polymerase chain reaction–restriction fragment length polymorphism of five chloroplast DNA fragments obtained from 344 individuals of A. digitata collected from 74 populations covering the entire extant distribution range of the species. Our goal was to reconstruct the phylogeographical history of the species and, if possible, to identify its centre of origin, which has been a subject of controversy for many decades. We identified five haplotypes whose distribution is clearly geographically structured. Using several species of Adansonia and of closely related genera as outgroups, the haplotypes showed a clear phylogeographical pattern of three groups. Two are phylogenetically related to the outgroup taxa, and are distributed in West Africa. The third group is substantially more differentiated genetically from outgroup species, and it corresponds to southern and eastern Africa, Arabia and the Indian Ocean islands, including Madagascar. According to our results, the tetraploid A. digitata , or its diploid progenitor, probably originated in West Africa and migrated subsequently throughout the tropical parts of that continent, and beyond, by natural and human-mediated terrestrial and overseas dispersal.     

Worldwide phylogeography of the blacktip shark (Carcharhinus limbatus) inferred from mitochondrial DNA reveals isolation of western Atlantic populations coupled with recent Pacific dispersal

Although many coastal shark species have widespread distributions, the genetic relatedness of worldwide populations has been examined for few species. The blacktip shark, (Carcharhinus limbatus), inhabits tropical and subtropical coastal waters throughout the world. In this study, we examined the genetic relationships of blacktip shark populations (n = 364 sharks) throughout the majority of the species' range using the entire mitochondrial control region (1067-1070 nucleotides). Two geographically distinct maternal lineages (western Atlantic, Gulf of Mexico, and Caribbean Sea clades, and eastern Atlantic, Indian, and Pacific Ocean clades) were identified and shallow population structure was detected throughout their geographic ranges. These findings indicate that a major population subdivision exists across the Atlantic Ocean, but not the Pacific Ocean. The historical dispersal of this widespread, coastal species may have been interrupted by the rise of the Isthmus of Panama. This scenario implies historical dispersal across the Pacific Ocean (supported by the recovery of the same common haplotype from the Philippines, Hawaii, and the Gulf of California reflecting recent/contemporary dispersal abilities) and an oceanic barrier to recent migration across the Atlantic. Genetic structure within the eastern Atlantic/Indo-Pacific (Phi(ST) = 0.612, P < 0.001) supports maternal philopatry throughout this area, expanding previous western Atlantic findings. Eastern Atlantic/Indo-Pacific C. limbatus control region haplotypes were paraphyletic to Carcharhinus tilstoni haplotypes in our maximum-parsimony analysis. The greater divergence of western Atlantic C. limbatus than C. tilstoni from eastern Atlantic/Indo-Pacific C. limbatus reflects the taxonomic uncertainty of western Atlantic C. limbatus.