Impacts of Long-Term Irrigation of Domestic Treated Wastewater on Soil Biogeochemistry and Bacterial Community Structure

Freshwater scarcity and regulations on wastewater disposal have necessitated the reuse of treated wastewater (TWW) for soil irrigation, which has several environmental and economic benefits. However, TWW irrigation can cause nutrient loading to the receiving environments. We assessed bacterial community structure and associated biogeochemical changes in soil plots irrigated with nitrate-rich TWW(referred to as pivots) for periods ranging from 13 to 30 years. Soil cores (0 to 40 cm) were collected in summer and winter from five irrigated pivots and three adjacently located nonirrigated plots. Total bacterial and denitrifier gene abundances were estimated by quantitative PCR (qPCR), and community structure was assessed by 454 massively parallel tag sequencing (MPTS) of small-subunit (SSU) rRNA genes along with terminal restriction fragment length polymorphism (T-RFLP) analysis of nirK, nirS, and nosZ functional genes responsible for denitrification of the TWW-associated nitrate. Soil physicochemical analyses showed that, regardless of the seasons, pH and moisture contents (MC) were higher in the irrigated (IR) pivots than in the nonirrigated (NIR) plots; organic matter (OM) and microbial biomass carbon (MBC) were higher as a function of season but not of irrigation treatment. MPTS analysis showed that TWW loading resulted in the following: (i) an increase in the relative abundance of Proteobacteria, especially Betaproteobacteria and Gammaproteobacteria; (ii) a decrease in the relative abundance of Actinobacteria; (iii) shifts in the communities of acidobacterial groups, along with a shift in the nirK and nirS denitrifier guilds as shown by T-RFLP analysis. Additionally, bacterial biomass estimated by genus/group-specific real-time qPCR analyses revealed that higher numbers of total bacteria, Acidobacteria, Actinobacteria, Alphaproteobacteria, and the nirS denitrifier guilds were present in the IR pivots than in the NIR plots. Identification of the nirK-containing microbiota as a proxy for the denitrifier community indicated that bacteria belonged to alphaproteobacteria from the Rhizobiaceae family within the agroecosystem studied. Multivariate statistical analyses further confirmed some of the above soil physicochemical and bacterial community structure changes as a function of long-term TWW application within this agroecosystem.