pyCapsid: identifying dominant dynamics and quasi-rigid mechanical units in protein shells.

SUMMARY: pyCapsid is a Python package developed to facilitate the characterization of the dynamics and quasi-rigid mechanical units of protein shells and other protein complexes. The package was developed in response to the rapid increase of high-resolution structures, particularly capsids of viruses, requiring multiscale biophysical analyses. Given a protein shell, pyCapsid generates the collective vibrations of its amino-acid residues, identifies quasi-rigid mechanical regions associated with the disassembly of the structure, and maps the results back to the input proteins for interpretation. pyCapsid summarizes the main results in a report that includes publication-quality figures. AVAILABILITY AND IMPLEMENTATION: pyCapsid's source code is available under MIT License on GitHub. It is compatible with Python 3.8-3.10 and has been deployed in two leading Python package-management systems, PIP and Conda. Installation instructions and tutorials are available in the online documentation and in the pyCapsid's YouTube playlist. In addition, a cloud-based implementation of pyCapsid is available as a Google Colab notebook. pyCapsid Colab does not require installation and generates the same report and outputs as the installable version. Users can post issues regarding pyCapsid in the repository's issues section.

Viral predation pressure on coral reefs.

BACKGROUND: Predation pressure and herbivory exert cascading effects on coral reef health and stability. However, the extent of these cascading effects can vary considerably across space and time. This variability is likely a result of the complex interactions between coral reefs' biotic and abiotic dimensions. A major biological component that has been poorly integrated into the reefs' trophic studies is the microbial community, despite its role in coral death and bleaching susceptibility. Viruses that infect bacteria can control microbial densities and may positively affect coral health by controlling microbialization. We hypothesize that viral predation of bacteria has analogous effects to the top-down pressure of macroorganisms on the trophic structure and reef health. RESULTS: Here, we investigated the relationships between live coral cover and viruses, bacteria, benthic algae, fish biomass, and water chemistry in 110 reefs spanning inhabited and uninhabited islands and atolls across the Pacific Ocean. Statistical learning showed that the abundance of turf algae, viruses, and bacteria, in that order, were the variables best predicting the variance in coral cover. While fish biomass was not a strong predictor of coral cover, the relationship between fish and corals became apparent when analyzed in the context of viral predation: high coral cover (> 50%) occurred on reefs with a combination of high predator fish biomass (sum of sharks and piscivores > 200 g m-2) and high virus-to-bacteria ratios (> 10), an indicator of viral predation pressure. However, these relationships were non-linear, with reefs at the higher and lower ends of the coral cover continuum displaying a narrow combination of abiotic and biotic variables, while reefs at intermediate coral cover showed a wider range of parameter combinations. CONCLUSIONS: The results presented here support the hypothesis that viral predation of bacteria is associated with high coral cover and, thus, coral health and stability. We propose that combined predation pressures from fishes and viruses control energy fluxes, inhibiting the detrimental accumulation of ecosystem energy in the microbial food web.

The landscape of lysogeny across microbial community density, diversity and energetics.

Lysogens are common at high bacterial densities, an observation that contrasts with the prevailing view of lysogeny as a low-density refugium strategy. Here, we review the mechanisms regulating lysogeny in complex communities and show that the additive effects of coinfections, diversity and host energic status yield a bimodal distribution of lysogeny as a function of microbial densities. At high cell densities (above 106 cells ml-1 or g-1 ) and low diversity, coinfections by two or more phages are frequent and excess energy availability stimulates inefficient metabolism. Both mechanisms favour phage integration and characterize the Piggyback-the-Winner dynamic. At low densities (below 105 cells ml-1 or g-1 ), starvation represses lytic genes and extends the time window for lysogenic commitment, resulting in a higher frequency of coinfections that cause integration. This pattern follows the predictions of the refugium hypothesis. At intermediary densities (between 105 and 106 cells ml-1 or g-1 ), encounter rates and efficient energy metabolism favour lysis. This may involve Kill-the-Winner lytic dynamics and induction. Based on these three regimes, we propose a framework wherein phage integration occurs more frequently at both ends of the host density gradient, with distinct underlying molecular mechanisms (coinfections and host metabolism) dominating at each extreme.

Genomic and ecological attributes of marine bacteriophages encoding bacterial virulence genes.

BACKGROUND: Bacteriophages encode genes that modify bacterial functions during infection. The acquisition of phage-encoded virulence genes is a major mechanism for the rise of bacterial pathogens. In coral reefs, high bacterial density and lysogeny has been proposed to exacerbate reef decline through the transfer of phage-encoded virulence genes. However, the functions and distribution of these genes in phage virions on the reef remain unknown. RESULTS: Here, over 28,000 assembled viral genomes from the free viral community in Atlantic and Pacific Ocean coral reefs were queried against a curated database of virulence genes. The diversity of virulence genes encoded in the viral genomes was tested for relationships with host taxonomy and bacterial density in the environment. These analyses showed that bacterial density predicted the profile of virulence genes encoded by phages. The Shannon diversity of virulence-encoding phages was negatively related with bacterial density, leading to dominance of fewer genes at high bacterial abundances. A statistical learning analysis showed that reefs with high microbial density were enriched in viruses encoding genes enabling bacterial recognition and invasion of metazoan epithelium. Over 60% of phages could not have their hosts identified due to limitations of host prediction tools; for those which hosts were identified, host taxonomy was not an indicator of the presence of virulence genes. CONCLUSIONS: This study described bacterial virulence factors encoded in the genomes of bacteriophages at the community level. The results showed that the increase in microbial densities that occurs during coral reef degradation is associated with a change in the genomic repertoire of bacteriophages, specifically in the diversity and distribution of bacterial virulence genes. This suggests that phages are implicated in the rise of pathogens in disturbed marine ecosystems.

Correlation among DNA Linker Length, Linker Histone Concentration, and Histone Tails in Chromatin.

Eukaryotic cells condense their genetic material in the nucleus in the form of chromatin, a macromolecular complex made of DNA and multiple proteins. The structure of chromatin is intimately connected to the regulation of all eukaryotic organisms, from amoebas to humans, but its organization remains largely unknown. The nucleosome repeat length (NRL) and the concentration of linker histones (ρLH) are two structural parameters that vary among cell types and cell cycles; the NRL is the number of DNA basepairs wound around each nucleosome core plus the number of basepairs linking successive nucleosomes. Recent studies have found a linear empirical relationship between the variation of these two properties for different cells, but its underlying mechanism remains elusive. Here we apply our established mesoscale chromatin model to explore the mechanisms responsible for this relationship, by investigating chromatin fibers as a function of NRL and ρLH combinations. We find that a threshold of linker histone concentration triggers the compaction of chromatin into well-formed 30-nm fibers; this critical value increases linearly with NRL, except for long NRLs, where the fibers remain disorganized. Remarkably, the interaction patterns between core histone tails and chromatin elements are highly sensitive to the NRL and ρLH combination, suggesting a molecular mechanism that could have a key role in regulating the structural state of the fibers in the cell. An estimate of the minimized work and volume associated with storage of chromatin fibers in the nucleus further suggests factors that could spontaneously regulate the NRL as a function of linker histone concentration. Both the tail interaction map and DNA packing considerations support the empirical NRL/ρLH relationship and offer a framework to interpret experiments for different chromatin conditions in the cell.

Dynamic condensation of linker histone C-terminal domain regulates chromatin structure.

The basic and intrinsically disordered C-terminal domain (CTD) of the linker histone (LH) is essential for chromatin compaction. However, its conformation upon nucleosome binding and its impact on chromatin organization remain unknown. Our mesoscale chromatin model with a flexible LH CTD captures a dynamic, salt-dependent condensation mechanism driven by charge neutralization between the LH and linker DNA. Namely, at low salt concentration, CTD condenses, but LH only interacts with the nucleosome and one linker DNA, resulting in a semi-open nucleosome configuration; at higher salt, LH interacts with the nucleosome and two linker DNAs, promoting stem formation and chromatin compaction. CTD charge reduction unfolds the domain and decondenses chromatin, a mechanism in consonance with reduced counterion screening in vitro and phosphorylated LH in vivo. Divalent ions counteract this decondensation effect by maintaining nucleosome stems and expelling the CTDs to the fiber exterior. Additionally, we explain that the CTD folding depends on the chromatin fiber size, and we show that the asymmetric structure of the LH globular head is responsible for the uneven interaction observed between the LH and the linker DNAs. All these mechanisms may impact epigenetic regulation and higher levels of chromatin folding.

Theoretical studies on assembly, physical stability and dynamics of viruses.

All matter has to obey the general laws of physics and living matter is not an exception. Viruses have not only learnt how to cope with them, but have managed to use them for their own survival. In this chapter we will review some of the exciting physics behind viruses and discuss simple physical models that can shed some light on different aspects of the viral life cycle and viral properties. In particular, we will focus on how the structure and shape of the capsid, its assembly and stability, and the entry and exit of viral particles and their genomes can be understood using fundamental physics theories.

Optimal architectures of elongated viruses.

Many viruses protect their genetic material by a closed elongated protein shell. Unlike spherical viruses, the structure of these prolates is not yet well understood, and only a few of them have been fully characterized. We present the results of a simple phenomenological model, which describes the remarkable structures of prolate or bacilliform viral shells. Surprisingly, we find that the special well-defined geometry of these elongated viruses arises just as a consequence of free-energy minimization of a generic interaction between the structural units of the capsid. Hemispherical T-number caps centered along the 5-, 3-, and 2-fold axes with hexagonally ordered cylindrical bodies are found to be local energy minima, thus justifying their occurrence as optimal viral structures. Moreover, closed elongated viruses show a sequence of magic numbers for the end-caps, leading to strict selection rules for the length and structure of the body as well as for the number of capsomers and proteins of the capsid. The model reproduces the architecture of spherical and bacilliform viruses, both in vivo and in vitro, and constitutes an important step towards understanding viral assembly and its potential control for biological and nanotechnological applications.

A novel stabilization mechanism accommodating genome length variation in evolutionarily related viral capsids.

Tailed bacteriophages are one of the most numerous and diverse group of viruses. They store their genome at quasi-crystalline densities in capsids built from multiple copies of proteins adopting the HK97-fold. The high density of the genome exerts an internal pressure, requiring a maturation process that reinforces their capsids. However, it is unclear how capsid stabilization strategies have adapted to accommodate the evolution of larger genomes in this virus group. Here we characterized a novel capsid reinforcement mechanism in two evolutionary-related actinobacteriophages that modifies the length of a stabilization protein to accommodate a larger genome while maintaining the same capsid size. We used cryo-EM to reveal that capsids contained split hexamers of HK97-fold proteins with a stabilization protein in the chasm. The observation of split hexamers in mature capsids was unprecedented, so we rationalized this result mathematically, discovering that icosahedral capsids can be formed by all split or skewed hexamers as long as their T-number is not a multiple of three. Our results suggest that analogous stabilization mechanisms can be present in other icosahedral capsids, and they provide a strategy for engineering capsids accommodating larger DNA cargoes as gene delivery systems.

The Promise and Pitfalls of Prophages.

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2.

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome. Impact statement: Bacteriophages play a crucial role in shaping microbial communities within the human gut. Among the most dominant bacteriophages in the human gut microbiome are Crassvirales phages, which infect Bacteroides. Despite being widely distributed, only a few Crassvirales genomes have been isolated, leading to a limited understanding of their biology, ecology, and evolution. This study isolated and characterized three novel Crassvirales genomes belonging to two different families, and three genera, but infecting one bacterial host, Bacteroides cellulosilyticus WH2. Notably, the observation confirmed the phages are not co-evolving with their bacterial hosts, rather have a shared ability to exploit similar features in their bacterial host. Additionally, the identification of a critical viral protein undergoing purifying selection and interacting with the bacterial receptors opens doors to targeted therapies against bacterial infections. Given Bacteroides role in polysaccharide degradation in the human gut, our findings advance our understanding of the phage-host interactions and could have important implications for the development of phage-based therapies. These discoveries may hold implications for improving gut health and metabolism to support overall well-being. Data summary: The genomes used in this research are available on Sequence Read Archive (SRA) within the project, PRJNA737576. Bacteroides cellulosilyticus WH2, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. ' frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11 are all available on GenBank with accessions NZ_CP072251.1 ( B. cellulosilyticus WH2), QQ198717 (Bc01), QQ198718 (Bc03), and QQ198719 (Bc11), and we are working on making the strains available through ATCC. The 3D protein structures for the three Crassvirales genomes are available to download at doi.org/10.25451/flinders.21946034.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2.

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

Relevance of capsid structure in the buckling and maturation of spherical viruses.

The shape and mechanical properties of viral capsids play an important role in several biological processes during the virus life cycle. In particular, to become infective, many viruses require a maturation stage where the capsid undergoes a buckling transition, from an initial spherical procapsid into a final icosahedral faceted shell. Here we study, using a minimal physical model, how the capsid shape and the buckling transition depend on the triangulation number T and the icosahedral class P of the virus structure. We find that, for small shells, capsids with P = 1 are most likely to produce polyhedral shapes that minimize their energy and accumulated stress, whereas viruses with P = 3 prefer to remain spherical. For big capsids, all shells are more stable adopting an icosahedral shape, in agreement with continuum elastic theory. Moreover, spherical viruses show a buckling transition to polyhedral shells under expansion, in consonance with virus maturation. The resulting icosahedral shell is mechanically stiffer, tolerates larger expansions and withstands higher internal pressures before failing, which could explain why some dsDNA viruses, which rely on the pressurization of their genetic material to facilitate the infection, undergo a buckling transition. We emphasize that the results are general and could also be applied to non-biological systems.

Physics of shell assembly: line tension, hole implosion, and closure catastrophe.

The self-assembly of perfectly ordered closed shells is a challenging process involved in many biological and nanoscale systems. However, most of the aspects that determine their formation are still unknown. Here we investigate the growth of shells by simulating the assembly of spherical structures made of N identical subunits. Remarkably, we show that the formation and energetics of partially assembled shells are dominated by an effective line-tension that can be described in simple thermodynamic terms. In addition, we unveil two mechanisms that can prevent the correct formation of defect-free structures: "hole implosion," which leads to a premature closure of the shell; and "closure catastrophe," which causes a dramatic production of structural disorder during the later stages of the growth of big shells.

Predicting the capsid architecture of phages from metagenomic data.

Tailed phages are viruses that infect bacteria and are the most abundant biological entities on Earth. Their ecological, evolutionary, and biogeochemical roles in the planet stem from their genomic diversity. Known tailed phage genomes range from 10 to 735 kilobase pairs thanks to the size variability of the protective protein capsids that store them. However, the role of tailed phage capsids' diversity in ecosystems is unclear. A fundamental gap is the difficulty of associating genomic information with viral capsids in the environment. To address this problem, here, we introduce a computational approach to predict the capsid architecture (T-number) of tailed phages using the sequence of a single gene-the major capsid protein. This approach relies on an allometric model that relates the genome length and capsid architecture of tailed phages. This allometric model was applied to isolated phage genomes to generate a library that associated major capsid proteins and putative capsid architectures. This library was used to train machine learning methods, and the most computationally scalable model investigated (random forest) was applied to human gut metagenomes. Compared to isolated phages, the analysis of gut data reveals a large abundance of mid-sized (T = 7) capsids, as expected, followed by a relatively large frequency of jumbo-like tailed phage capsids (T ≥ 25) and small capsids (T = 4) that have been under-sampled. We discussed how to increase the method's accuracy and how to extend the approach to other viruses. The computational pipeline introduced here opens the doors to monitor the ongoing evolution and selection of viral capsids across ecosystems.

The International Virus Bioinformatics Meeting 2022.

The International Virus Bioinformatics Meeting 2022 took place online, on 23-25 March 2022, and has attracted about 380 participants from all over the world. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The participants created a highly interactive scientific environment even without physical face-to-face interactions. This meeting is a focal point to gain an insight into the state-of-the-art of the virus bioinformatics research landscape and to interact with researchers in the forefront as well as aspiring young scientists. The meeting featured eight invited and 18 contributed talks in eight sessions on three days, as well as 52 posters, which were presented during three virtual poster sessions. The main topics were: SARS-CoV-2, viral emergence and surveillance, virus-host interactions, viral sequence analysis, virus identification and annotation, phages, and viral diversity. This report summarizes the main research findings and highlights presented at the meeting.

Space-filling and benthic competition on coral reefs.

Reef-building corals are ecosystem engineers that compete with other benthic organisms for space and resources. Corals harvest energy through their surface by photosynthesis and heterotrophic feeding, and they divert part of this energy to defend their outer colony perimeter against competitors. Here, we hypothesized that corals with a larger space-filling surface and smaller perimeters increase energy gain while reducing the exposure to competitors. This predicted an association between these two geometric properties of corals and the competitive outcome against other benthic organisms. To test the prediction, fifty coral colonies from the Caribbean island of Curaçao were rendered using digital 3D and 2D reconstructions. The surface areas, perimeters, box-counting dimensions (as a proxy of surface and perimeter space-filling), and other geometric properties were extracted and analyzed with respect to the percentage of the perimeter losing or winning against competitors based on the coral tissue apparent growth or damage. The increase in surface space-filling dimension was the only significant single indicator of coral winning outcomes, but the combination of surface space-filling dimension with perimeter length increased the statistical prediction of coral competition outcomes. Corals with larger surface space-filling dimensions (Ds > 2) and smaller perimeters displayed more winning outcomes, confirming the initial hypothesis. We propose that the space-filling property of coral surfaces complemented with other proxies of coral competitiveness, such as life history traits, will provide a more accurate quantitative characterization of coral competition outcomes on coral reefs. This framework also applies to other organisms or ecological systems that rely on complex surfaces to obtain energy for competition.

The structure of elongated viral capsids.

There are many viruses whose genetic material is protected by a closed elongated protein shell. Unlike spherical viruses, the structure and construction principles of these elongated capsids are not fully known. In this article, we have developed a general geometrical model to describe the structure of prolate or bacilliform capsids. We show that only a limited set of tubular architectures can be built closed by hemispherical icosahedral caps. In particular, the length and number of proteins adopt a very special set of discrete values dictated by the axial symmetry (fivefold, threefold, or twofold) and the triangulation number of the caps. The results are supported by experimental observations and simulations of simplified physical models. This work brings about a general classification of elongated viruses that will help to predict their structure, and to design viral cages with tailored geometrical properties for biomedical and nanotechnological applications.

Quantification of Lysogeny Caused by Phage Coinfections in Microbial Communities from Biophysical Principles.

Temperate phages can associate with their bacterial host to form a lysogen, often modifying the phenotype of the host. Lysogens are dominant in the microbially dense environment of the mammalian gut. This observation contrasts with the long-standing hypothesis of lysogeny being favored at low microbial densities, such as in oligotrophic marine environments. Here, we hypothesized that phage coinfections-a well-understood molecular mechanism of lysogenization-increase at high microbial abundances. To test this hypothesis, we developed a biophysical model of coinfection for marine and gut microbiomes. The model stochastically sampled ranges of phage and bacterial concentrations, adsorption rates, lysogenic commitment times, and community diversity from each environment. In 90% of the sampled marine communities, less than 10% of the bacteria were predicted to be lysogenized via coinfection. In contrast, 25% of the sampled gut communities displayed more than 25% of lysogenization. The probability of lysogenization in the gut was a consequence of the higher densities and higher adsorption rates. These results suggest that, on average, coinfections can form two trillion lysogens in the human gut every day. In marine microbiomes, which were characterized by lower densities and phage adsorption rates, lysogeny via coinfection was still possible for communities with long lysogenic commitment times. Our study indicates that different physical factors causing coinfections can reconcile the traditional view of lysogeny at poor host growth (long commitment times) and the recent Piggyback-the-Winner framework proposing that lysogeny is favored in rich environments (high densities and adsorption rates).IMPORTANCE The association of temperate phages and bacterial hosts during lysogeny manipulates microbial dynamics from the oceans to the human gut. Lysogeny is well studied in laboratory models, but its environmental drivers remain unclear. Here, we quantified the probability of lysogenization caused by phage coinfections, a well-known trigger of lysogeny, in marine and gut microbial environments. Coinfections were quantified by developing a biophysical model that incorporated the traits of viral and bacterial communities. Lysogenization via coinfection was more frequent in highly productive environments like the gut, due to higher microbial densities and higher phage adsorption rates. At low cell densities, lysogenization occurred in bacteria with long duplication times. These results bridge the molecular understanding of lysogeny with the ecology of complex microbial communities.

The Missing Tailed Phages: Prediction of Small Capsid Candidates.

Tailed phages are the most abundant and diverse group of viruses on the planet. Yet, the smallest tailed phages display relatively complex capsids and large genomes compared to other viruses. The lack of tailed phages forming the common icosahedral capsid architectures T = 1 and T = 3 is puzzling. Here, we extracted geometrical features from high-resolution tailed phage capsid reconstructions and built a statistical model based on physical principles to predict the capsid diameter and genome length of the missing small-tailed phage capsids. We applied the model to 3348 isolated tailed phage genomes and 1496 gut metagenome-assembled tailed phage genomes. Four isolated tailed phages were predicted to form T = 3 icosahedral capsids, and twenty-one metagenome-assembled tailed phages were predicted to form T < 3 capsids. The smallest capsid predicted was a T = 4/3 ≈ 1.33 architecture. No tailed phages were predicted to form the smallest icosahedral architecture, T = 1. We discuss the feasibility of the missing T = 1 tailed phage capsids and the implications of isolating and characterizing small-tailed phages for viral evolution and phage therapy.