Primary endosymbiosis events date to the later Proterozoic with cross-calibrated phylogenetic dating of duplicated ATPase proteins

Chloroplasts and mitochondria descended from bacterial ancestors, but the dating of these primary endosymbiosis events remains very uncertain, despite their importance for our understanding of the evolution of both bacteria and eukaryotes. All phylogenetic dating in the Proterozoic and before is difficult: Significant debates surround potential fossil calibration points based on the interpretation of the Precambrian microbial fossil record, and strict molecular clock methods cannot be expected to yield accurate dates over such vast timescales because of strong heterogeneity in rates. Even with more sophisticated relaxed-clock analyses, nodes that are distant from fossil calibrations will have a very high uncertainty in dating. However, endosymbiosis events and gene duplications provide some additional information that has never been exploited in dating; namely, that certain nodes on a gene tree must represent the same events, and thus must have the same or very similar dates, even if the exact date is uncertain. We devised techniques to exploit this information: cross-calibration, in which node date calibrations are reused across a phylogeny, and cross-bracing, in which node date calibrations are formally linked in a hierarchical Bayesian model. We apply these methods to proteins with ancient duplications that have remained associated and originated from plastid and mitochondrial endosymbionts: the α and β subunits of ATP synthase and its relatives, and the elongation factor thermo unstable. The methods yield reductions in dating uncertainty of 14–26% while only using date calibrations derived from phylogenetically unambiguous Phanerozoic fossils of multicellular plants and animals. Our results suggest that primary plastid endosymbiosis occurred ∼900 Mya and mitochondrial endosymbiosis occurred ∼1,200 Mya

Phototrophy and Carbon Fixation in Chlorobi Postdate the Rise of Oxygen

While most productivity on the surface of the Earth today is fueled by oxygenic photosynthesis, for much of Earth history it is thought that anoxygenic photosynthesis—using compounds like ferrous iron or sulfide as electron donors—drove most global carbon fixation. Anoxygenic photosynthesis is still performed by diverse bacteria in niche environments today. Of these, the Chlorobi (formerly green sulfur bacteria) are often interpreted as being particularly ancient and are frequently proposed to have fueled the biosphere during late Archean and early Paleoproterozoic time before the rise of oxygenic photosynthesis. Here, we perform comparative genomic, phylogenetic, and molecular clock analyses to determine the antiquity of the Chlorobi and their characteristic phenotypes. We show that contrary to common assumptions, the Chlorobi clade is relatively young, with anoxygenic phototrophy, carbon fixation via the rTCA pathway, and iron oxidation all significantly postdating the rise of oxygen ~2.3 billion years ago. The Chlorobi therefore could not have fueled the Archean biosphere, but instead represent a relatively young radiation of organisms which likely acquired the capacity for anoxygenic photosynthesis and other traits via horizontal gene transfer sometime after the evolution of oxygenic Cyanobacteria.1

Evolution of Phototrophy in the Chloroflexi Phylum Driven by Horizontal Gene Transfer

The evolutionary mechanisms behind the extant distribution of photosynthesis is a point of substantial contention. Hypotheses range from the presence of phototrophy in the last universal common ancestor and massive gene loss in most lineages, to a later origin in Cyanobacteria followed by extensive horizontal gene transfer into the extant phototrophic clades, with intermediate scenarios that incorporate aspects of both end-members. Here, we report draft genomes of 11 Chloroflexi: the phototrophic Chloroflexia isolate Kouleothrix aurantiaca as well as 10 genome bins recovered from metagenomic sequencing of microbial mats found in Japanese hot springs. Two of these metagenome bins encode photrophic reaction centers and several of these bins form a metabolically diverse, monophyletic clade sister to the Anaerolineae class that we term Candidatus Thermofonsia. Comparisons of organismal (based on conserved ribosomal) and phototrophy (reaction center and bacteriochlorophyll synthesis) protein phylogenies throughout the Chloroflexi demonstrate that two new lineages acquired phototrophy independently via horizontal gene transfer (HGT) from different ancestral donors within the classically phototrophic Chloroflexia class. These results illustrate a complex history of phototrophy within this group, with metabolic innovation tied to HGT. These observations do not support simple hypotheses for the evolution of photosynthesis that require massive character loss from many clades; rather, HGT appears to be the defining mechanic for the distribution of phototrophy in many of the extant clades in which it appears

Bayesian analysis of congruence of core genes in Prochlorococcus and Synechococcus and implications on horizontal gene transfer

It is often suggested that horizontal gene transfer is so ubiquitous in microbes that the concept of a phylogenetic tree representing the pattern of vertical inheritance is oversimplified or even positively misleading. “Universal proteins” have been used to infer the organismal phylogeny, but have been criticized as being only the “tree of one percent.” Currently, few options exist for those wishing to rigorously assess how well a universal protein phylogeny, based on a relative handful of well-conserved genes, represents the phylogenetic histories of hundreds of genes. Here, we address this problem by proposing a visualization method and a statistical test within a Bayesian framework. We use the genomes of marine cyanobacteria, a group thought to exhibit substantial amounts of HGT, as a test case. We take 379 orthologous gene families from 28 cyanobacteria genomes and estimate the Bayesian posterior distributions of trees – a “treecloud” – for each, as well as for a concatenated dataset based on putative “universal proteins.” We then calculate the average distance between trees within and between all treeclouds on various metrics and visualize this high-dimensional space with non-metric multidimensional scaling (NMMDS). We show that the tree space is strongly clustered and that the universal protein treecloud is statistically significantly closer to the center of this tree space than any individual gene treecloud. We apply several commonly-used tests for incongruence/HGT and show that they agree HGT is rare in this dataset, but make different choices about which genes were subject to HGT. Our results show that the question of the representativeness of the “tree of one percent” is a quantitative empirical question, and that the phylogenetic central tendency is a meaningful observation even if many individual genes disagree due to the various sources of incongruence

MetaPOAP: Presence or Absence of Metabolic Pathways in Metagenome-Assembled Genomes

Genome-resolved metagenomics allows the construction of draft microbial genomes from short-read shotgun metagenomics (Metagenome-Assembled Genomes, or MAGs); however, even high-quality MAGs are typically somewhat incomplete and contain a small amount of contaminant sequence, making accurate prediction of metabolic potential challenging. Here, we describe MetaPOAP, an algorithm for probabalistic assessment of the statistical likelihoods for the presence or absence of metabolic pathways in MAGs

Phototrophic Methane Oxidation in a Member of the Chloroflexi Phylum

Biological methane cycling plays an important role in Earth's climate and the global carbon cycle, with biological methane oxidation (methanotrophy) modulating methane release from numerous environments including soils, sediments, and water columns. Methanotrophy is typically coupled to aerobic respiration or anaerobically via the reduction of sulfate, nitrate, or metal oxides, and while the possibility of coupling methane oxidation to phototrophy (photomethanotrophy) has been proposed, no organism has ever been described that is capable of this metabolism. Here we described a new bacterial genome from a member of the Chloroflexi phylum--termed here Candidatus Chlorolinea photomethanotrophicum--with cooccurring methanotrophy and phototrophy pathways, suggesting a novel link between these two metabolisms. Recovered as a metagenome-assembled genome from microbial mats in an iron-rich hot spring in Japan, Ca. "C. photomethanotrophicum" forms a new lineage within the Chloroflexi phylum and expands the known metabolic diversity of this already diverse clade. Ca. "C. photomethanotrophicum" appears to be metabolically versatile, capable of phototrophy (via a Type 2 reaction center), aerobic respiration, nitrite reduction, oxidation of methane and carbon monoxide, and potentially carbon fixation via a novel pathway composed of hybridized components of the serine cycle and the 3-hydroxypropionate bicycle. The biochemical network of this organism is constructed from components from multiple organisms and pathways, further demonstrating the modular nature of metabolic machinery and the ecological and evolutionary importance of horizontal gene transfer in the establishment of novel pathways

Bacterial diversification through geological time

Numerous studies have estimated plant and animal diversification dynamics; however, no comparable rigorous estimates exist for bacteria—the most ancient and widespread form of life on Earth. Here, we analyse phylogenies comprising up to 448,112 bacterial lineages to reconstruct global bacterial diversification dynamics. To handle such large phylogenies, we developed methods based on the statistical properties of infinitely large trees. We further analysed sequencing data from 60 environmental studies to determine the fraction of extant bacterial diversity missing from the phylogenies—a crucial parameter for estimating speciation and extinction rates. We estimate that there are about 1.4–1.9 million extant bacterial lineages when lineages are defined by 99% similarity in the 16S ribosomal RNA gene, and that bacterial diversity has been continuously increasing over the past 1 billion years (Gyr). Recent bacterial extinction rates are estimated at 0.03–0.05 per lineage per million years (lineage^(–1) Myr^(–1)), and are only slightly below estimated recent bacterial speciation rates. Most bacterial lineages ever to have inhabited this planet are estimated to be extinct. Our findings disprove the notion that bacteria are unlikely to go extinct, and provide a valuable perspective on the evolutionary history of a domain of life with a sparse and cryptic fossil record

MetaPOAP: Presence or Absence of Metabolic Pathways in Metagenome-Assembled Genomes

Genome-resolved metagenomics allows the construction of draft microbial genomes from short-read shotgun metagenomics (Metagenome-Assembled Genomes, or MAGs); however, even high-quality MAGs are typically somewhat incomplete and contain a small amount of contaminant sequence, making accurate prediction of metabolic potential challenging. Here, we describe MetaPOAP, an algorithm for probabalistic assessment of the statistical likelihoods for the presence or absence of metabolic pathways in MAGs

Ice-Cap: A Method for Growing Arabidopsis and Tomato Plants in 96-well Plates for High-Throughput Genotyping

It is becoming common for plant scientists to develop projects that require the genotyping of large numbers of plants. The first step in any genotyping project is to collect a tissue sample from each individual plant. The traditional approach to this task is to sample plants one-at-a-time. If one wishes to genotype hundreds or thousands of individuals, however, using this strategy results in a significant bottleneck in the genotyping pipeline. The Ice-Cap method that we describe here provides a high-throughput solution to this challenge by allowing one scientist to collect tissue from several thousand seedlings in a single day 1,2. This level of throughput is made possible by the fact that tissue is harvested from plants 96-at-a-time, rather than one-at-a-time