Genome Evolution of the Obligate Endosymbiont Buchnera aphidicola

Nature should be the model for microbial sciences

Until recently, microbiologists have relied on cultures to understand the microbial world. As a result, model organisms have been the focus of research into understanding Bacteria and Archaea at a molecular level. Diversity surveys and metagenomic sequencing have revealed that these model species are often present in low abundance in the environment; instead, there are microbial taxa that are cosmopolitan in nature. Due to the numerical dominance of these microorganisms and the size of their habitats, these lineages comprise mind-boggling population sizes upward of 1028 cells on the planet. Many of these dominant groups have cultured representatives and have been shown to be involved in mediating key processes in nature. Given their importance and the increasing need to understand changes due to climate change, we propose that members of Nitrosophaerota (Nitrosopumilus maritimus), SAR11 (Pelagibacter ubique), Hadesarchaeia, Bathyarchaeia, and others become models in the future. Abundance should not be the only measure of a good model system; there are other organisms that are well suited to advance our understanding of ecology and evolution. For example, the most well-studied symbiotic bacteria, like Buchnera, Aliivibrio, and Rhizobium, should be models for understanding host-associations. Also, there are organisms that hold new insights into major transitions in the evolution of life on the planet like the Asgard Archaea (Heimdallarchaeia). Innovations in a variety of in situ techniques have enabled us to circumvent culturing when studying everything from genetics to physiology. Our deepest understanding of microbiology and its impact on the planet will come from studying these microbes in nature. Laboratory-based studies must be grounded in nature, not the other way around.

High-throughput transposon mutagenesis in the family Enterobacteriaceae reveals core essential genes and rapid turnover of essentiality

The Enterobacteriaceae are a scientifically and medically important clade of bacteria, containing the model organism Escherichia coli, as well as major human pathogens including Salmonella enterica and Klebsiella pneumoniae. Essential gene sets have been determined for several members of the Enterobacteriaceae, with the Keio E. coli single-gene deletion library often regarded as a gold standard. However, it remains unclear how gene essentiality varies between related strains and species. To investigate this, we have assembled a collection of 13 sequenced high-density transposon mutant libraries from five genera within the Enterobacteriaceae. We first assess several gene essentiality prediction approaches, investigate the effects of transposon density on essentiality prediction, and identify biases in transposon insertion sequencing data. Based on these investigations, we develop a new classifier for gene essentiality. Using this new classifier, we define a core essential genome in the Enterobacteriaceae of 201 universally essential genes. Despite the presence of a large cohort of variably essential genes, we find an absence of evidence for genus-specific essential genes. A clear example of this sporadic essentiality is given by the set of genes regulating the σE extracytoplasmic stress response, which appears to have independently acquired essentiality multiple times in the Enterobacteriaceae. Finally, we compare our essential gene sets to the natural experiment of gene loss in obligate insect endosymbionts that have emerged from within the Enterobacteriaceae. This isolates a remarkably small set of genes absolutely required for survival and identifies several instances of essential stress responses masked by redundancy in free-living bacteria.IMPORTANCEThe essential genome, that is the set of genes absolutely required to sustain life, is a core concept in genetics. Essential genes in bacteria serve as drug targets, put constraints on the engineering of biological chassis for technological or industrial purposes, and are key to constructing synthetic life. Despite decades of study, relatively little is known about how gene essentiality varies across related bacteria. In this study, we have collected gene essentiality data for 13 bacteria related to the model organism Escherichia coli, including several human pathogens, and investigated the conservation of essentiality. We find that approximately a third of the genes essential in any particular strain are non-essential in another related strain. Surprisingly, we do not find evidence for essential genes unique to specific genera; rather it appears a substantial fraction of the essential genome rapidly gains or loses essentiality during evolution. This suggests that essentiality is not an immutable characteristic but depends crucially on the genomic context. We illustrate this through a comparison of our essential genes in free-living bacteria to genes conserved in 34 insect endosymbionts with naturally reduced genomes, finding several cases where genes generally regarded as being important for specific stress responses appear to have become essential in endosymbionts due to a loss of functional redundancy in the genome.

'Drifting' Buchnera genomes track the microevolutionary trajectories of their aphid hosts

Evolution of Buchnera-aphid host symbioses is often studied among species at macroevolutionary scales. Investigations within species offer a different perspective about how eco-evolutionary processes shape patterns of genetic variation at microevolutionary scales. Our study leverages new and publicly available whole-genome sequencing data to study Buchnera-aphid host evolution in Myzus persicae, the peach potato aphid, a globally invasive and polyphagous pest. Across 43 different asexual, clonally reproducing isofemale strains, we examined patterns of genomic covariation between Buchnera and their aphid host and considered the distribution of mutations in protein-coding regions of the Buchnera genome. We found Buchnera polymorphisms within aphid strains, suggesting the presence of genetically different Buchnera strains within the same clonal lineage. Genetic distance between pairs of Buchnera samples was positively correlated to genetic distance between their aphid hosts, indicating shared evolutionary histories. However, there was no segregation of genetic variation for both M. persicae and Buchnera with plant host (Brassicaceae and non-tobacco Solanaceae) and no associations between genetic and geographic distance at global or regional spatial scales. Abundance patterns of non-synonymous mutations were similar to synonymous mutations in the Buchnera genome, and both mutation classes had similar site frequency spectra. We hypothesize that a predominance of neutral processes results in the Buchnera of M. persicae to simply 'drift' with the evolutionary trajectory of their aphid hosts. Our study presents a unique microevolutionary characterization of Buchnera-aphid host genomic covariation across multiple aphid clones. This provides a new perspective on the eco-evolutionary processes generating and maintaining polymorphisms in a major pest aphid species and its obligate primary endosymbiont.

Repeated horizontal acquisition of lagriamide-producing symbionts in Lagriinae beetles

Microbial symbionts associate with multicellular organisms on a continuum from facultative associations to mutual codependency. In some of the oldest intracellular symbioses there is exclusive vertical symbiont transmission, and co-diversification of symbiotic partners over millions of years. Such symbionts often undergo genome reduction due to low effective population sizes, frequent population bottlenecks, and reduced purifying selection. Here, we describe multiple independent acquisition events of closely related defensive symbionts followed by genome erosion in a group of Lagriinae beetles. Previous work in Lagria villosa revealed the dominant genome-eroded symbiont of the genus Burkholderia produces the antifungal compound lagriamide and protects the beetle's eggs and larvae from antagonistic fungi. Here, we use metagenomics to assemble 11 additional genomes of lagriamide-producing symbionts from seven different host species within Lagriinae from five countries, to unravel the evolutionary history of this symbiotic relationship. In each host species, we detected one dominant genome-eroded Burkholderia symbiont encoding the lagriamide biosynthetic gene cluster (BGC). Surprisingly, however, we did not find evidence for host-symbiont co-diversification, or for a monophyly of the lagriamide-producing symbionts. Instead, our analyses support at least four independent acquisition events of lagriamide-encoding symbionts and subsequent genome erosion in each of these lineages. By contrast, a clade of plant-associated relatives retained large genomes but secondarily lost the lagriamide BGC. In conclusion, our results reveal a dynamic evolutionary history with multiple independent symbiont acquisitions characterized by high degree of specificity. They highlight the importance of the specialized metabolite lagriamide for the establishment and maintenance of this defensive symbiosis.

Twenty-nine newly sequenced genomes and a comprehensive genome dataset for the insect endosymbiont Buchnera

Most phloem-feeding insects face nutritional deficiency and rely on their intracellular symbionts to provide nutrients, and most of endosymbiont genomes have undergone reduction. However, the study of genome reduction processes of endosymbionts has been constrained by the limited availability of genome data from different insect lineages. The obligate relationship between aphids and Buchnera aphidicola (hereafter Buchnera) makes them a classic model for studying insect-endosymbiont interaction. Here, we report 29 newly sequenced Buchnera genomes from 11 aphid subfamilies, and a comprehensive dataset based on 90 Buchnera genomes from 14 aphid subfamilies. The dataset shows a significant genomic difference of Buchnera among different aphid lineages. The dataset exhibits a more balanced distribution of Buchnera (from 14 aphid subfamilies) genome sizes, ranging from 400 kb to 600 kb, which can illustrate the genome reduction process of Buchnera. The new genome data provide valuable insights into the microevolutionary processes leading to genomic reduction of insect endosymbionts.

Comparative genomics of the primary endosymbiont Buchnera aphidicola in aphid hosts and their coevolutionary relationships

BACKGROUND

Coevolution between modern aphids and their primary obligate, bacterial endosymbiont, Buchnera aphidicola, has been previously reported at different classification levels based on molecular phylogenetic analyses. However, the Buchnera genome remains poorly understood within the Rhus gall aphids.

RESULTS

We assembled the complete genome of the endosymbiont Buchnera in 16 aphid samples, representing 13 species in all six genera of Rhus gall aphids by shotgun genome skimming method. We compared the newly assembled genomes with those from GenBank to comprehensively investigate patterns of coevolution between the bacteria Buchnera and their aphid hosts. Buchnera genomes were mostly collinear, and the pan-genome contained 684 genes, in which the core genome contained 256 genes with some lineages having large numbers of tandem gene duplications. There has been substantial gene-loss in each Buchnera lineage. We also reconstructed the phylogeny for Buchnera and their host aphids, respectively, using 72 complete genomes of Buchnera, along with the complete mitochondrial genomes and three nuclear genes of 31 corresponding host aphid accessions. The cophylogenetic test demonstrated significant coevolution between these two partner groups at individual, species, generic, and tribal levels.

CONCLUSIONS

Buchnera exhibits very high levels of genomic sequence divergence but relative stability in gene order. The relationship between the symbionts Buchnera and its aphid hosts shows a significant coevolutionary pattern and supports complexity of the obligate symbiotic relationship.

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.

Paleocene origin of a streamlined digestive symbiosis in leaf beetles

Timing the acquisition of a beneficial microbe relative to the evolutionary history of its host can shed light on the adaptive impact of a partnership. Here, we investigated the onset and molecular evolution of an obligate symbiosis between Cassidinae leaf beetles and Candidatus Stammera capleta, a γ-proteobacterium. Residing extracellularly within foregut symbiotic organs, Stammera upgrades the digestive physiology of its host by supplementing plant cell wall-degrading enzymes. We observe that Stammera is a shared symbiont across tortoise and hispine beetles that collectively comprise the Cassidinae subfamily, despite differences in their folivorous habits. In contrast to its transcriptional profile during vertical transmission, Stammera elevates the expression of genes encoding digestive enzymes while in the foregut symbiotic organs, matching the nutritional requirements of its host. Despite the widespread distribution of Stammera across Cassidinae beetles, symbiont acquisition during the Paleocene (∼62 mya) did not coincide with the origin of the subfamily. Early diverging lineages lack the symbiont and the specialized organs that house it. Reconstructing the ancestral state of host-beneficial factors revealed that Stammera encoded three digestive enzymes at the onset of symbiosis, including polygalacturonase-a pectinase that is universally shared. Although non-symbiotic cassidines encode polygalacturonase endogenously, their repertoire of plant cell wall-degrading enzymes is more limited compared with symbiotic beetles supplemented with digestive enzymes from Stammera. Highlighting the potential impact of a symbiotic condition and an upgraded metabolic potential, Stammera-harboring beetles exploit a greater variety of plants and are more speciose compared with non-symbiotic members of the Cassidinae.

Host-specific co-evolution likely driven by diet in Buchnera aphidicola

BACKGROUND

Russian wheat aphid (Diuraphis noxia Kurd.) is a severe pest to wheat, and even though resistance varieties are available to curb this pest, they are becoming obsolete with the development of new virulent aphid populations. Unlike many other aphids, D noxia only harbours a single endosymbiont, Buchnera aphidicola. Considering the importance of Buchnera, this study aimed to elucidate commonalities and dissimilarities between various hosts, to better understand its distinctiveness within its symbiotic relationship with D. noxia. To do so, the genome of the D. noxia's Buchnera was assembled and compared to those of other aphid species that feed on diverse host species.

RESULTS

The overall importance of several features such as gene length and percentage GC content was found to be critical for the maintenance of Buchnera genes when compared to their closest free-living relative, Escherichia coli. Buchnera protein coding genes were found to have percentage GC contents that tended towards a mean of ~ 26% which had strong correlation to their identity to their E. coli homologs. Several SNPs were identified between different aphid populations and multiple isolates of Buchnera were confirmed in single aphids.

CONCLUSIONS

Establishing the strong correlation of percentage GC content of protein coding genes and gene identity will allow for identifying which genes will be lost in the continually shrinking Buchnera genome. This is also the first report of a parthenogenically reproducing aphid that hosts multiple Buchnera strains in a single aphid, raising questions regarding the benefits of maintaining multiple strains. We also found preliminary evidence for post-transcriptional regulation of Buchnera genes in the form of polyadenylation.

Genome evolution of Buchnera aphidicola (Gammaproteobacteria): Insights into strand compositional asymmetry, codon usage bias, and phylogenetic implications

Taxa of Buchnera aphidicola (hereafter "Buchnera") are mutualistic intracellular symbionts of aphids, known for their remarkable biological traits such as genome reduction, strand compositional asymmetry, and symbiont-host coevolution. With the growing availability of genomic data, we performed a comprehensive analysis of 103 genomes of Buchnera strains from 12 host subfamilies, focusing on the genomic characterizations, codon usage patterns, and phylogenetic implications. Our findings revealed consistent features among all genomes, including small genome sizes, low GC contents, and gene losses. We also identified strong strand compositional asymmetries in all strains at the genome level. Further investigation suggested that mutation pressure may have played a crucial role in shaping codon usage of Buchnera. Moreover, the genomic asymmetries were reflected in asymmetric codon usage preferences within chromosomal genes. Notably, the levels of these asymmetries were varied among strains and were significantly influenced by the degrees of genome shrinkages. Lastly, our phylogenetic analyses presented an alternative topology of Aphididae, based on the Buchnera symbionts, providing robust confirmation of the paraphylies of Eriosomatinae, and Macrosiphini. Our objectives are to further understand the strand compositional asymmetry and codon usage bias of Buchnera taxa, and provide new perspectives for phylogenetic studies of Aphididae.

The nutritional dimension of facultative bacterial symbiosis in aphids: Current status and methodological considerations for future research

Aphids are valuable models for studying the functional diversity of bacterial symbiosis in insects. In addition to their ancestral obligate nutritional symbiont Buchnera aphidicola, these insects can host a myriad of so-called facultative symbionts. The diversity of these heritable bacterial associates is now well known, and some of the ecologically important traits associated with them have been well documented. Some twenty years ago, it was suggested that facultative symbionts could play an important role in aphid nutrition, notably by improving feeding performance on specific host plants, thus influencing the adaptation of these insects to host plants. However, the underlying mechanisms have never been elucidated, and the nutritional role that facultative symbionts might perform in aphids remains enigmatic. In this opinion piece, I put forward a series of arguments in support of the hypothesis that facultative symbionts play a central role in aphid nutrition and emphasize methodological considerations for testing this hypothesis in future work. In particular, I hypothesize that the metabolic capacities of B. aphidicola alone may not always be able to counterbalance the nutritional deficiencies of phloem sap. The association with one or several facultative symbionts with extensive metabolic capabilities would then be necessary to buffer the insect from host plant-derived nutrient deficiencies, thus enabling it to gain access to certain host plants.

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

SOCfinder: a genomic tool for identifying social genes in bacteria

Bacteria cooperate by working collaboratively to defend their colonies, share nutrients, and resist antibiotics. Nevertheless, our understanding of these remarkable behaviours primarily comes from studying a few well-characterized species. Consequently, there is a significant gap in our understanding of microbial social traits, particularly in natural environments. To address this gap, we can use bioinformatic tools to identify genes that control cooperative or otherwise social traits. Existing tools address this challenge through two approaches. One approach is to identify genes that encode extracellular proteins, which can provide benefits to neighbouring cells. An alternative approach is to predict gene function using annotation tools. However, these tools have several limitations. Not all extracellular proteins are cooperative, and not all cooperative behaviours are controlled by extracellular proteins. Furthermore, existing functional annotation methods frequently miss known cooperative genes. We introduce SOCfinder as a new tool to find bacterial genes that control cooperative or otherwise social traits. SOCfinder combines information from several methods, considering if a gene is likely to [1] code for an extracellular protein [2], have a cooperative functional annotation, or [3] be part of the biosynthesis of a cooperative secondary metabolite. We use data on two extensively-studied species (P. aeruginosa and B. subtilis) to show that SOCfinder is better at finding known cooperative genes than existing tools. We also use theory from population genetics to identify a signature of kin selection in SOCfinder cooperative genes, which is lacking in genes identified by existing tools. SOCfinder opens up a number of exciting directions for future research, and is available to download from https://github.com/lauriebelch/SOCfinder.

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.

Comparative Genomics Reveals Three Genetic Groups of the Whitefly Obligate Endosymbiont Candidatus Portiera aleyrodidarum

Maternally inherited obligate endosymbionts codiverge with their invertebrate hosts and reflect their host's evolutionary history. Whiteflies (Hemiptera: Aleyrodidae) harbor one obligate endosymbiont, Candidatus Portiera aleyrodidarum (hereafter Portiera). Portiera was anciently acquired by whitefly and has been coevolving with its host ever since. Uncovering the divergence of endosymbionts provides a fundamental basis for inspecting the coevolutionary processes between the bacteria and their hosts. To illustrate the divergence of Portiera lineages across different whitefly species, we sequenced the Portiera genome from Aleyrodes shizuokensis and conducted a comparative analysis on the basic features and gene evolution with bacterial genomes from five whitefly genera, namely Aleurodicus, Aleyrodes, Bemisia, Pealius, and Trialeurodes. The results indicated that Portiera from Bemisia possessed significantly larger genomes, fewer coding sequences (CDSs), and a lower coding density. Their gene arrangement differed notably from those of other genera. The phylogeny of the nine Portiera lineages resembled that of their hosts. Moreover, the lineages were classified into three distinct genetic groups based on the genetic distance, one from Aleurodicus (Aleurodicinae), one from Bemisia (Aleyrodinae), and another from Aleyrodes, Pealius, and Trialeurrodes (Aleyrodinae). Synonymous and nonsynonymous rate analyses, parity rule 2 plot analyses, neutrality plot analyses, and effective number of codons analyses supported the distinction of the three genetic groups. Our results indicated that Portiera from distant hosts exhibit distinct genomic contents, implying codivergence between hosts and their endosymbionts. This work will enhance our understanding of coevolution between hosts and their endosymbionts.

Phylosymbiosis: The Eco-Evolutionary Pattern of Insect-Symbiont Interactions

Insects harbor diverse assemblages of bacterial and fungal symbionts, which play crucial roles in host life history. Insects and their various symbionts represent a good model for studying host-microbe interactions. Phylosymbiosis is used to describe an eco-evolutionary pattern, providing a new cross-system trend in the research of host-associated microbiota. The phylosymbiosis pattern is characterized by a significant positive correlation between the host phylogeny and microbial community dissimilarities. Although host-symbiont interactions have been demonstrated in many insect groups, our knowledge of the prevalence and mechanisms of phylosymbiosis in insects is still limited. Here, we provide an order-by-order summary of the phylosymbiosis patterns in insects, including Blattodea, Coleoptera, Diptera, Hemiptera, Hymenoptera, and Lepidoptera. Then, we highlight the potential contributions of stochastic effects, evolutionary processes, and ecological filtering in shaping phylosymbiotic microbiota. Phylosymbiosis in insects can arise from a combination of stochastic and deterministic mechanisms, such as the dispersal limitations of microbes, codiversification between symbionts and hosts, and the filtering of phylogenetically conserved host traits (incl., host immune system, diet, and physiological characteristics).

Evolution of an Alternative Genetic Code in the Providencia Symbiont of the Hematophagous Leech Haementeria acuecueyetzin

Strict blood-feeding animals are confronted with a strong B-vitamin deficiency. Blood-feeding leeches from the Glossiphoniidae family, similarly to hematophagous insects, have evolved specialized organs called bacteriomes to harbor symbiotic bacteria. Leeches of the Haementeria genus have two pairs of globular bacteriomes attached to the esophagus which house intracellular "Candidatus Providencia siddallii" bacteria. Previous work analyzing a draft genome of the Providencia symbiont of the Mexican leech Haementeria officinalis showed that, in this species, the bacteria hold a reduced genome capable of synthesizing B vitamins. In this work, we aimed to expand our knowledge on the diversity and evolution of Providencia symbionts of Haementeria. For this purpose, we sequenced the symbiont genomes of three selected leech species. We found that all genomes are highly syntenic and have kept a stable genetic repertoire, mirroring ancient insect endosymbionts. Additionally, we found B-vitamin pathways to be conserved among these symbionts, pointing to a conserved symbiotic role. Lastly and most notably, we found that the symbiont of H. acuecueyetzin has evolved an alternative genetic code, affecting a portion of its proteome and showing evidence of a lineage-specific and likely intermediate stage of genetic code reassignment.

Phylloxera and aphids show distinct features of genome evolution despite similar reproductive modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared to the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with one sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared to their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Endogenous Plasmids and Chromosomal Genome Reduction in the Cardinium Endosymbiont of Dermatophagoides farinae

Cardinium bacteria are well known as endosymbionts that infect a wide range of arthropods and can manipulate host reproduction to promote their vertical transmission. As intracellular bacteria, Cardinium species undergo dramatic genome evolution, especially their chromosomal genome reduction. Although Cardinium plasmids have been reported to harbor important genes, the role of these plasmids in the genome evolution is yet to be fully understood. In this study, 2 genomes of Cardinium endosymbiont bacteria in astigmatic mites were de novo assembled, including the complete circular chromosomal genome of Cardinium sp. DF that was constructed in high quality using high-coverage long-read sequencing data. Intriguingly, 2 circular plasmids were assembled in Cardinium sp. DF and were identified to be endogenous for over 10 homologous genes shared with the chromosomal genome. Comparative genomics analysis illustrated an outline of the genome evolution of Cardinium bacteria, and the in-depth analysis of Cardinium sp. DF shed light on the multiple roles of endogenous plasmids in the molecular process of the chromosomal genome reduction. The endogenous plasmids of Cardinium sp. DF not only harbor massive homologous sequences that enable homologous recombination with the chromosome, but also can provide necessary functional proteins when the coding genes decayed in the chromosomal genome. IMPORTANCE As bacterial endosymbionts, Cardinium typically undergoes genome reduction, but the molecular process is still unclear, such as how plasmids get involved in chromosome reduction. Here, we de novo assembled 2 genomes of Cardinium in astigmatic mites, especially the chromosome of Cardinium sp. DF was assembled in a complete circular DNA using high-coverage long-read sequencing data. In the genome assembly of Cardinium sp. DF, 2 circular endogenous plasmids were identified to share at least 10 homologous genes with the chromosomal genome. In the comparative analysis, we identified a range of genes decayed in the chromosomal genome of Cardinium sp. DF but preserved in the 2 plasmids. Taken together with in-depth analyses, our results unveil that the endogenous plasmids harbor homologous sequences of chromosomal genome and can provide a structural basis of homologous recombination. Overall, this study reveals that endogenous plasmids participate in the ongoing chromosomal genome reduction of Cardinium sp. DF.

Loss of Flagella-Related Genes Enables a Nonflagellated, Fungal-Predating Bacterium To Strengthen the Synthesis of an Antifungal Weapon

Loss of flagellar genes causes a nonmotile phenotype. The genus Lysobacter consists of numerous environmentally ubiquitous, nonflagellated bacteria, including Lysobacter enzymogenes, an antifungal bacterium that is beneficial to plants. L. enzymogenes still has many flagellar genes on its genome, although this bacterium does not engage in flagella-driven motility. Here, we report that loss of certain flagellar genes allows L. enzymogenes to strengthen its evolutionarily gained capacity in fungal killing. To clarify why this bacterium loses flagellar genes during the evolutionary process, we cloned several representative flagellar genes from Xanthomonas oryzae, a flagellated, phylogenetically related species of Lysobacter, and introduced them individually into L. enzymogenes to mimic genomic reacquisition of lost flagellar genes. Heterogeneous expression of the three X. oryzae flagellar structural genes (Xo-motA, Xo-motB, Xo-fliE) and one flagellar regulatory gene (Xo-fleQ) remarkably weakened the bacterial capacity to kill fungal pathogens by impairing the synthesis of an antifungal weapon, known as the heat-stable antifungal factor (HSAF). We further investigated the underlying mechanism by selecting Xo-FleQ as the representative because it is a master transcription factor responsible for flagellar gene expression. Xo-FleQ inhibited the transcription of operon genes responsible for HSAF synthesis via direct binding of Xo-FleQ to the promoter region, thereby decreasing HSAF biosynthesis by L. enzymogenes. These observations suggest a possible genome and function coevolution event, in which an antifungal bacterium deletes certain flagellar genes in order to enhance its ability to kill fungi. IMPORTANCE It is generally recognized that flagellar genes are commonly responsible for the flagella-driven bacterial motility. Thus, finding nonflagellated bacteria partially or fully lost flagellar genes is not a surprise. However, the present study provides new insights into this common idea. We found that loss of either certain flagellar structural or regulatory genes (such as motA, motB, fliE, and fleQ) allows a nonflagellated, antifungal bacterium (L. enzymogenes) to stimulate its fungal-killing capacity, outlining a genome-function coevolution event, where an antifungal bacterium "smartly" designed its genome to "delete" crucial flagellar genes to coordinate flagellar loss and fungal predation. This unusual finding might trigger bacteriologists to reconsider previously ignored functions of the lost flagellar genes in any nonflagellated, pathogenic, or beneficial bacteria.

Secondary symbionts affect aphid fitness and the titer of primary symbiont

Bacterial symbionts associated with aphids are important for their ecological fitness. The corn leaf aphid, Rhopalosiphum maidis (Fitch), is one of the most damaging aphid pests on maize and has been reported to harbor Hamiltonella defensa and Regiella insecticola while the effects of the secondary symbionts (S-symbionts) on host ecology and primary symbiont Buchnera aphidicola remain unclear. Here, four aphid strains were established, two of which were collected from Langfang - Hebei Province, China, with similar symbiont pattern except for the presence of H. defensa. Two other aphid strains were collected from Nanning - Guangxi Province, China, with the same symbiont infection except for the presence of R. insecticola. Phylogenetic analysis and aphid genotyping indicated that the S-symbiont-infected and free aphid strains from the same location had identical genetic backgrounds. Aphid fitness measurement showed that aphid strain infected with H. defensa performed shortened developmental duration for 1^st instar and total nymph stages, reduced aphid survival rate, offspring, and longevity. While the developmental duration of H-infected strains was accelerated, and the adult weight was significantly higher compared to the H-free strain. Infection with R. insecticola did not affect the aphid's entire nymph stage duration and survival rate. As the H-strain does, aphids infected with R. insecticola also underwent a drop in offspring, along with marginally lower longevity. Unlike the H-infected strain, the R-infected strain performed delayed developmental duration and lower adult weight. The B. aphidicola titers of the H-infected strains showed a steep drop during the aphid 1^st to 3^rd instar stages, while the augmentation of B. aphidicola titers was found in the R-infected strain during the aphid 1^st to 3^rd instar. Our study investigated for the first time the effect of the S-symbionts on the ecology fitness and primary symbiont in R. maidis, indicating that infection with secondary symbionts leads to the modulation of aphid primary symbiont abundance, together inducing significant fitness costs on aphids with further impact on environmental adaptation and trophic interactions.

Complex host/symbiont integration of a multi-partner symbiotic system in the eusocial aphid Ceratovacuna japonica

Some hemipteran insects rely on multiple endosymbionts for essential nutrients. However, the evolution of multi-partner symbiotic systems is not well-established. Here, we report a co-obligate symbiosis in the eusocial aphid, Ceratovacuna japonica. 16S rRNA amplicon sequencing unveiled co-infection with a novel Arsenophonus sp. symbiont and Buchnera aphidicola, a common obligate endosymbiont in aphids. Both symbionts were housed within distinct bacteriocytes and were maternally transmitted. The Buchnera and Arsenophonus symbionts had streamlined genomes of 432,286 bp and 853,149 bp, respectively, and exhibited metabolic complementarity in riboflavin and peptidoglycan synthesis pathways. These anatomical and genomic properties were similar to those of independently evolved multi-partner symbiotic systems, such as Buchnera-Serratia in Lachninae and Periphyllus aphids, representing remarkable parallelism. Furthermore, symbiont populations and bacteriome morphology differed between reproductive and soldier castes. Our study provides the first example of co-obligate symbiosis in Hormaphidinae and gives insight into the evolutionary genetics of this complex system.

Microbiome of the Aphid Genus Dysaphis Börner (Hemiptera: Aphidinae) and Its Relation to Ant Attendance

Among mutualistic relationships of aphids with other organisms, there are two that seem to be of major importance: trophobiosis with ants and endosymbiosis of bacteria. While the former is well studied, the latter is the subject of an increasing amount of research constantly revealing new aspects of this symbiosis. Here, we studied the possible influence of ant attendance on the composition of aphid microbiota on primary and secondary hosts exploited by the aphid genus Dysaphis. The microbiome of 44 samples representing 12 aphid species was studied using an Illumina HiSeq 4000 with the V3-V4 region of 16S rRNA. The results showed a higher abundance of common facultative symbionts (Serratia, Regiella, Fukatsuia) in aphid species unattended by ants, but also on secondary hosts. However, in colonies attended by ants, the general species composition of bacterial symbionts was more rich in genera than in unattended colonies (Wolbachia, Gilliamella, Spiroplasma, Sphingomonas, Pelomonas). The results indicate a huge variability of facultative symbionts without clear correlation with ant attendance or aphid species. The possibility of multiple routes of bacterial infection mediated by ant-made environmental conditions is discussed.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

Host’s demand for essential amino acids is compensated by an extracellular bacterial symbiont in a hemipteran insect model

Plant sap is a nutritionally unbalanced diet that constitutes a challenge for insects that feed exclusively on it. Sap-sucking hemipteran insects generally overcome this challenge by harboring beneficial microorganisms in their specialized symbiotic organ, either intracellularly or extracellularly. Genomic information of these bacterial symbionts suggests that their primary role is to supply essential amino acids, but empirical evidence has been virtually limited to the intracellular symbiosis between aphids and Buchnera. Here we investigated the amino acid complementation by the extracellular symbiotic bacterium Ishikawaella harbored in the midgut symbiotic organ of the stinkbug Megacopta punctatissima. We evaluated amino acid compositions of the phloem sap of plants on which the insect feeds, as well as those of its hemolymph, whole body hydrolysate, and excreta. The results highlighted that the essential amino acids in the diet are apparently insufficient for the stinkbug development. Experimental symbiont removal caused severe shortfalls of some essential amino acids, including branched-chain and aromatic amino acids. In vitro culturing of the isolated symbiotic organ demonstrated that hemolymph-circulating metabolites, glutamine and trehalose, efficiently fuel the production of essential amino acids. Branched-chain amino acids and aromatic amino acids are the ones preferentially synthesized despite the symbiont’s synthetic capability of all essential amino acids. These results indicate that the symbiont-mediated amino acid compensation is quantitatively optimized in the stinkbug-Ishikawaella gut symbiotic association as in the aphid-Buchnera intracellular symbiotic association. The convergence of symbiont functions across distinct nutritional symbiotic systems provides insight into how host-symbiont interactions have been shaped over evolutionary time.

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.

A pluralistic view of holobionts in the context of process ontology

Developing precise definitions and fine categories is an important part of the scientific endeavour, enabling fidelity of transfers of knowledge and the progress of science. Currently, as a result of research on symbiotic microorganisms, science has been flooded with discoveries which appear to undermine many commonly accepted concepts and to introduce new ones that often require updated conceptualisations. One question currently being debated concerns whether or not a holobiont can be considered an organism. Based on which concept, physiology or evolutionary, of the organism is chosen, the verdict differs. We attempt here to show how a change in perspective, from that of substance ontology into that of process ontology, is capable of reconciling opposing positions within the existing discussion and enabling the implementation of conceptual pluralism.

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.

Impact of host demography and evolutionary history on endosymbiont molecular evolution: A test in carpenter ants (genus Camponotus) and their Blochmannia endosymbionts

Obligate endosymbioses are tight associations between symbionts and the hosts they live inside. Hosts and their associated obligate endosymbionts generally exhibit codiversification, which has been documented in taxonomically diverse insect lineages. Host demography (e.g., effective population sizes) may impact the demography of endosymbionts, which may lead to an association between host demography and the patterns and processes of endosymbiont molecular evolution. Here, we used whole-genome sequencing data for carpenter ants (Genus Camponotus; subgenera Camponotus and Tanaemyrmex) and their Blochmannia endosymbionts as our study system to address whether Camponotus demography shapes Blochmannia molecular evolution. Using whole-genome phylogenomics, we confirmed previous work identifying codiversification between carpenter ants and their Blochmannia endosymbionts. We found that Blochmannia genes have evolved at a pace ~30× faster than that of their hosts' molecular evolution and that these rates are positively associated with host rates of molecular evolution. Using multiple tests for selection in Blochmannia genes, we found signatures of positive selection and shifts in selection strength across the phylogeny. Host demography was associated with Blochmannia shifts toward increased selection strengths, but not associated with Blochmannia selection relaxation, positive selection, genetic drift rates, or genome size evolution. Mixed support for relationships between host effective population sizes and Blochmannia molecular evolution suggests weak or uncoupled relationships between host demography and Blochmannia population genomic processes. Finally, we found that Blochmannia genome size evolution was associated with genome-wide estimates of genetic drift and number of genes with relaxed selection pressures.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel's lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

Elucidation of host and symbiont contributions to peptidoglycan metabolism based on comparative genomics of eight aphid subfamilies and their Buchnera

Pea aphids (Acyrthosiphon pisum) are insects containing genes of bacterial origin with putative functions in peptidoglycan (PGN) metabolism. Of these, rlpA1-5, amiD, and ldcA are highly expressed in bacteriocytes, specialized aphid cells that harbor the obligate bacterial symbiont Buchnera aphidicola, required for amino acid supplementation of the host's nutrient-poor diet. Despite genome reduction associated with endosymbiosis, pea aphid Buchnera retains genes for the synthesis of PGN while Buchnera of many other aphid species partially or completely lack these genes. To explore the evolution of aphid horizontally-transferred genes (HTGs) and to elucidate how host and symbiont genes contribute to PGN production, we sequenced genomes from four deeply branching lineages, such that paired aphid and Buchnera genomes are now available for 17 species representing eight subfamilies. We identified all host and symbiont genes putatively involved in PGN metabolism. Phylogenetic analyses indicate that each HTG family was present in the aphid shared ancestor, but that each underwent a unique pattern of gene loss or duplication in descendant lineages. While four aphid rlpA gene subfamilies show no relation to symbiont PGN gene repertoire, the loss of aphid amiD and ldcA HTGs coincides with the loss of symbiont PGN metabolism genes. In particular, the coincident loss of host amiD and symbiont murCEF in tribe Aphidini, in contrast to tribe Macrosiphini, suggests either 1) functional linkage between these host and symbiont genes, or 2) Aphidini has lost functional PGN synthesis and other retained PGN pathway genes are non-functional. To test these hypotheses experimentally, we used cell-wall labeling methods involving a d-alanine probe and found that both Macrosiphini and Aphidini retain Buchnera PGN synthesis. Our results imply that compensatory adaptations can preserve PGN synthesis despite the loss of some genes considered essential for this pathway, highlighting the importance of the cell wall in these symbioses.

Recently introduced Wolbachia reduces bacterial species richness and reshapes bacterial community structure in Nilaparvata lugens

BACKGROUND

Wolbachia has been developed as an effective tool to suppress insect pests and arbovirus transmission. Recently, the brown planthopper Nilaparvata lugens, a serious agricultural pest, has been successfully transinfected with Wolbachia wStri strain from Laodelphax striatellus. However, before conducting the field experiments, the impacts of wStri on the bacterial microbiota in N. lugens and how it differs from native Wolbachia wLug strain have not been clarified.

RESULTS

Here, we found that wStri reduced bacterial diversity and shaped bacterial community structure more than wLug in both developmental stage and different adult tissues. Overall, the relative abundance of Wolbachia was negatively correlated with bacterial diversity, but the bacterial diversity gradually decreased only when the relative abundance of Wolbachia was higher than 60%. Further analysis found that wStri reduced species richness of other bacteria but not their evenness. wStri infection also affected many bacterial functions (e.g., amino acid metabolism & signaling and cellular processes) in the developmental stages, with a stronger effect than wLug in nymphs. Moreover, although Wolbachia occupied a high relative abundance in infected individuals, Acinetobacter was consistently a core part of microbiome.

CONCLUSION

These results showed the significant impacts of recently introduced wStri on bacterial microbiota in N. lugens, with the effects differing from native wLug. This study will aid in understanding the relationship between Wolbachia, its host and the host's microbiota, and provide a reference for future field experiments.

Aphid hologenomics: current status and future challenges

Aphids are important model organisms in ecological, developmental, and evolutionary studies of, for example, symbiosis, insect-plant interactions, pest management, and developmental polyphenism. Here, we review the recent progress made in the genomics of aphids and their symbionts: hologenomics. The reference genome of Acyrthosiphon pisum has been greatly improved, and chromosome-level assembly is now available. The genomes of over 20 aphid species have been sequenced, and comparative genomic analyses have revealed pervasive gene duplication and dynamic chromosomal rearrangements. Over 120 symbiont genomes (both obligate and facultative) have been sequenced, and modern deep-sequencing technologies have identified novel symbionts. The advances in hologenomics have helped to elucidate the dynamic evolution of facultative and co-obligate symbionts with the ancient obligate symbiont Buchnera.

Genomic evolution and adaptation of arthropod-associated Rickettsia

Rickettsia species are endosymbionts hosted by arthropods and are known to cause mild to fatal diseases in humans. Here, we analyse the evolution and diversity of 34 Rickettsia species using a pangenomic meta-analysis (80 genomes/41 plasmids). Phylogenomic trees showed that Rickettsia spp. diverged into two Spotted Fever groups, a Typhus group, a Canadensis group and a Bellii group, and may have inherited their plasmids from an ancestral plasmid that persisted in some strains or may have been lost by others. The results suggested that the ancestors of Rickettsia spp. might have infected Acari and/or Insecta and probably diverged by persisting inside and/or switching hosts. Pangenomic analysis revealed that the Rickettsia genus evolved through a strong interplay between genome degradation/reduction and/or expansion leading to possible distinct adaptive trajectories. The genus mainly shared evolutionary relationships with α-proteobacteria, and also with γ/β/δ-proteobacteria, cytophagia, actinobacteria, cyanobacteria, chlamydiia and viruses, suggesting lateral exchanges of several critical genes. These evolutionary processes have probably been orchestrated by an abundance of mobile genetic elements, especially in the Spotted Fever and Bellii groups. In this study, we provided a global evolutionary genomic view of the intracellular Rickettsia that may help our understanding of their diversity, adaptation and fitness.

Genome sequence of the English grain aphid, Sitobion avenae and its endosymbiont Buchnera aphidicola

The English grain aphid, Sitobion avenae, is a major agricultural pest of wheat, barley and oats, and one of the principal vectors of barley yellow dwarf virus leading to significant reductions in grain yield, annually. Emerging resistance to and increasing regulation of insecticides has resulted in limited options for their control. Using PacBio HiFi data, we have produced a high-quality draft assembly of the S. avenae genome; generating a primary assembly with a total assembly size of 475.7 Mb, and an alternate assembly with a total assembly size of 430.8 Mb. Our primary assembly was highly contiguous with only 326 contigs and a contig N50 of 15.95 Mb. Assembly completeness was estimated at 97.7% using BUSCO analysis and 31,007 and 29,037 protein-coding genes were predicted from the primary and alternate assemblies, respectively. This assembly, which is to our knowledge the first for an insecticide resistant clonal lineage of English grain aphid, will provide novel insight into the molecular and mechanistic determinants of resistance and will facilitate future research into mechanisms of viral transmission and aphid behavior.

The Phylosymbiosis Pattern Between the Fig Wasps of the Same Genus and Their Associated Microbiota

Microbial communities can be critical for many metazoans, which can lead to the observation of phylosymbiosis with phylogenetically related species sharing similar microbial communities. Most of the previous studies on phylosymbiosis were conducted across the host families or genera. However, it is unclear whether the phylosymbiosis signal is still prevalent at lower taxonomic levels. In this study, 54 individuals from six species of the fig wasp genus Ceratosolen (Hymenoptera: Agaonidae) collected from nine natural populations and their associated microbiota were investigated. The fig wasp species were morphologically identified and further determined by mitochondrial CO1 gene fragments and nuclear ITS2 sequences, and the V4 region of 16S rRNA gene was sequenced to analyze the bacterial communities. The results suggest a significant positive correlation between host genetic characteristics and microbial diversity characteristics, indicating the phylosymbiosis signal between the phylogeny of insect hosts and the associated microbiota in the lower classification level within a genus. Moreover, we found that the endosymbiotic Wolbachia carried by fig wasps led to a decrease in bacterial diversity of host-associated microbial communities. This study contributes to our understanding of the role of host phylogeny, as well as the role of endosymbionts in shaping the host-associated microbial community.

Genomic features underlying the evolutionary transitions of Apibacter to honey bee gut symbionts

The gut bacteria of honey bee recognized as a mutualistic partner with the insect host might have originated from a free-living or parasitic lifestyle. However, little is known about the genomic features underlying this lifestyle transition. Here we compared the genomes of bee gut bacteria Apibacter with their close relatives living in different lifestyles. We found that despite general reduction in the Apibacter genome, genes involved in amino acid synthesis and monosaccharide detoxification were retained, which is putatively beneficial to the host. Interestingly, the microaerobic Apibacter species specifically acquired genes encoding for the nitrate respiration (NAR). These together with nitrate transporter and enzymatic cofactor synthesis genes were found clustered in the genomes. The NAR system is also conserved in the cohabitating bee gut microbe Snodgrassella, although with a different structure. This convergence suggests a key role of respiratory nitrate reduction for microaerophilic microbiomes to colonize bee gut epithelium. Genes involved in lipid, histidine degradation were found partially or completely lost in Apibacter. Particularly, genes encoding for the conversion to the toxic intermediates in phenylacetate degradation, as well as other potential virulence factors, are specifically lost in Apibacter group. Antibiotic resistance genes are only sporadically distributed among Apibacter species, but are prevalent in their relatives, which may be related to the remotely living feature and less exposure to antibiotics of their bee hosts. Collectively, this study advanced our knowledge of genomic features specialized to bee gut symbionts.

Microbe Profile: Buchnera aphidicola: ancient aphid accomplice and endosymbiont exemplar

Buchnera aphidicola is an obligate endosymbiont of aphids that cannot be cultured outside of hosts. It exists as diverse strains in different aphid species, and phylogenetic reconstructions show that it has been maternally transmitted in aphids for >100 million years. B. aphidicola genomes are highly reduced and show conserved gene order and no gene acquisition, but encoded proteins undergo rapid evolution. Aphids depend on B. aphidicola for biosynthesis of essential amino acids and as an integral part of embryonic development. How B. aphidicola populations are regulated within hosts remains little known.

Horizontal-Acquisition of a Promiscuous Peptidoglycan-Recycling Enzyme Enables Aphids To Influence Symbiont Cell Wall Metabolism

During evolution, enzymes can undergo shifts in preferred substrates or in catalytic activities. An intriguing question is how enzyme function changes following horizontal gene transfer, especially for bacterial genes that have moved to animal genomes. Some insects have acquired genes that encode enzymes for the biosynthesis of bacterial cell wall components and that appear to function to support or control their obligate endosymbiotic bacteria. In aphids, the bacterial endosymbiont Buchnera aphidicola provides essential amino acids for aphid hosts but lacks most genes for remodeling of the bacterial cell wall. The aphid genome has acquired seven genes with putative functions in cell wall metabolism that are primarily expressed in the aphid cells harboring Buchnera. In analyses of aphid homogenates, we detected peptidoglycan (PGN) muropeptides indicative of the reactions of PGN hydrolases encoded by horizontally acquired aphid genes but not by Buchnera genes. We produced one such host enzyme, ApLdcA, and characterized its activity with both cell wall derived and synthetic PGN. Both ApLdcA and the homologous enzyme in Escherichia coli, which functions as an l,d-carboxypeptidase in the cytoplasmic PGN recycling pathway, exhibit turnover of PGN substrates containing stem pentapeptides and cross-linkages via l,d-endopeptidase activity, consistent with a potential role in cell wall remodeling. Our results suggest that ApLdcA derives its functions from the promiscuous activities of an ancestral LdcA enzyme, whose acquisition by the aphid genome may have enabled hosts to influence Buchnera cell wall metabolism as a means to control symbiont growth and division. IMPORTANCE Most enzymes are capable of performing biologically irrelevant side reactions. During evolution, promiscuous enzyme activities may acquire new biological roles, especially after horizontal gene transfer to new organisms. Pea aphids harbor obligate bacterial symbionts called Buchnera and encode horizontally acquired bacterial genes with putative roles in cell wall metabolism. Though Buchnera lacks cell wall endopeptidase genes, we found evidence of endopeptidase activity among peptidoglycan muropeptides purified from aphids. We characterized a multifunctional, aphid-encoded enzyme, ApLdcA, which displays l,d-endopeptidase activities considered promiscuous for the Escherichia coli homolog, for which these activities do not contribute to its native role in peptidoglycan recycling. These results exemplify the roles of enzyme promiscuity and horizontal gene transfer in enzyme evolution and demonstrate how aphids influence symbiont cell wall metabolism.

Genome Evolution of Filamentous Cyanobacterium Nostoc Species: From Facultative Symbiosis to Free Living

In contrast to obligate bacteria, facultative symbiotic bacteria are mainly characterized by genome enlargement. However, the underlying relationship of this feature with adaptations to various habitats remains unclear. In this study, we used the global genome data of Nostoc strains, including 10 novel genomes sequenced in this study and 26 genomes available from public databases, and analyzed their evolutionary history. The evolutionary boundary of the real clade of Nostoc species was identified and was found to be consistent with the results of polyphasic taxonomy. The initial ancestral species of Nostoc was demonstrated to be consistent with a facultative symbiotic population. Further analyses revealed that Nostoc strains tended to shift from facultative symbiosis to a free-living one, along with an increase in genome sizes during the dispersal of each exterior branch. Intracellular symbiosis was proved to be essentially related to Nostoc evolution, and the adaptation of its members to free-living environments was coupled with a large preference for gene acquisition involved in gene repair and recombination. These findings provided unique evidence of genomic mechanisms by which homologous microbes adapt to distinct life manners and external environments.

Social environment drives sex and age-specific variation in Drosophila melanogaster microbiome composition and predicted function

Social environments influence multiple traits of individuals including immunity, stress and ageing, often in sex-specific ways. The composition of the microbiome (the assemblage of symbiotic microorganisms within a host) is determined by environmental factors and the host's immune, endocrine and neural systems. The social environment could alter host microbiomes extrinsically by affecting transmission between individuals, probably promoting homogeneity in the microbiome of social partners. Alternatively, intrinsic effects arising from interactions between the microbiome and host physiology (the microbiota-gut-brain axis) could translate social stress into dysbiotic microbiomes, with consequences for host health. We investigated how manipulating social environments during larval and adult life-stages altered the microbiome composition of Drosophila melanogaster fruit flies. We used social contexts that particularly alter the development and lifespan of males, predicting that any intrinsic social effects on the microbiome would therefore be sex-specific. The presence of adult males during the larval stage significantly altered the microbiome of pupae of both sexes. In adults, same-sex grouping increased bacterial diversity in both sexes. Importantly, the microbiome community structure of males was more sensitive to social contact at older ages, an effect partially mitigated by housing focal males with young rather than coaged groups. Functional analyses suggest that these microbiome changes impact ageing and immune responses. This is consistent with the hypothesis that the substantial effects of the social environment on individual health are mediated through intrinsic effects on the microbiome, and provides a model for understanding the mechanistic basis of the microbiota-gut-brain axis.

Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids

Numerous animal lineages have maternally inherited symbionts that are required for host reproduction and growth. Endosymbionts also pose a risk to their hosts because of the mutational decay of their genomes through genetic drift or to selfish mutations that favor symbiont fitness over host fitness. One model for heritable endosymbiosis is the association of aphids with their obligate bacterial symbiont, Buchnera We experimentally established heteroplasmic pea aphid matrilines containing pairs of closely related Buchnera haplotypes and used deep sequencing of diagnostic markers to measure haplotype frequencies in successive host generations. These frequencies were used to estimate the effective population size of Buchnera within hosts (i.e., the transmission bottleneck size) and the extent of within-host selection. The within-host effective population size was in the range of 10 to 20, indicating a strong potential for genetic drift and fixation of deleterious mutations. Remarkably, closely related haplotypes were subject to strong within-host selection, with selection coefficients as high as 0.5 per aphid generation. In one case, the direction of selection depended on the thermal environment and went in the same direction as between-host selection. In another, a new mutant haplotype had a strong within-host advantage under both environments but had no discernible effect on host-level fitness under laboratory conditions. Thus, within-host selection can be strong, resulting in a rapid fixation of mutations with little impact on host-level fitness. Together, these results show that within-host selection can drive evolution of an obligate symbiont, accelerating sequence evolution.

'Candidatus Liberibacter asiaticus' Multimeric LotP Mediates Citrus sinensis Defense Response Activation

‘Candidatus Liberibacter asiaticus’ is known as the most pathogenic organism associated with citrus greening disease. Since its publicized emergence in Florida in 2005, ‘Ca. L. asiaticus’ remains unculturable. Currently, a limited number of potential disease effectors have been identified through in silico analysis. Therefore, these potential effectors remain poorly characterized and do not fully explain the complexity of symptoms observed in citrus trees infected with ‘Ca. L. asiaticus.’ LotP has been identified as a potential effector and have been partially characterized. This protein retains structural homology to the substrate binding domain of the Lon protease. LotP interacts with chaperones like GroEL, Hsp40, DnaJ, and ClpX and may exercise its biological role through interactions with different proteins involved in proteostasis networks. Here, we evaluate the interactome of LotP—revealing a new protein–protein interaction target (Lon-serine protease) and its effect on citrus plant tissue integrity. We found that via protein–protein interactions, LotP can enhance Lon protease activity, increasing the degradation rate of its specific targets. Infiltration of purified LotP strained citrus plant tissue causing photoinhibition and chlorosis after several days. Proteomics analysis of LotP tissues recovering after the infiltration revealed a large abundance of plant proteins associated with the stabilization and processing of mRNA transcripts, a subset of important transcription factors; and pathways associated with innate plant defense were highly expressed. Furthermore, interactions and substrate binding module of LotP suggest potential interactions with plant proteins, most likely proteases.

Can gene-inactivating mutations lead to evolutionary novelty?

Evolutionary novelty is difficult to define. It typically involves shifts in organismal or biochemical phenotypes that can be seen as qualitative as well as quantitative changes. In laboratory-based experimental evolution of novel phenotypes and the human domestication of crops, the majority of the mutations that lead to adaptation are loss-of-function mutations that impair or eliminate the function of genes rather than gain-of-function mutations that increase or qualitatively alter the function of proteins. Here, I speculate that easier access to loss-of-function mutations has led them to play a major role in the adaptive radiations that occur when populations have access to many unoccupied ecological niches. I discuss five possible objections to this claim: that genes can only survive if they confer benefits to the organisms that bear them, antagonistic pleiotropy, the importance of pre-existing genetic variation in populations, the danger that adaptation by breaking genes will, over long times, cause organisms to run out of genes, and the recessive nature of most loss-of-function mutations.

Genetic innovations in animal-microbe symbioses

Animal hosts have initiated myriad symbiotic associations with microorganisms and often have maintained these symbioses for millions of years, spanning drastic changes in ecological conditions and lifestyles. The establishment and persistence of these relationships require genetic innovations on the parts of both symbionts and hosts. The nature of symbiont innovations depends on their genetic population structure, categorized here as open, closed or mixed. These categories reflect modes of inter-host transmission that result in distinct genomic features, or genomic syndromes, in symbionts. Although less studied, hosts also innovate in order to preserve and control symbiotic partnerships. New capabilities to sequence host-associated microbial communities and to experimentally manipulate both hosts and symbionts are providing unprecedented insights into how genetic innovations arise under different symbiont population structures and how these innovations function to support symbiotic relationships.

Isolation of the Buchnera aphidicola flagellum basal body complexes from the Buchnera membrane

Buchnera aphidicola is an intracellular bacterial symbiont of aphids and maintains a small genome of only 600 kbps. Buchnera is thought to maintain only genes relevant to the symbiosis with its aphid host. Curiously, the Buchnera genome contains gene clusters coding for flagellum basal body structural proteins and for flagellum type III export machinery. These structures have been shown to be highly expressed and present in large numbers on Buchnera cells. No recognizable pathogenicity factors or secreted proteins have been identified in the Buchnera genome, and the relevance of this protein complex to the symbiosis is unknown. Here, we show isolation of Buchnera flagellum basal body proteins from the cellular membrane of Buchnera, confirming the enrichment of flagellum basal body proteins relative to other proteins in the Buchnera proteome. This will facilitate studies of the structure and function of the Buchnera flagellum structure, and its role in this model symbiosis.

Vertical Transmission at the Pathogen-Symbiont Interface: Serratia symbiotica and Aphids

Insects have evolved various mechanisms to reliably transmit their beneficial bacterial symbionts to the next generation. Sap-sucking insects, including aphids, transmit symbionts by endocytosis of the symbiont into cells of the early embryo within the mother’s body.

Many insects possess beneficial bacterial symbionts that occupy specialized host cells and are maternally transmitted. As a consequence of their host-restricted lifestyle, these symbionts often possess reduced genomes and cannot be cultured outside hosts, limiting their study. The bacterial species Serratia symbiotica was originally characterized as noncultured strains that live as mutualistic symbionts of aphids and are vertically transmitted through transovarial endocytosis within the mother’s body. More recently, culturable strains of S. symbiotica were discovered that retain a larger set of ancestral Serratia genes, are gut pathogens in aphid hosts, and are principally transmitted via a fecal-oral route. We find that these culturable strains, when injected into pea aphids, replicate in the hemolymph and are pathogenic. Unexpectedly, they are also capable of maternal transmission via transovarial endocytosis: using green fluorescent protein (GFP)-tagged strains, we observe that pathogenic S. symbiotica strains, but not Escherichia coli, are endocytosed into early embryos. Furthermore, pathogenic S. symbiotica strains are compartmentalized into specialized aphid cells in a fashion similar to that of mutualistic S. symbiotica strains during later stages of embryonic development. However, infected embryos do not appear to develop properly, and offspring infected by a transovarial route are not observed. Thus, cultured pathogenic strains of S. symbiotica have the latent capacity to transition to lifestyles as mutualistic symbionts of aphid hosts, but persistent vertical transmission is blocked by their pathogenicity. To transition into stably inherited symbionts, culturable S. symbiotica strains may need to adapt to regulate their titer, limit their pathogenicity, and/or provide benefits to aphids that outweigh their cost.

Evolutionary Rates are Correlated Between Buchnera Endosymbionts and the Mitochondrial Genomes of Their Aphid Hosts

The evolution of bacterial endosymbiont genomes is strongly influenced by host-driven selection. Factors affecting host genome evolution will potentially affect endosymbiont genomes in similar ways. One potential outcome is correlations in molecular rates between the genomes of the symbiotic partners. Recently, we presented the first evidence of such correlations between the mitochondrial genomes of cockroaches and the genomes of their endosymbiont (Blattabacterium cuenoti). Here we investigate whether similar patterns are found in additional host-symbiont partners. We use partial genome data from multiple strains of the bacterial endosymbionts Buchnera aphidicola and Sulcia muelleri, and the mitochondrial genomes of their sap-feeding insect hosts. Both endosymbionts show phylogenetic congruence with the mitochondria of their hosts, a result that is expected due to their identical mode of inheritance. We compared root-to-tip distances and branch lengths of phylogenetically independent species pairs. Both analyses showed a highly significant correlation of molecular rates between the genomes of Buchnera and the mitochondrial genomes of their hosts. A similar correlation was detected between Sulcia and their hosts, but was not statistically significant. Our results indicate that evolutionary rate correlations between hosts and long-term symbionts may be a widespread phenomenon.

Transmission of the sponge microbiome: moving towards a unified model

Sponges have co-evolved for millions of years alongside several types of microorganisms, which aside from participating in the animal's diet, are mostly symbionts. Since most of the genetic repertoire in the holobiont genome is provided by microbes, it is expected that the host-associated microbiome will be at least partially heritable. Sponges can therefore acquire their symbionts in different ways. Both vertical transmission (VT) and horizontal transmission (HT) have different advantages and disadvantages in the life cycle of these invertebrates. However, a third mode of transmission, called leaky vertical transmission or mixed mode of transmission (MMT), which incorporates both VT and HT modes, has gained relevance and seems to be the most robust model. In that regard, the aim of this review is to present the evolving knowledge on these main modes of transmission of the sponge microbiome. Our conclusions lead us to suggest that MMT may be more common for all sponges, with its frequency varying across the transmission spectrum between species and the environment. This hybrid model supports the stable and specific transmission of these microbial partners and reinforces their assistance in the resilience of sponges over the years.

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.