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

Accelerated Pseudogenization in the Ancient Endosymbionts of Giant Scale Insects

Symbiotic microorganisms are subject to a complex interplay of environmental and population-genetic pressures that drive their gene loss. Despite the widely held perception that ancient symbionts have stable genomes, even tiny genomes experience ongoing pseudogenization. Whether these tiny genomes also experience bursts of rapid gene loss is, however, less understood. Giant scale insects (Monophlebidae) feed on plant sap and rely on the symbiotic bacterium Walczuchella, which provides them with essential nutrients. When compared with other ancient symbionts with similar genome sizes, such as Karelsulcia, Walczuchella's genome was previously reported as unusually pseudogene-rich (10% of coding sequences). However, this result was based on only one genome assembly, raising questions about the assembly quality or a recent ecological shift such as co-symbiont acquisition driving the gene loss. Here, we generated six complete genomes of Walczuchella from three genera of giant scales, each with distinct co-symbiotic partners. We show that all the genomes are highly degraded, and particularly genes related to the cellular envelope and energy metabolism seem to be undergoing pseudogenization. Apart from general mechanisms driving genome reduction, such as the long-term intracellular lifestyle with transmission bottlenecks, we hypothesize that a more profound loss of DNA replication and repair genes, together with recent co-obligate symbiont acquisitions, likely contribute to the accelerated degradation of Walczuchella genomes. Our results highlight that even ancient symbionts with small genomes can experience significant bursts of gene loss when stochastic processes erase a gene that accelerates gene loss or when the selection pressure changes such as after co-symbiont acquisition.

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.

Chromosome-level genome assembly of the aster leafhopper (Macrosteles quadrilineatus) reveals the role of environment and microbial symbiosis in shaping pest insect genome evolution

Leafhoppers comprise over 20,000 plant-sap feeding species, many of which are important agricultural pests. Most species rely on two ancestral bacterial symbionts, Sulcia and Nasuia, for essential nutrition lacking in their phloem and xylem plant sap diets. To understand how pest leafhopper genomes evolve and are shaped by microbial symbioses, we completed a chromosomal-level assembly of the aster leafhopper's genome (ALF; Macrosteles quadrilineatus). We compared ALF's genome to three other pest leafhoppers, Nephotettix cincticeps, Homalodisca vitripennis, and Empoasca onukii, which have distinct ecologies and symbiotic relationships. Despite diverging ~155 million years ago, leafhoppers have high levels of chromosomal synteny and gene family conservation. Conserved genes include those involved in plant chemical detoxification, resistance to various insecticides, and defence against environmental stress. Positive selection acting upon these genes further points to ongoing adaptive evolution in response to agricultural environments. In relation to leafhoppers' general dependence on symbionts, species that retain the ancestral symbiont, Sulcia, displayed gene enrichment of metabolic processes in their genomes. Leafhoppers with both Sulcia and its ancient partner, Nasuia, showed genomic enrichment in genes related to microbial population regulation and immune responses. Finally, horizontally transferred genes (HTGs) associated with symbiont support of Sulcia and Nasuia are only observed in leafhoppers that maintain symbionts. In contrast, HTGs involved in non-symbiotic functions are conserved across all species. The high-quality ALF genome provides deep insights into how host ecology and symbioses shape genome evolution and a wealth of genetic resources for pest control targets.

Insect Bacteriocytes: Adaptation, Development, and Evolution

Bacteriocytes are host cells specialized to harbor symbionts in certain insect taxa. The adaptation, development, and evolution of bacteriocytes underlie insect symbiosis maintenance. Bacteriocytes carry enriched host genes of insect and bacterial origin whose transcription can be regulated by microRNAs, which are involved in host-symbiont metabolic interactions. Recognition proteins of peptidoglycan, the bacterial cell wall component, and autophagy regulate symbiont abundance in bacteriocytes. Horizontally transferred genes expressed in bacteriocytes influence the metabolism of symbiont peptidoglycan, which may affect the bacteriocyte immune response against symbionts. Bacteriocytes release or transport symbionts into ovaries for symbiont vertical transmission. Bacteriocyte development and death, regulated by transcriptional factors, are variable in different insect species. The evolutionary origin of insect bacteriocytes remains unclear. Future research should elucidate bacteriocyte cell biology, the molecular interplay between bacteriocyte metabolic and immune functions, the genetic basis of bacteriocyte origin, and the coordination between bacteriocyte function and host biology in diverse symbioses.

Nutrient supplementation by genome-eroded Burkholderia symbionts of scale insects

Hemipterans are known as hosts to bacterial or fungal symbionts that supplement their unbalanced diet with essential nutrients. Among them, scale insects (Coccomorpha) are characterized by a particularly large diversity of symbiotic systems. Here, using microscopic and genomic approaches, we functionally characterized the symbionts of two scale insects belonging to the Eriococcidae family, Acanthococcus aceris and Gossyparia spuria. These species host Burkholderia bacteria that are localized in the cytoplasm of the fat body cells. Metagenome sequencing revealed very similar and highly reduced genomes (<900KBp) with a low GC content (~38%), making them the smallest and most AT-biased Burkholderia genomes yet sequenced. In their eroded genomes, both symbionts retain biosynthetic pathways for the essential amino acids leucine, isoleucine, valine, threonine, lysine, arginine, histidine, phenylalanine, and precursors for the semi-essential amino acid tyrosine, as well as the cobalamin-dependent methionine synthase MetH. A tryptophan biosynthesis pathway is conserved in the symbiont of G. spuria, but appeared pseudogenized in A. aceris, suggesting differential availability of tryptophan in the two host species' diets. In addition to the pathways for essential amino acid biosynthesis, both symbionts maintain biosynthetic pathways for multiple cofactors, including riboflavin, cobalamin, thiamine, and folate. The localization of Burkholderia symbionts and their genome traits indicate that the symbiosis between Burkholderia and eriococcids is younger than other hemipteran symbioses, but is functionally convergent. Our results add to the emerging picture of dynamic symbiont replacements in sap-sucking Hemiptera and highlight Burkholderia as widespread and versatile intra- and extracellular symbionts of animals, plants, and fungi.

A carbon-nitrogen negative feedback loop underlies the repeated evolution of cnidarian-Symbiodiniaceae symbioses

Symbiotic associations with Symbiodiniaceae have evolved independently across a diverse range of cnidarian taxa including reef-building corals, sea anemones, and jellyfish, yet the molecular mechanisms underlying their regulation and repeated evolution are still elusive. Here, we show that despite their independent evolution, cnidarian hosts use the same carbon-nitrogen negative feedback loop to control symbiont proliferation. Symbiont-derived photosynthates are used to assimilate nitrogenous waste via glutamine synthetase-glutamate synthase-mediated amino acid biosynthesis in a carbon-dependent manner, which regulates the availability of nitrogen to the symbionts. Using nutrient supplementation experiments, we show that the provision of additional carbohydrates significantly reduces symbiont density while ammonium promotes symbiont proliferation. High-resolution metabolic analysis confirmed that all hosts co-incorporated glucose-derived 13C and ammonium-derived 15N via glutamine synthetase-glutamate synthase-mediated amino acid biosynthesis. Our results reveal a general carbon-nitrogen negative feedback loop underlying these symbioses and provide a parsimonious explanation for their repeated evolution.

Unveiling metabolic integration in psyllids and their nutritional endosymbionts through comparative transcriptomics analysis

Psyllids, a group of insects that feed on plant sap, have a symbiotic relationship with an endosymbiont called Carsonella. Carsonella synthesizes essential amino acids and vitamins for its psyllid host, but lacks certain genes required for this process, suggesting a compensatory role of psyllid host genes. To investigate this, gene expression was compared between two psyllid species, Bactericera cockerelli and Diaphorina citri, in specialized cells where Carsonella resides (bacteriomes). Collaborative psyllid genes, including horizontally transferred genes, showed patterns of conserved gene expression; however, species-specific patterns were also observed, suggesting differences in the nutritional metabolism between psyllid species. Also, the recycling of nitrogen in bacteriomes may primarily rely on glutamate dehydrogenase (GDH). Additionally, lineage-specific gene clusters were differentially expressed in B. cockerelli and D. citri bacteriomes and are highlighted here. These findings shed light on potential host adaptations for the regulation of this symbiosis due to host, microbiome, and environmental differences.

Complete genomes of mutualistic bacterial co-symbionts "Candidatus Sulcia muelleri" and "Candidatus Nasuia deltocephalinicola" of the rice green leafhopper Nephotettix cincticeps

The genomes of obligate bacterial co-symbionts of the green rice leafhopper Nephotettix cincticeps, which is notorious as an agricultural pest, were determined. The streamlined genomes of "Candidatus Sulcia muelleri" and "Candidatus Nasuia deltocephalinicola" exhibited complementary metabolic pathways for synthesizing essential nutrients that contribute to host adaptation.

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.

Variable organization of symbiont-containing tissue across planthoppers hosting different heritable endosymbionts

Sap-feeding hemipteran insects live in associations with diverse heritable symbiotic microorganisms (bacteria and fungi) that provide essential nutrients deficient in their hosts' diets. These symbionts typically reside in highly specialized organs called bacteriomes (with bacterial symbionts) or mycetomes (with fungal symbionts). The organization of these organs varies between insect clades that are ancestrally associated with different microbes. As these symbioses evolve and additional microorganisms complement or replace the ancient associates, the organization of the symbiont-containing tissue becomes even more variable. Planthoppers (Hemiptera: Fulgoromorpha) are ancestrally associated with bacterial symbionts Sulcia and Vidania, but in many of the planthopper lineages, these symbionts are now accompanied or have been replaced by other heritable bacteria (e.g., Sodalis, Arsenophonus, Purcelliella) or fungi. We know the identity of many of these microbes, but the symbiont distribution within the host tissues and the bacteriome organization have not been systematically studied using modern microscopy techniques. Here, we combine light, fluorescence, and transmission electron microscopy with phylogenomic data to compare symbiont tissue distributions and the bacteriome organization across planthoppers representing 15 families. We identify and describe seven primary types of symbiont localization and seven types of the organization of the bacteriome. We show that Sulcia and Vidania, when present, usually occupy distinct bacteriomes distributed within the body cavity. The more recently acquired gammaproteobacterial and fungal symbionts generally occupy separate groups of cells organized into distinct bacteriomes or mycetomes, distinct from those with Sulcia and Vidania. They can also be localized in the cytoplasm of fat body cells. Alphaproteobacterial symbionts colonize a wider range of host body habitats: Asaia-like symbionts often colonize the host gut lumen, whereas Wolbachia and Rickettsia are usually scattered across insect tissues and cell types, including cells containing other symbionts, bacteriome sheath, fat body cells, gut epithelium, as well as hemolymph. However, there are exceptions, including Gammaproteobacteria that share bacteriome with Vidania, or Alphaproteobacteria that colonize Sulcia cells. We discuss how planthopper symbiont localization correlates with their acquisition and replacement patterns and the symbionts' likely functions. We also discuss the evolutionary consequences, constraints, and significance of these findings.

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.

The genomic basis of evolutionary novelties in a leafhopper

Evolutionary innovations generate phenotypic and species diversity. Elucidating the genomic processes underlying such innovations is central to understanding biodiversity. In this study, we addressed the genomic basis of evolutionary novelties in the glassy-winged sharpshooter (Homalodisca vitripennis, GWSS), an agricultural pest. Prominent evolutionary innovations in leafhoppers include brochosomes, proteinaceous structures that are excreted and used to coat the body, and obligate symbiotic associations with two bacterial types that reside within cytoplasm of distinctive cell types. Using PacBio long-read sequencing and Dovetail Omni-C technology, we generated a chromosome-level genome assembly for the GWSS and then validated the assembly using flow cytometry and karyotyping. Additional transcriptomic and proteomic data were used to identify novel genes that underlie brochosome production. We found that brochosome-associated genes include novel gene families that have diversified through tandem duplications. We also identified the locations of genes involved in interactions with bacterial symbionts. Ancestors of the GWSS acquired bacterial genes through horizontal gene transfer (HGT), and these genes appear to contribute to symbiont support. Using a phylogenomics approach, we inferred HGT sources and timing. We found that some HGT events date to the common ancestor of the hemipteran suborder Auchenorrhyncha, representing some of the oldest known examples of HGT in animals. Overall, we show that evolutionary novelties in leafhoppers are generated by the combination of acquiring novel genes, produced both de novo and through tandem duplication, acquiring new symbiotic associations that enable use of novel diets and niches, and recruiting foreign genes to support symbionts and enhance herbivory.

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.

Substantial rearrangements, single nucleotide frameshift deletion and low diversity in mitogenome of Wolbachia-infected strepsipteran endoparasitoid in comparison to its tephritid hosts

Insect mitogenome organisation is highly conserved, yet, some insects, especially with parasitic life cycles, have rearranged mitogenomes. Furthermore, intraspecific mitochondrial diversity can be reduced by fitness-affecting bacterial endosymbionts like Wolbachia due to their maternal coinheritance with mitochondria. We have sequenced mitogenomes of the Wolbachia-infected endoparasitoid Dipterophagus daci (Strepsiptera: Halictophagidae) and four of its 22 known tephritid fruit fly host species using total genomic extracts of parasitised flies collected across > 700 km in Australia. This halictophagid mitogenome revealed extensive rearrangements relative to the four fly mitogenomes which exhibited the ancestral insect mitogenome pattern. Compared to the only four available other strepsipteran mitogenomes, the D. daci mitogenome had additional transpositions of one rRNA and two tRNA genes, and a single nucleotide frameshift deletion in nad5 requiring translational frameshifting or, alternatively, resulting in a large protein truncation. Dipterophagus daci displays an almost completely endoparasitic life cycle when compared to Strepsiptera that have maintained the ancestral state of free-living adults. Our results support the hypothesis that the transition to extreme endoparasitism evolved together with increased levels of mitogenome changes. Furthermore, intraspecific mitogenome diversity was substantially smaller in D. daci than the parasitised flies suggesting Wolbachia reduced mitochondrial diversity because of a role in D. daci fitness.

Microbial evolution and transitions along the parasite-mutualist continuum

Virtually all plants and animals, including humans, are home to symbiotic microorganisms. Symbiotic interactions can be neutral, harmful or have beneficial effects on the host organism. However, growing evidence suggests that microbial symbionts can evolve rapidly, resulting in drastic transitions along the parasite-mutualist continuum. In this Review, we integrate theoretical and empirical findings to discuss the mechanisms underpinning these evolutionary shifts, as well as the ecological drivers and why some host-microorganism interactions may be stuck at the end of the continuum. In addition to having biomedical consequences, understanding the dynamic life of microorganisms reveals how symbioses can shape an organism's biology and the entire community, particularly in a changing world.

Comparative Genomics Underlines Multiple Roles of Profftella, an Obligate Symbiont of Psyllids: Providing Toxins, Vitamins, and Carotenoids

The Asian citrus psyllid Diaphorina citri (Insecta: Hemiptera: Psylloidea), a serious pest of citrus species worldwide, harbors vertically transmitted intracellular mutualists, Candidatus Profftella armatura (Profftella_DC, Gammaproteobacteria: Burkholderiales) and Candidatus Carsonella ruddii (Carsonella_DC, Gammaproteobacteria: Oceanospirillales). Whereas Carsonella_DC is a typical nutritional symbiont, Profftella_DC is a unique defensive symbiont with organelle-like features, including intracellular localization within the host, perfect infection in host populations, vertical transmission over evolutionary time, and drastic genome reduction down to much less than 1 Mb. Large parts of the 460-kb genome of Profftella_DC are devoted to genes for synthesizing a polyketide toxin; diaphorin. To better understand the evolution of this unusual symbiont, the present study analyzed the genome of Profftella_Dco, a sister lineage to Profftella_DC, using Diaphorina cf. continua, a host psyllid congeneric with D. citri. The genome of coresiding Carsonella (Carsonella_Dco) was also analyzed. The analysis revealed nearly perfect synteny conservation in these genomes with their counterparts from D. citri. The substitution rate analysis further demonstrated genomic stability of Profftella which is comparable to that of Carsonella. Profftella_Dco and Profftella_DC shared all genes for the biosynthesis of diaphorin, hemolysin, riboflavin, biotin, and carotenoids, underlining multiple roles of Profftella, which may contribute to stabilizing symbiotic relationships with the host. However, acyl carrier proteins were extensively amplified in polyketide synthases DipP and DipT for diaphorin synthesis in Profftella_Dco. This level of acyl carrier protein augmentation, unprecedented in modular polyketide synthases of any known organism, is not thought to influence the polyketide structure but may improve the synthesis efficiency.

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.

Hawaiian Nysius Insects Rely on an Obligate Symbiont with a Reduced Genome That Retains a Discrete Nutritional Profile to Match Their Plant Seed Diet

Seed-feeding Nysius insects (Hemiptera: Lygaeidae) have a symbiotic association with distinct intracellular bacteria, “Candidatus Schneideria nysicola” (Gammaproteobacteria). Although many other hemipteran insect groups generally rely on bacterial symbionts that synthesize all ten essential amino acids lacking in their plant sap diets, the nutritional role of Schneideria in Nysius hosts that specialize on a more nutritionally complete seed-based diet has remained unknown. To determine the nutritional and functional capabilities of Schneideria, we sequenced the complete Schneideria genomes from three distantly related endemic Hawaiian Nysius seed bug species. The complete Schneideria genomes are highly conserved and perfectly syntenic among Hawaiian Nysius host species. Each circular chromosome is ∼0.57 Mb in size and encodes 537 protein-coding genes. They further exhibit a strong A + T nucleotide substitution bias with an average G + C nucleotide content of 29%. The predicted nutritional contribution of Schneideria includes four B vitamins and five of the ten essential amino acids that likely match its hosts’ seed-based diet. Disrupted and degraded genes in Schneideria suggests that Hawaiian lineages are undergoing continued gene losses observed in the smaller genomes of the other more ancient hemipteran symbionts.

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.

The Gut Microbiota of the Insect Infraorder Pentatomomorpha (Hemiptera: Heteroptera) for the Light of Ecology and Evolution

The stinkbugs of the infraorder Pentatomomorpha are a group of important plant sap-feeding insects, which host diverse microorganisms. Some are located in their complex morphological midgut compartments, while some within the specialized bacteriomes of insect hosts. This perpetuation of symbioses through host generations is reinforced via the diverse routes of vertical transmission or environmental acquisition of the symbionts. These symbiotic partners, reside either through the extracellular associations in midgut or intracellular associations in specialized cells, not only have contributed nutritional benefits to the insect hosts but also shaped their ecological and evolutionary basis. The stinkbugs and gut microbe symbioses present a valuable model that provides insights into symbiotic interactions between agricultural insects and microorganisms and may become potential agents for insect pest management.

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

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

Further evidence that mechanisms of host/symbiont integration are dissimilar in the maternal versus embryonic Acyrthosiphon pisum bacteriome

Background

Host/symbiont integration is a signature of evolutionarily ancient, obligate endosymbioses. However, little is known about the cellular and developmental mechanisms of host/symbiont integration at the molecular level. Many insects possess obligate bacterial endosymbionts that provide essential nutrients. To advance understanding of the developmental and metabolic integration of hosts and endosymbionts, we track the localization of a non-essential amino acid transporter, ApNEAAT1, across asexual embryogenesis in the aphid, Acyrthosiphon pisum. Previous work in adult bacteriomes revealed that ApNEAAT1 functions to exchange non-essential amino acids at the A. pisum/Buchnera aphidicola symbiotic interface. Driven by amino acid concentration gradients, ApNEAAT1 moves proline, serine, and alanine from A. pisum to Buchnera and cysteine from Buchnera to A. pisum. Here, we test the hypothesis that ApNEAAT1 is localized to the symbiotic interface during asexual embryogenesis.

Results

During A. pisum asexual embryogenesis, ApNEAAT1 does not localize to the symbiotic interface. We observed ApNEAAT1 localization to the maternal follicular epithelium, the germline, and, in late-stage embryos, to anterior neural structures and insect immune cells (hemocytes). We predict that ApNEAAT1 provisions non-essential amino acids to developing oocytes and embryos, as well as to the brain and related neural structures. Additionally, ApNEAAT1 may perform roles related to host immunity.

Conclusions

Our work provides further evidence that the embryonic and adult bacteriomes of asexual A. pisum are not equivalent. Future research is needed to elucidate the developmental time point at which the bacteriome reaches maturity.

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.

Endosymbionts of Metazoans Dwelling in the PACManus Hydrothermal Vent: Diversity and Potential Adaptive Features Revealed by Genome Analysis

Deep-sea hydrothermal vents are dominated by several invertebrate species. The establishment of symbiosis has long been thought to be the key to successful colonization by these sedentary species in such harsh environments. However, the relationships between symbiotic bacteria and their hosts and their role in environmental adaptations generally remain unclear. In this paper, we show that the distribution of three host species showed characteristic niche partitioning in the Manus Basin, giving us the opportunity to understand how they adapt to their particular habitats. This study also revealed three novel genomes of symbionts from the snails of A. boucheti. Combined with a data set on other ectosymbiont and free-living bacteria, genome comparisons for the snail endosymbionts pointed to several genetic traits that may have contributed to the lifestyle shift of Epsilonproteobacteria into the epithelial cells. These findings could increase our understanding of invertebrate-endosymbiont relationships in deep-sea ecosystems.

Deep-sea hydrothermal vent communities are dominated by invertebrates, namely, bathymodiolin mussels, siboglinid tubeworms, and provannid snails. Symbiosis is considered key to successful colonization by these sedentary species in such extreme environments. In the PACManus vent fields, snails, tubeworms, and mussels each colonized a niche with distinct geochemical characteristics. To better understand the metabolic potentials and genomic features contributing to host-environment adaptation, we compared the genomes of the symbionts of Bathymodiolus manusensis, Arcovestia ivanovi, and Alviniconcha boucheti sampled at PACManus, and we discuss their environmentally adaptive features. We found that B. manusensis and A. ivanovi are colonized by Gammaproteobacteria from distinct clades, whereas endosymbionts of B. manusensis feature high intraspecific heterogeneity with differing metabolic potentials. A. boucheti harbored three novel Epsilonproteobacteria symbionts, suggesting potential species-level diversity of snail symbionts. Genome comparisons revealed that the relative abundance of gene families related to low-pH homeostasis, metal resistance, oxidative stress resistance, environmental sensing/responses, and chemotaxis and motility was the highest in A. ivanovi’s symbiont, followed by symbionts of the vent-mouth-dwelling snail A. boucheti, and was relatively low in the symbiont of the vent-periphery-dwelling mussel B. manusensis, which is consistent with their environmental adaptations and host-symbiont interactions. Gene families classified as encoding host interaction/attachment, virulence factors/toxins, and eukaryotic-like proteins were most abundant in symbionts of mussels and least abundant in those of snails, indicating that these symbionts may differ in their host colonization strategies. Comparison of Epsilonproteobacteria symbionts to nonsymbionts demonstrated that the expanded gene families in symbionts were related to vitamin B₁₂ synthesis, toxin-antitoxin systems, methylation, and lipopolysaccharide biosynthesis, suggesting that these are vital to symbiont establishment and development in Epsilonproteobacteria.

IMPORTANCE Deep-sea hydrothermal vents are dominated by several invertebrate species. The establishment of symbiosis has long been thought to be the key to successful colonization by these sedentary species in such harsh environments. However, the relationships between symbiotic bacteria and their hosts and their role in environmental adaptations generally remain unclear. In this paper, we show that the distribution of three host species showed characteristic niche partitioning in the Manus Basin, giving us the opportunity to understand how they adapt to their particular habitats. This study also revealed three novel genomes of symbionts from the snails of A. boucheti. Combined with a data set on other ectosymbiont and free-living bacteria, genome comparisons for the snail endosymbionts pointed to several genetic traits that may have contributed to the lifestyle shift of Epsilonproteobacteria into the epithelial cells. These findings could increase our understanding of invertebrate-endosymbiont relationships in deep-sea ecosystems.

Ongoing Coevolution of Wolbachia and a Widespread Invasive Ant, Anoplolepis gracilipes

While Wolbachia are commonly found among arthropods, intraspecific infection rates can vary substantially across the geographic populations. Here we report nearly 100% prevalence of Wolbachia in the global populations of the yellow crazy ant, Anoplolepis gracilipes. To understand coevolutionary history between Wolbachia and A. gracilipes, we identified single nucleotide polymorphisms (SNPs) in Wolbachia from the ant across 12 geographical regions and compared the phylogeny of SNP-based Wolbachia to patterns of the ant's mitochondrial DNA (mtDNA) variation. Our results revealed a strong concordance between phylogenies of Wolbachia and host mtDNA, providing immediate evidence of co-divergence. Among eight identified SNP loci separating the genetic clusters of Wolbachia, seven loci are located in potential protein-coding genes, three of which being non-synonymous SNPs that may influence gene functions. We found a Wolbachia hypothetical protein gene with signature of positive selection. These findings jointly allow us to characterize Wolbachia-ant coevolution and also raise a question about mechanism(s) underlying maintenance of high prevalence of Wolbachia during the colonization of this invasive ant.

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

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

The Life of an Insect Endosymbiont from the Cradle to the Grave

Host-beneficial endosymbioses, which are formed when a microorganism takes up residence inside another cell and provides a fitness advantage to the host, have had a dramatic influence on the evolution of life. These intimate relationships have yielded the mitochondrion and the plastid (chloroplast) - the ancient organelles that in part define eukaryotic life - along with many more recent associations involving a wide variety of hosts and microbial partners. These relationships are often envisioned as stable associations that appear cooperative and persist for extremely long periods of time. But recent evidence suggests that this stable state is often born from turbulent and conflicting origins, and that the apparent stability of many beneficial endosymbiotic relationships - although certainly real in many cases - is not an inevitable outcome of these associations. Here we review how stable endosymbioses form, how they are maintained, and how they sometimes break down and are reborn. We focus on relationships formed by insects and their resident microorganisms because these symbioses have been the focus of significant empirical work over the last two decades. We review these relationships over five life stages: origin, birth, middle age, old age, and death.

Insect-Symbiont Gene Expression in the Midgut Bacteriocytes of a Blood-Sucking Parasite

Animals interact with a diverse array of both beneficial and detrimental microorganisms. In insects, these symbioses in many cases allow feeding on nutritionally unbalanced diets. It is, however, still not clear how are obligate symbioses maintained at the cellular level for up to several hundred million years. Exact mechanisms driving host-symbiont interactions are only understood for a handful of model species and data on blood-feeding hosts with intracellular bacteria are particularly scarce. Here, we analyzed interactions between an obligately blood-sucking parasite of sheep, the louse fly Melophagus ovinus, and its obligate endosymbiont, Arsenophonus melophagi. We assembled a reference transcriptome for the insect host and used dual RNA-Seq with five biological replicates to compare expression in the midgut cells specialized for housing symbiotic bacteria (bacteriocytes) to the rest of the gut (foregut-hindgut). We found strong evidence for the importance of zinc in the system likely caused by symbionts using zinc-dependent proteases when acquiring amino acids, and for different immunity mechanisms controlling the symbionts than in closely related tsetse flies. Our results show that cellular and nutritional interactions between this blood-sucking insect and its symbionts are less intimate than what was previously found in most plant-sap sucking insects. This finding is likely interconnected to several features observed in symbionts in blood-sucking arthropods, particularly their midgut intracellular localization, intracytoplasmic presence, less severe genome reduction, and relatively recent associations caused by frequent evolutionary losses and replacements.

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

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

Evolutionary rates are correlated between cockroach symbionts and mitochondrial genomes

Bacterial endosymbionts evolve under strong host-driven selection. Factors influencing host evolution might affect symbionts in similar ways, potentially leading to correlations between the molecular evolutionary rates of hosts and symbionts. Although there is evidence of rate correlations between mitochondrial and nuclear genes, similar investigations of hosts and symbionts are lacking. Here, we demonstrate a correlation in molecular rates between the genomes of an endosymbiont (Blattabacterium cuenoti) and the mitochondrial genomes of their hosts (cockroaches). We used partial genome data for multiple strains of B. cuenoti to compare phylogenetic relationships and evolutionary rates for 55 cockroach/symbiont pairs. The phylogenies inferred for B. cuenoti and the mitochondrial genomes of their hosts were largely congruent, as expected from their identical maternal and cytoplasmic mode of inheritance. We found a correlation between evolutionary rates of the two genomes, based on comparisons of root-to-tip distances and on comparisons of the branch lengths of phylogenetically independent species pairs. Our results underscore the profound effects that long-term symbiosis can have on the biology of each symbiotic partner.

Trends in Symbiont-Induced Host Cellular Differentiation

Bacteria participate in a wide diversity of symbiotic associations with eukaryotic hosts that require precise interactions for bacterial recognition and persistence. Most commonly, host-associated bacteria interfere with host gene expression to modulate the immune response to the infection. However, many of these bacteria also interfere with host cellular differentiation pathways to create a hospitable niche, resulting in the formation of novel cell types, tissues, and organs. In both of these situations, bacterial symbionts must interact with eukaryotic regulatory pathways. Here, we detail what is known about how bacterial symbionts, from pathogens to mutualists, control host cellular differentiation across the central dogma, from epigenetic chromatin modifications, to transcription and mRNA processing, to translation and protein modifications. We identify four main trends from this survey. First, mechanisms for controlling host gene expression appear to evolve from symbionts co-opting cross-talk between host signaling pathways. Second, symbiont regulatory capacity is constrained by the processes that drive reductive genome evolution in host-associated bacteria. Third, the regulatory mechanisms symbionts exhibit correlate with the cost/benefit nature of the association. And, fourth, symbiont mechanisms for interacting with host genetic regulatory elements are not bound by native bacterial capabilities. Using this knowledge, we explore how the ubiquitous intracellular Wolbachia symbiont of arthropods and nematodes may modulate host cellular differentiation to manipulate host reproduction. Our survey of the literature on how infection alters gene expression in Wolbachia and its hosts revealed that, despite their intermediate-sized genomes, different strains appear capable of a wide diversity of regulatory manipulations. Given this and Wolbachia's diversity of phenotypes and eukaryotic-like proteins, we expect that many symbiont-induced host differentiation mechanisms will be discovered in this system.

Changes in Aphid Host Plant Diet Influence the Small-RNA Expression Profiles of Its Obligate Nutritional Symbiont, Buchnera

Plants are a difficult food resource to use, and herbivorous insects have evolved a variety of mechanisms that allow them to fully exploit this poor nutritional resource. One such mechanism is the maintenance of bacterial symbionts that aid in host plant feeding and development. The majority of these intracellular symbionts have highly eroded genomes that lack many key regulatory genes; consequently, it is unclear if these symbionts can respond to changes in the insect's diet to facilitate host plant use. There is emerging evidence that symbionts with highly eroded genomes express small RNAs (sRNAs), some of which potentially regulate gene expression. In this study, we sought to determine if the reduced genome of the nutritional symbiont (Buchnera) in the pea aphid responds to changes in the aphid's host plant diet. Using transcriptome sequencing (RNA-seq), Buchnera sRNA expression profiles were characterized within two Buchnera life stages when pea aphids fed on either alfalfa or fava bean. Overall, this study demonstrates that Buchnera sRNA expression changes not only with life stage but also with changes in aphid host plant diet. Of the 321 sRNAs characterized in this study, 47% were previously identified and 22% showed evidence of conservation in two or more Buchnera taxa. Functionally, 13 differentially expressed sRNAs were predicted to target genes related to pathways involved in essential amino acid biosynthesis. Overall, results from this study reveal that host plant diet influences the expression of conserved and lineage-specific sRNAs in Buchnera and that these sRNAs display distinct host plant-specific expression profiles among biological replicates.IMPORTANCE In general, the genomes of intracellular bacterial symbionts are reduced compared to those of free-living relatives and lack many key regulatory genes. Many of these reduced genomes belong to obligate mutualists of insects that feed on a diet that is deficient in essential nutrients, such as essential amino acids. It is unclear if these symbionts respond with their host to changes in insect diet, because of their reduced regulatory capacity. Emerging evidence suggests that these symbionts express small RNAs (sRNAs) that regulate gene expression at the posttranscriptional level. Therefore, in this study, we sought to determine if the reduced genome of the nutritional symbiont Buchnera in the pea aphid responds to changes in the aphid's host plant diet. This study demonstrates for the first time that Buchnera sRNAs, some conserved in two or more Buchnera lineages, are differentially expressed when aphids feed on different plant species and potentially target genes within essential amino acid biosynthesis pathways.

Cicada Endosymbionts Have tRNAs That Are Correctly Processed Despite Having Genomes That Do Not Encode All of the tRNA Processing Machinery

The smallest bacterial genomes, in the range of about 0.1 to 0.5 million base pairs, are commonly found in the nutritional endosymbionts of insects. These tiny genomes are missing genes that encode proteins and RNAs required for the translation of mRNAs, one of the most highly conserved and important cellular processes. In this study, we found that the bacterial endosymbionts of cicadas have genomes which encode incomplete tRNA sets and lack genes required for tRNA processing. Nevertheless, we found that endosymbiont tRNAs are correctly processed at their 5′ and 3′ ends and, surprisingly, that mostly exist as tRNA halves. We hypothesize that the cicada host must supply its symbionts with these missing tRNA processing activities.

Gene loss and genome reduction are defining characteristics of endosymbiotic bacteria. The most highly reduced endosymbiont genomes have lost numerous essential genes related to core cellular processes such as replication, transcription, and translation. Computational gene predictions performed for the genomes of the two bacterial symbionts of the cicada Diceroprocta semicincta, “Candidatus Hodgkinia cicadicola” (Alphaproteobacteria) and “Ca. Sulcia muelleri” (Bacteroidetes), have found only 26 and 16 tRNA genes and 15 and 10 aminoacyl tRNA synthetase genes, respectively. Furthermore, the original “Ca. Hodgkinia cicadicola” genome annotation was missing several essential genes involved in tRNA processing, such as those encoding RNase P and CCA tRNA nucleotidyltransferase as well as several RNA editing enzymes required for tRNA maturation. How these cicada endosymbionts perform basic translation-related processes remains unknown. Here, by sequencing eukaryotic mRNAs and total small RNAs, we show that the limited tRNA set predicted by computational annotation of “Ca. Sulcia muelleri” and “Ca. Hodgkinia cicadicola” is likely correct. Furthermore, we show that despite the absence of genes encoding tRNA processing activities in the symbiont genomes, symbiont tRNAs have correctly processed 5′ and 3′ ends and seem to undergo nucleotide modification. Surprisingly, we found that most “Ca. Hodgkinia cicadicola” and “Ca. Sulcia muelleri” tRNAs exist as tRNA halves. We hypothesize that “Ca. Sulcia muelleri” and “Ca. Hodgkinia cicadicola” tRNAs function in bacterial translation but require host-encoded enzymes to do so.