Next Article in Journal
Comparison of the Virulence of Korean African Swine Fever Isolates from Pig Farms during 2019–2021
Next Article in Special Issue
Bayesian Phylogeographic Inference Suggests Japan as the Center for the Origin and Dissemination of Rice Stripe Virus
Previous Article in Journal
Epidemiology and Transmission Dynamics of Infectious Diseases and Control Measures
Previous Article in Special Issue
Epidemiology of Yam Viruses in Guadeloupe: Role of Cropping Practices and Seed-Tuber Supply
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Weed Hosts Represent an Important Reservoir of Turnip Yellows Virus and a Possible Source of Virus Introduction into Oilseed Rape Crop

1
Plant Virus and Vector Interactions, Crop Research Institute, Drnovská 507, Prague 16106, Czech Republic
2
Department of Plant Protection, Faculty of Agrobiology, Food and Natural Resources, Czech University of Life Sciences Prague, Kamýcká 129, Prague 16500, Czech Republic
3
Biomedical Research Center of the Slovak Academy of Sciences, Institute of Virology, Dúbravská Cesta 9, 84505 Bratislava, Slovakia
4
Faculty of Natural Sciences, University of Ss. Cyril and Methodius, Nám. J. Herdu 2, 91701 Trnava, Slovakia
5
Department of Microbiology and Virology, Comenius University in Bratislava, Ilkovičova 6, 84104 Bratislava, Slovakia
*
Authors to whom correspondence should be addressed.
Viruses 2022, 14(11), 2511; https://doi.org/10.3390/v14112511
Submission received: 21 October 2022 / Revised: 8 November 2022 / Accepted: 10 November 2022 / Published: 13 November 2022
(This article belongs to the Special Issue Plant Virus Epidemiology and Control 2022)

Abstract

:
Turnip yellows virus (TuYV) is one of the most important pathogens of oilseed rape worldwide. The virus has a large host range including many crop species (e.g., oilseed rape, pea, chickpea) and weeds from more than twenty plant families. Other than oilseed rape, we detected TuYV in many commonly grown weed species that share the fields and vegetation period together with canola crops in Czech and Slovak Republics. TuYV was detected by reverse-transcription polymerase chain reaction (RT-PCR) in at least 26 species including main crop hosts (oilseed rape), intercrops and weeds such as Amaranthus retroflexus, Atriplex patula (Amaranthaceae), Arctium lappa, Lactuca serriola, Taraxacum officinale, Tripleurospermum inodorum (Asteraceae), Phacelia tanacetifolia (Boraginaceae), Brassica napus, Capsella bursa–pastoris, Descurainia Sophia, Raphanus raphanistrum, Sinapis alba, Sisymbrium officinale, Thlaspi arvense (Brassicaceae), Silene alba, Stellaria media (Caryophyllaceae), Euphorbia helioscopia (Euphorbiaceae), Geranium rotundifolium (Geraniaceae), Lamium purpureum (Lamiaceae), Fumaria officinalis, Papaver rhoeas (Papaveraceae), Veronica persica (Plantaginaceae syn. Scrophulariaceae), Fallopia convolvulus (Polygonaceae), Solanum nigrum (Solanaceae), Urtica dioica (Urticaceae) and Viola arvensis (Violaceae). The detection of TuYV was further confirmed by RT-qPCR as well as Sanger sequencing of the PCR fragments. We discovered four new weed species as hosts of TuYV such as T. inodorum, S. alba, G. rotundifolium and E. helioscopia, representing their three respective plant families. The readthrough domain (RTD) gene sequence analysis of the Czech and Slovak TuYV isolates from oilseed rape and weed species showed similar within-group nucleotide divergence (7.1% and 5.6%, respectively) and the absence of geographical- or host-based phylogenetic clustering. The high-throughput sequencing of the P. rhoeas sample enabled the obtention of a nearly complete genome of TuYV and revealed the mixed infection of TuYV with turnip mosaic virus and cucumber mosaic virus. Our results thus show that weed species are an important TuYV reservoir and play a significant role in the spread and incidence of the disease in field crops such as oilseed rape.

1. Introduction

Oilseed rape (Brassica napus) is an important crop worldwide and has a long tradition in the European cropping system. Viral diseases are among the limiting biotic factors in oilseed rape cultivation. Four viral pathogens such as turnip yellows virus (TuYV, Polerovirus), turnip mosaic virus (TuMV, Potyvirus), broccoli necrotic yellows virus (BNYV, Cytorhabdovirus) and cauliflower mosaic virus (CaMV, Caulimovirus) are reported to be the main viruses infecting oilseed rape [1,2]. Among them, TuYV is the most important and causes significant losses in oilseed grain production. The high incidence of TuYV in oilseed rape crops could also have important consequences as it serves as a virus source for vegetable crops from Brassicaceae and other host families [3]. Yield reduction in plants showing severe virus symptoms due to TuYV has been estimated at 70 to 79% [3]. TuYV is transmitted by aphids in a persistent manner. The spread of the virus depends on the abundance and movement of the major aphid vectors, namely Myzus persicae and Brevicoryne brassicae [1], as well as virus reservoirs in weeds and volunteer host plants. TuYV has a broad host range, which includes, other than oilseed, a number of weeds [4] as well as economically important vegetable crops and intercrops [5]. Weed species serve not only as virus reservoirs but also as hosts for aphid vectors [6] and promote the spread of the virus to infect new plants and fields [7]. The weeds or cultivated hosts are also of concern for virus epidemiology because they overwinter and thus can constitute inoculum sources for early spring disease outbreaks in oilseed and vegetable crops [8,9,10]. Similarly, in many other cropping systems, weeds are the main reservoir of viruses infecting the field [2,11,12,13] and horticultural crops [14,15,16], especially insect-transmitted viruses [17,18].
TuYV belongs to the genus Polerovirus within the family Solemoviridae [19]. The virus was formerly known as beet western yellows virus (BWYV). TuYV has a single-stranded, positive-sense, approximately 5.6–6 kb RNA genome coding for at last six open reading frames (ORFs) [20]. ORFs 0, 1, and 2 are expressed from genomic RNA and ORFs 3, 4, and 5 are expressed from a subgenomic RNA. ORF 0 encodes a P0 protein that determines disease symptoms and host and acts as a silencing suppressor [21]. The P0 protein interacts with the RNA-induced silencing complex (RISC) and renders the complex inefficient by targeting the degradation of individual Argonaute 1 proteins before the RISC is formed [22]. ORF 1 and ORF 2 overlap and encode a replication-related protein that is individually expressed as a P1–2 fusion protein or a P1 protein [23]. The P2 protein encodes an RNA-dependent RNA polymerase (RdRp) and replicates viral genomic RNA using VPg as a primer template [24]. The major coat protein (CP, P3) is encoded by ORF 3, which is followed in the scaffold by ORF 5 (P5). P3–5 are expressed as a readthrough domain (RTD) protein and are required for aphid transmission and the efficient accumulation and spread of the virus in the plant [25,26,27]. An RTD is also the product of ORF 4 (P4), which is encoded within ORF 3 and is required for viral movement. A small ORF (P3a) is located in an intergenic region between ORF 2 and ORF 3. A viral genomic protein (VPg) is located at the 5′ end of the TuYV genome, as in all poleroviruses, and a 200-nucleotide-long non-coding region is located at the 3′ end [28].
A recent study has shown that TuYV isolates have a large diversity of genome sequences that contain few recombination events. The greatest molecular variability occurs in the RTD gene [10,29], which might be related to the different efficacy of the virus isolates in the aphids’ transmission [26,30].
In this work, we studied TuYV in oilseed rape crops, several weed species, and volunteer plants from different plant families (see Supplementary Table S1). These species are among the most important weeds in oilseed rape crops [31,32] and were present in large numbers in fields throughout the year, forming a “green belt” for the virus and vector. We detected TuYV in oilseed rape, 2 intercrop species and 23 weed species including 4 newly reported hosts such as T. inodorum, S. alba, G. rotundifolium and E. helioscopia. To study the genetic diversity of the virus in the Czech and Slovak republic, sequence analysis of selected isolates targeting the RTD domain of TuYV genome was performed. Furthermore, high-throughput sequencing (HTS) analysis of P. rhoeas with a multiple viral infection enabled us to determine the nearly complete genome of the TuYV isolate.

2. Materials and Methods

2.1. Sampling of Oilseed Rape and Weed Plants from Fields

Sampling was carried out randomly in oilseed rape fields in the vicinity of the Czech and Slovak Republic (Supplementary Table S1). The samples included Brassica napus (oilseed rape), Sinapis alba, Phacelia tanacetifolia (intercrop) and 38 weed species (including the most common for oilseed rape) from 17 plant families.

2.2. Detection of TuYV by RT-PCR

Plant samples were ground in liquid nitrogen, and RNA was extracted from 100 mg of leaf tissue. RNA was isolated using a Trizol kit (Thermo Scientific, Wilmington, DE, USA), according to the manufacturer’s instructions. The concentration and purity of the isolated RNA were then measured spectrophotometrically (NanoDrop 2000; Thermo Scientific, Wilmington, DE, USA). cDNA was prepared using 1.0 μg/μL of total RNA, a random hexamer primer and RevertAid RT (Thermo Scientific, Wilmington, DE, USA), according to the manufacturer’s instructions. The cDNA was amplified with the generic primer pair Luteovirus F/R as well as with the specific primers TuYV_5cpF/TuYVR-K2. The primers are listed in Supplementary Table S2 [33,34]. The 25 μL reaction mixture contained 1 μL cDNA, 1 μL each of 10 μmol forward and reverse primer, 12.5 μL 2× DreamTaq Green PCR Master Mix (Thermo Scientific, Wilmington, DE, USA) and 9.5 μL sterile nuclease-free water. PCR conditions were as follows: initial denaturation at 94 °C for 5 min, followed by 30 amplification cycles (94 °C, 45 s; 60 °C, 30 s; and 72 °C, 1 min), and a final extension at 72 °C for 10 min. The PCR product was then separated on a 1% agarose gel and stained with SYBR® Safe DNA Gel Stain (Invitrogen, Wilmington, DE, USA).

2.3. Analysis of TuYV Titre in Oilseed Rape and Weed Plants by RT-qPCR

Oilseed rape, intercrop and weed samples from different fields were analyzed for TuYV titres using RT-qPCR. RT-qPCR was performed using a LightCycler® 480 (Roche, Basel, Switzerland) with SYBR Green I as described in [35] and Dráb et al. [36]. Absolute quantification of viral RNA copies was performed using cDNA samples with a LightCycler® 480 Instrument II (Roche, Basel, Switzerland) in 384-well plates containing 12 μL of reaction solutions per well. The PCR mastermix in each well consisted of the primer pair mix (0.5 μL of TuYVF-K1 and TuYVR-K1 or TuYVF-K2 and TuYVR-K2 at a final concentration of 10 µM; the primers are listed in Supplementary Table S2), 6 μL of LightCycler® 480 SYBR Green I Master (Roche, Basel, Switzerland), and 4.5 μL of sterile nuclease-free water. Finally, 1 μL cDNA was added to this mixture to obtain a final volume of 12 μL. The thermal cycling protocol was as follows: 42 °C for 30 min, 95 °C for 10 min, followed by 40 cycles at 95 °C for 5 s, 60 °C for 30 s and 20 s at 72 °C. The fluorescence signal was measured at the end of each extension step at 72 °C. After amplification, a melting curve analysis was performed with a temperature gradient from 60 to 97 °C to confirm that only the specific products were amplified. Finally, the samples were cooled to 40 °C for 10 s. The fit points method in the software was used to determine the threshold for the cycle (Ct). TuYV copies in each plant sample were calculated using a standard curve generated by inserting the purified PCR fragment of the target sequence of TuYV (isolate No28dobrovizCZ, maintained in Chinese cabbage in Crop Research Institute) into the vector pGEM-T Easy Vector (Promega, Madison, WI, USA) and cloned into E. coli DH5a. The selected colony with confirmed insertion was cultured and the plasmid DNA was subsequently purified (Plasmid plus Midi Kit, Qiagen, Hilden, Germany). The amount of DNA was quantified by spectrophotometry (NanoDrop 2000; Thermo Scientific). To calculate the TuYV titre in the tested samples, a ten-fold standard dilution curve was constructed using the plasmid DNA of a positive TuYV sample. The plasmid quantity in μg was converted to pmol using the average molecular weight of a ribonucleotide (340 Da) and the number of bases of the transcript (Nb). The following mathematical formula was applied: pmol of dsDNA = μg of dsDNA) × 106 pg μg–1/660 pmol pg−1/Nb. Avogadro’s constant (6.023 × 1023 molecules mol–1) was used to estimate the number of transcripts [27]. Standards and samples were prepared in triplicate. The resulting Ct value of 35 cycles was set as the threshold for virus-positive samples; all samples with a value above 35 were classified as negative.

2.4. Sequence Analysis of Readthrough Domain (RTD) Gene Sequence

A genomic region encompassing the RTD from a group of 26 isolates of oilseed rape, intercrops and weeds was determined from four overlapping fragments amplified using RT-PCR and the primer pairs listed in Table S2. Amplification of cDNA of TuYV isolates was carried out using high-fidelity Phusion High-Fidelity DNA Polymerase (Thermo Scientific, Waltham, MA, USA). The reaction conditions were as follows: initial denaturation at 94 °C for 2 min; followed by 30 amplification cycles of 94 °C for 10 s, 58 °C for 20 s and 72 °C for 30 s; and a final extension at 72 °C for 10 min. The PCR amplicons were purified by using QIAquick PCR Purification Kit according to the manufacture’s instruction (Qiagen, Hilden, Germany). PCR amplicons were commercially sequenced by using forward and reverse primers (Eurofins, Val Fleuri, Luxembourg). Thereafter, the sequences were edited using the bioinformatics tools Reverse Complement (https://www.bioinformatics.org/sms/rev_comp.html; accessed on 1 January 2000) and CAP3 assembly (https://doua.prabi.fr/cgi-bin/run_cap3; accessed on 1 January 2022), and the resulting contigs of nucleotide sequences were analyzed by Blast (NCBI: https://www.ncbi.nlm.nih.gov/; accessed on 1 January 2021). The determined partial TuYV genome sequences were deposited in the GenBank database under accession numbers OP699021-OP699053 (listed in Supplementary Table S3).

2.5. High-Throughput Analysis of the PK2 P. rhoeas Sample

A common poppy (P. rhoeas) plant (labelled as PK2) found in a garden in Pezinok (western Slovakia) in 2017 and showing severe symptoms involving mosaics and leaf deformations was analyzed using HTS. Total RNAs from symptomatic leaves were extracted using the NucleoSpin RNA Plant kit (Macherey-Nagel, Duren, Germany) and ribosomal RNA was removed using the Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, USA). The rRNA-depleted total RNAs were used for double-stranded cDNA synthesis using the SuperScript II (Thermo Fisher Scientific, Waltham, MA, USA) kit. The sample was further processed with the transposon-based chemistry library preparation kit (Nextera XT, Illumina, San Diego, CA, USA) and the cDNA library was sequenced on the Illumina MiSeq platform (Illumina, San Diego, CA, USA, 300-bp paired-end sequencing). High-quality trimmed reads were used for de novo assembly and contigs were blasted to the viral genomes database (ftp://ftp.ncbi.nih.gov/genomes/Viruses/all.fna.tar.gz; accessed on 1 January 2022) using CLC Genomics Workbench 7.5 and Geneious v.8.1.9. The nearly complete TuYV genome sequence reported in this paper (missing 11 and 103 nt at the non-coding 5′ and 3′ extremities as compared to the reference NC_003743 genome) was deposited in the GenBank database under accession number OP562900.

2.6. Multiple Sequence Alignment and Phylogenetic Analysis

Multiple sequence alignment of nucleotide and amino acid sequences was per-formed using Clustal X (version 2.0) [37] with default parameters, and phylogenetic trees were generated using Clustal X and MEGA 7 [38] (www.megasoftware.net/; accessed on 1 January 2022) applying maximum likelihood method. The reliability of branches was inferred from bootstrap analysis of 1000 replicates. The analysis of TuYV diversity was performed using DnaSP software [39].

3. Results and Discussion

3.1. Detection of TuYV in Oilseed Rape and Weed Hosts

We detected TuYV in 160 out of 417 samples tested, using RT-PCR from different species from different plant families, including Amaranthaceae: Amaranthus retroflexus, Atriplex patula; Asteraceae: Arctium lappa, Lactuca serriola, Taraxacum officinale, Tripleurospermum inodorum; Boraginaceae: Phacelia tanacetifolia; Brassicaceae: Brassica napus, Capsella bursa–pastoris, Descurainia Sophia, Raphanus raphanistrum, Sinapis alba, Sisysimbrium officinale, Thlaspi arvense; Caryophyllaceae: Silene alba, Stellaria media; Euphorbiaceae: Euphorbia helioscopia; Geraniaceae: Geranium rotundifolium; Lamiaceae: Lamium purpureum; Papaveraceae: Fumaria officinalis, Papaver rhoeas; Plantaginaceae (Scrophulariaceae): Veronica persica; Polygonaceae: Fallopia convolvulus; Solanaceae: Solanum nigrum; Urticaceae: Urtica dioica; and Violaceae: Viola arvensis (Table S1, Figure 1 and Figure 2). As previously shown, the main hosts of TuYV in arable crops belong to families such as Brassicaceae, Fabaceae, Chenopodiaceae and Asteraceae [5,8,9,40]. Moreover, many weed species from different plant families serve as hosts for the virus, namely Amaranthaceae, Apiaceae, Asteraceae, Boraginaceae, Brassicacae, Caryophyllaceae, Cucurbitaceae, Euphorbiaceae, Fabaceae, Fumarianceae, Geraniaceae, Hydrophyllaceae, Malvaceae, Oxalidaceae, Lamiaceae, Papaveraceae, Polemoniaceae, Portulaceae, Primulaceae, Rhizophoraceae, Scrophulariaceae, Solanacecae, Tropaeolaceae, Urticaceae, and Violaceae [2,5,8,9,28,40,41,42]. However, most of these host data date back from the earlier reports related to the host range of BWYV. Recent data on the newly classified TuYV linked to the weed hosts or intercrop species are still limited. Our results confirm that in addition to oilseed rape and two intercrops (P. tanacetifolia and S. alba), at least 23 weed species serve as hosts for TuYV, including four new host species from Asteraceae (T. inodorum), Caryophyllaceae (S. alba), Geraniaceae (G. rotundifolium), and Euphorbiaceae (E. helioscopia) families. The host range of TuYV thus appears to be very wide among weed species, as it so far includes plant species from 24 families that are commonly grown in and within annual crops worldwide, such as brassicas, fabacious and sugar beet plants. The weed species in which TuYV has been detected in the Czech and Slovak Republics are either winter species that share a common growing season with oilseed rape (Figure 1 and Figure 2) or spring weeds or volunteer plants, which makes them a very effective reservoir for the virus and aphid vectors. The winter species such as V. arvensis, T. arvense, C. bursa-pastoris, D. sophia F. officinalis and G. rotundifolium are the most important weeds for oilseed rape [30,31,43]. The spices such as S. nigrum, A. retroflexus and A. patula are usually spring weeds of broad-row crops (e.g., sugar beet, cruciferous and other species of vegetables) [44] and apparently represent TuYV sources from spring to summer, transmitting the virus to winter weed species as well as crops, mainly by cabbage aphids (Brevicoryne brassicae). The green peach aphids (Myzus persicae) afterward transmit the virus to newly emerging winter crops in early autumn [2,9,45]. Species such as A. lappa, S. alba, and U. dioica mainly grow on the field margins as volunteer plants [31] and are additional reservoirs. In addition, biennial weeds such as A. lappa, V. persica (both annual or biennials) and some perennial spices of Conyza spp. [9] or G. rotundifolium can be a permanent reservoir for TuYV and provide an overwintering green-bridge for the disease cycle. Important reservoir hosts for TuYV are also the species used for intercropping green manures, such as S. alba or P. tanacetifolia, which usually cover the green gap between the post-harvest period and the newly emerging oilseed rape in the fields (Figure 1). The broad and diverse host range of TuYV that includes cultivated, volunteer and weed plants as well as plants grown in winter, spring, and summer, ensures that the virus is constantly maintained in the ecosystem and can be easily transmitted to field crops by aphids [8,46]. Apart from the fact that weeds often grow near or within fields and orchards, the cultivated varieties of some weeds can be used as ornamental plants. Moreover, weed seeds can lie dormant in the soil for a long period and are easy to grow, thus they are abundant in the environment [47] and are important components of disease control in crops.

3.2. TuYV Titre in Oilseed Rape and Weed Plants

RT-qPCR was performed using two primer pairs TuYVF-K1/TuYVR-K1 and TuYVF-K2 and TuYVR-K2. The Ct values were measured in triplicate and plotted against the known copy numbers of the standard sample. The standard curve covered a linear range of five orders of magnitude. For the TuYV K1 primers, the slope was −3.369, intersection = 17.73, r2 = 0.9996 and efficiency = 1.98; for the TuYV K2 primers, the slope was −3.707, intersection = 20.91, r2 = 0.9975 and efficiency = 1.89 (Figure 3A,B) (2.0 being a theoretical perfect efficiency, indicating the doubling of product with every cycle). The laboratory isolate of TuYV was used for the qPCR optimization and both primers in the RT-qPCR assay successfully amplified their respective targets (see Supplementary Table S2). The copies of transcripts are equivalent to the viral gene copies.
The RT-qPCR assay using the primer pair TuYVF-K1/TuYVR-K1 was used to detect the viral titre to confirm the RT-PCR-based detection in the field samples. The detection of the virus titre by RT-qPCR ranged from a few hundred to thousands of virus copies. In most of the weeds, the virus titre remained low, with the exception of A. lappa, D. sophia and V. persica, where titres similar to B. napus were detected (Figure 4).

3.3. The Diversity of TuYV Sequences from Crop and Non-Crop Host

To assess the genetic diversity of TuYV in the Czech and Slovak regions, we analyzed 26 TuYV isolates in a ca. 2040 bp-long genome portion encoding the P3–P5 readthrough domain (RTD), which is required for virus circulation, accumulation and persistence in the aphid vector [24,25,26]. The phylogenetic analysis revealed that TuYV isolates, including those available in Genbank, were split into three major groups (Figure 5). Most of the Czech and Slovak clustered in Group I, irrespective of their host origin. Only two Czech isolates from B. napus and two isolates from P. rhoeas were assigned to Group II, while no isolate was found belonging to Group III. Based on the phylogenetic analyses, no clear-cut clustering was noted among TuYV isolates in the examined genomic region. The calculation of the within mean group nucleotide distance (p-distance model) between Czech and Slovak B. napus (n = 11) and weed isolates (n = 15) showed similar values of 7.1% and 5.6%, respectively. The similar diversity within both populations and the same between-group distance (6.2%) may indicate the absence of particular selection pressure imposed on “culture” and “wild” TuYV isolates and the presumed mixing of both populations during the vegetation period due to vector transmission.

3.4. Virome Analysis of a Common Poppy (P. rhoeas) Plant

In an effort to further compare the genetic relatedness of a “weed” TuYV isolate to known TuYV diversity, the nearly complete genome of the PK2 isolate from additional common poppy (P. rhoeas) was determined from HTS data (available in Genbank under accession number OP562900). The HTS approach provides an unbiased access to the plant virome and viral diversity studies, overcoming the limits of conventional detection tools [48]. The analysis of the generated PK2 sequence data (ca. 3.9 millions of reads with a mean length of 144.8 nts) enabled us to obtain a large contig of ca. 5.3 kb, corresponding to TuYV. Subsequent rounds of mapping of individual reads against the TuYV genomes allowed the validation and reconstruction of the nearly complete genome (5477 nts, 677 reads mapped, with a mean coverage of 17.1×). Interestingly, the phylogenetic analysis placed this isolate in the same cluster as Pisum sativum isolates from the UK and Germany (Figure 6), further confirming the absence of host-related phylogenetic groups within TuYV. Similarly, the BLAST searches identified the P. sativum isolates OK030792, MN497803 and OK030785 as the closest PK2 relatives with more than 98% nt identity. Interestingly, RTD- (Figure 5) and nearly complete genome-based phylogenetic trees (Figure 6) display slightly different topology, probably due to recombination event(s) in the TuYV evolutionary history affecting different genome portions [28]. HTS revealed that besides TuYV, the PK2 plant was co-infected by two additional viruses, i.e., cucumber mosaic virus and turnip mosaic virus. Consequently, due to the mixed infection, the potential role of individual viruses in the observed symptomatology (Supplementary Figure S1) could not be determined. Complex viral infection of plants is a common phenomenon [49], possibly further enhancing the viral disease impact and complicating the phytosanitary management efforts.

4. Conclusions

We detected TuYV in oilseed rape fields in both crops and in various weeds, including winter and spring species as well as volunteer plants. These species are important annual weeds both winter and spring (or some species can be bi-annual e.g., A. lappa, T. inodorum, and perennial, e.g., S. media, G. rotundifolium) in central Europe. Some weed species share a common growing season with oilseed rape, and others are spring weeds and therefore serve as reservoirs for the virus and the aphid vectors (e.g., Brevicoryne brassicae, Myzus persicae) during both pre- and post-harvest periods throughout the vegetation. Similarly, intercrop species (S. alba or P. tanacetifolia) could play a particularly important role in the disease cycle of TuYV as they harbor both the virus and the aphid vectors and close the green gap between crops. A molecular analysis of the RTD coding region showed that there is no strict regional or host-specific clustering of TuYV isolates, suggesting a long-term diversification of the virus, easy gene flow between crop and wild plant communities, and inter-regional spread of the virus through infected plant material. As the virome analysis of the P. rhoeas plant shows, the infection of weeds can be complex and involve different viruses. Due to the reliance on insecticidal seed treatments, the occurrence of the virus becomes a threat to oilseeds or other host plants (e.g., Brassica spps., Pisum vulgaris, Cicer arietinum Lupinus albus Vicia faba, Lens culinaris, Spinacea oleracea, Lactuca sativa). Therefore, weed control as well as the control of aphids in intercrops can significantly reduce disease pressure and the incidence of TuYV in field crops, including oilseed rape.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14112511/s1, Table S1: The plant sample tested for turnip yellows virus; Table S2. Primers used in this study; Table S3: TuYV isolates from the Czech and Slovak Republics and reference isolates of the readthrough domain (RTD) gene sequence analyzed in this study; Figure S1: Common poppy PK2 plant naturally mixed-infected with TuYV, CMV and TuMV.

Author Contributions

Conceptualization, J.K.K. and L.S. methodology, L.S., G.A., E.I., J.T. and K.Š.; molecular analysis: L.S., E.I., G.A., M.G., K.Š. and J.K.K.; bioinformatics analysis, M.G.; writing—original draft preparation, L.S., J.K.K. and M.G.; writing—review and editing, J.K.K. and M.G.; visualization, L.S., E.I. and M.G.; supervision, M.G. and J.K.K.; project administration, J.K.K.; funding acquisition, J.K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the grant from the Technology Agency of the Czech Republic from project TM01000044. M.G. was supported by grant APVV-18-0005 from the Slovak Research and Development Agency.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data used in this study are publicly available on NCBI. A list of the accession numbers used is found in Supplementary Table S3.

Acknowledgments

We thank Markéta Vítámvásová and Hana Novotná (Plant Virus and Vector Interactions, Crop Research Institute, Czech Republic) for their excellent technical assistance. We also thank Jan Mikulka (Crop Research Institute, Czech Republic) for his help with the growing season determination of weed species.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Walsh, J.A.; Tomlinson, J.A. Viruses infecting winter oilseed rape (Brassica-napus ssp oleifera). Ann. Appl. Biol. 1985, 107, 485–495. [Google Scholar] [CrossRef]
  2. Jones, R.A.C.; Sharman, M.; Trebicki, P.; Maina, S.; Congdon, B.S. Virus diseases of cereal and oilseed crops in Australia: Current position and future challenges. Viruses 2021, 13, 2051. [Google Scholar] [CrossRef] [PubMed]
  3. Hardwick, N.V.; Davies, J.M.L.; Wright, D.M. The incidence of three diseases of winter oilseed rape in England and Wales in the 1991/92 and 1992/93 growing season. Plant Pathol. 1994, 43, 1045–1049. [Google Scholar] [CrossRef]
  4. Duffus, J.E. Host relationship of beet western yellows virus strains. Phytopathology 1964, 54, 736–738. [Google Scholar]
  5. Graichen, K.; Rabenstein, F. European isolates of Beet western yellows virus (BWYV) from oilseed rape (Brassica napus L. ssp. napus) are non-pathogenic on sugar beet (Beta vulgaris L. var. altissima) but represent isolates of Turnip yellows virus (TuYV). J. Plant Dis. Prot. 1996, 103, 233–245. [Google Scholar]
  6. Wallis, R.L. Some host plants of the green peach aphid and beet western yellows virus in the Pacific Northwest. J. Econ. Entomol. 1967, 60, 904–907. [Google Scholar] [CrossRef]
  7. Goyal, G.; Harsimran, K.G.; McSorley, R. Common Weed Hosts of Insect-Transmitted Viruses of Florida Vegetable Crops; ENY-863; University of Florida Institute of Food and Agricultural Sciences: Gainesville, FL, USA, 2022; Available online: https://edis.ifas.ufl.edu/publication/in931 (accessed on 29 April 2022).
  8. Stevens, M.; Smith, H.G.; Hallsworth, P.B. The host-range of Beet yellowing viruses among common arable weed species. Plant Pathol. 1994, 43, 579–588. [Google Scholar] [CrossRef]
  9. Stevens, M.; McGrann, G.; Clark, B.; Authority, H. Turnip yellows virus (syn Beet western yellows virus): An emerging threat to European oilseed rape production? HGCA Res. Rev. 2008, 69, 1–37. [Google Scholar]
  10. Freeman, A.; Aftab, M. Effective management of viruses in pulse crops in south eastern Australia should include management of weeds. Australas. Plant Pathol. 2011, 40, 430–441. [Google Scholar] [CrossRef]
  11. Chalupníková, J.; Kundu, J.K.; Singh, K.; Bartakova, P.; Beoni, E. Wheat streak mosaic virus: Incidence in field crops, potential reservoir within grass species and uptake in winter wheat cultivars. J. Integr. Agric. 2017, 16, 523–531. [Google Scholar] [CrossRef]
  12. Rashidi, M.; Cruzado, K.R.; Hutchinson, S.J.P.; Bosque-Pérez, A.N.; Marshall, M.J.; Rashed, A. Grassy weeds and corn as potential sources of barley yellow dwarf virus spread into winter wheat. Plant Dis. 2021, 105, 444–449. [Google Scholar] [CrossRef] [PubMed]
  13. Yazdkhasti, E.; Hopkins, J.R.; Kvarnheden, A. Reservoirs of plant virus disease: Occurrence of wheat dwarf virus and barley/cereal yellow dwarf viruses in Sweden. Plant Pathol. 2021, 70, 1552–1561. [Google Scholar] [CrossRef]
  14. Souza, T.A.; Macedo, M.A.; Albuquerque, L.C.; Inoue-Nagata, K.A. Host range and natural infection of tomato chlorosis virus in weeds collected in Central Brazil. Trop. Plant Pathol. 2020, 45, 84–90. [Google Scholar] [CrossRef]
  15. Korbecka-Glinka, G.; Przyby´s, M.; Feledyn-Szewczyk, B.A. Survey of five plant viruses in weeds and tobacco in Poland. Agronomy 2021, 11, 1667. [Google Scholar] [CrossRef]
  16. Hancinsky, R.; Mihalik, D.; Mrkvova, M.; Candresse, T.; Glasa, M. Plant viruses infecting Solanaceae family members in the cultivated and wild environments: A review. Plants 2020, 9, 667. [Google Scholar] [CrossRef] [PubMed]
  17. Maliano, M.R.; Macedo, M.A.; Rojas, M.R.; Gilbertson, R.L. Weed-infecting viruses in a tropical agroecosystem present different threats to crops and evolutionary histories. PLoS ONE 2021, 16, e0250066. [Google Scholar] [CrossRef]
  18. Chen, G.; Pan, H.; Xie, W.; Wang, S.; Wu, Q.; Fang, Y.; Shi, X.; Zhang, Y. Virus infection of a weed increases vector attraction to and vector fitness on the weed. Sci. Rep. 2013, 3, 2253. [Google Scholar] [CrossRef] [Green Version]
  19. Sõmera, M.; Fargette, D.; Hébrard, E.; Sarmiento, C. ICTV Report Consortium: ICTV Virus Taxonomy Profile: Solemoviridae. J. Gen. Virol. 2021, 102, 001707. [Google Scholar]
  20. Stevens, M.; Freeman, B.; Liu, H.Y.; Herrbach, E.; Lemaire, O. Beet poleroviruses: Close friends or distant relatives? Mol. Plant Pathol. 2005, 6, 1–9. [Google Scholar] [CrossRef] [Green Version]
  21. Pfeffer, S.; Dunoyer, P.; Heim, F.; Richards, K.E.; Jonard, G.; Ziegler-Graff, V. P0 of Beet western yellows virus is a suppressor of posttranscriptional gene silencing. J. Virol. 2002, 76, 6815–6824. [Google Scholar] [CrossRef] [Green Version]
  22. Bortolamiol, D.; Pazhouhandeh, M.; Marrocco, K.; Genschik, P.; Ziegler-Graff, V. The Polerovirus F box protein P0 targets ARGONAUTE1 to suppress RNA silencing. Curr. Biol. 2007, 17, 1615–1621. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Reutenauer, A.; Ziegler-Graff, V.; Lot, H.; Scheidecker, D.; Guilley, H.; Richards, K.; Jonard, G. Identification of Beet western yellows virus luteovirus genes implicated in viral replication and particle morphogenesis. Virology 1993, 195, 692–699. [Google Scholar] [CrossRef] [PubMed]
  24. Martin, R.R.; Keese, P.K.; Young, M.J.; Waterhouse, P.M.; Gerlach, W.L. Evolution and molecular biology of luteoviruses. Annu. Rev. Phytopathol. 1990, 28, 341–363. [Google Scholar] [CrossRef]
  25. Brault, V.; van den Heuvel, J.F.; Verbeek, M.; Ziegler-Graff, V.; Reutenauer, A.; Herrbach, E.; Garaud, J.C.; Guilley, H.; Richards, K.; Jonard, G. Aphid transmission of beet western yellows luteovirus requires the minor capsid read-through protein P74. EMBO J. 1995, 14, 650–659. [Google Scholar] [CrossRef]
  26. Van den Heuvel, J.F.J.M.; Bruyère, A.; Hogenhout, A.; Ziegler-Graff, V.; Brault, V.; Verbeek, M.; van der Wilk, F.; Richards, K. The N-terminal region of the luteovirus readthrough domain determines virus binding to Buchnera GroEL and is essential for virus persistence in the aphid. J. Virol. 1997, 71, 7258–7265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Peter, K.A.; Gildow, F.; Palukaitis, P.; Gray, S.M. The C terminus of the polerovirus P5 readthrough domain limits virus infection to the phloem. J. Virol. 2009, 83, 5419–5429. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Smirnova, E.; Firth, A.E.; Miller, W.A.; Scheidecker, D.; Brault, V.; Reinbold, C.; Rakotondrafara, A.M.; Chung, B.Y.W.; Ziegler-Graff, V. Discovery of a small non-AUG-initiated ORF in poleroviruses and luteoviruses that is required for long-distance movement. PLoS Pathog. 2015, 11, e1004868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Filardo, F.; Nancarrow, N.; Kehoe, M.; McTaggart, A.R.; Congdon, B.; Kumari, S.; Aftab, M.; Trębicki, P.; Rodoni, B.; Thomas, J.; et al. Genetic diversity and recombination between turnip yellows virus strains in Australia. Arch. Virol. 2021, 166, 813–829. [Google Scholar] [CrossRef]
  30. Brault, V.; Périgon, S.; Reinbold, C.; Erdinger, M.; Scheidecker, D.; Herrbach, E.; Richards, K.; Ziegler-Graff, V. The polerovirus minor capsid protein determines vector specificity and intestinal tropism in the aphid. J. Virol. 2005, 79, 9685–9693. [Google Scholar] [CrossRef] [Green Version]
  31. Hanzlik, K.; Gerowitt, B. Occurrence and distribution of important weed species in German winter oilseed rape fields. J. Plant Dis. Prot. 2012, 119, 107–120. [Google Scholar] [CrossRef]
  32. Krähmer, H.; Andreasen, C.; Economou-Antonaka, G.; Holec, J.; Kalivas, D.; Kolářová, M.; Novák, R.; Panozzog, S.; Pinke, G.; Salonen, J.; et al. Weed surveys and weed mapping in Europe: State of the art and future tasks. Crop Prot. 2020, 129, 105010. [Google Scholar] [CrossRef]
  33. Abraham, A.D.; Menzel, W.; Lesemann, D.E.; Varrelmann, M.; Vetten, H.J. Chickpea chlorotic stunt virus: A new polerovirus infecting cool-season food legumes in ethiopia. Phytopathology 2006, 96, 437–446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Singh, K.; Slavíková, L.; Kumar, J. Reakční Směs pro Kvantifikaci Viru Žloutenky Vodnice v Kulturních a Nekulturních Vzorcích Pomocí qPCR (Reaction Mixture for Quantification of Turnip Yellows Virus in Cultivated and Uncultivated Plant Samples by qPCR); Užitný Vzor (Utility Model) no. 35269; Industrial Property Office: Prague, Czech Republic, 2021. [Google Scholar]
  35. Jarošová, J.; Kundu, J.K. Validation of reference genes as internal control for studying viral infections in cereals by quantitative real-time RT-PCR. BMC Plant Biol. 2010, 10, 146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Dráb, T.; Svobodová, E.; Ripl, J.; Jarošová, J.; Rabenstein, F.; Melcher, U.; Kundu, J.K. SYBR Green I based RT-qPCR assays for the detection of RNA viruses of cereals and grasses. Crop Pasture Sci. 2014, 65, 1323–1328. [Google Scholar] [CrossRef]
  37. Larkin, M.A.; Blackshields, G.; Brown, N.P.; Chenna, R.; McGettigan, P.A.; McWilliam, H.; Valentin, F.; Wallace, I.M.; Wilm, A.; Lopez, R.; et al. Clustal W and Clustal X version 2.0. Bioinformatics 2007, 23, 2947–2948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  39. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  40. Walkey, D.G.A.; Pink, D.A.C. Studies on resistance to Beet western yellows virus in lettuce (Lactuca sativa) and the occurrence of field sources of the virus. Plant Pathol. 1990, 39, 141–155. [Google Scholar] [CrossRef]
  41. Pallett, D.W.; Thurston, M.I.; Cortina-Borja, M.; Edwards, M.L.; Alexander, M.; Mitchell, E.; Raybould, A.F.; Cooper, J.I. The incidence of viruses in wild Brassica rapa ssp. sylvestris in southern England. Ann. Appl. Biol. 2002, 141, 163–170. [Google Scholar] [CrossRef]
  42. Thurston, M.I.; Pallett, D.W.; Cortina-Borja, M.; Edwards, M.L.; Raybould, A.F.; Cooper, J.I. The incidence of viruses in wild Brassica nigra in Dorset (UK). Ann. Appl. Biol. 2001, 139, 277–284. [Google Scholar] [CrossRef]
  43. Coutts, B.A.; Hawkes, J.R.; Jones, R.A.C. Occurrence of Beet western yellows virus and its aphid vectors in over-summering broad-leafed weeds and volunteer crop plants in the grainbelt region of south-western Australia. Aust. J. Agric. Res. 2006, 57, 975–982. [Google Scholar] [CrossRef]
  44. Kolářová, M.; Tyšer, J.; Soukup, J. Survey about the weed occurrence on arable land in the Czech Republic. Sci. Agric. Bohem. 2013, 44, 63–69. [Google Scholar] [CrossRef]
  45. Jursík, M.; Holec, J.; Soukup, J.; Venclová, V. Competitive relationships between sugar beet and weeds in dependence on time of weed control. Plant Soil Environ. 2008, 54, 108–116. [Google Scholar] [CrossRef] [Green Version]
  46. Latham, L.J.; Smith, L.J.; Jones, R.A.C. Incidence of three viruses in vegetable brassica plantings and associated wild radish weeds in south-west Australia. Australas. Plant Pathol. 2003, 32, 387–391. [Google Scholar] [CrossRef]
  47. Primot, S.; Valantin-Morison, M.; Makowski, D. Predicting the risk of weed infestation in winter oilseed rape crops. Weed Res. 2005, 46, 22–33. [Google Scholar] [CrossRef]
  48. Maree, H.J.; Fox, A.; Al Rwahnih, M.; Boonham, N.; Candresse, T. Application of HTS for routine plant virus diagnostics: State of the art and challenges. Front. Plant Sci. 2018, 9, 1082. [Google Scholar] [CrossRef] [Green Version]
  49. Mrkvová, M.; Hančinský, R.; Predajňa, L.; Alaxin, P.; Achs, A.; Tomašechová, J.; Šoltys, K.; Mihálik, D.; Olmos, A.; Ruiz-García, A.B.; et al. High-throughput sequencing discloses the Cucumber mosaic virus (CMV) diversity in Slovakia and reveals new hosts of CMV from the Papaveraceae family. Plants 2022, 11, 1665. [Google Scholar] [CrossRef]
Figure 1. The growing season of oilseed rape and weeds covers the whole year. Most species are annual weeds, with the exception of A. lappa, T. arvense (* either annual or bi-annual) and S. media, G. rotundifolium (** either annual or perennial). (The blue color indicates the beginning of the growing season).
Figure 1. The growing season of oilseed rape and weeds covers the whole year. Most species are annual weeds, with the exception of A. lappa, T. arvense (* either annual or bi-annual) and S. media, G. rotundifolium (** either annual or perennial). (The blue color indicates the beginning of the growing season).
Viruses 14 02511 g001
Figure 2. Oilseed rape, intercrops and weed species infected by TuYV. (13): Brassica napus, (4): Arctium lappa, (5): Capsella bursa-pastoris, (6): Descurainia sophia, (7): Euphorbia helioscopia, (8): Fumaria officinalis, (9): Geranium rotundifolium, (10): Lamium purpureum, (11): Papaver rhoeas, (12): Phacelia tanacetifolia (intercrop), (13): Sinapis alba (intercrop), (14): Sisymbrium officinale, (15): Stellaria media, (16): Thlaspi arvense, (17): Tripleurospermum inodoorum, (18): Veronica persica, (19): Viola arvensis, (20): Lactuca seriola, (21): Amaranthus retroflexus, (22): Atriplex patula, (23): Silene alba, (24): Solanum nigrum, (25): Taraxacum officinalis, (26): Urtica dioica.
Figure 2. Oilseed rape, intercrops and weed species infected by TuYV. (13): Brassica napus, (4): Arctium lappa, (5): Capsella bursa-pastoris, (6): Descurainia sophia, (7): Euphorbia helioscopia, (8): Fumaria officinalis, (9): Geranium rotundifolium, (10): Lamium purpureum, (11): Papaver rhoeas, (12): Phacelia tanacetifolia (intercrop), (13): Sinapis alba (intercrop), (14): Sisymbrium officinale, (15): Stellaria media, (16): Thlaspi arvense, (17): Tripleurospermum inodoorum, (18): Veronica persica, (19): Viola arvensis, (20): Lactuca seriola, (21): Amaranthus retroflexus, (22): Atriplex patula, (23): Silene alba, (24): Solanum nigrum, (25): Taraxacum officinalis, (26): Urtica dioica.
Viruses 14 02511 g002
Figure 3. Standard curves for SYBR Green-I-based real-time RT-PCR amplification of standard TuYV DNA with specific primer sets (A) TuYV K1 and (B) TuYV K2. The amplification plots show the testing, in triplicate, of a ten-fold dilution series containing standard TuYV DNA from 2.1 × 102 to 2.1 × 104 template copies/reaction.
Figure 3. Standard curves for SYBR Green-I-based real-time RT-PCR amplification of standard TuYV DNA with specific primer sets (A) TuYV K1 and (B) TuYV K2. The amplification plots show the testing, in triplicate, of a ten-fold dilution series containing standard TuYV DNA from 2.1 × 102 to 2.1 × 104 template copies/reaction.
Viruses 14 02511 g003
Figure 4. (A). TuYV titre detected in RT-PCR-positive plants of oilseed rape, intercrop and weed species by RT-qPCR. (B). TuYV detection by RT-PCR using the primer pair TuYV_5cpF/TuYVR-K2, 797 bp (lanes 1–33, representing the same samples as in Figure A {from left to right}; lanes 34 and 35 positive and negative control, respectively; lane M GeneRuler 1 kb DNA ladder, Thermo fisher Scientific, Wilmington, DE, USA).
Figure 4. (A). TuYV titre detected in RT-PCR-positive plants of oilseed rape, intercrop and weed species by RT-qPCR. (B). TuYV detection by RT-PCR using the primer pair TuYV_5cpF/TuYVR-K2, 797 bp (lanes 1–33, representing the same samples as in Figure A {from left to right}; lanes 34 and 35 positive and negative control, respectively; lane M GeneRuler 1 kb DNA ladder, Thermo fisher Scientific, Wilmington, DE, USA).
Viruses 14 02511 g004
Figure 5. Maximum likelihood phylogenetic tree constructed from alignments of partial TuYV sequences of the genomic regions encompassing the RTD domain (nts 3483–5508 based on the reference NC_003743 genome). The Czech, Slovak and selected Genbank sequences of previously characterized isolates are identified by their accession number, isolate name and geographical location. Only bootstrap values ≥ 70% (1000 bootstrap resamplings) are indicated.
Figure 5. Maximum likelihood phylogenetic tree constructed from alignments of partial TuYV sequences of the genomic regions encompassing the RTD domain (nts 3483–5508 based on the reference NC_003743 genome). The Czech, Slovak and selected Genbank sequences of previously characterized isolates are identified by their accession number, isolate name and geographical location. Only bootstrap values ≥ 70% (1000 bootstrap resamplings) are indicated.
Viruses 14 02511 g005
Figure 6. (A). Phylogenetic tree generated on available nearly complete nucleotide genome sequences of TuYV isolates (n = 65). The phylogenetic analysis was inferred using maximum likelihood (ML) based on the General Time Reversible (GTR) model selected as the best-fit model of nucleotide substitution based on Bayesian information criterion (BIC) as implemented in MEGA X. Isolates are identified by their GenBank accession numbers and host (if known). The Slovak PK2 isolate sequenced in the present study is highlighted by an arrow. The scale bar indicates a genetic distance of 0.01. Bootstrap values higher than 70% (1000 bootstrap resamplings) are indicated. (B). Nucleotide diversity (π) plotted along the TuYV genomes with a window of 100 nt and a step of 25 nt. The blue dotted line indicates the mean nucleotide diversity (π = 9.2%) within TuYV genomes. In scale schematic representation of TuYV genome is shown below.
Figure 6. (A). Phylogenetic tree generated on available nearly complete nucleotide genome sequences of TuYV isolates (n = 65). The phylogenetic analysis was inferred using maximum likelihood (ML) based on the General Time Reversible (GTR) model selected as the best-fit model of nucleotide substitution based on Bayesian information criterion (BIC) as implemented in MEGA X. Isolates are identified by their GenBank accession numbers and host (if known). The Slovak PK2 isolate sequenced in the present study is highlighted by an arrow. The scale bar indicates a genetic distance of 0.01. Bootstrap values higher than 70% (1000 bootstrap resamplings) are indicated. (B). Nucleotide diversity (π) plotted along the TuYV genomes with a window of 100 nt and a step of 25 nt. The blue dotted line indicates the mean nucleotide diversity (π = 9.2%) within TuYV genomes. In scale schematic representation of TuYV genome is shown below.
Viruses 14 02511 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Slavíková, L.; Ibrahim, E.; Alquicer, G.; Tomašechová, J.; Šoltys, K.; Glasa, M.; Kundu, J.K. Weed Hosts Represent an Important Reservoir of Turnip Yellows Virus and a Possible Source of Virus Introduction into Oilseed Rape Crop. Viruses 2022, 14, 2511. https://doi.org/10.3390/v14112511

AMA Style

Slavíková L, Ibrahim E, Alquicer G, Tomašechová J, Šoltys K, Glasa M, Kundu JK. Weed Hosts Represent an Important Reservoir of Turnip Yellows Virus and a Possible Source of Virus Introduction into Oilseed Rape Crop. Viruses. 2022; 14(11):2511. https://doi.org/10.3390/v14112511

Chicago/Turabian Style

Slavíková, Lucie, Emad Ibrahim, Glenda Alquicer, Jana Tomašechová, Katarína Šoltys, Miroslav Glasa, and Jiban Kumar Kundu. 2022. "Weed Hosts Represent an Important Reservoir of Turnip Yellows Virus and a Possible Source of Virus Introduction into Oilseed Rape Crop" Viruses 14, no. 11: 2511. https://doi.org/10.3390/v14112511

APA Style

Slavíková, L., Ibrahim, E., Alquicer, G., Tomašechová, J., Šoltys, K., Glasa, M., & Kundu, J. K. (2022). Weed Hosts Represent an Important Reservoir of Turnip Yellows Virus and a Possible Source of Virus Introduction into Oilseed Rape Crop. Viruses, 14(11), 2511. https://doi.org/10.3390/v14112511

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

Article Metrics

Back to TopTop