Crop Science Grow Your Career with CSSA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


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 Similar articles in this journal
Right arrow Similar articles in ISI 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 ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by DeHaan, L. R.
Right arrow Articles by Wyse, D. L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by DeHaan, L. R.
Right arrow Articles by Wyse, D. L.
Agricola
Right arrow Articles by DeHaan, L. R.
Right arrow Articles by Wyse, D. L.
Related Collections
Right arrow Other Forage Crops
Right arrow Crop Genetics
Crop Science 43:402-408 (2003)
© 2003 Crop Science Society of America

PLANT GENETIC RESOURCES

Illinois Bundleflower Genetic Diversity Determined by AFLP Analysis

L. R. DeHaan*,a, N. J. Ehlkeb, C. C. Sheafferb, G. J. Muehlbauerb and D. L. Wyseb

a The Land Institute, 2440 E. Water Well Rd., Salina, KS 67401
b Dep. of Agronomy and Plant Genetics, Univ. of Minnesota, 411 Borlaug Hall, 1991 Upper Buford Circle, St. Paul, MN 55108

* Corresponding Author (dehaan{at}landinstitute.org)


    ABSTRACT
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Illinois bundleflower [Desmanthus illinoensis (Michx.) MacMillan] is an herbaceous perennial legume native to North America. Useful as a N2–fixing plant in warm-season grass pastures, Illinois bundleflower is also a promising perennial grain crop. Knowledge of the distribution of genetic variation in Illinois bundleflower would increase efficiency of germplasm preservation and expedite plant breeding progress. The objective of this experiment was to determine the distribution of genetic variation within and between Illinois bundleflower accessions by amplified fragment length polymorphism (AFLP) markers. Semi-automated fluorescence-based AFLP analysis was performed on three individuals from each of 50 accessions. We identified 222 markers, 159 of which were polymorphic. Within-accession diversity was low (HS = 0.013) compared with total gene diversity (HT = 0.086). Analysis of molecular variance (AMOVA) showed that 83% of the molecular variance was explained by two major clusters, supporting a previous phenotypic study that suggested the existence of two distinct races of Illinois bundleflower. The less commonly occurring of the two races contained greater genetic diversity and originated from a localized region in the south central USA. The collection of more accessions from throughout the USA and focused on the Southeast would likely increase the genetic variation available to plant breeders.

Abbreviations: AFLP, amplified fragment length polymorphism • AMOVA, analysis of molecular variance • bp, base pair • HS, within-accession diversity • HT, total gene diversity • RAPD, random amplified polymorphic DNA • UPGMA, unweighted pair-group method, arithmetic average


    INTRODUCTION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
ILLINOIS BUNDLEFLOWER is an herbaceous warm-season N2–fixing perennial legume native to North America and useful in warm-season grass pastures. Illinois bundleflower is an excellent candidate for use as a perennial grain crop because accessions collected from the wild can produce in excess of 3000 kg ha-1 of seed (Adjei and Pitman, 1993) and the crude protein content of the seed is about 380 g kg-1 (Kulakow et al., 1990).

The potential of Illinois bundleflower as a forage crop has been illustrated in several experiments. ‘Sabine’ Illinois bundleflower (Muncrief and Heizer, 1985) was readily established in a Texas kleingrass (Panicum coloratum L.) pasture (Dovel et al., 1990). In the second, third, and fourth years of the study, kleingrass plots interseeded with Illinois bundleflower yielded 43, 47, and 45% more forage dry matter than noninterseeded plots in the respective years. Comparison tests have been made between warm-season grasses interseeded with Illinois bundleflower and warm-season grass monocultures (Posler et al., 1993). In the study, mixtures including Illinois bundleflower generally had forage yield and crude protein content more than double that of grass alone.

Kulakow (1999) conducted an extensive evaluation of phenotypic variation in 141 Illinois bundleflower accessions collected primarily from the Great Plains. Substantial variation was observed among accessions for traits including plant height, seed yield, seed size, survival, shatter resistance, and maturity. On the basis of phenotype, this study showed that there are at least two distinct races of Illinois bundleflower. Although variation within accession was not documented, accessions appeared uniform. Kulakow (1999) stated that molecular markers would be helpful for estimating within- versus between-population genetic diversity.

The AFLP technique for DNA fingerprinting enables visualization of restriction fragments without knowledge of nucleotide sequence (Vos et al., 1995). Fluorescence labeling and semi-automated detection of fragments have increased the speed and accuracy of the AFLP technique (Huang and Sun, 1999). Because fluorescence AFLP allows researchers to obtain genetic fingerprints rapidly without prior sequence information, the technique has been widely used in plant species with which little or no molecular research has been done. Fluorescence AFLP has been used to analyze genetic diversity in hops (Humulus lupulus L.), bermudagrass (Cynodon spp.), Onopordum thistles, ryegrasses (Lolium spp.), azaleas (Rhododendron simsii Planch.), and mulberry (Morus spp.) (Hartl and Seefelder, 1998; Zhang et al., 1999; O'Hanlon et al., 1999; De Riek et al., 1999; Roldán-Ruiz et al., 2000; and Sharma et al., 2000).

Although phenotypic diversity of Illinois bundleflower has been studied, the genetic diversity within and between Illinois bundleflower accessions is unknown. Knowledge of genetic diversity within and among accessions would enable plant breeders to choose parental sources that will generate diverse populations for selection. Knowledge of the distribution of genetic diversity within the species would also allow for an optimal germplasm preservation strategy, with future sampling focused on regions of greatest diversity. Our objective was to determine the genetic variation within and between Illinois bundleflower accessions with AFLP molecular markers.


    MATERIALS AND METHODS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Plant Materials and DNA Extraction
Seed of 50 Illinois bundleflower accessions was obtained from throughout the central USA (Fig. 1). Fifteen accessions were from the University of Minnesota Native Perennial Legume Collection, 25 accessions were from The Land Institute, Salina, KS, and 10 accessions were from the National Plant Germplasm System (USDA, ARS, 2001). For purposes of this study, the accessions were assigned numbers from 1 to 50 (Table 1).



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 1. Geographic origin of the Illinois bundleflower accessions used for AFLP analysis. Circled accession numbers belong to Cluster 2 and noncircled accession numbers belong to Cluster 1, as determined by the unweighted pair-group method, arithmetic average (UPGMA) clustering procedure applied to 159 polymorphic AFLP markers.

 

View this table:
[in this window]
[in a new window]
 
Table 1. Illinois bundleflower accession information and cluster grouping.

 
Plants of Accessions 1 through 20 were grown in a regularly irrigated field environment for a separate study. Tissue from these accessions was collected on July 11 from randomly selected plants in the third growing season after establishment. Plants of Accessions 21 through 50 were grown in the greenhouse and tissue was harvested from these plants about 3 wk after planting.

About 0.1 g of tissue was collected from developing leaves of three plants from every accession. The tissue was sampled and immediately frozen in liquid N2. The samples were ground in liquid N2 with a glass pestle. Total genomic DNA was extracted by means of the Dneasy Plant Mini Kit1 (Qiagen Inc., Valencia, CA) because the kit provided a clean extraction despite the fact that many samples were taken from mature field-grown plants. The DNA was concentrated to a final volume of about 10 µL with a DNA SpeedVac DNA 100 (Savant Instruments, Holbrook, NY) and quantified with an Ultrospec 3000 UV/Visible Spectrophotometer (Pharmacia Biotech Ltd., Cambridge, England).

The AFLP fingerprinting was performed with the AFLP Plant Mapping Kit (PE Applied Biosystems, Foster City, CA). The enzymes EcoRI and MseI (New England Biolabs, Inc., Beverly, MA) were used to digest 0.5 µg genomic DNA from each sample. Preselective amplification was performed with EcoRI-A and MseI-C primers. Selective amplification was performed with two fluorescently labeled primer sets: EcoRI-AA-Joe/MseI-CTA and EcoRI-AA-Joe/MseI-CAT. All amplifications were performed with a GeneAmp PCR System 9700 (PE Applied Biosystems). The selective amplification step was performed on two subsamples of each amplification product to replicate a portion of the laboratory procedure.

The AFLP products were separated with an automated DNA sequencer (ABI Prism 377, PE Applied Biosystems). The GeneScan-500 ROX (PE Applied Biosystems) internal size standard was used to accurately size the amplified fragments. The electropherograms generated by the sequencer were interpreted with GeneScan software (PE Applied Biosystems). Genotyper software (PE Applied Biosystems) was then used to create a list of fragments detected in each lane by fragment size. Fragments sized from 50 to 500 base pairs (bp) with a peak height >50 in the electropherogram were retained for subsequent analysis. Peakmatcher software (DeHaan et al., 2002) was used to convert the list of fragments detected in each sample into a binary (1/0) table for the presence or absence of each fragment in each sample. Peakmatcher was set to retain markers with repeatability >90%.

A separate matrix of Nei's unbiased genetic distances (Nei, 1978) was calculated by POPGENE (Yeh and Boyle, 1997) using markers derived from each of the two AFLP primer combinations. The correlation between the two distance matrices was determined with the MXCOMP function of NTSYSpc (Rohlf, 2000), and the significance of the correlation was tested with a Mantel test of 1000 random permutations.

POPGENE software was used to calculate Nei's unbiased genetic distance between accessions with all markers, including monomorphic markers. Nei's unbiased genetic distance is an accurate estimate of the number of gene differences per locus when populations are small (Nei, 1978). The matrix of unbiased genetic distances between accessions was analyzed by NTSYSpc. Dendrograms were generated by the unweighted pair-group method, arithmetic average (UPGMA), and the ball-cluster method of Jardine et al. (1969). Ordination analysis via nonmetric multidimensional scaling was used to represent the relationships between the accessions in two dimensions by means of the MDSCALE module of NTSYSpc. To evaluate goodness of fit of the UPGMA dendrogram and ordination analysis to the original distance matrix, the cophenetic correlation was calculated by means of the MXCOMP function of NTSYSpc. Within-accession diversity (HS) and total gene diversity (HT) (Nei, 1973) were calculated within the species and within two major groups by POPGENE software. The calculations were performed with all markers, both monomorphic and polymorphic.

The partitioning of molecular variance within and among groups and accessions was calculated by the AMOVA technique (Excoffier et al., 1992) in ARLEQUIN software (Schneider et al., 2001). Three models were used: all accessions within two clusters; accessions within Cluster 1; and accessions within Cluster 2. All significance tests were calculated by performing 1023 permutations.


    RESULTS AND DISCUSSION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Two AFLP primer combinations amplified a total of 222 unique fragments (Table 2). Of these, 159 fragments were polymorphic (72%). The average repeatability of the AFLP fragments across two replications was 98%. The primer pair AA-Joe, CAT amplified 130 fragments and the primer pair AA-Joe, CTA amplified 92 fragments. The unbiased genetic distance matrices calculated using each primer pair were similar (r = 0.965, P = 0.001) indicating that the information generated by the two primer pairs was consistent.


View this table:
[in this window]
[in a new window]
 
Table 2. Illinois bundleflower AFLP markers obtained from two primer combinations.

 
The present study is the first to use Peakmatcher software (DeHaan et al., 2002) to automate the gel scoring process. We verified the results produced with Peakmatcher by crosschecking a subset of the data with visual scoring using Genographer (Benham et al., 1999). The results of the two methods were consistent, and Peakmatcher enabled us to score a total of 600 lanes containing on average 49 fragments in about 6 h.

Cluster and Ordination Analysis
The dendrogram calculated from Nei's unbiased genetic distances between populations by the UPGMA method revealed two distinct clusters (Fig. 2). The unbiased genetic distance between the two clusters is 0.23. Cluster 1 comprises 42 accessions originating from throughout the central USA (Fig. 1). Cluster 2 comprises eight accessions originating from a localized region east of 96° W and south of 35°06' N (Table 1).



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 2. Illinois bundleflower accessions clustered by the unweighted pair-group method, arithmetic average (UPGMA) procedure applied to 159 polymorphic AFLP markers.

 
The UPGMA-generated dendrogram (Fig. 2) is in close agreement with the matrix of unbiased genetic distances (cophenetic correlation = 0.98). The high cophenetic correlation is due primarily to the existence of two very distinct clusters. The UPGMA method will produce a hierarchical clustering even if there is little evidence for clustering in the original data. Therefore, we generated a dendrogram of ball clusters (Fig. 3). A ball cluster dendrogram has the advantage of only producing robust clusters, but it has the disadvantage of losing a readily interpreted scale. Each accession within a ball cluster is more similar to every other accession within the cluster than it is to any accession outside the cluster (Jardine et al., 1969).



View larger version (17K):
[in this window]
[in a new window]
 
Fig. 3. Ball cluster dendrogram of Illinois bundleflower accessions based on 159 polymorphic AFLP markers.

 
The ball cluster dendrogram (Fig. 3) reveals little clustering among the majority of the accessions. The clustering that is present is generally in agreement with geographic origin. For instance, Accessions 25, 37, and 38, which originated from Kansas, all cluster together. Exceptions are Accessions 14 and 16 from Illinois that clustered with Accessions 40, 42, and 47 from Oklahoma, and Accession 15 from Iowa that clustered with Accession 29 from Texas. All accessions obtained from latitude greater than 36° N were placed in the same third-level ball cluster, except Accessions 3 from Minnesota and 10 from South Dakota. Therefore, these two accessions may be important sources of genetic variation within northern germplasm.

Ordination analysis by nonmetric multidimensional scaling (Fig. 4) provided an excellent summary of the data in two dimensions (cophenetic correlation = 1.00). Ordination analysis represents large distances most accurately, while cluster analysis (Fig. 2 and Fig. 3) represents the fine structure of the data most accurately (Sneath and Sokal, 1973). In the ordination analysis, Accessions 44 and 26 are placed on the horizontal axis between the centers of the two clusters. Therefore, although Accessions 44 and 26 are more similar to Cluster 1 than Cluster 2, they somewhat resemble Cluster 2. Accession 44 resembles Cluster 2, probably because it was obtained from the geographical margin between the two clusters (Fig. 1). Accession 26 is isolated both genetically and geographically from the other accessions (Fig. 1 and Fig. 4). In a previous phenotypic study by Kulakow (1999), Accession 26 clustered with Accessions 24, 35, 36, and 39, which are found in Cluster 2 in this study. Therefore, Accession 26 is unique in its phenotypic resemblance to accessions in Cluster 2 and its genetic similarity to accessions in Cluster 1.



View larger version (11K):
[in this window]
[in a new window]
 
Fig. 4. Illinois bundleflower accessions plotted from nonmetric multidimensional scaling analysis of 159 polymorphic AFLP markers. The dashed line separates the two major clusters from unweighted pair-group method, arithmetic average (UPGMA) analysis (Fig. 2). The tightly grouped unlabeled points represent Accessions 1, 2, 6, 8, 9, 11, 12, 14, 15, 16, 17, 18, 20, 21, 22, 25, 29, 30, 37, 38, 40, 42, 43, and 47.

 
Gene Diversity and AMOVA
The total gene diversity across all accessions (HT) was 0.086 and the within-accession gene diversity (HS) was 0.013 (Table 3). Low HS relative to HT is consistent with the species being approximately 80% self-pollinating (DeHaan, unpublished data, 2000). The values of HT and HS for Cluster 2 exceed those of Cluster 1 by 61.3 and 62.5%, respectively, which indicates that Cluster 2 contains greater genetic diversity both among and within accessions. Although gene diversity estimates are commonly used to make comparisons between species, calculations derived from AFLP data should not be compared across studies because AFLP overestimates the number of loci and underestimates the number of alleles (Caicedo et al., 1999). Because of these factors and because fluorescence AFLP can detect many fragments that occur at low frequency, the low diversity estimates obtained may be due to the marker technology used rather than to low genetic diversity in the species. Thus relative, but not absolute, diversity estimates from this study should be considered reliable.


View this table:
[in this window]
[in a new window]
 
Table 3. Total gene diversity among and within (HT) and gene diversity within (HS) 50 accessions of Illinois bundleflower.

 
The AMOVA analysis showed that 83% of the molecular variance was explained by the two major clusters (Table 4). This result indicates that the majority of the genetic variation in the Illinois bundleflower accessions is found between two distinct types represented as Cluster 1 and Cluster 2. Kulakow (1999) described four clusters based on phenotypic data, but principal component analysis showed that Clusters 1, 2, and 3 formed one large group. Cluster 4 was distinctly different from the other three clusters, which led to the conclusion that there are two distinct races of Illinois bundleflower. The AFLP analysis lends support to this conclusion. Four of the five accessions that Kulakow's (1999) analysis placed Cluster 4 were included in Cluster 2 in the current study (Table 1). Accession 26 was the exception, although ordination analysis revealed that this accession was distinctly different from most other Cluster 1 accessions (Fig. 4).


View this table:
[in this window]
[in a new window]
 
Table 4. Analysis of molecular variance (AMOVA) of 159 polymorphic AFLP markers from 50 Illinois bundleflower accessions.

 
The AMOVA-derived estimate of population differentiation ({Phi}ST) is 0.905, which is high when compared with studies of self-pollinated species. Nybom and Bartish (2000) reported that, across eight studies using random amplified polymorphic DNAs (RAPDs), average {Phi}ST reported for selfing species was 0.70. The exceptionally large {Phi}ST estimated for Illinois bundleflower can be attributed mainly to the existence of two very distinct clusters. Allocation of molecular variance within and between accessions is similar for both clusters. About 45% of the molecular variance is among accessions and about 55% is within accessions (Table 4). This means that when the two major clusters are disregarded, about 45% of the genetic variation is between accessions.

Geographic Distribution of Genetic Variation
The mean geographic distance between accessions is 763 km in Cluster 1 and 411 km in Cluster 2. Although accessions in Cluster 2 originate from a more localized geography, they are more diverse. Because accessions in Cluster 2 originated from the southeastern range of the area represented in this study, additional genetic diversity might be obtained by collecting seed of plants from the southeastern USA. According to the PLANTS Database (USDA, NRCS, 2001), the range of Illinois bundleflower extends south and east to Alabama, Georgia, Florida, South Carolina, North Carolina, Virginia, Maryland, and Pennsylvania. Accessions from these states are not available through the National Plant Germplasm System (USDA, ARS, 2001), but obtaining accessions from this region could increase the genetic diversity available to plant breeders working with the species.

Implications for Plant Breeding Programs
Accessions in Cluster 1 are the more common type among the accessions used in this experiment (Table 1) and in a previous study of 141 accessions (Kulakow, 1999). Kulakow (1999) described many of the accessions in this cluster as being erect, high seed yielding, and having good survivability, while accessions belonging to Cluster 2 were described as important for their large seed size and vigorous seedling growth. Because polymorphisms between the two races are abundant, the phenotypic differences between the clusters are likely to be a result of diversity across the genome rather than because of a major gene effect. Therefore, accessions of Cluster 2 are likely to be important sources of variation for many quantitative traits, including seed and forage yield.

The germplasm collections presently maintained by the University of Minnesota, The Land Institute, and the National Plant Germplasm System (USDA, ARS, 2001) are weighted heavily toward the more genetically uniform Cluster 1 type and contain few accessions of the more genetically diverse Cluster 2 type of Illinois bundleflower. To maximize the genetic diversity available for use by plant breeders, additional collections need to be made at the periphery of Illinois bundleflower's range in the southern, eastern, and western United States. The use of genetically uniform seed sources such as Sabine Illinois bundleflower for widespread roadside revegetation could threaten the diversity within the species by contaminating or overwhelming it. Local ecotypes should be used for this purpose whenever possible.

We obtained valuable information about several specific accessions. Although Accession 44 was placed in Cluster 1, it was nearly intermediate between the two clusters of Illinois bundleflower. It may be useful to examine the phenotype of this accession because of its potential to contain a combination of useful traits from the two clusters. Kulakow (1999) found Accession 26 to be most similar phenotypically to accessions that were placed in Cluster 2 in the present study. He also described it as being an important source of shatter resistance and large seed size. Because it is genetically similar to the high seed yielding accessions of Cluster 1, shatter resistance might be readily introduced from this accession into Cluster 1 types without substantially reducing seed yield or survival. Accessions 3 and 10 are genetically distant from other accessions collected in the northern United States. Therefore, they may be important sources of diversity within winter-hardy germplasm. The existence of these unique accessions underscores the importance of preserving populations of Illinois bundleflower native to the northern region. Illinois bundleflower has become a rare species in Minnesota because of draining wetlands and widespread herbicide usage (Coffin and Pfannmuller, 1988).

We have documented low within-accession genetic diversity relative to total genetic diversity present in Illinois bundleflower. This pattern underscores the importance of evaluating large numbers of accessions rather than many individuals per accession when searching for valuable traits for plant breeding programs.


    ACKNOWLEDGMENTS
 
We gratefully acknowledge funding provided by The Land Institute, Salina, KS, and the University of Minnesota Graduate School. We are also grateful to Robb DeHaan and The Land Institute for obtaining many of the accessions used in this study.


    NOTES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Contribution of the Minnesota Agric. Exp. Stn.

1 Names are necessary to report factually on available data; however, the Univ. of Minnesota neither guarantees nor warrants the standard of the product, and the use of the name by the Univ. of Minnesota implies no approval of the product to the exclusion of others that may be suitable. Back

Received for publication October 8, 2001.


    REFERENCES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 




This article has been cited by other articles:


Home page
ANN BOT (LOND)Home page
E. L. PLOSCHUK, G. A. SLAFER, and D. A. RAVETTA
Reproductive Allocation of Biomass and Nitrogen in Annual and Perennial Lesquerella Crops
Ann. Bot., July 1, 2005; 96(1): 127 - 135.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
L. R. DeHaan, N. J. Ehlke, C. C. Sheaffer, R. L. DeHaan, and D. L. Wyse
Evaluation of Diversity among and within Accessions of Illinois Bundleflower
Crop Sci., July 1, 2003; 43(4): 1528 - 1537.
[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 Similar articles in this journal
Right arrow Similar articles in ISI 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 ISI Web of Science (4)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by DeHaan, L. R.
Right arrow Articles by Wyse, D. L.
Right arrow Search for Related Content
PubMed
Right arrow Articles by DeHaan, L. R.
Right arrow Articles by Wyse, D. L.
Agricola
Right arrow Articles by DeHaan, L. R.
Right arrow Articles by Wyse, D. L.
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