Skip to main content
  • Brief Report
  • Open access
  • Published:

Two +ssRNA mycoviruses cohabiting the fungal cultivar of leafcutter ants

Abstract

Leafcutter ants are dominant herbivores in the Neotropics and rely on a fungus (Leucoagaricus gongylophorus) to transform freshly gathered leaves into a source of nourishment rather than consuming the vegetation directly. Here we report two virus-like particles that were isolated from L. gongylophorus and observed using transmission electron microscopy. RNA sequencing identified two +ssRNA mycovirus strains, Leucoagaricus gongylophorus tymo-like virus 1 (LgTlV1) and Leucoagaricus gongylophorus magoulivirus 1 (LgMV1). Genome annotation of LgTlV1 (7401 nt) showed conserved domains for methyltransferase, endopeptidase, viral RNA helicase, and RNA-dependent RNA polymerase (RdRp). The smaller genome of LgMV1 (2636 nt) contains one open reading frame encoding an RdRp. While we hypothesize these mycoviruses function as symbionts in leafcutter farming systems, further study will be needed to test whether they are mutualists, commensals, or parasites.

Viruses infecting fungi (mycoviruses) appear to be common across the fungal kingdom [1]. Mycoviruses are obligately dependent on their fungal host since their inability to lyse or live outside of host cells precludes dispersal [1, 2]. Mycoviruses can be parasites that reduce host fitness, commensals that cause asymptomatic infections [3, 4], or mutualists that promote adaptive phenotypes in their hosts [5, 6]. Mycoviruses are thus likely a frontier in our understanding of fungal ecology and evolution.

Leafcutter ants depend on a fungus (Leucoagaricus gongylophorus) that they cultivate by provisioning with freshly cut plant fragments [7, 8]. In return, the fungus produces food reward structures for the ants, in the form of swollen hyphal cells called gongylidia. Each gongylidium contains a large vacuole that originates from an autophagic process and stores metabolites that go beyond nutritionally sustaining the ant farmers, to enzymatically breaking down freshly deposited plant material [9,10,11]. When leafcutter ants ingest gongylidia, they do not digest these enzymes but vector them in their fecal fluid onto the fungus garden to assist with the decomposition of newly deposited plant fragments [12, 13].

In the present study, we hypothesized that mycoviruses are common in leafcutter farming systems, given their ubiquitousness across the fungal kingdom, and could potentially have important functions. Other secondary microbial symbionts (e.g., bacteria) are present in the leafcutter farming system and promote disease resistance [14,15,16] or provide crucial nutritional services in ant guts [17] and in the fungus garden matrix [18]. In contrast, little is known about potential roles for mycoviruses.

Two seldom cited studies published nearly 40 years ago used transmission electron microscopy (TEM) to detect what were provisionally considered to be virus-like particles (VLPs) in L. gongylophorus isolated from a colony of the leafcutter ant Acromyrmex octospinosus [19, 20]. More recently, an in silico study using publicly available fungal transcriptomes identified potential viral sequences from the fungal cultivar isolated from an Acromyrmex echinatior colony (Ae332) [21, 22]. Here, we sought to first verify the presence of VLPs using TEM and then characterize the resident fungal virome using total RNA sequencing. Additionally, we sought to investigate whether leafcutter ants could potentially vector VLPs from ingested gongylidium cells back to the fungal cultivar in their fecal fluid. This could indicate that the mycoviruses can exist in an extracellular state and suggest further symbiotic functions in fungus farming systems as beneficial mutualists, costly parasites, or neutral commensals.

We first sought to visualize VLPs in axenic cultures of L. gongylophorus isolated from an Atta colombica colony (Ac2012-1) collected in a Panamanian rainforest (Table S1). We used a modified affinity chromatography method that was initially designed for isolating enveloped viruses by binding viral surface glycoconjugates to a lectin (concanavalin A), but that is also likely effective in binding other viral components in non-enveloped viruses [23, 24]. Briefly, samples were individually freeze-dried, ground in liquid nitrogen and then resuspended in Concanavalin A Sepharose 4B (Cytiva, MA, USA) buffer. We next filtered the samples, isolated VLPs with a modified version of the manufacturer’s instructions and prepared a negative control using only buffer (Supp. File S1). In addition, fecal fluid was collected from three worker ants of A. colombica (Ac2012-1) by pressing their abdomens with forceps and mixing the exuded liquid with 200 µL ddH2O. Isolated VLPs were prepared for TEM visualization on glow-discharged 200 mesh carbon-coated grids using 2% phosphotungstic acid (pH 7) as a negative stain. These samples were then examined on a CM100 BioTWIN (Phillips, Eindhoven, The Netherlands) transmission electron microscope at an accelerating voltage of 80 kV.

TEM images revealed two VLP phenotypes found in fungal isolate samples and in the fecal fluid of worker ants: (1) an icosahedral VLP of 30-nm in diameter (Fig. 1; VLP1), and (2) a bacilliform VLP 15-nm wide with lengths of ca. 30–150 nm (Fig. 1; VLP2). No VLPs were found in the negative control.

Fig. 1
figure 1

TEM images of negative stained virus-like particles (VLPs) extracted from L. gongylophorus fungal isolates farmed by A. colombica (Ac2012-1) and the fecal fluid of the worker ants from the same colony. VLP1 presents a 30-nm wide icosahedral geometry with a high similarity to known tymovirus virions [25]. VLP2 has a 15-nm wide bacilliform morphology and varies between 30 and 100 nm long like known virions of ourmiaviruses in the family Botourmiaviridae [26, 27]. The scale bar for all panels is 50 nm

We complemented this qualitative visualization with molecular and computational methods to identify the VLPs. We extracted total RNA from in vitro fungus cultures from the same A. colombica colony (Ac2012-1) with a RNeasy® Plant Mini Kit (Qiagen) using buffer RLC supplemented with DTT, following the manufacturer's protocol. Total RNA was processed for rRNA depletion, and next by 2 × 100 bp library construction using stranded cDNA (dUTP) synthesis. Library preparation and sequencing (DNBSEQ platform) were performed by BGI (Hong Kong). Low-quality reads were filtered out, and trimmed reads were provided for downstream analysis. The short reads were mapped to a high-quality draft genome from the same isolate of L. gongylophorus (Ac2012-1) [28] using bwa-mem2 v.2.2.1 [29]. A total of 1,323,206 unmapped read pairs were extracted with Samtools v.1.16 [30], and used for de novo RNA assembly with Trinity v.2.12.0 [31]. This generated 798 contiguous sequences.

We next used cd-hit v.4.8.1 [32] to cluster sequences with 95% similarity and filtered the dataset to remove sequences shorter than 700 nucleotides and thus limit the search for near-complete viral genomes based on the smaller genome segments in the NCBI plant and fungi virus genome database. The final dataset comprised 48 sequences ranging from 702 to 7401 nucleotides (Supp. File S2). We then identified open reading frames (ORFs) longer than 100 amino acids and translated them into amino acid sequences using Transdecoder v.5.5.0 (Haas, BJ. https://github.com/TransDecoder/TransDecoder) with default parameters (Supp. File S3). These assembled nucleotide and amino acid sequences were assessed using the Basic Local Alignment Search Tool (BLAST) against the NCBI Viral RefSeq and a custom database including all viral genomes known from fungi and plants, using BLASTn, BLASTx, and BLASTp, with both NCBI BLAST v.2.12.0+ [33] and Diamond v.2.0.13 [34] (Supp. File S1). After manually curating these sequences, we found only two viral candidate sequences of 7401 and 2636 nucleotides matching mycoviruses sequences from Tymovirales and Botourmiaviridae, respectively. These sequences are deposited on NCBI under the accessions OR914225 and OR914226, respectively. We further attempted to obtain the 5' and 3' termini of the two candidate viral sequences using 5' and 3' RACE, however, we were unsuccessful, and the terminal sequences of both viruses remain undetermined. However, the main ORFs were completely spanned by the total RNA sequencing, which allowed us to identify the main components of the viral genomes (see below).

We further mapped all reads to the viral sequences and used BEDTools v2.30.0 [35] to calculate the read count and read depth for each of them. We found a read depth of 102.17 on average for the Tymovirales sequence with read count of 7,615 reads while the Botourmiaviridae sequence had a read depth of 688.3 on average with 18,194 mapping reads (Supp. Fig. S1). Given that we used a dUTP method to sequence strand specific RNA reads, we used the mapped reads to calculate sequence polarity bias for both viral sequences. We found strong bias for positive strand in both sequences, with only two reads (0.02%) mapping to the negative strand in the Tymovirales sequence and 389 (2.13%) in the Botourmiaviridae, suggesting a low replication rate for both viruses. We then used BCFtools v1.20 [30] to call variant sites and quantify the polymorphisms in each sequence. This analysis showed that the Tymovirales sequence has 496 (6.7%) polymorphic sites, while the Botourmiaviridae has only 13 (0.49%).

To further validate the presence of the viruses, we attempted to amplify a short region (150–200 nt) of their genomes with PCR followed by Sanger sequencing. However, we only succeeded in amplifying and sequencing the sequence related to Botourmiaviridae, although we designed and tested 18 different primer pairs in multiple regions of the sequence related to Tymovirales (see Supp. File S1). While we were not able to validate the presence of the virus related to Tymovirales in our sample via amplification, we can assume that it maybe is present in L. gongylophorus, given that it was also found by Jo et al. [21] in three other samples previously sequenced by De Fine Licht et al. [22]. The high level of polymorphic sites, combined with the relatively low abundance of these sequences as indicated by the read depth and the low replication rate, likely contributed to the failure of RACE and PCR reactions.

We further assessed by BLAST the two mycovirus sequences against NCBI non-redundant viral sequences. The BLASTn analysis of the 7401 nt sequence showed high similarity to three sequences previously identified by Jo et al. [21] as Leucoagaricus tymovirus A (n = 1) and Leucoagaricus tymovirus B (n = 2) found in an isolate of L. gongylophorus. These sequences showed at least 95% nucleotide identity with our sequence. However, none of these three tymo-like viral sequences spanned the complete extension of our sequence (Supp. File S1, Fig. S2). The BLASTp alignment of the longest ORF showed a similar result with at least 98.2% amino acid identity to the polyprotein encoded by the same Leucoagaricus tymovirus A and B that also did not span completely our sequence (Table S2). This suggests that the three sequences found by Jo et al. [21]: (1) belong to the same viral species that we observed, (2) were more fragmented than the sequence we recovered (see Supp. File S1, Figs. S2 and S3) and (3) might be a common mycovirus across farming systems of different leafcutter ant genera as it has been found in L. gongylophorus from two different colonies: Atta colombica (this study) and Acromyrmex echinatior (Ae322) [21].

In contrast, BLASTn analysis of the 2636 nt sequence had a single hit to Leucoagaricus ourmiavirus F also identified from the same fungal isolate as above by Jo et al. [21] with 84% nucleotide identity, but in a very short (7%) region of our sequence. The BLASTp analysis of the longest ORF showed the highest amino acid similarity to the RNA-dependent RNA polymerase (RdRp) of Lentinula edodes magoulivirus virus 1 with 43% amino acid identity (Table S2). While this low similarity to known sequences suggests our viral sequence is an unclassified mycovirus, it still appears phylogenetically aligned with another virus found in leafcutter ant farming system.

To assess these conjectures, we next performed phylogenetic analyses for both putative viral genomes to confirm their placement among known mycoviruses. We used the amino acid sequence of the RdRp predicted from each of our genomes combined with their best hits from BLAST results as well as sequences accepted by the International Committee on Taxonomy of Viruses; (ICTV) related to Tymovirales or Botourmiaviridae. Local pairwise alignment was carried out using MAFFT v7.5 [36] and was followed by manual trimming with Aliview v1.8 [37]. We used IQ-TREE v 1.6.12 [38] to choose the best evolutionary model with the built-in ModelFinder [39] and then build the maximum likelihood trees using 1000 ultrafast bootstraps (UFBoot) [40]. Alignments are available as supplementary file S4.

We first provided support for the hypothesis that the putative tymo-like virus we sequenced belonged to the order Tymovirales since it grouped in a clade with, among others, Leucoagaricus tymovirus A and Leucoagaricus tymovirus B (Fig. 2A). Considering these sequences are highly similar to each other (Tables S2 and S3) and the genomic features discussed below, and considering the clade formed by them do not contain tymovirus species accepted by the ICTV (Fig. 2), we propose to reassign the strain name to “Leucoagaricus gongylophorus tymo-like virus 1” (LgTlV1) [41]. We next found support that the hypothesized magoulivirus classification was correct as our sequence grouped with maximum support (UFBoot = 100) in a clade with Lentinula edodes magoulivirus 1 (Fig. 3A). The phylogenetic tree further supported close relationships between this clade and viruses in the family Botourmiaviridae. Based on this evidence and the genomic structure that we describe below, we propose that our novel mycovirus be named “Leucoagaricus gongylophorus magoulivirus 1” (LgMV1).

Fig. 2
figure 2

Phylogenetic placement and genomic structure of Leucoagaricus gongylophorus tymo-like virus 1 (LgTlV1) A maximum-likelihood phylogenetic tree based on RdRp alignment (1755 aa) showing the position of LgTlV1 (green) within the order Tymovirales. Support values on the branches are ultrafast bootstrap values calculated with 1000 replicates (only shown UFBoot ≥ 50). The scale bar indicates 0.5 amino acid substitutions per site. B Genome structure with ORF positions and conserved sequence regions predicted

Fig. 3
figure 3

Phylogenetic placement and genomic structure of Leucoagaricus gongylophorus magoulivirus 1 (LgMV1). A Maximum-likelihood phylogenetic tree based on RdRp alignment (1748 aa) showing the position of LgMV1 (green) within the family Botourmiaviridae. Support values on the branches are ultrafast bootstraps calculated with 1000 replicates (only shown UFBoot ≥ 50). The scale bar indicates 0.5 amino acid substitutions per site. B LgMV1 contains a single ORF that is highly similar to known RdRp sequences

We next annotated the viral genomes. Three main positive-sense ORFs were identified within the LgTlV1 genome (Fig. 2B). We annotated these ORFs using pHMMER 2.41.2 [42], InterProScan 5 v5.57–90.0 [43], and cdBlast v3.18 [44]. The largest ORF (ORF1) was located in the second reading frame and spanned 5652 nucleotides. It was predicted to encode a polyprotein containing domains with similarities to methyltransferase, endopeptidase, viral RNA helicase, and RdRp. The two other ORFs were located in the third reading frame. The smaller (ORF2) spanned 810 nucleotides and partially overlapped with ORF1. The third ORF (ORF3) spanned 885 nucleotides. No similarities to known proteins were found for ORF2 or ORF3. The annotation complements our BLAST and phylogenetic results, highlighting the consistency with other genomic studies of tymo-like viruses in terms of genome size and structure [25, 26]. These results further support that our sequence data provide a near-complete genome of Leucoagaricus gongylophorus tymo-like virus 1 (LgTlV1) relative to the ones that were previously identified, such as those of Leucoagaricus tymovirus A and Leucoagaricus tymovirus B [21], which are small and fragmented sequences of LgTlV1.

We did not identify a coat protein for LgTlV1. To investigate the possibility of a functional coat protein, we used AlphaFold v2.3.2 [45] to predict the protein structures originating from ORF2 and ORF3 and compare with the coat proteins of the related and well-characterized tymoviruses turnip yellow mosaic virus (TYMV) and belladonna mottle virus (BDMV) [25] (Fig. S4). While ORF2 and ORF3 exhibited several disordered segments, ORF3 displayed a structural central domain with the presence of alpha-helices and a beta-barrel within its protein structure (Fig. S4) reminiscent of the structure of coat proteins of TYMV and BDMV. Further investigation of this structure is needed to confirm the presence of a coat protein in LgTlV1 [46]. Using the RNAfold webserver [47], we further predicted a tRNA-like secondary structure in the final 140 nucleotides of the LgTlV1 genome (Fig. 2F). Since this feature is commonly found in plant tymoviruses [25], it reinforces their similarities and the completeness of our assembled genome sequence.

In the genome of LgMV1, we identified a positive-sense ORF spanning 1839 nucleotides. This ORF encodes a protein closely resembling a viral +ssRNA RdRp (Fig. 3B). LgMV1 resembles other species of Botourmiaviridae where all except ourmiaviruses are not encapsidated and have monopartite genomes encoding RdRp [26, 48]. While bacilliform virions similar to VLP2 (Fig. 1) are associated with ourmiaviruses that encode a coat protein in a separate genome segment, we did not detect such a genomic segment.

We thus cannot exclude the possibility that the observed VLP1 and VLP2 belong to LgTlV1 and LgMV1, respectively, or that they are other viruses we failed to detect (or perhaps not virus particles).

While the present study confirms the presence of two mycoviruses cohabiting L. gongylophorus, further study will be needed to test their roles as symbionts in the leafcutter ant system. We assume the mycoviruses we detected are not usually lethal pathogens, since their host leafcutter colony: (1) was likely infected since it was collected in Panama, and (2) survived in the laboratory for nearly a decade. It is of course possible that viral loads fluctuate over time and that costly parasitic functions become more pronounced if the host farming system becomes immunocompromised.

We hypothesize that, since they lack an obvious dispersal mechanism, the mycoviruses we detected are strongly dependent on the continuous health of the fungus garden, and are thus either commensals or mutualists. Detection of intact viruses in ant fecal droplets suggests that ants ingest them when they consume gongylidia and then actively vector them back to the fungus as they prepare the plant substrate. Viral vectoring would add to the list of other fungus-derived metabolites actively vectored in an undigested fully functional form back to the fungal cultivar [10, 12, 13]. Further study can explore whether the adaptive absence of active proteases in leafcutter ant digestive systems also promotes viral stability.

We propose two potential benefits of such viral stability. First, as viruses are known to trigger autophagy in their hosts [49], the identified mycoviruses may trigger autophagic processes associated with gongylidia formation [11]. Second, mycoviruses may cause the fungus to produce secondary metabolites that in turn promote antimicrobial defense. While both scenarios would support the hypothesis that the virus is a mutualist, further experimentation will be needed to test whether and how viral infection persists due to the reinoculation of the virus through the fecal fluid of ants. Potential tests of this hypothesis could include virus-curing experiments, reinfection through the fecal fluid, and subsequent measures linking viral infection and L. gongylophorus performance.

The apparent occurrence of LgTlV1 across leafcutter ant fungi farmed by other leafcutter ant genera supports both the assumption of vertical viral transmission and the potential for long-term co-evolution with this virus since earlier periods of leafcutter ant evolution. It will be important to confirm the presence of these mycoviruses in: (1) L. gongylophorus fungi from field-collected leafcutter ant farming systems ranging across the biogeographic distribution of leafcutter ant species from Argentina to northern Texas, and (2) the fungal pellets carried by individual dispersing leafcutter ant queens. Our findings pave the way for testing these hypotheses about symbiotic functions within leafcutter ant farming systems and highlight that the study of mycoviruses is a nascent field with potentially major biological impact.

Availability of data and materials

Viral genome assemblies and RNASeq data have been uploaded to NCBI Genbank under accession numbers OR914225, OR914226, and SRR27078396 (BioProject PRJNA879936). All other data are available in the supplementary material.

References

  1. Ghabrial SA, Castón JR, Jiang D, Nibert ML, Suzuki N. 50-plus years of fungal viruses. Virology. 2015;479:356–68.

    Article  PubMed  Google Scholar 

  2. Yu X, Li B, Fu Y, Xie J, Cheng J, Ghabrial SA, et al. Extracellular transmission of a DNA mycovirus and its use as a natural fungicide. Proc Natl Acad Sci. 2013;110(4):1452–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Ghabrial SA. Origin, adaptation and evolutionary pathways of fungal viruses. Virus Genes. 1998;16(1):119–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Lerer V, Shlezinger N. Inseparable companions: Fungal viruses as regulators of fungal fitness and host adaptation. Front Cellul Infect Microbiol. 2022:1556.

  5. Márquez LM, Redman RS, Rodriguez RJ, Roossinck MJ. A virus in a fungus in a plant: three-way symbiosis required for thermal tolerance. Science. 2007;315(5811):513–5.

    Article  PubMed  Google Scholar 

  6. Redman RS, Sheehan KB, Stout RG, Rodriguez RJ, Henson JM. Thermotolerance generated by plant/fungal symbiosis. Science. 2002;298(5598):1581.

    Article  CAS  PubMed  Google Scholar 

  7. Hölldobler B, Wilson EO. The leafcutter ants: civilization by instinct. New York: WW Norton & Company; 2010.

    Google Scholar 

  8. Schultz TR, Brady SG. Major evolutionary transitions in ant agriculture. Proc Natl Acad Sci. 2008;105(14):5435–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Aylward FO, Khadempour L, Tremmel DM, McDonald BR, Nicora CD, Wu S, et al. Enrichment and broad representation of plant biomass-degrading enzymes in the specialized hyphal swellings of Leucoagaricus gongylophorus, the fungal symbiont of leaf-cutter ants. PLoS ONE. 2015;10(8): e0134752.

    Article  PubMed  PubMed Central  Google Scholar 

  10. De Fine Licht HH, Schiøtt M, Rogowska-Wrzesinska A, Nygaard S, Roepstorff P, Boomsma JJ. Laccase detoxification mediates the nutritional alliance between leaf-cutting ants and fungus-garden symbionts. Proc Natl Acad Sci. 2013;110(2):583–7.

    Article  PubMed  Google Scholar 

  11. Leal-Dutra CA, Yuen LM, Guedes BAM, Contreras-Serrano M, Marques PE, Shik JZ. Evidence that the domesticated fungus Leucoagaricus gongylophorus recycles its cytoplasmic contents as nutritional rewards to feed its leafcutter ant farmers. IMA Fungus. 2023;14(1):19.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Kooij PW, Rogowska-Wrzesinska A, Hoffmann D, Roepstorff P, Boomsma JJ, Schiøtt M. Leucoagaricus gongylophorus uses leaf-cutting ants to vector proteolytic enzymes towards new plant substrate. ISME J. 2014;8(5):1032–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Rønhede S, Boomsma JJ, Rosendahl S. Fungal enzymes transferred by leaf-cutting ants in their fungus gardens. Mycol Res. 2004;108(1):101–6.

    Article  PubMed  Google Scholar 

  14. Currie CR, Bot ANM, Boomsma JJ. Experimental evidence of a tripartite mutualism: bacteria protect ant fungus gardens from specialized parasites. Oikos. 2003;101(1):91–102.

    Article  Google Scholar 

  15. Heine D, Holmes NA, Worsley SF, Santos ACA, Innocent TM, Scherlach K, et al. Chemical warfare between leafcutter ant symbionts and a co-evolved pathogen. Nat Commun. 2018;9(1):2208.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Rodrigues A, Cable RN, Mueller UG, Bacci M, Pagnocca FC. Antagonistic interactions between garden yeasts and microfungal garden pathogens of leaf-cutting ants. Antonie Van Leeuwenhoek. 2009;96(3):331–42.

    Article  PubMed  Google Scholar 

  17. Sapountzis P, Zhukova M, Shik JZ, Schiott M, Boomsma JJ. Reconstructing the functions of endosymbiotic Mollicutes in fungus-growing ants. Elife. 2018;7: e39209.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Martiarena MJS, Deveau A, Montoya QV, Flórez LV, Rodrigues A. The hyphosphere of leaf-cutting ant cultivars is enriched with helper bacteria. Microb Ecol. 2023;86(3):1773–88.

    Article  CAS  PubMed  Google Scholar 

  19. Angeli-Papa J, Eymé J. Le champignon cultivé par la" fourmi-manioc", Acromyrmex octospinosus Reich en Guadeloupe; résultats préliminaires sur le mycélium en culturue pure et sur l’infrastructure des hyphes. CR Acad Sci. 1979;281:21–4.

    Google Scholar 

  20. Angeli-Papa J, Eymé J. Les champignons cultivés par les fourmis Attinæ - Évolution des structures cellulaires au cours du dévelopment. Annales des Sciences Naturelles, Botanique, Paris. 1985(13e série, Tome 7):103–29.

  21. Jo Y, Choi H, Chu H, Cho WK. Unveiling mycoviromes using fungal transcriptomes. Int J Mol Sci. 2022;23(18):10926.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. De Fine Licht HH, Boomsma JJ, Tunlid A. Symbiotic adaptations in the fungal cultivar of leaf-cutting ants. Nat Commun. 2014;5(1):5675.

    Article  PubMed  Google Scholar 

  23. Chen GY, Chen CY, Chang MDT, Matsuura Y, Hu YC. Concanavalin A affinity chromatography for efficient baculovirus purification. Biotechnol Prog. 2009;25(6):1669–77.

    Article  CAS  PubMed  Google Scholar 

  24. Mazalovska M, Kouokam JC. Lectins as promising therapeutics for the prevention and treatment of HIV and other potential coinfections. Biomed Res Int. 2018;2018:3750646.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Martelli GP, Sabanadzovic S, Abou-Ghanem Sabanadzovic N, Edwards MC, Dreher T. The family tymoviridae. Arch Virol. 2002;147(9):1837–46.

    Article  CAS  PubMed  Google Scholar 

  26. Ayllón MA, Vainio EJ. Chapter One - Mycoviruses as a part of the global virome: diversity, evolutionary links and lifestyle. In: Kielian M, Roossinck MJ, editors. Adv Virus Res, vol. 115. New York: Academic Press; 2023. p. 1–86.

    Google Scholar 

  27. Turina M, Hillman BI, Izadpanah K, Rastgou M, Rosa C, Consortium IR. ICTV virus taxonomy profile: ourmiavirus. J Gen Virol. 2017;98(2):129.

    Article  Google Scholar 

  28. Leal-Dutra CA, Vizueta J, Baril T, Kooij PW, Rødsgaard-Jørgensen A, Schantz-Conlon BH, Shik JZ. Genomic signatures of domestication in a fungus obligately farmed by leafcutter ants. bioRxiv. 2024. https://doi.org/10.1101/2024.07.02.601708.

    Article  Google Scholar 

  29. Vasimuddin M, Misra S, Li H, Aluru S, editors. Efficient architecture-aware acceleration of BWA-MEM for multicore systems. In: 2019 IEEE International parallel and distributed processing symposium (IPDPS); 2019 20–24 May 2019.

  30. Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, et al. Twelve years of SAMtools and BCFtools. GigaScience. 2021. https://doi.org/10.1093/gigascience/giab008.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Fu L, Niu B, Zhu Z, Wu S, Li W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 2012;28(23):3150–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinform. 2009;10(1):421.

    Article  Google Scholar 

  34. Buchfink B, Reuter K, Drost H-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat Methods. 2021;18(4):366–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6.

    Article  CAS  PubMed  Google Scholar 

  37. Larsson A. AliView: a fast and lightweight alignment viewer and editor for large datasets. Bioinformatics. 2014;30(22):3276–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35(2):518–22.

    Article  CAS  PubMed  Google Scholar 

  41. Zerbini FM, Siddell SG, Mushegian AR, Walker PJ, Lefkowitz EJ, Adriaenssens EM, et al. Differentiating between viruses and virus species by writing their names correctly. Arch Virol. 2022;167(4):1231–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucl Acids Res. 2011;39(Suppl_2):W29–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30(9):1236–40.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Lu S, Wang J, Chitsaz F, Derbyshire MK, Geer RC, Gonzales NR, et al. CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res. 2020;48(D1):D265–8.

    Article  CAS  PubMed  Google Scholar 

  45. Jumper J, Evans R, Pritzel A, Green T, Figurnov M, Ronneberger O, et al. Highly accurate protein structure prediction with AlphaFold. Nature. 2021;596(7873):583–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Li F, Jia M, Wang A. Hidden viral proteins: how powerful are they? PLoS Pathog. 2024;20(1): e1011905.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Gruber AR, Lorenz R, Bernhart SH, Neuböck R, Hofacker IL. The Vienna RNA Websuite. Nucl Acids Res. 2008;36(12):W70–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Ayllón MA, Turina M, Xie J, Nerva L, Marzano S-YL, Donaire L, et al. ICTV virus taxonomy profile: botourmiaviridae. J General Virol. 2020;101(5):454.

    Article  Google Scholar 

  49. Shi L, Wang J, Quan R, Yang F, Shang J, Chen B. CpATG8, a homolog of yeast autophagy protein ATG8, is required for pathogenesis and hypovirus accumulation in the chest blight fungus. Front Cell Infect Microbiol. 2019;9:222.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Sylvia Mathiasen for laboratory assistance, and Michael James Johnson for assistance with TEM. Data processing and analyses were performed using Computerome at the National Life Science Supercomputer at DTU (www.computerome.dk). TEM was performed at the Core Facility for Integrated Microscopy in the Faculty of Health and Medical Sciences at the University of Copenhagen. Zofia Magnowska and Nina Molin Høyland-Kroghsbo provided support with initial VLPs isolations. The Ministerio de Ambiente, Republica de Panama provided permits for leafcutter ant colony collection and exportation (SEX/A-31-12 and SEX/A-28-17).

Funding

Open access funding provided by Copenhagen University. This study was funded in part by the European Research Council (ELEVATE: ERC-2017-STG-757810), the Villum Foundation (VEX Grant No. 50281) and the Carlsberg Foundation (CF22-0664) to JZS. ERGRA and REM are CNPq Research fellows.

Author information

Authors and Affiliations

Authors

Contributions

ARJør, CLD, and JZS conceived the study. ARJør, ARJen and CLD performed laboratory experiments. CLD, ERGRA, and SFS performed genomic analyses and REM performed protein structure predictions. ARJør and CLD contributed equally to the manuscript draft with support from JZS. All authors provided feedback on and approved the final manuscript draft.

Corresponding author

Correspondence to Caio Ambrosio Leal-Dutra.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Supplementary file S1:

Supplementary methods, discussion, figures and tables.

Supplementary file S2:

Fasta file of total RNA assembly.

Supplementary file S3:

Amino acid sequences file of identified ORFs)from total RNA.

Supplementary file S4:

RdRp alignments used for phylogenetic reconstructions.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rødsgaard-Jørgensen, A., Leal-Dutra, C.A., de Santana, S.F. et al. Two +ssRNA mycoviruses cohabiting the fungal cultivar of leafcutter ants. Virol J 21, 211 (2024). https://doi.org/10.1186/s12985-024-02465-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-024-02465-0

Keywords