Probing mixed-genotype infections II: high multiplicity in natural infections of the trypanosomatid, Crithidia bombi, in its host, Bombus spp

Cultivation of monoxenous trypanosomatids: A minireview

Trypanosomatids are obligatory parasites, some of which are responsible for important human and animal diseases, but the vast majority of trypanosomatids are restricted to invertebrate hosts. Isolation and in vitro cultivation of trypanosomatids from insect hosts enable their description, characterization, and subsequently genetic and genomic studies. However, exact nutritional requirements are still unknown for most trypanosomatids and thus very few defined media are available. This mini review provides information about the role of different ingredients, recommendations and advice on essential supplements and important physicochemical parameters of culture media with the aim of facilitating first attempts to cultivate insect-infesting trypanosomatids, with a focus on monoxenous trypanosomatids.

Genetic variation and microbiota in bumble bees cross-infected by different strains of C. bombi

The bumblebee Bombus terrestris is commonly infected by a trypanosomatid gut parasite Crithidia bombi. This system shows a striking degree of genetic specificity where host genotypes are susceptible to different genotypes of parasite. To a degree, variation in host gene expression underlies these differences, however, the effects of standing genetic variation has not yet been explored. Here we report on an extensive experiment where workers of twenty colonies of B. terrestris were each infected by one of twenty strains of C. bombi. To elucidate the host's genetic bases of susceptibility to infection (measured as infection intensity), we used a low-coverage (~2 x) genome-wide association study (GWAS), based on angsd, and a standard high-coverage (~15x) GWAS (with a reduced set from a 8 x 8 interaction matrix, selected from the full set of twenty). The results from the low-coverage approach remained ambiguous. The high-coverage approach suggested potentially relevant genetic variation in cell surface and adhesion processes. In particular, mucin, a surface mucoglycoprotein, potentially affecting parasite binding to the host gut epithelia, emerged as a candidate. Sequencing the gut microbial community of the same bees showed that the abundance of bacterial taxa, such as Gilliamella, Snodgrassella, or Lactobacillus, differed between 'susceptible' and 'resistant' microbiota, in line with earlier studies. Our study suggests that the constitutive microbiota and binding processes at the cell surface are candidates to affect infection intensity after the first response (captured by gene expression) has run its course. We also note that a low-coverage approach may not be powerful enough to analyse such complex traits. Furthermore, testing large interactions matrices (as with the full 20 x 20 combinations) for the effect of interaction terms on infection intensity seems to blur the specific host x parasite interaction effects, likely because the outcome of an infection is a highly non-linear process dominated by variation in individually different pathways of host defence (immune) responses.

Effects of an alternative host on the prevalence and infection intensity of a bumble bee parasite

Several bee parasites are transmitted through flowers, and some of them can infect multiple host species. Given the shared use of flowers by bee species, parasites can potentially encounter multiple host species, which could affect the evolution of parasite virulence. We used the trypanosomatid parasite Crithidia bombi and its host, the common eastern bumble bee (Bombus impatiens), to explore the effect of infecting an alternative host, the alfalfa leaf-cutter bee (Megachile rotundata), on parasite infectivity and ability to replicate. We conducted a serial passage experiment on primary and alternative hosts, assessing infectivity and intensity of infection during five passes. Parasite cells from each pass through the alternative host were also used to infect a group of primary hosts. We found that serial passes through the alternative host increased infectivity, but there was no effect on intensity of infection. Interestingly, both the probability and intensity of infection on the primary host increased after serial passage through the alternative host. This increase in intensity of infection could be due to maladaptation after selection of new C. bombi strains has occurred in the alternative host. This study suggests that host switching has the potential to affect the adaptation of bee parasites to their hosts.

Why Are Queens Broodless? Failed Nest Initiation Not Linked to Parasites, Mating Status, or Ovary Development in Two Bumble Bee Species of Pyrobombus (Hymenoptera: Apidae: Bombus)

Bumble bees (Bombus [Hymenoptera: Apidae]) are important pollinators for agricultural crops, which has led to their commercial domestication. Despite their importance, little is known about the reproductive biology of bumble bees native to North America. The Hunt bumble bee (Bombus huntii Greene [Hymenoptera: Apidae]) and the Vosnesensky bumble bee (Bombus vosnesenskii Radoszkowski [Hymenoptera: Apidae] are native candidates for commercial production in western North America due to their efficacy in providing commercial pollination services. Availability of pollinators native to the region in which services would be provided would minimize the likelihood of introducing exotic species and spreading novel disease. Some parasites are known to affect bumble bee reproduction, but little is known about their prevalence in North America or how they affect queen success. Only 38% of wild-caught B. huntii and 51% wild-caught B. vosnesenskii queens collected between 2015 and 2017 initiated nests in the laboratory. Our objective was to identify causal factors leading to a queen's inability to oviposit. To address this, we dissected each broodless queen and diagnosed diseases, assessed mating status, and characterized ovary development. Nematodes, arthropods, and microorganisms were detected in both species. Overall, 20% of queens were infected by parasites, with higher rates in B. vosnesenskii. Over 95% of both species were mated, and over 88% had developed ovaries. This suggests that parasitism and mating status were not primary causes of broodlessness. Although some failure to nest can be attributed to assessed factors, additional research is needed to fully understand the challenges presented by captive rearing.

Insect trypanosomatids in Papua New Guinea: high endemism and diversity

The extreme biological diversity of Oceanian archipelagos has long stimulated research in ecology and evolution. However, parasitic protists in this geographic area remained neglected and no molecular analyses have been carried out to understand the evolutionary patterns and relationships with their hosts. Papua New Guinea (PNG) is a biodiversity hotspot containing over 5% of the world's biodiversity in less than 0.5% of the total land area. In the current work, we examined insect heteropteran hosts collected in PNG for the presence of trypanosomatid parasites. The diversity of insect flagellates was analysed, to our knowledge for the first time, east of Wallace's Line, one of the most distinct biogeographic boundaries of the world. Out of 907 investigated specimens from 138 species and 23 families of the true bugs collected in eight localities, 135 (15%) were infected by at least one trypanosomatid species. High species diversity of captured hosts correlated with high diversity of detected trypanosomatids. Of 46 trypanosomatid Typing Units documented in PNG, only eight were known from other geographic locations, while 38 TUs (~83%) have not been previously encountered. The widespread trypanosomatid TUs were found in both widely distributed and endemic/sub-endemic insects. Approximately one-third of the endemic trypanosomatid TUs were found in widely distributed hosts, while the remaining species were confined to endemic and sub-endemic insects. The TUs from PNG form clades with conspicuous host-parasite coevolutionary patterns, as well as those with a remarkable lack of this trait. In addition, our analysis revealed new members of the subfamilies Leishmaniinae and Strigomonadinae, potentially representing new genera of trypanosomatids.

Genomic Variation among Strains of Crithidia bombi and C. expoeki

In this study, we sequenced and analyzed the genomes of 40 strains, in addition to the already-reported two type strains, of two Crithidia species infecting bumblebees in Alaska and Central Europe and demonstrated that different strains of Crithidia bombi and C. expoeki vary considerably in terms of single nucleotide polymorphisms and gene copy number. Based on the genomic structure, phylogenetic analyses, and the pattern of copy number variation, we confirmed the status of C. expoeki as a separate species. The Alaskan populations appear to be clearly separated from those of Central Europe. This pattern fits a scenario of rapid host-parasite coevolution, where the selective advantage of a given parasite strain is only temporary. This study provides helpful insights into possible scenarios of selection and diversification of trypanosomatid parasites.IMPORTANCE A group of trypanosomatid flagellates includes several well-studied medically and economically important parasites of vertebrates and plants. Nevertheless, the vast majority of trypanosomatids infect only insects (mostly flies and true bugs) and, because of that, has attracted little research attention in the past. Of several hundred trypanosomatid species, only four can infect bees (honeybees and bumblebees). Because of such scarcity, these parasites are severely understudied. We analyzed whole-genome information for a total of 42 representatives of bee-infecting trypanosomatids collected in Central Europe and Alaska from a population genetics point of view. Our data shed light on the evolution, selection, and diversification in this important group of trypanosomatid parasites.

Large and finite sample properties of a maximum-likelihood estimator for multiplicity of infection

Reliable measures of transmission intensities can be incorporated into metrics for monitoring disease-control interventions. Genetic (molecular) measures like multiplicity of infection (MOI) have several advantages compared with traditional measures, e.g., R0. Here, we investigate the properties of a maximum-likelihood approach to estimate MOI and pathogen-lineage frequencies. By verifying regulatory conditions, we prove asymptotical unbiasedness, consistency and efficiency of the estimator. Finite sample properties concerning bias and variance are evaluated over a comprehensive parameter range by a systematic simulation study. Moreover, the estimator's sensitivity to model violations is studied. The estimator performs well for realistic sample sizes and parameter ranges. In particular, the lineage-frequency estimates are almost unbiased independently of sample size. The MOI estimate's bias vanishes with increasing sample size, but might be substantial if sample size is too small. The estimator's variance matrix agrees well with the Cramér-Rao lower bound, even for small sample size. The numerical and analytical results of this study can be used for study design. This is exemplified by a malaria data set from Venezuela. It is shown how the results can be used to determine the necessary sample size to achieve certain performance goals. An implementation of the likelihood method and a simulation algorithm for study design, implemented as an R script, is available as S1 File alongside a documentation (S2 File) and example data (S3 File).

Context-dependent medicinal effects of anabasine and infection-dependent toxicity in bumble bees

BACKGROUND

Floral phytochemicals are ubiquitous in nature, and can function both as antimicrobials and as insecticides. Although many phytochemicals act as toxins and deterrents to consumers, the same chemicals may counteract disease and be preferred by infected individuals. The roles of nectar and pollen phytochemicals in pollinator ecology and conservation are complex, with evidence for both toxicity and medicinal effects against parasites. However, it remains unclear how consistent the effects of phytochemicals are across different parasite lineages and environmental conditions, and whether pollinators actively self-medicate with these compounds when infected.

APPROACH

Here, we test effects of the nectar alkaloid anabasine, found in Nicotiana, on infection intensity, dietary preference, and survival and performance of bumble bees (Bombus impatiens). We examined variation in the effects of anabasine on infection with different lineages of the intestinal parasite Crithidia under pollen-fed and pollen-starved conditions.

RESULTS

We found that anabasine did not reduce infection intensity in individual bees infected with any of four Crithidia lineages that were tested in parallel, nor did anabasine reduce infection intensity in microcolonies of queenless workers. In addition, neither anabasine nor its isomer, nicotine, was preferred by infected bees in choice experiments, and infected bees consumed less anabasine than did uninfected bees under no-choice conditions. Furthermore, anabasine exacerbated the negative effects of infection on bee survival and microcolony performance. Anabasine reduced infection in only one experiment, in which bees were deprived of pollen and post-pupal contact with nestmates. In this experiment, anabasine had antiparasitic effects in bees from only two of four colonies, and infected bees exhibited reduced-rather than increased-phytochemical consumption relative to uninfected bees.

CONCLUSIONS

Variation in the effect of anabasine on infection suggests potential modulation of tritrophic interactions by both host genotype and environmental variables. Overall, our results demonstrate that Bombus impatiens prefer diets without nicotine and anabasine, and suggest that the medicinal effects and toxicity of anabasine may be context dependent. Future research should identify the specific environmental and genotypic factors that determine whether nectar phytochemicals have medicinal or deleterious effects on pollinators.

Immune response and gut microbial community structure in bumblebees after microbiota transplants

Microbial communities are a key component of host health. As the microbiota is initially 'foreign' to a host, the host's immune system should respond to its acquisition. Such variation in the response should relate not only to host genetic background, but also to differences in the beneficial properties of the microbiota. However, little is known about such interactions. Here, we investigate the gut microbiota of the bumblebee, Bombus terrestris, which has a protective function against the bee's natural trypanosome gut parasite, Crithidia bombi We transplanted 'resistant' and 'susceptible' microbiota into 'resistant' and 'susceptible' host backgrounds, and studied the activity of the host immune system. We found that bees from different resistance backgrounds receiving a microbiota differed in aspects of their immune response. At the same time, the elicited immune response also depended on the received microbiota's resistance phenotype. Furthermore, the microbial community composition differed between microbiota resistance phenotypes (resistant versus susceptible). Our results underline the complex feedback between the host's ability to potentially exert selection on the establishment of a microbial community and the influence of the microbial community on the host immune response in turn.

Evolution of resistance to single and combined floral phytochemicals by a bumble bee parasite

Repeated exposure to inhibitory compounds can drive the evolution of resistance, which weakens chemical defence against antagonists. Floral phytochemicals in nectar and pollen have antimicrobial properties that can ameliorate infection in pollinators, but evolved resistance among parasites could diminish the medicinal efficacy of phytochemicals. However, multicompound blends, which occur in nectar and pollen, present simultaneous chemical challenges that may slow resistance evolution. We assessed evolution of resistance by the common bumble bee gut parasite Crithidia bombi to two floral phytochemicals, singly and combined, over 6 weeks (~100 generations) of chronic exposure. Resistance of C. bombi increased under single and combined phytochemical exposure, without any associated costs of reduced growth under phytochemical-free conditions. After 6 weeks' exposure, phytochemical concentrations that initially inhibited growth by > 50%, and exceeded concentrations in floral nectar, had minimal effects on evolved parasite lines. Unexpectedly, the phytochemical combination did not impede resistance evolution compared to single compounds. These results demonstrate that repeated phytochemical exposure, which could occur in homogeneous floral landscapes or with therapeutic phytochemical treatment of managed hives, can cause rapid evolution of resistance in pollinator parasites. We discuss possible explanations for submaximal phytochemical resistance in natural populations. Evolved resistance could diminish the antiparasitic value of phytochemical ingestion, weakening an important natural defence against infection.

Insect antimicrobial peptides act synergistically to inhibit a trypanosome parasite

The innate immune system provides protection from infection by producing essential effector molecules, such as antimicrobial peptides (AMPs) that possess broad-spectrum activity. This is also the case for bumblebees, Bombus terrestris, when infected by the trypanosome, Crithidia bombi Furthermore, the expressed mixture of AMPs varies with host genetic background and infecting parasite strain (genotype). Here, we used the fact that clones of C. bombi can be cultivated and kept as strains in medium to test the effect of various combinations of AMPs on the growth rate of the parasite. In particular, we used pairwise combinations and a range of physiological concentrations of three AMPs, namely Abaecin, Defensin and Hymenoptaecin, synthetized from the respective genomic sequences. We found that these AMPs indeed suppress the growth of eight different strains of C. bombi, and that combinations of AMPs were typically more effective than the use of a single AMP alone. Furthermore, the most effective combinations were rarely those consisting of maximum concentrations. In addition, the AMP combination treatments revealed parasite strain specificity, such that strains varied in their sensitivity towards the same mixtures. Hence, variable expression of AMPs could be an alternative strategy to combat highly variable infections.This article is part of the themed issue 'Evolutionary ecology of arthropod antimicrobial peptides'.

Genome of Leptomonas pyrrhocoris: a high-quality reference for monoxenous trypanosomatids and new insights into evolution of Leishmania

Many high-quality genomes are available for dixenous (two hosts) trypanosomatid species of the genera Trypanosoma, Leishmania, and Phytomonas, but only fragmentary information is available for monoxenous (single-host) trypanosomatids. In trypanosomatids, monoxeny is ancestral to dixeny, thus it is anticipated that the genome sequences of the key monoxenous parasites will be instrumental for both understanding the origin of parasitism and the evolution of dixeny. Here, we present a high-quality genome for Leptomonas pyrrhocoris, which is closely related to the dixenous genus Leishmania. The L. pyrrhocoris genome (30.4 Mbp in 60 scaffolds) encodes 10,148 genes. Using the L. pyrrhocoris genome, we pinpointed genes gained in Leishmania. Among those genes, 20 genes with unknown function had expression patterns in the Leishmania mexicana life cycle suggesting their involvement in virulence. By combining differential expression data for L. mexicana, L. major and Leptomonas seymouri, we have identified several additional proteins potentially involved in virulence, including SpoU methylase and U3 small nucleolar ribonucleoprotein IMP3. The population genetics of L. pyrrhocoris was also addressed by sequencing thirteen strains of different geographic origin, allowing the identification of 1,318 genes under positive selection. This set of genes was significantly enriched in components of the cytoskeleton and the flagellum.

Possible Synergistic Effects of Thymol and Nicotine Against Crithidia bombi Parasitism in Bumble Bees

Floral nectar contains secondary compounds with antimicrobial properties that can affect not only plant-pollinator interactions, but also interactions between pollinators and their parasites. Although recent work has shown that consumption of plant secondary compounds can reduce pollinator parasite loads, little is known about the effects of dosage or compound combinations. We used the generalist pollinator Bombus impatiens and its obligate gut parasite Crithidia bombi to study the effects of nectar chemistry on host-parasite interactions. In two experiments we tested (1) whether the secondary compounds thymol and nicotine act synergistically to reduce parasitism, and (2) whether dietary thymol concentration affects parasite resistance. In both experiments, uninfected Bombus impatiens were inoculated with Crithidia and then fed particular diet treatments for 7 days, after which infection levels were assessed. In the synergism experiment, thymol and nicotine alone and in combination did not significantly affect parasite load or host mortality. However, the thymol-nicotine combination treatment reduced log-transformed parasite counts by 30% relative to the control group (P = 0.08). For the experiment in which we manipulated thymol concentration, we found no significant effect of any thymol concentration on Crithidia load, but moderate (2 ppm) thymol concentrations incurred a near-significant increase in mortality (P = 0.054). Our results tentatively suggest the value of a mixed diet for host immunity, yet contrast with research on the antimicrobial activity of dietary thymol and nicotine in vertebrate and other invertebrate systems. We suggest that future research evaluate genetic variation in Crithidia virulence, multi-strain competition, and Crithidia interactions with the gut microbe community that may mediate antimicrobial activities of secondary compounds.

Variable effects of nicotine, anabasine, and their interactions on parasitized bumble bees

Secondary metabolites in floral nectar have been shown to reduce parasite load in two common bumble bee species. Previous studies on the effects of nectar secondary metabolites on parasitized bees have focused on single compounds in isolation; however, in nature, bees are simultaneously exposed to multiple compounds. We tested for interactions between the effects of two alkaloids found in the nectar of Nicotiana spp. plants, nicotine and anabasine, on parasite load and mortality in bumble bees ( Bombus impatiens) infected with the intestinal parasite Crithidia bombi. Adult worker bees inoculated with C. bombi were fed nicotine and anabasine diet treatments in a factorial design, resulting in four nectar treatment combinations:  2 ppm nicotine, 5 ppm anabasine, 2ppm nicotine and 5 ppm anabasine together, or a control alkaloid-free solution. We conducted the experiment twice: first, with bees incubated under variable environmental conditions ('Variable'; temperatures varied from 10-35°C with ambient lighting); and second, under carefully controlled environmental conditions ('Stable'; 27°C incubator, constant darkness). In 'Variable', each alkaloid alone significantly decreased parasite loads, but this effect was not realized with the alkaloids in combination, suggesting an antagonistic interaction. Nicotine but not anabasine significantly increased mortality, and the two compounds had no interactive effects on mortality. In 'Stable', nicotine significantly increased parasite loads, the opposite of its effect in 'Variable'. While not significant, the relationship between anabasine and parasite loads was also positive. Interactive effects between the two alkaloids on parasite load were non-significant, but the pattern of antagonistic interaction was similar to that in the variable experiment. Neither alkaloid, nor their interaction, significantly affected mortality under controlled conditions. Our results do not indicate synergy between Nicotiana nectar alkaloids; however, they do suggest a complex interaction between secondary metabolites, parasites, and environmental variables, in which secondary metabolites can be either toxic or medicinal depending on context.

The distribution of parasite strains among hosts affects disease spread in a social insect

Social insects present highly interesting and experimentally amenable systems for the study of disease transmission because they naturally live in dense groups of frequently interacting individuals. Using experimental inoculations of five trypanosomatid strains into groups of its natural host, the bumblebee Bombus terrestris, we investigate the effects of the initial parasite strain distribution across group members on the establishment and transmission success of the different strains to new hosts. For a given number of parasite strains circulating within a host group, transmission to new hosts was increased when the strains were initially inoculated as mixed infections (as opposed to separate single infections), presumably because mixed infections generally favored fast replicating strains. In contrast, separate single infections reduced transmission at least in part through a precedence effect, whereby weak strains appeared to persist by making their host unavailable to superinfection. These results suggest that host groups could benefit from 'compartmentalizing' infections by different parasite strains across different group members, which might be achieved in social insects, for example, by division of labor.

A likelihood approach to estimate the number of co-infections

The number of co-infections of a pathogen (multiplicity of infection or MOI) is a relevant parameter in epidemiology as it relates to transmission intensity. Notably, such quantities can be built into a metric in the context of disease control and prevention. Having applications to malaria in mind, we develop here a maximum-likelihood (ML) framework to estimate the quantities of interest at low computational and no additional costs to study designs or data collection. We show how the ML estimate for the quantities of interest and corresponding confidence-regions are obtained from multiple genetic loci. Assuming specifically that infections are rare and independent events, the number of infections per host follows a conditional Poisson distribution. Under this assumption, we show that a unique ML estimate for the parameter () describing MOI exists which is found by a simple recursion. Moreover, we provide explicit formulas for asymptotic confidence intervals, and show that profile-likelihood-based confidence intervals exist, which are found by a simple two-dimensional recursion. Based on the confidence intervals we provide alternative statistical tests for the MOI parameter. Finally, we illustrate the methods on three malaria data sets. The statistical framework however is not limited to malaria.

Can host ecology and kin selection predict parasite virulence?

Parasite virulence, or the damage a parasite does to its host, is measured in terms of both host costs (reductions in host growth, reproduction and survival) and parasite benefits (increased transmission and parasite numbers) in the literature. Much work has shown that ecological and genetic factors can be strong selective forces in virulence evolution. This review uses kin selection theory to explore how variations in host ecological parameters impact the genetic relatedness of parasite populations and thus virulence. We provide a broad overview of virulence and population genetics studies and then draw connections to existing knowledge about natural parasite populations. The impact of host movement (transporting parasites) and host resistance (filtering parasites) on the genetic structure and virulence of parasite populations is explored, and empirical studies of these factors using Plasmodium and trematode systems are proposed.

Gene expression differences underlying genotype-by-genotype specificity in a host-parasite system

In many systems, host-parasite evolutionary dynamics have led to the emergence and maintenance of diverse parasite and host genotypes within the same population. Genotypes vary in key attributes: Parasite genotypes vary in ability to infect, host genotypes vary in susceptibility, and infection outcome is frequently the result of both parties' genotypic identities. These host-parasite genotype-by-genotype (GH × GP) interactions influence evolutionary and ecological dynamics in important ways. Interactions can be produced through genetic variation; however, here, we assess the role of variable gene expression as an additional source of GH × GP interactions. The bumblebee Bombus terrestris and its trypanosome gut parasite Crithidia bombi are a model system for host-parasite matching. Full-transcriptome sequencing of the bumblebee host revealed that different parasite genotypes indeed induce fundamentally different host expression responses and host genotypes vary in their responses to the infecting parasite genotype. It appears that broadly and successfully infecting parasite genotypes lead to reduced host immune gene expression relative to unexposed bees but induce the expression of genes responsible for controlling gene expression. Contrastingly, a poorly infecting parasite genotype induced the expression of immunologically important genes, including antimicrobial peptides. A targeted expression assay confirmed the transcriptome results and also revealed strong host genotype effects. In all, the expression of a number of genes depends on the host genotype and the parasite genotype and the interaction between both host and parasite genotypes. These results suggest that alongside sequence variation in coding immunological genes, variation that controls immune gene expression can also produce patterns of host-parasite specificity.

Probing mixed-genotype infections I: extraction and cloning of infections from hosts of the trypanosomatid Crithidia bombi

We here present an efficient, precise and reliable method to isolate and cultivate healthy and viable single Crithidia bombi cells from bumblebee faeces using flow cytometry. We report a precision of >99% in obtaining single trypanosomatid cells for further culture and analysis (“cloning”). In the study, we have investigated the use of different liquid media to cultivate C. bombi and present an optimal medium for obtaining viable clones from all tested, infected host donors. We show that this method can be applied to genotype a collection of clones from natural infections. Furthermore, we show how to cryo-preserve C. bombi cells to be revived to become infective clones after at least 4 years of storage. Considering the high prevalence of infections in natural populations, our method provides a powerful tool in studying the level and diversity of these infections, and thus enriches the current methodology for the studies of complex host-parasite interactions.