Transcriptome-Wide Identification of WRKY Transcription Factor and Functional Characterization of RgWRKY37 Involved in Acteoside Biosynthesis in Rehmannia glutinosa

WRKYs play important roles in plant metabolism, but their regulation mechanism in Rehmannia glutinosa remains elusive. In this study, 37 putative WRKY transcription factors (TFs) with complete WRKY domain from R. glutinosa transcriptome sequence data were identified. Based on their conserved domains and zinc finger motif, the R. glutinosa WRKY TFs were divided into five groups. Structural feature analysis shows that the 37 RgWRKY proteins contain WRKYGQK/GKK domains and a C2H2/C2HC-type zinc finger structure. To identify the function of RgWRKY members involved in acteoside biosynthesis, transcriptional profiles of 37 RgWRKYs in hairy roots under salicylic acid (SA), methyl jasmonate (MeJA), and hydrogen peroxide (H2O2) treatments were systematically established using RNA-seq analysis. Based on the correlationship between the expression levels of RgWRKY genes and acteoside content, RgWRKY7, RgWRKY23, RgWRKY34, RgWRKY35, and RgWRKY37 were suggested to be involved in acteoside biosynthesis in R. glutinosa, and RgWRKY37 was selected for gene functional research. Overexpression of RgWRKY37 increased the content of acteoside and total phenylethanoid glycosides (PhGs) in hairy roots and enhanced the transcript abundance of seven enzyme genes involved in the acteoside biosynthesis pathway. These results strongly suggest the involvement of the WRKY transcription factor in the regulation of acteoside biosynthesis.

Response of bitter and sweet Chenopodium quinoa varieties to cucumber mosaic virus: Transcriptome and small RNASeq perspective

Saponins are secondary metabolites with antiviral properties. Low saponin (sweet) varieties of quinoa (Chenopodium quinoa) have been developed because seeds high in saponins taste bitter. The aim of this study was to elucidate the role of saponin in resistance of quinoa to Cucumber mosaic virus (CMV). Differential gene expression was studied in time-series study of CMV infection. High-throughput transcriptome sequence data were obtained from 36 samples (3 varieties × +/- CMV × 1 or 4 days after inoculation × 3 replicates). Translation, lipid, nitrogen, amino acid metabolism, and mono- and sesquiterpenoid biosynthesis genes were upregulated in CMV infections. In ‘Red Head’ (bitter), CMV-induced systemic symptoms were concurrent with downregulation of a key saponin biosynthesis gene, TSARL1, four days after inoculation. In local lesion responses (sweet and semi-sweet), TSARL1 levels remained up-regulated. Known microRNAs (miRNA) (81) from 11 miR families and 876 predicted novel miRNAs were identified. Differentially expressed miRNA and short interfering RNA clusters (24nt) induced by CMV infection are predicted to target genomic and intergenic regions enriched in repetitive elements. This is the first report of integrated RNASeq and sRNASeq data in quinoa-virus interactions and provides comprehensive understanding of involved genes, non-coding regions, and biological pathways in virus resistance.

Integration of Transcriptome and Small RNA Sequencing to Decipher Molecular Interaction of Chenopodium quinoa Varieties with Cucumber Mosaic Virus

Saponins are secondary metabolites with antiviral properties. Low saponin (sweet) varieties of quinoa (Chenopodium quinoa) have been developed because seeds high in saponins taste bitter. The aim of this study was to elucidate the role of saponin in resistance of quinoa to Cucumber mosaic virus (CMV). Differential gene expression was studied in time-series study of CMV infection. High-throughput transcriptome sequence data were obtained from 36 samples (3 varieties × +/- CMV × 1 or 4 days after inoculation × 3 replicates). Translation, lipid, nitrogen, amino acid metabolism, and mono- and sesquiterpenoid biosynthesis genes were upregulated in CMV infections. In ‘Red Head’, (bitter), CMV-induced systemic symptoms were concurrent with downregulation of a key saponin biosynthesis gene, TSARL1, four days after inoculation. In local lesion responses (sweet and semi-sweet), TSARL1 levels remained up-regulated. Known microRNAs (miRNA) (81) from 11 miR families and 876 predicted novel miRNAs were identified. Differentially expressed miRNA and short interfering RNA clusters (24nt) induced by CMV infection are predicted to target genomic and intergenic regions enriched in repetitive elements. This is the first report of integrated RNASeq and sRNASeq data in quinoa-virus interactions and provides comprehensive understanding of involved genes, non-coding regions, and biological pathways in virus resistance.

Gene copy number is associated with phytochemistry in Cannabis sativa

Gene copy number (CN) variation is known to be important in nearly every species where it has been examined. Alterations in gene CN may provide a fast way of acquiring diversity, allowing rapid adaptation under strong selective pressures, and may also be a key component of standing genetic variation within species. Cannabis sativa plants produce a distinguishing set of secondary metabolites, the cannabinoids, many of which have medicinal utility. Two major cannabinoids—THCA (delta-9-tetrahydrocannabinolic acid) and CBDA (cannabidiolic acid)—are products of a three-step biochemical pathway. Using whole-genome shotgun sequence data for 69 Cannabis cultivars from diverse lineages within the species, we found that genes encoding the synthases in this pathway vary in CN. Transcriptome sequence data show that the cannabinoid paralogs are differentially expressed among lineages within the species. We also found that CN partially explains variation in cannabinoid content levels among Cannabis plants. Our results demonstrate that biosynthetic genes found at multiple points in the pathway could be useful for breeding purposes, and suggest that natural and artificial selection have shaped CN variation. Truncations in specific paralogs are associated with lack of production of particular cannabinoids, showing how phytochemical diversity can evolve through a complex combination of processes.

Discovery of novel representatives of bilaterian neuropeptide families and reconstruction of neuropeptide precursor evolution in ophiuroid echinoderms

Neuropeptides are a diverse class of intercellular signalling molecules that mediate neuronal regulation of many physiological and behavioural processes. Recent advances in genome/transcriptome sequencing are enabling identification of neuropeptide precursor proteins in species from a growing variety of animal taxa, providing new insights into the evolution of neuropeptide signalling. Here, detailed analysis of transcriptome sequence data from three brittle star species, Ophionotus victoriae , Amphiura filiformis and Ophiopsila aranea , has enabled the first comprehensive identification of neuropeptide precursors in the class Ophiuroidea of the phylum Echinodermata. Representatives of over 30 bilaterian neuropeptide precursor families were identified, some of which occur as paralogues. Furthermore, homologues of endothelin/CCHamide, eclosion hormone, neuropeptide-F/Y and nucleobinin/nesfatin were discovered here in a deuterostome/echinoderm for the first time. The majority of ophiuroid neuropeptide precursors contain a single copy of a neuropeptide, but several precursors comprise multiple copies of identical or non-identical, but structurally related, neuropeptides. Here, we performed an unprecedented investigation of the evolution of neuropeptide copy number over a period of approximately 270 Myr by analysing sequence data from over 50 ophiuroid species, with reference to a robust phylogeny. Our analysis indicates that the composition of neuropeptide ‘cocktails’ is functionally important, but with plasticity over long evolutionary time scales.

Bioinformatics analysis for L polymerase homology of arenaviruses

Many pathogenic viruses can transmit between humans and animals as zoonotic infectious often to cause dangerous disease with obvious clinical signs, sometimes globally. However, these outbreaks are generally only dealt with seriously when the viruses concerned have been confirmed and classified as zoonotic. Arenavirus haemorrhagic fever viruses can cause epidemics in some areas of the world. To assess the potential of similar viruses to become zoonotic in the future, bioinformatics of relatively conserved virus proteins, such as the L polymerase of arenaviruses, can be used to identify homology with viruses that might give future outbreaks. In this paper, new and archival metazoan transcriptome sequence data was used as a resource to identify similar sequences to known arenavirus sequences for the purpose of risk prediction. In essence, bioinformatics was utilized to provide a better understanding of the potential evolution and natural history of uncharacterized virus sequences in the GenBank database. Several matches were identified which, along with a reasoned approach to their phyla, was used to provide a likelihood score of their zoonotic potential.