Protozoan parasites of the Leishmania genus produce a spectrum of diseases in their mammalian hosts, including humans, ranging from self-limiting cutaneous lesions and tissue-destructive mucosal involvement to infection of the deep viscera that is fatal in the absence of treatment. These clinical outcomes are associated with an extraordinary diversity of Leishmania species and strains, with more than 20 species that are pathogenic to humans. All members of the genus have a dimorphic life cycle consisting of extracellular promastigotes that multiply asexually in the digestive tract of their sand fly vectors, and amastigotes that multiply asexually in the phagosomes of their mammalian host cells, primarily macrophages. Leishmania were long considered to be essentially clonal, with genetic diversity arising from gradual accumulation of somatic mutations (Tibayrenc et al., 1993).

The reproductive strategies of Leishmania are now known to include a cryptic sexual cycle that has been inferred from the analysis of hybrid genotypes observed in natural isolates (Bañuls et al., 1997; Dujardin et al., 1995; Nolder et al., 2007; Odiwuor et al., 2011; Ravel et al., 2006; Rogers et al., 2014) and directly demonstrated by the generation of hybrids between different Leishmania strains and species in the laboratory (Akopyants et al., 2009; Inbar et al., 2013; Romano et al., 2014; Sadlova et al., 2011). The latter involved mating between extracellular promastigotes stages developing in the sand fly and was achieved by coinfection of flies with two parental lines bearing different drug resistance markers, with subsequent selection of double-drug-resistant cells. Experimentally at least, the sexual cycle is nonobligatory, relatively rare, and confined to life-cycle stages present in the sand fly midgut. Based on whole-genome sequencing analyses, the allele inheritance patterns of experimental hybrids provide strong evidence that the system of genetic exchange in Leishmania is Mendelian and involves meiosis-like sexual recombination (Inbar et al., 2019). Homologues of meiosis-specific genes are found within the genomes of Leishmania and are expressed by stages in the sand fly, including the core meiotic genes SPO11, HOP1, and DMC1 involved in creating DNA double-strand breaks, homologous chromosome alignment, and recombination (Inbar et al., 2017; Ramesh et al., 2005). Homologues of the cell and nuclear fusion proteins HAP2 and GEX1, respectively, have also been found in ancestral trypanosomatid lineages and extant Leishmania species, providing indirect evidence that sex was already present in the last eukaryotic common ancestor (Speijer et al., 2015). Importantly, there has been no direct evidence that any of these genes or their products are involved in Leishmania genetic exchange, and due to the difficulty of making direct observations of the mating events in the sand fly, the precise nature of this reproductive process, including the identification of putative gametic cells, has been difficult to study.

We recently reported that Leishmania promastigotes, which can be grown readily in axenic culture, can form stable hybrids entirely in vitro (Louradour et al., 2020). While this result proved that mating competent cells can arise in culture, its broader applicability to study the biology of mating and to generate recombinant parasites for genetic linkage analysis was limited by the fact that the in vitro hybridization was confined to only two strains of Leishmania tropica and occurred at frequencies far lower than in the sand fly. In the current studies, we have investigated the conditions that might potentiate the mating success of Leishmania in vitro. It is common for organisms that are facultatively sexual to undergo sex in response to environmental stress, including conditions that produce DNA damage (Ram and Hadany, 2016; Schoustra et al., 2010). Thus, a primary function of sex is proposed to be an adaptation for DNA repair, for which a homologous, undamaged chromosome serves as a template (Bernstein et al., 2018). Amongst eukaryotic microbes, oxidative DNA-damaging conditions have been shown to promote sexual spore formation in the yeast Schizosaccharomyces pombe (Bernstein and Johns, 1989) and in the multicellular green alga Volvox carteri (Nedelcu et al., 2004). Oxidative stress also efficiently induced the expression of sexual pheromone precursors and same-sex mating in Candida albicans (Guan et al., 2019). DNA damage caused by X-irradiation induced meiotic recombination in the budding yeast Sacchromyces cerevisiae (Thorne and Byers, 1993) and in the nematode Caenorhabditis elegans (Dernburg et al., 1998). In the current studies, we have used DNA-damaging conditions to promote in vitro hybridization of Leishmania. Exposure to H₂O_2, methyl methanesulfonate (MMS), or to γ-radiation produced a remarkably enhanced efficiency in the in vitro hybridization of the L. tropica strains, and extended to other species and strains a capacity for in vitro hybridization. Single-cell RNA-seq comparisons of untreated and irradiated promastigotes revealed greatly expanded subpopulations amongst the irradiated cells that expressed homologues of genes involved in genetic exchange. We identified expression of the ancestral gamete fusogen HAP2 as a selectable marker for promastigotes able to hybridize in vitro.

Reproduction in Leishmania includes a cryptic sexual cycle that operates among the extracellular promastigote stages developing in the sand fly vector. The mating-competent cells, and in particular the existence of a gametic stage, have yet to be identified. Studying the biology of mating in Leishmania would be greatly facilitated if these events could be reproduced using cells from axenic cultures, and we have recently achieved the recovery of stable hybrids generated between cultured promastigotes entirely in vitro (Louradour et al., 2020). While these initial findings established that mating-competent forms can arise in axenic culture, the mating frequencies were extremely low, and only a single pairwise combination of L. tropica parental lines would successfully hybridize. In the current studies, we demonstrate that under various conditions of environmental stress, including exposure to H₂O₂, MMS, or γ-radiation, the in vitro hybridization capacity of the L. tropica strains could be markedly enhanced and extended to other L. tropica strains and multiple Leishmania species, including L. infantum, L. donovani, and L. braziliensis. The products of these cell fusion events were in each case full-genome, polyploid hybrids. By single-cell RNA-seq analysis, we could identify transcriptionally unique populations of promastigotes whose numbers were drastically increased in each parental line following exposure to γ-radiation, and for which specific meiotic gene homologues, including the ancestral gamete fusogen HAP2, were found to be upregulated. By generating reporter constructs for HAP2, we could select for promastigotes that could hybridize or not in vitro.

Our findings add to the list of facultatively sexual eukaryotes, including yeast and green alga, that exhibit more sexual events under stress conditions. As each of the stress conditions that we imposed are known to produce DNA damage, the results also support studies in varied model eukaryotes, suggesting that DNA repair is a central adaptive function of meiotic recombination (Bernstein et al., 2018; Ram and Hadany, 2016; Schoustra et al., 2010). Homologous chromosome pairing during meiosis provides an undamaged, template chromosome for DNA repair. The generation and repair of double-strand breaks is evidenced in our in vitro hybrids by the substantial number of chromosomes for which there were changes within the chromosome in the relative contribution of the parental alleles. A role for SPO11 protein in the induction of programmed DNA double-strand breaks is argued to be a universal feature of meiotic crossover events in eukaryotes (Keeney et al., 1997). However, meiotic recombination induced by X-irradiation in S. cerevisiae and in C. elegans was shown to occur even when SPO11 was absent or nonfunctional (Dernburg et al., 1998; Thorne and Byers, 1993). In our studies, expression of the homologous SPO11 gene in L. tropica was not upregulated in the irradiated, hybridization-competent cells.

In order to identify hybridization-competent cells present in Leishmania cultures, we performed single-cell RNA sequencing of L747 and MA37 cultured cells exposed or not to γ-radiation. To our knowledge, this is the first application of single-cell transcriptomics in Leishmania, although it has been applied to other kinetoplastid pathogens to reveal the heterogeneity of the parasite populations arising in their mammalian hosts or insect vectors (reviewed in Briggs et al., 2021). Our analyses revealed transcriptomic heterogeneity within and between the parental L. tropica lines in the promastigote cultures. The heterogeneity likely reflects cells in different phases of the cell cycle, as well as differences in the developmental stage progression of the cultured promastigotes. The extensive differences in the transcriptomes of cells comparing the parental lines, apparent in both the PCA and UMAP visualizations, are surprising given the genetic similarity between these L. tropica strains (Iantorno et al., 2017) and likely reflect their different culture histories.

More meaningful to the current studies is the cluster heterogeneity comparing the untreated versus irradiated cells, which for each parental line revealed the emergence of a transcriptionally unique population in the irradiated cultures. Probing the clusters for enriched transcripts related to meiotic processes identified homologues of HAP2, GEX1, and RAD51, involved respectively in cell fusion, nuclear fusion, and DNA repair, that were largely confined to the unique clusters that emerged in the irradiated cultures. Even though a large fraction of the transcriptional shift is likely explained as part of a more global cellular mechanism to recover from stress caused by the irradiation, the expansion of cells expressing HAP2 transcripts in treated samples, the co-regulation of other meiotic genes (a HAP2-like FusM homologue, GEX1, and RAD51), and validation in vitro using HAP2 reporters, all strongly suggest that these cells are associated with the enhanced hybridization frequency observed post-treatment. The possible role in genetic exchange of other genes coordinately expressed with HAP2 but without any described meiotic function will be investigated in future studies. Of note, RNA-binding proteins, such as RBP5 and ZC3H11, were found to be upregulated in these cells and could have novel regulatory functions during Leishmania mating. Critically, we did not find significant upregulation of other meiosis-specific markers, including SPO11, DMC1, and MND1. This suggests that cell fusion might be a response to DNA stress that is independent of a meiotic program, especially as the meiotic reduction steps were for the most part absent in the stress-induced hybrids. Nevertheless, a similar finding was recently reported in an scRNA-seq analysis of T. brucei stages in the tsetse fly salivary glands, in which the gamete cell cluster was defined based solely on the upregulation of HOP1 and HAP2 (Hutchinson et al., 2021). It is possible that some of these markers of gametogenesis are expressed within a narrow time frame, and by analyzing single-cell gene expression at a single time point post-irradiation, we may have missed their expression peak in what is likely a rare cell population. In addition, expression of several genes may not have been detected due to ‘dropout’ events or excessive zeros, which has often been described in this type of approach due to a lower sequencing depth when compared to bulk RNA-seq (Choi et al., 2020).

HAP2 is a type I transmembrane protein that displays the same three-dimensional fold as class II viral fusion proteins. It plays a key role in the membrane fusion of gametes or mating cells that was first identified in plants, and is broadly conserved from protists to invertebrates (Cole et al., 2014; Fédry et al., 2017; Johnson et al., 2004; Liu et al., 2008; Mori et al., 2006). HAP2 expression was recently reported in T. brucei, both in gametes and in meiotic intermediates (Peacock et al., 2021). In the current work, we generated reporter constructs for HAP2 in L747 and MA37 and could show enhanced expression or altered expression kinetics in the stressed cultures. By sorting the cells into mNG⁺ and mNG^- populations, we could confirm that HAP2 expression is a marker for hybridization-competent cells since hybridization failed unless the co-cultures included mNG⁺ selected cells from one or both of the parents. To the extent that HAP2 is itself required for hybridization, then the apparent requirement for HAP2 expression on just one of the fusing cells is consistent with a number of studies showing that it can function unilaterally; either on the mating type minus gamete in Chlamydomonas algae (Fédry et al., 2017; Liu et al., 2008) or on the male gamete in Plasmodium (Liu et al., 2008), the ciliated protozoan Tetrahymena thermophila (Cole et al., 2014), or the plant Arabidopsis thaliana (Johnson et al., 2004; Mori et al., 2006). It is interesting that the MA37 parent showed a relatively high frequency of HAP2-expressing cells even absent exposure to the exogenous stress conditions. This may explain why this L. tropica strain hybridizes with such high efficiency in flies (Inbar et al., 2019), and why we are able to obtain some in vitro hybrids even without exogenous stress so long as MA37 is a fusion partner.

Contrary to the in vitro LMA hybrids that were previously generated, which were diploid, triploid, and tetraploid (Louradour et al., 2020), the irradiation-facilitated LMA hybrids, as well as the other intra- and interspecies hybrids described in this report, are in their great majority tetraploid. These results strongly suggest fusion between the two diploid parental cells, and that the stress conditions that promoted the fusogenic potential of the cells may have produced unusual features leading to unreduced gamete formation, which is a common mechanism underlying polyploidization in plants (De Storme and Geelen, 2013). Attempts to condition the tetraploid hybrids to complete a putative meiotic program by passing them through flies or by re-exposure to stress conditions failed to reduce their ploidy, showing that their polyploid state is relatively stable. As tetraploid progeny are occasional products of hybridization events in flies, it is also possible that tetraploid intermediates, as described for some fungi (Heitman, 2010; Roman et al., 1955), are a normal part of the mating process in Leishmania, with the ploidy reduction steps somehow lacking in the in vitro protocol. Since in the related protist T. brucei, haploid gametic cells have been identified prior to fusion in cells localized to the salivary gland of the tsetse fly vector (Peacock et al., 2014), we favor a comparable sexual process can operate during Leishmania development in sand flies. The fact that a few triploid hybrids, also observed in flies, were generated by the in vitro protocol could mean that a haploid gamete was produced by at least one of the parents and may encourage additional manipulations of the culture conditions to promote a more complete and uniform sexual cycle in vitro.

The absence of random segregation and genome-wide recombinations associated with a reductional meiotic program will undermine the utility of the hybrids for use in SNP association studies. It is interesting, nonetheless, that the hybrids seem to have retained an ability for homologue pairing and recombination after fusion, as evidenced by the breaks within chromosomes in the relative contributions of the parental alleles. Since the maintenance of the tetraploid state (two genome complements from each parent) would be expected to mask recombinations generated by reciprocal crossover events typically associated with meiosis I, the data suggest noncrossovers or gene conversion as a mechanism of DNA repair that would result in gene exchange in only one of the homologs. The asymmetric parental contributions observed both within and between different chromosomes could also have arisen by chromosome loss or gain at the level of individual somy, which is well described for Leishmania genomes during vegetative growth (Dujardin et al., 2007; Sterkers et al., 2011). Aneuploidy, including mosaic aneuploidy, is a constitutive feature of the Leishmania genome (Sterkers et al., 2012). Since gene dosage can directly control the levels of gene expression in Leishmania and impact phenotypes (Iantorno et al., 2017; Prieto Barja et al., 2017), the chromosome-specific patterns of aneuploidy and recombination in the hybrids can be used for linkage studies. It is notable that the progeny that were recovered are in almost every case distinct from the other progeny clones generated in the same experimental cross with respect to the combinations of the relative parental contributions to individual chromosomes that they possess. The distinct haplotypes observed in each of the 17 hybrids generated between the L. donovani strains from India and Sri Lanka are a case in point. As these parental strains produce, respectively, visceral and cutaneous forms of disease in humans and in mice (Kariyawasam et al., 2018; Karunaweera et al., 2003), it may be possible to map the genes controlling the behavior of the progeny clones in mice by genome-wide association studies. At a minimum, the admixture of these genomes should permit assessment as to whether or not the cutaneous or visceral tropisms are expressed as a dominant trait.

The stress-facilitated, in vitro hybridization protocol described in this report removes the requirement for sand flies as the main constraint to generating large numbers of genomic hybrids between Leishmania species and strains, and offers a transformative approach for genetic analysis to understand the extraordinary phenotypic diversity of the genus.

The raw sequence data containing reads from the 51 WGS samples and 8 scRNA-seq samples sequenced are deposited in the SRA database with Accession numbers PRJNA756557 and PRJNA756571, respectively. The source data files for Figures 1 and 6 are provided. Summary statistics on the sequencing data are available in Supplementary files 6 and 7.

1. 1. Sacks D

   (2022) NCBI BioProject

   ID PRJNA756557. raw sequence data containing reads from the 51 WGS samples.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA756557
2. 1. Sacks D

   (2022) NCBI BioProject

   ID PRJNA756571. scRNA-seq samples sequenced.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA756571

1. 1. Akopyants NS
   2. Kimblin N
   3. Secundino N
   4. Patrick R
   5. Peters N
   6. Lawyer P
   7. Dobson DE
   8. Beverley SM
   9. Sacks DL

   (2009) Demonstration of genetic exchange during cyclical development of Leishmania in the sand fly vector

   Science 324:265–268.

   https://doi.org/10.1126/science.1169464

   • PubMed
   • Google Scholar
2. 1. Alves CL
   2. Repolês BM
   3. da Silva MS
   4. Mendes IC
   5. Marin PA
   6. Aguiar PHN
   7. Santos SS
   8. Franco GR
   9. Macedo AM
   10. Pena SDJ
   11. Andrade LO
   12. Guarneri AA
   13. Tahara EB
   14. Elias MC
   15. Machado CR

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

   Scientific Reports 8:13335.

   https://doi.org/10.1038/s41598-018-31541-z

   • PubMed
   • Google Scholar
3. 1. Bañuls AL
   2. Guerrini F
   3. Le Pont F
   4. Barrera C
   5. Espinel I
   6. Guderian R
   7. Echeverria R
   8. Tibayrenc M

   (1997) Evidence for hybridization by multilocus enzyme electrophoresis and random amplified polymorphic DNA between Leishmania braziliensis and Leishmania panamensis/guyanensis in Ecuador

   The Journal of Eukaryotic Microbiology 44:408–411.

   https://doi.org/10.1111/j.1550-7408.1997.tb05716.x

   • PubMed
   • Google Scholar
4. 1. Beneke T
   2. Madden R
   3. Makin L
   4. Valli J
   5. Sunter J
   6. Gluenz E

   (2017) A CRISPR Cas9 high-throughput genome editing toolkit for kinetoplastids

   Royal Society Open Science 4:170095.

   https://doi.org/10.1098/rsos.170095

   • PubMed
   • Google Scholar
5. 1. Bernstein C
   2. Johns V

   (1989) Sexual reproduction as a response to H2O2 damage in Schizosaccharomyces pombe

   Journal of Bacteriology 171:1893–1897.

   https://doi.org/10.1128/jb.171.4.1893-1897.1989

   • PubMed
   • Google Scholar
6. 1. Bernstein H
   2. Bernstein C
   3. Michod RE

   (2018) Sex in microbial pathogens

   Infection, Genetics and Evolution 57:8–25.

   https://doi.org/10.1016/j.meegid.2017.10.024

   • PubMed
   • Google Scholar
7. 1. Bolger AM
   2. Lohse M
   3. Usadel B

   (2014) Trimmomatic: a flexible trimmer for Illumina sequence data

   Bioinformatics 30:2114–2120.

   https://doi.org/10.1093/bioinformatics/btu170

   • PubMed
   • Google Scholar
8. 1. Briggs EM
   2. Warren FSL
   3. Matthews KR
   4. McCulloch R
   5. Otto TD

   (2021) Application of single-cell transcriptomics to kinetoplastid research

   Parasitology 148:1223–1236.

   https://doi.org/10.1017/S003118202100041X

   • PubMed
   • Google Scholar
9. 1. Chagas AC
   2. Oliveira F
   3. Debrabant A
   4. Valenzuela JG
   5. Ribeiro JMC
   6. Calvo E

   (2014) Lundep, a sand fly salivary endonuclease increases Leishmania parasite survival in neutrophils and inhibits XIIa contact activation in human plasma

   PLOS Pathogens 10:e1003923.

   https://doi.org/10.1371/journal.ppat.1003923

   • PubMed
   • Google Scholar
10. 1. Choi K
    2. Chen Y
    3. Skelly DA
    4. Churchill GA

    (2020) Bayesian model selection reveals biological origins of zero inflation in single-cell transcriptomics

    Genome Biology 21:183.

    https://doi.org/10.1186/s13059-020-02103-2

    • PubMed
    • Google Scholar
11. 1. Cole ES
    2. Cassidy-Hanley D
    3. Fricke Pinello J
    4. Zeng H
    5. Hsueh M
    6. Kolbin D
    7. Ozzello C
    8. Giddings T Jr
    9. Winey M
    10. Clark TG

    (2014) Function of the male-gamete-specific fusion protein HAP2 in a seven-sexed ciliate

    Current Biology 24:2168–2173.

    https://doi.org/10.1016/j.cub.2014.07.064

    • PubMed
    • Google Scholar
12. 1. De Pablos LM
    2. Ferreira TR
    3. Walrad PB

    (2016) Developmental differentiation in Leishmania lifecycle progression: post-transcriptional control conducts the orchestra

    Current Opinion in Microbiology 34:82–89.

    https://doi.org/10.1016/j.mib.2016.08.004

    • PubMed
    • Google Scholar
13. 1. De Storme N
    2. Geelen D

    (2013) Sexual polyploidization in plants--cytological mechanisms and molecular regulation

    The New Phytologist 198:670–684.

    https://doi.org/10.1111/nph.12184

    • PubMed
    • Google Scholar
14. 1. Dernburg AF
    2. McDonald K
    3. Moulder G
    4. Barstead R
    5. Dresser M
    6. Villeneuve AM

    (1998) Meiotic recombination in C. elegans initiates by a conserved mechanism and is dispensable for homologous chromosome synapsis

    Cell 94:387–398.

    https://doi.org/10.1016/s0092-8674(00)81481-6

    • PubMed
    • Google Scholar
15. 1. Dillon LAL
    2. Okrah K
    3. Hughitt VK
    4. Suresh R
    5. Li Y
    6. Fernandes MC
    7. Belew AT
    8. Corrada Bravo H
    9. Mosser DM
    10. El-Sayed NM

    (2015) Transcriptomic profiling of gene expression and RNA processing during Leishmania major differentiation

    Nucleic Acids Research 43:6799–6813.

    https://doi.org/10.1093/nar/gkv656

    • PubMed
    • Google Scholar
16. 1. Dujardin JC
    2. Bañuls AL
    3. Llanos-Cuentas A
    4. Alvarez E
    5. DeDoncker S
    6. Jacquet D
    7. Le Ray D
    8. Arevalo J
    9. Tibayrenc M

    (1995) Putative Leishmania hybrids in the Eastern Andean valley of Huanuco, Peru

    Acta Tropica 59:293–307.

    https://doi.org/10.1016/0001-706x(95)00094-u

    • PubMed
    • Google Scholar
17. 1. Dujardin JC
    2. De Doncker S
    3. Jacquet D
    4. Bañuls AL
    5. Balavoine M
    6. Van Bockstaele D
    7. Tibayrenc M
    8. Arevalo J
    9. Le Ray D

    (2007) Clonal propagation and the fast generation of karyotype diversity: An in vitro Leishmania model

    Parasitology 134:33–39.

    https://doi.org/10.1017/S0031182006001156

    • PubMed
    • Google Scholar
18. 1. Fédry J
    2. Liu Y
    3. Péhau-Arnaudet G
    4. Pei J
    5. Li W
    6. Tortorici MA
    7. Traincard F
    8. Meola A
    9. Bricogne G
    10. Grishin NV
    11. Snell WJ
    12. Rey FA
    13. Krey T

    (2017) The Ancient Gamete Fusogen HAP2 Is a Eukaryotic Class II Fusion Protein

    Cell 168:904–915.

    https://doi.org/10.1016/j.cell.2017.01.024

    • PubMed
    • Google Scholar
19. 1. Franssen SU
    2. Durrant C
    3. Stark O
    4. Moser B
    5. Downing T
    6. Imamura H
    7. Dujardin JC
    8. Sanders MJ
    9. Mauricio I
    10. Miles MA
    11. Schnur LF
    12. Jaffe CL
    13. Nasereddin A
    14. Schallig H
    15. Yeo M
    16. Bhattacharyya T
    17. Alam MZ
    18. Berriman M
    19. Wirth T
    20. Schönian G
    21. Cotton JA

    (2020) Global genome diversity of the Leishmania donovani complex

    eLife 9:e51243.

    https://doi.org/10.7554/eLife.51243

    • PubMed
    • Google Scholar
20. 1. Gilabert Carbajo C
    2. Cornell LJ
    3. Madbouly Y
    4. Lai Z
    5. Yates PA
    6. Tinti M
    7. Tiengwe C

    (2021) Novel aspects of iron homeostasis in pathogenic bloodstream form Trypanosoma brucei

    PLOS Pathogens 17:e1009696.

    https://doi.org/10.1371/journal.ppat.1009696

    • PubMed
    • Google Scholar
21. 1. Guan G
    2. Tao L
    3. Yue H
    4. Liang W
    5. Gong J
    6. Bing J
    7. Zheng Q
    8. Veri AO
    9. Fan S
    10. Robbins N
    11. Cowen LE
    12. Huang G

    (2019) Environment-induced same-sex mating in the yeast Candida albicans through the Hsf1-Hsp90 pathway

    PLOS Biology 17:e2006966.

    https://doi.org/10.1371/journal.pbio.2006966

    • PubMed
    • Google Scholar
22. 1. Hafemeister C
    2. Satija R

    (2019) Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression

    Genome Biology 20:296.

    https://doi.org/10.1186/s13059-019-1874-1

    • PubMed
    • Google Scholar
23. 1. Hao Y
    2. Hao S
    3. Andersen-Nissen E
    4. Mauck WM
    5. Zheng S
    6. Butler A
    7. Lee MJ
    8. Wilk AJ
    9. Darby C
    10. Zager M
    11. Hoffman P
    12. Stoeckius M
    13. Papalexi E
    14. Mimitou EP
    15. Jain J
    16. Srivastava A
    17. Stuart T
    18. Fleming LM
    19. Yeung B
    20. Rogers AJ
    21. McElrath JM
    22. Blish CA
    23. Gottardo R
    24. Smibert P
    25. Satija R

    (2021) Integrated analysis of multimodal single-cell data

    Cell 184:3573–3587.

    https://doi.org/10.1016/j.cell.2021.04.048

    • PubMed
    • Google Scholar
24. 1. Heitman J

    (2010) Evolution of eukaryotic microbial pathogens via covert sexual reproduction

    Cell Host & Microbe 8:86–99.

    https://doi.org/10.1016/j.chom.2010.06.011

    • PubMed
    • Google Scholar
25. 1. Hutchinson S
    2. Foulon S
    3. Crouzols A
    4. Menafra R
    5. Rotureau B
    6. Griffiths AD
    7. Bastin P

    (2021) The establishment of variant surface glycoprotein monoallelic expression revealed by single-cell RNA-seq of Trypanosoma brucei in the tsetse fly salivary glands

    PLOS Pathogens 17:e1009904.

    https://doi.org/10.1371/journal.ppat.1009904

    • PubMed
    • Google Scholar
26. 1. Iantorno SA
    2. Durrant C
    3. Khan A
    4. Sanders MJ
    5. Beverley SM
    6. Warren WC
    7. Berriman M
    8. Sacks DL
    9. Cotton JA
    10. Grigg ME

    (2017) Gene Expression in Leishmania Is Regulated Predominantly by Gene Dosage

    MBio 8:e01393-17.

    https://doi.org/10.1128/mBio.01393-17

    • PubMed
    • Google Scholar
27. 1. Inbar E
    2. Akopyants NS
    3. Charmoy M
    4. Romano A
    5. Lawyer P
    6. Elnaiem DEA
    7. Kauffmann F
    8. Barhoumi M
    9. Grigg M
    10. Owens K
    11. Fay M
    12. Dobson DE
    13. Shaik J
    14. Beverley SM
    15. Sacks D

    (2013) The mating competence of geographically diverse Leishmania major strains in their natural and unnatural sand fly vectors

    PLOS Genetics 9:e1003672.

    https://doi.org/10.1371/journal.pgen.1003672

    • PubMed
    • Google Scholar
28. 1. Inbar E
    2. Hughitt VK
    3. Dillon LAL
    4. Ghosh K
    5. El-Sayed NM
    6. Sacks DL

    (2017) The Transcriptome of Leishmania major Developmental Stages in Their Natural Sand Fly Vector

    MBio 8:e00029-17.

    https://doi.org/10.1128/mBio.00029-17

    • PubMed
    • Google Scholar
29. 1. Inbar E
    2. Shaik J
    3. Iantorno SA
    4. Romano A
    5. Nzelu CO
    6. Owens K
    7. Sanders MJ
    8. Dobson D
    9. Cotton JA
    10. Grigg ME
    11. Beverley SM
    12. Sacks D

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

    PLOS Genetics 15:e1008042.

    https://doi.org/10.1371/journal.pgen.1008042

    • PubMed
    • Google Scholar
30. 1. Johnson MA
    2. von Besser K
    3. Zhou Q
    4. Smith E
    5. Aux G
    6. Patton D
    7. Levin JZ
    8. Preuss D

    (2004) Arabidopsis hapless mutations define essential gametophytic functions

    Genetics 168:971–982.

    https://doi.org/10.1534/genetics.104.029447

    • PubMed
    • Google Scholar
31. 1. Jones NG
    2. Thomas EB
    3. Brown E
    4. Dickens NJ
    5. Hammarton TC
    6. Mottram JC

    (2014) Regulators of Trypanosoma brucei cell cycle progression and differentiation identified using a kinome-wide RNAi screen

    PLOS Pathogens 10:e1003886.

    https://doi.org/10.1371/journal.ppat.1003886

    • PubMed
    • Google Scholar
32. 1. Kariyawasam K
    2. Selvapandiyan A
    3. Siriwardana H
    4. Dube A
    5. Karunanayake P
    6. Senanayake S
    7. Dey R
    8. Gannavaram S
    9. Nakhasi HL
    10. Karunaweera ND

    (2018) Dermotropic Leishmania donovani in Sri Lanka: visceralizing potential in clinical and preclinical studies

    Parasitology 145:443–452.

    https://doi.org/10.1017/S003118201700169X

    • PubMed
    • Google Scholar
33. 1. Karunaweera ND
    2. Pratlong F
    3. Siriwardane HVYD
    4. Ihalamulla RL
    5. Dedet JP

    (2003) Sri Lankan cutaneous leishmaniasis is caused by Leishmania donovani zymodeme MON-37

    Transactions of the Royal Society of Tropical Medicine and Hygiene 97:380–381.

    https://doi.org/10.1016/s0035-9203(03)90061-7

    • PubMed
    • Google Scholar
34. 1. Keeney S
    2. Giroux CN
    3. Kleckner N

    (1997) Meiosis-specific DNA double-strand breaks are catalyzed by Spo11, a member of a widely conserved protein family

    Cell 88:375–384.

    https://doi.org/10.1016/s0092-8674(00)81876-0

    • PubMed
    • Google Scholar
35. 1. Kramer S
    2. Queiroz R
    3. Ellis L
    4. Hoheisel JD
    5. Clayton C
    6. Carrington M

    (2010) The RNA helicase DHH1 is central to the correct expression of many developmentally regulated mRNAs in trypanosomes

    Journal of Cell Science 123:699–711.

    https://doi.org/10.1242/jcs.058511

    • PubMed
    • Google Scholar
36. 1. Krzywinski M
    2. Schein J
    3. Birol I
    4. Connors J
    5. Gascoyne R
    6. Horsman D
    7. Jones SJ
    8. Marra MA

    (2009) Circos: an information aesthetic for comparative genomics

    Genome Research 19:1639–1645.

    https://doi.org/10.1101/gr.092759.109

    • PubMed
    • Google Scholar
37. 1. Liu Y
    2. Tewari R
    3. Ning J
    4. Blagborough AM
    5. Garbom S
    6. Pei J
    7. Grishin NV
    8. Steele RE
    9. Sinden RE
    10. Snell WJ
    11. Billker O

    (2008) The conserved plant sterility gene HAP2 functions after attachment of fusogenic membranes in Chlamydomonas and Plasmodium gametes

    Genes & Development 22:1051–1068.

    https://doi.org/10.1101/gad.1656508

    • PubMed
    • Google Scholar
38. 1. Louradour I
    2. Ferreira TR
    3. Ghosh K
    4. Shaik J
    5. Sacks D

    (2020) In Vitro Generation of Leishmania Hybrids

    Cell Reports 31:107507.

    https://doi.org/10.1016/j.celrep.2020.03.071

    • PubMed
    • Google Scholar
39. 1. Love MI
    2. Huber W
    3. Anders S

    (2014) Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2

    Genome Biology 15:550.

    https://doi.org/10.1186/s13059-014-0550-8

    • PubMed
    • Google Scholar
40. 1. Mori T
    2. Kuroiwa H
    3. Higashiyama T
    4. Kuroiwa T

    (2006) GENERATIVE CELL SPECIFIC 1 is essential for angiosperm fertilization

    Nature Cell Biology 8:64–71.

    https://doi.org/10.1038/ncb1345

    • PubMed
    • Google Scholar
41. 1. Nedelcu AM
    2. Marcu O
    3. Michod RE

    (2004) Sex as a response to oxidative stress: a twofold increase in cellular reactive oxygen species activates sex genes

    Proceedings. Biological Sciences 271:1591–1596.

    https://doi.org/10.1098/rspb.2004.2747

    • PubMed
    • Google Scholar
42. 1. Nolder D
    2. Roncal N
    3. Davies CR
    4. Llanos-Cuentas A
    5. Miles MA

    (2007) Multiple hybrid genotypes of Leishmania (viannia) in a focus of mucocutaneous Leishmaniasis

    The American Journal of Tropical Medicine and Hygiene 76:573–578.

    https://doi.org/10.4269/ajtmh.2007.76.573

    • PubMed
    • Google Scholar
43. 1. Odiwuor S
    2. De Doncker S
    3. Maes I
    4. Dujardin J-C
    5. Van der Auwera G

    (2011) Natural Leishmania donovani/Leishmania aethiopica hybrids identified from Ethiopia

    Infection, Genetics and Evolution 11:2113–2118.

    https://doi.org/10.1016/j.meegid.2011.04.026

    • PubMed
    • Google Scholar
44. 1. Peacock L
    2. Bailey M
    3. Carrington M
    4. Gibson W

    (2014) Meiosis and haploid gametes in the pathogen Trypanosoma brucei

    Current Biology 24:181–186.

    https://doi.org/10.1016/j.cub.2013.11.044

    • PubMed
    • Google Scholar
45. 1. Peacock L
    2. Kay C
    3. Farren C
    4. Bailey M
    5. Carrington M
    6. Gibson W

    (2021) Sequential production of gametes during meiosis in trypanosomes

    Communications Biology 4:555.

    https://doi.org/10.1038/s42003-021-02058-5

    • PubMed
    • Google Scholar
46. 1. Prieto Barja P
    2. Pescher P
    3. Bussotti G
    4. Dumetz F
    5. Imamura H
    6. Kedra D
    7. Domagalska M
    8. Chaumeau V
    9. Himmelbauer H
    10. Pages M
    11. Sterkers Y
    12. Dujardin J-C
    13. Notredame C
    14. Späth GF

    (2017) Haplotype selection as an adaptive mechanism in the protozoan pathogen Leishmania donovani

    Nature Ecology & Evolution 1:1961–1969.

    https://doi.org/10.1038/s41559-017-0361-x

    • PubMed
    • Google Scholar
47. 1. Ram Y
    2. Hadany L

    (2016) Condition-dependent sex: who does it, when and why?

    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 371:20150539.

    https://doi.org/10.1098/rstb.2015.0539

    • PubMed
    • Google Scholar
48. 1. Ramesh MA
    2. Malik S-B
    3. Logsdon JM Jr

    (2005) A phylogenomic inventory of meiotic genes; evidence for sex in Giardia and an early eukaryotic origin of meiosis

    Current Biology 15:185–191.

    https://doi.org/10.1016/j.cub.2005.01.003

    • PubMed
    • Google Scholar
49. 1. Ravel C
    2. Cortes S
    3. Pratlong F
    4. Morio F
    5. Dedet J-P
    6. Campino L

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

    International Journal for Parasitology 36:1383–1388.

    https://doi.org/10.1016/j.ijpara.2006.06.019

    • PubMed
    • Google Scholar
50. 1. Robinson JT
    2. Thorvaldsdóttir H
    3. Winckler W
    4. Guttman M
    5. Lander ES
    6. Getz G
    7. Mesirov JP

    (2011) Integrative genomics viewer

    Nature Biotechnology 29:24–26.

    https://doi.org/10.1038/nbt.1754

    • PubMed
    • Google Scholar
51. 1. Rogers MB
    2. Downing T
    3. Smith BA
    4. Imamura H
    5. Sanders M
    6. Svobodova M
    7. Volf P
    8. Berriman M
    9. Cotton JA
    10. Smith DF

    (2014) Genomic confirmation of hybridisation and recent inbreeding in a vector-isolated Leishmania population

    PLOS Genetics 10:e1004092.

    https://doi.org/10.1371/journal.pgen.1004092

    • PubMed
    • Google Scholar
52. 1. Roman H
    2. Phillips MM
    3. Sands SM

    (1955) Studies of Polyploid Saccharomyces. I. Tetraploid Segregation

    Genetics 40:546–561.

    https://doi.org/10.1093/genetics/40.4.546

    • PubMed
    • Google Scholar
53. 1. Romano A
    2. Inbar E
    3. Debrabant A
    4. Charmoy M
    5. Lawyer P
    6. Ribeiro-Gomes F
    7. Barhoumi M
    8. Grigg M
    9. Shaik J
    10. Dobson D
    11. Beverley SM
    12. Sacks DL

    (2014) Cross-species genetic exchange between visceral and cutaneous strains of Leishmania in the sand fly vector

    PNAS 111:16808–16813.

    https://doi.org/10.1073/pnas.1415109111

    • PubMed
    • Google Scholar
54. 1. Sadlova J
    2. Yeo M
    3. Seblova V
    4. Lewis MD
    5. Mauricio I
    6. Volf P
    7. Miles MA

    (2011) Visualisation of Leishmania donovani fluorescent hybrids during early stage development in the sand fly vector

    PLOS ONE 6:e19851.

    https://doi.org/10.1371/journal.pone.0019851

    • PubMed
    • Google Scholar
55. 1. Schoustra S
    2. Rundle HD
    3. Dali R
    4. Kassen R

    (2010) Fitness-associated sexual reproduction in a filamentous fungus

    Current Biology 20:1350–1355.

    https://doi.org/10.1016/j.cub.2010.05.060

    • PubMed
    • Google Scholar
56. 1. Shaik JS
    2. Dobson DE
    3. Sacks DL
    4. Beverley SM

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

    Genes 12:167.

    https://doi.org/10.3390/genes12020167

    • PubMed
    • Google Scholar
57. 1. Speijer D
    2. Lukeš J
    3. Eliáš M

    (2015) Sex is a ubiquitous, ancient, and inherent attribute of eukaryotic life

    PNAS 112:8827–8834.

    https://doi.org/10.1073/pnas.1501725112

    • PubMed
    • Google Scholar
58. 1. Sterkers Y
    2. Lachaud L
    3. Crobu L
    4. Bastien P
    5. Pagès M

    (2011) FISH analysis reveals aneuploidy and continual generation of chromosomal mosaicism in Leishmania major

    Cellular Microbiology 13:274–283.

    https://doi.org/10.1111/j.1462-5822.2010.01534.x

    • PubMed
    • Google Scholar
59. 1. Sterkers Y
    2. Lachaud L
    3. Bourgeois N
    4. Crobu L
    5. Bastien P
    6. Pagès M

    (2012) Novel insights into genome plasticity in Eukaryotes: mosaic aneuploidy in Leishmania

    Molecular Microbiology 86:15–23.

    https://doi.org/10.1111/j.1365-2958.2012.08185.x

    • PubMed
    • Google Scholar
60. 1. Thorne LW
    2. Byers B

    (1993) Stage-specific effects of X-irradiation on yeast meiosis

    Genetics 134:29–42.

    https://doi.org/10.1093/genetics/134.1.29

    • PubMed
    • Google Scholar
61. 1. Tibayrenc M
    2. Ben Abderrazak S
    3. Guerrini F
    4. Bañuls A

    (1993)

    Leishmania and the clonal theory of parasitic protozoa

    Archives de l’Institut Pasteur de Tunis 70:375–382.

    • Google Scholar
62. 1. Van den Broeck F
    2. Savill NJ
    3. Imamura H
    4. Sanders M
    5. Maes I
    6. Cooper S
    7. Mateus D
    8. Jara M
    9. Adaui V
    10. Arevalo J
    11. Llanos-Cuentas A
    12. Garcia L
    13. Cupolillo E
    14. Miles M
    15. Berriman M
    16. Schnaufer A
    17. Cotton JA
    18. Dujardin JC

    (2020) Ecological divergence and hybridization of Neotropical Leishmania parasites

    PNAS 117:25159–25168.

    https://doi.org/10.1073/pnas.1920136117

    • PubMed
    • Google Scholar

Author details

1. Isabelle Louradour

   Laboratory of Parasitic Diseases, National Institute of Allergy and Infectious Diseases, Bethesda, United States

   Present address

   Department of Parasites and Insect Vectors, Institut Pasteur/INSERM U1201, Paris, France

   Contribution

   Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing

   Contributed equally with

   Tiago Rodrigues Ferreira

   Competing interests

   No competing interests declared

2. Tiago Rodrigues Ferreira

   Laboratory of Parasitic Diseases, National Institute of Allergy and Infectious Diseases, Bethesda, United States

   Contribution

   Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing

   Contributed equally with

   Isabelle Louradour

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-4894-1515

3. Emma Duge

   Laboratory of Parasitic Diseases, National Institute of Allergy and Infectious Diseases, Bethesda, United States

   Contribution

   Formal analysis, Investigation, Methodology

   Competing interests

   No competing interests declared

4. Nadira Karunaweera

   Department of Parasitology, Faculty of Medicine, University of Colombo, Colombo, Sri Lanka

   Contribution

   Conceptualization, Resources

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0003-3985-1817

5. Andrea Paun

   Laboratory of Parasitic Diseases, National Institute of Allergy and Infectious Diseases, Bethesda, United States

   Contribution

   Investigation

   Competing interests

   No competing interests declared

6. David Sacks

   Laboratory of Parasitic Diseases, National Institute of Allergy and Infectious Diseases, Bethesda, United States

   Contribution

   Conceptualization, Formal analysis, Funding acquisition, Supervision, Writing – original draft, Writing – review and editing

   For correspondence

   dsacks@nih.gov

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-7557-3124

• 1,930

  views

• 308

  downloads

• 26

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.