Snowball Earth, population bottleneck and Prochlorococcus evolution

Physical constraints during Snowball Earth drive the evolution of multicellularity

Molecular and fossil evidence suggests that complex eukaryotic multicellularity evolved during the late Neoproterozoic era, coincident with Snowball Earth glaciations, where ice sheets covered most of the globe. During this period, environmental conditions-such as seawater temperature and the availability of photosynthetically active light in the oceans-likely changed dramatically. Such changes would have had significant effects on both resource availability and optimal phenotypes. Here, we construct and apply mechanistic models to explore (i) how environmental changes during Snowball Earth and biophysical constraints generated selective pressures, and (ii) how these pressures may have had differential effects on organisms with different forms of biological organization. By testing a series of alternative-and commonly debated-hypotheses, we demonstrate how multicellularity was likely acquired differently in eukaryotes and prokaryotes owing to selective differences on their size due to the biophysical and metabolic regimes they inhabit: decreasing temperatures and resource availability instigated by the onset of glaciations generated selective pressures towards smaller sizes in organisms in the diffusive regime and towards larger sizes in motile heterotrophs. These results suggest that changing environmental conditions during Snowball Earth glaciations gave multicellular eukaryotes an evolutionary advantage, paving the way for the complex multicellular lineages that followed.

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 10 k - 100 k 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.

A timeline of bacterial and archaeal diversification in the ocean

Microbial plankton play a central role in marine biogeochemical cycles, but the timing in which abundant lineages diversified into ocean environments remains unclear. Here, we reconstructed the timeline in which major clades of bacteria and archaea colonized the ocean using a high-resolution benchmarked phylogenetic tree that allows for simultaneous and direct comparison of the ages of multiple divergent lineages. Our findings show that the diversification of the most prevalent marine clades spans throughout a period of 2.2 Ga, with most clades colonizing the ocean during the last 800 million years. The oldest clades - SAR202, SAR324, Ca. Marinimicrobia, and Marine Group II - diversified around the time of the Great Oxidation Event, during which oxygen concentration increased but remained at microaerophilic levels throughout the Mid-Proterozoic, consistent with the prevalence of some clades within these groups in oxygen minimum zones today. We found the diversification of the prevalent heterotrophic marine clades SAR11, SAR116, SAR92, SAR86, and Roseobacter as well as the Marine Group I to occur near to the Neoproterozoic Oxygenation Event (0.8-0.4 Ga). The diversification of these clades is concomitant with an overall increase of oxygen and nutrients in the ocean at this time, as well as the diversification of eukaryotic algae, consistent with the previous hypothesis that the diversification of heterotrophic bacteria is linked to the emergence of large eukaryotic phytoplankton. The youngest clades correspond to the widespread phototrophic clades Prochlorococcus, Synechococcus, and Crocosphaera, whose diversification happened after the Phanerozoic Oxidation Event (0.45-0.4 Ga), in which oxygen concentrations had already reached their modern levels in the atmosphere and the ocean. Our work clarifies the timing at which abundant lineages of bacteria and archaea colonized the ocean, thereby providing key insights into the evolutionary history of lineages that comprise the majority of prokaryotic biomass in the modern ocean.

Genomic Legacies of Ancient Adaptation Illuminate GC-Content Evolution in Bacteria

Bacterial evolution is characterized by strong purifying selection as well as rapid adaptive evolution in changing environments. In this context, the genomic GC content (genomic GC) varies greatly but presents some level of phylogenetic stability, making it challenging to explain based on current hypotheses. To illuminate the evolutionary mechanisms of the genomic GC, we analyzed the base composition and functional inventory of 11,083 representative genomes. A phylogenetically constrained bimodal distribution of the genomic GC, which mainly originated from parallel divergences in the early evolution, was demonstrated. Such variation of the genomic GC can be well explained by DNA replication and repair (DRR), in which multiple pathways correlate with the genomic GC. Furthermore, the biased conservation of various stress-related genes, especially the DRR-related ones, implies distinct adaptive processes in the ancestral lineages of high- or low-GC clades which are likely induced by major environmental changes. Our findings support that the mutational biases resulting from these legacies of ancient adaptation have changed the course of adaptive evolution and generated great variation in the genomic GC. This highlights the importance of indirect effects of natural selection, which indicates a new model for bacterial evolution. IMPORTANCE GC content has been shown to be an important factor in microbial ecology and evolution, and the genomic GC of bacteria can be characterized by great intergenomic heterogeneity, high intragenomic homogeneity, and strong phylogenetic inertia, as well as being associated with the environment. Current hypotheses concerning direct selection or mutational biases cannot well explain these features simultaneously. Our findings of the genomic GC showing that ancient adaptations have transformed the DRR system and that the resulting mutational biases further contributed to a bimodal distribution of it offer a more reasonable scenario for the mechanism. This would imply that, when thinking about the evolution of life, diverse processes of adaptation exist, and combined effects of natural selection should be considered.

Assessing a Role of Genetic Drift for Deep-Time Evolutionary Events

Effective population size (N_e) determines the amount of genetic diversity and the fate of genetic variants in a species and thus is an essential parameter in evolutionary genetics. There are standard approaches to determine the N_e of evolving species. For example, the long-term N_e of an extant species is calculated based on its unbiased global mutation rate and the neutral genetic diversity of the species. However, approaches for inferring N_e of ancestral lineages are less known. Here, we introduce an evolutionary genetic statistic and an analytical procedure to assess the efficiency of natural selection for deep nodes by calculating rates of nonsynonymous nucleotide substitutions leading to radical (d_R) and conservative (d_C) amino acid replacements, respectively. Given that radical variants are more likely to be deleterious than conservative ones, an elevated d_R/d_C ratio in gene families across the genome means an accelerated genome-wide accumulation of the more deleterious type of mutations (i.e., radical variants), which indicates that natural selection is less efficient and genetic drift becomes more powerful. Earlier approaches that calculate d_R/d_C do not consider the impact of nucleotide composition (G+C content) on the d_R/d_C result, which is partially accounted for in more recent methods. Here, we use these methods to demonstrate that genetic drift may have driven the early evolution of Prochlorococcus, the most abundant carbon-fixing photosynthetic bacteria in the ocean.

A Computational Protocol for Dating the Evolution of Cyanobacteria

Cyanobacteria are known to play important roles in driving biological and geochemical innovations in ancient Earth. The origin of Cyanobacteria is the key to understanding these evolutionary events and thus has gained much interest to biologists and geobiologists. Recent development of the molecular dating approaches provides us an opportunity to assess the timeline of Cyanobacteria evolution based on relaxed clock models. The implementation of Bayesian phylogenetic approaches accommodates the uncertainties from different sources, such as fossil calibrations and topological structure of the phylogenomic tree, and provides us converged estimates of posterior mean ages. In this chapter, by taking Cyanobacteria as an example, we introduce a refined strategy to perform molecular dating analysis, as well as a practical method to evaluate the precision of dating analysis.