The Mating Competence of Geographically Diverse Leishmania major Strains in Their Natural and Unnatural Sand Fly Vectors

Detecting complex infections in trypanosomatids using whole genome sequencing

BACKGROUND

Trypanosomatid parasites are a group of protozoans that cause devastating diseases that disproportionately affect developing countries. These protozoans have developed several mechanisms for adaptation to survive in the mammalian host, such as extensive expansion of multigene families enrolled in host-parasite interaction, adaptation to invade and modulate host cells, and the presence of aneuploidy and polyploidy. Two mechanisms might result in "complex" isolates, with more than two haplotypes being present in a single sample: multiplicity of infections (MOI) and polyploidy. We have developed and validated a methodology to identify multiclonal infections and polyploidy using whole genome sequencing reads, based on fluctuations in allelic read depth in heterozygous positions, which can be easily implemented in experiments sequencing genomes from one sample to larger population surveys.

RESULTS

The methodology estimates the complexity index (CI) of an isolate, and compares real samples with simulated clonal infections at individual and populational level, excluding regions with somy and gene copy number variation. It was primarily validated with simulated MOI and known polyploid isolates respectively from Leishmania and Trypanosoma cruzi. Then, the approach was used to assess the complexity of infection using genome wide SNP data from 497 trypanosomatid samples from four clades, L. donovani/L. infantum, L. braziliensis, T. cruzi and T. brucei providing an overview of multiclonal infection and polyploidy in these cultured parasites. We show that our method robustly detects complex infections in samples with at least 25x coverage, 100 heterozygous SNPs and where 5-10% of the reads correspond to the secondary clone. We find that relatively small proportions (≤ 7%) of cultured trypanosomatid isolates are complex.

CONCLUSIONS

The method can accurately identify polyploid isolates, and can identify multiclonal infections in scenarios with sufficient genome read coverage. We pack our method in a single R script that requires only a standard variant call format (VCF) file to run ( https://github.com/jaumlrc/Complex-Infections ). Our analyses indicate that multiclonality and polyploidy do occur in all clades, but not very frequently in cultured trypanosomatids. We caution that our estimates are lower bounds due to the limitations of current laboratory and bioinformatic methods.

Comparative genomics of Leishmania donovani progeny from genetic crosses in two sand fly species and impact on the diversity of diagnostic and vaccine candidates

Sand fly transmitted Leishmania species are responsible for severe, wide ranging, visceral and cutaneous leishmaniases. Genetic exchange can occur among natural Leishmania populations and hybrids can now be produced experimentally, with limitations. Feeding Phlebotomus orientalis or Phlebotomus argentipes on two strains of Leishmania donovani yielded hybrid progeny, selected using double drug resistance and fluorescence markers. Fluorescence activated cell sorting of cultured clones derived from these hybrids indicated diploid progeny. Multilocus sequence typing of the clones showed hybridisation and nuclear heterozygosity, although with inheritance of single haplotypes in a kinetoplastid target. Comparative genomics showed diversity of clonal progeny between single chromosomes, and extraordinary heterozygosity across all 36 chromosomes. Diversity between progeny was seen for the HASPB antigen, which has been noted previously as having implications for design of a therapeutic vaccine. Genomic diversity seen among Leishmania strains and hybrid progeny is of great importance in understanding the epidemiology and control of leishmaniasis. As an outcome of this study we strongly recommend that wider biological archives of different Leishmania species from endemic regions should be established and made available for comparative genomics. However, in parallel, performance of genetic crosses and genomic comparisons should give fundamental insight into the specificity, diversity and limitations of candidate diagnostics, vaccines and drugs, for targeted control of leishmaniasis.

HOP1 and HAP2 are conserved components of the meiosis-related machinery required for successful mating in Leishmania

Whole genome analysis of Leishmania hybrids generated experimentally in sand flies supports a meiotic mechanism of genetic exchange, with Mendelian segregation of the nuclear genome. Here, we perform functional analyses through the generation of double drug-resistant hybrids in vitro and in vivo (during sand fly infections) to assess the importance of conserved meiosis-related genes in recombination and plasmogamy. We report that HOP1 and a HAP2-paralog (HAP2-2) are essential components of the Leishmania meiosis machinery and cell-to-cell fusion mechanism, respectively, since deletion of either gene in one or both parents significantly reduces or completely abrogates mating competence. These findings significantly advance our understanding of sexual reproduction in Leishmania, with likely relevance to other trypanosomatids, by formally demonstrating the involvement of a meiotic protein homolog and a distinct fusogen that mediates non-canonical, bilateral fusion in the hybridizing cells.

Leishmania genetic exchange is mediated by IgM natural antibodies

Host factors that mediate Leishmania genetic exchange are not well defined. Here we demonstrate that natural IgM (IgMn)1-4 antibodies mediate parasite genetic exchange by inducing the transient formation of a spherical parasite clump that promotes parasite fusion and hybrid formation. We establish that IgMn from Leishmania-free animals binds to the surface of Leishmania parasites to induce significant changes in the expression of parasite transcripts and proteins. Leishmania binding to IgMn is partially lost after glycosidase treatment, although parasite surface phosphoglycans, including lipophosphoglycan, are not required for IgMn-induced parasite clumping. Notably, the transient formation of parasite clumps is essential for Leishmania hybridization in vitro. In vivo, we observed a 12-fold increase in hybrid formation in sand flies provided a second blood meal containing IgMn compared with controls. Furthermore, the generation of recombinant progeny from mating hybrids and parental lines were only observed in sand flies provided with IgMn. Both in vitro and in vivo IgM-induced Leishmania crosses resulted in full genome hybrids that show equal patterns of biparental contribution. Leishmania co-option of a host natural antibody to facilitate mating in the insect vector establishes a new paradigm of parasite-host-vector interdependence that contributes to parasite diversity and fitness by promoting genetic exchange.

First evidence of experimental genetic hybridization between cutaneous and visceral strains of Leishmania donovani within its natural vector Phlebotomus argentipes

Leishmaniasis is a neglected tropical disease caused by protozoan parasites of genus Leishmania, and transmitted by different species of Phlebotomine sand flies. More than 20 species of Leishmania are known to cause disease in humans and other animals. Leishmania donovani species complex is known to have a vast diversity of clinical manifestations in humans, but underlying mechanisms for such diversity are yet unknown. Long believed to be strictly asexual, Leishmania have been shown to undergo a cryptic sexual cycle inside its sandfly vector. Natural populations of hybrid parasites have been associated with the rise of atypical clinical outcomes in the Indian subcontinent (ISC). However, formal demonstration of genetic crossing in the major endemic sandfly species in the ISC remain unexplored. Here, we investigated the ability of two distinct variants of L. donovani associated with strikingly different forms of the disease to undergo genetic exchange inside its natural vector, Phlebotomus argentipes. Clinical isolates of L. donovani either from a Sri Lankan cutaneous leishmaniasis (CL) patient or an Indian visceral leishmaniasis (VL) patient were genetically engineered to express different fluorescent proteins and drug-resistance markers and subsequently used as parental strains in experimental sandfly co-infection. After 8 days of infection, sand flies were dissected and midgut promastigotes were transferred into double drug-selective media. Two double drug-resistant, dual fluorescent hybrid cell lines were recovered, which after cloning and whole genome sequencing, were shown to be full genomic hybrids. This study provides the first evidence of L. donovani hybridization within its natural vector Ph. argentipes.

Natural infection with Leishmania (Mundinia) martiniquensis supports Culicoides peregrinus (Diptera: Ceratopogonidae) as a potential vector of leishmaniasis and characterization of a Crithidia sp. isolated from the midges

The prevalence of autochthonous leishmaniasis in Thailand is increasing but the natural vectors that are responsible for transmission remain unknown. Experimental in vivo infections in Culicoides spp. with Leishmania (Mundinia) martiniquensis and Leishmania (Mundinia) orientalis, the major causative pathogens in Thailand, have demonstrated that biting midges can act as competent vectors. Therefore, the isolation and detection of Leishmania and other trypanosomatids were performed in biting midges collected at a field site in an endemic area of leishmaniasis in Tha Ruea and a mixed farm of chickens, goats, and cattle in Khuan Phang, Nakhon Si Thammarat province, southern Thailand. Results showed that Culicoides peregrinus was the abundant species (>84%) found in both locations and only cow blood DNA was detected in engorged females. Microscopic examination revealed various forms of Leishmania promastigotes in the foregut of several C. peregrinus in the absence of bloodmeal remnants, indicating established infections. Molecular identification using ITS1 and 3'UTR HSP70 type I markers showed that the Leishmania parasites found in the midges were L. martiniquensis. The infection rate of L. martiniquensis in the collected flies was 2% in Tha Ruea and 6% in Khuan Phang, but no L. orientalis DNA or parasites were found. Additionally, organisms from two different clades of Crithidia, both possibly new species, were identified using SSU rRNA and gGAPDH genes. Choanomastigotes and promastigotes of both Crithidia spp. were observed in the hindgut of the dissected C. peregrinus. Interestingly, midges infected with both L. martiniquensis and Crithidia were found. Moreover, four strains of Crithidia from one of the clades were successfully isolated into culture. These parasites could grow at 37°C in the culture and infect BALB/c mice macrophages but no multiplication was observed, suggesting they are thermotolerant monoxenous trypanosomatids similar to Cr. thermophila. These findings provide the first evidence of natural infection of L. martiniquensis in C. peregrinus supporting it as a potential vector of L. martiniquensis.

Ethanolaminephosphate cytidylyltransferase is essential for survival, lipid homeostasis and stress tolerance in Leishmania major

Glycerophospholipids including phosphatidylethanolamine (PE) and phosphatidylcholine (PC) are vital components of biological membranes. Trypanosomatid parasites of the genus Leishmania can acquire PE and PC via de novo synthesis and the uptake/remodeling of host lipids. In this study, we investigated the ethanolaminephosphate cytidylyltransferase (EPCT) in Leishmania major, which is the causative agent for cutaneous leishmaniasis. EPCT is a key enzyme in the ethanolamine branch of the Kennedy pathway which is responsible for the de novo synthesis of PE. Our results demonstrate that L. major EPCT is a cytosolic protein capable of catalyzing the formation of CDP-ethanolamine from ethanolamine-phosphate and cytidine triphosphate. Genetic manipulation experiments indicate that EPCT is essential in both the promastigote and amastigote stages of L. major as the chromosomal null mutants cannot survive without the episomal expression of EPCT. This differs from our previous findings on the choline branch of the Kennedy pathway (responsible for PC synthesis) which is required only in promastigotes but not amastigotes. While episomal EPCT expression does not affect promastigote proliferation under normal conditions, it leads to reduced production of ethanolamine plasmalogen or plasmenylethanolamine, the dominant PE subtype in Leishmania. In addition, parasites with episomal EPCT exhibit heightened sensitivity to acidic pH and starvation stress, and significant reduction in virulence. In summary, our investigation demonstrates that proper regulation of EPCT expression is crucial for PE synthesis, stress response, and survival of Leishmania parasites throughout their life cycle.

Leishmania allelic selection during experimental sand fly infection correlates with mutational signatures of oxidative DNA damage

Trypanosomatid pathogens are transmitted by blood-feeding insects, causing devastating human infections. These parasites show important phenotypic shifts that often impact parasite pathogenicity, tissue tropism, or drug susceptibility. The evolutionary mechanisms that allow for the selection of such adaptive phenotypes remain only poorly investigated. Here, we use Leishmania donovani as a trypanosomatid model pathogen to assess parasite evolutionary adaptation during experimental sand fly infection. Comparing the genome of the parasites before and after sand fly infection revealed a strong population bottleneck effect as judged by allele frequency analysis. Apart from random genetic drift caused by the bottleneck effect, our analyses revealed haplotype and allelic changes during sand fly infection that seem under natural selection given their convergence between independent biological replicates. Our analyses further uncovered signature mutations of oxidative DNA damage in the parasite genomes after sand fly infection, suggesting that Leishmania suffers from oxidative stress inside the insect digestive tract. Our results propose a model of Leishmania genomic adaptation during sand fly infection, with oxidative DNA damage and DNA repair processes likely driving haplotype and allelic selection. The experimental and computational framework presented here provides a useful blueprint to assess evolutionary adaptation of other eukaryotic pathogens inside their insect vectors, such as Plasmodium spp, Trypanosoma brucei, and Trypanosoma cruzi.

Ethanolaminephosphate cytidyltransferase is essential for survival, lipid homeostasis and stress tolerance in Leishmania major

Glycerophospholipids including phosphatidylethanolamine (PE) and phosphatidylcholine (PC) are vital components of biological membranes. Trypanosomatid parasites of the genus Leishmania can acquire PE and PC via de novo synthesis and the uptake/remodeling of host lipids. In this study, we investigated the ethanolaminephosphate cytidyltransferase (EPCT) in Leishmania major , which is the causative agent for cutaneous leishmaniasis. EPCT is a key enzyme in the ethanolamine branch of the Kennedy pathway which is responsible for the de novo synthesis of PE. Our results demonstrate that L. major EPCT is a cytosolic protein capable of catalyzing the formation of CDP-ethanolamine from ethanolamine-phosphate and cytidine triphosphate. Genetic manipulation experiments indicate that EPCT is essential in both the promastigote and amastigote stages of L. major as the chromosomal null mutants cannot survive without the episomal expression of EPCT. This differs from our previous findings on the choline branch of the Kennedy pathway (responsible for PC synthesis) which is required only in promastigotes but not amastigotes. While episomal EPCT expression does not affect promastigote proliferation under normal conditions, it leads to reduced production of ethanolamine plasmalogen or plasmenylethanolamine, the dominant PE subtype in Leishmania . In addition, parasites with epsiomal EPCT exhibit heightened sensitivity to acidic pH and starvation stress, and significant reduction in virulence. In summary, our investigation demonstrates that proper regulation of EPCT expression is crucial for PE synthesis, stress response, and survival of Leishmania parasites throughout their life cycle.

AUTHOR SUMMARY

In nature, Leishmania parasites alternate between fast replicating, extracellular promastigotes in sand fly gut and slow growing, intracellular amastigotes in macrophages. Previous studies suggest that promastigotes acquire most of their lipids via de novo synthesis whereas amastigotes rely on the uptake and remodeling of host lipids. Here we investigated the function of ethanolaminephosphate cytidyltransferase (EPCT) which catalyzes a key step in the de novo synthesis of phosphatidylethanolamine (PE) in Leishmania major . Results showed that EPCT is indispensable for both promastigotes and amastigotes, indicating that de novo PE synthesis is still needed at certain capacity for the intracellular form of Leishmania parasites. In addition, elevated EPCT expression alters overall PE synthesis and compromises parasite’s tolerance to adverse conditions and is deleterious to the growth of intracellular amastigotes. These findings provide new insight into how Leishmania acquire essential phospholipids and how disturbance of lipid metabolism can impact parasite fitness.

Self-Hybridization in Leishmania major

Genetic exchange between different Leishmania strains in the sand fly vector has been experimentally demonstrated and is supported by population genetic studies. In nature, opportunities for Leishmania interstrain mating are restricted to flies biting multiply infected hosts or through multiple bites of different hosts. In contrast, self-mating could occur in any infected sand fly. By crossing two recombinant lines derived from the same Leishmania major strain, each expressing a different drug-resistance marker, self-hybridization in L. major was confirmed in a natural sand fly vector, Phlebotomus duboscqi, and in frequencies comparable to interstrain crosses. We provide the first high resolution, whole-genome sequencing analysis of large numbers of selfing progeny, their parents, and parental subclones. Genetic exchange consistent with classical meiosis is supported by the biallelic inheritance of the rare homozygous single nucleotide polymorphisms (SNPs) that arose by mutation during the generation of the parental clones. In contrast, heterozygous SNPs largely failed to be transmitted in Mendelian ratios for reasons not understood. SNPs that were heterozygous in both parents, however, recombined to produce homozygous alleles in some hybrids. For trisomic chromosomes present in both parents, transmittal to the progeny was only altered by self-hybridization, involving a gain or loss of somy in frequencies predicted by a meiotic process. Whole-genome polyploidization was also observed in the selfing progeny. Thus, self-hybridization in Leishmania, with its potential to occur in any infected sand fly, may be an important source of karyotype variation, loss of heterozygosity, and functional diversity. IMPORTANCE Leishmania are parasitic protozoa that cause a wide spectrum of diseases collectively known as the leishmaniases. Sexual reproduction in Leishmania has been proposed as an important source of genetic diversity and has been formally demonstrated to occur inside the sand fly vector midgut. Nevertheless, in the wild, opportunities for genetic exchange between different Leishmania species or strains are restricted by the capacity of different Leishmania strains to colonize the same sand fly. In this work, we report the first high resolution, whole-genome sequence analysis of intraclonal genetic exchange as a type of self-mating in Leishmania. Our data reveal that self-hybridization can occur with comparable frequency as interstrain mating under experimental lab conditions, leading to important genomic alterations that can potentially take place within every naturally infected sand fly.

Impact of collaborative actions of Leishmania (Viannia) braziliensis subpopulations on the infection profile

This study focuses on the role of the population structure of Leishmania spp. on the adaptive capacity of the parasite. Herein, we investigate the contribution of subpopulations of the L. (V.) braziliensis Thor strain (Thor03, Thor10 and Thor22) in the profile of murine macrophages infection. Infection assays were performed with binary combinations of these subpopulations at stationary phases. The initial interaction time showed major effects on the combination assays, as demonstrated by the significant increase in the infection rate at 5 h. Based on the endocytic index (EI), Thor10 (EI = 563.6) and Thor03 (EI = 497) showed a higher infection load compared to Thor22 (EI = 227.3). However, the EI decreased in Thor03 after 48 h (EI = 447) and 72 h (EI = 388.3) of infection, and showed changes in the infection level in all Thor10/Thor22 combinations. Assays with CellTrace CFSE-labelled Thor22 promastigotes indicated an increase (~1.5 fold) in infection by this subpopulation in the presence of Thor10 when compared to the infection profile of Thor03/Thor22 combinations in the same proportions. In addition, the potential of these subpopulations, alone or in binary combinations, to modulate the expression of cytokines and nitric oxide (NO) in vitro was investigated. Lower NO and tumour necrosis factor-α production levels were observed for all Thor10/Thor22 combinations at 24 h compared to these subpopulations alone. In contrast, Thor03/Thor22 combination assays increased IL-10 production at this time. Collectively, these results provide in vitro evidence on the potential of L. (V.) braziliensis population structure to play a relevant role in a host infection by this parasite.

Experimental Hybridization in Leishmania: Tools for the Study of Genetic Exchange

Despite major advances over the last decade in our understanding of Leishmania reproductive strategies, the sexual cycle in Leishmania has defied direct observation and remains poorly investigated due to experimental constraints. Here, we summarize the findings and conclusions drawn from genetic analysis of experimental hybrids generated in sand flies and highlight the recent advances in generating hybrids in vitro. The ability to hybridize between culture forms of different species and strains of Leishmania should invite more intensive investigation of the mechanisms underlying genetic exchange and provide a rich source of recombinant parasites for future genetic analyses.

Effect of hybridization on Lipophosphoglycan expression in Leishmania major

Leishmania major is the causative agent of cutaneous leishmaniasis (CL). It is one of the most studied Leishmania species not only during vector interaction, but also in the vertebrate host. Lipophosphoglycan (LPG) is the Leishmania multifunctional virulence factor during host-parasite interaction, whose polymorphisms are involved in the immunopathology of leishmaniasis. Although natural hybrids occur in nature, hybridization of L. major strains in the laboratory was successfully demonstrated. However, LPG expression in the hybrids remains unknown. LPGs from parental (Friedlin-Fn and Seidman-Sd) and hybrids (FnSd3, FnSd4A, FnSd4B and FnSd6F) were extracted, purified and their repeat units analyzed by immunoblotting and fluorophore-assisted carbohydrate electrophoresis (FACE). Parental strains have distinct profiles in LPG expression, and a mixed profile was observed for all hybrids. Variable levels of NO production by macrophages were detected after LPG exposure (parental and hybrids) and were strain specific. This article is protected by copyright. All rights reserved.

Signatures of hybridization in Trypanosoma brucei

Genetic exchange among disease-causing micro-organisms can generate progeny that combine different pathogenic traits. Though sexual reproduction has been described in trypanosomes, its impact on the epidemiology of Human African Trypanosomiasis (HAT) remains controversial. However, human infective and non-human infective strains of Trypanosoma brucei circulate in the same transmission cycles in HAT endemic areas in subsaharan Africa, providing the opportunity for mating during the developmental cycle in the tsetse fly vector. Here we investigated inheritance among progeny from a laboratory cross of T. brucei and then applied these insights to genomic analysis of field-collected isolates to identify signatures of past genetic exchange. Genomes of two parental and four hybrid progeny clones with a range of DNA contents were assembled and analysed by k-mer and single nucleotide polymorphism (SNP) frequencies to determine heterozygosity and chromosomal inheritance. Variant surface glycoprotein (VSG) genes and kinetoplast (mitochondrial) DNA maxi- and minicircles were extracted from each genome to examine how each of these components was inherited in the hybrid progeny. The same bioinformatic approaches were applied to an additional 37 genomes representing the diversity of T. brucei in subsaharan Africa and T. evansi. SNP analysis provided evidence of crossover events affecting all 11 pairs of megabase chromosomes and demonstrated that polyploid hybrids were formed post-meiotically and not by fusion of the parental diploid cells. VSGs and kinetoplast DNA minicircles were inherited biparentally, with approximately equal numbers from each parent, whereas maxicircles were inherited uniparentally. Extrapolation of these findings to field isolates allowed us to distinguish clonal descent from hybridization by comparing maxicircle genotype to VSG and minicircle repertoires. Discordance between maxicircle genotype and VSG and minicircle repertoires indicated inter-lineage hybridization. Significantly, some of the hybridization events we identified involved human infective and non-human infective trypanosomes circulating in the same geographic areas.

Axenic interspecies and intraclonal hybrid formation in Leishmania: Successful crossings between visceral and cutaneous strains

Diseases caused by trypanosomatids are serious public health concerns in low-income endemic countries. Leishmaniasis is presented in two main clinical forms, visceral leishmaniasis-caused by L. infantum and L. donovani-and cutaneous leishmaniasis-caused by many species, including L. major, L. tropica and L. braziliensis. As for certain other trypanosomatids, sexual reproduction has been confirmed in these parasites, and formation of hybrids can contribute to virulence, drug resistance or adaptation to the host immune system. In the present work, the capability of intraclonal and interspecies genetic exchange has been investigated using three parental strains: L. donovani, L. tropica and L. major, which have been engineered to express different fluorescent proteins and antibiotic resistance markers in order to facilitate the phenotypic selection of hybrid parasites after mating events. Stationary and exponential-phase promastigotes of each species were used, in in vitro experiments, some of them containing LULO cells (an embryonic cell line derived from Lutzomyia longipalpis). Several intraclonal hybrids were obtained with L. tropica as crossing progenitor, but not with L. donovani or L. major. In interspecies crossings, three L. donovani x L. major hybrids and two L. donovani x L. tropica hybrids were isolated, thereby demonstrating the feasibility to obtain in vitro hybrids of parental lines causing different tropism of leishmaniasis. Ploidy analysis revealed an increase in DNA content in all hybrids compared to the parental strains, and nuclear analysis showed that interspecies hybrids are complete hybrids, i.e. each of them showing at least one chromosomal set from each parental. Regarding kDNA inheritance, discrepancies were observed between maxi and minicircle heritage. Finally, phenotypic studies showed either intermediate phenotypes in terms of growth profiles, or a decreased in vitro infection capacity compared to the parental cells. To the best of our knowledge, this is the first time that in vitro interspecies outcrossing has been demonstrated between Leishmania species with different tropism, thus contributing to shed light on the mechanisms underlying sexual reproduction in these parasites.

High throughput single-cell genome sequencing gives insights into the generation and evolution of mosaic aneuploidy in Leishmania donovani

Leishmania, a unicellular eukaryotic parasite, is a unique model for aneuploidy and cellular heterogeneity, along with their potential role in adaptation to environmental stresses. Somy variation within clonal populations was previously explored in a small subset of chromosomes using fluorescence hybridization methods. This phenomenon, termed mosaic aneuploidy (MA), might have important evolutionary and functional implications but remains under-explored due to technological limitations. Here, we applied and validated a high throughput single-cell genome sequencing method to study for the first time the extent and dynamics of whole karyotype heterogeneity in two clonal populations of Leishmania promastigotes representing different stages of MA evolution in vitro. We found that drastic changes in karyotypes quickly emerge in a population stemming from an almost euploid founder cell. This possibly involves polyploidization/hybridization at an early stage of population expansion, followed by assorted ploidy reduction. During further stages of expansion, MA increases by moderate and gradual karyotypic alterations, affecting a defined subset of chromosomes. Our data provide the first complete characterization of MA in Leishmania and pave the way for further functional studies.

Stress conditions promote Leishmania hybridization in vitro marked by expression of the ancestral gamete fusogen HAP2 as revealed by single-cell RNA-seq

Leishmania are protozoan parasites transmitted by the bite of sand fly vectors producing a wide spectrum of diseases in their mammalian hosts. These diverse clinical outcomes are directly associated with parasite strain and species diversity. Although Leishmania reproduction is mainly clonal, a cryptic sexual cycle capable of producing hybrid genotypes has been inferred from population genetic studies and directly demonstrated by laboratory crosses. Experimentally, mating competence has been largely confined to promastigotes developing in the sand fly midgut. The ability to hybridize culture promastigotes in vitro has been limited so far to low-efficiency crosses between two Leishmania tropica strains, L747 and MA37, that mate with high efficiency in flies. Here, we show that exposure of promastigote cultures to DNA damage stress produces a remarkably enhanced efficiency of in vitro hybridization of the L. tropica strains and extends to other species, including Leishmania donovani, Leishmania infantum, and Leishmania braziliensis, a capacity to generate intra- and interspecific hybrids. Whole-genome sequencing and total DNA content analyses indicate that the hybrids are in each case full genome, mostly tetraploid hybrids. Single-cell RNA sequencing of the L747 and MA37 parental lines highlights the transcriptome heterogeneity of culture promastigotes and reveals discrete clusters that emerge post-irradiation in which genes potentially involved in genetic exchange are expressed, including the ancestral gamete fusogen HAP2. By generating reporter constructs for HAP2, we could select for promastigotes that could either hybridize or not in vitro. Overall, this work reveals that there are specific populations involved in Leishmania hybridization associated with a discernible transcriptomic signature, and that stress facilitated in vitro hybridization can be a transformative approach to generate large numbers of hybrid genotypes between diverse species and strains.

High genome plasticity and frequent genetic exchange in Leishmania tropica isolates from Afghanistan, Iran and Syria

Background

The kinetoplastid protozoan Leishmania tropica mainly causes cutaneous leishmaniasis in humans in the Middle East, and relapse or treatment failure after treatment are common in this area. L. tropica’s digenic life cycle includes distinct stages in the vector sandfly and the mammalian host. Sexual reproduction and genetic exchange appear to occur more frequently than in other Leishmania species. Understanding these processes is complicated by chromosome instability during cell division that yields aneuploidy, recombination and heterozygosity. This combination of rare recombination and aneuploid permits may reveal signs of hypothetical parasexual mating, where diploid cells fuse to form a transient tetraploid that undergoes chromosomal recombination and gradual chromosomal loss.

Methodology/principal findings

The genome-wide SNP diversity from 22 L. tropica isolates showed chromosome-specific runs of patchy heterozygosity and extensive chromosome copy number variation. All these isolates were collected during 2007–2017 in Sweden from patients infected in the Middle East and included isolates from a patient possessing two genetically distinct leishmaniasis infections three years apart with no evidence of re-infection. We found differing ancestries on the same chromosome (chr36) across multiple samples: matching the reference genome with few derived alleles, followed by blocks of heterozygous SNPs, and then by clusters of homozygous SNPs with specific recombination breakpoints at an inferred origin of replication. Other chromosomes had similar marked changes in heterozygosity at strand-switch regions separating polycistronic transcriptional units.

Conclusion/significance

These large-scale intra- and inter-chromosomal changes in diversity driven by recombination and aneuploidy suggest multiple mechanisms of cell reproduction and diversification in L. tropica, including mitotic, meiotic and parasexual processes. It underpins the need for more genomic surveillance of Leishmania, to detect emerging hybrids that could spread more widely and to better understand the association between genetic variation and treatment outcome. Furthering our understanding of Leishmania genome evolution and ancestry will aid better diagnostics and treatment for cutaneous leishmaniasis caused by L.tropica in the Middle East.

Cutaneous leishmaniasis is mainly caused by Leishmania tropica in the Middle East, where it is known for treatment failure and a need for prolonged and/or multiple treatments. Several factors affect the clinical presentation and treatment outcome, such as host genetic variability and specific immune response, as well as environmental factors and the vector species. Little is known about the parasite genome and its influence on treatment response. By analysing the genome of 22 isolates of L. tropica, we have revealed extensive genomic variation and a complex population structure with evidence of genetic exchange within and among the isolates, indicating a possible presence of sexual or parasexual mechanisms. Understanding the Leishmania genome better may improve future treatment and better understanding of treatment failure and relapse.

Leishmania and the Model of Predominant Clonal Evolution

As it is the case for other pathogenic microorganisms, the respective impact of clonality and genetic exchange on Leishmania natural populations has been the object of lively debates since the early 1980s. The predominant clonal evolution (PCE) model states that genetic exchange in these parasites’ natural populations may have a high relevance on an evolutionary scale, but is not sufficient to erase a persistent phylogenetic signal and the existence of bifurcating trees. Recent data based on high-resolution markers and genomic polymorphisms fully confirm the PCE model down to a microevolutionary level.

Balancing de novo synthesis and salvage of lipids by Leishmania amastigotes

Leishmania parasites replicate as flagellated, extracellular promastigotes in the sand fly vector and then differentiate into non-flagellated, intracellular amastigotes in the vertebrate host. Promastigotes rely on de novo synthesis to produce the majority of their lipids including glycerophospholipids, sterols and sphingolipids. In contrast, amastigotes acquire most of their lipids from the host although they retain some capacity for de novo synthesis. The switch from de novo synthesis to salvage reflects the transition of Leishmania from fast-replicating promastigotes to slow-growing, metabolically quiescent amastigotes. Future studies will reveal the uptake and remodeling of host lipids by amastigotes at the cellular and molecular levels. Blocking the lipid transfer from host to parasites may present a novel strategy to control Leishmania growth.

Serine proteases profiles of Leishmania (Viannia) braziliensis clinical isolates with distinct susceptibilities to antimony

Glucantime (Sb^V) is the first-line treatment against American Tegumentary Leishmaniasis. Resistance cases to this drug have been reported and related to host characteristics and parasite phenotypes. In this study, 12 Leishmania (Viannia) braziliensis isolates from patients that presented clinical cure (Responders—R) and relapse or therapeutic failure (Non-responders—NR) after treatment with antimony, were analyzed. These parasites were assessed by in vitro susceptibility to Sb^III and Sb^V, serine proteases activity measured with substrate (z-FR-AMC) and specific inhibitors (TLCK, AEBSF and PMSF). In vitro susceptibility of axenic amastigotes to Sb^III showed a significant difference between R and NR groups. The protease assays showed that TLCK inhibited almost 100% of activity in both axenic amastigotes and promastigotes while AEBSF inhibited around 70%, and PMSF showed lower inhibition of some isolates. Principal component and clustering analysis performed with these data yielded one homogeneous cluster with only NR isolates and three heterogeneous clusters with R and NR isolates. Additionally, differential expression of subtilisins (LbrM.13.0860 and LbrM.28.2570) and TXNPx (LbrM.15.1080) was evaluated in promastigotes and axenic amastigotes from both groups. The results showed a higher expression of LbrM.13.0860 and LbrM.15.1080 genes in axenic amastigotes, while LbrM.28.2570 gene had the lowest expression in all isolates, regardless of the parasite form. The data presented here show a phenotypic heterogeneity among the parasites, suggesting that exploration of in vitro phenotypes based on Sb^III and serine proteases profiles can aid in the characterization of L. (V.) braziliensis clinical isolates.

Reproduction in Trypanosomatids: Past and Present

Simple Summary

The reproduction of trypanosomatids is a fundamental issue for host–parasite interaction, and its biological importance lies in knowing how these species acquire new defense mechanisms against the countermeasures imposed by the host, which is consistent with the theory of the endless race or the Red Queen hypothesis for the existence of meiotic sex. Moreover, the way these species re-produce may also be at the origin of novel and more virulent clades and is relevant from a thera-peutic or vaccination point of view, as sex may contribute to increased tolerance and even to the rapid acquisition of drug resistance mechanisms. Kinetoplastids are single-celled organisms, many of them being responsible for important parasitic diseases, globally termed neglected diseases, which are endemic in low-income countries. Leishmaniasis, African (sleeping sickness) and American trypanosomiasis (Chagas disease) caused by trypanosomatids are among the most ne-glected tropical scourges related to poverty and poor health systems. The reproduction of these microorganisms has long been considered to be clonal due to population genetic observations. However, there is increasing evidence of true sex and genetic exchange events under laboratory conditions. We would like to highlight the importance of this topic in the field of host/parasite in-terplay, virulence, and drug resistance.

Abstract

Diseases caused by trypanosomatids (Sleeping sickness, Chagas disease, and leishmaniasis) are a serious public health concern in low-income endemic countries. These diseases are produced by single-celled parasites with a diploid genome (although aneuploidy is frequent) organized in pairs of non-condensable chromosomes. To explain the way they reproduce through the analysis of natural populations, the theory of strict clonal propagation of these microorganisms was taken as a rule at the beginning of the studies, since it partially justified their genomic stability. However, numerous experimental works provide evidence of sexual reproduction, thus explaining certain naturally occurring events that link the number of meiosis per mitosis and the frequency of mating. Recent techniques have demonstrated genetic exchange between individuals of the same species under laboratory conditions, as well as the expression of meiosis specific genes. The current debate focuses on the frequency of genomic recombination events and its impact on the natural parasite population structure. This paper reviews the results and techniques used to demonstrate the existence of sex in trypanosomatids, the inheritance of kinetoplast DNA (maxi- and minicircles), the impact of genetic exchange in these parasites, and how it can contribute to the phenotypic diversity of natural populations.

Molecular tools confirm natural Leishmania (Viannia) guyanensis/L. (V.) shawi hybrids causing cutaneous leishmaniasis in the Amazon region of Brazil

Seven isolates from patients with American cutaneous leishmaniasis in the Amazon region of Brazil were phenotypically suggestive of Leishmania (Viannia) guyanensis/L. (V.) shawi hybrids. In this work, two molecular targets were employed to check the hybrid identity of the putative hybrids. Heat shock protein 70 (hsp70) gene sequences were analyzed by three different polymerase chain reaction (PCR) approaches, and two different patterns of inherited hsp70 alleles were found. Three isolates presented heterozygous L. (V.) guyanensis/L. (V.) shawi patterns, and four presented homozygous hsp70 patterns involving only L. (V.) shawi alleles. The amplicon sequences confirmed the RFLP patterns. The high-resolution melting method detected variant heterozygous and homozygous profiles. Single-nucleotide polymorphism genotyping/cleaved amplified polymorphic site analysis suggested a higher contribution from L. (V.) guyanensis in hsp70 heterozygous hybrids. Additionally, PCR-RFLP analysis targeting the enzyme mannose phosphate isomerase (mpi) gene indicated heterozygous and homozygous cleavage patterns for L. (V.) shawi and L. (V.) guyanensis, corroborating the hsp70 findings. In this communication, we present molecular findings based on partial informative regions of the coding sequences of hsp70 and mpi as markers confirming that some of the parasite strains from the Brazilian Amazon region are indeed hybrids between L. (V.) guyanensis and L. (V.) shawi.

In Vitro Generation of Leishmania Hybrids

Protozoan parasites in the genus Leishmania produce a broad spectrum of diseases in their human hosts. The strain and species-specific genes controlling these diverse clinical outcomes have remained poorly tractable using reverse genetics approaches. A cryptic sexual cycle involving a meiotic-like process has been described in Leishmania but is so far confined to parasites growing in the sand fly vector. Here, we describe the reproducible in vitro generation of hybrid clones using axenic culture forms of Leishmania tropica promastigotes. Analysis of SNPs marker inheritance and whole-genome sequencing data indicate that the progeny clones are full genomic hybrids. The demonstration that mating-competent forms arise in culture should facilitate experimental study of the mating biology of Leishmania and the generation of large numbers of recombinant parasites for positional cloning of important genes.

The sexual side of parasitic protists

Much of the vast evolutionary landscape occupied by Eukaryotes is dominated by protists. Though parasitism has arisen in many lineages, there are three main groups of parasitic protists of relevance to human and livestock health: the Apicomplexa, including the malaria parasite Plasmodium and coccidian pathogens of livestock such as Eimeria; the excavate flagellates, encompassing a diverse range of protist pathogens including trypanosomes, Leishmania, Giardia and Trichomonas; and the Amoebozoa, including pathogenic amoebae such as Entamoeba. These three groups represent separate, deep branches of the eukaryote tree, underlining their divergent evolutionary histories. Here, I explore what is known about sex in these three main groups of parasitic protists.

De Novo Synthesis of Phosphatidylcholine Is Essential for the Promastigote But Not Amastigote Stage in Leishmania major

Phosphatidylcholine (PC) is the most abundant type of phospholipids in eukaryotes constituting ~30% of total lipids in Leishmania. PC synthesis mainly occurs via the choline branch of the Kennedy pathway (choline ⇒ choline-phosphate ⇒ CDP-choline ⇒ PC) and the N-methylation of phosphatidylethanolamine (PE). In addition, Leishmania parasites can acquire PC and other lipids from the host or culture medium. In this study, we assessed the function and essentiality of choline ethanolamine phosphotransferase (CEPT) in Leishmania major which is responsible for the final step of the de novo synthesis of PC and PE. Our data indicate that CEPT is localized in the endoplasmic reticulum and possesses the activity to generate PC from CDP-choline and diacylglycerol. Targeted deletion of CEPT is only possible in the presence of an episomal CEPT gene in the promastigote stage of L. major. These chromosomal null parasites require the episomal expression of CEPT to survive in culture, confirming its essentiality during the promastigote stage. In contrast, during in vivo infection of BALB/c mice, these chromosomal null parasites appeared to lose the episomal copy of CEPT while maintaining normal levels of virulence, replication and cellular PC. Therefore, while the de novo synthesis of PC/PE is indispensable for the proliferation of promastigotes, intracellular amastigotes appear to acquire most of their lipids through salvage and remodeling.

Leishmania Sexual Reproductive Strategies as Resolved through Computational Methods Designed for Aneuploid Genomes

A cryptic sexual reproductive cycle in Leishmania has been inferred through population genetic studies revealing the presence of hybrid genotypes in natural isolates, with attempts made to decipher sexual strategies by studying complex chromosomal inheritance patterns. A more informative approach is to study the products of controlled, laboratory-based experiments where known strains or species are crossed in the sand fly vector to generate hybrid progeny. These hybrids can be subsequently studied through high resolution sequencing technologies and software suites such as PAINT that disclose inheritance patterns including ploidies, parental chromosome contributions and recombinations, all of which can inform the sexual strategy. In this work, we discuss the computational methods in PAINT that can be used to interpret the sexual strategies adopted specifically by aneuploid organisms and summarize how PAINT has been applied to the analysis of experimental hybrids to reveal meiosis-like sexual recombination in Leishmania.

Study on the Occurrence of Genetic Exchange Among Parasites of the Leishmania mexicana Complex

In Leishmania, genetic exchange has been experimentally demonstrated to occur in the sand fly vector and in promastigote axenic cultures through a meiotic-like process. No evidence of genetic exchange in mammalian hosts have been reported so far, possibly due to the fact that the Leishmania species used in previous studies replicate within individual parasitophorous vacuoles. In the present work, we explored the possibility that residing in communal vacuoles may provide conditions favorable for genetic exchange for L. mexicana and L. amazonensis. Using promastigote lines of both species harboring integrated or episomal drug-resistance markers, we assessed whether genetic exchange can occur in axenic cultures, in infected macrophages as well as in infected mice. We obtained evidence of genetic exchange for L. amazonensis in both axenic promastigote cultures and infected macrophages. However, the resulting products of those putative genetic events were unstable as they did not sustain growth in subsequent sub-cultures, precluding further characterization.

Designing and developing a high-resolution melting technique for accurate identification of Leishmania species by targeting amino acid permease 3 and cytochrome oxidase II genes using real-time PCR and in silico genetic evaluation

Discrimination, accurate identification, and reliable techniques are required for accurate identification of Leishmania parasites. High-resolution melting (HRM) is recognized as an authentic and exact method. The main objective of this research was optimizing HRM analysis for detecting and screening Leishmania major, Leishmania tropica and mix infections. Thirty-six DNA samples of Leishmania parasite were prepared and analyzed. Two gene regions of amino acid permease 3 (AAP3) and cytochrome oxidase II (COII) were targeted and six pairs of specific new primers were designed. Bioinformatics analysis was employed to predict DNA temperature resolution for each species and compared with in-vitro results. The genetic diversity of the selected gene regions was analyzed using PCR-sequencing method and DnaSP 5.10.01 software. They were submitted in GenBank (KU680818- KU680821 and KY041643- KY041649). The haplotype diversity for both AAP3 and COII genes was 96% and 87%, respectively. Tajima's D index was 0.65 for AAP3 and 0.36 for COII. CLC Genomics Workbench 11 software predictions were significant and close to these findings. The designed primers could be able to identify at least two Leishmania species. Temperature variations in HRM technique separated Iranian Leishmania parasites of L. major, L. tropica and mix infections. The target genes and our modified HRM method proved this technique could be useful in both clinical and experimental settings. Also, it can be effective for detecting Leishmania parasites in different hosts such as humans, reservoir hosts and vectors. Indeed, HRM can be used as a technique in Leishmania identification as well as for ecological and epidemiological research.

Binding of Leishmania infantum Lipophosphoglycan to the Midgut Is Not Sufficient To Define Vector Competence in Lutzomyia longipalpis Sand Flies

It is well established that the presence of LPG is sufficient to define the vector competence of restrictive sand fly vectors with respect to Leishmania parasites. However, the permissiveness of other sand flies with respect to multiple Leishmania species suggests that other factors might define vector competence for these vectors. In this study, we investigated the underpinnings of Leishmania infantum survival and development in its natural vector, Lutzomyia longipalpis. We found that LPG-mediated midgut binding persists in late-stage parasites. This observation is of relevance for the understanding of vector-parasite molecular interactions and suggests that only a subset of infective metacyclic-stage parasites (metacyclics) lose their ability to attach to the midgut, with implications for parasite transmission dynamics. However, our data also demonstrate that LPG is not a determining factor in Leishmania infantum retention in the midgut of Lutzomyia longipalpis, a permissive vector. Rather, LPG appears to be more important in protecting some parasite strains from the toxic environment generated during blood meal digestion in the insect gut. Thus, the relevance of LPG in parasite development in permissive vectors appears to be a complex issue and should be investigated on a strain-specific basis.

The major surface lipophosphoglycan (LPG) of Leishmania parasites is critical to vector competence in restrictive sand fly vectors in mediating Leishmania attachment to the midgut epithelium, considered essential to parasite survival and development. However, the relevance of LPG for sand flies that harbor multiple species of Leishmania remains elusive. We tested binding of Leishmania infantum wild-type (WT), LPG-defective (Δlpg1 mutants), and add-back (Δlpg1 + LPG1) lines to sand fly midguts in vitro and their survival in Lutzomyia longipalpis sand flies in vivo. Le. infantum WT parasites attached to the Lu. longipalpis midgut in vitro, with late-stage parasites binding to midguts in significantly higher numbers than were seen with early-stage promastigotes. Δlpg1 mutants did not bind to Lu. longipalpis midguts, and this was rescued in the Δlpg1 + LPG1 lines, indicating that midgut binding is mediated by LPG. When Lu. longipalpis sand flies were infected with the Le. infantum WT or Le. infantum Δlpg1 or Le. infantum Δlpg1 + LPG1 line of the BH46 or BA262 strains, the BH46 Δlpg1 mutant, but not the BA262 Δlpg1 mutant, survived and grew to numbers similar to those seen with the WT and Δlpg1 + LPG1 lines. Exposure of BH46 and BA262 Δlpg1 mutants to blood-engorged midgut extracts led to mortality of the BA262 Δlpg1 but not the BH46 Δlpg1 parasites. These findings suggest that Le. infantum LPG protects parasites on a strain-specific basis early in infection, likely against toxic components of blood digestion, but that it is not necessary to prevent Le. infantum evacuation along with the feces in the permissive vector Lu. longipalpis.

IMPORTANCE It is well established that the presence of LPG is sufficient to define the vector competence of restrictive sand fly vectors with respect to Leishmania parasites. However, the permissiveness of other sand flies with respect to multiple Leishmania species suggests that other factors might define vector competence for these vectors. In this study, we investigated the underpinnings of Leishmania infantum survival and development in its natural vector, Lutzomyia longipalpis. We found that LPG-mediated midgut binding persists in late-stage parasites. This observation is of relevance for the understanding of vector-parasite molecular interactions and suggests that only a subset of infective metacyclic-stage parasites (metacyclics) lose their ability to attach to the midgut, with implications for parasite transmission dynamics. However, our data also demonstrate that LPG is not a determining factor in Leishmania infantum retention in the midgut of Lutzomyia longipalpis, a permissive vector. Rather, LPG appears to be more important in protecting some parasite strains from the toxic environment generated during blood meal digestion in the insect gut. Thus, the relevance of LPG in parasite development in permissive vectors appears to be a complex issue and should be investigated on a strain-specific basis.

A Way Straight-Forward for Leishmania Genetics

Genetic exchange between Leishmania parasites was demonstrated in sandflies over 10 years ago. Louradour et al. have shown in vitro hybridization of two Leishmania tropica isolates, with the potential to remove a major roadblock to using forward genetics in Leishmania, understanding Leishmania reproductive biology, and analyzing gene flow in natural populations.

Sand Fly Studies Predict Transmission Potential of Drug-resistant Leishmania

Leishmania parasites have the capacity to rapidly adapt to changing environments in their digenetic life cycle which alternates between a vertebrate and an invertebrate host. Emergence of resistance following drug exposure can evoke phenotypic alterations that affect several aspects of parasite fitness in both hosts. Current studies of the impact of resistance are mostly limited to interactions with the mammalian host and characterization of in vitro parasite growth and differentiation. Development in the vector and transmission capacity have been largely ignored. This review reflects on the impact of drug resistance on its spreading potential with specific focus on the use of the sand fly infection model to evaluate parasite development in the vector and the ensuing transmission potential of drug-resistant phenotypes.

Genomic analysis of natural intra-specific hybrids among Ethiopian isolates of Leishmania donovani

Parasites of the genus Leishmania (Kinetoplastida: Trypanosomatidae) cause widespread and devastating human diseases. Visceral leishmaniasis due to Leishmania donovani is endemic in Ethiopia where it has also been responsible for major epidemics. The presence of hybrid genotypes has been widely reported in surveys of natural populations, genetic variation reported in a number of Leishmania species, and the extant capacity for genetic exchange demonstrated in laboratory experiments. However, patterns of recombination and the evolutionary history of admixture that produced these hybrid populations remain unclear. Here, we use whole-genome sequence data to investigate Ethiopian L. donovani isolates previously characterized as hybrids by microsatellite and multi-locus sequencing. To date there is only one previous study on a natural population of Leishmania hybrids based on whole-genome sequences. We propose that these hybrids originate from recombination between two different lineages of Ethiopian L. donovani occurring in the same region. Patterns of inheritance are more complex than previously reported with multiple, apparently independent, origins from similar parents that include backcrossing with parental types. Analysis indicates that hybrids are representative of at least three different histories. Furthermore, isolates were highly polysomic at the level of chromosomes with differences between parasites recovered from a recrudescent infection from a previously treated individual. The results demonstrate that recombination is a significant feature of natural populations and contributes to the growing body of data that shows how recombination, and gene flow, shape natural populations of Leishmania.

Global genome diversity of the Leishmania donovani complex

Protozoan parasites of the Leishmania donovani complex - L. donovani and L. infantum - cause the fatal disease visceral leishmaniasis. We present the first comprehensive genome-wide global study, with 151 cultured field isolates representing most of the geographical distribution. L. donovani isolates separated into five groups that largely coincide with geographical origin but vary greatly in diversity. In contrast, the majority of L. infantum samples fell into one globally-distributed group with little diversity. This picture is complicated by several hybrid lineages. Identified genetic groups vary in heterozygosity and levels of linkage, suggesting different recombination histories. We characterise chromosome-specific patterns of aneuploidy and identified extensive structural variation, including known and suspected drug resistance loci. This study reveals greater genetic diversity than suggested by geographically-focused studies, provides a resource of genomic variation for future work and sets the scene for a new understanding of the evolution and genetics of the Leishmania donovani complex.

Experimental infections and co-infections with Leishmania braziliensis and Leishmania infantum in two sand fly species, Lutzomyia migonei and Lutzomyia longipalpis

Leishmaniases are neglected tropical diseases and Leishmania (Leishmania) infantum and Leishmania (Viannia) braziliensis are the most important causative agents of leishmaniases in the New World. These two parasite species may co-circulate in a given endemic area but their interactions in the vector have not been studied yet. We conducted experimental infections using both single infections and co-infections to compare the development of L. (L.) infantum (OGVL/mCherry) and L. (V.) braziliensis (XB29/GFP) in Lutzomyia longipalpis and Lutzomyia migonei. Parasite labelling by different fluorescein proteins enabled studying interspecific competition and localization of different parasite species during co-infections. Both Leishmania species completed their life cycle, producing infective forms in both sand fly species studied. The same happens in the co infections, demonstrating that the two parasites conclude their development and do not compete with each other. However, infections produced by L. (L.) infantum reached higher rates and grew more vigorously, as compared to L. (V.) braziliensis. In late-stage infections, L. (L.) infantum was present in all midgut regions, showing typical suprapylarian type of development, whereas L. (V.) braziliensis was concentrated in the hindgut and the abdominal midgut (peripylarian development). We concluded that both Lu. migonei and Lu. longipalpis are equally susceptible vectors for L. (L.) infantum, in laboratory colonies. In relation to L. (V.) braziliensis, Lu. migonei appears to be more susceptible to this parasite than Lu. longipalpis.

Global selective sweep of a highly inbred genome of the cattle parasite Neospora caninum

Neospora caninum, a cyst-forming apicomplexan parasite, is a leading cause of neuromuscular diseases in dogs as well as fetal abortion in cattle worldwide. The importance of the domestic and sylvatic life cycles of Neospora, and the role of vertical transmission in the expansion and transmission of infection in cattle, is not sufficiently understood. To elucidate the population genomics of Neospora, we genotyped 50 isolates collected worldwide from a wide range of hosts using 19 linked and unlinked genetic markers. Phylogenetic analysis and genetic distance indices resolved a single genotype of N. caninum Whole-genome sequencing of 7 isolates from 2 different continents identified high linkage disequilibrium, significant structural variation, but only limited polymorphism genome-wide, with only 5,766 biallelic single nucleotide polymorphisms (SNPs) total. Greater than half of these SNPs (∼3,000) clustered into 6 distinct haploblocks and each block possessed limited allelic diversity (with only 4 to 6 haplotypes resolved at each cluster). Importantly, the alleles at each haploblock had independently segregated across the strains sequenced, supporting a unisexual expansion model that is mosaic at 6 genomic blocks. Integrating seroprevalence data from African cattle, our data support a global selective sweep of a highly inbred livestock pathogen that originated within European dairy stock and expanded transcontinentally via unisexual mating and vertical transmission very recently, likely the result of human activities, including recurrent migration, domestication, and breed development of bovid and canid hosts within similar proximities.

Adaptation by Loss of Heterozygosity in Saccharomyces cerevisiae Clones Under Divergent Selection

Loss of heterozygosity (LOH) is observed during vegetative growth and reproduction of diploid genotypes through mitotic crossovers, aneuploidy caused by nondisjunction, and gene conversion. We aimed to test the role that LOH plays during adaptation of two highly heterozygous Saccharomyces cerevisiae genotypes to multiple environments over a short time span in the laboratory. We hypothesized that adaptation would be observed through parallel LOH events across replicate populations. Using genome resequencing of 70 clones, we found that LOH was widespread with 5.2 LOH events per clone after ∼500 generations. The most common mode of LOH was gene conversion (51%) followed by crossing over consistent with either break-induced replication or double Holliday junction resolution. There was no evidence that LOH involved nondisjunction of whole chromosomes. We observed parallel LOH in both an environment-specific and environment-independent manner. LOH largely involved recombining existing variation between the parental genotypes, but also was observed after de novo, presumably beneficial, mutations occurred in the presence of canavanine, a toxic analog of arginine. One highly parallel LOH event involved the ENA salt efflux pump locus on chromosome IV, which showed repeated LOH to the allele from the European parent, an allele originally derived by introgression from S. paradoxus Using CRISPR-engineered LOH we showed that the fitness advantage provided by this single LOH event was 27%. Overall, we found extensive evidence that LOH could be adaptive and is likely to be a greater source of initial variation than de novo mutation for rapid evolution of diploid genotypes.

Meiotic sex in Chagas disease parasite Trypanosoma cruzi

Genetic exchange enables parasites to rapidly transform disease phenotypes and exploit new host populations. Trypanosoma cruzi, the parasitic agent of Chagas disease and a public health concern throughout Latin America, has for decades been presumed to exchange genetic material rarely and without classic meiotic sex. We present compelling evidence from 45 genomes sequenced from southern Ecuador that T. cruzi in fact maintains truly sexual, panmictic groups that can occur alongside others that remain highly clonal after past hybridization events. These groups with divergent reproductive strategies appear genetically isolated despite possible co-occurrence in vectors and hosts. We propose biological explanations for the fine-scale disconnectivity we observe and discuss the epidemiological consequences of flexible reproductive modes. Our study reinvigorates the hunt for the site of genetic exchange in the T. cruzi life cycle, provides tools to define the genetic determinants of parasite virulence, and reforms longstanding theory on clonality in trypanosomatid parasites.

Comparative Proteomic Analysis of Murine Cutaneous Lesions Induced by Leishmania amazonensis or Leishmania major

Cutaneous leishmaniasisis is the most common clinical form of leishmaniasis and one of the most relevant neglected diseases. It is known that the progress of the disease is species specific and the host's immune response plays an important role in its outcome. However, the pathways that lead to parasite clearance or survival remain unknown. In this work, skin tissue from mice experimentally infected with L. amazonensis, one of the causative agents of cutaneous leishmaniasis in the Amazon region, L. major, another causative agent of cutaneous leishmaniasis in Africa, the Middle East, China, and India, or lipopolysaccharides from Escherichia coli as an inflammation model were investigated using label-free proteomics to unveil Leishmania-specific protein alterations. Proteomics is a powerful tool to investigate host-pathogen relationships to address biological questions. In this work, proteins from mice skin biopsies were identified and quantified using nano-LC coupled with tandem mass spectrometry analyses. Integrated Proteomics Pipeline was used for peptide/protein identification and quantification. Western blot was used for validation of protein quantification by mass spectrometry, and protein pathways were predicted using Ingenuity Pathway Analysis. In this proteomics study, several proteins were pointed out as hypothetical targets to guide future studies on Leishmania-specific modulation of proteins in the host. We identified hundreds of exclusively modulated proteins after Leishmania spp. infection and 17 proteins that were differentially modulated in the host after L. amazonensis or L. major infection.

The impact of refugees on leishmaniasis in Turkey: a new Syrian/Turkish Leishmania tropica population structure described by multilocus microsatellite typing (MLMT)

Turkey is one of the leishmaniasis endemic countries, and according to the recent reports, more than 45% of the cases were reported from the Southeastern part of Turkey. The disease is endemic in Syria with annually 25,000 cases, and it is emphasized by WHO that the actual number was estimated to be 2-5-fold higher than the reported numbers. Due to the civil war in Syria, more than seven million people were displaced and migrate to neighboring countries. The population structure of Leishmania tropica was investigated in the present study using clinical samples, which were obtained from Syrian patients residing in Turkey. Previously reported database was used to compare the results obtained in the present study. According to the multilocus microsatellite typing profiles, three populations (Şanlıurfa, Mediterranean, and Syrian/Turkish) were identified. Syrian/Turkish population, which is a new structure and identified for the first time in the present study, was comprised of clinical samples obtained from Syrian patients. The newly described population structure was homogeneous and solid comparing to previously identified population structures in Turkey. Further analyses revealed two sub-populations under the main Syrian/Turkish population structure. The findings of the present study revealed that the epidemiological status of leishmaniasis is more complicated than it is estimated. We believe that the data presented here will provide valuable information on the leishmaniasis epidemiology.

Label-Free Proteomic Analysis Reveals Parasite-Specific Protein Alterations in Macrophages Following Leishmania amazonensis, Leishmania major, or Leishmania infantum Infection

Leishmania is an obligate intracellular parasite known to modulate the host cell to survive and proliferate. However, the complexity of host-parasite interactions remains unclear. Also, the outcome of the disease has been recognized to be species-specific and dependent on the host's immune responses. Proteomics has emerged as a powerful tool to investigate the host-pathogen interface, allowing us to deepen our knowledge about infectious diseases. Quantification of the relative amount of proteins in a sample can be achieved using label-free proteomics, and for the first time, we have used it to quantify Leishmania-specific protein alterations in macrophages. Protein extracts were obtained and digested, and peptides were identified and quantified using nano-LC coupled with tandem mass spectrometry analyses. Protein expression was validated by Western blot analysis. Integrated Proteomics Pipeline was used for peptide/protein identification and for quantification and data processing. Ingenuity Pathway Analysis was used for network analysis. In this work, we investigated how this intracellular parasite modulates protein expression on a host macrophage by comparing three different Leishmania species- L. amazonensis, one of the causative agents of cutaneous disease in the Amazon region; L. major, another causative agent of cutaneous leishmaniasis in Africa, the Middle East, China, and India; L. infantum, the causative agent of visceral leishmaniasis affecting humans and dogs in Latin America-and lipopolysaccharide stimulated macrophages as an in vitro inflammation model. Our results revealed that Leishmania infection downregulates apoptosis pathways while upregulating the activation of phagocytes/leukocytes and lipid accumulation.

Identification of divergent Leishmania (Viannia) braziliensis ecotypes derived from a geographically restricted area through whole genome analysis

Leishmania braziliensis, the main etiological agent of cutaneous leishmaniasis (CL) in Latin America, is characterized by major differences in basic biology in comparison with better-known Leishmania species. It is also associated with a high phenotypic and possibly genetic diversity that need to be more adequately defined. Here we used whole genome sequences to evaluate the genetic diversity of ten L. braziliensis isolates from a CL endemic area from Northeastern Brazil, previously classified by Multi Locus Enzyme Electrophoresis (MLEE) into ten distinct zymodemes. These sequences were first mapped using the L. braziliensis M2904 reference genome followed by identification of Single Nucleotide Polymorphisms (SNPs). A substantial level of diversity was observed when compared with the reference genome, with SNP counts ranging from ~95,000 to ~131,000 for the different isolates. When the genome data was used to infer relationship between isolates, those belonging to zymodemes Z72/Z75, recovered from forested environments, were found to cluster separately from the others, generally associated with more urban environments. Among the remaining isolates, those from zymodemes Z74/Z106 were also found to form a separate group. Phylogenetic analyses were also performed using Multi-Locus Sequence Analysis from genes coding for four metabolic enzymes used for MLEE as well as the gene sequence coding for the Hsp70 heat shock protein. All 10 isolates were firmly identified as L. braziliensis, including the zymodeme Z26 isolate previously classified as Leishmania shawi, with the clustering into three groups confirmed. Aneuploidy was also investigated but found in general restricted to chromosome 31, with a single isolate, from zymodeme Z27, characterized by extra copies for other chromosomes. Noteworthy, both Z72 and Z75 isolates are characterized by a much reduced heterozygosity. Our data is consistent with the existence of distinct evolutionary groups in the restricted area sampled and a substantial genetic diversity within L. braziliensis.

Whole genome sequencing of experimental hybrids supports meiosis-like sexual recombination in Leishmania

Hybrid genotypes have been repeatedly described among natural isolates of Leishmania, and the recovery of experimental hybrids from sand flies co-infected with different strains or species of Leishmania has formally demonstrated that members of the genus possess the machinery for genetic exchange. As neither gamete stages nor cell fusion events have been directly observed during parasite development in the vector, we have relied on a classical genetic analysis to determine if Leishmania has a true sexual cycle. Here, we used whole genome sequencing to follow the chromosomal inheritance patterns of experimental hybrids generated within and between different strains of L. major and L. infantum. We also generated and sequenced the first experimental hybrids in L. tropica. We found that in each case the parental somy and allele contributions matched the inheritance patterns expected under meiosis 97–99% of the time. The hybrids were equivalent to F1 progeny, heterozygous throughout most of the genome for the markers that were homozygous and different between the parents. Rare, non-Mendelian patterns of chromosomal inheritance were observed, including a gain or loss of somy, and loss of heterozygosity, that likely arose during meiosis or during mitotic divisions of the progeny clones in the fly or culture. While the interspecies hybrids appeared to be sterile, the intraspecies hybrids were able to produce backcross and outcross progeny. Analysis of 5 backcross and outcross progeny clones generated from an L. major F1 hybrid, as well as 17 progeny clones generated from backcrosses involving a natural hybrid of L. tropica, revealed genome wide patterns of recombination, demonstrating that classical crossing over occurs at meiosis, and allowed us to construct the first physical and genetic maps in Leishmania. Altogether, the findings provide strong evidence for meiosis-like sexual recombination in Leishmania, presenting clear opportunities for forward genetic analysis and positional cloning of important genes.

Leishmania promastigotes are able to undergo genetic exchange during their growth and development in the sand fly vector, however, it is still not known if they have a true sexual cycle involving meiosis. Here, we used whole genome sequencing to follow the chromosomal inheritance patterns of 44 experimental hybrids generated between different strains of L. major, L. infantum, and L. tropica. In almost every case the number of chromosomes and the allele contributions from each parent matched the inheritance patterns expected under meiosis. Rare instances of hybrid chromosomes that did not fit Mendelian expectations were observed, including gain or loss of somy, and loss of heterozygosity. Strong evidence for a meiotic-like process was also obtained from the genome wide patterns of recombination observed in the offspring generated from backcrosses involving an experimental or natural hybrid, consistent with crossing over occurring between homologous chromosomes during meiosis. The frequency and position of the recombination breakpoints observed on each chromosome allowed us to construct the first physical and genetic maps in Leishmania. The results demonstrate that forward genetic approaches are possible for positional cloning of important genes.

Elucidating in vitro and in vivo phenotypic behaviour of L. infantum/L. major natural hybrids

The clinical manifestation and course of Leishmania infections depend on factors such as species, virulence and host-immunity. Although trypanosomatids are considered to have clonal propagation, genetic hybridization has produced successful natural hybrid lineages. Hybrids displaying strong selective advantages may have an impact on pathogenesis and the eco-epidemiology of leishmaniasis. Thus, characterization of phenotypic properties of Leishmania hybrids could bring significant insight into the biology, infectivity, pathogenicity and transmission dynamics of these atypical strains. The present study focuses on phenotypic features and survival capacity of Leishmania infantum/Leishmania major hybrid isolates as compared with representative putative parental species, L. infantum and L. major. In vitro assays (growth kinetics, susceptibility to different conditions) and in vivo infection (parasite detection and histopathological alterations) showed that hybrids present higher growth capacity and decreased susceptibility to reactive oxygen species. Furthermore, evaluation of infected spleen tissue suggests that hybrids induce a stronger immune reaction than their putative parents, leading to the development of white pulp hyperplasia in B-lymphocyte compartments. Overall, these hybrids have shown high plasticity in terms of their general behaviour within the different phenotypic parameters, suggesting that they might have acquired genetic features conferring different mechanisms to evade host cells.

The recombinase Rad51 plays a key role in events of genetic exchange in Trypanosoma cruzi

Detection of genetic exchange has been a limiting factor to deepen the knowledge on the mechanisms by which Trypanosoma cruzi is able to generate progeny and genetic diversity. Here we show that incorporation of halogenated thymidine analogues, followed by immunostaining, is a reliable method not only to detect T. cruzi fused-cell hybrids, but also to quantify their percentage in populations of this parasite. Through this approach, we were able to detect and quantify fused-cell hybrids of T. cruzi clones CL Brener and Y. Given the increased detection of fused-cell hybrids in naturally-occurring hybrid CL Brener strain, which displays increased levels of RAD51 and BRCA2 transcripts, we further investigated the role of Rad51 - a recombinase involved in homologous recombination - in the process of genetic exchange. We also verified that the detection of fused-cell hybrids in T. cruzi overexpressing RAD51 is increased when compared to wild-type cells, suggesting a key role for Rad51 either in the formation or in the stabilization of fused-cell hybrids in this organism.

Recent advances in trypanosomatid research: genome organization, expression, metabolism, taxonomy and evolution

Unicellular flagellates of the family Trypanosomatidae are obligatory parasites of invertebrates, vertebrates and plants. Dixenous species are aetiological agents of a number of diseases in humans, domestic animals and plants. Their monoxenous relatives are restricted to insects. Because of the high biological diversity, adaptability to dramatically different environmental conditions, and omnipresence, these protists have major impact on all biotic communities that still needs to be fully elucidated. In addition, as these organisms represent a highly divergent evolutionary lineage, they are strikingly different from the common 'model system' eukaryotes, such as some mammals, plants or fungi. A number of excellent reviews, published over the past decade, were dedicated to specialized topics from the areas of trypanosomatid molecular and cell biology, biochemistry, host-parasite relationships or other aspects of these fascinating organisms. However, there is a need for a more comprehensive review that summarizing recent advances in the studies of trypanosomatids in the last 30 years, a task, which we tried to accomplish with the current paper.

Gene and Chromosomal Copy Number Variations as an Adaptive Mechanism Towards a Parasitic Lifestyle in Trypanosomatids

Trypanosomatids are a group of kinetoplastid parasites including some of great public health importance, causing debilitating and life-long lasting diseases that affect more than 24 million people worldwide. Among the trypanosomatids, Trypanosoma cruzi, Trypanosoma brucei and species from the Leishmania genus are the most well studied parasites, due to their high prevalence in human infections. These parasites have an extreme genomic and phenotypic variability, with a massive expansion in the copy number of species-specific multigene families enrolled in host-parasite interactions that mediate cellular invasion and immune evasion processes. As most trypanosomatids are heteroxenous, and therefore their lifecycles involve the transition between different hosts, these parasites have developed several strategies to ensure a rapid adaptation to changing environments. Among these strategies, a rapid shift in the repertoire of expressed genes, genetic variability and genome plasticity are key mechanisms. Trypanosomatid genomes are organized into large directional gene clusters that are transcribed polycistronically, where genes derived from the same polycistron may have very distinct mRNA levels. This particular mode of transcription implies that the control of gene expression operates mainly at post-transcriptional level. In this sense, gene duplications/losses were already associated with changes in mRNA levels in these parasites. Gene duplications also allow the generation of sequence variability, as the newly formed copy can diverge without loss of function of the original copy. Recently, aneuploidies have been shown to occur in several Leishmania species and T. cruzi strains. Although aneuploidies are usually associated with debilitating phenotypes in superior eukaryotes, recent data shows that it could also provide increased fitness in stress conditions and generate drug resistance in unicellular eukaryotes. In this review, we will focus on gene and chromosomal copy number variations and their relevance to the evolution of trypanosomatid parasites.

Mannose receptor high, M2 dermal macrophages mediate nonhealing Leishmania major infection in a Th1 immune environment

The origin and functional specialization of dermal macrophages in cutaneous infections have been little studied. In this paper, we show that a strain of Leishmania major (L. major Seidman [LmSd]) that produces nonhealing cutaneous lesions in conventionally resistant C57BL/6 mice was more efficiently taken up by M2-polarized bone marrow (BM)-derived macrophages (BMDMs) in vitro and by mannose receptor (MR)^hi dermal macrophages in vivo compared with a healing strain (L. major Friedlin V1). Both in steady and in T helper type 1 (Th1) cell-driven inflammatory states, the MR^hi dermal macrophages showed M2 characteristics. The dermal macrophages were radio resistant and not replaced by monocytes or adult BM-derived cells during infection, but were locally maintained by IL-4 and IL-10. Notably, the favored infection of M2 BMDMs by LmSd in vitro was MR dependent, and genetic deletion of MR or selective depletion of MR^hi dermal macrophages by anti-CSF-1 receptor antibody reversed the nonhealing phenotype. We conclude that embryonic-derived, MR^hi dermal macrophages are permissive for parasite growth even in a strong Th1-immune environment, and the preferential infection of these cells plays a crucial role in the severity of cutaneous disease.

Gene Expression in Leishmania Is Regulated Predominantly by Gene Dosage

Leishmania tropica, a unicellular eukaryotic parasite present in North and East Africa, the Middle East, and the Indian subcontinent, has been linked to large outbreaks of cutaneous leishmaniasis in displaced populations in Iraq, Jordan, and Syria. Here, we report the genome sequence of this pathogen and 7,863 identified protein-coding genes, and we show that the majority of clinical isolates possess high levels of allelic diversity, genetic admixture, heterozygosity, and extensive aneuploidy. By utilizing paired genome-wide high-throughput DNA sequencing (DNA-seq) with RNA-seq, we found that gene dosage, at the level of individual genes or chromosomal "somy" (a general term covering disomy, trisomy, tetrasomy, etc.), accounted for greater than 85% of total gene expression variation in genes with a 2-fold or greater change in expression. High gene copy number variation (CNV) among membrane-bound transporters, a class of proteins previously implicated in drug resistance, was found for the most highly differentially expressed genes. Our results suggest that gene dosage is an adaptive trait that confers phenotypic plasticity among natural Leishmania populations by rapid down- or upregulation of transporter proteins to limit the effects of environmental stresses, such as drug selection.IMPORTANCELeishmania is a genus of unicellular eukaryotic parasites that is responsible for a spectrum of human diseases that range from cutaneous leishmaniasis (CL) and mucocutaneous leishmaniasis (MCL) to life-threatening visceral leishmaniasis (VL). Developmental and strain-specific gene expression is largely thought to be due to mRNA message stability or posttranscriptional regulatory networks for this species, whose genome is organized into polycistronic gene clusters in the absence of promoter-mediated regulation of transcription initiation of nuclear genes. Genetic hybridization has been demonstrated to yield dramatic structural genomic variation, but whether such changes in gene dosage impact gene expression has not been formally investigated. Here we show that the predominant mechanism determining transcript abundance differences (>85%) in Leishmania tropica is that of gene dosage at the level of individual genes or chromosomal somy.