The evolution of pathogenic trypanosomes

Mitochondrial DNAs provide insight into trypanosome phylogeny and molecular evolution

Background

Trypanosomes are single-celled eukaryotic parasites characterised by the unique biology of their mitochondrial DNA. African livestock trypanosomes impose a major burden on agriculture across sub-Saharan Africa, but are poorly understood compared to those that cause sleeping sickness and Chagas disease in humans. Here we explore the potential of the maxicircle, a component of trypanosome mitochondrial DNA to study the evolutionary history of trypanosomes.

Results

We used long-read sequencing to completely assemble maxicircle mitochondrial DNA from four previously uncharacterized African trypanosomes, and leveraged these assemblies to scaffold and assemble a further 103 trypanosome maxicircle gene coding regions from published short-read data. While synteny was largely conserved, there were repeated, independent losses of Complex I genes. Comparison of pre-edited and non-edited genes revealed the impact of RNA editing on nucleotide composition, with non-edited genes approaching the limits of GC loss. African tsetse-transmitted trypanosomes showed high levels of RNA editing compared to other trypanosomes. The gene coding regions of maxicircle mitochondrial DNAs were used to construct time-resolved phylogenetic trees, revealing deep divergence events among isolates of the pathogens Trypanosoma brucei and T. congolense.

Conclusions

Our data represents a new resource for experimental and evolutionary analyses of trypanosome phylogeny, molecular evolution and function. Molecular clock analyses yielded a timescale for trypanosome evolution congruent with major biogeographical events in Africa and revealed the recent emergence of Trypanosoma brucei gambiense and T. equiperdum, major human and animal pathogens.

The First Non-LRV RNA Virus in Leishmania

In this work, we describe the first Leishmania-infecting leishbunyavirus-the first virus other than Leishmania RNA virus (LRV) found in trypanosomatid parasites. Its host is Leishmania martiniquensis, a human pathogen causing infections with a wide range of manifestations from asymptomatic to severe visceral disease. This virus (LmarLBV1) possesses many characteristic features of leishbunyaviruses, such as tripartite organization of its RNA genome, with ORFs encoding RNA-dependent RNA polymerase, surface glycoprotein, and nucleoprotein on L, M, and S segments, respectively. Our phylogenetic analyses suggest that LmarLBV1 originated from leishbunyaviruses of monoxenous trypanosomatids and, probably, is a result of genomic re-assortment. The LmarLBV1 facilitates parasites' infectivity in vitro in primary murine macrophages model. The discovery of a virus in L. martiniquensis poses the question of whether it influences pathogenicity of this parasite in vivo, similarly to the LRV in other Leishmania species.

Comparative genomics of Leishmania (Mundinia)

Background

Trypanosomatids of the genus Leishmania are parasites of mammals or reptiles transmitted by bloodsucking dipterans. Many species of these flagellates cause important human diseases with clinical symptoms ranging from skin sores to life-threatening damage of visceral organs. The genus Leishmania contains four subgenera: Leishmania, Sauroleishmania, Viannia, and Mundinia. The last subgenus has been established recently and remains understudied, although Mundinia contains human-infecting species. In addition, it is interesting from the evolutionary viewpoint, representing the earliest branch within the genus and possibly with a different type of vector. Here we analyzed the genomes of L. (M.) martiniquensis, L. (M.) enriettii and L. (M.) macropodum to better understand the biology and evolution of these parasites.

Results

All three genomes analyzed were approximately of the same size (~ 30 Mb) and similar to that of L. (Sauroleishmania) tarentolae, but smaller than those of the members of subgenera Leishmania and Viannia, or the genus Endotrypanum (~ 32 Mb). This difference was explained by domination of gene losses over gains and contractions over expansions at the Mundinia node, although only a few of these genes could be identified. The analysis predicts significant changes in the Mundinia cell surface architecture, with the most important ones relating to losses of LPG-modifying side chain galactosyltransferases and arabinosyltransferases, as well as β-amastins. Among other important changes were gene family contractions for the oxygen-sensing adenylate cyclases and FYVE zinc finger-containing proteins.

Conclusions

We suggest that adaptation of Mundinia to different vectors and hosts has led to alternative host-parasite relationships and, thereby, made some proteins redundant. Thus, the evolution of genomes in the genus Leishmania and, in particular, in the subgenus Mundinia was mainly shaped by host (or vector) switches.

Development of Phytomonas lipae sp. n. (Kinetoplastea: Trypanosomatidae) in the true bug Coreus marginatus (Heteroptera: Coreidae) and insights into the evolution of life cycles in the genus Phytomonas

Here we described a new trypanosomatid species, Phytomonas lipae, parasitizing the dock bug Coreus marginatus based on axenic culture and in vivo material. Using light and electron microscopy we characterized the development of this flagellate in the intestine, hemolymph and salivary glands of its insect host. The intestinal promastigotes of Phytomonas lipae do not divide and occur only in the anterior part of the midgut. From there they pass into hemolymph, increasing in size, and then to salivary glands, where they actively proliferate without attachment to the host's epithelium and form infective endomastigotes. We conducted molecular phylogenetic analyses based on 18s rRNA, gGAPDH and HSP83 gene sequences, of which the third marker performed the best in terms of resolving phylogenetic relationships within the genus Phytomonas. Our inference demonstrated rather early origin of the lineage comprising the new species, right after that of P. oxycareni, which represents the earliest known branch within the Phytomonas clade. This allowed us to compare the development of P. lipae and three other Phytomonas spp. in their insect hosts and reconstruct the vectorial part of the life cycle of their common ancestor.

Phylogenetic Studies

Phylogenetics is an important component of the systems biology approach. Knowledge about evolution of the genus Leishmania is essential to understand various aspects of basic biology of these parasites, such as parasite-host or parasite-vector relationships, biogeography, or epidemiology. Here, we present a comprehensive guideline for performing phylogenetic studies based on DNA sequence data, but with principles that can be adapted to protein sequences or other molecular markers. It is presented as a compilation of the most commonly used genetic targets for phylogenetic studies of Leishmania, including their respective primers for amplification and references, as well as details of PCR assays. Guidelines are, then, presented to choose the best targets in relation to the types of samples under study. Finally, and importantly, instructions are given to obtain optimal sequences, alignments, and datasets for the subsequent data analysis and phylogenetic inference. Different bioinformatics methods and software for phylogenetic inference are presented and explained. This chapter aims to provide a compilation of methods and generic guidelines to conduct phylogenetics of Leishmania for nonspecialists.

Functional genomics of trypanosomatids

The decoding of the Tritryp reference genomes nearly 7 years ago provided a first peek into the biology of pathogenic trypanosomatids and a blueprint that has paved the way for genome-wide studies. Although 60-70% of the predicted protein coding genes in Trypanosoma brucei, Trypanosoma cruzi and Leishmania major remain unannotated, the functional genomics landscape is rapidly changing. Facilitated by the advent of next-generation sequencing technologies, improved structural and functional annotation and genes and their products are emerging. Information is also growing for the interactions between cellular components as transcriptomes, regulatory networks and metabolomes are characterized, ushering in a new era of systems biology. Simultaneously, the launch of comparative sequencing of multiple strains of kinetoplastids will finally lead to the investigation of a vast, yet to be explored, evolutionary and pathogenomic space.

Pathogenesis of chagas' disease: parasite persistence and autoimmunity

Acute Trypanosoma cruzi infections can be asymptomatic, but chronically infected individuals can die of Chagas' disease. The transfer of the parasite mitochondrial kinetoplast DNA (kDNA) minicircle to the genome of chagasic patients can explain the pathogenesis of the disease; in cases of Chagas' disease with evident cardiomyopathy, the kDNA minicircles integrate mainly into retrotransposons at several chromosomes, but the minicircles are also detected in coding regions of genes that regulate cell growth, differentiation, and immune responses. An accurate evaluation of the role played by the genotype alterations in the autoimmune rejection of self-tissues in Chagas' disease is achieved with the cross-kingdom chicken model system, which is refractory to T. cruzi infections. The inoculation of T. cruzi into embryonated eggs prior to incubation generates parasite-free chicks, which retain the kDNA minicircle sequence mainly in the macrochromosome coding genes. Crossbreeding transfers the kDNA mutations to the chicken progeny. The kDNA-mutated chickens develop severe cardiomyopathy in adult life and die of heart failure. The phenotyping of the lesions revealed that cytotoxic CD45, CD8(+) γδ, and CD8α(+) T lymphocytes carry out the rejection of the chicken heart. These results suggest that the inflammatory cardiomyopathy of Chagas' disease is a genetically driven autoimmune disease.

CRK9 contributes to regulation of mitosis and cytokinesis in the procyclic form of Trypanosoma brucei

Background

The Trypanosoma brucei cell cycle is regulated by combinations of cyclin/CRKs (cdc2 related kinases). Recently, two additional cyclins (CYC10, CYC11) and six new CRK (CRK7-12) homologues were identified in the T. brucei genome database [1,2].

Results

Individual RNAi knockdowns of these new proteins in the procyclic form of T. brucei showed no apparent phenotype except for the CRK9 depletion, which enriched the cells in G2/M phase. But a similar CRK9 knockdown in the bloodstream form caused no apparent phenotype. CRK9 lacks the typical PSTAIRE motif for cyclin binding and the phenylalanine "gatekeeper" but binds to cyclin B2 in vitro and localizes to the nucleus in both forms of T. brucei. CRK9-depleted procyclic-form generated no detectable anucleate cells, suggesting an inhibition of cytokinesis by CRK9 depletion as well. The knockdown enriched cells with one nucleus, one kinetoplast and two closely associated basal bodies with an average distance of 1.08 mm in between, which was shorter than the control value of 1.36 μm, and the cells became morphologically deformed and rounded with time.

Conclusion

CRK9 may play a role in mediating the segregation between the two kinetoplast/basal body pairs prior to cytokinetic initiation. Since such a segregation over a relatively significant distance is essential for cytokinetic initiation only in the procyclic but may not be in the bloodstream form, CRK9 could be specifically involved in regulating cytokinetic initiation in the procyclic form of T. brucei.

Evidence of recent interkingdom horizontal gene transfer between bacteria and Candida parapsilosis

BACKGROUND

To date very few incidences of interdomain gene transfer into fungi have been identified. Here, we used the emerging genome sequences of Candida albicans WO-1, Candida tropicalis, Candida parapsilosis, Clavispora lusitaniae, Pichia guilliermondii, and Lodderomyces elongisporus to identify recent interdomain HGT events. We refer to these as CTG species because they translate the CTG codon as serine rather than leucine, and share a recent common ancestor.

RESULTS

Phylogenetic and syntenic information infer that two C. parapsilosis genes originate from bacterial sources. One encodes a putative proline racemase (PR). Phylogenetic analysis also infers that there were independent transfers of bacterial PR enzymes into members of the Pezizomycotina, and protists. The second HGT gene in C. parapsilosis belongs to the phenazine F (PhzF) superfamily. Most CTG species also contain a fungal PhzF homolog. Our phylogeny suggests that the CTG homolog originated from an ancient HGT event, from a member of the proteobacteria. An analysis of synteny suggests that C. parapsilosis has lost the endogenous fungal form of PhzF, and subsequently reacquired it from a proteobacterial source. There is evidence that Schizosaccharomyces pombe and Basidiomycotina also obtained a PhzF homolog through HGT.

CONCLUSION

Our search revealed two instances of well-supported HGT from bacteria into the CTG clade, both specific to C. parapsilosis. Therefore, while recent interkingdom gene transfer has taken place in the CTG lineage, its occurrence is rare. However, our analysis will not detect ancient gene transfers, and we may have underestimated the global extent of HGT into CTG species.

Variables that modulate the spatial distribution of Trypanosoma cruzi and Trypanosoma evansi in the Brazilian Pantanal

An evaluation was made on how the landscape and cattle ranching affect the transmission cycles and the patterns of tripanosomatid infection (Trypanosoma cruzi and Trypanosoma evansi) of small wild mammals in the Pantanal. This region comprises a large natural environment with a multiplicity of habitats, wide variety of biodiversity besides the presence of livestock. T. cruzi and T. evansi infections were evaluated by parasitological and serological methods in one preserved and one cattle ranching area. The diversity of the small mammal fauna showed to be the same in the two studied areas, however, their relative abundance was different. Distinct enzootiological scenarios of both trypanosomatids could be observed. Transmission of T. cruzi occurred mainly in forested areas, in the two study areas, while T. evansi occurred dispersed among all habitats studied in the unpreserved area. The arboreal rodent Oecomys mamorae, the most abundant species in both areas, displayed high T. cruzi and T. evansi serum prevalence and parasitemias. Also, the caviomorph rodent Thrichomys pachyurus was shown to be an important host due to its expressive relative abundance, prevalence of infection by both trypanosomatid species and a broad range use of habitats. The role of the small mammal fauna in the transmission cycle of both trypanosomes species seems to be distinct according to land use since we found a broad range of T. evansi infected hosts in the preserved area in contrast to cattle ranching area and a half number of the rodents species infected with T. cruzi in unpreserved in comparison to protect area. The present study showed that cattle ranching in this study area did not enhance overall prevalence of T. cruzi infection among small wild mammals. Together with the observation that small mammals diversity in FA is similar to RN area suggest that ranching activity may also not necessarily conduct to biodiversity loss or risk of Chagas disease.

Comparative genomics of trypanosomatid parasitic protozoa

A comparison of gene content and genome architecture of Trypanosoma brucei, Trypanosoma cruzi, and Leishmania major, three related pathogens with different life cycles and disease pathology, revealed a conserved core proteome of about 6200 genes in large syntenic polycistronic gene clusters. Many species-specific genes, especially large surface antigen families, occur at nonsyntenic chromosome-internal and subtelomeric regions. Retroelements, structural RNAs, and gene family expansion are often associated with syntenic discontinuities that-along with gene divergence, acquisition and loss, and rearrangement within the syntenic regions-have shaped the genomes of each parasite. Contrary to recent reports, our analyses reveal no evidence that these species are descended from an ancestor that contained a photosynthetic endosymbiont.

Evaluation of the Presence of a Thapsigargin-Sensitive Calcium Store in Trypanosomatids Using Trypanosoma evansi as a Model

Ca2+ plays an important role in the regulation of several important activities in different trypanosomatids. These parasites possess a Ca2+ transport system in the endoplasmic reticulum (ER) involved in Ca2+ homeostasis, which has been reported to be insensitive to thapsigargin, a classical inhibitor of the sarcoplasmic-ER Ca2+ adenosine triphosphatase (ATPase) (SERCA) in most eukaryotic cells. However, currently there is a controversy regarding the existence of a thapsigargin-sensitive ER Ca2+ store in these parasites. Therefore, we decided to explore the effect of this inhibitor using different methodological approaches. First, we selected Trypanosoma evansi as a parasite model to warrant the homogeneity of the population because this parasite has only a single life cycle, i.e., bloodstream-form trypomastigotes. Second, we compared the thapsigargin effect on Ca2+ homeostasis by spectrophotometrical Ca2+ measurements using 3 different approaches: whole-cell populations, cells that have been permeabilized by treatment with digitonin, and intact single cells. Our results demonstrate that a low concentration of thapsigargin induces Ca2+ release from intracellular Ca2+ stores in this parasite, which can be observed independently of the method used. Furthermore, the addition of thapsigargin before or after nigericin did not abolish its effect, showing that thapsigargin acts specifically on the ER. In conclusion, our results indicate the presence of a nonmitochondrial thapsigargin-sensitive Ca2+ store in T. evansi.

Gene synteny and evolution of genome architecture in trypanosomatids

The trypanosomatid protozoa Trypanosoma brucei, Trypanosoma cruzi and Leishmania major are related human pathogens that cause markedly distinct diseases. Using information from genome sequencing projects currently underway, we have compared the sequences of large chromosomal fragments from each species. Despite high levels of divergence at the sequence level, these three species exhibit a striking conservation of gene order, suggesting that selection has maintained gene order among the trypanosomatids over hundreds of millions of years of evolution. The few sites of genome rearrangement between these species are marked by the presence of retrotransposon-like elements, suggesting that retrotransposons may have played an important role in shaping trypanosomatid genome organization. A degenerate retroelement was identified in L. major by examining the regions near breakage points of the synteny. This is the first such element found in L. major suggesting that retroelements were found in the common ancestor of all three species.

Parasitism, the diversity of life, and paleoparasitology

The parasite-host-environment system is dynamic, with several points of equilibrium. This makes it difficult to trace the thresholds between benefit and damage, and therefore, the definitions of commensalism, mutualism, and symbiosis become worthless. Therefore, the same concept of parasitism may encompass commensalism, mutualism, and symbiosis. Parasitism is essential for life. Life emerged as a consequence of parasitism at the molecular level, and intracellular parasitism created evolutive events that allowed species to diversify. An ecological and evolutive approach to the study of parasitism is presented here. Studies of the origin and evolution of parasitism have new perspectives with the development of molecular paleoparasitology, by which ancient parasite and host genomes can be recovered from disappeared populations. Molecular paleoparasitology points to host-parasite co-evolutive mechanisms of evolution traceable through genome retrospective studies.

Some current problems in the systematics of Trypanosomatids

Trypanosomatids have been traditionally allocated to a number of genera that were described based on morphological features and host range. Recently molecular studies have provided new data that has allowed a reexamination of the genera. While in some cases the molecular data has been in agreement with the morphological characters they have also reinforced existing doubts about some current generic divisions as well as raising new concerns. A revision of the trypanosomatid genera is required. Suggested features of such a revision would include: (1) The possible division of Trypanosoma into new genera to reflect the wide genetic diversity of this group; (2) The inclusion of Leishmania, Sauroleishmania and Endotrypanum within a single genus given their high genetic affinity; (3) The complete revision of the monogenetic typanosomatid genera to reflect monophyletic groups; (4) A more precise redescription of Phytomonas so as to only include the monophyletic plant flagellates.