Anvi'o: an advanced analysis and visualization platform for 'omics data

Seasonal shifts in Fe-acquisition strategies in Southern Ocean microbial communities revealed by metagenomics and autonomous sampling

Iron (Fe) governs the cycling of organic carbon in large parts of the Southern Ocean. The strategies of diverse microbes to acquire the different chemical forms of Fe under seasonally changing organic carbon regimes remain, however, poorly understood. Here, we report high-resolution seasonal metagenomic observations from the region off Kerguelen Island (Indian Sector of the Southern Ocean) where natural Fe-fertilization induces consecutive spring and summer phytoplankton blooms. Our data illustrate pronounced, but distinct seasonal patterns in the abundance of genes implicated in the transport of different forms of Fe and organic substrates, of siderophore biosynthesis and carbohydrate-active enzymes. The seasonal dynamics suggest a temporal decoupling in the prokaryotic requirements of Fe and organic carbon during the spring phytoplankton bloom and a concerted access to these resources after the summer bloom. Taxonomic assignments revealed differences in the prokaryotic groups harbouring genes of a given Fe-related category and pronounced seasonal successions were observed. Using MAGs we could decipher the respective Fe- and organic substrate-related genes of individual taxa assigned to abundant groups. The ecological strategies related to Fe-acquisition provide insights on how this element could shape microbial community composition with potential implications on organic matter transformations in the Southern Ocean.

Novel pelagiphage isolate Polarivirus skadi is a polar specialist that dominates SAR11-associated bacteriophage communities at high latitudes

The SAR11 clade are the most abundant members of surface marine bacterioplankton and a critical component of global biogeochemical cycles. Similarly, pelagiphages that infect SAR11 are ubiquitous and highly abundant in the oceans. Pelagiphages are predicted to shape SAR11 community structures and increase carbon turnover throughout the oceans. Yet, ecological drivers of host and niche specificity of pelagiphage populations are poorly understood. Here we report the global distribution of a novel pelagiphage called "Polarivirus skadi", which is the sole representative of a novel genus. P. skadi was isolated from the Western English Channel using a cold-water ecotype of SAR11 as bait. P. skadi is closely related to the globally dominant pelagiphage HTVC010P. Along with other HTVC010P-type viruses, P. skadi belongs to a distinct viral family within the order Caudovirales, for which we propose the name Ubiqueviridae. Metagenomic read recruitment identified P. skadi as one of the most abundant pelagiphages on Earth. P. skadi is a polar specialist, replacing HTVC010P at high latitudes. Experimental evaluation of P. skadi host range against cold- and warm-water SAR11 ecotypes supported cold-water specialism. Relative abundance of P. skadi in marine metagenomes correlated negatively with temperature, and positively with nutrients, available oxygen, and chlorophyll concentrations. In contrast, relative abundance of HTVC010P correlated negatively with oxygen and positively with salinity, with no significant correlation to temperature. The majority of other pelagiphages were scarce in most marine provinces, with a few representatives constrained to discrete ecological niches. Our results suggest that pelagiphage populations persist within a global viral seed bank, with environmental parameters and host availability selecting for a few ecotypes that dominate ocean viromes.

happi: a hierarchical approach to pangenomics inference

Recovering metagenome-assembled genomes (MAGs) from shotgun sequencing data is an increasingly common task in microbiome studies, as MAGs provide deeper insight into the functional potential of both culturable and non-culturable microorganisms. However, metagenome-assembled genomes vary in quality and may contain omissions and contamination. These errors present challenges for detecting genes and comparing gene enrichment across sample types. To address this, we propose happi, an approach to testing hypotheses about gene enrichment that accounts for genome quality. We illustrate the advantages of happi over existing approaches using published Saccharibacteria MAGs, Streptococcus thermophilus MAGs, and via simulation.

Potential for the anaerobic oxidation of benzene and naphthalene in thermophilic microorganisms from the Guaymas Basin

Unsubstituted aromatic hydrocarbons (UAHs) are recalcitrant molecules abundant in crude oil, which is accumulated in subsurface reservoirs and occasionally enters the marine environment through natural seepage or human-caused spillage. The challenging anaerobic degradation of UAHs by microorganisms, in particular under thermophilic conditions, is poorly understood. Here, we established benzene- and naphthalene-degrading cultures under sulfate-reducing conditions at 50°C and 70°C from Guaymas Basin sediments. We investigated the microorganisms in the enrichment cultures and their potential for UAH oxidation through short-read metagenome sequencing and analysis. Dependent on the combination of UAH and temperature, different microorganisms became enriched. A Thermoplasmatota archaeon was abundant in the benzene-degrading culture at 50°C, but catabolic pathways remained elusive, because the archaeon lacked most known genes for benzene degradation. Two novel species of Desulfatiglandales bacteria were strongly enriched in the benzene-degrading culture at 70°C and in the naphthalene-degrading culture at 50°C. Both bacteria encode almost complete pathways for UAH degradation and for downstream degradation. They likely activate benzene via methylation, and naphthalene via direct carboxylation, respectively. The two species constitute the first thermophilic UAH degraders of the Desulfatiglandales. In the naphthalene-degrading culture incubated at 70°C, a Dehalococcoidia bacterium became enriched, which encoded a partial pathway for UAH degradation. Comparison of enriched bacteria with related genomes from environmental samples indicated that pathways for benzene degradation are widely distributed, while thermophily and capacity for naphthalene activation are rare. Our study highlights the capacities of uncultured thermophilic microbes for UAH degradation in petroleum reservoirs and in contaminated environments.

Proposal to transfer Bacillus massiliigorillae to the genus Peribacillus as Peribacillus massiliigorillae comb. nov., and Bacillus sinesaloumensis to the genus Ferdinandcohnia as Ferdinandcohnia sinesaloumensis comb. nov

The present study was carried out to clarify the taxonomic position of Bacillus massiliigorillae and Bacillus sinesaloumensis. The 16S rRNA gene sequences extracted from the Bacillus sinesaloumensis Marseille-P3516T (FTOX00000000) and Bacillus massiliigorillae G2T (CAVL000000000) genomes showed 98.5 and 99.1% similarity with the type strains of Ferdinandcohnia humi and Peribacillus endoradicis, respectively. The amino acid identity (AAI) values of Bacillus sinesaloumensis Marseille-P3516T were higher with Ferdinandcohnia members, while Bacillus massiliigorillae G2T with Peribacillus members. In phylogenomic and phylogenetic trees, Bacillus sinesaloumensis Marseille-P3516T and Bacillus massiliigorillae G2T clade with members of the genera Ferdinandcohnia and Peribacillus, respectively. Based on the above results, we propose to transfer Bacillus massiliigorillae to the genus Peribacillus as Peribacillus massiliigorillae comb. nov., and Bacillus sinesaloumensis to the genus Ferdinandcohnia as Ferdinandcohnia sinesaloumensis comb. nov.

A glimpse of the paleome in endolithic microbial communities

BACKGROUND

The terrestrial subsurface is home to a significant proportion of the Earth's microbial biomass. Our understanding about terrestrial subsurface microbiomes is almost exclusively derived from groundwater and porous sediments mainly by using 16S rRNA gene surveys. To obtain more insights about biomass of consolidated rocks and the metabolic status of endolithic microbiomes, we investigated interbedded limestone and mudstone from the vadose zone, fractured aquifers, and deep aquitards.

RESULTS

By adapting methods from microbial archaeology and paleogenomics, we could recover sufficient DNA for downstream metagenomic analysis from seven rock specimens independent of porosity, lithology, and depth. Based on the extracted DNA, we estimated between 2.81 and 4.25 × 105 cells × g-1 rock. Analyzing DNA damage patterns revealed paleome signatures (genetic records of past microbial communities) for three rock specimens, all obtained from the vadose zone. DNA obtained from deep aquitards isolated from surface input was not affected by DNA decay indicating that water saturation and not flow is controlling subsurface microbial survival. Decoding the taxonomy and functional potential of paleome communities revealed increased abundances for sequences affiliated with chemolithoautotrophs and taxa such as Cand. Rokubacteria. We also found a broader metabolic potential in terms of aromatic hydrocarbon breakdown, suggesting a preferred utilization of sedimentary organic matter in the past.

CONCLUSIONS

Our study suggests that limestones function as archives for genetic records of past microbial communities including those sensitive to environmental stress at modern times, due to their specific conditions facilitating long-term DNA preservation. Video Abstract.

Treatment performance and microbial community structure in an aerobic granular sludge sequencing batch reactor amended with diclofenac, erythromycin, and gemfibrozil

This study characterizes the effects of three commonly detected pharmaceuticals-diclofenac, erythromycin, and gemfibrozil-on aerobic granular sludge. Approximately 150 μg/L of each pharmaceutical was fed in the influent to a sequencing batch reactor for 80 days, and the performance of the test reactor was compared with that of a control reactor. Wastewater treatment efficacy in the test reactor dropped by approximately 30-40%, and ammonia oxidation was particularly inhibited. The relative abundance of active Rhodocyclaceae, Nitrosomonadaceae, and Nitrospiraceae families declined throughout exposure, likely explaining reductions in wastewater treatment performance. Pharmaceuticals were temporarily removed in the first 12 days of the test via both sorption and degradation; both removal processes declined sharply thereafter. This study demonstrates that aerobic granular sludge may successfully remove pharmaceuticals in the short term, but long-term tests are necessary to confirm if pharmaceutical removal is sustainable.

ACR: metagenome-assembled prokaryotic and eukaryotic genome refinement tool

Microbial genome recovery from metagenomes can further explain microbial ecosystem structures, functions and dynamics. Thus, this study developed the Additional Clustering Refiner (ACR) to enhance high-purity prokaryotic and eukaryotic metagenome-assembled genome (MAGs) recovery. ACR refines low-quality MAGs by subjecting them to iterative k-means clustering predicated on contig abundance and increasing bin purity through validated universal marker genes. Synthetic and real-world metagenomic datasets, including short- and long-read sequences, evaluated ACR's effectiveness. The results demonstrated improved MAG purity and a significant increase in high- and medium-quality MAG recovery rates. In addition, ACR seamlessly integrates with various binning algorithms, augmenting their strengths without modifying core features. Furthermore, its multiple sequencing technology compatibilities expand its applicability. By efficiently recovering high-quality prokaryotic and eukaryotic genomes, ACR is a promising tool for deepening our understanding of microbial communities through genome-centric metagenomics.

Differences in fine arabinoxylan structures govern microbial selection and competition among human gut microbiota

Dietary fibers are known to modulate microbiome composition, but it is unclear to what extent minor fiber structural differences impact community assembly, microbial division of labor, and organismal metabolic responses. To test the hypothesis that fine linkage variations afford different ecological niches for distinct communities and metabolism, we employed a 7-day in vitro sequential batch fecal fermentation with four fecal inocula and measured responses using an integrated multi-omics approach. Two sorghum arabinoxylans (SAXs) were fermented, with one (RSAX) having slightly more complex branch linkages than the other (WSAX). Although there were minor glycoysl linkage differences, consortia on RSAX retained much higher species diversity (42 members) than on WSAX (18-23 members) with distinct species-level genomes and metabolic outcomes (e.g., higher short chain fatty acid production from RSAX and more lactic acid produced from WSAX). The major SAX-selected members were from genera of Bacteroides and Bifidobacterium and family Lachnospiraceae. Carbohydrate active enzyme (CAZyme) genes in metagenomes revealed broad AX-related hydrolytic potentials among key members; however, CAZyme genes enriched in different consortia displayed various catabolic domain fusions with diverse accessory motifs that differ among the two SAX types. These results suggest that fine polysaccharide structure exerts deterministic selection effect for distinct fermenting consortia.

Genome Analysis and Potential Ecological Functions of Members of the Genus Ensifer from Subsurface Environments and Description of Ensifer oleiphilus sp. nov

The current work deals with genomic analysis, possible ecological functions, and biotechnological potential of two bacterial strains, HO-A22T and SHC 2-14, isolated from unique subsurface environments, the Cheremukhovskoe oil field (Tatarstan, Russia) and nitrate- and radionuclide-contaminated groundwater (Tomsk region, Russia), respectively. New isolates were characterized using polyphasic taxonomy approaches and genomic analysis. The genomes of the strains HO-A22T and SHC 2-14 contain the genes involved in nitrate reduction, hydrocarbon degradation, extracellular polysaccharide synthesis, and heavy metal detoxification, confirming the potential for their application in various environmental biotechnologies. Genomic data were confirmed by cultivation studies. Both strains were found to be neutrophilic, chemoorganotrophic, facultatively anaerobic bacteria, growing at 15-33 °C and 0-1.6% NaCl (w/v). The 16S rRNA gene sequences of the strains were similar to those of the type strains of the genus Ensifer (99.0-100.0%). Nevertheless, genomic characteristics of strain HO-A22T were below the thresholds for species delineation: the calculated average nucleotide identity (ANI) values were 83.7-92.4% (<95%), and digital DNA-DNA hybridization (dDDH) values were within the range of 25.4-45.9% (<70%), which supported our conclusion that HO-A22T (=VKM B-3646T = KCTC 92427T) represented a novel species of the genus Ensifer, with the proposed name Ensifer oleiphilus sp. nov. Strain SHC 2-14 was assigned to the species 'Ensifer canadensis', which has not been validly published. This study expanded the knowledge about the phenotypic diversity among members of the genus Ensifer and its potential for the biotechnologies of oil recovery and radionuclide pollution treatment.

Sulfur-cycling chemolithoautotrophic microbial community dominates a cold, anoxic, hypersaline Arctic spring

BACKGROUND

Gypsum Hill Spring, located in Nunavut in the Canadian High Arctic, is a rare example of a cold saline spring arising through thick permafrost. It perennially discharges cold (~ 7 °C), hypersaline (7-8% salinity), anoxic (~ 0.04 ppm O2), and highly reducing (~ - 430 mV) brines rich in sulfate (2.2 g.L-1) and sulfide (9.5 ppm), making Gypsum Hill an analog to putative sulfate-rich briny habitats on extraterrestrial bodies such as Mars.

RESULTS

Genome-resolved metagenomics and metatranscriptomics were utilized to describe an active microbial community containing novel metagenome-assembled genomes and dominated by sulfur-cycling Desulfobacterota and Gammaproteobacteria. Sulfate reduction was dominated by hydrogen-oxidizing chemolithoautotrophic Desulfovibrionaceae sp. and was identified in phyla not typically associated with sulfate reduction in novel lineages of Spirochaetota and Bacteroidota. Highly abundant and active sulfur-reducing Desulfuromusa sp. highly transcribed non-coding RNAs associated with transcriptional regulation, showing potential evidence of putative metabolic flexibility in response to substrate availability. Despite low oxygen availability, sulfide oxidation was primarily attributed to aerobic chemolithoautotrophic Halothiobacillaceae. Low abundance and transcription of photoautotrophs indicated sulfur-based chemolithoautotrophy drives primary productivity even during periods of constant illumination.

CONCLUSIONS

We identified a rare surficial chemolithoautotrophic, sulfur-cycling microbial community active in a unique anoxic, cold, hypersaline Arctic spring. We detected Mars-relevant metabolisms including hydrogenotrophic sulfate reduction, sulfur reduction, and sulfide oxidation, which indicate the potential for microbial life in analogous S-rich brines on past and present Mars. Video Abstract.

Description of Tellurirhabdus bombi sp. nov., Isolated from Bumblebee

A Gram-stain-negative, aerobic, non-motile, and rod-shaped bacterium, designated IE-0392T, was isolated from a bumblebee. The 16S rRNA gene sequence (highest 16S rRNA gene sequence similarity with the type strain of Tellurirhabdus rosea (90.0%) and phylogenetic analysis suggest that strain IE-0392T was a member of the genus Tellurirhabdus. Strain IE-0392T optimally grew at 25 ℃ and pH 7.0. Menaquinone 7 (MK-7) was the only isoprenoid quinone present in strain IE-0392T. The major fatty acids (> 10%) of strain IE-0392T were iso-C15:0, C16:1 ω5c, and iso-C17:0 3-OH. The polar lipids of strain IE-0392T were phosphatidylethanolamine, phosphatidylserine, unidentified aminophospholipids, unidentified aminolipid, unidentified phospholipid, and unidentified lipids. The genomic DNA G + C content of strain IE-0392T was 48.8%. The amino acid identity (AAI) and the average nucleotide identity (ANI) values suggest that strain IE-0392T is a novel member of the genus Tellurirhabdus. The results suggest that strain IE-0392T represents a novel species of the genus Tellurirhabdus, for which the name Tellurirhabdus bombi sp. nov., is proposed. The type strain is IE-0392T (= GDMCC 1.2794T = JCM 35040T).

Genomes and transcriptomes help unravel the complex life cycle of the blastoclad fungus, Coelomomyces lativittatus, an obligate parasite of mosquitoes and microcrustaceans

Species of the phylum Blastocladiomycota, early-diverging zoosporic (flagellated) lineages of fungi, are vastly understudied. This phylum includes the genus Coelomomyces, which consists of more than 80 fungal species that are obligate parasites of arthropods. Known Coelomomyces species lack a complete asexual life cycle, instead surviving through an obligate heterecious alternation of generations life cycle. Despite their global distribution and interesting life cycle, little is known about the genomics of any Coelomomyces species. To address this, we generated three draft-level genomes and annotations for C. lativittatus representing its haploid meiospore, orange gamete, and amber gamete life stages. These draft genome assemblies ranged in size from 5002 to 5799 contigs, with a total length of 19.8-22.8 Mb and a mean of 7416 protein-coding genes. We then demonstrated the utility of these genomes by combining the draft annotations as a reference for analysis of C. lativittatus transcriptomes. We analyzed transcriptomes from across host-associated life stages, including infected larvae and excised mature sporangia from the mosquito Anopheles quadrimaculatus. We identified differentially expressed genes and enriched GO terms both across and within life stages and used these to make hypotheses about C. lativittatus biology. Generally, we found the C. lativittatus transcriptome to be a complex and dynamic expression landscape; GO terms related to metabolism and transport processes were enriched during infection and terms related to dispersal were enriched during sporulation. We further identified five high mobility group (HMG)-box genes in C. lativittatus, three belonging to clades with mating type (MAT) loci from other fungi, as well as four ortholog expansions in C. lativittatus compared with other fungi. The C. lativittatus genomes and transcriptomes reported here are a valuable resource and may be leveraged toward furthering understanding of the biology of these and other early-diverging fungal lineages.

Taxonomic considerations on Aerococcus urinae with proposal of subdivision into Aerococcus urinae, Aerococcus tenax sp. nov., Aerococcus mictus sp. nov., and Aerococcus loyolae sp. nov

Average nucleotide identity analysis, based on whole genome sequences of 115 strains previously identified as Aerococcus urinae, an emerging uropathogen, discriminates at least six unique genomic taxa. The whole genome analysis affords clearer species boundaries over 16S rRNA gene sequencing and traditional phenotypic approaches for the identification and phylogenetic organization of Aerococcus species. The newly described species can be differentiated by matrix-assisted laser desorption ionization time-of-flight analysis of protein signatures. We propose the emendation of the description of A. urinae (type strain ATCC 51268T = CCUG 34223T=NCFB 2893) and the names of Aerococcus tenax sp. nov. (ATCC TSD-302T = DSM 115700T = CCUG 76531T=NR-58630T), Aerococcus mictus sp. nov. (ATCC TSD-301T = DSM 115699T = CCUG 76532T=NR-58629T), and Aerococcus loyolae sp. nov. (ATCC TSD-300T = DSM 115698T = CCUG 76533T=NR-58628T) for three of the newly identified genomic taxa.

Predicting variable gene content in Escherichia coli using conserved genes

Having the ability to predict the protein-encoding gene content of an incomplete genome or metagenome-assembled genome is important for a variety of bioinformatic tasks. In this study, as a proof of concept, we built machine learning classifiers for predicting variable gene content in Escherichia coli genomes using only the nucleotide k-mers from a set of 100 conserved genes as features. Protein families were used to define orthologs, and a single classifier was built for predicting the presence or absence of each protein family occurring in 10%-90% of all E. coli genomes. The resulting set of 3,259 extreme gradient boosting classifiers had a per-genome average macro F1 score of 0.944 [0.943-0.945, 95% CI]. We show that the F1 scores are stable across multi-locus sequence types and that the trend can be recapitulated by sampling a smaller number of core genes or diverse input genomes. Surprisingly, the presence or absence of poorly annotated proteins, including "hypothetical proteins" was accurately predicted (F1 = 0.902 [0.898-0.906, 95% CI]). Models for proteins with horizontal gene transfer-related functions had slightly lower F1 scores but were still accurate (F1s = 0.895, 0.872, 0.824, and 0.841 for transposon, phage, plasmid, and antimicrobial resistance-related functions, respectively). Finally, using a holdout set of 419 diverse E. coli genomes that were isolated from freshwater environmental sources, we observed an average per-genome F1 score of 0.880 [0.876-0.883, 95% CI], demonstrating the extensibility of the models. Overall, this study provides a framework for predicting variable gene content using a limited amount of input sequence data. IMPORTANCE Having the ability to predict the protein-encoding gene content of a genome is important for assessing genome quality, binning genomes from shotgun metagenomic assemblies, and assessing risk due to the presence of antimicrobial resistance and other virulence genes. In this study, we built a set of binary classifiers for predicting the presence or absence of variable genes occurring in 10%-90% of all publicly available E. coli genomes. Overall, the results show that a large portion of the E. coli variable gene content can be predicted with high accuracy, including genes with functions relating to horizontal gene transfer. This study offers a strategy for predicting gene content using limited input sequence data.

Stratified microbial communities in Australia's only anchialine cave are taxonomically novel and drive chemotrophic energy production via coupled nitrogen-sulphur cycling

BACKGROUND

Anchialine environments, in which oceanic water mixes with freshwater in coastal aquifers, are characterised by stratified water columns with complex physicochemical profiles. These environments, also known as subterranean estuaries, support an abundance of endemic macro and microorganisms. There is now growing interest in characterising the metabolisms of anchialine microbial communities, which is essential for understanding how complex ecosystems are supported in extreme environments, and assessing their vulnerability to environmental change. However, the diversity of metabolic strategies that are utilised in anchialine ecosystems remains poorly understood.

RESULTS

Here, we employ shotgun metagenomics to elucidate the key microorganisms and their dominant metabolisms along a physicochemical profile in Bundera Sinkhole, the only known continental subterranean estuary in the Southern Hemisphere. Genome-resolved metagenomics suggests that the communities are largely represented by novel taxonomic lineages, with 75% of metagenome-assembled genomes assigned to entirely new or uncharacterised families. These diverse and novel taxa displayed depth-dependent metabolisms, reflecting distinct phases along dissolved oxygen and salinity gradients. In particular, the communities appear to drive nutrient feedback loops involving nitrification, nitrate ammonification, and sulphate cycling. Genomic analysis of the most highly abundant members in this system suggests that an important source of chemotrophic energy is generated via the metabolic coupling of nitrogen and sulphur cycling.

CONCLUSION

These findings substantially contribute to our understanding of the novel and specialised microbial communities in anchialine ecosystems, and highlight key chemosynthetic pathways that appear to be important in these energy-limited environments. Such knowledge is essential for the conservation of anchialine ecosystems, and sheds light on adaptive processes in extreme environments. Video Abstract.

Complete De Novo Assembly of Wolbachia Endosymbiont of Frankliniella intonsa

As an endosymbiont, Wolbachia exerts significant effects on the host, including on reproduction, immunity, and metabolism. However, the study of Wolbachia in Thysanopteran insects, such as flower thrips Frankliniella intonsa, remains limited. Here, we assembled a gap-free looped genome assembly of Wolbachia strain wFI in a length of 1,463,884 bp (GC content 33.80%), using Nanopore long reads and Illumina short reads. The annotation of wFI identified a total of 1838 protein-coding genes (including 85 pseudogenes), 3 ribosomal RNAs (rRNAs), 35 transfer RNAs (tRNAs), and 1 transfer-messenger RNA (tmRNA). Beyond this basic description, we identified mobile genetic elements, such as prophage and insertion sequences (ISs), which make up 17% of the entire wFI genome, as well as genes involved in riboflavin and biotin synthesis and metabolism. This research lays the foundation for understanding the nutritional mutualism between Wolbachia and flower thrips. It also serves as a valuable resource for future studies delving into the intricate interactions between Wolbachia and its host.

Microbial communities related to the sulfur cycle in the Sansha Yongle Blue Hole

The Sansha Yongle Blue Hole (SYBH), the deepest blue hole in the world, is an excellent habitat for revealing biogeochemical cycles in the anaerobic environment. However, how sulfur cycling is mediated by microorganisms in the SYBH hasn't been fully understood. In this study, the water layers of the SYBH were divided into oxic zone, hypoxic zone, anoxic zone I and II, and microbial-mediated sulfur cycling in the SYBH was comprehensively interpreted. The 16S rRNA genes/transcripts analyses showed that the microbial community structures associated with the sulfur cycling in each zone had distinctive features. Sulfur-oxidizing bacteria were mostly constituted by Gammaproteobacteria, Alphaproteobacteria, Campylobacterota, and Chlorobia above the anoxic zone I and sulfate-reducing bacteria were dominated by Desulfobacterota in anoxic zones. Metagenomic analyses showed that the sulfide-oxidation-related gene sqr and genes encoding the Sox system were mainly distributed in the anoxic zone I, while genes related to dissimilatory sulfate reduction and sulfur intermediate metabolite reduction were mainly distributed in the anoxic zone II, indicating different sulfur metabolic processes between these two zones. Moreover, sulfur-metabolism-related genes were identified in 81 metagenome-assembled genomes (MAGs), indicating a high diversity of microbial communities involved in sulfur cycling. Among them, three MAGs from the candidate phyla JdFR-76 and AABM5-125-24 with genes related to dissimilatory sulfate reduction exhibited distinctive metabolic features. Our results showed unique and novel microbial populations in the SYBH sulfur cycle correlated to the sharp redox gradients, revealing complex biogeochemical processes in this extreme environment. IMPORTANCE Oxygen-deficient regions in the global ocean are expanding rapidly and affect the growth, reproduction and ecological processes of marine organisms. The anaerobic water body of about 150 m in the Sansha Yongle Blue Hole (SYBH) provided a suitable environment to study the specific microbial metabolism in anaerobic seawater. Here, we found that the vertical distributions of the total and active communities of sulfur-oxidizing bacteria (SOB) and sulfate-reducing bacteria (SRB) were different in each water layer of the SYBH according to the dissolved oxygen content. Genes related to sulfur metabolism also showed distinct stratification characteristics. Furthermore, we have obtained diverse metagenome-assembled genomes, some of which exhibit special sulfur metabolic characteristics, especially candidate phyla JdFR-76 and AABM5-125-24 were identified as potential novel SRB. The results of this study will promote further understanding of the sulfur cycle in extreme environments, as well as the environmental adaptability of microorganisms in blue holes.

Exploring the genotypic and phenotypic differences distinguishing Lactobacillus jensenii and Lactobacillus mulieris

Lactobacillus crispatus, Lactobacillus gasseri, Lactobacillus iners, and Lactobacillus jensenii are dominant species of the urogenital microbiota. Prior studies suggest that these Lactobacillus species play a significant role in the urobiome of healthy females. In our prior genomic analysis of all publicly available L. jensenii and Lactobacillus mulieris genomes at the time (n = 43), we identified genes unique to these two closely related species. This motivated our further exploration here into their genotypic differences as well as into their phenotypic differences. First, we expanded genome sequence representatives of both species to 61 strains, including publicly available strains and nine new strains sequenced here. Genomic analyses conducted include phylogenetics of the core genome as well as biosynthetic gene cluster analysis and metabolic pathway analyses. Urinary strains of both species were assayed for their ability to utilize four simple carbohydrates. We found that L. jensenii strains can efficiently catabolize maltose, trehalose, and glucose, but not ribose, and L. mulieris strains can utilize maltose and glucose, but not trehalose and ribose. Metabolic pathway analysis clearly shows the lack of treB within L. mulieris strains, indicative of its inability to catabolize external sources of trehalose. While genotypic and phenotypic observations provide insight into the differences between these two species, we did not find any association with urinary symptom status. Through this genomic and phenotypic investigation, we identify markers that can be leveraged to clearly distinguish these two species in investigations of the female urogenital microbiota. IMPORTANCE We have expanded upon our prior genomic analysis of L. jensenii and L. mulieris strains, including nine new genome sequences. Our bioinformatic analysis finds that L. jensenii and L. mulieris cannot be distinguished by short-read 16S rRNA gene sequencing alone. Thus, to discriminate between these two species, future studies of the female urogenital microbiome should employ metagenomic sequencing and/or sequence species-specific genes, such as those identified here. Our bioinformatic examination also confirmed our prior observations of differences between the two species related to genes associated with carbohydrate utilization, which we tested here. We found that the transport and utilization of trehalose are key distinguishing traits of L. jensenii, which is further supported by our metabolic pathway analysis. In contrast with other urinary Lactobacillus species, we did not find strong evidence for either species, nor particular genotypes, to be associated with lower urinary tract symptoms (or the lack thereof).

Comparative Genomic Analyses of Virulence and Antimicrobial Resistance in Citrobacter werkmanii, an Emerging Opportunistic Pathogen

Citrobacter werkmanii is an emerging and opportunistic human pathogen found in developing countries and is a causative agent of wound, urinary tract, and blood infections. The present study conducted comparative genomic analyses of a C. werkmanii strain collection from diverse geographical locations and sources to identify the relevant virulence and antimicrobial resistance genes. Pangenome analyses divided the examined C. werkmanii strains into five distinct clades; the subsequent classification identified genes with functional roles in carbohydrate and general metabolism for the core genome and genes with a role in secretion, adherence, and the mobilome for the shell and cloud genomes. A maximum-likelihood phylogenetic tree with a heatmap, showing the virulence and antimicrobial genes' presence or absence, demonstrated the presence of genes with functional roles in secretion systems, adherence, enterobactin, and siderophore among the strains belonging to the different clades. C. werkmanii strains in clade V, predominantly from clinical sources, harbored genes implicated in type II and type Vb secretion systems as well as multidrug resistance to aminoglycoside, beta-lactamase, fluoroquinolone, phenicol, trimethoprim, macrolides, sulfonamide, and tetracycline. In summary, these comparative genomic analyses have demonstrated highly pathogenic and multidrug-resistant genetic profiles in C. werkmanii strains, indicating a virulence potential for this commensal and opportunistic human pathogen.

Seasonal patterns in microbial carbon and iron transporter expression in the Southern Ocean

BACKGROUND

Heterotrophic microbes in the Southern Ocean are challenged by the double constraint of low concentrations of organic carbon (C) and iron (Fe). These essential elements are tightly coupled in cellular processes; however, the prokaryotic requirements of C and Fe under varying environmental settings remain poorly studied. Here, we used a combination of metatranscriptomics and metaproteomics to identify prokaryotic membrane transporters for organic substrates and Fe in naturally iron-fertilized and high-nutrient, low-chlorophyll waters of the Southern Ocean during spring and late summer.

RESULTS

Pronounced differences in membrane transporter profiles between seasons were observed at both sites, both at the transcript and protein level. When specific compound classes were considered, the two approaches revealed different patterns. At the transcript level, seasonal patterns were only observed for subsets of genes belonging to each transporter category. At the protein level, membrane transporters of organic compounds were relatively more abundant in spring as compared to summer, while the opposite pattern was observed for Fe transporters. These observations suggest an enhanced requirement for organic C in early spring and for Fe in late summer. Mapping transcripts and proteins to 50 metagenomic-assembled genomes revealed distinct taxon-specific seasonal differences pointing to potentially opportunistic clades, such as Pseudomonadales and Nitrincolaceae, and groups with a more restricted repertoire of expressed transporters, such as Alphaproteobacteria and Flavobacteriaceae.

CONCLUSION

The combined investigations of C and Fe membrane transporters suggest seasonal changes in the microbial requirements of these elements under different productivity regimes. The taxon-specific acquisition strategies of different forms of C and Fe illustrate how diverse microbes could shape transcript and protein expression profiles at the community level at different seasons. Our results on the C- and Fe-related metabolic capabilities of microbial taxa provide new insights into their potential role in the cycling of C and Fe under varying nutrient regimes in the Southern Ocean. Video Abstract.

BinaRena: a dedicated interactive platform for human-guided exploration and binning of metagenomes

BACKGROUND

Exploring metagenomic contigs and "binning" them into metagenome-assembled genomes (MAGs) are essential for the delineation of functional and evolutionary guilds within microbial communities. Despite the advances in automated binning algorithms, their capabilities in recovering MAGs with accuracy and biological relevance are so far limited. Researchers often find that human involvement is necessary to achieve representative binning results. This manual process however is expertise demanding and labor intensive, and it deserves to be supported by software infrastructure.

RESULTS

We present BinaRena, a comprehensive and versatile graphic interface dedicated to aiding human operators to explore metagenome assemblies via customizable visualization and to associate contigs with bins. Contigs are rendered as an interactive scatter plot based on various data types, including sequence metrics, coverage profiles, taxonomic assignments, and functional annotations. Various contig-level operations are permitted, such as selection, masking, highlighting, focusing, and searching. Binning plans can be conveniently edited, inspected, and compared visually or using metrics including silhouette coefficient and adjusted Rand index. Completeness and contamination of user-selected contigs can be calculated in real time. In demonstration of BinaRena's usability, we show that it facilitated biological pattern discovery, hypothesis generation, and bin refinement in a complex tropical peatland metagenome. It enabled isolation of pathogenic genomes within closely related populations from the gut microbiota of diarrheal human subjects. It significantly improved overall binning quality after curating results of automated binners using a simulated marine dataset.

CONCLUSIONS

BinaRena is an installation-free, dependency-free, client-end web application that operates directly in any modern web browser, facilitating ease of deployment and accessibility for researchers of all skill levels. The program is hosted at https://github.com/qiyunlab/binarena , together with documentation, tutorials, example data, and a live demo. It effectively supports human researchers in intuitive interpretation and fine tuning of metagenomic data. Video Abstract.

Sargasso Sea Vibrio bacteria: Underexplored potential pathovars in a perturbed habitat

We fully sequenced the genomes of 16 Vibrio cultivars isolated from eel larvae, plastic marine debris (PMD), the pelagic brown macroalga Sargassum, and seawater samples collected from the Caribbean and Sargasso Seas of the North Atlantic Ocean. Annotation and mapping of these 16 bacterial genome sequences to a PMD-derived Vibrio metagenome-assembled genome created for this study showcased vertebrate pathogen genes closely-related to cholera and non-cholera pathovars. Phenotype testing of cultivars confirmed rapid biofilm formation, hemolytic, and lipophospholytic activities, consistent with pathogenic potential. Our study illustrates that open ocean vibrios represent a heretofore undescribed group of microbes, some representing potential new species, possessing an amalgam of pathogenic and low nutrient acquisition genes, reflecting their pelagic habitat and the substrates and hosts they colonize.

Genomic characterization and virulence of Streptococcus suis serotype 4 clonal complex 94 recovered from human and swine samples

Streptococcus suis is a zoonotic pathogen that causes invasive infections in humans and pigs. Herein, we performed genomic analysis of seven S. suis serotype 4 strains belonging to clonal complex (CC) 94 that were recovered from a human patient or from diseased and clinically healthy pigs. Genomic exploration and comparisons, as well as in vitro cytotoxicity tests, indicated that S. suis CC94 serotype 4 strains are potentially virulent. Genomic analysis revealed that all seven strains clustered within minimum core genome group 3 (MCG-3) and had a high number of virulence-associated genes similar to those of virulent serotype 2 strains. Cytotoxicity assays showed that both the human lung adenocarcinoma cell line and HeLa cells rapidly lost viability following incubation for 4 h with the strains at a concentration of 106 bacterial cells. The human serotype 4 strain (ID36054) decreased cell viability profoundly and similarly to the control serotype 2 strain P1/7. In addition, strain ST1689 (ID34572), isolated from a clinically healthy pig, presented similar behaviour in an adenocarcinoma cell line and HeLa cells. The antimicrobial resistance genes tet(O) and ermB that confer resistance to tetracyclines, macrolides, and lincosamides were commonly found in the strains. However, aminoglycoside and streptothricin resistance genes were found only in certain strains in this study. Our results indicate that S. suis CC94 serotype 4 strains are potentially pathogenic and virulent and should be monitored.

Novel complete methanogenic pathways in longitudinal genomic study of monogastric age-associated archaea

BACKGROUND

Archaea perform critical roles in the microbiome system, including utilizing hydrogen to allow for enhanced microbiome member growth and influencing overall host health. With the majority of microbiome research focusing on bacteria, the functions of archaea are largely still under investigation. Understanding methanogenic functions during the host lifetime will add to the limited knowledge on archaeal influence on gut and host health. In our study, we determined lifelong archaea dynamics, including detection and methanogenic functions, while assessing global, temporal and host distribution of our novel archaeal metagenome-assembled genomes (MAGs). We followed 7 monogastric swine throughout their life, from birth to adult (1-156 days of age), and collected feces at 22 time points. The samples underwent gDNA extraction, Illumina sequencing, bioinformatic quality and assembly processes, MAG taxonomic assignment and functional annotation. MAGs were utilized in downstream phylogenetic analysis for global, temporal and host distribution in addition to methanogenic functional potential determination.

RESULTS

We generated 1130 non-redundant MAGs, representing 588 unique taxa at the species level, with 8 classified as methanogenic archaea. The taxonomic classifications were as follows: orders Methanomassiliicoccales (5) and Methanobacteriales (3); genera UBA71 (3), Methanomethylophilus (1), MX-02 (1), and Methanobrevibacter (3). We recovered the first US swine Methanobrevibacter UBA71 sp006954425 and Methanobrevibacter gottschalkii MAGs. The Methanobacteriales MAGs were identified primarily during the young, preweaned host whereas Methanomassiliicoccales primarily in the adult host. Moreover, we identified our methanogens in metagenomic sequences from Chinese swine, US adult humans, Mexican adult humans, Swedish adult humans, and paleontological humans, indicating that methanogens span different hosts, geography and time. We determined complete metabolic pathways for all three methanogenic pathways: hydrogenotrophic, methylotrophic, and acetoclastic. This study provided the first evidence of acetoclastic methanogenesis in archaea of monogastric hosts which indicated a previously unknown capability for acetate utilization in methanogenesis for monogastric methanogens. Overall, we hypothesized that the age-associated detection patterns were due to differential substrate availability via the host diet and microbial metabolism, and that these methanogenic functions are likely crucial to methanogens across hosts. This study provided a comprehensive, genome-centric investigation of monogastric-associated methanogens which will further improve our understanding of microbiome development and functions.

Insights into the physiological and genomic characterization of three bacterial isolates from a highly alkaline, terrestrial serpentinizing system

The terrestrial serpentinite-hosted ecosystem known as "The Cedars" is home to a diverse microbial community persisting under highly alkaline (pH ~ 12) and reducing (Eh < -550 mV) conditions. This extreme environment presents particular difficulties for microbial life, and efforts to isolate microorganisms from The Cedars over the past decade have remained challenging. Herein, we report the initial physiological assessment and/or full genomic characterization of three isolates: Paenibacillus sp. Cedars ('Paeni-Cedars'), Alishewanella sp. BS5-314 ('Ali-BS5-314'), and Anaerobacillus sp. CMMVII ('Anaero-CMMVII'). Paeni-Cedars is a Gram-positive, rod-shaped, mesophilic facultative anaerobe that grows between pH 7-10 (minimum pH tested was 7), temperatures 20-40°C, and 0-3% NaCl concentration. The addition of 10-20 mM CaCl₂ enhanced growth, and iron reduction was observed in the following order, 2-line ferrihydrite > magnetite > serpentinite ~ chromite ~ hematite. Genome analysis identified genes for flavin-mediated iron reduction and synthesis of a bacillibactin-like, catechol-type siderophore. Ali-BS5-314 is a Gram-negative, rod-shaped, mesophilic, facultative anaerobic alkaliphile that grows between pH 10-12 and temperatures 10-40°C, with limited growth observed 1-5% NaCl. Nitrate is used as a terminal electron acceptor under anaerobic conditions, which was corroborated by genome analysis. The Ali-BS5-314 genome also includes genes for benzoate-like compound metabolism. Anaero-CMMVII remained difficult to cultivate for physiological studies; however, growth was observed between pH 9-12, with the addition of 0.01-1% yeast extract. Anaero-CMMVII is a probable oxygen-tolerant anaerobic alkaliphile with hydrogenotrophic respiration coupled with nitrate reduction, as determined by genome analysis. Based on single-copy genes, ANI, AAI and dDDH analyses, Paeni-Cedars and Ali-BS5-314 are related to other species (P. glucanolyticus and A. aestuarii, respectively), and Anaero-CMMVII represents a new species. The characterization of these three isolates demonstrate the range of ecophysiological adaptations and metabolisms present in serpentinite-hosted ecosystems, including mineral reduction, alkaliphily, and siderophore production.

Ultra-deep sequencing of Hadza hunter-gatherers recovers vanishing gut microbes

The gut microbiome modulates immune and metabolic health. Human microbiome data are biased toward industrialized populations, limiting our understanding of non-industrialized microbiomes. Here, we performed ultra-deep metagenomic sequencing on 351 fecal samples from the Hadza hunter-gatherers of Tanzania and comparative populations in Nepal and California. We recovered 91,662 genomes of bacteria, archaea, bacteriophages, and eukaryotes, 44% of which are absent from existing unified datasets. We identified 124 gut-resident species vanishing in industrialized populations and highlighted distinct aspects of the Hadza gut microbiome related to in situ replication rates, signatures of selection, and strain sharing. Industrialized gut microbes were found to be enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. This unparalleled view of the Hadza gut microbiome provides a valuable resource, expands our understanding of microbes capable of colonizing the human gut, and clarifies the extensive perturbation induced by the industrialized lifestyle.

Metapangenomic investigation provides insight into niche differentiation of methanogenic populations from the subsurface serpentinizing environment, Samail Ophiolite, Oman

Serpentinization reactions produce highly reduced waters that have hyperalkaline pH and that can have high concentrations of H₂ and CH₄. Putatively autotrophic methanogenic archaea have been identified in the subsurface waters of the Samail Ophiolite, Sultanate of Oman, though the strategies to overcome hyperalkaline pH and dissolved inorganic carbon limitation remain to be fully understood. Here, we recovered metagenome assembled genomes (MAGs) and applied a metapangenomic approach to three different Methanobacterium populations to assess habitat-specific functional gene distribution. A Type I population was identified in the fluids with neutral pH, while a Type II and "Mixed" population were identified in the most hyperalkaline fluids (pH 11.63). The core genome of all Methanobacterium populations highlighted potential DNA scavenging techniques to overcome phosphate or nitrogen limitation induced by environmental conditions. With particular emphasis on the Mixed and Type II population found in the most hyperalkaline fluids, the accessory genomes unique to each population reflected adaptation mechanisms suggesting lifestyles that minimize niche overlap. In addition to previously reported metabolic capability to utilize formate as an electron donor and generate intracellular CO₂, the Type II population possessed genes relevant to defense against antimicrobials and assimilating potential osmoprotectants to provide cellular stability. The accessory genome of the Mixed population was enriched in genes for multiple glycosyltransferases suggesting reduced energetic costs by adhering to mineral surfaces or to other microorganisms, and fostering a non-motile lifestyle. These results highlight the niche differentiation of distinct Methanobacterium populations to circumvent the challenges of serpentinization impacted fluids through coexistence strategies, supporting our ability to understand controls on methanogenic lifestyles and adaptations within the serpentinizing subsurface fluids of the Samail Ophiolite.

Ten quick tips for avoiding pitfalls in multi-omics data integration analyses

Data are the most important elements of bioinformatics: Computational analysis of bioinformatics data, in fact, can help researchers infer new knowledge about biology, chemistry, biophysics, and sometimes even medicine, influencing treatments and therapies for patients. Bioinformatics and high-throughput biological data coming from different sources can even be more helpful, because each of these different data chunks can provide alternative, complementary information about a specific biological phenomenon, similar to multiple photos of the same subject taken from different angles. In this context, the integration of bioinformatics and high-throughput biological data gets a pivotal role in running a successful bioinformatics study. In the last decades, data originating from proteomics, metabolomics, metagenomics, phenomics, transcriptomics, and epigenomics have been labelled -omics data, as a unique name to refer to them, and the integration of these omics data has gained importance in all biological areas. Even if this omics data integration is useful and relevant, due to its heterogeneity, it is not uncommon to make mistakes during the integration phases. We therefore decided to present these ten quick tips to perform an omics data integration correctly, avoiding common mistakes we experienced or noticed in published studies in the past. Even if we designed our ten guidelines for beginners, by using a simple language that (we hope) can be understood by anyone, we believe our ten recommendations should be taken into account by all the bioinformaticians performing omics data integration, including experts.

Genome-based reclassification of Bacillus acidicola, Bacillus pervagus and the genera Heyndrickxia, Margalitia and Weizmannia

In the present study, the taxonomic positions of Bacillus acidicola, Bacillus pervagus and members of the genera Heyndrickxia, Margalitia and Weizmannia were evaluated. The 16S rRNA gene sequence similarity between Bacillus acidicola DSM 14745^T, Bacillus pervagus DSM 23947^T and members of the genera Heyndrickxia and Margalitia were above the cut-off level (>95 %) for genus delineation. Amino acid identity (AAI) values and the results of phylogenomic analysis suggested that B. acidicola and the members of the genera Heyndrickxia, Margalitia and Weizmannia belong to the same genus. Furthermore, the AAI and phylogenomic results also differentiate B. pervagus from B. acidicola and the members of the genera Heyndrickxia, Margalitia and Weizmannia. Based on the results, we propose to transfer Bacillus acidicola, Margalitia and Weizmannia to the genus Heyndrickxia. We also propose the reclassification of B. pervagus into a new genus Oikeobacillus gen. nov., with the type species Oikeobacillus pervagus comb. nov.

The Western Lake Erie culture collection: A promising resource for evaluating the physiological and genetic diversity of Microcystis and its associated microbiome

Cyanobacterial harmful algal blooms (cyanoHABs) dominated by Microcystis spp. have significant public health and economic implications in freshwater bodies around the world. These blooms are capable of producing a variety of cyanotoxins, including microcystins, that affect fishing and tourism industries, human and environmental health, and access to drinking water. In this study, we isolated and sequenced the genomes of 21 primarily unialgal Microcystis cultures collected from western Lake Erie between 2017 and 2019. While some cultures isolated in different years have a high degree of genetic similarity (genomic Average Nucleotide Identity >99%), genomic data show that these cultures also represent much of the breadth of known Microcystis diversity in natural populations. Only five isolates contained all the genes required for microcystin biosynthesis while two isolates contained a previously described partial mcy operon. Microcystin production within cultures was also assessed using Enzyme-Linked Immunosorbent Assay (ELISA) and supported genomic results with high concentrations (up to 900 μg L⁻¹) in cultures with complete mcy operons and no or low toxin detected otherwise. These xenic cultures also contained a substantial diversity of bacteria associated with Microcystis, which has become increasingly recognized as an essential component of cyanoHAB community dynamics. These results highlight the genomic diversity among Microcystis strains and associated bacteria in Lake Erie, and their potential impacts on bloom development, toxin production, and toxin degradation. This culture collection significantly increases the availability of environmentally relevant Microcystis strains from temperate North America.

Arcobacteraceae comparative genome analysis demonstrates genome heterogeneity and reduction in species isolated from animals and associated with human illness

The Arcobacteraceae family groups Gram-negative bacterial species previously included in the family Campylobacteraceae. These species of which some are considered foodborne pathogens, have been isolated from different environmental niches and hosts. They have been isolated from various types of foods, though predominantly from food of animal origin, as well as from stool of humans with enteritis. Their different abilities to survive in different hosts and environments suggest an evolutionary pressure with consequent variation in their genome content. Moreover, their different physiological and genomic characteristics led to the recent proposal to create new genera within this family, which is however criticized due to the lack of discriminatory features and biological and clinical relevance. Aims of the present study were to assess the Arcobacteraceae pangenome, and to characterize existing similarities and differences in 20 validly described species. For this, analysis has been conducted on the genomes of the corresponding type strains obtained by Illumina sequencing, applying several bioinformatic tools. Results of the present study do not support the proposed division into different genera and revealed the presence of pangenome partitions with numbers comparable to other Gram-negative bacteria genera, such as Campylobacter. Different gene class compositions in animal and human-associated species are present, including a higher percentage of virulence-related gene classes such as cell motility genes. The adaptation to environmental and/or host conditions of some species was identified by the presence of specific genes. Furthermore, a division into pathogenic and non-pathogenic species is suggested, which can support future research on food safety and public health.

In situ cell division and mortality rates of SAR11, SAR86, Bacteroidetes, and Aurantivirga during phytoplankton blooms reveal differences in population controls

Net growth of microbial populations, that is, changes in abundances over time, can be studied using 16S rRNA fluorescence in situ hybridization (FISH). However, this approach does not differentiate between mortality and cell division rates. We used FISH-based image cytometry in combination with dilution culture experiments to study net growth, cell division, and mortality rates of four bacterial taxa over two distinct phytoplankton blooms: the oligotrophs SAR11 and SAR86, and the copiotrophic phylum Bacteroidetes, and its genus Aurantivirga. Cell volumes, ribosome content, and frequency of dividing cells (FDC) co-varied over time. Among the three, FDC was the most suitable predictor to calculate cell division rates for the selected taxa. The FDC-derived cell division rates for SAR86 of up to 0.8/day and Aurantivirga of up to 1.9/day differed, as expected for oligotrophs and copiotrophs. Surprisingly, SAR11 also reached high cell division rates of up to 1.9/day, even before the onset of phytoplankton blooms. For all four taxonomic groups, the abundance-derived net growth (-0.6 to 0.5/day) was about an order of magnitude lower than the cell division rates. Consequently, mortality rates were comparably high to cell division rates, indicating that about 90% of bacterial production is recycled without apparent time lag within 1 day. Our study shows that determining taxon-specific cell division rates complements omics-based tools and provides unprecedented clues on individual bacterial growth strategies including bottom-up and top-down controls. IMPORTANCE The growth of a microbial population is often calculated from their numerical abundance over time. However, this does not take cell division and mortality rates into account, which are important for deriving ecological processes like bottom-up and top-down control. In this study, we determined growth by numerical abundance and calibrated microscopy-based methods to determine the frequency of dividing cells and subsequently calculate taxon-specific cell division rates in situ. The cell division and mortality rates of two oligotrophic (SAR11 and SAR86) and two copiotrophic (Bacteroidetes and Aurantivirga) taxa during two spring phytoplankton blooms showed a tight coupling for all four taxa throughout the blooms without any temporal offset. Unexpectedly, SAR11 showed high cell division rates days before the bloom while cell abundances remained constant, which is indicative of strong top-down control. Microscopy remains the method of choice to understand ecological processes like top-down and bottom-up control on a cellular level.

An initial genomic blueprint of the healthy human oesophageal microbiome

Background

The oesophageal microbiome is thought to contribute to the pathogenesis of oesophageal cancer. However, investigations using culture and molecular barcodes have provided only a low-resolution view of this important microbial community. We therefore explored the potential of culturomics and metagenomic binning to generate a catalogue of reference genomes from the healthy human oesophageal microbiome, alongside a comparison set from saliva.

Results

Twenty-two distinct colonial morphotypes from healthy oesophageal samples were genome-sequenced. These fell into twelve species clusters, eleven of which represented previously defined species. Two isolates belonged to a novel species, which we have named Rothia gullae. We performed metagenomic binning of reads generated from UK samples from this study alongside reads generated from Australian samples in a recent study. Metagenomic binning generated 136 medium or high-quality metagenome-assembled genomes (MAGs). MAGs were assigned to 56 species clusters, eight representing novel Candidatus species, which we have named Ca. Granulicatella gullae, Ca. Streptococcus gullae, Ca. Nanosynbacter quadramensis, Ca. Nanosynbacter gullae, Ca. Nanosynbacter colneyensis, Ca. Nanosynbacter norwichensis, Ca. Nanosynococcus oralis and Ca. Haemophilus gullae. Five of these novel species belong to the recently described phylum Patescibacteria . Although members of the Patescibacteria are known to inhabit the oral cavity, this is the first report of their presence in the oesophagus. Eighteen of the metagenomic species were, until recently, identified only by hard-to-remember alphanumeric placeholder designations. Here we illustrate the utility of a set of recently published arbitrary Latinate species names in providing user-friendly taxonomic labels for microbiome analyses.Our non-redundant species catalogue contained 63 species derived from cultured isolates or MAGs. Mapping revealed that these species account for around half of the sequences in the oesophageal and saliva metagenomes. Although no species was present in all oesophageal samples, 60 species occurred in at least one oesophageal metagenome from either study, with 50 identified in both cohorts.

Conclusions

Recovery of genomes and discovery of new species represents an important step forward in our understanding of the oesophageal microbiome. The genes and genomes that we have released into the public domain will provide a base line for future comparative, mechanistic and intervention studies.

Characterization of two keystone taxa, sulfur-oxidizing, and nitrate-reducing bacteria, by tracking their role transitions in the benzo[a]pyrene degradative microbiome

BACKGROUND

Keystone taxa are drivers of microbiome structure and functioning, which may play critical roles in microbiome-level responses to recalcitrant pollution and are a key to bioremediation. However, the characterization and manipulation of such taxa is a major challenge due to the complexity of microbial communities and rapid turnover in both time and space. Here, microcosms were set up with benzo[a]-pyrene (BaP) and/or nitrate based on C-rich, S-rich, and N-limited mangrove sediments as reductive experimental models to trigger and track the turnover of keystone taxa to address this challenge.

RESULTS

Based on microbial co-occurrence network analysis, two keystone taxa, Sulfurovum and Sulfurimonas, were found to exhibit significant role transitions in different microcosms, where these two taxa played nonkeystone roles with neutral relationships in in situ mangrove sediments. However, Sulfurimonas transitioned to be keystone taxa in nitrate-replenished microcosms and formed a keystone guild with Thioalkalispira. Sulfurovum stood out in BaP-added microcosms and mutualized in a densely polycyclic aromatic hydrocarbon (PAH)-degrader-centric keystone guild with Novosphingobium and Robiginitalea, where 63.25% of added BaP was removed. Under the occurrence of nitrate and BaP, they simultaneously played roles as keystone taxa in their respective guilds but exhibited significant competition. Comparative genomics and metagenome-assembled genome (MAG) analysis was then performed to reveal the metabolic potential of those keystone taxa and to empirically deduce their functional role in keystone guilds. Sulfurimonas possesses a better sense system and motility, indicative of its aggressive role in nitrate acquisition and conversion; Sulfurovum exhibited a better ability for oxidation resistance and transporting nutrients and electrons. High-efficiency thermal asymmetric interlaced polymerase reaction (hiTAIL-PCR) and enhanced green fluorescent protein (eGFP)-labeling approaches were employed to capture and label the BaP key degrader to further experimentally verify the roles of keystone taxa Sulfurovum in the keystone guilds. Observations of the enhancement in reactive oxygen species (ROS) removal, cell growth, and degradation efficiency by co-culture of isolated keystone taxa strains experimentally demonstrated that Sulfurovum contributes to the BaP degradative microbiome against BaP toxicity.

CONCLUSIONS

Our findings suggest that the combined use of co-occurrence network analysis, comparative genomics, and co-culture of captured keystone taxa (3C-strategy) in microbial communities whose structure is strongly shaped by changing environmental factors can characterize keystone taxa roles in keystone guilds and may provide targets for manipulation to improve the function of the microbiome. Video Abstract.

Site-specialization of human oral Gemella species

Gemella species are core members of the human oral microbiome in healthy subjects and are regarded as commensals, although they can cause opportunistic infections. Our objective was to evaluate the site-specialization of Gemella species among various habitats within the mouth by combining pangenomics and metagenomics. With pangenomics, we identified genome relationships and categorized genes as core and accessory to each species. With metagenomics, we identified the primary oral habitat of individual genomes. Our results establish that the genomes of three species, G. haemolysans, G. sanguinis and G. morbillorum, are abundant and prevalent in human mouths at different oral sites: G. haemolysans on buccal mucosa and keratinized gingiva; G. sanguinis on tongue dorsum, throat, and tonsils; and G. morbillorum in dental plaque. The gene-level basis of site-specificity was investigated by identifying genes that were core to Gemella genomes at a specific oral site but absent from other Gemella genomes. The riboflavin biosynthesis pathway was present in G. haemolysans genomes associated with buccal mucosa but absent from the rest of the genomes. Overall, metapangenomics show that Gemella species have clear ecological preferences in the oral cavity of healthy humans and provides an approach to identifying gene-level drivers of site specificity.

Fusobacterium sphaericum sp. nov. , isolated from a human colon tumor, is prevalent in various human body sites and induces IL-8 secretion from colorectal cancer cells

Cancerous tissue is a largely unexplored microbial niche that provides a unique environment for the colonization and growth of specific bacterial communities, and with it, the opportunity to identify novel bacterial species. Here, we report distinct features of a novel Fusobacterium species, F. sphaericum sp. nov. ( Fs ), isolated from primary colon adenocarcinoma tissue. We acquire the complete, closed genome of this organism and phylogenetically confirm its classification into the Fusobacterium genus. Phenotypic and genomic analysis of Fs reveal that this novel organism is of coccoid shape, rare for Fusobacterium members, and has species-distinct gene content. Fs displays a metabolic profile and antibiotic resistance repertoire consistent with other Fusobacterium species. In vitro, Fs has adherent and immunomodulatory capabilities, as it intimately associates with human colon cancer epithelial cells and promotes IL-8 secretion. Analysis of the prevalence and abundance of Fs in ∼1,750 human metagenomic samples shows that it is a moderately prevalent member of the human oral cavity and stool. Intriguingly, analysis of ∼1,270 specimens from patients with colorectal cancer demonstrate that Fs is significantly enriched in colonic and tumor tissue as compared to mucosa or feces. Our study sheds light on a novel bacterial species that is prevalent within the human intestinal microbiota and whose role in human health and disease requires further investigation.

Comparative genomic analysis of Cohnella hashimotonis sp. nov. isolated from the International Space Station

A single strain from the family Paenibacillaceae was isolated from the wall behind the Waste Hygiene Compartment aboard the International Space Station (ISS) in April 2018, as part of the Microbial Tracking mission series. This strain was identified as a gram-positive, rod-shaped, oxidase-positive, catalase-negative motile bacterium in the genus Cohnella, designated as F6_2S_P_1T. The 16S sequence of the F6_2S_P_1T strain places it in a clade with C. rhizosphaerae and C. ginsengisoli, which were originally isolated from plant tissue or rhizosphere environments. The closest 16S and gyrB matches to strain F6_2S_P_1T are to C. rhizosphaerae with 98.84 and 93.99% sequence similarity, while a core single-copy gene phylogeny from all publicly available Cohnella genomes places it as more closely related to C. ginsengisoli. Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values to any described Cohnella species are <89 and <22%, respectively. The major fatty acids for strain F6_2S_P_1T are anteiso-C15:0 (51.7%), iso-C16:0 (23.1%), and iso-C15:0 (10.5%), and it is able to metabolize a wide range of carbon compounds. Given the results of the ANI and dDDH analyses, this ISS strain is a novel species within the genus Cohnella for which we propose the name Cohnella hashimotonis, with the type strain F6_2S_P_1T (=NRRL B-65657T and DSMZ 115098T). Because no closely related Cohnella genomes were available, this study generated the whole-genome sequences (WGSs) of the type strains for C. rhizosphaerae and C. ginsengisoli. Phylogenetic and pangenomic analysis reveals that F6_2S_P_1T, C. rhizosphaerae, and C. ginsengisoli, along with two uncharacterized Cohnella strains, possess a shared set of 332 gene clusters which are not shared with any other WGS of Cohnella species, and form a distinct clade branching off from C. nanjingensis. Functional traits were predicted for the genomes of strain F6_2S_P_1T and other members of this clade.

Inter-kingdom interactions and stability of methanogens revealed by machine-learning guided multi-omics analysis of industrial-scale biogas plants

Multi-omics analysis is a powerful tool for the detection and study of inter-kingdom interactions, such as those between bacterial and archaeal members of complex biogas-producing microbial communities. In the present study, the microbiomes of three industrial-scale biogas digesters, each fed with different substrates, were analysed using a machine-learning guided genome-centric metagenomics framework complemented with metatranscriptome data. This data permitted us to elucidate the relationship between abundant core methanogenic communities and their syntrophic bacterial partners. In total, we detected 297 high-quality, non-redundant metagenome-assembled genomes (nrMAGs). Moreover, the assembled 16 S rRNA gene profiles of these nrMAGs showed that the phylum Firmicutes possessed the highest copy number, while the representatives of the archaeal domain had the lowest. Further investigation of the three anaerobic microbial communities showed characteristic alterations over time but remained specific to each industrial-scale biogas plant. The relative abundance of various microorganisms as revealed by metagenome data was independent from corresponding metatranscriptome activity data. Archaea showed considerably higher activity than was expected from their abundance. We detected 51 nrMAGs that were present in all three biogas plant microbiomes with different abundances. The core microbiome correlated with the main chemical fermentation parameters, and no individual parameter emerged as a predominant shaper of community composition. Various interspecies H₂/electron transfer mechanisms were assigned to hydrogenotrophic methanogens in the biogas plants that ran on agricultural biomass and wastewater. Analysis of metatranscriptome data revealed that methanogenesis pathways were the most active of all main metabolic pathways.

Ecology of Endozoicomonadaceae in three coral genera across the Pacific Ocean

Health and resilience of the coral holobiont depend on diverse bacterial communities often dominated by key marine symbionts of the Endozoicomonadaceae family. The factors controlling their distribution and their functional diversity remain, however, poorly known. Here, we study the ecology of Endozoicomonadaceae at an ocean basin-scale by sampling specimens from three coral genera (Pocillopora, Porites, Millepora) on 99 reefs from 32 islands across the Pacific Ocean. The analysis of 2447 metabarcoding and 270 metagenomic samples reveals that each coral genus harbored a distinct new species of Endozoicomonadaceae. These species are composed of nine lineages that have distinct biogeographic patterns. The most common one, found in Pocillopora, appears to be a globally distributed symbiont with distinct metabolic capabilities, including the synthesis of amino acids and vitamins not produced by the host. The other lineages are structured partly by the host genetic lineage in Pocillopora and mainly by the geographic location in Porites. Millepora is more rarely associated to Endozoicomonadaceae. Our results show that different coral genera exhibit distinct strategies of host-Endozoicomonadaceae associations that are defined at the bacteria lineage level.

uBin: A manual refining tool for genomes from metagenomes

Resolving bacterial and archaeal genomes from metagenomes has revolutionized our understanding of Earth's biomes yet producing high-quality genomes from assembled fragments has been an ever-standing problem. While automated binning software and their combination produce prokaryotic bins in high throughput, their manual refinement has been slow, sometimes difficult or missing entirely facilitating error propagation in public databases. Here, we present uBin, a GUI-based, standalone bin refiner that runs on all major operating platforms and was additionally designed for educational purposes. When applied to the public CAMI dataset, refinement of bins using GC content, coverage and taxonomy was able to improve 78.9% of bins by decreasing their contamination. We also applied the bin refiner as a standalone binner to public metagenomes from the International Space Station and demonstrate the recovery of near-complete genomes, whose replication indices indicate the active proliferation of microbes in Earth's lower orbit. uBin is an easy to instal software for bin refinement, binning of simple metagenomes and communication of metagenomic results to other scientists and in classrooms. The software and its helper scripts are open source and available under https://github.com/ProbstLab/uBin.

Candidatus Alkanophaga archaea from Guaymas Basin hydrothermal vent sediment oxidize petroleum alkanes

Methanogenic and methanotrophic archaea produce and consume the greenhouse gas methane, respectively, using the reversible enzyme methyl-coenzyme M reductase (Mcr). Recently, Mcr variants that can activate multicarbon alkanes have been recovered from archaeal enrichment cultures. These enzymes, called alkyl-coenzyme M reductase (Acrs), are widespread in the environment but remain poorly understood. Here we produced anoxic cultures degrading mid-chain petroleum n-alkanes between pentane (C₅) and tetradecane (C₁₄) at 70 °C using oil-rich Guaymas Basin sediments. In these cultures, archaea of the genus Candidatus Alkanophaga activate the alkanes with Acrs and completely oxidize the alkyl groups to CO₂. Ca. Alkanophaga form a deep-branching sister clade to the methanotrophs ANME-1 and are closely related to the short-chain alkane oxidizers Ca. Syntrophoarchaeum. Incapable of sulfate reduction, Ca. Alkanophaga shuttle electrons released from alkane oxidation to the sulfate-reducing Ca. Thermodesulfobacterium syntrophicum. These syntrophic consortia are potential key players in petroleum degradation in heated oil reservoirs.

Epiphytic common core bacteria in the microbiomes of co-located green (Ulva), brown (Saccharina) and red (Grateloupia, Gelidium) macroalgae

BACKGROUND

Macroalgal epiphytic microbial communities constitute a rich resource for novel enzymes and compounds, but studies so far largely focused on tag-based microbial diversity analyses or limited metagenome sequencing of single macroalgal species.

RESULTS

We sampled epiphytic bacteria from specimens of Ulva sp. (green algae), Saccharina sp. (brown algae), Grateloupia sp. and Gelidium sp. (both red algae) together with seawater and sediment controls from a coastal reef in Weihai, China, during all seasons. Using 16S rRNA amplicon sequencing, we identified 14 core genera (consistently present on all macroalgae), and 14 dominant genera (consistently present on three of the macroalgae). Core genera represented ~ 0.7% of all genera, yet accounted for on average 51.1% of the bacterial abundances. Plate cultivation from all samples yielded 5,527 strains (macroalgae: 4,426) representing 1,235 species (685 potentially novel). Sequencing of selected strains yielded 820 non-redundant draft genomes (506 potentially novel), and sequencing of 23 sampled metagenomes yielded 1,619 metagenome-assembled genomes (MAGs), representing further 1,183 non-redundant genomes. 230 isolates and 153 genomes were obtained from the 28 core/dominant genera. We analyzed the genomic potential of phycosphere bacteria to degrade algal polysaccharides and to produce bioactive secondary metabolites. We predicted 4,451 polysaccharide utilization loci (PULs) and 8,810 biosynthetic gene clusters (BGCs). These were particularly prevalent in core/dominant genera.

CONCLUSIONS

Our metabolic annotations and analyses of MAGs and genomes provide new insights into novel species of phycosphere bacteria and their ecological niches for an improved understanding of the macroalgal phycosphere microbiome. Video Abstract.

Effect of Arsenic on EPS Synthesis, Biofilm Formation, and Plant Growth-Promoting Abilities of the Endophytes Pseudomonas PD9R and Rahnella laticis PD12R

Phytoremediation, a cost-effective, eco-friendly alternative to conventional remediation, could expand efforts to remediate arsenic-contaminated soils. As with other pollutants, the plant microbiome may improve phytoremediation outcomes for arsenic-contaminated sites. We used in vitro and in silico methods to compare the arsenic resistance mechanisms, synthesis of extracellular polymeric substances (EPS), biofilm formation, and plant growth-promoting abilities of the endophytes Pseudomonas sp. PD9R and Rahnella laticis PD12R. PD12R, which tolerates arsenate (As(V)) and arsenite (As(III)) to concentrations fivefold greater than PD9R, synthesizes high volumes of EPS in response to arsenic, and sequesters arsenic in the capsular EPS and cells. While arsenic exposure induced EPS synthesis in both strains, only PD12R continued to form biofilms at high As(III) and As(V) concentrations. The effects of endophyte inoculation on Arabidopsis growth varied by strain and As(V) concentration, and PD9R had positive effect on plants exposed to low levels of arsenic. Comparative genomic analyses exploring the EPS synthesis and arsenic resistance mechanisms against other Pseudomonas and Rahnella strains suggest that both strains possess atypical arsenic resistance mechanisms from other plant-associated strains, while the configuration of the EPS synthesis systems appeared to be more broadly distributed among plant- and non-plant-associated strains.

Expression of Microcystis Biosynthetic Gene Clusters in Natural Populations Suggests Temporally Dynamic Synthesis of Novel and Known Secondary Metabolites in Western Lake Erie

Microcystis spp. produce diverse secondary metabolites within freshwater cyanobacterial harmful algal blooms (cyanoHABs) around the world. In addition to the biosynthetic gene clusters (BGCs) encoding known compounds, Microcystis genomes harbor numerous BGCs of unknown function, indicating a poorly understood chemical repertoire. While recent studies show that Microcystis produces several metabolites in the lab and field, little work has focused on analyzing the abundance and expression of its broader suite of BGCs during cyanoHAB events. Here, we use metagenomic and metatranscriptomic approaches to track the relative abundance of Microcystis BGCs and their transcripts throughout the 2014 western Lake Erie cyanoHAB. The results indicate the presence of several transcriptionally active BGCs that are predicted to synthesize both known and novel secondary metabolites. The abundance and expression of these BGCs shifted throughout the bloom, with transcript abundance levels correlating with temperature, nitrate, and phosphorus concentrations and the abundance of co-occurring predatory and competitive eukaryotic microorganisms, suggesting the importance of both abiotic and biotic controls in regulating expression. This work highlights the need for understanding the chemical ecology and potential risks to human and environmental health posed by secondary metabolites that are produced but often unmonitored. It also indicates the prospects for identifying pharmaceutical-like molecules from cyanoHAB-derived BGCs. IMPORTANCE Microcystis spp. dominate cyanobacterial harmful algal blooms (cyanoHABs) worldwide and pose significant threats to water quality through the production of secondary metabolites, many of which are toxic. While the toxicity and biochemistry of microcystins and several other compounds have been studied, the broader suite of secondary metabolites produced by Microcystis remains poorly understood, leaving gaps in our understanding of their impacts on human and ecosystem health. We used community DNA and RNA sequences to track the diversity of genes encoding synthesis of secondary metabolites in natural Microcystis populations and assess patterns of transcription in western Lake Erie cyanoHABs. Our results reveal the presence of both known gene clusters that encode toxic secondary metabolites as well as novel ones that may encode cryptic compounds. This research highlights the need for targeted studies of the secondary metabolite diversity in western Lake Erie, a vital freshwater source to the United States and Canada.

Distribution and diversity of 'Tectomicrobia', a deep-branching uncultivated bacterial lineage harboring rich producers of bioactive metabolites

Genomic and functional analyses of bacterial sponge symbionts belonging to the uncultivated candidate genus 'Entotheonella' has revealed them as the prolific producers of bioactive compounds previously identified from their invertebrate hosts. These studies also suggested 'Entotheonella' as the first members of a new candidate phylum, 'Tectomicrobia'. Here we analyzed the phylogenetic structure and environmental distribution of this as-yet sparsely populated phylum-like lineage. The data show that 'Entotheonella' and other 'Tectomicrobia' are not restricted to marine habitats but widely distributed among terrestrial locations. The inferred phylogenetic trees suggest several intra-phylum lineages with diverse lifestyles. Of these, the previously described 'Entotheonella' lineage can be more accurately divided into at least three different candidate genera with the terrestrial 'Candidatus Prasianella', the largely terrestrial 'Candidatus Allonella', the 'Candidatus Thalassonella' comprising sponge-associated members, and the more widely distributed 'Candidatus Entotheonella'. Genomic characterization of 'Thalassonella' members from a range of sponge hosts did not suggest a role as providers of natural products, despite high genomic similarity to 'Entotheonella' regarding primary metabolism and implied lifestyle. In contrast, the analysis revealed a correlation between the revised 'Entotheonella' 16S rRNA gene phylogeny and a specific association with sponges and their natural products. This feature might serve as a discovery method to accelerate the identification of new chemically rich 'Entotheonella' variants, and led to the identification of the first 'Entotheonella' symbiont in a non-tetractinellid sponge, Psammocinia sp., indicating a wide host distribution of 'Entotheonella'-based chemical symbiosis.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

Cataloging the Phylogenetic Diversity of Human Bladder Bacterial Isolates

Although the human bladder is reported to harbor unique microbiota, our understanding of how these microbial communities interact with their human hosts is limited, mostly owing to the lack of isolates to test mechanistic hypotheses. Niche-specific bacterial collections and associated reference genome databases have been instrumental in expanding knowledge of the microbiota of other anatomical sites, e.g., the gut and oral cavity. To facilitate genomic, functional, and experimental analyses of the human bladder microbiota, here we present a bladder-specific bacterial reference collection comprised of 1134 genomes. These genomes were culled from bacterial isolates obtained by a metaculturomic method from bladder urine collected by transurethral catheterization. This bladder-specific bacterial reference collection includes 196 different species, including representatives of major aerobes and facultative anaerobes, as well as some anaerobes. It captures 72.2 % of the genera found when we reexamined previously published 16S rRNA gene sequencing of 392 adult female bladder urine samples. Comparative genomic analysis found that the taxonomies and functions of the bladder microbiota shared more similarities with the vaginal microbiota than the gut microbiota. Whole-genome phylogenetic and functional analyses of 186 bladder E. coli isolates and 387 gut E. coli isolates supports the hypothesis that phylogroup distribution and functions of E. coli strains differ dramatically between these two very different niches. This bladder-specific bacterial reference collection is a unique resource that will enable hypothesis-driven bladder microbiota research and comparison to isolates from other anatomical sites.

Draft genomes, phylogenomic reconstruction and comparative genome analysis of three Xenorhabdus strains isolated from soil-dwelling nematodes in Kenya

As a proven source of potent and selective antimicrobials, Xenorhabdus bacteria are important to an age plagued with difficult-to-treat microbial infections. Yet, only 27 species have been described to date. In this study, a novel Xenorhabdus species was discovered through genomic studies on three isolates from Kenyan soils. Soils in Western Kenya were surveyed for steinernematids and Steinernema isolates VH1 and BG5 were recovered from red volcanic loam soils from cultivated land in Vihiga and clay soils from riverine land in Bungoma respectively. From the two nematode isolates, Xenorhabdus sp. BG5 and Xenorhabdus sp. VH1 were isolated. The genomes of these two, plus that of X. griffiniae XN45 - this was previously isolated from Steinernema sp. scarpo that also originated from Kenyan soils - were sequenced and assembled. Nascent genome assemblies of the three isolates were of good quality with over 70 % of their proteome having known functions. These three isolates formed the X. griffiniae clade in a phylogenomic reconstruction of the genus. Their species were delineated using three overall genome relatedness indices: an unnamed species of the genus, Xenorhabdus sp. BG5, X. griffiniae VH1 and X. griffiniae XN45. A pangenome analysis of this clade revealed that over 70 % of species-specific genes encoded unknown functions. Transposases were linked to genomic islands in Xenorhabdus sp. BG5. Thus, overall genome-related indices sufficiently delineated species of two new Xenorhabdus isolates from Kenya, both of which were closely related to X. griffiniae . The functions encoded by most species-specific genes in the X. griffiniae clade remain unknown.

Assessing the efficacy of probiotics in augmenting bovine reproductive health: an integrated in vitro, in silico, and in vivo study

The aim of this study was to isolate and characterize bovine-vaginal probiotics genotypically and phenotypically using in silico and evaluate their in vivo performance in buffaloes with endometritis. For the in vitro isolation and characterization, vaginal swabs were collected from 34 cows and 17 buffaloes, and 709 primary bacterial isolates with probiotic activity were obtained using MRS agar media. Two isolates Lactiplantibacillus plantarum KUGBRC (LPKUGBRC) and Pediococcus pentosaceus GBRCKU (PPGBRCKU) demonstrated optimum in vitro probiotic activities as compared to Lacticaseibacillus rhamnosus GG including, acid production, secretion of fatty acids and exopolysaccharide, cell surface hydrophobicity, self-aggregating and co-aggregating capacity with pathogens, anti-microbial activity and bacteriocin-like compounds against pathogens Escherichia coli and Staphylococcus aureus in cell-free supernatant and absence of hemolytic activity. Their phenotypic capacity was confirmed by analyzing the whole genome sequencing data and identifying genes and pathways associated with probiotic properties. These probiotic isolates have shown no virulence genes were discovered in their genomic study. In vivo study of 92 buffaloes suffering from clinical endometritis with purulent cervico-vaginal mucus (CVM) were randomly allocated 40 × 10⁸ CFU/ml LPKUGBRC and PPGBRCKU and 40 ml Normal saline. The LPKUGBRC reduced the duration between administration of probiotic to induction of healthy estrus significantly. However, no effect was observed on pregnancy rate. These results suggest that LPKUGBRC and PPGBRCKU probiotic bacteria demonstrate probiotic efficiency and adaptability. Further sourced from the same niche as the targeted infection, they offer a distinct advantage in targeting the specific microbial population associated with endometritis. The findings of this study highlight the potential of LPKUGBRC and PPGBRCKU probiotics in treating endometritis and suggest further exploration of their clinical applications.