Group-specific primer and probe sets to detect methanogenic communities using quantitative real-time polymerase chain reaction

Evaluation of system performance and microbial communities of a temperature-phased anaerobic digestion system treating dairy manure: thermophilic digester operated at acidic pH

A temperature-phased anaerobic digestion system with the thermophilic digester acidified by acidogenesis products (referred to as AT-TPAD) was evaluated to treat high-strength dairy cattle manure at a 15-day retention time. Three temperatures (50, 55, and 60°C) were tested on the thermophilic digester, and 50°C was found to be the optimal temperature for overall performance of the AT-TPAD system, achieving 31% VS removal and 0.22 L methane/g VS fed. The mesophilic digester contributed significantly more to the overall system performance than the thermophilic digester. The thermophilic and the mesophilic digesters had different microbial communities under all conditions, and both microbial communities exhibited dynamic changes in response to different conditions. Certain microbial groups were found significantly correlated with the system performance. Methanosarcina was the most important methanogen genus of the AT-TPAD system and its population abundance was inversely correlated with high concentrations of volatile fatty acids (VFA).

Use of real-time QPCR in biokinetics and modeling of two different ammonia-oxidizing bacteria growing simultaneously

A real-time quantitative polymerase chain reaction (QPCR) was used to evaluate biokinetic coefficients of Nitrosomonas nitrosa and N. cryotolerans clusters growing simultaneously in a batch mode of ammonia oxidation. The mathematical models based on Monod equation were employed to describe the competitive relationship between these clusters and were fitted to experimental data to obtain biokinetic values. The maximum growth rates (μ(m)), half-saturation coefficients (K(S)), microbial yields (Y) and decay coefficients (k(d)) of N. nitrosa and N. cryotolerans were 1.77 and 1.21 day(-1), 23.25 and 23.06 mg N·L(-1), 16 × 10(8) and 1 × 10(8) copies·mg N(-1), 0.26 and 0.20 day(-1), respectively. The estimated coefficients were applied for modeling continuous operations at various hydraulic retention times (HRTs) with an influent ammonia concentration of 300 mg N·L(-1). Modeling results revealed that ammonia oxidation efficiencies were achieved 55-98 % at 0.8-10 days HRTs and that the system was predicted to be washed out at HRT of 0.7 days. Overall, use of QPCR for estimating biokinetic coefficients of the two AOB cluster growing simultaneously by use of ammonia were successful. This idea may open a new direction towards biokinetics of ammonia oxidation in which respirometry tests are usually employed.

Organic loading rate shock impact on operation and microbial communities in different anaerobic fixed-bed reactors

For the fixed-bed reactors in this experiment, during 40 days of stable operation and under different organic loading shocks, biogas production remained stable at 21 L, effluent pH remained between 6.8 and 7.5, and chemical oxygen demand (COD) removal efficiency and the biogas methane content were greater than 80% and 75%, respectively. The community was analyzed using denaturing gradient gel electrophoresis (DGGE), 16S rRNA gene clone library screening, and quantitative PCR. Findings revealed that bacteria and methanogenic archaea were typically dominant in the adhering sludge. Methanomicrobiales was identified in carbon fiber carriers, they were breeding slowly, and attached easily. The 16S rRNA gene concentration of methanogenic archaea was higher in the adhering sludge than in the deposited sludge. Our results indicated that the colonization of the microorganism played a very important role in the carbon fiber carriers, as well as in the improvement of sludge activity and the shock resistance of the reactor.

Anaerobic-aerobic sequencing bioreactors improve energy efficiency for treatment of personal care product industry wastes

Personal care product (PCP) industry liquid wastes contain shampoo residues, which are usually treated by aerobic activated sludge (AS). Unfortunately, AS is expensive for PCP wastes because of high aeration and energy demands, whereas potentially energy-positive anaerobic designs cannot meet effluent targets. Therefore, combined anaerobic-aerobic systems may be the best solution. Seven treatment systems were assessed in terms of energy and treatment performance for shampoo wastes, including one aerobic, three anaerobic (HUASB, AHR and AnCSTR) and three anaerobic-aerobic reactor designs. COD removals were highest in the HUASB-aerobic (87.9 ± 0.4%) and AHR-aerobic (86.8±0.5%) systems, which used 69.2% and 62.5% less energy than aerobic AS. However, actual methane production rates were low relative to theoretical in the UASB and AHR units (∼10% methane/COD removed) compared with the AnCSTR unit (∼70%). Anaerobic-aerobic sequence reactors show promise for treating shampoo wastes, but optimal designs depend upon whether methane production or COD removal is most important to operations.

High-rate iron-rich activated sludge as stabilizing agent for the anaerobic digestion of kitchen waste

Anaerobic digestion is a key technology in the bio-based economy and can be applied to convert a wide range of organic substrates into CH4 and CO2. Kitchen waste is a valuable substrate for anaerobic digestion, since it is an abundant source of organic matter. Yet, digestion of single kitchen waste often results in process failure. High-rate activated sludge or A-sludge is produced during the highly loaded first stage of the two-phase 'Adsorptions-Belebungsverfahren' or A/B activated sludge system for municipal wastewater treatment. In this specific case, the A-sludge was amended with FeSO4 to enhance phosphorous removal and coagulation during the water treatment step. This study therefore evaluated whether this Fe-rich A-sludge could be used to obtain stable methanation and higher methane production values during co-digestion with kitchen waste. It was revealed that Fe-rich A-sludge can be a suitable co-substrate for kitchen waste; i.e. methane production rate values of 1.15 ± 0.22 and 1.12 ± 0.28 L L(-1) d(-1) were obtained during mesophilic and thermophilic co-digestion respectively of a feed-mixture consisting of 15% KW and 85% A-sludge. The thermophilic process led to higher residual VFA concentrations, up to 2070 mg COD L(-1), and can therefore be considered less stable. Addition of micro- and macronutrients provided a more stable digestion of single kitchen waste, i.e. a methane production of 0.45 L L(-1) d(-1) was obtained in the micronutrient treatment compared to 0.30 L L(-1) d(-1) in the control treatment on day 61. Yet, methane production during single kitchen waste digestion still decreased toward the end of the experiment, despite the addition of micronutrients. Methane production rates were clearly influenced by the total numbers of archaea in the different reactors. This study showed that Fe-rich A-sludge and kitchen waste are suitable for co-digestion.

DNA-SIP reveals that Syntrophaceae play an important role in methanogenic hexadecane degradation

The methanogenic degradation of linear alkanes is a common process in oil-impacted environments. However, little is known about the key players involved in this process. Here, the hexadecane-degrading organisms in a methanogenic, hexadecane-degrading consortium designated M82 obtained from Shengli oilfield and maintained at 35°C for over 4 years, were identified by DNA-stable isotope probing with UL-¹³C-hexadecane, followed by density-resolved terminal restriction fragment length polymorphism (T-RFLP) analysis, cloning and phylogenetic analysis of 16S rRNA gene fragments. Compared to the fractions of the ¹²C treatment, the relative abundance of two phylotypes significantly increased in the heavy fractions of the ¹³C-hexadecane incubated microcosm. One belongs to a uncultured member of the bacterial family Syntrophaceae, which show 95–97% rRNA sequence identity with Smithella propionica, and the other is affiliated with Methanoculleus receptaculi (>99% sequence identity). The results of the present study prove the significant role of uncultured Syntrophaceae in degradation of hexadecane, probably through syntrophic interactions with hydrogenotrophic methanogens.

Monitoring methanogenic population dynamics in a full-scale anaerobic digester to facilitate operational management

Microbial populations in a full-scale anaerobic digester fed on food waste were monitored over an 18-month period using qPCR. The digester exhibited a highly dynamic environment in which methanogenic populations changed constantly in response to availability of substrates and inhibitors. The methanogenic population in the digester was dominated by Methanosaetaceae, suggesting that aceticlastic methanogenesis was the main route for the production of methane. Sudden losses (69%) in Methanosaetaceae were followed by a build-up of VFAs which were subsequently consumed when populations recovered. A build up of ammonium inhibited Methanosaetaceae and resulted in shifts from acetate to hydrogen utilization. Addition of trace elements and alkalinity when propionate levels were high stimulated microbial growth. Routine monitoring of microbial populations and VFAs provided valuable insights into the complex processes occurring within the digester and could be used to predict digester stability and facilitate digester optimization.

454 pyrosequencing analyses of bacterial and archaeal richness in 21 full-scale biogas digesters

The microbial community of 21 full-scale biogas reactors was examined using 454 pyrosequencing of 16S rRNA gene sequences. These reactors included seven (six mesophilic and one thermophilic) digesting sewage sludge (SS) and 14 (ten mesophilic and four thermophilic) codigesting (CD) various combinations of wastes from slaughterhouses, restaurants, households, etc. The pyrosequencing generated more than 160,000 sequences representing 11 phyla, 23 classes, and 95 genera of Bacteria and Archaea. The bacterial community was always both more abundant and more diverse than the archaeal community. At the phylum level, the foremost populations in the SS reactors included Actinobacteria, Proteobacteria, Chloroflexi, Spirochetes, and Euryarchaeota, while Firmicutes was the most prevalent in the CD reactors. The main bacterial class in all reactors was Clostridia. Acetoclastic methanogens were detected in the SS, but not in the CD reactors. Their absence suggests that methane formation from acetate takes place mainly via syntrophic acetate oxidation in the CD reactors. A principal component analysis of the communities at genus level revealed three clusters: SS reactors, mesophilic CD reactors (including one thermophilic CD and one SS), and thermophilic CD reactors. Thus, the microbial composition was mainly governed by the substrate differences and the process temperature.

Molecular Genetic Diversity and Quantitation of Methanogen in Ruminal Fluid of Buffalo (Bubalus bubalis) Fed Ration (Wheat Straw and Concentrate Mixture Diet)

High roughage diet causes more methane emissions; however, the total methanogen abundance is not influenced by roughage proportion. Technologies to reduce methane emissions are lacking, and development of inhibitors and vaccines that mitigate rumen-derived methane by targeting methanogens relies on present knowledge of the methanogens. In this work, we have investigated molecular diversity of rumen methanogens of Surti buffalo. DNA from rumen fluid was extracted, and 16S rRNA encoding genes were amplified using methanogen specific primer to generate 16S rDNA clone libraries. Seventy-six clones were randomly selected and analysed by RFLP resulting in 21 operational taxonomic units (OTUs). BLAST analysis with available sequences in database revealed sequences of 13 OTUs (55 clones) showing similarity with Methanomicrobium sp, 3 OTUs (15 clones) with Methanobrevibacter sp. The remaining 5 OTUs (6 clones) belonged to uncultured archaea. The phylogenetic analysis indicated that methanogenic communities found in the library were clustered in the order of Methanomicrobiales (18 OTUs) and Methanobacteriales (3 OTUs). The population of Methanomicrobiales, Methanobacteriales, and Methanococcales were also observed, accounting for 1.94%, 0.72%, and 0.47% of total archaea, respectively.

UASB performance and electron competition between methane-producing archaea and sulfate-reducing bacteria in treating sulfate-rich wastewater containing ethanol and acetate

To find an appropriate method for sulfate-rich wastewater containing ethanol and acetate with COD/sulfate ratio of 1, a UASB reactor was operated for more than 180 days. The influences of HRT (hydraulic retention time) and OLR (organic loading rate) on organics and sulfate removal, gas production, and electrons utilization were investigated. The sludge activity and microorganism composition were also determined. The results indicated that this system removed more than 80% of COD and 30% of sulfate with HRT above 6h and OLR below 12.3 gCOD/L d. Further HRT decrease caused volatile fatty acids accumulation and performance deterioration. Except at HRT of 2h, COD and electron flow were mostly utilized by methane-producing archaea (MPA), and methane yield remained in the range of 0.18-0.24 LCH4/gCOD. Methane was mainly generated by Methanosaeta concilii GP6 with acetate as substrate, whereas sulfate was mainly reduced by incomplete-oxidizing Desulfovibrio species with ethanol as substrate.

The human gut chip "HuGChip", an explorative phylogenetic microarray for determining gut microbiome diversity at family level

Evaluating the composition of the human gut microbiota greatly facilitates studies on its role in human pathophysiology, and is heavily reliant on culture-independent molecular methods. A microarray designated the Human Gut Chip (HuGChip) was developed to analyze and compare human gut microbiota samples. The PhylArray software was used to design specific and sensitive probes. The DNA chip was composed of 4,441 probes (2,442 specific and 1,919 explorative probes) targeting 66 bacterial families. A mock community composed of 16S rRNA gene sequences from intestinal species was used to define the threshold criteria to be used to analyze complex samples. This was then experimentally verified with three human faecal samples and results were compared (i) with pyrosequencing of the V4 hypervariable region of the 16S rRNA gene, (ii) metagenomic data, and (iii) qPCR analysis of three phyla. When compared at both the phylum and the family level, high Pearson's correlation coefficients were obtained between data from all methods. The HuGChip development and validation showed that it is not only able to assess the known human gut microbiota but could also detect unknown species with the explorative probes to reveal the large number of bacterial sequences not yet described in the human gut microbiota, overcoming the main inconvenience encountered when developing microarrays.

Elucidating the microbial resuscitation cascade in biological soil crusts following a simulated rain event

Biological soil crusts (biocrusts) are photosynthetic mats formed through an association of prokaryotic and eukaryotic microorganisms with soil particles. Biocrusts are found in virtually any terrestrial ecosystem where vascular plant coverage is abiotically limited, with drylands comprising the primary habitat for them. We studied the dynamics of the active bacterial community in two biocrusts from an arid and a hyperarid region in the Negev Desert, Israel, under light-oxic and dark-anoxic incubation conditions after simulated rainfall. We used H2(18)O for hydrating the crusts and analysed the bacterial community in the upper and lower parts of the biocrust using an RNA-stable isotope probing approach coupled with 454-pyrosequencing. In both biocrusts, two distinct bacterial communities developed under each incubation condition. The active anaerobic communities were initially dominated by members of the order Bacillales which were later replaced by Clostridiales. The aerobic communities on the other hand were dominated by Sphingobacteriales and several Alphaproteobacteria (Rhizobiales, Rhodobacterales, Rhodospirillales and Rubrobacteriales). Actinomycetales were the dominant bacterial order in the dry crusts but quickly collapsed and accounted for < 1% of the community by the end of the incubation. Our study shows that biocrusts host a diverse community whose members display complex interactions as they resuscitate from dormancy.

Quantitatively different, yet qualitatively alike: a meta-analysis of the mouse core gut microbiome with a view towards the human gut microbiome

Background

A number of human diseases such as obesity and diabetes are associated with changes or imbalances in the gut microbiota (GM). Laboratory mice are commonly used as experimental models for such disorders. The introduction and dynamic development of next generation sequencing techniques have enabled detailed mapping of the GM of both humans and animal models. Nevertheless there is still a significant knowledge gap regarding the human and mouse common GM core and thus the applicability of the latter as an animal model. The aim of the present study was to identify inter- and intra-individual differences and similarities between the GM composition of particular mouse strains and humans.

Methodology/Principal Findings

A total of 1509428 high quality tag-encoded partial 16S rRNA gene sequences determined using 454/FLX Titanium (Roche) pyro-sequencing reflecting the GM composition of 32 human samples from 16 individuals and 88 mouse samples from three laboratory mouse strains commonly used in diabetes research were analyzed using Principal Coordinate Analysis (PCoA), nonparametric multivariate analysis of similarity (ANOSIM) and alpha diversity measures. A reliable cutoff threshold for low abundant taxa estimated on the basis of the present study is recommended for similar trials.

Conclusions/Significance

Distinctive quantitative differences in the relative abundance of most taxonomic groups between the examined categories were found. All investigated mouse strains clustered separately, but with a range of shared features when compared to the human GM. However, both mouse fecal, caecal and human fecal samples shared to a large extent not only representatives of the same phyla, but also a substantial fraction of common genera, where the number of shared genera increased with sequencing depth. In conclusion, the GM of mice and humans is quantitatively different (in terms of abundance of specific phyla and species) but share a large qualitatively similar core.

Microbial community structure and dynamics in two-stage vs single-stage thermophilic anaerobic digestion of mixed swine slurry and market bio-waste

The microbial community of a thermophilic two-stage process was monitored during two-months operation and compared to a conventional single-stage process. Qualitative and quantitative microbial dynamics were analysed by Denaturing Gradient Gel Electrophoresis (DGGE) and real-time PCR techniques, respectively. The bacterial community was dominated by heat-shock resistant, spore-forming clostridia in the two-stage process, whereas a more diverse and dynamic community (Firmicutes, Bacteroidetes, Synergistes) was observed in the single-stage process. A significant evolution of bacterial community occurred over time in the acidogenic phase of the two-phase process with the selection of few dominant species associated to stable hydrogen production. The archaeal community, dominated by the acetoclastic Methanosarcinales in both methanogen reactors, showed a significant diversity change in the single-stage process after a period of adaptation to the feeding conditions, compared to a constant stability in the methanogenic reactor of the two-stage process. The more diverse and dynamic bacterial and archaeal community of single-stage process compared to the two-stage process accounted for the best degradation activity, and consequently the best performance, in this reactor. The microbiological perspective proved a useful tool for a better understanding and comparison of anaerobic digestion processes.

Microbial community signature of high-solid content methanogenic ecosystems

In this study, changes in bacterial and archaeal communities involved in anaerobic digestion processes operated with high solid contents were investigated. Batch tests were performed within a range of total solids (TS) of 10-35%. Between 10% and 25% TS, high methanogenic activity was observed and no overall specific structure of active bacterial communities was found. At 30% and 35%, methanogenesis was inhibited as a consequence of volatile fatty acids accumulation. Here, a specific bacterial signature was observed with three main dominant bacteria related to Clostridium sp., known for their ability to grow at low pH. Additionally, archaeal community was gradually impacted by TS content. Three archaeal community structures were observed with a gradual shift from Methanobacterium sp. to Methanosarcina sp., according to the TS content. Overall, several species were identified as biomarkers of methanogenesis inhibition, since bacterial and archaeal communities were highly specific at high TS contents.

Improving PCR efficiency for accurate quantification of 16S rRNA genes

Quantitative real-time PCR is a valuable tool for microbial ecologists. To obtain accurate absolute quantification it is essential that PCR efficiency for pure standards is close to amplification efficiency for test samples. Counter to normal expectation that PCR efficiency might be lower in environmental DNA, due to the presence of PCR inhibitors, we report the counterintuitive observation that PCR efficiency of pure standards can be lower than for environmental DNA. This can lead to overestimation of gene abundances if not corrected. SYBR green-based qPCR assays of 16S rRNA genes targeting Bacteria, Syntrophus and Smithella spp., Marinobacter spp., Methanomicrobiales, Methanosarcinaceae, and Methanosaetaceae in samples from methanogenic crude oil biodegradation enrichments were tested. In five out of the six assays, PCR efficiency was lower with pure standards than with environmental DNA samples. We developed a solution to this problem based on amending pure clone standards with a background of non-target environmental 16S rRNA genes which significantly improved PCR efficiency of standards in the qPCR assays that exhibited this phenomenon. Overall this method of qPCR standard preparation achieved a more reliable and robust quantification of 16S rRNA genes. We believe this may be a potentially common issue in microbial ecology that often goes unreported, as intuitively one would not expect standards to have poorer PCR efficiency than samples.

Archaeal and anaerobic methane oxidizer communities in the Sonora Margin cold seeps, Guaymas Basin (Gulf of California)

Cold seeps, located along the Sonora Margin transform fault in the Guaymas Basin, were extensively explored during the 'BIG' cruise in June 2010. They present a seafloor mosaic pattern consisting of different faunal assemblages and microbial mats. To investigate this mostly unknown cold and hydrocarbon-rich environment, geochemical and microbiological surveys of the sediments underlying two microbial mats and a surrounding macrofaunal habitat were analyzed in detail. The geochemical measurements suggest biogenic methane production and local advective sulfate-rich fluxes in the sediments. The distributions of archaeal communities, particularly those involved in the methane cycle, were investigated at different depths (surface to 18 cm below the sea floor (cmbsf)) using complementary molecular approaches, such as Automated method of Ribosomal Intergenic Spacer Analysis (ARISA), 16S rRNA libraries, fluorescence in situ hybridization and quantitative polymerase chain reaction with new specific primer sets targeting methanogenic and anaerobic methanotrophic lineages. Molecular results indicate that metabolically active archaeal communities were dominated by known clades of anaerobic methane oxidizers (archaeal anaerobic methanotroph (ANME)-1, -2 and -3), including a novel 'ANME-2c Sonora' lineage. ANME-2c were found to be dominant, metabolically active and physically associated with syntrophic Bacteria in sulfate-rich shallow sediment layers. In contrast, ANME-1 were more prevalent in the deepest sediment samples and presented a versatile behavior in terms of syntrophic association, depending on the sulfate concentration. ANME-3 were concentrated in small aggregates without bacterial partners in a restricted sediment horizon below the first centimetres. These niche specificities and syntrophic behaviors, depending on biological surface assemblages and environmental availability of electron donors, acceptors and carbon substrates, suggest that ANME could support alternative metabolic pathways than syntrophic anaerobic oxidation of methane.

Activated zeolite--suitable carriers for microorganisms in anaerobic digestion processes?

Plant cell wall structures represent a barrier in the biodegradation process to produce biogas for combustion and energy production. Consequently, approaches concerning a more efficient de-polymerisation of cellulose and hemicellulose to monomeric sugars are required. Here, we show that natural activated zeolites (i.e. trace metal activated zeolites) represent eminently suitable mineral microhabitats and potential carriers for immobilisation of microorganisms responsible for anaerobic hydrolysis of biopolymers stabilising related bacterial and methanogenic communities. A strategy for comprehensive analysis of immobilised anaerobic populations was developed that includes the visualisation of biofilm formation via scanning electron microscopy and confocal laser scanning microscopy, community and fingerprint analysis as well as enzyme activity and identification analyses. Using SDS polyacrylamide gel electrophoresis, hydrolytical active protein bands were traced by congo red staining. Liquid chromatography/mass spectroscopy revealed cellulolytical endo- and exoglucanase (exocellobiohydrolase) as well as hemicellulolytical xylanase/mannase after proteolytic digestion. Relations to hydrolytic/fermentative zeolite colonisers were obtained by using single-strand conformation polymorphism analysis (SSCP) based on amplification of bacterial and archaeal 16S rRNA fragments. Thereby, dominant colonisers were affiliated to the genera Clostridium, Pseudomonas and Methanoculleus. The specific immobilisation on natural zeolites with functional microbes already colonising naturally during the fermentation offers a strategy to systematically supply the biogas formation process responsive to population dynamics and process requirements.

Antigen dose and humoral immune response correspond with protection for inactivated infectious pancreatic necrosis virus vaccines in Atlantic salmon (Salmo salar L)

An enduring challenge in the vaccinology of infectious pancreatic necrosis virus (IPNV) is the lack of correlation between neutralizing antibodies and protection against mortality. To better understand the immunological basis of vaccine protection, an efficacy trial including Atlantic salmon (Salmo salar L.) vaccinated with a high antigen (HiAg) or low antigen (LoAg) dose vaccine was carried out in a cohabitation challenge model using the highly virulent Norwegian Sp strain NVI015. To pinpoint the immunological basis of vaccine protection, pathogenic mechanisms of IPNV were unraveled in control fish while obtaining feedback on mechanisms of protection in the vaccinated fish. During the incubation period, infection rates were highest in control fish, followed by the LoAg group with the lowest infections being in the HiAg group. Although both the liver and pancreas are target organs prone to tissue damage, infection in the liver was delayed until acute infection in most fish. A correlate of pathology determined as the cutoff threshold of viral copy numbers linked to tissue damage in target organs was estimated at ≥ 10^7.0, which corresponded with an increase in mortality. The kinetics of IFNα and Mx expression suggests that these genes can be used as biomarkers of IPNV infection progression. Mechanisms of vaccine protection involved reducing infection rates, preventing infection of the liver and reducing virus replication in target organs to levels below the correlate of pathology. Overall, the study shows that antigen dose corresponds with vaccine efficacy and that antibody levels can be used as a signature of protective immunity against pathological disease and mortality.

Microbial ecosystem and fermentation traits in the caecum of growing rabbits given diets varying in neutral detergent soluble and insoluble fibre levels

The effect of the level of neutral detergent fibre (NDF: 0.35, LI and 0.42, HI) and neutral detergent soluble fibre (NDSF: 0.14, LS and 0.17, HS) in the caecal ecosystem was studied in 24 weaned (28 days of age) rabbits, weighing 630 ± 80.2 g in a 2 × 2 factorial design. After 22 days, rabbits were slaughtered and their caecal contents sampled. The caecal pH (on average 6.2) and molar volatile fatty acids (VFA) proportions were not affected by dietary treatments, but total VFA concentration tended to be lower with NDF (84.7 vs. 74.1 mmol/l; P = 0.095). The amount of total bacteria tended (P = 0.075) to increase with NDSF, but only in diets with 0.35 NDF. The caecal proportions of Ruminococcus albus and Fibrobacter succinogenes were not affected by type or level of fibre, but Butyrivibrio fibrisolvens decreased (P = 0.055) with the NDF proportion in LS diets. Denaturing gradient gel electrophoresis (DGGE) analysis showed that bacterial communities clustered according to each combination of NDF and NDSF, but did not greatly differ among diets (similarity indexes between 0.67 and 0.70), nor biodiversity was affected (average Shannon and richness indexes 3.50 and 33.1; P > 0.10). Archaeal population revealed changes in the amount and composition that were particularly evident in HS diets, decreasing in concentration (from 4.37 to 4.12 log10 gene copy number/g) and biodiversity (Shannon index from 3.14 to 2.52 and richness index from 23.7 to 13.9) compared to LS. The type and level of dietary fibre had a minor impact on caecal fermentation traits or caecal bacterial community. However, the increase in NDSF from 0.14 to 0.17 reduced concentration and diversity of methanogenic archaea.

Comparative study of changes in reaction profile and microbial community structure in two anaerobic repeated-batch reactors started up with different seed sludges

Microbial community structure and dynamics were examined in two anaerobic reactors run in repeated-batch mode to treat whey permeate. Despite being started up using different seeding sources, the reactors showed generally similar reaction patterns and performances. During the repeated-batch operation for three cycles, the overall reaction rate increased with the increase in the initial population size of both bacteria and methanogens over cycles. Clostridium- and Methanospirillum-related microorganisms were likely the main acidogenic and methanogenic populations, respectively, in both reactors. Bacterial community structure shifted dynamically over cycles, while little change was observed in methanogen community structure throughout the operation. This means that the changes in bacterial community structure changes had little influence on the formation and evolution of methanogen community structure in the reactors. The increased methanogenesis rate with cycles seemed therefore more likely due to the effect of the increase in methanogen abundance rather than the alteration of community structure.

Lovastatin-enriched rice straw enhances biomass quality and suppresses ruminal methanogenesis

The primary objective of this study was to test the hypothesis that solid state fermentation (SSF) of agro-biomass (using rice straw as model); besides, breaking down its lignocellulose content to improve its nutritive values also produces lovastatin which could be used to suppress methanogenesis in the rumen ecosystem. Fermented rice straw (FRS) containing lovastatin after fermentation with Aspergillus terreus was used as substrate for growth study of rumen microorganisms using in vitro gas production method. In the first experiment, the extract from the FRS (FRSE) which contained lovastatin was evaluated for its efficacy for reduction in methane (CH4) production, microbial population, and activity in the rumen fluid. FRSE reduced total gas and CH4 productions (P < 0.01). It also reduced (P < 0.01) total methanogens population and increased the cellulolytic bacteria including Ruminococcus albus, Fibrobacter succinogenes (P < 0.01), and Ruminococcus flavefaciens (P < 0.05). Similarly, FRS reduced total gas and CH4 productions, methanogens population, but increased in vitro dry mater digestibility compared to the non-fermented rice straw. Lovastatin in the FRSE and the FRS significantly increased the expression of HMG-CoA reductase gene that produces HMG-CoA reductase, a key enzyme for cell membrane production in methanogenic Archaea.

Repeated pulse feeding induces functional stability in anaerobic digestion

Anaerobic digestion is an environmental key technology in the future bio-based economy. To achieve functional stability, a minimal microbial community diversity is required. This microbial community should also have a certain ‘elasticity’, i.e. the ability to rapidly adapt to suboptimal conditions or stress. In this study it was evaluated whether a higher degree of functional stability could be achieved by changing the feeding pattern, which can change the evenness, dynamics and richness of the bacterial community. The first reactor (CSTR_stable) was fed on daily basis, whereas the second reactor (CSTR_dynamic) was fed every 2 days. Average biogas production was 0.30 l CH₄ l⁻¹ day⁻¹ in both reactors, although daily variation was up to four times higher in the CSTR_dynamic compared with the CSTR_stable during the first 50 days. Bacterial analysis revealed that this CSTR_dynamic had a two times higher degree of bacterial community dynamics. The CSTR_dynamic also appeared to be more tolerant to an organic shock load of 8 g COD l⁻¹ and ammonium levels up to 8000 mg TAN l⁻¹. These results suggest that the regular application of a limited pulse of organic material and/or a variation in the substrate composition might promote higher functional stability in anaerobic digestion.

The role of anaerobic sludge recycle in improving anaerobic digester performance

Solids retention time (SRT) is a critical parameter for the performance of anaerobic digesters (AD) in wastewater treatment plants. AD SRT should increase when active biomass is input to the AD by recycling anaerobic sludge via the wastewater-treatment tanks, creating a hybrid aerobic/anaerobic system. When 85% of the flow through the AD was recycled in pilot-scale hybrid systems, the AD SRT increased by as much as 9-fold, compared to a parallel system without anaerobic-sludge recycle. Longer AD SRTs resulted in increased hydrolysis and methanogenesis in the AD: net solids yield decreased by 39-96% for overall and 23-94% in the AD alone, and AD methane yield increased 1.5- to 5.5-fold. Microbial community assays demonstrated higher, more consistent Archaea concentrations in all tanks in the wastewater-treatment system with anaerobic-sludge recycle. Thus, multiple lines of evidence support that AD-sludge recycle increased AD SRT, solids hydrolysis, and methane generation.

Nitrogen removal through different pathways in an aged refuse bioreactor treating mature landfill leachate

In this study, an aged refuse bioreactor was constructed to remove nitrogen in a mature landfill leachate. The nitrogen removal efficiency and the microbial community composition in the bioreactor were investigated. The results showed that the aged refuse bioreactor removed more than 90 % of total nitrogen in the leachate under the nitrogen loading rate (NLR) of 0.74 g/kg (vs) day, and the total nitrogen removal rate decreased to 62.2 % when NLR increased up to 2.03 g/kg (vs) day. Quantitative polymerase chain reaction results showed that the average cell number of ammonia-oxidizing bacteria in the bioreactor was 1.58 × 10(8) cells/g, which accounted for 0.41 % of total bacteria. The number of anammox bacteria in the reactor was 1.09 × 10(8) cells/g, which accounted for 0.27 % of total bacteria. Isotopic (15)N tracing experiments showed that nearly 10 % of nitrogen was removed by anammox. High-throughout 454 pyrosequencing revealed that the predominant bacteria in the bioreactor were Proteobacteria, Chloroflexi, Actinobacteria, Bacteroidetes, and Gemmatimonadetes, including various nitrifiers and denitrifiers with diverse heterotrophic and autotrophic metabolic pathways, supporting that nitrogen was removed through different pathways in this aged refuse bioreactor.

Enterohemorrhagic Escherichia coli O157:H7 survival in an in vitro model of the human large intestine and interactions with probiotic yeasts and resident microbiota

This is the first report on the fate of enterohemorrhagic Escherichia coli O157:H7 in simulated human colonic conditions. The pathogen was progressively eliminated from the bioreactor and did not modify the major populations of resident microbiota. The coadministration of the Saccharomyces cerevisiae CNCM I-3856 probiotic strain led to a significant increase in acetate production but did not reduce pathogen viability.

Methane production by treating vinasses from hydrous ethanol using a modified UASB reactor

BACKGROUND

A modified laboratory-scale upflow anaerobic sludge blanket (UASB) reactor was used to obtain methane by treating hydrous ethanol vinasse. Vinasses or stillage are waste materials with high organic loads, and a complex composition resulting from the process of alcohol distillation. They must initially be treated with anaerobic processes due to their high organic loads. Vinasses can be considered multipurpose waste for energy recovery and once treated they can be used in agriculture without the risk of polluting soil, underground water or crops. In this sense, treatment of vinasse combines the elimination of organic waste with the formation of methane. Biogas is considered as a promising renewable energy source. The aim of this study was to determine the optimum organic loading rate for operating a modified UASB reactor to treat vinasse generated in the production of hydrous ethanol from sugar cane molasses.

RESULTS

The study showed that chemical oxygen demand (COD) removal efficiency was 69% at an optimum organic loading rate (OLR) of 17.05 kg COD/m3-day, achieving a methane yield of 0.263 m3/kg CODadded and a biogas methane content of 84%. During this stage, effluent characterization presented lower values than the vinasse, except for potassium, sulfide and ammonia nitrogen. On the other hand, primers used to amplify the 16S-rDNA genes for the domains Archaea and Bacteria showed the presence of microorganisms which favor methane production at the optimum organic loading rate.

CONCLUSIONS

The modified UASB reactor proposed in this study provided a successful treatment of the vinasse obtained from hydrous ethanol production.Methanogen groups (Methanobacteriales and Methanosarcinales) detected by PCR during operational optimum OLR of the modified UASB reactor, favored methane production.

Bacterial profile from caecal contents and soft faeces in growing rabbits given diets differing in soluble and insoluble fibre levels

To verify if non-invasive collection of soft faeces (SF) from rabbits can be used as an index of bacterial biodiversity in caecal contents (CC), 24 weaned rabbits were given diets with low (LI) and high (HI) levels of insoluble fibre (neutral detergent fibre, NDF) and low (LS) and high (HS) levels of soluble fibre (neutral detergent soluble fibre, NDSF). After 21 days, animals were fitted with neck collars for SF collection. Two days later, animals were slaughtered and CC sampled. Total bacterial concentration quantified by real time PCR (log(10) ng DNA/mg DM) was higher in SF than CC (2.615 vs. 2.383). Among diets, in CC it was (P = 0.059) lowest in LILS diet, whereas in SF it decreased (P = 0.025) with the NDF level. DGGE profiles showed that structure of bacterial communities of SF was close to that of CC; however, similarity was higher in LI than HI diets (0.82 vs. 0.74). Diversity indexes in CC decreased with NDSF (P < 0.05), whereas the effect of NDF (P < 0.05) was also appreciated in SF. Soft faeces can be an alternative to surgery or slaughter techniques to monitor changes in caecal bacterial community; however, high dietary NDF may decrease similarity between both communities.

Gene number determination and genetic polymorphism of the gamma delta T cell co-receptor WC1 genes

Background

WC1 co-receptors belong to the scavenger receptor cysteine-rich (SRCR) superfamily and are encoded by a multi-gene family. Expression of particular WC1 genes defines functional subpopulations of WC1⁺ γδ T cells. We have previously identified partial or complete genomic sequences for thirteen different WC1 genes through annotation of the bovine genome Btau_3.1 build. We also identified two WC1 cDNA sequences from other cattle that did not correspond to sequences in the Btau_3.1 build. Their absence in the Btau_3.1 build may have reflected gaps in the genome assembly or polymorphisms among animals. Since the response of γδ T cells to bacterial challenge is determined by WC1 gene expression, it was critical to understand whether individual cattle or breeds differ in the number of WC1 genes or display polymorphisms.

Results

Real-time quantitative PCR using DNA from the animal whose genome was sequenced (“Dominette”) and sixteen other animals representing ten breeds of cattle, showed that the number of genes coding for WC1 co-receptors is thirteen. The complete coding sequences of those thirteen WC1 genes is presented, including the correction of an error in the WC1-2 gene due to mis-assembly in the Btau_3.1 build. All other cDNA sequences were found to agree with the previous annotation of complete or partial WC1 genes. PCR amplification and sequencing of the most variable N-terminal SRCR domain (domain 1 which has the SRCR “a” pattern) of each of the thirteen WC1 genes showed that the sequences are highly conserved among individuals and breeds. Of 160 sequences of domain 1 from three breeds of cattle, no additional sequences beyond the thirteen described WC1 genes were found. Analysis of the complete WC1 cDNA sequences indicated that the thirteen WC1 genes code for three distinct WC1 molecular forms.

Conclusion

The bovine WC1 multi-gene family is composed of thirteen genes coding for three structural forms whose sequences are highly conserved among individual cattle and breeds. The sequence diversity necessary for WC1 genes to function as a multi-genic pattern recognition receptor array is encoded in the genome, rather than generated by recombinatorial diversity or hypermutation.

Transfer of intestinal microbiota from lean donors increases insulin sensitivity in individuals with metabolic syndrome

Alterations in intestinal microbiota are associated with obesity and insulin resistance. We studied the effects of infusing intestinal microbiota from lean donors to male recipients with metabolic syndrome on the recipients' microbiota composition and glucose metabolism. Subjects were assigned randomly to groups that were given small intestinal infusions of allogenic or autologous microbiota. Six weeks after infusion of microbiota from lean donors, insulin sensitivity of recipients increased (median rate of glucose disappearance changed from 26.2 to 45.3 μmol/kg/min; P < .05) along with levels of butyrate-producing intestinal microbiota. Intestinal microbiota might be developed as therapeutic agents to increase insulin sensitivity in humans; www.trialregister.nl; registered at the Dutch Trial Register (NTR1776).

Methanogenic potential of formate in thermophilic anaerobic digestion

In the present study the methanogenic potential of formate (HCOO(-)) during thermophilic anaerobic digestion was investigated. After appropriate conditions for methanogenesis (HCOO(-) and inoculum concentration, pH and duration of incubation) were assessed, an experiment with initial 31 replicates was run. Diluted fermenter sludge was used as inoculum, and process parameters including the pH, quality and quantity of the produced biogas and the concentrations of volatile fatty acids and HCO(3) (-) were determined. Remarkably, after 5 days of incubation the highest CH(4) production was calculated for a HCOO(-) concentration of 200 mmol L(-1), a concentration, however, which might not occur in situ. During the phase of high CH(4) production HCOO(-) was degraded with a rate of 1.5 mmol L(-1) h(-1), and distinct changes of Gibbs free energy for several reactions were observed. Based on denaturing high-performance liquid chromatography, denaturing gradient gel electrophoresis, and additional subsequent sequencing approaches the hydrogenotrophic Methanothermobacter wolfeii was the dominant methanogen responsible for CH(4) production. Further confirmation was achieved due to the detection of autofluorescing rods with a size of up to ~3 µm, which were often arranged in pairs and chains. It was shown that even high concentrations of HCOO(-) are readily degraded, which might lead to an underestimation of both, the concentration and thus, the importance of HCOO(-) in anaerobic digestion.

Culture-dependent and -independent investigations of microbial diversity on urinary catheters

Catheter-associated urinary tract infection is caused by bacteria, which ascend the catheter along its external or internal surface to the bladder and subsequently develop into biofilms on the catheter and uroepithelium. Antibiotic-treated bacteria and bacteria residing in biofilm can be difficult to culture. In this study we used culture-based and 16S rRNA gene-based culture-independent methods (fingerprinting, cloning, and pyrosequencing) to determine the microbial diversity of biofilms on 24 urinary catheters. Most of the patients were catheterized for <30 days and had undergone recent antibiotic treatment. In addition, the corresponding urine samples for 16 patients were cultured. We found that gene analyses of the catheters were consistent with cultures of the corresponding urine samples for the presence of bacteria but sometimes discordant for the identity of the species. Cultures of catheter tips detected bacteria more frequently than urine cultures and gene analyses; coagulase-negative staphylococci were, in particular, cultured much more often from catheter tips, indicating potential contamination of the catheter tips during sampling. The external and internal surfaces of 19 catheters were separately analyzed by molecular methods, and discordant results were found in six catheters, suggesting that bacterial colonization intra- and extraluminally may be different. Molecular analyses showed that most of the species identified in this study were known uropathogens, and infected catheters were generally colonized by one to two species, probably due to antibiotic usage and short-term catheterization. In conclusion, our data showed that culture-independent molecular methods did not detect bacteria from urinary catheters more frequently than culture-based methods.

Role of bicarbonate as a pH buffer and electron sink in microbial dechlorination of chloroethenes

BACKGROUND

Buffering to achieve pH control is crucial for successful trichloroethene (TCE) anaerobic bioremediation. Bicarbonate (HCO3-) is the natural buffer in groundwater and the buffer of choice in the laboratory and at contaminated sites undergoing biological treatment with organohalide respiring microorganisms. However, HCO3- also serves as the electron acceptor for hydrogenotrophic methanogens and hydrogenotrophic homoacetogens, two microbial groups competing with organohalide respirers for hydrogen (H2). We studied the effect of HCO3- as a buffering agent and the effect of HCO3--consuming reactions in a range of concentrations (2.5-30 mM) with an initial pH of 7.5 in H2-fed TCE reductively dechlorinating communities containing Dehalococcoides, hydrogenotrophic methanogens, and hydrogenotrophic homoacetogens.

RESULTS

Rate differences in TCE dechlorination were observed as a result of added varying HCO3- concentrations due to H2-fed electrons channeled towards methanogenesis and homoacetogenesis and pH increases (up to 8.7) from biological HCO3- consumption. Significantly faster dechlorination rates were noted at all HCO3- concentrations tested when the pH buffering was improved by providing 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) as an additional buffer. Electron balances and quantitative PCR revealed that methanogenesis was the main electron sink when the initial HCO3- concentrations were 2.5 and 5 mM, while homoacetogenesis was the dominant process and sink when 10 and 30 mM HCO3- were provided initially.

CONCLUSIONS

Our study reveals that HCO3- is an important variable for bioremediation of chloroethenes as it has a prominent role as an electron acceptor for methanogenesis and homoacetogenesis. It also illustrates the changes in rates and extent of reductive dechlorination resulting from the combined effect of electron donor competition stimulated by HCO3- and the changes in pH exerted by methanogens and homoacetogens.

Microbial community structure and dynamics in anaerobic fluidized-bed and granular sludge-bed reactors: influence of operational temperature and reactor configuration

Methanogenic community structure and dynamics were investigated in two different, replicated anaerobic wastewater treatment reactor configurations [inverted fluidized bed (IFB) and expanded granular sludge bed (EGSB)] treating synthetic dairy wastewater, during operating temperature transitions from 37°C to 25°C, and from 25°C to 15°C, over a 430-day trial. Non-metric multidimensional scaling (NMS) and moving-window analyses, based on quantitative real-time PCR data, along with denaturing gradient gel electrophoresis (DGGE) profiling, demonstrated that the methanogenic communities developed in a different manner in these reactor configurations. A comparable level of performance was recorded for both systems at 37°C and 25°C, but a more dynamic and diverse microbial community in the IFB reactors supported better stability and adaptative capacity towards low temperature operation. The emergence and maintenance of particular bacterial genotypes (phylum Firmicutes and Bacteroidetes) was associated with efficient protein hydrolysis in the IFB, while protein hydrolysis was inefficient in the EGSB. A significant community shift from a Methanobacteriales and Methanosaetaceae towards a Methanomicrobiales-predominated community was demonstrated during operation at 15°C in both reactor configurations.

Stratification of Archaea in the deep sediments of a freshwater meromictic lake: vertical shift from methanogenic to uncultured archaeal lineages

As for lineages of known methanogens, several lineages of uncultured archaea were recurrently retrieved in freshwater sediments. However, knowledge is missing about how these lineages might be affected and structured according to depth. In the present study, the vertical changes of archaeal communities were characterized in the deep sediment of the freshwater meromictic Lake Pavin. For that purpose, an integrated molecular approach was performed to gain information on the structure, composition, abundance and vertical stratification of archaeal communities thriving in anoxic freshwater sediments along a gradient of sediments encompassing 130 years of sedimentation. Huge changes occurred in the structure and composition of archaeal assemblages along the sediment core. Methanogenic taxa (i.e. Methanosaeta and Methanomicrobiales) were progressively replaced by uncultured archaeal lineages (i.e. Marine Benthic Group-D (MBG-D) and Miscellaneous Crenarchaeal Group (MCG)) which are suspected to be involved in the methane cycle.

Selection for Cu-tolerant bacterial communities with altered composition, but unaltered richness, via long-term Cu exposure

Toxic metal pollution affects the composition and metal tolerance of soil bacterial communities. However, there is virtually no knowledge concerning the responses of members of specific bacterial taxa (e.g., phyla or classes) to metal toxicity, and contradictory results have been obtained regarding the impact of metals on operational taxonomic unit (OTU) richness. We used tag-coded pyrosequencing of the 16S rRNA gene to elucidate the impacts of copper (Cu) on bacterial community composition and diversity within a well-described Cu gradient (20 to 3,537 μg g(-1)) stemming from industrial contamination with CuSO(4) more than 85 years ago. DNA sequence information was linked to analysis of pollution-induced community tolerance (PICT) to Cu, as determined by the [(3)H]leucine incorporation technique, and to chemical characterization of the soil. PICT was significantly correlated to bioavailable Cu, as determined by the results seen with a Cu-specific bioluminescent biosensor strain, demonstrating a specific community response to Cu. The relative abundances of members of several phyla or candidate phyla, including the Proteobacteria, Bacteroidetes, Verrumicrobia, Chloroflexi, WS3, and Planctomycetes, decreased with increasing bioavailable Cu, while members of the dominant phylum, the Actinobacteria, showed no response and members of the Acidobacteria showed a marked increase in abundance. Interestingly, changes in the relative abundances of classes frequently deviated from the responses of the phyla to which they belong. Despite the apparent Cu impacts on Cu resistance and community structure, bioavailable Cu levels did not show any correlation to bacterial OTU richness (97% similarity level). Our report highlights several bacterial taxa responding to Cu and thereby provides new guidelines for future studies aiming to explore the bacterial domain for members of metal-responding taxa.

Early life treatment with vancomycin propagates Akkermansia muciniphila and reduces diabetes incidence in the NOD mouse

AIMS/HYPOTHESIS

Increasing evidence suggests that environmental factors changing the normal colonisation pattern in the gut strongly influence the risk of developing autoimmune diabetes. The aim of this study was to investigate, both during infancy and adulthood, whether treatment with vancomycin, a glycopeptide antibiotic specifically directed against Gram-positive bacteria, could influence immune homeostasis and the development of diabetic symptoms in the NOD mouse model for diabetes.

METHODS

Accordingly, one group of mice received vancomycin from birth until weaning (day 28), while another group received vancomycin from 8 weeks of age until onset of diabetes. Pyrosequencing of the gut microbiota and flow cytometry of intestinal immune cells was used to investigate the effect of vancomycin treatment.

RESULTS

At the end of the study, the cumulative diabetes incidence was found to be significantly lower for the neonatally treated group compared with the untreated group, whereas the insulitis score and blood glucose levels were significantly lower for the mice treated as adults compared with the other groups. Mucosal inflammation was investigated by intracellular cytokine staining of the small intestinal lymphocytes, which displayed an increase in cluster of differentiation (CD)4(+) T cells producing pro-inflammatory cytokines in the neonatally treated mice. Furthermore, bacteriological examination of the gut microbiota composition by pyrosequencing revealed that vancomycin depleted many major genera of Gram-positive and Gram-negative microbes while, interestingly, one single species, Akkermansia muciniphila, became dominant.

CONCLUSIONS/INTERPRETATION

The early postnatal period is a critical time for microbial protection from type 1 diabetes and it is suggested that the mucolytic bacterium A. muciniphila plays a protective role in autoimmune diabetes development, particularly during infancy.

SELECTION OF INTERNAL CONTROL GENE FOR EXPRESSION STUDIES IN PORPHYRA HAITANENSIS (RHODOPHYTA) AT DIFFERENT LIFE-HISTORY STAGES(1)

Accurate gene quantification depends on the use of an appropriate internal control gene, which should be verified before its use for normalizing data. Housekeeping genes, which are expressed at relatively constant levels, are generally regarded as candidate internal control genes. To determine the ideal internal control for gene expression profiles for Porphyra haitanensis T. J. Chang et B. F. Zheng (Bangiales, Rhodophyta) at different life-history stages, we used absolute quantification to assess the expression levels of six housekeeping genes (18S ribosomal RNA, 30S ribosomal protein, glyceraldehyde-3-phosphate dehydrogenase, elongation factor 3, alpha-tubulin, and beta-tubulin) at the sporophyte and gametophyte stages. Housekeeping genes were selected by comparing the differences of observed copy numbers in sporophytes and in gametophytes. TubB (beta-tubulin) was found to be the optimal internal control gene, because it showed the smallest difference of gene expression. Compared with TubB, other housekeeping genes had greater variation of expression to different degrees.

Effect of the rumen ciliates Entodinium caudatum, Epidinium ecaudatum and Eudiplodinium maggii, and combinations thereof, on ruminal fermentation and total tract digestion in sheep

The quantitative importance of individual ciliate species and their interaction in the rumen is still unclear. The present study was performed to test whether there are species differences in the influence on ruminal fermentation in vivo and if combinations of ciliates act additive in that respect. Six adult wethers fed a hay-concentrate diet were defaunated, then refaunated either with Entodinium caudatum (EC), Epidinium ecaudatum (EE) or Eudiplodinium maggii (EM) alone, then progressively with all possible species combinations. Feed, faeces, urine, ruminal fluid and gas were sampled for eight days always after at least 21 days of adaptation. With a linear mixed model, accounting for the 2 x 2 x 2 full factorial study design, mean marginal effect sizes, i.e., the magnitude of change in variables as caused by the presence of each ciliate species or of combinations of them, were estimated. The apparent digestibility of organic matter and neutral detergent fibre remained unaffected. The apparent N digestibility increased by 0.054 with EM (0.716 with defaunation). Ruminal ammonia increased by 1.6, 4.0 and 8.7 mmol/l in the presence of EM, EC and EE, respectively, compared to defaunation (6.9 mmol/l). In the EM + EE combination, ruminal ammonia was lower than would have been expected from an additive effect. With EE, total short-chain fatty acids increased by 23 mmol/l (100 mmol/l with defaunation), but not when EE was combined with EM. The acetate-to-propionate ratio decreased by 0.73 units in the presence of EE (4.0 with defaunation), but only when EE was the sole ciliate species in the rumen. In the presence of any ciliate species, the 16S rDNA copies of total Bacteria and major fibrolytic species decreased to 0.52- and 0.22-fold values, respectively of that found without protozoa. Total Archaea were unaffected; however, Methanobacteriales copies increased 1.44-fold with EC. The CH4-to-CO2 ratio of ruminal gas decreased by 0.036 with EM and 0.051 with EE (0.454 with defaunation). In conclusion, individual ciliates affected ruminal fermentation differently and, when different species were combined, sometimes in a non-additive manner. From the ciliates investigated, EE affected ruminal fermentation most and might play a dominant role in mixed ciliate populations.