Assessing Cat Flea Microbiomes in Northern and Southern California by 16S rRNA Next-Generation Sequencing

Unique microbial communities of parasitic fleas on wild animals from the Qinghai-Tibet Plateau

Fleas, one of the most significant ectoparasites, play a crucial role as vectors in spreading zoonotic diseases globally. The Qinghai Province, as part of the Qinghai-Tibet Plateau, is one of the provinces in China with the largest number of flea species. In this study, we characterized the microbial communities of eighty-five adult fleas, belonging to nineteen species within four families (Ceratophyllidae, Ctenophthalmidae, Leptopsyllidae, and Pulicidae). We identified a total of 1162 unique operational taxonomic units at the genus level, with flea-borne pathogens such as Wolbachia, Bartonella, Rickettsia being the members of top abundant taxa. Except for comparison between Ctenophthalmidae and Leptopsyllidae families, the analyses of both alpha- and beta- diversity indicators suggested that bacterial diversity varied among flea families. This could be attributed to flea phylogeny, which also influenced by their geographical sites and animal hosts. Results of Linear discriminant analysis effect size (LEfSe) indicated that 29 genera in Ceratophylloidea, 11 genera in Ctenophthalmidae, 15 genera in Leptopsyllidae, and 22 genera in Pulicidae were significantly responsible for explaining the differences among the four flea families (linear discriminant analysis score > 2, P < 0.05). Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt2) analyses showed that the functional pathways varied significantly across flea families, which was supported by the significant correlation between the functional pathways and the microbial communities.

Cat Flea Coinfection with Rickettsia felis and Rickettsia typhi

Purpose: Flea-borne rickettsioses, collectively referred to as a term for etiological agents Rickettsia felis, Rickettsia typhi, and RFLOs (R. felis-like organisms), has become a public health concern around the world, specifically in the United States. Due to a shared arthropod vector (the cat flea) and clinical signs, discriminating between Rickettsia species has proven difficult. While the effects of microbial coinfections in the vector can result in antagonistic or synergistic interrelationships, subsequently altering potential human exposure and disease, the impact of bacterial interactions within flea populations remains poorly defined. Methods: In this study, in vitro and in vivo systems were utilized to assess rickettsial interactions in arthropods. Results: Coinfection of both R. felis and R. typhi within a tick-derived cell line indicated that the two species could infect the same cell, but distinct growth kinetics led to reduced R. felis growth over time, regardless of infection order. Sequential flea coinfections revealed the vector could acquire both Rickettsia spp. and sustain coinfection for up to 2 weeks, but rickettsial loads in coinfected fleas and feces were altered during coinfection. Conclusion: Altered rickettsial loads during coinfection suggest R. felis and R. typhi interactions may enhance the transmission potential of either agent. Thus, this study provides a functional foundation to disentangle transmission events propelled by complex interspecies relationships during vector coinfections.

Vector biology of the cat flea Ctenocephalides felis

Ctenocephalides felis, the cat flea, is among the most prevalent and widely dispersed vectors worldwide. Unfortunately, research on C. felis and associated pathogens (Bartonella and Rickettsia spp.) lags behind that of other vectors and vector-borne pathogens. Therefore, we aimed to review fundamental aspects of C. felis as a vector (behavior, epidemiology, phylogenetics, immunology, and microbiome composition) with an emphasis on key techniques and research avenues employed in other vector species. Future laboratory C. felis experimental infections with Bartonella, Rickettsia, and Wolbachia species/strains should examine the vector-pathogen interface utilizing contemporary visualization, transcriptomic, and gene-editing techniques. Further environmental sampling will inform the range and prevalence of C. felis and associated pathogens, improving the accuracy of vector and pathogen modeling to improve infection/infestation risk assessment and diagnostic recommendations.

Documenting the microbiome diversity and distribution in selected fleas from South Africa with an emphasis on the cat flea, Ctenocephalides f. felis

The factors that influence parasite associated bacterial microbial diversity and the geographic distributions of bacteria are not fully understood. In an effort to gain a deeper understanding of the relationship between the bacterial diversity of Ctenocephalides fleas and host species and the external environment, we conducted a metagenetic analysis of 107 flea samples collected from 8 distinct sampling sites in South Africa. Pooled DNA samples mostly comprising of 2 or 3 individuals sampled from the same host, and belonging to the same genetic cluster, were sequenced using the Ion PGM™ Hi-Q™ Kit and the Ion 316™ Chip v2. Differences were detected in the microbiome compositions between Ctenocephalides felis, Ctenocephalides canis and Ctenocephalides connatus. Although based on a small sample, C. connatus occurring on wildlife harboured a higher bacterial richness when compared to C. felis on domestic animals. Intraspecific differences in the microbial OTU diversity were detected within C. f. felis that occurred on domestic cats and dogs. Different genetic lineages of C. f. felis were similar in microbial compositions but some differences exist in the presence or absence of rare bacteria. Rickettsia and Bartonella OTU's identified in South African cat fleas differ from those identified in the USA and Australia. Intraspecific microbial compositions also differ across geographic sampling sites. Generalized dissimilarity modelling showed that temperature and humidity are potentially important environmental factors explaining the pattern obtained.

Complete De Novo Assembly of Wolbachia Endosymbiont of Frankliniella intonsa

As an endosymbiont, Wolbachia exerts significant effects on the host, including on reproduction, immunity, and metabolism. However, the study of Wolbachia in Thysanopteran insects, such as flower thrips Frankliniella intonsa, remains limited. Here, we assembled a gap-free looped genome assembly of Wolbachia strain wFI in a length of 1,463,884 bp (GC content 33.80%), using Nanopore long reads and Illumina short reads. The annotation of wFI identified a total of 1838 protein-coding genes (including 85 pseudogenes), 3 ribosomal RNAs (rRNAs), 35 transfer RNAs (tRNAs), and 1 transfer-messenger RNA (tmRNA). Beyond this basic description, we identified mobile genetic elements, such as prophage and insertion sequences (ISs), which make up 17% of the entire wFI genome, as well as genes involved in riboflavin and biotin synthesis and metabolism. This research lays the foundation for understanding the nutritional mutualism between Wolbachia and flower thrips. It also serves as a valuable resource for future studies delving into the intricate interactions between Wolbachia and its host.

Metagenomics of culture isolates and insect tissue illuminate the evolution of Wolbachia, Rickettsia and Bartonella symbionts in Ctenocephalides spp. fleas

While fleas are often perceived simply as a biting nuisance and a cause of allergic dermatitis, they represent important disease vectors worldwide, especially for bacterial zoonoses such as plague (transmitted by rodent fleas) and some of the rickettsioses and bartonelloses. The cosmopolitan cat (Ctenocephalides felis ) and dog (Ctenocephalides canis ) fleas, as well as Ctenocephalides orientis (restricted to tropical and subtropical Asia), breed in human dwellings and are vectors of cat-scratch fever (caused by Bartonella spp.) and Rickettsia spp., including Rickettsia felis (agent of flea-borne spotted fever) and Rickettsia asembonensis , a suspected pathogen. These Rickettsia spp. are members of a phylogenetic clade known as the ‘transitional group’, which includes both human pathogens and arthropod-specific endosymbionts. The relatively depauperate flea microbiome can also contain other endosymbionts, including a diverse range of Wolbachia strains. Here, we present circularized genome assemblies for two C. orientis -derived pathogens (Bartonella clarridgeiae and R. asembonensis ) from Malaysia, a novel Wolbachia strain (w Cori), and the C. orientis mitochondrion; all were obtained by direct metagenomic sequencing of flea tissues. Moreover, we isolated two Wolbachia strains from Malaysian C. felis into tick cell culture and recovered circularized genome assemblies for both, one of which (w CfeF) is newly sequenced. We demonstrate that the three Wolbachia strains are representatives of different major clades (‘supergroups’), two of which appear to be flea-specific. These Wolbachia genomes exhibit unique combinations of features associated with reproductive parasitism or mutualism, including prophage WO, cytoplasmic incompatibility factors and the biotin operon of obligate intracellular microbes. The first circularized assembly for R. asembonensis includes a plasmid with a markedly different structure and gene content compared to the published plasmid; moreover, this novel plasmid was also detected in cat flea metagenomes from the USA. Analysis of loci under positive selection in the transitional group revealed genes involved in host–pathogen interactions that may facilitate host switching. Finally, the first B. clarridgeiae genome from Asia exhibited large-scale genome stability compared to isolates from other continents, except for SNPs in regions predicted to mediate interactions with the vertebrate host. These findings highlight the paucity of data on the genomic diversity of Ctenocephalides -associated bacteria and raise questions regarding how interactions between members of the flea microbiome might influence vector competence.

Symbioses shape feeding niches and diversification across insects

For over 300 million years, insects have relied on symbiotic microbes for nutrition and defence. However, it is unclear whether specific ecological conditions have repeatedly favoured the evolution of symbioses, and how this has influenced insect diversification. Here, using data on 1,850 microbe-insect symbioses across 402 insect families, we found that symbionts have allowed insects to specialize on a range of nutrient-imbalanced diets, including phloem, blood and wood. Across diets, the only limiting nutrient consistently associated with the evolution of obligate symbiosis was B vitamins. The shift to new diets, facilitated by symbionts, had mixed consequences for insect diversification. In some cases, such as herbivory, it resulted in spectacular species proliferation. In other niches, such as strict blood feeding, diversification has been severely constrained. Symbioses therefore appear to solve widespread nutrient deficiencies for insects, but the consequences for insect diversification depend on the feeding niche that is invaded.

Complete bacterial profile and potential pathogens of cat fleas Ctenocephalides felis

Fleas are important ectoparasites and vectors associated with a wide range of pathogenic diseases, posing threats to public health concerns, especially cat fleas that spread worldwide. Understanding the microbial components is essential due to cat fleas are capable of transmitting pathogens to humans, causing diseases like plague and murine typhus. In the present study, metagenomic next-generation sequencing was applied to obtain the complete microbiota and related functions in the gut of Ctenocephalides felis. A total of 1,870 species was taxonomically recognized including 1,407 bacteria, 365 eukaryotes, 69 viruses, and 29 archaea. Proteobacteria was the dominant phylum among the six samples. Pathogens Rickettsia felis, Acinetobacter baumannii, Coxiella burnetii, and Anaplasma phagocytophilum were taxonomically identified and had high abundances in all samples. The resistance gene MexD was predominant in microbial communities of all cat fleas. We also performed epidemiological surveys of pathogens R. felis, A. baumannii, C. burnetii, and A. phagocytophilum among 165 cat fleas collected from seven provinces in China, while only the DNAs of R. felis (38/165, 23.03%) and C. burnetii (2/165, 1.21%) were obtained. The data provide new insight and understanding of flea intestinal microbiota and provided novel information for preventing and controlling fleas and their transmitted diseases.

Feeding on a Bartonella henselae Infected Host Triggers Temporary Changes in the Ctenocephalides felis Microbiome

The effect of Bartonella henselae on the microbiome of its vector, Ctenocephalides felis (the cat flea) is largely unknown, as the majority of C. felis microbiome studies have utilized wild-caught pooled fleas. We surveyed the microbiome of laboratory-origin C. felis fed on B. henselae-infected cats for 24 h or 9 days to identify changes to microbiome diversity and microbe prevalence compared to unfed fleas, and fleas fed on uninfected cats. Utilizing Next Generation Sequencing (NGS) on the Illumina platform, we documented an increase in microbial diversity in C. felis fed on Bartonella-infected cats for 24 h. These changes returned to baseline (unfed fleas or fleas fed on uninfected cats) after 9 days on the host. Increased diversity in the C. felis microbiome when fed on B. henselae-infected cats may be related to the mammalian, flea, or endosymbiont response. Poor B. henselae acquisition was documented with only one of four infected flea pools having B. henselae detected by NGS. We hypothesize this is due to the use of adult fleas, flea genetic variation, or lack of co-feeding with B. henselae-infected fleas. Future studies are necessary to fully characterize the effect of endosymbionts and C. felis diversity on B. henselae acquisition.

Identification of microbial taxa present in Ctenocephalides felis (cat flea) reveals widespread co-infection and associations with vector phylogeny

BACKGROUND

Ctenocephalides felis, the cat flea, is the most common ectoparasite of cats and dogs worldwide. As a cause of flea allergy dermatitis and a vector for two genera of zoonotic pathogens (Bartonella and Rickettsia spp.), the effect of the C. felis microbiome on pathogen transmission and vector survival is of substantial medical importance to both human and veterinary medicine. The aim of this study was to assay the pathogenic and commensal eubacterial microbial communities of individual C. felis from multiple geographic locations and analyze these findings by location, qPCR pathogen prevalence, and flea genetic diversity.

METHODS

16S Next Generation Sequencing (NGS) was utilized to sequence the microbiome of fleas collected from free-roaming cats, and the cox1 gene was used for flea phylogenetic analysis. NGS data were analyzed for 168 individual fleas from seven locations within the US and UK. Given inconsistency in the genera historically reported to constitute the C. felis microbiome, we utilized the decontam prevalence method followed by literature review to separate contaminants from true microbiome members.

RESULTS

NGS identified a single dominant and cosmopolitan amplicon sequence variant (ASV) from Rickettsia and Wolbachia while identifying one dominant Bartonella clarridgeiae and one dominant Bartonella henselae/Bartonella koehlerae ASV. Multiple less common ASVs from these genera were detected within restricted geographical ranges. Co-detection of two or more genera (Bartonella, Rickettsia, and/or Wolbachia) or multiple ASVs from a single genus in a single flea was common. Achromobacter, Peptoniphilus, and Rhodococcus were identified as additional candidate members of the C. felis microbiome on the basis of decontam analysis and literature review. Ctenocephalides felis phylogenetic diversity as assessed by the cox1 gene fell within currently characterized clades while identifying seven novel haplotypes. NGS sensitivity and specificity for Bartonella and Rickettsia spp. DNA detection were compared to targeted qPCR.

CONCLUSIONS

Our findings confirm the widespread coinfection of fleas with multiple bacterial genera and strains, proposing three additional microbiome members. The presence of minor Bartonella, Rickettsia, and Wolbachia ASVs was found to vary by location and flea haplotype. These findings have important implications for flea-borne pathogen transmission and control.

Complete de novo assembly of Wolbachia endosymbiont of Diaphorina citri Kuwayama (Hemiptera: Liviidae) using long-read genome sequencing

Wolbachia, a gram-negative \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{\alpha }$$\end{document}α-proteobacterium, is an endosymbiont found in some arthropods and nematodes. Diaphorina citri Kuwayama, the vector of ‘Candidatus Liberibacter asiaticus’ (CLas), are naturally infected with a strain of Wolbachia (wDi), which has been shown to colocalize with the bacteria pathogens CLas, the pathogen associated with huanglongbing (HLB) disease of citrus. The relationship between wDi and CLas is poorly understood in part because the complete genome of wDi has not been available. Using high-quality long-read PacBio circular consensus sequences, we present the largest complete circular wDi genome among supergroup-B members. The assembled circular chromosome is 1.52 megabases with 95.7% genome completeness with contamination of 1.45%, as assessed by checkM. We identified Insertion Sequences (ISs) and prophage genes scattered throughout the genomes. The proteins were annotated using Pfam, eggNOG, and COG that assigned unique domains and functions. The wDi genome was compared with previously sequenced Wolbachia genomes using pangenome and phylogenetic analyses. The availability of a complete circular chromosome of wDi will facilitate understanding of its role within the insect vector, which may assist in developing tools for disease management. This information also provides a baseline for understanding phylogenetic relationships among Wolbachia of other insect vectors.

Data analysis workflow for the detection of canine vector-borne pathogens using 16 S rRNA Next-Generation Sequencing

BACKGROUND

Vector-borne diseases (VBDs) impact both human and veterinary medicine and pose special public health challenges. The main bacterial vector-borne pathogens (VBPs) of importance in veterinary medicine include Anaplasma spp., Bartonella spp., Ehrlichia spp., and Spotted Fever Group Rickettsia. Taxon-targeted PCR assays are the current gold standard for VBP diagnostics but limitations on the detection of genetically diverse organisms support a novel approach for broader detection of VBPs. We present a methodology for genetic characterization of VBPs using Next-Generation Sequencing (NGS) and computational approaches. A major advantage of NGS is the ability to detect multiple organisms present in the same clinical sample in an unsupervised (i.e. non-targeted) and semi-quantitative way. The Standard Operating Procedure (SOP) presented here combines industry-standard microbiome analysis tools with our ad-hoc bioinformatic scripts to form a complete analysis pipeline accessible to veterinary scientists and freely available for download and use at https://github.com/eltonjrv/microbiome.westernu/tree/SOP .

RESULTS

We tested and validated our SOP by mimicking single, double, and triple infections in genomic canine DNA using serial dilutions of plasmids containing the entire 16 S rRNA gene sequence of (A) phagocytophilum, (B) v. berkhoffii, and E. canis. NGS with broad-range 16 S rRNA primers followed by our bioinformatics SOP was capable of detecting these pathogens in biological replicates of different dilutions. These results illustrate the ability of NGS to detect and genetically characterize multi-infections with different amounts of pathogens in a single sample.

CONCLUSIONS

Bloodborne microbiomics & metagenomics approaches may help expand the molecular diagnostic toolbox in veterinary and human medicine. In this paper, we present both in vitro and in silico detailed protocols that can be combined into a single workflow that may provide a significant improvement in VBP diagnostics and also facilitate future applications of microbiome research in veterinary medicine.

Changes in Bacterial Diversity, Composition and Interactions During the Development of the Seabird Tick Ornithodoros maritimus (Argasidae)

Characterising within-host microbial interactions is essential to understand the drivers that shape these interactions and their consequences for host ecology and evolution. Here, we examined the bacterial microbiota hosted by the seabird soft tick Ornithodoros maritimus (Argasidae) in order to uncover bacterial interactions within ticks and how these interactions change over tick development. Bacterial communities were characterised through next-generation sequencing of the V3-V4 hypervariable region of the bacterial 16S ribosomal RNA gene. Bacterial co-occurrence and co-exclusion were determined by analysing networks generated from the metagenomic data obtained at each life stage. Overall, the microbiota of O. maritimus was dominated by four bacterial genera, namely Coxiella, Rickettsia, Brevibacterium and Arsenophonus, representing almost 60% of the reads. Bacterial diversity increased over tick development, and adult male ticks showed higher diversity than did adult female ticks. Bacterial networks showed that co-occurrence was more frequent than co-exclusion and highlighted substantial shifts across tick life stages; interaction networks changed from one stage to the next with a steady increase in the number of interactions through development. Although many bacterial interactions appeared unstable across life stages, some were maintained throughout development and were found in both sexes, such as Coxiella and Arsenophonus. Our data support the existence of a few stable interactions in O. maritimus ticks, on top of which bacterial taxa accumulate from hosts and/or the environment during development. We propose that stable associations delineate core microbial interactions, which are likely to be responsible for key biological functions.

Molecular Evidence of a Broad Range of Pathogenic Bacteria in Ctenocephalides spp.: Should We Re-Examine the Role of Fleas in the Transmission of Pathogens?

The internal microbiome of common cat and dog fleas was studied for DNA evidence of pathogenic bacteria. Fleas were grouped in pools by parasitized animal. DNA was extracted and investigated with 16S metagenomics for medically relevant (MR) bacteria, based on the definitions of the International Statistical Classification of Diseases and Related Health Problems (WHO). The MR bacterial species totaled 40, were found in 60% of flea-pools (N = 100), and included Acinetobacterbaumannii, Bacteroidesfragilis, Clostridiumperfringens, Enterococcusfaecalis, E. mundtii, Fusobacteriumnucleatum, Haemophilusaegyptius, Kingellakingae, Klebsiellapneumoniae, Leptotrichiabuccalis, L. hofstadii, Moraxellalacunata, Pasteurellamultocida, Propionibacteriumacnes, P. propionicum, Proteusmirabilis, Pseudomonasaeruginosa, Rickettsiaaustralis, R. hoogstraalii, Salmonellaenterica, and various Bartonella, Staphylococcus, and Streptococcus species. B. henselae (p = 0.004) and B. clarridgeiae (p = 0.006) occurred more frequently in fleas from cats, whereas Rickettsiahoogstraalii (p = 0.031) and Propionibacteriumacnes (p = 0.029) had a preference in fleas from stray animals. Most of the discovered MR species can form biofilm, and human exposure may theoretically occur through the flea-host interface. The fitness of these pathogenic bacteria to cause infection and the potential role of fleas in the transmission of a broad range of diseases should be further investigated.

Rigorous quantification of bacterial richness in ticks: A case study

Kwan et al. (2017) published an informative study comparing results obtained by next-generation sequencing (NGS) of mean bacterial genera richness among different life stages, male and female adults, and rearing conditions (field vs. laboratory) for Ixodes pacificus. The current paper examines Kwan et al. (2017) as a case study to provide guidance on statistical design and analysis for estimation of richness, derived from next generation sequencing technology, of the bacterial microbiome in field-collected I. pacificus. Suggestions are provided to further strengthen quantification of microbiome richness in studies in ticks, with focus on sampling design. In-depth treatment is provided of the relative merits of estimating mean richness versus median richness. Research on microbiome diversity in ticks can be made quantitatively rigorous; although, more research on methods is needed.

Evolution of Wolbachia mutualism and reproductive parasitism: insight from two novel strains that co-infect cat fleas

Wolbachiae are obligate intracellular bacteria that infect arthropods and certain nematodes. Usually maternally inherited, they may provision nutrients to (mutualism) or alter sexual biology of (reproductive parasitism) their invertebrate hosts. We report the assembly of closed genomes for two novel wolbachiae, wCfeT and wCfeJ, found co-infecting cat fleas (Ctenocephalides felis) of the Elward Laboratory colony (Soquel, CA, USA). wCfeT is basal to nearly all described Wolbachia supergroups, while wCfeJ is related to supergroups C, D and F. Both genomes contain laterally transferred genes that inform on the evolution of Wolbachia host associations. wCfeT carries the Biotin synthesis Operon of Obligate intracellular Microbes (BOOM); our analyses reveal five independent acquisitions of BOOM across the Wolbachia tree, indicating parallel evolution towards mutualism. Alternately, wCfeJ harbors a toxin-antidote operon analogous to the wPip cinAB operon recently characterized as an inducer of cytoplasmic incompatibility (CI) in flies. wCfeJ cinB and three adjacent genes are collectively similar to large modular toxins encoded in CI-like operons of certain Wolbachia strains and Rickettsia species, signifying that CI toxins streamline by fission of large modular toxins. Remarkably, the C. felis genome itself contains two CI-like antidote genes, divergent from wCfeJ cinA, revealing episodic reproductive parasitism in cat fleas and evidencing mobility of CI loci independent of WO-phage. Additional screening revealed predominant co-infection (wCfeT/wCfeJ) amongst C. felis colonies, though fleas in wild populations mostly harbor wCfeT alone. Collectively, genomes of wCfeT, wCfeJ, and their cat flea host supply instances of lateral gene transfers that could drive transitions between parasitism and mutualism.

A chromosome-level assembly of the cat flea genome uncovers rampant gene duplication and genome size plasticity

BACKGROUND

Fleas (Insecta: Siphonaptera) are small flightless parasites of birds and mammals; their blood-feeding can transmit many serious pathogens (i.e., the etiological agents of bubonic plague, endemic and murine typhus). The lack of flea genome assemblies has hindered research, especially comparisons to other disease vectors. Accordingly, we sequenced the genome of the cat flea, Ctenocephalides felis, an insect with substantial human health and veterinary importance across the globe.

RESULTS

By combining Illumina and PacBio sequencing of DNA derived from multiple inbred female fleas with Hi-C scaffolding techniques, we generated a chromosome-level genome assembly for C. felis. Unexpectedly, our assembly revealed extensive gene duplication across the entire genome, exemplified by ~ 38% of protein-coding genes with two or more copies and over 4000 tRNA genes. A broad range of genome size determinations (433-551 Mb) for individual fleas sampled across different populations supports the widespread presence of fluctuating copy number variation (CNV) in C. felis. Similarly, broad genome sizes were also calculated for individuals of Xenopsylla cheopis (Oriental rat flea), indicating that this remarkable "genome-in-flux" phenomenon could be a siphonapteran-wide trait. Finally, from the C. felis sequence reads, we also generated closed genomes for two novel strains of Wolbachia, one parasitic and one symbiotic, found to co-infect individual fleas.

CONCLUSION

Rampant CNV in C. felis has dire implications for gene-targeting pest control measures and stands to complicate standard normalization procedures utilized in comparative transcriptomics analysis. Coupled with co-infection by novel Wolbachia endosymbionts-potential tools for blocking pathogen transmission-these oddities highlight a unique and underappreciated disease vector.