Microbiome Announcement | Open | Published:
The PAMPA datasets: a metagenomic survey of microbial communities in Argentinean pampean soils
Microbiomevolume 1, Article number: 21 (2013)
Soil is among the most diverse and complex environments in the world. Soil microorganisms play an essential role in biogeochemical cycles and affect plant growth and crop production. However, our knowledge of the relationship between species-assemblies and soil ecosystem processes is still very limited. The aim of this study was to generate a comprehensive metagenomic survey to evaluate the effect of high-input agricultural practices on soil microbial communities.
We collected soil samples from three different areas in the Argentinean Pampean region under three different types of land uses and two soil sources (bulk and rhizospheric). We extracted total DNA from all samples and also synthetized cDNA from rhizospheric samples. Using 454-FLX technology, we generated 112 16S ribosomal DNA and 14 16S ribosomal RNA amplicon libraries totaling 1.3 M reads and 36 shotgun metagenome libraries totaling 17.8 million reads (7.7 GB). Our preliminary results suggested that water availability could be the primary driver that defined microbial assemblages over land use and soil source. However, when water was not a limiting resource (annual precipitation >800 mm) land use was a primary driver.
This was the first metagenomic study of soil conducted in Argentina and our datasets are among the few large soil datasets publicly available. The detailed analysis of these data will provide a step forward in our understanding of how soil microbiomes respond to high-input agricultural systems, and they will serve as a useful comparison with other soil metagenomic studies worldwide.
The Argentine Pampas is a plain area of 60 million ha. Because of its large expanse and high yields, it is one of the most productive areas for grain crop production in the world . Indeed, 90% of the pampean surface is currently used for high-input agricultural purposes. Argentina is currently the third and fourth world producer of soybean and maize, respectively . This production is mostly concentrated in the pampean region.
Since 1980, agriculture has rapidly expanded in the region, replacing grasslands, with the widespread adoption of limited tillage systems, particularly no-till with crop rotation . These practices have been reported to preserve surface water, prevent soil erosion and return nutrients to soil [4–6]. However, concerns remain regarding the impact of these practices on soil quality, microbial diversity and community assemblages.
Changes in microbial communities throughout the Argentine Pampas are poorly reported. Most studies have focused on the tillage effects on microbial biomass or specific microbial activities such as the utilization of specific substrates, extracellular enzyme production, or mineralization [7–9]. Other studies have focused on well-studied and particular bacterial taxa rather than the microbial community structure itself [10, 11]. Studies conducted with an ecological approach have usually focused on the individual effects of land use such as the application of herbicides [12, 13]. In such cases, community variability was assessed using classical fingerprinting techniques (such as RFLP and DGGE), which only capture the most dominant species in the environment [14, 15]. In this regard, classical approaches are inadequate for describing highly diverse soil microbial communities.
High-throughput sequencing (HTS) has opened a new era for environmental microbial studies as large amounts of genetic information can be obtained without culturing. Some recent studies have used amplicon and shotgun metagenome pyrosequencing to characterize soil microbial communities worldwide [16–20]. These strategies have allowed a more exhaustive characterization of community patterns, composition and metabolic capabilities, and continue to change our understanding of the microbial world. To date, however, HTS approaches have not been employed in Argentina as a means to compare tillage systems and evaluate land use effects on soil microbial communities.
In this study, we examined the impact of agricultural management on soil microbial communities. To do so, we collected soil samples from sites under three different types of land use (conventional tillage, no till and no agriculture), at each of five different locations in the Argentine Pampas region (Figure 1). From these samples, we generated amplicon and shotgun metagenome libraries, which were subsequently sequenced using 454-FLX pyrosequencing. Together these data compose the designated PAMPA datasets.
Soil samples were obtained at five different sites in the Argentinean Pampas located in three isohyet regions (Figure 1): three production fields in the rolling pampas (La Estrella: LE, La Negrita: LN, Criadero Klein: CK, wet weather, 1,000 to 1,200 annual mm) and two experimental stations, at Balcarce (Ba, semi-wet, 800 to 1,000 annual mm) and Anguil (An, semi-arid, 600 to 800 annual mm). At each experimental station, soils were collected from three plots, with three different types of land use: conventional tillage (CT), no till (NT) and soils with no agricultural (NA) management. Bulk soil was obtained from all plots included in this study. In addition, wheat rhizospheric soil was also obtained from the Anguil CT and NT plots. Only one sampling campaign was performed at each site, except at the La Estrella production field in the rolling pampas where there were six sampling time points over a year. At least two independent soil samples from each plot and land use site were collected, resulting in a total of 30 samples for Anguil station, 20 for Balcarce station and 62 for the rolling pampas region (see Additional files 1 and 2 for a detailed description of sampling strategy and sample processing). Total DNA was prepared from all soil samples. In addition, total cDNA was also prepared from Anguil rhizospheric samples. Amplicon sequencing libraries were constructed by PCR amplification of the V4 variable region in the 16S rRNA gene. Shotgun metagenome libraries were also constructed from one genomic DNA (gDNA) (and one cDNA, when available) sample obtained from each plot (see Additional files 1 and 2 for further details). Amplicon and shotgun libraries were sequenced using 454-FLX-Titanium chemistry. Raw data processing was performed following standard procedures suggested by the manufacturer.
We obtained a total of 19,325,913 reads and 7,740,811,541 bases from 30 samples by 454-FLX shotgun metagenome sequencing and 1,051,470 16S ribosomal DNA and ribosomal RNA (rDNA/rRNA) reads from 126 samples by amplicon sequencing. The metatranscriptomic shotgun libraries were excluded from the analysis due to the low number of reads recovered after rRNA trimming (more than ten fold below other samples). The amplicon dataset was analyzed using QIIME v1.5 software package . Shotgun metagenome datasets were annotated by BLAST against the NCBI database and subsequent results imported into MEGAN  for further analysis. Numerical and statistical analyses were performed using the METAGENassist software  and the R packages ‘BiodiversityR’ and ‘Vegan’ (R Development Core Team) (see Additional file 1).
To rule out possible contaminants from non-microbial species, such as plant, human or any other allochthonous DNA, in our metagenome shotgun libraries, a taxonomy assignment of all reads was assessed. We performed peptide prediction using FragGeneScan  followed by BlastP annotation against the NCBI Database. The Blast output was analyzed using MEGAN . The results showed that 95% of the classified sequences were identified as Bacteria, 1% as Eukarya and 0.6% as Archaea, whereas the remaining 3.4% of sequences could not be classified above the cellular organism level (data not shown). Within the Eukarya, 42% of reads were classified as Viridiplantae (plants), 27% as Fungi, 12% as Metazoa, 6% as diatoms and 13% to other groups or could not be classified (data not shown). Plant sequences are likely to be from decomposing material. These results suggest that contamination with allochthonous DNA is minimal or nonexistent as we could not identify any genetic material from unexpected species in the soils (for example, humans).
We found that geographic-specific differences, possibly associated with water availability, were evident in the 16S rRNA amplicon analysis of 103 soil communities (23 samples were excluded from the preliminary analysis due to differences in sequencing depth and other biases, see Additional files 1 and 2). The semi-arid soils (An) harbored communities that clustered separately from the wet (LE, LN, CK) and semi-wet (Ba) soil microbial communities (analysis of similarity: ANOSIM = 0.672, P < 0.001, Figure 2A, Additional file 3: Figure S1). This observation could be explained by the very different environmental conditions in both areas: the eastern area (wet and semi-wet) is humid and fertile with fine-textured soils that are rich in organic matter, while the western area is semi-arid with shallow coarse-textured soils with low levels of organic matter. We used Bioenv analysis (see Additional file 1 for further details of the analysis) to test which soil properties best explained the variation in microbial community structure. We found that clay, organic matter content, pH and salinity were the most influential variables (Mantel test: r = 0.6209, P = 0.001).
Differences in microbial communities within the semi-arid region (An) were largely determined by soil source, that is rhizospheric compared to bulk soil (ANOSIM = 0.5614, P < 0.001, Figure 2A, Additional file 3: Figure S1). In addition, rhizospheric samples clustered separately depending on the type of genetic material amplified (ANOSIM = 0.5169, P = 0.001, Figure 2B, Additional file 3: Figure S1). At the DNA level, active, inactive and even dead microorganisms were detected, that is, all the microbes present in the sample. However, at the RNA level, only metabolically active microorganisms were detected due to their high rates of rRNA expression. Our results show that rhizospheric microbial signatures detected by 16S rDNA are clearly distinct from those detected by 16S rRNA, suggesting that bacterial activity was not necessarily correlated with bacterial abundance.
Land use was another important driver that defined microbial community assemblages. Bulk soil samples clustered separately depending on land use (ANOSIM: Anguil = 0.3954, P = 0.017; Balcarce = 0.3795, P = 0.001; rolling pampas = 0.2072, P = 0.01, Additional file 3: Figure S1). Moreover, samples collected from soils under different tillage systems at the two experimental stations (Ba, An) also clustered separately in the analysis (ANOSIM: Balcarce = 0.5476, P = 0.001; Anguil = 0.2652, P = 0.001, Additional file 3: Figure S1). These results suggest that different microbial communities were selected under each type of soil management.
The evaluation of metabolic categories using the shotgun metagenome libraries also showed that semi-arid western locations were different from wet and semi-wet eastern sites (ANOSIM = 0.2806, P < 0.001). Therefore, we propose that water availability is probably the primary driver that shapes microbial communities (Figure 2C, Additional file 3: Figure S1). There was also clear separation by soil source in western semi-arid samples (ANOSIM = 0.6688, P < 0.001, Figure 2C, Additional file 3: Figure S1). In addition, bulk soil samples clustered separately according to tillage system in An and Ba (ANOSIM: Balcarce = 0.5391, P = 0.01; Anguil = 0.2346, P = 0.02, Additional file 3: Figure S1). However, the latter observation was less defined for rhizospheric samples, suggesting that other conditions, such as plant phenotype and exudates, could determine bacterial populations in rhizospheric communities. The soil properties that best explained the functional variation between samples for shotgun sequencing analysis were silt, organic matter, nitrogen content, pH and salinity (Mantel test: r = 0.2771, P = 0.002).
Even though additional work is required, preliminary results indicated that differences in microbial communities were largely defined by the variables considered, for example, water availability, geographic location, soil source, genetic material amplified and land use or tillage system. However, this was not always observed at the functional metagenomic level, since some samples showed patterns different from those in amplicon analysis (Additional file 3: Figure S1). Differences between the amplicon and shotgun analyses could be due to the fact that the 16S rDNA/rRNA operational taxonomic unit (OTU) analysis was performed by clustering sequences based on similarity, while the metagenomic analysis was based on sequence annotation, constrained by SEED database size, its limited number of categories and their ambiguity in sequence identity. Nevertheless, we could not rule out the possibility that very different microbial species could have similar metabolisms, thus minimizing the differences at metabolic level.
The present project represents the first large-scale metagenomic study of soils in Argentina that explores the link between agricultural management and soil microbiome. The resulting PAMPA datasets are among the few available soil metagenomic datasets based on high-throughput sequencing  and, here we presented a preliminary analysis of our data. While more detailed analysis will be needed to test the ideas presented in this paper, results so far have shed considerable light on the largely unknown soil micro-ecosystem of the Argentine Pampas. We showed that the soil microbiome changes primarily because of water availability and agricultural land use, and that these changes are also linked to different tillage systems (no-till or conventional tillage).
Additional analysis of the PAMPA datasets will continue to expand our knowledge of soil microbiome composition and function. Future efforts will be directed at identifying particular species and metabolisms associated with each tillage system in each geographic region and enriched by the rhizosphere. In addition, the PAMPA datasets can also be used in future worldwide soil metagenomic projects for comparative purposes. Additional experimental and sequencing efforts will be needed to describe in detail the root-associated microorganisms for different crops in different conditions. Understanding soil microbial dynamics and identifying specific plant-interacting microbes will be important steps towards improving current agricultural and soil sustainability practices.
Availability of supporting data
All data are publicly available and can be accessed through the Bioproject PRJNA178180 or directly by the NCBI Sequence Read Archive (SRA) under the accession numbers SRA058523 and SRA056866 (Additional file 2: Table S1 for detailed information). Additional information to that presented in this paper will be available from the Soil Genetic Network (SoilGeNe) website .
analysis of similarity
no till farming
operational taxonomic unit
principal component analysis
polymerase chain reaction
Satorre EH, Slafer GA, In Wheat: Wheat production systems of the Pampas. 1999, New York: The Haworth Press, Inc: Ecology and Physiology of Yield Determination, 333-348.
FAOSTAT: 2013, [http://faostat.fao.org/] data sourced January 2013
Kovalevski E, García F: Norte de argentina bajo siembra directa. 2007, 1: 159-172.
Phillips RE, Thomas GW, Blevins RL, Frye WW, Phillips SH: No-tillage agriculture. Science. 1980, 208: 1108-1113. 10.1126/science.208.4448.1108.
Gebhardt MR, Daniel TC, Schweizer EE, Allmaras RR: Conservation tillage. Science. 1985, 230: 625-630. 10.1126/science.230.4726.625.
Montgomery DR: Soil erosion and agricultural sustainability. Proc Nat Acad Sci USA. 2007, 104: 13268-13272. 10.1073/pnas.0611508104.
Alvarez R, Santanatoglia OJ, García R: Soil respiration and carbon inputs from crops in a wheat-soyabean rotation under different tillage systems. Soil Use Manage. 1995, 11: 45-50. 10.1111/j.1475-2743.1995.tb00495.x.
Gomez E, Bisaro V, Conti M: Potential C-source utilization patterns of bacterial communities as influenced by clearing and land use in a vertic soil of Argentina. Applied Soil Ecol. 2000, 15: 273-281. 10.1016/S0929-1393(00)00078-0.
Aon M, Cabello M, Sarena D, Colaneri A, Franco M, Burgos J, Cortassa S: I. Spatio-temporal patterns of soil microbial and enzymatic activities in an agricultural soil. Applied Soil Ecol. 2001, 18: 239-254. 10.1016/S0929-1393(01)00153-6.
Nievas F, Bogino P, Nocelli N, Giordano W: Genotypic analysis of isolated peanut-nodulating rhizobial strains reveals differences among populations obtained from soils with different cropping histories. Applied Soil Ecol. 2012, 53: 74-82.
Agaras B, Wall LG, Valverde C: Specific enumeration and analysis of the community structure of culturable pseudomonads in agricultural soils under no-till management in Argentina. Applied Soil Ecol. 2012, 61: 305-319.
Zabaloy MC, Garland JL, Gómez MA: An integrated approach to evaluate the impacts of the herbicides glyphosate, 2,4-D and metsulfuron-methyl on soil microbial communities in the Pampas region, Argentina. Applied Soil Ecol. 2008, 40: 1-12. 10.1016/j.apsoil.2008.02.004.
Zabaloy MC, Gómez E, Garland JL, Gómez MA: Assessment of microbial community function and structure in soil microcosms exposed to glyphosate. Applied Soil Ecol. 2012, 61: 333-339.
Smalla K, Oros-Sichler M, Milling A, Heuer H, Baumgarte S, Becker R, Neuber G, Kropf S, Ulrich A, Tebbe CC: Bacterial diversity of soils assessed by DGGE, T-RFLP and SSCP fingerprints of PCR-amplified 16S rRNA gene fragments: do the different methods provide similar results?. J Microbiol Methods. 2007, 69: 470-479. 10.1016/j.mimet.2007.02.014.
Muyzer G, Smalla K: Application of denaturing gradient gel electrophoresis (DGGE) and temperature gradient gel electrophoresis (TGGE) in microbial ecology. Antonie van Leeuwenhoek. 1998, 73: 127-141. 10.1023/A:1000669317571.
Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, Arrieta JM, Herndl GJ: Microbial diversity in the deep sea and the underexplored ‘rare biosphere’. Proc Nat Acad Sci USA. 2006, 103: 12115-12120. 10.1073/pnas.0605127103.
Vogel TM, Simonet P, Jansson JK, Hirsch PR, Tiedje JM, Van Elsas JD, Bailey MJ, Nalin R, Philippot L: TerraGenome: a consortium for the sequencing of a soil metagenome. Nat Rev Microbiol. 2009, 7: 252-252. 10.1038/nrmicro2119.
Yin C, Jones KL, Peterson DE, Garrett KA, Hulbert SH, Paulitz TC: Members of soil bacterial communities sensitive to tillage and crop rotation. Soil Biol Biochem. 2010, 42: 2111-2118. 10.1016/j.soilbio.2010.08.006.
Roesch LFW, Fulthorpe RR, Riva A, Casella G, Hadwin AKM, Kent AD, Daroub SH, Camargo FAO, Farmerie WG, Triplett EW: Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 2007, 1: 283-290.
Delmont TO, Prestat E, Keegan KP, Faubladier M, Clark IM, Pelletier E, Hirsch PR, Meyer F, Gilbert JA, Le Paslier D, Simonet P, Vogel TM: Structure, fluctuation and magnitude of a natural grassland soil metagenome. ISME J. 2012, 6: 1677-1687. 10.1038/ismej.2011.197.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Peña AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nature Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303.
Huson DH, Auch AF, Qi J, Schuster SC: MEGAN analysis of metagenomic data. Genome Res. 2007, 17: 377-386. 10.1101/gr.5969107.
Arndt D, Xia J, Liu Y, Zhou Y, Guo AC, Cruz JA, Sinelnikov I, Budwill K, Nesbø CL, Wishart DS: METAGENassist: a comprehensive web server for comparative metagenomics. Nucleic Acids Res. 2012, 40: W88-W95. 10.1093/nar/gks497.
Rho M, Tang H, Ye Y: FragGeneScan: predicting genes in short and error-prone reads. Nucleic Acids Res. 2010, 38: e191-10.1093/nar/gkq747.
Soil Genetic Network: [http://soilgene.net]
Funding for this work was provided by the Agencia Nacional de Promoción Científica y Tecnológica of Argentina through PAE37164, “Asociación para el estudio y desarrollo de la Biología del Suelo”. We acknowledge Alejandro Arrieta, Franco Puccio, Florencia Barbarich, Gonzalo Berhongaray, Marcelo Soria and Olga Correa for their help during sampling and for helpful discussions and advice. We also would like to thank Fernando Lopez for the administrative work that made this project possible.
The authors declare that they have no competing interests.
NR, BC and RA did the sampling. NR, BC and MR processed the samples in the lab and did the sequencing using the 454-FLX. SR performed the bioinformatic processing of sequence data. NR did the analysis, generated preliminary results and uploaded the sequences to the Sequence Read Archive (SRA). BC did some of the statistical analysis. MPV supervised all the work done. NR, BC and MPV participated in the writing of this manuscript. AM supervised all groups involved in this project. All authors participated in the experimental design, discussions about data interpretation and manuscript supervision. All authors read and approved the final manuscript.