Unusual Metabolism and Hypervariation in the Genome of a Gracilibacterium (BD1-5) from an Oil-Degrading Community

Bioinformatic discovery of type 11 secretion system (T11SS) cargo across the Proteobacteria

Type 11 secretion systems (T11SS) are broadly distributed amongst Proteobacteria, with more than 3,000 T11SS family outer membrane proteins (OMPs) comprising ten major sequence similarity network clusters. Of these, only seven, all from animal-associated cluster 1, have been experimentally verified as secretins of cargo, including adhesins, haemophores and metal-binding proteins. To identify novel cargo of a more diverse set of T11SS, we identified gene families co-occurring in gene neighbourhoods with either cluster 1 or marine microbe-associated cluster 3 T11SS OMP genes. We developed bioinformatic controls to ensure that perceived co-occurrences are specific to T11SS, and not general to OMPs. We found that both cluster 1 and cluster 3 T11SS OMPs frequently co-occur with single-carbon metabolism and nucleotide synthesis pathways, but that only cluster 1 T11SS OMPs had significant co-occurrence with metal and haem pathways, as well as with mobile genetic islands, potentially indicating the diversified function of this cluster. Cluster 1 T11SS co-occurrences included 2,556 predicted cargo proteins, unified by the presence of a C-terminal β-barrel domain, which fall into 141 predicted UniRef50 clusters and approximately ten different architectures: four similar to known cargo and six uncharacterized types. We experimentally demonstrate T11SS-dependent secretion of an uncharacterized cargo type with homology to plasmin-sensitive protein. Unexpectedly, genes encoding marine cluster 3 T11SS OMPs only rarely co-occurred with the C-terminal β-barrel domain and instead frequently co-occurred with DUF1194-containing genes. Overall, our results show that with sufficiently large-scale and controlled genomic data, T11SS-dependent cargo proteins can be accurately predicted.

Metaproteogenomics resolution of a high-CO2 aquifer community reveals a complex cellular adaptation of groundwater Gracilibacteria to a host-dependent lifestyle

BACKGROUND

Bacteria of the candidate phyla radiation (CPR), constituting about 25% of the bacterial biodiversity, are characterized by small cell size and patchy genomes without complete key metabolic pathways, suggesting a symbiotic lifestyle. Gracilibacteria (BD1-5), which are part of the CPR branch, possess alternate coded genomes and have not yet been cultivated. The lifestyle of Gracilibacteria, their temporal dynamics, and activity in natural ecosystems, particularly in groundwater, has remained largely unexplored. Here, we aimed to investigate Gracilibacteria activity in situ and to discern their lifestyle based on expressed genes, using the metaproteogenome of Gracilibacteria as a function of time in the cold-water geyser Wallender Born in the Volcanic Eifel region in Germany.

RESULTS

We coupled genome-resolved metagenomics and metaproteomics to investigate a cold-water geyser microbial community enriched in Gracilibacteria across a 12-day time-series. Groundwater was collected and sequentially filtered to fraction CPR and other bacteria. Based on 725 Gbps of metagenomic data, 1129 different ribosomal protein S3 marker genes, and 751 high-quality genomes (123 population genomes after dereplication), we identified dominant bacteria belonging to Gallionellales and Gracilibacteria along with keystone microbes, which were low in genomic abundance but substantially contributing to proteomic abundance. Seven high-quality Gracilibacteria genomes showed typical limitations, such as limited amino acid or nucleotide synthesis, in their central metabolism but no co-occurrence with potential hosts. The genomes of these Gracilibacteria were encoded for a high number of proteins involved in cell to cell interaction, supporting the previously surmised host-dependent lifestyle, e.g., type IV and type II secretion system subunits, transporters, and features related to cell motility, which were also detected on protein level.

CONCLUSIONS

We here identified microbial keystone taxa in a high-CO2 aquifer, and revealed microbial dynamics of Gracilibacteria. Although Gracilibacteria in this ecosystem did not appear to target specific organisms in this ecosystem due to lack of co-occurrence despite enrichment on 0.2-µm filter fraction, we provide proteomic evidence for the complex machinery behind the host-dependent lifestyle of groundwater Gracilibacteria. Video Abstract.

Refocusing the microbial rare biosphere concept through a functional lens

The influential concept of the rare biosphere in microbial ecology has underscored the importance of taxa occurring at low abundances yet potentially playing key roles in communities and ecosystems. Here, we refocus the concept of rare biosphere through a functional trait-based lens and provide a framework to characterize microbial functional rarity, a combination of numerical scarcity across space or time and trait distinctiveness. We demonstrate how this novel interpretation of the rare biosphere, rooted in microbial functions, can enhance our mechanistic understanding of microbial community structure. It also sheds light on functionally distinct microbes, directing conservation efforts towards taxa harboring rare yet ecologically crucial functions.

Interpretable metric learning in comparative metagenomics: The adaptive Haar-like distance

Random forests have emerged as a promising tool in comparative metagenomics because they can predict environmental characteristics based on microbial composition in datasets where β-diversity metrics fall short of revealing meaningful relationships between samples. Nevertheless, despite this efficacy, they lack biological insight in tandem with their predictions, potentially hindering scientific advancement. To overcome this limitation, we leverage a geometric characterization of random forests to introduce a data-driven phylogenetic β-diversity metric, the adaptive Haar-like distance. This new metric assigns a weight to each internal node (i.e., split or bifurcation) of a reference phylogeny, indicating the relative importance of that node in discerning environmental samples based on their microbial composition. Alongside this, a weighted nearest-neighbors classifier, constructed using the adaptive metric, can be used as a proxy for the random forest while maintaining accuracy on par with that of the original forest and another state-of-the-art classifier, CoDaCoRe. As shown in datasets from diverse microbial environments, however, the new metric and classifier significantly enhance the biological interpretability and visualization of high-dimensional metagenomic samples.

Metabolic implications for predatory and parasitic bacterial lineages in activated sludge wastewater treatment systems

Deciphering unclear microbial interactions is key to improving biological wastewater treatment processes. Microbial predation and parasitism in wastewater treatment ecosystems are unexplored survival strategies that have long been known and have recently attracted attention because these interspecies interactions may contribute to the reduction of excess sludge. Here, microbial community profiling of 600 activated sludge samples taken from six industrial and one municipal wastewater treatment processes (WWTPs) was conducted. To identify the shared lineages in the WWTPs, the shared microbial constituents were defined as the family level taxa that had ≥ 0.1% average relative abundance and detected in all processes. The microbial community analysis assigned 106 families as the shared microbial constituents in the WWTPs. Correlation analysis showed that 98 of the 106 shared families were significantly correlated with total carbon (TC) and/or total nitrogen (TN) concentrations, suggesting that they may contribute to wastewater remediation. Most possible predatory or parasitic bacteria belonging to the phyla Bdellovibrionota, Myxococcota, and Candidatus Patescibacteria were found to be the shared families and negatively correlated with TC/TN; thus, they were frequently present in the WWTPs and could be involved in the removal of carbon/nitrogen derived from cell components. Shotgun metagenome-resolved metabolic reconstructions indicated that gene homologs associated with predation or parasitism are conserved in the Bdellovibrionota, Myxococcota, and Ca. Patescibacteria genomes (e.g., host interaction (hit) locus, Tad-like secretion complexes, and type IV pilus assembly proteins). This study provides insights into the complex microbial interactions potentially linked to the reduction of excess sludge biomass in these processes.

Infant Saliva Microbiome Activity Modulates Nutritional Impacts on Neurodevelopment

Neurodevelopment is influenced by complex interactions between environmental factors, including social determinants of health (SDOH), nutrition, and even the microbiome. This longitudinal cohort study of 142 infants tested the hypothesis that microbial activity modulates the effects of nutrition on neurodevelopment. Salivary microbiome activity was measured at 6 months using RNA sequencing. Infant nutrition was assessed longitudinally with the Infant Feeding Practices survey. The primary outcome was presence/absence of neurodevelopmental delay (NDD) at 18 months on the Survey of Wellbeing in Young Children. A logistic regression model employing two microbial factors, one nutritional factor, and two SDOH accounted for 33.3% of the variance between neurodevelopmental groups (p < 0.001, AIC = 77.7). NDD was associated with Hispanic ethnicity (OR 18.1, 2.36-139.3; p = 0.003), no fish consumption (OR 10.6, 2.0-54.1; p = 0.003), and increased Candidatus Gracilibacteria activity (OR 1.43, 1.00-2.07; p = 0.007). Home built after 1977 (OR 0.02, 0.001-0.53; p = 0.004) and Chlorobi activity (OR 0.76, 0.62-0.93, p = 0.001) were associated with reduced risk of NDD. Microbial alpha diversity modulated the effect of fish consumption on NDD (X2 = 5.7, p = 0.017). These data suggest the benefits of fish consumption for neurodevelopment may be mediated by microbial diversity. Confirmation in a larger, randomized trial is required.

Using Genomics to Decipher the Enigmatic Properties and Survival Adaptation of Candidate Phyla Radiation

Microbial ecology is a critical field for understanding the composition, diversity, and functions of microorganisms in various environmental and health-related processes. The discovery of Candidate Phyla Radiation (CPR) through culture-independent methods has introduced a new division of microbes characterized by a symbiotic/parasitic lifestyle, small cell size, and small genome. Despite being poorly understood, CPRs have garnered significant attention in recent years due to their widespread detection in a variety of environmental and clinical samples. These microorganisms have been found to exhibit a high degree of genetic diversity compared to other microbes. Several studies have shed light on their potential importance in global biogeochemical cycles and their impact on various human activities. In this review, we provide a systematic overview of the discovery of CPRs. We then focus on describing how the genomic characteristics of CPRs have helped them interact with and adapt to other microbes in different ecological niches. Future works should focus on discovering the metabolic capacities of CPRs and, if possible, isolating them to obtain a better understanding of these microorganisms.

Responses of Free-Living Planktonic Bacterial Communities to Experimental Acidification and Warming

Climate change driven by human activities encompasses the increase in atmospheric CO₂ concentration and sea-surface temperature. Little is known regarding the synergistic effects of these phenomena on bacterial communities in oligotrophic marine ecosystems that are expected to be particularly vulnerable. Here, we studied bacterial community composition changes based on 16S rRNA sequencing at two fractions (0.1-0.2 and >0.2 μm) during a 10- day fully factorial mesocosm experiment in the eastern Mediterranean where the pH decreased by ~0.3 units and temperature increased by ~3 °C to project possible future changes in surface waters. The bacterial community experienced significant taxonomic differences driven by the combined effect of time and treatment; a community shift one day after the manipulations was noticed, followed by a similar state between all mesocosms at the third day, and mild shifts later on, which were remarkable mainly under sole acidification. The abundance of Synechococcus increased in response to warming, while the SAR11 clade immediately benefited from the combined acidification and warming. The effect of the acidification itself had a more persistent impact on community composition. This study highlights the importance of studying climate change consequences on ecosystem functioning both separately and simultaneously, considering the ambient environmental parameters.

Global patterns of diversity and metabolism of microbial communities in deep-sea hydrothermal vent deposits

BACKGROUND

When deep-sea hydrothermal fluids mix with cold oxygenated fluids, minerals precipitate out of solution and form hydrothermal deposits. These actively venting deep-sea hydrothermal deposits support a rich diversity of thermophilic microorganisms which are involved in a range of carbon, sulfur, nitrogen, and hydrogen metabolisms. Global patterns of thermophilic microbial diversity in deep-sea hydrothermal ecosystems have illustrated the strong connectivity between geological processes and microbial colonization, but little is known about the genomic diversity and physiological potential of these novel taxa. Here we explore this genomic diversity in 42 metagenomes from four deep-sea hydrothermal vent fields and a deep-sea volcano collected from 2004 to 2018 and document their potential implications in biogeochemical cycles.

RESULTS

Our dataset represents 3635 metagenome-assembled genomes encompassing 511 novel and recently identified genera from deep-sea hydrothermal settings. Some of the novel bacterial (107) and archaeal genera (30) that were recently reported from the deep-sea Brothers volcano were also detected at the deep-sea hydrothermal vent fields, while 99 bacterial and 54 archaeal genera were endemic to the deep-sea Brothers volcano deposits. We report some of the first examples of medium- (≥ 50% complete, ≤ 10% contaminated) to high-quality (> 90% complete, < 5% contaminated) MAGs from phyla and families never previously identified, or poorly sampled, from deep-sea hydrothermal environments. We greatly expand the novel diversity of Thermoproteia, Patescibacteria (Candidate Phyla Radiation, CPR), and Chloroflexota found at deep-sea hydrothermal vents and identify a small sampling of two potentially novel phyla, designated JALSQH01 and JALWCF01. Metabolic pathway analysis of metagenomes provides insights into the prevalent carbon, nitrogen, sulfur, and hydrogen metabolic processes across all sites and illustrates sulfur and nitrogen metabolic "handoffs" in community interactions. We confirm that Campylobacteria and Gammaproteobacteria occupy similar ecological guilds but their prevalence in a particular site is driven by shifts in the geochemical environment.

CONCLUSION

Our study of globally distributed hydrothermal vent deposits provides a significant expansion of microbial genomic diversity associated with hydrothermal vent deposits and highlights the metabolic adaptation of taxonomic guilds. Collectively, our results illustrate the importance of comparative biodiversity studies in establishing patterns of shared phylogenetic diversity and physiological ecology, while providing many targets for enrichment and cultivation of novel and endemic taxa. Video Abstract.

The Microbiome Structure of the Symbiosis between the Desert Truffle Terfezia boudieri and Its Host Plant Helianthemum sessiliflorum

The desert truffle Terfezia boudieri is an ascomycete fungus that forms ect-endomycorrhiza in the roots of plants belonging to Cistaceae. The fungus forms hypogeous edible fruit bodies, appreciated as gourmet food. Truffles and host plants are colonized by various microbes, which may contribute to their development. However, the diversity and composition of the bacterial community under field conditions in the Negev desert are still unknown. The overall goal of this research was to identify the rhizosphere microbial community supporting the establishment of a symbiotic association between T. boudieri and Helianthemum sessiliflorum. The bacterial community was characterized by fruiting bodies, mycorrhized roots, and rhizosphere soil. Based on next-generation sequencing meta-analyses of the 16S rRNA gene, we discovered diverse bacterial communities of fruit bodies that differed from those found in the roots and rhizosphere. Families of Proteobacteria, Planctomycetes, and Actinobacteria were present in all four samples. Alpha diversity analysis revealed that the rhizosphere and roots contain significantly higher bacterial species numbers compared to the fruit. Additionally, ANOSIM and PCoA provided a comparative analysis of the bacterial taxa associated with fruiting bodies, roots, and rhizosphere. The core microbiome described consists of groups whose biological role triggers important traits supporting plant growth and fruit body development.

Investigation of microbial metabolisms in an extremely high pH marine-like terrestrial serpentinizing system: Ney Springs

Ney Springs, a continental serpentinizing spring in northern California, has an exceptionally high reported pH (12.4) for a naturally occurring water source. With high conductivity fluids, it is geochemically more akin to marine serpentinizing systems than other terrestrial locations. Our geochemical analyses also revealed high sulfide concentrations (544 mg/L) and methane emissions (83% volume gas content) relative to other serpentinizing systems. Thermodynamic calculations were used to investigate the potential for substrates resulting from serpentinization to fuel microbial life, and were found to support the energetic feasibility of sulfate reduction, anaerobic methane oxidation, denitrification, and anaerobic sulfide oxidation within this system. Assessment of the microbial community via 16S rRNA taxonomic gene surveys and metagenome sequencing revealed a community composition dominated by poorly characterized members of the Izemoplasmatales and Clostridiales. The genomes of these dominant taxa point to a fermentative lifestyle, though other highly complete (>90%) metagenome assembled genomes support the potential for organisms to perform sulfate reduction, sulfur disproportionation and/or sulfur oxidation (aerobic and anaerobic). Two chemolithoheterotrophs identified in the metagenome, a Halomonas sp. and a Rhodobacteraceae sp., were isolated and shown to oxidize thiosulfate and were capable of growth in conditions up to pH 12.4. Despite being characteristic products of serpentinization reactions, little evidence was seen for hydrogen and methane utilization in the Ney Springs microbial community. Hydrogen is not highly abundant and could be consumed prior to reaching the spring community. Other metabolic strategies may be outcompeted by more energetically favorable heterotrophic or fermentation reactions, or even inhibited by other compounds in the spring such as ammonia. The unique geochemistry of Ney Springs provides an opportunity to study how local geology interacts with serpentinized fluids, while its microbial community can better inform us of the metabolic strategies employed in hyperalkaline environments.

Integrating high-throughput sequencing and metabolomics to investigate the stereoselective responses of soil microorganisms to chiral fungicide cis-epoxiconazole

The use of the chiral triazole fungicide cis-epoxiconazole in agricultural production continues to increase; however, little is known about the stereoselective and toxic responses of soil microorganisms to cis-epoxiconazole in the soil microenvironment. High-throughput sequencing and metabolomics were integrated to investigate the stereoselective response of soil microbial community structure, metabolic profile to cis-epoxiconazole exposure, and the correlation between the microbiomes and different metabolites. Soil microbial community structure and soil metabolic profile were significantly altered and exhibited significant enantioselectivity. The alpha diversity (Chao, Shannon, and Simpson diversity) of bacterial and fungus was not significantly affected, whereas the beta diversity (Bray-Curtis dissimilarity and PLS-DA) of bacterial and fungus was significantly altered in treatment of cis-epoxiconazole and its enantiomers (p-value < 0.05). The variation in bacterial and fungus community structure was the highest under (+)-enantiomer exposure, followed by exposure to racemate and (-)-enantiomer. Soil metabolomic analysis revealed that exposure to high or low doses of cis-epoxiconazole and its enantiomers resulted in different degrees of reprogramming of the soil metabolic pool. The 39 significantly changed metabolites mainly included small molecular organic acids, amino acids and their intermediates, and purine and adenosine intermediates. Six metabolic pathways were significantly disrupted. Different correlation patterns were observed between the significantly altered metabolites and microbes (p-value < 0.05) by Pearson correlation-based analysis. In conclusion, as xenobiotic pollutant, epoxiconazole altered the structure and metabolism of soil microorganisms with significant stereoselectivity mainly driven by 2R, 3S-(+)-cis-epoxiconazole. This study provided a more robust assessment of the risks of epoxiconazole exposure to soil microorganisms. Given the importance of the soil environment in agricultural production, characterization of the soil microbiome and metabolome can provide new insights into the ecological risks posed by exposure to the chiral triazole pesticide cis-epoxiconazole and its enantiomers.

Metabolic Potential of the Superphylum Patescibacteria Reconstructed from Activated Sludge Samples from a Municipal Wastewater Treatment Plant

Patescibacteria are widely distributed in various environments and often detected in activated sludge. However, limited information is currently available on their phylogeny, morphology, and ecophysiological role in activated sludge or interactions with other microorganisms. In the present study, we identified microorganisms that interacted with Patescibacteria in activated sludge via a correlation ana-lysis using the 16S rRNA gene, and predicted the metabolic potential of Patescibacteria using a metagenomic ana-lysis. The metagenome-assembled genomes of Patescibacteria consisted of three Saccharimonadia, three Parcubacteria, and one Gracilibacteria, and showed a strong positive correlation of relative abundance with Chitinophagales. Metabolic predictions from ten recovered patescibacterial and five Chitinophagales metagenome-assembled genomes supported mutualistic interactions between a member of Saccharimonadia and Chitinophagales via N-acetylglucosamine, between a member of Parcubacteria and Chitinophagales via nitrogen compounds related to denitrification, and between Gracilibacteria and Chitinophagales via phospholipids in activated sludge. The present results indicate that various interactions between Patescibacteria and Chitinophagales are important for the survival of Patescibacteria in activated sludge ecosystems.

Bacterial Necromass Is Rapidly Metabolized by Heterotrophic Bacteria and Supports Multiple Trophic Levels of the Groundwater Microbiome

Pristine groundwater is a highly stable environment with microbes adapted to dark, oligotrophic conditions. Input events like heavy rainfalls can introduce the excess particulate organic matter, including surface-derived microorganisms, thereby disturbing the groundwater microbiome. Some surface-derived bacteria will not survive this translocation, leading to an input of necromass to the groundwater. Here, we investigated the effects of necromass addition to the microbial community in fractured bedrock groundwater, using groundwater mesocosms as model systems. We followed the uptake of ¹³C-labeled necromass by the bacterial and eukaryotic groundwater community quantitatively and over time using a complementary protein-stable and DNA-stable isotope probing approach. Necromass was rapidly depleted in the mesocosms within 4 days, accompanied by a strong decrease in Shannon diversity and a 10-fold increase in bacterial 16S rRNA gene copy numbers. Species of Flavobacterium, Massilia, Rheinheimera, Rhodoferax, and Undibacterium dominated the microbial community within 2 days and were identified as key players in necromass degradation, based on a ¹³C incorporation of >90% in their peptides. Their proteomes comprised various proteins for uptake and transport functions and amino acid metabolization. After 4 and 8 days, the autotrophic and mixotrophic taxa Nitrosomonas, Limnohabitans, Paucibacter, and Acidovorax increased in abundance with a ¹³C incorporation between 0.5% and 23%. Likewise, eukaryotes assimilated necromass-derived carbon either directly or indirectly. Our data point toward a fast and exclusive uptake of labeled necromass by a few specialists followed by a concerted action of groundwater microorganisms, including autotrophs presumably fueled by released, reduced nitrogen and sulfur compounds generated during necromass degradation.

IMPORTANCE Subsurface microbiomes provide essential ecosystem services, like the generation of drinking water. These ecosystems are devoid of light-driven primary production, and microbial life is adapted to the resulting oligotrophic conditions. Modern groundwater is most vulnerable to anthropogenic and climatic impacts. Heavy rainfalls, which will increase with climate change, can result in high surface inputs into shallow aquifers by percolation or lateral flow. These inputs include terrestrial organic matter and surface-derived microbes that are not all capable to flourish in aquatic subsurface habitats. Here, we investigated the response of groundwater mesocosms to the addition of bacterial necromass, simulating event-driven surface input. We found that the groundwater microbiome responds with a rapid bloom of only a few primary degraders, followed by the activation of typical groundwater autotrophs and mixotrophs, as well as eukaryotes. Our results suggest that this multiphase strategy is essential to maintain the balance of the groundwater microbiome to provide ecosystem services.

Candidate Phyla Radiation, an Underappreciated Division of the Human Microbiome, and Its Impact on Health and Disease

Candidate phyla radiation (CPR) is an emerging division of the bacterial domain within the human microbiota. Still poorly known, these microorganisms were first described in the environment in 1981 as "ultramicrobacteria" with a cell volume under 0.1 μm³ and were first associated with the human oral microbiota in 2007. The evolution of technology has been paramount for the study of CPR within the human microbiota. In fact, since these ultramicrobacteria have yet to be axenically cultured despite ongoing efforts, progress in imaging technology has allowed their observation and morphological description. Although their genomic abilities and taxonomy are still being studied, great strides have been made regarding their taxonomic classification, as well as their lifestyle. In addition, advancements in next-generation sequencing and the continued development of bioinformatics tools have allowed their detection as commensals in different human habitats, including the oral cavity and gastrointestinal and genital tracts, thus highlighting CPR as a nonnegligible part of the human microbiota with an impact on physiological settings. Conversely, several pathologies present dysbiosis affecting CPR levels, including inflammatory, mucosal, and infectious diseases. In this exhaustive review of the literature, we provide a historical perspective on the study of CPR, an overview of the methods available to study these organisms and a description of their taxonomy and lifestyle. In addition, their distribution in the human microbiome is presented in both homeostatic and dysbiotic settings. Future efforts should focus on developing cocultures and, if possible, axenic cultures to obtain isolates and therefore genomes that would provide a better understanding of these ultramicrobacteria, the importance of which in the human microbiome is undeniable.

Unifying the known and unknown microbial coding sequence space

Genes of unknown function are among the biggest challenges in molecular biology, especially in microbial systems, where 40-60% of the predicted genes are unknown. Despite previous attempts, systematic approaches to include the unknown fraction into analytical workflows are still lacking. Here, we present a conceptual framework, its translation into the computational workflow AGNOSTOS and a demonstration on how we can bridge the known-unknown gap in genomes and metagenomes. By analyzing 415,971,742 genes predicted from 1749 metagenomes and 28,941 bacterial and archaeal genomes, we quantify the extent of the unknown fraction, its diversity, and its relevance across multiple organisms and environments. The unknown sequence space is exceptionally diverse, phylogenetically more conserved than the known fraction and predominantly taxonomically restricted at the species level. From the 71 M genes identified to be of unknown function, we compiled a collection of 283,874 lineage-specific genes of unknown function for Cand. Patescibacteria (also known as Candidate Phyla Radiation, CPR), which provides a significant resource to expand our understanding of their unusual biology. Finally, by identifying a target gene of unknown function for antibiotic resistance, we demonstrate how we can enable the generation of hypotheses that can be used to augment experimental data.

Microbial Community Structure and Ecological Networks during Simulation of Diatom Sinking

Microbial-mediated utilization of particulate organic matter (POM) during its downward transport from the surface to the deep ocean constitutes a critical component of the global ocean carbon cycle. However, it remains unclear as to how high hydrostatic pressure (HHP) and low temperature (LT) with the sinking particles affects community structure and network interactions of the particle-attached microorganisms (PAM) and those free-living microorganisms (FLM) in the surrounding water. In this study, we investigated microbial succession and network interactions in experiments simulating POM sinking in the ocean. Diatom-derived ¹³C- and ¹²C-labeled POM were used to incubate surface water microbial communities from the East China Sea (ECS) under pressure (temperature) of 0.1 (25 °C), 20 (4 °C), and 40 (4 °C) MPa (megapascal). Our results show that the diversity and species richness of the PAM and FLM communities decreased significantly with HHP and LT. Microbial community analysis indicated an increase in the relative abundance of Bacteroidetes at high pressure (40 MPa), mostly at the expense of Gammaproteobacteria, Alphaproteobacteria, and Gracilibacteria at atmospheric pressure. Hydrostatic pressure and temperature affected lifestyle preferences between particle-attached (PA) and free-living (FL) microbes. Ecological network analysis showed that HHP and LT enhanced microbial network interactions and resulted in higher vulnerability to networks of the PAM communities and more resilience of those of the FLM communities. Most interestingly, the PAM communities occupied most of the module hubs of the networks, whereas the FLM communities mainly served as connectors of the modules, suggesting their different ecological roles of the two groups of microbes. These results provided novel insights into how HHP and LT affected microbial community dynamics, ecological networks during POM sinking, and the implications for carbon cycling in the ocean.

Ecological Interactions of Cyanobacteria and Heterotrophs Enhances the Robustness of Cyanobacterial Consortium for Carbon Sequestration

Lack of robustness is a major barrier to foster a sustainable cyanobacterial biotechnology. Use of cyanobacterial consortium increases biodiversity, which provides functional redundancy and prevents invading species from disrupting the production ecosystem. Here we characterized a cyanobacterial consortium enriched from microbial mats of alkaline soda lakes in BC, Canada, at high pH and alkalinity. This consortium has been grown in open laboratory culture for 4 years without crashes. Using shotgun metagenomic sequencing, 29 heterotrophic metagenome-assembled-genomes (MAGs) were retrieved and were assigned to Bacteroidota, Alphaproteobacteria, Gammaproteobacteria, Verrucomicrobiota, Patescibacteria, Planctomycetota, and Archaea. In combination with metaproteomics, the overall stability of the consortium was determined under different cultivation conditions. Genome information from each heterotrophic population was investigated for six ecological niches created by cyanobacterial metabolism and one niche for phototrophy. Genome-resolved metaproteomics with stable isotope probing using 13C-bicarbonate (protein/SIP) showed tight coupling of carbon transfer from cyanobacteria to the heterotrophic populations, specially Wenzhouxiangella. The community structure was compared to a previously described consortium of a closely related cyanobacteria, which indicated that the results may be generalized. Productivity losses associated with heterotrophic metabolism were relatively small compared to other losses during photosynthesis.

Dominance of Arcobacter in the white filaments from the thermal sulfidic spring of Fetida Cave (Apulia, southern Italy)

The thermal spring of Fetida Cave, a still active sulfuric acid cave opening at sea level and located in Santa Cesarea Terme, southeastern Salento (Apulia region, Southern Italy) hosts abundant floating white filaments. The white filaments were mainly composed of sulfur crystals surrounded by microbial mass of the phyla Epsilonbacteraeota, Proteobacteria, Bacteroidetes, and Patescibacteria. The most abundant genus in the white filaments collected from the waters in the innermost part of the cave dominated by sulfidic exhalations was Arcobacter. This abundance can be related to the higher concentration of sulfide dissolved in water, and low oxygen and pH values. Conversely, lower Arcobacter abundances were obtained in the filaments collected in the entrance and middle part of the cave, where sulfidic water mixes with seawater, as the cave is subjected to tides and the mixing of fresh (continental) with marine water. The geochemical analysis of water and atmospheric gases confirmed these environmental constraints. In fact, the highest concentrations of H₂S in the air and water were recorded closest to the spring upwelling in the innermost part of the cave, and the lowest ones near the cave entrance. The metabolic versatility of Arcobacter might provide a competitive advantage in the colonization of water bodies characterized by high sulfide, low oxygen, and dynamic fluid movement.

Cultivation of a vampire: 'Candidatus Absconditicoccus praedator'

Halorhodospira halophila, one of the most-xerophilic halophiles, inhabits biophysically stressful and energetically expensive, salt-saturated alkaline brines. Here, we report an additional stress factor that is biotic: a diminutive Candidate-Phyla-Radiation bacterium, that we named 'Ca. Absconditicoccus praedator' M39-6, which predates H. halophila M39-5, an obligately photosynthetic, anaerobic purple-sulfur bacterium. We cultivated this association (isolated from the hypersaline alkaline Lake Hotontyn Nur, Mongolia) and characterized their biology. 'Ca. Absconditicoccus praedator' is the first stably cultivated species from the candidate class-level lineage Gracilibacteria (order-level lineage Absconditabacterales). Its closed-and-curated genome lacks genes for the glycolytic, pentose phosphate- and Entner-Doudoroff pathways which would generate energy/reducing equivalents and produce central carbon currencies. Therefore, 'Ca. Absconditicoccus praedator' is dependent on host-derived building blocks for nucleic acid-, protein-, and peptidoglycan synthesis. It shares traits with (the uncultured) 'Ca. Vampirococcus lugosii', which is also of the Gracilibacteria lineage. These are obligate parasitic lifestyle, feeding on photosynthetic anoxygenic Gammaproteobacteria, and absorption of host cytoplasm. Commonalities in their genomic composition and structure suggest that the entire Absconditabacterales lineage consists of predatory species which act to cull the populations of their respective host bacteria. Cultivation of vampire : host associations can shed light on unresolved aspects of their metabolism and ecosystem dynamics at life-limiting extremes.

Patterns of Gene Content and Co-occurrence Constrain the Evolutionary Path toward Animal Association in Candidate Phyla Radiation Bacteria

Candidate Phyla Radiation (CPR) bacteria are small, likely episymbiotic organisms found across Earth's ecosystems. Despite their prevalence, the distribution of CPR lineages across habitats and the genomic signatures of transitions among these habitats remain unclear. Here, we expand the genome inventory for Absconditabacteria (SR1), Gracilibacteria, and Saccharibacteria (TM7), CPR bacteria known to occur in both animal-associated and environmental microbiomes, and investigate variation in gene content with habitat of origin. By overlaying phylogeny with habitat information, we show that bacteria from these three lineages have undergone multiple transitions from environmental habitats into animal microbiomes. Based on co-occurrence analyses of hundreds of metagenomes, we extend the prior suggestion that certain Saccharibacteria have broad bacterial host ranges and constrain possible host relationships for Absconditabacteria and Gracilibacteria. Full-proteome analyses show that animal-associated Saccharibacteria have smaller gene repertoires than their environmental counterparts and are enriched in numerous protein families, including those likely functioning in amino acid metabolism, phage defense, and detoxification of peroxide. In contrast, some freshwater Saccharibacteria encode a putative rhodopsin. For protein families exhibiting the clearest patterns of differential habitat distribution, we compared protein and species phylogenies to estimate the incidence of lateral gene transfer and genomic loss occurring over the species tree. These analyses suggest that habitat transitions were likely not accompanied by large transfer or loss events but rather were associated with continuous proteome remodeling. Thus, we speculate that CPR habitat transitions were driven largely by availability of suitable host taxa and were reinforced by acquisition and loss of some capacities. IMPORTANCE Studying the genetic differences between related microorganisms from different environment types can indicate factors associated with their movement among habitats. This is particularly interesting for bacteria from the Candidate Phyla Radiation because their minimal metabolic capabilities require associations with microbial hosts. We found that shifts of Absconditabacteria, Gracilibacteria, and Saccharibacteria between environmental ecosystems and mammalian mouths/guts probably did not involve major episodes of gene gain and loss; rather, gradual genomic change likely followed habitat migration. The results inform our understanding of how little-known microorganisms establish in the human microbiota where they may ultimately impact health.

Deciphering Symbiotic Interactions of "Candidatus Aenigmarchaeota" with Inferred Horizontal Gene Transfers and Co-occurrence Networks

"Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota") represents one of the earliest proposed evolutionary branches within the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, their ecological roles and potential host-symbiont interactions are still poorly understood. Here, eight metagenome-assembled genomes (MAGs) were reconstructed from hot spring ecosystems, and further in-depth comparative and evolutionary genomic analyses were conducted on these MAGs and other genomes downloaded from public databases. Although with limited metabolic capacities, we reported that "Ca. Aenigmarchaeota" in thermal environments harbor more genes related to carbohydrate metabolism than "Ca. Aenigmarchaeota" in nonthermal environments. Evolutionary analyses suggested that members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota contribute substantially to the niche expansion of "Ca. Aenigmarchaeota" via horizontal gene transfer (HGT), especially genes related to virus defense and stress responses. Based on co-occurrence network results and recent genetic exchanges among community members, we conjectured that "Ca. Aenigmarchaeota" may be symbionts associated with one MAG affiliated with the genus Pyrobaculum, though host specificity might be wide and variable across different "Ca. Aenigmarchaeota" organisms. This study provides significant insight into possible DPANN-host interactions and ecological roles of "Ca. Aenigmarchaeota." IMPORTANCE Recent advances in sequencing technology promoted the blowout discovery of super tiny microbes in the Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota (DPANN) superphylum. However, the unculturable properties of the majority of microbes impeded our investigation of their behavior and symbiotic lifestyle in the corresponding community. By integrating horizontal gene transfer (HGT) detection and co-occurrence network analysis on "Candidatus Aenigmarchaeota" ("Ca. Aenigmarchaeota"), we made one of the first attempts to infer their putative interaction partners and further decipher the potential functional and genetic interactions between the symbionts. We revealed that HGTs contributed by members from the Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota (TACK) superphylum and Euryarchaeota conferred "Ca. Aenigmarchaeota" with the ability to survive under different environmental stresses, such as virus infection, high temperature, and oxidative stress. This study demonstrates that the interaction partners might be inferable by applying informatics analyses on metagenomic sequencing data.

Untapped Potential of Moving Bed Biofilm Reactors with Different Biocarrier Types for Bilge Water Treatment: A Laboratory-Scale Study

Two labscale aerobic moving bed biofilm reactor (MBBR) systems, with a different type of biocarrier in each (K3 and Mutag BioChip), were operated in parallel for the treatment of real saline bilge water. During the operation, different stress conditions were applied in order to evaluate the performance of the systems: organic/hydraulic load shock (chemical oxygen demand (COD): 9 g L−1; hydraulic retention time (HRT): 48–72 h) and salinity shock (salinity: 40 ppt). At the same time, the microbiome in the biofilm and suspended biomass was monitored through 16S rRNA gene analysis in order to describe the changes in the microbial community. The dominant classes were Alphaproteobacteria (families Rhodospirillaceae and Rhodobacteraceae) and Bacteroidia (family Lentimicrobiaceae), being recorded at high relative abundance in all MBBRs. The structure of the biofilm was examined and visualized with scanning electron microscopy (SEM) analysis. Both systems exhibited competent performance, reaching up to 86% removal of COD under high organic loading conditions (COD: 9 g L−1). In the system in which K3 biocarriers were used, the attached and suspended biomass demonstrated a similar trend regarding the changes observed in the microbial communities. In the bioreactor filled with K3 biocarriers, higher concentration of biomass was observed. Biofilm developed on Mutag BioChip biocarriers presented lower biodiversity, while the few species identified in the raw wastewater were not dominant in the bioreactors. Through energy-dispersive X-ray (EDX) analysis of the biofilm, the presence of calcium carbonate was discovered, indicating that biomineralization occurred.

Procaryotic Diversity and Hydrogenotrophic Methanogenesis in an Alkaline Spring (La Crouen, New Caledonia)

(1) Background: The geothermal spring of La Crouen (New Caledonia) discharges warm (42 °C) alkaline water (pH~9) enriched in dissolved nitrogen with traces of methane, but its microbial diversity has not yet been studied. (2) Methods: Cultivation-dependent and -independent methods (e.g., Illumina sequencing and quantitative PCR based on 16S rRNA gene) were used to describe the prokaryotic diversity of this spring. (3) Results: Prokaryotes were mainly represented by Proteobacteria (57% on average), followed by Cyanobacteria, Chlorofexi, and Candidatus Gracilibacteria (GN02/BD1-5) (each > 5%). Both potential aerobes and anaerobes, as well as mesophilic and thermophilic microorganisms, were identified. Some of them had previously been detected in continental hyperalkaline springs found in serpentinizing environments (The Cedars, Samail, Voltri, and Zambales ophiolites). Gammaproteobacteria, Ca. Gracilibacteria and Thermotogae were significantly more abundant in spring water than in sediments. Potential chemolithotrophs mainly included beta- and gammaproteobacterial genera of sulfate-reducers (Ca. Desulfobacillus), methylotrophs (Methyloversatilis), sulfur-oxidizers (Thiofaba, Thiovirga), or hydrogen-oxidizers (Hydrogenophaga). Methanogens (Methanobacteriales and Methanosarcinales) were the dominant Archaea, as found in serpentinization-driven and deep subsurface ecosystems. A novel alkaliphilic hydrogenotrophic methanogen (strain CAN) belonging to the genus Methanobacterium was isolated, suggesting that hydrogenotrophic methanogenesis occurs at La Crouen.

Niche partitioning in the Rimicaris exoculata holobiont: the case of the first symbiotic Zetaproteobacteria

BACKGROUND

Free-living and symbiotic chemosynthetic microbial communities support primary production and higher trophic levels in deep-sea hydrothermal vents. The shrimp Rimicaris exoculata, which dominates animal communities along the Mid-Atlantic Ridge, houses a complex bacterial community in its enlarged cephalothorax. The dominant bacteria present are from the taxonomic groups Campylobacteria, Desulfobulbia (formerly Deltaproteobacteria), Alphaproteobacteria, Gammaproteobacteria, and some recently discovered iron oxyhydroxide-coated Zetaproteobacteria. This epibiotic consortium uses iron, sulfide, methane, and hydrogen as energy sources. Here, we generated shotgun metagenomes from Rimicaris exoculata cephalothoracic epibiotic communities to reconstruct and investigate symbiotic genomes. We collected specimens from three geochemically contrasted vent fields, TAG, Rainbow, and Snake Pit, to unravel the specificity, variability, and adaptation of Rimicaris-microbe associations.

RESULTS

Our data enabled us to reconstruct 49 metagenome-assembled genomes (MAGs) from the TAG and Rainbow vent fields, including 16 with more than 90% completion and less than 5% contamination based on single copy core genes. These MAGs belonged to the dominant Campylobacteria, Desulfobulbia, Thiotrichaceae, and some novel candidate phyla radiation (CPR) lineages. In addition, most importantly, two MAGs in our collection were affiliated to Zetaproteobacteria and had no close relatives (average nucleotide identity ANI < 77% with the closest relative Ghiorsea bivora isolated from TAG, and 88% with each other), suggesting potential novel species. Genes for Calvin-Benson Bassham (CBB) carbon fixation, iron, and sulfur oxidation, as well as nitrate reduction, occurred in both MAGs. However, genes for hydrogen oxidation and multicopper oxidases occurred in one MAG only, suggesting shared and specific potential functions for these two novel Zetaproteobacteria symbiotic lineages. Overall, we observed highly similar symbionts co-existing in a single shrimp at both the basaltic TAG and ultramafic Rainbow vent sites. Nevertheless, further examination of the seeming functional redundancy among these epibionts revealed important differences.

CONCLUSION

These data highlight microniche partitioning in the Rimicaris holobiont and support recent studies showing that functional diversity enables multiple symbiont strains to coexist in animals colonizing hydrothermal vents. Video Abstract.

Schizasterid Heart Urchins Host Microorganisms in a Digestive Symbiosis of Mesozoic Origin

Because of their lifestyles, abundance, and feeding habits, infaunal marine deposit feeders have a significant impact on the ocean floor. As these animals also ingest microorganisms associated with their sediment and seawater diet, their digestive tract usually contains a diverse array of bacteria. However, while most of these microorganisms are transients, some may become part of a resident gut microbiome, in particular when sheltered from the main flow of digesta in specialized gut compartments. Here, we provide an in-depth analysis of the structure and contents of the intestinal caecum (IC), a hindgut diverticulum found exclusively in schizasterid heart urchins (Echinoidea: Spatangoida: Schizasteridae). Based on specimens of Brisaster townsendi, in addition to various other schizasterid taxa, our structural characterization of the IC shows that the organ is a highly specialized gut compartment with unique structural properties. Next generation sequencing shows that the IC contains a microbial population composed predominantly of Bacteroidales, Desulfobacterales, and Spirochaetales. The microbiome of this gut compartment is significantly different in composition and lower in diversity than the microbial population in the sediment-filled main digestive tract. Inferences on the function and evolution of the IC and its microbiome suggest that this symbiosis plays a distinct role in host nutrition and that it evolved at least 66 million years ago during the final phase of the Mesozoic.

The rise of diversity in metabolic platforms across the Candidate Phyla Radiation

BACKGROUND

A unifying feature of the bacterial Candidate Phyla Radiation (CPR) is a limited and highly variable repertoire of biosynthetic capabilities. However, the distribution of metabolic traits across the CPR and the evolutionary processes underlying them are incompletely resolved.

RESULTS

Here, we selected ~ 1000 genomes of CPR bacteria from diverse environments to construct a robust internal phylogeny that was consistent across two unlinked marker sets. Mapping of glycolysis, the pentose phosphate pathway, and pyruvate metabolism onto the tree showed that some components of these pathways are sparsely distributed and that similarity between metabolic platforms is only partially predicted by phylogenetic relationships. To evaluate the extent to which gene loss and lateral gene transfer have shaped trait distribution, we analyzed the patchiness of gene presence in a phylogenetic context, examined the phylogenetic depth of clades with shared traits, and compared the reference tree topology with those of specific metabolic proteins. While the central glycolytic pathway in CPR is widely conserved and has likely been shaped primarily by vertical transmission, there is evidence for both gene loss and transfer especially in steps that convert glucose into fructose 1,6-bisphosphate and glycerate 3P into pyruvate. Additionally, the distribution of Group 3 and Group 4-related NiFe hydrogenases is patchy and suggests multiple events of ancient gene transfer.

CONCLUSIONS

We infer that patterns of gene gain and loss in CPR, including acquisition of accessory traits in independent transfer events, could have been driven by shifts in host-derived resources and led to sparse but varied genetic inventories.

Role of the gut microbiome in mediating standard metabolic rate after dietary shifts in the viviparous cockroach, Diploptera punctata

Diet may be a significant determinant of insect gut microbiome composition. However, the extent to which dietary shifts shape both the composition and relevant functions of insect gut microbiomes, and ultimately impact host energy balance (i.e. metabolic phenotype), is not well understood. We investigated the impacts of diet switching on Diploptera punctata females maintained on a dog food (DF) diet relative to those fed a comparatively sub-optimal cellulose-amended dog food (CADF) diet for 4 weeks. After this period, dietary shift resulted in a significantly higher average mass-specific standard metabolic rate (SMR) in CADF-fed females compared with DF-fed females. We also uncovered significant 13C-enrichment in DF-fed insect samples relative to CADF-fed insect samples and lowered bacterial essential amino acid (EAA) provisioning in CADF-fed samples. Differences in SMR and EAA provisioning were not accompanied by significant differences in overall microbiome composition between the two groups. However, cellulolytic and nitrogen-fixing bacterial families dominant in wild omnivorous cockroaches and wood-feeding termites were significantly enriched in CADF-fed females than in DF-fed females, at the end of the study. We propose that these changes in microbiome composition after dietary shifts are associated with changes in EAA provisioning and possibly SMR. Further studies are needed to comprehensively understand the relative importance of gut microbial functions among the complexity of factors known to underscore SMR responses in insects under varying dietary conditions.

Accurate and complete genomes from metagenomes

Genomes are an integral component of the biological information about an organism; thus, the more complete the genome, the more informative it is. Historically, bacterial and archaeal genomes were reconstructed from pure (monoclonal) cultures, and the first reported sequences were manually curated to completion. However, the bottleneck imposed by the requirement for isolates precluded genomic insights for the vast majority of microbial life. Shotgun sequencing of microbial communities, referred to initially as community genomics and subsequently as genome-resolved metagenomics, can circumvent this limitation by obtaining metagenome-assembled genomes (MAGs); but gaps, local assembly errors, chimeras, and contamination by fragments from other genomes limit the value of these genomes. Here, we discuss genome curation to improve and, in some cases, achieve complete (circularized, no gaps) MAGs (CMAGs). To date, few CMAGs have been generated, although notably some are from very complex systems such as soil and sediment. Through analysis of about 7000 published complete bacterial isolate genomes, we verify the value of cumulative GC skew in combination with other metrics to establish bacterial genome sequence accuracy. The analysis of cumulative GC skew identified potential misassemblies in some reference genomes of isolated bacteria and the repeat sequences that likely gave rise to them. We discuss methods that could be implemented in bioinformatic approaches for curation to ensure that metabolic and evolutionary analyses can be based on very high-quality genomes.

Complete Genome of a Member of a New Bacterial Lineage in the Microgenomates Group Reveals an Unusual Nucleotide Composition Disparity Between Two Strands of DNA and Limited Metabolic Potential

The candidate phyla radiation is a large monophyletic lineage comprising unculturable bacterial taxa with small cell and genome sizes, mostly known from genomes obtained from environmental sources without cultivation. Here, we present the closed complete genome of a member of the superphylum Microgenomates obtained from the metagenome of a deep subsurface thermal aquifer. Phylogenetic analysis indicates that the new bacterium, designated Ch65, represents a novel phylum-level lineage within the Microgenomates group, sibling to the candidate phylum Collierbacteria. The Ch65 genome has a highly unusual nucleotide composition with one strand of highly enriched in cytosine versus guanine throughout the whole length. Such nucleotide composition asymmetry, also detected in the members of Ca. Collierbacteria and Ca. Beckwithbacteria, suggests that most of the Ch65 chromosome is replicated in one direction. A genome analysis predicted that the Ch65 bacterium has fermentative metabolism and could produce acetate and lactate. It lacks respiratory capacity, as well as complete pathways for the biosynthesis of lipids, amino acids, and nucleotides. The Embden–Meyerhof glycolytic pathway and nonoxidative pentose phosphate pathway are mostly complete, although glucokinase, 6-phosphofructokinase, and transaldolase were not found. The Ch65 bacterium lacks secreted glycoside hydrolases and conventional transporters for importing sugars and amino acids. Overall, the metabolic predictions imply that Ch65 adopts the lifestyle of a symbiont/parasite, or a scavenger, obtaining resources from the lysed microbial biomass. We propose the provisional taxonomic assignment ‘Candidatus Chazhemtobacterium aquaticus’, genus ‘Chazhemtobacterium‘, family ‘Chazhemtobacteraceae‘ in the Microgenomates group.

Characterization of the Rumen Microbiota and Volatile Fatty Acid Profiles of Weaned Goat Kids under Shrub-Grassland Grazing and Indoor Feeding

In this study, we conducted comparative analyses to characterize the rumen microbiota and volatile fatty acid (VFA) profiles of weaned Nanjiang Yellow goat kids under shrub-grassland grazing (GR), shrub-grassland grazing and supplementary feeding (SF), and indoor feeding (IF) systems. We observed significant differences (p < 0.05) in the concentrations of total VFA and the proportions of acetate and butyrate in the rumen fluid among the three groups, whereas the proportions of propionate and the acetate/propionate ratio did not differ substantially. Alpha diversity of the rumen bacterial and archaeal populations in the GR and SF kids was significantly higher (p < 0.05) than that in the IF goat kids, and significant differences (p < 0.05) in similarity were observed in the comparisons of GR vs. IF and SF vs. IF. The most predominant bacterial phyla were Bacteroidetes and Firmicutes across the three groups, and the archaeal community was mainly composed of Euryarchaeota. At the genus and species levels, the cellulose-degrading bacteria, including Lachnospiraceae, Ruminococcaceae and Butyrivibrio fibrisolvens, were abundant in the GR and SF groups. Furthermore, 27 bacterial and 11 unique archaeal taxa, such as Lachnospiraceae, Butyrivibrio fibrisolvens, and Methanobrevibacter ruminantium, were identified as biomarkers, and showed significantly different (p < 0.05) abundances among the three groups. Significant Spearman correlations (p < 0.05), between the abundances of several microbial biomarkers and the concentrations of VFAs, were further observed. In summary, our results demonstrated that the adaptation to grazing required more rumen bacterial populations due to complex forage types in shrub-grassland, although the rumen fermentation pattern did not change substantially among the three feeding systems. Some microbial taxa could be used as biomarkers for different feeding systems, particularly cellulose-degrading bacteria associated with grazing.