Widespread endogenization of giant viruses shapes genomes of green algae.

Endogenous viral elements (EVEs)-viruses that have integrated their genomes into those of their hosts-are prevalent in eukaryotes and have an important role in genome evolution1,2. The vast majority of EVEs that have been identified to date are small genomic regions comprising a few genes2, but recent evidence suggests that some large double-stranded DNA viruses may also endogenize into the genome of the host1. Nucleocytoplasmic large DNA viruses (NCLDVs) have recently become of great interest owing to their large genomes and complex evolutionary origins3-6, but it is not yet known whether they are a prominent component of eukaryotic EVEs. Here we report the widespread endogenization of NCLDVs in diverse green algae; these giant EVEs reached sizes greater than 1 million base pairs and contained as many as around 10% of the total open reading frames in some genomes, substantially increasing the scale of known viral genes in eukaryotic genomes. These endogenized elements often shared genes with host genomic loci and contained numerous spliceosomal introns and large duplications, suggesting tight assimilation into host genomes. NCLDVs contain large and mosaic genomes with genes derived from multiple sources, and their endogenization represents an underappreciated conduit of new genetic material into eukaryotic lineages that can substantially impact genome composition.

The genomic and phylogenetic analysis of Marseillevirus cajuinensis raises questions about the evolution of Marseilleviridae lineages and their taxonomical organization.

Marseilleviruses (MsV) are a group of viruses that compose the Marseilleviridae family within the Nucleocytoviricota phylum. They have been found in different samples, mainly in freshwater. MsV are classically organized into five phylogenetic lineages (A/B/C/D/E), but the current taxonomy does not fully represent all the diversity of the MsV lineages. Here, we describe a novel strain isolated from a Brazilian saltwater sample named Marseillevirus cajuinensis. Based on genomics and phylogenetic analyses, M. cajuinensis exhibits a 380,653-bp genome that encodes 515 open reading frames. Additionally, M. cajuinensis encodes a transfer RNA, a feature that is rarely described for Marseilleviridae. Phylogeny suggests that M. cajuinensis forms a divergent branch within the MsV lineage A. Furthermore, our analysis suggests that the common ancestor for the five classical lineages of MsV diversified into three major groups. The organization of MsV into three main groups is reinforced by a comprehensive analysis of clusters of orthologous groups, sequence identities, and evolutionary distances considering several MsV isolates. Taken together, our results highlight the importance of discovering new viruses to expand the knowledge about known viruses that belong to the same lineages or families. This work proposes a new perspective on the Marseilleviridae lineages organization that could be helpful to a future update in the taxonomy of the Marseilleviridae family. IMPORTANCE: Marseilleviridae is a family of viruses whose members were mostly isolated from freshwater samples. In this work, we describe the first Marseillevirus isolated from saltwater samples, which we called Marseillevirus cajuinensis. Most of M. cajuinensis genomic features are comparable to other Marseilleviridae members, such as its high number of unknown proteins. On the other hand, M. cajuinensis encodes a transfer RNA, which is a gene category involved in protein translation that is rarely described in this viral family. Additionally, our phylogenetic analyses suggested the existence of, at least, three major Marseilleviridae groups. These observations provide a new perspective on Marseilleviridae lineages organization, which will be valuable in future updates to the taxonomy of the family since the current official classification does not capture all the Marseilleviridae known diversity.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales.

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel's lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

Gene duplication as a major force driving the genome expansion in some giant viruses.

IMPORTANCE: Giant viruses are noteworthy not only due to their enormous particles but also because of their gigantic genomes. In this context, a fundamental question has persisted: how did these genomes evolve? Here we present the discovery of cedratvirus pambiensis, featuring the largest genome ever described for a cedratvirus. Our data suggest that the larger size of the genome can be attributed to an unprecedented number of duplicated genes. Further investigation of this phenomenon in other viruses has illuminated gene duplication as a key evolutionary mechanism driving genome expansion in diverse giant viruses. Although gene duplication has been described as a recurrent event in cellular organisms, our data highlights its potential as a pivotal event in the evolution of gigantic viral genomes.

A phylogenomic framework for charting the diversity and evolution of giant viruses.

Large DNA viruses of the phylum Nucleocytoviricota have recently emerged as important members of ecosystems around the globe that challenge traditional views of viral complexity. Numerous members of this phylum that cannot be classified within established families have recently been reported, and there is presently a strong need for a robust phylogenomic and taxonomic framework for these viruses. Here, we report a comprehensive phylogenomic analysis of the Nucleocytoviricota, present a set of giant virus orthologous groups (GVOGs) together with a benchmarked reference phylogeny, and delineate a hierarchical taxonomy within this phylum. We show that the majority of Nucleocytoviricota diversity can be partitioned into 6 orders, 32 families, and 344 genera, substantially expanding the number of currently recognized taxonomic ranks for these viruses. We integrate our results within a taxonomy that has been adopted for all viruses to establish a unifying framework for the study of Nucleocytoviricota diversity, evolution, and environmental distribution.

Taxonomic update for giant viruses in the order Imitervirales (phylum Nucleocytoviricota).

Large DNA viruses in the phylum Nucleocytoviricota, sometimes referred to as "giant viruses" owing to their large genomes and virions, have been the subject of burgeoning interest over the last decade. Here, we describe recently adopted taxonomic updates for giant viruses within the order Imitervirales. The families Allomimiviridae, Mesomimiviridae, and Schizomimiviridae have been created to accommodate the increasing diversity of mimivirus relatives that have sometimes been referred to in the literature as "extended Mimiviridae". In addition, the subfamilies Aliimimivirinae, Megamimivirinae, and Klosneuvirinae have been established to refer to subgroups of the Mimiviridae. Binomial names have also been adopted for all recognized species in the order. For example, Acanthamoeba polyphaga mimivirus is now classified in the species Mimivirus bradfordmassiliense.

Phylogenetic Signal, Congruence, and Uncertainty across Bacteria and Archaea.

Reconstruction of the Tree of Life is a central goal in biology. Although numerous novel phyla of bacteria and archaea have recently been discovered, inconsistent phylogenetic relationships are routinely reported, and many inter-phylum and inter-domain evolutionary relationships remain unclear. Here, we benchmark different marker genes often used in constructing multidomain phylogenetic trees of bacteria and archaea and present a set of marker genes that perform best for multidomain trees constructed from concatenated alignments. We use recently-developed Tree Certainty metrics to assess the confidence of our results and to obviate the complications of traditional bootstrap-based metrics. Given the vastly disparate number of genomes available for different phyla of bacteria and archaea, we also assessed the impact of taxon sampling on multidomain tree construction. Our results demonstrate that biases between the representation of different taxonomic groups can dramatically impact the topology of resulting trees. Inspection of our highest-quality tree supports the division of most bacteria into Terrabacteria and Gracilicutes, with Thermatogota and Synergistota branching earlier from these superphyla. This tree also supports the inclusion of the Patescibacteria within the Terrabacteria as a sister group to the Chloroflexota instead of as a basal-branching lineage. For the Archaea, our tree supports three monophyletic lineages (DPANN, Euryarchaeota, and TACK/Asgard), although we note the basal placement of the DPANN may still represent an artifact caused by biased sequence composition. Our findings provide a robust and standardized framework for multidomain phylogenetic reconstruction that can be used to evaluate inter-phylum relationships and assess uncertainty in conflicting topologies of the Tree of Life.

Historical forest disturbance mediates soil microbial community responses to drought.

Despite the abundance of studies demonstrating the effects of drought on soil microbial communities, the role of land use legacies in mediating these drought effects is unclear. To assess historical land use influences on microbial drought responses, we conducted a drought-rewetting experiment in soils from two adjacent and currently forested watersheds with distinct land use histories: an undisturbed 'reference' site and a 'disturbed' site that was clear-cut and converted to agriculture ~60 years prior. We incubated intact soil cores at either constant moisture or under a drought-rewet treatment and characterized bacterial and fungal communities using amplicon sequencing throughout the experiment. Bacterial alpha diversity decreased following drought-rewetting while fungal diversity increased. Bacterial beta diversity also changed markedly following drought-rewetting, especially in historically disturbed soils, while fungal beta diversity exhibited little response. Additionally, bacterial beta diversity in disturbed soils recovered less from drought-rewetting compared with reference soils. Disturbed soil communities also exhibited notable reductions in nitrifying taxa, increases in putative r-selected bacteria, and reductions in network connectivity following drought-rewetting. Overall, our study reveals historical land use to be important in mediating responses of soil bacterial communities to drought, which will influence the ecosystem-scale trajectories of these environments under ongoing and future climate change.

Diel cycling and long-term persistence of viruses in the ocean's euphotic zone.

Viruses are fundamental components of marine microbial communities that significantly influence oceanic productivity, biogeochemistry, and ecosystem processes. Despite their importance, the temporal activities and dynamics of viral assemblages in natural settings remain largely unexplored. Here we report the transcriptional activities and variability of dominant dsDNA viruses in the open ocean's euphotic zone over daily and seasonal timescales. While dsDNA viruses exhibited some fluctuation in abundance in both cellular and viral size fractions, the viral assemblage was remarkably stable, with the most abundant viral types persisting over many days. More extended time series indicated that long-term persistence (>1 y) was the rule for most dsDNA viruses observed, suggesting that both core viral genomes as well as viral community structure were conserved over interannual periods. Viral gene transcription in host cell assemblages revealed diel cycling among many different viral types. Most notably, an afternoon peak in cyanophage transcriptional activity coincided with a peak in Prochlorococcus DNA replication, indicating coordinated diurnal coupling of virus and host reproduction. In aggregate, our analyses suggested a tightly synchronized diel coupling of viral and cellular replication cycles in both photoautotrophic and heterotrophic bacterial hosts. A surprising consequence of these findings is that diel cycles in the ocean's photic zone appear to be universal organizing principles that shape ecosystem dynamics, ecological interactions, and biogeochemical cycling of both cellular and acellular community components.

Microbial community transcriptional networks are conserved in three domains at ocean basin scales.

Planktonic microbial communities in the ocean are typically dominated by several cosmopolitan clades of Bacteria, Archaea, and Eukarya characterized by their ribosomal RNA gene phylogenies and genomic features. Although the environments these communities inhabit range from coastal to open ocean waters, how the biological dynamics vary between such disparate habitats is not well known. To gain insight into the differential activities of microbial populations inhabiting different oceanic provinces we compared the daily metatranscriptome profiles of related microbial populations inhabiting surface waters of both a coastal California upwelling region (CC) as well as the oligotrophic North Pacific Subtropical Gyre (NPSG). Transcriptional networks revealed that the dominant photoautotrophic microbes in each environment (Ostreococcus in CC, Prochlorococcus in NPSG) were central determinants of overall community transcriptome dynamics. Furthermore, heterotrophic bacterial clades common to both ecosystems (SAR11, SAR116, SAR86, SAR406, and Roseobacter) displayed conserved, genome-wide inter- and intrataxon transcriptional patterns and diel cycles. Populations of SAR11 and SAR86 clades in particular exhibited tightly coordinated transcriptional patterns in both coastal and pelagic ecosystems, suggesting that specific biological interactions between these groups are widespread in nature. Our results identify common diurnally oscillating behaviors among diverse planktonic microbial species regardless of habitat, suggesting that highly conserved temporally phased biotic interactions are ubiquitous among planktonic microbial communities worldwide.

Single-cell RNA-seq of the rare virosphere reveals the native hosts of giant viruses in the marine environment.

Giant viruses (phylum Nucleocytoviricota) are globally distributed in aquatic ecosystems. They play fundamental roles as evolutionary drivers of eukaryotic plankton and regulators of global biogeochemical cycles. However, we lack knowledge about their native hosts, hindering our understanding of their life cycle and ecological importance. In the present study, we applied a single-cell RNA sequencing (scRNA-seq) approach to samples collected during an induced algal bloom, which enabled pairing active giant viruses with their native protist hosts. We detected hundreds of single cells from multiple host lineages infected by diverse giant viruses. These host cells included members of the algal groups Chrysophycae and Prymnesiophycae, as well as heterotrophic flagellates in the class Katablepharidaceae. Katablepharids were infected with a rare Imitervirales-07 giant virus lineage expressing a large repertoire of cell-fate regulation genes. Analysis of the temporal dynamics of these host-virus interactions revealed an important role for the Imitervirales-07 in controlling the population size of the host Katablepharid population. Our results demonstrate that scRNA-seq can be used to identify previously undescribed host-virus interactions and study their ecological importance and impact.

Mountain pine beetles colonizing historical and naive host trees are associated with a bacterial community highly enriched in genes contributing to terpene metabolism.

The mountain pine beetle, Dendroctonus ponderosae, is a subcortical herbivore native to western North America that can kill healthy conifers by overcoming host tree defenses, which consist largely of high terpene concentrations. The mechanisms by which these beetles contend with toxic compounds are not well understood. Here, we explore a component of the hypothesis that beetle-associated bacterial symbionts contribute to the ability of D. ponderosae to overcome tree defenses by assisting with terpene detoxification. Such symbionts may facilitate host tree transitions during range expansions currently being driven by climate change. For example, this insect has recently breached the historical geophysical barrier of the Canadian Rocky Mountains, providing access to näive tree hosts and unprecedented connectivity to eastern forests. We use culture-independent techniques to describe the bacterial community associated with D. ponderosae beetles and their galleries from their historical host, Pinus contorta, and their more recent host, hybrid P. contorta-Pinus banksiana. We show that these communities are enriched with genes involved in terpene degradation compared with other plant biomass-processing microbial communities. These pine beetle microbial communities are dominated by members of the genera Pseudomonas, Rahnella, Serratia, and Burkholderia, and the majority of genes involved in terpene degradation belong to these genera. Our work provides the first metagenome of bacterial communities associated with a bark beetle and is consistent with a potential microbial contribution to detoxification of tree defenses needed to survive the subcortical environment.

Comparison of 26 sphingomonad genomes reveals diverse environmental adaptations and biodegradative capabilities.

Sphingomonads comprise a physiologically versatile group within the Alphaproteobacteria that includes strains of interest for biotechnology, human health, and environmental nutrient cycling. In this study, we compared 26 sphingomonad genome sequences to gain insight into their ecology, metabolic versatility, and environmental adaptations. Our multilocus phylogenetic and average amino acid identity (AAI) analyses confirm that Sphingomonas, Sphingobium, Sphingopyxis, and Novosphingobium are well-resolved monophyletic groups with the exception of Sphingomonas sp. strain SKA58, which we propose belongs to the genus Sphingobium. Our pan-genomic analysis of sphingomonads reveals numerous species-specific open reading frames (ORFs) but few signatures of genus-specific cores. The organization and coding potential of the sphingomonad genomes appear to be highly variable, and plasmid-mediated gene transfer and chromosome-plasmid recombination, together with prophage- and transposon-mediated rearrangements, appear to play prominent roles in the genome evolution of this group. We find that many of the sphingomonad genomes encode numerous oxygenases and glycoside hydrolases, which are likely responsible for their ability to degrade various recalcitrant aromatic compounds and polysaccharides, respectively. Many of these enzymes are encoded on megaplasmids, suggesting that they may be readily transferred between species. We also identified enzymes putatively used for the catabolism of sulfonate and nitroaromatic compounds in many of the genomes, suggesting that plant-based compounds or chemical contaminants may be sources of nitrogen and sulfur. Many of these sphingomonads appear to be adapted to oligotrophic environments, but several contain genomic features indicative of host associations. Our work provides a basis for understanding the ecological strategies employed by sphingomonads and their role in environmental nutrient cycling.

Leucoagaricus gongylophorus produces diverse enzymes for the degradation of recalcitrant plant polymers in leaf-cutter ant fungus gardens.

Plants represent a large reservoir of organic carbon comprised primarily of recalcitrant polymers that most metazoans are unable to deconstruct. Many herbivores gain access to nutrients in this material indirectly by associating with microbial symbionts, and leaf-cutter ants are a paradigmatic example. These ants use fresh foliar biomass as manure to cultivate gardens composed primarily of Leucoagaricus gongylophorus, a basidiomycetous fungus that produces specialized hyphal swellings that serve as a food source for the host ant colony. Although leaf-cutter ants are conspicuous herbivores that contribute substantially to carbon turnover in Neotropical ecosystems, the process through which plant biomass is degraded in their fungus gardens is not well understood. Here we present the first draft genome of L. gongylophorus, and, using genomic and metaproteomic tools, we investigate its role in lignocellulose degradation in the gardens of both Atta cephalotes and Acromyrmex echinatior leaf-cutter ants. We show that L. gongylophorus produces a diversity of lignocellulases in ant gardens and is likely the primary driver of plant biomass degradation in these ecosystems. We also show that this fungus produces distinct sets of lignocellulases throughout the different stages of biomass degradation, including numerous cellulases and laccases that likely play an important role in lignocellulose degradation. Our study provides a detailed analysis of plant biomass degradation in leaf-cutter ant fungus gardens and insight into the enzymes underlying the symbiosis between these dominant herbivores and their obligate fungal cultivar.

An insect herbivore microbiome with high plant biomass-degrading capacity.

Herbivores can gain indirect access to recalcitrant carbon present in plant cell walls through symbiotic associations with lignocellulolytic microbes. A paradigmatic example is the leaf-cutter ant (Tribe: Attini), which uses fresh leaves to cultivate a fungus for food in specialized gardens. Using a combination of sugar composition analyses, metagenomics, and whole-genome sequencing, we reveal that the fungus garden microbiome of leaf-cutter ants is composed of a diverse community of bacteria with high plant biomass-degrading capacity. Comparison of this microbiome's predicted carbohydrate-degrading enzyme profile with other metagenomes shows closest similarity to the bovine rumen, indicating evolutionary convergence of plant biomass degrading potential between two important herbivorous animals. Genomic and physiological characterization of two dominant bacteria in the fungus garden microbiome provides evidence of their capacity to degrade cellulose. Given the recent interest in cellulosic biofuels, understanding how large-scale and rapid plant biomass degradation occurs in a highly evolved insect herbivore is of particular relevance for bioenergy.

Microbiology: The curious case of the mysterious mirusvirus.

Long-read sequencing of a marine stramenopile genome yields a trove of insights into protist genomics and solves a 50-year-old viral mystery.

Automated classification of giant virus genomes using a random forest model built on trademark protein families.

Viruses of the phylum Nucleocytoviricota, often referred to as "giant viruses," are prevalent in various environments around the globe and play significant roles in shaping eukaryotic diversity and activities in global ecosystems. Given the extensive phylogenetic diversity within this viral group and the highly complex composition of their genomes, taxonomic classification of giant viruses, particularly incomplete metagenome-assembled genomes (MAGs) can present a considerable challenge. Here we developed TIGTOG (Taxonomic Information of Giant viruses using Trademark Orthologous Groups), a machine learning-based approach to predict the taxonomic classification of novel giant virus MAGs based on profiles of protein family content. We applied a random forest algorithm to a training set of 1,531 quality-checked, phylogenetically diverse Nucleocytoviricota genomes using pre-selected sets of giant virus orthologous groups (GVOGs). The classification models were predictive of viral taxonomic assignments with a cross-validation accuracy of 99.6% to the order level and 97.3% to the family level. We found that no individual GVOGs or genome features significantly influenced the algorithm's performance or the models' predictions, indicating that classification predictions were based on a comprehensive genomic signature, which reduced the necessity of a fixed set of marker genes for taxonomic assigning purposes. Our classification models were validated with an independent test set of 823 giant virus genomes with varied genomic completeness and taxonomy and demonstrated an accuracy of 98.6% and 95.9% to the order and family level, respectively. Our results indicate that protein family profiles can be used to accurately classify large DNA viruses at different taxonomic levels and provide a fast and accurate method for the classification of giant viruses. This approach could easily be adapted to other viral groups.

Assessing the biogeography of marine giant viruses in four oceanic transects.

Viruses of the phylum Nucleocytoviricota are ubiquitous in ocean waters and play important roles in shaping the dynamics of marine ecosystems. In this study, we leveraged the bioGEOTRACES metagenomic dataset collected across the Atlantic and Pacific Oceans to investigate the biogeography of these viruses in marine environments. We identified 330 viral genomes, including 212 in the order Imitervirales and 54 in the order Algavirales . We found that most viruses appeared to be prevalent in shallow waters (<150 meters), and that viruses of the Mesomimiviridae ( Imitervirales ) and Prasinoviridae ( Algavirales ) are by far the most abundant and diverse groups in our survey. Five mesomimiviruses and one prasinovirus are particularly widespread in oligotrophic waters; annotation of these genomes revealed common stress response systems, photosynthesis-associated genes, and oxidative stress modulation that may be key to their broad distribution in the pelagic ocean. We identified a latitudinal pattern in viral diversity in one cruise that traversed the North and South Atlantic Ocean, with viral diversity peaking at high latitudes of the northern hemisphere. Community analyses revealed three distinct Nucleocytoviricota communities across latitudes, categorized by latitudinal distance towards the equator. Our results contribute to the understanding of the biogeography of these viruses in marine systems.