Epigenetics has been proposed as crucial in shaping plant phenotypic plasticity, putatively explaining the rapid and reversible alterations in gene expression in response to environmental changes. This fine-tuning of gene expression can be achieved through DNA methylation, histone modifications and chromatin remodeling . Small non-coding RNAs are ubiquitous and adjustable repressors of gene expression across a broad group of eukaryotic species and are directly involved in controlling, in a sequence specific manner, multiple epigenetic phenomena such as RNA-directed DNA methylation and chromatin remodeling and might play a role in genotype by environment interactions. In plants, small ncRNAs are typically 20–24 nt long RNA molecules and participate in a wide series of biological processes controlling gene expression via transcriptional and post-transcriptional regulation . Moreover, small RNAs have been recently shown to play an important role in plants environmental plasticity . Fruit maturation, the process that starts with fruit-set and ends with fruit ripening , has been largely investigated in fleshy fruits such as tomato and grapevine. These studies highlighted, among others, plastic planters the vast transcriptomic reprogramming underlying the berry ripening process , the extensive plasticity of berry maturation in the context of a changing environment , and the epigenetic regulatory network which contributes to adjust gene expression to internal and external stimuli .
In particular, small RNAs, and especially microRNAs , are involved, among others, in those biological processes governing fruit ripening . In this work, we assessed the role of small ncRNAs in the plasticity of grapevine berries development, by employing next-generation sequencing. We focused on two cultivars of Vitis vinifera, Cabernet Sauvignon, and Sangiovese, collecting berries at four different developmental stages in three Italian vineyards, diversely located. First, we described the general landscape of small RNAs originated from hotspots present along the genome, examining their accumulation according to cultivars, environments and developmental stages. Subsequently, we analyzed miRNAs, identifying known and novel miRNA candidates and their distribution profiles in the various samples. Based on the in silico prediction of their targets, we suggest the potential involvement of this class of small RNAs in GxE interactions. The results obtained provide insights into the complex molecular machinery that connects the genotype and the environment.RNA extraction was performed as described in Kullan et al. . Briefly, total RNA was extracted from 200 mg of ground berries pericarp tissue using 1 ml of Plant RNA Isolation Reagent following manufacturer’s recommendations. The small RNA fraction was isolated from the total RNA using the mirPremier R microRNA Isolation kit and dissolved in DEPC water. All the steps suggested in the technical bulletin for small RNA isolation of plant tissues were followed except the “Filter Lysate” step, which was omitted.
The quality and quantity of small RNAs were evaluated by a NanoDrop 1000 spectrometer , and their integrity assessed by an Agilent 2100 Bioanalyzer using a small RNA chip according to the manufacturer’s instructions. Small RNA libraries were prepared using the TruSeq Small RNA Sample Preparation Kit , following all manufacturers’ instructions. Forty-eight bar-coded small RNA libraries were constructed starting from 50 ng of small RNAs. The quality of each library was assessed using an Agilent DNA 1000 chip for the Agilent 2100 Bioanalyzer. Libraries were grouped in pools with six libraries each . The pools of libraries were sequenced on an Illumina Hiseq 2000 at IGA Technology Services . The sequencing data were submitted to GEO–NCBI under the accession number GSE85611.In order to investigate whether the overall distribution and accumulation of small RNA is affected by the interaction between different V. vinifera genotypes [Cabernet Sauvignon and Sangiovese ] and environments [Bolgheri , Montalcino and Riccione ], we investigated the regions in the grapevine genome from where a high number of small RNAs were being produced , by applying a proximity-based pipeline to group and quantify clusters of small RNAs as described by Lee et al. . The nuclear grapevine genome was divided in 972,413 adjacent, non-overlapping, fixed-size windows or clusters. To determine the small RNA cluster abundance, we summed the hits-normalized-abundance values of all the small RNAs mapping to each of the 500 bp clusters, for each library .
To reduce the number of false positives, we considered a cluster as expressed when the cluster abundance was greater than the threshold for a given library, eliminating regions where few small RNAs were generated, possibly by chance. Libraries from bunch closure, representing green berries, and 19 ◦Brix representing ripened berries, where used in this analysis. From the 972,413 clusters covering the whole grapevine genome, 4408 were identified as expressed in at least one sample. As showed in Figure 1, CS-derived libraries have a higher number of expressed clusters when compared to SG-derived libraries of the same developmental stage and from the same vineyard. The exceptions were the Sangiovese green berries collected in Riccione and Sangiovese ripened berries collected in Montalcino, which have a higher number of expressed clusters than the respective CS ones. The two cultivars show a completely different small RNA profile across environments. When Cabernet berries were green, a higher number of sRNA-generating regions were found active in Bolgheri than in Montalcino and Riccione. Differently, ripened berries had the highest number of sRNA producing regions expressed in Riccione, while Bolgheri and Montalcino show a similar level of expressed clusters . Sangiovese green berries instead show the highest number of active sRNA-generating regions in Riccione, and this number is twice the number found in Bolgheri and Montalcino that is similar. Ripened berries collected in Montalcino and Riccione show almost the same high level of sRNA-generating clusters, whereas those collected in Bogheri present a lower number . We also noted that when cultivated in Bolgheri, neither Cabernet Sauvignon or Sangiovese change dramatically the number of expressed clusters during ripening, while in Riccione Cabernet Sauvignon shows a 2-fold increase of sRNAproducing clusters, which is not observed in Sangiovese. Next, the small RNA-generating clusters were characterized on the basis of the genomic regions where they map, i.e., genic, intergenic and transposable elements. In general, when the berries were green, the numbers of sRNA-generating loci located in genic and intergenic regions were roughly equal in all environments and for both cultivars, except for Sangiovese berries collected in Riccione, plastic nursery plant pot which show a slight intergenic disposition of sRNA-producing regions . Differently, in ripened berries on average 65% of the sRNA-generating loci were in genic regions, indicating a strong genic disposition of the sRNA-producing clusters . The shift of sRNA-producing clusters from intergenic to mostly genic is more pronounced in Cabernet Sauvignon berries collected in Riccione, with an increase of approximately 20% of expressed clusters in genic regions when berries pass from the green to the ripened stage. When comparing the clusters abundance among libraries, we found that 462 clusters were expressed in all libraries. The remaining 3946 expressed clusters were either shared among groups of libraries or specific to unique libraries. Interestingly, 1335 of the 4408 expressed clusters were specific to Riccione-derived libraries . The other two environments showed a much lower percentage of specific clusters, 263 and 140 in Bolgheri and Montalcino respectively . Comparing the expressed clusters between cultivars or developmental stages, we did not observe a similar discrepancy of specific clusters toward one cultivar or developmental stage; roughly the same proportion of specific clusters was found for each cultivar and for each developmental stage .
Among the 1335 specific clusters of Riccione, 605 were specific to Cabernet Sauvignon ripened berries of and 499 to Sangiovese green berries. Other smaller groups of expressed clusters were identified as specific to one cultivar, one developmental stage or also one cultivar in a specific developmental stage. When comparing the expressed clusters with the presence of transposable elements annotated in the grapevine genome , we noticed that approximately 23% of the sRNA-generating regions were TE-associated. Sangiovese green berries from Riccione have the highest proportion of TE-associated expressed clusters, while Cabernet Sauvignon ripened berries also from Riccione show the lowest proportion of TE associated expressed clusters. Sangiovese berries have the highest percentage of expressed clusters located in TE when cultivated in Riccione, compared to the other two vineyards. Interestingly, Cabernet Sauvignon berries show the lowest proportion of TE-associated clusters when growing in Riccione , independently from the berry stage. In all the libraries, Long Terminal Repeat retrotransposons were the most represented TE. More specifically, the gypsy family was the LTR class associated with the highest number of sRNA hotspots. The other classes of TE associated with the sRNA-generating regions can be visualized in Figure 3B.To determine the global relationship of small RNA-producing loci in the different environments, cultivars and developmental stages, we performed a hierarchical clustering analysis. As showed in Figure 4, the libraries clearly clustered according to the developmental stage and cultivar and not according to the environments. Ripened and green berries had their profile of sRNA-generating loci clearly distinguished from each other. Inside each branch of green and ripened samples, Cabernet Sauvignon and Sangiovese were also well separated, indicating that, the cultivar and the stage of development in which the berries were sampled modulate the outline of sRNA-producing loci more than the environment. Notwithstanding the evidence that developmental stage and variety have the strongest effect in terms of distinguishing samples clustering, we were interested to verify the environmental influence on small RNA loci expression in the two cultivars. Thus, for each sRNA-generating cluster we calculated the ratio between cluster abundance in Cabernet Sauvignon and Sangiovese in each environment and developmental stage, thereby revealing the genomic regions with regulated clusters, considering a 2-fold change threshold, a minimum abundance of 5 HNA in each library and a minimum sum of abundance of 30 HNA . Figure 5 shows how different environments affect the production of small RNAs. In Bolgheri, regardless the developmental stage there were many clusters with a very high abundance level in Cabernet Sauvignon . In Montalcino and even more in Riccione we also observed differences between the expressions of clusters in the two cultivars, with ripened and green berries showing an almost opposite profile in terms of number of clusters more expressed in Cabernet Sauvignon or Sangiovese. When the berries were green, in Montalcino Cabernet Sauvignon shows the highest number of up-regulated clusters, while in Riccione, Sangiovese has the highest number of up-regulated clusters. The opposite behavior was noticed in ripened berries, with Sangiovese having the highest number of up-regulated clusters in Montalcino and Cabernet Sauvignon in Riccione . Notably, we observed a small percentage of regulated clusters exhibiting at least a 10-fold higher abundance of small RNA in Cabernet Sauvignon or Sangiovese when compared to each other . An examination of those clusters showed that a substantial difference could exist between the cultivars, depending on the vineyard and the developmental stage. For example, in Riccione, a cluster matching a locus encoding a BURP domain-containing protein showed a fold change of 390 when comparing green berries of Sangiovese with Cabernet Sauvignon. Similarly, the majority of the highly differentially expressed clusters showed a similar profile: strong bias toward 21-nt sRNAs and a low strand bias. These findings suggest that these small RNAs might be the productof RDR polymerase activity rather than degradation products of mRNAs.We applied a pipeline adapted from Jeong et al. and Zhai et al. to identify annotated vvi-miRNAs, their variants, novel species-specific candidates and, when possible, the complementary 3p or 5p sequences. Starting from 25,437,525 distinct sequences from all the 48 libraries, the first filter of the pipeline removed sequences matching t/r/sn/snoRNAs as well as those that did not meet the threshold of 30 TP4M in at least one library or, conversely, that mapped in more than 20 loci of the grapevine genome . Only sequences 18–26-nt in length were retained. Overall, 27,332 sequences, including 56 known vvi-miRNAs, passed through this first filter and were subsequently analyzed by a modified version of miREAP as described by Jeong et al. . miREAP identified 1819 miRNA precursors producing 1108 unique miRNA candidates, including 47 known vvi-miRNA. Next, the sequences were submitted to the third filter to evaluate the single-strand and abundance bias retrieving only one or two most abundant miRNA sequence for each precursor previously identified. A total of 150 unique miRNA corresponding to 209 precursors were identified as candidate miRNAs.