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.

SymbiQuant: A Machine Learning Object Detection Tool for Polyploid Independent Estimates of Endosymbiont Population Size

Quantifying the size of endosymbiont populations is challenging because endosymbionts are typically difficult or impossible to culture and commonly polyploid. Current approaches to estimating endosymbiont population sizes include quantitative PCR (qPCR) targeting endosymbiont genomic DNA and flow-cytometry. While qPCR captures genome copy number data, it does not capture the number of bacterial cells in polyploid endosymbiont populations. In contrast, flow cytometry can capture accurate estimates of whole host-level endosymbiont population size, but it is not readily able to capture data at the level of endosymbiotic host cells. To complement these existing approaches for estimating endosymbiont population size, we designed and implemented an object detection/segmentation tool for counting the number of endosymbiont cells in micrographs of host tissues. The tool, called SymbiQuant, which makes use of recent advances in deep neural networks includes a graphic user interface that allows for human curation of tool output. We trained SymbiQuant for use in the model aphid/Buchnera endosymbiosis and studied Buchnera population dynamics and phenotype over aphid postembryonic development. We show that SymbiQuant returns accurate counts of endosymbionts, and readily captures Buchnera phenotype. By replacing our training data with data composed of annotated microscopy images from other models of endosymbiosis, SymbiQuant has the potential for broad application. Our tool, which is available on GitHub, adds to the repertoire of methods researchers can use to study endosymbiosis at the organismal, genome, and now endosymbiotic host tissue or cell levels.

Global patterns in genomic diversity underpinning the evolution of insecticide resistance in the aphid crop pest Myzus persicae

The aphid Myzus persicae is a destructive agricultural pest that displays an exceptional ability to develop resistance to both natural and synthetic insecticides. To investigate the evolution of resistance in this species we generated a chromosome-scale genome assembly and living panel of >110 fully sequenced globally sampled clonal lines. Our analyses reveal a remarkable diversity of resistance mutations segregating in global populations of M. persicae. We show that the emergence and spread of these mechanisms is influenced by host-plant associations, uncovering the widespread co-option of a host-plant adaptation that also offers resistance against synthetic insecticides. We identify both the repeated evolution of independent resistance mutations at the same locus, and multiple instances of the evolution of novel resistance mechanisms against key insecticides. Our findings provide fundamental insights into the genomic responses of global insect populations to strong selective forces, and hold practical relevance for the control of pests and parasites.

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 genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest

BACKGROUND

Although native to North America, the invasion of the aphid-like grape phylloxera Daktulosphaira vitifoliae across the globe altered the course of grape cultivation. For the past 150 years, viticulture relied on grafting-resistant North American Vitis species as rootstocks, thereby limiting genetic stocks tolerant to other stressors such as pathogens and climate change. Limited understanding of the insect genetics resulted in successive outbreaks across the globe when rootstocks failed. Here we report the 294-Mb genome of D. vitifoliae as a basic tool to understand host plant manipulation, nutritional endosymbiosis, and enhance global viticulture.

RESULTS

Using a combination of genome, RNA, and population resequencing, we found grape phylloxera showed high duplication rates since its common ancestor with aphids, but similarity in most metabolic genes, despite lacking obligate nutritional symbioses and feeding from parenchyma. Similarly, no enrichment occurred in development genes in relation to viviparity. However, phylloxera evolved > 2700 unique genes that resemble putative effectors and are active during feeding. Population sequencing revealed the global invasion began from the upper Mississippi River in North America, spread to Europe and from there to the rest of the world.

CONCLUSIONS

The grape phylloxera genome reveals genetic architecture relative to the evolution of nutritional endosymbiosis, viviparity, and herbivory. The extraordinary expansion in effector genes also suggests novel adaptations to plant feeding and how insects induce complex plant phenotypes, for instance galls. Finally, our understanding of the origin of this invasive species and its genome provide genetics resources to alleviate rootstock bottlenecks restricting the advancement of viticulture.

microRNA regulation in an ancient obligate endosymbiosis

Although many insects are associated with obligate bacterial endosymbionts, the mechanisms by which these host/endosymbiont associations are regulated remain mysterious. While microRNAs (miRNAs) have been recently identified as regulators of host/microbe interactions, including host/pathogen and host/facultative endosymbiont interactions, the role miRNAs may play in mediating host/obligate endosymbiont interactions is virtually unknown. Here, we identified conserved miRNAs that potentially mediate symbiotic interactions between aphids and their obligate endosymbiont, Buchnera aphidicola. Using small RNA sequence data from Myzus persicae and Acyrthosiphon pisum, we annotated 93 M. persicae and 89 A. pisum miRNAs, among which 69 were shared. We found 14 miRNAs that were either highly expressed in aphid bacteriome, the Buchnera-housing tissue, or differentially expressed in bacteriome vs. gut, a non-Buchnera-housing tissue. Strikingly, 10 of these 14 miRNAs have been implicated previously in other host/microbe interaction studies. Investigating the interaction networks of these miRNAs using a custom computational pipeline, we identified 103 miRNA::mRNA interactions shared between M. persicae and A. pisum. Functional annotation of the shared mRNA targets revealed only two over-represented cluster of orthologous group categories: amino acid transport and metabolism, and signal transduction mechanisms. Our work supports a role for miRNAs in mediating host/symbiont interactions between aphids and their obligate endosymbiont Buchnera. In addition, our results highlight the probable importance of signal transduction mechanisms to host/endosymbiont coevolution.

Rapid transcriptional plasticity of duplicated gene clusters enables a clonally reproducing aphid to colonise diverse plant species

Background

The prevailing paradigm of host-parasite evolution is that arms races lead to increasing specialisation via genetic adaptation. Insect herbivores are no exception and the majority have evolved to colonise a small number of closely related host species. Remarkably, the green peach aphid, Myzus persicae, colonises plant species across 40 families and single M. persicae clonal lineages can colonise distantly related plants. This remarkable ability makes M. persicae a highly destructive pest of many important crop species.

Results

To investigate the exceptional phenotypic plasticity of M. persicae, we sequenced the M. persicae genome and assessed how one clonal lineage responds to host plant species of different families. We show that genetically identical individuals are able to colonise distantly related host species through the differential regulation of genes belonging to aphid-expanded gene families. Multigene clusters collectively upregulate in single aphids within two days upon host switch. Furthermore, we demonstrate the functional significance of this rapid transcriptional change using RNA interference (RNAi)-mediated knock-down of genes belonging to the cathepsin B gene family. Knock-down of cathepsin B genes reduced aphid fitness, but only on the host that induced upregulation of these genes.

Conclusions

Previous research has focused on the role of genetic adaptation of parasites to their hosts. Here we show that the generalist aphid pest M. persicae is able to colonise diverse host plant species in the absence of genetic specialisation. This is achieved through rapid transcriptional plasticity of genes that have duplicated during aphid evolution.

Electronic supplementary material

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

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.

Proton-dependent glutamine uptake by aphid bacteriocyte amino acid transporter ApGLNT1

Aphids house large populations of the gammaproteobacterial symbiont Buchnera aphidicola in specialized bacteriocyte cells. The combined biosynthetic capability of the holobiont (Acyrthosiphon pisum and Buchnera) is sufficient for biosynthesis of all twenty protein coding amino acids, including amino acids that animals alone cannot synthesize; and that are present at low concentrations in A. pisum's plant phloem sap diet. Collaborative holobiont amino acid biosynthesis depends on glutamine import into bacteriocytes, which serves as a nitrogen-rich amino donor for biosynthesis of other amino acids. Recently, we characterized A. pisum glutamine transporter 1 (ApGLNT1), a member of the amino acid/auxin permease family, as the dominant bacteriocyte plasma membrane glutamine transporter. Here we show ApGLNT1 to be structurally and functionally related to mammalian proton-dependent amino acid transporters (PATs 1-4). Using functional expression in Xenopus laevis oocytes, combined with two-electrode voltage clamp electrophysiology we demonstrate that ApGLNT1 is electrogenic and that glutamine induces large inward currents. ApGLNT1 glutamine induced currents are dependent on external glutamine concentration, proton (H+) gradient across the membrane, and membrane potential. Based on these transport properties, ApGLNT1-mediated glutamine uptake into A. pisum bacteriocytes can be regulated by changes in either proton gradients across the plasma membrane or membrane potential.

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

Background

Mutualistic obligate endosymbioses shape the evolution of endosymbiont genomes, but their impact on host genomes remains unclear. Insects of the sub-order Sternorrhyncha (Hemiptera) depend on bacterial endosymbionts for essential amino acids present at low abundances in their phloem-based diet. This obligate dependency has been proposed to explain why multiple amino acid transporter genes are maintained in the genomes of the insect hosts. We implemented phylogenetic comparative methods to test whether amino acid transporters have proliferated in sternorrhynchan genomes at rates grater than expected by chance.

Results

By applying a series of methods to reconcile gene and species trees, inferring the size of gene families in ancestral lineages, and simulating the null process of birth and death in multi-gene families, we uncovered a 10-fold increase in duplication rate in the AAAP family of amino acid transporters within Sternorrhyncha. This gene family expansion was unmatched in other closely related clades lacking endosymbionts that provide essential amino acids.

Conclusions

Our findings support the influence of obligate endosymbioses on host genome evolution by both inferring significant expansions of gene families involved in symbiotic interactions, and discovering increases in the rate of duplication associated with multiple emergences of obligate symbiosis in Sternorrhyncha.

Electronic supplementary material

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

A substrate ambiguous enzyme facilitates genome reduction in an intracellular symbiont

BACKGROUND

Genome evolution in intracellular microbial symbionts is characterized by gene loss, generating some of the smallest and most gene-poor genomes known. As a result of gene loss these genomes commonly contain metabolic pathways that are fragmented relative to their free-living relatives. The evolutionary retention of fragmented metabolic pathways in the gene-poor genomes of endosymbionts suggests that they are functional. However, it is not always clear how they maintain functionality. To date, the fragmented metabolic pathways of endosymbionts have been shown to maintain functionality through complementation by host genes, complementation by genes of another endosymbiont and complementation by genes in host genomes that have been horizontally acquired from a microbial source that is not the endosymbiont. Here, we demonstrate a fourth mechanism.

RESULTS

We investigate the evolutionary retention of a fragmented pathway for the essential nutrient pantothenate (vitamin B5) in the pea aphid, Acyrthosiphon pisum endosymbiosis with Buchnera aphidicola. Using quantitative analysis of gene expression we present evidence for complementation of the Buchnera pantothenate biosynthesis pathway by host genes. Further, using complementation assays in an Escherichia coli mutant we demonstrate functional replacement of a pantothenate biosynthesis enzyme, 2-dehydropantoate 2-reductase (E.C. 1.1.1.169), by an endosymbiont gene, ilvC, encoding a substrate ambiguous enzyme.

CONCLUSIONS

Earlier studies have speculated that missing enzyme steps in fragmented endosymbiont metabolic pathways are completed by adaptable endosymbiont enzymes from other pathways. Here, we experimentally demonstrate completion of a fragmented endosymbiont vitamin biosynthesis pathway by recruitment of a substrate ambiguous enzyme from another pathway. In addition, this work extends host/symbiont metabolic collaboration in the aphid/Buchnera symbiosis from amino acid metabolism to include vitamin biosynthesis.

Dynamic recruitment of amino acid transporters to the insect/symbiont interface

Symbiosis is well known to influence bacterial symbiont genome evolution and has recently been shown to shape eukaryotic host genomes. Intriguing patterns of host genome evolution, including remarkable numbers of gene duplications, have been observed in the pea aphid, a sap-feeding insect that relies on a bacterial endosymbiont for amino acid provisioning. Previously, we proposed that gene duplication has been important for the evolution of symbiosis based on aphid-specific gene duplication in amino acid transporters (AATs), with some paralogs highly expressed in the cells housing symbionts (bacteriocytes). Here, we use a comparative approach to test the role of gene duplication in enabling recruitment of AATs to bacteriocytes. Using genomic and transcriptomic data, we annotate AATs from sap-feeding and non sap-feeding insects and find that, like aphids, AAT gene families have undergone independent large-scale gene duplications in three of four additional sap-feeding insects. RNA-seq differential expression data indicate that, like aphids, the sap-feeding citrus mealybug possesses several lineage-specific bacteriocyte-enriched paralogs. Further, differential expression data combined with quantitative PCR support independent evolution of bacteriocyte enrichment in sap-feeding insect AATs. Although these data indicate that gene duplication is not necessary to initiate host/symbiont amino acid exchange, they support a role for gene duplication in enabling AATs to mediate novel host/symbiont interactions broadly in the sap-feeding suborder Sternorrhyncha. In combination with recent studies on other symbiotic systems, gene duplication is emerging as a general pattern in host genome evolution.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Horizontal gene transfer from diverse bacteria to an insect genome enables a tripartite nested mealybug symbiosis

The smallest reported bacterial genome belongs to Tremblaya princeps, a symbiont of Planococcus citri mealybugs (PCIT). Tremblaya PCIT not only has a 139 kb genome, but possesses its own bacterial endosymbiont, Moranella endobia. Genome and transcriptome sequencing, including genome sequencing from a Tremblaya lineage lacking intracellular bacteria, reveals that the extreme genomic degeneracy of Tremblaya PCIT likely resulted from acquiring Moranella as an endosymbiont. In addition, at least 22 expressed horizontally transferred genes from multiple diverse bacteria to the mealybug genome likely complement missing symbiont genes. However, none of these horizontally transferred genes are from Tremblaya, showing that genome reduction in this symbiont has not been enabled by gene transfer to the host nucleus. Our results thus indicate that the functioning of this three-way symbiosis is dependent on genes from at least six lineages of organisms and reveal a path to intimate endosymbiosis distinct from that followed by organelles.

Novel male-biased expression in paralogs of the aphid slimfast nutrient amino acid transporter expansion

Background

A major goal of molecular evolutionary biology is to understand the fate and consequences of duplicated genes. In this context, aphids are intriguing because the newly sequenced pea aphid genome harbors an extraordinary number of lineage-specific gene duplications relative to other insect genomes. Though many of their duplicated genes may be involved in their complex life cycle, duplications in nutrient amino acid transporters appear to be associated rather with their essential amino acid poor diet and the intracellular symbiosis aphids rely on to compensate for dietary deficits. Past work has shown that some duplicated amino acid transporters are highly expressed in the specialized cells housing the symbionts, including a paralog of an aphid-specific expansion homologous to the Drosophila gene slimfast. Previous data provide evidence that these bacteriocyte-expressed transporters mediate amino acid exchange between aphids and their symbionts.

Results

We report that some nutrient amino acid transporters show male-biased expression. Male-biased expression characterizes three paralogs in the aphid-specific slimfast expansion, and the male-biased expression is conserved across two aphid species for at least two paralogs. One of the male-biased paralogs has additionally experienced an accelerated rate of non-synonymous substitutions.

Conclusions

This is the first study to document male-biased slimfast expression. Our data suggest that the male-biased aphid slimfast paralogs diverged from their ancestral function to fill a functional role in males. Furthermore, our results provide evidence that members of the slimfast expansion are maintained in the aphid genome not only for the previously hypothesized role in mediating amino acid exchange between the symbiotic partners, but also for sex-specific roles.

The power of paired genomes

Species interactions are fundamental to ecology. Classic studies of competition, predation, parasitism and mutualism between macroscopic organisms have provided a foundation for the discipline, but many of the most important and intimate ecological interactions are microscopic in scale. These microscopic interactions include those occurring between eukaryotic hosts and their microbial symbionts. Such symbioses, ubiquitous in nature, provide experimental challenges because the partners often cannot live outside the symbiosis. With respect to the symbionts, this precludes utilizing classical microbiological and genetic techniques that require in vitro cultivation. Genomics, however, has rapidly changed the study of symbioses. In this issue of Molecular Ecology, MacDonald et al. (2011), coupling symbiont whole-genome sequencing, experimental studies and metabolic modelling, provide novel insights into one of the best-studied symbioses, that between aphids and their obligate, nutrient-provisioning, intracellular bacteria, Buchnera aphidicola (Fig. 1). MacDonald and colleagues assessed variation in the ability of aphid–Buchnera pairs to thrive on artificial diets missing different amino acids. As shown previously (e.g. Wilkinson & Douglas 2003), aphid–Buchnera pairs can differ in their requirements for external sources of essential amino acids. Such phenotypic variation could result from differences in Buchnera’s amino acid biosynthetic capabilities or in the ability of aphids to interact with their symbionts. Whole-genome sequencing of the Buchnera genomes from four aphid lines with alternate nutritional phenotypes revealed that the environmental nutrients required by the aphid–Buchnera pairs could not be explained by sequence variation in the symbionts. Instead, a novel metabolic modelling approach suggested that much of the variation in nutritional phenotype could be explained by host variation in the capacity to provide necessary nutrient precursors to their symbionts. MacDonald et al.’s work complements a recent study by Vogel & Moran (2011), who through crossing experiments investigating the inheritance of a nutritional phenotype associated with a frameshift mutation in a Buchnera amino acid biosynthesis gene powerfully demonstrated that different host genotypes paired with the same symbiont genome could exhibit substantially different nutritional requirements.† Thus, while there is little doubt that Buchnera are evolutionarily central to the nutritional ecology of aphids, the current work by MacDonald et al. (2011) together with that of Vogel & Moran (2011) surprisingly demonstrates host dominance in defining and controlling the ecological niche of this particular symbiosis.

Aphids alter host-plant nitrogen isotope fractionation

Plant sap-feeding insects and blood-feeding parasites are frequently depleted in (15)N relative to their diet. Unfortunately, most fluid-feeder/host nitrogen stable-isotope studies simply report stable-isotope signatures, but few attempt to elucidate the mechanism of isotopic trophic depletion. Here we address this deficit by investigating the nitrogen stable-isotope dynamics of a fluid-feeding herbivore-host plant system: the green peach aphid, Myzus persicae, feeding on multiple brassicaceous host plants. M. persicae was consistently more than 6‰ depleted in (15)N relative to their hosts, although aphid colonized plants were 1.5‰ to 2.0‰ enriched in (15)N relative to uncolonized control plants. Isotopic depletion of aphids relative to hosts was strongly related to host nitrogen content. We tested whether the concomitant aphid (15)N depletion and host (15)N enrichment was coupled by isotopic mass balance and determined that aphid (15)N depletion and host (15)N enrichment are uncoupled processes. We hypothesized that colonized plants would have higher nitrate reductase activity than uncolonized plants because previous studies had demonstrated that high nitrate reductase activity under substrate-limiting conditions can result in increased plant δ(15)N values. Consistent with our hypothesis, nitrate reductase activity in colonized plants was twice that of uncolonized plants. This study offers two important insights that are likely applicable to understanding nitrogen dynamics in fluid-feeder/host systems. First, isotopic separation of aphid and host depends on nitrogen availability. Second, aphid colonization alters host nitrogen metabolism and subsequently host nitrogen stable-isotope signature. Notably, this work establishes a metabolic framework for future hypothesis-driven studies focused on aphid manipulation of host nitrogen metabolism.

Genome expansion and differential expression of amino acid transporters at the aphid/Buchnera symbiotic interface

In insects, some of the most ecologically important symbioses are nutritional symbioses that provide hosts with novel traits and thereby facilitate exploitation of otherwise inaccessible niches. One such symbiosis is the ancient obligate intracellular symbiosis of aphids with the γ-proteobacteria, Buchnera aphidicola. Although the nutritional basis of the aphid/Buchnera symbiosis is well understood, the processes and structures that mediate the intimate interactions of symbiotic partners remain uncharacterized. Here, using a de novo approach, we characterize the complement of 40 amino acid polyamine organocation (APC) superfamily member amino acid transporters (AATs) encoded in the genome of the pea aphid, Acyrthosiphon pisum. We find that the A. pisum APC superfamily is characterized by extensive gene duplications such that A. pisum has more APC superfamily transporters than other fully sequenced insects, including a ten paralog aphid-specific expansion of the APC transporter slimfast. Detailed expression analysis of 17 transporters selected on the basis of their phylogenetic relationship to five AATs identified in an earlier bacteriocyte expressed sequence tag study distinguished a subset of eight transporters that have been recruited for amino acid transport in bacteriocyte cells at the symbiotic interface. These eight transporters include transporters that are highly expressed and/or highly enriched in bacteriocytes and intriguingly, the four AATs that show bacteriocyte-enriched expression are all members of gene family expansions, whereas three of the four that are highly expressed but not enriched in bacteriocytes retain one-to-one orthology with transporters in other genomes. Finally, analysis of evolutionary rates within the large A. pisum slimfast expansion demonstrated increased rates of molecular evolution coinciding with two major shifts in expression: 1) a loss of gut expression and possibly a gain of bacteriocyte expression and 2) loss of expression in all surveyed tissues in asexual females. Taken together, our characterization of nutrient AATs at the aphid/Buchnera symbiotic interface provides the first examination of the processes and structures operating at the interface of an obligate intracellular insect nutritional symbiosis, offering unique insight into the types of genomic change that likely facilitated evolutionary maintenance of the symbiosis.

The invasion route for an insect pest species: the tobacco aphid in the New World

Biological invasions are rapid evolutionary events in which populations are usually subject to a founder event during introduction followed by rapid adaptation to the new environment. Molecular tools and Bayesian approaches have shown their utility in exploring different evolutionary scenarios regarding the invasion routes of introduced species. We examined the situation for the tobacco aphid, Myzus persicae nicotianae, a recently introduced aphid species in Chile. Using seven microsatellite loci and approximate Bayesian computation, we studied populations of the tobacco aphid sampled from several American and European countries, identifying the most likely source populations and tracking the route of introduction to Chile. Our population genetic data are consistent with available historical information, pointing to an introduction route of the tobacco aphid from Europe and/or from other putative populations (e.g. Asia) with subsequent introduction through North America to South America. Evidence of multiple introductions to North America from different genetic pools, with successive loss of genetic diversity from Europe towards North America and a strong bottleneck during the southward introduction to South America, was also found. Additionally, we examined the special case of a widespread multilocus genotype that was found in all American countries examined. This case provides further evidence for the existence of highly successful genotypes or 'superclones' in asexually reproducing organisms.

Genomic resources for Myzus persicae: EST sequencing, SNP identification, and microarray design

Background

The green peach aphid, Myzus persicae (Sulzer), is a world-wide insect pest capable of infesting more than 40 plant families, including many crop species. However, despite the significant damage inflicted by M. persicae in agricultural systems through direct feeding damage and by its ability to transmit plant viruses, limited genomic information is available for this species.

Results

Sequencing of 16 M. persicae cDNA libraries generated 26,669 expressed sequence tags (ESTs). Aphids for library construction were raised on Arabidopsis thaliana, Nicotiana benthamiana, Brassica oleracea, B. napus, and Physalis floridana (with and without Potato leafroll virus infection). The M. persicae cDNA libraries include ones made from sexual and asexual whole aphids, guts, heads, and salivary glands. In silico comparison of cDNA libraries identified aphid genes with tissue-specific expression patterns, and gene expression that is induced by feeding on Nicotiana benthamiana. Furthermore, 2423 genes that are novel to science and potentially aphid-specific were identified. Comparison of cDNA data from three aphid lineages identified single nucleotide polymorphisms that can be used as genetic markers and, in some cases, may represent functional differences in the protein products. In particular, non-conservative amino acid substitutions in a highly expressed gut protease may be of adaptive significance for M. persicae feeding on different host plants. The Agilent eArray platform was used to design an M. persicae oligonucleotide microarray representing over 10,000 unique genes.

Conclusion

New genomic resources have been developed for M. persicae, an agriculturally important insect pest. These include previously unknown sequence data, a collection of expressed genes, molecular markers, and a DNA microarray that can be used to study aphid gene expression. These resources will help elucidate the adaptations that allow M. persicae to develop compatible interactions with its host plants, complementing ongoing work illuminating plant molecular responses to phloem-feeding insects.

The genetic outcomes of sex and recombination in long-term functionally parthenogenetic lineages of Australian Sitobion aphids

The typical life cycle of an aphid is cyclical parthenogenesis which involves the alternation of sexual and asexual reproduction. However, aphid life cycles, even within a species, can encompass everything on a continuum from obligate sexuality, through facultative sexuality to obligate asexuality. Loss of the sexual cycle in aphids is frequently associated with the introduction of a new pest and can occur for a number of environmental and genetic reasons. Here we investigate loss of sexual function in Sitobion aphids in Australia. Specifically, we aimed to determine whether an absence of sexual reproduction in Australian Sitobion results from genetic loss of sexual function or environmental constraints in the introduced range. We addressed our aims by performing a series of breeding experiments. We found that some lineages have genetically lost sexual function while others retain sexual function and appear environmentally constrained to asexuality. Further, in our crosses, using autosomal and X-linked microsatellite markers, we identified processes deviating from normal Mendelian segregation. We observed strong deviations in X chromosome transmission through the sexual cycle. Additionally, when progeny genotypes were examined across multiple loci simultaneously we found that some multilocus genotypes are significantly over-represented in the sample and that levels of heterozygosity were much higher than expected at almost all loci. This study demonstrates that strong biases in the transmission of X chromosomes through the sexual cycle are likely to be widespread in aphids. The mechanisms underlying these patterns are not clear. We discuss several possible alternatives, including mutation accumulation during periods of functional asexuality and genetic imprinting.

Microsatellites reveal male recombination and neo-sex chromosome formation in Scaptodrosophila hibisci (Drosophilidae)

In drosophilid flies, male recombination and neo-sex chromosome formation are rare. Following the genotyping of full-sib families with 20 microsatellite markers and subsequent cytological work, we found evidence of both male recombination and neo-sex chromosome formation in Scaptodrosophila hibisci. As far as we are aware, this is the first report of male recombination and neo-sex chromosome formation co-occurring in a drosophilid fly. Two autosomal loci, Sh29c and Sh90, showed aberrant segregation of male parental alleles. We describe how an autosomal fission followed by fusion of one of the autosomal fragments to the Y chromosome to create a Y1Y2X1X2/X1X1X2X2 sex determination system provides the most parsimonious explanation of the patterns we observe. Male recombination was observed in three families, including autosomal linkage groups and the Y1/X2 linkage group. In addition to the X1 linkage group, two autosomal linkage groups were identified.

Lack of detectable genetic recombination on the X chromosome during the parthenogenetic production of female and male aphids

We used polymorphic microsatellite markers to look for recombination during parthenogenetic oogenesis between the X chromosomes of aphids of the tribe Macrosiphini. We examined the X chromosome because it comprises approximately 25 % of the genome and previous cytological observations of chromosome pairing and nucleolar organizer (NOR) heteromorphism suggest recombination, although the same is not true for autosomes. A total of 564 parthenogenetic females of Myzus clones with three distinct reproductive modes (cyclical parthenogenesis, obligate parthenogenesis and obligate parthenogenesis with male production) were genotyped at three informative X-linked loci. Also, parthenogenetically produced males from clones encompassing the full range of male-producing reproductive strategies were genotyped. These included 391 Myzus persicae males that were genotyped at three X-linked loci and 538 males from Sitobion clones that were genotyped at five informative X-linked loci. Our results show no departure from clonality in parthenogenetic generations of aphids of the tribe Macrosiphini: no recombinant genotypes were observed in parthenogenetically produced males or females.

Segregation of autosomes during spermatogenesis in the peach-potato aphid (Myzus persicae) (Sulzer) (Hemiptera: Aphididae)

Most aphids are cyclic parthenogens, so are ideal models in studies of the mechanisms and consequences of sex and recombination. However, owing to a shortage of physical and genetic markers, there have been few studies of the most fundamental genetic processes in these organisms. For example, it is not known whether autosomal segregation during male spermatogenesis is in Mendelian proportions: we address that question here. The aphid Myzus persicae has a typical karyotype of 2n = 12 in females (XX), while males are XO (2n = 11). During male meiosis, only the spermatocytes with an X chromosome are viable. We hypothesized that assortment of autosomes might be non-random because chromosomal imprinting leading to elimination of the paternal autosomes is seen in the closely related coccoids. In other aphid models, we have observed segregation distortions at single microsatellite loci (Wilson, 2000). Such distortions may have nothing to do with 'selfish' behaviour, but may be caused by mutation accumulation causing fitness differentials. Thus single-locus distortions might be predicted to be more likely to be detected via the male lines of clones that have lost the ability to reproduce sexually (male-producing obligate parthenogenesis (androcyclic)). Using microsatellites we show that genetic imprinting or selfish autosome behaviour does not occur in male M. persicae. Generally, loci segregated in Mendelian proportions in both sexes of cyclically parthenogenetic (holocyclic) clones. However, in androcyclic clones, segregation distortions consistently involved the same two autosomes. This is consistent with linkage of markers to deleterious mutations associated with a loss of sexual reproduction.

Microsatellite variation in cyclically parthenogenetic populations of Myzus persicae in south-eastern Australia

We examined the population structure of the introduced aphid, Myzus persicae collected mainly from its primary host, Prunus persica, in south-east Australia. Myzus persicae has been present in Australia since at least 1893. Samples were collected in the spring of 1998 from two mainland and three Tasmanian localities and isofemale lines were established in the laboratory. The reproductive mode (life cycle), karyotype and 17-locus microsatellite genotype of each clone were determined. All populations showed significant population differentiation (F(ST) 0.058-0.202) even over small geographic distances (<50 km). All clones were karyotypically normal except for a subset of clones from one site that was exposed to the carbamate insecticide, Pirimor, the week prior to sampling. Those clones were heterozygous for an autosomal 1,3 translocation frequently associated in M. persicae with insecticide resistance. In contrast to other loci and despite being on different chromosomes, loci myz2(A) and M55(A) showed general and significant linkage disequilibrium. These loci may be affected by epistatic selection. We discuss the observed high clonal diversity, moderate but significant population differentiation, general conformance to Hardy-Weinberg equilibria and low linkage disequilibria with particular focus on the global population biology of M. persicae.