Legionella Becoming a Mutualist: Adaptive Processes Shaping the Genome of Symbiont in the Louse Polyplax serrata

Discovery of a novel symbiotic lineage associated with a hematophagous leech from the genus Haementeria

Similarly to other strict blood feeders, leeches from the Haementeria genus (Hirudinida: Glossiphoniidae) have established a symbiotic association with bacteria harbored intracellularly in esophageal bacteriomes. Previous genome sequence analyses of these endosymbionts revealed co-divergence with their hosts, a strong genome reduction, and a simplified metabolism largely dedicated to the production of B vitamins, which are nutrients lacking from a blood diet. 'Candidatus Providencia siddallii' has been identified as the obligate nutritional endosymbiont of a monophyletic clade of Mexican and South American Haementeria spp. However, the Haementeria genus includes a sister clade of congeners from Central and South America, where the presence or absence of the aforementioned symbiont taxon remains unknown. In this work, we report on a novel bacterial endosymbiont found in a representative from this Haementeria clade. We found that this symbiont lineage has evolved from within the Pluralibacter genus, known mainly from clinical but also environmental strains. Similarly to Ca. Providencia siddallii, the Haementeria-associated Pluralibacter symbiont displays clear signs of genome reduction, accompanied by an A+T-biased sequence composition. Genomic analysis of its metabolic potential revealed a retention of pathways related to B vitamin biosynthesis, supporting its role as a nutritional endosymbiont. Finally, comparative genomics of both Haementeria symbiont lineages suggests that an ancient Providencia symbiont was likely replaced by the novel Pluralibacter one, thus constituting the first reported case of nutritional symbiont replacement in a leech without morphological changes in the bacteriome.

IMPORTANCE

Obligate symbiotic associations with a nutritional base have likely evolved more than once in strict blood-feeding leeches. Unlike those symbioses found in hematophagous arthropods, the nature, identity, and evolutionary history of these remains poorly studied. In this work, we further explored obligate nutritional associations between Haementeria leeches and their microbial symbionts, which led to the unexpected discovery of a novel symbiosis with a member of the Pluralibacter genus. When compared to Providencia siddallii, an obligate nutritional symbiont of other Haementeria leeches, this novel bacterial symbiont shows convergent retention of the metabolic pathways involved in B vitamin biosynthesis. Moreover, the genomic characteristics of this Pluralibacter symbiont suggest a more recent association than that of Pr. siddallii and Haementeria. We conclude that the once-thought stable associations between blood-feeding Glossiphoniidae and their symbionts (i.e., one bacteriome structure, one symbiont lineage) can break down, mirroring symbiont turnover observed in various arthropod lineages.

Highly Resolved Genomes of Two Closely Related Lineages of the Rodent Louse Polyplax serrata with Different Host Specificities

Sucking lice of the parvorder Anoplura are permanent ectoparasites with specific lifestyle and highly derived features. Currently, genomic data are only available for a single species, the human louse Pediculus humanus. Here, we present genomes of two distinct lineages, with different host spectra, of a rodent louse Polyplax serrata. Genomes of these ecologically different lineages are closely similar in gene content and display a conserved order of genes, with the exception of a single translocation. Compared with P. humanus, the P. serrata genomes are noticeably larger (139 vs. 111 Mbp) and encode a higher number of genes. Similar to P. humanus, they are reduced in sensory-related categories such as vision and olfaction. Utilizing genome-wide data, we perform phylogenetic reconstruction and evolutionary dating of the P. serrata lineages. Obtained estimates reveal their relatively deep divergence (∼6.5 Mya), comparable with the split between the human and chimpanzee lice P. humanus and Pediculus schaeffi. This supports the view that the P. serrata lineages are likely to represent two cryptic species with different host spectra. Historical demographies show glaciation-related population size (Ne) reduction, but recent restoration of Ne was seen only in the less host-specific lineage. Together with the louse genomes, we analyze genomes of their bacterial symbiont Legionella polyplacis and evaluate their potential complementarity in synthesis of amino acids and B vitamins. We show that both systems, Polyplax/Legionella and Pediculus/Riesia, display almost identical patterns, with symbionts involved in synthesis of B vitamins but not amino acids.

Arsenophonus symbiosis with louse flies: multiple origins, coevolutionary dynamics, and metabolic significance

IMPORTANCE

Insects that live exclusively on vertebrate blood utilize symbiotic bacteria as a source of essential compounds, e.g., B vitamins. In louse flies, the most frequent symbiont originated in genus Arsenophonus, known from a wide range of insects. Here, we analyze genomic traits, phylogenetic origins, and metabolic capacities of 11 Arsenophonus strains associated with louse flies. We show that in louse flies, Arsenophonus established symbiosis in at least four independent events, reaching different stages of symbiogenesis. This allowed for comparative genomic analysis, including convergence of metabolic capacities. The significance of the results is twofold. First, based on a comparison of independently originated Arsenophonus symbioses, it determines the importance of individual B vitamins for the insect host. This expands our theoretical insight into insect-bacteria symbiosis. The second outcome is of methodological significance. We show that the comparative approach reveals artifacts that would be difficult to identify based on a single-genome analysis.

Microbiomes of Blood-Feeding Triatomines in the Context of Their Predatory Relatives and the Environment

The importance of gut microbiomes has become generally recognized in vector biology. This study addresses microbiome signatures in North American Triatoma species of public health significance (vectors of Trypanosoma cruzi) linked to their blood-feeding strategy and the natural habitat. To place the Triatoma-associated microbiomes within a complex evolutionary and ecological context, we sampled sympatric Triatoma populations, related predatory reduviids, unrelated ticks, and environmental material from vertebrate nests where these arthropods reside. Along with five Triatoma species, we have characterized microbiomes of five reduviids (Stenolemoides arizonensis, Ploiaria hirticornis, Zelus longipes, and two Reduvius species), a single soft tick species, Ornithodoros turicata, and environmental microbiomes from selected sites in Arizona, Texas, Florida, and Georgia. The microbiomes of predatory reduviids lack a shared core microbiota. As in triatomines, microbiome dissimilarities among species correlate with dominance of a single bacterial taxon. These include Rickettsia, Lactobacillus, "Candidatus Midichloria," and Zymobacter, which are often accompanied by known symbiotic genera, i.e., Wolbachia, "Candidatus Lariskella," Asaia, Gilliamella, and Burkholderia. We have further identified a compositional convergence of the analyzed microbiomes in regard to the host phylogenetic distance in both blood-feeding and predatory reduviids. While the microbiomes of the two reduviid species from the Emesinae family reflect their close relationship, the microbiomes of all Triatoma species repeatedly form a distinct monophyletic cluster highlighting their phylosymbiosis. Furthermore, based on environmental microbiome profiles and blood meal analysis, we propose three epidemiologically relevant and mutually interrelated bacterial sources for Triatoma microbiomes, i.e., host abiotic environment, host skin microbiome, and pathogens circulating in host blood. IMPORTANCE This study places microbiomes of blood-feeding North American Triatoma vectors (Reduviidae) into a broader evolutionary and ecological context provided by related predatory assassin bugs (Reduviidae), another unrelated vector species (soft tick Ornithodoros turicata), and the environment these arthropods coinhabit. For both vectors, microbiome analyses suggest three interrelated sources of bacteria, i.e., the microbiome of vertebrate nests as their natural habitat, the vertebrate skin microbiome, and the pathobiome circulating in vertebrate blood. Despite an apparent influx of environment-associated bacteria into the arthropod microbiomes, Triatoma microbiomes retain their specificity, forming a distinct cluster that significantly differs from both predatory relatives and ecologically comparable ticks. Similarly, within the related predatory Reduviidae, we found the host phylogenetic distance to underlie microbiome similarities.

Supergroup F Wolbachia with extremely reduced genome: transition to obligate insect symbionts

BACKGROUND

Wolbachia belong to highly abundant bacteria which are frequently found in invertebrate microbiomes and manifest by a broad spectrum of lifestyles from parasitism to mutualism. Wolbachia supergroup F is a particularly interesting clade as it gave rise to symbionts of both arthropods and nematodes, and some of its members are obligate mutualists. Investigations on evolutionary transitions among the different symbiotic stages have been hampered by a lack of the known diversity and genomic data for the supergroup F members.

RESULTS

Based on amplicon screening, short- and long-read WGS approaches, and laser confocal microscopy, we characterize five new supergroup F Wolbachia strains from four chewing lice species. These strains reached different evolutionary stages and represent two remarkably different types of symbiont genomes. Three of the genomes resemble other known members of Wolbachia F supergroup, while the other two show typical signs of ongoing gene inactivation and removal (genome size, coding density, low number of pseudogenes). Particularly, wMeur1, a symbiont fixed in microbiomes of Menacanthus eurysternus across four continents, possesses a highly reduced genome of 733,850 bp. The horizontally acquired capacity for pantothenate synthesis and localization in specialized bacteriocytes suggest its obligate nutritional role.

CONCLUSIONS

The genome of wMeur1 strain, from the M. eurysternus microbiome, represents the smallest currently known Wolbachia genome and the first example of Wolbachia which has completed genomic streamlining as known from the typical obligate symbionts. This points out that despite the large amount and great diversity of the known Wolbachia strains, evolutionary potential of these bacteria still remains underexplored. The diversity of the four chewing lice microbiomes indicates that this vast parasitic group may provide suitable models for further investigations. Video Abstract.

Genomic Approaches to Uncovering the Coevolutionary History of Parasitic Lice

Simple Summary

New sequencing technologies have now made it possible to sequence entire genomes for a diversity of life on earth. Parasites comprise nearly half of all species. Lice are one important group of parasites of birds and mammals, including humans. Genome sequencing approaches have been applied to this group of parasites to uncover patterns of diversification. These patterns can be compared to the patterns of diversification in their hosts. Key findings from these studies have revealed that parasitic lice likely originated on birds and then switched to mammals multiple times. Within groups of birds and mammals, the evolutionary trees of lice match those for mammal hosts more than those for birds. Genomic approaches have also revealed that individual birds and mammals harbor distinct populations of lice. Thus, these new techniques allow for the study of patterns of diversification at a wide variety of scales.

Abstract

Next-generation sequencing technologies are revolutionizing the fields of genomics, phylogenetics, and population genetics. These new genomic approaches have been extensively applied to a major group of parasites, the lice (Insecta: Phthiraptera) of birds and mammals. Two louse genomes have been assembled and annotated to date, and these have opened up new resources for the study of louse biology. Whole genome sequencing has been used to assemble large phylogenomic datasets for lice, incorporating sequences of thousands of genes. These datasets have provided highly supported trees at all taxonomic levels, ranging from relationships among the major groups of lice to those among closely related species. Such approaches have also been applied at the population scale in lice, revealing patterns of population subdivision and inbreeding. Finally, whole genome sequence datasets can also be used for additional study beyond that of the louse nuclear genome, such as in the study of mitochondrial genome fragmentation or endosymbiont function.

Endosymbiotic bacteria of the boar louse Haematopinus apri (Insecta: Phthiraptera: Anoplura)

Insects exclusively feeding on vertebrate blood are usually dependent on symbiotic bacteria for provisioning of B vitamins. Among them, sucking lice are prominent in that their symbiotic bacteria as well as their symbiotic organs exhibit striking diversity. Here we investigated the bacterial diversity associated with the boar louse Haematopinus apri in comparison with the hog louse Haematopinus suis. Amplicon sequencing analysis identified the primary endosymbiont predominantly detected from all populations of H. apri with some minor secondary bacterial associates. Sequencing and phylogenetic analysis of bacterial 16S rRNA gene confirmed that the endosymbionts of the boar louse H. apri, the hog louse H. suis and the cattle louse Haematopinus eurysternus form a distinct clade in the Gammaproteobacteria. The endosymbiont clade of Haematopinus spp. was phylogenetically distinct from the primary endosymbionts of other louse lineages. Fluorescence in situ hybridization visualized the endosymbiont localization within midgut epithelium, ovarial ampulla and posterior oocyte of H. apri, which were substantially the same as the endosymbiont localization previously described in H. suis and H. eurysternus. Mitochondrial haplotype analysis revealed that, although the domestic pig was derived from the wild boar over the past 8,000 years of human history, the populations of H. apri constituted a distinct sister clade to the populations of H. suis. Based on these results, we discussed possible evolutionary trajectories of the boar louse, the hog louse and their endosymbionts in the context of swine domestication. We proposed ‘Candidatus Haematopinicola symbiotica’ for the distinct clade of the endosymbionts of Haematopinus spp.

Lightella neohaematopini: A new lineage of highly reduced endosymbionts coevolving with chipmunk lice of the genus Neohaematopinus

Sucking lice (Anoplura) are known to have established symbiotic associations multiple times with different groups of bacteria as diverse as Enterobacteriales, Legionellales, and Neisseriales. This diversity, together with absence of a common coevolving symbiont (such as Buchnera, in aphids), indicates that sucking lice underwent a series of symbiont acquisitions, losses, and replacements. To better understand evolution and significance of louse symbionts, genomic and phylogenetic data are needed from a broader taxonomic diversity of lice and their symbiotic bacteria. In this study, we extend the known spectrum of the louse symbionts with a new lineage associated with Neohaematopinus pacificus, a louse species that commonly parasitizes North American chipmunks. The recent coevolutionary analysis showed that rather than a single species, these lice form a cluster of unique phylogenetic lineages specific to separate chipmunk species (or group of closely related species). Using metagenomic assemblies, we show that the lice harbor a bacterium which mirrors their phylogeny and displays traits typical for obligate mutualists. Phylogenetic analyses place this bacterium within Enterobacteriaceae on a long branch related to another louse symbiont, “Candidatus Puchtella pedicinophila.” We propose for this symbiotic lineage the name “Candidatus Lightella neohaematopini.” Based on the reconstruction of metabolic pathways, we suggest that like other louse symbionts, L. neohaematopini provides its host with at least some B vitamins. In addition, several samples harbored another symbiotic bacterium phylogenetically affiliated with the Neisseriales-related symbionts described previously from the lice Polyplax serrata and Hoplopleura acanthopus. Characterizing these bacteria further extend the known diversity of the symbiotic associations in lice and show unique complexity and dynamics of the system.

Development and validation of an LC-MS/MS method for determination of B vitamins and some its derivatives in whole blood

Obligate symbiotic bacteria associated with the insects feeding exclusively on vertebrate blood are supposed to complement B vitamins presumably lacking in their diet. Recent genomic analyses revealed considerable differences in biosynthetic capacities across different symbionts, suggesting that levels of B vitamins may vary across different vertebrate hosts. However, a rigorous determination of B vitamins content in blood of various vertebrates has not yet been approached. A reliable analytical method focused on B vitamin complex in blood can provide valuable informative background and understanding of general principles of insect symbiosis. In this work, a chromatographic separation of eight B vitamins (thiamine, riboflavin, niacin, pantothenic acid, pyridoxine, biotin, folic acid, and cyanocobalamine), four B vitamin derivatives (niacinamide, pyridoxal-5-phosphate, 4-pyridoxic acid, and tetrahydrofolic acid), and 3 stable isotope labelled internal standards was developed. Detection was carried out using dual-pressure linear ion trap mass spectrometer in FullScan MS/MS and SIM mode. Except for vitamin B9 (tetrahydrofolic acid), the instrument quantitation limits of all analytes were ranging from 0.42 to 5.0 μg/L, correlation coefficients from 0.9997 to 1.0000, and QC coefficients from 0.53 to 3.2%. Optimization of whole blood sample preparation step was focused especially on evaluation of two types of protein-precipitation agents: trichloroacetic acid and zinc sulphate in methanol. The best results were obtained for zinc sulphate in methanol, but only nine analytes were successfully validated. Accuracy of the procedure using this protein-precipitating agent was ranging from 89 to 120%, precision from 0.5 to 13%, and process efficiency from 65 to 108%. The content of B vitamins in whole blood samples from human and various vertebrates is presented as an application example of this newly developed method.

Cryptic Genes for Interbacterial Antagonism Distinguish Rickettsia Species Infecting Blacklegged Ticks From Other Rickettsia Pathogens

Background

The genus Rickettsia (Alphaproteobacteria: Rickettsiales) encompasses numerous obligate intracellular species with predominantly ciliate and arthropod hosts. Notable species are pathogens transmitted to mammals by blood-feeding arthropods. Mammalian pathogenicity evolved from basal, non-pathogenic host-associations; however, some non-pathogens are closely related to pathogens. One such species, Rickettsia buchneri, is prevalent in the blacklegged tick, Ixodes scapularis. While I. scapularis transmits several pathogens to humans, it does not transmit Rickettsia pathogens. We hypothesize that R. buchneri established a mutualism with I. scapularis, blocking tick superinfection with Rickettsia pathogens.

Methods

To improve estimates for assessing R. buchneri infection frequency in blacklegged tick populations, we used comparative genomics to identify an R. buchneri gene (REIS_1424) not present in other Rickettsia species present throughout the I. scapularis geographic range. Bioinformatic and phylogenomics approaches were employed to propose a function for the hypothetical protein (263 aa) encoded by REIS_1424.

Results

REIS_1424 has few analogs in other Rickettsiales genomes and greatest similarity to non-Proteobacteria proteins. This cohort of proteins varies greatly in size and domain composition, possessing characteristics of Recombination hotspot (Rhs) and contact dependent growth inhibition (CDI) toxins, with similarity limited to proximal C-termini (~145 aa). This domain was named CDI-like/Rhs-like C-terminal toxin (CRCT). As such proteins are often found as toxin-antidote (TA) modules, we interrogated REIS_1423 (151 aa) as a putative antidote. Indeed, REIS_1423 is similar to proteins encoded upstream of CRCT domain-containing proteins. Accordingly, we named these proteins CDI-like/Rhs-like C-terminal toxin antidotes (CRCA). R. buchneri expressed both REIS_1423 and REIS_1424 in tick cell culture, and PCR assays showed specificity for R. buchneri over other rickettsiae and utility for positive detection in three tick populations. Finally, phylogenomics analyses uncovered divergent CRCT/CRCA modules in varying states of conservation; however, only R. buchneri and related Tamurae/Ixodes Group rickettsiae carry complete TA modules.

Conclusion

We hypothesize that Rickettsia CRCT/CRCA modules circulate in the Rickettsia mobile gene pool, arming rickettsiae for battle over arthropod colonization. While its functional significance remains to be tested, R. buchneri CRCT/CRCA serves as a marker to positively identify infection and begin deciphering the role this endosymbiont plays in the biology of the blacklegged tick.

Ecological Contacts and Host Specificity Promote Replacement of Nutritional Endosymbionts in Ticks

Symbiosis with vitamin-provisioning microbes is essential for the nutrition of animals with some specialized feeding habits. While coevolution favors the interdependence between symbiotic partners, their associations are not necessarily stable: Recently acquired symbionts can replace ancestral symbionts. In this study, we demonstrate successful replacement by Francisella-like endosymbionts (-LE), a group of B-vitamin-provisioning endosymbionts, across tick communities driven by horizontal transfers. Using a broad collection of Francisella-LE-infected tick species, we determined the diversity of Francisella-LE haplotypes through a multi-locus strain typing approach and further characterized their phylogenetic relationships and their association with biological traits of their tick hosts. The patterns observed showed that Francisella-LE commonly transfer through similar ecological networks and geographic distributions shared among different tick species and, in certain cases, through preferential shuffling across congeneric tick species. Altogether, these findings reveal the importance of geographic, ecological, and phylogenetic proximity in shaping the replacement pattern in which new nutritional symbioses are initiated.

Host-adaptation in Legionellales is 1.9 Ga, coincident with eukaryogenesis

Bacteria adapting to living in a host cell caused the most salient events in the evolution of eukaryotes, namely the seminal fusion with an archaeon, and the emergence of both mitochondrion and chloroplast. A bacterial clade that may hold the key to understanding these events is the deep-branching gammaproteobacterial order Legionellales-containing among others Coxiella and Legionella-of which all known members grow inside eukaryotic cells. Here, by analyzing 35 novel Legionellales genomes mainly acquired through metagenomics, we show that this group is much more diverse than previously thought, and that key host-adaptation events took place very early in its evolution. Crucial virulence factors like the Type IVB secretion (Dot/Icm) system and two shared effector proteins were gained in the last Legionellales common ancestor (LLCA). Many metabolic gene families were lost in LLCA and its immediate descendants, including functions directly and indirectly related to molybdenum metabolism. On the other hand, genome sizes increased in the ancestors of the Legionella genus. We estimate that LLCA lived circa 1.89 Ga ago, probably predating the last eukaryotic common ancestor (LECA) by circa 0.4-1.0 Ga. These elements strongly indicate that host-adaptation arose only once in Legionellales, and that these bacteria were using advanced molecular machinery to exploit and manipulate host cells early in eukaryogenesis.

A dual endosymbiosis supports nutritional adaptation to hematophagy in the invasive tick Hyalomma marginatum

Many animals are dependent on microbial partners that provide essential nutrients lacking from their diet. Ticks, whose diet consists exclusively on vertebrate blood, rely on maternally inherited bacterial symbionts to supply B vitamins. While previously studied tick species consistently harbor a single lineage of those nutritional symbionts, we evidence here that the invasive tick Hyalomma marginatum harbors a unique dual-partner nutritional system between an ancestral symbiont, Francisella, and a more recently acquired symbiont, Midichloria. Using metagenomics, we show that Francisella exhibits extensive genome erosion that endangers the nutritional symbiotic interactions. Its genome includes folate and riboflavin biosynthesis pathways but deprived functional biotin biosynthesis on account of massive pseudogenization. Co-symbiosis compensates this deficiency since the Midichloria genome encompasses an intact biotin operon, which was primarily acquired via lateral gene transfer from unrelated intracellular bacteria commonly infecting arthropods. Thus, in H. marginatum, a mosaic of co-evolved symbionts incorporating gene combinations of distant phylogenetic origins emerged to prevent the collapse of an ancestral nutritional symbiosis. Such dual endosymbiosis was never reported in other blood feeders but was recently documented in agricultural pests feeding on plant sap, suggesting that it may be a key mechanism for advanced adaptation of arthropods to specialized diets.

A Rickettsiella Endosymbiont Is a Potential Source of Essential B-Vitamins for the Poultry Red Mite, Dermanyssus gallinae

Many obligate blood-sucking arthropods rely on symbiotic bacteria to provision essential B vitamins that are either missing or at sub-optimal levels in their nutritionally challenging blood diet. The poultry red mite Dermanyssus gallinae, an obligate blood-feeding ectoparasite, is a serious threat to the hen egg industry. Poultry red mite infestation has a major impact on hen health and welfare and causes a significant reduction in both egg quality and production. Thus far, the identity and biological role of nutrient provisioning bacterial mutualists from D. gallinae are little understood. Here, we demonstrate that an obligate intracellular bacterium of the Rickettsiella genus is detected in D. gallinae mites collected from 63 sites (from 15 countries) across Europe. In addition, we report the genome sequence of Rickettsiella from D. gallinae (Rickettsiella - D. gallinae endosymbiont; Rickettsiella DGE). Rickettsiella DGE has a circular 1.89Mbp genome that encodes 1,973 proteins. Phylogenetic analysis confirms the placement of Rickettsiella DGE within the Rickettsiella genus, related to a facultative endosymbiont from the pea aphid and Coxiella-like endosymbionts (CLEs) from blood feeding ticks. Analysis of the Rickettsiella DGE genome reveals that many protein-coding sequences are either pseudogenized or lost, but Rickettsiella DGE has retained several B vitamin biosynthesis pathways, suggesting the importance of these pathways in evolution of a nutritional symbiosis with D. gallinae. In silico metabolic pathway reconstruction revealed that Rickettsiella DGE is unable to synthesize protein amino acids and, therefore, amino acids are potentially provisioned by the host. In contrast, Rickettsiella DGE retains biosynthetic pathways for B vitamins: thiamine (vitamin B1) via the salvage pathway; riboflavin (vitamin B2) and pyridoxine (vitamin B6) and the cofactors: flavin adenine dinucleotide (FAD) and coenzyme A (CoA) that likely provision these nutrients to the host.

Patterns of Microbiome Variation Among Infrapopulations of Permanent Bloodsucking Parasites

While interspecific variation in microbiome composition can often be readily explained by factors such as host species identity, there is still limited knowledge of how microbiomes vary at scales lower than the species level (e.g., between individuals or populations). Here, we evaluated variation in microbiome composition of individual parasites among infrapopulations (i.e., populations of parasites of the same species living on a single host individual). To address this question, we used genome-resolved and shotgun metagenomic data of 17 infrapopulations (balanced design) of the permanent, bloodsucking seal louse Echinophthirius horridus sampled from individual Saimaa ringed seals Pusa hispida saimensis. Both genome-resolved and read-based metagenomic classification approaches consistently show that parasite infrapopulation identity is a significant factor that explains both qualitative and quantitative patterns of microbiome variation at the intraspecific level. This study contributes to the general understanding of the factors driving patterns of intraspecific variation in microbiome composition, especially of bloodsucking parasites, and has implications for understanding how well-known processes occurring at higher taxonomic levels, such as phylosymbiosis, might arise in these systems.

Antimicrobials Affect the Fat Body Microbiome and Increase the Brown Planthopper Mortality

Symbionts in the abdomen fat body of brown planthopper (BPH) play an important role in the growth and reproduction of their host, Nilaparvata lugens Stål (Hemiptera: Delphacidae). Thus, controlling BPH infection on rice by inhibiting symbionts with antimicrobials is feasible. However, the effect of antimicrobials on the microbiome in the fat body and the relationship between microbial community and mortality have not been fully elucidated. A decrease in the total number of yeast-like symbiotes in the fat body and elevated mortality were observed after exposure to toyocamycin, tebuconazole, and zhongshengmycin. Additionally, we found that the antimicrobials reduced bacterial diversity and increased fungal diversity in the fat body and altered the bacterial and fungal community structure. Although the total absolute abundance of bacteria and fungi decreased after antimicrobial exposure, the absolute abundance of Serratia increased, indicating that Serratia, which was the most dominant in the fat body, is an important symbiont involved in resistance to antimicrobials. After antimicrobial exposure, seven genera, which probably participated in the nutrition and development function of the host, were totally eliminated from the fat body. Overall, our study enriches the knowledge of microbiomes in the fat body of BPH under antimicrobial treatment and the disturbance of symbionts would be further used to help other pesticides to control pests.

A new symbiotic lineage related to Neisseria and Snodgrassella arises from the dynamic and diverse microbiomes in sucking lice

The phylogenetic diversity of symbiotic bacteria in sucking lice suggests that lice have a complex history of symbiont acquisition, loss, and replacement throughout their evolution. These processes have resulted in the establishment of different, phylogenetically distant bacteria as obligate mutualists in different louse groups. By combining metagenomics and amplicon screening across several populations of three louse species (members of the genera Polyplax and Hoplopleura) we describe a novel louse symbiont lineage related to Neisseria and Snodgrassella, and show its independent origin in the two louse genera. While the genomes of these symbionts are highly similar, their respective distributions and status within lice microbiomes indicate that they have different functions and history. In Hoplopleura acanthopus, the Neisseriaceae-related bacterium is a dominant obligate symbiont present across several host populations. In contrast, the Polyplax microbiomes are dominated by the obligate symbiont Legionella polyplacis, with the Neisseriaceae-related bacterium co-occurring only in some samples and with much lower abundance. The results thus support the view that compared to other exclusively blood feeding insects, Anoplura possess a unique capacity to acquire symbionts from diverse groups of bacteria.

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.

"Parasite turnover zone" at secondary contact: A new pattern in host-parasite population genetics

We describe here a new pattern of population genetic structure in a host-parasite system that can arise after secondary contact of previously isolated populations. Due to different generation times, and therefore different tempos of molecular evolution, the host and parasite populations reach different degrees of genetic differentiation during their separation (e.g., in refugia). Consequently, upon secondary contact, the host populations are able to re-establish a single panmictic population across the area of contact, while the parasite populations stop their dispersal at the secondary contact zone and create a narrow hybrid zone. From the host's perspective, the parasite's hybrid zone functions on a microevolutionary scale as a "parasite turnover zone": while the hosts are passing from area A to area B, their parasites turn genetically from the area A genotypes to the area B genotypes. We demonstrate this novel pattern with a model composed of Apodemus mice and Polyplax lice by comparing maternally inherited markers (complete mitochondrial genomes, and complete genomes of the vertically transmitted symbiont Legionella polyplacis) with single nucleotide polymorphisms derived from louse genomic data. We discuss the circumstances that may lead to this pattern and possible reasons why it has been overlooked in studies of host-parasite population genetics.

Convergence of Nutritional Symbioses in Obligate Blood Feeders

Symbiosis with intracellular or gut bacteria is essential for the nutrition of animals with an obligate blood-feeding habit. Divergent bacterial lineages have independently evolved functional interactions with obligate blood feeders, but all converge to an analogous biochemical feature: the provisioning of B vitamins. Although symbionts and blood feeders coevolved interdependently for millions of years we stress that their associations are not necessarily stable. Ancestral symbionts can be replaced by recently acquired bacteria with similar biochemical features, a dynamic that emerges through a combination of phylogenetic and ecological constraints. Specifically, we highlight the lateral transfer of a streamlined biotin (B₇ vitamin) operon, and conjecture that its extensive spread across bacterial lineages may drive the emergence of novel nutritional symbioses with blood feeders.

Evaluation of an ensemble-based distance statistic for clustering MLST datasets using epidemiologically defined clusters of cyclosporiasis

Outbreaks of cyclosporiasis, a food-borne illness caused by the coccidian parasite Cyclospora cayetanensis have increased in the USA in recent years, with approximately 2300 laboratory-confirmed cases reported in 2018. Genotyping tools are needed to inform epidemiological investigations, yet genotyping Cyclospora has proven challenging due to its sexual reproductive cycle which produces complex infections characterized by high genetic heterogeneity. We used targeted amplicon deep sequencing and a recently described ensemble-based distance statistic that accommodates heterogeneous (mixed) genotypes and specimens with partial genotyping data, to genotype and cluster 648 C. cayetanensis samples submitted to CDC in 2018. The performance of the ensemble was assessed by comparing ensemble-identified genetic clusters to analogous clusters identified independently based on common food exposures. Using these epidemiologic clusters as a gold standard, the ensemble facilitated genetic clustering with 93.8% sensitivity and 99.7% specificity. Hence, we anticipate that this procedure will greatly complement epidemiologic investigations of cyclosporiasis.

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.

Insect-Symbiont Gene Expression in the Midgut Bacteriocytes of a Blood-Sucking Parasite

Animals interact with a diverse array of both beneficial and detrimental microorganisms. In insects, these symbioses in many cases allow feeding on nutritionally unbalanced diets. It is, however, still not clear how are obligate symbioses maintained at the cellular level for up to several hundred million years. Exact mechanisms driving host-symbiont interactions are only understood for a handful of model species and data on blood-feeding hosts with intracellular bacteria are particularly scarce. Here, we analyzed interactions between an obligately blood-sucking parasite of sheep, the louse fly Melophagus ovinus, and its obligate endosymbiont, Arsenophonus melophagi. We assembled a reference transcriptome for the insect host and used dual RNA-Seq with five biological replicates to compare expression in the midgut cells specialized for housing symbiotic bacteria (bacteriocytes) to the rest of the gut (foregut-hindgut). We found strong evidence for the importance of zinc in the system likely caused by symbionts using zinc-dependent proteases when acquiring amino acids, and for different immunity mechanisms controlling the symbionts than in closely related tsetse flies. Our results show that cellular and nutritional interactions between this blood-sucking insect and its symbionts are less intimate than what was previously found in most plant-sap sucking insects. This finding is likely interconnected to several features observed in symbionts in blood-sucking arthropods, particularly their midgut intracellular localization, intracytoplasmic presence, less severe genome reduction, and relatively recent associations caused by frequent evolutionary losses and replacements.

Serial horizontal transfer of vitamin-biosynthetic genes enables the establishment of new nutritional symbionts in aphids' di-symbiotic systems

Many insects depend on obligate mutualistic bacteria to provide essential nutrients lacking from their diet. Most aphids, whose diet consists of phloem, rely on the bacterial endosymbiont Buchnera aphidicola to supply essential amino acids and B vitamins. However, in some aphid species, provision of these nutrients is partitioned between Buchnera and a younger bacterial partner, whose identity varies across aphid lineages. Little is known about the origin and the evolutionary stability of these di-symbiotic systems. It is also unclear whether the novel symbionts merely compensate for losses in Buchnera or carry new nutritional functions. Using whole-genome endosymbiont sequences of nine Cinara aphids that harbour an Erwinia-related symbiont to complement Buchnera, we show that the Erwinia association arose from a single event of symbiont lifestyle shift, from a free-living to an obligate intracellular one. This event resulted in drastic genome reduction, long-term genome stasis, and co-divergence with aphids. Fluorescence in situ hybridisation reveals that Erwinia inhabits its own bacteriocytes near Buchnera's. Altogether these results depict a scenario for the establishment of Erwinia as an obligate symbiont that mirrors Buchnera's. Additionally, we found that the Erwinia vitamin-biosynthetic genes not only compensate for Buchnera's deficiencies, but also provide a new nutritional function; whose genes have been horizontally acquired from a Sodalis-related bacterium. A subset of these genes have been subsequently transferred to a new Hamiltonella co-obligate symbiont in one specific Cinara lineage. These results show that the establishment and dynamics of multi-partner endosymbioses can be mediated by lateral gene transfers between co-ocurring symbionts.

The all-intracellular order Legionellales is unexpectedly diverse, globally distributed and lowly abundant

Legionellales is an order of the Gammaproteobacteria, only composed of host-adapted, intracellular bacteria, including the accidental human pathogens Legionella pneumophila and Coxiella burnetii. Although the diversity in terms of lifestyle is large across the order, only a few genera have been sequenced, owing to the difficulty to grow intracellular bacteria in pure culture. In particular, we know little about their global distribution and abundance.

Here, we analyze 16/18S rDNA amplicons both from tens of thousands of published studies and from two separate sampling campaigns in and around ponds and in a silver mine. We demonstrate that the diversity of the order is much larger than previously thought, with over 450 uncultured genera. We show that Legionellales are found in about half of the samples from freshwater, soil and marine environments and quasi-ubiquitous in man-made environments. Their abundance is low, typically 0.1%, with few samples up to 1%. Most Legionellales OTUs are globally distributed, while many do not belong to a previously identified species.

This study sheds a new light on the ubiquity and diversity of one major group of host-adapted bacteria. It also emphasizes the need to use metagenomics to better understand the role of host-adapted bacteria in all environments.

The Importance of Revisiting Legionellales Diversity

Bacteria of the order Legionellales, such as Legionella pneumophila, the agent of Legionnaires' disease, and Coxiella burnetii, the agent of Q fever, are widely recognized as human pathogens. While our view of the Legionellales is often limited to clinical isolates, ecological surveys are continually uncovering new members of the Legionellales that do not fall into the recognized pathogenic species. Here we emphasize that most of these Legionellales are nonpathogenic forms that have evolved symbiotic lifestyles with nonvertebrate hosts. The diversity of nonpathogenic forms remains, however, largely underexplored. We conjecture that its characterization, once contrasted with the data on pathogenic species, will reveal novel highlights on the mechanisms underlying lifestyle transitions of intracellular bacteria, including the emergence of pathogenesis and mutualism, transmission routes, and host specificity.

Host specificity driving genetic structure and diversity in ectoparasite populations: Coevolutionary patterns in Apodemus mice and their lice

A degree of host specificity, manifested by the processes of host-parasite cospeciations and host switches, is assumed to be a major determinant of parasites' evolution. To understand these patterns and formulate appropriate ecological hypotheses, we need better insight into the coevolutionary processes at the intraspecific level, including the maintenance of genetic diversity and population structure of parasites and their hosts. Here, we address these questions by analyzing large-scale molecular data on the louse Polyplax serrata and its hosts, mice of the genus Apodemus, across a broad range of European localities. Using mitochondrial DNA sequences and microsatellite data, we demonstrate the general genetic correspondence of the Apodemus/Polyplax system to the scenario of the postglacial recolonization of Europe, but we also show several striking discrepancies. Among the most interesting are the evolution of different degrees of host specificity in closely related louse lineages in sympatry, or decoupled population structures of the host and parasites in central Europe. We also find strong support for the prediction that parasites with narrower host specificity possess a lower level of genetic diversity and a deeper pattern of interpopulation structure as a result of limited dispersal and smaller effective population size.

Host-symbiont-pathogen interactions in blood-feeding parasites: nutrition, immune cross-talk and gene exchange

Animals are common hosts of mutualistic, commensal and pathogenic microorganisms. Blood-feeding parasites feed on a diet that is nutritionally unbalanced and thus often rely on symbionts to supplement essential nutrients. However, they are also of medical importance as they can be infected by pathogens such as bacteria, protists or viruses that take advantage of the blood-feeding nutritional strategy for own transmission. Since blood-feeding evolved multiple times independently in diverse animals, it showcases a gradient of host-microbe interactions. While some parasitic lineages are possibly asymbiotic and manage to supplement their diet from other food sources, other lineages are either loosely associated with extracellular gut symbionts or harbour intracellular obligate symbionts that are essential for the host development and reproduction. What is perhaps even more diverse are the pathogenic lineages that infect blood-feeding parasites. This microbial diversity not only puts the host into a complicated situation - distinguishing between microorganisms that can greatly decrease or increase its fitness - but also increases opportunity for horizontal gene transfer to occur in this environment. In this review, I first introduce this diversity of mutualistic and pathogenic microorganisms associated with blood-feeding animals and then focus on patterns in their interactions, particularly nutrition, immune cross-talk and gene exchange.

Pangenomic Definition of Prokaryotic Species and the Phylogenetic Structure of Prochlorococcus spp

The pangenome is the collection of all groups of orthologous genes (OGGs) from a set of genomes. We apply the pangenome analysis to propose a definition of prokaryotic species based on identification of lineage-specific gene sets. While being similar to the classical biological definition based on allele flow, it does not rely on DNA similarity levels and does not require analysis of homologous recombination. Hence this definition is relatively objective and independent of arbitrary thresholds. A systematic analysis of 110 accepted species with the largest numbers of sequenced strains yields results largely consistent with the existing nomenclature. However, it has revealed that abundant marine cyanobacteria Prochlorococcus marinus should be divided into two species. As a control we have confirmed the paraphyletic origin of Yersinia pseudotuberculosis (with embedded, monophyletic Y. pestis) and Burkholderia pseudomallei (with B. mallei). We also demonstrate that by our definition and in accordance with recent studies Escherichia coli and Shigella spp. are one species.