Evaluating the potential impact of proton carriers on syntrophic propionate oxidation

Anaerobic propionic acid degradation relies on interspecies electron transfer (IET) between propionate oxidisers and electron acceptor microorganisms, via either molecular hydrogen, formate or direct transfers. We evaluated the possibility of stimulating direct IET, hence enhancing propionate oxidation, by increasing availability of proton carriers to decrease solution resistance and reduce pH gradients. Phosphate was used as a proton carrying anion, and chloride as control ion together with potassium as counter ion. Propionic acid consumption in anaerobic granules was assessed in a square factorial design with ratios (1:0, 2:1, 1:1, 1:2 and 0:1) of total phosphate (TP) to Cl−, at 1X, 10X, and 30X native conductivity (1.5 mS.cm−1). Maximum specific uptake rate, half saturation, and time delay were estimated using model-based analysis. Community profiles were analysed by fluorescent in situ hybridisation and 16S rRNA gene pyrosequencing. The strongest performance was at balanced (1:1) ratios at 10X conductivity where presumptive propionate oxidisers namely Syntrophobacter and Candidatus Cloacamonas were more abundant. There was a shift from Methanobacteriales at high phosphate, to Methanosaeta at low TP:Cl ratios and low conductivity. A lack of response to TP, and low percentage of presumptive electroactive organisms suggested that DIET was not favoured under the current experimental conditions.

In a methanogenic environment, in the absence of an external electron acceptor such as nitrate or sulfate, anaerobic oxidation of organic acids to acetate and hydrogen (H 2 ) is only thermodynamically feasible where microorganisms associate to allow electron sinking to a methanogen 1 . This process is known as syntrophy and involves the transfer of electrons from one microorganism (electron donor) to the other (electron acceptor) by interspecies electron transfer (IET) 2,3 .
Propionate is a key substrate oxidised via syntrophy under methanogenic conditions, with two distinct reactions combining to generate catabolic energy necessary for the survival of the partners involved. Known syntrophic propionate oxidisers include Syntrophobacter wolinii 4 , Syntrophobacter sulfatireducens 5 , Syntrophobacter fumaroxidans 6 , Pelotomaculum schinkii 7 and Smithella propionica 8 . Strain Candidatus Cloacamonas acidaminovorans has also been linked to syntrophic propionate oxidation 9,10 . The net microbial reaction is oxidation of propionate to acetate and hydrogen with energy generated through substrate level phosphorylation on the propionyl-CoA (Eq. 1) 1,11,12 . The butyrate is further oxidised to acetate and H 2 13,14 , which results in a net equation equivalent to Eq 1. In order to enable oxidation reactions in either case, electrons are utilised by methanogens to reduce carbon dioxide (CO 2 ) to methane (CH 4 ). Hydrogenotrophic methanogenesis reduces CO 2 to CH 4 . Most H 2 utilising methanogens do not have cytochromes 15,16 , which suggests that energy is generated mostly through a sodium (cation) motive force (SMF) 17 .
The most commonly proposed mechanism for electron transfer is mediated IET (MIET) with H 2 and/or formate as the electron carriers 18,19 . However other oxidative organisms fulfilling similar metabolic niches on other substrates have been observed to engage in direct electron transfer as detailed below.
Co-cultures of Geobacter sp engineered to study interspecies H 2 transfer have been shown to engage in direct IET instead which was metabolically more favourable, with growth observed within 21 days of incubation instead of 7 months 20 . DIET, via membrane-bound proteins (OmcS) 21 or through conductive pili, was observed in environmental samples during the degradation of ethanol where Geobacter sp, capable of DIET, were the dominant oxidisers and Methanosaeta sp the electron acceptor microorganism 21,22 .
DIET partnerships have not been widely observed in natural communities possibly because existing culture conditions favour MIET and difficulty in transferring DIET engaged cultures 23 . Advantages of DIET vs MIET are that molecular diffusional limitations (and gradients) are avoided and translation of electrons to a molecular carrier are not required. However, as a potential limitation, DIET involves the direct transport of electrons between cells and also requires the migration of a charged ions to maintain electroneutrality 24 . Charge can be transferred either as protons, or as saline anions and cations, and generally, the latter (anions and cations) will dominate, due to the very low concentrations of protons at neutral pH conditions. Migration of saline ions will result in a pH gradient, with a low pH at the DIET oxidiser, and high pH at the DIET reducer. This can result in inhibition by pH changes and decreased thermodynamic viability, due to charge accumulation. Under such conditions, increasing proton mobility outside the cells would increase proton availability.
Buffers increase proton transport rate by molecular diffusion of the buffer, particularly in its associated form (allowing pH gradient neutralisation). In addition, the presence of the buffer will mitigate pH gradients caused by transport of saline ions due to the buffering effect of the weak base. As early as 1979, Engasser and Wilhem-Despres proved the ability of buffers, at a higher concentration, to shuttle protons, which could help improve processes such as ion exchange 25 . Torres et al. (2008) have shown that phosphate buffer facilitates proton transfer from anode-respiring electroactive biofilm to bulk liquid in single chamber microbial fuel cells, translating into higher current densities, while Nagarajan et al. (2013) have predicted a direct influence of proton mobility on the energetic efficiency of DIET. It was noted that total salinity will also enhance DIET, through a decrease in solution resistance, but this will also subject the microbial community to saline inhibition 24,26,27 .
Therefore, the presence of buffers, particularly phosphate, which have pKa values near circumneutral values (6)(7)(8), should enhance DIET, by reducing pH gradients, and enhancing proton mobility. To test this hypothesis (that DIET and hence activity is promoted by proton motility), this paper evaluates propionate utilisation capacity under varying TP:Cl − (TP = total phosphate) ratio (5 levels) and varying total conductivity (3 levels) in a square factorial design. Potassium was used as counter ion to minimise inhibition effects. Experiments were designed to assess whether either solution salinity or proton motility separately have an influence on propionate utilisation rate, and by inference, direct interspecies electron transfer. The inoculum used consisted of anaerobic granules from anaerobic digester treating brewery waste, which are known to be dominated by IET 14 .

Results
Specific Methane Activity tests were initially carried out on intact anaerobic granules to establish the influence of crushing the granules on the overall experiment (Fig. S1). A small delay was noted at the onset of the experiment when the physical property was modified, but the exponential curves are comparable.
Since increase in biomass has previously been related to increase in microbial growth 28 , the changes in VS were used as a quantitative assumption for any increase or decrease in the microbial population (Table S2). With an initial VS value of 2.27 ± 0.04 g L −1 , an overall increase in biomass was noted after SMA in most experiment. Figure 2 reports the maximum specific rate (r max ) of methane production for each assay. The r max values increased with changes in TP:Cl − ratio, with the highest values of 4.24 (± 0.27) gCOD L −1 day −1 and 4.28 (± 0.25) gCOD L −1 day −1 found with ratios of 1:1 and 2:1 respectively. The highest r max of all the tested batches was at 10X conductivity with a 1:1 ratio, at 4.67 (± 0.28) gCOD L −1 day −1 . Lower values were recorded at higher Cl − concentration, for samples at 1:2 and 0:1 ratio. At 30X conductivity, the r max value decreased significantly across all the samples, with the highest value at only 2.70 (± 0.32) gCOD L −1 day −1 with equimolar PO 4 3− and Cl − (Fig. 2). According to correlation analysis the variation in r max across treatments was positively correlated with variability in conductivity (r = − 0.809, p = 0.008). High conductivities were also significantly correlated with time delay (r = 0.774, p = 0.001) and had a significant impact on changes in biomass (r = 0.543, p = 0.036) with, also, a negative correlation with the rate of the reaction (r = 0.682, p = 0.043).
Scientific RepoRts | 5:18364 | DOI: 10.1038/srep18364 Microbial community response. Both FISH and pyrosequencing performed on the initial inoculum showed that bacteria were relatively low in abundance, compared to the archaea (Fig. 3). Pyrosequencing results initially identified Syntrophomonas (4.6%), Bacteroidales (3.9%) and the class of SHA-114 (2.6%) as the most prevalent bacteria in the community dominated by three main genera of archaea as being Methanobacterium, Methanosaeta and Methanolinea (Fig. 4 & Table S3).
This work focused specifically on presumptive syntrophic propionate oxidisers. While, very few cells from the Syntrophobacter species were identified with FISH (determined with mixed Syn 835 with EUB mix-Cy3 in Fig. 3), the pyrosequencing analysis identified two known propionate oxidisers belonging to the Syntrophobacter genus and Candidatus Cloacamonas genus initially at 1.3% and 0.7% relative abundance, respectively. As seen in Fig. 4, while there were no obvious changes in the abundance of Syntrophobacter at 1X conductivity, an increase was observed at 10X conductivity with percentages of 3.1%, 3.6% and 2.6% at 1:0, 1:1 and 0:1 TP to Cl − ratio respectively. In comparison, Candidatus Cloacamonas was more susceptible to conditions with TP only and Cl − only. At higher conductivities (10X and 30X), presumptive propionate oxidisers (see above) were positively correlated with samples that had TP and Cl − ratio of 1:1 as further confirmed by PCA across PC2 (Fig. 5).
The only bacteria known to participate in DIET belonged to the Geobacteraceae family and were originally in quite low abundance (0.6%), with an increase in relative abundance noted in assays at 1X conductivity and a more substantial increase of 2.13% with Cl − only at 10X conductivity (Fig. 4).
Another feature demonstrated in Fig. 5 is that low conductivity samples, and high-conductivity chloride only samples cluster (negative on PC1), while higher conductivity samples with phosphate are more widely distributed and are generally positive in PC1. This is mainly driven by a shift from dominant OTU Methanobacteriales in the inoculum and high phosphate samples, towards dominant OTU Methanosaeta in low conductivity and low phosphate samples. Specifically, the population shift averaged 46% to 13% for Methanobacteriales, corresponding   to an average increase for Methanosaeta from 30% to 60%. Variation in high-phosphate samples is mainly caused by a variation in bacterial community.

Discussion
Propionic acid utilisation with consequent methane production is an indication of syntrophic IET. Kinetic parameters indicate that the best results were obtained with a balanced TP:Cl − ratio in all conductivities tested and that activity decreased as the balance shifted to phosphate or chloride, particularly at higher conductivities.
Overall, correlation analysis indicates that conductivity is the major factor influencing the system, including changes in community profile (also related to availability of phosphate), delay (which is an indicator of short term adaptation), and rate. The correlation analysis tests for linear relationships and indicates no linear relationship between TP:Cl − ratio and activity. However, results clearly suggest a non-linear relationship, with an optimum occurring at a balanced TP:Cl − ratio, suggesting that the reaction is most effective in the presence of both anions.
Despite previous research showing that low pH conditions can have a negative impact during anaerobic oxidation of organics 29,30 , in our study the system managed to recover and this was attributed to the presence of both  buffers and potassium ions (K + ). Specifically, buffers mitigate pH decreases caused by accumulation of organic acids, while potassium (after adaptation) enables a degree of pH tolerance 31 .
Mechanistically, an increase in buffer and/or salt concentration, results in an accumulation of K + , in order to maintain cell turgor pressure 32 . At high conductivity, in the presence of Na + , inhibition is stronger if K + is limiting. This is known as K + starvation and is translated in initial growth, reaching a quick stationary phase as the level of K + ions in the media decreases 32,33 , a condition shown in batches treated with NaCl ( Fig S2).
A more diverse microbial population was observed in those treatments that resulted in a higher growth rate, namely at 1:0 and 1:1 TP to Cl − ratio at 10X conductivity and 1:1 ratio at 30X conductivity. This was confirmed by both rarefaction measurement (Fig. S5) and the Shannon index (Table S4). Two populations of hydrogenotrophic methanogens belonging to the genus Methanobacterium and Methanolinea were identified. Our findings support the initial sensitivity of Methanobacterium species, particularly at higher Cl − to TP ratio, over the acetoclastic Methanosaeta and H 2 utilising Methanolinea, phylogenetically related to Methanolinea tarda (Fig. S3) of the Methanoregulaceae family, capable of accepting electrons from both H 2 and formate for methane production 34 .
The two main propionate oxidisers identified belonged to the genus Syntrophobacter and Candidatus Cloacamonas, with both abundant at 10X conductivity and 1:1 TP to Cl − ratio, where the highest r max was recorded. So far, neither has been linked to DIET. A literature search in the genome of these species was carried out. Since the genome of Syntrophobacter wolinii, phylogenetically linked to our Syntrophobacter cluster (Fig. S4), is not yet available, Syntrophobacter fumaroxidans 35 was used as the closest relative strain. While membrane-bound cytochrome involved in ion-translocating ferredoxin:NADH oxidoreductase (Sfum_2694-99) and other transmembrane cytochromes linked to proton transfer were characterised, no pili-type genes or cytochromes that could be linked to DIET was detected in the genome of Syntrophobacter fumaroxidans 35 . As for the Candidatus Cloacamonas species similar to the acidaminovorans strain (Fig. S4), potential links to IET have previously been drawn 13 . Its genome does contain putative membrane proteins and genes encoding for type IV pili such as PilB, PilC as well as PilT 13,36 but no link to DIET has been reported yet.
The key species that could potentially be involved in DIET belongs to the Geobacteraceae family and was identified as Geobacter hephaestius (Fig. S4). Members of this family can use a variety of carbon sources including propionate 37 . Growth of Geobacter was affected by increasing conductivity, with its abundance more significant at higher Cl − to TP ratio, a condition where both Syntrophobacter and Candidatus Cloacamonas sp were sensitive.
At this stage and under the current experimental conditions, the results obtained , including lack of response to improved cation motility (TP compared to Cl − controls), and no substantial population of known electroactive organisms, does not support the hypothesis that DIET is a mode for propionate oxidation in this community, or alternatively, that DIET is not controlled by proton mobility. Further investigation under acclimatised conditions and in a known electroactive system are both areas which could further clarify controlling mechanisms.

Methods
Inoculum. Anaerobic granules grown on high levels of organics, including propionic acid 38 , were obtained from an active anaerobic digester (Fosters Breweries, Yatala, Queensland). These were separated into 20 g lots and enclosed into clean zip lock bags and crushed using a mortar and pestle. Granules were crushed (10-100 microns) in order to better visualise communities during FISH and avoid artefacts due to the deep biofilms, and performance was investigated against intact granules. All tests were carried out in triplicate to allow for the calculation of 95% confidence intervals in means (two-tailed t-test).

Chemical analysis
Biogas analysis. Methane (CH 4 ), hydrogen (H 2 ) and carbon dioxide (CO 2 ) were analysed at regular intervals by Gas Chromatography (GC) analysis. Using a 10 ml glass syringe, 5 ml of sampled headspace was analysed in a gas chromatograph (Perkin Elmer, Autosystem) with a thermo conductivity detector (GC-TCD) at 100 °C. a 99.99% purity Nitrogen (N 2 ) gas was used as the carrier through a column (Alltech, 8011/2) of 50 m to 320 μ m diameter. Results are reported as gCOD L −1 of CH 4 produced over time. For every mole of propionic acid in the batch at initial conditions 2.75 mole of CH 4 is expected with 0.75 moles produced from hydrogenotrophic methanogenesis and 2 moles from acetoclastic methanogenesis (derived from Eq. 1).
VFA analysis. Liquid samples, (3-4 ml), were taken at the same intervals as the gas samples and were filter-sterilised using 0.22 μ m pore size Millipore TM membrane cartridge filters. A 1:4 dilution was applied to the first 4 samples while the other samples remained undiluted. A 1% v/v of sodium azide (10% stock solution) was added to each vial before analysis by an HP Agilent high pressure liquid chromatography (HPLC) system, adapted with an inline degasser and a variable wavelength detector (VWD), at a flow rate of 0.0.8 ml min −1 and detection wavelength at Scientific RepoRts | 5:18364 | DOI: 10.1038/srep18364 210 nm. In an attempt to achieve better separation, this analysis was carried out using a Carbo Sep Coregel-87H Guard Column, size 200 × 300 mm (Part No: ICE-99-9861) at a temperature of 65 °C.

16S rRNA gene pyrosequencing
Total DNA was extracted from the granules prior and at the end of the batches using the PowerSoil ® DNA Isolation Kit (MB 12888-MO BIO Laboratories, Inc.) with the following procedure; 2 g of sample was added to the PowerBead tubes. The samples were mechanically disrupted by bead beating at high speed for 60 seconds (Biopsec mini bead beater TM ) followed by vortexing for another 60 seconds. DNA concentration and quality were determined by nanodrop (ND-1000) before being used for 16s rRNA gene amplification through polymerase chain reaction (PCR). Amplification of the 16S rRNA gene region was verified using primer set of pyroL926F (5′ -aaactyaaakgaattgacgg-3′ ) and pyroL1392R (5′ -acgggcggtgtgtrc-3′ ) before submission to the Australian Centre for Ecogenomics (ACE, University of Queensland) for gene pyrotag sequencing on a 454 sequencing technology platform (Roche, USA) with the same primers.
Analysis of pyrosequencing data. Pyrosequencing data was grouped in operational taxonomic units (OTUs) at 97% similarity and aligned with the 16S rRNA identified sequences in the Greengenes database via the Quantitative Insights Into Microbial Ecology (QIIME) software package. These were further refined using the ACACIA software 48 .

Principal Component Analysis (PCA).
Normalised OTU values were Hellinger transformed to enhance subdominants 24 . A principal component analysis (PCA), based on eigenvector multivariate analysis, was done using the vegan functions RDA in RStudio (R-project for statistical computing-R Foundation) and a Biplot used to visualise data.
Phylogenetic analysis of pyrosequencing data. Cluster identity at species level was assessed using phylogenetic analysis. The 16S rRNA sequence from the 10 dominant clusters were aligned in MEGA 5.10 (molecular evolutionary genetics analysis -www.megasoftware.net) imported into the ARB software. These were first