PacBio Hi-Fi genome assembly of Sipha maydis, a model for the study of multipartite mutualism in insects.

Dependence on multiple nutritional endosymbionts has evolved repeatedly in insects feeding on unbalanced diets. However, reference genomes for species hosting multi-symbiotic nutritional systems are lacking, even though they are essential for deciphering the processes governing cooperative life between insects and anatomically integrated symbionts. The cereal aphid Sipha maydis is a promising model for addressing these issues, as it has evolved a nutritional dependence on two bacterial endosymbionts that complement each other. In this study, we used PacBio High fidelity (HiFi) long-read sequencing to generate a highly contiguous genome assembly of S. maydis with a length of 410 Mb, 3,570 contigs with a contig N50 length of 187 kb, and BUSCO completeness of 95.5%. We identified 117 Mb of repetitive sequences, accounting for 29% of the genome assembly, and predicted 24,453 protein-coding genes, of which 2,541 were predicted enzymes included in an integrated metabolic network with the two aphid-associated endosymbionts. These resources provide valuable genetic and metabolic information for understanding the evolution and functioning of multi-symbiotic systems in insects.

MicroAnnot: A Dedicated Workflow for Accurate Microsporidian Genome Annotation.

With nearly 1700 species, Microsporidia represent a group of obligate intracellular eukaryotes with veterinary, economic and medical impacts. To help understand the biological functions of these microorganisms, complete genome sequencing is routinely used. Nevertheless, the proper prediction of their gene catalogue is challenging due to their taxon-specific evolutionary features. As innovative genome annotation strategies are needed to obtain a representative snapshot of the overall lifestyle of these parasites, the MicroAnnot tool, a dedicated workflow for microsporidian sequence annotation using data from curated databases of accurately annotated microsporidian genes, has been developed. Furthermore, specific modules have been implemented to perform small gene (<300 bp) and transposable element identification. Finally, functional annotation was performed using the signature-based InterProScan software. MicroAnnot's accuracy has been verified by the re-annotation of four microsporidian genomes for which structural annotation had previously been validated. With its comparative approach and transcriptional signal identification method, MicroAnnot provides an accurate prediction of translation initiation sites, an efficient identification of transposable elements, as well as high specificity and sensitivity for microsporidian genes, including those under 300 bp.

Coordination of host and endosymbiont gene expression governs endosymbiont growth and elimination in the cereal weevil Sitophilus spp.

BACKGROUND: Insects living in nutritionally poor environments often establish long-term relationships with intracellular bacteria that supplement their diets and improve their adaptive and invasive powers. Even though these symbiotic associations have been extensively studied on physiological, ecological, and evolutionary levels, few studies have focused on the molecular dialogue between host and endosymbionts to identify genes and pathways involved in endosymbiosis control and dynamics throughout host development. RESULTS: We simultaneously analyzed host and endosymbiont gene expression during the life cycle of the cereal weevil Sitophilus oryzae, from larval stages to adults, with a particular emphasis on emerging adults where the endosymbiont Sodalis pierantonius experiences a contrasted growth-climax-elimination dynamics. We unraveled a constant arms race in which different biological functions are intertwined and coregulated across both partners. These include immunity, metabolism, metal control, apoptosis, and bacterial stress response. CONCLUSIONS: The study of these tightly regulated functions, which are at the center of symbiotic regulations, provides evidence on how hosts and bacteria finely tune their gene expression and respond to different physiological challenges constrained by insect development in a nutritionally limited ecological niche. Video Abstract.

Differential responses of Ceratitis capitata to infection by the entomopathogenic fungus Purpureocillium lilacinum.

The medfly Ceratitis capitata is one of the most damaging fruit pests with quarantine significance due to its extremely wide host range. The use of entomopathogenic fungi constitutes a promising approach with potential applications in integrated pest management. Furthermore, developing insect control methods can involve the use of fungal machinery to cause metabolic disruption, which may increase its effectiveness by impairing insect development. Insect species, including C. capitata, relies on reproduction potential, nutrient reserves, metabolic activities, and immune response for survival. Accordingly, the purpose of this study was to investigate the impacts of the entomopathogenic fungus Purpureocillium lilacinum on C. capitata pre-mortality. The medfly V8 strain was subjected to laboratory bioassays, which consisted on determining the virulence of P. lilacinum on the medfly. Purpureocillium lilacinum was applied on abdominal topical of 5-day-old males and females. Following the fungal inoculation, we have confirmed (i) a significant increase in tissue sugar content, (ii) a significant decrease in carbohydrase activities, digestive glycosyl hydrolase, and proteinase activities in whole midguts of treated flies, (iii) the antimicrobial peptides (AMPs) genes expression profile was significantly influenced by fly gender, fly status (virgin, mature, and mated), and time after infection, but infection itself had no discernible impact on the AMPs for the genes that were examined. This study provides the first insight into how P. lilacinum could affect C. capitata physiological mechanisms and provides the foundation for considering P. lilacinum as a novel, promising biocontrol agent.

Endosymbiont-containing germarium transcriptional survey in a cereal weevil depicts downregulation of immune effectors at the onset of sexual maturity.

Insects often establish long-term relationships with intracellular symbiotic bacteria, i.e., endosymbionts, that provide them with essential nutrients such as amino acids and vitamins. Endosymbionts are typically confined within specialized host cells called bacteriocytes that may form an organ, the bacteriome. Compartmentalization within host cells is paramount for protecting the endosymbionts and also avoiding chronic activation of the host immune system. In the cereal weevil Sitophilus oryzae, bacteriomes are present as a single organ at the larval foregut-midgut junction, and in adults, at the apex of midgut mesenteric caeca and at the apex of the four ovarioles. While the adult midgut endosymbionts experience a drastic proliferation during early adulthood followed by complete elimination through apoptosis and autophagy, ovarian endosymbionts are maintained throughout the weevil lifetime by unknown mechanisms. Bacteria present in ovarian bacteriomes are thought to be involved in the maternal transmission of endosymbionts through infection of the female germline, but the exact mode of transmission is not fully understood. Here, we show that endosymbionts are able to colonize the germarium in one-week-old females, pinpointing a potential infection route of oocytes. To identify potential immune regulators of ovarian endosymbionts, we have analyzed the transcriptomes of the ovarian bacteriomes through young adult development, from one-day-old adults to sexually mature ones. In contrast with midgut bacteriomes, immune effectors are downregulated in ovarian bacteriomes at the onset of sexual maturation. We hypothesize that relaxation of endosymbiont control by antimicrobial peptides might allow bacterial migration and potential oocyte infection, ensuring endosymbiont transmission.

GReNaDIne: A Data-Driven Python Library to Infer Gene Regulatory Networks from Gene Expression Data.

Context: Inferring gene regulatory networks (GRN) from high-throughput gene expression data is a challenging task for which different strategies have been developed. Nevertheless, no ever-winning method exists, and each method has its advantages, intrinsic biases, and application domains. Thus, in order to analyze a dataset, users should be able to test different techniques and choose the most appropriate one. This step can be particularly difficult and time consuming, since most methods' implementations are made available independently, possibly in different programming languages. The implementation of an open-source library containing different inference methods within a common framework is expected to be a valuable toolkit for the systems biology community. Results: In this work, we introduce GReNaDIne (Gene Regulatory Network Data-driven Inference), a Python package that implements 18 machine learning data-driven gene regulatory network inference methods. It also includes eight generalist preprocessing techniques, suitable for both RNA-seq and microarray dataset analysis, as well as four normalization techniques dedicated to RNA-seq. In addition, this package implements the possibility to combine the results of different inference tools to form robust and efficient ensembles. This package has been successfully assessed under the DREAM5 challenge benchmark dataset. The open-source GReNaDIne Python package is made freely available in a dedicated GitLab repository, as well as in the official third-party software repository PyPI Python Package Index. The latest documentation on the GReNaDIne library is also available at Read the Docs, an open-source software documentation hosting platform. Contribution: The GReNaDIne tool represents a technological contribution to the field of systems biology. This package can be used to infer gene regulatory networks from high-throughput gene expression data using different algorithms within the same framework. In order to analyze their datasets, users can apply a battery of preprocessing and postprocessing tools and choose the most adapted inference method from the GReNaDIne library and even combine the output of different methods to obtain more robust results. The results format provided by GReNaDIne is compatible with well-known complementary refinement tools such as PYSCENIC.

Aphid BCR4 Structure and Activity Uncover a New Defensin Peptide Superfamily.

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

Efficient compartmentalization in insect bacteriomes protects symbiotic bacteria from host immune system.

BACKGROUND: Many insects house symbiotic intracellular bacteria (endosymbionts) that provide them with essential nutrients, thus promoting the usage of nutrient-poor habitats. Endosymbiont seclusion within host specialized cells, called bacteriocytes, often organized in a dedicated organ, the bacteriome, is crucial in protecting them from host immune defenses while avoiding chronic host immune activation. Previous evidence obtained in the cereal weevil Sitophilus oryzae has shown that bacteriome immunity is activated against invading pathogens, suggesting endosymbionts might be targeted and impacted by immune effectors during an immune challenge. To pinpoint any molecular determinants associated with such challenges, we conducted a dual transcriptomic analysis of S. oryzae's bacteriome subjected to immunogenic peptidoglycan fragments. RESULTS: We show that upon immune challenge, the bacteriome actively participates in the innate immune response via induction of antimicrobial peptides (AMPs). Surprisingly, endosymbionts do not undergo any transcriptomic changes, indicating that this potential threat goes unnoticed. Immunohistochemistry showed that TCT-induced AMPs are located outside the bacteriome, excluding direct contact with the endosymbionts. CONCLUSIONS: This work demonstrates that endosymbiont protection during an immune challenge is mainly achieved by efficient confinement within bacteriomes, which provides physical separation between host systemic response and endosymbionts. Video Abstract.

The Di-Symbiotic Systems in the Aphids Sipha maydis and Periphyllus lyropictus Provide a Contrasting Picture of Recent Co-Obligate Nutritional Endosymbiosis in Aphids.

Dependence on multiple nutritional bacterial symbionts forming a metabolic unit has repeatedly evolved in many insect species that feed on nutritionally unbalanced diets such as plant sap. This is the case for aphids of the subfamilies Lachninae and Chaitophorinae, which have evolved di-symbiotic systems in which the ancient obligate nutritional symbiont Buchnera aphidicola is metabolically complemented by an additional nutritional symbiont acquired more recently. Deciphering how different symbionts integrate both metabolically and anatomically in such systems is crucial to understanding how complex nutritional symbiotic systems function and evolve. In this study, we sequenced and analyzed the genomes of the symbionts B. aphidicola and Serratia symbiotica associated with the Chaitophorinae aphids Sipha maydis and Periphyllus lyropictus. Our results show that, in these two species, B. aphidicola and S. symbiotica complement each other metabolically (and their hosts) for the biosynthesis of essential amino acids and vitamins, but with distinct metabolic reactions supported by each symbiont depending on the host species. Furthermore, the S. symbiotica symbiont associated with S. maydis appears to be strictly compartmentalized into the specialized host cells housing symbionts in aphids, the bacteriocytes, whereas the S. symbiotica symbiont associated with P. lyropictus exhibits a highly invasive phenotype, presumably because it is capable of expressing a larger set of virulence factors, including a complete flagellum for bacterial motility. Such contrasting levels of metabolic and anatomical integration for two S. symbiotica symbionts that were recently acquired as nutritional co-obligate partners reflect distinct coevolutionary processes specific to each association.

The transposable element-rich genome of the cereal pest Sitophilus oryzae.

BACKGROUND: The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions. RESULTS: We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis. CONCLUSIONS: Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Isolation of Insect Bacteriocytes as a Platform for Transcriptomic Analyses.

Over the past few decades, various techniques have been developed and optimized for the accurate measurement of RNA abundance in cells or tissues. These methods have been instrumental in gaining insight in complex systems such as host-symbiont associations. The pea aphid model has recently emerged as a powerful and experimentally tractable system for studying symbiotic relationships and it is the subject of a growing number of molecular studies. Nevertheless, the lack of standardized protocols for the collection of bacteriocytes, the specialized host cells harboring the symbionts, has limited its use. This chapter provides a simple, step-by-step dissection protocol for the rapid isolation of aphid bacteriocytes. We then describe an adapted protocol for efficient extraction and purification of bacteriocyte RNA that can be used for most downstream transcriptomic analyses.

Evolutionary novelty in the apoptotic pathway of aphids.

Apoptosis, a conserved form of programmed cell death, shows interspecies differences that may reflect evolutionary diversification and adaptation, a notion that remains largely untested. Among insects, the most speciose animal group, the apoptotic pathway has only been fully characterized in Drosophila melanogaster, and apoptosis-related proteins have been studied in a few other dipteran and lepidopteran species. Here, we studied the apoptotic pathway in the aphid Acyrthosiphon pisum, an insect pest belonging to the Hemiptera, an earlier-diverging and distantly related order. We combined phylogenetic analyses and conserved domain identification to annotate the apoptotic pathway in A. pisum and found low caspase diversity and a large expansion of its inhibitory part, with 28 inhibitors of apoptosis (IAPs). We analyzed the spatiotemporal expression of a selected set of pea aphid IAPs and showed that they are differentially expressed in different life stages and tissues, suggesting functional diversification. Five IAPs are specifically induced in bacteriocytes, the specialized cells housing symbiotic bacteria, during their cell death. We demonstrated the antiapoptotic role of these five IAPs using heterologous expression in a tractable in vivo model, the Drosophila melanogaster developing eye. Interestingly, IAPs with the strongest antiapoptotic potential contain two BIR and two RING domains, a domain association that has not been observed in any other species. We finally analyzed all available aphid genomes and found that they all show large IAP expansion, with new combinations of protein domains, suggestive of evolutionarily novel aphid-specific functions.

Correction to: The genome sequence of the grape phylloxera provides insights into the evolution, adaptation, and invasion routes of an iconic pest.

An amendment to this paper has been published and can be accessed via the original article.

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.

Spatial and morphological reorganization of endosymbiosis during metamorphosis accommodates adult metabolic requirements in a weevil.

Bacterial intracellular symbiosis (endosymbiosis) is widespread in nature and impacts many biological processes. In holometabolous symbiotic insects, metamorphosis entails a complete and abrupt internal reorganization that creates a constraint for endosymbiont transmission from larvae to adults. To assess how endosymbiosis copes-and potentially evolves-throughout this major host-tissue reorganization, we used the association between the cereal weevil Sitophilus oryzae and the bacterium Sodalis pierantonius as a model system. S. pierantonius are contained inside specialized host cells, the bacteriocytes, that group into an organ, the bacteriome. Cereal weevils require metabolic inputs from their endosymbiont, particularly during adult cuticle synthesis, when endosymbiont load increases dramatically. By combining dual RNA-sequencing analyses and cell imaging, we show that the larval bacteriome dissociates at the onset of metamorphosis and releases bacteriocytes that undergo endosymbiosis-dependent transcriptomic changes affecting cell motility, cell adhesion, and cytoskeleton organization. Remarkably, bacteriocytes turn into spindle cells and migrate along the midgut epithelium, thereby conveying endosymbionts to midgut sites where future mesenteric caeca will develop. Concomitantly, endosymbiont genes encoding a type III secretion system and a flagellum apparatus are transiently up-regulated while endosymbionts infect putative stem cells and enter their nuclei. Infected cells then turn into new differentiated bacteriocytes and form multiple new bacteriomes in adults. These findings show that endosymbiosis reorganization in a holometabolous insect relies on a synchronized host-symbiont molecular and cellular "choreography" and illustrates an adaptive feature that promotes bacteriome multiplication to match increased metabolic requirements in emerging adults.

Drosophila-associated bacteria differentially shape the nutritional requirements of their host during juvenile growth.

The interplay between nutrition and the microbial communities colonizing the gastrointestinal tract (i.e., gut microbiota) determines juvenile growth trajectory. Nutritional deficiencies trigger developmental delays, and an immature gut microbiota is a hallmark of pathologies related to childhood undernutrition. However, how host-associated bacteria modulate the impact of nutrition on juvenile growth remains elusive. Here, using gnotobiotic Drosophila melanogaster larvae independently associated with Acetobacter pomorumWJL (ApWJL) and Lactobacillus plantarumNC8 (LpNC8), 2 model Drosophila-associated bacteria, we performed a large-scale, systematic nutritional screen based on larval growth in 40 different and precisely controlled nutritional environments. We combined these results with genome-based metabolic network reconstruction to define the biosynthetic capacities of Drosophila germ-free (GF) larvae and its 2 bacterial partners. We first established that ApWJL and LpNC8 differentially fulfill the nutritional requirements of the ex-GF larvae and parsed such difference down to individual amino acids, vitamins, other micronutrients, and trace metals. We found that Drosophila-associated bacteria not only fortify the host's diet with essential nutrients but, in specific instances, functionally compensate for host auxotrophies by either providing a metabolic intermediate or nutrient derivative to the host or by uptaking, concentrating, and delivering contaminant traces of micronutrients. Our systematic work reveals that beyond the molecular dialogue engaged between the host and its bacterial partners, Drosophila and its associated bacteria establish an integrated nutritional network relying on nutrient provision and utilization.

Weevil pgrp-lb prevents endosymbiont TCT dissemination and chronic host systemic immune activation.

Long-term intracellular symbiosis (or endosymbiosis) is widely distributed across invertebrates and is recognized as a major driving force in evolution. However, the maintenance of immune homeostasis in organisms chronically infected with mutualistic bacteria is a challenging task, and little is known about the molecular processes that limit endosymbiont immunogenicity and host inflammation. Here, we investigated peptidoglycan recognition protein (PGRP)-encoding genes in the cereal weevil Sitophilus zeamais's association with Sodalis pierantonius endosymbiont. We discovered that weevil pgrp-lb generates three transcripts via alternative splicing and differential regulation. A secreted isoform is expressed in insect tissues under pathogenic conditions through activation of the PGRP-LC receptor of the immune deficiency pathway. In addition, cytosolic and transmembrane isoforms are permanently produced within endosymbiont-bearing organ, the bacteriome, in a PGRP-LC-independent manner. Bacteriome isoforms specifically cleave the tracheal cytotoxin (TCT), a peptidoglycan monomer released by endosymbionts. pgrp-lb silencing by RNAi results in TCT escape from the bacteriome to other insect tissues, where it chronically activates the host systemic immunity through PGRP-LC. While such immune deregulations did not impact endosymbiont load, they did negatively affect host physiology, as attested by a diminished sexual maturation of adult weevils. Whereas pgrp-lb was first described in pathogenic interactions, this work shows that, in an endosymbiosis context, specific bacteriome isoforms have evolved, allowing endosymbiont TCT scavenging and preventing chronic endosymbiont-induced immune responses, thus promoting host homeostasis.

Identification by Tn-seq of Dickeya dadantii genes required for survival in chicory plants.

The identification of the virulence factors of plant-pathogenic bacteria has relied on the testing of individual mutants on plants, a time-consuming process. Transposon sequencing (Tn-seq) is a very powerful method for the identification of the genes required for bacterial growth in their host. We used this method in a soft-rot pathogenic bacterium to identify the genes required for the multiplication of Dickeya dadantii in chicory. About 100 genes were identified showing decreased or increased fitness in the plant. Most had no previously attributed role in plant-bacterium interactions. Following our screening, in planta competition assays confirmed that the uridine monophosphate biosynthesis pathway and the purine biosynthesis pathway were essential to the survival of D. dadantii in the plant, as the mutants ∆carA, ∆purF, ∆purL, ∆guaB and ∆pyrE were unable to survive in the plant in contrast with the wild-type (WT) bacterium. This study also demonstrated that the biosynthetic pathways of leucine, cysteine and lysine were essential for bacterial survival in the plant and that RsmC and GcpA were important in the regulation of the infection process, as the mutants ∆rsmC and ∆gcpA were hypervirulent. Finally, our study showed that D. dadantii flagellin was glycosylated and that this modification conferred fitness to the bacterium during plant infection. Assay by this method of the large collections of environmental pathogenic strains now available will allow an easy and rapid identification of new virulence factors.

Bacteriocyte Reprogramming to Cope With Nutritional Stress in a Phloem Sap Feeding Hemipteran, the Pea Aphid Acyrthosiphon pisum.

Nutritional symbioses play a central role in the ability of insects to thrive on unbalanced diets and in ensuring their evolutionary success. A genomic model for nutritional symbiosis comprises the hemipteran Acyrthosiphon pisum, and the gamma-3-proteobacterium, Buchnera aphidicola, with genomes encoding highly integrated metabolic pathways. A. pisum feeds exclusively on plant phloem sap, a nutritionally unbalanced diet highly variable in composition, thus raising the question of how this symbiotic system responds to nutritional stress. We addressed this by combining transcriptomic, phenotypic and life history trait analyses to determine the organismal impact of deprivation of tyrosine and phenylalanine. These two aromatic amino acids are essential for aphid development, are synthesized in a metabolic pathway for which the aphid host and the endosymbiont are interdependent, and their concentration can be highly variable in plant phloem sap. We found that this nutritional challenge does not have major phenotypic effects on the pea aphid, except for a limited weight reduction and a 2-day delay in onset of nymph laying. Transcriptomic analyses through aphid development showed a prominent response in bacteriocytes (the core symbiotic tissue which houses the symbionts), but not in gut, thus highlighting the role of bacteriocytes as major modulators of this homeostasis. This response does not involve a direct regulation of tyrosine and phenylalanine biosynthetic pathway and transporter genes. Instead, we observed an extensive transcriptional reprogramming of the bacteriocyte with a rapid down-regulation of genes encoding sugar transporters and genes required for sugar metabolism. Consistently, we observed continued overexpression of the A. pisum homolog of RRAD, a small GTPase implicated in repressing aerobic glycolysis. In addition, we found increased transcription of genes involved in proliferation, cell size control and signaling. We experimentally confirmed the significance of these gene expression changes detecting an increase in bacteriocyte number and cell size in vivo under tyrosine and phenylalanine depletion. Our results support a central role of bacteriocytes in the aphid response to amino acid deprivation: their transcriptional and cellular responses fine-tune host physiology providing the host insect with an effective way to cope with the challenges posed by the variability in composition of phloem sap.

Bacteriocyte cell death in the pea aphid/Buchnera symbiotic system.

Symbiotic associations play a pivotal role in multicellular life by facilitating acquisition of new traits and expanding the ecological capabilities of organisms. In insects that are obligatorily dependent on intracellular bacterial symbionts, novel host cells (bacteriocytes) or organs (bacteriomes) have evolved for harboring beneficial microbial partners. The processes regulating the cellular life cycle of these endosymbiont-bearing cells, such as the cell-death mechanisms controlling their fate and elimination in response to host physiology, are fundamental questions in the biology of symbiosis. Here we report the discovery of a cell-death process involved in the degeneration of bacteriocytes in the hemipteran insect Acyrthosiphon pisum This process is activated progressively throughout aphid adulthood and exhibits morphological features distinct from known cell-death pathways. By combining electron microscopy, immunohistochemistry, and molecular analyses, we demonstrated that the initial event of bacteriocyte cell death is the cytoplasmic accumulation of nonautophagic vacuoles, followed by a sequence of cellular stress responses including the formation of autophagosomes in intervacuolar spaces, activation of reactive oxygen species, and Buchnera endosymbiont degradation by the lysosomal system. We showed that this multistep cell-death process originates from the endoplasmic reticulum, an organelle exhibiting a unique reticular network organization spread throughout the entire cytoplasm and surrounding Buchnera aphidicola endosymbionts. Our findings provide insights into the cellular and molecular processes that coordinate eukaryotic host and endosymbiont homeostasis and death in a symbiotic system and shed light on previously unknown aspects of bacteriocyte biological functioning.