DNA Polymorphism and Assessments of Genetic Relationships in genus Zoysia Based on Simple Sequence Repeat Markers
DNA Polymorphism and Assessments of Genetic Relationships in genus Zoysia Based on Simple Sequence Repeat Markers
Journal of Life Science. 2015. Mar, 25(3): 256-262
Copyright © 2015, Korean Society of Life Science
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : October 08, 2014
  • Accepted : March 09, 2015
  • Published : March 30, 2015
Export by style
Cited by
About the Authors
만규 허

The genetic variability of four species of the genus Zoysia collected from South Korea was analyzed using an inter-simple sequence repeat (ISSR) marker system. Polymerase chain reactions (PCR) with eight ISSR primers generated 86 amplicons, 76 (87.1%) of which were polymorphisms. The polymorphism information content (PIC) value of the ISSR marker system was 0.848. The percentage of polymorphic loci ( P p ) ranged from 41.2% to 44.7%. Nei’s gene diversity ( H ) ranged from 0.149 to 0.186, with an average overall value of 0.170. The mean of Shannon’s information index ( I ) value was 0.250. Total genetic diversity values ( H T ) varied between 0.356 (ISSR-1) and 0.418 (ISSR-16), for an average overall polymorphic loci of 0.345. Interlocus variation in within-species genetic diversity ( H S ) was low (0.170). On a per-locus basis, the proportion of total genetic variation due to differences among species (G ST ) was 0.601. This indicated that about 60.1% of the total variation was among species. Thus, about 39.9 of genetic variation was within species. The estimate of gene flow, based on G ST , was very low among species of the genus Zoysia ( N m = 0.332). The phylogenic tree showed three distinct groups: Z. macrostachya and Z. tenuifolia clades and other species were formed the separated clusters. In conclusion, the ISSR assay was useful for detecting genetic variation in the genus Zoysia , and its discriminatory power was comparable to that of other genotyping tools.
The genus Zoysia Willd. (Family Poaceae , subfamily Chloridoideae , tribe Zoysiece ) includes about ten species [7] and is commonly known as zoysiagrass ( Zoysia sp.) or lawn. Zoysiagrass is found throughout East Asia (China, Korea and Japan) and Australasia, grown especially in warm regions [17 , 27] .
Wild species belonging genus Zoysia have been survived though long-term natural selection [13] , and thereby can tolerate wide variations in temperature, sunlight, and water. Zoysiagrass is widely used for turfgass and forage grass in as Korea and other countries in East Asia. Zoysia grasses can diminish soil erosion on slopes, and are excellent at repelling weeds throughout the year. Some types of Zoysia are available commercially as making lawn ground, park grass, and garden grass. In addition, they are used on golf courses to create fairways and teeing areas after several hybridization processes.
DNA markers have numerous applications in plant molecular genetic research [6] . The two most common uses of DNA markers have been the assessment of genetic diversity within plant germplasm and evolutionary relationship within or among species [12 , 19] . One major use of DNA techniques is to reveal genetic diversity within and between populations or species.
Inter simple sequence repeat (ISSR) technique is a polymerase chain reaction (PCR) based technique, reported by Zietkiewicz et al. [34] . ISSR technique involves PCR amplification of DNA using a single primer composed of a microsatellite sequence. ISSR technique involves amplification of DNA segments between two identical microsatellite repeat regions oriented in opposite direction using primers designed from microsatellite core regions [29] . The ISSR has mild technical difficulty, good reproducibility and reasonable cost, permitting its use for genetic studies of population. PCR products were resolved by performing agarose gel electrophoresis. Compared to arbitrary markers such as random amplified polymorphic DNA (RAPD), ISSR markers are highly reproducible due to the use of longer primers.
Considering the potentials of the ISSR marker based genetic diversity analysis, the present study aimed to evaluate the extent of genetic diversity and phylogenetic relationships among four species of genus Zoysia collected from Korea.
Materials and Methods
- Sample materials
Four species within genus Zoysia Willd. occur in Korea. Zoysia japonica Steud., Z. macrostachya Franch. et Savat., Z. sinica Hance, and Z. tenuifolia Willd. were selected. Plant materials collected from four large natural populations per species. To avoid including individuals from the same lineage by clonal reproduction, the distance between the selected individuals was about 5.0 m. Arundinella hirta (Thunb.) Koidz. was used as an outgroup species in this study.
- Genomic DNA isolation and ISSR analysis
The total genomic DNA was extracted from the leaf tissues using the plant DNA Zol Reagent (Life Technologies Inc., Grand Island, New York, USA) according to the manufacturer’s protocol. DNA quantity was checked using a mini fluorometer (TKO 100 Mini-Fluorometer, Hoefer Scientific Instruments).
ISSR primers (University of British Columbia, Canada) synthesized by Sigma Aldrich Inc., were used for the polymorphism survey ( Table 1 ).
Lists of decamer oligonucleotide utilized as primers, their sequences, and associated bands for Korean genusZoysia
PPT Slide
Lager Image
Lists of decamer oligonucleotide utilized as primers, their sequences, and associated bands for Korean genus Zoysia
PCR (polymerase chain reaction) amplification was carried out in 25 μl reaction containing 25 ng of genomic DNA, 1 μM of primer, 10 mM of dNTPs (2.5 mM each), 50 mM KCl, 10 mM Tris HCl (pH 8.3), 2.5 mM MgCl2, 0.3 units (U) of Taq DNA polymerase (Promega, USA). PCR conditions were programmed for initial denaturation at 94℃ for 5 min, 36 cycles of 1 min denaturation at 94℃, 45 sec annealing at 42-50℃, 2 min extension at 72℃, and final extension for 10 min at 72℃. The annealing temperature for PCR amplification was maintained based on the specificity of the primer pair used. PCR amplified products of ISSR primers were subject to horizontal gel electrophoresis using 2.0% agarose gel and stained with ethidium bromide. The bands on gels were documented using Alpha Imager 1200TM (Alpha Innotech Co., USA). All the reactions were repeated twice and only reproducible bands were scored for analyses.
- Data analysis
The unambiguous ISSR bands were scored visually as present (1) or absent (0) and the binary data matrix was constructed. Several standard genetic parameters, the percentage of polymorphic loci ( P p ), mean number of alleles per locus ( A ), effective number of alleles per locus ( A E ), Nei’s [21] gene diversity ( H ), and Shannon’s Information index ( I ) were estimated using the computer program, POPGENE ver. 1.31 [33] . Polymorphism information content (PIC) value was calculated using the formula PIC, 1 - Σ pi 2 , where pi is the frequency of the i th allele [25] .
The estimation of genetic similarity (GS) between genotypes was based on the probability that an amplified fragment from one individual will also be present in another [23] . GS was converted to genetic distance (1-GS) [15] .
Species differentiation analyses have been assessed within and among according to Nei's gene diversity formulae ( H T , H S , and G ST ) [22] . The G ST coefficient corresponds to the relative amount of differentiation among populations or species. Furthermore, gene flow ( N m ) between the pairs of species was calculated from G ST values by Nm = 0.5(1/ G ST -1) [20] .
Shannon’s index of genotypic diversity ( H O ) for ISSR was estimated to quantity the degree of within species diversity following the formula [3] : H O = –Σ pi log pi , where pi is the frequency of a particular phenotype i .
H O can be calculated and compared for different populations [24] . Let
H POP = 1/n H O
be the average diversity over the n different populations and let
H SP =–Σp log p
be the diversity of species calculated from the phenotypic frequencies p in all the species considered together [24] . Then the proportion of diversity within species, H POP / H SP , can compared with that between species ( H SP - H POP )/ H SP .
- Cluster analyses
A phenetic relationship was constructed by the neighbor joining (NJ) method using the NEIGHBOR program in MEGA5 [28] .
For analysis of variability of genus Zoysia , eight primers were used for studying the ISSR banding patterns across the entire samples. The eight primers generated a total of 85 consistently scorable bands, 76 of which were polymorphic (89.4% polymorphism) ( Table 1 ). The average number of amplification products was 10.6 per primer; the maximum was 15 with (GT) 8 YT (ISSR-12), whereas the minimum was 8 with (GA) 8 CG (ISSR-5).
Polymorphism information content (PIC) for ISSR primers ranged from 0.791 to 0.896 with an average of 0.848 per primer ( Table 1 ).
In a simple measure of inter-populations variability i.e. the percentage of polymorphic bands ( P p ), Z. japonica exhibited the highest variation (47.1%) and Z. tenuifolia the lowest (41.2%) ( Table 2 ). The average number of alleles per locus ( A ) was 1.439 across species, varying from 1.412 to 1.471. The effective numbers of alleles per locus ( A E ) at the lowest species and the highest species level were 1.253 and 1.334, respectively. The mean genetic diversity within species was 0.170. Overall, Z. sinica and Z. japonica exhibited high variation among species. Z. macrostachya and Z. tenuifolia was shown the low genetic variation.
Codes of countries and measurements of genetic variation for genusZoysia
PPT Slide
Lager Image
Codes of countries and measurements of genetic variation for genus Zoysia
Total genetic diversity values ( H T ) varied between 0.356 (ISSR-1) and 0.418 (ISSR-16), for an average over all polymorphic loci of 0.345 ( Table 3 ). Interlocus variation in the within-species genetic diversity ( H S ) was low (0.170). On a per-locus basis, the proportion of total genetic variation due to differences among species ( G ST ) ranged from 0.273 for ISSR-1 to 0.660 for ISSR-6, with a mean of 0.506. This indicated indicated that about 50.6% of the total variation was among species. Thus, about genetic variation (40.4%) resided within species. The estimate of gene flow, based on G ST , was very low among species of genus Zoysia ( N m = 0.488). Values of genetic distance (D) were <0.470 ( Table 4 ). Genetic identity values among pairs of populations ranged from 0.625 to 0.813. Genetic distances between species were high. There was not shown significant difference among four species.
Estimates of genetic diversity of genusZoysiain Korea
PPT Slide
Lager Image
Total genetic diversity (HT), genetic diversity within populations (HS), the proportion of total genetic diversity partitioned among populations (GST), and gene flow (Nm).
Genetic identity (above diagonal) and genetic distances (below diagonal) of genusZoysiabased on ISSR
PPT Slide
Lager Image
Genetic identity (above diagonal) and genetic distances (below diagonal) of genus Zoysia based on ISSR
An assessment of the proportion of diversity present within species, H POP / H SP , indicated that about 71.0% the total genetic diversity was among species. Thus, the other portion of genetic variation (29.0%) resided within species ( Table 5 ).
Partitioning of the genetic diversity into within and among genusZoysiain Korea
PPT Slide
Lager Image
Partitioning of the genetic diversity into within and among genus Zoysia in Korea
Clustering of four species of genus Zoysia , using the NJ algorithm, was performed based on the matrix of calculated distances ( Fig. 1 ). Four species of genus Zoysia were well separated each other. The phylogenic tree showed three distinct groups; Z. macrostachya and Z. tenuifolia clade and the other species.
PPT Slide
Lager Image
The tree for Korean genus Zoysia based on ISSR analysis using MEGA5.
Genetic diversity of Korean Zoysia species is comparable with other species, although the use of different methods (e.g., isozyme [co-dominant marker] and RAPD [dominant marker], ISSR [sequence based marker system) may preclude meaningful direct comparisons.
In ISSR analysis, nine species belonging to genus Zoysia maintain moderate or higher than mean level of genetic diversity compared with other plant species [10] . For example, the percentage of polymorphic loci at the species level for Zoysia is 43.9%, which is similar to species with temperature- zone distributions (48.5%), species with a sexual and asexual reproduction mode (43.8%), and those with a short-lived perennial herbaceous (41.3%) (Hamrick and Godt 1989). Its genetic diversity of 0.170 is higher than that for temperate-zone species (0.146), species with a sexual and asexual reproduction mode (0.138), and those with a short-lived perennial herbaceous (0.116) [10] .
The results of other research for Zoysia were similar to the results of this study. Choi et al. [5] analyzed 93 native zoysiagrass lines collected from the southern and western coastal regions of South Korea. They also estimated genetic diversity of zoysiagrass using RAPD [4] . Weng et al. [30] also analyzed genetic diversity of 131 clones of Taiwan Zoysia species by using isozyme and RAPD. Li and Tong [16] reported genetic diversity of 105 plants from seven Z. sinica populations in different regions of China using ten RAPD primers. About 4.8, 30, and 70% of the genetic variation was detected among groups, populations, and within populations, respectively. In ISSR literature, polymorphism and genetic diversity ( H ) within species in Chinese zoysiagrass were 96.7% and 0.250, respectively [31] .
The relatively high level of genetic variation found in genus Zoysia is consistent with several aspects of its biology. First, the particular breeding system that a plant possesses is an important determinant of variability at both the species and population levels. Flowering plants exhibit a spectacular diversity in reproductive systems, and this can have important effects on the amount and structuring of genetic variation within and among populations [9] . Species of the genus Zoysia are rarely reproduced sexually and pertetuates theyselves clonally by rhizomes [1] . Whereas sexual reproduction might initially act to enhance genetic variation, asexual reproduction can maintain this enhanced variability [2] . Second, a perennial and/or long-lived species generally maintains relatively higher levels of variation than do annuals [18] . Because individuals of genus Zoysia are so longlived, opportunities for the accumulation of mutations should be high [14] .
Genetic differentiation among populations is principally a function of natural selection, genetic drift, and gene flow among populations via pollen and seed dispersal [26] . Of the total variation observed in Zoysia species about 20.0-27.7% was due to differences among species ( G ST = 0.601, H VAR )/ H SP = 0.710). Predominantly wind-pollinated outcrossing species have on an average less than 10 % of the genetic variation between populations [11] . This high level of genetic differentiation also suggests that gene flow among species is low ( Nm = 0.332), though no chromosomal barrier of interspecific crossing [8 , 32] .
Anderson S. J. 2000 Taxonomy of Zoysia (Poaceae): Morphological and molecular variation. Ph.D. dissertation Texas A&M Univ. Texas, USA
Bayer R. J. 1990 Patterns of clonal diversity in the Antennaria rosea (Asteraceae) polyploid agamic complex Am. J. Bot. 77 1313 - 1319    DOI : 10.2307/2444591
Bowman K. D. , Hutcheson K. , Odum E. P. , Shenton L. R. 1971 Comments on the distribution of indices of diversity Stat. Ecol. 3 315 - 359
Choi J. S. , Ahn B. J. , Yang G. M. 1997 Classification of zoysiagrasses (Zoysia spp.) native to the southwest coastal regions of Korea using RAPDs J. Kor. Soc. Hort. Sci. 38 789 - 795
Choi J. S. , Ahn B. J. , Yang G. M. 1997 Distribution of native zoysiagrasses (Zoysia spp.) in the south and west coastal regions of Korea and classification using morphological characteristics J. Kor. Soc. Hort. Sci. 38 399 - 407
Collard B. C. Y. , Mackill D. J. 2009 Start codon targeted (SCoT) polymorphism: a simple, novel DNA marker technique for generating gene-targeted markers in plants Plant Mol. Biol. Rep. 27 86 - 93    DOI : 10.1007/s11105-008-0060-5
Engelke M. , Anderson S. , Casler M. D. , Cuncan R. R. 2003 Turfgrass Biology, Genetics, and Breeding John Wiley & Sons Hoboken, NJ Zoysiagrasses (Zoysia spp.) 271 - 286
Forbes I. 1952 Chromosome numbers and hybrids in Zoysia Agron. J. 44 194 - 199    DOI : 10.2134/agronj1952.00021962004400040008x
Glemin S. , Bazin E. , Charlesworth D. 2006 Impact of mating systems on patterns of sequence polymorphism in flowering plants Proc. Biol. Sci. 273 3011 - 3019    DOI : 10.1098/rspb.2006.3657
Hamrick J. L. , Godt M. J. W. , Brown A. H. D. , Clegg M. T. , Kahler A. L. , Weir B. S. 1989 Plant Population Genetics, Breeding and Genetic Resources Sinauer Sunderland, MA Allozyme diversity in plant species 304 - 319
Hamrick J. L. , Godt M. J. W. , Sherman-Broyles S. L. 1992 Factors influencing levels of genetic diversity in woody plant species New For. 6 95 - 124    DOI : 10.1007/BF00120641
Hu J. , Vick B. A. 2003 Target region amplification polymorphism: a novel marker technique for plant genotyping Plant Mol. Biol. Rep. 21 289 - 294    DOI : 10.1007/BF02772804
Jin H. , Han L. B. 2004 Progress on genetic diversity of Zoysia japonica Steud J. Bejing. For. Uni. 26 91 - 95
Ledig F. T. , Soule M. E. 1986 Conservation Biology Sinauer Sunderland, MA Heterozygosity, heterosis, and fitness in outbreeding plants 77 - 104
Le Thierry d’Enneequin M. , Poupance B. , Starr A. 2000 Assessment of genetic relationships between Setaria italica and its wild relative S. viridis using AFLP markers Theor. Appl. Genet. 100 1061 - 1066    DOI : 10.1007/s001220051387
Li Y. , Tong H. Y. 2004 Genetic differentiation in Zoysia sinica populations revealed by RAPD markers Guihaia 24 345 - 349
Loch D. S. , Simon B. K. , Poulter R. E. 2005 Taxonomy, distribution, and ecology of Zoysia macrantha Desv., an Australian native species with turf breeding potential Int. Turfgrass Soc. Res. J. 10 593 - 599
Loveless M. D. , Hamrick J. L. 1984 Ecological determinants of genetic structure in plant populations Ann. Rev. Ecol. Syst. 15 65 - 95    DOI : 10.1146/
Lubberstedt T. , Melchinger A. E. , Duble C. , Vuylsteke M. , Kuiper M. 2000 Relationships among early Europe maize inbreds: IV. Genetic diversity revealed with AFLP markers and comparison with RFLP, RAPD, and pedigree data Crop Sci. 40 783 - 791    DOI : 10.2135/cropsci2000.403783x
McDermott J. M. , McDonald B. A. 1993 Gene flow in plant pathosystems Ann. Rev. Phytopathy 31 353 - 373    DOI : 10.1146/
Nei M. 1973 Analysis of gene diversity in subdivided populations Proc. Natl. Acad. Sci. USA 73 3321 - 3323
Nei M. 1977 F-statistics and analysis of gene diversity in subdivided populations Ann. Human Genet. 41 225 - 233    DOI : 10.1111/j.1469-1809.1977.tb01918.x
Nei M. , Li W. H. 1979 Mathematical model for studying genetical variation in terms of restriction endonucleases Proc. Natl. Acad. Sci. USA 74 5267 - 5273
Paul S. P. , Wachira F. N. , Powell W. , Waugh R. 1997 Diversity and genetic differentiation among populations of Indian and Kenyan tea (Camellia sinensis (L.) O. Kuntze) revealed by AFLP markers Theor. Appl. Genet. 94 255 - 263    DOI : 10.1007/s001220050408
Shete S. , Tiwari H. , Elston R. C. 2000 On estimating the heterozygosity and polymorphism information content value Theor. Pop. Biol. 57 265 - 271    DOI : 10.1006/tpbi.2000.1452
Slatkin M. 1985 Rare alleles as indicators of gene flow Evolution 39 53 - 65    DOI : 10.2307/2408516
Stewart A. 2005 The potential for domestication and seed propagation of native New Zealand grasses for turf Royal New Zealand Institute of Horticulture Conference 2003: Greening the City-Bringing Biodiversity Back into the Urban Environment Christchurch, New Zealand October 21-24 277 - 284
Tamura K. , Peterson D. , Peterson N. , Stecher G. , Nei M. , Kumar S. 2011 MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods Mol. Biol. Evol. 28 2731 - 2739    DOI : 10.1093/molbev/msr121
Vijayan K. 2005 Inter simple sequence repeat (ISSR) Ppolymorphism and its application in mulberry genome analysis Int. J. Indust. Entomol. 10 79 - 86
Weng J. H. , Fan M. J. , Lin C. Y. , Liu Y. H. , Huang S. Y. 2007 Genetic variation of Zoysia as revealed by random amplified polymorphic DNA (RAPD) and isozyme pattern Plant. Prod. Sci. 10 80 - 85    DOI : 10.1626/pps.10.80
Xie Y. , Liu L. , Fu J. , Li H. 2012 Genetic diversity in Chinese natural zoysiagrass based on inter-simple sequence repeat (ISSR) analysis Afr. J. Biotechol. 11 7659 - 7669
Yaneshita M. , Kaneko S. , Sasakuma T. 1999 Allotetraploidy of Zoysia species with 2n=40 based on a RFLP genetic map Theor. Appl. Genet. 98 751 - 756    DOI : 10.1007/s001220051131
Yeh F. C. , Yang R. C. , Boyle T. 1999 POPGENE Version 1.31, Microsoft Windows-based Freeware for Population Genetic Analysis University of Alberta Alberta
Zietkiewicz E. , Rafalski A. , Labuda D. 1994 Genome fingerprinting by simple sequence repeat (SSR)-anchored polymerase chain reaction amplification Genomics 20 176 - 183    DOI : 10.1006/geno.1994.1151