For RNA-seq analysis, a total of 8 RNA-seq libraries were generated, comprising four biological replicates from each of the two vineyards . The libraries were constructed as previously described using the NEBNext Ultra II RNA Library Prep Kit for Illumina . Subsequently, these libraries were pooled in equal amounts and subjected to paired-end 150-base sequencing on two lanes of the NovaSeq 6000 platform at the Novogene Co., Ltd .Illumina sequencing of the multiplexed RNA-seq libraries resulted in 8 FASTQ files containing sequences, and the dataprocessing followed the methods described in our previous work . In summary, the quality of reads was assessed using FASTQ before and after trimming with Trimmomatic v0.39 . Subsequently, the trimmed reads were quantified using Salmon in non-alignment based mode to estimate transcript abundance . The transcripts were mapped to the Vitis transcriptome file “Vvinifera_457_v2.1.transcript_primaryTranscript Only.fa” extracted from Phytozome database , resulting in a mapping rate higher than 61.9% . To identify differentially expressed genes between V7 and V9 at the sampling point, we utilized the DESeq2 and EdgeR packages with default parameters . For convenience, the DEGs generated by both DESeq2 and EdgeR pipelines, large pot with drainage with a threshold of PFDR<0.05 and log2fold change > 1.5 or < –1.5, were considered as being expressed . For the analysis of Gene Ontology terms, we employed the g:Profiler website with the g:SCS multiple testing correction method, using a significance threshold of 0.05 .
Finally, to visualize the consensus result, the Web-based tool Venny was used .Co-expression network modules were constructed using the variance stabilizing transformation values and the R package WGCNA . Before analyzing the data, lowly expressed genes among all sample types were removed by DESeq2, and the remaining non-lowly expressed genes of the 8 samples were used in module construction. The co-expression modules were obtained using the default settings, except that the soft threshold power was set to 9, TOMType was set to signed, minModuleSize was set to 30, mergeCutHeight was set to 0.25, and scale-free topology fit index was set to 0.8 . A module eigengene value, which summarizes the expression profile of a given module as the first principal component, was calculated and used to evaluate the association of modules with berry biochemical characteristics of V7-berries and V9-berries at the fifth sampling time . The resultant final WGCNA matrix had 42 modules with 17,553 genes. The module membership and gene significance values were calculated, subsequently the intramodular hub genes were identified . Despite the premium fruit quality of the variety, in some cases, an undesirable taste was observed under certain unknown circumstances. To gain comprehensive insights into the development of the occasional berry astringency of Scarlet Royal and understand the underlying mechanism of this phenomenon, berries were investigated at two contrasting vineyards , both following the same commercial cultural practices. However,leaf petioles analysis of grapes from both vineyards showed considerable differences in nutrient levels, especially in the primary macronutrients .
During both seasons, the amount of nitrogen in the form of nitrate in LP-V9 was roughly 2 to 3 times higher than the normal levels, in contrast to its counterpart in LP-V7, which slightly accumulated more or less N. Similarly, LP-V9 contained higher percentages of phosphorus and potassium compared to LP-V7 . Conversely, the amounts of secondary macronutrients, calcium and magnesium , in LP-V7 were within the normal range but greater than LP-V9, which showed Mg deficiency in the first year only. Regarding the micronutrients, their levels were mainly within or around the normal range at both vineyards and during both seasons, with some differences . For example, zinc was slightly higher in LP-V9, especially in the first year. On the contrary, manganese and chlorine were roughly 2 times higher in V7 . Similarly, soil analysis shoed a higher level of nitrogen, potassium and magnesium . However, no significant difference was observed in all other soil macro and micronutrients. During the two seasons of the study, we determined the total marketable yield and the number of clusters in both vineyards. Our data revealed a higher yield in V7 compared to V9 in 2016 and 2017, respectively. The lower yield in V9 can likely be attributed to the smaller number of clusters in V9 compared to V7 during 2016 and 2017. To monitor the changes in the biochemical composition of Scarlet Royal berries, V7 and V9 berries were periodically sampled at six time points from veraison until the end of the season . The obtained data showed that berry polyphenols exhibited discernible patterns in both vineyards, most notably during the ripening stage . Of special interest were the tannin compounds, which widely affect organoleptic properties such as astringency and bitterness . Our data showed that berries from both V7 and V9 vineyards maintained lower levels of tannin from veraison up to the middle of August . Subsequently, a significant gradual increment of tannin took place. However, only V9-berries showed consistent accumulation of tannin over the two studied seasons compared to V7-berries, where the significant induction occurred only during the first season.
It is worth noting that the levels of tannin were lower in both vineyards during the second year compared to the first season. Nevertheless, they were more pronounced in V9-berries compared to V7-berries, with roughly 2- to 4.5-fold increases by the end of the harvesting time during the two seasons, respectively . The patterns of catechin and quercetin glycosides were inconsistent during both seasons, particularly within V7-berries . During the first year, for instance, the levels of catechin were similar in both vineyards, showing a dramatic increase only by the end of the season . In contrast, during the second year, such induction of catechin was exclusively restricted to V9-berries, starting from time S3 . For quercetin glycosides, V7-berries exhibited significantly higher amounts at early stages during both seasons relative to V9-berries . However, subsequent amounts were comparable in both vineyards during the first season only , but not in the second one, where V7-berries showed a significant drop at the last sample S6 . Interestingly, the levels of quercetin glycosides were roughly equal at the last V9-berries sample between both seasons despite such inconsistency. For total anthocyanins , the levels in early samples were comparable in both vineyards and seasons . Afterwards, their pattern started to vary between V7 and V9 within the same season, as well as from the first season to the second, as the nutrient amounts fluctuated as well . Nevertheless, TAC accumulation was positively correlated with the progress of ripening in V7-berries, but not V9-berries. To further confirm our data, we measured these phenolic compounds for the third time in mid-September of the next year . Overall, the results showed that the patterns of tannins and TAC were reciprocally inverted between V7-berries and V9-berries as ripening advanced. In addition, both catechin and quercetin glycosides most likely followed the pattern of tannins despite their seasonal fluctuations. To further distinguish V7-berries and V9-berries and assess their astringency development, a panel test was performed using samples at three commercial harvest times . A group of 12 nontechnical panelists scored berry astringency on a scale from 1 to 7, square pot where 1 is extremely low and 7 is extremely high. The panelists were trained using samples from contrasting standard varieties, including Flame Seedless and Crimson as non-astringent and Vintage Red known for its astringent taste . The results showed that V7-berries exhibited lower intensity of astringency compared to V9-berries . As ripening proceeded, astringency levels increased in V9-berries, but decreased in V7-berries. Moreover, we collected samples from clusters with various astringent taste and measured its tannins content. We were able to determine that the threshold level of tannins that causes the Scarlet Royal astringency taste is around 400 mg/L . Taking into account the levels of polyphenol compounds and the taste panel data together , it is evident that astringency development is positively associated with tannins’ accumulation throughout the ripening process of V9-berries. Nevertheless, organoleptic analysis revealed a significant difference in the berries of the two vineyards, particularly in terms of total soluble solids and titratable acidity . Notably, V9 berries exhibited higher titratable acidity and lower total soluble solids, especially in the later stages . It’s worth noting that the weight of V9 berries is also higher than that of V7 .To better understand the molecular events associated with the induction of tannins and astringency upon ripening, the berry transcriptome profile was analyzed in both V7-berries and V9- berries at the late commercial harvest date . Following the quality and quantity check, extracted RNA from quadruplicate samples was deeply sequenced . Of the 19.7 to 24.4 million high-quality clean reads per replicate, 61.9% to 66.1% were mapped against the V. vinifera transcriptome . Hierarchical clustering of the RNAseq data showed explicit changes in the berry transcriptome profile between V7- berries and V9-berries . The Principal Component Analysis showed high consistency among biological replicates . Samples were mainly separated along the first component , which was responsible for 97% of the variance, and was definitely associated with the site ofcultivation; V7 and V9. In contrast, the second component was trivial, accounting for only 1% of the variance and was probably attributed to experimental error. Such results were expected, as berry samples came from the same cultivar, Scarlet Royal , and the only difference between them was the vineyard locations.
To identify the differentially expressed genes in V7- berries and V9-berries at this specific time within the ripening window, the RNAseq data were analyzed using two different Bioconductor packages, DESeq2, and EdgeR . Subsequently, the DEGs with FDR < 0.05 and log2fold change > 1.5 or < –1.5 generated by both pipelines were considered . The pairwise comparison between berry transcriptomes resulted in 2134 DEGs, with 1514 up-regulated and 620 down-regulated . The data manifested the impact of the cultivation site on the transcriptional reprogramming of a large number of genes that ultimately affect berry quality. Most apparently, at the V9 vineyard, where roughly 2.5-fold higher number of berry transcripts were upregulated compared to V7 . Subsequently, the enrichment of Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways were analyzed among the up- and down-regulated DEGs using the Vitis vinifera Ensembl GeneID . Among the significantly enriched GO terms, the up-regulated transcripts in V9-berries exhibited high enrichment in the molecular function GO terms for quercetin 3-O-glucosyltransferase activity and quercetin 7-Oglucosyltransferase activity . Additionally, the V9-berries induced DEGs were highly enriched in the biological process GO terms for the jasmonic acid signaling pathway and cellular response , Lphenylalanine metabolic process , L-phenylalanine biosynthetic process , and nitrogen compound metabolic process . Similarly, these DEGs were highly enriched in the KEGG pathways for the biosynthesis of secondary metabolites and phenylpropanoid biosynthesis . On the other hand, the down-regulated transcripts in V9-berries showed substantial augmentation in the MF GO terms for hormone binding , abscisic acid binding , and potassium ion transmembrane transporter activity . Correspondingly, the BP GO terms for hormone-mediated signaling pathway and response , auxin-activated signaling, cellular response, and homeostasis , abscisic acid-activated signaling, response, and cellular response , response to strigolactone , potassium ion transmembrane transport , and potassium ion transport , as well as the KEGG pathways for plant hormone signal transduction , brassinosteroid biosynthesis , and carotenoid biosynthesis were highly enriched in the down-regulated genes of V9-berries . Overall, the transcriptome analysis pointed out the substantial changes in transcript abundance that coordinate and reflect the observed induction of tannins/astringency during the maturation and ripening of V9-berries compared to the V7-berries .To elucidate which fundamental processes were altered during tannins/astringency induction within berries, the Weighted Gene CoExpression Network Analysis was applied to construct coexpression networks. Forty-two modules were identified based on pairwise correlations among the 17553 non-lowly expressed genes . Subsequently, the biochemical data from both V7-berries and V9-berries were correlated to the WGCNA modules, and only 2 modules, M21 and M30, displayed substantial correlations with berry polyphenols, containing 5349 and 4559 genes, respectively . The M21 module was positively linked with TAC , but negatively associated with tannins, catechin, and quercetin glycosides . On the contrary, the M30 module exhibited a positive correlation with tannins, catechin, and quercetin glycosides , but was negatively linked with TAC .