The grapem1355 candidate maps on chromosome 3, exactly on thefirst exon of its target , in a region where another two isoforms of the same gene are located . The last target of this miRNA candidate, codes a cinnamyl alcohol dehydrogenase known to be involved in the lignin biosynthesis . Other novel vvi-miRNA candidates seem to be involved in cell proliferation and in chloroplasts-related functions . Furthermore, for the new vvi-miRC482b candidate, besides the already known involvement of this miRNA family with disease resistance also predicted here, one predicted target encodes an anthocyanin 5-aromatic acyltransferase-like protein known to be involved in the biosynthesis of anthocyanin in different species . As for the conserved known vvi-miRNAs, most of the well-established examples of miR-targets, such as miR156-SPB, miR166-HD-ZIP, miR171-GRAS, miR172-AP2, confirmed in several plant species and already predicted in grapevine, blueberry containers were also predicted here. We studied miRNA profile of accumulation in the different samples.
Using their normalized abundance , i.e., their relative cloning frequency, we set an empirical cut off value equal to at least 10 TP4M in both biological replicates to consider a miRNA as expressed in a given library. Also, a miRNA was considered specific when it was expressed in one or more libraries of a unique cultivar, unique environment or unique developmental stage. According to our established cut off, 175 miRNAs were classified as expressed in at least one of our libraries . The libraries constructed from Sangiovese berries at bunch closure collected in Bolgheri showed only 24 expressed miRNAs . For all the other libraries, expressed miRNAs ranged from 76 to 148 . We found very few miRNAs specific to a given condition. The number of specific miRNAs for each cultivar, developmental stage and environment is reported in Figures 8A–C, respectively. Thirty-nine vvi-miRNAs were highly expressed in almost all libraries [21 ubiquitous plus 18 expressed in all libraries except in Bol_SG_bc ], whereas other miRNAs had different accumulation patterns. The normalized expression values of miRNAs were subjected to hierarchical clustering and represented in a heat map . To examine the relatedness among cultivars, environments and developmental stages, we generated a correlation dendrogram . The dendrogram shows, as already suggested by the heatmaps, that a fundamental dichotomy emerges between ripened and green berries. The most evident pattern of expression is observed when comparing different developmental stages, and confirm previous observation of miRNA modulation during fruit ripening .
For example, some members of the vvi-miRNA156 family were highly expressed in all ripened berries, but weakly or not expressed in green berries. Differently, vvi-miR396a-3p and vvi-miR396b-3p showed the opposite profile. Similarly, vvi-miR172d, vvi-miR166b-5p, vvi-miR166f-5p, and vvi-miR396d-5p were highly expressed in green berries but weakly expressed in ripened berries and the members of the vvi-miR319 family showed a gradient of decreasing abundance from pea size to harvest.To gain statistical evidence of miRNA differential expression driven by the environment and/or genotype, we made pairwise comparisons, keeping constant the developmental stage, and evaluating the miRNA modulation among vineyards or between cultivars . The analyses reveal that some miRNAs are differentially expressed between the two genotypes grown in the same environment, but also that a number of miRNAs are modulated by the environment. In particular the number of differentially expressed miRNAs is higher in ripened berries , while no miRNAs are differentially expressed at bunch closure stage . In details, 14 reads are differentially expressed at pea size stage, in at least one comparison, corresponding to 6 distinct miRNA families; 27 reads are modulated at 19 ◦Brix stage, corresponding to 12 miRNA families and 35 reads are differentially expressed in berries at harvest, corresponding to 12 miRNA families. It is worth noting that 4 of the 6 families modulated in the berries at pea size, are still present among the miRNAs differentially expressed in the berries sampled at 19 ◦Brix and at harvest , even though not always in the same comparisons. Some of the modulated miRNAs, both novel and known are intriguingly connected to berry development and secondary metabolism, even though most of the modulated families are still uncharacterized, or with targets not clearly involved in berry ripening and development, and deserve further studies to fully understand their biological roles. Using high throughput sequencing coupled with robust bio-informatics pipelines we analyzed small RNAs derived from the berries of Cabernet Sauvignon and Sangiovese, grown sideby-side in three vineyards, representative of different grapevine cultivation areas in Italy . We obtained nearly 750 MB reads comprising a significant proportion of small RNAs.
The size distribution profiles of our libraries were in general consistent with previous reports in berry grapevine, where the 21-nt class was more abundant than the 24-nt class . Our analysis revealed dynamic features of the regulatory network mediated by miRNAs and other small RNAs, at the basis of genotype-environment interactions. Plants evolved a series of pathways that generate small RNAs of different sizes with dedicated functions . Although the various small RNA classes have been intensively studied, we are still far from understanding how many small RNA pathways exist, and how they are connected . Additionally, new classes of small non-coding RNAs continue to be discovered and many studies demonstrate a substantial redundancy and cross-talk between known small RNA pathways . Estimating the exact percentage of the plant genome covered by small RNA-generating loci still remains a challenge. By applying static cluster analysis, we investigated small RNA abundances across the genome, identifying 4408 small RNAs producing hotspots. We analyzed their expression in different cultivars, environments and developmental stages, highlighting that the majority of the considered small RNA producing regions was modulated in different conditions. This suggests a strong influence of small RNAs in the response to environment in grapevine berries. Only 462 small RNA-generating loci, corresponding to about 10% of the total, were expressed in all the analyzed libraries, possibly involved in essential biological pathways. Comparing the two cultivars, we observed, with few exceptions, that Cabernet Sauvignon berries have a higher number of expressed sRNA-generating loci than Sangiovese berries when collected in the same conditions . Considering the fact that small RNAs are implicated in the regulation of gene expression in several processes , the higher number of small RNAs expressed in Cabernet Sauvignon compared to Sangiovese berries may reflect a buffering effect of small RNAs influencing grapevine response to diverse growing environments. We believe that these characteristics may have contributed to the wide diffusion of Cabernet Sauvignon, allowing its wide cultivation in almost all wine producing countries. This is not the case for Sangiovese whose cultivation is more restricted. It is worth noting that Sangiovese is considered a very unsettled grapevine cultivar , showing a wide range of variability in response to year, clone and bunch exposure . Differently, Cabernet Sauvignon is a cultivars showing less inter-annual differences in terms, for example, best indoor plant pots of concentration of secondary metabolites . To better evaluate varietal differences in response to the environment, we calculated the CS/SG ratio for the small RNA producing hotspots in the three vineyards. A region on chromosome 4 showed a 390-fold change in the small RNA abundance, when comparing Cabernet vs. Sangiovese . Most of the reads produced in this region are 21 nt long and are also phased in intervals of 21 nt from both strands, typical of a phased locus . The gene in this locus, also known as VvRD22g, encodes a BURP domain containing protein, involved in an ABA-mediated abiotic stress response, which persists still after long periods of stress . The small RNAs profile suggests that the locus is regulated by phased siRNAs similarly to the mechanisms already described for PPR, NB-LRR, and MYB gene families . This is a clear example of GxE interactions since the BURP domain gene modulates phased siRNAs production in the two cultivars only when grown in Riccione. When removing the threshold of minimum cluster abundance set to 5 HNA, in the CS/SG ratio, a high number of clusters with fold change greater than 50 was found, where one of the libraries has 0 HNA and the other any number greater than 30 HNA. This fact suggests a very strong modulation of the expression of small RNAs between the two cultivars, which is more or less pronounced depending on the vineyard where the berries were cultivated. A similar situation was observed comparing the expression level of small RNAs between reciprocal hybrids of Solanum lycopersicum and S. pimpinellifolium .
The ripening process of grapevine berries is highly affected by the environment and we observed the impact of the environment on the ripening process in the expression of small RNAs. The most relevant observation is that Riccione is very peculiar in relation to the activation of sRNA hotspots, as indicated by the high number of Riccione specific clusters and by the extreme modification it induces in the CS/SG ratio : in Riccione in fact this ratio decreases in green berries and increases in ripened berries, and this is not observed in any other vineyard; in addition to this the already discussed example of BURP domain gene, is observed in Riccione, as well. Riccione is the most diverse environment when compared to Montalcino and Bolgheri. Riccione is located at the Adriatic coast and has a temperate sub-littoral climate, while Montalcino and Bolgheri are both located in Tuscany with typically Mediterranean climate. Moreover, both cultivars show a peculiar profile of small RNA loci during berries ripening, in Riccione. The expression of small RNA loci in Cabernet Sauvignon berries drastically changed during development, especially when collected in Riccione , not only in the number of active loci but also in the different genic or intergenic disposition: ripened berries have a 2.6-fold increase in small RNA loci active in genic regions. Differently, when Sangiovese is grown in Riccione, there is a very high number of small RNA loci active in green berries, mainly associated to transposable elements that remains almost stable during development although the proportion of intergenic loci is reduced. Sangiovese berries collected in Montalcino show a 2.5- fold increase of small RNA producing loci during development. Differences during berry development between the cultivars may explain their different behavior in different environments, and the characteristics of each vineyard may favor one or other variety according to their demands. For example, Sangiovese needs a long growing season with sufficient warmth to fully ripen . Consequently, cooler environments will require a reprograming of Sangiovese gene expression in order to achieve ripening. Other factors such as composition of soil, level of humidity, photo period and density of cultivation may be exerting the same influence on the ripening of the berries triggering the activation of different small RNA loci.Applying a conservative pipeline to the analysis of our 48 small RNA libraries, we recognized 89 known and annotated grapevine miRNAs. In addition, when compared to previous reports in grapevine we identified 7 completely novel miRNAs plus 26 homologous to other plant species, but novel to grapevine. This is a remarkable number considering the stringency of our pipeline and that our study is based only on four developmental stages of berries. The outline of miRNA accumulation across samples is different from that of sRNA-producing loci. While the expression of sRNA-generating regions allows distinguishing very well between ripened and green berries and also between cultivars , the accumulation of miRNAs shows a clear distinction only between ripened and green berries, and when the berries were green, we observe a further dichotomy separating the two cultivars and the two green developmental stages. The same pattern of miRNA accumulation among green and ripened berries of grapevine was observed when we described the miRNA expression atlas of Vitis vinifera . Comparing the distribution of miRNAs expressed throughout our samples, we found a set of 39 miRNAs ubiquitous or nearly ubiquitous to all the libraries, and very few miRNAs specific of a cultivar, vineyard or developmental stage. All these 39 miRNAs belong to known vvi-miRNA families. With few exceptions, the same set of miRNAs was also found expressed in all the small RNA libraries constructed with different tissues of the grapevine cv. Corvina , where the population of expressed miRNAs appears highly variable apart from a well-defined group of miRNAs, probably related to the basal metabolism. These findings are also consistent with previous report in grapevine where a small number of known tissue-specific miRNAs was described .