Molecular and biological characteristics of a peach latent mosaic viroid PC isolate in peach from China: base mutations in hairpin stems and implications for symptomatology.

Peach latent mosaic viroid (PLMVd) infects peach trees in China and induces a conspicuous albino phenotype (peach calico, PC) that is closely associated with variants containing a 12-14 nucleotide hairpin insertion capped by a U-rich loop. Initially, PC disease distribution was limited to parts of Italy, and it was first detected in the field in China in 2019. To explore the molecular and biological characteristics of PLMVd PC isolates in peach in China, we conducted a comprehensive analysis of disease phenotype development, and investigated the data-associated pathogenicity and in vivo dynamics of Chinese isolate PC-A2 using slash-inoculated into GF-305 peach seedlings. Inoculated seedlings displayed PC symptoms much earlier following topping treatment, and PLMVd infectivity was further assessed using bioassay and semiquantitative RT-PCR experiments. Evolutionary analysis showed that the PC isolate and its progeny variants clustered into a single phylogroup distinct from reference PC-C40 isolates from Italy, and PC-K1 and PC-K2 from South Korea. Some PC-A2 progeny variants from green leaves of PC-expressing seedlings showed unbalanced point mutations in hairpin stems compared with the PC-C40 reference sequence, and constituted a new stem insertion type. The results reveal associations between the recessive phenotypes of peach albino symptoms and base variation in hairpin stem insertions relative to the PC-C40/chloroplastic heat shock protein 90 reference sequence.

Analyses of virus/viroid communities in nectarine trees by next-generation sequencing and insight into viral synergisms implication in host disease symptoms.

We analyzed virus and viroid communities in five individual trees of two nectarine cultivars with different disease phenotypes using next-generation sequencing technology. Different viral communities were found in different cultivars and individual trees. A total of eight viruses and one viroid in five families were identified in a single tree. To our knowledge, this is the first report showing that the most-frequently identified viral and viroid species co-infect a single individual peach tree, and is also the first report of peach virus D infecting Prunus in China. Combining analyses of genetic variation and sRNA data for co-infecting viruses/viroid in individual trees revealed for the first time that viral synergisms involving a few virus genera in the Betaflexiviridae, Closteroviridae, and Luteoviridae families play a role in determining disease symptoms. Evolutionary analysis of one of the most dominant peach pathogens, peach latent mosaic viroid (PLMVd), shows that the PLMVd sequences recovered from symptomatic and asymptomatic nectarine leaves did not all cluster together, and intra-isolate divergent sequence variants co-infected individual trees. Our study provides insight into the role that mixed viral/viroid communities infecting nectarine play in host symptom development, and will be important in further studies of epidemiological features of host-pathogen interactions.

Complete Genome Sequence of a Divergent Isolate of Cherry Virus A from Prunus avium in China.

Here, we report the complete genome sequence of a divergent cherry virus A (CVA) isolate (ChYT56) from Prunus avium in China. The genome nucleotide sequence has low identity (80.7%) with a CVA from P. avium (GenBank accession number FN691959) and high identity (97%) with a CVA from P. armeniaca (GenBank accession number LC125634).

Complete nucleotide sequence of a new virus, peach chlorotic leaf spot virus, isolated from flat peach in China.

Contigs with sequence homologies to apple chlorotic leaf spot virus (ACLSV) were identified by high-throughput sequencing analysis in three peach samples. Complete genomic sequences of RP19-1 and RP19-2 of the virus consisted of 7,466 and 7,465 nucleotides (nts), respectively, excluding the poly (A) tails. They shared the highest identity with Ta Tao 5, but lower than 70% of sequence similarity with other ACLSV isolates. Furthermore, phylogenetic analysis revealed that these two isolates clustered with Ta Tao 5, which is distinct from other ACLSV isolates. According to the criteria for species demarcation within the genus Trichovirus, these two isolates as well as Ta Tao 5 should be recognized as a new virus species, tentatively named "Peach chlorotic leaf spot virus".

Further insight into genetic variation and haplotype diversity of Cherry virus A from China.

Cherry virus A (CVA) infection appears to be prevalent in cherry plantations worldwide. In this study, the diversity of CVA isolates from 31 cherry samples collected from different orchards around Bohai Bay in northeastern China was analyzed. The complete genome of one of these isolates, ChYT52, was found to be 7,434 nt in length excluding the poly (A) tail. It shares between 79.9-98.7% identity with CVA genome sequences in GenBank, while its RdRp core is more divergent (79.1-90.7% nt identity), likely as a consequence of a recombination event. Phylogenetic analysis of ChYT52 genome with CVA genomes in Genbank resulted in at least 7 major clusters plus additional 5 isolates alone at the end of long branches suggesting the existence of further phylogroups diversity in CVA. The genetic diversity of Chinese CVA isolates from 31 samples and GenBank sequences were analyzed in three genomic regions that correspond to the coat protein, the RNA-dependent RNA polymerase core region, and the movement protein genes. With few exceptions likely representing further recombination impact, the trees various trees are largely congruent, indicating that each region provides valuable phylogenetic information. In all cases, the majority of the Chinese CVA isolates clustering in phylogroup I, together with the X82547 reference sequence from Germany. Statistically significant negative values were obtained for Tajima's D in the three genes for phylogroup I, suggesting that it may be undergoing a period of expansion. There was considerable haplotype diversity in the individual samples and more than half samples contained genetically diverse haplotypes belonging to different phylogroups. In addition, a number of statistically significant recombination events were detected in CVA genomes or in the partial genomic sequences indicating an important contribution of recombination to CVA evolution. This work provides a foundation for elucidation of the epidemiological characteristics and evolutionary history of CVA populations.

Global Transcriptomic Changes Induced by Infection of Cucumber (Cucumis sativus L.) with Mild and Severe Variants of Hop Stunt Viroid.

Fifteen years after transfer to hops, hop stunt viroid-grapevine (HSVd-g) was replaced by HSVd-hop (HSVd-h), a sequence variant that contains changes at five different positions. HSVd-g54 is a laboratory mutant derived from HSVd-g that differs from its progenitor by a single G to A substitution at position 54. While infection by HSVd-h induces only mild stunting in cucumber (Cucumis sativus L.), HSVd-g54 induces much more severe symptoms in this indicator host. Comparison of transcriptome profiles of cucumber infected with HSVd-h or HSVd-g54 with those of mock-inoculated controls obtained by whole transcriptome shotgun sequencing revealed that many genes related to photosynthesis were down-regulated following infection. In contrast, genes encoding RNA-dependent RNA polymerase 1 (CsRDR1), especially CsRDR1c1 and CsRDR1c2, as well as those related to basal defense responses were up-regulated. Expression of genes associated with phytohormone signaling pathways were also altered, indicating that viroid infection initiates a complex array of changes in the host transcriptome. HSVd-g54 induced an earlier and stronger response than HSVd-h, and further examination of these differences will contribute to a better understanding of the mechanisms that determine viroid pathogenicity.

Identification and characterization of a novel geminivirus with a monopartite genome infecting apple trees.

A novel circular DNA virus sequence has been identified through next-generation sequencing and in silico assembly of small RNAs of 21-24 nt from an apple tree grown in China. The virus genome was cloned using two independent approaches and sequenced. With a size of 2932 nt, it showed the same genomic structure and conserved origin of replication reported for members of the family Geminiviridae. However, the low nucleotide and amino acid sequence identity with known geminiviruses indicated that it was a novel virus, for which the provisional name apple geminivirus (AGV) is proposed. Rolling circle amplification followed by RFLP analyses indicated that AGV was a virus with a monopartite DNA genome. This result was in line with bioassays showing that the cloned viral genome was infectious in several herbaceous plants (Nicotiana bethamiana, Nicotiana glutinosa and Solanum lycopersicum), thus confirming it was complete and biologically active, although no symptoms were observed in these experimental hosts. AGV genome structure and phylogenetic analyses did not support the inclusion of this novel species in any of the established genera in the family Geminiviridae. A survey of 165 apple trees grown in four Chinese provinces showed a prevalence of 7.2% for AGV, confirming its presence in several cultivars and geographical areas in China, although no obvious relationship between virus infection and specific symptoms was found.

Identification and molecular characterization of a novel monopartite geminivirus associated with mulberry mosaic dwarf disease.

High-throughput sequencing of small RNAs allowed the identification of a novel DNA virus in a Chinese mulberry tree affected by a disease showing mosaic and dwarfing symptoms. Rolling-circle amplification and PCR with specific primers, followed by sequencing of eleven independent full-length clones, showed that this virus has a monopartite circular DNA genome (∼ 2.95 kb) containing ORFs in both polarity strands, as reported previously for geminiviruses. A field survey showed the close association of the virus with diseased mulberries, so we tentatively named the virus mulberry mosaic dwarf-associated virus (MMDaV). The MMDaV genome codes for five and two putative proteins in the virion-sense and in the complementary-sense strands, respectively. Although three MMDaV virion-sense putative proteins did not share sequence homology with any protein in the databases, functional domains [coiled-coil and transmembrane (TM) domains] were identified in two of them. In addition, the protein containing a TM domain was encoded by an ORF located in a similar genomic position in MMDaV and in several other geminiviruses. As reported for members of the genera Mastrevirus and Becurtovirus, MMDaV replication-associated proteins are expressed through the alternative splicing of an intron, which was shown to be functional in vivo. A similar intron was found in the genome of citrus chlorotic dwarf-associated virus (CCDaV), a divergent geminivirus found recently in citrus. On the basis of pairwise comparisons and phylogenetic analyses, CCDaV and MMDaV appear to be closely related to each other, thus supporting their inclusion in a putative novel genus in the family Geminiviridae.

Discovery of replicating circular RNAs by RNA-seq and computational algorithms.

Replicating circular RNAs are independent plant pathogens known as viroids, or act to modulate the pathogenesis of plant and animal viruses as their satellite RNAs. The rate of discovery of these subviral pathogens was low over the past 40 years because the classical approaches are technical demanding and time-consuming. We previously described an approach for homology-independent discovery of replicating circular RNAs by analysing the total small RNA populations from samples of diseased tissues with a computational program known as progressive filtering of overlapping small RNAs (PFOR). However, PFOR written in PERL language is extremely slow and is unable to discover those subviral pathogens that do not trigger in vivo accumulation of extensively overlapping small RNAs. Moreover, PFOR is yet to identify a new viroid capable of initiating independent infection. Here we report the development of PFOR2 that adopted parallel programming in the C++ language and was 3 to 8 times faster than PFOR. A new computational program was further developed and incorporated into PFOR2 to allow the identification of circular RNAs by deep sequencing of long RNAs instead of small RNAs. PFOR2 analysis of the small RNA libraries from grapevine and apple plants led to the discovery of Grapevine latent viroid (GLVd) and Apple hammerhead viroid-like RNA (AHVd-like RNA), respectively. GLVd was proposed as a new species in the genus Apscaviroid, because it contained the typical structural elements found in this group of viroids and initiated independent infection in grapevine seedlings. AHVd-like RNA encoded a biologically active hammerhead ribozyme in both polarities, and was not specifically associated with any of the viruses found in apple plants. We propose that these computational algorithms have the potential to discover novel circular RNAs in plants, invertebrates and vertebrates regardless of whether they replicate and/or induce the in vivo accumulation of small RNAs.

Simultaneous detection and identification of four viruses infecting pepino by multiplex RT-PCR.

Potato virus M (PVM), pepino mosaic virus (PepMV), tomato mosaic virus (ToMV), and potato virus S (PVS) infect pepino and cause serious crop losses. In this study, a multiplex RT-PCR method was developed for simultaneous detection and differentiation of PVS, ToMV, PepMV and PVM. The method was highly reliable and sensitive; validation was accomplished by testing pepino samples collected from different regions of China. In this survey, PVM, ToMV and PVS were detected in 37.0 %, 31.0 % and 5.5 % of samples tested, respectively, confirming the widespread occurrence of these three viruses in China. PepMV was not detected in any of the samples, which indicated that this virus may not be prevalent in China. The results suggest that the new multiplex RT-PCR method has potential to be used routinely for surveys of pepino for virus infection.

Complete nucleotide sequences of two isolates of cherry green ring mottle virus from peach (Prunus persica) in China.

Two complete nucleotide sequences of cherry green ring mottle virus (CGRMV) isolated from peach in Hebei (Hs10) and Fujian (F9) Provinces, China, were determined. Five open reading frames (ORFs) were found in the genomes of both isolates. The F9 and Hs10 isolates shared 82.2 % and 83.4-94.4 % nucleotide sequence identity, respectively, with two CGRMV isolates from cherry. Analysis of the nucleotide and amino acid sequences from the five ORFs of both isolates showed that Hs10 shares the greatest sequence identity with P1A (GenBank AJ291761) from cherry. Phylogenetic analysis indicated that CGRMV isolates from peach and cherry are closely related to members of the genus Foveavirus.

Selection and heritability of resistance to Bacillus thuringiensis subsp kurstaki and transgenic cotton in Helicoverpa armigera (Lepidoptera: Noctuidae).

Compared with an unselected susceptible population, a cotton bollworm, Helicoverpa armigera (Hübner), population selected for 22 generations with transgenic cotton leaves (modified Cry1A) in the laboratory developed 11.0-fold resistance to Cry1Ac (one single-protein product MVPII). Resistance to Bacillus thuringiensis Berliner subsp kurstaki (Btk) was selected for 22 generations with a 5.2-fold increase in LC50. The estimated realized heritabilities (h2) of resistance for transgenic-cotton- and Btk-selected populations were 0.1008 and 0.2341, respectively. This reflects the higher phenotypic variation in response to Cry1Ac in the transgenic-cotton-selected population. This variation may have been caused by differences in protein toxin levels expressed in different growth stages of the transgenic cotton. Because of the different slopes of the probit regression lines between Cry1Ac and Btk, the estimated realized h2 cannot be used visually to compare resistance development to Cry1Ac and Btk in H armigera. Thus, the response quotient (Q) of resistance was also estimated. The Q values of resistance for transgenic-cotton- and Btk-selected populations were 0.0763 and 0.0836, respectively. This showed that the rate of resistance development would be similar in both selection populations. This result indicates that the selection of resistance using transgenic cotton is different from that selected using the single toxin. Resistance risk to transgenic cotton and Btk in field populations was assessed assuming different pressures of selection by using the estimated h2. Assuming the h2 of resistance in a field population was half of the estimated h2, and the population received prolonged and uniform exposure to transgenic cotton or Btk causing >70% mortality in each generation, we predicted that resistance would increase 10-fold after <23 generations for Cry1Ac in transgenic cotton-selected-populations and after <21 generations for Btk in Btk-selected populations. Cross-resistance would be expected after <48 generations for Btk in transgenic-cotton-selected populations and after <21 generations for Cry1Ac in Btk-selected population. The results show that the potential to evolve resistance is similar in both transgenic-cotton- and Btk-selected populations, but that cross-resistance development to Btk is slower in transgenic-cotton-selected populations than cross-resistance development to Cry1Ac in Btk-selected populations.