Cell biology of the trypanosome genome. Microbiol Mol Biol Rev. 2010;74:552-569. [PMID: 21119017 DOI: 10.1128/mmbr.00024-10]

Probing the binding and antiparasitic efficacy of azobenzene G-quadruplex ligands to investigate G4 ligand design

Novel strategies against parasitic infections are of great importance. Here, we describe a G4 DNA ligand with subnanomolar antiparasitic activity against T. brucei and a remarkable selectivity index (IC50 MRC-5/T. brucei) of 2285-fold. We also correlate the impact of small structural changes to G4 binding activity and antiparasitic activity.

DNA lesions that block transcription induce the death of Trypanosoma cruzi via ATR activation, which is dependent on the presence of R-loops

Trypanosoma cruzi is the etiological agent of Chagas disease and a peculiar eukaryote with unique biological characteristics. DNA damage can block RNA polymerase, activating transcription-coupled nucleotide excision repair (TC-NER), a DNA repair pathway specialized in lesions that compromise transcription. If transcriptional stress is unresolved, arrested RNA polymerase can activate programmed cell death. Nonetheless, how this parasite modulates these processes is unknown. Here, we demonstrate that T. cruzi cell death after UV irradiation, a genotoxic agent that generates lesions resolved by TC-NER, depends on active transcription and is signaled mainly by an apoptotic-like pathway. Pre-treated parasites with α-amanitin, a selective RNA polymerase II inhibitor, become resistant to such cell death. Similarly, the gamma pre-irradiated cells are more resistant to UV when the transcription processes are absent. The Cockayne Syndrome B protein (CSB) recognizes blocked RNA polymerase and can initiate TC-NER. Curiously, CSB overexpression increases parasites' cell death shortly after UV exposure. On the other hand, at the same time after irradiation, the single-knockout CSB cells show resistance to the same treatment. UV-induced fast death is signalized by the exposition of phosphatidylserine to the outer layer of the membrane, indicating a cell death mainly by an apoptotic-like pathway. Furthermore, such death is suppressed in WT parasites pre-treated with inhibitors of ataxia telangiectasia and Rad3-related (ATR), a key DDR kinase. Signaling for UV radiation death may be related to R-loops since the overexpression of genes associated with the resolution of these structures suppress it. Together, results suggest that transcription blockage triggered by UV radiation activates an ATR-dependent apoptosis-like mechanism in T. cruzi, with the participation of CSB protein in this process.

Exploring the functionality and conservation of Alba proteins in Trypanosoma cruzi: A focus on biological diversity and RNA binding ability

Trypanosoma cruzi is the causative agent of Chagas disease, as well as a trypanosomatid parasite with a complex biological cycle that requires precise mechanisms for regulating gene expression. In Trypanosomatidae, gene regulation occurs mainly at the mRNA level through the recognition of cis elements by RNA-binding proteins (RBPs). Alba family members are ubiquitous DNA/RNA-binding proteins with representatives in trypanosomatid parasites functionally related to gene expression regulation. Although T. cruzi possesses two groups of Alba proteins (Alba1/2 and Alba30/40), their functional role remains poorly understood. Thus, herein, a characterization of T. cruzi Alba (TcAlba) proteins was undertaken. Physicochemical, structural, and phylogenetic analysis of TcAlba showed features compatible with RBPs, such as hydrophilicity, RBP domains/motifs, and evolutionary conservation of the Alba-domain, mainly regarding other trypanosomatid Alba. However, in silico RNA interaction analysis of T. cruzi Alba proteins showed that TcAlba30/40 proteins, but not TcAlba1/2, would directly interact with the assayed RNA molecules, suggesting that these two groups of TcAlba proteins have different targets. Given the marked differences existing between both T. cruzi Alba groups (TcAlba1/2 and TcAlba30/40), regarding sequence divergence, RNA binding potential, and life-cycle expression patterns, we suggest that they would be involved in different biological processes.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

Expectations of duplicate gene retention under the gene duplicability hypothesis

BACKGROUND

Gene duplication is an important process in evolution. What causes some genes to be retained after duplication and others to be lost is a process not well understood. The most prevalent theory is the gene duplicability hypothesis, that something about the function and number of interacting partners (number of subunits of protein complex, etc.), determines whether copies have more opportunity to be retained for long evolutionary periods. Some genes are also more susceptible to dosage balance effects following WGD events, making them more likely to be retained for longer periods of time. One would expect these processes that affect the retention of duplicate copies to affect the conditional probability ratio after consecutive whole genome duplication events. The probability that a gene will be retained after a second whole genome duplication event (WGD2), given that it was retained after the first whole genome duplication event (WGD1) versus the probability a gene will be retained after WGD2, given it was lost after WGD1 defines the probability ratio that is calculated.

RESULTS

Since duplicate gene retention is a time heterogeneous process, the time between the events (t1) and the time since the most recent event (t2) are relevant factors in calculating the expectation for observation in any genome. Here, we use a survival analysis framework to predict the probability ratio for genomes with different values of t1 and t2 under the gene duplicability hypothesis, that some genes are more susceptible to selectable functional shifts, some more susceptible to dosage compensation, and others only drifting. We also predict the probability ratio with different values of t1 and t2 under the mutational opportunity hypothesis, that probability of retention for certain genes changes in subsequent events depending upon how they were previously retained. These models are nested such that the mutational opportunity model encompasses the gene duplicability model with shifting duplicability over time. Here we present a formalization of the gene duplicability and mutational opportunity hypotheses to characterize evolutionary dynamics and explanatory power in a recently developed statistical framework.

CONCLUSIONS

This work presents expectations of the gene duplicability and mutational opportunity hypotheses over time under different sets of assumptions. This expectation will enable formal testing of processes leading to duplicate gene retention.

Immunoprecipitation of RNA-DNA hybrid interacting proteins in Trypanosoma brucei reveals conserved and novel activities, including in the control of surface antigen expression needed for immune evasion by antigenic variation

RNA-DNA hybrids are epigenetic features of genomes that provide a diverse and growing range of activities. Understanding of these functions has been informed by characterising the proteins that interact with the hybrids, but all such analyses have so far focused on mammals, meaning it is unclear if a similar spectrum of RNA-DNA hybrid interactors is found in other eukaryotes. The African trypanosome is a single-cell eukaryotic parasite of the Discoba grouping and displays substantial divergence in several aspects of core biology from its mammalian host. Here, we show that DNA-RNA hybrid immunoprecipitation coupled with mass spectrometry recovers 602 putative interactors in T. brucei mammal- and insect-infective cells, some providing activities also found in mammals and some lineage-specific. We demonstrate that loss of three factors, two putative helicases and a RAD51 paralogue, alters T. brucei nuclear RNA-DNA hybrid and DNA damage levels. Moreover, loss of each factor affects the operation of the parasite immune survival mechanism of antigenic variation. Thus, our work reveals the broad range of activities contributed by RNA-DNA hybrids to T. brucei biology, including new functions in host immune evasion as well as activities likely fundamental to eukaryotic genome function.

Targeting trypanosomes: how chemogenomics and artificial intelligence can guide drug discovery

Trypanosomatids are protozoan parasites that cause human and animal neglected diseases. Despite global efforts, effective treatments are still much needed. Phenotypic screens have provided several chemical leads for drug discovery, but the mechanism of action for many of these chemicals is currently unknown. Recently, chemogenomic screens assessing the susceptibility or resistance of parasites carrying genome-wide modifications started to define the mechanism of action of drugs at large scale. In this review, we discuss how genomics is being used for drug discovery in trypanosomatids, how integration of chemical and genomics data from these and other organisms has guided prioritisations of candidate therapeutic targets and additional chemical starting points, and how these data can fuel the expansion of drug discovery pipelines into the era of artificial intelligence.

Nucleolar targeting in an early-branching eukaryote suggests a general mechanism for ribosome protein sorting

The compartmentalised eukaryotic cell demands accurate targeting of proteins to the organelles in which they function, whether membrane-bound (like the nucleus) or non-membrane-bound (like the nucleolus). Nucleolar targeting relies on positively charged localisation signals and has received rejuvenated interest since the widespread recognition of liquid-liquid phase separation (LLPS) as a mechanism contributing to nucleolus formation. Here, we exploit a new genome-wide analysis of protein localisation in the early-branching eukaryote Trypanosoma brucei to analyse general nucleolar protein properties. T. brucei nucleolar proteins have similar properties to those in common model eukaryotes, specifically basic amino acids. Using protein truncations and addition of candidate targeting sequences to proteins, we show both homopolymer runs and distributed basic amino acids give nucleolar partition, further aided by a nuclear localisation signal (NLS). These findings are consistent with phase separation models of nucleolar formation and physical protein properties being a major contributing mechanism for eukaryotic nucleolar targeting, conserved from the last eukaryotic common ancestor. Importantly, cytoplasmic ribosome proteins, unlike mitochondrial ribosome proteins, have more basic residues - pointing to adaptation of physicochemical properties to assist segregation.

Trypanosoma brucei histones are heavily modified with combinatorial post-translational modifications and mark Pol II transcription start regions with hyperacetylated H2A

Trypanosomes diverged from the main eukaryotic lineage about 600 million years ago, and display some unusual genomic and epigenetic properties that provide valuable insight into the early processes employed by eukaryotic ancestors to regulate chromatin-mediated functions. We analysed Trypanosoma brucei core histones by high mass accuracy middle-down mass spectrometry to map core histone post-translational modifications (PTMs) and elucidate cis-histone combinatorial PTMs (cPTMs). T. brucei histones are heavily modified and display intricate cPTMs patterns, with numerous hypermodified cPTMs that could contribute to the formation of non-repressive euchromatic states. The Trypanosoma brucei H2A C-terminal tail is hyperacetylated, containing up to five acetylated lysine residues. MNase-ChIP-seq revealed a striking enrichment of hyperacetylated H2A at Pol II transcription start regions, and showed that H2A histones that are hyperacetylated in different combinations localised to different genomic regions, suggesting distinct epigenetic functions. Our genomics and proteomics data provide insight into the complex epigenetic mechanisms used by this parasite to regulate a genome that lacks the transcriptional control mechanisms found in later-branched eukaryotes. The findings further demonstrate the complexity of epigenetic mechanisms that were probably shared with the last eukaryotic common ancestor.

Stage-specific transcription activator ESB1 regulates monoallelic antigen expression in Trypanosoma brucei

Variant surface glycoprotein (VSG) coats bloodstream form Trypanosoma brucei parasites, and monoallelic VSG expression underpins the antigenic variation necessary for pathogenicity. One of thousands of VSG genes is transcribed by RNA polymerase I in a singular nuclear structure called the expression site body (ESB), but how monoallelic VSG transcription is achieved remains unclear. Using a localization screen of 153 proteins we found one, ESB-specific protein 1 (ESB1), that localized only to the ESB and is expressed only in VSG-expressing life cycle stages. ESB1 associates with DNA near the active VSG promoter and is necessary for VSG expression, with overexpression activating inactive VSG promoters. Mechanistically, ESB1 is necessary for recruitment of a subset of ESB components, including RNA polymerase I, revealing that the ESB has separately assembled subdomains. Because many trypanosomatid parasites have divergent ESB1 orthologues yet do not undergo antigenic variation, ESB1 probably represents an important class of transcription regulators.

The genomic basis of host and vector specificity in non-pathogenic trypanosomatids

Trypanosoma theileri, a non-pathogenic parasite of bovines, has a predicted surface protein architecture that likely aids survival in its mammalian host. Their surface proteins are encoded by genes which account for ∼10% of their genome. A non-pathogenic parasite of sheep, Trypanosoma melophagium, is transmitted by the sheep ked and is closely related to T. theileri. To explore host and vector specificity between these species, we sequenced the T. melophagium genome and transcriptome and an annotated draft genome was assembled. T. melophagium was compared to 43 kinetoplastid genomes, including T. theileri. T. melophagium and T. theileri have an AT biased genome, the greatest bias of publicly available trypanosomatids. This trend may result from selection acting to decrease the genomic nucleotide cost. The T. melophagium genome is 6.3Mb smaller than T. theileri and large families of proteins, characteristic of the predicted surface of T. theileri, were found to be absent or greatly reduced in T. melophagium. Instead, T. melophagium has modestly expanded protein families associated with the avoidance of complement-mediated lysis. We propose that the contrasting genomic features of these species is linked to their mode of transmission from their insect vector to their mammalian host. This article has an associated First Person interview with the first author of the paper.

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.

Identification of sequence-specific promoters driving polycistronic transcription initiation by RNA polymerase II in trypanosomes

Protein-coding genes in trypanosomes occur in polycistronic transcription units (PTUs). How RNA polymerase II (Pol II) initiates transcription of PTUs has not been resolved; the current model favors chromatin modifications inducing transcription rather than sequence-specific promoters. Here, we uncover core promoters by functional characterization of Pol II peaks identified by chromatin immunoprecipitation sequencing (ChIP-seq). Two distinct promoters are located between divergent PTUs, each driving unidirectional transcription. Detailed analysis identifies a 75-bp promoter that is necessary and sufficient to drive full reporter expression and contains functional motifs. Analysis of further promoters suggests transcription initiation is regulated and promoters are either focused or dispersed. In contrast to the previous model of unregulated and promoter-independent transcription initiation, we find that sequence-specific promoters determine the initiation of Pol II transcription of protein-coding genes PTUs. These findings in Trypanosoma brucei suggest that in addition of chromatin modifications, promoter motifs-based regulation of gene expression is deeply conserved among eukaryotes.

The SPARC complex defines RNAPII promoters in Trypanosoma brucei

Kinetoplastids are a highly divergent lineage of eukaryotes with unusual mechanisms for regulating gene expression. We previously surveyed 65 putative chromatin factors in the kinetoplastid Trypanosoma brucei. Our analyses revealed that the predicted histone methyltransferase SET27 and the Chromodomain protein CRD1 are tightly concentrated at RNAPII transcription start regions (TSRs). Here, we report that SET27 and CRD1, together with four previously uncharacterized constituents, form the SET27 promoter-associated regulatory complex (SPARC), which is specifically enriched at TSRs. SET27 loss leads to aberrant RNAPII recruitment to promoter sites, accumulation of polyadenylated transcripts upstream of normal transcription start sites, and conversion of some normally unidirectional promoters to bidirectional promoters. Transcriptome analysis in the absence of SET27 revealed upregulated mRNA expression in the vicinity of SPARC peaks within the main body of chromosomes in addition to derepression of genes encoding variant surface glycoproteins (VSGs) located in subtelomeric regions. These analyses uncover a novel chromatin-associated complex required to establish accurate promoter position and directionality.

A resource for improved predictions of Trypanosoma and Leishmania protein three-dimensional structure

AlphaFold2 and RoseTTAfold represent a transformative advance for predicting protein structure. They are able to make very high-quality predictions given a high-quality alignment of the protein sequence with related proteins. These predictions are now readily available via the AlphaFold database of predicted structures and AlphaFold or RoseTTAfold Colaboratory notebooks for custom predictions. However, predictions for some species tend to be lower confidence than model organisms. Problematic species include Trypanosoma cruzi and Leishmania infantum: important unicellular eukaryotic human parasites in an early-branching eukaryotic lineage. The cause appears to be due to poor sampling of this branch of life (Discoba) in the protein sequences databases used for the AlphaFold database and ColabFold. Here, by comprehensively gathering openly available protein sequence data for Discoba species, significant improvements to AlphaFold2 protein structure prediction over the AlphaFold database and ColabFold are demonstrated. This is made available as an easy-to-use tool for the parasitology community in the form of Colaboratory notebooks for generating multiple sequence alignments and AlphaFold2 predictions of protein structure for Trypanosoma, Leishmania and related species.

A nuclear enterprise: zooming in on nuclear organization and gene expression control in the African trypanosome

African trypanosomes are early divergent protozoan parasites responsible for high mortality and morbidity as well as a great economic burden among the world's poorest populations. Trypanosomes undergo antigenic variation in their mammalian hosts, a highly sophisticated immune evasion mechanism. Their nuclear organization and mechanisms for gene expression control present several conventional features but also a number of striking differences to the mammalian counterparts. Some of these unorthodox characteristics, such as lack of controlled transcription initiation or enhancer sequences, render their monogenic antigen transcription, which is critical for successful antigenic variation, even more enigmatic. Recent technological developments have advanced our understanding of nuclear organization and gene expression control in trypanosomes, opening novel research avenues. This review is focused on Trypanosoma brucei nuclear organization and how it impacts gene expression, with an emphasis on antigen expression. It highlights several dedicated sub-nuclear bodies that compartmentalize specific functions, whilst outlining similarities and differences to more complex eukaryotes. Notably, understanding the mechanisms underpinning antigen as well as general gene expression control is of great importance, as it might help designing effective control strategies against these organisms.

TBP and SNAP50 transcription factors bind specifically to the Pr77 promoter sequence from trypanosomatid non-LTR retrotransposons

Background

Trypanosomatid genomes are colonized by active and inactive mobile DNA elements, such as LINE, SINE-like, SIDER and DIRE retrotransposons. These elements all share a 77-nucleotide-long sequence at their 5′ ends, known as Pr77, which activates transcription, thereby generating abundant unspliced and translatable transcripts. However, transcription factors that mediates this process have still not been reported.

Methods

TATA-binding protein (TBP) and small nuclear RNA-activating protein 50 kDa (SNAP50) recombinant proteins and specific antibodies raised against them were generated. Protein capture assay, electrophoretic mobility-shift assays (EMSA) and EMSA competition assays carried out using these proteins and nuclear proteins of the parasite together to specific DNA sequences used as probes allowed detecting direct interaction of these transcription factors to Pr77 sequence.

Results

This study identified TBP and SNAP50 as part of the DNA-protein complex formed by the Pr77 promoter sequence and nuclear proteins of Trypanosoma cruzi. TBP establishes direct and specific contact with the Pr77 sequence, where the DPE and DPE downstream regions are docking sites with preferential binding. TBP binds cooperatively (Hill coefficient = 1.67) to Pr77 and to both strands of the Pr77 sequence, while the conformation of this highly structured sequence is not involved in TBP binding. Direct binding of SNAP50 to the Pr77 sequence is weak and may be mediated by protein–protein interactions through other trypanosomatid nuclear proteins.

Conclusions

Identification of the transcription factors that mediate Pr77 transcription may help to elucidate how these retrotransposons are mobilized within the trypanosomatid genomes and their roles in gene regulation processes in this human parasite.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04803-5.

Nuclear tRNA export in trypanosomes: a journey full of twists and turns guided by tRNA modifications

Transfer RNAs play a key role in protein synthesis. Following transcription, tRNAs are extensively processed prior to their departure from the nucleus to become fully functional during translation. This includes removal of 5′ leaders and 3′ trailers by a specific endo- and/or exonuclease, 3′ CCA tail addition, posttranscriptional modifications and in some cases intron removal. In this minireview, the critical factors of nuclear tRNA trafficking are described based on studies in classical models such as yeast and human cell lines. In addition, recent findings and identification of novel regulatory loops of nuclear tRNA trafficking in trypanosomes are discussed with emphasis on tRNA modifications. The comparison between the representatives of opisthokonts and excavates serves here to understand the evolutionary conservation and diversity of nuclear tRNA export mechanisms.

Testing the CRISPR-Cas9 and glmS ribozyme systems in Leishmania tarentolae

Leishmania parasites include important pathogens and model organisms and are even used for the production of recombinant proteins. However, functional genomics and the characterization of essential genes are often limited in Leishmania because of low-throughput technologies for gene disruption or tagging and the absence of components for RNA interference. Here, we tested the T7 RNA polymerase-dependent CRISPR-Cas9 system by Beneke et al. and the glmS ribozyme-based knock-down system in the model parasite Leishmania tarentolae. We successfully deleted two reference genes encoding the flagellar motility factor Pf16 and the salvage-pathway enzyme adenine phosphoribosyltransferase, resulting in immotile and drug-resistant parasites, respectively. In contrast, we were unable to disrupt the gene encoding the mitochondrial flavoprotein Erv. Cultivation of L. tarentolae in standard BHI medium resulted in a constitutive down-regulation of an episomal mCherry-glmS reporter by 40 to 60%. For inducible knock-downs, we evaluated the growth of L. tarentolae in alternative media and identified supplemented MEM, IMDM and McCoy's 5A medium as candidates. Cultivation in supplemented MEM allowed an inducible, glucosamine concentration-dependent down-regulation of the episomal mCherry-glmS reporter by more than 70%. However, chromosomal glmS-tagging of the genes encoding Pf16, adenine phosphoribosyltransferase or Erv did not reveal a knock-down phenotype. Our data demonstrate the suitability of the CRISPR-Cas9 system for the disruption and tagging of genes in L. tarentolae as well as the limitations of the glmS system, which was restricted to moderate efficiencies for episomal knock-downs and caused no detectable phenotype for chromosomal knock-downs.

Trypanosomatid selenophosphate synthetase structure, function and interaction with selenocysteine lyase

Eukaryotes from the Excavata superphylum have been used as models to study the evolution of cellular molecular processes. Strikingly, human parasites of the Trypanosomatidae family (T. brucei, T. cruzi and L. major) conserve the complex machinery responsible for selenocysteine biosynthesis and incorporation in selenoproteins (SELENOK/SelK, SELENOT/SelT and SELENOTryp/SelTryp), although these proteins do not seem to be essential for parasite viability under laboratory controlled conditions. Selenophosphate synthetase (SEPHS/SPS) plays an indispensable role in selenium metabolism, being responsible for catalyzing the formation of selenophosphate, the biological selenium donor for selenocysteine synthesis. We solved the crystal structure of the L. major selenophosphate synthetase and confirmed that its dimeric organization is functionally important throughout the domains of life. We also demonstrated its interaction with selenocysteine lyase (SCLY) and showed that it is not present in other stable assemblies involved in the selenocysteine pathway, namely the phosphoseryl-tRNA^Sec kinase (PSTK)-Sec-tRNA^Sec synthase (SEPSECS) complex and the tRNA^Sec-specific elongation factor (eEFSec) complex. Endoplasmic reticulum stress with dithiothreitol (DTT) or tunicamycin upon selenophosphate synthetase ablation in procyclic T. brucei cells led to a growth defect. On the other hand, only DTT presented a negative effect in bloodstream T. brucei expressing selenophosphate synthetase-RNAi. Furthermore, selenoprotein T (SELENOT) was dispensable for both forms of the parasite. Together, our data suggest a role for the T. brucei selenophosphate synthetase in the regulation of the parasite’s ER stress response.

Selenium is both a toxic compound and a micronutrient. As a micronutrient, it participates in the synthesis of specific proteins, selenoproteins, as the amino acid selenocysteine. The synthesis of selenocysteine is present in organisms ranging from bacteria to humans. The protist parasites of the Trypanosomatidae family, that cause major tropical diseases, conserve the complex machinery responsible for selenocysteine biosynthesis and incorporation in selenoproteins. However, this pathway has been considered dispensable for the parasitic protist cells. This has intrigued us, and lead to question that if maintained in the cell it should be under selective pressure and therefore be necessary. Also, extensive and dynamic protein-protein interactions must happen to deliver selenium-containing intermediates along the pathway in order to warrant efficient usage of biological selenium in the cell. In this study we have investigated the molecular interactions of different proteins involved in selenocysteine synthesis and its putative involvement in the endoplasmic reticulum redox homeostasis.

Sleeping Sickness: A Tale of Two Clocks

Sleeping sickness is caused by a eukaryotic unicellular parasite known to infect wild animals, cattle, and humans. It causes a fatal disease that disrupts many rhythmic physiological processes, including daily rhythms of hormonal secretion, temperature regulation, and sleep, all of which are under circadian (24-h) control. In this review, we summarize research on sleeping sickness parasite biology and the impact it has on host health. We also consider the possible evolutionary advantages of sleep and circadian deregulation for the parasite.

Replication origin location might contribute to genetic variability in Trypanosoma cruzi

Background

DNA replication in trypanosomatids operates in a uniquely challenging environment, since most of their genomes are constitutively transcribed. Trypanosoma cruzi, the etiological agent of Chagas disease, presents high variability in both chromosomes size and copy number among strains, though the underlying mechanisms are unknown.

Results

Here we have mapped sites of DNA replication initiation across the T. cruzi genome using Marker Frequency Analysis, which has previously only been deployed in two related trypanosomatids. The putative origins identified in T. cruzi show a notable enrichment of GC content, a preferential position at subtelomeric regions, coinciding with genes transcribed towards the telomeres, and a pronounced enrichment within coding DNA sequences, most notably in genes from the Dispersed Gene Family 1 (DGF-1).

Conclusions

These findings suggest a scenario where collisions between DNA replication and transcription are frequent, leading to increased genetic variability, as seen by the increase SNP levels at chromosome subtelomeres and in DGF-1 genes containing putative origins.

Novel structural CYP51 mutation in Trypanosoma cruzi associated with multidrug resistance to CYP51 inhibitors and reduced infectivity

Ergosterol biosynthesis inhibitors, such as posaconazole and ravuconazole, have been proposed as drug candidates for Chagas disease, a neglected infectious tropical disease caused by the protozoan parasite Trypanosoma cruzi. To understand better the mechanism of action and resistance to these inhibitors, a clone of the T. cruzi Y strain was cultured under intermittent and increasing concentrations of ravuconazole until phenotypic stability was achieved. The ravuconazole-selected clone exhibited loss in fitness in vitro when compared to the wild-type parental clone, as observed in reduced invasion capacity and slowed population growth in both mammalian and insect stages of the parasite. In drug activity assays, the resistant clone was above 300-fold more tolerant to ravuconazole than the sensitive parental clone, when the half-maximum effective concentration (EC₅₀) was considered. The resistant clones also showed reduced virulence in vivo, when compared to parental sensitive clones. Cross-resistance to posaconazole and other CYP51 inhibitors, but not to other antichagasic drugs that act independently of CYP51, such as benznidazole and nifurtimox, was also observed. A novel amino acid residue change, T297M, was found in the TcCYP51 gene in the resistant but not in the sensitive clones. The structural effects of the T297M, and of the previously described P355S residue changes, were modelled to understand their impact on interaction with CYP51 inhibitors.

•

A ravuconazole-resistant T. cruzi clone presented reduced in vitro and in vivo fitness.

•

The ravuconazole-resistant clone presented cross-resistance to other CYP51 inhibitors.

•

There was no cross-resistance to benznidazole and nifurtimox.

•

Resistance is associated with a novel structural mutation in the TcCYP51 protein.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL’s homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

The general mRNA exporters Mex67 and Mtr2 play distinct roles in nuclear export of tRNAs in Trypanosoma brucei

Transfer RNAs (tRNAs) are central players in protein synthesis, which in Eukarya need to be delivered from the nucleus to the cytoplasm by specific transport receptors, most of which belong to the evolutionarily conserved beta-importin family. Based on the available literature, we identified two candidates, Xpo-t and Xpo-5 for tRNA export in Trypanosoma brucei. However, down-regulation of expression of these genes did not disrupt the export of tRNAs to the cytoplasm. In search of alternative pathways, we tested the mRNA export complex Mex67-Mtr2, for a role in tRNA nuclear export, as described previously in yeast. Down-regulation of either exporter affected the subcellular distribution of tRNAs. However, contrary to yeast, TbMex67 and TbMtr2 accumulated different subsets of tRNAs in the nucleus. While TbMtr2 perturbed the export of all the tRNAs tested, silencing of TbMex67, led to the nuclear accumulation of tRNAs that are typically modified with queuosine. In turn, inhibition of tRNA nuclear export also affected the levels of queuosine modification in tRNAs. Taken together, the results presented demonstrate the dynamic nature of tRNA trafficking in T. brucei and its potential impact not only on the availability of tRNAs for protein synthesis but also on their modification status.

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.

African trypanosomes

African trypanosomes cause human African trypanosomiasis and animal African trypanosomiasis. They are transmitted by tsetse flies in sub-Saharan Africa. Although most famous for their mechanisms of immune evasion by antigenic variation, there have been recent important studies that illuminate important aspects of the biology of these parasites both in their mammalian host and during passage through their tsetse fly vector. This Primer overviews current research themes focused on these parasites and discusses how these biological insights and the development of new technologies to interrogate gene function are being used in the search for new approaches to control the parasite. The new insights into the biology of trypanosomes in their host and vector highlight that we are in a ‘golden age’ of discovery for these fascinating parasites.

Involvement in surface antigen expression by a moonlighting FG-repeat nucleoporin in trypanosomes

The nuclear pore complex is an ancient component of the eukaryotic cell. We show here that an FG nucleoporin, TbNup53b, in trypanosomes has an association with the splicing machinery and roles in gene expression, indicating that moonlighting roles for nucleoporins are highly ancient and present in the earliest eukaryotes.

Components of the nuclear periphery coordinate a multitude of activities, including macromolecular transport, cell-cycle progression, and chromatin organization. Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport, mRNA processing, and transcriptional regulation, and NPC components can define regions of high transcriptional activity in some organisms at the nuclear periphery and nucleoplasm. Lineage-specific features underpin several core nuclear functions and in trypanosomatids, which branched very early from other eukaryotes, unique protein components constitute the lamina, kinetochores, and parts of the NPCs. Here we describe a phenylalanine-glycine (FG)-repeat nucleoporin, TbNup53b, that has dual localizations within the nucleoplasm and NPC. In addition to association with nucleoporins, TbNup53b interacts with a known trans-splicing component, TSR1, and has a role in controlling expression of surface proteins including the nucleolar periphery-located, procyclin genes. Significantly, while several nucleoporins are implicated in intranuclear transcriptional regulation in metazoa, TbNup53b appears orthologous to components of the yeast/human Nup49/Nup58 complex, for which no transcriptional functions are known. These data suggest that FG-Nups are frequently co-opted to transcriptional functions during evolution and extend the presence of FG-repeat nucleoporin control of gene expression to trypanosomes, suggesting that this is a widespread and ancient eukaryotic feature, as well as underscoring once more flexibility within nucleoporin function.

Genome-wide mapping reveals conserved and diverged R-loop activities in the unusual genetic landscape of the African trypanosome genome

R-loops are stable RNA-DNA hybrids that have been implicated in transcription initiation and termination, as well as in telomere maintenance, chromatin formation, and genome replication and instability. RNA Polymerase (Pol) II transcription in the protozoan parasite Trypanosoma brucei is highly unusual: virtually all genes are co-transcribed from multigene transcription units, with mRNAs generated by linked trans-splicing and polyadenylation, and transcription initiation sites display no conserved promoter motifs. Here, we describe the genome-wide distribution of R-loops in wild type mammal-infective T. brucei and in mutants lacking RNase H1, revealing both conserved and diverged functions. Conserved localization was found at centromeres, rRNA genes and retrotransposon-associated genes. RNA Pol II transcription initiation sites also displayed R-loops, suggesting a broadly conserved role despite the lack of promoter conservation or transcription initiation regulation. However, the most abundant sites of R-loop enrichment were within the regions between coding sequences of the multigene transcription units, where the hybrids coincide with sites of polyadenylation and nucleosome-depletion. Thus, instead of functioning in transcription termination the most widespread localization of R-loops in T. brucei suggests a novel correlation with pre-mRNA processing. Finally, we find little evidence for correlation between R-loop localization and mapped sites of DNA replication initiation.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

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.

RNA-seq transcriptional profiling of Leishmania amazonensis reveals an arginase-dependent gene expression regulation

Background

Leishmania is a protozoan parasite that alternates its life cycle between the sand-fly vector and the mammalian host. This alternation involves environmental changes and leads the parasite to dynamic modifications in morphology, metabolism, cellular signaling and regulation of gene expression to allow for a rapid adaptation to new conditions. The L-arginine pathway in L. amazonensis is important during the parasite life cycle and interferes in the establishment and maintenance of the infection in mammalian macrophages. Host arginase is an immune-regulatory enzyme that can reduce the production of nitric oxide by activated macrophages, directing the availability of L-arginine to the polyamine pathway, resulting in parasite replication. In this work, we performed transcriptional profiling to identify differentially expressed genes in L. amazonensis wild-type (La-WT) versus L. amazonensis arginase knockout (La-arg^-) promastigotes and axenic amastigotes.

Methodology/Principal findings

A total of 8253 transcripts were identified in La-WT and La-arg^- promastigotes and axenic amastigotes, about 60% of them codifying hypothetical proteins and 443 novel transcripts, which did not match any previously annotated genes. Our RNA-seq data revealed that 85% of genes were constitutively expressed. The comparison of transcriptome and metabolome data showed lower levels of arginase and higher levels of glutamate-5-kinase in La-WT axenic amastigotes compared to promastigotes. The absence of arginase activity in promastigotes increased the levels of pyrroline 5-carboxylate reductase, but decreased the levels of arginosuccinate synthase, pyrroline 5-carboxylate dehydrogenase, acetylornithine deacetylase and spermidine synthase transcripts levels. These observations can explain previous metabolomic data pointing to the increase of L-arginine, citrulline and L-glutamate and reduction of aspartate, proline, ornithine and putrescine. Altogether, these results indicate that arginase activity is important in Leishmania gene expression modulation during differentiation and adaptation to environmental changes. Here, we confirmed this hypothesis with the identification of differential gene expression of the enzymes involved in biosynthesis of amino acids, arginine and proline metabolism and arginine biosynthesis.

Conclusions/Significance

All data provided information about the transcriptomic profiling and the expression levels of La-WT and La-arg^- promastigotes and axenic amastigotes. These findings revealed the importance of arginase in parasite survival and differentiation, and indicated the existence of a coordinated response in the absence of arginase activity related to arginine and polyamine pathways.

Leishmania are auxotrophic for many essential nutrients, including amino acids. In this way, the parasite needs to uptake the amino acids from the environment. The uptake of amino acids is mediated by amino acid transporters that are unique for Leishmania. As part of polyamine pathway, the arginase converts L-arginine to ornithine and furthermore to putrescine, products which are essential for parasite growth. On the other hand, the absence of arginase activity could alter the metabolism of the parasite to surpass the external signals during the life cycle and the fate of infection. The transcriptional profiling of La-WT and La-arg^- promastigotes and axenic amastigotes revealed 8253 transcripts, 60% encoding hypothetical proteins and 443 novel transcripts. In addition, our data revealed that 85% of the genes were constitutively expressed. Among the 15% (1268 genes) of the differentially expressed genes, we identified genes up- and down-regulated comparing the transcript abundance from different life cycle stages of the parasite and in the presence or absence of arginase. We also combined the transcriptional with metabolic profile that revealed a proportional correlation between enzyme and metabolites in the polyamine pathway. The differentiation of promastigotes to amastigotes alters the expression of enzymes from polyamines biosynthesis, which modulates ornithine, L-glutamate, proline and putrescine levels. In addition, the absence of arginase activity increased the levels of L-arginine, citrulline and L-glutamate and decreased the levels of aspartate, proline, ornithine and putrescine in promastigotes by differential modulation of genes involved in its metabolism. Altogether these data provided additional insights into how Leishmania is able to modulate its biological functions in the presence or absence of arginase activity to survive during environmental changes.

Genomic analysis of Isometamidium Chloride resistance in Trypanosoma congolense

Isometamidium Chloride (ISM) is one of the principal drugs used to counteract Trypanosoma congolense infection in livestock, both as a prophylactic as well as a curative treatment. However, numerous cases of ISM resistance have been reported in different African regions, representing a significant constraint in the battle against Animal African Trypanosomiasis. In order to identify genetic signatures associated with ISM resistance in T. congolense, the sensitive strain MSOROM7 was selected for induction of ISM resistance in a murine host. Administered ISM concentrations in immune-suppressed mice were gradually increased from 0.001 mg/kg to 1 mg/kg, the maximal dose used in livestock. As a result, three independent MSOROM7 lines acquired full resistance to this concentration after five months of induction, and retained this full resistant phenotype following a six months period without drug pressure. In contrast, parasites did not acquire ISM resistance in immune-competent animals, even after more than two years under ISM pressure, suggesting that the development of full ISM resistance is strongly enhanced when the host immune response is compromised. Genomic analyses comparing the ISM resistant lines with the parental sensitive line identified shifts in read depth at heterozygous loci in genes coding for different transporters and transmembrane products, and several of these shifts were also found within natural ISM resistant isolates. These findings suggested that the transport and accumulation of ISM inside the resistant parasites may be modified, which was confirmed by flow cytometry and ex vivo ISM uptake assays that showed a decrease in the accumulation of ISM in the resistant parasites.

Cell cycle stage-specific transcriptional activation of cyclins mediated by HAT2-dependent H4K10 acetylation of promoters in Leishmania donovani

Chromatin modifications affect several processes. In investigating the Leishmania donovani histone acetyltransferase HAT2, using in vitro biochemical assays and HAT2-heterozygous genomic knockout we found the constitutively nuclear HAT2 acetylated histone H4K10 in vitro and in vivo. HAT2 was essential. HAT2-depleted cells displayed growth and cell cycle defects, and poor survival in host cells. Real time PCR and DNA microarray analyses, as well as rescue experiments, revealed that downregulation of cyclins CYC4 and CYC9 were responsible for S phase and G2/M defects of HAT2-depleted cells respectively. Leishmania genes are arranged in unidirectional clusters, and clustered genes are coordinately transcribed as long polycistronic units, typically from divergent strand switch regions (dSSRs) which initiate transcription bidirectionally on opposite strands. In investigating the mechanism by which CYC4 and CYC9 expression levels are reduced in HAT2-depleted cells without other genes in their polycistronic transcription units being coordinately downregulated, we found using reporter assays that CYC4 and CYC9 have their own specific promoters. Chromatin immunoprecipitation assays with H4acetylK10 antibodies and real time PCR analyses of RNA suggested these gene-specific promoters were activated in cell cycle-dependent manner. Nuclear run-on analyses confirmed that CYC4 and CYC9 were transcriptionally activated from their own promoters at specific cell cycle stages. Thus, there are two tiers of gene regulation. Transcription of polycistronic units primarily initiates at dSSRs, and this most likely occurs constitutively. A subset of genes have their own promoters, at least some of which are activated in a cell-cycle dependent manner. This second tier of regulation is more sensitive to H4K10 acetylation levels, resulting in downregulation of expression in HAT2-depleted cells. This report presents the first data pointing to cell cycle-specific activation of promoters in trypanosomatids, thus uncovering new facets of gene regulation in this parasite family.

The proteome and transcriptome of the infectious metacyclic form of Trypanosoma brucei define quiescent cells primed for mammalian invasion

The infectious metacyclic forms of Trypanosoma brucei result from a complex development in the tsetse fly vector. When they infect mammals, they cause African sleeping sickness in humans. Due to scarcity of biological material and difficulties of the tsetse fly as an experimental system, very limited information is available concerning the gene expression profile of metacyclic forms. We used an in vitro system based on expressing the RNA binding protein 6 to obtain infectious metacyclics and determined their protein and mRNA repertoires by mass-spectrometry (MS) based proteomics and mRNA sequencing (RNA-Seq) in comparison to non-infectious procyclic trypanosomes. We showed that metacyclics are quiescent cells, and propose this influences the choice of a monocistronic variant surface glycoprotein expression site. Metacyclics have a largely bloodstream-form type transcriptome, and thus are programmed to translate a bloodstream-form type proteome upon entry into the mammalian host and resumption of cell division. Genes encoding cell surface components showed the largest changes between procyclics and metacyclics, observed at both the transcript and protein levels. Genes encoding metabolic enzymes exhibited expression in metacyclics with features of both procyclic and bloodstream forms, suggesting that this intermediate-type metabolism is dictated by the availability of nutrients in the tsetse fly vector.

Well-positioned nucleosomes punctuate polycistronic pol II transcription units and flank silent VSG gene arrays in Trypanosoma brucei

Background

The compaction of DNA in chromatin in eukaryotes allowed the expansion of genome size and coincided with significant evolutionary diversification. However, chromatin generally represses DNA function, and mechanisms coevolved to regulate chromatin structure and its impact on DNA. This included the selection of specific nucleosome positions to modulate accessibility to the DNA molecule. Trypanosoma brucei, a member of the Excavates supergroup, falls in an ancient evolutionary branch of eukaryotes and provides valuable insight into the organization of chromatin in early genomes.

Results

We have mapped nucleosome positions in T. brucei and identified important differences compared to other eukaryotes: The RNA polymerase II initiation regions in T. brucei do not exhibit pronounced nucleosome depletion, and show little evidence for defined −1 and +1 nucleosomes. In contrast, a well-positioned nucleosome is present directly on the splice acceptor sites within the polycistronic transcription units. The RNA polyadenylation sites were depleted of nucleosomes, with a single well-positioned nucleosome present immediately downstream of the predicted sites. The regions flanking the silent variant surface glycoprotein (VSG) gene cassettes showed extensive arrays of well-positioned nucleosomes, which may repress cryptic transcription initiation. The silent VSG genes themselves exhibited a less regular nucleosomal pattern in both bloodstream and procyclic form trypanosomes. The DNA replication origins, when present within silent VSG gene cassettes, displayed a defined nucleosomal organization compared with replication origins in other chromosomal core regions.

Conclusions

Our results indicate that some organizational features of chromatin are evolutionarily ancient, and may already have been present in the last eukaryotic common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0121-9) contains supplementary material, which is available to authorized users.

Transcriptomic analysis reveals metabolic switches and surface remodeling as key processes for stage transition in Trypanosoma cruzi

American trypanosomiasis is a chronic and endemic disease which affects millions of people. Trypanosoma cruzi, its causative agent, has a life cycle that involves complex morphological and functional transitions, as well as a variety of environmental conditions. This requires a tight regulation of gene expression, which is achieved mainly by post-transcriptional regulation. In this work we conducted an RNAseq analysis of the three major life cycle stages of T. cruzi: amastigotes, epimastigotes and trypomastigotes. This analysis allowed us to delineate specific transcriptomic profiling for each stage, and also to identify those biological processes of major relevance in each state. Stage specific expression profiling evidenced the plasticity of T. cruzi to adapt quickly to different conditions, with particular focus on membrane remodeling and metabolic shifts along the life cycle. Epimastigotes, which replicate in the gut of insect vectors, showed higher expression of genes related to energy metabolism, mainly Krebs cycle, respiratory chain and oxidative phosphorylation related genes, and anabolism related genes associated to nucleotide and steroid biosynthesis; also, a general down-regulation of surface glycoprotein coding genes was seen at this stage. Trypomastigotes, living extracellularly in the bloodstream of mammals, express a plethora of surface proteins and signaling genes involved in invasion and evasion of immune response. Amastigotes mostly express membrane transporters and genes involved in regulation of cell cycle, and also express a specific subset of surface glycoprotein coding genes. In addition, these results allowed us to improve the annotation of the Dm28c genome, identifying new ORFs and set the stage for construction of networks of co-expression, which can give clues about coded proteins of unknown functions.

Lineage-specific proteins essential for endocytosis in trypanosomes

Clathrin-mediated endocytosis (CME) is the most evolutionarily ancient endocytic mechanism known, and in many lineages the sole mechanism for internalisation. Significantly, in mammalian cells CME is responsible for the vast bulk of endocytic flux and has likely undergone multiple adaptations to accommodate specific requirements by individual species. In African trypanosomes, we previously demonstrated that CME is independent of the AP-2 adaptor protein complex, that orthologues to many of the animal and fungal CME protein cohort are absent, and that a novel, trypanosome-restricted protein cohort interacts with clathrin and drives CME. Here, we used a novel cryomilling affinity isolation strategy to preserve transient low-affinity interactions, giving the most comprehensive trypanosome clathrin interactome to date. We identified the trypanosome AP-1 complex, Trypanosoma brucei (Tb)EpsinR, several endosomal SNAREs plus orthologues of SMAP and the AP-2 associated kinase AAK1 as interacting with clathrin. Novel lineage-specific proteins were identified, which we designate TbCAP80 and TbCAP141. Their depletion produced extensive defects in endocytosis and endomembrane system organisation, revealing a novel molecular pathway subtending an early-branching and highly divergent form of CME, which is conserved and likely functionally important across the kinetoplastid parasites.

Summary: Endocytosis is a vital process in most cells, and here we identify important proteins required for this process in trypanosomes. Significantly, these are unique and not present in animals, fungi or plants.

Comprehensive Identification of mRNA-Binding Proteins of Leishmania donovani by Interactome Capture

Leishmania are unicellular eukaryotes responsible for leishmaniasis in humans. Like other trypanosomatids, leishmania regulate protein coding gene expression almost exclusively at the post-transcriptional level with the help of RNA binding proteins (RBPs). Due to the presence of polycystronic transcription units, leishmania do not regulate RNA polymerase II-dependent transcription initiation. Recent evidence suggests that the main control points in gene expression are mRNA degradation and translation. Protein-RNA interactions are involved in every aspect of RNA biology, such as mRNA splicing, polyadenylation, localization, degradation, and translation. A detailed picture of these interactions would likely prove to be highly informative in understanding leishmania biology and virulence. We developed a strategy involving covalent UV cross-linking of RBPs to mRNA in vivo, followed by interactome capture using oligo(dT) magnetic beads to define comprehensively the mRNA interactome of growing L. donovani amastigotes. The protein mass spectrometry analysis of captured proteins identified 79 mRNA interacting proteins which withstood very stringent washing conditions. Strikingly, we found that 49 of these mRNA interacting proteins had no orthologs or homologs in the human genome. Consequently, these may represent high quality candidates for selective drug targeting leading to novel therapeutics. These results show that this unbiased, systematic strategy has the promise to be applicable to study the mRNA interactome during various biological settings such as metabolic changes, stress (low pH environment, oxidative stress and nutrient deprivation) or drug treatment.

Genome wide occurrence and insertion preferences of INGI/RIME and SLACS CRE transposable elements in Trypanosoma brucei

Retrotransposons play significant role in genome remodelling of T. brucei and about 5% of its genome consists of retrotransposons including INGI/RIME elements. INGI is one of the dispersed repetitive elements in T. brucei genome which is found distributed throughout all the chromosomes. SLACS (Spliced Leader Associated Conserved Sequence) however, is more conserved in its nature and lacks the typical poly-distributional pattern of LINE like transposons. We have found total 589 copies of these TEs with only 17.06 % (104 copies out of 589) copies with both ends intact thus showing a majority of truncated copies. Complete SLACS CRE were found only on chromosome 9, whereas, complete INGI/RIME were more ubiquitously distributed. The 50 - 70 bp upstream flanking sequence of these elements shows suitable nucleotide biophysical properties to favor transposition.

Bromodomains in Protozoan Parasites: Evolution, Function, and Opportunities for Drug Development

Parasitic infections remain one of the most pressing global health concerns of our day, affecting billions of people and producing unsustainable economic burdens. The rise of drug-resistant parasites has created an urgent need to study their biology in hopes of uncovering new potential drug targets. It has been established that disrupting gene expression by interfering with lysine acetylation is detrimental to survival of apicomplexan (Toxoplasma gondii and Plasmodium spp.) and kinetoplastid (Leishmania spp. and Trypanosoma spp.) parasites. As "readers" of lysine acetylation, bromodomain proteins have emerged as key gene expression regulators and a promising new class of drug target. Here we review recent studies that demonstrate the essential roles played by bromodomain-containing proteins in parasite viability, invasion, and stage switching and present work showing the efficacy of bromodomain inhibitors as novel antiparasitic agents. In addition, we performed a phylogenetic analysis of bromodomain proteins in representative pathogens, some of which possess unique features that may be specific to parasite processes and useful in future drug development.

Trypanosome outer kinetochore proteins suggest conservation of chromosome segregation machinery across eukaryotes

The kinetochore complex is conserved across many eukaryotes, but the protozoan lineage Kinetoplastida builds kinetochores from components without apparent homology to models. D’Archivio and Wickstead describe a new family of proteins with homology to outer kinetochore proteins Ndc80 and Nuf2 that defines the outer kinetochore of trypanosomes, suggesting that all eukaryotes have divergent versions of a universal kinetochore machine.

Kinetochores are multiprotein complexes that couple eukaryotic chromosomes to the mitotic spindle to ensure proper segregation. The model for kinetochore assembly is conserved between humans and yeast, and homologues of several components are widely distributed in eukaryotes, but key components are absent in some lineages. The recent discovery in a lineage of protozoa called kinetoplastids of unconventional kinetochores with no apparent homology to model organisms suggests that more than one system for eukaryotic chromosome segregation may exist. In this study, we report a new family of proteins distantly related to outer kinetochore proteins Ndc80 and Nuf2. The family member in kinetoplastids, KKT-interacting protein 1 (KKIP1), associates with the kinetochore, and its depletion causes severe defects in karyokinesis, loss of individual chromosomes, and gross defects in spindle assembly or stability. Immunopurification of KKIP1 from stabilized kinetochores identifies six further components, which form part of a trypanosome outer kinetochore complex. These findings suggest that kinetochores in organisms such as kinetoplastids are built from a divergent, but not ancestrally distinct, set of components and that Ndc80/Nuf2-like proteins are universal in eukaryotic division.

The Major Protein Arginine Methyltransferase in Trypanosoma brucei Functions as an Enzyme-Prozyme Complex

Prozymes are catalytically inactive enzyme paralogs that dramatically stimulate the function of weakly active enzymes through complex formation. The two prozymes described to date reside in the polyamine biosynthesis pathway of the human parasite Trypanosoma brucei, an early branching eukaryote that lacks transcriptional regulation and regulates its proteome through posttranscriptional and posttranslational means. Arginine methylation is a common posttranslational modification in eukaryotes catalyzed by protein arginine methyltransferases (PRMTs) that are typically thought to function as homodimers. We demonstrate that a major T. brucei PRMT, TbPRMT1, functions as a heterotetrameric enzyme-prozyme pair. The inactive PRMT paralog, TbPRMT1^PRO, is essential for catalytic activity of the TbPRMT1^ENZ subunit. Mutational analysis definitively demonstrates that TbPRMT1^ENZ is the cofactor-binding subunit and carries all catalytic activity of the complex. Our results are the first demonstration of an obligate heteromeric PRMT, and they suggest that enzyme-prozyme organization is expanded in trypanosomes as a posttranslational means of enzyme regulation.

Nuclear DNA replication and repair in parasites of the genus Leishmania: Exploiting differences to develop innovative therapeutic approaches

Leishmaniasis is a common tropical disease that affects mainly poor people in underdeveloped and developing countries. This largely neglected infection is caused by Leishmania spp, a parasite from the Trypanosomatidae family. This parasitic disease has different clinical manifestations, ranging from localized cutaneous to more harmful visceral forms. The main limitations of the current treatments are their high cost, toxicity, lack of specificity, and long duration. Efforts to improve treatments are necessary to deal with this infectious disease. Many approved drugs to combat diseases as diverse as cancer, bacterial, or viral infections take advantage of specific features of the causing agent or of the disease. Recent evidence indicates that the specific characteristics of the Trypanosomatidae replication and repair machineries could be used as possible targets for the development of new treatments. Here, we review in detail the molecular mechanisms of DNA replication and repair regulation in trypanosomatids of the genus Leishmania and the drugs that could be useful against this disease.

Co-dependence between trypanosome nuclear lamina components in nuclear stability and control of gene expression

The nuclear lamina is a filamentous structure subtending the nuclear envelope and required for chromatin organization, transcriptional regulation and maintaining nuclear structure. The trypanosomatid coiled-coil NUP-1 protein is a lamina component functionally analogous to lamins, the major lamina proteins of metazoa. There is little evidence for shared ancestry, suggesting the presence of a distinct lamina system in trypanosomes. To find additional trypanosomatid lamina components we identified NUP-1 interacting proteins by affinity capture and mass-spectrometry. Multiple components of the nuclear pore complex (NPC) and a second coiled-coil protein, which we termed NUP-2, were found. NUP-2 has a punctate distribution at the nuclear periphery throughout the cell cycle and is in close proximity to NUP-1, the NPCs and telomeric chromosomal regions. RNAi-mediated silencing of NUP-2 leads to severe proliferation defects, gross alterations to nuclear structure, chromosomal organization and nuclear envelope architecture. Further, transcription is altered at telomere-proximal variant surface glycoprotein (VSG) expression sites (ESs), suggesting a role in controlling ES expression, although NUP-2 silencing does not increase VSG switching. Transcriptome analysis suggests specific alterations to Pol I-dependent transcription. NUP-1 is mislocalized in NUP-2 knockdown cells and vice versa, implying that NUP-1 and NUP-2 form a co-dependent network and identifying NUP-2 as a second trypanosomatid nuclear lamina component.

Comparative Transcriptome Profiling of Human Foreskin Fibroblasts Infected with the Sylvio and Y Strains of Trypanosoma cruzi

Trypanosoma cruzi, the causative agent of Chagas Disease, is phylogeneticaly distributed into nearly identical genetic strains which show divergent clinical presentations including differences in rates of cardiomyopathy in humans, different vector species and transmission cycles, differential congenital transmission in a mouse model, and differing immune and heart inflammation response in dogs. The population structure of these strains divides into two groups, which are geographically and clinically distinct. The aim of this study was to compare the transcriptome of two strains of T. cruzi, Sylvio vs. Y, to identify differences in expression that could account for clinical and biochemical differences. We collected and sequenced RNA from T. cruzi-infected and control Human Foreskin Fibroblasts at three timepoints. Differential expression analysis identified gene expression different timepoints in Sylvio infections, and between Sylvio and Y infections in both parasite and host. The Sylvio strain parasite and the host response to Sylvio infection largely mirrored the host-pathogen interaction seen in our previous Y strain work. IL-8 was more highly expressed in Sylvio-infected HFFs than in Y-infected HFFs.

The nuclear envelope and gene organization in parasitic protozoa: Specializations associated with disease

The parasitic protozoa Trypanosoma brucei and Plasmodium falciparum are lethal human parasites that have developed elegant strategies of immune evasion by antigenic variation. Despite the vast evolutionary distance between the two taxa, both parasites employ strict monoallelic expression of their membrane proteins, variant surface glycoproteins in Trypanosomes and the var, rif and stevor genes in Plasmodium, in order to evade their host's immune system. Additionally, both telomeric location and epigenetic controls are prominent features of these membrane proteins. As such, telomeres, chromatin structure and nuclear organization all contribute to control of gene expression and immune evasion. Here, we discuss the importance of epigenetics and sub-nuclear context for the survival of these disease-causing parasites.