Pyrosequencing Investigation into the Bacterial Community in Permafrost Soils along the China-Russia Crude Oil Pipeline (CRCOP)

The China-Russia Crude Oil Pipeline (CRCOP) goes through 441 km permafrost soils in northeastern China. The bioremediation in case of oil spills is a major concern. So far, little is known about the indigenous bacteria inhabiting in the permafrost soils along the pipeline. A pilot 454 pyrosequencing analysis on the communities from four selected sites which possess high environment risk along the CRCOP is herein presented. The results reveal an immense bacterial diversity than previously anticipated. A total of 14448 OTUs with 84834 reads are identified, which could be assigned into 39 different phyla, and 223 families or 386 genera. Only five phyla sustain a mean OTU abundance more than 5% in all the samples, but they altogether account for 85.08% of total reads. Proteobacteria accounts for 41.65% of the total OTUs or 45% of the reads across all samples, and its proportion generally increases with soil depth, but OTUs numerically decline. Among Proteobacteria, the abundance of Beta-, Alpha-, Delta- and Gamma- subdivisions average to 38.7% (2331 OTUs), 37.5% (2257 OTUs), 10.35% (616 OTUs), and 6.21% (374 OTUs), respectively. Acidobacteria (esp. Acidobacteriaceae), Actinobacteria (esp. Intrasporangiaceae), Bacteroidetes (esp. Sphingobacteria and Flavobacteria) and Chloroflexi (esp. Anaerolineaceae) are also very common, accounting for 8.56% (1237 OTUs), 7.86% (1136 OTUs); 7.35% (1063 OTUs) and 8.27% (1195 OTUs) of total libraries, respectively. The ordination analysis indicates that bacteria communities in the upper active layer cluster together (similar), while bacterial consortia from the lower active layer and permafrost table scatter (less similar). The abundance of Rhodococcus (12 OTUs), Pseudomonas (71 OTUs) and Sphingomonas (87 OTUs) is even less (<0.01%). This effort to profile the background diversity may set the first stage for better evaluating the bacterial dynamics in response to accidental oil spills.     

石油集输系统中微生物群落结构研究

采用16SrRNA基因克隆一变性梯度凝胶电泳分析方法研究了石油集输系统原油和油田产水中的微生物群落结构。变性梯度凝胶电泳图谱显示：油田产水中微生物群落远比原油中的菌群丰富。所有的油田水样和原油样本中都存在与Ochrobactrum sp．和Stenotrophomonas sp．相关的细菌;原油样本中检测出与Burkholderia sp．、Brevundimonas sp．和Propionibacterium sp．相关的细菌,而这些细菌在油田水样中未检出;在油田水样中检出与Hippea sp．、Acidovorax sp．、Arcobacter sp．、Pseudomonas sp．、Thiomicrospira sp．、Brevibacterium sp．、Tissierella sp．和Peptostreptococcus sp．相关的细菌,而这些细菌在原油样本中未检出。用古细菌特异性引物进行检测发现在油田水样中存在与Methanomicrobials和Methanosarcinales相关的产甲烷菌,而这些细菌在原油样本中未检出。在石油集输过程中,油田水样和原油中微生物群落的相似性分别为83．3％和88．2％,说明微生物群落结构较为稳定。     

The Dynamic Change of Microbial Communities in Crude Oil-Contaminated Soils from Oil Fields in China

To study the biodegradability of microbial communities in crude oil contamination, crude oil-contaminated soil samples from different areas of China were collected. Using polyphasic approach, this study explored the dynamic change of the microbial communities during natural accumulation in oil field and how the constructed bioremediation systems reshape the composition of microbial communities. The abundance of oil-degrading microbes was highest when oil content was 3–8%. This oil content is potentially optimal for oil degrading bacteria proliferation. During a ～12 months natural accumulation, the quantity of oil-degrading microbes increased from 10⁵ to 10⁸ cells/g of soil. A typical sample of Liaohe (LH, oil-contaminated site near Liaohe River, Liaoning Province, China) was remediated for 50 days to investigate the dynamic change of microbial communities. The average FDA (a fluorescein diacetate approach) activities reached 0.25 abs/hr·g dry soil in the artificially enhanced repair system, 32% higher than the 0.19 abs/hr·g dry soil in natural circumstances. The abundance of oil-degrading microbes increased steadily from 0.001 to 0.068. During remediation treatment, oil content in the soil sample was reduced from 6.0% to 3.7%. GC–MS analysis indicated up to 67% utilization of C₁₀–C₂₀ normal paraffin hydrocarbons, the typical compounds that undergo microbial degradation.     

Phylogenetic Comparisons of Bacterial Communities from Serpentine and Nonserpentine Soils

I present the results of a culture-independent survey of soil bacterial communities from serpentine soils and adjacent nonserpentine comparator soils using a variety of newly developed phylogenetically based statistical tools. The study design included site-based replication of the serpentine-to-nonserpentine community comparison over a regional scale (~100 km) in Northern California and Southern Oregon by producing 16S rRNA clone libraries from pairs of samples taken on either side of the serepentine-nonserpentine edaphic boundary at three geographical sites. At the division level, the serpentine and nonserpentine communities were similar to each other and to previous data from forest soils. Comparisons of both richness and Shannon diversity produced no significant differences between any of the libraries, but the vast majority of phylogenetically based tests were significant, even with only 50 sequences per library. These results suggest that most samples were distinct, consisting of a collection of lineages generally not found in other samples. The pattern of results showed that serpentine communities tended to be more similar to each other than they were to nonserpentine communities, and these differences were at a lower taxonomic scale. Comparisons of two nonserpentine communities generally showed differences, and some results suggest that the geographical site may control community composition as well. These results show the power of phylogenetic tests to discern differences between 16S rRNA libraries compared to tests that discard DNA data to bin sequences into operational taxonomic units, and they stress the importance of replication at larger scales for inferences regarding microbial biogeography.     

Prepared Bed Land Treatment of Soils Containing Diesel and Crude Oil Hydrocarbons

An ex situ, field-scale, prepared bed land treatment unit (LTU) was used to bio-remediate soils containing petroleum hydrocarbons. Two soils were treated in side-by-side units to compare performance: (1) a clayey silt containing crude oil hydrocarbons from releases 30 to 40 years ago and (2) a silty sand containing diesel fuel hydrocarbons from a leak about three years prior to the bioremediation. The effectiveness of the bioremediation in the LTU was evaluated over a period of 18 months. The results indicated that: (1) prepared bed bioremediation reduced the hydrocarbon concentration, mobility, and relative toxicity in the soil with the diesel fuel, and (2) chemical bioavailability appeared to limit bioremediation of the soil containing the crude oil hydrocarbons. Although the soils containing the crude oil hydrocarbons contained an average of 10,000?mg TPH/kg dry soil, these soils had limited hydrocarbon availability, nontoxic conditions, and low potential for chemical migration. For the soils containing the diesel fuel, active prepared bed bioremediation of about 15 weeks was adequate to reach an environmentally acceptable endpoint. At that time, there was little further TPH loss, no Microtox^TM toxicity, and limited hydrocarbon mobility.     

Soil Type Is the Primary Determinant of the Composition of the Total and Active Bacterial Communities in Arable Soils

Degradation of agricultural land and the resulting loss of soil biodiversity and productivity are of great concern. Land-use management practices can be used to ameliorate such degradation. The soil bacterial communities at three separate arable farms in eastern England, with different farm management practices, were investigated by using a polyphasic approach combining traditional soil analyses, physiological analysis, and nucleic acid profiling. Organic farming did not necessarily result in elevated organic matter levels; instead, a strong association with increased nitrate availability was apparent. Ordination of the physiological (BIOLOG) data separated the soil bacterial communities into two clusters, determined by soil type. Denaturing gradient gel electrophoresis and terminal restriction fragment length polymorphism analyses of 16S ribosomal DNA identified three bacterial communities largely on the basis of soil type but with discrimination for pea cropping. Five fields from geographically distinct soils, with different cropping regimens, produced highly similar profiles. The active communities (16S rRNA) were further discriminated by farm location and, to some degree, by land-use practices. The results of this investigation indicated that soil type was the key factor determining bacterial community composition in these arable soils. Leguminous crops on particular soil types had a positive effect upon organic matter levels and resulted in small changes in the active bacterial population. The active population was therefore more indicative of short-term management changes.     

Diversity and Structure of Bacterial Chemolithotrophic Communities in Pine Forest and Agroecosystem Soils

Obligate lithotrophs (e.g., ammonia oxidizers) and facultative lithotrophs (e.g., CO and hydrogen oxidizers) collectively comprise a phylogenetically diverse functional group that contributes significantly to carbon and nitrogen cycles in soils and plays important roles in trace gas dynamics (e.g., carbon monoxide and nitrous and nitric oxides) that affect tropospheric chemistry and radiative forcing. In spite of their diverse physiologies, facultative and obligate lithotrophs typically possess the Calvin-Benson-Bassham cycle enzyme, ribulose-1,5-bisphosphate carboxylase/oxygenase (rubisCO). In an effort designed to understand the structure of lithotrophic communities in soil, genomic DNA extracts from surface (0 to 2 cm) and subsurface (5 to 7 cm) soils have been obtained from two sites in a Georgia agroecosystem (peanut and cotton plots) and an unmanaged pine stand (>50 years old). The extracts have been used in PCR amplifications of the cbbL gene for the rubisCO large subunit protein. cbbL PCR products were cloned, sequenced, and subjected to phylogenetic and statistical analyses. Numerous novel lineages affiliated with the form IC clade (one of four form I rubisCO clades), which is typified by facultative lithotrophs, comprised lithotrophic communities from all soils. One of the form IC clone sequences clustered with a form IC clade of ammonia-oxidizing Nitrosospira. Distinct assemblages were obtained from each of the sites and from surface and subsurface soils. The results suggest that lithotrophic populations respond differentially to plant type and land use, perhaps forming characteristic associations. The paucity of clone sequences attributed to ammonia-oxidizing bacteria indicates that even though ammonia oxidation occurs in the various soils, the relevant populations are small compared to those of facultative lithotrophs.     

Spatial and Temporal Variation of Archaeal,Bacterial and Fungal Communities in Agricultural Soils

Background

Soil microbial communities are in constant change at many different temporal and spatial scales. However, the importance of these changes to the turnover of the soil microbial communities has been rarely studied simultaneously in space and time.

Methodology/Principal Findings

In this study, we explored the temporal and spatial responses of soil bacterial, archaeal and fungal β-diversities to abiotic parameters. Taking into account data from a 3-year sampling period, we analyzed the abundances and community structures of Archaea, Bacteria and Fungi along with key soil chemical parameters. We questioned how these abiotic variables influence the turnover of bacterial, archaeal and fungal communities and how they impact the long-term patterns of changes of the aforementioned soil communities. Interestingly, we found that the bacterial and fungal β-diversities are quite stable over time, whereas archaeal diversity showed significantly higher fluctuations. These fluctuations were reflected in temporal turnover caused by soil management through addition of N-fertilizers.

Conclusions

Our study showed that management practices applied to agricultural soils might not significantly affect the bacterial and fungal communities, but cause slow and long-term changes in the abundance and structure of the archaeal community. Moreover, the results suggest that, to different extents, abiotic and biotic factors determine the community assembly of archaeal, bacterial and fungal communities.     

Leaching of BTEX from Aged Crude Oil Contaminated Model Soils: Experimental and Modeling Results

It is generally assumed that soil properties such as organic matter content, porosity, and mineral surface area have a significant effect on the bioavailability and leachability of aged petroleum hydrocarbons. In order to test this hypothesis, nine model soils or sorbents (i.e., fine and coarse quartz sand, montmorillonite and kaolinite clay, peat, 60Å and 150Å silica gel, a loam soil, and non-porous glass beads) were spiked with a crude oil, aged for 27 months in the laboratory, and transferred to glass columns for the performance of continuous flow leaching experiments. The column effluents were periodically sampled for 43 days and analyzed for BTEX. A one-dimensional flow model for predicting the dissolution and dispersion of individual hydrocarbons from a multi-component NAPL such as crude oil was used to fit the leaching data (i.e., the BTEX concentration versus time curves) by adjusting the equilibrium oil-leachate partitioning coefficient (K ^ol ) for each respective hydrocarbon. The Peclet number, which is a measure of dispersion and a required modeling parameter, was measured in separate chloride tracer experiments for each soil column.

Results demonstrate that soil properties did not significantly affect the leaching kinetics of BTEX from the columns. Instead, BTEX leaching curves could be successfully fitted with the one-dimensional NAPL dissolution flow model for all sorbents with the exception of montmorillonite clay. The fitting parameter K ^ol for each hydrocarbon was found to be similar to the K ^ol values that were independently measured for the same crude oil by Rixey et al. (Journal of Hazardous Materials B, 65: 137–156, 1999 Rixey, W. G., Garg, S. and Nie, Y. 1999. Comparison of the fixed-bed and batch leaching characteristics of aromatic compounds in residually trapped crude oils and oily wastes. J. Hazard. Mat. B, 64: 137–156. [CSA][CROSSREF] [Google Scholar]). In addition, the fitted K ^ol values were very similar for BTEX leaching from aged compared to freshly spiked loam soil. These findings indicate that leaching of BTEX in the aged soils that are contaminated with crude oil at the high concentrations commonly found in the environment (i.e., > 20,000 mg/kg) was not affected by soil properties or aging but rather was governed by the equilibrium dissolution of these hydrocarbons from the crude oil NAPL that is coating the soil particles.     

Assessment of Bacterial Communities in Auriferous and Non-Auriferous Soils Using Genetic and Functional Fingerprinting

Understanding the microbial processes affecting the mobility of Au is important in the development of biogeochemical models describing the formation of secondary anomalies and Au grains in soils and deeper regolith materials. This study characterizes bacterial activity in auriferous soils that is linked to the microbially mediated solubilization of Au, as a result of production and consumption of free amino acids, which can form stable complexes with Au. Through the application of 16S rDNA fingerprinting and community level physiological profiling (CLPP), concurrently with Au mobility data, microcosm experiments have demonstrated the role that mobile Au plays in determining the structure and function of bacterial communities in auriferous soils. The bacterial community of auriferous soils displayed genetic differences compared to non-auriferous (background) soils associated with the appearance of Methylocella sp., Arthrobacter sp. and Bacillus sp., as well as functional differences in the utilization of D-Cellobiose, L-Serine, L-Phenylalanine, L-Arginine and N-Acetyl-D-Glucosamine. These results suggest that soil bacterial communities are linked to biogeochemical Au cycling, and that microbial fingerprinting analyses may be used as a screening tool in Au exploration to differentiate auriferous from background terrains.     

Response of Antarctic Soil Bacterial Assemblages to Contamination by Diesel Fuel and Crude Oil

Abstract The effects of diesel fuel and ``Arabian light' crude oil addition on Antarctic bacterial assemblages were studied in four contaminated soils during 1 year in the Terre Adelie land area. Monthly sampling allowed a regular survey of the bacterial changes occurring in the contaminated soils. All samples were analyzed for total bacteria, heterotrophic culturable microbiota, and hydrocarbon-utilizing microbiota. Crude oil contamination induced an initial increase of all bacterial parameters in all contaminated soils. Diesel oil contamination had a more complex effect. Hydrocarbon degrading bacterial abundance increases occurred after diesel oil addition. In contrast, general heterotrophic bacterial abundance could significantly decrease in the same conditions. In all cases the stimulatory effects of oil addition disappeared after several months of contamination. Received: 13 April 1999; Accepted: 24 February 2000; Online Publication: 29 May 2000     

Detection of a Reproducible,Single-Member Shift in Soil Bacterial Communities Exposed to Low Levels of Hydrogen

Soil is exposed to hydrogen when symbiotic rhizobia in legume root nodules cannot recycle the hydrogen that is generated during nitrogen fixation. The hydrogen emitted is most likely taken up by free-living soil bacteria that use hydrogen as an energy source, though the bacteria that do this in situ remain unclear. In this study, we investigated the effect of hydrogen exposure on the bacteria of two different soils in a microcosm setup designed to simulate hydrogen-emitting root nodules. Although the size and overall composition of the soil bacterial community did not significantly alter after hydrogen exposure, one ribotype increased in relative abundance within each soil. This single-ribotype shift was identified by generating multiple terminal restriction fragment length polymorphism (T-RFLP) profiles of 16S rRNA genes from each soil sample, with gene sequence confirmation to identify terminal restriction fragments. The increased abundance of a single ribotype after hydrogen exposure, within an otherwise similar community, was found in replicate samples taken from each microcosm and was reproducible across replicate experiments. Similarly, only one member of the soil bacterial community increased in abundance in response to hydrogen exposure in soil surrounding the root nodules of field-grown soybean (Glycine max). The ribotypes that increased after hydrogen exposure in each soil system tested were all from known hydrogen-oxidizing lineages within the order Actinomycetales. We suggest that soil actinomycetes are important utilizers of hydrogen at relevant concentrations in soil and could be key contributors to soil''s function as a sink in the global hydrogen cycle.Soil is the major sink in the global hydrogen cycle and accounts for approximately 75 to 80% of uptake from the atmosphere (7, 11). Soil is such a strong sink that the atmospheric mixing ratio of molecular hydrogen, H₂, is hemispherically asymmetric because of the greater landmass in the Northern Hemisphere (11). Many nitrogen-fixing bacteria that form symbiotic relationships with legume plants cannot recycle the H₂ that is generated during N₂ fixation (2, 13). Most of the H₂ emitted from legume root nodules is taken up by the surrounding soil, within a few centimeters of the nodule surface, and is not released to the atmosphere (20). Although the H₂ emitted by the rhizobial symbionts costs the legume approximately 5% of its daily photosynthate and “represents a significant investment by the plant” (9), there is growing evidence to suggest that soil exposed to H₂ is beneficial to plant growth, separate from the benefits derived from N₂ fixation (8, 10, 28). Previously, La Favre and Focht have hypothesized that “the hydrogen which is evolved during N₂ fixation represents an additional energy input into the plant-soil ecosystem… since metabolism of H₂ by chemolithotrophic bacteria results in an input of fixed carbon to the system” (20). A number of studies have found that when H₂ is taken up by soil, net CO₂ fixation occurs at the rate of 7 to 8 nmol CO₂ per g of soil per h (22, 34). For a legume fixing 200 kg of atmospheric N₂ per hectare, over 200,000 liters of H₂ could be released into the legume''s rhizosphere over the duration of the growing season and CO₂ fixation could result in an extra 25 kg of soil carbon fixed per hectare (9, 10, 28).Many bacteria isolated from soil are able to utilize H₂ as an energy source (2, 5-7, 21), and free-living bacteria are most likely responsible for the H₂ uptake observed by soil surrounding legume roots (22). Adding a bacterial energy source, such as H₂, could affect the microbial population size, as has been observed previously (34), but more specific shifts within the bacterial community may occur if just the microorganisms able to utilize the energy source multiply. Their activity could also have downstream consequences specifically on other members of the community. Most H₂-oxidizing cultures have required enrichment with concentrations of H₂ that are not environmentally relevant and therefore cannot be assumed to be carrying out H₂ oxidation at much lower, naturally occurring concentrations (5-7). Recent surveys of microbes present in soil samples, via their nucleic acids, have revealed many novel bacterial inhabitants that have been little studied and thus may also be contributing to the repertoire of bacterial soil processes, such as H₂ uptake (16). A recent study into the effect of H₂ on soil bacteria focused on a few groups of H₂-oxidizing, autotrophic bacteria and thus ignored many other H₂ utilizers potentially present in soil (34).There are now many ways of characterizing the entire microbial community in environmental samples, either via their entire genomic content, though metagenomic analysis of soil is difficult at present, or via analysis of the lineages present according to 16S rRNA gene sequences, or ribotypes (36). A recent study comparing high-throughput pyrosequencing of 16S rRNA genes and an easily accessible profiling method, known as terminal restriction fragment length polymorphism (T-RFLP), found the simpler profiles were appropriate for comparing the dominant ribotypes in multiple samples (24). Although T-RFLP profiles only provide a simplified snapshot of the dominant members in microbial communities, compared to the deeper analyses provided by microarrays or high-throughput sequencing technologies, T-RFLP profiling is a cost-effective, reproducible, and robust method of “fingerprinting” many soil samples rapidly and efficiently (14, 24, 25, 32).In this study, we chose to assess the dominant members of the soil bacterial community via T-RFLP profiles of ribotypes present in H₂-treated and control soils to avoid a narrow focus on well-studied H₂ oxidizers. We investigated the bacterial community structure in two different soils, utilizing a microcosm setup with concentrations of H₂ calculated to occur in the rhizosphere of N₂-fixing legumes, to determine whether common responses to H₂ exposure could be predicted from soils that differ by climate, edaphic characteristics, and starting communities. Soil in microcosms has previously been shown to have similar H₂ uptake properties to soil close to H₂-emitting legume nodules (9), but we also complemented our plant-free microcosm work with an examination of soil collected from the root systems of field-grown soybean (Glycine max (L.) Merr.).     

Relative Roles of Deterministic and Stochastic Processes in Driving the Vertical Distribution of Bacterial Communities in a Permafrost Core from the Qinghai-Tibet Plateau,China

Understanding the processes that influence the structure of biotic communities is one of the major ecological topics, and both stochastic and deterministic processes are expected to be at work simultaneously in most communities. Here, we investigated the vertical distribution patterns of bacterial communities in a 10-m-long soil core taken within permafrost of the Qinghai-Tibet Plateau. To get a better understanding of the forces that govern these patterns, we examined the diversity and structure of bacterial communities, and the change in community composition along the vertical distance (spatial turnover) from both taxonomic and phylogenetic perspectives. Measures of taxonomic and phylogenetic beta diversity revealed that bacterial community composition changed continuously along the soil core, and showed a vertical distance-decay relationship. Multiple stepwise regression analysis suggested that bacterial alpha diversity and phylogenetic structure were strongly correlated with soil conductivity and pH but weakly correlated with depth. There was evidence that deterministic and stochastic processes collectively drived bacterial vertically-structured pattern. Bacterial communities in five soil horizons (two originated from the active layer and three from permafrost) of the permafrost core were phylogenetically random, indicator of stochastic processes. However, we found a stronger effect of deterministic processes related to soil pH, conductivity, and organic carbon content that were structuring the bacterial communities. We therefore conclude that the vertical distribution of bacterial communities was governed primarily by deterministic ecological selection, although stochastic processes were also at work. Furthermore, the strong impact of environmental conditions (for example, soil physicochemical parameters and seasonal freeze-thaw cycles) on these communities underlines the sensitivity of permafrost microorganisms to climate change and potentially subsequent permafrost thaw.     

Impact of Logging and Forest Conversion to Oil Palm Plantations on Soil Bacterial Communities in Borneo

Tropical forests are being rapidly altered by logging and cleared for agriculture. Understanding the effects of these land use changes on soil bacteria, which constitute a large proportion of total biodiversity and perform important ecosystem functions, is a major conservation frontier. Here we studied the effects of logging history and forest conversion to oil palm plantations in Sabah, Borneo, on the soil bacterial community. We used paired-end Illumina sequencing of the 16S rRNA gene, V3 region, to compare the bacterial communities in primary, once-logged, and twice-logged forest and land converted to oil palm plantations. Bacteria were grouped into operational taxonomic units (OTUs) at the 97% similarity level, and OTU richness and local-scale α-diversity showed no difference between the various forest types and oil palm plantations. Focusing on the turnover of bacteria across space, true β-diversity was higher in oil palm plantation soil than in forest soil, whereas community dissimilarity-based metrics of β-diversity were only marginally different between habitats, suggesting that at large scales, oil palm plantation soil could have higher overall γ-diversity than forest soil, driven by a slightly more heterogeneous community across space. Clearance of primary and logged forest for oil palm plantations did, however, significantly impact the composition of soil bacterial communities, reflecting in part the loss of some forest bacteria, whereas primary and logged forests did not differ in composition. Overall, our results suggest that the soil bacteria of tropical forest are to some extent resilient or resistant to logging but that the impacts of forest conversion to oil palm plantations are more severe.     

Responses of Bacterial Communities in Arable Soils in a Rice-Wheat Cropping System to Different Fertilizer Regimes and Sampling Times

Soil physicochemical properties, soil microbial biomass and bacterial community structures in a rice-wheat cropping system subjected to different fertilizer regimes were investigated in two seasons (June and October). All fertilizer regimes increased the soil microbial biomass carbon and nitrogen. Both fertilizer regime and time had a significant effect on soil physicochemical properties and bacterial community structure. The combined application of inorganic fertilizer and manure organic-inorganic fertilizer significantly enhanced the bacterial diversity in both seasons. The bacterial communities across all samples were dominated by Proteobacteria, Acidobacteria and Chloroflexi at the phylum level. Permutational multivariate analysis confirmed that both fertilizer treatment and season were significant factors in the variation of the composition of the bacterial community. Hierarchical cluster analysis based on Bray-Curtis distances further revealed that bacterial communities were separated primarily by season. The effect of fertilizer treatment is significant (P = 0.005) and accounts for 7.43% of the total variation in bacterial community. Soil nutrients (e.g., available K, total N, total P and organic matter) rather than pH showed significant correlation with the majority of abundant taxa. In conclusion, both fertilizer treatment and seasonal changes affect soil properties, microbial biomass and bacterial community structure. The application of NPK plus manure organic-inorganic fertilizer may be a sound fertilizer practice for sustainable food production.