Amino acid transporter expansions associated with the evolution of obligate endosymbiosis in sap-feeding insects (Hemiptera: sternorrhyncha)

The role of insect gut microbiota in host fitness, detoxification and nutrient supplementation

Insects are incredibly diverse, ubiquitous and have successfully flourished out of the dynamic and often unpredictable nature of evolutionary processes. The resident microbiome has accompanied the physical and biological adaptations that enable their continued survival and proliferation in a wide array of environments. The host insect and microbiome's bidirectional relationship exhibits their capability to influence each other's physiology, behavior and characteristics. Insects are reported to rely directly on the microbial community to break down complex food, adapt to nutrient-deficit environments, protect themselves from natural adversaries and control the expression of social behavior. High-throughput metagenomic approaches have enhanced the potential for determining the abundance, composition, diversity and functional activities of microbial fauna associated with insect hosts, enabling in-depth investigation into insect-microbe interactions. We undertook a review of some of the major advances in the field of metagenomics, focusing on insect-microbe interaction, diversity and composition of resident microbiota, the functional capability of endosymbionts and discussions on different symbiotic relationships. The review aims to be a valuable resource on insect gut symbiotic microbiota by providing a comprehensive understanding of how insect gut symbionts systematically perform a range of functions, viz., insecticide degradation, nutritional support and immune fitness. A thorough understanding of manipulating specific gut symbionts may aid in developing advanced insect-associated research to attain health and design strategies for pest management.

Functional Analysis of Amino Acid Transporter Genes ACYPI000536 and ACYPI004320 in Acyrthosiphon pisum

In recent years, pea aphids have become major pests of alfalfa. Our previous study found that "Gannong 5" is a highly aphid-resistant alfalfa variety and that "Lie Renhe" is a susceptible one. The average field susceptibility index of "Gannong 5" was 31.31, and the average field susceptibility index of "Lie Renhe" was 80.34. The uptake and balance of amino acids in insects are usually dependent on amino acid transporters. RT-qPCR was used to detect the relative expression levels of seven amino acid transporter differential genes in the different instar pea aphids fed on resistant and susceptible alfalfa varieties after 24 h, and two key genes were selected. When pea aphids fed on "Gannong 5", the expression of ACYPI004320 was significantly higher than that in pea aphids fed on "Lie Renhe"; however, the expression of ACYPI000536 was significantly lower than that in pea aphids fed on "Lie Renhe". Afterward, the RNA interference with pea aphid ACYPI000536 and ACYPI004320 genes was performed using a plant-mediated method, and gene function was verified via liquid chromatography-mass spectrometry and pea aphid sensitivity to aphid-resistant and susceptible alfalfa varieties. The results showed that the down-regulation of the ACYPI000536 gene expression led to an increase in the histidine and lysine contents in pea aphids, which, in turn, led to an increase in mortality when pea aphids fed on the susceptible variety "Lie Renhe". The down-regulation of the ACYPI004320 gene expression led to an increase in phenylalanine content in pea aphids, which, in turn, led to a decrease in mortality when pea aphids fed on the resistant variety "Gannong 5".

Co-option of a conserved host glutamine transporter facilitates aphid/Buchnera metabolic integration

Organisms across the tree of life colonize novel environments by partnering with bacterial symbionts. These symbioses are characterized by intimate integration of host/endosymbiont biology at multiple levels, including metabolically. Metabolic integration is particularly important for sap-feeding insects and their symbionts, which supplement nutritionally unbalanced host diets. Many studies reveal parallel evolution of host/endosymbiont metabolic complementarity in amino acid biosynthesis, raising questions about how amino acid metabolism is regulated, how regulatory mechanisms evolve, and the extent to which similar mechanisms evolve in different systems. In the aphid/Buchnera symbiosis, the transporter ApGLNT1 (Acyrthosiphon pisum glutamine transporter 1) supplies glutamine, an amino donor in transamination reactions, to bacteriocytes (where Buchnera reside) and is competitively inhibited by Buchnera-supplied arginine-consistent with a role regulating amino acid metabolism given host demand for Buchnera-produced amino acids. We examined how ApGLNT1 evolved a regulatory role by functionally characterizing orthologs in insects with and without endosymbionts. ApGLNT1 orthologs are functionally similar, and orthology searches coupled with homology modeling revealed that GLNT1 is ancient and structurally conserved across insects. Our results indicate that the ApGLNT1 symbiotic regulatory role is derived from its ancestral role and, in aphids, is likely facilitated by loss of arginine biosynthesis through the urea cycle. Given consistent loss of host arginine biosynthesis and retention of endosymbiont arginine supply, we hypothesize that GLNT1 is a general mechanism regulating amino acid metabolism in sap-feeding insects. This work fills a gap, highlighting the broad importance of co-option of ancestral proteins to novel contexts in the evolution of host/symbiont systems.

Genes possibly related to symbiosis in early life stages of Acropora tenuis inoculated with Symbiodinium microadriaticum

Due to the ecological importance of mutualism between reef-building corals and symbiotic algae (Family Symbiodiniaceae), various transcriptomic studies on coral-algal symbiosis have been performed; however, molecular mechanisms, especially genes essential to initiate and maintain these symbioses remain unknown. We investigated transcriptomic responses of Acropora tenuis to inoculation with the native algal symbiont, Symbiodinium microadriaticum, during early life stages, and identified possible symbiosis-related genes. Genes involved in immune regulation, protection against oxidative stress, and metabolic interactions between partners are particularly important for symbiosis during Acropora early life stages. In addition, molecular phylogenetic analysis revealed that some possible symbiosis-related genes originated by gene duplication in the Acropora lineage, suggesting that gene duplication may have been the driving force to establish stable mutualism in Acropora, and that symbiotic molecular mechanisms may vary among coral lineages.

Genome-Wide Analysis of the Amino Acid Auxin Permease (AAAP) Gene Family and Identification of an AAAP Gene Associated with the Growth and Reproduction of the Brown Planthopper, Nilaparvata lugens (Stål)

Amino acids play a vital role in several biological processes in organisms and are mainly acquired through diet by most insects. The amino acid auxin permease (AAAP) transporter family is an important amino acid transporter gene family in insects for the transportation of amino acids into and out of cells across the plasma membrane. Here, we identified 21 putative AAAP family members in the genome of the brown planthopper (BPH), Nilaparvata lugens, a devastating pest that feeds only on the phloem sap of rice plants. Molecular characteristic analysis indicated large variations in protein features and amino acid sequences among the predicted AAAP family members in BPH. Phylogenetic analysis clustered these AAAP transporters into three subgroups, with the members in the same group sharing a similar pattern of conserved motif distribution. Through ortholog gene recognition and spatiotemporal gene expression analysis, the AAAP gene NlAAAP07, which was predicted to regulate BPH larval growth and female fecundity, was identified. RNA interference (RNAi)-mediated suppression of NlAAAP07 significantly postponed the duration of 3rd instar nymphs developing into adults from 7.4 days to 9.0 days, and decreased the oviposition amount and egg hatching rate of females by 30.7% and 11.0%, respectively. Our results provide a foundation for further functional analysis of AAAP transporters in BPH.

Genome-Wide Identification and Characterization of Amino Acid Polyamine Organocation Transporter Family Genes Reveal Their Role in Fecundity Regulation in a Brown Planthopper Species (Nilaparvata lugens)

The brown planthopper (BPH), Nilaparvata lugens Stål (Hemiptera:Delphacidae), is one of the most destructive pests of rice worldwide. As a sap-feeding insect, the BPH is incapable of synthesizing several amino acids which are essential for normal growth and development. Therefore, the insects have to acquire these amino acids from dietary sources or their endosymbionts, in which amino acid transporters (AATs) play a crucial role by enabling the movement of amino acids into and out of insect cells. In this study, a common amino acid transporter gene family of amino acid/polyamine/organocation (APC) was identified in BPHs and analyzed. Based on a homology search and conserved functional domain recognition, 20 putative APC transporters were identified in the BPH genome. Molecular trait analysis showed that the verified BPH APC family members were highly variable in protein features, conserved motif distribution patterns, and exon/intron organization. Phylogenetic analysis of five hemipteran species revealed an evolutionary pattern of interfamily conservation and lineage-specific expansion of this gene family. Moreover, stage- and tissue-specific expression analysis revealed diverse expression patterns in the 20 BPH APC transporter genes. Lastly, a potential BPH fecundity regulatory gene of NlAPC09 was identified and shown to participate in the fecundity regulation through the use of quantitative polymerase chain reaction (qPCR) and RNA inference experiments. Our results provide a basis for further functional investigations of APC transporters in BPH.

Genome and transcriptome analysis of the mealybug Maconellicoccus hirsutus: Correlation with its unique phenotypes

Mealybugs are aggressive pests with world-wide distribution and are suitable for the study of different phenomena like genomic imprinting and epigenetics. Genomic approaches facilitate these studies in absence of robust genetics in this system. We sequenced, de novo assembled, annotated Maconellicoccus hirsutus genome. We carried out comparative genomics it with four mealybug and eight other insect species, to identify expanded, specific and contracted gene classes that relate to pesticide and desiccation resistance. We identified horizontally transferred genes adding to the mutualism between the mealybug and its endosymbionts. Male and female transcriptome analysis indicates differential expression of metabolic pathway genes correlating with their physiology and the genes for sexual dimorphism. The significantly lower expression of endosymbiont genes in males relates to the depletion of endosymbionts in males during development.

Comparative genomics of Klebsiella michiganensis BD177 and related members of Klebsiella sp. reveal the symbiotic relationship with Bactrocera dorsalis

Background

Bactrocera dorsalis is a destructive polyphagous and highly invasive insect pest of tropical and subtropical species of fruit and vegetable crops. The sterile insect technique (SIT) has been used for decades to control insect pests of agricultural, veterinary, and human health importance. Irradiation of pupae in SIT can reduce the ecological fitness of the sterile insects. Our previous study has shown that a gut bacterial strain BD177 that could restore ecological fitness by promoting host food intake and metabolic activities.

Results

Using long-read sequence technologies, we assembled the complete genome of K. michiganensis BD177 strain. The complete genome of K. michiganensis BD177 comprises one circular chromosome and four plasmids with a GC content of 55.03%. The pan-genome analysis was performed on 119 genomes (strain BD177 genome and 118 out of 128 published Klebsiella sp. genomes since ten were discarded). The pan-genome includes a total of 49305 gene clusters, a small number of 858 core genes, and a high number of accessory (10566) genes. Pan-genome and average nucleotide identity (ANI) analysis showed that BD177 is more similar to the type strain K. michiganensis DSM2544, while away from the type strain K. oxytoca ATCC13182. Comparative genome analysis with 21 K. oxytoca and 12 K. michiganensis strains, identified 213 unique genes, several of them related to amino acid metabolism, metabolism of cofactors and vitamins, and xenobiotics biodegradation and metabolism in BD177 genome.

Conclusions

Phylogenomics analysis reclassified strain BD177 as a member of the species K. michiganensis. Comparative genome analysis suggested that K. michiganensis BD177 has the strain-specific ability to provide three essential amino acids (phenylalanine, tryptophan and methionine) and two vitamins B (folate and riboflavin) to B. dorsalis. The clear classification status of BD177 strain and identification of unique genetic characteristics may contribute to expanding our understanding of the symbiotic relationship of gut microbiota and B. dorsalis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-020-00945-0.

Whole-Genome Transcriptome Analyses of Native Symbionts Reveal Host Coral Genomic Novelties for Establishing Coral-Algae Symbioses

Reef-building corals and photosynthetic, endosymbiotic algae of the family Symbiodiniaceae establish mutualistic relationships that are fundamental to coral biology, enabling coral reefs to support a vast diversity of marine species. Although numerous types of Symbiodiniaceae occur in coral reef environments, Acropora corals select specific types in early life stages. In order to study molecular mechanisms of coral–algal symbioses occurring in nature, we performed whole-genome transcriptomic analyses of Acropora tenuis larvae inoculated with Symbiodinium microadriaticum strains isolated from an Acropora recruit. In order to identify genes specifically involved in symbioses with native symbionts in early life stages, we also investigated transcriptomic responses of Acropora larvae exposed to closely related, nonsymbiotic, and occasionally symbiotic Symbiodinium strains. We found that the number of differentially expressed genes was largest when larvae acquired native symbionts. Repertoires of differentially expressed genes indicated that corals reduced amino acid, sugar, and lipid metabolism, such that metabolic enzymes performing these functions were derived primarily from S. microadriaticum rather than from A. tenuis. Upregulated gene expression of transporters for those metabolites occurred only when coral larvae acquired their natural symbionts, suggesting active utilization of native symbionts by host corals. We also discovered that in Acropora, genes for sugar and amino acid transporters, prosaposin-like, and Notch ligand-like, were upregulated only in response to native symbionts, and included tandemly duplicated genes. Gene duplications in coral genomes may have been essential to establish genomic novelties for coral–algae symbiosis.

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

Background

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

Results

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

Conclusions

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

The Identification and Evolutionary Trends of the Solute Carrier Superfamily in Arthropods

The solute carrier (SLC) transporter superfamily comprises an ancient and ubiquitous group of proteins capable of translocating a range of nutrients, endogenous molecules, and xenobiotics. Although the group has been the subject of intense investigation in both bacteria and mammals, its systematic identification in arthropods has not yet been undertaken. Here, we present a genome-wide identification of all 66 human SLC families in 174 arthropod species. A pipeline (SLC_id) was constructed to identify and group SLCs using a combination of hidden Markov model and BLAST searches followed by filtering based on polypeptide length and the number of transmembrane domains. Comparative analysis of the number of transporters in each family across diverse arthropod lineages was accomplished using one-way analysis of variance (ANOVA) and the Computational Analysis of gene Family Evolution (CAFE). These results suggested that many SLC families have undergone expansions or contractions in particular evolutionary lineages. Notably, the sugar transporting SLC2 family was significantly larger in insects compared with arachnids. This difference may have been complemented by a rapid expansion of the SLC60 family in arachnids which also acts on dietary sugars. Furthermore, the SLC33 family underwent a recent and drastic expansion in aphids, although the biological relevance of this expansion was not possible to infer. Information on specific SLC transporter families across arthropod species can be accessed through an R shiny web application at http://chrysalida.imbb.forth.gr : 3838/Arthropod_SLC_Database/. The present study greatly facilitates further investigation of the diverse group of SLC transporters in arthropods.

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

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

Gene Family Evolution in the Pea Aphid Based on Chromosome-Level Genome Assembly

Genome structural variations, including duplications, deletions, insertions, and inversions, are central in the evolution of eukaryotic genomes. However, structural variations present challenges for high-quality genome assembly, hampering efforts to understand the evolution of gene families and genome architecture. An example is the genome of the pea aphid (Acyrthosiphon pisum) for which the current assembly is composed of thousands of short scaffolds, many of which are known to be misassembled. Here, we present an improved version of the A. pisum genome based on the use of two long-range proximity ligation methods. The new assembly contains four long scaffolds (40-170 Mb), corresponding to the three autosomes and the X chromosome of A. pisum, and encompassing 86% of the new assembly. Assembly accuracy is supported by several quality assessments. Using this assembly, we identify the chromosomal locations and relative ages of duplication events, and the locations of horizontally acquired genes. The improved assembly illuminates the mode of gene family evolution by providing proximity information between paralogs. By estimating nucleotide polymorphism and coverage depth from resequencing data, we determined that many short scaffolds not assembling to chromosomes represent hemizygous regions, which are especially frequent on the highly repetitive X chromosome. Aligning the X-linked aphicarus region, responsible for male wing dimorphism, to the new assembly revealed a 50-kb deletion that cosegregates with the winged male phenotype in some clones. These results show that long-range scaffolding methods can substantially improve assemblies of repetitive genomes and facilitate study of gene family evolution and structural variation.

Protocols for the Molecular Evolutionary Analysis of Membrane Protein Gene Duplicates

Gene duplication is an important process in the evolution of gene content in eukaryotic genomes. Understanding when gene duplicates contribute new molecular functions to genomes through molecular adaptation is one important goal in comparative genomics. In large gene families, however, characterizing adaptation and neofunctionalization across species is challenging, as models have traditionally quantified the timing of duplications without considering underlying gene trees. This protocol combines multiple approaches to detect adaptation in protein duplicates at a phylogenetic scale. We include a description of models for gene tree-species tree reconciliation that enable different types of inference, as well as a practical guide to their use. Although simulation-based approaches successfully detect shifts in the rate of duplication/retention, the conflation between the duplication and retention processes, the distinct trajectories of duplicates under non-, sub-, and neofunctionalization, as well as dosage effects offer hitherto unexplored analytical avenues. We introduce mathematical descriptions of these probabilities and offer a road map to computational implementation whose starting point is parsimony reconciliation. Sequence evolution information based on the ratio of nonsynonymous to synonymous nucleotide substitution rates (dN/dS) can be combined with duplicate survival probabilities to better predict the emergence of new molecular functions in retained duplicates. Together, these methods enable characterization of potentially adaptive candidate duplicates whose neofunctionalization may contribute to phenotypic divergence across species.

Divergent evolution peaks under intermediate population bottlenecks during bacterial experimental evolution

There is growing evidence that parallel molecular evolution is common, but its causes remain poorly understood. Demographic parameters such as population bottlenecks are predicted to be major determinants of parallelism. Here, we test the hypothesis that bottleneck intensity shapes parallel evolution by elucidating the genomic basis of adaptation to antibiotic-supplemented media in hundreds of populations of the bacterium Pseudomonas fluorescens Pf0-1. As expected, bottlenecking decreased the rate of phenotypic and molecular adaptation. Surprisingly, bottlenecking had no impact on the likelihood of parallel adaptive molecular evolution at a genome-wide scale. However, bottlenecking had a profound impact on the genes involved in antibiotic resistance. Specifically, under either intense or weak bottlenecking, resistance predominantly evolved by strongly beneficial mutations which provide high levels of antibiotic resistance. In contrast with intermediate bottlenecking regimes, resistance evolved by a greater diversity of genetic mechanisms, significantly reducing the observed levels of parallel genetic evolution. Our results demonstrate that population bottlenecking can be a major predictor of parallel evolution, but precisely how may be more complex than many simple theoretical predictions.

Genome-wide Identification and Expression Analysis of Amino Acid Transporters in the Whitefly, Bemisia tabaci (Gennadius)

The whitefly (Bemisia tabaci) is a cosmopolitan and devastating pest of agricultural crops and ornamentals. B. tabaci causes extensive damage by feeding on phloem and by transmitting plant viruses. Like many other organisms, insects depend on amino acid transporters (AATs) to transport amino acids into and out of its cells. We present a genome-wide and transcriptome-wide investigation of the following two families of AATs in B. tabaci biotype B: amino acid/auxin permease (AAAP) and amino acid/polyamine/organocation (APC). A total of 14 putative APCs and 25 putative AAAPs were identified, and a 10-paralog B. tabaci-specific expansion of AAAPs was found by maximum likelihood phylogeny. Detailed gene structure information revealed that 9 members of the B. tabaci-specific AAAP family expansion closely situated on a same scaffold. Expression profiling of the B. tabaci B APC and AAAP genes as affected by stage and plant host showed diverse expression patterns. The analysis of evolutionary rates indicated that purifying selection can explain the B. tabaci-specific AAAP expansion. RNA interference (RNAi)-mediated suppression of two AAAP genes (BtAAAP15 and BtAAAP21) significantly increased the mortality of B. tabaci B adults. The results provide a foundation for future functional analysis of APC and AAAP genes in B. tabaci.

Annotation of the Asian Citrus Psyllid Genome Reveals a Reduced Innate Immune System

Citrus production worldwide is currently facing significant losses due to citrus greening disease, also known as Huanglongbing. The citrus greening bacteria, Candidatus Liberibacter asiaticus (CLas), is a persistent propagative pathogen transmitted by the Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Liviidae). Hemipterans characterized to date lack a number of insect immune genes, including those associated with the Imd pathway targeting Gram-negative bacteria. The D. citri draft genome was used to characterize the immune defense genes present in D. citri. Predicted mRNAs identified by screening the published D. citri annotated draft genome were manually searched using a custom database of immune genes from previously annotated insect genomes. Toll and JAK/STAT pathways, general defense genes Dual oxidase, Nitric oxide synthase, prophenoloxidase, and cellular immune defense genes were present in D. citri. In contrast, D. citri lacked genes for the Imd pathway, most antimicrobial peptides, 1,3-β-glucan recognition proteins (GNBPs), and complete peptidoglycan recognition proteins. These data suggest that D. citri has a reduced immune capability similar to that observed in A. pisum, P. humanus, and R. prolixus. The absence of immune system genes from the D. citri genome may facilitate CLas infections, and is possibly compensated for by their relationship with their microbial endosymbionts.

Amino acid transporters implicated in endocytosis of Buchnera during symbiont transmission in the pea aphid

Background

Many insects host their obligate, maternally transmitted symbiotic bacteria in specialized cells called bacteriocytes. One of the best-studied insect nutritional endosymbioses is that of the aphid and its endosymbiont, Buchnera aphidicola. Aphids and Buchnera are metabolically and developmentally integrated, but the molecular mechanisms underlying Buchnera transmission and coordination with aphid development remain largely unknown. Previous work using electron microscopy to study aphid asexual embryogenesis has revealed that Buchnera transmission involves exocytosis from a maternal bacteriocyte followed by endocytotic uptake by a blastula. While the importance of exo- and endocytic cellular processes for symbiont transmission is clear, the molecular mechanisms that regulate these processes are not known. Here, we shed light on the molecular mechanisms that regulate Buchnera transmission and developmental integration.

Results

We present the developmental atlas of ACYPI000536 and ACYPI008904 mRNAs during asexual embryogenesis in the pea aphid, Acyrthosiphon pisum. Immediately before Buchnera invasion, transcripts of both genes were detected by whole-mount in situ hybridization in the posterior syncytial nuclei of late blastula embryos. Following Buchnera invasion, expression of both genes was identified in the region occupied by Buchnera throughout embryogenesis. Notably during Buchnera migration, expression of both genes was not concomitant with the entirety of the bacterial mass but rather expression colocalized with Buchnera in the anterior region of the bacterial mass. In addition, we found that ACYPI000536 was expressed in nuclei at the leading edge of the bacterial mass, joining the bacterial mass in subsequent developmental stages. Finally, quantitative reverse transcription real-time PCR suggested that early in development both transcripts were maternally provisioned to embryos.

Conclusions

We venture that ACYPI000536 and ACYPI008904 function as nutrient sensors at the site of symbiont invasion to facilitate TOR-pathway-mediated endocytosis of Buchnera by the aphid blastula. Our data support earlier reports of bacteriocyte determination involving a two-step recruitment process but suggest that the second wave of recruitment occurs earlier than previously described. Finally, our work highlights that bacteriocyte-enriched amino acid transporter paralogs have additionally been retained to play novel developmental roles in both symbiont recruitment and bacteriome development.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-016-0061-7) contains supplementary material, which is available to authorized users.

Gene Family Expansions in Aphids Maintained by Endosymbiotic and Nonsymbiotic Traits

Facilitating the evolution of new gene functions, gene duplication is a major mechanism driving evolutionary innovation. Gene family expansions relevant to host/symbiont interactions are increasingly being discovered in eukaryotes that host endosymbiotic microbes. Such discoveries entice speculation that gene duplication facilitates the evolution of novel, endosymbiotic relationships. Here, using a comparative transcriptomic approach combined with differential gene expression analysis, we investigate the importance of endosymbiosis in retention of amino acid transporter paralogs in aphid genomes. To pinpoint the timing of amino acid transporter duplications we inferred gene phylogenies for five aphid species and three outgroups. We found that while some duplications arose in the aphid common ancestor concurrent with endosymbiont acquisition, others predate aphid divergence from related insects without intracellular symbionts, and still others appeared during aphid diversification. Interestingly, several aphid-specific paralogs have conserved enriched expression in bacteriocytes, the insect cells that host primary symbionts. Conserved bacteriocyte enrichment suggests that the transporters were recruited to the aphid/endosymbiont interface in the aphid common ancestor, consistent with a role for gene duplication in facilitating the evolution of endosymbiosis in aphids. In contrast, the temporal variability of amino acid transporter duplication indicates that endosymbiosis is not the only trait driving selection for retention of amino acid transporter paralogs in sap-feeding insects. This study cautions against simplistic interpretations of the role of gene family expansion in the evolution of novel host/symbiont interactions by further highlighting that multiple complex factors maintain gene family paralogs in the genomes of eukaryotes that host endosymbiotic microbes.

Ontogenetic differences in localization of glutamine transporter ApGLNT1 in the pea aphid demonstrate that mechanisms of host/symbiont integration are not similar in the maternal versus embryonic bacteriome

Background

Obligate intracellular symbionts of insects are metabolically and developmentally integrated with their hosts. Typically, reproduction fails in many insect nutritional endosymbioses when host insects are cured of their bacterial symbionts, and yet remarkably little is known about the processes that developmentally integrate host and symbiont. Here in the best studied insect obligate intracellular symbiosis, that of the pea aphid, Acyrthosiphon pisum, with the gammaproteobacterium Buchnera aphidicola, we tracked the expression and localization of amino acid transporter ApGLNT1 gene products during asexual embryogenesis. Recently being characterized as a glutamine transporter, ApGLNT1 has been proposed to be a key regulator of amino acid biosynthesis in A. pisum bacteriocytes. To determine when this important mediator of the symbiosis becomes expressed in aphid embryonic bacteriocytes, we applied whole-mount in situ hybridization and fluorescent immunostaining with a specific anti-ApGLNT1 antibody to detect the temporal and spatial expression of ApGLNT1 gene products during asexual embryogenesis.

Results

During embryogenesis, ApGLNT1 mRNA and protein localize to the follicular epithelium that surrounds parthenogenetic viviparous embryos, where we speculate that it functions to supply developing embryos with glutamine from maternal hemolymph. Unexpectedly, in the embryonic bacteriome ApGLNT1 protein does not localize to the membrane of bacteriocytes, a pattern that leads us to conclude that the regulation of amino acid metabolism in the embryonic bacteriome mechanistically differs from that in the maternal bacteriome. Paralleling our earlier report of punctate cytoplasmic localization of ApGLNT1 in maternal bacteriocytes, we find ApGLNT1 protein localizing as cytoplasmic puncta throughout development in association with Buchnera.

Conclusions

Our work that documents ontogenetic shifts in the localization of ApGLNT1 protein in the host bacteriome demonstrates that maternal and embryonic bacteriomes are not equivalent. Significantly, the persistent punctate cytoplasmic localization of ApGLNT1 in association with Buchnera in embryos prior to bacteriocyte formation and later in both embryonic and maternal bacteriomes suggests that ApGLNT1 plays multiple roles in this symbiosis, roles that include amino acid transport and possibly nutrient sensing.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0038-y) contains supplementary material, which is available to authorized users.

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.