Lily plants exhibiting virus-like symptoms of leaf yellowing, twisting and brownish necrotic spots were collected, and next-generation sequencing of small RNAs was conducted to identify the associated viruses. Cucumber mosaic virus, lily symptomless virus and a hitherto unrecorded potyvirus, tentatively named "lily yellow mosaic virus" (LYMV), were detected. The genomic RNA of LYMV was 9811 nt in length, encoding a large polyprotein of 3,124 amino acids with a predicted M of 353.3 kDa. BLAST analysis showed that LYMV shared a high degree of amino acid sequence identity with Thunberg fritillary mosaic virus (55%), bean yellow mosaic virus (52%), clover yellow vein virus (51%), leek yellow stripe virus (51%), and lily mottle virus (52%), and these viruses clustered together in a phylogenetic tree.

To identify the viruses in tree peony plants associated with the symptoms of yellowing, leaf rolling, stunted growth, and decline, high-throughput sequencing of small RNA and mRNA was conducted from a single symptomatic plant. Bioinformatic analyses and reconstruction of viral genomes indicated mixed viral infections involving cycas necrotic stunt virus (CNSV), apple stem grooving virus (ASGV), lychnis mottle virus (LycMoV), grapevine line pattern virus (GLPV), and three new viruses designated as peony yellowing-associated citrivirus (PYaCV, Citrivirus in Betaflexiviridae), peony betaflexivirus 1 (PeV1, unclassified in Betaflexiviridae), and peony leafroll-associated virus (PLRaV, Ampelovirus in Closteroviridae). PYaCV was 8,666 nucleaotides (nt) in length, comprising three open reading frames (ORFs) and shared 63.8–75.9% nucleotide sequence identity with citrus leaf blotch virus (CLBV) isolates. However, the ORF encoding the replication-associated protein (REP) shared 57% and 52% sequence identities at the nt and amino acid (aa) level, respectively, with other reported CLBV isolates, which were below the criterion for species classification within the family Betaflexiviridae. Recombination analysis identified putative recombination sites in PYaCV, which originated from CLBV. PeV1, only identified from the transcriptome data, was 8,124 nt in length with five ORFs encoding the REP (ORF1), triple gene block (TGB, ORF2–4) and coat protein (CP, ORF5) proteins. Phylogenetic analysis and sequence comparison showed that PeV1 clustered with an unassigned member, the garlic yellow mosaic-associated virus (GYMaV) within the Betaflexiviridae family, into a separate clade. Partial genome sequence analysis of PLRaV (12,545 nt) showed it contained seven ORFs encoding the partial polyprotein 1a, the RNA-dependent RNA polymerase (RdRp), two small hydrophobic proteins p11 and p6, HSP70h, p55, and a CP duplicate, which shared low aa sequence identity with Closteroviridae family members. Phylogenetic analysis based on the aa sequences of RdRp or HSP70h indicated that PLRaV clustered with grapevine leafroll-associated virus 1 (GLRaV-1) and GLRaV-13 in the Ampelovirus genus. Field investigation confirmed the wide distribution of these viruses, causing mixed infections of peony plants in Beijing.