|
|
||||||||
Dep. of Plant Pathology, Univ. of Wisconsin, 1630 Linden Dr., Madison, WI 53706
* Corresponding author (jung{at}plantpath.wisc.edu).
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: MDS, multidimensional scaling plot RAPD, random amplified polymorphic DNA
| INTRODUCTION |
|---|
|
|
|---|
The above advantages, along with improvement of intraspecific hybridization techniques (Pepin and Funk, 1971) have caused an increased interest in Kentucky bluegrass (Bonos et al., 2000a), resulting in a rapid increase in the number of registered Kentucky bluegrass cultivars. These new varieties are in addition to the many genotypes resulting from selection of promising, usually apomictic plants from old fields (Pepin and Funk, 1971) and other natural populations (Bonos et al., 2000b). Further, this species, and hence the genotypes within it, are inherently diverse because of adaptation to many climates and locations over many years (Burt and Christians, 1990).
The large, increasing number of Kentucky bluegrass cultivars can make registration of new cultivars difficult. In Europe, new Brassica cultivars must be tested for distinctness from previous cultivars, and uniformity and stability in their traits (Lombard et al., 2000); this is done to protect the breeder's intellectual property and verify that new cultivars have sufficient difference from previously existing ones (Lombard et al., 2000). This testing is tedious, expensive, and the results are sensitive to environmental influences on trait expression.
In the USA, numerous parameters are listed in the USDA's Plant Variety Protection Application Form for Kentucky bluegrass (Bonos et al., 2000b). These parameters include morphological traits such as plant height, panicle height, flag leaf dimensions, subtending leaf dimensions, and rhizome spread and length. These traits have also been used to classify Kentucky bluegrass cultivars into either seven (Murphy et al., 1997; Bonos et al., 2000b) or 12 types (Bonos et al., 2000a). The different morphological types are listed in Table 1, with the exception of the "Julia" type. The two main differences between the grouping systems are that (i) the "Compact" type in the seven group system has been divided into "Compact," "Compact-America," and "Compact-Midnight" in the revised 12 group system, and (ii) three new types have been added ("CELA," "Julia," and "Shamrock"). The "Other" type contains cultivars with characteristics intermediate between two or more types (Bonos et al., 2000b), while the "Unknown" listing contains cultivars for which the type name could not be determined in this study. The morphological data used to classify cultivars into these types are supplemented by data on disease susceptibility, including stripe smut [caused by Ustilago striiformis (Westend.) Niessl.] and leaf spot [caused by Drechslera poae (Baudys.) Shoem.] (Murphy et al., 1997). The use of morphological and disease reaction traits to distinguish varieties is time-consuming, environmentally dependent, and is not always sufficient to reliably distinguish the new and existing varieties.
|
Molecular markers, particularly RAPDs, have recently been applied to turfgrass research (Chai and Sticklen, 1998). In Kentucky bluegrass, markers have been successfully applied to characterization of the USDA germplasm collection (Johnson et al., 2002), identifying cultivar succession in swards (Lickfeldt et al., 2002), and determining the genetic origin of non-identical progeny from facultatively apomictic parents (Huff and Bara, 1993). Further, the utility of RAPD markers in discriminating among multiple Kentucky bluegrass cultivars has been demonstrated (Huff, 2001).
Therefore, this evidence suggests RAPD markers will provide an objective, rapid method of developing unique genetic profiles of each genotype. The use of RAPD marker data, combined with the known morphologically based classification system (Bonos et al., 2000b), can assist in developing cultivar blends that actually contain different genotypes, while providing the desired mix of phenotypes. Because cultivars belonging to a given morphological trait-based type presumably have similar morphology, the morphological trait-based grouping system is useful for deciding which cultivars to include in blends. Therefore, it would be useful to know whether or not cultivars in a given morphological type are genetically related.
In addition, Kentucky bluegrass is facultatively apomictic (Pepin and Funk, 1971), meaning that some sexual seed production is possible and not all seeds produced will be genetically identical to the parent plant. The proportion of such apomictic seeds produced varies between cultivars (Porceddu et al., 2002). Thus, rapid molecular verification of apomixis (as opposed to progeny testing) will facilitate production of genetically stable turf material (seed or sod) not subject to changes in appearance over time due to changes in the genetic composition of the turf stand.
Furthermore, marker-derived genetic diversity analyses will facilitate the use of hybrid vigor in developing new cultivars. Applying molecular markers to determine genetic distinctiveness in new germplasm such as exotic plant introductions (PIs) or interspecific hybrids, especially if they are found to be genetically distant from currently available commercial cultivars, will be very useful in broadening the genetic base for breeders, providing opportunities for development of new, unique cultivars. Molecular markers applied in this way have significant economic potential.
This study used RAPD markers to study the genetic relationships of 123 Kentucky bluegrass genotypes. The specific objectives of this study were (i) to assess the amount of genetic divergence between Kentucky bluegrass commercial cultivars, plant introductions (PIs), experimental Kentucky bluegrass breeding lines, and hybrids between Kentucky bluegrass and Texas bluegrass, (ii) assess the correlation between trait-based classification and marker-based grouping, and (iii) determine the amount of intracultivar and seed source variability, which can help assess the degree of apomixis.
| MATERIALS AND METHODS |
|---|
|
|
|---|
RAPD
For each genotype, DNA was individually extracted from three different plants by grinding tissue in PEX (potassium ethyl xanthogenate) buffer with a Bio101 (Carlsbad, CA) Fastprep machine, which pulverizes tissue by rapidly vibrating a ceramic bead within a 2 mL screw-top microcentrifuge tube containing the tissue. The DNA was precipitated with 6:1 ethanol and 7.5 M ammonium acetate, treated with RNase, and finally precipitated with 10:1 ethanol and 3 M sodium acetate as described in Scheef et al. (2003).
The 12 RAPD primers used in this study (Operon primers A13, G19, M14, O15, P8, Y5, Y9, AE4, AE18, AF20, AG10, and AG14) consisted of 10 base-pair random sequences (Operon Technologies, Alameda, CA), chosen for ability to produce bright, polymorphic bands. Reaction mixtures, 10 µL total, were 50 mM TrisHCl, pH 8.5, 10 mM KCl, 2 mM MgCl2, 500 mg/mL bovine serum albumin (BSA), 0.01% (v/v) xylene cyanole, 1.5% (v/v) Ficoll 400, 20 ng plant DNA, and 0.6 U Taq polymerase (Scheef et al., 2003).
All PCR reactions were run in an MJ Research (Waltham, MA) PTC-100 Programmable Thermal Cycler. Reaction conditions consisted of one cycle of 91°C for 1 min, 42°C for 15 s, and 72°C for 1 min 10 s, followed by 38 cycles of 91°C for 15 s, 42°C for 15 s, and 72°C for 1 min 10 s, and then cooling to 4°C at the end of cycling.
Reaction products were electrophoresed on 1.5% (w/v) agarose gels stained with ethidium bromide and photographed with an instant camera under UV light. Eighty-five bright, reproducible bands were scored and confirmed by two different people for presence or absence.
Data Analysis
Genetic distances between genotypes were calculated as the complement to the simple matching coefficient (Gower, 1972), with the complement to the simple matching coefficient equaling 1 minus (simple matching coefficient). This is equal to the number of discordant bands between two genotypes divided by the sum of the number of discordant and concordant bands, as in Beebe et al. (1995). Genetic distances were calculated for all genotypes (each cultivar or seed source times three seedling replicates of each genotype). This was done to allow calculation of genetic distances between seedling replicates of cultivars or seed sources as well, by measuring genetic distance between the three replicates of each genotype.
These genetic distance matrices were then used to create a multidimensional scaling (MDS) plot using the Kruskal scaling option in Systat version 5.2 (Wilkinson et al., 1992). This type of plot graphically represents the matrices in the form of x and y coordinates (Wilkinson et al., 1992). The MDS coordinates are quantitative variables that contain weighted information from each of the 85 markers and, as such, hypotheses regarding among-group variability can be tested by analysis of variance techniques. The MDS procedure thus is very similar to principal components analysis (PCA), as each MDS coordinate explains a proportion of the total variance in the data. As more MDS coordinates are added, the proportion of variance explained decreases. An important difference from PCA is that MDS can often fit an appropriate model in fewer dimensions (Wilkinson et al., 1992). For example, Nienhuis et al. (1993) noted that MDS analysis provided a clearer separation of Brassica entries into two dimensions than did PCA.
Both MDS coordinates were analyzed by mixed models analysis, using the Statistical Analysis System, SAS (Littel et al., 1996). The repeated measures option of SAS was used to model heterogeneous variances among germplasm sources by a trial-and-error system and the Akaike's and Bayesian Information Criteria to choose the best model. Germplasm sources were considered to be a fixed effect and genotypes within sources were a random effect. Mixed models analyses were based on 15 germplasm sources: 12 morphological-trait-based cultivar types, plant introductions, experimental breeding lines, and interspecific hybrids. This was done to formalize the visual comparisons of the marker-based MDS groupings of the reported 12-group trait-based classification system (Bonos et al., 2000a), plant introductions, experimental breeding lines, and interspecific hybrids. Mixed models analysis based on means of MDS coordinates has been previously used by Casler et al. (2003) to detect sources of variation in MDS coordinate values among creeping bentgrass clones.
In addition, analysis of molecular variance (AMOVA, Schneider et al., 2000) was conducted by Arlequin version 2.000 software. This was done to determine the overall similarity of the genetic distance matrices obtained from the three seedling replicates, mainly as an error-checking measure.
| RESULTS |
|---|
|
|
|---|
|
When mixed models analysis using MDS-derived x and y coordinates of the 123 entries separated into each of the four above-mentioned germplasm sources was conducted, statistical significance was detected (p < 0.0001 for both coordinates). Looking at the group means separately, the mean of the PIs appears different by both coordinates from the mean of the cultivars (Table 2). An unequal variance t test using SAS shows that the means of both the x and y coordinates are significantly different between the cultivars and PIs (p < 0.0001 for both coordinates). Also, the mean of the interspecific hybrids was different from the mean of the cultivars by their y coordinates (p < 0.0001) but not by their x coordinates.
|
|
|
As an example of this variability, Fig. 3 shows the RAPD profile based on Operon primer AF20 of the three replicate seedlings of cultivars Nassau and Chateau. Replicates 2 and 3 of cultivar Nassau had nearly identical banding patterns, while the two bands marked with arrows were amplified only in replicate 1. Thus seedling replicate 1 does not appear to be apomictically derived. On the other hand, all three replicates of cultivar Chateau have the same amplification pattern. Further, the amount of variability was not constant over the four major gene pools in this study (Table 4). The PIs had the highest level of variability, with the cultivars and interspecific hybrids intermediate, and the experimental breeding lines with the least within seed source variability.
|
|
| DISCUSSION |
|---|
|
|
|---|
Within the cultivars themselves, the distribution over the MDS plot was not uniform. Two groups of cultivars were noted, which upon comparison with the morphological trait-based classification systems were found to correspond with the BVMG type and the Compact-Midnight type. The BVMG type is interesting since it seems to form two groups in the MDS plot (Fig. 2), and as such little statistical significance was found between this type and others. The Compact-Midnight type, on the other hand, showed significant differences in the mean MDS y coordinate with several of the other morphological types. Since the y coordinate explains 37.2% of the genetic distance variance, this provides statistical evidence to support the grouping pattern seen in the MDS plot.
This correspondence suggests that most cultivars falling into these two groups are likely to be very similar genetically. No other correspondences between morphologically based cultivar groupings and molecular marker genetic distances were noted, except possibly for the Common type, which would suggest that the cultivars within these other groups are not necessarily genetically related. This way, they can be used for blended plantings with similar turf characteristics while still retaining genetic diversity within the planting.
These results generally agree with those of Huff (2001), in which the genetic relationship of 93 Kentucky bluegrass entries was studied using RAPDs. In that study, four out of five Compact-Midnight members were in the same dendrogram clade, while the fifth one, Explorer, clustered near the BVMG members as it did in this study. Further, all four BVMG cultivars clustered together, following the pattern of this study, and six of the eleven experimental lines tested in that study clustered in the same clade as BVMG, again similar to the results of this study, even though different experimental lines were used. Thus, evidence for this partial correspondence between RAPD-based grouping and morphological trait-based classification of cultivars is provided by two different studies. Further, since some of the cultivar groups have very few entries in this study, it is possible more correspondence between the two grouping methods may be found if a larger sample size is tested in the future.
When the three seedling replicates of each genotype were compared against each other to calculate their genetic distance, a rather high amount of genetic variability within cultivar or seed source was found for some of them, with values as high as 0.5. The mean within-source distance for all the genotypes was 0.22, or an average of 22% of RAPD bands being different. This value appears high for a facultatively apomictic species, although the frequency of apomixis can vary from close to zero to nearly 100% (Porceddu et al., 2002). Further, occasional mitotic aberrations in meristems can change chromosome numbers within a plant (Porceddu et al., 2002), and presumably other genetic characters, and these changes can then be perpetuated by apomixis. So, it is possible that these differences reflect the known variation in the amount of reproduction by apomixis in this species.
Also, this estimate is likely to be biased upward slightly because of possible low frequency errors in scoring or misamplification of some of the RAPD bands. Skroch and Nienhuis (1995) reported a scoring error value of 2%, so a similar scoring error value can probably be expected in this study as well. Additionally, if fewer than 85 markers had been included in the analysis, it is possible that this genetic distance estimate would have been lower, because of possible undersampling of genetic differences between the three replicates (Skroch et al., 1992). However this seems unlikely since sampling efficiency, as measured by reduction in coefficient of variation, is near its maximum at a marker number of 85 (Skroch et al., 1992).
It should be noted that the sample size of three replicates for each genotype may be low. For example, Huff (2001) used 96 seedling replicates of two cultivars, Baron and Unique, to estimate the number of off-type plants in these two cultivars. However, the present study used three seedling replicates of all 123 cultivars and seed sources, so that comparisons of RAPD variability among replicates of different germplasm sources were possible.
Examination of Fig. 3, however, would suggest technical factors likely do not explain all of the observed variation. In this figure the second and third replicate of cultivar Nassau and all three replicates of cultivar Chateau appear identical as would be expected for apomictic reproduction. However, it appears that replicate 1 of cultivar Nassau is an off-type. It is also possible that this difference is caused by seed contamination, but these cultivars are well spaced from each other as well as from the rest of the cultivars in the MDS plot, which is based on all 85 RAPD markers. Another possibility is that it is the result of a sexually produced seed, either cross- or self-fertilized (Huff and Bara, 1993).
In contrast, there is also evidence for genetic uniformity of seeds of a given cultivar between seedlots. Specifically, a fourth seedling of cultivar Award, from a different seedlot than the other three (listed in Table 1 as Award-B), was included in this study. Compared with the three replicates of Award, the RAPD banding patterns were 3.5, 8.2, and 8.2% different. Thus, even though some cultivars had variability between seedling replicates, this cultivar seemed fairly uniform even between seedlots.
Although further experimentation would be needed to determine the origins of the off-types, their frequency may increase as the age of the seed production field increases, due to a gradual buildup of an initially low proportion of sexually produced plants, which may or may not be derived from different cultivars. The age of the seed harvest fields from where the seeds originated in this study is not known. However, the experimental breeding lines had a lower value for within seed source variability than the cultivars, and since these genotypes are still in development rather than in commercial production, it is likely that the seed fields are kept more pure than for the cultivars. The PIs have the highest value for within seed source variability, suggesting a lower level of apomixis in these genotypes.
The information presented in this study provides support for the utility of RAPD marker-based genetic relationships for turfgrass genetics. Furthermore, the genetic information will be a valuable resource for turf breeders and managers, particularly in recommending cultivars to include in blends, choosing new, exotic genotypes to incorporate into breeding programs, and possibly providing a preliminary indication of variable levels of among-seedling genetic variability over different germplasm sources in Kentucky bluegrass.
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
Received for publication February 23, 2003.
| REFERENCES |
|---|
|
|
|---|
Related articles in Crop Science:
This article has been cited by other articles:
![]() |
T. D. Eaton, J. Curley, R. C. Williamson, and G. Jung Determination of the Level of Variation in Polyploidy among Kentucky Bluegrass Cultivars by Means of Flow Cytometry Crop Sci., November 1, 2004; 44(6): 2168 - 2174. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||