Maxicircle (mitochondrial) genome sequence (partial) of Leishmania major: gene content, arrangement and composition compared with Leishmania tarentolae

We report 8420 bp of DNA sequence data from the maxicircle (mitochondrial) genome of Leishmania major (MHOM/SU/73/5ASKH), a much larger portion of this genome than has been reported previously from any Leishmania species infecting humans. This region contains 10 partial and complete genes: 5 protein-encoding genes (COII, COIII, ND1, ND7 and Cyt b); two ribosomal RNA subunits (12S and 9S) and three unidentified open reading frames (MURF1, MURF4 (ATPase6) and MURF5), as in the lizard-infecting species L. tarentolae. The genes from L. major exhibit 85-87% identity with those of L. tarentolae at the nucleotide level and 71-94% identity at the amino acid level. Most differences between sequences from the two species are transversions. The gene order and arrangement within the maxicircle of L. major are similar to those in L. tarentolae, but base composition and codon usage differ between the species. Codons assigned for initiation for protein-coding genes available for comparison are similar in five genes in the two species. Pre-editing was identified in some of the protein-coding genes. Short intergenic non-coding regions are also present in L. major as they are in L. tarentolae. Intergenic regions between 9S rRNA and MURF5, MURF1 and ND1 genes are G+C rich and considered to be extensive RNA editing regions. The RNA editing process is likely to be conserved in similar pattern in L. major as in L. tarentolae.     

Comparison of the maxicircle (mitochondrial) genomes of Leishmania tarentolae and Trypanosoma brucei at the level of nucleotide sequence

The entire 16.7-kilobase (kb) transcribed region of the Leishmania tarentolae maxicircle was compared to the entire 15-kb transcribed region of the Trypanosoma brucei maxicircle at the nucleotide sequence level by dot matrix analysis and by alignments of individual genes. The L. tarentolae NADH dehydrogenase subunit 1 (ND1) gene was identified in a newly obtained 2.9-kb sequence. All but two regions which flank the cytochrome b gene are highly conserved in both species. One 3.1-kb region in L. tarentolae that contains the cytochrome oxidase subunit III (COIII) gene and several open reading frames corresponds to a 2-kb sequence in T. brucei with limited sequence homology that lacks the COIII gene. Another 0.6-kb region that comprises an unidentified open reading frame (open reading frame 12) in L. tarentolae is substituted by a nonhomologous 0.4-kb open reading frame in T. brucei. A short intergenic region between the ND1 gene and the maxicircle unidentified reading frame 1 gene shows limited sequence homology, and the regions between the ND4 and ND5 genes and between the COI and ND4 genes are not conserved. All of the intergenic regions share G + C richness and a similar pattern of G versus C strand bias. 1.8 kb of the L. tarentolae divergent region (DV) and around 3 kb of the T. brucei DV were also obtained. The T. brucei DV sequences were not homologous to the L. tarentolae DV sequence but were organized in a similar fashion with tandem repeats of varying complexity.     

Genome sequencing of the lizard parasite Leishmania tarentolae reveals loss of genes associated to the intracellular stage of human pathogenic species

The Leishmania tarentolae Parrot-TarII strain genome sequence was resolved to an average 16-fold mean coverage by next-generation DNA sequencing technologies. This is the first non-pathogenic to humans kinetoplastid protozoan genome to be described thus providing an opportunity for comparison with the completed genomes of pathogenic Leishmania species. A high synteny was observed between all sequenced Leishmania species. A limited number of chromosomal regions diverged between L. tarentolae and L. infantum, while remaining syntenic to L. major. Globally, >90% of the L. tarentolae gene content was shared with the other Leishmania species. We identified 95 predicted coding sequences unique to L. tarentolae and 250 genes that were absent from L. tarentolae. Interestingly, many of the latter genes were expressed in the intracellular amastigote stage of pathogenic species. In addition, genes coding for products involved in antioxidant defence or participating in vesicular-mediated protein transport were underrepresented in L. tarentolae. In contrast to other Leishmania genomes, two gene families were expanded in L. tarentolae, namely the zinc metallo-peptidase surface glycoprotein GP63 and the promastigote surface antigen PSA31C. Overall, L. tarentolae's gene content appears better adapted to the promastigote insect stage rather than the amastigote mammalian stage.     

Direct visualisation of RNA editing within a Leishmania tarentolae mitochondrial extract

The coding sequence within several mitochondrial mRNAs of the trypanosomatid protozoa is created through editing by the precise insertion and deletion of U nucleotides. The biochemical characterisation of the editing reaction in the Leishmania genus of the trypanosomatids has been hindered by the lack of a direct in vitro assay. We describe here the first direct assay for the detection of guide RNA-directed editing mediated by a mitochondrial extract prepared from two independent isolates of Leishmania tarentolae. The assay enabled the editing activity within a L. tarentolae mitochondrial extract to be significantly enriched and will facilitate the characterisation of the editing reaction. The results suggest that the difficulty in establishing an assay for the L. tarentolae reaction was not simply a result of the catalytic machinery being limiting but rather reflected the presence of constraints on both the guide RNA and mRNA sequences.     

Identification of a gp63 surface glycoprotein in Leishmania tarentolae

The promastigote stage of most if not all Leishmania species possesses an abundant surface glycoprotein of 63 kDa (gp63) that has protease activity. We show that the lizard parasite Leishmania tarentolae appears to lack the surface protease activity. L. tarentolae does, however, possess an approximately 63-kDa molecule that is antigenically cross-reactive with the L. major gp63. Additionally, the genome of L. tarentolae contains sequences that hybridise at high stringency to a L. major gp63 gene probe.     

Sequences of six genes and several open reading frames in the kinetoplast maxicircle DNA of Leishmania tarentolae

The DNA sequence of approximately 80% of the transcribed region of the kinetoplast maxicircle DNA of Leishmania tarentolae was obtained, and structural genes were localized by comparison of the translated amino acid sequences with those of known mitochondrial genes from other organisms. By this method, the genes for cytochrome oxidase subunits I, II, and III, cytochrome b, and human mitochondrial unidentified reading frames 4 and 5 were identified. By comparing the amino acid sequences of the putative L. tarentolae genes with those of known genes, we conclude that TGA codes for tryptophan, as in most other mitochondrial systems. This is the only apparent change from the universal genetic code. The six identified structural genes show various degrees of divergence from the homologous genes in other species, with cytochrome oxidase subunit I being the most conserved and cytochrome oxidase subunit III being the least conserved. A comparison of the cytochrome b genes from L. tarentolae and Trypanosoma brucei showed that the ratio of transversions to transitions is 1:1, suggesting that these species diverged from each other more than 80 X 10(6) years ago. Several as yet unidentified open reading frames were also present in the maxicircle sequence. These data confirm that maxicircle DNA has a coding potential which typifies other mitochondrial systems.     

In vitro RNA editing-like activity in a mitochondrial extract from Leishmania tarentolae.

A mitochondrial extract from Leishmania tarentolae directs the incorporation of uridylate (U) residues within the pre-edited domain of synthetic cytochrome b (CYb) and NADH dehydrogenase subunit 7 mRNA. This has several characteristics of an in vitro RNA editing activity, but no direct evidence for involvement of guide RNAs was obtained. Inhibition by micrococcal nuclease suggests a requirement for some type of endogenous RNA. The limitation of internal U-incorporation to the pre-edited region in the CYb mRNA and the inhibition by deletion or substitution of both mRNA anchor sequences for CYb gRNA-I and -II could be consistent either with a gRNA-mediated process or a secondary structure-mediated process. A low level of incorporation of [alpha-32P]CTP occurs at the same sites as UTP. Internal U-incorporation activity is selectively inhibited by heterologous RNAs, suggesting an involvement of low affinity RNA-binding proteins which can be competed by the added RNA.     

Editing of the grapevine mitochondrial cytochrome b mRNA and molecular modeling of the protein

Cytochrome b (COB), the central catalytic subunit of ubiquinol cytochrome c reductase, is a component of the transmembrane electron transfer chain that generates proton motive force. Some plant COB mRNAs are processed by RNA editing, which changes the gene coding sequence. This report presents the sequences of the grapevine (Vitis vinifera L.) mitochondrial gene for apocytochrome b (cob), the edited mRNA and the deduced protein. Grapevine COB is 393 amino acids long and is 98% identical to homologs in rapeseed, Arabidopsis thaliana and Oenothera sp. Twenty-one C-U editing sites were identified in the grapevine cob mRNA, resulting in 20 amino acid changes. These changes increase the overall hydrophobicity of the protein and result in a more conserved protein. Molecular modeling of grapevine COB shows that residues changed by RNA editing fit the secondary structure characteristic of an integral membrane protein. This is the first complete mitochondrial gene reported for grapevine. Novel RNA editing sites were identified in grapevine cob, which have not been previously reported for other plants.     

Sequence and structural requirements for optimal guide RNA-directed insertional editing within Leishmania tarentolae

The coding sequence of several mitochondrial mRNAs of the trypanosomatid family of protozoa is created by the guide RNA-directed insertion and deletion of uridylates (Us). Selection-amplification was used to explore the sequence and structure of the guide RNA and mRNA required for efficient insertional editing within a mitochondrial extract prepared from Leishmania tarentolae. This study identifies several novel features of the editing reaction in addition to several that are consistent with the previous mutagenesis and phylogenetic analysis of the reaction in Trypanosoma brucei, a distantly related trypanosomatid. Specifically, there is a strong bias against cytidines 5' of the editing sites and guanosines immediately 3' of guiding nucleotides. U insertions are directed both 5' and 3' of a genomically encoded U, which was previously assumed not to occur. Base pairing immediately flanking an editing site can significantly stimulate the editing reaction and affect the reaction fidelity but is not essential. Likewise, single-stranded RNA in the region upstream of the editing site, not necessarily immediately adjacent, can facilitate editing but is also not essential. The editing of an RNA containing many of the optimal features is linear with increasing quantities of extract permitting specific activity measurements to be made that are not possible with previously described T. brucei and L. tarentolae assays. The reaction catalyzed by the L. tarentolae extract can be highly accurate, which does not support a proposed model for editing that was based largely on the inaccuracy of an earlier in vitro reaction.     

Trypanosoma brucei ATPase subunit 6 mRNA bound to gA6-14 forms a conserved three-helical structure

T. brucei survival relies on the expression of mitochondrial genes, most of which require RNA editing to become translatable. In trypanosomes, RNA editing involves the insertion and deletion of uridylates, a developmentally regulated process directed by guide RNAs (gRNAs) and catalyzed by the editosome, a complex of proteins. The pathway for mRNA/gRNA complex formation and assembly with the editosome is still unknown. Work from our laboratory has suggested that distinct mRNA/gRNA complexes anneal to form a conserved core structure that may be important for editosome assembly. The secondary structure for the apocytochrome b (CYb) pair has been previously determined and is consistant with our model of a three-helical structure. Here, we used cross-linking and solution structure probing experiments to determine the structure of the ATPase subunit 6 (A6) mRNA hybridized to its cognate gA6-14 gRNA in different stages of editing. Our results indicate that both unedited and partially edited A6/gA6-14 pairs fold into a three-helical structure similar to the previously characterized CYb/gCYb-558 pair. These results lead us to conclude that at least two mRNA/gRNA pairs with distinct editing sites and distinct primary sequences fold to a three-helical secondary configuration that persists through the first few editing events.     

Trypanosome mitochondrial 3' terminal uridylyl transferase (TUTase): the key enzyme in U-insertion/deletion RNA editing

A 3' terminal RNA uridylyltransferase was purified from mitochondria of Leishmania tarentolae and the gene cloned and expressed from this species and from Trypanosoma brucei. The enzyme is specific for 3' U-addition in the presence of Mg(2+). TUTase is present in vivo in at least two stable configurations: one contains a approximately 500 kDa TUTase oligomer and the other a approximately 700 kDa TUTase complex. Anti-TUTase antiserum specifically coprecipitates a small portion of the p45 and p50 RNA ligases and approximately 40% of the guide RNAs. Inhibition of TUTase expression in procyclic T. brucei by RNAi downregulates RNA editing and appears to affect parasite viability.