Estimation of micro-biota in the Upo wetland using eukaryotic barcode molecular markers
Estimation of micro-biota in the Upo wetland using eukaryotic barcode molecular markers
Journal of Ecology and Environment. 2011. Sep, 34(3): 323-331
Copyright ©2011, The Ecological Society of Korea
This is an Open Access article distributed under the terms of the CreativeCommons Attribution Non-Commercial License( permits unrestricted non-commercial use,distribution, and reproduction in any medium, provided the original workis properly cited.
  • Received : April 28, 2011
  • Accepted : July 10, 2011
  • Published : September 01, 2011
Export by style
Cited by
About the Authors
Hyun Chul, Park
School of Biological Sciences, Seoul National University, Seoul 151-747, Korea
Chang Hwan, Bae
National Institute of Biological Resources, Environmental Research Complex, Incheon 404-708, Korea
Jumin, Jun
Department of Life Sciences, Ewha Womans University, Seoul 120-750, Korea
Myounghai, Kwak
National Institute of Biological Resources, Environmental Research Complex, Incheon 404-708, Korea

Biodiversity and the community composition of micro-eukaryotic organisms were investigated in the Upo wetland in Korea using molecular analysis. Molecular identification was performed using cytochrome oxidase I (COI) and small subunit ribosomal DNA (SSU rDNA). The genomic DNA was isolated directly from soil samples. The COI and SSU rDNA regions were amplified using universal primers and then sequenced after cloning. In a similarity search of the obtained sequences with BLAST in the Genbank database, the closely related sequences from NCBI were used to identify the amplified sequences. A total of six eukaryotic groups (Annelida, Arthropoda, Rotifera, Chlorophyta, Bacillariophyta, and Stramenopiles) with COI and six groups (Annelida, Arthropoda, Rotifera, Alveolata, Fungi, and Apicomplexa) with SSU rDNA genes were determined in the Upo wetland. Among 38 taxa in 20 genera, which are closely related to the amplified sequences, 10 genera (50%) were newly reported in Korea and five genera (25%) were shown to be distributed in the Upo wetland. This approach is applicable to the development of an efficient method for monitoring biodiversity without traditional taxonomic processes and is expected to produce more accurate results in depositing molecular barcode data in the near future.
Micro eukaryotes such as small invertebrates, larvae and eggs, and plankton are very abundant in the environment and play important roles in supporting the ecosystem(Diez et al. 2001). However, they are often neglected in most biodiversity surveys due to difficulty of identification,a lack of specialists to identify such organisms,time-consuming soil screening processes, and a lack of taxonomic treatment of small organisms (Lefranc et al.2005). While most biota investigation has been generally performed using a taxonomic method based on morphological analysis, species identification systems using molecular analysis have been established for the identification of microbes, as a tool for microbial community biology, and for the classification of microbial population data from environmental samples (Olsen et al. 1986, Jahn et al. 2007, Zwolinski 2007).
Recently, DNA-based identification (DNA barcoding)has been started to the animal and plant kingdoms (Hebert et al. 2004, Barrett and Hebert 2005, Kress and Erickson 2007). Although vast sequence data for DNA barcoding regions in animals and plants have been deposited, DNA barcoding is mainly limited to the identification of a single sample. While molecular species identification has been popularly adopted to microbial community analysis for environmental samples (Sogin et al. 2006, Huber et al. 2007), it has been seldom applied to the identification of eukaryotes from environmental samples (Creer et al. 2010). Several researchers have studied the biodiversity of micro eukaryotes in extreme environments such as Antarctic soil and tropical lakes using small subunit ribosomal DNA (SSU rDNA) (Lawley et al. 2004, Lefranc et al. 2005, Fell et al. 2006).
In the present study, we describe a molecular survey of organisms from soil samples taken from the Upo wetland in Korea using two genes, cytochome oxidase I (COI) and SSU rDNA. As the COI gene is a mitochondrial gene, it has a universal primer for amplification and a better phylogenetic signal than other mitochondrial genes in most animal phyla (Folmer et al. 1994). For these reasons, COI is utilized as a DNA barcode for the discrimination of species (Hebert et al. 2003). However, since the molecular community estimation of environmental samples has been performed using 18S SSU rDNA (Savin et al. 2004, Behnke et al. 2006, Medinger et al. 2010), we also analyzed SSU rDNA for comparison with the COI results. Moreover, SSU rDNA is widely applied in phylogenetic analyses, constituting a vast portion of all available sequences and covering diverse taxa. Thus, we expected that SSU rDNA could compensate for the limitation of COI, which has been less studied compared to other taxa except metazoan animals. Through this study, we identified numerically dominant groups in the study area as well as the possible existence of previously unknown organisms.
- Study site
The Upo wetland is the largest natural lowland swamp in Korea and was formed from ground subsidence around 140 million years ago ( Fig. 1 ). The Upo wetland was designated as an ecological conservation area by the Ministry of Environment in 1997 and registered as a Ramsar wetland. It has an average water depth of 1 m (Baek 1988) and has maintained an average temperature of 13°C for 30 years. The annual rainfall measures 1,233 mm and mainly occurs during the summer period (Kim et al. 2004). Water comes from the Topyeong stream but is directly influenced by a flood from the adjacent Nakdong River (Do et al. 2007). The biota of the Upo region has been investigated by several research groups who have focused on macro invertebrates, which are diverse aquatic plant and vertebrates (amphibian and fish). However, micro eukaryotes have been rarely studied, except for phytoplankton (Kim 2001, Bae et al. 2004, Nakdong River Basin Environmental Office 2006).
- Sample collection
Soil samples were collected on 2 March in 2009 (35°32?49.51˝ N, 128°24?40.01˝ E) ( Fig. 1 ). The collecting site was submerged under water, and it contained many suspended plant particles. Using a shovel at a depth of 10 cm near the water’s edge, around 2 kg of soil from the top layer, including the surface layer, suspended particles, and water, was collected into a plastic bag. The sample was transferred to the laboratory and kept fresh in an ice box. Three 50 mL aliquots of the soil samples were then placed into sterile Falcon tubes using sterile spatulas and stored at -70°C in a deep freezer before genomic DNA extraction.
- DNA extraction and sequencing
Genomic DNA from each of the three aliquots was isolated using the PowerMax™ Soil DNA Isolation Kit (Mo Bio, West Carlsbad, CA, USA) according to the manufacturer’s instructions. PCR amplification was performed using a premix taq of Top-Taq™ PreMix (2×) (Takara, Tokyo, Japan) with 10 ng of genomic DNA and 0.2 μM primers in a final volume of 50 μL. The COI region was amplified with LCO1490 (5′-GGTCAACAAATCATAAAGATATTGG-3) and HCO2198 (5′-TAAACTTCAGGGTGACCAAAAAATCA-3) (Former et al. 1994), whereas the SSU rDNA region was amplified with EukA (5′-AACCTGGTTGATCCTGCCAGT-3) and EukB (5-TGATCCTTCTGCAGGTTCACCTAC-3) (Jeong et al. 2005). PCR were composed of the following steps: pre?denaturation (10 min at 94°C), 35 cycles of denaturation (1 min at 94°C), annealing for 2 min at 51°C for both COI and SSU rDNA, and extension (2 min at 72°C) and final extension for 7 min at 72°C using a Mastercycler (Eppendorf, Hamburg, Germany). PCR products were visualized on 1% agarose gel and purified with a Core-One™ PCR purification kit (CoreBio, Seoul, Korea). TA cloning was performed using the pGEM-T easy vector system (Promega, Madison, WI, USA) by Macrogen Inc. (Seoul, Korea). The forty eight clones for each marker were sequenced with forward
Lager Image
Study site in the Upo wetland. (a) Map of Korea with sampling site indicated. (b) Map of Upo swamp area. Circled dot: site where soil samples were collected. (c) Overview of the study site.
(SP6) and reverse (T7) primers using an ABI 3730XL SEQUENCER (Applied Biosystems, Foster City, CA, USA).
- Data analyses
After sequence editing with SEQUENCHER ver. 4.8, sequence alignment was initially performed using Clustal X 1.81 (Thompson et al. 1997). To obtain species information for each operational taxonomic unit (OTU), BLAST search was performed using the blastn program against the NCBI database to retrieve the taxonomic classification of the best hit. In the case of low quality hits with query coverage < 90%, identity < 75% for COI and query coverage < 90%, identity < 87% for SSU rDNA, the sequences were discarded. As a neighbor-joining tree was constructed using MEGA ver. 4.0 from the aligned sequences and the sequences of the best hit (Saitou and Nei, 1987, Tamura et al. 2007), phylogeny groups for each sequence were assigned, referred from the taxonomic information of the best hits. The obtained OTUs were compared with the biota reports on the Upo wetland (Bae et al. 2004, Nakdong River Basin Environmental Office 2006), the list of animals in Korea (The Korean Society of Systematic Zoology 1997), and the check list of insects from Korea (The Entomological Society of Korea and Korean Society of Applied Entomology 1994).
Forty-seven and 44 sequences for COI and SSU rDNA, respectively, were identified. All obtained sequences were deposited in NCBI (GU070858-GU070938). The sizes of COI and SSU rDNA were about 700 and 1,700 base pairs (bp), respectively. After four COI and six SSU rDNA sequences were removed from the data matrix according to the cut-off, a total of 38 closest relatives in 20 genera (7 groups) were identified using blast search ( Table 1 ). Interestingly, only five genera among the 22 total genera have been previously found in the Upo wetland.
Forty-three COI sequences had high similarity with nine species in three animal phyla (Annelida, Arthropoda, and Rotifera), 11 species in a stramenopiles group, and two species in one green algae group (Viridiplantae: Chlorophyta). Especially, stramenopiles groups (Oomycetes and Bacillariophyta) were widely represented compared to the animal groups in the phylogenetic tree ( Fig. 2 a). Bacillariophyta (Diatom) is known to be common in freshwater systems. In the Upo wetland area, 36 diatom taxa have been previously reported (Nakdong River Basin Environmental Office 2006). The other Stamenopiles group is Oomycetes, known as water mold, and is usually pathogenic to plants. Since the Upo wetland is surrounded by a rice paddy and crop field and contains various wetland plants, it is not surprising that Oomycetes OTUs were dominantly amplified. However, the biodiversity of this group has not been investigated in the Upo wetland, and blast results did not show high similarity between obtained OTUs and their best hit sequences (below 90%) ( Table 1 ), suggesting that further analysis is required to examine the distribution of oomycetes in Upo.
COI sequences from the genus Paratanytarsus (Arthropod: Chironomidae) and genus Dero (Annelid: Naididae) were matched from Genbank sequences with high identity scores of 100% and 98%, respectively ( Table 1 ). However, these genera have not yet been reported in the Korean peninsula (The Entomological Society of Korea and Korean Society of Applied Entomology 1994, The Korean Society of Systematic Zoology 1997). In general, COI variations within species of Arthropoda and Annelida were reported to occur at a rate of less than 3% (Quan et al. 2001, Erseus and Kvist 2007). This highly suggests the existence of new or unrecorded Paratanytarsus and Dero in the Upo study area.
Thirty-eight SSU rDNA sequences were identified; eight from three animal phyla (Annelida, Arthropoda, and Rotifera), five species in Ciliophora, one Apicomplexa, one Ascomycota, and three unidentified eukaryote sequences. Using SSU rDNA primers, sequences closely related to Copepoda and Arthropoda were frequently obtained (37%) ( Fig. 2 b). Further, various genera of Ciliophora were obtained using SSU rDNA primers. In a molecular survey of the fresh-water system, Ciliophora were shown to be the dominant taxa using SSU rDNA primers (Behnke et al. 2006, Medlinger et al. 2010).
Biodiversity assessment is the first step to understanding the characteristics of ecosystems in a study area. However, a lack of experts and the difficulty of specimen isolation hinder the detection of small eukaryotes by biodiversity survey. Thus, we examined alternative molecular methods for screening biodiversity at a study site. In this study, we showed that approximately three fourths of all genera in the Upo area were previously unreported ( Table 1 ), suggesting a large number of species that have been ignored in biodiversity screenings. Especially, COI sequences, which match genus Paratanytarsus and genus Dero with high similarity, indicates their existence in Upo area, even though it is never been recognized in Upo area as well as in Korean peninsula.
Hebert et al. (2003) reported the identification of COI sequences with a good resolution in seven major animal phyla: Platyhelminthes, Nematoda, Mollusca, Annelida, Echinodermata, Chordata, and Arthropod. Based on this, a huge number of COI sequences from animals were deposited to develop Barcode of Life Data Systems (BOLD, ) as well as the NCBI genbank database. Our results also show that the COI region was applicable to the identification of organisms in the phyla Arthropoda and Annelida from environmental samples. However, COI alone could not contribute to the identification of all organisms due to a lack of sequence information in certain taxa and the absence of target genes in some eukaryotes (Scicluna et al. 2006). Protist groups like stramenopiles have been identified using ribosomal DNAs such as SSU rDNA and LSU rDNA rather than COI (Diez et al. 2001, Fell et al. 2006, Scicluna et al. 2006).
A total of 15 species including fungi groups were searched with 38 SSU rDNA sequences ( Table 1 ). The same animal phyla as a COI tree have been identified, and alveolate groups, a major line of protest, were found in SSU rDNA sequences ( Fig. 2 ). In the amplification using SSU rDNA primers, the Stamenopile and Viridiplantae groups were not identified. Although the genera Stentor and Eucyclops are common in freshwater systems in Korea, they have not been reported in Upo. This research also revealed that a molecular approach based on barcode regions is able to facilitate the discovery of unreported taxa.
The animal phyla Annelida, Arthropoda, and Rotifera were identified at similar proportions with both markers, whereas Chlorophyta, Stamenopiles, Bacillariophyta were only amplified with COI marker and Ciliophora,
List of taxa identified in this study compared with those reported in the Upo wetlandCOI, cytochrome oxidase I; SSU, small subunit. *A species of the corresponding genus has been reported in the Upo wetland.
Lager Image
List of taxa identified in this study compared with those reported in the Upo wetland COI, cytochrome oxidase I; SSU, small subunit. *A species of the corresponding genus has been reported in the Upo wetland.
Lager Image
Phylogenetic tree constructed using (a) cytochrome oxidase I and (b) small subunit ribosomal rDNA gene sequences. Bootstrap support (1000 replications) ≻ 50% is indicated.
Lager Image
Lager Image
Proportion of groups identified using cytochrome oxidase I (COI) marker and small subunit ribosomal rDNA (SSU rDNA) markers.
Apicomplexa, Ascomycota with SSU rDNA marker ( Fig. 3 ). However, we could not determine whether or not the proportions of obtained sequences reflect the number of species/taxa at the study site. Firstly, the total volume of genomic DNA varies according to the species since the number of cells can differ from single cell organisms such as protest to relatively larger multicellular organisms. Moreover, copies of especially ribosomal DNA vary among taxonomic groups. Thus, larger species with more cells and with higher copy numbers of genes would be more frequently detected compared to others (Potvin and Lovejoy 2009). Lastly, the number of final amplified sequence copies is dependent on PCR reaction efficiency, which is mainly affected by primer affinity. Unfortunately, there is no perfect primer to amplify specific target regions from eukaryotic taxa. This means that diversity may be easily underestimated. However, this molecular estimation method would be more efficient than traditional groups in terms of covering diverse taxonomic groups, since it is impossible to gather experts to cover all taxonomic groups in many cases.
This study therefore demonstrates the possibility of molecular barcoding systems to survey biodiversity in a specific environment, especially for identifying micro flora and fauna. To effectively apply this method for biodiversity surveying as well as biodiversity monitoring, the following questions should be considered. Firstly, amplified multiple sequences were separated by cloning steps before sequencing. However, due to the high expense of the cloning step, it would be somewhat difficult to represent all of the species of a region with a limited number of sequences. When performing a new DNA sequencing method such as pyrosequencing to compensate for these drawbacks, a quantitative approach based on the number of DNA sequences should be implemented (Valentini et al. 2008). Secondly, we totally depended on the Genbank database using the blastn program to identify species (GenBank, ). Although more sequence data have been deposited in a public database, the quality of sequence data is not entirely reliable due to imperfections such as misidentification and taxonomic problems (Harris 2003, Valentini et al. 2008). Thus, it is important to complete the reliable reference database for DNA barcoding such as BOLD and the Ribosomal Database Project (RDP, ). Moreover, if more gene regions such as internal transcribed spacer and large subunit ribosomal DNA are analyzed in parallel, more abundant and diverse species will be found.
This work was supported by a research grant from the National Institute of Biological Resources, Republic of Korea. The authors thank to Young-Mi Lee, Jung-Hyun An, and Sun-Nam Park for technical support and Hyun-Kyung Oh for helpful comments.
Bae YJ , Jo SI , Hoang DH , Lee HG , Na KB 2004 Biodiversity and community composition of benthic macroinvertebrates from Upo wetlands in Korea Korean J Environ Ecol KISTI1.1003/JNL.JAKO200411922098352 18 75 - 91
Baek SY 1988 Ecological studies on the Upo wetland Chung-Ang University Seoul Korea
Barrett RDH , Hebert PDN 2005 Identifying spiders through DNA barcodes Can J Zool 83 481 - 491    DOI : 10.1139/z05-024
Behnke A , Bunge J , Barger K , Breiner HW , Alla V , Stoeck T 2006 Microeukaryote community patterns along an O2/H2S gradient in a supersulfidic anoxic fjord (Framvaren Norway) Appl Environ Microbiol 72 3626 - 3636    DOI : 10.1128/AEM.72.5.3626-3636.2006
Creer S , Fonseca VG , Porazinska DL , Giblin-Davis RM , Sung W , Power DM , Packer M , Carvalho GR , Blaxter ML , Lambshead PJD , Thomas WK 2010 Ultrasequencing of the meiofaunal biosphere: practice pitfalls and promises Mol Ecol 19 (Suppl 1) 4 - 20    DOI : 10.1111/j.1365-294X.2009.04473.x
Diez B , Pedros-Alio C , Massana R 2001 Study of genetic diversity of eukaryotic picoplankton in different oceanic regions by small-subunit rRNA gene cloning and sequencing Appl Environ Microbiol 67 2932 - 2941    DOI : 10.1128/AEM.67.7.2932-2941.2001
Do Y , Jang MH , Kim DK , Joo GJ 2007 Change of carabid beetle (Coleoptera Carabidae) diversity and species composition after flooding events in Woopo wetlands Korean J Limnol KISTI1.1003/JNL.JAKO200709905951465 40 346 - 351
Erseus C , Kvist S 2007 COI variation in Scandinavian marine species of Tubificoides (Annelida: Clitellata: Tubificidae J Mar Biol Assoc U K 87 1121 - 1126
Fell JW , Scorzetti G , Connell L , Craig S 2006 Biodiversity of micro-eukaryotes in antarctic dry valley soils with <5% soil moisture Soil Biol Biochem 38 3107 - 3119    DOI : 10.1016/j.soilbio.2006.01.014
Folmer O , Black M , Hoeh W , Lutz R , Vrijenhoek R 1994 DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates Mol Mar Biol Biotechnol 3 294 - 299
Harris DJ 2003 Can you bank on genbank? Trends Ecol Evol 18 317 - 319    DOI : 10.1016/S0169-5347(03)00150-2
Hebert PDN , Cywinska A , Ball SL , deWaard JR 2003 Biological identifications through DNA barcodes Proc R Soc Lond B 270 313 - 321    DOI : 10.1098/rspb.2002.2218
Hebert PDN , Stoeckle MY , Zemlak TS , Francis CM 2004 Identification of birds through DNA barcodes PLoS Biol 2 e312 -    DOI : 10.1371/journal.pbio.0020312
Huber JA , Welch DBM , Morrison HG , Huse SM , Neal PR , Butterfield DA , Sogin ML 2007 Microbial population structures in the deep marine biosphere Science 318 97 - 100    DOI : 10.1126/science.1146689
Jahn R , Zetzsche H , Reinhardt R , Gemeinholzer B 2007 Diatoms and DNA barcoding: A pilot study on an environmental sample BGBM Press Berlin Proceedings of the 1st Central European Diatom Meeting. 2007 Mar 23-25 Berlin-Dahlem.
Jeong HJ , Kim JS , Park JY , Kim JH , Kim S , Lee I , Lee SH , Ha JH , Yih WH 2005 Stoeckeria algicida n. gen. n. sp. (Dinophyceae) from the coastal waters off southern Korea: morphology and small subunit ribosomal DNA gene sequence J Eukaryot Microbiol 52 382 - 390
Kim GY , Lee CW , Joo GJ 2004 The evaluation of early growth pattern of miscanthus sacchariflorus after cutting and burning in the Woopo wetland Korean J Limnol KISTI1.1003/JNL.JAKO200418317182871 37 255 - 262
Kim HS 2001 Seasonal changes of phytoplankton community in the Woopo and Mokpo swamp Korean J Limnol KISTI1.1003/JNL.JAKO200118317177235 34 90 - 97
Kress WJ , Erickson DL 2007 A two-locus global DNA barcode for land plants: the coding rbcL gene complements the non-codingtrnH-psbAspacer region PLoS One 2 e508 -    DOI : 10.1371/journal.pone.0000508
Lawley B , Ripley S , Bridge P , Convey P 2004 Molecular analysis of geographic patterns of eukaryotic diversity in antarctic soils Appl Environ Microbiol 70 5963 - 5972    DOI : 10.1128/AEM.70.10.5963-5972.2004
Lefranc M , Thenot A , Lepere C , Debroas D 2005 Genetic diversity of small eukaryotes in lakes differing by their trophic status Appl Environ Microbiol 71 5935 - 5942    DOI : 10.1128/AEM.71.10.5935-5942.2005
Medinger R , Nolte V , Pandey RV , Jost S , Ottenwalder B , Schlotterer C , Boenigk J 2010 Diversity in a hidden world: potential and limitation of next-generation sequencing for surveys of molecular diversity of eukaryotic microorganisms Mol Ecol 19 (Suppl 1) 32 - 40    DOI : 10.1111/j.1365-294X.2009.04478.x
Nakdong River Basin Environmental Office 2006 A Report for Biota of Lake (Upo) Nakdong River Basin Environmental Office Changwon in Korean
Olsen GJ , Lane DJ , Giovannoni SJ , Pace NR , Stahl DA 1986 Microbial ecology and evolution: a ribosomal RNA approach Annu Rev Microbiol 40 337 - 365
Potvin M , Lovejoy C 2009 PCR-based diversity estimates of artificial and environmental 18S rRNA gene libraries J Eukaryot Microbiol 56 174 - 181    DOI : 10.1111/j.1550-7408.2008.00386.x
Quan J , Lu XM , Zhuang Z , Dai J , Deng J , Zhang YP 2001 Low genetic variation of Penaeus chinensis as revealed by mitochondrial COI and 16S rRNA gene sequences Biochem Genet 39 279 - 284    DOI : 10.1023/A:1010234816514
Savin MC , Martin JL , LeGresley M , Giewat M , Ronney-Varga J 2004 Plankton diversity in the bay of fundy as measured by morphological and molecular methods Microb Ecol 48 51 - 65    DOI : 10.1007/s00248-003-1033-8
Saitou N , Nei M 1987 The neighbor-joining method: a new method for reconstructing phylogenetic trees Mol Biol Evol 4 406 - 425
Scicluna SM , Tawari B , Clark CG 2006 DNA barcoding of Blastocystis Protist 157 77 - 85    DOI : 10.1016/j.protis.2005.12.001
Sogin ML , Morrison HG , Huber JA , Welch DM , Huse SM , Neal PR , Arrieta JM , Herndl GJ 2006 Microbial diversity in the deep sea and the underexplored “rare biosphere” Proc Natl Acad Sci U S A 103 12115 - 12120    DOI : 10.1073/pnas.0605127103
Tamura K , Dudley J , Nei M , Kumar S 2007 MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0 Mol Biol Evol 24 1596 - 1599    DOI : 10.1093/molbev/msm092
The Entomological Society of Korea and Korean Society of Applied Entomology 1994 Check List of Insects from Korea Konkuk University Press Seoul
The Korean Society of Systematic Zoology 1997 List of Animals in Korea (Excluding Insects) Academy Publishing Co. Seoul
Thompson JD , Gibson TJ , Plewniak F , Jeanmougin F , Higgins DG 1997 The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools Nucleic Acids Res 25 4876 - 4882    DOI : 10.1093/nar/25.24.4876
Valentini A , Pompanon F , Taberlet P 2008 DNA barcoding for ecologists Trends Ecol Evol 24 110 - 117
Zwolinski MD 2007 DNA sequencing: strategies for soil microbiology Soil Sci Soc Am J 71 592 - 600    DOI : 10.2136/sssaj2006.0125