Physical and transcriptional analysis of the Trypanosoma brucei genome reveals a typical eukaryotic arrangement with close interspersionof RNA polymerase II- and III-transcribed genes

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.

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.

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.

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.

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.

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.

Nuclear architecture, genome and chromatin organisation in Trypanosoma brucei

The nucleus of the human pathogen Trypanosoma brucei not only has unusual chromosomal composition, characterised by the presence of megabase, intermediate and minichromosomes, but also chromosome and gene organisation that is unique amongst eukaryotes. Here I provide an overview of current knowledge of nuclear structure, chromatin organisation and chromosome dynamics during interphase and mitosis. New technologies such as chromatin immunoprecipitation, in combination with new generation sequencing and proteomic analysis of subnuclear fractions, have led to novel insights into the organisation of the nucleus and chromatin. In particular, we are beginning to understand how universal mechanisms of chromatin modifications and nuclear position effects are deployed for parasite-specific functions and are centrally involved in genomic organisation and transcriptional regulation. These advances also have a major impact on progress in understanding the molecular basis of antigenic variation.

The transcriptome of the human pathogen Trypanosoma brucei at single-nucleotide resolution

The genome of Trypanosoma brucei, the causative agent of African trypanosomiasis, was published five years ago, yet identification of all genes and their transcripts remains to be accomplished. Annotation is challenged by the organization of genes transcribed by RNA polymerase II (Pol II) into long unidirectional gene clusters with no knowledge of how transcription is initiated. Here we report a single-nucleotide resolution genomic map of the T. brucei transcriptome, adding 1,114 new transcripts, including 103 non-coding RNAs, confirming and correcting many of the annotated features and revealing an extensive heterogeneity of 5′ and 3′ ends. Some of the new transcripts encode polypeptides that are either conserved in T. cruzi and Leishmania major or were previously detected in mass spectrometry analyses. High-throughput RNA sequencing (RNA-Seq) was sensitive enough to detect transcripts at putative Pol II transcription initiation sites. Our results, as well as recent data from the literature, indicate that transcription initiation is not solely restricted to regions at the beginning of gene clusters, but may occur at internal sites. We also provide evidence that transcription at all putative initiation sites in T. brucei is bidirectional, a recently recognized fundamental property of eukaryotic promoters. Our results have implications for gene expression patterns in other important human pathogens with similar genome organization (Trypanosoma cruzi, Leishmania sp.) and revealed heterogeneity in pre-mRNA processing that could potentially contribute to the survival and success of the parasite population in the insect vector and the mammalian host.

Identifying genes essential for survival in the host is fundamental to unraveling the biology of human pathogens and understanding mechanisms of pathogenesis. The protozoan parasite Trypanosoma brucei causes devastating diseases in humans and animals in sub-Saharan Africa, and the publication in 2005 of the genome sequence provided the first glance at the coding potential of this organism. Although at present there is a catalogue of predicted protein coding genes, the challenge remains to identify all authentic genes, including their boundaries. We used next generation RNA sequencing (RNA-Seq) to map transcribed regions and RNA polymerase II transcription initiation sites on a genome-wide scale. This approach allowed us to improve and correct the current annotation, to reveal a widespread heterogeneity of RNA processing sites (trans-splicing and polyadenylation) and to estimate that most genes are expressed at levels corresponding to 1 to 10 mRNAs per cell. Our data indicate that different transcript forms representing the same gene are present stochastically within the mRNA population. This unanticipated scenario may contribute to determining gene expression landscapes to adapt to different environments in the parasite life cycle.

U3 snoRNA genes are multi-copy and frequently linked to U5 snRNA genes in Euglena gracilis

Background

U3 snoRNA is a box C/D small nucleolar RNA (snoRNA) involved in the processing events that liberate 18S rRNA from the ribosomal RNA precursor (pre-rRNA). Although U3 snoRNA is present in all eukaryotic organisms, most investigations of it have focused on fungi (particularly yeasts), animals and plants. Relatively little is known about U3 snoRNA and its gene(s) in the phylogenetically broad assemblage of protists (mostly unicellular eukaryotes). In the euglenozoon Euglena gracilis, a distant relative of the kinetoplastid protozoa, Southern analysis had previously revealed at least 13 bands hybridizing with U3 snoRNA, suggesting the existence of multiple copies of U3 snoRNA genes.

Results

Through screening of a λ genomic library and PCR amplification, we recovered 14 U3 snoRNA gene variants, defined by sequence heterogeneities that are mostly located in the U3 3'-stem-loop domain. We identified three different genomic arrangements of Euglena U3 snoRNA genes: i) stand-alone, ii) linked to tRNA^Arg genes, and iii) linked to a U5 snRNA gene. In arrangement ii), the U3 snoRNA gene is positioned upstream of two identical tRNA^Arg genes that are convergently transcribed relative to the U3 gene. This scenario is reminiscent of a U3 snoRNA-tRNA gene linkage previously described in trypanosomatids. We document here twelve different U3 snoRNA-U5 snRNA gene arrangements in Euglena; in each case, the U3 gene is linked to a downstream and convergently oriented U5 gene, with the intergenic region differing in length and sequence among the variants.

Conclusion

The multiple U3 snoRNA-U5 snRNA gene linkages, which cluster into distinct families based on sequence similarities within the intergenic spacer, presumably arose by genome, chromosome, and/or locus duplications. We discuss possible reasons for the existence of the unusually large number of U3 snoRNA genes in the Euglena genome. Variability in the signal intensities of the multiple Southern hybridization bands raises the possibility that Euglena contains a naturally aneuploid chromosome complement.

Four histone variants mark the boundaries of polycistronic transcription units in Trypanosoma brucei

Unusually for a eukaryote, genes transcribed by RNA polymerase II (pol II) in Trypanosoma brucei are arranged in polycistronic transcription units. With one exception, no pol II promoter motifs have been identified, and how transcription is initiated remains an enigma. T. brucei has four histone variants: H2AZ, H2BV, H3V, and H4V. Using chromatin immunoprecipitation (ChIP) and sequencing (ChIP-seq) to examine the genome-wide distribution of chromatin components, we show that histones H4K10ac, H2AZ, H2BV, and the bromodomain factor BDF3 are enriched up to 300-fold at probable pol II transcription start sites (TSSs). We also show that nucleosomes containing H2AZ and H2BV are less stable than canonical nucleosomes. Our analysis also identifies >60 unexpected TSS candidates and reveals the presence of long guanine runs at probable TSSs. Apparently unique to trypanosomes, additional histone variants H3V and H4V are enriched at probable pol II transcription termination sites. Our findings suggest that histone modifications and histone variants play crucial roles in transcription initiation and termination in trypanosomes and that destabilization of nucleosomes by histone variants is an evolutionarily ancient and general mechanism of transcription initiation, demonstrated in an organism in which general pol II transcription factors have been elusive.

Cloning and molecular characterization of Trypanosoma cruzi U2, U4, U5, and U6 small nuclear RNAs

Small nuclear RNAs (snRNAs) are important factors in the functioning of eukaryotic cells that form several small complexes with proteins; these ribonucleoprotein particles (U snRNPs) have an essential role in the pre-mRNA processing, particularly in splicing, catalyzed by spliceosomes, large RNA-protein complexes composed of various snRNPs. Even though they are well defined in mammals, snRNPs are still not totally characterized in certain trypanosomatids as Trypanosoma cruzi. For this reason we subjected snRNAs (U2, U4, U5, and U6) from T. cruzi epimastigotes to molecular characterization by polymerase chain reaction (PCR) and reverse transcription-PCR. These amplified sequences were cloned, sequenced, and compared with those other of trypanosomatids. Among these snRNAs, U5 was less conserved and U6 the most conserved. Their respective secondary structures were predicted and compared with known T. brucei structures. In addition, the copy number of each snRNA in the T. cruzi genome was characterized by Southern blotting.

Trypanosoma brucei encodes a bifunctional capping enzyme essential for cap 4 formation on the spliced leader RNA

The 5' end of kinetoplastid mRNA possesses a hypermethylated cap 4 structure, which is derived from standard m7GpppN (cap 0) with additional methylations at seven sites within the first four nucleosides on the spliced leader RNA. In addition to TbCe1 guanylyltransferase and TbCmt1 (guanine N-7) methyltransferase, Trypanosoma brucei encodes a second cap 0 forming enzyme. TbCgm1 (T. brucei cap guanylyltransferase-methyltransferase) is a novel bifunctional capping enzyme consisting of an amino-terminal guanylyltransferase domain and a carboxyl-terminal methyltransferase domain. Recombinant TbCgm1 transfers the GMP to spliced leader RNA (SL RNA) via a covalent enzyme-GMP intermediate, and methylates the guanine N-7 position of the GpppN-terminated RNA to form cap 0 structure. The two domains can function autonomously in vitro. TbCGM1 is essential for parasite growth. Silencing of TbCGM1 by RNA interference increased the abundance of uncapped SL RNA and lead to accumulation of hypomethylated SL RNA. In contrast, silencing of TbCE1 and TbCMT1 did not affect parasite growth or SL RNA capping. We conclude that TbCgm1 specifically cap SL RNA, and cap 0 is a prerequisite for subsequent methylation events leading to the formation of mature SL RNA.

Transcription initiation and termination on Leishmania major chromosome 3

Genome projects involving Leishmania and other trypanosomatids have revealed that most genes in these organisms are organized into large clusters of genes on the same DNA strand. We have previously shown that transcription of the entire Leishmania major Friedlin (LmjF) chromosome 1 (chr1) initiates bidirectionally between two divergent gene clusters. Here, we analyze transcription of LmjF chr3, which contains two convergent clusters of 67 and 30 genes, separated by a tRNA gene, with a single divergent protein-coding gene located close to the "left" telomere. Nuclear run-on analyses indicate that specific transcription of chr3 initiates bidirectionally between the single subtelomeric gene and the adjacent 67-gene cluster, close to the "right" telomere upstream of the 30-gene cluster, and upstream of the tRNA gene. Transcription on both strands terminates within the tRNA-gene region. Transient-transfection studies support the role of the tRNA-gene region as a transcription terminator for RNA polymerase II (Pol II) and Pol III, and also for Pol I.

Transcription in kinetoplastid protozoa: why be normal?

Transcription in the kinetoplastid protozoa shows substantial variation from the paradigms of eukaryotic gene expression, including polycistronic transcription, a paucity of RNA polymerase (RNAP) II promoters, no qualitative regulated transcription initiation for most protein-coding genes, transcription of some protein-coding genes by RNAP I, an exclusive subnuclear location for VSG transcription, the dependence of small nuclear RNA gene transcription on an upstream tRNA gene, and the synthesis of mitochondrial tRNAs in the nucleus. Here, we present a broad overview of what is known about transcription in the kinetoplastids and what has yet to be determined.

Delineation of the regulated Variant Surface Glycoprotein gene expression site domain of Trypanosoma brucei

The African trypanosome Trypanosoma brucei is protected in the bloodstream of the mammalian host by a dense Variant Surface Glycoprotein (VSG) coat. Although an individual cell has hundreds of VSG genes, the active VSG is transcribed in a mutually exclusive fashion from one of about twenty telomeric VSG expression sites. Expression sites are regulated domains flanked by 50 bp repeat arrays and extensive tracts of repetitive elements. We have integrated exogenous rDNA and expression site promoters upstream of the 50 bp repeats of the VO2 VSG expression site. Transcription from both types of exogenous promoter is downregulated and comparable to promoters targeted into the VSG Basic Copy arrays. We show that the upstream exogenous rDNA promoter escapes VSG expression site control, as switching the downstream VO2 VSG expression site on and off does not affect its activity. Therefore, the 50 bp repeat arrays appear to be the boundary of the regulated expression site domain.

The architecture of variant surface glycoprotein gene expression sites in Trypanosoma brucei

Trypanosoma brucei evades the immune system by switching between Variant Surface Glycoprotein (VSG) genes. The active VSG gene is transcribed in one of approximately 20 telomeric expression sites (ESs). It has been postulated that ES polymorphism plays a role in host adaptation. To gain more insight into ES architecture, we have determined the complete sequence of Bacterial Artificial Chromosomes (BACs) containing DNA from three ESs and their flanking regions. There was variation in the order and number of ES-associated genes (ESAGs). ESAGs 6 and 7, encoding transferrin receptor subunits, are the only ESAGs with functional copies in every ES that has been sequenced until now. A BAC clone containing the VO2 ES sequences comprised approximately half of a 330 kb 'intermediate' chromosome. The extensive similarity between this intermediate chromosome and the left telomere of T. brucei 927 chromosome I, suggests that this previously uncharacterised intermediate size class of chromosomes could have arisen from breakage of megabase chromosomes. Unexpected conservation of sequences, including pseudogenes, indicates that the multiple ESs could have arisen through a relatively recent amplification of a single ES.

tRNAs in Trypanosoma brucei: genomic organization, expression, and mitochondrial import

The mitochondrial genome of Trypanosoma brucei does not encode tRNAs. Consequently, all mitochondrial tRNAs are imported from the cytosol and originate from nucleus-encoded genes. Analysis of all currently available T. brucei sequences revealed that its genome carries 50 tRNA genes representing 40 different isoacceptors. The identified set is expected to be nearly complete since all but four codons are accounted for. The number of tRNA genes in T. brucei is very low for a eukaryote and lower than those of many prokaryotes. Using quantitative Northern analysis we have determined the absolute abundance in the cell and the mitochondrion of a group of 15 tRNAs specific for 12 amino acids. Except for the initiator type tRNA(Met), which is cytosol specific, the cytosolic and the mitochondrial sets of tRNAs were qualitatively identical. However, the extent of mitochondrial localization was variable for the different tRNAs, ranging from 1 to 7.5% per cell. Finally, by using transgenic cell lines in combination with quantitative Northern analysis it was shown that import of tRNA(Leu)(CAA) is independent of its 5'-genomic context, suggesting that the in vivo import substrate corresponds to the mature, fully processed tRNA.

Two novel RNA binding proteins from Trypanosoma brucei are associated with 5S rRNA

We have previously reported the identification of two closely related RNA binding proteins from Trypanosoma brucei which we have termed p34 and p37. The predicted primary structures of the two proteins are highly homologous with one major difference, an 18-amino-acid insert in the N-terminal region of p37. These two proteins have been localized to the nucleus based on immunofluorescence microscopy. To gain insight into their function, we have utilized UV crosslinking, coimmunoprecipitation, and sucrose density gradients to identify T. brucei RNA species that associate with p34 and p37. These experiments have demonstrated a specific interaction of both p34 and p37 with the 5S ribosomal RNA and indicate that other RNA species are unlikely to be specifically bound. This suggests a role for p34 and p37 in the import and/or assembly pathway of T. brucei 5S rRNA in ribosome biogenesis.

Recent developments from the Leishmania genome project

A first generation cosmid contig map of the Leishmania major Friedlin genome has been constructed, and genomic sequencing is well underway. Chromosome 1 (Chr1) and Chr3 have been completely sequenced, and Chr4 is virtually complete. Sequencing of several other chromosomes is in progress and the complete genome sequence may be available as soon as 2003. More than 600 completely sequenced new genes have been identified, representing approximately 8% of the total gene complement (approximately 8,600 genes) of Leishmania. Notably, a large proportion (approximately 69%) of the genes remain unclassified, with 40% of these being potentially Leishmania- (or kinetoplastid-) specific. Most interestingly, the genes are organized into large (>100-300 kb) polycistronic clusters of adjacent genes on the same DNA strand. Chr1 contains two such clusters organized in a 'divergent' manner, whereas Chr3 contains two 'convergent' clusters, with a single 'divergent' gene at one telomere, with the two large clusters separated by a tRNA gene. Statistical analyses of Chr1 show that the 'divergent junction' region between the two polycistronic gene clusters may be a candidate for an origin of DNA replication.