0737-B2

The population Structure of Gmelin Larch in The Central Evenkia.

A. Ya. Larionova, N. V. Yakhneva 1


Abstract

On the basis of 17 genes, coding allozyme diversity MDH, IDH, G-6-PD, ME, 6-PGD, GDH, LAP, SKDH, EST, AAT, the values of principal parameters of genetic variation were determinated in three stands of Gmelin larch (Larix gmelinii (Rupr.) Rupr) occurring near the settlement Tura (Evenkia). It was established that 50.98% of gene loci assayed were polymorphic at 95% and 78.43% at 99% polymorphism criterions.

The mean number of alleles per locus, the mean observed and expected heterozygosities, effective number of alleles per locus were 2.000, 0.0915, 0.1246 and 1.179, respectively. More than 98% of total genetic variation resides within stand and only 1.66% among stands. The mean value of Nei's genetic distance (D) between stands of L. gmelinii was 0.0035. The data obtained point to a very low differentiation in the Evenkia Gmelin larch stands.


Intruduction

Larix gmelinii is the main forest-forming species in the Central Evenkia. The most characteristic feature of Gmelin larch is its location in permafrost region (Pozdnyakov, 1975). Gmelin larch grows in the northern regions of Krasnoyarsk Territory and Irkutsk district, including the Putoran Mountains, Nizhnyaya Tunguska and the higher parts of Podkamennay Tunguska River basins. (Milyutin, Vishnevetskaia, 1995). In the Evenkia this species is dominating. The share of Gmelin larch in of forest cover territory is 81% (Abaimov et al., 2000). Genetic diversity and differentiation of Gmelin larch forests in the Evenkia up to the present time remains pooly studied. The purpose of our study was to investigate gene diversity, population structure and differentiation of Gmelin larch from Central Evenkia by starch-gel electrophoresis.

Material and Methods.

Seed material for electrophoretic analysis was collected from 76 trees in three stands of Gmelin larch in Central Evenkia near the settlement Tura. The stand I is situated at right bank of Kochechum River (64˚ 19' N, 100˚ 07' E). The stand II is located at left bank of Kochechum River (64˚ 19' N, 100˚ 13' E). The stand III is placed in the of Nizhnyaya Tunguska River valley (64˚ 17' N, 100˚ 14' E).

Individual trees were genotyped using 8 to 50 megagametophytes. The megagametophytes were sampled randomly from each tree a set extracted from minimum of 20 cones. The gomogenates were electrophoresed in horizontal chamber on 13-14 % starch gels. For electrophoresis of enzymes two buffer systems were used: A) tris-citrate, pH 6.2 (Adams, Joly, 1980), B) tris-citrate, pH 8.5 / LiOH-borate, pH 8.1 (Ridgway et al., 1970). Enzyme extraction, electrophoresis, histochemical staining followed standart methods (Brewer, 1970; Vallejos, 1983; Goncharenko et al., 1989) with insignificant modifications. The enzymes assayed, their abbreviations, the buffer systems used are given in Table 1.

To estimate level of genetic variation in the populations studied the values of principal parameters (P - percentage of polymorphic loci, A - number of alleles per locus, Ho - observed and He - expected heterozygosities, ne - effective number of alleles) were determinated. The percentage of polymorphic loci was estimated at 95% (the frequency of the most common allele was not greater than 95%) and 99% polymorphism criterions. The mean number of alleles per locus was calculated for all the alleles revealed. To measure gene diversity in natural stands of Gmelin larch the parameters of Wright's F-statistics were used (Guries, Ledig, 1982). The genetic distance coefficient (D) among the stands was estimated by Nei's method (Nei, 1972). The effective number of alleles was calculated, using following formula:

ne=1/(1-He).

Results and Discussion

The electrophoretic analysis revealed 40 allelic variants at 17 loci. The most variable loci in the Evenkia stands of Gmelin larch were Est-1, Mdh-3 and G-6-pd. The Mdh-4, Skdh, Aat-1, Aat-2, Aat-3, Idh, Me-2, 6-Pgd loci showed an intermediate level of variability. Five loci (Lap-1, Lap-2, Mdh-2, Me-1, Gdh) were less polymorphic and one locus (Mdh-1) was monomorphic in all stands (Table 2).

On the basis of allelic frequencies for 17 loci the principal parameters of genetic variation were computed (Table 3). The values of all parameters vary among the stands. The most high level of genetic variation was revealed in the oldest stand III, situated in the Nizhnyaya Tunguska River valley. On average in L. gmelinii stands from Evenkia 50.98% of the loci assayed were polymorphic at 95% and 78.43% at 99% polymorphism criterions, the mean number of alleles per locus was 2.000, the average observed and expected heterozygosities, effective number of alleles were 0.0915 and 0.1246, 1.179 respectively. In all stands studied of Gmelin larch the mean values of observed heterozygosity were lesser than those expected under Hardy-Weinberg conditions.

Obtained values of observed heterozygosity were lower than values of Ho established by Semerikov et al. (1999) at 15 loci in the populations of L. gmelinii from Eastern Siberia and than those determinated by Potenko, Rasumov (1996) at 21 loci in two populations from Khabarovsk Territory (Far East). The values of expected heterozygosity, mean number of alleles per locus, effective number of alleles were over the range of values established in other populations studied of this species.

Population structure of L. gmelinii in the Evenkia was analysed by means of Wright's F-statistics (Guries, Ledig, 1982). Data on the correlation between uniting gametes within stands (Fis), among stands (Fst), for the species as a whole (Fit) were calculated on the basis of analysis for 17 loci (Table 4).

The Fis values ranged from - 0.0780 at Aat-2 to +0.5218 at Est-1 and averaged 0.2658. Obtained value is positive. This indicates to a deficiency of heterozygotes within L. gmelinii stands studied relative to Hardy-Weinberg expections. As a whole in L. gmelinii from Evenkia deficiency of heterozygotes was 0.2781. Value of Fst for polymorphic loci ranged from 0.0022 (Lap-1) to 0.0267 (Aat-3). The average value of Fst for all the 17 loci assayed was equal to 0.0166 (1.66%). This means that more than 98% of the total genetic diversity resides within stands of L. gmelinii from Evenkia and only 1.66% among stands. Fst value point to low extent genetic differentiation between the stands studied.

Genetic similarity among L. gmelinii stands in Central Evenkia was quantified using Nei's genetic distance coefficient (Nei, 1972). The values D obtained for three stands are listed in Table 5.

It is obvious that genetic differences between stands are insignificant. The values of D range from 0.0025 to 0.0037 and averaged - 0.0035. These values D are typical for stands within one population and for very closely related geographically connected populations of conifers. The most genetic similarity (D=0.0025) was revealed between stands I, situated at right bank of Kochechum River, and III, placed in the of Nizhnyaya Tunguska River valley.

References

Abaimov A.P., Zyryanova O.A., Prokushkin S.G., Koike T., Matsuura Y., 2000. Forest Ecosystems of the Cryolithic Zone of Siberia; Regional Features, Mechanisms of Stability and Pyrogenic Changes. Eurasian J. For. Res., 1: 1 - 10.

Adams W.T., Joly R.I., 1980. Genetics of Allozyme Variants in Loblolly Pine. Heredity, 71: 33 - 40.

Brewer G.J., 1970. Introduction to isozyme techniques. N.Y.-L. Academ press: 186 p.

Goncharenko G.G., Padutov V.E., Potenko V.V., 1989. Guide to conifer species research by isozymes methods. Gomel.: 164 p. (In Russian).

Guries R.P., Ledig F.T., 1982. Genetic diversity and population structure in Pitch pine (Pinus rigida Mill.). Evolution, 36: 387 - 402.

Milyutin L.I., Vishnevetskaia K.D., 1995. Larch and larch forests of Siberia. Ecol. and Manag. of Larix Forests: A Look Ahead, Ogden, USA: 50 - 53.

Nei M., 1972. Genetic distance between populations. Amer. Natur., 106: 283 - 291.

Potenko V.V., Rasumov P.N., 1996. Genetic variability and population structure in Gmelin larch at territory of Khabarovsk region. Lesovedenie, 5: 11 - 18. (In Russian).

Pozdnyakov L.K., 1975. Larix dahurica. Moscow, Russia, Nauka: 312 p. (In Russian).

Ridgway G.J., Sherburne S.W., Lewis R.D., 1970. Polymorphism in the Esterases of Atlantic Harring. Trans. Am. Fish. Soc., 99: 147 - 151.

Semerikov V.L, Matveev A.V., 1995. Investigation of geneticvariability in Siberian larch (Larix sibirica Ldb.) at isozyme loci. Genetika, 31: 1107 - 1113.(In Russian)

Semerikov V.L., Semerikov L.F., Lascoux M., 1999. Intra and interspecific allozyme variability in Eurasian Larix Mill. species. The Genetical Society of Great Britain. Heredity, 82: 193 - 204.

Shurchal A.V., Podogas A.V., Semerikov V.L., Zhivotovskii L.А., 1989. Allozyme polymorphism in Siberian larch (Larix sibirica Ledeb.). Genetika, 25: 1899 - 1901. (In Russian).

Vallejos G.E., 1983. Enzyme activity staining. Isozymes in plant genetics and breeding. Amsterdam: Elsevier Sci. Publ. B.V.: 469 - 516

Table 1. Enzymes, their Abbreviations (Abbr.), Enzyme Commission numbers (EC) and buffer systems used for electrophoresis.

Enzyme

Abbr.

EC

Buffer

Malate dehydrogenase

MDH

1.1.1.37.

A

Glutamate dehydrogenase

GDH

1.4.1.2.

B

Isocitrate dehydrogenase

IDH

1.1.1.42.

A

6-Phosphogluconate dehydrogenase

6-PGD

1.1.1.44.

A

Glucose-6-phosphate dehydrogenase

G-6-PD

1.1.1.49.

A

Malic enzyme

ME

1.1.1.40.

B

Leucine aminopeptidase

LAP

3.4.11.1.

B

Esterase

EST

3.1.1.1.

B

Aspartate aminotransferase

AAT

2.6.1.1.

B

Shikimate dehydrogenase

SKDH

1.1.1.25

A

Table 2. Allele frequencies for 17 loci in three stands of Gmelin larch in Evenkia

Locus

Allele

Stand

I

II

III

Gdh

1
2

0.9583
0.0417

0.9800
0.0200

1.0000
-

Idh

1
2

0.0625
0.9375

-
1.000

0.0741
0.9259

G-6-pd

1
2

0.8542
0.1458

0.9400
0.0600

0.8148
0.1852

6-Pgd

1
2
3

0.0625
0.9375
-

0.1200
0.8800
-

0.0370
0.9444
0.0185

Skdh

1
2

0.0208
0.9792

0.0600
0.9400

0.0185
0.9815

Mdh-1

1

1.0000

1.0000

1.0000

Mdh-2

1
2

0.9792
0.0208

1.0000
-

0.9815
0.0185

Mdh-3

1
2

0.1458
0.8542

0.3000
0.7000

0.2037
0.7963

Mdh-4

1
2
3
n

-
0.9583
-
0.0417

0.0400
0.9400
-
0.0200

0.0185
0.9074
0.0370
0.0370

Me-1

1
2

-
1.0000

-
1.0000

0.0370
0.9630

Me-2

1
2

0.0833
0.9167

0/0400
0.9600

0.0556
0.9444

Lap-1

1
n

0.9792
0.0208

0.9600
0.0400

0.9630
0.0370

Lap-1

1
n

1.0000
-

1.0000
-

0.9630
0.0370

Aat-1

1
2
3

-
1.0000
-

0.0400
0.9400
0.0200

0.0185
0.9259
0.0556

Aat-2

1
2

0.0833
0.9167

0.0200
0.9800

0.0741`
0.9259

Aat-3

1
2
3

0.0208
0.9375
0.0417

-
0.9400
0.0600

0.0185
0.8333
0.1482

Est-1

1
2
3
4

0.2708
-
0.6875
0.0417

0.3600
0.0600
0.5800
-

0.3704
0.0370
0.5926
-

Table 3. Genetic variation in natural stands of L. gmelinii.

Stand

P95, %

P99, %

A

Ho

He

ne

I
II
III

47.06
47.06
58.82

76.47
70.59
88.24

1.882
1.882
2.235

0.0858
0.0776
0.1111

0.1113
0.1153
0.1473

1.150
1.177
1.211

Mean

50.98

78.43

2.000

0.0915

0.1246

1.179

Table 4. Estimates of Fis, Fit, Fst- in Gmelin larch

Locus

Fis

Fit

Fst

Mdh-1
Mdh-2
Mdh-3
Mdh-4
Est-1
Skdh
G-6pd
Lap-1
Lap-2
Aat-1
Aat-2
Aat-3
Idh
6-Pgd
Me-1
Me-2
Gdh

-
-0.0201
0.1670
-0.0564
0.5218
-0.0461
0.4077
-0.0361
-0.0385
0.2582
-0.0780
0.4818
0.2538
0.2728
-0.0385
0.1519
-0.0359

-
-0.0133
0.1869
-0.0470
0.5269
-0.0343
0.4219
-0.0337
-0.0125
0.2739
-0.0629
0.4957
0.2719
0.2833
-0.0125
0.1568
-0.0210

0.0000
0.0067
0.0238
0.0089
0.0105
0.0113
0.0241
0.0022
0.0250
0.0211
0.0140
0.0267
0.0244
0.0144
0.0250
0.0057
0.0144

Average over 17 loci

0.2658

0.2781

0.0166

Table 5. Estimates of Nei's genetic distance coefficient (D), based upon data from 17 loci

Stand

I

II

III

I

-

   

II

0.0042

-

 

III

0.0025

0.0037

-



1 V.N. Sukachev Institute of Forest, SB RAS, Akademgorodok, Krasnoyarsk, Russia, 660036
Phone: +7 (3912) 433686,
Fax: +7 (391-2) 43-36-86,
e-mail: [email protected]