Abstract
The variation at 15 allozyme loci was examined in 32 populations of Larix taxa across their range in northern Eurasia. Fifteen populations of L. sibirica Ledeb., six populations of L. gmelinii Rupr., three populations of L. olgensis A. Henry, two populations of L. decidua Mill., and one population each of L. cajanderi Mayr, L. kaempferi (Lamb.) Carr. (=L. leptolepis (Sieb. et Zucc.) Endl.), L. kamtschatica (Rupr.) Carr., L. czekanovskii Szaf., L. amurensis Kolesn. and L. ochotensis Kolesn. were analysed. Interpopulational genetic differentiation within taxa was estimated with Wright's FST and varied from 2.6 to 8.1% according to taxon. Highly significant allelic differentiation was detected with Fisher's exact test. Genetic distances among L. sibirica populations were associated primarily with longitude. Allozyme frequencies changed gradually from the Urals eastward to Siberia, and eastern populations of L. sibirica are genetically more similar to L. olgensis than to western L. sibirica. Far-eastern populations of Larix species, morphologically similar to L. gmelinii and defined as L. cajanderi, L. amurensis and L. ochotensis, appear to be genetically close to both L. gmelinii and to L. olgensis, suggesting that they originated through introgressive hybridization between these two species.
Similar content being viewed by others
Introduction
The present distribution of genetic variation in larch in Eurasia may, in large part, reflect the way the species spread and interacted over the postglacial recolonization period. Unlike northern Europe which was entirely covered by ice at the time of the last glacial maximum (18 kya), northern Siberia and the Russian far east were only locally glaciated (Bennett, 1997). This was recently confirmed indirectly by molecular studies in Arctic lemmings (Dicrostonyx); using mitochondrial DNA, Fedorov et al. (in press) revealed a complex phylogeographical pattern reflecting the past subdivisions of the present range by ice sheets. They also showed that the formerly glaciated Verchoyansky Mountains, on the eastern shore of the Lena river, acted as a strong barrier to gene flow. In comparison, postglacial recolonization of Siberia by plant species is still poorly understood. Glacial refuges and routes of recolonization have been partly documented for a few forest species: Abies sibirica Ledeb., Picea obovata Ledeb., Pinus sibirica Du Tour. and Pinus sylvestris L. All of these species were absent from the northern part of the current range during the last glaciation. In contrast, four Larix fossils from northern Siberia were radiocarbon-dated to 13, 15, 16 and 21 kya, respectively (Kremenetski, 1994), thus suggesting the presence of Larix species during the last glaciation, or at least shortly before or after its maximum. On the other hand, palaeobotanical data indicate a radical transformation of plant cover during the late Pleistocene and Holocene in Siberia (Neishtadt, 1957; Bojarskaja & Malajeva, 1967; Belova, 1985).
Introgression has played an important part in plant evolution as, for instance, in the recolonization of Europe by oaks after the last glaciation (Dumolin-Lapégue, 1998). Larch species can easily be hybridized, and morphological and genetic studies indicate that hybridization is widespread in natural populations. The natural hybrid complex between L. sibirica Ledeb. and L. gmelinii Rupr. occurs along a long belt running from the Taymir Peninsula down to Mongolia (Kruklis & Milyutin, 1977). Another possible case of larch introgressive hybridization can be found in the Russian far east (Dylis, 1961; Bobrov, 1972), where hybrid origins are suspected for many of the larch species described previously (L. amurensis Dylis, L. maritima Sukacz., L. komarovii Kolesn., L. middendorfii Kolesn., L. ochotensis Kolesn. and L. lubarskii Sukacz.). Because morphological differences between larch taxa are small, the taxonomy of Larix cannot be resolved on the basis of morphological traits only and there is currently no common view on the taxonomic status of some geographical varieties, and hence on the total number of species (Bobrov, 1972). Unsurprisingly, there is little congruence among phylogenetic trees based on morphological traits (Sukaczew, 1924; Kolesnikov, 1946; Dylis, 1947, 1961; Lepage & Basinger, 1995) and phylogenies based on molecular data (Kisanuki et al., 1995; Qian et al., 1995). Introgression and the population genetic structure of Eurasian larch species remain poorly known (but see Semerikov & Matveev, 1995; Shigapov et al., 1998, on L. sibirica; Maier, 1992 on L. decidua Mill.).
In the present study, allozyme diversity within and between Eurasian Larix species of the section Pauciseriales will be examined. More specifically, we will (i) compare genetic distances within and between species at allozyme loci with the differentiation pattern based on morphological traits, and (ii) relate genetic divergence within and between species to what is known of Larix fossil history and the extent of glaciation in eastern Siberia. If Larix species survived at high latitudes during the last glaciation, as suggested by the fossil record, no clear north– south cline in genetic variation is expected, and northern populations might contain as much genetic variation as southern ones. Furthermore, large glaciated areas seem to have existed (in the Verchojansky and Putorana mountains in the Taymir Peninsula) and acted as barriers to longitudinal migration (Fedorov et al., 1998).
Materials and methods
The species investigated, the number of populations per species and their locations are given in Table 1 and Fig. 1. The identification of the species was based on characteristic morphological traits and previously published information on the species distributions (Kolesnikov, 1946; Bobrov, 1972; Dylis, 1947, 1961). For most populations, needle tissues were analysed. When only seeds were available, seedlings were grown under artificial light for 2–4 weeks and the whole seedling used for analysis. In some populations, seeds from half-sib families were available, whereas in others the number of trees from which the seeds originated was unknown.
Allozymes
Eleven enzyme systems representing 15 protein loci were examined: glutamic-oxaloacetic transaminase (GOT, EC 2.6.1.1, two loci), glutamate dehydrogenase (GDH, EC 1.4.1.2), isocitrate dehydrogenase (IDH, EC 1.1.1.42), diaphorase (DIA, EC 1.8.1.4), phosphoglucoisomerase (PGI, EC 5.3.1.9, two loci), 6-phosphogluconate dehydrogenase (6PGD, EC 1.1.1.44), shikimate dehydrogenase (SKDH, EC 1.1.1.25), glucose-6-phosphate dehydrogenase (G6PDH, EC 1.1.1.49), phosphoglucomutase (PGM, EC 5.4.2.2, two loci), superoxide dismutase (SOD, EC 1.15.1.1, two loci) and fluorescent esterase (EST, EC 3.1.1.1).
Protein extraction and polyacrylamide gel electrophoresis were carried out according to Shurkhal et al. (1992) and Semerikov & Matveev (1995). 150 mg of needles and 150 mg of insoluble PVP were ground with liquid nitrogen in a mortar, and mixed with extraction buffer. The extraction buffer was composed of 1 M sucrose, 5.7 mM L-ascorbic acid, 8.3 mM DL-cysteine, 0.02 M dithiothreitol, and 1.5 mM aminocapronic acid dissolved in electrode buffer diluted 1:1.7. 1 mL Tween-80 was added to 100 mL of this solution and after 2–14 h at 4°C the mixture of extraction buffer, ground needles and PVP was filtered through nylon filters. The supernatant was centrifuged after adding a small amount of CCl4 in order to improve the centrifugation of the supernatant from small tissue particles. Extracts from seedlings were prepared by grinding needles in 0.15 mL extraction buffer and centrifugation with CCl4. Electrophoresis in 7.0% polyacrylamide gel in the Tris-EDTA-borate system was conducted. The electrode buffer (pH 8.0) was the following: 116 mM Tris, 3.5 mM EDTA, 161 mM boric acid. The gel buffer (pH 8.6) was the following: 118 mM Tris, 3.5 mM EDTA, 118 mM boric acid. Histochemical staining was carried out using standard methods (Harris & Hopkinson, 1976).
Data analysis
Data were analysed using GENEPOP (version 2; Raymond & Rousset, 1995b), FSTAT (version 1.2; Goudet, 1995), PHYLIP (Felsenstein, 1993) and NTSYS-pc (Rohlf, 1988). One individual per half-sib family was used in analyses of Hardy–Weinberg and linkage disequilibrium. All individuals were used in cluster analysis and ordination (similar results were obtained when one individual per half-sib family was used).
Hardy–Weinberg expectations
The fit of genotypic distributions to Hardy–Weinberg expectations was tested by the exact test proposed by Haldane (1954). The overall significance for each locus was estimated by Fisher's combined probability test (Fisher, 1954). According to this test, if P-values are obtained for each locus separately under the null hypothesis, then
is distributed according to a χ2 distribution with n d.f., where n is the number of loci (Sokal & Rohlf, 1995). FIS values, where FIS is the correlation between two uniting gametes within a subpopulation, were estimated according to Weir & Cockerham (1984). Heterozygote deficiencies or excesses were tested using an exact test (Rousset & Raymond, 1995).
Linkage disequilibrium
For each population, the nonrandom association between pairs of loci, or linkage disequilibrium, was tested using Fisher's exact test on contingency tables. Contingency tables are created for all pairs of loci in each population and an unbiased estimate of the exact probability obtained using a Markov Chain Monte Carlo method (Raymond & Rousset, 1995a). Each test is unaffected by a potential departure from Hardy–Weinberg, because each contingency table considers the genotypic composition, not the allelic composition. For each pair of loci, a global measure was obtained by averaging across populations and a global test was obtained through Fisher's combined test.
Population differentiation
Genetic differentiation between populations or groups of populations was tested for each locus separately using Fisher's exact test on contingency tables. As for linkage disequilibrium, a Markov Chain Monte Carlo method permits the attainment of an unbiased estimate of the exact probability (Raymond & Rousset, 1995a). Wright's F-statistics, FIS, FIT and FST, were estimated according to Weir & Cockerham (1984) and a 95% confidence interval was estimated by bootstrapping over loci. FIS and FIT are the correlations between two uniting gametes relative to the subpopulation and relative to the total population, respectively, and FST is the correlation between two gametes drawn at random from each subpopulation and measures the degree of genetic differentiation of subpopulations (Nei, 1987). Only statistically independent loci which did not depart from Hardy–Weinberg proportions were retained. Finally, genetic distances between populations were estimated with Nei's genetic distances (PHYLIP; Felsenstein, 1993). Ordination of populations was performed by Multidimensional Scaling (procedure MDSCALE in NTSYS-pc). A Mantel test (Sokal & Rohlf, 1995) was used to test for a relationship between genetic and geographical distances. The distance between populations was estimated using simple trigonometric formulae under the assumption that the Earth is a sphere with radius 6360 km.
Migration and isolation by distance
The number of effective migrants (Nm) between populations was estimated using the relationship Nm=(1/FST−1)/4 that holds for an island model (Wright, 1969) and by the private alleles method (Slatkin, 1985). Finally, isolation by distance was analysed according to Slatkin (1993).
Results
Polymorphism
Mean heterozygosity (Hexp) ranged from 8.8% (L. olgensis; Olga Bay) to 19.5% (Table 2; the allozyme frequency table is available from the first author). There were very few rare alleles, the only case being found in L. decidua, where Skdh-82 and Skdh-71 were present with average frequencies of 60% and 2%, respectively, whereas they were absent in other populations. It should be noted that in seedling tissues the products of these alleles displayed much lower activity (‘seminull alleles’) than other alleles. A null-allele could only be recorded in the homozygous state, as only diploid individuals were examined. Such a homozygote was observed only once in the progeny of a tree from population 11 (Yada); the seedling extract showed no activity for the analysed Skdh locus and was considered as a null-allele homozygote. Additional analysis of megagametophytes from this parent tree confirmed that the parental tree was a null-allele heterozygote.
Statistical independence among loci
Overall, 1199 pairs of loci were analysed. Among them only 50 pairs showed a statistically significant linkage disequilibrium at the 5% level. Thus, globally, independence between loci cannot be rejected.
Hardy–Weinberg equilibrium
A total of 255 exact tests of Hardy–Weinberg proportions were carried out with the Markov Chain Monte Carlo method for the 13 polymorphic loci and 32 populations. Statistically significant departures (P<0.05) were detected in 19 cases: the loci Got-B, Pgm-B and Sod-A were involved in one case each, G6pdh, 6Pgd and Skdh in two cases each, Pgm-A in three cases, and Est in seven cases. Most (18 out of 19) of these departures were heterozygote deficits. Overall, only the Est locus significantly departed from Hardy–Weinberg equilibrium (χ262=126.6, P<0.0001 for Fisher's method). Among the populations only the Magadan population (L. ochotensis) and the population from Chadan (L. sibirica) had significant deviations over all loci.
Within-species differentiation
The levels of interpopulation differentiation estimated as FST are given in Table 3. The number of effective migrants Nm, estimated according to the equation Nm=(1/FST−1)/4, was 2.91 for L. sibirica and 11.65 for L. gmelinii. The same parameter computed using Slatkin's ‘private alleles’ method was 5.57 for L. sibirica and 9.87 for L. gmelinii. There was no clear isolation by distance among L. sibirica populations (log(M)=−0.185 log(distance)+1.068, where M=(1/FST−1)/4).
The L. sibirica populations split into two groups (Fig. 2). The western group consists of the populations from the Urals and western Siberia. The eastern one includes populations located in eastern Siberia and in the north of western Siberia. The frequencies of alleles Got-A-112, Dia-85, 6Pgd-109, Sod-119 and Est-108 give the strongest differentiation between these two groups. Some of the L. sibirica populations (populations 4, 10, 7, 31, 32), characterized by intermediate frequencies of these alleles, occur between these two groups on the scatter-plot. There is a clear dependence between Nei's genetic distance between any two populations and their longitudinal distance (P=0.0094, Mantel test), but no significant relationship with latitude (P=0.6482, Mantel test; Fig. 2). Moreover, the higher the latitude, the more narrow the zone of intermediate populations. For instance, in the northern part of the Siberian larch range, population 3 which is classified as ‘western’ and population 11 which is classified as ‘eastern’ are only 260 km apart. In contrast, in southern Siberia, the allele frequencies of populations ranging from the Altai Mountains up to Tuva (populations 7, 31, 30) are intermediate between those characterizing western and eastern populations.
Between-species differentiation
Larix amurensis, L. cajanderi, L. olgensis and L. gmelinii are very similar genetically (Table 4; Figs. 2 and 3). The frequency of allele Pgi-B-79 in population 29 (Magadan) of L. ochotensis was about 50%, whereas it was less than 13% in all other populations. As illustrated by the UPGMA dendrogram (Fig. 3) and ordination of the populations using multidimensional scaling (Fig. 2), Nei's genetic distances and classical taxonomy give congruent groupings. The genetic difference between L. ochotensis (Magadan population) on the one hand, and L. olgensis, L. gmelinii and L. cajanderi on the other hand, is mainly caused by frequency differences of allele Pgi-B-79. Finally, genetic differences were more pronounced between L. kamtschatica and L. kaempferi (Table 4), than between L. kamtschatica and L. cajanderi, L. ochotensis and L. amurensis.
Discussion
Three main conclusions can be drawn from our study. First, the classification of species according to Nei's genetic distances is congruent with classical taxonomy. For instance, allozymes reveal the difference between eastern and western (L. sukaczevii Dyl.) populations of L. sibirica. Secondly, whereas there is a marked dependence of the genetic distance between any two populations of L. sibirica and the longitudinal distance separating them, there is no significant relationship with latitude. Thirdly, there is a decrease in genetic variation in marginal populations.
Genetic variation at the margin
As in Douglas-fir (Li & Adams, 1989) and Jeffrey pine (Furnier & Adams, 1986), marginal populations of larch were genetically depauperate. For example, the populations at Olga Bay and at the Polar tree line (populations 11 and 19) had low heterozygosity. These populations are generally very small and isolated and the narrowness of the ecological niche occupied by larch at Olga Bay probably results from competition from the Manchurian broad-leaved flora and historical events (climatic fluctuation, anthropogenic factors).
Population differentiation
The FST values obtained in our study were of the same order of magnitude as those observed in other conifer species (Picea abies Karst., 5.2% (Lagercrantz & Ryman, 1990); Picea maritima (Mill.) B.S.P., 5.9% (Yeh et al., 1986); Pinus sibirica du Tour, 1.6% (Krutovsky et al., 1989); Pinus sylvestris L., 2.5% (Semerikov et al., 1993)). In L. occidentalis, estimates vary from 8.6% (Fins & Seeb, 1986) to 5% (Cheliak et al., 1988; overall range). The majority of Larix species of section Pauciseriales exhibit little genetic differentiation, as do other conifers which occupy continuous and recently recolonized areas. Substantial gene flow and history (colonization from a single refuge) are possible causes of the observed lack of genetic differentiation. Based on morphological variability, Dylis (1947) concluded that L. sibirica includes two species, a western one, which he named L. sukaczewii, and an eastern one, L. sibirica, with a border lying approximately along the Ob and Irtish Rivers (see Fig. 1). Allozymes reveal that there are two genetically differentiated groups in the western and eastern parts of the L. sibirica range, but allozyme data do not correspond exactly to Dylis's hypothesis because in the southern part of the range, allele frequencies change gradually and eastern and western populations cannot clearly be distinguished. This population differentiation pattern may be the result of the recolonization of northern Eurasia by L. sibirica during the late Pleistocene and Holocene. There are strong suggestions that forests remained only in restricted refuges during the last glacial maxima. Some of the refugia were preserved in southern Siberia, south of the Urals and in northern Kazakhstan. So, at the end of glaciation, forest species extended into the north of the Urals and northern Siberia, preserving the allozyme allele frequencies of the refugia. The genetic similarity of larches from the Yamal and Taz Peninsulas (populations 8, 9 and 11) with those from the Upper Yenisey (populations 30 and 31) could, for instance, be connected with the spread of the forests along the Yenisey and Taz rivers. Fossil records also suggest that, unlike other tree species, larch glacial refugia might never have completely disappeared from higher latitudes during the last glaciation (Kremenetski, 1994), and the present genetic pattern might reflect postglacial expansion from these surviving populations. However, the results reported here are more consistent with the latter because a more pronounced longitudinal, rather than latitudinal gradient was observed. Furthermore, had the southern populations been glacial refugia one might expect to find a markedly higher number of alleles in these populations, which was not the case. On the other hand, differentiation between southern and northern populations is not pronounced and neither of the two alternatives can be unambiguously rejected. This issue may be clarified by the use of DNA markers, for instance cpDNA or mtDNA sequences, that give information on allele genealogies.
Hybridization and species differentiation
In the Russian far east only two larch species, L. gmelinii and L. olgensis, can be easily identified. Both of them have characteristic morphological traits and occupy well-defined areas. Outside these areas, morphologically variable populations occur, in which features of L. gmelinii, L. olgensis and even L. kaempferi can be found. Characteristic features of mature cones (flat, marginally jagged and nonpubescent cone scales) of some far-eastern larch species (L. cajanderi, L. ochotensis, L. amurensis and L. kamtschatica) are close to those of L. gmelinii, but other traits (large cones, numerous and sometimes convex or marginally recurved cone scales) resemble more closely those found in either L. olgensis or L. kaempferi. Therefore a hypothesis of introgressive hybridization of larches in the far east is quite plausible (Dylis, 1961; Bobrov, 1972). The introgression zone between L. gmelinii and L. sibirica probably dates from the Holocene or late Pleistocene, at least to the north of Lake Baikal. This is supported by the narrowness of this zone and the co-occurrence in this region of the marginal areas of many plant species, e.g. Chosenia arbutifolia (Pall.) A. Skvorts., Pinus pumila R.G.L., Betula middendorfii Trautv. et Mey, B. exilis Sukacz, Populus suaveolens Fish.; Sokolov et al. (1977). In contrast, the introgressive zone between L. gmelinii and L. olgensis is wide, and possibly originated in the early Pleistocene when L. gmelinii first occurred (Bobrov, 1972).
Larix olgensis appears to be genetically more similar to eastern populations of L. sibirica (D=0.020) than to L. gmelinii (D=0.028), suggesting that the L. olgensis – L. sibirica divergence took place more recently than that of L. olgensis and L. gmelinii, although the difference is not very large. However, taking into consideration the similarity of L. olgensis and L. sibirica for the main taxonomic morphological traits (Bobrov, 1972), this inference may be regarded as quite plausible. The oldest larch fossil records are from the early Pleistocene in north-eastern Siberia (Dylis, 1961). It is generally supposed that after the emergence of L. gmelinii, it expanded to the south and south-west and forced out both L. sibirica and L. olgensis, as they were less adapted to the new, more continental, climatic conditions of the Pleistocene (Dylis, 1961). As a result, the common L. sibirica – L. olgensis range was disrupted. The genetic closeness of these two species might therefore be a consequence of the relative short period of time that elapsed since this event, although additional investigations are needed to clarify this.
The larger genetic differences between L. kamtschatica and L. kaempferi (Table 4) compared to those between L. kamtschatica and L. cajanderi, L. ochotensis and L. amurensis suggest that L. kamtschatica originated from continental populations rather than from L. kaempferi. A clearer picture emerges from this allozyme study than from the cpDNA study carried out by Qian et al. (1995). These authors suggest that it might be a consequence of the importance of hybridization, cpDNA, which is paternally inherited in larch, crossing barrier species more readily than nuclear DNA. Whether this is actually the case remains to be shown.
References
Belova, V. A. (1985). Vegetation and Climate of South of East Siberia in Later Cenozoic. Nauka, Novosibirsk (in Russian).
Bennett, K. D. (1997). Evolution and Ecology. The Pace of Life. Cambridge University Press, Cambridge.
Bobrov, E. G. (1972). History and Systematics of Larix. Nauka, Leningrad (in Russian).
Bojarskaja, T. D. and Malajeva, E. M. (1967). Development of Vegetation in Siberia and Far East in Quaternary. Nauka, Moscow (in Russian).
Cheliak, W. M., Wang, J. and Pitel, J. A. (1988). Population structure and genetic diversity in tamarack L. laricina (Du Roi) K. Koch. Can J Forest Res, 18: 1318–1324.
Dumolin-Lapégue, S. (1998). Variabilité chloroplastique et mitochondriale chez les chênes blancs européens à deux échelles spatiales. Ph.D. Thesis, Université de Paris-Sud.
Dylis, N. W. (1947). Siberian Larch. Moscow Society of Nature Investigators, Moscow (in Russian).
Dylis, N. W. (1961). Larch of Eastern Siberia and Far East. Academy of Sciences of USSR, Moscow (in Russian).
Fedorov, V., Fredga, L. and Jarrell, G. H. (1998). Mitochondrial DNA variation and the evolutionary history of chromosome races of collared lemmings (Dicrostonyx) in the Eurasian Arctic. J Evol Biol. in press.
Felsenstein, J. (1993). PHYLIP (Phylogeny Inference Package). Version 3.5c. Distributed by the author, Department of Genetics, University of Washington, Seattle.
Fins, L. and Seeb, L. W. (1986). Genetic variation in allozymes of western larch (Larix occidentalis). Can J Forest Res, 16: 1013–1018.
Fisher, R. A. (1954). Statistical Methods for Research Workers. 12th edn. Oliver & Boyd, Edinburgh.
Furnier, G. R. and Adams, W. T. (1986). Geographic patterns of allozyme variation in Jeffrey pine. Am J Bot, 73: 1009–1015.
Goudet, J. (1995). FSTAT, version 1.2. A computer program to calculate F-statistics. J Hered, 86: 485–486.
Haldane, J. B. S. (1954). An exact test for randomness of mating. J Genet, 52: 631–635.
Harris, H. and Hopkinson, D. A. (1976). Handbook of Enzyme Electrophoresis. North-Holland, Amsterdam.
Kisanuki, H., Kurahashi, A., Kato, H., Terauchi, R., Kawano, S., Ide, Y. and Watanabe, S. (1995). Interspecific relationship of genus Larix inferred from the RFLPs of chloroplast DNA. J Jap For Soc, 77: 83–85.
Kolesnikov, B. P. (1946). Contribution to the taxonomy and the historical development of larches (section Pauciseriales Patschke). In: Materials on the History of the Flora and Vegetation of the USSR vol. 2, pp. 321–364. Nauka, Moscow (in Russian).
Kremenetski, C. V. (1994). Holocene history of ranges of main coniferous trees in Siberia. In: Short Periodical and Sharp Landscape and Climatic Fluctuations During Last 15 000 Years pp. 160–210. Institute of Geography of Russian Academy of Science, Moscow (in Russian).
Kruklis, M. V. and Milyutin, L. I. (1977). Larix czekanovskii. Nauka, Moscow (in Russian).
Krutovsky, K. V., Politov, D. V., Altukhov, YR. P., Milutin, L. I., Kuznetsova, G. V., Iroshnikov, A. I. et al (1989). Genetic variability in Siberian Cedar pine. IV. Genetic diversity and differentiation between populations. Genetika, 25: 2009–2032 (in Russian).
Lagercrantz, U. and Ryman, N. (1990). Genetic structure of Norway spruce (Picea abies): concordance of morphological and allozymic variation. Evolution, 44: 38–53.
Lepage, B. A. and Basinger, J. F. (1995). The evolutionary history of the genus Larix (Pinaceae). In: Ecology and Management of Larix Forests: A Look Ahead, pp. 19–29. Proceedings of an International Symposium. October, 5–9, 1992, Whitefish, MT, U.S.A.
Li, P. and Adams, W. T. (1989). Ring-wide patterns of allozyme variation in Douglas-fir (Pseudotsuga menziesii). Can J Forest Res, 19: 149–161.
Maier, J. (1992). Genetic variability in European larch (Larix decidua Mill.). Ann Sci For, 49: 39–47.
Nei, M. (1987). Molecular Evolutionary Genetics. Columbia University Press, New York.
Neishtadt, M. I. (1957). History of Forests and Paleogeography of USSR in Holocene. Academy of Sciences of USSR, Moscow (in Russian).
Qian, T., Ennos, R. A. and Helgason, T. (1995). Genetic relationships among larch species based on analysis of restriction fragment variation for chloroplast DNA. Can J Forest Res, 25: 1197–1202.
Raymond, M. and Rousset, F. (1995a). An exact test for population differentiation. Evolution, 49: 1283–1286.
Raymond, M. and Rousset, F. (1995b). GENEPOP (version 1.2). A population genetics software for exact tests and ecumenism. J Hered, 86: 248–249.
Rohlf, F. J. (1988). NTSYS-PC. Numerical taxonomy and multivariate analysis system. Exeter Publishing, New York.
Rousset, F. and Raymond, M. (1995). Testing heterozygote excess and deficiency. Genetics, 140: 1413–1419.
Semerikov, V. L. and Matveev, A. V. (1995). Investigation of genetic variation of allozyme loci in Siberian larch. Russian J Genet, 31: 944–949.
Semerikov, V. L., Podogas, A. V. and Shurkhal, A. V. (1993). Variability of allozyme loci in populations of common pine. Russian J Ecol, 24: 14–20. (in Russian).
Shigapov, Z. KH., Putenikhin, V. P., Shigapova, A. I. and Urazbakhtina, K. A. (1998). Genetic structure of the Ural populations of Larix sukaczewii Dyl. Genetika, 34: 65–74. (in Russian).
Shurkhal, A. V., Podogas, A. V. and Zhivotovsky, L. A. (1992). Allozyme differentiation in the genus Pinus. Silvae Genet, 41: 105–109.
Slatkin, M. (1985). Rare alleles as indicators of gene flow. Evolution, 39: 53–65.
Slatkin, M. (1993). Isolation by distance in equilibrium and non-equilibrium populations. Evolution, 47: 264–279.
Sokal, R. R. and Rohlf, F. J. (1995). Biometry. 3rd edn. W.H. Freeman and Co., New York.
Sokolov, S. J., Svazeva, O. A. and Kubli, V. A. (1977). Areas of Trees and Shrubs in USSR. vol. 1. Nauka, Leningrad (in Russian).
Sukaczew, V. N. (1924). On the Evolutionary History of Larix. People's Commissariat of Agriculture, Leningrad (in Russian).
Weir, B. S. and Cockerham, C. C. (1984). Estimating F-statistics for the analysis of population structure. Evolution, 38: 1358–1370.
Wright, S. (1969). Evolution and the Genetics of Populations, vol. 2, The Theory of Gene Frequencies University of Chicago Press, Chicago.
Yeh, F. C., Khalil, M. A. K., El-Kassaby, Y. A. and Trust, D. C. (1986). Allozyme variation in Picea maritima from Newfoundland: genetic diversity, population structure and analysis of differentiation. Can J Forest Res, 16: 713–720.
Acknowledgements
Financial support for this research was provided, partly, by the Russian Fundamental Research Foundation, grant no. 94–04–1282-a and Foundation ‘Bioraznoobrasie’. We thank Dr S. Shiyatov, Dr S. Zalesov, Dr V. Tarakanov, L.I. Milutin and Dr D. Malokvasov for seed collections. We are also grateful to Prof. N.V. Glotov, Dr E.A. Gileva and Dr V. Fedorov for useful discussions while preparing the manuscript, Dr Y. I. Manko for help in collection of larch species from the far east and for useful discussion of some of the results, and A. Podogas and I. Sarapultsev for collaboration in the field. Finally, we thank two anonymous referees for helpful comments on the manuscript.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Semerikov, V., Semerikov, L. & Lascoux, M. Intra- and interspecific allozyme variability in Eurasian Larix Mill. species. Heredity 82, 193–204 (1999). https://doi.org/10.1038/sj.hdy.6884710
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/sj.hdy.6884710
Keywords
This article is cited by
-
High gene flow and complex treeline dynamics of Larix Mill. stands on the Taymyr Peninsula (north-central Siberia) revealed by nuclear microsatellites
Tree Genetics & Genomes (2018)
-
Cytogenetic analysis of reciprocal hybrids and their parents between Larix leptolepis and Larix gmelinii: implications for identifying hybrids
Tree Genetics & Genomes (2010)
-
Phylogeography of Larix sukaczewii Dyl. and Larix sibirica L. inferred from nucleotide variation of nuclear genes
Tree Genetics & Genomes (2008)
-
Enhanced hybridization rates in a Larix gmelinii var. japonica × L. kaempferi interspecific seed orchard with a single maternal clone revealed by cytoplasmic DNA markers
Tree Genetics & Genomes (2008)
-
Effect of fluoride pollution on genetic variability ofLarix gmelinii(Pinaceae) in East Siberia
Journal of Forest Research (2007)