Characterization of grapevine microR164 and its target genes

NAC61 regulates late- and post-ripening osmotic, oxidative, and biotic stress responses in grapevine

During late- and post-ripening stages, grape berry undergoes profound biochemical and physiological changes whose molecular control is poorly understood. Here, we report the role of NAC61, a grapevine NAC transcription factor, in regulating different processes involved in berry ripening progression. NAC61 is highly expressed during post-harvest berry dehydration and its expression pattern is closely related to sugar concentration. The ectopic expression of NAC61 in Nicotiana benthamiana leaves resulted in low stomatal conductance, high leaf temperature, tissue collapse and a higher relative water content. Transcriptome analysis of grapevine leaves transiently overexpressing NAC61 and DNA affinity purification and sequencing analyses allowed us to narrow down a list of NAC61-regulated genes. Direct regulation of the stilbene synthase regulator MYB14, the osmotic stress-related gene DHN1b, the Botrytis cinerea susceptibility gene WRKY52, and NAC61 itself was validated. We also demonstrate that NAC61 interacts with NAC60, a proposed master regulator of grapevine organ maturation, in the activation of MYB14 and NAC61 expression. Overall, our findings establish NAC61 as a key player in a regulatory network that governs stilbenoid metabolism and osmotic, oxidative, and biotic stress responses that are the hallmark of late- and post-ripening grape stages.

Dissecting the effect of soil on plant phenology and berry transcriptional plasticity in two Italian grapevine varieties (Vitis vinifera L.)

Grapevine embodies a fascinating species as regards phenotypic plasticity and genotype-per-environment interactions. The terroir, namely the set of agri-environmental factors to which a variety is subjected, can influence the phenotype at the physiological, molecular, and biochemical level, representing an important phenomenon connected to the typicality of productions. We investigated the determinants of plasticity by conducting a field-experiment where all terroir variables, except soil, were kept as constant as possible. We isolated the effect of soils collected from different areas, on phenology, physiology, and transcriptional responses of skin and flesh of a red and a white variety of great economic value: Corvina and Glera. Molecular results, together with physio-phenological parameters, suggest a specific effect of soil on grapevine plastic response, highlighting a higher transcriptional plasticity of Glera in respect to Corvina and a marked response of skin compared to flesh. Using a novel statistical approach, we identified clusters of plastic genes subjected to the specific influence of soil. These findings could represent an issue of applicative value, posing the basis for targeted agricultural practices to enhance the desired characteristics for any soil/cultivar combination, to improve vineyards management for a better resource usage and to valorize vineyards uniqueness maximizing the terroir-effect.

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.

VviNAC33 promotes organ de-greening and represses vegetative growth during the vegetative-to-mature phase transition in grapevine

Plants undergo several developmental transitions during their life cycle. In grapevine, a perennial woody fruit crop, the transition from vegetative/green-to-mature/woody growth involves transcriptomic reprogramming orchestrated by a small group of genes encoding regulators, but the underlying molecular mechanisms are not fully understood. We investigated the function of the transcriptional regulator VviNAC33 by generating and characterizing transgenic overexpressing grapevine lines and a chimeric repressor, and by exploring its putative targets through a DNA affinity purification sequencing (DAP-seq) approach combined with transcriptomic data. We demonstrated that VviNAC33 induces leaf de-greening, inhibits organ growth and directly activates the expression of STAY-GREEN PROTEIN 1 (SGR1), which is involved in Chl and photosystem degradation, and AUTOPHAGY 8f (ATG8f), which is involved in the maturation of autophagosomes. Furthermore, we show that VviNAC33 directly inhibits AUXIN EFFLUX FACILITATOR PIN1, RopGEF1 and ATP SYNTHASE GAMMA CHAIN 1T (ATPC1), which are involved in photosystem II integrity and activity. Our results show that VviNAC33 plays a major role in terminating photosynthetic activity and organ growth as part of a regulatory network governing the vegetative-to-mature phase transition.

Transcriptomic and biochemical investigations support the role of rootstock-scion interaction in grapevine berry quality

Background

In viticulture, rootstock genotype plays a critical role to improve scion physiology, berry quality and to adapt grapevine (Vitis vinifera L.) to different environmental conditions. This study aimed at investigating the effect of two different rootstocks (1103 Paulsen - P - and Mgt 101–14 - M) in comparison with not grafted plants - NGC - on transcriptome (RNA-seq and small RNA-seq) and chemical composition of berry skin in Pinot noir, and exploring the influence of rootstock-scion interaction on grape quality. Berry samples, collected at veraison and maturity, were investigated at transcriptional and biochemical levels to depict the impact of rootstock on berry maturation.

Results

RNA- and miRNA-seq analyses highlighted that, at veraison, the transcriptomes of the berry skin are extremely similar, while variations associated with the different rootstocks become evident at maturity, suggesting a greater diversification at transcriptional level towards the end of the ripening process. In the experimental design, resembling standard agronomic growth conditions, the vines grafted on the two different rootstocks do not show a high degree of diversity. In general, the few genes differentially expressed at veraison were linked to photosynthesis, putatively because of a ripening delay in not grafted vines, while at maturity the differentially expressed genes were mainly involved in the synthesis and transport of phenylpropanoids (e.g. flavonoids), cell wall loosening, and stress response. These results were supported by some differences in berry phenolic composition detected between grafted and not grafted plants, in particular in resveratrol derivatives accumulation.

Conclusions

Transcriptomic and biochemical data demonstrate a stronger impact of 1103 Paulsen rootstock than Mgt 101–14 or not grafted plants on ripening processes related to the secondary metabolite accumulations in berry skin tissue. Interestingly, the MYB14 gene, involved in the feedback regulation of resveratrol biosynthesis was up-regulated in 1103 Paulsen thus supporting a putative greater accumulation of stilbenes in mature berries.

The Critical Role of miRNAs in Regulation of Flowering Time and Flower Development

Flowering is an important biological process for plants that ensures reproductive success. The onset of flowering needs to be coordinated with an appropriate time of year, which requires tight control of gene expression acting in concert to form a regulatory network. MicroRNAs (miRNAs) are non-coding RNAs known as master modulators of gene expression at the post-transcriptional level. Many different miRNA families are involved in flowering-related processes such as the induction of floral competence, floral patterning, and the development of floral organs. This review highlights the diverse roles of miRNAs in controlling the flowering process and flower development, in combination with potential biotechnological applications for miRNAs implicated in flower regulation.

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.

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.

MicroRNA profiling analysis of developing berries for 'Kyoho' and its early-ripening mutant during berry ripening

BACKGROUND

'Fengzao' is an early-ripening bud mutant of 'Kyoho', which matures nearly 30 days earlier than 'Kyoho'. To gain a better understanding of the regulatory role of miRNAs in early-ripening of grape berry, high-throughput sequencing approach and quantitative RT-PCR validation were employed to identify miRNAs at the genome-wide level and profile the expression patterns of the miRNAs during berry development in 'Kyho' and 'Fengzao', respectively.

RESULTS

Nine independent small RNA libraries were constructed and sequenced in two varieties from key berry development stages. A total of 108 known miRNAs and 61 novel miRNAs were identified. Among that, 159 miRNAs identified in 'Fengzao' all completely expressed in 'Kyoho' and there were 10 miRNAs specifically expressed in 'Kyoho'. The expression profiles of known and novel miRNAs were quite similar between two varieties. As the major differentially expressed miRNAs, novel_144, vvi-miR3626-3p and vvi-miR3626-5p only expressed in 'Kyoho', vvi-miR399b and vvi-miR399e were down-regulated in 'Fengzao', while vvi-miR477b-3p up-regulated in 'Fengzao'. According to the expression analysis and previous reports, miR169-NF-Y subunit, miR398-CSD, miR3626-RNA helicase, miR399- phosphate transporter and miR477-GRAS transcription factor were selected as the candidates for further investigations of miRNA regulation role in the early-ripening of grape. The qRT-PCR analyses validated the contrasting expression patterns for these miRNAs and their target genes.

CONCLUSIONS

The miRNAome of the grape berry development of 'Kyoho', and its early-ripening bud mutant, 'Fengzao' were compared by high-throughput sequencing. The expression pattern of several key miRNAs and their target genes during grape berry development and ripening stages was examined. Our results provide valuable basis towards understanding the regulatory mechanisms of early-ripening of grape berry.

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.

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.

Characterization of miR061 and its target genes in grapevine responding to exogenous gibberellic acid

MicroRNAs (miRNAs), as an important growth regulator, are also involved in gibberellic acid (GA) signaling, revealing much relationship between miRNAs and GA in various plant responses. Grape is highly sensitive to GA₃, which plays a significant regulatory role in regulation of flower development, berry expansion, berry set, berry ripening, and seedlessness induction; further, it was found that grapevine miR061 (VvmiR061) is a GA₃ responsive miRNA. In this study, grapevine REV (VvREV) and HOX32 (VvHOX32), two target genes of VvmiR061, were predicted, verified, and cloned; homologous conservation was analyzed in various plants. The expression profiles of both VvmiR061 and its target genes (VvREV and VvHOX32) under GA₃ treatment were detected by qRT-PCR during grapevine flower and berry development. Results revealed that GA₃ treatment has upregulated the transcription of VvREV and VvHOX32, while it downregulated the expression of VvmiR061. The function of VvmiR061 in cleaving target genes VvREV and VvHOX32 was diminished by GA₃ treatment during flower developmental process. The results of this study exhibited the importance of VvmiR061 in regulating flower development and GA₃ signaling pathway and also contributed some to the knowledge of small RNA-mediated regulation in grape.

A characterization of grapevine of GRAS domain transcription factor gene family

GRAS domain genes are a group of important plant-specific transcription factors that have been reported to be involved in plant development. In order to know the roles of GRAS genes in grapevine, a widely cultivated fruit crop, the study on grapevine GRAS (VvGRAS) was carried out, and from which, 43 were identified from 12× assemble grapevine genomic sequences. Further, the genomic structures, synteny, phylogeny, expression profiles in different tissues of these genes, and their roles in response to stress were investigated. Among the genes, two potential target genes (VvSCL15 and VvSCL22) for VvmiR171 were experimentally verified by PPM-RACE and RLM-RACE, in that not only the cleavage sites of miR171 on the target mRNA were mapped but also the cleaved fragments and their expressing patterns were detected. Transgenic Arabidopsis plants over expression VvSCL15 showed lower tolerance to drought and salt treatments.

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.

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.

Transporters, chaperones, and P-type ATPases controlling grapevine copper homeostasis

With more copper and copper-containing compounds used as bactericides and fungicides in viticulture, copper homeostasis in grapevine (Vitis) has become one of the serious environmental crises with great risk. To better understand the regulation of Cu homeostasis in grapevine, grapevine seedlings cultured in vitro with different levels of Cu were utilized to investigate the tolerance mechanisms of grapevine responding to copper availability at physiological and molecular levels. The results indicated that Cu contents in roots and leaves arose with increasing levels of Cu application. With copper concentration increasing, malondialdehyde (MDA) content increased in roots and leaves and the activities of superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) increased to protect the plant itself from damage. The expression patterns of 19 genes, encoding transporters, chaperones, and P-type ATPases involved in copper homeostasis in root and leaf of grapevine seedling under various levels of Cu(2+) were further analyzed. The expression patterns indicated that CTr1, CTr2, and CTr8 transporters were significantly upregulated in response both to Cu excess and deficiency. ZIP2 was downregulated in response to Cu excess and upregulated under Cu-deficient conditions, while ZIP4 had an opposite expression pattern under similar conditions. The expression of chaperones and P-type ATPases in response to Cu availability in grapevine were also briefly studied.

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.

Identification of chilling stress-responsive tomato microRNAs and their target genes by high-throughput sequencing and degradome analysis

BACKGROUND

MicroRNAs (miRNAs) are a class of noncoding small RNAs (sRNAs) that are 20-24 nucleotides (nt) in length. Extensive studies have indicated that miRNAs play versatile roles in plants, functioning in processes such as growth, development and stress responses. Chilling is a common abiotic stress that seriously affects plants growth and development. Recently, chilling-responsive miRNAs have been detected in several plant species. However, little is known about the miRNAs in the model plant tomato. 'LA1777' (Solanum habrochaites) has been shown to survive chilling stress due to its various characteristics.

RESULTS

Here, two small RNA libraries and two degradome libraries were produced from chilling-treated (CT) and non-chilling-treated (NT) leaves of S. habrochaites seedlings. Following high-throughput sequencing and filtering, 161 conserved and 236 novel miRNAs were identified in the two libraries. Of these miRNAs, 192 increased in the response to chilling stress while 205 decreased. Furthermore, the target genes of the miRNAs were predicted using a degradome sequencing approach. It was found that 62 target genes were cleaved by 42 conserved miRNAs, while nine target genes were cleaved by nine novel miRNAs. Additionally, nine miRNAs and six target genes were validated by quantitative real-time PCR (qRT-PCR). Target gene functional analysis showed that most target genes played positive roles in the chilling response, primarily by regulating the expression of anti-stress proteins, antioxidant enzyme and genes involved in cell wall formation.

CONCLUSIONS

Tomato is an important model plant for basic biological research. In this study, numerous conserved and novel miRNAs involved in the chilling response were identified using high-throughput sequencing, and the target genes were analyzed by degradome sequencing. The work helps identify chilling-responsive miRNAs in tomato and increases the number of identified miRNAs involved in chilling stress. Furthermore, the work provides a foundation for further study of the regulation of miRNAs in the plant response to chilling stress.

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.

De novo transcriptome assembly, gene annotation, marker development, and miRNA potential target genes validation under abiotic stresses in Oenanthe javanica

Oenanthe javanica is an aquatic perennial herb with known medicinal properties and an edible vegetable with high vitamin and mineral content. The understanding of the biology of O. javanica is limited by the absence of information on its genome, transcriptome, and small RNA. In this study, transcriptome sequencing and small RNA sequencing were performed to annotate function genes, develop SSR markers and analyze potential target genes of miRNAs in O. javanica. All reads with total nucleotides number of 1,440,321,408 bp were assembled into 58,072 transcripts and 40,208 unigenes. A total of 1,233 SSRs were identified from O. javanica. Generated unigenes were aligned against seven databases and annotated with functions. A total of 29 potential targets were predicted. Expression of 10 miRNAs and their corresponding target genes under abiotic stresses (heat, cold, salinity, and drought) was validated. All ten miRNAs were confirmed to response to abiotic stresses. A pair of miRNA and its target gene was found. This study can serve as a valuable resource for future studies on O. javanica, which may focus on novel gene discovery, SSR development, gene mapping, and miRNA-affected processes and pathways. This can promote the development of the useful medicinal properties of O. javanica in medical science.

Advances in identification and validation of plant microRNAs and their target genes

Developments in the field of molecular biology and genetics, such as microarray, gene transfer and discovery of small regulatory RNAs, have led to significant advances in plant biotechnology. Among the small RNAs, microRNAs (miRNAs) have elicited much interest as key post-transcriptional regulators in eukaryotic gene expression. Advances in genome and transcriptome sequencing of plants have facilitated the generation of a huge wealth of sequence information that can find much use in the discovery of novel miRNAs and their target genes. In this review, we present an overview of the developments in the strategies and methods used to identify and study miRNAs, their target genes and the mechanisms by which these miRNAs interact with their target genes since the discovery of the first miRNA. The approaches discussed include both reverse and forward genetics. We observed that despite the availability of advanced methods, certain limitations ranging from the cost of materials, equipment and personnel to the availability of genome sequences for many plant species present a number of challenges for the development and utilization of modern scientific methods for the elucidation and development of miRNAs in many important plant species.

Genomic resources in fruit plants: an assessment of current status

The availability of many genomic resources such as genome sequences, functional genomics resources including microarrays and RNA-seq, sufficient numbers of molecular markers, express sequence tags (ESTs) and high-density genetic maps is causing a rapid acceleration of genetics and genomic research of many fruit plants. This is leading to an increase in our knowledge of the genes that are linked to many horticultural and agronomically important traits. Recently, some progress has also been made on the identification and functional analysis of miRNAs in some fruit plants. This is one of the most active research fields in plant sciences. The last decade has witnessed development of genomic resources in many fruit plants such as apple, banana, citrus, grapes, papaya, pears, strawberry etc.; however, many of them are still not being exploited. Furthermore, owing to lack of resources, infrastructure and research facilities in many lesser-developed countries, development of genomic resources in many underutilized or less-studied fruit crops, which grow in these countries, is limited. Thus, research emphasis should be given to those fruit crops for which genomic resources are relatively scarce. The development of genomic databases of these less-studied fruit crops will enable biotechnologists to identify target genes that underlie key horticultural and agronomical traits. This review presents an overview of the current status of the development of genomic resources in fruit plants with the main emphasis being on genome sequencing, EST resources, functional genomics resources including microarray and RNA-seq, identification of quantitative trait loci and construction of genetic maps as well as efforts made on the identification and functional analysis of miRNAs in fruit plants.

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

MicroRNAs (miRNAs) play critical regulatory roles mainly through cleaving their target mRNAs or repressing gene translation during plant development. Grapevines are among the most economically important fruit crops with available whole genome sequences. Studies on grapevine miRNAs (Vv-miRNAs) are also widely available. However, studies on the regulation mode of Vv-miRNAs on their target mRNAs during grapevine development have not been studied well, especially at the transcriptome-wide level. Here, six small RNA and mRNA libraries from various grapevine tissues were constructed for Illumina and Degradome sequencing. Subsequently, we systematically analyzed the spatiotemporal variations in the regulation of the target genes of regulation of Vv-miRNAs. In total, 242 known and 132 novel Vv-miRNAs and 193 target mRNAs were identified, including 103 target mRNAs for known and 90 target mRNAs for novel miRNAs, were validated in one or more of the tissues examined. More than 50 % of novel miRNAs were expressed exclusively in the flowers and berries, where they cleaved their target genes in a tissue-specific manner, especially, the breadth of their cleavage sites in flower tissues. Moreover, six novel miRNAs in berries responded to exogenous gibberellin and/or ethylene under a quantitative real time RT-PCR analysis, which confirmed their regulatory functions during berry development. Up to 93.6 % of the known miRNAs were highly conserved in various tissues, where their expression levels exhibited dynamic variations during grapevine development. Significantly, some Vv-miRNA families had one key member that acted as the main regulator of their target genes during grapevine development.

Characterization of regulatory mechanism of Poncirus trifoliata microRNAs on their target genes with an integrated strategy of newly developed PPM-RACE and RLM-RACE

MicroRNAs (miRNAs) play an important role in post-transcriptional gene regulation that involved various biological and metabolic processes. Many extensive studies have been done in model plant species, to discover miRNAs' regulating expression of their target genes and analyze their functions. But, the function of Poncirus trifoliata miRNAs has not been properly investigated. In this study, we employed the RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-RACE) and the newly developed method called poly (A) polymerase-mediated 3' rapid amplification of cDNA ends (PPM-RACE), which mapped the cleavage site of target mRNAs and detected expression patterns of cleaved fragments that could in turn indicate the regulatory functions of the miRNAs on their target genes. Furthermore, the spatiotemporal expression levels of target genes were analyzed by qRT-PCR, with exhibiting different expression trends from their corresponding miRNAs, thus indicating the cleavage mode of miRNAs on their target genes. The expression patterns of miRNAs, their target mRNAs and cleaved target mRNAs in different organs of juvenile and adult trifoliate orange were studied. The results showed that the expression of miRNAs and their target mRNAs was in a trade-off trend. When the miRNA expression was high, its corresponding target mRNA expression was low, while the cleaved target mRNA expression was high; when the miRNA expression was low, its target mRNA expression was high, while the expression of cleaved target mRNAs follows that of the miRNA. The validation of the cleavage site of target mRNAs and the detection of expression patterns of cleaved fragments can further broaden the knowledge of small RNA-mediated regulation in P. trifoliate.

Characterization of target mRNAs for grapevine microRNAs with an integrated strategy of modified RLM-RACE, newly developed PPM-RACE and qPCRs

MicroRNAs (miRNAs) regulate target gene expression by mediating target gene cleavage or inhibition of translation at transcriptional and post-transcriptional levels in higher plants. Until now, many grapevine microRNAs (Vv-miRNAs) have been identified and quite a number of miRNA target genes were also verified by various analysis. However, global interaction of miRNAs with their target genes still remained to perform more research. We reported experimental validation of a number of miRNA target genes in table grapevine that had been previously identified by bioinformatics in our earlier studies. To verify more predicted target genes of Vv-miRNAs and elucidate the modes by which these Vv-miRNAs work on their target genes, 31 unverified potential target genes for 18 Vv-miRNAs were experimentally verified by a new integrated strategy employing a modified 5'-RLM-RACE (RNA ligase-mediated 5' rapid amplification of cDNA ends), 3'-PPM-RACE (poly(A) polymerase-mediated 3' rapid amplification of cDNA ends) and qRT-PCRs of cleavage products. The results showed that these Vv-miRNAs negatively regulated expression of their target messenger RNAs (mRNAs) through guiding corresponding target mRNA cleavage, of which about 94.4% Vv-miRNAs cleaved their target mRNAs mainly at the tenth nucleotide of 5'-end of miRNAs. Expression levels of both miRNAs and their target mRNAs in eight tissues exhibited inverse relationships, and expressions both of cleaved targets and miRNAs indicated a cleavage mode of Vv-miRNAs on their target genes. Our results confirm the importance of Vv-miRNAs in grapevine growth and development, and suggest more study on Vv-miRNAs and targets can enrich the knowledge of miRNA mediated-regulation in grapevine.

Evolutionary conservation of microRNA regulatory programs in plant flower development

MicroRNAs (miRNAs) are post-transcriptional regulators of growth and development in both plants and animals. Flowering is critical for the reproduction of angiosperms. Flower development entails the transition from vegetative growth to reproductive growth, floral organ initiation, and the development of floral organs. These developmental processes are genetically regulated by miRNAs, which participate in complex genetic networks of flower development. A survey of the literature shows that miRNAs, their specific targets, and the regulatory programs in which they participate are conserved throughout the plant kingdom. This review summarizes the role of miRNAs and their targets in the regulation of gene expression during the floral developmental phase, which includes the floral transition stage, followed by floral patterning, and then the development of floral organs. The conservation patterns observed in each component of the miRNA regulatory system suggest that these miRNAs play important roles in the evolution of flower development.

Genome-wide analysis of 24-nt siRNAs dynamic variations during rice superior and inferior grain filling

24 nt-siRNAs are the most abundant small interfering RNAs in rice grains aside from microRNAs. To investigate the roles that 24 nt-siRNAs played in the poor grain filling of rice inferior grains, dynamic variations of 24 nt-siRNAs in inferior grains were compared with those of superior grains by using small RNA deep sequencing technology. The results showed that 24 nt-siRNAs derived from multiple regions of rice genome, and the maintenance of the two strands of 24 nt-siRNA duplex was a non-random process. The amounts of 24 nt-siRNAs declined with the process of grain filling in both superior and inferior grains, but 24 nt-siRNAs in inferior grains was much higher than that of superior grains in each period we sampled. Bioinformatics prediction indicated that 24 nt-siRNAs targeted on more genes involved in most of the known KEGG rice pathways, such as the starch and sucrose biosynthesis pathway. Combined with digital gene expression profiling of target genes, 24 nt-siRNAs mapped on the antisense strands of exons were specifically investigated, but the abundance of 24 nt-siRNAs did not show negative correlations with their corresponding target genes. The results indicated that 24 nt-siRNAs were not involved in down-regulation of target genes. The potential biological meanings for this inconsistency were probably the results of methylation directed gene expression activation, or competition for small RNA stability methylation.