GToTree: a user-friendly workflow for phylogenomics

Sulfoquinovose is a select nutrient of prominent bacteria and a source of hydrogen sulfide in the human gut

Responses of the microbiota to diet are highly personalized but mechanistically not well understood because many metabolic capabilities and interactions of human gut microorganisms are unknown. Here we show that sulfoquinovose (SQ), a sulfonated monosaccharide omnipresent in green vegetables, is a selective yet relevant substrate for few but ubiquitous bacteria in the human gut. In human feces and in defined co-culture, Eubacterium rectale and Bilophila wadsworthia used recently identified pathways to cooperatively catabolize SQ with 2,3-dihydroxypropane-1-sulfonate as a transient intermediate to hydrogen sulfide (H₂S), a key intestinal metabolite with disparate effects on host health. SQ-degradation capability is encoded in almost half of E. rectale genomes but otherwise sparsely distributed among microbial species in the human intestine. However, re-analysis of fecal metatranscriptome datasets of four human cohorts showed that SQ degradation (mostly from E. rectale and Faecalibacterium prausnitzii) and H₂S production (mostly from B. wadsworthia) pathways were expressed abundantly across various health states, demonstrating that these microbial functions are core attributes of the human gut. The discovery of green-diet-derived SQ as an exclusive microbial nutrient and an additional source of H₂S in the human gut highlights the role of individual dietary compounds and organosulfur metabolism on microbial activity and has implications for precision editing of the gut microbiota by dietary and prebiotic interventions.

Methylobacterium ajmalii sp. nov., Isolated From the International Space Station

Four strains belonging to the family of Methylobacteriaceae were isolated from different locations on the International Space Station (ISS) across two consecutive flights. Of these, three were identified as Gram-negative, rod-shaped, catalase-positive, oxidase-positive, motile bacteria, designated as IF7SW-B2^T, IIF1SW-B5, and IIF4SW-B5, whereas the fourth was identified as Methylorubrum rhodesianum. The sequence similarity of these three ISS strains, designated as IF7SW-B2^T, IIF1SW-B5, and IIF4SW-B5, was <99.4% for 16S rRNA genes and <97.3% for gyrB gene, with the closest being Methylobacterium indicum SE2.11^T. Furthermore, the multi-locus sequence analysis placed these three ISS strains in the same clade of M. indicum. The average nucleotide identity (ANI) values of these three ISS strains were <93% and digital DNA-DNA hybridization (dDDH) values were <46.4% with any described Methylobacterium species. Based on the ANI and dDDH analyses, these three ISS strains were considered as novel species belonging to the genus Methylobacterium. The three ISS strains showed 100% ANI similarity and dDDH values with each other, indicating that these three ISS strains, isolated during various flights and from different locations, belong to the same species. These three ISS strains were found to grow optimally at temperatures from 25 to 30°C, pH 6.0 to 8.0, and NaCl 0 to 1%. Phenotypically, these three ISS strains resemble M. aquaticum and M. terrae since they assimilate similar sugars as sole carbon substrate when compared to other Methylobacterium species. Fatty acid analysis showed that the major fatty acid produced by the ISS strains are C_18:1−ω7c and C_18:1−ω6c. The predominant quinone was ubiquinone 10, and the major polar lipids were diphosphatidylglycerol, phosphatidylcholine, phosphatidylethanolamine, phosphatidylglycerol, and an unidentified lipid. Therefore, based on genomic, phylogenetic, biochemical, and fatty acid analyses, strains IF7SW-B2^T, IIF1SW-B5, and IIF4SW-B5, are assigned to a novel species within the genus Methylobacterium, and the name Methylobacterium ajmalii sp. nov. is proposed. The type strain is IF7SW-B2^T (NRRL B-65601^T and LMG 32165^T).

Description of Candidatus Mesopelagibacter carboxydoxydans and Candidatus Anoxipelagibacter denitrificans: Nitrate-reducing SAR11 genera that dominate mesopelagic and anoxic marine zones

The diverse and ubiquitous members of the SAR11 lineage (Alphaproteobacteria) represent up to 30-40% of the surface and mesopelagic oceanic microbial communities. However, the molecular and ecological mechanisms that differentiate closely related, yet distinct, SAR11 members that often co-occur under similar environmental conditions remain speculative. Recently, two mesopelagic and oxygen minimum zone (OMZ)-associated subclades of SAR11 (Ic and IIa.A) were described using single-cell amplified genomes (SAGs) linked to nitrate reduction in OMZs. In this current study, the collection of genomes belonging to these two subclades was expanded with thirteen new metagenome-assembled genomes (MAGs), thus providing a more detailed phylogenetic and functional characterization of these subclades. Gene content-based predictions of metabolic functions revealed similarities in central carbon metabolism between subclades Ic and IIa.A and surface SAR11 clades, with small variations in central pathways. These variations included more versatile sulfur assimilation pathways, as well as a previously predicted capacity for nitrate reduction that conferred unique versatility on mesopelagic-adapted clades compared to their surface counterparts. Finally, consistent with previously reported abundances of carbon monoxide (CO) in surface and mesopelagic waters, subclades Ia (surface) and Ic (mesopelagic) have the genetic potential to oxidize carbon monoxide (CO), presumably taking advantage of this abundant compound as an electron donor. Based on genomic analyses, environmental distribution and metabolic reconstruction, we propose two new SAR11 genera, Ca. Mesopelagibacter carboxydoxydans (subclade Ic) and Ca. Anoxipelagibacter denitrificans (subclade IIa.A), which represent members of the mesopelagic and OMZ-adapted SAR11 clades.

Reference-guided metagenomics reveals genome-level evidence of potential microbial transmission from the ISS environment to an astronaut's microbiome

Monitoring microbial communities aboard the International Space Station (ISS) is essential to maintaining astronaut health and the integrity of life-support systems. Using assembled genomes of ISS-derived microbial isolates as references, recruiting metagenomic reads from an astronaut's nasal microbiome revealed no recruitment to a Staphylococcus aureus isolate from samples before launch, yet systematic recruitment across the genome when sampled after 3 months aboard the ISS, with a median percent identity of 100%. This suggests that either a highly similar S. aureus population colonized the astronaut's nasal microbiome while the astronaut was aboard the ISS or that it may have been below detection before spaceflight, instead supporting a shift in community composition. This work highlights the value in generating genomic libraries of microbes from built-environments such as the ISS and demonstrates one way such data can be integrated with metagenomics to facilitate the tracking and monitoring of astronaut microbiomes and health.

•

Understanding built-environment microbiomes is critical for human space exploration

•

Reference-guided metagenomics is a powerful approach for monitoring microbiomes

•

We show potential microbial colonization of an astronaut microbiome while aboard the ISS

Arenibacter amylolyticus sp. nov., an amylase-producing bacterium of the family Flavobacteriaceae isolated from marine water in India

A novel Gram-stain-negative, curved rod-shaped, 0.5-0.7 µm wide and 3.0-10.0 µm long, non-motile bacterium, designated strain AK53^T, was isolated from a 5 m depth water sample collected from the Bay of Bengal, Visakhapatnam, India. Colonies on marine agar were circular, small, dark orange, shiny, smooth, translucent, flat, with an entire margin. The major fatty acids included iso-C_15 : 0, iso-C_15 : 0 3OH, anteiso-C_15 : 0, iso-C_15 : 1 G, iso-C_17 : 0 3OH and summed feature 3 (C_16 : 1  ω7c and/or C_16 : 1  ω6c and/or iso-C_15 : 0-2OH). Polar lipids included phosphatidylethanolamine and five unidentified lipids. The DNA G+C content of the strain AK53^T was found to be 40.8 mol%. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain AK53^T was closely related to Arenibacter latericius KMM 426^T and Arenibacter certesii KMM3941^T (pair-wise sequence similarity of 99.17 and 98.89 %, respectively), forming a distinct branch within the genus Arenibacter and clustering with A. latericius. Strain AK53^T shared average nucleotide identity (ANIb, based on blast) of 78.07 and 77.44 % with A. latericius JCM 13508^T and A. certesii JCM 13507^T, respectively. Based on the observed phenotypic, chemotaxonomic characteristics and phylogenetic analysis, strain AK53^T is described in this study as representing a novel species in the genus Arenibacter, for which the name Arenibacter amylolyticus sp. nov. is proposed. The type strain of Arenibacter amylolyticus is AK53^T (=MTCC 12004^T= JCM 19206^T=KCTC 62553^T).

Integrated approach of whole-genome analysis, toxicological evaluation and life cycle assessment for pyrene biodegradation by a psychrophilic strain, Shewanella sp. ISTPL2

Polycyclic aromatic hydrocarbons (PAHs) such as pyrene are universal contaminants existing in the environment which have known cancer-causing and mutagenic characteristics. A psychrophilic bacterial strain Shewanella sp. ISTPL2 was isolated from the sediment sample collected from the Pangong lake, Jammu & Kashmir, India. In our previous study, the pyrene degradation potential of the ISTPL2 strain was studied in both mineral salt media as well as in soil artificially spiked with different concentrations of pyrene. Whole-genome sequencing of ISTPL2 strain in the current study highlighted the key genes of pyrene metabolism, including alcohol dehydrogenase and ring hydroxylating dioxygenase alpha-subunit. Pyrene cytotoxicity was evaluated on HepG2, a human hepato-carcinoma cell line. The cytotoxicity of the organic extract decreased with the increasing duration of bacterial treatment. To develop a more sustainable biodegradation approach, the potential impacts were evaluated for human health and ecosystem using life-cycle assessment (LCA) following the ReCiPe methodology for the considered PAH. The results implemented that global warming potential (GWP) had the highest impact, whereas both ecotoxicity and human toxicity had least from this study.

Metagenomic Profile of Microbial Communities in a Drinking Water Storage Tank Sediment after Sequential Exposure to Monochloramine, Free Chlorine, and Monochloramine

Sediment accumulation in drinking water storage facilities may lead to water quality degradation, including biological growth and disinfectant decay. The current research evaluated the microbiome present in a sediment after sequential exposure to monochloramine, free chlorine, and monochloramine. Chemical profiles within the sediment based on microelectrodes showed evidence of nitrification, and monochloramine slowly penetrated the sediment but was not measurable at lower depths. A metagenomic approach was used to characterize the microbial communities and functional potential of top (0-1 cm) and bottom (1-2 cm) layers in sediment cores. Differential abundance analysis revealed both an enrichment and depletion associated with depth of microbial populations. We assembled 30 metagenome-assembled genomes (MAGs) representing bacterial and archaeal microorganisms. Most metabolic functions were represented in both layers, suggesting the capability of the microbiomes to respond to environmental fluctuations. However, niche-specific abundance differences were identified in biotransformation processes (e.g., nitrogen). Metagenome-level analyses indicated that nitrification and denitrification can potentially occur simultaneously in the sediments, but the exact location of their occurrence within the sediment will depend on the localized physicochemical conditions. Even though monochloramine was maintained in the bulk water there was limited penetration into the sediment, and the microbial community remained functionally diverse and active.

Differential prevalence and host-association of antimicrobial resistance traits in disinfected and non-disinfected drinking water systems

Antimicrobial resistance (AMR) in drinking water has received less attention than its counterparts in the urban water cycle. While culture-based techniques or gene-centric PCR have been used to probe the impact of treatment approaches (e.g., disinfection) on AMR in drinking water, to our knowledge there is no systematic comparison of AMR trait distribution and prevalence between disinfected and disinfectant residual-free drinking water systems. We used metagenomics to assess the associations between disinfectant residuals and AMR prevalence and its host association in full-scale drinking water distribution systems (DWDSs) with and without disinfectant residuals. While the differences in AMR profiles between DWDSs were associated with the presence or absence of disinfectant, they were also associated with overall water chemistry and more importantly with microbial community structure. AMR genes and mechanisms differentially abundant in disinfected systems were primarily associated with nontuberculous mycobacteria (NTM). Finally, evaluation of metagenome assembled genomes (MAGs) also suggests that NTM possessing AMR genes conferring intrinsic resistance to key antibiotics were prevalent in disinfected systems, whereas such NTM genomes were not detected in disinfectant residual free DWDSs. Altogether, our findings provide insights into the drinking water resistome and its association with potential opportunistic pathogens, particularly in systems with disinfectant residual.

High-Quality Genome-Scale Models From Error-Prone, Long-Read Assemblies

Advances in nanopore-based sequencing techniques have enabled rapid characterization of genomes and transcriptomes. An emerging application of this sequencing technology is point-of-care characterization of pathogenic bacteria. However, genome assessments alone are unable to provide a complete understanding of the pathogenic phenotype. Genome-scale metabolic reconstruction and analysis is a bottom-up Systems Biology technique that has elucidated the phenotypic nuances of antimicrobial resistant (AMR) bacteria and other human pathogens. Combining these genome-scale models (GEMs) with point-of-care nanopore sequencing is a promising strategy for combating the emerging health challenge of AMR pathogens. However, the sequencing errors inherent to the nanopore technique may negatively affect the quality, and therefore the utility, of GEMs reconstructed from nanopore assemblies. Here we describe and validate a workflow for rapid construction of GEMs from nanopore (MinION) derived assemblies. Benchmarking the pipeline against a high-quality reference GEM of Escherichia coli K-12 resulted in nanopore-derived models that were >99% complete even at sequencing depths of less than 10× coverage. Applying the pipeline to clinical isolates of pathogenic bacteria resulted in strain-specific GEMs that identified canonical AMR genome content and enabled simulations of strain-specific microbial growth. Additionally, we show that treating the sequencing run as a mock metagenome did not degrade the quality of models derived from metagenome assemblies. Taken together, this study demonstrates that combining nanopore sequencing with GEM construction pipelines enables rapid, in situ characterization of microbial metabolism.

The saccharibacterium TM7x elicits differential responses across its host range

Host range is a fundamental component of symbiotic interactions, yet it remains poorly characterized for the prevalent yet enigmatic subcategory of bacteria/bacteria symbioses. The recently characterized obligate bacterial epibiont Candidatus Nanosynbacter lyticus TM7x with its bacterial host Actinomyces odontolyticus XH001 offers an ideal system to study such a novel relationship. In this study, the host range of TM7x was investigated by coculturing TM7x with various related Actinomyces strains and characterizing their growth dynamics from initial infection through subsequent co-passages. Of the twenty-seven tested Actinomyces, thirteen strains, including XH001, could host TM7x, and further classified into "permissive" and "nonpermissive" based on their varying initial responses to TM7x. Ten permissive strains exhibited growth/crash/recovery phases following TM7x infection, with crash timing and extent dependent on initial TM7x dosage. Meanwhile, three nonpermissive strains hosted TM7x without a growth-crash phase despite high TM7x dosage. The physical association of TM7x with all hosts, including nonpermissive strains, was confirmed by microscopy. Comparative genomic analyses revealed distinguishing genomic features between permissive and nonpermissive hosts. Our results expand the concept of host range beyond a binary to a wider spectrum, and the varying susceptibility of Actinomyces strains to TM7x underscores how small genetic differences between hosts can underly divergent selective trajectories.

Metabolic Contributions of an Alphaproteobacterial Endosymbiont in the Apicomplexan Cardiosporidium cionae

Apicomplexa is a diverse protistan phylum composed almost exclusively of metazoan-infecting parasites, including the causative agents of malaria, cryptosporidiosis, and toxoplasmosis. A single apicomplexan genus, Nephromyces, was described in 2010 as a mutualist partner to its tunicate host. Here we present genomic and transcriptomic data from the parasitic sister species to this mutualist, Cardiosporidium cionae, and its associated bacterial endosymbiont. Cardiosporidium cionae and Nephromyces both infect tunicate hosts, localize to similar organs within these hosts, and maintain bacterial endosymbionts. Though many other protists are known to harbor bacterial endosymbionts, these associations are completely unknown in Apicomplexa outside of the Nephromycidae clade. Our data indicate that a vertically transmitted α-proteobacteria has been retained in each lineage since Nephromyces and Cardiosporidium diverged. This α-proteobacterial endosymbiont has highly reduced metabolic capabilities, but contributes the essential amino acid lysine, and essential cofactor lipoic acid to C. cionae. This partnership likely reduces resource competition with the tunicate host. However, our data indicate that the contribution of the single α-proteobacterial endosymbiont in C. cionae is minimal compared to the three taxa of endosymbionts present in the Nephromyces system, and is a potential explanation for the virulence disparity between these lineages.

Complete Genome Sequence of Microbacterium sp. Strain Nx66, Isolated from Waters Contaminated with Petrochemicals in El Saf-Saf Valley, Algeria

Microbacterium sp. strain Nx66 was isolated from waters contaminated by petrochemical effluents collected in Algeria. Its genome was sequenced using Illumina MiSeq (2 × 150-bp read pairs) and Oxford Nanopore (long reads) technologies and was assembled using Unicycler. It is composed of one chromosome of 3.42 Mb and one plasmid of 34.22 kb.

Microbacterium sp. strain Nx66 was isolated from waters contaminated by petrochemical effluents collected in Algeria. Its genome was sequenced using Illumina MiSeq (2 × 150-bp read pairs) and Oxford Nanopore (long reads) technologies and was assembled using Unicycler. It is composed of one chromosome of 3.42 Mb and one plasmid of 34.22 kb.

Anoxygenic photosynthesis and iron-sulfur metabolic potential of Chlorobia populations from seasonally anoxic Boreal Shield lakes

Aquatic environments with high levels of dissolved ferrous iron and low levels of sulfate serve as an important systems for exploring biogeochemical processes relevant to the early Earth. Boreal Shield lakes, which number in the tens of millions globally, commonly develop seasonally anoxic waters that become iron rich and sulfate poor, yet the iron-sulfur microbiology of these systems has been poorly examined. Here we use genome-resolved metagenomics and enrichment cultivation to explore the metabolic diversity and ecology of anoxygenic photosynthesis and iron/sulfur cycling in the anoxic water columns of three Boreal Shield lakes. We recovered four high-completeness and low-contamination draft genome bins assigned to the class Chlorobia (formerly phylum Chlorobi) from environmental metagenome data and enriched two novel sulfide-oxidizing species, also from the Chlorobia. The sequenced genomes of both enriched species, including the novel "Candidatus Chlorobium canadense", encoded the cyc2 gene that is associated with photoferrotrophy among cultured Chlorobia members, along with genes for phototrophic sulfide oxidation. One environmental genome bin also encoded cyc2. Despite the presence of cyc2 in the corresponding draft genome, we were unable to induce photoferrotrophy in "Ca. Chlorobium canadense". Genomic potential for phototrophic sulfide oxidation was more commonly detected than cyc2 among environmental genome bins of Chlorobia, and metagenome and cultivation data suggested the potential for cryptic sulfur cycling to fuel sulfide-based growth. Overall, our results provide an important basis for further probing the functional role of cyc2 and indicate that anoxygenic photoautotrophs in Boreal Shield lakes could have underexplored photophysiology pertinent to understanding Earth's early microbial communities.

Genetic examination of the marine bacterium Pseudoalteromonas luteoviolacea and effects of its metamorphosis-inducing factors

Pseudoalteromonas luteoviolacea is a globally distributed marine bacterium that stimulates the metamorphosis of marine animal larvae, an important bacteria-animal interaction that can promote the recruitment of animals to benthic ecosystems. Recently, different P. luteoviolacea isolates have been shown to produce two stimulatory factors that can induce tubeworm and coral metamorphosis; Metamorphosis-Associated Contractile structures (MACs) and tetrabromopyrrole (TBP) respectively. However, it remains unclear what proportion of P. luteoviolacea isolates possess the genes encoding MACs, and what phenotypic effect MACs and TBP have on other larval species. Here, we show that 9 of 19 sequenced P. luteoviolacea genomes genetically encode both MACs and TBP. While P. luteoviolacea biofilms producing MACs stimulate the metamorphosis of the tubeworm Hydroides elegans, TBP biosynthesis genes had no effect under the conditions tested. Although MACs are lethal to larvae of the cnidarian Hydractinia symbiologicarpus, P. luteoviolacea mutants unable to produce MACs are capable of stimulating metamorphosis. Our findings reveal a hidden complexity of interactions between a single bacterial species, the factors it produces and two species of larvae belonging to different phyla.

Comparative Genome Analysis of Bifidobacterium longum subsp. infantis Strains Reveals Variation in Human Milk Oligosaccharide Utilization Genes among Commercial Probiotics

Dysbiosis is associated with acute and long-term consequences for neonates. Probiotics can be effective in limiting the growth of bacteria associated with dysbiosis and promoting the healthy development of the infant microbiome. Given its adaptation to the infant gut, and promising data from animal and in vitro models, Bifidobacterium longum subsp. infantis is an attractive candidate for use in infant probiotics. However, strain-level differences in the ability of commercialized strains to utilize human milk oligosaccharides (HMOs) may have implications in the performance of strains in the infant gut. In this study, we characterized twelve B. infantis probiotic strains and identified two main variants in one of the HMO utilization gene clusters. Some strains possessed the full repertoire of HMO utilization genes (H5-positive strains), while H5-negative strains lack an ABC-type transporter known to bind core HMO structures. H5-positive strains achieved significantly superior growth on lacto-N-tetraose and lacto-N-neotetraose. In vitro, H5-positive strains had a significant fitness advantage over H5-negative strains, which was also observed in vivo in breastfed infants. This work provides evidence of the functional implications of genetic differences among B. infantis strains and highlights that genotype and HMO utilization phenotype should be considered when selecting a strain for probiotic use in infants.

Unlocking the Genomic Taxonomy of the Prochlorococcus Collective

Prochlorococcus is the most abundant photosynthetic prokaryote on our planet. The extensive ecological literature on the Prochlorococcus collective (PC) is based on the assumption that it comprises one single genus comprising the species Prochlorococcus marinus, containing itself a collective of ecotypes. Ecologists adopt the distributed genome hypothesis of an open pan-genome to explain the observed genomic diversity and evolution patterns of the ecotypes within PC. Novel genomic data for the PC prompted us to revisit this group, applying the current methods used in genomic taxonomy. As a result, we were able to distinguish the five genera: Prochlorococcus, Eurycolium, Prolificoccus, Thaumococcus, and Riococcus. The novel genera have distinct genomic and ecological attributes.

Intra-species variation within Lactobacillus rhamnosus correlates to beneficial or harmful outcomes: lessons from the oral cavity

Background

The origin of most of the Lactobacillus rhamnosus genome sequences lodged in NCBI can be traced to food and faecal isolates followed by blood and tissue sites but with minimal representation from oral and vaginal isolates. However, on the L. rhamnosus phylogenetic tree no apparent clade is linked to the origin of isolation or to the relevant clinical source, except for a distinct clade exclusively shared by L. rhamnosus isolates from early stages of dental pulp infection (LRHMDP2 and LRHMDP3) and from bronchoalveolar lavage (699_LRHA and 708_LRHA) from a critical care patient. These L. rhamnosus strains, LRHMDP2, LRHMDP3, 699_LRHA and 708_LRHA isolated from different continents, display closest genome neighbour gapped identity of 99.95%. The aim of this study was to define a potentially unique complement of genes of clinical relevance shared between these L. rhamnosus clinical isolates in comparison to probiotic L. rhamnosus strains.

Results

In this analysis we used orthologous protein identification tools such as ProteinOrtho followed by tblastn alignments to identify a novel tyrosine protein phosphatase (wzb)-tyrosine-protein kinase modulator EpsC (wzd)- synteny exopolysaccharide (EPS) cluster. This EPS cluster was specifically conserved in a clade of 5 clinical isolates containing the four L. rhamnosus clinical isolates noted above and Lactobacillus spp. HMSC077C11, a clinical isolate from a neck abscess. The EPS cluster was shared with only two other strains, L. rhamnosus BPL5 and BPL15, which formed a distant clade on the L. rhamnosus phylogenetic tree, with a closest genome neighbour gapped identity of 97.51% with L. rhamnosus LRHMDP2 and LRHMDP3.

Exclusivity of this EPS cluster (from those identified before) was defined by five EPS genes, which were specifically conserved between the clade of 5 clinical isolates and L. rhamnosus BPL5 and BPL15 when compared to the remaining L. rhamnosus strains. Comparative genome analysis between the clade of 5 clinical isolates and L. rhamnosus BPL5 and BPL15 showed a set of 58 potentially unique genes characteristic of the clade of 5.

Conclusion

The potentially unique functional protein orthologs associated with the clade of 5 clinical isolates may provide understanding of fitness under selective pressure.

Benchmarking microbial growth rate predictions from metagenomes

Growth rates are central to understanding microbial interactions and community dynamics. Metagenomic growth estimators have been developed, specifically codon usage bias (CUB) for maximum growth rates and “peak-to-trough ratio” (PTR) for in situ rates. Both were originally tested with pure cultures, but natural populations are more heterogeneous, especially in individual cell histories pertinent to PTR. To test these methods, we compared predictors with observed growth rates of freshly collected marine prokaryotes in unamended seawater. We prefiltered and diluted samples to remove grazers and greatly reduce virus infection, so net growth approximated gross growth. We sampled over 44 h for abundances and metagenomes, generating 101 metagenome-assembled genomes (MAGs), including Actinobacteria, Verrucomicrobia, SAR406, MGII archaea, etc. We tracked each MAG population by cell-abundance-normalized read recruitment, finding growth rates of 0 to 5.99 per day, the first reported rates for several groups, and used these rates as benchmarks. PTR, calculated by three methods, rarely correlated to growth (r ~−0.26–0.08), except for rapidly growing γ-Proteobacteria (r ~0.63–0.92), while CUB correlated moderately well to observed maximum growth rates (r = 0.57). This suggests that current PTR approaches poorly predict actual growth of most marine bacterial populations, but maximum growth rates can be approximated from genomic characteristics.

CeMbio - The Caenorhabditis elegans Microbiome Resource

The study of microbiomes by sequencing has revealed a plethora of correlations between microbial community composition and various life-history characteristics of the corresponding host species. However, inferring causation from correlation is often hampered by the sheer compositional complexity of microbiomes, even in simple organisms. Synthetic communities offer an effective approach to infer cause-effect relationships in host-microbiome systems. Yet the available communities suffer from several drawbacks, such as artificial (thus non-natural) choice of microbes, microbe-host mismatch (e.g., human microbes in gnotobiotic mice), or hosts lacking genetic tractability. Here we introduce CeMbio, a simplified natural Caenorhabditis elegans microbiota derived from our previous meta-analysis of the natural microbiome of this nematode. The CeMbio resource is amenable to all strengths of the C. elegans model system, strains included are readily culturable, they all colonize the worm gut individually, and comprise a robust community that distinctly affects nematode life-history. Several tools have additionally been developed for the CeMbio strains, including diagnostic PCR primers, completely sequenced genomes, and metabolic network models. With CeMbio, we provide a versatile resource and toolbox for the in-depth dissection of naturally relevant host-microbiome interactions in C. elegans.

A new genomic taxonomy system for the Synechococcus collective

Cyanobacteria of the genus Synechococcus are major contributors to global primary productivity and are found in a wide range of aquatic ecosystems. This Synechococcus collective (SC) is metabolically diverse, with some lineages thriving in polar and nutrient-rich locations and others in tropical or riverine waters. Although many studies have discussed the ecology and evolution of the SC, there is a paucity of knowledge on its taxonomic structure. Thus, we present a new taxonomic classification framework for the SC based on recent advances in microbial genomic taxonomy. Phylogenomic analyses of 1085 cyanobacterial genomes demonstrate that organisms classified as Synechococcus are polyphyletic at the order rank. The SC is classified into 15 genera, which are placed into five distinct orders within the phylum Cyanobacteria: (i) Synechococcales (Cyanobium, Inmanicoccus, Lacustricoccus gen. Nov., Parasynechococcus, Pseudosynechococcus, Regnicoccus, Synechospongium gen. nov., Synechococcus and Vulcanococcus); (ii) Cyanobacteriales (Limnothrix); (iii) Leptococcales (Brevicoccus and Leptococcus); (iv) Thermosynechococcales (Stenotopis and Thermosynechococcus) and (v) Neosynechococcales (Neosynechococcus). The newly proposed classification is consistent with habitat distribution patterns (seawater, freshwater, brackish and thermal environments) and reflects the ecological and evolutionary relationships of the SC.

Draft Genome Sequences of Bacillus glennii V44-8, Bacillus saganii V47-23a, Bacillus sp. Strain V59.32b, Bacillus sp. Strain MER_TA_151, and Paenibacillus sp. Strain MER_111, Isolated from Cleanrooms Where the Viking and Mars Exploration Rover Spacecraft Were Assembled

We report the draft genome sequences of Bacillus glennii V44-8, Bacillus saganii V47-23a, and Bacillus sp. strain V59.32b, isolated from the Viking spacecraft assembly cleanroom, and Bacillus sp. strain MER_TA_151 and Paenibacillus sp. strain MER_111, isolated from the Mars Exploration Rover (MER) assembly cleanroom.

Single-Cell Genomics of Novel Actinobacteria With the Wood-Ljungdahl Pathway Discovered in a Serpentinizing System

Serpentinite-hosted systems represent modern-day analogs of early Earth environments. In these systems, water-rock interactions generate highly alkaline and reducing fluids that can contain hydrogen, methane, and low-molecular-weight hydrocarbons-potent reductants capable of fueling microbial metabolism. In this study, we investigated the microbiota of Hakuba Happo hot springs (∼50°C; pH∼10.5-11), located in Nagano (Japan), which are impacted by the serpentinization process. Analysis of the 16S rRNA gene amplicon sequences revealed that the bacterial community comprises Nitrospirae (47%), "Parcubacteria" (19%), Deinococcus-Thermus (16%), and Actinobacteria (9%), among others. Notably, only 57 amplicon sequence variants (ASV) were detected, and fifteen of these accounted for 90% of the amplicons. Among the abundant ASVs, an early-branching, uncultivated actinobacterial clade identified as RBG-16-55-12 in the SILVA database was detected. Ten single-cell genomes (average pairwise nucleotide identity: 0.98-1.00; estimated completeness: 33-93%; estimated genome size: ∼2.3 Mb) that affiliated with this clade were obtained. Taxonomic classification using single copy genes indicates that the genomes belong to the actinobacterial class-level clade UBA1414 in the Genome Taxonomy Database. Based on metabolic pathway predictions, these actinobacteria are anaerobes, capable of glycolysis, dissimilatory nitrate reduction and CO₂ fixation via the Wood-Ljungdahl (WL) pathway. Several other genomes within UBA1414 and two related class-level clades also encode the WL pathway, which has not yet been reported for the Actinobacteria phylum. For the Hakuba actinobacterium, the energy metabolism related to the WL pathway is likely supported by a combination of the Rnf complex, group 3b and 3d [NiFe]-hydrogenases, [FeFe]-hydrogenases, and V-type (H⁺/Na⁺ pump) ATPase. The genomes also harbor a form IV ribulose 1,5-bisphosphate carboxylase/oxygenase (RubisCO) complex, also known as a RubisCO-like protein, and contain signatures of interactions with viruses, including clustered regularly interspaced short palindromic repeat (CRISPR) regions and several phage integrases. This is the first report and detailed genome analysis of a bacterium within the Actinobacteria phylum capable of utilizing the WL pathway. The Hakuba actinobacterium is a member of the clade UBA1414/RBG-16-55-12, formerly within the group "OPB41." We propose to name this bacterium 'Candidatus Hakubanella thermoalkaliphilus.'

Precise phylogenetic analysis of microbial isolates and genomes from metagenomes using PhyloPhlAn 3.0

Microbial genomes are available at an ever-increasing pace, as cultivation and sequencing become cheaper and obtaining metagenome-assembled genomes (MAGs) becomes more effective. Phylogenetic placement methods to contextualize hundreds of thousands of genomes must thus be efficiently scalable and sensitive from closely related strains to divergent phyla. We present PhyloPhlAn 3.0, an accurate, rapid, and easy-to-use method for large-scale microbial genome characterization and phylogenetic analysis at multiple levels of resolution. PhyloPhlAn 3.0 can assign genomes from isolate sequencing or MAGs to species-level genome bins built from >230,000 publically available sequences. For individual clades of interest, it reconstructs strain-level phylogenies from among the closest species using clade-specific maximally informative markers. At the other extreme of resolution, it scales to large phylogenies comprising >17,000 microbial species. Examples including Staphylococcus aureus isolates, gut metagenomes, and meta-analyses demonstrate the ability of PhyloPhlAn 3.0 to support genomic and metagenomic analyses.

Honey-bee-associated prokaryotic viral communities reveal wide viral diversity and a profound metabolic coding potential

This study uses viral-like particle purification and subsequent unbiased genome sequencing to identify prokaryotic viruses associated with Apis mellifera. Interestingly, bacteriophages found in honey bees show a high diversity and span different viral taxa. This diversity sharply contrasts with the state-of-the-art knowledge on the relatively simple bee bacterial microbiome. The identification of multiple auxiliary metabolic genes suggests that these bacteriophages possess the coding potential to intervene in essential microbial pathways related to health and possibly also to disease. This study sheds light on a neglected part of the bee microbiota and opens avenues of in vivo research on the interaction of bacteriophages with their bacterial host, which likely has strongly underappreciated consequences on bee health.

Honey bees (Apis mellifera) produce an enormous economic value through their pollination activities and play a central role in the biodiversity of entire ecosystems. Recent efforts have revealed the substantial influence that the gut microbiota exert on bee development, food digestion, and homeostasis in general. In this study, deep sequencing was used to characterize prokaryotic viral communities associated with honey bees, which was a blind spot in research up until now. The vast majority of the prokaryotic viral populations are novel at the genus level, and most of the encoded proteins comprise unknown functions. Nevertheless, genomes of bacteriophages were predicted to infect nearly every major bee-gut bacterium, and functional annotation and auxiliary metabolic gene discovery imply the potential to influence microbial metabolism. Furthermore, undiscovered genes involved in the synthesis of secondary metabolic biosynthetic gene clusters reflect a wealth of previously untapped enzymatic resources hidden in the bee bacteriophage community.

Parallel Evolution in the Integration of a Co-obligate Aphid Symbiosis

Insects evolve dependence-often extreme-on microbes for nutrition. This includes cases in which insects harbor multiple endosymbionts that function collectively as a metabolic unit [1-5]. How do these dependences originate [6], and is there a predictable sequence of events leading to the integration of new symbionts? While co-obligate symbioses, in which hosts rely on multiple nutrient-provisioning symbionts, have evolved numerous times across sap-feeding insects, there is only one known case in aphids, involving Buchnera aphidicola and Serratia symbiotica in the Lachninae subfamily [7-9]. Here, we identify three additional independent transitions to the same co-obligate symbiosis in different aphids. Comparing recent and ancient associations allow us to investigate intermediate stages of metabolic and anatomical integration of Serratia. We find that these uniquely replicated evolutionary events support the idea that co-obligate associations initiate in a predictable manner-through parallel evolutionary processes. Specifically, we show how the repeated losses of the riboflavin and peptidoglycan pathways in Buchnera lead to dependence on Serratia. We then provide evidence of a stepwise process of symbiont integration, whereby dependence evolves first. Then, essential amino acid pathways are lost (at ∼30-60 mya), which coincides with the increased anatomical integration of the companion symbiont. Finally, we demonstrate that dependence can evolve ahead of specialized structures (e.g., bacteriocytes), and in one case with no direct nutritional basis. More generally, our results suggest the energetic costs of synthesizing nutrients may provide a unified explanation for the sequence of gene losses that occur during the evolution of co-obligate symbiosis.

Whole Genome Sequencing and Comparative Genomics of Two Nematicidal Bacillus Strains Reveals a Wide Range of Possible Virulence Factors

Bacillus firmus nematicidal bacterial strains are used to control plant parasitic nematode infestation of crops in agricultural production. Proteases are presumed to be the primary nematode virulence factors in nematicidal B. firmus degrading the nematode cuticle and other organs. We determined and compared the whole genome sequences of two nematicidal strains. Comparative genomics with a particular focus on possible virulence determinants revealed a wider range of possible virulence factors in a B. firmus isolate from a commercial bionematicide and a wild type Bacillus sp. isolate with nematicidal activity. The resulting 4.6 Mb B. firmus I-1582 and 5.3 Mb Bacillus sp. ZZV12-4809 genome assemblies contain respectively 18 and 19 homologs to nematode-virulent proteases, two nematode-virulent chitinase homologs in ZZV12-4809 and 28 and 36 secondary metabolite biosynthetic clusters, projected to encode antibiotics, small peptides, toxins and siderophores. The results of this study point to the genetic capability of B. firmus and related species for nematode virulence through a range of direct and indirect mechanisms.

Genomic and phenotypic characterization of Burkholderia isolates from the potable water system of the International Space Station

The opportunistic pathogens Burkholderia cepacia and Burkholderia contaminans, both genomovars of the Burkholderia cepacia complex (BCC), are frequently cultured from the potable water dispenser (PWD) of the International Space Station (ISS). Here, we sequenced the genomes and conducted phenotypic assays to characterize these Burkholderia isolates. All recovered isolates of the two species fall within monophyletic clades based on phylogenomic trees of conserved single-copy core genes. Within species, the ISS-derived isolates all demonstrate greater than 99% average nucleotide identity (with 95-99% of genomes aligning) and share around 90% of the identified gene clusters from a pangenomic analysis-suggesting that the two groups are each composed of highly similar genomic lineages and their members may have all stemmed from the same two founding populations. The differences that can be observed between the recovered isolates at the pangenomic level are primarily located within putative plasmids. Phenotypically, macrophage intracellularization and lysis occurred at generally similar rates between all ISS-derived isolates, as well as with their respective type-terrestrial strain references. All ISS-derived isolates exhibited antibiotic sensitivity similar to that of the terrestrial reference strains, and minimal differences between isolates were observed. With a few exceptions, biofilm formation rates were generally consistent across each species. And lastly, though isolation date does not necessarily provide any insight into how long a given isolate had been aboard the ISS, none of the assayed physiology correlated with either date of isolation or distances based on nucleotide variation. Overall, we find that while the populations of Burkholderia present in the ISS PWS each maintain virulence, they are likely are not more virulent than those that might be encountered on planet and remain susceptible to clinically used antibiotics.

Draft Genome Sequence of Cyclobacterium marinum Strain Atlantic-IS, Isolated from the Atlantic Slope off the Coast of Virginia, USA

Here, we report a draft genome sequence for Cyclobacterium marinum strain Atlantic-IS, isolated from the Atlantic slope off the coast of Virginia. The whole-genome sequence will help us understand its adaptive metabolic responses to diverse C sources in low-nutrient environments.

Isolation and Characterization of Akhmeta Virus from Wild-Caught Rodents (Apodemus spp.) in Georgia

Akhmeta virus is a unique Orthopoxvirus that was described in 2013 from the country of Georgia. This paper presents the first isolation of this virus from small mammal (Rodentia; Apodemus spp.) samples and the molecular characterization of those isolates. The identification of the virus in small mammals is an essential component to understanding the natural history of this virus and its transmission to human populations and could guide public health interventions in Georgia. Akhmeta virus genomes harbor evidence suggestive of recombination with a variety of other orthopoxviruses; this has implications for the evolution of orthopoxviruses, their ability to infect mammalian hosts, and their ability to adapt to novel host species.

In 2013, a novel orthopoxvirus was detected in skin lesions of two cattle herders from the Kakheti region of Georgia (country); this virus was named Akhmeta virus. Subsequent investigation of these cases revealed that small mammals in the area had serological evidence of orthopoxvirus infections, suggesting their involvement in the maintenance of these viruses in nature. In October 2015, we began a longitudinal study assessing the natural history of orthopoxviruses in Georgia. As part of this effort, we trapped small mammals near Akhmeta (n = 176) and Gudauri (n = 110). Here, we describe the isolation and molecular characterization of Akhmeta virus from lesion material and pooled heart and lung samples collected from five wood mice (Apodemus uralensis and Apodemus flavicollis) in these two locations. The genomes of Akhmeta virus obtained from rodents group into 2 clades: one clade represented by viruses isolated from A. uralensis samples, and one clade represented by viruses isolated from A. flavicollis samples. These genomes also display several presumptive recombination events for which gene truncation and identity have been examined.

IMPORTANCE Akhmeta virus is a unique Orthopoxvirus that was described in 2013 from the country of Georgia. This paper presents the first isolation of this virus from small mammal (Rodentia; Apodemus spp.) samples and the molecular characterization of those isolates. The identification of the virus in small mammals is an essential component to understanding the natural history of this virus and its transmission to human populations and could guide public health interventions in Georgia. Akhmeta virus genomes harbor evidence suggestive of recombination with a variety of other orthopoxviruses; this has implications for the evolution of orthopoxviruses, their ability to infect mammalian hosts, and their ability to adapt to novel host species.

Draft Genome Sequences of Five Proteobacteria Isolated from Lechuguilla Cave, New Mexico, USA, and Insights into Taxonomy and Quorum Sensing

Genomic resources remain scarce for bacteria isolated from oligotrophic caves. We sequenced the genomes of five Proteobacteria isolated from Lechuguilla Cave in New Mexico. Genome-based phylogeny indicates that each strain belongs to a distinct genus. Two Rhizobiaceae isolates possess genomic potential for the biosynthesis of acyl-homoserine lactone.

Genomic resources remain scarce for bacteria isolated from oligotrophic caves. We sequenced the genomes of five Proteobacteria isolated from Lechuguilla Cave in New Mexico. Genome-based phylogeny indicates that each strain belongs to a distinct genus. Two Rhizobiaceae isolates possess genomic potential for the biosynthesis of acyl-homoserine lactone.

Applications and Considerations of GToTree: A User-Friendly Workflow for Phylogenomics

Phylogenomics is the practice of attempting to infer evolutionary relationships at a genome-level. This is becoming a standard step in the characterization of newly recovered genomes and to direct/constrain further research; yet the process from start to finish of building a de novo phylogenomic tree that is specific to the organisms of interest can still be computationally intractable for many biologists. GToTree is a recently published user-friendly workflow for phylogenomics intended to give more researchers the capability to generate phylogenomic trees to help guide their work. This commentary describes two common applications where GToTree can be helpful and then discusses some things to consider when using the program.

Mycolicibacterium nivoides sp. nov isolated from a peat bog

A fast-growing, non-chromogenic, acid-fast-staining bacterium (DL90T) was isolated from a peat bog in northern Minnesota. On the basis of 16S rRNA gene sequence similarity (99.8 % identity with Mycolicibacterium septicum and 98 % with Mycolicibacterium peregrinum) and chemotaxonomic data (fatty acid content), strain DL90T represents a member of the genus Mycolicibacterium. Physiological tests (growth curves, biofilm formation, antibiotic sensitivity, colony morphologies and heat tolerance) and biochemical analysis (arylsulfatase activity and fatty acid profiles) distinguish DL90T from its closest relative M. septicum. Phylogenomic reconstruction of the 'Fortuitium-Vaccae' clade, digital DNA-DNA hybridization (DDH) values of 61 %, and average nucleotide identity (ANI) values of approximately 95 % indicate that DL90T is likely to be diverged from M. septicum. Thus, we propose that DL90T represents a novel species, given the name Mycolicibacterium nivoides with the type strain being isolate DL90T (=JCM 32796T=NCCB 100660T).

Survey of Toxin⁻Antitoxin Systems in Erwinia amylovora Reveals Insights into Diversity and Functional Specificity

Toxin–antitoxin (TA) systems are diverse genetic modules with demonstrated roles in plasmid stability, stress management, biofilm formation and antibiotic persistence. However, relatively little is known about their functional significance in plant pathogens. In this study we characterize type II and IV TA systems in the economically important plant pathogen Erwinia amylovora. Hidden Markov Model (HMM) and BLAST-based programs were used to predict the identity and distribution of putative TA systems among sequenced genomes of E. amylovora and other plant-associated Erwinia spp. Of six conserved TA systems tested for function from E. amylovora, three (CbtA/CbeA, ParE/RHH and Doc/PhD) were validated as functional. CbtA was toxic to E. amylovora, but not to Escherichia coli. While the E. coli homolog of CbtA elicits the formation of lemon-shaped cells upon overexpression and targets cytoskeletal proteins FtsZ and MreB, E. amylovora CbtA led to cell elongation and did not interact with these cytoskeletal proteins. Phylogenetic analysis revealed that E. amylovora CbtA belongs to a distinct clade from the CbtA of pathogenic E. coli. This study expands the repertoire of experimentally validated TA systems in plant pathogenic bacteria, and suggests that the E. amylovora homolog of CbtA is functionally distinct from that of E. coli.