Genome reduction and co-evolution between the primary and secondary bacterial symbionts of psyllids

Complete genome of the mutualistic symbiont "Candidatus Carsonella ruddii" from a Japanese island strain of the Asian citrus psyllid Diaphorina citri

The complete genome, 173,958 bp in size, of "Candidatus Carsonella ruddii" DC-OKEB1, an obligate bacterial endosymbiont of the Asian citrus psyllid Diaphorina citri, was determined. The genome sequence provides valuable information for comparative and evolutionary aspects of the intimate insect-microbe mutualism.

A dual insect symbiont and plant pathogen improves insect host fitness under arginine limitation

Some facultative bacterial symbionts are known to benefit insects, but nutritional advantages are rare among these non-obligate symbionts. Here, we demonstrate that the facultative symbiont Candidatus Liberibacter psyllaurous enhances the fitness of its psyllid insect host, Bactericera cockerelli, by providing nutritional benefits. L. psyllaurous, an unculturable pathogen of solanaceous crops, also establishes a close relationship with its insect vector, B. cockerelli, increasing in titer during insect development, vertically transmitting through eggs, and colonizing various tissues, including the bacteriome, which houses the obligate nutritional symbiont, Carsonella. Carsonella supplies essential amino acids to its insect host but has gaps in some of its essential amino acid pathways that the psyllid complements with its own genes, many of which have been acquired through horizontal gene transfer (HGT) from bacteria. Our findings reveal that L. psyllaurous increases psyllid fitness on plants by reducing developmental time and increasing adult weight. In addition, through metagenomic sequencing, we reveal that L. psyllaurous maintains complete pathways for synthesizing the essential amino acids arginine, lysine, and threonine, unlike the psyllid's other resident microbiota, Carsonella, and two co-occurring Wolbachia strains. RNA sequencing reveals the downregulation of a HGT collaborative psyllid gene (ASL), which indicates a reduced demand for arginine supplied by Carsonella when the psyllid is infected with L. psyllaurous. Notably, artificial diet assays show that L. psyllaurous enhances psyllid fitness on an arginine-deplete diet. These results corroborate the role of L. psyllaurous as a beneficial insect symbiont, contributing to the nutrition of its insect host.IMPORTANCEUnlike obligate symbionts that are permanently associated with their hosts, facultative symbionts rarely show direct nutritional contributions, especially under nutrient-limited conditions. This study demonstrates, for the first time, that Candidatus Liberibacter psyllaurous, a facultative symbiont and a plant pathogen, enhances the fitness of its Bactericera cockerelli host by supplying an essential nutrient arginine that is lacking in the plant sap diet. Our findings reveal how facultative symbionts can play a vital role in helping their insect hosts adapt to nutrient-limited environments. This work provides new insights into the dynamic interactions between insect hosts, their symbiotic microbes, and their shared ecological niches, broadening our understanding of symbiosis and its role in shaping adaptation and survival.

Origin and function of beneficial bacterial symbioses in insects

Beneficial bacterial symbionts are widespread in insects and affect the fitness of their hosts by contributing to nutrition, digestion, detoxification, communication or protection from abiotic stressors or natural enemies. Decades of research have formed our understanding of the identity, localization and functional benefits of insect symbionts, and the increasing availability of genome sequences spanning a diversity of pathogens and beneficial bacteria now enables comparative approaches of their metabolic features and their phylogenetic affiliations, shedding new light on the origin and function of beneficial symbioses in insects. In this Review, we explore the symbionts' metabolic traits that can provide benefits to insect hosts and discuss the evolutionary paths to the formation of host-beneficial symbiotic associations. Phylogenetic analyses and molecular studies reveal that extracellular symbioses colonizing cuticular organs or the digestive tract evolved from a broad diversity of bacterial partners, whereas intracellular beneficial symbionts appear to be restricted to a limited number of lineages within the Gram-negative bacteria and probably originated from parasitic ancestors. To unravel the general principles underlying host-symbiont interactions and recapitulate the early evolutionary steps leading towards beneficial symbioses, future efforts should aim to establish more symbiotic systems that are amenable to genetic manipulation and experimental evolution.

Metrics of Genomic Complexity in the Evolution of Bacterial Endosymbiosis

Endosymbiosis can be considered a regressive or degenerative evolutionary process characterized at the genomic level by genome erosion and degeneration due to high mutational pressure toward AT (adenine and thymine) bases. The genomic and biological complexity of endosymbionts must be lower than that of the free-living bacteria from which they evolved. In the present work, we contrasted whether two proposed metrics for measuring genomic complexity in both types of bacteria, GS and BB, reflect their complexity, expecting higher values in free-living bacteria than in endosymbionts. On the other hand, we endeavored to delve into the factors that contribute to the reduction in metric values in endosymbionts, as well as their eventual relationship with six genomic parameters associated with functionality. This study aimed to test the robustness of these proposed metrics in a well-known biological scenario, such as the endosymbiosis process.

Genomes of the Bacterial Endosymbionts of Carrot Psyllid Trioza apicalis Suggest Complementary Biosynthetic Capabilities

Carrot psyllid Trioza apicalis is a serious pest of cultivated carrot and also a vector of the plant pathogen 'Candidatus Liberibacter solanacearum' (Lso). To find out whether T. apicalis harbours other species of bacteria that might affect the Lso infection rate, the bacterial communities and metagenome in T. apicalis were studied. Lso haplotype C was detected in a third of the psyllids sampled, at different relative amounts. Surprisingly, T. apicalis was found to harbour only one secondary endosymbiont, a previously unknown species of gamma proteobacterium endosymbiont (Gpe), beside the primary endosymbiont 'Candidatus Carsonella ruddii' (CCr). The relative abundancies of these two endosymbionts were approximately equal. The genomes of CCr, Gpe and Lso were assembled from a T. apicalis metagenome sample. Based on the 16S rRNA gene, the closest relative of Gpe of T. apicalis could be a secondary endosymbiont of Trioza magnoliae. The 253.171 kb Gpe genome contains all the tRNA and rRNA genes and most of the protein-coding genes required for DNA replication, transcription and translation, but it lacks most of the genes for amino acid biosynthesis. Gpe has no genes encoding cell wall peptidoglycan synthesis, suggesting it has no cell wall, and could thus live as an intracellular endosymbiont. Like the CCr of other psyllids, CCr of T. apicalis retains a broad amino acid biosynthetic capacity, whilst lacking many genes required for DNA replication and repair and for transcription and translation. These findings suggest that these two endosymbionts of T. apicalis are complementary in their biosynthetic capabilities.

Detection of Arsenophonus in Glycaspis brimblecombei (Hemiptera: Aphalaridae) populations in Brazil

Eucalyptus is the most intensively managed tree genus in the world. Different factors, including damage by insect pests, affect its growth and productivity. Among these pests is Glycaspis brimblecombei Moore (Hemiptera: Aphalaridae), an exotic insect of Australian origin. The evolutionary success of this insect depends on symbiotic associations with microorganisms. The influence of these microorganisms on insect pests and their natural enemies is important for integrated management tactics. Within this context, this work aimed to detect Arsenophonus in populations of G. brimblecombei in Brazil. Eucalyptus branches infested with G. brimblecombei nymphs were collected in commercial eucalyptus plantations in six Brazilian states. Specimens of this pest were sampled soon after emergence and frozen for molecular analysis. The genomic DNA of G. brimblecombei adults from each population was extracted and used to detect the endosymbiont Arsenophonus by polymerase chain reaction (PCR) employing specific primers that target its 23 S rRNA gene. This endosymbiont was identified in all of the studied G. brimblecombei populations. This is the first report on the association between Arsenophonus and G. brimblecombei in Brazil.

Shrinking in the dark: Parallel endosymbiont genome erosions are associated with repeated host transitions to an underground life

Microbial symbioses have had profound impacts on the evolution of animals. Conversely, changes in host biology may impact the evolutionary trajectory of symbionts themselves. Blattabacterium cuenoti is present in almost all cockroach species and enables hosts to subsist on a nutrient-poor diet. To investigate if host biology has impacted Blattabacterium at the genomic level, we sequenced and analyzed 25 genomes from Australian soil-burrowing cockroaches (Blaberidae: Panesthiinae), which have undergone at least seven separate subterranean, subsocial transitions from above-ground, wood-feeding ancestors. We find at least three independent instances of genome erosion have occurred in Blattabacterium strains exclusive to Australian soil-burrowing cockroaches. These shrinkages have involved the repeated inactivation of genes involved in amino acid biosynthesis and nitrogen recycling, the core role of Blattabacterium in the host-symbiont relationship. The most drastic of these erosions have occurred in hosts thought to have transitioned underground the earliest relative to other lineages, further suggestive of a link between gene loss in Blattabacterium and the burrowing behavior of hosts. As Blattabacterium is unable to fulfill its core function in certain host lineages, these findings suggest soil-burrowing cockroaches must acquire these nutrients from novel sources. Our study represents one of the first cases, to our knowledge, of parallel host adaptations leading to concomitant parallelism in their mutualistic symbionts, further underscoring the intimate relationship between these two partners.

Chromosome Structural Rearrangements in Invasive Haplodiploid Ambrosia Beetles Revealed by the Genomes of Euwallacea fornicatus (Eichhoff) and Euwallacea similis (Ferrari) (Coleoptera, Curculionidae, Scolytinae)

Bark and ambrosia beetles are among the most ecologically and economically damaging introduced plant pests worldwide. Life history traits including polyphagy, haplodiploidy, inbreeding polygyny, and symbiosis with fungi contribute to their dispersal and impact. Species vary in their interactions with host trees, with many attacking stressed or recently dead trees, such as the globally distributed Euwallacea similis (Ferrari). Other species, like the Polyphagous Shot Hole Borer Euwallacea fornicatus (Eichhoff), can attack over 680 host plants and is causing considerable economic damage in several countries. Despite their notoriety, publicly accessible genomic resources for Euwallacea Hopkins species are scarce, hampering our understanding of their invasive capabilities as well as modern control measures, surveillance, and management. Using a combination of long and short read sequencing platforms, we assembled and annotated high quality (BUSCO > 98% complete) pseudo-chromosome-level genomes for these species. Comparative macrosynteny analysis identified an increased number of pseudo-chromosome scaffolds in the haplodiploid inbreeding species of Euwallacea compared to diploid outbred species, due to fission events. This suggests that life history traits can impact chromosome structure. Further, the genome of E. fornicatus had a higher relative proportion of repetitive elements, up to 17% more, than E. similis. Metagenomic assembly pipelines identified microbiota associated with both species including Fusarium fungal symbionts and a novel Wolbachia strain. These novel genomes of haplodiploid inbreeding species will contribute to the understanding of how life history traits are related to their evolution and to the management of these invasive pests.

Rapid Loss of Nutritional Symbionts in an Endemic Hawaiian Herbivore Radiation Is Associated with Plant Galling Habit

Insect herbivores frequently cospeciate with symbionts that enable them to survive on nutritionally unbalanced diets. While ancient symbiont gain and loss events have been pivotal for insect diversification and feeding niche specialization, evidence of recent events is scarce. We examine the recent loss of nutritional symbionts (in as little as 1 MY) in sap-feeding Pariaconus, an endemic Hawaiian insect genus that has undergone adaptive radiation, evolving various galling and free-living ecologies on a single host-plant species, Metrosideros polymorpha within the last ∼5 MY. Using 16S rRNA sequencing, we investigated the bacterial microbiomes of 19 Pariaconus species and identified distinct symbiont profiles associated with specific host-plant ecologies. Phylogenetic analyses and metagenomic reconstructions revealed significant differences in microbial diversity and functions among psyllids with different host-plant ecologies. Within a few millions of years, Pariaconus species convergently evolved the closed-gall habit twice. This shift to enclosed galls coincided with the loss of the Morganella-like symbiont that provides the essential amino acid arginine to free-living and open-gall sister species. After the Pariaconus lineage left Kauai and colonized younger islands, both open- and closed-gall species lost the Dickeya-like symbiont. This symbiont is crucial for synthesizing essential amino acids (phenylalanine, tyrosine, and lysine) as well as B vitamins in free-living species. The recurrent loss of these symbionts in galling species reinforces evidence that galls are nutrient sinks and, combined with the rapidity of the evolutionary timeline, highlights the dynamic role of insect-symbiont relationships during the diversification of feeding ecologies. We propose new Candidatus names for the novel Morganella-like and Dickeya-like symbionts.

Creating and leveraging bespoke large-scale knowledge graphs for comparative genomics and multi-omics drug discovery with SocialGene

The rapid expansion of multi-omics data has transformed biological research, offering unprecedented opportunities to explore complex genomic relationships across diverse organisms. However, the vast volume and heterogeneity of these datasets presents significant challenges for analyses. Here we introduce SocialGene, a comprehensive software suite designed to collect, analyze, and organize multi-omics data into structured knowledge graphs, with the ability to handle small projects to repository-scale analyses. Originally developed to enhance genome mining for natural product drug discovery, SocialGene has been effective across various applications, including functional genomics, evolutionary studies, and systems biology. SocialGene's concerted Python and Nextflow libraries streamline data ingestion, manipulation, aggregation, and analysis, culminating in a custom Neo4j database. The software not only facilitates the exploration of genomic synteny but also provides a foundational knowledge graph supporting the integration of additional diverse datasets and the development of advanced search engines and analyses. This manuscript introduces some of SocialGene's capabilities through brief case studies including targeted genome mining for drug discovery, accelerated searches for similar and distantly related biosynthetic gene clusters in biobank-available organisms, integration of chemical and analytical data, and more. SocialGene is free, open-source, MIT-licensed, designed for adaptability and extension, and available from github.com/socialgene.

Twenty-nine newly sequenced genomes and a comprehensive genome dataset for the insect endosymbiont Buchnera

Most phloem-feeding insects face nutritional deficiency and rely on their intracellular symbionts to provide nutrients, and most of endosymbiont genomes have undergone reduction. However, the study of genome reduction processes of endosymbionts has been constrained by the limited availability of genome data from different insect lineages. The obligate relationship between aphids and Buchnera aphidicola (hereafter Buchnera) makes them a classic model for studying insect-endosymbiont interaction. Here, we report 29 newly sequenced Buchnera genomes from 11 aphid subfamilies, and a comprehensive dataset based on 90 Buchnera genomes from 14 aphid subfamilies. The dataset shows a significant genomic difference of Buchnera among different aphid lineages. The dataset exhibits a more balanced distribution of Buchnera (from 14 aphid subfamilies) genome sizes, ranging from 400 kb to 600 kb, which can illustrate the genome reduction process of Buchnera. The new genome data provide valuable insights into the microevolutionary processes leading to genomic reduction of insect endosymbionts.

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.

Host genetics and larval host plant modulate microbiome structure and evolution underlying the intimate insect-microbe-plant interactions in Parnassius species on the Qinghai-Tibet Plateau

Insects harbor a remarkable diversity of gut microbiomes critical for host survival, health, and fitness, but the mechanism of this structured symbiotic community remains poorly known, especially for the insect group consisting of many closely related species that inhabit the Qinghai-Tibet Plateau. Here, we firstly analyzed population-level 16S rRNA microbial dataset, comprising 11 Parnassius species covering 5 subgenera, from 14 populations mostly sampled in mountainous regions across northwestern-to-southeastern China, and meanwhile clarified the relative importance of multiple factors on gut microbial community structure and evolution. Our findings indicated that both host genetics and larval host plant modulated gut microbial diversity and community structure. Moreover, the effect analysis of host genetics and larval diet on gut microbiomes showed that host genetics played a critical role in governing the gut microbial beta diversity and the symbiotic community structure, while larval host plant remarkably influenced the functional evolution of gut microbiomes. These findings of the intimate insect-microbe-plant interactions jointly provide some new insights into the correlation among the host genetic background, larval host plant, the structure and evolution of gut microbiome, as well as the mechanisms of high-altitude adaptation in closely related species of this alpine butterfly group.

The Hidden Secrets of Psylloidea: Biology, Behavior, Symbionts, and Ecology

Psyllids constitute a diverse group of sap-feeding Sternorrhyncha that were relatively obscure until it was discovered that a handful of species transmit bacterial plant pathogens. Yet the superfamily Psylloidea is much richer than the sum of its crop-associated vectors, with over 4,000 described species exhibiting diverse life histories and host exploitation strategies. A growing body of research is uncovering fascinating insights into psyllid evolution, biology, behavior, and species interactions. This work has revealed commonalities and differences with better-studied Sternorrhyncha, as well as unique evolutionary patterns of lineage divergence and host use. We are also learning how psyllid evolution and foraging ecology underlie life history traits and the roles of psyllids in communities. At finer scales, we are untangling the web of symbionts across the psyllid family tree, linking symbiont and psyllid lineages, and revealing mechanisms underlying reciprocal exchange between symbiont and host. In this review, we synthesize and summarize key advances within these areas with a focus on free-living (nongalling) Psylloidea.

Ultrastructure of the bacteriome and bacterial symbionts in the Asian citrus psyllid, Diaphorina citri

IMPORTANCE

Omics analyses suggested a mutually indispensable tripartite association among the host D. citri and organelle-like bacteriome associates, Carsonella and Profftella, which are vertically transmitted through host generations. This relationship is based on the metabolic complementarity among these organisms, which is partly enabled by horizontal gene transfer between partners. However, little was known about the fine morphology of the symbionts and the bacteriome, the interface among these organisms. As a first step to address this issue, the present study performed transmission electron microscopy, which revealed previously unrecognized ultrastructures, including aggregations of ribosomes in Carsonella, numerous tubes and occasional protrusions of Profftella, apparently degrading Profftella, and host organelles with different abundance and morphology in distinct cell types. These findings provide insights into the behaviors of the symbionts and host cells to maintain the symbiotic relationship in D. citri.

Highly Reduced Complementary Genomes of Dual Bacterial Symbionts in the Mulberry Psyllid Anomoneura mori

The genomes of obligately host-restricted bacteria suffer from accumulating mildly deleterious mutations, resulting in marked size reductions. Psyllids (Hemiptera) are phloem sap-sucking insects with a specialized organ called the bacteriome, which typically harbors two vertically transmitted bacterial symbionts: the primary symbiont "Candidatus Carsonella ruddii" (Gammaproteobacteria) and a secondary symbiont that is phylogenetically diverse among psyllid lineages. The genomes of several Carsonella lineages were revealed to be markedly reduced (158-174‍ ‍kb), AT-rich (14.0-17.9% GC), and structurally conserved with similar gene inventories devoted to synthesizing essential amino acids that are scarce in the phloem sap. However, limited genomic information is currently available on secondary symbionts. Therefore, the present study investigated the genomes of the bacteriome-associated dual symbionts, Secondary_AM (Gammaproteobacteria) and Carsonella_AM, in the mulberry psyllid Anomoneura mori (Psyllidae). The results obtained revealed that the Secondary_AM genome is as small and AT-rich (229,822 bp, 17.3% GC) as those of Carsonella lineages, including Carsonella_AM (169,120 bp, 16.2% GC), implying that Secondary_AM is an evolutionarily ancient obligate mutualist, as is Carsonella. Phylogenomic ana-lyses showed that Secondary_AM is sister to "Candidatus Psyllophila symbiotica" of Cacopsylla spp. (Psyllidae), the genomes of which were recently reported (221-237‍ ‍kb, 17.3-18.6% GC). The Secondary_AM and Psyllophila genomes showed highly conserved synteny, sharing all genes for complementing the incomplete tryptophan biosynthetic pathway of Carsonella and those for synthesizing B vitamins. However, sulfur assimilation and carotenoid-synthesizing genes were only retained in Secondary_AM and Psyllophila, respectively, indicating ongoing gene silencing. Average nucleotide identity, gene ortholog similarity, genome-wide synteny, and substitution rates suggest that the Secondary_AM/Psyllophila genomes are more labile than Carsonella genomes.

Insect Microbial Symbionts: Ecology, Interactions, and Biological Significance

The guts of insect pests are typical habitats for microbial colonization and the presence of bacterial species inside the gut confers several potential advantages to the insects. These gut bacteria are located symbiotically inside the digestive tracts of insects and help in food digestion, phytotoxin breakdown, and pesticide detoxification. Different shapes and chemical assets of insect gastrointestinal tracts have a significant impact on the structure and makeup of the microbial population. The number of microbial communities inside the gastrointestinal system differs owing to the varying shape and chemical composition of digestive tracts. Due to their short generation times and rapid evolutionary rates, insect gut bacteria can develop numerous metabolic pathways and can adapt to diverse ecological niches. In addition, despite hindering insecticide management programs, they still have several biotechnological uses, including industrial, clinical, and environmental uses. This review discusses the prevalent bacterial species associated with insect guts, their mode of symbiotic interaction, their role in insecticide resistance, and various other biological significance, along with knowledge gaps and future perspectives. The practical consequences of the gut microbiome and its interaction with the insect host may lead to encountering the mechanisms behind the evolution of pesticide resistance in insects.

Division of labor within psyllids: metagenomics reveals an ancient dual endosymbiosis with metabolic complementarity in the genus Cacopsylla

IMPORTANCE

Heritable beneficial bacterial endosymbionts have been crucial for the evolutionary success of numerous insects by enabling the exploitation of nutritionally limited food sources. Herein, we describe a previously unknown dual endosymbiosis in the psyllid genus Cacopsylla, consisting of the primary endosymbiont "Candidatus Carsonella ruddii" and a co-occurring Enterobacteriaceae bacterium for which we propose the name "Candidatus Psyllophila symbiotica." Its localization within the bacteriome and its small genome size confirm that Psyllophila is a co-primary endosymbiont widespread within the genus Cacopsylla. Despite its highly eroded genome, Psyllophila perfectly complements the tryptophan biosynthesis pathway that is incomplete in the co-occurring Carsonella. Moreover, the genome of Psyllophila is almost as small as Carsonella's, suggesting an ancient dual endosymbiosis that has now reached a precarious stage where any additional gene loss would make the system collapse. Hence, our results shed light on the dynamic interactions of psyllids and their endosymbionts over evolutionary time.

Unveiling metabolic integration in psyllids and their nutritional endosymbionts through comparative transcriptomics analysis

Psyllids, a group of insects that feed on plant sap, have a symbiotic relationship with an endosymbiont called Carsonella. Carsonella synthesizes essential amino acids and vitamins for its psyllid host, but lacks certain genes required for this process, suggesting a compensatory role of psyllid host genes. To investigate this, gene expression was compared between two psyllid species, Bactericera cockerelli and Diaphorina citri, in specialized cells where Carsonella resides (bacteriomes). Collaborative psyllid genes, including horizontally transferred genes, showed patterns of conserved gene expression; however, species-specific patterns were also observed, suggesting differences in the nutritional metabolism between psyllid species. Also, the recycling of nitrogen in bacteriomes may primarily rely on glutamate dehydrogenase (GDH). Additionally, lineage-specific gene clusters were differentially expressed in B. cockerelli and D. citri bacteriomes and are highlighted here. These findings shed light on potential host adaptations for the regulation of this symbiosis due to host, microbiome, and environmental differences.

Seasonal wild dance of dual endosymbionts in the pear psyllid Cacopsylla pyricola (Hemiptera: Psylloidea)

Most sap-feeding insects maintain obligate relationships with endosymbiotic bacteria that provide their hosts with essential nutrients. However, knowledge about the dynamics of endosymbiont titers across seasons in natural host populations is scarce. Here, we used quantitative PCR to investigate the seasonal dynamics of the dual endosymbionts "Candidatus Carsonella ruddii" and "Ca. Psyllophila symbiotica" in a natural population of the pear psyllid Cacopsylla pyricola (Hemiptera: Psylloidea: Psyllidae). Psyllid individuals were collected across an entire year, covering both summer and overwintering generations. Immatures harboured the highest titers of both endosymbionts, while the lowest endosymbiont density was observed in males. The density of Carsonella remained high and relatively stable across the vegetative period of the pear trees, but significantly dropped during the non-vegetative period, overlapping with C. pyricola's reproductive diapause. In contrast, the titer of Psyllophila was consistently higher than Carsonella's and exhibited fluctuations throughout the sampling year, which might be related to host age. Despite a tightly integrated metabolic complementarity between Carsonella and Psyllophila, our findings highlight differences in their density dynamics throughout the year, that might be linked to their metabolic roles at different life stages of the host.

Variation of endosymbiont and citrus tristeza virus (CTV) titers in the Huanglongbing insect vector, Diaphorina citri, on CTV-infected plants

"Candidatus Liberibacter asiaticus" (CLas) is a notorious agent that causes Citrus Huanglongbing (HLB), which is transmitted by Diaphorina citri (D. citri). We recently found that the acquisition and transmission of CLas by D. citri was facilitated by Citrus tristeza virus (CTV), a widely distributed virus in the field. In this study, we further studied whether different CTV strains manipulate the host preference of D. citri, and whether endosymbionts variation is related to CTV strains in D. citri. The results showed that the non-viruliferous D. citri preferred to select the shoots infected with CTV, without strain differences was observed in the selection. However, the viruliferous D. citri prefered to select the mixed strain that is similar to the field's. Furthermore, D. citri effectively acquired the CTV within 2-12 h depending on the strains of the virus. The persistence period of CTV in D. citri was longer than 24 days, without reduction of the CTV titers being observed. These results provide a foundation for understanding the transmission mode of D. citri on CTV. During the process of CTV acquisition and persistence, the titers of main endosymbionts in D. citri showed similar variation trend, but their relative titers were different at different time points. The titers of the "Candidatus Profftella armatura" and CTV tended to be positively correlated, and the titers of Wolbachia and "Candidatus Carsonella ruddii" were mostly negatively related with titers of CT31. These results showed the relationship among D. citri, endosymbionts, and CTV and provided useful information for further research on the interactions between D. citri and CLas, which may benefit the development of approaches for the prevention of CLas transmission and control of citrus HLB.

Comparative genomics of Mollicutes-related endobacteria supports a late invasion into Mucoromycota fungi

Diverse members of early-diverging Mucoromycota, including mycorrhizal taxa and soil-associated Mortierellaceae, are known to harbor Mollicutes-related endobacteria (MRE). It has been hypothesized that MRE were acquired by a common ancestor and transmitted vertically. Alternatively, MRE endosymbionts could have invaded after the divergence of Mucoromycota lineages and subsequently spread to new hosts horizontally. To better understand the evolutionary history of MRE symbionts, we generated and analyzed four complete MRE genomes from two Mortierellaceae genera: Linnemannia (MRE-L) and Benniella (MRE-B). These genomes include the smallest known of fungal endosymbionts and showed signals of a tight relationship with hosts including a reduced functional capacity and genes transferred from fungal hosts to MRE. Phylogenetic reconstruction including nine MRE from mycorrhizal fungi revealed that MRE-B genomes are more closely related to MRE from Glomeromycotina than MRE-L from the same host family. We posit that reductions in genome size, GC content, pseudogene content, and repeat content in MRE-L may reflect a longer-term relationship with their fungal hosts. These data indicate Linnemannia and Benniella MRE were likely acquired independently after their fungal hosts diverged from a common ancestor. This work expands upon foundational knowledge on minimal genomes and provides insights into the evolution of bacterial endosymbionts.

Variability of Inverted Repeats in All Available Genomes of Bacteria

Noncanonical secondary structures in nucleic acids have been studied intensively in recent years. Important biological roles of cruciform structures formed by inverted repeats (IRs) have been demonstrated in diverse organisms, including humans. Using Palindrome analyser, we analyzed IRs in all accessible bacterial genome sequences to determine their frequencies, lengths, and localizations. IR sequences were identified in all species, but their frequencies differed significantly across various evolutionary groups. We detected 242,373,717 IRs in all 1,565 bacterial genomes. The highest mean IR frequency was detected in the Tenericutes (61.89 IRs/kbp) and the lowest mean frequency was found in the Alphaproteobacteria (27.08 IRs/kbp). IRs were abundant near genes and around regulatory, tRNA, transfer-messenger RNA (tmRNA), and rRNA regions, pointing to the importance of IRs in such basic cellular processes as genome maintenance, DNA replication, and transcription. Moreover, we found that organisms with high IR frequencies were more likely to be endosymbiotic, antibiotic producing, or pathogenic. On the other hand, those with low IR frequencies were far more likely to be thermophilic. This first comprehensive analysis of IRs in all available bacterial genomes demonstrates their genomic ubiquity, nonrandom distribution, and enrichment in genomic regulatory regions. IMPORTANCE Our manuscript reports for the first time a complete analysis of inverted repeats in all fully sequenced bacterial genomes. Thanks to the availability of unique computational resources, we were able to statistically evaluate the presence and localization of these important regulatory sequences in bacterial genomes. This work revealed a strong abundance of these sequences in regulatory regions and provides researchers with a valuable tool for their manipulation.

16S rRNA Gene Sequencing of Six Psyllid Species of the Family Carsidaridae Identified Various Bacteria Including Symbiopectobacterium

Psyllids (Hemiptera: Sternorrhyncha: Psylloidea) are plant sap-sucking insects that are closely associated with various microbes. To obtain a more detailed understanding of the ecological and evolutionary behaviors of microbes in Psylloidea, the bacterial populations of six psyllid species, belonging to the family Carsidaridae, were analyzed using high-throughput amplicon sequencing of the 16S rRNA gene. The majority of the secondary symbionts identified in the present study were gammaproteobacteria, particularly those of the order Enterobacterales, including Arsenophonus and Sodalis, which are lineages found in a wide variety of insect hosts. Additionally, Symbiopectobacterium, another Enterobacterales lineage, which has recently been recognized and increasingly shown to be vertically transmitted and mutualistic in various invertebrates, was identified for the first time in Psylloidea. This lineage is closely related to Pectobacterium spp., which are plant pathogens, but forms a distinct clade exhibiting no pathogenicity to plants. Non-Enterobacterales gammaproteobacteria found in the present study were Acinetobacter, Pseudomonas (both Pseudomonadales), Delftia, Comamonas (both Burkholderiales), and Xanthomonas (Xanthomonadales), a putative plant pathogen. Regarding alphaproteobacteria, three Wolbachia (Rickettsiales) lineages belonging to supergroup B, the major group in insect lineages, were detected in four psyllid species. In addition, a Wolbachia lineage of supergroup O, a minor group recently found for the first time in Psylloidea, was detected in one psyllid species. These results suggest the pervasive transfer of bacterial symbionts among animals and plants, providing deeper insights into the evolution of the interactions among these organisms.

Insect phylogeny structures the bacterial communities in the microbiome of psyllids (Hemiptera: Psylloidea) in Aotearoa New Zealand

The bacterial microbiome of psyllids has been studied for decades, with a strong focus on the primary and secondary endosymbionts capable of providing essential amino acids for the insects' diet and therefore playing a key role in the insects' ability to radiate on novel plant hosts. Here, we combine metabarcoding analysis of the bacterial communities hosted by psyllids with a multi-gene phylogenetic analysis of the insect hosts to determine what factors influence the bacterial diversity of the psyllids' microbiomes, especially in the context of the dispersal and evolutionary radiation of these insects in Aotearoa New Zealand. Using multi-gene phylogenetics with COI, 18S and EF-1α sequences from 102 psyllid species, we confirmed for the first time monophyly for all the six genera of native/endemic Aotearoa New Zealand psyllids, with indications that they derive from at least six dispersal events to the country. This also revealed that, after its ancestral arrival, the genus Powellia has radiated onto a larger and more diverse range of plants than either Psylla or Ctenarytaina, which is uncommon amongst monophyletic psyllids globally. DNA metabarcoding of the bacterial 16S gene here represents the largest dataset analysed to date from psyllids, including 246 individuals from 73 species. This provides novel evidence that bacterial diversity across psyllid species is strongly associated with psyllid phylogenetic structure, and to a lesser degree to their host plant association and geographic distribution. Furthermore, while the strongest co-phylogenetic signals were derived from the primary and secondary symbionts, a signal of phylosymbiosis was still retained among the remaining taxa of the bacterial microbiome, suggesting potential vertical transmission of bacterial lineages previously unknown to have symbiotic roles.

Microbiome of pear psyllids: A tale about closely related species sharing their endosymbionts

Psyllids are phloem-feeding insects that can transmit plant pathogens such as phytoplasmas, intracellular bacteria causing numerous plant diseases worldwide. Their microbiomes are essential for insect physiology and may also influence the capacity of vectors to transmit pathogens. Using 16S rRNA gene metabarcoding, we compared the microbiomes of three sympatric psyllid species associated with pear trees in Central Europe. All three species are able to transmit 'Candidatus Phytoplasma pyri', albeit with different efficiencies. Our results revealed potential relationships between insect biology and microbiome composition that varied during psyllid ontogeny and between generations in Cacopsylla pyri and C. pyricola, as well as between localities in C. pyri. In contrast, no variations related to psyllid life cycle and geography were detected in C. pyrisuga. In addition to the primary endosymbiont Carsonella ruddii, we detected another highly abundant endosymbiont (unclassified Enterobacteriaceae). C. pyri and C. pyricola shared the same taxon of Enterobacteriaceae which is related to endosymbionts harboured by other psyllid species from various families. In contrast, C. pyrisuga carried a different Enterobacteriaceae taxon related to the genus Sodalis. Our study provides new insights into host-symbiont interactions in psyllids and highlights the importance of host biology and geography in shaping microbiome structure.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

Investigating the microbial community of Cacopsylla spp. as potential factor in vector competence of phytoplasma

Phytoplasmas are obligatory intracellular bacteria that colonize the phloem of many plant species and cause hundreds of plant diseases worldwide. In nature, phytoplasmas are primarily transmitted by hemipteran vectors. While all phloem-feeding insects could in principle transmit phytoplasmas, only a limited number of species have been confirmed as vectors. Knowledge about factors that might determine the vector capacity is currently scarce. Here, we characterized the microbiomes of vector and non-vector species of apple proliferation (AP) phytoplasma 'Candidatus Phytoplasma mali' to investigate their potential role in the vector capacity of the host. We performed high-throughput 16S rRNA metabarcoding of the two principal AP-vectors Cacopsylla picta and Cacopsylla melanoneura and eight Cacopsylla species, which are not AP-vectors but co-occur in apple orchards. The microbiomes of all species are dominated by Carsonella, the primary endosymbiont of psyllids and a second uncharacterized Enterobacteriaceae endosymbiont. Each Cacopsylla species harboured a species-specific phylotype of both symbionts. Moreover, we investigated differences between the microbiomes of AP-vector versus non-vector species and identified the predominant endosymbionts but also Wolbachia and several minor taxa as potential indicator species. Our study highlights the importance of considering the microbiome in future investigations of potential factors influencing host vector competence. We investigated the potential role of symbiotic bacteria in the acquisition and transmission of phytoplasma. By comparing the two main psyillid vector species of Apple proliferation (AP) phytoplasma and eight co-occurring species, which are not able to vector AP-phytoplasma, we found differences in the microbial communities of AP-vector and non-vector species, which appear to be driven by the predominant symbionts in both vector species and Wolbachia and several minor taxa in the non-vector species. In contrast, infection with AP-phytoplasma did not affect microbiome composition in both vector species. Our study provides new insights into the endosymbiont diversity of Cacopsylla spp. and highlights the importance of considering the microbiome when investigating potential factors influencing host vector competence.

Convergent evolution of a labile nutritional symbiosis in ants

Ants are among the most successful organisms on Earth. It has been suggested that forming symbioses with nutrient-supplementing microbes may have contributed to their success, by allowing ants to invade otherwise inaccessible niches. However, it is unclear whether ants have evolved symbioses repeatedly to overcome the same nutrient limitations. Here, we address this question by comparing the independently evolved symbioses in Camponotus, Plagiolepis, Formica and Cardiocondyla ants. Our analysis reveals the only metabolic function consistently retained in all of the symbiont genomes is the capacity to synthesise tyrosine. We also show that in certain multi-queen lineages that have co-diversified with their symbiont for millions of years, only a fraction of queens carry the symbiont, suggesting ants differ in their colony-level reliance on symbiont-derived resources. Our results imply that symbioses can arise to solve common problems, but hosts may differ in their dependence on symbionts, highlighting the evolutionary forces influencing the persistence of long-term endosymbiotic mutualisms.

A complex interplay of evolutionary forces continues to shape ancient co-occurring symbiont genomes

Many insects depend on ancient associations with intracellular bacteria for essential nutrition. The genomes of these bacteria are often highly reduced. Although drift is a major driver of symbiont evolution, other evolutionary forces continue to influence them. To understand how ongoing molecular evolution and gene loss shape symbiont genomes, we sequenced two of the most ancient symbionts known, Sulcia and Nasuia, from 20 Hawaiian Nesophrosyne leafhoppers. We leveraged the parallel divergence of Nesophrosyne lineages throughout Hawaii as a natural experimental framework. Sulcia and Nasuia experience ongoing-but divergent-gene loss, often in a convergent fashion. Although some genes are under relaxed selection, purifying and positive selection are also important drivers of genome evolution, particularly in maintaining certain nutritional and cellular functions. Our results further demonstrate that symbionts experience dramatically different evolutionary environments, as evidenced by the finding that Sulcia and Nasuia have one of the slowest and fastest rates of molecular evolution known.

Chromosomal-level assembly of Bactericera cockerelli reveals rampant gene family expansions impacting genome structure, function and insect-microbe-plant-interactions

Lineage specific expansions and gene duplications are some of the most important sources of evolutionary novelty in eukaryotes. Although not as prevalent in eukaryotes compared to bacteria, horizontal gene transfer events can also result in key adaptations for insects, especially for those involved in insect-microbe interactions. In this study we assemble the first chromosomal assembly of the psyllid Bactericera cockerelli and reveal that the B. cockerelli genome has experienced significantly more gene expansion events compared to other Hemipteran representatives with fully sequenced genomes. We also reveal that B. cockerelli's genome is the largest psyllid genome (567 Mb) sequenced to date and is ~15% larger than the other two psyllid species genomes sequenced (Pachypsylla venusta and Diaphorina citri). Structurally, B. cockerelli appears to have an additional chromosome compared to the distantly related psyllid species P. venusta due to a previous chromosomal fission or fusion event. The increase in genome size and dynamic nature of the B. cockerelli genome may largely be contributed to the widespread expansion of type I and type II repeat elements that are rampant across all of B. cockerelli's. chromosomes. These repeat elements are distributed near equally in both euchromatic and heterochromatic regions. Furthermore, significant gene family expansions and gene duplications were uncovered for genes that are expected to be important in its adaptation to insect-plant and microbe interactions, which include transcription factors, proteases, odorant receptors, and horizontally transferred genes that are involved in the nutritional symbioses with their long-term nutritional endosymbiont Carsonella.

Divergent Host-Microbe Interaction and Pathogenesis Proteins Detected in Recently Identified Liberibacter Species

Candidatus (Ca.) Liberibacter taxa are economically important bacterial plant pathogens that are not culturable; however, genome-enabled insights can help us develop a deeper understanding of their host-microbe interactions and evolution. The draft genome of a recently identified Liberibacter taxa, Ca. Liberibacter capsica, was curated and annotated here with a total draft genome size of 1.1 MB with 1,036 proteins, which is comparable to other Liberibacter species with complete genomes. A total of 459 orthologous clusters were identified among Ca. L. capsica, Ca. L. asiaticus, Ca. L. psyllaurous, Ca. L. americanus, Ca. L. africanus, and L. crescens, and these genes within these clusters consisted of housekeeping and environmental response functions. We estimated the rates of molecular evolution for each of the 443 one-to-one ortholog clusters and found that all Ca. L. capsica orthologous pairs were under purifying selection when the synonymous substitutions per synonymous site (dS) were not saturated. These results suggest that these genes are largely maintaining their conserved functions. We also identified the most divergent single-copy orthologous proteins in Ca. L. capsica by analyzing the ortholog pairs that represented the highest nonsynonymous substitutions per nonsynonymous site (dN) values for each pairwise comparison. From these analyses, we found that 21 proteins which are known to be involved in pathogenesis and host-microbe interactions, including the Tad pilus complex, were consistently divergent between Ca. L. capsica and the majority of other Liberibacter species. These results further our understanding of the evolutionary genetics of Ca. L. capsica and, more broadly, the evolution of Liberibacter.

IMPORTANCE“Candidatus” (Ca.) Liberibacter taxa are economically important plant pathogens vectored by insects; however, these host-dependent bacterial taxa are extremely difficult to study because they are unculturable. Recently, we identified a new Ca. Liberibacter lineage (Ca. Liberibacter capsica) from a rare insect metagenomic sample. In this current study, we report that the draft genome of Ca. Liberibacter capsica is similar in genome size and protein content compared to the other Ca. Liberibacter taxa. We provide evidence that many of their shared genes, which encode housekeeping and environmental response functions, are evolving under purifying selection, suggesting that these genes are maintaining similar functions. Our study also identifies 21 proteins that are rapidly evolving amino acid changes in Ca. Liberibacter capsica compared to the majority of other Liberibacter taxa. Many of these proteins represent key genes involved in Liberibacter-host interactions and pathogenesis and are valuable candidate genes for future studies.

Developmental Integration of Endosymbionts in Insects

In endosymbiosis, two independently existing entities are inextricably intertwined such that they behave as a single unit. For multicellular hosts, the endosymbiont must be integrated within the host developmental genetic network to maintain the relationship. Developmental integration requires innovations in cell type, gene function, gene regulation, and metabolism. These innovations are contingent upon the existing ecological interactions and may evolve mutual interdependence. Recent studies have taken significant steps toward characterizing the proximate mechanisms underlying interdependence. However, the study of developmental integration is only in its early stages of investigation. Here, we review the literature on mutualistic endosymbiosis to explore how unicellular endosymbionts developmentally integrate into their multicellular hosts with emphasis on insects as a model. Exploration of this process will help gain a more complete understanding of endosymbiosis. This will pave the way for a better understanding of the endosymbiotic theory of evolution in the future.

Extreme Polyploidy of Carsonella , an Organelle-Like Bacterium with a Drastically Reduced Genome

Mitochondria and plastids are endosymbiotic organelles in eukaryotic cells and are derived from free-living bacteria. They have many highly reduced genomes from which numerous genes have been transferred to the host nucleus.

The Minimal Translation Machinery: What We Can Learn From Naturally and Experimentally Reduced Genomes

The current theoretical proposals of minimal genomes have not attempted to outline the essential machinery for proper translation in cells. Here, we present a proposal of a minimal translation machinery based on (1) a comparative analysis of bacterial genomes of insects’ endosymbionts using a machine learning classification algorithm, (2) the empiric genomic information obtained from Mycoplasma mycoides JCVI-syn3.0 the first minimal bacterial genome obtained by design and synthesis, and (3) a detailed functional analysis of the candidate genes based on essentiality according to the DEG database (Escherichia coli and Bacillus subtilis) and the literature. This proposed minimal translational machinery is composed by 142 genes which must be present in any synthetic prokaryotic cell designed for biotechnological purposes, 76.8% of which are shared with JCVI-syn3.0. Eight additional genes were manually included in the proposal for a proper and efficient translation.

Transitional genomes and nutritional role reversals identified for dual symbionts of adelgids (Aphidoidea: Adelgidae)

Many plant-sap-feeding insects have maintained a single, obligate, nutritional symbiont over the long history of their lineage. This senior symbiont may be joined by one or more junior symbionts that compensate for gaps in function incurred through genome-degradative forces. Adelgids are sap-sucking insects that feed solely on conifer trees and follow complex life cycles in which the diet fluctuates in nutrient levels. Adelgids are unusual in that both senior and junior symbionts appear to have been replaced repeatedly over their evolutionary history. Genomes can provide clues to understanding symbiont replacements, but only the dual symbionts of hemlock adelgids have been examined thus far. Here, we sequence and compare genomes of four additional dual-symbiont pairs in adelgids. We show that these symbionts are nutritional partners originating from diverse bacterial lineages and exhibiting wide variation in general genome characteristics. Although dual symbionts cooperate to produce nutrients, the balance of contributions varies widely across pairs, and total genome contents reflect a range of ages and degrees of degradation. Most symbionts appear to be in transitional states of genome reduction. Our findings support a hypothesis of periodic symbiont turnover driven by fluctuating selection for nutritional provisioning related to gains and losses of complex life cycles in their hosts.

Bacterial Endosymbionts of Bactericera maculipennis and Three Mitochondrial Haplotypes of B. cockerelli (Hemiptera: Psylloidea: Triozidae)

Insects harbor bacterial endosymbionts that provide their hosts with nutritional benefit or with protection against natural enemies, plant defenses, insecticides, or abiotic stresses. We used directed sequencing of 16S rDNA to identify and compare endosymbionts of Bactericera maculipennis (Crawford) and the western, central, and northwestern haplotypes of B. cockerelli (Šulc) (Hemiptera: Psylloidea: Triozidae). Both species are native to North America, are known to harbor the plant pathogen 'Candidatus Liberibacter solanacearum' and develop on shared host plants within the Convolvulaceae. The Old-World species Heterotrioza chenopodii (Reuter) (Psylloidea: Triozidae), now found in North America, was included as an outgroup. 16S sequencing confirmed that both Bactericera species harbor 'Candidatus Liberibacter solanacearum' and revealed that both species harbor unique strains of Wolbachia and Sodalis. However, the presence of Wolbachia and Sodalis varied among haplotypes of B. cockerelli. The central and western haplotypes harbored the same strains of Wolbachia, which was confirmed by Sanger sequencing of the wsp and ftsZ genes. Wolbachia was also detected in very low abundance from the northwestern haplotype by high-throughput sequencing of 16S but was not detected from this haplotype by PCR screening. The northwestern and central haplotypes also harbored Sodalis, which was not detected in the western haplotype. Heterotrioza chenopodii harbored an entirely different community of potential endosymbionts compared with the Bactericera spp. that included Rickettsia and an unidentified bacterium in the Enterobacteriaceae. Results of this study provide a foundation for further research on the interactions between psyllids and their bacterial endosymbionts.

Microbiome analyses of 12 psyllid species of the family Psyllidae identified various bacteria including Fukatsuia and Serratia symbiotica, known as secondary symbionts of aphids

Background

Psyllids (Hemiptera: Psylloidea) comprise a group of plant sap-sucking insects that includes important agricultural pests. They have close associations not only with plant pathogens, but also with various microbes, including obligate mutualists and facultative symbionts. Recent studies are revealing that interactions among such bacterial populations are important for psyllid biology and host plant pathology. In the present study, to obtain further insight into the ecological and evolutionary behaviors of bacteria in Psylloidea, we analyzed the microbiomes of 12 psyllid species belonging to the family Psyllidae (11 from Psyllinae and one from Macrocorsinae), using high-throughput amplicon sequencing of the 16S rRNA gene.

Results

The analysis showed that all 12 psyllids have the primary symbiont, Candidatus Carsonella ruddii (Gammaproteobacteria: Oceanospirillales), and at least one secondary symbiont. The majority of the secondary symbionts were gammaproteobacteria, especially those of the family Enterobacteriaceae (order: Enterobacteriales). Among them, symbionts belonging to “endosymbionts3”, which is a genus-level monophyletic group assigned by the SILVA rRNA database, were the most prevalent and were found in 9 of 11 Psyllinae species. Ca. Fukatsuia symbiotica and Serratia symbiotica, which were recognized only as secondary symbionts of aphids, were also identified. In addition to other Enterobacteriaceae bacteria, including Arsenophonus, Sodalis, and “endosymbionts2”, which is another genus-level clade, Pseudomonas (Pseudomonadales: Pseudomonadaceae) and Diplorickettsia (Diplorickettsiales: Diplorickettsiaceae) were identified. Regarding Alphaproteobacteria, the potential plant pathogen Ca. Liberibacter europaeus (Rhizobiales: Rhizobiaceae) was detected for the first time in Anomoneura mori (Psyllinae), a mulberry pest. Wolbachia (Rickettsiales: Anaplasmataceae) and Rickettsia (Rickettsiales: Rickettsiaceae), plausible host reproduction manipulators that are potential tools to control pest insects, were also detected.

Conclusions

The present study identified various bacterial symbionts including previously unexpected lineages in psyllids, suggesting considerable interspecific transfer of arthropod symbionts. The findings provide deeper insights into the evolution of interactions among insects, bacteria, and plants, which may be exploited to facilitate the control of pest psyllids in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02429-2.

Evolutionarily recent dual obligatory symbiosis among adelgids indicates a transition between fungus- and insect-associated lifestyles

Adelgids (Insecta: Hemiptera: Adelgidae) form a small group of insects but harbor a surprisingly diverse set of bacteriocyte-associated endosymbionts, which suggest multiple replacement and acquisition of symbionts over evolutionary time. Specific pairs of symbionts have been associated with adelgid lineages specialized on different secondary host conifers. Using a metagenomic approach, we investigated the symbiosis of the Adelges laricis/Adelges tardus species complex containing betaproteobacterial ("Candidatus Vallotia tarda") and gammaproteobacterial ("Candidatus Profftia tarda") symbionts. Genomic characteristics and metabolic pathway reconstructions revealed that Vallotia and Profftia are evolutionary young endosymbionts, which complement each other's role in essential amino acid production. Phylogenomic analyses and a high level of genomic synteny indicate an origin of the betaproteobacterial symbiont from endosymbionts of Rhizopus fungi. This evolutionary transition was accompanied with substantial loss of functions related to transcription regulation, secondary metabolite production, bacterial defense mechanisms, host infection, and manipulation. The transition from fungus to insect endosymbionts extends our current framework about evolutionary trajectories of host-associated microbes.

High-resolution Microbiome Analyses of Nine Psyllid Species of the Family Triozidae Identified Previously Unrecognized but Major Bacterial Populations, including Liberibacter and Wolbachia of Supergroup O

Psyllids (Hemiptera: Sternorrhyncha: Psylloidea) are plant sap-sucking insects that include important agricultural pests. To obtain insights into the ecological and evolutionary behaviors of microbes, including plant pathogens, in Psylloidea, high-resolution ana-lyses of the microbiomes of nine psyllid species belonging to the family Triozidae were performed using high-throughput amplicon sequencing of the 16S rRNA gene. Analyses identified various bacterial populations, showing that all nine psyllids have at least one secondary symbiont, along with the primary symbiont "Candidatus Carsonella ruddii" (Gammaproteobacteria: Oceanospirillales: Halomonadaceae). The majority of the secondary symbionts were gammaproteobacteria, particularly those of the order Enterobacterales, which included Arsenophonus and Serratia symbiotica, a bacterium formerly recognized only as a secondary symbiont of aphids (Hemiptera: Sternorrhyncha: Aphidoidea). The non-Enterobacterales gammaproteobacteria identified in the present study were Diplorickettsia (Diplorickettsiales: Diplorickettsiaceae), a potential human pathogen, and Carnimonas (Oceanospirillales: Halomonadaceae), a lineage detected for the first time in Psylloidea. Regarding alphaproteobacteria, the potential plant pathogen "Ca. Liberibacter europaeus" (Rhizobiales: Rhizobiaceae) was detected for the first time in Epitrioza yasumatsui, which feeds on the Japanese silverberry Elaeagnus umbellata (Elaeagnaceae), an aggressive invasive plant in the United States and Europe. Besides the detection of Wolbachia (Rickettsiales: Anaplasmataceae) of supergroup B in three psyllid species, a lineage belonging to supergroup O was identified for the first time in Psylloidea. These results suggest the rampant transfer of bacterial symbionts among animals and plants, thereby providing deeper insights into the evolution of interkingdom interactions among multicellular organisms and bacteria, which will facilitate the control of pest psyllids.

Uncovering Symbionts Across the Psyllid Tree of Life and the Discovery of a New Liberibacter Species, "Candidatus" Liberibacter capsica

Sap-feeding insects in the order Hemiptera associate with obligate endosymbionts that are required for survival and facultative endosymbionts that can potentially modify resistance to stress, enemies, development, and reproduction. In the superfamily Psylloidea, the jumping plant lice (psyllids), less is known about the diversity and prevalence of their endosymbionts compared to other sap-feeding pests such as aphids (Aphididae). To address this knowledge gap, using 16S rRNA sequencing we identify symbionts across divergent psyllid host lineages from around the world. Taking advantage of a new comprehensive phylogenomic analyses of Psylloidea, we included psyllid samples from 44 species of 35 genera of five families, collected from 11 international locations for this study. Across psyllid lineages, a total of 91 OTUs were recovered, predominantly of the Enterobacteriaceae (68%). The diversity of endosymbionts harbored by each psyllid species was low with an average of approximately 3 OTUs. Two clades of endosymbionts (clade 1 and 2), belonging to Enterobacteriaceae, were identified that appear to be long term endosymbionts of the psyllid families Triozidae and Psyllidae, respectively. We also conducted high throughput metagenomic sequencing on three Ca. Liberibacter infected psyllid species (Russelliana capsici, Trichochermes walkeri, and Macrohomotoma gladiata), initially identified from 16S rRNA sequencing, to obtain more genomic information on these putative Liberibacter plant pathogens. The phylogenomic analyses from these data identified a new Ca. Liberibacter species, Candidatus Liberibacter capsica, that is a potential pathogen of solanaceous crops. This new species shares a distant ancestor with Ca. L. americanus, which occurs in the same range as R. capsici in South America. We also detected the first association between a psyllid specializing on woody hosts and the Liberibacter species Ca. L. psyllaurous, which is a globally distributed pathogen of herbaceous crop hosts in the Solanaceae. Finally, we detected a potential association between a psyllid pest of figs (M. gladiata) and a Ca. Liberibacter related to Ca. L. asiaticus, which causes severe disease in citrus. Our findings reveal a wider diversity of associations between facultative symbionts and psyllids than previously reported and suggest numerous avenues for future work to clarify novel associations of ecological, evolutionary, and pathogenic interest.

Comparative Genomics Underlines Multiple Roles of Profftella, an Obligate Symbiont of Psyllids: Providing Toxins, Vitamins, and Carotenoids

The Asian citrus psyllid Diaphorina citri (Insecta: Hemiptera: Psylloidea), a serious pest of citrus species worldwide, harbors vertically transmitted intracellular mutualists, Candidatus Profftella armatura (Profftella_DC, Gammaproteobacteria: Burkholderiales) and Candidatus Carsonella ruddii (Carsonella_DC, Gammaproteobacteria: Oceanospirillales). Whereas Carsonella_DC is a typical nutritional symbiont, Profftella_DC is a unique defensive symbiont with organelle-like features, including intracellular localization within the host, perfect infection in host populations, vertical transmission over evolutionary time, and drastic genome reduction down to much less than 1 Mb. Large parts of the 460-kb genome of Profftella_DC are devoted to genes for synthesizing a polyketide toxin; diaphorin. To better understand the evolution of this unusual symbiont, the present study analyzed the genome of Profftella_Dco, a sister lineage to Profftella_DC, using Diaphorina cf. continua, a host psyllid congeneric with D. citri. The genome of coresiding Carsonella (Carsonella_Dco) was also analyzed. The analysis revealed nearly perfect synteny conservation in these genomes with their counterparts from D. citri. The substitution rate analysis further demonstrated genomic stability of Profftella which is comparable to that of Carsonella. Profftella_Dco and Profftella_DC shared all genes for the biosynthesis of diaphorin, hemolysin, riboflavin, biotin, and carotenoids, underlining multiple roles of Profftella, which may contribute to stabilizing symbiotic relationships with the host. However, acyl carrier proteins were extensively amplified in polyketide synthases DipP and DipT for diaphorin synthesis in Profftella_Dco. This level of acyl carrier protein augmentation, unprecedented in modular polyketide synthases of any known organism, is not thought to influence the polyketide structure but may improve the synthesis efficiency.

The Evolution of Interdependence in a Four-Way Mealybug Symbiosis

Mealybugs are insects that maintain intracellular bacterial symbionts to supplement their nutrient-poor plant sap diets. Some mealybugs have a single betaproteobacterial endosymbiont, a Candidatus Tremblaya species (hereafter Tremblaya) that alone provides the insect with its required nutrients. Other mealybugs have two nutritional endosymbionts that together provision these same nutrients, where Tremblaya has gained a gammaproteobacterial partner that resides in its cytoplasm. Previous work had established that Pseudococcus longispinus mealybugs maintain not one but two species of gammaproteobacterial endosymbionts along with Tremblaya. Preliminary genomic analyses suggested that these two gammaproteobacterial endosymbionts have large genomes with features consistent with a relatively recent origin as insect endosymbionts, but the patterns of genomic complementarity between members of the symbiosis and their relative cellular locations were unknown. Here, using long-read sequencing and various types of microscopy, we show that the two gammaproteobacterial symbionts of P. longispinus are mixed together within Tremblaya cells, and that their genomes are somewhat reduced in size compared with their closest nonendosymbiotic relatives. Both gammaproteobacterial genomes contain thousands of pseudogenes, consistent with a relatively recent shift from a free-living to an endosymbiotic lifestyle. Biosynthetic pathways of key metabolites are partitioned in complex interdependent patterns among the two gammaproteobacterial genomes, the Tremblaya genome, and horizontally acquired bacterial genes that are encoded on the mealybug nuclear genome. Although these two gammaproteobacterial endosymbionts have been acquired recently in evolutionary time, they have already evolved codependencies with each other, Tremblaya, and their insect host.

RNA polymerases in strict endosymbiont bacteria with extreme genome reduction show distinct erosions that might result in limited and differential promoter recognition

Strict endosymbiont bacteria present high degree genome reduction, retain smaller proteins, and in some instances, lack complete functional domains compared to free-living counterparts. Until now, the mechanisms underlying these genetic reductions are not well understood. In this study, the conservation of RNA polymerases, the essential machinery for gene expression, is analyzed in endosymbiont bacteria with extreme genome reductions. We analyzed the RNA polymerase subunits to identify and define domains, subdomains, and specific amino acids involved in precise biological functions known in Escherichia coli. We also perform phylogenetic analysis and three-dimensional models over four lineages of endosymbiotic proteobacteria with the smallest genomes known to date: Candidatus Hodgkinia cicadicola, Candidatus Tremblaya phenacola, Candidatus Tremblaya Princeps, Candidatus Nasuia deltocephalinicola, and Candidatus Carsonella ruddii. We found that some Hodgkinia strains do not encode for the RNA polymerase α subunit. The rest encode genes for α, β, β', and σ subunits to form the RNA polymerase. However, 16% shorter, on average, respect their orthologous in E. coli. In the α subunit, the amino-terminal domain is the most conserved. Regarding the β and β' subunits, both the catalytic core and the assembly domains are the most conserved. However, they showed compensatory amino acid substitutions to adapt to changes in the σ subunit. Precisely, the most erosive diversity occurs within the σ subunit. We identified broad amino acid substitution even in those recognizing and binding to the -10-box promoter element. In an overall conceptual image, the RNA polymerase from Candidatus Nasuia conserved the highest similarity with Escherichia coli RNA polymerase and their σ70. It might be recognizing the two main promoter elements (-10 and -35) and the two promoter accessory elements (-10 extended and UP-element). In Candidatus Carsonella, the RNA polymerase could recognize all the promoter elements except the -10-box extended. In Candidatus Tremblaya and Hodgkinia, due to the α carboxyl-terminal domain absence, they might not recognize the UP-promoter element. We also identified the lack of the β flap-tip helix domain in most Hodgkinia's that suggests the inability to bind the -35-box promoter element.

Natural selection shapes maintenance of orthologous sRNAs in divergent host-restricted bacterial genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.

At the Gate of Mutualism: Identification of Genomic Traits Predisposing to Insect-Bacterial Symbiosis in Pathogenic Strains of the Aphid Symbiont Serratia symbiotica

Mutualistic associations between insects and heritable bacterial symbionts are ubiquitous in nature. The aphid symbiont Serratia symbiotica is a valuable candidate for studying the evolution of bacterial symbiosis in insects because it includes a wide diversity of strains that reflect the diverse relationships in which bacteria can be engaged with insects, from pathogenic interactions to obligate intracellular mutualism. The recent discovery of culturable strains, which are hypothesized to resemble the ancestors of intracellular strains, provide an opportunity to study the mechanisms underlying bacterial symbiosis in its early stages. In this study, we analyzed the genomes of three of these culturable strains that are pathogenic to aphid hosts, and performed comparative genomic analyses including mutualistic host-dependent strains. All three genomes are larger than those of the host-restricted S. symbiotica strains described so far, and show significant enrichment in pseudogenes and mobile elements, suggesting that these three pathogenic strains are in the early stages of the adaptation to their host. Compared to their intracellular mutualistic relatives, the three strains harbor a greater diversity of genes coding for virulence factors and metabolic pathways, suggesting that they are likely adapted to infect new hosts and are a potential source of metabolic innovation for insects. The presence in their genomes of secondary metabolism gene clusters associated with the production of antimicrobial compounds and phytotoxins supports the hypothesis that S. symbiotia symbionts evolved from plant-associated strains and that plants may serve as intermediate hosts. Mutualistic associations between insects and bacteria are the result of independent transitions to endosymbiosis initiated by the acquisition of environmental progenitors. In this context, the genomes of free-living S. symbiotica strains provide a rare opportunity to study the inventory of genes held by bacterial associates of insects that are at the gateway to a host-dependent lifestyle.

Genome Reduction and Secondary Metabolism of the Marine Sponge-Associated Cyanobacterium Leptothoe

Sponges form symbiotic relationships with diverse and abundant microbial communities. Cyanobacteria are among the most important members of the microbial communities that are associated with sponges. Here, we performed a genus-wide comparative genomic analysis of the newly described marine benthic cyanobacterial genus Leptothoe (Synechococcales). We obtained draft genomes from Le. kymatousa TAU-MAC 1615 and Le. spongobia TAU-MAC 1115, isolated from marine sponges. We identified five additional Leptothoe genomes, host-associated or free-living, using a phylogenomic approach, and the comparison of all genomes showed that the sponge-associated strains display features of a symbiotic lifestyle. Le. kymatousa and Le. spongobia have undergone genome reduction; they harbored considerably fewer genes encoding for (i) cofactors, vitamins, prosthetic groups, pigments, proteins, and amino acid biosynthesis; (ii) DNA repair; (iii) antioxidant enzymes; and (iv) biosynthesis of capsular and extracellular polysaccharides. They have also lost several genes related to chemotaxis and motility. Eukaryotic-like proteins, such as ankyrin repeats, playing important roles in sponge-symbiont interactions, were identified in sponge-associated Leptothoe genomes. The sponge-associated Leptothoe stains harbored biosynthetic gene clusters encoding novel natural products despite genome reduction. Comparisons of the biosynthetic capacities of Leptothoe with chemically rich cyanobacteria revealed that Leptothoe is another promising marine cyanobacterium for the biosynthesis of novel natural products.

Virtual 2-D map of the fungal proteome

The molecular weight and isoelectric point (pI) of the proteins plays important role in the cell. Depending upon the shape, size, and charge, protein provides its functional role in different parts of the cell. Therefore, understanding to the knowledge of their molecular weight and charges is (pI) is very important. Therefore, we conducted a proteome-wide analysis of protein sequences of 689 fungal species (7.15 million protein sequences) and construct a virtual 2-D map of the fungal proteome. The analysis of the constructed map revealed the presence of a bimodal distribution of fungal proteomes. The molecular mass of individual fungal proteins ranged from 0.202 to 2546.166 kDa and the predicted isoelectric point (pI) ranged from 1.85 to 13.759 while average molecular weight of fungal proteome was 50.98 kDa. A non-ribosomal peptide synthase (RFU80400.1) found in Trichoderma arundinaceum was identified as the largest protein in the fungal kingdom. The collective fungal proteome is dominated by the presence of acidic rather than basic pI proteins and Leu is the most abundant amino acid while Cys is the least abundant amino acid. Aspergillus ustus encodes the highest percentage (76.62%) of acidic pI proteins while Nosema ceranae was found to encode the highest percentage (66.15%) of basic pI proteins. Selenocysteine and pyrrolysine amino acids were not found in any of the analysed fungal proteomes. Although the molecular weight and pI of the protein are of enormous important to understand their functional roles, the amino acid compositions of the fungal protein will enable us to understand the synonymous codon usage in the fungal kingdom. The small peptides identified during the study can provide additional biotechnological implication.

Pantothenate mediates the coordination of whitefly and symbiont fitness

Intracellular symbionts in insects often have reduced genomes. Host acquisition of genes from bacteria is an important adaptation that supports symbionts. However, the function of horizontally transferred genes in insect symbiosis remains largely unclear. The primary symbiont Portiera housed in bacteriocytes lacks pantothenate synthesis genes: panB and panC, which is presumably complemented by a fused gene panB-panC (hereafter panBC) horizontally transferred from bacteria in Bemisia tabaci MEAM1. We found panBC in many laboratory cultures, and species of B. tabaci shares a common evolutionary origin. We demonstrated that complementation with whitefly panBC rescued E. coli pantothenate gene knockout mutants. Portiera elimination decreased the pantothenate level and PanBC abundance in bacteriocytes, and reduced whitefly survival and fecundity. Silencing PanBC decreased the Portiera titer, reduced the pantothenate level, and decreased whitefly survival and fecundity. Supplementation with pantothenate restored the symbiont titer, PanBC level, and fitness of RNAi whiteflies. These data suggest that pantothenate synthesis requires cooperation and coordination of whitefly PanBC expression and Portiera. This host-symbiont co-regulation was mediated by the pantothenate level. Our findings demonstrated that pantothenate production, by the cooperation of a horizontally acquired, fused bacteria gene and Portiera, facilitates the coordination of whitefly and symbiont fitness. Thus, this study extends our understanding on the basis of complex host-symbiont interactions.