The genome sequence of Trypanosoma brucei gambiense, causative agent of chronic human african trypanosomiasis

Insights into Peptidyl-Prolyl cis-trans Isomerases from Clinically Important Protozoans: From Structure to Potential Biotechnological Applications

Peptidyl-prolyl cis/trans isomerases (PPIases) are present in a wide variety of microorganisms, including protozoan parasites such as Trypanosoma cruzi, Trypanosoma brucei, Trichomonas vaginalis, Leishmania major, Leishmania donovani, Plasmodium falciparum, Plasmodium vivax, Entamoeba histolytica, Giardia intestinalis, Cryptosporidium parvum, and Cryptosporidium hominis, all of which cause important neglected diseases. PPIases are classified as cyclophilins, FKBPs, or parvulins and play crucial roles in catalyzing the cis-trans isomerization of the peptide bond preceding a proline residue. This activity assists in correct protein folding. However, experimentally, the biological structure-function characterization of PPIases from these protozoan parasites has been poorly addressed. The recombinant production of these enzymes is highly relevant for this ongoing research. Thus, this review explores the structural diversity, functions, recombinant production, activity, and inhibition of protozoan PPIases. We also highlight their potential as biotechnological tools for the in vitro refolding of other recombinant proteins from these parasites. These applications are invaluable for the development of diagnostic and therapeutic tools.

Unveiling the peptidases of parasites from the office chair - The endothelin-converting enzyme case study

The emergence of high-throughput methodologies such as next-generation sequencing and proteomics has necessitated significant advancements in biological databases and bioinformatic tools, therefore reshaping the landscape of research into parasitic peptidases. In this review we outline the development of these resources along the -omics technologies and their transformative impact on the field. Apart from extensive summary of general and specific databases and tools, we provide a general pipeline on how to use these resources effectively to identify candidate peptidases from these large datasets and how to gain as much information about them as possible without leaving the office chair. This pipeline is then applied in an illustrative case study on the endothelin-converting enzyme 1 homologue from Schistosoma mansoni and attempts to highlight the contemporary capabilities of bioinformatics. The case study demonstrate how such approach can aid to hypothesize enzyme functions and interactions through computational analysis alone effectively and emphasizes how such virtual investigations can guide and optimize subsequent wet lab experiments therefore potentially saving precious time and resources. Finally, by showing what can be achieved without traditional wet laboratory methods, this review provides a compelling narrative on the use of bioinformatics to bridge the gap between big data and practical research applications, highlighting the key role of these technologies in furthering our understanding of parasitic diseases.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Tackling Sleeping Sickness: Current and Promising Therapeutics and Treatment Strategies

Human African trypanosomiasis is a neglected tropical disease caused by the extracellular protozoan parasite Trypanosoma brucei, and targeted for eradication by 2030. The COVID-19 pandemic contributed to the lengthening of the proposed time frame for eliminating human African trypanosomiasis as control programs were interrupted. Armed with extensive antigenic variation and the depletion of the B cell population during an infectious cycle, attempts to develop a vaccine have remained unachievable. With the absence of a vaccine, control of the disease has relied heavily on intensive screening measures and the use of drugs. The chemotherapeutics previously available for disease management were plagued by issues such as toxicity, resistance, and difficulty in administration. The approval of the latest and first oral drug, fexinidazole, is a major chemotherapeutic achievement for the treatment of human African trypanosomiasis in the past few decades. Timely and accurate diagnosis is essential for effective treatment, while poor compliance and resistance remain outstanding challenges. Drug discovery is on-going, and herein we review the recent advances in anti-trypanosomal drug discovery, including novel potential drug targets. The numerous challenges associated with disease eradication will also be addressed.

VSGs Expressed during Natural T. b. gambiense Infection Exhibit Extensive Sequence Divergence and a Subspecies-Specific Bias towards Type B N-Terminal Domains

Trypanosoma brucei gambiense is the primary causative agent of human African trypanosomiasis (HAT), a vector-borne disease endemic to West and Central Africa. The extracellular parasite evades antibody recognition within the host bloodstream by altering its variant surface glycoprotein (VSG) coat through a process of antigenic variation. The serological tests that are widely used to screen for HAT use VSG as one of the target antigens. However, the VSGs expressed during human infection have not been characterized. Here, we use VSG sequencing (VSG-seq) to analyze the VSGs expressed in the blood of patients infected with T. b. gambiense and compared them to VSG expression in Trypanosoma brucei rhodesiense infections in humans as well as Trypanosoma brucei brucei infections in mice. The 44 VSGs expressed during T. b. gambiense infection revealed a striking bias toward expression of type B N termini (82% of detected VSGs). This bias is specific to T. b. gambiense, which is unique among T. brucei subspecies in its chronic clinical presentation and anthroponotic nature. The expressed T. b. gambiense VSGs also share very little similarity to sequences from 36 T. b. gambiense whole-genome sequencing data sets, particularly in areas of the VSG protein exposed to host antibodies, suggesting the antigen repertoire is under strong selective pressure to diversify. Overall, this work demonstrates new features of antigenic variation in T. brucei gambiense and highlights the importance of understanding VSG repertoires in nature. IMPORTANCE Human African trypanosomiasis is a neglected tropical disease primarily caused by the extracellular parasite Trypanosoma brucei gambiense. To avoid elimination by the host, these parasites repeatedly replace their variant surface glycoprotein (VSG) coat. Despite the important role of VSGs in prolonging infection, VSG expression during human infections is poorly understood. A better understanding of natural VSG gene expression dynamics can clarify the mechanisms that T. brucei uses to alter its VSG coat. We analyzed the expressed VSGs detected in the blood of patients with trypanosomiasis. Our findings indicate that there are features of antigenic variation unique to human-infective T. brucei subspecies and that natural VSG repertoires may vary more than previously expected.

Genome-Wide Screening for Pathogenic Proteins and microRNAs Associated with Parasite-Host Interactions in Trypanosoma brucei

Tsetse flies are a type of blood-sucking insect living in diverse locations in sub-Saharan Africa. These insects can transmit the unicellular parasite Trypanosoma brucei (T. brucei) which causes African trypanosomiasis in mammals. There remain huge unmet needs for prevention, early detection, and effective treatments for this disease. Currently, few studies have investigated the molecular mechanisms of parasite-host interactions underlying African trypanosomiasis, mainly due to a lack of understanding of the T. brucei genome. In this study, we dissected the genomic and transcriptomic profiles of T. brucei by annotating the genome and analyzing the gene expression. We found about 5% of T. brucei proteins in the human proteome, while more than 80% of T. brucei protein in other trypanosomes. Sequence alignment analysis showed that 142 protein homologs were shared among T. brucei and mammalian genomes. We identified several novel proteins with pathogenic potential supported by their molecular functions in T. brucei, including 24 RNA-binding proteins and six variant surface glycoproteins. In addition, 26 novel microRNAs were characterized, among which five miRNAs were not found in the mammalian genomes. Topology analysis of the miRNA-gene network revealed three genes (RPS27A, UBA52 and GAPDH) involved in the regulation of critical pathways related to the development of African trypanosomiasis. In conclusion, our work opens a new door to understanding the parasite-host interaction mechanisms by resolving the genome and transcriptome of T. brucei.

Fold-change of chromatin condensation in yeast is a conserved property

During mitosis, chromatin is condensed and organized into mitotic chromosomes. Condensation is critical for genome stability and dynamics, yet the degree of condensation is significantly different between multicellular and single-cell eukaryotes. What is less clear is whether there is a minimum degree of chromosome condensation in unicellular eukaryotes. Here, we exploited two-photon microscopy to analyze chromatin condensation in live and fixed cells, enabling studies of some organisms that are not readily amenable to genetic modification. This includes the yeasts Saccharomyces cerevisiae, Schizosaccharomyces pombe, Kluyveromyces lactis, and Candida albicans, as well as a protist Trypanosoma brucei. We found that mitotic chromosomes in this range of species are condensed about 1.5-fold relative to interphase chromatin. In addition, we used two-photon microscopy to reveal that chromatin reorganization in interphase human hepatoma cells infected by the hepatitis C virus is decondensed compared to uninfected cells, which correlates with the previously reported viral-induced changes in chromatin dynamics. This work demonstrates the power of two-photon microscopy to analyze chromatin in a broad range of cell types and conditions, including non-model single-cell eukaryotes. We suggest that similar condensation levels are an evolutionarily conserved property in unicellular eukaryotes and important for proper chromosome segregation. Furthermore, this provides new insights into the process of chromatin condensation during mitosis in unicellular organisms as well as the response of human cells to viral infection.

In silico analysis of the HSP90 chaperone system from the African trypanosome, Trypanosoma brucei

African trypanosomiasis is a neglected tropical disease caused by Trypanosoma brucei (T. brucei) and spread by the tsetse fly in sub-Saharan Africa. The trypanosome relies on heat shock proteins for survival in the insect vector and mammalian host. Heat shock protein 90 (HSP90) plays a crucial role in the stress response at the cellular level. Inhibition of its interactions with chaperones and co-chaperones is being explored as a potential therapeutic target for numerous diseases. This study provides an in silico overview of HSP90 and its co-chaperones in both T. brucei brucei and T. brucei gambiense in relation to human and other trypanosomal species, including non-parasitic Bodo saltans and the insect infecting Crithidia fasciculata. A structural analysis of T. brucei HSP90 revealed differences in the orientation of the linker and C-terminal domain in comparison to human HSP90. Phylogenetic analysis displayed the T. brucei HSP90 proteins clustering into three distinct groups based on subcellular localizations, namely, cytosol, mitochondria, and endoplasmic reticulum. Syntenic analysis of cytosolic HSP90 genes revealed that T. b. brucei encoded for 10 tandem copies, while T. b. gambiense encoded for three tandem copies; Leishmania major (L. major) had the highest gene copy number with 17 tandem copies. The updated information on HSP90 from recently published proteomics on T. brucei was examined for different life cycle stages and subcellular localizations. The results show a difference between T. b. brucei and T. b. gambiense with T. b. brucei encoding a total of twelve putative HSP90 genes, while T. b. gambiense encodes five HSP90 genes. Eighteen putative co-chaperones were identified with one notable absence being cell division cycle 37 (Cdc37). These results provide an updated framework on approaching HSP90 and its interactions as drug targets in the African trypanosome.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

Transcriptomic profiling of Trypanosoma congolense mouthpart parasites from naturally infected flies

Background

Animal African trypanosomiasis, or nagana, is a veterinary disease caused by African trypanosomes transmitted by tsetse flies. In Africa, Trypanosoma congolense is one of the most pathogenic and prevalent causes of nagana in livestock, resulting in high animal morbidity and mortality and extensive production losses. In the tsetse fly, parasites colonise the midgut and eventually reach the mouthparts, from where they can be transmitted as the fly feeds on vertebrate hosts such as cattle. Despite the extreme importance of mouthpart-form parasites for disease transmission, very few global expression profile studies have been conducted in these parasite forms.

Methods

Here, we collected tsetse flies from the Shimba Hills National Reserve, a wildlife area in southeast Kenya, diagnosed T. congolense infections, and sequenced the transcriptomes of the T. congolense parasites colonising the mouthparts of the flies.

Results

We found little correlation between mouthpart parasites from natural and experimental fly infections. Furthermore, we performed differential gene expression analysis between mouthpart and bloodstream parasite forms and identified several surface-expressed genes and 152 novel hypothetical proteins differentially expressed in mouthpart parasites. Finally, we profiled variant antigen expression and observed that a variant surface glycoprotein (VSG) transcript belonging to T. congolense phylotype 8 (i.e. TcIL3000.A.H_000381200), previously observed to be enriched in metacyclic transcriptomes, was present in all wild-caught mouthpart samples as well as bloodstream-form parasites, suggestive of constitutive expression.

Conclusion

Our study provides transcriptomes of trypanosome parasites from naturally infected tsetse flies and suggests that a phylotype 8 VSG gene is constitutively expressed in metacyclic- and bloodstream-form parasites at the population level.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05258-y.

An Update on African Trypanocide Pharmaceutics and Resistance

African trypanosomiasis is associated with Trypanosoma evansi, T. vivax, T. congolense, and T. brucei pathogens in African animal trypanosomiasis (AAT) while T. b gambiense and T. b rhodesiense are responsible for chronic and acute human African trypanosomiasis (HAT), respectively. Suramin sodium suppresses ATP generation during the glycolytic pathway and is ineffective against T. vivax and T. congolense infections. Resistance to suramin is associated with pathogen altered transport proteins. Melarsoprol binds irreversibly with pyruvate kinase protein sulfhydryl groups and neutralizes enzymes which interrupts the trypanosome ATP generation. Melarsoprol resistance is associated with the adenine-adenosine transporter, P2, due to point mutations within this transporter. Eflornithine is used in combination with nifurtimox. Resistance to eflornithine is caused by the deletion or mutation of TbAAT6 gene which encodes the transmembrane amino acid transporter that delivers eflornithine into the cell, thus loss of transporter protein results in eflornithine resistance. Nifurtimox alone is regarded as a poor trypanocide, however, it is effective in melarsoprol-resistant gHAT patients. Resistance is associated with loss of a single copy of the genes encoding for nitroreductase enzymes. Fexinidazole is recommended for first-stage and non-severe second-stage illnesses in gHAT and resistance is associated with trypanosome bacterial nitroreductases which reduce fexinidazole. In AAT, quinapyramine sulfate interferes with DNA synthesis and suppression of cytoplasmic ribosomal activity in the mitochondria. Quinapyramine sulfate resistance is due to variations in the potential of the parasite's mitochondrial membrane. Pentamidines create cross-links between two adenines at 4–5 pairs apart in adenine-thymine-rich portions of Trypanosoma DNA. It also suppresses type II topoisomerase in the mitochondria of Trypanosoma parasites. Pentamidine resistance is due to loss of mitochondria transport proteins P2 and HAPT1. Diamidines are most effective against Trypanosome brucei group and act via the P2/TbAT1 transporters. Diminazene aceturate resistance is due to mutations that alter the activity of P2, TeDR40 (T. b. evansi). Isometamidium chloride is primarily employed in the early stages of trypanosomiasis and resistance is associated with diminazene resistance. Phenanthridine (homidium bromide, also known as ethidium bromide) acts by a breakdown of the kinetoplast network and homidium resistance is comparable to isometamidium. In humans, the development of resistance and adverse side effects against monotherapies has led to the adoption of nifurtimox-eflornithine combination therapy. Current efforts to develop new prodrug combinations of nifurtimox and eflornithine and nitroimidazole fexinidazole as well as benzoxaborole SCYX-7158 (AN5568) for HAT are in progress while little comparable progress has been done for the development of novel therapies to address trypanocide resistance in AAT.

Single nucleotide polymorphisms and copy-number variations in the Trypanosoma brucei repeat (TBR) sequence can be used to enhance amplification and genotyping of Trypanozoon strains

The Trypanosoma brucei repeat (TBR) is a tandem repeat sequence present on the Trypanozoon minichromosomes. Here, we report that the TBR sequence is not as homogenous as previously believed. BLAST analysis of the available T. brucei genomes reveals various TBR sequences of 177 bp and 176 bp in length, which can be sorted into two TBR groups based on a few key single nucleotide polymorphisms. Conventional and quantitative PCR with primers matched to consensus sequences that target either TBR group show substantial copy-number variations in the TBR repertoire within a collection of 77 Trypanozoon strains. We developed the qTBR, a novel PCR consisting of three primers and two probes, to simultaneously amplify target sequences from each of the two TBR groups into one single qPCR reaction. This dual probe setup offers increased analytical sensitivity for the molecular detection of all Trypanozoon taxa, in particular for T.b. gambiense and T. evansi, when compared to existing TBR PCRs. By combining the qTBR with 18S rDNA amplification as an internal standard, the relative copy-number of each TBR target sequence can be calculated and plotted, allowing for further classification of strains into TBR genotypes associated with East, West or Central Africa. Thus, the qTBR takes advantage of the single-nucleotide polymorphisms and copy number variations in the TBR sequences to enhance amplification and genotyping of all Trypanozoon strains, making it a promising tool for prevalence studies of African trypanosomiasis in both humans and animals.

Allele-specific assembly of a eukaryotic genome corrects apparent frameshifts and reveals a lack of nonsense-mediated mRNA decay

To date, most reference genomes represent a mosaic consensus sequence in which the homologous chromosomes are collapsed into one sequence. This approach produces sequence artefacts and impedes analyses of allele-specific mechanisms. Here, we report an allele-specific genome assembly of the diploid parasite Trypanosoma brucei and reveal allelic variants affecting gene expression. Using long-read sequencing and chromosome conformation capture data, we could assign 99.5% of all heterozygote variants to a specific homologous chromosome and build a 66 Mb long allele-specific genome assembly. The phasing of haplotypes allowed us to resolve hundreds of artefacts present in the previous mosaic consensus assembly. In addition, it revealed allelic recombination events, visible as regions of low allelic heterozygosity, enabling the lineage tracing of T. brucei isolates. Interestingly, analyses of transcriptome and translatome data of genes with allele-specific premature termination codons point to the absence of a nonsense-mediated decay mechanism in trypanosomes. Taken together, this study delivers a reference quality allele-specific genome assembly of T. brucei and demonstrates the importance of such assemblies for the study of gene expression control. We expect the new genome assembly will increase the awareness of allele-specific phenomena and provide a platform to investigate them.

Genomics of Trypanosomatidae: Where We Stand and What Needs to Be Done?

Trypanosomatids are easy to cultivate and they are (in many cases) amenable to genetic manipulation. Genome sequencing has become a standard tool routinely used in the study of these flagellates. In this review, we summarize the current state of the field and our vision of what needs to be done in order to achieve a more comprehensive picture of trypanosomatid evolution. This will also help to illuminate the lineage-specific proteins and pathways, which can be used as potential targets in treating diseases caused by these parasites.

Evolution of the variant surface glycoprotein family in African trypanosomes

An intriguing and remarkable feature of African trypanosomes is their antigenic variation system, mediated by the variant surface glycoprotein (VSG) family and fundamental to both immune evasion and disease epidemiology within host populations. Recent studies have revealed that the VSG repertoire has a complex evolutionary history. Sequence diversity, genomic organization, and expression patterns are species-specific, which may explain other variations in parasite virulence and disease pathology. Evidence also shows that we may be underestimating the extent to what VSGs are repurposed beyond their roles as variant antigens, establishing a need to examine VSG functionality more deeply. Here, we review sequence variation within the VSG gene family, and highlight the many opportunities to explore their likely diverse contributions to parasite survival.

Molecular detection of Trypanosoma evansi in dogs from India and Southeast Asia

Trypanosoma evansi, the causative agent of surra, is a hemoflagellate protozoan mechanically transmitted by hematophagous flies, mainly in tropical and subtropical regions. This protozoan affects several mammalian hosts, including dogs, which are highly susceptible to the infection. To investigate the occurrence of T. evansi in dogs, a total of 672 DNA samples from India (n = 228), Indonesia (n = 57), Malaysia (n = 45), the Philippines (n = 103), Thailand (n = 120), and Vietnam (n = 119) were screened by using species-specific conventional PCR. Of the tested dogs, 10 (1.5%) scored positive to T. evansi. In particular, positive samples were detected in canine blood samples collected from India (n = 4; 1.8%), Indonesia (n = 4; 7%), and Malaysia (n = 2; 4.4%). All tested samples from the Philippines, Thailand and Vietnam were negative. Nucleotide sequence analysis revealed a high variation (i.e. from 0.4% to 6.2%) among the RoTat 1.2 variant surface glycoprotein (vsg) gene. Although the number of sequences included in this analysis is relatively small, this nucleotide variation may indicate the divergence of T. evansi RoTat 1.2 vsg gene among different strains. The high incidence of T. evansi previously reported in cattle and buffaloes in India and Southeast Asia suggests that these animals are the main source of infection to dogs.

Molecular identification of Trypanosoma brucei gambiense in naturally infected pigs, dogs and small ruminants confirms domestic animals as potential reservoirs for sleeping sickness in Chad

Human African trypanosomiasis (HAT) has been targeted for zero transmission to humans by 2030. Animal reservoirs of gambiense-HAT could jeopardize these elimination goals. This study was undertaken to identify potential host reservoirs for Trypanosoma brucei gambiense by detecting its natural infections in domestic animals of Chadian HAT foci. Blood samples were collected from 267 goats, 181 sheep, 154 dogs, and 67 pigs. Rapid diagnostic test (RDT) and capillary tube centrifugation (CTC) were performed to search for trypanosomes. DNA was extracted from the buffy coat, and trypanosomes of the subgenus Trypanozoon as well as T. b. gambiense were identified by PCR. Of 669 blood samples, 19.4% were positive by RDT and 9.0% by CTC. PCR revealed 150 animals (22.4%) with trypanosomes belonging to Trypanozoon, including 18 (12%) T. b. gambiense. This trypanosome was found in all investigated animal species and all HAT foci. Between animal species or villages, no significant differences were observed in the number of animals harboring T. b. gambiense DNA. Pigs, dogs, sheep and goats appeared to be potential reservoir hosts for T. b. gambiense in Chad. The identification of T. b. gambiense in all animal species of all HAT foci suggests that these animals should be considered when designing new control strategies for sustainable elimination of HAT. Investigations aiming to decrypt their specific role in each epidemiological setting are important to achieve zero transmission of HAT.

Draft Genome Assembly of the Freshwater Apex Predator Wels Catfish (Silurus glanis) Using Linked-Read Sequencing

The wels catfish (Silurus glanis) is one of the largest freshwater fish species in the world. This top predator plays a key role in ecosystem stability, and represents an iconic trophy-fish for recreational fishermen. S. glanis is also a highly valued species for its high-quality boneless flesh, and has been cultivated for over 100 years in Eastern and Central Europe. The interest in rearing S. glanis continues to grow; the aquaculture production of this species has almost doubled during the last decade. However, despite its high ecological, cultural and economic importance, the available genomic resources for S. glanis are very limited. To fulfill this gap we report a de novo assembly and annotation of the whole genome sequence of a female S. glanis. The linked-read based technology with 10X Genomics Chromium chemistry and Supernova assembler produced a highly continuous draft genome of S. glanis: ∼0.8Gb assembly (scaffold N₅₀ = 3.2 Mb; longest individual scaffold = 13.9 Mb; BUSCO completeness = 84.2%), which included 313.3 Mb of putative repeated sequences. In total, 21,316 protein-coding genes were predicted, of which 96% were annotated functionally from either sequence homology or protein signature searches. The highly continuous genome assembly will be an invaluable resource for aquaculture genomics, genetics, conservation, and breeding research of S. glanis.

Trypanosome infections in naturally infected horses and donkeys of three active sleeping sickness foci in the south of Chad

Background

Equine trypanosomiases are complex infectious diseases with overlapping clinical signs defined by their mode of transmission. Despite their economic impacts, these diseases have been neglected by the scientific community, the veterinary authorities and regulatory organizations. To fill the observed knowledge gap, we undertook the identification of different trypanosome species and subspecies naturally infecting horses and donkeys within the Chadian sleeping sickness focus. The objective of the study was to investigate the potential role of these domestic animals as reservoirs of the human-infective Trypanosoma brucei gambiense.

Method

Blood samples were collected from 155 donkeys and 131 horses in three human African trypanosomiasis (HAT) foci in Chad. Rapid diagnostic test (RDT) and capillary tube centrifugation (CTC) test were used to search for trypanosome infections. DNA was extracted from each blood sample and different trypanosome species and subspecies were identified with molecular tools.

Results

From 286 blood samples collected, 54 (18.9%) and 36 (12.6%) were positive for RDT and CTC, respectively. PCR revealed 101 (35.3%) animals with trypanosome infections. The Cohen’s kappa coefficient used to evaluate the concordance between the diagnostic methods were low; ranging from 0.09 ± 0.05 to 0.48 ± 0.07. Trypanosomes of the subgenus Trypanozoon were the most prevalent (29.4%), followed by T. congolense forest (11.5%), Trypanosoma congolense savannah (4.9%) and Trypanosoma vivax (4.5%). Two donkeys and one horse from the Maro HAT focus were found with T. b. gambiense infections. No significant differences were observed in the infection rates of different trypanosomes between animal species and HAT foci.

Conclusions

This study revealed several trypanosome species and subspecies in donkeys and horses, highlighting the existence of AAT in HAT foci in Chad. The identification of T. b. gambiense in donkeys and horses suggests considering these animals as potential reservoir for HAT in Chad. The presence of both human-infective and human non-infective trypanosomes species highlights the need for developing joint control strategies for HAT and AAT.

P-type transport ATPases in Leishmania and Trypanosoma

P-type ATPases are critical to the maintenance and regulation of cellular ion homeostasis and membrane lipid asymmetry due to their ability to move ions and phospholipids against a concentration gradient by utilizing the energy of ATP hydrolysis. P-type ATPases are particularly relevant in human pathogenic trypanosomatids which are exposed to abrupt and dramatic changes in their external environment during their life cycles. This review describes the complete inventory of ion-motive, P-type ATPase genes in the human pathogenic Trypanosomatidae; eight Leishmania species (L. aethiopica, L. braziliensis, L. donovani, L. infantum, L. major, L. mexicana, L. panamensis, L. tropica), Trypanosoma cruzi and three Trypanosoma brucei subspecies (Trypanosoma brucei brucei TREU927, Trypanosoma brucei Lister strain 427, Trypanosoma brucei gambiense DAL972). The P-type ATPase complement in these trypanosomatids includes the P_1B (metal pumps), P_2A (SERCA, sarcoplasmic-endoplasmic reticulum calcium ATPases), P_2B (PMCA, plasma membrane calcium ATPases), P_2D (Na⁺ pumps), P_3A (H⁺ pumps), P₄ (aminophospholipid translocators), and P_5B (no assigned specificity) subfamilies. These subfamilies represent the P-type ATPase transport functions necessary for survival in the Trypanosomatidae as P-type ATPases for each of these seven subfamilies are found in all Leishmania and Trypanosoma species included in this analysis. These P-type ATPase subfamilies are correlated with current molecular and biochemical knowledge of their function in trypanosomatid growth, adaptation, infectivity, and survival.

Transcriptional and genomic parallels between the monoxenous parasite Herpetomonas muscarum and Leishmania

Trypanosomatid parasites are causative agents of important human and animal diseases such as sleeping sickness and leishmaniasis. Most trypanosomatids are transmitted to their mammalian hosts by insects, often belonging to Diptera (or true flies). These are called dixenous trypanosomatids since they infect two different hosts, in contrast to those that infect just insects (monoxenous). However, it is still unclear whether dixenous and monoxenous trypanosomatids interact similarly with their insect host, as fly-monoxenous trypanosomatid interaction systems are rarely reported and under-studied-despite being common in nature. Here we present the genome of monoxenous trypanosomatid Herpetomonas muscarum and discuss its transcriptome during in vitro culture and during infection of its natural insect host Drosophila melanogaster. The H. muscarum genome is broadly syntenic with that of human parasite Leishmania major. We also found strong similarities between the H. muscarum transcriptome during fruit fly infection, and those of Leishmania during sand fly infections. Overall this suggests Drosophila-Herpetomonas is a suitable model for less accessible insect-trypanosomatid host-parasite systems such as sand fly-Leishmania.

Mutualist-Provisioned Resources Impact Vector Competency

Parasites elicit several physiological changes in their host to enhance transmission. Little is known about the functional association between parasitism and microbiota-provisioned resources typically dedicated to animal hosts and how these goods may be rerouted to optimize parasite development. This study is the first to identify a specific symbiont-generated metabolite that impacts insect vector competence by facilitating parasite establishment and, thus, eventual transmission. Specifically, we demonstrate that the tsetse fly obligate mutualist Wigglesworthia provisions folate (vitamin B₉) that pathogenic African trypanosomes exploit in an effort to successfully establish an infection in the vector’s MG. This process is essential for the parasite to complete its life cycle and be transmitted to a new vertebrate host. Disrupting metabolic contributions provided by the microbiota of arthropod disease vectors may fuel future innovative control strategies while also offering minimal nontarget effects.

Many symbionts supplement their host’s diet with essential nutrients. However, whether these nutrients also enhance parasitism is unknown. In this study, we investigated whether folate (vitamin B₉) production by the tsetse fly (Glossina spp.) essential mutualist, Wigglesworthia, aids auxotrophic African trypanosomes in completing their life cycle within this obligate vector. We show that the expression of Wigglesworthia folate biosynthesis genes changes with the progression of trypanosome infection within tsetse. The disruption of Wigglesworthia folate production caused a reduction in the percentage of flies that housed midgut (MG) trypanosome infections. However, decreased folate did not prevent MG trypanosomes from migrating to and establishing an infection in the fly’s salivary glands, thus suggesting that nutrient requirements vary throughout the trypanosome life cycle. We further substantiated that trypanosomes rely on symbiont-generated folate by feeding this vitamin to Glossina brevipalpis, which exhibits low trypanosome vector competency and houses Wigglesworthia incapable of producing folate. Folate-supplemented G. brevipalpis flies were significantly more susceptible to trypanosome infection, further demonstrating that this vitamin facilitates parasite infection establishment. Our cumulative results provide evidence that Wigglesworthia provides a key metabolite (folate) that is “hijacked” by trypanosomes to enhance their infectivity, thus indirectly impacting tsetse species vector competency. Parasite dependence on symbiont-derived micronutrients, which likely also occurs in other arthropod vectors, represents a relationship that may be exploited to reduce disease transmission.

The Hsp70/J-protein machinery of the African trypanosome, Trypanosoma brucei

The etiological agent of the neglected tropical disease African trypanosomiasis, Trypanosoma brucei, possesses an expanded and diverse repertoire of heat shock proteins, which have been implicated in cytoprotection, differentiation, as well as progression and transmission of the disease. Hsp70 plays a crucial role in proteostasis, and inhibition of its interactions with co-chaperones is emerging as a potential therapeutic target for numerous diseases. In light of genome annotations and the release of the genome sequence of the human infective subspecies, an updated and current in silico overview of the Hsp70/J-protein machinery in both T. brucei brucei and T. brucei gambiense was conducted. Functional, structural, and evolutionary analyses of the T. brucei Hsp70 and J-protein families were performed. The Hsp70 and J-proteins from humans and selected kinetoplastid parasites were used to assist in identifying proteins from T. brucei, as well as the prediction of potential Hsp70-J-protein partnerships. The Hsp70 and J-proteins were mined from numerous genome-wide proteomics studies, which included different lifecycle stages and subcellular localisations. In this study, 12 putative Hsp70 proteins and 67 putative J-proteins were identified to be encoded on the genomes of both T. brucei subspecies. Interestingly there are 6 type III J-proteins that possess tetratricopeptide repeat-containing (TPR) motifs. Overall, it is envisioned that the results of this study will provide a future context for studying the biology of the African trypanosome and evaluating Hsp70 and J-protein interactions as potential drug targets.

Diversity and Evolution of Sensor Histidine Kinases in Eukaryotes

Histidine kinases (HKs) are primary sensor proteins that act in cell signaling pathways generically referred to as "two-component systems" (TCSs). TCSs are among the most widely distributed transduction systems used by both prokaryotic and eukaryotic organisms to detect and respond to a broad range of environmental cues. The structure and distribution of HK proteins are now well documented in prokaryotes, but information is still fragmentary for eukaryotes. Here, we have taken advantage of recent genomic resources to explore the structural diversity and the phylogenetic distribution of HKs in the prominent eukaryotic supergroups. Searches of the genomes of 67 eukaryotic species spread evenly throughout the phylogenetic tree of life identified 748 predicted HK proteins. Independent phylogenetic analyses of predicted HK proteins were carried out for each of the major eukaryotic supergroups. This allowed most of the compiled sequences to be categorized into previously described HK groups. Beyond the phylogenetic analysis of eukaryotic HKs, this study revealed some interesting findings: 1) characterization of some previously undescribed eukaryotic HK groups with predicted functions putatively related to physiological traits; 2) discovery of HK groups that were previously believed to be restricted to a single kingdom in additional supergroups, and 3) indications that some evolutionary paths have led to the appearance, transfer, duplication, and loss of HK genes in some phylogenetic lineages. This study provides an unprecedented overview of the structure and distribution of HKs in the Eukaryota and represents a first step toward deciphering the evolution of TCS signaling in living organisms.

Highly Continuous Genome Assembly of Eurasian Perch (Perca fluviatilis) Using Linked-Read Sequencing

The Eurasian perch (Perca fluviatilis) is the most common fish of the Percidae family and is widely distributed across Eurasia. Perch is a popular target for professional and recreational fisheries, and a promising freshwater aquaculture species in Europe. However, despite its high ecological, economical and societal importance, the available genomic resources for P. fluviatilis are rather limited. In this work, we report de novo assembly and annotation of the whole genome sequence of perch. The linked-read based technology with 10X Genomics Chromium chemistry and Supernova assembler produced a draft perch genome ∼1.0 Gbp assembly (scaffold N₅₀ = 6.3 Mb; the longest individual scaffold of 29.3 Mb; BUSCO completeness of 88.0%), which included 281.6 Mb of putative repeated sequences. The perch genome assembly presented here, generated from small amount of starting material (0.75 ng) and a single linked-read library, is highly continuous and considerably more complete than the currently available draft of P. fluviatilis genome. A total of 23,397 protein-coding genes were predicted, 23,171 (99%) of which were annotated functionally from either sequence homology or protein signature searches. Linked-read technology enables fast, accurate and cost-effective de novo assembly of large non-model eukaryote genomes. The highly continuous assembly of the Eurasian perch genome presented in this study will be an invaluable resource for a range of genetic, ecological, physiological, ecotoxicological, functional and comparative genomic studies in perch and other fish species of the Percidae family.

Evaluation of multiplex PCR assay for detection of Babesia spp, Ehrlichia canis and Trypanosoma evansi in dogs

A multiplex PCR test was evaluated to detect the DNA of three important dog haemoparasites by comparing with singular PCR counterpart on clinical blood samples of dogs in and around Chennai, Tamil Nadu, India. Initial screening of samples was done by microscopic examination of peripheral blood smear and singular PCR and those found exclusively positive for Babesia spp, Ehrlichia canis and Trypanosoma evansi and concurrent infections were used to standardize multiplex PCR. Amplicons of 619 bp, 377 bp and 227 bp corresponding to Babesia spp (18S rRNA gene), E. canis (VirB9 gene), and T.evansi (VSG gene) respectively were amplified, without any non-specific amplification. The laboratory sensitivity (91.7% to 100%) and specificity (100%) of the multiplex PCR were calculated using 'true positive' and 'true negative' dog blood samples obtained in the initial screening process. Clinical blood samples from 287 dogs were screened using singular PCR and multiplex PCR tests for the presence of genome of Babesia spp, E. canis and T. evansi. The multiplex PCR was found to have high level of diagnostic specificity (97.5%-100%) in the detection of all three dog blood parasites and high level of diagnostic sensitivity (95%) in the detection of T. evansi from field level clinical blood samples compared to the singular PCR. However, the diagnostic sensitivity of the multiplex PCR was found to be low to moderate (40.45%-66.7%) in detection of Babesia spp and E. canis from field level clinical blood samples. The strength of agreement between singular and multiplex PCR assays was 'moderate' (0.445), 'good' (0.708) and 'very good' (0.968) in detection of DNA of Babesia spp, E. canis and T. evansi. The multiplex PCR was found to be 10 fold less sensitive in comparison with the singular PCR counterpart.

Blood feeding tsetse flies as hosts and vectors of mammals-pre-adapted African Trypanosoma: current and expected research directions

Research on the zoo-anthropophilic blood feeding tsetse flies' biology conducted, by different teams, in laboratory settings and at the level of the ecosystems- where also co-perpetuate African Trypanosoma- has allowed to unveil and characterize key features of tsetse flies' bacterial symbionts on which rely both (a) the perpetuation of the tsetse fly populations and (b) the completion of the developmental program of the African Trypanosoma. Transcriptomic analyses have already provided much information on tsetse fly genes as well as on genes of the fly symbiotic partners Sodalis glossinidius and Wigglesworthia, which account for the successful onset or not of the African Trypanosoma developmental program. In parallel, identification of the non- symbiotic bacterial communities hosted in the tsetse fly gut has recently been initiated: are briefly introduced those bacteria genera and species common to tsetse flies collected from distinct ecosystems, that could be further studied as potential biologicals preventing the onset of the African Trypanosoma developmental program. Finally, future work will need to concentrate on how to render tsetse flies refractory, and the best means to disseminate them in the field in order to establish an overall refractory fly population.

The Structure of a Conserved Telomeric Region Associated with Variant Antigen Loci in the Blood Parasite Trypanosoma congolense

African trypanosomiasis is a vector-borne disease of humans and livestock caused by African trypanosomes (Trypanosoma spp.). Survival in the vertebrate bloodstream depends on antigenic variation of Variant Surface Glycoproteins (VSGs) coating the parasite surface. In T. brucei, a model for antigenic variation, monoallelic VSG expression originates from dedicated VSG expression sites (VES). Trypanosoma brucei VES have a conserved structure consisting of a telomeric VSG locus downstream of unique, repeat sequences, and an independent promoter. Additional protein-coding sequences, known as “Expression Site Associated Genes (ESAGs)”, are also often present and are implicated in diverse, bloodstream-stage functions. Trypanosoma congolense is a related veterinary pathogen, also displaying VSG-mediated antigenic variation. A T. congolense VES has not been described, making it unclear if regulation of VSG expression is conserved between species. Here, we describe a conserved telomeric region associated with VSG loci from long-read DNA sequencing of two T. congolense strains, which consists of a distal repeat, conserved noncoding elements and other genes besides the VSG; although these are not orthologous to T. brucei ESAGs. Most conserved telomeric regions are associated with accessory minichromosomes, but the same structure may also be associated with megabase chromosomes. We propose that this region represents the T. congolense VES, and through comparison with T. brucei, we discuss the parallel evolution of antigenic switching mechanisms, and unique adaptation of the T. brucei VES for developmental regulation of bloodstream-stage genes. Hence, we provide a basis for understanding antigenic switching in T. congolense and the origins of the African trypanosome VES.

Variant antigen repertoires in Trypanosoma congolense populations and experimental infections can be profiled from deep sequence data using universal protein motifs

African trypanosomes are vector-borne hemoparasites of humans and animals. In the mammal, parasites evade the immune response through antigenic variation. Periodic switching of the variant surface glycoprotein (VSG) coat covering their cell surface allows sequential expansion of serologically distinct parasite clones. Trypanosome genomes contain many hundreds of VSG genes, subject to rapid changes in nucleotide sequence, copy number, and chromosomal position. Thus, analyzing, or even quantifying, VSG diversity over space and time presents an enormous challenge to conventional techniques. Indeed, previous population genomic studies have overlooked this vital aspect of pathogen biology for lack of analytical tools. Here we present a method for analyzing population-scale VSG diversity in Trypanosoma congolense from deep sequencing data. Previously, we suggested that T. congolense VSGs segregate into defined “phylotypes” that do not recombine. In our data set comprising 41 T. congolense genome sequences from across Africa, these phylotypes are universal and exhaustive. Screening sequence contigs with diagnostic protein motifs accurately quantifies relative phylotype frequencies, providing a metric of VSG diversity, called the “variant antigen profile.” We applied our metric to VSG expression in the tsetse fly, showing that certain, rare VSG phylotypes may be preferentially expressed in infective, metacyclic-stage parasites. Hence, variant antigen profiling accurately and rapidly determines the T. congolense VSG gene and transcript repertoire from sequence data, without need for manual curation or highly contiguous sequences. It offers a tractable approach to measuring VSG diversity across strains and during infections, which is imperative to understanding the host–parasite interaction at population and individual scales.

Tsetse fly evolution, genetics and the trypanosomiases - A review

This reviews work published since 2007. Relative efforts devoted to the agents of African trypanosomiasis and their tsetse fly vectors are given by the numbers of PubMed accessions. In the last 10 years PubMed citations number 3457 for Trypanosoma brucei and 769 for Glossina. The development of simple sequence repeats and single nucleotide polymorphisms afford much higher resolution of Glossina and Trypanosoma population structures than heretofore. Even greater resolution is offered by partial and whole genome sequencing. Reproduction in T. brucei sensu lato is principally clonal although genetic recombination in tsetse salivary glands has been demonstrated in T. b. brucei and T. b. rhodesiense but not in T. b. gambiense. In the past decade most genetic attention was given to the chief human African trypanosomiasis vectors in subgenus Nemorhina e.g., Glossina f. fuscipes, G. p. palpalis, and G. p. gambiense. The chief interest in Nemorhina population genetics seemed to be finding vector populations sufficiently isolated to enable efficient and long-lasting suppression. To this end estimates were made of gene flow, derived from F_ST and its analogues, and Ne, the size of a hypothetical population equivalent to that under study. Genetic drift was greater, gene flow and Ne typically lesser in savannah inhabiting tsetse (subgenus Glossina) than in riverine forms (Nemorhina). Population stabilities were examined by sequential sampling and genotypic analysis of nuclear and mitochondrial genomes in both groups and found to be stable. Gene frequencies estimated in sequential samplings differed by drift and allowed estimates of effective population numbers that were greater for Nemorhina spp than Glossina spp. Prospects are examined of genetic methods of vector control. The tsetse long generation time (c. 50 d) is a major contraindication to any suggested genetic method of tsetse population manipulation. Ecological and modelling research convincingly show that conventional methods of targeted insecticide applications and traps/targets can achieve cost-effective reduction in tsetse densities.

A spatial genetics approach to inform vector control of tsetse flies (Glossina fuscipes fuscipes) in Northern Uganda

Tsetse flies (genus Glossina) are the only vector for the parasitic trypanosomes responsible for sleeping sickness and nagana across sub-Saharan Africa. In Uganda, the tsetse fly Glossina fuscipes fuscipes is responsible for transmission of the parasite in 90% of sleeping sickness cases, and co-occurrence of both forms of human-infective trypanosomes makes vector control a priority. We use population genetic data from 38 samples from northern Uganda in a novel methodological pipeline that integrates genetic data, remotely sensed environmental data, and hundreds of field-survey observations. This methodological pipeline identifies isolated habitat by first identifying environmental parameters correlated with genetic differentiation, second, predicting spatial connectivity using field-survey observations and the most predictive environmental parameter(s), and third, overlaying the connectivity surface onto a habitat suitability map. Results from this pipeline indicated that net photosynthesis was the strongest predictor of genetic differentiation in G. f. fuscipes in northern Uganda. The resulting connectivity surface identified a large area of well-connected habitat in northwestern Uganda, and twenty-four isolated patches on the northeastern margin of the G. f. fuscipes distribution. We tested this novel methodological pipeline by completing an ad hoc sample and genetic screen of G. f. fuscipes samples from a model-predicted isolated patch, and evaluated whether the ad hoc sample was in fact as genetically isolated as predicted. Results indicated that genetic isolation of the ad hoc sample was as genetically isolated as predicted, with differentiation well above estimates made in samples from within well-connected habitat separated by similar geographic distances. This work has important practical implications for the control of tsetse and other disease vectors, because it provides a way to identify isolated populations where it will be safer and easier to implement vector control and that should be prioritized as study sites during the development and improvement of vector control methods.

Genome-Wide SNP Analysis Reveals Distinct Origins of Trypanosoma evansi and Trypanosoma equiperdum

Trypanosomes cause a variety of diseases in man and domestic animals in Africa, Latin America, and Asia. In the Trypanozoon subgenus, Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense cause human African trypanosomiasis, whereas Trypanosoma brucei brucei, Trypanosoma evansi, and Trypanosoma equiperdum are responsible for nagana, surra, and dourine in domestic animals, respectively. The genetic relationships between T. evansi and T. equiperdum and other Trypanozoon species remain unclear because the majority of phylogenetic analyses has been based on only a few genes. In this study, we have conducted a phylogenetic analysis based on genome-wide SNP analysis comprising 56 genomes from the Trypanozoon subgenus. Our data reveal that T. equiperdum has emerged at least once in Eastern Africa and T. evansi at two independent occasions in Western Africa. The genomes within the T. equiperdum and T. evansi monophyletic clusters show extremely little variation, probably due to the clonal spread linked to the independence from tsetse flies for their transmission.

UDP-glycosyltransferase genes in trypanosomatid genomes have diversified independently to meet the distinct developmental needs of parasite adaptations

BACKGROUND

Trypanosomatid parasites such as Trypanosoma spp. and Leishmania spp. are a major source of infectious disease in humans and domestic animals worldwide. Fundamental to the host-parasite interactions of these potent pathogens are their cell surfaces, which are highly decorated with glycosylated proteins and other macromolecules. Trypanosomatid genomes contain large multi-copy gene families encoding UDP-dependent glycosyltransferases (UGTs), the primary role of which is cell-surface decoration. Here we report a phylogenetic analysis of UGTs from diverse trypanosomatid genomes, the aim of which was to understand the origin and evolution of their diversity.

RESULTS

By combining phylogenetics with analyses of recombination, and selection, we compared UGT repertoire, genomic context and sequence evolution across 19 trypanosomatids. We identified a UGT lineage present in stercorarian trypanosomes and a free-living kinetoplastid Bodo saltans that likely represents the ancestral state of this gene family. The phylogeny of parasite-specific genes shows that UGTs repertoire in Leishmaniinae and salivarian trypanosomes has expanded independently and with distinct evolutionary dynamics. In the former, the ancestral UGT repertoire was organised in a tandem array from which sporadic transpositions to telomeric regions occurred, allowing expansion most likely through telomeric exchange. In the latter, the ancestral UGT repertoire was comprised of seven subtelomeric lineages, two of which have greatly expanded potentially by gene transposition between these dynamic regions of the genome.

CONCLUSIONS

The phylogeny of UGTs confirms that they represent a substantial parasite-specific innovation, which has diversified independently in the distinct trypanosomatid lineages. Nonetheless, developmental regulation has been a strong driver of UGTs diversification in both African trypanosomes and Leishmania.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

The genomes of Crithidia bombi and C. expoeki, common parasites of bumblebees

Trypanosomatids (Trypanosomatidae, Kinetoplastida) are flagellated protozoa containing many parasites of medical or agricultural importance. Among those, Crithidia bombi and C. expoeki, are common parasites in bumble bees around the world, and phylogenetically close to Leishmania and Leptomonas. They have a simple and direct life cycle with one host, and partially castrate the founding queens greatly reducing their fitness. Here, we report the nuclear genome sequences of one clone of each species, extracted from a field-collected infection. Using a combination of Roche 454 FLX Titanium, Pacific Biosciences PacBio RS, and Illumina GA2 instruments for C. bombi, and PacBio for C. expoeki, we could produce high-quality and well resolved sequences. We find that these genomes are around 32 and 34 MB, with 7,808 and 7,851 annotated genes for C. bombi and C. expoeki, respectively-which is somewhat less than reported from other trypanosomatids, with few introns, and organized in polycistronic units. A large fraction of genes received plausible functional support in comparison primarily with Leishmania and Trypanosoma. Comparing the annotated genes of the two species with those of six other trypanosomatids (C. fasciculata, L. pyrrhocoris, L. seymouri, B. ayalai, L. major, and T. brucei) shows similar gene repertoires and many orthologs. Similar to other trypanosomatids, we also find signs of concerted evolution in genes putatively involved in the interaction with the host, a high degree of synteny between C. bombi and C. expoeki, and considerable overlap with several other species in the set. A total of 86 orthologous gene groups show signatures of positive selection in the branch leading to the two Crithidia under study, mostly of unknown function. As an example, we examined the initiating glycosylation pathway of surface components in C. bombi, finding it deviates from most other eukaryotes and also from other kinetoplastids, which may indicate rapid evolution in the extracellular matrix that is involved in interactions with the host. Bumble bees are important pollinators and Crithidia-infections are suspected to cause substantial selection pressure on their host populations. These newly sequenced genomes provide tools that should help better understand host-parasite interactions in these pollinator pathogens.

Serological tests for gambiense human African trypanosomiasis detect antibodies in cattle

Background

Serological tests for gambiense human African trypanosomiasis (gHAT) detect antibodies to antigens on the cell surface of bloodstream trypanosomes. As trypanosomes that cause animal African trypanosomiasis (AAT) also express related antigens, we have evaluated two rapid diagnostic tests (RDTs) on cattle in trypanosomiasis endemic and non-endemic regions, to determine whether gHAT serological tests could also be used to screen for AAT.

Methods

Two RDTs, 1G RDT, made with native antigens, and p2G RDT, made with recombinant antigens, were tested on 121 cattle in a trypanosomiasis-free region, and on 312 cattle from a rhodesiense HAT and AAT endemic region. A subset of samples from the endemic region were also tested with two immune trypanolysis (TL) tests. The sensitivity of the tests was estimated by evaluating the result of the RDT on samples that were positive by both microscopy and internal transcribed spacer (ITS) PCR, whilst specificity was the result of the RDT on samples that were negative by ITS PCR and microscopy, and others from the non-endemic region.

Results

The specificity of the p2G RDT on cattle from the non-endemic region was 97.5% (95% CI: 93.0–99.2%), compared to only 57.9% (95% CI: 48.9–66.3%) for 1G RDT. The specificities of 1G RDT, p2G RDT and TL on endemic control cattle were 14.6% (95% CI: 9.7–21.5%), 22.6% (95% CI: 16.4–30.3%) and 68.3% (95% CI: 59.6–75.9%), respectively. The sensitivities of the tests on trypanosome positive samples were 85.1% (95% CI: 79.1–89.7%), 89.1% (95% CI: 83.7–93.0%) and 59.3% (95% CI: 51.8–66.4%), respectively. Among the same samples, 51.7% were positive by both TL and the 1G RDT.

Conclusions

These serological tests detect cross-reacting antibodies in cattle. The p2G RDT based on recombinant antigens had a high specificity in a non-endemic region, while the 1G RDT had a lower specificity, suggesting cross-reactivity with other pathogens.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2487-8) contains supplementary material, which is available to authorized users.

Genomic analyses of African Trypanozoon strains to assess evolutionary relationships and identify markers for strain identification

African trypanosomes of the sub-genus Trypanozoon) are eukaryotic parasitesthat cause disease in either humans or livestock. The development of genomic resources can be of great use to those interested in studying and controlling the spread of these trypanosomes. Here we present a large comparative analysis of Trypanozoon whole genomes, 83 in total, including human and animal infective African trypanosomes: 21 T. brucei brucei, 22 T. b. gambiense, 35 T. b. rhodesiense and 4 T. evansi strains, of which 21 were from Uganda. We constructed a maximum likelihood phylogeny based on 162,210 single nucleotide polymorphisms (SNPs.) The three Trypanosoma brucei sub-species and Trypanosoma evansi are not monophyletic, confirming earlier studies that indicated high similarity among Trypanosoma “sub-species”. We also used discriminant analysis of principal components (DAPC) on the same set of SNPs, identifying seven genetic clusters. These clusters do not correspond well with existing taxonomic classifications, in agreement with the phylogenetic analysis. Geographic origin is reflected in both the phylogeny and clustering analysis. Finally, we used sparse linear discriminant analysis to rank SNPs by their informativeness in differentiating the strains in our data set. As few as 84 SNPs can completely distinguish the strains used in our study, and discriminant analysis was still able to detect genetic structure using as few as 10 SNPs. Our results reinforce earlier results of high genetic similarity between the African Trypanozoon. Despite this, a small subset of SNPs can be used to identify genetic markers that can be used for strain identification or other epidemiological investigations.

Trypanosomes are a major health threat to the people and livestock of Sub-Saharan Africa. Building genomic resources and understanding the genetic structure of these parasites will aid researchers trying to control their spread. To this end, we compared the genomes from 83 trypanosome strains, identifying 162,210 single nucleotide polymorphisms (SNPs) between them. Our analysis shows high genetic similarity between the trypanosomes, and confirms earlier results indicating that the traditional taxonomic classifications do not correspond well with genetic data. Further, we demonstrate that, despite the high genetic similarity, each strain in the study can be distinguished using as few as 84 SNPs, suggesting that a small number of SNPs can be useful for tracking and classifying populations of African trypanosomes.

Epigenetic Regulation of Virulence Gene Expression in Parasitic Protozoa

Protozoan parasites colonize numerous metazoan hosts and insect vectors through their life cycles, with the need to respond quickly and reversibly while encountering diverse and often hostile ecological niches. To succeed, parasites must also persist within individuals until transmission between hosts is achieved. Several parasitic protozoa cause a huge burden of disease in humans and livestock, and here we focus on the parasites that cause malaria and African trypanosomiasis. Efforts to understand how these pathogens adapt to survive in varied host environments, cause disease, and transmit between hosts have revealed a wealth of epigenetic phenomena. Epigenetic switching mechanisms appear to be ideally suited for the regulation of clonal antigenic variation underlying successful parasitism. We review the molecular players and complex mechanistic layers that mediate the epigenetic regulation of virulence gene expression. Understanding epigenetic processes will aid the development of antiparasitic therapeutics.

A co-evolutionary arms race: trypanosomes shaping the human genome, humans shaping the trypanosome genome

Trypanosoma brucei is the causative agent of African sleeping sickness in humans and one of several pathogens that cause the related veterinary disease Nagana. A complex co-evolution has occurred between these parasites and primates that led to the emergence of trypanosome-specific defences and counter-measures. The first line of defence in humans and several other catarrhine primates is the trypanolytic protein apolipoprotein-L1 (APOL1) found within two serum protein complexes, trypanosome lytic factor 1 and 2 (TLF-1 and TLF-2). Two sub-species of T. brucei have evolved specific mechanisms to overcome this innate resistance, Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense. In T. b. rhodesiense, the presence of the serum resistance associated (SRA) gene, a truncated variable surface glycoprotein (VSG), is sufficient to confer resistance to lysis. The resistance mechanism of T. b. gambiense is more complex, involving multiple components: reduction in binding affinity of a receptor for TLF, increased cysteine protease activity and the presence of the truncated VSG, T. b. gambiense-specific glycoprotein (TgsGP). In a striking example of co-evolution, evidence is emerging that primates are responding to challenge by T. b. gambiense and T. b. rhodesiense, with several populations of humans and primates displaying resistance to infection by these two sub-species.

Genetic and structural study of DNA-directed RNA polymerase II of Trypanosoma brucei, towards the designing of novel antiparasitic agents

Trypanosoma brucei brucei (TBB) belongs to the unicellular parasitic protozoa organisms, specifically to the Trypanosoma genus of the Trypanosomatidae class. A variety of different vertebrate species can be infected by TBB, including humans and animals. Under particular conditions, the TBB can be hosted by wild and domestic animals; therefore, an important reservoir of infection always remains available to transmit through tsetse flies. Although the TBB parasite is one of the leading causes of death in the most underdeveloped countries, to date there is neither vaccination available nor any drug against TBB infection. The subunit RPB1 of the TBB DNA-directed RNA polymerase II (DdRpII) constitutes an ideal target for the design of novel inhibitors, since it is instrumental role is vital for the parasite’s survival, proliferation, and transmission. A major goal of the described study is to provide insights for novel anti-TBB agents via a state-of-the-art drug discovery approach of the TBB DdRpII RPB1. In an attempt to understand the function and action mechanisms of this parasite enzyme related to its molecular structure, an in-depth evolutionary study has been conducted in parallel to the in silico molecular designing of the 3D enzyme model, based on state-of-the-art comparative modelling and molecular dynamics techniques. Based on the evolutionary studies results nine new invariant, first-time reported, highly conserved regions have been identified within the DdRpII family enzymes. Consequently, those patches have been examined both at the sequence and structural level and have been evaluated in regard to their pharmacological targeting appropriateness. Finally, the pharmacophore elucidation study enabled us to virtually in silico screen hundreds of compounds and evaluate their interaction capabilities with the enzyme. It was found that a series of chlorine-rich set of compounds were the optimal inhibitors for the TBB DdRpII RPB1 enzyme. All-in-all, herein we present a series of new sites on the TBB DdRpII RPB1 of high pharmacological interest, alongside the construction of the 3D model of the enzyme and the suggestion of a new in silico pharmacophore model for fast screening of potential inhibiting agents.

Leishmania infections: Molecular targets and diagnosis

Progress in the diagnosis of leishmaniases depends on the development of effective methods and the discovery of suitable biomarkers. We propose firstly an update classification of Leishmania species and their synonymies. We demonstrate a global map highlighting the geography of known endemic Leishmania species pathogenic to humans. We summarize a complete list of techniques currently in use and discuss their advantages and limitations. The available data highlights the benefits of molecular markers in terms of their sensitivity and specificity to quantify variation from the subgeneric level to species complexes, (sub) species within complexes, and individual populations and infection foci. Each DNA-based detection method is supplied with a comprehensive description of markers and primers and proposal for a classification based on the role of each target and primer in the detection, identification and quantification of leishmaniasis infection. We outline a genome-wide map of genes informative for diagnosis that have been used for Leishmania genotyping. Furthermore, we propose a classification method based on the suitability of well-studied molecular markers for typing the 21 known Leishmania species pathogenic to humans. This can be applied to newly discovered species and to hybrid strains originating from inter-species crosses. Developing more effective and sensitive diagnostic methods and biomarkers is vital for enhancing Leishmania infection control programs.

Euglena gracilis Genome and Transcriptome: Organelles, Nuclear Genome Assembly Strategies and Initial Features

Euglena gracilis is a major component of the aquatic ecosystem and together with closely related species, is ubiquitous worldwide. Euglenoids are an important group of protists, possessing a secondarily acquired plastid and are relatives to the Kinetoplastidae, which themselves have global impact as disease agents. To understand the biology of E. gracilis, as well as to provide further insight into the evolution and origins of the Kinetoplastidae, we embarked on sequencing the nuclear genome; the plastid and mitochondrial genomes are already in the public domain. Earlier studies suggested an extensive nuclear DNA content, with likely a high degree of repetitive sequence, together with significant extrachromosomal elements. To produce a list of coding sequences we have combined transcriptome data from both published and new sources, as well as embarked on de novo sequencing using a combination of 454, Illumina paired end libraries and long PacBio reads. Preliminary analysis suggests a surprisingly large genome approaching 2 Gbp, with a highly fragmented architecture and extensive repeat composition. Over 80% of the RNAseq reads from E. gracilis maps to the assembled genome sequence, which is comparable with the well assembled genomes of T. brucei and T. cruzi. In order to achieve this level of assembly we employed multiple informatics pipelines, which are discussed here. Finally, as a preliminary view of the genome architecture, we discuss the tubulin and calmodulin genes, which highlight potential novel splicing mechanisms.

Grandeur Alliances: Symbiont Metabolic Integration and Obligate Arthropod Hematophagy

Several arthropod taxa live exclusively on vertebrate blood. This food source lacks essential metabolites required for the maintenance of metabolic homeostasis, and as such, these arthropods have formed symbioses with nutrient-supplementing microbes that facilitate their host's 'hematophagous' feeding ecology. Herein we highlight metabolic contributions of bacterial symbionts that reside within tsetse flies, bed bugs, lice, reduviid bugs, and ticks, with specific emphasis on B vitamin and cofactor biosynthesis. Importantly, these arthropods can transmit pathogens of medical and veterinary relevance and/or cause infestations that induce psychological and dermatological distress. Microbial metabolites, and the biochemical pathways that generate them, can serve as specific targets of novel control mechanisms aimed at disrupting the metabolism of hematophagous arthropods, thus combatting pest invasion and vector-borne pathogen transmission.

Apolipoprotein L1 Variant Associated with Increased Susceptibility to Trypanosome Infection

African trypanosomes, except Trypanosoma brucei gambiense and Trypanosoma brucei rhodesiense, which cause human African trypanosomiasis, are lysed by the human serum protein apolipoprotein L1 (ApoL1). These two subspecies can resist human ApoL1 because they express the serum resistance proteins T. b. gambiense glycoprotein (TgsGP) and serum resistance-associated protein (SRA), respectively. Whereas in T. b. rhodesiense, SRA is necessary and sufficient to inhibit ApoL1, in T. b. gambiense, TgsGP cannot protect against high ApoL1 uptake, so different additional mechanisms contribute to limit this uptake. Here we report a complex interplay between trypanosomes and an ApoL1 variant, revealing important insights into innate human immunity against these parasites. Using whole-genome sequencing, we characterized an atypical T. b. gambiense infection in a patient in Ghana. We show that the infecting trypanosome has diverged from the classical T. b. gambiense strains and lacks the TgsGP defense mechanism against human serum. By sequencing the ApoL1 gene of the patient and subsequent in vitro mutagenesis experiments, we demonstrate that a homozygous missense substitution (N264K) in the membrane-addressing domain of this ApoL1 variant knocks down the trypanolytic activity, allowing the trypanosome to avoid ApoL1-mediated immunity.

Most African trypanosomes are lysed by the ApoL1 protein in human serum. Only the subspecies Trypanosoma b. gambiense and T. b. rhodesiense can resist lysis by ApoL1 because they express specific serum resistance proteins. We here report a complex interplay between trypanosomes and an ApoL1 variant characterized by a homozygous missense substitution (N264K) in the domain that we hypothesize interacts with the endolysosomal membranes of trypanosomes. The N264K substitution knocks down the lytic activity of ApoL1 against T. b. gambiense strains lacking the TgsGP defense mechanism and against T. b. rhodesiense if N264K is accompanied by additional substitutions in the SRA-interacting domain. Our data suggest that populations with high frequencies of the homozygous N264K ApoL1 variant may be at increased risk of contracting human African trypanosomiasis.

DNA Recombination Strategies During Antigenic Variation in the African Trypanosome

Survival of the African trypanosome in its mammalian hosts has led to the evolution of antigenic variation, a process for evasion of adaptive immunity that has independently evolved in many other viral, bacterial and eukaryotic pathogens. The essential features of trypanosome antigenic variation have been understood for many years and comprise a dense, protective Variant Surface Glycoprotein (VSG) coat, which can be changed by recombination-based and transcription-based processes that focus on telomeric VSG gene transcription sites. However, it is only recently that the scale of this process has been truly appreciated. Genome sequencing of Trypanosoma brucei has revealed a massive archive of >1000 VSG genes, the huge majority of which are functionally impaired but are used to generate far greater numbers of VSG coats through segmental gene conversion. This chapter will discuss the implications of such VSG diversity for immune evasion by antigenic variation, and will consider how this expressed diversity can arise, drawing on a growing body of work that has begun to examine the proteins and sequences through which VSG switching is catalyzed. Most studies of trypanosome antigenic variation have focused on T. brucei, the causative agent of human sleeping sickness. Other work has begun to look at antigenic variation in animal-infective trypanosomes, and we will compare the findings that are emerging, as well as consider how antigenic variation relates to the dynamics of host-trypanosome interaction.

De Novo Genome Assembly Shows Genome Wide Similarity between Trypanosoma brucei brucei and Trypanosoma brucei rhodesiense

Background

Trypanosoma brucei is a eukaryotic pathogen which causes African trypanosomiasis. It is notable for its variant surface glycoprotein (VSG) coat, which undergoes antigenic variation enabled by a large suite of VSG pseudogenes, allowing for persistent evasion of host adaptive immunity. While Trypanosoma brucei rhodesiense (Tbr) and T. b gambiense (Tbg) are human infective, related T. b. brucei (Tbb) is cleared by human sera. A single gene, the Serum Resistance Associated (SRA) gene, confers Tbr its human infectivity phenotype. Potential genetic recombination of this gene between Tbr and non-human infective Tbb strains has significant epidemiological consequences for Human African Trypanosomiasis outbreaks.

Results

Using long and short read whole genome sequencing, we generated a hybrid de novo assembly of a Tbr strain, producing 4,210 scaffolds totaling approximately 38.8 megabases, which comprise a significant proportion of the Tbr genome, and thus represents a valuable tool for a comparative genomics analyses among human and non-human infective T. brucei and future complete genome assembly. We detected 5,970 putative genes, of which two, an alcohol oxidoreductase and a pentatricopeptide repeat-containing protein, were members of gene families common to all T. brucei subspecies, but variants specific to the Tbr strain sequenced in this study. Our findings confirmed the extremely high level of genomic similarity between the two parasite subspecies found in other studies.

Conclusions

We confirm at the whole genome level high similarity between the two Tbb and Tbr strains studied. The discovery of extremely minor genomic differentiation between Tbb and Tbr suggests that the transference of the SRA gene via genetic recombination could potentially result in novel human infective strains, thus all genetic backgrounds of T. brucei should be considered potentially human infective in regions where Tbr is prevalent.

Whole genome sequencing shows sleeping sickness relapse is due to parasite regrowth and not reinfection

The trypanosome Trypanosoma brucei gambiense (Tbg) is a cause of human African trypanosomiasis (HAT) endemic to many parts of sub-Saharan Africa. The disease is almost invariably fatal if untreated and there is no vaccine, which makes monitoring and managing drug resistance highly relevant. A recent study of HAT cases from the Democratic Republic of the Congo reported a high incidence of relapses in patients treated with melarsoprol. Of the 19 Tbg strains isolated from patients enrolled in this study, four pairs were obtained from the same patient before treatment and after relapse. We used whole genome sequencing to investigate whether these patients were infected with a new strain, or if the original strain had regrown to pathogenic levels. Clustering analysis of 5938 single nucleotide polymorphisms supports the hypothesis of regrowth of the original strain, as we found that strains isolated before and after treatment from the same patient were more similar to each other than to other isolates. We also identified 23 novel genes that could affect melarsoprol sensitivity, representing a promising new set of targets for future functional studies. This work exemplifies the utility of using evolutionary approaches to provide novel insights and tools for disease control.

Histone H3 Variant Regulates RNA Polymerase II Transcription Termination and Dual Strand Transcription of siRNA Loci in Trypanosoma brucei

Base J, β-D-glucosyl-hydroxymethyluracil, is a chromatin modification of thymine in the nuclear DNA of flagellated protozoa of the order Kinetoplastida. In Trypanosoma brucei, J is enriched, along with histone H3 variant (H3.V), at sites involved in RNA Polymerase (RNAP) II termination and telomeric sites involved in regulating variant surface glycoprotein gene (VSG) transcription by RNAP I. Reduction of J in T. brucei indicated a role of J in the regulation of RNAP II termination, where the loss of J at specific sites within polycistronic gene clusters led to read-through transcription and increased expression of downstream genes. We now demonstrate that the loss of H3.V leads to similar defects in RNAP II termination within gene clusters and increased expression of downstream genes. Gene derepression is intensified upon the subsequent loss of J in the H3.V knockout. mRNA-seq indicates gene derepression includes VSG genes within the silent RNAP I transcribed telomeric gene clusters, suggesting an important role for H3.V in telomeric gene repression and antigenic variation. Furthermore, the loss of H3.V at regions of overlapping transcription at the end of convergent gene clusters leads to increased nascent RNA and siRNA production. Our results suggest base J and H3.V can act independently as well as synergistically to regulate transcription termination and expression of coding and non-coding RNAs in T. brucei, depending on chromatin context (and transcribing polymerase). As such these studies provide the first direct evidence for histone H3.V negatively influencing transcription elongation to promote termination.