A survey of Leishmania braziliensis genome by shotgun sequencing

LeishDB: a database of coding gene annotation and non-coding RNAs in Leishmania braziliensis

Database URL

www.leishdb.com.

The putative Leishmania telomerase RNA (LeishTER) undergoes trans-splicing and contains a conserved template sequence

Telomerase RNAs (TERs) are highly divergent between species, varying in size and sequence composition. Here, we identify a candidate for the telomerase RNA component of Leishmania genus, which includes species that cause leishmaniasis, a neglected tropical disease. Merging a thorough computational screening combined with RNA-seq evidence, we mapped a non-coding RNA gene localized in a syntenic locus on chromosome 25 of five Leishmania species that shares partial synteny with both Trypanosoma brucei TER locus and a putative TER candidate-containing locus of Crithidia fasciculata. Using target-driven molecular biology approaches, we detected a ∼2,100 nt transcript (LeishTER) that contains a 5′ spliced leader (SL) cap, a putative 3′ polyA tail and a predicted C/D box snoRNA domain. LeishTER is expressed at similar levels in the logarithmic and stationary growth phases of promastigote forms. A 5′SL capped LeishTER co-immunoprecipitated and co-localized with the telomerase protein component (TERT) in a cell cycle-dependent manner. Prediction of its secondary structure strongly suggests the existence of a bona fide single-stranded template sequence and a conserved C[U/C]GUCA motif-containing helix II, representing the template boundary element. This study paves the way for further investigations on the biogenesis of parasite TERT ribonucleoproteins (RNPs) and its role in parasite telomere biology.

The Comparative Genomics and Phylogenomics of Leishmania amazonensis Parasite

Leishmaniasis is an infectious disease caused by Leishmania species. Leishmania amazonensis is a New World Leishmania species belonging to the Mexicana complex, which is able to cause all types of leishmaniasis infections. The L. amazonensis reference strain MHOM/BR/1973/M2269 was sequenced identifying 8,802 codifying sequences (CDS), most of them of hypothetical function. Comparative analysis using six Leishmania species showed a core set of 7,016 orthologs. L. amazonensis and Leishmania mexicana share the largest number of distinct orthologs, while Leishmania braziliensis presented the largest number of inparalogs. Additionally, phylogenomic analysis confirmed the taxonomic position for L. amazonensis within the “Mexicana complex”, reinforcing understanding of the split of New and Old World Leishmania. Potential non-homologous isofunctional enzymes (NISE) were identified between L. amazonensis and Homo sapiens that could provide new drug targets for development.

Identification of putative potassium channel homologues in pathogenic protozoa

K⁺ channels play a vital homeostatic role in cells and abnormal activity of these channels can dramatically alter cell function and survival, suggesting that they might be attractive drug targets in pathogenic organisms. Pathogenic protozoa lead to diseases such as malaria, leishmaniasis, trypanosomiasis and dysentery that are responsible for millions of deaths each year worldwide. The genomes of many protozoan parasites have recently been sequenced, allowing rational design of targeted therapies. We analyzed the genomes of pathogenic protozoa and show the existence within them of genes encoding putative homologues of K⁺ channels. These protozoan K⁺ channel homologues represent novel targets for anti-parasitic drugs. Differences in the sequences and diversity of human and parasite proteins may allow pathogen-specific targeting of these K⁺ channel homologues.

Trypanosomatid comparative genomics: Contributions to the study of parasite biology and different parasitic diseases

In 2005, draft sequences of the genomes of Trypanosoma brucei, Trypanosoma cruzi and Leishmania major, also known as the Tri-Tryp genomes, were published. These protozoan parasites are the causative agents of three distinct insect-borne diseases, namely sleeping sickness, Chagas disease and leishmaniasis, all with a worldwide distribution. Despite the large estimated evolutionary distance among them, a conserved core of ~6,200 trypanosomatid genes was found among the Tri-Tryp genomes. Extensive analysis of these genomic sequences has greatly increased our understanding of the biology of these parasites and their host-parasite interactions. In this article, we review the recent advances in the comparative genomics of these three species. This analysis also includes data on additional sequences derived from other trypanosmatid species, as well as recent data on gene expression and functional genomics. In addition to facilitating the identification of key parasite molecules that may provide a better understanding of these complex diseases, genome studies offer a rich source of new information that can be used to define potential new drug targets and vaccine candidates for controlling these parasitic infections.

Identification of intracellular and plasma membrane calcium channel homologues in pathogenic parasites

Ca²⁺ channels regulate many crucial processes within cells and their abnormal activity can be damaging to cell survival, suggesting that they might represent attractive therapeutic targets in pathogenic organisms. Parasitic diseases such as malaria, leishmaniasis, trypanosomiasis and schistosomiasis are responsible for millions of deaths each year worldwide. The genomes of many pathogenic parasites have recently been sequenced, opening the way for rational design of targeted therapies. We analyzed genomes of pathogenic protozoan parasites as well as the genome of Schistosoma mansoni, and show the existence within them of genes encoding homologues of mammalian intracellular Ca²⁺ release channels: inositol 1,4,5-trisphosphate receptors (IP₃Rs), ryanodine receptors (RyRs), two-pore Ca²⁺ channels (TPCs) and intracellular transient receptor potential (Trp) channels. The genomes of Trypanosoma, Leishmania and S. mansoni parasites encode IP₃R/RyR and Trp channel homologues, and that of S. mansoni additionally encodes a TPC homologue. In contrast, apicomplexan parasites lack genes encoding IP₃R/RyR homologues and possess only genes encoding TPC and Trp channel homologues (Toxoplasma gondii) or Trp channel homologues alone. The genomes of parasites also encode homologues of mammalian Ca²⁺influx channels, including voltage-gated Ca²⁺ channels and plasma membrane Trp channels. The genome of S. mansoni also encodes Orai Ca²⁺ channel and STIM Ca²⁺ sensor homologues, suggesting that store-operated Ca²⁺ entry may occur in this parasite. Many anti-parasitic agents alter parasite Ca²⁺ homeostasis and some are known modulators of mammalian Ca²⁺ channels, suggesting that parasite Ca²⁺ channel homologues might be the targets of some current anti-parasitic drugs. Differences between human and parasite Ca²⁺ channels suggest that pathogen-specific targeting of these channels may be an attractive therapeutic prospect.

Evidence of RNA editing in Leishmania braziliensis promastigotes

RNA editing in trypanosomatids is an elaborate form of post-transcriptional processing that inserts and deletes uridines in many mitochondrial pre-mRNAs, providing the genetic information needed to create functional transcripts. The process has been extensively analyzed in Trypanosoma brucei, Crithidia fasciculata, and Leishmania tarentolae. However, few data exist on this mechanism in pathogenic Leishmania species. Here, we show evidence that this process also operates in Leishmania braziliensis, being the first time that RNA editing has been described in a species of the Viannia subgenus. A partially edited transcript corresponding to the NADH dehydrogenase subunit 8 (ND8) gene was identified in L. braziliensis promastigotes. Sequence analysis allowed the identification of the maxicircle-encoded cryptogene, which shows a high degree of sequence conservation with the corresponding cryptogenes in other Leishmania species. Although an edition pattern could be postulated for the ND8 transcripts in L. braziliensis, attempts to isolate completely edited transcripts by RT-PCR were not fruitful; instead, many transcripts with partial and unexpected editing patterns were isolated. This data, together with our inability to detect full-size transcripts by Northern blotting in promastigotes of L. braziliensis, led us to the suggestion that the strain used in this study (M2904) lacks of critical RNA guides for a complete edition of ND8 transcripts.

ReRep: computational detection of repetitive sequences in genome survey sequences (GSS)

Background

Genome survey sequences (GSS) offer a preliminary global view of a genome since, unlike ESTs, they cover coding as well as non-coding DNA and include repetitive regions of the genome. A more precise estimation of the nature, quantity and variability of repetitive sequences very early in a genome sequencing project is of considerable importance, as such data strongly influence the estimation of genome coverage, library quality and progress in scaffold construction. Also, the elimination of repetitive sequences from the initial assembly process is important to avoid errors and unnecessary complexity. Repetitive sequences are also of interest in a variety of other studies, for instance as molecular markers.

Results

We designed and implemented a straightforward pipeline called ReRep, which combines bioinformatics tools for identifying repetitive structures in a GSS dataset. In a case study, we first applied the pipeline to a set of 970 GSSs, sequenced in our laboratory from the human pathogen Leishmania braziliensis, the causative agent of leishmaniosis, an important public health problem in Brazil. We also verified the applicability of ReRep to new sequencing technologies using a set of 454-reads of an Escheria coli. The behaviour of several parameters in the algorithm is evaluated and suggestions are made for tuning of the analysis.

Conclusion

The ReRep approach for identification of repetitive elements in GSS datasets proved to be straightforward and efficient. Several potential repetitive sequences were found in a L. braziliensis GSS dataset generated in our laboratory, and further validated by the analysis of a more complete genomic dataset from the EMBL and Sanger Centre databases. ReRep also identified most of the E. coli K12 repeats prior to assembly in an example dataset obtained by automated sequencing using 454 technology. The parameters controlling the algorithm behaved consistently and may be tuned to the properties of the dataset, in particular to the length of sequencing reads and the genome coverage. ReRep is freely available for academic use at .

Proteome analysis of Leishmania (Viannia) braziliensis by two-dimensional gel electrophoresis and mass spectrometry

Leishmania (Viannia) braziliensis, a protozoan parasite widespread in the New World, is responsible for the infection of different mammal orders, including humans. This species is considered to be a major etiological agent of American cutaneous leishmaniasis. A proteomic study was carried out to identify proteins expressed by L. (V.) braziliensis. One hundred and one spots representing 75 protein entries were identified by MALDI-TOF-TOF. Isoelectric point values estimated by gel electrophoresis matched closely with predicted values, although some discrepancies existed suggesting that post-translational protein modifications may be common in L. braziliensis. Moreover, 20 hypothetical proteins were experimentally identified. Identified proteins were classified into 15 groups according to biological process. Among the proteins identified, approximately 40% have not been previously reported in a proteomic map of Leishmania. In addition, a number of potential virulence factors and drug targets were identified in this protein map, including some proteins associated with the metastatic phenotype. This study describes the first compilation of a proteomic reference map for L. braziliensis (pI 4-7, M(r) 10-130 kDa) and provides a very useful tool for comparative studies of strains isolated from patients presenting different clinical manifestations of leishmaniasis as well as a potential tool to identify markers for clinical diagnosis, therapeutics, and prognosis.

Comparative genomics: from genotype to disease phenotype in the leishmaniases

Recent progress in sequencing the genomes of several Leishmania species, causative agents of cutaneous, mucocutaneous and visceral leishmaniasis, is revealing unusual features of potential relevance to parasite virulence and pathogenesis in the host. While the genomes of Leishmania major, Leishmania braziliensis and Leishmania infantum are highly similar in content and organisation, species-specific genes and mechanisms distinguish one from another. In particular, the presence of retrotransposons and the components of a putative RNA interference machinery in L. braziliensis suggest the potential for both greater diversity and more tractable experimentation in this Leishmania Viannia species.

Comparative genomic analysis of three Leishmania species that cause diverse human disease

Leishmania parasites cause a broad spectrum of clinical disease. Here we report the sequencing of the genomes of two species of Leishmania: Leishmania infantum and Leishmania braziliensis. The comparison of these sequences with the published genome of Leishmania major reveals marked conservation of synteny and identifies only approximately 200 genes with a differential distribution between the three species. L. braziliensis, contrary to Leishmania species examined so far, possesses components of a putative RNA-mediated interference pathway, telomere-associated transposable elements and spliced leader-associated SLACS retrotransposons. We show that pseudogene formation and gene loss are the principal forces shaping the different genomes. Genes that are differentially distributed between the species encode proteins implicated in host-pathogen interactions and parasite survival in the macrophage.

The use of Tn5 transposable elements in a gene trapping strategy for the protozoan Leishmania

The use of transposable elements as a gene-trapping strategy is a powerful tool for gene discovery. Herein we describe the development of a transposable system, based on the bacterial Tn5 transposon, which has been used successfully in Leishmania braziliensis. The transposon carries the neomycin phosphotransferase gene, which is expressed only when inserted in-frame with a Leishmania gene present in the target DNA. Four cosmid clones from a L. braziliensis genomic library were used as targets in transposition reactions and four insertional libraries were constructed and transfected in L. braziliensis. Clones resistant to G418 were selected and analysed by immunofluorescence in order to identify the subcellular localisation of the protein coded by the trapped gene. A definitive subcellular localisation for neomycin phosphotransferase/targeted protein fusion was not obtained in any of the four Leishmania clones investigated. However, the constructed transposable element is highly efficient considering the frequency of insertion in large targets and is therefore a useful tool for functional genetic studies in Leishmania. Our data confirm the utility of the Tn5 transposon system for insertion of sequencing priming sites into target DNA. Furthermore, the high frequency of insertion and even distribution are important in studying genomic regions bearing long and polymorphic repetitive sequences.

Organization of H locus conserved repeats in Leishmania (Viannia) braziliensis correlates with lack of gene amplification and drug resistance

Resistance to antimonials is a major problem when treating visceral leishmaniasis in India and has already been described for New World parasites. Clinical response to meglumine antimoniate in patients infected with parasites of the Viannia sub-genus can be widely variable, suggesting the presence of mechanisms of drug resistance. In this work, we have compared L. major and L. braziliensis mutants selected in different drugs. The cross-resistance profiles of some cell lines resembled those of mutants bearing H locus amplicons. However, amplified episomal molecules were exclusively detected in L. major mutants. The analysis of the L. braziliensis H region revealed a strong conservation of gene synteny. The typical intergenic repeats that are believed to mediate the amplification of the H locus in species of the Leishmania sub-genus are partially conserved in the Viannia species. The conservation of these non-coding elements in equivalent positions in both species is indicative of their relevance within this locus. The absence of amplicons in L. braziliensis suggests that this species may not favour extra-chromosomal gene amplification as a source of phenotypic heterogeneity and fitness maintenance in changing environments.

Genomic organization and expression of the HSP70 locus in New and Old World Leishmania species

Heat shock is believed to be a developmental inductor of differentiation in Leishmania. Furthermore, heat shock genes are extensively studied as gene models to decipher mechanisms of gene regulation in kinetoplastids. Here, we describe the organization and expression of the HSP70 loci in representative Leishmania species (L. infantum, L. major, L. tropica, L. mexicana, L. amazonensis and L. braziliensis). With the exception of L. braziliensis, the organization of the HSP70 loci was found to be well conserved among the other Leishmania species. Two types of genes, HSP70-I and HSP70-II, were found to be present in these Leishmania species except for L. braziliensis that lacks HSP70-II gene. Polymorphisms in the HSP70 locus allow the differentiation of the Old and New World species within the subgenus Leishmania. A notable discrepancy between our data and those of the L. major genome database in relation to the gene copy number composing the L. major HSP70 locus was revealed. The temperature-dependent accumulation of the HSP70-I mRNAs is also conserved among the different Leishmania species with the exception of L. braziliensis. In spite of these differences, analysis of the HSP70 synthesis indicated that the HSP70 mRNAs are also preferentially translated during heat shock in L. braziliensis.

Characterization of LST-R533: uncovering a novel repetitive element in Leishmania

We have previously isolated and sequenced a novel repetitive element, now named LST-R533, which is present in four different regions of one extremity of Leishmania major chromosome 20. The repeats are polymorphic in size, ranging from 367 to 533 bp and contain an internal 81 bp sequence with highly conserved segments (14-81 bp long) dispersed throughout the parasite's genome. These sequences were not found in coding regions of any predicted gene in L. major Friedlin genome, but are part of untranslated regions of some Leishmania transcripts. Analysis of the 81 bp sequence revealed significant degrees of identity with retrotransposons described in several other organisms. The presence of the sequence in other species from genus Leishmania was determined by Southern hybridisation and DNA sequencing. This analysis indicated the conservation of the 81-nucleotide element in all the Leishmania species evaluated. No sequences corresponding to LST-R533 or the 81 bp element were found on either Trypanosoma brucei or Trypanosoma cruzi databanks.

Protozoan genomes: gene identification and annotation

The draft sequence of several complete protozoan genomes is now available and genome projects are ongoing for a number of other species. Different strategies are being implemented to identify and annotate protein coding and RNA genes in these genomes, as well as study their genomic architecture. Since the genomes vary greatly in size, GC-content, nucleotide composition, and degree of repetitiveness, genome structure is often a factor in choosing the methodology utilised for annotation. In addition, the approach taken is dictated, to a greater or lesser extent, by the particular reasons for carrying out genome-wide analyses and the level of funding available for projects. Nevertheless, these projects have provided a plethora of material that will aid in understanding the biology and evolution of these parasites, as well as identifying new targets that can be used to design urgently required drug treatments for the diseases they cause.