Guosong Wen

A novel caulimovirus was identified from diseased leaves of Dregea volubilis exhibiting yellowing and vein-associated chlorosis in Yuanjiang County, Yunnan Province, China. The virus was tentatively named Dregea volubilis caulimovirus 1 (DVCaV1). The complete genome sequence of DVCaV1, determined by de novo assembly of high-throughput sequencing data, comprises 8,160 bp of circular double-stranded DNA containing two intergenic regions and seven open reading frames (ORFs). These ORFs encode (in order) a movement protein (MP), an aphid transmission factor (ATF), a virion-associated protein (VAP), a coat protein (CP), a polymerase polyprotein (Pol, containing protease, reverse transcriptase, and RNase H domains), a transactivator/viroplasmin (TAV) protein, and a hypothetical protein of unknown function. Sequence comparisons revealed the highest nucleotide similarity with strawberry vein banding virus (SVBV; NC_001725). Phylogenetic analysis confirmed DVCaV1 as a member of the genus Caulimovirus, with SVBV as its closest known relative. According to current ICTV species demarcation criteria for the genus Caulimovirus (host range and > 20% nucleotide sequence divergence in the polymerase region), DVCaV1 represents a novel species. This is, to our knowledge, the first report of a caulimovirus detected in naturally symptomatic Dregea volubilis.

Polygonati Rhizoma and Polygonati odorati Rhizoma, known as “Huangjing” and “Yuzhu” in China, are medicinal Polygonatum species resources with top-grade medical and edible properties. The chloroplast (cp) genome has been used to study species diversity, evolution, and breeding of species for applications in genetic engineering. Codon usage bias (CUB), a common and complex natural phenomenon, is essential for studies of codon optimization of exogenous genes, genetic engineering, and molecular evolution. However, the CUB of medicinal Polygonatum species chloroplast genomes has not been systematically studied. In our study, a detailed analysis of CUB was performed in the medicinal Polygonatum species chloroplast genomes. We investigated the codon bias of 204 plastid protein-coding genes (PCGs) in 4 medicinal Polygonatum species using CodonW and CUSP online software. Through the analysis of the codon bias index, we found that the medicinal Polygonatum species chloroplast genomes had weak codon usage bias. In addition, our results also showed a high preference for AT bases in medicinal Polygonatum species chloroplast genomes, and the preference to use AT-ending codons was observed in these species chloroplast genomes. The neutrality plot, ENC plot, PR2-Bias plot, and correspondence analysis showed that compared with mutation pressure, natural selection was the most important factor of CUB. Based on the comparative analysis of high-frequency codons and high expression codons, we also determined the 10-11 optimal codons of investigative medicinal Polygonatum species. Furthermore, the result of RSCU-based cluster analysis showed that the genetic relationship between different medicinal Polygonatum species could be well reflected. This study provided an essential understanding of CUB and evolution in the medicinal Polygonatum species chloroplast genomes.

The complete genome sequence of a novel badnavirus, tentatively named “fatsia badnavirus 1” (FaBV1, OM540428), was identified in Fatsia japonica. The infected plant displayed virus-like symptoms on leaves, including yellowing and chlorosis. The genome of FaBV1 is 7313 bp in length and similar in size and organization to other members of the genus Badnavirus (family Caulimoviridae), containing four open reading frames (ORFs), three of which are found in all known badnaviruses, and the other of which is only present in some badnaviruses. The virus has the genome characteristics of badnaviruses, including a tRNA^Met binding site (5’-TCTGAATTTATAGCGCTA-3’) and two cysteine-rich domains (C-X-C-2X-C-4X-H-4X-C and C-2X-C-11X-C-2X-C-4X-C-2X-C). Pairwise sequence comparisons of the RT+RNase H region indicated that FaBV1 shares 61.4-71.2% nucleotide (nt) sequence identity with other known badnaviruses, which is below the threshold (80% nt sequence identity in the RT+RNase H region) used for species demarcation in the genus Badnavirus. Phylogenetic analysis revealed that FaBV1, ivy ringspot-associated virus (IRSaV, MN850490.1), and cacao mild mosaic virus (CMMV, KX276640.1) together form a separate clade within the genus Badnavirus, suggesting that FaBV1 is a new member of the genus Badnavirus in the family Caulimoviridae. To our knowledge, this is the first report of a badnavirus infecting F. japonica.

Lily virus X (LVX) is a positive-sense ssRNA virus belonging to the genus Potexvirus in the family Alphaflexiviridae. LVX is known to infect plants of the genera Lilium and Tricyrtis in the family Liliacea. LVX was first reported in an asymptomatic lily (Lilium formosanum) from England (Stone, 1980), but has been shown to infect plants in the Netherlands (Chen et al. 2005), the United States (Jordan et al. 2008) and Japan (Nijo et al. 2018). To date, the complete genomes of two LVX isolates from the Netherlands and Japan have been reported. Paris polyphylla var. yunnanensis, known as Dianchonglou in China, is a perennial plant of the family Melanthiaceae (formerly belonging to the family Trillium). In China, its rhizome is commonly used as an antispasmodic agent for stroke and cancer treatment (Chang et al. 2017). From 2019 to 2022, leaf mottle and shrinkage which are typical symptoms of viral infections were observed on the leaves of P. polyphylla var. yunnanensis plants in Dianchonglou fields in Qujing, Yunnan. Disease incidence ranged from 19% to 45% across 5 fields (90 plants per field) in Qujing. To identify the possible viral pathogen(s) associated with the disease, the mirVanaTM miRNA isolation Kit was used to extract total RNA was from a mixed sample pool of 5 symptomatic leaf samples collected from the 5 fields. RNA sequencing library was constructed using TruSeqTM RNA sample preparation kit. Sequencing on the Illumina HiSeqTM 2500 platform (Illumina, USA) with 125-bp paired-end reads yielded 23,077,786 raw reads. 22,534,100 clean reads were obtained by removing reads of low quality and poly-N using Trimmomatic software (Bolger et al. 2014). By utilizing the paired-end splicing method in Trinity software (Grabherr et al. 2011) the the raw reads were De novo assembled into 184,596 contigs, of which 303 were related to viruses, including Paris mosaic necrosis virus (PMNV), Pear alphapartitivirus (PAPV), Dahlia mosaic virus (DMV), and Lily virus X (LVX). BLASTn analysis revealed that 12 contigs (lengths ranging from 344 nt to 5,981 nt, query cover 6% to 99%) were most similar (57.32% to 91.67% nt identities) to the genome sequences of LVX, suggesting a possible infection of LVX in the plants. To confirm the result, a full-length genomic sequence of LVX was obtained by reverse transcription polymerase chain reaction (RT-PCR) using specific primers designed based on the sequence of the assembled contigs. The PCR products were cloned into pGEM-T vector (Promega Corporation, USA) and sequenced using the Sanger method (Sangon Biotech, Shanghai, China). The obtained full-length genomic sequence of the LVX isolate (LVX-PP, accession number OM100017) was 5,981 nt in length. BLASTp analysis demonstrated that the putative Rep and CP of LVX-PP shared 76.27% to 81.05% and 80.81% to 81.82% aa sequence similarities with that of other LVX isolates, respectively. Maximum-likelihood phylogenetic trees inferred from the Rep and CP aa sequences showed that LVX-PP clustered closely with LVX isolates. The leaf samples were further analyzed using a lily virus X (LVX) ELISA kit (DEIAPV181, Creative Diagnostics, U.S.A.). Healthy P. polyphylla var. yunnanensis leaves were taken as a negative control and buffer solution as a blank control. The results showed a positive reaction for all five symptomatic plants (OD = 1.259 ± 0.007) relative to the negative (OD = 0.099) and blank (OD = 0.073) controls. These results indicate that LVX can infect P. polyphylla var. yunnanensis. To our knowledge, this is the first report that LVX has been detected in P. polyphylla var. yunnannensis. This study will serve as an important reference for the study of the host range of LVX. Further studies will be required to determine how LVX spreads between P. polyphylla var. yunnannensis and other host plants.

Background: Polygonatum kingianum has been widely used as a traditional Chinese medicine as well as a healthy food. Because of its highly variable morphology, this medicinal plant is often difficult to distinguish from other related verticillate leaf types of the Polygonatum species. The contaminants in P. kingianum products not only decrease the products’ quality but also threaten consumer safety, seriously inhibiting the industrial application of P. kingianum.Methods: Nine complete chloroplast (cp) genomes of six verticillate leaf types of the Polygonatum species were de novo assembled and systematically analyzed.Results: The total lengths of newly sequenced cp genomes ranged from 155,437 to 155,977 bp, including 86/87 protein-coding, 38 tRNA, and 8 rRNA genes, which all exhibited well-conserved genomic structures and gene orders. The differences in the IR/SC (inverted repeats/single-copy) boundary regions and simple sequence repeats were detected among the verticillate leaf types of the Polygonatum cp genomes. Comparative cp genomes analyses revealed that a higher similarity was conserved in the IR regions than in the SC regions. In addition, 11 divergent hotspot regions were selected, providing potential molecular markers for the identification of the Polygonatum species with verticillate leaf types. Phylogenetic analysis indicated that, as a super barcode, plastids realized a fast and efficient identification that clearly characterized the relationships within the verticillate leaf types of the Polygonatum species. In brief, our results not only enrich the data on the cp genomes of the genus Polygonatum but also provide references for the P. kingianum germplasm resource protection, herbal cultivation, and drug production.Conclusion: This study not only accurately identifies P. kingianum species, but also provides valuable information for the development of molecular markers and phylogenetic analyses of the Polygonatum species with verticillate leaf types.

The complete genome sequence of a putative novel closterovirus, tentatively named "Dregea volubilis closterovirus 1" (DvCV1, GenBank accession no. MZ779122), infecting Dregea volubilis in China was determined using high-throughput sequencing (HTS). The complete genome sequence of DvCV1 consists of 16,165 nucleotides (nt) and contains nine ORFs. The genome structure of DvCV1 is typical of members of the genus Closterovirus. Complete genome sequence analysis showed that DvCV1 shares 41.4–48.4% nucleotide sequence identity with other known closteroviruses. The putative RNA-dependent RNA polymerase (RdRp), heat shock protein 70-like protein (HSP70h), and coat protein (CP) of DvCV1 share 46.80-62.65%, 31.06–51.80%, and 28.34–37.37% amino acid sequence identity, respectively, with the RdRp, HSP70h and CP of other closteroviruses. Phylogenetic analysis based on HSP70h aa sequences placed DvCV1 alongside other members of the genus Closterovirus in the family Closteroviridae. These results suggest that DvCV1 is a new member of the genus Closterovirus. This is the first report of a closterovirus infecting D. volubilis.

Through high-throughput sequencing, a novel citlodavirus, tentatively named “Myrica rubra citlodavirus 1” (MRV1, accession no. OP374189), was isolated from the leaves of Myrica rubra in Yunnan exhibiting narrow deformity of leaf tips, shrinkage, and chlorosis along the veins. The complete genome sequence was determined and analyzed using cloning and Sanger sequencing. MRV1 is a single-stranded circular non-enveloped DNA virus with a genome size of 3775 nucleotides and contains six open reading frames (ORFs). The virion-sense genome strand encodes a coat protein (CP, nt 750-1,493, 247 aa), two hypothetical movement proteins (V3, nt 382–666, 94 aa; and V2, nt 461–895, 144 aa), and one movement protein (MP, nt 1,527-2,438, 303 aa). The complementary strand of the genome encodes two replication proteins (RepA, nt 3,712-2,834, 292 aa; Rep, nt 2,867-2,553, 104 aa). The MRV1 genome contains the stem-loop motif 5′-TAATATTAC-3′, which is a highly conserved nonanucleotide motif found in the origin of virion-strand replication in geminiviruses. Genome sequence alignment analysis showed that citrus chlorotic dwarf associated virus (CCDaV, accession no. JQ920490) shared the highest nucleotide sequence similarity (66.10% identity) with MRV1. Phylogenetic analysis showed that CCDaV is the closest known relative of MRV1, and that these viruses clustered in a single branch within a clade consisting of citlodaviruses. These results indicate that MRV1 should be regarded as a new species of the genus Citlodavirus in the family Geminiviridae.

Group Ⅲ WRKY transcription factors (TFs) play pivotal roles in responding to the diverse abiotic stress and secondary metabolism of plants. However, the evolution and function of WRKY66 remains unclear. Here, WRKY66 homologs were traced back to the origin of terrestrial plants and found to have been subjected to both motifs’ gain and loss, and purifying selection. A phylogenetic analysis showed that 145 WRKY66 genes could be divided into three main clades (Clade A–C). The substitution rate tests indicated that the WRKY66 lineage was significantly different from others. A sequence analysis displayed that the WRKY66 homologs had conserved WRKY and C2HC motifs with higher proportions of crucial amino acid residues in the average abundance. The AtWRKY66 is a nuclear protein, salt- and ABA- inducible transcription activator. Simultaneously, under salt stress and ABA treatments, the superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT) activities, as well as the seed germination rates of Atwrky66-knockdown plants generated by the clustered, regularly interspaced, short palindromic repeats/CRISPR-associated 9 (CRISPR/Cas9) system, were all lower than those of wild type (WT) plants, but the relative electrolyte leakage (REL) was higher, indicating the increased sensitivities of the knockdown plants to the salt stress and ABA treatments. Moreover, RNA-seq and qRT-PCR analyses revealed that several regulatory genes in the ABA-mediated signaling pathway involved in stress response of the knockdown plants were significantly regulated, being evidenced by the more moderate expressions of the genes. Therefore, the AtWRKY66 likely acts as a positive regulator in the salt stress response, which may be involved in an ABA-mediated signaling pathway.

A novel double-stranded RNA (dsRNA) virus, tentatively named “Paris alphapartitivirus 1” (ParAPV1, OL960006-OL960007), was detected in Paris polyphylla var. yunnanensis plants exhibiting leaf chlorosis and shrinkage symptoms in Yunnan. Its complete genome sequence was determined using Illumina and Sanger sequencing. ParAPV1 has a bipartite genome that consists of dsRNA1 (1,917 bp) encoding the viral RNA-dependent RNA polymerase (RdRp), and dsRNA2 (1,818 bp) encoding the putative coat protein (CP). Sequence comparisons showed that the RdRp and CP of ParAPV1 are most similar to those of pear alphapartitivirus (PpPV2), with 69.97% and 54.21% amino acid sequence identities respectively. Phylogenetic analysis of the RdRp amino acid sequences of ParAPV1 and other partitiviruses showed that ParAPV1 cluster with viruses in a clade containing alphapartitiviruses, and that its closest known relatives are PpPV2 (BBA66577) and rose partitivirus (RoPV, ANQ45203S). Taken together, these results suggest that ParAPV1 should be regarded as a new member of genus Alphapartitivirus in the family Partitiviridae. This is the first report of a partitivirus infecting P. polyphylla var. yunnanensis.