The Complete Chromosomal Organization of the Reference Strain of the Leishmania Genome Project, L. major `Friedlin'

Making the Most of Its Short Reads: A Bioinformatics Workflow for Analysing the Short-Read-Only Data of Leishmania orientalis (Formerly Named Leishmania siamensis) Isolate PCM2 in Thailand

Simple Summary

Leishmaniasis is a parasitic disease caused by flagellated protozoa of the genus Leishmania. Multiple genome sequencing platforms have been employed to complete Leishmania genomes at the expense of high cost. This study proposes an integrative bioinformatic workflow for assembling only the short-read data of Leishmania orientalis isolate PCM2 from Thailand and produce an acceptable-quality genome for further genomic analysis. This workflow gives extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Abstract

Background: Leishmania orientalis (formerly named Leishmania siamensis) has been neglected for years in Thailand. The genomic study of L. orientalis has gained much attention recently after the release of the first high-quality reference genome of the isolate LSCM4. The integrative approach of multiple sequencing platforms for whole-genome sequencing has proven effective at the expense of considerably expensive costs. This study presents a preliminary bioinformatic workflow including the use of multi-step de novo assembly coupled with the reference-based assembly method to produce high-quality genomic drafts from the short-read Illumina sequence data of L. orientalis isolate PCM2. Results: The integrating multi-step de novo assembly by MEGAHIT and SPAdes with the reference-based method using the L. enriettii genome and salvaging the unmapped reads resulted in the 30.27 Mb genomic draft of L. orientalis isolate PCM2 with 3367 contigs and 8887 predicted genes. The results from the integrated approach showed the best integrity, coverage, and contig alignment when compared to the genome of L. orientalis isolate LSCM4 collected from the northern province of Thailand. Similar patterns of gene ratios and frequency were observed from the GO biological process annotation. Fifty GO terms were assigned to the assembled genomes, and 23 of these (accounting for 61.6% of the annotated genes) showed higher gene counts and ratios when results from our workflow were compared to those of the LSCM4 isolate. Conclusions: These results indicated that our proposed bioinformatic workflow produced an acceptable-quality genome of L. orientalis strain PCM2 for functional genomic analysis, maximising the usage of the short-read data. This workflow would give extensive information required for identifying strain-specific markers and virulence-associated genes useful for drug and vaccine development before a more exhaustive and expensive investigation.

Phenotypical and genotypical differences among Leishmania (Leishmania) amazonensis isolates that caused different clinical frames in humans and dogs: A systematic review

Leishmania (Leishmania) amazonensis is an important etiological agent of American cutaneous leishmaniasis (ACL) in Brazil. The species causes a large spectrum of clinical manifestations in humans and dogs, ranging from cutaneous, cutaneous diffuse, mucocutaneous, and visceral involvement, however, the factors that drive the development of different disease forms by the same species are not yet fully known. In the present work, it was systematically reviewed the studies addressing phenotypic and genotypic characteristics of Leishmania (L.) amazonensis isolates causing cutaneous and visceral clinical frames in humans and dogs, comparing the results observed. For this, four research databases were searched for the following keywords: (Leishmania amazonensis AND visceral leishmaniasis) AND (tropism OR virulence OR visceralization OR adaptations OR mutation OR clinical presentation OR resistance OR survival OR wide spectrum). The results revealed that the complexity disease seems to involve the combination of genetic factors of the parasite (as modifications in molecules related to the virulence and metabolism) and also of the host's immune background and status. Nonetheless, the exact mechanism that leads to different clinical manifestations between strains of the same species is still uncertain and future studies must be developed to better elucidate this phenomenon.

Comparative genomics of Tunisian Leishmania major isolates causing human cutaneous leishmaniasis with contrasting clinical severity

Zoonotic cutaneous leishmaniasis caused by Leishmania (L.) major parasites affects urban and suburban areas in the center and south of Tunisia where the disease is endemo-epidemic. Several cases were reported in human patients for which infection due to L. major induced lesions with a broad range of severity. However, very little is known about the mechanisms underlying this diversity. Our hypothesis is that parasite genomic variability could, in addition to the host immunological background, contribute to the intra-species clinical variability observed in patients and explain the lesion size differences observed in the experimental model. Based on several epidemiological, in vivo and in vitro experiments, we focused on two clinical isolates showing contrasted severity in patients and BALB/c experimental mice model. We used DNA-seq as a high-throughput technology to facilitate the identification of genetic variants with discriminating potential between both isolates. Our results demonstrate that various levels of heterogeneity could be found between both L. major isolates in terms of chromosome or gene copy number variation (CNV), and that the intra-species divergence could surprisingly be related to single nucleotide polymorphisms (SNPs) and Insertion/Deletion (InDels) events. Interestingly, we particularly focused here on genes affected by both types of variants and correlated them with the observed gene CNV. Whether these differences are sufficient to explain the severity in patients is obviously still open to debate, but we do believe that additional layers of -omic information is needed to complement the genomic screen in order to draw a more complete map of severity determinants.

Parasexuality and mosaic aneuploidy in Leishmania: alternative genetics

Reproduction as identical or similar organisms in most biological systems depends on the extreme accuracy of the mitotic (and meiotic) mechanisms involved in the transmission of the genetic material to the two daughter cells. Character recombination and genotype diversification are ensured by the alternation between haploidy and diploidy, which corresponds to the most predominant model in sexually reproducing organisms. In Leishmania, the unique association of high levels of automixis and of constitutive 'mosaic aneuploidy' unexpectedly does not lead to loss of heterozygosity but constitutes an alternative for genotype recombination, hence a source of adaptability.

Constitutive mosaic aneuploidy is a unique genetic feature widespread in the Leishmania genus

Using fluorescence in situ hybridization, we determined the ploidy of four species of Leishmania: Leishmania infantum, Leishmania donovani, Leishmania tropica and Leishmania amazonensis. We found that each cell in a strain possesses a combination of mono-, di- and trisomies for all chromosomes; ploidy patterns were different among all strains/species. These results extend those we previously described in Leishmania major, demonstrating that mosaic aneuploidy is a genetic feature widespread to the Leishmania genus. In addition to the genetic consequences induced by this mosaicism, the apparent absence of alternation between haploid/diploid stages questions the modality of genetic exchange in Leishmania sp.

Novel insights into genome plasticity in Eukaryotes: mosaic aneuploidy in Leishmania

Leishmania are unicellular eukaryotes that have many markedly original molecular features compared with other uni- or multicellular eukaryotes like yeasts or mammals. Genome plasticity in this parasite has been the subject of many publications, and has been associated with drug resistance or adaptability. Aneuploidy has been suspected by several authors and it is now confirmed using state-of-the-art technologies such as high-throughput DNA sequencing. The analysis of genome contents at the single cell level using fluorescence in situ hybridization (FISH) has brought a new light on the genome organization: within a cell population, every chromosome, in every cell, may be present in at least two ploidy states (being either monosomic, disomic or trisomic), and the chromosomal content varies greatly from cell to cell, thus generating a constitutive intra-strain genomic heterogeneity, here termed 'mosaic aneuploidy'. Mosaic aneuploidy deeply affects the genetics of these organisms, leading, for example, to an extreme degree of intra-strain genomic diversity, as well as to a clearance of heterozygous cells in the population without however affecting genetic heterogeneity. Second, mosaic aneuploidy might be considered as a powerful strategy evolved by the parasite for adapting to modifications of environment conditions as well as for the emergence of drug resistance. On the whole, mosaic aneuploidy may be considered as a novel mechanism for generating phenotypic diversity driven by genomic plasticity.

Leishmania RNA virus: when the host pays the toll

The presence of an RNA virus in a South American subgenus of the Leishmania parasite, L. (Viannia), was detected several decades ago but its role in leishmanial virulence and metastasis was only recently described. In Leishmania guyanensis, the nucleic acid of Leishmania RNA virus (LRV1) acts as a potent innate immunogen, eliciting a hyper-inflammatory immune response through toll-like receptor 3 (TLR3). The resultant inflammatory cascade has been shown to increase disease severity, parasite persistence, and perhaps even resistance to anti-leishmanial drugs. Curiously, LRVs were found mostly in clinical isolates prone to infectious metastasis in both their human source and experimental animal model, suggesting an association between the viral hyperpathogen and metastatic complications such as mucocutaneous leishmaniasis (MCL). MCL presents as chronic secondary lesions in the mucosa of the mouth and nose, debilitatingly inflamed and notoriously refractory to treatment. Immunologically, this outcome has many of the same hallmarks associated with the reaction to LRV: production of type 1 interferons, bias toward a chronic Th1 inflammatory state and an impaired ability of host cells to eliminate parasites through oxidative stress. More intriguing, is that the risk of developing MCL is found almost exclusively in infections of the L. (Viannia) subtype, further indication that leishmanial metastasis is caused, at least in part, by a parasitic component. LRV present in this subgenus may contribute to the destructive inflammation of metastatic disease either by acting in concert with other intrinsic "metastatic factors" or by independently preying on host TLR3 hypersensitivity. Because LRV amplifies parasite virulence, its presence may provide a unique target for diagnostic and clinical intervention of metastatic leishmaniasis. Taking examples from other members of the Totiviridae virus family, this paper reviews the benefits and costs of endosymbiosis, specifically for the maintenance of LRV infection in Leishmania parasites, which is often at the expense of its human host.

FISH analysis reveals aneuploidy and continual generation of chromosomal mosaicism in Leishmania major

The protozoan parasite Leishmania is generally considered to be diploid, although a few chromosomes have been described as aneuploid. Using fluorescence in situ hybridization (FISH), we determined the number of homologous chromosomes per individual cell in L. major (i) during interphase and (ii) during mitosis. We show that, in Leishmania, aneuploidy appears to be the rule, as it affects all the chromosomes that we studied. Moreover, every chromosome was observed in at least two ploidy states, among monosomic, disomic or trisomic, in the cell population. This variable chromosomal ploidy among individual cells generates intra-strain heterogeneity, here precisely chromosomal mosaicism. We also show that this mosaicism, hence chromosome ploidy distribution, is variable among clones and strains. Finally, when we examined dividing nuclei, we found a surprisingly high rate of asymmetric chromosome allotments, showing that the transmission of genetic material during mitosis is highly unstable in this 'divergent' eukaryote: this leads to continual generation of chromosomal mosaicism. Using these results, we propose a model for the occurrence and persistence of this mosaicism. We discuss the implications of this additional unique feature of Leishmania for its biology and genetics, in particular as a novel genetic mechanism to generate phenotypic variability from genomic plasticity.

Regulation of gene expression in protozoa parasites

Infections with protozoa parasites are associated with high burdens of morbidity and mortality across the developing world. Despite extensive efforts to control the transmission of these parasites, the spread of populations resistant to drugs and the lack of effective vaccines against them contribute to their persistence as major public health problems. Parasites should perform a strict control on the expression of genes involved in their pathogenicity, differentiation, immune evasion, or drug resistance, and the comprehension of the mechanisms implicated in that control could help to develop novel therapeutic strategies. However, until now these mechanisms are poorly understood in protozoa. Recent investigations into gene expression in protozoa parasites suggest that they possess many of the canonical machineries employed by higher eukaryotes for the control of gene expression at transcriptional, posttranscriptional, and epigenetic levels, but they also contain exclusive mechanisms. Here, we review the current understanding about the regulation of gene expression in Plasmodium sp., Trypanosomatids, Entamoeba histolytica and Trichomonas vaginalis.

Demonstration of genetic exchange during cyclical development of Leishmania in the sand fly vector

Genetic exchange has not been shown to be a mechanism underlying the extensive diversity of Leishmania parasites. We report here evidence that the invertebrate stages of Leishmania are capable of having a sexual cycle consistent with a meiotic process like that described for African trypanosomes. Hybrid progeny were generated that bore full genomic complements from both parents, but kinetoplast DNA maxicircles from one parent. Mating occurred only in the sand fly vector, and hybrids were transmitted to the mammalian host by sand fly bite. Genetic exchange likely contributes to phenotypic diversity in natural populations, and analysis of hybrid progeny will be useful for positional cloning of the genes controlling traits such as virulence, tissue tropism, and drug resistance.

Leishmania (Viannia) braziliensis transfectants overexpressing the miniexon gene lose virulence in vivo

The miniexon gene has a central role in the processing of polycistronic pre-mRNA of kinetoplastids. It is added to the 5' extremity of each mRNA, supplying the 5'-capped structure to the molecule. Previous studies in Leishmania (Leishmania) major showed that the overexpression of the miniexon array attenuates the virulence of the parasite in in vivo assays. The results presented here extend those findings to Viannia subgenus. Leishmania (Viannia) braziliensis was transfected with a cosmid harboring a tandem array of one hundred miniexon gene copies and then characterized by Northern blot analysis. The overexpression of the exogenous gene was confirmed and its effect on the virulence of L. (V.) braziliensis was investigated in hamsters. In BALB/c mice we could not detect parasites during the course of 15 weeks of infection. In addition, hamsters infected with transfectants overexpressing the miniexon gene exhibited only a minor footpad swelling of late onset and failed to develop progressive lesion, these attenuated parasites could be recovered from the inoculation site 1 year after infection. The persistence of parasites in the host indicates that a stable line overexpressing the miniexon may be tested as live vaccine against leishmaniasis.

First complete chromosomal organization of a protozoan plant parasite (Phytomonas spp.)

Phytomonas spp. are members of the family Trypanosomatidae that parasitize plants and may cause lethal diseases in crops such as Coffee Phloem necrosis, Hartrot in coconut, and Marchitez sorpresiva in oil palm. In this study, the molecular karyotype of 6 isolates from latex plants has been entirely elucidated by pulsed-field gel electrophoresis and DNA hybridization. Twenty-one chromosomal linkage groups constituting heterologous chromosomes and sizing between 0.3 and 3 Mb could be physically defined by the use of 75 DNA markers (sequence-tagged sites and genes). From these data, the genome size can be estimated at 25.5 (+/-2) Mb. The physical linkage groups were consistently conserved in all strains examined. Moreover, the finding of several pairs of different-sized homologous chromosomes strongly suggest diploidy for this organism. The definition of the complete molecular karyotype of Phytomonas represents an essential primary step toward sequencing the genome of this parasite of economical importance.

Cytogenetic evidence for diversity of two nuclei within a single diplomonad cell of Giardia

Giardia intestinalis is an ancient protist that causes the most commonly reported human diarrheal disease of parasitic origin worldwide. An intriguing feature of the Giardia cell is the presence of two morphologically similar nuclei, generally considered equivalent, in spite of the fact that their karyotypes are unknown. We found that within a single cell, the two nuclei differ both in the number and the size of chromosomes and that representatives of two major genetic groups of G. intestinalis possess different karyotypes. Odd chromosome numbers indicate aneuploidy of Giardia nuclei, and their stable occurrence is suggestive of a long-term asexuality. A semi-open type of Giardia mitosis excludes a chromosome interfusion between the nuclei. Differences in karyotype and DNA content, and cell cycle-dependent asynchrony are indicative of diversity of the two Giardia nuclei.

Cell cycle-dependent expression regulation by the proteasome pathway and characterization of the nuclear targeting signal of a Leishmania major Kin-13 kinesin

The LmjF01.0030 gene of Leishmania major Friedlin, annotated as 'MCAK-like', was confirmed as a kinesin with an internally located motor domain and termed LmjKIN13-1. Both the native form of the protein and a green fluorescent protein (GFP)-fused recombinant version were shown to be exclusively intranuclear, and, more specifically, to localize to the spindle and spindle poles. Cell cycle-dependent regulation of the protein levels was demonstrated using synchronized Leishmania cells: LmjKIN13-1 was highly abundant in the G2+M phase and present at very low levels after mitosis. Altogether, these features suggest that this protein participates in mitosis. The construction of systematic deletion mutants allowed the localization of the primary sequence regions responsible for nuclear targeting on the one hand, and for cell cycle-dependent variations on the other hand. A 42-amino-acid region of the carboxy(C)-terminal domain mediates nuclear import and could be defined as an atypical nuclear localization signal. Protein level regulation during the cell cycle was shown to also depend upon the C-terminal domain, where apparently redundant degradation signals are present. Putative degradation signals appear to be present on both sides and inside the nuclear localization signal. Further experiments strongly suggest a role for the ubiquitin/proteasome pathway in this cell cycle-dependent regulation. These data underline the importance of post-translational regulation of protein abundance in this ancestral eukaryote where transcriptional regulation seems to be rare or near absent.

Ploidy changes associated with disruption of two adjacent genes on Leishmania major chromosome 1

Leishmania major Friedlin (LmjF) is a kinetoplastid protozoan whose genomic sequence has been recently elucidated. About 60% of the identified genes do not have a known function, and many are trypanosomatid-specific. Here we characterise two adjacent genes from LmjF chromosome 1 (chr1): LmjF01.0750, which encodes a predicted protein with a serine/threonine protein kinase motif and LmjF01.0760, which encodes a product with no similarity to other known proteins. Orthologues of both genes are present in Trypanosoma cruzi, but neither occur in Trypanosoma brucei. We have mapped polyadenylation and spliced-leader acceptor sites for both genes, and show that they differ between Leishmania species. Attempts to generate null mutants of LmjF01.0750 by homologous recombination were unsuccessful and led to the apparent triploidy of the entire genome, suggesting that it is an essential gene. Interestingly, at least two copies of LmjF01.0750 are required for cell survival. Further evidence of genome plasticity in Leishmania was provided by changes in chr1 copy number that occurred during in vitro growth of wild-type LmjF promastigotes and following replacement of a single copy of LmjF01.0760.

Shotgun optical mapping of the entire Leishmania major Friedlin genome

Leishmania is a group of protozoan parasites which causes a broad spectrum of diseases resulting in widespread human suffering and death, as well as economic loss from the infection of some domestic animals and wildlife. To further understand the fundamental genomic architecture of this parasite, and to accelerate the on-going sequencing project, a whole-genome XbaI restriction map was constructed using the optical mapping system. This map supplemented traditional physical maps that were generated by fingerprinting and hybridization of cosmid and P1 clone libraries. Thirty-six optical map contigs were constructed for the corresponding known 36 chromosomes of the Leishmania major Friedlin genome. The chromosome sizes ranged from 326.9 to 2821.3 kb, with a total genome size of 34.7 Mb; the average XbaI restriction fragment was 25.3 kb, and ranged from 15.7 to 77.8 kb on a per chromosomes basis. Comparison between the optical maps and the in silico maps of sequence drawn from completed, nearly finished, or large sequence contigs showed that optical maps served several useful functions within the path to create finished sequence by: guiding aspects of the sequence assembly, identifying misassemblies, detection of cosmid or PAC clones misplacements to chromosomes, and validation of sequence stemming from varying degrees of finishing. Our results also showed the potential use of optical maps as a means to detect and characterize map segmental duplication within genomes.

Genomic organization and gene expression in a chromosomal region of Leishmania major

Little is known about the relation between the genome organization and gene expression in Leishmania. Bioinformatic analysis can be used to predict genes and find homologies with known proteins. A model was proposed, in which genes are organized into large clusters and transcribed from only one strand, in the form of large polycistronic primary transcripts. To verify the validity of this model, we studied gene expression at the transcriptional, post-transcriptional and translational levels in a unique locus of 34kb located on chr27 and represented by cosmid L979. Sequence analysis revealed 115 ORFs on either DNA strand. Using computer programs developed for Leishmania genes, only nine of these ORFs, localized on the same strand, were predicted to code for proteins, some of which show homologies with known proteins. Additionally, one pseudogene, was identified. We verified the biological relevance of these predictions. mRNAs from nine predicted genes and proteins from seven were detected. Nuclear run-on analyses confirmed that the top strand is transcribed by RNA polymerase II and suggested that there is no polymerase entry site. Low levels of transcription were detected in regions of the bottom strand and stable transcripts were identified for four ORFs on this strand not predicted to be protein-coding. In conclusion, the transcriptional organization of the Leishmania genome is complex, raising the possibility that computer predictions may not be comprehensive.

Leishmania major chromosome 3 contains two long convergent polycistronic gene clusters separated by a tRNA gene

Leishmania parasites (order Kinetoplastida, family Trypanosomatidae) cause a spectrum of human diseases ranging from asymptomatic to lethal. The approximately 33.6 Mb genome is distributed among 36 chromosome pairs that range in size from approximately 0.3 to 2.8 Mb. The complete nucleotide sequence of Leishmania major Friedlin chromosome 1 revealed 79 protein-coding genes organized into two divergent polycistronic gene clusters with the mRNAs transcribed towards the telomeres. We report here the complete nucleotide sequence of chromosome 3 (384 518 bp) and an analysis revealing 95 putative protein-coding ORFs. The ORFs are primarily organized into two large convergent polycistronic gene clusters (i.e. transcribed from the telomeres). In addition, a single gene at the left end is transcribed divergently towards the telomere, and a tRNA gene separates the two convergent gene clusters. Numerous genes have been identified, including those for metabolic enzymes, kinases, transporters, ribosomal proteins, spliceosome components, helicases, an RNA-binding protein and a DNA primase subunit.

The switch region on Leishmania major chromosome 1 is not required for mitotic stability or gene expression, but appears to be essential

The Leishmania genome project reference strain, Leishmania major Friedlin, is trisomic for chromosome 1. The complete sequence of this chromosome has revealed that the genes are grouped into two large clusters of the polycistronic type, each borne by one DNA strand and located on each side of a 1.6-kb sequence often termed the switch region. Several hypotheses concerning the role of this switch region have been put forward (region of initiation of transcription for both gene clusters, origin of replication or centromeric sequence). In the present study, we have deleted this region on the three copies of chromosome 1 by sequential targeted replacements. The absence of the switch region did not alter the mitotic stability of the three deleted chromosomes. This region therefore does not appear necessary for chromosomal replication or segregation. However, during the third targeting round, which aimed at knocking out the last switch region, a fourth copy of chromosome 1 that retained this region appeared in all clones analysed. This suggests that the persistence of this switch region is necessary for parasite survival. We then showed that the presence/absence of the switch region did not act upon the expression of a resistance marker gene inserted beforehand into the left gene cluster of the same chromosomal molecule. This result suggests that the presence of this 1.6-kb sequence is not necessary for the expression of all genes on chromosome 1.

Mitotic stability of a coding DNA sequence-free version of Leishmania major chromosome 1 generated by targeted chromosome fragmentation

The deletion of a 260-kb segment containing all the coding DNA sequences (CDS) of chromosome 1 of Leishmania major Friedlin strain was performed through homologous recombination during a transfection experiment. This allowed the selection of a mutant clone containing a linear extra chromosome sizing 155 kb (XC155). The structure of XC155 was determined by restriction analysis and DNA cloning and sequencing of the gel-purified chromosome: it is made of a 'mirror' inverted duplication of the 'right' end of chromosome 1a (approximately 25 kb at each end), and in its central part of a complex tandem amplification of the linearized transfection vector containing the hygromycin resistance gene (over approximately 105 kb). No sequence of the coding region of chromosome 1 (including the 1.6-kb 'switch' region) was found. By contrast, XC155 contains two large (approximately 13 kb) clusters of tandemly repeated subtelomeric sequences (272-bp 'satellite' DNA) as well as telomeric hexamer repeats. This extra chromosome was found to be mitotically stable after >150 generations without selective pressure in vitro. Two sequence elements are considered which may have an effect on mitotic stability and participate to centromeric function in this extra chromosome: the amplification of the input vector and the 272-bp 'satellite' DNA bound by telomeric repeats.

Genetic and structural heterogeneity of proteophosphoglycans in Leishmania

Proteophosphoglycans (PPG) are a large family of extensively glycosylated proteins with some unusual and unique features. The ppg gene family is conserved in at least three Leishmania species and localises to chromosome 35. Previous studies using standard discontinuous SDS-PAGE have been incapable of resolving PPG heterogeneity with most material failing to enter the resolving gel. We have exploited a continuous electrophoretic system, which allows for the first time the separation and characterisation of a low molecular weight population of PPG polypeptides. We provide evidence of surface expressed and developmentally regulated forms. Among those, we identify for the first time the previously described membrane-bound PPG and a form of filamentous fPPG, which is altered, or absent in two of the three L. major isolates examined.

Charaterization of Leishmania major Friedlin telomeric terminus

Here we have characterized Leishmania major (Friedlin) telomeric terminus (the very end) using recombinants obtained by a vector-adaptor cloning protocol. As in L. donovani, the last nine nucleotides of L. major terminus are 5'-GGTTAGGGT-OH 3', differing from Trypanosoma cruzi and T. brucei terminus 5'GGGTTAGGG-OH 3', thus indicating that these sequences are genus specific. We have also made a comparative analysis between L. major and L. donovani telomere-associated sequences, and described a novel non-repeated telomeric associated sequence common to L. major low molecular weight chromosomal bands.

Genomic organization and functional characterization of the Leishmania major Friedlin ribosomal RNA gene locus

The sequence and gene organization of the ribosomal RNA (rRNA) genes of Leishmania major Friedlin (LmjF) were determined. Interestingly, the rDNA repeat unit contained a duplicated 526 bp fragment at the 3' end of the unit with two copies of the LSUepsilon rRNA gene. Our results suggested the presence of only approximately 24 copies of the rRNA unit per diploid genome in LmjF. Repetitive elements (IGSRE) of 63 bp occurred in the intergenic spacer (IGS) between the LSUepsilon and the SSU rRNA genes. Among the different rDNA units, the region containing the IGSRE fluctuated in length from approximately 1.3 to approximately 18 kb. The transcription initiation site (TIS) of the rRNA unit was localized by primer extension to 1043 bp upstream of the SSU gene and 184 bp downstream of the IGSRE. Sequence comparison among several species of Leishmania showed a high degree of conservation around the TIS. Moreover, the IGSRE also showed considerable similarity between Leishmania species. In transient transfection assays, a fragment containing the TIS directed a 164- to 178-fold increase in luciferase activity over the no-insert control, indicating the presence of a promoter within this 391 bp fragment. The LmjF promoter region was also functional in other species of Leishmania. Nuclear run-on analyses demonstrated that only the rRNA-coding strand is transcribed, downstream of this RNA polymerase I (pol I) promoter. These experiments also suggested that transcription terminates upstream of the IGSRE.

Effect of large targeted deletions on the mitotic stability of an extra chromosome mediating drug resistance in Leishmania

A mitotically stable linear extra chromosome obtained in a Leishmania donovani strain rendered mycophenolic acid-resistant has been physically mapped. This 290-kb chromosome has an inverted duplicated structure around a central inversion region, and is derived from a conservative amplification event of a approximately 140-kb subtelomeric end of chromosome 19. Large-sized targeted deletions of the central region were performed through homologous recombination using three specific transfection vectors. The size of the extra chromosome was thus successfully reduced from 290 to 260, 200 and 120 kb respectively. The mitotic stability of these chromosomes was then analysed in drug-free cultures over >140 days. Results differed according to the deletion created. By contrast with the smallest deletion the two largest deletions altered mitotic stability, leading to progressive loss of the size-reduced chromosomes with similar kinetics in both mutants. The 30-kb region common to both deletions may therefore be considered as involved in mitotic stability. A 44-kb contig covering this region could be assembled and sequenced. The analysis of this sequence did not reveal any sequence elements typical of centromeric DNA. By contrast, its enrichment in homopolymer tracts suggests that this region might contain an origin of replication.

Characterisation of three chromosomal ends of Leishmania major reveals transcriptional activity across arrays of reiterated and unique sequences

The 36 chromosomes of the parasite Leishmania major range in size from 200 kb to approximately 2.5 Mb and variation between homologues seems to be restricted to the telomeric and subtelomeric regions. We have isolated three cosmids carrying the telomere hexameric repeat and assigned them to the extreme location of chromosomes 3, 7 and 20. When considering the distribution of repetitive sequences, Southern analysis of the three chromosomal ends indicated the existence of at least two classes of chromosomal extremities: one of them is composed almost exclusively of unique sequences and the other is characterised by patches of both reiterated and unique sequences. We devised a transfection-based strategy that allowed the determination of a map of transcripts in each of the regions examined. Sequencing of the chromosome 20 cosmid revealed the existence of a novel class of reiterated sequence, LST-R378, and 10 ORFs drawing a map of putative genes compatible with the map of transcripts.

A random sequencing approach for the analysis of the Trypanosoma cruzi genome: general structure, large gene and repetitive DNA families, and gene discovery

A random sequence survey of the genome of Trypanosoma cruzi, the agent of Chagas disease, was performed and 11,459 genomic sequences were obtained, resulting in approximately 4.3 Mb of readable sequences or approximately 10% of the parasite haploid genome. The estimated total GC content was 50.9%, with a high representation of A and T di- and trinucleotide repeats. Out of the estimated 5000 parasite genes, 947 putative new genes were identified. Another 1723 sequences corresponded to genes detected previously in T. cruzi through expression sequence tag analysis. 7735 sequences had no matches in the database, but the presence of open reading frames that passed Fickett's test suggests that some might contain coding DNA. The survey was highly redundant, with approximately 35% of the sequences included in a few large sequence families. Some of them code for protein families present in dozens of copies, including proteins essential for parasite survival and retrotransposons. Other sequence families include repetitive DNA present in thousands of copies per haploid genome. Some families in the latter group are new, parasite-specific, repetitive DNAs. These results suggest that T. cruzi could constitute an interesting model to analyze gene and genome evolution due to its plasticity in terms of sequence amplification and divergence. Additional information can be found at http://www.iib.unsam.edu.ar/tcruzi.gss. html.

Chromosome structure and sequence organization between pathogenic and non-pathogenic Leishmania spp

We have used a chromosome fragmentation strategy based on systematic genomic insertions of the rare cutting yeast I-SceI endonuclease to assess structure and sequence organization of homologous chromosomes between evolutionary divergent pathogenic and non-pathogenic Leishmania species. This method was combined to physical mapping and hybridization studies using a number of specific chromosomal markers as probes. Our studies have concentrated on two different chromosomes of Leishmania major (L. major), L. donovani and L. infantum and of the non-pathogenic species L. tarentolae. Specific chromosome fragmentation events at the level of multiple I-SccI genomic integrations indicated that very similar distances separated internal genomic sequences between homologous chromosomes and that distances from chromosome ends were more variable. The order and orientation of genes along the homologous chromosomes were also conserved between species. With only few exceptions, genome organization between pathogenic and non-pathogenic Leishmania species was found to be highly conserved. Genomic comparison of pathogenic and non-pathogenic species may be useful for depicting regions involved in species-specific related pathologies.

A Random Sequencing Approach for the Analysis of the Trypanosoma cruzi Genome: General Structure, Large Gene and Repetitive DNA Families, and Gene Discovery

A random sequence survey of the genome of Trypanosoma cruzi, the agent of Chagas disease, was performed and 11,459 genomic sequences were obtained, resulting in ∼4.3 Mb of readable sequences or ∼10% of the parasite haploid genome. The estimated total GC content was 50.9%, with a high representation of A and T di- and trinucleotide repeats. Out of the estimated 5000 parasite genes, 947 putative new genes were identified. Another 1723 sequences corresponded to genes detected previously in T. cruzi through expression sequence tag analysis. 7735 sequences had no matches in the database, but the presence of open reading frames that passed Fickett's test suggests that some might contain coding DNA. The survey was highly redundant, with ∼35% of the sequences included in a few large sequence families. Some of them code for protein families present in dozens of copies, including proteins essential for parasite survival and retrotransposons. Other sequence families include repetitive DNA present in thousands of copies per haploid genome. Some families in the latter group are new, parasite-specific, repetitive DNAs. These results suggest that T. cruzi could constitute an interesting model to analyze gene and genome evolution due to its plasticity in terms of sequence amplification and divergence. Additional information can be found at http://www.iib.unsam.edu.ar/tcruzi.gss.html.

[The sequence data described in this paper have been submitted to the dbGSS database under the following GenBank accession nos.:AQ443439–AQ443513, AQ443743–AQ445667, AQ902981–AQ911366,AZ049857–AZ051184, and AZ302116–AZ302563.]

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.

The size difference between leishmania major friedlin chromosome one homologues is localized to sub-telomeric repeats at one chromosomal end

Leishmania species are members of the evolutionarily ancient protozoan order Kinetoplastidae and are important human pathogens. The Leishmania genome is relatively small (approximately 34 Mbp) and is distributed among 36 chromosome pairs, ranging in size from 0.3 to 2.5 Mbp. The smallest chromosome of Leishmania major Friedlin, chrl, consists of three homologues which differ in size by approximately 29 kb. Previous sequence and Southern analyses of all three homologues reveal a conserved chromosomal core, consisting of coding and adjacent 'non-informational' sequence. Here we show the size difference between homologues is largely restricted to variation in both the number and content of several sub-telomeric repetitive elements localized on one chromosomal end. These repetitive elements also occur on other chromosomes, but some are more dispersed in the Leishmania genome than others.

The Leishmania genome comes of Age

The Leishmania Genome Network (LGN) was born in Rio de Janeiro, Brazil in 1994. In the short period that has elapsed since then, the LGN has focused solely on the acquisition of the resources, and hence data, that have enabled a rational approach to genomic sequencing of the reference strain, Leishmania major Friedlin. This has now been achieved. In this review, Alasdair Ivens and Jennie Blackwell, secretary and chairman of the LGN, respectively, re-examine the approaches that were adopted, comment on some of the interesting data that have been obtained and introduce some genome-wide approaches that will facilitate functional studies of the parasite.

Altered ligands reveal limited plasticity in the T cell response to a pathogenic epitope

Experimental leishmaniasis offers a well characterized model of T helper type 1 cell (Th1)-mediated control of infection by an intracellular organism. Susceptible BALB/c mice aberrantly develop Th2 cells in response to infection and are unable to control parasite dissemination. The early CD4(+) T cell response in these mice is oligoclonal and reflects the expansion of Vbeta4/ Valpha8-bearing T cells in response to a single epitope from the parasite Leishmania homologue of mammalian RACK1 (LACK) antigen. Interleukin 4 (IL-4) generated by these cells is believed to direct the subsequent Th2 response. We used T cells from T cell receptor-transgenic mice expressing such a Vbeta4/Valpha8 receptor to characterize altered peptide ligands with similar affinity for I-Ad. Such altered ligands failed to activate IL-4 production from transgenic LACK-specific T cells or following injection into BALB/c mice. Pretreatment of susceptible mice with altered peptide ligands substantially altered the course of subsequent infection. The ability to confer a healer phenotype on otherwise susceptible mice using altered peptides that differed by a single amino acid suggests limited diversity in the endogenous T cell repertoire recognizing this antigen.

Conserved linkage groups associated with large-scale chromosomal rearrangements between Old World and New World Leishmania genomes

The genus Leishmania can be taxonomically separated into three main groups: the Old World subgenus L. (Leishmania), the New World subgenus L. (Leishmania) and the New World subgenus L. (Viannia). The haploid genome of Old World Leishmania species has been shown to contain 36 chromosomes defined as physical linkage groups; the latter were found entirely conserved across species. In the present study, we tried to verify whether this conservation of the genome structure extends to the New World species of Leishmania. 300 loci were explored by hybridization on optimized pulsed field gel electrophoresis separations of the chromosomes of polymorphic strains of the six main pathogenic Leishmania species of the New World. When comparing these New World karyotypes with their Old World counterparts, 32 out of 36 linkage groups were found conserved among all species. Four chromosomal rearrangements were found. All species belonging to the L. (Viannia) subgenus were characterized by the presence (i) of a short sequence exchange between chromosomes 26 and 35, and (ii) more importantly, of a fused version of chromosomes 20 and 34 which are separated in all Old World species. 69 additional markers were isolated from a plasmid library specifically constructed from the rearranged chromosomes 20+34 in an attempt to detect mechanisms other than a fusion or breakage: only two markers out of 40 did not belong to the linkage groups 20 and 34. On the other hand, all strains belonging to the New World subgenus L. (Leishmania) were characterized by two different chromosomal rearrangements of the same type (fusion/breakage) as above as compared with Old World species: chromosomes 8+29 and 20+36. Consequently, these two groups of species have 35 and 34 heterologous chromosomes, respectively. Overall, these results show that large-scale chromosomal rearrangements occurred during the evolution of the genus Leishmania, and that the three main groups of pathogenic species are characterized by different chromosome numbers. Nevertheless, translocations seem particularly rare, and the conservation of the major linkage groups should be an essential feature for the compared genetics between species of this parasite.