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

Understanding the genetic complexity of Leishmania infantum in the Americas: a focus on 3'NT/NU gene deletion

Visceral Leishmaniasis in the Americas is primarily associated with Leishmania (Leishmania) infantum. This parasite is non-native and was imported during the colonisation era. The constitutive instability of the Leishmania genome allows this parasite to express flexibility in adapting to environmental fluctuations and different selective pressures, such as those the parasite faced when arrived in the New World. Therefore, genetic diversity is expected among the populations of L. infantum in the Americas, despite the bottle neck of importation route. Indeed, subpopulation of strains of L. infantum carrying a homozygous deletion in the genome was detected exclusively in the continent. These strains are more spread across Brazilian territory to the detriment of the non-deleted; the locus includes four genes, two of which encode the enzyme ecto-3'-nucleotidase/nuclease (3'NT/NU), a virulence factor in L. infantum. In this review, we highlight the sub estimated genetic complexity of L. infantum populations in Brazil, addressing the biological importance of the 3'NT/NU enzyme and the possible phenotypic impacts of its deletion, pointing out how it may configure an adaptive strategy for L. infantum. Finally, we raise the discussion of how the genome of L. infantum might be shaped in a unique way under the ecological conditions of Brazil.

Chromosome-level genome assembly of Eimeria tenella at the single-oocyst level

BACKGROUND

Eimeria are obligate protozoan parasites, and more than 1,500 species have been reported. However, Eimeria genomes lag behind many other eukaryotes since obtaining many oocysts is difficult due to a lack of sustainable in vitro culture, highly repetitive sequences, and mixed species infections. To address this challenge, we used whole-genome amplification of a single oocyst followed by long-read sequencing and obtained a chromosome-level genome of Eimeria tenella.

RESULTS

The assembled genome was 52.13 Mb long, encompassing 15 chromosomes and 46.94% repeat sequences. In total, 7,296 protein-coding genes were predicted, exhibiting high completeness, with 92.00% single-copy BUSCO genes. To the best of our knowledge, this is the first chromosome-level assembly of E. tenella using a combination of single-oocyst whole-genome amplification and long-read sequencing. Comparative genomic and transcriptome analyses confirmed evolutionary relationship and supported estimates of divergence time of apicomplexan parasites and identified AP2 and Myb gene families that may play indispensable roles in regulating the growth and development of E. tenella.

CONCLUSION

This high-quality genome assembly and the established sequencing strategy provide valuable community resources for comparative genomic and evolutionary analyses of the Eimeria clade. Additionally, our study also provides a valuable resource for exploring the roles of AP2 and Myb transcription factor genes in regulating the development of Eimeria parasites.

The SET29 and SET7 proteins of Leishmania donovani exercise non-redundant convergent as well as collaborative functions in moderating the parasite's response to oxidative stress

SET proteins are lysine methyltransferases. In investigating Leishmania donovani SET29, we found depletion of LdSET29 by two-thirds did not affect promastigote growth, nor alter the parasite's response to UV-induced or HU-induced stress, but made it more tolerant to H2O2-induced oxidizing environment. The deviant response to oxidative stress was coupled to lowered accumulation of reactive oxygen species, which was linked to enhanced scavenging activity. The set29 mutants' response to H2O2 exposure was similar to that of set7 mutants, prompting an investigation into genetic and physical interactions between the two proteins. While neither protein could rescue the aberrant phenotypes of the other set mutant, the two proteins interacted physically in vitro and in vivo. Transcriptome analyses revealed that neither protein regulated global gene expression, but LdSET7 controlled transcript levels of a limited number of genes, including several peroxidases. In working towards identifying targets through which SET7/SET29 mediate the cell's response to an oxidative milieu, we found HSP60/CNP60 and TCP1 to be possible candidates. LdHSP60 has earlier been implicated in the regulation of the response of virulent promastigotes to H2O2 exposure, and LdTCP1 has previously been found to have a protective effect against oxidative stress. set7 and set29 mutants survived more proficiently in host macrophages as well. The data suggest an alliance between LdSET29 and LdSET7 in mounting the parasite's response to oxidative stress, each protein playing its own distinctive role. They ensure the parasite not only establishes infection but also maintains the balance with host cells to enable the persistence of infection.

Nanopore sequencing reveals that DNA replication compartmentalisation dictates genome stability and instability in Trypanosoma brucei

The Trypanosoma brucei genome is structurally complex. Eleven megabase-sized chromosomes each comprise a transcribed core flanked by silent subtelomeres, housing thousands of Variant Surface Glycoprotein (VSG) genes. Additionally, hundreds of sub-megabase chromosomes contain 177 bp repeats of unknown function, and VSG transcription sites localise to many telomeres. DNA replication dynamics have only been described in the megabase chromosome cores, and in the single active VSG transcription site. Using a Nanopore genome assembly, we show that megabase chromosome subtelomeres display a paucity of replication initiation events relative to the core, correlating with increased instability. In addition, replication of the active VSG transcription site is shown to originate from the telomere, likely causing targeted VSG recombination. Lastly, we provide evidence that the 177 bp repeats act as conserved DNA replication origins, explaining submegabase chromosome stability. Compartmentalized DNA replication therefore explains how T. brucei balances stable genome transmission with localised instability driving immune evasion.

Evolutionary genomics of Leishmania braziliensis across the neotropical realm

The Neotropical realm, one of the most biodiverse regions on Earth, houses a broad range of zoonoses that pose serious public health threats. Protozoan parasites of the Leishmania (Viannia) braziliensis clade cause zoonotic leishmaniasis in Latin America with clinical symptoms ranging from simple cutaneous to destructive, disfiguring mucosal lesions. We present the first comprehensive genome-wide continental study including 257 cultivated isolates representing most of the geographical distribution of this major human pathogen. The L. braziliensis clade is genetically highly heterogeneous, consisting of divergent parasite groups that are associated with different environments and vary greatly in diversity. Apart from several small ecologically isolated groups with little diversity, our sampling identifies two major parasite groups, one associated with the Amazon and the other with the Atlantic Forest biomes. These groups show different recombination histories, as suggested by high levels of heterozygosity and effective population sizes in the Amazonian group in contrast to high levels of linkage and clonality in the Atlantic group. We argue that these differences are linked to strong eco-epidemiological differences between the two regions. In contrast to geographically focused studies, our study provides a broad understanding of the molecular epidemiology of zoonotic parasites circulating in tropical America.

Diverse Genome Structures among Eukaryotes May Have Arisen in Response to Genetic Conflict

In contrast to the typified view of genome cycling only between haploidy and diploidy, there is evidence from across the tree of life of genome dynamics that alter both copy number (i.e. ploidy) and chromosome complements. Here, we highlight examples of such processes, including endoreplication, aneuploidy, inheritance of extrachromosomal DNA, and chromatin extrusion. Synthesizing data on eukaryotic genome dynamics in diverse extant lineages suggests the possibility that such processes were present before the last eukaryotic common ancestor. While present in some prokaryotes, these features appear exaggerated in eukaryotes where they are regulated by eukaryote-specific innovations including the nucleus, complex cytoskeleton, and synaptonemal complex. Based on these observations, we propose a model by which genome conflict drove the transformation of genomes during eukaryogenesis: from the origin of eukaryotes (i.e. first eukaryotic common ancestor) through the evolution of last eukaryotic common ancestor.

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.

Comprehensive analysis of flavohemoprotein copy number variation in Giardia intestinalis: exploring links to metronidazole resistance

BACKGROUND

Giardiasis, caused by the protozoan parasite Giardia intestinalis, often presents a treatment challenge, particularly in terms of resistance to metronidazole. Despite extensive research, markers for metronidazole resistance have not yet been identified.

METHODS

This study analysed 28 clinical samples of G. intestinalis from sub-assemblage AII, characterised by varying responses to metronidazole treatment. We focussed on copy number variation (CNV) of the multi-copy flavohemoprotein gene, analysed using digital polymerase chain reaction (dPCR) and next generation sequencing (NGS). Additionally, chromosomal ploidy was tested in 18 of these samples. Flavohemoprotein CNV was also assessed in 17 samples from other sub-assemblages.

RESULTS

Analyses revealed variable CNVs of the flavohemoprotein gene among the isolates, with no correlation to clinical metronidazole resistance. Discrepancies in CNVs detected from NGS data were attributed to biases linked to the whole genome amplification. However, dPCR helped to clarify these discrepancies by providing more consistent CNV data. Significant differences in flavohemoprotein CNVs were observed across different G. intestinalis sub-assemblages. Notably, Giardia exhibits a propensity for aneuploidy, contributing to genomic variability within and between sub-assemblages.

CONCLUSIONS

The complexity of the clinical metronidazole resistance in Giardia is influenced by multiple genetic factors, including CNVs and aneuploidy. No significant differences in the CNV of the flavohemoprotein gene between isolates from metronidazole-resistant and metronidazole-sensitive cases of giardiasis were found, underscoring the need for further research to identify reliable genetic markers for resistance. We demonstrate that dPCR and NGS are robust methods for analysing CNVs and provide cross-validating results, highlighting their utility in the genetic analyses of this parasite.

Ancestral aneuploidy and stable chromosomal duplication resulting in differential genome structure and gene expression control in trypanosomatid parasites

Aneuploidy is widely observed in both unicellular and multicellular eukaryotes, usually associated with adaptation to stress conditions. Chromosomal duplication stability is a tradeoff between the fitness cost of having unbalanced gene copies and the potential fitness gained from increased dosage of specific advantageous genes. Trypanosomatids, a family of protozoans that include species that cause neglected tropical diseases, are a relevant group to study aneuploidies. Their life cycle has several stressors that could select for different patterns of chromosomal duplications and/or losses, and their nearly universal use of polycistronic transcription increases their reliance on gene expansion/contraction, as well as post-transcriptional control as mechanisms for gene expression regulation. By evaluating the data from 866 isolates covering seven trypanosomatid genera, we have revealed that aneuploidy tolerance is an ancestral characteristic of trypanosomatids but has a reduced occurrence in a specific monophyletic clade that has undergone large genomic reorganization and chromosomal fusions. We have also identified an ancient chromosomal duplication that was maintained across these parasite's speciation, named collectively as the trypanosomatid ancestral supernumerary chromosome (TASC). TASC has most genes in the same coding strand, is expressed as a disomic chromosome (even having four copies), and has increased potential for functional variation, but it purges highly deleterious mutations more efficiently than other chromosomes. The evidence of stringent control over gene expression in this chromosome suggests that these parasites have adapted to mitigate the fitness cost associated with this ancient chromosomal duplication.

Comparing the Yield of Recombinant Human Factor VII Protein Expressed by the rDNA-Promoter with the CMV-Promoter in Iranian Lizard Leishmania

Background

Iranian Lizard Leishmania (I.L.L) is a nonpathogenic Leishmania strain. Due to its advantages, several recombinant proteins have been produced in this host. However, I.L.L shows a lower yield of recombinant protein expression compared to other commercial hosts. Considering the role of protease enzymes in protein digestion, we selected cysteine protease B (CPB) to investigate its impact on recombinant protein yield in I.L.L.

Methods

we generated gene knockouts by utilizing homologous recombination (HR) and CRISPR methods. To assess the efficacy of the designed construct, we compared the yield of recombinant human factor VII (rhFVII) production between cells transfected with the pLEXSY-hyg2-FVII vector and the CMV-promoter-based construct (pF7cmvneo).

Results

The knockout of a single CPB gene allele through the HR method or the complete knockout of all alleles through the CRISPR method led to cell death. This outcome suggests that even the deletion of a single CPB gene allele diminishes the protein to a level insufficient for the survival of I.L.L, indicating a critical dependency on the presence of this protein for the organism's viability. rhFVII exhibited a greater expression yield with the pLEXSY construct compared to the pF7cmvneo construct in I.L.L. The lower expression rate of pF7cmvneo may be influenced by epigenetic factors related to the CPC gene or the RNA polymerase used for the expression of that promoter.

Conclusion

Therefore, considering alternative integration targets for CMV-promoter-based constructs and incorporating UTR sequences of I.L.L high-expression proteins in the vector may enhance recombinant protein expression rates.

Drug Discovery for Cutaneous Leishmaniasis: A Review of Developments in the Past 15 Years

Leishmaniasis is a group of vector-borne, parasitic diseases caused by over 20 species of the protozoan Leishmania spp. The three major disease classifications, cutaneous, visceral, and mucocutaneous, have a range of clinical manifestations from self-healing skin lesions to hepatosplenomegaly and mucosal membrane damage to fatality. As a neglected tropical disease, leishmaniasis represents a major international health challenge, with nearly 350 million people living at risk of infection a year. The current chemotherapeutics used to treat leishmaniasis have harsh side effects, prolonged and costly treatment regimens, as well as emerging drug resistance, and are predominantly used for the treatment of visceral leishmaniasis. There is an undeniable need for the identification and development of novel chemotherapeutics targeting cutaneous leishmaniasis (CL), largely ignored by concerted drug development efforts. CL is mostly non-lethal and the most common presentation of this disease, with nearly 1 million new cases reported annually. Recognizing this unaddressed need, substantial yet fragmented progress in early drug discovery efforts for CL has occurred in the past 15 years and was outlined in this review. However, further work needs to be carried out to advance early discovery candidates towards the clinic. Importantly, there is a paucity of investment in the translation and development of therapies for CL, limiting the emergence of viable solutions to deal with this serious and complex international health problem.

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.

The adaptive roles of aneuploidy and polyclonality in Leishmania in response to environmental stress

Aneuploidy is generally considered harmful, but in some microorganisms, it can act as an adaptive mechanism against environmental stress. Here, we use Leishmania-a protozoan parasite with remarkable genome plasticity-to study the early steps of aneuploidy evolution under high drug pressure (using antimony or miltefosine as stressors). By combining single-cell genomics, lineage tracing with cellular barcodes, and longitudinal genome characterization, we reveal that aneuploidy changes under antimony pressure result from polyclonal selection of pre-existing karyotypes, complemented by further and rapid de novo alterations in chromosome copy number along evolution. In the case of miltefosine, early parasite adaptation is associated with independent point mutations in a miltefosine transporter gene, while aneuploidy changes only emerge later, upon exposure to increased drug levels. Therefore, polyclonality and genome plasticity are hallmarks of parasite adaptation, but the scenario of aneuploidy dynamics depends on the nature and strength of the environmental stress as well as on the existence of other pre-adaptive mechanisms.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae

BACKGROUND

Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group.

RESULTS

We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy.

CONCLUSIONS

The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

Drug resistance in Leishmania: does it really matter?

Treatment failure (TF) jeopardizes the management of parasitic diseases, including leishmaniasis. From the parasite's point of view, drug resistance (DR) is generally considered as central to TF. However, the link between TF and DR, as measured by in vitro drug susceptibility assays, is unclear, some studies revealing an association between treatment outcome and drug susceptibility, others not. Here we address three fundamental questions aiming to shed light on these ambiguities. First, are the right assays being used to measure DR? Second, are the parasites studied, which are generally those that adapt to in vitro culture, actually appropriate? Finally, are other parasite factors - such as the development of quiescent forms that are recalcitrant to drugs - responsible for TF without DR?

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Life in plastic, it's fantastic! How Leishmania exploit genome instability to shape gene expression

Leishmania are kinetoplastid pathogens that cause leishmaniasis, a debilitating and potentially life-threatening infection if untreated. Unusually, Leishmania regulate their gene expression largely post-transcriptionally due to the arrangement of their coding genes into polycistronic transcription units that may contain 100s of functionally unrelated genes. Yet, Leishmania are capable of rapid and responsive changes in gene expression to challenging environments, often instead correlating with dynamic changes in their genome composition, ranging from chromosome and gene copy number variations to the generation of extrachromosomal DNA and the accumulation of point mutations. Typically, such events indicate genome instability in other eukaryotes, coinciding with genetic abnormalities, but for Leishmania, exploiting these products of genome instability can provide selectable substrates to catalyse necessary gene expression changes by modifying gene copy number. Unorthodox DNA replication, DNA repair, replication stress factors and DNA repeats are recognised in Leishmania as contributors to this intrinsic instability, but how Leishmania regulate genome plasticity to enhance fitness whilst limiting toxic under- or over-expression of co-amplified and co-transcribed genes is unclear. Herein, we focus on fresh, and detailed insights that improve our understanding of genome plasticity in Leishmania. Furthermore, we discuss emerging models and factors that potentially circumvent regulatory issues arising from polycistronic transcription. Lastly, we highlight key gaps in our understanding of Leishmania genome plasticity and discuss future studies to define, in higher resolution, these complex regulatory interactions.

Genome diversity of Leishmania aethiopica

Leishmania aethiopica is a zoonotic Old World parasite transmitted by Phlebotomine sand flies and causing cutaneous leishmaniasis in Ethiopia and Kenya. Despite a range of clinical manifestations and a high prevalence of treatment failure, L. aethiopica is one of the most neglected species of the Leishmania genus in terms of scientific attention. Here, we explored the genome diversity of L. aethiopica by analyzing the genomes of twenty isolates from Ethiopia. Phylogenomic analyses identified two strains as interspecific hybrids involving L. aethiopica as one parent and L. donovani and L. tropica respectively as the other parent. High levels of genome-wide heterozygosity suggest that these two hybrids are equivalent to F1 progeny that propagated mitotically since the initial hybridization event. Analyses of allelic read depths further revealed that the L. aethiopica - L. tropica hybrid was diploid and the L. aethiopica - L. donovani hybrid was triploid, as has been described for other interspecific Leishmania hybrids. When focusing on L. aethiopica, we show that this species is genetically highly diverse and consists of both asexually evolving strains and groups of recombining parasites. A remarkable observation is that some L. aethiopica strains showed an extensive loss of heterozygosity across large regions of the nuclear genome, which likely arose from gene conversion/mitotic recombination. Hence, our prospection of L. aethiopica genomics revealed new insights into the genomic consequences of both meiotic and mitotic recombination in Leishmania.

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.

Paving the Way: Contributions of Big Data to Apicomplexan and Kinetoplastid Research

In the age of big data an important question is how to ensure we make the most out of the resources we generate. In this review, we discuss the major methods used in Apicomplexan and Kinetoplastid research to produce big datasets and advance our understanding of Plasmodium, Toxoplasma, Cryptosporidium, Trypanosoma and Leishmania biology. We debate the benefits and limitations of the current technologies, and propose future advancements that may be key to improving our use of these techniques. Finally, we consider the difficulties the field faces when trying to make the most of the abundance of data that has already been, and will continue to be, generated.

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.