Molecular karyotype analysis in Leishmania

Intrachromosomal amplification, locus deletion and point mutation in the aquaglyceroporin AQP1 gene in antimony resistant Leishmania (Viannia) guyanensis

BACKGROUND

Antimony resistance complicates the treatment of infections caused by the parasite Leishmania.

METHODOLOGY/PRINCIPAL FINDINGS

Using next generation sequencing, we sequenced the genome of four independent Leishmania guyanensis antimony-resistant (SbR) mutants and found different chromosomal alterations including aneuploidy, intrachromosomal gene amplification and gene deletion. A segment covering 30 genes on chromosome 19 was amplified intrachromosomally in three of the four mutants. The gene coding for the multidrug resistance associated protein A involved in antimony resistance was also amplified in the four mutants, most likely through chromosomal translocation. All mutants also displayed a reduced accumulation of antimony mainly due to genomic alterations at the level of the subtelomeric region of chromosome 31 harboring the gene coding for the aquaglyceroporin 1 (LgAQP1). Resistance involved the loss of LgAQP1 through subtelomeric deletions in three mutants. Interestingly, the fourth mutant harbored a single G133D point mutation in LgAQP1 whose role in resistance was functionality confirmed through drug sensitivity and antimony accumulation assays. In contrast to the Leishmania subspecies that resort to extrachromosomal amplification, the Viannia strains studied here used intrachromosomal amplification and locus deletion.

CONCLUSIONS/SIGNIFICANCE

This is the first report of a naturally occurred point mutation in AQP1 in antimony resistant parasites.

Karyotype variability in KP1(+) and KP1(-) strains of Trypanosoma rangeli isolated in Brazil and Colombia

In the present study, the molecular karyotypes of 12 KP1(+) and KP1(-) Trypanosoma rangeli strains were determined and 10 different molecular markers were hybridized to the chromosomes of the parasite, including seven obtained from T. rangeli [ubiquitin hydrolase (UH), a predicted serine/threonine protein kinase (STK), hexose transporter, hypothetical protein, three anonymous sequences] and three from Trypanosoma cruzi [ubiquitin-conjugating enzyme E2 (UBE2), ribosomal RNA methyltransferase (rRNAmtr), proteasome non-ATPase regulatory subunit 6 (PSMD6)]. Despite intraspecific variation, analysis of the karyotype profiles permitted the division of the T. rangeli strains into two groups coinciding with the KP1(+) and KP1(-) genotypes. Southern blot hybridization showed that, except for the hexose transporter probe, all other probes produced distinct patterns able to differentiate the KP1(+) and KP1(-) genotypes. The UH, STK and An-1A04 probes exclusively hybridized to the chromosomes of KP1(+) strains and can be used as markers of this group. In addition, the UBE2, rRNAmtr and PSMD6 markers, which are present in a conserved region in all trypanosomatid species sequenced so far, co-hybridized to the same T. rangeli chromosomal bands, suggesting the occurrence of gene synteny in these species. The finding of distinct molecular karyotypes in KP1(+) and KP1(-) strains of T. rangeli is noteworthy and might be used as a new approach to the study of genetic variability in this parasite. Together with the Southern blot hybridization results, these findings demonstrate that differences at the kDNA level might be associated with variations in nuclear DNA.

Evolutionary consequences of a large duplication event in Trypanosoma brucei: chromosomes 4 and 8 are partial duplicons

Background

Gene order along the genome sequence of the human parasite Trypanosoma brucei provides evidence for a 0.5 Mb duplication, comprising the 3' regions of chromosomes 4 and 8. Here, the principal aim was to examine the contribution made by this duplication event to the T. brucei genome sequence, emphasising the consequences for gene content and the evolutionary change subsequently experienced by paralogous gene copies. The duplicated region may be browsed online at

Results

Comparisons of trypanosomatid genomes demonstrated widespread gene loss from each duplicon, but also showed that 47% of duplicated genes were retained on both chromosomes as paralogous loci. Secreted and surface-expressed genes were over-represented among retained paralogs, reflecting a bias towards important factors at the host-parasite interface, and consistent with a dosage-balance hypothesis. Genetic divergence in both coding and regulatory regions of retained paralogs was bimodal, with a deficit in moderately divergent paralogs; in particular, non-coding sequences were either conserved or entirely remodelled. The conserved paralogs included examples of remarkable sequence conservation, but also considerable divergence of both coding and regulatory regions. Sequence divergence typically displayed strong negative selection; but several features, such as asymmetric evolutionary rates, positively-selected codons and other non-neutral substitutions, suggested that divergence of some paralogs was driven by functional change. The absence of orthologs to retained paralogs in T. congolense indicated that the duplication event was specific to T. brucei.

Conclusion

The duplication of this chromosomal region doubled the dosage of many genes. Rather than creating 'more of the same', these results show that paralogs were structurally modified according to various evolutionary trajectories. The retention of paralogs, and subsequent elaboration of both their primary structures and regulatory regions, strongly suggests that this duplication was a seminal development, stimulating functional innovation and fundamentally altering the genetic repertoire of T. brucei relative to other trypanosomatids.

Identification of a DNA fragment that increases mitotic stability of episomal linear DNAs in Leishmania major

The centromere is a specialized region of eukaryotic chromosomes, the site of kinetochore formation, spindle attachment and regulation of chromosome segregation during mitotic and meiotic cell divisions. To identify sequences which increase mitotic stability and/or act as potential centromeres in Leishmania major, we first generated libraries of Leishmania linear artificial chromosomes (LACs) bearing 30 kb inserts of randomly selected genomic DNAs. These were introduced into parasites, and then their stability was assessed following a period of 10 passages of growth in the absence of selective pressure. Approximately 80% of the 108 transfectants tested lost their LACs promptly and only 20% of the recombinants were retained; of these six showed strong but partial stability (maintained in 30-46% of cells). Mapping and sequencing of one clone (cSC10), which confers the highest degree of maintenance, revealed the presence of a sequence that was found within another stable episome, and which is dispersed in the genome of L. major. The implications of these data to the possible mechanisms of chromosomal maintenance are discussed.

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.

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.

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.

Leishmania major Friedlin chromosome 1 has an unusual distribution of protein-coding genes

Leishmania are evolutionarily ancient protozoans (Kinetoplastidae) and important human pathogens that cause a spectrum of diseases ranging from the asymptomatic to the lethal. The Leishmania genome is relatively small [ approximately 34 megabases (Mb)], lacks substantial repetitive DNA, and is distributed among 36 chromosomes pairs ranging in size from 0.3 Mb to 2.5 Mb, making it a useful candidate for complete genome sequence determination. We report here the nucleotide sequence of the smallest chromosome, chr1. The sequence of chr1 has a 257-kilobase region that is densely packed with 79 protein-coding genes. This region is flanked by telomeric and subtelomeric repetitive elements that vary in number and content among the chr1 homologs, resulting in an approximately 27.5-kilobase size difference. Strikingly, the first 29 genes are all encoded on one DNA strand, whereas the remaining 50 genes are encoded on the opposite strand. Based on the gene density of chr1, we predict a total of approximately 9,800 genes in Leishmania, of which 40% may encode unknown proteins.

Genomics and the biology of parasites

Despite the advances of modern medicine, the threat of chronic illness, disfigurement, or death that can result from parasitic infection still affects the majority of the world population, retarding economic development. For most parasitic diseases, current therapeutics often leave much to be desired in terms of administration regime, toxicity, or effectiveness and potential vaccines are a long way from market. Our best prospects for identifying new targets for drug, vaccine, and diagnostics development and for dissecting the biological basis of drug resistance, antigenic diversity, infectivity and pathology lie in parasite genome analysis, and international mapping and gene discovery initiatives are under way for a variety of protozoan and helminth parasites. These are far from ideal experimental organisms, and the influence of biological and genomic characteristics on experimental approaches is discussed, progress is reviewed and future prospects are examined.

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.

A physical map of the Leishmania major Friedlin genome

An extensive physical map of the Leishmania major Friedlin genome has been assembled by the combination of fingerprint analysis of a shuttle vector cosmid library and probe hybridization. The integrated data obtained for 9004 fingerprinted clones and 974 probes have placed 91.2% of the 33.58-Mb genome into contigs representing each of the 36 chromosomes. This first-generation map has already provided a suitable framework for both high-throughput DNA sequencing and functional studies of the L. major parasite.

Parasite genome analysis. A global map of the Leishmania major genome: prelude to genomic sequencing

In 1994, the World Health Organization (TDR) launched a new strategic initiative in parasite genome analysis, establishing international genome networks for filariae, Schistosoma, Leishmania, Trypanosoma brucei and T. cruzi. For Leishmania, a number of different but complementary approaches have been adopted by members of the Leishmania Genome Network. Our laboratory has been using cosmid clone fingerprinting to produce a physical map of the genome. Progress towards the completion of an integrated physical and biological map of L. major, and the preparations for genomic sequencing, are described.

Selection against the dihydrofolate reductase-thymidylate synthase (DHFR-TS) locus as a probe of genetic alterations in Leishmania major

The genome of the trypanosomatid protozoan genus Leishmania has been shown to undergo a number of changes relevant to drug resistance and virulence, such as gene amplification, chromosomal rearrangement, and variation in ploidy. Experimental approaches to the study of genomic changes have in some cases been limited by the fact that Leishmania cells are asexual diploids, as are some other trypanosomatids, pathogenic fungi, and cultured mammalian cells. Here we report upon a system which permits the measurement of several types of genomic change occurring at the dihydrofolate reductase-thymidylate synthase (DHFR-TS) locus. First, we show that DHFR-TS can function as a positive/negative marker. We used selection against DHFR-TS on a heterozygous line (+/HYG) to generate colonies exhibiting both loss of heterozygosity and structural mutations in DHFR-TS, permitting the first measurement of mutation frequencies in this parasite. Loss of heterozygosity occurred at a frequency ranging from 10(-4) to 10(-6) and was elevated 24-fold by treatment with gamma-irradiation, while the frequency of other events was less than 10(-6) and was increased more than 1,000-fold by nitrosoguanidine treatment. The frequency of loss of heterozygosity relative to other processes such as mutation and gene replacement has important implications for genetic variability in natural Leishmania populations and the generation of both targeted and random mutations. We also developed a protocol for null targeting of diploid cells, in which transfection of a DHFR-TS deletion construct into Leishmania cells followed by negative selection yielded parasites lacking DHFR-TS or foreign sequences. The null-targeting method can be applied to any diploid cell, at any locus for which a negative selection exists. Such marker-free auxotrophic Leishmania cells show potential as an attenuated vaccine, and the methods developed here provide a new approach for manipulating and characterizing the plasticity of the Leishmania genome.

The Leishmania genome comprises 36 chromosomes conserved across widely divergent human pathogenic species

All the physical linkage groups constituting the genome of Leishmania infantum have been identified for the first time by hybridization of specific DNA probes to pulsed field gradient-separated chromosomes. The numerous co-migrating chromosomes were individualised using the distinctive size polymorphisms which occur among strains of the L. infantum/L. donovani complex as a tool. A total of 244 probes, consisting of 41 known genes, 66 expressed sequence tags (ESTs) and 137 anonymous DNA sequences, were assigned to a specific linkage group. We show that this genome comprises 36 chromosomes ranging in size from 0.35 to -3 Mb. This information enabled us to compare the genome structure of L. infantum with those of the three other main Leishmania species that infect man in the Old World, L. major, L. tropica and L. aethiopica. The linkage groups were consistently conserved in all species examined. This result is in striking contrast to the large genetic distances that separate these species and suggests that conservation of the chromosome structure may be critical for this human pathogen. Finally, the high density of markers obtained during the present study (with a mean of 1 marker/130 kb) will speed up the construction of a detailed physical map that would facilitate the genetic analysis of this parasite, for which no classical genetics is available.

Conservation among Old World Leishmania species of six physical linkage groups defined in Leishmania infantum small chromosomes

We have characterised 49 DNA probes specific for each of the six smallest chromosomes in Leishmania infantum and have examined the allocation of these probes in the molecular karyotypes of the other Old World Leishmania species Leishmania donovani, Leishmania major, Leishmania tropica and Leishmania aethiopica. These 49 probes define 6 physical linkage groups in the molecular karyotypes of various strains of L. infantum. 40 of these probes hybridise in the other Old World Leishmania species and show a remarkably conserved linkage pattern. No interchromosomal exchange nor fusion could be detected. Thus, in spite of the chromosomal size polymorphisms, the general structure of the genome seems to be conserved in the six smallest chromosomes among Old World Leishmania species. This structural genomic homogeneity should be helpful for mapping studies of any Old World Leishmania genomes.

Multiple forms of chromosome I, II and V in a restricted population of Leishmania infantum contrasting with monomorphism in individual strains suggest haploidy or automixy

We have resolved the molecular karyotypes of 22 Leishmania infantum strains isolated between 1980 and 1988 in a restricted geographic area and belonging to zymodemes MON-11, -29 and -33. Three strains were isolated from sandflies and all the others from human cutaneous lesions. A high degree of karyotypic homology is observed among these strains, contrasting with the highly polymorphic MON-1 strains isolated in the same area. We have analysed the time-dependent evolution of size variants of chromosomes I to V, each identified by chromosome-specific DNA probes. More evidence is given for the role of subtelomeric regions in chromosomal size variation in Leishmania for both chromosomes I and II. At the population level, the chromosomes I, II and V are present in respectively 8, 4 and 3 distinct sizes. Furthermore, and despite the small size of the sample, various combinations were observed among these different chromosomal forms. These results could be explained by the occurrence of a high rate of recurrent mutations or of genetic exchange. In contrast, only one chromosomal form was observed in individual karyotypes for the chromosomes I-V. These results could tally with the hypothesis of a haploid organisation for these chromosomes and strains, or, in the frame of a diploid organisation, with the hypothesis of a predominantly automictic sexuality giving rise to 2 identical forms of the homologues in the same strain.