Abiotic and biotic stressors causing equivalent mortality induce highly variable transcriptional responses in the soybean aphid

Exposure to heat stress leads to striking clone-specific nymph deformity in pea aphid

Climatic changes, such as heatwaves, pose unprecedented challenges for insects, as escalated temperatures above the thermal optimum alter insect reproductive strategies and energy metabolism. While thermal stress responses have been reported in different insect species, thermo-induced developmental abnormalities in phloem-feeding pests are largely unknown. In this laboratory study, we raised two groups of first instar nymphs belonging to two clones of the pea aphid Acyrthosiphon pisum, on fava beans Vicia faba. The instars developed and then asexually reproduced under constant exposure to a sub-lethal heatwave (27°C) for 14 days. Most mothers survived but their progenies showed abnormalities, as stillbirths and appendageless or weak nymphs with folded appendages were delivered. Clone N116 produced more deceased and appendageless embryos, contrary to N127, which produced fewer dead and more malformed premature embryos. Interestingly, the expression of the HSP70 and HSP83 genes differed in mothers between the clones. Moreover, noticeable changes in metabolism, e.g., lipids, were also detected and that differed in response to stress. Deformed offspring production after heat exposure may be due to heat injury and differential HSP gene expression, but may also be indicative of a conflict between maternal and offspring fitness. Reproductive altruism might have occurred to ensure some of the genetically identical daughters survive. This is because maintaining homeostasis and complete embryogenesis could not be simultaneously fulfilled due to the high costs of stress. Our findings shine new light on pea aphid responses to heatwaves and merit further examination across different lineages and species.

Assessment of avian health status: suitability and constraints of the Zoetis VetScan VS2 blood analyser for ecological and evolutionary studies

Biochemical analyses of blood can decipher physiological conditions of living animals and unravel mechanistic underpinnings of life-history strategies and trade-offs. Yet, researchers in ecology and evolution often face constraints in which methods to apply, not least due to blood volume restrictions or field settings. Here, we test the suitability of a portable biochemical analyser (Zoetis VetScan VS2) for ecological and evolutionary studies that may help solve those problems. Using as little as 80 µl of whole-bird blood from free-living Jackdaws (Corvus monedula) and captive Zebra Finches (Taeniopygia guttata), we show that eight (out of 10) blood analytes show high repeatability after short-term storage (approximately 2 h) and six after 12 h storage time. Handling stress had a clear impact on all except two analytes by 16 min after catching. Finally, six analytes showed consistency within individuals over a period of 30 days, and three even showed individual consistency over a year. Taken together, we conclude that the VetScan VS2 captures biologically relevant variation in blood analytes using just 80 µl of whole blood and, thus, provides valuable physiological measurements of (small) birds sampled in semi-field and field conditions.

Bacteriocyte plasticity in pea aphids facing amino acid stress or starvation during development

An important contributing factor to the evolutionary success of insects is nutritional association with microbial symbionts, which provide the host insects with nutrients lacking in their unbalanced diets. These symbionts are often compartmentalized in specialized cells of the host, the bacteriocytes. Even though bacteriocytes were first described more than a century ago, few studies have explored their dynamics throughout the insect life cycle and in response to environmental stressors. Here, we use the Buchnera aphidicola/pea aphid symbiotic system to study how bacteriocytes are regulated in response to nutritional stress throughout aphid development. Using artificial diets, we analyzed the effects of depletion or excess of phenylalanine or leucine, two amino acids essential for aphid growth and whose biosynthetic pathways are shared between the host and the symbiont. Bacteriocytes responded dynamically to those treatments, while other tissues showed no obvious morphological change. Amino acid depletion resulted in an increase in bacteriocyte numbers, with the extent of the increase depending on the amino acid, while excess either caused a decrease (for leucine) or an increase (for phenylalanine). Only a limited impact on survival and fecundity was observed, suggesting that the adjustment in bacteriocyte (and symbiont) numbers is sufficient to withstand these nutritional challenges. We also studied the impact of more extreme conditions by exposing aphids to a 24 h starvation period at the beginning of nymphal development. This led to a dramatic drop in aphid survival and fecundity and a significant developmental delay. Again, bacteriocytes responded dynamically, with a considerable decrease in number and size, correlated with a decrease in the number of symbionts, which were prematurely degraded by the lysosomal system. This study shows how bacteriocyte dynamics is integrated in the physiology of insects and highlights the high plasticity of these cells.

Evolved transcriptional responses and their trade-offs after long-term adaptation of Bemisia tabaci to a marginally-suitable host

Although generalist insect herbivores can migrate and rapidly adapt to a broad range of host plants, they can face significant difficulties when accidentally migrating to novel and marginally-suitable hosts. What happens, both in performance and gene expression regulation, if these marginally-suitable hosts must be used for multiple generations before migration to a suitable host can take place, largely remains unknown. In this study, we established multigenerational colonies of the whitefly Bemisia tabaci, a generalist phloem-feeding species, adapted to a marginally-suitable host (habanero pepper) or an optimal host (cotton). We used reciprocal host tests to estimate the differences in performance of the populations on both hosts under optimal (30 oC) and mild-stressful (24 oC) temperature conditions, and documented the associated transcriptomic changes. The habanero pepper-adapted population greatly improved its performance on habanero pepper but did not reach its performance level on cotton, the original host. It also showed reduced performance on cotton, relative to the non-adapted population, and an antagonistic effect of the lower-temperature stressor. The transcriptomic data revealed that most of the expression changes, associated with long-term adaptation to habanero pepper, can be categorized as "evolved" with no initial plastic response. Three molecular functions dominated: enhanced formation of cuticle structural constituents, enhanced activity of oxidation-reduction processes involved in neutralization of phytotoxins and reduced production of proteins from the cathepsin B family. Taken together, these findings indicate that generalist insects can adapt to novel host plants by modifying the expression of a relatively small set of specific molecular functions.

Potato leafroll virus reduces Buchnera aphidocola titer and alters vector transcriptome responses

Viruses in the Luteoviridae family, such as Potato leafroll virus (PLRV), are transmitted by aphids in a circulative and nonpropagative mode. This means the virions enter the aphid body through the gut when they feed from infected plants and then the virions circulate through the hemolymph to enter the salivary glands before being released into the saliva. Although these viruses do not replicate in their insect vectors, previous studies have demonstrated viruliferous aphid behavior is altered and the obligate symbiont of aphids, Buchnera aphidocola, may be involved in transmission. Here we provide the transcriptome of green peach aphids (Myzus persicae) carrying PLRV and virus-free control aphids using Illumina sequencing. Over 150 million paired-end reads were obtained through Illumina sequencing, with an average of 19 million reads per library. The comparative analysis identified 134 differentially expressed genes (DEGs) between the M. persicae transcriptomes, including 64 and 70 genes that were up- and down-regulated in aphids carrying PLRV, respectively. Using functional classification in the GO databases, 80 of the DEGs were assigned to 391 functional subcategories at category level 2. The most highly up-regulated genes in aphids carrying PLRV were cytochrome p450s, genes related to cuticle production, and genes related to development, while genes related to heat shock proteins, histones, and histone modification were the most down-regulated. PLRV aphids had reduced Buchnera titer and lower abundance of several Buchnera transcripts related to stress responses and metabolism. These results suggest carrying PLRV may reduce both aphid and Buchnera genes in response to stress. This work provides valuable basis for further investigation into the complicated mechanisms of circulative and nonpropagative transmission.

Whole genome sequence of the soybean aphid, Aphis glycines

Aphids are emerging as model organisms for both basic and applied research. Of the 5,000 estimated species, only three aphids have published whole genome sequences: the pea aphid Acyrthosiphon pisum, the Russian wheat aphid, Diuraphis noxia, and the green peach aphid, Myzus persicae. We present the whole genome sequence of a fourth aphid, the soybean aphid (Aphis glycines), which is an extreme specialist and an important invasive pest of soybean (Glycine max). The availability of genomic resources is important to establish effective and sustainable pest control, as well as to expand our understanding of aphid evolution. We generated a 302.9 Mbp draft genome assembly for Ap. glycines using a hybrid sequencing approach. This assembly shows high completeness with 19,182 predicted genes, 92% of known Ap. glycines transcripts mapping to contigs, and substantial continuity with a scaffold N50 of 174,505 bp. The assembly represents 95.5% of the predicted genome size of 317.1 Mbp based on flow cytometry. Ap. glycines contains the smallest known aphid genome to date, based on updated genome sizes for 19 aphid species. The repetitive DNA content of the Ap. glycines genome assembly (81.6 Mbp or 26.94% of the 302.9 Mbp assembly) shows a reduction in the number of classified transposable elements compared to Ac. pisum, and likely contributes to the small estimated genome size. We include comparative analyses of gene families related to host-specificity (cytochrome P450's and effectors), which may be important in Ap. glycines evolution. This Ap. glycines draft genome sequence will provide a resource for the study of aphid genome evolution, their interaction with host plants, and candidate genes for novel insect control methods.

Phenotypic and Transcriptomic Response of the Grasshopper Oedaleus asiaticus (Orthoptera: Acrididae) to Toxic Rutin

Rutin, a widely distributed phytochemical flavonoid, can be used to control insect pests. In this study, we studied the growth performance of the grasshopper Oedaleus asiaticus Bey-Bienko given xenobiotic rutin using feeding experiments and transcriptomic analysis. O. asiaticus had reduced body size, lower survival rate, and reduced growth performance when fed with xenobiotic rutin. Rutin-fed nymphs had large variation in gene expression profiles, with a total of 308 genes significantly upregulated and 287 genes downregulated. The upregulated genes were significantly enriched in stress resistance-, immune-, and detoxification-related biological processes and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Downregulated genes mainly involved cuticle biosynthesis and nutrition metabolism-related pathways. The quantitative real-time PCR (qRT-PCR) analysis of 15 candidate genes also produced results consistent with the transcriptome data. These results suggested that grasshoppers’ capacity for biosynthesis and nutrition metabolism decreased, and stress resistance and metabolized capacity to toxic substances were significantly induced when O. asiaticus was fed on xenobiotic rutin. Rutin, as a phytotoxin, had detrimental effects and induced changes in gene expression profiles for O. asiaticus. This study can provide a molecular basis and offer future opportunities for the development of rutin-related insecticides and their application to grasshopper control.

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.

Divergent Switchgrass Cultivars Modify Cereal Aphid Transcriptomes

Schizaphis graminum Rondani (Hemiptera: Aphididae) and Sipha flava Forbes (Hemiptera: Aphididae) are two common pests of bioenergy grasses. Despite the fact that they are both considered generalists, they differ in their ability to colonize Panicum virgatum cultivars. For example, S. flava colonizes both P. virgatum cv. Summer and P. virgatum cv. Kanlow whereas S. graminum can only colonize Summer. To study the molecular responses of these aphids to these two switchgrass cultivars, we generated de novo transcriptome assemblies and compared the expression profiles of aphids feeding on both cultivars to profiles associated with feeding on a highly susceptible sorghum host and a starvation treatment. Transcriptome assemblies yielded 8,428 and 8,866 high-quality unigenes for S. graminum and S. flava, respectively. Overall, S. graminum responded strongly to all three treatments after 12 h with an upregulation of unigenes coding for detoxification enzymes while major transcriptional changes were not observed in S. flava until 24 h. Additionally, while the two aphids responded to the switchgrass feeding treatment by downregulating unigenes linked to growth and development, their responses to Summer and Kanlow diverged significantly. Schizaphis graminum upregulated more unigenes coding for stress-responsive enzymes in the Summer treatment compared to S. flava; however, many of these unigenes were actually downregulated in the Kanlow treatment. In contrast, S. flava appeared capable of overcoming host defenses by upregulating a larger number of unigenes coding for detoxification enzymes in the Kanlow treatment. Overall, these findings are consistent with previous studies on the interactions of these two cereal aphids to divergent switchgrass hosts.

Trichoderma atroviride P1 Colonization of Tomato Plants Enhances Both Direct and Indirect Defense Barriers Against Insects

Numerous microbial root symbionts are known to induce different levels of enhanced plant protection against a variety of pathogens. However, more recent studies have demonstrated that beneficial microbes are able to induce plant systemic resistance that confers some degree of protection against insects. Here, we report how treatments with the fungal biocontrol agent Trichoderma atroviride strain P1 in tomato plants induce responses that affect pest insects with different feeding habits: the noctuid moth Spodoptera littoralis (Boisduval) and the aphid Macrosiphum euphorbiae (Thomas). We observed that the tomato plant-Trichoderma P1 interaction had a negative impact on the development of moth larvae and on aphid longevity. These effects were attributed to a plant response induced by Trichoderma that was associated with transcriptional changes of a wide array of defense-related genes. While the impact on aphids could be related to the up-regulation of genes involved in the oxidative burst reaction, which occur early in the defense reaction, the negative performance of moth larvae was associated with the enhanced expression of genes encoding for protective enzymes (i.e., Proteinase inhibitor I (PI), Threonine deaminase, Leucine aminopeptidase A1, Arginase 2, and Polyphenol oxidase) that are activated downstream in the defense cascade. In addition, Trichoderma P1 produced alterations in plant metabolic pathways leading to the production and release of volatile organic compounds (VOCs) that are involved in the attraction of the aphid parasitoid Aphidius ervi, thus reinforcing the indirect plant defense barriers. Our findings, along with the evidence available in the literature, indicate that the outcome of the tripartite interaction among plant, Trichoderma, and pests is highly specific and only a comprehensive approach, integrating both insect phenotypic changes and plant transcriptomic alterations, can allow a reliable prediction of its potential for plant protection.

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.

Learning to starve: impacts of food limitation beyond the stress period

Starvation is common among wild animal populations, and many individuals experience repeated bouts of starvation over the course of their lives. Although much information has been gained through laboratory studies of acute starvation, little is known about how starvation affects an animal once food is again available (i.e. during the refeeding and recovery phases). Many animals exhibit a curious phenomenon – some seem to ‘get better’ at starving following exposure to one or more starvation events – by this we mean that they exhibit potentially adaptive responses, including reduced rates of mass loss, reduced metabolic rates, and lower costs of digestion. During subsequent refeedings they may also exhibit improved digestive efficiency and more rapid mass gain. Importantly, these responses can last until the next starvation bout or even be inherited and expressed in the subsequent generation. Currently, however, little is known about the molecular regulation and physiological mechanisms underlying these changes. Here, we identify areas of research that can fill in the most pressing knowledge gaps. In particular, we highlight how recently refined techniques (e.g. stable isotope tracers, quantitative magnetic resonance and thermal measurement) as well as next-generation sequencing approaches (e.g. RNA-seq, proteomics and holobiome sequencing) can address specific starvation-focused questions. We also describe outstanding unknowns ripe for future research regarding the timing and severity of starvation, and concerning the persistence of these responses and their interactions with other ecological stressors.

Transcriptional profile and differential fitness in a specialist milkweed insect across host plants varying in toxicity

Interactions between plants and herbivorous insects have been models for theories of specialization and co-evolution for over a century. Phytochemicals govern many aspects of these interactions and have fostered the evolution of adaptations by insects to tolerate or even specialize on plant defensive chemistry. While genomic approaches are providing new insights into the genes and mechanisms insect specialists employ to tolerate plant secondary metabolites, open questions remain about the evolution and conservation of insect counterdefences, how insects respond to the diversity defences mounted by their host plants, and the costs and benefits of resistance and tolerance to plant defences in natural ecological communities. Using a milkweed-specialist aphid (Aphis nerii) model, we test the effects of host plant species with increased toxicity, likely driven primarily by increased secondary metabolites, on aphid life history traits and whole-body gene expression. We show that more toxic plant species have a negative effect on aphid development and lifetime fecundity. When feeding on more toxic host plants with higher levels of secondary metabolites, aphids regulate a narrow, targeted set of genes, including those involved in canonical detoxification processes (e.g., cytochrome P450s, hydrolases, UDP-glucuronosyltransferases and ABC transporters). These results indicate that A. nerii marshal a variety of metabolic detoxification mechanisms to circumvent milkweed toxicity and facilitate host plant specialization, yet, despite these detoxification mechanisms, aphids experience reduced fitness when feeding on more toxic host plants. Disentangling how specialist insects respond to challenging host plants is a pivotal step in understanding the evolution of specialized diet breadths.

Diet alters performance and transcription patterns in Oedaleus asiaticus (Orthoptera: Acrididae) grasshoppers

We reared Oedaleus asiaticus grasshoppers under four different single-plant diets to examine the relationships among diet, performance, stress, and transcription patterns. Grasshoppers fed only Artemisia frigida (Asteraceae) were stressed, as indicated by their lower growth, size, development, and survival, in comparison to grasshoppers fed on any of three grasses, Cleistogenes squarrosa, Leymus chinensis, or Stipa krylovii (all Poaceae). We then used transcriptome analysis to examine how gene expression levels in O. asiaticus were altered by feeding on these diets. Nymphs fed A. frigida had the largest variation in gene expression profiles with a total of 299 genes significantly up- or down-regulated compared to those feeding on the three grasses: down-regulated genes included those involved in cuticle biosynthesis, DNA replication, biosynthesis and metabolism of nutrition. The up-regulated genes included stress-resistant and detoxifying enzymes. GO and KEGG enrichment analysis also showed that feeding on A. frigida could down-regulate biosynthesis and metabolism related pathways, and up-regulate stress-resistant and detoxification terms and pathways. Our results show that diet significantly altered gene-expression, and that unfavorable, stressful diets induce more transcriptional changes than favorable diets. Altered gene-expression represents phenotypic plasticity, and many such changes appear to be evolved, adaptive responses. The ease and regularity by which individuals shift phenotypes via altered transcription suggests that populations consist not of similar, fixed phenotypes, but of a collection of ever-changing, divergent phenotypes.

Molecular Ecological Basis of Grasshopper (Oedaleus asiaticus) Phenotypic Plasticity under Environmental Selection

While ecological adaptation in insects can be reflected by plasticity of phenotype, determining the causes and molecular mechanisms for phenotypic plasticity (PP) remains a crucial and still difficult question in ecology, especially where control of insect pests is involved. Oedaleus asiaticus is one of the most dominant pests in the Inner Mongolia steppe and represents an excellent system to study phenotypic plasticity. To better understand ecological factors affecting grasshopper phenotypic plasticity and its molecular control, we conducted a full transcriptional screening of O. asiaticus grasshoppers reared in four different grassland patches in Inner Mongolia. Grasshoppers showed different degrees of PP associated with unique gene expressions and different habitat plant community compositions. Grasshopper performance variables were susceptible to habitat environment conditions and closely associated with plant architectures. Intriguingly, eco-transcriptome analysis revealed five potential candidate genes playing important roles in grasshopper performance, with gene expression closely relating to PP and plant community factors. By linking the grasshopper performances to gene profiles and ecological factors using canonical regression, we first demonstrated the eco-transcriptomic architecture (ETA) of grasshopper phenotypic traits (ETAGPTs). ETAGPTs revealed plant food type, plant density, coverage, and height were the main ecological factors influencing PP, while insect cuticle protein (ICP), negative elongation factor A (NELFA), and lactase-phlorizin hydrolase (LCT) were the key genes associated with PP. Our study gives a clear picture of gene-environment interaction in the formation and maintenance of PP and enriches our understanding of the transcriptional events underlying molecular control of rapid phenotypic plasticity associated with environmental variability. The findings of this study may also provide new targets for pest control and highlight the significance of ecological management practice on grassland conservation.

Biology, physiology and gene expression of grasshopper Oedaleus asiaticus exposed to diet stress from plant secondary compounds

We studied the role of plant primary and secondary metabolites in mediating plant-insect interactions by conducting a no-choice single-plant species field experiment to compare the suitability, enzyme activities, and gene expression of Oedaleus asiaticus grasshoppers feeding on four host and non-host plants with different chemical traits. O. asiaticus growth showed a positive relationship to food nutrition content and a negative relationship to secondary compounds content. Grasshopper amylase, chymotrypsin, and lipase activities were positively related to food starch, crude protein, and lipid content, respectively. Activity of cytochrome P450s, glutathione-S-transferase, and carboxylesterase were positively related to levels of secondary plant compounds. Gene expression of UDP-glucuronosyltransferase 2C1, cytochrome P450 6K1 were also positively related to secondary compounds content in the diet. Grasshoppers feeding on Artemisia frigida, a species with low nutrient content and a high level of secondary compounds, had reduced growth and digestive enzyme activity. They also had higher detoxification enzyme activity and gene expression compared to grasshoppers feeding on the grasses Cleistogenes squarrosa, Leymus chinensis, or Stipa krylovii. These results illustrated Oedaleus asiaticus adaptive responses to diet stress resulting from toxic chemicals, and support the hypothesis that nutritious food benefits insect growth, but plant secondary compounds are detrimental for insect growth.

The genomewide transcriptional response underlying the pea aphid wing polyphenism

Phenotypic plasticity is a key life history strategy used by many plants and animals living in heterogeneous environments. A multitude of studies have investigated the costs and limits of plasticity, as well as the conditions under which it evolves. Much less well understood are the molecular genetic mechanisms that enable an organism to sense its environment and respond in a plastic manner. The pea aphid wing polyphenism is a compelling laboratory model to study these mechanisms. In this polyphenism, environmental stressors like high density cause asexual, viviparous adult female aphids to change the development of their embryos from wingless to winged morphs. The life history trade-offs between the two morphs have been intensively studied, but the molecular mechanisms underlying this process remain largely unknown. We therefore performed a genomewide study of the maternal transcriptome at two time points with and without a crowding stress to discover the maternal molecular changes that lead to the development of winged vs. wingless offspring. We observed significant transcriptional changes in genes associated with odorant binding, neurotransmitter transport, hormonal activity and chromatin remodelling in the maternal transcriptome. We also found that titres of serotonin, dopamine and octopamine were higher in solitary compared to crowded aphids. We use these results to posit a model for how maternal signals inform a developing embryo to be winged or wingless. Our findings add significant insights into the identity of the molecular mechanisms that underlie environmentally induced morph determination and suggest a possible role for biogenic amine regulation in polyphenisms generally.

Stress-induced changes in abundance differ among obligate and facultative endosymbionts of the soybean aphid

Bacterial endosymbionts can drive evolutionary novelty by conferring adaptive benefits under adverse environmental conditions. Among aphid species there is growing evidence that symbionts influence tolerance to various forms of stress. However, the extent to which stress inflicted on the aphid host has cascading effects on symbiont community dynamics remains poorly understood. Here we simultaneously quantified the effect of host-plant induced and xenobiotic stress on soybean aphid (Aphis glycines) fitness and relative abundance of its three bacterial symbionts. Exposure to soybean defensive stress (Rag1 gene) and a neurotoxic insecticide (thiamethoxam) substantially reduced aphid composite fitness (survival × reproduction) by 74 ± 10% and 92 ± 2%, respectively, which in turn induced distinctive changes in the endosymbiont microbiota. When challenged by host-plant defenses a 1.4-fold reduction in abundance of the obligate symbiont Buchnera was observed across four aphid clonal lines. Among facultative symbionts of Rag1-stressed aphids, Wolbachia abundance increased twofold and Arsenophonus decreased 1.5-fold. A similar pattern was observed under xenobiotic stress, with Buchnera and Arsenophonus titers decreasing (1.3-fold) and Wolbachia increasing (1.5-fold). Furthermore, variation in aphid virulence to Rag1 was positively correlated with changes in Arsenophonus titers, but not Wolbachia or Buchnera. A single Arsenophonus multi-locus genotype was found among aphid clonal lines, indicating strain diversity is not primarily responsible for correlated host-symbiont stress levels. Overall, our results demonstrate the nature of aphid symbioses can significantly affect the outcome of interactions under stress and suggests general changes in the microbiome can occur across multiple stress types.