Multiple origins of endosymbiosis within the Enterobacteriaceae (γ-Proteobacteria): convergence of complex phylogenetic approaches

Discovery of a novel symbiotic lineage associated with a hematophagous leech from the genus Haementeria

Similarly to other strict blood feeders, leeches from the Haementeria genus (Hirudinida: Glossiphoniidae) have established a symbiotic association with bacteria harbored intracellularly in esophageal bacteriomes. Previous genome sequence analyses of these endosymbionts revealed co-divergence with their hosts, a strong genome reduction, and a simplified metabolism largely dedicated to the production of B vitamins, which are nutrients lacking from a blood diet. 'Candidatus Providencia siddallii' has been identified as the obligate nutritional endosymbiont of a monophyletic clade of Mexican and South American Haementeria spp. However, the Haementeria genus includes a sister clade of congeners from Central and South America, where the presence or absence of the aforementioned symbiont taxon remains unknown. In this work, we report on a novel bacterial endosymbiont found in a representative from this Haementeria clade. We found that this symbiont lineage has evolved from within the Pluralibacter genus, known mainly from clinical but also environmental strains. Similarly to Ca. Providencia siddallii, the Haementeria-associated Pluralibacter symbiont displays clear signs of genome reduction, accompanied by an A+T-biased sequence composition. Genomic analysis of its metabolic potential revealed a retention of pathways related to B vitamin biosynthesis, supporting its role as a nutritional endosymbiont. Finally, comparative genomics of both Haementeria symbiont lineages suggests that an ancient Providencia symbiont was likely replaced by the novel Pluralibacter one, thus constituting the first reported case of nutritional symbiont replacement in a leech without morphological changes in the bacteriome.

IMPORTANCE

Obligate symbiotic associations with a nutritional base have likely evolved more than once in strict blood-feeding leeches. Unlike those symbioses found in hematophagous arthropods, the nature, identity, and evolutionary history of these remains poorly studied. In this work, we further explored obligate nutritional associations between Haementeria leeches and their microbial symbionts, which led to the unexpected discovery of a novel symbiosis with a member of the Pluralibacter genus. When compared to Providencia siddallii, an obligate nutritional symbiont of other Haementeria leeches, this novel bacterial symbiont shows convergent retention of the metabolic pathways involved in B vitamin biosynthesis. Moreover, the genomic characteristics of this Pluralibacter symbiont suggest a more recent association than that of Pr. siddallii and Haementeria. We conclude that the once-thought stable associations between blood-feeding Glossiphoniidae and their symbionts (i.e., one bacteriome structure, one symbiont lineage) can break down, mirroring symbiont turnover observed in various arthropod lineages.

The evolution of entomopathogeny in nematodes

Understanding how parasites evolved is crucial to understand the host and parasite interaction. The evolution of entomopathogenesis in rhabditid nematodes has traditionally been thought to have occurred twice within the phylum Nematoda: in Steinernematidae and Heterorhabditidae families, which are associated with the entomopathogenic bacteria Xenorhabdus and Photorhabdus, respectively. However, nematodes from other families that are associated with entomopathogenic bacteria have not been considered to meet the criteria for "entomopathogenic nematodes." The evolution of parasitism in nematodes suggests that ecological and evolutionary properties shared by families in the order Rhabditida favor the convergent evolution of the entomopathogenic trait in lineages with diverse lifestyles, such as saprotrophs, phoretic, and necromenic nematodes. For this reason, this paper proposes expanding the term "entomopathogenic nematode" considering the diverse modes of this attribute within Rhabditida. Despite studies are required to test the authenticity of the entomopathogenic trait in the reported species, they are valuable links that represent the early stages of specialized lineages to entomopathogenic lifestyle. An ecological and evolutionary exploration of these nematodes has the potential to deepen our comprehension of the evolution of entomopathogenesis as a convergent trait spanning across the Nematoda.

Incongruence in the phylogenomics era

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Convergent evolution of a labile nutritional symbiosis in ants

Ants are among the most successful organisms on Earth. It has been suggested that forming symbioses with nutrient-supplementing microbes may have contributed to their success, by allowing ants to invade otherwise inaccessible niches. However, it is unclear whether ants have evolved symbioses repeatedly to overcome the same nutrient limitations. Here, we address this question by comparing the independently evolved symbioses in Camponotus, Plagiolepis, Formica and Cardiocondyla ants. Our analysis reveals the only metabolic function consistently retained in all of the symbiont genomes is the capacity to synthesise tyrosine. We also show that in certain multi-queen lineages that have co-diversified with their symbiont for millions of years, only a fraction of queens carry the symbiont, suggesting ants differ in their colony-level reliance on symbiont-derived resources. Our results imply that symbioses can arise to solve common problems, but hosts may differ in their dependence on symbionts, highlighting the evolutionary forces influencing the persistence of long-term endosymbiotic mutualisms.

Lightella neohaematopini: A new lineage of highly reduced endosymbionts coevolving with chipmunk lice of the genus Neohaematopinus

Sucking lice (Anoplura) are known to have established symbiotic associations multiple times with different groups of bacteria as diverse as Enterobacteriales, Legionellales, and Neisseriales. This diversity, together with absence of a common coevolving symbiont (such as Buchnera, in aphids), indicates that sucking lice underwent a series of symbiont acquisitions, losses, and replacements. To better understand evolution and significance of louse symbionts, genomic and phylogenetic data are needed from a broader taxonomic diversity of lice and their symbiotic bacteria. In this study, we extend the known spectrum of the louse symbionts with a new lineage associated with Neohaematopinus pacificus, a louse species that commonly parasitizes North American chipmunks. The recent coevolutionary analysis showed that rather than a single species, these lice form a cluster of unique phylogenetic lineages specific to separate chipmunk species (or group of closely related species). Using metagenomic assemblies, we show that the lice harbor a bacterium which mirrors their phylogeny and displays traits typical for obligate mutualists. Phylogenetic analyses place this bacterium within Enterobacteriaceae on a long branch related to another louse symbiont, “Candidatus Puchtella pedicinophila.” We propose for this symbiotic lineage the name “Candidatus Lightella neohaematopini.” Based on the reconstruction of metabolic pathways, we suggest that like other louse symbionts, L. neohaematopini provides its host with at least some B vitamins. In addition, several samples harbored another symbiotic bacterium phylogenetically affiliated with the Neisseriales-related symbionts described previously from the lice Polyplax serrata and Hoplopleura acanthopus. Characterizing these bacteria further extend the known diversity of the symbiotic associations in lice and show unique complexity and dynamics of the system.

Annotation-free delineation of prokaryotic homology groups

Phylogenomic studies of prokaryotic taxa often assume conserved marker genes are homologous across their length. However, processes such as horizontal gene transfer or gene duplication and loss may disrupt this homology by recombining only parts of genes, causing gene fission or fusion. We show using simulation that it is necessary to delineate homology groups in a set of bacterial genomes without relying on gene annotations to define the boundaries of homologous regions. To solve this problem, we have developed a graph-based algorithm to partition a set of bacterial genomes into Maximal Homologous Groups of sequences (MHGs) where each MHG is a maximal set of maximum-length sequences which are homologous across the entire sequence alignment. We applied our algorithm to a dataset of 19 Enterobacteriaceae species and found that MHGs cover much greater proportions of genomes than markers and, relatedly, are less biased in terms of the functions of the genes they cover. We zoomed in on the correlation between each individual marker and their overlapping MHGs, and show that few phylogenetic splits supported by the markers are supported by the MHGs while many marker-supported splits are contradicted by the MHGs. A comparison of the species tree inferred from marker genes with the species tree inferred from MHGs suggests that the increased bias and lack of genome coverage by markers causes incorrect inferences as to the overall relationship between bacterial taxa.

Assuming genes to be the basic evolutionary unit has been commonplace in bacterial genomics. For example, when quantifying the extent of horizontal gene transfer it is common to infer gene trees and reconcile them against a species tree to account for recombination-based processes. We have developed a new method which challenges this assumption by identifying contiguous regions of true homology without regards to gene boundaries and applied it to Enterobacteriaceae, a family of bacteria containing several important human pathogens. Our results show that genes are composed of distinct homologous regions with conflicting phylogenetic histories. We further demonstrate that failing to take account of this conflict, together with the functional biases we show exist among single-copy marker genes, significantly changes the consensus evolutionary tree of Enterobacteriaceae.

Transitional genomes and nutritional role reversals identified for dual symbionts of adelgids (Aphidoidea: Adelgidae)

Many plant-sap-feeding insects have maintained a single, obligate, nutritional symbiont over the long history of their lineage. This senior symbiont may be joined by one or more junior symbionts that compensate for gaps in function incurred through genome-degradative forces. Adelgids are sap-sucking insects that feed solely on conifer trees and follow complex life cycles in which the diet fluctuates in nutrient levels. Adelgids are unusual in that both senior and junior symbionts appear to have been replaced repeatedly over their evolutionary history. Genomes can provide clues to understanding symbiont replacements, but only the dual symbionts of hemlock adelgids have been examined thus far. Here, we sequence and compare genomes of four additional dual-symbiont pairs in adelgids. We show that these symbionts are nutritional partners originating from diverse bacterial lineages and exhibiting wide variation in general genome characteristics. Although dual symbionts cooperate to produce nutrients, the balance of contributions varies widely across pairs, and total genome contents reflect a range of ages and degrees of degradation. Most symbionts appear to be in transitional states of genome reduction. Our findings support a hypothesis of periodic symbiont turnover driven by fluctuating selection for nutritional provisioning related to gains and losses of complex life cycles in their hosts.

Pervasiveness of the symbiont Serratia symbiotica in the aphid natural environment: distribution, diversity and evolution at a multitrophic level

Symbioses are significant drivers of insect evolutionary ecology. Despite recent findings that these associations can emerge from environmentally derived bacterial precursors, there is still little information on how these potential progenitors of insect symbionts circulate in trophic systems. Serratia symbiotica represents a valuable model for deciphering evolutionary scenarios of bacterial acquisition by insects, as its diversity includes gut-associated strains that retained the ability to live independently of their hosts, representing a potential reservoir for symbioses emergence. Here, we conducted a field study to examine the distribution and diversity of S. symbiotica found in aphid populations, and in different compartments of their surrounding environment. Twenty % of aphids colonies were infected with S. symbiotica, including a wide diversity of strains with varied tissue tropism corresponding to different lifestyle. We also showed that the prevalence of S. symbiotica is influenced by seasonal temperatures. We found that S. symbiotica was present in non-aphid species and in host plants, and that its prevalence in these samples was higher when associated aphid colonies were infected. Furthermore, phylogenetic analyses suggest the existence of horizontal transfers between the different trophic levels. These results provide a new picture of the pervasiveness of an insect symbiont in nature.

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.

The Genomics and Cell Biology of Host-Beneficial Intracellular Infections

Microbes gain access to eukaryotic cells as food for bacteria-grazing protists, for host protection by microbe-killing immune cells, or for microbial benefit when pathogens enter host cells to replicate. But microbes can also gain access to a host cell and become an important-often required-beneficial partner. The oldest beneficial microbial infections are the ancient eukaryotic organelles now called the mitochondrion and plastid. But numerous other host-beneficial intracellular infections occur throughout eukaryotes. Here I review the genomics and cell biology of these interactions with a focus on intracellular bacteria. The genomes of host-beneficial intracellular bacteria have features that span a previously unfilled gap between pathogens and organelles. Host cell adaptations to allow the intracellular persistence of beneficial bacteria are found along with evidence for the microbial manipulation of host cells, but the cellular mechanisms of beneficial bacterial infections are not well understood. Expected final online publication date for the Annual Review of Cell and Developmental Biology, Volume 37 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

At the Gate of Mutualism: Identification of Genomic Traits Predisposing to Insect-Bacterial Symbiosis in Pathogenic Strains of the Aphid Symbiont Serratia symbiotica

Mutualistic associations between insects and heritable bacterial symbionts are ubiquitous in nature. The aphid symbiont Serratia symbiotica is a valuable candidate for studying the evolution of bacterial symbiosis in insects because it includes a wide diversity of strains that reflect the diverse relationships in which bacteria can be engaged with insects, from pathogenic interactions to obligate intracellular mutualism. The recent discovery of culturable strains, which are hypothesized to resemble the ancestors of intracellular strains, provide an opportunity to study the mechanisms underlying bacterial symbiosis in its early stages. In this study, we analyzed the genomes of three of these culturable strains that are pathogenic to aphid hosts, and performed comparative genomic analyses including mutualistic host-dependent strains. All three genomes are larger than those of the host-restricted S. symbiotica strains described so far, and show significant enrichment in pseudogenes and mobile elements, suggesting that these three pathogenic strains are in the early stages of the adaptation to their host. Compared to their intracellular mutualistic relatives, the three strains harbor a greater diversity of genes coding for virulence factors and metabolic pathways, suggesting that they are likely adapted to infect new hosts and are a potential source of metabolic innovation for insects. The presence in their genomes of secondary metabolism gene clusters associated with the production of antimicrobial compounds and phytotoxins supports the hypothesis that S. symbiotia symbionts evolved from plant-associated strains and that plants may serve as intermediate hosts. Mutualistic associations between insects and bacteria are the result of independent transitions to endosymbiosis initiated by the acquisition of environmental progenitors. In this context, the genomes of free-living S. symbiotica strains provide a rare opportunity to study the inventory of genes held by bacterial associates of insects that are at the gateway to a host-dependent lifestyle.

Influential Insider: Wolbachia, an Intracellular Symbiont, Manipulates Bacterial Diversity in Its Insect Host

Facultative intracellular symbionts like the α-proteobacteria Wolbachia influence their insect host phenotype but little is known about how much they affect their host microbiota. Here, we quantified the impact of Wolbachia infection on the bacterial community of the cabbage root fly Delia radicum by comparing the microbiota of Wolbachia-free and infected adult flies of both sexes. We used high-throughput DNA sequencing (Illumina MiSeq, 16S rRNA, V5-V7 region) and performed a community and a network analysis. In both sexes, Wolbachia infection significantly decreased the diversity of D. radicum bacterial communities and modified their structure and composition by reducing abundance in some taxa but increasing it in others. Infection by Wolbachia was negatively correlated to 8 bacteria genera (Erwinia was the most impacted), and positively correlated to Providencia and Serratia. We suggest that Wolbachia might antagonize Erwinia for being entomopathogenic (and potentially intracellular), but would favor Providencia and Serratia because they might protect the host against chemical plant defenses. Although they might seem prisoners in a cell, endocellular symbionts can impact the whole microbiota of their host, hence its extended phenotype, which provides them with a way to interact with the outside world.

Multiple concurrent and convergent stages of genome reduction in bacterial symbionts across a stink bug family

Nutritional symbioses between bacteria and insects are prevalent and diverse, allowing insects to expand their feeding strategies and niches. A common consequence of long-term associations is a considerable reduction in symbiont genome size likely influenced by the radical shift in selective pressures as a result of the less variable environment within the host. While several of these cases can be found across distinct insect species, most examples provide a limited view of a single or few stages of the process of genome reduction. Stink bugs (Pentatomidae) contain inherited gamma-proteobacterial symbionts in a modified organ in their midgut and are an example of a long-term nutritional symbiosis, but multiple cases of new symbiont acquisition throughout the history of the family have been described. We sequenced the genomes of 11 symbionts of stink bugs with sizes that ranged from equal to those of their free-living relatives to less than 20%. Comparative genomics of these and previously sequenced symbionts revealed initial stages of genome reduction including an initial pseudogenization before genome reduction, followed by multiple stages of progressive degeneration of existing metabolic pathways likely to impact host interactions such as cell wall component biosynthesis. Amino acid biosynthesis pathways were retained in a similar manner as in other nutritional symbionts. Stink bug symbionts display convergent genome reduction events showing progressive changes from a free-living bacterium to a host-dependent symbiont. This system can therefore be used to study convergent genome evolution of symbiosis at a scale not previously available.

A Codon Model for Associating Phenotypic Traits with Altered Selective Patterns of Sequence Evolution

Detecting the signature of selection in coding sequences and associating it with shifts in phenotypic states can unveil genes underlying complex traits. Of the various signatures of selection exhibited at the molecular level, changes in the pattern of selection at protein-coding genes have been of main interest. To this end, phylogenetic branch-site codon models are routinely applied to detect changes in selective patterns along specific branches of the phylogeny. Many of these methods rely on a prespecified partition of the phylogeny to branch categories, thus treating the course of trait evolution as fully resolved and assuming that phenotypic transitions have occurred only at speciation events. Here, we present TraitRELAX, a new phylogenetic model that alleviates these strong assumptions by explicitly accounting for the uncertainty in the evolution of both trait and coding sequences. This joint statistical framework enables the detection of changes in selection intensity upon repeated trait transitions. We evaluated the performance of TraitRELAX using simulations and then applied it to two case studies. Using TraitRELAX, we found an intensification of selection in the primate SEMG2 gene in polygynandrous species compared to species of other mating forms, as well as changes in the intensity of purifying selection operating on sixteen bacterial genes upon transitioning from a free-living to an endosymbiotic lifestyle.[Evolutionary selection; intensification; $\gamma $-proteobacteria; genotype-phenotype; relaxation; SEMG2.].

Convergence of Nutritional Symbioses in Obligate Blood Feeders

Symbiosis with intracellular or gut bacteria is essential for the nutrition of animals with an obligate blood-feeding habit. Divergent bacterial lineages have independently evolved functional interactions with obligate blood feeders, but all converge to an analogous biochemical feature: the provisioning of B vitamins. Although symbionts and blood feeders coevolved interdependently for millions of years we stress that their associations are not necessarily stable. Ancestral symbionts can be replaced by recently acquired bacteria with similar biochemical features, a dynamic that emerges through a combination of phylogenetic and ecological constraints. Specifically, we highlight the lateral transfer of a streamlined biotin (B₇ vitamin) operon, and conjecture that its extensive spread across bacterial lineages may drive the emergence of novel nutritional symbioses with blood feeders.

Reduced Genome of the Gut Symbiotic Bacterium "Candidatus Benitsuchiphilus tojoi" Provides Insight Into Its Possible Roles in Ecology and Adaptation of the Host Insect

Diverse animals, including insects, harbor microbial symbionts within their gut, body cavity, or cells. The subsocial parastrachiid stinkbug Parastrachia japonensis is well-known for its peculiar ecological and behavioral traits, including its prolonged non-feeding diapause period and maternal care of eggs/nymphs in an underground nest. P. japonensis harbors a specific bacterial symbiont within the gut cavity extracellularly, which is vertically inherited through maternal excretion of symbiont-containing white mucus. Thus far, biological roles of the symbiont in the host lifecycle has been little understood. Here we sequenced the genome of the uncultivable gut symbiont “Candidatus Benitsuchiphilus tojoi.” The symbiont has an 804 kb circular chromosome encoding 606 proteins and a 14.5 kb plasmid encoding 13 proteins. Phylogenetic analysis indicated that the bacterium is closely related to other obligate insect symbionts belonging to the Gammaproteobacteria, including Buchnera of aphids and Blochmannia of ants, and the most closely related to Ishikawaella, an extracellular gut symbiont of plataspid stinkbugs. These data suggested that the symbiont genome has evolved like highly reduced gamma-proteobacterial symbiont genomes reported from a variety of insects. The presence of genes involved in biosynthesis pathways for amino acids, vitamins, and cofactors in the genome implicated the symbiont as a nutritional mutualist, supplementing essential nutrients to the host. Interestingly, the symbiont’s plasmid encoded genes for thiamine and carotenoid synthesis pathways, suggesting the possibility of additional functions of the symbiont for protecting the host against oxidative stress and DNA damage. Finally, possible involvement of the symbiont in uric acid metabolism during diapause is discussed.

Evolution and taxonomy of nematode-associated entomopathogenic bacteria of the genera Xenorhabdus and Photorhabdus: an overview

Entomopathogenic bacteria from the genera Photorhabdus and Xenorhabdus are closely related Gram-negative bacilli from the family Enterobacteriaceae (γ-Proteobacteria). They establish obligate mutualistic associations with soil nematodes from the genera Steinernema and Heterorhabditis to facilitate insect pathogenesis. The research of these two bacterial genera is focused mainly on their unique interactions with two different animal hosts, i.e. nematodes and insects. So far, studies of the mutualistic bacteria of nematodes collected from around the world have contributed to an increase in the number of the described Xenorhabdus and Photorhabdus species. Recently, the classification system of entomopatogenic nematode microsymbionts has undergone profound revision and now 26 species of the genus Xenorhabdus and 19 species of the genus Photorhabdus have been identified. Despite their similar life style and close phylogenetic origin, Photorhabdus and Xenorhabdus bacterial species differ significantly in e.g. the nematode host range, symbiotic strategies for parasite success, and arrays of released antibiotics and insecticidal toxins. As the knowledge of the diversity of entomopathogenic nematode microsymbionts helps to enable the use thereof, assessment of the phylogenetic relationships of these astounding bacterial genera is now a major challenge for researchers. The present article summarizes the main information on the taxonomy and evolutionary history of Xenorhabdus and Photorhabdus, entomopathogenic nematode symbionts.

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.

Legionella Becoming a Mutualist: Adaptive Processes Shaping the Genome of Symbiont in the Louse Polyplax serrata

Legionellaceae are intracellular bacteria known as important human pathogens. In the environment, they are mainly found in biofilms associated with amoebas. In contrast to the gammaproteobacterial family Enterobacteriaceae, which established a broad spectrum of symbioses with many insect taxa, the only instance of legionella-like symbiont has been reported from lice of the genus Polyplax. Here, we sequenced the complete genome of this symbiont and compared its main characteristics to other Legionella species and insect symbionts. Based on rigorous multigene phylogenetic analyses, we confirm this bacterium as a member of the genus Legionella and propose the name Candidatus Legionella polyplacis, sp.n. We show that the genome of Ca. Legionella polyplacis underwent massive degeneration, including considerable size reduction (529.746 bp, 484 protein coding genes) and a severe decrease in GC content (23%). We identify several possible constraints underlying the evolution of this bacterium. On one hand, Ca. Legionella polyplacis and the louse symbionts Riesia and Puchtella experienced convergent evolution, perhaps due to adaptation to similar hosts. On the other hand, some metabolic differences are likely to reflect different phylogenetic positions of the symbionts and hence availability of particular metabolic function in the ancestor. This is exemplified by different arrangements of thiamine metabolism in Ca. Legionella polyplacis and Riesia. Finally, horizontal gene transfer is shown to play a significant role in the adaptive and diversification process. Particularly, we show that Ca. L. polyplacis horizontally acquired a complete biotin operon (bioADCHFB) that likely assisted this bacterium when becoming an obligate mutualist.

Comparative genomics of a quadripartite symbiosis in a planthopper host reveals the origins and rearranged nutritional responsibilities of anciently diverged bacterial lineages

Insects in the Auchenorrhyncha (Hemiptera: Suborder) established nutritional symbioses with bacteria approximately 300 million years ago (MYA). The suborder split early during its diversification (~ 250 MYA) into the Fulgoroidea (planthoppers) and Cicadomorpha (leafhoppers and cicadas). The two lineages share some symbionts, including Sulcia and possibly a Betaproteobacteria that collaboratively provide their hosts with 10 essential amino acids (EAA). Some hosts harbour three bacteria, as is common among planthoppers. However, genomic studies are currently restricted to the dual-bacterial symbioses found in Cicadomorpha, leaving the origins and functions of these more complex symbioses unclear. To address these questions, we sequenced the genomes and performed phylogenomic analyses of 'Candidatus Sulcia muelleri' (Bacteroidetes), 'Ca. Vidania fulgoroideae' (Betaproteobacteria) and 'Ca. Purcelliella pentastirinorum' (Gammaproteobacteria) from a planthopper (Cixiidae: Oliarus). In contrast to the Cicadomorpha, nutritional synthesis responsibilities are rearranged between the cixiid symbionts. Although Sulcia has a highly conserved genome across the Auchenorrhyncha, in the cixiids it is greatly reduced and provides only three EAAs. Vidania contributes the remaining seven EAAs. Phylogenomic results suggest that it represents an ancient symbiont lineage paired with Sulcia throughout the Auchenorrhyncha. Finally, Purcelliella was recently acquired from plant-insect associated bacteria (Pantoea-Erwinia) to provide B vitamins and metabolic support to its degenerate partners.

Metagenomic analysis reveals changes of the Drosophila suzukii microbiota in the newly colonized regions

The spotted wing drosophila, Drosophila suzukii (Matsumura) (Diptera: Drosophilidae) is a highly polyphagous pest of a wide variety of wild or cultivated berry and stone fruit. Originating from Southeast Asia, it has recently invaded a wide range of regions in Europe and North America. It is well known that insect microbiotas may significantly influence several aspects of the host biology and play an important role in invasive species introduction into new areas. However, in spite of the great economic importance of D. suzukii, a limited attention has been given so far to its microbiota. In this study, we present the first in-depth characterization of gut bacterial diversity from field (native and invasive range) and lab-reared populations of this insect. The gut bacterial communities of field insects were dominated, regardless of their origin, by 2 families of the phylum Proteobacteria: Acetobacteraceae and Enterobacteriaceae, while Firmicutes, mainly represented by the family Staphylococcaceae, prevailed in lab-reared population. Locality was the most significant factor in shaping the microbiota of wild flies. Moreover, a negative correlation between diversity and abundance of Enterobacteriaceae and the time elapsed since the establishment of D. suzukii in a new region was observed. Altogether our results indicate that habitat, food resources as well as the colonization phase of a new region contribute to shape the bacterial communities of the invasive species which, in turn, by evolving more quickly, could influence host adaptation in a new environment.

The All-Rounder Sodalis: A New Bacteriome-Associated Endosymbiont of the Lygaeoid Bug Henestaris halophilus (Heteroptera: Henestarinae) and a Critical Examination of Its Evolution

Hemipteran insects are well-known in their ability to establish symbiotic relationships with bacteria. Among them, heteropteran insects present an array of symbiotic systems, ranging from the most common gut crypt symbiosis to the more restricted bacteriome-associated endosymbiosis, which have only been detected in members of the superfamily Lygaeoidea and the family Cimicidae so far. Genomic data of heteropteran endosymbionts are scarce and have merely been analyzed from the Wolbachia endosymbiont in bed bug and a few gut crypt-associated symbionts in pentatomoid bugs. In this study, we present the first detailed genomic analysis of a bacteriome-associated endosymbiont of a phytophagous heteropteran, present in the seed bug Henestaris halophilus (Hemiptera: Heteroptera: Lygaeoidea). Using phylogenomics and genomics approaches, we have assigned the newly characterized endosymbiont to the Sodalis genus, named as Candidatus Sodalis baculum sp. nov. strain kilmister. In addition, our findings support the reunification of the Sodalis genus, currently divided into six different genera. We have also conducted comparative analyses between 15 Sodalis species that present different genome sizes and symbiotic relationships. These analyses suggest that Ca. Sodalis baculum is a mutualistic endosymbiont capable of supplying the amino acids tyrosine, lysine, and some cofactors to its host. It has a small genome with pseudogenes but no mobile elements, which indicates middle-stage reductive evolution. Most of the genes in Ca. Sodalis baculum are likely to be evolving under purifying selection with several signals pointing to the retention of the lysine/tyrosine biosynthetic pathways compared with other Sodalis.

A Freeloader? The Highly Eroded Yet Large Genome of the Serratia symbiotica Symbiont of Cinara strobi

Genome reduction is pervasive among maternally inherited bacterial endosymbionts. This genome reduction can eventually lead to serious deterioration of essential metabolic pathways, thus rendering an obligate endosymbiont unable to provide essential nutrients to its host. This loss of essential pathways can lead to either symbiont complementation (sharing of the nutrient production with a novel co-obligate symbiont) or symbiont replacement (complete takeover of nutrient production by the novel symbiont). However, the process by which these two evolutionary events happen remains somewhat enigmatic by the lack of examples of intermediate stages of this process. Cinara aphids (Hemiptera: Aphididae) typically harbor two obligate bacterial symbionts: Buchnera and Serratia symbiotica. However, the latter has been replaced by different bacterial taxa in specific lineages, and thus species within this aphid lineage could provide important clues into the process of symbiont replacement. In the present study, using 16S rRNA high-throughput amplicon sequencing, we determined that the aphid Cinara strobi harbors not two, but three fixed bacterial symbionts: Buchnera aphidicola, a Sodalis sp., and S. symbiotica. Through genome assembly and genome-based metabolic inference, we have found that only the first two symbionts (Buchnera and Sodalis) actually contribute to the hosts' supply of essential nutrients while S. symbiotica has become unable to contribute towards this task. We found that S. symbiotica has a rather large and highly eroded genome which codes only for a few proteins and displays extensive pseudogenization. Thus, we propose an ongoing symbiont replacement within C. strobi, in which a once "competent" S. symbiotica does no longer contribute towards the beneficial association. These results suggest that in dual symbiotic systems, when a substitute cosymbiont is available, genome deterioration can precede genome reduction and a symbiont can be maintained despite the apparent lack of benefit to its host.

Comparative genomics and the nature of placozoan species

Placozoans are a phylum of nonbilaterian marine animals currently represented by a single described species, Trichoplax adhaerens, Schulze 1883. Placozoans arguably show the simplest animal morphology, which is identical among isolates collected worldwide, despite an apparently sizeable genetic diversity within the phylum. Here, we use a comparative genomics approach for a deeper appreciation of the structure and causes of the deeply diverging lineages in the Placozoa. We generated a high-quality draft genome of the genetic lineage H13 isolated from Hong Kong and compared it to the distantly related T. adhaerens. We uncovered substantial structural differences between the two genomes that point to a deep genomic separation and provide support that adaptation by gene duplication is likely a crucial mechanism in placozoan speciation. We further provide genetic evidence for reproductively isolated species and suggest a genus-level difference of H13 to T. adhaerens, justifying the designation of H13 as a new species, Hoilungia hongkongensis nov. gen., nov. spec., now the second described placozoan species and the first in a new genus. Our multilevel comparative genomics approach is, therefore, likely to prove valuable for species distinctions in other cryptic microscopic animal groups that lack diagnostic morphological characters, such as some nematodes, copepods, rotifers, or mites.

Partnering With a Pest: Genomes of Hemlock Woolly Adelgid Symbionts Reveal Atypical Nutritional Provisioning Patterns in Dual-Obligate Bacteria

Nutritional bacterial symbionts enhance the diets of sap-feeding insects with amino acids and vitamins missing from their diets. In many lineages, an ancestral senior symbiont is joined by a younger junior symbiont. To date, an emergent pattern is that senior symbionts supply a majority of amino acids, and junior symbionts supply a minority. Similar to other hemipterans, adelgids harbor obligate symbionts, but have higher diversity of bacterial associates, suggesting a history of symbiont turnover. The metabolic roles of dual symbionts in adelgids and their contributions to the consortium are largely unexplored. Here, we investigate the symbionts of Adelges tsugae, the hemlock woolly adelgid (HWA), an invasive species introduced from Japan to the eastern United States, where it kills hemlock trees. The response of hemlocks to HWA feeding has aspects of a defensive reaction against pathogens, and some have speculated that symbionts may be involved. We sequenced the genomes of "Ca. Annandia adelgestsuga" and "Ca. Pseudomonas adelgestsugas" symbionts to detail their metabolic capabilities, infer ages of relationship, and search for effectors of plant defenses. We also tested the relationship of "Ca. Annandia" to symbionts of other insects. We find that both symbionts provide nutrients, but in more balanced proportions than dual symbionts of other hemipterans. The lesser contributions of the senior "Ca. Annandia" support our hypothesis for symbiont replacements in adelgids. Phylogenomic results were ambiguous regarding the position of "Ca. Annandia". We found no obvious effectors of plant defenses related to insect virulence, but hypothetical proteins in symbionts are unknown players.

Compositional and mutational rate heterogeneity in mitochondrial genomes and its effect on the phylogenetic inferences of Cimicomorpha (Hemiptera: Heteroptera)

Background

Mitochondrial genome (mt-genome) data can potentially return artefactual relationships in the higher-level phylogenetic inference of insects due to the biases of accelerated substitution rates and compositional heterogeneity. Previous studies based on mt-genome data alone showed a paraphyly of Cimicomorpha (Insecta, Hemiptera) due to the positions of the families Tingidae and Reduviidae rather than the monophyly that was supported based on morphological characters, morphological and molecular combined data and large scale molecular datasets. Various strategies have been proposed to ameliorate the effects of potential mt-genome biases, including dense taxon sampling, removal of third codon positions or purine-pyrimidine coding and the use of site-heterogeneous models. In this study, we sequenced the mt-genomes of five additional Tingidae species and discussed the compositional and mutational rate heterogeneity in mt-genomes and its effect on the phylogenetic inferences of Cimicomorpha by implementing the bias-reduction strategies mentioned above.

Results

Heterogeneity in nucleotide composition and mutational biases were found in mt protein-coding genes, and the third codon exhibited high levels of saturation. Dense taxon sampling of Tingidae and Reduviidae and the other common strategies mentioned above were insufficient to recover the monophyly of the well-established group Cimicomorpha. When the sites with weak phylogenetic signals in the dataset were removed, the remaining dataset of mt-genomes can support the monophyly of Cimicomorpha; this support demonstrates that mt-genomes possess strong phylogenetic signals for the inference of higher-level phylogeny of this group. Comparison of the ratio of the removal of amino acids for each PCG showed that ATP8 has the highest ratio while CO1 has the lowest. This pattern is largely congruent with the evolutionary rate of 13 PCGs that ATP8 represents the highest evolutionary rate, whereas CO1 appears to be the lowest. Notably, the value of Ka/Ks ratios of all PCGs is less than 1, indicating that these genes are likely evolving under purifying selection.

Conclusions

Our results demonstrate that mt-genomes have sites with strong phylogenetic signals for the inference of higher-level phylogeny of Cimicomorpha. Consequently, bioinformatic approaches to removing sites with weak phylogenetic signals in mt-genome without relying on an a priori tree topology would greatly improve this field.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4650-9) contains supplementary material, which is available to authorized users.

Ancient symbiosis confers desiccation resistance to stored grain pest beetles

Microbial symbionts of insects provide a range of ecological traits to their hosts that are beneficial in the context of biotic interactions. However, little is known about insect symbiont-mediated adaptation to the abiotic environment, for example, temperature and humidity. Here, we report on an ancient clade of intracellular, bacteriome-located Bacteroidetes symbionts that are associated with grain and wood pest beetles of the phylogenetically distant families Silvanidae and Bostrichidae. In the saw-toothed grain beetle Oryzaephilus surinamensis, we demonstrate that the symbionts affect cuticle thickness, melanization and hydrocarbon profile, enhancing desiccation resistance and thereby strongly improving fitness under dry conditions. Together with earlier observations on symbiont contributions to cuticle biosynthesis in weevils, our findings indicate that convergent acquisitions of bacterial mutualists represented key adaptations enabling diverse pest beetle groups to survive and proliferate under the low ambient humidity that characterizes dry grain storage facilities.

Insect-bacteria parallel evolution in multiple-co-obligate-aphid association: a case in Lachninae (Hemiptera: Aphididae)

Parallel phylogenies between aphid and its obligate symbiont Buchnera are hot topics which always focused on aphid lower taxonomic levels. Symbionts in the subfamily Lachninae are special. Buchnera in many lachnine species has undergone functional and genome size reduction that was replaced by other co-obligate symbionts. In this study, we constructed the phylogenetic relationships of Lachninae with a combined dataset of five genes sequenced from Buchnera to estimate the effects of a dual symbiotic system in the aphid-Buchnera cospeciation association. The phylogeny of Buchnera in Lachninae was well-resolved in the combined dataset. Each of the genera formed strongly supported monophyletic groups, with the exception of the genus Cinara. The phylogeny based on sequences from Buchnera was divided into five tribes according to the clades of the Lachninae hosts tree, with the phylogenies of Buchnera and Lachninae being generally congruent. These results first provided evidence of parallel evolution at the aphid subfamily level comprehensively and supported the view that topological congruence between the phylogenies of Buchnera and Lachninae would not be interfered with the other co-obligate symbionts, such as Sarretia, in aphid-entosymbiont association. These results also provided new insight in understanding host-plant coevolution in lachnine lineages.

Effects of contaminants of emerging concern on Megaselia scalaris (Lowe, Diptera: Phoridae) and its microbial community

Drought, rising temperatures, and expanding human populations are increasing water demands. Many countries are extending potable water supplies by irrigating crops with wastewater. Unfortunately, wastewater contains biologically active, long-lived pharmaceuticals, even after treatment. Run-off from farms and wastewater treatment plant overflows contribute high concentrations of pharmaceuticals to the environment. This study assessed the effects of common pharmaceuticals on a cosmopolitan saprophagous insect, Megaselia scalaris (Diptera: Phoridae). Larvae were reared on artificial diets spiked with contaminants of emerging concern (CECs) at environmentally relevant concentrations. Female flies showed no oviposition preference for treated or untreated diets. Larvae exposed to caffeine in diets showed increased mortality, and larvae fed antibiotics and hormones showed signs of slowed development, especially in females. The normal sex ratio observed in M. scalaris from control diets was affected by exposure to caffeine and pharmaceutical mixture treatments. There was an overall effect of treatment on the flies’ microbial communities; notably, caffeine fed insects displayed higher microbial variability. Eight bacterial families accounted for approximately 95% of the total microbes in diet and insects. Our results suggest that CECs at environmentally relevant concentrations can affect the biology and microbial communities of an insect of ecological and medical importance.

Convergent patterns in the evolution of mealybug symbioses involving different intrabacterial symbionts

Mealybugs (Insecta: Hemiptera: Pseudococcidae) maintain obligatory relationships with bacterial symbionts, which provide essential nutrients to their insect hosts. Most pseudococcinae mealybugs harbor a unique symbiosis setup with enlarged betaproteobacterial symbionts (‘Candidatus Tremblaya princeps'), which themselves contain gammaproteobacterial symbionts. Here we investigated the symbiosis of the manna mealybug, Trabutina mannipara, using a metagenomic approach. Phylogenetic analyses revealed that the intrabacterial symbiont of T. mannipara represents a novel lineage within the Gammaproteobacteria, for which we propose the tentative name ‘Candidatus Trabutinella endobia'. Combining our results with previous data available for the nested symbiosis of the citrus mealybug Planococcus citri, we show that synthesis of essential amino acids and vitamins and translation-related functions partition between the symbiotic partners in a highly similar manner in the two systems, despite the distinct evolutionary origin of the intrabacterial symbionts. Bacterial genes found in both mealybug genomes and complementing missing functions in both symbioses were likely integrated in ancestral mealybugs before T. mannipara and P. citri diversified. The high level of correspondence between the two mealybug systems and their highly intertwined metabolic pathways are unprecedented. Our work contributes to a better understanding of the only known intracellular symbiosis between two bacteria and suggests that the evolution of this unique symbiosis included the replacement of intrabacterial symbionts in ancestral mealybugs.

Reinventing the Wheel and Making It Round Again: Evolutionary Convergence in Buchnera-Serratia Symbiotic Consortia between the Distantly Related Lachninae Aphids Tuberolachnus salignus and Cinara cedri

Virtually all aphids (Aphididae) harbor Buchnera aphidicola as an obligate endosymbiont to compensate nutritional deficiencies arising from their phloem diet. Many species within the Lachninae subfamily seem to be consistently associated also with Serratia symbiotica We have previously shown that both Cinara (Cinara) cedri and Cinara (Cupressobium) tujafilina (Lachninae: Eulachnini tribe) have indeed established co-obligate associations with both Buchnera and S. symbiotica However, while Buchnera genomes of both Cinara species are similar, genome degradation differs greatly between the two S. symbiotica strains. To gain insight into the essentiality and degree of integration of S. symbiotica within the Lachninae, we sequenced the genome of both Buchnera and S. symbiotica endosymbionts from the distantly related aphid Tuberolachnus salignus (Lachninae: Tuberolachnini tribe). We found a striking level of similarity between the endosymbiotic system of this aphid and that of C. cedri In both aphid hosts, S. symbiotica possesses a highly reduced genome and is found exclusively intracellularly inside bacteriocytes. Interestingly, T. salignus' endosymbionts present the same tryptophan biosynthetic metabolic complementation as C. cedri's, which is not present in C. tujafilina's. Moreover, we corroborate the riboflavin-biosynthetic-role take-over/rescue by S. symbiotica in T. salignus, and therefore, provide further evidence for the previously proposed establishment of a secondary co-obligate endosymbiont in the common ancestor of the Lachninae aphids. Finally, we propose that the putative convergent split of the tryptophan biosynthetic role between Buchnera and S. symbiotica could be behind the establishment of S. symbiotica as an obligate intracellular symbiont and the triggering of further genome degradation.

Capturing the Phylogeny of Holometabola with Mitochondrial Genome Data and Bayesian Site-Heterogeneous Mixture Models

After decades of debate, a mostly satisfactory resolution of relationships among the 11 recognized holometabolan orders of insects has been reached based on nuclear genes, resolving one of the most substantial branches of the tree-of-life, but the relationships are still not well established with mitochondrial genome data. The main reasons have been the absence of sufficient data in several orders and lack of appropriate phylogenetic methods that avoid the systematic errors from compositional and mutational biases in insect mitochondrial genomes. In this study, we assembled the richest taxon sampling of Holometabola to date (199 species in 11 orders), and analyzed both nucleotide and amino acid data sets using several methods. We find the standard Bayesian inference and maximum-likelihood analyses were strongly affected by systematic biases, but the site-heterogeneous mixture model implemented in PhyloBayes avoided the false grouping of unrelated taxa exhibiting similar base composition and accelerated evolutionary rate. The inclusion of rRNA genes and removal of fast-evolving sites with the observed variability sorting method for identifying sites deviating from the mean rates improved the phylogenetic inferences under a site-heterogeneous model, correctly recovering most deep branches of the Holometabola phylogeny. We suggest that the use of mitochondrial genome data for resolving deep phylogenetic relationships requires an assessment of the potential impact of substitutional saturation and compositional biases through data deletion strategies and by using site-heterogeneous mixture models. Our study suggests a practical approach for how to use densely sampled mitochondrial genome data in phylogenetic analyses.

Loss of Conserved Noncoding RNAs in Genomes of Bacterial Endosymbionts

The genomes of intracellular symbiotic or pathogenic bacteria, such as of Buchnera, Mycoplasma, and Rickettsia, are typically smaller compared with their free-living counterparts. Here we showed that noncoding RNA (ncRNA) families, which are conserved in free-living bacteria, frequently could not be detected by computational methods in the small genomes. Statistical tests demonstrated that their absence is not an artifact of low GC content or small deletions in these small genomes, and thus it was indicative of an independent loss of ncRNAs in different endosymbiotic lineages. By analyzing the synteny (conservation of gene order) between the reduced and nonreduced genomes, we revealed instances of protein-coding genes that were preserved in the reduced genomes but lost cis-regulatory elements. We found that the loss of cis-regulatory ncRNA sequences, which regulate the expression of cognate protein-coding genes, is characterized by the reduction of secondary structure formation propensity, GC content, and length of the corresponding genomic regions.

The Distant Siblings-A Phylogenomic Roadmap Illuminates the Origins of Extant Diversity in Fungal Aromatic Polyketide Biosynthesis

In recent years, the influx of newly sequenced fungal genomes has enabled sampling of secondary metabolite biosynthesis on an unprecedented scale. However, explanations of extant diversity which take into account both large-scale phylogeny reconstructions and knowledge gained from multiple genome projects are still lacking. We analyzed the evolutionary sources of genetic diversity in aromatic polyketide biosynthesis in over 100 model fungal genomes. By reconciling the history of over 400 nonreducing polyketide synthases (NR-PKSs) with corresponding species history, we demonstrate that extant fungal NR-PKSs are clades of distant siblings, originating from a burst of duplications in early Pezizomycotina and thinned by extensive losses. The capability of higher fungi to biosynthesize the simplest precursor molecule (orsellinic acid) is highlighted as an ancestral trait underlying biosynthesis of aromatic compounds. This base activity was modified during early evolution of filamentous fungi, toward divergent reaction schemes associated with biosynthesis of, for example, aflatoxins and fusarubins (C4–C9 cyclization) or various anthraquinone derivatives (C6–C11 cyclization). The functional plasticity is further shown to have been supplemented by modularization of domain architecture into discrete pieces (conserved splice junctions within product template domain), as well as tight linkage of key accessory enzyme families and divergence in employed transcriptional factors. Although the majority of discord between species and gene history is explained by ancient duplications, this landscape has been altered by more recent duplications, as well as multiple horizontal gene transfers. The 25 detected transfers include previously undescribed events leading to emergence of, for example, fusarubin biosynthesis in Fusarium genus. Both the underlying data and the results of present analysis (including alternative scenarios revealed by sampling multiple reconciliation optima) are maintained as a freely available web-based resource: http://cropnet.pl/metasites/sekmet/nrpks_2014.

Solving a Bloody Mess: B-Vitamin Independent Metabolic Convergence among Gammaproteobacterial Obligate Endosymbionts from Blood-Feeding Arthropods and the Leech Haementeria officinalis

Endosymbiosis is a common phenomenon in nature, especially between bacteria and insects, whose typically unbalanced diets are usually complemented by their obligate endosymbionts. While much interest and focus has been directed toward phloem-feeders like aphids and mealybugs, blood-feeders such as the Lone star tick (Amblyomma americanum), Glossina flies, and the human body louse (Pediculus humanus corporis) depend on obligate endosymbionts which complement their B-vitamin-deficient diets, and thus are required for growth and survival. Glossiphoniid leeches have also been found to harbor distinct endosymbionts housed in specialized organs. Here, we present the genome of the bacterial endosymbiont from Haementeria officinalis, first of a glossiphoniid leech. This as-yet-unnamed endosymbiont belongs to the Gammaproteobacteria, has a pleomorphic shape and is restricted to bacteriocytes. For this bacterial endosymbiont, we propose the name Candidatus Providencia siddallii. This symbiont possesses a highly reduced genome with high A+T content and a reduced set of metabolic capabilities, all of which are common characteristics of ancient obligate endosymbionts of arthropods. Its genome has retained many pathways related to the biosynthesis of B-vitamins, pointing toward a role in supplementing the blood-restricted diet of its host. Through comparative genomics against the endosymbionts of A. americanum, Glossina flies, and P. humanus corporis, we were able to detect a high degree of metabolic convergence among these four very distantly related endosymbiotic bacteria.

Avoidance and Potential Remedy Solutions of Chimeras in Reconstructing the Phylogeny of Aphids Using the 16S rRNA Gene of Buchnera: A Case in Lachninae (Hemiptera)

It is known that PCR amplification of highly homologous genes from complex DNA mixtures can generate a significant proportion of chimeric sequences. The 16S rRNA gene is not only widely used in estimating the species diversity of endosymbionts in aphids but also used to explore the co-diversification of aphids and their endosymbionts. Thus, chimeric sequences may lead to the discovery of non-existent endosymbiont species and mislead Buchnera-based phylogenetic analysis that lead to false conclusions. In this study, a high probability (6.49%) of chimeric sequence occurrence was found in the amplified 16S rRNA gene sequences of endosymbionts from aphid species in the subfamily Lachninae. These chimeras are hybrid products of multiple parent sequences from the dominant species of endosymbionts in each corresponding host. It is difficult to identify the chimeric sequences of a new or unidentified species due to the high variability of their main parent, Buchnera aphidicola, and because the chimeric sequences can confuse the phylogenetic analysis of 16S rRNA gene sequences. These chimeras present a challenge to Buchnera-based phylogenetic research in aphids. Thus, our study strongly suggests that using appropriate methods to detect chimeric 16S rRNA sequences may avoid some false conclusions in endosymbiont-based aphid research.

Draft Genomes, Phylogenetic Reconstruction, and Comparative Genomics of Two Novel Cohabiting Bacterial Symbionts Isolated from Frankliniella occidentalis

Obligate bacterial symbionts are widespread in many invertebrates, where they are often confined to specialized host cells and are transmitted directly from mother to progeny. Increasing numbers of these bacteria are being characterized but questions remain about their population structure and evolution. Here we take a comparative genomics approach to investigate two prominent bacterial symbionts (BFo1 and BFo2) isolated from geographically separated populations of western flower thrips, Frankliniella occidentalis. Our multifaceted approach to classifying these symbionts includes concatenated multilocus sequence analysis (MLSA) phylogenies, ribosomal multilocus sequence typing (rMLST), construction of whole-genome phylogenies, and in-depth genomic comparisons. We showed that the BFo1 genome clusters more closely to species in the genus Erwinia, and is a putative close relative to Erwinia aphidicola. BFo1 is also likely to have shared a common ancestor with Erwinia pyrifoliae/Erwinia amylovora and the nonpathogenic Erwinia tasmaniensis and genetic traits similar to Erwinia billingiae. The BFo1 genome contained virulence factors found in the genus Erwinia but represented a divergent lineage. In contrast, we showed that BFo2 belongs within the Enterobacteriales but does not group closely with any currently known bacterial species. Concatenated MLSA phylogenies indicate that it may have shared a common ancestor to the Erwinia and Pantoea genera, and based on the clustering of rMLST genes, it was most closely related to Pantoea ananatis but represented a divergent lineage. We reconstructed a core genome of a putative common ancestor of Erwinia and Pantoea and compared this with the genomes of BFo bacteria. BFo2 possessed none of the virulence determinants that were omnipresent in the Erwinia and Pantoea genera. Taken together, these data are consistent with BFo2 representing a highly novel species that maybe related to known Pantoea.

A novel intracellular mutualistic bacterium in the invasive ant Cardiocondyla obscurior

The evolution of eukaryotic organisms is often strongly influenced by microbial symbionts that confer novel traits to their hosts. Here we describe the intracellular Enterobacteriaceae symbiont of the invasive ant Cardiocondyla obscurior, 'Candidatus Westeberhardia cardiocondylae'. Upon metamorphosis, Westeberhardia is found in gut-associated bacteriomes that deteriorate following eclosion. Only queens maintain Westeberhardia in the ovarian nurse cells from where the symbionts are transmitted to late-stage oocytes during nurse cell depletion. Functional analyses of the streamlined genome of Westeberhardia (533 kb, 23.41% GC content) indicate that neither vitamins nor essential amino acids are provided for the host. However, the genome encodes for an almost complete shikimate pathway leading to 4-hydroxyphenylpyruvate, which could be converted into tyrosine by the host. Taken together with increasing titers of Westeberhardia during pupal stage, this suggests a contribution of Westeberhardia to cuticle formation. Despite a widespread occurrence of Westeberhardia across host populations, one ant lineage was found to be naturally symbiont-free, pointing to the loss of an otherwise prevalent endosymbiont. This study yields insights into a novel intracellular mutualist that could play a role in the invasive success of C. obscurior.

Candidatus Frankia Datiscae Dg1, the Actinobacterial Microsymbiont of Datisca glomerata, Expresses the Canonical nod Genes nodABC in Symbiosis with Its Host Plant

Frankia strains are nitrogen-fixing soil actinobacteria that can form root symbioses with actinorhizal plants. Phylogenetically, symbiotic frankiae can be divided into three clusters, and this division also corresponds to host specificity groups. The strains of cluster II which form symbioses with actinorhizal Rosales and Cucurbitales, thus displaying a broad host range, show suprisingly low genetic diversity and to date can not be cultured. The genome of the first representative of this cluster, Candidatus Frankia datiscae Dg1 (Dg1), a microsymbiont of Datisca glomerata, was recently sequenced. A phylogenetic analysis of 50 different housekeeping genes of Dg1 and three published Frankia genomes showed that cluster II is basal among the symbiotic Frankia clusters. Detailed analysis showed that nodules of D. glomerata, independent of the origin of the inoculum, contain several closely related cluster II Frankia operational taxonomic units. Actinorhizal plants and legumes both belong to the nitrogen-fixing plant clade, and bacterial signaling in both groups involves the common symbiotic pathway also used by arbuscular mycorrhizal fungi. However, so far, no molecules resembling rhizobial Nod factors could be isolated from Frankia cultures. Alone among Frankia genomes available to date, the genome of Dg1 contains the canonical nod genes nodA, nodB and nodC known from rhizobia, and these genes are arranged in two operons which are expressed in D. glomerata nodules. Furthermore, Frankia Dg1 nodC was able to partially complement a Rhizobium leguminosarum A34 nodC::Tn5 mutant. Phylogenetic analysis showed that Dg1 Nod proteins are positioned at the root of both α- and β-rhizobial NodABC proteins. NodA-like acyl transferases were found across the phylum Actinobacteria, but among Proteobacteria only in nodulators. Taken together, our evidence indicates an Actinobacterial origin of rhizobial Nod factors.

The Recent Evolution of a Maternally-Inherited Endosymbiont of Ticks Led to the Emergence of the Q Fever Pathogen, Coxiella burnetii

Q fever is a highly infectious disease with a worldwide distribution. Its causative agent, the intracellular bacterium Coxiella burnetii, infects a variety of vertebrate species, including humans. Its evolutionary origin remains almost entirely unknown and uncertainty persists regarding the identity and lifestyle of its ancestors. A few tick species were recently found to harbor maternally-inherited Coxiella-like organisms engaged in symbiotic interactions, but their relationships to the Q fever pathogen remain unclear. Here, we extensively sampled ticks, identifying new and atypical Coxiella strains from 40 of 58 examined species, and used this data to infer the evolutionary processes leading to the emergence of C. burnetii. Phylogenetic analyses of multi-locus typing and whole-genome sequencing data revealed that Coxiella-like organisms represent an ancient and monophyletic group allied to ticks. Remarkably, all known C. burnetii strains originate within this group and are the descendants of a Coxiella-like progenitor hosted by ticks. Using both colony-reared and field-collected gravid females, we further establish the presence of highly efficient maternal transmission of these Coxiella-like organisms in four examined tick species, a pattern coherent with an endosymbiotic lifestyle. Our laboratory culture assays also showed that these Coxiella-like organisms were not amenable to culture in the vertebrate cell environment, suggesting different metabolic requirements compared to C. burnetii. Altogether, this corpus of data demonstrates that C. burnetii recently evolved from an inherited symbiont of ticks which succeeded in infecting vertebrate cells, likely by the acquisition of novel virulence factors.

How virulent infectious diseases emerge from non-pathogenic organisms is a challenging question. Here, we address this evolutionary issue in the case of Q fever. Its causative agent, the intracellular bacterium Coxiella burnetii, is extremely infectious to humans and a variety of animals. However, uncertainty persists regarding its evolutionary origin, including the identity and lifestyle of its ancestors. In this article, we show that C. burnetii arose from a rare evolutionary transformation of a maternally-inherited endosymbiont of ticks into a specialized and virulent pathogen of vertebrates. While arthropod symbionts are typically transmitted maternally and thought not to be infectious to vertebrates, we establish here that one Coxiella symbiont has evolved the necessary adaptations to exploit the vertebrate cell, leading to the emergence of Q fever.

Endosymbiont evolution: predictions from theory and surprises from genomes

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

Genome evolution in an ancient bacteria-ant symbiosis: parallel gene loss among Blochmannia spanning the origin of the ant tribe Camponotini

Stable associations between bacterial endosymbionts and insect hosts provide opportunities to explore genome evolution in the context of established mutualisms and assess the roles of selection and genetic drift across host lineages and habitats. Blochmannia, obligate endosymbionts of ants of the tribe Camponotini, have coevolved with their ant hosts for ∼40 MY. To investigate early events in Blochmannia genome evolution across this ant host tribe, we sequenced Blochmannia from two divergent host lineages, Colobopsis obliquus and Polyrhachis turneri, and compared them with four published genomes from Blochmannia of Camponotus sensu stricto. Reconstructed gene content of the last common ancestor (LCA) of these six Blochmannia genomes is reduced (690 protein coding genes), consistent with rapid gene loss soon after establishment of the symbiosis. Differential gene loss among Blochmannia lineages has affected cellular functions and metabolic pathways, including DNA replication and repair, vitamin biosynthesis and membrane proteins. Blochmannia of P. turneri (i.e., B. turneri) encodes an intact DnaA chromosomal replication initiation protein, demonstrating that loss of dnaA was not essential for establishment of the symbiosis. Based on gene content, B. obliquus and B. turneri are unable to provision hosts with riboflavin. Of the six sequenced Blochmannia, B. obliquus is the earliest diverging lineage (i.e., the sister group of other Blochmannia sampled) and encodes the fewest protein-coding genes and the most pseudogenes. We identified 55 genes involved in parallel gene loss, including glutamine synthetase, which may participate in nitrogen recycling. Pathways for biosynthesis of coenzyme A, terpenoids and riboflavin were lost in multiple lineages, suggesting relaxed selection on the pathway after inactivation of one component. Analysis of Illumina read datasets did not detect evidence of plasmids encoding missing functions, nor the presence of coresident symbionts other than Wolbachia. Although gene order is strictly conserved in four Blochmannia of Camponotus sensu stricto, comparisons with deeply divergent lineages revealed inversions in eight genomic regions, indicating ongoing recombination despite ancestral loss of recA. In sum, the addition of two Blochmannia genomes of divergent host lineages enables reconstruction of early events in evolution of this symbiosis and suggests that Blochmannia lineages may experience distinct, host-associated selective pressures. Understanding how evolutionary forces shape genome reduction in this system may help to clarify forces driving gene loss in other bacteria, including intracellular pathogens.

Comparative genomic and phylogenetic analyses of Gammaproteobacterial glg genes traced the origin of the Escherichia coli glycogen glgBXCAP operon to the last common ancestor of the sister orders Enterobacteriales and Pasteurellales

Production of branched α-glucan, glycogen-like polymers is widely spread in the Bacteria domain. The glycogen pathway of synthesis and degradation has been fairly well characterized in the model enterobacterial species Escherichia coli (order Enterobacteriales, class Gammaproteobacteria), in which the cognate genes (branching enzyme glgB, debranching enzyme glgX, ADP-glucose pyrophosphorylase glgC, glycogen synthase glgA, and glycogen phosphorylase glgP) are clustered in a glgBXCAP operon arrangement. However, the evolutionary origin of this particular arrangement and of its constituent genes is unknown. Here, by using 265 complete gammaproteobacterial genomes we have carried out a comparative analysis of the presence, copy number and arrangement of glg genes in all lineages of the Gammaproteobacteria. These analyses revealed large variations in glg gene presence, copy number and arrangements among different gammaproteobacterial lineages. However, the glgBXCAP arrangement was remarkably conserved in all glg-possessing species of the orders Enterobacteriales and Pasteurellales (the E/P group). Subsequent phylogenetic analyses of glg genes present in the Gammaproteobacteria and in other main bacterial groups indicated that glg genes have undergone a complex evolutionary history in which horizontal gene transfer may have played an important role. These analyses also revealed that the E/P glgBXCAP genes (a) share a common evolutionary origin, (b) were vertically transmitted within the E/P group, and (c) are closely related to glg genes of some phylogenetically distant betaproteobacterial species. The overall data allowed tracing the origin of the E. coli glgBXCAP operon to the last common ancestor of the E/P group, and also to uncover a likely glgBXCAP transfer event from the E/P group to particular lineages of the Betaproteobacteria.

RELAX: detecting relaxed selection in a phylogenetic framework

Relaxation of selective strength, manifested as a reduction in the efficiency or intensity of natural selection, can drive evolutionary innovation and presage lineage extinction or loss of function. Mechanisms through which selection can be relaxed range from the removal of an existing selective constraint to a reduction in effective population size. Standard methods for estimating the strength and extent of purifying or positive selection from molecular sequence data are not suitable for detecting relaxed selection, because they lack power and can mistake an increase in the intensity of positive selection for relaxation of both purifying and positive selection. Here, we present a general hypothesis testing framework (RELAX) for detecting relaxed selection in a codon-based phylogenetic framework. Given two subsets of branches in a phylogeny, RELAX can determine whether selective strength was relaxed or intensified in one of these subsets relative to the other. We establish the validity of our test via simulations and show that it can distinguish between increased positive selection and a relaxation of selective strength. We also demonstrate the power of RELAX in a variety of biological scenarios where relaxation of selection has been hypothesized or demonstrated previously. We find that obligate and facultative γ-proteobacteria endosymbionts of insects are under relaxed selection compared with their free-living relatives and obligate endosymbionts are under relaxed selection compared with facultative endosymbionts. Selective strength is also relaxed in asexual Daphnia pulex lineages, compared with sexual lineages. Endogenous, nonfunctional, bornavirus-like elements are found to be under relaxed selection compared with exogenous Borna viruses. Finally, selection on the short-wavelength sensitive, SWS1, opsin genes in echolocating and nonecholocating bats is relaxed only in lineages in which this gene underwent pseudogenization; however, selection on the functional medium/long-wavelength sensitive opsin, M/LWS1, is found to be relaxed in all echolocating bats compared with nonecholocating bats.

Phylogenomic analyses uncover origin and spread of the Wolbachia pandemic

Of all obligate intracellular bacteria, Wolbachia is probably the most common. In general, Wolbachia are either widespread, opportunistic reproductive parasites of arthropods or essential mutualists in a single group of filarial nematodes, including many species of medical significance. To date, a robust phylogenetic backbone of Wolbachia is lacking and consequently, many Wolbachia-related phenomena cannot be discussed in a broader evolutionary context. Here we present the first comprehensive phylogenomic analysis of Wolbachia supergroup relationships based on new whole-genome-shotgun data. Our results suggest that Wolbachia has switched between its two major host groups at least twice. The ability of some arthropod-infecting Wolbachia to universally infect and to adapt to a broad range of hosts quickly is restricted to a single monophyletic lineage (containing supergroups A and B). Thus, the currently observable pandemic has likely a single evolutionary origin and is unique within the radiation of Wolbachia strains.