GToTree: a user-friendly workflow for phylogenomics

Unveiling microbial guilds and symbiotic relationships in Antarctic sponge microbiomes

Marine sponges host diverse microbial communities. Although we know many of its ecological patterns, a deeper understanding of the polar sponge holobiont is still needed. We combine high-throughput sequencing of ribosomal genes, including the largest taxonomic repertoire of Antarctic sponge species analyzed to date, functional metagenomics, and metagenome-assembled genomes (MAGs). Our findings show that sponges harbor more exclusive bacterial and archaeal communities than seawater, while microbial eukaryotes are mostly shared. Furthermore, bacteria in Antarctic sponge holobionts establish more cooperative interactions than in sponge holobionts from other environments. The bacterial classes that established more positive relations were Bacteroidia, Gamma- and Alphaproteobacteria. Antarctic sponge microbiomes contain microbial guilds that encompass ammonia-oxidizing archaea, ammonia-oxidizing bacteria, nitrite-oxidizing bacteria, and sulfur-oxidizing bacteria. The retrieved MAGs showed a high level of novelty and streamlining signals and belong to the most abundant members of the main microbial guilds in the Antarctic sponge holobiont. Moreover, the genomes of these symbiotic bacteria contain highly abundant functions related to their adaptation to the cold environment, vitamin production, and symbiotic lifestyle, helping the holobiont survive in this extreme environment.

Uncultivated DPANN archaea are ubiquitous inhabitants of global oxygen-deficient zones with diverse metabolic potential

Archaea belonging to the DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota) superphylum have been found in an expanding number of environments and perform a variety of biogeochemical roles, including contributing to carbon, sulfur, and nitrogen cycling. Generally characterized by ultrasmall cell sizes and reduced genomes, DPANN archaea may form mutualistic, commensal, or parasitic interactions with various archaeal and bacterial hosts, influencing the ecology and functioning of microbial communities. While DPANN archaea reportedly comprise a sizeable fraction of the archaeal community within marine oxygen-deficient zone (ODZ) water columns, little is known about their metabolic capabilities in these ecosystems. We report 33 novel metagenome-assembled genomes (MAGs) belonging to the DPANN phyla Nanoarchaeota, Pacearchaeota, Woesearchaeota, Undinarchaeota, Iainarchaeota, and SpSt-1190 from pelagic ODZs in the Eastern Tropical North Pacific and the Arabian Sea. We find these archaea to be permanent, stable residents of all three major ODZs only within anoxic depths, comprising up to 1% of the total microbial community and up to 25%-50% of archaea as estimated from read mapping to MAGs. ODZ DPANN appear to be capable of diverse metabolic functions, including fermentation, organic carbon scavenging, and the cycling of sulfur, hydrogen, and methane. Within a majority of ODZ DPANN, we identify a gene homologous to nitrous oxide reductase. Modeling analyses indicate the feasibility of a nitrous oxide reduction metabolism for host-attached symbionts, and the small genome sizes and reduced metabolic capabilities of most DPANN MAGs suggest host-associated lifestyles within ODZs.

IMPORTANCE

Archaea from the DPANN (Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, and Nanohaloarchaeota) superphylum have diverse metabolic capabilities and participate in multiple biogeochemical cycles. While metagenomics and enrichments have revealed that many DPANN are characterized by ultrasmall genomes, few biosynthetic genes, and episymbiotic lifestyles, much remains unknown about their biology. We report 33 new DPANN metagenome-assembled genomes originating from the three global marine oxygen-deficient zones (ODZs), the first from these regions. We survey DPANN abundance and distribution within the ODZ water column, investigate their biosynthetic capabilities, and report potential roles in the cycling of organic carbon, methane, and nitrogen. We test the hypothesis that nitrous oxide reductases found within several ODZ DPANN genomes may enable ultrasmall episymbionts to serve as nitrous oxide consumers when attached to a host nitrous oxide producer. Our results indicate DPANN archaea as ubiquitous residents within the anoxic core of ODZs with the potential to produce or consume key compounds.

Biofilm colonization and succession in a full-scale partial nitritation-anammox moving bed biofilm reactor

BACKGROUND

Partial nitritation-anammox (PNA) is a biological nitrogen removal process commonly used in wastewater treatment plants for the treatment of warm and nitrogen-rich sludge liquor from anaerobic digestion, often referred to as sidestream wastewater. In these systems, biofilms are frequently used to retain biomass with aerobic ammonia-oxidizing bacteria (AOB) and anammox bacteria, which together convert ammonium to nitrogen gas. Little is known about how these biofilm communities develop, and whether knowledge about the assembly of biofilms in natural communities can be applied to PNA biofilms.

RESULTS

We followed the start-up of a full-scale PNA moving bed biofilm reactor for 175 days using shotgun metagenomics. Environmental filtering likely restricted initial biofilm colonization, resulting in low phylogenetic diversity, with the initial microbial community comprised mainly of Proteobacteria. Facilitative priority effects allowed further biofilm colonization, with the growth of initial aerobic colonizers promoting the arrival and growth of anaerobic taxa like methanogens and anammox bacteria. Among the early colonizers were known 'oligotrophic' ammonia oxidizers including comammox Nitrospira and Nitrosomonas cluster 6a AOB. Increasing the nitrogen load in the bioreactor allowed colonization by 'copiotrophic' Nitrosomonas cluster 7 AOB and resulted in the exclusion of the initial ammonia- and nitrite oxidizers.

CONCLUSIONS

We show that complex dynamic processes occur in PNA microbial communities before a stable bioreactor process is achieved. The results of this study not only contribute to our knowledge about biofilm assembly and PNA bioreactor start-up but could also help guide strategies for the successful implementation of PNA bioreactors. Video Abstract.

Anoxygenic phototroph of the Chloroflexota uses a type I reaction centre

Scientific exploration of phototrophic bacteria over nearly 200 years has revealed large phylogenetic gaps between known phototrophic groups that limit understanding of how phototrophy evolved and diversified1,2. Here, through Boreal Shield lake water incubations, we cultivated an anoxygenic phototrophic bacterium from a previously unknown order within the Chloroflexota phylum that represents a highly novel transition form in the evolution of photosynthesis. Unlike all other known phototrophs, this bacterium uses a type I reaction centre (RCI) for light energy conversion yet belongs to the same bacterial phylum as organisms that use a type II reaction centre (RCII) for phototrophy. Using physiological, phylogenomic and environmental metatranscriptomic data, we demonstrate active RCI-utilizing metabolism by the strain alongside usage of chlorosomes3 and bacteriochlorophylls4 related to those of RCII-utilizing Chloroflexota members. Despite using different reaction centres, our phylogenomic data provide strong evidence that RCI-utilizing and RCII-utilizing Chloroflexia members inherited phototrophy from a most recent common phototrophic ancestor. The Chloroflexota phylum preserves an evolutionary record of the use of contrasting phototrophic modes among genetically related bacteria, giving new context for exploring the diversification of phototrophy on Earth.

Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals

Cycads are known to host symbiotic cyanobacteria, including Nostocales species, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition to Nostoc spp., symbiotic cyanobacteria of the genus Aulosira were identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includes Aulosira spp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.

Natural variation of immune epitopes reveals intrabacterial antagonism

Plants and animals detect biomolecules termed Microbe-Associated Molecular Patterns (MAMPs) and induce immunity. Agricultural production is severely impacted by pathogens which can be controlled by transferring immune receptors. However, most studies use a single MAMP epitope and the impact of diverse multi-copy MAMPs on immune induction is unknown. Here we characterized the epitope landscape from five proteinaceous MAMPs across 4,228 plant-associated bacterial genomes. Despite the diversity sampled, natural variation was constrained and experimentally testable. Immune perception in both Arabidopsis and tomato depended on both epitope sequence and copy number variation. For example, Elongation Factor Tu is predominantly single copy and 92% of its epitopes are immunogenic. Conversely, 99.9% of bacterial genomes contain multiple Cold Shock Proteins and 46% carry a non-immunogenic form. We uncovered a new mechanism for immune evasion, intrabacterial antagonism, where a non-immunogenic Cold Shock Protein blocks perception of immunogenic forms encoded in the same genome. These data will lay the foundation for immune receptor deployment and engineering based on natural variation.

Increased environmental microbial diversity reduces the disease risk of a mosquitocidal pathogen

IMPORTANCE

The host-specific microbiotas of animals can both reduce and increase disease risks from pathogens. In contrast, how environmental microbial communities affect pathogens is largely unexplored. Aquatic habitats are of interest because water enables environmental microbes to readily interact with animal pathogens. Here, we focused on mosquitoes, which are important disease vectors as terrestrial adults but are strictly aquatic as larvae. We identified a pathogen of mosquito larvae from the field as a strain of Chromobacterium haemolyticum. Comparative genomic analyses and functional assays indicate this strain and other Chromobacterium are mosquitocidal but are also opportunistic pathogens of other animals. We also identify a critical role for diversity of the environmental microbiota in disease risk. Our study characterizes both the virulence mechanisms of a pathogen and the role of the environmental microbiota in disease risk to an aquatic animal of significant importance to human health.

Assessing the metabolism, phylogenomic, and taxonomic classification of the halophilic genus Halarchaeum

In this study, a genomic approach was employed to evaluate the metabolic potentials and taxonomic classification of the halophilic genus Halarchaeum. Genomic analysis revealed that Halarchaeum members exhibit a predilection for amino acids as their primary energy source in high-salinity environments over carbohydrates. Genome analysis unveiled the presence of crucial genes associated with metabolic pathways, including the Embden-Meyerhof pathway, semi-phosphorylative Entner-Doudoroff pathway, and the urea cycle. Furthermore, the genomic analysis indicated that Halarchaeum members employ diverse mechanisms for osmotic regulation (encompassing both salt-in and salt-out strategies). Halarchaeum members also encode genes to alleviate acid and heat stress. The average nucleotide identity value between Halarchaeum solikamskense and Halarchaeum nitratireducens exceeded the established threshold (95%-96%) for defining distinct species. This high similarity suggests a close relationship between these two species, prompting the proposal to reclassify Halarchaeum solikamskense as a heterotypic synonym of Halarchaeum nitratireducens. The results of this study contribute to our knowledge of taxonomic classification and shed light on the adaptive strategies employed by Halarchaeum species in their specific ecological niches.

Unexpected carbon utilization activity of sulfate-reducing microorganisms in temperate and permanently cold marine sediments

Significant amounts of organic carbon in marine sediments are degraded, coupled with sulfate reduction. However, the actual carbon and energy sources used in situ have not been assigned to each group of diverse sulfate-reducing microorganisms (SRM) owing to the microbial and environmental complexity in sediments. Here, we probed microbial activity in temperate and permanently cold marine sediments by using potential SRM substrates, organic fermentation products at very low concentrations (15-30 μM), with RNA-based stable isotope probing. Unexpectedly, SRM were involved only to a minor degree in organic fermentation product mineralization, whereas metal-reducing microbes were dominant. Contrastingly, distinct SRM strongly assimilated 13C-DIC (dissolved inorganic carbon) with H2 as the electron donor. Our study suggests that canonical SRM prefer autotrophic lifestyle, with hydrogen as the electron donor, while metal-reducing microorganisms are involved in heterotrophic organic matter turnover, and thus regulate carbon fluxes in an unexpected way in marine sediments.

Isolation and biogeography of the oligotrophic ocean diazotroph, Crocosphaera waterburyi nov. sp

Marine N2-fixing cyanobacteria, including the unicellular genus Crocosphaera, are considered keystone species in marine food webs. Crocosphaera are globally distributed and provide new sources of nitrogen and carbon, which fuel oligotrophic microbial communities and upper trophic levels. Despite their ecosystem importance, only one pelagic, oligotrophic, phycoerythrin-rich species, Crocosphaera watsonii, has ever been identified and characterized as widespread. Herein, we present a new species, named Crocosphaera waterburyi, enriched from the North Pacific Ocean. C. waterburyi was found to be phenotypically and genotypically distinct from C. watsonii, active in situ, distributed globally, and preferred warmer temperatures in culture and the ocean. Additionally, C. waterburyi was detectable in 150- and 4000-meter sediment export traps, had a relatively larger biovolume than C. watsonii, and appeared to aggregate in the environment and laboratory culture. Therefore, it represents an additional, previously unknown link between atmospheric CO2 and N2 gas and deep ocean carbon and nitrogen export and sequestration.

Denitrification genotypes of endospore-forming Bacillota

Denitrification is a key metabolic process in the global nitrogen cycle and is performed by taxonomically diverse microorganisms. Despite the widespread importance of this metabolism, challenges remain in identifying denitrifying populations and predicting their metabolic end-products based on their genotype. Here, genome-resolved metagenomics was used to explore the denitrification genotype of Bacillota enriched in nitrate-amended high temperature incubations with confirmed N2O and N2 production. A set of 12 hidden Markov models (HMMs) was created to target the diversity of denitrification genes in members of the phylum Bacillota. Genomic potential for complete denitrification was found in five metagenome-assembled genomes from nitrate-amended enrichments, including two novel members of the Brevibacillaceae family. Genomes of complete denitrifiers encode N2O reductase gene clusters with clade II-type nosZ and often include multiple variants of the nitric oxide reductase gene. The HMM set applied to all genomes of Bacillota from the Genome Taxonomy Database identified 17 genera inferred to contain complete denitrifiers based on their gene content. Among complete denitrifiers it was common for three distinct nitric oxide reductases to be present (qNOR, bNOR, and sNOR) that may reflect the metabolic adaptability of Bacillota in environments with variable redox conditions.

Metagenome-assembled genomes of deep-sea sediments: changes in microbial functional potential lag behind redox transitions

Hadal sediments are hotspots of microbial activity in the deep sea and exhibit strong biogeochemical gradients. But although these gradients are widely assumed to exert selective forces on hadal microbial communities, the actual relationship between biogeochemistry, functional traits, and microbial community structure remains poorly understood. We tested whether the biogeochemical conditions in hadal sediments select for microbes based on their genomic capacity for respiration and carbohydrate utilization via a metagenomic analysis of over 153 samples from the Atacama Trench region (max. depth = 8085 m). The obtained 1357 non-redundant microbial genomes were affiliated with about one-third of all known microbial phyla, with more than half belonging to unknown genera. This indicated that the capability to withstand extreme hydrostatic pressure is a phylogenetically widespread trait and that hadal sediments are inhabited by diverse microbial lineages. Although community composition changed gradually over sediment depth, these changes were not driven by selection for respiratory or carbohydrate degradation capability in the oxic and nitrogenous zones, except in the case of anammox bacteria and nitrifying archaea. However, selection based on respiration and carbohydrate degradation capacity did structure the communities of the ferruginous zone, where aerobic and nitrogen respiring microbes declined exponentially (half-life = 125-419 years) and were replaced by subsurface communities. These results highlight a delayed response of microbial community composition to selective pressure imposed by redox zonation and indicated that gradual changes in microbial composition are shaped by the high-resilience and slow growth of microbes in the seafloor.

Phylogenetics and environmental distribution of nitric oxide-forming nitrite reductases reveal their distinct functional and ecological roles

The two evolutionarily unrelated nitric oxide-producing nitrite reductases, NirK and NirS, are best known for their redundant role in denitrification. They are also often found in organisms that do not perform denitrification. To assess the functional roles of the two enzymes and to address the sequence and structural variation within each, we reconstructed robust phylogenies of both proteins with sequences recovered from 6973 isolate and metagenome-assembled genomes and identified 32 well-supported clades of structurally distinct protein lineages. We then inferred the potential niche of each clade by considering other functional genes of the organisms carrying them as well as the relative abundances of each nir gene in 4082 environmental metagenomes across diverse aquatic, terrestrial, host-associated, and engineered biomes. We demonstrate that Nir phylogenies recapitulate ecology distinctly from the corresponding organismal phylogeny. While some clades of the nitrite reductase were equally prevalent across biomes, others had more restricted ranges. Nitrifiers make up a sizeable proportion of the nitrite-reducing community, especially for NirK in marine waters and dry soils. Furthermore, the two reductases showed distinct associations with genes involved in oxidizing and reducing other compounds, indicating that the NirS and NirK activities may be linked to different elemental cycles. Accordingly, the relative abundance and diversity of NirS versus NirK vary between biomes. Our results show the divergent ecological roles NirK and NirS-encoding organisms may play in the environment and provide a phylogenetic framework to distinguish the traits associated with organisms encoding the different lineages of nitrite reductases.

Succession of microbial community composition and secondary metabolism during marine biofilm development

In nature, secondary metabolites mediate interactions between microorganisms residing in complex microbial communities. However, the degree to which community dynamics can be linked to secondary metabolite potential remains largely unknown. In this study, we address the relationship between community succession and secondary metabolism variation. We used 16S and 18S rRNA gene and adenylation domain amplicon sequencing, genome-resolved metagenomics, and untargeted metabolomics to track the taxons, biosynthetic gene clusters, and metabolome dynamics in situ of microorganisms during marine biofilm succession over 113 days. Two phases were identified during the community succession, with a clear shift around Day 29, where the alkaloid secondary metabolites, pseudanes, were also detected. The microbial secondary metabolite potential changed between the phases, and only a few community members, including Myxococotta spp., were responsible for the majority of the biosynthetic gene cluster potential in the early succession phase. In the late phase, bryozoans and benthic copepods were detected, and the microbial nonribosomal peptide potential drastically decreased in association with a reduction in the relative abundance of the prolific secondary metabolite producers. Conclusively, this study provides evidence that the early succession of the marine biofilm community favors prokaryotes with high nonribosomal peptide synthetase potential. In contrast, the late succession is dominated by multicellular eukaryotes and a reduction in bacterial nonribosomal peptide synthetase potential.

High quality Bathyarchaeia MAGs from lignocellulose-impacted environments elucidate metabolism and evolutionary mechanisms

The archaeal class Bathyarchaeia is widely and abundantly distributed in anoxic habitats. Metagenomic studies have suggested that they are mixotrophic, capable of CO2 fixation and heterotrophic growth, and involved in acetogenesis and lignin degradation. We analyzed 35 Bathyarchaeia metagenome-assembled genomes (MAGs), including the first complete circularized MAG (cMAG) of the Bathy-6 subgroup, from the metagenomes of three full-scale pulp and paper mill anaerobic digesters and three laboratory methanogenic enrichment cultures maintained on pre-treated poplar. Thirty-three MAGs belong to the Bathy-6, lineage while two are from the Bathy-8 lineage. In our previous analysis of the microbial community in the pulp mill digesters, Bathyarchaeia were abundant and positively correlated to hydrogenotrophic and methylotrophic methanogenesis. Several factors likely contribute to the success of the Bathy-6 lineage compared to Bathy-8 in the reactors. The Bathy-6 genomes are larger than those of Bathy-8 and have more genes involved in lignocellulose degradation, including carbohydrate-active enzymes not present in the Bathy-8. Bathy-6 also shares the Bathyarchaeal O-demethylase system recently identified in Bathy-8. All the Bathy-6 MAGs had numerous membrane-associated pyrroloquinoline quinone-domain proteins that we suggest are involved in lignin modification or degradation, together with Radical-S-adenosylmethionine (SAM) and Rieske domain proteins, and AA2, AA3, and AA6-family oxidoreductases. We also identified a complete B12 synthesis pathway and a complete nitrogenase gene locus. Finally, comparative genomic analyses revealed that Bathyarchaeia genomes are dynamic and have interacted with other organisms in their environments through gene transfer to expand their gene repertoire.

Age, metabolisms, and potential origin of dominant anammox bacteria in the global oxygen-deficient zones

Anammox bacteria inhabiting oxygen-deficient zones (ODZs) are a major functional group mediating fixed nitrogen loss in the global ocean. However, many basic questions regarding the diversity, broad metabolisms, origin, and adaptive mechanisms of ODZ anammox bacteria remain unaddressed. Here we report two novel metagenome-assembled genomes of anammox bacteria affiliated with the Scalindua genus, which represent most, if not all, of the anammox bacteria in the global ODZs. Metagenomic read-recruiting and comparison with historical data show that they are ubiquitously present in all three major ODZs. Beyond the core anammox metabolism, both organisms contain cyanase, and the more dominant one encodes a urease, indicating most ODZ anammox bacteria can utilize cyanate and urea in addition to ammonium. Molecular clock analysis suggests that the evolutionary radiation of these bacteria into ODZs occurred no earlier than 310 million years ago, ~1 billion years after the emergence of the earliest modern-type ODZs. Different strains of the ODZ Scalindua species are also found in benthic sediments, and the first ODZ Scalindua is likely derived from the benthos. Compared to benthic strains of the same clade, ODZ Scalindua uniquely encodes genes for urea utilization but has lost genes related to growth arrest, flagellum synthesis, and chemotaxis, presumably for adaptation to thrive in the global ODZ waters. Our findings expand the known metabolisms and evolutionary history of the bacteria controlling the global nitrogen budget.

Meta-omics reveals role of photosynthesis in microbially induced carbonate precipitation at a CO2-rich geyser

Microbially induced carbonate precipitation (MICP) is a natural process with potential biotechnological applications to address both carbon sequestration and sustainable construction needs. However, our understanding of the microbial processes involved in MICP is limited to a few well-researched pathways such as ureolytic hydrolysis. To expand our knowledge of MICP, we conducted an omics-based study on sedimentary communities from travertine around the CO2-driven Crystal Geyser near Green River, Utah. Using metagenomics and metaproteomics, we identified the community members and potential metabolic pathways involved in MICP. We found variations in microbial community composition between the two sites we sampled, but Rhodobacterales were consistently the most abundant order, including both chemoheterotrophs and anoxygenic phototrophs. We also identified several highly abundant genera of Cyanobacteriales. The dominance of these community members across both sites and the abundant presence of photosynthesis-related proteins suggest that photosynthesis could play a role in MICP at Crystal Geyser. We also found abundant bacterial proteins involved in phosphorous starvation response at both sites suggesting that P-limitation shapes both composition and function of the microbial community driving MICP.

Complete genome sequence of a winter season Vibrio facilitates discovery of a novel subclade of cold-adapted species in the albus clade

In temperate marine climate zones, seasonal changes in water temperature contribute to distinct populations of warm- and cold-water vibrios. We report here the complete genome sequence (BUSCO score=94.8) of the novel strain Vibrio sp. VB16 isolated in late winter from the intertidal zone near Virginia Beach, Virginia, USA with the ability to form colonies at 4 °C. The 5.2 Mbp genome is composed of a large (3.6 Mbp) and small (1.6 Mbp) chromosome. Based on paired average nucleotide identity (ANI), average amino acid identity (AAI) and digital DNA-DNA hybridization (dDDH), V. sp. VB16 is the same species as V. sp. UBA2437 from a North Sea tidal flat and is closely related to V. sp. DW001 from Antarctic sea ice. Our phylogenomic and bioinformatic analyses placed VB16, UBA2437 and DW001 into a cold-tolerant subclade within the albus clade, along with two non-cold-tolerant subclades. Orthovenn analysis indicated that VB16 and its other albus clade members shared 1544 gene orthologue clusters, including clusters for biosynthesis of polar flagella and tight adhesion pili that predict multiple lifestyles, either free-living or as an opportunistic pathogen within a marine eukaryotic host. The cold-tolerant subclade shared 552 orthologue proteins, including genes known to promote survival in cold or freezing temperatures, such as the eicosapentaenoic acid biosynthetic gene cluster, syp exopolysaccharide gene cluster and novel giant proteins with ice-binding domains. This subclade represents a group of psychrotolerant or 'moderate psychrophile' winter season Vibrio species. The discovery of this subclade opens the door for experimental work on the physiological features, virulence potential and ecological importance of this subclade.

A novel UV-resistant bacterium Sphingomonas endolithica sp. nov., and genomic analysis, isolated from the north slope of Mount Everest

Gram-stain-negative, aerobic, rod-shaped, non-motile bacterium strain ZFBP2030T was isolated from a rock on the North slope of Mount Everest. This strain contained a unique ubiquinone-10 (Q-10) as a predominant respiratory quinone. Among the tested fatty acids, the strain contained summed feature 8, C14:0 2OH, and C16:0, as major cellular fatty acids. The polar lipid profile contained phosphatidyl glycerol, phosphatidyl ethanolamine, three unidentified phospholipids, two unidentified aminolipids, and six unidentified lipids. The cell-wall peptidoglycan was a meso-diaminopimelic acid, and cell-wall sugars were ribose and galactose. Phylogenetic analyses based on 16S rRNA gene sequence revealed that strain ZFBP2030T was a member of the genus Sphingomonas, exhibiting high sequence similarity to the 16S rRNA gene sequences of Sphingomonas aliaeris DH-S5T (97.9%), Sphingomonas alpina DSM 22537T (97.3%) and Sphingomonas hylomeconis CCTCC AB 2013304T (97.0%). The 16S rRNA gene sequence similarity between ZFBP2030T and other typical strains was less than 97.0%. The average amino acid identity values, average nucleotide identity, and digital DNA-DNA hybridization values between strain ZFBP2030T and its highest sequence similarity strains were 56.9-79.9%, 65.1-82.2%, and 19.3-25.8%, respectively. The whole-genome size of the novel strain ZFBP2030T was 4.1 Mbp, annotated with 3838 protein-coding genes and 54 RNA genes. Moreover, DNA G + C content was 64.7 mol%. Stress-related functions predicted in the subsystem classification of the strain ZFBP2030T genome included osmotic, oxidative, cold/heat shock, detoxification, and periplasmic stress responses. The overall results of this study clearly showed that strain ZFBP2030T is a novel species of the genus Sphingomonas, for which the name Sphingomonas endolithica sp. nov. is proposed. The type of strain is ZFBP2030T (= EE 013T = GDMCC 1.3123T = JCM 35386T).

Distribution and genomic variation of ammonia-oxidizing archaea in abyssal and hadal surface sediments

Ammonia-oxidizing archaea of the phylum Thaumarchaeota play a central role in the biogeochemical cycling of nitrogen in benthic sediments, at the interface between pelagic and subsurface ecosystems. However, our understanding of their niche separation and of the processes controlling their population structure in hadal and abyssal surface sediments is still limited. Here, we reconstructed 47 AOA metagenome-assembled genomes (MAGs) from surface sediments of the Atacama and Kermadec trench systems. They formed deep-sea-specific groups within the family Nitrosopumilaceae and were assigned to six amoA gene-based clades. MAGs from different clades had distinct distribution patterns along oxygen-ammonium counter gradients in surface sediments. At the species level, MAGs thus seemed to form different ecotypes and follow deterministic niche-based distributions. In contrast, intraspecific population structure, defined by patterns of Single Nucleotide Variants (SNV), seemed to reflect more complex contributions of both deterministic and stochastic processes. Firstly, the bathymetric range had a strong effect on population structure, with distinct populations in abyssal plains and hadal trenches. Then, hadal populations were clearly separated by trench system, suggesting a strong isolation-by-topography effect, whereas abyssal populations were rather controlled by sediment depth or geographic distances, depending on the clade considered. Interestingly, genetic variability between samples was lowest in sediment layers where the mean MAG coverage was highest, highlighting the importance of selective pressure linked with each AOA clade's ecological niche. Overall, our results show that deep-sea AOA genome distributions seem to follow both deterministic and stochastic processes, depending on the genomic variability scale considered.

Taxonomic distribution of metabolic functions in bacteria associated with Trichodesmium consortia

IMPORTANCE

Colonies of the cyanobacteria Trichodesmium act as a biological hotspot for the usage and recycling of key resources such as C, N, P, and Fe within an otherwise oligotrophic environment. While Trichodesmium colonies are known to interact and support a unique community of algae and particle-associated microbes, our understanding of the taxa that populate these colonies and the gene functions they encode is still limited. Characterizing the taxa and adaptive strategies that influence consortium physiology and its concomitant biogeochemistry is critical in a future ocean predicted to have increasingly resource-depleted regions.

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton

IMPORTANCE

Single-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

Persistent enrichment of multidrug resistant Klebsiella in oral and nasal communities during long-term starvation

The human oral and nasal cavities can act as reservoirs for opportunistic pathogens capable of causing acute infection. These microbes asymptomatically colonize the human oral and nasal cavities which facilitates transmission within human populations via the environment, and they routinely possess a clinically-significant antibiotic-resistance genes. Among these opportunistic pathogens, the Klebsiella genus stands out as a notable example, with its members frequently linked to nosocomial infections and multidrug resistance. As with many colonizing opportunistic pathogens, how Klebsiella transitions from an asymptomatic colonizer to a pathogen remains unclear. Here, we explored a possible explanation by investigating the ability of oral and nasal Klebsiella to outcompete their native microbial community members under in vitro starvation conditions, which could be analogous to external hospital environments. When Klebsiella was present within a healthy human oral or nasal sample, the bacterial community composition shifted dramatically under starvation conditions and typically became dominated by Klebsiella. Furthermore, introducing K. pneumoniae exogenously into a native microbial community lacking K. pneumoniae, even at low inoculum, led to repeated dominance under starvation. K.pneumoniae strains isolated from healthy individuals' oral and nasal cavities also exhibited resistance to multiple classes of antibiotics and were genetically similar to clinical and gut isolates. In addition, we found that in the absence of Klebsiella, other understudied opportunistic pathogens, such as Peptostreptococcus, dominate under starvation conditions. Our findings establish an environmental circumstance that allows for the outgrowth of Klebsiella and other opportunistic pathogens. The ability to outcompete other commensal bacteria and to persist under harsh environmental conditions may contribute to the colonization-to-infection transition of these opportunistic pathogens.

Endophytic bacteria associated with wild-type banana seed (Musa balbisiana): whole genome sequencing

We present the whole-genome sequences of five endophytic bacteria isolated from Musa balbisiana seeds. These strains represent five different genera: Bacillus, Brachybacterium, Enterobacter, Enterococcus, and Pantoea. Among these, three genera (Bacillus, Pantoea, and Enterobacter) were previously recognized for their antagonistic effects against Fusarium wilt, a highly destructive disease that affects banana plants.

TADA: taxonomy-aware dataset aggregator

SUMMARY

The profusion of sequenced genomes across the bacterial and archeal domains offers unprecedented possibilities for phylogenetic and comparative genomic analyses. In general, phylogenetic reconstruction is improved by the use of more data. However, including all available data is (i) not computationally tractable, and (ii) prone to biases, as the abundance of genomes is very unequally distributed over the biological diversity. Thus, in most cases, subsampling taxa to build a phylogeny is necessary. Currently, though, there is no available software to perform that handily. Here we present TADA, a taxonomic-aware dataset selection workflow that allows sampling across user-defined portions of the prokaryotic diversity with variable granularity, while setting constraints on genome quality and balance between branches.

AVAILABILITY AND IMPLEMENTATION

TADA is implemented as a snakemake workflow and is freely available at https://github.com/emilhaegglund/TADA.

A metagenomic view of novel microbial and metabolic diversity found within the deep terrestrial biosphere at DeMMO: A microbial observatory in South Dakota, USA

The deep terrestrial subsurface is a large and diverse microbial habitat and vast repository of biomass. However, in relation to its size and physical heterogeneity we have limited understanding of taxonomic and metabolic diversity in this realm. Here we present a detailed metagenomic analysis of samples from the Deep Mine Microbial Observatory (DeMMO) spanning depths from the surface to 1.5 km into the crust. From eight geochemically and spatially distinct fluid samples we reconstructed ~600 partial to near-complete metagenome-assembled genomes (MAGs), representing 50 distinct phyla and including 18 candidate phyla. These novel clades include members of the candidate phyla radiation, two new MAGs from OLB16, a phylum originally identified in DeMMO fluids and for which only one other MAG is currently available, and new MAGs from the Eisenbacteria, Omnitrophota, and Edwardsbacteria. We find that microbes spanning this expansive phylogenetic diversity and physical subsurface space gain a competitive edge by maintaining a wide variety of functional pathways, are often capable of numerous dissimilatory energy metabolisms and poised to take advantage of nutrients as they become available in isolated fracture fluids. Our results support and expand on emerging themes of tight nutrient cycling and genomic plasticity in deep subsurface biosphere taxa.

Hadza Prevotella require diet-derived microbiota-accessible carbohydrates to persist in mice

Industrialization has transformed the gut microbiota, reducing the prevalence of Prevotella relative to Bacteroides. Here, we isolate Bacteroides and Prevotella strains from the microbiota of Hadza hunter-gatherers in Tanzania, a population with high levels of Prevotella. We demonstrate that plant-derived microbiota-accessible carbohydrates (MACs) are required for persistence of Prevotella copri but not Bacteroides thetaiotaomicron in vivo. Differences in carbohydrate metabolism gene content, expression, and in vitro growth reveal that Hadza Prevotella strains specialize in degrading plant carbohydrates, while Hadza Bacteroides isolates use both plant and host-derived carbohydrates, a difference mirrored in Bacteroides from non-Hadza populations. When competing directly, P. copri requires plant-derived MACs to maintain colonization in the presence of B. thetaiotaomicron, as a no-MAC diet eliminates P. copri colonization. Prevotella's reliance on plant-derived MACs and Bacteroides' ability to use host mucus carbohydrates could explain the reduced prevalence of Prevotella in populations consuming a low-MAC, industrialized diet.

Microbial methane cycling in a landfill on a decadal time scale

Landfills generate outsized environmental footprints due to microbial degradation of organic matter in municipal solid waste, which produces the potent greenhouse gas methane. With global solid waste production predicted to increase substantially in the next few decades, there is a pressing need to better understand the temporal dynamics of biogeochemical processes that control methane cycling in landfills. Here, we use metagenomic approaches to characterize microbial methane cycling in waste that was landfilled over 39 years. Our analyses indicate that newer waste supports more diverse communities with similar composition compared to older waste, which contains lower diversity and more varied communities. Older waste contains primarily autotrophic organisms with versatile redox metabolisms, whereas newer waste is dominated by anaerobic fermenters. Methane-producing microbes are more abundant, diverse, and metabolically versatile in new waste compared to old waste. Our findings indicate that predictive models for methane emission in landfills overlook methane oxidation in the absence of oxygen, as well as certain microbial lineages that can potentially contribute to methane sinks in diverse habitats.

A genome catalog of the early-life human skin microbiome

BACKGROUND

Metagenome-assembled genomes have greatly expanded the reference genomes for skin microbiome. However, the current reference genomes are largely based on samples from adults in North America and lack representation from infants and individuals from other continents.

RESULTS

Here we use deep shotgun metagenomic sequencing to profile the skin microbiota of 215 infants at age 2-3 months and 12 months who are part of the VITALITY trial in Australia as well as 67 maternally matched samples. Based on the infant samples, we present the Early-Life Skin Genomes (ELSG) catalog, comprising 9483 prokaryotic genomes from 1056 species, 206 fungal genomes from 13 species, and 39 eukaryotic viral sequences. This genome catalog substantially expands the diversity of species previously known to comprise human skin microbiome and improves the classification rate of sequenced data by 21%. The protein catalog derived from these genomes provides insights into the functional elements such as defense mechanisms that distinguish early-life skin microbiome. We also find evidence for microbial sharing at the community, bacterial species, and strain levels between mothers and infants.

CONCLUSIONS

Overall, the ELSG catalog uncovers the skin microbiome of a previously underrepresented age group and population and provides a comprehensive view of human skin microbiome diversity, function, and development in early life.

Phylogenomics, phenotypic, and functional traits of five novel (Earth-derived) bacterial species isolated from the International Space Station and their prevalence in metagenomes

With the advent of long-term human habitation in space and on the moon, understanding how the built environment microbiome of space habitats differs from Earth habitats, and how microbes survive, proliferate and spread in space conditions, is becoming more important. The microbial tracking mission series has been monitoring the microbiome of the International Space Station (ISS) for almost a decade. During this mission series, six unique strains of Gram-stain-positive bacteria, including two spore-forming and three non-spore-forming species, were isolated from the environmental surfaces of the ISS. The analysis of their 16S rRNA gene sequences revealed > 99% similarities with previously described bacterial species. To further explore their phylogenetic affiliation, whole genome sequencing was undertaken. For all strains, the gyrB gene exhibited < 93% similarity with closely related species, which proved effective in categorizing these ISS strains as novel species. Average nucleotide identity and digital DNA-DNA hybridization values, when compared to any known bacterial species, were < 94% and <50% respectively for all species described here. Traditional biochemical tests, fatty acid profiling, polar lipid, and cell wall composition analyses were performed to generate phenotypic characterization of these ISS strains. A study of the shotgun metagenomic reads from the ISS samples, from which the novel species were isolated, showed that only 0.1% of the total reads mapped to the novel species, supporting the idea that these novel species are rare in the ISS environments. In-depth annotation of the genomes unveiled a variety of genes linked to amino acid and derivative synthesis, carbohydrate metabolism, cofactors, vitamins, prosthetic groups, pigments, and protein metabolism. Further analysis of these ISS-isolated organisms revealed that, on average, they contain 46 genes associated with virulence, disease, and defense. The main predicted functions of these genes are: conferring resistance to antibiotics and toxic compounds, and enabling invasion and intracellular resistance. After conducting antiSMASH analysis, it was found that there are roughly 16 cluster types across the six strains, including β-lactone and type III polyketide synthase (T3PKS) clusters. Based on these multi-faceted taxonomic methods, it was concluded that these six ISS strains represent five novel species, which we propose to name as follows: Arthrobacter burdickii IIF3SC-B10T (= NRRL B-65660T = DSM 115933T), Leifsonia virtsii F6_8S_P_1AT (= NRRL B-65661T = DSM 115931T), Leifsonia williamsii F6_8S_P_1BT (= NRRL B-65662T = DSM 115932T), Paenibacillus vandeheii F6_3S_P_1CT (= NRRL B-65663T = DSM 115940T), and Sporosarcina highlanderae F6_3S_P_2T (= NRRL B-65664T = DSM 115943T). Identifying and characterizing the genomes and phenotypes of novel microbes found in space habitats, like those explored in this study, is integral for expanding our genomic databases of space-relevant microbes. This approach offers the only reliable method to determine species composition, track microbial dispersion, and anticipate potential threats to human health from monitoring microbes on the surfaces and equipment within space habitats. By unraveling these microbial mysteries, we take a crucial step towards ensuring the safety and success of future space missions.

Desulfoferula mesophilus gen. nov. sp. nov., a mesophilic sulfate-reducing bacterium isolated from a brackish lake sediment

A novel sulfate-reducing bacterium (strain 12FAKT) was isolated from sediment sampled from a brackish lake in Japan. Respiratory growth was observed with formate and pyruvate as an electron donor. Sulfate, thiosulfate, elemental sulfur and dimethyl sulfoxide were utilized as an electron acceptor. The isolate grew over a temperature range of 18-42 °C (optimum 35-37 °C), a NaCl concentration range of 50-450 mM (optimum 150-300 mM) and a pH range of 6.6-7.5. The 12FAKT genome consists of a circular chromosome with a length of 4.5 Mbp and G + C content of 63.6%. Based on 16S rRNA gene sequence similarity, the closest cultured relative was Desulfarculus baarsii 2st14T (92.2%). Genome-based phylogenetic analysis placed strain 12FAKT within the family Desulfarculaceae but did not affiliate the strain with any existing genus. Taken together, we propose a novel species of a novel genus, Desulfoferula mesophilus gen. nov. sp. nov. with the type strain 12FAKT (= DSM 115219T = JCM 39399T).

Genes and genome-resolved metagenomics reveal the microbial functional make up of Amazon peatlands under geochemical gradients

The Pastaza-Marañón Foreland Basin (PMFB) holds the most extensive tropical peatland area in South America. PMFB peatlands store ~7.07 Gt of organic carbon interacting with multiple microbial heterotrophic, methanogenic, and other aerobic/anaerobic respirations. Little is understood about the contribution of distinct microbial community members inhabiting tropical peatlands. Here, we studied the metagenomes of three geochemically distinct peatlands spanning minerotrophic, mixed, and ombrotrophic conditions. Using gene- and genome-centric approaches, we evaluate the functional potential of the underlying microbial communities. Abundance analyses show significant differences in C, N, P, and S acquisition genes. Furthermore, community interactions mediated by toxin-antitoxin and CRISPR-Cas systems were enriched in oligotrophic soils, suggesting that non-metabolic interactions may exert additional controls in low-nutrient environments. Additionally, we reconstructed 519 metagenome-assembled genomes spanning 28 phyla. Our analyses detail key differences across the geochemical gradient in the predicted microbial populations involved in degradation of organic matter, and the cycling of N and S. Notably, we observed differences in the nitric oxide (NO) reduction strategies between sites with high and low N2 O fluxes and found phyla putatively capable of both NO and sulfate reduction. Our findings detail how gene abundances and microbial populations are influenced by geochemical differences in tropical peatlands.

Microbial-enrichment method enables high-throughput metagenomic characterization from host-rich samples

Host-microbe interactions have been linked to health and disease states through the use of microbial taxonomic profiling, mostly via 16S ribosomal RNA gene sequencing. However, many mechanistic insights remain elusive, in part because studying the genomes of microbes associated with mammalian tissue is difficult due to the high ratio of host to microbial DNA in such samples. Here we describe a microbial-enrichment method (MEM), which we demonstrate on a wide range of sample types, including saliva, stool, intestinal scrapings, and intestinal mucosal biopsies. MEM enabled high-throughput characterization of microbial metagenomes from human intestinal biopsies by reducing host DNA more than 1,000-fold with minimal microbial community changes (roughly 90% of taxa had no significant differences between MEM-treated and untreated control groups). Shotgun sequencing of MEM-treated human intestinal biopsies enabled characterization of both high- and low-abundance microbial taxa, pathways and genes longitudinally along the gastrointestinal tract. We report the construction of metagenome-assembled genomes directly from human intestinal biopsies for bacteria and archaea at relative abundances as low as 1%. Analysis of metagenome-assembled genomes reveals distinct subpopulation structures between the small and large intestine for some taxa. MEM opens a path for the microbiome field to acquire deeper insights into host-microbe interactions by enabling in-depth characterization of host-tissue-associated microbial communities.

Metagenome-assembled genomes from High Arctic glaciers highlight the vulnerability of glacier-associated microbiota and their activities to habitat loss

The rapid warming of the Arctic is threatening the demise of its glaciers and their associated ecosystems. Therefore, there is an urgent need to explore and understand the diversity of genomes resident within glacial ecosystems endangered by human-induced climate change. In this study we use genome-resolved metagenomics to explore the taxonomic and functional diversity of different habitats within glacier-occupied catchments. Comparing different habitats within such catchments offers a natural experiment for understanding the effects of changing habitat extent or even loss upon Arctic microbiota. Through binning and annotation of metagenome-assembled genomes (MAGs) we describe the spatial differences in taxon distribution and their implications for glacier-associated biogeochemical cycling. Multiple taxa associated with carbon cycling included organisms with the potential for carbon monoxide oxidation. Meanwhile, nitrogen fixation was mediated by a single taxon, although diverse taxa contribute to other nitrogen conversions. Genes for sulphur oxidation were prevalent within MAGs implying the potential capacity for sulphur cycling. Finally, we focused on cyanobacterial MAGs, and those within cryoconite, a biodiverse microbe-mineral granular aggregate responsible for darkening glacier surfaces. Although the metagenome-assembled genome of Phormidesmis priestleyi, the cyanobacterium responsible for forming Arctic cryoconite was represented with high coverage, evidence for the biosynthesis of multiple vitamins and co-factors was absent from its MAG. Our results indicate the potential for cross-feeding to sustain P. priestleyi within granular cryoconite. Taken together, genome-resolved metagenomics reveals the vulnerability of glacier-associated microbiota to the deletion of glacial habitats through the rapid warming of the Arctic.

EasyCGTree: a pipeline for prokaryotic phylogenomic analysis based on core gene sets

BACKGROUND

Genome-scale phylogenetic analysis based on core gene sets is routinely used in microbiological research. However, the techniques are still not approachable for individuals with little bioinformatics experience. Here, we present EasyCGTree, a user-friendly and cross-platform pipeline to reconstruct genome-scale maximum-likehood (ML) phylogenetic tree using supermatrix (SM) and supertree (ST) approaches.

RESULTS

EasyCGTree was implemented in Perl programming languages and was built using a collection of published reputable programs. All the programs were precompiled as standalone executable files and contained in the EasyCGTree package. It can run after installing Perl language environment. Several profile hidden Markov models (HMMs) of core gene sets were prepared in advance to construct a profile HMM database (PHD) that was enclosed in the package and available for homolog searching. Customized gene sets can also be used to build profile HMM and added to the PHD via EasyCGTree. Taking 43 genomes of the genus Paracoccus as the testing data set, consensus (a variant of the typical SM), SM, and ST trees were inferred via EasyCGTree successfully, and the SM trees were compared with those inferred via the pipelines UBCG and bcgTree, using the metrics of cophenetic correlation coefficients (CCC) and Robinson-Foulds distance (topological distance). The results suggested that EasyCGTree can infer SM trees with nearly identical topology (distance < 0.1) and accuracy (CCC > 0.99) to those of trees inferred with the two pipelines.

CONCLUSIONS

EasyCGTree is an all-in-one automatic pipeline from input data to phylogenomic tree with guaranteed accuracy, and is much easier to install and use than the reference pipelines. In addition, ST is implemented in EasyCGTree conveniently and can be used to explore prokaryotic evolutionary signals from a different perspective. The EasyCGTree version 4 is freely available for Linux and Windows users at Github ( https://github.com/zdf1987/EasyCGTree4 ).

Phycobacter azelaicus gen. nov. sp. nov., a diatom symbiont isolated from the phycosphere of Asterionellopsis glacialis

A diatom-associated bacterium, designated as strain F10T, was isolated from a pure culture of the pennate diatom Asterionellopsis glacialis A3 and has since been used to characterize molecular mechanisms of symbiosis between phytoplankton and bacteria, including interactions using diatom-derived azelaic acid. Its origin from a hypersaline environment, combined with its capacity for quorum sensing, biofilm formation, and potential for dimethylsulfoniopropionate methylation/cleavage, suggest it is within the family Roseobacteraceae. Initial phylogenetic analysis of the 16S rRNA gene sequence placed this isolate within the Phaeobacter genus, but recent genomic and phylogenomic analyses show strain F10T is a separate lineage diverging from the genus Pseudophaeobacter. The genomic DNA G+C content is 60.0 mol%. The predominant respiratory quinone is Q-10. The major fatty acids are C18 : 1 ω7c and C16 : 0. Strain F10T also contains C10 : 03-OH and the furan-containing fatty acid 10,13-epoxy-11-methyl-octadecadienoate (9-(3-methyl-5-pentylfuran-2-yl)nonanoic acid). The major polar lipids are diphosphatidylglycerol, phosphatidylethanolamine and phosphatidylglycerol. Based on genomic, phylogenomic, phenotypic and chemotaxonomic characterizations, strain F10T represents a novel genus and species with the proposed name, Phycobacter azelaicus gen. nov. sp. nov. The type strain is F10T (=NCMA B37T=NCIMB 15470T=NRIC 2002T).

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.

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).

Conserved functions of prion candidates suggest a primeval role of protein self-templating

Amyloid-based prions have simple structures, a wide phylogenetic distribution, and a plethora of functions in contemporary organisms, suggesting they may be an ancient phenomenon. However, this hypothesis has yet to be addressed with a systematic, computational, and experimental approach. Here we present a framework to help guide future experimental verification of candidate prions with conserved functions to understand their role in the early stages of evolution and potentially in the origins of life. We identified candidate prions in all high-quality proteomes available in UniProt computationally, assessed their phylogenomic distributions, and analyzed candidate-prion functional annotations. Of the 27 980 560 proteins scanned, 228 561 were identified as candidate prions (~0.82%). Among these candidates, there were 84 Gene Ontology (GO) terms conserved across the three domains of life. We found that candidate prions with a possible role in adaptation were particularly well-represented within this group. We discuss unifying features of candidate prions to elucidate the primeval roles of prions and their associated functions. Candidate prions annotated as transcription factors, DNA binding, and kinases are particularly well suited to generating diverse responses to changes in their environment and could allow for adaptation and population expansion into more diverse environments. We hypothesized that a relationship between these functions and candidate prions could be evolutionarily ancient, even if individual prion domains themselves are not evolutionarily conserved. Candidate prions annotated with these universally occurring functions potentially represent the oldest extant prions on Earth and are therefore excellent experimental targets.

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.

Dynamic genetic adaptation of Bacteroides thetaiotaomicron during murine gut colonization

To understand how a bacterium ultimately succeeds or fails in adapting to a new host, it is essential to assess the temporal dynamics of its fitness over the course of colonization. Here, we introduce a human-derived commensal organism, Bacteroides thetaiotaomicron (Bt), into the guts of germ-free mice to determine whether and how the genetic requirements for colonization shift over time. Combining a high-throughput functional genetics assay and transcriptomics, we find that gene usage changes drastically during the first days of colonization, shifting from high expression of amino acid biosynthesis genes to broad upregulation of diverse polysaccharide utilization loci. Within the first week, metabolism becomes centered around utilization of a predominant dietary oligosaccharide, and these changes are largely sustained through 6 weeks of colonization. Spontaneous mutations in wild-type Bt also evolve around this locus. These findings highlight the importance of considering temporal colonization dynamics in developing more effective microbiome-based therapies.

Revealing the Host-Dependent Nature of an Engineered Genetic Inverter in Concordance with Physiology

Broad-host-range synthetic biology is an emerging frontier that aims to expand our current engineerable domain of microbial hosts for biodesign applications. As more novel species are brought to "model status," synthetic biologists are discovering that identically engineered genetic circuits can exhibit different performances depending on the organism it operates within, an observation referred to as the "chassis effect." It remains a major challenge to uncover which genome-encoded and biological determinants will underpin chassis effects that govern the performance of engineered genetic devices. In this study, we compared model and novel bacterial hosts to ask whether phylogenomic relatedness or similarity in host physiology is a better predictor of genetic circuit performance. This was accomplished using a comparative framework based on multivariate statistical approaches to systematically demonstrate the chassis effect and characterize the performance dynamics of a genetic inverter circuit operating within 6 Gammaproteobacteria. Our results solidify the notion that genetic devices are strongly impacted by the host context. Furthermore, we formally determined that hosts exhibiting more similar metrics of growth and molecular physiology also exhibit more similar performance of the genetic inverter, indicating that specific bacterial physiology underpins measurable chassis effects. The result of this study contributes to the field of broad-host-range synthetic biology by lending increased predictive power to the implementation of genetic devices in less-established microbial hosts.

Vitamin interdependencies predicted by metagenomics-informed network analyses and validated in microbial community microcosms

Metagenomic or metabarcoding data are often used to predict microbial interactions in complex communities, but these predictions are rarely explored experimentally. Here, we use an organism abundance correlation network to investigate factors that control community organization in mine tailings-derived laboratory microbial consortia grown under dozens of conditions. The network is overlaid with metagenomic information about functional capacities to generate testable hypotheses. We develop a metric to predict the importance of each node within its local network environments relative to correlated vitamin auxotrophs, and predict that a Variovorax species is a hub as an important source of thiamine. Quantification of thiamine during the growth of Variovorax in minimal media show high levels of thiamine production, up to 100 mg/L. A few of the correlated thiamine auxotrophs are predicted to produce pantothenate, which we show is required for growth of Variovorax, supporting that a subset of vitamin-dependent interactions are mutualistic. A Cryptococcus yeast produces the B-vitamin pantothenate, and co-culturing with Variovorax leads to a 90-130-fold fitness increase for both organisms. Our study demonstrates the predictive power of metagenome-informed, microbial consortia-based network analyses for identifying microbial interactions that underpin the structure and functioning of microbial communities.

Thiovibrio frasassiensis gen. nov., sp. nov., an autotrophic, elemental sulphur disproportionating bacterium isolated from sulphidic karst sediment, and proposal of Thiovibrionaceae fam. nov

A novel, autotrophic, mesophilic bacterium, strain RS19-109T, was isolated from sulphidic stream sediments in the Frasassi Caves, Italy. The cells of this strain grew chemolithoautotrophically under anaerobic conditions while disproportionating elemental sulphur (S0) and thiosulphate, but not sulphite with bicarbonate/CO2 as a carbon source. Autotrophic growth was also observed with molecular hydrogen as an electron donor, and S0, sulphate, thiosulphate, nitrate and ferric iron as electron acceptors. Oxygen was not used as an electron acceptor and sulphide was not used as an electron donor. Weak growth was observed with sulphate as an electron acceptor and organic carbon as an electron donor and carbon source. The strain also showed weak growth by fermentation of tryptone. It grew at pH 5.5–7.5 (optimum, pH 7.0), 4–35 °C (optimum, 30 °C) and between 0–1.7 % NaCl. Strain RS19-109T was found to be phylogenetically distinct based on 16S rRNA gene sequence similarity (89.2 %) to its closest relative, Desulfurivibrio alkaliphilus AHT2T. The draft genome sequence for strain RS19-109T had average nucleotide identity, average amino acid identity and in silico DNA–DNA hybridization values of 72.2, 63.0 and 18.3 %, respectively, compared with the genome sequence of D. alkaliphilus AHT2T. On the basis of its physiological and genomic properties, strain RS19-109T is proposed as the type strain of a novel species of a novel genus, Thiovibrio frasassiensis gen. nov., sp. nov. A novel family, Thiovibrionaceae fam. nov., is proposed to accommodate Thiovibrio within the order Desulfobulbales . Strain RS19-109T has been deposited at the DSMZ-German Collection of Microorganisms and Cell Cultures (=DSM 115074T) and the American Type Culture Collection (=ATCC TSD-325T).

Mining-impacted rice paddies select for Archaeal methylators and reveal a putative (Archaeal) regulator of mercury methylation

Methylmercury (MeHg) is a microbially produced neurotoxin derived from inorganic mercury (Hg), which accumulation in rice represents a major health concern to humans. However, the microbial control of MeHg dynamics in the environment remains elusive. Here, leveraging three rice paddy fields with distinct concentrations of Hg (Total Hg (THg): 0.21-513 mg kg-1 dry wt. soil; MeHg: 1.21-6.82 ng g-1 dry wt. soil), we resorted to metagenomics to determine the microbial determinants involved in MeHg production under contrasted contamination settings. We show that Hg methylating Archaea, along with methane-cycling genes, were enriched in severely contaminated paddy soils. Metagenome-resolved Genomes of novel putative Hg methylators belonging to Nitrospinota (UBA7883), with poorly resolved taxonomy despite high completeness, showed evidence of facultative anaerobic metabolism and adaptations to fluctuating redox potential. Furthermore, we found evidence of environmental filtering effects that influenced the phylogenies of not only hgcA genes under different THg concentrations, but also of two housekeeping genes, rpoB and glnA, highlighting the need for further experimental validation of whether THg drives the evolution of hgcAB. Finally, assessment of the genomic environment surrounding hgcAB suggests that this gene pair may be regulated by an archaeal toxin-antitoxin (TA) system, instead of the more frequently found arsR-like genes in bacterial methylators. This suggests the presence of distinct hgcAB regulation systems in bacteria and archaea. Our results support the emerging role of Archaea in MeHg cycling under mining-impacted environments and shed light on the differential control of the expression of genes involved in MeHg formation between Archaea and Bacteria.

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.

The evolution and spread of sulfur cycling enzymes reflect the redox state of the early Earth

The biogeochemical sulfur cycle plays a central role in fueling microbial metabolisms, regulating the Earth's redox state, and affecting climate. However, geochemical reconstructions of the ancient sulfur cycle are confounded by ambiguous isotopic signals. We use phylogenetic reconciliation to ascertain the timing of ancient sulfur cycling gene events across the tree of life. Our results suggest that metabolisms using sulfide oxidation emerged in the Archean, but those involving thiosulfate emerged only after the Great Oxidation Event. Our data reveal that observed geochemical signatures resulted not from the expansion of a single type of organism but were instead associated with genomic innovation across the biosphere. Moreover, our results provide the first indication of organic sulfur cycling from the Mid-Proterozoic onwards, with implications for climate regulation and atmospheric biosignatures. Overall, our results provide insights into how the biological sulfur cycle evolved in tandem with the redox state of the early Earth.

Phylogenetic affiliations and genomic characterization of novel bacterial species and their abundance in the International Space Station

Background

With the advent of long-term human habitation in space and on the moon, understanding how the built environment microbiome of space habitats differs from Earth habits, and how microbes survive, proliferate and spread in space conditions, is coming more and more important. The Microbial Tracking mission series has been monitoring the microbiome of the International Space Station (ISS) for almost a decade. During this mission series, six unique strains of Gram-positive bacteria, including two spore-forming and three non-spore-forming species, were isolated from the environmental surfaces of the International Space Station (ISS).

Results

The analysis of their 16S rRNA gene sequences revealed <99% similarities with previously described bacterial species. To further explore their phylogenetic affiliation, whole genome sequencing (WGS) was undertaken. For all strains, the gyrB gene exhibited <93% similarity with closely related species, which proved effective in categorizing these ISS strains as novel species. Average ucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values, when compared to any known bacterial species, were less than <94% and 50% respectively for all species described here. Traditional biochemical tests, fatty acid profiling, polar lipid, and cell wall composition analyses were performed to generate phenotypic characterization of these ISS strains. A study of the shotgun metagenomic reads from the ISS samples, from which the novel species were isolated, showed that only 0.1% of the total reads mapped to the novel species, supporting the idea that these novel species are rare in the ISS environments. In-depth annotation of the genomes unveiled a variety of genes linked to amino acid and derivative synthesis, carbohydrate metabolism, cofactors, vitamins, prosthetic groups, pigments, and protein metabolism. Further analysis of these ISS-isolated organisms revealed that, on average, they contain 46 genes associated with virulence, disease, and defense. The main predicted functions of these genes are: conferring resistance to antibiotics and toxic compounds, and enabling invasion and intracellular resistance. After conducting antiSMASH analysis, it was found that there are roughly 16 cluster types across the six strains, including β-lactone and type III polyketide synthase (T3PKS) clusters.

Conclusions

Based on these multi-faceted taxonomic methods, it was concluded that these six ISS strains represent five novel species, which we propose to name as follows: Arthrobacter burdickii IIF3SC-B10T (=NRRL B-65660T), Leifsonia virtsii, F6_8S_P_1AT (=NRRL B-65661T), Leifsonia williamsii, F6_8S_P_1BT (=NRRL B- 65662T and DSMZ 115932T), Paenibacillus vandeheii, F6_3S_P_1CT(=NRRL B-65663T and DSMZ 115940T), and Sporosarcina highlanderae F6_3S_P_2 T(=NRRL B-65664T and DSMZ 115943T). Identifying and characterizing the genomes and phenotypes of novel microbes found in space habitats, like those explored in this study, is integral for expanding our genomic databases of space-relevant microbes. This approach offers the only reliable method to determine species composition, track microbial dispersion, and anticipate potential threats to human health from monitoring microbes on the surfaces and equipment within space habitats. By unraveling these microbial mysteries, we take a crucial step towards ensuring the safety and success of future space missions.

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.