Contribution of parasite and host genotype to immunopathology of schistosome infections.

BACKGROUND: The role of pathogen genotype in determining disease severity and immunopathology has been studied intensively in microbial pathogens including bacteria, fungi, protozoa and viruses but is poorly understood in parasitic helminths. The medically important blood fluke Schistosoma mansoni is an excellent model system to study the impact of helminth genetic variation on immunopathology. Our laboratory has demonstrated that laboratory schistosome populations differ in sporocyst growth and cercarial production in the intermediate snail host and worm establishment and fecundity in the vertebrate host. Here, we (i) investigate the hypothesis that schistosome genotype plays a significant role in immunopathology and related parasite life history traits in the vertebrate mouse host and (ii) quantify the relative impact of parasite and host genetics on infection outcomes. METHODS: We infected BALB/c and C57BL/6 mice with four different laboratory schistosome populations from Africa and the Americas. We quantified disease progression in the vertebrate host by measuring body weight and complete blood count (CBC) with differential over a 12-week infection period. On sacrifice, we assessed parasitological (egg and worm counts, fecundity), immunopathological (organ measurements and histopathology) and immunological (CBC with differential and cytokine profiles) characteristics to determine the impact of parasite and host genetics. RESULTS: We found significant variation between parasite populations in worm numbers, fecundity, liver and intestine egg counts, liver and spleen weight, and fibrotic area but not in granuloma size. Variation in organ weight was explained by egg burden and intrinsic parasite factors independent of egg burden. We found significant variation between infected mouse lines in cytokine levels (IFN-γ, TNF-α), eosinophils, lymphocytes and monocyte counts. CONCLUSIONS: This study showed that both parasite and host genotype impact the outcome of infection. While host genotype explains most of the variation in immunological traits, parasite genotype explains most of the variation in parasitological traits, and both host and parasite genotypes impact immunopathology outcomes.

Decoding bee cleptoparasitism through comparative transcriptomics of Coelioxoides waltheriae and its host Tetrapedia diversipes.

Cleptoparasitism, also known as brood parasitism, is a widespread strategy among bee species in which the parasite lays eggs into the nests of the host species. Even though this behavior has significant ecological implications for the dynamics of several species, little is known about the molecular pathways associated with cleptoparasitism. To shed some light on this issue, we used gene expression data to perform a comparative analysis between two solitary neotropical bees: Coelioxoides waltheriae, an obligate parasite, and their specific host Tetrapedia diversipes. We found that ortholog genes involved in signal transduction, sensory perception, learning, and memory formation were differentially expressed between the cleptoparasite and the host. We hypothesize that these genes and their associated molecular pathways are engaged in cleptoparasitism-related processes and, hence, are appealing subjects for further investigation into functional and evolutionary aspects of cleptoparasitism in bees.

Characterizing genetic variation on the Z chromosome in Schistosoma japonicum reveals host-parasite co-evolution.

BACKGROUND: Schistosomiasis is a neglected tropical disease that afflicts millions of people worldwide; it is caused by Schistosoma, the only dioecious flukes with ZW systems. Schistosoma japonicum is endemic to Asia; the Z chromosome of S. japonicum comprises one-quarter of the entire genome. Detection of positive selection using resequencing data to understand adaptive evolution has been applied to a variety of pathogens, including S. japonicum. However, the contribution of the Z chromosome to evolution and adaptation is often neglected. METHODS: We obtained 1,077,526 high-quality SNPs on the Z chromosome in 72 S. japonicum using re-sequencing data publicly. To examine the faster Z effect, we compared the sequence divergence of S. japonicum with two closely related species, Schistosoma haematobium and S. mansoni. Genetic diversity was compared between the Z chromosome and autosomes in S. japonicum by calculating the nucleotide diversity (π) and Dxy values. Population structure was also assessed based on PCA and structure analysis. Besides, we employed multiple methods including Tajima's D, FST, iHS, XP-EHH, and CMS to detect positive selection signals on the Z chromosome. Further RNAi knockdown experiments were performed to investigate the potential biological functions of the candidate genes. RESULTS: Our study found that the Z chromosome of S. japonicum showed faster evolution and more pronounced genetic divergence than autosomes, although the effect may be smaller than the variation among genes. Compared with autosomes, the Z chromosome in S. japonicum had a more pronounced genetic divergence of sub-populations. Notably, we identified a set of candidate genes associated with host-parasite co-evolution. In particular, LCAT exhibited significant selection signals within the Taiwan population. Further RNA interference experiments suggested that LCAT is necessary for S. japonicum survival and propagation in the definitive host. In addition, we identified several genes related to the specificity of the intermediate host in the C-M population, including Rab6 and VCP, which are involved in adaptive immune evasion to the host. CONCLUSIONS: Our study provides valuable insights into the adaptive evolution of the Z chromosome in S. japonicum and further advances our understanding of the co-evolution of this medically important parasite and its hosts.

Physiological Measurements and Transcriptomics Reveal the Fitness Costs of Monochamus saltuarius to Bursaphelenchus xylophilus.

The pine wood nematode (PWN) uses several Monochamus species as vehicles, through a temporary hitchhiking process known as phoresy, enabling it to access new host plant resources. Monochamus saltuarius acts as a new and major vector of the PWN in Northeastern China, showing lower PWN carrying capacity and a shorter transmission cycle compared to established vectors. The apparently altered symbiotic relationship offers an interesting area for researching the costs and adaptions involved in nematode-beetle, a specialized phoresy. We analyzed the response and fitness costs of M. saltuarius through physiological measurements and transcriptomics. The PWN exerted adverse repercussions on the growth and development of M. saltuarius. The PWN accelerated larval development into pupae, while beetle adults carrying the PWN exhibited an elevated abnormality rate and mortality, and reduced starvation resistance. During the pupal stage, the expression of growth-related genes, including ecdysone-inducible genes (E74EA), cuticle proteins, and chitin genes (CHTs), markedly increased. Meanwhile, the induced immune response, mainly by the IMD and Toll signaling pathways, could be a contributing factor to adult abnormality and mortality. Adult gonads and trachea exhibited enrichment in pathways related to fatty acid elongation, biosynthesis, and metabolism. FASN, ELOVL, and SCD possibly contributed to resistance against PWN. Our research indicated that phoretic interactions between vector beetles and PWN vary throughout the vector's lifespan, particularly before and after entry into the trachea. This study highlighted the fitness costs of immunity and metabolism on the vector beetle, indicating the adaptation mechanisms and evolutionary trade-offs to PWN.

Polymorphism in the aggressive mimicry lure of the parasitic freshwater mussel Lampsilis fasciola.

Unionoid freshwater mussels (Bivalvia: Unionidae) are free-living apart from a brief, obligately parasitic, larval stage that infects fish hosts, and gravid female mussels have evolved a spectrum of strategies to infect fish hosts with their larvae. In many North American species, this involves displaying a mantle lure: a pigmented fleshy extension that acts as an aggressive mimic of a host fish prey, thereby eliciting a feeding response that results in host infection. The mantle lure of Lampsilis fasciola is of particular interest because it is apparently polymorphic, with two distinct primary lure phenotypes. One, described as "darter-like", has "eyespots", a mottled body coloration, prominent marginal extensions, and a distinct "tail". The other, described as "worm-like", lacks those features and has an orange and black coloration. We investigated this phenomenon using genomics, captive rearing, biogeographic, and behavioral analyses. Within-brood lure variation and within-population phylogenomic (ddRAD-seq) analyses of individuals bearing different lures confirmed that this phenomenon is a true polymorphism. The relative abundance of the two morphs appears stable over ecological timeframes: the ratio of the two lure phenotypes in a River Raisin (MI) population in 2017 was consistent with that of museum samples collected at the same site six decades earlier. Within the River Raisin, four main "darter-like" lure motifs visually approximated four co-occurring darter species (Etheostoma blennioides, E. exile, E. microperca, and Percina maculata), and the "worm-like" lure resembled a widespread common leech, Macrobdella decora. Darters and leeches are typical prey of Micropterus dolomieui (smallmouth bass), the primary fish host of L. fasciola. In situ field recordings of the L. fasciola "darter" and "leech" lure display behaviors, and the lure display of co-occurring congener L. cardium, were captured. Despite having putative models in distinct phyla, both L. fasciola lure morphs have largely similar display behaviors that differ significantly from that of sympatric L. cardium individuals. Some minor differences in the behavior between the two L. fasciola morphs were observed, but we found no clear evidence for a behavioral component of the polymorphism given the criteria measured. Discovery of discrete within-brood inheritance of the lure polymorphism implies potential control by a single genetic locus and identifies L. fasciola as a promising study system to identify regulatory genes controlling a key adaptive trait of freshwater mussels.

New insights into Trypanosoma cruzi genetic diversity, and its influence on parasite biology and clinical outcomes.

Chagas disease, caused by Trypanosoma cruzi, remains a serious public health problem worldwide. The parasite was subdivided into six distinct genetic groups, called "discrete typing units" (DTUs), from TcI to TcVI. Several studies have indicated that the heterogeneity of T. cruzi species directly affects the diversity of clinical manifestations of Chagas disease, control, diagnosis performance, and susceptibility to treatment. Thus, this review aims to describe how T. cruzi genetic diversity influences the biology of the parasite and/or clinical parameters in humans. Regarding the geographic dispersion of T. cruzi, evident differences were observed in the distribution of DTUs in distinct areas. For example, TcII is the main DTU detected in Brazilian patients from the central and southeastern regions, where there are also registers of TcVI as a secondary T. cruzi DTU. An important aspect observed in previous studies is that the genetic variability of T. cruzi can impact parasite infectivity, reproduction, and differentiation in the vectors. It has been proposed that T. cruzi DTU influences the host immune response and affects disease progression. Genetic aspects of the parasite play an important role in determining which host tissues will be infected, thus heavily influencing Chagas disease's pathogenesis. Several teams have investigated the correlation between T. cruzi DTU and the reactivation of Chagas disease. In agreement with these data, it is reasonable to suppose that the immunological condition of the patient, whether or not associated with the reactivation of the T. cruzi infection and the parasite strain, may have an important role in the pathogenesis of Chagas disease. In this context, understanding the genetics of T. cruzi and its biological and clinical implications will provide new knowledge that may contribute to additional strategies in the diagnosis and clinical outcome follow-up of patients with Chagas disease, in addition to the reactivation of immunocompromised patients infected with T. cruzi.

Combining a transcriptomic approach and a targeted metabolomics approach for deciphering the molecular bases of compatibility phenotype in the snail Biomphalaria glabrata toward Schistosoma mansoni.

Biomphalaria glabrata is a freshwater snail and the obligatory intermediate host of Schistosoma mansoni parasite, the etiologic agent of intestinal Schistosomiasis, in South America and Caribbean. Interestingly in such host-parasite interactions, compatibility varies between populations, strains or individuals. This observed compatibility polymorphism is based on a complex molecular-matching-phenotype, the molecular bases of which have been investigated in numerous studies, notably by comparing between different strains or geographical isolates or clonal selected snail lines. Herein we propose to decipher the constitutive molecular support of this interaction in selected non-clonal resistant and susceptible snail strain originating from the same natural population from Brazil and thus having the same genetic background. Thanks to a global RNAseq transcriptomic approach on whole snail, we identified a total of 328 differentially expressed genes between resistant and susceptible phenotypes among which 129 were up-regulated and 199 down-regulated. Metabolomic studies were used to corroborate the RNAseq results. The activation of immune genes and specific metabolic pathways in resistant snails might provide them with the capacity to better respond to parasite infection.

Dispersal stabilizes coupled ecological and evolutionary dynamics in a host-parasitoid system.

When ecological and evolutionary dynamics occur on comparable timescales, persistence of the ensuing eco-evolutionary dynamics requires both ecological and evolutionary stability. This unites key questions in ecology and evolution: How do species coexist, and what maintains genetic variation in a population? In this work, we investigated a host-parasitoid system in which pea aphid hosts rapidly evolve resistance to Aphidius ervi parasitoids. Field data and mathematical simulations showed that heterogeneity in parasitoid dispersal can generate variation in parasitism-mediated selection on hosts through time and space. Experiments showed how evolutionary trade-offs plus moderate host dispersal across this selection mosaic cause host-parasitoid coexistence and maintenance of genetic variation in host resistance. Our results show how dispersal can stabilize both the ecological and evolutionary components of eco-evolutionary dynamics.

Manipulating multi-level selection in a fungal entomopathogen reveals social conflicts and a method for improving biocontrol traits.

Changes in parasite virulence are commonly expected to lead to trade-offs in other life history traits that can affect fitness. Understanding these trade-offs is particularly important if we want to manipulate the virulence of microbial biological control agents. Theoretically, selection across different spatial scales, i.e. between- and within-hosts, shapes these trade-offs. However, trade-offs are also dependent on parasite biology. Despite their applied importance the evolution of virulence in fungal parasites is poorly understood: virulence can be unstable in culture and commonly fails to increase in simple passage experiments. We hypothesized that manipulating selection intensity at different scales would reveal virulence trade-offs in a fungal pathogen of aphids, Akanthomyces muscarius. Starting with a genetically diverse stock we selected for speed of kill, parasite yield or infectivity by manipulating competition within and between hosts and between-populations of hosts over 7 rounds of infection. We characterized ancestral and evolved lineages by whole genome sequencing and by measuring virulence, growth rate, sporulation and fitness. While several lineages showed increases in virulence, we saw none of the trade-offs commonly found in obligately-killing parasites. Phenotypically similar lineages within treatments often shared multiple single-nucleotide variants, indicating strong convergent evolution. The most dramatic phenotypic changes were in timing of sporulation and spore production in vitro. We found that early sporulation led to reduced competitive fitness but could increase yield of spores on media, a trade-off characteristic of social conflict. Notably, the selection regime with strongest between-population competition and lowest genetic diversity produced the most consistent shift to early sporulation, as predicted by social evolution theory. Multi-level selection therefore revealed social interactions novel to fungi and showed that these biocontrol agents have the genomic flexibility to improve multiple traits-virulence and spore production-that are often in conflict in other parasites.

Genetic factors underlying host resistance to Rhipicephalus microplus tick infestation in Braford cattle: a systems biology perspective.

Approximately 80% of the world's cattle are raised in regions with a high risk of tick-borne diseases, resulting in significant economic losses due to parasitism by Rhipicephalus (Boophilus) microplus. However, the lack of a systemic biology approach hampers a comprehensive understanding of tick-host interactions that mediate tick resistance phenotypes. Here, we conducted a genome-wide association study (GWAS) of 2933 Braford cattle and found 340 single-nucleotide polymorphisms (SNPs) associated with tick counts. Gene expression analyses were performed on skin samples obtained from previously tick-exposed heifers with extremely high or low estimated breeding values for R. microplus counts. Evaluations were performed both before and after artificial infestation with ticks. Differentially expressed genes were found within 1-Mb windows centered at significant SNPs from GWAS. A total of 330 genes were related to the breakdown of homeostasis that was induced by larval attachment to bovine skin. Enrichment analysis pointed to a key role of proteolysis and signal transduction via JAK/STAT, NFKB and WNT/beta catenin signaling pathways. Integrative analysis on matrixEQTL revealed two cis-eQTLs and four significant SNPs in the genes peptidyl arginine deiminase type IV (PADI4) and LOC11449251. The integration of genomic data from QTL maps and transcriptome analyses has identified a set of twelve key genes that show significant associations with tick loads. These genes could be key candidates to improve the accuracy of genomic predictions for tick resistance in Braford cattle.

High specificity of symbiont-conferred resistance in an aphid-parasitoid field community.

Host-parasite coevolution is mediated by genetic interactions between the antagonists and may lead to reciprocal adaptation. In the black bean aphid, Aphis fabae fabae, resistance to parasitoids can be conferred by the heritable bacterial endosymbiont Hamiltonella defensa. H. defensa has been shown to be variably protective against different parasitoid species, and different genotypes of the black bean aphid's main parasitoid Lysiphlebus fabarum. However, these results were obtained using haphazard combinations of laboratory-reared insect lines with different origins, making it unclear how representative they are of natural, locally (co)adapted communities. We therefore comprehensively sampled the parasitoids of a natural A. f. fabae population and measured the ability of the five most abundant species to parasitize aphids carrying the locally prevalent H. defensa haplotypes. H. defensa provided resistance only against the dominant parasitoid L. fabarum (70% of all parasitoids), but not against less abundant parasitoids, and resistance to L. fabarum acted in a genotype-specific manner (G × G interactions between H. defensa and L. fabarum). These results confirm that strong species- and genotype-specificity of symbiont-conferred resistance is indeed a hallmark of wild A. f. fabae populations, and they are consistent with symbiont-mediated adaptation of aphids to the parasitoids posing the highest risk.

The evolution of constitutively active humoral immune defenses in Drosophila populations under high parasite pressure.

Both constitutive and inducible immune mechanisms are employed by hosts for defense against infection. Constitutive immunity allows for a faster response, but it comes with an associated cost that is always present. This trade-off between speed and fitness costs leads to the theoretical prediction that constitutive immunity will be favored where parasite exposure is frequent. We selected populations of Drosophila melanogaster under high parasite pressure from the parasitoid wasp Leptopilina boulardi. With RNA sequencing, we found the evolution of resistance in these populations was associated with them developing constitutively active humoral immunity, mediated by the larval fat body. Furthermore, these evolved populations were also able to induce gene expression in response to infection to a greater level, which indicates an overall more activated humoral immune response to parasitization. The anti-parasitoid immune response also relies on the JAK/STAT signaling pathway being activated in muscles following infection, and this induced response was only seen in populations that had evolved under high parasite pressure. We found that the cytokine Upd3, which induces this JAK/STAT response, is being expressed by immature lamellocytes. Furthermore, these immune cells became constitutively present when populations evolved resistance, potentially explaining why they gained the ability to activate JAK/STAT signaling. Thus, under intense parasitism, populations evolved resistance by increasing both constitutive and induced immune defenses, and there is likely an interplay between these two forms of immunity.

Evaluating evidence for co-geography in the Anopheles-Plasmodium host-parasite system.

The often tight association between parasites and their hosts means that under certain scenarios, the evolutionary histories of the two species can become closely coupled both through time and across space. Using spatial genetic inference, we identify a potential signal of common dispersal patterns in the Anopheles gambiae and Plasmodium falciparum host-parasite system as seen through a between-species correlation of the differences between geographic sampling location and geographic location predicted from the genome. This correlation may be due to coupled dispersal dynamics between host and parasite but may also reflect statistical artifacts due to uneven spatial distribution of sampling locations. Using continuous-space population genetics simulations, we investigate the degree to which uneven distribution of sampling locations leads to bias in prediction of spatial location from genetic data and implement methods to counter this effect. We demonstrate that while algorithmic bias presents a problem in inference from spatio-genetic data, the correlation structure between A. gambiae and P. falciparum predictions cannot be attributed to spatial bias alone and is thus likely a genetic signal of co-dispersal in a host-parasite system.

Separation of evolutionary timescales in coevolving species.

Many coevolutionary processes, including host-parasite and host-symbiont interactions, involve one species or trait which evolves much faster than the other. Whether or not a coevolutionary trajectory converges depends on the relative rates of evolutionary change in the two species, and so current adaptive dynamics approaches generally either determine convergence stability by considering arbitrary (often comparable) rates of evolutionary change or else rely on necessary or sufficient conditions for convergence stability. We propose a method for determining convergence stability in the case where one species is expected to evolve much faster than the other. This requires a second separation of timescales, which assumes that the faster evolving species will reach its evolutionary equilibrium (if one exists) before a new mutation arises in the more slowly evolving species. This method, which is likely to be a reasonable approximation for many coevolving species, both provides straightforward conditions for convergence stability and is less computationally expensive than traditional analysis of coevolution models, as it reduces the trait space from a two-dimensional plane to a one-dimensional manifold. In this paper, we present the theory underlying this new separation of timescales and provide examples of how it could be used to determine coevolutionary outcomes from models.

Rapid bacteria-phage coevolution drives the emergence of multiscale networks.

Interactions between species catalyze the evolution of multiscale ecological networks, including both nested and modular elements that regulate the function of diverse communities. One common assumption is that such complex pattern formation requires spatial isolation or long evolutionary timescales. We show that multiscale network structure can evolve rapidly under simple ecological conditions without spatial structure. In just 21 days of laboratory coevolution, Escherichia coli and bacteriophage Φ21 coevolve and diversify to form elaborate cross-infection networks. By measuring ~10,000 phage-bacteria infections and testing the genetic basis of interactions, we identify the mechanisms that create each component of the multiscale pattern. Our results demonstrate how multiscale networks evolve in parasite-host systems, illustrating Darwin's idea that simple adaptive processes can generate entangled banks of ecological interactions.

Lack of host specialization despite selective host use in brood parasitic cuckoo catfish.

Host-parasite dynamics involve coevolutionary arms races, which may lead to host specialization and ensuing diversification. Our general understanding of the evolution of host specialization in brood parasites is compromised by a restricted focus on bird and insect lineages. The cuckoo catfish (Synodontis multipunctatus) is an obligate parasite of parental care of mouthbrooding cichlids in Lake Tanganyika. Given the ecological and taxonomic diversity of mouthbrooding cichlids in the lake, we hypothesized the existence of sympatric host-specific lineages in the cuckoo catfish. In a sample of 779 broods from 20 cichlid species, we found four species parasitized by cuckoo catfish (with prevalence of parasitism of 2%-18%). All parasitized cichlids were from the tribe Tropheini, maternal mouthbrooders that spawn over a substrate (rather than in open water). Phylogenetic analysis based on genomic (ddRAD sequencing) and mitochondrial (Dloop) data from cuckoo catfish embryos showed an absence of host-specific lineages. This was corroborated by analyses of genetic structure and co-ancestry matrix. Within host species, parasitism was not associated with any individual characteristic we recorded (parent size, water depth), but was costly as parasitized parents carried smaller clutches of their own offspring. We conclude that the cuckoo catfish is an intermediate generalist and discuss costs, benefits and constraints of host specialization in this species and brood parasites in general.

Leishmania genetic exchange is mediated by IgM natural antibodies.

Host factors that mediate Leishmania genetic exchange are not well defined. Here we demonstrate that natural IgM (IgMn)1-4 antibodies mediate parasite genetic exchange by inducing the transient formation of a spherical parasite clump that promotes parasite fusion and hybrid formation. We establish that IgMn from Leishmania-free animals binds to the surface of Leishmania parasites to induce significant changes in the expression of parasite transcripts and proteins. Leishmania binding to IgMn is partially lost after glycosidase treatment, although parasite surface phosphoglycans, including lipophosphoglycan, are not required for IgMn-induced parasite clumping. Notably, the transient formation of parasite clumps is essential for Leishmania hybridization in vitro. In vivo, we observed a 12-fold increase in hybrid formation in sand flies provided a second blood meal containing IgMn compared with controls. Furthermore, the generation of recombinant progeny from mating hybrids and parental lines were only observed in sand flies provided with IgMn. Both in vitro and in vivo IgM-induced Leishmania crosses resulted in full genome hybrids that show equal patterns of biparental contribution. Leishmania co-option of a host natural antibody to facilitate mating in the insect vector establishes a new paradigm of parasite-host-vector interdependence that contributes to parasite diversity and fitness by promoting genetic exchange.

Massive horizontal gene transfer and the evolution of nematomorph-driven behavioral manipulation of mantids.

To complete their life cycle, a wide range of parasites must manipulate the behavior of their hosts.1 This manipulation is a well-known example of the "extended phenotype,2" where genes in one organism have phenotypic effects on another organism. Recent studies have explored the parasite genes responsible for such manipulation of host behavior, including the potential molecular mechanisms.3,4 However, little is known about how parasites have acquired the genes involved in manipulating phylogenetically distinct hosts.4 In a fascinating example of the extended phenotype, nematomorph parasites have evolved the ability to induce their terrestrial insect hosts to enter bodies of water, where the parasite then reproduces. Here, we comprehensively analyzed nematomorphs and their mantid hosts, focusing on the transcriptomic changes associated with host manipulations and sequence similarity between host and parasite genes to test molecular mimicry. The nematomorph's transcriptome changed during host manipulation, whereas no distinct changes were found in mantids. We then discovered numerous possible host-derived genes in nematomorphs, and these genes were frequently up-regulated during host manipulation. Our findings suggest a possible general role of horizontal gene transfer (HGT) in the molecular mechanisms of host manipulation, as well as in the genome evolution of manipulative parasites. The evidence of HGT between multicellular eukaryotes remains scarce but is increasing and, therefore, elucidating its mechanisms will advance our understanding of the enduring influence of HGT on the evolution of the web of life.

Long live the host! Proteomic analysis reveals possible strategies for parasitic manipulation of its social host.

Parasites with complex life cycles often manipulate the phenotype of their intermediate hosts to increase the probability of transmission to their definitive hosts. Infection with Anomotaenia brevis, a cestode that uses Temnothorax nylanderi ants as intermediate hosts, leads to a multiple-fold extension of host lifespan and to changes in behaviour, morphology and colouration. The mechanisms behind these changes are unknown, as is whether the increased longevity is achieved through parasite manipulation. Here, we demonstrate that the parasite releases proteins into its host with functions that might explain the observed changes. These parasitic proteins make up a substantial portion of the proteome of the hosts' haemolymph, and thioredoxin peroxidase and superoxide dismutase, two antioxidants, exhibited the highest abundances among them. The largest part of the secreted proteins could not be annotated, indicating they are either novel or severely altered during recent coevolution to function in host manipulation. We also detected shifts in the hosts' proteome with infection, in particular an overabundance of vitellogenin-like A in infected ants, a protein that regulates division of labour in Temnothorax ants, which could explain the observed behavioural changes. Our results thus suggest two different strategies that might be employed by this parasite to manipulate its host: secreting proteins with immediate influence on the host's phenotype and altering the host's translational activity. Our findings highlight the intricate molecular interplay required to influence the phenotype of a host and point to potential signalling pathways and genes involved in parasite-host communication.

Balanophora genomes display massively convergent evolution with other extreme holoparasites and provide novel insights into parasite-host interactions.

Parasitic plants have evolved to be subtly or severely dependent on host plants to complete their life cycle. To provide new insights into the biology of parasitic plants in general, we assembled genomes for members of the sandalwood order Santalales, including a stem hemiparasite (Scurrula) and two highly modified root holoparasites (Balanophora) that possess chimaeric host-parasite tubers. Comprehensive genome comparisons reveal that hemiparasitic Scurrula has experienced a relatively minor degree of gene loss compared with autotrophic plants, consistent with its moderate degree of parasitism. Nonetheless, patterns of gene loss appear to be substantially divergent across distantly related lineages of hemiparasites. In contrast, Balanophora has experienced substantial gene loss for the same sets of genes as an independently evolved holoparasite lineage, the endoparasitic Sapria (Malpighiales), and the two holoparasite lineages experienced convergent contraction of large gene families through loss of paralogues. This unprecedented convergence supports the idea that despite their extreme and strikingly divergent life histories and morphology, the evolution of these and other holoparasitic lineages can be shaped by highly predictable modes of genome reduction. We observe substantial evidence of relaxed selection in retained genes for both hemi- and holoparasitic species. Transcriptome data also document unusual and novel interactions between Balanophora and host plants at the host-parasite tuber interface tissues, with evidence of mRNA exchange, substantial and active hormone exchange and immune responses in parasite and host.