Substrate characterization
The characteristics of the substrates rice straw and swine manure, are depicted in detail in Table 1 according to the physicochemical parameters evaluated. It was found that the percentage of TS in rice straw was higher than that of swine manure, making it a good carbon source for the digestion of lignocellulosic biomass [12]. The opposite occured with ammonium nitrogen where there was a higher contribution from swine manure, which is an essential nutrient for microorganisms. Therefore, the co-digestion of these substrates could effectively contribute to a better balance in the total C/N ratio, in addition to the contribution of micro and macronutrients, mainly K+, Na+, Ca2+, and Fe2+. Nitrogen deficiency could lead to failure of anaerobic digestion, but at high concentrations (> 3.5 g L− 1) can also lead to inhibition of methanogenic archaeal species [40]. The alkalinity ratio and pH values are in the range previously reported for swine manure [41, 42].
Table 1
Substrate characteristics before anaerobic digestion. Standard deviation in brackets (n = 3)
Parameter | Swine manure | Rice straw |
*TS (%) | 1.30 (0.06) | 89.95 (0.22) |
*VS (% TS) | 74.08 (0.80) | 80.55 (0.51) |
Ammonium Nitrogen (g·L− 1) | 2.07 (0.02) | 0.08 (0.03) |
Total alkalinity (g CaCO3·L− 1) | 3.1 (0.15) | n/d* |
Alkalinity ratio (α) | 0.53 (0.04) | n/d |
pH | 7.62 (0.06) | n/d |
C/N ratio | 1:19 | 48:1 |
*TS: total solids, VS: volatile solids, n/d = not determined |
The results shown in Table 1 are consistent with previous studies, that worked with the same substrates [4, 5, 40, 43]. The characteristics of the rice straw substrate, in terms of VS content, demonstrate the importance of the use of this biomass for energy purposes, despite its C/N ratio (48:1) being a little high. Furthermore, swine manure is considered a good substrate for anaerobic digestion, not only due to the variety of nutrients but also because of the high buffering capacity, which favors both bacterial and methanogenic trophic groups [44, 45]. Its combination with rice straw, creates a favorable balance for the C/N ratio, and allows the retention of the buffering capacity since the degradation of the straw can lead to the accumulation of volatile organic acids.
Evaluation of anaerobic digestion of rice straw
At start-up, the reactors showed similar methane yield values (without significant differences) during the first 10 days of the experimental stage (Fig. 1). Subsequently, a slight increase in methane yield values was observed in the bioaugmented reactor compared to the control, with increases on average of 50 units of normalized liters of methane per kg of volatile solid added. After 55 days, starting at an organic load of 1.24 gVSL− 1d− 1, this difference became noticeable.
A greater instability was observed in the control reactor decreasing the methane yield, probably due to the effects of the low degradability of the rice straw, influenced by the steric effect of non-degraded straw. In contrast, in the bioaugmented reactor, the methane production increased, obtaining values that doubled those obtained in the control reactor and increasing by 40% (v/v) the methane yield. The statistical analysis with the nonparametric test of multiple comparisons of related samples Mann Whitney (p < 0.05) corroborated these results, where significant differences existed between the control and both treatments, but not when the comparison was made between treatments.
The increase in methane yield obtained from the bioaugmentation strategy demonstrated that this could be a promising method to increase the efficiency of lignocellulosic biomass degradation. These results are in agreement with studies related to the impact of bioaugmentation as a bioremediation technique to increase methane potentials from cellulolytic wastes [18, 46, 47] and specifically with the use of microbial consortia to optimize bioethanol and methane production [23, 48–51].
During the degradation process, the pH ranged between 5.6 and 6.4 in both reactors. This may have been due to the addition of the lignocellulosic substrate rice straw which has a high C/N ratio and can accumulate high levels of volatile fatty acids (VFAs), resulting in a low pH [52]. However, Sun et al, [53] showed that pH might not be reliable enough as a monitoring indicator if the AD process is designed to treat high-strength organic waste, such as manure with much higher buffer capacity. With the increase in volumetric organic load, ammonium nitrogen levels also increased (Fig. 2) due to the addition of swine manure, showing a jump from weeks 9 and 10 (50–65 days of operation) in both reactors, being even higher in the bioaugmented reactor. Despite this, the highest methane yields were obtained in the bioaugmented reactor at this time, which is evidence that the high buffering capacity and the high concentrations of NH4+-N did not greatly affect the methanogenic archaea, although at higher levels (> 2.6 g L− 1) it can be inhibitory [40, 54]. Another indicator was the ratio between partial and total alkalinity (α), which should range between 0.5 and 0.9 [9]; in this case, values between 0.63 and 0.88 were obtained, which demonstrated the stability of the anaerobic process. The ratio of volatile fatty acids to total alkalinity (VFAs/TA) decreased during reactor start-up and remained in the optimal range for this parameter (0.1–0.35) [55] from weeks 3 to 10 (Fig. 2). After this time, VFAs/TA continued to decrease, which was closely related to alkaline pH values and reflected the acid-base balance of the system [56]. Some authors suggest that it may also be influenced by overloading conditions [9, 53, 57] or the effect of high ammonium concentrations that have their effect with accumulated VFAs, which would cause a decrease in TA. Alkalinity due to bicarbonates does not undergo any such change [56].
Molecular characterization of microbial consortium used
The taxonomic distribution up to the species level from the Bacterial Domain was detected in the microbial consortium through the metagenomic study. According to the consulted literature, there are few studies reporting similar metagenomic analyses in microbial consortia [58–61]. The Phyla Proteobacteria (71.8%), Firmicutes (21.55%), and Bacteroidetes (6.59%) predominated (Fig. 3). The genera Pseudomonas (17.74%), Lactobacillus (17.65%), Azotobacter (13.16%), and Arcobacter (8.15%), were identified with an abundance greater than 1%. Species of the genera Clostridium, Alcaligenes, Novispirillum, Bacteroides, Parabacteroides, Desulfotomaculum, Sedimentibacter, Dechlorospirillum, Aeromonas, Lysinibacillus and Acinetobacter were also identified, but in lower abundance. Of the total consortium community, 7% of species were unidentified, indicating that they were probably species that have not been previously detected in this type of sample by the Ion Torrent technique.
The metagenomic analysis of the microbial consortium is in agreement with the indices used to study the structure and organization of microbial community. A richness of 10 species was identified in the consortium related to the diversity according to Shannon's index (1.97), however, they had an equitable distribution in the community expressed in Pielou's equity index (0.85) since the value is close to 1.
Some of the genera identified in the microbial consortium (Lactobacillus, Lysinibacillus, Clostridium, Ruminococcus, Desulfotomaculum, Sedimentibacter, Pseudomonas, Acinetobacter, Arcobacter, Dechlorospirillum) match with those identified in a previous study of molecular characterization of efficient microorganisms (EM16) and microorganisms isolated from soil with agricultural residues [28], which were used for the development of this microbial consortium. In the previous study, the main metabolic pathways of the most abundant species were described. The presence of hydrolytic and cellulolytic species in the consortium demonstrates that it could positively influence the biodegradability of agricultural residues.
Effect of bioaugmentation on the microbial community dynamics in the reactors
Bacterial dynamics and diversity were analyzed in the reactors at different intervals (initial, 30 and 60 days). More than 95% of the community was identified to the species level through the Ion Torrent metagenomic technique. Of special interest for the analysis of the community in the reactors are those that are represented with a relative abundance higher than 1%, especially those predominant at 60 days, since they turned out to be the most stable in time and those that support the greatest methane formation.
In the control reactor, 31 operational taxonomic units (OTUs) with a relative abundance greater than 1% were found in the initial sample. At 30 and 60 days, 29 and 20 OTUs were obtained, respectively. This indicates that the number of dominant species decreased during reactor operation. The diversity determined using the Shannon index (H) was 3.30 at the initial time, and decreased as the days of operation passed to values of 2.76 at 60 days. The Pielou's index (H/Hmax) in this reactor showed values of 0.96 from the initial moment to 0.92 at 60 days, which elucidated the evenness of the community, which was close to 1.
According to Jaccard's similarity index (J'), 46% of the species found at the initial time were also detected at 30 days, but of these, only 31% were detected at 60 days in the control reactor. A comparison of the samples at 30 and 60 days showed that 63% of the species were similar. Of the microbial groups identified with the highest abundance in this reactor (> 1%), 12 were detected at all three sampling times. However, at 30 days, 10 families were detected that were not detected at the initial time, and at 60 days only 1 family was identified that was not detected either at the initial time or at 30 days. These microbial dynamics suggest adaptation and specialization of the anaerobic sludge as the organic load of the reactor increased over time.
In the bioaugmented reactor, 27 OTUs with a relative abundance above 1% were found in the initial sample. At 30 and 60 days, 30 and 21 OTUs were obtained, respectively. In this reactor at 30 days, the diversity of species increased slightly with respect to the initial sample, but as in the control reactor, the diversity decreased during reactor operating days. According to Shannon's index, this reactor showed values of 3.14 at the initial moment to 2.80 at 60 days, and the Pielou evenness index value of 0.95 from the initial moment to 0.94 at 60 days, which was also close to 1.
According to Jaccard's similarity index (J'), 63% of the species found at the initial time were also detected at 30 days, but of these only 27% were detected at 60 days in the bioaugmented reactor. Comparing the samples at 30 and 60 days, it can be said that 32% of the species were similar (lower similarity index than in the control reactor), which may have been because in this reactor the microbial community was exposed to more significant variability due to the addition of the consortium. Of the microbial groups identified with the highest abundance in this reactor (> 1%), 10 were detected at all three sampling times. However, at 30 days, 7 were detected that were not detected at the initial time and at 60 days, and 4 were detected that were not with a relative abundance greater than 1% at the initial time or at 30 days. It is possible that certain species found initially in the inoculum and substrates may not have been favored by increases in organic load and others that were better adapted to this substrate may have benefited. As in the control reactor, these microbial dynamics suggest a better specialization of the community.
Different microbial populations existed in each of the stages of anaerobic digestion, each with specific metabolic pathways. As observed in the results, there were also differences in microbial diversity and dynamics between different reactor operation times. Therefore, it can be stated that there was a close correlation between the operating time and the microbial diversity of the reactors. In the bioaugmented reactor, with the specialization of the community, the diversity of species decreased, but the abundance of those capable of surviving and degrading these lignocellulosic substrates increased, which had a positive effect on methane production, corroborating the higher yield obtained in the bioaugmented reactor to the control. Similar results have been obtained by other authors [62–65].
Analysis of the bacterial community in the reactors
In the control reactor at the initial time, the following families were detected in the highest abundance: Clostridiaceae (genus Clostridium, 5.72%, Clostridium butyricum, 1.98%, Clostridium sp., 1.08%), Chlorobiaceae (genus Chlorobium, 5.49%), Syntrophaceae (4.13%), Ignavibacteriaceae (3.97%), Synergistaceae (3.73%), Dehalococcoidaceae (3.14%), Prolixibacteraceae (3.12%), Erysipelotrichaceae (Turicibacter sanguinis, 2.87%), Planctomycetaceae (2.67%), Syntrophorhabdaceae (genus Syntrophorhabdus, 2.53%), Marinilabiliaceae (2.51%), Geobacteraceae (2.43%), Cytophagaceae (2.42%), Eubacteriaceae (2.25%), and Anaerolineaceae (2.04%) (Fig. 4).
In the sample taken at 30 days the families Prolixibacteraceae (4.6%), Porphyromonadaceae (Parabacteroides chartae, 4.39%), Synergistaceae (3.61%), Clostridiaceae (genus Clostridium 3, 25% and Clostridium butyricum, 1.43%), Cytophagaceae (3.39%), Dehalococcoidaceae (2.61%), Syntrophaceae (2.55%), Marinilabiliaceae (2.44%), Flavobacteriaceae (2.34%) and Eubacteriaceae (2.08%) were found. In addition, others that were not detected in the initial sample were found, such as Porphyromonadaceae (2.78%), Ruminococcaceae (2.57%), Moraxellaceae (genus Acinetobacter, 2.37%), Oxalobacteraceae (genus Undibacterium, 2.09), Bacteroidaceae (Bacteroides sp., 1.06% and Bacteroides graminisolvens, 1.48%), Paenibacillaceae (1.06%) and Prevotellaceae (1.02%) (Fig. 4).
At 60 days in the control reactor the families Clostridiaceae (8.7%), Ruminococcaceae (6.62%) Prolixibacteraceae (5.84%), Moraxellaceae (genus Acinetobacter, 5.73%), Synergistaceae (3.99%), Cytophagaceae (3.56%) Prevotellaceae (3, 29%), Porphyromonadaceae (3.28%), Marinilabiliaceae (2.52%), Porphyromonadaceae (Parabacteroides chartae, 2.4%), Flavobacteriaceae (2.06%), Dehalococcoidaceae (1.59%) and the previously undetected family Nocardioidaceae (1.31%) were found (Fig. 4).
In each sample evaluated from the control reactor there was a percentage of species that were not detected by the Ion Torrent metagenomic technique, which increased with time (1%, 3% and 5%, at the beginning, 30 and 60 days respectively), indicating that they are probably sequences not previously reported in the NCBI database, which inhabit this type of ecosystems.
In the bioaugmented reactor, at start up, the Clostridiaceae family was detected in highest abundance (17.25%). Other families in high proportion included; Erysipelotrichaceae 11.41% (genus Turicibacter, 5.34%, and Turicibacter sanguinis 6.07%); Ignavibacteriaceae (7.71%); Peptostreptococcaceae (6.27%); Syntrophorhabdaceae (6.21%) (genus Syntrophorhabdus, 3.41%); Syntrophaceae (4.55%); Synergistaceae (3.93%); Marinilabiliaceae (3.86%); Cytophagaceae (3.56%); Eubacteriaceae (3.17%); Geobacteraceae (3.06%), Dehalococcoidaceae (2.35%) and Prolixibacteraceae (2.32%) (Fig. 5). Among these families identified at the initial time it is interesting that some of them were also identified in the microbial consortium, such as: Porphyromonadaceae, Prevotellaceae, Cytophagaceae, Flavobacteriaceae, Sphingobacteriaceae (phylum Bacteroidetes), Bacillaceae, Clostridiaceae, Eubacteriaceae, Lachnospiraceae, Peptococcaceae, Peptostreptococcaceae (phylum Firmicutes), and Alcaligenaceae (phylum Proteobacteria).
At 30 days in the bioaugmented bioreactor, the Clostridiaceae families predominated with 16.04% (genus Clostridium, 5.45%, Clostridium butyricum, 3.82%, Clostridium disporicum, 1.09%, Clostridium sp., 1.49%), Erysipelotrichaceae 9.89% (genus Turicibacter, 4.19%, and Turicibacter sanguinis, 4.75%), Prolixibacteraceae (8.88%), Syntrophorhabdaceae 5.3% (genus Syntrophorhabdus, 2.89%), Porphyromonadaceae (4.91%), Bacteroidaceae 4, 13% (Bacteroides graminisolvens, 1.01%), Peptostreptococcaceae (3.87%), Cytophagaceae (3.25%), Ignavibacteriaceae (2.98%), Ruminocococcaceae (2.83%), Marinilabiliaceae (2.59%), Synergistaceae (2.57%), Flavobacteriaceae (2.33%) and Eubacteriaceae (2.23%) (Fig. 5). The bacterial community at 30 days, had change in composition from the sample taken at the beginning of the experiment. Species that probably came from the microbial consortium were favored, increasing their abundance, such as the families Bacteroidaceae, Porphyromonadaceae, Prevotellaceae, Flavobacteriaceae, Bacillaceae, Paenibacillaceae, Ruminococcaceae, Rhodospirillaceae, Campylobacteraceae, Aeromonadaceae, Moraxellaceae, and Pseudomonadaceae. However, others decreased in abundance such as Clostridiaceae, Eubacteriaceae, and Peptostreptococcaceae. This decrease was not very significant however, because at 30 days these three families represented 22.14% of the total community.
At 60 days the families Clostridiaceae (9.76%), Prolixibacteraceae (8.58%), Prevotellaceae (7.21%), Spirochaetaceae (6.51%), Bacteroidaceae (5.52%), Porphyromonadaceae (5.05%), Ruminococcaceae (5%), Leptospiraceae (4, 33%), Cytophagaceae (3.95%), Syntrophaceae (3.89%), Marinilabiliaceae (2.97%), Erysipelotrichaceae (2.48%) Pseudanabaenaceae (2.37%), Pseudomonadaceae 2.15%, and Flavobacteriaceae (2.11%) were more abundant in the bioaugmented bioreactor (Fig. 5). Many of these families were also identified at baseline, at 30 days and in the microbial consortium. Similar behavior was observed in relation to the sample taken at 30 days, where many of the families from the consortium increased their abundance.
If the two reactors are compared, it can be noted that some of the families that are in the bioaugmented reactor were not present in the control reactor because they came from the addition of the consortium, among these were Porphyromonadaceae, Prevotellaceae, Sphingobacteriaceae, Bacillaceae, Lactobacillaceae, Eubacteriaceae, Lachnospiraceae, Peptococcaceae, Ruminocococcaceae, Peptostreptococcaceae, Acetobacteraceae, Alcaligenaceae, Campylobacteraceae, Aeromonadaceae, Moraxellaceae, and Pseudomonadaceae. In this reactor as in the control, in the 3 sampling moments, some species were not identified by this technique. In this case unknown species were less than 1% in the initial sample, 2% at 30 days, and 3% at 60 days.
In both reactors, bacterial species of the phylum Firmicutes predominated, which is similar to the observations of other authors [9, 66, 67]. Cellulose degrading species were dominant, mainly from the Clostridiaceae family. Other families detected in both reactors with a relative abundance greater than 1% predominated up to 60 days, which could be inferred to be part of the basic structure of the community.
As mentioned above, some of the families from the bioaugmented reactor were also part of the consortium community, so the addition of the consortium to the reactor contributed to strengthening its microbial population. The fact that after 60 days, species from the consortium were maintained in this reactor suggests that they easily adapted to the environmental and nutritional conditions inside the reactor and compete with the species that came from the inoculum or substrates. Therefore, it is likely that they increased the capacity to degrade the lignocellulosic material found in the rice straw, which in turn supported the increase in methane production.
Methanogenic archaea analysis
Real-time amplification of genes coding for the four groups of methanogenic archaea of interest revealed that all groups (acetoclastic and hydrogenotrophic methanogens) were detected in both the control and bioaugmented reactors. However, species of the order Methanosarcinales and the family Methanosaetacea predominated in both reactors (Fig. 6). In the control reactor, no significant differences were found among the three time samples analyzed, where similar profiles of methanogenic species were shown. In the bioaugmented reactor, the total number of methanogenic archaea decreased slightly from the initial time until 60 days.
In the reactors, the methanosarcins and methanosaetas, are favored [68]. Members of this group are frequently isolated from anaerobic reactors under mesophilic conditions, which suggests that the acetoclastic metabolic pathway is favored for methane production [69, 70]. In addition, the presence of species belonging to the orders Methanobacteriales and Methanomicrobiales, although in lower abundance, suggests that the anaerobic process was balanced, and in turn allowed the abundance of species of the order Methanosarcinales responsible for the increase in methane yield [68, 71].
This fact contrasts with Kim et al., [72], who detected a predominance of the orders Methanobacteriales and Methanomicrobiales in batch reactors with swine manure. However, these results are in agreement with those obtained by Jiménez-Hernández et al., [73] and Padmasiri et al., [70], in which these groups of archaea predominate under mesophilic conditions. In general, Metanosarcina ssp. has a higher maximum growth rate but a lower affinity to acetate than Metanosaeta ssp. [71, 74]. Therefore, it is to be expected that when acetate is limited, the Methanosaetaceae family predominates. The coexistence of different groups of methanogens in the same ecosystem has been reported, as they are key players in methane production and preserve the functional stability of the ecosystem [75–77]. However, substrate availability and/or concentration and environmental changes seem to be crucial factors regulating the change of methanogenic communities [76, 78].
The microbial diversity and dynamics of both bacteria and archaea can support the positive effect that the addition of the enriched microbial consortium had on the anaerobic digestion process of rice straw residues. This allowed an increase in the utilization of the methane potential contained in this substrate even when the organic load was increased. This is evidence of the importance of the degradative capacity when microorganisms (bioaugmentation) are applied to anaerobic reactors either because the metabolic pathways of degradation of the species present or those that are added are favored.