Sequence elements in both the intergenic space and the 3' untranslated region of the Crithidia fasciculata KAP3 gene are required for cell cycle regulation of KAP3 mRNA

Gene co-expression network analysis of Trypanosoma brucei in tsetse fly vector

Background

Trypanosoma brucei species are motile protozoan parasites that are cyclically transmitted by tsetse fly (genus Glossina) causing human sleeping sickness and nagana in livestock in sub-Saharan Africa. African trypanosomes display digenetic life cycle stages in the tsetse fly vector and in their mammalian host. Experimental work on insect-stage trypanosomes is challenging because of the difficulty in setting up successful in vitro cultures. Therefore, there is limited knowledge on the trypanosome biology during its development in the tsetse fly. Consequently, this limits the development of new strategies for blocking parasite transmission in the tsetse fly.

Methods

In this study, RNA-Seq data of insect-stage trypanosomes were used to construct a T. brucei gene co-expression network using the weighted gene co-expression analysis (WGCNA) method. The study identified significant enriched modules for genes that play key roles during the parasite’s development in tsetse fly. Furthermore, potential 3′ untranslated region (UTR) regulatory elements for genes that clustered in the same module were identified using the Finding Informative Regulatory Elements (FIRE) tool.

Results

A fraction of gene modules (12 out of 27 modules) in the constructed network were found to be enriched in functional roles associated with the cell division, protein biosynthesis, mitochondrion, and cell surface. Additionally, 12 hub genes encoding proteins such as RNA-binding protein 6 (RBP6), arginine kinase 1 (AK1), brucei alanine-rich protein (BARP), among others, were identified for the 12 significantly enriched gene modules. In addition, the potential regulatory elements located in the 3′ untranslated regions of genes within the same module were predicted.

Conclusions

The constructed gene co-expression network provides a useful resource for network-based data mining to identify candidate genes for functional studies. This will enhance understanding of the molecular mechanisms that underlie important biological processes during parasite’s development in tsetse fly. Ultimately, these findings will be key in the identification of potential molecular targets for disease control.

Transcriptome-wide analysis of the Trypanosoma cruzi proliferative cycle identifies the periodically expressed mRNAs and their multiple levels of control

Trypanosoma cruzi is the protozoan parasite causing American trypanosomiasis or Chagas disease, a neglected parasitosis with important human health impact in Latin America. The efficacy of current therapy is limited, and its toxicity is high. Since parasite proliferation is a fundamental target for rational drug design, we sought to progress into its understanding by applying a genome-wide approach. Treating a TcI linage strain with hydroxyurea, we isolated epimastigotes in late G1, S and G2/M cell cycle stages at 70% purity. The sequencing of each phase identified 305 stage-specific transcripts (1.5-fold change, p≤0.01), coding for conserved cell cycle regulated proteins and numerous proteins whose cell cycle dependence has not been recognized before. Comparisons with the parasite T. brucei and the human host reveal important differences. The meta-analysis of T. cruzi transcriptomic and ribonomic data indicates that cell cycle regulated mRNAs are subject to sub-cellular compartmentalization. Compositional and structural biases of these genes- including CAI, GC content, UTR length, and polycistron position- may contribute to their regulation. To discover nucleotide motifs responsible for the co-regulation of cell cycle regulated genes, we looked for overrepresented motifs at their UTRs and found a variant of the cell cycle sequence motif at the 3' UTR of most of the S and G2 stage genes. We additionally identified hairpin structures at the 5' UTRs of a high proportion of the transcripts, suggesting that periodic gene expression might also rely on translation initiation in T. cruzi. In summary, we report a comprehensive list of T. cruzi cell cycle regulated genes, including many previously unstudied proteins, we show evidence favoring a multi-step control of their expression, and we identify mRNA motifs that may mediate their regulation. Our results provide novel information of the T. cruzi proliferative proteins and the integrated levels of their gene expression control.

Transcriptome-wide analysis of trypanosome mRNA decay reveals complex degradation kinetics and suggests a role for co-transcriptional degradation in determining mRNA levels

African trypanosomes are an excellent system for quantitative modelling of post-transcriptional mRNA control. Transcription is constitutive and polycistronic; individual mRNAs are excised by trans splicing and polyadenylation. We here measure mRNA decay kinetics in two life cycle stages, bloodstream and procyclic forms, by transcription inhibition and RNASeq. Messenger RNAs with short half-lives tend to show initial fast degradation, followed by a slower phase; they are often stabilized by depletion of the 5′–3′ exoribonuclease XRNA. Many longer-lived mRNAs show initial slow degradation followed by rapid destruction: we suggest that the slow phase reflects gradual deadenylation. Developmentally regulated mRNAs often show regulated decay, and switch their decay pattern. Rates of mRNA decay are good predictors of steady state levels for short mRNAs, but mRNAs longer than 3 kb show unexpectedly low abundances. Modelling shows that variations in splicing and polyadenylation rates can contribute to steady-state mRNA levels, but this is completely dependent on competition between processing and co-transcriptional mRNA precursor destruction.

The cell cycle regulated transcriptome of Trypanosoma brucei

Progression of the eukaryotic cell cycle requires the regulation of hundreds of genes to ensure that they are expressed at the required times. Integral to cell cycle progression in yeast and animal cells are temporally controlled, progressive waves of transcription mediated by cell cycle-regulated transcription factors. However, in the kinetoplastids, a group of early-branching eukaryotes including many important pathogens, transcriptional regulation is almost completely absent, raising questions about the extent of cell-cycle regulation in these organisms and the mechanisms whereby regulation is achieved. Here, we analyse gene expression over the Trypanosoma brucei cell cycle, measuring changes in mRNA abundance on a transcriptome-wide scale. We developed a “double-cut” elutriation procedure to select unperturbed, highly synchronous cell populations from log-phase cultures, and compared this to synchronization by starvation. Transcriptome profiling over the cell cycle revealed the regulation of at least 430 genes. While only a minority were homologous to known cell cycle regulated transcripts in yeast or human, their functions correlated with the cellular processes occurring at the time of peak expression. We searched for potential target sites of RNA-binding proteins in these transcripts, which might earmark them for selective degradation or stabilization. Over-represented sequence motifs were found in several co-regulated transcript groups and were conserved in other kinetoplastids. Furthermore, we found evidence for cell-cycle regulation of a flagellar protein regulon with a highly conserved sequence motif, bearing similarity to consensus PUF-protein binding motifs. RNA sequence motifs that are functional in cell-cycle regulation were more widespread than previously expected and conserved within kinetoplastids. These findings highlight the central importance of post-transcriptional regulation in the proliferation of parasitic kinetoplastids.

Cell biology of the trypanosome genome

Trypanosomes are a group of protozoan eukaryotes, many of which are major parasites of humans and livestock. The genomes of trypanosomes and their modes of gene expression differ in several important aspects from those of other eukaryotic model organisms. Protein-coding genes are organized in large directional gene clusters on a genome-wide scale, and their polycistronic transcription is not generally regulated at initiation. Transcripts from these polycistrons are processed by global trans-splicing of pre-mRNA. Furthermore, in African trypanosomes, some protein-coding genes are transcribed by a multifunctional RNA polymerase I from a specialized extranucleolar compartment. The primary DNA sequence of the trypanosome genomes and their cellular organization have usually been treated as separate entities. However, it is becoming increasingly clear that in order to understand how a genome functions in a living cell, we will need to unravel how the one-dimensional genomic sequence and its trans-acting factors are arranged in the three-dimensional space of the eukaryotic nucleus. Understanding this cell biology of the genome will be crucial if we are to elucidate the genetic control mechanisms of parasitism. Here, we integrate the concepts of nuclear architecture, deduced largely from studies of yeast and mammalian nuclei, with recent developments in our knowledge of the trypanosome genome, gene expression, and nuclear organization. We also compare this nuclear organization to those in other systems in order to shed light on the evolution of nuclear architecture in eukaryotes.

Trans-acting proteins regulating mRNA maturation, stability and translation in trypanosomatids

In trypanosomatids, alterations in gene expression in response to intrinsic or extrinsic signals are achieved through post-transcriptional mechanisms. In the last 20 years, research has concentrated on defining the responsible cis-elements in the untranslated regions of several regulated mRNAs. More recently, the focus has shifted towards the identification of RNA-binding proteins that act as trans-acting factors. Trypanosomatids have a large number of predicted RNA-binding proteins of which the vast majority have no orthologues in other eukaryotes. Several RNA-binding proteins have been shown to bind and/or regulate the expression of a group of mRNAs that code for functionally related proteins, indicating the possible presence of co-regulated mRNA cohorts.

Trypanosoma brucei PUF9 regulates mRNAs for proteins involved in replicative processes over the cell cycle

Many genes that are required at specific points in the cell cycle exhibit cell cycle–dependent expression. In the early-diverging model eukaryote and important human pathogen Trypanosoma brucei, regulation of gene expression in the cell cycle and other processes is almost entirely post-transcriptional. Here, we show that the T. brucei RNA-binding protein PUF9 stabilizes certain transcripts during S-phase. Target transcripts of PUF9—LIGKA, PNT1 and PNT2—were identified by affinity purification with TAP-tagged PUF9. RNAi against PUF9 caused an accumulation of cells in G2/M phase and unexpectedly destabilized the PUF9 target mRNAs, despite the fact that most known Puf-domain proteins promote degradation of their target mRNAs. The levels of the PUF9-regulated transcripts were cell cycle dependent, peaking in mid- to late- S-phase, and this effect was abolished when PUF9 was targeted by RNAi. The sequence UUGUACC was over-represented in the 3′ UTRs of PUF9 targets; a point mutation in this motif abolished PUF9-dependent stabilization of a reporter transcript carrying the PNT1 3′ UTR. LIGKA is involved in replication of the kinetoplast, and here we show that PNT1 is also kinetoplast-associated and its over-expression causes kinetoplast-related defects, while PNT2 is localized to the nucleus in G1 phase and redistributes to the mitotic spindle during mitosis. PUF9 targets may constitute a post-transcriptional regulon, encoding proteins involved in temporally coordinated replicative processes in early G2 phase.

The unicellular protozoan Trypanosoma brucei is the causative agent of African sleeping sickness, responsible for over 100,000 deaths annually, and is related to other important pathogens (e.g. Leishmania major and Trypanosoma cruzi). Unusually, these organisms do not regulate their genes by changing the rate at which they are copied into RNA, but by changing the rate of RNA destruction or the rate of translation into protein. We identified an RNA-binding protein, PUF9, responsible for the accumulation of several RNA molecules at a specific time point in the cell division cycle, just after DNA replication. Correspondingly, the proteins encoded by these RNAs appear to function in the division of various cellular structures at this time point or shortly afterwards. Two of them facilitate replication of the kinetoplast (an organelle containing the mitochondrial DNA) while another was found in the mitotic spindle. Their temporal co-expression may stem from another unusual feature of trypanosomes: only one copy of the kinetoplast (and several other organelles) are present per cell, their replication being coordinated with cell division. Indeed, PUF9 may be important in the control of organelle copy-number because suppression of PUF9 resulted in cells with too many kinetoplasts, flagella, or nuclei.

Identification of new kinetoplast DNA replication proteins in trypanosomatids based on predicted S-phase expression and mitochondrial targeting

Trypanosomatid parasites contain an unusual form of mitochondrial DNA (kinetoplast DNA [kDNA]) consisting of a catenated network of several thousand minicircles and a smaller number of maxicircles. Many of the proteins involved in the replication and division of kDNA are likely to have no counterparts in other organisms and would not be identified by similarity to known replication proteins in other organisms. A new kDNA replication protein conserved in kinetoplastids has been identified based on the presence of posttranscriptional regulatory sequences associated with S-phase gene expression and predicted mitochondrial targeting. The Leishmania major protein P105 (LmP105) and Trypanosoma brucei protein P93 (TbP93) localize to antipodal sites flanking the kDNA disk, where several other replication proteins and nascent minicircles have been localized. Like some of these kDNA replication proteins, the LmP105 protein is only present at the antipodal sites during S phase. RNA interference (RNAi) of TbP93 expression resulted in a cessation of cell growth and the loss of kDNA. Nicked/gapped forms of minicircles, the products of minicircle replication, were preferentially lost from the population of free minicircles during RNAi, suggesting involvement of TbP93 in minicircle replication. This approach should allow the identification of other novel proteins involved in the duplication of kDNA.

Evolution and pathology in chagas disease--a review

Trypanosoma cruzi acute infections often go unperceived, but one third of chronically infected individuals die of Chagas disease, showing diverse manifestations affecting the heart, intestines, and nervous systems. A common denominator of pathology in Chagas disease is the minimal rejection unit, whereby parasite-free target host cells are destroyed by immune system mononuclear effectors cells infiltrates. Another key feature stemming from T. cruzi infection is the integration of kDNA minicircles into the vertebrate host genome; horizontal transfer of the parasite DNA can undergo vertical transmission to the progeny of mammals and birds. kDNA integration-induced mutations can enter multiple loci in diverse chromosomes, generating new genes, pseudo genes and knock-outs, and resulting in genomic shuffling and remodeling over time. As a result of the juxtaposition of kDNA insertions with host open reading frames, novel chimeric products may be generated. Germ line transmission of kDNA-mutations determined the appearance of lesions in birds that are indistinguishable from those seen in Chagas disease patients. The production of tissue lesions showing typical minimal rejection units in birds' refractory to T. cruzi infection is consistent with the hypothesis that autoimmunity, likely triggered by integration-induced phenotypic alterations, plays a major role in the pathogenesis of Chagas disease.

Cell cycle-dependent localization and properties of a second mitochondrial DNA ligase in Crithidia fasciculata

The mitochondrial DNA in kinetoplastid protozoa is contained in a single highly condensed structure consisting of thousands of minicircles and approximately 25 maxicircles. The disk-shaped structure is termed kinetoplast DNA (kDNA) and is located in the mitochondrial matrix near the basal body. We have previously identified a mitochondrial DNA ligase (LIG kbeta) in the trypanosomatid Crithidia fasciculata that localizes to antipodal sites flanking the kDNA disk where several other replication proteins are localized. We describe here a second mitochondrial DNA ligase (LIG kalpha). LIG kalpha localizes to the kinetoplast primarily in cells that have completed mitosis and contain either a dividing kinetoplast or two newly divided kinetoplasts. Essentially all dividing or newly divided kinetoplasts show localization of LIG kalpha. The ligase is present on both faces of the kDNA disk and at a high level in the kinetoflagellar zone of the mitochondrial matrix. Cells containing a single nucleus show localization of the LIG kalpha to the kDNA but at a much lower frequency. The mRNA level of LIG kalpha varies during the cell cycle out of phase with that of LIG kbeta. LIG kalpha transcript levels are maximal during the phase when cells contain two nuclei, whereas LIG kbeta transcript levels are maximal during S phase. The LIG kalpha protein decays with a half-life of 100 min in the absence of protein synthesis. The periodic expression of the LIG kalpha transcript and the instability of the LIG kalpha protein suggest a possible role of the ligase in regulating minicircle replication.

Assigning functions to genes: identification of S-phase expressed genes in Leishmania major based on post-transcriptional control elements

Assigning functions to genes is one of the major challenges of the post-genomic era. Usually, functions are assigned based on similarity of the coding sequences to sequences of known genes, or by identification of transcriptional cis-regulatory elements that are known to be associated with specific pathways or conditions. In trypanosomatids, where regulation of gene expression takes place mainly at the post-transcriptional level, new approaches for function assignment are needed. Here we demonstrate the identification of novel S-phase expressed genes in Leishmania major, based on a post-transcriptional control element that was recognized in Crithidia fasciculata as involved in the cell cycle-dependent expression of several nuclear and mitochondrial S-phase expressed genes. Hypothesizing that a similar regulatory mechanism is manifested in L.major, we have applied a computational search for similar control elements in the genome of L.major. Our computational scan yielded 132 genes, of which 33% are homologues of known DNA metabolism genes and 63% lack any annotation. Experimental testing of seven of these genes revealed that their mRNAs cycle throughout the cell cycle, reaching a maximum level during S-phase or just prior to it. It is suggested that screening for post-transcriptional control elements associated with a specific function provides an efficient method for assigning functions to trypanosomatid genes.

A novel strategy to identify the location of necessary and sufficient cis-acting regulatory mRNA elements in trypanosomes

Expression of nearly all protein coding genes in trypanosomes is regulated post-transcriptionally, predominantly at the level of mRNA half-life. The identification of cis-acting elements involved in mRNA stability has been hindered by a lack of ability to screen for loss-of-regulation mutants. The method described in this article allows the region containing the necessary and sufficient elements within a mRNA to be identified and uses antibiotic resistance genes as both selectable markers and reporters. In the case of unstable mRNAs, the strategy can be extended by performing a screen for spontaneous loss-of-function mutants in regulatory parts of a mRNA. The method was validated by using the GPI-PLC mRNA, which is unstable in procyclic form trypanosomes and showed that the 3'UTR of the GPI-PLC mRNA contains all elements required for developmentally regulated instability. Loss-of-instability mutants all contained deletions within the 2300-nucleotide-long 3'UTR, and their analysis showed that a deletion including the last 800 nt of the gene stabilized the mRNA. The method is nonpresumptive, allows far more rapid screening for cis-elements than existing procedures, and has the advantage of identifying functional mutants. It is applicable to all eukaryotes using polycistronic transcription.

The 3'-untranslated region of cytochrome oxidase II mRNA functions in RNA editing of African trypanosomes exclusively as a cis guide RNA

RNA editing in trypanosomes is a post-transcriptional process responsible for correcting the coding sequences of many mitochondrial mRNAs. Uridines are specifically added or deleted from mRNA by an enzymatic cascade in which a pre-edited mRNA is specifically cleaved, uridines are added or removed, and the corrected mRNA is ligated. The process is directed by RNA molecules, termed guide RNAs (gRNA). The ability of this class of small, noncoding RNA to function in RNA editing is essential for these organisms. Typically, gRNAs are transcribed independent of the their cognate mRNA and anneal to form a binary RNA complex . An exception for this process may be cytochrome oxidase subunit II (COII) mRNA since a gene encoding a trans acting gRNA has not been identified. Using an in vitro editing assay we find that the 3' UTR of COII, indeed, functions as a guide for both the site and number of uridines added to the coding region of the COII mRNA. We further show that the guiding sequence within the COII 3' UTR can only function in COII editing when contiguous with the editing substrate, indicating that the 3' UTR of COII lacks sequence or structure information necessary to function as a trans-acting gRNA. While other RNAs have been shown to "guide" RNA processing reactions, our discovery that the COII 3' UTR directs editing of its cognate mRNA in cis, is a unique function for a 3' UTR. The findings described here have led us to propose a new model for the evolution of gRNAs in kinetoplastids.

Mitochondrial histone-like DNA-binding proteins are essential for normal cell growth and mitochondrial function in Crithidia fasciculata

The Crithidia fasciculata KAP2 and KAP3 proteins are closely related kinetoplast-specific histone-like DNA-binding proteins. The KAP2 and KAP3 genes are 46% identical and are arranged in tandem on the chromosomal DNA. Disruption of both alleles of either gene alone shows no detectable phenotype. However, replacement of both copies of the sequence encoding the entire KAP2 and KAP3 locus increases maxicircle mRNA levels two- to fourfold. These double-knockout cells are viable but grow extremely slowly, have reduced respiration and very abnormal cell morphologies, and accumulate numerous large vacuoles. The extreme phenotype of these mutant cells suggests an important role for the KAP2 and KAP3 proteins in mitochondrial metabolism and cell growth.