Systemic analysis of the symbiotic function of Buchnera aphidicola, the primary endosymbiont of the pea aphid Acyrthosiphon pisum

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Aphelinus nigritus Howard (Hymenoptera: Aphelinidae) Preference for Sorghum Aphid, Melanaphis sorghi (Theobald, 1904) (Hemiptera: Aphididae), Honeydew Is Stronger in Johnson Grass, Sorghum halepense, Than in Grain Sorghum, Sorghum bicolor

How aphid parasitoids of recent invasive species interact with their hosts can affect the feasibility of biological control. In this study, we focus on a recent invasive pest of US grain sorghum, Sorghum bicolor, the sorghum aphid (SA), Melanaphis sorghi. Understanding this pest's ecology in the grain sorghum agroecosystem is critical to develop effective control strategies. As parasitoids often use aphid honeydew as a sugar resource, and honeydew is known to mediate parasitoid-aphid interactions, we investigated the ability of SA honeydew to retain the parasitoid Aphelinus nigritus. Since SAs in the US have multiple plant hosts, and host-plant diet can modulate parasitoid retention (a major component in host foraging), we measured SA honeydew sugar, organic acid, and amino acid profiles, then assessed via retention time A. nigritus preference for honeydew produced on grain sorghum or Johnson grass, Sorghum halepense. Compared to a water control, A. nigritus spent more time on SA honeydew produced on either host plant. Despite similar honeydew profiles from both plant species, A. nigritus preferred honeydew produced on Johnson grass. Our results suggest the potential for SA honeydew to facilitate augmentation strategies aimed at maintaining A. nigritus on Johnson grass to suppress SAs before grain sorghum is planted.

Evolution and ontogeny of bacteriocytes in insects

The ontogenetic origins of the bacteriocytes, which are cells that harbour bacterial intracellular endosymbionts in multicellular animals, are unknown. During embryonic development, a series of morphological and transcriptional changes determine the fate of distinct cell types. The ontogeny of bacteriocytes is intimately linked with the evolutionary transition of endosymbionts from an extracellular to an intracellular environment, which in turn is linked to the diet of the host insect. Here we review the evolution and development of bacteriocytes in insects. We first classify the endosymbiotic occupants of bacteriocytes, highlighting the complex challenges they pose to the host. Then, we recall the historical account of the discovery of bacteriocytes. We then summarize the molecular interactions between the endosymbiont and the host. In addition, we illustrate the genetic contexts in which the bacteriocytes develop, with examples of the genetic changes in the hosts and endosymbionts, during specific endosymbiotic associations. We finally address the evolutionary origin as well as the putative ontogenetic or developmental source of bacteriocytes in insects.

Aphid BCR4 Structure and Activity Uncover a New Defensin Peptide Superfamily

Aphids (Hemiptera: Aphidoidea) are among the most detrimental insects for agricultural plants, and their management is a great challenge in agronomical research. A new class of proteins, called Bacteriocyte-specific Cysteine-Rich (BCR) peptides, provides an alternative to chemical insecticides for pest control. BCRs were initially identified in the pea aphid Acyrthosiphon pisum. They are small disulfide bond-rich proteins expressed exclusively in aphid bacteriocytes, the insect cells that host intracellular symbiotic bacteria. Here, we show that one of the A. pisum BCRs, BCR4, displays prominent insecticidal activity against the pea aphid, impairing insect survival and nymphal growth, providing evidence for its potential use as a new biopesticide. Our comparative genomics and phylogenetic analyses indicate that BCRs are restricted to the aphid lineage. The 3D structure of BCR4 reveals that this peptide belongs to an as-yet-unknown structural class of peptides and defines a new superfamily of defensins.

Microbiome Structure of the Aphid Myzus persicae (Sulzer) Is Shaped by Different Solanaceae Plant Diets

Myzus persicae (Sulzer) is an important insect pest in agriculture that has a very broad host range. Previous research has shown that the microbiota of insects has implications for their growth, development, and environmental adaptation. So far, there is little detailed knowledge about the factors that influence and shape the microbiota of aphids. In the present study, we aimed to investigate diet-induced changes in the microbiome of M. persicae using high-throughput sequencing of bacterial 16S ribosomal RNA gene fragments in combination with molecular and microbiological experiments. The transfer of aphids to different plants from the Solanaceae family resulted in a substantial decrease in the abundance of the primary symbiont Buchnera. In parallel, a substantial increase in the abundance of Pseudomonas was observed; it accounted for up to 69.4% of the bacterial community in M. persicae guts and the attached bacteriocytes. In addition, we observed negative effects on aphid population dynamics when they were transferred to pepper plants (Capsicum annuum L.). The microbiome of this treatment group showed a significantly lower increase in the abundance of Pseudomonas when compared with the other Solanaceae plant diets, which might be related to the adaptability of the host to this diet. Molecular quantifications of bacterial genera that were substantially affected by the different diets were implemented as an additional verification of the microbiome-based observations. Complementary experiments with bacteria isolated from aphids that were fed with different plants indicated that nicotine-tolerant strains occur in Solanaceae-fed specimens, but they were not restricted to them. Overall, our mechanistic approach conducted under controlled conditions provided strong indications that the aphid microbiome shows responses to different plant diets. This knowledge could be used in the future to develop environmentally friendly methods for the control of insect pests in agriculture.

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.

Aphid infestation differently affects the defences of nitrate-fed and nitrogen-fixing Medicago truncatula and alters symbiotic nitrogen fixation

Legumes can meet their nitrogen requirements through root nodule symbiosis, which could also trigger plant systemic resistance against pests. The pea aphid Acyrthosiphon pisum, a legume pest, can harbour different facultative symbionts (FS) influencing various traits of their hosts. It is therefore worth determining if and how the symbionts of the plant and the aphid modulate their interaction. We used different pea aphid lines without FS or with a single one (Hamiltonella defensa, Regiella insecticola, Serratia symbiotica) to infest Medicago truncatula plants inoculated with Sinorhizobium meliloti (symbiotic nitrogen fixation, SNF) or supplemented with nitrate (non-inoculated, NI). The growth of SNF and NI plants was reduced by aphid infestation, while aphid weight (but not survival) was lowered on SNF compared to NI plants. Aphids strongly affected the plant nitrogen fixation depending on their symbiotic status, suggesting indirect relationships between aphid- and plant-associated microbes. Finally, all aphid lines triggered expression of Pathogenesis-Related Protein 1 (PR1) and Proteinase Inhibitor (PI), respective markers for salicylic and jasmonic pathways, in SNF plants, compared to only PR1 in NI plants. We demonstrate that the plant symbiotic status influences plant-aphid interactions while that of the aphid can modulate the amplitude of the plant's defence response.

Characterization of the bacterial communities of psyllids associated with Rutaceae in Bhutan by high throughput sequencing

Background

Several plant-pathogenic bacteria are transmitted by insect vector species that often also act as hosts. In this interface, these bacteria encounter plant endophytic, insect endosymbiotic and other microbes. Here, we used high throughput sequencing to examine the bacterial communities of five different psyllids associated with citrus and related plants of Rutaceae in Bhutan: Diaphorina citri, Diaphorina communis, Cornopsylla rotundiconis, Cacopsylla heterogena and an unidentified Cacopsylla sp.

Results

The microbiomes of the psyllids largely comprised their obligate P-endosymbiont ‘Candidatus Carsonella ruddii’, and one or two S-endosymbionts that are fixed and specific to each lineage. In addition, all contained Wolbachia strains; the Bhutanese accessions of D. citri were dominated by a Wolbachia strain first found in American isolates of D. citri, while D. communis accessions were dominated by the Wolbachia strain, wDi, first detected in D. citri from China. The S-endosymbionts from the five psyllids grouped with those from other psyllid taxa; all D. citri and D. communis individuals contained sequences matching ‘Candidatus Profftella armatura’ that has previously only been reported from other Diaphorina species, and the remaining psyllid species contained OTUs related to unclassified Enterobacteriaceae. The plant pathogenic ‘Candidatus Liberibacter asiaticus’ was found in D. citri but not in D. communis. Furthermore, an unidentified ‘Candidatus Liberibacter sp.’ occurred at low abundance in both Co. rotundiconis and the unidentified Cacopsylla sp. sampled from Zanthoxylum sp.; the status of this new liberibacter as a plant pathogen and its potential plant hosts are currently unknown. The bacterial communities of Co. rotundiconis also contained a range of OTUs with similarities to bacteria previously found in samples taken from various environmental sources.

Conclusions

The bacterial microbiota detected in these Bhutanese psyllids support the trends that have been seen in previous studies: psyllids have microbiomes largely comprising their obligate P-endosymbiont and one or two S-endosymbionts. In addition, the association with plant pathogens has been demonstrated, with the detection of liberibacters in a known host, D. citri, and identification of a putative new species of liberibacter in Co. rotundiconis and Cacopsylla sp.

Parallel Evolution in the Integration of a Co-obligate Aphid Symbiosis

Insects evolve dependence-often extreme-on microbes for nutrition. This includes cases in which insects harbor multiple endosymbionts that function collectively as a metabolic unit [1-5]. How do these dependences originate [6], and is there a predictable sequence of events leading to the integration of new symbionts? While co-obligate symbioses, in which hosts rely on multiple nutrient-provisioning symbionts, have evolved numerous times across sap-feeding insects, there is only one known case in aphids, involving Buchnera aphidicola and Serratia symbiotica in the Lachninae subfamily [7-9]. Here, we identify three additional independent transitions to the same co-obligate symbiosis in different aphids. Comparing recent and ancient associations allow us to investigate intermediate stages of metabolic and anatomical integration of Serratia. We find that these uniquely replicated evolutionary events support the idea that co-obligate associations initiate in a predictable manner-through parallel evolutionary processes. Specifically, we show how the repeated losses of the riboflavin and peptidoglycan pathways in Buchnera lead to dependence on Serratia. We then provide evidence of a stepwise process of symbiont integration, whereby dependence evolves first. Then, essential amino acid pathways are lost (at ∼30-60 mya), which coincides with the increased anatomical integration of the companion symbiont. Finally, we demonstrate that dependence can evolve ahead of specialized structures (e.g., bacteriocytes), and in one case with no direct nutritional basis. More generally, our results suggest the energetic costs of synthesizing nutrients may provide a unified explanation for the sequence of gene losses that occur during the evolution of co-obligate symbiosis.

Tissue- and Population-Level Microbiome Analysis of the Wasp Spider Argiope bruennichi Identified a Novel Dominant Bacterial Symbiont

Many ecological and evolutionary processes in animals depend upon microbial symbioses. In spiders, the role of the microbiome in these processes remains mostly unknown. We compared the microbiome between populations, individuals, and tissue types of a range-expanding spider, using 16S rRNA gene sequencing. Our study is one of the first to go beyond targeting known endosymbionts in spiders and characterizes the total microbiome across different body compartments (leg, prosoma, hemolymph, book lungs, ovaries, silk glands, midgut, and fecal pellets). Overall, the microbiome differed significantly between populations and individuals, but not between tissue types. The microbiome of the wasp spider Argiope bruennichi features a novel dominant bacterial symbiont, which is abundant in every tissue type in spiders from geographically distinct populations and that is also present in offspring. The novel symbiont is affiliated with the Tenericutes, but has low sequence identity (<85%) to all previously named taxa, suggesting that the novel symbiont represents a new bacterial clade. Its presence in offspring implies that it is vertically transmitted. Our results shed light on the processes that shape microbiome differentiation in this species and raise several questions about the implications of the novel dominant bacterial symbiont on the biology and evolution of its host.

The reproductive system of the male and oviparous female of a model organism-the pea aphid, Acyrthosiphon pisum (Hemiptera, Aphididae)

The structure of the reproductive system of the sexual generation-males and oviparous females-of the pea aphid Acyrthosiphon pisum (Harris) (Hemiptera, Aphididae), a serious pest of cultivated plants of Fabaceae, was investigated. For the first time we describe the morphology, histology and ultrastructure of the reproductive system in both morphs of the sexual generation of aphids within one species, using light and fluorescent microscopy, as well as transmission and scanning electron microscopy. The results revealed that males have testes composed of three follicles fused by the upper ends of the vasa efferentia, the vasa deferentia run independently, the accessory glands are asymmetric and the ejaculatory duct shortened. Oviparous females have ovaries composed of seven ovarioles each. The lateral oviducts join to a short common oviduct connected with the unpaired spermatheca and paired accessory glands. Yolky eggs with an aggregation of symbiotic bacteria at the posterior pole are produced. Histologically, the components of genital tracts are broadly similar: the epithelial cells of the walls of the vasa deferentia and accessory glands of the male and oviparous female have secretory functions which correlate with the age of the studied morphs. We also found symbiotic bacteria within the vasa deferentia epithelial cells in males and within the cells of the lateral oviducts of females. Because the pea aphid is listed among the 14 species that are of the greatest economic importance, our results will be useful for managing aphid populations, protecting plants and ensuring global food security.

Reconstructing the functions of endosymbiotic Mollicutes in fungus-growing ants

Mollicutes, a widespread class of bacteria associated with animals and plants, were recently identified as abundant abdominal endosymbionts in healthy workers of attine fungus-farming leaf-cutting ants. We obtained draft genomes of the two most common strains harbored by Panamanian fungus-growing ants. Reconstructions of their functional significance showed that they are independently acquired symbionts, most likely to decompose excess arginine consistent with the farmed fungal cultivars providing this nitrogen-rich amino-acid in variable quantities. Across the attine lineages, the relative abundances of the two Mollicutes strains are associated with the substrate types that foraging workers offer to fungus gardens. One of the symbionts is specific to the leaf-cutting ants and has special genomic machinery to catabolize citrate/glucose into acetate, which appears to deliver direct metabolic energy to the ant workers. Unlike other Mollicutes associated with insect hosts, both attine ant strains have complete phage-defense systems, underlining that they are actively maintained as mutualistic symbionts.

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.

The Origin of Mucosal Immunity: Lessons from the Holobiont Hydra

Historically, mucosal immunity—i.e., the portion of the immune system that protects an organism’s various mucous membranes from invasion by potentially pathogenic microbes—has been studied in single-cell epithelia in the gastrointestinal and upper respiratory tracts of vertebrates. Phylogenetically, mucosal surfaces appeared for the first time about 560 million years ago in members of the phylum Cnidaria. There are remarkable similarities and shared functions of mucosal immunity in vertebrates and innate immunity in cnidarians, such as Hydra species. Here, we propose a common origin for both systems and review observations that indicate that the ultimately simple holobiont Hydra provides both a new perspective on the relationship between bacteria and animal cells and a new prism for viewing the emergence and evolution of epithelial tissue-based innate immunity. In addition, recent breakthroughs in our understanding of immune responses in Hydra polyps reared under defined short-term gnotobiotic conditions open up the potential of Hydra as an animal research model for the study of common mucosal disorders.

Direct flow cytometry measurements reveal a fine-tuning of symbiotic cell dynamics according to the host developmental needs in aphid symbiosis

Endosymbiotic associations constitute a driving force in the ecological and evolutionary diversification of metazoan organisms. Little is known about whether and how symbiotic cells are coordinated according to host physiology. Here, we use the nutritional symbiosis between the insect pest, Acyrthosiphon pisum, and its obligate symbiont, Buchnera aphidicola, as a model system. We have developed a novel approach for unculturable bacteria, based on flow cytometry, and used this method to estimate the absolute numbers of symbionts at key stages of aphid life. The endosymbiont population increases exponentially throughout nymphal development, showing a growing rate which has never been characterized by indirect molecular techniques. Using histology and imaging techniques, we have shown that the endosymbiont-bearing cells (bacteriocytes) increase significantly in number and size during the nymphal development, and clustering in the insect abdomen. Once adulthood is reached and the laying period has begun, the dynamics of symbiont and host cells is reversed: the number of endosymbionts decreases progressively and the bacteriocyte structure degenerates during insect aging. In summary, these results show a coordination of the cellular dynamics between bacteriocytes and primary symbionts and reveal a fine-tuning of aphid symbiotic cells to the nutritional demand imposed by the host physiology throughout development.

A metagenomic approach from aphid's hemolymph sheds light on the potential roles of co-existing endosymbionts

BACKGROUND

Aphids are known to live in symbiosis with specific bacteria, called endosymbionts which can be classified as obligate or accessory. Buchnera aphidicola is generally the only obligatory symbiont present in aphids, supplying essential nutrients that are missing in the plants phloem to its host. Pentalonia nigronervosa is the main vector of the banana bunchy top virus, one of the most damageable viruses in banana. This aphid is carrying two symbionts: B. aphidicola (BPn) and Wolbachia sp. (wPn). The high occurrence of Wolbachia in the banana aphid raises questions about the role it plays in this insect. The goal of this study was to go further in the understanding of the role played by the two symbionts in P. nigronervosa. To do so, microinjection tests were made to see the effect of wPn elimination on the host, and then, high-throughput sequencing of the haemolymph was used to analyze the gene content of the symbionts.

RESULTS

We observed that the elimination of wPn systematically led to the death of aphids, suggesting that the bacterium could play a mutualistic role. In addition, we identify and annotate 587 and 250 genes for wPn and BPn, respectively, through high-throughput sequencing. Analysis of these genes suggests that the two bacteria are working together for the production of several essential nutrients. The most striking cases are for lysin and riboflavin which are usually provided by B. aphidicola alone to the host. In the banana aphid, the genes involved in the production pathways of these metabolites are shared between the two bacteria making them both essential for the survival of the aphid host.

CONCLUSIONS

Our results suggest that a co-obligatory symbiosis between B. aphidicola and Wolbachia occurs in the banana aphid, the two bacteria acting together to supply essential nutrients to the host. This is, to our knowledge, the first time Wolbachia is reported to play an essential role in aphids.

Adenine methylation may contribute to endosymbiont selection in a clonal aphid population

BACKGROUND

The pea aphid Acyrthosiphon pisum has two modes of reproduction: parthenogenetic during the spring and summer and sexual in autumn. This ability to alternate between reproductive modes and the emergence of clonal populations under favorable conditions make this organism an interesting model for genetic and epigenetic studies. The pea aphid hosts different types of endosymbiotic bacteria within bacteriocytes which help the aphids survive and adapt to new environmental conditions and habitats. The obligate endosymbiont Buchnera aphidicola has a drastically reduced and stable genome, whereas facultative endosymbionts such as Regiella insecticola have large and dynamic genomes due to phages, mobile elements and high levels of genetic recombination. In previous work, selection toward cold adaptation resulted in the appearance of parthenogenetic A. pisum individuals characterized by heavier weights and remarkable green pigmentation.

RESULTS

Six adenine-methylated DNA fragments were isolated from genomic DNA (gDNA) extracted from the cold-induced green variant of A. pisum using deoxyadenosine methylase (Dam) by digesting the gDNA with the restriction enzymes DpnI and DpnII, which recognize the methylated and unmethylated GATC sites, respectively. The six resultant fragments did not match any sequence in the A. pisum or Buchnera genomes, implying that they came from facultative endosymbionts. The A1 fragment encoding a putative transposase and the A6 fragment encoding a putative helicase were selected for further comparison between the two A. pisum variants (green and orange) based on Dam analysis followed by PCR amplification. An association between adenine methylation and the two A. pisum variants was demonstrated by higher adenine methylation levels on both genes in the green variant as compared to the orange one.

CONCLUSION

Temperature selection may affect the secondary endosymbiont and the sensitive Dam involved in the survival and adaptation of aphids to cold temperatures. There is a high degree of adenine methylation at the GATC sites of the endosymbiont genes at 8°C, an effect that disappears at 22°C. We suggest that endosymbionts can be modified or selected to increase host fitness under unfavorable climatic conditions, and that the phenotype of the newly adapted aphids can be inherited.

New insight into the RNA interference response against cathepsin-L gene in the pea aphid, Acyrthosiphon pisum: molting or gut phenotypes specifically induced by injection or feeding treatments

RNA interference (RNAi) has been widely and successfully used for gene inactivation in insects, including aphids, where dsRNA administration can be performed either by feeding or microinjection. However, several aspects related to the aphid response to RNAi, as well as the influence of the administration method on tissue response, or the mixed success to observe phenotypes specific to the gene targeted, are still unclear in this insect group. In the present study, we made the first direct comparison of two administration methods (injection or feeding) for delivery of dsRNA targeting the cathepsin-L gene in the pea aphid, Acyrthosiphon pisum. In order to maximize the possibility of discovering specific phenotypes, the effect of the treatment was analyzed in single individual aphids at the level of five body compartments: the bacteriocytes, the gut, the embryonic chains, the head and the remaining body carcass. Our analysis revealed that gene expression knockdown effect in each single body compartment was dependent on the administration method used, and allowed us to discover new functions for the cathepsin-L gene in aphids. Injection of cathepsin-L dsRNA was much more effective on carcass and head, inducing body morphology alterations, and suggesting a novel role of this gene in the molting of these insects. Administration by feeding provoked cathepsin-L knockdown in the gut and specific gut epithelial cell alteration, therefore allowing a better characterization of tissue specific role of this gene in aphids.

Tyrosine pathway regulation is host-mediated in the pea aphid symbiosis during late embryonic and early larval development

BACKGROUND

Nutritional symbioses play a central role in insects' adaptation to specialized diets and in their evolutionary success. The obligatory symbiosis between the pea aphid, Acyrthosiphon pisum, and the bacterium, Buchnera aphidicola, is no exception as it enables this important agricultural pest insect to develop on a diet exclusively based on plant phloem sap. The symbiotic bacteria provide the host with essential amino acids lacking in its diet but necessary for the rapid embryonic growth seen in the parthenogenetic viviparous reproduction of aphids. The aphid furnishes, in exchange, non-essential amino acids and other important metabolites. Understanding the regulations acting on this integrated metabolic system during the development of this insect is essential in elucidating aphid biology.

RESULTS

We used a microarray-based approach to analyse gene expression in the late embryonic and the early larval stages of the pea aphid, characterizing, for the first time, the transcriptional profiles in these developmental phases. Our analyses allowed us to identify key genes in the phenylalanine, tyrosine and dopamine pathways and we identified ACYPI004243, one of the four genes encoding for the aspartate transaminase (E.C. 2.6.1.1), as specifically regulated during development. Indeed, the tyrosine biosynthetic pathway is crucial for the symbiotic metabolism as it is shared between the two partners, all the precursors being produced by B. aphidicola. Our microarray data are supported by HPLC amino acid analyses demonstrating an accumulation of tyrosine at the same developmental stages, with an up-regulation of the tyrosine biosynthetic genes. Tyrosine is also essential for the synthesis of cuticular proteins and it is an important precursor for cuticle maturation: together with the up-regulation of tyrosine biosynthesis, we observed an up-regulation of cuticular genes expression. We were also able to identify some amino acid transporter genes which are essential for the switch over to the late embryonic stages in pea aphid development.

CONCLUSIONS

Our data show that, in the development of A. pisum, a specific host gene set regulates the biosynthetic pathways of amino acids, demonstrating how the regulation of gene expression enables an insect to control the production of metabolites crucial for its own development and symbiotic metabolism.

Genomic analysis of the regulatory elements and links with intrinsic DNA structural properties in the shrunken genome of Buchnera

Background

Buchnera aphidicola is an obligate symbiotic bacterium, associated with most of the aphididae, whose genome has drastically shrunk during intracellular evolution. Gene regulation in Buchnera has been a matter of controversy in recent years as the combination of genomic information with the experimental results has been contradictory, refuting or arguing in favour of a functional and responsive transcription regulation in Buchnera.

The goal of this study was to describe the gene transcription regulation capabilities of Buchnera based on the inventory of cis- and trans-regulators encoded in the genomes of five strains from different aphids (Acyrthosiphon pisum, Schizaphis graminum, Baizongia pistacea, Cinara cedri and Cinara tujafilina), as well as on the characterisation of some intrinsic structural properties of the DNA molecule in these bacteria.

Results

Interaction graph analysis shows that gene neighbourhoods are conserved between E. coli and Buchnera in structures called transcriptons, interactons and metabolons, indicating that selective pressures have acted on the evolution of transcriptional, protein-protein interaction and metabolic networks in Buchnera. The transcriptional regulatory network in Buchnera is composed of a few general DNA-topological regulators (Nucleoid Associated Proteins and topoisomerases), with the quasi-absence of any specific ones (except for multifunctional enzymes with a known gene expression regulatory role in Escherichia coli, such as AlaS, PepA and BolA, and the uncharacterized hypothetical regulators YchA and YrbA). The relative positioning of regulatory genes along the chromosome of Buchnera seems to have conserved its ancestral state, despite the genome erosion. Sigma-70 promoters with canonical thermodynamic sequence profiles were detected upstream of about 94% of the CDS of Buchnera in the different aphids. Based on Stress-Induced Duplex Destabilization (SIDD) measurements, unstable σ⁷⁰ promoters were found specifically associated with the regulator and transporter genes.

Conclusions

This genomic analysis provides supporting evidence of a selection of functional regulatory structures and it has enabled us to propose hypotheses concerning possible links between these regulatory elements and the DNA-topology (i.e., supercoiling, curvature, flexibility and base-pair stability) in the regulation of gene expression in the shrunken genome of Buchnera.

Knowing your friends: invertebrate innate immunity fosters beneficial bacterial symbioses

The innate immune system is present in all animals and is a crucial first line of defence against pathogens. However, animals also harbour large numbers of beneficial microorganisms that can be housed in the digestive tract, in specialized organs or on tissue surfaces. Although invertebrates lack conventional antibody-based immunity, they are capable of eliminating pathogens and, perhaps more importantly, discriminating them from other microorganisms. This Review examines the interactions between the innate immune systems of several model invertebrates and the symbionts of these organisms, and addresses the central question of how these long-lived and specific associations are established and maintained.

A genomic reappraisal of symbiotic function in the aphid/Buchnera symbiosis: reduced transporter sets and variable membrane organisations

Buchnera aphidicola is an obligate symbiotic bacterium that sustains the physiology of aphids by complementing their exclusive phloem sap diet. In this study, we reappraised the transport function of different Buchnera strains, from the aphids Acyrthosiphon pisum, Schizaphis graminum, Baizongia pistaciae and Cinara cedri, using the re-annotation of their transmembrane proteins coupled with an exploration of their metabolic networks. Although metabolic analyses revealed high interdependencies between the host and the bacteria, we demonstrate here that transport in Buchnera is assured by low transporter diversity, when compared to free-living bacteria, being mostly based on a few general transporters, some of which probably have lost their substrate specificity. Moreover, in the four strains studied, an astonishing lack of inner-membrane importers was observed. In Buchnera, the transport function has been shaped by the distinct selective constraints occurring in the Aphididae lineages. Buchnera from A. pisum and S. graminum have a three-membraned system and similar sets of transporters corresponding to most compound classes. Transmission electronic microscopic observations and confocal microscopic analysis of intracellular pH fields revealed that Buchnera does not show any of the typical structures and properties observed in integrated organelles. Buchnera from B. pistaciae seem to possess a unique double membrane system and has, accordingly, lost all of its outer-membrane integral proteins. Lastly, Buchnera from C. cedri revealed an extremely poor repertoire of transporters, with almost no ATP-driven active transport left, despite the clear persistence of the ancestral three-membraned system.

Detection and characterization of Wolbachia infections in natural populations of aphids: is the hidden diversity fully unraveled?

Aphids are a serious threat to agriculture, despite being a rather small group of insects. The about 4,000 species worldwide engage in highly interesting and complex relationships with their microbial fauna. One of the key symbionts in arthropods is Wolbachia, an α-Proteobacterium implicated in many important biological processes and believed to be a potential tool for biological control. Aphids were thought not to harbour Wolbachia; however, current data suggest that its presence in aphids has been missed, probably due to the low titre of the infection and/or to the high divergence of the Wolbachia strains of aphids. The goal of the present study is to map the Wolbachia infection status of natural aphids populations, along with the characterization of the detected Wolbachia strains. Out of 425 samples from Spain, Portugal, Greece, Israel and Iran, 37 were found to be infected. Our results, based mainly on 16S rRNA gene sequencing, indicate the presence of two new Wolbachia supergroups prevailing in aphids, along with some strains belonging either to supergroup B or to supergroup A.

Genome economization in the endosymbiont of the wood roach Cryptocercus punctulatus due to drastic loss of amino acid synthesis capabilities

Cockroaches (Blattaria: Dictyoptera) harbor the endosymbiont Blattabacterium sp. in their abdominal fat body. This endosymbiont is involved in nitrogen recycling and amino acid provision to its host. In this study, the genome of Blattabacterium sp. of Cryptocercus punctulatus (BCpu) was sequenced and compared with those of the symbionts of Blattella germanica and Periplaneta americana, BBge and BPam, respectively. The BCpu genome consists of a chromosome of 605.7 kb and a plasmid of 3.8 kb and is therefore approximately 31 kb smaller than the other two aforementioned genomes. The size reduction is due to the loss of 55 genes, 23 of which belong to biosynthetic pathways for amino acids. The pathways for the production of tryptophan, leucine, isoleucine/threonine/valine, methionine, and cysteine have been completely lost. Additionally, the genes for the enzymes catalyzing the last steps of arginine and lysine biosynthesis, argH and lysA, were found to be missing and pseudogenized, respectively. These gene losses render BCpu auxotrophic for nine amino acids more than those corresponding to BBge and BPam. BCpu has also lost capacities for sulfate reduction, production of heme groups, as well as genes for several other unlinked metabolic processes, and genes present in BBge and BPam in duplicates. Amino acids and cofactors that are not synthesized by BCpu are either produced in abundance by hindgut microbiota or are provisioned via a copious diet of dampwood colonized by putrefying microbiota, supplying host and Blattabacterium symbiont with the necessary nutrients and thus permitting genome economization of BCpu.

Metabolic modeling of endosymbiont genome reduction on a temporal scale

This study explores the order in which individual metabolic genes are lost in an in silico evolutionary process leading from the metabolic network of Eschericia coli to that of the genome-reduced endosymbiont Buchnera aphidicola.

Simulating the reductive evolutionary process under several growth conditions, a remarkable correlation between in silico and phylogenetically reconstructed gene loss time is obtained.

A gene's k-robustness (its depth of backups) is prime determinant of its loss time.

In silico gene loss time is a better predictor of their actual loss times than genomic features and network properties.

Simulating the reductive evolutionary process by the loss of large blocks followed by single-gene deletions, as known to occur in evolution, yields a remarkable correspondence with the phylogenetic reconstruction and the block loss reported in the literature.

An open fundamental challenge in Systems Biology is whether a genome-scale model can predict patterns of genome evolution by realistically accounting for the associated biochemical constraints. In this study, we explore the order in which individual genes are lost in an in silico evolutionary process, leading from the metabolic network of Eschericia coli to that of the endosymbiont Buchnera aphidicola.

To evaluate the in silico gene loss time, we repeated the reductive evolutionary process introduced by Pál et al (2006), denoting the in silico deletion time of a gene in a single run of the reductive evolutionary process as the number of genes deleted before its own deletion occurred. By comparing the in silico evaluations of the gene loss time to that obtained by a phylogenetic reconstruction (Figure 1), we could evaluate the ability of an in silico process to predict temporal patterns of genome reduction. Applying this procedure on a literature-based viable media, we obtained a mean Spearman's correlation of 0.46 (53% of the maximal correlation, empirical P-value <9.9e−4) between in silico and phylogenetically reconstructed loss times. In order to provide an upper bound on evolutionary necessity stemming from metabolic constraints, we searched the space of potential growth media and biomass functions via a simulated annealing search algorithm aimed at identifying an environment/biomass function that maximizes the target correlation between in silico and reconstructed loss times. Simulating the reductive evolutionary process under the growth conditions and biomass function obtained in this process, we managed to improve the correlation between in silico and reconstructed loss times to a mean Spearman's correlation of 0.54 (63% of the maximal correlation, empirical P-value <9.9e−4, Figure 3).

Examining the dependency of the predicted loss time of each gene on its intrinsic network-level properties we find a very strong inverse Spearman's correlation of −0.84 (empirical P-value <9.9e−4) between the order of gene loss predicted in silico and the k-robustness levels of the genes, the latter denoting the depth of their functional backups in the network (Deutscher et al, 2006). Moreover, in order to examine whether the relative loss time of a gene is influenced by its functional dependencies with other genes, we performed a flux-coupling analysis and identified pairs of reactions whose activities asymmetrically depend on each other, i.e., are directionally coupled (Burgard et al, 2004). We find that genes encoding reactions whose activity is needed for activating the other reaction (and not vice versa) have a tendency to be lost later, as one would expect (binomial P-value <1e−14).

To assess the scale of these results, we examined as a control how well genomic features and network properties predict the phylogenetically reconstructed gene loss times. We examined the dependency of the latter on several factors that are known be inversely correlated with the propensity of a gene to be lost (Brinza et al, 2009; Delmotte et al, 2006; Tamames et al, 2007), including the genes' mRNA levels, tAI values (Covert et al, 2004; Reis et al, 2004; Sharp and Li, 1987; Tuller et al, 2010a) and the number of partners the gene products have in a protein–protein interaction network. Remarkably, these genomic features yield considerably lower Spearman's correlation than that obtained by the in silico simulations. Moreover, multiply regressing the loss times from the phylogenetic reconstruction on the in silico gene loss time predictions and the genomic and network variables, we found that the (normalized) coefficient of the in silico predictions in the regression is much higher than those of the genomic features, further testifying to the considerable independent predictive power of the metabolic model.

Finally, simulating the evolutionary process as large block deletions at first followed by single-gene deletions as is thought to occur in evolution (Moran and Mira, 2001; van Ham et al, 2003), a remarkable correspondence with the phylogenetic reconstruction was found. Namely, we find that after a certain amount of genes are deleted from the genome, no further block deletions can occur due to the increasing density of essential genes. Notably, the maximum amount of genes that can be deleted in blocks (i.e., until no more blocks can be deleted) corresponds to the number of genes appearing in our phylogenetic reconstruction from the LCA (last common ancestor of Buchnera and E. coli) to the LCSA (last common symbiotic ancestor, nodes 1–3 in Figure 1A), as described in the literature.

A fundamental challenge in Systems Biology is whether a cell-scale metabolic model can predict patterns of genome evolution by realistically accounting for associated biochemical constraints. Here, we study the order in which genes are lost in an in silico evolutionary process, leading from the metabolic network of Eschericia coli to that of the endosymbiont Buchnera aphidicola. We examine how this order correlates with the order by which the genes were actually lost, as estimated from a phylogenetic reconstruction. By optimizing this correlation across the space of potential growth and biomass conditions, we compute an upper bound estimate on the model's prediction accuracy (R=0.54). The model's network-based predictive ability outperforms predictions obtained using genomic features of individual genes, reflecting the effect of selection imposed by metabolic stoichiometric constraints. Thus, while the timing of gene loss might be expected to be a completely stochastic evolutionary process, remarkably, we find that metabolic considerations, on their own, make a marked 40% contribution to determining when such losses occur.

Multimodal dynamic response of the Buchnera aphidicola pLeu plasmid to variations in leucine demand of its host, the pea aphid Acyrthosiphon pisum

Aphids, important agricultural pests, can grow and reproduce thanks to their intimate symbiosis with the γ-proteobacterium Buchnera aphidicola that furnishes them with essential amino acids lacking in their phloem sap diet. To study how B. aphidicola, with its reduced genome containing very few transcriptional regulators, responds to variations in the metabolic requirements of its host, we concentrated on the leucine metabolic pathway. We show that leucine is a limiting factor for aphid growth and it displays a stimulatory feeding effect. Our metabolic analyses demonstrate that symbiotic aphids are able to respond to leucine starvation or excess by modulating the neosynthesis of this amino acid. At a molecular level, this response involves an early important transcriptional regulation (after 12 h of treatment) followed by a moderate change in the pLeu plasmid copy number. Both responses are no longer apparent after 7 days of treatment. These experimental data are discussed in the light of a re-annotation of the pLeu plasmid regulatory elements. Taken together, our data show that the response of B. aphidicola to the leucine demand of its host is multimodal and dynamically regulated, providing new insights concerning the genetic regulation capabilities of this bacterium in relation to its symbiotic functions.

Metaorganisms as the new frontier

Because it appears that almost all organisms are part of an interdependent metaorganism, an understanding of the underlying host-microbe species associations, and of evolution and molecular underpinnings, has become the new frontier in zoology. The availability of novel high-throughput sequencing methods, together with the conceptual understanding that advances mostly originate at the intersection of traditional disciplinary boundaries, enable biologists to dissect the mechanisms that control the interdependent associations of species. In this review article, we outline some of the issues in inter-species interactions, present two case studies illuminating the necessity of interfacial research when addressing complex and fundamental zoological problems, and show that an interdisciplinary approach that seeks to understand co-evolved multi-species relationships will connect genomes, phenotypes, ecosystems and the evolutionary forces that have shaped them. We hope that this article inspires other collaborations of a similar nature on the diverse landscape commonly referred to as "zoology".

Genomic revelations of a mutualism: the pea aphid and its obligate bacterial symbiont

The symbiosis of the pea aphid Acyrthosphion pisum with the bacterium Buchnera aphidicola APS represents the best-studied insect obligate symbiosis. Here we present a refined picture of this symbiosis by linking pre-genomic observations to new genomic data that includes the complete genomes of the eukaryotic and prokaryotic symbiotic partners. In doing so, we address four issues central to understanding the patterns and processes operating at the A. pisum/Buchnera APS interface. These four issues include: (1) lateral gene transfer, (2) host immunity, (3) symbiotic metabolism, and (4) regulation.

Horizontal gene transfer between bacteria and animals

Horizontal gene transfer is increasingly described between bacteria and animals. Such transfers that are vertically inherited have the potential to influence the evolution of animals. One classic example is the transfer of DNA from mitochondria and chloroplasts to the nucleus after the acquisition of these organelles by eukaryotes. Even today, many of the described instances of bacteria-to-animal transfer occur as part of intimate relationships such as those of endosymbionts and their invertebrate hosts, particularly insects and nematodes, while numerous transfers are also found in asexual animals. Both of these observations are consistent with modern evolutionary theory, in particular the serial endosymbiotic theory and Muller's ratchet. Although it is tempting to suggest that these particular lifestyles promote horizontal gene transfer, it is difficult to ascertain given the nonrandom sampling of animal genome sequencing projects and the lack of a systematic analysis of animal genomes for such transfers.

Host preference between symbiotic and aposymbiotic Aphis fabae, by the aphid parasitoid, Lysiphlebus ambiguus

Few empirical studies have directly explored the association between Buchnera aphidicola (Enterobacteriales: Enterobacteriaceae), the primary endosymbiont of aphids, and the life history strategies of aphid parasitoids. A series of paired-choice experiments were conducted to explore the preference of the parasitoid Lysiphlebus ambiguus Halliday (Hymenoptera: Aphididae) for symbiotic and aposymbiotic Aphis fabae Scopoli (Hemiptera: Aphididae) and the suitability of these hosts for parasitoid development. When given a choice between symbiotic and aposymbiotic aphids of the same instar, the parasitoid significantly preferred symbiotic over aposymbiotic aphids only during the later instars (L(4) and adult). The suitability of aposymbiotic aphids for parasitoid development was equal to that of symbiotic aphids in terms of survivorship and sex ratio, but was significantly lower than that of symbiotic aphids for L(4) and adult instars in development rate and/or female adult size. When given a choice between similar-sized symbiotic L(2) and aposymbiotic L(4) aphids, the parasitoid preferred the former. No significant differences in preference or host suitability were demonstrated when the parasitoid was given a choice between different instars of aposymbiotic aphids. While parasitoid lifetime fecundity increased with aphid instar at the time of oviposition, there was no significant influence of previous development from symbiotic versus aposymbiotic aphids. These results suggest that while L. ambiguus can discriminate between symbiotic and aposymbiotic A. fabae during later instars and when the aphids are of a similar size, the primary endosymbiont is not needed for successful parasitoid development; and its absence only compromises parasitoid growth reared from later instar aposymbiotic host.

Graph-based analysis of the metabolic exchanges between two co-resident intracellular symbionts, Baumannia cicadellinicola and Sulcia muelleri, with their insect host, Homalodisca coagulata

Endosymbiotic bacteria from different species can live inside cells of the same eukaryotic organism. Metabolic exchanges occur between host and bacteria but also between different endocytobionts. Since a complete genome annotation is available for both, we built the metabolic network of two endosymbiotic bacteria, Sulcia muelleri and Baumannia cicadellinicola, that live inside specific cells of the sharpshooter Homalodisca coagulata and studied the metabolic exchanges involving transfers of carbon atoms between the three. We automatically determined the set of metabolites potentially exogenously acquired (seeds) for both metabolic networks. We show that the number of seeds needed by both bacteria in the carbon metabolism is extremely reduced. Moreover, only three seeds are common to both metabolic networks, indicating that the complementarity of the two metabolisms is not only manifested in the metabolic capabilities of each bacterium, but also by their different use of the same environment. Furthermore, our results show that the carbon metabolism of S. muelleri may be completely independent of the metabolic network of B. cicadellinicola. On the contrary, the carbon metabolism of the latter appears dependent on the metabolism of S. muelleri, at least for two essential amino acids, threonine and lysine. Next, in order to define which subsets of seeds (precursor sets) are sufficient to produce the metabolites involved in a symbiotic function, we used a graph-based method, PITUFO, that we recently developed. Our results highly refine our knowledge about the complementarity between the metabolisms of the two bacteria and their host. We thus indicate seeds that appear obligatory in the synthesis of metabolites are involved in the symbiotic function. Our results suggest both B. cicadellinicola and S. muelleri may be completely independent of the metabolites provided by the co-resident endocytobiont to produce the carbon backbone of the metabolites provided to the symbiotic system (., thr and lys are only exploited by B. cicadellinicola to produce its proteins).

Fitness and freezing: vector biology and human health

Microbes transmitted to mammals by arthropods contend with many factors that could impede survival. To survive, host fitness with infection must outweigh costs. In this issue of the JCI, Neelakanta et al. demonstrate that ticks infected with Anaplasma phagocytophilum show enhanced fitness against freezing injury owing to induced expression of tick "antifreeze glycoprotein." This allows A. phagocytophilum to successfully propagate and survive to cause disease in nonnatural hosts, such as humans. How an intracellular microbe with a small genome subverts host cell function for survival provides insight into the control of some cellular function programs and underscores how vector biology can have an impact on human health.