Bot. Bull. Acad. Sin. (2003) 44: 13-18

Fu et al. — RAPD analysis for genetic diversity in Changium smyrnioides

RAPD analysis for genetic diversity in Changium smyrnioides (Apiaceae), an endangered plant1

Chengxin Fu*, Yingxiong Qiu, and Hanghui Kong

Lab of Plant Systematic Evolution and Biodiversity, College of Life Sciences, Zhejiang University, Hangzhou 310029, P.R. China

(Received January 31, 2002; Accepted September 2, 2002)

Abstract. Random Amplified Polymorphic DNA (RAPD) markers were used to measure genetic diversity of Changium smyrnioides Wolff (Apiaceae), a genus endemic to eastern China and an endangered medicinal plant, collected from five populations along the Yangzi River. A total of 92 amplified bands were scored from the 13 RAPD primers, and a mean of 7.1 amplified bands per primer and 69% (64 bands) percentages of polymorphic bands (PPB) was found. The Shannon's index was used to partition genetic diversity. Genetic diversity estimates indicated that 51.2% of total diversity was among populations and 48.8% within populations. The species shows higher genetic diversity between populations than other endangered plants. In this study, the RAPD results of genetic diversity are similar to those of previous studies employing allozymes (58% among populations and 42% within populations). The RAPDs provide a useful tool for assessing genetic diversity of rare, endemic species and for resolving relationships among populations. The results show that the genetic diversity of this species is high, possibly allowing it to more easily adapt to environmental variations. A pattern of greater variation among rather than within populations was observed with both markers. This distributive pattern of genetic variation of C. smyrnioides populations provides important baseline data for conservation and collection strategies for this species. It is suggested that most of the populations should be sampled in ex-situation protection so as to retain as much genetic diversity as possible.

Keywords: Allozymes; Changium smyrnioides; Genetic diversity; RAPD.

Introduction

Changium smyrnioides Wolff, a species of monotypic genus in the family Apiaceae, is an endangered herbaceous perennial plant in China. It is endemic to eastern China, and is distributed only in fragments in Zhejiang, Anhui, Jiangsu, and Hubei Provinces (see Figure 1). The species is not only of important medicinal value; it is also of interest to plant systematists. The drugs made from this species promote the production of body fluid to quench thirst, nourish "yin" to moisten the lung, soothe the throat, and clear away toxic material that may cause skin infection. As there might be some internal reasons, such as questions of reproductive biology and population ecology, this species has been assigned an endangered status and is listed in the Red Data Book as an endangered plant by the government of China (Zhang, 1991).

The species was named C. smyrnioides by Wolff (1924) based on a specimen from Wuxin, Zhejiang. Various aspects of the genus have been discussed in earlier papers (Qiu et al., 2000; Qiu and Fu, 2001), but it should be emphasized that natural populations of C. smyrnioides are very small and rare. Because allozymes have sometimes

proven insufficiently variable for assessing genetic diversity within and among populations of rare endemic plants, attention has turned to more variable regions of the genome. A more recently employed approach in plant systematics and population biology is random amplified polymorphic DNA (RAPD) markers, a PCR-based (polymerase chain reaction) technique.

1State Key Basic Research and Development Plan of China, G2000046806, and project of NNSF of Zhejiang, China, No. 395064.

*Corresponding author. Tel: +86-571-86971576; Fax: 86-571-86971634; E-mail: fucxsmi@mail.hz.zj.cn

Figure 1. The distribution area and sampling locations of Changium smyrnioides.


Botanical Bulletin of Academia Sinica, Vol. 44, 2003

Although earlier studies (Pan et al., 1983, 1995; Zhang et al., 1983; Li et al., 1987; Sheh and Su, 1987; Su and Sheh, 1990; Cheng et al., 1993) have reported some results of research on chemical, cytological, pollen morphological, ecological, geographical, and karyotypic characters of the species; neither the genetic diversity nor the divergence were clear. The purpose of the present study is to assess genetic diversity and divergence within and among populations of the species using RAPD markers and to compare the results with those obtained from allozymes. Another important aim is to provide genetic data and a theoretical basis for protection of the species.

Materials and Methods

Population Sampling

Total DNA was extracted from plants collected in the field. Leaves were dried in plastic bags with silica gel until extracted in the laboratory at Zhejiang University, Hangzhou, Zhejiang, China. A total of 35 individuals from five populations (Ningbo, Langya, Nanjing, Anqing, Hangzhou) were included in the study. The distance between plants collected was at least 10 m to increase the possibility of detecting the variation potential within each population. Details on material and vouchers are given in Table 1.

DNA Extraction and Amplification

The CTAB (hexadecyltri-methylammonium bromide) method of Dellaporta et al. (1983) was used to extract total DNA. Leaf material was powdered in liquid nitrogen, mixed with 2 ml extraction buffer [1.4 M NaCl, 100 mM Tris-HCl (pH 8.0), 20 mM EDTA, 2% CTAB, and 1% 2-mercaptoethanol] at 65°C, and incubated at 65°C for 45 min with slow shaking every 10 min. Proteins were extracted twice with 2 ml chloroform: iso-amylalcohol (24:1), then centrifuged at 13,000 rpm for 2 min. RNase (10 µg/ml) was added to the supernatants for 30 min at 37°C, and the mix

ture was centrifuged at 13,000 rpm for 2 min again. The sediment was washed twice in 70% ethanol, vacuum-dried, and resuspended in 100 µl Tris-EDTA buffer (1 mM Tris-HCl and 0.1 mM EDTA, pH 8.0).

Amplification of genomic DNA was made on a Perkin Elmer DNA Cycler (BIOMETRA, Germany), using the arbitrary decamers. The 13 primers were selected from the "S" RAPD primer kit purchased from Sangon (Shanghai, China) (Table 2); these primers included S14, S15, S16, S21, S23, S24, S31, S34, S35, S38, S39, S44 and S45. Amplifications of genomic DNA were performed in 25-µl reaction volumes containing 1.2 units of Taq polymerase (Sangon, Shanghai, China), 10 mM Tris-HCl (pH 9.0), 25 mM KCl, 2 mM MgCl2, 0.2 mM of each dNTP, 24 ng each of random primer and 40 ng of template DNA. The cycle program included an initial 75 sec denaturation at 94°C, followed by 45 cycles of 15 sec at 94°C, 30 sec at 42°C and 75 sec at 72°C, with a final extension at 72°C for 7 min. RAPD fragments were separated electrophoretically on 1.5% agarose gels in 1X TBE buffer, stained with ethidium bromide, and photographed on a UV transilluminator using a digital camera. DNA from each plant was amplified with the same primer more than once, and the banding patterns were compared.

Optimization of RAPD Protocol

Because the RAPD-PCR technology is sensitive to changes in experimental parameters, a total of 50 primers were initially screened against ten plants selected from all populations. The effects of magnesium, template DNA concentrations, pH values, and length of the denaturation stage of the amplification were all examined. When trying to optimize annealing temperatures, we ran the test reactions at 37°C, 39°C and 42°C. The decamer primers can be clearly amplified at 42°C. A subset of 13 primers for further analysis were based on the following criteria: (i) consistent, strong amplification products, and (ii) production of uniform, reproducible fragments between replicate PCRs.


Fu et al. — RAPD analysis for genetic diversity in Changium smyrnioides

Data Analysis

Fragment sizes were designated as amplified bands, and bands were shared as diallelic characters (present = 1, absent = 0). Those bands amplifying in each instance were scored and included in the analyses. The number of multi-locus genotypes (unique arrays of amplified bands) was calculated for each population. RAPD locus diversity was calculated with the Shannon-Weaver information statistic, employing the Brilliouin formula for eliminating the bias of finite sample size (Yun et al., 1998; Peet, 1974; Whitkus et al., 1998). The formula is Ho=-SPilog2 (Pi) (Pi for frequency of the i band, Ho for diversity), which estimates the genetic diversity within and among populations of the species.

Results

A total of 92 bands were scored for the 13 RAPD primers for a range from 4 to 10, corresponding to an average of 7.1 bands per primer, and 69% (64 bands) of these were

polymorphic. Percentages of polymorphic bands (PPB) for each primer ranged from 43% to 89%. PPB within populations for each primer ranged from 0.10 of S23 in NJ to 0.78 of S44 in NB and LY, which can be seen in Table 3. Mean RAPD similarities for individuals of different populations ranged from 0.22 to 0.89, with an average of 0.31 for all pairwise comparisons. Populations HZ and LY exhibit the highest level of variability (PPB=0.43), and population NJ exhibits the lowest (PPB=0.28).

Shannon's index of phenotypic diversity was used to class diversity into components within and among populations. Table 4 shows that the diversity within populations related to different primers and different populations. Primer S44 detected the highest genetic diversity within these populations while primer S34 detected the lowest. As for different populations, the populations HZ and LY showed greater average variation (1.246 and 1.215) than the other populations while the population NJ on average showed less variation (0.465) than the others.


Botanical Bulletin of Academia Sinica, Vol. 44, 2003

found that the two species of Dendroseris were with the highest allozyme diversity and also exhibited the highest RAPD diversity, and one species showed low diversity with both markers.

In the studies, the level and distribution of genetic diversity detected by RAPDs are in agreement with previous results by allozymes (Qiu et al., 2000). The allozyme study found significantly lower gene diversity within populations of NJ and AQ compared with the other three populations. Likewise, the results from this RAPD analysis show that the NJ and AQ populations are of lower genetic polymorphism at the population levels. The RAPD genetic diversity within and among populations was also similar to the results from allozyme analysis, which showed about 51.2% among and 48.8% within populations in the former (Table 5) and 57.8% among and 42.2% within population in the latter. This is only slightly lower than the diversity seen with allozymes. Considering the data from both markers, we can conclude that genetic diversity among populations is higher than within populations of Changium smyrnioides. It was found from field investigations that the distributive pattern of the species is becoming island-like as a result of human activity and forest destruction. Results from RAPD and allozymes indicate that genetic drift might have occurred among the studied populations, thereby producing population differentiation. Genetic distance among populations of Changium smyrnioides at allozyme loci ranged from 0.386 for HZ and AQ to 0.006 for NJ and LY. It was found from field observations that the NJ population was a small, artificial group that was planted there from the Zijing Mountain near the LY population decades ago. This might explain why the genetic identity between them is high (0.994).

The percentage of polymorphic bands (69%) of RAPD in the species was higher than in other endangered plants, e.g. Lactoris fernandeziana (Lactoridaceae) 24.5% (Brauner et al., 1992), Cathaya argyrophylla 32% (Wang et al., 1996), Paeonia suffruticosa 22.5% and P. rockii 27.6% (Pei et al., 1995), and Dacydium pierrei 33.3% (Su et al., 1999). This shows that the species' genetic diversity is not low, and it should be able to fit the environmental variation. Accordingly, it is a main reason that the human's activities to damage their habitats and dig them excessively for medicine make the species' population decrease in size and their habitat be island-like in distribution. The populations HZ and LY, with higher genetic diversity, are located in the West Lake Mountains and in the Forest Park of Langya Mountain, respectively, and both sites enjoy good environmental protection. With a larger area of population, the probability of crossing among the individuals increases, which results in the retention of genetic

Usually, Shannon's index of mean phenotypic diversity can be divided into two types: within populations and among populations (Table 5). In brief, a little more than half of the total variation (51.2%) is presented among the populations, and 48.8% of the variation is presented within the populations. Of course different primers had different percentages. For example, the variation among populations detected by primer S24 was only 18.8% whereas primer S34 detected 73.5%.

Discussion

Allozymes' value is sometimes limited for rare endangered plants because of their limited genetic variation although they have been widely and effectively used to analyze genetic diversity within and among plant populations (Crawford et al., 1994; Hamrick and Godt, 1996). RAPD markers have been employed as an alternative (Brauner et al., 1992; Dawson et al., 1993; Pei et al., 1995; Su et al., 1999; Wolfe and Liston, 1998; Yun et al., 1998; Esselman et al., 2000). Data from RAPD analysis has indicated that their diversity is usually similar to or greater than diversity from allozymes in plant species (Esselman et al., 1999, 2000). Comparing allozyme and RAPD diversity in rare endangered plants, the latter often prove more variable than the former. For example, no allozyme variation was found in some species of Dendroseris although RAPD diversity was detected (Esselman et al., 2000).

Both the Shannon-Weaver index (King and Schaal, 1989) and the Nei index (Nei, 1972) can be used to estimate the genetic diversity of populations, but the data analyzed by the Nei index need strict dominant and recessive allelic frequency. The data from RAPDs do not depend on these criteria. Estimating Gst using RAPDs has been problematical due to their dominance, and analytical methods usually rely on knowledge of the selfing rate or assume Hardy-Weinberg equilibrium. This assumption does not hold when populations exhibit fixed heterozygosity, so an alternative method, the Shannon-Weaver index, was used to partition the genetic diversity. However, this index is, to some extent, in accord with the analysis of RAPD data. Wei et al. (1999) used the two indices to analyze RAPD data from Caragana and found that both were consistent in showing the variation of populations. Do allozymes and RAPD markers provide similar estimates of the genetic diversity for rare endangered plants? As discussed above, we know that RAPD markers are inherited as dominants although the banding patterns cannot be used to analyze gene diversity. This is in contrast to the analysis of allozyme data, in which band frequencies can be directly interpreted as allelic frequencies. Esselman et al. (2000)


Fu et al. — RAPD analysis for genetic diversity in Changium smyrnioides

variation. Though many individuals of the species have been recorded in the Lao Mountains of Jiangsu Province, they have disappeared gradually along with environmental changes in their habitat.

Only with a firm grasp of the genetic structure and the diversity grade of populations in rare and endangered species can we make an efficacious measurement and strategy of protecting them. For a species with limited gene flow and over 50% variation among populations, it is necessary to collect samples from at least six populations in order to conserve 95% of the genetic diversity of the species. For a species with only 20% variation among populations, the samples taken from two populations are enough to get the same results above (Hamrick et al., 1991; Pei et al., 1995; Yun et al., 1998). The 51.2% genetic differentiation coefficient of Changium smyrnioides from RAPD (Table 5) suggests that the species is of a higher genetic diversity among populations than other endangered plants. The damage to their habitats is a main reason they are so rare. It is, therefore, a good strategy to protect more of their habitats. The results of RAPD and isozyme also show that we need to take individuals from more different populations if we are to construct an artificial conservation area so as to preserve their diversity for the future.

Acknowledgements. The authors thank the Institute of Botany, Jiangsu, Naijing for permission to observe specimens; Mr. Tai Zhe-Ming, Jiangxi for doing portions of the field work; and Dr. Kameron, New York Botanical Garden for valuable comments on the manuscript.

Literature Cited

Brauner, S., D.J. Crawford, and T.F. Stuessy. 1992. Ribosomal DNA and RAPD variation in the rare plant family Lactoridaceae. Amer. J. Bot. 79: 1436-1439.

Cheng, Q., Z.Y. Hang, and S.X. Zong. 1993. Ecogeograpgical distribution, utilization and conservation on Changium smyrnioides. J. Chin. Materia Med. 18: 327-329.

Crawford, D.J., D.W. Haines, M.B. Cosner, D. Wiens, and P. Lopez. 1994. Lactoris fernandeziana on the Juan Fernandez Islands: allozyme uniformity and field observations. Conserv. Biol. 8: 277-280.

Dawson, K., K.J. Chalmers, and R. Waugh. 1993. Detection and analysis of genetic variation in Hordeum spontaneum populations from Israwl using RAPD markers. Mol. Ecol. 2: 151-159.

Dellaporta, S.L., J. Wood, and J.B. Hicks. 1983. A plant DNA minipreparation (version II). Plant Mol. Biol. Rep. 1: 19-22.

Esselman, E.J., D.J. Crawford, S. Brauner, T.F. Stuessy, G.J. Anderson, and M.O. Silva. 2000. RAPD marker diversity within and divergence among species of Dendroseris (Asteraceae: Lactuceae). Amer. J. Bot. 87: 591-596.

Esselman, E.J., L. Jianqiang, D.J. Crawford, J.L. Windus, and A.D. Wolfe. 1999. Clonal diversity in the rare Calamagrostis porteri ssp. insperata (Poaceae): comparative results for allozymes and random amplified polymorphic DNA and intersimple sequence repeat markers. Mol. Ecol. 8: 443-453.

Hamrick, J.L., H.J.W. Godt, and K.A. Murawski. 1991. Correlations between species traits and allozyme diversity implications for conservation biology. In D. A. Falk and K. E. Holsinger (eds.), Genetics and Conservation of Rare Plants. Oxford Univ. Press, New York, Oxford, pp. 75-86.

Hamrick, J.L. and H.J. W. Godt. 1996. Conservation genetics of endemic plant species. In J. C. Avise and J. L. Hamrick (eds.), Conservation Genetics: Case Studies from Nature. Chapman and Hall, New York, New York, USA, pp. 281-304.

King, L.M. and B.A. Schaal. 1989. Ribosomal DNA variation and distribution in Rudbckia missouriensis. Evolution 43:1117-1119.

Nei, M. 1972. Genetic distance between populations. Amer. Nat. 106: 283-292.

Li, X., J.W. Cheng, and J. Sun. 1987. Isolation, identification and quantitative analysis of water soluble active principle in the roots of Changium smyrnioides. Nat. Prod. Research Develop. 7: 1-5.

Pan, Z.H., Z.J. Wu, and X.T. Liu. 1983. Karyotype analysis of Changium smyronioides. In S. A. He (ed.), Bulletin of Nanjing Botanical Garden. Jiangsu Sci. and Tech. Press, Nanjing, pp. 33-35.

Pan, Z.H., M.L. She, X.T. Liu, and X.M. Yao. 1995. On karyotypes and geographical distribution of endemic genera in Umbelliferae from China. J. Plant Resources Environ. 4: 1-8.

Pei, Y.L., Y.P. Zou, Z. Yin, X.Q. Wang, Z.X. Zhang, and D.Y. Hong. 1995. Preliminary report of RAPD analysis in Paeonia suffruticosa subsp. spontanea and P. rockii. Acta Phytotaxon. Sin. 33: 350-356.

Peet, R.K. 1974. The measurement of species diversity. Annu. Rev. Ecol. Sys. 5: 285-307.

Qiu, Y.X., A.J. Huang, and C.X. Fu. 2000. Studies on genetic diversity in Changium smyrnioides Wolff (Umbelliferae). Acta Phytotaxon. Sin. 38: 111-120.

Qiu, Y.X. and C.X. Fu. 2001. Studies on the endangerment mechanism of and conservation strategies for Changium smyrnioides. Conserv. Biol. (Beijing) 9: 151-156.

Sheh, M.L. and P. Su. 1987. The floristic analysis of endemic genera in Chinese Umbellifferae. In S. A. He (ed.), Bulletin of the Nanjing Botanical Garden Mem. Sun Yat Sen. Jiangsu Sci. and Tech. Press, Nanjing, pp. 14-23.

Su, P. and M.L. Sheh. 1990. Study on pollen morphology of endemic genera in Chinese Umbellifferae. In S. A. He (ed.), Bulletin of the Nanjing Botanical Garden Mem. Sun Yat Sen Jiangsu Sci. and Tech. Press, Nanjing, pp. 33-36.

Su, Y.J., T. Wang, and C. Huang. 1999. RAPD analysis of different population of Dacydium pierrei. Acta Sci. Nat. Univ. Sunyatseni 38: 99-101.

Yun, R., M. Zhong, H.X. Wang, W. Wei, Z.A. Hu, and Y.Q. Qian. 1998. Study on DNA diversity of Liaodong population at Dongling mountain region. Acta Bot. Sin. 40: 169-175.

Wang, X.Q., Y.P. Zou, D.M. Zhang, and D.Y. Hong. 1996. RAPD analysis for genetic polymorphism in Cathaya argyrophylla. Science in China (C) 26: 437-441.

Wei, W., H.X. Wang, Z.A. Hu, M. Zhong, R. Yun, and Y.Q. Qian. 1999. Primary studies on molecular ecology of Caragana spp. populations distribution over Mao Wusu Sandy grassland from RAPD data. Acta Ecol. Sin. 19: 17-


Botanical Bulletin of Academia Sinica, Vol. 44, 2003

21.

Whitkus, R., M.D.L. Cruz, L. Mota-Bravo, and A. Gomez-Pompa. 1998. Genetic diversity and relationships of cacao (Theobroma cacao L.) in southern Mexico. Theor. Appl. Genet. 96: 621-627.

Wolfe, A.D. and A. Liston. 1998. Contributions of PCR-based methods to plant systematics and evolutionary biology. In D. E. Soltis, P. S. Soltis, and J. J. Doyle (eds.), Molecular Systematics of Plants II, DNA Sequencing. Kluwer Aca

demic Publishers, New York, New York, USA, pp. 43-86.

Wolff, H. 1924. In Fedde Repert. Sp. 19: 315.

Zhang, W.P. (eds.). 1991. Conservation and Study of Endangered Plant. Chinese Environment Science Press, Beijing, pp. 187-188.

Zhang, Z.H., Z.J. Wu, and X.T. Liu. 1983. Karyotypic analysis of Changium Smyrnioides. In S. A. He (ed.), Bulletin of the Nanjing Botanical Garden Mem. Sun Yat Sen. Jiangsu Sci. and Tech. Press, Nanjing, pp. 33-35.