We report a high-quality hybrid genome assembly and annotation for a Leishmania (Viannia) guyanensis line (M4147) expressing a firefly luciferase reporter gene and bearing the totivirus LRV1 strongly implicated in parasite hypervirulence. This reference genome should facilitate the study of infectivity and pathogenesis of cutaneous and severe leishmaniasis.

RNA interference (RNAi) is a biological process in which RNA molecules are involved in sequence-specific suppression of gene expression, via small RNA triggers derived from double-stranded RNA that can target specific genes; it is a natural process that plays a role in both the regulation of protein synthesis and in immunity. Discovery of RNAi by Fire and Mello in 1998 had a profound impact on unraveling novel aspects of eukaryotic biology. RNA interference (RNAi) has proven to be an immensely useful tool for studying gene function and validation of potential drug targets in almost all organisms. A great advance in parasitic protozoa was achieved by the experimental demonstration of RNAi in Trypanosoma brucei, and in other protists such as Leishmania braziliensis, Entamoeba histolytica and Giardia lamblia/intestinalis. These organisms exhibit numerous differences beyond the core ‘dicer’ and ‘slicer’ activities, thereby expanding knowledge of the evolutionary diversification of this pathway in eukaryotes. When present, RNAi has led to new technologies for engineering powerful and facile knockdowns in gene expression, revolutionizing biomedical research and opening clinical potentialities. In this review, we discuss the distribution of RNAi pathways, their biological roles, and experimental applications in protozoan parasites.

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.

We report the high-quality draft assemblies and gene annotations for 13 species and/or strains of the protozoan parasite genera Leishmania , Endotrypanum , and Crithidia , which span the phylogenetic diversity of the subfamily Leishmaniinae within the kinetoplastid order of the phylum Euglenazoa. These resources will support studies on the origins of parasitism.

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.

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.

Knowledge of viral diversity is expanding greatly, but many lineages remain underexplored. We surveyed RNA viruses in 52 cultured monoxenous relatives of the human parasite Leishmania (Crithidia and Leptomonas), as well as plant-infecting PhytomonasLeptomonas pyrrhocoris was a hotbed for viral discovery, carrying a virus (Leptomonas pyrrhocoris ostravirus 1) with a highly divergent RNA-dependent RNA polymerase missed by conventional BLAST searches, an emergent clade of tombus-like viruses, and an example of viral endogenization. A deep-branching clade of trypanosomatid narnaviruses was found, notable as Leptomonas seymouri bearing Narna-like virus 1 (LepseyNLV1) have been reported in cultures recovered from patients with visceral leishmaniasis. A deep-branching trypanosomatid viral lineage showing strong affinities to bunyaviruses was termed "Leishbunyavirus" (LBV) and judged sufficiently distinct to warrant assignment within a proposed family termed "Leishbunyaviridae" Numerous relatives of trypanosomatid viruses were found in insect metatranscriptomic surveys, which likely arise from trypanosomatid microbiota. Despite extensive sampling we found no relatives of the totivirus Leishmaniavirus (LRV1/2), implying that it was acquired at about the same time the Leishmania became able to parasitize vertebrates. As viruses were found in over a quarter of isolates tested, many more are likely to be found in the >600 unsurveyed trypanosomatid species. Viral loss was occasionally observed in culture, providing potentially isogenic virus-free lines enabling studies probing the biological role of trypanosomatid viruses. These data shed important insights on the emergence of viruses within an important trypanosomatid clade relevant to human disease.

In the absence of a copyright statement, users should assume that standard copyright protection applies, unless the article contains an explicit statement to the contrary. In case of doubt, contact the journal publisher to verify the copyright status of an article. Title: Viral discovery and diversity in trypanosomatid protozoa with a focus on relatives of the human parasite Leishmania. Authors: Grybchuk D., Akopyants N.S., Kostygov A.Y., Konovalovas A., Lye L.F., Dobson D.E., Zangger H., Fasel N., Butenko A. , Frolov A.O., Votýpka J., d'AvilaLevy C.M., Kulich P., Moravcová J., Plevka P., Rogozin I.B., Serva S., Lukeš J., Beverley S.M., Yurchenko V.

We report here the sequences for all three segments of a novel RNA virus (LepmorLBV1) from the insect trypanosomatid parasite Leptomonas moramango. This virus belongs to a newly discovered group of bunyavirus-like elements termed Leishbunyaviruses (LBV), the first discovered from protists related to arboviruses infecting humans.

Genome sequences were determined for a novel RNA virus, Leptomonas seymouri Narna-like virus 1 (LepseyNLV1). A 2.9-kb segment encodes an RNA-dependent RNA polymerase (RdRp), while a smaller 1.5-kb segment showed no database search matches. This is the first report of bisegmented Narnaviridae from insect trypanosomatids.

ABSTRACT We describe here a new RNA virus (PserNV1) from the plant protist parasite Phytomonas serpens (family Trypanosomatidae , Kinetoplastida , supergroup Excavata ). The properties of PserNV1 permit assignment to the genus Narnavirus ( Narnaviridae ), the first reported from a host other than fungi or oomycetes.

Background: The Leishmania (Viannia) braziliensis complex is responsible for most cases of New World tegumentary leishmaniasis. This complex includes two closely related species but with different geographic distribution and disease phenotypes, L. (V.) peruviana and L. (V.) braziliensis. However, the genetic basis of these differences is not well understood and the status of L. (V.) peruviana as distinct species has been questioned by some.Here we sequenced the genomes of two L. (V.) peruviana isolates (LEM1537 and PAB-4377) using Illumina high throughput sequencing and performed comparative analyses against the L. (V.) braziliensis M2904 reference genome. Comparisons were focused on the detection of Single Nucleotide Polymorphisms (SNPs), insertions and deletions (INDELs), aneuploidy and gene copy number variations.Results: We found 94,070 variants shared by both L. (V.) peruviana isolates (144,079 in PAB-4377 and 136,946 in LEM1537) against the L. (V.) braziliensis M2904 reference genome while only 26,853 variants separated both L. (V.) peruviana genomes. Analysis in coding sequences detected 26,750 SNPs and 1,513 indels shared by both L. (V.) peruviana isolates against L. (V.) braziliensis M2904 and revealed two L. (V.) braziliensis pseudogenes that are likely to have coding potential in L. (V.) peruviana. Chromosomal read density and allele frequency profiling showed a heterogeneous pattern of aneuploidy with an overall disomic tendency in both L. (V.) peruviana isolates, in contrast with a trisomic pattern in the L. (V.) braziliensis M2904 reference.Read depth analysis allowed us to detect more than 368 gene expansions and 14 expanded gene arrays in L. (V.) peruviana, and the likely absence of expanded amastin gene arrays.Conclusions: The greater numbers of interspecific SNP/indel differences between L. (V.) peruviana and L. (V.) braziliensis and the presence of different gene and chromosome copy number variations support the classification of both organisms as closely related but distinct species.The extensive nucleotide polymorphisms and differences in gene and chromosome copy numbers in L. (V.) peruviana suggests the possibility that these may contribute to some of the unique features of its biology, including a lower pathology and lack of mucosal development.

Cutaneous and mucosal leishmaniasis, caused in South America by Leishmania braziliensis, is difficult to cure by chemotherapy (primarily pentavalent antimonials [Sb(V)]). Treatment failure does not correlate well with resistance in vitro, and the factors responsible for treatment failure in patients are not well understood. Many isolates of L. braziliensis (>25%) contain a double-stranded RNA virus named Leishmaniavirus 1 (LRV1), which has also been reported in Leishmania guyanensis, for which an association with increased pathology, metastasis, and parasite replication was found in murine models. Here we probed the relationship of LRV1 to drug treatment success and disease in 97 L. braziliensis-infected patients from Peru and Bolivia. In vitro cultures were established, parasites were typed as L. braziliensis, and the presence of LRV1 was determined by reverse transcription-polymerase chain reaction, followed by sequence analysis. LRV1 was associated significantly with an increased risk of treatment failure (odds ratio, 3.99; P = .04). There was no significant association with intrinsic Sb(V) resistance among parasites, suggesting that treatment failure arises from LRV1-mediated effects on host metabolism and/or parasite survival. The association of LRV1 with clinical drug treatment failure could serve to guide more-effective treatment of tegumentary disease caused by L. braziliensis.

Genetic exchange between Leishmania major strains during their development in the sand fly vector has been experimentally shown. To investigate the possibility of genetic exchange between different Leishmania species, a cutaneous strain of L. major and a visceral strain of Leishmania infantum, each bearing a different drug-resistant marker, were used to coinfect Lutzomyia longipalpis sand flies. Eleven double-drug-resistant progeny clones, each the product of an independent mating event, were generated and submitted to genotype and phenotype analyses. The analysis of multiple allelic markers across the genome suggested that each progeny clone inherited at least one full set of chromosomes from each parent, with loss of heterozygosity at some loci, and uniparental retention of maxicircle kinetoplast DNA. Hybrids with DNA contents of approximately 2n, 3n, and 4n were observed. In vivo studies revealed clear differences in the ability of the hybrids to produce pathology in the skin or to disseminate to and grow in the viscera, suggesting polymorphisms and differential inheritance of the gene(s) controlling these traits. The studies, to our knowledge, represent the first experimental confirmation of cross-species mating in Leishmania, opening the way toward genetic linkage analysis of important traits and providing strong evidence that genetic exchange is responsible for the generation of the mixed-species genotypes observed in natural populations.

Infection with Leishmania parasites causes mainly cutaneous lesions at the site of the sand fly bite. Inflammatory metastatic forms have been reported with Leishmania species such as L. braziliensis, guyanensis and aethiopica. Little is known about the factors underlying such exacerbated clinical presentations. Leishmania RNA virus (LRV) is mainly found within South American Leishmania braziliensis and guyanensis. In a mouse model of L. guyanensis infection, its presence is responsible for an hyper-inflammatory response driven by the recognition of the viral dsRNA genome by the host Toll-like Receptor 3 leading to an exacerbation of the disease. In one instance, LRV was reported outside of South America, namely in the L. major ASKH strain from Turkmenistan, suggesting that LRV appeared before the divergence of Leishmania subgenera. LRV presence inside Leishmania parasites could be one of the factors implicated in disease severity, providing rationale for LRV screening in L. aethiopica.A new LRV member was identified in four L. aethiopica strains (LRV-Lae). Three LRV-Lae genomes were sequenced and compared to L. guyanensis LRV1 and L. major LRV2. LRV-Lae more closely resembled LRV2. Despite their similar genomic organization, a notable difference was observed in the region where the capsid protein and viral polymerase open reading frames overlap, with a unique -1 situation in LRV-Lae. In vitro infection of murine macrophages showed that LRV-Lae induced a TLR3-dependent inflammatory response as previously observed for LRV1.In this study, we report the presence of an immunogenic dsRNA virus in L. aethiopica human isolates. This is the first observation of LRV in Africa, and together with the unique description of LRV2 in Turkmenistan, it confirmed that LRV was present before the divergence of the L. (Leishmania) and (Viannia) subgenera. The potential implication of LRV-Lae on disease severity due to L. aethiopica infections is discussed.

Invertebrate stages of Leishmania are capable of genetic exchange during their extracellular growth and development in the sand fly vector. Here we explore two variables: the ability of diverse L. major strains from across its natural range to undergo mating in pairwise tests; and the timing of the appearance of hybrids and their developmental stage associations within both natural (Phlebotomus duboscqi) and unnatural (Lutzomyia longipalpis) sand fly vectors. Following co-infection of flies with parental lines bearing independent drug markers, doubly-drug resistant hybrid progeny were selected, from which 96 clonal lines were analyzed for DNA content and genotyped for parent alleles at 4-6 unlinked nuclear loci as well as the maxicircle DNA. As seen previously, the majority of hybrids showed '2n' DNA contents, but with a significant number of '3n' and one '4n' offspring. In the natural vector, 97% of the nuclear loci showed both parental alleles; however, 3% (4/150) showed only one parental allele. In the unnatural vector, the frequency of uniparental inheritance rose to 10% (27/275). We attribute this to loss of heterozygosity after mating, most likely arising from aneuploidy which is both common and temporally variable in Leishmania. As seen previously, only uniparental inheritance of maxicircle kDNA was observed. Hybrids were recovered at similar efficiencies in all pairwise crosses tested, suggesting that L. major lacks detectable 'mating types' that limit free genetic exchange. In the natural vector, comparisons of the timing of hybrid formation with the presence of developmental stages suggest nectomonads as the most likely sexually competent stage, with hybrids emerging well before the first appearance of metacyclic promastigotes. These studies provide an important perspective on the prevalence of genetic exchange in natural populations of L. major and a guide for experimental studies to understand the biology of mating.

Leishmaniasis is a spectral disease caused by the Trypanosomatid protozoan parasite Leishmania . One manifestation is a disfiguring mucocutaneous form (MCL) involving dissemination to nasopharyngeal areas and tissue destruction, caused by species from the South American Leishmania subgenus Viannia . We reported previously that metastatic L. guyanensis bear high levels of a novel dsRNA totivirus (LRV1), while non‐metastatic lines have low levels or lack LRV1. LRV1‐infected Leishmania were associated with a hyper‐inflammatory response dependent on the host Toll‐like receptor 3 (TLR3) signaling. Thus, LRV in metastatic parasites subverts the host immune response and promotes dissemination. This is the first parasite factor implicated causally in phenotypes relevant to metastatic mucocutaneous disease, with LRV acting as a ‘parasite within a parasite’. The closely related species L. braziliensis is responsible for the great majority of human MCL. Thus we are surveying L. braziliensis strains to establish the distribution and potential for clinical association(s) of LRVs in parasites taken from humans. Several novel LRVs were identified; these are currently being characterized by virus and genome sequencing, classical virological approaches, and infectivity tests of the LRV‐bearing L. braziliensis in animal models.

Phlebotomine sand flies that transmit the protozoan parasite Leishmania differ greatly in their ability to support different parasite species or strains in the laboratory: while some show considerable selectivity, others are more permissive. In "selective" sand flies, Leishmania binding and survival in the fly midgut typically depends upon the abundant promastigote surface adhesin lipophosphoglycan (LPG), which exhibits species- and strain-specific modifications of the dominant phosphoglycan (PG) repeat units. For the "selective" fly Phlebotomus papatasi PpapJ, side chain galactosyl-modifications (scGal) of PG repeats play key roles in parasite binding. We probed the specificity and properties of this scGal-LPG PAMP (Pathogen Associated Molecular Pattern) through studies of natural isolates exhibiting a wide range of galactosylation patterns, and of a panel of isogenic L. major engineered to express similar scGal-LPG diversity by transfection of SCG-encoded β1,3-galactosyltransferases with different activities. Surprisingly, both 'poly-scGal' and 'null-scGal' lines survived poorly relative to PpapJ-sympatric L. major FV1 and other 'mono-scGal' lines. However, survival of all lines was equivalent in P. duboscqi, which naturally transmit L. major strains bearing 'null-scGal'-LPG PAMPs. We then asked whether scGal-LPG-mediated interactions were sufficient for PpapJ midgut survival by engineering Leishmania donovani, which normally express unsubstituted LPG, to express a 'PpapJ-optimal' scGal-LPG PAMP. Unexpectedly, these "L. major FV1-cloaked" L. donovani-SCG lines remained unable to survive within PpapJ flies. These studies establish that midgut survival of L. major in PpapJ flies is exquisitely sensitive to the scGal-LPG PAMP, requiring a specific 'mono-scGal' pattern. However, failure of 'mono-scGal' L. donovani-SCG lines to survive in selective PpapJ flies suggests a requirement for an additional, as yet unidentified L. major-specific parasite factor(s). The interplay of the LPG PAMP and additional factor(s) with sand fly midgut receptors may determine whether a given sand fly host is "selective" or "permissive", with important consequences to both disease transmission and the natural co-evolution of sand flies and Leishmania.