The complete genome of a broad bean wilt virus 1 (BBWV1; Fabavirus alphaviciae) isolate (BBWV1-IR) was recovered from an uncultivated Plantago lanceolata plant in Iran by viral metagenomics and validated by RT-PCR. RNA1 (5,779 nucleotides) contains a single ORF encoding replication-associated proteins (Pro-Co, HEL, VPg, Pro, RdRp), while RNA2 (3,414 nucleotides) harbors two overlapping ORFs encoding the large and small coat proteins (LCP and SCP) and two additional proteins (VP47 and VP37). Comparative analyses revealed that BBWV1-IR shared 81–92% and 79.5–83% nucleotide identity in RNA1 and RNA2, respectively, with global isolates. No intragenic recombination was detected; however, reassortment analysis identified three distinct events, including one involving BBWV1-IR, whose RNA1 segment likely originated from Austrian (major) and UK (minor) parental lineages. ORF1 and ORF2a showed substantial variability and high haplotype diversity, with VP37 displaying the greatest nucleotide diversity. Evolutionary analyses indicated that BBWV1 genes were predominantly shaped by negative selection, with essential replication proteins (HEL and Pro) under strong purifying pressure, while VP47 and VP37 experienced more relaxed constraints. A few codons in ORF1 and ORF2a were under episodic positive selection. Phylogenetic analysis clustered 19 non-recombinant isolates into two major clades (A and B), with BBWV1-IR positioned in subclade I of clade A alongside geographically distant isolates, reflecting human-mediated long-distance dispersal. Incongruent clustering of ORF1 and ORF2a in several isolates supports RNA segment reassortment as a key driver of novel variant emergence. Collectively, these findings highlight the roles of mutation, selection, reassortment, and gene flow in shaping BBWV1 evolution, exemplified by the Iranian isolate.

Keywords