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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Transcriptional profiling of Klebsiella pneumoniae defines signatures for planktonic, sessile and biofilm-dispersed cells.

BACKGROUND: Surface-associated communities of bacteria, known as biofilms, play a critical role in the persistence and dissemination of bacteria in various environments. Biofilm development is a sequential dynamic process from an initial bacterial adhesion to a three-dimensional structure formation, and a subsequent bacterial dispersion. Transitions between these different modes of growth are governed by complex and partially known molecular pathways. RESULTS: Using RNA-seq technology, our work provided an exhaustive overview of the transcriptomic behavior of the opportunistic pathogen Klebsiella pneumoniae derived from free-living, biofilm and biofilm-dispersed states. For each of these conditions, the combined use of Z-scores and principal component analysis provided a clear illustration of distinct expression profiles. In particular, biofilm-dispersed cells appeared as a unique stage in the bacteria lifecycle, different from both planktonic and sessile states. The K-means cluster analysis showed clusters of Coding DNA Sequences (CDS) and non-coding RNA (ncRNA) genes differentially transcribed between conditions. Most of them included dominant functional classes, emphasizing the transcriptional changes occurring in the course of K. pneumoniae lifestyle transitions. Furthermore, analysis of the whole transcriptome allowed the selection of an overall of 40 transcriptional signature genes for the five bacterial physiological states. CONCLUSIONS: This transcriptional study provides additional clues to understand the key molecular mechanisms involved in the transition between biofilm and the free-living lifestyles, which represents an important challenge to control both beneficial and harmful biofilm. Moreover, this exhaustive study identified physiological state specific transcriptomic reference dataset useful for the research community.

Weevil endosymbiont dynamics is associated with a clamping of immunity.

BACKGROUND: Insects subsisting on nutritionally unbalanced diets have evolved long-term mutualistic relationships with intracellular symbiotic bacteria (endosymbionts). The endosymbiont population load undergoes changes along with insect development. In the cereal weevil Sitophilus oryzae, the midgut endosymbionts Sodalis pierantonius drastically multiply following adult metamorphosis and rapidly decline until total elimination when the insect achieves its cuticle synthesis. Whilst symbiont load was shown to timely meet insect metabolic needs, little is known about the host molecular and immune processes underlying this dynamics. METHODS: We performed RNA sequencing analysis on weevil midguts at three representative phases of the endosymbiont dynamics (i.e. increase, climax and decrease). To screen genes which transcriptional changes are specifically related to symbiont dynamics and not to the intrinsic development of the midgut, we further have monitored by RT-qPCR sixteen gene transcript levels in symbiotic and artificially non-symbiotic (aposymbiotic) weevils. We also localized the endosymbionts during the elimination process by fluorescence microscopy. RESULTS: Functional analysis of the host differentially expressed genes by RNA sequencing showed that the main transcriptional changes occur during endosymbiont growth phase and affect cell proliferation, apoptosis, autophagy, phagocytosis, and metabolism of fatty acids and nucleic acids. We also showed that symbiont dynamics alters the expression of several genes involved in insect development. Our results strengthened the implication of apoptosis and autophagy processes in symbiont elimination and recycling. Remarkably, apart from the coleoptericin A that is known to target endosymbionts and controls their division and location, no gene coding antimicrobial peptide was upregulated during the symbiont growth and elimination phases. CONCLUSION: We show that endosymbiont dynamics parallels numerous transcriptional changes in weevil developing adults and affects several biological processes, including metabolism and development. It also triggers cell apoptosis, autophagy and gut epithelial cell swelling and delamination. Strikingly, immunity is repressed during the whole process, presumably avoiding tissue inflammation and allowing insects to optimize nutrient recovery from recycled endosymbiont.

Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine.

BACKGROUND: A seventh order of methanogens, the Methanomassiliicoccales, has been identified in diverse anaerobic environments including the gastrointestinal tracts (GIT) of humans and other animals and may contribute significantly to methane emission and global warming. Methanomassiliicoccales are phylogenetically distant from all other orders of methanogens and belong to a large evolutionary branch composed by lineages of non-methanogenic archaea such as Thermoplasmatales, the Deep Hydrothermal Vent Euryarchaeota-2 (DHVE-2, Aciduliprofundum boonei) and the Marine Group-II (MG-II). To better understand this new order and its relationship to other archaea, we manually curated and extensively compared the genome sequences of three Methanomassiliicoccales representatives derived from human GIT microbiota, "Candidatus Methanomethylophilus alvus", "Candidatus Methanomassiliicoccus intestinalis" and Methanomassiliicoccus luminyensis. RESULTS: Comparative analyses revealed atypical features, such as the scattering of the ribosomal RNA genes in the genome and the absence of eukaryotic-like histone gene otherwise present in most of Euryarchaeota genomes. Previously identified in Thermoplasmatales genomes, these features are presently extended to several completely sequenced genomes of this large evolutionary branch, including MG-II and DHVE2. The three Methanomassiliicoccales genomes share a unique composition of genes involved in energy conservation suggesting an original combination of two main energy conservation processes previously described in other methanogens. They also display substantial differences with each other, such as their codon usage, the nature and origin of their CRISPRs systems and the genes possibly involved in particular environmental adaptations. The genome of M. luminyensis encodes several features to thrive in soil and sediment conditions suggesting its larger environmental distribution than GIT. Conversely, "Ca. M. alvus" and "Ca. M. intestinalis" do not present these features and could be more restricted and specialized on GIT. Prediction of the amber codon usage, either as a termination signal of translation or coding for pyrrolysine revealed contrasted patterns among the three genomes and suggests a different handling of the Pyl-encoding capacity. CONCLUSIONS: This study represents the first insights into the genomic organization and metabolic traits of the seventh order of methanogens. It suggests contrasted evolutionary history among the three analyzed Methanomassiliicoccales representatives and provides information on conserved characteristics among the overall methanogens and among Thermoplasmata.

Large scale explorative oligonucleotide probe selection for thousands of genetic groups on a computing grid: application to phylogenetic probe design using a curated small subunit ribosomal RNA gene database.

Phylogenetic Oligonucleotide Arrays (POAs) were recently adapted for studying the huge microbial communities in a flexible and easy-to-use way. POA coupled with the use of explorative probes to detect the unknown part is now one of the most powerful approaches for a better understanding of microbial community functioning. However, the selection of probes remains a very difficult task. The rapid growth of environmental databases has led to an exponential increase of data to be managed for an efficient design. Consequently, the use of high performance computing facilities is mandatory. In this paper, we present an efficient parallelization method to select known and explorative oligonucleotide probes at large scale using computing grids. We implemented a software that generates and monitors thousands of jobs over the European Computing Grid Infrastructure (EGI). We also developed a new algorithm for the construction of a high-quality curated phylogenetic database to avoid erroneous design due to bad sequence affiliation. We present here the performance and statistics of our method on real biological datasets based on a phylogenetic prokaryotic database at the genus level and a complete design of about 20,000 probes for 2,069 genera of prokaryotes.

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.

KASpOD--a web service for highly specific and explorative oligonucleotide design.

SUMMARY: KASpOD is a web service dedicated to the design of signature sequences using a k-mer-based algorithm. Such highly specific and explorative oligonucleotides are then suitable for various goals, including Phylogenetic Oligonucleotide Arrays. AVAILABILITY: http://g2im.u-clermont1.fr/kaspod. CONTACT: eric.peyretaillade@udamail.fr SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

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.

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.