Genome-Inferred Correspondence between Phylogeny and Metabolic Traits in the Wild Drosophila Gut Microbiome.

Annotated genome sequences provide valuable insight into the functional capabilities of members of microbial communities. Nevertheless, most studies on the microbiome in animal guts use metagenomic data, hampering the assignment of genes to specific microbial taxa. Here, we make use of the readily culturable bacterial communities in the gut of the fruit fly Drosophila melanogaster to obtain draft genome sequences for 96 isolates from wild flies. These include 81 new de novo assembled genomes, assigned to three orders (Enterobacterales, Lactobacillales, and Rhodospirillales) with 80% of strains identified to species level using average nucleotide identity and phylogenomic reconstruction. Based on annotations by the RAST pipeline, among-isolate variation in metabolic function partitioned strongly by bacterial order, particularly by amino acid metabolism (Rhodospirillales), fermentation, and nucleotide metabolism (Lactobacillales) and arginine, urea, and polyamine metabolism (Enterobacterales). Seven bacterial species, comprising 2-3 species in each order, were well-represented among the isolates and included ≥5 strains, permitting analysis of metabolic functions in the accessory genome (i.e., genes not present in every strain). Overall, the metabolic function in the accessory genome partitioned by bacterial order. Two species, Gluconobacter cerinus (Rhodospirillales) and Lactiplantibacillus plantarum (Lactobacillales) had large accessory genomes, and metabolic functions were dominated by amino acid metabolism (G. cerinus) and carbohydrate metabolism (La. plantarum). The patterns of variation in metabolic capabilities at multiple phylogenetic scales provide the basis for future studies of the ecological and evolutionary processes shaping the diversity of microorganisms associated with natural populations of Drosophila.

Predicted Metabolic Function of the Gut Microbiota of Drosophila melanogaster.

An important goal for many nutrition-based microbiome studies is to identify the metabolic function of microbes in complex microbial communities and their impact on host physiology. This research can be confounded by poorly understood effects of community composition and host diet on the metabolic traits of individual taxa. Here, we investigated these multiway interactions by constructing and analyzing metabolic models comprising every combination of five bacterial members of the Drosophila gut microbiome (from single taxa to the five-member community of Acetobacter and Lactobacillus species) under three nutrient regimes. We show that the metabolic function of Drosophila gut bacteria is dynamic, influenced by community composition, and responsive to dietary modulation. Furthermore, we show that ecological interactions such as competition and mutualism identified from the growth patterns of gut bacteria are underlain by a diversity of metabolic interactions, and show that the bacteria tend to compete for amino acids and B vitamins more frequently than for carbon sources. Our results reveal that, in addition to fermentation products such as acetate, intermediates of the tricarboxylic acid (TCA) cycle, including 2-oxoglutarate and succinate, are produced at high flux and cross-fed between bacterial taxa, suggesting important roles for TCA cycle intermediates in modulating Drosophila gut microbe interactions and the potential to influence host traits. These metabolic models provide specific predictions of the patterns of ecological and metabolic interactions among gut bacteria under different nutrient regimes, with potentially important consequences for overall community metabolic function and nutritional interactions with the host.IMPORTANCE Drosophila is an important model for microbiome research partly because of the low complexity of its mostly culturable gut microbiota. Our current understanding of how Drosophila interacts with its gut microbes and how these interactions influence host traits derives almost entirely from empirical studies that focus on individual microbial taxa or classes of metabolites. These studies have failed to capture fully the complexity of metabolic interactions that occur between host and microbe. To overcome this limitation, we reconstructed and analyzed 31 metabolic models for every combination of the five principal bacterial taxa in the gut microbiome of Drosophila This revealed that metabolic interactions between Drosophila gut bacterial taxa are highly dynamic and influenced by cooccurring bacteria and nutrient availability. Our results generate testable hypotheses about among-microbe ecological interactions in the Drosophila gut and the diversity of metabolites available to influence host traits.

Non-Target Effects of dsRNA Molecules in Hemipteran Insects.

Insect pest control by RNA interference (RNAi)-mediated gene expression knockdown can be undermined by many factors, including small sequence differences between double-stranded RNA (dsRNA) and the target gene. It can also be compromised by effects that are independent of the dsRNA sequence on non-target organisms (known as sequence-non-specific effects). This study investigated the species-specificity of RNAi in plant sap-feeding hemipteran pests. We first demonstrated sequence-non-specific suppression of aphid feeding by dsRNA at dietary concentrations ≥0.5 µg µL-1. Then we quantified the expression of NUC (nuclease) genes in insects administered homologous dsRNA (with perfect sequence identity to the target species) or heterologous dsRNA (generated against a related gene of non-identical sequence in a different insect species). For the aphids Acyrthosiphon pisum and Myzus persicae, significantly reduced NUC expression was obtained with the homologous but not heterologous dsRNA at 0.2 µg µL-1, despite high dsNUC sequence identity. Follow-up experiments demonstrated significantly reduced expression of NUC genes in the whitefly Bemisia tabaci and mealybug Planococcus maritimus administered homologous dsNUCs, but not heterologous aphid dsNUCs. Our demonstration of inefficient expression knockdown by heterologous dsRNA in these insects suggests that maximal dsRNA sequence identity is required for RNAi targeting of related pest species, and that heterologous dsRNAs at appropriate concentrations may not be a major risk to non-target sap-feeding hemipterans.

On the evolutionary origins of host-microbe associations.

Many microorganisms with high prevalence in host populations are beneficial to the host and maintained by specialized transmission mechanisms. Although microbial promotion of host fitness and specificity of the associations undoubtedly enhance microbial prevalence, it is an open question whether these symbiotic traits are also a prerequisite for the evolutionary origin of prevalent microbial taxa. To address this issue, we investigate how processes without positive microbial effects on host fitness or host choice can influence the prevalence of certain microbes in a host population. Specifically, we develop a theoretical model to assess the conditions under which particular microbes can become enriched in animal hosts even when they are not providing a specific benefit to a particular host. We find increased prevalence of specific microbes in a host when both show some overlap in their lifecycles, and especially when both share dispersal routes across a patchy habitat distribution. Our results emphasize that host enrichment per se is not a reliable indicator of beneficial host-microbe interactions. The resulting increase in time spent associated with a host may nevertheless give rise to new selection conditions, which can favor microbial adaptations toward a host-associated lifestyle, and, thus, it could be the foundation for subsequent evolution of mutually beneficial coevolved symbioses.

The impact of the gut microbiome on memory and sleep in Drosophila.

The gut microbiome has been proposed to influence diverse behavioral traits of animals, although the experimental evidence is limited and often contradictory. Here, we made use of the tractability of Drosophila melanogaster for both behavioral analyses and microbiome studies to test how elimination of microorganisms affects a number of behavioral traits. Relative to conventional flies (i.e. with unaltered microbiome), microbiologically sterile (axenic) flies displayed a moderate reduction in memory performance in olfactory appetitive conditioning and courtship assays. The microbiological status of the flies had a small or no effect on anxiety-like behavior (centrophobism) or circadian rhythmicity of locomotor activity, but axenic flies tended to sleep for longer and displayed reduced sleep rebound after sleep deprivation. These last two effects were robust for most tests conducted on both wild-type Canton S and w1118 strains, as well for tests using an isogenized panel of flies with mutations in the period gene, which causes altered circadian rhythmicity. Interestingly, the effect of absence of microbiota on a few behavioral features, most notably instantaneous locomotor activity speed, varied among wild-type strains. Taken together, our findings demonstrate that the microbiome can have subtle but significant effects on specific aspects of Drosophila behavior, some of which are dependent on genetic background.

Evaluation of RNA Interference for Control of the Grape Mealybug Pseudococcus maritimus (Hemiptera: Pseudococcidae).

The grape mealybug Pseudococcus maritimus (Ehrhorn, 1900) (Hemiptera: Pseudococcidae) is a significant pest of grapevines (Vitis spp.) and a vector of disease-causing grape viruses, linked to its feeding on phloem sap. The management of this pest is constrained by the lack of naturally occurring resistance traits in Vitis. Here, we obtained proof of concept that RNA interference (RNAi) using double-stranded RNA (dsRNA) molecules against essential genes for phloem sap feeding can depress insect survival. The genes of interest code for an aquaporin (AQP) and a sucrase (SUC) that are required for osmoregulation in related phloem sap-feeding hemipteran insects (aphids and whiteflies). In parallel, we investigated the grape mealybug genes coding non-specific nucleases (NUC), which reduce RNAi efficacy by degrading administered dsRNA. Homologs of AQP and SUC with experimentally validated function in aphids, together with NUC, were identified in the published transcriptome of the citrus mealybug Planococcus citri by phylogenetic analysis, and sequences of the candidate genes were obtained for Ps. maritimus by PCR with degenerate primers. Using this first sequence information for Ps. maritimus, dsRNA was prepared and administered to the insects via an artificial diet. The treatment comprising dsRNA against AQP, SUC and NUC significantly increased insect mortality over three days, relative to dsRNA-free controls. The dsRNA constructs for AQP and NUC were predicted, from sequence analysis to have some activity against other mealybugs, but none of the three dsRNA constructs have predicted activity against aphids. This study provides the basis to develop in planta RNAi strategies against Ps. maritimus and other mealybug pests of grapevines.

How gut microbiome interactions affect nutritional traits of Drosophila melanogaster.

Most research on the impact of the gut microbiome on animal nutrition is designed to identify the effects of single microbial taxa and single metabolites of microbial origin, without considering the potentially complex network of interactions among co-occurring microorganisms. Here, we investigated how different microbial associations and their fermentation products affect host nutrition, using Drosophila melanogaster colonized with three gut microorganisms (the bacteria Acetobacter fabarum and Lactobacillus brevis, and the yeast Hanseniaspora uvarum) in all seven possible combinations. Some microbial effects on host traits could be attributed to single taxa (e.g. yeast-mediated reduction of insect development time), while other effects were sex specific and driven by among-microbe interactions (e.g. male lipid content determined by interactions between the yeast and both bacteria). Parallel analysis of nutritional indices of microbe-free flies administered different microbial fermentation products (acetic acid, acetoin, ethanol and lactic acid) revealed a single consistent effect: that the lipid content of both male and female flies is reduced by acetic acid. This effect was recapitulated in male flies colonized with both yeast and A. fabarum, but not for any microbial treatment in females or males with other microbial complements. These data suggest that the effect of microbial fermentation products on host nutritional status is strongly context dependent, with respect to both the combination of associated microorganisms and host sex. Taken together, our findings demonstrate that among-microbe interactions can play a critically important role in determining the physiological outcome of host-microbiome interactions in Drosophila and, likely, in other animal hosts.

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.

Housing microbial symbionts: evolutionary origins and diversification of symbiotic organs in animals.

In many animal hosts, microbial symbionts are housed within specialized structures known as symbiotic organs, but the evolutionary origins of these structures have rarely been investigated. Here, I adopt an evolutionary developmental (evo-devo) approach, specifically to apply knowledge of the development of symbiotic organs to gain insights into their evolutionary origins and diversification. In particular, host genetic changes associated with evolution of symbiotic organs can be inferred from studies to identify the host genes that orchestrate the development of symbiotic organs, recognizing that microbial products may also play a key role in triggering the developmental programme in some associations. These studies may also reveal whether higher animal taxonomic groups (order, class, phylum, etc.) possess a common genetic regulatory network for symbiosis that is latent in taxa lacking symbiotic organs, and activated at the origination of symbiosis in different host lineages. In this way, apparent instances of convergent evolution of symbiotic organs may be homologous in terms of a common genetic blueprint for symbiosis. Advances in genetic technologies, including reverse genetic tools and genome editing, will facilitate the application of evo-devo approaches to investigate the evolution of symbiotic organs in animals. This article is part of the theme issue 'The role of the microbiome in host evolution'.

The role of the microbiome in host evolution.

In the last decade, we have witnessed a major paradigm shift in the life sciences: the recognition that the microbiome, i.e. the set of microorganisms associated with healthy animals (including humans) and plants, plays a crucial role in the sustained health and fitness of its host. Enabled by rapid advances in sequencing technologies and analytical methods, substantial advances have been achieved in both identifying the microbial taxa and understanding the relationship between microbiome composition and host phenotype. These breakthroughs are leading to novel strategies for improved human and animal health, enhanced crop yield and nutritional quality, and the control of various pests and disease agents. This article is part of the theme issue 'The role of the microbiome in host evolution'.

B-vitamin nutrition in the pea aphid-Buchnera symbiosis.

Various insects that utilize vitamin-deficient diets derive a supplementary supply of these micronutrients from their symbiotic microorganisms. Here, we tested the inference from genome annotation that the symbiotic bacterium Buchnera aphidicola in the pea aphid Acyrthosiphon pisum provides the insect with vitamins B2 and B5 but no other B-vitamins. Contrary to expectation, aphid survival over five days of larval development on artificial diets individually lacking each B-vitamin not synthesized by Buchnera was not significantly reduced, despite significantly lower carcass B1, B3, B6 and B7 concentrations in the aphids on diets lacking each of these B-vitamins than on the vitamin-complete diet. Aphid survival was, however, significantly reduced on diet containing low concentrations (≤0.2 mM) or no pantothenate (B5). Complementary transcriptome analysis revealed low abundance of the sense-transcript, but high abundance of the antisense transcript, of the Buchnera gene panC encoding the enzyme mediating the terminal reaction in pantothenate synthesis. We hypothesize that metabolic constraints or antisense transcripts may reduce Buchnera-mediated production of pantothenate, resulting in poor aphid performance on pantothenate-free diets. The discrepancy between predictions from genome data and empirical data illustrates the need for physiological study to test functional inferences made from genome annotations.

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.

Impact of Facultative Bacteria on the Metabolic Function of an Obligate Insect-Bacterial Symbiosis.

Beneficial microorganisms associated with animals derive their nutritional requirements entirely from the animal host, but the impact of these microorganisms on host metabolism is largely unknown. The focus of this study was the experimentally tractable tripartite symbiosis between the pea aphid Acyrthosiphon pisum, its obligate intracellular bacterial symbiont Buchnera, and the facultative bacterium Hamiltonella which is localized primarily to the aphid hemolymph (blood). Metabolome experiments on, first, multiple aphid genotypes that naturally bear or lack Hamiltonella and, second, one aphid genotype from which Hamiltonella was experimentally eliminated revealed no significant effects of Hamiltonella on aphid metabolite profiles, indicating that Hamiltonella does not cause major reconfiguration of host metabolism. However, the titer of just one metabolite, 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR), displayed near-significant enrichment in Hamiltonella-positive aphids in both metabolome experiments. AICAR is a by-product of biosynthesis of the essential amino acid histidine in Buchnera and, hence, an index of histidine biosynthetic rates, suggesting that Buchnera-mediated histidine production is elevated in Hamiltonella-bearing aphids. Consistent with this prediction, aphids fed on [13C]histidine yielded a significantly elevated 12C/13C ratio of histidine in Hamiltonella-bearing aphids, indicative of increased (∼25%) histidine synthesized de novo by Buchnera However, in silico analysis predicted an increase of only 0.8% in Buchnera histidine synthesis in Hamiltonella-bearing aphids. We hypothesize that Hamiltonella imposes increased host demand for histidine, possibly for heightened immune-related functions. These results demonstrate that facultative bacteria can alter the dynamics of host metabolic interactions with co-occurring microorganisms, even when the overall metabolic homeostasis of the host is not substantially perturbed.IMPORTANCE Although microbial colonization of the internal tissues of animals generally causes septicemia and death, various animals are persistently associated with benign or beneficial microorganisms in their blood or internal organs. The metabolic consequences of these persistent associations for the animal host are largely unknown. Our research on the facultative bacterium Hamiltonella, localized primarily to the hemolymph of pea aphids, demonstrated that although Hamiltonella imposed no major reconfiguration of the aphid metabolome, it did alter the metabolic relations between the aphid and its obligate intracellular symbiont, Buchnera Specifically, Buchnera produced more histidine in Hamiltonella-positive aphids to support both Hamiltonella demand for histidine and Hamiltonella-induced increase in host demand. This study demonstrates how microorganisms associated with internal tissues of animals can influence specific aspects of metabolic interactions between the animal host and co-occurring microorganisms.

Candidate genetic determinants of intraspecific variation in pea aphid susceptibility to RNA interference.

RNA interference (RNAi) plays a key role in insect defense against viruses and transposable elements, and it is being applied as an experimental tool and for insect pest control. However, RNAi efficiency is highly variable for some insects, notably the pea aphid Acyrthosiphon pisum. In this study, we used natural variation in RNAi susceptibility of pea aphids to identify genes that influence RNAi efficiency. Susceptibility to orally-delivered dsRNA against the gut aquaporin gene AQP1 (ds-AQP1) varied widely across a panel of 83 pea aphid genotypes, from zero to total mortality. Genome-wide association between aphid performance on ds-AQP1 supplemented diet and aphid genetic variants yielded 103 significantly associated single nucleotide polymorphisms (SNPs), including variants in 55 genes, at the 10-4 probability cut-off. When ds-AQP1 was co-administered with dsRNA against six candidate genes, aphid mortality was reduced for three (50%) genes: the orthologs of the Drosophila genes trachealess (CG42865), headcase (CG15532) and a gene coding a peritrophin-A domain (CG8192), indicating that these genes function to promote RNAi efficiency against AQP1 in the pea aphid. Aphid susceptibility (quantified as mortality) to ds-AQP1 was correlated with RNAi against a further gene, snakeskin with essential gut function unrelated to AQP1, for some but not all aphid genotypes tested, suggesting that the determinants of RNAi efficiency may be partly gene-specific. This study demonstrates high levels of natural variation in susceptibility to RNAi and demonstrates the value of harnessing this variation to identify genes influencing RNAi efficiency.

Syntrophic splitting of central carbon metabolism in host cells bearing functionally different symbiotic bacteria.

Insects feeding on the nutrient-poor diet of xylem plant sap generally bear two microbial symbionts that are localized to different organs (bacteriomes) and provide complementary sets of essential amino acids (EAAs). Here, we investigate the metabolic basis for the apparent paradox that xylem-feeding insects are under intense selection for metabolic efficiency but incur the cost of maintaining two symbionts for functions mediated by one symbiont in other associations. Using stable isotope analysis of central carbon metabolism and metabolic modeling, we provide evidence that the bacteriomes of the spittlebug Clastoptera proteus display high rates of aerobic glycolysis, with syntrophic splitting of glucose oxidation. Specifically, our data suggest that one bacteriome (containing the bacterium Sulcia, which synthesizes seven EAAs) predominantly processes glucose glycolytically, producing pyruvate and lactate, and the exported pyruvate and lactate is assimilated by the second bacteriome (containing the bacterium Zinderia, which synthesizes three energetically costly EAAs) and channeled through the TCA cycle for energy generation by oxidative phosphorylation. We, furthermore, calculate that this metabolic arrangement supports the high ATP demand in Zinderia bacteriomes for Zinderia-mediated synthesis of energy-intensive EAAs. We predict that metabolite cross-feeding among host cells may be widespread in animal-microbe symbioses utilizing low-nutrient diets.

The microbial exometabolome: ecological resource and architect of microbial communities.

All microorganisms release many metabolites, collectively known as the exometabolome. The resultant multi-way cross-feeding of metabolites among microorganisms distributes resources, thereby increasing total biomass of the microbial community, and promotes the recruitment and persistence of phylogenetically and functionally diverse taxa in microbial communities. Metabolite transfer can also select for evolutionary diversification, yielding multiple closely related but functionally distinct strains. Depending on starting conditions, the evolved strains may be auxotrophs requiring metabolic outputs from producer cells or, alternatively, display loss of complementary reactions in metabolic pathways, with increased metabolic efficiency. Metabolite cross-feeding is widespread in many microbial communities associated with animals and plants, including the animal gut microbiome, and these metabolic interactions can yield products valuable to the host. However, metabolite exchange between pairs of intracellular microbial taxa that share the same host cell or organ can be very limited compared to pairs of free-living microorganisms, perhaps as a consequence of host controls over the metabolic function of intracellular microorganisms. Priorities for future research include the development of tools for improved quantification of metabolite exchange in complex communities and greater integration of the roles of metabolic cross-feeding and other ecological processes, including priority effects and antagonistic interactions, in shaping microbial communities. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Functional traits of the gut microbiome correlated with host lipid content in a natural population of Drosophila melanogaster.

Most research on the nutritional significance of the gut microbiome is conducted on laboratory animals, and its relevance to wild animals is largely unknown. This study investigated the microbiome correlates of lipid content in individual wild fruit flies, Drosophila melanogaster. Lipid content varied 3.6-fold among the flies and was significantly correlated with the abundance of gut-derived bacterial DNA sequences that were assigned to genes contributing to 16 KEGG pathways. These included genes encoding sugar transporters and enzymes in glycolysis/gluconeogenesis, potentially promoting sugar consumption by the gut microbiome and, thereby, a lean fly phenotype. Furthermore, the lipid content of wild flies was significantly lower than laboratory flies, indicating that, as for some mammalian models, certain laboratory protocols might be obesogenic for Drosophila. This study demonstrates the value of research on natural populations to identify candidate microbial genes that influence ecologically important host traits.

Host and symbiont genetic determinants of nutritional phenotype in a natural population of the pea aphid.

A defining feature of the nutritional ecology of plant sap-feeding insects is that the dietary deficit of essential amino acids (EAAs) in plant sap is supplemented by EAA-provisioning microbial symbionts in the insect. Here, we demonstrated substantial variation in the nutritional phenotype of 208 genotypes of the pea aphid Acyrthosiphon pisum collected from a natural population. Specifically, the genotypes varied in performance (larval growth rates) on four test diets lacking the EAAs arginine, histidine and methionine or aromatic EAAs (phenylalanine and tryptophan), relative to the diet containing all EAAs. These data indicate that EAA supply from the symbiotic bacteria Buchnera can meet total aphid nutritional demand for only a subset of the EAA/aphid genotype combinations. We then correlated single nucleotide polymorphisms (SNPs) identified in the aphid and Buchnera genomes by reduced genome sequencing against aphid performance for each EAA deletion diet. This yielded significant associations between performance on the histidine-free diet and Buchnera SNPs, including metabolism genes predicted to influence histidine biosynthesis. Aphid genetic correlates of performance were obtained for all four deletion diets, with associations on the arginine-free diet and aromatic-free diets dominated by genes functioning in the regulation of metabolic and cellular processes. The specific aphid genes associated with performance on different EAA deletion diets are largely nonoverlapping, indicating some independence in the regulatory circuits determining aphid phenotype for the different EAAs. This study demonstrates how variation in the phenotype of associations collected from natural populations can be applied to elucidate the genetic basis of ecologically important traits in systems intractable to traditional forward/reverse genetic techniques.

Patterns of host cell inheritance in the bacterial symbiosis of whiteflies.

Whiteflies possess bacterial symbionts Candidatus Portiera aleyrodidium that are housed in specialized cells called bacteriocytes and are faithfully transmitted via the ovary to insect offspring. In one whitefly species studied previously, Bemisia tabaci MEAM1, transmission is mediated by somatic inheritance of bacteriocytes, with a single bacteriocyte transferred to each oocyte and persisting through embryogenesis to the next generation. Here, we investigate the mode of bacteriocyte transmission in two whitefly species, B. tabaci MED, the sister species of MEAM1, and the phylogenetically distant species Trialeurodes vaporariorum. Microsatellite analysis supported by microscopical studies demonstrates that B. tabaci MED bacteriocytes are genetically different from other somatic cells and persist through embryogenesis, as for MEAM1, but T. vaporariorum bacteriocytes are genetically identical to other somatic cells of the insect, likely mediated by the degradation of maternal bacteriocytes in the embryo. These two alternative modes of transmission provide a first demonstration among insect symbioses that the cellular processes underlying vertical transmission of bacterial symbionts can diversify among related host species associated with a single lineage of symbiotic bacteria.