Genes and chromosomes of Leishmania infantum

Gene Annotation and Transcriptome Delineation on a De Novo Genome Assembly for the Reference Leishmania major Friedlin Strain

Leishmania major is the main causative agent of cutaneous leishmaniasis in humans. The Friedlin strain of this species (LmjF) was chosen when a multi-laboratory consortium undertook the objective of deciphering the first genome sequence for a parasite of the genus Leishmania. The objective was successfully attained in 2005, and this represented a milestone for Leishmania molecular biology studies around the world. Although the LmjF genome sequence was done following a shotgun strategy and using classical Sanger sequencing, the results were excellent, and this genome assembly served as the reference for subsequent genome assemblies in other Leishmania species. Here, we present a new assembly for the genome of this strain (named LMJFC for clarity), generated by the combination of two high throughput sequencing platforms, Illumina short-read sequencing and PacBio Single Molecular Real-Time (SMRT) sequencing, which provides long-read sequences. Apart from resolving uncertain nucleotide positions, several genomic regions were reorganized and a more precise composition of tandemly repeated gene loci was attained. Additionally, the genome annotation was improved by adding 542 genes and more accurate coding-sequences defined for around two hundred genes, based on the transcriptome delimitation also carried out in this work. As a result, we are providing gene models (including untranslated regions and introns) for 11,238 genes. Genomic information ultimately determines the biology of every organism; therefore, our understanding of molecular mechanisms will depend on the availability of precise genome sequences and accurate gene annotations. In this regard, this work is providing an improved genome sequence and updated transcriptome annotations for the reference L. major Friedlin strain.

Dipeptidyl peptidase III as a DNA marker to investigate epidemiology and taxonomy of Old World Leishmania species

BACKGROUND

Dipeptidyl peptidase III (DPPIII) member of M49 peptidase family is a zinc-dependent metallopeptidase that cleaves dipeptides sequentially from the N-terminus of its substrates. In Leishmania, DPPIII, was reported with other peptidases to play a significant role in parasites' growth and survival. In a previous study, we used a coding sequence annotated as DPPIII to develop and evaluate a PCR assay that is specific to dermotropic Old World (OW) Leishmania species. Thus, our objective was to further assess use of this gene for Leishmania species identification and for phylogeny, and thus for diagnostic and molecular epidemiology studies of Old World Leishmania species.

METHODOLOGY

Orthologous DDPIII genes were searched in all Leishmania genomes and aligned to design PCR primers and identify relevant restriction enzymes. A PCR assays was developed and seventy-two Leishmania fragment sequences were analyzed using MEGA X genetics software to infer evolution and phylogenetic relationships of studied species and strains. A PCR-RFLP scheme was also designed and tested on 58 OW Leishmania strains belonging to 8 Leishmania species and evaluated on 75 human clinical skin samples.

FINDINGS

Sequence analysis showed 478 variable sites (302 being parsimony informative). Test of natural selection (dN-dS) (-0.164, SE = 0.013) inferred a negative selection, characteristic of essential genes, corroborating the DPPIII importance for parasite survival. Inter- and intra-specific genetic diversity was used to develop universal amplification of a 662bp fragment. Sequence analyses and phylogenies confirmed occurrence of 6 clusters congruent to L. major, L. tropica, L. aethiopica, L. arabica, L. turanica, L. tarentolae species, and one to the L. infantum and L. donovani species complex. A PCR-RFLP algorithm for Leishmania species identification was designed using double digestions with HaeIII and KpnI and with SacI and PvuII endonucleases. Overall, this PCR-RFLP yielded distinct profiles for each of the species L. major, L. tropica, L. aethiopica, L. arabica and L. turanica and the L. (Sauroleishmania) L. tarentolae. The species L. donovani, and L. infantum shared the same profile except for strains of Indian origin. When tested on clinical samples, the DPPIII PCR showed sensitivities of 82.22% when compared to direct examination and was able to identify 84.78% of the positive samples.

CONCLUSION

The study demonstrates that DPPIII gene is suitable to detect and identify Leishmania species and to complement other molecular methods for leishmaniases diagnosis and epidemiology. Thus, it can contribute to evidence-based disease control and surveillance.

Complete assembly of the Leishmania donovani (HU3 strain) genome and transcriptome annotation

Leishmania donovani is a unicellular parasite that causes visceral leishmaniasis, a fatal disease in humans. In this study, a complete assembly of the genome of L. donovani is provided. Apart from being the first published genome of this strain (HU3), this constitutes the best assembly for an L. donovani genome attained to date. The use of a combination of sequencing platforms enabled to assemble, without any sequence gap, the 36 chromosomes for this species. Additionally, based on this assembly and using RNA-seq reads derived from poly-A + RNA, the transcriptome for this species, not yet available, was delineated. Alternative SL addition sites and heterogeneity in the poly-A addition sites were commonly observed for most of the genes. After a complete annotation of the transcriptome, 2,410 novel transcripts were defined. Additionally, the relative expression for all transcripts present in the promastigote stage was determined. Events of cis-splicing have been documented to occur during the maturation of the transcripts derived from genes LDHU3_07.0430 and LDHU3_29.3990. The complete genome assembly and the availability of the gene models (including annotation of untranslated regions) are important pieces to understand how differential gene expression occurs in this pathogen, and to decipher phenotypic peculiarities like tissue tropism, clinical disease, and drug susceptibility.

Leishmania infections: Molecular targets and diagnosis

Progress in the diagnosis of leishmaniases depends on the development of effective methods and the discovery of suitable biomarkers. We propose firstly an update classification of Leishmania species and their synonymies. We demonstrate a global map highlighting the geography of known endemic Leishmania species pathogenic to humans. We summarize a complete list of techniques currently in use and discuss their advantages and limitations. The available data highlights the benefits of molecular markers in terms of their sensitivity and specificity to quantify variation from the subgeneric level to species complexes, (sub) species within complexes, and individual populations and infection foci. Each DNA-based detection method is supplied with a comprehensive description of markers and primers and proposal for a classification based on the role of each target and primer in the detection, identification and quantification of leishmaniasis infection. We outline a genome-wide map of genes informative for diagnosis that have been used for Leishmania genotyping. Furthermore, we propose a classification method based on the suitability of well-studied molecular markers for typing the 21 known Leishmania species pathogenic to humans. This can be applied to newly discovered species and to hybrid strains originating from inter-species crosses. Developing more effective and sensitive diagnostic methods and biomarkers is vital for enhancing Leishmania infection control programs.

Leishmania-specific surface antigens show sub-genus sequence variation and immune recognition

Background

A family of hydrophilic acylated surface (HASP) proteins, containing extensive and variant amino acid repeats, is expressed at the plasma membrane in infective extracellular (metacyclic) and intracellular (amastigote) stages of Old World Leishmania species. While HASPs are antigenic in the host and can induce protective immune responses, the biological functions of these Leishmania-specific proteins remain unresolved. Previous genome analysis has suggested that parasites of the sub-genus Leishmania (Viannia) have lost HASP genes from their genomes.

Methods/Principal Findings

We have used molecular and cellular methods to analyse HASP expression in New World Leishmania mexicana complex species and show that, unlike in L. major, these proteins are expressed predominantly following differentiation into amastigotes within macrophages. Further genome analysis has revealed that the L. (Viannia) species, L. (V.) braziliensis, does express HASP-like proteins of low amino acid similarity but with similar biochemical characteristics, from genes present on a region of chromosome 23 that is syntenic with the HASP/SHERP locus in Old World Leishmania species and the L. (L.) mexicana complex. A related gene is also present in Leptomonas seymouri and this may represent the ancestral copy of these Leishmania-genus specific sequences. The L. braziliensis HASP-like proteins (named the orthologous (o) HASPs) are predominantly expressed on the plasma membrane in amastigotes and are recognised by immune sera taken from 4 out of 6 leishmaniasis patients tested in an endemic region of Brazil. Analysis of the repetitive domains of the oHASPs has shown considerable genetic variation in parasite isolates taken from the same patients, suggesting that antigenic change may play a role in immune recognition of this protein family.

Conclusions/Significance

These findings confirm that antigenic hydrophilic acylated proteins are expressed from genes in the same chromosomal region in species across the genus Leishmania. These proteins are surface-exposed on amastigotes (although L. (L.) major parasites also express HASPB on the metacyclic plasma membrane). The central repetitive domains of the HASPs are highly variant in their amino acid sequences, both within and between species, consistent with a role in immune recognition in the host.

Single-celled Leishmania parasites, transmitted by sand flies, infect humans and other mammals in many tropical and sub-tropical regions, giving rise to a spectrum of diseases called the leishmaniases. Species of parasite within the Leishmania genus can be divided into two groups (referred to as sub-genera) that are separated by up to 100 million years of evolution yet are highly related at the genome level. Our research is focused on identifying gene differences between these sub-genera that may identify proteins that impact on the transmission and pathogenicity of different Leishmania species. Here we report the presence of a highly-variant genomic locus (OHL) that was previously described as absent in parasites of the L. (Viannia) subgenus (on the basis of lack of key genes) but is present and well-characterised (as the LmcDNA16 locus) in all members of the alternative subgenus, L. (Leishmania). We demonstrate that the proteins encoded within the LmcDNA16 and OHL loci are similar in their structure and surface localisation in mammalian-infective amastigotes, despite significant differences in their DNA sequences. Most importantly, we demonstrate that the OHL locus proteins, like the HASP proteins from the LmcDNA16 locus, contain highly variable amino acid repeats that are antigenic in man and may therefore contribute to future vaccine development.

Analysis of gene order conservation in eukaryotes identifies transcriptionally and functionally linked genes

The order of genes in eukaryotes is not entirely random. Studies of gene order conservation are important to understand genome evolution and to reveal mechanisms why certain neighboring genes are more difficult to separate during evolution. Here, genome-wide gene order information was compiled for 64 species, representing a wide variety of eukaryotic phyla. This information is presented in a browser where gene order may be displayed and compared between species. Factors related to non-random gene order in eukaryotes were examined by considering pairs of neighboring genes. The evolutionary conservation of gene pairs was studied with respect to relative transcriptional direction, intergenic distance and functional relationship as inferred by gene ontology. The results show that among gene pairs that are conserved the divergently and co-directionally transcribed genes are much more common than those that are convergently transcribed. Furthermore, highly conserved pairs, in particular those of fungi, are characterized by a short intergenic distance. Finally, gene pairs of metazoa and fungi that are evolutionary conserved and that are divergently transcribed are much more likely to be related by function as compared to poorly conserved gene pairs. One example is the ribosomal protein gene pair L13/S16, which is unusual as it occurs both in fungi and alveolates. A specific functional relationship between these two proteins is also suggested by the fact that they are part of the same operon in both eubacteria and archaea. In conclusion, factors associated with non-random gene order in eukaryotes include relative gene orientation, intergenic distance and functional relationships. It seems likely that certain pairs of genes are conserved because the genes involved have a transcriptional and/or functional relationship. The results also indicate that studies of gene order conservation aid in identifying genes that are related in terms of transcriptional control.

The translational efficiencies of the two Leishmania infantum HSP70 mRNAs, differing in their 3'-untranslated regions, are affected by shifts in the temperature of growth through different mechanisms

Exposure of Leishmania promastigotes to the temperature of their mammalian hosts induces a typical heat-shock response. In Leishmania infantum, HSP70 is encoded by two types of genes that differ in their 3'-untranslated regions (3'-UTRs). Previously, we have shown that specific transcripts for each gene are present in promastigotes growing at normal temperature (26 degrees C), but only transcripts with 3'-UTR-type I (3'-UTRI) accumulate in a temperature-dependent manner. Here, we have investigated the translational efficiencies of both types of HSP70 transcripts at the different temperatures that the parasite encounters in the insect (26 degrees C, normal temperature) or in the mammalian host (heat-shock temperatures). Interestingly, 3'-UTRI-bearing transcripts (HSP70-I) were found associated with ribosomes in promastigotes at normal and heat-shock temperatures, whereas the HSP70-II transcripts appear to be preferentially translated at heat-shock temperatures but not at 26 degrees C. We have analyzed the function of these UTRs in the translational control by use of plasmid constructs in which the CAT reporter gene was flanked by UTRs of the HSP70 genes. Unexpectedly, it was found that CAT transcripts with 3'-UTRII bind to ribosomes at 26 degrees C, and, indeed, the CAT protein is synthesized. A valid conclusion of these experiments was that both types of 3'-UTRs are essential for translation of HSP70 mRNAs at heat shock temperatures, although the 3'-UTRII is more efficient during severe heat shock (39 degrees C). In addition, these results suggest that sequence region other than the 3'-UTR of HSP70-II gene is involved in the translational silent state of HSP70-II transcripts at 26 degrees C. Finally, a null mutant has been created by targeted disruption of both HSP70-II alleles. Remarkably, the deltaHSP70 mutant synthesizes HSP70 at a lower rate than the wild-type parasites. Overall, our data suggest that the biological function of the HSP70-II gene is to top up HSP70 levels under conditions of stress.