Managing Maedi-Visna in sheep flocks infected by Small Ruminant Lentivirus (SRLV) by limiting production losses and reducing transmission could be a more cost-effective approach than eliminating the disease. We quantified the progressive physiological costs of SRLV infection by analysing performance measures as readily accessible criteria of disease severity in data from a five-year study of a Polish flock with bi-annual tests of serostatus. SRLV infection had no effect on ewes’ pre-natal reproductive traits but caused increased mortality and reduced growth rate in their lambs, and reduced ewes’ productive lifespan and output. These performance deficits were negligible soon after seroconversion but increased significantly with time since seroconversion, such that they became substantial in older ewes with longer durations of infection. We argue that infectiousness is likely to be correlated with disease severity and therefore tends to increase with the duration of infection, implying that the oldest animals pose the greatest infection risk. Selective culling of under-performing older ewes would limit ongoing production losses and may also reduce disease transmission in affected sheep flocks.

Immunity to severe malaria is acquired quickly, operates independently of pathogen load, and represents a highly effective form of disease tolerance. The mechanism that underpins tolerance remains unknown. We used a human rechallenge model of falciparum malaria in which healthy adult volunteers were infected three times over a 12 mo period to track the development of disease tolerance in real-time. We found that parasitemia triggered a hardwired innate immune response that led to systemic inflammation, pyrexia, and hallmark symptoms of clinical malaria across the first three infections of life. In contrast, a single infection was sufficient to reprogram T cell activation and reduce the number and diversity of effector cells upon rechallenge. Crucially, this did not silence stem-like memory cells but instead prevented the generation of cytotoxic effectors associated with autoinflammatory disease. Tolerized hosts were thus able to prevent collateral tissue damage in the absence of antiparasite immunity.

ABSTRACT Malaria infections represent an iconic example of developmental synchrony, where periodic fevers can result when the population of parasites develops synchronously within host red blood cells. The level of synchrony appears to vary across individual hosts and across parasite species and strains, variation that—once quantified—can illuminate the ecological and evolutionary drivers of synchrony. Yet current approaches for quantifying synchrony in parasites are either biased by population dynamics or unsuitable when population growth rates vary through time, features ubiquitous to parasite populations in vitro and in vivo . Here we develop an approach to estimate synchrony that accounts for population dynamics, including changing population growth rates, and validate it with simulated time series data encompassing a range of synchrony levels in two different host-parasite systems: malaria infections of mice and human malaria parasites in vitro . This new method accurately quantifies developmental synchrony from per capita growth rates using obtainable abundance data even with realistic sampling noise, without the need to sort parasites into developmental stages. Our approach enables variability in developmental schedules to be disentangled from even extreme variation in population dynamics, providing a comparative metric of developmental synchrony.

Unknown factors regulate mitochondrial U-insertion/deletion (U-indel) RNA editing in procyclic-form (PCF) and bloodstream-form (BSF) T. brucei. This editing, directed by anti-sense gRNAs, creates canonical protein-encoding mRNAs and may developmentally control respiration. Canonical editing by gRNAs that specify protein-encoding mRNA sequences occurs amid massive non-canonical editing of unclear sources and biological significance. We found PCF-specific repression at a major early checkpoint in mRNA ND7, involving helicase KREH2-dependent opposite modulation of canonical and non-canonical ‘terminator’ gRNA utilization. Terminator-programmed editing derails canonical editing and installs proposed repressive structure in 30% of the ND7 transcriptome. BSF-to-PCF differentiation in vitro recreated this negative control. Remarkably, KREH2-RNAi knockdown relieved repression and increased editing progression by reverting canonical/terminator gRNA utilization. ND7 transcripts lacking early terminator-directed editing in PCF exhibited similar negative editing control along the mRNA sequence, suggesting global modulation of gRNA utilization fidelity. The terminator is a ‘moonlighting’ gRNA also associated with mRNA COX3 canonical editing, so the gRNA transcriptome seems multifunctional. Thus, KREH2 is the first identified repressor in developmental editing control. This and our prior work support a model whereby KREH2 activates or represses editing in a stage and substrate-specific manner. KREH2’s novel dual role tunes mitochondrial gene expression in either direction during development.

The relative importance of maternal and horizontal transmission of small ruminant lentivirus (SRLV), the causative organism in maedi-visna, is poorly understood. Review of the literature shows that maternal transmission is inefficient, infecting only about 10-25 % of the lambs of infected ewes. Theory proves that maternal transmission alone cannot achieve the rates of transmission that would be required to start or maintain an outbreak. Maternal and horizontal transmission are additive in effect, and we use modelling to show that maternal transmission does not amplify or enhance prevalence in the presence of horizontal transmission. Taking steps to avoid maternal transmission by rearing lambs without infected maternal colostrum does have a role in producing a clean flock, but has no significance for the control of a disease outbreak if the conditions for horizontal transmission are present. Efforts to prevent disease by reducing the spread of SRLV must be focussed on minimising horizontal transmission.

U-insertion/deletion (U-indel) RNA editing in trypanosome mitochondria is directed by guide RNAs (gRNAs). This editing may developmentally control respiration in bloodstream forms (BSF) and insect procyclic forms (PCF). Holo-editosomes include the accessory RNA Editing Substrate Binding Complex (RESC) and RNA Editing Helicase 2 Complex (REH2C), but the specific proteins controlling differential editing remain unknown. Also, RNA editing appears highly error prone because most U-indels do not match the canonical pattern. However, despite extensive non-canonical editing of unknown functions, accurate canonical editing is required for normal cell growth. In PCF, REH2C controls editing fidelity in RESC-bound mRNAs. Here, we report that KREH2, a REH2C-associated helicase, developmentally controls programmed non-canonical editing, including an abundant 3' element in ATPase subunit 6 (A6) mRNA. The 3' element sequence is directed by a proposed novel regulatory gRNA. In PCF, KREH2 RNAi-knockdown up-regulates the 3' element, which establishes a stable structure hindering element removal by canonical initiator-gRNA-directed editing. In BSF, KREH2-knockdown does not up-regulate the 3' element but reduces its high abundance. Thus, KREH2 differentially controls extensive non-canonical editing and associated RNA structure via a novel regulatory gRNA, potentially hijacking factors as a 'molecular sponge'. Furthermore, this gRNA is bifunctional, serving in canonical CR4 mRNA editing whilst installing a structural element in A6 mRNA.

The development of Trypanosoma brucei in its mammalian host is marked by a distinct morphological change as replicative “slender” forms differentiate into cell cycle arrested “stumpy” forms in a quorum-sensing-dependent manner. Although stumpy forms dominate chronic infections at the population level, the proportion of replicative parasites at the individual cell level and the irreversibility of arrest in the bloodstream are unclear. Here, we experimentally demonstrate that developmental cell cycle arrest is definitively irreversible in acute and chronic infections in mice. Furthermore, analysis of replicative capacity and single-cell transcriptome profiling reveal a temporal hierarchy, whereby cell cycle arrest and appearance of a reversible stumpy-like transcriptome precede irreversible commitment and morphological change. Unexpectedly, we show that proliferating parasites are exceptionally scarce in the blood after infections are established. This challenges the ability of bloodstream trypanosomes to sustain infection by proliferation or antigenic variation, these parasites instead being overwhelmingly adapted for transmission.

The World Health Organization targeted Trypanosoma brucei gambiense (Tbg) human African trypanosomiasis for elimination of transmission by 2030. Sensitive molecular markers that specifically detect Tbg type 1 (Tbg1) parasites will be important tools to assist in reaching this goal. We aim at improving molecular diagnosis of Tbg1 infections by targeting the abundant mitochondrial minicircles within the kinetoplast of these parasites. Using Next-Generation Sequencing of total cellular DNA extracts, we assembled and annotated the kinetoplast genome and investigated minicircle sequence diversity in 38 animal- and human-infective trypanosome strains. Computational analyses recognized a total of 241 Minicircle Sequence Classes as Tbg1-specific, of which three were shared by the 18 studied Tbg1 strains. We developed a minicircle-based assay that is applicable on animals and as specific as the TgsGP-based assay, the current golden standard for molecular detection of Tbg1. The median copy number of the targeted minicircle was equal to eight, suggesting our minicircle-based assay may be used for the sensitive detection of Tbg1 parasites. Annotation of the targeted minicircle sequence indicated that it encodes genes essential for the survival of the parasite and will thus likely be preserved in natural Tbg1 populations, the latter ensuring the reliability of our novel diagnostic assay.

Mitochondrial DNA of protists of order Kinetoplastida comprises thousands of interlinked circular molecules arranged in a network. There are two types of molecules called minicircles and maxicircles. Minicircles encode guide RNA (gRNA) genes whose transcripts mediate post-transcriptional editing of maxicircle encoded genes. Minicircles are diverse. The human sleeping sickness parasite Trypanosoma brucei has one of the most diverse sets of minicircle classes of all studied trypanosomatids with hundreds of different classes, each encoding one to four genes mainly within cassettes framed by 18 bp inverted repeats. A third of cassettes have no identifiable gRNA genes even though their sequence structures are similar to cassettes with identifiable genes. Only recently have almost all minicircle classes for some subspecies and isolates of T. brucei been sequenced and annotated with corresponding verification of gRNA expression by small-RNA transcriptome data. These data sets provide a rich resource for understanding the structure of minicircle classes, cassettes and gRNA genes and their transcription. Here, we provide a statistical description of the functionality, expression status, structure and sequence of gRNA genes in a differentiation-competent, laboratory-adapted strain of T. brucei We obtain a clearer definition of what is a gRNA gene. Our analysis supports the idea that many, if not all, cassettes without an identifiable gRNA gene contain decaying remnants of once functional gRNA genes. Finally, we report several new, unexplained discoveries such as the association between cassette position on the minicircle and gene expression and functionality, and the association between gene initiation sequence and anchor position.

Immunity to severe malaria is acquired quickly, operates independently of pathogen load and represents a highly effective form of disease tolerance. The mechanism that underpins tolerance in human malaria remains unknown. We developed a re-challenge model of falciparum malaria in which healthy naive adult volunteers were infected three times over a 12-month period to track the development of disease tolerance in real-time. We found that parasitaemia triggered a hardwired emergency myeloid response that led to systemic inflammation, pyrexia and hallmark symptoms of clinical malaria across the first three infections of life. In contrast, CD4 + T cell activation was quickly modified to reduce the number and diversity of effector cells upon re-challenge. Crucially, this did not silence critical helper T cell functions but instead prevented the generation of cytotoxic effectors associated with autoinflammatory disease. Tolerised hosts were thus able to prevent collateral tissue damage and injury. Host control of T cell activation can therefore be established after a single infection and in the absence of anti-parasite immunity. And furthermore, this rapid host adaptation can protect vital organs to minimise the harm caused by systemic inflammation and sequestration.

Persistent pathogens have evolved to avoid elimination by the mammalian immune system including mechanisms to evade complement. Infections with African trypanosomes can persist for years and cause human and animal disease throughout sub-Saharan Africa. It is not known how trypanosomes limit the action of the alternative complement pathway. Here we identify an African trypanosome receptor for mammalian factor H, a negative regulator of the alternative pathway. Structural studies show how the receptor binds ligand, leaving inhibitory domains of factor H free to inactivate complement C3b deposited on the trypanosome surface. Receptor expression is highest in developmental stages transmitted to the tsetse fly vector and those exposed to blood meals in the tsetse gut. Receptor gene deletion reduced tsetse infection, identifying this receptor as a virulence factor for transmission. This demonstrates how a pathogen evolved a molecular mechanism to increase transmission to an insect vector by exploitation of a mammalian complement regulator.

Maedi-visna (MV) is a complex lentiviral disease syndrome characterised by long immunological and clinical latencies and chronic progressive inflammatory pathology. Incurable at the individual level, it is widespread in most sheep-keeping countries, and is a cause of lost production and poor animal welfare. Culling seropositive animals is the main means of control, but it might be possible to manage virus transmission effectively if its epidemiology was better quantified. We derive a mathematical epidemiological model of the temporal distributions of seroconversion probabilities and estimate susceptibility, transmission rate and latencies in three serological datasets. We demonstrate the existence of epidemiological latency, which has not explicitly been recognised in the SRLV literaure. This time delay between infection and infectiousness apparently exceeds the delay between infection and seroconversion. Poor body condition was associated with more rapid seroconversion, but not with a higher probability of infection. We estimate transmission rates amongst housed sheep to be at about 1,000 times faster than when sheep were at grass, when transmission was negligible. Maternal transmission has only a small role in transmission, because lambs from infected ewes have a low probability of being infected directly by them, and only a small proportion of lambs need be retained to maintain flock size. Our results show that MV is overwhelmingly a disease of housing, where sheep are kept in close proximity. Prevalence of MV is likely to double each year from an initial low incidence in housed flocks penned in typically-sized groups of sheep (c. 50) for even a few days per year. Ewes kept entirely at grass are unlikely to experience transmission frequently enough for MV to persist, and pre-existing infection should die out as older ewes are replaced, thereby essentially curing the flock.

The tropical Andes are an important natural laboratory to understand speciation in many taxa. Here we examined the evolutionary history of parasites of the Leishmania braziliensis species complex based on whole-genome sequencing of 67 isolates from 47 localities in Peru. We first show the origin of Andean Leishmania as a clade of near-clonal lineages that diverged from admixed Amazonian ancestors, accompanied by a significant reduction in genome diversity and large structural variations implicated in host-parasite interactions. Within the Andean species, patterns of population structure were strongly associated with biogeographical origin. Molecular clock and ecological niche modeling suggested that the history of diversification of the Andean lineages is limited to the Late Pleistocene and intimately associated with habitat contractions driven by climate change. These results suggest that changes in forestation over the past 150,000 y have influenced speciation and diversity of these Neotropical parasites. Second, genome-scale analyses provided evidence of meiotic-like recombination between Andean and Amazonian Leishmania species, resulting in full-genome hybrids. The mitochondrial genome of these hybrids consisted of homogeneous uniparental maxicircles, but minicircles originated from both parental species. We further show that mitochondrial minicircles-but not maxicircles-show a similar evolutionary pattern to the nuclear genome, suggesting that compatibility between nuclear-encoded mitochondrial genes and minicircle-encoded guide RNA genes is essential to maintain efficient respiration. By comparing full nuclear and mitochondrial genome ancestries, our data expand our appreciation on the genetic consequences of diversification and hybridization in parasitic protozoa.

additional data associated with Nucleic Acids Research submission by Cooper et al.

Malaria infection is often accompanied by periodic fevers, triggered by synchronous cycles of parasite replication within the host. The degree of synchrony in parasite development influences the efficacy of drugs and immune defenses and is therefore relevant to host health and infectiousness. Synchrony is thought to vary over the course of infection and across different host-parasite genotype or species combinations, but the evolutionary significance - if any - of this diversity remains elusive. Standardized methods are lacking, but the most common metric for quantifying synchrony is the percentage of parasites in a particular developmental stage. We use a heuristic model to show that this metric is often unacceptably biased. Methodological challenges must be addressed to characterize diverse patterns of synchrony and their consequences for disease severity and spread.

Kinetoplastids are protists defined by one of the most complex mitochondrial genomes in nature, the kinetoplast. In the sleeping sickness parasite Trypanosoma brucei, the kinetoplast is a chain mail-like network of two types of interlocked DNA molecules: a few dozen ∼23-kb maxicircles (homologs of the mitochondrial genome of other eukaryotes) and thousands of ∼1-kb minicircles. Maxicircles encode components of respiratory chain complexes and the mitoribosome. Several maxicircle-encoded mRNAs undergo extensive post-transcriptional RNA editing via addition and deletion of uridines. The process is mediated by hundreds of species of minicircle-encoded guide RNAs (gRNAs), but the precise number of minicircle classes and gRNA genes was unknown. Here we present the first essentially complete assembly and annotation of the kinetoplast genome of T. brucei. We have identified 391 minicircles, encoding not only ∼930 predicted 'canonical' gRNA genes that cover nearly all known editing events (accessible via the web at http://hank.bio.ed.ac.uk), but also ∼370 'non-canonical' gRNA genes of unknown function. Small RNA transcriptome data confirmed expression of the majority of both categories of gRNAs. Finally, we have used our data set to refine definitions for minicircle structure and to explore dynamics of minicircle copy numbers.

The ability of malaria (Plasmodium) parasites to adjust investment into sexual transmission stages versus asexually replicating stages is well known, but plasticity in other traits underpinning the replication rate of asexual stages in the blood has received less attention. Such traits include burst size (the number of merozoites produced per schizont), the duration of the asexual cycle, and invasion preference for different ages of red blood cell (RBC). Here, plasticity [environment (E) effects] and genetic variation [genotype (G) effects] in traits relating to asexual replication rate are examined for 4 genotypes of the rodent malaria parasite Plasmodium chabaudi. An experiment tested whether asexual dynamics differ between parasites infecting control versus anaemic hosts, and whether variation in replication rate can be explained by differences in burst size, asexual cycle, and invasion rates. The within-host environment affected each trait to different extents but generally had similar impacts across genotypes. The dynamics of asexual densities exhibited a genotype by environment effect (G×E), in which one of the genotypes increased replication rate more than the others in anaemic hosts. Burst size and cycle duration varied between the genotypes (G), while burst size increased and cycle duration became longer in anaemic hosts (E). Variation in invasion rates of differently aged RBCs was not explained by environmental or genetic effects. Plasticity in burst size and genotype are the only traits making significant contributions to the increase in asexual densities observed in anaemic hosts, together explaining 46.4% of the variation in replication rate. That host anaemia induces several species of malaria parasites to alter conversion rate is well documented. Here, previously unknown plasticity in other traits underpinning asexual replication is revealed. These findings contribute to mounting evidence that malaria parasites deploy a suite of sophisticated strategies to maximize fitness by coping with, or exploiting the opportunities provided by, the variable within-host conditions experienced during infections. That genetic variation and genotype by environment interactions also shape these traits highlights their evolutionary potential. Asexual replication rate is a major determinant of virulence and so, understanding the evolution of virulence requires knowledge of the ecological (within-host environment) and genetic drivers of variation among parasites.

Infection can dramatically alter behavioural and physiological traits as hosts become sick and subsequently return to health. Such “sickness behaviours” include disrupted circadian rhythms in both locomotor activity and body temperature. Host sickness behaviours vary in pathogen species-specific manners but the influence of pathogen intraspecific variation is rarely studied. We examine how infection with the murine malaria parasite, Plasmodium chabaudi , shapes sickness in terms of parasite genotype-specific effects on host circadian rhythms. We reveal that circadian rhythms in host locomotor activity patterns and body temperature become differentially disrupted and in parasite genotype-specific manners. Locomotor activity and body temperature in combination provide more sensitive measures of health than commonly used virulence metrics for malaria (e.g. anaemia). Moreover, patterns of host disruption cannot be explained simply by variation in replication rate across parasite genotypes or the severity of anaemia each parasite genotype causes. It is well known that disruption to circadian rhythms is associated with non-infectious diseases, including cancer, type 2 diabetes, and obesity. Our results reveal that disruption of host circadian rhythms is a genetically variable virulence trait of pathogens with implications for host health and disease tolerance.

The sleeping sickness parasite Trypanosoma brucei has a complex life cycle, alternating between a mammalian host and the tsetse fly vector. A tightly controlled developmental programme ensures parasite transmission between hosts as well as survival within them and involves strict regulation of mitochondrial activities. In the glucose-rich bloodstream, the replicative ‘slender’ stage is thought to produce ATP exclusively via glycolysis and uses the mitochondrial F1FO-ATP synthase as an ATP hydrolysis-driven proton pump to generate the mitochondrial membrane potential (ΔΨm). The ‘procyclic’ stage in the glucose-poor tsetse midgut depends on mitochondrial catabolism of amino acids for energy production, which involves oxidative phosphorylation with ATP production via the F1FO-ATP synthase. Both modes of the F1FO enzyme critically depend on FO subunit a, which is encoded in the parasite’s mitochondrial DNA (kinetoplast or kDNA). Comparatively little is known about mitochondrial function and the role of kDNA in non-replicative ‘stumpy’ bloodstream forms, a developmental stage essential for disease transmission. Here we show that the L262P mutation in the nuclear-encoded F1 subunit γ that permits survival of ‘slender’ bloodstream forms lacking kDNA (‘akinetoplastic’ forms), via FO-independent generation of ΔΨm, also permits their differentiation into stumpy forms. However, these akinetoplastic stumpy cells lack a ΔΨm and have a reduced lifespan in vitro and in mice, which significantly alters the within-host dynamics of the parasite. We further show that generation of ΔΨm in stumpy parasites and their ability to use α-ketoglutarate to sustain viability depend on F1-ATPase activity. Surprisingly, however, loss of ΔΨm does not reduce stumpy life span. We conclude that the L262P γ subunit mutation does not enable FO-independent generation of ΔΨm in stumpy cells, most likely as a consequence of mitochondrial ATP production in these cells. In addition, kDNA-encoded genes other than FO subunit a are important for stumpy form viability.