Regulation of transcription in a reduced bacterial genome: nutrient-provisioning genes of the obligate symbiont Buchnera aphidicola

Dynamic response of essential amino acid biosynthesis in Buchnera aphidicola to supplement sub-optimal host nutrition

The endosymbiotic bacterium Buchnera aphidicola allows its host Acyrthosiphon pisum to utilise a nutritionally limited phloem sap diet without significant mortality by providing essential amino acids (EAAs), which it biosynthesises de novo via complex pathways consisting of multiple enzymes. Previous studies have reported how non-essential amino acids (NEAAs) provided by the host are utilised by B. aphidicola, along with how genes within the biosynthetic pathways respond to amino acid deficiency. Although the effect on B. aphidicola gene expression upon the removal of a single EAA and multiple NEAAs from the A. pisum diet has been reported, little is known about the effects of the complete simultaneous removal of multiple EAAs, especially branched-chain amino acids (BCAAs). To investigate this, A. pisum was provided with amino acid deficient diets ilv- (lacking isoleucine, leucine, valine) or thra- (lacking threonine, methionine, lysine). Due to their involvement in the production of several amino acids, the expression of genes ilvC, ilvD (both involved in isoleucine, leucine and valine biosynthesis) and thrA (involved in threonine, methionine and lysine biosynthesis) was analysed and the expression of trpC (involved in tryptophan biosynthesis) was used as a control. Survival was reduced significantly when A. pisum was reared on ilv- or thra- (P < 0.001 and P = 0.000 respectively) compared to optimal artificial diet and was significantly lower on ilv- (P < 0.001) than thra-. This is likely attributed to the EAAs absent from ilv- being required at higher concentrations for aphid growth, than those EAAs absent from thra-. Expression of ilvC and ilvD were upregulated 2.49- and 2.08-fold (respectively) and thrA expression increased 2.35- and 2.12-fold when A. pisum was reared on ilv- and thra- (respectively). The surprisingly large upregulation of thrA when reared on ilv- is likely due to threonine being an intermediate in isoleucine biosynthesis. Expression of trpC was not affected by rearing on either of the two amino acid deficient diets. To our knowledge this study has shown, for the first time, how genes within the biosynthetic pathways of an endosymbiont respond to the simultaneous complete omission of multiple EAAs as well as all three BCAAs (leucine, isoleucine, valine), from the host diet.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Coordination of host and endosymbiont gene expression governs endosymbiont growth and elimination in the cereal weevil Sitophilus spp

BACKGROUND

Insects living in nutritionally poor environments often establish long-term relationships with intracellular bacteria that supplement their diets and improve their adaptive and invasive powers. Even though these symbiotic associations have been extensively studied on physiological, ecological, and evolutionary levels, few studies have focused on the molecular dialogue between host and endosymbionts to identify genes and pathways involved in endosymbiosis control and dynamics throughout host development.

RESULTS

We simultaneously analyzed host and endosymbiont gene expression during the life cycle of the cereal weevil Sitophilus oryzae, from larval stages to adults, with a particular emphasis on emerging adults where the endosymbiont Sodalis pierantonius experiences a contrasted growth-climax-elimination dynamics. We unraveled a constant arms race in which different biological functions are intertwined and coregulated across both partners. These include immunity, metabolism, metal control, apoptosis, and bacterial stress response.

CONCLUSIONS

The study of these tightly regulated functions, which are at the center of symbiotic regulations, provides evidence on how hosts and bacteria finely tune their gene expression and respond to different physiological challenges constrained by insect development in a nutritionally limited ecological niche. Video Abstract.

Arsenophonus Interacts with Buchnera to Improve Growth Performance of Aphids under Amino Acid Stress

Amino acids play a crucial role in the growth and development of insects. Aphids cannot ingest enough amino acids in plant phloem to meet their requirements, and therefore, they are mainly dependent on the obligate symbiont Buchnera aphidicola to synthesize essential amino acids. Besides Buchnera, aphids may harbor another facultative symbiont, Arsenophonus, which alters the requirement of the cotton-melon aphid Aphis gossypii for amino acid. However, it is unclear how Arsenophonus regulates the requirement. Here, we found that Arsenophonus ameliorated growth performance of A. gossypii on an amino acid-deficient diet. A deficiency in lysine (Lys) or methionine (Met) led to changes in the abundance of Arsenophonus. Arsenophonus suppressed the abundance of Buchnera when aphids were fed a normal amino acid diet, but this suppression was eliminated or reversed when aphids were on a Lys- or Met-deficient diet. The relative abundance of Arsenophonus was positively correlated with that of Buchnera, but neither of them was correlated with the body weight of aphids. The relative expression levels of Lys and Met synthase genes of Buchnera were affected by the interaction between Arsenophonus infections and Buchnera abundance, especially in aphids reared on a Lys- or Met-deficient diet. Arsenophonus coexisted with Buchnera in bacteriocytes, which strengthens the interaction. IMPORTANCE The obligate symbiont Buchnera can synthesize amino acids for aphids. In this study, we found that a facultative symbiont, Arsenophonus, can help improve aphids' growth performance under amino acid deficiency stress by changing the relative abundance of Buchnera and the expression levels of amino acid synthase genes. This study highlights the interaction between Arsenophonus and Buchnera to ameliorate aphid growth under amino acid stress.

No Transcriptional Compensation for Extreme Gene Dosage Imbalance in Fragmented Bacterial Endosymbionts of Cicadas

Bacteria that form long-term intracellular associations with host cells lose many genes, a process that often results in tiny, gene-dense, and stable genomes. Paradoxically, the some of the same evolutionary processes that drive genome reduction and simplification may also cause genome expansion and complexification. A bacterial endosymbiont of cicadas, Hodgkinia cicadicola, exemplifies this paradox. In many cicada species, a single Hodgkinia lineage with a tiny, gene-dense genome has split into several interdependent cell and genome lineages. Each new Hodgkinia lineage encodes a unique subset of the ancestral unsplit genome in a complementary way, such that the collective gene contents of all lineages match the total found in the ancestral single genome. This splitting creates genetically distinct Hodgkinia cells that must function together to carry out basic cellular processes. It also creates a gene dosage problem where some genes are encoded by only a small fraction of cells while others are much more abundant. Here, by sequencing DNA and RNA of Hodgkinia from different cicada species with different amounts of splitting-along with its structurally stable, unsplit partner endosymbiont Sulcia muelleri-we show that Hodgkinia does not transcriptionally compensate to rescue the wildly unbalanced gene and genome ratios that result from lineage splitting. We also find that Hodgkinia has a reduced capacity for basic transcriptional control independent of the splitting process. Our findings reveal another layer of degeneration further pushing the limits of canonical molecular and cell biology in Hodgkinia and may partially explain its propensity to go extinct through symbiont replacement.

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

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

Single mutation makes Escherichia coli an insect mutualist

Microorganisms often live in symbiosis with their hosts, and some are considered mutualists, where all species involved benefit from the interaction. How free-living microorganisms have evolved to become mutualists is unclear. Here we report an experimental system in which non-symbiotic Escherichia coli evolves into an insect mutualist. The stinkbug Plautia stali is typically associated with its essential gut symbiont, Pantoea sp., which colonizes a specialized symbiotic organ. When sterilized newborn nymphs were infected with E. coli rather than Pantoea sp., only a few insects survived, in which E. coli exhibited specific localization to the symbiotic organ and vertical transmission to the offspring. Through transgenerational maintenance with P. stali, several hypermutating E. coli lines independently evolved to support the host's high adult emergence and improved body colour; these were called 'mutualistic' E. coli. These mutants exhibited slower bacterial growth, smaller size, loss of flagellar motility and lack of an extracellular matrix. Transcriptomic and genomic analyses of 'mutualistic' E. coli lines revealed independent mutations that disrupted the carbon catabolite repression global transcriptional regulator system. Each mutation reproduced the mutualistic phenotypes when introduced into wild-type E. coli, confirming that single carbon catabolite repression mutations can make E. coli an insect mutualist. These findings provide an experimental system for future work on host-microbe symbioses and may explain why microbial mutualisms are omnipresent in nature.

Roles of Bacterial Symbionts in Transmission of Plant Virus by Hemipteran Vectors

The majority of plant viruses are transmitted by hemipteran insects. Bacterial symbionts in hemipteran hosts have a significant impact on the host life, physiology and ecology. Recently, the involvement of bacterial symbionts in hemipteran vector-virus and vector-plant interactions has been documented. Thus, the exploitation and manipulation of bacterial symbionts have great potential for plant viral disease control. Herein, we review the studies performed on the impact of symbiotic bacteria on plant virus transmission, including insect-bacterial symbiont associations, the role of these bacterial symbionts in viral acquisition, stability and release during viral circulation in insect bodies, and in viral vertical transmission. Besides, we prospect further studies aimed to understand tripartite interactions of the virus-symbiotic microorganisms-insect vector.

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.

Insect-host control of obligate, intracellular symbiont density

Many insects rely on intracellular bacterial symbionts to supplement their specialized diets with micronutrients. Using data from diverse and well-studied insect systems, we propose three lines of evidence suggesting that hosts have tight control over the density of their obligate, intracellular bacterial partners. First, empirical studies have demonstrated that the within-host symbiont density varies depending on the nutritional and developmental requirements of the host. Second, symbiont genomes are highly reduced and have limited capacity for self-replication or transcriptional regulation. Third, several mechanisms exist for hosts to tolerate, regulate and remove symbionts including physical compartmentalization and autophagy. We then consider whether such regulation is adaptive, by discussing the relationship between symbiont density and host fitness. We discuss current limitations of empirical studies for exploring fitness effects in host-symbiont relationships, and emphasize the potential for using mathematical models to formalize evolutionary hypotheses and to generate testable predictions for future work.

RNA polymerases in strict endosymbiont bacteria with extreme genome reduction show distinct erosions that might result in limited and differential promoter recognition

Strict endosymbiont bacteria present high degree genome reduction, retain smaller proteins, and in some instances, lack complete functional domains compared to free-living counterparts. Until now, the mechanisms underlying these genetic reductions are not well understood. In this study, the conservation of RNA polymerases, the essential machinery for gene expression, is analyzed in endosymbiont bacteria with extreme genome reductions. We analyzed the RNA polymerase subunits to identify and define domains, subdomains, and specific amino acids involved in precise biological functions known in Escherichia coli. We also perform phylogenetic analysis and three-dimensional models over four lineages of endosymbiotic proteobacteria with the smallest genomes known to date: Candidatus Hodgkinia cicadicola, Candidatus Tremblaya phenacola, Candidatus Tremblaya Princeps, Candidatus Nasuia deltocephalinicola, and Candidatus Carsonella ruddii. We found that some Hodgkinia strains do not encode for the RNA polymerase α subunit. The rest encode genes for α, β, β', and σ subunits to form the RNA polymerase. However, 16% shorter, on average, respect their orthologous in E. coli. In the α subunit, the amino-terminal domain is the most conserved. Regarding the β and β' subunits, both the catalytic core and the assembly domains are the most conserved. However, they showed compensatory amino acid substitutions to adapt to changes in the σ subunit. Precisely, the most erosive diversity occurs within the σ subunit. We identified broad amino acid substitution even in those recognizing and binding to the -10-box promoter element. In an overall conceptual image, the RNA polymerase from Candidatus Nasuia conserved the highest similarity with Escherichia coli RNA polymerase and their σ70. It might be recognizing the two main promoter elements (-10 and -35) and the two promoter accessory elements (-10 extended and UP-element). In Candidatus Carsonella, the RNA polymerase could recognize all the promoter elements except the -10-box extended. In Candidatus Tremblaya and Hodgkinia, due to the α carboxyl-terminal domain absence, they might not recognize the UP-promoter element. We also identified the lack of the β flap-tip helix domain in most Hodgkinia's that suggests the inability to bind the -35-box promoter element.

Host-Endosymbiont Genome Integration in a Deep-Sea Chemosymbiotic Clam

Endosymbiosis with chemosynthetic bacteria has enabled many deep-sea invertebrates to thrive at hydrothermal vents and cold seeps, but most previous studies on this mutualism have focused on the bacteria only. Vesicomyid clams dominate global deep-sea chemosynthesis-based ecosystems. They differ from most deep-sea symbiotic animals in passing their symbionts from parent to offspring, enabling intricate coevolution between the host and the symbiont. Here, we sequenced the genomes of the clam Archivesica marissinica (Bivalvia: Vesicomyidae) and its bacterial symbiont to understand the genomic/metabolic integration behind this symbiosis. At 1.52 Gb, the clam genome encodes 28 genes horizontally transferred from bacteria, a large number of pseudogenes and transposable elements whose massive expansion corresponded to the timing of the rise and subsequent divergence of symbiont-bearing vesicomyids. The genome exhibits gene family expansion in cellular processes that likely facilitate chemoautotrophy, including gas delivery to support energy and carbon production, metabolite exchange with the symbiont, and regulation of the bacteriocyte population. Contraction in cellulase genes is likely adaptive to the shift from phytoplankton-derived to bacteria-based food. It also shows contraction in bacterial recognition gene families, indicative of suppressed immune response to the endosymbiont. The gammaproteobacterium endosymbiont has a reduced genome of 1.03 Mb but retains complete pathways for sulfur oxidation, carbon fixation, and biosynthesis of 20 common amino acids, indicating the host’s high dependence on the symbiont for nutrition. Overall, the host–symbiont genomes show not only tight metabolic complementarity but also distinct signatures of coevolution allowing the vesicomyids to thrive in chemosynthesis-based ecosystems.

Triple RNA-Seq characterizes aphid gene expression in response to infection with unequally virulent strains of the endosymbiont Hamiltonella defensa

Background

Secondary endosymbionts of aphids provide benefits to their hosts, but also impose costs such as reduced lifespan and reproductive output. The aphid Aphis fabae is host to different strains of the secondary endosymbiont Hamiltonella defensa, which encode different putative toxins. These strains have very different phenotypes: They reach different densities in the host, and the costs and benefits (protection against parasitoid wasps) they confer to the host vary strongly.

Results

We used RNA-Seq to generate hypotheses on why four of these strains inflict such different costs to A. fabae. We found different H. defensa strains to cause strain-specific changes in aphid gene expression, but little effect of H. defensa on gene expression of the primary endosymbiont, Buchnera aphidicola. The highly costly and over-replicating H. defensa strain H85 was associated with strongly reduced aphid expression of hemocytin, a marker of hemocytes in Drosophila. The closely related strain H15 was associated with downregulation of ubiquitin-related modifier 1, which is related to nutrient-sensing and oxidative stress in other organisms. Strain H402 was associated with strong differential regulation of a set of hypothetical proteins, the majority of which were only differentially regulated in presence of H402.

Conclusions

Overall, our results suggest that costs of different strains of H. defensa are likely caused by different mechanisms, and that these costs are imposed by interacting with the host rather than the host’s obligatory endosymbiont B. aphidicola.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07742-8.

Enhanced Symbiotic Characteristics in Bacterial Genomes with the Disruption of rRNA Operon

Simple Summary

Exploring the genomic changes that organisms have undergone to adapt to their specific environment is one of the most important processes in ecology and evolutionary biology. Here, we found that almost all rRNA operon-unlinked bacteria are symbiotic bacteria, which could be evidence of specific selective pressures in symbionts like genome reduction. This is meaningful and suggests that not only does the copy number variation of the rRNA operon sensitively respond to the bacterial lifestyle, but structural modification can also strongly reflect adaptation to the surrounding environmental conditions.

Abstract

Ribosomal RNA is an indispensable molecule in living organisms that plays an essential role in protein synthesis. Especially in bacteria, 16S, 23S, and 5S rRNAs are usually co-transcribed as operons. Despite the positive effects of rRNA co-transcription on growth and reproduction rate, a recent study revealed that bacteria with unlinked rRNA operons are more widespread than expected. However, it is still unclear why the rRNA operon is broken. Here, we explored rRNA operon linkage status in 15,898 bacterial genomes and investigated whether they have common features or lifestyles; 574 genomes were found to have unlinked rRNA operons and tended to be phylogenetically conserved. Most of them were symbionts and showed enhanced symbiotic genomic features such as reduced genome size and high adenine–thymine (AT) content. In an eggNOG-mapper analysis, they were also found to have significantly fewer genes than rRNA operon-linked bacteria in the “transcription” and “energy production and conversion in metabolism” categories. These genomes also tend to decrease RNases related to the synthesis of ribosomes and tRNA processing. Based on these results, the disruption of the rRNA operon seems to be one of the tendencies associated with the characteristics of bacteria requiring a low dynamic range.

Spatial and morphological reorganization of endosymbiosis during metamorphosis accommodates adult metabolic requirements in a weevil

Bacterial intracellular symbiosis (endosymbiosis) is widespread in nature and impacts many biological processes. In holometabolous symbiotic insects, metamorphosis entails a complete and abrupt internal reorganization that creates a constraint for endosymbiont transmission from larvae to adults. To assess how endosymbiosis copes-and potentially evolves-throughout this major host-tissue reorganization, we used the association between the cereal weevil Sitophilus oryzae and the bacterium Sodalis pierantonius as a model system. S. pierantonius are contained inside specialized host cells, the bacteriocytes, that group into an organ, the bacteriome. Cereal weevils require metabolic inputs from their endosymbiont, particularly during adult cuticle synthesis, when endosymbiont load increases dramatically. By combining dual RNA-sequencing analyses and cell imaging, we show that the larval bacteriome dissociates at the onset of metamorphosis and releases bacteriocytes that undergo endosymbiosis-dependent transcriptomic changes affecting cell motility, cell adhesion, and cytoskeleton organization. Remarkably, bacteriocytes turn into spindle cells and migrate along the midgut epithelium, thereby conveying endosymbionts to midgut sites where future mesenteric caeca will develop. Concomitantly, endosymbiont genes encoding a type III secretion system and a flagellum apparatus are transiently up-regulated while endosymbionts infect putative stem cells and enter their nuclei. Infected cells then turn into new differentiated bacteriocytes and form multiple new bacteriomes in adults. These findings show that endosymbiosis reorganization in a holometabolous insect relies on a synchronized host-symbiont molecular and cellular "choreography" and illustrates an adaptive feature that promotes bacteriome multiplication to match increased metabolic requirements in emerging adults.

Performance and Transcriptional Response of the Green Peach Aphid Myzus persicae to the Restriction of Dietary Amino Acids

Free amino acids in the phloem sap are the dominant nitrogen source for aphids, but their availability is usually poor. Although some studies have explored the effect of dietary amino acid restriction on aphid performance, little is known about the molecular basis of these effects. Here, we examined the performance and transcriptome of the green peach aphid, Myzus persicae, fed a standard diet (Control diet) or a diet containing 50% of the total amino acids of the Control diet (Half diet). Aphid weight and fecundity were significantly reduced in the Half diet group. Transcriptomic analysis showed that a total of 1460 genes were differentially expressed between the groups were fed on the two diets, which many of them were associated with nutrient and energy metabolism. When feeding on the Half diet, aphids upregulated genes associated with the amino acid biosynthetic pathway (predominantly amino acid biosynthesis genes and some amino acid transporter genes) as well as the cysteine and serine protease genes. Furthermore, these aphids displayed increased expression of genes associated with glycolysis, which could generate intermediates for de novo amino acid biosynthesis. Consistent with this, elevated glucose levels were observed in aphids in the Half diet group. Additionally, the expression levels of several genes associated with hormonal signaling pathway were altered. Several genes related to juvenile hormone and insulin-like peptide (ILP) signaling were downregulated, including Krüppel homolog 1 (Kr-h1) and insulin-like peptide 5 (Ilp5), respectively. In contrast, several genes related to ecdysone signaling were upregulated including broad-complex core protein (Br-c) and shade (Shd). Despite their poor performances, M. persicae adapted to dietary restriction of amino acids, through upregulation of genes involved in amino acid biosynthesis, glycolysis, and protein degradation, as well as by altering the expression level of genes involved in hormone signaling pathways.

Coordination of host and symbiont gene expression reveals a metabolic tug-of-war between aphids and Buchnera

Symbioses between animals and microbes are often described as mutualistic, but are subject to tradeoffs that may manifest as shifts in host and symbiont metabolism, cellular processes, or symbiont density. In pea aphids, the bacterial symbiont Buchnera is confined to specialized aphid cells called bacteriocytes, where it produces essential amino acids needed by hosts. This relationship is dynamic; Buchnera titer varies within individual aphids and among different clonal aphid lineages, and is affected by environmental and host genetic factors. We examined how host genotypic variation relates to host and symbiont function among seven aphid clones differing in Buchnera titer. We found that bacteriocyte gene expression varies among individual aphids and among aphid clones, and that Buchnera gene expression changes in response. By comparing hosts with low and high Buchnera titer, we found that aphids and Buchnera oppositely regulate genes underlying amino acid biosynthesis and cell growth. In high-titer hosts, both bacteriocytes and symbionts show elevated expression of genes underlying energy metabolism. Several eukaryotic cell signaling pathways are differentially expressed in bacteriocytes of low- versus high-titer hosts: Cell-growth pathways are up-regulated in low-titer genotypes, while membrane trafficking, lysosomal processes, and mechanistic target of rapamycin (mTOR) and cytokine pathways are up-regulated in high-titer genotypes. Specific Buchnera functions are up-regulated within different bacteriocyte environments, with genes underlying flagellar body secretion and flagellar assembly overexpressed in low- and high-titer hosts, respectively. Overall, our results reveal allowances and demands made by both host and symbiont engaged in a metabolic "tug-of-war."

Genome analysis of sponge symbiont 'Candidatus Halichondribacter symbioticus' shows genomic adaptation to a host-dependent lifestyle

The marine sponge Halichondria panicea inhabits coastal areas around the globe and is a widely studied sponge species in terms of its biology, yet the ecological functions of its dominant bacterial symbiont 'Candidatus Halichondribacter symbioticus' remain unknown. Here, we present the draft genome of 'Ca. H. symbioticus' HS1 (2.8 Mbp, ca. 87.6% genome coverage) recovered from the sponge metagenome of H. panicea in order to study functions and symbiotic interactions at the genome level. Functional genome comparison of HS1 against closely related free-living seawater bacteria revealed a reduction of genes associated with carbohydrate transport and transcription regulation, pointing towards a limited carbohydrate metabolism, and static transcriptional dynamics reminiscent of other bacterial symbionts. In addition, HS1 was enriched in sponge symbiont specific gene families related to host-symbiont interactions and defence. Similarity in the functional gene repertoire between HS1 and a phylogenetically more distant symbiont in the marine sponge Aplysina aerophoba, based on COG category distribution, suggest a convergent evolution of symbiont specific traits and general metabolic features. This warrants further investigation into convergent genomic evolution of symbionts across different sponge species and habitats.

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.

microRNA-92a regulates the expression of aphid bacteriocyte-specific secreted protein 1

OBJECTIVE

Aphids harbor a nutritional obligate endosymbiont in specialized cells called bacteriocytes, which aggregate to form an organ known as the bacteriome. Aphid bacteriomes display distinct gene expression profiles that facilitate the symbiotic relationship. Currently, the mechanisms that regulate these patterns of gene expression are unknown. Recently using computational pipelines, we identified miRNAs that are conserved in expression in the bacteriomes of two aphid species and proposed that they function as important regulators of bacteriocyte gene expression. Here using a dual luciferase assay in mouse NIH/3T3 cell culture, we aimed to experimentally validate the computationally predicted interaction between Myzus persicae miR-92a and the predicted target region of M. persicae bacteriocyte-specific secreted protein 1 (SP1) mRNA.

RESULTS

In the dual luciferase assay, miR-92a interacted with the SP1 target region resulting in a significant downregulation of the luciferase signal. Our results demonstrate that miR-92a interacts with SP1 to alter expression in a heterologous expression system, thereby supporting our earlier assertion that miRNAs are regulators of the aphid/Buchnera symbiotic interaction.

Conservation of transcriptional elements in the obligate symbiont of the whitefly Bemisia tabaci

Background

Bacterial symbiosis is widespread in arthropods, especially in insects. Some of the symbionts undergo a long-term co-evolution with the host, resulting in massive genome decay. One particular consequence of genome decay is thought to be the elimination of transcriptional elements within both the coding region and intergenic sequences. In the whitefly Bemisia tabaci species complex, the obligate symbiont Candidatus Portiera aleyrodidarum is of vital importance in nutrient provision, and yet little is known about the regulatory capacities of it.

Methods

Portiera genomes of two whitefly species in China were sequenced and assembled. Gene content of these two Portiera genomes was predicted, and then subjected to Kyoto Encyclopedia of Genes and Genomes pathway analysis. Together with two other Portiera genomes from whitefly species available previously, four Portiera genomes were utilized to investigate regulatory capacities of Portiera, focusing on transcriptional elements, including genes related with transcription and functional elements within the intergenic spacers.

Results

Comparative analyses of the four Portiera genomes of whitefly B. tabaci indicate that the obligate symbionts Portiera is similar in different species of whiteflies, in terms of general genome features and possible functions in the biosynthesis of essential amino acids. The screening of transcriptional factors suggests compromised ability of Portiera to regulate the essential amino acid biosynthesis pathways. Meanwhile, thermal tolerance ability of Portiera is indicated with the detection of a σ³² factor, as well as two predicted σ³² binding sites. Within intergenic spacers, functional elements are predicted, including 37 Shine-Dalgarno sequences and 34 putative small RNAs.

Exaggeration and cooption of innate immunity for social defense

Social insects often exhibit striking altruistic behaviors, of which the most spectacular ones may be self-destructive defensive behaviors called autothysis, "self-explosion," or "suicidal bombing." In the social aphid Nipponaphis monzeni, when enemies damage their plant-made nest called the gall, soldier nymphs erupt to discharge a large amount of body fluid, mix the secretion with their legs, and skillfully plaster it over the plant injury. Dozens of soldiers come out, erupt, mix, and plaster, and the gall breach is promptly sealed with the coagulated body fluid. What molecular and cellular mechanisms underlie the self-sacrificing nest repair with body fluid for the insect society? Here we demonstrate that the body cavity of soldier nymphs is full of highly differentiated large hemocytes that contain huge amounts of lipid droplets and phenoloxidase (PO), whereas their hemolymph accumulates huge amounts of tyrosine and a unique repeat-containing protein (RCP). Upon breakage of the gall, soldiers gather around the breach and massively discharge the body fluid. The large hemocytes rupture and release lipid droplets, which promptly form a lipidic clot, and, concurrently, activated PO converts tyrosine to reactive quinones, which cross-link RCP and other macromolecules to physically reinforce the clot to seal the gall breach. Here, soldiers' humoral and cellular immune mechanisms for wound sealing are extremely up-regulated and utilized for colony defense, which provides a striking case of direct evolutionary connection between individual immunity and social immunity and highlights the importance of exaggeration and cooption of preexisting traits to create evolutionary novelties.

Bacterial communities of Aphis gossypii and Myzus persicae (Hemiptera: Aphididae) from pepper crops (Capsicum sp.)

Insects harbor a wide variety of microorganisms that form complex and changing communities and play an important role in the biology and evolution of their hosts. Aphids have been used as model organisms to study microorganism-insect interactions. Almost all aphids are infected with the obligate endosymbiont Buchnera aphidicola and can host different bacteria that allow them to acquire traits of agronomic importance, such as resistance to high temperatures and/or defense against natural enemies. However, the bacterial communities of most aphid species remain poorly characterized. In this study, we used high-throughput DNA sequencing to characterize the bacterial communities of Aphis gossypii and Myzus persicae from two cultivable pepper species, Capsicum frutescens (Tabasco variety) and C. annuum (Cayenne variety), in four localities of southwestern Colombia. In addition, we evaluated the dynamics of A. gossypii-associated microorganisms on a seasonal basis. Our results show that the bacterial communities of A. gossypii and M. persicae are dominated by the primary endosymbiont B. aphidicola, while the presence of the facultative symbiont Arsenophonus sp. was only detected in one A. gossypii population from cayenne pepper. In addition to these two known symbionts, eight bacterial OTUs were identified that presented a frequency of 1% or more in at least one of the analyzed populations. The results show that the bacterial communities of aphids associated with pepper crops appears to be structured according to the host aphid species and the geographical location, while no differences were observed in the diversity of bacteria between host plants. Finally, the diversity and abundance of the A. gossypii bacterial community was variable among the four sampling points evaluated over the year and showed a relation with the aphid’s population dynamics. This study represents the first approach to the knowledge of the bacterial community present in chili pepper aphids from Colombia. Nevertheless, more in-depth studies, including replicates, are required to confirm the patterns observed in the microbial communities of aphids from pepper crops.

Nephromyces Encodes a Urate Metabolism Pathway and Predicted Peroxisomes, Demonstrating That These Are Not Ancient Losses of Apicomplexans

The phylum Apicomplexa is a quintessentially parasitic lineage, whose members infect a broad range of animals. One exception to this may be the apicomplexan genus Nephromyces, which has been described as having a mutualistic relationship with its host. Here we analyze transcriptome data from Nephromyces and its parasitic sister taxon, Cardiosporidium, revealing an ancestral purine degradation pathway thought to have been lost early in apicomplexan evolution. The predicted localization of many of the purine degradation enzymes to peroxisomes, and the in silico identification of a full set of peroxisome proteins, indicates that loss of both features in other apicomplexans occurred multiple times. The degradation of purines is thought to play a key role in the unusual relationship between Nephromyces and its host. Transcriptome data confirm previous biochemical results of a functional pathway for the utilization of uric acid as a primary nitrogen source for this unusual apicomplexan.

Bacteriocyte Reprogramming to Cope With Nutritional Stress in a Phloem Sap Feeding Hemipteran, the Pea Aphid Acyrthosiphon pisum

Nutritional symbioses play a central role in the ability of insects to thrive on unbalanced diets and in ensuring their evolutionary success. A genomic model for nutritional symbiosis comprises the hemipteran Acyrthosiphon pisum, and the gamma-3-proteobacterium, Buchnera aphidicola, with genomes encoding highly integrated metabolic pathways. A. pisum feeds exclusively on plant phloem sap, a nutritionally unbalanced diet highly variable in composition, thus raising the question of how this symbiotic system responds to nutritional stress. We addressed this by combining transcriptomic, phenotypic and life history trait analyses to determine the organismal impact of deprivation of tyrosine and phenylalanine. These two aromatic amino acids are essential for aphid development, are synthesized in a metabolic pathway for which the aphid host and the endosymbiont are interdependent, and their concentration can be highly variable in plant phloem sap. We found that this nutritional challenge does not have major phenotypic effects on the pea aphid, except for a limited weight reduction and a 2-day delay in onset of nymph laying. Transcriptomic analyses through aphid development showed a prominent response in bacteriocytes (the core symbiotic tissue which houses the symbionts), but not in gut, thus highlighting the role of bacteriocytes as major modulators of this homeostasis. This response does not involve a direct regulation of tyrosine and phenylalanine biosynthetic pathway and transporter genes. Instead, we observed an extensive transcriptional reprogramming of the bacteriocyte with a rapid down-regulation of genes encoding sugar transporters and genes required for sugar metabolism. Consistently, we observed continued overexpression of the A. pisum homolog of RRAD, a small GTPase implicated in repressing aerobic glycolysis. In addition, we found increased transcription of genes involved in proliferation, cell size control and signaling. We experimentally confirmed the significance of these gene expression changes detecting an increase in bacteriocyte number and cell size in vivo under tyrosine and phenylalanine depletion. Our results support a central role of bacteriocytes in the aphid response to amino acid deprivation: their transcriptional and cellular responses fine-tune host physiology providing the host insect with an effective way to cope with the challenges posed by the variability in composition of phloem sap.

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.

Feminizing Wolbachia influence microbiota composition in the terrestrial isopod Armadillidium vulgare

Wolbachia are widespread heritable endosymbionts of arthropods notorious for their profound effects on host fitness as well as for providing protection against viruses and eukaryotic parasites, indicating that they can interact with other microorganisms sharing the same host environment. Using the terrestrial isopod crustacean Armadillidium vulgare, its highly diverse microbiota (>200 bacterial genera) and its three feminizing Wolbachia strains (wVulC, wVulM, wVulP) as a model system, the present study demonstrates that Wolbachia can even influence the composition of a diverse bacterial community under both laboratory and natural conditions. While host origin is the major determinant of the taxonomic composition of the microbiota in A. vulgare, Wolbachia infection affected both the presence and, more importantly, the abundance of many bacterial taxa within each host population, possibly due to competitive interactions. Moreover, different Wolbachia strains had different impacts on microbiota composition. As such, infection with wVulC affected a higher number of taxa than infection with wVulM, possibly due to intrinsic differences in virulence and titer between these two strains. In conclusion, this study shows that heritable endosymbionts such as Wolbachia can act as biotic factors shaping the microbiota of arthropods, with as yet unknown consequences on host fitness.

A sRNA in a reduced mutualistic symbiont genome regulates its own gene expression

Similar to other nutritional endosymbionts that are obligate for host survival, the mutualistic aphid endosymbiont, Buchnera, has a highly reduced genome with few regulatory elements. Until recently, it was thought that aphid hosts were primarily responsible for regulating their symbiotic relationship. However, we recently revealed that Buchnera displays differential protein regulation, but not mRNA expression. We also identified a number of conserved small RNAs (sRNAs) that are expressed among Buchnera taxa. In this study, we investigate whether differential protein regulation in Buchnera is the result of post-transcriptional gene regulation via sRNAs. We characterize the sRNA profile of two Buchnera life stages: (i) when Buchnera is transitioning from an extracellular proliferating state in aphid embryos and (ii) when Buchnera is in an intracellular nonproliferating state in aphid bacteriocytes (specialized symbiont cells). Overall, we identified 90 differentially expressed sRNAs, 97% of which were upregulated in aphid embryos. Of these sRNAs, the majority were predicted to be involved in the regulation of various metabolic processes, including arginine biosynthesis. Using a heterologous dual expression vector, we reveal for the first time that a Buchnera antisense sRNA can post-transcriptionally interact with its cognate Buchnera coding sequence, carB, a gene involved in arginine biosynthesis. These results corroborate our in vivo RNAseq and proteomic data, where the candidate antisense sRNA carB and the protein CarB are significantly upregulated in aphid embryos. Overall, we demonstrate that Buchnera may regulate gene expression independently from its host by utilizing sRNAs.

The evolutionary development of plant-feeding insects and their nutritional endosymbionts

Herbivorous insects have evolved diverse mechanisms enabling them to feed on plants with suboptimal nutrient availability. Low nutrient availability negatively impacts insect herbivore development and fitness. To overcome this obstacle numerous insect lineages have evolved intimate associations with nutritional endosymbionts. This is especially true for insects that specialize on nitrogen-poor substrates, as these insects are highly dependent on intracellular symbionts to provide nitrogen lacking in their insect host's diet. Emerging evidence in these systems suggest that the symbiont's and/or the insect's biosynthetic pathways are dynamically regulated throughout the insect's development to potentially cope with the insect's changing nutritional demands. In this review, we evaluate the evolutionary development of symbiotic insect cells (bacteriocytes) by comparing and contrasting genes and mechanisms involved in maintaining and regulating the nutritional symbiosis throughout insect development in a diversity of insect herbivore-endosymbiont associations. With new advances in genome sequencing and functional genomics, we evaluate to what extent nutritional symbioses are shaped by (i) the regulation of symbiont titer, (ii) the regulation of insect symbiosis genes, and (iii) the regulation of symbiont genes. We discuss how important these mechanisms are for the biosynthesis of essential amino acids and vitamins across insect life stages in divergent insect-symbiont systems. We conclude by suggesting future directions of research to further elucidate the evolutionary development of bacteriocytes and the impact of these nutritional symbioses on insect-plant interactions.

Genome-Wide Transcriptional Dynamics in the Companion Bacterial Symbionts of the Glassy-Winged Sharpshooter (Cicadellidae: Homalodisca vitripennis) Reveal Differential Gene Expression in Bacteria Occupying Multiple Host Organs

The agricultural pest known as the glassy-winged sharpshooter (GWSS) or Homalodisca vitripennis (Hemiptera: Cicadellidae) harbors two bacterial symbionts, “Candidatus Sulcia muelleri” and “Ca. Baumannia cicadellinicola,” which provide the 10 essential amino acids (EAAs) that are limited in the host plant-sap diet. Although they differ in origin and symbiotic age, both bacteria have experienced extensive genome degradation resulting from their ancient restriction to specialized host organs (bacteriomes) that provide cellular support and ensure vertical transmission. GWSS bacteriomes are of different origins and distinctly colored red and yellow. While Sulcia occupies the yellow bacteriome, Baumannia inhabits both. Aside from genomic predictions, little is currently known about the cellular functions of these bacterial symbionts, particularly whether Baumannia in different bacteriomes perform different roles in the symbiosis. To address these questions, we conducted a replicated, strand-specific RNA-seq experiment to assay global gene expression patterns in Sulcia and Baumannia. Despite differences in genomic capabilities, the symbionts exhibit similar profiles of their most highly expressed genes, including those involved in nutrition synthesis and protein stability (chaperonins dnaK and groESL) that likely aid impaired proteins. Baumannia populations in separate bacteriomes differentially express genes enriched in essential nutrient synthesis, including EAAs (histidine and methionine) and B vitamins (biotin and thiamine). Patterns of differential gene expression further reveal complexity in methionine synthesis. Baumannia’s capability to differentially express genes is unusual, as ancient symbionts lose the capability to independently regulate transcription. Combined with previous microscopy, our results suggest that the GWSS may rely on distinct Baumannia populations for essential nutrition and vertical transmission.

Effects of different temperature regimes on survival of Diaphorina citri and its endosymbiotic bacterial communities

The Asian citrus psyllid, Diaphorina citri, is a major pest of citrus and vector of citrus greening (huanglongbing) in Asian. In our field-collected psyllid samples, we discovered that Fuzhou (China) and Faisalabad (Pakistan), populations harbored an obligate primary endosymbiont Candidatus Carsonella (gen. nov.) with a single species, Candidatus Carsonella ruddii (sp. nov.) and a secondary endosymbiont, Wolbachia surface proteins (WSP) which are intracellular endosymbionts residing in the bacteriomes. Responses of these symbionts to different temperatures were examined and their host survival assessed. Diagnostic PCR assays showed that the endosymbionts infection rates were not significantly reduced in both D. citri populations after 24 h exposure to cold or heat treatments. Although quantitative PCR assays showed significant reduction of WSP relative densities at 40°C for 24 h, a substantial decrease occurred as the exposure duration increased beyond 3 days. Under the same temperature regimes, Ca. C. ruddii density was initially less affected during the first exposure day, but rapidly reduced at 3-5 days compared to WSP. However, the mortality of the psyllids increased rapidly as exposure time to heat treatment increased. The responses of the two symbionts to unfavorable temperature regimes highlight the complex host-symbionts interactions between D. citri and its associated endosymbionts.

Chance and necessity in the genome evolution of endosymbiotic bacteria of insects

An open question in evolutionary biology is how does the selection–drift balance determine the fates of biological interactions. We searched for signatures of selection and drift in genomes of five endosymbiotic bacterial groups known to evolve under strong genetic drift. Although most genes in endosymbiotic bacteria showed evidence of relaxed purifying selection, many genes in these bacteria exhibited stronger selective constraints than their orthologs in free-living bacterial relatives. Remarkably, most of these highly constrained genes had no role in the host–symbiont interactions but were involved in either buffering the deleterious consequences of drift or other host-unrelated functions, suggesting that they have either acquired new roles or their role became more central in endosymbiotic bacteria. Experimental evolution of Escherichia coli under strong genetic drift revealed remarkable similarities in the mutational spectrum, genome reduction patterns and gene losses to endosymbiotic bacteria of insects. Interestingly, the transcriptome of the experimentally evolved lines showed a generalized deregulation of the genome that affected genes encoding proteins involved in mutational buffering, regulation and amino acid biosynthesis, patterns identical to those found in endosymbiotic bacteria. Our results indicate that drift has shaped endosymbiotic associations through a change in the functional landscape of bacterial genes and that the host had only a small role in such a shift.

The regulation of antimicrobial peptide resistance in the transition to insect symbiosis

Many bacteria utilize two-component systems consisting of a sensor kinase and a transcriptional response regulator to detect environmental signals and modulate gene expression for adaptation. The response regulator PhoP and its cognate sensor kinase PhoQ compose a two-component system known for its role in responding to low levels of Mg²⁺ , Ca²⁺ , pH and to the presence of antimicrobial peptides and activating the expression of genes involved in adaptation to host association. Compared with their free-living relatives, mutualistic insect symbiotic bacteria inhabit a static environment where the requirement for sensory functions is expected to be relaxed. The insect symbiont, Sodalis glossinidius, requires PhoP to resist killing by host derived antimicrobial peptides. However, the S. glossinidius PhoQ was found to be insensitive to Mg²⁺ , Ca²⁺ and pH. Here they show that Sodalis praecaptivus, a close non host-associated relative of S. glossinidius, utilizes a magnesium sensing PhoP-PhoQ and an uncharacterized MarR-like transcriptional regulator (Sant_4061) to control antimicrobial peptide resistance in vitro. While the inactivation of phoP, phoQ or Sant_4061 completely retards the growth of S. praecaptivus in the presence of an antimicrobial peptide in vitro, inactivation of both phoP and Sant_4061 is necessary to abrogate growth of this bacterium in an insect host.

Single sample resolution of rare microbial dark matter in a marine invertebrate metagenome

Direct, untargeted sequencing of environmental samples (metagenomics) and de novo genome assembly enable the study of uncultured and phylogenetically divergent organisms. However, separating individual genomes from a mixed community has often relied on the differential-coverage analysis of multiple, deeply sequenced samples. In the metagenomic investigation of the marine bryozoan Bugula neritina, we uncovered seven bacterial genomes associated with a single B. neritina individual that appeared to be transient associates, two of which were unique to one individual and undetectable using certain “universal” 16S rRNA primers and probes. We recovered high quality genome assemblies for several rare instances of “microbial dark matter,” or phylogenetically divergent bacteria lacking genomes in reference databases, from a single tissue sample that was not subjected to any physical or chemical pre-treatment. One of these rare, divergent organisms has a small (593 kbp), poorly annotated genome with low GC content (20.9%) and a 16S rRNA gene with just 65% sequence similarity to the closest reference sequence. Our findings illustrate the importance of sampling strategy and de novo assembly of metagenomic reads to understand the extent and function of bacterial biodiversity.

The Terrestrial Isopod Microbiome: An All-in-One Toolbox for Animal-Microbe Interactions of Ecological Relevance

Bacterial symbionts represent essential drivers of arthropod ecology and evolution, influencing host traits such as nutrition, reproduction, immunity, and speciation. However, the majority of work on arthropod microbiota has been conducted in insects and more studies in non-model species across different ecological niches will be needed to complete our understanding of host–microbiota interactions. In this review, we present terrestrial isopod crustaceans as an emerging model organism to investigate symbiotic associations with potential relevance to ecosystem functioning. Terrestrial isopods comprise a group of crustaceans that have evolved a terrestrial lifestyle and represent keystone species in terrestrial ecosystems, contributing to the decomposition of organic matter and regulating the microbial food web. Since their nutrition is based on plant detritus, it has long been suspected that bacterial symbionts located in the digestive tissues might play an important role in host nutrition via the provisioning of digestive enzymes, thereby enabling the utilization of recalcitrant food compounds (e.g., cellulose or lignins). If this were the case, then (i) the acquisition of these bacteria might have been an important evolutionary prerequisite for the colonization of land by isopods, and (ii) these bacterial symbionts would directly mediate the role of their hosts in ecosystem functioning. Several bacterial symbionts have indeed been discovered in the midgut caeca of terrestrial isopods and some of them might be specific to this group of animals (i.e., Candidatus Hepatoplasma crinochetorum, Candidatus Hepatincola porcellionum, and Rhabdochlamydia porcellionis), while others are well-known intracellular pathogens (Rickettsiella spp.) or reproductive parasites (Wolbachia sp.). Moreover, a recent investigation of the microbiota in Armadillidium vulgare has revealed that this species harbors a highly diverse bacterial community which varies between host populations, suggesting an important share of environmental microbes in the host-associated microbiota. In this review, we synthesize our current knowledge on the terrestrial isopod microbiome and identify future directions to (i) fully understand the functional roles of particular bacteria (both intracellular or intestinal symbionts and environmental gut passengers), and (ii) whether and how the host-associated microbiota could influence the performance of terrestrial isopods as keystone species in soil ecosystems.

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.

Host Plant Determines the Population Size of an Obligate Symbiont (Buchnera aphidicola) in Aphids

Buchnera aphidicolais an obligate endosymbiont that provides aphids with several essential nutrients. Though much is known about aphid-Buchnera interactions, the effect of the host plant on Buchnera population size remains unclear. Here we used quantitative PCR (qPCR) techniques to explore the effects of the host plant on Buchnera densities in the cotton-melon aphid, Aphis gossypii Buchneratiters were significantly higher in populations that had been reared on cucumber for over 10 years than in populations maintained on cotton for a similar length of time. Aphids collected in the wild from hibiscus and zucchini harbored more Buchnera symbionts than those collected from cucumber and cotton. The effect of aphid genotype on the population size of Buchnera depended on the host plant upon which they fed. When aphids from populations maintained on cucumber or cotton were transferred to novel host plants, host survival and Buchnera population size fluctuated markedly for the first two generations before becoming relatively stable in the third and later generations. Host plant extracts from cucumber, pumpkin, zucchini, and cowpea added to artificial diets led to a significant increase in Buchnera titers in the aphids from the population reared on cotton, while plant extracts from cotton and zucchini led to a decrease in Buchnera titers in the aphids reared on cucumber. Gossypol, a secondary metabolite from cotton, suppressed Buchnera populations in populations from both cotton and cucumber, while cucurbitacin from cucurbit plants led to higher densities. Together, the results suggest that host plants influence Buchnera population processes and that this may provide phenotypic plasticity in host plant use for clonal aphids.

Host origin and tissue microhabitat shaping the microbiota of the terrestrial isopod Armadillidium vulgare

We present the first in-depth investigation of the host-associated microbiota of the terrestrial isopod crustacean Armadillidium vulgare. This species is an important decomposer of organic matter in terrestrial ecosystems and a major model organism for arthropod-Wolbachia symbioses due to its well-characterized association with feminizing Wolbachia 16S rRNA gene pyrotags were used to characterize its bacterial microbiota at multiple levels: (i) in individuals from laboratory lineages and field populations and (ii) in various host tissues. This integrative approach allowed us to reveal an unexpectedly high bacterial diversity, placing this species in the same league as termites in terms of symbiotic diversity. Interestingly, both animal groups belong to the same ecological guild in terrestrial ecosystems. While Wolbachia represented the predominant taxon in infected individuals, it was not the only major player. Together, the most abundant taxa represented a large scope of symbiotic interactions, including bacterial pathogens, a reproductive parasite (Wolbachia) and potential nutritional symbionts. Furthermore, we demonstrate that individuals from different populations harboured distinct bacterial communities, indicating a strong link between the host-associated microbiota and environmental bacteria, possibly due to terrestrial isopod nutritional ecology. Overall, this work highlights the need for more studies of host-microbiota interactions and bacterial diversity in non-insect arthropods.

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.

Dynamics of Wolbachia pipientis Gene Expression Across the Drosophila melanogaster Life Cycle

Symbiotic interactions between microbes and their multicellular hosts have manifold biological consequences. To better understand how bacteria maintain symbiotic associations with animal hosts, we analyzed genome-wide gene expression for the endosymbiotic α-proteobacteria Wolbachia pipientis across the entire life cycle of Drosophila melanogaster. We found that the majority of Wolbachia genes are expressed stably across the D. melanogaster life cycle, but that 7.8% of Wolbachia genes exhibit robust stage- or sex-specific expression differences when studied in the whole-organism context. Differentially-expressed Wolbachia genes are typically up-regulated after Drosophila embryogenesis and include many bacterial membrane, secretion system, and ankyrin repeat-containing proteins. Sex-biased genes are often organized as small operons of uncharacterized genes and are mainly up-regulated in adult Drosophila males in an age-dependent manner. We also systematically investigated expression levels of previously-reported candidate genes thought to be involved in host-microbe interaction, including those in the WO-A and WO-B prophages and in the Octomom region, which has been implicated in regulating bacterial titer and pathogenicity. Our work provides comprehensive insight into the developmental dynamics of gene expression for a widespread endosymbiont in its natural host context, and shows that public gene expression data harbor rich resources to probe the functional basis of the Wolbachia-Drosophila symbiosis and annotate the transcriptional outputs of the Wolbachia genome.

Matching the supply of bacterial nutrients to the nutritional demand of the animal host

Various animals derive nutrients from symbiotic microorganisms with much-reduced genomes, but it is unknown whether, and how, the supply of these nutrients is regulated. Here, we demonstrate that the production of essential amino acids (EAAs) by the bacterium Buchnera aphidicola in the pea aphid Acyrthosiphon pisum is elevated when aphids are reared on diets from which that EAA are omitted, demonstrating that Buchnera scale EAA production to host demand. Quantitative proteomics of bacteriocytes (host cells bearing Buchnera) revealed that these metabolic changes are not accompanied by significant change in Buchnera or host proteins, suggesting that EAA production is regulated post-translationally. Bacteriocytes in aphids reared on diet lacking the EAA methionine had elevated concentrations of both methionine and the precursor cystathionine, indicating that methionine production is promoted by precursor supply and is not subject to feedback inhibition by methionine. Furthermore, methionine production by isolated Buchnera increased with increasing cystathionine concentration. We propose that Buchnera metabolism is poised for EAA production at certain maximal rates, and the realized release rate is determined by precursor supply from the host. The incidence of host regulation of symbiont nutritional function via supply of key nutritional inputs in other symbioses remains to be investigated.

Whole Genome Sequence of the Soybean Aphid Endosymbiont Buchnera aphidicola and Genetic Differentiation among Biotype-Specific Strains

Endosymbiosis with microorganisms is common in insects, with more than 10% of species requiring the metabolic capabilities of intracellular bacteria for their nutrient acquisition. Aphids harbor an obligate mutualism with the vertically transferred endosymbiont, Buchnera aphidicola, which produces key nutrients lacking in the aphid's phloem-based diet that are necessary for normal development and reproduction. It is thought that, in some groups of insects, bacterial symbionts may play key roles in biotype evolution against host-plant resistance. The genome of Buchnera has been sequenced in several aphid strains but little genomic data is currently available for the soybean aphid (Aphis glycines), one of the most important pests of soybean in North America. In this study, DNA sequencing was used to assemble and annotate the genome sequence of the Buchnera A. glycines strain and to reconstruct phylogenetic relationships among different strains. In addition, we identified several fixed Buchnera SNPs between Aphis glycines biotypes that were avirulent or virulent to a soybean aphid resistance gene (Rag1). The results of this study describe the genetic and evolutionary relationships of the Buchnera A. glycines strain, and begin to define the roles of an aphid symbiont in host-plant resistance.

Shifts in Buchnera aphidicola density in soybean aphids (Aphis glycines) feeding on virus-infected soybean

Vertically transmitted bacterial symbionts are common in arthropods. Aphids undergo an obligate symbiosis with Buchnera aphidicola, which provides essential amino acids to its host and contributes directly to nymph growth and reproduction. We previously found that newly adult Aphis glycines feeding on soybean infected with the beetle-transmitted Bean pod mottle virus (BPMV) had significantly reduced fecundity. We hypothesized that the reduced fecundity was attributable to detrimental impacts of the virus on the aphid microbiome, namely Buchnera. To test this, mRNA sequencing and quantitative real-time PCR were used to assay Buchnera transcript abundance and titre in A. glycines feeding on Soybean mosaic virus-infected, BPMV-infected, and healthy soybean for up to 14 days. Our results indicated that Buchnera density was lower and ultimately suppressed in aphids feeding on virus-infected soybean. While the decreased Buchnera titre may be associated with reduced aphid fecundity, additional mechanisms are probably involved. The present report begins to describe how interactions among insects, plants, and plant pathogens influence endosymbiont population dynamics.

Signatures of host/symbiont genome coevolution in insect nutritional endosymbioses

The role of symbiosis in bacterial symbiont genome evolution is well understood, yet the ways that symbiosis shapes host genomes or more particularly, host/symbiont genome coevolution in the holobiont is only now being revealed. Here, we identify three coevolutionary signatures that characterize holobiont genomes. The first signature, host/symbiont collaboration, arises when completion of essential pathways requires host/endosymbiont genome complementarity. Metabolic collaboration has evolved numerous times in the pathways of amino acid and vitamin biosynthesis. Here, we highlight collaboration in branched-chain amino acid and pantothenate (vitamin B5) biosynthesis. The second coevolutionary signature is acquisition, referring to the observation that holobiont genomes acquire novel genetic material through various means, including gene duplication, lateral gene transfer from bacteria that are not their current obligate symbionts, and full or partial endosymbiont replacement. The third signature, constraint, introduces the idea that holobiont genome evolution is constrained by the processes governing symbiont genome evolution. In addition, we propose that collaboration is constrained by the expression profile of the cell lineage from which endosymbiont-containing host cells, called bacteriocytes, are derived. In particular, we propose that such differences in bacteriocyte cell lineage may explain differences in patterns of host/endosymbiont metabolic collaboration between the sap-feeding suborders Sternorrhyncha and Auchenorrhynca. Finally, we review recent studies at the frontier of symbiosis research that are applying functional genomic approaches to characterization of the developmental and cellular mechanisms of host/endosymbiont integration, work that heralds a new era in symbiosis research.

Widespread expression of conserved small RNAs in small symbiont genomes

Genome architecture of a microbe markedly changes when it transitions from a free-living lifestyle to an obligate symbiotic association within eukaryotic cells. These symbiont genomes experience numerous rearrangements and massive gene loss, which is expected to radically alter gene regulatory networks compared with those of free-living relatives. As such, it remains unclear whether and how these small symbiont genomes regulate gene expression. Here, using a label-free mass-spec quantification approach we found that differential protein regulation occurs in Buchnera, a model symbiont with a reduced genome, when it transitions between two distinct life stages. However, differential mRNA expression could not be detected between Buchnera life stages, despite the presence of a small number of putative transcriptional regulators. Instead a comparative analysis of small RNA expression profiles among five divergent Buchnera lineages, spanning a variety of Buchnera life stages, reveals 140 novel intergenic and antisense small RNAs and 517 untranslated regions that were significantly expressed, some of which have been conserved for ∼65 million years. In addition, the majority of these small RNAs exhibit both sequence covariation and thermodynamic stability, indicators of a potential structural RNA role. Together, these data suggest that gene regulation at the post-transcriptional level may be important in Buchnera. This is the first study to empirically identify Buchnera small RNAs, and we propose that these novel small RNAs may facilitate post-transcriptional regulation through translational inhibition/activation, and/or transcript stability. Ultimately, post-transcriptional regulation may shape metabolic complementation between Buchnera and its aphid host, thus impacting the animal's ecology and evolution.

Settling down: the genome of Serratia symbiotica from the aphid Cinara tujafilina zooms in on the process of accommodation to a cooperative intracellular life

Particularly interesting cases of mutualistic endosymbioses come from the establishment of co-obligate associations of more than one species of endosymbiotic bacteria. Throughout symbiotic accommodation from a free-living bacterium, passing through a facultative stage and ending as an obligate intracellular one, the symbiont experiences massive genomic losses and phenotypic adjustments. Here, we scrutinized the changes in the coevolution of Serratia symbiotica and Buchnera aphidicola endosymbionts in aphids, paying particular attention to the transformations undergone by S. symbiotica to become an obligate endosymbiont. Although it is already known that S. symbiotica is facultative in Acyrthosiphon pisum, in Cinara cedri it has established a co-obligate endosymbiotic consortium along with B. aphidicola to fulfill the aphid’s nutritional requirements. The state of this association in C. tujafilina, an aphid belonging to the same subfamily (Lachninae) that C. cedri, remained unknown. Here, we report the genome of S. symbiotica strain SCt-VLC from the aphid C. tujafilina. While being phylogenetically and genomically very closely related to the facultative endosymbiont S. symbiotica from the aphid A. pisum, it shows a variety of metabolic, genetic, and architectural features, which point toward this endosymbiont being one step closer to an obligate intracellular one. We also describe in depth the process of genome rearrangements suffered by S. symbiotica and the role mobile elements play in gene inactivations. Finally, we postulate the supply to the host of the essential riboflavin (vitamin B₂) as key to the establishment of S. symbiotica as a co-obligate endosymbiont in the aphids belonging to the subfamily Lachninane.

Comparative analysis of genome sequences from four strains of the Buchnera aphidicola Mp endosymbion of the green peach aphid, Myzus persicae

BACKGROUND

Myzus persicae, the green peach aphid, is a polyphagous herbivore that feeds from hundreds of species of mostly dicot crop plants. Like other phloem-feeding aphids, M. persicae rely on the endosymbiotic bacterium, Buchnera aphidicola (Buchnera Mp), for biosynthesis of essential amino acids and other nutrients that are not sufficiently abundant in their phloem sap diet. Tobacco-specialized M. persicae are typically red and somewhat distinct from other lineages of this species. To determine whether the endosymbiotic bacteria of M. persicae could play a role in tobacco adaptation, we sequenced the Buchnera Mp genomes from two tobacco-adapted and two non-tobacco M. persicae lineages.

RESULTS

With a genome size of 643.5 kb and 579 predicted genes, Buchnera Mp is the largest Buchnera genome sequenced to date. No differences in gene content were found between the four sequenced Buchnera Mp strains. Compared to Buchnera APS from the well-studied pea aphid, Acyrthosiphon pisum, Buchnera Mp has 21 additional genes. These include genes encoding five enzymes required for biosynthesis of the modified nucleoside queosine, the heme pathway enzyme uroporphyrinogen III synthase, and asparaginase. Asparaginase, which is also encoded by the genome of the aphid host, may allow Buchnera Mp to synthesize essential amino acids from asparagine, a relatively abundant phloem amino acid.

CONCLUSIONS

Together our results indicate that the obligate intracellular symbiont Buchnera aphidicola does not contribute to the adaptation of Myzus persicae to feeding on tobacco.

Aphid amino acid transporter regulates glutamine supply to intracellular bacterial symbionts

Endosymbiotic associations have played a major role in evolution. However, the molecular basis for the biochemical interdependence of these associations remains poorly understood. The aphid-Buchnera endosymbiosis provides a powerful system to elucidate how these symbioses are regulated. In aphids, the supply of essential amino acids depends on an ancient nutritional symbiotic association with the gamma-proteobacterium Buchnera aphidicola. Buchnera cells are densely packed in specialized aphid bacteriocyte cells. Here we confirm that five putative amino acid transporters are highly expressed and/or highly enriched in Acyrthosiphon pisum bacteriocyte tissues. When expressed in Xenopus laevis oocytes, two bacteriocyte amino acid transporters displayed significant levels of glutamine uptake, with transporter ACYPI001018, LOC100159667 (named here as Acyrthosiphon pisum glutamine transporter 1, ApGLNT1) functioning as the most active glutamine transporter. Transporter ApGLNT1 has narrow substrate selectivity, with high glutamine and low arginine transport capacity. Notably, ApGLNT1 has high binding affinity for arginine, and arginine acts as a competitive inhibitor for glutamine transport. Using immunocytochemistry, we show that ApGLNT1 is localized predominantly to the bacteriocyte plasma membrane, a location consistent with the transport of glutamine from A. pisum hemolymph to the bacteriocyte cytoplasm. On the basis of functional transport data and localization, we propose a substrate feedback inhibition model in which the accumulation of the essential amino acid arginine in A. pisum hemolymph reduces the transport of the precursor glutamine into bacteriocytes, thereby regulating amino acid biosynthesis in the bacteriocyte. Structural similarities in the arrangement of hosts and symbionts across endosymbiotic systems suggest that substrate feedback inhibition may be mechanistically important in other endosymbioses.