Crop Science
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online 1 January 2005
Published in Crop Sci 45:27-35 (2005)
© 2005 Crop Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Related articles in Crop Science
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Diaby, M.
Right arrow Articles by Casler, M. D.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Diaby, M.
Right arrow Articles by Casler, M. D.
Agricola
Right arrow Articles by Diaby, M.
Right arrow Articles by Casler, M. D.
Related Collections
Right arrow Other Forage Crops
Right arrow Crop Genetics

RAPD Marker Variation among Divergent Selections for Fiber Concentration in Smooth Bromegrass

M. Diabya and M. D. Caslerb,*

a Dep. of Agronomy, Univ. of Wisconsin-Madison, WI 53706-1597 USA
b USDA-ARS, U.S. Dairy Forage Research Center, 1925 Linden Dr. West, Madison, WI 53706-1108

* Corresponding author (mdcasler{at}wisc.edu).


    ABSTRACT
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Neutral detergent fiber (NDF) is considered the laboratory measure most closely correlated with voluntary intake of forages by ruminant livestock. The objectives of this study were to create smooth bromegrass (Bromus inermis Leyss.) populations divergent for NDF concentration in four smooth bromegrass germplasm pools and to identify random amplified polymorphic DNA (RAPD) marker changes in these divergently selected populations. Two cycles of divergent phenotypic selection led to significant linear responses in NDF among selection cycles for all germplasm pools. Within-population variation for RAPD markers was large, reflecting the outcrossing reproduction and the complex inheritance of smooth bromegrass. However, analysis of molecular variance revealed significant genetic differences among selected populations. Analysis of genetic distances showed that each cycle of selection created additional divergence, both among cycles and among germplasm pools. Up to 15 RAPD markers were associated with selection in each population, but only one marker was consistently associated with selection for NDF across all four germplasm pools (linear, homogeneous slopes, no drift). Some of the RAPD markers appear to have utility for marker selection or marker-assisted selection (MAS) to modify NDF concentration.

Abbreviations: IVDMD, in vitro dry matter digestibility • MAS, marker-assisted selection • NDF, neutral detergent fiber • PCR, polymerase chain reaction • QTL, quantitative trait loci • RAPD, random amplified polymorphic DNA


    INTRODUCTION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
DIGESTIBILITY, feed consumption, and energetic efficiency are the three components of forage nutritional value theoretically and/or empirically related to animal performance (Raymond, 1969). Digestibility, predicted simply and effectively by standardized in vitro laboratory techniques, has received the most attention by forage grass breeders. Quantitative genetic changes in in vitro dry matter digestibility (IVDMD) are generally due to changes in chemical composition of individual plant parts (Buxton and Casler, 1993) and only rarely due to changes in plant morphology (Kephart et al., 1989) or relative maturity (Buxton and Casler, 1993).

A forage sample consists of a cell-contents fraction soluble in neutral detergent and an insoluble cell-wall fraction (neutral detergent fiber, NDF), largely containing cellulose, hemicellulose, pectin, and lignin. The concentration of NDF is the single most effective laboratory predictor of an animal's ability to consume a forage ad libitum (Van Soest, 1994). However, it is not clear whether the mechanism of intake is regulated by negative feedback from ruminal tract distension due to the high bulk density of high-fiber forages (physiological regulation) or hormonal signals indicating satiety, a fulfilling of current nutritional requirements (metabolic regulation) (Van Soest, 1994). The concentration of NDF is the most common selection criterion to improve the intake potential of forage crops. Also, preference or palatability has been found to be associated with low fiber concentration and high digestibility (Falkner and Casler, 1998; Gangstad, 1964).

Genetic studies and germplasm research through intensive selection and breeding efforts have provided a solid scientific basis for improving forage quality of smooth bromegrass. However, selection designs have been inefficient in estimating independently the effects of some nutrients or constituents on forage nutritional value or separating negative associations between some forage quality traits and forage yield or disease resistance (Casler, 2001; Casler and Vogel, 1999). Therefore, a long-term divergent selection program for NDF was undertaken initially to create more variability in smooth bromegrass populations, and to establish relatively unconfounded populations differing in NDF concentration capable of addressing some important genetic and breeding questions.

Over time and selection cycles, population phenotypes change as a result of changes in the frequencies of favorable alleles. Allele frequencies may change as a result of selection or drift. Selection responses will occur for alleles under direct selection pressure and those within their linkage blocks. Although most molecular markers are considered to be selectively neutral, changes in marker frequencies associated with changes in population performance have been reported (Keithley and Bulfield, 1993; Ollivier et al., 1997; Stuber and Moll, 1972; Stuber et al., 1980). Consistent frequency changes in molecular markers, without effects of drift, may be a mechanism to identity linkages between markers and quantitative trait loci (QTL), potentially improving selection efficiency.

The objectives of this study were (i) to create divergent populations for NDF concentration and (ii) to identify RAPD marker changes during two cycles of divergent selection for NDF concentration, for eventual application of MAS to improve the efficiency and cost-effectiveness of phenotypic selection for reduced NDF concentration.


    MATERIALS AND METHODS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Germplasm
Divergent selection for NDF concentration was conducted in four germplasm pools representing three different levels of germplasm improvement: the cultivars Alpha and Lincoln and the two experimental germplasm pools WB19e and WB88S. The four germplasm pools are unrelated to each other, with the exception that Alpha and a population selected for high digestibility out of Lincoln were two of the four parents of WB19e. Lincoln and WB88S were derived directly from unrelated natural populations. The original parentage of Alpha and WB19e did not contain any germplasm related to WB88S. WB19e is a Syn2 generation strain cross among four high-IVDMD populations from Wisconsin (B8HD and ‘Badger’) and Nebraska (NE-BI-1 and NE-L-D), contributed by K.P. Vogel (USDA-ARS), and WB88S-Alt is a Syn1 strain cross among five accessions collected in the Altai Mountains of south-central Russia on a 1988 expedition (USDA-ARS, 1990).

Selection for Neutral Detergent Fiber—Cycle 1
Cycle 1 was completed by Casler (2002). Seed of the four germplasm pools was planted in early spring 1991 in the greenhouse. Three hundred 70-d-old seedlings for each population were transplanted to the field in a spaced-plant nursery on 1.0-m centers on a Plano silt loam (fine-silty, mixed, mesic Typic Argiudoll) in May 1992 at Arlington Agricultural Research Station, WI (Casler, 2002). Four germplasm pools were planted in separate selection nurseries with 4 m between germplasm pools. Plants within each population were arranged in ten blocks of 30 plants each, with a spacing of 0.9 m between all adjacent plants. Weeds were controlled by pre-emergence herbicide applications (Falkner and Casler, 1998). Plants were fertilized with 112 kg N ha–1 in late June.

A leaf tissue sample was clipped from each plant at a vegetative growth in mid-August 1992, with a stubble height of 10 cm. Tissue samples were placed in paper bags and dried at 55°C. Dried samples were ground through a 1-mm screen of Wiley-type mill and reground through a 1-mm screen of cyclone mill. Two independent subsets of each sample were scanned on a near-infrared reflectance spectrophotometer (NIRS). Two random plants from each block of each population comprised a stratified random subset of 80 plants that was subjected to wet-laboratory analysis. Concentration of NDF was determined using the procedure of Van Soest et al. (1991), omitting the {alpha}-amylase step. Data on NDF of the 80-plant subset were used to calibrate the NIRS for prediction of the entire set of 2400 scanned samples (four germplasm pools x 300 plants x two scanned subsets). Means over the two scanned subsets of each sample were computed before selection. Some field and laboratory variability was removed from the estimates of plant phenotypes by processing samples in numerical order and computing t scores to adjust for differences among block means (Casler, 1992): tj = (Xij Mj)/sj, where tj is the t score for NDF of the ijth plant, Xij the raw datum for the ijth plant, Mj the mean of the jth block, and sj the standard deviation of the jth block.

Ten plants within each population were selected in spring 1993 from the selection nursery as the most extreme plants in two categories: high NDF and low NDF. Each group of ten selected plants were transplanted into a four-replicate crossing block using 60 cm plant spacing. During an 8- to 12-d period of anthesis, pollen from plants was wind-dispersed among selected plants in each crossing block. The crossing blocks were isolated at the Arlington Agricultural Research Station by at least 100 m from each other and from other smooth bromegrass to avoid contamination. The crossing blocks were fertilized in spring at the rate of 56 kg N ha–1.

One cycle of divergent selection for NDF concentration in four germplasm pools (WB19e, Alpha, WB88S, and Lincoln) was completed with the harvest of seed from individual clones in each crossing block in August 1993 (Casler, 2002). Seed for each clone was individually harvested, dried, cleaned, and weighed. Balanced bulks of seed were made for each population. Selected germplasm pools were identified as C-1 (low NDF) and C+1 (high NDF).

Selection for Neutral Detergent Fiber—Cycle 2
Seed from the C-1 and C+1 cycles of all four germplasm pools generated in the first cycle, was germinated in January 1996 in the greenhouse. Three hundred fifty seedlings from each population were transplanted to the field at Arlington on a Plano silt loam (fine-silty, mixed, mesic Typic Argiudoll) in May 1996 in a spaced-plant nursery on 1.0-m centers. Eight populations were planted in separate selection nurseries with 4m between populations, and plants within each population were arranged in ten blocks of 35 plants each, with a spacing of 1 m between all adjacent plants. Plants were mowed twice during the establishment year. Weeds were controlled with herbicide (Falkner and Casler, 1998). Plants were fertilized with 112 kg N ha–1 in late June and sampled for laboratory analysis.

Neutral detergent fiber was determined on a random sample of approximately 20 tillers clipped from each of the 350 plants of each population in May 1997 at a 9-cm stubble height. Plants were approximately 20 to 30 cm tall and consisted entirely of leaves. Thirty-five plants were harvested from each block and always processed together throughout the drying, grinding, and laboratory processes. Plant samples were placed in perforated paper bags, dried at 55°C, ground through a 1-mm screen in a Wiley-type mill, reground through a 1-mm screen in a cyclone mill to improve their particle size uniformity, and scanned by NIRS. An 80-plant calibration set was chosen, analyzed for NDF concentration, and used to predict NDF for the remaining samples as described for Cycle 1. The only change from Cycle 1 was the method of choosing the calibration set, which was based on a cluster analysis of the samples, using their individual wavelength reflectance data (Shenk and Westerhaus, 1991).

Ten plants within each population were selected in spring 1998 from the selection nursery as the most extreme plants in two categories: high NDF and low NDF. Each group of 10 selected plants were transplanted into a four-replicate crossing block using 60 cm plant spacing. The crossing blocks were isolated at the Arlington Agricultural Research Station by at least 100 m from each other and from other smooth bromegrass to avoid contamination. The crossing blocks were fertilized in spring at the rate of 56 kg N ha–1.

Two cycles of divergent selection for NDF concentration in four germplasm pools (WB19e, Alpha, WB88S, and Lincoln) were completed with the harvest of seed for the second cycle of selection in July 1999. Seed for each clone was individually harvested, dried, cleaned, and weighed. Balanced bulks of seeds were made for each population. Selected populations were identified as C-2 (low NDF) and C+2 (high NDF).

Field Experiment
Seedlings of 20 populations (C-2, C-1, C0, C+1, and C+2 for each of the four germplasm pools) were transplanted in the field in May 2000. The field experiment was established at Arlington, WI, [43°20' N, 89°23' W; Plano silt loam (fine-silty, mixed, mesic Typic Argiudoll)] and Marshfield, WI, [44°40' N, 91°53' W; Withee silt loam (fine-loamy, mixed, frigid, Aeric Glossoboralf)]. The experimental design was a randomized complete block with four replicates at each location. Plot size was 1 x 3 m with 30 plants per plot and 30-cm centers. Weeds were controlled with herbicide (Falkner and Casler, 1998). Field plots at each location were fertilized with 112 kg N ha–1.

In 2001, field plots at Arlington and Marshfield were fertilized with 112 kg N ha–1 in April and June. Forage samples were hand-clipped at a 9-cm stubble height in May (immediately after heading) and August (vegetative growth stage, with 3 to 4 nodes per sterile culm), by taking a small sample of tillers from within each of thirty 0.3 x 0.3-m grids. Samples were dried at 55°C and ground to pass through a 1.0-mm screen in a wiley-type mill and a cyclone mill. Samples were analyzed for NDF concentration by NIRS as described for the Cycle-2 selection phase (80 samples in the calibration set).

DNA Isolation
Seed from the C-2, C-1, C0, C+1, and C+2 cycles divergently selected for NDF from all four germplasm pools was germinated in the greenhouse in 1999. A random sample of seedlings from the 20 populations were maintained in the greenhouse for RAPD marker analysis. The RAPD marker data and NDF data were collected on different random samples of plants from the 20 populations because we felt it was important to save the RAPD-genotyped plants for use in future studies. Because of the rhizomatous nature of smooth bromegrass, this would have been impossible with plants transplanted to the field.

DNA was extracted from 14 to 25 individual seedlings from each population as described by Skroch and Nienhuis (1995). Approximately 0.5 to 0.75 g of fresh tissue was harvested and ground with 500 µL of potassium ethyl xanthogenate (PEX) (Sigma-Aldrich, St. Louis, MO) at maximum speed of 5.0 m s–1 for 40 s with the Bio 101 (Vista, CA) Savant FP120 Fast PrepTM. After grinding, tissue was transferred to centrifuge tubes and allowed to incubate for 30 min in a 65°C water bath. After organic and aqueous phases of the extraction mixture were separated by centrifugation (Eppendorf 5415C microfuge), nucleic acids were precipitated by adding a 6:1 mixture of 95% (v/v) ethanol and 7.5 M ammonium acetate. After removing RNA (with 100 µg/mL RNase A for 1 h at 37°C) and any remaining debris, DNA was reprecipitated by the addition of 10:1 solution of ethanol and 3 M sodium acetate. After a 70% (v/v) ethanol wash and pelleting, DNA was hydrated in TE buffer (1 mM Tris, pH = 8.0, 0.1 mM EDTA, pH = 8.0). DNA concentrations were quantified in a logical numerical order with a Hoefer Scientific TKO-100 Fluorometer (Amersham Pharmacia Biotech, Piscataway, NJ).

RAPD Reactions
RAPD reactions were performed as described in Johns et al. (1997) in an M.J. Research, Inc. (Wattham, MA) PTC-100 Programmable Thermal Controller. Cycling temperature settings were 91°C for denaturation, 42°C for annealing, and 72°C for elongation. In the first cycle, cycling times were 60 s for denaturation, 15 s for annealing, and 70 s for elongation. During the subsequent 39 cycles, denaturation was set for 15 s, annealing for 15 s, and elongation for 70 s.

Polymerase chain reaction (PCR) amplifications were performed in a final reaction volume of 10 µL, containing the following reaction buffer: 50 mM Tris, pH 8.5, 20 mM KCl, 2 mM MgCl2, 500 µg/mL of bovine serum albumin (BSA), 2.5% (w/v) ficoll 400, and 0.02% (w/v) xylene cyanol. Reactant concentrations were 100 µM dNTPs (deoxy nucleotide triphosphates) (Promega, Madison, WI), 2 ng/µL of DNA template, 0.4 µM of decamer primer Operon Technologies, Inc. (Alameda, CA) and University of British Columbia, (Vancouver, BC, Canada), and 0.6 unit (5 units/µL) of Taq DNA polymerase (Promega, Madison, WI). All RAPD reaction products were electrophoresed in 20 x 25 cm, 1.5% (w/v) agarose gels in 1x TBE (Tris, Boric Acid, EDTA) buffer. Gels were run for 2 h at 300 V in Gibco/BRL Life Technologies (Invitrogen, Carlsbad, CA) H4 gel apparatus, stained with ethidium bromide, and illuminated by UV light and subsequently photographed with Polaroid 667 film.

Primer Screening
To enhance the probability of detecting a relationship between RAPD marker frequencies and NDF selection responses, primers were prescreened for a correlation between band frequencies and NDF concentration. Seventeen RAPD primers used by Diaby and Casler (2003) were evaluated by correlation analysis with NDF data from 27 smooth bromegrass populations, obtained from Casler et al., 2000). The four original populations included in this study were also included in the study of Casler et al. (2000). Seven primers had a relatively high frequency of polymorphic bands that had frequencies correlated with NDF concentration of the 27 populations (0.35 < |r| < 0.65). These primers (A09, A18, AE12, AF06, AF14, and AG14 from Operon Technologies Inc., Alameda, CA, and UBC318 from University of British Columbia, Vancouver, BC, Canada) were used on the 433 plants of this study, generating 83 polymorphic RAPD markers. The frequency of each marker (scored as "1" for presence of the band or "0" for absence of the band for each individual plant and stored as a binary matrix) was computed for each of the 20 populations.

Data Collection and Statistical Analysis
Field-plot NDF values over two harvests and two locations were analyzed by general linear models analysis of variance (SAS, 1999). Replicates and locations were assumed to be random effects, while populations, cycles, and harvests were fixed effects. Linear regression for mean divergent selection responses and contrasts for the linear effect of selection cycles (cycles-linear) were computed in the four germplasm pools. Significance of linear regressions were tested by contrasts in the analysis of variance.

From the 83 polymorphic RAPD markers transformed into a binary matrix, a pairwise Jaccard similarity coefficient matrix was computed by NTSYS-PC 2.01 (Rohlf, 1997) on all individuals across populations. The similarity matrix Sij (similarity between individual plants i and j) was converted to a Euclidean distance matrix by the elementwise formula: (1 – Sij)0.5. Euclidean distances, converted from Jaccard similarity coefficients, were used as the measure of genetic distance between all individuals.

Analysis of molecular variance (AMOVA; Excoffier et al., 1992; Schneider et al., 1997) was performed on all individuals, partitioning the Euclidean distance matrix into three sources of variation: among populations, among cycles within populations, and plants within populations and cycles. Variance components were estimated by equating AMOVA mean squares to their expectations. Variance components were tested by nonparametric permutation tests (Schneider et al., 1997). Cluster analysis, based on the unweighted pair-group method of arithmetic averages (UPGMA; SAS, 1999) was used to construct a distance dendrogram for the 20 populations.

Divergent selection for NDF utilized only additive genetic variation for NDF, and linear regressions of phenotype or marker frequencies across selection cycles are measures of additive effects (slopes) and additive genetic variation (variance due to regression) (Falconer and Mackey, 1996). Moreover, selection for NDF concentration was practiced in divergent directions, under the same effective population sizes and selection intensities. Therefore, the effects of drift and selection are largely independent across selection cycles. The effect of drift, measured as asymmetry of selection responses about the base population mean (Casler, 1999; Clayton et al., 1957; Falconer, 1953), is orthogonal to the linear effect of selection in this selection scheme. Nevertheless, to increase the rigor of the marker selection process, statistical procedures developed by Wilson (1980) were used to test whether fluctuations in frequencies from cycle to cycle could be attributed to genetic drift alone or whether divergent selection had occurred. Linear trends were tested by chi-square tests incorporating sampling variation due to restricted population size and restricted number of plants sampled for genotypic analysis. If significant deviations from random drift were detected, and significant directional changes as represented by linear trends were observed, then evidence would be judged sufficient to suggest that the particular marker locus may be linked to one or more loci affecting NDF.

Finally, RAPD marker regressions on cycle number were tested for homogeneity across germplasm pools, by contrasts within a generalized linear model (GENMOD) for binomial data (SAS, 1999). The LOGISTIC procedure (SAS, 1999) was used to fit a linear regression model for binary data to selected markers by maximum likelihood estimation.


    RESULTS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Consistent and significant linear responses were found for NDF concentration for all four populations divergently selected for NDF (Fig. 1) . Two cycles of divergent selection for NDF concentration led to consistent responses averaged over both harvests and locations. All four germplasm pools showed significant (P < 0.05) responses to selection, with responses ranging from 6.6 to 8.6 g kg–1 cycle–1 and R2 values ranging from 0.79 to 0.91.


Figure 1
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 1. Linear regressions of mean neutral detergent fiber (NDF) concentration on selection cycle for four smooth bromegrass populations divergently selected for NDF. Regressions statistics for Alpha, WB19e, Lincoln, and WB88S, respectively, were b = 7.2 ± 1.5, 7.0 ± 2.1, 6.6 ± 1.6, and 8.6 ± 1.4 g kg–1 cycle–1; R2 = 0.83, 0.79, 0.91, and 0.88.

 
Analysis of molecular variance (AMOVA) showed all groups of smooth bromegrass germplasm to have high within-population genetic variation, ranging from 94.2 to 95.0% of the total (Table 1). Intergroup variation, among the four smooth bromegrass germplasm pools, was low ranging from 0.6% for Alpha vs. WB19e to 3.7% for Alpha vs. WB88S, reflecting the breeding history of the different germplasm pools (Table 2). Alpha and a population selected from Lincoln were two of the four parents of WB19e, while WB88S is unrelated to the other three populations. The intergroup variation for WB19e vs. Alpha and WB19e vs. Lincoln was not significant, indicating that WB19e is very similar to both Alpha and Lincoln.


View this table:
[in this window]
[in a new window]

 
Table 1. Partitioning of total RAPD marker variation of pairwise comparisons into between-population and within-population components for populations divergently selected for neutral detergent fiber (NDF) concentration in the four germplasm pools of smooth bromegrass by analysis of molecular variance (AMOVA).

 

View this table:
[in this window]
[in a new window]

 
Table 2. Partitioning of total RAPD marker variation of pairwise comparisons among groups, among populations within groups and within-population components for contrasting smooth bromegrass germplasm groups by analysis of molecular variance (AMOVA).

 
The clustering of populations, based on their genetic distances, partially reflected selection history, but showed little correspondence with pedigree (Fig. 2) . WB88S should be the most distant of the four original populations, but it did not appear to be particularly unique. Consecutive selection cycles clustered together in several cases, and opposite selection cycles (high vs. low NDF) typically did not cluster together, demonstrating a relationship between genetic distances and selection divergence. This is illustrated in a general increase in genetic distances with each incremental selection event for each of the four germplasm pools (Fig. 3) . The linear regression of genetic distance on the numerical distance between selection cycles was significant (P < 0.0001), linear (R2 = 0.71), and homogeneous across germplasm pools (D = 0.004 + 0.024N, where D = genetic distance and N = numerical distance between paired cycles). Furthermore, genetic distances between the four germplasm pools (Alpha, WB19e, Lincoln, and WB88S) increased with selection in either direction, increasing the genetic diversity among germplasm pools for both high-NDF and low-NDF selections (Fig. 4) .


Figure 2
View larger version (15K):
[in this window]
[in a new window]

 
Fig. 2. Cluster dendrogram of 20 smooth bromegrass populations based on 83 RAPD marker band frequencies. Populations are abbreviated as follows: A = Alpha, B = WB19e, L = Lincoln, S = WB88S, 0 = original population, 1 = Cycle 1, 2 = Cycle 2, "+" = high NDF, and "-" = low NDF.

 

Figure 3
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 3. Scatterplot of genetic distances between paired smooth bromegrass populations that differ by one, two, three, or four selection cycles for divergent NDF concentration (enumerative difference between cycles). The linear regression equation was D = 0.004 + 0.024N, where N = enumerative difference between selection cycles, R2 = 0.71, P < 0.0001.

 

Figure 4
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 4. Scatterplot of genetic distances between paired smooth bromegrass populations within each of five divergent-NDF selection cycles, including the original populations (Cycle 0). The quadratic regression equation was D = 0.0579 – 0.0078C + 0.0064C2, R2 = 0.43, P = 0.0005.

 
The 83 polymorphic RAPD markers were divided into four categories on the basis of the significance of drift and selection effects (Table 3). Differential P values were used to discriminate markers for drift and selection because Type II errors are more serious for drift, whereas Type I errors are more serious for selection. On the basis of P < 0.01, 11 to 18% of the markers had significant linear selection responses, after accounting for the effects of drift. Most of these markers also had a significant drift effect. Significant drift, but nonsignificant selection was the largest category for each of the four germplasm pools, reflecting the importance of drift in determining changes in RAPD markers with selection.


View this table:
[in this window]
[in a new window]

 
Table 3. Number of RAPD markers, out of a total of 83, that were characterized by specific P values for the effect of random genetic drift or the effect of selection after accounting for drift.

 
The 83 polymorphic RAPD markers were also divided into four categories by generalized linear model and logistic regression across divergent selection cycles for NDF concentration in the four smooth bromegrass germplasm pools (data not shown). Fifteen markers had a significant linear main effect for cycles and a cycle-linear x germplasm pool interaction (P < 0.05). Thirty markers had nonsignificant linear main effect for cycles and significant cycle-linear x germplasm pool interaction (P < 0.05). For another 30 markers, neither the linear main effect for cycles or the cycle-linear x germplasm pool interaction were significant. The remaining eight RAPD bands had significant linear main effect for cycles and no cycle-linear x germplasm pool interaction (Table 4), suggesting a consistent association between those specific markers and potential QTL for NDF in the four smooth bromegrass germplasm pools. These eight markers were derived from only three primers.


View this table:
[in this window]
[in a new window]

 
Table 4. P values for polymorphic RAPD bands demonstrating high discriminating power between divergent selection cycles for NDF concentration of four smooth bromegrass populations.

 
Among the eight polymorphic bands found to have shifted in frequency through selection cycles, five tended to have negative slope and three (AG14.1300, AG14.0825, and UBC318.0325) tended to have positive slope across selection cycles (Table 4). Despite the lack of cycle-linear x germplasm pool interaction for these eight markers, logistic regressions were often inconsistent across populations. In some cases (e.g., AE12.0975 and AE12. 0650), regression coefficients were sufficiently small that random variability caused changes in the sign of the regression coefficient. Other markers had sufficiently large responses that variation among populations was relatively unimportant.

At least one germplasm pool had a significant drift effect for seven of the eight markers (Table 5). Often the significance of drift caused the chi-square test for linear selection response to be nonsignificant. Only one marker, AG14.0825, had nonsignificant drift effects (P > 0.10) and significant linear selection response (P < 0.05) for all four germplasm pools. This marker also had the most consistent and linear logistic regressions. The linearity, consistency, and goodness-of-fit of these selection responses is demonstrated in Fig. 5 .


View this table:
[in this window]
[in a new window]

 
Table 5. P values for random genetic drift and linear selection response and logistic regressions of eight polymorphic RAPD bands demonstrating statistically similar slope as a function of cycle and direction of selection for neutral detergent fiber (NDF) concentration in four smooth bromegrass germplasm pools.

 

Figure 5
View larger version (23K):
[in this window]
[in a new window]

 
Fig. 5. Logistic regressions of RAPD marker band AG14.0825 frequencies (Pij) as a function of cycle and direction of selection for neutral detergent fiber (NDF) concentration in four smooth bromegrass populations. The common-slope logistic regression equation was LN[Pij/(1 – Pij)] = 1.8603 – 0.4487di1 – 1.0557di2 – 0.6267di3 + 0.4620i, where dij = 0 (absence) or 1 (presence) for each population; j = 1 for Alpha, 2 for WB19e, or 3 for Lincoln; and i = cycle number.

 

    DISCUSSION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Significant linear selection responses for all four smooth bromegrass germplasm pools divergently selected for NDF support previous studies reporting genetic progress for NDF concentration and other forage quality traits in smooth bromegrass (Casler, 1999; Casler et al., 2000) and other forage species (Surprenant et al., 1988; Wolf et al., 1993). These results confirm the heritable nature of NDF concentration and that genetic changes in NDF concentration are relatively insensitive to genotype x environment interactions. The latter was indicated by the fact that selection responses were consistent at the two evaluation locations and that selection was evaluated in different environments from the original selection nurseries.

The large within-population variation observed for all germplasm groups reflected the outcrossing mode of reproduction and probably the complex genome organization of smooth bromegrass (Armstrong, 1979). The variation partitioned between and within populations is dependent on the material under study and on the breeding system of the species. The patterns of variation observed in this study were similar to those found in several outcrossing grass species, including smooth and meadow bromegrass (Bromus riparius Rehm.) (Ferdinandez et al., 2001), blue grama (Bouteloua gracilis [H.B.K.] Lag. ex. Steud.) (Phan, 2000), buffalograss (Buchloë dactyloides [Nutt.] Engelm) (Peakall et al., 1995), and perennial ryegrass (Lolium perenne L.) (Huff, 1997) where within-population variation was much higher than between-population variation.

Despite the large amount of within-population variation, variation among cycles was significant, indicating genetic differentiation among populations divergently selected for NDF for each smooth bromegrass germplasm pool. The high genetic similarity of WB19e with Alpha and Lincoln as revealed by AMOVA was not surprising, because Alpha and Lincoln represent half of the pedigree of WB19e. The relatively large intergroup variation values for all the WB88S comparisons confirm the unique nature of this germplasm pool relative to the other three germplasm pools in this study. WB88S is a wild germplasm and represents the meadow climatype of smooth bromegrass, while the other three germplasms are cultivated forms of the steppe climatype of smooth bromegrass. The meadow and steppe climatypes of smooth bromegrass are phenotypically and genetically distinct from each other (Casler et al., 2000; Diaby and Casler, 2003).

Genetic changes due to selection were confirmed by genetic distance analyses. The tendency of high-NDF or low-NDF populations to cluster together, particularly Cycle-1 and Cycle-2 populations within common germplasm pools, suggested that rapid changes occurred as a result of selection and/or drift. Furthermore, these pairs of consecutive cycles within the cluster dendrogram suggested that some RAPD markers were linked to QTL for NDF and that changes in RAPD marker frequency, in some cases, were a direct result of selection on linkage blocks. Each additional cycle of selection resulted in an average 0.024 increase in genetic distance between pairs of cycles. Conversely, the relatively large genetic distances between germplasm pools, regardless of the cycle of selection, suggested that these linkage blocks were rarely collinear across germplasm pools. Furthermore, the increase in genetic distances in both the high-NDF and low-NDF directions suggested that the four germplasm pools were more divergent after selection than before selection. Differential linkage disequilibria across the four germplasm pools is the most likely explanation for this observation. Selection appears to have acted on different markers in each germplasm pool, largely based on differential linkage relationships with QTL for NDF.

Development of reliable linkage maps in most forage crops (complex polyploids with polysomic inheritance) is more complicated than in simple diploids. Furthermore, because of differences in linkage relationships and allelic composition, linkage maps have limited usefulness across different populations and germplasms (Lübberstedt et al., 1997; 1998). Identification of individual markers associated with high-NDF vs. low-NDF plants within an array of germplasms could facilitate MAS in highly heterogeneous populations of cross-pollinated smooth bromegrass, without the need for a linkage map.

Phan (2000) found four RAPD bands from four different primers that distinguish between original, ecovar, and cultivar populations of blue grama on the basis of differences in the band frequencies among the assayed populations. In the present study, among the eight most discriminating bands across selection cycles, only one marker, AG14.0825, had R2 values sufficiently high for all four germplasm pools that it might serve as a selectable marker for NDF concentration in all four germplasms.

This study confirms the potentially strong effect of drift due to restricted effective population size (Kahler, 1983; Guse et al., 1988; Falconer, 1953). However, linear responses for marker AG14.0825, after accounting for drift, were significant across all four germplasm pools of smooth bromegrass. These significant linear trends observed across all germplasm pools through cycles could be attributed to divergent selection for NDF concentration. Furthermore, the linear selection response, due to change in allele frequency of this marker, accounted for 86 to 96% of the variation among cycles.

The remaining markers showed some selection responses, but were inconsistent or had weak linear trends across germplasm pools, implying that they were specific for a subset of germplasm pools or there may be unknown factors accounting for the variability observed for some markers. The majority of markers showing inconsistency in their linear trend across all germplasm pools suggest the presence of variable levels of association between markers and QTL for NDF in the four germplasm pools of smooth bromegrass. Germplasm pools may have differed in the number of linkage blocks dragged from one cycle to another during selection for NDF concentration. Chromosome blocks containing QTL for NDF were variable among germplasm pool genomes, reflecting probably different forms of alleles in the complex biochemical pathways that contribute to NDF accumulation.

In a randomly mating population, only polymorphisms with extremely tight linkage to a locus with phenotypic effects are likely to demonstrate significant marker-QTL associations (Falconer and Mackay, 1996; Weir, 1996; Labate et al., 2000). In this study, the logistic regression revealed evidence for different levels of association between markers and potential QTL for NDF concentration across four germplasm pools of smooth bromegrass. Neutral detergent fiber appears to be under the control of many effective factors in the smooth bromegrass genome. The NDF fraction consists of several compounds—lignin, cellulose, hemicellulose, proteins, pectins, and minerals—with complex chemical composition and biosynthetic pathways (Casler, 2001; Boudet et al., 1995). Changes in most of these compounds could lead to changes in NDF concentration. Moreover, the germplasm pools of smooth bromegrass used in the divergent selection experiment have different genetic backgrounds, contributing to variable linkage disequilibrium states and variable allele frequencies. Identification of QTL for ADF (acid detergent fiber, closely related to NDF) was strongly dependent on the tester used to create maize (Zea mays L.) testcrosses, showing almost no correspondence in QTL across testers (Lübberstedt et al., 1997).

This is the first report of association of a quantitative character with a molecular marker gene locus in smooth bromegrass. A statistically significant association of QTL for NDF with the presence of a marker indicates either that the marker locus resides within a QTL for NDF (pleiotropy) or that the marker locus is linked with one or more QTL for NDF. The strength of the association of a marker with a QTL depends on the degree of linkage and of the size of the effects of the QTL alleles. The consistency of linear responses to marker AG14.0825 among germplasm pools and the lack of drift effects for this marker suggest that it may be pleiotropic with a QTL for NDF concentration, although this is relatively improbable.

In conclusion, we identified several dominant PCR-based markers potentially detecting QTL for NDF concentration in selected populations of smooth bromegrass. Although RAPD markers can be applied in high throughput analyses required by most breeding programs, they are dominant markers and their reaction sensitivity may be an impediment where methodological conditions cannot be repeated or when the genotype scoring process is subject to error. The conversion of RAPD bands to sequence characterized amplified repeats (SCAR) markers will provide more specific and stable markers. Using a combination of SCAR markers, developed from this study, and codominant markers to identify linkage relationships, it should be possible to create interval maps of putative QTL for NDF within each germplasm pool and to utilize these QTL in a marker-selection or MAS program.


    ACKNOWLEDGMENTS
 
The authors would like to express their thanks to Drs. A.T. Phan, J.G. Coors, and M.J. Havey for helpful suggestions; M.E. Sass; R. Kethireddypally; P.M. Crump and T.J. Tabone for providing technical support.


    NOTES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Mention of a trademarked or proprietary product does not constitute a guarantee or warranty by the USDA-ARS and does not imply its approval over other suitable products.

Received for publication October 8, 2003.


    REFERENCES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 


Related articles in Crop Science:

THIS ISSUE IN CROP SCIENCE

Crop Science 2005 45: xi. [Full Text]  



This article has been cited by other articles:


Home page
Crop Sci.Home page
P. Sulima, J. A. Przyborowski, and D. Zaluski
RAPD Markers Reveal Genetic Diversity in Salix purpurea L.
Crop Sci., May 11, 2009; 49(3): 857 - 863.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
M. D. Casler and M. Diaby
Positive Genetic Correlation between Forage Yield and Fiber of Smooth Bromegrass
Crop Sci., November 24, 2008; 48(6): 2153 - 2158.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
M. D. Casler, C. A. Stendal, L. Kapich, and K. P. Vogel
Genetic Diversity, Plant Adaptation Regions, and Gene Pools for Switchgrass
Crop Sci., November 7, 2007; 47(6): 2261 - 2273.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
C. Stendal and M. D. Casler
Maximizing Efficiency of Recurrent Phenotypic Selection for Neutral Detergent Fiber Concentration in Smooth Bromegrass
Crop Sci., January 24, 2006; 46(1): 297 - 302.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
C. Stendal, M. D. Casler, and G. Jung
Marker-Assisted Selection for Neutral Detergent Fiber in Smooth Bromegrass
Crop Sci., January 24, 2006; 46(1): 303 - 311.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Related articles in Crop Science
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (3)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Diaby, M.
Right arrow Articles by Casler, M. D.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Diaby, M.
Right arrow Articles by Casler, M. D.
Agricola
Right arrow Articles by Diaby, M.
Right arrow Articles by Casler, M. D.
Related Collections
Right arrow Other Forage Crops
Right arrow Crop Genetics


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