The postglacial history of Scots pine (Pinus sylvestris L.) in western Europe: evidence from mitochondrial DNA variation

The geographical structure of mitochondrial (mt)DNA variants (mitotypes) was investigated in 38 western European populations of Scots pine Pinus sylvestris using restriction fragment length polymorphism (RFLP) analysis of total DNA and a homologous cox1 probe. Three major mitotypes (designated a, b and d ) were detected. Within Spain all three major mitotypes were found, gene diversity was high, H_T = 0.586, and this diversity was distributed predominantly among rather than within populations (F_ST(M) = 0.813 for the seven Spanish populations). Mitotype d was present only in the most southerly population from the Sierra Nevada . Elsewhere in Europe, populations showed little or no mtDNA diversity within regions, but there were marked differences between regions. Italian populations were fixed for mitotype b ; populations from northern France, Germany, Poland, Russia and southern Sweden were fixed for mitotype a ; while populations in northern Fennoscandia were fixed for mitotype b . The isolated Scottish populations were predominantly of mitotype a , but mitotype b was present in three of the 20 populations scored. In Scotland, UK gene diversity (H_T = 0.120) and genetic differentiation among populations (F_ST(M) = 0.37) was much lower than in Spain. When interpreted in the light of complementary data from pollen analysis and nuclear genetic markers, the results suggest that present-day populations of P. sylvestris in western Europe have been derived from at least three different sources after glaciation.     

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.     

Mitochondrial and chloroplast phylogeography of Picea crassifolia Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands

The disjunct distribution of forests in the Qinghai-Tibetan Plateau (QTP) and adjacent Helan Shan and Daqing Shan highlands provides an excellent model to examine vegetation shifts, glacial refugia and gene flow of key species in this complex landscape region in response to past climatic oscillations and human disturbance. In this study, we examined maternally inherited mitochondrial DNA (nad1 intron b/c and nad5 intron 1) and paternally inherited chloroplast DNA (trnC-trnD) sequence variation within a dominant forest species, Picea crassifolia Kom. We recovered nine mitotypes and two chlorotypes in a survey of 442 individuals from 32 populations sampled throughout the species' range. Significant mitochondrial DNA population subdivision was detected (G(ST) = 0.512; N(ST) = 0.679), suggesting low levels of recurrent gene flow through seeds among populations and significant phylogeographical structure (N(ST) > GST, P < 0.05). Plateau haplotypes differed in sequence from those in the adjacent highlands, suggesting a long period of allopatric fragmentation between the species in the two regions and the presence of independent refugia in each region during Quaternary glaciations. On the QTP platform, all but one of the disjunct populations surveyed were fixed for the same mitotype, while most populations at the plateau edge contained more than one haplotype with the mitotype that was fixed in plateau platform populations always present at high frequency. This distribution pattern suggests that present-day disjunct populations on the QTP platform experienced a common recolonization history. The same phylogeographical pattern, however, was not detected for paternally inherited chloroplast DNA haplotypes. Two chlorotypes were distributed throughout the range of the species with little geographical population differentiation (G(ST) = N(ST) = 0.093). This provides evidence for highly efficient pollen-mediated gene flow among isolated forest patches, both within and between the QTP and adjacent highland populations. A lack of isolation to pollen-mediated gene flow between forests on the QTP and adjacent highlands is surprising given that the Tengger Desert has been a geographical barrier between these two regions for approximately the last 1.8 million years.     

Polymorphism of mitochondrial DNA and infection with symbiotic cytoplasmic bacterium Wolbachia pipientis in mosquitoes of the Culex pipiens complex from Russia

A total of 208 mosquitoes of the Culex pipiens complex from 15 basement and terrestrial populations collected in different regions of the European part of Russia and Siberia were examined by genetic methods. Among these, two major mitotypes, M and P, were identified. These mitotypes differed by six substitutions in the 246-bp mitochondrial DNA cytochrome oxidase I gene fragment examined. Urban basement mosquito ecotype C. pipiens (form molestus) were characterized by the presence of mitotype M and infection with the endosymbiotic bacteria of the genus Wolbachia. Mosquitoes of the C. pipiens complex inhabiting opened biotopes harbored mitotype P, or its variety, mitotype P1, and were not infected with Wolbachia. Thus, in natural conditions marked linkage disequilibrium between cytoplasmic elements, mitochondrial DNA and Wolbachia, can be observed. Similarity of mitotypes in form molestus mosquito from different geographical localities favors the hypothesis on the common ancestry of urban mosquitoes.     

Chloroplast DNA markers reveal a geographical divide across Argentinean southern beech Nothofagus nervosa (Phil.) Dim. et Mil. distribution area

  Nothofagus nervosa is one of the most important species of the temperate forests of southern South America. On the eastern slope of the Andes Mountains (Argentina) it has a very small and narrow natural distribution area which follows the valleys of the numerous west-east lake basins. Re-colonisation after the last glaciation is assumed to have originated from refugia located in the western pacific coast or even in free interglacial patches of the Andes Mountains. Sixteen pairs of primers were used to amplify 16 non-coding regions of chloroplast DNA in 11 populations from over the entire distribution area in Argentina. After restriction analysis two polymorphic fragments were found which defined two haplotypes and allowed a differentiation among populations. A clear geographic divide occurs separating populations distributed south and north of a west-east mountain chain. This suggests that after the last glaciation the species spread out from at least two different refugia. Hypotheses about the possible location of these refugia are discussed. Received: 28 October 1997 / Accepted: 19 March 1998     

Morphometric variation in the Nearctic collared lemming (Dicrostonyx)

A. H. Macpherson suggested that much of the current geographic diversity in Canadian Arctic mammals resulted from isolation in refugia during the Wisconsin glacial stage. This study evaluates the refugium hypothesis, insofar as it applies to Nearctic Dicrostonyx , by means of a statistical analysis of geographic variation in 13 skull characters. Overall, geographic variation among samples is not significant, although D. hudsonius and D. unalascensis are geographically and morphologically distinct. Some variation in skull shape is correlated with winter temperature. Partitioning tests on other measures of shape variation indicate some discontinuities consistent with the refugial hypothesis. Discrete samples reflect possible refugial populations in northern North America, Eastern Beringia and two southern periglacial refugia, one in eastern North America and at least one in western North America.     

DNA footprints of European hedgehogs, Erinaceus europaeus and E. concolor: Pleistocene refugia, postglacial expansion and colonization routes

European hedgehogs, Erinaceus europaeus and E. concolor, are among the many European plant and animal taxa that have been subjected to cyclical restriction to glacial refugia and interglacial expansion. An analysis of 95 mitotypes, comprising partial cytochrome b and control region sequences, shows deep divergence between the two hedgehog species. Three europaeus and two concolor clades are clearly identified and are consistent with previously identified refugia for Europe: the Iberian peninsula, Italy, and the Balkans. The degree of mitochondrial divergence among these clades suggests pre-Pleistocene separation of the refugial populations. In contrast, analysis of two nuclear introns clearly separates the two concolor clades, as in the mitochondrial data, but cannot discriminate the three europaeus clades. This discrepancy between nuclear and mitochondrial data is attributed to historical differences in the refugial population size of europaeus and concolor. The geographical distribution of mitotypes is analysed using nested clade analysis. This method, by including unobserved ('missing') mitotypes, can identify mitotype groupings that remain undetected in conventional analyses. However, the application of nested clade analysis to the study of refugial populations may be hampered by such factors as the loss of haplotypes from the refugial areas by repeated contractions of the population and the recent time scale of colonization relative to mutation rate.     

Patterns of genetic variation in Pinus chiapensis,a threatened Mexican pine,detected by RAPD and mitochondrial DNA RFLP markers

Pinus chiapensis (Pinaceae) is a large conifer, endemic to central and southern Mexico and north-western Guatemala. In order to assess the extent of genetic variation within and between populations of this species, samples were obtained from throughout the natural range and analysed using random amplified polymorphic DNA (RAPD) and mtDNA RFLPs markers. Probes for the CoxI mitochondrial gene enabled two mitotypes to be observed. Populations from the eastern and western limit of the range of the species were fixed for one mitotype ('A'), whereas two populations distributed near the centre of the range were fixed for another ('B'). When the samples were screened with eight 10-mer RAPD primers, a total of 12 polymorphic bands were detected. The proportion of polymorphic bands was unusually low (24.5%) compared with other tree species. AMOVA analysis indicated that a significant proportion of the variation (P < 0.002) was distributed between populations; the extent of population differentiation detected (Phi(st) = 0.226; G(ST ) = 0.194) was exceptionally high for a pine species. Pair-wise comparison of Phi(st) values derived from AMOVA indicated that populations were significantly (P < 0.05) different from each other in virtually every case. These results are interpreted in the context of the evolutionary history of the species, and the implications for its in- and ex situ conservation are discussed.     

Genetic divergence and phylogeographic relationships among european perch (Perca fluviatilis) populations reflect glacial refugia and postglacial colonization

We used the widely distributed freshwater fish, perch (Perca fluviatilis), to investigate the postglacial colonization routes of freshwater fishes in Europe. Genetic variability within and among drainages was assessed using mitochondrial DNA (mtDNA) D-loop sequencing and RAPD markers from 55 populations all over Europe as well as one Siberian population. High level of structuring for both markers was observed among drainages and regions, while little differentiation was seen within drainages and regions. Phylogeographic relationships among European perch were determined from the distribution of 35 mtDNA haplotypes detected in the samples. In addition to a distinct southern European group, which includes a Greek and a southern Danubian population, three major groups of perch are observed: the western European drainages, the eastern European drainages including the Siberian population, and Norwegian populations from northern Norway, and western side of Oslofjord. Our data suggest that present perch populations in western and northern Europe were colonized from three main refugia, located in southeastern, northeastern and western Europe. In support of this, nested cladistic analysis of mtDNA clade and nested clade distances suggested historical range expansion as the main factor determining geographical distribution of haplotypes. The Baltic Sea has been colonized from all three refugia, and northeastern Europe harbours descendants from both eastern European refugia. In the upper part of the Danube lineages from the western European and the southern European refugia meet. The southern European refugium probably did not contribute to the recolonization of other western and northern European drainages after the last glaciation. However, phylogenetic analyses suggest that the southern European mtDNA lineage is the most ancient, and therefore likely to be the founder of all present perch lineages. The colonization routes used by perch probably also apply to other freshwater species with similar distribution patterns.     

Genetic differentiation among populations of Cicerbita alpina (L.) Wallroth (asteraceae) in the western Alps

In this study we analyzed the genetic population structure of the hygrophilous tall-herb Cicerbita alpina in the western Alps because this group of mountain plants is underrepresented in the biogeographical literature. AFLP (amplified fragment length polymorphism) fingerprints of 40 samples were analyzed from four populations situated in a transect from the southwestern Alps to the eastern part of the western Alps and one population from the Black Forest outside the Alps. Two genetic groups can be distinguished. The first group (A) comprises the populations from the northern and eastern parts of the western Alps, and the second group (B) comprises the populations from the southwestern Alps and the Black Forest. Group A originates most likely from at least one refugium in the southern piedmont regions of the Alps. This result provides molecular evidence for a humid climate at the southern margin of the Alps during the Würm glaciation. Group B originates presumably from western or northern direction and we discuss two possible scenarios for the colonization of the Alps, i.e. (1) long-distance dispersal from southwestern refugia and (2) colonization from nearby refugia in the western and/or northern Alpine forelands. The study demonstrates that the target species harbours considerable genetic diversity, even on a regional scale, and therefore is a suitable model for phylogeographic research.     

Polymorphism of mitochondrial DNA and infection with symbiotic cytoplasmic bacterium <Emphasis Type="Italic">Wolbachia pipientis</Emphasis> in mosquitoes of the <Emphasis Type="Italic">Culex pipiens</Emphasis> (Diptera,Culicidae) complex from Russia

A total of 208 mosquitoes of the Culex pipiens complex from 15 basement and terrestrial populations collected in different regions of the European part of Russia and Siberia were examined by genetic methods. Among these, two major mitotypes, M and P, were identified. These mitotypes differed by six substitutions in the 246-bp mitochondrial DNA cytochrome oxidase I gene fragment examined. Urban basement mosquitoes C. pipiens (form molestus) were characterized by the presence of mitotype M and infection with the endosymbiotic bacteria of the genus Wolbachia. Mosquitoes of the C. pipiens complex inhabiting opened biotopes harbored mitotype P, or its variety, mitotype P1, and were not infected with Wolbachia. Thus, in natural conditions marked linkage disequilibrium between cytoplasmic elements, mitochondrial DNA and Wolbachia, can be observed. Similarity of mitotypes in form molestus mosquito from different geographical localities favors the hypothesis on the common ancestry of urban mosquitoes.Translated from Genetika, Vol. 41, No. 3, 2005, pp. 320–325.Original Russian Text Copyright © 2005 by Shaikevich, Vinogradova, Platonov, Karan, Zakharov.     

Late Pleistocene divergence between eastern and western populations of wood ducks (Aix sponsa) inferred by the 'isolation with migration' coalescent method

During the Late Pleistocene, glaciers sundered many species into multiple glacial refugia where populations diverged in allopatry. Although deeply divergent mitochondrial DNA (mtDNA) lineages often reflect the number of refugia occupied, it is unlikely that populations that split during the recent Wisconsin glaciations will have reached reciprocal monophyly. We examined mtDNA control region sequences from eastern and western populations of wood ducks (Aix sponsa) to determine whether their current, disjunct distribution is consistent with the occupancy of two glacial refugia. We used the 'isolation with migration' coalescent method (im) to simultaneously estimate effective population sizes, maternal gene flow, and time since divergence. We found 24 unique haplotypes, none of which were shared between the eastern and western populations, but we did not find diagnostic monophyletic lineages suggestive of long-term isolation in multiple glacial refugia. However, a high Phi ST (0.31) indicates that eastern and western populations are well differentiated in mtDNA, and results from im suggest that these populations have been diverging, without extensive gene flow, for 10,000 to 124,000 years. Results from im further suggest that these populations most likely split about 34,000 years ago, and this time of divergence is consistent with the occupancy of multiple glacial refugia during the Late Wisconsin glaciation. Eastern wood ducks are characterized by high genetic diversity, a large effective population size, and a recent population expansion, while western wood ducks have much less genetic diversity, a smaller population size, and have not undergone a recent population expansion.     

MULTIPLE ORIGINS OF GENDER-ASSOCIATED MITOCHONDRIAL DNA LINEAGES IN BIVALVES (MOLLUSCA: BIVALVIA)

Previous studies have shown that marine mussels (genus Mytilus) and a freshwater mussel (Pyganodon grandis) contain two distinct gender-associated mitotypes, which is a characteristic feature of the phenomenon of doubly uniparental inheritance (DUI) of mitochondrial DNA (mtDNA). Here we present evidence for the presence of distinct male (M) and female (F) mitotypes in three other bivalve species, the mytilid Geukensia demissa, and the unionid species P. fragilis and Fusconaia flava. Nucleotide sequences of a segment of the COI gene from the M and F mitotypes from each of the three mytilid species (M. edulis, M. trossulus, G. demissa) and three unionid species (P. grandis, P. fragilis, F. flava) were used for phylogenetic analysis. The analysis suggests three independent origins of M and F mitotypes for the six species examined; one for the three unionid species, one for the two Mytilus species, and one for Geukensia. The first of these F/M divergence events, while of uncertain age, predates the divergence of the two unionid genera and is likely older than either of the two F/M divergence events in the mytilid taxa. The most parsimonious explanation of multiple F/M divergence events is that they represent independent origins of DUI. Another possibility is that, in a given taxon, an F or M mitotype assumes the role of the opposite mitotype (by virtue of a mechanism that remains to be clarified) and subsequently was fixed within its new gender. The fixation of a mtDNA lineage derived from a mitotype of switched function would reset the divergence of the gender-associated lineages to zero, thereby mimicking a de novo split of F and M lineages from a preexisting mtDNA genome that was not gender specific. Further broad-scale taxonomic studies of the occurrence of distinct M and F mitotypes may allow for the evaluation of the latter hypothesis.     

Mitochondrial DNA of pre‐last glacial maximum red deer from NW Spain suggests a more complex phylogeographical history for the species

The major climatic oscillations that characterized the Quaternary had a great influence on the evolution and distribution of several species. During cold periods, the distribution of temperate‐adapted species became fragmented with many surviving in southern refugia (Iberian, Italian, and Balkan Peninsulas). Red deer was one of the species that contracted its original range to southern refugia. Currently, two main lineages have been described for the species: western and eastern. We have analyzed fossils pre‐dating the last glacial maximum (LGM) from Liñares cave (NW Spain) that belongs to the peripheral range of the western clade, and fossils from the Danish Holocene belonging to the central part of the same clade. Phylogenetic analyses place our samples in the western clade. However, some specimens from Liñares represent an early split in the tree along with other pre‐LGM western samples from previous studies. Despite low bootstrap values in the Bayesian phylogenies, haplotype networks connect these foreign haplotypes to the eastern clade. We suggest a mixed phylogeographical model to explain this pattern with range expansions from the east during the expansion phase after the cold periods in marine isotope stage 3. We find slight isolation by distance in post‐LGM populations that could be a consequence of the recolonization from southern refugia after the LGM.     

Comparative Restriction Enzyme Analysis of Mitochondrial DNA in Several Populations of Lake Resident Chars of the Genus Salvelinus

Restriction enzyme analysis was employed in studying the mitochondrial DNA (mtDNA) ATPase6/ND4L region in several Northeast Asian populations of resident lake char of the genus Salvelinus. On evidence of mitotypes, genetic similarity was assumed for populations of neiva (Ueginskoe Lake), lake resident char from the northern coast of the Sea of Okhotsk (Mak-Mak Lake and Elekchanskie lakes), and dolly varden. Mitotype AAAA proved to be common for these populations. Lake char populations of the Juliet and Maxi lakes (the basin of the Kolyma River) had mitotype DBAA, which is similar to mitotype DBAB observed earlier in Taranetz char from Chukotka. The divergence between mitotypes AAAA and DBAA was estimated at 0.3%. Different origins were assumed for the lake resident char populations from the basins of the Sea of Okhotsk and of the Kolyma River, the former originating in the Pacific and the latter, in the Arctic basins.     

The influence of late pleistocene mountain glaciations on the genetic differentiation of long-tailed ground squirrel (<Emphasis Type="Italic">Urocitellus undulatus</Emphasis>)

Long-tailed ground squirrel (Urocitellus undulatus) is a polytypic species with a wide distribution from the Tien Shan to the Amur River region. Previously, considerable genetic differentiation between eastern and western populations of this species was demonstrated. Moreover, the greatest differences were observed in the western part of the range located in Central Asia, the region that was subjected to repeated glaciations in the past and represents one of the centers of the ground squirrel secondary diversification. The analysis of polymorphism of the mitochondrial DNA control region was carried out on long-tailed ground squirrels living in the northern part of Central Asia, on the territory of the Altai Mountains (45 individuals from 23 localities). The presence of two genetically differentiated (7.7% differences) and geographically separated lineages (western and eastern) was revealed. The data obtained disprove the hypothesis on unidirectional, from west to east, colonization of the Altai Mountains after the end of the last glacial maximum and show the two pathways of the ground squirrel colonization of the Altai, from both western and eastern refugia.     

Maternal lineages of Pinus densata,a diploid hybrid

Previous morphological, allozyme and chloroplast DNA data have suggested that Pinus densata originated through hybridization between P. tabuliformis and P. yunnanensis. In the present study, sequence and restriction site analyses of maternally inherited mitochondrial nad1 intron were used to detect variation patterns in 19 populations of P. tabuliformis, P. yunnanensis and P. densata. A total of three mitotypes (A, B, C) were detected. All but one of the populations of P. yunnanensis possessed mitotype B while all populations of P. tabuliformis had mitotype A. Pinus densata populations, on the other hand, harboured both mitotypes A and B, which are characteristic of P. tabuliformis and P. yunnanensis, respectively. This result gives strong additional evidence supporting the hybrid origin of this diploid pine. The distribution of mitotypes indicated very different mating compositions and evolutionary history among P. densata populations. It seems that local founder populations and backcrosses may have played important roles in the early establishment of P. densata populations. The uplift of the Tibetan Plateau had a significant impact on the distribution of maternal lineages of P. densata populations.     

Phylogeography and genetic structuring of European nine-spined sticklebacks (Pungitius pungitius)-mitochondrial DNA evidence

As a consequence of colonisation from different glacial refugia, many northern European taxa are split into distinct western and eastern lineages. However, as for the nine-spined stickleback (Pungitius pungitius), the exact location of the contact zone between lineages often remains poorly known. We assessed the genetic differentiation and diversity in the nine-spined stickleback within Europe using 1037 base pairs of cytochrome b sequence for 320 individuals from 57 locations, including pond, lake, river, and coastal habitats. Our main aims were (i) to locate the contact zone between the previously recognized western and eastern lineages, (ii) investigate latitudinal patterns in genetic diversity, (iii) compare genetic diversity among different habitat types, and (iv) date the known split between eastern and western lineages. The data revealed the split between eastern and western to be located across the Danish Straits and roughly following the Norway/Sweden border to the North. Reference sites from Canada form their own clades, and one of the Canadian sites was found to have a haplotype common to the Eastern European lineage, possibly representing an ancestral polymorphism. The split between the two European clades was dated to approximately 1.48 million years ago (Mya), and between Canada and Europe to approximately 1.62 Mya. After controlling for habitat effects, nucleotide (but not haplotype) diversity across populations decreased with increasing latitude. Coastal populations showed significantly higher haplotype diversity (but not nucleotide diversity) than pond populations, but there were no detectable differences in haplotype diversity among different freshwater habitat types (viz. river, lake and pond populations), or between coastal and lake/river populations. Sequences were found to cluster according to their geographic proximity, rather than by habitat type, and all habitat types were found within each major clade, implying that colonisation and adaptation between the coastal and freshwater environments in different regions must have occurred in parallel.