Key Transport and Ammonia Recycling Genes Involved in Aphid Symbiosis Respond to Host-Plant Specialization

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

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

Comparative transcriptomics of aphid species that diverged > 22 MYA reveals genes that are important for the maintenance of their symbiosis

Most plant-sap feeding insects have obligate relationships with maternally transmitted bacteria. Aphids require their nutritional endosymbiont, Buchnera aphidicola, for the production of essential amino acids. Such endosymbionts are harbored inside of specialized insect cells called bacteriocytes. Here, we use comparative transcriptomics of bacteriocytes between two recently diverged aphid species, Myzus persicae and Acyrthosiphon pisum, to identify key genes that are important for the maintenance of their nutritional mutualism. The majority of genes with conserved expression profiles in M. persicae and A. pisum are for orthologs previously identified in A. pisum to be important for the symbiosis. However, asparaginase which produces aspartate from asparagine was significantly up-regulated only in A. pisum bacteriocytes, potentially because Buchnera of M. persicae encodes its own asparaginase enzyme unlike Buchnera of A. pisum resulting in Buchnera of A. pisum to be dependent on its aphid host for aspartate. One-to-one orthologs that explained the most amount of variation for bacteriocyte specific mRNA expression for both species includes a collaborative gene for methionine biosynthesis, multiple transporters, a horizontally transmitted gene, and secreted proteins. Finally, we highlight species-specific gene clusters which may contribute to host adaptations and/or accommodations in gene regulation to changes in the symbiont or the symbiosis.

Parasitization of Aphis gossypii Glover by Binodoxys communis Gahan Causes Shifts in the Ovarian Bacterial Microbiota

BACKGROUND

Aphis gossypii Glover is an important agricultural pest distributed worldwide. Binodoxys communis Gahan is the main parasitoid wasp of A. gossypii. Previous studies have shown that parasitization causes reduced egg production in A. gossypii, but the effects of parasitism on the symbiotic bacteria in the host ovaries are unknown.

RESULTS

In this study, we analyzed the microbial communities in the ovaries of A. gossypii without and after parasitization. Whether parasitized or not, Buchnera was the dominant genus of symbiotic bacteria in the ovaries, followed by facultative symbionts including Arsenophonus, Pseudomonas, and Acinetobacter. The relative abundance of Buchnera in the aphid ovary increased after parasitization for 1 d in both third-instar nymph and adult stages, but decreased after parasitization for 3 d. The shifts in the relative abundance of Arsenophonus in both stages were the same as those observed for Buchnera. In addition, the relative abundance of Serratia remarkably decreased after parasitization for 1 d and increased after parasitization for 3 d. A functional predictive analysis of the control and parasitized ovary microbiomes revealed that pathways primarily enriched in parasitization were "amino acid transport and metabolism" and "energy production and conversion." Finally, RT-qPCR analysis was performed on Buchnera, Arsenophonus, and Serratia. The results of RT-qPCR were the same as the results of 16S rDNA sequencing.

CONCLUSIONS

These results provide a framework for investigating shifts in the microbial communities in host ovaries, which may be responsible for reduced egg production in aphids. These findings also broaden our understanding of the interactions among aphids, parasitoid wasps, and endosymbionts.

Proton gradient mediates hemolymph trehalose influx into aphid bacteriocytes

Aphids harbor proteobacterial endosymbionts such as Buchnera aphidicola housed in specialized bacteriocytes derived from host cells. The endosymbiont Buchnera supplies essential amino acids such as arginine to the host cells and, in turn, obtains sugars needed for its survival from the hemolymph. The mechanism of sugar supply in aphid bacteriocytes has been rarely studied. It also remains unclear how Buchnera acquires its carbon source. The hemolymph sugars in Acyrthosiphon pisum are composed of the disaccharide trehalose containing two glucose molecules. Here, we report for the first time that trehalose is transported and used as a potential carbon source by Buchnera across the bacteriocyte plasma membrane via trehalose transporters. The current study characterized the bacteriocyte trehalose transporter Ap_ST11 (LOC100159441) using the Xenopus oocyte expression system. The Ap_ST11 transporter was found to be proton-dependent with a K_m value ≥700 mM. We re-examined the hemolymph trehalose at 217.8 mM using a fluorescent trehalose sensor. The bacteriocytes did not obtain trehalose by facilitated diffusion along the gradient across cellular membranes. These findings suggest that trehalose influx into the bacteriocytes depends on the extracellular proton-driven secondary electrochemical transporter.

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.

The genetic adaptations of Toxoptera aurantii facilitated its rapid multiple plant hosts dispersal and invasion

Toxoptera aurantii Boyer de Fonscolombe (Hemiptera: Aphididae) can attack many plant hosts, including tea (Camellia sinensis L.), citrus (Citrus spp.), lychee (Litchi chinensis Sonn.), banana (Musa spp.), and pineapple (Ananas comasus L.) among others. It is a widely distributed hexapod and one of the most destructive pests in tea plantations, causing enormous economic losses in tea production each year. A high-quality reference genome is important to study the phylogenetics and evolution of T. aurantii because its genome is highly heterozygous and repetitive. We obtained a de novo genome assembly of T. aurantii at the chromosome level using a combination of long Nanopore reads from sequencing with high-throughput chromosome conformation capture technology. When finally assembled, the genome was 318.95 Mb on four chromosomes with a 15.19 Mb scaffold N50. A total of 12,162 genes encoded proteins, while there were 22.01% repetitive sequences that totaled 67.73 Mb. Phylogenetic analyses revealed that T. aurantii and Aphis gossypii parted ways approximately 7.6 million years ago (Mya). We used a combination of long-read single-molecule sequencing with Hi-C-based chromatin interaction maps that resulted in a reference chromosomal level reference genome of T. aurantii that was high quality. Our results will enable the exploration of the genetics behind the special biological features of T. aurantii and also provide a source of data that should be useful to compare the compare genome among the Hemiptera.

Symbiosis maintenance in the facultative coral, Oculina arbuscula, relies on nitrogen cycling, cell cycle modulation, and immunity

Symbiosis with unicellular algae in the family Symbiodiniaceae is common across tropical marine invertebrates. Reef-building corals offer a clear example of cellular dysfunction leading to a dysbiosis that disrupts entire ecosystems in a process termed coral bleaching. Due to their obligate symbiotic relationship, understanding the molecular underpinnings that sustain this symbiosis in tropical reef-building corals is challenging, as any aposymbiotic state is inherently coupled with severe physiological stress. Here, we leverage the subtropical, facultatively symbiotic and calcifying coral Oculina arbuscula to investigate gene expression differences between aposymbiotic and symbiotic branches within the same colonies under baseline conditions. We further compare gene ontology (GO) and KOG enrichment in gene expression patterns from O. arbuscula with prior work in the sea anemone Exaiptasia pallida (Aiptasia) and the salamander Ambystoma maculatum-both of which exhibit endophotosymbiosis with unicellular algae. We identify nitrogen cycling, cell cycle control, and immune responses as key pathways involved in the maintenance of symbiosis under baseline conditions. Understanding the mechanisms that sustain a healthy symbiosis between corals and Symbiodiniaceae algae is of urgent importance given the vulnerability of these partnerships to changing environmental conditions and their role in the continued functioning of critical and highly diverse marine ecosystems.

Chromosome-level genome assemblies of two cotton-melon aphid Aphis gossypii biotypes unveil mechanisms of host adaption

The cotton-melon aphid Aphis gossypii is a sap-sucking insect that is considered a serious global pest. The species is distributed over a large geographical range and uses a wide variety of hosts, with some populations being specialized to attack different plant species. Here, we provide de novo chromosome-level genome assemblies of a cotton specialist population (Hap1) and a cucurbit specialist population (Hap3). We achieved this by using a combination of third-generation sequencing platforms, namely Illumina and Hi-C sequencing technologies. We were able to anchor a total of 334.89 Mb (scaffold N50 of 89.13 Mb) and 359.95 Mb (scaffold N50 of 68.88 Mb) to four chromosomes for Hap1 and Hap3, respectively. Moreover, our results showed that the X-chromosome of Hap3 (113.01 Mb) was significantly longer than that of Hap1 (100.26 Mb), with a high level of sequence conservation between the aphid species. We also report variation in the number of protein-coding genes and repeat sequences between Hap1 and Hap3. In particular, olfactory and gustatory receptor genes underwent a high level of gene duplication and expansion events in A. gossypii, including between Hap1 and Hap3. Moreover, we identified two glutathione S-transferase genes which underwent single gene duplications in Hap3, and tandem duplication and inversion events affecting the cytochrome P450 monooxygenase between Hap1 and Hap3, all of which include the CYP3 family. Our results illustrate the variance in the genomic composition of two specialized A. gossypii populations and provide a helpful resource for the study of aphid population evolution, host adaption and insecticide resistance.

The boom and bust of the aphid's essential amino acid metabolism across nymphal development

Within long-term symbioses, animals integrate their physiology and development with their symbiont. In a model nutritional mutualism, aphids harbor the endosymbiont, Buchnera, within specialized bacteriocyte cells. Buchnera synthesizes essential amino acids (EAAs) and vitamins for their host, which are lacking from the aphid's plant sap diet. It is unclear if the aphid host differentially expresses aphid EAA metabolism pathways and genes that collaborate with Buchnera for the production of EAA and vitamins throughout nymphal development when feeding on plants. It is also unclear if aphid bacteriocytes are differentially methylated throughout aphid development as DNA methylation may play a role in gene regulation. By analyzing aphid gene expression, we determined that the bacteriocyte is metabolically more active in metabolizing Buchnera's EAAs and vitamins early in nymphal development compared to intermediate or later immature and adult lifestages. The largest changes in aphid bacteriocyte gene expression, especially for aphid genes that collaborate with Buchnera, occurred during the 3rd to 4th instar transition. During this transition, there is a huge shift in the bacteriocyte from a high energy "nutrient-consuming state" to a "recovery and growth state" where patterning and signaling genes and pathways are upregulated and differentially methylated, and de novo methylation is reduced as evidenced by homogenous DNA methylation profiles after the 2nd instar. Moreover, bacteriocyte number increased and Buchnera's titer decreased throughout aphid nymphal development. These data suggest in combination that bacteriocytes of older nymphal and adult lifestages depend less on the nutritional symbiosis compared to early nymphal lifestages.

An eco-systems biology approach for modeling tritrophic networks reveals the influence of dietary amino acids on symbiont dynamics of Bemisia tabaci

Metabolic conversions allow organisms to produce essential metabolites from the available nutrients in an environment, frequently requiring metabolic exchanges among co-inhabiting organisms. Here, we applied genomic-based simulations for exploring tri-trophic interactions among the sap-feeding insect whitefly (Bemisia tabaci), its host-plants, and symbiotic bacteria. The simplicity of this ecosystem allows capturing the interacting organisms (based on genomic data) and the environmental content (based on metabolomics data). Simulations explored the metabolic capacities of insect-symbiont combinations under environments representing natural phloem. Predictions were correlated with experimental data on the dynamics of symbionts under different diets. Simulation outcomes depict a puzzle of three-layer origins (plant-insect-symbionts) for the source of essential metabolites across habitats and stratify interactions enabling the whitefly to feed on diverse hosts. In parallel to simulations, natural and artificial feeding experiments provide supporting evidence for an environment-based effect on symbiont dynamics. Based on simulations, a decrease in the relative abundance of a symbiont can be associated with a loss of fitness advantage due to an environmental excess in amino-acids whose production in a deprived environment used to depend on the symbiont. The study demonstrates that genomic-based predictions can bridge environment and community dynamics and guide the design of symbiont manipulation strategies.

MetNet: A two-level approach to reconstructing and comparing metabolic networks

Metabolic pathway comparison and interaction between different species can detect important information for drug engineering and medical science. In the literature, proposals for reconstructing and comparing metabolic networks present two main problems: network reconstruction requires usually human intervention to integrate information from different sources and, in metabolic comparison, the size of the networks leads to a challenging computational problem. We propose to automatically reconstruct a metabolic network on the basis of KEGG database information. Our proposal relies on a two-level representation of the huge metabolic network: the first level is graph-based and depicts pathways as nodes and relations between pathways as edges; the second level represents each metabolic pathway in terms of its reactions content. The two-level representation complies with the KEGG database, which decomposes the metabolism of all the different organisms into “reference” pathways in a standardised way. On the basis of this two-level representation, we introduce some similarity measures for both levels. They allow for both a local comparison, pathway by pathway, and a global comparison of the entire metabolism. We developed a tool, MetNet, that implements the proposed methodology. MetNet makes it possible to automatically reconstruct the metabolic network of two organisms selected in KEGG and to compare their two networks both quantitatively and visually. We validate our methodology by presenting some experiments performed with MetNet.

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."

Gene Expression and Diet Breadth in Plant-Feeding Insects: Summarizing Trends

Transcriptomic studies lend insights into the role of transcriptional plasticity in adaptation and specialization. Recently, there has been growing interest in understanding the relationship between variation in herbivorous insect gene expression and the evolution of diet breadth. We review the studies that have emerged on insect gene expression and host plant use, and outline the questions and approaches in the field. Many candidate genes underlying herbivory and specialization have been identified, and a few key studies demonstrate increased transcriptional plasticity associated with generalist compared with specialist species. Addressing the roles that transcriptional variation plays in insect diet breadth will have important implications for our understanding of the evolution of specialization and the genetic and environmental factors that govern insect-plant interactions.

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.

The Effects of Different Diets and Transgenerational Stress on Acyrthosiphon pisum Development

Despite the fact that sap-feeding hemipterans are major agricultural pests, little is known about the pea aphid’s (Acyrthosiphon pisum) nymphal development, compared to other insect models. Given our limited understanding of A. pisum nymphal development and variability in the naming/timing of its developmental events between different environmental conditions and studies, here, we address developmental knowledge gaps by elucidating how diet impacts A. pisum nymphal development for the LSR1 strain when it develops on its universal host plant (Vicia faba), isolated leaves, and artificial diet. Moreover, we test how plant age and transgenerational stressors, such as overcrowding and low plant vigor, can affect nymphal development. We also validate a morphological method to quickly confirm the life stage of each nymphal instar within a mixed population. Overall, we found extremely high variation in the timing of developmental events and a significant delay in nymphal (~5–25-h/instar) and pre-reproductive adult (~40-h) development when reared on isolated leaves and artificial diets, compared to intact host plants. Also, delays in development were observed when reared on older host plants (~9–17-h/event, post 2nd instar) or when previous generations were exposed to overcrowding on host plants (~20-h delay in nymph laying) compared to controls.

Trading amino acids at the aphid-Buchnera symbiotic interface

Plant sap-feeding insects thrive despite feeding exclusively on a diet lacking in essential amino acids. This nutritional deficit is countered through endosymbiotic relationships with microbial symbionts. Nonessential amino acids, vital for microbial symbionts, are utilized by symbiont metabolic pathways and yield essential amino acids required by their eukaryotic hosts. Symbionts are completely dependent on their host to meet nutritional requirements. The endosymbionts are surrounded individually by host-derived symbiosomal membranes and are housed within specialized host bacteriocyte cells. The transport capabilities of the symbiosomal membrane remain unknown. Here, we identify a transport system that mediates a crucial step in this metabolic complementarity: a transporter capable of transporting nonessential amino acids across the symbiosomal membrane of the pea aphid Acyrthosiphon pisum.

Plant sap-feeding insects are widespread, having evolved to occupy diverse environmental niches despite exclusive feeding on an impoverished diet lacking in essential amino acids and vitamins. Success depends exquisitely on their symbiotic relationships with microbial symbionts housed within specialized eukaryotic bacteriocyte cells. Each bacteriocyte is packed with symbionts that are individually surrounded by a host-derived symbiosomal membrane representing the absolute host–symbiont interface. The symbiosomal membrane must be a dynamic and selectively permeable structure to enable bidirectional and differential movement of essential nutrients, metabolites, and biosynthetic intermediates, vital for growth and survival of host and symbiont. However, despite this crucial role, the molecular basis of membrane transport across the symbiosomal membrane remains unresolved in all bacteriocyte-containing insects. A transport protein was immunolocalized to the symbiosomal membrane separating the pea aphid Acyrthosiphon pisum from its intracellular symbiont Buchnera aphidicola. The transporter, A. pisum nonessential amino acid transporter 1, or ApNEAAT1 (gene: ACYPI008971), was characterized functionally following heterologous expression in Xenopus oocytes, and mediates both inward and outward transport of small dipolar amino acids (serine, proline, cysteine, alanine, glycine). Electroneutral ApNEAAT1 transport is driven by amino acid concentration gradients and is not coupled to transmembrane ion gradients. Previous metabolite profiling of hemolymph and bacteriocyte, alongside metabolic pathway analysis in host and symbiont, enable prediction of a physiological role for ApNEAAT1 in bidirectional host–symbiont amino acid transfer, supplying both host and symbiont with indispensable nutrients and biosynthetic precursors to facilitate metabolic complementarity.

The chemoreceptors and odorant binding proteins of the soybean and pea aphids

We examined the genome of the soybean aphid, Aphis glycines, and an updated genome assembly of the pea aphid, Acyrthosiphon pisum, for members of the three major families of chemoreceptors, the Odorant Receptors (ORs), Gustatory Receptors (GRs) and Ionotropic Receptors (IRs), as well as the Odorant Binding Proteins (OBPs). The soybean aphid has 47 ORs, 61 GRs, 19 IRs, and 10 OBPs, compared with 87 ORs, 78 Grs, 19 IRs, and 18 OBPs in the pea aphid, with variable numbers of pseudogenes in the OR and GR families. Phylogenetic analysis reveals that while all of the IRs are simple orthologs between these two species, the OR, GR, and OBP families in the pea aphid have experienced major expansions of particular gene lineages and fewer losses of gene lineages. This imbalance in birth-and-death of chemosensory genes has led to the larger pea aphid gene repertoire, which might be related to the broader host range of pea aphids versus the specialization of soybean aphids on a single summer host plant. Examination of the expression levels of these chemosensory genes in parthenogenetic and sexual females and males of pea aphids revealed multiple genes that are differentially expressed in sexual females or males and might be involved in reproductive biology. Examination of the soybean aphid genes in parthenogenetic females under multiple stressors revealed multiple genes whose expression levels changed with heat or starvation stress, the latter potentially important in finding new food sources.