Viruses are the most abundant and diverse biological entities on Earth and have been found in virtually all host species studied to date, although understanding of the virosphere is still very limited (Paez-Espino et al. 2016; Shi et al. 2016). Characterizing and classifying this enormous diversity has been a difficult task. Originally, this task was limited to viruses that could be isolated in culture cells and induce cell death or alterations. In recent decades, however, high-throughput sequencing technologies have enabled an enormous expansion of knowledge about the viral component of the microbiome, the so-called viromes, which are based on the high-resolution viral sequences recovery (Wolf et al., 2018).
Despite increasing efforts to characterize the virome of different organisms and environments, our current understanding is biased for chordates, arthropods and plants. Platyhelminthes, a phylum composed of free-living and parasitic flatworms, is a poorly studied group in terms of its virome, even though several studies have described viruses on certain species of planarians, Schistosoma and Taenia (Jones MK and Whittington 1992; Justine and Bonami 1993; Noury-Sraïria et al 1995; Crespo-González et al 2008; Rebrikov et al, 2002). More recently, comprehensive studies analyzing platyhelminth genomes produced by NSG have found more than a hundred viruses infecting 50 species of free-living and 65 species of parasitic flatworms, including viruses belonging to the Picornavirales, Jingchuvirales, Maretellivirales, Flaviviridae, Nyamiviridae, Ghabrivirales, Rhabdoviridae and Bunyavirales families (Dheilly et al, 2022; Zhang et al 2023).
To our knowledge, no studies have been conducted to characterize viruses of the Catenulida, which are a more basal group of platyhelminthes (Telford et al., 2003; Wallberg et al., 2007). Catenulida are small free-living flatworms ranging in size from 0.5 to 2.5 mm that reproduce asexually by paratomy, a process in which the anterior body structures form in the mid-body regions, resulting in segmented organisms (Rosa et al., 2015). We have maintained a clonal culture of S. leucops (Catenulida) for more than a decade. This strain has a flexible relationship with its microbiome that may be associated with phenotypic plasticity (Rosa and Loreto 2023) and is capable of expressing genes found in its microbiome or eDNA (Rosa and Loreto 2019).
In the present study, we characterized a complete genome of a (+) ssRNA virus from the Picornaviridae family that is maintained for at least five years in a S. leucops colony, suggesting persistent and long-term infection. Phylogenetic analysis revealed that this virus is a novel virus closely related to the virus found in Biomphalaria. Based on the long-term association of this virus with S. leucops and the apparent lack of deleterious host effects, it is likely that this virus and its host engage in a mutualistic or commensalistic interaction.
Stenostomum leucops specimens were collected in 2009 in a pond in Santa Maria, Brazil (53°17′W; 29°28′S). An isoline was first established and the worms have been kept in reconstituted water at 28°C and fed with milk powder, following the details described by Rosa et al. (2015).
To gain deeper insights into this non-model species, we performed a transcriptomic characterization. Total RNA was isolated using Trizol (Invitrogen®) according to the manufacturer's guidelines. RNA quality was checked using a Nanodrop spectrophotometer (LabTech, USA) and RNA integrity was assessed using a 2100 Bioanalyzer (Agilent Technologies, USA). For library preparation, mRNA was enriched using the Dynabeads® mRNA purification kit (Invitrogen) and library preparation was performed using the Total Ion RNA-Seq v2 kit (Thermo Fisher Scientific). Sequencing was performed with the IonTorrent S5 sequencer (Thermo Fisher Scientific) using the Ion 540 TM Kit-OT2 and the Ion 540 TM Chip.
Viral sequences were also searched for in the genome sequences. For Oxford Nanopore long-read sequencing, DNA from approximately 30,000 worms was extracted using the protocol described by Sassi et al. (2005). Sequencing was performed on a PromethION-Oxford Nanopore platform from GenOne, RJ, Brazil. The genome and RNAseq sequences were deposited in GenBank (BioProject PRJNA1037460; accession SAMN38172808 and the S_leucopus_picornavirus sequence was deposited with acession PP818844).
For Illumina sequencing, genomic DNA was isolated from a pool of 100 individuals according to a previously described protocol (Oliveira et al 2009). A genomic DNA library was prepared from 200 ng of DNA using the Illumina TruSeq Nano DNA Kit. For sequencing, a MiSeq Reagent Kit v3 (600 cycles) was used to generate paired-end reads. Sequencing was performed using an Illumina MiSeq platform from Unidade de Genômica Computacional Darcy Fontoura de Almeida/LNCC/Brazil (BAM accession number: SAMN12500763).
To analyse the presence of the picornavirus genome found by RNAseq analysis, RT-PCR was performed. RNA was isolated using Quick-Zol (Ludwig Biotecnologia/Brazil), quantified using a QuBit 3.0 fluorometer (Thermo Fisher Scientific) and reverse transcribed using the High-Capacity cDNA Reverse Transcription Kit (Thermo Fisher Scientific) according to the manufacturer's recommendations. Two primer pairs were used: 1) F_Ls_pic_sit1: CACTGGTGCCAAAATCTG; R_Ls_pic_sit1: GTAGCTGGACGATACATAC; 2) F_Ls_pic_sit2: CGCACTTGGCATAGTACAC; R_Ls_pic_sit2: GATGTCAGGTTCGACTGAC; amplifying a 468bp and 386 bp fragment, respectively.
Viral search on the assembled transcriptome was performed using a pipeline described by Silva et al (2023). Briefly, pairwise alignment using Blastp from DIAMOND (Buchfink et al., 2015) was performed using the assembled transcripts against the viral database from NCBI (protein database NCB--txid10239- using the E-Value of 1E-3 and using the more-sensitive parameter. After was removed false positives matches running a blast with the contigs that passed the threshold blasting it against the NCBI non-redundant protein bank. The sequences that came back as viral positive and among the top 5 searches were used for the next step. Only viral transcripts larger than 600 bp were used to obtain the taxonomy information and to carry out the phylogenetic reconstruction.
The complete virus sequence and the RdRp amino acid sequence of the viruses were used for separate phylogenetic analysis together with the other Picornavirus found in Platyhelminthes described by Dheilly et al, (2022). To these sequences we also added some sequences showing high similarity to S. leucops virus found in the NR NCBI database. The accession number of the sequences used can be found at the Supplementary table 1).
The sequences were aligned using Clustal W and ML phylogeny was reconstructed with IQ-Tree (Minh et al, 2020 ) using a substitution model: GTR+F, and 10000 bootstrap replicates.
For estimates of Evolutionary Divergence between virus sequences, the number of base differences per site between sequences are conducted in MEGA X (Kumar et al, 2018) using gamma distribution. and bootstrap procedure (1000 replicates).
Analysis of the transcriptome of S. leucops revealed the complete sequence of a novel RNA virus. A search for a similar sequence in the genome of S. leucops revealed that the viral sequence is not present in the genome, indicating that it is a bona fide RNA virus. Blast searches in the NCBI nr and virus databank showed that the most similar sequences correspond to polyprotein YP_009342320.1 [Biomphalaria virus 2] and polyprotein AVA16916.1 [Pittsburgh sewage-associated virus 1], with an overall nucleotide divergence between these viruses of approximately 0.35 to 0.38, depending on the region of the virus (Table 1). These are (+) ssRNA viruses from the Picornaviridae family.
The complete sequence of S. leucops picornavirus has 8548 bp with three ORFs (Fig. 1). The first ORF encodes a polyprotein ranging from nucleotide 766 to 4425. This polyprotein corresponds to an RNA helicase, a peptidase C3 and a peptidase C3G. The second ORF extends from nucleotide 4401 to 7358 and encodes an RNA-dependent RNA polymerase (RdRp) and an rhv-like capsid protein. There is an overlap of 24 base pairs on the two ORFs. The third ORF expands from 4372 to 8412 nucleotides and encodes a CRPV capsid protein that belongs to a family of proteins found in positive stranded ssRNA. In addition, there is a 16 bp overlap between the second and third ORFs
Since there is a five-year gap between the time when the RNA extraction for transcriptome analysis was performed and the virome analyses using this data, we decided to develop a PCR assay to perform RT-PCR analysis in the S. leucops culture to see if the virus is still found in the S. leucops colony and has been maintained over this time period. As can be seen in Figure 2, the amplicons were produced at the expected size. We can show that this virus has been maintained for at least 5 years in the S. leucops isolates cultured in the laboratory.
Phylogenetic analyses were performed using the nucleotide sequence of the virus and the amino acid sequence of the viral gene RdRp, which is the most conserved sequence in these viral genomes. As can be seen in Figures 3 and 4, the topology of the trees is very similar. The sequences used gave rise to two main clades, which we refer to as clades A and B. Clade A includes several viruses present in flatworms and from other sources such as Pittsburgh sewage virus, Biophalaria virus and a water virus identified in a metagenomic analysis. S. leucops picornavirus_1 is grouped with these viruses from other sources. Group B includes mainly viruses found in Tricladida and also some viruses found in Trematoda. The sequences of Prostheceraeus (Policladida) and of Dugesia japonica and Procotyla (Tricladida) showed some differences in the positions of the phylogenetic analysis of nucleotides and proteins, sometimes outside the two main clades.
Previous studies have shown that Picornaviridae are the most representative viruses in free-living flatworms (Dheilly et al, 2022). Stenostomum leucops is a free-living flatworm from the order Catenulida, a basal group of platyhelminthes (Telford et al., 2003; Wallberg et al., 2007). The new virus described here (Sleuc_Pic_1), which harbours S. leucops, also belongs to the Picornaviridae family. Although the picornaviruses of the Platyhelminthes apparently form a monophyletic group, the picornavirus of S. leucops described here belongs to a group with other picornaviruses that do not occur in flatworms but harbour Biomphalaria (Mollusca) and were found in metagenomic analyses (water and sewage). Since Catenulida is a cosmopolitan group (Larson et al. 2008) that lives in the aquatic environment, similar viruses found in the metagenomic analyses are likely to be from Catenulida. However, the presence of viruses of the same clade in Biomphalaria (Mollusca) may also indicate that this particular group of picornaviruses may be more widespread in limnic organisms.
One point to emphasize is that in the genome of Sleuc_Pic_1 the ORFs show an overlap, resulting in three separate polyproteins CDSs. In general, picornaviruses have a single long ORF that encodes the polyproteins. However, there are cases where the polyproteins are encoded by multiple ORFs, similar to the virus we sequenced from S. leucops, for example, the picornavirus found in platyhelminths (accession numbers BK059711.1; BK059741.1; BK059754.1; and BK059682.1).
High-throughput sequencing has significantly improved the discovery of new viruses and increased our knowledge of viromes. However, this approach little inform regarding the interactions between the virus and its hosts. S. leucops reproduces asexually, which theoretically allows the virus to be passed from cell to cell and maintained over generations. We have shown that Sleuc_Pic_1 persisted for at least five years in the worms of a laboratory colony and thus represents a persistent infection. We found no adverse effects on S. leucops associated with the presence of this virus. Similarly, Burrows et al. (2020) have shown that a double-stranded RNA (dsRNA) totivirus-like virus is maintained as a persistent infection in an asexual strain of Schmidtea mediterranea planaria that has a mutualistic or commensal relationship with its host. These examples raise the question of whether other asexual flatworms also maintain viruses with these ecological interactions. A limitation of our study is that we did not examine stressful conditions or other factors that may activate or overactivate Sleuc_Pic_1 virus to observe whether it can cause harmful infection in worms. However, these questions can be addressed in future studies.