Genome Streamlining Research Articles

Genome streamlining: effect of mutation rate and population size on genome size reduction.

Genome streamlining, i.e. genome size reduction, is observed in bacteria with very different life traits, including endosymbiotic bacteria and several marine bacteria, raising the question of its evolutionary origin. None of the hypotheses proposed in the literature is firmly established, mainly due to the many confounding factors related to the diverse habitats of species with streamlined genomes. Computational models may help overcome these difficulties and rigorously test hypotheses. In this work, we used Aevol, a platform designed to study the evolution of genome architecture, to test two main hypotheses: that an increase in population size (N) or mutation rate (μ) could cause genome reduction. In our experiments, both conditions lead to streamlining but have very different resulting genome structures. Under increased population sizes, genomes lose a significant fraction of non-coding sequences but maintain their coding size, resulting in densely packed genomes (akin to streamlined marine bacteria genomes). By contrast, under an increased mutation rate, genomes lose both coding and non-coding sequences (akin to endosymbiotic bacteria genomes). Hence, both factors lead to an overall reduction in genome size, but the coding density of the genome appears to be determined by N × μ. Thus, a broad range of genome size and density can be achieved by different combinations of N and μ. Our results suggest that genome size and coding density are determined by the interplay between selection for phenotypic adaptation and selection for robustness.

Genome streamlining of Pseudomonas putida B6-2 for bioremediation.

Microbial transformation is a favored approach for environmental remediation. However, the effectiveness of microbial remediation has been limited by the lack of chassis cells with satisfactory contaminant degradation performance. Pseudomonas putida B6-2, with a wide substrate spectrum and high solvent tolerance, is a chassis strain with great potential for application in environmental remediation. Here, guided by bioinformatic analyses and genome-scale metabolic model (GEM) predictions, we successfully optimized P. putida B6-2 by rationally reducing its nonessential genetic components and generating a more robust genome-streamlined strain, P. putida BGR4. Several improvements were observed compared with the original P. putida B6-2 strain, including a 1.4 × 105-fold increase in electroporation efficiency, an 8.3-fold increase in conjugation efficiency, improved glycerol utilization capability, and increased phenol utilization after heterologous expression of the phenol monooxygenase encoded by dmpKLMNOP. Additionally, P. putida BGR4 exhibited enhanced tolerance to several stressors, including starvation, oxidative stress, and DNA damage. Transcriptomic analysis revealed that genome streamlining led to the upregulation of genes involved in the "carbon metabolism" and "tricarboxylic acid cycle" pathways in P. putida BGR4, which likely contributed to the superior phenotype of P. putida BGR4 in terms of carbon source utilization and contaminant degradation capabilities. Furthermore, the absence of four prophages was identified as a potential cause of the enhanced stress resistance observed in P. putida BGR4. Overall, we developed a combined genome-streamlining strategy involving bioinformatic analyses and GEM predictions and generated a more robust chassis strain, P. putida BGR4, which expands the repertoire of chassis cells for environmental remediation.IMPORTANCEDespite the development of many chassis cells, there is still a lack of robust chassis cells with satisfactory contaminant degradation performance. Targeted genome streamlining is an effective way to provide powerful chassis cells. However, genome streamlining does not always lead to the improved phenotypes of genome-streamlined chassis cells. In this research, a novel procedure that combined bioinformatic analyses and GEM predictions was proposed to guide genome streamlining and predict the effects of genome streamlining. This genome streamlining procedure was successfully applied to Pseudomonas putida B6-2, which was a chassis cell with great potential for application in environmental remediation and resulted in the generation of a more robust chassis cell, P. putida BGR4, thereby providing a superior chassis cell for efficient and sustainable environmental remediation and a valuable framework for guiding the genome streamlining of strains for other applications.

Ubiquitous genome streamlined Acidobacteriota in freshwater environments.

Acidobacteriota are abundant in soil, peatlands, and sediments, but their ecology in freshwater environments remains understudied. UBA12189, an Acidobacteriota genus, is an uncultivated, genome-streamlined lineage with a small genome size found in aquatic environments where detailed genomic analyses are lacking. Here, we analyzed 66 MAGs of UBA12189 (including one complete genome) from freshwater lakes and rivers in Europe, North America, and Asia. UBA12189 has small genome sizes (<1.4 Mbp), low GC content, and a highly diverse pangenome. In freshwater lakes, this bacterial lineage is abundant from the surface waters (epilimnion) down to a 300-m depth (hypolimnion). UBA12189 appears to be free-living from CARD-FISH analysis. When compared to other genome-streamlined bacteria such as Nanopelagicales and Methylopumilus, genome reduction has caused UBA12189 to have a more limited metabolic repertoire in carbon, sulfur, and nitrogen metabolisms, limited numbers of membrane transporters, as well as a higher degree of auxotrophy for various amino acids, vitamins, and reduced sulfur. Despite having reduced genomes, UBA12189 encodes proteorhodopsin, complete biosynthesis pathways for heme and vitamin K2, cbb3-type cytochrome c oxidases, and heme-requiring enzymes. These genes may give a selective advantage during the genome streamlining process. We propose the new genus Acidiparvus, with two new species named "A. lacustris" and "A. fluvialis". Acidiparvus is the first described genome-streamlined lineage under the phylum Acidobacteriota, which is a free-living, slow-growing scavenger in freshwater environments.

Diagnosis and mitigation of the systemic impact of genome reduction in Escherichia coli DGF-298.

Genomic streamlining can be employed in model organisms to reduce complexity and enhance strain predictability. One of the most striking examples is the bacterial strain Escherichia coli DGF-298, notable for having over one-third of its genome deleted. However, such extensive genome modifications raise the question of how similar this simplified cell remains when compared with its parent, and what are the possible unintended consequences of this simplification. In this study, we used metabolic modeling along with iModulon-based transcriptomic analysis in different growth conditions to assess the impact of genome reduction on metabolism and gene regulation. We observed little impact of genomic reduction on the regulatory network of E. coli DGF-298 and identified a potential metabolic bottleneck leading to the constitutive activity of the SoxS iModulon. We then leveraged the model's predictions to successfully restore SoxS activity to the basal level.

Evolution recovers the fitness of Acinetobacter baylyi strains with large deletions through mutations in deletion-specific targets and global post-transcriptional regulators.

Organelles and endosymbionts have naturally evolved dramatically reduced genome sizes compared to their free-living ancestors. Synthetic biologists have purposefully engineered streamlined microbial genomes to create more efficient cellular chassis and define the minimal components of cellular life. During natural or engineered genome streamlining, deletion of many non-essential genes in combination often reduces bacterial fitness for idiosyncratic or unknown reasons. We investigated how and to what extent laboratory evolution could overcome these defects in six variants of the transposon-free Acinetobacter baylyi strain ADP1-ISx that each had a deletion of a different 22- to 42-kilobase region and two strains with larger deletions of 70 and 293 kilobases. We evolved replicate populations of ADP1-ISx and each deletion strain for ~300 generations in a chemically defined minimal medium or a complex medium and sequenced the genomes of endpoint clonal isolates. Fitness increased in all cases that were examined except for two ancestors that each failed to improve in one of the two environments. Mutations affecting nine protein-coding genes and two small RNAs were significantly associated with one of the two environments or with certain deletion ancestors. The global post-transcriptional regulators rnd (ribonuclease D), csrA (RNA-binding carbon storage regulator), and hfq (RNA-binding protein and chaperone) were frequently mutated across all strains, though the incidence and effects of these mutations on gene function and bacterial fitness varied with the ancestral deletion and evolution environment. Mutations in this regulatory network likely compensate for how an earlier deletion of a transposon in the ADP1-ISx ancestor of all the deletion strains restored csrA function. More generally, our results demonstrate that fitness lost during genome streamlining can usually be regained rapidly through laboratory evolution and that recovery tends to occur through a combination of deletion-specific compensation and global regulatory adjustments.

New insights into decoding the lifestyle of endophytic Fusarium lateritium Fl617 via comparing genomes

Fungal-plant interactions have persisted for 460 million years, and almost all terrestrial plants on Earth have endophytic fungi. However, the mechanism of symbiosis between endophytic fungi and host plants has been inconclusive. In this dissertation, we used a strain of endophytic Fusarium lateritium (Fl617), which was found in the previous stage to promote disease resistance in tomato, and selected the pathogenic Fusarium oxysporum Fo4287 and endophytic Fusarium oxysporum Fo47, which are in the same host and the closest relatives of Fl617, to carry out a comparative genomics analysis of the three systems and to provide a new perspective for the elucidation of the special lifestyle of the fungal endophytes. We found that endophytic F. lateritium has a smaller genome, fewer clusters and genes associated with pathogenicity, and fewer plant cell wall degrading enzymes (PCWDEs). There were also relatively fewer secondary metabolisms and typical Fusarium spp. toxins, and a lack of the key Fusarium spp. pathogenicity factor, secreted in xylem (SIX), but the endophytic fungi may be more sophisticated in their regulation of the colonization process. It is hypothesized that the endophytic fungi may have maintained their symbiosis with plants due to the relatively homogeneous microenvironment in plants for a long period of time, considering only plant interactions and discarding the relevant pathogenicity factors, and that their endophytic evolutionary tendency may tend to be genome streamlining and to enhance the fineness of the regulation of plant interactions, thus maintaining their symbiotic status with plants.

Eco-evolutionary strategies for relieving carbon limitation under salt stress differ across microbial clades

With the continuous expansion of saline soils under climate change, understanding the eco-evolutionary tradeoff between the microbial mitigation of carbon limitation and the maintenance of functional traits in saline soils represents a significant knowledge gap in predicting future soil health and ecological function. Through shotgun metagenomic sequencing of coastal soils along a salinity gradient, we show contrasting eco-evolutionary directions of soil bacteria and archaea that manifest in changes to genome size and the functional potential of the soil microbiome. In salt environments with high carbon requirements, bacteria exhibit reduced genome sizes associated with a depletion of metabolic genes, while archaea display larger genomes and enrichment of salt-resistance, metabolic, and carbon-acquisition genes. This suggests that bacteria conserve energy through genome streamlining when facing salt stress, while archaea invest in carbon-acquisition pathways to broaden their resource usage. These findings suggest divergent directions in eco-evolutionary adaptations to soil saline stress amongst microbial clades and serve as a foundation for understanding the response of soil microbiomes to escalating climate change.

Assembly and comparative analyses of the Geosiphon pyriformis metagenome.

Geosiphon pyriformis, a representative of the fungal sub-phylum Glomeromycotina, is unique in its endosymbiosis with cyanobacteria within a fungal cell. This symbiotic relationship occurs in bladders containing nuclei of G. pyriformis, Mollicutes-like bacterial endosymbionts (MRE), and photosynthetically active and dividing cells of Nostoc punctiforme. Recent genome analyses have shed light on the biology of G. pyriformis, but the genome content and biology of its endosymbionts remain unexplored. To fill this gap, we gathered and examined metagenomic data from the bladders of G. pyriformis, where N. punctiforme and MRE are located. This ensures that our analyses are focused on the organs directly involved in the symbiosis. By comparing this data with the genetic information of related cyanobacteria and MREs from other species of Arbuscular Mycorrhizal Fungi, we aimed to reveal the genetic content of these organisms and understand how they interact at a genetic level to establish a symbiotic relationship. Our analyses uncovered significant gene expansions in the Nostoc endosymbiont, particularly in mobile elements and genes potentially involved in xenobiotic degradation. We also confirmed that the MRE of Glomeromycotina are monophyletic and possess a highly streamlined genome. These genomes show dramatic differences in both structure and content, including the presence of enzymes involved in environmental sensing and stress response.

Genome streamlining in Parcubacteria transitioning from soil to groundwater

BackgroundTo better understand the influence of habitat on the genetic content of bacteria, with a focus on members of Candidate Phyla Radiation (CPR) bacteria, we studied the effects of transitioning from soil via seepage waters to groundwater on genomic composition of ultra-small Parcubacteria, the dominating CPR class in seepage waters, using genome resolved metagenomics.ResultsBacterial metagenome-assembled genomes (MAGs), (318 total, 32 of Parcubacteria) were generated from seepage waters and compared directly to groundwater counterparts. The estimated average genome sizes of members of major phyla Proteobacteria, Bacteroidota and Cand. Patescibacteria (Candidate Phyla Radiation – CPR bacteria) were significantly higher in soil-seepage water as compared to their groundwater counterparts. Seepage water Parcubacteria (Paceibacteria) exhibited 1.18-fold greater mean genome size and 2-fold lower mean proportion of pseudogenes than those in groundwater. Bacteroidota and Proteobacteria also showed a similar trend of reduced genomes in groundwater compared to seepage. While exploring gene loss and adaptive gains in closely related CPR lineages in groundwater, we identified a membrane protein, and a lipoglycopeptide resistance gene unique to a seepage Parcubacterium genome. A nitrite reductase gene was also identified and was unique to the groundwater Parcubacteria genomes, likely acquired from other planktonic microbes via horizontal gene transfer.ConclusionsOverall, our data suggest that bacteria in seepage waters, including ultra-small Parcubacteria, have significantly larger genomes and higher metabolic enrichment than their groundwater counterparts, highlighting possible genome streamlining of the latter in response to habitat selection in an oligotrophic environment.

Endosymbioses Have Shaped the Evolution of Biological Diversity and Complexity Time and Time Again.

Life on Earth comprises prokaryotes and a broad assemblage of endosymbioses. The pages of Molecular Biology and Evolution and Genome Biology and Evolution have provided an essential window into how these endosymbiotic interactions have evolved and shaped biological diversity. Here, we provide a current perspective on this knowledge by drawing on decades of revelatory research published in Molecular Biology and Evolution and Genome Biology and Evolution, and insights from the field at large. The accumulated work illustrates how endosymbioses provide hosts with novel phenotypes that allow them to transition between adaptive landscapes to access environmental resources. Such endosymbiotic relationships have shaped and reshaped life on Earth. The early serial establishment of mitochondria and chloroplasts through endosymbioses permitted massive upscaling of cellular energetics, multicellularity, and terrestrial planetary greening. These endosymbioses are also the foundation upon which all later ones are built, including everything from land-plant endosymbioses with fungi and bacteria to nutritional endosymbioses found in invertebrate animals. Common evolutionary mechanisms have shaped this broad range of interactions. Endosymbionts generally experience adaptive and stochastic genome streamlining, the extent of which depends on several key factors (e.g. mode of transmission). Hosts, in contrast, adapt complex mechanisms of resource exchange, cellular integration and regulation, and genetic support mechanisms to prop up degraded symbionts. However, there are significant differences between endosymbiotic interactions not only in how partners have evolved with each other but also in the scope of their influence on biological diversity. These differences are important considerations for predicting how endosymbioses will persist and adapt to a changing planet.

Host specificity of plant‐associated bacteria is negatively associated with genome size and host abundance along a latitudinal gradient

AbstractHost specialization plays a critical role in the ecology and evolution of plant–microbe symbiosis. Theory predicts that host specialization is associated with microbial genome streamlining and is influenced by the abundance of host species, both of which can vary across latitudes, leading to a latitudinal gradient in host specificity. Here, we quantified the host specificity and composition of plant–bacteria symbioses on leaves across 329 tree species spanning a latitudinal gradient. Our analysis revealed a predominance of host‐specialized leaf bacteria. The degree of host specificity was negatively correlated with bacterial genome size and the local abundance of host plants. Additionally, we found an increased host specificity at lower latitudes, aligning with the high prevalence of small bacterial genomes and rare host species in the tropics. These findings underscore the importance of genome streamlining and host abundance in the evolution of host specificity in plant‐associated bacteria along the latitudinal gradient.

Genome-resolved metagenomics reveals the effect of nutrient availability on bacterial genomic properties across 44 European freshwater lakes.

Understanding intricate microbial interactions in the environment is crucial. This is especially true for the relationships between nutrients and bacteria, as phosphorus, nitrogen and organic carbon availability are known to influence bacterial population dynamics. It has been suggested that low nutrient conditions prompt the evolutionary process of genome streamlining. This process helps conserve scarce nutrients and allows for proliferation. Genome streamlining is associated with genomic properties such as %GC content, genes encoding sigma factors, percent coding regions, gene redundancy, and functional shifts in processes like cell motility and ATP binding cassette transporters, among others. The current study aims to unveil the impact of nutrition on the genome size, %GC content, and functional properties of pelagic freshwater bacteria. We do this at finer taxonomic resolutions for many metagenomically characterized communities. Our study confirms the interplay of trophic level and genomic properties. It also highlights that different nutrient types, particularly phosphorus and nitrogen, impact these properties differently. We observed a covariation of functional traits with genome size. Larger genomes exhibit enriched pathways for motility, environmental interaction, and regulatory genes. ABC transporter genes reflect the availability of nutrients in the environment, with small genomes presumably relying more on metabolites from other organisms. We also discuss the distinct strategies different phyla adopt to adapt to oligotrophic environments. The findings contribute to our understanding of genomic adaptations within complex microbial communities.

Active microbial communities facilitate carbon turnover in brine pools found in the deep Southeastern Mediterranean Sea

Discharge of gas-rich brines fuels productive chemosynthetic ecosystems in the deep sea. In these salty, methanic and sulfidic brines, microbial communities adapt to specific niches along the physicochemical gradients. However, the molecular mechanisms that underpin these adaptations are not fully known. Using metagenomics, we investigated the dense (∼106 cell ml−1) microbial communities that occupy small deep-sea brine pools found in the Southeastern Mediterranean Sea (1150 m water depth, ∼22 °C, ∼60 PSU salinity, sulfide, methane, ammonia reaching millimolar levels, and oxygen usually depleted), reaching high productivity rates of 685 μg C L−1 d−1 ex-situ. We curated 266 metagenome-assembled genomes of bacteria and archaea from the several pools and adjacent sediment-water interface, highlighting the dominance of a single Sulfurimonas, which likely fuels its autotrophy using sulfide oxidation or inorganic sulfur disproportionation. This lineage may be dominant in its niche due to genome streamlining, limiting its metabolic repertoire, particularly by using a single variant of sulfide: quinone oxidoreductase. These primary producers co-exist with ANME-2c archaea that catalyze the anaerobic oxidation of methane. Other lineages can degrade the necromass aerobically (Halomonas and Alcanivorax), or anaerobically through fermentation of macromolecules (e.g., Caldatribacteriota, Bipolaricaulia, Chloroflexota, etc). These low-abundance organisms likely support the autotrophs, providing energy-rich H2, and vital organics such as vitamin B12.

Genome-scale community modelling reveals conserved metabolic cross-feedings in epipelagic bacterioplankton communities

Marine microorganisms form complex communities of interacting organisms that influence central ecosystem functions in the ocean such as primary production and nutrient cycling. Identifying the mechanisms controlling their assembly and activities is a major challenge in microbial ecology. Here, we integrated Tara Oceans meta-omics data to predict genome-scale community interactions within prokaryotic assemblages in the euphotic ocean. A global genome-resolved co-activity network revealed a significant number of inter-lineage associations across diverse phylogenetic distances. Identified co-active communities include species displaying smaller genomes but encoding a higher potential for quorum sensing, biofilm formation, and secondary metabolism. Community metabolic modelling reveals a higher potential for interaction within co-active communities and points towards conserved metabolic cross-feedings, in particular of specific amino acids and group B vitamins. Our integrated ecological and metabolic modelling approach suggests that genome streamlining and metabolic auxotrophies may act as joint mechanisms shaping bacterioplankton community assembly in the global ocean surface.

A redox switch allows binding of Fe(II) and Fe(III) ions in the cyanobacterial iron-binding protein FutA from Prochlorococcus

The marine cyanobacterium Prochlorococcus is a main contributor to global photosynthesis, whilst being limited by iron availability. Cyanobacterial genomes generally encode two different types of FutA iron-binding proteins: periplasmic FutA2 ABC transporter subunits bind Fe(III), while cytosolic FutA1 binds Fe(II). Owing to their small size and their economized genome Prochlorococcus ecotypes typically possess a single futA gene. How the encoded FutA protein might bind different Fe oxidation states was previously unknown. Here, we use structural biology techniques at room temperature to probe the dynamic behavior of FutA. Neutron diffraction confirmed four negatively charged tyrosinates, that together with a neutral water molecule coordinate iron in trigonal bipyramidal geometry. Positioning of the positively charged Arg103 side chain in the second coordination shell yields an overall charge-neutral Fe(III) binding state in structures determined by neutron diffraction and serial femtosecond crystallography. Conventional rotation X-ray crystallography using a home source revealed X-ray-induced photoreduction of the iron center with observation of the Fe(II) binding state; here, an additional positioning of the Arg203 side chain in the second coordination shell maintained an overall charge neutral Fe(II) binding site. Dose series using serial synchrotron crystallography and an XFEL X-ray pump-probe approach capture the transition between Fe(III) and Fe(II) states, revealing how Arg203 operates as a switch to accommodate the different iron oxidation states. This switching ability of the Prochlorococcus FutA protein may reflect ecological adaptation by genome streamlining and loss of specialized FutA proteins.

Genome streamlining to improve performance of a fast-growing cyanobacterium Synechococcus elongatus UTEX 2973.

Cyanobacteria are photosynthetic organisms that have garnered significant recognition as potential hosts for sustainable bioproduction. However, their complex regulatory networks pose significant challenges to major metabolic engineering efforts, thereby limiting their feasibility as production hosts. Genome streamlining has been demonstrated to be a successful approach for improving productivity and fitness in heterotrophs but is yet to be explored to its full potential in phototrophs. Here, we present the systematic reduction of the genome of the cyanobacterium exhibiting the fastest exponential growth, Synechococcus elongatus UTEX 2973. This work, the first of its kind in a photoautotroph, involved an iterative process using state-of-the-art genome-editing technology guided by experimental analysis and computational tools. CRISPR-Cas3 enabled large, progressive deletions of predicted dispensable regions and aided in the identification of essential genes. The large deletions were combined to obtain a strain with 55-kb genome reduction. The strains with streamlined genome showed improvement in growth (up to 23%) and productivity (by 22.7%) as compared to the wild type (WT). This streamlining strategy not only has the potential to develop cyanobacterial strains with improved growth and productivity traits but can also facilitate a better understanding of their genome-to-phenome relationships.IMPORTANCEGenome streamlining is an evolutionary strategy used by natural living systems to dispense unnecessary genes from their genome as a mechanism to adapt and evolve. While this strategy has been successfully borrowed to develop synthetic heterotrophic microbial systems with desired phenotype, it has not been extensively explored in photoautotrophs. Genome streamlining strategy incorporates both computational predictions to identify the dispensable regions and experimental validation using genome-editing tool, and in this study, we have employed a modified strategy with the goal to minimize the genome size to an extent that allows optimal cellular fitness under specified conditions. Our strategy has explored a novel genome-editing tool in photoautotrophs, which, unlike other existing tools, enables large, spontaneous optimal deletions from the genome. Our findings demonstrate the effectiveness of this modified strategy in obtaining strains with streamlined genome, exhibiting improved fitness and productivity.

PoMeLo: a systematic computational approach to predicting metabolic loss in pathogen genomes

BackgroundGenome streamlining, the process by which genomes become smaller and encode fewer genes over time, is a common phenomenon among pathogenic bacteria. This reduction is driven by selection for minimized energy expenditure in a nutrient-rich environment. As pathogens evolve to become more reliant on the host, metabolic genes and resulting capabilities are lost in favor of siphoning metabolites from the host. Characterizing genome streamlining, gene loss, and metabolic pathway degradation can be useful in assessing pathogen dependency on host metabolism and identifying potential targets for host-directed therapeutics.ResultsPoMeLo (Predictor of Metabolic Loss) is a novel evolutionary genomics-guided computational approach for identifying metabolic gaps in the genomes of pathogenic bacteria. PoMeLo leverages a centralized public database of high-quality genomes and annotations and allows the user to compare an unlimited number of genomes across individual genes and pathways. PoMeLo runs locally using user-friendly prompts in a matter of minutes and generates tabular and visual outputs for users to compare predicted metabolic capacity between groups of bacteria and individual species. Each pathway is assigned a Predicted Metabolic Loss (PML) score to assess the magnitude of genome streamlining. Optionally, PoMeLo places the results in an evolutionary context by including phylogenetic relationships in visual outputs. It can also initially compute phylogenetically-weighted mean genome sizes to identify genome streamlining events. Here, we describe PoMeLo and demonstrate its use in identifying metabolic gaps in genomes of pathogenic Treponema species.ConclusionsPoMeLo represents an advance over existing methods for identifying metabolic gaps in genomic data, allowing comparison across large numbers of genomes and placing the resulting data in a phylogenetic context. PoMeLo is freely available for academic and non-academic use at https://github.com/czbiohub-sf/pomelo.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size.

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10k - 100k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Functional annotation of a divergent genome using sequence and structure-based similarity

BackgroundMicrosporidia are a large taxon of intracellular pathogens characterized by extraordinarily streamlined genomes with unusually high sequence divergence and many species-specific adaptations. These unique factors pose challenges for traditional genome annotation methods based on sequence similarity. As a result, many of the microsporidian genomes sequenced to date contain numerous genes of unknown function. Recent innovations in rapid and accurate structure prediction and comparison, together with the growing amount of data in structural databases, provide new opportunities to assist in the functional annotation of newly sequenced genomes.ResultsIn this study, we established a workflow that combines sequence and structure-based functional gene annotation approaches employing a ChimeraX plugin named ANNOTEX (Annotation Extension for ChimeraX), allowing for visual inspection and manual curation. We employed this workflow on a high-quality telomere-to-telomere sequenced tetraploid genome of Vairimorpha necatrix. First, the 3080 predicted protein-coding DNA sequences, of which 89% were confirmed with RNA sequencing data, were used as input. Next, ColabFold was used to create protein structure predictions, followed by a Foldseek search for structural matching to the PDB and AlphaFold databases. The subsequent manual curation, using sequence and structure-based hits, increased the accuracy and quality of the functional genome annotation compared to results using only traditional annotation tools. Our workflow resulted in a comprehensive description of the V. necatrix genome, along with a structural summary of the most prevalent protein groups, such as the ricin B lectin family. In addition, and to test our tool, we identified the functions of several previously uncharacterized Encephalitozoon cuniculi genes.ConclusionWe provide a new functional annotation tool for divergent organisms and employ it on a newly sequenced, high-quality microsporidian genome to shed light on this uncharacterized intracellular pathogen of Lepidoptera. The addition of a structure-based annotation approach can serve as a valuable template for studying other microsporidian or similarly divergent species.

Evolving understanding of rumen methanogen ecophysiology.

Production of methane by methanogenic archaea, or methanogens, in the rumen of ruminants is a thermodynamic necessity for microbial conversion of feed to volatile fatty acids, which are essential nutrients for the animals. On the other hand, methane is a greenhouse gas and its production causes energy loss for the animal. Accordingly, there are ongoing efforts toward developing effective strategies for mitigating methane emissions from ruminant livestock that require a detailed understanding of the diversity and ecophysiology of rumen methanogens. Rumen methanogens evolved from free-living autotrophic ancestors through genome streamlining involving gene loss and acquisition. The process yielded an oligotrophic lifestyle, and metabolically efficient and ecologically adapted descendants. This specialization poses serious challenges to the efforts of obtaining axenic cultures of rumen methanogens, and consequently, the information on their physiological properties remains in most part inferred from those of their non-rumen representatives. This review presents the current knowledge of rumen methanogens and their metabolic contributions to enteric methane production. It also identifies the respective critical gaps that need to be filled for aiding the efforts to mitigate methane emission from livestock operations and at the same time increasing the productivity in this critical agriculture sector.