1. Introduction
===============

*Cynodon dactylon* (L.) Pers. (common bermudagrass) is the most important member of the genus *Cynodon* because of its widespread distribution in warmer parts of the world and its use as livestock herbage and turf \[[@b1-ijms-13-00383]\]. Harlan and de Wet describe the taxon as the ubiquitous, cosmopolitan weed of the world, containing an enormous variation ranging from small, fine turfgrasses used as golf course putting green turf to robust types grown for pasture or hay \[[@b2-ijms-13-00383]\].

The center of diversity for some races of *Cynodon dactylon* (L.) Pers. was reported to be Turkey, Iran, Afghanistan and the western part of Pakistan, which could be potentially rich in the desirable genes and genotypes \[[@b1-ijms-13-00383],[@b2-ijms-13-00383]\]. Characterizing its genetic diversity for germplasms is an essential step in selection and breeding of this grass.

The commonly used polymerase chain reaction (PCR) based DNA marker systems are random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP) and more recently simple sequence repeats (SSRs) or microsatellites \[[@b3-ijms-13-00383]\]. The major limitations of these methods are low reproducibility of RAPD, high cost of AFLP and the need to know the flanking sequences to develop species specific primers for SSR polymorphism \[[@b4-ijms-13-00383]\]. Markers such as inter-simple sequence repeat (ISSR) \[[@b5-ijms-13-00383]\] are widely used in genetic diversity studies because they need no prior DNA sequence information, development costs are low, and laboratory procedures can easily be transferred to any plant species \[[@b6-ijms-13-00383]\]. ISSR is a technique that overcomes most of these limitations \[[@b5-ijms-13-00383],[@b7-ijms-13-00383]\]. This technology has been used to DNA fingerprint a wide range of crops \[[@b5-ijms-13-00383],[@b7-ijms-13-00383],[@b8-ijms-13-00383]--[@b10-ijms-13-00383]\] and to understand ploidy complex and the geographic origin of some plant species \[[@b11-ijms-13-00383],[@b12-ijms-13-00383]--[@b16-ijms-13-00383]\]. DNA profiling techniques that have been successfully used in assessing relatedness of *Cynodon* accessions includes DNA amplification fingerprinting (DAF) \[[@b17-ijms-13-00383]\], RAPD \[[@b18-ijms-13-00383],[@b19-ijms-13-00383]\], AFLP \[[@b20-ijms-13-00383],[@b21-ijms-13-00383]\], ISSR \[[@b22-ijms-13-00383]\], SSR \[[@b23-ijms-13-00383]\] and chloroplast specific simple sequence repeat length polymorphism (CpSSRLP) \[[@b20-ijms-13-00383]\].

In this study, we reported the feasibility of the ISSR-based PCR as an interesting approach in providing accurate molecular markers for investigating the genetic diversity among bermudagrass accessions.

2. Results
==========

2.1. ISSR Amplification
-----------------------

The fourteen selected ISSR primers produced 389 bands with an average of 27.78 bands per primer, of which 313 (80.5%) were polymorphic with an average of 22.36 bands per primer. The different ISSR primers amplified the number of bands from 20 (ISSR-14) to 40 (ISSR-11) having a range of 400 to 1200 bp fragments size. The number of polymorphic fragments detected by each primer ranged from 16 to 27. The most polymorphism was shown by ISSR-5, which showed 92.8% polymorphism. The mean of polymorphism information content (PIC) value for primers was 0.328 which ranged from 0.205 (ISSR-5) to 0.433 (ISSR-15) ([Table 1](#t1-ijms-13-00383){ref-type="table"}). [Figure 1](#f1-ijms-13-00383){ref-type="fig"} shows the banding pattern of genotypes generated by the primer ISSR-7.

2.2. Genetic Variation and Similarity among Genotypes
-----------------------------------------------------

The average values of observed number of alleles (Na), effective number of alleles (Ne), Nei's gene diversity (He) and Shannon's information index (I) for all primers were 2, 1.56, 0.33 and 0.56, respectively.

The similarity matrix of the 27 genotypes from ISSR data using Jaccard's coefficient showed that the genetic relationships of genotypes were different varying between 0.057 and 0.524 (data not shown). The lowest similarity was revealed between "18-Gr" and "27-II" genotypes (0.057), whereas the highest similarity was between "3-Gs" and "6-Gs" genotypes (0.524). The latest two common bermudagrass accessions, "3-Gs" and "6-Gs", are tetraploid and are collected from the same origin ([Table 2](#t2-ijms-13-00383){ref-type="table"}) with some similar morphological characteristics such as length and width of leaf and diameter of rhizome. The very close genetic relationship of these accessions suggests that they are a clonal propagule of a single plant.

2.3. Cluster Analysis
---------------------

The data obtained from ISSR analysis of 27 bermudagrass genotypes was subjected to cluster analysis. In order to recognize the best clustering and similarity coefficient methods and the cophenetic correlation coefficient, a measure of the correlation between cophenetic matrix constructed from tree and similarity matrix, was calculated for each method combination. From different methods, the highest value (*r* = 0.828) was observed for the un-weighted paired group method with arithmetic average (UPGMA) clustering method based on Jaccard's similarity coefficient, suggesting that the cluster analysis represents the similarity matrix.

The UPGMA dendrogram grouped the 27 accessions into six clusters considered as A, B, C, D, E and F with a similarity coefficient of 0.23 ([Figure 2](#f2-ijms-13-00383){ref-type="fig"}). Seven tetraploid accessions were grouped together as cluster A of which three ("1-Lr", "3-Gs", "6-Gs") and four ("11-Ns", "13-Ns", "15-Bs", "58-Chr") of them were from Charmahal Bakhtiari and Isfahan provinces, respectively. These accessions have some similar morphological characteristics such as leaf and rhizome color. Genetic similarity coefficient values for accessions in cluster A ranged from 0.198 to 0.524.

The major cluster B included the largest number of accessions (11 tetraploides) that had some similar morphological characteristics, such as leaf width, length of flower, stem and leaf, and stolon color. The grouping of these accessions from the same province (Isfahan) may reflect their having a common geographic origin. Genetic similarity coefficient values in this group were from 0.128 to 0.458.

Cluster C embraced two diploid accessions ("17-GN1" and "88-Khl") and one triploid accession ("44-II"). These accessions have also similar morphological characteristics such as width of leaf and stolon color.

Four triploid hybrids (*Cynodon dactylon* × *Cynodon transvaalensis*), Tifgreen, Tifdwarf, Tifway and Midlawn were separated from all accessions and grouped together as cluster D.

Two tetraploid accessions, "19-Gcr" and "18-Gr", collected from Gilan (Annual rainfall 1220 mm) were separated from all accessions in single clusters as labeled E and F, respectively in [Figure 2](#f2-ijms-13-00383){ref-type="fig"}.

Principal coordinate analysis (PCoA) based on genetic similarity metrics was used to visualize the genetic relationships among accessions. The first three eigenvectors accounted for 25.61% of the total molecular variation. Since the original data are not highly correlated in PCoA, the first few PCs do not explain much of the original variation. Assessment of genetic relationships on the basis of the first three PCs could lead to misleading interpretations and therefore analysis of genetic relationships among accessions should be based on cluster analysis as well as the optimal number of PCs that explain the maximum amount of the original data.

3. Discussion
=============

The results obtained in our study showed that ISSR markers measured sufficient polymorphism for DNA typing and can be successfully used in bermudagrass germplasm characterization and for fingerprinting purposes. Similar results were obtained in several studies using other molecular markers such as RAPD \[[@b11-ijms-13-00383],[@b18-ijms-13-00383]\].

The polymorphic information content (PIC) and the polymorphism rate (*P*) were used to measure the genetic diversity in bermudagrass accessions. High, medium or low polymorphism is in accordance with PIC \> 0.5, 0.5 \> PIC \> 0.25 and PIC \< 0.25, respectively \[[@b24-ijms-13-00383],[@b25-ijms-13-00383]\]. Moreover, the mean value of the PIC obtained in this study was 0.328, indicating that the primers could develop medium polymorphism which is useful for genetic variation of genotypes studied in this research.

The average values of Na, Ne, He and I indicated that the background genetic data of germplasm accessions should be considered for integrated application in the breeding programs for germplasm improvement. It should be noted that, as a measure of genetic variation, the effective number of alleles (Ne) or expected heterozygosity (He) is more appropriate than the actual number of alleles (Na), since the latter depend on the sample size \[[@b26-ijms-13-00383]\].

Bootstrapping was effectively utilized for estimating the statistical support to the internal branches of the tree. The internal tree branches that have \>70% bootstrap are likely to be correct at the 95% level. However, a high bootstrap percentage still does not guarantee that long branch attractions have not biased the results. Also, in many cases, the overall tree structure provides better information than a particular branch. Wherever a priori information regarding genetic relationships is available, it is preferable to apply it in the interpretation \[[@b27-ijms-13-00383]\]. Therefore, a priori information about accessions was used in branching the constructed tree.

Hybridization is one of the most common methods used to create variation in bermudagrass \[[@b28-ijms-13-00383]\]. Most of the bermudagrass cultivars are interspecific hybrids between *Cynodon dactylon* and *Cynodon transvaalensis*. In this study, we included four hybrids, which were divided into two subgroups on the UPGMA tree. Among these hybrids, "Tifgreen" and its putative spontaneous somatic mutant (Tifdwraf) were placed in two separate subgroups, whereas, "Tifway" and "Midlawn" were placed in the same subgroup. Zhang *et al*. and Caetano-Anolles used AFLP and DAF techniques to effectively differentiate "Tifgreen" and the putative somatic mutants \[[@b29-ijms-13-00383],[@b30-ijms-13-00383]\]. They showed that "Tifgreen" and "Tifdwarf" formed separate groups that are consistent with the results of this study. Triploidy forms were not indicated for *Cynodon dactylon* in the taxonomic revision of the genus as listed by Harlan *et al*. \[[@b1-ijms-13-00383]\] probably because their collection contained none \[[@b11-ijms-13-00383],[@b31-ijms-13-00383]\].

Cluster analysis was able to cluster bermudagrass accessions based on geographical origin, which is in agreement with previous observations of Li *et al*. in Chinese natural bermudagrass (*Cynodon dactylon*) germplasm using ISSR markers \[[@b22-ijms-13-00383]\]. This showed that there was an incomplete direct relationship between the origin of the accessions and the molecular clusters. The reasons might be due to genetic overlap occurring in the bermudagrass accessions from two different regions, open pollination behavior of *Cynodon* plants and exchange of germplasm resources from different regions. This constructed dendrogram may be able to discriminate between genetically different chromosome numbers of *Cynodon dactylon* as reported by Anderson *et al*. \[[@b17-ijms-13-00383]\], Etemadi *et al*. \[[@b18-ijms-13-00383]\] and Gulsen *et al*. \[[@b11-ijms-13-00383]\].

4. Experimental Section
=======================

4.1. Plant Materials and DNA Extraction
---------------------------------------

Twenty seven bermudagrass genotypes were included in this study. Among them, 23 genotypes were common bermudagrass types (*Cynodon dactylon*) and four were triploid hybrids of *Cynodon dactylon × Cynodon transvaalensis* (Tifdwarf, Tifgreen, Tifway and Midlawn). Stolons of twenty-three accessions of *Cynodon dactylon* (L.) Pers. had been collected from Isfahan, Charmahal Bakhtiari, Gilan and Mazandaran, provinces of Iran ([Table 2](#t2-ijms-13-00383){ref-type="table"}). Each accession was grown in 20 cm diameter pots under green-house conditions. The green-house was maintained at 25 ± 1 °C. Bermudagrass DNA samples were isolated from fresh leaf tissue according to Dellaporta *et al*. \[[@b32-ijms-13-00383]\].

4.2. ISSR Assays
----------------

Fourteen random ISSR primers (Operan technologies, Inc, Japan) were used to amplify DNA fragments of the selected materials. These primers were selected based on preliminary screening, producing a higher level of polymorphism and reproducible fragment patterns ([Table 1](#t1-ijms-13-00383){ref-type="table"}). ISSR reactions were conducted in a volume of 25 μL consisting of 2.5 μL of 10× PCR buffer, 1.5 mM MgCl~2~, 1 U of *Taq* DNA polymerase (Roch Co. Germany), 200 μM of dNTPs, 1.25 μM of each primer, and template DNA of approximately 50 ng.

Amplifications were carried out in a Thermal Cycler (T-1 Theroblock, Germany) with the following PCR program: 3 min of initial denaturing at 94 °C, 40 cycles of 94 °C for 45 s, 1 min for annealing at the primer-specific melting temperature, and 72 °C for 4 min; followed by a final extension of 10 min at 72 °C. The amplified fragments were separated on a 1.5% agarose gels in 1× TBE buffer running at 65 V constant for 1.5--2 h and then stained by ethidium bromide (1 μg·mL^−1^) and photographed under UV light in a gel documentation system (Villber Lourmat, France)

4.3. Data Analysis
------------------

All clearly detectable ISSR product bands were scored as either presence (1), absence (0) or ambiguous (9) of each band for all accessions and the matrix of ISSRs data was assembled. All amplifications were repeated at least twice and only reproducible and well-defined bands were considered for analysis.

Polymorphism information content (PIC) was calculated using the formula PIC = 1 − ∑*pi*^2^, where *pi* is the frequency of the *i*th allele \[[@b33-ijms-13-00383]\]. Jaccard's similarity coefficient values for each pairwise comparison between accessions were calculated and a similarity coefficient matrix was constructed. This matrix was subjected to an un-weighted paired group method with arithmetic average (UPGMA) \[[@b34-ijms-13-00383]\] to generate a dendrogram using software NTSYS-pc Version 2.02 \[[@b35-ijms-13-00383]\]. The MXCOMP subroutine was used to calculate a cophenetic correlation coefficient between the similarity matrix and respective dendrogram derived matrix to measure suitability-of-fit \[[@b36-ijms-13-00383]\]. Winboot software \[[@b37-ijms-13-00383]\] was used for a bootstrap analysis with 1000 replicates to obtain the confidence of branches of the cluster tree. The frequency of occurrence of each marker in each genotype was computed, to render a matrix of 27 accessions by ISSR markers. These matrices were subjected to principal coordinate analysis (PCoA) to show the relationship between genotypes.

Observed number of alleles (Na), effective number of alleles (Ne), Nei's gene diversity (He) \[[@b38-ijms-13-00383]\] and Shannon's information index (I) \[[@b39-ijms-13-00383]\] were estimated for total genotypes using POPGENE software version 1.32 \[[@b40-ijms-13-00383]\].

5. Conclusion
=============

In summary, our results indicate that the level of polymorphism among bermudagrass accessions is appreciably high in Iran. The results of this study also suggest that the ISSR marker offers a powerful means to analyze the genetic diversity among accessions. Cluster analysis using the UPGMA method and Jaccard's similarity coefficient (*r* = 0.828) grouped the accessions into six main clusters according to some degree to geographical origin, their chromosome number and some morphological characteristics. The enlargement of the number of primer pairs and accessions would have provided more useful information to assist the classification of the bermudograss germplasm.

![Bands profiles generated by the primer ISSR-7 with the sequence 5′-(GA)8T-3′, marker 100(M), Tifdwarf, Tifway, Tifgreen, Midlawn, 1-Lr, 3-Gs, 6-Gs, 11-Ns, 13-Ns, 15-Bs, 17-GN1, 18-Gr, 19-Gcr, 24-II, 27-II, 36-Tr, 39-II, 44-II, 51-Chg, 52-Chg, 53-Chr, 55-Chg, 56-Chr, 58-Chr, 59-Chf, 78-Ts, 88-Khl.](ijms-13-00383f1){#f1-ijms-13-00383}

![UPGMA (un-weighted paired group method with arithmetic average) dendrogram depicting patterns of genetic diversity for 23 *Cynodon dactylon* (L.) and Tifdwarf, Tifway, Tifgreen, Midlawn using 313 ISSR markers.](ijms-13-00383f2){#f2-ijms-13-00383}

###### 

List of the inter-simple sequence repeat (ISSR) primers used in the present study, polymorphism number, rates and polymorphism information content (PIC) values.

  Primer name       Nucleotide sequence   Tm (°C)   No. of fragments amplified   No. of polymorphic fragments   Percentage of poymorphism   PIC value
  ----------------- --------------------- --------- ---------------------------- ------------------------------ --------------------------- -----------
  ISSR-1            5′-(AC)8 G-3′         54/9      23                           19                             82.6                        0.313
  ISSR-2            5′-(AC)8T-3′          49/2      29                           25                             86.2                        0.266
  ISSR-3            5′-(AC)8C-3′          53/3      25                           20                             80                          0.291
  ISSR-4            5′-(AG)8C-3′          46/8      24                           20                             83.3                        0.280
  ISSR-5            5′-(AG)8T-3′          42/5      28                           26                             92.8                        0.205
  ISSR-6            5′-(GA)8C-3′          43/3      25                           22                             88                          0.283
  ISSR-7            5′-(GA)8T-3′          42/9      36                           26                             72.6                        0.333
  ISSR-8            5′-(CT)8G-3′          44/9      30                           20                             66.7                        0.346
  ISSR-11           5′-(AG)8GT-3′         47/1      40                           27                             67.5                        0.333
  ISSR-14           5′-(AC)8GG-3′         57        20                           16                             80                          0.375
  ISSR-15           5′-(AC)8GA-3′         53/7      25                           20                             80                          0.433
  RAMP-TAG          5′-T(AG)9-3′          56/1      29                           24                             82.7                        0.428
  RAMP-GAC          5′-G(AC)9-3′          53/9      25                           22                             88                          0.391
  LK7               5′-CCA(CT)8-3′        56/1      30                           26                             86.7                        0.321
  **Total means**   **-**                 **-**     **389**                      **313**                        **80.5**                    **0.328**

###### 

List of the accessions of *Cynodon dactylon* (L.) Pers. used in the present study.

  Accession name   Collected site (province-city)   Altitude (m)   Annual Rainfall (mm)   Main daily temperature (°C)   Choromosome No.
  ---------------- -------------------------------- -------------- ---------------------- ----------------------------- -----------------
  1-Lr             Charmahal Bakhtiari-Lordegan     1900           551                    17/5                          36
  3-Gs             Charmahal Bakhtiari-Borujen      2200           245                    12/8                          36
  6-Gs             Charmahal Bakhtiari-Borujen      1930           245                    12/8                          36
  11-Ns            Isfahan-Natanz                   1800           198                    16/4                          36
  13-Ns            Isfahan-Natanz                   1800           198                    16/4                          36
  15-Bs            Isfahan-Natanz                   1300           198                    16/4                          36
  17-GN1           Mazandaran-Neka                  15             895                    16/8                          18
  18-Gr            Gilan-Chaboksar                  15-            1220                   15/9                          36
  19-Gcr           Gilan-Chaboksar                  15-            1220                   15/9                          36
  24-II            Isfahan-Isfahan                  2060           122                    16/2                          36
  27-II            Isfahan-Isfahan                  1577           122                    16/2                          36
  36-Tr            Isfahan-Isfahan                  1580           122                    16/2                          36
  39-II            Isfahan-Isfahan                  1573           122                    16/2                          36
  44-II            Isfahan-Isfahan                  1570           122                    16/2                          27
  51-Chg           Isfahan-Chadegan                 2150           317                    10/2                          36
  52-Chg           Isfahan-Chadegan                 2150           317                    10/2                          36
  53-Chr           Isfahan-Chadegan                 2150           317                    10/2                          36
  55-Chg           Isfahan-Chadegan                 2180           317                    10/2                          36
  56-Chr           Isfahan-Chadegan                 2180           317                    10/2                          36
  58-Chr           Isfahan-Chadegan                 2170           317                    10/2                          36
  59-Chf           Isfahan-Chadegan                 2160           317                    10/2                          36
  78-Ts            Isfahan-Tiran                    1900           175                    14/5                          36
  88-Khl           Isfahan-Khomeini shahr           1560           130                    15/8                          18
