Comparative Transcriptomic Analysis of Genes in the 20-Hydroxyecdysone Biosynthesis in the Fern Microsorum scolopendria towards Challenges with Foliar Application of Chitosan.

Microsorum scolopendria is an important medicinal plant that belongs to the Polypodiaceae family. In this study, we analyzed the effects of foliar spraying of chitosan on growth promotion and 20-hydroxyecdysone (20E) production in M. scolopendria. Treatment with chitosan at a concentration of 50 mg/L in both young and mature sterile fronds induced the highest increase in the amount of accumulated 20E. Using RNA sequencing, we identified 3552 differentially expressed genes (DEGs) in response to chitosan treatment. The identified DEGs were associated with 236 metabolic pathways. We identified several DEGs involved in the terpenoid and steroid biosynthetic pathways that might be associated with secondary metabolite 20E biosynthesis. Eight upregulated genes involved in cholesterol and phytosterol biosynthetic pathway, five upregulated genes related to the methylerythritol 4-phosphate (MEP) and mevalonate (MVA) pathways, and several DEGs that are members of cytochrome P450s and ABC transporters were identified. Quantitative real-time RT-PCR confirmed the results of RNA-sequencing. Taken together, we showed that chitosan treatment increased plant dry weight and 20E accumulation in M. scolopendria. RNA-sequencing and DEG analyses revealed key enzymes that might be related to the production of the secondary metabolite 20E in M. scolopendria.

Transcriptional Analysis of the Differences between ToLCNDV-India and ToLCNDV-ES Leading to Contrary Symptom Development in Cucumber.

Tomato leaf curl New Delhi virus-ES (ToLCNDV-ES), a high threat to cucurbits in the Mediterranean Basin, is listed as a different strain from the Asian ToLCNDV isolates. In this study, the infectivity of two clones previously isolated from Italy and Pakistan were compared in cucumbers, which resulted in the opposite symptom appearance. The swapping subgenome was processed; however, the mechanisms related to the disease phenotype remain unclear. To identify the disease-associated genes that could contribute to symptom development under the two ToLCNDV infections, the transcriptomes of ToLCNDV-infected and mock-inoculated cucumber plants were compared 21 days postinoculation. The number of differentially expressed genes in ToLCNDV-India-infected plants was 10 times higher than in ToLCNDV-ES-infected samples. The gene ontology (GO) and pathway enrichment were analyzed using the Cucurbits Genomics Database. The flavonoid pathway-related genes were upregulated in ToLCNDV-ES, but some were downregulated in ToLCNDV-India infection, suggesting their role in resistance to the two ToLCNDV infections. The relative expression levels of the selected candidate genes were validated by qRT-PCR under two ToLCNDV-infected conditions. Our results reveal the different infectivity of the two ToLCNDVs in cucumber and also provide primary information based on RNA-seq for further analysis related to different ToLCNDV infections.

Comparative Analysis of Water Extracts from Roselle (Hibiscus sabdariffa L.) Plants and Callus Cells: Constituents, Effects on Human Skin Cells, and Transcriptome Profiles.

Roselle (Hibiscus sabdariffa L.) is a plant that has traditionally been used in various food and beverage products. Here, we investigated the potential of water extracts derived from Roselle leaves and callus cells for cosmetic and pharmaceutical purposes. We generated calluses from Roselle leaves and produced two different water extracts through heat extraction, which we named Hibiscus sabdariffa plant extract (HSPE) and Hibiscus sabdariffa callus extract (HSCE). HPLC analysis showed that the two extracts have different components, with nucleic acids and metabolites such as phenylalanine and tryptophan being the most common components in both extracts. In vitro assays demonstrated that HSCE has strong anti-melanogenic effects and functions for skin barrier and antioxidant activity. Transcriptome profiling of human skin cells treated with HSPE and HSCE showed significant differences, with HSPE having more effects on human skin cells. Up-regulated genes by HSPE function in angiogenesis, the oxidation-reduction process, and glycolysis, while up-regulated genes by HSCE encode ribosome proteins and IFI6, functioning in the healing of radiation-injured skin cells. Therefore, we suggest that the two extracts from Roselle should be applied differently for cosmetics and pharmaceutical purposes. Our findings demonstrate the potential of Roselle extracts as a natural source for skincare products.

Virus classification based on in-depth sequence analyses and development of demarcation criteria using the Betaflexiviridae as a case study.

Currently, many viruses are classified based on their genome organization and nucleotide/amino acid sequence identities of their capsid and replication-associated proteins. Although biological traits such as vector specificities and host range are also considered, this later information is scarce for the majority of recently identified viruses, characterized only from genomic sequences. Accordingly, genomic sequences and derived information are being frequently used as the major, if not only, criteria for virus classification and this calls for a full review of the process. Herein, we critically addressed current issues concerning classification of viruses in the family Betaflexiviridae in the era of high-throughput sequencing and propose an updated set of demarcation criteria based on a process involving pairwise identity analyses and phylogenetics. The proposed framework has been designed to solve the majority of current conundrums in taxonomy and to facilitate future virus classification. Finally, the analyses performed herein, alongside the proposed approaches, could be used as a blueprint for virus classification at-large.

De Novo Assembly and Annotation of the Vaginal Metatranscriptome Associated with Bacterial Vaginosis.

The vaginal microbiome plays an important role in women's health and disease. Here we reanalyzed 40 vaginal transcriptomes from a previous study of de novo assembly (metaT-Assembly) followed by functional annotation. We identified 286,293 contigs and further assigned them to 25 phyla, 209 genera, and 339 species. Lactobacillus iners and Lactobacillus crispatus dominated the microbiome of non-bacterial vaginosis (BV) samples, while a complex of microbiota was identified from BV-associated samples. The metaT-Assembly identified a higher number of bacterial species than the 16S rRNA amplicon and metaT-Kraken methods. However, metaT-Assembly and metaT-Kraken exhibited similar major bacterial composition at the species level. Binning of metatranscriptome data resulted in 176 bins from major known bacteria and several unidentified bacteria in the vagina. Functional analyses based on Clusters of Orthologous Genes (COGs) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways suggested that a higher number of transcripts were expressed by the microbiome complex in the BV-associated samples than in non-BV-associated samples. The KEGG pathway analysis with an individual bacterial genome identified specific functions of the identified bacterial genome. Taken together, we demonstrated that the metaT-Assembly approach is an efficient tool to understand the dynamic microbial communities and their functional roles associated with the human vagina.

Viromes of 15 Pepper (Capsicum annuum L.) Cultivars.

Pepper (Capsicum annuum L.) plants produce berry fruits that are used as spices. Here, we examined the viromes of 15 pepper cultivars through RNA sequencing. We obtained 1,325 virus-associated contigs derived from 8 virus species. Bean broad wilt virus 2 (BBWV2) and cucumber mosaic virus (CMV) were identified as the major viruses infecting pepper plants, followed by potato virus Y, bell pepper endornavirus, and hot pepper endornavirus. The proportion of viral reads in each transcriptome ranged from 0.04% to 24.5%. BBWV2 was the dominant virus in seven cultivars, whereas CMV was dominant in five cultivars. All the bell pepper cultivars showed severe viral disease symptoms, whereas the commercially developed hot pepper cultivars were asymptomatic or had mild symptoms. In addition, 111 complete viral segments were obtained from 7 viruses. Based on the obtained viral genomes, the genetic relationship between the identified viruses and quasispecies of BBWV2 and CMV in each pepper plant was determined. Newly designed primers for nine viruses confirmed the results of RNA sequencing. Taken together, this study, for the first time, provides a comprehensive overview of viromes in 15 major pepper cultivars through RNA sequencing.

Unveiling Mycoviromes Using Fungal Transcriptomes.

Viruses infecting fungi are referred to as mycoviruses. Here, we carried out in silico mycovirome studies using public fungal transcriptomes mostly derived from mRNA libraries. We identified 468 virus-associated contigs assigned to 5 orders, 21 families, 26 genera, and 88 species. We assembled 120 viral genomes with diverse RNA and DNA genomes. The phylogenetic tree and genome organization unveiled the possible host origin of mycovirus species and diversity of their genome structures. Most identified mycoviruses originated from fungi; however, some mycoviruses had strong phylogenetic relationships with those from insects and plants. The viral abundance and mutation frequency of mycoviruses were very low; however, the compositions and populations of mycoviruses were very complex. Although coinfection of diverse mycoviruses in the fungi was common in our study, most mycoviromes had a dominant virus species. The compositions and populations of mycoviruses were more complex than we expected. Viromes of Monilinia species revealed that there were strong deviations in the composition of viruses and viral abundance among samples. Viromes of Gigaspora species showed that the chemical strigolactone might promote virus replication and mutations, while symbiosis with endobacteria might suppress virus replication and mutations. This study revealed the diversity and host distribution of mycoviruses.

Evolution and Phylogeny of Soybean Mosaic Virus Based on 143 Complete Genomes.

Soybean mosaic virus (SMV) of the genus Potyvirus is an important virus in cultivated soybeans. Here, we obtained 7 SMV genomes from soybean germplasms using RNA sequencing and conducted a comprehensive evolutionary and phylogenetic study of 143 SMV genomes derived from 10 plant species and 12 countries. The phylogenetic tree we constructed using coding DNA sequences revealed the existence of nine clades of SMV isolates/strains. Recombination analysis revealed 76 recombinant events and 141 recombinants in total. Clades 1 and 3 contain the most common SMV pathotypes, including G1 through G7, which are distributed worldwide. Clade 2 includes several Chinese SMV pathotypes. The SMV isolates were further divided into two groups. The SMV isolates in the first group, including clades 8 and 9, were identified from Pinellia and Atractylodes species, whereas those in the second group (clades 1 through 7) were mostly found in cultivated soybeans. The SMV polyprotein undergoes positive selection, whereas most mature proteins, except for the P1 protein, undergo negative selection. The P1 protein of SMV isolates in group 1 may be highly correlated with host adaptation. This study provides strong evidence that recombination and plant hosts are powerful forces driving the genetic diversity of the SMV genome.

De Novo Transcriptome Assembly of Two Microsorum Fern Species Identifies Enzymes Required for Two Upstream Pathways of Phytoecdysteroids.

Microsorum species produce a high amount of phytoecdysteroids (PEs), which are widely used in traditional medicine in the Pacific islands. The PEs in two different Microsorum species, M. punctatum (MP) and M. scolopendria (MS), were examined using high-performance liquid chromatography (HPLC). In particular, MS produces a high amount of 20-hydroxyecdysone, which is the main active compound in PEs. To identify genes for PE biosynthesis, we generated reference transcriptomes from sterile frond tissues using the NovaSeq 6000 system. De novo transcriptome assembly after deleting contaminants resulted in 57,252 and 54,618 clean transcripts for MP and MS, respectively. The clean Microsorum transcripts for each species were annotated according to gene ontology terms, UniProt pathways, and the clusters of the orthologous group protein database using the MEGAN6 and Sma3s programs. In total, 1852 and 1980 transcription factors were identified for MP and MS, respectively. We obtained transcripts encoding for 38 and 32 enzymes for MP and MS, respectively, potentially involved in mevalonate and sterol biosynthetic pathways, which produce precursors for PE biosynthesis. Phylogenetic analyses revealed many redundant and unique enzymes between the two species. Overall, this study provides two Microsorum reference transcriptomes that might be useful for further studies regarding PE biosynthesis in Microsorum species.

Using RNA-Sequencing Data to Examine Tissue-Specific Garlic Microbiomes.

Garlic (Allium sativum) is a perennial bulbous plant. Due to its clonal propagation, various diseases threaten the yield and quality of garlic. In this study, we conducted in silico analysis to identify microorganisms, bacteria, fungi, and viruses in six different tissues using garlic RNA-sequencing data. The number of identified microbial species was the highest in inflorescences, followed by flowers and bulb cloves. With the Kraken2 tool, 57% of identified microbial reads were assigned to bacteria and 41% were assigned to viruses. Fungi only made up 1% of microbial reads. At the species level, Streptomyces lividans was the most dominant bacteria while Fusarium pseudograminearum was the most abundant fungi. Several allexiviruses were identified. Of them, the most abundant virus was garlic virus C followed by shallot virus X. We obtained a total of 14 viral genome sequences for four allexiviruses. As we expected, the microbial community varied depending on the tissue types, although there was a dominant microorganism in each tissue. In addition, we found that Kraken2 was a very powerful and efficient tool for the bacteria using RNA-sequencing data with some limitations for virome study.

Comparative Study of Metagenomics and Metatranscriptomics to Reveal Microbiomes in Overwintering Pepper Fruits.

Red pepper (Capsicum annuum, L.), is one of the most important spice plants in Korea. Overwintering pepper fruits are a reservoir of various microbial pepper diseases. Here, we conducted metagenomics (DNA sequencing) and metatranscriptomics (RNA sequencing) using samples collected from three different fields. We compared two different library types and three different analytical methods for the identification of microbiomes in overwintering pepper fruits. Our results demonstrated that DNA sequencing might be useful for the identification of bacteria and DNA viruses such as bacteriophages, while mRNA sequencing might be beneficial for the identification of fungi and RNA viruses. Among three analytical methods, KRAKEN2 with raw data reads (KRAKEN2_R) might be superior for the identification of microbial species to other analytical methods. However, some microbial species with a low number of reads were wrongly assigned at the species level by KRAKEN2_R. Moreover, we found that the databases for bacteria and viruses were better established as compared to the fungal database with limited genome data. In summary, we carefully suggest that different library types and analytical methods with proper databases should be applied for the purpose of microbiome study.

First report of Cherry virus F infecting Japanese plum in Korea.

Cherry virus F (CVF) is a tentative member of the genus Fabavirus in the family Secoviridae, consisting of two RNA segments (Koloniuk et al. 2018). To date, CVF has been documented in only sweet cherry (Prunus avium) in the Czech Republic (Koloniuk et al. 2018), Canada, and Greece. In May 2014, we collected leaf samples from four symptomatic (leaf spots and dapple fruits) and two asymptomatic Japanese plum cultivars (Sun and Gadam) grown in an orchard in Hoengseong, South Korea, to identify viruses and viroids infecting plum trees. Total RNA from individual plum trees was extracted using two commercial kits: Fruit-mate for RNA Purification Kit (Takara, Shiga, Japan) and RNeasy Plant Mini Kit (Qiagen, Hilden, Germany). We generated six mRNA libraries from the six different plum cultivars for RNA-sequencing using the TruSeq RNA Library Preparation Kit v2 (Illumina, CA, U.S.A.) as described previously (Jo et al. 2017). The mRNA libraries were paired-end (2 X 100 bp) sequenced with a HiSeq 2000 system (Macrogen, Seoul, Korea). The raw sequence reads were de novo assembled by Trinity program v. 2.8.6, with default parameters (Haas et al. 2013). The assembled contigs were subjected to BLASTX search against the non-redundant protein database in NCBI. Of the two asymptomatic cultivars, the transcriptome of asymptomatic plum cv. Gadam contained five contigs specific to CVF. Two and three contigs were specific to CVF RNA1 (2,571 reads, coverage 42.15%) and RNA2 (2,025 reads, coverage 53.04%), respectively. The size of these five contigs ranged from 241 to 5,986 bp. Contigs of 5,986 and 3,867 bp in length, referred to as CVF isolate Gadam RNA1 (GenBank MN896996) and RNA2 (GenBank MN896995), respectively, were subjected to BLASTP search against NCBI's non-redundant protein database. The results showed that the polyprotein sequences of RNA1 and RNA2 shared 95.3% and 93.11% amino acid identities with isolates SwC-H_1a from the Czech Republic (GenBank acc. no. AWB36326) and Stac-3B_c8 from Canada (AZZ10055), respectively. To confirm the infection of CVF in cv. Gadam, RT-PCR was conducted using CVF RNA1-specific primers designed based on the CVF reference genome sequences (MH998210 and MH998216), including 5'-CCACCAAATAGGCAAGAGGTCAC-3' (position 3190-3212) and 5'-CACAATCACCATCAATGGTCTCTGC-3' (position 3742-3766), and CVF RNA2-specific primers, including 5'-CTGCTTTATGATGCTAGACATCAAGATG-3' (position 1015-1042) and 5'-ACAATAGGCATGCTCATCTCAACCTC-3' (position 1594-1619). We amplified 577-bp RNA1-specific and 605-bp RNA2-specific amplicons that were cloned and then performed Sanger sequencing. Sequencing of the cloned amplicons for isolate Gadam RNA1 (GenBank MN896993) and RNA2 (GenBank MN896994) revealed values of 99.48% and 99.17% nucleotide identity to that of RNA1 and RNA2 determined by high-throughput sequencing, respectively. Additionally, we tested five plants for each of the six plum cultivars grown in the same orchard. The detection of CVF was carried out through PCR using the primers and protocol described above. Of the 30 trees, CVF was detected in three trees of cv. Gadam by both primer pairs. To our knowledge, this is the first report of CVF infecting Japanese plum and the first report of the virus in Korea. However, its prevalence in other Prunus species, including apricot, European plum, and peach, should be further elucidated.

Identification of Viruses and Viroids Infecting Tomato and Pepper Plants in Vietnam by Metatranscriptomics.

Tomato (Lycopersicum esculentum L.) and pepper (Capsicum annuum L.) plants belonging to the family Solanaceae are cultivated worldwide. The rapid development of next-generation sequencing (NGS) technology facilitates the identification of viruses and viroids infecting plants. In this study, we carried out metatranscriptomics using RNA sequencing followed by bioinformatics analyses to identify viruses and viroids infecting tomato and pepper plants in Vietnam. We prepared a total of 16 libraries, including eight tomato and eight pepper libraries derived from different geographical regions in Vietnam. We identified a total of 602 virus-associated contigs, which were assigned to 18 different virus species belonging to nine different viral genera. We identified 13 different viruses and two viroids infecting tomato plants and 12 viruses and two viroids infecting pepper plants with viruses as dominantly observed pathogens. Our results showed that multiple infection of different viral pathogens was common in both plants. Moreover, geographical region and host plant were two major factors to determine viral populations. Taken together, our results provide the comprehensive overview of viral pathogens infecting two important plants in the family Solanaceae grown in Vietnam.

Differential Contribution of RNA Interference Components in Response to Distinct Fusarium graminearum Virus Infections.

The mechanisms of RNA interference (RNAi) as a defense response against viruses remain unclear in many plant-pathogenic fungi. In this study, we used reverse genetics and virus-derived small RNA profiling to investigate the contributions of RNAi components to the antiviral response against Fusarium graminearum viruses 1 to 3 (FgV1, -2, and -3). Real-time reverse transcription-quantitative PCR (qRT-PCR) indicated that infection of Fusarium graminearum by FgV1, -2, or -3 differentially induces the gene expression of RNAi components in F. graminearum Transcripts of the DICER-2 and AGO-1 genes of F. graminearum (FgDICER-2 and FgAGO-1) accumulated at lower levels following FgV1 infection than following FgV2 or FgV3 infection. We constructed gene disruption and overexpression mutants for each of the Argonaute and dicer genes and for two RNA-dependent RNA polymerase (RdRP) genes and generated virus-infected strains of each mutant. Interestingly, mycelial growth was significantly faster for the FgV1-infected FgAGO-1 overexpression mutant than for the FgV1-infected wild type, while neither FgV2 nor FgV3 infection altered the colony morphology of the gene deletion and overexpression mutants. FgV1 RNA accumulation was significantly decreased in the FgAGO-1 overexpression mutant. Furthermore, the levels of induction of FgAGO-1, FgDICER-2, and some of the FgRdRP genes caused by FgV2 and FgV3 infection were similar to those caused by hairpin RNA-induced gene silencing. Using small RNA sequencing analysis, we documented different patterns of virus-derived small interfering RNA (vsiRNA) production in strains infected with FgV1, -2, and -3. Our results suggest that the Argonaute protein encoded by FgAGO-1 is required for RNAi in F. graminearum, that FgAGO-1 induction differs in response to FgV1, -2, and -3, and that FgAGO-1 might contribute to the accumulation of vsiRNAs in FgV1-infected F. graminearumIMPORTANCE To increase our understanding of how RNAi components in Fusarium graminearum react to mycovirus infections, we characterized the role(s) of RNAi components involved in the antiviral defense response against Fusarium graminearum viruses (FgVs). We observed differences in the levels of induction of RNA silencing-related genes, including FgDICER-2 and FgAGO-1, in response to infection by three different FgVs. FgAGO-1 can efficiently induce a robust RNAi response against FgV1 infection, but FgDICER genes might be relatively redundant to FgAGO-1 with respect to antiviral defense. However, the contribution of this gene in the response to the other FgV infections might be small. Compared to previous studies of Cryphonectria parasitica, which showed dicer-like protein 2 and Argonaute-like protein 2 to be important in antiviral RNA silencing, our results showed that F. graminearum developed a more complex and robust RNA silencing system against mycoviruses and that FgDICER-1 and FgDICER-2 and FgAGO-1 and FgAGO-2 had redundant roles in antiviral RNA silencing.

Transcriptome profiles of tomato plants after neutron irradiation and infection with TYLCV.

Ionizing radiation is ubiquitous in the environment and can cause mutagenesis in living organisms. In this study, we examined the effects of neutron irradiation on tomato plants. Neutron irradiation decreased tomato germination rates, but most irradiated tomato plants did not show any significant phenotype. However, tomato mutants infected by Tomato yellow leaf curl virus (TYLCV) displayed resistance against TYLCV compared to the wild type (WT), which showed disease symptoms. RNA-Seq data demonstrated that the expression profiles of eight tomato mutants were significantly different from that of the WT. The transcriptomes obtained from presoaked seeds were highly altered compared to those of dry seeds. Increased irradiation time resulted in severe changes in the tomato transcriptome; however, different neutron irradiation intensities affected the expressions of different sets of genes. A high number of single-nucleotide polymorphisms in tomato transcriptomes suggest that neutron irradiation strongly impacts plant transcriptomes. The transition/transversion values among mutants were almost constant and were lower than that of the non-irradiated sample (WT), suggesting that neutron irradiation caused an effect. Taken together, this is the first report showing the effects of neutron irradiation on tomato plants by transcriptome analyses.

RNA viromes of the oriental hybrid lily cultivar "Sorbonne".

BACKGROUND: The lily is a perennial flowering plant belonging to the genus Lilium in the family Liliaceae. Most cultivated lily plants are propagated by bulbs. Therefore, numerous lily bulbs are frequently infected by diverse viruses causing viral diseases. To date, no study has examined the viromes of plants of one type with identical genetic backgrounds collected from different geographical regions. RESULTS: Here, we examined different viromes of the lily cultivar "Sorbonne" using 172 gigabytes of transcriptome data composed of 23 libraries from four different projects for the cultivar "Sorbonne." We identified 396 virus-associated contigs from all but one library. We identified six different viruses, including Plantago asiatica mosaic virus (PlAMV), Cucumber mosaic virus (CMV), Lily symptomless virus (LSV), Tulip virus X (TVX), Lily mottle virus (LMoV), and Tobacco rattle virus (TRV). Of them, PlAMV was the most common virus infecting the lily. Scale and flower samples possessed a high number of virus-associated reads. We assembled 32 nearly complete genomes for the six identified viruses possessing the polyadenylate tails. Genomes of all six viruses were highly conserved in the lily cultivar "Sorbonne" based on mutation analysis. We identified defective RNAs from LSV, TVX, and PlAMV localized in the triple gene block region. Phylogenetic analyses showed that virus genomes are highly correlated with geographical regions and host plants. CONCLUSIONS: We conducted comprehensive virome analyses of a single lily cultivar, "Sorbonne," using transcriptome data. Our results shed light on an array of lily virome-associated topics, including virus identification, the dominant virus, virus accumulation in different plant tissues, virus genome assembly, virus mutation, identification of defective RNAs, and phylogenetic relationships of identified viruses. Taken together, we provide very useful methods and valuable results that can be applied in other virome-associated studies.

Identification of viral domains integrated into Arabidopsis proteome.

Horizontal gene transfer (HGT) contributes to the genome evolution of living organisms. In particular, several recent studies provide convincing data on the integration of viral sequences into diverse organisms. Here, we identified 101 viral domains integrated into the model plant Arabidopsis proteome. Functional analysis based on gene ontology (GO) terms indicates that viral domains in the Arabidopsis proteome were involved in various stress responses with binding functions. Protein interaction networks support the strong protein interactions of viral domains with other Arabidopsis proteins. A proteome-wide analysis gave a comprehensive evolutionary view of viral domains integrated into 41 plant proteomes, revealing the specific and conserved integration of viral domains into plant proteomes. Phylogenetic analyses revealed the possible HGT between viral domains and plant proteomes. Our results provide an overview of the integration of viral domains into plant proteomes and their possible functional roles associated with plant defense mechanisms.

The pepper virome: natural co-infection of diverse viruses and their quasispecies.

BACKGROUND: The co-infection of diverse viruses in a host plant is common; however, little is known about viral populations and their quasispecies in the host. RESULTS: Here, we report the first pepper viromes that were co-infected by different types of viral genomes. The pepper viromes are dominated by geminivirus DNA-A followed by a novel carlavirus referred to as Pepper virus A. The two pepper cultivars share similar viral populations and replications. However, the quasispecies for double-stranded RNA virus and two satellite DNAs were heterogeneous and homogenous in susceptible and resistant cultivars, respectively, indicating the quasispecies of an individual virus depends on the host. CONCLUSIONS: Taken together, we provide the first evidence that the host plant resistant to viruses has an unrevealed antiviral system, affecting viral quasispecies, not replication.

Integrated analyses using RNA-Seq data reveal viral genomes, single nucleotide variations, the phylogenetic relationship, and recombination for Apple stem grooving virus.

BACKGROUND: Next-generation sequencing (NGS) provides many possibilities for plant virology research. In this study, we performed integrated analyses using plant transcriptome data for plant virus identification using Apple stem grooving virus (ASGV) as an exemplar virus. We used 15 publicly available transcriptome libraries from three different studies, two mRNA-Seq studies and a small RNA-Seq study. RESULTS: We de novo assembled nearly complete genomes of ASGV isolates Fuji and Cuiguan from apple and pear transcriptomes, respectively, and identified single nucleotide variations (SNVs) of ASGV within the transcriptomes. We demonstrated the application of NGS raw data to confirm viral infections in the plant transcriptomes. In addition, we compared the usability of two de novo assemblers, Trinity and Velvet, for virus identification and genome assembly. A phylogenetic tree revealed that ASGV and Citrus tatter leaf virus (CTLV) are the same virus, which was divided into two clades. Recombination analyses identified six recombination events from 21 viral genomes. CONCLUSIONS: Taken together, our in silico analyses using NGS data provide a successful application of plant transcriptomes to reveal extensive information associated with viral genome assembly, SNVs, phylogenetic relationships, and genetic recombination.

Comparative analysis of chrysanthemum transcriptome in response to three RNA viruses: Cucumber mosaic virus, Tomato spotted wilt virus and Potato virus X.

The chrysanthemum is one of popular flowers in the world and a host for several viruses. So far, molecular interaction studies between the chrysanthemum and viruses are limited. In this study, we carried out a transcriptome analysis of chrysanthemum in response to three different viruses including Cucumber mosaic virus (CMV), Tomato spotted wilt virus (TSWV) and Potato virus X (PVX). A chrysanthemum 135K microarray derived from expressed sequence tags was successfully applied for the expression profiles of the chrysanthemum at early stage of virus infection. Finally, we identified a total of 125, 70 and 124 differentially expressed genes (DEGs) for CMV, TSWV and PVX, respectively. Many DEGs were virus specific; however, 33 DEGs were commonly regulated by three viruses. Gene ontology (GO) enrichment analysis identified a total of 132 GO terms, and of them, six GO terms related stress response and MCM complex were commonly identified for three viruses. Several genes functioning in stress response such as chitin response and ethylene mediated signaling pathway were up-regulated indicating their involvement in establishment of host immune system. In particular, TSWV infection significantly down-regulated genes related to DNA metabolic process including DNA replication, chromatin organization, histone modification and cytokinesis, and they are mostly targeted to nucleosome and MCM complex. Taken together, our comparative transcriptome analysis revealed several genes related to hormone mediated viral stress response and DNA modification. The identified chrysanthemums genes could be good candidates for further functional study associated with resistant to various plant viruses.