L. (Apiaceae) is a perennial emergent plant designated as an endangered species in South Korea. According to the former records, only four natural habitats remain in South Korea. A former study suggested that three of four populations (Pyeongchang: PC, Hoengseong: HS, Gunsan: GS) would be classified as different ecotypes based on their different morphological characteristics and life cycle under different environmental conditions. To evaluate this suggestion, we estimated genetic diversity in each population and distance among three populations by random amplification of polymorphic DNA. Seven random primers generated a total of 61 different banding positions, 36 (59%) of them were polymorphic. Nei’s gene diversity and the Shannon diversity index increased in the order of PC < HS < GS, which is the same order of population size. In the two-dimensional (2D) plot of first two principal components in principal component analysis with the presence of 61 loci, individuals could be grouped as three populations easily (proportion of variance = 0.6125). Nei’s genetic distance for the three populations showed the same tendency with the geographical distance within three populations. And it is also similar to the result of discriminant analysis with the morphological or life-cycle factors from the previous study. From the results, we concluded that three different populations of
should be classified as ecotypes based on not only morphology and phenology but genetic differences in terms of diversity and distance as well.
L. (Apiaceae) is a perennial emergent plant that has toxic substance (cicutoxin) in its roots and tubers (Panter et al. 1988). It disperses by its overwintering rootstocks or buoyant fruits (Mulligan and Munro 1981). Natural populations of
were found at north and central Europe, northern Asia, northern America (
a, Mulligan and Munro 1981, Lee 2003). South Korea is a southern marginal area of distribution, population size of
has decreased (Shin et al. 2013). Under the experimental condition with elevated CO
concentration and temperature according to the general climate change scenario, the timing of flowering, flower stalk forming and seed formation were earlier than the control condition and germination rate of self-pollinated seeds were higher at experimental condition (Hong 2012). Therefore, longterm climate change seems to be able to affect the phenology of
. Several countries (e.g., Japan, Hungary) have designated
as a protected species (Tatar 2010, Japanese Wildlife Research Association 2012), and
also has been designated as a domestic endangered species in Korea (Korea National Arboretum 2012).
Former studies about
mainly have concentrated
Distribution map of Cicuta virosa. (a)Worldwide occurrence of C. virosa populations(■ = occurrence) and (b) three C. virosa populationsfor sampling in Korea. Biodiversity occurrencedata (a) published by Global BiodiversityInformation Facility (2013).
on toxicology or pharmacology of cicutoxin (Sarreveld and Hope 1975, Strauß et al. 1996).
requires wet substrate and this requirement limits its habitat range such as streams, ditches, lakes, rivers or marshy areas (Panter et al. 1988). In South Korea, only four natural habitats of
were remained and their water and soil conditions were similar to the optimal environmental range of their accompanying species. However, only the water level was different from optimal range of others (Shin et al. 2013). Three of four populations in Korea (except the smallest and rapidly decreasing population) were suggested to be classified as different ecotypes based on their different morphological characteristics and phenology under different environmental conditions (Shin and Kim 2013).
Understanding of the genetic structure of vulnerable species is important to formulate sustainable conservation strategy (Ellstrand and Elam 1993, Lee et al. 2004, Kim et al. 2008). Several molecular analysis tools, such as DNA sequencing, random amplification of polymorphic DNA (RAPD), inter simple sequence repeat (ISSR), and amplified fragment length polymorphism (AFLP) has been utilized in ecological studies (Beebee and Graham 2008). To study the genetic difference among three populations, we used the RAPD method that has been used for investigating the genetic diversity and structure of populations and that is independent to the former knowledge of their DNA sequence (Williams et al. 1990, Newbury and Ford-Lloyd 1993, Beebee and Graham 2008). This method has been used for understanding the formation and maintenance of natural plant populations (Bonnin et al. 1996, Min et al. 2012), intraspecific phylogenic studies (Kim et al. 2008, Ok and Yoo 2012), and genetic difference among the varieties or ecotypes (Virk et al. 1995, Ray and Roy 2009, Arghavani et al. 2010).
The aim of this study is to assess the genetic diversity within three
populations in Korea and the differences among them using RAPD markers and to understand the relationship between the genetic difference and ecotypic differences of Shin and Kim (2013). The suggestion to classify the three
populations as ecotypes would be able to be supported with the result of this study of the genetic difference.
MATERIALS AND METHODS
- Sampling sites and populations description
Only four natural habitats of
were reported in South Korea and sampling was conducted at three of the four habitats (
b). Pyeongchang (PC) was a wet meadow that located at 793 m above sea level, and accompanying species with
. At Hoengseong (HS), two subpopulations (12-year-old abandoned paddy field and streamlet) were located with about 700 m distance.
was a major accompanying species at the abandoned paddy field and
at the streamlet. Gunsan (GS) was a reservoir which located at relatively low latitude than PC and HS.
has been lived on the floating mat with
. The major accompanying species at GS was
. Population sizes in PC, the abandoned paddy field at HS, streamlet at HS and GS were <40, 170, 70, >500, respectively (see Shin et al. 2013 for detailed description). Population in Gangneung was excluded for sampling because of its small population size (<10).
- Seed sampling and germination
At the year of 2012, after the flowering time of each three populations, seeds of the senescent
individuals were collected randomly from more than 5 individuals. Seeds were stored at 4℃ until the germination. Collected seeds were germinated in a growth chamber, on the wet filter paper in petri dishes. Each cycle of light and
Random primer sequences and number of amplified RAPD fragments*Number in the parenthesis is the number of polymorphic bands for each primer.
Random primer sequences and number of amplified RAPD fragments *Number in the parenthesis is the number of polymorphic bands for each primer.
temperature condition was kept in 12 hours interval, light (at 30℃) and dark (at 15℃) for a week.
After about a week from the germination, five seedlings were randomly collected in each of three populations. Each sample size of them was determined in previous studies about the distinction among ecotypes using RAPD markers (Milella et al. 2006, Ray and Roy 2009, Arghavani et al. 2010). Whole genomic DNA extraction was followed with DNeasy Plant Mini Kit (Qiagen, Hilden, Germany) according to manufacturer’s protocol.
PCR for RAPD analysis was conducted with AccuPower PCR premix (Bioneer, Daejeon, South Korea). Total reaction volume of PCR was 20 μL, composed of 1 μL of DNA extracted from each seedling, 1 μL of primer (10 pmol), and 18 μL of deionized water. Seven of 33 random primers (Bioneer) which showed the polymorphism in electrophoresis were selected for RAPD analysis (
). Reaction cycle of PCR began with the 2 min of denaturation at 94℃. 40 cycles were followed with 45 s of denaturation at 94℃, 45 s of annealing at the T
of each primer, and 1 min 30 s of elongation at 72℃. Lastly, 2 min of the final extension step at 72℃ was followed.
Electrophoresis of amplified fragments was conducted in 1% agarose gel stained with GelRed™ (Biotium, Heyward, CA, USA). Under the UV light, bands were observed and recorded by binary character matrix with 1 (presence) and 0 (absence).
- Data analysis
Genetic diversity within the population (number of polymorphic loci, Nei’s gene diversity: h, Shannon index: i), genetic distance (Nei’s genetic distance) among populations and unweighted pair-group method with arithmetic averages (UPGMA) dendrogram were calculated by Popgen32 with the binary matrix data (Yeh and Boyle 1997). Principal component analysis (PCA) with presence and absence of each band was conducted by R ver. 2.15.0 (R Development Core Team 2013).
RESULT AND DISCUSSION
- Genetic diversity among individuals of three populations
61 band loci were detected from seven random primers and 36 loci (59%) of them were shown to be polymorphic. 10 (28%) of 36 polymorphic loci appeared only one or two populations. DNA sequence of these specific loci could be able to design the population-specific markers (Ray and Roy 2009).
In three populations, average of Nei’s gene diversity (h) was 0.0747, and it of Shannon index (i) was 0.1115. GS was shown to be the highest gene diversity (
In GS, the population size of
in 2011 over 500 individuals was largest in three populations. Population sizes of it at abandoned paddy field and Narrow streamlet of HS were about 170 and 70, respectively. PC has the smallest size, less than 40 individuals (Shin et al. 2013). Genetic diversity of the randomly selected seeds from three
populations also seems to be affected by their population size at natural habitat (Ellstrand and Elam 1993, Frankham 1996, Fischer et al. 2000, Kim et al. 2008).
From the result of PCA with band presence and absence of 61 loci, 5 principal components had 83% of total variance and proportion of variance of major 9 principal components was 95%. On the 2D plot of major two principal components (PC1, PC2), three populations were clearly separated and grouped within each population (
). This result implies that three populations could be grouped based on the RAPD-PCR results.
Number of polymorphic loci and percentage, Nei’s gene diversity (h), and Shannon index (i) in three populations*Total does not mean the sum of polymorphic loci numbers of three populations but the total loci number which considers overlapped loci among populations.
Number of polymorphic loci and percentage, Nei’s gene diversity (h), and Shannon index (i) in three populations *Total does not mean the sum of polymorphic loci numbers of three populations but the total loci number which considers overlapped loci among populations.
2D plot of principal component analysis (PCA) of 15 seedlingswith presence and absence of 61 band loci. 15 individuals are plotted atthe first two principal components.
UPGMA dendrogram based on Nei’s genetic distance for thethree populations of Cicuta virosa. Nei’s genetic distances were PC-HS =0.1999, PC-GS = 0.2331, and HS-GS = 0.2618.
- Genetic distance among three populations
Genetic distance among three populations based on Nei’s genetic distance was the highest between HS and GS (0.2618). The shortest distance was between PC and HS (0.1999), and the distance between PC and GS was 0.2331. GS seems to be genetically apart from PC and HS in the UPGMA dendrogram from the RAPD markers (
PC and HS are located in Gangwon province, the eastern part of Korean peninsula, whereas GS is located in Jeollabuk province, adjacent to the western sea of Korea and the reservoir is located on relatively lower latitude than PC and HS. Straight-line distances between PC and HS, PC and GS, HS and GS are about 50 km, 261 km and 218 km, respectively. Genetic distance among populations could be increased similar to the increasing of the geographic distance when genetic flow among populations and genetic drift in each population make equilibrium (Hutchinson and Templeton 1999). Based on the genetic and geographic distance among three populations, GS seems they had been separated earlier from the common ancestor than the other two populations.
Ecotypic characteristics such as seed length and weight, flowering time, and sensitivity about the moisture stress were significantly different among three populations (Shin and Kim 2013). However, such as root:shoot (R:S) ratio or relative growth rate (RGR) in length in GS was significantly different from HS and PC, whereas the two populations (HS and PC) were not significantly different from each other. In addition, as a result of a 2D plot of the discriminant analysis, GS was relatively away from HS and PC (Shin and Kim 2013). This result seems to be affected from the genetic distance among three populations.
Based on the RAPD analysis, genetic diversity in each population seems to be affected from their population size and some polymorphic loci appeared only one or two populations. In a result of PCA with the band presence, three populations were separated from each other and grouped within each population on 2D plot. In addition, genetic distance among three populations was associated to their geographic distance and it also corresponds to the Shin and Kim’s study (2013) about the morphological and phenological characteristics. In conclusion, we suggest it is possible to distinguish the
individuals into three populations by the result of RAPD-PCR and this support the suggestion of classifying three populations as different ecotype.
This study was supported by the Korea Ministry of Environment as “The Eco-Innovation project, 416-111-010” and by the Center for Aquatic Ecosystem Restoration (CAER) of the Eco-STAR project from the Ministry of Environment, Republic of Korea (MOE;08-Ⅲ-12).
Genetic diversity in ecotypes of two Agropyron species using RAPD markers.
Res J Environ Sci
DOI : 10.3923/rjes.2010.50.56
An Introduction to Molecular Ecology
Oxford University Press
High level of polymorphism and spatial structure in a selfing plant species, Medicago truncatula (Leguminosae), shown using RAPD markers.
Am J Bot
DOI : 10.2307/2446261
van Kleunen M
RAPD variation among and within small and large population of the rare clonal plant Ranunculus reptans (Ranunculaceae).
Am J Bot
DOI : 10.2307/2656649
Global Biodiversity Information Facility
GBIF data portal.
Effects of elevated CO2concentration and temperature on the phenology, growth response and reproductive ecology of Cicuta virosa, endangered plant in Korea.
Correlation of pairwisegenetic and geographic distance measures: inferring therelative influences of gene flow and drift on the distributionof genetic variability.
DOI : 10.2307/2640449
Japanese Wildlife Research Association
In: Search system of Japanese RED data.
Genetic diversity and populationstructure of endangered Isoetes coreana in SouthKorea based on RAPD analysis.
DOI : 10.1016/j.aquabot.2008.02.004
Korea National Arboretum
Korea biodiversity information system.
Colored Flora of Korea.
Genetic diversityand its implications in the conservation of endangeredZostera japonica in Korea.
J Plant Biol
DOI : 10.1007/BF03030519
Relationships between an Italianstrawberry ecotype and its ancestor using RAPD markers.
Genet Resour Crop Ev
DOI : 10.1007/s10722-005-1405-7
Assessment of genetic diversityof Typha angustifolia in the development of cattailstands.
J Ecol Field Biol
DOI : 10.5141/JEFB.2012.004
The biology of Canadianweeds: 48. Cicuta maculata L., C. douglasii (DC.) Coult.& Rose and C. virosa L.
Can J Plant Sci
DOI : 10.4141/cjps81-012
The use of RAPD for assessingvariation in plants.
Plant Growth Reg
DOI : 10.1007/BF00144581
Habitats ecological characteristics ofAsplenium scolopendrium L. and its RAPD Analysis.
Korean J Plant Res
DOI : 10.7732/kjpr.2012.25.6.719
Toxicoses in livestock from the hemlocks (Conium and Cicuta spp.).
J Anim Sci
R Development Core Team
R: a language and environment for statistical computing.
R Foundation for Statistical Computing
Genetic diversity of Amaranthus speciesfrom the Indo-Gangetic plains revealed by RAPD analysisleading to the development of ecotype-specific SCARmarker.
DOI : 10.1093/jhered/esn102
Comparison of environmentalcharacteristics at Cicuta virosa habitats, an endangeredspecies in South Korea.
J Ecol Environ
DOI : 10.5141/ecoenv.2013.003
Cicutoxin from Cicuta virosa: a new and potentpotassium channel blocker in T Lymphocytes.
Biochem Biophys Res Commun
DOI : 10.1006/bbrc.1996.0233
Useof RAPD for the study of diversity within plant germplasmcollections.
DOI : 10.1038/hdy.1995.25
DNA polymorphisms amplified by arbitrary primersare useful as genetic markers.
Nucleic Acids Res
DOI : 10.1093/nar/18.22.6531
Population genetic analysis of codominantand dominant markers and quantitative traits.
Belg J Bot