Conserved untranslated regions of multipartite viruses: Natural markers of novel viral genomic components and tags of viral evolution.

Viruses with split genomes are classified as being either segmented or multipartite based on whether their genomic segments occur within a single virion or across different virions. Despite variations in number and sequence during evolution, the genomic segments of many viruses are conserved within the untranslated regions (UTRs). In this study, we present a methodology that combines RNA sequencing with iterative BLASTn of UTRs (https://github.com/qq371260/Iterative-blast-v.1.0) to identify new viral genomic segments. Some novel multipartite-like viruses related to the phylum Kitrinoviricota were annotated using sequencing data from field plant samples and public databases. We identified potentially plant-infecting jingmen-related viruses (order Amarillovirales) and jivi-related viruses (order Martellivirales) with at least six genomic components. The number of RNA molecules associated with a genome of the novel viruses in the families Closteroviridae, Kitaviridae, and Virgaviridae within the order Martellivirales reached five. Several of these viruses seem to represent new taxa at the subgenus, genus, and family levels. The diversity of novel genomic components and the multiple duplication of proteins or protein domains within single or multiple genomic components emphasize the evolutionary roles of genetic recombination (horizontal gene transfer), reassortment, and deletion. The relatively conserved UTRs at the genome level might explain the relationships between monopartite and multipartite viruses, as well as how subviral agents such as defective RNAs and satellite viruses can coexist with their helper viruses.

Effects of RNA silencing during antagonism between citrus exocortis viroid and citrus bark cracking viroid in Etrog citron (Citrus medica).

Citrus exocortis viroid (CEVd) and citrus bark cracking viroid (CBCVd) are two important viroids that infect citrus plants and frequently occur as mixed infections in orchards. However, the mechanism of antagonism between the two viroids in mixed infections remains unclear. The CEVd/CBCVd-citron system and small RNA sequencing (sRNA-seq) were used to study the antagonism. When CBCVd was inoculated before CEVd, the CEVd titre was significantly reduced and the symptoms were attenuated. Viroid-derived sRNAs (vd-sRNAs) from CEVd and CBCVd were predominantly 21-nucleotide (nt) and 22-nt in length and had similar 5' base biases. Homologous sequences of the two viroids in the terminal right (TR) region are rich in vd-sRNAs, and the high frequency vd-sRNAs selected from the CBCVd TR region can be used to degrade the transcripts of CEVd in vivo directly. These results suggest that RNA silencing may play an important role in the antagonism of the two viroids, thus deepening our understanding of the molecular interaction of long noncoding RNAs in woody plants.

Identification and molecular characterization of a novel cytorhabdovirus from rose plants (Rosa chinensis Jacq.).

A novel negative-sense single-stranded RNA virus, tentatively named "rose-associated cytorhabdovirus" (RaCV), was identified by high-throughput sequencing. RaCV is 16,067 nucleotides in length and contains eight open reading frames (ORFs 1-8) encoding a nucleocapsid protein (N), a putative phosphoprotein (P), a putative P3 protein (P3), a putative P4 protein (P4), a putative matrix protein (M), a glycoprotein (G), a putative P7 protein (P7), and an RNA-dependent RNA polymerase (L), respectively. The coding genes are flanked by a 3' leader sequence (228 nt) and a 5' trailer sequence (251 nt) and are separated by conserved intergenic junctions (3'-AUUCUUUUUG(N)nCUN-5'). Phylogenetic analysis showed that RaCV clustered with yerba mate virus A (YmVA) within the cytorhabdovirus clade, and it exhibited low a degree of nt sequence similarity (<40% identity) to other rhabdoviruses. Amino acid sequence comparisons between the putative proteins of RaCV and the corresponding proteins of other cytorhabdoviruses showed that the sequence identity levels were far below the species demarcation cutoff of 80% for cytorhabdoviruses. These results suggest that RaCV should be classified as a new member of the genus Cytorhabdovirus.

Identification and Characterization of Two Novel Geminiviruses Associated with Paper Mulberry (Broussonetia papyrifera) Leaf Curl Disease.

Paper mulberry (Broussonetia papyrifera) is a perennial woody plant used as source material for Cai Lun paper making, in traditional Chinese medicine, and as livestock feed. To identify the presence of viruses in paper mulberry plants affected by a disease with leaf curl symptoms, high-throughput sequencing of total RNA was performed. Analysis of transcriptome libraries allowed the reconstruction of two geminivirus-like genomes. Rolling-circle amplification and PCR with back-to-back primers confirmed the presence of two geminiviruses with monopartite genomes in these plants, with the names paper mulberry leaf curl virus 1 and 2 (PMLCV-1 and PMLCV-2) proposed. The genomes of PMLCV-1 (3,056 nt) and PMLCV-2 (3,757 to 3,763 nt) encode six proteins, with the V4 protein of PMLCV-1 and the V3 proteins of both viruses having low similarities to any known protein in databases. Alternative splicing of an intron, akin to that of mastre-, becurto-, capula-, and grabloviruses, was identified by small RNA (sRNA)-seq and RNA-seq reads mapping to PMLCV-1 and PMLCV-2 antisense transcripts. Phylogenetic analyses and pairwise comparisons showed that PMLCV-1 and PMLCV-2 are most closely related to, but distinct from, two unassigned geminiviruses, citrus chlorotic dwarf associated virus and mulberry mosaic dwarf associated virus, suggesting that they are two new members of the family Geminiviridae. Field investigation confirmed the close association of the two viruses with leaf curl symptoms in paper mulberry plants and that coinfection can aggravate the symptoms.

Molecular characterization of a novel cytorhabdovirus associated with paper mulberry mosaic disease.

A novel cytorhabdovirus, tentatively named "paper mulberry mosaic-associated virus" (PMuMaV), was discovered and identified by transcriptome sequencing (RNA-seq), small RNA sequencing (sRNA-seq), and RT-PCR amplification. The whole-genome sequence of PMuMaV is 13,736 nucleotides (nt) in length and contains six open reading frames (ORFs) encoding a nucleocapsid protein (N), a phosphoprotein (P), a putative movement protein (P3), a matrix protein (M), a glycoprotein (G), and an RNA-dependent RNA polymerase (L). The coding sequences are flanked by a 194-nt leader and a 370-nt trailer sequence at the 3' terminus and 5' terminus, respectively. Pairwise sequence comparisons showed that PMuMaV is related to northern cereal mosaic virus (NCMV, 38.97%), barley yellow striate mosaic virus (BYSMV, 38.86%), and maize yellow striate virus (MYSV, 38.76%), and phylogenetic analysis also placed these viruses together into the same branch, thus suggesting that PMuMaV is a member of a new species in the genus Cytorhabdovirus.

Complete nucleotide sequence of loquat virus A, a member of the family Betaflexiviridae with a novel genome organization.

Analysis of a loquat tree with leaf curl symptoms by deep sequencing revealed a novel virus with a single-stranded RNA genome, for which the name "loquat virus A" (LoVA) was proposed. The complete genome sequence comprised 7553 nucleotides (nt) and an additional poly(A) tail at the 3' terminus. Sequence comparisons of LoVA showed moderate similarity to cherry virus A (CVA), currant virus A (CuVA), and mume virus A (MuVA), which are members of the genus Capillovirus in the family Betaflexiviridae. Phylogenetic analysis of the full-genome nt sequence and replicase-like protein supported the placement of LoVA within the genus Capillovirus. However, it has a distinct genome organization, differing from recognized capilloviruses, as it contains three open reading frames (ORFs), with the coat protein (CP) expressed separately from the replication-associated protein (RP) rather than being encoded in the same ORF. This indicates that LoVA is a novel member of the genus Capillovirus in the family Betaflexiviridae with a distinct genomic organization.

Global Transcriptomic Analysis Reveals Insights into the Response of 'Etrog' Citron (Citrus medica L.) to Citrus Exocortis Viroid Infection.

Citrus exocortis viroid (CEVd) is the causal agent of citrus exocortis disease. We employed CEVd-infected 'Etrog' citron as a system to study the feedback regulation mechanism using transcriptome analysis in this study. Three months after CEVd infection, the transcriptome of fresh leaves was analyzed, and 1530 differentially expressed genes were detected. The replication of CEVd in citron induced upregulation of genes encoding key proteins that were involved in the RNA silencing pathway such as Dicer-like 2, RNA-dependent RNA polymerase 1, argonaute 2, argonaute 7, and silencing defective 3, as well as those genes encoding proteins that are related to basic defense responses. Many genes involved in secondary metabolite biosynthesis and chitinase activity were upregulated, whereas other genes related to cell wall and phytohormone signal transduction were downregulated. Moreover, genes encoding disease resistance proteins, pathogenicity-related proteins, and heat shock cognate 70 kDa proteins were also upregulated in response to CEVd infection. These results suggest that basic defense and RNA silencing mechanisms are activated by CEVd infection, and this information improves our understanding of the pathogenesis of viroids in woody plants.

The mitochondrial genome of booklouse, Liposcelis sculptilis (Psocoptera: Liposcelididae) and the evolutionary timescale of Liposcelis.

Bilateral animals are featured by an extremely compact mitochondrial (mt) genome with 37 genes on a single circular chromosome. To date, the complete mt genome has only been determined for four species of Liposcelis, a genus with economic importance, including L. entomophila, L. decolor, L. bostrychophila, and L. paeta. They belong to A, B, or D group of Liposcelis, respectively. Unlike most bilateral animals, L. bostrychophila, L. entomophila and L. paeta have a bitipartite mt genome with genes on two chromosomes. However, the mt genome of L. decolor has the typical mt chromosome of bilateral animals. Here, we sequenced the mt genome of L. sculptilis, and identified 35 genes, which were on a single chromosome. The mt genome fragmentation is not shared by the D group of Liposcelis and the single chromosome of L. sculptilis differed from those of booklice known in gene content and gene arrangement. We inferred that different evolutionary patterns and rate existed in Liposcelis. Further, we reconstructed the evolutionary history of 21 psocodean taxa with phylogenetic analyses, which suggested that Liposcelididae and Phthiraptera have evolved 134 Ma and the sucking lice diversified in the Late Cretaceous.