Endosymbiont evolution: predictions from theory and surprises from genomes

A new lysine biosynthetic enzyme from a bacterial endosymbiont shaped by genetic drift and genome reduction

The effect of population bottlenecks and genome reduction on enzyme function is poorly understood. Candidatus Liberibacter solanacearum is a bacterium with a reduced genome that is transmitted vertically to the egg of an infected psyllid-a population bottleneck that imposes genetic drift and is predicted to affect protein structure and function. Here, we define the function of Ca. L. solanacearum dihydrodipicolinate synthase (CLsoDHDPS), which catalyzes the committed branchpoint reaction in diaminopimelate and lysine biosynthesis. We demonstrate that CLsoDHDPS is expressed in Ca. L. solanacearum and expression is increased ~2-fold in the insect host compared to in planta. CLsoDHDPS has decreased thermal stability and increased aggregation propensity, implying mutations have destabilized the enzyme but are compensated for through elevated chaperone expression and a stabilized oligomeric state. CLsoDHDPS uses a ternary-complex kinetic mechanism, which is to date unique among DHDPS enzymes, has unusually low catalytic ability, but an unusually high substrate affinity. Structural studies demonstrate that the active site is more open, and the structure of CLsoDHDPS with both pyruvate and the substrate analogue succinic-semialdehyde reveals that the product is both structurally and energetically different and therefore evolution has in this case fashioned a new enzyme. Our study suggests the effects of genome reduction and genetic drift on the function of essential enzymes and provides insights on bacteria-host co-evolutionary associations. We propose that bacteria with endosymbiotic lifestyles present a rich vein of interesting enzymes useful for understanding enzyme function and/or informing protein engineering efforts.

Genome reduction and horizontal gene transfer in the evolution of Endomicrobia-rise and fall of an intracellular symbiosis with termite gut flagellates

Bacterial endosymbionts of eukaryotic hosts typically experience massive genome reduction, but the underlying evolutionary processes are often obscured by the lack of free-living relatives. Endomicrobia, a family-level lineage of host-associated bacteria in the phylum Elusimicrobiota that comprises both free-living representatives and endosymbionts of termite gut flagellates, are an excellent model to study evolution of intracellular symbionts. We reconstructed 67 metagenome-assembled genomes (MAGs) of Endomicrobiaceae among more than 1,700 MAGs from the gut microbiota of a wide range of termites. Phylogenomic analysis confirmed a sister position of representatives from termites and ruminants, and allowed to propose eight new genera in the radiation of Endomicrobiaceae. Comparative genome analysis documented progressive genome erosion in the new genus Endomicrobiellum, which comprises all flagellate endosymbionts characterized to date. Massive gene losses were accompanied by the acquisition of new functions by horizontal gene transfer, which led to a shift from a glucose-based energy metabolism to one based on sugar phosphates. The breakdown of glycolysis and many anabolic pathways for amino acids and cofactors in several subgroups was compensated by the independent acquisition of new uptake systems, including an ATP/ADP antiporter, from other gut microbiota. The putative donors are mostly flagellate endosymbionts from other bacterial phyla, including several, hitherto unknown lineages of uncultured Alphaproteobacteria, documenting the importance of horizontal gene transfer in the convergent evolution of these intracellular symbioses. The loss of almost all biosynthetic capacities in some lineages of Endomicrobiellum suggests that their originally mutualistic relationship with flagellates is on its decline.IMPORTANCEUnicellular eukaryotes are frequently colonized by bacterial and archaeal symbionts. A prominent example are the cellulolytic gut flagellates of termites, which harbor diverse but host-specific bacterial symbionts that occur exclusively in termite guts. One of these lineages, the so-called Endomicrobia, comprises both free-living and endosymbiotic representatives, which offers the unique opportunity to study the evolutionary processes underpinning the transition from a free-living to an intracellular lifestyle. Our results revealed a progressive gene loss in energy metabolism and biosynthetic pathways, compensated by the acquisition of new functions via horizontal gene transfer from other gut bacteria, and suggest the eventual breakdown of an initially mutualistic symbiosis. Evidence for convergent evolution of unrelated endosymbionts reflects adaptations to the intracellular environment of termite gut flagellates.

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.

Genome evolution in intracellular parasites: Microsporidia and Apicomplexa

Microsporidia and Apicomplexa are eukaryotic, single-celled, intracellular parasites with huge public health and economic importance. Typically, these parasites are studied separately, emphasizing their uniqueness and diversity. In this review, we explore the huge amount of genomic data that has recently become available for the two groups. We compare and contrast their genome evolution and discuss how their transitions to intracellular life may have shaped it. In particular, we explore genome reduction and compaction, genome expansion and ploidy, gene shuffling and rearrangements, and the evolution of centromeres and telomeres.

Exploring the landscape of symbiotic diversity and distribution in unicellular ciliated protists

BACKGROUND

The eukaryotic-bacterial symbiotic system plays an important role in various physiological, developmental, and evolutionary processes. However, our current understanding is largely limited to multicellular eukaryotes without adequate consideration of diverse unicellular protists, including ciliates.

RESULTS

To investigate the bacterial profiles associated with unicellular organisms, we collected 246 ciliate samples spanning the entire Ciliophora phylum and conducted single-cell based metagenome sequencing. This effort has yielded the most extensive collection of bacteria linked to unicellular protists to date. From this dataset, we identified 883 bacterial species capable of cohabiting with ciliates, unveiling the genomes of 116 novel bacterial cohabitants along with 7 novel archaeal cohabitants. Highlighting the intimate relationship between ciliates and their cohabitants, our study unveiled that over 90% of ciliates coexist with bacteria, with individual hosts fostering symbiotic relationships with multiple bacteria concurrently, resulting in the observation of seven distinct symbiotic patterns among bacteria. Our exploration of symbiotic mechanisms revealed the impact of host digestion on the intracellular diversity of cohabitants. Additionally, we identified the presence of eukaryotic-like proteins in bacteria as a potential contributing factor to their resistance against host digestion, thereby expanding their potential host range.

CONCLUSIONS

As the first large-scale analysis of prokaryotic associations with ciliate protists, this study provides a valuable resource for future research on eukaryotic-bacterial symbioses. Video Abstract.

How do bacterial endosymbionts work with so few genes?

The move from a free-living environment to a long-term residence inside a host eukaryotic cell has profound effects on bacterial function. While endosymbioses are found in many eukaryotes, from protists to plants to animals, the bacteria that form these host-beneficial relationships are even more diverse. Endosymbiont genomes can become radically smaller than their free-living relatives, and their few remaining genes show extreme compositional biases. The details of how these reduced and divergent gene sets work, and how they interact with their host cell, remain mysterious. This Unsolved Mystery reviews how genome reduction alters endosymbiont biology and highlights a "tipping point" where the loss of the ability to build a cell envelope coincides with a marked erosion of translation-related genes.

Facultative lifestyle drives diversity of coral algal symbionts

The photosynthetic symbionts of corals sustain biodiverse reefs in nutrient-poor, tropical waters. Recent genomic data illuminate the evolution of coral symbionts under genome size constraints and suggest that retention of the facultative lifestyle, widespread among these algae, confers a selective advantage when compared with a strict symbiotic existence. We posit that the coral symbiosis is analogous to a 'bioreactor' that selects winner genotypes and allows them to rise to high numbers in a sheltered habitat prior to release by the coral host. Our observations lead to a novel hypothesis, the 'stepping-stone model', which predicts that local adaptation under both the symbiotic and free-living stages, in a stepwise fashion, accelerates coral alga diversity and the origin of endemic strains and species.

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.

Comparative genomics of a vertically transmitted thiotrophic bacterial ectosymbiont and its close free-living relative

Thiotrophic symbioses between sulphur-oxidizing bacteria and various unicellular and metazoan eukaryotes are widespread in reducing marine environments. The giant colonial ciliate Zoothamnium niveum, however, is the only host of thioautotrophic symbionts that has been cultivated along with its symbiont, the vertically transmitted ectosymbiont Candidatus Thiobius zoothamnicola (short Thiobius). Because theoretical predictions posit a smaller genome in vertically transmitted endosymbionts compared to free-living relatives, we investigated whether this is true also for an ectosymbiont. We used metagenomics to recover the high-quality draft genome of this bacterial symbiont. For comparison we have also sequenced a closely related free-living cultured but not formally described strain Milos ODIII6 (short ODIII6). We then performed comparative genomics to assess the functional capabilities at gene, metabolic pathway and trait level. 16S rRNA gene trees and average amino acid identity confirmed the close phylogenetic relationship of both bacteria. Indeed, Thiobius has about a third smaller genome than its free-living relative ODIII6, with reduced metabolic capabilities and fewer functional traits. The functional capabilities of Thiobius were a subset of those of the more versatile ODIII6, which possessed additional genes for oxygen, sulphur and hydrogen utilization and for the acquisition of phosphorus illustrating features that may be adaptive for the unstable environmental conditions at hydrothermal vents. In contrast, Thiobius possesses genes potentially enabling it to utilize lactate and acetate heterotrophically, compounds that may be provided as byproducts by the host. The present study illustrates the effect of strict host-dependence of a bacterial ectosymbiont on genome evolution and host adaptation.

Convergent reductive evolution and host adaptation in Mycoavidus bacterial endosymbionts of Mortierellaceae fungi

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.

Contrasting outcomes of genome reduction in mikrocytids and microsporidians

BACKGROUND

Intracellular symbionts often undergo genome reduction, losing both coding and non-coding DNA in a process that ultimately produces small, gene-dense genomes with few genes. Among eukaryotes, an extreme example is found in microsporidians, which are anaerobic, obligate intracellular parasites related to fungi that have the smallest nuclear genomes known (except for the relic nucleomorphs of some secondary plastids). Mikrocytids are superficially similar to microsporidians: they are also small, reduced, obligate parasites; however, as they belong to a very different branch of the tree of eukaryotes, the rhizarians, such similarities must have evolved in parallel. Since little genomic data are available from mikrocytids, we assembled a draft genome of the type species, Mikrocytos mackini, and compared the genomic architecture and content of microsporidians and mikrocytids to identify common characteristics of reduction and possible convergent evolution.

RESULTS

At the coarsest level, the genome of M. mackini does not exhibit signs of extreme genome reduction; at 49.7 Mbp with 14,372 genes, the assembly is much larger and gene-rich than those of microsporidians. However, much of the genomic sequence and most (8075) of the protein-coding genes code for transposons, and may not contribute much of functional relevance to the parasite. Indeed, the energy and carbon metabolism of M. mackini share several similarities with those of microsporidians. Overall, the predicted proteome involved in cellular functions is quite reduced and gene sequences are extremely divergent. Microsporidians and mikrocytids also share highly reduced spliceosomes that have retained a strikingly similar subset of proteins despite having reduced independently. In contrast, the spliceosomal introns in mikrocytids are very different from those of microsporidians in that they are numerous, conserved in sequence, and constrained to an exceptionally narrow size range (all 16 or 17 nucleotides long) at the shortest extreme of known intron lengths.

CONCLUSIONS

Nuclear genome reduction has taken place many times and has proceeded along different routes in different lineages. Mikrocytids show a mix of similarities and differences with other extreme cases, including uncoupling the actual size of a genome with its functional reduction.

From Pathogens to Cancer: Are Cancer Cells Evolved Mitochondrial Super Cells?

Life is based on a highly specific combination of atoms, metabolism, and genetics which eventually reflects the chemistry of the Universe which is composed of hydrogen, oxygen, nitrogen, sulfur, phosphorus, and carbon. The interaction of atomic, metabolic, and genetic cycles results in the organization and de-organization of chemical information of that which we consider as living entities, including cancer cells. In order to approach the problem of the origin of cancer it is therefore reasonable to start from the assumption that the sub-molecular level, the atomic structure, should be the considered starting point on which metabolism, genetics, and external insults eventually emanate. Second, it is crucial to characterize which of the entities and parts composing human cells may live a separate life; certainly, this theoretical standpoint would consider mitochondria, an organelle of "bacteria" origin embedded in conditions favorable for the onset of both. This organelle has not only been tolerated by immunity but has also been placed as a central regulator of cell defense. Virus, bacteria, and mitochondria are also similar in the light of genetic and metabolic elements; they share not only equivalent DNA and RNA features but also many basic biological activities. Thus, it is important to finalize that once the cellular integrity has been constantly broken down, the mitochondria like any other virus or bacteria return to their original autonomy to simply survive. The Warburg's law that states the ability of cancers to ferment glucose in the presence of oxygen, indicates mitochondria respiration abnormalities may be the underlying cause of this transformation towards super cancer cells. Though genetic events play a key part in altering biochemical metabolism, inducing aerobic glycolysis, this is not enough to impair mitochondrial function since mitochondrial biogenesis and quality control are constantly upregulated in cancers. While some cancers have mutations in the nuclear-encoded mitochondrial tricarboxylic acid (TCA) cycle, enzymes that produce oncogenic metabolites, there is also a bio-physic pathway for pathogenic mitochondrial genome mutations. The atomic level of all biological activities can be considered the very beginning, marked by the electron abnormal behavior that consequently affects DNA of both cells and mitochondria. Whilst the cell's nucleus DNA after a certain number of errors and defection tends to gradually switch off, the mitochondria DNA starts adopting several escape strategies, switching-on a few important genes that belong back at their original roots as independent beings. The ability to adopt this survival trick, by becoming completely immune to current life-threatening events, is probably the beginning of a differentiation process towards a "super-power cell", the cancer cells that remind many pathogens, including virus, bacteria, and fungi. Thus, here, we present a hypothesis regarding those changes that first begin at the mitochondria atomic level to steadily involve molecular, tissue and organ levels in response to the virus or bacteria constant insults that drive a mitochondria itself to become an "immortal cancer cell". Improved insights into this interplay between these pathogens and mitochondria progression may disclose newly epistemological paradigms as well as innovative procedures in targeting cancer cell progressive invasion.

The Influence of the Selection at the Amino Acid Level on Synonymous Codon Usage from the Viewpoint of Alternative Genetic Codes

Synonymous codon usage can be influenced by mutations and/or selection, e.g., for speed of protein translation and correct folding. However, this codon bias can also be affected by a general selection at the amino acid level due to differences in the acceptance of the loss and generation of these codons. To assess the importance of this effect, we constructed a mutation-selection model model, in which we generated almost 90,000 stationary nucleotide distributions produced by mutational processes and applied a selection based on differences in physicochemical properties of amino acids. Under these conditions, we calculated the usage of fourfold degenerated (4FD) codons and compared it with the usage characteristic of the pure mutations. We considered both the standard genetic code (SGC) and alternative genetic codes (AGCs). The analyses showed that a majority of AGCs produced a greater 4FD codon bias than the SGC. The mutations producing more thymine or adenine than guanine and cytosine increased the differences in usage. On the other hand, the mutational pressures generating a lot of cytosine or guanine with a low content of adenine and thymine decreased this bias because the nucleotide content of most 4FD codons stayed in the compositional equilibrium with these pressures. The comparison of the theoretical results with those for real protein coding sequences showed that the influence of selection at the amino acid level on the synonymous codon usage cannot be neglected. The analyses indicate that the effect of amino acid selection cannot be disregarded and that it can interfere with other selection factors influencing codon usage, especially in AT-rich genomes, in which AGCs are usually used.

Evolutionary jumps in bacterial GC content

Genomic GC (Guanine-Cytosine) content is a fundamental molecular trait linked with many key genomic features such as codon and amino acid use. Across bacteria, GC content is surprisingly diverse and has been studied for many decades; yet its evolution remains incompletely understood. Since it is difficult to observe GC content evolve on laboratory time scales, phylogenetic comparative approaches are instrumental; but this dimension is rarely studied systematically in the case of bacterial GC content. We applied phylogenetic comparative models to analyze GC content evolution in multiple bacterial groups across 2 major bacterial phyla. We find that GC content diversifies via a combination of gradual evolution and evolutionary “jumps.” Surprisingly, unlike prior reports that solely focused on reductions in GC, we found a comparable number of jumps with both increased and decreased GC content. Overall, many of the identified jumps occur in lineages beyond the well-studied peculiar examples of endosymbiotic and AT-rich marine bacteria and do not support the predicted role of oxygen dependence. Our analysis of rapid and large shifts in GC content thus identifies new clades and novel contexts to further understand the ecological and evolutionary drivers of this important genomic trait.

A complex interplay of evolutionary forces continues to shape ancient co-occurring symbiont genomes

Many insects depend on ancient associations with intracellular bacteria for essential nutrition. The genomes of these bacteria are often highly reduced. Although drift is a major driver of symbiont evolution, other evolutionary forces continue to influence them. To understand how ongoing molecular evolution and gene loss shape symbiont genomes, we sequenced two of the most ancient symbionts known, Sulcia and Nasuia, from 20 Hawaiian Nesophrosyne leafhoppers. We leveraged the parallel divergence of Nesophrosyne lineages throughout Hawaii as a natural experimental framework. Sulcia and Nasuia experience ongoing-but divergent-gene loss, often in a convergent fashion. Although some genes are under relaxed selection, purifying and positive selection are also important drivers of genome evolution, particularly in maintaining certain nutritional and cellular functions. Our results further demonstrate that symbionts experience dramatically different evolutionary environments, as evidenced by the finding that Sulcia and Nasuia have one of the slowest and fastest rates of molecular evolution known.

The Di-Symbiotic Systems in the Aphids Sipha maydis and Periphyllus lyropictus Provide a Contrasting Picture of Recent Co-Obligate Nutritional Endosymbiosis in Aphids

Dependence on multiple nutritional bacterial symbionts forming a metabolic unit has repeatedly evolved in many insect species that feed on nutritionally unbalanced diets such as plant sap. This is the case for aphids of the subfamilies Lachninae and Chaitophorinae, which have evolved di-symbiotic systems in which the ancient obligate nutritional symbiont Buchnera aphidicola is metabolically complemented by an additional nutritional symbiont acquired more recently. Deciphering how different symbionts integrate both metabolically and anatomically in such systems is crucial to understanding how complex nutritional symbiotic systems function and evolve. In this study, we sequenced and analyzed the genomes of the symbionts B. aphidicola and Serratia symbiotica associated with the Chaitophorinae aphids Sipha maydis and Periphyllus lyropictus. Our results show that, in these two species, B. aphidicola and S. symbiotica complement each other metabolically (and their hosts) for the biosynthesis of essential amino acids and vitamins, but with distinct metabolic reactions supported by each symbiont depending on the host species. Furthermore, the S. symbiotica symbiont associated with S. maydis appears to be strictly compartmentalized into the specialized host cells housing symbionts in aphids, the bacteriocytes, whereas the S. symbiotica symbiont associated with P. lyropictus exhibits a highly invasive phenotype, presumably because it is capable of expressing a larger set of virulence factors, including a complete flagellum for bacterial motility. Such contrasting levels of metabolic and anatomical integration for two S. symbiotica symbionts that were recently acquired as nutritional co-obligate partners reflect distinct coevolutionary processes specific to each association.

Beyond specialization: re-examining routes of host influence on symbiont evolution

Our understanding of host influence on microbial evolution has focused on symbiont specialization and the genomic streamlining that often accompanies it. However, a vast diversity of symbiotic lineages facultatively interact with hosts or associate with multiple hosts. Yet, there are no clear expectations for how host association influences the niche of these symbionts or their evolution. Here, we discuss how weak or variable selection on microbial symbiotic associations, horizontal transmission, and low costs of adaptation to novel host habitats are predicted to promote the expansion or maintenance of microbial niches. This broad perspective will aid in developing better and more general predictions for evolution in microbial symbioses.

Unraveling coevolutionary dynamics using ecological genomics

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

Developmental Integration of Endosymbionts in Insects

In endosymbiosis, two independently existing entities are inextricably intertwined such that they behave as a single unit. For multicellular hosts, the endosymbiont must be integrated within the host developmental genetic network to maintain the relationship. Developmental integration requires innovations in cell type, gene function, gene regulation, and metabolism. These innovations are contingent upon the existing ecological interactions and may evolve mutual interdependence. Recent studies have taken significant steps toward characterizing the proximate mechanisms underlying interdependence. However, the study of developmental integration is only in its early stages of investigation. Here, we review the literature on mutualistic endosymbiosis to explore how unicellular endosymbionts developmentally integrate into their multicellular hosts with emphasis on insects as a model. Exploration of this process will help gain a more complete understanding of endosymbiosis. This will pave the way for a better understanding of the endosymbiotic theory of evolution in the future.

Prochlorococcus have low global mutation rate and small effective population size

Prochlorococcus are the most abundant free-living photosynthetic carbon-fixing organisms in the ocean. Prochlorococcus show small genome sizes, low genomic G+C content, reduced DNA repair gene pool and fast evolutionary rates, which are typical features of endosymbiotic bacteria. Nevertheless, their evolutionary mechanisms are believed to be different. Evolution of endosymbiotic bacteria is dominated by genetic drift owing to repeated population bottlenecks, whereas Prochlorococcus are postulated to have extremely large effective population sizes (N_e) and thus drift has rarely been considered. However, accurately extrapolating N_e requires measuring an unbiased global mutation rate through mutation accumulation, which is challenging for Prochlorococcus. Here, we managed this experiment over 1,065 days using Prochlorococcus marinus AS9601, sequenced genomes of 141 mutant lines and determined its mutation rate to be 3.50 × 10^-10 per site per generation. Extrapolating N_e additionally requires identifying population boundaries, which we defined using PopCOGenT and over 400 genomes related to AS9601. Accordingly, we calculated its N_e to be 1.68 × 10⁷, which is only reasonably greater than that of endosymbiotic bacteria but surprisingly smaller than that of many free-living bacteria extrapolated using the same approach. Our results therefore suggest that genetic drift is a key driver of Prochlorococcus evolution.

Tightly Constrained Genome Reduction and Relaxation of Purifying Selection during Secondary Plastid Endosymbiosis

Endosymbiosis, the establishment of a former free-living prokaryotic or eukaryotic cell as an organelle inside a host cell, can dramatically alter the genomic architecture of the endosymbiont. Plastids or chloroplasts, the light-harvesting organelle of photosynthetic eukaryotes, are excellent models to study this phenomenon because plastid origin has occurred multiple times in evolution. Here, we investigate the genomic signature of molecular processes acting through secondary plastid endosymbiosis-the origination of a new plastid from a free-living eukaryotic alga. We used phylogenetic comparative methods to study gene loss and changes in selective regimes on plastid genomes, focusing on green algae that have given rise to three independent lineages with secondary plastids (euglenophytes, chlorarachniophytes, and Lepidodinium). Our results show an overall increase in gene loss associated with secondary endosymbiosis, but this loss is tightly constrained by the retention of genes essential for plastid function. The data show that secondary plastids have experienced temporary relaxation of purifying selection during secondary endosymbiosis. However, this process is tightly constrained, with selection relaxed only relative to the background in primary plastids. Purifying selection remains strong in absolute terms even during the endosymbiosis events. Selection intensity rebounds to pre-endosymbiosis levels following endosymbiosis events, demonstrating the changes in selection efficiency during different origin phases of secondary plastids. Independent endosymbiosis events in the euglenophytes, chlorarachniophytes, and Lepidodinium differ in their degree of relaxation of selection, highlighting the different evolutionary contexts of these events. This study reveals the selection-drift interplay during secondary endosymbiosis and evolutionary parallels during organellogenesis.

Evolutionarily recent dual obligatory symbiosis among adelgids indicates a transition between fungus- and insect-associated lifestyles

Adelgids (Insecta: Hemiptera: Adelgidae) form a small group of insects but harbor a surprisingly diverse set of bacteriocyte-associated endosymbionts, which suggest multiple replacement and acquisition of symbionts over evolutionary time. Specific pairs of symbionts have been associated with adelgid lineages specialized on different secondary host conifers. Using a metagenomic approach, we investigated the symbiosis of the Adelges laricis/Adelges tardus species complex containing betaproteobacterial ("Candidatus Vallotia tarda") and gammaproteobacterial ("Candidatus Profftia tarda") symbionts. Genomic characteristics and metabolic pathway reconstructions revealed that Vallotia and Profftia are evolutionary young endosymbionts, which complement each other's role in essential amino acid production. Phylogenomic analyses and a high level of genomic synteny indicate an origin of the betaproteobacterial symbiont from endosymbionts of Rhizopus fungi. This evolutionary transition was accompanied with substantial loss of functions related to transcription regulation, secondary metabolite production, bacterial defense mechanisms, host infection, and manipulation. The transition from fungus to insect endosymbionts extends our current framework about evolutionary trajectories of host-associated microbes.

Reconstitution and Transmission of Gut Microbiomes and Their Genes between Generations

Microbiomes are transmitted between generations by a variety of different vertical and/or horizontal modes, including vegetative reproduction (vertical), via female germ cells (vertical), coprophagy and regurgitation (vertical and horizontal), physical contact starting at birth (vertical and horizontal), breast-feeding (vertical), and via the environment (horizontal). Analyses of vertical transmission can result in false negatives (failure to detect rare microbes) and false positives (strain variants). In humans, offspring receive most of their initial gut microbiota vertically from mothers during birth, via breast-feeding and close contact. Horizontal transmission is common in marine organisms and involves selectivity in determining which environmental microbes can colonize the organism's microbiome. The following arguments are put forth concerning accurate microbial transmission: First, the transmission may be of functions, not necessarily of species; second, horizontal transmission may be as accurate as vertical transmission; third, detection techniques may fail to detect rare microbes; lastly, microbiomes develop and reach maturity with their hosts. In spite of the great variation in means of transmission discussed in this paper, microbiomes and their functions are transferred from one generation of holobionts to the next with fidelity. This provides a strong basis for each holobiont to be considered a unique biological entity and a level of selection in evolution, largely maintaining the uniqueness of the entity and conserving the species from one generation to the next.

Genome-Resolved Metagenomic Analyses Reveal the Presence of a Putative Bacterial Endosymbiont in an Avian Nasal Mite (Rhinonyssidae; Mesostigmata)

Rhinonyssidae (Mesostigmata) is a family of nasal mites only found in birds. All species are hematophagous endoparasites, which may damage the nasal cavities of birds, and also could be potential reservoirs or vectors of other infections. However, the role of members of Rhinonyssidae as disease vectors in wild bird populations remains uninvestigated, with studies of the microbiomes of Rhinonyssidae being almost non-existent. In the nasal mite (Tinaminyssus melloi) from rock doves (Columba livia), a previous study found evidence of a highly abundant putatively endosymbiotic bacteria from Class Alphaproteobacteria. Here, we expanded the sample size of this species (two different hosts- ten nasal mites from two independent samples per host), incorporated contamination controls, and increased sequencing depth in shotgun sequencing and genome-resolved metagenomic analyses. Our goal was to increase the information regarding this mite species and its putative endosymbiont. We obtained a metagenome assembled genome (MAG) that was estimated to be 98.1% complete and containing only 0.9% possible contamination. Moreover, the MAG has characteristics typical of endosymbionts (namely, small genome size an AT bias). Overall, our results support the presence of a potential endosymbiont, which is the first described for avian nasal mites to date, and improve the overall understanding of the microbiota inhabiting these mites.

The Role of Gene Duplication in the Divergence of Enzyme Function: A Comparative Approach

Gene duplication is a crucial process involved in the appearance of new genes and functions. It is thought to have played a major role in the growth of enzyme families and the expansion of metabolism at the biosphere's dawn and in recent times. Here, we analyzed paralogous enzyme content within each of the seven enzymatic classes for a representative sample of prokaryotes by a comparative approach. We found a high ratio of paralogs for three enzymatic classes: oxidoreductases, isomerases, and translocases, and within each of them, most of the paralogs belong to only a few subclasses. Our results suggest an intricate scenario for the evolution of prokaryotic enzymes, involving different fates for duplicated enzymes fixed in the genome, where around 20-40% of prokaryotic enzymes have paralogs. Intracellular organisms have a lesser ratio of duplicated enzymes, whereas free-living enzymes show the highest ratios. We also found that phylogenetically close phyla and some unrelated but with the same lifestyle share similar genomic and biochemical traits, which ultimately support the idea that gene duplication is associated with environmental adaptation.

Endosymbiont Capture, a Repeated Process of Endosymbiont Transfer with Replacement in Trypanosomatids Angomonas spp

Trypanosomatids of the subfamily Strigomonadinae bear permanent intracellular bacterial symbionts acquired by the common ancestor of these flagellates. However, the cospeciation pattern inherent to such relationships was revealed to be broken upon the description of Angomonas ambiguus, which is sister to A. desouzai, but bears an endosymbiont genetically close to that of A. deanei. Based on phylogenetic inferences, it was proposed that the bacterium from A. deanei had been horizontally transferred to A. ambiguus. Here, we sequenced the bacterial genomes from two A. ambiguus isolates, including a new one from Papua New Guinea, and compared them with the published genome of the A. deanei endosymbiont, revealing differences below the interspecific level. Our phylogenetic analyses confirmed that the endosymbionts of A. ambiguus were obtained from A. deanei and, in addition, demonstrated that this occurred more than once. We propose that coinfection of the same blowfly host and the phylogenetic relatedness of the trypanosomatids facilitate such transitions, whereas the drastic difference in the occurrence of the two trypanosomatid species determines the observed direction of this process. This phenomenon is analogous to organelle (mitochondrion/plastid) capture described in multicellular organisms and, thereafter, we name it endosymbiont capture.

Relaxed selection underlies genome erosion in socially parasitic ant species

Inquiline ants are highly specialized and obligate social parasites that infiltrate and exploit colonies of closely related species. They have evolved many times convergently, are often evolutionarily young lineages, and are almost invariably rare. Focusing on the leaf-cutting ant genus Acromyrmex, we compared genomes of three inquiline social parasites with their free-living, closely-related hosts. The social parasite genomes show distinct signatures of erosion compared to the host lineages, as a consequence of relaxed selective constraints on traits associated with cooperative ant colony life and of inquilines having very small effective population sizes. We find parallel gene losses, particularly in olfactory receptors, consistent with inquiline species having highly reduced social behavioral repertoires. Many of the genomic changes that we uncover resemble those observed in the genomes of obligate non-social parasites and intracellular endosymbionts that branched off into highly specialized, host-dependent niches.

Patterns of Microbiome Variation Among Infrapopulations of Permanent Bloodsucking Parasites

While interspecific variation in microbiome composition can often be readily explained by factors such as host species identity, there is still limited knowledge of how microbiomes vary at scales lower than the species level (e.g., between individuals or populations). Here, we evaluated variation in microbiome composition of individual parasites among infrapopulations (i.e., populations of parasites of the same species living on a single host individual). To address this question, we used genome-resolved and shotgun metagenomic data of 17 infrapopulations (balanced design) of the permanent, bloodsucking seal louse Echinophthirius horridus sampled from individual Saimaa ringed seals Pusa hispida saimensis. Both genome-resolved and read-based metagenomic classification approaches consistently show that parasite infrapopulation identity is a significant factor that explains both qualitative and quantitative patterns of microbiome variation at the intraspecific level. This study contributes to the general understanding of the factors driving patterns of intraspecific variation in microbiome composition, especially of bloodsucking parasites, and has implications for understanding how well-known processes occurring at higher taxonomic levels, such as phylosymbiosis, might arise in these systems.

Finding Needles in Haystacks and Inferring Their Function: Challenges and Successes in Beneficial Symbiosis Research

Symbioses between hosts and beneficial microbes are key drivers of biological innovation and diversity. While a range of systems have emerged that provide foundational insights into how symbioses function and evolve, we still have a limited understanding of the vast diversity of organisms that engage in such interactions. Recent advances in molecular tools, theory, and interdisciplinary approaches now permit researchers to expand our knowledge and to press forward the frontiers of symbiosis research. As described in a recent issue of mSystems, Myers and colleagues (K. N. Myers, D. Conn, and A. M. V. Brown, mSystems, 6:e01048-20, 2021, https://doi.org/10.1128/mSystems.01048-20) conducted a genome skimming approach to understand the role of obligate beneficial symbionts in plant-parasitic dagger nematodes. Nematodes are extraordinarily abundant and key players in ecosystem function and health. However, they are difficult to harness in the lab. The approach used by Myers et al. ameliorates these challenges to illustrate a relatively complete picture of a poorly understood beneficial symbiosis.

Interactions of host-associated multispecies bacterial communities

The oral microbiome comprises microbial communities colonizing biotic (epithelia, mucosa) and abiotic (enamel) surfaces. Different communities are associated with health (eg, immune development, pathogen resistance) and disease (eg, tooth loss and periodontal disease). Like any other host-associated microbiome, colonization and persistence of both beneficial and dysbiotic oral microbiomes are dictated by successful utilization of available nutrients and defense against host and competitor assaults. This chapter will explore these general features of microbe-host interactions through the lens of symbiotic (mutualistic and antagonistic/pathogenic) associations with nonmammalian animals. Investigations in such systems across a broad taxonomic range have revealed conserved mechanisms and processes that underlie the complex associations among microbes and between microbes and hosts.

Mutational Pressure Drives Differential Genome Conservation in Two Bacterial Endosymbionts of Sap-Feeding Insects

Compared with free-living bacteria, endosymbionts of sap-feeding insects have tiny and rapidly evolving genomes. Increased genetic drift, high mutation rates, and relaxed selection associated with host control of key cellular functions all likely contribute to genome decay. Phylogenetic comparisons have revealed massive variation in endosymbiont evolutionary rate, but such methods make it difficult to partition the effects of mutation versus selection. For example, the ancestor of Auchenorrhynchan insects contained two obligate endosymbionts, Sulcia and a betaproteobacterium (BetaSymb; called Nasuia in leafhoppers) that exhibit divergent rates of sequence evolution and different propensities for loss and replacement in the ensuing ∼300 Ma. Here, we use the auchenorrhynchan leafhopper Macrosteles sp. nr. severini, which retains both of the ancestral endosymbionts, to test the hypothesis that differences in evolutionary rate are driven by differential mutagenesis. We used a high-fidelity technique known as duplex sequencing to measure and compare low-frequency variants in each endosymbiont. Our direct detection of de novode novo mutations reveals that the rapidly evolving endosymbiont (Nasuia) has a much higher frequency of single-nucleotide variants than the more stable endosymbiont (Sulcia) and a mutation spectrum that is potentially even more AT-biased than implied by the 83.1% AT content of its genome. We show that indels are common in both endosymbionts but differ substantially in length and distribution around repetitive regions. Our results suggest that differences in long-term rates of sequence evolution in Sulcia versus BetaSymb, and perhaps the contrasting degrees of stability of their relationships with the host, are driven by differences in mutagenesis.

Portiera Gets Wild: Genome Instability Provides Insights into the Evolution of Both Whiteflies and Their Endosymbionts

Whiteflies (Hemiptera: Sternorrhyncha: Aleyrodidae) are a superfamily of small phloem-feeding insects. They rely on their primary endosymbionts "Candidatus Portiera aleyrodidarum" to produce essential amino acids not present in their diet. Portiera has been codiverging with whiteflies since their origin and therefore reflects its host's evolutionary history. Like in most primary endosymbionts, the genome of Portiera stays stable across the Aleyrodidae superfamily after millions of years of codivergence. However, Portiera of the whitefly Bemisia tabaci has lost the ancestral genome order, reflecting a rare event in the endosymbiont evolution: the appearance of genome instability. To gain a better understanding of Portiera genome evolution, identify the time point in which genome instability appeared and contribute to the reconstruction of whitefly phylogeny, we developed a new phylogenetic framework. It targeted five Portiera genes and determined the presence of the DNA polymerase proofreading subunit (dnaQ) gene, previously associated with genome instability, and two alternative gene rearrangements. Our results indicated that Portiera gene sequences provide a robust tool for studying intergenera phylogenetic relationships in whiteflies. Using these new framework, we found that whitefly species from the Singhiella, Aleurolobus, and Bemisia genera form a monophyletic tribe, the Aleurolobini, and that their Portiera exhibit genome instability. This instability likely arose once in the common ancestor of the Aleurolobini tribe (at least 70 Ma), drawing a link between the appearance of genome instability in Portiera and the switch from multibacteriocyte to a single-bacteriocyte mode of inheritance in this tribe.

Horizontal transmission and recombination maintain forever young bacterial symbiont genomes

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inefficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosymbiont genomes are occasionally horizontally transmitted and are much larger than their terrestrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and substantial functional genetic variation.

A simple stochastic model describing genomic evolution over time of GC content in microbial symbionts

An organism's genomic base composition is usually summarized by its AT or GC content due to Chargaff's parity laws. Variation in prokaryotic GC content can be substantial between taxa but is generally small within microbial genomes. This variation has been found to correlate with both phylogeny and environmental factors. Since novel single-nucleotide polymorphisms (SNPs) within genomes are at least partially linked to the environment through natural selection, SNP GC content can be considered a compound measure of an organism's environmental influences, lifestyle, phylogeny as well as other more or less random processes. While there are several models describing genomic GC content few, if any, consider AT/GC mutation rates subjected to random perturbations. We present a mathematical model that describes how GC content in microbial genomes evolves over time as a function of the AT → GC and GC → AT mutation rates with Gaussian white noise disturbances. The model, which is suited specifically to non-recombining vertically transmitted prokaryotic symbionts, suggests that small differences in the AT/GC mutation rates can lead to profound differences in outcome due to the ensuing stochastic process. In other words, the model indicates that time to extinction could be a consequence of the mutation rate trajectory on which the symbiont embarked early on in its evolutionary history.

The Similarities between Human Mitochondria and Bacteria in the Context of Structure, Genome, and Base Excision Repair System

Mitochondria emerged from bacterial ancestors during endosymbiosis and are crucial for cellular processes such as energy production and homeostasis, stress responses, cell survival, and more. They are the site of aerobic respiration and adenosine triphosphate (ATP) production in eukaryotes. However, oxidative phosphorylation (OXPHOS) is also the source of reactive oxygen species (ROS), which are both important and dangerous for the cell. Human mitochondria contain mitochondrial DNA (mtDNA), and its integrity may be endangered by the action of ROS. Fortunately, human mitochondria have repair mechanisms that allow protecting mtDNA and repairing lesions that may contribute to the occurrence of mutations. Mutagenesis of the mitochondrial genome may manifest in the form of pathological states such as mitochondrial, neurodegenerative, and/or cardiovascular diseases, premature aging, and cancer. The review describes the mitochondrial structure, genome, and the main mitochondrial repair mechanism (base excision repair (BER)) of oxidative lesions in the context of common features between human mitochondria and bacteria. The authors present a holistic view of the similarities of mitochondria and bacteria to show that bacteria may be an interesting experimental model for studying mitochondrial diseases, especially those where the mechanism of DNA repair is impaired.

Ancestral niche separation and evolutionary rate differentiation between sister marine flavobacteria lineages

Marine flavobacteria are specialists for polysaccharide degradation. They dominate in habitats enriched with polysaccharides, but are also prevalent in pelagic environments where polysaccharides are less available. These niches are likely occupied by distinct lineages, but evolutionary processes underlying their niche differentiation remain elusive. Here, genomic analyses and physiological assays indicate that the sister flavobacteria lineages Leeuwenhoekiella and Nonlabens likely explore polysaccharide-rich macroalgae and polysaccharide-poor pelagic niches respectively. Phylogenomic analyses inferred that the niche separation likely occurred anciently and coincided with increased sequence evolutionary rate in Nonlabens compared with Leeuwenhoekiella. Further analyses ruled out the known mechanisms likely driving evolutionary rate acceleration, including reduced selection efficiency, decreased generation time and increased mutation rate. In particular, the mutation rates were determined using an unbiased experimental method, which measures the present-day populations and may not reflect ancestral populations. These data collectively lead to a new hypothesis that an ancestral and transient mutation rate increase resulted in evolutionary rate increase in Nonlabens. This hypothesis was supported by inferring that gains and losses of genes involved in SOS response, a mechanism known to drive transiently increased mutation rate, coincided with evolutionary rate acceleration. Our analyses highlight the evolutionary mechanisms underlying niche differentiation of flavobacteria lineages.

Chlamydia psittaci in fulmars on the Faroe Islands: a causative link to South American psittacines eight decades after a severe epidemic

A psittacosis epidemic linked to fulmar hunting occurred on the Faroe Islands in the 1930s. This study investigates a plausible explanation to the 20% human mortality in this outbreak. Phylogenetic analysis showed that Chlamydia psittaci isolated from fulmars were closely related to the highly virulent 6BC strains from psittacines and are compatible with an acquisition by fulmars of an ancestor of the 6BC clade in the 1930s. This supports the hypothesis that the outbreak on the Faroe Islands started after naïve fulmars acquired C. psittaci from infected dead parrots thrown overboard when shipped to Europe in the 1930s.

Symbiont replacements reset the co-evolutionary relationship between insects and their heritable bacteria

Auchenorrhynchan insects (Hemiptera) generally depend on two bacterial symbionts for nutrition. These bacteria experience extreme genome reduction and loss of essential cell functions that require direct host support, or the replacement of failing symbionts with more capable ones. However, it remains unclear how hosts adapt to integrate symbionts into their systems, particularly when they are replaced. Here, we comparatively investigated the evolution of host-support mechanisms in the glassy-winged sharpshooter, Homalodisca vitripennis (GWSS), and the aster leafhopper, Macrosteles quadrilineatus (ALF). ALF harbors the ancestral co-symbionts of the Auchenorrhyncha that have tiny genomes, Sulcia (190 kb) and Nasuia (112 kb). In GWSS, Sulcia retains an expanded genome (245 kb), but Nasuia was replaced by the more capable Baumannia (686 kb). To support their symbionts, GWSS and ALF have evolved novel mechanisms via horizontal gene transfer, gene duplication, and co-option of mitochondrial support genes. However, GWSS has fewer support systems targeting essential bacterial processes. In particular, although both hosts use ancestral mechanisms to support Sulcia, GWSS does not encode all of the same support genes required to sustain Sulcia-ALF or Nasuia. Moreover, GWSS support of Baumannia is far more limited and tailored to its expanded capabilities. Our results demonstrate how symbiont replacements shape host genomes and the co-evolutionary process.

Horizontally transmitted symbiont populations in deep-sea mussels are genetically isolated

Eukaryotes are habitats for bacterial organisms where the host colonization and dispersal among individual hosts have consequences for the bacterial ecology and evolution. Vertical symbiont transmission leads to geographic isolation of the microbial population and consequently to genetic isolation of microbiotas from individual hosts. In contrast, the extent of geographic and genetic isolation of horizontally transmitted microbiota is poorly characterized. Here we show that chemosynthetic symbionts of individual Bathymodiolus brooksi mussels constitute genetically isolated subpopulations. The reconstruction of core genome-wide strains from high-resolution metagenomes revealed distinct phylogenetic clades. Nucleotide diversity and strain composition vary along the mussel life span and individual hosts show a high degree of genetic isolation. Our results suggest that the uptake of environmental bacteria is a restricted process in B. brooksi, where self-infection of the gill tissue results in serial founder effects during symbiont evolution. We conclude that bacterial colonization dynamics over the host life cycle is thus an important determinant of population structure and genome evolution of horizontally transmitted symbionts.

Comparative genomics reveals a novel genetic organization of the sad cluster in the sulfonamide-degrader 'Candidatus Leucobacter sulfamidivorax' strain GP

Background

Microbial communities recurrently establish metabolic associations resulting in increased fitness and ability to perform complex tasks, such as xenobiotic degradation. In a previous study, we have described a sulfonamide-degrading consortium consisting of a novel low-abundant actinobacterium, named strain GP, and Achromobacter denitrificans PR1. However, we found that strain GP was unable to grow independently and could not be further purified.

Results

Previous studies suggested that strain GP might represent a new putative species within the Leucobacter genus (16S rRNA gene similarity < 97%). In this study, we found that average nucleotide identity (ANI) with other Leucobacter spp. ranged between 76.8 and 82.1%, further corroborating the affiliation of strain GP to a new provisional species. The average amino acid identity (AAI) and percentage of conserved genes (POCP) values were near the lower edge of the genus delimitation thresholds (65 and 55%, respectively). Phylogenetic analysis of core genes between strain GP and Leucobacter spp. corroborated these findings. Comparative genomic analysis indicates that strain GP may have lost genes related to tetrapyrrole biosynthesis and thiol transporters, both crucial for the correct assembly of cytochromes and aerobic growth. However, supplying exogenous heme and catalase was insufficient to abolish the dependent phenotype. The actinobacterium harbors at least two copies of a novel genetic element containing a sulfonamide monooxygenase (sadA) flanked by a single IS1380 family transposase. Additionally, two homologs of sadB (4-aminophenol monooxygenase) were identified in the metagenome-assembled draft genome of strain GP, but these were not located in the vicinity of sadA nor of mobile or integrative elements.

Conclusions

Comparative genomics of the genus Leucobacter suggested the absence of some genes encoding for important metabolic traits in strain GP. Nevertheless, although media and culture conditions were tailored to supply its potential metabolic needs, these conditions were insufficient to isolate the PR1-dependent actinobacterium further. This study gives important insights regarding strain GP metabolism; however, gene expression and functional studies are necessary to characterize and further isolate strain GP. Based on our data, we propose to classify strain GP in a provisional new species within the genus Leucobacter, ‘Candidatus Leucobacter sulfamidivorax‘.

Identification of novel Legionella genes required for endosymbiosis in Paramecium based on comparative genome analysis with Holospora spp

The relationship between Legionella and protist hosts has a huge impact when considering the infectious risk in humans because it facilitates the long-term replication and survival of Legionella in the environment. The ciliate Paramecium is considered to be a protist host for Legionella in natural environments, but the details of their endosymbiosis are largely unknown. In this study, we determined candidate Legionella pneumophila genes that are likely to be involved in the establishment of endosymbiosis in Paramecium caudatum by comparing the genomes of Legionella spp. and Holospora spp. that are obligate endosymbiotic bacteria in Paramecium spp. Among the candidate genes, each single deletion mutant for five genes (lpg0492, lpg0522, lpg0523, lpg2141 and lpg2398) failed to establish endosymbiosis in P. caudatum despite showing intracellular growth in human macrophages. The mutants exhibited no characteristic changes in terms of their morphology, multiplication rate or capacity for modulating the phagosomes in which they were contained, but their resistance to lysozyme decreased significantly. This study provides insights into novel factors required by L. pneumophila for endosymbiosis in P. caudatum, and suggests that endosymbiotic organisms within conspecific hosts may have shared genes related to effective endosymbiosis establishment.

Cladogenesis and Genomic Streamlining in Extracellular Endosymbionts of Tropical Stink Bugs

Phytophagous stink bugs are globally distributed and many harbor vertically inherited bacterial symbionts that are extracellular, yet little is known about how the symbiont’s genomes have evolved under this transmission strategy. Genome reduction is common in insect intracellular symbionts but limited genome sampling of the extracellular symbionts of distantly related stink bugs has precluded inferring patterns of extracellular symbiont genome evolution. To address this knowledge gap, we completely sequenced the genomes of the uncultivable bacterial symbionts of four neotropical stink bugs of the Edessa genus. Phylogenetic and comparative analyses indicated that the symbionts form a clade within the Pantoea genus and their genomes are highly reduced (∼0.8 Mb). Furthermore, genome synteny analysis and a jackknife approach for phylogenetic reconstruction, which corrected for long branch attraction artifacts, indicated that the Edessa symbionts were the result of a single symbiotic event that was distinct from the symbiosis event giving rise to Candidatus “Pantoea carbekii,” the extracellular symbiont of the invasive pentatomid stink bug, Halyomorpha halys. Metabolic functions inferred from the Edessa symbiont genomes suggests a shift in genomic composition characteristic of its lifestyle in that they retained many host-supportive functions while undergoing dramatic gene loss and establishing a stable relationship with their host insects. Given the undersampled nature of extracellular insect symbionts, this study is the first comparative analysis of these symbiont genomes from four distinct Edessa stink bug species. Finally, we propose the candidate name “Candidatus Pantoea edessiphila” for the species of these symbionts with strain designations according to their host species.

The First Purification of Functional Proteins from the Unculturable, Genome-Reduced, Bottlenecked α-Proteobacterium 'Candidatus Liberibacter solanacearum'

'Candidatus Liberibacter solanacearum' is an unculturable α-proteobacterium that is the causal agent of zebra chip disease of potato-a major problem in potato-growing areas, because it affects growth and yield. Developing effective treatments for 'Ca. L. solanacearum' has been hampered by the difficulty in functionally characterizing the proteins of this organism, largely because they are not easily expressed and purified in standard expression systems. 'Ca. L. solanacearum' has a reduced genome and its proteins are predicted to be prone to instability and aggregation. Among intracellular-dwelling bacteria, chaperone proteins are conserved and overexpressed to buffer against problems in protein folding. We mimicked this approach for expressing and purifying 'Ca. L. solanacearum' proteins in Escherichia coli by coexpressing them with chaperones. Neither of the representative 'Ca. L. solanacearum' enzymes, dihydrodipicolinate synthase (key in lysine biosynthesis) and pyruvate kinase (involved in glycolysis), were overexpressed in standard E. coli expression plasmids or strains. However, soluble dihydrodipicolinate synthase was successfully coexpressed with GroEL/GroES, while soluble pyruvate kinase was successfully coexpressed with either GroEL/GroES, dnaK/dnaJ/grpE, or a trigger factor. Both enzymes, believed to be key proteins for the organism, were purified by a combination of affinity chromatography and size-exclusion chromatography. Additionally, both 'Ca. L. solanacearum' enzymes are active and have the canonical tetrameric oligomeric structure in solution, consistent with other bacterial orthologs. This is the first study to successfully isolate and functionally characterize proteins from 'Ca. L. solanacearum'. Thus, we provide a general strategy for characterizing its proteins, enabling new research and drug discovery programs to study and manage the pathogenicity of the organism.

Selective advantages favour high genomic AT-contents in intracellular elements

Extrachromosomal genetic elements such as bacterial endosymbionts and plasmids generally exhibit AT-contents that are increased relative to their hosts’ DNA. The AT-bias of endosymbiotic genomes is commonly explained by neutral evolutionary processes such as a mutational bias towards increased A+T. Here we show experimentally that an increased AT-content of host-dependent elements can be selectively favoured. Manipulating the nucleotide composition of bacterial cells by introducing A+T-rich or G+C-rich plasmids, we demonstrate that cells containing GC-rich plasmids are less fit than cells containing AT-rich plasmids. Moreover, the cost of GC-rich elements could be compensated by providing precursors of G+C, but not of A+T, thus linking the observed fitness effects to the cytoplasmic availability of nucleotides. Accordingly, introducing AT-rich and GC-rich plasmids into other bacterial species with different genomic GC-contents revealed that the costs of G+C-rich plasmids decreased with an increasing GC-content of their host’s genomic DNA. Taken together, our work identifies selection as a strong evolutionary force that drives the genomes of intracellular genetic elements toward higher A+T contents.

Genomes of endosymbiotic bacteria are commonly more AT-rich than the ones of their free-living relatives. Interestingly, genomes of other intracellular elements like plasmids or bacteriophages also tend to be richer in AT than the genomes of their hosts. The AT-bias of endosymbiotic genomes is commonly explained by neutral evolutionary processes. However, since A+T nucleotides are both more abundant and energetically less expensive than G+C nucleotides, an alternative explanation is that selective advantages drive the nucleotide composition of intracellular elements. Here we provide strong experimental evidence that intracellular elements, whose genome is more AT-rich than the genome of the host, are selectively favoured on the host level. Thus, our results emphasize the importance of selection for shaping the DNA base composition of extrachromosomal genetic elements.

Evolutionary dynamics of the chromatophore genome in three photosynthetic Paulinella species

The thecate amoeba Paulinella is a valuable model for understanding plastid organellogenesis because this lineage has independently gained plastids (termed chromatophores) of alpha-cyanobacterial provenance. Plastid primary endosymbiosis in Paulinella occurred relatively recently (90–140 million years ago, Mya), whereas the origin of the canonical Archaeplastida plastid occurred >1,500 Mya. Therefore, these two events provide independent perspectives on plastid formation on vastly different timescales. Here we generated the complete chromatophore genome sequence from P. longichromatophora (979,356 bp, GC-content = 38.8%, 915 predicted genes) and P. micropora NZ27 (977,190 bp, GC-content = 39.9%, 911 predicted genes) and compared these data to that from existing chromatophore genomes. Our analysis suggests that when a basal split occurred among photosynthetic Paulinella species ca. 60 Mya, only 35% of the ancestral orthologous gene families from the cyanobacterial endosymbiont remained in chromatophore DNA. Following major gene losses during the early stages of endosymbiosis, this process slowed down significantly, resulting in a conserved gene content across extant taxa. Chromatophore genes faced relaxed selection when compared to homologs in free-living alpha-cyanobacteria, likely reflecting the homogeneous intracellular environment of the Paulinella host. Comparison of nucleotide substitution and insertion/deletion events among different P. micropora strains demonstrates that increases in AT-content and genome reduction are ongoing and dynamic processes in chromatophore evolution.

Fueled by methane: deep-sea sponges from asphalt seeps gain their nutrition from methane-oxidizing symbionts

Sponges host a remarkable diversity of microbial symbionts, however, the benefit their microbes provide is rarely understood. Here, we describe two new sponge species from deep-sea asphalt seeps and show that they live in a nutritional symbiosis with methane-oxidizing (MOX) bacteria. Metagenomics and imaging analyses revealed unusually high amounts of MOX symbionts in hosts from a group previously assumed to have low microbial abundances. These symbionts belonged to the Marine Methylotrophic Group 2 clade. They are host-specific and likely vertically transmitted, based on their presence in sponge embryos and streamlined genomes, which lacked genes typical of related free-living MOX. Moreover, genes known to play a role in host–symbiont interactions, such as those that encode eukaryote-like proteins, were abundant and expressed. Methane assimilation by the symbionts was one of the most highly expressed metabolic pathways in the sponges. Molecular and stable carbon isotope patterns of lipids confirmed that methane-derived carbon was incorporated into the hosts. Our results revealed that two species of sponges, although distantly related, independently established highly specific, nutritional symbioses with two closely related methanotrophs. This convergence in symbiont acquisition underscores the strong selective advantage for these sponges in harboring MOX bacteria in the food-limited deep sea.

Trophic Specialization Results in Genomic Reduction in Free-Living Marine Idiomarina Bacteria

The streamlining hypothesis is usually used to explain the genomic reduction events in free-living bacteria like SAR11. However, we find that the genomic reduction phenomenon in the bacterial genus Idiomarina is different from that in SAR11. Therefore, we propose a new hypothesis to explain genomic reduction in this genus based on trophic specialization that could result in genomic reduction, which would be not uncommon in nature. Not only can the trophic specialization hypothesis explain the genomic reduction in the genus Idiomarina, but it also sheds new light on our understanding of the genomic reduction processes in other free-living bacterial lineages.

The streamlining hypothesis is generally used to explain the genomic reduction events related to the small genome size of free-living bacteria like marine bacteria SAR11. However, our current understanding of the correlation between bacterial genome size and environmental adaptation relies on too few species. It is still unclear whether there are other paths leading to genomic reduction in free-living bacteria. The genome size of marine free-living bacteria of the genus Idiomarina belonging to the order Alteromonadales (Gammaproteobacteria) is much smaller than the size of related genomes from bacteria in the same order. Comparative genomic and physiological analyses showed that the genomic reduction pattern in this genus is different from that of the classical SAR11 lineage. Genomic reduction reconstruction and substrate utilization profile showed that Idiomarina spp. lost a large number of genes related to carbohydrate utilization, and instead they specialized on using proteinaceous resources. Here we propose a new hypothesis to explain genomic reduction in this genus; we propose that trophic specialization increasing the metabolic efficiency for using one kind of substrate but reducing the substrate utilization spectrum could result in bacterial genomic reduction, which would be not uncommon in nature. This hypothesis was further tested in another free-living genus, Kangiella, which also shows dramatic genomic reduction. These findings highlight that trophic specialization is potentially an important path leading to genomic reduction in some marine free-living bacteria, which is distinct from the classical lineages like SAR11.

Evolution of host support for two ancient bacterial symbionts with differentially degraded genomes in a leafhopper host

Plant sap-feeding insects (Hemiptera) rely on bacterial symbionts for nutrition absent in their diets. These bacteria experience extreme genome reduction and require genetic resources from their hosts, particularly for basic cellular processes other than nutrition synthesis. The host-derived mechanisms that complete these processes have remained poorly understood. It is also unclear how hosts meet the distinct needs of multiple bacterial partners with differentially degraded genomes. To address these questions, we investigated the cell-specific gene-expression patterns in the symbiotic organs of the aster leafhopper (ALF), Macrosteles quadrilineatus (Cicadellidae). ALF harbors two intracellular symbionts that have two of the smallest known bacterial genomes: Nasuia (112 kb) and Sulcia (190 kb). Symbionts are segregated into distinct host cell types (bacteriocytes) and vary widely in their basic cellular capabilities. ALF differentially expresses thousands of genes between the bacteriocyte types to meet the functional needs of each symbiont, including the provisioning of metabolites and support of cellular processes. For example, the host highly expresses genes in the bacteriocytes that likely complement gene losses in nucleic acid synthesis, DNA repair mechanisms, transcription, and translation. Such genes are required to function in the bacterial cytosol. Many host genes comprising these support mechanisms are derived from the evolution of novel functional traits via horizontally transferred genes, reassigned mitochondrial support genes, and gene duplications with bacteriocyte-specific expression. Comparison across other hemipteran lineages reveals that hosts generally support the incomplete symbiont cellular processes, but the origins of these support mechanisms are generally specific to the host-symbiont system.