A bioinformatics approach to identifying Wolbachia infections in arthropods

Indexing and searching petabase-scale nucleotide resources

Searching vast and rapidly growing nucleotide content in resources, such as runs in the Sequence Read Archive and assemblies for whole-genome shotgun sequencing projects in GenBank, is currently impractical for most researchers. Here we present Pebblescout, a tool that navigates such content by providing indexing and search capabilities. Indexing uses dense sampling of the sequences in the resource. Search finds subjects (runs or assemblies) that have short sequence matches to a user query, with well-defined guarantees and ranks them using informativeness of the matches. We illustrate the functionality of Pebblescout by creating eight databases that index over 3.7 petabases. The web service of Pebblescout can be reached at https://pebblescout.ncbi.nlm.nih.gov . We show that for a wide range of query lengths, Pebblescout provides a data-driven way for finding relevant subsets of large nucleotide resources, reducing the effort for downstream analysis substantially. We also show that Pebblescout results compare favorably to MetaGraph and Sourmash.

Bacterial community and genome analysis of cytoplasmic incompatibility-inducing Wolbachia in American serpentine leafminer, Liriomyza trifolii

Liriomyza trifolii, an agricultural pest, is occasionally infected by Wolbachia. A Wolbachia strain present in Liriomyza trifolii is associated with cytoplasmic incompatibility (CI) effects, leading to the death of embryos resulting from incompatible crosses between antibiotic-treated or naturally Wolbachia-free strain females and Wolbachia-infected males. In this study, high-throughput sequencing of hypervariable rRNA genes was employed to characterize the bacterial community in Wolbachia-infected L. trifolii without antibiotic treatment. The analysis revealed that Wolbachia dominates the bacterial community in L. trifolii, with minor presence of Acinetobacter, Pseudomonas, and Limnobacter. To elucidate the genetic basis of the CI phenotype, metagenomic sequencing was also conducted to assemble the genome of the Wolbachia strain. The draft-genome of the Wolbachia strain wLtri was 1.35 Mbp with 34% GC content and contained 1,487 predicted genes. Notably, within the wLtri genome, there are three distinct types of cytoplasmic incompatibility factor (cif) genes: Type I, Type III, and Type V cifA;B. These genes are likely responsible for inducing the strong cytoplasmic incompatibility observed in L. trifolii.

Screening Host Genomic Data for Wolbachia Infections

Less than a decade ago, the production of Wolbachia genomic assemblies was tedious, time-consuming, and expensive. The production of Wolbachia genomic DNA free of contamination from host DNA, as required for Wolbachia-targeted sequencing, was then only possible after the amplification and extraction of a large amount of clonal Wolbachia DNA. However, as an endosymbiotic bacterium, Wolbachia does not grow outside the host cell environment, and large-scale recovery of the bacteria required mass rearing of their host, preferably clones of a single individual to avoid strain genetic diversity, or amplification of cell cultures infected with a single Wolbachia strain. Bacterial DNA could be separated from host DNA based on genomic size. Nowadays, the production of full Wolbachia genomes does not require the physical isolation of the bacterial strains from their respective hosts, and the bacterium is often sequenced as a by-catch of host genomic projects. Here, we provide a step-by-step protocol to (1) identify whether host genome projects contain reads from associated Wolbachia and (2) isolate/retrieve the Wolbachia reads from the rest of the sequenced material. We hope this simple protocol will support many projects aiming at studying diverse Wolbachia genome assemblies.

Geographical survey of the mycobiome and microbiome of Southern California glassy-winged sharpshooters

The glassy-winged sharpshooter, Homalodisca vitripennis Germar, is an invasive xylem-feeding leafhopper with a devastating economic impact on California agriculture through transmission of the plant pathogen, Xylella fastidiosa. While studies have focused on X. fastidiosa or known symbionts of H. vitripennis, little work has been done at the scale of the microbiome (the bacterial community) or mycobiome (the fungal community). Here, we characterize the mycobiome and the microbiome of H. vitripennis across Southern California and explore correlations with captivity and host insecticide resistance status. Using high-throughput sequencing of the ribosomal internal transcribed spacer 1 region and the 16S rRNA gene to profile the mycobiome and microbiome, respectively, we found that while the H. vitripennis mycobiome significantly varied across Southern California, the microbiome did not. We also observed a significant difference in both the mycobiome and microbiome between captive and wild H. vitripennis. Finally, we found that the mycobiome, but not the microbiome, was correlated with insecticide resistance status in wild H. vitripennis. This study serves as a foundational look at the H. vitripennis mycobiome and microbiome across Southern California. Future work should explore the putative link between microbes and insecticide resistance status and investigate whether microbial communities should be considered in H. vitripennis management practices. IMPORTANCE The glassy-winged sharpshooter is an invasive leafhopper that feeds on the xylem of plants and transmits the devastating pathogen, Xylella fastidiosa, resulting in significant economic damage to California's agricultural system. While studies have focused on this pathogen or obligate symbionts of the glassy-winged sharpshooter, there is limited knowledge of the bacterial and fungal communities that make up its microbiome and mycobiome. To address this knowledge gap, we explored the composition of the mycobiome and the microbiome of the glassy-winged sharpshooter across Southern California and identified differences associated with geography, captivity, and host insecticide resistance status. Understanding sources of variation in the microbial communities associated with the glassy-winged sharpshooter is an important consideration for developing management strategies to control this invasive insect. This study is a first step toward understanding the role microbes may play in the glassy-winged sharpshooter's resistance to insecticides.

Deep data mining reveals variable abundance and distribution of microbial reproductive manipulators within and among diverse host species

Bacterial symbionts that manipulate the reproduction of their hosts are important factors in invertebrate ecology and evolution, and are being leveraged for host biological control. Infection prevalence restricts which biological control strategies are possible and is thought to be strongly influenced by the density of symbiont infection within hosts, termed titer. Current methods to estimate infection prevalence and symbiont titers are low-throughput, biased towards sampling infected species, and rarely measure titer. Here we develop a data mining approach to estimate symbiont infection frequencies within host species and titers within host tissues. We applied this approach to screen ~32,000 publicly available sequence samples from the most common symbiont host taxa, discovering 2,083 arthropod and 119 nematode infected samples. From these data, we estimated that Wolbachia infects approximately 44% of all arthropod and 34% of all nematode species, while other reproductive manipulators only infect 1-8% of arthropod and nematode species. Although relative titers within hosts were highly variable within and between arthropod species, a combination of arthropod host species and Wolbachia strain explained approximately 36% of variation in Wolbachia titer across the dataset. To explore potential mechanisms for host control of symbiont titer, we leveraged population genomic data from the model system Drosophila melanogaster. In this host, we found a number of SNPs associated with titer in candidate genes potentially relevant to host interactions with Wolbachia. Our study demonstrates that data mining is a powerful tool to detect bacterial infections and quantify infection intensities, thus opening an array of previously inaccessible data for further analysis in host-symbiont evolution.

Genomic Assessment of the Contribution of the Wolbachia Endosymbiont of Eurosta solidaginis to Gall Induction

We explored the genome of the Wolbachia strain, wEsol, symbiotic with the plant-gall-inducing fly Eurosta solidaginis with the goal of determining if wEsol contributes to gall induction by its insect host. Gall induction by insects has been hypothesized to involve the secretion of the phytohormones cytokinin and auxin and/or proteinaceous effectors to stimulate cell division and growth in the host plant. We sequenced the metagenome of E. solidaginis and wEsol and assembled and annotated the genome of wEsol. The wEsol genome has an assembled length of 1.66 Mbp and contains 1878 protein-coding genes. The wEsol genome is replete with proteins encoded by mobile genetic elements and shows evidence of seven different prophages. We also detected evidence of multiple small insertions of wEsol genes into the genome of the host insect. Our characterization of the genome of wEsol indicates that it is compromised in the synthesis of dimethylallyl pyrophosphate (DMAPP) and S-adenosyl L-methionine (SAM), which are precursors required for the synthesis of cytokinins and methylthiolated cytokinins. wEsol is also incapable of synthesizing tryptophan, and its genome contains no enzymes in any of the known pathways for the synthesis of indole-3-acetic acid (IAA) from tryptophan. wEsol must steal DMAPP and L-methionine from its host and therefore is unlikely to provide cytokinin and auxin to its insect host for use in gall induction. Furthermore, in spite of its large repertoire of predicted Type IV secreted effector proteins, these effectors are more likely to contribute to the acquisition of nutrients and the manipulation of the host's cellular environment to contribute to growth and reproduction of wEsol than to aid E. solidaginis in manipulating its host plant. Combined with earlier work that shows that wEsol is absent from the salivary glands of E. solidaginis, our results suggest that wEsol does not contribute to gall induction by its host.

The toxins of vertically transmitted Spiroplasma

Vertically transmitted (VT) microbial symbionts play a vital role in the evolution of their insect hosts. A longstanding question in symbiont research is what genes help promote long-term stability of vertically transmitted lifestyles. Symbiont success in insect hosts is due in part to expression of beneficial or manipulative phenotypes that favor symbiont persistence in host populations. In Spiroplasma, these phenotypes have been linked to toxin and virulence domains among a few related strains. However, these domains also appear frequently in phylogenetically distant Spiroplasma, and little is known about their distribution across the Spiroplasma genus. In this study, we present the complete genome sequence of the Spiroplasma symbiont of Drosophila atripex, a non-manipulating member of the Ixodetis clade of Spiroplasma, for which genomic data are still limited. We perform a genus-wide comparative analysis of toxin domains implicated in defensive and reproductive phenotypes. From 12 VT and 31 non-VT Spiroplasma genomes, ribosome-inactivating proteins (RIPs), OTU-like cysteine proteases (OTUs), ankyrins, and ETX/MTX2 domains show high propensity for VT Spiroplasma compared to non-VT Spiroplasma. Specifically, OTU and ankyrin domains can be found only in VT-Spiroplasma, and RIP domains are found in all VT Spiroplasma and three non-VT Spiroplasma. These domains are frequently associated with Spiroplasma plasmids, suggesting a possible mechanism for dispersal and maintenance among heritable strains. Searching insect genome assemblies available on public databases uncovered uncharacterized Spiroplasma genomes from which we identified several spaid-like genes encoding RIP, OTU, and ankyrin domains, suggesting functional interactions among those domain types. Our results suggest a conserved core of symbiont domains play an important role in the evolution and persistence of VT Spiroplasma in insects.

Phylogenomic analysis of Wolbachia genomes from the Darwin Tree of Life biodiversity genomics project

The Darwin Tree of Life (DToL) project aims to sequence all described terrestrial and aquatic eukaryotic species found in Britain and Ireland. Reference genome sequences are generated from single individuals for each target species. In addition to the target genome, sequenced samples often contain genetic material from microbiomes, endosymbionts, parasites, and other cobionts. Wolbachia endosymbiotic bacteria are found in a diversity of terrestrial arthropods and nematodes, with supergroups A and B the most common in insects. We identified and assembled 110 complete Wolbachia genomes from 93 host species spanning 92 families by filtering data from 368 insect species generated by the DToL project. From 15 infected species, we assembled more than one Wolbachia genome, including cases where individuals carried simultaneous supergroup A and B infections. Different insect orders had distinct patterns of infection, with Lepidopteran hosts mostly infected with supergroup B, while infections in Diptera and Hymenoptera were dominated by A-type Wolbachia. Other than these large-scale order-level associations, host and Wolbachia phylogenies revealed no (or very limited) cophylogeny. This points to the occurrence of frequent host switching events, including between insect orders, in the evolutionary history of the Wolbachia pandemic. While supergroup A and B genomes had distinct GC% and GC skew, and B genomes had a larger core gene set and tended to be longer, it was the abundance of copies of bacteriophage WO who was a strong determinant of Wolbachia genome size. Mining raw genome data generated for reference genome assemblies is a robust way of identifying and analysing cobiont genomes and giving greater ecological context for their hosts.

A continental-scale survey of Wolbachia infections in blue butterflies reveals evidence of interspecific transfer and invasion dynamics

Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infection dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infections in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia presence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean = 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quantify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridization between butterfly lineages and complex dynamics between Wolbachia strains.

Genome sequencing and comparative analysis of Wolbachia strain wAlbA reveals Wolbachia-associated plasmids are common

Wolbachia are widespread maternally-transmitted bacteria of arthropods that often spread by manipulating their host's reproduction through cytoplasmic incompatibility (CI). Their invasive potential is currently being harnessed in field trials aiming to control mosquito-borne diseases. Wolbachia genomes commonly harbour prophage regions encoding the cif genes which confer their ability to induce CI. Recently, a plasmid-like element was discovered in wPip, a Wolbachia strain infecting Culex mosquitoes; however, it is unclear how common such extra-chromosomal elements are in Wolbachia. Here we sequenced the complete genome of wAlbA, a strain of the symbiont found in Aedes albopictus, after eliminating the co-infecting and higher density wAlbB strain that previously made sequencing of wAlbA challenging. We show that wAlbA is associated with two new plasmids and identified additional Wolbachia plasmids and related chromosomal islands in over 20% of publicly available Wolbachia genome datasets. These plasmids encode a variety of accessory genes, including several phage-like DNA packaging genes as well as genes potentially contributing to host-symbiont interactions. In particular, we recovered divergent homologues of the cif genes in both Wolbachia- and Rickettsia-associated plasmids. Our results indicate that plasmids are common in Wolbachia and raise fundamental questions around their role in symbiosis. In addition, our comparative analysis provides useful information for the future development of genetic tools to manipulate and study Wolbachia symbionts.

One's trash is someone else's treasure: sequence read archives from Lepidoptera genomes provide material for genome reconstruction of their endosymbionts

BACKGROUND

Maternally inherited bacterial symbionts are extremely widespread in insects. They owe their success to their ability to promote their own transmission through various manipulations of their hosts' life-histories. Many symbionts however very often go undetected. Consequently, we have only a restricted idea of the true symbiont diversity in insects, which may hinder our understanding of even bigger questions in the field such as the evolution or establishment of symbiosis.

RESULTS

In this study, we screened publicly available Lepidoptera genomic material for two of the most common insect endosymbionts, namely Wolbachia and Spiroplasma, in 1904 entries, encompassing 106 distinct species. We compared the performance of two screening software, Kraken2 and MetaPhlAn2, to identify the bacterial infections and using a baiting approach we reconstruct endosymbiont genome assemblies. Of the 106 species screened, 20 (19%) and nine (8.5%) were found to be infected with either Wolbachia or Spiroplasma, respectively. Construction of partial symbiotic genomes and phylogenetic analyses suggested the Wolbachia strains from the supergroup B were the most prevalent type of symbionts, while Spiroplasma infections were scarce in the Lepidoptera species screened here.

CONCLUSIONS

Our results indicate that many of the host-symbiont associations remain largely unexplored, with the majority of associations we identify never being recorded before. This highlights the usefulness of public databases to explore the hidden diversity of symbiotic entities, allowing the development of hypotheses regarding host-symbiont associations. The ever-expanding genomic databases provide a diverse databank from which one can characterize and explore the true diversity of symbiotic entities.

Diaphorin, a Polyketide Produced by a Bacterial Symbiont of the Asian Citrus Psyllid, Inhibits the Growth and Cell Division of Bacillus subtilis but Promotes the Growth and Metabolic Activity of Escherichia coli

Diaphorin is a polyketide produced by “Candidatus Profftella armatura” (Gammaproteobacteria: Burkholderiales), an obligate symbiont of a notorious agricultural pest, the Asian citrus psyllid Diaphorina citri (Hemiptera: Psyllidae). Diaphorin belongs to the pederin family of bioactive agents found in various host-symbiont systems, including beetles, lichens, and sponges, harboring phylogenetically diverse bacterial producers. Previous studies showed that diaphorin, which is present in D. citri at concentrations of 2 to 20 mM, has inhibitory effects on various eukaryotes, including the natural enemies of D. citri. However, little is known about its effects on prokaryotic organisms. To address this issue, the present study assessed the biological activities of diaphorin on two model prokaryotes, Escherichia coli (Gammaproteobacteria: Enterobacterales) and Bacillus subtilis (Firmicutes: Bacilli). Their growth and morphological features were analyzed using spectrophotometry, optical microscopy followed by image analysis, and transmission electron microscopy. The metabolic activity of E. coli was further assessed using the β-galactosidase assay. The results revealed that physiological concentrations of diaphorin inhibit the growth and cell division of B. subtilis but promote the growth and metabolic activity of E. coli. This finding implies that diaphorin functions as a defensive agent of the holobiont (host plus symbionts) against some bacterial lineages but is metabolically beneficial for others, which potentially include obligate symbionts of D. citri.

IMPORTANCE Certain secondary metabolites, including antibiotics, evolve to mediate interactions among organisms. These molecules have distinct spectra for microorganisms and are often more effective against Gram-positive bacteria than Gram-negative ones. However, it is rare that a single molecule has completely opposite activities on distinct bacterial lineages. The present study revealed that a secondary metabolite synthesized by an organelle-like bacterial symbiont of psyllids inhibits the growth of Gram-positive Bacillus subtilis but promotes the growth of Gram-negative Escherichia coli. This finding not only provides insights into the evolution of microbiomes in animal hosts but also may potentially be exploited to promote the effectiveness of industrial material production by microorganisms.

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.

Molecular detection and characterization of the endosymbiont Wolbachia in the European hedgehog flea, Archaeopsylla erinacei

Wolbachia, the endosymbiont of arthropods and onchocercid nematodes is present in many medically important insect species, being also considered for the indirect control of parasitic ones. Archaeopsylla erinacei is a flea species infesting hedgehogs acting as vector of Rickettsia felis, Bartonella henselae, and Rickettsia helvetica, thus having public health relevance. The Wolbachia surface protein (wsp) and 16S rRNA genes were used to determine the presence, prevalence and molecular typing of Wolbachia in this flea species collected in two regions of southern Italy. Of the 45 fleas tested (n = 16 males, 35.6%; n = 29 females, 64.4%), 43 (95.6%; 95% CI: 84.8-99.2) scored positive for Wolbachia, of which 15 (33.3%) and 28 (62.2%) were males and females, respectively. The sex-wise prevalence of this endosymbiont was almost equal in both sexes (males 93.8%; 95% CI: 69.5-99.7; females 96.7%; 95% CI: 83.1-99.8). Single locus sequence analysis (SLST) of Wolbachia revealed two sequence types for 16S rRNA gene, named as wAr_15227 and wAr_15234, which came from two different areas, equally distributed in male and female fleas, whilst only one sequence type was identified for wsp gene. The phylogenetic analysis placed the two 16S rRNA sequence types in paraphyletic clades belonging to the supergroup A and B, respectively. Whilst, the tree of wsp gene clustered the corresponding sequence in the same clade including those of Wolbachia supergroup A. In MLST analyses, both Wolbachia sequence types clustered in a monophyletic clade with Drosophila nikananu (wNik) and Drosophila sturtevanti (wStv) from supergroup A. ClonalFrame analysis revealed a recombination event in the wAr_15234 strain which came from Apulia region. Scientific knowledge of the presence/prevalence of Wolbachia among medically important fleas, may contribute to develop an alternative biological method for the vector control.

Mito-nuclear discordance at a mimicry color transition zone in bumble bee Bombus melanopygus

As hybrid zones exhibit selective patterns of gene flow between otherwise distinct lineages, they can be especially valuable for informing processes of microevolution and speciation. The bumble bee, Bombus melanopygus, displays two distinct color forms generated by Müllerian mimicry: a northern "Rocky Mountain'' color form with ferruginous mid-abdominal segments (B. m. melanopygus) and a southern "Pacific'' form with black mid-abdominal segments (B. m. edwardsii). These morphs meet in a mimetic transition zone in northern California and southern Oregon that is more narrow and transitions further west than comimetic bumble bee species. To understand the historical formation of this mimicry zone, we assessed color distribution data for B. melanopygus from the last 100 years. We then examined gene flow among the color forms in the transition zone by comparing sequences from mitochondrial COI barcode sequences, color-controlling loci, and the rest of the nuclear genome. These data support two geographically distinct mitochondrial haplogroups aligned to the ancestrally ferruginous and black forms that meet within the color transition zone. This clustering is also supported by the nuclear genome, which, while showing strong admixture across individuals, distinguishes individuals most by their mitochondrial haplotype, followed by geography. These data suggest the two lineages most likely were historically isolated, acquired fixed color differences, and then came into secondary contact with ongoing gene flow. The transition zone, however, exhibits asymmetries: mitochondrial haplotypes transition further south than color pattern, and both transition over shorter distances in the south. This system thus demonstrates alternative patterns of gene flow that occur in contact zones, presenting another example of mito-nuclear discordance. Discordant gene flow is inferred to most likely be driven by a combination of mimetic selection, dominance effects, and assortative mating.

Next-Generation Sequencing and Its Impacts on Entomological Research in Ecology and Evolution

The advent of NGS-based methods has been profoundly transforming entomological research. Through continual development and improvement of different methods and sequencing platforms, NGS has promoted mass elucidation of partial or whole genetic materials associated with beneficial insects, pests (of agriculture, forestry and animal, and human health), and species of conservation concern, helping to unravel ecological and evolutionary mechanisms and characterizing survival, trophic interactions, and dispersal. It is shifting the scale of biodiversity and environmental analyses from individuals and biodiversity indicator species to the large-scale study of communities and ecosystems using bulk samples of species or a mixed "soup" of environmental DNA. As the NGS-based methods have become more affordable, complexity demystified, and specificity and sensitivity proven, their use in entomological research has spread widely. This article presents several examples on how NGS-based methods have been used in entomology to provide incentives to apply them when appropriate and to open our minds to the expected advances in entomology that are yet to come.

Improved draft reference genome for the Glassy-winged Sharpshooter (Homalodisca vitripennis), a vector for Pierce's disease

Homalodisca vitripennis (Hemiptera: Cicadellidae), known as the glassy-winged sharpshooter, is a xylem feeding leafhopper and an important agricultural pest as a vector of Xylella fastidiosa, which causes Pierce’s disease in grapes and a variety of other scorch diseases. The current H. vitripennis reference genome from the Baylor College of Medicine's i5k pilot project is a 1.4-Gb assembly with 110,000 scaffolds, which still has significant gaps making identification of genes difficult. To improve on this effort, we used a combination of Oxford Nanopore long-read sequencing technology combined with Illumina sequencing reads to generate a better assembly and first-pass annotation of the whole genome sequence of a wild-caught Californian (Tulare County) individual of H. vitripennis. The improved reference genome assembly for H. vitripennis is 1.93-Gb in length (21,254 scaffolds, N50 = 650 Mb, BUSCO completeness = 94.3%), with 33.06% of the genome masked as repetitive. In total, 108,762 gene models were predicted including 98,296 protein-coding genes and 10,466 tRNA genes. As an additional community resource, we identified 27 orthologous candidate genes of interest for future experimental work including phenotypic marker genes like white. Furthermore, as part of the assembly process, we generated four endosymbiont metagenome-assembled genomes, including a high-quality near complete 1.7-Mb Wolbachia sp. genome (1 scaffold, CheckM completeness = 99.4%). The improved genome assembly and annotation for H. vitripennis, curated set of candidate genes, and endosymbiont MAGs will be invaluable resources for future research of H. vitripennis.

Two New Strains of Wolbachia Affecting Natural Avocado Thrips

Wolbachia is an obligate intracellular bacterium with a high frequency of infection and a continental distribution in arthropods and nematodes. This endosymbiont can induce various reproductive phenotypes in their hosts and has been previously found naturally in several pests including thrips (Thripidae). These insects cause physical fruit damage and economic losses in avocado. The presence of Wolbachia was evaluated for the first time in avocado thrips populations of Frankliniella sp. and Scirtothrips hansoni sp.n. from eastern Antioquia. DNA from adult thrips individuals was used to assess the detection of Wolbachia by amplifying a fragment (600 bp) of the Wolbachia major surface protein (wsp) gene. Results confirmed the presence of two new Wolbachia strains in these two thrips species, with a higher percentage of natural infection in S. hansoni sp.n. The first Wolbachia species was found in Frankliniella sp. and belongs to supergroup A and the second was detected in S. hansoni sp.n. and is part of supergroup B. Wolbachia was more frequently found in females (32.73%), and only found in one male. Analysis of phylogenetic relationships, suggests that the two new Wolbachia sequences (wFran: Frankliniella and wShan: Scirtothrips hansoni) detected here represent two new groups for this endosymbiont. The haplotype network shows the presence of two possible haplotypes for each strain. Future studies to evaluate the possible use of Wolbachia as a control agent in avocado thrips are necessary.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12088-021-00951-5.

Evolutionary transition to XY sex chromosomes associated with Y-linked duplication of a male hormone gene in a terrestrial isopod

Sex chromosomes are highly variable in some taxonomic groups, but the evolutionary mechanisms underlying this diversity are not well understood. In terrestrial isopod crustaceans, evolutionary turnovers in sex chromosomes are frequent, possibly caused by Wolbachia, a vertically-transmitted endosymbiont causing male-to-female sex reversal. Here, we use surgical manipulations and genetic crosses, plus genome sequencing, to examine sex chromosomes in the terrestrial isopod Trachelipus rathkei. Although an earlier cytogenetics study suggested a ZZ/ZW sex chromosome system in this species, we surprisingly find multiple lines of evidence that in our study population, sex is determined by an XX/XY system. Consistent with a recent evolutionary origin for this XX/XY system, the putative male-specific region of the genome is small. The genome shows evidence of Y-linked duplications of the gene encoding the androgenic gland hormone, a major component of male sexual differentiation in isopods. Our analyses also uncover sequences horizontally acquired from past Wolbachia infections, consistent with the hypothesis that Wolbachia may have interfered with the evolution of sex determination in T. rathkei. Overall, these results provide evidence for the co-occurrence of multiple sex chromosome systems within T. rathkei, further highlighting the relevance of terrestrial isopods as models for the study of sex chromosome evolution.

Wolbachia: endosymbiont of onchocercid nematodes and their vectors

Background

Wolbachia is an obligate intracellular maternally transmitted, gram-negative bacterium which forms a spectrum of endosymbiotic relationships from parasitism to obligatory mutualism in a wide range of arthropods and onchocercid nematodes, respectively. In arthropods Wolbachia produces reproductive manipulations such as male killing, feminization, parthenogenesis and cytoplasmic incompatibility for its propagation and provides an additional fitness benefit for the host to protect against pathogens, whilst in onchocercid nematodes, apart from the mutual metabolic dependence, this bacterium is involved in moulting, embryogenesis, growth and survival of the host.

Methods

This review details the molecular data of Wolbachia and its effect on host biology, immunity, ecology and evolution, reproduction, endosymbiont-based treatment and control strategies exploited for filariasis. Relevant peer-reviewed scientic papers available in various authenticated scientific data bases were considered while writing the review.

Conclusions

The information presented provides an overview on Wolbachia biology and its use in the control and/or treatment of vectors, onchocercid nematodes and viral diseases of medical and veterinary importance. This offers the development of new approaches for the control of a variety of vector-borne diseases.

Graphic Abstract

Stable high-density and maternally inherited Wolbachia infections in Anopheles moucheti and Anopheles demeilloni mosquitoes

Wolbachia, a widespread bacterium that can reduce pathogen transmission in mosquitoes, has recently been reported to be present in Anopheles (An.) species. In wild populations of the An. gambiae complex, the primary vectors of Plasmodium malaria in Sub-Saharan Africa, Wolbachia DNA sequences at low density and infection frequencies have been detected. As the majority of studies have used highly sensitive nested PCR as the only method of detection, more robust evidence is required to determine whether Wolbachia strains are established as endosymbionts in Anopheles species. Here, we describe high-density Wolbachia infections in geographically diverse populations of An. moucheti and An. demeilloni. Fluorescent in situ hybridization localized a heavy infection in the ovaries of An. moucheti, and maternal transmission was observed. Genome sequencing of both Wolbachia strains obtained genome depths and coverages comparable to those of other known infections. Notably, homologs of cytoplasmic incompatibility factor (cif) genes were present, indicating that these strains possess the capacity to induce the cytoplasmic incompatibility phenotype, which allows Wolbachia to spread through host populations. These strains should be further investigated as candidates for use in Wolbachia biocontrol strategies in Anopheles aiming to reduce the transmission of malaria.

•

High-density Wolbachia strains found in An. moucheti and An. demeilloni mosquitoes

•

Infections are visualized in the ovaries, and maternal transmission was observed

•

Sequencing at depths and coverages comparable to other known Wolbachia strains

•

Homologs of cytoplasmic incompatibility factor genes are present in both genomes

Life and Death of Selfish Genes: Comparative Genomics Reveals the Dynamic Evolution of Cytoplasmic Incompatibility

Cytoplasmic incompatibility is a selfish reproductive manipulation induced by the endosymbiont Wolbachia in arthropods. In males Wolbachia modifies sperm, leading to embryonic mortality in crosses with Wolbachia-free females. In females, Wolbachia rescues the cross and allows development to proceed normally. This provides a reproductive advantage to infected females, allowing the maternally transmitted symbiont to spread rapidly through host populations. We identified homologs of the genes underlying this phenotype, cifA and cifB, in 52 of 71 new and published Wolbachia genome sequences. They are strongly associated with cytoplasmic incompatibility. There are up to seven copies of the genes in each genome, and phylogenetic analysis shows that Wolbachia frequently acquires new copies due to pervasive horizontal transfer between strains. In many cases, the genes have subsequently acquired loss-of-function mutations to become pseudogenes. As predicted by theory, this tends to occur first in cifB, whose sole function is to modify sperm, and then in cifA, which is required to rescue the cross in females. Although cif genes recombine, recombination is largely restricted to closely related homologs. This is predicted under a model of coevolution between sperm modification and embryonic rescue, where recombination between distantly related pairs of genes would create a self-incompatible strain. Together, these patterns of gene gain, loss, and recombination support evolutionary models of cytoplasmic incompatibility.

Comparative analysis reveals the expansion of mitochondrial DNA control region containing unusually high G-C tandem repeat arrays in Nasonia vitripennis

Insect mitochondrial DNA (mtDNA) ranges from 14 to 19 kbp, and the size difference is attributed to the AT-rich control region. Jewel wasps have a parasitoid lifestyle, which may affect mitochondria function and evolution. We sequenced, assembled, and annotated mitochondrial genomes in Nasonia and outgroup species. Gene composition and order are conserved within Nasonia, but they differ from other parasitoids by two large inversion events that were not reported before. We observed a much higher substitution rate relative to the nuclear genome and mitochondrial introgression between N. giraulti and N. oneida, which is consistent with previous studies. Most strikingly, N. vitripennis mtDNA has an extremely long control region (7665 bp), containing twenty-nine 217 bp tandem repeats and can fold into a super-cruciform structure. In contrast to tandem repeats commonly found in other mitochondria, these high-copy repeats are highly conserved (98.7% sequence identity), much longer in length (approximately 8 Kb), extremely GC-rich (50.7%), and CpG-rich (percent CpG 19.4% vs. 1.1% in coding region), resulting in a 23 kbp mtDNA beyond the typical size range in insects. These N. vitripennis-specific mitochondrial repeats are not related to any known sequences in insect mitochondria. Their evolutionary origin and functional consequences warrant further investigations.

Large scale genome reconstructions illuminate Wolbachia evolution

Wolbachia is an iconic example of a successful intracellular bacterium. Despite its importance as a manipulator of invertebrate biology, its evolutionary dynamics have been poorly studied from a genomic viewpoint. To expand the number of Wolbachia genomes, we screen over 30,000 publicly available shotgun DNA sequencing samples from 500 hosts. By assembling over 1000 Wolbachia genomes, we provide a substantial increase in host representation. Our phylogenies based on both core-genome and gene content provide a robust reference for future studies, support new strains in model organisms, and reveal recent horizontal transfers amongst distantly related hosts. We find various instances of gene function gains and losses in different super-groups and in cytoplasmic incompatibility inducing strains. Our Wolbachia-host co-phylogenies indicate that horizontal transmission is widespread at the host intraspecific level and that there is no support for a general Wolbachia-mitochondrial synchronous divergence.

By greatly expanding the number of assembled genomes for Wolbachia (a group of intracellular bacteria) and constructing robust phylogenies, this study finds strong rate heterogeneity among Wolbachiapopulations and no support for synchronous divergence between Wolbachia and host mitochondria.

16S rRNA Sequencing Detected Profftella, Liberibacter, Wolbachia, and Diplorickettsia from Relatives of the Asian Citrus Psyllid

The Asian citrus psyllid Diaphorina citri (Hemiptera: Psylloidea) is a serious pest of citrus species worldwide because it transmits Candidatus Liberibacter spp. (Alphaproteobacteria: Rhizobiales), the causative agents of the incurable citrus disease, huanglongbing or greening disease. Diaphorina citri possesses a specialized organ called a bacteriome, which harbors vertically transmitted intracellular mutualists, Ca. Carsonella ruddii (Gammaproteobacteria: Oceanospirillales) and Ca. Profftella armatura (Gammaproteobacteria: Betaproteobacteriales). Whereas Carsonella is a typical nutritional symbiont, Profftella is an unprecedented type of toxin-producing defensive symbiont, unusually sharing organelle-like features with nutritional symbionts. Additionally, many D. citri strains are infected with Wolbachia, which manipulate reproduction in various arthropod hosts. In the present study, in an effort to obtain insights into the evolution of symbioses between Diaphorina and bacteria, microbiomes of psyllids closely related to D. citri were investigated. Bacterial populations of Diaphorina cf. continua and Diaphorina lycii were analyzed using Illumina sequencing of 16S rRNA gene amplicons and compared with data obtained from D. citri. The analysis revealed that all three Diaphorina spp. harbor Profftella as well as Carsonella lineages, implying that Profftella is widespread within the genus Diaphorina. Moreover, the analysis identified Ca. Liberibacter europaeus and Diplorickettsia sp. (Gammaproteobacteria: Diplorickettsiales) in D. cf. continua, and a total of four Wolbachia (Alphaproteobacteria: Rickettsiales) lineages in the three psyllid species. These results provide deeper insights into the interactions among insects, bacteria, and plants, which would eventually help to better manage horticulture.

Molecular characterization of Cardinium, Rickettsia, Spiroplasma and Wolbachia in mite species from citrus orchards

Tetranychidae spider mites are considered key citrus pests in some production areas, especially Tetranychus urticae Koch. Over the past decades, pesticide overuse seems to have promoted T. urticae population selection in citrus orchards. However, the microbiota has also been pointed out as a plausible explanation for population structure or plant host specialisation observed in several arthropod species. In this work, we have determined the incidence of Cardinium, Rickettsia, Spiroplasma and Wolbachia as representatives of major distorter bacteria genera in Aplonobia histricina (Berlese), Eutetranychus banksi (McGregor), Eutetranychus orientalis (Klein), Panonychus citri (McGregor), Tetranychus evansi Baker and Pritchard, Tetranychus turkestani Ugarov and Nikolskii, and T. urticae populations from Spanish citrus orchards. Only Wolbachia was detected by PCR. The multilocus alignment approach and phylogenetic inference indicated that all detected Wolbachia belong to supergroup B. The deep analysis of each 16S rDNA, ftsZ and wsp gene sequences allowed identifying several phylogenetically different Wolbachia sequences. It probably indicates the presence of several different races or strains, all of them belonging to supergroup B. The wsp sequence typing analysis unveiled the presence of the two already identified alleles (61 and 370) and allowed to contribute with five new alleles, supporting the presence of different but related B-races in the studied mite populations. The results are discussed and related to T. urticae population structure, previously observed in Spanish citrus orchards.

Host-associated microbiomes are predicted by immune system complexity and climate

BACKGROUND

Host-associated microbiomes, the microorganisms occurring inside and on host surfaces, influence evolutionary, immunological, and ecological processes. Interactions between host and microbiome affect metabolism and contribute to host adaptation to changing environments. Meta-analyses of host-associated bacterial communities have the potential to elucidate global-scale patterns of microbial community structure and function. It is possible that host surface-associated (external) microbiomes respond more strongly to variations in environmental factors, whereas internal microbiomes are more tightly linked to host factors.

RESULTS

Here, we use the dataset from the Earth Microbiome Project and accumulate data from 50 additional studies totaling 654 host species and over 15,000 samples to examine global-scale patterns of bacterial diversity and function. We analyze microbiomes from non-captive hosts sampled from natural habitats and find patterns with bioclimate and geophysical factors, as well as land use, host phylogeny, and trophic level/diet. Specifically, external microbiomes are best explained by variations in mean daily temperature range and precipitation seasonality. In contrast, internal microbiomes are best explained by host factors such as phylogeny/immune complexity and trophic level/diet, plus climate.

CONCLUSIONS

Internal microbiomes are predominantly associated with top-down effects, while climatic factors are stronger determinants of microbiomes on host external surfaces. Host immunity may act on microbiome diversity through top-down regulation analogous to predators in non-microbial ecosystems. Noting gaps in geographic and host sampling, this combined dataset represents a global baseline available for interrogation by future microbial ecology studies.

Complete Genome Assemblies for Three Variants of the Wolbachia Endosymbiont of Drosophila melanogaster

Here, we report genome assemblies for three strains of Wolbachia pipientis, assembled from unenriched, unfiltered long-read shotgun sequencing data of geographically distinct strains of Drosophila melanogaster. Our simple methodology can be applied to long-read data sets of other Wolbachia-infected species with limited Wolbachia-host lateral gene transfers to produce complete assemblies for this important model symbiont.

Here, we report genome assemblies for three strains of Wolbachia pipientis, assembled from unenriched, unfiltered long-read shotgun sequencing data of geographically distinct strains of Drosophila melanogaster. Our simple methodology can be applied to long-read data sets of other Wolbachia-infected species with limited Wolbachia-host lateral gene transfers to produce complete assemblies for this important model symbiont.

Is Anopheles gambiae a Natural Host of Wolbachia?

Wolbachia (Alphaproteobacteria, Rickettsiales) is an intraovarially transmitted symbiont of insects able to exert striking phenotypes, including reproductive manipulations and pathogen blocking. These phenotypes make Wolbachia a promising tool to combat mosquito-borne diseases. Although Wolbachia is present in the majority of terrestrial arthropods, including many disease vectors, it was considered absent from Anopheles gambiae mosquitos, the main vectors of malaria in sub-Saharan Africa. In 2014, Wolbachia sequences were detected in A. gambiae samples collected in Burkina Faso. Subsequently, similar evidence came from collections all over Africa, revealing a high Wolbachia 16S rRNA sequence diversity, low abundance, and a lack of congruence between host and symbiont phylogenies. Here, we reanalyze and discuss recent evidence on the presence of Wolbachia sequences in A. gambiae. We find that although detected at increasing frequencies, the unusual properties of these Wolbachia sequences render them insufficient to diagnose natural infections in A. gambiae Future studies should focus on uncovering the origin of Wolbachia sequence variants in Anopheles and seeking sequence-independent evidence for this new symbiosis. Understanding the ecology of Anopheles mosquitos and their interactions with Wolbachia will be key in designing successful, integrative approaches to limit malaria spread. Although the prospect of using Wolbachia to fight malaria is intriguing, the newly discovered strains do not bring it closer to realization.IMPORTANCEAnopheles gambiae mosquitos are the main vectors of malaria, threatening around half of the world's population. The bacterial symbiont Wolbachia can interfere with disease transmission by other important insect vectors, but until recently, it was thought to be absent from natural A. gambiae populations. Here, we critically analyze the genomic, metagenomic, PCR, imaging, and phenotypic data presented in support of the presence of natural Wolbachia infections in A. gambiae We find that they are insufficient to diagnose Wolbachia infections and argue for the need of obtaining robust data confirming basic Wolbachia characteristics in this system. Determining the Wolbachia infection status of Anopheles is critical due to its potential to influence Anopheles population structure and Plasmodium transmission.

Combined effects of elevated CO2 concentration and Wolbachia on Hylyphantes graminicola (Araneae: Linyphiidae)

The increasing concentration of carbon dioxide in atmosphere is not only a major cause of global warming, but it also adversely affects the ecological diversity of invertebrates. This study was conducted to evaluate the effect of elevated CO2 concentration (ambient, 400 ppm and high, 800 ppm) and Wolbachia (Wolbachia-infected, W+ and Wolbachia-uninfected, W-) on Hylyphantes graminicola. The total survival rate, developmental duration, carapace width and length, body weight, sex ratio, net reproductive rate, nutrition content, and enzyme activity in H. graminicola were examined under four treatments: W- 400 ppm, W- 800 ppm, W+ 400 ppm, and W+ 800 ppm. Results showed that Wolbachia-infected spiders had significantly decreased the total developmental duration. Different instars showed variations up to some extent, but no obvious effect was found under elevated CO2 concentration. Total survival rate, sex ratio, and net reproductive rate were not affected by elevated CO2 concentration or Wolbachia infection. The carapace width of Wolbachia-uninfected spiders decreased significantly under elevated CO2 concentration, while the width, length and weight were not significantly affected in Wolbachia-infected spiders reared at ambient CO2 concentration. The levels of protein, specific activities of peroxidase, and amylase were significantly increased under elevated CO2 concentration or Wolbachia-infected spiders, while the total amino content was only increased in Wolbachia-infected spiders. Thus, our current finding suggested that elevated CO2 concentration and Wolbachia enhance nutrient contents and enzyme activity of H. graminicola and decrease development duration hence explore the interactive effects of factors which were responsible for reproduction regulation, but it also gives a theoretical direction for spider's protection in such a dynamic environment. Increased activities of enzymes and nutrients caused by Wolbachia infection aids for better survival of H. graminicola under stress.

Unexpected Diversity of Wolbachia Associated with Bactrocera dorsalis (Diptera: Tephritidae) in Africa

Bactrocera dorsalis (Hendel) is an important pest of fruit-bearing plants in many countries worldwide. In Africa, this pest has spread rapidly and has become widely established since the first invasion report in 2003. Wolbachia is a vertically transmitted endosymbiont that can significantly influence aspects of the biology and, in particular, the reproduction of its host. In this study, we screened B. dorsalis specimens collected from several locations in Africa between 2005 and 2017 for Wolbachia using a PCR-based assay to target the Wolbachia surface protein wsp. Of the 357 individuals tested, 10 were positive for Wolbachia using the wsp assay. We identified four strains of Wolbachia infecting two B. dorsalis mitochondrial haplotypes. We found no strict association between the infecting strain and host haplotype, with one strain being present in two different host haplotypes. All the detected strains belonged to Super Group B Wolbachia and did not match any strains reported previously in B. dorsalis in Asia. These findings indicate that diverse Wolbachia infections are present in invasive populations of B. dorsalis.

Genomic Advances in Avian Malaria Research

Haemosporidian parasites causing malaria-like diseases in birds are globally distributed and have been associated with reduced host fitness and mortality in susceptible bird species. This group of parasites has not only enabled a greater understanding of host specificity, virulence, and parasite dispersal, but has also been crucial in restructuring the evolutionary history of apicomplexans. Despite their importance, genomic resources of avian haemosporidians have proved difficult to obtain, and they have, as a result, been lagging behind the congeneric Plasmodium species infecting mammals. In this review, I discuss recent genomic advances in the field of avian malaria research, and outline outstanding questions that will become possible to investigate with the continued successful efforts to generate avian haemosporidian genomic data.