Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Viral discovery and diversity in trypanosomatid protozoa with a focus on relatives of the human parasite Leishmania

Danyil Grybchuk, Natalia S. Akopyants, Alexei Y. Kostygov, Aleksandras Konovalovas, Lon-Fye Lye, Deborah E. Dobson, Haroun Zangger, Nicolas Fasel, Anzhelika Butenko, Alexander O. Frolov, Jan Votýpka, Claudia M. d’Avila-Levy, Pavel Kulich, Jana Moravcová, Pavel Plevka, Igor B. Rogozin, View ORCID ProfileSaulius Serva, Julius Lukeš, View ORCID ProfileStephen M. Beverley, and Vyacheslav Yurchenko
  1. aLife Science Research Centre, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic;
  2. bDepartment of Molecular Microbiology, Washington University School of Medicine, Saint Louis, MO 63110;
  3. cZoological Institute of the Russian Academy of Sciences, St. Petersburg, 199034, Russia;
  4. dDepartment of Biochemistry and Molecular Biology, Institute of Biosciences, Life Sciences Center, Vilnius University, Vilnius 10257, Lithuania;
  5. eDepartment of Biochemistry, University of Lausanne, 1066 Epalinges, Switzerland;
  6. fDepartment of Parasitology, Faculty of Science, Charles University, 128 44 Prague, Czech Republic;
  7. gBiology Centre, Institute of Parasitology, Czech Academy of Sciences, 370 05 České Budějovice, Czech Republic;
  8. hColeção de Protozoários, Laboratório de Estudos Integrados em Protozoologia, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, 21040-360 Rio de Janeiro, Brazil;
  9. iVeterinary Research Institute, 621 00 Brno, Czech Republic;
  10. jCentral European Institute of Technology – Masaryk University, 625 00 Brno, Czech Republic;
  11. kNational Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894;
  12. lDepartment of Chemistry and Bioengineering, Faculty of Fundamental Sciences, Vilnius Gediminas Technical University, Vilnius 10223, Lithuania;
  13. mUniversity of South Bohemia, Faculty of Sciences, 370 05 České Budějovice, Czech Republic;
  14. nInstitute of Environmental Technologies, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic

See allHide authors and affiliations

PNAS January 16, 2018 115 (3) E506-E515; first published December 28, 2017; https://doi.org/10.1073/pnas.1717806115
Danyil Grybchuk
aLife Science Research Centre, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Natalia S. Akopyants
bDepartment of Molecular Microbiology, Washington University School of Medicine, Saint Louis, MO 63110;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexei Y. Kostygov
aLife Science Research Centre, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic;
cZoological Institute of the Russian Academy of Sciences, St. Petersburg, 199034, Russia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Aleksandras Konovalovas
dDepartment of Biochemistry and Molecular Biology, Institute of Biosciences, Life Sciences Center, Vilnius University, Vilnius 10257, Lithuania;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Lon-Fye Lye
bDepartment of Molecular Microbiology, Washington University School of Medicine, Saint Louis, MO 63110;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Deborah E. Dobson
bDepartment of Molecular Microbiology, Washington University School of Medicine, Saint Louis, MO 63110;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Haroun Zangger
eDepartment of Biochemistry, University of Lausanne, 1066 Epalinges, Switzerland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nicolas Fasel
eDepartment of Biochemistry, University of Lausanne, 1066 Epalinges, Switzerland;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anzhelika Butenko
aLife Science Research Centre, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alexander O. Frolov
cZoological Institute of the Russian Academy of Sciences, St. Petersburg, 199034, Russia;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jan Votýpka
fDepartment of Parasitology, Faculty of Science, Charles University, 128 44 Prague, Czech Republic;
gBiology Centre, Institute of Parasitology, Czech Academy of Sciences, 370 05 České Budějovice, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Claudia M. d’Avila-Levy
hColeção de Protozoários, Laboratório de Estudos Integrados em Protozoologia, Instituto Oswaldo Cruz, Fundação Oswaldo Cruz, 21040-360 Rio de Janeiro, Brazil;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pavel Kulich
iVeterinary Research Institute, 621 00 Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jana Moravcová
jCentral European Institute of Technology – Masaryk University, 625 00 Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pavel Plevka
jCentral European Institute of Technology – Masaryk University, 625 00 Brno, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Igor B. Rogozin
kNational Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Bethesda, MD 20894;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Saulius Serva
dDepartment of Biochemistry and Molecular Biology, Institute of Biosciences, Life Sciences Center, Vilnius University, Vilnius 10257, Lithuania;
lDepartment of Chemistry and Bioengineering, Faculty of Fundamental Sciences, Vilnius Gediminas Technical University, Vilnius 10223, Lithuania;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Saulius Serva
Julius Lukeš
gBiology Centre, Institute of Parasitology, Czech Academy of Sciences, 370 05 České Budějovice, Czech Republic;
mUniversity of South Bohemia, Faculty of Sciences, 370 05 České Budějovice, Czech Republic;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Stephen M. Beverley
bDepartment of Molecular Microbiology, Washington University School of Medicine, Saint Louis, MO 63110;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Stephen M. Beverley
  • For correspondence: stephen.beverley@wustl.edu vyacheslav.yurchenko@osu.cz
Vyacheslav Yurchenko
aLife Science Research Centre, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic;
gBiology Centre, Institute of Parasitology, Czech Academy of Sciences, 370 05 České Budějovice, Czech Republic;
nInstitute of Environmental Technologies, Faculty of Science, University of Ostrava, 710 00 Ostrava, Czech Republic
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: stephen.beverley@wustl.edu vyacheslav.yurchenko@osu.cz
  1. Contributed by Stephen M. Beverley, November 26, 2017 (sent for review October 11, 2017; reviewed by Said A. Ghabrial, Max L. Nibert, and Larry Simpson)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Largely overlooked, the viruses of protists have started to attract more attention. Several viruses of the family Totiviridae are currently implicated in the increased pathogenicity of parasitic protozoa such as Leishmania to vertebrate hosts. We conducted a broad survey of RNA viruses within trypanosomatids, one of the iconic groups of protists. These revealed several previously unidentified viral taxa including one designated “Leishbunyaviridae” and a highly divergent virus termed “Leptomonas pyrrhocoris ostravirus 1.” Our studies provide important information on the origins as well as the diversity and distribution of viruses within a group of protists related to the human parasite Leishmania.

Abstract

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 Phytomonas. Leptomonas 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.

  • Trypanosomatidae
  • Bunyavirales
  • persistent virus infection
  • coinfection
  • coevolution

The ability of viruses to infect virtually any cellular life form on Earth contributes to their immense diversity. While many eukaryotic groups have been probed for the viral presence, the full diversity of viruses remains to be explored (1). Especially promising is the investigation of RNA viruses in simple eukaryotes such as fungi, green algae, diatoms, slime molds, oomycetes, dinoflagellates, apicomplexans, kinetoplastids, diplomonads, and trichomonads (2⇓–4). While originally considered to be little more than evolutionary curiosities, these viruses have started to attract more attention as their important biological roles are now emerging. For example, Cryphonectria hypovirus 1 plays a key role in limiting pathogenicity to its fungal hosts, with applications toward biological control (5), and several viruses of the family Totiviridae have been implicated in the increased pathogenicity of parasitic protozoa to vertebrate hosts (6, 7).

Most studies reporting unicellular eukaryotic viruses arose from fortuitous discovery of virus-like particles (VLPs) or abundant discrete RNA segments rather than from systematic searches often termed “virus hunting.” Here we present a broad survey of RNA viruses within trypanosomatids, one of the iconic groups of protists. Members of the family Trypanosomatidae exhibit strikingly unusual molecular and biochemical traits (8⇓⇓⇓–12). Several species cause widespread severe illnesses, such as sleeping sickness, Chagas disease, and kala-azar in humans (13). Monoxenous (with one host) parasites of invertebrates (primarily insects) were ancestors of these dixenous (with two hosts) pathogens and still represent the majority of trypanosomatid lineages (14, 15). Phylogenetic analysis of the Trypanosomatidae has shown convincingly that the transition from a monoxenous to a dixenous state occurred at least three times, giving rise to the genera Trypanosoma and Leishmania (both parasites of vertebrates), as well as plant-dwelling Phytomonas (16).

VLPs were reported from a number of trypanosomatid species including Endotrypanum schaudinni, Leishmania hertigi [now classified as Paraleishmania hertigi (17)], Phytomonas spp., Crithidia pragensis, Leptomonas seymouri, Angomonas desouzai, and others (18⇓⇓⇓⇓–23). The molecular era in the research of trypanosomatid viruses began with the pioneering studies of those found in South American Leishmania spp. including Leishmania RNA virus 1 (LRV1) from Leishmania guyanensis and Leishmania braziliensis (24, 25), and an unrelated RNA virus in Phytomonas (21). The biological significance of these lays fallow until the finding that LRV1 was associated with increased disease pathology, parasite numbers, and immune response in animal models (6, 26⇓⇓–29). Subsequent studies provided evidence linking LRV1 to the severity of human leishmaniasis, including acute pathology and drug-treatment failures (30⇓⇓–33), although data relating the viral presence to the chronic mucocutaneous leishmaniasis are mixed (32, 34⇓–36).

Recently, molecular descriptions have been made for the viruses from several additional trypanosomatid species. Among them were a bunyavirus-like virus of Leptomonas moramango (37) as well as narnavirus-like viruses of Leptomonas seymouri (38) and the dixenous plant pathogen Phytomonas serpens (39). Provocatively, Leptomonas seymouri has been recovered from cultures from visceral leishmaniasis patients infected with Leishmania donovani, and many of such Leptomonas seymouri strains bear NLV1 (40). Thus, there appears to be considerable unexplored viral diversity in trypanosomatids, the study of which may contribute to our understanding of the biology of trypanosomatids and their insect and/or plant hosts as well as the origins of viruses in Leishmania.

Results and Discussion

Screening of Trypanosomatid Isolates.

We surveyed 52 isolates including 44 belonging to the genera Crithidia and Leptomonas (subfamily Leishmaniinae), as well as eight belonging to Phytomonas spp. These originated from diverse insect or plant hosts and geographic regions (Table S1). Total RNA from these isolates was digested with S1 nuclease, removing most cellular RNAs, after which the remaining dsRNA arising from dsRNA viruses or replicative intermediates of ssRNA viruses could then be sensitively detected by gel electrophoresis (see Figs. 1A, 3A, and 4A) (41). From this analysis, 11 Leishmaniinae and three Phytomonas spp. exhibited dsRNA bands, while the remainder appeared to lack them (Table 1). Most RNA segments were sequenced, and the sequences of those encoding viral RNA-dependent RNA polymerase (RDRP) were used to assign affiliations to the known viral families.

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Tombus-like virus from Leptomonas pyrrhocoris. (A) Agarose gel electrophoresis of S1-digested total RNAs from strains H10 (lane 1), F19 (lane 2), and F165 (lane 3). LeppyrTLV1 segments are labeled “RNA-T1” and “RNA-T2” on the right and marked by green dots; LeppyrOV1 segments are marked by red dots. The left lane shows a 1-kb DNA ladder. (B) Genome structure of LeppyrTLV1. ORFs for different predicted proteins are shown in different colors. A 127-nt stem-loop is found within the predicted N-terminal region of ORF2. (C) Sequence of the ORF1/2 overlap region including a putative slippery sequence (yellow). The RDRP domain is predicted to start from the ACC coding for threonine as previously reported for the UUUUUA slippery sequence (43). (D) Maximum likelihood phylogenetic tree based on RDRP amino acid sequences. Host taxa are shown by symbols defined in the key for hosts. Numbers at the branches indicate Bayesian posterior probability and maximum likelihood bootstrap supports, respectively; those having a Bayesian posterior probability value of 1.0 and maximum likelihood bootstrap support of 100% are marked with black circles. (The scale bar indicates the number of substitutions per site.) The tree was rooted with the sequences of Nodaviridae. Abbreviations and GenBank accession numbers are given in Tables S2–S4.

View this table:
  • View inline
  • View popup
Table 1.

Virus-positive trypanosomatid isolates

RNA Viruses of Leptomonas pyrrhocoris.

Three of 18 isolates of Leptomonas pyrrhocoris (H10, F165, and F19) originating from various locations worldwide (42) exhibited viral dsRNA bands (Table S1). All three bore two common RNAs of 3.5 and 2.2 kb, termed “RNA-T1” and “RNA-T2” (marked by green dots in Fig. 1A), and two (H10 and F19) contained six additional bands termed “RNAs O1–O6” (marked by red dots in Fig. 1A). Sequence analysis of all RNA segments from H10 and F165 suggested the presence of two viruses. The first was distantly related to Tombusviridae. It comprised RNAs T1 and T2 and was named “Leptomonas pyrrhocoris tombus-like virus 1” (hereafter, “LeppyrTLV1”). The second virus comprising RNAs O1–O6 could not be associated with any of known viral groups and was named “Leptomonas pyrrhocoris ostravirus 1” after the city of Ostrava, where it was discovered (hereafter “LeppyrOV1”). PCR tests confirmed the presence/absence of assignments made by S1 nuclease analysis (Table 1).

LeppyrTLV1.

The sequences of segments T1 and T2 in the strains H10 and F165 were highly similar (96.7 and 97.05% nucleotide identity, respectively). RNA-T1 contained two overlapping ORFs with predicted proteins of 850 and 515 aa (Fig. 1 B and C). For ORF1, a BLAST search in the National Center for Biotechnology Information (NCBI) nonredundant protein database did not yield any hits. The ORF2 showed a clear homology to viral RDRP (cd01699 in the NCBI Conserved Domain Database, CDD) with closest relationships to positive-strand RNA viruses of the Tombusviridae/Nodaviridae group (1). The two ORFs showed an overlap of 880 nt (Fig. 1C). A putative slippery sequence, UUUUUUA, was found 6 nt into the overlap, followed by a 127-nt hairpin 6 nt further. Both elements are typical of the −1 ribosomal frameshift of various viruses (43⇓–45). These data suggest that the RDRP of this virus arises through the synthesis of an N-terminal frameshifted protein. While typical Tombusviridae encode RDRPs translated as a C-terminal extension of an upstream ORF by stop-codon read-through (46), several examples of −1 ribosomal frameshifting have been reported recently (1, 47, 48). RNA-T2 encoded a single ORF (ORF3) with a predicted protein of 455 aa (Fig. 1B), for which no homologs were identified in BLAST database searches.

Neither RNA T1 nor T2 exhibited conserved terminal sequences, which are also absent in both Tombusviridae and Nodaviridae (49, 50). Typically tombusviruses are monopartite, and the members of the related family Nodaviridae have two segments (49, 51). However, recent studies have shown remarkable variation within both groups (1).

Phylogenetic reconstruction using RDRP sequences placed LeppyrTLV1 within a clade distantly related to Tombusviridae, which usually infect plants (Fig. 1D). This clade includes viruses from invertebrates including parasitic nematodes, terrestrial myriapods, bivalves, cephalopods, freshwater crustaceans, and gastropods (Fig. 1D and Table S2) (1). Pyrrhocoris apterus, the firebug host of Leptomonas pyrrhocoris, is known to feed on the corpses of invertebrates (52), suggesting this as a possible route of acquisition.

Endogenous viral element related to LeppyrTLV1.

BLAST searches against the genome assembly of Leptomonas pyrrhocoris H10 (53) revealed that the ORF H10_02_0010 at the rightmost end of the chromosome 2 is homologous to the LeppyrTLV1 RDRP (Fig. S1A). Similar to the RNA-T1 of LeppyrTLV1, an overlapping ORF, H10_02_0020, was found immediately upstream. The overlap contained a potential slippery sequence, GGGAAAU, although we did not detect a stem-loop element thereafter (Fig. S1). The ORF H10_02_0010 and the LeppyrTLV1 RDRP shared 38% overall amino acid identity, including conservation of key RDRP motifs (Fig. S1D). Whole-transcriptome data for Leptomonas pyrrhocoris (53) confirmed transcription of both ORFs. No homology was detected between the ORF1 of LeppyrTLV1 and the predicted chromosomal protein H10_02_0020. We considered the two ORFs of the chromosome 2 as an endogenous viral element (EVE) related to LeppyrTLV1 and named it “LeppyrTLV-EVE1.”

PCR tests with primers specific to LeppyrTLV-EVE1 RDRP revealed its presence in four additional European isolates (P59, LP, PP1, and PP2), all of whose sequences were identical (Tables S3 and S4). In contrast, this EVE1 region differed by 180 nt substitutions (and 84 indels) from the corresponding part of the LeppyrTLV1 RDRP, while the TLV1 RDRP sequences of strains H10 and F165 differed by only seven nucleotide substitutions. The similarity between EVE1 and TLV1 suggests that a TLV1-like RNA was captured via reverse transcription and integration into the Leptomonas pyrrhocoris genome. EVEs occur frequently in evolution and are thought to be mediated primarily by reverse transcriptases encoded in host retroposons (54⇓–56). Indeed, a number of telomere-associated transposable element (TATE) and spliced leader-associated (SLAC) retroelements have been identified in the Leptomonas pyrrhocoris genome (53), including one located immediately upstream of the LeppyrTLV-EVE1 (Fig. S1A). The high level of sequence divergence with LeppyrTLV1 points to a relatively ancient origin of EVE1, perhaps predating the dispersal of Leptomonas pyrrhocoris across Europe (42).

LeppyrOV1.

The six RNAs O1–O6 of strains H10 and F19 (Fig. 1A) were initially viewed as “satellite” RNAs of LeppyrTLV1. However, several observations suggested that they comprise separate virus. First, unlike TLV1 RNAs T1 and T2, the termini of RNAs O1–O6 share common sequences: AAAGAAAAAA at the 5′ and ATGAGTTT at the 3′ ends (defined in the presumptive protein-coding strand orientation) (Fig. 2A). Conserved terminal sequences are known to participate in the replication of viruses and often are defining features of viral families (57). Second, in all strains the ratio of RNAs T1 to T2 was relatively constant; the same was true for RNAs O1–O6. However, the overall ratio of both RNA groups was substantially different.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Leptomonas pyrrhocoris ostravirus 1, a unique virus from Leptomonas pyrrhocoris. (A) Genome structure of LeppyrOV1, showing shared terminal sequences and a single ORF per segment (squiggle marks an incompletely sequenced end). The location of an RDRP domain predicted on RNA-O3 by CDD search, PHYRE2, and HHPred software is shown. (B) Multiple alignments of the LeppyrOV1 putative RDRP with those of Picorna-, Flavi-, and Caliciviridae. Identical residues are shown in red; similar residues are shown in blue. Amino acid motifs, typically found in viral RDRPs, are highlighted in yellow.

Segments O1–O6 each contained a single ORF, and conventional BLAST searches did not yield any homologs for the corresponding hypothetical proteins. However, search algorithms focused on both structural and sequence homology revealed a putative RDRP motif within the predicted 1,315-aa protein within segment O3 (Fig. 2A), albeit with modest statistical support (NCBI CDD, amino acids 767–870, e = 0.89; PHYRE 2, amino acids 684–874, confidence = 56%; HHPRED, amino acids 693–874, confidence = 89.7%) (58, 59). Within this region we identified conserved viral RDRP motifs responsible for catalytic activity and ribonucleotide selectivity (Fig. 2B) (60⇓–62). Analysis of the base frequencies of codon third positions of the viral ORFs showed significant differences between TLV1 and OV1 and a greater degree between these and the nuclear genome of Leptomonas pyrrhocoris (Table S5).

Thus, we conclude that RNAs O1–O6 comprise a previously undescribed virus, Leptomonas pyrrhocoris ostravirus 1 (LeppyrOV1). As yet, we have not found a trypanosomatid strain containing this virus alone, which would firmly establish its independence from LeppyrTLV1. Further studies are required to address the functional relationships between the six segments of this virus and significance of its co-occurrence with LeppyrTLV1.

A Bunyavirus-Like Genus, “Leishbunyavirus.”

Six isolates showed the presence of dsRNAs related to previously described viruses of Leptomonas moramango (37). LepmorLBV1a and b showed features characteristic of many other bunyaviruses, including a trisegmented genome, terminal “panhandle” repeats, and sequence relatedness of the predicted RDRP and nucleocapsid proteins, and were thus assigned as species within the genus Leishbunyavirus (LBV) (37). We confirmed the presence of LepmorLBV1a and 1b in Leptomonas moramango, as well as LBV1s in the dixenous phytopathogenic Phytomonas sp. TCC231 (PTCCLBV1) and four species of Crithidia: Crithidia otongatchiensis (CotoLBV1), Crithidia abscondita (CabsLBV1), Crithidia sp. G15 (CG15LBV1), and Crithidia sp. ZM (CZMLBV1) (Fig. 3A and Table 1). PCR tests with primers complementary to the conserved regions of LBV1 RDRPs showed the presence of these viruses in Crithidia sp. C4 and Crithidia pragensis as well (Table 1). These LBV1-positive strains showed three dsRNAs, except PTCCLBV1 which exhibited only two (Table 1). We sequenced all segments of CotoLBV1 and CabsLBV1 and the largest segment (completely or partially) of the others (Table 1 and Tables S3 and S4).

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Leishbunyaviruses. (A) Viral dsRNA from Crithidia otongatchiensis (lane 1), Crithidia abscondita (lane 2), Crithidia sp. G15 (lane 3), and Crithidia sp. ZM (lane 4). The left lane shows a 1-kb DNA ladder. (B) Genome structure of LBVs. The sizes of segments and their various features (except for terminal complementary sequences) are shown in proportion. EN, endonuclease domain. Orange, teal, and yellow labels in the M segment stand for the signal peptide, glycosylation site(s), and transmembrane domain, respectively. (C) Negative-stain transmission electron micrographs of the virus particle isolated from Crithidia otongatchiensis. (Scale bar, 100 nm.) (D) Maximum likelihood phylogenetic tree based on RDRP amino acid sequences. Numbers at the branches indicate Bayesian posterior probability and maximum likelihood bootstrap supports, respectively; those having a Bayesian posterior probability value of 1.0 and maximum likelihood bootstrap support of 100% are marked with black circles. (The scale bar indicates the number of substitutions per site.) The tree was rooted at the midpoint. Abbreviations and GenBank accession numbers are given in Tables S2–S4. The pictograms describing viral hosts are as in Fig. 1.

Sequence features and coding potential of LBV1s.

Prototypic bunyaviruses bear three RNA segments, termed “large: (L, 7–12 kb), “medium” (M, 3.2–4.9 kb), and “small” (S, 1–3 kb), encoding RDRP, envelope glycoproteins, and nucleocapsid, respectively (57, 63). The corresponding segments in LBV1s were considerably shorter: 6–6.3 kb (L), 1.0–1.9 kb (M), and 0.7–1.0 kb (S) (Table 1). Within each completely sequenced LBV1 segment we identified a single large ORF (Fig. 3B), in contrast to many bunyaviruses, which can encode multiple ORFs on the M and S segments (57). Bunyaviral RNA segments are typically flanked by panhandle inverted repeats mediating key steps of virus replication, transcription, and translation (64). Although the methods used here did not invariably yield full-length sequences, we were often able to identify panhandles in all fragments. In the L and M segments, we identified the sequence ACACAAAG at the 5′ end (as defined by the viral sense orientation) and the complementary sequence CTTTGTGT at the 3′ end. These terminal eight nucleotides are typically found in all viruses belonging to the family Phenuiviridae (Table S6).

The ORFs in all completely sequenced L segments encoded putative proteins of ∼2,000 aa homologous to the bunyaviral RDRP domain (amino acids 600–1,200; pfam04196, e-values <10−10; NCBI CDD). The region between amino acids 85 and 150 showed homology to the bunyaviral endonuclease domain, with conserved key residues involved in Mn2+ ion coordination and phosphodiester bond cleavage (Fig. S2A). Both the RDRP and endonuclease domains are essential for bunyaviral replication (65).

Various database searches (BLAST/CDD, PHYRE2, and HHpred) with the predicted proteins from the four sequenced M segments returned no hits. However, the predicted M proteins displayed a signal peptide as well as varying numbers of transmembrane domains, ranging from two in CabsLBV1, one in CotoLBV1, and zero or one in LepmorLBV1 (CCTOP, TMpred, and TMHMM algorithms), and N-glycosylation sites in CabsLBV1 and CotoLBV1 (Fig. 3B). These analyses suggest that LBV1s, much like other bunyaviruses, are able to exploit cellular secretory systems for glycoprotein synthesis and virion assembly (66, 67). Indeed, purified LBV1 virions visualized by negatively stained transmission electron microscopy displayed the typical envelope with surface projections or spikes spread evenly along its surface (Fig. 3C).

The predicted S segment proteins did not yield compelling BLASTP hits. However, PHYRE structural homology searches showed the similarity of those proteins from CabsLBV1 and LepmorLBV1b to the nucleocapsid proteins of Toscana- and Punta Toro viruses (57.2–89.3% confidence). Alignment of the predicted leishbunyaviral S proteins with the nucleocapsid proteins of other bunyaviruses revealed several universally conserved amino acid motifs (Fig. S2B). Hence, we concluded that LBV1 S segments encode nucleocapsid proteins.

Phylogenetic analysis of LBV1s suggests classification as a family Leishbunyaviridae within the Bunyavirales.

RDRP-based phylogenetic trees showed that the LBV1s formed a well-supported clade separate from other major Bunyavirales groups (Fig. 3D). The closest family was the Phenuiviridae, consistent with the similarities noted earlier in the terminal panhandle elements (Table S6). Many Phenuiviridae have been reported from insects and other arthropods, and viruses within the genus Phlebovirus are transmitted by the same sand fly species as Leishmania (68). However, our data show that LBV1s are far more ancient (Fig. 3D). As the divergence of the LBV1-containing clade from other bunyaviral families is comparable to or greater than other bunyavirus interfamilial divergences, we propose that this clade be recognized as a family, termed Leishbunyaviridae.

Identification of LBVs within metatranscriptomic viral surveys.

Interestingly, BLAST searches with trypanosomatid LBV1 RDRPs identified several hits in the sequences from metatranscriptomic virus-hunting surveys. These included Huangshi Humpbacked Fly virus (HHFV), Wuhan Spider virus (WSV) (69), Hubei bunya-like virus 5 (HBLV5) from a mix of dipterans, Hubei bunya-like virus 6 (HBLV6) from a horse leech (1), and two from honey bees—Apis bunyavirus 1 (ABV1) and Duke bunyavirus (DuBV) (70). On the reconstructed phylogenetic tree all these viruses from metatranscriptomes intermingled with the trypanosomatid LBVs, with high statistical support (Fig. 3D).

Recently, it was proposed that bunyaviruses originated within insects (71, 72), and one explanation for the interdigitation observed here could be multiple transitions of these viruses between arthropods and trypanosomatids. An alternate model is that the insect metatranscriptomic leishbunyaviruses arose not from the insects themselves but from their associated microbiota (73), given that trypanosomatids are well-known parasites of arthropods (14, 74). Thus, we searched the LBV-containing metatranscriptomic sequence read archives (SRAs) for trypanosomatid signatures, a challenging task given the relatively low number of viral reads in these pooled datasets. Nonetheless, BLASTN searches of assembled contigs revealed several abundant trypanosomatid transcripts, such as 18S rRNA or paraflagellar rod proteins (Tables S3 and S4), in the HHFV-, WSV-, HBLV5-, and ABV1-containing SRAs (data for DuBV were not available). Indeed, phylogenetic analysis of these putative transcripts revealed affinities to various trypanosomatids. Based on these data, ABV1 could speculatively be associated with Lotmaria passim (subfamily Leishmaniinae), HHFV and HBLV5 with subfamily Strigomonadinae, and HBLV6 and WSV with the genera Trypanosoma and Herpetomonas, respectively (Tables S2–S4). While the co-occurrence of reads for leishbunyaviruses and trypanosomatids in the metatranscriptomic read sets is not definitive proof that these flagellates actually contained viruses, we consider this a plausible explanation.

These findings provide support for the model postulating a trypanosomatid microbiota origin of LBVs emerging from the metatranscriptome datasets. If borne out, this suggests that, instead of multiple origins from insects, trypanosomatid LBV1s may have originated less frequently and perhaps only once. Consistent with the latter, a significant, albeit imperfect, level of phylogenetic congruency can be seen between trypanosomatid LBVs and nuclear genome phylogenies (Fig. S3). However, the possibility of multiple acquisitions of LBVs by trypanosomatids from insects or other trypanosomatids cannot be formally excluded, given that trypanosomatid LBV1s bear hallmarks of infectious bunyaviruses and reports of mixed trypanosomatid infections (75⇓⇓–78). Currently we favor a model with a single transition of an ancestral insect virus to a trypanosomatid, but further investigations will be required to rigorously establish this hypothesis.

Narnaviridae.

In Leptomonas seymouri and two isolates of Phytomonas serpens we documented the presence of dsRNA (fragments of 2.9 + 1.5 kb and 3.8 kb, respectively) (Fig. 4 A and B and Table 1), in agreement with previous findings (38⇓–40) that these species bear Leptomonas seymouri narna-like virus 1 (LepseyNLV1) and Phytomonas serpens narnavirus 1 (PserNV1).

Fig. 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 4.

Narnaviruses of trypanosomatids. (A, Upper) Viral dsRNA in two subcultures of Phytomonas serpens (isolate 9T). Lane 1 is from the University of California, Riverside and lane 2 is from the Institute of Parasitology in Budweis, Czech Republic. The lane marked “M” is a 1-kb DNA ladder. (Lower) Total RNA was used as a loading control. (B, Upper) Viral dsRNA in four subcultures of Leptomonas seymouri ATCC30220. Lanes 1 and 4 are the original ATCC culture; lanes 2 and 3 are substrains 2003WT and 294-1993VB (Rutgers University), respectively; lane 5 is a culture from Zoological Institute of the Russian Academy of Sciences. Leishmania guyanensis (Lgy, strain M4147) bearing the 5.3-kb LRV1 served as a positive control. The lane marked “M” is a 1-kb DNA ladder. (Lower) Total RNA was used as a loading control for virus-negative substrains. (C) Genome structure of LepseyNLV1 and PserNV1. ORFs for different proteins are shown in different colors. The stem-loop structures and terminal complementary sequences are indicated. Squiggles mark incompletely sequenced ends. (D) Maximum likelihood phylogenetic tree of Narnaviridae based on RDRP amino acid sequences. LepseyNLV1 and PserNV1 are indicated by a trypanosomatid symbol. Numbers at the branches indicate Bayesian posterior probability and maximum likelihood bootstrap supports, respectively; those having a Bayesian posterior probability value of 1.0 and maximum likelihood bootstrap support value of 100%, respectively are marked with black circles. (The scale bar indicates number of substitutions per site.) The tree was rooted with the sequences of Leviviridae. Abbreviations and GenBank accession numbers are given in Tables S2–S4. The pictograms are as in Fig. 1.

Sequence features of trypanosomatid narnaviruses.

The genome of PserNV1 was monosegmented (Fig. 4A, lane 2), and its RNA contained a single ORF for RDRP. In contrast, LepseyNLV1 displayed a bipartite organization (Fig. 4B, lanes 1, 4, and 5) with RNA1 encoding RDRP and RNA2 comprising two overlapping ORFs with no homologs identified in database searches. The region of overlap displayed several structural features associated with +1 ribosomal frameshifting, including a hairpin preceded by a slippery sequence (Fig. S4C), suggesting these two ORFs may be expressed as a fusion protein.

PserNV1 termini were determined by ligating an adapter followed by sequencing across the adapter–virus junction. They revealed features common for Narnaviridae (79): short terminal complementary sequences 5′-ACGC ... GCGT-3′ and putative subterminal hairpin structures (Fig. 4C and Fig. S4A). Intriguingly, the very 5′ end of the viral RNA showed similarity to the spliced leader (SL) of Phytomonas serpens (GenBank X87137). The SL is a 39-nt capped sequence added to the 5′ end of every trypanosomatid mRNA by trans-splicing (80, 81). However, the PserNV1 SL-related sequence lacked the first five nucleotides and had three internal mismatches (Fig. S4B), rendering it unlikely to be functional based on current knowledge of SL function (82). Thus, in the past the PserNV1 may have “snatched” the host’s SL, substituting it for the original terminus. In LepseyNLV1 we did not determine the terminal sequences explicitly; however, typical narnaviral subterminal hairpins were predicted in the RNA2 assembly (Fig. 4C).

Phylogeny and evolutionary origins.

RDRP-based phylogenetic reconstruction showed LepseyNLV1 and PserNV1 to be the closest relatives, forming a well-supported clade along with prototypical narnaviruses—Saccharomyces cerevisiae 20S and 23S viruses (Fig. 4 B and C)—as well as the oomycete-infecting Phytophthora infestans RNA virus 4 (PiRV4) (83). Interestingly, we identified a metatranscriptomic virus from the fly Teleopsis dalmanni (84), whose transcriptome assembly also contained two contigs (GBBP01074304 and GBBP01074305) corresponding to trypanosomatid 18S rRNA genes. We were not able to closely associate these contigs with known trypanosomatid sequences, suggesting that this virus may belong to a yet-uncharacterized lineage. Thus, similar to leishbunyaviruses, the insect metatranscriptomic Narnavirus may have arisen from its trypanosomatid microbiota.

As inferred earlier, Ourmiavirus and ourmia-like viruses (family Ourmiaviridae) clustered preferentially with Narnavirus, while Mitovirus (another genus of Narnaviridae) was sister to the clade comprising those three groups (85, 86). Previous studies suggested that narnaviruses were ancestral parasites of fungi, which later switched to other organisms (87, 88). Yeasts represent a normal component of insect’s intestine, where they could encounter trypanosomatids (89, 90).

While narnaviruses are typically monosegmented, Ourmiaviridae typically contain several segments (87, 88). LepseyNLV1 with its two segments exhibits an independently evolved genome organization being intermediate between those of Narnaviridae and Ourmiaviridae. While definitive evidence that the two segments represent a single virus is lacking, we consider such an association likely, given that both segments are maintained or lost in parallel as described below.

Viral Stability.

In the course of our studies we observed that upon in vitro cultivation some viruses could be occasionally lost. PserNV1 was originally found in the 9T strain from the Czech Republic; however, the same strain maintained elsewhere lacked it (Fig. 4A, lanes 1 and 2). Similarly, while LepseyNLV1 occurred in the ATCC30220 isolate, it was absent in the same strain and in a transfectant derivative obtained from another source (Fig. 4B, lanes 1–3). However, it persisted during continuous cultivation (∼300 passages) in the Zoological Institute of the Russian Academy of Sciences (Fig. 4B, lanes 4 and 5). For leishbunyaviruses, we noticed a gradual decrease of viral dsRNA levels in Crithidia otongatchiensis over 6 mo of continuous cultivation and their disappearance from Crithidia pragensis and Leptomonas moramango after 2 wk of passaging (although low levels could be detected by RT-qPCR). However, no changes in dsRNA abundance were seen for CabsLBV1, CZMLBV1, CG15LBV1, or PTCCLBV1.

Several nonexclusive mechanisms might explain these observations. Some viruses may be intrinsically unstable or lost because the selective pressures on their trypanosomatid hosts may differ in vitro and in vivo. Alternatively, the culture may be heterogeneous in terms of viral presence, and virus-free cells may outcompete their infected counterparts. Our data collectively suggest that caution is warranted when interpreting viral absence in cultured parasites. Serendipitously, virus-free derivatives may serve as isogenic tools for probing potential roles for viruses in parasite biology, as for Leishmania guyanensis LRV1 (6, 91). Indeed, the coincidental loss of LepseyNLV1 RNA1 and RNA2 provides some support for a functional association (Fig. 4B).

Conclusions

Here, we conducted a survey of RNA viruses in two groups of Trypanosomatidae: insect-restricted (monoxenous) relatives of Leishmania (Crithidia and Leptomonas, subfamily Leishmaniinae) and plant-infecting Phytomonas. This greatly expanded the known diversity of RNA viruses in these flagellates, showing that trypanosomatids can be infected by various unrelated viruses: Totiviridae, Narnaviridae, Bunyavirales, tombus-like viruses, and a previously unknown virus that was termed “Ostravirus” and is currently defined by LeppyrOV1, whose RDRP was so divergent that it escaped generic BLAST searches. We also documented EVE formation in trypanosomatids (LeppyrTLV-EVE1), presumably enabled by the activity of the endogenous retroposons.

One interesting question is whether the trypanosomatid viruses can be shed and infect other parasites. Current data suggest that LRV1, like the great majority of other Totiviridae, is not shed (92). Narnaviruses, by virtue of lacking either a capsid or an envelope, are transmitted only vertically or during mating (4). However, the presence of an extra segment in LepseyNLV1 (Fig. 4B) might be associated with transmission, as in related ourmiaviruses (85). Similarly, the two Leptomonas pyrrhocoris viruses OV1 and TLV1 have sufficient coding capacity for transmission. Last, for several trypanosomatid LBV1s we visualized the presence of enveloped virions bearing surface proteins (Fig. 3C), the hallmarks of infectious bunyaviruses. These fascinating questions will be addressed in future studies.

Phylogenetic relationships of relevant trypanosomatid taxa permit a broader view of the origins and evolution of their viruses (Fig. 5). First, Leptomonas pyrrhocoris appears to be a hotbed for viral discovery, with two previously unreported viruses (LeppyrTLV1 and LeppyrOV1) and the presence of an EVE. Secondly, narnaviruses, LBV1s, and LRV1/2s appear to be distributed over the trypanosomatid phylogenetic tree in a patchy manner, with many seemingly virus-free lineages interspersed with ones bearing diverse viruses (Fig. 5). This poses a number of challenges. If one postulates the presence of virus in the common ancestor of a particular group (marked by arrows in Fig. 5), viral loss must have occurred independently in a great many subsequent taxa. Alternatively, if one assumes the common ancestor to be virus-free, independent viral acquisitions must have occurred. The chances of this are speculative at best, perhaps being more likely for the viruses showing increased likelihood for infectivity (LBV1s and, conceivably, OV1, TLV1, and NLV1). Superimposed upon or alternative to this is the possibility of viral exchange via infectious shedding during coinfections, as mixed trypanosomatid infections are quite frequent (see above). Importantly, these latter two processes would be expected to further blur signs of virus–parasite coevolution. Thus, it is remarkable that for LRV1/2 (93) and, to some extent, for LBVs (Fig. S3), the phylogenetic trees for the parasite and their viruses show significant congruency. This suggests that there must be some constraints on the horizontal viral transmission, if present, especially among kingdoms.

Fig. 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 5.

Overview of trypanosomatid relationships and viruses. The evolutionary tree shows the maximum likelihood phylogenetic tree of trypanosomatids reconstructed using 18S rRNA and gGAPDH genes, over which the absence or presence of viruses is marked (see graphical legend). Arrows denote hypothetical acquisition of viruses under the assumption of a single origin in the common ancestor. Maximal bootstrap supports are marked by filled circles, and bootstrap supports over 70% are denoted by open circles. (The scale bar indicates the number of substitutions per site.)

Notably, in our survey we did not find any LRV-related Totiviridae (Fig. 5 and Table 1), although numerous Leishmaniinae were tested. This suggests that these viruses were acquired upon the involvement of vertebrates into the life cycle of Leishmaniinae. Given the elevated pathogenicity of LRV1-bearing Leishmania to the vertebrate host, viral acquisition could be viewed as beneficial for the parasite, if one equates pathogenicity with an increased evolutionary fitness. However, most Leishmania except Viannia and a handful of Leishmania major or Leishmania aethiopica isolates lack LRV1/2 (94). This implies that if LRV1/2 presence was indeed beneficial early in evolution, it became less important in modern lineages and/or was replaced by other mechanisms contributing to virulence, such as type I IFN induction (29).

Variation in the RNAi pathway may contribute to the observed patchiness in viral distribution, as this pathway acts as an antiviral defense mechanism in many species (95). In agreement with this, the RNAi pathway (believed to be ancestral to all eukaryotes) is absent in Phytomonas spp., Leptomonas seymouri, and LRV2-bearing Leishmania aethiopica and Leishmania major (23, 96, 97). The RNAi pathway may be especially important for narnaviruses, which are presumably defenseless because of the lack of capsids. However, LRV1-containing Leishmania guyanensis and Leishmania braziliensis have a highly active RNAi pathway (97), and accordingly LRV1 has mastered the ability to coexist in the face of RNAi attack, although under some circumstances RNAi can lead to its elimination (27). In addition, numerous Crithidia and Leptomonas spp. retain the RNAi pathway (96). It is thus possible that RNAi plays only a weak role in the evolutionary distribution of trypanosomatid viruses.

Several studies have established a role for trypanosomatid viruses in the vertebrate host (6, 27⇓–29, 40, 98). Our studies now suggest that potential role(s) for trypanosomatid viruses in parasite biology within their insect hosts should be considered. While LRV1 and other Totiviridae have been implicated in vertebrate pathogenicity (6, 7, 35), there are no direct data concerning the influence of this virus on the relationships between Leishmania and sand flies. Given that Toll-like receptors were first discovered in insects (99), and TLR3 specifically was implicated in LRV1 pathogenicity (6), this possibility clearly merits attention. Alternatively, viruses may invade and persist as mere parasitic elements rather than providing any advantage to their trypanosomatid hosts. Resolution of these questions may benefit from the serendipitous identification of virus-free isolates of Phytomonas serpens and Leptomonas seymouri and their use in studies assessing potential functional roles.

Previously unidentified viruses were found in considerable numbers in the species/isolates tested (Table 1). The actual diversity of trypanosomatids is not known, but at least 600 species have been already described (100). In addition, the example of Leptomonas pyrrhocoris with its multiple isolates showing variation in viral presence and composition, illustrates another level of diversity. Indeed, as noted in Fig. 5, there are several trypanosomatid lineages for which VLPs have been reported but not studied by modern molecular methods. Furthermore, in various invertebrate metatranscriptomes we found several viruses possibly originating from their trypanosomatid microbiota. Such metatranscriptomes may also provide important new information about the diversity of trypanosomatids themselves. Taken together these findings suggest that a great number of viruses remain to be found in this important group of parasites.

Methods

Isolation of Viral RNA and Primary Screening.

Total RNA was isolated from trypanosomatid cultures using the TRI reagent (MRC Inc.) as described previously (101). For primary screening, 50 µg of total RNA from each sample were treated with RNase-free DNase I (New England Biolabs) and nuclease S1 from Aspergillus oryzae (Sigma-Aldrich) (41). The resulting dsRNA was resolved on 0.8% agarose gel and stained with ethidium bromide. For preparative isolation 400 µg of total RNA from virus-positive cultures were digested with DNase I, followed by ssRNA precipitation by LTS solution (2 M LiCl, 150 mM NaCl, 15 mM Tris HCl, pH 8.0) at 4 °C overnight as described previously (102). The ssRNA fraction was removed by centrifugation for 30 min at 20,000 × g at 4 °C, and dsRNA was precipitated by EtOH and visualized as above. Individual dsRNA bands were gel-purified using Zymoclean Gel RNA Recovery Kit (Zymo Research).

Viral dsRNA Amplification, Cloning, and Sequencing.

Gel-extracted dsRNA was polyadenylated at both 3′ ends using Escherichia coli Poly(A) Polymerase (New England Biolabs) and then purified on a PCR clean-up column (Thermo Fisher Scientific) according to the manufacturer’s protocol. Next, polyadenylated dsRNA was reverse-transcribed using the Transcriptor First Strand cDNA Synthesis Kit (Roche), and an anchored-oligo (dT) primer QD2-T20 5′-ggcaattaaccctcactatagaattcgttcgatctttttttttttttttttttt-3′ (modified from ref. 103). To prevent renaturation of the cRNA strands, DMSO was added to a final concentration of 7.5%, and the residual RNA was lysed with 0.1 M NaOH. The cDNA was then reannealed at 65 °C for 90 min followed by gradual cooling to 4 °C (102). The cDNA was purified on PCR clean-up columns and amplified using Phusion High-Fidelity DNA Polymerase (Thermo Fisher Scientific) with primer QD2, 5′-tcactatagaattcgttcgatc-3′, that anneals to the fragment introduced by QD2-T20. The first PCR step that included end repair (72 °C for 5 min) was followed by the manufacturer’s recommended cycling conditions: 98 °C for 10 s, 55 °C for 30 s, and 72 °C for 40 s/kb. Obtained PCR products were cloned into the pTZ57R vector (Thermo Fisher Scientific) and sequenced by primer walking. For the analysis described above, we were unable to obtain enough dsRNA from Crithidia sp. C4 and Crithidia pragensis. In these cases, the partial RDRP gene (∼900 bp) was amplified using degenerate primers designed to amplify known LBV1s (LeiBunyaF: 5′-ttykcvacnttcaagaaragcac-3′ and LeiBunyaR: 5′-ccagartcatcwgadgadaccat-3′) (ref. 37 and this work), and the products were cloned into the pTZ57R vector and sequenced. To assess the presence of the Leptomonas pyrrhocoris RNA virus, total cDNA of all Leptomonas pyrrhocoris isolates (both positive and negative as judged by gel-based assay) was amplified with primers LpTLV1F 5′-ttactcctataacggggca-3′ and LpTLV1R 5′-taaaggagcgaattctgct-3′ specific to the RDRP region (∼300 bp) of this virus and directly sequenced. Similarly, the occurrence of integrated virus in these isolates was checked by amplification using primers LpIVF 5′-cctatgcggatgcactcaa-3′ and LpIVR 5′-cttgtgcattttctatccaag-3′. PCR Primers M200 5′-atggctccvvtcaargtwggmat-3′ and M201 5′-takccccactcrttrtcrtacca-3′ for the glycosomal GAPDH (gGAPDH) gene were used as an internal positive control (104). Additional methods for cultivation of trypanosomatids, phylogenetic, genomic and transcriptomic analyses, and negative-stain transmission electron microscopy can be found in the Supporting Information.

Acknowledgments

We thank members of our groups and D. Wang and S. Krishnamurthy (Washington University School of Medicine in St. Louis) for discussion and D. Maslov, L. Simpson, and V. Bellofatto for providing Phytomonas and Leptomonas strains. This work was supported by NIH Grants R01-AI029646 and R56-AI099364 (to S.M.B.), the Intramural Research Program of the National Library of Medicine at the NIH (I.B.R.), Grant Agency of Czech Republic Awards 17-10656S (to V.Y. and J.V.) and 16-18699S (to J.L. and V.Y.), Moravskoslezský kraj Research Initiative DT01-021358 (to V.Y. and A.Y.K.), the Russian Foundation for Basic Research Project no. 15-29-02734 and the State Assignment for the Zoological Institute AAAA-A17-117030310322-3 (to A.O.F.), and the European Cooperation in Science and Technology (COST) action CM1307 (to J.L.). Work in the V.Y. laboratory is financially supported by the Ministry of Education, Youth and Sports of the Czech Republic in the National Feasibility Program I, project LO1208 TEWEP. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Footnotes

  • ↵1N.S.A. and A.Y.K. contributed equally to this work.

  • ↵2S.M.B. and V.Y. contributed equally to this work.

  • ↵3To whom correspondence may be addressed. Email: stephen.beverley{at}wustl.edu or vyacheslav.yurchenko{at}osu.cz.
  • Author contributions: D.G., N.S.A., A.Y.K., N.F., P.P., S.S., S.M.B., and V.Y. designed research; D.G., N.S.A., A.Y.K., A.K., L.-F.L., D.E.D., H.Z., P.K., and J.M. performed research; A.O.F., J.V., C.M.d.-L., and J.L. contributed new reagents/analytic tools; D.G., N.S.A., A.Y.K., A.B., I.B.R., S.M.B., and V.Y. analyzed data; and D.G., A.Y.K., S.M.B., and V.Y. wrote the paper.

  • Reviewers: S.A.G., University of Kentucky; M.L.N., Harvard Medical School; and L.S., University of California, Los Angeles.

  • The authors declare no conflict of interest.

  • Data deposition: The deposition and accession numbers for the sequences reported in this work are summarized in Table S3. Metatranscriptomic contigs assembled from single-read archive depositions, which cannot be deposited in GenBank because they have not been experimentally verified, are listed in Table S4.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1717806115/-/DCSupplemental.

Published under the PNAS license.

References

  1. ↵
    1. Shi M, et al.
    (2016) Redefining the invertebrate RNA virosphere. Nature 540:539–543.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Wang AL,
    2. Wang CC
    (1991) Viruses of the protozoa. Annu Rev Microbiol 45:251–263.
    OpenUrlCrossRefPubMed
  3. ↵
    1. Banik GR,
    2. Stark D,
    3. Rashid H,
    4. Ellis JT
    (2014) Recent advances in molecular biology of parasitic viruses. Infect Disord Drug Targets 14:155–167.
    OpenUrl
  4. ↵
    1. Ghabrial SA,
    2. Castón JR,
    3. Jiang D,
    4. Nibert ML,
    5. Suzuki N
    (2015) 50-plus years of fungal viruses. Virology 479-480:356–368.
    OpenUrl
  5. ↵
    1. Chen B,
    2. Geletka LM,
    3. Nuss DL
    (2000) Using chimeric hypoviruses to fine-tune the interaction between a pathogenic fungus and its plant host. J Virol 74:7562–7567.
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Ives A, et al.
    (2011) Leishmania RNA virus controls the severity of mucocutaneous leishmaniasis. Science 331:775–778.
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Fichorova RN, et al.
    (2012) Endobiont viruses sensed by the human host–Beyond conventional antiparasitic therapy. PLoS One 7:e48418.
    OpenUrlCrossRefPubMed
  8. ↵
    1. Opperdoes FR,
    2. Butenko A,
    3. Flegontov P,
    4. Yurchenko V,
    5. Lukeš J
    (2016) Comparative metabolism of free-living Bodo saltans and parasitic trypanosomatids. J Eukaryot Microbiol 63:657–678.
    OpenUrlCrossRef
  9. ↵
    1. Stuart K,
    2. Panigrahi AK
    (2002) RNA editing: Complexity and complications. Mol Microbiol 45:591–596.
    OpenUrlCrossRefPubMed
  10. ↵
    1. Záhonová K,
    2. Kostygov AY,
    3. Ševčíková T,
    4. Yurchenko V,
    5. Eliáš M
    (2016) An unprecedented non-canonical nuclear genetic code with all three termination codons reassigned as sense codons. Curr Biol 26:2364–2369.
    OpenUrlCrossRefPubMed
  11. ↵
    1. Simpson L
    (1998) A base called J. Proc Natl Acad Sci USA 95:2037–2038.
    OpenUrlFREE Full Text
  12. ↵
    1. Michaeli S
    (2011) Trans-splicing in trypanosomes: Machinery and its impact on the parasite transcriptome. Future Microbiol 6:459–474.
    OpenUrlCrossRefPubMed
  13. ↵
    1. Nussbaum K,
    2. Honek J,
    3. Cadmus CM,
    4. Efferth T
    (2010) Trypanosomatid parasites causing neglected diseases. Curr Med Chem 17:1594–1617.
    OpenUrlCrossRefPubMed
  14. ↵
    1. Maslov DA,
    2. Votýpka J,
    3. Yurchenko V,
    4. Lukeš J
    (2013) Diversity and phylogeny of insect trypanosomatids: All that is hidden shall be revealed. Trends Parasitol 29:43–52.
    OpenUrlCrossRefPubMed
  15. ↵
    1. Flegontov P, et al.
    (2013) Paratrypanosoma is a novel early-branching trypanosomatid. Curr Biol 23:1787–1793.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Lukeš J,
    2. Skalický T,
    3. Týč J,
    4. Votýpka J,
    5. Yurchenko V
    (2014) Evolution of parasitism in kinetoplastid flagellates. Mol Biochem Parasitol 195:115–122.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Kostygov AY,
    2. Yurchenko V
    (2017) Revised classification of the subfamily Leishmaniinae (Trypanosomatidae). Folia Parasitol 64:020.
    OpenUrlCrossRef
  18. ↵
    1. Molyneux DH
    (1974) Virus-like particles in Leishmania parasites. Nature 249:588–589.
    OpenUrlCrossRefPubMed
  19. ↵
    1. Croft SL,
    2. Chance ML,
    3. Gardener PJ
    (1980) Ultrastructural and biochemical characterization of stocks of Endotrypanum. Ann Trop Med Parasitol 74:585–589.
    OpenUrlPubMed
  20. ↵
    1. Soares MJ,
    2. Motta MC,
    3. de Souza W
    (1989) Bacterium-like endosymbiont and virus-like particles in the trypanosomatid Crithidia desouzai. Microbios Lett 41:137–141.
    OpenUrl
  21. ↵
    1. Marche S,
    2. Roth C,
    3. Manohar SK,
    4. Dollet M,
    5. Baltz T
    (1993) RNA virus-like particles in pathogenic plant trypanosomatids. Mol Biochem Parasitol 57:261–267.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Yurchenko V, et al.
    (2014) Ultrastructure and molecular phylogeny of four new species of monoxenous trypanosomatids from flies (Diptera: Brachycera) with redefinition of the genus Wallaceina. Folia Parasitol 61:97–112.
    OpenUrlCrossRefPubMed
  23. ↵
    1. Kraeva N, et al.
    (2015) Leptomonas seymouri: Adaptations to the dixenous life cycle analyzed by genome sequencing, transcriptome profiling and co-infection with Leishmania donovani. PLoS Pathog 11:e1005127.
    OpenUrlCrossRef
  24. ↵
    1. Widmer G,
    2. Comeau AM,
    3. Furlong DB,
    4. Wirth DF,
    5. Patterson JL
    (1989) Characterization of a RNA virus from the parasite Leishmania. Proc Natl Acad Sci USA 86:5979–5982.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Tarr PI, et al.
    (1988) LR1: A candidate RNA virus of Leishmania. Proc Natl Acad Sci USA 85:9572–9575.
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Eren RO, et al.
    (2016) Mammalian innate immune response to a Leishmania-resident RNA Virus increases macrophage survival to promote parasite persistence. Cell Host Microbe 20:318–328.
    OpenUrlCrossRefPubMed
  27. ↵
    1. Brettmann EA, et al.
    (2016) Tilting the balance between RNA interference and replication eradicates Leishmania RNA virus 1 and mitigates the inflammatory response. Proc Natl Acad Sci USA 113:11998–12005.
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Zangger H, et al.
    (2014) Leishmania aethiopica field isolates bearing an endosymbiontic dsRNA virus induce pro-inflammatory cytokine response. PLoS Negl Trop Dis 8:e2836.
    OpenUrlCrossRefPubMed
  29. ↵
    1. Rossi M, et al.
    (2017) Type I interferons induced by endogenous or exogenous viral infections promote metastasis and relapse of leishmaniasis. Proc Natl Acad Sci USA 114:4987–4992.
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Hartley MA,
    2. Drexler S,
    3. Ronet C,
    4. Beverley SM,
    5. Fasel N
    (2014) The immunological, environmental, and phylogenetic perpetrators of metastatic leishmaniasis. Trends Parasitol 30:412–422.
    OpenUrlCrossRefPubMed
  31. ↵
    1. Ronet C,
    2. Beverley SM,
    3. Fasel N
    (2011) Muco-cutaneous leishmaniasis in the New World: The ultimate subversion. Virulence 2:547–552.
    OpenUrlCrossRefPubMed
  32. ↵
    1. Adaui V, et al.
    (2016) Association of the endobiont double-stranded RNA virus LRV1 with treatment failure for human leishmaniasis caused by Leishmania braziliensis in Peru and Bolivia. J Infect Dis 213:112–121.
    OpenUrlCrossRefPubMed
  33. ↵
    1. Bourreau E, et al.
    (2016) Presence of Leishmania RNA Virus 1 in Leishmania guyanensis increases the risk of first-line treatment failure and symptomatic relapse. J Infect Dis 213:105–111.
    OpenUrlCrossRefPubMed
  34. ↵
    1. Ito MM, et al.
    (2015) Correlation between presence of Leishmania RNA virus 1 and clinical characteristics of nasal mucosal leishmaniosis. Rev Bras Otorrinolaringol (Engl Ed) 81:533–540.
    OpenUrl
  35. ↵
    1. Cantanhêde LM, et al.
    (2015) Further evidence of an association between the presence of Leishmania RNA Virus 1 and the mucosal manifestations in tegumentary leishmaniasis patients. PLoS Negl Trop Dis 9:e0004079.
    OpenUrlCrossRefPubMed
  36. ↵
    1. Pereira LdeO, et al.
    (2013) Severity of tegumentary leishmaniasis is not exclusively associated with Leishmania RNA virus 1 infection in Brazil. Mem Inst Oswaldo Cruz 108:665–667.
    OpenUrlCrossRefPubMed
  37. ↵
    1. Akopyants NS,
    2. Lye LF,
    3. Dobson DE,
    4. Lukeš J,
    5. Beverley SM
    (2016) A novel Bunyavirus-like virus of trypanosomatid protist parasites. Genome Announc 4:e00715–e00716.
    OpenUrlCrossRef
  38. ↵
    1. Lye LF,
    2. Akopyants NS,
    3. Dobson DE,
    4. Beverley SM
    (2016) A Narnavirus-like element from the trypanosomatid protozoan parasite Leptomonas seymouri. Genome Announc 4:e00713–e00716.
    OpenUrlCrossRef
  39. ↵
    1. Akopyants NS,
    2. Lye LF,
    3. Dobson DE,
    4. Lukeš J,
    5. Beverley SM
    (2016) A Narnavirus in the trypanosomatid protist plant pathogen Phytomonas serpens. Genome Announc 4:e00711–e00716.
    OpenUrlCrossRef
  40. ↵
    1. Sukla S,
    2. Roy S,
    3. Sundar S,
    4. Biswas S
    (2017) Leptomonas seymouri narna-like virus 1 and not leishmaniaviruses detected in kala-azar samples from India. Arch Virol 162:3827–3835.
    OpenUrlCrossRef
  41. ↵
    1. Beiting DP, et al.
    (2014) Differential induction of TLR3-dependent innate immune signaling by closely related parasite species. PLoS One 9:e88398.
    OpenUrlCrossRefPubMed
  42. ↵
    1. Votýpka J, et al.
    (2012) Cosmopolitan distribution of a trypanosomatid Leptomonas pyrrhocoris. Protist 163:616–631.
    OpenUrlCrossRefPubMed
  43. ↵
    1. Jacks T,
    2. Madhani HD,
    3. Masiarz FR,
    4. Varmus HE
    (1988) Signals for ribosomal frameshifting in the Rous sarcoma virus gag-pol region. Cell 55:447–458.
    OpenUrlCrossRefPubMed
  44. ↵
    1. Firth AE,
    2. Chung BY,
    3. Fleeton MN,
    4. Atkins JF
    (2008) Discovery of frameshifting in Alphavirus 6K resolves a 20-year enigma. Virol J 5:108.
    OpenUrlCrossRefPubMed
  45. ↵
    1. Atkins JF,
    2. Loughran G,
    3. Bhatt PR,
    4. Firth AE,
    5. Baranov PV
    (2016) Ribosomal frameshifting and transcriptional slippage: From genetic steganography and cryptography to adventitious use. Nucleic Acids Res 44:7007–7078.
    OpenUrlCrossRefPubMed
  46. ↵
    1. Dreher TW,
    2. Miller WA
    (2006) Translational control in positive strand RNA plant viruses. Virology 344:185–197.
    OpenUrlCrossRefPubMed
  47. ↵
    1. Kim KH,
    2. Lommel SA
    (1998) Sequence element required for efficient -1 ribosomal frameshifting in red clover necrotic mosaic dianthovirus. Virology 250:50–59.
    OpenUrlCrossRefPubMed
  48. ↵
    1. Olivier V, et al.
    (2008) Molecular characterisation and phylogenetic analysis of Chronic bee paralysis virus, a honey bee virus. Virus Res 132:59–68.
    OpenUrlCrossRefPubMed
  49. ↵
    1. Johnson KN,
    2. Johnson KL,
    3. Dasgupta R,
    4. Gratsch T,
    5. Ball LA
    (2001) Comparisons among the larger genome segments of six nodaviruses and their encoded RNA replicases. J Gen Virol 82:1855–1866.
    OpenUrlPubMed
  50. ↵
    1. Na H,
    2. White KA
    (2006) Structure and prevalence of replication silencer-3′ terminus RNA interactions in Tombusviridae. Virology 345:305–316.
    OpenUrlCrossRefPubMed
  51. ↵
    1. Fauquet CM,
    2. Mayo MA,
    3. Maniloff J,
    4. Desselberger U,
    5. Ball LA
    1. Schneemann A,
    2. Ball LA,
    3. Delsert C,
    4. Johnson JE,
    5. Nishizawa T
    (2005) Family Nodaviridae. Virus Taxonomy: Eighth Report of the International Committee on Taxonomy of Viruses, eds Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA (Elsevier, New York), pp 865–872.
  52. ↵
    1. Socha R
    (1993) Pyrrhocoris apterus (Heteroptera) - an experimental model species: A review. Eur J Entomol 90:241–286.
    OpenUrl
  53. ↵
    1. Flegontov P, et al.
    (2016) Genome of Leptomonas pyrrhocoris: A high-quality reference for monoxenous trypanosomatids and new insights into evolution of Leishmania. Sci Rep 6:23704.
    OpenUrl
  54. ↵
    1. Feschotte C,
    2. Gilbert C
    (2012) Endogenous viruses: Insights into viral evolution and impact on host biology. Nat Rev Genet 13:283–296.
    OpenUrlCrossRefPubMed
  55. ↵
    1. Bruenn JA,
    2. Warner BE,
    3. Yerramsetty P
    (2015) Widespread mitovirus sequences in plant genomes. PeerJ 3:e876.
    OpenUrlCrossRefPubMed
  56. ↵
    1. Geuking MB, et al.
    (2009) Recombination of retrotransposon and exogenous RNA virus results in nonretroviral cDNA integration. Science 323:393–396.
    OpenUrlAbstract/FREE Full Text
  57. ↵
    1. Viruses ICoTo
    (2012) Virus Taxonomy: Classification and Nomenclature of Viruses: Ninth Report of the International Committee on Taxonomy of Viruses (Academic, London), p 1327.
  58. ↵
    1. Kelley LA,
    2. Mezulis S,
    3. Yates CM,
    4. Wass MN,
    5. Sternberg MJ
    (2015) The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc 10:845–858.
    OpenUrlCrossRefPubMed
  59. ↵
    1. Soding J,
    2. Biegert A,
    3. Lupas AN
    (2005) The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res 33:W244–W248.
    OpenUrlCrossRefPubMed
  60. ↵
    1. Routhier E,
    2. Bruenn JA
    (1998) Functions of conserved motifs in the RNA-dependent RNA polymerase of a yeast double-stranded RNA virus. J Virol 72:4427–4429.
    OpenUrlAbstract/FREE Full Text
  61. ↵
    1. te Velthuis AJW
    (2014) Common and unique features of viral RNA-dependent polymerases. Cell Mol Life Sci 71:4403–4420.
    OpenUrlCrossRefPubMed
  62. ↵
    1. Gong P,
    2. Peersen OB
    (2010) Structural basis for active site closure by the poliovirus RNA-dependent RNA polymerase. Proc Natl Acad Sci USA 107:22505–22510.
    OpenUrlAbstract/FREE Full Text
  63. ↵
    1. Bishop DH
    (1986) Ambisense RNA viruses: Positive and negative polarities combined in RNA virus genomes. Microbiol Sci 3:183–187.
    OpenUrlPubMed
  64. ↵
    1. Elliott RM
    (1990) Molecular biology of the Bunyaviridae. J Gen Virol 71:501–522.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Reguera J,
    2. Weber F,
    3. Cusack S
    (2010) Bunyaviridae RNA polymerases (L-protein) have an N-terminal, influenza-like endonuclease domain, essential for viral cap-dependent transcription. PLoS Pathog 6:e1001101.
    OpenUrlCrossRefPubMed
  66. ↵
    1. Sanchez AJ,
    2. Vincent MJ,
    3. Nichol ST
    (2002) Characterization of the glycoproteins of Crimean-Congo hemorrhagic fever virus. J Virol 76:7263–7275.
    OpenUrlAbstract/FREE Full Text
  67. ↵
    1. Cifuentes-Muñoz N,
    2. Salazar-Quiroz N,
    3. Tischler ND
    (2014) Hantavirus Gn and Gc envelope glycoproteins: Key structural units for virus cell entry and virus assembly. Viruses 6:1801–1822.
    OpenUrlCrossRefPubMed
  68. ↵
    1. Moriconi M, et al.
    (2017) Phlebotomine sand fly-borne pathogens in the Mediterranean Basin: Human leishmaniasis and phlebovirus infections. PLoS Negl Trop Dis 11:e0005660.
    OpenUrl
  69. ↵
    1. Li CX, et al.
    (2015) Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. eLife 4:e05378.
    OpenUrlCrossRef
  70. ↵
    1. Remnant EJ, et al.
    (2017) A diverse range of novel RNA viruses in geographically distinct honey bee populations. J Virol 91:e00158–17.
    OpenUrl
  71. ↵
    1. Dudas G,
    2. Obbard DJ
    (2015) Are arthropods at the heart of virus evolution? eLife 4:e06837.
    OpenUrl
  72. ↵
    1. Marklewitz M,
    2. Zirkel F,
    3. Kurth A,
    4. Drosten C,
    5. Junglen S
    (2015) Evolutionary and phenotypic analysis of live virus isolates suggests arthropod origin of a pathogenic RNA virus family. Proc Natl Acad Sci USA 112:7536–7541.
    OpenUrlAbstract/FREE Full Text
  73. ↵
    1. Franzosa EA, et al.
    (2015) Sequencing and beyond: Integrating molecular ‘omics’ for microbial community profiling. Nat Rev Microbiol 13:360–372.
    OpenUrlCrossRefPubMed
  74. ↵
    1. d’Avila-Levy CM, et al.
    (2015) Exploring the environmental diversity of kinetoplastid flagellates in the high-throughput DNA sequencing era. Mem Inst Oswaldo Cruz 110:956–965.
    OpenUrlCrossRefPubMed
  75. ↵
    1. Votýpka J, et al.
    (2014) Kentomonas gen. n., a new genus of endosymbiont-containing trypanosomatids of Strigomonadinae subfam. n. Protist 165:825–838.
    OpenUrlCrossRefPubMed
  76. ↵
    1. Kostygov AY,
    2. Grybchuk-Ieremenko A,
    3. Malysheva MN,
    4. Frolov AO,
    5. Yurchenko V
    (2014) Molecular revision of the genus Wallaceina. Protist 165:594–604.
    OpenUrlCrossRefPubMed
  77. ↵
    1. Yurchenko VY,
    2. Lukeš J,
    3. Jirků M,
    4. Maslov DA
    (2009) Selective recovery of the cultivation-prone components from mixed trypanosomatid infections: A case of several novel species isolated from Neotropical Heteroptera. Int J Syst Evol Microbiol 59:893–909.
    OpenUrlCrossRefPubMed
  78. ↵
    1. Votýpka J, et al.
    (2010) Probing into the diversity of trypanosomatid flagellates parasitizing insect hosts in South-West China reveals both endemism and global dispersal. Mol Phylogenet Evol 54:243–253.
    OpenUrlCrossRefPubMed
  79. ↵
    1. Rodríguez-Cousiño N,
    2. Solórzano A,
    3. Fujimura T,
    4. Esteban R
    (1998) Yeast positive-stranded virus-like RNA replicons. 20 S and 23 S RNA terminal nucleotide sequences and 3′ end secondary structures resemble those of RNA coliphages. J Biol Chem 273:20363–20371.
    OpenUrlAbstract/FREE Full Text
  80. ↵
    1. Murphy WJ,
    2. Watkins KP,
    3. Agabian N
    (1986) Identification of a novel Y branch structure as an intermediate in trypanosome mRNA processing: Evidence for trans splicing. Cell 47:517–525.
    OpenUrlCrossRefPubMed
  81. ↵
    1. Günzl A
    (2010) The pre-mRNA splicing machinery of trypanosomes: Complex or simplified? Eukaryot Cell 9:1159–1170.
    OpenUrlAbstract/FREE Full Text
  82. ↵
    1. Mandelboim M,
    2. Estraño CL,
    3. Tschudi C,
    4. Ullu E,
    5. Michaeli S
    (2002) On the role of exon and intron sequences in trans-splicing utilization and cap 4 modification of the trypanosomatid Leptomonas collosoma SL RNA. J Biol Chem 277:35210–35218.
    OpenUrlAbstract/FREE Full Text
  83. ↵
    1. Cai G,
    2. Myers K,
    3. Fry WE,
    4. Hillman BI
    (2012) A member of the virus family Narnaviridae from the plant pathogenic oomycete Phytophthora infestans. Arch Virol 157:165–169.
    OpenUrlCrossRefPubMed
  84. ↵
    1. Reinhardt JA, et al.
    (2014) Meiotic drive impacts expression and evolution of x-linked genes in stalk-eyed flies. PLoS Genet 10:e1004362.
    OpenUrlCrossRefPubMed
  85. ↵
    1. Rastgou M, et al.
    (2009) Molecular characterization of the plant virus genus Ourmiavirus and evidence of inter-kingdom reassortment of viral genome segments as its possible route of origin. J Gen Virol 90:2525–2535.
    OpenUrlCrossRefPubMed
  86. ↵
    1. Hillman BI,
    2. Cai G
    (2013) The family narnaviridae: Simplest of RNA viruses. Adv Virus Res 86:149–176.
    OpenUrlCrossRefPubMed
  87. ↵
    1. Koonin EV,
    2. Dolja VV,
    3. Krupovic M
    (2015) Origins and evolution of viruses of eukaryotes: The ultimate modularity. Virology 479–480:2–25.
    OpenUrl
  88. ↵
    1. Dolja VV,
    2. Krupovic M
    (2013) Accelerating expansion of the viral universe. Curr Opin Virol 3:542–545.
    OpenUrl
  89. ↵
    1. Araújo JP,
    2. Hughes DP
    (2016) Diversity of entomopathogenic fungi: Which groups conquered the insect body? Adv Genet 94:1–39.
    OpenUrlCrossRef
  90. ↵
    1. Mitchell PL
    (2004) Heteroptera as vectors of plant pathogens. Neotrop Entomol 33:519–545.
    OpenUrlCrossRef
  91. ↵
    1. Ro YT,
    2. Scheffter SM,
    3. Patterson JL
    (1997) Hygromycin B resistance mediates elimination of Leishmania virus from persistently infected parasites. J Virol 71:8991–8998.
    OpenUrlAbstract/FREE Full Text
  92. ↵
    1. Patterson JL
    (1993) The current status of Leishmania RNA virus 1. Parasitol Today 9:135–136.
    OpenUrlPubMed
  93. ↵
    1. Widmer G,
    2. Dooley S
    (1995) Phylogenetic analysis of Leishmania RNA virus and Leishmania suggests ancient virus-parasite association. Nucleic Acids Res 23:2300–2304.
    OpenUrlCrossRefPubMed
  94. ↵
    1. Hartley MA,
    2. Ronet C,
    3. Zangger H,
    4. Beverley SM,
    5. Fasel N
    (2012) Leishmania RNA virus: When the host pays the toll. Front Cell Infect Microbiol 2:99.
    OpenUrlPubMed
  95. ↵
    1. Obbard DJ,
    2. Gordon KH,
    3. Buck AH,
    4. Jiggins FM
    (2009) The evolution of RNAi as a defence against viruses and transposable elements. Philos Trans R Soc Lond B Biol Sci 364:99–115.
    OpenUrlCrossRefPubMed
  96. ↵
    1. Matveyev AV, et al.
    (2017) The evolutionary loss of RNAi key determinants in kinetoplastids as a multiple sporadic phenomenon. J Mol Evol 84:104–115.
    OpenUrl
  97. ↵
    1. Lye LF, et al.
    (2010) Retention and loss of RNA interference pathways in trypanosomatid protozoans. PLoS Pathog 6:e1001161.
    OpenUrlCrossRefPubMed
  98. ↵
    1. Ghosh S,
    2. Banerjee P,
    3. Sarkar A,
    4. Datta S,
    5. Chatterjee M
    (2012) Coinfection of Leptomonas seymouri and Leishmania donovani in Indian leishmaniasis. J Clin Microbiol 50:2774–2778.
    OpenUrlAbstract/FREE Full Text
  99. ↵
    1. Krautz R,
    2. Arefin B,
    3. Theopold U
    (2014) Damage signals in the insect immune response. Front Plant Sci 5:342.
    OpenUrlCrossRefPubMed
  100. ↵
    1. Podlipaev SA
    (2000) Insect trypanosomatids: The need to know more. Mem Inst Oswaldo Cruz 95:517–522.
    OpenUrlCrossRefPubMed
  101. ↵
    1. Chomczyński P,
    2. Sacchi N
    (1987) Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 162:156–159.
    OpenUrlCrossRefPubMed
  102. ↵
    1. Potgieter AC, et al.
    (2009) Improved strategies for sequence-independent amplification and sequencing of viral double-stranded RNA genomes. J Gen Virol 90:1423–1432.
    OpenUrlCrossRefPubMed
  103. ↵
    1. Liu H, et al.
    (2012) Evolutionary genomics of mycovirus-related dsRNA viruses reveals cross-family horizontal gene transfer and evolution of diverse viral lineages. BMC Evol Biol 12:91.
    OpenUrlCrossRefPubMed
  104. ↵
    1. Maslov DA,
    2. Yurchenko VY,
    3. Jirků M,
    4. Lukeš J
    (2010) Two new species of trypanosomatid parasites isolated from Heteroptera in Costa Rica. J Eukaryot Microbiol 57:177–188.
    OpenUrlCrossRefPubMed
    1. Yurchenko V,
    2. Lukeš J,
    3. Tesarová M,
    4. Jirků M,
    5. Maslov DA
    (2008) Morphological discordance of the new trypanosomatid species phylogenetically associated with the genus Crithidia. Protist 159:99–114.
    OpenUrlCrossRefPubMed
    1. Katoh K,
    2. Standley DM
    (2013) MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol Biol Evol 30:772–780.
    OpenUrlCrossRefPubMed
    1. Capella-Gutiérrez S,
    2. Silla-Martínez JM,
    3. Gabaldón T
    (2009) trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973.
    OpenUrlCrossRefPubMed
    1. Nguyen LT,
    2. Schmidt HA,
    3. von Haeseler A,
    4. Minh BQ
    (2015) IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274.
    OpenUrlCrossRefPubMed
    1. Ronquist F, et al.
    (2012) MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol 61:539–542.
    OpenUrlCrossRefPubMed
    1. Chistyakova LV,
    2. Kostygov AY,
    3. Kornilova OA,
    4. Yurchenko V
    (2014) Reisolation and redescription of Balantidium duodeni Stein, 1867 (Litostomatea, Trichostomatia). Parasitol Res 113:4207–4215.
    OpenUrlCrossRefPubMed
    1. Kalyaanamoorthy S,
    2. Minh BQ,
    3. Wong TKF,
    4. von Haeseler A,
    5. Jermiin LS
    (2017) ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat Methods 14:587–589.
    OpenUrlCrossRefPubMed
    1. Huang X,
    2. Madan A
    (1999) CAP3: A DNA sequence assembly program. Genome Res 9:868–877.
    OpenUrlAbstract/FREE Full Text
    1. Mullapudi E, et al.
    (2016) Structure and genome release mechanism of the human cardiovirus Saffold virus 3. J Virol 90:7628–7639.
    OpenUrlAbstract/FREE Full Text
    1. Owczarzy R, et al.
    (2008) IDT SciTools: A suite for analysis and design of nucleic acid oligomers. Nucleic Acids Res 36:W163–W169.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Viral discovery and diversity in trypanosomatid protozoa with a focus on relatives of the human parasite Leishmania
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Viral discovery in trypanosomatids
Danyil Grybchuk, Natalia S. Akopyants, Alexei Y. Kostygov, Aleksandras Konovalovas, Lon-Fye Lye, Deborah E. Dobson, Haroun Zangger, Nicolas Fasel, Anzhelika Butenko, Alexander O. Frolov, Jan Votýpka, Claudia M. d’Avila-Levy, Pavel Kulich, Jana Moravcová, Pavel Plevka, Igor B. Rogozin, Saulius Serva, Julius Lukeš, Stephen M. Beverley, Vyacheslav Yurchenko
Proceedings of the National Academy of Sciences Jan 2018, 115 (3) E506-E515; DOI: 10.1073/pnas.1717806115

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Viral discovery in trypanosomatids
Danyil Grybchuk, Natalia S. Akopyants, Alexei Y. Kostygov, Aleksandras Konovalovas, Lon-Fye Lye, Deborah E. Dobson, Haroun Zangger, Nicolas Fasel, Anzhelika Butenko, Alexander O. Frolov, Jan Votýpka, Claudia M. d’Avila-Levy, Pavel Kulich, Jana Moravcová, Pavel Plevka, Igor B. Rogozin, Saulius Serva, Julius Lukeš, Stephen M. Beverley, Vyacheslav Yurchenko
Proceedings of the National Academy of Sciences Jan 2018, 115 (3) E506-E515; DOI: 10.1073/pnas.1717806115
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Microbiology
Proceedings of the National Academy of Sciences: 115 (3)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results and Discussion
    • Conclusions
    • Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Water from a faucet fills a glass.
News Feature: How “forever chemicals” might impair the immune system
Researchers are exploring whether these ubiquitous fluorinated molecules might worsen infections or hamper vaccine effectiveness.
Image credit: Shutterstock/Dmitry Naumov.
Reflection of clouds in the still waters of Mono Lake in California.
Inner Workings: Making headway with the mysteries of life’s origins
Recent experiments and simulations are starting to answer some fundamental questions about how life came to be.
Image credit: Shutterstock/Radoslaw Lecyk.
Cave in coastal Kenya with tree growing in the middle.
Journal Club: Small, sharp blades mark shift from Middle to Later Stone Age in coastal Kenya
Archaeologists have long tried to define the transition between the two time periods.
Image credit: Ceri Shipton.
Mouse fibroblast cells. Electron bifurcation reactions keep mammalian cells alive.
Exploring electron bifurcation
Jonathon Yuly, David Beratan, and Peng Zhang investigate how electron bifurcation reactions work.
Listen
Past PodcastsSubscribe
Panda bear hanging in a tree
How horse manure helps giant pandas tolerate cold
A study finds that giant pandas roll in horse manure to increase their cold tolerance.
Image credit: Fuwen Wei.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Cozzarelli Prize
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490