Complete genome sequence of fig leaf mottle-associated virus 2

Archives of Virology (2025) 170:72 https://doi.org/10.1007/s00705-025-06262-0 ANNOTATED SEQUENCE RECORD Complete genome sequence of fig leaf mottle‑associated virus 2 Rachelle Bester1,2 · Shannon Goodchild1 · Hans J. Maree1,2 Received: 13 December 2024 / Accepted: 19 January 2025 / Published online: 11 March 2025 © The Author(s) 2025 Abstract The complete genome sequence of fig leaf mottle-associated virus 2 (FLMaV2), a positive-sense single-stranded RNA virus with a genome length of 16,925 nucleotides, was determined via RT-PCR and rapid amplification of cDNA ends (RACE). Its genome organization resembles those of little cherry virus 2 (LChV2) and yam asymptomatic virus 1 (YaV1), both of which are currently classified as members of the genus Ampelovirus but differ from typical ampeloviruses. Phylogenetic analysis based on amino acid sequences of the RNA-dependent RNA polymerase and coat protein of members of the family Closteroviridae showed that FLMaV2 clustered with LChV2 and YaV1. Because of their unusual genome organization, FLMaV2, LChV2, and YaV1 might represent a distinct genus within the family Closteroviridae. The family Closteroviridae includes non-enveloped, filamentous viruses with a 13- to 19-kb positive-sense RNA genome that is either translated directly to make polyproteins or transcribed to make subgenomic RNAs [1]. The family includes seven genera [2], one of which is the genus Ampelovirus, which has 13 classified species (https://ictv. global/t axono my). Ampeloviruses share a conserved genome organization, featuring a two-gene module (ORF1a–ORF1b) encoding replication proteins and a five-gene module (small hydrophobic protein, heat shock protein 70 homologue [HSP70h], ~60-kDa protein, coat protein [CP], and minor coat protein [CPm]) [3, 4]. Fig leaf mottle-associated virus 2 (FLMaV2), first identified in Algeria [5] and later reported in Syria, Egypt, Lebanon, Tunisia, Albania, Turkey, and South Africa [6, 7], has been proposed to be an ampelovirus but remains unclassified due to the lack of a complete genome sequence. Only a partial 5868-nt sequence (FJ473383) is available, and a complete genome sequence is essential for its accurate classification and characterization. In July 2021, fig leaves with mottling and mosaic patterns from three garden trees in the Western Cape were analysed Handling Editor: Elvira Fiallo-Olivé. * Hans J. Maree 1 Department of Genetics, Stellenbosch University, Private Bag X1, Matieland 7602, South Africa 2 Citrus Research International, PO Box 2201, Matieland 7602, South Africa using high-throughput sequencing (HTS) as described previously [7], yielding a 15,206-nt contig with a high degree of sequence similarity to FLMaV2 (FJ473383.1). This draft genome sequence was used to design primers to produce overlapping amplicons (Supplementary Information S1). In 2023, additional fig leaf samples from the same garden were collected, and RNA was extracted [8, 9] and screened for FLMaV2 using primer set 16 (Supplementary Information S1). Complementary DNA (cDNA) was synthesized from total RNA using random hexamer primers (Promega) and Maxima Reverse Transcriptase (Thermo Fisher Scientific), following the manufacturer’s instructions. PCR was then performed using Phusion HF polymerase (New England Biolabs) and genome-specific forward and reverse primers (IDT) (Supplementary Information S1), following the manufacturer’s instructions. DNA samples containing the amplicons of interest were sent to the Central Analytical Facility (CAF) of Stellenbosch University for bidirectional Sanger sequencing using the same primers that were used for amplification. Five of the eight trees were positive for FLMaV2 by RT-PCR using primer set 16, and pairwise comparisons of the Sanger-sequenced amplicons from the different tree samples revealed a 96.45–100% identity to the HTS draft genome sequence. The sample with the most divergent sequence of FLMaV2 compared to the draft genome sequence was selected for genome confirmation. RNA from sample F8D3 was amplified using primer sets 1-19 and Phusion HF polymerase (New England Biolabs) (Supplementary Information S1). The amplicons were sequenced by the Sanger method as described above. To amplify the Vol.:(0123456789) 72 Page 2 of 5 5’ end of the genome, a SMARTer rapid amplification of cDNA ends (RACE) 5’ kit (Takara Bio) was used according to the manufacturer’s instructions. Plasmid DNA from eight clones was sent to CAF at Stellenbosch University for bidirectional sequencing using M13 primers. For the 3’ end, total RNA was polyadenylated using E. coli poly(A) polymerase (New England Biolabs) according to manufacturer’s instructions and used for reverse transcription with an oligo d(T) primer (primer 21, Supplementary Information S1) to generate cDNA. The cDNA was subjected to thermal cycling using Phusion HF polymerase, an oligo d(T) primer, and a specific forward primer (primer 22, Supplementary Information S1). A DNA sample containing the amplicon of interest was sent to CAF at Stellenbosch University for Sanger sequencing using primer 22. The complete genome sequence of isolate F8D3 (16925 nt) was assembled from the overlapping amplicons using CLC Main Workbench 7 (QIAGEN) (PQ727362). ORF Finder (https://www.ncbi.nlm.nih.gov/ orffinder/) was used to predict open reading frames (ORFs). The partial FLMaV2 sequence that had been submitted previously to NCBI GenBank (FJ473383.1) shared 83.85% nucleotide (nt) sequence identity with the corresponding region of the complete FLMaV2 genome sequence from this study (PQ727362), suggesting that the South African FLMaV2 isolate is a divergent variant. This was further supported by analysis of a 1179-bp amplicon spanning the ORFs encoding the HSP70 and HSP90 homologues. Sample F8D3 exhibited the greatest genetic divergence, with 96.62% nt sequence identity to the other four FLMaV2 sequences, which were 100% identical over 1064 nt. F8D3 shared 82.07% nt sequence identity with FJ473383.1, while the other four sequences exhibited 81.13% identity. The FLMaV2 isolate from this study has the same genome organization as FJ473383.1 (Fig. 1). The 5’ UTR is 1005 nt in length and is followed by nine putative ORFs in the positive orientation and a 3’ UTR of 231 nt. This organization resembles that of little cherry virus 2 (LChV2) and yam asymptomatic virus 1 (YaV1), which are classified as ampeloviruses [10, 11] but are distinct from other members of this genus. For FLMaV2, LChV2, and YaV1, the L-Pro and AlkB domains in ORF1a were not detected, and several 3’ ORFs (e.g., p21, p20A, p20B, p4, and p7) are absent. Notably, the conserved ampelovirus five-gene module (p5, Hsp70h, p55, CP, CPm) is disrupted by an additional p21 R. Bester et al. gene before the CP gene in these viruses (Fig. 1). A similar-sized ORF (p18) that was in the -1 reading frame of both YaV1 and LChV2 (Fig. 1) could not be identified in FLMaV2. However an ORF from position 383 to 511 in the (...truncated)