- Open Access
Genotypic and phenotypic evaluation of off-type grasses in hybrid Bermudagrass [Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy] putting greens using genotyping-by-sequencing and morphological characterization
Hereditas volume 155, Article number: 8 (2018)
Interspecific hybrid bermudagrass [Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy] is one of the most widely used grasses on golf courses, with cultivars derived from ‘Tifgreen’ or ‘Tifdwarf’ particularly used for putting greens. Many bermudagrass cultivars established for putting greens can be genetically unstable and lead to the occurrence of undesirable off-type grasses that vary in phenotype. The objective of this research was to genetically and phenotypically differentiate off-type grasses and hybrid cultivars. Beginning in 2013, off-type and desirable hybrid bermudagrass samples were collected from golf course putting greens in the southeastern United States and genetically and phenotypically characterized using genotyping-by-sequencing and morphology.
Genotyping-by-sequencing determined that 11% (5) of off-type and desirable samples from putting greens were genetically divergent from standard cultivars such as Champion, MiniVerde, Tifdwarf, TifEagle, and Tifgreen. In addition, genotyping-by-sequencing was unable to genetically distinguish all standard cultivars from one another due to their similar origin and clonal propagation; however, over 90,000 potentially informative nucleotide variants were identified among the triploid hybrid cultivars.
Although few genetic differences were found in this research, samples harvested from golf course putting greens had variable morphology and were clustered into three distinct phenotypic groups. The majority of off-type grasses in hybrid bermudagrass putting greens were genetically similar with variable morphological traits. Off-type grasses within golf course putting greens have the potential to compromise putting surface functionality and aesthetics.
The economic impact of the golf industry in the United States (U.S.) in 2011 was estimated to be $176.8 billion with a contribution of approximately 1.98 million jobs (SRI International; http://wearegolf.org/economy/impact). Interspecific hybrid bermudagrasses (Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy) are some of the most widely utilized grasses on golf courses throughout tropical, subtropical, and temperate climates due to their desirable turfgrass characteristics of texture, color, and stress tolerance . Putting greens are a vital aspect of a golf course and in 2007 hybrid bermudagrasses were grown on 80% of putting green acreage in the southern U.S. .
‘Tifgreen’ was one of the first interspecific hybrids developed for putting green use [6, 21]. Soon after its commercial release, ‘Tifdwarf’ was selected from a somatic mutation in a ‘Tifgreen’ establishment [7, 8]. Despite being genetically unstable [9, 10], ‘Tifgreen’ and ‘Tifdwarf’ are used on putting surfaces; however, superior mutations have been released as new “ultradwarf” cultivars . Ultradwarf bermudagrass cultivars are classified due to their more diminutive morphology (i.e., internode length and leaf length and width) compared to ‘Tifdwarf’. ‘Champion’ and ‘MiniVerde’ were selected from somatic mutations in established ‘Tifdwarf’ plantings [5, 23], whereas ‘TifEagle’ was a putative mutant from radiation-induced ‘Tifgreen’ or ‘Tifway II’ rhizome [18, 20, 41]. In all cases, “ultradwarf” cultivars such as ‘Champion’, ‘MiniVerde’, and ‘TifEagle’ were once identified as off-types.
Morphological characterization has been considered a traditional method of studying turfgrass classification and diversity [3, 22, 25, 38]. Morphological characteristics such as internode length, leaf length, leaf width, and stolon diameter are of particular interest for classification of bermudagrasses [25, 26, 38] because differences in morphology can differentiate off-types from desirable cultivars [9, 10]. Additionally, Roche and Loch  stated that morphological characterization could provide useful information to further research for adaptation and management of different hybrid bermudagrasses. Researchers have used morphological characterizations to compare hybrid bermudagrass cultivars within the ‘Tifgreen’ family [31, 37]; however, morphological inconsistencies suggest that molecular techniques are also needed to accurately evaluate hybrid bermudagrass diversity .
Several methods have been used to explore genotypic differences among off-type grasses and hybrid bermudagrass cultivars. DNA amplification fingerprinting (DAF) and signatures from amplification profiles identified contaminant off-types not related to the ‘Tifgreen’ family, but could not distinguish mutant off-types within the ‘Tifgreen’ family . Amplified fragment length polymorphisms (AFLPs) studies determined the genetic diversity among several ‘Tifgreen’-derived cultivars. Studies utilizing AFLPs grouped ‘Tifgreen’, ‘Tifdwarf’, ‘TifEagle’, and ‘Champion’ into the same genetic cluster despite the grasses having differing phenotypic characteristics [11, 12, 41]. Attempts to use various simple-sequence repeats (SSRs) to identify ‘Tifgreen’-derived hybrid bermudagrass cultivars met with limited success as well [19, 24]. While SSRs identified ‘TifEagle’ from other ‘Tifgreen’-derived cultivars and identified polymorphisms unique to ‘Tifdwarf’ and ‘MiniVerde’ [19, 24], SSRs are not able to readily distinguish all ‘Tifgreen’-derived hybrid bermudagrass cultivars from one another nor have they been able to identify weedy off-type grasses from standard cultivars used on golf course putting greens. Fig. 1 displays the aesthetic and functional disruptions that off-type grasses in bermudagrass putting greens can create.
Genotyping-by-sequencing (GBS) is a high-throughput, next generation sequencing method capable of generating large numbers of single nucleotide polymorphism (SNPs) from species with high diversity . GBS offers several advantages over other molecular marker techniques including the amount data generated and the price per sample [13, 15]. Moreover, GBS allows analysis of a species (e.g., hybrid bermudagrass) for which a complete reference genome sequence is not available. Fiedler et al.  identified over 4600 high-quality SNPs in switchgrass (Panicum virgatum) using an early draft genome assembly with only half of the assembled DNA contigs scaffolded. Poland et al.  used GBS to map over 34,000 SNPs for the Oregon Wolfe Barley (Hordeum vulgare) reference population and 20,000 SNPs for the Synthetic W9784xOpata85 wheat (Triticum aestivum) reference population, both of which lacked a complete reference genome sequence. Based on the robustness of the technique and successful use in other grasses without a complete reference genome [13, 15, 32, 33], we hypothesize that GBS may be able to identify genetic variation among off-types and hybrid bermudagrasses used on putting greens. Therefore, our objectives were to explore the genetic and the phenotypic variation among off-type grasses sampled from hybrid bermudagrass putting greens using GBS and morphological characterization.
Single replicates of desirable and off-type hybrid bermudagrass samples were harvested in 2013 from putting greens on golf courses in Alabama, Arkansas, Florida, Mississippi, South Carolina, and Tennessee (Table 1). The greenkeeper at each golf course determined samples that were desirable from those that were off-type grasses. Samples were harvested with a 7.5 cm diameter tubular plugger (Turf Tec International; Tallahassee, FL, USA) and established using one three node stolon planted in a 64 cm2 pot filled with a peat moss based growing medium (Pro-Mix BX Mycorrhizae; Premier Horticulture, Inc.; Quakertown, PA, USA) in a glasshouse environment at the University of Tennessee (Knoxville, TN, USA; 35.5°N, −83.5°W). Plants were maintained with 24 kg N ha−1 wk.−1 of a water-soluble complete fertilizer (20 N-8.7P-16.6 K; Southern Agriculture; Hendersonville, NC, USA), irrigated to promote active growth, and insecticides (abamectin 0.01 kg ai ha−1, Avid 0.15EC, Syngenta; pymetrozine 0.35 kg ai ha−1, Endeavor, Syngenta) were applied on a preventive basis.
Ploidy levels were confirmed for each sample included in the GBS assay using flow cytometry at the University of Georgia Coastal Plain Experiment Station (Tifton, GA). Fresh leaf tissue was isolated from samples and chopped using a razor in 300 μL of LB01-lysis buffer (15 mM Tris, 2 mM Na2EDTA, 0.5 mM spermine-4HCl, 80 mM KCl, 20 mM NaCl, 0.1% v/v Triton X-100 pH 7.5 and 16 mM ß-mercaptoethanol) to release nuclei. Each bermudagrass sample was combined with Sorghum bicolor cv. BTx623 for a standard genome size comparison. Samples were passed through a 30-μm filter (CellTrics; Partec; Munster, Germany) and then 150 μL of RNase and propidium iodide solution (PI/RNase Staining Buffer, BD Biosciences, San Jose, CA, USA) was added. Samples were incubated on ice for 15 min and analyzed on an Accuri C6 flow cytometer (BD Biosciences; San Jose, CA, USA). Gating was set by the selection of objects that exhibited a strong correlation between the FL2 and FL3 signals using a flow rate of 14 μL per minute and a minimum cell count of 10,000. The mean FL2-A peaks from the signals were determined for S. bicolar and each unknown hybrid bermudagrass sample using Accuri C6 software (BD Biosciences; San Jose, CA, USA). These mean FL2-A values were then used with S. bicolar genome size (1.67 pg/2C) to calculate the genome size of each unknown hybrid bermudagrass sample .
Plant material and DNA isolation
Desirable and off-type samples labeled S1 to S47 in Table 1 were included in GBS due to the expense of the analysis. Hybrid bermudagrass cultivars [Champion (CH1–6), MiniVerde (MV1–6), Tifdwarf (TD1–6), TifEagle (TE1–6), Tifgreen (TG1–6), and Tifway (TW1–6)] and progenitor species [C. dactylon (TA1–3 and TB1–3) and C. transvaalensis (DA1–3 and DB1–3)] were used as standards in the GBS analysis. Only three biological replicates of the progenitor species were included also due to the expense of the analysis. Plant material for these standard entries was obtained from the University of Georgia Coastal Plain Experiment Station (Tifton, GA, USA).
For all samples, plant genomic DNA was isolated from actively growing leaf tissue on a single stolon using the Qiagen DNeasy Plant Mini Kit (Qiagen; Valencia, California, USA) according to the manufacturer’s protocol. DNA concentration was quantified using an intercalating dye (Quant-iT™ PicoGreen dsDNA Asasy Kit; Life Technologies; Carlsbad, CA, USA). DNA working solutions for the GBS protocol had a total volume of 30 μL and a concentration ranging from 50 to 105 ng μL−1.
Genotyping-by-sequencing was conducted at the Cornell University Institute for Biotechnology (Ithaca, NY) using the protocol described by Elshire et al. . ApeKI restriction enzyme was selected based on optimization trials for the GBS digestion to maximize the number of sampled genomic loci . Libraries for next-generation sequencing were constructed from DNA samples and multiplexed using Illumina HiSeq 2500 and then Illumina NextSeq 500 to increase read coverage and depth.
The combined Illumina data sets were initially analyzed with the UNEAK pipeline of the Tassel software package [4, 17]. One of the limitations of the UNEAK pipeline is that its nucleotide variant calling algorithm relies on a diploid model. Bermudagrasses sequenced in our GBS analysis included diploid (2n = 2× = 18), triploid (2n = 3× = 27), and tetraploid (2n = 4× = 36) samples; therefore, an alternative approach was used to call variants. Sequence tags with a predicted variant were extracted from the topm.bin libraries generated during the Tag-Pair-Export phase using the UNEAK Binary-to-text-plugin; then raw reads were mapped with Bowtie2 v2.2.7 to these tags generated from the UNEAK pipeline as a pseudo-reference .
The haplotype-based variants caller, Freebayes v1.0.2–15, was used to call variants for each set of samples with the same ploidy level with the correct ploidy level specified with parameter p . The sorting, indexing, and merging of alignment files was performed with the SAMtools v1.3 package . Multidimensional scaling (MDS) plots were generated from these variants using PLINK v1.9 to illustrate the variation among samples. The bermudagrass samples S19, S28, S30, S32, and S44 were not included because they had less than one million raw-reads . The individual samples for each triploid cultivar were pooled to increase the read depth for each cultivar. The read alignment files were pooled using SAMtools v1.3 package  and then the FreeBayes method of determining variants was utilized again. The pooled data was then used in a custom Python script to determine loci that differed between at least two cultivars (github.com/statonlab/UDBG_Informative_SNPs/blob/master/find_informative_SNPs.py). Any loci with three different genotypes (homozygous for the reference allele, homozygous for the alternate allele, or heterozygous) for at least two cultivars were flagged. Loci were not flagged if heterozygosity differed between or within the two subgenomes in triploid cultivars (i.e., 0/0/1 or 0/1/1). Lastly, each cultivar was filtered to use only individual variants with a read depth greater than or equal to 40. All raw read data has been submitted to NCBI under BioProject accession PRJNA353769.
Bermudagrass samples labeled S1 to S62, except for S13, S14, S19, S20, and S43 were used in phenotypic evaluation (Table 1). Bermudagrass samples S13, S14, S19, S20, and S43 were excluded from phenotypic evaluation due to their harvest date. The greenkeeper at each golf course determined samples that were desirable from those that were off-type grasses. Samples were harvested with a 7.5 cm diameter tubular plugger (Turf Tec International; Tallahassee, FL, USA) and established during the summer of 2013 using one three node stolon planted in a 64 cm2 pot filled with a peat moss based growing medium (Pro-Mix BX Mycorrhizae; Premier Horticulture, Inc.; Quakertown, PA, USA) in a glasshouse environment at the University of Tennessee (Knoxville, TN, USA; 35.5°N, −83.5°W). Plants were maintained with 24 kg N ha−1 wk.−1 of a water-soluble complete fertilizer (20 N-8.7P-16.6 K; Southern Agriculture; Hendersonville, NC, USA), irrigated to promote active growth, and insecticides (abamectin 0.01 kg ai ha−1, Avid 0.15EC, Syngenta; pymetrozine 0.35 kg ai ha−1, Endeavor, Syngenta) were applied on a preventive basis.
For phenotypic evaluation, one three-node stolon of each sample was established in four, 64 cm2 pots filled with a peat moss based growing medium (Pro-Mix BX Mycorrhizae; Premier Horticulture, Inc.; Quakertown, PA, USA) on 7 April 2014. The stolon length at transplanting of the 52 selections ranged from 3.4 to 11.3 cm. The plants were maintained as previously described but regular clipping was ceased 2 weeks prior to evaluation.
Morphological measurements and statistical analysis
Phenotypic evaluation of off-type and desirable samples was conducted by measuring plant morphological characteristics via methods outlined by Roche and Loch . Five parameters were assessed and included internode length and stolon diameter, leaf length and width, and the leaf length:width ratio (LWR). Measurements were made between the third and fourth node and on the outer leaf from the third node using digital calipers (Digimatic Caliper, Model No.CD-6″ CX, Mitutoyo Corporation, Kawasaki-shi, Kanagawa, Japan). The experiment was a completely randomized design with pots replicated four times and morphology measured on three stolons per pot. Morphology was assessed on 3 June 2014 and repeated again on 25 June 2014.
All morphological data describing desirable and off-type hybrid bermudagrass samples were analyzed using cluster analysis in SAS Enterprise Guide (Version 6.1, SAS Institute, Cary, NC, USA). K-means clustering algorithm was used to partition the data set into a user-defined number of clusters . Three clusters were determined based on the cubic clustering criterion and the frequency of observations in each cluster . Cluster means and standard deviations for each morphological measurement were then graphed in Prism (Prism 6 for Mac OS X; GraphPad Software, Inc.) to determine statistical differences among cluster means.
Over 878 million reads were generated through Illumina sequencing, with 271 million from the Illumina HiSeq platform and an additional 606 million from the same libraries run on the Illumina NextSeq platform. Of the 878 million reads, 757 million were determined to be of high quality and were able to be demultiplexed and assigned to individual plant samples. After excluding five samples with less than a million reads, the remaining individual samples had a range of 1.3 (sample S44) million reads to 13.9 million reads (sample S22). The variants yielded from GBS analysis include single nucleotide variants, multiple nucleotide variants, and indels. An initial 1,088,920 total variants were identified with an average read depth of 4.9 sequences for the triploid bermudagrass (C. dactylon x C. transvaalensis) samples. GBS analysis identified 347,512 total variants with an average read depth of 9.5 per individual for the two diploid, C. transvaalensis selections. For the tetraploid, C. dactylon samples, 587,053 total variants were identified with an average sequence read depth of 7.4 (Fig. 2). Only 136,205 variants were shared among diploid, triploid, and tetraploid species; therefore, the remaining variants are fixed in at least one species (Fig. 2).
The majority of samples harvested from golf courses clustered with the hybrid bermudagrass cultivars Champion, MiniVerde, Tifdwarf, TifEagle, and Tifgreen in the MDS plot (Fig. 3). The clustering suggested that samples from golf courses were genetically similar to those hybrid bermudagrass cultivars typically established on golf course putting greens (Fig. 3). Of the 47 unknown samples, only five (~11%) were genetically divergent from the standard cultivars (S4, S16, S31, S33, and S45), as illustrated by the MDS plot (Fig. 3). Hybrid bermudagrass cultivars Champion, MiniVerde, Tifdwarf, TifEagle, and Tifgreen were genetically similar to one another; however, GBS separated these cultivars from ‘Tifway’ hybrid bermudagrass (Fig. 3). Pooling individual cultivar samples yielded a higher average read depth of 31 per variant site per cultivar. Using the pooled data, 675,578 loci were identified as different between at least two cultivars. The majority of these genotype differences were only able to differentiate cultivars within the ‘Tifgreen’-cultivar family from those with different lineage (i.e., ‘Tifway’) (Table 2). Table 2 exhibits the number of nucleotide variants with different genotypes between each pair of triploid hybrid bermudagrass cultivars. Variants were included if they were homozygous for the reference allele in one cultivar and homozygous for the alternate allele in the other cultivar or if they were homozygous in one cultivar and heterozygous in the other. ‘Tifway’ has a much larger number of variants compared to the other five cultivars (29,614 variants when compared to ‘MiniVerde’ to 37,802 when compared to ‘Tifdwarf’). The five other cultivars are more similar to each other, with the highest number of identified variants of 4476 between ‘Tifdwarf’ and ‘Tifgreen’.
The pooled data from variants among triploid cultivars still encompassed individual loci with both very low and very high individual read depths. Low read depth could miss heterozygotes, whereas high read depth could indicate a repetitive region instead of an individual locus. To mitigate read depth issues, a further filter was applied to identify only the most robust variants with a sequence read depth of at least 40, but no more than 100. The upper limit was set to filter reads originating from repetitive elements where detected variation is not likely to be from a single locus. Filtering using these read depths yielded 93,188 variants between at least two genotypes (Table 2).
The MDS plot revealed clear clustering of the diploid (DA1–3 and DB1–3) and tetraploid progenitor species (TA1–3 and TB1–3) apart from the standard hybrid bermudagrass cultivars (CH1–6, MV1–6, TD1–6, TE1–6, TG1–6, and TW1–6) and the majority of samples harvested from putting greens (Fig. 3a). The two progenitor species also clustered separately from one another with the exception of DA2 due to possible contamination during DNA isolation. The clustering demonstrated that GBS was effective for distinguishing diploid, triploid, and tetraploid bermudagrasses. The ability to distinguish among bermudagrass species is likely due to the large number of unshared variants (Fig. 2).
The K-means cluster algorithm yielded three clusters containing 14, 26, and 12 samples, respectively. Cluster one contained nine off-type and five desirable samples, cluster two had 12 off-types and 14 desirables, and cluster three had eight off-types and four desirable samples. The cluster analysis overall expected R 2 was 0.61 with a cubic clustering criterion of −19.36. Cluster means and standard deviations for each morphological assessment are presented in Fig. 4. Internode length, leaf length, and LWR were the only statistically different morphological parameters among clusters (Fig. 4). A representative hybrid bermudagrass sample from each cluster is illustrated in Fig. 5.
The average internode length for grasses in cluster one (34.6 mm) was significantly longer than the grasses within clusters two (21.9 mm) and three (24.7 mm) (Fig. 4). Grasses in cluster three had significantly longer leaves than those grouped in clusters one and two (Fig. 4). The leaf length mean for cluster three was 29.8 mm, compared to 14.9 and 9.9 mm for clusters one and two, respectively (Fig. 4). This relationship was also present in LWR among clusters (Fig. 3). Stolon diameter ranged from 0.7 to 0.8 mm and leaf width ranged from 2.0 to 2.2 mm, with no statistical differences present among clusters for either parameter.
Caetano-Anollés  and Caetano-Anollés et al.  revealed that eight of 16 off-types were genetically divergent from standard cultivars using DAF, leading researchers to conclude that off-types that were not genetically distinct, but were the result of somatic mutations within ‘Tifgreen’ and ‘Tifdwarf’. The inability of GBS, as well as other molecular marker techniques, to distinguish off-types from hybrid bermudagrass cultivars used on putting greens could be the result of aneuploidy within the ‘Tifgreen’-cultivar family (B.M. Schwartz, unpublished data, 2016; ). Reasor et al.  suggested that aneuploidy could have resulted during the original hybridization of ‘Tifgreen’ or through intense putting green management techniques. It is expected that some variant locations are not going to be sampled by random chance due to the sparse nature of the GBS analysis. This is a limitation of GBS because it cannot determine presence/absence or copy number variations for individual locations that are needed to determine aneuploidy . There were a total of 93,188 variants shared between at least two genotypes (Table 2); however, our experiment only included six biological replications of each standard cultivar from a single geographic location (Tifton, GA). Additional research and replication of this study with more samples will be needed to ascertain which variants can be used to identify standard hybrid bermudagrass cultivars, specifically hybrid ultradwarf cultivars, from one another.
Despite also being a triploid hybrid, ‘Tifway’ bermudagrass has been genetically distinguished from ‘Tifgreen’-cultivar family using SSRs [20, 24, 40] and AFLPs [12, 41]. Arumuganathan et al.  reported that ‘Tifway’ had less nuclear DNA content (1.37 ± 0.01 pg/2C) than ‘Tifgreen’ (1.61 ± 0.00 pg/2C) despite having the same number of chromosomes. Furthermore, Reasor et al.  hypothesized that this difference in DNA content could also be a result of aneuploidy in the ‘Tifgreen’-cultivar family of hybrid bermudagrass, which includes the hybrid ultradwarf cultivars. This difference in DNA content could also aid in genetic identification between hybrid bermudagrasses using GBS. In addition to the identification of ‘Tifway’ from other triploid hybrids, the ability of GBS to distinguish diploid, triploid, and tetraploid bermudagrasses align with previous efforts to genetically identify these grasses from one another using AFLPs [12, 41] and SSRs [20, 24, 40].
It is not clear why the majority of grasses included in our experiment exhibited variable morphological characteristics while being similar in genotype. The majority of bermudagrass cultivars established on golf course putting greens were selected from other bermudagrass cultivars . The off-type grasses studied in this experiment were also selected from existing cultivars. Differential gene expression driven by epigenetic mechanisms such as DNA methylation, histone modification, and small RNA expression may also be a possible explanation for the genetic similarities among hybrid bermudagrass samples varying in phenotype . Multiple genes control important turfgrass traits and gene expression can be greatly influenced by environment or management practices . Golf course putting greens are intensely managed surfaces subjected to daily mowing (often at heights of cut ≤3 mm), annual aerification and cultivation, as well as treatment with plant growth regulators and silica sand topdressing on a weekly basis. Any of these practices or other environmental influences could cause lasting epigenetic effects that result in the up or down regulation of genes associated with hybrid bermudagrass phenotypic characteristics; however, no research has been conducted on this possibility. Studying changes in gene expression as a result of these maintenance practices could benefit researchers and industry practitioners to better understand how putting green management could potentially lead to the occurrence of phenotypically different off-type grasses in hybrid bermudagrass putting surfaces and generate new hypotheses into how these changes are induced.
The internode length of grasses in this experiment align with Magni et al.  who reported an internode length range of 15 to 34 mm on ultradwarf hybrid bermudagrass cultivars used on putting greens. However, Roche and Loch  reported internode lengths of 9.4 to 12.5 mm for hybrid bermudagrasses used on putting greens. In our experiment, mean internode length for each cluster was in the uppermost half of the internode length range reported by Magni et al.  and greater than the range measured by Roche and Loch . Internode lengths measured in this experiment varied greatly among desirable and off-type grasses as well as grasses measured in other experiments. This is an indication of the amount of phenotypic variability that can occur in individual putting greens as well as from golf course-to-golf course and cultivar-to-cultivar. Differences in internode length within the same putting surface can lead to decreased turfgrass density and reductions in putting surface quality and playability . Similar to the internode length data, leaf length values (and subsequently LWR values) documented in our experiment were far greater than those reported by Roche and Loch . Stolon diameter and leaf width values were similar to those reported by Roche and Loch , but less than those reported by Magni et al. .
Off-type grasses reported to have phenotypic differences from standard hybrid bermudagrass cultivars were sampled from golf course putting greens and subjected to GBS and morphological characterization under controlled growth conditions. Genotyping-by-sequencing only distinguished five off-type grasses from standard hybrid bermudagrass cultivars. In addition, GBS failed to completely distinguish standard hybrid bermudagrass cultivars from one another, including ‘Champion’, ‘MiniVerde’, ‘Tifdwarf’, ‘TifEagle’, and ‘Tifgreen’. These results are not unexpected given their common origin. The final bioinformatics analysis did yield 93,188 variants that offer the potential to be useful in distinguishing standard cultivars from one another; however, additional research beyond the scope if this project would be needed to determine which ones are diagnostic. GBS was successful in determining triploid hybrid bermudagrass cultivars from two diploid and tetraploid progenitor samples. Additionally, GBS was also successful in determining triploid hybrid bermudagrass cultivars with lineage to ‘Tifgreen’ from those not developed from ‘Tifgreen’ (e.g., ‘Tifway’). Morphological characteristics varied among sampled grasses that allowed them to be clustered into three distinct phenotypic groups varying predominately in internode and leaf length.
Arumuganathan K, Tallury SP, Fraser ML, Bruneau AH, Qu R. Nuclear DNA content of thirteen turfgrass species by flow cytometry. Crop Sci. 1999;39:1518–21.
Beard JB. Turf management for golf courses 2nd ed. Chelsea: Wiley; 2002.
Bonos SA, Plumley KA, Meyer WA. Ploidy determination in Agrostis using flow cytometry and morphological traits. Crop Sci. 2002;42:192–6.
Bradbury PJ, Zhang Z, Kroon DE, Terry Casstevens M, Ramdoss Y, Buckler ES. TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5.
Brown RM, Brown MA, Brown SD (1997) ‘champion’ dwarf’ hybrid bermudagrass. US Plant patent 9888 date issued: 6 may 1997.
Burton GW. ‘Tifgreen’ (Tifton 328) bermudagrass for golf greens. USGA Green Section Record. 1964;2:11–3.
Burton GW. ‘Tifdwarf’ bermudagrass. Crop Sci. 1966;6:94.
Burton GW, Elsner JE. ‘Tifdwarf’-a new bermudagrass for golf greens. USGA Green Section Record. 1965;2:8–9.
Caetano-Anollés G. Genetic instability of bermudagrass (Cynodon) cultivars ‘Tifgreen’ and ‘Tifdwarf’ detected by DAF and ASAP analysis of accessions and off-types. Euphytica. 1998;101:165–73.
Caetano-Anollés G, Callahan LM, Gresshoff PM. The origin of bermudagrass (Cynodon) off-types inferred by DNA amplification fingerprinting. Crop Sci. 1997;37:81–7.
Capo-chichi LJ, Goatley JM, Philley W, Krans J, Davis D, Kato A, Santen E. Dinitroaniline-induced genetic changes in bermudagrass. Crop Sci. 2005;45:1504–10.
Chen Z, Wang ML, Waltz C, Raymer P. Genetic diversity of warm-season turfgrass: seashore paspalum, bermudagrass, and zoysiagrass revealed by AFLPs. Floriculture Ornamental Biotechnol. 2009;3:20–4.
Elshire RJ, Glaubitz JC, Sun G, Poland JA, Kawamoto K, Buckler ES, Mitchell SE. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011;6:E19379. doi:10.1371/journal.pone.0019379.
Fei S. Recent progresses on turfgrass molecular genetics and biotechnology. Acta Hortic. 2008;783:247–60.
Fiedler JD, Lanzatella C, Okada M, Jenkins J, Schmutz J, Tobias CM. High-density single nucleotide polymorphism linkage map of lowland switchgrass using genotyping-by-sequencing. Plant Genome. 2015;8. doi:10.3835/plantgenome2014.10.0065.
Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. Cornell University Library; 2012. https://arxiv.org/abs/1207.3907.
Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, Buckler ES. TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline. PLoS One. 2014;9(2):e90346.
Hanna WW, Elsner JE. Registration of ‘TifEagle’ bermudagrass. Crop Sci. 1999;39:1258.
Harris-Shultz KR, Schwartz BM, Brady JA. Identification of simple sequence repeat markers that differentiate bermudagrass cultivars derived from ‘Tifgreen’. J Amer Soc Hort Sci. 2011;136:211–8.
Harris-Shultz KR, Schwartz BM, Hanna WW, Brady JA. Development, linkage mapping, and use of microsatellites in bermudagrass. J Amer Soc Hort Sci. 2010;135:511–20.
Hein MA. Registration of varieties and strains of bermudagrass, III. (Cynodon dactylon (L.) Pers.). Agron J. 1961;53:276.
Hitchcock AS. Manual of the grasses of the United States. 2nd ed. New York: Dover Publications; 1971.
Kaerwer HE, Kaerwer B (2001) Bermudagrass ‘P-18’. United States plant patent. 12084P3 date issued: 11 September.
Kamps TL, Williams NR, Ortega VM, Chamusco KC, Harris-Shultz KR, Scully BT, Chase CD. DNA polymorphisms at bermudagrass microsatellite loci and their use in genotype fingerprinting. Crop Sci. 2011;51:1122–31.
Kang SY, Lee GJ, Lim KB, Lee HJ, Park IS, Chung SJ, Kim JB, Kim DS, Rhee HK. Genetic diversity among Korean bermudagrass (Cynodon spp.) ecotypes characterized by morphological, cytological, and molecular approaches. Mol Cell. 2008;25:163–71.
Kenworthy KE, Taliaferro CM, Carver BF, Martin DL, Anderson JA. Genetic variation in Cynodon transvaalensis Burtt-Davy. Crop Sci. 2006;46:2376–81.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Lyman GT, Throssell CS, Johnson ME, Stacey GA, Brown CD. Golf course profile describes turfgrass, landscape, and environmental stewardship features. Applied Turfgrass Science. 2007; doi:10.1094/ATS-2007-1107-01-RS.
MacQueen JB. Some methods for classification and analysis of multivariate observations. Proceedings of the fifth symposium on math, statistics, and probability. Berkeley: University of California Press; 1967. p. 281–97.
Magni S, Gaetani M, Caturegli L, Leto C, Tuttolomodo T, La Bella S, Virga G, Ntoulas N, Volterrani M. Phenotypic traits and establishment speed of 44 turf bermudagrass accessions. Acta Agric Scand Sect B Soil Plant Sci. 2014;64:722–33. http://dx.doi.org/10.1080/09064710.2014.955524
Poland JA, Brown PJ, Sorrells ME, Jannink JL. Development of high-density genetic maps for barley and wheat using a novel two enzyme genotyping-by sequencing approach. PLoS One. 2012;7:E32253. doi:10.1371/journal.pone.0032253.
Poland JA, Rife TW. Genotyping-by-sequencing for plant breeding and genetics. Plant Genome. 2012;5:92–105. doi:10.3853/plantgenome2012.05.0005.
Price HJ, Dillon SL, Hodnett G, Rooney WL, Ross L, Johnston JS. Genome evolution in the genus Sorghum (Poaceae). Ann Bot. 2005;95:219–27. doi:10.1093/aob/mci015.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Amer J Human Genet. 2007;81:559–75.
Reasor EH, Brosnan JT, Trigiano RN, Elsner JE, Henry GM, Schwartz BM. The genetic and phenotypic variability of interspecific hybrid bermudagrasses (Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy) used on golf course putting greens. Planta. 2016;244:761–73. doi:10.1007/s00425-016-2571-8.
Roche MB, Loch DS. Morphological and development comparisons of seven greens quality hybrid bermudagrass [Cynodon dactylon (L.) Pers. X C. transvaalensis Burtt-Davy] cultivars. Intl Turfgrass Soc Res J. 2005;10:627–34.
Romani M, Piano E, Carroni AM, Pecetti L. Evaluation of native bermudagrass (Cynodon dactylon) germplasm from Italy for the selection of adapted turfgrass cultivars. Acta Hort. 2004;661:381–6.
Sahu PP, Pandey G, Sharma N, Puranik S, Muthamilarason M, Prasad M. Epigentic mechanisms of plant stress responses and adaptation. Plant Cell Rep. 2013;32:1151–9. doi:10.1007/s00299-013-1462-x.
Wang Z, Wu Y, Martin DL, Gao H, Samuels T, Tan C. Identification of vegetatively propagated turf bermudagrass cultivars using simple sequence repeat markers. Crop Sci. 2010;50:2103–11.
Zhang LH, Ozias-Akins P, Kochert G, Kresovich S, Dean R, Hanna WW. Differentiation of bermudagrass (Cynodon spp.) genotypes by AFLP analyses. Theor Appl Genet. 1999;98:895–902.
The authors would like to thank Qi Sun, Jeffe Dunne, Sarah Boggess, Anne Hatmaker, Monil Mehta, Laura Poplawski, Sujata Agarwal, Javier Vargas, Tyler Campbell, Jimmy Greenway, Greg Breeden, Daniel Farnsworth, Trevor Mills, Mitchell Riffey, Kelly Arnholt, Jason Burris, Amanda Webb, Rebecca Grantham, Gerald Henry, and John Sorochan for assistance with this project. Mention of trade names or commercial products in this article is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the University of Tennessee Institute of Agriculture or the U.S. Department of Agriculture. USDA is an equal opportunity provider and employer.
The University of Tennessee Institute of Agriculture and the Memphis Area chapter of the Golf Course Superintendents Association of America provided funding for this project.
Availability of data and materials
All raw sequence data for bermudagrass samples have been submitted to NCBI under BioProject accession PRJNA353769.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Reasor, E.H., Brosnan, J.T., Staton, M.E. et al. Genotypic and phenotypic evaluation of off-type grasses in hybrid Bermudagrass [Cynodon dactylon (L.) Pers. x C. transvaalensis Burtt-Davy] putting greens using genotyping-by-sequencing and morphological characterization. Hereditas 155, 8 (2018). https://doi.org/10.1186/s41065-017-0043-3
- Putting greens