Transcriptional Changes in Damask Rose Suspension Cell Culture Revealed by RNA Sequencing

Damask roses (Rosa x damascena) are widely used in cosmetics and pharmaceutics. Here, we established an in vitro suspension cell culture for calli derived from damask rose petals. We analyzed rose suspension cell transcriptomes obtained at two different time points by RNA sequencing to reveal transcriptional changes during rose suspension cell culture. Of the 580 coding RNAs (1.3%) highly expressed in the suspension rose cells, 68 encoded cell wall-associated proteins. However, most RNAs encoded by the chloroplasts and mitochondria are not expressed. Many highly expressed coding RNAs are involved in translation, catalyzing peptide synthesis in ribosomes. Moreover, the amide metabolic process producing naturally occurring alkaloids was the most abundant metabolic process during the propagation of rose suspension cells. During rose cell propagation, coding RNAs involved in the stress response were upregulated at an early stage, while coding RNAs associated with detoxification and transmembrane transport were upregulated at the late stage. We used transcriptome analyses to reveal important biological processes and molecular mechanisms during rose suspension cell culture. Most non-coding (nc) RNAs were not expressed in the rose suspension cells, but a few ncRNAs with unknown functions were highly expressed. The expression of ncRNAs and their target coding RNAs was highly correlated. Taken together, we revealed significant biological processes and molecular mechanisms occurring during rose suspension cell culture using transcriptome analyses.

Exploring Tomato Fruit Viromes through Transcriptome Data Analysis

This study delves into the complex landscape of viral infections in tomatoes (Solanum lycopersicum) using available transcriptome data. We conducted a virome analysis, revealing 219 viral contigs linked to four distinct viruses: tomato chlorosis virus (ToCV), southern tomato virus (STV), tomato yellow leaf curl virus (TYLCV), and cucumber mosaic virus (CMV). Among these, ToCV predominated in contig count, followed by STV, TYLCV, and CMV. A notable finding was the prevalence of coinfections, emphasizing the concurrent presence of multiple viruses in tomato plants. Despite generally low viral levels in fruit transcriptomes, STV emerged as the primary virus based on viral read count. We delved deeper into viral abundance and the contributions of RNA segments to replication. While initially focused on studying the impact of sound treatment on tomato fruit transcriptomes, the unexpected viral presence underscores the importance of considering viruses in plant research. Geographical variations in virome communities hint at potential forensic applications. Phylogenetic analysis provided insights into viral origins and genetic diversity, enhancing our understanding of the Korean tomato virome. In conclusion, this study advances our knowledge of the tomato virome, stressing the need for robust pest control in greenhouse-grown tomatoes and offering insights into virus management and crop protection.

Reappraisal of an Adequate Negative Margin Following Breast Conservation Therapy in Young Patients with Invasive Breast Cancer

PURPOSE/OBJECTIVE(S)

This study aimed to examine the relationship between resection margin (RM) width and ipsilateral breast tumor recurrence (IBTR). The current definition of an adequate negative margin was validated in young breast cancer patients (<40 years).

MATERIALS/METHODS

This study included 4,042 patients who underwent breast conservation therapy between 2006 and 2012. There were 595 women <40 years (14.7%). Re-excision to achieve a wider RM was not considered in patients with no ink on tumor. Systemic treatments were administered in 3,985 patients (98.6%). The impact of RM width on the incidence of IBTR was assessed using narrowly specified intervals (positive, >0 to <1 mm, 1 to <2 mm, 2 to <3 mm, and ≥3 mm).

RESULTS

A significant relationship was observed between RM width and incidence of IBTR. In young women, the 10-year cumulative incidence rates of IBTR were 21.6% for positive RM, 13.9% for >0 to <1 mm, 10.9% for 1 to <2 mm, 5.7% for 2 to <3 mm, and 5.8% for ≥3 mm (p <0.01). The corresponding incidence rates of IBTR in women ≥40 years were 8.6%, 8.0%, 4.4%, 2.7%, and 2.1%, respectively (p <0.01). The difference in the incidence of IBTR between each RM width interval increased over time, and the magnitude of the increment was greater in young women. RM width was the only risk factor for IBTR in young women.

CONCLUSION

The relationship between RM width and IBTR was definite in young women. Compared with women ≥40 years, young women might have long-term benefits from a wider RM than no ink on tumor.

Investigating Variability in Viral Presence and Abundance across Soybean Seed Development Stages Using Transcriptome Analysis

Plant transcriptomes offer a valuable resource for studying viral communities (viromes). In this study, we explore how plant transcriptome data can be applied to virome research. We analyzed 40 soybean transcriptomes across different growth stages and identified six viruses: broad bean wilt virus 2 (BBWV2), brassica yellow virus (BrYV), beet western yellow virus (BWYV), cucumber mosaic virus (CMV), milk vetch dwarf virus (MDV), and soybean mosaic virus (SMV). SMV was the predominant virus in both Glycine max (GM) and Glycine soja (GS) cultivars. Our analysis confirmed its abundance in both, while BBWV2 and CMV were more prevalent in GS than GM. The viral proportions varied across developmental stages, peaking in open flowers. Comparing viral abundance measured by viral reads and fragments per kilobase of transcript per million (FPKM) values revealed insights. SMV showed similar FPKM values in GM and GS, but BBWV2 and CMV displayed higher FPKM proportions in GS. Notably, the differences in viral abundance between GM and GS were generally insignificant based on the FPKM values across developmental stages, except for the apical bud stage in four GM cultivars. We also detected MDV, a multi-segmented virus, in two GM samples, with variable proportions of its segments. In conclusion, our study demonstrates the potential of plant transcriptomes for virome research, highlighting their strengths and limitations.

In Silico Virome Analysis of Chinese Narcissus Transcriptomes Reveals Diverse Virus Species and Genetic Diversity at Different Flower Development Stages

Viromes of Chinese narcissus flowers were explored using transcriptome data from 20 samples collected at different flower development stages. Quality controlled raw data underwent de novo assembly, resulting in 5893 viral contigs that matched the seven virus species. The most abundant viruses were narcissus common latent virus (NCLV), narcissus yellow stripe virus (NYSV), and narcissus mottling-associated virus (NMaV). As flower development stages advanced, white tepal plants showed an increase in the proportion of viral reads, while the variation in viral proportion among yellow tepal plants was relatively small. Narcissus degeneration virus (NDV) dominated the white tepal samples, whereas NDV and NYSV prevailed in the yellow tepal samples. Potyviruses, particularly NDV, are the primary infectious viruses. De novo assembly generated viral contigs for five viruses, yielding complete genomes for NCLV, NDV, narcissus late season yellow virus (NLSYV), and NYSV. Phylogenetic analysis revealed genetic diversity, with distinct NCLV, NMaV, NDV, NLSYV, and NYSV groups. This study provides valuable insights into the viromes and genetic diversity of viruses in Chinese narcissus flowers.

Immune-checkpoint proteins, cytokines, and microbiome impact on patients with cervical insufficiency and preterm birth

Background

Microenvironmental factors, including microbe-induced inflammation and immune-checkpoint proteins that modulate immune cells have been associated with both cervical insufficiency and preterm delivery. These factors are incompletely understood. This study aimed to explore and compare interactions among microbiome and inflammatory factors, such as cytokines and immune-checkpoint proteins, in patients with cervical insufficiency and preterm birth. In particular, factors related to predicting preterm birth were identified and the performance of the combination of these factors was evaluated.

Methods

A total of 220 swab samples from 110 pregnant women, prospectively recruited at the High-Risk Maternal Neonatal Intensive Care Center, were collected between February 2020 and March 2021. This study included 63 patients with cervical insufficiency receiving cerclage and 47 control participants. Endo- and exocervical swabs and fluids were collected simultaneously. Shotgun metagenomic sequencing for the microbiome and the measurement of 34 immune-checkpoint proteins and inflammatory cytokines were performed.

Results

First, we demonstrated that immune-checkpoint proteins, the key immune-regulatory molecules, could be measured in endocervical and exocervical samples. Secondly, we identified significantly different microenvironments in cervical insufficiency and preterm birth, with precise cervical locations, to provide information about practically useful cervical locations in clinical settings. Finally, the presence of Moraxella osloensis (odds ratio = 14.785; P = 0.037) and chemokine CC motif ligand 2 levels higher than 73 pg/mL (odds ratio = 40.049; P = 0.005) in endocervical samples were associated with preterm birth. Combining M. osloensis and chemokine CC motif ligand 2 yielded excellent performance for predicting preterm birth (area under the receiver operating characteristic curve = 0.846, 95% confidence interval = 0.733-0.925).

Conclusion

Multiple relationships between microbiomes, immune-checkpoint proteins, and inflammatory cytokines in the cervical microenvironment were identified. We focus on these factors to aid in the comprehensive understanding and therapeutic modulation of local microbial and immunologic compositions for the management of cervical insufficiency and preterm birth.

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.

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.

Phylogenetic and Phylodynamic Analyses of Soybean Mosaic Virus Using 305 Coat Protein Gene Sequences

Soybean mosaic virus (SMV) of the family Potyviridae is the most devastating virus that infects soybean plants. In this study, we obtained 83 SMV coat protein (CP) sequences from seven provinces in Korea using RT-PCR and Sanger sequencing. Phylogenetic and haplotype analyses revealed eight groups of 83 SMV isolates and a network of 50 SMV haplotypes in Korea. The phylogenetic tree using 305 SMV CP sequences available worldwide revealed 12 clades that were further divided into two groups according to the plant hosts. Recombination rarely occurred in the CP sequences, while negative selection was dominant in the SMV CP sequences. Genetic diversity analyses revealed that plant species had a greater impact on the genetic diversity of SMV CP sequences than geographical origin or location. SMV isolates identified from Pinellia species in China showed the highest genetic diversity. Phylodynamic analysis showed that the SMV isolates between the two Pinellia species diverged in the year 1248. Since the divergence of the first SMV isolate from Glycine max in 1486, major clades for SMV isolates infecting Glycine species seem to have diverged from 1791 to 1886. Taken together, we provide a comprehensive overview of the genetic diversity and divergence of SMV CP sequences.

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.

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.

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.

Analysis of the transcriptomic, metabolomic, and gene regulatory responses to Puccinia sorghi in maize

Common rust, caused by Puccinia sorghi, is a widespread and destructive disease of maize. The Rp1-D gene confers resistance to the P. sorghi IN2 isolate, mediating a hypersensitive cell death response (HR). To identify differentially expressed genes (DEGs) and metabolites associated with the compatible (susceptible) interaction and with Rp1-D-mediated resistance in maize, we performed transcriptomics and targeted metabolome analyses of P. sorghi IN2-infected leaves from the near-isogenic lines H95 and H95:Rp1-D, which differed for the presence of Rp1-D. We observed up-regulation of genes involved in the defence response and secondary metabolism, including the phenylpropanoid, flavonoid, and terpenoid pathways. Metabolome analyses confirmed that intermediates from several transcriptionally up-regulated pathways accumulated during the defence response. We identified a common response in H95:Rp1-D and H95 with an additional H95:Rp1-D-specific resistance response observed at early time points at both transcriptional and metabolic levels. To better understand the mechanisms underlying Rp1-D-mediated resistance, we inferred gene regulatory networks occurring in response to P. sorghi infection. A number of transcription factors including WRKY53, BHLH124, NKD1, BZIP84, and MYB100 were identified as potentially important signalling hubs in the resistance-specific response. Overall, this study provides a novel and multifaceted understanding of the maize susceptible and resistance-specific responses to P. sorghi.

Gene expression profile of human follicle dermal papilla cells in response to Camellia japonica phytoplacenta extract

Camellia japonica L. is a flowering tree with several medicinal and cosmetic applications. Here, we investigated the efficacy of C. japonica placenta extract (CJPE) as a potential therapeutic agent for promotion of hair growth and scalp health by using various in vitro and in vivo assays. Moreover, we performed transcriptome analysis to examine the relative expression of human follicle dermal papilla cells (HFDPC) in response to CJPE by RNA-sequencing (RNA-seq). In vitro assays revealed upregulation of the expression of hair growth marker genes in HFDPC after CJPE treatment. Moreover, in vivo clinical tests with 42 adult female participants showed that a solution containing 0.5% CJPE increased the moisture content of the scalp and decreased the scalp's sebum content, dead scalp keratin, and erythema. Furthermore, RNA-seq analysis revealed key genes in HFDPC which are associated with CJPE. Interestingly, genes associated with lipid metabolism and cholesterol efflux were upregulated. Genes upregulated by CJPE are associated with several hormones, including parathyroid, adrenocorticotropic hormone, α-melanocyte-stimulating hormone (alpha-MSH), and norepinephrine, which are involved in hair follicle biology. Furthermore, some upregulated genes are associated with the regulation of axon guidance. In contrast, many genes downregulated by CJPE are associated with structural components of the cytoskeleton. In addition, CJPE suppressed genes associated with muscle structure and development. Taken together, this study provides extensive evidence that CJPE may have potential as a therapeutic agent for scalp treatment and hair growth promotion.

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.

Comparative Microbiome Study of Mummified Peach Fruits by Metagenomics and Metatranscriptomics

The dried peach fruits clinging to peach trees or lying on the ground nearby are known as mummified peach fruits. Here, we examined the microbiome communities of three different mummified peach fruits from the nectarine cultivar "Hahong" by DNA- and RNA-sequencing. We found the dominance of Monilinia fructigena followed by Sclerotinia borealis, S. sclerotiorum, and Botrytis cinerea in the mummified peach fruits. Moreover, we found a high number of Proteobacteria, including Frateuria aurantia, Neoasaia chiangmaiensis, Robbsia andropogonis, and Ewingella Americana. Furthermore, we identified several viruses and viroids. Bacteriophages were identified by DNA- and RNA-sequencing, while viruses and viroids with RNA genomes were identified by only RNA-sequencing. Moreover, we identified a novel mycovirus referred to as Monilinia umbra-like virus 1 (MULV1) from M. fructigena. Our results revealed the co-inhabitance of fungi and bacteria in the mummified peach fruits, although dominant microorganisms were present. RNA-sequencing revealed that several fungal and bacterial genes were actively transcribed. Comparative analyses suggested that RNA-sequencing provides more detailed information on microbial communities; however, combining DNA- and RNA-sequencing results increased the diversity of microorganisms, suggesting the importance of databases and analysis tools for microbiome studies. Taken together, our study provides a comprehensive overview of microbial communities in mummified peach fruits by DNA shotgun sequencing and RNA-sequencing.

Identification of viruses infecting six plum cultivars in Korea by RNA-sequencing

BACKGROUND

Plums are a kind of stone fruit, a category that includes peaches, cherries, apricots, and almonds. In Korea, Japanese plum trees are usually cultivated as they best suit the climate. To date, there have been few studies in Korea on viruses infecting plum trees compared to those infecting peach trees.

METHODS

To identify viruses and viroids infecting plum trees, we collected leaf samples from six different plum cultivars and subjected them to RNA-sequencing (RNA-seq). Six different plum transcriptomes were de novo assembled using the Trinity assembler followed by BLAST searching against a viral reference database.

RESULTS

We identified hop stunt viroid (HSVd) and six viruses, including apple chlorotic leaf spot virus (ACLSV), little cherry virus-1 (LChV-1), peach virus D (PeVD), peach leaf pitting-associated virus (PLPaV), plum bark necrosis stem pitting-associated virus (PBNSPaV), and prunus necrotic ringspot virus (PNRSV), from six plum cultivars by RNA-seq. RT-PCR confirmed the infection of HSVd and three viruses-ACLSV, PBNSPaV, and PNRSV-in plum trees. However, RT-PCR demonstrated that plum trees in this study were not infected by LChV-1, PeVD, or PLPaV. It is likely that the three viruses LChV-1, PeVD, and PLPaV as identified by RNA-seq were contaminants from other peach libraries caused by index misassignment, which suggests that careful confirmation by other methods should be carried out in next-generation sequencing (NGS)-based virus identification. Taken together, we identified a viroid and three viruses infecting plum trees in Korea.

Anti-Aging Effects of Leontopodium alpinum (Edelweiss) Callus Culture Extract Through Transcriptome Profiling

Edelweiss (Leontopodium Alpinum) in the family Asteraceae is a wildflower that grows in rocky limestone places. Here, we investigated the efficacy of edelweiss callus culture extract (Leontopodium Alpinum callus culture extract; LACCE) using multiple assays from in vitro to in vivo as well as transcriptome profiling. Several in vitro assay results showed the strong antioxidant activity of LACCE in response to UVB treatment. Moreover, LACCE suppressed inflammation and wrinkling; however, moisturizing activity was increased by LACCE. The clinical test in vivo demonstrated that constant application of LACCE on the face and skin tissues improved anti-periorbital wrinkles, skin elasticity, dermal density, and skin thickness compared with the placebo. The RNA-Sequencing results showed at least 16.56% of human genes were expressed in keratinocyte cells. LACCE up-regulated genes encoding several KRT proteins; DDIT4, BNIP3, and IGFBP3 were involved in the positive regulation of the developmental process, programmed cell death, keratinization, and cornification forming skin barriers, which provide many advantages in the human skin. By contrast, down-regulated genes were stress-responsive genes, including metal, oxidation, wounding, hypoxia, and virus infection, suggesting LACCE did not cause any harmful stress on the skin. Our comprehensive study demonstrated LACCE is a promising agent for anti-aging cosmetics.

Sweet potato viromes in eight different geographical regions in Korea and two different cultivars

The sweet potato in the family Convolvulaceae is a dicotyledonous perennial plant. Here, we conducted a comprehensive sweet potato virome study using 10 different libraries from eight regions in Korea and two different sweet potato cultivars by RNA-Sequencing. Comprehensive bioinformatics analyses revealed 10 different virus species infecting sweet potato. Moreover, we identified two novel viruses infecting sweet potato referred to as Sweet potato virus E (SPVE) in the genus Potyvirus and Sweet potato virus F (SPVF) in the genus Carlavirus. Of the identified viruses, Sweet potato feathery mottle virus (SPFMV) was the dominant virus followed by Sweet potato virus C (SPVC) and SPVE in Korea. We obtained a total of 30 viral genomes for eight viruses. Our phylogenetic analyses showed many potyvirus isolates are highly correlated with geographical regions. However, two isolates of SPFMV and a single isolate of Sweet potato virus G (SPVG) were genetically distant from other known isolates. The mutation rate was the highest in SPFMV followed by SPVC and SPVG. Two different sweet potato cultivars, Beni Haruka and Hogammi, were infected by seven and five viruses, respectively. Taken together, we provide a complete list of viruses infecting sweet potato in Korea and diagnostic methods.

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.

Genome-wide identification of long non-coding RNAs in tomato plants irradiated by neutrons followed by infection with Tomato yellow leaf curl virus

Long non-coding RNAs (lncRNAs) play an important role in regulating many biological processes. In this study, tomato seeds were first irradiated by neutrons. Eight tomato mutants were then selected and infected by Tomato yellow leaf curl virus (TYLCV). RNA sequencing followed by bioinformatics analyses identified 1,563 tomato lncRNAs. About half of the lncRNAs were derived from intergenic regions, whereas antisense lncRNAs accounted for 35%. There were fewer lncRNAs identified in our study than in other studies identifying tomato lncRNAs. Functional classification of 794 lncRNAs associated with tomato genes showed that many lncRNAs were associated with binding functions required for interactions with other molecules and localized in the cytosol and membrane. In addition, we identified 19 up-regulated and 11 down-regulated tomato lncRNAs by comparing TYLCV infected plants to non-infected plants using previously published data. Based on these results, the lncRNAs identified in this study provide important resources for characterization of tomato lncRNAs in response to TYLCV infection.

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.

Barley RNA viromes in six different geographical regions in Korea

Barley is a kind of cereal grass belonging to the family Poaceae. To examine viruses infecting winter barley in Korea, we carried out a comprehensive study of barley RNA viromes using next-generation sequencing (NGS). A total of 110 barley leaf samples from 17 geographical locations were collected. NGS followed by extensive bioinformatics analyses revealed six different barley viromes: Barley yellow mosaic virus (BaYMV), Barley mild mosaic virus (BaMMV), Barley yellow dwarf virus (BYDV), Hordeum vulgare endornavirus (HvEV), and Barley virus G (BVG). BaYMV and HvEV were identified in all libraries, while other viruses were identified in some specific library. Based on the number of virus-associated reads, BaYMV was a dominant virus infecting winter barley in Korea causing yellow disease symptoms. We obtained nearly complete genomes of six BaYMV isolates and two BaMMV isolates. Phylogenetic analyses indicate that BaYMV and BaMMV were largely grouped based on geographical regions such as Asia and Europe. Single nucleotide polymorphisms analyses suggested that most BaYMV and BaMMV showed strong genetic variations; however, BaYMV isolate Jeonju and BaMMV isolate Gunsan exhibited a few and no SNPs, respectively, suggesting low level of genetic variation. Taken together, this is the first study of barley RNA viromes in Korea.

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.

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.

De novo Genome Assembly and Single Nucleotide Variations for Soybean Mosaic Virus Using Soybean Seed Transcriptome Data

Soybean is the most important legume crop in the world. Several diseases in soybean lead to serious yield losses in major soybean-producing countries. Moreover, soybean can be infected by diverse viruses. Recently, we carried out a large-scale screening to identify viruses infecting soybean using available soybean transcriptome data. Of the screened transcriptomes, a soybean transcriptome for soybean seed development analysis contains several virus-associated sequences. In this study, we identified five viruses, including soybean mosaic virus (SMV), infecting soybean by de novo transcriptome assembly followed by blast search. We assembled a nearly complete consensus genome sequence of SMV China using transcriptome data. Based on phylogenetic analysis, the consensus genome sequence of SMV China was closely related to SMV isolates from South Korea. We examined single nucleotide variations (SNVs) for SMVs in the soybean seed transcriptome revealing 780 SNVs, which were evenly distributed on the SMV genome. Four SNVs, C-U, U-C, A-G, and G-A, were frequently identified. This result demonstrated the quasispecies variation of the SMV genome. Taken together, this study carried out bioinformatics analyses to identify viruses using soybean transcriptome data. In addition, we demonstrated the application of soybean transcriptome data for virus genome assembly and SNV analysis.

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.

Identification of residues or motif(s) of the rice stripe virus NS3 protein required for self-interaction and for silencing suppressor activity

Rice stripe virus (RSV) is an important pathogen of rice. The RSV genome consists of four single-stranded RNA segments that encode seven viral proteins. A previous report found that NS3 is a viral suppressor of RNA silencing and self interacts. Using a model that predicts protein structure, we identified amino acid residues or motifs, including four α-helix motifs, required for NS3 self-interaction. We then used yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays to study the interactions between full-length NS3 and its truncated and alanine substitution mutants. Y2H and BiFC results showed that the N-terminal region of NS3 is essential for self-interaction. All α-helix deletion mutants and substitution mutants lost the ability to self-interact. To identify the relationship between NS3 self-interaction and silencing suppressor activity, we used a GFP silencing system in Nicotiana benthamiana with Agrobacterium-mediated transient overexpression of each mutated NS3 protein. All of the deletion and the α-helix substitution mutants that had lost the ability to self-interact also lost their silencing suppressor ability. The substitution of amino acids with alanine at positions 70-75, 76-83, and 173-177, however, resulted in mutants that were able to self-interact but were unable to function as silencing suppressors. These results suggest that RSV requires NS3 self-interaction to suppress RNA silencing and to thereby counter host defenses.

Sequence variability of Chrysanthemum stunt viroid in different chrysanthemum cultivars

Viroids are the smallest infectious agents, and their genomes consist of a short single strand of RNA that does not encode any protein. Chrysanthemum stunt viroid (CSVd), a member of the family Pospiviroidae, causes chrysanthemum stunt disease. Here, we report the genomic variations of CSVd to understand the sequence variability of CSVd in different chrysanthemum cultivars. We randomly sampled 36 different chrysanthemum cultivars and examined the infection of CSVd in each cultivar by reverse transcription polymerase chain reaction (RT-PCR). Eleven cultivars were infected by CSVd. Cloning followed by Sanger sequencing successfully identified a total of 271 CSVd genomes derived from 12 plants from 11 cultivars. They were further classified into 105 CSVd variants. Each single chrysanthemum plant had a different set of CSVd variants. Moreover, different single plants from the same cultivar had different sets of CSVd variants but identical consensus genome sequences. A phylogenetic tree using 12 consensus genome sequences revealed three groups of CSVd genomes, while six different groups were defined by the phylogenetic analysis using 105 variants. Based on the consensus CSVd genome, by combining all variant sequences, we identified 99 single-nucleotide variations (SNVs) as well as three nucleotide positions showing high mutation rates. Although 99 SNVs were identified, most CSVd genomes in this study were derived from variant 1, which is identical to known CSVd SK1 showing pathogenicity.

Time-Course Small RNA Profiling Reveals Rice miRNAs and Their Target Genes in Response to Rice Stripe Virus Infection

It has been known that many microRNAs (miRNAs) are involved in the regulation for the plant development and defense mechanism by regulating the expression of the target gene. Several previous studies has demonstrated functional roles of miRNAs in antiviral defense mechanisms. In this study, we employed high-throughput sequencing technology to identify rice miRNAs upon rice stripe virus (RSV) infection at three different time points. Six libraries from mock and RSV-infected samples were subjected for small RNA sequencing. Bioinformatic analyses revealed 374 known miRNAs and 19 novel miRNAs. Expression of most identified miRNAs was not dramatically changed at 3 days post infection (dpi) and 7 dpi by RSV infection. However, many numbers of miRNAs were up-regulated in mock and RSV-infected samples at 15 dpi by RSV infection. Moreover, expression profiles of identified miRNAs revealed that only few numbers of miRNAs were strongly regulated by RSV infection. In addition, 15 resistance genes were targets of six miRNAs suggesting that those identified miRNAs and 15 NBS-LRR resistance genes might be involved in RSV infection. Taken together, our results provide novel insight into the dynamic expression profiles of rice miRNAs upon RSV infection and clues for the understanding of the regulatory roles of miRNAs via time-course.

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2994-6) contains supplementary material, which is available to authorized users.

Development of Virus-Induced Gene Expression and Silencing Vector Derived from Grapevine Algerian Latent Virus

Grapevine Algerian latent virus (GALV) is a member of the genus Tombusvirus in the Tombusviridae and infects not only woody perennial grapevine plant but also herbaceous Nicotiana benthamiana plant. In this study, we developed GALV-based gene expression and virus-induced gene silencing (VIGS) vectors in N. benthamiana. The GALV coat protein deletion vector, pGMG, was applied to express the reporter gene, green fluorescence protein (GFP), but the expression of GFP was not detected due to the necrotic cell death on the infiltrated leaves. The p19 silencing suppressor of GALV was engineered to inactivate its expression and GFP was successfully expressed with unrelated silencing suppressor, HC-Pro, from soybean mosaic virus. The pGMG vector was used to knock down magnesium chelatase (ChlH) gene in N. benthamaina and the silencing phenotype was clearly observed on systemic leaves. Altogether, the GALV-derived vector is expected to be an attractive tool for useful gene expression and VIGS vectors in grapevine as well as N. benthamiana.

Clinical impact of coexisting retinopathy and vascular calcification on chronic kidney disease progression and cardiovascular events

BACKGROUND AND AIMS

Retinopathy and vascular calcification (VC) are representative markers of microvascular and macrovascular dysfunction in patients with chronic kidney disease (CKD). However, their relationship and combined effects on clinical outcomes remain undetermined.

METHODS AND RESULTS

We included 523 patients with nondialysis-dependent CKD stage 3-5 who had been examined with fundus photography for diabetic or hypertensive retinopathy. Simple radiographs were analyzed for the presence of VC. The clinical significance of VC of the abdominal aorta and iliofemoral artery (apVC) and retinopathy was evaluated in terms of the rate of renal function decline and composite of any cardiovascular event or death. CKD patients with retinopathy showed higher prevalence of apVC than those without retinopathy (25.6% vs. 12.5%, P < 0.001).The presence of retinopathy was independently associated with apVC (OR 2.13, 95% CI 1.31, 3.49). In multivariate analysis, compared with subjects with neither apVC nor retinopathy, the coexistence of both apVC and retinopathy were independently associated with rapid renal function decline (β = -1.51; 95% CI -2.40, -0.61), whereas apVC or retinopathy alone were not. Compared with subjects with neither apVC nor retinopathy, the HRs for composite end points were 1.05 (95% CI 0.48, 2.27), 1.79 (95% CI 1.14, 2.80), and 2.07 (95% CI 1.17, 3.67) for patients with apVC only, those with retinopathy only, and those with both apVC and retinopathy, respectively.

CONCLUSION

The coexistence of VC and retinopathy was independently associated with CKD progression and cardiovascular events or deaths, and its combined effect was stronger than any separate condition.

Two homologous host proteins interact with potato virus X RNAs and CPs and affect viral replication and movement

Because viruses encode only a small number of proteins, all steps of virus infection rely on specific interactions between viruses and hosts. We previously screened several Nicotiana benthamiana (Nb) proteins that interact with the stem-loop 1 (SL1) RNA structure located at the 5' end of the potato virus X (PVX) genome. In this study, we characterized two of these proteins (NbCPIP2a and NbCPIP2b), which are homologous and are induced upon PVX infection. Electrophoretic mobility shift assay confirmed that both proteins bind to either SL1(+) or SL1(-) RNAs of PVX. The two proteins also interact with the PVX capsid protein (CP) in planta. Overexpression of NbCPIP2a positively regulated systemic movement of PVX in N. benthamiana, whereas NbCPIP2b overexpression did not affect systemic movement of PVX. Transient overexpression and silencing experiments demonstrated that NbCPIP2a and NbCPIP2b are positive regulators of PVX replication and that the effect on replication was greater for NbCPIP2a than for NbCPIP2b. Although these two host proteins are associated with plasma membranes, PVX infection did not affect their subcellular localization. Taken together, these results indicate that NbCPIP2a and NbCPIP2b specifically bind to PVX SL1 RNAs as well as to CP and enhance PVX replication and movement.

De novo transcriptome assembly of two Vigna angularis varieties collected from Korea

The adzuki bean (Vigna angularis), a member of the family Fabaceae, is widely grown in Asia, from East Asia to the Himalayas. The adzuki bean is known as an ingredient that adds sweetness to diverse desserts made in Eastern Asian countries. Libraries prepared from two V. angularis varieties referred to as Taejin Black and Taejin Red were paired-end sequenced using the Illumina HiSeq 2000 system. The raw data in this study can be available in NCBI SRA database with accession numbers of SRR3406660 and SRR3406553. After de novo transcriptome assembly using Trinity, we obtained 324,219 and 280,056 transcripts from Taejin Black and Taejin Red, respectively. We predicted a total of 238,321 proteins and 179,519 proteins for Taejin Black and Taejin Red, respectively, by the TransDecoder program. We carried out BLASTP on the predicted proteins against the Swiss-Prot protein sequence database to predict the putative functions of identified proteins. Taken together, we provide transcriptomes of two adzuki bean varieties by RNA-Seq, which might be usefully applied to generate molecular markers.

De novo transcriptome assembly of Sorghum bicolor variety Taejin

Sorghum (Sorghum bicolor), also known as great millet, is one of the most popular cultivated grass species in the world. Sorghum is frequently consumed as food for humans and animals as well as used for ethanol production. In this study, we conducted de novo transcriptome assembly for sorghum variety Taejin by next-generation sequencing, obtaining 8.748 GB of raw data. The raw data in this study can be available in NCBI SRA database with accession number of SRX1715644. Using the Trinity program, we identified 222,161 transcripts from sorghum variety Taejin. We further predicted coding regions within the assembled transcripts by the TransDecoder program, resulting in a total of 148,531 proteins. We carried out BLASTP against the Swiss-Prot protein sequence database to annotate the functions of the identified proteins. To our knowledge, this is the first transcriptome data for a sorghum variety derived from Korea, and it can be usefully applied to the generation of genetic markers.

De novo transcriptome assembly of Setatria italica variety Taejin

Foxtail millet (Setaria italica) belonging to the family Poaceae is an important millet that is widely cultivated in East Asia. Of the cultivated millets, the foxtail millet has the longest history and is one of the main food crops in South India and China. Moreover, foxtail millet is a model plant system for biofuel generation utilizing the C₄ photosynthetic pathway. In this study, we carried out de novo transcriptome assembly for the foxtail millet variety Taejin collected from Korea using next-generation sequencing. We obtained a total of 8.676 GB raw data by paired-end sequencing. The raw data in this study can be available in NCBI SRA database with accession number of SRR3406552. The Trinity program was used to de novo assemble 145,332 transcripts. Using the TransDecoder program, we predicted 82,925 putative proteins. BLASTP was performed against the Swiss-Prot protein sequence database to annotate the functions of identified proteins, resulting in 20,555 potentially novel proteins. Taken together, this study provides transcriptome data for the foxtail millet variety Taejin by RNA-Seq.

De novo transcriptome assembly of two different Prunus mume cultivars

Prunus mume, belonging to the Prunus genus, is an Asian tree, and its common names are Chinese plum and Japanese plum. P. mume are cultivated for fruit production as well as ornamental purposes. In this study, we conducted de novo transcriptome assembly for two selected P. mume cultivars referred to as Takada and Wallyoung (commercially important cultivars for fruit production and ornamental trees, respectively) by RNA-sequencing. We obtained 9.14 GB and 9.48 GB sequence data from Takada and Wallyoung (NCBI accession numbers: SRX1187101 and SRX1187169), respectively. De novo transcriptome assembly identified 130,989 and 116,941 transcripts for Takada and Wallyoung, respectively. In addition, we identified 96,681 and 91,429 proteins from Takada and Wallyoung, respectively, by TransDecoder program. We performed BLASTP against the NCBI non-redundant (nr) datasets to annotate identified proteins. This study provides transcriptomes and proteomes for two different P. mume cultivars, which might be useful for comparative transcriptome analyses and assist development of genetic markers.

De novo transcriptome assembly of two different Prunus salicina cultivars

Plum is a globally grown stone fruit and can be divided into several species. In particular, the Prunus salicina, which is native to China, is widely grown in many fruit orchards in Korea and Japan, as well as the United States and Australia. The transcriptome data for Prunus salicina has not been reported to our knowledge. In this study, we performed de novo transcriptome assembly for two selected P. salicina cultivars referred to as Akihime and Formosa (commercially important plum cultivars in Korea) using next generation sequencing. We obtained a total of 9.04 GB and 8.68 GB raw data from Akihime and Formosa, respectively. De novo transcriptome assembly using Trinity revealed 155,169 and 160,186 transcripts for Akihime and Formosa. Next, we identified 121,278 and 116,544 proteins from Akihime and Formosa using TransDecoder. We performed BLASTP against the NCBI non-redundant (nr) dataset to annotate proteins. Taken together, this is the first transcriptome data for P. salicina to our knowledge.

Time-Course RNA-Seq Analysis Reveals Transcriptional Changes in Rice Plants Triggered by Rice stripe virus Infection

Rice stripe virus (RSV) has become a major pathogen of rice. To determine how the rice transcriptome is modified in response to RSV infection, we used RNA-Seq to perform a genome-wide gene expression analysis of a susceptible rice cultivar. The transcriptomes of RSV-infected samples were compared to those of mock-treated samples at 3, 7, and 15 days post-infection (dpi). From 8 to 11% of the genes were differentially expressed (>2-fold difference in expression) in RSV-infected vs. noninfected rice. Among them, 532 genes were differentially expressed at all three time points. Surprisingly, 37.6% of the 532 genes are related to transposons. Gene ontology enrichment analysis revealed that many chloroplast genes were down-regulated in infected plants at 3 and 15 dpi. Expression of genes associated with cell differentiation and flowering was significantly down-regulated in infected plants at 15 dpi. In contrast, most of the up-regulated genes in infected plants concern the cell wall, plasma membrane, and vacuole and are known to function in various metabolic pathways and stress responses. In addition, transcripts of diverse transcription factors gradually accumulated in infected plants with increasing infection time. We also confirmed that the expression of gene subsets (including NBS-LRR domain-containing genes, receptor-like kinase genes, and genes involving RNA silencing) was changed by RSV infection. Taken together, we demonstrated that down-regulation of genes related to photosynthesis and flowering was strongly associated with disease symptoms caused by RSV and that up-regulation of genes involved in metabolic pathways, stress responses, and transcription was related to host defense mechanisms.

Proteomic Analysis to Identify Tightly-Bound Cell Wall Protein in Rice Calli

Rice is a model plant widely used for basic and applied research programs. Plant cell wall proteins play key roles in a broad range of biological processes. However, presently, knowledge on the rice cell wall proteome is rudimentary in nature. In the present study, the tightly-bound cell wall proteome of rice callus cultured cells using sequential extraction protocols was developed using mass spectrometry and bioinformatics methods, leading to the identification of 1568 candidate proteins. Based on bioinformatics analyses, 389 classical rice cell wall proteins, possessing a signal peptide, and 334 putative non-classical cell wall proteins, lacking a signal peptide, were identified. By combining previously established rice cell wall protein databases with current data for the classical rice cell wall proteins, a comprehensive rice cell wall proteome, comprised of 496 proteins, was constructed. A comparative analysis of the rice and Arabidopsis cell wall proteomes revealed a high level of homology, suggesting a predominant conservation between monocot and eudicot cell wall proteins. This study importantly increased information on cell wall proteins, which serves for future functional analyses of these identified rice cell wall proteins.

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.

Bacterial communities in the phylloplane of Prunus species

Bacterial populations in the phylloplane of four different Prunus species were investigated by 16 S rRNA pyrosequencing. Bioinformatic analysis identified an average of 510 operational taxonomic units belonging to 159 genera in 76 families. The two genera, Sphingomonas and Methylobacterium, were dominant in the phylloplane of four Prunus species. Twenty three genera were commonly identified in the four Prunus species, indicating a high level of bacterial diversity dependent on the plant species. Our study based on 16 S rRNA sequencing reveals the complexity of bacterial diversity in the phylloplane of Prunus species in detail.

Type 2C protein phosphatase is a key regulator of antiviral extreme resistance limiting virus spread

Effector-triggered immunity (ETI) is an active immune response triggered by interactions between host resistance proteins and their cognate effectors. Although ETI is often associated with the hypersensitive response (HR), various R genes mediate an HR-independent process known as extreme resistance (ER). In the soybean-Soybean mosaic virus (SMV) pathosystem, the strain-specific CI protein of SMV functions as an effector of Rsv3-mediated ER. In this study, we used the soybean (Rsv3)-SMV (CI) pathosystem to gain insight into the molecular signaling pathway involved in ER. We used genome-wide transcriptome analysis to identify a subset of the type 2C protein phophatase (PP2C) genes that are specifically up-regulated in Rsv3-mediated ER. Gain-of-function analysis of the most significantly expressed soybean PP2C gene, GmPP2C3a, showed that ABA-induced GmPP2C3a functions as a key regulator of Rsv3-mediated ER. Our results further suggest that the primary mechanism of ER against viruses is the inhibition of viral cell-to-cell movement by callose deposition in an ABA signaling-dependent manner.

Current Insights into Research on Rice stripe virus

Rice stripe virus (RSV) is one of the most destructive viruses of rice, and greatly reduces rice production in China, Japan, and Korea, where mostly japonica cultivars of rice are grown. RSV is transmitted by the small brown plant-hopper (SBPH) in a persistent and circulative-propagative manner. Several methods have been developed for detection of RSV, which is composed of four single-stranded RNAs that encode seven proteins. Genome sequence data and comparative phylogenetic analysis have been used to identify the origin and diversity of RSV isolates. Several rice varieties resistant to RSV have been selected and QTL analysis and fine mapping have been intensively performed to map RSV resistance loci or genes. RSV genes have been used to generate several genetically modified transgenic rice plants with RSV resistance. Recently, genome-wide transcriptome analyses and deep sequencing have been used to identify mRNAs and small RNAs involved in RSV infection; several rice host factors that interact with RSV proteins have also been identified. In this article, we review the current statues of RSV research and propose integrated approaches for the study of interactions among RSV, rice, and the SBPH.