Journal of Biochemical and Microbial Toxicology
Open Access

Liriodendron tulipifera is widely distributed in Eastern America, from southern Canada to central Florida peninsula, and from Atlanta coast to the Mississippi River basin. Its sister-species, L. chinense , is regarded as an endangered species. In contrast, L. tulipifera is a pioneer and dominant tree species with strong adaption in forest community in Southeast US. Genetic diversity is the basis and premise of biodiversity, hence, one specie with low genetic diversity would not evolve or adapt to the changing environments. Furthermore, genetic diversity is a decisive factor of successful response of one specie to human disturbance, and the degree of intra specific genetic variation also has key effect on its potentiality of evolution. Since genetic diversity can be inspected with various DNA markers, this study analyzes the level of genetic diversity of L. tulipifera germplasm by using SSR markers.

Plant material

All germplasm plants (Table 1) were conserved at Xiashu forest farm affiliated to Nanjing Forestry University. In the summer of 2009~2010, 2~3 fresh leaves of each plant were collected randomly and put into the ice box rapidly after they were stored in a valve bag. The fresh leaves were brought back and kept in a refrigerator under the condition of -70°C prior to DNA extraction.

Table 1: Forty genotypes (germplasm) from five provenances of L. tulipifera were taken as plant materials

Extraction and purification of total DNA

In this experiment, the improved CTAB lysis-silica beads assay was used, in order to extract the genomic DNA of L. tulipifera [1].

SSR-PCR amplification reaction

10 μl system of SSR-PCR amplification reaction (Table 2) and touch-down PCR amplification reaction process (Table 3) were used in this experiment.

Table 2: System of SSR-PCR amplification. Note: dNTP and magnesium 2+ to oscillation blender before use.

Table 3: Procedure of PCR amplification.

Polyacrylamide gel electrophoresis of PCR products

All polyacrylamide gel electrophoresis have been carried out according to Almasizadehyaghuti et al.

Primer selection

Screening criteria of polymorphism primers followed in this experiment are as follows: amplification products should have distinct banding pattern on electropherogram of 8% polyacrylamide denaturing gel electrophoresis, various allelic variations should exist among samples, and the eliminated loci of primer should only have 1~2 bp differences. After preliminary screening, the primers which have no amplification products or undesirable polymorphisms of amplification products should be eliminated, moreover, the primers which have clear amplification bands and obvious polymorphism were selected. Finally, 22 SSR primers with obvious banding pattern, good repeatability and high stability were screened and used for this experiment.

Data management and statistical analysis

Band interpretation: SSR is co-dominant marker; the band with consistent electrophoretic mobility in the amplification product of the same primer has homology. The band data should be numbered as increasing sequences and the length of band with A, B, C, D, E…. for instance, the lowest is A, and the top most is D or E or others, the result of band reading is AB, BC, etc., besides, if one sample only has one band, the result of band reading is AA, BB, and so on.

Data analysis

Analysis of genetic diversity: The genetic parameter of populations can be obtained by using POPGENE 1.32 [2], including average number of alleles (A), effective number of allele (Ne), Kimura and Crow observed heterozygosity (Ho), expected heterozygosity (He), [3,4], Shannon diversity index (I), [5], Nei diversity index [6,7] and genetic distance (D).

Analysis of genetic relationship: Clustering analysis is conducted by utilizing software NTSYS-pc 2.10_e and the clustering figure should be established with the method of UMPGA [8].

Analysis of genetic diversity on different SSR loci

40 plants of L. tulipifera , 112 alleles are tested totally through 26 pairs SSR primers. The number of alleles (A) on each locus changes from two to six, the average value is 4.3077 (Table 2).

The effective number of allele (Ne) changes from 1.4369 to 4.2781, the average value is 2.9714. The maximum of Shannon diversity index (I) is 1.5162, the minimum is 0.5850 and the average value is 1.1842. The variation range of observed heterozygosity (Ho) is 0.0750~0.8718, the variation range of expected heterozygosity (He) is 0.3079~0.7759, with an average expected heterozygosity of 0.6541. Nei expected heterozygosity (Nei) ranges between 0.3041~0.7662, indicating that L. tulipifera maintains high genetic variation (Figure 1 and Table 1).

Figure 1: Amplification results of primer LT446.

Fixation index (F) measures the degree of deviation between actual proportion of genotype and expected proportion of Hardy-Weinberg theory (Table 4). F<0 or F>0 means that heterozygosis is excessive or homozygosis is excessive respectively. When random mating of individuals in a population cannot be achieved, its actual proportion of genotype will deviate to some extent from Hardy-Weinberg equilibrium inevitably. As for LT102, LT297, LT488, and LT921 on 26 loci, their F<0, indicating that the heterozygosis are excessive on these four loci, and as for the rest loci, their F>0, showing insufficient heterozygote on the rest loci.

Table 4: Genetic diversity of 26 SSR loci within 5 provenances of L. tulipifera. Note: An average number of alleles; Ne Effective number of alleles; I Shannon’s Information index; Ho observed heterozygosity; He average expected heterozygosity. F Fixation index; Nei Nei’s expected heterozygosity.

Genetic variation among five provenances of L. tulipifera

From Table 5, it can be found that, among five provenances, the effective number of alleles (Ne) of provenance LYS is the highest (2.4031), while provenance ZZY is the lowest (1.9382). The variation range of Nei diversity index ranges from 0.4330 (ZZY) to 0.5643 (LYS); the maximum of Shannon diversity index (I) is 0.9532 (LYS), and the minimum is 0.7009 (ZZY). Generally, LYS has the highest level of genetic diversity, while ZZY the lowest genetic diversity. As for the fixation index, the degree of deviation of NK from equilibrium is largest, although it has maximal fixation index (F=0.4187), however, the degree of deviation of ZZY from equilibrium is smallest, although it has minimal fixation index (F=0.0623).

Table 5: The genetic diversity among 5 provenances of L. tulipifera. Note: An average number of alleles; Ne Effective number of alleles; I Shannon’s Information index; Ho observed heterozygosity; He average expected heterozygosity. F Fixation index; Nei Nei’s expected heterozygosity.

Genetic structure among five provenances of L. tulipifera

Fis and Fit represent the degree of deviation of single population and total population from Hardy-Weinberg equilibrium. As for five L. tulipifera provenances, Fis=0.1965, the average Fit=0.3789, which shows that heterozygote is insufficient not only in total level but also in individuals of each provenance. Among five L. tulipifera provenances, gene flow Nm=0.8513 and genetic differentiation coefficient Fst=0.2270, which means that 22.70% of genetic variation exist among provenances, and 77.30% of genetic variation exist within provenance.

Genetic relationship of five L. tulipifera provenances

From Table 6, it can be found that the variation range of genetic identity (I) of these five L. tulipifera provenances is 0.5333~0.8135, and the variation range of genetic distance (D) is 0.2064~0.6286, with an average genetic distance of 0.4574, thus it can be deduced that these five L. tulipifera provenances have great genetic variation, and genetic relationship is relative remote. Among which, genetic identity (I) of MSL and BK is highest, which is 0.8135, showing that genetic relationship of these two populations is much closer than other ones. However, genetic identity (I) of NK and ZZY is the lowest (0.5333), indicating that genetic relationship of these two populations is much remoter than other ones.

Table 6: Nei's genetic identity and genetic distance in 5 provenances of L. tulipifera. Nei's genetic identity (above diagonal) and genetic distance (below diagonal).

According to genetic distance and genetic similarity calculated by Nei method, the clustering figure of five L. tulipifera provenances can be established with the method of UMPGA. As shown in Figure 2, when the threshold value is 0.51, these five L. tulipifera provenances can be divided into two categories, among which, ZZY and LYS belong to one category and the rest three provenances belong to the other category. Meanwhile, when the threshold value is 0.47, these three provenances can be divided in two subgroups further: namely, BK and MSL belong to one category, and NK constitutes one category alone.

Figure 2: Dendrogram of L. tulipifera in 5 provenances.

Genetic diversity of five L. tulipifera provenances

Genetic diversity is the result of long-term evolution of the species and represents the potential of evolution of species, which is the premise for the survival and development of populations. In order to survive and develop in an unstable environment, the species are obliged to generate and accumulate genetic variation. On the basis of the results of this research, the effective number of allele of five L. tulipifera provenances is 2.9714, the variation range of Nei expected heterozygosity (Nei) is 0.3041~0.7662, which show that L. tulipifera has high level of genetic variation. It can further declare that L.tulipifera is one pioneer tree species with high adaption and wide distribution range, moreover, continues geographical distribution increases the pollen flow and seed dispersal among populations. Meanwhile, the large gene flow among populations also results in high level of allelic richness and genetic diversity of L.tulipifera.

We would like to thank all professors’ Huogen Li and Qiang Zhuge lab crews and students that helped us to carry out experiments and prepare manuscript.