High-throughput sequencing, characterization and detection of new and conserved cucumber miRNAs

Role of microRNAs and their putative mechanism in regulating potato (Solanum tuberosum L.) life cycle and response to various environmental stresses

The exponentially increasing population and the demand for food is inextricably linked. This has shifted global attention to improving crop plant traits to meet global food demands. Potato (Solanum tuberosum L.) is a major non-grain food crop that is grown all over the world. Currently, some of the major global potato research work focuses on the significance of microRNAs (miRNAs) in potato. miRNAs are a type of non-coding RNAs that regulate the gene expression of their target mRNA genes by cleavage and/or their translational inhibition. This suggests an essential role of miRNAs in a multitude of plant biological processes, including maintenance of genome integrity, plant growth, development and maturation, and initiation of responses to various stress conditions. Therefore, engineering miRNAs to generate stress-resistant varieties of potato may result in high yield and improved nutritional qualities. In this review, we discuss the potato miRNAs specifically known to play an essential role in the various stages of the potato life cycle, conferring stress-resistant characteristics, and modifying gene expression. This review highlights the significance of the miRNA machinery in plants, especially potato, encouraging further research into engineering miRNAs to boost crop yields and tolerance towards stress.

Genome-Wide Identification, Quantification, and Validation of Differentially Expressed miRNAs in Eggplant (Solanum melongena L.) Based on Their Response to Ralstonia solanacearum Infection

MicroRNAs (miRNAs), a type of short noncoding RNA molecule (21-23 nucleotides), mediate repressive gene regulation through RNA silencing at the posttranscriptional level and play an important role in the defense response to abiotic and biotic stresses. miRNAs of the plant system have been studied in model crops for their diverse regulatory role while less is known about their significance in other plants whose genome and transcriptome data are scarce in the database, including eggplant (Solanum melongena L.). In the present study, a next-generation sequencing platform was used for the sequencing of miRNA, and real-time quantitative PCR for miRNAs was used to validate the gene expression patterns of miRNAs in Solanum melongena plantlets infected with the bacterial wilt-causing pathogen Ralstonia solanacearum (R. solanacearum). Sequence analyses showed the presence of 375 miRNAs belonging to 29 conserved families. The miR414 is highly conserved miRNA across the plant system while miR5658 and miR5021 were found exclusively in Arabidopsis thaliana surprisingly, these miRNAs were found in eggplants too. The most abundant families were miR5658 and miR414. Ppt-miR414, hvu-miR444b, stu-miR8020, and sly miR5303 were upregulated in Pusa purple long (PPL) (susceptible) at 48 h postinfection, followed by a decline after 96 h postinfection. A similar trend was obtained in ath-miR414, stu-mir5303h, alymiR847-5p, far-miR1134, ath-miR5021, ath-miR5658, osa-miR2873c, lja-miR7530, stu-miR7997c, and gra-miR8741 but at very low levels after infection in the susceptible variety, indicating their negative role in the suppression of host immunity. On the other hand, osa-miR2873c was found to be slightly increased after 96 hpi from 48 hpi. Most of the miRNAs under study showed relatively lower expression in the resistant variety Arka Nidhi after infection than in the susceptible variety. These results shed light on a deeper regulatory role of miRNAs and their targets in regulation of the plant response to bacterial infection. The present experiment and their results suggested that the higher expression of miRNA leads to a decline in host mRNA and thus shows susceptibility.

Parasitic plant small RNA analyses unveil parasite-specific signatures of microRNA retention, loss, and gain

Parasitism is a successful life strategy that has evolved independently in several families of vascular plants. The genera Cuscuta and Orobanche represent examples of the two profoundly different groups of parasites: one parasitizing host shoots and the other infecting host roots. In this study, we sequenced and described the overall repertoire of small RNAs from Cuscuta campestris and Orobanche aegyptiaca. We showed that C. campestris contains a number of novel microRNAs (miRNAs) in addition to a conspicuous retention of miRNAs that are typically lacking in other Solanales, while several typically conserved miRNAs seem to have become obsolete in the parasite. One new miRNA appears to be derived from a horizontal gene transfer event. The exploratory analysis of the miRNA population (exploratory due to the absence of a full genomic sequence for reference) from the root parasitic O. aegyptiaca also revealed a loss of a number of miRNAs compared to photosynthetic species from the same order. In summary, our study shows partly similar evolutionary signatures in the RNA silencing machinery in both parasites. Our data bear proof for the dynamism of this regulatory mechanism in parasitic plants.

Exogenous Application of dsRNA for the Control of Viruses in Cucurbits

The recurrent emergence of viral diseases in intensive horticultural crops requires alternative control strategies. The topical application of double-stranded RNA (dsRNA) molecules homologous to pathogens has been proposed as a tool for virus control in plants. These dsRNAs induce the silencing mechanism, the RNA interference (RNAi), that degrades homologous dsRNAs. Cucumber green mottle mosaic virus (CGMMV) represents a serious threat to cucurbit crops. Since genetic resistance to the virus is not yet available in commercial varieties, we aimed to control this virus by RNAi. For this purpose, we obtained constructions both for expressing dsRNA in bacteria to treat cucumber plants by topical application and for agroinoculation in experiments done in the growth chamber. Besides, greenhouse tests were performed in spring and in summer when plants were challenged with the virus, and differences in several parameters were investigated, including the severity of symptoms, dry weight, total height, virus accumulation, and virus-derived small interfering RNAs (vsiRNAs). Spraying of plants with dsRNA reduced significatively CGMMV symptoms in the plants in growth chamber tests. Agroinfiltration experiments done under identical conditions were also effective in limiting the progress of CGMMV disease. In the greenhouse assay performed in spring, symptoms were significantly reduced in dsRNA-sprayed plants, and the development of the plants improved with respect to non-treated plants. Virus titers and vsiRNAs were clearly reduced in dsRNA-treated plants. The effect of protection of the dsRNA was less evident in the greenhouse assay carried out in the summer. Besides, we investigated the mobility of long (ds)RNA derived from spraying or agroinfiltrated dsRNA and found that it could be detected in local, close distal, and far distal points from the site of application. VsiRNAs were also detected in local and distal points and the differences in accumulation were compared. In parallel, we investigated the capacity of dsRNAs derived from genes of tomato leaf curl New Delhi virus (ToLCNDV), another economically important virus in cucurbits, to limit the disease in zucchini, both by agroinfiltration or by direct spraying, but found no protective effect. In view of the results, the topical application of dsRNAs is postulated as a promising strategy for CGMMV control in the cucumber.

miRNA Profiling and Its Role in Multi-Omics Regulatory Networks Connected with Somaclonal Variation in Cucumber (Cucumis sativus L.)

The role of miRNAs in connection with the phenomenon of somaclonal variation, which occurs during plant in vitro culture, remains uncertain. This study aims to investigate the possible role of miRNAs in multi-omics regulatory pathways in cucumber somaclonal lines. For this purpose, we performed sRNA sequencing (sRNA-seq) from cucumber fruit samples identified 8, 10 and 44 miRNAs that are differentially expressed between somaclones (S1, S2, S3 lines) and the reference B10 line of Cucumis sativus. For miRNA identification, we use ShortStack software designed to filter miRNAs from sRNAs according to specific program criteria. The identification of predicted in-silico targets revealed 2,886 mRNAs encoded by 644 genes. The functional annotation of miRNA's target genes and gene ontology classification revealed their association with metabolic processes, response to stress, multicellular organism development, biosynthetic process and catalytic activity. We checked with bioinformatic analyses for possible interactions at the level of target proteins, differentially expressed genes (DEGs) and genes affected by genomic polymorphisms. We assume that miRNAs can indirectly influence molecular networks and play a role in many different regulatory pathways, leading to somaclonal variation. This regulation is supposed to occur through the process of the target gene cleavage or translation inhibition, which in turn affects the proteome, as we have shown in the example of molecular networks. This is a new approach combining levels from DNA-seq through mRNA-seq, sRNA-seq and in silico PPI in the area of plants' somaclonal variation.

Genome-wide miRNA analysis and integrated network for flavonoid biosynthesis in Osmanthus fragrans

Background

Osmanthus fragrans is an important economical plant containing multiple secondary metabolites including flavonoids and anthocyanins. During the past years, the roles of miRNAs in regulating the biosynthesis of secondary metabolites in plants have been widely investigated. However, few studies on miRNA expression profiles and the potential roles in regulating flavonoid biosynthesis have been reported in O. fragrans.

Results

In this study, we used high-throughput sequencing technology to analyze the expression profiles of miRNAs in leaf and flower tissues of O. fragrans. As a result, 106 conserved miRNAs distributed in 47 families and 88 novel miRNAs were identified. Further analysis showed there were 133 miRNAs differentially expressed in leaves and flowers. Additionally, the potential target genes of miRNAs as well as the related metabolic pathways were predicted. In the end, flavonoid content was measured in flower and leaf tissues and potential role of miR858 in regulating flavonoid synthesis was illustrated in O. fragrans.

Conclusions

This study not only provided the genome-wide miRNA profiles in the flower and leaf tissue of O. fragrans, but also investigated the potential regulatory role of miR858a in flavonoid synthesis in O. fragrans. The results specifically indicated the connection of miRNAs to the regulation of secondary metabolite biosynthesis in non-model economical plant.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07439-y.

Assessment of ZnO-NPs toxicity in maize: An integrative microRNAomic approach

Rapid expansion of nanotechnology and indiscriminate discharge of metal oxide nanoparticles (NPs) into the environment pose a serious hazard to the ecological receptors including plants. To better understand the role of miRNAs in ZnO-NPs stress adaptation, two small RNA libraries were prepared from control and ZnO-NPs (800 ppm, <50 nm particle size) stressed maize leaves. Meager performance of ZnO-NPs treated seedlings was associated with elevated tissue zinc accumulation, enhanced ROS generation, loss of root cell viability, increased foliar MDA content, decrease in chlorophyll and carotenoids contents. Deep sequencing identified 3 (2 known and 1 novel) up- and 77 (73 known and 4 novel) down-regulated miRNAs from ZnO-NPs challenged leaves. GO analysis reveals that potential targets of ZnO-NPs responsive miRNAs regulate diverse biological processes viz. plant growth and development (miR159f-3p, zma_18), ROS homeostasis (miR156b, miR166l), heavy metal transport and detoxification (miR444a, miR167c-3p), photosynthesis (miR171b) etc. Up-regulation of SCARECROW 6 in ZnO-NPs treated leaves might be responsible for suppression of chlorophyll biosynthesis leading to yellowing of leaves. miR156b.1 mediated up-regulation of CALLOSE SYNTHASE also does not give much protection against ZnO-NPs treatment. Taken together, the findings shed light on the miRNA-guided stress regulatory networks involved in plant adaptive responses to ZnO-NPs stress.

Deep-sequencing of Solanum commersonii small RNA libraries reveals riboregulators involved in cold stress response

Among wild species used in potato breeding, Solanum commersonii displays the highest tolerance to low temperatures under both acclimated (ACC) and non-acclimated (NACC) conditions. It is also the first wild potato relative with a known whole genome sequence. Recent studies have shown that abiotic stresses induce changes in the expression of many small non-coding RNA (sncRNA). We determined the small non-coding RNA (sncRNAome) of two clones of S. commersonii contrasting in their cold response phenotypes via smRNAseq. Differential analysis provided evidence that expression of several miRNAs changed in response to cold stress conditions. Conserved miR408a and miR408b changed their expression under NACC conditions, whereas miR156 and miR169 were differentially expressed only under ACC conditions. We also report changes in tasiRNA and secondary siRNA expression under both stress conditions. Our results reveal possible roles of sncRNA in the regulatory networks associated with tolerance to low temperatures and provide useful information for a more strategic use of genomic resources in potato breeding.

Small RNA sequencing identifies cucumber miRNA roles in waterlogging-triggered adventitious root primordia formation

The formation of adventitious roots (ARs) is a key morphological adaptation of cucumber (Cucumis sativus L.) to waterlogging stress. MicroRNAs (miRNAs) constitute a group of non-coding small RNAs (sRNA) that play crucial roles in regulating diverse biological processes, including waterlogging acclimation. However, which specific miRNAs and how they are involved in waterlogging-triggered de novo AR primordia formation are not fully known. Here, Illumina sRNA sequencing was applied to sequence six sRNA libraries generated from the waterlogging-tolerant cucumber Zaoer-N after 48 h of waterlogging and the control. A total of 358 cucumber miRNAs, 312 known and 46 novel, were obtained. Among them, 23 were differentially expressed, with 10 and 13 being up- and downregulated, respectively. A qPCR expression study confirmed that the identified differentially expressed miRNAs were credible. A total of 657 putative miRNA target genes were predicted for the 23 miRNAs using an in silico approach. A gene ontology enrichment analysis revealed that target genes functioning in cell redox homeostasis, cytoskeleton, photosynthesis and cell growth were over-represented. In total, 58 of the 657 target genes showed inverse expression patterns compared with their respective miRNAs through a combined analysis of sRNA- and RNA-sequencing-based transcriptome datasets using the same experimental design. The target gene annotation included a peroxidase, a GDSL esterases/lipase and two heavy metal-associated isoprenylated plant proteins. Our results provide an important framework for understanding the unique miRNA patterns seen in responses to waterlogging and the miRNA-mediated formation of de novo AR primordia in cucumber.

Functional Identification of Corynespora cassiicola-Responsive miRNAs and Their Targets in Cucumber

Target leaf spot (TLS), which is caused by Corynespora cassiicola (C. cassiicola), is one of the most important diseases in cucumber (Cucumis sativus L.). Our previous research identified several C. cassiicola-responsive miRNAs in cucumber by high-throughput sequencing, including two known miRNAs and two novel miRNAs. The target genes of these miRNAs were related to secondary metabolism. In this study, we verified the interaction between these miRNAs and target genes by histochemical staining and fluorescence quantitative assays of GUS. We transiently expressed the candidate miRNAs and target genes in cucumber cotyledons to investigate the resistance to C. cassiicola. Transient expression of miR164d, miR396b, Novel-miR1, and Novel-miR7 in cucumber resulted in decreased resistance to C. cassiicola, while transient expression of NAC (inhibited by miR164d), APE (inhibited by miR396b), 4CL (inhibited by Novel-miR1), and PAL (inhibited by Novel-miR7) led to enhanced resistance to C. cassiicola. In addition, overexpression of 4CL and PAL downregulated lignin synthesis, and overexpression of Novel-miR1 and Novel-miR7 also downregulated lignin synthesis, indicating that the regulation of 4CL and PAL by Novel-miR1 and Novel-miR7 could affect lignin content. The tobacco rattle virus (TRV) induced short tandem target mimic (STTM)-miRNA silencing vector was successfully constructed, and target miRNAs were successfully silenced. The identification of disease resistance and lignin content showed that silencing candidate miRNAs could improve cucumber resistance to C. cassiicola.

High-throughput sequencing reveals microRNAs in response to heat stress in the head kidney of rainbow trout (Oncorhynchus mykiss)

Recently, the research of animal microRNAs (miRNAs) has attracted wide attention for its regulatory effect in the development process and the response to abiotic stresses. Rainbow trout is a commercially and cold water fish species, and usually encounters heat stress, which affects its growth and leads to a huge economic loss. But there were few investigations about the roles of miRNAs in heat stress in rainbow trout. In this study, miRNAs of rainbow trout which were involved in heat stress were identified by high-throughput sequencing of six small RNA libraries from head kidney tissues under control (18 °C) and heat-treated (24 °C) conditions. A total of 392 conserved miRNAs and 989 novel miRNAs were identified, of which 78 miRNAs were expressed in different response to heat stress. Ten of these miRNAs were further validated by quantitative real-time PCR. In addition to, including 393 negative correlation miRNA-target gene pairs, several important regulatory pathways were involved in heat stress of the potential target genes, including protein processing in endoplasmic reticulum, NOD-like receptor signaling pathway, and phagosome. Our data significantly advance understanding of heat stress regulatory mechanism of miRNA in the head kidney of rainbow trout, which provide a useful resource for the cultivation of rainbow trout.

Identification of cucumber circular RNAs responsive to salt stress

BACKGROUND

Circular RNAs (circRNAs) are 3'-5' head-to-tail covalently closed non-coding RNA that have been proved to play essential roles in many cellular and developmental processes. However, no information relate to cucumber circRNAs is available currently, especially under salt stress condition.

RESULTS

In this study, we sequenced circRNAs in cucumber and a total of 2787 were identified, with 1934 in root and 44 in leaf being differentially regulated under salt stress. Characteristics analysis of these circRNAs revealed following features: most of them are exon circRNAs (79.51%) and they prefer to arise from middle exon(s) of parent genes (2035/2516); moreover, most of circularization events (88.3%) use non-canonical-GT/AG splicing signals; last but not least, pairing-driven circularization is not the major way to generate cucumber circRNAs since very few circRNAs (18) contain sufficient flanking complementary sequences. Annotation and enrichment analysis of both parental genes and target mRNAs were launched to uncover the functions of differentially expressed circRNAs induced by salt stress. The results showed that circRNAs may be paly roles in salt stress response by mediating transcription, signal transcription, cell cycle, metabolism adaptation, and ion homeostasis related pathways. Moreover, circRNAs may function to regulate proline metabolisms through regulating associated biosynthesis and degradation genes.

CONCLUSIONS

The present study identified large number of cucumber circRNAs and function annotation revealed their possible biological roles in response to salt stress. Our findings will lay a solid foundation for further structure and function studies of cucumber circRNAs.

Signatures of B-cell receptor diversity in B lymphocytes following Epstein-Barr virus transformation

Epstein-Barr virus (EBV) is a widespread human virus that establishes latent infection, potentially leading to tumors, hematological disorders, and other severe diseases. EBV infections are associated with diverse symptoms and affect various organs; therefore, early diagnosis and treatment are crucial. B cell receptor (BCR) repertoires of B cell surface immunoglobulins have been widely studied for their association with various infectious diseases. However, the specific genetic changes that modulate the BCR repertoires after an EBV infection are still poorly understood. In this study, we employed high-throughput sequencing (HTS) to investigate the diversity of BCR repertoires in an EBV-transformed lymphoblastic cell line (LCL). Compared with the noninfected control B cell line, the LCL exhibited a decrease in overall BCR diversity but displayed an increase in the expansion of some dominant rearrangements such as IGHV4-31/IGHJ4, IGHV4-59/IGHJ4, IGHV5-51/IGHJ3, and IGHV3-74/IGHJ3. A higher frequency of occurrence of these rearrangement types was confirmed in patients with EBV infection. Interestingly, the IGHV3-74 rearrangement was only detected in EBV-infected children, suggesting that our experimental observations were not coincidental. In addition, we identified a highly dominant consensus motif, CAR(xRx)YGSG(xYx)FD, in complementarity-determining region 3 (CDR3) sequences of the heavy chain in the LCL. Our findings demonstrated the utility of HTS technology for studying the variations in signature motifs of the BCR repertoires after EBV infection. We propose that the analysis of BCR repertoire sequences represents a promising method for diagnosing early EBV infections and developing novel antibody- and vaccine-based therapies against such infections.

Expression profiling and regulatory network of cucumber microRNAs and their putative target genes in response to cucumber green mottle mosaic virus infection

Cucumber green mottle mosaic virus (CGMMV) is an important pathogen of cucumber (Cucumis sativus). The molecular mechanisms mediating host-pathogen interactions are likely to be strongly influenced by microRNAs (miRNAs), which are known to regulate gene expression during the disease cycle. This study focused on 14 miRNAs (miR159, miR169, miR172, miR838, miR854, miR5658, csa-miRn1-3p, csa-miRn2-3p, csa-miRn3-3p, csa-miRn4-5p, csa-miRn5-5p, csa-miRn6-3p, csa-miRn7-5p and csa-miRn8-3p) and their target genes. The data collected was used to construct a regulatory network of miRNAs and target genes associated with cucumber-CGMMV interactions, which identified 608 potential target genes associated with all of the miRNAs except csa-miRn7-5p. Five of the miRNAs (miR159, miR838, miR854, miR5658 and csa-miRn6-3p) were found to be mutually linked by target genes, while another eight (miR169, miR172, csa-miRn1-3p, csa-miRn2-3p, csa-miRn3-3p, csa-miRn4-5p, csa-miRn5-5p and csa-miRn8-3p) formed subnetworks that did not display any connectivity with other miRNAs or their target genes. Reverse transcription quantitative real-time PCR (RT-qPCR) was used to analyze the expression levels of the different miRNAs and their putative target genes in leaf, stem and root samples of cucumber over a 42-day period after inoculation with CGMMV. A positive correlation was found between some of the miRNAs and their respective target genes, although for most, the response varied greatly depending on the time point, indicating that additional factors are likely to be involved in the interaction between cucumber miRNAs and their target genes. Several miRNAs, including miR159 and csa-miRn6-3p, were linked to target genes that have been associated with plant responses to disease. A model linking miRNAs, their targets and downstream biological processes is proposed to indicate the roles of these miRNAs in the cucumber-CGMMV pathosystem.

Identification of novel miRNAs and their target genes from Populus szechuanica infected with Melampsora larici-populina

Two novel miRNAs were selected from a pre-constructed RNA library of Populus szechuanica infected with the foliar rust fungus Melampsora larici-populina in order to detect the genes regulated as targets of the miRNAs novel_mir_11 and novel_mir_357. The novel miRNAs were identified from P. szechuanica using stem-loop methods and their precursors were able to fold into a complete stem loop structure. The predicted target genes of the novel miRNAs were verified with RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). The full-length sequences of target genes, RPM1 and RPS2/5, in P. szechuanica were obtained through rapid amplification of cDNA ends (RACE) and officially named PsRPM1 and PsRPS2/5. These genes contain nucleotide binding site-leucine-rich repeats (NBS-LRR) domains typical of resistance genes. The expression levels of miRNAs and their target genes in different periods post infection were analysed with quantitative real-time PCR (qRT-PCR). After infection with the foliar rust fungus, the expression levels of the novel miRNAs and their target genes were dynamic. Both novel_mir_11 and novel_mir_357 negatively regulated the expression of their target genes. In this study, the regulatory effects of two novel miRNAs through their target genes were characterized to provide further mechanistic information regarding the interaction between Populus and a foliar rust fungus. Results of this study improve our understanding of the defence response mechanisms of Populus and will stimulate future work to characterize strategies to prevent and control Populus diseases.

MicroRNAs and their targets in cucumber shoot apices in response to temperature and photoperiod

BACKGROUND

The cucumber is one of the most important vegetables worldwide and is used as a research model for study of phloem transport, sex determination and temperature-photoperiod physiology. The shoot apex is the most important plant tissue in which the cell fate and organ meristems have been determined. In this study, a series of whole-genome small RNA, degradome and transcriptome analyses were performed on cucumber shoot apical tissues treated with high vs. low temperature and long vs. short photoperiod.

RESULTS

A total of 164 known miRNAs derived from 68 families and 203 novel miRNAs from 182 families were identified. Their 4611 targets were predicted using psRobot and TargetFinder, amongst which 349 were validated by degradome sequencing. Fourteen targets of six miRNAs were differentially expressed between the treatments. A total of eight known and 16 novel miRNAs were affected by temperature and photoperiod. Functional annotations revealed that "Plant hormone signal transduction" pathway was significantly over-represented in the miRNA targets. The miR156/157/SBP-Boxes and novel-mir153/ethylene-responsive transcription factor/senescence-related protein/aminotransferase/acyl-CoA thioesterase are the two most credible miRNA/targets combinations modulating the plant's responsive processes to the temperature-photoperiod changes. Moreover, the newly evolved, cucumber-specific novel miRNA (novel-mir153) was found to target 2087 mRNAs by prediction and has 232 targets proven by degradome analysis, accounting for 45.26-58.88% of the total miRNA targets in this plant. This is the largest sum of genes targeted by a single miRNA to the best of our knowledge.

CONCLUSIONS

These results contribute to a better understanding of the miRNAs mediating plant adaptation to combinations of temperature and photoperiod and sheds light on the recent evolution of new miRNAs in cucumber.

V. AM, Solanke AU, Singh NK, Sharma TR, Gaikwad K. Genome-wide discovery of tissue-specific miRNAs in clusterbean (Cyamopsis tetragonoloba) indicates their association with galactomannan biosynthesis

Owing to the presence of 80% soluble dietary fibre, high protein content and high value gum, clusterbean (Cyamopsis tetragonoloba) has recently emerged as an economically important legume. The developing clusterbean seeds accumulate 90% galactomannans in the endosperm and, therefore, can be used as a model crop to understand galactomannan biosynthesis and its regulation. miRNAs are tiny master regulators of their corresponding target genes, resulting in variations in the amounts of their metabolic end products. To understand the role of these regulators in galactomannan biosynthesis regulation, small RNA libraries were prepared and sequenced from five tissues of clusterbean genotype RGC-936, and miRanalyzer and DSAP programs were used to identify conserved miRNAs and novel small RNAs. A total of 187 known and 171 novel miRNAs were found to be differentially expressed, of which 10 miRNAs were validated. A complicated network topology and 35% sharing of the target mRNAs between known and novel miRNAs suggest random evolution of novel miRNAs. The gene ontology (GO) annotation of potential target genes revealed the genes coding for signalling and carbohydrate metabolism (50.10%), kinases and other enzymes (20.75%), transcription factors (10.20%), transporters (8.35%) and other targets (10.6%). Two novel unigenes were annotated as ManS (mannosyltransferase/mannan synthase) and UGE (UDP- D-glucose 4-epimerase) and validated as targets for three novel miRNAs, that is Ct-miR3130, Ct-miR3135 and Ct-miR3157. Our findings reveal that these novel miRNAs could play an important role in the regulation of the galactomannan pathway in C. tetragonoloba and possibly other galactomannan-producing species.

Transcriptome and miRNA analyses of the response to Corynespora cassiicola in cucumber

Cucumber (Cucumis sativus L.) target leaf spot (TLS), which is caused by the fungus Corynespora cassiicola (C. cassiicola), seriously endangers the production of cucumber. In this assay, we performed comprehensive sequencing of the transcriptome and microRNAs (miRNAs) of a resistant cucumber (Jinyou 38) during C. cassiicola inoculation using the Illumina NextSeq 500 platform. The possible genes related to the response to C. cassiicola were associated with plant hormones, transcription factors, primary metabolism, Ca²⁺ signaling pathways, secondary metabolism and defense genes. In total, 150 target genes of these differentially expressed miRNAs were predicted by the bioinformatic analysis. By analyzing the function of the target genes, several candidate miRNAs that may be related to the response to C. cassiicola stress were selected. We also predicted 7 novel miRNAs and predicted their target genes. Moreover, the expression patterns of the candidate genes and miRNAs were tested by quantitative real-time RT-PCR. According to the analysis, genes and miRNAs associated with secondary metabolism, particularly the phenylpropanoid biosynthesis pathway, may play a major role in the resistance to C. cassiicola stress in cucumber. These results offer a foundation for future studies exploring the mechanism and key genes of resistance to cucumber TLS.

Identification of cold stress responsive microRNAs in two winter turnip rape (Brassica rapa L.) by high throughput sequencing

BACKGROUND

Low temperature is a major abiotic stress affecting the production of rapeseed in China by impeding plant growth and development. A comprehensive knowledge of small-RNA expression pattern in Brassica rapa under cold stress could improve our knowledge of microRNA-mediated stress responses.

RESULTS

A total of 353 cold-responsive miRNAs, 84 putative novel and 269 conserved miRNAs, were identified from the leaves and roots of two winter turnip rape varieties 'Longyou 7' (cold-tolerant) and 'Tianyou 4' (cold-sensitive), which were stressed under - 4 °C for 8 h. Eight conserved (miR166h-3p-1, miR398b-3p, miR398b-3p-1, miR408d, miR156a-5p, miR396h, miR845a-1, miR166u) and two novel miRNAs (Bra-novel-miR3153-5p and Bra-novel-miR3172-5p) were differentially expressed in leaves of 'Longyou 7' under cold stress. Bra-novel-miR3936-5p was up-regulated in roots of 'Longyou 7' under cold stress. Four and five conserved miRNAs were differentially expressed in leaves and roots of 'Tianyou 4' after cold stress. Besides, we found two conserved miRNAs (miR319e and miR166m-2) were down-regulated in non-stressed roots of 'Longyou 7' compared with 'Tianyou 4'. After cold stress, we found two and eight miRNAs were differentially expressed in leaves and roots of 'Longyou 7' compared with 'Tianyou 4'. The differentially expressed miRNAs between two cultivars under cold stress include novel miRNAs and the members of the miR166 and miR319 families. A total of 211 target genes for 15 known miRNAs and two novel miRNAs were predicted by bioinformatic analysis, mainly involved in metabolic processes and stress responses. Five differentially expressed miRNAs and predicted target genes were confirmed by quantitative reverse transcription PCR, and the expressional changes of target genes were negatively correlated to differentially expressed miRNAs. Our data indicated that some candidate miRNAs (e.g., miR166e, miR319, and Bra-novel-miR3936-5p) may play important roles in plant response to cold stress.

CONCLUSIONS

Our work indicates that miRNA and putative target genes mediated metabolic processes and stress responses are significant to cold tolerance in B. rapa.

Identifying optimal reference genes for the normalization of microRNA expression in cucumber under viral stress

Cucumber green mottle mosaic virus (CGMMV) is an economically important pathogen and causes significant reduction of both yield and quality of cucumber (Cucumis sativus). Currently, there were no satisfied strategies for controlling the disease. A better understanding of microRNA (miRNA) expression related to the regulation of plant-virus interactions and virus resistance would be of great assistance when developing control strategies for CGMMV. However, accurate expression analysis is highly dependent on robust and reliable reference gene used as an internal control for normalization of miRNA expression. Most commonly used reference genes involved in CGMMV-infected cucumber are not universally expressed depending on tissue types and stages of plant development. It is therefore crucial to identify suitable reference genes in investigating the role of miRNA expression. In this study, seven reference genes, including Actin, Tubulin, EF-1α, 18S rRNA, Ubiquitin, GAPDH and Cyclophilin, were evaluated for the most accurate results in analyses using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Gene expression was assayed on cucumber leaves, stems and roots that were collected at different days post inoculation with CGMMV. The expression data were analyzed using algorithms including delta-Ct, geNorm, NormFinder, and BestKeeper as well as the comparative tool RefFinder. The reference genes were subsequently validated using miR159. The results showed that EF-1α and GAPDH were the most reliable reference genes for normalizing miRNA expression in leaf, root and stem samples, while Ubiquitin and EF-1α were the most suitable combination overall.

Identification and characterization of microRNAs in tree peony during chilling induced dormancy release by high-throughput sequencing

Tree peony, one of the most valuable horticultural and medicinal plants in the world, has to go through winter to break dormancy. Growing studies from molecular aspects on dormancy release process have been reported, but inadequate study has been done on miRNA-guided regulation in tree peony. In this study, high-throughput sequencing was employed to identify and characterize miRNAs in three libraries (6 d, 18 d and 24 d chilling treatments). There were 7,122, 10,076 and 9,097 unique miRNA sequences belonging to 52, 87 and 68 miRNA families, respectively. A total of 32 conserved miRNAs and 17 putative novel miRNAs were identified during dormancy release. There were 771 unigenes as potential targets of 62 miRNA families. Total 112 known miRNAs were differentially expressed, of which 55 miRNAs were shared among three libraries and 28 miRNAs were only found in 18 d chilling duration library. The expression patterns of 15 conserved miRNAs were validated and classified into four types by RT-qPCR. Combining with our microarray data under same treatments, five miRNAs (miR156k, miR159a, miR167a, miR169a and miR172a) were inversely correlated to those of their target genes. Our results would provide new molecular basis about dormancy release in tree peony.

Identification and comparative analysis of microRNAs from tomato varieties showing contrasting response to ToLCV infections

Increasing incidence of viral infections in crop plants adversely affects their growth and yield. Tomato (Solanum lycopersicum) is considered to be a favorite host for viruses with over 50 species of begomoviruses naturally infecting this crop. Tomato leaf curl virus (ToLCV) is among the most widespread and devastating begomoviruses affecting tomato production. microRNAs (miRs) have been established as key regulators of gene expression and plant development. The miR pathways are disturbed during infection by viruses. Thus, comprehension of regulatory miR networks is crucial in understanding the effect of viral pathogenicity. To identify key miRs involved in ToLCV infection, a high throughput approach involving next generation sequencing was employed. Healthy and infected leaf tissues of two tomato varieties, differing in their susceptibility to ToLCV infection were analyzed. NGS data analysis followed by computational predictions, led to identification of 91 known miRs, 15 novel homologs and 53 novel miRs covering two different varieties of tomato, susceptible (Pusa Ruby) and tolerant (LA1777) to ToLCV infection. The cleaved targets of these miRs were identified using online available degradome libraries from leaf, flower and fruit of tomato and showed their involvement in various biological pathways through KEGG Orthology. With detailed comparative profiling of expression pattern of these miRs, we could associate the specific miRs with the resistant and infected genotypes. This study depicted that in depth analysis of miR expression patterns and their functions will help in identification of molecules that can be used for manipulation of gene expression to increase crop production and developing resistance against diseases.

Evolution of microRNA827 targeting in the plant kingdom

Unlike most ancient microRNAs, which conservatively target homologous genes across species, microRNA827 (miR827) targets two different types of SPX (SYG1/PHO81/XPR1)-domain-containing genes, NITROGEN LIMITATION ADAPTATION (NLA) and PHOSPHATE TRANSPORTER 5 (PHT5), in Arabidopsis thaliana and Oryza sativa to regulate phosphate (Pi) transport and storage, respectively. However, how miR827 shifted its target preference and its evolutionary history are unknown. Based on target prediction analysis, we found that in most angiosperms, miR827 conservatively targets PHT5 homologs, but in Brassicaceae and Cleomaceae it preferentially targets NLA homologs, and we provide evidence for the transition of target preference during Brassicales evolution. Intriguingly, we found a lineage-specific loss of the miR827-regulatory module in legumes. Analysis of miR827-mediated cleavage efficiency and the expression of PHT5 in A. thaliana indicated that accumulation of mutations in the target site and the exclusion of the target site by alternative transcriptional initiation eliminated PHT5 targeting by miR827. Here, we identified a transition of miR827 target preference during plant evolution and revealed the uniqueness of miR827-mediated regulation among conserved plant miRNAs. Despite the change in its target preference, upregulation of miR827 by Pi starvation and its role in regulating cellular Pi homeostasis were retained.

Comparative analysis of miRNA and mRNA abundance in determinate cucumber by high-throughput sequencing

Determinate cucumber is a type of short vines, fewer nodes, and terminal flowers, it is suitable for high-density planting and available harvesting in field cultivation, whereas the indeterminate cucumber is preferred to cultivate in greenhouses. However, many biotic or abiotic stresses could lead indeterminate cucumber to be determinate in greenhouse cultivation, which may decrease yield and fruit quality. Therefore, it is urgent and essential to investigate the key factors forming determinate and terminal flowering in cucumber. In this study, two close background inbred lines were selected and conducted the miRNA and mRNA high throughput sequencing. Interestingly, ethylene-associated miRNAs and mRNAs were intensively obtained, indicating that the plant hormone ethylene is a key factor impacting determinate and terminal flowering in cucumber. The ethylene metabolites analysis showed that significant higher ethylene was observed in determinate line than that in the indeterminate line. The RT-qPCR validation of ethylene related miRNAs Cas-miR172, Cas-miR396, and Cas-miR414 and their target mRNAs showed a significant negative correlation. These data suggested that ethylene-associated miRNAs might affect determinate and terminal flower phenotypes by regulating their target genes expression. This study not only provides a potential molecular mechanism for determinate formation in cucumber but also establishes a method to demonstrate important physiological processes through the comprehensive association of miRNA and mRNA high-throughput sequencing.

Mapping the microRNA Expression Profiles in Glyoxalase Over-expressing Salinity Tolerant Rice

In the recent years, glyoxalase pathway has been an active area of research in both human and plants. This pathway is reported to confer stress tolerance in plants, by modulating the glutathione homeostasis to achieve detoxification of a potent cytotoxic and mutagenic compound, methylglyoxal. The microRNAs (miRNAs) are also reported to play significant role in stress tolerance for plants. However, the cross-talk of miRNAs with the metabolism regulated by glyoxalase in the salinity-tolerance is unexplored. We therefore investigated whether expression profiles of miRNAs are altered in response to glyoxalase overexpression, and if any of these are also responsible for modulating the stress responses of plants. In this study, the Next Generation Sequencing (NGS) was employed to profile miRNA expression levels from glyoxalase overexpressing transgenic lines. The associated targets of differentially expressed miRNAs were predicted and their functional annotation was carried out using Gene Ontology (GO) and KEGG Orthology (KO), which showed their involvement in several crucial biological pathways. The analysis of NGS datasets also identified other isoforms or isomiRs of selected miRNAs, which may have an active role in developing tolerance against salt stress. Different aspects of miRNA modifications were also studied in glyoxalase overexpressing lines.

Uncovering Male Fertility Transition Responsive miRNA in a Wheat Photo-Thermosensitive Genic Male Sterile Line by Deep Sequencing and Degradome Analysis

MicroRNAs (miRNAs) are endogenous small RNAs which play important negative regulatory roles at both the transcriptional and post-transcriptional levels in plants. Wheat is the most commonly cultivated plant species worldwide. In this study, RNA-seq analysis was used to examine the expression profiles of miRNA in the spikelets of photo-thermosenisitive genic male sterile (PTGMS) wheat line BS366 during male fertility transition. Through mapping on their corresponding precursors, 917-7,762 novel miRNAs were found in six libraries. Six novel miRNAs were selected for examination of their secondary structures and confirmation by stem-loop RT-PCR. In a differential expression analysis, 20, 22, and 58 known miRNAs exhibited significant differential expression between developmental stages 1 (secondary sporogenous cells had formed), 2 (all cells layers were present and mitosis had ceased), and 3 (meiotic division stage), respectively, of fertile and sterile plants. Some of these differential expressed miRNAs, such as tae-miR156, tae-miR164, tae-miR171, and tae-miR172, were shown to be associated with their targets. These targets were previously reported to be related to pollen development and/or male sterility, indicating that these miRNAs and their targets may be involved in the regulation of male fertility transition in the PTGMS wheat line BS366. Furthermore, target genes of miRNA cleavage sites were validated by degradome sequencing. In this study, a possible signal model for the miRNA-mediated signaling pathway during the process of male fertility transition in the PTGMS wheat line BS366 was developed. This study provides a new perspective for understanding the roles of miRNAs in male fertility in PTGMS lines of wheat.

Nano-based delivery of RNAi in cancer therapy

Background

RNA interference (RNAi), a newly developed method in which RNA molecules inhibit gene expression, has recently received considerable research attention. In the development of RNAi-based therapies, nanoparticles, which have distinctive size effects along with facile modification strategies and are capable of mediating effective RNAi with targeting potential, are attracting extensive interest.

Objective

This review presents an overview of the mechanisms of RNAi molecules in gene therapy and the different nanoparticles used to deliver RNAi molecules; briefly describes the current uses of RNAi in cancer therapy along with the nano-based delivery of RNA molecules in previous studies; and highlights some other carriers that have been applied in clinical settings. Finally, we discuss the nano-based delivery of RNAi therapeutics in preclinical development, including the current status and limitations of anti-cancer treatment.

Conclusion

With the growing number of RNAi therapeutics entering the clinical phase, various nanocarriers are expected to play important roles in the delivery of RNAi molecules for cancer therapeutics.

Genome-wide identification of Hami melon miRNAs with putative roles during fruit development

MicroRNAs represent a family of small endogenous, non-coding RNAs that play critical regulatory roles in plant growth, development, and environmental stress responses. Hami melon is famous for its attractive flavor and excellent nutritional value, however, the mechanisms underlying the fruit development and ripening remains largely unknown. Here, we performed small RNA sequencing to investigate the roles of miRNAs during Hami melon fruit development. Two batches of flesh samples were collected at four fruit development stages. Small RNA sequencing yielded a total of 54,553,424 raw reads from eight libraries. 113 conserved miRNAs belonging to 30 miRNA families and nine novel miRNAs comprising nine miRNA families were identified. The expression of 42 conserved miRNAs and three Hami melon-specific miRNAs significantly changed during fruit development. Furthermore, 484 and 124 melon genes were predicted as putative targets of 29 conserved and nine Hami melon-specific miRNA families, respectively. GO enrichment analysis were performed on target genes, "transcription, DNA-dependent", "rRNA processing", "oxidation reduction", "signal transduction", "regulation of transcription, DNA-dependent", and "metabolic process" were the over-represented biological process terms. Cleavage sites of six target genes were validated using 5' RACE. Our results present a comprehensive set of identification and characterization of Hami melon fruit miRNAs and their potential targets, which provide valuable basis towards understanding the regulatory mechanisms in programmed process of normal Hami fruit development and ripening. Specific miRNAs could be selected for further research and applications in breeding practices.

Genome-wide analysis of microRNA targeting impacted by SNPs in cucumber genome

Background

microRNAs (miRNAs) are endogenous small RNAs that play important regulatory functions in plant development. Genetic variations in miRNAs sequences or their target-binding sites (microRNA-target interaction sites) can alter miRNA targets in animal and human. Whether these single nucleotide polymorphisms (SNPs) in plant are functional have not yet been determined.

Results

In this study, we constructed leaf, root, and stem-derived small libraries of cucumber (Cucumis sativus) line 9930 (cultivated China-group cucumber) and C. sativus var. hardwickii (wild India group cucumber). A total of 22 conserved miRNA families, nine less-conserved miRNA families, and 49 cucumber-specific miRNAs were identified in both line 9930 and hardwickii. We employed cucumber resequencing data to perform a genome-wide scan for SNPs in cucumber miRNA-target interaction sites, including miRNA mature sequences and miRNA-target binding sites. As a result, we identified a total of 19 SNPs in mature miRNA sequences and 113 SNPs in miRNA-target binding sites with the potential to affect miRNA-target interactions. Furthermore, we experimentally confirmed that these SNPs produced 14 9930-unique targets mRNAs and 15 hardwickii-unique targets mRNA for cucumber miRNAs. This is the first experimental validation of SNPs in miRNA-target interaction sites affecting miRNA-target binding in plants.

Conclusions

Our results indicate that SNPs can alter miRNA function and produce unique miRNA targets in cultivated and wild cucumbers. Therefore, miRNA-related SNPs may have played important in events that led to the agronomic differences between domestic and wild cucumber.

Electronic supplementary material

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

Deciphering microRNAs and Their Associated Hairpin Precursors in a Non-Model Plant, Abelmoschus esculentus

MicroRNAs (miRNAs) are crucial regulatory RNAs, originated from hairpin precursors. For the past decade, researchers have been focusing extensively on miRNA profiles in various plants. However, there have been few studies on the global profiling of precursor miRNAs (pre-miRNAs), even in model plants. Here, for the first time in a non-model plant—Abelmoschus esculentus with negligible genome information—we are reporting the global profiling to characterize the miRNAs and their associated pre-miRNAs by applying a next generation sequencing approach. Preliminarily, we performed small RNA (sRNA) sequencing with five biological replicates of leaf samples to attain 207,285,863 reads; data analysis using miRPlant revealed 128 known and 845 novel miRNA candidates. With the objective of seizing their associated hairpin precursors, we accomplished pre-miRNA sequencing to attain 83,269,844 reads. The paired end reads are merged and adaptor trimmed, and the resulting 40–241 nt (nucleotide) sequences were picked out for analysis by using perl scripts from the miRGrep tool and an in-house built shell script for Minimum Fold Energy Index (MFEI) calculation. Applying the stringent criteria of the Dicer cleavage pattern and the perfect stem loop structure, precursors for 57 known miRNAs of 15 families and 18 novel miRNAs were revealed. Quantitative Real Time (qRT) PCR was performed to determine the expression of selected miRNAs.

The miRNAome of Catharanthus roseus: identification, expression analysis, and potential roles of microRNAs in regulation of terpenoid indole alkaloid biosynthesis

MicroRNAs (miRNAs) regulate numerous crucial biological processes in plants. However, information is limited on their involvement in the biosynthesis of specialized metabolites in plants, including Catharanthus roseus that produces a number of pharmaceutically valuable, bioactive terpenoid indole alkaloids (TIAs). Using small RNA-sequencing, we identified 181 conserved and 173 novel miRNAs (cro-miRNAs) in C. roseus seedlings. Genome-wide expression analysis revealed that a set of cro-miRNAs are differentially regulated in response to methyl jasmonate (MeJA). In silico target prediction identified 519 potential cro-miRNA targets that include several auxin response factors (ARFs). The presence of cleaved transcripts of miRNA-targeted ARFs in C. roseus cells was confirmed by Poly(A) Polymerase-Mediated Rapid Amplification of cDNA Ends (PPM-RACE). We showed that auxin (indole acetic acid, IAA) repressed the expression of key TIA pathway genes in C. roseus seedlings. Moreover, we demonstrated that a miRNA-regulated ARF, CrARF16, binds to the promoters of key TIA pathway genes and repress their expression. The C. roseus miRNAome reported here provides a comprehensive account of the cro-miRNA populations, as well as their abundance and expression profiles in response to MeJA. In addition, our findings underscore the importance of miRNAs in posttranscriptional control of the biosynthesis of specialized metabolites.

Genome-Wide Identification and Characterization of MicroRNAs and Target Genes in Lonicera japonica

MiRNAs function in post-transcriptional regulation of gene expression and play very important roles in plant development. Lonicera japonica is one of the important medicinal plants in China. However, few studies on the discovery of conserved and novel miRNAs from L. japonica were reported. In this study, we employed deep sequencing technology to identify miRNAs in leaf and flower tissues of L. japonica. A total of 22.97 million clean reads from flower and leaf tissues were obtained, which generated 146 conserved miRNAs distributed in 20 families and 110 novel miRNAs. Accordingly, 72 differentially expressed miRNAs (P≤0.001) between leaves and flowers and their potential target genes were identified and validated. The qRT-PCR validation showed that majority of the differentially expressed miRNAs showed significant tissue-specific expression in L. japonica. Furthermore, the miRNA-mRNA and mRNA-mRNA regulatory networks were constructed using Cytoscape software. Taken together, this study identified a large number of miRNAs and target genes in L. japonica, which not only provides the first global miRNA expression profiles, but also sheds light on functional genomics research on L. japonica in the future.

Characterization of miRNAs responsive to exogenous ethylene in grapevine berries at whole genome level

MicroRNAs (miRNAs) are critical regulators of various biological and metabolic processes of plants. Numerous miRNAs and their functions have been identified and analyzed in many plants. However, till now, the involvement of miRNAs in the response of grapevine berries to ethylene has not been reported yet. Here, Solexa technology was employed to deeply sequence small RNA libraries constructed from grapevine berries treated with and without ethylene. A total of 124 known and 78 novel miRNAs were identified. Among these miRNAs, 162 miRNAs were clearly responsive to ethylene, with 55 downregulated, 59 upregulated, and 14 unchanged miRNAs detected only in the control. The other 35 miRNAs responsive to ethylene were induced by ethylene and detected only in the ethylene-treated grapevine materials. Expression analysis of 27 conserved and 26 novel miRNAs revealed that 13 conserved and 18 novel ones were regulated by ethylene during the whole development of grapevine berries. High-throughput sequencing and qRT-PCR assays revealed consistent results on the expression results of ethylene-responsive miRNAs. Moreover, 90 target genes for 34 novel miRNAs were predicted, most of which were involved in responses to various stresses, especially like exogenous ethylene treatment. The identified miRNAs may be mainly involved in grapevine berry development and response to various environmental conditions.

Identification of conserved and novel microRNAs in Porphyridium purpureum via deep sequencing and bioinformatics

Background

Porphyridium purpureum has been utilized in important industrial and pharmaceutical fields. The identification of microRNAs (miRNAs) in this unique species is of great importance: such identification can help fill gaps in the small RNA (sRNA) studies of this organism and help to elucidate essential biological processes and their regulation mechanisms in this special micro alga.

Results

In this study, 254 high-confidence miRNAs (203 conserved miRNAs and 51 novel miRNAs) were identified by sRNA deep sequencing (sRNA-seq) combined with bioinformatics. A total of 235 putative miRNA families were predicted, including 192 conserved families and 43 species-specific families. The conservation and diversity of predicted miRNA families were analysed in different plant species. Both the 100 % northern blot validation rate (VR) of four randomly selected miRNAs and the results of stem-loop quantitative real time RT-PCR (qRT-PCR) assays of 25 randomly selected miRNAs demonstrated that the majority of the miRNAs identified in this study are credible. A total of 14,958 and 2184 genes were predicted to be targeted by the 186 conserved and 41 novel miRNAs. Gene ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that some target genes likely provide valuable references for further understanding of vital functions in P. purpureum. In addition, a cytoscape network will provide some clues for research into the complex biological processes that occur in this unique alga.

Conclusions

We first identified a large set of conserved and novel miRNAs in P. purpureum. The characteristic and validation analysis on miRNAs demonstrated authenticity of identification data. Functional annotation of target genes and metabolic pathways they involved in illuminated the direction for further utilization and development this micro alga based on its unique properties.

Electronic supplementary material

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

Identification of Ice Plant (Mesembryanthemum crystallinum L.) MicroRNAs Using RNA-Seq and Their Putative Roles in High Salinity Responses in Seedlings

The halophyte Mesembryanthemum crystallinum (common or crystalline ice plant) is a useful model for studying molecular mechanisms of salt tolerance. The morphology, physiology, metabolism, and gene expression of ice plant have been studied and large-scale analyses of gene expression profiling have drawn an outline of salt tolerance in ice plant. A rapid root growth to a sudden increase in salinity was observed in ice plant seedlings. Using a fluorescent dye to detect Na(+), we found that ice plant roots respond to an increased flux of Na(+) by either secreting or storing Na(+) in specialized cells. High-throughput sequencing was used to identify small RNA profiles in 3-day-old seedlings treated with or without 200 mM NaCl. In total, 135 conserved miRNAs belonging to 21 families were found. The hairpin precursor of 19 conserved mcr-miRNAs and 12 novel mcr-miRNAs were identified. After 6 h of salt stress, the expression of most mcr-miRNAs showed decreased relative abundance, whereas the expression of their corresponding target genes showed increased mRNA relative abundance. The cognate target genes are involved in a broad range of biological processes: transcription factors that regulate growth and development, enzymes that catalyze miRNA biogenesis for the most conserved mcr-miRNA, and proteins that are involved in ion homeostasis and drought-stress responses for some novel mcr-miRNAs. Analyses of the functions of target genes revealed that cellular processes, including growth and development, metabolism, and ion transport activity are likely to be enhanced in roots under salt stress. The expression of eleven conserved miRNAs and two novel miRNAs were correlated reciprocally with predicted targets within hours after salt stress exposure. Several conserved miRNAs have been known to regulate root elongation, root apical meristem activity, and lateral root formation. Based upon the expression pattern of miRNA and target genes in combination with the observation of Na(+) distribution, ice plant likely responds to increased salinity by using Na(+) as an osmoticum for cell expansion and guard cell opening. Excessive Na(+) could either be secreted through the root epidermis or stored in specialized leaf epidermal cells. These responses are regulated in part at the miRNA-mediated post-transcriptional level.

Microarray and genetic analysis reveals that csa-miR159b plays a critical role in abscisic acid-mediated heat tolerance in grafted cucumber plants

Root-shoot communication plays a vital role in plant growth, development and adaptation to environmental stimuli. Grafting-induced stress tolerance is associated with the induction of plentiful stress-related genes and proteins; the mechanism involved, however, remains obscure. Here, we show that the enhanced tolerance against heat stress in cucumber plants with luffa as rootstock was accompanied with an increased accumulation of abscisic acid (ABA), down-regulation of a subset of microRNAs (miRNAs) but up-regulation of their target genes and CsHSP70 accumulation in the shoots. Significantly, luffa rootstock and foliar application of ABA both down-regulated csa-miR159b and up-regulated its target mRNAs CsGAMYB1 and CsMYB29-like and CsHSP70 accumulation in cucumber, while ectopic expression of csa-miR159b led to decreased heat tolerance, AtMYB33 transcript and AtHSP70 accumulation in Arabidopsis plants. Taken together, our results suggest that root-originated signals such as ABA could alter miRNAs in the shoots, which have a major role in the post-transcriptional regulation of the stress-responsive genes.

Genome-wide identification of microRNAs in pomegranate (Punica granatum L.) by high-throughput sequencing

BACKGROUND

MicroRNAs (miRNAs), a class of small non-coding endogenous RNAs that regulate gene expression post-transcriptionally, play multiple key roles in plant growth and development and in biotic and abiotic stress response. Knowledge and roles of miRNAs in pomegranate fruit development have not been explored.

RESULTS

Pomegranate, which accumulates a large amount of anthocyanins in skin and arils, is valuable to human health, mainly because of its antioxidant properties. In this study, we developed a small RNA library from pooled RNA samples from young seedlings to mature fruits and identified both conserved and pomegranate-specific miRNA from 29,948,480 high-quality reads. For the pool of 15- to 30-nt small RNAs, ~50 % were 24 nt. The miR157 family was the most abundant, followed by miR156, miR166, and miR168, with variants within each family. The base bias at the first position from the 5' end had a strong preference for U for most 18- to 26-nt sRNAs but a preference for A for 18-nt sRNAs. In addition, for all 24-nt sRNAs, the nucleotide U was preferred (97 %) in the first position. Stem-loop RT-qPCR was used to validate the expression of the predominant miRNAs and novel miRNAs in leaves, male and female flowers, and multiple fruit developmental stages; miR156, miR156a, miR159a, miR159b, and miR319b were upregulated during the later stages of fruit development. Higher expression of miR156 in later fruit developmental may positively regulate anthocyanin biosynthesis by reducing SPL transcription factor. Novel miRNAs showed variation in expression among different tissues. These novel miRNAs targeted different transcription factors and hormone related regulators. Gene ontology and KEGG pathway analyses revealed predominant metabolic processes and catalytic activities, important for fruit development. In addition, KEGG pathway analyses revealed the involvement of miRNAs in ascorbate and linolenic acid, starch and sucrose metabolism; RNA transport; plant hormone signaling pathways; and circadian clock.

CONCLUSION

Our first and preliminary report of miRNAs will provide information on the synthesis of biochemical compounds of pomegranate for future research. The functions of the targets of the novel miRNAs need further investigation.

Identification and Characterization of miRNAs in Chondrus crispus by High-Throughput Sequencing and Bioinformatics Analysis

Chondrus crispus, an economically and medicinally important red alga, is a medicinally active substance and important for anti-tumor research. In this study, 117 C. crispus miRNAs (108 conserved and 9 novel) were identified from 2,416,181 small-RNA reads using high-throughput sequencing and bioinformatics methods. According to the BLAST search against the miRBase database, these miRNAs belonged to 110 miRNA families. Sequence alignment combined with homology searching revealed both the conservation and diversity of predicted potential miRNA families in different plant species. Four and 19 randomly selected miRNAs were validated by northern blotting and stem-loop quantitative real-time reverse transcription polymerase chain reaction detection, respectively. The validation rates (75% and 94.7%) demonstrated that most of the identified miRNAs could be credible. A total of 160 potential target genes were predicted and functionally annotated by Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes analysis. We also analyzed the interrelationship of miRNAs, miRNA-target genes and target genes in C. crispus by constructing a Cytoscape network. The 117 miRNAs identified in our study should supply large quantities of information that will be important for red algae small RNA research.

Transcriptome and Small RNAome Dynamics during a Resistant and Susceptible Interaction between Cucumber and Downy Mildew

Cucumber ( L.) downy mildew, caused by the obligate oomycete pathogen (Berk. and Curt.) Rostov., is the primary factor limiting cucumber production. Although sources of resistance have been identified, such as plant introduction line PI 197088, the genes and processes involved in mediating resistance are still unknown. In the current study, we conducted a comprehensive transcriptome and small RNAome analysis of a resistant (PI 197088) and susceptible ('Vlaspik') cucumber during a time course of infection using Illumina sequencing. We identified significantly differentially expressed (DE) genes within and between resistant and susceptible cucumber leaves over a time course of infection. Weighted gene correlation network analyses (WGCNA) created coexpression modules containing genes with unique expression patterns between Vlaspik and PI 197088. Recurring data trends indicated that resistance to cucumber downy mildew is associated with earlier response to the pathogen, hormone signaling, and regulation of nutrient supply. Candidate resistance genes were identified from multiple transcriptome analyses and literature support. Additionally, parallel sequencing of small RNAs (sRNAs) from cucumber and during the infection time course was used to identify and quantify novel and existing microRNA (miRNA) in both species. Predicted miRNA targets of cucumber transcripts suggest a complex interconnectedness of gene expression regulation in this plant-pathogen system. This work bioinformatically uncovered gene expression patterns involved in the mediation of or response to resistance. Herein, we provide the foundation for future work to validate candidate resistance genes and miRNA-based regulation proposed in this study.

Dehydration-responsive miRNAs in foxtail millet: genome-wide identification, characterization and expression profiling

A set of novel and known dehydration-responsive miRNAs have been identified in foxtail millet. These findings provide new insights into understanding the functional role of miRNAs and their respective targets in regulating plant response to dehydration stress. MicroRNAs perform significant regulatory roles in growth, development and stress response of plants. Though the miRNA-mediated gene regulatory networks under dehydration stress remain largely unexplored in plant including foxtail millet (Setaria italica), which is a natural abiotic stress tolerant crop. To find out the dehydration-responsive miRNAs at the global level, four small RNA libraries were constructed from control and dehydration stress treated seedlings of two foxtail millet cultivars showing contrasting tolerance behavior towards dehydration stress. Using Illumina sequencing technology, 55 known and 136 novel miRNAs were identified, representing 22 and 48 miRNA families, respectively. Eighteen known and 33 novel miRNAs were differentially expressed during dehydration stress. After the stress treatment, 32 dehydration-responsive miRNAs were up-regulated in tolerant cultivar and 22 miRNAs were down-regulated in sensitive cultivar, suggesting that miRNA-mediated molecular regulation might play important roles in providing contrasting characteristics to these cultivars. Predicted targets of identified miRNAs were found to encode various transcription factors and functional enzymes, indicating their involvement in broad spectrum regulatory functions and biological processes. Further, differential expression patterns of seven known miRNAs were validated by northern blot and expression of ten novel dehydration-responsive miRNAs were confirmed by SL-qRT PCR. Differential expression behavior of five miRNA-target genes was verified under dehydration stress treatment and two of them also validated by RLM RACE. Overall, the present study highlights the importance of dehydration stress-associated post-transcriptional regulation governed by miRNAs and their targets in a naturally stress-tolerant model crop.

Unveiling the Micronome of Cassava (Manihot esculenta Crantz)

MicroRNAs (miRNAs) are an important class of endogenous non-coding single-stranded small RNAs (21-24 nt in length), which serve as post-transcriptional negative regulators of gene expression in plants. Despite the economic importance of Manihot esculenta Crantz (cassava) only 153 putative cassava miRNAs (from multiple germplasm) are available to date in miRBase (Version 21), and identification of a number of miRNAs from the cassava EST database have been limited to comparisons with Arabidopsis. In this study, mature sequences of all known plant miRNAs were used as a query for homologous searches against cassava EST and GSS databases, and additional identification of novel and conserved miRNAs were gleaned from next generation sequencing (NGS) of two cassava landraces (T200 from southern Africa and TME3 from West Africa) at three different stages post explant transplantation and acclimatization. EST and GSS derived data revealed 259 and 32 miRNAs in cassava, and one of the miRNA families (miR2118) from previous studies has not been reported in cassava. NGS data collectively displayed expression of 289 conserved miRNAs in leaf tissue, of which 230 had not been reported previously. Of the 289 conserved miRNAs identified in T200 and TME3, 208 were isomiRs. Thirty-nine novel cassava-specific miRNAs of low abundance, belonging to 29 families, were identified. Thirty-eight (98.6%) of the putative new miRNAs identified by NGS have not been previously reported in cassava. Several miRNA targets were identified in T200 and TME3, highlighting differential temporal miRNA expression between the two cassava landraces. This study contributes to the expanding knowledge base of the micronome of this important crop.

Next generation sequencing and omics in cucumber (Cucumis sativus L.) breeding directed research

In the post-genomic era the availability of genomic tools and resources is leading us to novel generation methods in plant breeding, as they facilitate the study of the genotype and its relationship with the phenotype, in particular for complex traits. In this study we have mainly concentrated on the Cucumis sativus and (but much less) Cucurbitaceae family several important vegetable crops. There are many reports on research conducted in Cucurbitaceae plant breeding programs on the ripening process, phloem transport, disease resistance, cold tolerance and fruit quality traits. This paper presents the role played by new omic technologies in the creation of knowledge on the mechanisms of the formation of the breeding features. The analysis of NGS (NGS-next generation sequencing) data allows the discovery of new genes and regulatory sequences, their positions, and makes available large collections of molecular markers. Genome-wide expression studies provide breeders with an understanding of the molecular basis of complex traits. Firstly a high density map should be created for the reference genome, then each re-sequencing data could be mapped and new markers brought out into breeding populations. The paper also presents methods that could be used in the future for the creation of variability and genomic modification of the species in question. It has been shown also the state and usefulness in breeding the chloroplastomic and mitochondriomic study.

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.

Identification of MicroRNAs in Response to Different Day Lengths in Soybean Using High-Throughput Sequencing and qRT-PCR

MicroRNAs (miRNAs) are short, non-coding single-strand RNA molecules that play important roles in plant growth, development and stress responses. Flowering time affects the seed yield and quality of soybean. However, the miRNAs involved in the regulation of flowering time in soybean have not been reported until recently. Here, high-throughput sequencing and qRT-PCR were used to identify miRNAs involved in soybean photoperiodic pathways. The first trifoliate leaves of soybean that receive the signal of light treatment were used to construct six libraries (0, 8, and 16 h under short-day (SD) treatment and 0, 8, and 16 h under long-day (LD) treatment). The libraries were sequenced using Illumina Solexa. A total of 318 known plant miRNAs belonging to 163 miRNA families and 81 novel predicted miRNAs were identified. Among these, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h, showed differences in abundance between LD and SD treatments. Furthermore, the results of GO and KEGG analyses indicated that most of the miRNA targets were transcription factors. Seven miRNAs at 0 h, 23 miRNAs (including four novel predicted miRNAs) at 8 h, 16 miRNAs (including one novel predicted miRNA) at 16 h and miRNA targets were selected for qRT-PCR analysis to assess the accuracy of the sequencing and target prediction. The results indicated that the expression patterns of the selected miRNAs and miRNA targets showed no differences between the qRT-PCR and sequencing results. In addition, 23 miRNAs at 0 h, 65 miRNAs at 8 h and 83 miRNAs at 16 h responded to day length changes in soybean, including six novel predicted miRNAs at 8 h and six novel predicted miRNAs at 16 h. These results provided an important molecular basis to understand the regulation of flowering time through photoperiodic pathways in soybean.

High-Throughput Sequencing Identifies Novel and Conserved Cucumber (Cucumis sativus L.) microRNAs in Response to Cucumber Green Mottle Mosaic Virus Infection

Seedlings of Cucumis sativus L. (cv. 'Zhongnong 16') were artificially inoculated with Cucumber green mottle mosaic virus (CGMMV) at the three-true-leaf stage. Leaf and flower samples were collected at different time points post-inoculation (10, 30 and 50 d), and processed by high throughput sequencing analysis to identify candidate miRNA sequences. Bioinformatic analysis using screening criteria, and secondary structure prediction, indicated that 8 novel and 23 known miRNAs (including 15 miRNAs described for the first time in vivo) were produced by cucumber plants in response to CGMMV infection. Moreover, gene expression profiles (p-value <0.01) validated the expression of 3 of the novel miRNAs and 3 of the putative candidate miRNAs and identified a further 82 conserved miRNAs in CGMMV-infected cucumbers. Gene ontology (GO) analysis revealed that the predicted target genes of these 88 miRNAs, which were screened using the psRNATarget and miRanda algorithms, were involved in three functional categories: 2265 in molecular function, 1362 as cellular components and 276 in biological process. The subsequent Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the predicted target genes were frequently involved in metabolic processes (166 pathways) and genetic information processes (40 pathways) and to a lesser degree the biosynthesis of secondary metabolites (12 pathways). These results could provide useful clues to help elucidate host-pathogen interactions in CGMMV and cucumber, as well as for the screening of resistance genes.

High-throughput sequencing of small RNAs and anatomical characteristics associated with leaf development in celery

MicroRNAs (miRNAs) exhibit diverse and important roles in plant growth, development, and stress responses and regulate gene expression at the post-transcriptional level. Knowledge about the diversity of miRNAs and their roles in leaf development in celery remains unknown. To elucidate the roles of miRNAs in celery leaf development, we identified leaf development-related miRNAs through high-throughput sequencing. Small RNA libraries were constructed using leaves from three stages (10, 20, and 30 cm) of celery cv.'Ventura' and then subjected to high-throughput sequencing and bioinformatics analysis. At Stage 1, Stage 2, and Stage 3 of 'Ventura', a total of 333, 329, and 344 conserved miRNAs (belonging to 35, 35, and 32 families, respectively) were identified. A total of 131 miRNAs were identified as novel in 'Ventura'. Potential miRNA target genes were predicted and annotated using the eggNOG, GO, and KEGG databases to explore gene functions. The abundance of five conserved miRNAs and their corresponding potential target genes were validated. Expression profiles of novel potential miRNAs were also detected. Anatomical characteristics of the leaf blades and petioles at three leaf stages were further analyzed. This study contributes to our understanding on the functions and molecular regulatory mechanisms of miRNAs in celery leaf development.

Identification of fruit related microRNAs in cucumber (Cucumis sativus L.) using high-throughput sequencing technology

MicroRNAs (miRNAs) are approximately 21 nt noncoding RNAs that influence the phenotypes of different species through the post-transcriptional regulation of gene expression. Although many miRNAs have been identified in a few model plants, less is known about miRNAs specific to cucumber (Cucumis sativus L.). In this study, two libraries of cucumber RNA, one based on fruit samples and another based on mixed samples from leaves, stems, and roots, were prepared for deep-sequencing. A total of 110 sequences were matched to known miRNAs in 47 families, while 56 sequences in 46 families are newly identified in cucumber. Of these, 77 known and 44 new miRNAs were differentially expressed, with a fold-change of at least 2 and p-value < 0.05. In addition, we predicted the potential targets of known and new miRNAs. The identification and characterization of known and new miRNAs will enable us to better understand the role of these miRNAs in the formation of cucumber fruit.

High Throughput Sequencing of Small RNAs in the Two Cucurbita Germplasm with Different Sodium Accumulation Patterns Identifies Novel MicroRNAs Involved in Salt Stress Response

MicroRNAs (miRNAs), a class of small non-coding RNAs, recognize their mRNA targets based on perfect sequence complementarity. MiRNAs lead to broader changes in gene expression after plants are exposed to stress. High-throughput sequencing is an effective method to identify and profile small RNA populations in non-model plants under salt stresses, significantly improving our knowledge regarding miRNA functions in salt tolerance. Cucurbits are sensitive to soil salinity, and the Cucurbita genus is used as the rootstock of other cucurbits to enhance salt tolerance. Several cucurbit crops have been used for miRNA sequencing but salt stress-related miRNAs in cucurbit species have not been reported. In this study, we subjected two Cucurbita germplasm, namely, N12 (Cucurbita. maxima Duch.) and N15 (Cucurbita. moschata Duch.), with different sodium accumulation patterns, to Illumina sequencing to determine small RNA populations in root tissues after 4 h of salt treatment and control. A total of 21,548,326 and 19,394,108 reads were generated from the control and salt-treated N12 root tissues, respectively. By contrast, 19,108,240 and 20,546,052 reads were obtained from the control and salt-treated N15 root tissues, respectively. Fifty-eight conserved miRNA families and 33 novel miRNAs were identified in the two Cucurbita germplasm. Seven miRNAs (six conserved miRNAs and one novel miRNAs) were up-regulated in salt-treated N12 and N15 samples. Most target genes of differentially expressed novel miRNAs were transcription factors and salt stress-responsive proteins, including dehydration-induced protein, cation/H⁺ antiporter 18, and CBL-interacting serine/threonine-protein kinase. The differential expression of miRNAs between the two Cucurbita germplasm under salt stress conditions and their target genes demonstrated that novel miRNAs play an important role in the response of the two Cucurbita germplasm to salt stress. The present study initially explored small RNAs in the response of pumpkin to salt stress, and provided valuable information on novel miRNAs and their target genes in Cucurbita.