The Leishmania donovani complex: genotypes of five metabolic enzymes (ICD, ME, MPI, G6PDH, and FH), new targets for multilocus sequence typing

Flagellates of the Leishmania donovani complex are causative agents of human cutaneous and visceral leishmaniasis. The complex is comprised of L. donovani, Leishmania infantum and Leishmania archibaldi, although the latter is not now considered to be a valid species. Morphological distinction of Leishmania species is impractical, so biochemical, immunological and DNA-based criteria were introduced. Multilocus enzyme electrophoresis (MLEE) is the present gold standard. We have sequenced the genes encoding five metabolic enzymes used for MLEE, both to resolve the DNA diversity underlying isoenzyme mobility differences and to explore the potential of these targets for higher resolution PCR-based multilocus sequence typing. The genes sequenced were isocitrate dehydrogenase, malic enzyme, mannose phosphate isomerase, glucose-6-phosphate dehydrogenase, and fumarate hydratase, for 17 strains of L. infantum, seven strains of L. donovani, and three strains of L. archibaldi. Protein mobilities predicted from amino acid sequences did not always accord precisely with reported MLEE profiles. A high number of heterozygous sites was detected. Heterozygosity was particularly frequent in some strains and indirectly supported the presence of genetic exchange in Leishmania. Phylogenetic analysis of a concatenated alignment based on a total of 263 kb protein-coding sequences showed strong correlation of genotype with geographical origin. Europe and Africa appear to represent independent evolutionary centres.     

Strain typing in Leishmania donovani by using sequence-confirmed amplified region analysis

To differentiate strains of Leishmania donovani, allelic markers at the DNA level were developed by sequence-confirmed amplified region analysis (SCAR). Homologous fragments from different strains of L. donovani were amplified by PCR using random primers and subsequently screened for single-strand conformation polymorphisms. Direct sequencing revealed 55 sequence polymorphisms in eight co-dominant DNA markers; 38 of them were single point mutations. Heterozygosity was evident for 69% and fixed heterozygosity for 25% of all polymorphisms. At most polymorphic sites one of the segregation genotypes was missing. Nineteen unique multilocus genotypes were identified among 29 strains of L. donovani. One genotype was represented by eight Sudanese strains; also two strains from Sudan as well as two strains from Kenya, respectively, shared identical genotypes. All other strains had individual multilocus genotypes. Calculation of genetic distances showed a correlation between multilocus genotypes and the geographical origin of these strains. African strains were found in one well-supported cluster with Kenyan and Sudanese strains clearly separated. SCAR markers seem to represent a random sample of neutral genetic variation present in natural populations. They are co-dominant because they can detect all possible allele combinations in a diploid organism and may, therefore, be very useful for population genetic analysis in Leishmania.     

First report of genetic hybrids between two very divergent Leishmania species: Leishmania infantum and Leishmania major

The genus Leishmania includes many pathogenic species which are genetically very distant. The possibility of genetic exchange between different strains is still an important and debated question. Very few genetic hybrids (i.e., offspring of genetically dissimilar species) have been described in Leishmania. In this study, we report the first example of genetic hybrids occurring between two divergent Leishmania species, Leishmania infantum and Leishmania major. These two species have distinct geographical distributions and are transmitted by different vector species to different mammalian reservoir hosts. These hybrid strains were isolated in Portugal from immunocompromised patients and characterized by molecular and isoenzymatic techniques. These approaches showed that these chimeric strains probably contained the complete genome of both L. major and L. infantum. We believe this is the first report of genetic hybrids between such phylogenetically and epidemiologically distant species of Leishmania. This raises questions about the frequency of such cross-species genetic exchange in natural conditions, modalities of hybrid transmission, their long term maintenance as well as the consequences of these transfers on phenotypes such as drug resistance or pathogenicity.     

Leishmania chagasi/infantum: further investigations on Leishmania tropisms in atypical cutaneous and visceral leishmaniasis foci in Central America

In Central America, apparently genetically identical Leishmania chagasi/infantum parasites cause cutaneous (CL) and visceral leishmaniasis (VL), the latter being more frequent in young children. The present study investigated if there were pathology-related differences in virulence between Honduran CL and VL strains using Mediterranean L. infantum strains as a reference. Macrophage infectivity and serum sensitivity, properties thought to be associated with virulence, were similar between CL and VL strains from both regions. Attention focused on the genome organisation of genes for two candidate virulence factors: Leishmania mitogen activated protein kinase (LMPK) and cysteine proteinase b (Cpb). Interestingly, the Mediterranean strains exhibited restriction enzyme polymorphisms associated with tropism for both LMPK and Cpb genes whereas no differences were observed for the Honduran strains. We also report relative genetic homogeneity of the Honduran strains as compared to the Mediterranean strains and discuss it in terms of the probable origin for the Central American L. chagasi/infantum.     

Multiplex single nucleotide polymorphism genotyping by adapter ligation-mediated allele-specific amplification

An improved approach for increasing the multiplex level of single nucleotide polymorphism (SNP) typing by adapter ligation-mediated allele-specific amplification (ALM-ASA) has been developed. Based on an adapter ligation, each reaction requires n allele-specific primers plus an adapter-specific primer that is common for all SNPs. Thus, only n+1 primers are used for an n-plex PCR amplification. The specificity of ALM-ASA was increased by a special design of the adapter structure and PCR suppression. Given that the genetic polymorphisms in the liver enzyme cytochrome P450 CYP2D6 (debrisoquine 4-hydroxylase) have profound effects on responses of individuals to a particular drug, we selected 17 SNPs in the CYP2D6 gene as an example for the multiplex SNP typing. Without extensive optimization, we successfully typed 17-plex SNPs in the CYP2D6 gene by ALM-ASA. The results for genotyping 70 different genome samples by the 17-plex ALM-ASA were completely consistent with those obtained by both Sanger's sequencing and PCR restriction fragment length polymorphism (PCR-RFLP) analysis. ALM-ASA is a potential method for SNP typing at an ultra-low cost because of a high multiplex level and a simple optimization step for PCR. High-throughput SNP typing could be readily realized by coupling ALM-ASA with a well-developed automation device for sample processing.     

Multilocus sequence and microsatellite identification of intra-specific hybrids and ancestor-like donors among natural Ethiopian isolates of Leishmania donovani

Protozoan parasites of the genus Leishmania (Kinetoplastida: Trypanosomatidae) cause widespread and devastating human diseases. Visceral leishmaniasis is endemic in Ethiopia where it has also been responsible for fatal epidemics. It is postulated that genetic exchange in Leishmania has implications for heterosis (hybrid vigour), spread of virulent strains, resistance to chemotherapeutics, and exploitation of different hosts and vectors. Here we analyse 11 natural Ethiopian Leishmania donovani isolates consisting of four putative hybrids, seven parent-like isolates and over 90 derived biological clones. We apply a novel combination of high resolution multilocus microsatellite typing (five loci) and multilocus sequence typing (four loci) that together distinguish parent-like and hybrid L. donovani strains. Results indicate that the four isolates (and their associated biological clones) are genetic hybrids, not the results of mixed infections, each possessing heterozygous markers consistent with inheritance of divergent alleles from genetically distinct Ethiopian L. donovani lineages. The allelic profiles of the putative hybrids may have arisen from a single hybridisation event followed by inbreeding or gene conversion, or alternatively from two or more hybridisation events. Mitochondrial sequencing showed uniparental maxicircle inheritance for all of the hybrids, each possessing a single mitochondrial genotype. Fluorescence activated cell sorting analysis of DNA content demonstrated that all hybrids and their associated clones were diploid. Together the data imply that intra-specific genetic exchange is a recurrent feature of natural L. donovani populations, with substantial implications for the phyloepidemiology of Leishmania.     

Microsatellite analysis of Toxoplasma gondii shows considerable polymorphism structured into two main clonal groups.

Previous studies on Toxoplasma gondii population structure, based essentially on multilocus restriction fragment length polymorphism analysis or on multilocus enzyme electrophoresis, indicated that T. gondii comprises three clonal lineages. These studies showed a weak polymorphism of the markers (2-4 alleles by locus). In this study, we used eight microsatellite markers to type 84 independent isolates from humans and animals. Two microsatellite markers were present in the introns of two genes, one coding for beta-tubulin and the other for myosin A, and six were found in expressed sequence tags. With 3-16 alleles detected, these markers can be considered as the most discriminating multilocus single-copy markers available for typing T. gondii isolates. This high discriminatory power of microsatellites made it possible to detect mixed infections and epidemiologically related isolates. Evolutionary genetic analyses of diversity show that the T. gondii population structure consists of only two clonal lineages that can be equated to discrete typing units, but there is some evidence of occasional genetic exchange that could explain why one of these discrete typing units is less clearly individualised than the other.     

Plasmodium vivax: microsatellite analysis of multiple-clone infections

We used mixtures of genomic DNA from two genetically distinct isolates from Brazil, 42M and 312M, to investigate how accurately 12-locus microsatellite typing describes the overall genetic diversity and characterizes multilocus haplotypes in multiple-clone Plasmodium vivax infections. We found varying PCR amplification efficiencies of microsatellite alleles; for example, from the same 1:1 mixture of 42M and 312M DNA we amplified predominantly 312M-type alleles at 10 loci and 42M-type alleles at 2 loci. All microsatellite alleles were accurately scored in 1:0.5 and 1:0.25 312M:42M DNA mixtures, even when minor peak heights did not meet previously suggested criteria for minor allele detection in multiple-clone infections. Relative proportions of major and minor alleles were unaffected by multiple displacement amplification of template DNA prior to PCR-based microsatellite typing. Although microsatellite typing may detect minor alleles in clone mixtures, amplification biases may lead to inaccurate assignment of predominant haplotypes in multiple-clone P. vivax infections.     

Population structure of Tunisian Leishmania infantum and evidence for the existence of hybrids and gene flow between genetically different populations

Twenty-seven strains of Leishmania infantum from north and central Tunisia belonging to the three main MON zymodemes (the MON-typing system is based on multilocus enzyme electrophoresis (MLEE) of 15 enzymes) found in this country (MON-1, MON-24 and MON-80) and representing different pathologies (visceral, cutaneous and canine leishmaniasis) have been studied to understand the genetic polymorphism within this species. Intraspecific variation could be detected in L. infantum by the use of 14 hypervariable microsatellite markers. In addition to microsatellite repeat length variation, a high degree of allelic heterozygosity has been observed among the strains investigated, suggestive of sexual recombination within L. infantum groups. The two major clusters found by using Bayesian statistics as well as distance analysis are consistent with the classification based on isoenzymes, dividing Tunisian L. infantum into MON-1 and MON-24/MON-80. Moreover, the existence of hybrid strains between the MON-1 and the non-MON-1 populations has been shown and verified by analysis of clones of one of these strains. Substructure analysis discriminated four groups of L. infantum. The major MON-1 cluster split into two groups, one comprising only Tunisian strains and the second both Tunisian and European strains. The major MON-24 cluster was subdivided into two groups with geographical and clinical feature correlations: a dermotropic group of strains mainly from the north, and a viscerotropic group of strains from the centre of Tunisia. The four viscerotropic hybrid strains all originated from central Tunisia and were typed by MLEE as MON-24 or MON-80. To our knowledge, this is the first report describing relationships between clinical picture and population substructure of L. infantum MON-24 based on genotype data, as well as the existence of hybrids between zymodemes MON-1 and MON-24/MON-80, and proving one of these hybrid strains by molecular analysis of the parent strain and its clones.     

Comparison of four methods, including semi-automated rep-PCR, for the typing of vancomycin-resistant Enterococcus faecium

We have assessed the performance of semi-automated rep-PCR (Diversilab®) and multilocus sequence typing (MLST) in comparison to pulsed-field gel electrophoresis (PFGE) for typing a collection of 29 epidemiologically characterized vancomycin-resistant Enterococcus faecium (VRE). Sixteen strains that harbored the Tn1546 element were typed by PCR mapping. The discriminative power of the typing methods was calculated by the Simpson's index of diversity, and the concordance between methods was evaluated by the Kendall's coefficient of concordance. Semi-automated rep-PCR appeared as discriminative as PFGE and was further compared with PFGE for typing 67 VRE isolated during a hospital outbreak. Rep-PCR appeared to be more discriminative than PFGE for this second set of strains. Reproducibility of DiversiLab® was also tested against 35 selected isolates. Only three showed less than 97% similarity, indicating high reproducibility at this level of discrimination. In conclusion, semi-automated rep-PCR is a useful tool for rapid screening of VRE isolates during an outbreak, although cost of the system may be limiting for routine implementation. PFGE, which remains the reference method, should be used for confirmation and evaluation of the genetic relatedness of epidemic isolates.     

Multilocus sequence typing of Salmonella strains by high-throughput sequencing of selectively amplified target genes

Rapid development of next generation sequencing (NGS) technologies in recent years has made whole genome sequencing of bacterial genomes widely accessible. However, it is often unnecessary or not feasible to sequence the whole genome for most applications of genetic analyses in bacteria. Selectively capturing defined genomic regions followed by NGS analysis could be a promising approach for high-resolution molecular typing of a large set of strains. In this study, we describe a novel and straightforward PCR-based target-capturing method, hairpin-primed multiplex amplification (HPMA), which allows for simultaneous amplification of numerous target genes. To test the feasibility of NGS-based strain typing using HPMA, 20 target gene sequences were simultaneously amplified with barcode tagging in each of 41 Salmonella strains. The amplicons were then pooled and analyzed by 454 pyrosequencing. Analysis of the sequence data, as an extension of multilocus sequence typing (MLST), demonstrated the utility and potential of this novel typing method, MLST-seq, as a high-resolution strain typing method. With the rapidly increasing sequencing capacity of NGS, MLST-seq or its variations using different target enrichment methods can be expected to become a high-resolution typing method in the near future for high-throughput analysis of a large collection of bacterial strains.     

Genetic diversity and population structure of Escherichia coli isolated from freshwater beaches

Escherichia coli is an important member of the gastrointestinal tract of humans and warm-blooded animals (primary habitat). In the external environment outside the host (secondary habitat), it is often considered to be only a transient member of the microbiota found in water and soil, although recent evidence suggests that some strains can persist in temperate soils and freshwater beaches. Here we quantified the population genetic structure of E. coli from a longitudinal collection of environmental strains isolated from six freshwater beaches along Lake Huron and the St. Clair River in Michigan. Multilocus enzyme electrophoresis (MLEE) and multilocus sequence typing (MLST) revealed extensive genetic diversity among 185 E. coli isolates with an average of 40 alleles per locus. Despite evidence for extensive recombination generating new alleles and genotypic diversity, several genotypes marked by distinct MLEE and MLST profiles were repeatedly recovered from separate sites at different times. A PCR-based phylogrouping technique showed that the persistent, naturalized E. coli belonged to the B1 group. These results support the hypothesis that persistent genotypes have an adaptive advantage in the secondary habitat outside the host.     

Biochemical genetic relationship of Thailand and Hawaii isolates of Parastrongylus cantonensis (Nematoda: Angiostrongylidae)

The genetic relationship of the Thailand and Hawaii isolates (strains) of the rat lungworm Parastrongylus (=Angiostrongylus) cantonensis was investigated by gene–enzyme systems using vertical polyacrylamide slab gel electrophoresis. Six gene–enzyme systems were successfully determined, with each being represented by two presumptive loci. Glucose phosphate dehydrogenase, glucose phosphate isomerase, lactate dehydrogenase, malate dehydrogenase and malic enzyme were monomorphic at both loci, and the respective bands exhibited similar mobility in both isolates implying absence of genetic variation.     

New Insights on Taxonomy, Phylogeny and Population Genetics of Leishmania (Viannia) Parasites Based on Multilocus Sequence Analysis

The Leishmania genus comprises up to 35 species, some with status still under discussion. The multilocus sequence typing (MLST)—extensively used for bacteria—has been proposed for pathogenic trypanosomatids. For Leishmania, however, a detailed analysis and revision on the taxonomy is still required. We have partially sequenced four housekeeping genes—glucose-6-phosphate dehydrogenase (G6PD), 6-phosphogluconate dehydrogenase (6PGD), mannose phosphate isomerase (MPI) and isocitrate dehydrogenase (ICD)—from 96 Leishmania (Viannia) strains and assessed their discriminatory typing capacity. The fragments had different degrees of diversity, and are thus suitable to be used in combination for intra- and inter-specific inferences. Species-specific single nucleotide polymorphisms were detected, but not for all species; ambiguous sites indicating heterozygosis were observed, as well as the putative homozygous donor. A large number of haplotypes were detected for each marker; for 6PGD a possible ancestral allele for L. (Viannia) was found. Maximum parsimony-based haplotype networks were built. Strains of different species, as identified by multilocus enzyme electrophoresis (MLEE), formed separated clusters in each network, with exceptions. NeighborNet of concatenated sequences confirmed species-specific clusters, suggesting recombination occurring in L. braziliensis and L. guyanensis. Phylogenetic analysis indicates L. lainsoni and L. naiffi as the most divergent species and does not support L. shawi as a distinct species, placing it in the L. guyanensis cluster. BURST analysis resulted in six clonal complexes (CC), corresponding to distinct species. The L. braziliensis strains evaluated correspond to one widely geographically distributed CC and another restricted to one endemic area. This study demonstrates the value of systematic multilocus sequence analysis (MLSA) for determining intra- and inter-species relationships and presents an approach to validate the species status of some entities. Furthermore, it contributes to the phylogeny of L. (Viannia) and might be helpful for epidemiological and population genetics analysis based on haplotype/diplotype determinations and inferences.     

Electrophoretic Transfer from Polyacrylamide Gel to Nitrocellulose Sheets, a New Method To Characterize Multilocus Enzyme Genotypes of Klebsiella Strains

A new method for multilocus enzyme electrophoresis, based on electrophoretic transfers to nitrocellulose after polyacrylamide-agarose gel electrophoresis was explored. Electrophoretic separation was performed on 1-mm-thick slab gels with 6-μl samples of bacterial extracts and was followed by serial 5-min consecutive transfers. The transferability of 19 metabolic enzymes of Klebsiella strains was studied and allowed the simultaneous examination of one enzyme in the separation gel and at least five enzymes on nitrocellulose sheets. The resolution of enzyme bands was increased on nitrocellulose; thus, well-separated bands were recorded for nucleoside phosphorylase, peptidase, and phosphoglucose isomerase whereas their mobility variants could not be clearly distinguished in the separation gel because of stain diffusion. The study of genetic relationships of 42 strains of Klebsiella pneumoniae and 24 strains of Klebsiella oxytoca demonstrated the reliability of the method, since clustering analysis of electrophoretic types, based on electrophoretic polymorphism of 10 metabolic enzymes, showed two main clusters well correlated with the two species. The 57 electrophoretic types described confirm the usefulness of the method for the study of genetic relationships between closely related strains.     

Leishmania donovani: Genetic diversity of isolates from Sudan characterized by PCR-based RAPD

Drug unresponsiveness in patients with visceral leishmaniasis (VL) is a problem in many endemic areas. This study aimed to determine genetic diversity of Leishmania donovani isolates from a VL endemic area in Sudan as a possible explanation for drug unresponsiveness in some patients. Thirty clinically stibogluconate (SSG)-sensitive isolates were made SSG-unresponsive in vitro by gradually increasing SSG concentrations. The sensitive isolates and their SSG-unresponsive counterparts were typed using mini-circle kDNA and categorized using PCR-RAPD. All the isolates were typed as L. donovani, the resulting PCR-RAPD characterization of the SSG-sensitive isolates gave three distinct primary genotypes while, the SSG-unresponsive isolates showed only a single band. L. donovani isolates from eastern Sudan are diverse; this probably resulted from emergence of new L. donovani strains during epidemics due to the pressure of widespread use of antimonials.In this communication the possible role of isolates diversity in antimonial unresponsiveness and the in vitro changing PCR-RAPD band pattern in SSG-unresponsive strains were discussed.     

Extensive intragenic recombination and patterns of linkage disequilibrium at the CSN3 locus in European rabbit

Kappa-casein (CSN3) plays an important role in stabilising the Ca-sensitive caseins in the micelle. The European rabbit (Oryctolagus cuniculus) CSN3 has previously been shown to possess two alleles (A and B), which differ deeply in their intronic regions (indels of 100 and 1550 nucleotides in introns 1 and 4, respectively). Furthermore, a correlation between several reproductive performance traits and the different alleles was described. However, all these data were exclusively collected in rabbit domestic breeds, preventing a deeper understanding of the extensive polymorphism observed in the CSN3 gene. Additionally, the techniques available for the typing of both indel polymorphisms were until now not suitable for large-scale studies. In this report, we describe a simple, PCR-based typing method to distinguish rabbit CSN3 alleles. We analyse both ancient wild rabbit populations from the Iberian Peninsula and France, and the more recently derived English wild rabbits and domestic stocks. A new allele (C) showing another major indel (250 bp) in intron 1 was found, but exclusively detected in Iberian wild rabbits. In addition, our survey revealed the occurrence of new haplotypes in wild populations, suggesting that intragenic recombination is important in creating genetic diversity at this locus. This easy and low cost single-step PCR-based method results in an improvement over previous described techniques, can be easily set up in a routine molecular laboratory and would probably be a valuable tool in the management of rabbit domestic breeds.     

Discrimination of Bacillus anthracis and closely related microorganisms by analysis of 16S and 23S rRNA with oligonucleotide microarray

Analysis of 16S rRNA sequences is a commonly used method for the identification and discrimination of microorganisms. However, the high similarity of 16S and 23S rRNA sequences of Bacillus cereus group organisms (up to 99-100%) and repeatedly failed attempts to develop molecular typing systems that would use DNA sequences to discriminate between species within this group have resulted in several suggestions to consider B. cereus and B. thuringiensis, or these two species together with B. anthracis, as one species. Recently, we divided the B. cereus group into seven subgroups, Anthracis, Cereus A and B, Thuringiensis A and B, and Mycoides A and B, based on 16S rRNA, 23S rRNA and gyrB gene sequences and identified subgroup-specific makers in each of these three genes. Here we for the first time demonstrated discrimination of these seven subgroups, including subgroup Anthracis, with a 3D gel element microarray of oligonucleotide probes targeting 16S and 23S rRNA markers. This is the first microarray enabled identification of B. anthracis and discrimination of these seven subgroups in pure cell cultures and in environmental samples using rRNA sequences. The microarray bearing perfect match/mismatch (p/mm) probe pairs was specific enough to discriminate single nucleotide polymorphisms (SNPs) and was able to identify targeted organisms in 5min. We also demonstrated the ability of the microarray to determine subgroup affiliations for B. cereus group isolates without rRNA sequencing. Correlation of these seven subgroups with groupings based on multilocus sequence typing (MLST), fluorescent amplified fragment length polymorphism analysis (AFLP) and multilocus enzyme electrophoresis (MME) analysis of a wide spectrum of different genes, and the demonstration of subgroup-specific differences in toxin profiles, psychrotolerance, and the ability to harbor some plasmids, suggest that these seven subgroups are not based solely on neutral genomic polymorphisms, but instead reflect differences in both the genotypes and phenotypes of the B. cereus group organisms.