Rugang Yu

Genome-wide identification and characterization of cadmium-responsive microRNAs and their target genes in radish (Raphanus sativus L.) roots

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that play vital regulatory roles in plant growth, development, and environmental stress responses. Cadmium (Cd) is a non-essential heavy metal that is highly toxic to living organisms. To date, a number of conserved and non-conserved miRNAs have been identified to be involved in response to Cd stress in some plant species. However, the miRNA-mediated gene regulatory networks responsive to Cd stress in radish (Raphanus sativus L.) remain largely unexplored. To dissect Cd-responsive miRNAs and their targets systematically at the global level, two small RNA libraries were constructed from Cd-treated and Cd-free roots of radish seedlings. Using Solexa sequencing technology, 93 conserved and 16 non-conserved miRNAs (representing 26 miRNA families) and 28 novel miRNAs (representing 22 miRNA families) were identified. In all, 15 known and eight novel miRNA families were significantly differently regulated under Cd stress. The expression patterns of a set of Cd-responsive miRNAs were validated by quantitative real-time PCR. Based on the radish mRNA transcriptome, 18 and 71 targets for novel and known miRNA families, respectively, were identified by the degradome sequencing approach. Furthermore, a few target transcripts including phytochelatin synthase 1 (PCS1), iron transporter protein, and ABC transporter protein were involved in plant response to Cd stress. This study represents the first transcriptome-based analysis of miRNAs and their targets responsive to Cd stress in radish roots. These findings could provide valuable information for functional characterization of miRNAs and their targets in regulatory networks responsive to Cd stress in radish.

Identification of novel and salt-responsive miRNAs to explore miRNA-mediated regulatory network of salt stress response in radish (Raphanus sativus L.)

BackgroundSalt stress is one of the most representative abiotic stresses that severely affect plant growth and development. MicroRNAs (miRNAs) are well known for their significant involvement in plant responses to abiotic stresses. Although miRNAs implicated in salt stress response have been widely reported in numerous plant species, their regulatory roles in the adaptive response to salt stress in radish (Raphanus sativus L.), an important root vegetable crop worldwide, remain largely unknown.ResultsSolexa sequencing of two sRNA libraries from NaCl-free (CK) and NaCl-treated (Na200) radish roots were performed for systematical identification of salt-responsive miRNAs and their expression profiling in radish. Totally, 136 known miRNAs (representing 43 miRNA families) and 68 potential novel miRNAs (belonging to 51 miRNA families) were identified. Of these miRNAs, 49 known and 22 novel miRNAs were differentially expressed under salt stress. Target prediction and annotation indicated that these miRNAs exerted a role by regulating specific stress-responsive genes, such as squamosa promoter binding-like proteins (SPLs), auxin response factors (ARFs), nuclear transcription factor Y (NF-Y) and superoxide dismutase [Cu-Zn] (CSD1). Further functional analysis suggested that these target genes were mainly implicated in signal perception and transduction, regulation of ion homeostasis, basic metabolic processes, secondary stress responses, as well as modulation of attenuated plant growth and development under salt stress. Additionally, the expression patterns of ten miRNAs and five corresponding target genes were validated by reverse-transcription quantitative PCR (RT-qPCR).ConclusionsWith the sRNA sequencing, salt-responsive miRNAs and their target genes in radish were comprehensively identified. The results provide novel insight into complex miRNA-mediated regulatory network of salt stress response in radish, and facilitate further dissection of molecular mechanism underlying plant adaptive response to salt stress in root vegetable crops.

De novo sequencing of root transcriptome reveals complex cadmium-responsive regulatory networks in radish(Raphanus sativus L.)

Cadmium (Cd) is a nonessential metallic trace element that poses potential chronic toxicity to living organisms. To date, little is known about the Cd-responsive regulatory network in root vegetable crops including radish. In this study, 31,015 unigenes representing 66,552 assembled unique transcripts were isolated from radish root under Cd stress based on de novo transcriptome assembly. In all, 1496 differentially expressed genes (DEGs) consisted of 3579 transcripts were identified from Cd-free (CK) and Cd-treated (Cd200) libraries. Gene Ontology and pathway enrichment analysis indicated that the up- and down-regulated DEGs were predominately involved in glucosinolate biosynthesis as well as cysteine and methionine-related pathways, respectively. RT-qPCR showed that the expression profiles of DEGs were in consistent with results from RNA-Seq analysis. Several candidate genes encoding phytochelatin synthase (PCS), metallothioneins (MTs), glutathione (GSH), zinc iron permease (ZIPs) and ABC transporter were responsible for Cd uptake, accumulation, translocation and detoxification in radish. The schematic model of DEGs and microRNAs-involved in Cd-responsive regulatory network was proposed. This study represents a first comprehensive transcriptome-based characterization of Cd-responsive DEGs in radish. These results could provide fundamental insight into complex Cd-responsive regulatory networks and facilitate further genetic manipulation of Cd accumulation in root vegetable crops.

Transcriptome-based gene profiling provides novel insights into the characteristics of radish root response to Cr stress with next-generation sequencing

Radish (Raphanus sativus L.) is an important worldwide root vegetable crop with high nutrient values and is adversely affected by non-essential heavy metals including chromium (Cr). Little is known about the molecular mechanism underlying Cr stress response in radish. In this study, RNA-Seq technique was employed to identify differentially expressed genes (DEGs) under Cr stress. Based on de novo transcriptome assembly, there were 30,676 unigenes representing 60,881 transcripts isolated from radish root under Cr stress. Differential gene analysis revealed that 2985 uingenes were significantly differentially expressed between Cr-free (CK) and Cr-treated (Cr600) libraries, among which 1424 were up-regulated and 1561 down-regulated. Gene ontology (GO) analysis revealed that these DEGs were mainly involved in primary metabolic process, response to abiotic stimulus, cellular metabolic process and small molecule metabolic process. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that the DEGs were mainly involved in protein processing in endoplasmic reticulum, starch and sucrose metabolism, amino acid metabolism, glutathione metabolism, drug and xenobiotics by cytochrome P450 metabolism. RT-qPCR analysis showed that the expression patterns of 12 randomly selected DEGs were highly accordant with the results from RNA-seq. Furthermore, many candidate genes including signaling protein kinases, transcription factors and metal transporters, chelate compound biosynthesis and antioxidant system, were involved in defense and detoxification mechanisms of Cr stress response regulatory networks. These results would provide novel insight into molecular mechanism underlying plant responsiveness to Cr stress and facilitate further genetic manipulation on Cr uptake and accumulation in radish.

Genome-Wide Identification of Embryogenesis-Associated microRNAs in Radish (Raphanus sativus L.) by High-Throughput Sequencing

AbstractmicroRNAs (miRNAs), endogenous non-coding RNAs of approximately 21–24 nucleotides, are important regulators of transcriptional and post-transcriptional gene expression. These regulators play a key role during plant growth and development, including embryogenesis, which is crucial to the life cycle of most plant species. However, although embryogenesis-associated miRNAs have been isolated in a few species, the diversity of these regulatory miRNAs remains largely unexplored, especially in radish. In this study, two small RNA libraries were constructed from radish ovules before and after fertilization. Both libraries were sequenced by next generation sequencing (NGS) technology. This analysis identified 144 conserved and 34 non-conserved miRNAs (representing 43 known miRNA families) and 38 novel miRNAs (representing 28 miRNA families). Comparative analysis revealed that 29 known and 10 novel miRNA families were differentially expressed during embryogenesis. QRT-PCR analysis confirmed miRNA expression patterns and revealed tissue-specific and/or developmental stage-dependent expression for some miRNAs. Moreover, potential target predictions indicated that most of these targets were transcription factors involved in regulating plant growth, development and metabolism. Notably, target transcripts such as squamosa promoter-binding protein, auxin response factor and agamous-like MADS-box protein participated in radish embryogenesis.

Transcriptome-Wide Characterization of Novel and Heat-Stress-Responsive microRNAs in Radish (Raphanus Sativus L.) Using Next-Generation Sequencing

AbstractmicroRNAs (miRNAs) are a class of single-stranded endogenous non-coding RNAs that play critical roles in plant growth, development, and environmental stress responses. Temperature is one of the major physical parameters disturbing cellular homeostasis and causing leaf etiolation in plants. Previous studies have reported that several conserved and novel miRNAs were responsive to heat stress in plants. However, the characterization of miRNAs responsive to heat stress in radish remains poorly understood. To better understand miRNAs and their target genes under heat stress, two small RNA libraries were constructed from heat-treated (Heat24) and heat-untreated (CK) radish roots. Using Solexa system, totally, 26 known and 19 novel miRNAs were identified as differentially expressed under heat stress. Expression patterns of a set of heat-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). Furthermore, 422 sliced targets for 25 known miRNAs were identified by degradome sequencing technology, and most of the identified targets are involved in multiple biological processes including transcriptional regulation and response to biotic and abiotic stresses. Moreover, some miRNAs and their corresponding targets, which are related to the accumulation of heat stress transcription factors and heat shock proteins, played important roles in thermo-tolerance in radish. These findings could enhance the understanding of molecular mechanisms underlying miRNAs and their targets in regulating plant responses to heat stress.

De novo Taproot Transcriptome Sequencing and Analysis of Major Genes Involved in Sucrose Metabolism in Radish (Raphanus sativus L.).

Radish (Raphanus sativus L.) is an important annual or biennial root vegetable crop. The fleshy taproot comprises the main edible portion of the plant with high nutrition and medical value. Molecular biology study of radish begun rather later, and lacks sufficient transcriptomic and genomic data in pubic databases for understanding of the molecular mechanism during the radish taproot formation. To develop a comprehensive overview of the ‘NAU-YH’ root transcriptome, a cDNA library, prepared from three equally mixed RNA of taproots at different developmental stages including pre-cortex splitting stage, cortex splitting stage, and expanding stage was sequenced using high-throughput Illumina RNA sequencing. From approximately 51 million clean reads, a total of 70,168 unigenes with a total length of 50.28 Mb, an average length of 717 bp and a N50 of 994 bp were obtained. In total, 63,991 (about 91.20% of the assembled unigenes) unigenes were successfully annotated in five public databases including NR, GO, COG, KEGG, and Nt. GO analysis revealed that the majority of these unigenes were predominately involved in basic physiological and metabolic processes, catalytic, binding, and cellular process. In addition, a total of 103 unigenes encoding eight enzymes involved in the sucrose metabolism related pathways were also identified by KEGG pathway analysis. Sucrose synthase (29 unigenes), invertase (17 unigenes), sucrose-phosphate synthase (16 unigenes), fructokinase (17 unigenes), and hexokinase (11 unigenes) ranked top five in these eight key enzymes. From which, two genes (RsSuSy1, RsSPS1) were validated by T-A cloning and sequenced, while the expression of six unigenes were profiled with RT-qPCR analysis. These results would be served as an important public reference platform to identify the related key genes during taproot thickening and facilitate the dissection of molecular mechanisms underlying taproot formation in radish.

Transcriptome-wide analysis of chromium-stress responsive microRNAs to explore miRNA-mediated regulatory networks in radish (Raphanus sativus L.)

MicroRNAs (miRNAs) are small noncoding RNAs that play pivotal roles in plant growth, development and stress response. Chromium (Cr) is one of common environmental contaminants possessing potential health hazards to living organisms. To date, little is known about the regulatory roles of miRNAs in response to Cr stress in radish. To systematically identify Cr-responsive miRNAs and their targets in radish, two sRNA libraries derived from Cr-free (CK) and Cr-treated (Cr200) roots were constructed. With Solexa sequencing, 81 known and 72 novel miRNAs were identified, from which 54 known and 16 novel miRNAs were significantly differentially expressed under Cr stress. Several target genes for Cr-responsive miRNAs encode different transcription factor (TF) families, including SPLs, MYBs, ERFs and bZIPs, might regulate corresponding HM-related transcriptional processes in plants. Notably, a few key responsive enzymes or proteins, including HMA, YSL1 and ABC transporter protein were involved in Cr uptake and homeostasis process. Furthermore, the expression patterns of some Cr-responsive miRNAs and their targets were validated by RT-qPCR. This study represents the first characterization of Cr-responsive miRNAs and their targets in radish. The outcomes of this study could provide novel insights into miRNA-mediated regulatory mechanisms underlying plant response to Cr stress in root vegetable crops.

Transcriptome profiling of root microRNAs reveals novel insights into taproot thickening in radish ( Raphanus sativus L.)

BackgroundRadish (Raphanus sativus L.) is an economically important root vegetable crop, and the taproot-thickening process is the most critical period for the final productivity and quality formation. MicroRNAs (miRNAs) are a family of non-coding small RNAs that play an important regulatory function in plant growth and development. However, the characterization of miRNAs and their roles in regulating radish taproot growth and thickening remain largely unexplored. A Solexa high-throughput sequencing technology was used to identify key miRNAs involved in taproot thickening in radish.ResultsThree small RNA libraries from ‘NAU-YH’ taproot collected at pre-cortex splitting stage, cortex splitting stage and expanding stage were constructed. In all, 175 known and 107 potential novel miRNAs were discovered, from which 85 known and 13 novel miRNAs were found to be significantly differentially expressed during taproot thickening. Furthermore, totally 191 target genes were identified for the differentially expressed miRNAs. These target genes were annotated as transcription factors and other functional proteins, which were involved in various biological functions including plant growth and development, metabolism, cell organization and biogenesis, signal sensing and transduction, and plant defense response. RT-qPCR analysis validated miRNA expression patterns for five miRNAs and their corresponding target genes.ConclusionsThe small RNA populations of radish taproot at different thickening stages were firstly identified by Solexa sequencing. Totally 98 differentially expressed miRNAs identified from three taproot libraries might play important regulatory roles in taproot thickening. Their targets encoding transcription factors and other functional proteins including NF-YA2, ILR1, bHLH74, XTH16, CEL41 and EXPA9 were involved in radish taproot thickening. These results could provide new insights into the regulatory roles of miRNAs during the taproot thickening and facilitate genetic improvement of taproot in radish.

Transcriptome Profiling of Taproot Reveals Complex Regulatory Networks during Taproot Thickening in Radish (Raphanus sativus L.).

Radish (Raphanus sativus L.) is one of the most important vegetable crops worldwide. Taproot thickening represents a critical developmental period that determines yield and quality in radish life cycle. To isolate differentially expressed genes (DGEs) involved in radish taproot thickening process and explore the molecular mechanism underlying taproot development, three cDNA libraries from radish taproot collected at pre-cortex splitting stage (L1), cortex splitting stage (L2), and expanding stage (L3) were constructed and sequenced by RNA-Seq technology. More than seven million clean reads were obtained from the three libraries, from which 4,717,617 (L1, 65.35%), 4,809,588 (L2, 68.24%) and 4,973,745 (L3, 69.45%) reads were matched to the radish reference genes, respectively. A total of 85,939 transcripts were generated from three libraries, from which 10,450, 12,325, and 7392 differentially expressed transcripts (DETs) were detected in L1 vs. L2, L1 vs. L3, and L2 vs. L3 comparisons, respectively. Gene Ontology and pathway analysis showed that many DEGs, including EXPA9, Cyclin, CaM, Syntaxin, MADS-box, SAUR, and CalS were involved in cell events, cell wall modification, regulation of plant hormone levels, signal transduction and metabolisms, which may relate to taproot thickening. Furthermore, the integrated analysis of mRNA-miRNA revealed that 43 miRNAs and 92 genes formed 114 miRNA-target mRNA pairs were co-expressed, and three miRNA-target regulatory networks of taproot were constructed from different libraries. Finally, the expression patterns of 16 selected genes were confirmed using RT-qPCR analysis. A hypothetical model of genetic regulatory network associated with taproot thickening in radish was put forward. The taproot formation of radish is mainly attributed to cell differentiation, division and expansion, which are regulated and promoted by certain specific signal transduction pathways and metabolism processes. These results could provide new insights into the complex molecular mechanism underlying taproot thickening and facilitate genetic improvement of taproot in radish.