Histone acetylations mark origins of polycistronic transcription in Leishmania major

Characterization of Tau95 led to the identification of a four-subunit TFIIIC complex in trypanosomatid parasites

RNA polymerase III (RNAP III) synthetizes small essential non-coding RNA molecules such as tRNAs and 5S rRNA. In yeast and vertebrates, RNAP III needs general transcription factors TFIIIA, TFIIIB, and TFIIIC to initiate transcription. TFIIIC, composed of six subunits, binds to internal promoter elements in RNAP III-dependent genes. Limited information is available about RNAP III transcription in the trypanosomatid protozoa Trypanosoma brucei and Leishmania major, which diverged early from the eukaryotic lineage. Analyses of the first published draft of the trypanosomatid genome sequences failed to recognize orthologs of any of the TFIIIC subunits, suggesting that this transcription factor is absent in these parasites. However, a putative TFIIIC subunit was recently annotated in the databases. Here we characterize this subunit in T. brucei and L. major and demonstrate that it corresponds to Tau95. In silico analyses showed that both proteins possess the typical Tau95 sequences: the DNA binding region and the dimerization domain. As anticipated for a transcription factor, Tau95 localized to the nucleus in insect forms of both parasites. Chromatin immunoprecipitation (ChIP) assays demonstrated that Tau95 binds to tRNA and U2 snRNA genes in T. brucei. Remarkably, by performing tandem affinity purifications we identified orthologs of TFIIIC subunits Tau55, Tau131, and Tau138 in T. brucei and L. major. Thus, contrary to what was assumed, trypanosomatid parasites do possess a TFIIIC complex. Other putative interacting partners of Tau95 were identified in T. brucei and L. major. KEY POINTS: • A four-subunit TFIIIC complex is present in T. brucei and L. major • TbTau95 associates with tRNA and U2 snRNA genes • Putative interacting partners of Tau95 might include some RNAP II regulators.

Mono-allelic epigenetic regulation of bi-directional silencing of RNA Polymerase II polycistronic transcription initiation in Trypanosoma brucei

Unique for a eukaryote, protein-coding genes in trypanosomes are arranged in polycistronic units (PTUs). This genome arrangement has led to a model where Pol II transcription of PTUs is unregulated and that changes in gene expression are entirely post-transcriptional. Trypanosoma brucei brucei is unable to infect humans because of its susceptibility to an innate immune complex, trypanosome lytic factor (TLF) in the circulation of humans. The initial step in TLF mediated lysis of T.b.brucei requires high affinity haptoglobin/hemoglobin receptor (HpHbR) binding. Here we demonstrate that by in vitro selection with TLF, resistance is obtained in a stepwise process correlating with loss of HpHbR expression at an allelic level. RNA-seq, Pol II ChIP and run-on analysis indicate HpHbR silencing is at the transcriptional level, where loss of Pol II binding at the promoter region specifically shuts down transcription of the HpHbR containing gene cluster and the adjacent opposing gene cluster. Reversible transcriptional silencing of the divergent PTUs correlates with DNA base J modification of the shared promoter region. Therefore, we show that epigenetic mechanisms, including base J modification, are involved in regulating gene expression via Pol II transcription initiation of gene clusters in a mono-allelic fashion. These findings suggest epigenetic chromatin-based regulation of gene expression is deeply conserved among eukaryotes, including early divergent eukaryotes that rely on polycistronic transcription.

IMPORTANCE

The single-cell parasite Trypanosoma brucei causes lethal diseases in both humans and livestock. T. brucei undergoes multiple developmental changes to adapt in different environments during its digenetic life cycle. With protein-coding genes organized as polycistronic transcription and apparent absence of promoter-mediated regulation of transcription initiation, it is believed that developmental gene regulation in trypanosomes is essentially post-transcriptional. In this study, we found reversible Pol II transcriptional silencing of two adjacent polycistronic gene arrays that correlates with the novel DNA base J modification of the shared promoter region. Our findings support epigenetic regulation of Pol II transcription initiation as a viable mechanism of gene expression control in T. brucei . This has implications for our understanding how trypanosomes utilize polycistronic genome organization to regulate gene expression during its life cycle.

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.

Protein phosphatase PP1 regulation of RNA polymerase II transcription termination and allelic exclusion of VSG genes in trypanosomes

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, complicating the analysis of PP1 function in termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei and dephosphorylation of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated RPB1 accompanied by readthrough transcription and aberrant transcription of the chromosome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Clinical polymorphism of zoonotic cutaneous leishmaniasis: combination of the clinical and the parasitological diagnosis

Zoonotic cutaneous leishmaniasis (ZCL) is a neglected tropical disease caused by Leishmania (L.) major. This zoonosis is characterized by a broad-spectrum clinical polymorphism and may be underestimated and poorly treated since it is a simulator of various dermatoses. The aim of our study was to analyze the clinical polymorphism of patients with ZCL. A total of 142 patients with confirmed CL based on the microscopic examination of skin lesion biopsies were included in this study. Molecular typing of Leishmania species revealed that all patients were infected with L. major. In total, 14 clinical forms were observed. Six were typical and eight were atypical. The typical ZCL forms are grouped as follows: papular (26.76%), ulcero-crusted (26.05%), ulcerated (13.38%), impetiginous (9.86%), nodular (9.15%), and papulo-nodular (5.63%) lesions. In atypical ZCL forms, we described erythematous (2.81%), erysipeloid (1.4%), sporotrichoid, (1.4%), keratotic (0.7%) lupoid (0.7%), lichenoid (0.7%), psoriasiform (0.7%), and zosteriform (0.7%) lesions. Here, the lichenoid and the keratotic forms caused by L. major were reported for the first time in Tunisia. These findings will help physicians to be aware of the unusual lesions of ZCL that could be confused with other dermatological diseases. For this reason, it will be necessary to improve the diagnosis of CL especially in endemic areas. Such large clinical polymorphism caused by L. major may be the result of a complex association between the vector microbiota, the parasite, and the host immune state, and further studies should be carried out in order to reveal the mechanisms involved in clinical polymorphism of ZCL.

Leishmania PNUTS discriminates between PP1 catalytic subunits through an RVxF-ΦΦ-F motif and polymorphisms in the PP1 C-tail and catalytic domain

Phosphoprotein phosphatase 1 (PP1) associates with specific regulatory subunits to achieve, among other functions, substrate selectivity. Among the eight PP1 isotypes in Leishmania, PP1-8e associates with the regulatory protein PNUTS along with the structural factors JBP3 and Wdr82 in the PJW/PP1 complex that modulates RNA polymerase II (pol II) phosphorylation and transcription termination. Little is known regarding interactions involved in PJW/PP1 complex formation, including how PP1-8e is the selective isotype associated with PNUTS. Here, we show that PNUTS uses an established RVxF-ΦΦ-F motif to bind the PP1 catalytic domain with similar interfacial interactions as mammalian PP1-PNUTS and noncanonical motifs. These atypical interactions involve residues within the PP1-8e catalytic domain and N and C terminus for isoform-specific regulator binding. This work advances our understanding of PP1 isoform selectivity and reveals key roles of PP1 residues in regulator binding. We also explore the role of PNUTS as a scaffold protein for the complex by identifying the C-terminal region involved in binding JBP3 and Wdr82 and impact of PNUTS on the stability of complex components and function in pol II transcription in vivo. Taken together, these studies provide a potential mechanism where multiple motifs within PNUTS are used combinatorially to tune binding affinity to PP1, and the C terminus for JBP3 and Wdr82 association, in the Leishmania PJW/PP1 complex. Overall, our data provide insights in the formation of the PJW/PP1 complex involved in regulating pol II transcription in divergent protozoans where little is understood.

Bromodomain Factor 5 as a Target for Antileishmanial Drug Discovery

Leishmaniases are a collection of neglected tropical diseases caused by kinetoplastid parasites in the genus Leishmania. Current chemotherapies are severely limited, and the need for new antileishmanials is of pressing international importance. Bromodomains are epigenetic reader domains that have shown promising therapeutic potential for cancer therapy and may also present an attractive target to treat parasitic diseases. Here, we investigate Leishmania donovani bromodomain factor 5 (LdBDF5) as a target for antileishmanial drug discovery. LdBDF5 contains a pair of bromodomains (BD5.1 and BD5.2) in an N-terminal tandem repeat. We purified recombinant bromodomains of L. donovani BDF5 and determined the structure of BD5.2 by X-ray crystallography. Using a histone peptide microarray and fluorescence polarization assay, we identified binding interactions of LdBDF5 bromodomains with acetylated peptides derived from histones H2B and H4. In orthogonal biophysical assays including thermal shift assays, fluorescence polarization, and NMR, we showed that BDF5 bromodomains bind to human bromodomain inhibitors SGC-CBP30, bromosporine, and I-BRD9; moreover, SGC-CBP30 exhibited activity against Leishmania promastigotes in cell viability assays. These findings exemplify the potential BDF5 holds as a possible drug target in Leishmania and provide a foundation for the future development of optimized antileishmanial compounds targeting this epigenetic reader protein.

Protein Phosphatase PP1 Regulation of Pol II Phosphorylation is Linked to Transcription Termination and Allelic Exclusion of VSG Genes and TERRA in Trypanosomes

The genomes of Leishmania and trypanosomes are organized into polycistronic transcription units flanked by a modified DNA base J involved in promoting RNA polymerase II (Pol II) termination. We recently characterized a Leishmania complex containing a J-binding protein, PP1 protein phosphatase 1, and PP1 regulatory protein (PNUTS) that controls transcription termination potentially via dephosphorylation of Pol II by PP1. While T. brucei contains eight PP1 isoforms, none purified with the PNUTS complex, suggesting a unique PP1-independent mechanism of termination. We now demonstrate that the PP1-binding motif of TbPNUTS is required for function in termination in vivo and that TbPP1-1 modulates Pol II termination in T. brucei involving dephosphorylation of the C-terminal domain of the large subunit of Pol II. PP1-1 knock-down results in increased cellular levels of phosphorylated large subunit of Pol II accompanied by readthrough transcription and pervasive transcription of the entire genome by Pol II, including Pol I transcribed loci that are typically silent, such as telomeric VSG expression sites involved in antigenic variation and production of TERRA RNA. These results provide important insights into the mechanism underlying Pol II transcription termination in primitive eukaryotes that rely on polycistronic transcription and maintain allelic exclusion of VSG genes.

Knockout of protein phosphatase 1 in Leishmania major reveals its role during RNA polymerase II transcription termination

The genomes of kinetoplastids are organized into polycistronic transcription units that are flanked by a modified DNA base (base J, beta-D-glucosyl-hydroxymethyluracil). Previous work established a role of base J in promoting RNA polymerase II (Pol II) termination in Leishmania major and Trypanosoma brucei. We recently identified a PJW/PP1 complex in Leishmania containing a J-binding protein (JBP3), PP1 phosphatase 1, PP1 interactive-regulatory protein (PNUTS) and Wdr82. Analyses suggested the complex regulates transcription termination by recruitment to termination sites via JBP3-base J interactions and dephosphorylation of proteins, including Pol II, by PP1. However, we never addressed the role of PP1, the sole catalytic component, in Pol II transcription termination. We now demonstrate that deletion of the PP1 component of the PJW/PP1 complex in L. major, PP1-8e, leads to readthrough transcription at the 3'-end of polycistronic gene arrays. We show PP1-8e has in vitro phosphatase activity that is lost upon mutation of a key catalytic residue and associates with PNUTS via the conserved RVxF motif. Additionally, purified PJW complex with associated PP1-8e, but not complex lacking PP1-8e, led to dephosphorylation of Pol II, suggesting a direct role of PNUTS/PP1 holoenzymes in regulating transcription termination via dephosphorylating Pol II in the nucleus.

Advances in Protozoan Epigenetic Targets and Their Inhibitors for the Development of New Potential Drugs

Protozoan parasite diseases cause significant mortality and morbidity worldwide. Factors such as climate change, extreme poverty, migration, and a lack of life opportunities lead to the propagation of diseases classified as tropical or non-endemic. Although there are several drugs to combat parasitic diseases, strains resistant to routinely used drugs have been reported. In addition, many first-line drugs have adverse effects ranging from mild to severe, including potential carcinogenic effects. Therefore, new lead compounds are needed to combat these parasites. Although little has been studied regarding the epigenetic mechanisms in lower eukaryotes, it is believed that epigenetics plays an essential role in vital aspects of the organism, from controlling the life cycle to the expression of genes involved in pathogenicity. Therefore, using epigenetic targets to combat these parasites is foreseen as an area with great potential for development. This review summarizes the main known epigenetic mechanisms and their potential as therapeutics for a group of medically important protozoal parasites. Different epigenetic mechanisms are discussed, highlighting those that can be used for drug repositioning, such as histone post-translational modifications (HPTMs). Exclusive parasite targets are also emphasized, including the base J and DNA 6 mA. These two categories have the greatest potential for developing drugs to treat or eradicate these diseases.

The paradigm of intracellular parasite survival and drug resistance in leishmanial parasite through genome plasticity and epigenetics: Perception and future perspective

Leishmania is an intracellular, zoonotic, kinetoplastid eukaryote with more than 1.2 million cases all over the world. The leishmanial chromosomes are divided into polymorphic chromosomal ends, conserved central domains, and antigen-encoding genes found in telomere-proximal regions. The genome flexibility of chromosomal ends of the leishmanial parasite is known to cause drug resistance and intracellular survival through the evasion of host defense mechanisms. Therefore, in this review, we discuss the plasticity of Leishmania genome organization which is the primary cause of drug resistance and parasite survival. Moreover, we have not only elucidated the causes of such genome plasticity which includes aneuploidy, epigenetic factors, copy number variation (CNV), and post-translation modification (PTM) but also highlighted their impact on drug resistance and parasite survival.

Localization of Epigenetic Markers in Leishmania Chromatin

Eukaryotes use histone variants and post-translation modifications (PTMs), as well as DNA base modifications, to regulate DNA replication/repair, chromosome condensation, and gene expression. Despite the unusual organization of their protein-coding genes into large polycistronic transcription units (PTUs), trypanosomatid parasites also employ a “histone code” to control these processes, but the details of this epigenetic code are poorly understood. Here, we present the results of experiments designed to elucidate the distribution of histone variants and PTMs over the chromatin landscape of Leishmania tarentolae. These experiments show that two histone variants (H2A.Z and H2B.V) and three histone H3 PTMs (H3K4me3, H3K16ac, and H3K76me3) are enriched at transcription start sites (TSSs); while a histone variant (H3.V) and the trypanosomatid-specific hyper-modified DNA base J are located at transcription termination sites (TTSs). Reduced nucleosome density was observed at all TTSs and TSSs for RNA genes transcribed by RNA polymerases I (RNAPI) or RNAPIII; as well as (to a lesser extent) at TSSs for the PTUs transcribed by RNAPII. Several PTMs (H3K4me3, H3K16ac H3K20me2 and H3K36me3) and base J were enriched at centromeres, while H3K50ac was specifically associated with the periphery of these centromeric sequences. These findings significantly expand our knowledge of the epigenetic markers associated with transcription, DNA replication and/or chromosome segregation in these early diverging eukaryotes and will hopefully lay the groundwork for future studies to elucidate how they control these fundamental processes.

Bromodomain factor 5 is an essential regulator of transcription in Leishmania

Leishmania are unicellular parasites that cause human and animal diseases. Like other kinetoplastids, they possess large transcriptional start regions (TSRs) which are defined by histone variants and histone lysine acetylation. Cellular interpretation of these chromatin marks is not well understood. Eight bromodomain factors, the reader modules for acetyl-lysine, are found across Leishmania genomes. Using L. mexicana, Cas9-driven gene deletions indicate that BDF1–5 are essential for promastigotes. Dimerisable, split Cre recombinase (DiCre)-inducible gene deletion of BDF5 show it is essential for both promastigotes and murine infection. ChIP-seq identifies BDF5 as enriched at TSRs. XL-BioID proximity proteomics shows the BDF5 landscape is enriched for BDFs, HAT2, proteins involved in transcriptional activity, and RNA processing; revealing a Conserved Regulators of Kinetoplastid Transcription (CRKT) Complex. Inducible deletion of BDF5 causes global reduction in RNA polymerase II transcription. Our results indicate the requirement of Leishmania to interpret histone acetylation marks through the bromodomain-enriched CRKT complex for normal gene expression and cellular viability.

Leishmania use large (5–10 kb) transcriptional start regions, where the chromatin is highly enriched for acetylated histones, to drive the expression of polycistronic gene arrays. Here the authors show bromodomain-containing protein BDF5 is enriched at transcriptional start sites and its depletion leads to cell death in vitro and in murine infections, and they identify its interactors.

Essential Bromodomain TcBDF2 as a Drug Target against Chagas Disease

Trypanosoma cruzi is a unicellular parasite that causes Chagas disease, which is endemic in the American continent but also worldwide, distributed by migratory movements. A striking feature of trypanosomatids is the polycistronic transcription associated with post-transcriptional mechanisms that regulate the levels of translatable mRNA. In this context, epigenetic regulatory mechanisms have been revealed to be of great importance, since they are the only ones that would control the access of RNA polymerases to chromatin. Bromodomains are epigenetic protein readers that recognize and specifically bind to acetylated lysine residues, mostly at histone proteins. There are seven coding sequences for BD-containing proteins in trypanosomatids, named TcBDF1 to TcBDF7, and a putative new protein containing a bromodomain was recently described. Using the Tet-regulated overexpression plasmid pTcINDEX-GW and CRISPR/Cas9 genome editing, we were able to demonstrate the essentiality of TcBDF2 in T. cruzi. This bromodomain is located in the nucleus, through a bipartite nuclear localization signal. TcBDF2 was shown to be important for host cell invasion, amastigote replication, and differentiation from amastigotes to trypomastigotes. Overexpression of TcBDF2 diminished epimastigote replication. Also, some processes involved in pathogenesis were altered in these parasites, such as infection of mammalian cells, replication of amastigotes, and the number of trypomastigotes released from host cells. In in vitro studies, TcBDF2 was also able to bind inhibitors showing a specificity profile different from that of the previously characterized TcBDF3. These results point to TcBDF2 as a druggable target against T. cruzi.

Impact of Genetic Diversity and Genome Plasticity of Leishmania spp. in Treatment and the Search for Novel Chemotherapeutic Targets

Leishmaniasis is one of the major public health concerns in Latin America, Africa, Asia, and Europe. The absence of vaccines for human use and the lack of effective vector control programs make chemotherapy the main strategy to control all forms of the disease. However, the high toxicity of available drugs, limited choice of therapeutic agents, and occurrence of drug-resistant parasite strains are the main challenges related to chemotherapy. Currently, only a small number of drugs are available for leishmaniasis treatment, including pentavalent antimonials (Sb^V), amphotericin B and its formulations, miltefosine, paromomycin sulphate, and pentamidine isethionate. In addition to drug toxicity, therapeutic failure of leishmaniasis is a serious concern. The occurrence of drug-resistant parasites is one of the causes of therapeutic failure and is closely related to the diversity of parasites in this genus. Owing to the enormous plasticity of the genome, resistance can occur by altering different metabolic pathways, demonstrating that resistance mechanisms are multifactorial and extremely complex. Genetic variability and genome plasticity cause not only the available drugs to have limitations, but also make the search for new drugs challenging. Here, we examined the biological characteristics of parasites that hinder drug discovery.

High genome plasticity and frequent genetic exchange in Leishmania tropica isolates from Afghanistan, Iran and Syria

Background

The kinetoplastid protozoan Leishmania tropica mainly causes cutaneous leishmaniasis in humans in the Middle East, and relapse or treatment failure after treatment are common in this area. L. tropica’s digenic life cycle includes distinct stages in the vector sandfly and the mammalian host. Sexual reproduction and genetic exchange appear to occur more frequently than in other Leishmania species. Understanding these processes is complicated by chromosome instability during cell division that yields aneuploidy, recombination and heterozygosity. This combination of rare recombination and aneuploid permits may reveal signs of hypothetical parasexual mating, where diploid cells fuse to form a transient tetraploid that undergoes chromosomal recombination and gradual chromosomal loss.

Methodology/principal findings

The genome-wide SNP diversity from 22 L. tropica isolates showed chromosome-specific runs of patchy heterozygosity and extensive chromosome copy number variation. All these isolates were collected during 2007–2017 in Sweden from patients infected in the Middle East and included isolates from a patient possessing two genetically distinct leishmaniasis infections three years apart with no evidence of re-infection. We found differing ancestries on the same chromosome (chr36) across multiple samples: matching the reference genome with few derived alleles, followed by blocks of heterozygous SNPs, and then by clusters of homozygous SNPs with specific recombination breakpoints at an inferred origin of replication. Other chromosomes had similar marked changes in heterozygosity at strand-switch regions separating polycistronic transcriptional units.

Conclusion/significance

These large-scale intra- and inter-chromosomal changes in diversity driven by recombination and aneuploidy suggest multiple mechanisms of cell reproduction and diversification in L. tropica, including mitotic, meiotic and parasexual processes. It underpins the need for more genomic surveillance of Leishmania, to detect emerging hybrids that could spread more widely and to better understand the association between genetic variation and treatment outcome. Furthering our understanding of Leishmania genome evolution and ancestry will aid better diagnostics and treatment for cutaneous leishmaniasis caused by L.tropica in the Middle East.

Cutaneous leishmaniasis is mainly caused by Leishmania tropica in the Middle East, where it is known for treatment failure and a need for prolonged and/or multiple treatments. Several factors affect the clinical presentation and treatment outcome, such as host genetic variability and specific immune response, as well as environmental factors and the vector species. Little is known about the parasite genome and its influence on treatment response. By analysing the genome of 22 isolates of L. tropica, we have revealed extensive genomic variation and a complex population structure with evidence of genetic exchange within and among the isolates, indicating a possible presence of sexual or parasexual mechanisms. Understanding the Leishmania genome better may improve future treatment and better understanding of treatment failure and relapse.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

The thermal proteome stability profile of Trypanosoma cruzi in epimastigote and trypomastigote life stages

Trypanosoma cruzi is a flagellate protozoa being the etiological agent of Chagas disease, a neglected tropical disease, which still poses a public health problem worldwide. The intricate molecular changes during T. cruzi-host interaction have been explored using different largescale omics techniques. However, protein stability is largely unknown. Thermal proteome profiling (TPP) methodology has the potential to characterize proteome-wide stability highlighting key proteins during T. cruzi infection and life stage transition from the invertebrate to the mammalian host. In the present work, T. cruzi epimastigotes and trypomastigotes cell lysates were subjected to TPP workflow and analyzed by quantitative large-scale mass spectrometry-based proteomics to fit a melting profile for each protein. A total of 2884 proteins were identified and associated to 1741 melting curves being 1370 in trypomastigotes (Tm_AVG 53.53 °C) and 1279 in epimastigotes (Tm_AVG 50.89 °C). A total of 453 proteins were identified with statistically different melting profiles between the two life stages. Proteins associated to pathogenesis and intracellular transport had regulated melting temperatures. Membrane and glycosylated proteins had a higher average Tm in trypomastigotes compared to epimastigotes. This study represents the first large-scale comparison of parasite protein stability between life stages. SIGNIFICANCE: Trypanosoma cruzi, a unicellular flagellate parasite, is the etiological agent of Chagas disease, endemic in South America and affecting more that 7 million people worldwide. There is an intense research to identify novel chemotherapeutic and diagnostic targets of Chagas disease. Proteomic approaches have helped in elucidating the quantitative proteome and PTMs changes of T. cruzi during life cycle transition and upon different biotic and abiotic stimuli. However, a comprehensive knowledge of the protein-protein interaction and protein conformation is still missing. In order to fill this gap, this manuscript elucidates the T. cruzi Y strain proteome-wide thermal stability map in the epimastigote and trypomastigote life stages. Comparison between life stages showed a higher average melting temperature stability for trypomastigotes than epimastigotes indicating a host temperature adaptation. Both presented a selective thermal stability shift for cellular compartments, molecular functions and biological processes based on the T. cruzi life stage. Membrane and glycosylated proteins presented a higher thermal stability in trypomastigotes when compared to the epimastigotes.

Life Cycle Stage-Specific Accessibility of Leishmania donovani Chromatin at Transcription Start Regions

Leishmania donovani is a parasitic protist that causes the lethal Kala-azar fever in India and East Africa. Gene expression in Leishmania is regulated by gene copy number variation and inducible translation while RNA synthesis initiates at a small number of sites per chromosome and proceeds through polycistronic transcription units, precluding a gene-specific regulation (C. Clayton and M. Shapira, Mol Biochem Parasitol 156:93–101, 2007, https://doi.org/10.1016/j.molbiopara.2007.07.007). Here, we analyze the dynamics of chromatin structure in both life cycle stages of the parasite and find evidence for an additional, epigenetic gene regulation pathway in this early branching eukaryote. The assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis (J. D. Buenrostro, P. G. Giresi, L. C. Zaba, H. Y. Chang, and W. J. Greenleaf, Nat Methods 10:1213–1218, 2013, https://doi.org/10.1038/nmeth.2688) predominantly shows euchromatin at transcription start regions in fast-growing promastigotes, but mostly heterochromatin in the slowly proliferating amastigotes, the mammalian stage, reflecting a previously shown increase of histone synthesis in the latter stage.

IMPORTANCELeishmania parasites are important pathogens with a global impact and cause poverty-related illness and death. They are devoid of classic cis- and trans-acting transcription regulators but use regulated translation and gene copy number variations to adapt to hosts and environments. In this work, we show that transcription start regions present as open euchromatin in fast-growing insect stages but as less-accessible heterochromatin in the slowly proliferating amastigote stage, indicating an epigenetic control of gene accessibility in this early branching eukaryotic pathogen. This finding should stimulate renewed interest in the control of RNA synthesis in Leishmania and related parasites.

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.

An updated map of Trypanosoma cruzi histone post-translational modifications

In humans and other eukaryotes, histone post-translational modifications (hPTMs) play an essential role in the epigenetic control of gene expression. In trypanosomatid parasites, conversely, gene regulation occurs mainly at the post-transcriptional level. However, our group has recently shown that hPTMs are abundant and varied in Trypanosoma cruzi, the etiological agent of Chagas Disease, signaling for possible conserved epigenetic functions. Here, we applied an optimized mass spectrometry-based proteomic workflow to provide a high-confidence comprehensive map of hPTMs, distributed in all canonical, variant and linker histones of T. cruzi. Our work expands the number of known T. cruzi hPTMs by almost 2-fold, representing the largest dataset of hPTMs available to any trypanosomatid to date, and can be used as a basis for functional studies on the dynamic regulation of chromatin by epigenetic mechanisms and the selection of candidates for the development of epigenetic drugs against trypanosomatids.

Measurement(s)	histone_modification
Technology Type(s)	mass spectrometry • nanoflow liquid chromatography-tandem mass spectrometry • Data-Dependent Acquisition
Sample Characteristic - Organism	Trypanosoma cruzi

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13491165

Chromatin-Associated Protein Complexes Link DNA Base J and Transcription Termination in Leishmania

Unlike most other eukaryotes, Leishmania and other trypanosomatid protozoa have largely eschewed transcriptional control of gene expression, relying instead on posttranscriptional regulation of mRNAs derived from polycistronic transcription units (PTUs). In these parasites, a novel modified nucleotide base (β-d-glucopyranosyloxymethyluracil) known as J plays a critical role in ensuring that transcription termination occurs only at the end of each PTU, rather than at the polyadenylation sites of individual genes. To further understand the biology of J-associated processes, we used tandem affinity purification (TAP) tagging and mass spectrometry to reveal proteins that interact with the glucosyltransferase performing the final step in J synthesis. These studies identified four proteins reminiscent of subunits in the PTW/PP1 complex that controls transcription termination in higher eukaryotes. Moreover, bioinformatic analyses identified the DNA-binding subunit of Leishmania PTW/PP1 as a novel J-binding protein (JBP3), which is also part of another complex containing proteins with domains suggestive of a role in chromatin modification/remodeling. Additionally, JBP3 associates (albeit transiently and/or indirectly) with the trypanosomatid equivalent of the PAF1 complex involved in the regulation of transcription in other eukaryotes. The downregulation of JBP3 expression levels in Leishmania resulted in a substantial increase in transcriptional readthrough at the 3′ end of most PTUs. We propose that JBP3 recruits one or more of these complexes to the J-containing regions at the end of PTUs, where they halt the progression of the RNA polymerase. This decoupling of transcription termination from the splicing of individual genes enables the parasites’ unique reliance on polycistronic transcription and posttranscriptional regulation of gene expression.

IMPORTANCELeishmania parasites cause a variety of serious human diseases, with no effective vaccine and emerging resistance to current drug therapy. We have previously shown that a novel DNA base called J is critical for transcription termination at the ends of the polycistronic gene clusters that are a hallmark of Leishmania and related trypanosomatids. Here, we describe a new J-binding protein (JBP3) associated with three different protein complexes that are reminiscent of those involved in the control of transcription in other eukaryotes. However, the parasite complexes have been reprogrammed to regulate transcription and gene expression in trypanosomatids differently than in the mammalian hosts, providing new opportunities to develop novel chemotherapeutic treatments against these important pathogens.

Participation of TFIIIB Subunit Brf1 in Transcription Regulation in the Human Pathogen Leishmania major

In yeast and higher eukaryotes, transcription factor TFIIIB is required for accurate initiation of transcription by RNA Polymerase III (Pol III), which synthesizes transfer RNAs (tRNAs), 5S ribosomal RNA (rRNA), and other essential RNA molecules. TFIIIB is composed of three subunits: B double prime 1 (Bdp1), TATA-binding protein (TBP), and TFIIB-related factor 1 (Brf1). Here, we report the molecular characterization of Brf1 in Leishmania major (LmBrf1), a parasitic protozoan that shows distinctive transcription characteristics, including the apparent absence of Pol III general transcription factors TFIIIA and TFIIIC. Although single-knockout parasites of LmBrf1 were obtained, attempts to generate LmBrf1-null mutants were unsuccessful, which suggests that LmBrf1 is essential in promastigotes of L. major. Notably, Northern blot analyses showed that the half-lives of the messenger RNAs (mRNAs) from LmBrf1 and other components of the Pol III transcription machinery (Bdp1 and Pol III subunit RPC1) are very similar (~40 min). Stabilization of these transcripts was observed in stationary-phase parasites. Chromatin immunoprecipitation (ChIP) experiments showed that LmBrf1 binds to tRNA, small nuclear RNA (snRNA), and 5S rRNA genes. Unexpectedly, the results also indicated that LmBrf1 associates to the promoter region of the 18S rRNA genes and to three Pol II-dependent regions here analyzed. Tandem affinity purification and mass spectrometry analyses allowed the identification of a putative TFIIIC subunit. Moreover, several proteins involved in transcription by all three RNA polymerases co-purified with the tagged version of LmBrf1.

Nuclear mRNA maturation and mRNA export control: from trypanosomes to opisthokonts

The passage of mRNAs through the nuclear pores into the cytoplasm is essential in all eukaryotes. For regulation, mRNA export is tightly connected to the full machinery of nuclear mRNA processing, starting at transcription. Export competence of pre-mRNAs gradually increases by both transient and permanent interactions with multiple RNA processing and export factors. mRNA export is best understood in opisthokonts, with limited knowledge in plants and protozoa. Here, I review and compare nuclear mRNA processing and export between opisthokonts and Trypanosoma brucei. The parasite has many unusual features in nuclear mRNA processing, such as polycistronic transcription and trans-splicing. It lacks several nuclear complexes and nuclear-pore-associated proteins that in opisthokonts play major roles in mRNA export. As a consequence, trypanosome mRNA export control is not tight and export can even start co-transcriptionally. Whether trypanosomes regulate mRNA export at all, or whether leakage of immature mRNA to the cytoplasm is kept to a low level by a fast kinetics of mRNA processing remains to be investigated. mRNA export had to be present in the last common ancestor of eukaryotes. Trypanosomes are evolutionary very distant from opisthokonts and a comparison helps understanding the evolution of mRNA export.

Nucleosome landscape reflects phenotypic differences in Trypanosoma cruzi life forms

Trypanosoma cruzi alternates between replicative and nonreplicative life forms, accompanied by a shift in global transcription levels and by changes in the nuclear architecture, the chromatin proteome and histone posttranslational modifications. To gain further insights into the epigenetic regulation that accompanies life form changes, we performed genome-wide high-resolution nucleosome mapping using two T. cruzi life forms (epimastigotes and cellular trypomastigotes). By combining a powerful pipeline that allowed us to faithfully compare nucleosome positioning and occupancy, more than 125 thousand nucleosomes were mapped, and approximately 20% of them differed between replicative and nonreplicative forms. The nonreplicative forms have less dynamic nucleosomes, possibly reflecting their lower global transcription levels and DNA replication arrest. However, dynamic nucleosomes are enriched at nonreplicative regulatory transcription initiation regions and at multigenic family members, which are associated with infective-stage and virulence factors. Strikingly, dynamic nucleosome regions are associated with GO terms related to nuclear division, translation, gene regulation and metabolism and, notably, associated with transcripts with different expression levels among life forms. Finally, the nucleosome landscape reflects the steady-state transcription expression: more abundant genes have a more deeply nucleosome-depleted region at putative 5' splice sites, likely associated with trans-splicing efficiency. Taken together, our results indicate that chromatin architecture, defined primarily by nucleosome positioning and occupancy, reflects the phenotypic differences found among T. cruzi life forms despite the lack of a canonical transcriptional control context.

Trypanosoma Cruzi Genome: Organization, Multi-Gene Families, Transcription, and Biological Implications

Chagas disease caused by the parasite Trypanosoma cruzi affects millions of people. Although its first genome dates from 2005, its complexity hindered a complete assembly and annotation. However, the new sequencing methods have improved genome annotation of some strains elucidating the broad genetic diversity and complexity of this parasite. Here, we reviewed the genomic structure and regulation, the genetic diversity, and the analysis of the principal multi-gene families of the recent genomes for several strains. The telomeric and sub-telomeric regions are sites with high recombination events, the genome displays two different compartments, the core and the disruptive, and the genome plasticity seems to play a key role in the survival and the infection process. Trypanosoma cruzi (T. cruzi) genome is composed mainly of multi-gene families as the trans-sialidases, mucins, and mucin-associated surface proteins. Trans-sialidases are the most abundant genes in the genome and show an important role in the effectiveness of the infection and the parasite survival. Mucins and MASPs are also important glycosylated proteins of the surface of the parasite that play a major biological role in both insect and mammal-dwelling stages. Altogether, these studies confirm the complexity of T. cruzi genome revealing relevant concepts to better understand Chagas disease.

Genome duplication in Leishmania major relies on persistent subtelomeric DNA replication

DNA replication is needed to duplicate a cell’s genome in S phase and segregate it during cell division. Previous work in Leishmania detected DNA replication initiation at just a single region in each chromosome, an organisation predicted to be insufficient for complete genome duplication within S phase. Here, we show that acetylated histone H3 (AcH3), base J and a kinetochore factor co-localise in each chromosome at only a single locus, which corresponds with previously mapped DNA replication initiation regions and is demarcated by localised G/T skew and G4 patterns. In addition, we describe previously undetected subtelomeric DNA replication in G2/M and G1-phase-enriched cells. Finally, we show that subtelomeric DNA replication, unlike chromosome-internal DNA replication, is sensitive to hydroxyurea and dependent on 9-1-1 activity. These findings indicate that Leishmania’s genome duplication programme employs subtelomeric DNA replication initiation, possibly extending beyond S phase, to support predominantly chromosome-internal DNA replication initiation within S phase.

The Experimental Proteome of Leishmania infantum Promastigote and Its Usefulness for Improving Gene Annotations

Leishmania infantum causes visceral leishmaniasis (kala-azar), the most severe form of leishmaniasis, which is lethal if untreated. A few years ago, the re-sequencing and de novo assembling of the L. infantum (JPCM5 strain) genome was accomplished, and now we aimed to describe and characterize the experimental proteome of this species. In this work, we performed a proteomic analysis from axenic cultured promastigotes and carried out a detailed comparison with other Leishmania experimental proteomes published to date. We identified 2352 proteins based on a search of mass spectrometry data against a database built from the six-frame translated genome sequence of L. infantum. We detected many proteins belonging to organelles such as glycosomes, mitochondria, or flagellum, as well as many metabolic enzymes and many putative RNA binding proteins and molecular chaperones. Moreover, we listed some proteins presenting post-translational modifications, such as phosphorylations, acetylations, and methylations. On the other hand, the identification of peptides mapping to genomic regions previously annotated as non-coding allowed for the correction of annotations, leading to the N-terminal extension of protein sequences and the uncovering of eight novel protein-coding genes. The alliance of proteomics, genomics, and transcriptomics has resulted in a powerful combination for improving the annotation of the L. infantum reference genome.

Histone Modifications and Other Facets of Epigenetic Regulation in Trypanosomatids: Leaving Their Mark

Histone posttranslational modifications (PTMs) modulate several eukaryotic cellular processes, including transcription, replication, and repair. Vast arrays of modifications have been identified in conventional eukaryotes over the last 20 to 25 years. While initial studies uncovered these primarily on histone tails, multiple modifications were subsequently found on the central globular domains as well. Histones are evolutionarily conserved across eukaryotes, and a large number of their PTMs and the functional relevance of these PTMs are largely conserved.

Histone posttranslational modifications (PTMs) modulate several eukaryotic cellular processes, including transcription, replication, and repair. Vast arrays of modifications have been identified in conventional eukaryotes over the last 20 to 25 years. While initial studies uncovered these primarily on histone tails, multiple modifications were subsequently found on the central globular domains as well. Histones are evolutionarily conserved across eukaryotes, and a large number of their PTMs and the functional relevance of these PTMs are largely conserved. Trypanosomatids, however, are early diverging eukaryotes. Although possessing all four canonical histones as well as several variants, their sequences diverge from those of other eukaryotes, particularly in the tails. Consequently, the modifications they carry also vary. Initial analyses almost 15 years ago suggested that trypanosomatids possessed a smaller collection of histone modifications. However, exhaustive high resolution mass spectrometry analyses in the last few years have overturned this belief, and it is now evident that the “histone code” proposed by Allis and coworkers in the early years of this century is as complex in these organisms as in other eukaryotes. Trypanosomatids cause several diseases, and the members of this group of organisms have varied lifestyles, evolving diverse mechanisms to evade the host immune system, some of which have been found to be principally controlled by epigenetic mechanisms. This minireview aims to acquaint the reader with the impact of histone PTMs on trypanosomatid cellular processes, as well as other facets of trypanosomatid epigenetic regulation, including the influence of three-dimensional (3D) genome architecture, and discusses avenues for future investigations.

The highly diverse TATA box-binding proteins among protists: A review

Transcription is the first step of gene expression regulation and is a fundamental mechanism for establishing the viability and development of a cell. The TATA box-binding protein (TBP) interaction with a TATA box in a promoter is one of the best studied mechanisms in transcription initiation. TBP is a transcription factor that is highly conserved from archaea to humans and is essential for the transcription initiated by each of the three RNA polymerases. In addition, the discovery of TBP-related factor 1 (TRF1) and other factors related to TBP shed light on the variability among transcription initiation complexes, thus demonstrating that the compositions of these complexes are, in fact, more complicated than originally believed. Despite these facts, the majority of studies on transcription have been performed on animal, plant and fungal cells, which serve as canonical models, and information regarding protist cells is relatively scarce. The aim of this work is to review the diversity of the TBPs that have been documented in protists and describe some of the specific features that differentiate them from their counterparts in higher eukaryotes.

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.

An RNA polymerase II-associated TFIIF-like complex is indispensable for SL RNA gene transcription in Trypanosoma brucei

Trypanosomes are protistan parasites that diverged early in evolution from most eukaryotes. Their streamlined genomes are packed with arrays of tandemly linked genes that are transcribed polycistronically by RNA polymerase (pol) II. Individual mRNAs are processed from pre-mRNA by spliced leader (SL) trans splicing and polyadenylation. While there is no strong evidence that general transcription factors are needed for transcription initiation at these gene arrays, a RNA pol II transcription pre-initiation complex (PIC) is formed on promoters of SLRNA genes, which encode the small nuclear SL RNA, the SL donor in trans splicing. The factors that form the PIC are extremely divergent orthologues of the small nuclear RNA-activating complex, TBP, TFIIA, TFIIB, TFIIH, TFIIE and Mediator. Here, we functionally characterized a heterodimeric complex of unannotated, nuclear proteins that interacts with RNA pol II and is essential for PIC formation, SL RNA synthesis in vivo, SLRNA transcription in vitro, and parasite viability. These functional attributes suggest that the factor represents TFIIF although the amino acid sequences are too divergent to firmly make this conclusion. This work strongly indicates that early-diverged trypanosomes have orthologues of each and every general transcription factor, requiring them for the synthesis of SL RNA.

Nucleolar Structure and Function in Trypanosomatid Protozoa

The nucleolus is the conspicuous nuclear body where ribosomal RNA genes are transcribed by RNA polymerase I, pre-ribosomal RNA is processed, and ribosomal subunits are assembled. Other important functions have been attributed to the nucleolus over the years. Here we review the current knowledge about the structure and function of the nucleolus in the trypanosomatid parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania ssp., which represent one of the earliest branching lineages among the eukaryotes. These protozoan parasites present a single nucleolus that is preserved throughout the closed nuclear division, and that seems to lack fibrillar centers. Trypanosomatids possess a relatively low number of rRNA genes, which encode rRNA molecules that contain large expansion segments, including several that are trypanosomatid-specific. Notably, the large subunit rRNA (28S-type) is fragmented into two large and four small rRNA species. Hence, compared to other organisms, the rRNA primary transcript requires additional processing steps in trypanosomatids. Accordingly, this group of parasites contains the highest number ever reported of snoRNAs that participate in rRNA processing. The number of modified rRNA nucleotides in trypanosomatids is also higher than in other organisms. Regarding the structure and biogenesis of the ribosomes, recent cryo-electron microscopy analyses have revealed several trypanosomatid-specific features that are discussed here. Additional functions of the nucleolus in trypanosomatids are also reviewed.

Analysis by RNA-seq of transcriptomic changes elicited by heat shock in Leishmania major

Besides their medical relevance, Leishmania is an adequate model for studying post-transcriptional mechanisms of gene expression. In this microorganism, mRNA degradation/stabilization mechanisms together with translational control and post-translational modifications of proteins are the major drivers of gene expression. Leishmania parasites develop as promastigotes in sandflies and as amastigotes in mammalians, and during host transmission, the parasite experiences a sudden temperature increase. Here, changes in the transcriptome of Leishmania major promastigotes after a moderate heat shock were analysed by RNA-seq. Several of the up-regulated transcripts code for heat shock proteins, other for proteins previously reported to be amastigote-specific and many for hypothetical proteins. Many of the transcripts experiencing a decrease in their steady-state levels code for transporters, proteins involved in RNA metabolism or translational factors. In addition, putative long noncoding RNAs were identified among the differentially expressed transcripts. Finally, temperature-dependent changes in the selection of the spliced leader addition sites were inferred from the RNA-seq data, and particular cases were further validated by RT-PCR and Northern blotting. This study provides new insights into the post-transcriptional mechanisms by which Leishmania modulate gene expression.

Leishmania-Host Interactions-An Epigenetic Paradigm

Leishmaniasis is one of the major neglected tropical diseases, for which no vaccines exist. Chemotherapy is hampered by limited efficacy coupled with development of resistance and other side effects. Leishmania parasites elude the host defensive mechanisms by modulating their surface proteins as well as dampening the host's immune responses. The parasites use the conventional RNA polymerases peculiarly under different environmental cues or pressures such as the host's milieu or the drugs. The mechanisms that restructure post-translational modifications are poorly understood but altered epigenetic histone modifications are believed to be instrumental in influencing the chromatin remodeling in the parasite. Interestingly, the parasite also modulates gene expression of the hosts, thereby hijacking or dampening the host immune response. Epigenetic factor such as DNA methylation of cytosine residues has been incriminated in silencing of macrophage-specific genes responsible for defense against these parasites. Although there is dearth of information regarding the epigenetic alterations-mediated pathogenesis in these parasites and the host, the unique epigenetic marks may represent targets for potential anti-leishmanial drug candidates. This review circumscribes the epigenetic changes during Leishmania infection, and the epigenetic modifications they enforce upon the host cells to ensure a safe haven. The non-coding micro RNAs as post-transcriptional regulators and correlates of wound healing and toll-like receptor signaling, as well as prognostic biomarkers of therapeutic failure and healing time are also explored. Finally, we highlight the recent advances on how the epigenetic perturbations may impact leishmaniasis vaccine development as biomarkers of safety and immunogenicity.

UTRme: A Scoring-Based Tool to Annotate Untranslated Regions in Trypanosomatid Genomes

Most signals involved in post-transcriptional regulatory networks are located in the untranslated regions (UTRs) of the mRNAs. Therefore, to deepen our understanding of gene expression regulation, delimitation of these regions with high accuracy is needed. The trypanosomatid lineage includes a variety of parasitic protozoans causing a significant worldwide burden on human health. Given their peculiar mechanisms of gene expression, these organisms depend on post-transcriptional regulation as the main level of gene expression control. In this context, the definition of the UTR regions becomes of key importance. We have developed UTR-mini-exon (UTRme), a graphical user interface (GUI) stand-alone application to identify and annotate 5′ and 3′ UTR regions in a highly accurate way. UTRme implements a multiple scoring system tailored to address the issue of false positive UTR assignment that frequently arise because of the characteristics of the intergenic regions. Even though it was developed for trypanosomatids, the tool can be used to predict 3′ sites in any eukaryote and 5′ UTRs in any organism where trans-splicing occurs (such as the model organism C. elegans). UTRme offers a way for non-bioinformaticians to precisely determine UTRs from transcriptomic data. The tool is freely available via the conda and github repositories.

ChIP-ping the branches of the tree: functional genomics and the evolution of eukaryotic gene regulation

Advances in the methods for detecting protein-DNA interactions have played a key role in determining the directions of research into the mechanisms of transcriptional regulation. The most recent major technological transformation happened a decade ago, with the move from using tiling arrays [chromatin immunoprecipitation (ChIP)-on-Chip] to high-throughput sequencing (ChIP-seq) as a readout for ChIP assays. In addition to the numerous other ways in which it is superior to arrays, by eliminating the need to design and manufacture them, sequencing also opened the door to carrying out comparative analyses of genome-wide transcription factor occupancy across species and studying chromatin biology in previously less accessible model and nonmodel organisms, thus allowing us to understand the evolution and diversity of regulatory mechanisms in unprecedented detail. Here, we review the biological insights obtained from such studies in recent years and discuss anticipated future developments in the field.

Epigenetic Regulation of Transcription in Trypanosomatid Protozoa

The Trypanosomatid family includes flagellated parasites that cause fatal human diseases. Remarkably, protein-coding genes in these organisms are positioned in long tandem arrays that are transcribed polycistronically. However, the knowledge about regulation of transcription initiation and termination in trypanosomatids is scarce. The importance of epigenetic regulation in these processes has become evident in the last years, as distinctive histone modifications and histone variants have been found in transcription initiation and termination regions. Moreover, multiple chromatin-related proteins have been identified and characterized in trypanosomatids, including histone-modifying enzymes, effector complexes, chromatin-remodelling enzymes and histone chaperones. Notably, base J, a modified thymine residue present in the nuclear DNA of trypanosomatids, has been implicated in transcriptional regulation. Here we review the current knowledge on epigenetic control of transcription by all three RNA polymerases in this group of early-diverged eukaryotes.

Intrinsic DNA curvature in trypanosomes

Background

Trypanosoma cruzi and Trypanosoma brucei are protozoan parasites causing Chagas disease and African sleeping sickness, displaying unique features of cellular and molecular biology. Remarkably, no canonical signals for RNA polymerase II promoters, which drive protein coding genes transcription, have been identified so far. The secondary structure of DNA has long been recognized as a signal in biological processes and more recently, its involvement in transcription initiation in Leishmania was proposed. In order to study whether this feature is conserved in trypanosomatids, we undertook a genome wide search for intrinsic DNA curvature in T. cruzi and T. brucei.

Results

Using a region integrated intrinsic curvature (RIIC) scoring that we previously developed, a non-random distribution of sequence-dependent curvature was observed. High RIIC scores were found to be significantly correlated with transcription start sites in T. cruzi, which have been mapped in divergent switch regions, whereas in T. brucei, the high RIIC scores correlated with sites that have been involved not only in RNA polymerase II initiation but also in termination. In addition, we observed regions with high RIIC score presenting in-phase tracts of Adenines, in the subtelomeric regions of the T. brucei chromosomes that harbor the variable surface glycoproteins genes.

Conclusions

In both T. cruzi and T. brucei genomes, a link between DNA conformational signals and gene expression was found. High sequence dependent curvature is associated with transcriptional regulation regions. High intrinsic curvature also occurs at the T. brucei chromosome subtelomeric regions where the recombination processes involved in the evasion of the immune host system take place. These findings underscore the relevance of indirect DNA readout in these ancient eukaryotes.

Electronic supplementary material

The online version of this article (10.1186/s13104-017-2908-y) contains supplementary material, which is available to authorized users.

An essential domain of an early-diverged RNA polymerase II functions to accurately decode a primitive chromatin landscape

A unique feature of RNA polymerase II (RNA pol II) is its long C-terminal extension, called the carboxy-terminal domain (CTD). The well-studied eukaryotes possess a tandemly repeated 7-amino-acid sequence, called the canonical CTD, which orchestrates various steps in mRNA synthesis. Many eukaryotes possess a CTD devoid of repeats, appropriately called a non-canonical CTD, which performs completely unknown functions. Trypanosoma brucei, the etiologic agent of African Sleeping Sickness, deploys an RNA pol II that contains a non-canonical CTD to accomplish an unusual transcriptional program; all protein-coding genes are transcribed as part of a polygenic precursor mRNA (pre-mRNA) that is initiated within a several-kilobase-long region, called the transcription start site (TSS), which is upstream of the first protein-coding gene in the polygenic array. In this report, we show that the non-canonical CTD of T. brucei RNA pol II is important for normal protein-coding gene expression, likely directing RNA pol II to the TSSs within the genome. Our work reveals the presence of a primordial CTD code within eukarya and indicates that proper recognition of the chromatin landscape is a central function of this RNA pol II-distinguishing domain.

Evidence for regulated expression of Telomeric Repeat-containing RNAs (TERRA) in parasitic trypanosomatids

The Telomeric Repeat-containing RNAs (TERRA) participate in the homeostasis of telomeres in higher eukaryotes. Here, we investigated the expression of TERRA in Leishmania spp. and Trypanosoma brucei and found evidences for its expression as a specific RNA class. The trypanosomatid TERRA are heterogeneous in size and partially polyadenylated. The levels of TERRA transcripts appear to be modulated through the life cycle in both trypanosomatids investigated, suggesting that TERRA play a stage-specific role in the life cycle of these early-branching eukaryotes.

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.

Identification of the centromeres of Leishmania major: revealing the hidden pieces

Leishmania affects millions of people worldwide. Its genome undergoes constitutive mosaic aneuploidy, a type of genomic plasticity that may serve as an adaptive strategy to survive distinct host environments. We previously found high rates of asymmetric chromosome allotments during mitosis that lead to the generation of such ploidy. However, the underlying molecular events remain elusive. Centromeres and kinetochores most likely play a key role in this process, yet their identification has failed using classical methods. Our analysis of the unconventional kinetochore complex recently discovered in Trypanosoma brucei (KKTs) leads to the identification of a Leishmania KKT gene candidate (LmKKT1). The GFP-tagged LmKKT1 displays "kinetochore-like" dynamics of intranuclear localization throughout the cell cycle. By ChIP-Seq assay, one major peak per chromosome is revealed, covering a region of 4 ±2 kb. We find two largely conserved motifs mapping to 14 of 36 chromosomes while a higher density of retroposons are observed in 27 of 36 centromeres. The identification of centromeres and of a kinetochore component of Leishmania chromosomes opens avenues to explore their role in mosaic aneuploidy.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

HAT2 mediates histone H4K4 acetylation and affects micrococcal nuclease sensitivity of chromatin in Leishmania donovani

Histone post-translational modifications (PTMs) such as acetylation and methylation are known to affect chromatin higher order structures. Primary targets of these modifications include basic residues present at N-terminus tail region of core histones. Four histone acetyltransferase (HAT) genes have been identified in trypanosomatids. HAT1, HAT3 and HAT4 of Leishmania donovani have been partially characterized. However, there is no report about HAT2 of Leishmania donovani. Lysine residues present on the N-terminal tail of Leishmania donovani histone H4 are conserved in other trypanosomatids and humans. PTMs of lysines modulate various functions at chromatin level. The four histone acetyltransferases encoded in Leishmania genome were over-expressed to analyse their functional activity. All four HATs were found actively acetylating core histones H3/H4. Similar to L. donovani HAT3 and HAT4, HAT2 was found to be a member of MYST family protein and have SAS2 type domain. Over-expression of HAT2 significantly increases acetylation of H4K4. To analyse the effect of HAT2 over-expression on chromatin accessibility, micrococcal nuclease digestion assay was performed. MNase digestion resulted in a higher proportion of the mononucleosomes and dinucleosomes in HAT2 over-expressing cells as compared to WT L. donovani cells. Acetylation of lysine-4 neutralizes the amino terminal region of histone H4. This weakens its interaction with neighbouring nucleosomes and the linker DNA. HAT2 over-expression in L. donovani resulted in highly accessible chromatin suggesting chromatin decondensation. HAT2 may have an important role to play in global regulation of transcription in L. donovani. Better understanding of these epigenetic determinants of parasite would help in designing novel therapeutic strategies.

The genome of Leishmania adleri from a mammalian host highlights chromosome fission in Sauroleishmania

Control of pathogens arising from humans, livestock and wild animals can be enhanced by genome-based investigation. Phylogenetically classifying and optimal construction of these genomes using short sequence reads are key to this process. We examined the mammal-infecting unicellular parasite Leishmania adleri belonging to the lizard-infecting Sauroleishmania subgenus. L. adleri has been associated with cutaneous disease in humans, but can be asymptomatic in wild animals. We sequenced, assembled and investigated the L. adleri genome isolated from an asymptomatic Ethiopian rodent (MARV/ET/75/HO174) and verified it as L. adleri by comparison with other Sauroleishmania species. Chromosome-level scaffolding was achieved by combining reference-guided with de novo assembly followed by extensive improvement steps to produce a final draft genome with contiguity comparable with other references. L. tarentolae and L. major genome annotation was transferred and these gene models were manually verified and improved. This first high-quality draft Leishmania adleri reference genome is also the first Sauroleishmania genome from a non-reptilian host. Comparison of the L. adleri HO174 genome with those of L. tarentolae Parrot-TarII and lizard-infecting L. adleri RLAT/KE/1957/SKINK-7 showed extensive gene amplifications, pervasive aneuploidy, and fission of chromosomes 30 and 36. There was little genetic differentiation between L. adleri extracted from mammals and reptiles, highlighting challenges for leishmaniasis surveillance.

Genome-wide mapping of 5-hydroxymethyluracil in the eukaryote parasite Leishmania

BACKGROUND

5-Hydroxymethyluracil (5hmU) is a thymine base modification found in the genomes of a diverse range of organisms. To explore the functional importance of 5hmU, we develop a method for the genome-wide mapping of 5hmU-modified loci based on a chemical tagging strategy for the hydroxymethyl group.

RESULTS

We apply the method to generate genome-wide maps of 5hmU in the parasitic protozoan Leishmania sp. In this genus, another thymine modification, 5-(β-glucopyranosyl) hydroxymethyluracil (base J), plays a key role during transcription. To elucidate the relationship between 5hmU and base J, we also map base J loci by introducing a chemical tagging strategy for the glucopyranoside residue. Observed 5hmU peaks are highly consistent among technical replicates, confirming the robustness of the method. 5hmU is enriched in strand switch regions, telomeric regions, and intergenic regions. Over 90% of 5hmU-enriched loci overlapped with base J-enriched loci, which occurs mostly within strand switch regions. We also identify loci comprising 5hmU but not base J, which are enriched with motifs consisting of a stretch of thymine bases.

CONCLUSIONS

By chemically detecting 5hmU we present a method to provide a genome-wide map of this modification, which will help address the emerging interest in the role of 5hmU. This method will also be applicable to other organisms bearing 5hmU.