Characterization of a Novel Emaravirus Affecting Ash Species ( Fraxinus spp.) in Europe

: We identiﬁed a novel virus in diseased European ash ( Fraxinus excelsior) and manna ash ( F. ornus ) trees exhibiting chlorotic ringspots, mottle and leaf deformation such as curling and shoestring symptoms. High-throughput sequencing (HTS, Illumina RNASeq) of total RNA isolated from diseased leaf material in combination with RT-PCR-based ampliﬁcation techniques and Sanger sequencing determined ﬁve complete genome segments, each encoding a single open reading frame. Sequence analyses of RNA1–RNA5 revealed a genome organization typical for emaraviruses, i.e., (i) conserved and complementary terminal 5 (cid:48) and 3 (cid:48) termini of each genome segment (ii) proteins showing signiﬁcant homologies to the RNA-dependent RNA polymerase (RdRP) encoded by RNA1, the glycoprotein precursor (GPP) encoded by RNA2, the viral nucleocapsid protein (N, RNA3), the movement protein (MP, RNA4), and a protein of 26 kDA (P26, RNA5) highly similar to proteins of unknown function encoded by other emaraviruses. Furthermore, we identiﬁed spherical particles (double-membrane bodies, DMB) of different sizes (70–80 nm in diameter) which are typical for emaraviruses exclusively in virus-infected leaf tissue exhibiting mottle and leaf deformation. Sequence comparison and phylogenetic analyses conﬁrmed the identiﬁed novel virus as a new member of the genus Emaravirus . We established a species-speciﬁc RT-PCR detection protocol and could associate the observed disease symptoms with the infection of the novel emaravirus in F. excelsior and F. ornus . Therefore, we propose the name ash shoestring-associated emaravirus (ASaV). Investigation of ASaV-infected sample trees originating from different locations in Switzerland, Germany, Italy and Sweden provided a wide geographical distribution of the virus in affected ash species. To our knowledge, this is the ﬁrst conﬁrmation of an emaravirus affecting ash tree species with shoestring symptoms of leaves in Europe.


Introduction
Although virus-like symptoms have long been observed in ash [1,2], no viral agents affecting Fraxinus spp. have to date been studied in detail. The list of plant viruses known to infect ash comprises of generalist species including arabis mosaic virus (ArMV), cherry leaf roll virus (CLRV), tobacco necrosis virus (TNV), tobacco mosaic virus (TMV), tobacco rattle virus (TRV), tomato ringspot virus (ToRSV) and tobacco ringspot virus (TRSV) [3,4]. Our investigations of Fraxinus spp. were instigated by current interest in ash dieback, a disease that has decimated populations of European ash (Fraxinus excelsior L.) across the continent in recent years [5]. Although the disease is primarily thought to be caused by the pathogenic fungus Hymenoscyphus fraxineus, the role of other biotic and abiotic stress factors, including viruses, in disease development remains to be investigated.
Ash (Fraxinus spp.; Family: Oleaceae) is a group of broad-leaved, mostly deciduous trees, with both the European ash (F. excelsior) and the manna ash (F. ornus) being native to Europe. Not only do both species have important roles in primary or secondary succession stages in forest ecosystems, they also have a decorative, ornamental use in urban environments [6][7][8]. Manna ash trees, which are naturally found in the Mediterranean region, are being increasingly trialed in towns and cities due to their drought tolerance in the face of climate change [9,10]. Manna ash trees also appear to be less susceptible to ash dieback [11].
The causal agents of virus-suspect symptoms in forestry and urban green spaces were for a long time unknown, their identification only being made possible with advances in molecular biology over the last decades. Although the focus of much research on plant viral pathogens has concerned important agricultural and horticultural crops including fruit trees [3], the increasing application of high throughput sequencing (HTS) technologies is revolutionizing the diagnosis and characterization of plant viruses [12,13]. It is clear that a large uncharacterized forest tree virome is yet to be discovered [14]. HTS technologies have already been successfully used in the identification of a number of novel viruses in trees [15] and also in agricultural crops [16]. For instance, the virome of leaf roll diseased birch trees (Betula spp.) was studied revealing several novel viruses belonging to different genera [17,18]. Likewise, novel viruses of the genus Emaravirus were discovered in pedunculate oak (Quercus robur) [19,20] in sycamore maple (Acer pseudoplatanus) [21] and in Eurasian aspen (Populus tremula) [22] by Illumina RNASeq technology, applying total RNA extracted from diseased leaf material.
The genus Emaravirus (Order: Bunyavirales; Family: Fimoviridae) is comprised of viruses with a single-stranded (ss), segmented RNA genome in negative-sense orientation, which are typically transmitted by eriophyid gall mites [23,24]. Other characteristics of emaraviruses include conserved, complementary terminal nucleotide sequences, enabling the circularization of genome segments and the formation of a panhandle structure [25], which is thought to be necessary for transcription and replication [26]. Spherical double membrane-bound bodies (DMBs) found in infected cells of host plants are also thought to represent the enveloped virus particles [24]. Since the description of the first member species European mountain ash ringspot-associated emaravirus (EMARaV) and the establishment of the genus Emaravirus [24,26], many other viruses have been assigned to the genus, currently numbering 21 [27]. Most emaraviruses have a narrow host range. Nevertheless, many of them cause economic losses in their respective natural host including important crops, or they could be associated with characteristic disease symptoms in herbaceous plants as well as in deciduous tree species [24]. However, prior to this study, no emaravirus species has been identified infecting Fraxinus spp.
We chose to apply HTS for the identification of viruses associated with observed symptoms in diseased Fraxinus trees. These investigations were complemented by standard, conventional methods of RT-PCR, cloning and Sanger sequencing to obtain the complete sequence of the identified viruses, to develop a reliable detection method and to associate the observed disease symptoms with the novel ash pathogen. Here, we report on the identification of a novel emaravirus in Fraxinus spp. and provide the first information towards an in-depth biological characterization of the virus including genome structure, phylogenetic analyses and a first assessment of genetic variability. Such information not only contributes to our knowledge on emaraviruses, but is also essential for determining the significance of the virus, evaluating ecological and economic impacts and developing possible management strategies in the future.

Plant Material and RNA Extraction
Following observations of suspected virus symptoms in F. excelsior in Switzerland, other symptomatic trees, including manna ash (Fraxinus ornus L.), were found at different geographical locations across Europe. These included trees in both managed semi-natural ecosystems and urban environments, independent of latitude and elevation. Affected ash exhibited chlorotic ringspots, mottle and leaf deformations, ranging from curling to a severe narrowing of the leaflets termed shoestring. Cauliflower galls-growths caused by the gall mite Aceria fraxinivora were also recorded on some of the symptomatic trees. A selection of symptoms documented in F. excelsior from Basel, Switzerland in 2016 and F. ornus from Hamburg, Germany in 2018 are shown in Figure 1.
x FOR PEER REVIEW 3 of 21

Plant Material and RNA Extraction
Following observations of suspected virus symptoms in F. excelsior in Switzerland, other symptomatic trees, including manna ash (Fraxinus ornus L.), were found at different geographical locations across Europe. These included trees in both managed semi-natural ecosystems and urban environments, independent of latitude and elevation. Affected ash exhibited chlorotic ringspots, mottle and leaf deformations, ranging from curling to a severe narrowing of the leaflets termed shoestring. Cauliflower galls-growths caused by the gall mite Aceria fraxinivora were also recorded on some of the symptomatic trees. A selection of symptoms documented in F. excelsior from Basel, Switzerland in 2016 and F. ornus from Hamburg, Germany in 2018 are shown in Figure 1. For HTS analyses, total RNA was extracted as described by Mielke and Mühlbach 2007 [26], excluding the phenol/chloroform step, from leaf material collected in 2016 from two Fraxinus excelsior trees (E55262 and E55268, Basel, Switzerland), displaying virus-like disease symptoms thought to be caused by virus infection. Leaf samples displayed chlorotic ringspots, lines and leaf deformations including a severe narrowing of the leaflets, termed shoestring. RNA was also extracted from a second set of pooled leaf samples (E57321-E57325, E57328-E57330, E57335-E57337, E57339) originating from 12 Fraxinus ornus trees (Hamburg, Germany, sampled in 2018) displaying identical symptoms as the ash trees in Basel. The same total RNA extraction method was also applied to the following samples collected from  For HTS analyses, total RNA was extracted as described by Mielke and Mühlbach 2007 [26], excluding the phenol/chloroform step, from leaf material collected in 2016 from two Fraxinus excelsior trees (E55262 and E55268, Basel, Switzerland), displaying viruslike disease symptoms thought to be caused by virus infection. Leaf samples displayed chlorotic ringspots, lines and leaf deformations including a severe narrowing of the leaflets, termed shoestring. RNA was also extracted from a second set of pooled leaf samples (E57321-E57325, E57328-E57330, E57335-E57337, E57339) originating from 12 Fraxinus ornus trees (Hamburg, Germany, sampled in 2018) displaying identical symptoms as the ash trees in Basel. The same total RNA extraction method was also applied to the following samples collected from locations across Europe between 2016  Hamburg, Germany, sampled 2019). Samples were collected from trees in urban locations, with both symptomatic and asymptomatic leaves being collected for analysis. Eleven F. ornus samples (E57504-E57514) were also collected in 2018 in Hamburg from symptomatic and asymptomatic trees at a local tree nursery and the same methodology was applied.

HTS of Leaf Material from Diseased Fraxinus spp.
Equimolar ratios of total RNA extracted from F. excelsior samples E55262 and E55268 were pooled and 60 µg was subjected to cDNA-synthesis with random hexamer primer (Biolegio B.V., Nijmegen, The Netherlands) as described previously by Bandte et al. [20] and von Bargen et al. [28]. The second sample pool was generated in the same way using total RNA extracted from leaf samples of 12 F. ornus trees (E57321-E57325, E57328-E57330, E57335-E57337, E57339). Both sample pools were further prepared for high-throughput Illumina RNA sequencing (HTS) as outlined in the paper by von Bargen et al. [28]. A standard paired-end 125-bp library (Nextera XT) was prepared from 140 ng of ds cDNA generated from the F. excelsior and 410 ng of ds cDNA of the F. ornus sample pool, respectively, by the sequencing company BaseClear B.V. (Leiden, The Netherlands) and sequenced on a HiSeq2500 machine (Illumina, Illumina Solutions, Berlin, Germany). Quality trimming of raw data and de novo assembly into scaffolds was performed according to an in-house protocol of Baseclear. This protocol included generation of FASTQ sequence files with the Illumina Casava pipeline version 1.8.3 (Illumina Solutions, Berlin, Germany). Reads containing PhiX control signal were removed and reads containing partial adapter sequences up to a minimum read length of 50 bp were clipped. De novo assembly was carried out with CLC Genomics Workbench version 9.5.1 (Qiagen GmbH, Hilden, Germany) after trimming off low-quality bases of FASTQ sequences. Misassembles and nucleotide disagreements between Illumina generated data and the contig sequences were corrected with Pilon version 1.11 [29]. Scaffolds with a minimum length of 300 nucleotides were generated with the SSPACE Premium scaffolder version 2.3 [30] and gapped regions were closed when possible, using GapFiller version 1.10 [31].
Viral sequences in the datasets were identified by blastx search of the assembled sequence scaffolds against the non-redundant protein sequence database provided by the National Center for Bioinformatic Information (NCBI) as described previously [22]. All scaffolds relating to plant viruses were aligned to reference genomes in BioEdit, version 7.0.5.3 [32] using the integrated Clustal-W option.

Determination of the Complete Genome of the Novel Virus
Preliminary analysis of HTS sequence data identified scaffolds relating to five different genomic RNAs of virus species belonging to the genus Emaravirus. In order to determine the complete sequence of the novel virus, full-length RT-PCR was performed, applying the PDAP213 primer [33] using total RNA prepared from F. excelsior leaves (E55270) displaying the shoestring symptom. Full-length PCR products were cloned and sequenced as described in von Bargen et al. [28] and Rehanek et al. [19] with at least two identical fragments being required for sequence confirmation. As amplification of the complete RNA1 could not be achieved, eight primer pairs producing overlapping fragments spanning the complete genome segment were constructed based on original HTS scaffolds (Table S1, Supplementary Information). PCR products of the viral RNA1 of sample E55270 were generated by RT-PCR. Following cDNA synthesis from 1-3 µg of total RNA using random hexamer oligonucleotides (Biolegio GmbH, Nijmegen, The Netherlands), RT-PCR was performed in a 50 µL reaction mixture of 37.5 µL MilliQ-H2O, 10 µL 5 × HiFi Reaction buffer, 1 µL dNTPs (10 mM each), 0.5 µL forward and reverse primers (50 µM), 1 unit (U) Cycling conditions were as follows: initial denaturation at 94 • C for 2 min, followed by 30 to 35 cycles of denaturation at 94 • C for 20 s, annealing at 55 • C for 20 s and elongation at 72 • C for 15-30 s (adjusted for fragment length according to manufactures instructions). Following their purification using the MSB Spin PCRapace kit (Invitek Molecular, Berlin, Germany), generated products were sequenced from both directions with fragment-specific primers (Macrogen Europe B.V., Amsterdam, The Netherlands).
A modified RT-PCR protocol, a tailing PCR, was employed to confirm terminal ends of the RNA1 and 2. Following cDNA synthesis as described above, enzymatic degradation of residual RNA was achieved by adding 3 units (U) of RNase H (Thermo Fisher Scientific, Waltham, MA, USA) to 40 µL of the cDNA and incubating the samples at 37 • C for 30 min. A poly-A (adenine) tail was added to 3 cDNA ends with the help of 1 U of terminal transferase (Thermo Fisher Scientific, Waltham, MA, USA). The poly-A cDNA was then purified (MSB Spin PCRapace kit, Invitek Molecular, Berlin, Germany) and used as a template for PCR, utilizing two virus-specific primers in the form of a semi-nested PCR, allowing the amplification of a specific terminal end of a particular genome segment. Primer sets and PCR conditions for the tailing PCR with velocity polymerase are provided in the Supplementary Information (Table S1). PCR products of expected size were cloned in E. coli (XL1blue, MRF , Thermo Fisher Scientific, Waltham, MA, USA) using the CloneJET PCR cloning kit (Thermo Fisher Scientific, Waltham, MA, USA) and 2-5 clones were sequenced (Macrogen Europe B.V., Amsterdam, The Netherlands) to determine the terminal ends of viral RNA1 and RNA2.
All HTS-and RT-PCR-generated sequence data were edited, aligned and translated in BioEdit, version 7.0.5.3 [32]. The single open reading frame (ORF) of each RNA segment was located and translated using the NCBI ORF finder (https://www.ncbi.nlm.nih.gov/ orffinder (accessed on 21 March 2021)). A comparison of the RNA1-5 encoded proteins of the novel virus with corresponding proteins of classified and putative members of the genus Emaravirus was performed using the Sequence Identity Matrix function of BioEdit [32]. Phylogenetic analyses were performed using MUSCLE [34] to align the RNA1 and RNA3 encoded amino acid sequences of the novel virus identified in F. excelsior (E55270) and F. ornus (HTS data) with corresponding sequences of selected putative and classified emaraviruses. Phylogenetic trees were constructed using neighbor-joining methodology in Geneious prime 2019.1.1 (Biomatters Ltd., Auckland, New Zealand).

Detection of the Novel Virus and Association with Disease Symptoms in Fraxinus spp.
Six specific primer pairs (Table 1) targeting the five genome segments of the novel virus discovered in diseased ash were designed using HTS sequence information from the original pooled sample (E55262 and E55268). The specificity of each segment-specific primer set was evaluated using total RNA extracted from symptomatic and asymptomatic F. excelsior leaves and leaf material of four other deciduous tree species, known to be infected with putative emaraviruses: common oak ringspot-associated virus (CORaV)-infected Quercus robur [19,20], maple mottle-associated virus (MaMaV)-infected Acer pseudoplatanus [21], fig mosaic virus (FMV)-infected Ficus carica [35] and aspen mosaic-associated virus (AsMaV)infected Populus tremula [22].
Following total RNA extraction, first strand cDNA synthesis was performed using random hexamer oligonucleotides. Two µL of the primer (50 µM) was added to 1 µg of RNA and topped up with to 12.   Table 1. Segment-specific primer sets for the detection of the novel virus identified in diseased Fraxinus spp.

Genome Segment
Primer Fifty-two F. excelsior and 57 F. ornus leaf samples, collected during the vegetation periods between 2016 and 2019, were investigated for the presence of the novel virus. Samples originated from urban locations in Germany, Northern Italy, Southern Sweden and Switzerland as listed in the Section 2.1. Suspected virus symptoms were documented and all samples were tested by RT-PCR, applying the virus segment-specific primer sets (Table 1). Twelve samples from both F. excelsior and F. ornus displayed no observable symptoms of virus infection.

Variability Analysis of Virus Variants
Genetic variability in the novel virus was investigated using partial PCR products generated during the study. Fifty-two RNA3 (18 from F. excelsior and 34 from F. ornus), 17 RNA4 (10 F. excelsior and 7 F. ornus) and 49 RNA5 (17 F. excelsior and 32 F. ornus) PCR amplicons were purified (MSB Spin PCRapace kit, Invitek Molecular, Berlin, Germany), bidirectionally sequenced (Macrogen Europe B.V., Amsterdam, The Netherlands) and aligned in BioEdit. Sequences were trimmed to be in-frame and translated. Partial RNA3-5 sequences were 297 nt (99 aa), 453 nt (151 aa) and 357 nt (119 aa) in length, respectively, with all being found within the coding region of the respective genome segments. Pairwise sequence comparisons were performed in MEGA-X, version 10.0.5 [37] using the distance estimation function, applying the maximum composite likelihood model to nucleotide sequences and the equal input model to amino acid sequences. The standard error of the mean genetic distance was calculated for each set of values.

Electron Microscopic Studies
Leaves from a virus-infected ash tree (F. excelsior, E57075) exhibiting chlorotic line pattern and mottle accompanied by leaf deformation such as curling and shoestring were compared with asymptomatic leaf tissue from a healthy ash tree (E57079) sampled in Crans Montana in September 2017. Specimens were kept in sterile water until sample preparation. For this purpose, fragments of symptomatic or asymptomatic petiole and leaf blade were pre-fixed using a paraformaldehyde (2%) glutaraldehyde (3%) solution at pH 7.0 (0.07 M PIPES buffer) for 3 h at room temperature and descending atmospheric pressure. To ensure the infiltration of the resin, samples were incubated for 2 h under vacuum (400 bars). Polymerization of the resin took place at 60 • C for 48 h. Infiltrated samples were then thin cut (0.08 µm) and stained with 2% uranyl acetate followed by lead citrate, according to Reynolds [38]. Thin sections were observed with a transmission electron microscope CM10 (Philips, Hillsboro, OR, USA) with a Mega View II camera.

Identification and Assembly of the Complete Genome of a Novel Virus in Shoestring-Diseased Fraxinus spp.
In the HTS datasets obtained by RNAseq from pooled leaf samples of two different ash species (F. excelsior and F. ornus) exhibiting chlorotic ringspots and lines, vein yellowing and leaf deformations such as curling and shoestring ( Figure 1), we could identify several scaffolds relating to five different genome segments of the plant infecting emaraviruses ( Table 2). The first HTS dataset (194 MB), originating from two diseased F. excelsior trees (Basel, Switzerland), yielded a total of 1,657,082 reads (1,199,175 mapped reads) leading to 599 assembled scaffolds with a maximum length of 7048 nucleotides. Six scaffolds produced significant hits in blastx to the RNA1 of Actinidia emaravirus 2 (AcV-2), Pigeonpea sterility mosaic emaravirus-2 (PPSMV-2), and FMV, respectively. Two scaffolds related to RNA2 of PPSMV-2, one scaffold showed 99% sequence identity to the coding region of RNA3 of the putative novel emaravirus pea-associated virus, while two other scaffolds exhibited sequence identities of 77% to the RNA4-encoded movement protein (MP) of FMV and 46% with the P6 of AcV-2, respectively. All identified scaffolds relating to emaraviruses were represented by 2.57% of the total reads in the HTS pool from F. excelsior. The second HTS dataset was generated from approximately three times the amount of ds cDNA of pooled symptomatic leaf material of 12 diseased F. ornus trees (Hamburg, Germany) as compared to the first RNASeq run from the F. excelsior sample pool. This allowed a deeper sequencing yielding in 896 MB (6,633,156 total reads, 5,930,967 mapped reads) assembled into 10,568 scaffolds with a maximum length of 37,836 nucleotides. As shown in Table 2, we found five scaffolds relating to emaravirus-encoded protein sequences encoded by RNA1-RNA4, and RNA6 of the aforementioned emaraviruses, respectively. Although the overall percentage of total reads (0.23%) of these five scaffolds is lower than in the first HTS dataset, the mean coverage of each scaffold was considerably higher in the HTS dataset from F. ornus. Therefore, it was possible to assemble five nearly complete genome segments of a novel virus relating to emaraviruses. Sequence comparison at the nucleotide level of overlapping regions of all obtained scaffolds identified in F. ornus and F. excelsior showed high sequence identities between 95.3-97.1% (RNA1), 95.0-96.2% (RNA2), 97.4% (RNA3), 95.0% (RNA4) and 96.8% (RNA5), respectively, thus confirming that the same virus is present in both datasets derived from European ash and manna ash, respectively (data not shown). Identified scaffolds originating from the F. excelsior and F. ornus HTS datasets, which relate to emaraviral sequences, were deposited in the European Nucleotide Archive (ENA) and are available under the study number: PRJEB46740 (Accession numbers OU466864-OU466879). Scaffolds producing significant hits in blastx with sequences from other viral genera with plant infecting viruses were not found in either dataset. Table 2. Assembled scaffolds of the novel virus relating to emaravirus RNA segments identified by HTS in Fraxinus excelsior sample pool (n = 2; E55262, E55268) and Fraxinus ornus sample pool (n = 12; E57321, E57322, E57323-E57325, E57328-E57330, E57335-E57337, E57339). Identified scaffolds are presented with size in nucleotides (nt), mean coverage and percentage of mapped reads. Best matches in blastx search including query coverage, E value and sequence identity are given for each scaffold. The sample E55270 from F. excelsior was analyzed further to obtain the complete sequence of all five genome segments of the novel virus identified in Fraxinus. The complete RNA1 was determined by assembling eight overlapping sequenced RT-PCR products, generated by eight specific primer sets (Table S1, Supplementary Information). The complete sequences of the RNA2-5 were determined by cloning and sequencing of full-length PCR products. No further RNA segments were amplified by this method, suggesting that the five RNAs represent the complete genome of the novel virus. The 234 3 terminal nucleotides of the RNA1 were confirmed by a tailing PCR. Likewise, the 13 nucleotides of the 3 terminal of RNA2 could be confirmed, which were determined by the PDAP213 primer sequence. At least two individual clones were analyzed for sequence confirmation.

Genome Organization of the Novel Virus Identified in Fraxinus spp. and Comparison with Members of the Genus Emaravirus
The novel virus identified in F. excelsior (E55270) contained five RNA genome segments, each encoding a single open reading frame (ORF). Protein sequences encoded by ORFs of RNA1-RNA4 showed sequence identities and contained conserved domains characteristic for the viral replicase (RdRp, RNA1), a glycoprotein precursor (GPP, RNA2), the nucleocapsid protein (N, RNA3), and the movement protein (MP, RNA4) to certain members of the genus Emaravirus (Table 3 and Table S2, Supplementary Material). The RNA1 was found to be 7031 nucleotides (nt) long, with a single ORF (6891 nucleotides, 2297 aa) between 45 and 6935 nt in the positive sense (mRNA) encoding the viral RdRp. The 268 kDa protein exhibited a highest sequence identity of 77.9% to the RdRp of AcV-2. The RNA2 is 2243 nt in length, with a single ORF (1926 nt, 642 nt) between 57 and 1982 nt, encoding a 73 kDa glycoprotein precursor (GPP). The GPP was found to have the highest similarity to the pigeon pea sterility mosaic virus -2 (PPSMV-2) GPP (58.2%). The RNA3 is 1483 nt long, with a single ORF (948 nt, 316 aa) between 104 nt and 1051 nt, encoding a 35 kDa nucleocapsid protein (N). The N protein is 99% identical to the complementary P3 protein of the recently identified pea-associated emaravirus (PaEV) [39]. The RNA4 is 1518 nt in length with a single ORF (1086 nt, 362 aa) between 83 and 1168 nt, thought to encode the movement protein (MP). The 41 kDa protein was found to be 76.1% identical to the MP of FMV. The RNA5 is 1333 nt in length, with a single ORF (696 nt, 232 aa) between 70 and 765 nt. The function of the 26 kDa protein remains unknown. However, blastx revealed similarities to other undetermined proteins, of classified emaraviruses, including AcV-2 P6 (42.4%) and AsMaV (34.4%).
The 5 and 3 termini of all five identified genome segments are conserved and show high complementary to each other. Therefore, genome segments may be capable of circularization by forming a panhandle structure through hybridization of the terminal sequences. However, the length of the terminal complementarity region varies in the different genome segments. The first 13 nucleotides seem to be strictly conserved in all RNAs, while complementarity of RNA1 extends to the 20 5 and 3 terminal nucleotides. The terminal complementarity regions of RNA2, RNA4 and RNA5 are 19 nt long, while in RNA3 only a stretch of 17 terminal nts are fully complementary to each other ( Figure S1, Supplementary  Information).
Due to the current species demarcation criteria established by the ICTV [23] of more than 25% sequence divergence in RNA1-RNA3 encoded proteins of emaraviruses (RdRP, GPP, N), we consider the novel virus identified in diseased Fraxinus to be a distinct species within the genus. Due to the prominent leaf-narrowing symptom frequently observed, the name ash shoestring-associated virus (ASaV) is proposed. The complete genome organization of ASaV is shown in Figure 2 and the complete RNA1-5 sequences were deposited in the ENA database under the accession numbers OU466880-OU466884, respectively.
RdRp of AcV-2. The RNA2 is 2243 nt in length, with a single ORF (1926 nt, 642 nt) be 57 and 1982 nt, encoding a 73 kDa glycoprotein precursor (GPP). The GPP was fo have the highest similarity to the pigeon pea sterility mosaic virus -2 (PPSMV-2 (58.2%). The RNA3 is 1483 nt long, with a single ORF (948 nt, 316 aa) between 104 1051 nt, encoding a 35 kDa nucleocapsid protein (N). The N protein is 99% identica complementary P3 protein of the recently identified pea-associated emaravirus [39]. The RNA4 is 1518 nt in length with a single ORF (1086 nt, 362 aa) between 1168 nt, thought to encode the movement protein (MP). The 41 kDa protein was fo be 76.1% identical to the MP of FMV. The RNA5 is 1333 nt in length, with a sing (696 nt, 232 aa) between 70 and 765 nt. The function of the 26 kDa protein re unknown. However, blastx revealed similarities to other undetermined prote classified emaraviruses, including AcV-2 P6 (42.4%) and AsMaV (34.4%).
The 5′and 3′termini of all five identified genome segments are conserved and high complementary to each other. Therefore, genome segments may be capa circularization by forming a panhandle structure through hybridization of the te sequences. However, the length of the terminal complementarity region varies different genome segments. The first 13 nucleotides seem to be strictly conserved RNAs, while complementarity of RNA1 extends to the 20 5′and 3′terminal nucleotid terminal complementarity regions of RNA2, RNA4 and RNA5 are 19 nt long, while in only a stretch of 17 terminal nts are fully complementary to each other (Figu Supplementary Information).
Due to the current species demarcation criteria established by the ICTV [23] o than 25% sequence divergence in RNA1-RNA3 encoded proteins of emaraviruses GPP, N), we consider the novel virus identified in diseased Fraxinus to be a distinct within the genus. Due to the prominent leaf-narrowing symptom frequently observ name ash shoestring-associated virus (ASaV) is proposed. The complete g organization of ASaV is shown in Figure 2 and the complete RNA 1-5 sequence deposited in the ENA database under the accession numbers OU466880-OU4 respectively.  Phylogenetic trees based on amino acid (aa) alignments of the RNA1 encoded RdRp and the RNA3 encoded N protein were calculated for the ASaV variants identified in F. excelsior (E55270) and F. ornus (HTS data) and selected putative and established emaraviruses [27]. The phylogenetic relationships between the selected viruses, shown in Figure 3, exhibit four distinct subgroups A, B, C and D. Viruses in subgroups A and B are more closely related, sharing a common origin. The subgroups C and D are more distantly related. The RNA1 and RNA3 predicted proteins of ASaV were found to consistently cluster in subgroup A, together with those of classified emaraviruses including AcV-2, FMV and PPSMV-2. The trees also confirm the high sequences identities between the ASaV variants identified in F. excelsior and F. ornus, with both occupying a single terminal clade. The close phylogenetic relationships depicted therefore support the proposal of ASaV being a new, distinct member of the genus Emaravirus.
Phylogenetic trees based on amino acid (aa) alignments of the RNA 1 encoded RdRp and the RNA 3 encoded N protein were calculated for the ASaV variants identified in F. excelsior (E55270) and F. ornus (HTS data) and selected putative and established emaraviruses [27]. The phylogenetic relationships between the selected viruses, shown in Figure 3, exhibit four distinct subgroups A, B, C and D. Viruses in subgroups A and B are more closely related, sharing a common origin. The subgroups C and D are more distantly related. The RNA 1 and RNA 3 predicted proteins of ASaV were found to consistently cluster in subgroup A, together with those of classified emaraviruses including AcV-2, FMV and PPSMV-2. The trees also confirm the high sequences identities between the ASaV variants identified in F. excelsior and F. ornus, with both occupying a single terminal clade. The close phylogenetic relationships depicted therefore support the proposal of ASaV being a new, distinct member of the genus Emaravirus.

Establishment of a Virus Specific RT-PCR Detection
We tested six segment-specific primer sets (Table 1) targeting the RNA1-5 to detect the putative novel emaravirus ASaV by RT-PCR. The suitability and specificity of the detection method was evaluated using total RNA extracted from 11 symptomatic and four asymptomatic F. excelsior leaves and four additional leaf samples infected by the putative and established emaraviruses, i.e., CORaV-infected oak, MaMaV-infected maple, AsMaVinfected aspen, and FMV-infected Figure

Establishment of a Virus Specific RT-PCR Detection
We tested six segment-specific primer sets (Table 1) targeting the RNA1-5 to detect the putative novel emaravirus ASaV by RT-PCR. The suitability and specificity of the detection method was evaluated using total RNA extracted from 11 symptomatic and four asymptomatic F. excelsior leaves and four additional leaf samples infected by the putative and established emaraviruses, i.e., CORaV-infected oak, MaMaV-infected maple, AsMaVinfected aspen, and FMV-infected Figure

Geographical Distribution of ASaV and Association with Disease Symptoms
Following its discovery in 2016, further infections of ASaV in F. excelsior were confirmed in Switzerland (Basel) and at new locations discovered in Crans Montana (Switzerland), Southern Sweden (Dalby Söderskog, Sigtuna, Bromma, Vreta, Krusenstern Herrgard, Hällekis Blomberg and Västerplana), and Germany (Berlin, Hamburg). Infections in the second host species, F. ornus, were also found in Germany (Hamburg), including trees at a local tree nursery and Northern Italy (Bolzano). Sample location, host species and ASaV detection are summarized in Figure 5 and Table 4.
Fraxinus samples collected during the vegetation periods from 2016 to 2019 displayed various suspected virus symptoms, including chlorotic ringspots, line pattern, checking and a range of leaf deformations, including curling and a narrowing of the leaflets termed shoestring. Figure 1 shows a selection of documented symptoms. In total, 52 F. excelsior and 57 F. ornus samples were tested by RT-PCR applying ASaV specific primer sets. The RNA 3-5 primer sets were found to be most reliable. ASaV infections were confirmed by Sanger sequencing of PCR products (data not shown). Of 41 symptomatic F. excelsior samples, 37 (90%) were found to be infected with ASaV. ASaV was also detected in 44 of 45 (98%) symptomatic F. ornus samples (Table 4).
ASaV was also detectable in 6 of the 12 (50%) asymptomatic F. excelsior samples and 3 of the 12 (25%) samples from F. ornus with no observable disease symptoms. Our results show a strong correlation between suspected virus symptoms and ASaV infection in Fraxinus spp. ASaV was detected in 94% of the symptomatic samples overall, including all 21 samples displaying the shoestring symptom. The current known geographical range of ASaV extends across Northern and central Europe, from Southern Sweden in the north to Northern Italy in the south ( Figure 5).

Geographical Distribution of ASaV and Association with Disease Symptoms
Following its discovery in 2016, further infections of ASaV in F. excelsior were confirmed in Switzerland (Basel) and at new locations discovered in Crans Montana (Switzerland), Southern Sweden (Dalby Söderskog, Sigtuna, Bromma, Vreta, Krusenstern Herrgard, Hällekis Blomberg and Västerplana), and Germany (Berlin, Hamburg). Infections in the second host species, F. ornus, were also found in Germany (Hamburg), including trees at a local tree nursery and Northern Italy (Bolzano). Sample location, host species and ASaV detection are summarized in Figure 5 and Table 4.   Fraxinus samples collected during the vegetation periods from 2016 to 2019 displayed various suspected virus symptoms, including chlorotic ringspots, line pattern, checking and a range of leaf deformations, including curling and a narrowing of the leaflets termed shoestring. Figure 1 shows a selection of documented symptoms. In total, 52 F. excelsior and 57 F. ornus samples were tested by RT-PCR applying ASaV specific primer sets. The RNA3-5 primer sets were found to be most reliable. ASaV infections were confirmed by Sanger sequencing of PCR products (data not shown). Of 41 symptomatic F. excelsior samples, 37 (90%) were found to be infected with ASaV. ASaV was also detected in 44 of 45 (98%) symptomatic F. ornus samples (Table 4). ASaV was also detectable in 6 of the 12 (50%) asymptomatic F. excelsior samples and 3 of the 12 (25%) samples from F. ornus with no observable disease symptoms. Our results show a strong correlation between suspected virus symptoms and ASaV infection in Fraxinus spp. ASaV was detected in 94% of the symptomatic samples overall, including all 21 samples displaying the shoestring symptom. The current known geographical range of ASaV extends across Northern and central Europe, from Southern Sweden in the north to Northern Italy in the south ( Figure 5).

Genetic Variability of ASaV Based on Partial Sequences
A total of 118 partial ASaV PCR amplicons (RNA3 = 52, RNA4 = 17, RNA5 = 49) from infected F. excelsior and F. ornus were analyzed by pairwise sequence comparisons to assess genetic variability. All partial sequences lie within the coding region of the respective genome segment. Genetic distances calculated in MEGA-X are listed in Table 5. Table 5. Comparison of genetic distances of ASaV RNA3, RNA4 and RNA5 partial nucleotide and amino acid sequences. Pairwise comparisons conducted on all partial sequences for each genome segment (F. excelsior, F. ornus) and for sequences originating from each host species separately (F. excelsior or F. ornus). 0.039 (standard error 0.001-0.003) and 0-0.04 (standard error 0-0.002) at the nt and aa level, respectively. The results of this first variability analysis indicated there is no correlation between genetic variability and geographical location or ASaV host species.

Visualization of Double-Membrane Bodies in Virus-Infected Tissues
In ultrathin sections of symptomatic ash leaf tissue, we found double-membrane bodies (DMBs), with an approximate diameter of 70-80 nm, which were associated with the smooth endoplasmatic reticulum located in the cytoplasm of mesophyll cells ( Figure 6, white arrows). These DMBs are typically found in plant cells infected by emaraviruses [24] and we could observe them only in ASaV-infected leaf material but never in non-symptomatic ash leaves. This further supports the close association of the novel putative emaravirus ASaV with the chlorotic ringspots, curling and shoestring symptoms of leaves of diseased ash trees.
Number of sequences analyzed; Fragment length excluding primers; Genetic distance calculated by maximum composite likelihood model implemented in Mega-X; 4 Genetic distance calculated by equal input model implemented in Mega-X; 5 Standard error of the mean.
Results demonstrate a high level of conservation of partial sequences at both the nucleotide (nt) and amino acid (aa) level, with mean genetic distances ranging from 0.025-0.039 (standard error 0.001-0.003) and 0-0.04 (standard error 0-0.002) at the nt and aa level, respectively. The results of this first variability analysis indicated there is no correlation between genetic variability and geographical location or ASaV host species.

Visualization of Double-Membrane Bodies in Virus-Infected Tissues
In ultrathin sections of symptomatic ash leaf tissue, we found double-membrane bodies (DMBs), with an approximate diameter of 70-80 nm, which were associated with the smooth endoplasmatic reticulum located in the cytoplasm of mesophyll cells ( Figure  6, white arrows). These DMBs are typically found in plant cells infected by emaraviruses [24] and we could observe them only in ASaV-infected leaf material but never in nonsymptomatic ash leaves. This further supports the close association of the novel putative emaravirus ASaV with the chlorotic ringspots, curling and shoestring symptoms of leaves of diseased ash trees.

Discussion
We identified a novel plant virus infecting European ash (Fraxinus excelsior) and manna ash (F. ornus) across Europe. Infected trees displayed symptoms including chlorotic ringspots, mottle, line pattern and leaf deformations of curling and shoestring. The virus is proposed as being a new member of the plant virus genus Emaravirus (Order: Bunyavirales; Family: Fimoviridae) [24], fulfilling established criteria for classification to the genus

Discussion
We identified a novel plant virus infecting European ash (Fraxinus excelsior) and manna ash (F. ornus) across Europe. We identified five RNA genome segments (RNA1-5) of the novel virus in both affected ash species by HTS in combination with full-length PCR and Sanger sequencing and could attribute them to the putative emaravirus ASaV. The sequence divergence of the five ASaV encoded proteins to the analogous proteins of other emaraviruses known to date clearly demonstrates the novel virus ASaV as being a separate species. Despite the fact that the divergence of ASaV encoded proteins is less than 25%, when compared to the divergence of respective proteins of different members of the genus, this applies only to one of the respective protein sequences (RdRp, GPP, N), which are considered as relevant for species demarcation on a molecular level [23]. When comparing the ASaV encoded protein sequences to AcV2, only the RdRp shows identity above 75%, while the other four encoded proteins (GPP, N, MP, P26) exhibit diversities above 25% to the respective homologs. In the case of the putative N protein of ASaV, the aa sequence is nearly identical with the P3 of PaEV, suggesting that the emaravirus identified in pea may be the same as the novel virus we identified in ash. However, as RdRp, GPP, and MP sequences of PaEV are not available in the databases to date, this cannot be ascertained. In any case, the putative ASaV GPP and P26 demonstrate a diversity greater than 25%, when compared to the analogous protein of any other emaravirus. Therefore, the species demarcation criteria established for the genus Emaravirus [23] are fulfilled. Whether these five RNA segments represent the complete genome of ASaV remains to be clarified. Although all emaraviruses comprise of at least four core RNA segments, the number of additional genome segments is known to vary [23,33] with the 10 of perilla mosaic-associated virus (PerMV) being the highest number known to date [40]. However, as the five RNA segments attributed to AsaV were identified in two HTS datasets from pooled F. excelsior and F. ornus and also additionally by full-length RT-PCR, our results strongly suggest that we determined the complete genome of AsaV.
The fifth RNA5 encoded protein of AsaV was found to be homologous to a group of similar-sized, non-structural emaravirus proteins, designated P26-P28 according to their size (26-28 kDa), as described by von Bargen et al. [22]. The function and importance of these analogous proteins in the lifecycle of these emaraviruses, however, remains elusive. A specific interaction with host or vector proteins, or a role in virus pathogenicity has previously been speculated for such unknown proteins [33,41]. Viral suppressors of host RNA silencing mechanisms, for example, are well documented in other plant viruses [42,43], including the emaravirus High Plains wheat mosaic virus (HPWMoV), in which the P8 was reported to have suppressor activity [44,45]. However, no significant homology of the ASaV P26 was found to this protein.
Phylogenetic analyses based on amino acid alignments of the RdRp and N proteins of ASaV and selected emaraviruses revealed cladograms with four distinct subgroups (A, B, C and D), whereby ASaV consistently clustered in subgroup A. This not only corroborates the cumulative findings of previous emaravirus studies [21,22], but also confirms ASaV to be a new distinct emaravirus species, being to the best of our knowledge the first to have been detected in Fraxinus spp. Surprisingly, the N protein (RNA3) of ASaV was found to be nearly identical (99%) to that of the recently discovered pea-associated emaravirus (PaEV) which infects green peas (Pisum sativum L.) [39]. This raises the question as to whether both viruses are in fact the same species. With the exception of mechanical transmission of emaraviruses occurring under experimental conditions [46][47][48], emaraviruses have largely been found to infect a single host or relatively few, closely related host species. Although a natural infection of both ash and peas would therefore seem unlikely, it would raise new epidemiological questions, if this were the case. However, as only the RNA3 of PaEV is currently available in the databases for analysis, the characterization of further PaEV genome segments is necessary to determine the phylogenetic relationship of both viruses.
3 terminal ends of the RNA1 and RNA2 were confirmed by tailing PCR as being the same as those generated by the genus specific PDAP213 primer [33]. The remaining terminal ends of ASaV are still to be confirmed. Possible explanations for the difficulties experienced in amplifying particular termini by this method may include an unequal distribution of virus particles in the host plant, or low concentrations of RNA complementary to the genomic viral strands. Tatineni et al. [49] found that the concentrations of negative sense RNA were up to 20 times higher than those of sense viral RNA in HPW-MoV infected samples. Rapid amplification of cDNA ends (RACE) [50,51] or inverse PCR methods [52] could be used as alternative methods in future studies to confirm remaining terminal sequences.
During this study, a successful RT-PCR-based method was established for the detection of ASaV in Fraxinus spp., applying virus-specific primer sets targeting the five RNA segments. All primer pairs are located within the coding regions of each genome segment and comparable detection results for ASaV were achieved for all sets. The existence of more divergent ASaV variants, which may not be detected by these primer sets, cannot, however, be excluded due to the limited sequence information currently available. Tatineni et al. (2014) [49] demonstrated the existence of HPWMoV variants differing significantly in the identity of only a single genome segment. To avoid the possibility of false negative results, samples should therefore be tested by applying primer sets for at least two different genome segments, as suggested by McGavin et al. [47]. However, as a 5% maximum nucleotide diversity was found between two ASaV variants, each identified in a different Fraxinus species, this suggests that the overall sequence variability of ASaV is not expected to be high.
A close association between the presence of AsaV and the described host symptoms could be demonstrated in F. excelsior and F. ornus, indicating the virus to be the probable causal agent. Symptoms of chlorotic ringspots, mottle and line pattern are typical of many other emaravirus infections, including rowan (Sorbus aucuparia L.) infected with European mountain ash ringspot-associated virus (EMARaV) [3,26]. Many emaravirusinfected plants are also known to suffer a general decline over time, often leading to their eventual death [3,41,53]. Although several ash trees displaying premature leaf loss and branch dieback were observed during field excursions, these symptoms were most likely attributed to the fungus Hymenoscyphus fraxineus responsible for the ash dieback disease. This remains to be determined for ASaV-infected ash.
Interestingly, a range of leaf deformations, from curling to a severe narrowing of the leaflets known as shoestring, were also found to be associated with the presence of the virus in ash. Although symptom severity is likely a function of the developmental stage of the host at the time of infection and length of infection, Kallinen et al. [54] reported that symptom severity in rowan was not correlated with the EMARaV titer, as might be expected if this were the case. Symptom expression in the infected rowan was also found to vary over consecutive seasons. Virus strain and virulence, tolerance of the host plant and even differences in environmental conditions are factors likely to affect symptom development. Despite the absence of further sequences of known viruses in HTS datasets, Koch postulates establishing a direct causal link between ASaV and the disease symptoms have yet to be fulfilled. The difficulty of isolating the necessary infectious emaravirus particles to achieve this is well known [24]. The use of infectious full-length ASaV clones or graft transmission studies, as outlined by von Bargen et al. [22] and Führling and Büttner [55], would provide conclusive evidence of ASaV as the cause of the disease in previously healthy host plants. ASaV was also detected in a number of trees displaying no observable symptoms, thus confirming a latent infection of F. excelsior and F. ornus. Latent emaravirus infections have been documented for a number of other emaraviruses, including rowan infected with EMARaV [54] and pear infected with pear chlorotic leaf spot associated virus (PCLSaV) [56]. It is currently unclear whether the infected ash lacking disease symptoms were simply sampled at an early stage of infection and may have exhibited ASaV-related symptoms at other timepoints over the vegetation period or whether they remained latent, with disease development possibly only being induced during the influence of additional abiotic or biotic stress factors, for instance fungal infections.
ASaV infections were found to be widespread across Europe, ranging from Southern Sweden in the north to Bolzano, Northern Italy in the south. The native habitat of both host Fraxinus species, however, extends much further than the investigated area covered in this study, which is particularly true for F. ornus, which is native to the Mediterranean region [6]. These preliminary findings suggest however that ASaV is likely to be present throughout its host's range. As other emaraviruses have been found to infect more than one host species, the possibility of ASaV infecting additional Fraxinus spp., including nonnative ornamental species such as the red ash (F. pennsylvanica), is possible. EMARaV, for example, is known to infect several closely related species and interspecies hybrids of the family Rosaceae [26,28,[57][58][59]. A European-wide ASaV monitoring would most probably reveal a more extensive geographical distribution and additional host species of the novel virus ASaV.
Preliminary sequence data based on partial RNA3, 4 and 5 sequences, suggest a high level of genetic identity of ASaV at both the nucleotide and the amino acid level. There also appears to be no correlation between genetic variability and either host species or geographical location of trees sampled. Although more extensive sequence data are required from future studies to confirm this observation, these first results are in line with the theory that RNA virus populations are generally genetically stable despite having a high potential for genetic diversity [60]. Biological properties of viruses such as host type and mode of transmission are thought to play an important role in shaping genetic diversity, with purifying selection pressure constraining the variability of proteins required for interactions between virus, host and vector [54,60,61]. High sequence identities between emaravirus variants at different locations or between different hosts have been documented in several studies. An identity of over 97.1% for EMARaV was reported by Grimova et al. [58] in different, closely-related hosts. With respect to geographical location, Roßbach et al. [62], von   [61] and Kallinen et al. [54] all reported a high level of sequence identity of EMARaV at locations across Northern Europe and Russia. Roßbach et al. [62] for example, found 95.3-100% nucleotide identities in the RNA3 coding region (944 nt) of 36 EMARaV variants across Northern Europe. Interestingly, the authors also identified a small number of EMARaV variants with notably lower sequence identities, an observation also documented in rose rosette virus (RRV) [53]. The evolutionary significance of comparably few divergent variants in a virus population remains to be investigated. Although other studies have demonstrated genetically distinct emaravirus populations, highly similar variants were also found at geographically distant locations [63,64]. Although three distinct FMV populations were identified, based on the analysis of partial RNA3 fragments (873 nt) of 53 virus variants worldwide, highly similar virus variants at geographically distant regions were also found [64]. It was suggested that the trading of plants, leading to the long-distance migration of particular virus variants, may give rise to such observations [64]. The detection of ASaV at a tree nursery in our study indicates that the trading of ash trees could play an important role in determining genetic variability, but also in aiding the spread of the virus in general.
The increasing use of HTS technologies is revolutionizing the diagnosis and characterization of plant pathogens, particularly viruses [12,13]. Despite the new opportunities created by HTS for detecting novel viruses, an in-depth biological characterization remains essential for evaluating the significance and impact of these new species. This has led to a bottle neck in the evaluation of new species, and therefore a new standardized framework is required to solve this problem [65]. In our study, we report the discovery of a novel emaravirus in diseased Fraxinus spp. with the proposed name ash shoestring-associated virus (ASaV) and provide the first important information to genome structure, phylogenetic placement and genetic variability. The establishment of a reliable RT-PCR-based detection method achieved in this study and its refinement in subsequent investigations will not only help unravel the epidemiology of ASaV, but will also aid the development of potential management strategies in the future, such as the testing and certification of plant material, ensuring healthy, virus-free ash for planting. A focus on virus infection will become particularly important in the selection of healthy trees for urban and forest plantings, with tolerance to the changing climatic conditions due to global warming. We have shown that ASaV-infected ash trees are present in several European countries. They have been found in rural stands and forests such as the Swiss mountainous regions near Basel and Crans Montana but also have been propagated in tree nurseries and planted as street trees in urban sites in Germany and Sweden. Therefore, a European-wide ASaV monitoring in the future would most probably reveal a more extensive range as well as additional host species of the novel virus, thus allowing a better estimation of the impact of the identified novel virus and its role in the observed disease of investigated ash species. Although our observations of cauliflower gall growths on symptomatic ash strongly suggest Aceria fraxinivora to be a likely candidate for the natural transmission vector of ASaV, this also remains to be confirmed.
Viruses are most likely to have an important role in the health of plants and forests [3], and a better understanding of their interaction with other stress factors affecting plants is essential. A loss of vitality and a predisposition of trees to other diseases resulting from viral infection is not only relevant for ASaV when selecting healthy ash trees, suitable for changing climatic conditions, but also with regard to the ash dieback disease, currently causing the decline in Fraxinus across Europe.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/f12111574/s1, Figure S1: Nucleotide sequences of the terminal complementarity regions of the five RNA genome segments of the novel emaravirus (E55270) infecting ash species are shown. Strictly conserved 13 terminal nucleotides (represented by the PDAP213 primer) are shown in bold. Sequences confirmed by tailing PCR are underlined. Table S1: Primer sets for the complete amplification of the ASaV RNA1 and tailing PCR for ASaV RNA1 and RNA2. Table S2: Accession numbers of putative and classified emaraviruse sequences used for amino acid sequence comparisons (Table 3)   Data Availability Statement: Identified scaffolds originating from the F. excelsior and F. ornus HTS datasets, which relate to emaraviral sequences as well as assembled sequences of the complete genome segments of the novel emaravirus, were deposited in the European Nucleotide Archive (ENA) and are available under the study number: PRJEB46740 (Accession numbers OU466864-OU466879).