First report of apple rubbery wood virus 1 in apple in China

More than 30 viral and subviral pathogens infect apple (Malus domestica, an important fruit crop in China) trees and rootstocks, posing a threat to its production. With advances in diagnostic technologies, new viruses including apple rubbery wood virus 1 (ARWV-1), apple rubbery wood virus 2 (ARWV-2), apple luteovirus 1 (ALV), and citrus virus A (CiVA) have been detected (Beatriz et al. 2018; Rott et al. 2018; Hu et al. 2021). ARWV-1 (family Phenuiviridae) is a negative-sense single-stranded RNA virus with three RNA segments (large [L], medium [M], and small [S]). It causes apple rubbery wood disease (Rott et al. 2018) and is found in apple rootstocks, causing leaf yellowing and mottle symptoms in Korea (Lim et al. 2018). To determine virus prevalence in apple trees in China, 200 apple leaf and shoot samples were collected from orchards in Hebei (n = 26), Liaoning (n = 40), Shandong (n = 100), Yunnan (n = 25), Shanxi (n = 4) and Inner Mongolia (5) in 2020. Total RNA was extracted from the shoot phloem or leaf tissues (Hu et al., 2015) and subjected to reverse transcription (RT)-PCR to detect apple chlorotic leaf spot virus (ACLSV), apple stem pitting virus (ASPV), apple stem grooving virus (ASGV), apple necrotic mosaic virus (ApNMV), apple scar skin viroid (ASSVd), ARWV-2, ARWV-1, ALVand CiVA using primers specific to respective viruses (Supplementary Table 1). The prevalence of ACLSV, ASPV, ASGV, ApNMV, ASSVd, ARWV-2, ARWV-1, ALV and CiVA was found to be 75.5%, 85.5%, 86.0%, 43.0%, 4.0%, 48.5%, 10.5%, 0% and 0%, respectively (Supplementary Table 2). Among the 21 positive samples for ARWV-1, three, five and 13 samples were from Hebei, Liaoning and Shandong, respectively. Five ARWV-1-positive samples (cultivars Xinhongjiangjun, Xiangfu-1, Xiangfu-2 and Tianhong) showed leaf mosaic symptoms. To confirm the RT-PCR assay, the projected ARWV-1 amplicons from cvs. Xiangfu-1 and Tianhong were cloned into the pMD18-T vector (Takara, Dalian, China), and three clones of each sample were sequenced. BLASTn analyses demonstrated that the sequences (accession nos. MW507810-MW507811) shared 96.9%-98.9% identity withARWV-1 sequences (MH714536, MF062127, and MF062138) in GenBank. An lncRNA library was prepared for high-throughput sequencing (HTS) with the Illumina HiSeq platform using Xiangfu-1 RNA. A total of 71,613,294 reads were obtained. De novo assembly of the reads revealed 135 viral sequence contigs of ACLSV, ASGV, ASPV, ApNMV, ARWV-1, and ARWV-2. The sequences of co