Strong genome-wide selection early in the evolution of Prochlorococcus resulted in a reduced genome through the loss of a large number of small effect genes

Genomic Insights on the Carbon-Negative Workhorse: Systematical Comparative Genomic Analysis on 56 Synechococcus Strains

Synechococcus, a type of ancient photosynthetic cyanobacteria, is crucial in modern carbon-negative synthetic biology due to its potential for producing bioenergy and high-value products. With its high biomass, fast growth rate, and established genetic manipulation tools, Synechococcus has become a research focus in recent years. Abundant germplasm resources have been accumulated from various habitats, including temperature and salinity conditions relevant to industrialization. In this study, a comprehensive analysis of complete genomes of the 56 Synechococcus strains currently available in public databases was performed, clarifying genetic relationships, the adaptability of Synechococcus to the environment, and its reflection at the genomic level. This was carried out via pan-genome analysis and a detailed comparison of the functional gene groups. The results revealed an open-genome pattern, with 275 core genes and variable genome sizes within these strains. The KEGG annotation and orthology composition comparisons unveiled that the cold and thermophile strains have 32 and 84 unique KO functional units in their shared core gene functional units, respectively. Each KO functional unit reflects unique gene families and pathways. In terms of salt tolerance and comparative genomics, there are 65 unique KO functional units in freshwater-adapted strains and 154 in strictly marine strains. By delving into these aspects, our understanding of the metabolic potential of Synechococcus was deepened, promoting the development and industrial application of cyanobacterial biotechnology.

Regulatory and metabolic adaptations in the nitrogen assimilation of marine picocyanobacteria

Prochlorococcus and Synechococcus are the two most abundant photosynthetic organisms on Earth, with a strong influence on the biogeochemical carbon and nitrogen cycles. Early reports demonstrated the streamlining of regulatory mechanisms in nitrogen metabolism and the removal of genes not strictly essential. The availability of a large series of genomes, and the utilization of latest generation molecular techniques have allowed elucidating the main mechanisms developed by marine picocyanobacteria to adapt to the environments where they thrive, with a particular interest in the strains inhabiting oligotrophic oceans. Given that nitrogen is often limited in those environments, a series of studies have explored the strategies utilized by Prochlorococcus and Synechococcus to exploit the low concentrations of nitrogen-containing molecules available in large areas of the oceans. These strategies include the reduction in the GC and the cellular protein contents; the utilization of truncated proteins; a reduced average amount of N in the proteome; the development of metabolic mechanisms to perceive and utilize nanomolar nitrate concentrations; and the reduced responsiveness of key molecular regulatory systems such as NtcA to 2-oxoglutarate. These findings are in sharp contrast with the large body of knowledge obtained in freshwater cyanobacteria. We will outline the main discoveries, stressing their relevance to the ecological success of these important microorganisms.

Long-Read-Resolved, Ecosystem-Wide Exploration of Nucleotide and Structural Microdiversity of Lake Bacterioplankton Genomes

Reconstruction of metagenome-assembled genomes (MAGs) has become a fundamental approach in microbial ecology. However, a MAG is hardly complete and overlooks genomic microdiversity because metagenomic assembly fails to resolve microvariants among closely related genotypes. Aiming at understanding the universal factors that drive or constrain prokaryotic genome diversification, we performed an ecosystem-wide high-resolution metagenomic exploration of microdiversity by combining spatiotemporal (2 depths × 12 months) sampling from a pelagic freshwater system, high-quality MAG reconstruction using long- and short-read metagenomic sequences, and profiling of single nucleotide variants (SNVs) and structural variants (SVs) through mapping of short and long reads to the MAGs, respectively. We reconstructed 575 MAGs, including 29 circular assemblies, providing high-quality reference genomes of freshwater bacterioplankton. Read mapping against these MAGs identified 100 to 101,781 SNVs/Mb and 0 to 305 insertions, 0 to 467 deletions, 0 to 41 duplications, and 0 to 6 inversions for each MAG. Nonsynonymous SNVs were accumulated in genes potentially involved in cell surface structural modification to evade phage recognition. Most (80.2%) deletions overlapped with a gene coding region, and genes of prokaryotic defense systems were most frequently (>8% of the genes) overlapped with a deletion. Some such deletions exhibited a monthly shift in their allele frequency, suggesting a rapid turnover of genotypes in response to phage predation. MAGs with extremely low microdiversity were either rare or opportunistic bloomers, suggesting that population persistency is key to their genomic diversification. The results concluded that prokaryotic genomic diversification is driven primarily by viral load and constrained by a population bottleneck. IMPORTANCE Identifying intraspecies genomic diversity (microdiversity) is crucial to understanding microbial ecology and evolution. However, microdiversity among environmental assemblages is not well investigated, because most microbes are difficult to culture. In this study, we performed cultivation-independent exploration of bacterial genomic microdiversity in a lake ecosystem using a combination of short- and long-read metagenomic analyses. The results revealed the broad spectrum of genomic microdiversity among the diverse bacterial species in the ecosystem, which has been overlooked by conventional approaches. Our ecosystem-wide exploration further allowed comparative analysis among the genomes and genes and revealed factors behind microbial genomic diversification, namely, that diversification is driven primarily by resistance against viral infection and constrained by the population size.

To Die or Not to Die—Regulated Cell Death and Survival in Cyanobacteria

Regulated cell death (RCD) is central to the development, integrity, and functionality of multicellular organisms. In the last decade, evidence has accumulated that RCD is a universal phenomenon in all life domains. Cyanobacteria are of specific interest due to their importance in aquatic and terrestrial habitats and their role as primary producers in global nutrient cycling. Current knowledge on cyanobacterial RCD is based mainly on biochemical and morphological observations, often by methods directly transferred from vertebrate research and with limited understanding of the molecular genetic basis. However, the metabolism of different cyanobacteria groups relies on photosynthesis and nitrogen fixation, whereas mitochondria are the central executioner of cell death in vertebrates. Moreover, cyanobacteria chosen as biological models in RCD studies are mainly colonial or filamentous multicellular organisms. On the other hand, unicellular cyanobacteria have regulated programs of cellular survival (RCS) such as chlorosis and post-chlorosis resuscitation. The co-existence of different genetically regulated programs in cyanobacterial populations may have been a top engine in life diversification. Development of cyanobacteria-specific methods for identification and characterization of RCD and wider use of single-cell analysis combined with intelligent image-based cell sorting and metagenomics would shed more light on the underlying molecular mechanisms and help us to address the complex colonial interactions during these events. In this review, we focus on the functional implications of RCD in cyanobacterial communities.

Exploring the Diversity of the Thioredoxin Systems in Cyanobacteria

Cyanobacteria evolved the ability to perform oxygenic photosynthesis using light energy to reduce CO₂ from electrons extracted from water and form nutrients. These organisms also developed light-dependent redox regulation through the Trx system, formed by thioredoxins (Trxs) and thioredoxin reductases (TRs). Trxs are thiol-disulfide oxidoreductases that serve as reducing substrates for target enzymes involved in numerous processes such as photosynthetic CO₂ fixation and stress responses. We focus on the evolutionary diversity of Trx systems in cyanobacteria and discuss their phylogenetic relationships. The study shows that most cyanobacteria contain at least one copy of each identified Trx, and TrxA is the only one present in all genomes analyzed. Ferredoxin thioredoxin reductase (FTR) is present in all groups except Gloeobacter and Prochlorococcus, where there is a ferredoxin flavin-thioredoxin reductase (FFTR). Our data suggest that both TRs may have coexisted in ancestral cyanobacteria together with other evolutionarily related proteins such as NTRC or DDOR, probably used against oxidative stress. Phylogenetic studies indicate that they have different evolutionary histories. As cyanobacteria diversified to occupy new habitats, some of these proteins were gradually lost in some groups. Finally, we also review the physiological relevance of redox regulation in cyanobacteria through the study of target enzymes.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Titanium dioxide-mediated fatty acids promoted carotenoid synthesis in Phaffia rhodozyma PR106 analyzed whole genome

Carotenoid, as good colorant and antioxidant, is widely used in the fields of food, medicine and feed. The whole genome of P. rhodozyma PR106 strain with 228.77 mg/L carotenoid (mainly included astaxanthin, β-carotene and lycopene) yield was sequenced, and the genome size was 16.18 Mb, the GC content was 47%. The genetic evolution analysis indicated that PR106 greatly changed in evolution process, and closely related to P. rhodozyma CBS7918. Under 500 mg/L titanium dioxide (TiO₂) stress, carotenoid yield of PR106 was 2.15 times that of the control for 48 h, and was 305.12 mg/L in PR106 to 72 h, interestingly, the yield of oleate, linoleate and α-linolenate also increased significantly among 51 fatty acids by targeted metabolomics analysis. TiO₂ promoted carotenoid synthesis of PR106 by forming astaxanthin esters to reduce the feedback inhibition of carotenoid synthesis. These results provided a theoretical basis for carotenoid production and development using P. rhodozyma.

Comparative Genomics Reveals Insights into Induction of Violacein Biosynthesis and Adaptive Evolution in Janthinobacterium

Violacein has different bioactive properties conferring distinct selective advantages, such as defense from predation and interspecific competition. Adaptation of Janthinobacterium to diverse habitats likely leads to variation in violacein production among phylogenetically closely related species inhabiting different environments, yet genomic mechanisms and the influence of adaptive evolution underpinning violacein biosynthesis in Janthinobacterium are not clear. In this study, we performed genome sequencing, comparative genomic analysis, and phenotypic characterization to investigate genomic factors regulating violacein production in nine Janthinobacterium strains, including a type strain from soil and eight strains we isolated from terrestrial subsurface sediment and groundwater. Results show that although all nine Janthinobacterium strains are phylogenetically closely related and contain genes essential for violacein biosynthesis, they vary in carbon usage and violacein production. Sediment and groundwater strains are weak violacein producers and possess far fewer secondary metabolite biosynthesis genes, indicating genome adaptation compared to soil strains. Further examination suggests that quorum sensing (QS) may play an important role in regulating violacein in Janthinobacterium: the strains exhibiting strong potential in violacein production possess both N-acyl-homoserine lactone (AHL) QS and Janthinobacterium QS (JQS) systems in their genomes, while weaker violacein-producing strains harbor only the JQS system. Preliminary tests of spent media of two Janthinobacterium strains possessing both AHL QS and JQS systems support the potential role of AHLs in inducing violacein production in Janthinobacterium. Overall, results from this study reveal potential genomic mechanisms involved in violacein biosynthesis in Janthinobacterium and provide insights into evolution of Janthinobacterium for adaptation to oligotrophic terrestrial subsurface environment.

IMPORTANCE Phylogenetically closely related bacteria can thrive in diverse environmental habitats due to adaptive evolution. Genomic changes resulting from adaptive evolution lead to variations in cellular function, metabolism, and secondary metabolite biosynthesis. The most well-known secondary metabolite produced by Janthinobacterium is the purple-violet pigment violacein. To date, the mechanisms of induction of violacein biosynthesis in Janthinobacterium is not clear. Comparative genome analysis of closely related Janthinobacterium strains isolated from different environmental habitats not only reveals potential mechanisms involved in induction of violacein production by Janthinobacterium but also provides insights into the survival strategy of Janthinobacterium for adaptation to oligotrophic terrestrial subsurface environment.

Selfishness driving reductive evolution shapes interdependent patterns in spatially structured microbial communities

Microbes release a wide variety of metabolites to the environment that benefit the whole population, called public goods. Public goods sharing drives adaptive function loss, and allows the rise of metabolic cross-feeding. However, how public goods sharing governs the succession of communities over evolutionary time scales remains unclear. To resolve this issue, we constructed an individual-based model, where an autonomous population that possessed functions to produce three essential public goods, was allowed to randomly lose functions. Simulations revealed that function loss genotypes could evolve from the autonomous ancestor, driven by the selfish public production trade-off at the individual level. These genotypes could then automatically develop to three possible types of interdependent patterns: complete functional division, one-way dependency, and asymmetric functional complementation, which were influenced by function cost and function redundancy. In addition, we found random evolutionary events, i.e., the priority and the relative spatial positioning of genotype emergence, are also important in governing community assembly. Moreover, communities occupied by interdependent patterns exhibited better resistance to environmental perturbation, suggesting such patterns are selectively favored. Our work integrates ecological interactions with evolution dynamics, providing a new perspective to explain how reductive evolution shapes microbial interdependencies and governs the succession of communities.

Evolution of Ycf54-independent chlorophyll biosynthesis in cyanobacteria

Chlorophylls (Chls) are essential cofactors for photosynthesis. One of the least understood steps of Chl biosynthesis is formation of the fifth (E) ring, where the red substrate, magnesium protoporphyrin IX monomethyl ester, is converted to the green product, 3,8-divinyl protochlorophyllide a In oxygenic phototrophs, this reaction is catalyzed by an oxygen-dependent cyclase, consisting of a catalytic subunit (AcsF/CycI) and an auxiliary protein, Ycf54. Deletion of Ycf54 impairs cyclase activity and results in severe Chl deficiency, but its exact role is not clear. Here, we used a Δycf54 mutant of the model cyanobacterium Synechocystis sp. PCC 6803 to generate suppressor mutations that restore normal levels of Chl. Sequencing Δycf54 revertants identified a single D219G amino acid substitution in CycI and frameshifts in slr1916, which encodes a putative esterase. Introduction of these mutations to the original Δycf54 mutant validated the suppressor effect, especially in combination. However, comprehensive analysis of the Δycf54 suppressor strains revealed that the D219G-substituted CycI is only partially active and its accumulation is misregulated, suggesting that Ycf54 controls both the level and activity of CycI. We also show that Slr1916 has Chl dephytylase activity in vitro and its inactivation up-regulates the entire Chl biosynthetic pathway, resulting in improved cyclase activity. Finally, large-scale bioinformatic analysis indicates that our laboratory evolution of Ycf54-independent CycI mimics natural evolution of AcsF in low-light-adapted ecotypes of the oceanic cyanobacteria Prochlorococcus, which lack Ycf54, providing insight into the evolutionary history of the cyclase enzyme.

Ecological Genomics and Evolution of Trichoderma reesei

The filamentous fungus Trichoderma reesei (Hypocreales, Ascomycota) is an efficient industrial cell factory for the production of cellulolytic enzymes used for biofuel and other applications. Therefore, researches addressing T. reesei are relatively advanced compared to other Trichoderma spp. because of the significant bulk of available knowledge, multiple genomic data, and gene manipulation techniques. However, the established role of T. reesei in industry has resulted in a frequently biased understanding of the biology of this fungus. Thus, the recent studies unexpectedly show that the superior cellulolytic activity of T. reesei and other Trichoderma species evolved due to multiple lateral gene transfer events, while the innate ability to parasitize other fungi (mycoparasitism) was maintained in the genus, including T. reesei. In this chapter, we will follow the concept of ecological genomics and describe the ecology, distribution, and evolution of T. reesei, as well as critically discuss several common misconceptions that originate from the success of this species in applied sciences and industry.

Combined Methylome, Transcriptome and Proteome Analyses Document Rapid Acclimatization of a Bacterium to Environmental Changes

Polynucleobacter asymbioticus strain QLW-P1DMWA-1^T represents a group of highly successful heterotrophic ultramicrobacteria that is frequently very abundant (up to 70% of total bacterioplankton) in freshwater habitats across all seven continents. This strain was originally isolated from a shallow Alpine pond characterized by rapid changes in water temperature and elevated UV radiation due to its location at an altitude of 1300 m. To elucidate the strain’s adjustment to fluctuating environmental conditions, we recorded changes occurring in its transcriptomic and proteomic profiles under contrasting experimental conditions by simulating thermal conditions in winter and summer as well as high UV irradiation. To analyze the potential connection between gene expression and regulation via methyl group modification of the genome, we also analyzed its methylome. The methylation pattern differed between the three treatments, pointing to its potential role in differential gene expression. An adaptive process due to evolutionary pressure in the genus was deduced by calculating the ratios of non-synonymous to synonymous substitution rates for 20 Polynucleobacter spp. genomes obtained from geographically diverse isolates. The results indicate purifying selection.

Presence of toxin-antitoxin systems in picocyanobacteria and their ecological implications

Picocyanobacteria (mainly Synechococcus and Prochlorococcus) contribute significantly to ocean's primary production. Toxin-Antitoxin (TA) systems present in bacteria and archaea are known to regulate cell growth in response to environmental stresses. However, little is known about the presence of TA systems in picocyanobacteria. This study investigated complete genomes of Synechococcus and Prochlorococcus to understand the prevalence of TA systems in picocyanobacteria. Using the TAfinder software, Type II TA systems were predicted in 27 of 33 (81%) Synechococcus strains, but none of 38 Prochlorococcus strains contain TA genes. Synechococcus strains with larger genomes tend to contain more putative type II TA systems. The number of TA pairs varies from 0 to 42 in Synechococcus strains isolated from various environments. A linear correlation between the genome size and the number of putative TA systems in both coastal and freshwater Synechococcus was established. In general, open ocean Synechococcus contain no or few TA systems, while coastal and freshwater Synechococcus contain more TA systems. The type II TA systems inhibit microbial translation via ribonucleases and allow cells to enter the "dormant" stage in adverse environments. Inheritance of TA genes in freshwater and coastal Synechococcus could confer a recoverable persister mechanism important to survive in variable environments.

Prochlorococcus phage ferredoxin: structural characterization and electron transfer to cyanobacterial sulfite reductases

Marine cyanobacteria are infected by phages whose genomes encode ferredoxin (Fd) electron carriers. These Fds are thought to redirect the energy harvested from light to phage-encoded oxidoreductases that enhance viral fitness, but it is unclear how the biophysical properties and partner specificities of phage Fds relate to those of photosynthetic organisms. Here, results of a bioinformatics analysis using a sequence similarity network revealed that phage Fds are most closely related to cyanobacterial Fds that transfer electrons from photosystems to oxidoreductases involved in nutrient assimilation. Structural analysis of myovirus P-SSM2 Fd (pssm2-Fd), which infects the cyanobacterium Prochlorococcus marinus, revealed high levels of similarity to cyanobacterial Fds (root mean square deviations of ≤0.5 Å). Additionally, pssm2-Fd exhibited a low midpoint reduction potential (-336 mV versus a standard hydrogen electrode), similar to other photosynthetic Fds, although it had lower thermostability (T_m = 28 °C) than did many other Fds. When expressed in an Escherichia coli strain deficient in sulfite assimilation, pssm2-Fd complemented bacterial growth when coexpressed with a P. marinus sulfite reductase, revealing that pssm2-Fd can transfer electrons to a host protein involved in nutrient assimilation. The high levels of structural similarity with cyanobacterial Fds and reactivity with a host sulfite reductase suggest that phage Fds evolved to transfer electrons to cyanobacterially encoded oxidoreductases.

Strong Purifying Selection Is Associated with Genome Streamlining in Epipelagic Marinimicrobia

Marine microorganisms inhabiting nutrient-depleted waters play critical roles in global biogeochemical cycles due to their abundance and broad distribution. Many of these microbes share similar genomic features including small genome size, low % G + C content, short intergenic regions, and low nitrogen content in encoded amino acid residue side chains (N-ARSC), but the evolutionary drivers of these characteristics are unclear. Here, we compared the strength of purifying selection across the Marinimicrobia, a candidate phylum which encompasses a broad range of phylogenetic groups with disparate genomic features, by estimating the ratio of nonsynonymous and synonymous substitutions (dN/dS) in conserved marker genes. Our analysis reveals that epipelagic Marinimicrobia that exhibit features consistent with genome streamlining have significantly lower dN/dS values when compared with their mesopelagic counterparts. We also found a significant positive correlation between median dN/dS values and % G + C content, N-ARSC, and intergenic region length. We did not identify a significant correlation between dN/dS ratios and estimated genome size, suggesting the strength of selection is not a primary factor shaping genome size in this group. Our findings are generally consistent with genome streamlining theory, which postulates that many genomic features of abundant epipelagic bacteria are the result of adaptation to oligotrophic nutrient conditions. Our results are also in agreement with previous findings that genome streamlining is common in epipelagic waters, suggesting that microbes inhabiting this region of the ocean have been shaped by strong selection together with prevalent nutritional constraints characteristic of this environment.

An Expanded Ribosomal Phylogeny of Cyanobacteria Supports a Deep Placement of Plastids

The phylum Cyanobacteria includes free-living bacteria and plastids, the descendants of cyanobacteria that were engulfed by the ancestral lineage of the major photosynthetic eukaryotic group Archaeplastida. Endosymbiotic events that followed this primary endosymbiosis spread plastids across diverse eukaryotic groups. The remnants of the ancestral cyanobacterial genome present in all modern plastids, enable the placement of plastids within Cyanobacteria using sequence-based phylogenetic analyses. To date, such phylogenetic studies have produced conflicting results and two competing hypotheses: (1) plastids diverge relatively recently in cyanobacterial evolution and are most closely related to nitrogen-fixing cyanobacteria, or (2) plastids diverge early in the evolutionary history of cyanobacteria, before the divergence of most cyanobacterial lineages. Here, we use phylogenetic analysis of ribosomal proteins from an expanded data set of cyanobacterial and representative plastid genomes to infer a deep placement for the divergence of the plastid ancestor lineage. We recover plastids as sister to Gloeomargarita and show that the group diverges from other cyanobacterial groups before Pseudanabaena, a previously unreported placement. The tree topologies and phylogenetic distances in our study have implications for future molecular clock studies that aim to model accurate divergence times, especially with respect to groups containing fossil calibrations. The newly sequenced cyanobacterial groups included here will also enable the use of novel cyanobacterial microfossil calibrations.

panX: pan-genome analysis and exploration

Horizontal transfer, gene loss, and duplication result in dynamic bacterial genomes shaped by a complex mixture of different modes of evolution. Closely related strains can differ in the presence or absence of many genes, and the total number of distinct genes found in a set of related isolates—the pan-genome—is often many times larger than the genome of individual isolates. We have developed a pipeline that efficiently identifies orthologous gene clusters in the pan-genome. This pipeline is coupled to a powerful yet easy-to-use web-based visualization for interactive exploration of the pan-genome. The visualization consists of connected components that allow rapid filtering and searching of genes and inspection of their evolutionary history. For each gene cluster, panX displays an alignment, a phylogenetic tree, maps mutations within that cluster to the branches of the tree and infers gain and loss of genes on the core-genome phylogeny. PanX is available at pangenome.de. Custom pan-genomes can be visualized either using a web server or by serving panX locally as a browser-based application.

Evolution and comparative genomics of the most common Trichoderma species

BACKGROUND

The growing importance of the ubiquitous fungal genus Trichoderma (Hypocreales, Ascomycota) requires understanding of its biology and evolution. Many Trichoderma species are used as biofertilizers and biofungicides and T. reesei is the model organism for industrial production of cellulolytic enzymes. In addition, some highly opportunistic species devastate mushroom farms and can become pathogens of humans. A comparative analysis of the first three whole genomes revealed mycoparasitism as the innate feature of Trichoderma. However, the evolution of these traits is not yet understood.

RESULTS

We selected 12 most commonly occurring Trichoderma species and studied the evolution of their genome sequences. Trichoderma evolved in the time of the Cretaceous-Palaeogene extinction event 66 (±15) mya, but the formation of extant sections (Longibrachiatum, Trichoderma) or clades (Harzianum/Virens) happened in Oligocene. The evolution of the Harzianum clade and section Trichoderma was accompanied by significant gene gain, but the ancestor of section Longibrachiatum experienced rapid gene loss. The highest number of genes gained encoded ankyrins, HET domain proteins and transcription factors. We also identified the Trichoderma core genome, completely curated its annotation, investigated several gene families in detail and compared the results to those of other fungi. Eighty percent of those genes for which a function could be predicted were also found in other fungi, but only 67% of those without a predictable function.

CONCLUSIONS

Our study presents a time scaled pattern of genome evolution in 12 Trichoderma species from three phylogenetically distant clades/sections and a comprehensive analysis of their genes. The data offer insights in the evolution of a mycoparasite towards a generalist.

Insights Into the Evolution of Picocyanobacteria and Phycoerythrin Genes (mpeBA and cpeBA)

Marine picocyanobacteria, Prochlorococcus and Synechococcus, substantially contribute to marine primary production and have been the subject of extensive ecological and genomic studies. Little is known about their close relatives from freshwater and non-marine environments. Phylogenomic analyses (using 136 proteins) provide strong support for the monophyly of a clade of non-marine picocyanobacteria consisting of Cyanobium, Synechococcus and marine Sub-cluster 5.2; this clade itself is sister to marine Synechococcus and Prochlorococcus. The most basal lineage within the Syn/Pro clade, Sub-Cluster 5.3, includes marine and freshwater strains. Relaxed molecular clock (SSU, LSU) analyses show that while ancestors of the Syn/Pro clade date as far back as the end of the Pre-Cambrian, modern crown groups evolved during the Carboniferous and Triassic. Comparative genomic analyses reveal novel gene cluster arrangements involved in phycobilisome (PBS) metabolism in freshwater strains. Whilst PBS genes in marine Synechococcus are mostly found in one type of phycoerythrin (PE) rich gene cluster (Type III), strains from non-marine habitats, so far, appear to be more diverse both in terms of pigment content and gene arrangement, likely reflecting a wider range of habitats. Our phylogenetic analyses show that the PE genes (mpeBA) evolved via a duplication of the cpeBA genes in an ancestor of the marine and non-marine picocyanobacteria and of the symbiotic strains Synechococcus spongiarum. A 'primitive' Type III-like ancestor containing cpeBA and mpeBA had thus evolved prior to the divergence of the Syn/Pro clade and S. spongiarum. During the diversification of Synechococcus lineages, losses of mpeBA genes may explain the emergence of pigment cluster Types I, II, IIB, and III in both marine and non-marine habitats, with few lateral gene transfer events in specific taxa.

Ongoing Transposon-Mediated Genome Reduction in the Luminous Bacterial Symbionts of Deep-Sea Ceratioid Anglerfishes

Diverse marine fish and squid form symbiotic associations with extracellular bioluminescent bacteria. These symbionts are typically free-living bacteria with large genomes, but one known lineage of symbionts has undergone genomic reduction and evolution of host dependence. It is not known why distinct evolutionary trajectories have occurred among different luminous symbionts, and not all known lineages previously had genome sequences available. In order to better understand patterns of evolution across diverse bioluminescent symbionts, we de novo sequenced the genomes of bacteria from a poorly studied interaction, the extracellular symbionts from the "lures" of deep-sea ceratioid anglerfishes. Deep-sea anglerfish symbiont genomes are reduced in size by about 50% compared to free-living relatives. They show a striking convergence of genome reduction and loss of metabolic capabilities with a distinct lineage of obligately host-dependent luminous symbionts. These losses include reductions in amino acid synthesis pathways and abilities to utilize diverse sugars. However, the symbiont genomes have retained a number of categories of genes predicted to be useful only outside the host, such as those involved in chemotaxis and motility, suggesting that they may persist in the environment. These genomes contain very high numbers of pseudogenes and show massive expansions of transposable elements, with transposases accounting for 28 and 31% of coding sequences in the symbiont genomes. Transposon expansions appear to have occurred at different times in each symbiont lineage, indicating either independent evolutions of reduction or symbiont replacement. These results suggest ongoing genomic reduction in extracellular luminous symbionts that is facilitated by transposon proliferations.IMPORTANCE Many female deep-sea anglerfishes possess a "lure" containing luminous bacterial symbionts. Here we show that unlike most luminous symbionts, these bacteria are undergoing an evolutionary transition toward small genomes with limited metabolic capabilities. Comparative analyses of the symbiont genomes indicate that this transition is ongoing and facilitated by transposon expansions. This transition may have occurred independently in different symbiont lineages, although it is unclear why. Genomic reduction is common in bacteria that only live within host cells but less common in bacteria that, like anglerfish symbionts, live outside host cells. Since multiple evolutions of genomic reduction have occurred convergently in luminous bacteria, they make a useful system with which to understand patterns of genome evolution in extracellular symbionts. This work demonstrates that ecological factors other than an intracellular lifestyle can lead to dramatic gene loss and evolutionary changes and that transposon expansions may play important roles in this process.

The future of genomics in polar and alpine cyanobacteria

In recent years, genomic analyses have arisen as an exciting way of investigating the functional capacity and environmental adaptations of numerous micro-organisms of global relevance, including cyanobacteria. In the extreme cold of Arctic, Antarctic and alpine environments, cyanobacteria are of fundamental ecological importance as primary producers and ecosystem engineers. While their role in biogeochemical cycles is well appreciated, little is known about the genomic makeup of polar and alpine cyanobacteria. In this article, we present ways that genomic techniques might be used to further our understanding of cyanobacteria in cold environments in terms of their evolution and ecology. Existing examples from other environments (e.g. marine/hot springs) are used to discuss how methods developed there might be used to investigate specific questions in the cryosphere. Phylogenomics, comparative genomics and population genomics are identified as methods for understanding the evolution and biogeography of polar and alpine cyanobacteria. Transcriptomics will allow us to investigate gene expression under extreme environmental conditions, and metagenomics can be used to complement tradition amplicon-based methods of community profiling. Finally, new techniques such as single cell genomics and metagenome assembled genomes will also help to expand our understanding of polar and alpine cyanobacteria that cannot readily be cultured.

Pangenomic Definition of Prokaryotic Species and the Phylogenetic Structure of Prochlorococcus spp

The pangenome is the collection of all groups of orthologous genes (OGGs) from a set of genomes. We apply the pangenome analysis to propose a definition of prokaryotic species based on identification of lineage-specific gene sets. While being similar to the classical biological definition based on allele flow, it does not rely on DNA similarity levels and does not require analysis of homologous recombination. Hence this definition is relatively objective and independent of arbitrary thresholds. A systematic analysis of 110 accepted species with the largest numbers of sequenced strains yields results largely consistent with the existing nomenclature. However, it has revealed that abundant marine cyanobacteria Prochlorococcus marinus should be divided into two species. As a control we have confirmed the paraphyletic origin of Yersinia pseudotuberculosis (with embedded, monophyletic Y. pestis) and Burkholderia pseudomallei (with B. mallei). We also demonstrate that by our definition and in accordance with recent studies Escherichia coli and Shigella spp. are one species.

In it for the long haul: evolutionary consequences of persistent endosymbiosis

Phylogenetically independent bacterial lineages have undergone a profound lifestyle shift: from a free-living to obligately host-associated existence. Among these lineages, intracellular bacterial mutualists of insects are among the most intimate, constrained symbioses known. These obligate endosymbionts exhibit severe gene loss and apparent genome deterioration. Evolutionary theory provides a basis to link their unusual genomic features with shifts in fundamental mechanisms - selection, genetic drift, mutation, and recombination. This mini-review highlights recent comparative and experimental research of processes shaping ongoing diversification within these ancient associations. Recent work supports clear contributions of stochastic processes, including genetic drift and exceptionally strong mutational pressure, toward degenerative evolution. Despite possible compensatory mechanisms, genome degradation may constrain how persistent endosymbionts (and their hosts) respond to environmental fluctuations.

Fundamental differences in diversity and genomic population structure between Atlantic and Pacific Prochlorococcus

The Atlantic and Pacific Oceans represent different biogeochemical regimes in which the abundant marine cyanobacterium Prochlorococcus thrives. We have shown that Prochlorococcus populations in the Atlantic are composed of hundreds of genomically, and likely ecologically, distinct coexisting subpopulations with distinct genomic backbones. Here we ask if differences in the ecology and selection pressures between the Atlantic and Pacific are reflected in the diversity and genomic composition of their indigenous Prochlorococcus populations. We applied large-scale single-cell genomics and compared the cell-by-cell genomic composition of wild populations of co-occurring cells from samples from Station ALOHA off Hawaii, and from Bermuda Atlantic Time Series Station off Bermuda. We reveal fundamental differences in diversity and genomic structure of populations between the sites. The Pacific populations are more diverse than those in the Atlantic, composed of significantly more coexisting subpopulations and lacking dominant subpopulations. Prochlorococcus from the two sites seem to be composed of mostly non-overlapping distinct sets of subpopulations with different genomic backbones-likely reflecting different sets of ocean-specific micro-niches. Furthermore, phylogenetically closely related strains carry ocean-associated nutrient acquisition genes likely reflecting differences in major selection pressures between the oceans. This differential selection, along with geographic separation, clearly has a significant role in shaping these populations.

The effect of extrinsic mortality on genome size evolution in prokaryotes

Mortality has a significant role in prokaryotic ecology and evolution, yet the impact of variations in extrinsic mortality on prokaryotic genome evolution has received little attention. We used both mathematical and agent-based models to reveal how variations in extrinsic mortality affect prokaryotic genome evolution. Our results suggest that the genome size of bacteria increases with increased mortality. A high extrinsic mortality increases the pool of free resources and shortens life expectancy, which selects for faster reproduction, a phenotype we called 'scramblers'. This phenotype is realised by the expansion of gene families involved in nutrient acquisition and metabolism. In contrast, a low mortality rate increases an individual's life expectancy, which results in natural selection favouring tolerance to starvation when conditions are unfavourable. This leads to the evolution of small, streamlined genomes ('stayers'). Our models predict that large genomes, gene family expansion and horizontal gene transfer should be observed in prokaryotes occupying ecosystems exposed to high abiotic stress, as well as those under strong predator- and/or pathogen-mediated selection. A comparison of genome size of cyanobacteria in relatively stable marine versus more turbulent freshwater environments corroborates our predictions, although other factors between these environments could also be responsible.

In silico screening for candidate chassis strains of free fatty acid-producing cyanobacteria

BACKGROUND

Finding a source from which high-energy-density biofuels can be derived at an industrial scale has become an urgent challenge for renewable energy production. Some microorganisms can produce free fatty acids (FFA) as precursors towards such high-energy-density biofuels. In particular, photosynthetic cyanobacteria are capable of directly converting carbon dioxide into FFA. However, current engineered strains need several rounds of engineering to reach the level of production of FFA to be commercially viable; thus new chassis strains that require less engineering are needed. Although more than 120 cyanobacterial genomes are sequenced, the natural potential of these strains for FFA production and excretion has not been systematically estimated.

RESULTS

Here we present the FFA SC (FFASC), an in silico screening method that evaluates the potential for FFA production and excretion of cyanobacterial strains based on their proteomes. A literature search allowed for the compilation of 64 proteins, most of which influence FFA production and a few of which affect FFA excretion. The proteins are classified into 49 orthologous groups (OGs) that helped create rules used in the scoring/ranking of algorithms developed to estimate the potential for FFA production and excretion of an organism. Among 125 cyanobacterial strains, FFASC identified 20 candidate chassis strains that rank in their FFA producing and excreting potential above the specifically engineered reference strain, Synechococcus sp. PCC 7002. We further show that the top ranked cyanobacterial strains are unicellular and primarily include Prochlorococcus (order Prochlorales) and marine Synechococcus (order Chroococcales) that cluster phylogenetically. Moreover, two principal categories of enzymes were shown to influence FFA production the most: those ensuring precursor availability for the biosynthesis of lipids, and those involved in handling the oxidative stress associated to FFA synthesis.

CONCLUSION

To our knowledge FFASC is the first in silico method to screen cyanobacteria proteomes for their potential to produce and excrete FFA, as well as the first attempt to parameterize the criteria derived from genetic characteristics that are favorable/non-favorable for this purpose. Thus, FFASC helps focus experimental evaluation only on the most promising cyanobacteria.

Global genetic capacity for mixotrophy in marine picocyanobacteria

The assimilation of organic nutrients by autotrophs, a form of mixotrophy, has been demonstrated in the globally abundant marine picocyanobacterial genera Prochlorococcus and Synechococcus. However, the range of compounds used and the distribution of organic compound uptake genes within picocyanobacteria are unknown. Here we analyze genomic and metagenomic data from around the world to determine the extent and distribution of mixotrophy in these phototrophs. Analysis of 49 Prochlorococcus and 18 Synechococcus isolate genomes reveals that all have the transporters necessary to take up amino acids, peptides and sugars. However, the number and type of transporters and associated catabolic genes differ between different phylogenetic groups, with low-light IV Prochlorococcus, and 5.1B, 5.2 and 5.3 Synechococcus strains having the largest number. Metagenomic data from 68 stations from the Tara Oceans expedition indicate that the genetic potential for mixotrophy in picocyanobacteria is globally distributed and differs between clades. Phylogenetic analyses indicate gradual organic nutrient transporter gene loss from the low-light IV to the high-light II Prochlorococcus. The phylogenetic differences in genetic capacity for mixotrophy, combined with the ubiquity of picocyanobacterial organic compound uptake genes suggests that mixotrophy has a more central role in picocyanobacterial ecology than was previously thought.

Comparative Genomics of DNA Recombination and Repair in Cyanobacteria: Biotechnological Implications

Cyanobacteria are fascinating photosynthetic prokaryotes that are regarded as the ancestors of the plant chloroplast; the purveyors of oxygen and biomass for the food chain; and promising cell factories for an environmentally friendly production of chemicals. In colonizing most waters and soils of our planet, cyanobacteria are inevitably challenged by environmental stresses that generate DNA damages. Furthermore, many strains engineered for biotechnological purposes can use DNA recombination to stop synthesizing the biotechnological product. Hence, it is important to study DNA recombination and repair in cyanobacteria for both basic and applied research. This review reports what is known in a few widely studied model cyanobacteria and what can be inferred by mining the sequenced genomes of morphologically and physiologically diverse strains. We show that cyanobacteria possess many E. coli-like DNA recombination and repair genes, and possibly other genes not yet identified. E. coli-homolog genes are unevenly distributed in cyanobacteria, in agreement with their wide genome diversity. Many genes are extremely well conserved in cyanobacteria (mutMS, radA, recA, recFO, recG, recN, ruvABC, ssb, and uvrABCD), even in small genomes, suggesting that they encode the core DNA repair process. In addition to these core genes, the marine Prochlorococcus and Synechococcus strains harbor recBCD (DNA recombination), umuCD (mutational DNA replication), as well as the key SOS genes lexA (regulation of the SOS system) and sulA (postponing of cell division until completion of DNA reparation). Hence, these strains could possess an E. coli-type SOS system. In contrast, several cyanobacteria endowed with larger genomes lack typical SOS genes. For examples, the two studied Gloeobacter strains lack alkB, lexA, and sulA; and Synechococcus PCC7942 has neither lexA nor recCD. Furthermore, the Synechocystis PCC6803 lexA product does not regulate DNA repair genes. Collectively, these findings indicate that not all cyanobacteria have an E. coli-type SOS system. Also interestingly, several cyanobacteria possess multiple copies of E. coli-like DNA repair genes, such as Acaryochloris marina MBIC11017 (2 alkB, 3 ogt, 7 recA, 3 recD, 2 ssb, 3 umuC, 4 umuD, and 8 xerC), Cyanothece ATCC51142 (2 lexA and 4 ruvC), and Nostoc PCC7120 (2 ssb and 3 xerC).

Lack of Overt Genome Reduction in the Bryostatin-Producing Bryozoan Symbiont "Candidatus Endobugula sertula"

The uncultured bacterial symbiont "Candidatus Endobugula sertula" is known to produce cytotoxic compounds called bryostatins, which protect the larvae of its host, Bugula neritina The symbiont has never been successfully cultured, and it was thought that its genome might be significantly reduced. Here, we took a shotgun metagenomics and metatranscriptomics approach to assemble and characterize the genome of "Ca Endobugula sertula." We found that it had specific metabolic deficiencies in the biosynthesis of certain amino acids but few other signs of genome degradation, such as small size, abundant pseudogenes, and low coding density. We also identified homologs to genes associated with insect pathogenesis in other gammaproteobacteria, and these genes may be involved in host-symbiont interactions and vertical transmission. Metatranscriptomics revealed that these genes were highly expressed in a reproductive host, along with bry genes for the biosynthesis of bryostatins. We identified two new putative bry genes fragmented from the main bry operon, accounting for previously missing enzymatic functions in the pathway. We also determined that a gene previously assigned to the pathway, bryS, is not expressed in reproductive tissue, suggesting that it is not involved in the production of bryostatins. Our findings suggest that "Ca Endobugula sertula" may be able to live outside the host if its metabolic deficiencies are alleviated by medium components, which is consistent with recent findings that it may be possible for "Ca Endobugula sertula" to be transmitted horizontally.

IMPORTANCE

The bryostatins are potent protein kinase C activators that have been evaluated in clinical trials for a number of indications, including cancer and Alzheimer's disease. There is, therefore, considerable interest in securing a renewable supply of these compounds, which is currently only possible through aquaculture of Bugula neritina and total chemical synthesis. However, these approaches are labor-intensive and low-yielding and thus preclude the use of bryostatins as a viable therapeutic agent. Our genome assembly and transcriptome analysis for "Ca Endobugula sertula" shed light on the metabolism of this symbiont, potentially aiding isolation and culturing efforts. Our identification of additional bry genes may also facilitate efforts to express the complete pathway heterologously.

Losing Complexity: The Role of Simplification in Macroevolution

Macroevolutionary patterns can be produced by combinations of diverse and even oppositional dynamics. A growing body of data indicates that secondary simplifications of molecular and cellular structures are common. Some major diversifications in eukaryotes have occurred because of loss and minimalisation; numerous episodes in prokaryote evolution have likewise been driven by the reduction of structure. After examining a range of examples of secondary simplification and its consequences across the tree of life, we address how macroevolutionary explanations might incorporate simplification as well as complexification, and adaptive as well as nonadaptive dynamics.

Endosymbiont evolution: predictions from theory and surprises from genomes

Genome data have created new opportunities to untangle evolutionary processes shaping microbial variation. Among bacteria, long-term mutualists of insects represent the smallest and (typically) most AT-rich genomes. Evolutionary theory provides a context to predict how an endosymbiotic lifestyle may alter fundamental evolutionary processes--mutation, selection, genetic drift, and recombination--and thus contribute to extreme genomic outcomes. These predictions can then be explored by comparing evolutionary rates, genome size and stability, and base compositional biases across endosymbiotic and free-living bacteria. Recent surprises from such comparisons include genome reduction among uncultured, free-living species. Some studies suggest that selection generally drives this streamlining, while drift drives genome reduction in endosymbionts; however, this remains an hypothesis requiring additional data. Unexpected evidence of selection acting on endosymbiont GC content hints that even weak selection may be effective in some long-term mutualists. Moving forward, intraspecific analysis offers a promising approach to distinguish underlying mechanisms, by testing the null hypothesis of neutrality and by quantifying mutational spectra. Such analyses may clarify whether endosymbionts and free-living bacteria occupy distinct evolutionary trajectories or, alternatively, represent varied outcomes of similar underlying forces.

Genomes of diverse isolates of the marine cyanobacterium Prochlorococcus

The marine cyanobacterium Prochlorococcus is the numerically dominant photosynthetic organism in the oligotrophic oceans, and a model system in marine microbial ecology. Here we report 27 new whole genome sequences (2 complete and closed; 25 of draft quality) of cultured isolates, representing five major phylogenetic clades of Prochlorococcus. The sequenced strains were isolated from diverse regions of the oceans, facilitating studies of the drivers of microbial diversity—both in the lab and in the field. To improve the utility of these genomes for comparative genomics, we also define pre-computed clusters of orthologous groups of proteins (COGs), indicating how genes are distributed among these and other publicly available Prochlorococcus genomes. These data represent a significant expansion of Prochlorococcus reference genomes that are useful for numerous applications in microbial ecology, evolution and oceanography.