A Distinct Arabidopsis Latent Virus 1 Isolate Was Found in Wild Brassica hirta Plants and Bees, Suggesting the Potential Involvement of Pollinators in Virus Spread

During our search for aphid-pathogenic viruses, a comovirus was isolated from wild asymptomatic Brassica hirta (white mustard) plants harboring a dense population of Brevicoryne brassicae aphids. The transmission-electron-microscopy visualization of purified virions revealed icosahedral particles. The virus was mechanically transmitted to plants belonging to Brassicaceae, Solanaceae, Amaranthaceae, and Fabaceae families, showing unique ringspot symptoms only on B. rapa var. perviridis plants. The complete viral genome, comprised of two RNA segments, was sequenced. RNA1 and RNA2 contained 5921 and 3457 nucleotides, respectively, excluding the 3′ terminal poly-adenylated tails. RNA1 and RNA2 each had one open-reading frame encoding a polyprotein of 1850 and 1050 amino acids, respectively. The deduced amino acids at the Pro-Pol region, delineated between a conserved CG motif of 3C-like proteinase and a GDD motif of RNA-dependent RNA polymerase, shared a 96.5% and 90% identity with the newly identified Apis mellifera-associated comovirus and Arabidopsis latent virus 1 (ArLV1), respectively. Because ArLV1 was identified early in 2018, the B. hirta comovirus was designated as ArLV1-IL-Bh. A high-throughput-sequencing-analyses of the extracted RNA from managed honeybees and three abundant wild bee genera, mining bees, long-horned bees, and masked bees, sampled while co-foraging in a Mediterranean ecosystem, allowed the assembly of ArLV1-IL-Bh, suggesting pollinators’ involvement in comovirus spread in weeds.


Introduction
The Comovirus genus belongs to the Comovirinae subfamily in the Secoviridae family, which contains a wide variety of members [1].Comoviruses have a narrow host range that belongs to the Leguminosae family, showing characteristic mosaic and mottling symptoms.Viruses belonging to this genus are non-enveloped with icosahedral morphology containing a bipartite positive sense RNA genome [2].The polyprotein of RNA1 and the small polyprotein of RNA2 are processed into five and three cleavage products, respectively [3,4].
A 3C-like proteinase encoded by RNA1 has unique cleavage sites for polyprotein processing.RNA1 is necessary for viral replication and includes a protease cofactor, a helicase, a viral protein-genome-linked (VPg), a 3C-like proteinase, and an RNA-dependent RNA polymerase (RdRp).The N-terminal protease cofactor of RNA1 assists in processing the RNA2 polyprotein [5].RNA2 encodes the movement protein and large and small capsid proteins required for viral movement [6].Comoviruses are transmitted by beetles [7] and mechanical means [2].A recent study on Arabidopsis latent virus 1 (ArLV1), a comovirus infecting Arabidopsis thaliana plants of the Brassicaceae family, has shown a highly efficient seed transmission of the virus [8].
Mechanically transmitted viruses are sometimes transmitted by non-host vectors, primarily via mechanical adherence.Bee colonies were considered vectors of plant viruses [9][10][11][12][13][14][15], and metagenomics studies of bee populations were found to be beneficial in the identification of plant viruses in various ecosystems [16][17][18][19][20][21].The ability to share viruses between managed and wild bees [22,23] may lead to a wide plant virus spread within a specific geographic region [18].For our study of plant virus spread in wild vegetation, we have sampled the managed honeybees (Apis mellifera) and three abundant wild bee genera, mining bees (Andrena), long-horned bees (Eucera), and masked bees (Hylaeus), while co-foraging in a Mediterranean shrub land in central Israel.These four bee genera differ in their nesting, diet, pollen-collecting organs, sociality, body size, and seasonal activity [24].The honeybee A. mellifera is a generalist forager, and compared to most of the wild bees in our system, they are long-distance foragers and are active during most of the year in the examined region.They live in large colonies with a high level of nest-mate interactions.Unlike A. mellifera, some mining bees (Andrena), long-horned bees (Eucera), and masked bees (Hylaeus) have narrower diets [25][26][27]; they are short-distance foragers [28,29] and active for short periods limited to a few weeks per year.In the current study of plant virus spread, we have recently found that all the sampled genera harbored a comovirus, which showed similarity with a previously identified comovirus found in asymptomatic Brassica hirta (white mustard) plants in Mediterranean shrubland ecosystems in central Israel.The virus found in mustard showed a high similarity to the recently described latent ArLV1.This study emphasizes the likely involvement of wild and managed bees in the spread and persistence of plant viruses in wild vegetation, allowing a possible spread to adjacent economically important or research-essential crops.

Wild Plant Collections
During the spring of 2011, asymptomatic B. hirta plants, heavily infested with the cabbage aphid (Brevicoryne brassicae), were collected from a no-till farming area for virus identification and characterization.The collected B. hirta plants, found to harbor a comovirus, were used as source material to establish a mustard comovirus culture.During the winter of 2022, asymptomatic B. hirta plants were collected from fields near commercial hives, and the presence of the mustard comovirus was confirmed by RT-PCR using the primer sets F10-R8 and F15-R13 for RNA1 and RNA2, respectively (Table 1) (see below).

Viral Particle Purification and Transmission Electron Microscopy Visualization
Virus purification was conducted as described previously [30] with some modifications.Briefly, 100 g of mustard leaves were homogenized in 200 mL of 0.067 M sodium phosphate buffer (pH 7.2) and 50 mL of 0.1 M ascorbic acid.The homogenate was centrifuged at 9700× g for 10 min at 4 • C. The separated supernatant was mixed with 10% chloroform and shaken for 10 min at 4 • C, followed by a pH adjustment to 5.3.The mixture was centrifuged at 200,000× g for 2 h in a fixed-angle (Beckman Ti 35) rotor.The pellets were suspended in a cold 0.01 M sodium phosphate buffer (pH 7.0).Transmission electron microscopy (TEM) was used to identify the purified viral particles.For TEM visualization, a purified virion sample (3.5 µL) was applied onto 300 mesh carbon-coated copper TEM grids for 30 s. Excess fluids were blotted, and after a wash with distilled water, the grids were stained with 2% uranyl acetate and visualized using a Tecnai G 2 , FEI-Philips (Philips, Eindhoven, The Netherlands).

Viral RNA Extractions from Purified Virions
Viral RNA extractions from viral particles were conducted as described previously [31,32].In brief, purified virions were subjected to RNase-free DNase I digestion (Promega; Madison, WI, USA) (1 h at 37 • C) followed by Proteinase K treatment (1 h at 37 • C).The viral RNA was purified using acid-phenol chloroform (Ambion/Applied Biosystems; Austin, TX, USA).Viral RNA precipitation was conducted overnight at −20 • C in the presence of glycogen (Fermentas-Thermo Fisher Scientific; Burlington, ON, Canada), 0.1 M sodium acetate, and isopropanol.The obtained viral RNA was then washed twice with 75% ethanol, air-dried for 10 min, and suspended in 40 µL of 0.01 M Tris-borate EDTA (TBE) buffer.

Double-Stranded (ds) cDNA Synthesis
Viral RNA samples from purified virions served as a template for cDNA synthesis using random hexamer and oligo dT 17 VN primers with a Maxima reverse transcriptase (Fermentas).The cDNA served as a template for the second strand synthesis using the Universal RiboClone cDNA Synthesis System (Promega, Madison, WI, USA), according to the manufacturer's instructions.The ds-cDNA was purified using the Zymoclean kit (Zymo Research, Irvine, CA, USA), and the obtained double-stranded fragments were cloned as a library into a commercial pUC19/SmaI (Fermentas).The library was then transformed into DH5α competent cells, and insert-positive colonies were cultured in LB media, including antibiotics.The plasmids were extracted using a plasmid extraction kit (Bioneer, Daejeon, Republic of Korea) and sequenced by Sanger sequencing (HyLabs, Rehovot, Israel).
Based on the obtained sequences, the primer pairs were designed (Table 1) and used for RT-PCR to obtain the unidentified genome segments.The 5 ′ end of the genome was identified using the RACE strategy on cDNA derived from viral RNA extractions, while the 3' end was sequenced using the oligo dT 17 VN primer.

Viral RNA Extractions from Plants
The Viral RNA extraction kit (Bioneer, Daejeon, Republic of Korea) was used to extract viral RNA from plant leaves.A sample of 0.5 g of plant tissue was ground in the presence of the supplied RNA extraction buffer, and RNA was extracted according to the manufacturer's instructions.

RT-PCR Amplifications
Primers were designed based on sequences obtained from the ds-cDNA sequencing to amplify the viral genome.The cDNA served as a template for PCR using JMR polymerase (JMR, Kent, UK) or Advantage 2 Polymerase Mix (Clontech-Takara Bio, Mountain View, CA, USA) and specific primer sets (Table 1).Amplicons were separated on a 1% Agarose gel (HyLabs, Rehovot, Israel), extracted with a Zymoclean Gel DNA recovery kit (Zymo Research, Irvine, CA, USA), and Sanger sequenced (HyLabs).The primer sets F10, R8, and F15, R13 (Table 1) were used for the detection of the mustard comovirus RNA1 and RNA2, respectively.For the detection of Actin, the primer set was as follows: F 5 ′ ATGCCAACACTGTCCTTTCTGG 3 ′ and R 5 ′ GACCCACCAATCCATACGGA 3 ′ .

Viral Genome Assembly
Nucleotide sequences obtained from the ds-cDNA library were analyzed using the NCBI database and Basic Local Alignment Search Tool (BLAST, NIH, Bethesda, MD, USA).Genome sequence assembly was generated by DNAMAN (Lynnon BioSoft; Montreal, QC, Canada) and SnapGene (GSL Biotech LLC, Boston, MA), and the reference genomes of the turnip ringspot virus (TuRSV) (accession numbers FJ712026 and FJ712027) [33], radish mosaic virus (RaMV) (accession numbers AB295643 and AB295644), and ArLV1 (accession numbers MH899120 and MH899121).The deduced amino acid sequence was obtained using the ORFfinder from NCBI and SnapGene software version 6.0.4.

Characterization of the Viral Coat Protein
The molecular weight (MW) of the viral coat protein (CP) subunits was estimated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE).The purified virion preparation was fractionated using a 12% polyacrylamide gel containing 0.4% SDS, as described by Laemmli [34].The resulting protein bands were visualized by staining the gel with Coomassie Brilliant Blue (Sigma-Aldrich, Burlington, MA, USA).

Bee Sampling
Bee sampling was conducted between February and May 2021 at five shrubland sites in the Judean foothills, a Mediterranean agroecosystem in central Israel (Table S1).Each site was sampled 5-7 times, at 7-10 day intervals.Sampling included two sets of observations on bees' flower visits, 30 min each, during peak bee activity hours.Subsequently, 15 individuals from the most abundant bee genera were netted.These included managed honeybees (A.mellifera) in all sampling days (n = 540), as well as solitary mining bees (Andrena spp.; n = 376), long-horned bees (Eucera spp.; n = 208), and masked bees (Hylaeus spp.; n = 109) (Table S1).The captured bees were immobilized on ice, identified at the genus level, and subsequently kept on dry ice and stored in a −80 • C freezer until processed.

Local Flora Characterization
At the beginning of every sampling day, the local flora available for forage was recorded by counting flowering units for every species from ten 1-square-meter quadrates that were located on two transects in the sites.

Viral RNA Extractions from Bees
Total RNA was extracted from every bee individually using the standard phenol and guanidinium isothiocyanate (TriReagent, Cincinnata, OH, USA) protocol.Briefly, each bee was homogenized in 400 µL of BioTri (BioLab Ltd., Jerusalem, Israel), followed by 10 min incubation.A total of 80 µL of chloroform was added to each tube, followed by vortex, 5 min incubation, and 15 min centrifugation.The upper phase was mixed in a ratio of 1:1 (v:v) with 9 M LiCl in isopropanol, followed by 1 h incubation at −20 • C for precipitation.Following an 8 min centrifugation, the pellet was washed with 75% ethanol, and elution was performed with 40 µL DDW.The RNA from all the bees was pooled according to genus and enriched for viral RNA using the Viral RNA Extraction kit (Bioneer, Daejeon, Republic of Korea) according to the manufacturer's instructions.The presence of the mustard comovirus was tested by RT-PCR followed by Sanger sequencing as detailed above.

High Throughput Sequencing (HTS) Analyses
In order to diagnose viruses that are both poly-adenylated and non-poly-adenylated, the RNA pools from the four bee genera were in vitro poly-adenylated and sequenced using Illumina Hiseq2500 (50 cycles) platform (Technion Genome Center, Haifa, Israel).Low-quality sequences were filtered and trimmed using Trimmomatic version 0.39 [35].Clean reads were scanned for matched viral sequences using VirusDetect software version 1.7 [36].For mapping the reads, VirusDetect software (version 1.7) employed a pipeline that combined de novo assembly using Velvet software version 1.1.07[37] with mapping to plant viral references from Genbank utilizing the Burrows-Wheeler Alignment tool (BWA) [38].In parallel, we also used trinity assembler version v2.13.2 [39]; the assembled contigs were then searched for plant virus sequences using diamond Blastx [40] against the NCBI non-redundant protein database.The depth coverage of the viral contigs was calculated using Bowtie2 alignment [41] and Samtools version 1.7 [42].

Statistical Analysis
To determine the tendency of bees from the dominant genera to forage on brassica plants, we calculated, for every sampling day, the difference between the proportional abundance of the brassica plants in the floral community and the proportion of visits to brassica plants by each bee genus.We included in this analysis only days in which bees from the examined genera visited brassica flowers.Thus, the sample size was in accordance with the number of those days.Due to the non-normal distribution of proportional data, we used the Wilcoxon signed-rank test to compare these values with zero.A significant difference for a genus indicates its tendency to visit or avoid brassica plants relative to their proportional abundance.The masked bees were omitted from this analysis because their activity on brassica plants was too low to allow a meaningful inference.All statistical analyses were performed in R version 4.3.1.

A Phylogenetic Tree Analysis
A phylogenetic tree analysis was performed based on the Comovirinae subfamily Pro-Pol deduced amino acid region delineated between the conserved CG motif of the 3C-like proteinase and the GDD motif of the RdRp.Sequences were aligned by multiple sequence alignment using Muscle.The grapevine fabavirus (GFabV) (accession number KX241482) served as an out-group.The tree was constructed using MEGA software version 6 based on the maximum likelihood method with the parameter of 1000 bootstraps.

A Comovirus Was Identified in Asymptomatic Wild B. hirta Plants
Wild B. hirta plants harboring a dense population of B. brassicae aphids were analyzed for viral infection.The TEM visualizations of virions prepared from the plants showed icosahedral particles (Figure 1a).SDS-PAGE analysis of the virion preparation showed two potential CPs of ~20 kDa and ~40 kDa, characteristic of the Comovirus genus (Figure 1b).In parallel, several molecular methods were carried out to characterize the entire viral genome sequence.The ds-cDNA procedure was conducted on virion RNA extractions to identify the viral genus.The obtained sequences within the viral genome shared a high similarity with the Comovirus genus, which comprises two RNA molecules, RNA1 and RNA2 [43].The ds-cDNA sequences were aligned and mapped to RNA1 and RNA2 molecules using the reference viruses TuRSV (accession numbers FJ712026 and FJ712027) and RaMV (accession numbers AB295643 and AB295644) (Figure 1c,d).Based on the obtained sequences, primer pairs were designed and used for RT-PCR to attain the unidentified genome segments (Table 1).The 5 ′ end of the genome was identified using the RACE strategy on cDNA derived from the viral RNA extractions, while the 3' end was sequenced using the oligo dT 17 VN primer (Figure 1c,d).

Genome Organization
The assembled genome was analyzed using ORFfinder (NCBI) to find the open reading frames (ORFs) and to annotate putative proteins.The ORFs were then identified using the SnapGene software version 6.0.4 and aligned with ArLV1 (accession numbers

Genome Organization
The assembled genome was analyzed using ORFfinder (NCBI) to find the open reading frames (ORFs) and to annotate putative proteins.The ORFs were then identified using the SnapGene software version 6.0.4 and aligned with ArLV1 (accession numbers MH899120 and MH899121).Similar to other comoviruses [2], the mustard isolate showed one ORF within the RNA1, encoding five proteins: a Pro-co, a helicase, a viral protein genome-linked (VPg), a 3C-like proteinase, and an RdRp (i.e., Pol) (Figure 1c).The RNA2 is also comprised of one ORF encoding three known proteins: a movement protein (MP), a large capsid protein (CPL), and a small capsid protein (CPS) (Figure 1d).
The nucleotide sequence of RNA1, comprised of 5921 nucleotides excluding the polyA tail, showed an 85% similarity with the newly identified Apis mellifera-associated comovirus (AmCV) (A) and (B) (accession numbers OP972917 and OP972918), and a 76% similarity with ArLV1 (accession number MH899120).The nucleotide sequence of RNA2, comprised of 3457 nucleotides excluding the polyA tail, showed an 84% and 83% similarity with AmCV (A) and (B), respectively (accession numbers OP972919 and OP972920), but no similarity was found with ArLV1 (accession number MH899121).The deduced amino acid sequence at the Pro-Pol region, encompassing the domain between the conserved CG motif of the proteinase and the conserved GDD motif of the polymerase, shared a 96.5% and 90.3% identity with AmCV and ArLV1, respectively.The deduced amino acid sequence of the CPs shared a 93% and 82% identity with AmCV and ArLV1, respectively.As ArLV1 was already submitted to the GenBank database in 2018, and according to the International Committee on Taxonomy of Viruses (ICTV), the species demarcation criteria for the conserved Pro-Pol region and the CPs were 80% and 75% identical, respectively [2]; the Israeli B. hirta isolate was designated as ArLV1-IL-Bh and was deposited to GenBank (accession numbers OR840696 and OR840697).
In RNA1 at the Pro-Co region, the conserved amino acids 124 F, 151 W, and 180 E [44] were identified.At the helicase region, the nucleoside triphosphate binding motif had the conserved amino acids 479 GKSRVGKT 486 and the conserved DD preceded by hydrophobic amino acids 526 ILIDD 530 [45].The VPg had the conserved amino acids at positions 909 K, 914 D, 918 Y, 922 N, and 927 R [46].The 3C-like cysteine proteinase showed the catalytic cysteine at the 1100 CG motif [47,48].The RdRp showed the conserved amino acids 1423 DYSSFDGLLSK 1433 and 1484 SGFPLTVICNS 1494 , and a GDD with the preceding hydrophobic residue 1531 YGDD 1534 and 1584 FLKR 1587 [49].In RNA2, the conserved amino acids identified at the MP region were 69 P, 129 G, and 153 D (counting from the second initiation codon at nt 432: GAAATGG, which is placed in a favorable context compared to the first initiation codon at nt 120: ACTATGT) [2,50].

A Host Range Analysis
The field-collected asymptomatic B. hirta plants served as a source of virus inoculum.The virus was easily transmitted to test plants by mechanical sap inoculation.A broad host range that included members of four plant families was analyzed by RT-PCR (Table 2).Like the wild B. hirta plants, most plants that were positive for ArLV1-IL-Bh in RT-PCR analysis were asymptomatic.Unique symptoms of ringspot were observed on B. perviridis plants (Figure 2a,b), whereas infected N. benthamiana plants showed the classic comovirus symptoms of mosaic and leaf distortions (Figure 2c).

ArLV1-IL-Bh Identified in Managed and Wild Bee Populations
B. hirta (white mustard) and other Brassicaceae family species are indigenous to the Israeli wild vegetation, attracting a wide range of insect pests and pollinators such as wild bees and honeybees.Recently collected B. hirta plants infected with ArLV1-IL-Bh (Figure 2d) were found in several locations in proximity to managed A. mellifera bees and wild bee populations, which were sampled and subjected to RNA extractions and RT-PCR tests for the presence of ArLV1-IL-Bh and HTS analyses (Figure 2e, Table 3).The managed A. mellifera bees revealed, by HTS and de novo assembly followed by BLAST analyses, the presence of an A. mellifera comovirus isolate, which showed 99-100% sequence identity with the RNA1 of ArLV1-IL-Bh and a contig of 3378 nucleotides showing a 99.8% sequence identity with RNA2 of ArLV1-IL-Bh (Table 4).The Israeli A. mellifera comovirus showed a ~90% similarity with each RNA1 of ArLV1 and Zymoseptoria comovirus A and a 96% similarity with RNA1 of AmCV.The large contig covering RNA2 showed an ~80% similarity with each RNA2 of ArLV1 and Zymoseptoria comovirus A and a 92% similarity with RNA2 of AmCV.We, therefore, nominated the Israeli A. mellifera comovirus as ArLV1-IL-Am and deposited it to GenBank (accession numbers OR840694 and OR840695).Following the alignments of HTS de novo assembled contigs, the ArLV1-IL-Bh large genome segments were also identified in the wild bee populations, mining bees (Andrena), long-horned bees (Eucera) and masked bees (Hylaeus), and showed the highest similarity with the bee comovirus contigs compared to ArLV1, Zymoseptoria comovirus A, and AmCV (Table 4, Figure 3).Following the Sanger sequencing of ArLV1-IL-Bh, the obtained complete genome sequence served as a reference template for the assembly of the HTS reads in order to re- Following the Sanger sequencing of ArLV1-IL-Bh, the obtained complete genome sequence served as a reference template for the assembly of the HTS reads in order to re-assess the genome integrity of the plant comovirus in the bee samples.The analysis revealed a complete coverage of both RNA1 and RNA2 of ArLV1-IL-Bh in all analyzed bee populations (Figure 4).The highest depth in the reads, mapped to ArLV1-IL-Bh in A. mellifera compared to the wild bees (Figure 4g,h vs. Figure 4a-f), apparently reflected the high number of clean reads obtained from the honeybee sample compared to the wild bees while all bee samples showed a similar percentage of number of reads mapped to ArLV1-IL-Bh genome (Table 3).

Bee Tendency to Visit Brassica Flowers
We asked whether the presence of ArLV1-IL-Bh in all bee samples was positively correlated with the bees' tendency to visit brassica flowers.The median of the difference between the proportion of visits to brassica plants by each bee genus and the relative abundance of brassica plants in the floral community in the site was 0.2 (with an interquartile range (IQR) of 0.108-0.251)and 0.5 (IQR = 0.365-0.612)for honeybees (A.mellifera) and mining bees (Andrena), respectively (Figure 5).Both were significantly different from zero (Wilcoxon signed-rank test; p = 0.0004 and 0.007, respectively), indicating an attraction to brassica plants.Conversely, the median delta for long-horned bees (Eucera) was −0.02 (IQR = −0.106-0.077)and not significantly different from zero (Wilcoxon signed-rank test; p = 1) (Figure 5).

Bee Tendency to Visit Brassica Flowers
We asked whether the presence of ArLV1-IL-Bh in all bee samples was positively correlated with the bees' tendency to visit brassica flowers.The median of the difference between the proportion of visits to brassica plants by each bee genus and the relative abundance of brassica plants in the floral community in the site was 0.2 (with an interquartile range (IQR) of 0.108-0.251)and 0.5 (IQR = 0.365-0.612)for honeybees (A.mellifera) and mining bees (Andrena), respectively (Figure 5).Both were significantly different from zero (Wilcoxon signed-rank test; p = 0.0004 and 0.007, respectively), indicating an attraction to brassica plants.Conversely, the median delta for long-horned bees (Eucera) was −0.02 (IQR = −0.106-0.077)and not significantly different from zero (Wilcoxon signed-rank test; p = 1) (Figure 5).

A Phylogenetic Tree Analysis
A phylogenetic tree was constructed based on the alignment of amino acid sequences of the conserved Pro-Pol region.A high similarity was observed between the ArLV1-IL-Bh and the comovirus found in sampled honeybees ArLV1-IL-Am.Both ArLV1-IL isolates were clustered in a clade with AmCV, while ArLV1 and Zymoseptoria comovirus A sequences were clustered in the nearest separate clade.The Israeli strains ArLV1-IL-Bh and ArLV1-IL-Am were separated from AmCV in a sub-cluster.A cluster that included ArLV1-IL isolates, AmCV, ArLV1, and Zymoseptoria comovirus A was separated from other clusters, including TuRSV.The grapevine fabavirus (GFabV) served as an outgroup (Figure 6).

A Phylogenetic Tree Analysis
A phylogenetic tree was constructed based on the alignment of amino acid sequences of the conserved Pro-Pol region.A high similarity was observed between the ArLV1-IL-Bh and the comovirus found in sampled honeybees ArLV1-IL-Am.Both ArLV1-IL isolates were clustered in a clade with AmCV, while ArLV1 and Zymoseptoria comovirus A sequences were clustered in the nearest separate clade.The Israeli strains ArLV1-IL-Bh and ArLV1-IL-Am were separated from AmCV in a sub-cluster.A cluster that included ArLV1-IL isolates, AmCV, ArLV1, and Zymoseptoria comovirus A was separated from other clusters, including TuRSV.The grapevine fabavirus (GFabV) served as an outgroup (Figure 6).

Discussion
We have identified a strain of comovirus in field-collected wild B. hirta plants.The characterization of the mustard comovirus strain showed high similarity to AmCV (accession numbers OP972917-20) and ArLV1 (accession numbers MH899120-21), the latter submitted to GenBank early in 2018; according to the ICTV taxonomy criteria, the isolate was designated as ArLV1-IL-Bh.The alignments of the putative amino acid sequence of ArLV1-IL-Bh at the Pro-Pol region showed a 96.5% and 90.3% similarity to AmCV and ArLV1, respectively.The putative amino acids of the CPs of ArLV1-IL-Bh shared a 93% and 82% similarity with AmCV and ArLV1, respectively.The putative cleavage sites of the 3C-like proteinase were identical between the ArLV1-IL-Bh and AmCV.However, when compared to ArLV1, the putative sites identified at the Pro-Co/Helicase domain differed.In ArLV1, the putative cleavage site was 308 VAQ/AGP 313 , while in ArLV1-IL-Bh and AmCV, the cleavage site was 308 VAQ/SGP 313 .Other putative 3C-like proteinase cleavage sites in the polyprotein of RNA1 and the small polyprotein of RNA2 were similar between ArLV1 and ArLV1-IL-Bh.
We have described unique symptoms of ArLV1-IL-Bh upon infection of the mustard plant B. rapa var.perviridis.Symptom manifestations of ringspot were associated with the infection.Ringspots are known characteristic symptoms of TuRSV and viruses belonging to the Fabavirus genus in the Comovirinae subfamily.Fabaviruses are transmitted by mechanical sap inoculations and aphids in a non-persistent manner and have a broad host range [2,51].However, the phylogenetic tree analysis showed a separate clustering of ArLV1-IL isolates and TuRSV when grapevine fabavirus (GFabV) served as an outgroup (Figure 6).Host range analyses revealed that, similar to ArLV1, the Brassicaceae and Solanaceae family members were infected by the ArLV1-IL-Bh, and we added members of the Amaranthaceae and Fabaceae families to the ArLV1-IL-Bh host range (Table 2).
There are several modes of comovirus transmission that were reported, including mechanical means [2], beetles [7], and seed transmission [8].In the current study, we have shown that ArLV1-IL-Bh can be mechanically transmitted between plants, but furthermore, we have identified ArLV1-IL-Bh in managed honeybees and three wild bee genera, which are phylogenetically distant [24].A whole-genome coverage of ArLV1-IL-Bh was demonstrated in all four bee libraries (Figure 4).
The reads mapped to ArLV1-IL-Bh encompassed both RNA1 and RNA2 in the four tested bee genera (Figure 4).Among the bees, the mapped reads showed a similar percentage of reads that were mapped to the ArLV1-IL-Bh reference genome, with 0.4-0.44%and 0.34-0.37% of the reads being mapped to RNA1 and RNA2, respectively (Table 3).The minor differences in the abundance of the two segments in the bee libraries (Figure 4, Table 3) were previously reported for numerous A. thaliana accessions hosting ArLV1 [8].This phenomenon may result from separate encapsidations of the bipartite genome in the host plants, which are reported to be host-specific in other systems [52][53][54].
Low percentage values of the mapped reads are in accordance with a non-replicating virus in a non-host vector where the viral particles externally adhered to the bees [10,12,13].Furthermore, the percentage of viral RNA reads did not seem to co-vary with the tendencies of the studied genera's foraging behavior on brassica plants (Figure 5).The low virus loads of ArLV1-IL-Bh are in accordance with passive adherence via floral visitation, suggesting that the quantity of externally adhering virus particles saturates rapidly at lower visitation frequencies compared to the proportion of brassica flowers in the field.In addition, this study focused on bees that significantly differ in their morphological traits, including body size, tongue length, and relative hair cover density, as well as foraging behavior, different pollen collection mechanisms, nesting habits, and more [24].These differences might affect the efficiency of plant virus vectoring by bees.These results could also suggest that the comovirus ArLV1-IL-Bh was shared between the tested genera in the specific geographic region, as was previously demonstrated in viruses infecting arthropods [22,23].Alternatively, the shared source of ArLV1-IL-Bh acquisition might be a plant unrelated to the Brassicaceae family that hosts the comovirus.
The managed and wild bees may be potential non-host vectors of ArLV1-IL-Bh and presumably ArLV1 in the wild, implying a role of bee populations in comovirus spread in wild vegetation, as suggested for AmCV in managed honeybees [20].In addition, the seed transmission of ArLV1 [8] suggests viral presence in pollen grains.This transmission mode supports the possibility that bee pollinators would transmit the comoviruses by virus adherence to their body parts, as was observed with other mechanically transmitted viruses [10,12,13].However, further studies are necessary to establish the transmission of ArLV1-IL-Bh via bee foraging behavior.A broad spectrum view of this possible transmission mode of ArLV1-IL-Bh opens a question regarding the significance of the possible contribution of the managed and wild bees to comovirus spread compared to the reported seed transmission of the virus [8].
It has been proposed that a low virulence could allow plant-pathogen mutualistic benefits by conferring tolerance toward abiotic stress or cross-protection against co-infecting pathogenic viruses [55][56][57].This might be the case for ArLV1-IL-Bh, which caused distinct symptoms in the B. rapa var.perviridis mustard plants with no significant impact on the host plant.In the natural habitats, the B. hirta plants infected by ArLV1-IL-Bh harbored a dense population of the aphid B. brassicae but were asymptomatic.Whether ArLV1-IL-Bh protects brassica plants from aphid-transmitted viruses is unknown.

Conclusions
This study described a possible contribution of managed and wild pollinators to the long-term adaptation of plants with a seed-borne latent comovirus.Our results demonstrate the complexity of these interactions when multiple factors are involved in the plant virus's natural spread.Modulating the characteristics of virus-host plant biological interactions could shed light on trade-offs for both plants and viruses in nature.

Figure 1 .
Figure 1.Transmission electron microscopy visualization and genome organization of a comovirus identified in B. hirta plants harboring a dense aphid population.(a) A representative micrograph showing icosahedral particles.(b) A Coomassie brilliant blue stained gel showing predicted ~20 kDa and ~40 kDa CPs.M, molecular size marker; V, virions.(c,d) Sequencing and genome organization of the mustard comovirus comprised of (c) RNA1 and (d) RNA2 genome segments.

Figure 1 .
Figure 1.Transmission electron microscopy visualization and genome organization of a comovirus identified in B. hirta plants harboring a dense aphid population.(a) A representative micrograph

Plants 2024 , 18 Figure 3 .
Figure 3.A graphical genome presentation of the ArLV1-IL-Bh, sequence assembly, and alignment with the four bee species' identified comovirus contigs.(a-c) RNA1 and (d-f) RNA2.(a,d) represent alignments of HTS de novo assembled contigs from wild and managed bees on the genome of ArLV1-IL-Bh.(b,e) represent the genome amplicons generated from B. hitra-infected plants sequenced by Sanger.(c,f) represent the genome size of both RNA particles, where the dark gray color represents the 3' and 5' UTRs, and the black color represents the polyprotein coding region.(g) studied bee genera frontal and lateral view (left and right, respectively).

Figure 3 .
Figure 3.A graphical genome presentation of the ArLV1-IL-Bh, sequence assembly, and alignment with the four bee species' identified comovirus contigs.(a-c) RNA1 and (d-f) RNA2.(a,d) represent

Figure 5 .
Figure 5.The distributions of differences between the proportion of bee visits on brassica flowers (crucifers) and the relative abundance of brassica plants on each sampling day.Positive and negative values indicate attraction to and avoidance of brassica plants, respectively.The dashed line represents visitation rates that are expected under random encounters."•" a black dot represents an outlier value.

Figure 5 .
Figure 5.The distributions of differences between the proportion of bee visits on brassica flowers (crucifers) and the relative abundance of brassica plants on each sampling day.Positive and negative values indicate attraction to and avoidance of brassica plants, respectively.The dashed line represents visitation rates that are expected under random encounters."•" a black dot represents an outlier value.

Figure 6 .
Figure 6.A phylogenetic tree based on Comovirinae subfamily Pro-Pol region.The following viruses and virus isolates, along with the ArLV1-IL-Bh (OR840696) and ArLV1-IL Am (OR840694), were included:

Table 1 .
Primers are used for sequencing and detection of the viral genome.

Table 2 .
A host range analysis for ArLV1-IL-Bh.

Table 3 .
HTS read the results of four sampled bee populations.

Table 4 .
BLAST analysis of HTS de novo assembled contigs derived from managed and wild bees' pooled samples.