Transcriptome-wide analysis of dynamic variations in regulation modes of grapevine microRNAs on their target genes during grapevine development

Construction of a miRNA Signature Using Support Vector Machine to Identify Microsatellite Instability Status and Prognosis in Gastric Cancer

Background. The specific role and prognostic value of DNA repair and replication-associated miRNAs in gastric cancer (GC) have not been clearly elucidated. Therefore, comprehensive analysis of miRNAs in GC is crucial for proposing therapeutic strategies and survival prediction. Methods. Firstly, clinical information and transcriptome data of TCGA-GC were downloaded from the database. In the entire cohort, we performed differential analysis in all miRNAs and support vector machine (SVM) was used to eliminate redundant miRNAs. Subsequently, we combined survival data and cox regression analysis to construct a miRNA signature in the training cohort. In addition, we used PCA, Kaplan-Meier, and ROC analysis to explore the prognosis value of risk score in the training and testing cohort. It is worth noting that multiple algorithms were used to evaluate difference of immune microenvironment (TME), microsatellite instability (MSI), tumor mutational burden (TMB), and immunotherapy in different risk groups. Finally, we investigated the potential mechanism about miRNA signature. Results. We constructed miRNA signature based on the following 4 miRNAs: hsa-miR-139-5p, hsa-miR-139-3p, hsa-miR-146b-5p, and hsa-miR-181a-3p. Univariate and multivariate Cox regression analyses suggested that risk score is a risk factor and an independent prognostic factor in GC patients. The AUC value of ROC analysis showed a robust prediction accuracy in each cohort. Moreover, significant differences in immune functions, immune cell content, immune checkpoint, MSI status, and TMB score were excavated in different groups distinguished by risk score. Finally, based on the above four miRNA target genes, we revealed that the signature was enriched in DNA repair and replication. Conclusion. We have developed a robust risk-formula based on 4 miRNAs that provides accurate risk stratification and prognostic prediction for GC patients. In addition, different risk subgroups may potentially guide the choice of targeted therapy.

Stress responses and epigenomic instability mark the loss of somatic embryogenesis competence in grapevine

Somatic embryogenesis (SE) represents the most appropriate tool for next-generation breeding methods in woody plants such as grapevine (Vitis vinifera L.). However, in this species, the SE competence is strongly genotype-dependent and the molecular basis of this phenomenon is poorly understood. We explored the genetic and epigenetic basis of SE in grapevine by profiling the transcriptome, epigenome, and small RNAome of undifferentiated, embryogenic, and non-embryogenic callus tissues derived from two genotypes differing in competence for SE, Sangiovese and Cabernet Sauvignon. During the successful formation of embryonic callus, we observed the upregulation of epigenetic-related transcripts and short interfering RNAs in association with DNA hypermethylation at transposable elements in both varieties. Nevertheless, the switch to nonembryonic development matched the incomplete reinforcement of transposon silencing, and the evidence of such effect was more apparent in the recalcitrant Cabernet Sauvignon. Transcriptomic differences between the two genotypes were maximized already at early stage of culture where the recalcitrant variety expressed a broad panel of genes related to stress responses and secondary metabolism. Our data provide a different angle on the SE molecular dynamics that can be exploited to leverage SE as a biotechnological tool for fruit crop breeding.

Identification of miRNAs-mediated seed and stone-hardening regulatory networks and their signal pathway of GA-induced seedless berries in grapevine (V. vinifera L.)

BACKGROUND

Stone-hardening stage is crucial to the development of grape seed and berry quality. A significant body of evidence supports the important roles of MicroRNAs in grape-berry development, but their specific molecular functions during grape stone-hardening stage remain unclear.

RESULTS

Here, a total of 161 conserved and 85 species-specific miRNAs/miRNAs* (precursor) were identified in grape berries at stone-hardening stage using Solexa sequencing. Amongst them, 30 VvmiRNAs were stone-hardening stage-specific, whereas 52 exhibited differential expression profiles during berry development, potentially participating in the modulation of berry development as verified by their expression patterns. GO and KEGG pathway analysis showed that 13 VvmiRNAs might be involved in the regulation of embryo development, another 11 in lignin and cellulose biosynthesis, and also 28 in the modulation of hormone signaling, sugar, and proline metabolism. Furthermore, the target genes for 4 novel VvmiRNAs related to berry development were validated using RNA Ligase-Mediated (RLM)-RACE and Poly(A) Polymerase-Mediated (PPM)-RACE methods, and their cleavage mainly occurred at the 9th-11th sites from the 5' ends of miRNAs at their binding regions. In view of the regulatory roles of GA in seed embryo development and stone-hardening in grape, we investigated the expression modes of VvmiRNAs and their target genes during GA-induced grape seedless-berry development, and we validated that GA induced the expression of VvmiR31-3p and VvmiR8-5p to negatively regulate the expression levels of CAFFEOYL COENZYME A-3-O-METHYLTRANSFERASE (VvCCoAOMT), and DDB1-CUL4 ASSOCIATED FACTOR1 (VvDCAF1). The series of changes might repress grape stone hardening and embryo development, which might be a potential key molecular mechanism in GA-induced grape seedless-berry development. Finally, a schematic model of miRNA-mediated grape seed and stone-hardening development was proposed.

CONCLUSION

This work identified 30 stone-hardening stage-specific VvmiRNAs and 52 significant differential expression ones, and preliminary interpreted the potential molecular mechanism of GA-induced grape parthenocarpy. GA negatively manipulate the expression of VvCCoAOMT and VvDCAF1 by up-regulation the expression of VvmiR31-3p and VvmiR8-5p, thereby repressing seed stone and embryo development to produce grape seedless berries.

Gastric cancer-associated microRNA expression signatures: integrated bioinformatics analysis, validation, and clinical significance

Background

Gastric cancer (GC) is one of the common gastrointestinal malignancy worldwide and exhibits a poor prognosis. Increasing studies have indicated that microRNAs play critical roles in the cancer progression and have shown great potential as useful biomarkers. The search for potential diagnostic and prognostic biomarkers of gastric cancer (GC) with integrated bioinformatics analyses has been undertaken in previous studies.

Methods

In this study, the robust rank aggregation (RRA) method was used to perform an integrated analysis of differentially expressed miRNAs (DEMs) from five microarray datasets in the Gene Expression Omnibus (GEO) database to find robust biomarkers for GC. Ultimately, seven miRNAs were filtered from fourteen primary miRNAs using the validation set of The Cancer Genome Atlas (TCGA) database. Based on these results, diagnostic and survival analyses were performed, and logistic regression and Cox regression were used to determine the clinicopathological characteristics of the DEM expression and overall survival.

Results

Nine eligible miRNA datasets related to GC were selected from the GEO database for integrated analysis in this study. Diagnostic analysis implied that these miRNAs could be regarded as promising candidate diagnostic biomarkers in GC tissues, but whether the results of the tissue analysis are consistent with those of peripheral blood analysis requires further validation. The logistic regression indicated that the ectopic expression of these DEMs was relevant to the histological type, anatomical region, and pathological grade of GC. However, the survival and Cox regression analyses suggested that the poor prognosis of GC patients was not strongly dependent on the ectopic expression of the seven miRNAs, but rather, a poor prognosis was associated with age, metastasis, and histological grade.

Conclusions

Based on the results presented in this study it can be concluded that these miRNAs (miR-455-3p, miR-135b-5p, let-7a-3p, miR-195-5p, miR-204-5p, miR-149-5p, and miR-143-3p) might be potential biomarkers for the early diagnosis of GC patients, but this finding should be regarded with caution. A large-scale, prospective, and multicenter cohort study should be performed.

Integrative Analysis of Gene Expression and miRNAs Reveal Biological Pathways Associated with Bud Paradormancy and Endodormancy in Grapevine

Transition of grapevine buds from paradormancy to endodormancy is coordinated by changes in gene expression, phytohormones, transcription factors, and other molecular regulators, but the mechanisms involved in transcriptional and post-transcriptional regulation of dormancy stages are not well delineated. To identify potential regulatory targets, an integrative analysis of differential gene expression profiles and their inverse relationships with miRNA abundance was performed in paradormant (long day (LD) 15 h) or endodormant (short day (SD), 13 h) Vitis riparia buds. There were 400 up- and 936 downregulated differentially expressed genes in SD relative to LD budsGene set and gene ontology enrichment analysis indicated that hormone signaling and cell cycling genes were downregulated in SD relative to LD buds. miRNA abundance and inverse expression analyses of miRNA target genes indicated increased abundance of miRNAs that negatively regulate genes involved with cell cycle and meristem development in endodormant buds and miRNAs targeting starch metabolism related genes in paradormant buds. Analysis of interactions between abundant miRNAs and transcription factors identified a network with coinciding regulation of cell cycle and epigenetic regulation related genes in SD buds. This network provides evidence for cross regulation occurring between miRNA and transcription factors both upstream and downstream of MYB3R1.

Characterization of the regulation mechanism of grapevine microRNA172 family members during flower development

BACKGROUND

Grapevine (Vitis vinifera L.), which has important nutritional values and health benefits, is one of the most economically important fruit crops cultivated worldwide. Several studies showed a large number of microRNAs (VvmiRNAs) involved in the modulation of grape growth and development, and many VvmiRNA families have multiple members. However, the way by which various members from the same miRNA family work is unclear, particularly in grapes.

RESULTS

In this study, an important conserved VvmiR172 family (VvmiR172s) and their targets were set as a good example for elucidating the interaction degree, mechanism, and spatio-temporal traits of diverse members from the same miRNA family. miR-RACE and Stem-loop RT-PCR were used to identify the spatio-temporal expressions of various members of VvmiR172s; together with RLM-RACE, PPM-RACE, Western blot, transgenic technologies, their interaction degree, and regulation mechanism were further validated. The expression of VvmiR172c was significantly higher than that of VvmiR172a, b, and d and showed a positive correlation with the abundance of VvAP2 cleavage products. These findings indicated that VvmiR172c might be one of the main action factors of the VvmiR172 family in flower development. The ability of VvmiR172c to cleave target genes differed due to divergence in complementary degree with VvAP2 and expression levels of various members. In VvmiR172 transgenic lines, we observed that 35S::VvmiR172c resulted in the earliest and abundant flowering, indicating the strong function of VvmiR172c. In contrast, the non-significant phenotypic changes were detected in the VvAP2 transgenic lines. The qRT-PCR and Western bolt results demonstrated that VvmiR172c plays a major role in targeting VvAP2.

CONCLUSIONS

VvmiR172 up-regulated the expression of NtFT and decreased the expression of NtFLC. The up/down regulation of VvmiR172c was the most pronounced. The functions of four VvmiR172 members in grape differed, and miR172c had the strongest regulation on AP2.

Characterization on the conservation and diversification of miRNA156 gene family from lower to higher plant species based on phylogenetic analysis at the whole genomic level

miRNA156 family members (miR156s) participate in regulating the transition of plant vegetative and reproductive growth, flower development, and formation of berry skin color by negatively modulating their target gene SPLs. However, the evolution and functional diversification of miR156s in plants remain elusive. Phylogenetic analysis on 310 miR156s from 51 plant species on miRBase 21.0 showed that only miR156a could be conserved in the 51 plant species, but their sequences exhibited variation; another set of miR156s, such as miR156m/n/o/p/q/r/s/t/u/v/w/x/y/z, was identified only in certain special plant species (Glycine max and Malus); also, all base variations in the sequences of 310 miR156s occurred within one miR156 seed sequence, "TGACAGAAGAGAGTGAGCAC," and the changed base sites were mainly located at the 11th and 14th bases from the 5' end of the miR156 seed sequence, in which some base variations of miR156s resulted in a difference in miR156 targeting modes; by contrast, miR156 precursor sequences are highly divergent across diverse species. Similarly, cis-regulatory motifs on the promoter sequence of MIR156s in various plants also exhibited significant discrepancy. The intragenic MIR156 genes overlapped their target SBP genes, thereby suggesting that some microRNAs (miRNAs) originate from duplication of target genes. These traits might be the reasons of the conservation and diversification of miR156 gene family. This study identified the conserved seed sequence "TGACAGAAGAGAGTGAGCAC," and the sequence variation characterization, of miR156 family evolution, also investigated the varied traits of their promoters, precursors, and mature sequences in sequence evolutions and found some miRNAs might originate from duplication of target genes. Our findings will contribute to our understanding of the functional diversification of miRNAs and the interactions of miRNA/target pairs based on the evolutionary history of miRNA genes.

VvmiR160s/VvARFs interaction and their spatio-temporal expression/cleavage products during GA-induced grape parthenocarpy

BACKGROUND

Grape (Vitis vinifera) is highly sensitive to gibberellin (GA), which effectively induce grape parthenocarpy. Studies showed that miR160s and their target AUXIN RESPONSIVE FACTOR (ARF) responding hormones are indispensable for various aspects of plant growth and development, but their functions in GA-induced grape parthenocarpy remain elusive.

RESULTS

In this study, the morphological changes during flower development in response to GA treatments were examined in the 'Rosario Bianco' cultivar. The precise sequences of VvmiR160a/b/c/d/e and their VvARF10/16/17 target genes were cloned, sequenced and characterized. The phylogenetic relationship and intron-exon structure of VvARFs and other ARF family members derived from different species were investigated. All VvmiR160s (except VvmiR160b) and VvARF10/16/17 had the common cis-elements responsive to GA, which support their function in GA-mediated grape parthenocarpy. The cleavage role of VvmiR160s-mediated VvARF10/16/17 was verified in grape flowers. Moreover, spatio-temporal expression analysis demonstrated that among VvmiR160 family, VvmiR160a/b/c highly expressed at late stage of flower/berry development, while VvARF10/16/17showed a reverse expression trend. Interestingly, GA exhibited a long-term effect through inducing the expression of VvmiR160a/b/c/e to increase their cleavage product accumulations from 5 to 9 days after treatment, but GA enhanced the expressions of VvARF10/16/17 only at short term. Pearson correlation analysis based on expression data revealed a negative correlation between VvmiR160a/b/c and VvARF10/16/17 in flowers not berries during GA-induced grape parthenocarpy.

CONCLUSIONS

This work demonstrated that the negative regulation of VvARF10/16/17 expression by VvmiR160a/b/c as key regulatory factors is critical for GA-mediated grape parthenocarpy, and provide significant implications for molecular breeding of high-quality seedless berry.

De novo origination of MIRNAs through generation of short inverted repeats in target genes

MIRNA (MIR) gene origin and early evolutionary processes, such as hairpin precursor sequence origination, promoter activity acquirement and the sequence of these two processes, are fundamental and fascinating subjects. Three models, including inverted gene duplication, spontaneous evolution and transposon transposition, have been proposed for de novo origination of hairpin precursor sequence. However, these models still open to discussion. In addition, de novo origination of MIR gene promoters has not been well investigated. Here, I systematically investigated the origin of evolutionarily young polyphenol oxidase gene (PPO)-targeting MIRs, including MIR1444, MIR058 and MIR12112, and a genomic region termed AasPPO-as-hp, which contained a hairpin-forming sequence. I found that MIR058 precursors and the hairpin-forming sequence of AasPPO-as-hp originated in an ancient PPO gene through forming short inverted repeats. Palindromic-like sequences and imperfect inverted repeats in the ancient PPO gene contributed to initiate the generation of short inverted repeats probably by causing errors during DNA duplication. Analysis of MIR058 and AasPPO-as-hp promoters showed that they originated in the 3'-flanking region of the ancient PPO gene. Promoter activities were gained by insertion of a CAAT-box and multiple-copper-response element (CuRE)-containing miniature inverted-repeat transposable element (MITE) in the upstream of AT-rich TATA-box-like sequence. Gain of promoter activities occurred before hairpin-forming sequence origination. Sequence comparison of MIR1444, MIR058 and MIR12112 promoters showed frequent birth and death of CuREs, indicating copper could be vital for the origination and evolution of PPO-targeting MIRs. Based on the evidence obtained, a novel model for plant MIR origination and evolution is proposed.

The Role of UV-B light on Small RNA Activity During Grapevine Berry Development

We explored the effects of ultraviolet B radiation (UV-B) on the developmental dynamics of microRNAs and phased small-interfering-RNA (phasi-RNAs)-producing loci by sequencing small RNAs in vegetative and reproductive organs of grapevine (Vitis vinifera L.). In particular, we tested different UV-B conditions in in vitro-grown plantlets (high-fluence exposition) and in berries from field-grown (radiation filtering) and greenhouse-grown (low- and high-fluence expositions) adult plants throughout fruit development and ripening. The functional significance of the observed UV-coordinated miRNA responses was supported by degradome evidences of ARGONAUTE (AGO)-programmed slicing of mRNAs. Co-expression patterns of the up-regulated miRNAs miR156, miR482, miR530, and miR828 with cognate target gene expressions in response to high-fluence UV-B was tested by q-RT-PCR. The observed UV-response relationships were also interrogated against two published UV-stress and developmental transcriptome datasets. Together, the dynamics observed between miRNAs and targets suggest that changes in target abundance are mediated transcriptionally and, in some cases, modulated post-transcriptionally by miRNAs. Despite the major changes in target abundance are being controlled primarily by those developmental effects that are similar between treatments, we show evidence for novel miRNA-regulatory networks in grape. A model is proposed where high-fluence UV-B increases miR168 and miR530 that target ARGONAUTE 1 (AGO1) and a Plus-3 domain mRNA, respectively, while decreasing miR403 that targets AGO2, thereby coordinating post-transcriptional gene silencing activities by different AGOs. Up-regulation of miR3627/4376 could facilitate anthocyanin accumulation by antagonizing a calcium effector, whereas miR395 and miR399, induced by micronutrient deficiencies known to trigger anthocyanin accumulation, respond positively to UV-B radiation. Finally, increases in the abundance of an anthocyanin-regulatory MYB-bHLH-WD40 complex elucidated in Arabidopsis, mediated by UV-B-induced changes in miR156/miR535, could contribute to the observed up-regulation of miR828. In turn, miR828 would regulate the AtMYB113-ortologues MYBA5, A6 and A7 (and thereby anthocyanins) via a widely conserved and previously validated auto-regulatory loop involving miR828 and phasi TAS4abc RNAs.

Identification of copper (Cu) stress-responsive grapevine microRNAs and their target genes by high-throughput sequencing

MicroRNAs (miRNAs) are a class of single-stranded non-coding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play important roles in plant growth, development and stress responses. With more copper (Cu) and copper containing compounds used as bactericides and fungicides in plants, Cu stress has become one of the serious environmental problems that affect plant growth and development. In order to uncover the hidden response mechanisms of Cu stress, two small RNA libraries were constructed from Cu-treated and water-treated (Control) leaves of 'Summer Black' grapevine. Following high-throughput sequencing and filtering, a total of 158 known and 98 putative novel miRNAs were identified in the two libraries. Among these, 100 known and 47 novel miRNAs were identified as differentially expressed under Cu stress. Subsequently, the expression patterns of nine Cu-responsive miRNAs were validated by quantitative real-time PCR (qRT-PCR). There existed some consistency in expression levels of Cu-responsive miRNAs between high throughput sequencing and qRT-PCR assays. The targets prediction of miRNAs indicates that miRNA may regulate some transcription factors, including AP2, SBP, NAC, MYB and ARF during Cu stress. The target genes for two known and two novel miRNAs showed specific cleavage sites at the 10th and/or 11th nucleotide from the 5'-end of the miRNA corresponding to their miRNA complementary sequences. The findings will lay the foundation for exploring the role of the regulation of miRNAs in response to Cu stress and provide valuable gene information for breeding some Cu-tolerant grapevine cultivars.

miRVIT: A Novel miRNA Database and Its Application to Uncover Vitis Responses to Flavescence dorée Infection

Micro(mi)RNAs play crucial roles in plant developmental processes and in defense responses to biotic and abiotic stresses. In the last years, many works on small RNAs in grapevine (Vitis spp.) were published, and several conserved and putative novel grapevine-specific miRNAs were identified. In order to reorganize the high quantity of available data, we produced "miRVIT," the first database of all novel grapevine miRNA candidates characterized so far, and still not deposited in miRBase. To this aim, each miRNA accession was renamed, repositioned in the last version of the grapevine genome, and compared with all the novel and conserved miRNAs detected in grapevine. Conserved and novel miRNAs cataloged in miRVIT were then used for analyzing Vitis vinifera plants infected by Flavescence dorée (FD), one of the most severe phytoplasma diseases affecting grapevine. The analysis of small RNAs from healthy, recovered (plants showing spontaneous and stable remission of symptoms), and FD-infected "Barbera" grapevines showed that FD altered the expression profiles of several miRNAs, including those involved in cell development and photosynthesis, jasmonate signaling, and disease resistance response. The application of miRVIT in a biological context confirmed the effectiveness of the followed approach, especially for the identification of novel miRNA candidates in grapevine. miRVIT database is available at http://mirvit.ipsp.cnr.it. Highlights: The application of the newly produced database of grapevine novel miRNAs to the analysis of plants infected by Flavescence dorée reveals key roles of miRNAs in photosynthesis and jasmonate signaling.

Spatio-temporal expression of miRNA159 family members and their GAMYB target gene during the modulation of gibberellin-induced grapevine parthenocarpy

Grapevine, Vitis vinifera, is an important economic fruit crop that is highly sensitive to gibberellin (GA), and the exogenous application of GA can efficiently induce grapevine parthenocarpy. However, the molecular mechanisms underlying this process remain elusive. In this study, morphological changes during flower development in response to GA treatments were examined in the 'Zuijinxiang' cultivar. To obtain insights into the roles of miRNA159s in GA-induced grapevine parthenocarpy, VvmiR159a, VvmiR159b, VvmiR159c, and their target gene VvGAMYB were isolated, sequenced and characterized. Spatial-temporal expression analyses showed that VvmiR159c exhibited the highest expression levels at 4 d before flowering, followed by a gradual decrease, while VvGAMYB displayed an opposite pattern of expression with the lowest expression at the corresponding stage in response to GA treatment. A cleavage interaction between VvmiR159s and VvGAMYB and variations of their cleavage roles were confirmed in grapevine floral development. In addition, the potential roles of VvmiR159s in GA signaling were investigated through DELLA-protein repressors, indicating that GA-DELLA (SLR1)-VvmiR159c-VvGAMYB is the key signaling regulatory module in grapevine. Our findings provide novel insights into the GA-responsive roles of VvmiR159s in modulating grapevine floral development, which have important implications for the molecular breeding of high-quality seedless grapevine berry.

Characterization of Vv-miR156: Vv-SPL pairs involved in the modulation of grape berry development and ripening

SPL is a plant-specific transcription factor family. Many researchers reported that SPL members targeted by miR156s could play crucial roles in the modulation of plant growth and development. Although there are similar reports on grapes, till now little is known about grape berry development and ripening. To gain more insight into how grape miR156s (Vv-miR156s) modulated the above given processes of grape berries by mediating their target gene Vv-SPLs, here we identified the precise sequences of Vv-miR156s in 'Giant Rose' grape berries, predicted their potential targets, and revealed that the matching degree of various Vv-miR156: Vv-SPL pairs exhibited some discrepancy, implying the divergence of their interaction. Subsequently, we also discovered similar motifs such as ABRE, CGTCA and ERE, which are more specific to berry development and ripening, within the promoters of both Vv-MIR156s and Vv-SPLs. With berry development and ripening, meanwhile, Vv-miR156a, b/c/d, e and f/g/i exhibited an overall increasing expression trend, while their targets showed opposite trends at the corresponding stages. Additionally, exogenous ABA and NAA application promoted or curbed the expression of Vv-miR156s to some extent, before grape berry ripening stage. The cleavage products, sites and frequencies of Vv-miR156a, b/c/d, e, f/g/i and their respective targets (Vv-SPL2, 9, 10, 16) during grape berry development and ripening process were validated by our developed PPM-RACE and modified RLM-RACE together with qRT-PCR, which demonstrated that Vv-miR156s can be involved in the modulation of grape berry development and ripening process by mediating the expression of Vv-SPL2, 9, 10, 16. Our findings lay an important foundation for further recognizing their functions in grape berries, and enrich the knowledge of the regulatory mechanism of miRNA-mediated grape berry development and ripening.

Comparative Analysis of Fruit Ripening-Related miRNAs and Their Targets in Blueberry Using Small RNA and Degradome Sequencing

MicroRNAs (miRNAs) play vital roles in the regulation of fruit development and ripening. Blueberry is an important small berry fruit crop with economical and nutritional value. However, nothing is known about the miRNAs and their targets involved in blueberry fruit ripening. In this study, using high-throughput sequencing of small RNAs, 84 known miRNAs belonging to 28 families and 16 novel miRNAs were identified in white fruit (WF) and blue fruit (BF) libraries, which represent fruit ripening onset and in progress, respectively. Among them, 41 miRNAs were shown to be differentially expressed during fruit maturation, and 16 miRNAs representing 16 families were further chosen to validate the sRNA sequencing data by stem-loop qRT-PCR. Meanwhile, 178 targets were identified for 41 known and 7 novel miRNAs in WF and BF libraries using degradome sequencing, and targets of miR160 were validated using RLM-RACE (RNA Ligase-Mediated (RLM)-Rapid Amplification of cDNA Ends) approach. Moreover, the expression patterns of 6 miRNAs and their targets were examined during fruit development and ripening. Finally, integrative analysis of miRNAs and their targets revealed a complex miRNA-mRNA regulatory network involving a wide variety of biological processes. The findings will facilitate future investigations of the miRNA-mediated mechanisms that regulate fruit development and ripening in blueberry.

Ectopic expression of CSD1 and CSD2 targeting genes of miR398 in grapevine is associated with oxidative stress tolerance

MicroRNAs (miRNAs) are endogenous small RNAs of -21 nucleotides that play an important role in diverse plant physiological processes at the post-transcriptional level by directing mRNA cleavage or translational inhibition. Previous studies have indicated that down-regulation of miR398 in response to oxidative stress allows up-regulation of the two target genes, cytosolic CSD1 and chloroplastic CSD2 (copper/zinc superoxide dismutase), resulting in protecting the plants to tolerate oxidative stress. In this study, we provide evidence that grapevine miR398 (Vv-miR398), by regulating the expression of its target genes, VvCSD1 and VvCSD2, mediates responses of grapevine to copper (Cu) stress which have been magnified due to increase in Cu-containing pesticide application. The expression of Vv-miR398 was inhibited by different concentrations of Cu stress; on the other hand, there was a steady increase in the activity of VvCSD1 and VvCSD2 genes. The function of VvCSD1 and VvCSD2 under Cu stress was thoroughly examined by overexpressing the use of the VvCSD1 and VvCSD2 in transgenic tobacco (Nicotiana tabacum). We found that both the overexpressed transgenic lines had lower Cu sensitivity and higher Cu tolerance compared with the wild type. In addition, lower levels of ROS and higher levels of SOD activities were accumulated in the transgenic lines in comparison with the wild type under the higher Cu conditions. Furthermore, these transgenic tobacco lines also recorded a higher UV and salt tolerance than the WT plants. These results suggested that overexpressing the VvCSDs will enhance the ROS-scavenging systems and protect the plant against more oxidative damage. Also, more investigations in this line are needed that would provide significant improvements in our understanding the resistance of fruit crops to environmental stress.

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.

Functional Roles of microRNAs in Agronomically Important Plants-Potential as Targets for Crop Improvement and Protection

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have recently emerged as important regulators of gene expression, mainly through cleavage and/or translation inhibition of the target mRNAs during or after transcription. miRNAs play important roles by regulating a multitude of biological processes in plants which include maintenance of genome integrity, development, metabolism, and adaptive responses toward environmental stresses. The increasing population of the world and their food demands requires focused efforts for the improvement of crop plants to ensure sustainable food production. Manipulation of mRNA transcript abundance via miRNA control provides a unique strategy for modulating differential plant gene expression and miRNAs are thus emerging as the next generation targets for genetic engineering for improvement of the agronomic properties of crops. However, a deeper understanding of its potential and the mechanisms involved will facilitate the design of suitable strategies to obtain the desirable traits with minimum trade-offs in the modified crops. In this regard, this review highlights the diverse roles of conserved and newly identified miRNAs in various food and industrial crops and recent advances made in the uses of miRNAs to improve plants of agronomically importance so as to significantly enhance crop yields and increase tolerance to various environmental stress agents of biotic-or abiotic origin.

Genome-wide characterisation and expression profile of the grapevine ATL ubiquitin ligase family reveal biotic and abiotic stress-responsive and development-related members

The Arabidopsis Tóxicos en Levadura (ATL) protein family is a class of E3 ubiquitin ligases with a characteristic RING-H2 Zn-finger structure that mediates diverse physiological processes and stress responses in plants. We carried out a genome-wide survey of grapevine (Vitis vinifera L.) ATL genes and retrieved 96 sequences containing the canonical ATL RING-H2 domain. We analysed their genomic organisation, gene structure and evolution, protein domains and phylogenetic relationships. Clustering revealed several clades, as already reported in Arabidopsis thaliana and rice (Oryza sativa), with an expanded subgroup of grapevine-specific genes. Most of the grapevine ATL genes lacked introns and were scattered among the 19 chromosomes, with a high level of duplication retention. Expression profiling revealed that some ATL genes are expressed specifically during early or late development and may participate in the juvenile to mature plant transition, whereas others may play a role in pathogen and/or abiotic stress responses, making them key candidates for further functional analysis. Our data offer the first genome-wide overview and annotation of the grapevine ATL family, and provide a basis for investigating the roles of specific family members in grapevine physiology and stress responses, as well as potential biotechnological applications.

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.

Molecular phenology in plants: in natura systems biology for the comprehensive understanding of seasonal responses under natural environments

Phenology refers to the study of seasonal schedules of organisms. Molecular phenology is defined here as the study of the seasonal patterns of organisms captured by molecular biology techniques. The history of molecular phenology is reviewed briefly in relation to advances in the quantification technology of gene expression. High-resolution molecular phenology (HMP) data have enabled us to study phenology with an approach of in natura systems biology. I review recent analyses of FLOWERING LOCUS C (FLC), a temperature-responsive repressor of flowering, along the six steps in the typical flow of in natura systems biology. The extensive studies of the regulation of FLC have made this example a successful case in which a comprehensive understanding of gene functions has been progressing. The FLC-mediated long-term memory of past temperatures creates time lags with other seasonal signals, such as photoperiod and short-term temperature. Major signals that control flowering time have a phase lag between them under natural conditions, and hypothetical phase lag calendars are proposed as mechanisms of season detection in plants. Transcriptomic HMP brings a novel strategy to the study of molecular phenology, because it provides a comprehensive representation of plant functions. I discuss future perspectives of molecular phenology from the standpoints of molecular biology, evolutionary biology and ecology.

Novel functional microRNAs from virus-free and infected Vitis vinifera plants under water stress

MicroRNAs (miRNAs) are small non-coding RNAs that regulate the post-transcriptional control of several pathway intermediates, thus playing pivotal roles in plant growth, development and response to biotic and abiotic stresses. In recent years, the grapevine genome release, small(s)-RNAseq and degradome-RNAseq together has allowed the discovery and characterisation of many miRNA species, thus rendering the discovery of additional miRNAs difficult and uncertain. Taking advantage of the miRNA responsiveness to stresses and the availability of virus-free Vitis vinifera plants and those infected only by a latent virus, we have analysed grapevines subjected to drought in greenhouse conditions. The sRNA-seq and other sequence-specific molecular analyses have allowed us to characterise conserved miRNA expression profiles in association with specific eco-physiological parameters. In addition, we here report 12 novel grapevine-specific miRNA candidates and describe their expression profile. We show that latent viral infection can influence the miRNA profiles of V. vinifera in response to drought. Moreover, study of eco-physiological parameters showed that photosynthetic rate, stomatal conductance and hydraulic resistance to water transport were significantly influenced by drought and viral infection. Although no unequivocal cause–effect explanation could be attributed to each miRNA target, their contribution to the drought response is discussed.

The Influence of Genotype and Environment on Small RNA Profiles in Grapevine Berry

Understanding the molecular mechanisms involved in the interaction between the genetic composition and the environment is crucial for modern viticulture. We approached this issue by focusing on the small RNA transcriptome in grapevine berries of the two varieties Cabernet Sauvignon and Sangiovese, growing in adjacent vineyards in three different environments. Four different developmental stages were studied and a total of 48 libraries of small RNAs were produced and sequenced. Using a proximity-based pipeline, we determined the general landscape of small RNAs accumulation in grapevine berries. We also investigated the presence of known and novel miRNAs and analyzed their accumulation profile. The results showed that the distribution of small RNA-producing loci is variable between the two cultivars, and that the level of variation depends on the vineyard. Differently, the profile of miRNA accumulation mainly depends on the developmental stage. The vineyard in Riccione maximizes the differences between the varieties, promoting the production of more than 1000 specific small RNA loci and modulating their expression depending on the cultivar and the maturation stage. In total, 89 known vvi-miRNAs and 33 novel vvi-miRNA candidates were identified in our samples, many of them showing the accumulation profile modulated by at least one of the factors studied. The in silico prediction of miRNA targets suggests their involvement in berry development and in secondary metabolites accumulation such as anthocyanins and polyphenols.

Cs-miR156 is involved in the nitrogen form regulation of catechins accumulation in tea plant (Camellia sinensis L.)

The nitrogen source affects the growth of tea plants and regulates the accumulation of catechins in the leaves. In this report, we assessed the influences of NH4(+) and NO3(-) on plant growth, catechins accumulation and associated gene expression. Compared with the preferential nitrogen source NH4(+), when NO3(-) was supplied as the sole nitrogen source, tea plants showed similar symptoms with the nitrogen-free treatments and showed lower nitrogen, free amino acid accumulation, chlorophyll content and biomass gain, indicating NO3(-) was not efficiently used by these plants. However, the total shoot catechins content was significantly higher for NO3(-) treatments than that for NH4(+) treatment or combined NH4(+)+NO3(-) treatment, suggesting that, in addition to its influence on plant growth, the nitrogen form regulated the accumulation of catechins in tea. The expression of catechins biosynthesis-related genes was associated with the regulation of catechins accumulation and composition changes mediated by nitrogen form. PAL, CHS, CHI, and DFR genes exhibited higher expression levels in plants supplied with NO3(-), in which the transcript level of DFR in the shoots was significantly correlated with the catechins content. In the end, we identified a new function for the Cs-miR156, which was drastically induced through NH4(+). Moreover, a potential mechanism of the Cs-miR156 pathway in regulating catechins biosynthesis in tea plants has been suggested, with particular respect to nitrogen forms. Cs-miR156 might repress the expression of the target gene SPL to regulate the DFR gene, which plays a vital role in catechins biosynthesis.

Determination of the precise sequences of computationally predicted miRNAs in Citrus reticulata by miR-RACE and characterization of the related target genes using RLM-RACE

MicroRNAs play vital roles in various biological and metabolic processes by regulating the expression of their target genes in model plants. Since there are limited reports on miRNAs in Citrus reticulata (Crt-miRNAs), the determination of precise sequences of miRNAs is essential to further analyze the functions of miRNAs in Citrus reticulata. Here, miR-RACE, a recently developed technique for determination of the potential miRNAs computationally, was employed to identify the precise sequences of Crt-miRNAs. Tissue- and development-specific expression of nine miRNAs were identified by quantitative RT-PCR in the leaves, stems, flowers and fruits Subsequently, 10 potential target genes were predicated for the eight Crt-miRNAs, most of which were transcription factors and disease resistance proteins. Four target genes were experimentally validated by Poly (A) polymerase-mediated 3′ rapid amplification of cDNA ends and RNA ligase-mediated 5′ rapid amplification of cDNA ends (PPM-RACE and RLM-RACE). Our findings showed that regulatory miRNAs in C. reticulata may play a key role in regulating growth, development, and response to disease. Future work is required to study the functions of miRNAs and their targets of C. reticulata.

miRVine: a microRNA expression atlas of grapevine based on small RNA sequencing

Background

miRNAs are the most abundant class of small non-coding RNAs, and they are involved in post-transcriptional regulations, playing a crucial role in the refinement of genetic programming during plant development. Here we present a comprehensive picture of miRNA regulation in Vitis vinifera L. plant during its complete life cycle. Furthering our knowledge about the post-transcriptional regulation of plant development is fundamental to understand the biology of such an important crop.

Results

We analyzed 70 small RNA libraries, prepared from berries, inflorescences, tendrils, buds, carpels, stamens and other samples at different developmental stages. One-hundred and ten known and 175 novel miRNAs have been identified and a wide grapevine expression atlas has been described. The distribution of miRNA abundance reveals that 22 novel miRNAs are specific to stamen, and two of them are, interestingly, involved in ethylene biosynthesis, while only few miRNAs are highly specific to other organs. Thirty-eight miRNAs are present in all our samples, suggesting a role in key regulatory circuit. On the basis of miRNAs abundance and distribution across samples and on the estimated correlation, we suggest that miRNA expression define organ identity. We performed target prediction analysis and focused on miRNA expression analysis in berries and inflorescence during their development, providing an initial functional description of the identified miRNAs.

Conclusions

Our findings represent a very extensive miRNA expression atlas in grapevine, allowing the definition of how the spatio-temporal distribution of miRNAs defines organ identity. We describe miRNAs abundance in specific tissues not previously described in grapevine and contribute to future targeted functional analyses. Finally, we present a deep characterization of miRNA involvement in berry and inflorescence development, suggesting a role for miRNA-driven hormonal regulation.

Electronic supplementary material

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

miRNA expression during prickly pear cactus fruit development

miRNAs are a class of small non-coding RNAs that regulate gene expression. They are involved in the control of many developmental processes, including fruit development. The increasing amount of information on miRNAs, on their expression, abundance, and conservation between various species, provides a new opportunity to study the role of miRNAs in non-model plant species. In this work, we used a combination of Northern blot and tissue print hybridization analysis to identify conserved miRNAs expressed during prickly pear cactus (Opuntia ficus indica) fruit development. Comparative profiling detected the expression of 34 miRNAs, which were clustered in three different groups that were associated with the different phases of fruit development. Variation in the level of miRNA expression was observed. Gradual expression increase of several miRNAs was observed during fruit development, including miR164. miR164 was selected for stem-loop RT-PCR and for a detailed spatial-temporal expression analysis. At early floral stages, miR164 was mainly localized in meristematic tissues, boundaries and fusion zones, while it was more homogenously expressed in fruit tissues. Our results provide the first evidence of miRNA expression in the prickly pear cactus and provide the basis for future research on miRNAs in Opuntia. Moreover, our analyses suggest that miR164 plays different roles during prickly pear cactus fruit development.

MicroRNA-based biotechnology for plant improvement

MicroRNAs (miRNAs) are an extensive class of newly discovered endogenous small RNAs, which negatively regulate gene expression at the post-transcription levels. As the application of next-generation deep sequencing and advanced bioinformatics, the miRNA-related study has been expended to non-model plant species and the number of identified miRNAs has dramatically increased in the past years. miRNAs play a critical role in almost all biological and metabolic processes, and provide a unique strategy for plant improvement. Here, we first briefly review the discovery, history, and biogenesis of miRNAs, then focus more on the application of miRNAs on plant breeding and the future directions. Increased plant biomass through controlling plant development and phase change has been one achievement for miRNA-based biotechnology; plant tolerance to abiotic and biotic stress was also significantly enhanced by regulating the expression of an individual miRNA. Both endogenous and artificial miRNAs may serve as important tools for plant improvement.

Integrated network analysis identifies fight-club nodes as a class of hubs encompassing key putative switch genes that induce major transcriptome reprogramming during grapevine development

We developed an approach that integrates different network-based methods to analyze the correlation network arising from large-scale gene expression data. By studying grapevine (Vitis vinifera) and tomato (Solanum lycopersicum) gene expression atlases and a grapevine berry transcriptomic data set during the transition from immature to mature growth, we identified a category named "fight-club hubs" characterized by a marked negative correlation with the expression profiles of neighboring genes in the network. A special subset named "switch genes" was identified, with the additional property of many significant negative correlations outside their own group in the network. Switch genes are involved in multiple processes and include transcription factors that may be considered master regulators of the previously reported transcriptome remodeling that marks the developmental shift from immature to mature growth. All switch genes, expressed at low levels in vegetative/green tissues, showed a significant increase in mature/woody organs, suggesting a potential regulatory role during the developmental transition. Finally, our analysis of tomato gene expression data sets showed that wild-type switch genes are downregulated in ripening-deficient mutants. The identification of known master regulators of tomato fruit maturation suggests our method is suitable for the detection of key regulators of organ development in different fleshy fruit crops.

Non-coding RNAs and gastric cancer

Non-coding RNAs (ncRNAs) play key roles in development, proliferation, differentiation and apoptosis. Altered ncRNA expression is associated with gastric cancer occurrence, invasion, and metastasis. Moreover, aberrant expression of microRNAs (miRNAs) is significantly related to gastric cancer tumor stage, size, differentiation and metastasis. MiRNAs interrupt cellular signaling pathways, inhibit the activity of tumor suppressor genes, and affect the cell cycle in gastric cancer cells. Some miRNAs, including miR-21, miR-106a and miR-421, could be potential markers for the diagnosis of gastric cancer. Long non-coding RNAs (lncRNAs), a new research hotspot among cancer-associated ncRNAs, play important roles in epigenetic, transcriptional and post-transcriptional regulation. Several gastric cancer-associated lncRNAs, such as CCAT1, GACAT1, H19, and SUMO1P3, have been explored. In addition, Piwi-interacting RNAs, another type of small ncRNA that is recognized by gastroenterologists, are involved in gastric carcinogenesis, and piR-651/823 represents an efficient diagnostic biomarker of gastric cancer that can be detected in the blood and gastric juice. Small interfering RNAs also function in post-transcriptional regulation in gastric cancer and might be useful in gastric cancer treatment.

Grapevine microRNAs responsive to exogenous gibberellin

Background

MicroRNAs (miRNAs), involving in various biological and metabolic processes, have been discovered and analyzed in quite a number of plants species, such as Arabidopsis, rice and other plants. However, there have been few reports about grapevine miRNAs in response to gibberelline (GA₃).

Results

Solexa technology was used to sequence small RNA libraries constructed from grapevine berries treated with GA₃ and the control. A total of 122 known and 90 novel grapevine miRNAs (Vvi-miRNAs) were identified. Totally, 137 ones were found to be clearly responsive to GA₃, among which 58 were down-regulated, 51 were up-regulated, 21 could only be detected in the control, and seven were only detected in the treatment. Subsequently, we found that 28 of them were differentially regulated by GA₃, with 12 conserved and 16 novel Vvi-miRNAs, based on the analysis of qRT-PCR essays. There existed some consistency in expression levels of GA₃-responsive Vvi-miRNAs between high throughput sequencing and qRT-PCR essays. In addition, 117 target genes for 29 novel miRNAs were predicted.

Conclusions

Deep sequencing of short RNAs from grapevine berries treated with GA3 and the control identified 137 GA₃-responsive miRNAs, among which 28 exhibited different expression profiles of response to GA₃ in the diverse developmental stages of grapevine berries. These identified Vvi-miRNAs might be involved in the grapevine berry development and response to environmental stresses.