Two thymidine hydroxylases differentially regulate the formation of glucosylated DNA at regions flanking polymerase II polycistronic transcription units throughout the genome of Trypanosoma brucei

JBP1 and JBP3 have conserved structures but different affinity to base-J

Base-J (β-D-glucopyranosyloxymethyluracil) is an unusual kinetoplastid-specific DNA modification, recognized by base-J containing DNA (J-DNA) binding proteins JBP1 and JBP3. Recognition of J-DNA by both JBP1 and JBP3 takes place by a conserved J-DNA binding domain (JDBD). Here we show that JDBD-JBP3 has about 1,000-fold weaker affinity to base-J than JDBD-JBP1 and discriminates between J-DNA and unmodified DNA with a factor ∼5, whereas JDBD-JBP1 discriminates with a factor ∼10,000. Comparison of the crystal structures of JDBD-JBP3 we present here, with that of the previously characterized JDBD-JBP1, shows a flexible α5-helix that lacks a positively charged patch in JBP3. Mutations removing this positive charge in JDBD-JBP1, resulted in decreased binding affinity relative to wild-type JDBD-JBP1, indicating this patch is involved in DNA binding. We suggest that the α5-helix might rearrange upon JBP1 binding to J-DNA stabilizing the complex. This work contributes to our understanding of how JBPs bind to this unique DNA modification, which may contribute to identifying potential drug targets to end the base-J dependent parasite life cycle.

Mono-allelic epigenetic regulation of polycistronic transcription initiation by RNA polymerase II in Trypanosoma brucei

Unique for a eukaryote, protein-coding genes in trypanosomes are arranged in polycistronic transcription units (PTUs). This genome arrangement has led to a model where Pol II transcription of PTUs is unregulated and 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. Base J function in establishing transcriptional silencing, rather than maintenance, is suggested by the maintenance of PTU silencing following the inhibition of J-biosynthesis and subsequent loss of the modified DNA base. Therefore, we show that epigenetic mechanisms exist to regulate 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.IMPORTANCEThe 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 correlate 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.

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.

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.

Epigenetic regulation of TERRA transcription and metacyclogenesis by base J in Leishmania major

The hyper-modified DNA base J helps control termination of Pol II transcription at polycistronic transcription units (PTUs) in T. brucei and L. major , allowing epigenetic control of gene expression. The Telomere Repeat-containing RNA (TERRA) is synthesized in T. brucei by Pol I readthrough transcription of a telomeric PTU. While little is understood regarding TERRA synthesis and function, the hyper-modified DNA base J is highly enriched at telomeres in L. major promastigotes. We now show that TERRA is synthesized by Pol II in L. major and loss of base J leads to increased TERRA. For at least one site, the increased TERRA is by Pol II readthrough transcription from an adjacent PTU. Furthermore, Pol II readthrough defects and increased TERRA correlate with increased differentiation of promastigotes to the infectious metacyclic life stage and decreased cell viability. These results help explain the essential nature of base J in Leishmania and provide insight regarding epigenetic control of coding and non-coding RNA expression and parasite development during the life cycle of L. major .

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

BACKGROUND

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties are originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools.

RESULTS

In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

CONCLUSIONS

Our results provide the first window into the 5-hmU and chromatin accessibility landscapes in dinoflagellates.

Epigenetic marks or not? The discovery of novel DNA modifications in eukaryotes

DNA modifications add another layer of complexity to the eukaryotic genome to regulate gene expression, playing critical roles as epigenetic marks. In eukaryotes, the study of DNA epigenetic modifications has been confined to 5mC and its derivatives for decades. However, rapid developing approaches have witnessed the expansion of DNA modification reservoirs during the past several years, including the identification of 6mA, 5gmC, 4mC, and 4acC in diverse organisms. However, whether these DNA modifications function as epigenetic marks requires careful consideration. In this review, we try to present a panorama of all the DNA epigenetic modifications in eukaryotes, emphasizing recent breakthroughs in the identification of novel DNA modifications. The characterization of their roles in transcriptional regulation as potential epigenetic marks is summarized. More importantly, the pathways for generating or eliminating these DNA modifications, as well as the proteins involved are comprehensively dissected. Furthermore, we briefly discuss the potential challenges and perspectives, which should be taken into account while investigating novel DNA modifications.

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.

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.

Genome-wide distribution of 5-hydroxymethyluracil and chromatin accessibility in the Breviolum minutum genome

In dinoflagellates, a unique and extremely divergent genomic and nuclear organization has evolved. The highly unusual features of dinoflagellate nuclei and genomes include permanently condensed liquid crystalline chromosomes, primarily packaged by proteins other than histones, genes organized in very long unidirectional gene arrays, a general absence of transcriptional regulation, high abundance of the otherwise very rare DNA modification 5-hydroxymethyluracil (5-hmU), and many others. While most of these fascinating properties were originally identified in the 1970s and 1980s, they have not yet been investigated using modern genomic tools. In this work, we address some of the outstanding questions regarding dinoflagellate genome organization by mapping the genome-wide distribution of 5-hmU (using both immunoprecipitation-based and basepair-resolution chemical mapping approaches) and of chromatin accessibility in the genome of the Symbiodiniaceae dinoflagellate Breviolum minutum. We find that the 5-hmU modification is preferentially enriched over certain classes of repetitive elements, often coincides with the boundaries between gene arrays, and is generally correlated with decreased chromatin accessibility, the latter otherwise being largely uniform along the genome. We discuss the potential roles of 5-hmU in the functional organization of dinoflagellate genomes and its relationship to the transcriptional landscape of gene arrays.

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

PP1 phosphatases lack substrate specificity and associate with specific regulatory subunits to achieve 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 non-canonical 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-termini for independent binding of JBP3 and Wdr82, 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.

Distant sequence regions of JBP1 contribute to J-DNA binding

Base-J (β-D-glucopyranosyloxymethyluracil) is a modified DNA nucleotide that replaces 1% of thymine in kinetoplastid flagellates. The biosynthesis and maintenance of base-J depends on the base-J-binding protein 1 (JBP1) that has a thymidine hydroxylase domain and a J-DNA-binding domain (JDBD). How the thymidine hydroxylase domain synergizes with the JDBD to hydroxylate thymine in specific genomic sites, maintaining base-J during semi-conservative DNA replication, remains unclear. Here, we present a crystal structure of the JDBD including a previously disordered DNA-contacting loop and use it as starting point for molecular dynamics simulations and computational docking studies to propose recognition models for JDBD binding to J-DNA. These models guided mutagenesis experiments, providing additional data for docking, which reveals a binding mode for JDBD onto J-DNA. This model, together with the crystallographic structure of the TET2 JBP1-homologue in complex with DNA and the AlphaFold model of full-length JBP1, allowed us to hypothesize that the flexible JBP1 N-terminus contributes to DNA-binding, which we confirmed experimentally. Α high-resolution JBP1:J-DNA complex, which must involve conformational changes, would however need to be determined experimentally to further understand this unique underlying molecular mechanism that ensures replication of epigenetic information.

Navigating the boundaries between metabolism and epigenetics in trypanosomes

Epigenetic marks enable cells to acquire new biological features that favor their adaptation to environmental changes. These marks are chemical modifications on chromatin-associated proteins and nucleic acids that lead to changes in the chromatin landscape and may eventually affect gene expression. The chemical tags of these epigenetic marks are comprised of intermediate cellular metabolites. The number of discovered associations between metabolism and epigenetics has increased, revealing how environment influences gene regulation and phenotype diversity. This connection is relevant to all organisms but underappreciated in digenetic parasites, which must adapt to different environments as they progress through their life cycles. This review speculates and proposes associations between epigenetics and metabolism in trypanosomes, which are protozoan parasites that cause human and livestock diseases.

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.

G-Quadruplexes as Key Transcriptional Regulators in Neglected Trypanosomatid Parasites

G-quadruplexes (G4s) are nucleic acid secondary structures that have been linked to the functional regulation of eukaryotic organisms. G4s have been extensively characterised in humans and emerging evidence suggests that they might also be biologically relevant for human pathogens. This indicates that G4s might represent a novel class of therapeutic targets for tackling infectious diseases. Bioinformatic studies revealed a high prevalence of putative quadruplex-forming sequences (PQSs) in the genome of protozoans, which highlights their potential roles in regulating vital processes of these parasites, including DNA transcription and replication. In this work, we focus on the neglected trypanosomatid parasites, Trypanosoma and Leishmania spp., which cause debilitating and deadly diseases across the poorest populations worldwide. We review three examples where G4-formation might be key to modulate transcriptional activity in trypanosomatids, providing an overview of experimental approaches that can be used to exploit the regulatory roles and relevance of these structures to fight parasitic infections.

Kinetic Studies on the 2-Oxoglutarate/Fe(II)-Dependent Nucleic Acid Modifying Enzymes from the AlkB and TET Families

Nucleic acid methylations are important genetic and epigenetic biomarkers. The formation and removal of these markers is related to either methylation or demethylation. In this review, we focus on the demethylation or oxidative modification that is mediated by the 2-oxoglutarate (2-OG)/Fe(II)-dependent AlkB/TET family enzymes. In the catalytic process, most enzymes oxidize 2-OG to succinate, in the meantime oxidizing methyl to hydroxymethyl, leaving formaldehyde and generating demethylated base. The AlkB enzyme from Escherichia coli has nine human homologs (ALKBH1-8 and FTO) and the TET family includes three members, TET1 to 3. Among them, some enzymes have been carefully studied, but for certain enzymes, few studies have been carried out. This review focuses on the kinetic properties of those 2-OG/Fe(II)-dependent enzymes and their alkyl substrates. We also provide some discussions on the future directions of this field.

Behind Base J: The Roles of JBP1 and JBP2 on Trypanosomatids

β-D-glucopyranosyloxymethiluracil (Base J) is a modified thymidine base found in kinetoplastids and some related organisms. Interestingly, Base J distribution into the genome can vary depending on the organism and its life stage. Base J is reported to be found mostly at telomeric repeats, on inactive variant surface glycoproteins (VSG’s) expression sites (e.g., T. brucei), in RNA polymerase II termination sites and sub-telomeric regions (e.g., Leishmania). This hypermodified nucleotide is synthesized in two steps with the participation of two distinct thymidine hydroxylases, J-binding protein 1 and 2 (JBP1 and JBP2, respectively) and a β-glucosyl transferase. A third J-binding protein, named JBP3, was recently identified as part of a multimeric complex. Although its structural similarities with JBP1, it seems not to be involved in J biosynthesis but to play roles in gene expression regulation in trypanosomatids. Over the years, with the characterization of JBP1 and JBP2 mutant lines, Base J functions have been targeted and shone a light on that matter, showing genus-specific features. This review aims to explore Base J’s reported participation as a regulator of RNA polymerase II transcription termination and to summarize the functional and structural characteristics and similarities of the remarkable JBP proteins in pathogenic trypanosomatids.

The RNA-binding protein RBP33 dampens non-productive transcription in trypanosomes

In-depth analysis of the transcriptomes of several model organisms has revealed that genomes are pervasively transcribed, giving rise to an abundance of non-canonical and mainly antisense RNA polymerase II-derived transcripts that are produced from almost any genomic context. Pervasive RNAs are degraded by surveillance mechanisms, but the repertoire of proteins that control the fate of these non-productive transcripts is still incomplete. Trypanosomes are single-celled eukaryotes that show constitutive RNA polymerase II transcription and in which initiation and termination of transcription occur at a limited number of sites per chromosome. It is not known whether pervasive transcription exists in organisms with unregulated RNA polymerase II activity, and which factors could be involved in the process. We show here that depletion of RBP33 results in overexpression of ∼40% of all annotated genes in the genome, with a marked accumulation of sense and antisense transcripts derived from silenced regions. RBP33 loss does not result in a significant increase in chromatin accessibility. Finally, we have found that transcripts that increase in abundance upon RBP33 knockdown are significantly more stable in RBP33-depleted trypanosomes, and that the exosome complex is responsible for their degradation. Our results provide strong evidence that RBP33 dampens non-productive transcription in trypanosomes.

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.

Open chromatin analysis in Trypanosoma cruzi life forms highlights critical differences in genomic compartments and developmental regulation at tDNA loci

BACKGROUND

Genomic organization and gene expression regulation in trypanosomes are remarkable because protein-coding genes are organized into codirectional gene clusters with unrelated functions. Moreover, there is no dedicated promoter for each gene, resulting in polycistronic gene transcription, with posttranscriptional control playing a major role. Nonetheless, these parasites harbor epigenetic modifications at critical regulatory genome features that dynamically change among parasite stages, which are not fully understood.

RESULTS

Here, we investigated the impact of chromatin changes in a scenario commanded by posttranscriptional control exploring the parasite Trypanosoma cruzi and its differentiation program using FAIRE-seq approach supported by transmission electron microscopy. We identified differences in T. cruzi genome compartments, putative transcriptional start regions, and virulence factors. In addition, we also detected a developmental chromatin regulation at tRNA loci (tDNA), which could be linked to the intense chromatin remodeling and/or the translation regulatory mechanism required for parasite differentiation. We further integrated the open chromatin profile with public transcriptomic and MNase-seq datasets. Strikingly, a positive correlation was observed between active chromatin and steady-state transcription levels.

CONCLUSION

Taken together, our results indicate that chromatin changes reflect the unusual gene expression regulation of trypanosomes and the differences among parasite developmental stages, even in the context of a lack of canonical transcriptional control of protein-coding genes.

H2B.V demarcates divergent strand-switch regions, some tDNA loci, and genome compartments in Trypanosoma cruzi and affects parasite differentiation and host cell invasion

Histone variants play a crucial role in chromatin structure organization and gene expression. Trypanosomatids have an unusual H2B variant (H2B.V) that is known to dimerize with the variant H2A.Z generating unstable nucleosomes. Previously, we found that H2B.V protein is enriched in tissue-derived trypomastigote (TCT) life forms, a nonreplicative stage of Trypanosoma cruzi, suggesting that this variant may contribute to the differences in chromatin structure and global transcription rates observed among parasite life forms. Here, we performed the first genome-wide profiling of histone localization in T. cruzi using epimastigotes and TCT life forms, and we found that H2B.V was preferentially located at the edges of divergent transcriptional strand switch regions, which encompass putative transcriptional start regions; at some tDNA loci; and between the conserved and disrupted genome compartments, mainly at trans-sialidase, mucin and MASP genes. Remarkably, the chromatin of TCT forms was depleted of H2B.V-enriched peaks in comparison to epimastigote forms. Interactome assays indicated that H2B.V associated specifically with H2A.Z, bromodomain factor 2, nucleolar proteins and a histone chaperone, among others. Parasites expressing reduced H2B.V levels were associated with higher rates of parasite differentiation and mammalian cell infectivity. Taken together, H2B.V demarcates critical genomic regions and associates with regulatory chromatin proteins, suggesting a scenario wherein local chromatin structures associated with parasite differentiation and invasion are regulated during the parasite life cycle.

A systematic analysis of Trypanosoma brucei chromatin factors identifies novel protein interaction networks associated with sites of transcription initiation and termination

Nucleosomes composed of histones are the fundamental units around which DNA is wrapped to form chromatin. Transcriptionally active euchromatin or repressive heterochromatin is regulated in part by the addition or removal of histone post-translational modifications (PTMs) by "writer" and "eraser" enzymes, respectively. Nucleosomal PTMs are recognized by a variety of "reader" proteins that alter gene expression accordingly. The histone tails of the evolutionarily divergent eukaryotic parasite Trypanosoma brucei have atypical sequences and PTMs distinct from those often considered universally conserved. Here we identify 65 predicted readers, writers, and erasers of histone acetylation and methylation encoded in the T. brucei genome and, by epitope tagging, systemically localize 60 of them in the parasite's bloodstream form. ChIP-seq shows that 15 candidate proteins associate with regions of RNAPII transcription initiation. Eight other proteins show a distinct distribution with specific peaks at a subset of RNAPII transcription termination regions marked by RNAPIII-transcribed tRNA and snRNA genes. Proteomic analyses identify distinct protein interaction networks comprising known chromatin regulators and novel trypanosome-specific components. Notably, several SET- and Bromo-domain protein networks suggest parallels to RNAPII promoter-associated complexes in conventional eukaryotes. Further, we identify likely components of TbSWR1 and TbNuA4 complexes whose enrichment coincides with the SWR1-C exchange substrate H2A.Z at RNAPII transcription start regions. The systematic approach used provides details of the composition and organization of the chromatin regulatory machinery in T. brucei and establishes a route to explore divergence from eukaryotic norms in an evolutionarily ancient but experimentally accessible eukaryote.

Genetic Interaction Between Site-Specific Epigenetic Marks and Roles of H4v in Transcription Termination in Trypanosoma brucei

In Trypanosoma brucei, genes are assembled in polycistronic transcription units (PTUs). Boundaries of PTUs are designated transcription start sites and transcription termination sites (TTSs). Messenger RNAs are generated by trans-splicing and polyadenylation of precursor RNAs, and regulatory information in the 3' un-translated region (UTR), rather than promoter activity/sequence-specific transcription factors, controls mRNA levels. Given this peculiar genome structure, special strategies must be utilized to control transcription in T. brucei. TTSs are deposition sites for three non-essential chromatin factors-two of non-canonical histone variants (H3v and H4v) and a DNA modification (base J, which is a hydroxyl-glucosyl dT). This association generated the hypothesis that these three chromatin marks define a transcription termination site in T. brucei. Using a panel of null mutants lacking H3v, H4v, and base J, here I show that H4v is a major sign for transcription termination at TTSs. While having a secondary function at TTSs, H3v is important for monoallelic transcription of telomeric antigen genes. The simultaneous absence of both histone variants leads to proliferation and replication defects, which are exacerbated by the J absence, accompanied by accumulation of sub-G1 population. Thus, I propose that the coordinated actions of H3v, H4v, and J provide compensatory mechanisms for each other in chromatin organization, transcription, replication, and cell-cycle progression.

Quantification and mapping of DNA modifications

Apart from the four canonical nucleobases, DNA molecules carry a number of natural modifications. Substantial evidence shows that DNA modifications can regulate diverse biological processes. Dynamic and reversible modifications of DNA are critical for cell differentiation and development. Dysregulation of DNA modifications is closely related to many human diseases. The research of DNA modifications is a rapidly expanding area and has been significantly stimulated by the innovations of analytical methods. With the recent advances in methods and techniques, a series of new DNA modifications have been discovered in the genomes of prokaryotes and eukaryotes. Deciphering the biological roles of DNA modifications depends on the sensitive detection, accurate quantification, and genome-wide mapping of modifications in genomic DNA. This review provides an overview of the recent advances in analytical methods and techniques for both the quantification and genome-wide mapping of natural DNA modifications. We discuss the principles, advantages, and limitations of these developed methods. It is anticipated that new methods and techniques will resolve the current challenges in this burgeoning research field and expedite the elucidation of the functions of DNA modifications.

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.

Strain-specific genome evolution in Trypanosoma cruzi, the agent of Chagas disease

The protozoan Trypanosoma cruzi almost invariably establishes life-long infections in humans and other mammals, despite the development of potent host immune responses that constrain parasite numbers. The consistent, decades-long persistence of T. cruzi in human hosts arises at least in part from the remarkable level of genetic diversity in multiple families of genes encoding the primary target antigens of anti-parasite immune responses. However, the highly repetitive nature of the genome-largely a result of these same extensive families of genes-have prevented a full understanding of the extent of gene diversity and its maintenance in T. cruzi. In this study, we have combined long-read sequencing and proximity ligation mapping to generate very high-quality assemblies of two T. cruzi strains representing the apparent ancestral lineages of the species. These assemblies reveal not only the full repertoire of the members of large gene families in the two strains, demonstrating extreme diversity within and between isolates, but also provide evidence of the processes that generate and maintain that diversity, including extensive gene amplification, dispersion of copies throughout the genome and diversification via recombination and in situ mutations. Gene amplification events also yield significant copy number variations in a substantial number of genes presumably not required for or involved in immune evasion, thus forming a second level of strain-dependent variation in this species. The extreme genome flexibility evident in T. cruzi also appears to create unique challenges with respect to preserving core genome functions and gene expression that sets this species apart from related kinetoplastids.

Identification of a novel base J binding protein complex involved in RNA polymerase II transcription termination in trypanosomes

Base J, β-D-glucosyl-hydroxymethyluracil, is a modification of thymine DNA base involved in RNA Polymerase (Pol) II transcription termination in kinetoplastid protozoa. Little is understood regarding how specific thymine residues are targeted for J-modification or the mechanism of J regulated transcription termination. To identify proteins involved in J-synthesis, we expressed a tagged version of the J-glucosyltransferase (JGT) in Leishmania tarentolae, and identified four co-purified proteins by mass spectrometry: protein phosphatase (PP1), a homolog of Wdr82, a potential PP1 regulatory protein (PNUTS) and a protein containing a J-DNA binding domain (named JBP3). Gel shift studies indicate JBP3 is a J-DNA binding protein. Reciprocal tagging, co-IP and sucrose gradient analyses indicate PP1, JGT, JBP3, Wdr82 and PNUTS form a multimeric complex in kinetoplastids, similar to the mammalian PTW/PP1 complex involved in transcription termination via PP1 mediated dephosphorylation of Pol II. Using RNAi and analysis of Pol II termination by RNA-seq and RT-PCR, we demonstrate that ablation of PNUTS, JBP3 and Wdr82 lead to defects in Pol II termination at the 3'-end of polycistronic gene arrays in Trypanosoma brucei. Mutants also contain increased antisense RNA levels upstream of transcription start sites, suggesting an additional role of the complex in regulating termination of bi-directional transcription. In addition, PNUTS loss causes derepression of silent Variant Surface Glycoprotein genes involved in host immune evasion. Our results suggest a novel mechanistic link between base J and Pol II polycistronic transcription termination in kinetoplastids.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

Trypanosoma brucei ribonuclease H2A is an essential R-loop processing enzyme whose loss causes DNA damage during transcription initiation and antigenic variation

Ribonucleotides represent a threat to DNA genome stability and transmission. Two types of Ribonuclease H (RNase H) excise ribonucleotides when they form part of the DNA strand, or hydrolyse RNA when it base-pairs with DNA in structures termed R-loops. Loss of either RNase H is lethal in mammals, whereas yeast survives the absence of both enzymes. RNase H1 loss is tolerated by the parasite Trypanosoma brucei but no work has examined the function of RNase H2. Here we show that loss of T. brucei RNase H2 (TbRH2A) leads to growth and cell cycle arrest that is concomitant with accumulation of nuclear damage at sites of RNA polymerase (Pol) II transcription initiation, revealing a novel and critical role for RNase H2. Differential gene expression analysis reveals limited overall changes in RNA levels for RNA Pol II genes after TbRH2A loss, but increased perturbation of nucleotide metabolic genes. Finally, we show that TbRH2A loss causes R-loop and DNA damage accumulation in telomeric RNA Pol I transcription sites, also leading to altered gene expression. Thus, we demonstrate separation of function between two nuclear T. brucei RNase H enzymes during RNA Pol II transcription, but overlap in function during RNA Pol I-mediated gene expression during host immune evasion.

The domain architecture of the protozoan protein J-DNA-binding protein 1 suggests synergy between base J DNA binding and thymidine hydroxylase activity

J-DNA-binding protein 1 (JBP1) contributes to the biosynthesis and maintenance of base J (β-d-glucosyl-hydroxymethyluracil), an epigenetic modification of thymidine (T) confined to pathogenic protozoa such as Trypanosoma and Leishmania JBP1 has two known functional domains: an N-terminal T hydroxylase (TH) homologous to the 5-methylcytosine hydroxylase domain in TET proteins and a J-DNA-binding domain (JDBD) that resides in the middle of JBP1. Here, we show that removing JDBD from JBP1 results in a soluble protein (Δ-JDBD) with the N- and C-terminal regions tightly associated together in a well-ordered structure. We found that this Δ-JDBD domain retains TH activity in vitro but displays a 15-fold lower apparent rate of hydroxylation compared with JBP1. Small-angle X-ray scattering (SAXS) experiments on JBP1 and JDBD in the presence or absence of J-DNA and on Δ-JDBD enabled us to generate low-resolution three-dimensional models. We conclude that Δ-JDBD, and not the N-terminal region of JBP1 alone, is a distinct folding unit. Our SAXS-based model supports the notion that binding of JDBD specifically to J-DNA can facilitate T hydroxylation 12-14 bp downstream on the complementary strand of the J-recognition site. We postulate that insertion of the JDBD module into the Δ-JDBD scaffold during evolution provided a mechanism that synergized J recognition and T hydroxylation, ensuring inheritance of base J in specific sequence patterns following DNA replication in kinetoplastid parasites.

Regulation of gene expression in trypanosomatids: living with polycistronic transcription

In trypanosomes, RNA polymerase II transcription is polycistronic and individual mRNAs are excised by trans-splicing and polyadenylation. The lack of individual gene transcription control is compensated by control of mRNA processing, translation and degradation. Although the basic mechanisms of mRNA decay and translation are evolutionarily conserved, there are also unique aspects, such as the existence of six cap-binding translation initiation factor homologues, a novel decapping enzyme and an mRNA stabilizing complex that is recruited by RNA-binding proteins. High-throughput analyses have identified nearly a hundred regulatory mRNA-binding proteins, making trypanosomes valuable as a model system to investigate post-transcriptional regulation.

Characterization and structure determination of a llama-derived nanobody targeting the J-base binding protein 1

J-base binding protein 1 (JBP1) contributes to the biosynthesis and maintenance of base J (β-D-glucosylhydroxymethyluracil), a modification of thymidine confined to some protozoa. Camelid (llama) single-domain antibody fragments (nanobodies) targeting JBP1 were produced for use as crystallization chaperones. Surface plasmon resonance screening identified Nb6 as a strong binder, recognizing JBP1 with a 1:1 stoichiometry and high affinity (Kd = 30 nM). Crystallization trials of JBP1 in complex with Nb6 yielded crystals that diffracted to 1.47 Å resolution. However, the dimensions of the asymmetric unit and molecular replacement with a nanobody structure clearly showed that the crystals of the expected complex with JBP1 were of the nanobody alone. Nb6 crystallizes in space group P31 with two molecules in the asymmetric unit; its crystal structure was refined to a final resolution of 1.64 Å. Ensemble refinement suggests that in the ligand-free state one of the complementarity-determining regions (CDRs) is flexible, while the other two adopt well defined conformations.

Fluorescence-based assay for accurate measurement of transcriptional activity in trypanosomatid parasites

Trypanosomatid parasites are causative agents of neglected human diseases. Their lineage diverged early from the common eukaryotic ancestor, and they evolved singular mechanisms of gene expression that are crucial for their survival. Studies on unusual and essential molecular pathways lead to new drug targets. In this respect, assays to analyze transcriptional activity will provide useful information to identify essential and specific factors. However, the current methods are laborious and do not provide global and accurate measures. For this purpose, a previously reported radiolabeling in vitro nascent mRNA methodology was used to establish an alternative fluorescent-based assay that is able to precisely quantify nascent mRNA using both flow cytometry and a high-content image system. The method allowed accurate and global measurements in Trypanosoma brucei, a representative species of trypanosomatid parasites. We obtained data demonstrating that approximately 70% of parasites from a population under normal growth conditions displayed mRNA transcriptional activity, whilst the treatment with α-amanitin (75 µg/ml) inhibited the polymerase II activity. The adaptation of the method also allowed the analyses of the transcriptional activity during the cell cycle. Therefore, the methodology described herein contributes to obtaining precise measurements of transcriptional rates using multiparametric analysis. This alternative method can facilitate investigations of genetic and biochemical processes in trypanosome parasites and consequently provide additional information related to new treatment or prophylaxis strategies involving these important human parasites.

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.

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.

A method for the efficient and selective identification of 5-hydroxymethyluracil in genomic DNA

Recently, 5-hydroxymethyluracil (5hmU) was identified in mammalian genomic DNA as an oxidative product of thymine by the ten-eleven translocation (TET) proteins. While the biological role of this modification remains unclear, identifying its genomic location will assist in elucidating function. Here we present a rapid and robust method to selectively tag and enrich genomic regions containing 5hmU. This method involves the selective glucosylation of 5hmU residues by the base J glucosyltransferase from trypanosomes creating glucosylhydroxymethyluracil (base J). The base J can then be efficiently and selectively pulled down by antibodies against base J or by J-binding protein 1. DNA that is enriched is suitable for analysis by quantitative PCR or sequencing. We utilized this tagging reaction to provide proof of concept for the enrichment of 5hmU containing DNA from a pool that contains modified and unmodified DNA. Furthermore, we demonstrate that the base J pull-down assay identifies 5hmU at specific regions of the trypanosome genome involved in transcriptional repression. The method described here will allow for a greater understanding of the functional role and dynamics of 5hmU in biology.

Base J represses genes at the end of polycistronic gene clusters in Leishmania major by promoting RNAP II termination

The genomes of kinetoplastids are organized into polycistronic gene clusters that are flanked by the modified DNA base J. Previous work has established a role of base J in promoting RNA polymerase II termination in Leishmania spp. where the loss of J leads to termination defects and transcription into adjacent gene clusters. It remains unclear whether these termination defects affect gene expression and whether read through transcription is detrimental to cell growth, thus explaining the essential nature of J. We now demonstrate that reduction of base J at specific sites within polycistronic gene clusters in L. major leads to read through transcription and increased expression of downstream genes in the cluster. Interestingly, subsequent transcription into the opposing polycistronic gene cluster does not lead to downregulation of sense mRNAs. These findings indicate a conserved role for J regulating transcription termination and expression of genes within polycistronic gene clusters in trypanosomatids. In contrast to the expectations often attributed to opposing transcription, the essential nature of J in Leishmania spp. is related to its role in gene repression rather than preventing transcriptional interference resulting from read through and dual strand transcription.

Base J glucosyltransferase does not regulate the sequence specificity of J synthesis in trypanosomatid telomeric DNA

Telomeric DNA of trypanosomatids possesses a modified thymine base, called base J, that is synthesized in a two-step process; the base is hydroxylated by a thymidine hydroxylase forming hydroxymethyluracil (hmU) and a glucose moiety is then attached by the J-associated glucosyltransferase (JGT). To examine the importance of JGT in modifiying specific thymine in DNA, we used a Leishmania episome system to demonstrate that the telomeric repeat (GGGTTA) stimulates J synthesis in vivo while mutant telomeric sequences (GGGTTT, GGGATT, and GGGAAA) do not. Utilizing an in vitro GT assay we find that JGT can glycosylate hmU within any sequence with no significant change in Km or kcat, even mutant telomeric sequences that are unable to be J-modified in vivo. The data suggests that JGT possesses no DNA sequence specificity in vitro, lending support to the hypothesis that the specificity of base J synthesis is not at the level of the JGT reaction.

Base J and H3.V Regulate Transcriptional Termination in Trypanosoma brucei

Trypanosoma brucei is a protozoan parasite that lacks many transcription factors found in other eukaryotes, such as those whose binding demarcates enhancers. T. brucei retains histone variants and modifications, however, and it is hypothesized that it relies on epigenetic marks to define transcription-related boundaries. The histone H3 variant (H3.V) and an alternate nucleotide, base J (ß-D-glucosyl-hydroxymethyluracil), are two chromatin marks found at both transcription termination sites (TTSs) and telomeres. Here, we report that the absence of both base J and H3.V result in transcription readthrough and the appearance of antisense transcripts near TTSs. Additionally, we find that maintaining the transcriptional silencing of pol I-transcribed telomeric Variant Surface Glycoprotein (VSG) genes appears to be dependent on deposition of H3.V alone. Our study reveals that gene expression depends on different epigenetic cues depending on chromosomal location and on the transcribing polymerase. This work provides insight into how these signals may have evolved into the more nuanced and fine-tuned gene regulatory mechanisms observed in other model systems.

In eukaryotes, DNA is wrapped around histones to form chromatin. Modifications on the DNA itself, or on the canonical histones around which the DNA is wrapped, can lead to changes in gene expression. Alternate histones are also used to regulate gene expression. The African trypanosome, the causative agent of African sleeping sickness, transcribes its genes in long arrays called polycistronic transcription units (PTUs). In this study, we describe how the African Trypanosome uses two chromatin marks to regulate genes that lie close to the end of the PTU or close to the ends of chromosomes, called telomeres. One mark, base J, is on the DNA itself, while the other utilized mark is an alternate histone protein called H3.V. In the absence of these marks, there is an increase in antisense RNA that corresponds to genes that lie near the ends of the PTUs, and an increase in sense RNA for genes near telomeres. Since trypanosomes are evolutionarily distant from humans, these findings shed light on how gene expression mechanisms that are modulated by chromatin marks may have evolved to produce the complex gene regulatory networks found in our own tissues.

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

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

The response of trypanosomes and other eukaryotes to ER stress and the spliced leader RNA silencing (SLS) pathway in Trypanosoma brucei

The unfolded protein response (UPR) is induced when the quality control machinery of the cell is overloaded with unfolded proteins or when one of the functions of the endoplasmic reticulum (ER) is perturbed. Here, I describe UPR in yeast and mammals, and compare it to what we know about pathogenic fungi and the parasitic protozoans from the order kinetoplastida, focusing on the novel pathway the spliced leader silencing (SLS) in Trypanosoma brucei. Trypanosomes lack conventional transcription regulation, and thus, lack most of the UPR machinery present in other eukaryotes. Trypanosome genes are transcribed in polycistronic units that are processed by trans-splicing and polyadenylation. In trans-splicing, which is essential for processing of each mRNA, an exon known as the spliced leader (SL) is added to all mRNAs from a small RNA, the SL RNA. Under severe ER stress, T. brucei elicits the SLS pathway. In SLS, the transcription of the SL RNA gene is extinguished, and the entire transcription complex dissociates from the SL RNA promoter. Induction of SLS is mediated by an ER-associated kinase (PK3) that migrates to the nucleus, where it phosphorylates the TATA-binding protein (TRF4), leading shut-off of SL RNA transcription. As a result, trans-splicing is inhibited and the parasites activate a programmed cell death (PCD) pathway. Despite the ability to sense the ER stress, the different eukaryotes, especially unicellular parasites and pathogenic fungi, developed a variety of unique and different ways to sense and adjust to this stress in a manner different from their host.

The TET/JBP Family of Nucleic Acid Base-Modifying 2-Oxoglutarate and Iron-Dependent Dioxygenases

The TET/JBP family of enzymes includes 2-oxoglutarate- and Fe(ii)-dependent dioxygenases that oxidize 5-methylpyrimidines in nucleic acids. They include euglenozoan JBP enzymes that catalyse the first step in the biosynthesis of the hypermodified thymine, base J, and metazoan TET enzymes that generate oxidized 5-methylcytosines (hydroxy-, formyl- and carboxymethylcytosine) in DNA. Recent studies suggest that these modified bases function as epigenetic marks and/or as potential intermediates for DNA demethylation during resetting of epigenetic 5mC marks upon zygote formation and in primordial germ cell development. Studies in mammalian models also point to an important role for these enzymes in haematopoiesis, tumour suppression, cell differentiation and neural behavioural adaptation. The TET/JBP family has undergone extensive gene expansion in fungi, such as mushrooms, in conjunction with a novel class of transposons and might play a role in genomic plasticity and speciation. Certain versions from stramenopiles and chlorophytes are likely to modify RNA and often show fusions to other RNA-modifying enzymatic domains. The ultimate origin of the TET/JBP family lies in bacteriophages where the enzymes are likely to catalyse formation of modified bases with key roles in DNA packaging and evasion of host restriction.

25 years of African trypanosome research: From description to molecular dissection and new drug discovery

The Molecular Parasitology conference was first held at the Marine Biological laboratory, Woods Hole, USA 25 years ago. Since that first meeting, the conference has evolved and expanded but has remained the showcase for the latest research developments in molecular parasitology. In this perspective, I reflect on the scientific discoveries focussed on African trypanosomes (Trypanosoma brucei spp.) that have occurred since the inaugural MPM meeting and discuss the current and future status of research on these parasites.

Defining the sequence requirements for the positioning of base J in DNA using SMRT sequencing

Base J (β-D-glucosyl-hydroxymethyluracil) replaces 1% of T in the Leishmania genome and is only found in telomeric repeats (99%) and in regions where transcription starts and stops. This highly restricted distribution must be co-determined by the thymidine hydroxylases (JBP1 and JBP2) that catalyze the initial step in J synthesis. To determine the DNA sequences recognized by JBP1/2, we used SMRT sequencing of DNA segments inserted into plasmids grown in Leishmania tarentolae. We show that SMRT sequencing recognizes base J in DNA. Leishmania DNA segments that normally contain J also picked up J when present in the plasmid, whereas control sequences did not. Even a segment of only 10 telomeric (GGGTTA) repeats was modified in the plasmid. We show that J modification usually occurs at pairs of Ts on opposite DNA strands, separated by 12 nucleotides. Modifications occur near G-rich sequences capable of forming G-quadruplexes and JBP2 is needed, as it does not occur in JBP2-null cells. We propose a model whereby de novo J insertion is mediated by JBP2. JBP1 then binds to J and hydroxylates another T 13 bp downstream (but not upstream) on the complementary strand, allowing JBP1 to maintain existing J following DNA replication.

Suppression of subtelomeric VSG switching by Trypanosoma brucei TRF requires its TTAGGG repeat-binding activity

Trypanosoma brucei causes human African trypanosomiasis and regularly switches its major surface antigen, VSG, in the bloodstream of its mammalian host to evade the host immune response. VSGs are expressed exclusively from subtelomeric loci, and we have previously shown that telomere proteins TbTIF2 and TbRAP1 play important roles in VSG switching and VSG silencing regulation, respectively. We now discover that the telomere duplex DNA-binding factor, TbTRF, also plays a critical role in VSG switching regulation, as a transient depletion of TbTRF leads to significantly more VSG switching events. We solved the NMR structure of the DNA-binding Myb domain of TbTRF, which folds into a canonical helix-loop-helix structure that is conserved to the Myb domains of mammalian TRF proteins. The TbTRF Myb domain tolerates well the bulky J base in T. brucei telomere DNA, and the DNA-binding affinity of TbTRF is not affected by the presence of J both in vitro and in vivo. In addition, we find that point mutations in TbTRF Myb that significantly reduced its in vivo telomere DNA-binding affinity also led to significantly increased VSG switching frequencies, indicating that the telomere DNA-binding activity is critical for TbTRF's role in VSG switching regulation.

Quantitative mass spectrometry-based analysis of β-D-glucosyl-5-hydroxymethyluracil in genomic DNA of Trypanosoma brucei

β-D-glucosyl-5-hydroxymethyluracil (base J) is a hyper-modified nucleobase found in the nuclear DNA of kinetoplastid parasites. With replacement of a fraction of thymine in DNA, J is localized primarily in telomeric regions of all organisms carrying this modified base. The biosynthesis of J occurs in two putative steps: first, a specific thymine in DNA is recognized and converted into 5-hydroxymethyluracil (5-HmU) by J-binding proteins (JBP1 and JBP2); a glucosyl transferase (GT) subsequently glucosylates the 5-HmU to yield J. Although several recent studies revealed the roles of internal J in regulating transcription in kinetoplastids, functions of telomeric J and proteins involved in J synthesis remain elusive. Assessing the functions of base J and understanding fully its biosynthesis necessitate the measurement of its level in cells and organisms. In this study, we reported a reversed-phase HPLC coupled with tandem mass spectrometry (LC-MS/MS) method, together with the use of a surrogate internal standard (β-D-glucosyl-5-hydroxymethyl-2'-deoxycytidine, 5-gHmdC), for the accurate detection of β-D-glucosyl-5-hydroxymethyl-2'-deoxyuridine (dJ) in Trypanosoma brucei DNA. For comparison, we also measured the level of the precursor for dJ synthesis [i.e. 5-hydroxymethyl-2'-deoxyuridine (5-HmdU)]. We found that base J was not detectable in the JBP-null cells whereas it replaced approximately 0.5% thymine in wild-type cells, which was accompanied with a markedly decreased level of 5-HmdU in JBP1/JBP2-null strain relative to the wild-type strain. These results provided direct evidence supporting that JBP proteins play an important role in oxidizing thymidine to form 5-HmdU, which facilitated the generation of dJ. This is the first report about the application of LC-MS/MS for the quantification of base J. The analytical method built a solid foundation for dissecting the molecular mechanisms of J biosynthesis and assessing the biological functions of base J in the future.ᅟ