GToTree: a user-friendly workflow for phylogenomics

Identification of microbial metabolic functional guilds from large genomic datasets

Heterotrophic microbes play an important role in the Earth System as key drivers of major biogeochemical cycles. Specifically, the consumption rate of organic matter is set by the interaction between diverse microbial communities and the chemical and physical environment in which they reside. Modeling these dynamics requires reducing the complexity of microbial communities and linking directly with biogeochemical functions. Microbial metabolic functional guilds provide one approach for reducing microbial complexity and incorporating microbial biogeochemical functions into models. However, we lack a way to identify these guilds. In this study, we present a method for defining metabolic functional guilds from annotated genomes, which are derived from both uncultured and cultured organisms. This method utilizes an Aspect Bernoulli (AB) model and was tested on three large genomic datasets with 1,733-3,840 genomes each. Ecologically relevant microbial metabolic functional guilds were identified including guilds related to DMSP degradation, dissimilatory nitrate reduction to ammonia, and motile copiotrophy. This method presents a way to generate hypotheses about functions co-occurring within individual microbes without relying on cultured representatives. Applying the concept of metabolic functional guilds to environmental samples will provide new insight into the role that heterotrophic microbial communities play in setting rates of carbon cycling.

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.

Infrequent intrahousehold transmission of Clostridioides difficile between pet owners and their pets

Companion animals have been shown to carry Clostridioides difficile strains that are similar or identical to strains found in people, and a small number of studies have shown that pets carry genetically identical C. difficile isolates as their owners, suggesting inter-species transmission. However, the directionality of transmission is ultimately unknown, and the frequency with which animals acquire C. difficile following their owners' infection is unclear. The goal of this study was to assess how often pets belonging to people with C. difficile infection carry genetically related C. difficile isolates. We enrolled pet owners from two medical institutions (University of Pennsylvania Health System (UPHS) and The Ohio State University Wexner Medical Center (OSUWMC)) who had diarrhoea with or without positive C. difficile assays and tested their faeces and their pets' faeces for C. difficile using both anaerobic culture and PCR assays. When microorganisms were obtained from both the owner and pet and had the same toxin profile or ribotype, isolates underwent genomic sequencing. Faecal samples were obtained from a total of 59 humans, 72 dogs and 9 cats, representing 47 complete households (i.e. where a sample was available from the owner and at least one pet). Of these, C. difficile was detected in 30 humans, 10 dogs and 0 cats. There were only two households where C. difficile was detected in both the owner and pet. In one of these households, the C. difficile isolates were of different toxin profiles/ribotypes (A+/B+ / RT 499 from the owner, A-/B- / RT PR22386 from the dog). In the other household, the isolates were genetically identical (one SNP difference). Interestingly, the dog from this household had recently received a course of antibiotics (cefpodoxime and metronidazole). Our findings suggest that inter-species transmission of C. difficile occurs infrequently in households with human C. difficile infections.

Targeted IS-element sequencing uncovers transposition dynamics during selective pressure in enterococci

Insertion sequences (IS) are simple transposons implicated in the genome evolution of diverse pathogenic bacterial species. Enterococci have emerged as important human intestinal pathogens with newly adapted virulence potential and antibiotic resistance. These genetic features arose in tandem with large-scale genome evolution mediated by mobile elements. Pathoadaptation in enterococci is thought to be mediated in part by the IS element IS256 through gene inactivation and recombination events. However, the regulation of IS256 and the mechanisms controlling its activation are not well understood. Here, we adapt an IS256-specfic deep sequencing method to describe how chronic lytic phage infection drives widespread diversification of IS256 in E. faecalis and how antibiotic exposure is associated with IS256 diversification in E. faecium during a clinical human infection. We show through comparative genomics that IS256 is primarily found in hospital-adapted enterococcal isolates. Analyses of IS256 transposase gene levels reveal that IS256 mobility is regulated at the transcriptional level by multiple mechanisms in E. faecalis, indicating tight control of IS256 activation in the absence of selective pressure. Our findings reveal that stressors such as phages and antibiotic exposure drives rapid genome-scale transposition in the enterococci. IS256 diversification can therefore explain how selective pressures mediate evolution of the enterococcal genome, ultimately leading to the emergence of dominant nosocomial lineages that threaten human health.

Physiological changes during torpor favor association with Endozoicomonas endosymbionts in the urochordate Botrylloides leachii

Environmental perturbations evoke down-regulation of metabolism in some multicellular organisms, leading to dormancy, or torpor. Colonies of the urochordate Botrylloides leachii enter torpor in response to changes in seawater temperature and may survive for months as small vasculature remnants that lack feeding and reproductive organs but possess torpor-specific microbiota. Upon returning to milder conditions, the colonies rapidly restore their original morphology, cytology and functionality while harboring re-occurring microbiota, a phenomenon that has not been described in detail to date. Here we investigated the stability of B. leachii microbiome and its functionality in active and dormant colonies, using microscopy, qPCR, in situ hybridization, genomics and transcriptomics. A novel lineage of Endozoicomonas, proposed here as Candidatus Endozoicomonas endoleachii, was dominant in torpor animals (53-79% read abundance), and potentially occupied specific hemocytes found only in torpid animals. Functional analysis of the metagenome-assembled genome and genome-targeted transcriptomics revealed that Endozoicomonas can use various cellular substrates, like amino acids and sugars, potentially producing biotin and thiamine, but also expressing various features involved in autocatalytic symbiosis. Our study suggests that the microbiome can be linked to the metabolic and physiological states of the host, B. leachii, introducing a model organism for the study of symbioses during drastic physiological changes, such as torpor.

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.

Metabolic handoffs between multiple symbionts may benefit the deep-sea bathymodioline mussels

Bathymodioline mussels rely on thiotrophic and/or methanotrophic chemosynthetic symbionts for nutrition, yet, secondary heterotrophic symbionts are often present and play an unknown role in the fitness of the organism. The bathymodioline Idas mussels that thrive in gas seeps and on sunken wood in the Mediterranean Sea and the Atlantic Ocean, host at least six symbiont lineages that often co-occur. These lineages include the primary symbionts chemosynthetic methane- and sulfur-oxidizing gammaproteobacteria, and the secondary symbionts, Methylophagaceae, Nitrincolaceae and Flavobacteriaceae, whose physiology and metabolism are obscure. Little is known about if and how these symbionts interact or exchange metabolites. Here we curated metagenome-assembled genomes of Idas modiolaeformis symbionts and used genome-centered metatranscriptomics and metaproteomics to assess key symbiont functions. The Methylophagaceae symbiont is a methylotrophic autotroph, as it encoded and expressed the ribulose monophosphate and Calvin-Benson-Bassham cycle enzymes, particularly RuBisCO. The Nitrincolaceae ASP10-02a symbiont likely fuels its metabolism with nitrogen-rich macromolecules and may provide the holobiont with vitamin B12. The Urechidicola (Flavobacteriaceae) symbionts likely degrade glycans and may remove NO. Our findings indicate that these flexible associations allow for expanding the range of substrates and environmental niches, via new metabolic functions and handoffs.

Novel biological aqua crust enhances in situ metal(loid) bioremediation driven by phototrophic/diazotrophic biofilm

BACKGROUND

Understanding the ecological and environmental functions of phototrophic biofilms in the biological crust is crucial for improving metal(loid) (e.g. Cd, As) bioremediation in mining ecosystems. In this study, in combination with metal(loid) monitoring and metagenomic analysis, we systematically evaluated the effect of biofilm in a novel biological aqua crust (biogenic aqua crust-BAC) on in situ metal(loid) bioremediation of a representative Pb/Zn tailing pond.

RESULTS

We observed strong accumulation of potentially bioavailable metal(loid)s and visible phototrophic biofilms in the BAC. Furthermore, dominating taxa Leptolyngbyaceae (10.2-10.4%, Cyanobacteria) and Cytophagales (12.3-22.1%, Bacteroidota) were enriched in biofilm. Along with predominant heterotrophs (e.g. Cytophagales sp.) as well as diazotrophs (e.g. Hyphomonadaceae sp.), autotrophs/diazotrophs (e.g. Leptolyngbyaceae sp.) in phototrophic biofilm enriched the genes encoding extracellular peptidase (e.g. family S9, S1), CAZymes (e.g. CBM50, GT2) and biofilm formation (e.g. OmpR, CRP and LuxS), thus enhancing the capacity of nutrient accumulation and metal(loid) bioremediation in BAC system.

CONCLUSIONS

Our study demonstrated that a phototrophic/diazotrophic biofilm constitutes the structured communities containing specific autotrophs (e.g. Leptolyngbyaceae sp.) and heterotrophs (e.g. Cytophagales sp.), which effectively control metal(loid) and nutrient input using solar energy in aquatic environments. Elucidation of the mechanisms of biofilm formation coupled with metal(loid) immobilization in BAC expands the fundamental understanding of the geochemical fate of metal(loid)s, which may be harnessed to enhance in situ metal(loid) bioremediation in the aquatic ecosystem of the mining area. Video Abstract.

Transcriptionally active nitrogen fixation and biosynthesis of diverse secondary metabolites by Dolichospermum and Aphanizomenon-like Cyanobacteria in western Lake Erie Microcystis blooms

Cyanobacterial harmful algal blooms (cyanoHABs) in the western basin of Lake Erie are dominated by microcystin producing Microcystis spp., but other cyanobacterial taxa that coexist in these communities may play important roles in production of toxins and shaping bloom dynamics and community function. In this study, we used metagenomic and metatranscriptomic data from the 2014 western Lake Erie cyanoHAB to explore the genetic diversity and biosynthetic potential of cyanobacteria belonging to the Anabaena, Dolichospermum, Aphanizomenon (ADA) clade. We reconstructed two near-complete metagenome-assembled genomes from two distinct ADA clade species, each containing biosynthetic gene clusters that encode novel and known secondary metabolites, including those with toxic and/or known taste and odor properties, that were transcriptionally active. However, neither ADA metagenome-assembled genome contained genes encoding guanitoxins, anatoxins, or saxitoxins, which are known to be produced by ADA. The ADA cyanobacteria accounted for most of the metagenomic and metatranscriptomic reads from nitrogen fixation genes, suggesting they were the dominant N-fixers at the times and stations sampled. Despite their relatively low abundance, our results highlight the possibility that ADA taxa could influence the water quality and ecology of Microcystis blooms, although the extent of these impacts remains to be quantified.

Functional annotation and importance of marine bacterial transporters of plankton exometabolites

Metabolite exchange within marine microbial communities transfers carbon and other major elements through global cycles and forms the basis of microbial interactions. Yet lack of gene annotations and concern about the quality of existing ones remain major impediments to revealing currencies of carbon flux. We employed an arrayed mutant library of the marine bacterium Ruegeria pomeroyi DSS-3 to experimentally annotate substrates of organic compound transporter systems, using mutant growth and compound drawdown analyses to link transporters to their cognate substrates. Mutant experiments verified substrates for thirteen R. pomeroyi transporters. Four were previously hypothesized based on gene expression data (taurine, glucose/xylose, isethionate, and cadaverine/putrescine/spermidine); five were previously hypothesized based on homology to experimentally annotated transporters in other bacteria (citrate, glycerol, N-acetylglucosamine, fumarate/malate/succinate, and dimethylsulfoniopropionate); and four had no previous annotations (thymidine, carnitine, cysteate, and 3-hydroxybutyrate). These bring the total number of experimentally-verified organic carbon influx transporters to 18 of 126 in the R. pomeroyi genome. In a longitudinal study of a coastal phytoplankton bloom, expression patterns of the experimentally annotated transporters linked them to different stages of the bloom, and also led to the hypothesis that citrate and 3-hydroxybutyrate were among the most highly available bacterial substrates. Improved functional annotation of the gatekeepers of organic carbon uptake is critical for deciphering carbon flux and fate in microbial ecosystems.

Expression of macromolecular organic nitrogen degrading enzymes identifies potential mediators of soil organic N availability to an annual grass

Nitrogen (N) is frequently limiting to plant growth, in part because most soil N is present as polymeric organic compounds that are not readily taken up by plants. Microbial depolymerization of these large macromolecular N-substrates gradually releases available inorganic N. While many studies have researched and modeled controls on soil organic matter formation and bulk N mineralization, the ecological-spatial, temporal and phylogenetic-patterns underlying organic N degradation remain unclear. We analyzed 48 time-resolved metatranscriptomes and quantified N-depolymerization gene expression to resolve differential expression by soil habitat and time in specific taxonomic groups and gene-based guilds. We observed much higher expression of extracellular serine-type proteases than other extracellular N-degrading enzymes, with protease expression of predatory bacteria declining with time and other taxonomic patterns driven by the presence (Gammaproteobacteria) or absence (Thermoproteota) of live roots and root detritus (Deltaproteobacteria and Fungi). The primary chitinase chit1 gene was more highly expressed by eukaryotes near root detritus, suggesting predation of fungi. In some lineages, increased gene expression over time suggests increased competitiveness with rhizosphere age (Chloroflexi). Phylotypes from some genera had protease expression patterns that could benefit plant N nutrition, for example, we identified a Janthinobacterium phylotype and two Burkholderiales that depolymerize organic N near young roots and a Rhizobacter with elevated protease levels near mature roots. These taxon-resolved gene expression results provide an ecological read-out of microbial interactions and controls on N dynamics in specific soil microhabitats and could be used to target potential plant N bioaugmentation strategies.

Genome-resolving metagenomics reveals wild western capercaillies (Tetrao urogallus) as avian hosts for antibiotic-resistance bacteria and their interactions with the gut-virome community

The gut microbiome is a critical component of avian health, influencing nutrient uptake and immune functions. While the gut microbiomes of agriculturally important birds have been studied, the microbiomes of wild birds still need to be explored. Filling this knowledge gap could have implications for the microbial rewilding of captive birds and managing avian hosts for antibiotic-resistant bacteria (ARB). Using genome-resolved metagenomics, we recovered 112 metagenome-assembled genomes (MAGs) from the faeces of wild and captive western capercaillies (Tetrao urogallus) (n = 8). Comparisons of bacterial diversity between the wild and captive capercaillies suggest that the reduced diversity in the captive individual could be due to differences in diet. This was further substantiated through the analyses of 517,657 clusters of orthologous groups (COGs), which revealed that gene functions related to amino acids and carbohydrate metabolisms were more abundant in wild capercaillies. Metagenomics mining of resistome identified 751 antibiotic resistance genes (ARGs), of which 40.7 % were specific to wild capercaillies suggesting that capercaillies could be potential reservoirs for hosting ARG-associated bacteria. Additionally, the core resistome shared between wild and captive capercaillies indicates that birds can acquire these ARG-associated bacteria naturally from the environment (43.1 % of ARGs). The association of 26 MAGs with 120 ARGs and 378 virus operational taxonomic units (vOTUs) also suggests a possible interplay between these elements, where putative phages could have roles in modulating the gut microbiota of avian hosts. These findings can have important implications for conservation and human health, such as avian gut microbiota rewilding, identifying the emerging threats or opportunities due to phage-microbe interactions, and monitoring the potential spread of ARG-associated bacteria from wild avian populations.

Evolutionary investigations of the biosynthetic diversity in the skin microbiome using lsaBGC

Bacterial secondary metabolites, synthesized by enzymes encoded in biosynthetic gene clusters (BGCs), can underlie microbiome homeostasis and serve as commercialized products, which have historically been mined from a select group of taxa. While evolutionary approaches have proven beneficial for prioritizing BGCs for experimental characterization efforts to uncover new natural products, dedicated bioinformatics tools designed for comparative and evolutionary analysis of BGCs within focal taxa are limited. We thus developed l ineage s pecific a nalysis of BGCs (lsa BGC; https://github.com/Kalan-Lab/lsaBGC ) to aid exploration of microdiversity and evolutionary trends across homologous groupings of BGCs, gene cluster families (GCFs), in any bacterial taxa of interest. lsa BGC enables rapid and direct identification of GCFs in genomes, calculates evolutionary statistics and conservation for BGC genes, and builds a framework to allow for base resolution mining of novel variants through metagenomic exploration. Through application of the suite to four genera commonly found in skin microbiomes, we uncover new insights into the evolution and diversity of their BGCs. We show that the BGC of the virulence-associated carotenoid staphyloxanthin in Staphylococcus aureus is ubiquitous across the genus Staphylococcus . While one GCF encoding the biosynthesis of staphyloxanthin showcases evidence for plasmid-mediated horizontal gene transfer (HGT) between species, another GCF appears to be transmitted vertically amongst a sub-clade of skin-associated Staphylococcus . Further, the latter GCF, which is well conserved in S. aureus , has been lost in most Staphylococcus epidermidis , which is the most common Staphylococcus species on human skin and is also regarded as a commensal. We also identify thousands of novel single-nucleotide variants (SNVs) within BGCs from the Corynebacterium tuberculostearicum sp. complex, a narrow, multi-species clade that features the most prevalent Corynebacterium in healthy skin microbiomes. Although novel SNVs were approximately 10 times as likely to correspond to synonymous changes when located in the top five percentile of conserved sites, lsa BGC identified SNVs that defied this trend and are predicted to underlie amino acid changes within functionally key enzymatic domains. Ultimately, beyond supporting evolutionary investigations of BGCs, lsa BGC also provides important functionalities to aid efforts for the discovery or directed modification of natural products.

Co-Occurrence and Cooperation between Comammox and Anammox Bacteria in a Full-Scale Attached Growth Municipal Wastewater Treatment Process

Cooperation between comammox and anammox bacteria for nitrogen removal has been recently reported in laboratory-scale systems, including synthetic community constructs; however, there are no reports of full-scale municipal wastewater treatment systems with such cooperation. Here, we report intrinsic and extant kinetics as well as genome-resolved community characterization of a full-scale integrated fixed film activated sludge (IFAS) system where comammox and anammox bacteria co-occur and appear to drive nitrogen loss. Intrinsic batch kinetic assays indicated that majority of the aerobic ammonia oxidation was driven by comammox bacteria (1.75 ± 0.08 mg-N/g TS-h) in the attached growth phase, with minimal contribution by ammonia-oxidizing bacteria. Interestingly, a portion of total inorganic nitrogen (∼8%) was consistently lost during these aerobic assays. Aerobic nitrite oxidation assays eliminated the possibility of denitrification as a cause of nitrogen loss, while anaerobic ammonia oxidation assays resulted in rates consistent with anammox stoichiometry. Full-scale experiments at different dissolved oxygen (DO = 2 - 6 mg/L) setpoints indicated persistent nitrogen loss that was partly sensitive to DO concentrations. Genome-resolved metagenomics confirmed the high abundance (relative abundance 6.53 ± 0.34%) of two Brocadia-like anammox populations, while comammox bacteria within the Ca. Nitrospira nitrosa cluster were lower in abundance (0.37 ± 0.03%) and Nitrosomonas-like ammonia oxidizers were even lower (0.12 ± 0.02%). Collectively, our study reports for the first time the co-occurrence and cooperation of comammox and anammox bacteria in a full-scale municipal wastewater treatment system.

Whole-Genome Sequence of Endophytic Bacteria Associated with Poison Ivy Vine (Toxicodendron radicans)

Here, we report the genome assemblies of 11 endophytic bacteria, isolated from poison ivy vine (Toxicodendron radicans). Five species belonging to the genus Pseudomonas, two species of Curtobacterium, one strain of Pantoea agglomerans, and one species from the Bacillus, Cellulomonas, and Enterobacter genera were isolated from the interior tissue of poison ivy.

Bacterial aerobic methane cycling by the marine sponge-associated microbiome

BACKGROUND

Methanotrophy by the sponge-hosted microbiome has been mainly reported in the ecological context of deep-sea hydrocarbon seep niches where methane is either produced geothermically or via anaerobic methanogenic archaea inhabiting the sulfate-depleted sediments. However, methane-oxidizing bacteria from the candidate phylum Binatota have recently been described and shown to be present in oxic shallow-water marine sponges, where sources of methane remain undescribed.

RESULTS

Here, using an integrative -omics approach, we provide evidence for sponge-hosted bacterial methane synthesis occurring in fully oxygenated shallow-water habitats. Specifically, we suggest methane generation occurs via at least two independent pathways involving methylamine and methylphosphonate transformations that, concomitantly to aerobic methane production, generate bioavailable nitrogen and phosphate, respectively. Methylphosphonate may be sourced from seawater continuously filtered by the sponge host. Methylamines may also be externally sourced or, alternatively, generated by a multi-step metabolic process where carnitine, derived from sponge cell debris, is transformed to methylamine by different sponge-hosted microbial lineages. Finally, methanotrophs specialized in pigment production, affiliated to the phylum Binatota, may provide a photoprotective function, closing a previously undescribed C1-metabolic loop that involves both the sponge host and specific members of the associated microbial community.

CONCLUSION

Given the global distribution of this ancient animal lineage and their remarkable water filtration activity, sponge-hosted methane cycling may affect methane supersaturation in oxic coastal environments. Depending on the net balance between methane production and consumption, sponges may serve as marine sources or sinks of this potent greenhouse gas. Video Abstract.

Intracytoplasmic-membrane development in alphaproteobacteria involves the homolog of the mitochondrial crista-developing protein Mic60

Mitochondrial cristae expand the surface area of respiratory membranes and ultimately allow for the evolutionary scaling of respiration with cell volume across eukaryotes. The discovery of Mic60 homologs among alphaproteobacteria, the closest extant relatives of mitochondria, suggested that cristae might have evolved from bacterial intracytoplasmic membranes (ICMs). Here, we investigated the predicted structure and function of alphaproteobacterial Mic60, and a protein encoded by an adjacent gene Orf52, in two distantly related purple alphaproteobacteria, Rhodobacter sphaeroides and Rhodopseudomonas palustris. In addition, we assessed the potential physical interactors of Mic60 and Orf52 in R. sphaeroides. We show that the three α helices of mitochondrial Mic60's mitofilin domain, as well as its adjacent membrane-binding amphipathic helix, are present in alphaproteobacterial Mic60. The disruption of Mic60 and Orf52 caused photoheterotrophic growth defects, which are most severe under low light conditions, and both their disruption and overexpression led to enlarged ICMs in both studied alphaproteobacteria. We also found that alphaproteobacterial Mic60 physically interacts with BamA, the homolog of Sam50, one of the main physical interactors of eukaryotic Mic60. This interaction, responsible for making contact sites at mitochondrial envelopes, has been conserved in modern alphaproteobacteria despite more than a billion years of evolutionary divergence. Our results suggest a role for Mic60 in photosynthetic ICM development and contact site formation at alphaproteobacterial envelopes. Overall, we provide support for the hypothesis that mitochondrial cristae evolved from alphaproteobacterial ICMs and have therefore improved our understanding of the nature of the mitochondrial ancestor.

Genomic Characterization and Wetland Occurrence of a Novel Campylobacter Isolate from Canada Geese

Populations of resident, non-migratory Canada geese are rapidly increasing. Canada geese are known to transmit viral and bacterial diseases, posing a possible threat to human health. The most prevalent pathogens vectored by geese are Campylobacter species, yet the current understanding of the identity and virulence of these pathogens is limited. In our previous study, we observed a high prevalence of Campylobacter spp. in the Banklick Creek wetland—a constructed treatment wetland (CTW) located in northern KY (USA) used to understand sources of fecal contamination originating from humans and waterfowl frequenting the area. To identify the types of Campylobacter spp. found contaminating the CTW, we performed genetic analyses of Campylobacter 16s ribosomal RNA amplified from CTW water samples and collected fecal material from birds frequenting those areas. Our results showed a high occurrence of a Campylobacter canadensis-like clade from the sampling sites. Whole-genome sequence analyses of an isolate from Canada goose fecal material, called MG1, were used to confirm the identity of the CTW isolates. Further, we examined the phylogenomic position, virulence gene content, and antimicrobial resistance gene profile of MG1. Lastly, we developed an MG1-specific real-time PCR assay and confirmed the presence of MG1 in Canada goose fecal samples surrounding the CTW. Our findings reveal that the Canada goose-vectored Campylobacter sp. MG1 is a novel isolate compared to C. canadensis that possesses possible zoonotic potential, which may be of human health concern.

Nitrogen-metabolising microorganism analysis in rapid sand filters from drinking water treatment plant

Sand filters (SFs) are common treatment processes for nitrogen pollutant removal in drinking water treatment plants (DWTPs). However, the mechanisms on the nitrogen-cycling role of SFs are still unclear. In this study, 16S rRNA gene amplicon sequencing was used to characterise the diversity and composition of the bacterial community in SFs from DWTPs. Additionally, metagenomics approach was used to determine the functional microorganisms involved in nitrogen cycle in SFs. Our results showed that Pseudomonadota, Acidobacteria, Nitrospirae and Chloroflexi dominated in SFs. Subsequently, 85 high-quality metagenome-assembled genomes (MAGs) were retrieved from metagenome datasets of selected SFs involving nitrification, assimilatory nitrogen reduction, denitrification and anaerobic ammonia oxidation (anammox) processes. Read mapping to reference genomes of Nitrospira and the phylogenetic tree of the ammonia monooxygenase subunit A gene, amoA, suggested that Nitrospira is abundantly found in SFs. Furthermore, according to their genetic content, a nitrogen metabolic model in SFs was proposed using representative MAGs and pure culture isolate. Quantitative real-time polymerase chain reaction (qPCR) showed that ammonia-oxidising bacteria (AOB) and archaea (AOA), and complete ammonia oxidisers (comammox) were ubiquitous in the SFs, with the abundance of comammox being higher than that of AOA and AOB. Moreover, we identified a bacterial strain with a high NO₃-N removal rate as Pseudomonas sp. DW-5, which could be applied in the bioremediation of micro-polluted drinking water sources. Our study provides insights into functional nitrogen-metabolising microbes in SFs of DWTPs.

CrAssphage May Be Viable Markers of Contamination in Pristine and Contaminated River Water

Viruses are the most biologically abundant entities and may be ideal indicators of fecal pollutants in water. Anthropogenic activities have triggered drastic ecosystem changes in rivers, leading to substantial shifts in chemical and biological attributes. Here, we evaluate the viability of using the presence of crAssphage as indicators of fecal contamination in South African rivers. Shotgun analysis revealed diverse crAssphage viruses in these rivers, which are impacted by chemical and biological pollution. Overall, the diversity and relative abundances of these viruses was higher in contaminated sites compared to pristine locations. In contrast to fecal coliform counts, crAssphage sequences were detected in pristine rivers, supporting the assertion that the afore mentioned marker may be a more accurate indicator of fecal contamination. Our data demonstrate the presence of diverse putative hosts which includes members of the phyla Bacteroidota, Pseudomonadota, Verrucomicrobiota, and Bacillota. Phylogenetic analysis revealed novel subfamilies, suggesting that rivers potentially harbor distinct and uncharacterized clades of crAssphage. These data provide the first insights regarding the diversity, distribution, and functional roles of crAssphage in rivers. Taken together, the results support the potential application of crAssphage as viable markers for water quality monitoring. IMPORTANCE Rivers support substantial populations and provide important ecosystem services. Despite the application of fecal coliform tests and other markers, we lack rapid and reproducible approaches for determining fecal contamination in rivers. Waterborne viral outbreaks have been reported even after fecal indicator bacteria (FIB) were suggested to be absent or below regulated levels of coliforms. This indicates a need to develop and apply improved indicators of pollutants in aquatic ecosystems. Here, we evaluate the viability of crAssphage as indicators of fecal contamination in two South African rivers. We assess the abundance, distribution, and diversity of these viruses in sites that had been predicted pristine or contaminated by FIB analysis. We show that crAssphage are ideal and sensitive markers for fecal contamination and describe novel clades of crAss-like phages. Known crAss-like subfamilies were unrepresented in our data, suggesting that the diversity of these viruses may reflect geographic locality and dependence.

Importance of mobile genetic element immunity in numerically abundant Trichodesmium clades

The colony-forming cyanobacteria Trichodesmium spp. are considered one of the most important nitrogen-fixing genera in the warm, low nutrient ocean. Despite this central biogeochemical role, many questions about their evolution, physiology, and trophic interactions remain unanswered. To address these questions, we describe Trichodesmium pangenomic potential via significantly improved genomic assemblies from two isolates and 15 new >50% complete Trichodesmium metagenome-assembled genomes from hand-picked, Trichodesmium colonies spanning the Atlantic Ocean. Phylogenomics identified ~four N2 fixing clades of Trichodesmium across the transect, with T. thiebautii dominating the colony-specific reads. Pangenomic analyses showed that all T. thiebautii MAGs are enriched in COG defense mechanisms and encode a vertically inherited Type III-B Clustered Regularly Interspaced Short Palindromic Repeats and associated protein-based immunity system (CRISPR-Cas). Surprisingly, this CRISPR-Cas system was absent in all T. erythraeum genomes, vertically inherited by T. thiebautii, and correlated with increased signatures of horizontal gene transfer. Additionally, the system was expressed in metaproteomic and transcriptomic datasets and CRISPR spacer sequences with 100% identical hits to field-assembled, putative phage genome fragments were identified. While the currently CO2-limited T. erythraeum is expected to be a 'winner' of anthropogenic climate change, their genomic dearth of known phage resistance mechanisms, compared to T. thiebautii, could put this outcome in question. Thus, the clear demarcation of T. thiebautii maintaining CRISPR-Cas systems, while T. erythraeum does not, identifies Trichodesmium as an ecologically important CRISPR-Cas model system, and highlights the need for more research on phage-Trichodesmium interactions.

Seasonal activities of the phyllosphere microbiome of perennial crops

Understanding the interactions between plants and microorganisms can inform microbiome management to enhance crop productivity and resilience to stress. Here, we apply a genome-centric approach to identify ecologically important leaf microbiome members on replicated plots of field-grown switchgrass and miscanthus, and to quantify their activities over two growing seasons for switchgrass. We use metagenome and metatranscriptome sequencing and curate 40 medium- and high-quality metagenome-assembled-genomes (MAGs). We find that classes represented by these MAGs (Actinomycetia, Alpha- and Gamma- Proteobacteria, and Bacteroidota) are active in the late season, and upregulate transcripts for short-chain dehydrogenase, molybdopterin oxidoreductase, and polyketide cyclase. Stress-associated pathways are expressed for most MAGs, suggesting engagement with the host environment. We also detect seasonally activated biosynthetic pathways for terpenes and various non-ribosomal peptide pathways that are poorly annotated. Our findings support that leaf-associated bacterial populations are seasonally dynamic and responsive to host cues.

Co-diversification of an intestinal Mycoplasma and its salmonid host

Understanding the evolutionary relationships between a host and its intestinal resident bacteria can transform how we understand adaptive phenotypic traits. The interplay between hosts and their resident bacteria inevitably affects the intestinal environment and, thereby, the living conditions of both the host and the microbiota. Thereby this co-existence likely influences the fitness of both bacteria and host. Whether this co-existence leads to evolutionary co-diversification in animals is largely unexplored, mainly due to the complexity of the environment and microbial communities and the often low host selection. We present the gut metagenome from wild Atlantic salmon (Salmo salar), a new wild organism model with an intestinal microbiota of low complexity and a well-described population structure, making it well-suited for investigating co-evolution. Our data reveal a strong host selection of a core gut microbiota dominated by a single Mycoplasma species. We found a clear co-diversification between the population structure of Atlantic salmon and nucleotide variability of the intestinal Mycoplasma populations conforming to expectations from co-evolution between host and resident bacteria. Our results show that the stable microbiota of Atlantic salmon has evolved with its salmonid host populations while potentially providing adaptive traits to the salmon host populations, including defence mechanisms, biosynthesis of essential amino acids, and metabolism of B vitamins. We highlight Atlantic salmon as a novel model for studying co-evolution between vertebrate hosts and their resident bacteria.

Site Specialization of Human Oral Veillonella Species

Veillonella species are abundant members of the human oral microbiome with multiple interspecies commensal relationships. Examining the distribution patterns of Veillonella species across the oral cavity is fundamental to understanding their oral ecology. In this study, we used a combination of pangenomic analysis and oral metagenomic information to clarify Veillonella taxonomy and to test the site specialist hypothesis for the Veillonella genus, which contends that most oral bacterial species are adapted to live at specific oral sites. Using isolate genome sequences combined with shotgun metagenomic sequence data, we showed that Veillonella species have clear, differential site specificity: Veillonella parvula showed strong preference for supra- and subgingival plaque, while closely related V. dispar, as well as more distantly related V. atypica, preferred the tongue dorsum, tonsils, throat, and hard palate. In addition, the provisionally named Veillonella sp. Human Microbial Taxon 780 showed strong site specificity for keratinized gingiva. Using comparative genomic analysis, we identified genes associated with thiamine biosynthesis and the reductive pentose phosphate cycle that may enable Veillonella species to occupy their respective habitats. IMPORTANCE Understanding the microbial ecology of the mouth is fundamental for understanding human physiology. In this study, metapangenomics demonstrated that different Veillonella species have clear ecological preferences in the oral cavity of healthy humans, validating the site specialist hypothesis. Furthermore, the gene pool of different Veillonella species was found to be reflective of their ecology, illuminating the potential role of vitamins and carbohydrates in determining Veillonella distribution patterns and interspecies interactions.

Comparative Genomic and Metagenomic Investigations of the Corynebacterium tuberculostearicum Species Complex Reveals Potential Mechanisms Underlying Associations To Skin Health and Disease

Corynebacterium are a diverse genus and dominant member of the human skin microbiome. Recently, we reported that the most prevalent Corynebacterium species found on skin, including Corynebacterium tuberculostearicum and Corynebacterium kefirresidentii, comprise a narrow species complex despite the diversity of the genus. Here, we apply high-resolution phylogenomics and comparative genomics to describe the structure of the C. tuberculostearicum species complex and highlight genetic traits which are enriched or depleted in it relative to other Corynebacterium. Through metagenomic investigations, we also find that individual species within the complex can associate with specific body sites. Finally, we discover that one species from the complex, C. kefirresidentii, increases in relative abundance during atopic dermatitis flares, and show that most genomes of this species encode a colocalized set of putative virulence genes. IMPORTANCE Corynebacterium are commonly found bacteria on the human skin. In this study, we perform comparative genomics to gain insight into genetic traits which differentiate a phylogenetically related group of Corynebacterium, the Corynebacterium tuberculostearicum species complex, that includes the most prevalent species from the genus in skin microbiomes. After resolving the presence of distinct species within the complex, we applied metagenomic analysis to uncover biogeographic associations of individual species within the complex with specific body sites and discovered that one species, commonly found in the nares of individuals, increases in abundance across multiple body sites during atopic dermatitis flares.

Whole Genome-Based Characterization of Multidrug Resistant Enterobacter and Klebsiella aerogenes Isolates from Lebanon

Enterobacter spp. and Klebsiella aerogenes are rod-shaped Gram-negative opportunistic pathogens. This study aimed at the molecular and genomic characterization of multidrug resistant Enterobacter spp. and K. aerogenes isolates recovered from hospitalized patients in a tertiary care hospital in Lebanon. A total of 59 Enterobacter spp. clinical isolates consisting of 41 carbapenem-resistant and 18 susceptible by Etest were included in this study. Genotypic identification through whole-genome sequencing (WGS) was performed and confirmed in silico. Resistance and plasmid profiles were studied using ResFinder4.0 and Plasmid-Finder2.1. Multilocus sequence typing (MLST) was used to determine the isolates' clonality. Using the average nucleotide identity (ANI) we identified and confirmed that 47 (80%) isolates were E. hormaechei, 11 (18%) were Klebsiella aerogenes and 1 (2%) was an E. cloacae. Carbapenem-resistance was detected among 41 isolates all showing an MIC₉₀ of ≥ 32 μg/mL for ertapenem, imipenem, and meropenem. bla_NDM-1 (58.5%), bla_ACT-16 (54%), and bla_OXA-1 (54%) were the most common detected β-lactamases, while bla_CTX-M-15 (68%) was the main detected extended-spectrum β-lactamase (ESBL) encoding gene. Chromosomal ampC, carbapenemase encoding genes, and porin modifications were among the detected carbapenem resistance determinants. The carbapenemase encoding genes were linked to three well-defined plasmid Inc groups, IncFII/IncFIB, IncX3, and IncL. MLST typing revealed the diversity within the studied isolates, with ST114 being the most common among the studied E. hormaechei.: The spread of carbapenem-resistant isolates in clinical settings in Lebanon is a serious challenge. Screening and continuous monitoring through WGS analysis could effectively limit the dissemination of drug-resistant isolates in hospitalized patients. IMPORTANCE Drug resistance is an increasing global public health threat that involves most disease-causing organisms and antimicrobial drugs. Drug-resistant organisms spread in health care settings, and resistance to multiple drugs is common. Our study demonstrated the mechanisms leading to resistance against the last resort antimicrobial agents among members of the Enterobacteriaceae family. The spread of carbapenem-resistant bacteria in clinical settings is a serious challenge. Screening and continuous monitoring could effectively limit the dissemination of drug-resistant isolates in hospitalized patients.

Phocaeicola oris sp. nov., an anaerobic bacterium isolated from the saliva of a patient with oral squamous cell carcinoma

A Gram-stain-negative or -positive, strictly anaerobic, non-spore-forming and pleomorphic bacterium (designated 14-104^T) was isolated from the saliva sample of a patient with oral squamous cell carcinoma. It was an acid-tolerant neutralophilic mesophile, growing at between 20 and 40 °C (with optimum growth at 30 °C) and pH between pH 3.0 and 7.0 (with optimum growth at pH 6.0-7.0). It contained anteiso-C_15 : 0 and C_15 : 0 as the major fatty acids. The genome size of strain 14-104^T was 2.98 Mbp, and the G+C content was 39.6 mol%. It shared <87 % 16S rRNA sequence similarity, <71 % orthologous average nucleotide identity, <76 % average amino acid identity and <68 %% of conserved proteins with its closest relative, Phocaeicola abscessus CCUG 55929^T. Reconstruction of phylogenetic and phylogenomic trees revealed that strain 14-104^T and P. abscessus CCUG 55929^T were clustered as a distinct clade without any other terminal node. The phylogenetic and phylogenomic analyses along with physiological and chemotaxonomic data indicated that strain 14-104^T represents a novel species in the genus Phocaeicola, for which the name Phocaeicola oris sp. nov. is proposed. The type strain is 14-104^T (=BCRC 81305^T= NBRC 115041^T).

Genome-centered metagenomics illuminates adaptations of core members to a partial Nitritation-Anammox bioreactor under periodic microaeration

The partial nitritation-anaerobic ammonium oxidation (anammox; PN-A) process has been considered a sustainable method for wastewater ammonium removal, with recent attempts to treat low-strength wastewater. However, how microbes adapt to the alternate microaerobic-anoxic operation of the process when treating low ammonium concentrations remains poorly understood. In this study, we applied a metagenomic approach to determine the genomic contents of core members in a PN-A reactor treating inorganic ammonium wastewater at loading as low as 0.0192 kg-N/m³/day. The metabolic traits of metagenome-assembled genomes from 18 core species were analyzed. Taxonomically diverse ammonia oxidizers, including two Nitrosomonas species, a comammox Nitrospira species, a novel Chloroflexota-related species, and two anammox bacteria, Ca. Brocadia and Ca. Jettenia, accounted for the PN-A reactions. The characteristics of a series of genes encoding class II ribonucleotide reductase, high-affinity bd-type terminal oxidase, and diverse antioxidant enzymes revealed that comammox Nitrospira has a superior adaptation ability over the competitors, which may confer the privileged partnership with anammox bacteria in the PN-A reactor. This finding is supported by the long-term monitoring experiment, showing the predominance of the comammox Nitrospira in the ammonia-oxidizing community. Metagenomic analysis of seven heterotrophs suggested that nitrate reduction is a common capability in potentially using endogenous carbohydrates and peptides to enhance nitrogen removals. The prevalence of class II ribonucleotide reductase and antioxidant enzymes genes may grant the adaptation to cyclically microaerobic/anoxic environments. The predominant heterotroph is affiliated with Chloroflexota; its genome encodes complete pathways for synthesizing vitamin B6 and methionine. By contrast, other than the two growth factors, Nitrospira and anammox bacteria are complementary to produce various vitamins and amino acids. Besides, the novel Chloroflexota-related ammonia oxidizer lacks corresponding genes for detoxifying the reactive oxygen species and thus requires the aid of co-existing members to alleviate oxidative stress. The analysis results forecast the exchanges of substrates and nutrients as well as the collective alleviation of oxidative stress among the core populations. The new findings of the genomic features and predicted microbial interplay shed light on microbial adaptation to intermittent microaeration specific to the PN-A reactor, which may aid in improving its application to low-strength ammonium wastewater.

Novel and unusual genes for nitrogen and metal cycling in Planctomycetota- and KSB1-affiliated metagenome-assembled genomes reconstructed from a marine subsea tunnel

The Oslofjord subsea road tunnel is a unique environment in which the typically anoxic marine deep subsurface is exposed to oxygen. Concrete biodeterioration and steel corrosion in the tunnel have been linked to the growth of iron- and manganese-oxidizing biofilms in areas of saline water seepage. Surprisingly, previous 16S rRNA gene surveys of biofilm samples revealed microbial communities dominated by sequences affiliated with nitrogen-cycling microorganisms. This study aimed to identify microbial genomes with metabolic potential for novel nitrogen- and metal-cycling reactions, representing biofilm microorganisms that could link these cycles and play a role in concrete biodeterioration. We reconstructed 33 abundant, novel metagenome-assembled genomes (MAGs) affiliated with the phylum Planctomycetota and the candidate phylum KSB1. We identified novel and unusual genes and gene clusters in these MAGs related to anaerobic ammonium oxidation, nitrite oxidation, and other nitrogen-cycling reactions. Additionally, 26 of 33 MAGs also had the potential for iron, manganese, and arsenite cycling, suggesting that bacteria represented by these genomes might couple these reactions. Our results expand the diversity of microorganisms putatively involved in nitrogen and metal cycling, and contribute to our understanding of potential biofilm impacts on built infrastructure.

The ESKAPE mobilome contributes to the spread of antimicrobial resistance and CRISPR-mediated conflict between mobile genetic elements

Mobile genetic elements (MGEs) mediate the shuffling of genes among organisms. They contribute to the spread of virulence and antibiotic resistance (AMR) genes in human pathogens, such as the particularly problematic group of ESKAPE pathogens. Here, we performed the first systematic analysis of MGEs, including plasmids, prophages, and integrative and conjugative/mobilizable elements (ICEs/IMEs), across all ESKAPE pathogens. We found that different MGE types are asymmetrically distributed across these pathogens, and that most horizontal gene transfer (HGT) events are restricted by phylum or genus. We show that the MGEs proteome is involved in diverse functional processes and distinguish widespread proteins within the ESKAPE context. Moreover, anti-CRISPRs and AMR genes are overrepresented in the ESKAPE mobilome. Our results also underscore species-specific trends shaping the number of MGEs, AMR, and virulence genes across pairs of conspecific ESKAPE genomes with and without CRISPR-Cas systems. Finally, we observed that CRISPR spacers found on prophages, ICEs/IMEs, and plasmids have different targeting biases: while plasmid and prophage CRISPRs almost exclusively target other plasmids and prophages, respectively, ICEs/IMEs CRISPRs preferentially target prophages. Overall, our study highlights the general importance of the ESKAPE mobilome in contributing to the spread of AMR and mediating conflict among MGEs.

Genetic and physiological insights into the diazotrophic activity of a non-cyanobacterial marine diazotroph

Nitrogen (N₂ ) fixation, or diazotrophy, supports a large part of primary production in oceans. Culture-independent approaches highlighted the presence in abundance of marine non-cyanobacterial diazotrophs (NCD), but their ecophysiology remains elusive, mostly because of the low number of isolated NCD and because of the lack of available genetic tools for these isolates. Here, a dual genetic and functional approach allowed unveiling the ecophysiology of a marine NCD affiliated to the species Vibrio diazotrophicus. Physiological characterization of the first marine NCD mutant obtained so far was performed using a soft-gellan assay, demonstrating that a ΔnifH mutant is not able to grow in nitrogen-free media. Furthermore, we demonstrated that V. diazotrophicus produces a thick biofilm under diazotrophic conditions, suggesting biofilm production as an adaptive response of this NCD to cope with the inhibition of nitrogen fixation by molecular oxygen. Finally, the genomic signature of V. diazotrophicus is essentially absent from metagenomic data of Tara Ocean expeditions, despite having been isolated from various marine environments. We think that the genetically tractable V. diazotrophicus strain used in this study may serve as an ideal model to study the ecophysiology of these overlooked procaryotic group.

Adaptations of endolithic communities to abrupt environmental changes in a hyper-arid desert

The adaptation mechanisms of microbial communities to natural perturbations remain unexplored, particularly in extreme environments. The extremophilic communities of halite (NaCl) nodules from the hyper-arid core of the Atacama Desert are self-sustained and represent a unique opportunity to study functional adaptations and community dynamics with changing environmental conditions. We transplanted halite nodules to different sites in the desert and investigated how their taxonomic, cellular, and biochemical changes correlated with water availability, using environmental data modeling and metagenomic analyses. Salt-in strategists, mainly represented by haloarchaea, significantly increased in relative abundance at sites characterized by extreme dryness, multiple wet/dry cycles, and colder conditions. The functional analysis of metagenome-assembled genomes (MAGs) revealed site-specific enrichments in archaeal MAGs encoding for the uptake of various compatible solutes and for glycerol utilization. These findings suggest that opportunistic salt-in strategists took over the halite communities at the driest sites. They most likely benefited from compounds newly released in the environment by the death of microorganisms least adapted to the new conditions. The observed changes were consistent with the need to maximize cellular bioenergetics when confronted with lower water availability and higher salinity, providing valuable information on microbial community adaptations and resilience to climate change.

Hyperthermophilic endospores germinate and metabolize organic carbon in sediments heated to 80°C

Cold surface sediments host a seedbank of functionally diverse thermophilic bacteria. These thermophiles are present as endospores, which are widely dispersed in aquatic environments. Here, we investigated the functional potential of endospore populations in cold surface sediments heated to 80°C. Microbial production of acetate was observed at 80°C and could be enhanced by supplying additional organic carbon substrates. Comparison of 16S rRNA gene amplicon libraries from 80°C enrichments to sediments heated to lower temperatures (50-70°C) showed that temperature selects for distinct populations of endospore-forming bacteria. Whereas sulfate-reducing thermophiles were enriched in 50-70°C incubations, 80°C exceeds their thermal tolerance and selects for hyperthermophilic organotrophic bacteria that are similarly detected in amplicon libraries from sediments heated to 90°C. Genome-resolved metagenomics revealed novel carbon cycling members of Symbiobacteriales, Thermosediminibacteraceae, Thermanaeromonas and Calditerricola with the genomic potential for the degradation of carbohydrates, sugars, amino acids and nucleotides. Endospores of thermophilic bacteria are deposited on seabed sediments worldwide where they remain dormant as they are buried in the accumulating sediments. Our results suggest that endospore populations could be activated by temperature increases encountered during burial and show the potential for organotrophic metabolic activity contributing to acetate generation in deep hot sediments.

Phylogenetic and functional diverse ANME-1 thrive in Arctic hydrothermal vents

The methane-rich areas, the Loki's Castle vent field and the Jan Mayen vent field at the Arctic Mid Ocean Ridge (AMOR), host abundant niches for anaerobic methane-oxidizers, which are predominantly filled by members of the ANME-1. In this study, we used a metagenomic-based approach that revealed the presence of phylogenetic and functional different ANME-1 subgroups at AMOR, with heterogeneous distribution. Based on a common analysis of ANME-1 genomes from AMOR and other geographic locations, we observed that AMOR subgroups clustered with a vent-specific ANME-1 group that occurs solely at vents, and with a generalist ANME-1 group, with a mixed environmental origin. Generalist ANME-1 are enriched in genes coding for stress response and defense strategies, suggesting functional diversity among AMOR subgroups. ANME-1 encode a conserved energy metabolism, indicating strong adaptation to sulfate-methane-rich sediments in marine systems, which does not however prevent global dispersion. A deep branching family named Ca. Veteromethanophagaceae was identified. The basal position of vent-related ANME-1 in phylogenomic trees suggests that ANME-1 originated at hydrothermal vents. The heterogeneous and variable physicochemical conditions present in diffuse venting areas of hydrothermal fields could have favored the diversification of ANME-1 into lineages that can tolerate geochemical and environmental variations.

VEBA: a modular end-to-end suite for in silico recovery, clustering, and analysis of prokaryotic, microeukaryotic, and viral genomes from metagenomes

BACKGROUND

With the advent of metagenomics, the importance of microorganisms and how their interactions are relevant to ecosystem resilience, sustainability, and human health has become evident. Cataloging and preserving biodiversity is paramount not only for the Earth's natural systems but also for discovering solutions to challenges that we face as a growing civilization. Metagenomics pertains to the in silico study of all microorganisms within an ecological community in situ, however, many software suites recover only prokaryotes and have limited to no support for viruses and eukaryotes.

RESULTS

In this study, we introduce the Viral Eukaryotic Bacterial Archaeal (VEBA) open-source software suite developed to recover genomes from all domains. To our knowledge, VEBA is the first end-to-end metagenomics suite that can directly recover, quality assess, and classify prokaryotic, eukaryotic, and viral genomes from metagenomes. VEBA implements a novel iterative binning procedure and hybrid sample-specific/multi-sample framework that yields more genomes than any existing methodology alone. VEBA includes a consensus microeukaryotic database containing proteins from existing databases to optimize microeukaryotic gene modeling and taxonomic classification. VEBA also provides a unique clustering-based dereplication strategy allowing for sample-specific genomes and genes to be directly compared across non-overlapping biological samples. Finally, VEBA is the only pipeline that automates the detection of candidate phyla radiation bacteria and implements the appropriate genome quality assessments. VEBA's capabilities are demonstrated by reanalyzing 3 existing public datasets which recovered a total of 948 MAGs (458 prokaryotic, 8 eukaryotic, and 482 viral) including several uncharacterized organisms and organisms with no public genome representatives.

CONCLUSIONS

The VEBA software suite allows for the in silico recovery of microorganisms from all domains of life by integrating cutting edge algorithms in novel ways. VEBA fully integrates both end-to-end and task-specific metagenomic analysis in a modular architecture that minimizes dependencies and maximizes productivity. The contributions of VEBA to the metagenomics community includes seamless end-to-end metagenomics analysis but also provides users with the flexibility to perform specific analytical tasks. VEBA allows for the automation of several metagenomics steps and shows that new information can be recovered from existing datasets.

Oral microbiome diversity: The curious case of Corynebacterium sp. isolation

Oral microbiome sequencing efforts revealed the presence of hundreds of different microbes. Interindividual differences at strain and species resolution suggest that microbiome diversity could lead to mechanistically distinct gene regulation as well as species-related differences in phenotypes. Commonly, gene regulation and related phenotypes are studied in a few selected strains of a particular species with conclusions that are mostly generalized. The aim of this study was to isolate several species of Corynebacterium using an established protocol that led to the previous isolation of C. durum. Characterization of C. durum interspecies interactions revealed a specific mechanism for chain elongation in Streptococcus sanguinis that was the result of corynebacterial fatty acid production and secretion. While the protocol was successfully applied to isolate what we presumed to be additional Corynebacterium based on several phenotypic traits that seem to be identical to C. durum, genome sequencing of the newly isolated strains placed them closer to Actinomyces. Both Corynebacterium and Actinomyces are suborders of the Actinobacteridae and related species. Our study suggests to take several comprehensive strategies into consideration when taxonomically identifying closely related microorganisms. Furthermore, it seems to be important to test common core phenotypes in bacterial ecology to understand the behavior of specific groups of microbes, rather than simply relying upon genome sequence homology to establish relationships in the microbiome.

Integrating phylogenetic and functional data in microbiome studies

MOTIVATION

Microbiome functional data are frequently analyzed to identify associations between microbial functions (e.g. genes) and sample groups of interest. However, it is challenging to distinguish between different possible explanations for variation in community-wide functional profiles by considering functions alone. To help address this problem, we have developed POMS, a package that implements multiple phylogeny-aware frameworks to more robustly identify enriched functions.

RESULTS

The key contribution is an extended balance-tree workflow that incorporates functional and taxonomic information to identify functions that are consistently enriched in sample groups across independent taxonomic lineages. Our package also includes a workflow for running phylogenetic regression. Based on simulated data we demonstrate that these approaches more accurately identify gene families that confer a selective advantage compared with commonly used tools. We also show that POMS in particular can identify enriched functions in real-world metagenomics datasets that are potential targets of strong selection on multiple members of the microbiome.

AVAILABILITY AND IMPLEMENTATION

These workflows are freely available in the POMS R package at https://github.com/gavinmdouglas/POMS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparative Genomic Analysis of Agarolytic Flavobacterium faecale WV33T

Flavobacteria are widely dispersed in a variety of environments and produce various polysaccharide-degrading enzymes. Here, we report the complete genome of Flavobacterium faecale WV33^T, an agar-degrading bacterium isolated from the stools of Antarctic penguins. The sequenced genome of F. faecale WV33^T represents a single circular chromosome (4,621,116 bp, 35.2% G + C content), containing 3984 coding DNA sequences and 85 RNA-coding genes. The genome of F. faecale WV33^T contains 154 genes that encode carbohydrate-active enzymes (CAZymes). Among the CAZymes, seven putative genes encoding agarases have been identified in the genome. Transcriptional analysis revealed that the expression of these putative agarases was significantly enhanced by the presence of agar in the culture medium, suggesting that these proteins are involved in agar hydrolysis. Pangenome analysis revealed that the genomes of the 27 Flavobacterium type strains, including F. faecale WV33^T, tend to be very plastic, and Flavobacterium strains are unique species with a tiny core genome and a large non-core region. The average nucleotide identity and phylogenomic analysis of the 27 Flavobacterium-type strains showed that F. faecale WV33^T was positioned in a unique clade in the evolutionary tree.

Convergent evolution of a labile nutritional symbiosis in ants

Ants are among the most successful organisms on Earth. It has been suggested that forming symbioses with nutrient-supplementing microbes may have contributed to their success, by allowing ants to invade otherwise inaccessible niches. However, it is unclear whether ants have evolved symbioses repeatedly to overcome the same nutrient limitations. Here, we address this question by comparing the independently evolved symbioses in Camponotus, Plagiolepis, Formica and Cardiocondyla ants. Our analysis reveals the only metabolic function consistently retained in all of the symbiont genomes is the capacity to synthesise tyrosine. We also show that in certain multi-queen lineages that have co-diversified with their symbiont for millions of years, only a fraction of queens carry the symbiont, suggesting ants differ in their colony-level reliance on symbiont-derived resources. Our results imply that symbioses can arise to solve common problems, but hosts may differ in their dependence on symbionts, highlighting the evolutionary forces influencing the persistence of long-term endosymbiotic mutualisms.

Plant Growth-Promoting Activity of Bacteria Isolated from Asian Rice (Oryza sativa L.) Depends on Rice Genotype

Asian rice is one of the most important crops because it is a staple food for almost half of the world's population. To have production of rice keep pace with a growing world population, it is anticipated that the use of fertilizers will also need to increase, which may cause environmental damage through runoff impacts. An alternative strategy to increase crop yield is the use of plant growth-promoting bacteria. Thousands of microbial species can exist in association with plant roots and shoots, and some are critical to the plant's survival. We isolated 140 bacteria from two distantly related rice accessions and investigated whether their impact on the growth of four different rice accessions. The bacterial isolates were screened for their ability to solubilize phosphate, a known plant growth-promoting characteristic, and 25 isolates were selected for further analysis. These 25 phosphate-solubilizing isolates were also able to produce other potentially growth-promoting factors. Five of the most promising bacterial isolates were chosen for whole-genome sequencing. Four of these bacteria, isolates related to Pseudomonas mosselii, a Microvirga sp., Paenibacillus rigui, and Paenibacillus graminis, improved root and shoot growth in a rice genotype-dependent manner. This indicates that while bacteria have several known plant growth-promoting functions, their effects on growth parameters are rice genotype dependent and suggest a close relationship between plants and their microbial partners. IMPORTANCE In this study, endophytic bacterial isolates from roots and shoots of two distantly related rice accessions were characterized phenotypically and genotypically. From the isolated bacterial species, five of the most promising plant growth-promoting bacteria were selected to test their abilities to enhance growth of the four rice accessions. Interestingly, plant growth enhancement was both bacterial isolate specific and plant genotype specific. However, the positive interactions between plant and bacteria could not easily be predicted because rice growth-promoting bacteria isolated from their host plants did not necessarily stimulate growth of their own host.

Metagenome-Assembled Genomes of Bacteria Associated with Massospora cicadina Fungal Plugs from Infected Brood VIII Periodical Cicadas

We report six metagenome-assembled genomes (MAGs) associated with Massospora cicadina strain MCPNR19 (ARSEF 14555), an obligate entomopathogenic fungus of periodical cicadas. The MAGs include representatives of Pantoea, Pseudomonas, Lactococcus, and one potential new Chryseobacterium species. Future research is needed to resolve the ecology of these MAGs and determine whether they represent symbionts or contaminants.

Alkalibacter rhizosphaerae sp. nov., a CO-utilizing bacterium isolated from tidal flat sediment, and emended description of the genus Alkalibacter

A novel anaerobic, rod-shaped, non-motile bacterium, designated strain ES005T, was isolated from tidal flat sediments near the rhizosphere of Phragmites australis at Eulsukdo Island, Republic of Korea. A polyphasic approach revealed that cells of the strain were Gram-stain-positive, catalase- and oxidase-negative, non-spore-forming rods. Phylogenetic analyses based on 16S rRNA gene sequences revealed that strain ES005T belonged to the family Eubacteriaceae , class Clostridia and showed the highest sequence similarity to Alkalibacter mobili s (97.52 %) and followed by Alkalibacter saccharofermentans Z-79820T (96.72%). The OrthoANI value between strain ES005T and A. mobilis was 69.67 %. Strain ES005T grew optimally at 33–37 °C, at pH 6.0–7.0 and in the presence of 1–2 % (w/v) NaCl. Growth in 12.5 % CO atmosphere was observed. Acetate and formate were end products of fructose fermentation and growth on CO. The major cellular fatty acids of strain ES005T were C14 : 0 (39.1 %) and C16 : 0 (26.6 %). The major polar lipids were diphoshatidylgycerol, phosphatidylglycerol and three unidentified phospholipids. The DNA G+C content of strain ES005T was 46.9 mol%. Based on the phenotypic, phylogenetic, genomic and chemotaxonomic features of the isloate, strain ES005T represents a novel species, for which the name Alkalibacter rhizosphaerae sp. nov. is proposed. The type strain is ES005T (=KCTC 25246T=JCM 34530T)

The microbiome of a bacterivorous marine choanoflagellate contains a resource-demanding obligate bacterial associate

Microbial predators such as choanoflagellates are key players in ocean food webs. Choanoflagellates, which are the closest unicellular relatives of animals, consume bacteria and also exhibit marked biological transitions triggered by bacterial compounds, yet their native microbiomes remain uncharacterized. Here we report the discovery of a ubiquitous, uncultured bacterial lineage we name Candidatus Comchoanobacterales ord. nov., related to the human pathogen Coxiella and physically associated with the uncultured marine choanoflagellate Bicosta minor. We analyse complete ‘Comchoano’ genomes acquired after sorting single Bicosta cells, finding signatures of obligate host-dependence, including reduction of pathways encoding glycolysis, membrane components, amino acids and B-vitamins. Comchoano encode the necessary apparatus to import energy and other compounds from the host, proteins for host-cell associations and a type IV secretion system closest to Coxiella’s that is expressed in Pacific Ocean metatranscriptomes. Interactions between choanoflagellates and their microbiota could reshape the direction of energy and resource flow attributed to microbial predators, adding complexity and nuance to marine food webs.

Choanoflagellates are the closest living unicellular relatives of animals and are important bacterivorous predators in the ocean. Here the authors show that the microbiome of this predator includes an obligate, host resource-dependent bacterial associate.

Dethiosulfovibrio faecalis sp. nov., a novel proteolytic, non-sulphur-reducing bacterium isolated from a marine aquaculture solid waste bioreactor

A new Dethiosulfovibrio strain, designated F2BT, was isolated from an anaerobic digester for treating solid waste from a marine recirculating aquaculture system. The motile, Gram-negative, non-spore-forming curved rods were 2–7 µm long and 1 µm in diameter. Growth occurred at temperatures ranging from 20 to 40 °C with a maximum rate of growth at 30 °C. The pH range for growth was pH 6.0–8.0, with a maximum rate of growth at pH 7.5. This isolate was halotolerant growing in NaCl concentrations ranging from 0 to 1.6 M with a maximum rate of growth at 0.4 M. Similarly to the five described Dethiosulfovibrio species, this obligate anaerobe isolate was fermentative, capable of utilizing peptides, amino acids and some organic acids for growth, but unlike described strains in the genus did not reduce thiosulphate or elemental sulphur to hydrogen sulphide during fermentation of organic substrates. The G+C content of 55 mol% is similar to the described Dethiosulfovibrio species. The average nucleotide identity analysis between whole genome sequences showed less than 93.15% sequence similarity between strain F2BT and the five other described Dethiosulfovibrio species. Differences in the physiological and phylogenetic characteristics between the new strain and other Dethiosulfovibrio specied indicate that F2BT represents a novel species of this genus and the epithet Dethiosulfovibrio faecalis sp. nov. is proposed. The type strain is F2BT (=DSM 112078T=KCTC25378T).

Cloacimonadota metabolisms include adaptations in engineered environments that are reflected in the evolutionary history of the phylum

Phylum Cloacimonadota (previously Cloacimonetes, WWE1) is an understudied bacterial lineage frequently associated with engineered and wastewater systems. Cloacimonadota members were abundant and diverse in metagenomic datasets from a municipal landfill, prompting an examination of phylogenetic relationships, metabolic diversity, and pangenomic dynamics across the phylum, based on the 30 publicly available genomes and 24 new metagenome-assembled genomes (MAGs) from landfill samples. We found that Cloacimonadota have distinct evolutionary histories associated with engineered versus natural environments and identified genomic features and metabolic strategies that correlate to habitat of origin. Metabolic reconstructions for MAGs predict an anaerobic, acetogenic, and mixed fermentative and flavin-bifurcation-based anaerobic respiratory lifestyle for the majority of Cloacimonadota surveyed. Genomes from engineered ecosystems encode a suite of genes not typically found in genomes from natural environments including acetate kinase, genes for cysteine degradation to pyruvate, increased diversity of carbon utilization enzymes, and different mechanisms for generating membrane potential and ATP synthesis. This phylum-level examination also clarifies the distribution of functions previously observed for members of the phylum, where propionate oxidation and reverse TCA cycles are not common components of Cloacimonadota metabolism.

Metagenome-Assembled Genomes of Bacterial Symbionts Associated with Insecticide-Resistant and -Susceptible Individuals of the Glassy-Winged Sharpshooter (Homalodisca vitripennis)

The role of microbes in insecticide resistance is an emerging question. Here, we describe six metagenome-assembled genomes (MAGs) associated with the glassy-winged sharpshooter (Homalodisca vitripennis [Germar, 1821]) (Hemiptera, Cicadellidae). MAGs representing the obligate symbionts Candidatus Sulcia muelleri and Candidatus Baumannia cicadellinicola and the facultative symbiont Wolbachia were obtained from imidacloprid-resistant and imidacloprid-susceptible sharpshooters.

Genome of a novel Sediminibacterium discovered in association with two species of freshwater cyanobacteria from streams in Southern California

Here, we report the discovery of a novel Sediminibacterium sequenced from laboratory cultures of freshwater stream cyanobacteria from sites in Southern California, grown in BG11 medium. Our genome-wide analyses reveal a highly contiguous and complete genome (97% BUSCO) that is placed within sediminibacterial clades in phylogenomic analyses. Functional annotation indicates the presence of genes that could be involved in mutualistic/commensal relationship with associated cyanobacterial hosts.

Genomic diversity of genus Limosilactobacillus

The genus Limosilactobacillus (formerly Lactobacillus) contains multiple species considered to be adapted to vertebrates, yet their genomic diversity has not been explored. In this study, we performed comparative genomic analysis of Limosilactobacillus (22 species; 332 genomes) isolated from different niches, further focusing on human strains (11 species; 74 genomes) and their adaptation features to specific body sites. Phylogenomic analysis of Limosilactobacillus showed misidentification of some strains deposited in public databases and existence of putative novel Limosilactobacillus species. The pangenome analysis revealed a remarkable genomic diversity (only 1.3 % of gene clusters are shared), and we did not observe a strong association of the accessory genome with different niches. The pangenome of Limosilactobacillus reuteri and Limosilactobacillus fermentum was open, suggesting that acquisition of genes is still occurring. Although most Limosilactobacillus were predicted as antibiotic susceptible (83%), acquired antibiotic-resistance genes were common in L. reuteri from food-producing animals. Genes related to lactic acid isoform production (>95 %) and putative bacteriocins (70.2%) were identified in most Limosilactobacillus strains, while prophages (55.4%) and CRISPR-Cas systems (32.0%) were less prevalent. Among strains from human sources, several metabolic pathways were predicted as conserved and completed. Their accessory genome was highly variable and did not cluster according to different human body sites, with some exceptions (urogenital Limosilactobacillus vaginalis, Limosilactobacillus portuensis, Limosilactobacillus urinaemulieris and Limosilactobacillus coleohominis or gastrointestinal Limosilactobacillus mucosae). Moreover, we identified 12 Kyoto Encyclopedia of Genes and Genomes (KEGG) orthologues that were significantly enriched in strains from particular body sites. We concluded that evolution of the highly diverse Limosilactobacillus is complex and not always related to niche or human body site origin.

High molecular weight DNA extraction strategies for long-read sequencing of complex metagenomes

By offering extremely long contiguous characterization of individual DNA molecules, rapidly emerging long-read sequencing strategies offer comprehensive insights into the organization of genetic information in genomes and metagenomes. However, successful long-read sequencing experiments demand high concentrations of highly purified DNA of high molecular weight (HMW), which limits the utility of established DNA extraction kits designed for short-read sequencing. The challenges associated with input DNA quality intensify further when working with complex environmental samples of low microbial biomass, which requires new protocols that are tailored to study metagenomes with long-read sequencing. Here, we use human tongue scrapings to benchmark six HMW DNA extraction strategies that are based on commercially available kits, phenol-chloroform (PC) extraction and agarose encasement followed by agarase digestion. A typical end goal of HMW DNA extractions is to obtain the longest possible reads during sequencing, which is often achieved by PC extractions, as demonstrated in sequencing of cultured cells. Yet our analyses that consider overall read-size distribution, assembly performance and the number of circularized elements found in sequencing results suggest that column-based kits with enzyme supplementation, rather than PC methods, may be more appropriate for long-read sequencing of metagenomes.

A closed Candidatus Odinarchaeum chromosome exposes Asgard archaeal viruses

Asgard archaea have recently been identified as the closest archaeal relatives of eukaryotes. Their ecology, and particularly their virome, remain enigmatic. We reassembled and closed the chromosome of Candidatus Odinarchaeum yellowstonii LCB_4, through long-range PCR, revealing CRISPR spacers targeting viral contigs. We found related viruses in the genomes of diverse prokaryotes from geothermal environments, including other Asgard archaea. These viruses open research avenues into the ecology and evolution of Asgard archaea.