Multicellular magnetotactic bacteria are genetically heterogeneous consortia with metabolically differentiated cells.

Consortia of multicellular magnetotactic bacteria (MMB) are currently the only known example of bacteria without a unicellular stage in their life cycle. Because of their recalcitrance to cultivation, most previous studies of MMB have been limited to microscopic observations. To study the biology of these unique organisms in more detail, we use multiple culture-independent approaches to analyze the genomics and physiology of MMB consortia at single-cell resolution. We separately sequenced the metagenomes of 22 individual MMB consortia, representing 8 new species, and quantified the genetic diversity within each MMB consortium. This revealed that, counter to conventional views, cells within MMB consortia are not clonal. Single consortia metagenomes were then used to reconstruct the species-specific metabolic potential and infer the physiological capabilities of MMB. To validate genomic predictions, we performed stable isotope probing (SIP) experiments and interrogated MMB consortia using fluorescence in situ hybridization (FISH) combined with nanoscale secondary ion mass spectrometry (NanoSIMS). By coupling FISH with bioorthogonal noncanonical amino acid tagging (BONCAT), we explored their in situ activity as well as variation of protein synthesis within cells. We demonstrate that MMB consortia are mixotrophic sulfate reducers and that they exhibit metabolic differentiation between individual cells, suggesting that MMB consortia are more complex than previously thought. These findings expand our understanding of MMB diversity, ecology, genomics, and physiology, as well as offer insights into the mechanisms underpinning the multicellular nature of their unique lifestyle.

Giant virus diversity and host interactions through global metagenomics.

Our current knowledge about nucleocytoplasmic large DNA viruses (NCLDVs) is largely derived from viral isolates that are co-cultivated with protists and algae. Here we reconstructed 2,074 NCLDV genomes from sampling sites across the globe by building on the rapidly increasing amount of publicly available metagenome data. This led to an 11-fold increase in phylogenetic diversity and a parallel 10-fold expansion in functional diversity. Analysis of 58,023 major capsid proteins from large and giant viruses using metagenomic data revealed the global distribution patterns and cosmopolitan nature of these viruses. The discovered viral genomes encoded a wide range of proteins with putative roles in photosynthesis and diverse substrate transport processes, indicating that host reprogramming is probably a common strategy in the NCLDVs. Furthermore, inferences of horizontal gene transfer connected viral lineages to diverse eukaryotic hosts. We anticipate that the global diversity of NCLDVs that we describe here will establish giant viruses-which are associated with most major eukaryotic lineages-as important players in ecosystems across Earth's biomes.

Novel Z-DNA binding domains in giant viruses.

Z-nucleic acid structures play vital roles in cellular processes and have implications in innate immunity due to their recognition by Zα domains containing proteins (Z-DNA/Z-RNA binding proteins, ZBPs). Although Zα domains have been identified in six proteins, including viral E3L, ORF112, and I73R, as well as, cellular ADAR1, ZBP1, and PKZ, their prevalence across living organisms remains largely unexplored. In this study, we introduce a computational approach to predict Zα domains, leading to the revelation of previously unidentified Zα domain-containing proteins in eukaryotic organisms, including non-metazoan species. Our findings encompass the discovery of new ZBPs in previously unexplored giant viruses, members of the Nucleocytoviricota phylum. Through experimental validation, we confirm the Zα functionality of select proteins, establishing their capability to induce the B-to-Z conversion. Additionally, we identify Zα-like domains within bacterial proteins. While these domains share certain features with Zα domains, they lack the ability to bind to Z-nucleic acids or facilitate the B-to-Z DNA conversion. Our findings significantly expand the ZBP family across a wide spectrum of organisms and raise intriguing questions about the evolutionary origins of Zα-containing proteins. Moreover, our study offers fresh perspectives on the functional significance of Zα domains in virus sensing and innate immunity and opens avenues for exploring hitherto undiscovered functions of ZBPs.

IMG/VR v4: an expanded database of uncultivated virus genomes within a framework of extensive functional, taxonomic, and ecological metadata.

Viruses are widely recognized as critical members of all microbiomes. Metagenomics enables large-scale exploration of the global virosphere, progressively revealing the extensive genomic diversity of viruses on Earth and highlighting the myriad of ways by which viruses impact biological processes. IMG/VR provides access to the largest collection of viral sequences obtained from (meta)genomes, along with functional annotation and rich metadata. A web interface enables users to efficiently browse and search viruses based on genome features and/or sequence similarity. Here, we present the fourth version of IMG/VR, composed of >15 million virus genomes and genome fragments, a ≈6-fold increase in size compared to the previous version. These clustered into 8.7 million viral operational taxonomic units, including 231 408 with at least one high-quality representative. Viral sequences in IMG/VR are now systematically identified from genomes, metagenomes, and metatranscriptomes using a new detection approach (geNomad), and IMG standard annotation are complemented with genome quality estimation using CheckV, taxonomic classification reflecting the latest taxonomic standards, and microbial host taxonomy prediction. IMG/VR v4 is available at https://img.jgi.doe.gov/vr, and the underlying data are available to download at https://genome.jgi.doe.gov/portal/IMG_VR.

IMG/VR v3: an integrated ecological and evolutionary framework for interrogating genomes of uncultivated viruses.

Viruses are integral components of all ecosystems and microbiomes on Earth. Through pervasive infections of their cellular hosts, viruses can reshape microbial community structure and drive global nutrient cycling. Over the past decade, viral sequences identified from genomes and metagenomes have provided an unprecedented view of viral genome diversity in nature. Since 2016, the IMG/VR database has provided access to the largest collection of viral sequences obtained from (meta)genomes. Here, we present the third version of IMG/VR, composed of 18 373 cultivated and 2 314 329 uncultivated viral genomes (UViGs), nearly tripling the total number of sequences compared to the previous version. These clustered into 935 362 viral Operational Taxonomic Units (vOTUs), including 188 930 with two or more members. UViGs in IMG/VR are now reported as single viral contigs, integrated proviruses or genome bins, and are annotated with a new standardized pipeline including genome quality estimation using CheckV, taxonomic classification reflecting the latest ICTV update, and expanded host taxonomy prediction. The new IMG/VR interface enables users to efficiently browse, search, and select UViGs based on genome features and/or sequence similarity. IMG/VR v3 is available at https://img.jgi.doe.gov/vr, and the underlying data are available to download at https://genome.jgi.doe.gov/portal/IMG_VR.

Yaravirus: A novel 80-nm virus infecting Acanthamoeba castellanii.

Here we report the discovery of Yaravirus, a lineage of amoebal virus with a puzzling origin and evolution. Yaravirus presents 80-nm-sized particles and a 44,924-bp dsDNA genome encoding for 74 predicted proteins. Yaravirus genome annotation showed that none of its genes matched with sequences of known organisms at the nucleotide level; at the amino acid level, six predicted proteins had distant matches in the nr database. Complimentary prediction of three-dimensional structures indicated possible function of 17 proteins in total. Furthermore, we were not able to retrieve viral genomes closely related to Yaravirus in 8,535 publicly available metagenomes spanning diverse habitats around the globe. The Yaravirus genome also contained six types of tRNAs that did not match commonly used codons. Proteomics revealed that Yaravirus particles contain 26 viral proteins, one of which potentially representing a divergent major capsid protein (MCP) with a predicted double jelly-roll domain. Structure-guided phylogeny of MCP suggests that Yaravirus groups together with the MCPs of Pleurochrysis endemic viruses. Yaravirus expands our knowledge of the diversity of DNA viruses. The phylogenetic distance between Yaravirus and all other viruses highlights our still preliminary assessment of the genomic diversity of eukaryotic viruses, reinforcing the need for the isolation of new viruses of protists.

Hydrogenotrophic methanogenesis in archaeal phylum Verstraetearchaeota reveals the shared ancestry of all methanogens.

Methanogenic archaea are major contributors to the global carbon cycle and were long thought to belong exclusively to the euryarchaeal phylum. Discovery of the methanogenesis gene cluster methyl-coenzyme M reductase (Mcr) in the Bathyarchaeota, and thereafter the Verstraetearchaeota, led to a paradigm shift, pushing back the evolutionary origin of methanogenesis to predate that of the Euryarchaeota. The methylotrophic methanogenesis found in the non-Euryarchaota distinguished itself from the predominantly hydrogenotrophic methanogens found in euryarchaeal orders as the former do not couple methanogenesis to carbon fixation through the reductive acetyl-CoA [Wood-Ljungdahl pathway (WLP)], which was interpreted as evidence for independent evolution of the two methanogenesis pathways. Here, we report the discovery of a complete and divergent hydrogenotrophic methanogenesis pathway in a thermophilic order of the Verstraetearchaeota, which we have named Candidatus Methanohydrogenales, as well as the presence of the WLP in the crenarchaeal order Desulfurococcales. Our findings support the ancient origin of hydrogenotrophic methanogenesis, suggest that methylotrophic methanogenesis might be a later adaptation of specific orders, and provide insight into how the transition from hydrogenotrophic to methylotrophic methanogenesis might have occurred.

The trajectory of microbial single-cell sequencing.

Over the past decade, it has become nearly routine to sequence genomes of individual microbial cells directly isolated from environmental samples ranging from deep-sea hydrothermal vents to insect guts, providing a powerful complement to shotgun metagenomics in microbial community studies. In this review, we address the technical aspects and challenges of single-cell genome sequencing and discuss some of the scientific endeavors that it has enabled. Specifically, we highlight newly added leaves and branches in the genomic tree of bacterial and archaeal life and illustrate the unique and exciting advantages that single-cell genomics offers over metagenomics, both now and in the near future.

Coexistence of novel gammaproteobacterial and Arsenophonus symbionts in the scale insect Greenisca brachypodii (Hemiptera, Coccomorpha: Eriococcidae).

Scale insects are commonly associated with obligate, intracellular microorganisms which play important roles in complementing their hosts with essential nutrients. Here we characterized the symbiotic system of Greenisca brachypodii, a member of the family Eriococcidae. Histological and ultrastructural analyses have indicated that G. brachypodii is stably associated with coccoid and rod-shaped bacteria. Phylogenetic analyses have revealed that the coccoid bacteria represent a sister group to the secondary symbiont of the mealybug Melanococcus albizziae, whereas the rod-shaped symbionts are close relatives of Arsenophonus symbionts in insects - to our knowledge, this is the first report of the presence of Arsenophonus bacterium in scale insects. As a comparison of 16S and 23S rRNA genes sequences of the G. brachypodii coccoid symbiont with other gammaprotebacterial sequences showed only low similarity (∼90%), we propose the name 'Candidatus Kotejella greeniscae' for its tentative classification. Both symbionts are transovarially transmitted from one generation to the next. The infection takes place in the neck region of the ovariole. The bacteria migrate between follicular cells, as well as through the cytoplasm of those cells to the perivitelline space, where they form a characteristic 'symbiont ball'. Our findings provide evidence for a polyphyletic origin of symbionts of Eriococcidae.

Genome-based evolutionary history of Pseudomonas spp.

Pseudomonas is a large and diverse genus of Gammaproteobacteria. To provide a framework for discovery of evolutionary and taxonomic relationships of these bacteria, we compared the genomes of type strains of 163 species and 3 additional subspecies of Pseudomonas, including 118 genomes sequenced herein. A maximum likelihood phylogeny of the 166 type strains based on protein sequences of 100 single-copy orthologous genes revealed thirteen groups of Pseudomonas, composed of two to sixty three species each. Pairwise average nucleotide identities and alignment fractions were calculated for the data set of the 166 type strains and 1224 genomes of Pseudomonas available in public databases. Results revealed that 394 of the 1224 genomes were distinct from any type strain, suggesting that the type strains represent only a fraction of the genomic diversity of the genus. The core genome of Pseudomonas was determined to contain 794 genes conferring primarily housekeeping functions. The results of this study provide a phylogenetic framework for future studies aiming to resolve the classification and phylogenetic relationships, identify new gene functions and phenotypes, and explore the ecological and metabolic potential of the Pseudomonas spp.

A Rickettsiales symbiont of amoebae with ancient features.

The Rickettsiae comprise intracellular bacterial symbionts and pathogens infecting diverse eukaryotes. Here, we provide a detailed characterization of 'Candidatus Jidaibacter acanthamoeba', a rickettsial symbiont of Acanthamoeba. The bacterium establishes the infection in its amoeba host within 2 h where it replicates within vacuoles. Higher bacterial loads and accelerated spread of infection at elevated temperatures were observed. The infection had a negative impact on host growth rate, although no increased levels of host cell lysis were seen. Phylogenomic analysis identified this bacterium as member of the Midichloriaceae. Its 2.4 Mb genome represents the largest among Rickettsiales and is characterized by a moderate degree of pseudogenization and a high coding density. We found an unusually large number of genes encoding proteins with eukaryotic-like domains such as ankyrins, leucine-rich repeats and tetratricopeptide repeats, which likely function in host interaction. There are a total of three divergent, independently acquired type IV secretion systems, and 35 flagellar genes representing the most complete set found in an obligate intracellular Alphaproteobacterium. The deeply branching phylogenetic position of 'Candidatus Jidaibacter acanthamoeba' together with its ancient features place it closely to the rickettsial ancestor and helps to better understand the transition from a free-living to an intracellular lifestyle.

Chlamydial seasonal dynamics and isolation of 'Candidatus†Neptunochlamydia vexilliferae' from a Tyrrhenian coastal lake.

The Chlamydiae are a phylum of obligate intracellular bacteria comprising important human and animal pathogens, yet their occurrence in the environment, their phylogenetic diversity and their host range has been largely underestimated. We investigated the seasonality of environmental chlamydiae in a Tyrrhenian coastal lake. By catalysed reporter deposition fluorescence in situ hybridization, we quantified the small planktonic cells and detected a peak in the abundance of environmental chlamydiae in early autumn with up to 5.9 × 10(4) cells ml(-1) . Super-resolution microscopy improved the visualization and quantification of these bacteria and enabled the detection of pleomorphic chlamydial cells in their protist host directly in an environmental sample. To isolate environmental chlamydiae together with their host, we applied a high-throughput limited dilution approach and successfully recovered a Vexillifera sp., strain harbouring chlamydiae (93% 16S rRNA sequence identity to Simkania negevensis), tentatively named 'Candidatus Neptunochlamydia vexilliferae'. Transmission electron microscopy in combination with fluorescence in situ hybridization was used to prove the intracellular location of these bacteria representing the first strain of marine chlamydiae stably maintained alongside with their host in a laboratory culture. Taken together, this study contributes to a better understanding of the distribution and diversity of environmental chlamydiae in previously neglected marine environments.

Prediction of microbial phenotypes based on comparative genomics.

The accessibility of almost complete genome sequences of uncultivable microbial species from metagenomes necessitates computational methods predicting microbial phenotypes solely based on genomic data. Here we investigate how comparative genomics can be utilized for the prediction of microbial phenotypes. The PICA framework facilitates application and comparison of different machine learning techniques for phenotypic trait prediction. We have improved and extended PICA's support vector machine plug-in and suggest its applicability to large-scale genome databases and incomplete genome sequences. We have demonstrated the stability of the predictive power for phenotypic traits, not perturbed by the rapid growth of genome databases. A new software tool facilitates the in-depth analysis of phenotype models, which associate expected and unexpected protein functions with particular traits. Most of the traits can be reliably predicted in only 60-70% complete genomes. We have established a new phenotypic model that predicts intracellular microorganisms. Thereby we could demonstrate that also independently evolved phenotypic traits, characterized by genome reduction, can be reliably predicted based on comparative genomics. Our results suggest that the extended PICA framework can be used to automatically annotate phenotypes in near-complete microbial genome sequences, as generated in large numbers in current metagenomics studies.

Adaptation strategies of giant viruses to low-temperature marine ecosystems.

Microbes in marine ecosystems have evolved their gene content to thrive successfully in the cold. Although this process has been reasonably well studied in bacteria and selected eukaryotes, less is known about the impact of cold environments on the genomes of viruses that infect eukaryotes. Here, we analyzed cold adaptations in giant viruses (Nucleocytoviricota and Mirusviricota) from austral marine environments and compared them with their Arctic and temperate counterparts. We recovered giant virus metagenome-assembled genomes (98 Nucleocytoviricota and 12 Mirusviricota MAGs) from 61 newly sequenced metagenomes and metaviromes from sub-Antarctic Patagonian fjords and Antarctic seawater samples. When analyzing our data set alongside Antarctic and Arctic giant viruses MAGs already deposited in the Global Ocean Eukaryotic Viral (GOEV) database, we found that Antarctic and Arctic giant viruses predominantly inhabit sub-10°C environments, featuring a high proportion of unique phylotypes in each ecosystem. In contrast, giant viruses in Patagonian fjords were subject to broader temperature ranges and showed a lower degree of endemicity. However, despite differences in their distribution, giant viruses inhabiting low-temperature marine ecosystems evolved genomic cold-adaptation strategies that led to changes in genetic functions and amino acid frequencies that ultimately affect both gene content and protein structure. Such changes seem to be absent in their mesophilic counterparts. The uniqueness of these cold-adapted marine giant viruses may now be threatened by climate change, leading to a potential reduction in their biodiversity.

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake.

BACKGROUND: Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria. RESULTS: Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota. CONCLUSIONS: We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

Distinct and rich assemblages of giant viruses in Arctic and Antarctic lakes.

Giant viruses (GVs) are key players in ecosystem functioning, biogeochemistry, and eukaryotic genome evolution. GV diversity and abundance in aquatic systems can exceed that of prokaryotes, but their diversity and ecology in lakes, especially polar ones, remain poorly understood. We conducted a comprehensive survey and meta-analysis of GV diversity across 20 lakes, spanning polar to temperate regions, combining our extensive lake metagenome database from the Canadian Arctic and subarctic with publicly available datasets. Leveraging a novel GV genome identification tool, we identified 3304 GV metagenome-assembled genomes, revealing lakes as untapped GV reservoirs. Phylogenomic analysis highlighted their dispersion across all Nucleocytoviricota orders. Strong GV population endemism emerged between lakes from similar regions and biomes (Antarctic and Arctic), but a polar/temperate barrier in lacustrine GV populations and differences in their gene content could be observed. Our study establishes a robust genomic reference for future investigations into lacustrine GV ecology in fast changing polar environments.

Comparative genomics of a vertically transmitted thiotrophic bacterial ectosymbiont and its close free-living relative.

Thiotrophic symbioses between sulphur-oxidizing bacteria and various unicellular and metazoan eukaryotes are widespread in reducing marine environments. The giant colonial ciliate Zoothamnium niveum, however, is the only host of thioautotrophic symbionts that has been cultivated along with its symbiont, the vertically transmitted ectosymbiont Candidatus Thiobius zoothamnicola (short Thiobius). Because theoretical predictions posit a smaller genome in vertically transmitted endosymbionts compared to free-living relatives, we investigated whether this is true also for an ectosymbiont. We used metagenomics to recover the high-quality draft genome of this bacterial symbiont. For comparison we have also sequenced a closely related free-living cultured but not formally described strain Milos ODIII6 (short ODIII6). We then performed comparative genomics to assess the functional capabilities at gene, metabolic pathway and trait level. 16S rRNA gene trees and average amino acid identity confirmed the close phylogenetic relationship of both bacteria. Indeed, Thiobius has about a third smaller genome than its free-living relative ODIII6, with reduced metabolic capabilities and fewer functional traits. The functional capabilities of Thiobius were a subset of those of the more versatile ODIII6, which possessed additional genes for oxygen, sulphur and hydrogen utilization and for the acquisition of phosphorus illustrating features that may be adaptive for the unstable environmental conditions at hydrothermal vents. In contrast, Thiobius possesses genes potentially enabling it to utilize lactate and acetate heterotrophically, compounds that may be provided as byproducts by the host. The present study illustrates the effect of strict host-dependence of a bacterial ectosymbiont on genome evolution and host adaptation.

A giant virus infecting the amoeboflagellate Naegleria.

Giant viruses (Nucleocytoviricota) are significant lethality agents of various eukaryotic hosts. Although metagenomics indicates their ubiquitous distribution, available giant virus isolates are restricted to a very small number of protist and algal hosts. Here we report on the first viral isolate that replicates in the amoeboflagellate Naegleria. This genus comprises the notorious human pathogen Naegleria fowleri, the causative agent of the rare but fatal primary amoebic meningoencephalitis. We have elucidated the structure and infection cycle of this giant virus, Catovirus naegleriensis (a.k.a. Naegleriavirus, NiV), and show its unique adaptations to its Naegleria host using fluorescence in situ hybridization, electron microscopy, genomics, and proteomics. Naegleriavirus is only the fourth isolate of the highly diverse subfamily Klosneuvirinae, and like its relatives the NiV genome contains a large number of translation genes, but lacks transfer RNAs (tRNAs). NiV has acquired genes from its Naegleria host, which code for heat shock proteins and apoptosis inhibiting factors, presumably for host interactions. Notably, NiV infection was lethal to all Naegleria species tested, including the human pathogen N. fowleri. This study expands our experimental framework for investigating giant viruses and may help to better understand the basic biology of the human pathogen N. fowleri.

Phylogenomics and genetic analysis of solvent-producing Clostridium species.

The genus Clostridium is a large and diverse group within the Bacillota (formerly Firmicutes), whose members can encode useful complex traits such as solvent production, gas-fermentation, and lignocellulose breakdown. We describe 270 genome sequences of solventogenic clostridia from a comprehensive industrial strain collection assembled by Professor David Jones that includes 194 C. beijerinckii, 57 C. saccharobutylicum, 4 C. saccharoperbutylacetonicum, 5 C. butyricum, 7 C. acetobutylicum, and 3 C. tetanomorphum genomes. We report methods, analyses and characterization for phylogeny, key attributes, core biosynthetic genes, secondary metabolites, plasmids, prophage/CRISPR diversity, cellulosomes and quorum sensing for the 6 species. The expanded genomic data described here will facilitate engineering of solvent-producing clostridia as well as non-model microorganisms with innately desirable traits. Sequences could be applied in conventional platform biocatalysts such as yeast or Escherichia coli for enhanced chemical production. Recently, gene sequences from this collection were used to engineer Clostridium autoethanogenum, a gas-fermenting autotrophic acetogen, for continuous acetone or isopropanol production, as well as butanol, butanoic acid, hexanol and hexanoic acid production.

Solvent-Producing Clostridia Revisited.

The review provides an overview of the current status of the solvent-producing clostridia. The origin and development of industrial clostridial species, as well as the history of the industrial Acetone Butanol Ethanol fermentation process, is reexamined, and the recent resurgence of interest in the production of biobutanol is reviewed. Over 300 fully sequenced genomes for solvent-producing and closely related clostridial species are currently available in public databases. These include 270 genomes sourced from the David Jones culture collection. These genomes were allocated arbitrary DJ codes, and a conversion table to identify the species and strains has now been provided. The expanded genomic database facilitated new comparative genomic and phylogenetic analysis. A synopsis of the common features, molecular taxonomy, and phylogeny of solvent-producing clostridia and the application of comparative phylogenomics are evaluated. A survey and analysis of resident prophages in solvent-producing clostridia are discussed, and the discovery, occurrence, and role of novel R-type tailocins are reported. Prophage genomes with R-type tailocin-like features were detected in all 12 species investigated. The widespread occurrence of tailocins in Gram-negative species is well documented; this survey has indicated that they may also be widespread in clostridia.