Next Article in Journal
Molecular Analysis of the Avian H7 Influenza Viruses Circulating in South Korea during 2018–2019: Evolutionary Significance and Associated Zoonotic Threats
Next Article in Special Issue
Composition of Eukaryotic Viruses and Bacteriophages in Individuals with Acute Gastroenteritis
Previous Article in Journal
CRISPR Screening: Molecular Tools for Studying Virus–Host Interactions
Previous Article in Special Issue
Discovery and Genetic Characterization of Novel Paramyxoviruses Related to the Genus Henipavirus in Crocidura Species in the Republic of Korea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 11
10.3390/v13112259
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Virome Analysis of Signal Crayfish (Pacifastacus leniusculus) along Its Invasion Range Reveals Diverse and Divergent RNA Viruses

by
Katarina Bačnik
1,2,
Denis Kutnjak
1,*,
Silvija Černi
3,
Ana Bielen
4 and
Sandra Hudina
3,*
1
Department of Biotechnology and Systems Biology, National Institute of Biology, 1000 Ljubljana, Slovenia
2
Jozef Stefan International Postgraduate School, 1000 Ljubljana, Slovenia
3
Department of Biology, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia
4
Department of Biochemical Engineering, Faculty of Food Technology and Biotechnology, University of Zagreb, 10000 Zagreb, Croatia
*
Authors to whom correspondence should be addressed.
Viruses 2021, 13(11), 2259; https://doi.org/10.3390/v13112259
Submission received: 1 October 2021 / Revised: 29 October 2021 / Accepted: 4 November 2021 / Published: 11 November 2021
(This article belongs to the Special Issue Applications of Next-Generation Sequencing in Virus Discovery)
Browse Figures
Versions Notes

Abstract

:
Crayfish are a keystone species of freshwater ecosystems and a successful invasive species. However, their pathogens, including viruses, remain understudied. The aim of this study was to analyze the virome of the invasive signal crayfish (Pacifastacus leniusculus) and to elucidate the potential differences in viral composition and abundance along its invasion range in the Korana River, Croatia. By the high-throughput sequencing of ribosomal RNA, depleted total RNA isolated from the crayfish hepatopancreas, and subsequent sequence data analysis, we identified novel and divergent RNA viruses, including signal crayfish-associated reo-like, hepe-like, toti-like, and picorna-like viruses, phylogenetically related to viruses previously associated with crustacean hosts. The patterns of reads abundance and calculated nucleotide diversities of the detected viral sequences varied along the invasion range. This could indicate the possible influence of different factors and processes on signal crayfish virome composition: e.g., the differences in signal crayfish population density, the non-random dispersal of host individuals from the core to the invasion fronts, and the transfer of viruses from the native co-occurring and phylogenetically related crayfish species. The study reveals a high, previously undiscovered diversity of divergent RNA viruses associated with signal crayfish, and sets foundations for understanding the potential risk of virus transmissions as a result of this invader’s dispersal.

1. Introduction

Until recently, knowledge of invertebrate viruses was limited mostly to viral pathogens causing high mortalities and arboviruses (arthropod-borne viruses), which are vectored by arthropods and cause disease in humans and other vertebrate species [1]. This view has changed, with metagenomics revealing remarkable levels of RNA virus diversity in invertebrates [2,3,4,5], and, in this perspective, the detection of invertebrate disease-causing viruses became the exception rather than the rule [1]. In addition to transforming our understanding of virus diversity, virus metagenomics provides insights into different aspects of virus evolution [6], including the role of invertebrates and their ecological interactions with other organisms in the evolution of RNA viruses [5].
For understanding the diversity of invertebrate viruses, arthropods, including crustaceans, may be of particular importance, as they account for over 80% of the total animal diversity [7], are globally abundant, and often live in extremely large and dense populations [6].
Crayfish, belonging to a diverse order of decapod crustaceans, are a keystone species of freshwater ecosystems, and their populations can reach high densities. They are an important component of freshwater food webs and ecosystem engineers due to their larger size and longer life span compared to other benthic macroinvertebrates, as well as due to their bioturbation and burrowing activity [8]. Crayfish are also among the most successful invasive alien species (IAS) in freshwater ecosystems [9]. Throughout history, crayfish have been frequently translocated for ornamental or aquaculture purposes [10] and, after being intentionally or unintentionally released into the wild, have often rapidly expanded their range and exerted a high number of documented negative impacts [10,11]. On the other hand, native crayfish populations in the wild are increasingly endangered on the global scale [9,10,12,13], and invasive crayfish and their pathogens have been identified as one of the most prominent factors contributing to the decline of native species’ populations [14,15]. Some of the pathogens introduced by invasive crayfish are significantly more virulent towards the new (i.e., native) hosts than toward IAS hosts [16,17]. Thus, invasive crayfish species have a high potential for introducing and spreading the emerging diseases that may incur significant economic and ecological losses [18,19]. However, the existing literature on crayfish diseases is biased mostly towards the crayfish plague pathogen, the oomycete Aphanomyces astaci, while other disease agents, including fungi, bacteria, and viruses, are significantly understudied [20].
Most of the knowledge on crustacean pathogenic viruses comes from the study of viral diseases in aquaculture [20]. For instance, white spot syndrome virus causes the white spot disease that has mostly been described due to its symptoms in aquaculturally relevant species such as Cherax quadricarinatus [21], although it displays a low degree of host specificity, and is therefore also an important pathogen in crayfish populations in the wild [22]. Despite their large potential impact on aquaculture and populations in the wild, viruses are the least studied group of crayfish pathogens [20], and relatively few have been formally characterized and classified [23]. However, with the increasing availability of high-throughput sequencing (HTS) technologies and metagenomics approaches, which can reveal all of the viruses (i.e., virome) associated with an individual or an environmental sample, comprehensive descriptions and a revolution in our understanding of the viral diversity of crustacean species are occurring [23]. For instance, this approach was recently employed to discover bunya-like brown spot virus, a new negative, single stranded RNA virus belonging to the Phenuiviridae family that was associated with a massive disease outbreak in the population of the endangered white-clawed crayfish Austropotamobius pallipes [24].
The aim of our study was to analyze the virome of the signal crayfish, Pacifastacus leniusculus, (Dana, 1852), one of the most successful freshwater invaders in Europe. The signal crayfish is present in over 29 EU countries [25], reaches very high dispersal rates within European watercourses [26], and is listed as IAS of Union concern (EU Regulation No. 1143/2014 on invasive alien species). We have investigated the signal crayfish population in a recently invaded Korana river, Croatia, where signal crayfish range expansion is well monitored [27,28,29], and occurs in both upstream and downstream directions [27]. Recently established populations at invasion fronts have 6–11 times lower crayfish abundance in comparison to longer established populations in invasion cores, and also contain co-occurring native congener (narrow-clawed crayfish, Pontastacus leptodactylus, Eschscholtz, 1823) [27]. Host population density, with increased contact rates among individuals in highly dense populations, is an important factor affecting rates of virus transmission [30,31].
In this study, we sampled the signal crayfish individuals at four locations along its invasion range, including upstream and downstream invasion cores and invasion fronts. We sequenced and analyzed signal crayfish samples to: (1) identify putative viral sequences in the signal crayfish hepatopancreas, as this is the tissue most often associated with viral infections in freshwater crayfish [20]; (2) describe the phylogenetic relationships of the most abundant novel signal crayfish-associated viruses; (3) discuss the potential differences in virome composition and the abundance of viral sequences along the invasion range. Our results reveal a high, previously undiscovered diversity of viruses in signal crayfish, with putative signal crayfish-associated reo-like, hepe-like, toti-like, and picorna-like viruses, representing divergent sequences most similar to the viruses previously associated with crayfish hosts.

2. Materials and Methods

2.1. Study Area and Sample Collection

Research was carried out in the continental part of Croatia, in the Korana River. Here, the signal crayfish was illegally introduced in the lower section of the river, and is spreading both upstream and downstream [28], with its invasion range currently stretching along 33 km [27]. Significant differences in crayfish abundance occur along the signal crayfish invasion range, i.e., between the longer established populations in the center of its distribution (invasion core) compared to the recently established populations at distribution edges (invasion fronts). Invasion front locations have up to 11 times lower crayfish abundance compared to the invasion core [27]. Also, at invasion fronts, interspecific populations of the signal crayfish and the native narrow-clawed crayfish (Pantastacus leptodactylus, Eschscholtz, 1823) are present, with 2.5–12.9 higher abundances of the native P. leptodactylus compared to the signal crayfish [27]. Populations in the invasion cores comprise only signal crayfish, which have gradually displaced the native P. leptodactylus [27,29], as depicted in Figure 1.
Crayfish were collected in the period of increased crayfish activity (early autumn 2018) using baited LiNi traps [32]. A total of 110 signal crayfish individuals were captured from four invasion range endpoints (Figure 1), which were previously [27] identified as: the upstream invasion front (UF), upstream invasion core (UC), downstream invasion core (DC), and downstream invasion front (DF). For each location, we calculated the catch per unit effort (CPUE) that equals the number of crayfish captured per each trap per trapping night, and is a frequently used measure of relative crayfish abundance (Figure 1). The captured indigenous crayfish species were returned to the same location where they were caught, while signal crayfish individuals from all four sampling locations were taken on ice to the laboratory where hepatopancreas samples were randomly collected from 25 individuals from each location. For each individual, the complete organ was removed from the body, placed in a sterile Petri dish, and carefully chopped into smaller pieces using a sterile scalpel. The cut pieces of the organ were stored at −80 °C in an RNA stabilizing agent (RNA later; Sigma Aldrich, MO, USA) until RNA extraction.

2.2. RNA Extraction and Sequencing

The total RNA was extracted from hepatopancreas tissue using the RNeasy Lipid Tissue Mini Kit (Qiagen CA, USA). Before RNA extraction, equal amounts of tissue (30 mg) of each of the 25 individuals collected at the same location were pooled together, resulting in four composite pools corresponding to four sampling locations. After homogenizing tissue in liquid nitrogen, RNA extraction, followed by on-column DNase digestion step (RNase-Free DNase Set, Qiagen), was performed as recommended by the manufacturer. The eluted RNA from all four samples was sent for ribosomal RNA depletion (Illumina Ribo-Zero Plus rRNA Depletion Kit), sequencing library preparation (NEBNext Ultra RNA Library Prep Kit), and shotgun sequencing (Illumina HiSeqX, 2 × 150 bp) to CD Genomics, USA.

2.3. Bioinformatic Analysis for Detection of Viral Sequences

Sequencing datasets were quality checked, filtered, and analyzed as described below. The overview of the bioinformatic analyses performed is schematized in Supplementary material, Figure S1. The remains of the sequencing adaptors were trimmed, and resulting reads were filtered by using a quality filter (Limit = 0.01; no ambiguous nucleotide allowed) and by size (reads shorter than 25 bp were discarded) in CLC Genomics Workbench 20 (Qiagen Bioinformatics, CA, USA). Trimmed, size, and quality filtered reads were exported from the CLC Genomics Workbench and compared for similarity with a complete NCBI nr database (2020) using Diamond blastx [33] with default parameters, followed by the taxonomic classification of the reads and visualization of the results using Megan (Metagenome Analyzer, version 6.20.19, Tübingen, Germany) [34].
To identify longer virus-like sequences, reads were used for de novo assembly using SPAdes (version 3.14.0, St. Petersburg, Russia) [35], and the de novo constructed contigs were queried for similarity on protein level using Diamond blastx [33], followed by the taxonomic classification of the reads and visualization of the results using Megan [34]. The contig sequences that were classified as viral in all four analyzed samples were imported into the CLC Genomics Workbench, where additional steps of de novo assembly of viral contigs from all four samples was performed, resulting in a further assembly of overlapping sequences from different samples. The contigs that were not classified as viral by Diamond-Megan analysis were further analyzed using conserved domain search against Pfam (v32) database, and additional virus-like sequences detected using this analysis were added to the list of virus-like sequences from Diamond-Megan analysis.
The contigs that were identified as virus-like sequences according to the described analysis (Figure S1) were further investigated using blastn (NCBI-nt, April 2021) analysis to identify any possible misclassified sequences. After the additional manual inspection of all the similarity searches, we classified virus-like sequences into different categories (Table S2). Non-viral hits, such as host sequences misclassified as viral sequences (similar to viruses from Nudiviridae family) and possible endogenized virus-like sequences (e.g., members of Retroviridae and Adintoviridae families), as well as viral sequences belonging to bacteriophages (e.g., members of Caudovirales order, Leviviridae family), plant infecting viruses (e.g., members of Solemoviridae and Bromoviridae families), and sequences identified as laboratory contamination (e.g., Cryphonectria hypovirus) were excluded from further analysis. The final selection of virus-like contigs included only sequences of potential invertebrate viruses longer than 300 nts. The differences between the abundance of reads corresponding to selected virus-like contigs detected in individual samples from different locations were identified by mapping the reads to virus-like contig sequences in CLC Genomics Workbench 20 (length fraction = 0.90, similarity fraction = 0.90), and the percentage of mapped reads was calculated (Table S1). In order to normalize samples according to the number of reads before mapping them to contig sequences, random subsampling was done and all samples were normalized to 97,563,944 reads, which was the minimum number of reads obtained among the samples. Decimal logarithm transformed average read coverage values were visualized as a heatmap, and used for clustering the locations and viruses using the heatmap.2 function in a gplots package in R 4.0.5 [36]
All sequencing data that support the findings of this study are linked under the accession number PRJNA754774 in the NCBI BioProject database. The sequences of selected viral contigs identified in this study are available in Table S2 and have been deposited in the GenBank under the accession numbers OK317706-OK317734.

2.4. PCR Detection of Selected Viral Contigs in Hepatopancreas Samples

Based on the presence of RNA-dependent RNA polymerase (RdRp) domain, contig length, and average read coverage, we selected four viral contigs for further analysis; virus-like contig 4 (signal crayfish-associated reo-like virus 1), virus-like contig 139 (signal crayfish-associated hepe-like virus 1), virus-like contig 1 (signal crayfish-associated picorna-like virus 1), and virus-like contig 141 (signal crayfish-associated toti-like virus 1). To confirm their presence in the samples, primer pairs (Table S3) amplifying the fragments of newly detected sequences were constructed using Primer-BLAST [37]. cDNA was synthesized from 2 µL of total RNA in a reaction volume of 20 µL by random priming using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, CA, USA) as described by the manufacturer, using the following reaction conditions: 25 °C for 10 min, 37 °C for 120 min, and 85 °C for 5 min. Target fragments were amplified in a total volume of 25 µL using 1 µL of cDNA template, 2 mM MgCl2, 0.2 mM dNTP mixture, 0.5 µM each primer, and 0.0625 U of G2 polymerase (GoTaq Flexi DNA Polymerase, Promega, WI, USA) and following cycling conditions: 95 °C for 2 min, 40 cycles of 95 °C for 30 s, 58 °C for 30 s, and 72 °C for 30 s; followed by 5 min of final extension at 72 °C. Amplicons were detected after electrophoresis in 1% agarose gel stained with SYBR Safe DNA Gel Stain (Invitrogen, CA, USA).

2.5. Phylogenetic Analyses

For selected virus-like contigs 4, 139, 1, and 141, the most similar homologous sequences with corresponding RdRp sequences were retrieved from the NCBI using blastp search against NCBI nr database (May 2021), and used in the phylogenetic analysis. The RdRp protein sequences of recognized viral species belonging to reo-like, hepe-like, picorna-like and toti-like virus groups selected according to previous studies [3] were additionally included in the phylogenetic analysis from the GenBank database. RdRp amino acid sequences were aligned using the L-INS-I algorithm of the MAFFT program [38], and alignments were trimmed manually and with TrimAI (automated 1 mode) [39]. Alignments were visualized in MEGAX [40]. Sequences included into the alignments with corresponding GenBank accession numbers are given in Alignments S1–S4. Maximum likelihood trees were inferred using PhyML (v. 3.00, Montpellier, France) [41] available in the Phylemon 2.0 web server [42] employing the Le Gascuel (LG + I + G) amino acid replacement model, selected by IQTree ModelFinder [43], using 1000 bootstrap replicates.

2.6. Analysis of Nucleotide Diversities

To compare the possible differences in the diversities of populations of selected viruses detected at all sampling locations, we determined the nucleotide diversities of populations of three signal crayfish-associated viruses, which were present at all sampling locations with an average read coverage higher than 40. For these analyses, read data was additionally filtered. First, reads with low nucleotide quality (<20) were discarded using a FASTQ Quality Filter (FASTX-toolkit). Reads from individual samples were then mapped to selected viral contigs (length fraction = 0.90, similarity fraction = 0.90). Before mapping, the read datasets were subsampled to amount for the same read coverage of individual contig sequence in different samples (higher read coverage was normalized by using a lower number of reads per sample). Single-nucleotide polymorphism (SNP) calling was performed using a Low Frequency Variant Detection tool in the CLC Genomics Workbench 20 (required significance = 1%, minimum coverage = 10, minimum count = 2, minimum frequency = 1%). SNP tables were used as an input for SNPGenie software [44], which was used to calculate diversity indices for each viral contig—location combination. Calculated nucleotide diversity (π) values, which represent the mean number of pairwise differences per nucleotide in a population of mapped reads, were plotted and compared according to the location.

3. Results and Discussion

3.1. Novel Viral Sequences Identified in Hepatopancreas of Signal Crayfish

3.1.1. Overview of the Newly Identified Viral Sequences in Hepatopancreas of Signal Crayfish

To identify viral sequences, we performed sequencing of ribosomal RNA depleted total RNA resulting in high numbers (100–164 M) of reads per sample (Table 1). First, in a reads-based protein similarity search, a small fraction of reads (at most, 0.018%) was classified as viral sequences (Table 1). Read sequences similar to different RNA viruses (classified as Riboviria) were the most abundant viral sequences in all four samples. Our analysis of de novo assembled contigs, presented in Table S2 together with blast resulting similarities and average read coverage values, supported the observed presence of sequences similar to different unclassified RNA viruses. Results revealed the presence of diverse virus-like sequences, with selected virus-like contigs (longer than 300 nts) representing sequences associated with invertebrate hosts shown in Table 2.
Sequences belonging to potential invertebrate viruses exhibited relatively low levels of similarity to known viral sequences, with sequence identities ranging from 25% to 84% according to blastx results. Many virus-like contigs representing novel, highly divergent RNA viral sequences were most similar to reo-like, hepe-like, toti-like, picorna-like, tombus-like, chu-like, and partiti-like viruses. A small fraction of viral reads and the dominance of previously unknown RNA viruses in aquatic arthropod hosts samples is in line with previous studies [3,5,45]. Many of the virus-like sequences identified in this study had the highest identities with unclassified RNA viruses discovered as part of a large-scale metagenomics study of invertebrate viromes [3], where similarities with viruses described in different crustacean hosts were also observed.
Some viral contig sequences (virus-like contigs 139, 141) discovered in this study likely correspond to full-length genomes, given that the length and genomic organization is similar to the annotated closest hits from the blastx analysis. High average read coverage values were obtained for some virus-like contigs (Figure 2), which may further support the assumption that they represent crayfish infecting viruses. However, confidently assigning host associations based on metagenomics studies and bioinformatic analyses remains a challenge [1], and therefore, some virus-like contigs identified here may be associated with other organisms that are present within the studied host [46].

3.1.2. Phylogenetic Relationships and Genome Organization of Selected Newly Identified Viral Sequences

Based on the presence of the RdRp domain and contig length, we selected four distinct viral contigs for further analyses. They represent putative novel virus species belonging to reo-like (virus-like contig 4), hepe-like (virus-like contig 139), toti-like (virus-like contig 141), and picorna-like (virus-like contig 1) virus clades. These viral sequences showed similarities to viruses previously described in crustacean hosts, suggesting that they were likely infecting signal crayfish rather than being associated with the environment or co-infecting microorganisms. The presence of signal crayfish-associated reo-like virus 1, hepe-like virus 1, and toti-like virus 1 sequences was confirmed in all four analyzed samples by obtaining the products of the expected size after RT-PCR (Figure 3), while the presence of signal crayfish-associated picorna-like virus 1 was confirmed in the UF sample only, in line with the sequencing results (Figure 2).
In our study, signal crayfish-associated reo-like virus 1 (virus-like contig 4, 4234 nts) had the highest average read coverage in all samples (Figure 2). A variety of reoviruses have been described previously in crustacean hosts [23]. Using blastx, we observed a distant similarity to the RdRp partial coding sequence of Cherax quadricarinatus reovirus, which was recently associated with hepatopancreatic samples in redclaw crayfish (Cherax quadricarinatus) from Australia [47]. This relationship was also confirmed by phylogenetic analysis (Figure 4). Cherax quadricarinatus reovirus causes the necrosis of hepatopancreocytes and inflammatory cells in hepatopancreatic tubules [48]. Notably, we have recorded a similar condition, the acute necrotizing hepatopancreatitis, along the whole invasion range of the signal crayfish in the Korana River, which is currently classified as idiopathic [49]. Thus, future studies should examine the recorded signal crayfish-associated reo-like virus 1 as a potential causative agent of this disease.
Signal crayfish-associated hepe-like virus 1 (virus-like contig 139, 10,400 nts), detected in all samples, is represented by near complete genome sequence with four open reading frames (ORF), including ORF1—nonstructural polyprotein with RdRp domain (Figure 5). In blastx analysis, this contig had the highest (44%) amino acid identity to sequence of Beihai hepe-like virus 4, identified from the blue swimmer crab (deposited sequence is 11,635 nts long) and octopus (deposited sequence is 12,648 nts long) in China [3]. Consistent with the blastx search, phylogenetic analysis showed clustering of the detected virus sequence with Beihai hepe-like virus 4 and other hepe-like viruses from different crustacean hosts (Figure 5). Aside from reports in large-scale metagenomics studies [3], novel hepe-like virus has recently been characterized from the cephalothoraxes of the economically important giant freshwater prawn Macrobrachium rosenbergii with growth retardation, however, it could not be determined as the causative agent of the disease [50].
The signal crayfish-associated toti-like virus 1 (virus-like contig 141, 8576 nts), containing four ORFs, was clustered with the Wenzhou crab virus 5 (deposited sequence is 8691 nts long) in a phylogenetic analysis (Figure 6). It had a 41% amino acid identity in blastx analysis with the same virus, which was associated with the Charybdis crab host [3] and was detected in all of our samples. Similarities to other viruses, especially toti-like viruses, from different insect hosts were also detected using blastx. Other viruses from toti-like virus clade, with no blastx similarity to sequence identified in this study, and not included in phylogenetic analysis, are known to be associated with crayfish; infectious myonecrosis virus (IMNV), the causative agent of infectious myonecrosis in the Pacific white shrimp (Litopenaeus vannamei) [51], and the Cherax giardiavirus-like virus (CGV) in freshwater crayfish (Cherax quadricarinatus) [52]. However, IMNV targets mesoderm-derived tissues, like muscles, connective tissue, and hemocytes, and does not replicate in enteric tissues like hepatopancreas [53]. CGV was, however, reported by histology in hepatopancreas, similar to the toti-like virus identified in this study, but it has been neither sequenced nor isolated [52,53].
Distinct contigs similar to previously described picorna-like viruses, a loosely defined broad group of viruses, were detected mostly in UF, with the highest blastx similarities to different crustacean picorna-like viruses (Table 2), such as Changjiang picorna-like virus 6, Beihai picorna-like virus 99, and Wenzhou picorna-like virus 38 [3]. The longest (4587 nts) picorna-like virus contig (virus-like contig 1, signal crayfish-associated picorna-like virus 1) with the highest read coverage had 26% amino acid identity in the blastx analysis with the Beihai picorna-like virus 99 (deposited sequence is 6503 nts long) isolated from a hermit crab in China [3]. Based on phylogenetic analysis and genome organization (Figure 7), with ORF1 nonstructural polyprotein comprising RdRp domain, signal crayfish-associated picorna-like virus 1 represents a single-sequence branch placed in the cluster that contains viruses of different invertebrate hosts, including crustacean viruses. In the same cluster, Taura syndrome virus is a known disease-causing crustacean virus that has spread rapidly throughout shrimp farming regions globally [23], and was reported to induce significant changes in the hepatopancreas transcriptome of the Pacific white shrimp [54]. Besides most similar viruses included in our phylogenetic analysis, mud crab dicistrovirus has been described previously as pathogenic to the mud crab (Scylla paramamosain) [55]. Other sequences distantly similar to picorna-like viruses (virus-like contigs 10, 11, 15, 9, 5, 13, 116, 2, 36) were also identified with low average read coverage, and could not be further assembled into longer contigs.
Further viral contigs potentially associated with signal crayfish, often detected in only one sample and with low average read coverage, belonged to tombus-like, chu-like, partiti-like, sobemo-like, and narna-like virus groups. Sequences similar to tombus-like viruses were previously associated with plant hosts, however, numerous tombus-related viruses have recently been described in non-plant hosts, such as marine and terrestrial invertebrates. Based on the phylogenetic analyses, it was suggested that this group of viruses was primarily associated with invertebrates, and was later horizontally transferred from aquatic invertebrates to plant hosts [3,6]. Apart from HTS based studies, tombus-like viruses have not been widely studied as infectious agents in crustacean hosts, except for the mud crab tombus-like virus that was associated with the sleeping disease in crabs, but only as a part of co-infection with mud crab dicistrovirus [56]. In this study, distinct tombus-like virus sequences were detected in signal crayfish hepatopancreas. Signal crayfish-associated tombus-like virus 1 (virus-like contig 66, 4504 nts) with a near complete genome sequence found in the DC location was most similar to the Caledonia beadlet anemone tombus-like virus 1 (41% identity on amino acid level) isolated from the anemone species Actinia equine [57]. Signal crayfish-associated tombus-like virus 2 (virus-like contig 84, 2981 nts) was detected in the DF location and was also most similar to Caledonia beadlet anemone tombus-like virus 1 (39% identity on amino acid level) [57]. These two tombus-like virus contigs had an average read coverage higher than ten, in contrast to lower read coverage of additional tombus-like virus contigs (virus-like contigs 55, 35, 169) (1425, 665, 301 nts respectively) detected in the UC location with the highest blastx-based identities with different invertebrate and plant infecting tombus-like viruses. Low read coverage tombus-like contigs may represent plant-infecting viruses, derived from food sources, while those with higher read coverage are more likely to be true crayfish viruses.
At all locations, different contigs similar to the Beihai hermit crab virus 3 [3] were identified. Virus-like contig 140 (signal crayfish-associated chu-like virus 1, 2216 nts) was the longest, and had a 27% amino acid sequence identity with hypothetical protein 2. Virus-like contig 145 (signal crayfish-associated chu-like virus 4, 1009 bp) had a 30% amino acid sequence identity to RdRp. The Beihai hermit crab virus 3 belongs to the Mivirus genus from the recently discovered negative-strand RNA Chuviridae family, characterized by HTS analysis only [3,31,58]. Diverse chuvirus-derived endogenous viral elements were also detected in mosquito genomes [59]. Additional contigs (virus-like contigs 65, 7, 146, 222) exhibited similarity with the Beihai hermit crab virus 3 hypothetical protein 2 or RdRp, but could not be further assembled into longer contigs, and did not contain full-length RdRp domain. In the UF location we have also detected a signal crayfish-associated partiti-like virus 1 (virus-like contig 27, 1188 nts) with a 56% amino acid identity to Caledonia partiti-like virus isolated from the breadcrumb sponge (Halichondria panacea) [57]. According to the blastx results, signal crayfish-associated partiti-like virus 1 was also similar to other crustacean-associated viruses (Wenling partiti-like virus 13, Wenling partiti-like virus 11) [3] and may with relatively high average read coverage (28) represent a crayfish-infecting virus.

3.2. Abundance of Viral Sequences and Nucleotide Diversity of Selected Virus Populations along the Invasion Range

We have compared the abundance of viral sequences as well as the nucleotide diversity of selected virus populations at different sampling locations in order to obtain insights into the possible patterns of variations along the invasion range. Based on the observed patterns, we could postulate hypotheses about how the observed patterns could be connected with ecological conditions at different sample sites, which represent a baseline for future research work on the virus ecology of crayfish and/or invasive alien species.
In our study, an abundance of virus-like reads, distinct virus-like contigs, and sequences mapping to signal crayfish-associated virus-like contigs varied between locations along the invasion range, with virus-like reads constituting between 0.0016% of the total reads in UF and 0.0174% in DC (Table 1). When comparing the two upstream and two downstream samples, the ones at the invasion cores had higher numbers of virus-like reads detected in comparison to those from the corresponding front sites. In addition, higher read coverage of viral sequences belonging to two out of three of the most represented newly identified signal crayfish-associated viruses was observed in the invasion cores compared to invasion fronts (Figure 2). At the same time, higher crayfish abundances were observed at the invasion cores (Figure 1).
The differences in signal crayfish abundance along the invasion range could influence the abundance of virus sequences detected, according to the classical epidemiological theory that links larger populations with higher contact rates with an increased likelihood of viral transmission [60,61,62]. The UF location differed the most from other locations in several aspects. The lowest percentage of virus-like reads was detected, and also lower reads coverages for signal crayfish-associated reo-like virus 1 (virus-like contig 4), crayfish-associated hepe-like virus 1 (virus-like contig 139), and toti-like virus 1 (virus-like contig 141) were found for this locality (Figure 2). Also, the lowest crayfish abundance was estimated at this locality, which probably resulted in lower encounter rates between individuals. Additionally, in the UF, we have recorded a lower severity of the acute necrotizing hepatopancreatitis [49]. This correlates with above mentioned lower read abundance at this location, and implies that some of the detected putative novel viruses in this study may be associated with the observed pathological changes in hepatopancreas. To test this hypothesis, further studies would be needed to associate individual crayfish samples of a known health status with the presence of putative novel viruses discovered in this study.
Despite the lower total number of virus-like reads in the UF sample, the highest number of distinct putative virus-like contigs similar to unclassified RNA viruses (Table S2) was found here, and picorna-like virus contigs were found predominantly in this sample. According to previous research [63,64,65], a decrease in viral diversification is expected in less dense groups, such as invasion fronts. In contrast, we observed the highest number of distinct viral contigs in the UF location, a site with the lowest crayfish abundance. However, at invasion fronts, the population of the native P. leptodactylus is also present with 2.5–12.9 higher abundances than the signal crayfish [27], potentially contributing to the inter-specific transmission of viruses. The observed differences in the number of total virus-like reads between the samples could also partially result from methodological limitations, such as difficulties in detecting highly divergent viral sequences that are considerably different from the reference genomes of known viruses [66].
For the three viruses, which were present at all four locations with relatively high average read coverage, we determined the diversities of their populations by calculating nucleotide diversities. This enabled us to estimate how the diversity of each discovered viral species population might be changing along its host invasion range. For all analyzed viruses, we observed differences in nucleotide diversity between populations at the invasion cores and invasion fronts. We observed two contrasting patterns. For signal crayfish-associated reo-like virus 1, the nucleotide diversity was decreased at the invasion fronts, while for signal crayfish-associated hepe-like virus 1 and toti-like virus 1, the nucleotide diversity was higher at the invasion fronts (Figure 8).
Observed differences in nucleotide diversities of viral populations of the three investigated viruses from samples collected along the signal crayfish invasion range might result from different processes, and several hypotheses can be put forward to explain the observed patterns. For instance, higher nucleotide diversity might indicate a higher prevalence of a particular virus, e.g., reo-like virus in the core sites. However, since we were analyzing composite hepatopancreas samples, we are limited in our estimations of the prevalence of individual viruses. Furthermore, we speculate that the estimated nucleotide diversity (π) values potentially reflect transmission dynamics of analyzed viruses, with a higher core nucleotide diversity of signal crayfish-associated reo-like virus 1 suggesting the transmission of the virus from longer established core populations with high signal crayfish abundance to recently established front populations of low signal crayfish abundance. The signal crayfish-associated reo-like virus 1 might be an endemic virus of signal crayfish, introduced to the new habitat along with its host, which could be confirmed/rejected by detecting the signal crayfish-associated reo-like virus 1 sequence in signal crayfish from other geographic regions. On the contrary, a higher nucleotide diversity of signal crayfish-associated toti-like virus 1 and hepe-like virus 1 in samples from invasion fronts could potentially reflect the introduction of these viruses from other host populations (e.g., from the native P. leptodactylus with higher density at the fronts) to populations of low signal crayfish abundance at invasion fronts. Successful inter-specific transmission relies on physical opportunities for transmission, such as occupying the same environment at the same time [67]. At invasion cores, the signal crayfish have displaced the native host, which would, under the proposed scenario, explain the decrease in nucleotide diversity for signal crayfish-associated toti-like virus 1 and hepe-like virus 1. However, such host associations should be further investigated, since very few of the crustacean viruses have been tested for infectivity in species other than the original host [68]; this hypothesis could be tested by analyzing native crayfish species for the presence of signal crayfish-associated toti-like virus 1 and hepe-like virus 1. Another possible hypothesis for the observed patterns might be a non-random dispersal of the host individuals during the invasion process, as demonstrated for a number of invasive species (i.e., [69]).
Finally, the results regarding the comparison of different locations can be confounded by other factors, which we were not able to elucidate using our experimental design (sequencing of samples pooled by location), and further analysis of individual samples would be needed to advance the understanding and test the postulated hypotheses.

4. Conclusions

  • We reported, for the first time, the virome of signal crayfish hepatopancreas tissue and found a high diversity of novel divergent viral sequences most similar to different unclassified RNA viruses.
  • We identified putative novel RNA viruses, including near complete genome sequence of signal crayfish-associated hepe-like virus 1 and toti-like virus 1, and the partial genomes of signal crayfish-associated reo-like virus 1 and picorna-like viruses. We identified additional tombus-like, partiti-like, and chu-like virus sequences potentially representing novel crayfish viruses. This pioneer study represents a baseline for the future research of a signal crayfish virome, e.g., to confirm the association of novel viruses with signal crayfish host, and to investigate their potential involvement in the observed necrotizing hepatopancreatitis.
  • We speculate that the differences in the signal crayfish population density along the invasion range, non-random dispersal, and possibilities of inter-specific viral transmissions may have an effect on the diversity and abundance of signal crayfish-associated viral sequences. Different hypotheses can be postulated to explain these patterns, and this study represents a baseline for the further research of virus transmission dynamics as a result of the invader’s fast dispersal, including inter-species transmission between the signal crayfish as an invader and P. leptodactylus as a co-occurring and phylogenetically related native species.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13112259/s1, Figure S1. Schematic representation of different steps of the bioinformatic analysis of sequencing reads (black arrows) and contigs (green arrows) to identify virus-like sequences in the signal crayfish hepatopancreas samples; Table S1. Selected putative invertebrate virus-like contigs (>300 nts) representing signal crayfish associated viruses identified in this study with their contig length (nts) and percentage of mapped reads resulting from mapping of reads from individual samples to virus-like contigs; Table S2. (A) Virus-like contigs from different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front) identified in Diamond analysis together with average coverage values, Diamond classification, blastn similarity search results (NCBI-nt, April 2021) and sequences of individual virus-like contigs. (B) virus-like contigs from different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front) identified using pfam domain search, where the remaining contigs not classified as viral by Diamond (963721) were translated and compared with the entire Pfam database, together with average coverage values, Pfam search results and blastn (NCBI-nt, April 2021) and blastx (NCBI-nr, June 2021) similarity search results; Table S3. Sequences of primers used for PCR amplification of selected virus-like contigs representing signal crayfish associated viruses identified in this stud; Alignment S1. Sequence alignments (.fasta) with virus names and corresponding GenBank accession numbers of viruses used for phylogenetic analysis of signal crayfish associated reo-like virus 1; Alignment S2. Sequence alignments (.fasta) with virus names and corresponding GenBank accession numbers of viruses used for phylogenetic analysis of signal crayfish associated hepe-like virus 1; Alignment S3. Sequence alignments (.fasta) with virus names and corresponding GenBank accession numbers of viruses used for phylogenetic analysis of signal crayfish associated toti-like virus 1; Alignment S4. Sequence alignments (.fasta) with virus names and corresponding GenBank accession numbers of viruses used for phylogenetic analysis of signal crayfish associated picorna-like virus 1.

Author Contributions

Conceptualization, S.H. and S.Č.; methodology, D.K., K.B., S.Č. and S.H.; formal analysis, K.B.; investigation, S.H., A.B. and S.Č.; data curation, K.B.; writing, K.B. and S.H.; original draft preparation, K.B. and S.H.; writing—review and editing, K.B., A.B., D.K., S.H. and S.Č.; visualization, K.B., D.K., S.H. and S.Č.; supervision, D.K. and S.H.; project administration, S.H. and D.K.; funding acquisition, S.H. and D.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Croatian Science Foundation grant (HRZZ UIP 2017-05-1720) and Slovenian research agency (research core funding No P4-0407 and program for young researchers awarded to K.B.).

Institutional Review Board Statement

Ethical approval was not required for this study since no institutional, national, or EU ethical guidelines exist for crayfish. However, the dissection procedure (rapid cut of individuals’ nerve cord from the thorax to the end of abdomen) was performed according to available guidelines for the humane killing of crayfish [70].

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequencing data that support the findings of this study are linked under the accession number PRJNA754774 in the NCBI BioProject database (https://www.ncbi.nlm.nih.gov/sra, accessed on 9 November 2021). The sequences of selected viral contigs identified in this study are available in Table S2 and have been deposited in the GenBank under the accession numbers OK317706-OK317734.

Acknowledgments

We thank Paula Dragičević, Ivana Maguire, Dora Pavić, and Leona Lovrenčić for their help in crayfish sampling, and to Paula Dragičević for her help with RNA extraction.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Y.Z.; Shi, M.; Holmes, E.C. Using Metagenomics to Characterize an Expanding Virosphere. Cell 2018, 172, 1168–1172. [Google Scholar] [CrossRef]
  2. Li, C.X.; Shi, M.; Tian, J.H.; Lin, X.D.; Kang, Y.J.; Chen, L.J.; Qin, X.C.; Xu, J.; Holmes, E.C.; Zhang, Y.Z. Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. Elife 2015, 2015, 1–26. [Google Scholar] [CrossRef]
  3. Shi, M.; Lin, X.D.; Tian, J.H.; Chen, L.J.; Chen, X.; Li, C.X.; Qin, X.C.; Li, J.; Cao, J.P.; Eden, J.S.; et al. Redefining the invertebrate RNA virosphere. Nature 2016, 540, 539–543. [Google Scholar] [CrossRef]
  4. Wu, H.; Pang, R.; Cheng, T.; Xue, L.; Zeng, H.; Lei, T.; Chen, M.; Wu, S.; Ding, Y.; Zhang, J.; et al. Abundant and Diverse RNA Viruses in Insects Revealed by RNA-Seq Analysis: Ecological and Evolutionary Implications. Msystems 2020, 5, e00039-20. [Google Scholar] [CrossRef] [PubMed]
  5. Chang, T.; Hirai, J.; Hunt, B.P.V.; Suttle, C.A. Arthropods and the evolution of RNA viruses. bioRxiv 2021. [Google Scholar] [CrossRef]
  6. Dolja, V.V.; Koonin, E.V. Metagenomics reshapes the concepts of RNA virus evolution by revealing extensive horizontal virus transfer. Virus Res. 2018, 244, 36–52. [Google Scholar] [CrossRef]
  7. Zhang, Z.-Q. Phylum Arthropoda. In Animal Biodiversity: An Outline of Higher-Level Classification and Survey of Taxonomic Richness; Zootaxa: Auckland, New Zealand, 2013; Volume 3703, p. 17. [Google Scholar]
  8. Reynolds, J.; Souty-Grosset, C.; Richardson, A. Ecological roles of crayfish in freshwater and terrestrial habitats. Freshw. Crayfish 2013, 19, 197–218. [Google Scholar] [CrossRef]
  9. Holdich, D.M.; Reynolds, J.D.; Souty-Grosset, C.; Sibley, P.J. A review of the ever increasing threat to European crayfish from non-indigenous crayfish species. Knowl. Manag. Aquat. Ecosyst. 2009, 2009, 394–395. [Google Scholar] [CrossRef] [Green Version]
  10. Lodge, D.M.; Deines, A.; Gherardi, F.; Yeo, D.C.J.; Arcella, T.; Baldridge, A.K.; Barnes, M.A.; Lindsay Chadderton, W.; Feder, J.L.; Gantz, C.A.; et al. Global introductions of crayfishes: Evaluating the impact of species invasions on ecosystem services. Annu. Rev. Ecol. Evol. Syst. 2012, 43, 449–472. [Google Scholar] [CrossRef]
  11. Twardochleb, L.A.; Olden, J.D.; Larson, E.R. A global meta-analysis of the ecological impacts of nonnative crayfish. Freshw. Sci. 2013, 32, 1367–1382. [Google Scholar] [CrossRef]
  12. Souty-Grosset, C.; Grandjean, F.; Gouin, N. Conservation and management of native crayfish populations. Freshw. Crayfish 2004, 14, 1–20. [Google Scholar]
  13. Richman, N.I.; Böhm, M.; Adams, S.B.; Alvarez, F.; Bergey, E.A.; Bunn, J.J.S.; Burnham, Q.; Cordeiro, J.; Coughran, J.; Crandall, K.A.; et al. Multiple drivers of decline in the global status of freshwater crayfish (Decapoda: Astacidea). Philos. Trans. R. Soc. B Biol. Sci. 2015, 370, 1–11. [Google Scholar] [CrossRef] [PubMed]
  14. Crowl, T.A.; Crist, T.O.; Parmenter, R.R.; Belovsky, G.; Lugo, A.E. The spread of invasive species and infectious disease as drivers of ecosystem change. Front. Ecol. Environ. 2008, 6, 238–246. [Google Scholar] [CrossRef]
  15. Van der Veer, G.; Nentwig, W. Environmental and economic impact assessment of alien and invasive fish species in Europe using the generic impact scoring system. Ecol. Freshw. Fish 2015, 24, 646–656. [Google Scholar] [CrossRef]
  16. Unestam, T.; Weiss, D.W. The host-parasite relationship between freshwater crayfish and the crayfish disease fungus Aphanomyces astaci: Responses to infection by a susceptible and a resistant species. J. Gen. Microbiol. 1970, 60, 77–90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Martín-Torrijos, L.; Martínez-Ríos, M.; Casabella-Herrero, G.; Adams, S.B.; Jackson, C.R.; Diéguez-Uribeondo, J. Tracing the origin of the crayfish plague pathogen, Aphanomyces astaci, to the Southeastern United States. Sci. Rep. 2021, 11, 1–11. [Google Scholar] [CrossRef]
  18. Zhang, Q.Y.; Gui, J.F. Virus genomes and virus-host interactions in aquaculture animals. Sci. China Life Sci. 2015, 58, 156–169. [Google Scholar] [CrossRef] [Green Version]
  19. Geoghegan, J.L.; Di Giallonardo, F.; Cousins, K.; Shi, M.; Williamson, J.E.; Holmes, E.C. Hidden diversity and evolution of viruses in market fish. Virus Evol. 2018, 4, 1–11. [Google Scholar] [CrossRef]
  20. Dragičević, P.; Bielen, A.; Petrić, I.; Hudina, S. Microbial pathogens of freshwater crayfish: A critical review and systematization of the existing data with directions for future research. J. Fish Dis. 2020, 44, 221–247. [Google Scholar] [CrossRef]
  21. Sánchez-Paz, A. White spot syndrome virus: An overview on an emergent concern. Vet. Res. 2010, 41, 43. [Google Scholar] [CrossRef] [Green Version]
  22. Longshaw, M. Parasites, Commensals, Pathogens and Diseases of Crayfish. In Biology and Ecology of Crayfish; Stebbing, M., Longshaw, P., Eds.; Taylor & Francis Group, LLC: Boca Raton, FL, USA, 2016; pp. 171–250. [Google Scholar]
  23. Bateman, K.S.; Stentiford, G.D. A taxonomic review of viruses infecting crustaceans with an emphasis on wild hosts. J. Invertebr. Pathol. 2017, 147, 86–110. [Google Scholar] [CrossRef]
  24. Grandjean, F.; Gilbert, C.; Razafimafondy, F.; Vucić, M.; Delaunay, C.; Gindre, P.; Bouchard, J.; Raimond, M.; Moumen, B. A new bunya-like virus associated with mass mortality of white-clawed crayfish in the wild. Virology 2019, 533, 115–124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kouba, A.; Petrusek, A.; Kozák, P. Continental-wide distribution of crayfish species in Europe: Update and maps. Knowl. Manag. Aquat. Ecosyst. 2014, 5, 413. [Google Scholar] [CrossRef]
  26. Hudina, S.; Faller, M.; Lucić, A.; Klobučar, G.; Maguire, I. Distribution and dispersal of two invasive crayfish species in the Drava River basin, Croatia. Knowl. Manag. Aquat. Ecosyst. 2009, 9, 394–395. [Google Scholar] [CrossRef] [Green Version]
  27. Dragičević, P.; Faller, M.; Kutleša, P.; Hudina, S. Update on the signal crayfish, pacifastacus leniusculus (Dana, 1852) range expansion in Croatia: A 10-year report. BioInvasions Rec. 2020, 9, 793–807. [Google Scholar] [CrossRef]
  28. Hudina, S.; Žganec, K.; Lucić, A.; Trgovčić, K.; Maguire, I. Recent invasion of the karstic river systems in croatia through illegal introductions of the signal crayfish. Freshw. Crayfish 2013, 19, 21–27. [Google Scholar] [CrossRef]
  29. Hudina, S.; Kutleša, P.; Trgovčić, K.; Duplić, A. Dynamics of range expansion of the signal crayfish (Pacifastacus leniusculus) in a recently invaded region in Croatia. Aquat. Invasions 2017, 12, 67–75. [Google Scholar] [CrossRef]
  30. Anderson, R.M.; May, R.M. Coevolution of Hosts and Parasites. Parasitology 1982, 85, 411–426. [Google Scholar] [CrossRef]
  31. Costa, V.A.; Mifsud, J.C.O.; Gilligan, D.; Williamson, J.E.; Holmes, E.C.; Geoghegan, J.L. Metagenomic sequencing reveals a lack of virus exchange between native and invasive freshwater fish across the Murray-Darling Basin, Australia. Virus Evol. 2021, 7, 1–15. [Google Scholar] [CrossRef] [PubMed]
  32. Westman, K. A new folding trap model which prevents crayfish from escaping. Freshw. Crayfish 1978, 4, 235–242. [Google Scholar]
  33. Buchfink, B.; Xie, C.; Huson, D.H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 2015, 12, 59–60. [Google Scholar] [CrossRef] [PubMed]
  34. Huson, D.H.; Beier, S.; Flade, I.; Gorska, A.; El-Hadidi, M.; Mitra, S.; Ruscheweyh, H.J.; Tappu, R. MEGAN Community Edition—Interactive Exploration and Analysis of Large-Scale Microbiome Sequencing Data. PLoS Comput. Biol. 2016, 12, e1004957. [Google Scholar] [CrossRef] [Green Version]
  35. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. R Core Team. A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. 2021. Available online: https://www.R-project.org (accessed on 9 November 2021).
  37. Ye, J.; Coulouris, G.; Zaretskaya, I.; Cutcutache, I.; Rozen, S.; Madden, T.L. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinform. 2012, 13, 134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  39. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef]
  40. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  41. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [Green Version]
  42. Sánchez, R.; Serra, F.; Tárraga, J.; Medina, I.; Carbonell, J.; Pulido, L.; De María, A.; Capella-Gutíerrez, S.; Huerta-Cepas, J.; Gabaldón, T.; et al. Phylemon 2.0: A suite of web-tools for molecular evolution, phylogenetics, phylogenomics and hypotheses testing. Nucleic Acids Res. 2011, 39, 470–474. [Google Scholar] [CrossRef] [Green Version]
  43. Nguyen, L.T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef]
  44. Nelson, C.W.; Moncla, L.H.; Hughes, A.L. SNPGenie: Estimating evolutionary parameters to detect natural selection using pooled next-generation sequencing data. Bioinformatics 2015, 31, 3709–3711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wolf, Y.I.; Silas, S.; Wang, Y.; Wu, S.; Bocek, M.; Kazlauskas, D.; Krupovic, M.; Fire, A.; Dolja, V.V.; Koonin, E.V. Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome. Nat. Microbiol. 2020, 5, 1262–1270. [Google Scholar] [CrossRef] [PubMed]
  46. Xu, C.; Sun, X.; Taylor, A.; Jiao, C.; Xu, Y.; Cai, X.; Wang, X.; Ge, C.; Pan, G.; Wang, Q.; et al. Diversity, Distribution, and Evolution of Tomato Viruses in China Uncovered by Small RNA Sequencing. J. Virol. 2017, 91, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Hayakijkosol, O.; Jaroenram, W.; Owens, L.; Elliman, J. Reverse transcription polymerase chain reaction (RT-PCR) detection for Australian Cherax reovirus from redclaw crayfish (Cherax quadricarinatus). Aquaculture 2021, 530, 735881. [Google Scholar] [CrossRef]
  48. Hayakijkosol, O.; Owens, L. Investigation into the pathogenicity of reovirus to juvenile Cherax quadricarinatus. Aquaculture 2011, 316, 1–5. [Google Scholar] [CrossRef]
  49. Bekavac, A.; Dragičević, P.; Dragun, Z.; Beck, A.; Maguire, I.; Ivanković, D.; Fiket, Ž.; Gračan, R.; Hudina, S. Disturbance in invasion? Idiopathic necrotizing hepatopancreatitis in the signal crayfish Pacifastacus leniusculus (Dana, 1852) in Croatia. J. Fish Dis. 2021, 1–16. [Google Scholar] [CrossRef]
  50. Dong, X.; Hu, T.; Liu, Q.; Li, C.; Sun, Y.; Wang, Y.; Shi, W.; Zhao, Q.; Huang, J. A Novel Hepe-Like Virus from Farmed Giant Freshwater Prawn Macrobrachium rosenbergii. Viruses 2020, 12, 323. [Google Scholar] [CrossRef] [Green Version]
  51. Lightner, D.V.; Pantoja, C.R.; Poulos, B.T.; Tang, K.F.J.; Redman, R.M.; Andrade, T.P.; Bonami, J.R. University Infectious Myonecrosis- New Disease in Pacific White Shrimp. Glob. Aquac. Advocate 2004, 7, 85. [Google Scholar]
  52. Edgerton, B.; Owens, L.; Glasson, B.; De Beer, S. Description of a small dsRNA virus from freshwater crayfish Cherax quadricarinatus. Dis. Aquat. Organ. 1994, 18, 63–69. [Google Scholar] [CrossRef]
  53. Sunarto, A.; Naim, S. Totiviruses of Crustaceans. In Aquaculture Virology; Elsevier Inc.: Amsterdam, The Netherlands, 2016; pp. 425–439. [Google Scholar]
  54. Zeng, D.; Chen, X.; Xie, D.; Zhao, Y.; Yang, C.; Li, Y.; Ma, N.; Peng, M.; Yang, Q.; Liao, Z.; et al. Transcriptome Analysis of Pacific White Shrimp (Litopenaeus vannamei) Hepatopancreas in Response to Taura Syndrome Virus (TSV) Experimental Infection. PLoS ONE 2013, 8, e57515. [Google Scholar] [CrossRef]
  55. Guo, Z.X.; He, J.G.; Xu, H.D.; Weng, S.P. Pathogenicity and complete genome sequence analysis of the mud crab dicistrovirus-1. Virus Res. 2013, 171, 8–14. [Google Scholar] [CrossRef] [PubMed]
  56. Gao, Y.; Liu, S.; Huang, J.; Wang, Q.; Li, K.; He, J.; He, J.; Weng, S.; Zhang, Q. Cryo-electron Microscopy Structures of Novel Viruses from Mud Crab Scylla paramamosain with Multiple Infections. J. Virol. 2019, 93, e02255-18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Waldron, F.M.; Stone, G.N.; Obbard, D.J. Metagenomic sequencing suggests a diversity of RNA interference-like responses to viruses across multicellular eukaryotes. PLoS Genet. 2020, 16, e1007533. [Google Scholar] [CrossRef]
  58. Argenta, F.F.; Hepojoki, J.; Smura, T.; Szirovicza, L.; Hammerschmitt, M.E.; Driemeier, D.; Kipar, A.; Hetzel, U. Identification of Reptarenaviruses, Hartmaniviruses, and a Novel Chuvirus in Captive Native Brazilian Boa Constrictors with Boid Inclusion Body Disease. J. Virol. 2020, 94, e00001-20. [Google Scholar] [CrossRef] [PubMed]
  59. Dezordi, F.Z.; dos Santos Vasconcelos, C.R.; Rezende, A.M.; Wallau, G.L. In and Outs of Chuviridae Endogenous Viral Elements: Origin of a Potentially New Retrovirus and Signature of Ancient and Ongoing Arms Race in Mosquito Genomes. Front. Genet. 2020, 11, 1291. [Google Scholar] [CrossRef] [PubMed]
  60. Anderson, R.M.; May, R.M. Population biology of infectious diseases: Part I. Nature 1979, 280, 361–367. [Google Scholar] [CrossRef]
  61. Heesterbeek, J.A.P. A brief history of R0 and a recipe for its calculation. Acta Biotheor. 2002, 50, 189–204. [Google Scholar] [CrossRef]
  62. Lloyd-Smith, J.O.; Schreiber, S.J.; Kopp, P.E.; Getz, W.M. Superspreading and the effect of individual variation on disease emergence. Nature 2005, 438, 355–359. [Google Scholar] [CrossRef]
  63. Lindenfors, P.; Nunn, C.L.; Jones, K.E.; Cunningham, A.A.; Sechrest, W.; Gittleman, J.L. Parasite species richness in carnivores: Effects of host body mass, latitude, geographical range and population density. Glob. Ecol. Biogeogr. 2007, 16, 496–509. [Google Scholar] [CrossRef]
  64. Gay, N.; Olival, K.J.; Bumrungsri, S.; Siriaroonrat, B.; Bourgarel, M.; Morand, S. Parasite and viral species richness of Southeast Asian bats: Fragmentation of area distribution matters. Int. J. Parasitol. Parasites Wildl. 2014, 3, 161–170. [Google Scholar] [CrossRef] [Green Version]
  65. Webber, Q.M.R.; Fletcher, Q.E.; Willis, C.K.R. Viral Richness is Positively Related to Group Size, but Not Mating System, in Bats. Ecohealth 2017, 14, 652–661. [Google Scholar] [CrossRef]
  66. Rose, R.; Constantinides, B.; Tapinos, A.; Robertson, D.L.; Prosperi, M. Challenges in the analysis of viral metagenomes. Virus Evol. 2016, 2, vew022. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Pedersen, A.B.; Davies, T.J. Cross-species pathogen transmission and disease emergence in primates. Ecohealth 2009, 6, 496–508. [Google Scholar] [CrossRef] [PubMed]
  68. Zhao, M.; dos Santos Tavares, C.P.; Schott, E.J. Diversity and classification of reoviruses in crustaceans: A proposal. J. Invertebr. Pathol. 2021, 182, 107568. [Google Scholar] [CrossRef] [PubMed]
  69. Cote, J.; Brodin, T.; Fogarty, S.; Sih, A. Non-random dispersal mediates invader impacts on the invertebrate community. J. Anim. Ecol. 2017, 86, 1298–1307. [Google Scholar] [CrossRef] [Green Version]
  70. Conte, F.; Voslarova, E.; Vecerek, V.; Elwood, R.W.; Coluccio, P.; Pugliese, M.; Passantino, A. Humane slaughter of edible decapod crustaceans. Animals 2021, 11, 1089. [Google Scholar] [CrossRef]
Viruses 13 02259 g001 550
Figure 1. Position of sampling locations and differences in crayfish (Pacifastacus leniusculus and Pontastacus leptodactylus) presence and abundance (CPUE; catch per unit effort) along the signal crayfish invasion range in the Korana River in 2018. Sampling was performed at both upstream (UF) and downstream (DF) invasion fronts and upstream (UC) and downstream (DC) invasion cores.
Figure 1. Position of sampling locations and differences in crayfish (Pacifastacus leniusculus and Pontastacus leptodactylus) presence and abundance (CPUE; catch per unit effort) along the signal crayfish invasion range in the Korana River in 2018. Sampling was performed at both upstream (UF) and downstream (DF) invasion fronts and upstream (UC) and downstream (DC) invasion cores.
Viruses 13 02259 g001
Viruses 13 02259 g002 550
Figure 2. Heatmap showing average read coverage resulting from mapping of normalized read datasets from individual locations to selected putative invertebrate virus-like contigs together with the dendrograms clustering the sampling locations (top) and viral contigs (left) according to average read coverage values. Heatmap and dendrograms were constructed using decimal logarithm transformed values, however, non-transformed average read coverage numbers are plotted on the heatmap.
Figure 2. Heatmap showing average read coverage resulting from mapping of normalized read datasets from individual locations to selected putative invertebrate virus-like contigs together with the dendrograms clustering the sampling locations (top) and viral contigs (left) according to average read coverage values. Heatmap and dendrograms were constructed using decimal logarithm transformed values, however, non-transformed average read coverage numbers are plotted on the heatmap.
Viruses 13 02259 g002
Viruses 13 02259 g003 550
Figure 3. Gel electrophoresis of PCR amplicons originating from newly identified signal crayfish-associated viruses. Four primer pairs were used for amplifying app. 200 nts long regions of signal crayfish-associated reo-like virus 1 (V1), signal crayfish-associated hepe-like virus 1 (V2), signal crayfish-associated toti-like virus 1 (V3), and signal crayfish-associated picorna-like virus 1 (V4) genomes. Samples from different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front) were analyzed. Unmarked wells represent no template controls of each amplification reaction.
Figure 3. Gel electrophoresis of PCR amplicons originating from newly identified signal crayfish-associated viruses. Four primer pairs were used for amplifying app. 200 nts long regions of signal crayfish-associated reo-like virus 1 (V1), signal crayfish-associated hepe-like virus 1 (V2), signal crayfish-associated toti-like virus 1 (V3), and signal crayfish-associated picorna-like virus 1 (V4) genomes. Samples from different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front) were analyzed. Unmarked wells represent no template controls of each amplification reaction.
Viruses 13 02259 g003
Viruses 13 02259 g004 550
Figure 4. Phylogenetic relationships and genome organization of signal crayfish-associated reo-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated reo-like virus 1 (virus-like contig 4) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site (B) The predicted partial genome organization of the novel signal crayfish-associated reo-like virus 1, with positions and length of the open reading frames (ORF) indicated with corresponding arrow length.
Figure 4. Phylogenetic relationships and genome organization of signal crayfish-associated reo-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated reo-like virus 1 (virus-like contig 4) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site (B) The predicted partial genome organization of the novel signal crayfish-associated reo-like virus 1, with positions and length of the open reading frames (ORF) indicated with corresponding arrow length.
Viruses 13 02259 g004
Viruses 13 02259 g005 550
Figure 5. Phylogenetic relationships and genome organization of signal crayfish-associated hepe-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated hepe-like virus 1 (virus-like contig 139) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site. (B) The predicted genome organization of the novel signal crayfish-associated hepe-like virus 1, with positions and lengths of the open reading frames (ORF) indicated with corresponding arrow lengths.
Figure 5. Phylogenetic relationships and genome organization of signal crayfish-associated hepe-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated hepe-like virus 1 (virus-like contig 139) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site. (B) The predicted genome organization of the novel signal crayfish-associated hepe-like virus 1, with positions and lengths of the open reading frames (ORF) indicated with corresponding arrow lengths.
Viruses 13 02259 g005
Viruses 13 02259 g006 550
Figure 6. Phylogenetic relationships and genome organization of signal crayfish-associated toti-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated toti-like virus 1 (virus-like contig 141) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site. (B) The predicted genome organization of the novel signal crayfish-associated toti-like virus 1 with positions and lengths of the open reading frames (ORF) indicated with corresponding arrow lengths.
Figure 6. Phylogenetic relationships and genome organization of signal crayfish-associated toti-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated toti-like virus 1 (virus-like contig 141) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site. (B) The predicted genome organization of the novel signal crayfish-associated toti-like virus 1 with positions and lengths of the open reading frames (ORF) indicated with corresponding arrow lengths.
Viruses 13 02259 g006
Viruses 13 02259 g007 550
Figure 7. Phylogenetic relationships and genome organization of signal crayfish-associated picorna-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated picorna-like virus 1 (virus-like contig 1) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site. (B) The predicted partial genome organization of the novel signal crayfish-associated picorna-like virus 1, with positions and length of the open reading frame (ORF) indicated with corresponding arrow length.
Figure 7. Phylogenetic relationships and genome organization of signal crayfish-associated picorna-like virus 1. (A) Phylogenetic tree built with maximum likelihood approach, based on the alignment of the conserved segment of RNA-dependent RNA polymerase domain of signal crayfish-associated picorna-like virus 1 (virus-like contig 1) and representative selected sequences of phylogenetically related viruses. The numbers on the branches represent bootstrap support values (>50% shown), the branch length represents the average number of amino acid substitutions per site. (B) The predicted partial genome organization of the novel signal crayfish-associated picorna-like virus 1, with positions and length of the open reading frame (ORF) indicated with corresponding arrow length.
Viruses 13 02259 g007
Viruses 13 02259 g008 550
Figure 8. Nucleotide diversity values (π) calculated for signal crayfish-associated hepe-like virus 1, signal crayfish-associated toti-like virus 1, and signal crayfish-associated reo-like virus 1 at different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front) along the invasion range of the signal crayfish in the Korana River, Croatia.
Figure 8. Nucleotide diversity values (π) calculated for signal crayfish-associated hepe-like virus 1, signal crayfish-associated toti-like virus 1, and signal crayfish-associated reo-like virus 1 at different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front) along the invasion range of the signal crayfish in the Korana River, Croatia.
Viruses 13 02259 g008
Table
Table 1. Taxonomic classification of reads obtained from different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front).
Table 1. Taxonomic classification of reads obtained from different locations (UF—upstream front, UC—upstream core, DC—downstream core, DF—downstream front).
Sample NameRaw ReadsReads after TrimmingAccession Number (SRA *)% of Classified Reads (Diamond)% of Reads Class. as Viral (Diamond)% of Reads Class. as Riboviria (Diamond)No. of Viral Read Class. as Riboviria (Diamond)
UF164,043,698160,643,914SAMN2080094115.470.001590.000941511
UC149,357,758146,166,391SAMN2080094219.060.008950.0081611,929
DC100,131,42097,563,944SAMN2080094323.240.017370.0162015,802
DF136,134,398132,904,626SAMN2080094415.290.009610.0082110,905
* short reads archive (SRA).
Table
Table 2. Selected putative invertebrate virus-like contigs (>300 nts) representing signal crayfish-associated viruses identified in this study with their contig length (nts) and closest match in the blastx (NCBI-nr) analysis.
Table 2. Selected putative invertebrate virus-like contigs (>300 nts) representing signal crayfish-associated viruses identified in this study with their contig length (nts) and closest match in the blastx (NCBI-nr) analysis.
Contig NameVirus NameGenBank AccesionContig LengthBlastx Results
Closest Protein Hit Name and Speciese-ValueIdentities *
Virus-like contig 4Signal crayfish associated reo-like virus 1OK3177064234RdRp [Cherax quadricarinatus reovirus]2.74 × 10−20106/422 (25%)
Virus-like contig 139Signal crayfish associated hepe-like virus 1OK31770710,400hypothetical protein [Beihai hepe-like virus 4]2.64 × 10−117213/480 (44%)
Virus-like contig 141Signal crayfish associated toti-like virus 1OK3177088576hypothetical protein 4 [Wenzhou crab virus 5]<1.00 × 10−250353/869 (41%)
Virus-like contig 1Signal crayfish associated picorna-like virus 1OK3177114587hypothetical protein [Beihai picorna-like virus 99]1.28 × 10−86324/1251 (26%)
Virus-like contig 10Signal crayfish associated picorna-like virus 2OK317712531hypothetical protein [Beihai picorna-like virus 99]2.75 × 10−3979/175 (45%)
Virus-like contig 11Signal crayfish associated picorna-like virus 3OK317713320hypothetical protein [Beihai picorna-like virus 99]4.97 × 10−5180/105 (76%)
Virus-like contig 15Signal crayfish associated picorna-like virus 4OK317714898hypothetical protein 1 [Changjiang picorna-like virus 6]1.67 × 10−3785/151 (56%)
Virus-like contig 9Signal crayfish associated picorna-like virus 5OK317715771hypothetical protein [Beihai picorna-like virus 99] Y1.30 × 10−96170/257 (66%)
Virus-like contig 5Signal crayfish associated picorna-like virus 6OK317716714hypothetical protein [Wenzhou picorna-like virus 38]1.56 × 10−1962/152 (41%)
Virus-like contig 13Signal crayfish associated picorna-like virus 7OK3177171894hypothetical protein [Beihai sesarmid crab virus 2]9.08 × 10−99226/642 (35%)
Virus-like contig 66Signal crayfish associated tombus-like virus 1OK3177184504replicase [Caledonia beadlet anemone tombus-like virus 1]<1.00 × 10−250352/862 (41%)
Virus-like contig 84Signal crayfish associated tombus-like virus 2OK3177192981replicase [Caledonia beadlet anemone tombus-like virus 1]3.00 × 10−114246/628 (39%)
Virus-like contig 55Signal crayfish associated tombus-like virus 3OK3177201425RdRp [Riboviria sp.]2.00 × 10−92150/299 (50%)
Virus-like contig 35Signal crayfish associated tombus-like virus 4OK317721665RdRp [Riboviria sp.]4.00 × 10−3276/161 (47%)
Virus-like contig 24Signal crayfish associated tombus-like virus 5OK317722538hypothetical protein 1 [Hubei tombus-like virus 16]3.98 × 10−0443/99 (43%)
Virus-like contig 169Signal crayfish associated tombus-like virus 6OK317723301hypothetical protein 2 [Hubei unio douglasiae virus 2]2.00 × 10−4269/82 (84%)
Virus-like contig 140Signal crayfish associated chu-like virus 1OK3177242216hypothetical protein 2 [Beihai hermit crab virus 3]5.71 × 10−56161/588 (27%)
Virus-like contig 65Signal crayfish associated chu-like virus 2OK317725746hypothetical protein 2 [Beihai hermit crab virus 3]6.45 × 10−3569/167 (41%)
Virus-like contig 7Signal crayfish associated chu-like virus 3OK317726493RdRp [Beihai hermit crab virus 3]3.98 × 10−4844/54 (81%)
Virus-like contig 145Signal crayfish associated chu-like virus 4OK3177271009RdRp [Beihai hermit crab virus 3]1.63E × 10−60120/219 (55%)
Virus-like contig 146Signal crayfish associated chu-like virus 5OK317728418RdRp [Beihai hermit crab virus 3]1.00E × 10−2754/98 (55%)
Virus-like contig 222Signal crayfish associated chu-like virus 6OK317729356RdRp [Beihai hermit crab virus 3]2.00E × 10−4580/111 (72%)
Virus-like contig 27Signal crayfish associated partiti-like virus 1OK3177301188RdRp [Caledonia partiti-like virus]1.00E × 10−142221/394 (56%)
Virus-like contig 116Signal crayfish associated picorna-like virus 8OK317731734polyprotein [Picornaviridae sp.]3.70E × 10−1042/88 (48%)
Virus-like contig 2Signal crayfish associated picorna-like virus 9OK317732412hypothetical protein 1 [Picornavirales sp.]1.05E × 10−2445/100 (45%)
Virus-like contig 36Signal crayfish associated picorna-like virus 10OK317733457RdRp [Picornavirales sp.]4.32 × 10−1660/102 (58%)
Virus-like contig 83Signal crayfish associated sobemo-like virus 1OK317734431hypothetical protein [Hubei sobemo-like virus 43]1.00 × 10−5084/143 (59%)
Virus-like contig 147Signal crayfish associated sobemo-like virus 2OK317709512hypothetical protein 1 [Beihai sobemo-like virus 17]1.00E × 10−1044/128 (34%)
Virus-like contig 30Signal crayfish associated narna-like virus 1OK317710378RdRp [Beihai narna-like virus 18]5.60 × 10−1848/128 (38%)
* number of matching amino acids/number of total amino acids of the closest protein hit (percentage of amino acid identity).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bačnik, K.; Kutnjak, D.; Černi, S.; Bielen, A.; Hudina, S. Virome Analysis of Signal Crayfish (Pacifastacus leniusculus) along Its Invasion Range Reveals Diverse and Divergent RNA Viruses. Viruses 2021, 13, 2259. https://doi.org/10.3390/v13112259

AMA Style

Bačnik K, Kutnjak D, Černi S, Bielen A, Hudina S. Virome Analysis of Signal Crayfish (Pacifastacus leniusculus) along Its Invasion Range Reveals Diverse and Divergent RNA Viruses. Viruses. 2021; 13(11):2259. https://doi.org/10.3390/v13112259

Chicago/Turabian Style

Bačnik, Katarina, Denis Kutnjak, Silvija Černi, Ana Bielen, and Sandra Hudina. 2021. "Virome Analysis of Signal Crayfish (Pacifastacus leniusculus) along Its Invasion Range Reveals Diverse and Divergent RNA Viruses" Viruses 13, no. 11: 2259. https://doi.org/10.3390/v13112259

APA Style

Bačnik, K., Kutnjak, D., Černi, S., Bielen, A., & Hudina, S. (2021). Virome Analysis of Signal Crayfish (Pacifastacus leniusculus) along Its Invasion Range Reveals Diverse and Divergent RNA Viruses. Viruses, 13(11), 2259. https://doi.org/10.3390/v13112259

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop