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.

Exploring large-scale gene coexpression networks in peach (Prunus persica L.): a new tool for predicting gene function

Peach is a model for Prunus genetics and genomics, however, identifying and validating genes associated to peach breeding traits is a complex task. A gene coexpression network (GCN) capable of capturing stable gene-gene relationships would help researchers overcome the intrinsic limitations of peach genetics and genomics approaches and outline future research opportunities. In this study, we created four GCNs from 604 Illumina RNA-Seq libraries. We evaluated the performance of every GCN in predicting functional annotations using an algorithm based on the 'guilty-by-association' principle. The GCN with the best performance was COO300, encompassing 21 956 genes. To validate its performance predicting gene function, we performed two case studies. In case study 1, we used two genes involved in fruit flesh softening: the endopolygalacturonases PpPG21 and PpPG22. Genes coexpressing with both genes were extracted and referred to as melting flesh (MF) network. Finally, we performed an enrichment analysis of MF network and compared the results with the current knowledge regarding peach fruit softening. The MF network mostly included genes involved in cell wall expansion and remodeling, and with expressions triggered by ripening-related phytohormones, such as ethylene, auxin, and methyl jasmonate. In case study 2, we explored potential targets of the anthocyanin regulator PpMYB10.1 by comparing its gene-centered coexpression network with that of its grapevine orthologues, identifying a common regulatory network. These results validated COO300 as a powerful tool for peach and Prunus research. This network, renamed as PeachGCN v1.0, and the scripts required to perform a function prediction analysis are available at https://github.com/felipecobos/PeachGCN.

The serine-glycine-one-carbon metabolic network orchestrates changes in nitrogen and sulfur metabolism and shapes plant development

L-serine (Ser) and L-glycine (Gly) are critically important for the overall functioning of primary metabolism. We investigated the interaction of the phosphorylated pathway of Ser biosynthesis (PPSB) with the photorespiration-associated glycolate pathway of Ser biosynthesis (GPSB) using Arabidopsis thaliana PPSB-deficient lines, GPSB-deficient mutants, and crosses of PPSB with GPSB mutants. PPSB-deficient lines mainly showed retarded primary root growth. Mutation of the photorespiratory enzyme Ser-hydroxymethyltransferase 1 (SHMT1) in a PPSB-deficient background resumed primary root growth and induced a change in the plant metabolic pattern between roots and shoots. Grafting experiments demonstrated that metabolic changes in shoots were responsible for the changes in double mutant development. PPSB disruption led to a reduction in nitrogen (N) and sulfur (S) contents in shoots and a general transcriptional response to nutrient deficiency. Disruption of SHMT1 boosted the Gly flux out of the photorespiratory cycle, which increased the levels of the one-carbon (1C) metabolite 5,10-methylene-tetrahydrofolate and S-adenosylmethionine. Furthermore, disrupting SHMT1 reverted the transcriptional response to N and S deprivation and increased N and S contents in shoots of PPSB-deficient lines. Our work provides genetic evidence of the biological relevance of the Ser-Gly-1C metabolic network in N and S metabolism and in interorgan metabolic homeostasis.

MYB24 orchestrates terpene and flavonol metabolism as light responses to anthocyanin depletion in variegated grape berries

Variegation is a rare type of mosaicism not fully studied in plants, especially fruits. We examined red and white sections of grape (Vitis vinifera cv. 'Béquignol') variegated berries and found that accumulation of products from branches of the phenylpropanoid and isoprenoid pathways showed an opposite tendency. Light-responsive flavonol and monoterpene levels increased in anthocyanin-depleted areas in correlation with increasing MYB24 expression. Cistrome analysis suggested that MYB24 binds to the promoters of 22 terpene synthase (TPS) genes, as well as 32 photosynthesis/light-related genes, including carotenoid pathway members, the flavonol regulator HY5 HOMOLOGUE (HYH), and other radiation response genes. Indeed, TPS35, TPS09, the carotenoid isomerase gene CRTISO2, and HYH were activated in the presence of MYB24 and MYC2. We suggest that MYB24 modulates ultraviolet and high-intensity visible light stress responses that include terpene and flavonol synthesis and potentially affects carotenoids. The MYB24 regulatory network is developmentally triggered after the onset of berry ripening, while the absence of anthocyanin sunscreens accelerates its activation, likely in a dose-dependent manner due to increased radiation exposure. Anthocyanins and flavonols in variegated berry skins act as effective sunscreens but for different wavelength ranges. The expression patterns of stress marker genes in red and white sections of 'Béquignol' berries strongly suggest that MYB24 promotes light stress amelioration but only partly succeeds during late ripening.

Aggregated gene co-expression networks predict transcription factor regulatory landscapes in grapevine

Gene co-expression networks (GCNs) have not been extensively studied in non-model plants. However, the rapid accumulation of transcriptome datasets in certain species represents an opportunity to explore underutilized network aggregation approaches. In fact, aggregated GCNs (aggGCNs) highlight robust co-expression interactions and improve functional connectivity. We applied and evaluated two different aggregation methods on public grapevine RNA-Seq datasets from three different tissues (leaf, berry, and 'all organs'). Our results show that co-occurrence-based aggregation generally yielded the best-performing networks. We applied aggGCNs to study several transcription factor gene families, showing their capacity for detecting both already-described and novel regulatory relationships between R2R3-MYBs, bHLH/MYC, and multiple specialized metabolic pathways. Specifically, transcription factor gene- and pathway-centered network analyses successfully ascertained the previously established role of VviMYBPA1 in controlling the accumulation of proanthocyanidins while providing insights into its novel role as a regulator of p-coumaroyl-CoA biosynthesis as well as the shikimate and aromatic amino acid pathways. This network was validated using DNA affinity purification sequencing data, demonstrating that co-expression networks of transcriptional activators can serve as a proxy of gene regulatory networks. This study presents an open repository to reproduce networks in other crops and a GCN application within the Vitviz platform, a user-friendly tool for exploring co-expression relationships.

The complete reference genome for grapevine (Vitis vinifera L.) genetics and breeding

Grapevine is one of the most economically important crops worldwide. However, the previous versions of the grapevine reference genome tipically consist of thousands of fragments with missing centromeres and telomeres, limiting the accessibility of the repetitive sequences, the centromeric and telomeric regions, and the study of inheritance of important agronomic traits in these regions. Here, we assembled a telomere-to-telomere (T2T) gap-free reference genome for the cultivar PN40024 using PacBio HiFi long reads. The T2T reference genome (PN_T2T) is 69 Mb longer with 9018 more genes identified than the 12X.v0 version. We annotated 67% repetitive sequences, 19 centromeres and 36 telomeres, and incorporated gene annotations of previous versions into the PN_T2T assembly. We detected a total of 377 gene clusters, which showed associations with complex traits, such as aroma and disease resistance. Even though PN40024 derives from nine generations of selfing, we still found nine genomic hotspots of heterozygous sites associated with biological processes, such as the oxidation-reduction process and protein phosphorylation. The fully annotated complete reference genome therefore constitutes an important resource for grapevine genetic studies and breeding programs.

An improved reference of the grapevine genome reasserts the origin of the PN40024 highly-homozygous genotype

The genome sequence of the diploid and highly homozygous V. vinifera genotype PN40024 serves as the reference for many grapevine studies. Despite several improvements to the PN40024 genome assembly, its current version PN12X.v2 is quite fragmented and only represents the haploid state of the genome with mixed haplotypes. In fact, being nearly homozygous, this genome contains several heterozygous regions that are yet to be resolved. Taking the opportunity of improvements that long-read sequencing technologies offer to fully discriminate haplotype sequences, an improved version of the reference, called PN40024.v4, was generated. Through incorporating long genomic sequencing reads to the assembly, the continuity of the 12X.v2 scaffolds was highly increased with a total number decreasing from 2,059 to 640 and a reduction in N bases of 88%. Additionally, the full alternative haplotype sequence was built for the first time, the chromosome anchoring was improved and the number of unplaced scaffolds was reduced by half. To obtain a high-quality gene annotation that outperforms previous versions, a liftover approach was complemented with an optimized annotation workflow for Vitis. Integration of the gene reference catalogue and its manual curation have also assisted in improving the annotation, while defining the most reliable estimation of 35,230 genes to date. Finally, we demonstrated that PN40024 resulted from nine selfings of cv. 'Helfensteiner' (cross of cv. 'Pinot noir' and 'Schiava grossa') instead of a single 'Pinot noir'. These advances will help maintain the PN40024 genome as a gold-standard reference, also contributing towards the eventual elaboration of the grapevine pangenome.

The transcription factor VviNAC60 regulates senescence- and ripening-related processes in grapevine

Grapevine (Vitis vinifera L.) is one of the most widely cultivated fruit crops because the winemaking industry has huge economic relevance worldwide. Uncovering the molecular mechanisms controlling the developmental progression of plant organs will prove essential for maintaining high-quality grapes, expressly in the context of climate change, which impairs the ripening process. Through a deep inspection of transcriptomic data, we identified VviNAC60, a member of the NAC transcription factor family, as a putative regulator of grapevine organ maturation. We explored VviNAC60 binding landscapes through DNA affinity purification followed by sequencing and compared bound genes with transcriptomics datasets from grapevine plants stably and transiently overexpressing VviNAC60 to define a set of high-confidence targets. Among these, we identified key molecular markers associated with organ senescence and fruit ripening. Physiological, metabolic, and promoter activation analyses showed that VviNAC60 induces chlorophyll degradation and anthocyanin accumulation through the up-regulation of STAY-GREEN PROTEIN 1 (VviSGR1) and VviMYBA1, respectively, with the latter being up-regulated through a VviNAC60-VviNAC03 regulatory complex. Despite sharing a closer phylogenetic relationship with senescence-related homologues to the NAC transcription factor AtNAP, VviNAC60 complemented the non-ripening(nor) mutant phenotype in tomato (Solanum lycopersicum), suggesting a dual role as an orchestrator of both ripening- and senescence-related processes. Our data support VviNAC60 as a regulator of processes initiated in the grapevine vegetative- to mature-phase organ transition and therefore as a potential target for enhancing the environmental resilience of grapevine by fine-tuning the duration of the vegetative phase.

Peptidomics Methods Applied to the Study of Flower Development

Understanding the global and dynamic nature of plant developmental processes requires not only the study of the transcriptome, but also of the proteome, including its largely uncharacterized peptidome fraction. Recent advances in proteomics and high-throughput analyses of translating RNAs (ribosome profiling) have begun to address this issue, evidencing the existence of novel, uncharacterized, and possibly functional peptides. To validate the accumulation in tissues of sORF-encoded polypeptides (SEPs), the basic setup of proteomic analyses (i.e., LC-MS/MS) can be followed. However, the detection of peptides that are small (up to ~100 aa, 6-7 kDa) and novel (i.e., not annotated in reference databases) presents specific challenges that need to be addressed both experimentally and with computational biology resources. Several methods have been developed in recent years to isolate and identify peptides from plant tissues. In this chapter, we outline two different peptide extraction protocols and the subsequent peptide identification by mass spectrometry using the database search or the de novo identification methods.

Multi-Omics Methods Applied to Flower Development

Developmental processes in multicellular organisms depend on the proficiency of cells to orchestrate different gene expression programs. Over the past years, several studies of reproductive organ development have considered genomic analyses of transcription factors and global gene expression changes, modeling complex gene regulatory networks. Nevertheless, the dynamic view of developmental processes requires, as well, the study of the proteome in its expression, complexity, and relationship with the transcriptome. In this chapter, we describe a dual extraction method-for protein and RNA-for the characterization of genome expression at proteome level and its correlation to transcript expression data. We also present a shotgun proteomic procedure (LC-MS/MS) followed by a pipeline for the imputation of missing values in mass spectrometry results.

Transcriptomic and metabolomic integration as a resource in grapevine to study fruit metabolite quality traits

Transcriptomics and metabolomics are methodologies being increasingly chosen to perform molecular studies in grapevine (Vitis vinifera L.), focusing either on plant and fruit development or on interaction with abiotic or biotic factors. Currently, the integration of these approaches has become of utmost relevance when studying key plant physiological and metabolic processes. The results from these analyses can undoubtedly be incorporated in breeding programs whereby genes associated with better fruit quality (e.g., those enhancing the accumulation of health-promoting compounds) or with stress resistance (e.g., those regulating beneficial responses to environmental transition) can be used as selection markers in crop improvement programs. Despite the vast amount of data being generated, integrative transcriptome/metabolome meta-analyses (i.e., the joint analysis of several studies) have not yet been fully accomplished in this species, mainly due to particular specificities of metabolomic studies, such as differences in data acquisition (i.e., different compounds being investigated), unappropriated and unstandardized metadata, or simply no deposition of data in public repositories. These meta-analyses require a high computational capacity for data mining a priori, but they also need appropriate tools to explore and visualize the integrated results. This perspective article explores the universe of omics studies conducted in V. vinifera, focusing on fruit-transcriptome and metabolome analyses as leading approaches to understand berry physiology, secondary metabolism, and quality. Moreover, we show how omics data can be integrated in a simple format and offered to the research community as a web resource, giving the chance to inspect potential gene-to-gene and gene-to-metabolite relationships that can later be tested in hypothesis-driven research. In the frame of the activities promoted by the COST Action CA17111 INTEGRAPE, we present the first grapevine transcriptomic and metabolomic integrated database (TransMetaDb) developed within the Vitis Visualization (VitViz) platform (https://tomsbiolab.com/vitviz). This tool also enables the user to conduct and explore meta-analyses utilizing different experiments, therefore hopefully motivating the community to generate Findable, Accessible, Interoperable and Reusable (F.A.I.R.) data to be included in the future.

Direct regulation of shikimate, early phenylpropanoid, and stilbenoid pathways by Subgroup 2 R2R3-MYBs in grapevine

The stilbenoid pathway is responsible for the production of resveratrol in grapevine (Vitis vinifera L.). A few transcription factors (TFs) have been identified as regulators of this pathway but the extent of this control has not been deeply studied. Here we show how DNA affinity purification sequencing (DAP-Seq) allows for the genome-wide TF-binding site interrogation in grape. We obtained 5190 and 4443 binding events assigned to 4041 and 3626 genes for MYB14 and MYB15, respectively (approximately 40% of peaks located within −10 kb of transcription start sites). DAP-Seq of MYB14/MYB15 was combined with aggregate gene co-expression networks (GCNs) built from more than 1400 transcriptomic datasets from leaves, fruits, and flowers to narrow down bound genes to a set of high confidence targets. The analysis of MYB14, MYB15, and MYB13, a third uncharacterized member of Subgroup 2 (S2), showed that in addition to the few previously known stilbene synthase (STS) targets, these regulators bind to 30 of 47 STS family genes. Moreover, all three MYBs bind to several PAL, C4H, and 4CL genes, in addition to shikimate pathway genes, the WRKY03 stilbenoid co-regulator and resveratrol-modifying gene candidates among which ROMT2-3 were validated enzymatically. A high proportion of DAP-Seq bound genes were induced in the activated transcriptomes of transient MYB15-overexpressing grapevine leaves, validating our methodological approach for delimiting TF targets. Overall, Subgroup 2 R2R3-MYBs appear to play a key role in binding and directly regulating several primary and secondary metabolic steps leading to an increased flux towards stilbenoid production. The integration of DAP-Seq and reciprocal GCNs offers a rapid framework for gene function characterization using genome-wide approaches in the context of non-model plant species and stands up as a valid first approach for identifying gene regulatory networks of specialized metabolism.

A COMPASS for VESPUCCI: A FAIR Way to Explore the Grapevine Transcriptomic Landscape

Successfully integrating transcriptomic experiments is a challenging task with the ultimate goal of analyzing gene expression data in the broader context of all available measurements, all from a single point of access. In its second major release VESPUCCI, the integrated database of gene expression data for grapevine, has been updated to be FAIR-compliant, employing standards and created with open-source technologies. It includes all public grapevine gene expression experiments from both microarray and RNA-seq platforms. Transcriptomic data can be accessed in multiple ways through the newly developed COMPASS GraphQL interface, while the expression values are normalized using different methodologies to flexibly satisfy different analysis requirements. Sample annotations are manually curated and use standard formats and ontologies. The updated version of VESPUCCI provides easy querying and analyzing of integrated grapevine gene expression (meta)data and can be seamlessly embedded in any analysis workflow or tools. VESPUCCI is freely accessible and offers several ways of interaction, depending on the specific goals and purposes and/or user expertise; an overview can be found at https://vespucci.readthedocs.io/.

The Grape Gene Reference Catalogue as a Standard Resource for Gene Selection and Genetic Improvement

Effective crop improvement, whether through selective breeding or biotech strategies, is largely dependent on the cumulative knowledge of a species' pangenome and its containing genes. Acquiring this knowledge is specially challenging in grapevine, one of the oldest fruit crops grown worldwide, which is known to have more than 30,000 genes. Well-established research communities studying model organisms have created and maintained, through public and private funds, a diverse range of online tools and databases serving as repositories of genomes and gene function data. The lack of such resources for the non-model, but economically important, Vitis vinifera species has driven the need for a standardised collection of genes within the grapevine community. In an effort led by the Integrape COST Action CA17111, we have recently developed the first grape gene reference catalogue, where genes are ascribed to functional data, including their accession identifiers from different genome-annotation versions (https://integrape.eu/resources/genes-genomes/). We present and discuss this gene repository together with a validation-level scheme based on varied supporting evidence found in current literature. The catalogue structure and online submission form provided permits community curation. Finally, we present the Gene Cards tool, developed within the Vitis Visualization (VitViz) platform, to visualize the data collected in the catalogue and link gene function with tissue-specific expression derived from public transcriptomic data. This perspective article aims to present these resources to the community as well as highlight their potential use, in particular for plant-breeding applications.

Vitis OneGenE: A Causality-Based Approach to Generate Gene Networks in Vitis vinifera Sheds Light on the Laccase and Dirigent Gene Families

The abundance of transcriptomic data and the development of causal inference methods have paved the way for gene network analyses in grapevine. Vitis OneGenE is a transcriptomic data mining tool that finds direct correlations between genes, thus producing association networks. As a proof of concept, the stilbene synthase gene regulatory network obtained with OneGenE has been compared with published co-expression analysis and experimental data, including cistrome data for MYB stilbenoid regulators. As a case study, the two secondary metabolism pathways of stilbenoids and lignin synthesis were explored. Several isoforms of laccase, peroxidase, and dirigent protein genes, putatively involved in the final oxidative oligomerization steps, were identified as specifically belonging to either one of these pathways. Manual curation of the predicted sequences exploiting the last available genome assembly, and the integration of phylogenetic and OneGenE analyses, identified a group of laccases exclusively present in grapevine and related to stilbenoids. Here we show how network analysis by OneGenE can accelerate knowledge discovery by suggesting new candidates for functional characterization and application in breeding programs.

Metabolite analysis reveals distinct spatio-temporal accumulation of anthocyanins in two teinturier variants of cv. 'Gamay' grapevines (Vitis vinifera L.)

White-fleshed grape cv. 'Gamay' and its two teinturier variants presented distinct spatial-temporal accumulation of anthocyanins, with uncoupled accumulation of sugars and anthocyanins in 'Gamay Fréaux'. In most red grape cultivars, anthocyanins accumulate exclusively in the berry skin, while 'teinturier' cultivars also accumulate anthocyanins in the pulp. Here, we investigated the teinturier cvs. 'Gamay de Bouze' and 'Gamay Fréaux' (two somatic variants of the white-fleshed cv. 'Gamay') through metabolic and transcript analysis to clarify whether these two somatic variants have the same anthocyanin accumulation pattern in the skin and pulp, and whether primary metabolites are also affected. The skin of the three cultivars and the pulp of 'Gamay de Bouze' begun to accumulate anthocyanins at the onset of berry ripening. However, the pulp of 'Gamay Fréaux' exhibited a distinct anthocyanin accumulation pattern, starting as early as fruit set with very low level of sugars. The highest level of anthocyanins was found in 'Gamay Fréaux' skin, followed by 'Gamay de Bouze' and 'Gamay'. Consistently, the transcript abundance of genes involved in anthocyanin biosynthesis were in line with the anthocyanin levels in the three cultivars. Despite no evident differences in pulp sugar content, the concentration of glucose and fructose in the skin of 'Gamay Fréaux' was only half of those in the skin of 'Gamay' and 'Gamay de Bouze' throughout all berry ripening, suggesting an uncoupled accumulation of sugars and anthocyanins in 'Gamay Fréaux'. The study provides a comprehensive view of metabolic consequences in grape somatic variants and the three almost isogenic genotypes can serve as ideal reagents to further uncover the mechanisms underlying the linkage between sugar and anthocyanin accumulation.

Identification of ABA-Mediated Genetic and Metabolic Responses to Soil Flooding in Tomato (Solanum lycopersicum L. Mill)

Soil flooding is a compound abiotic stress that alters soil properties and limits atmospheric gas diffusion (O₂ and CO₂) to the roots. The involvement of abscisic acid (ABA) in the regulation of soil flooding-specific genetic and metabolic responses has been scarcely studied despite its key importance as regulator in other abiotic stress conditions. To attain this objective, wild type and ABA-deficient tomatoes were subjected to short-term (24 h) soil waterlogging. After this period, gas exchange parameters were reduced in the wild type but not in ABA-deficient plants that always had higher E and g _s . Transcript and metabolite alterations were more intense in waterlogged tissues, with genotype-specific variations. Waterlogging reduced the ABA levels in the roots while inducing PYR/PYL/RCAR ABA receptors and ABA-dependent transcription factor transcripts, of which induction was less pronounced in the ABA-deficient genotype. Ethylene/O₂-dependent genetic responses (ERFVIIs, plant anoxia survival responses, and genes involved in the N-degron pathway) were induced in hypoxic tissues independently of the genotype. Interestingly, genes encoding a nitrate reductase and a phytoglobin involved in NO biosynthesis and scavenging and ERFVII stability were induced in waterlogged tissues, but to a lower extent in ABA-deficient tomato. At the metabolic level, flooding-induced accumulation of Ala was enhanced in ABA-deficient lines following a differential accumulation of Glu and Asp in both hypoxic and aerated tissues, supporting their involvement as sources of oxalacetate to feed the tricarboxylic acid cycle in waterlogged tissues and constituting a potential advantage upon long periods of soil waterlogging. The promoter analysis of upregulated genes indicated that the production of oxalacetate from Asp via Asp oxidase, energy processes such as acetyl-CoA, ATP, and starch biosynthesis, and the lignification process were likely subjected to ABA regulation. Taken together, these data indicate that ABA depletion in waterlogged tissues acts as a positive signal, inducing several specific genetic and metabolic responses to soil flooding.

Salinity impairs photosynthetic capacity and enhances carotenoid-related gene expression and biosynthesis in tomato (Solanum lycopersicum L. cv. Micro-Tom)

Carotenoids are essential components of the photosynthetic antenna and reaction center complexes, being also responsible for antioxidant defense, coloration, and many other functions in multiple plant tissues. In tomato, salinity negatively affects the development of vegetative organs and productivity, but according to previous studies it might also increase fruit color and taste, improving its quality, which is a current agricultural challenge. The fruit quality parameters that are increased by salinity are cultivar-specific and include carotenoid, sugar, and organic acid contents. However, the relationship between vegetative and reproductive organs and response to salinity is still poorly understood. Considering this, Solanum lycopersicum cv. Micro-Tom plants were grown in the absence of salt supplementation as well as with increasing concentrations of NaCl for 14 weeks, evaluating plant performance from vegetative to reproductive stages. In response to salinity, plants showed a significant reduction in net photosynthesis, stomatal conductance, PSII quantum yield, and electron transport rate, in addition to an increase in non-photochemical quenching. In line with these responses the number of tomato clusters decreased, and smaller fruits with higher soluble solids content were obtained. Mature-green fruits also displayed a salt-dependent higher induction in the expression of PSY1, PDS, ZDS, and LYCB, key genes of the carotenoid biosynthesis pathway, in correlation with increased lycopene, lutein, β-carotene, and violaxanthin levels. These results suggest a key relationship between photosynthetic plant response and yield, involving impaired photosynthetic capacity, increased carotenoid-related gene expression, and carotenoid biosynthesis.

Synthetic conversion of leaf chloroplasts into carotenoid-rich plastids reveals mechanistic basis of natural chromoplast development

Plastids, the defining organelles of plant cells, undergo physiological and morphological changes to fulfill distinct biological functions. In particular, the differentiation of chloroplasts into chromoplasts results in an enhanced storage capacity for carotenoids with industrial and nutritional value such as beta-carotene (provitamin A). Here, we show that synthetically inducing a burst in the production of phytoene, the first committed intermediate of the carotenoid pathway, elicits an artificial chloroplast-to-chromoplast differentiation in leaves. Phytoene overproduction initially interferes with photosynthesis, acting as a metabolic threshold switch mechanism that weakens chloroplast identity. In a second stage, phytoene conversion into downstream carotenoids is required for the differentiation of chromoplasts, a process that involves a concurrent reprogramming of nuclear gene expression and plastid morphology for improved carotenoid storage. We hence demonstrate that loss of photosynthetic competence and enhanced production of carotenoids are not just consequences but requirements for chloroplasts to differentiate into chromoplasts.

Comprehending and improving cannabis specialized metabolism in the systems biology era

Cannabis sativa is a source of food, fiber and specialized metabolites such as cannabinoids, with psychoactive and pharmacological effects. Due to its expanding and increasingly-accepted use in medicine, cannabis cultivation is acquiring more importance and less social stigma. Humans initiated different domestication episodes whose later spread gave rise to a plethora of landrace cultivars. At present, breeders cross germplasms from different gene pools depending on their specific use. The fiber (hemp) and drug (marijuana) types of C. sativa differ in their cannabinoid chemical composition phenotype (chemotype) and also in the accumulation of terpenoid compounds that constitute a strain's particular flavor and scent. Cannabinoids are isoprenylated polyketides among which cannabidiolic acid (CBDA) and (-)-trans-Δ⁹-tetrahydrocannabinol acid (THCA) have been well-documented for their many effects on humans. Here, we review the most studied specialized metabolic pathways in C. sativa, showing how terpenes and cannabinoids share both part of the isoprenoid pathway and the same biosynthetic compartmentalization (i.e. glandular trichomes of leaves and flowers). We enlist the several studies that have deciphered these pathways in this species including physical and genetic maps, QTL analyses and localization and enzymatic studies of cannabinoid and terpene synthases. In addition, new comparative modeling of cannabinoid synthases and phylogenetic trees are presented. We describe the genome sequencing initiatives of several accessions with the concomitant generation of next-generation genome maps and transcriptomic data. Very recently, proteomic characterizations and systems biology approaches such as those applying network theory or the integration of multi-omics data have increased the knowledge on gene function, enzyme diversity and metabolite content in C. sativa. In this revision we drift through the history, present and future of cannabis research and on how second- and third-generation sequencing technologies are bringing light to the field of cannabis specialized metabolism. We also discuss different biotechnological approaches for producing cannabinoids in engineered microorganisms.

Genetic analysis of a white-to-red berry skin color reversion and its transcriptomic and metabolic consequences in grapevine (Vitis vinifera cv. 'Moscatel Galego')

BACKGROUND

Somatic mutations occurring within meristems of vegetative propagation material have had a major role in increasing the genetic diversity of the domesticated grapevine (Vitis vinifera subsp. vinifera). The most well studied somatic variation in this species is the one affecting fruit pigmentation, leading to a plethora of different berry skin colors. Color depletion and reversion are often observed in the field. In this study we analyzed the origin of a novel white-to-red skin color reversion and studied its possible metabolic and transcriptomic consequences on cv. 'Muscat à Petits Grains Blancs' (synonym cv. 'Moscatel Galego Branco'), a member of the large family of Muscats.

RESULTS

The mild red-skinned variant (cv. 'Muscat à Petits Grains Rouge', synonym cv. 'Moscatel Galego Roxo'), characterized by a preferential accumulation of di-hydroxylated anthocyanins, showed in heterozygosis a partially-excised Gret1 retrotransposon in the promoter region of the MYBA1 anthocyanin regulator, while MYBA2 was still in homozygosis for its non-functional allele. Through metabolic (anthocyanin, resveratrol and piceid quantifications) and transcriptomic (RNA-Seq) analyses, we show that within a near-isogenic background, the transcriptomic consequences of color reversion are largely associated to diminished light/UV-B responses probably as a consequence of the augment of metabolic sunscreens (i.e. anthocyanins).

CONCLUSIONS

We propose that the reduced activity of the flavonoid tri-hydroxylated sub-branch and decreased anthocyanin synthesis and modification (e.g. methylation and acylation) are the potential causes for the mild red-skinned coloration in the pigmented revertant. The observed positive relation between anthocyanins and stilbenes could be attributable to an increased influx of phenylpropanoid intermediaries due to the replenished activity of MYBA1, an effect yet to be demonstrated in other somatic variants.

Combinatorial Regulation of Stilbene Synthase Genes by WRKY and MYB Transcription Factors in Grapevine (Vitis vinifera L.)

Stilbene synthase (STS) is the key enzyme leading to the biosynthesis of resveratrol. Recently we reported two R2R3-MYB transcription factor (TF) genes that regulate the stilbene biosynthetic pathway in grapevine: VviMYB14 and VviMYB15. These genes are strongly co-expressed with STS genes under a range of stress and developmental conditions, in agreement with the specific activation of STS promoters by these TFs. Genome-wide gene co-expression analysis using two separate transcriptome compendia based on microarray and RNA sequencing data revealed that WRKY TFs were the top TF family correlated with STS genes. On the basis of correlation frequency, four WRKY genes, namely VviWRKY03, VviWRKY24, VviWRKY43 and VviWRKY53, were further shortlisted and functionally validated. Expression analyses under both unstressed and stressed conditions, together with promoter-luciferase reporter assays, suggested different hierarchies for these TFs in the regulation of the stilbene biosynthetic pathway. In particular, VviWRKY24 seems to act as a singular effector in the activation of the VviSTS29 promoter, while VviWRKY03 acts through a combinatorial effect with VviMYB14, suggesting that these two regulators may interact at the protein level as previously reported in other species.

The GARP/MYB-related grape transcription factor AQUILO improves cold tolerance and promotes the accumulation of raffinose family oligosaccharides

Grapevine (Vitis vinifera L.) is a widely cultivated fruit crop whose growth and productivity are greatly affected by low temperatures. On the other hand, wild Vitis species represent valuable genetic resources of natural stress tolerance. We have isolated and characterized a MYB-like gene encoding a putative GARP-type transcription factor from Amur grape (V. amurensis) designated as VaAQUILO. AQUILO (AQ) is induced by cold in both V. amurensis and V. vinifera, and its overexpression results in significantly improved tolerance to cold both in transgenic Arabidopsis and in Amur grape calli. In Arabidopsis, the ectopic expression of VaAQ increased antioxidant enzyme activities and up-regulated reactive oxygen species- (ROS) scavenging-related genes. Comparative mRNA sequencing profiling of 35S:VaAQ Arabidopsis plants suggests that this transcription factor is related to phosphate homeostasis like their Arabidopsis closest homologues: AtHRS1 and AtHHO2. However, when a cold stress is imposed, AQ is tightly associated with the cold-responsive pathway and with the raffinose family oligosaccharides (RFOs), as observed by the up-regulation of galactinol synthase (GoLS) and raffinose synthase genes. Gene co-expression network (GCN) and cis-regulatory element (CRE) analyses in grapevine indicated AQ as potentially regulating VvGoLS genes. Increased RFO content was confirmed in both transgenic Arabidopsis and Amur grape calli overexpressing VaAQ. Taken together, our results imply that AQ improves cold tolerance through promoting the accumulation of osmoprotectants.

A group of grapevine MYBA transcription factors located in chromosome 14 control anthocyanin synthesis in vegetative organs with different specificities compared with the berry color locus

Grapevine organs accumulate anthocyanins in a cultivar-specific and environmentally induced manner. The MYBA1-A2 genes within the berry color locus in chromosome 2 represent the major genetic determinants of fruit color. The simultaneous occurrence of transposon insertions and point mutations in these genes is responsible for most white-skinned phenotypes; however, the red pigmentation found in vegetative organs suggests the presence of additional regulators. This work describes a genomic region of chromosome 14 containing three closely related R2R3-MYB genes, named MYBA5, MYBA6 and MYBA7. Ectopic expression of the latter two genes in grapevine hairy roots promoted anthocyanin accumulation without affecting other phenylpropanoids. Transcriptomic profiling of hairy roots expressing MYBA1, MYBA6 and MYBA7 showed that these regulators share the activation of late biosynthetic and modification/transport-related genes, but differ in the activation of the FLAVONOID-3'5'-HYDROXYLASE (F3'5'H) family. An alternatively spliced MYBA6 variant was incapable of activating anthocyanin synthesis, however, because of the lack of an MYC1 interaction domain. MYBA1, MYBA6.1 and MYBA7 activated the promoters of UDP-GLUCOSE:FLAVONOID 3-O-GLUCOSYLTRANSFERASE (UFGT) and ANTHOCYANIN 3-O-GLUCOSIDE-6″-O-ACYLTRANSFERASE (3AT), but only MYBA1 induced F3'5'H in concordance with the low proportion of tri-hydroxylated anthocyanins found in MYBA6-A7 hairy roots. This putative new color locus is related to the red/cyanidic pigmentation of vegetative organs in black- and white-skinned cultivars, and forms part of the UV-B radiation response pathway orchestrated by ELONGATED HYPOCOTYL 5 (HY5). These results demonstrate the involvement of additional anthocyanin regulators in grapevine and suggest an evolutionary divergence between the two grape color loci for controlling additional targets of the flavonoid pathway.

Constructing Integrated Networks for Identifying New Secondary Metabolic Pathway Regulators in Grapevine: Recent Applications and Future Opportunities

Representing large biological data as networks is becoming increasingly adopted for predicting gene function while elucidating the multifaceted organization of life processes. In grapevine (Vitis vinifera L.), network analyses have been mostly adopted to contribute to the understanding of the regulatory mechanisms that control berry composition. Whereas, some studies have used gene co-expression networks to find common pathways and putative targets for transcription factors related to development and metabolism, others have defined networks of primary and secondary metabolites for characterizing the main metabolic differences between cultivars throughout fruit ripening. Lately, proteomic-related networks and those integrating genome-wide analyses of promoter regulatory elements have also been generated. The integration of all these data in multilayered networks allows building complex maps of molecular regulation and interaction. This perspective article describes the currently available network data and related resources for grapevine. With the aim of illustrating data integration approaches into network construction and analysis in grapevine, we searched for berry-specific regulators of the phenylpropanoid pathway. We generated a composite network consisting of overlaying maps of co-expression between structural and transcription factor genes, integrated with the presence of promoter cis-binding elements, microRNAs, and long non-coding RNAs (lncRNA). This approach revealed new uncharacterized transcription factors together with several microRNAs potentially regulating different steps of the phenylpropanoid pathway, and one particular lncRNA compromising the expression of nine stilbene synthase (STS) genes located in chromosome 10. Application of network-based approaches into multi-omics data will continue providing supplementary resources to address important questions regarding grapevine fruit quality and composition.

Transcriptome-Wide Identification of Novel UV-B- and Light Modulated Flavonol Pathway Genes Controlled by VviMYBF1

Flavonols constitute a group of flavonoids with important photoprotective roles in plants. In addition, flavonol content and composition greatly influences fruit quality. We previously demonstrated that the grapevine R2R3-MYB transcription factor (TF) VviMYBF1 promotes flavonol accumulation by inducing the expression of flavonol synthase (VviFLS1/VviFLS4), a key step of the initial flavonol pathway. Despite this, gene networks underlying flavonol modification in grapevine including both structural and regulatory genes remain poorly understood. In order to identify flavonol modifying genes and TFs acting downstream of VviMYBF1 a microarray-based transcriptome analysis was performed on grapevine hairy roots ectopically expressing VviMYBF1 or a Green Fluorescent Protein as control. VviFLS1 was induced in VviMYBF1 transgenic roots and glycosylated flavonols accumulated significantly compared with control lines. Among the differentially expressed genes, potential flavonol-modifying enzymes with predicted rhamnosyltransferase (e.g., RhaT1) or glycosyltransferase (e.g., GT3) activities were identified. In addition, important TFs of the MYB and bZIP families such as the proanthocyanidin regulator VviMYBPA1 and the UV-B light responsive HY5 homolog VviHYH were significantly altered in their expression pattern by overexpression of VviMYBF1. Co-temporal expression analysis demonstrated positive correlation of VviMYBF1 with VviFLS1, VviGT3, and VviRhaT1 during berry development and in fruits ripened with different light and UV-B radiation conditions at field. These results show that VviMYBF1 overexpression led to the identification of novel genes of the flavonol pathway and that the flavonol modifying machinery can be influenced by agricultural practices to optimize flavonol composition in grapes.

The photomorphogenic factors UV-B RECEPTOR 1, ELONGATED HYPOCOTYL 5, and HY5 HOMOLOGUE are part of the UV-B signalling pathway in grapevine and mediate flavonol accumulation in response to the environment

Grapevine (Vitis vinifera L.) is a species well known for its adaptation to radiation. However, photomorphogenic factors related to UV-B responses have not been molecularly characterized. We cloned and studied the role of UV-B RECEPTOR (UVR1), ELONGATED HYPOCOTYL 5 (HY5), and HY5 HOMOLOGUE (HYH) from V. vinifera We performed gene functional characterizations, generated co-expression networks, and tested them in different environmental conditions. These genes complemented the Arabidopsis uvr8 and hy5 mutants in morphological and secondary metabolic responses to radiation. We combined microarray and RNA sequencing (RNA-seq) data with promoter inspections to identify HY5 and HYH putative target genes and their DNA binding preferences. Despite sharing a large set of common co-expressed genes, we found different hierarchies for HY5 and HYH depending on the organ and stress condition, reflecting both co-operative and partially redundant roles. New candidate UV-B gene markers were supported by the presence of HY5-binding sites. These included a set of flavonol-related genes that were up-regulated in a HY5 transient expression assay. We irradiated in vitro plantlets and fruits from old potted vines with high and low UV-B exposures and followed the accumulation of flavonols and changes in gene expression in comparison with non-irradiated conditions. UVR1, HY5, and HYH expression varied with organ, developmental stage, and type of radiation. Surprisingly, UVR1 expression was modulated by shading and temperature in berries, but not by UV-B radiation. We propose that the UV-B response machinery favours berry flavonol accumulation through the activation of HY5 and HYH at different developmental stages at both high and low UV-B exposures.

Transcriptomic and Metabolomic Networks in the Grape Berry Illustrate That it Takes More Than Flavonoids to Fight Against Ultraviolet Radiation

Plants are constantly challenged by environmental fluctuations. In response, they have developed a wide range of morphological and biochemical adaptations committed to ameliorate the effects of abiotic stress. When exposed to higher solar radiation levels, plants activate the synthesis of a large set of enzymes and secondary metabolites as part of a complex sunscreen and antioxidant defense mechanism. Grapevine (Vitis vinifera L.) has become a widely used system for studying adaptive responses to this type of stress since changes in berry composition, positively influenced by increased ultraviolet (UV) radiation levels, improve the quality of wines subsequently produced. Despite the fact that most of the attention has been directed toward the synthesis of flavonoids, recent transcriptomic and metabolomic studies have shown that stilbenoids and isoprenoids (e.g., terpenes and carotenoids) are also an important part of the grape UV-response machinery. This minireview focuses on the latest findings referring to the metabolic responses of grapes to UV radiation and proposes a model for its transcriptional control. Depending on the berry developmental stage and the type of radiation (i.e., irradiance level, exposure length), increased UV levels activate different metabolic pathways through the activity of master regulators belonging to the basic Leucine Zipper Domain (bZIP) and R2R3-MYB transcription factor families. This transcriptional control is influenced by the interaction of other environmental factors such as light, temperature or soil water availability. In grapevine, phenylpropanoids are part of, but are not the whole story, in the fight against radiation damage.

A systems-oriented analysis of the grapevine R2R3-MYB transcription factor family uncovers new insights into the regulation of stilbene accumulation

R2R3-MYB transcription factors (TFs) belong to a large and functionally diverse protein superfamily in plants. In this study, we explore the evolution and function of this family in grapevine (Vitis vinifera L.), a high-value fruit crop. We identified and manually curated 134 genes using RNA-Seq data, and named them systematically according to the Super-Nomenclature Committee. We identified novel genes, splicing variants and grapevine/woody-specific duplicated subgroups, suggesting possible neo- and sub-functionalization events. Regulatory network analysis ascribed biological functions to uncharacterized genes and validated those of known genes (e.g. secondary cell wall biogenesis and flavonoid biosynthesis). A comprehensive analysis of different MYB binding motifs in the promoters of co-expressed genes predicted grape R2R3-MYB binding preferences and supported evidence for putative downstream targets. Enrichment of cis-regulatory motifs for diverse TFs reinforced the notion of transcriptional coordination and interaction between MYBs and other regulators. Analysis of the network of Subgroup 2 showed that the resveratrol-related VviMYB14 and VviMYB15 share common co-expressed STILBENE SYNTHASE genes with the uncharacterized VviMYB13. These regulators have distinct expression patterns within organs and in response to biotic and abiotic stresses, suggesting a pivotal role of VviMYB13 in regulating stilbene accumulation in vegetative tissues and under biotic stress conditions.

Genome-wide analyses for dissecting gene regulatory networks in the shoot apical meristem

Shoot apical meristem activity is controlled by complex regulatory networks in which components such as transcription factors, miRNAs, small peptides, hormones, enzymes and epigenetic marks all participate. Many key genes that determine the inherent characteristics of the shoot apical meristem have been identified through genetic approaches. Recent advances in genome-wide studies generating extensive transcriptomic and DNA-binding datasets have increased our understanding of the interactions within the regulatory networks that control the activity of the meristem, identifying new regulators and uncovering connections between previously unlinked network components. In this review, we focus on recent studies that illustrate the contribution of whole genome analyses to understand meristem function.

The phenylpropanoid pathway is controlled at different branches by a set of R2R3-MYB C2 repressors in grapevine

Because of the vast range of functions that phenylpropanoids possess, their synthesis requires precise spatiotemporal coordination throughout plant development and in response to the environment. The accumulation of these secondary metabolites is transcriptionally controlled by positive and negative regulators from the MYB and basic helix-loop-helix protein families. We characterized four grapevine (Vitis vinifera) R2R3-MYB proteins from the C2 repressor motif clade, all of which harbor the ethylene response factor-associated amphiphilic repression domain but differ in the presence of an additional TLLLFR repression motif found in the strong flavonoid repressor Arabidopsis (Arabidopsis thaliana) AtMYBL2. Constitutive expression of VvMYB4a and VvMYB4b in petunia (Petunia hybrida) repressed general phenylpropanoid biosynthetic genes and selectively reduced the amount of small-weight phenolic compounds. Conversely, transgenic petunia lines expressing VvMYBC2-L1 and VvMYBC2-L3 showed a severe reduction in petal anthocyanins and seed proanthocyanidins together with a higher pH of crude petal extracts. The distinct function of these regulators was further confirmed by transient expression in tobacco (Nicotiana benthamiana) leaves and grapevine plantlets. Finally, VvMYBC2-L3 was ectopically expressed in grapevine hairy roots, showing a reduction in proanthocyanidin content together with the down-regulation of structural and regulatory genes of the flavonoid pathway as revealed by a transcriptomic analysis. The physiological role of these repressors was inferred by combining the results of the functional analyses and their expression patterns in grapevine during development and in response to ultraviolet B radiation. Our results indicate that VvMYB4a and VvMYB4b may play a key role in negatively regulating the synthesis of small-weight phenolic compounds, whereas VvMYBC2-L1 and VvMYBC2-L3 may additionally fine tune flavonoid levels, balancing the inductive effects of transcriptional activators.

Inspection of the grapevine BURP superfamily highlights an expansion of RD22 genes with distinctive expression features in berry development and ABA-mediated stress responses

The RESPONSIVE TO DEHYDRATION 22 (RD22) gene is a molecular link between abscisic acid (ABA) signalling and abiotic stress responses. Its expression has been used as a reliable ABA early response marker. In Arabidopsis, the single copy RD22 gene possesses a BURP domain also located at the C-terminus of USP embryonic proteins and the beta subunit of polygalacturonases. In grapevine, a RD22 gene has been identified but putative paralogs are also found in the grape genome, possibly forming a large RD22 family in this species. In this work, we searched for annotations containing BURP domains in the Vitis vinifera genome. Nineteen proteins were defined by a comparative analysis between the two genome predictions and RNA-Seq data. These sequences were compared to other plant BURPs identified in previous genome surveys allowing us to reconceive group classifications based on phylogenetic relationships and protein motif occurrence. We observed a lineage-specific evolution of the RD22 family, with the biggest expansion in grapevine and poplar. In contrast, rice, sorghum and maize presented highly expanded monocot-specific groups. The Vitis RD22 group may have expanded from segmental duplications as most of its members are confined to a region in chromosome 4. The inspection of transcriptomic data revealed variable expression of BURP genes in vegetative and reproductive organs. Many genes were induced in specific tissues or by abiotic and biotic stresses. Three RD22 genes were further studied showing that they responded oppositely to ABA and to stress conditions. Our results show that the inclusion of RNA-Seq data is essential while describing gene families and improving gene annotations. Robust phylogenetic analyses including all BURP members from other sequenced species helped us redefine previous relationships that were erroneously established. This work provides additional evidence for RD22 genes serving as marker genes for different organs or stresses in grapevine.

Dynamics of chromatin accessibility and gene regulation by MADS-domain transcription factors in flower development

Background

Development of eukaryotic organisms is controlled by transcription factors that trigger specific and global changes in gene expression programs. In plants, MADS-domain transcription factors act as master regulators of developmental switches and organ specification. However, the mechanisms by which these factors dynamically regulate the expression of their target genes at different developmental stages are still poorly understood.

Results

We characterized the relationship of chromatin accessibility, gene expression, and DNA binding of two MADS-domain proteins at different stages of Arabidopsis flower development. Dynamic changes in APETALA1 and SEPALLATA3 DNA binding correlated with changes in gene expression, and many of the target genes could be associated with the developmental stage in which they are transcriptionally controlled. We also observe dynamic changes in chromatin accessibility during flower development. Remarkably, DNA binding of APETALA1 and SEPALLATA3 is largely independent of the accessibility status of their binding regions and it can precede increases in DNA accessibility. These results suggest that APETALA1 and SEPALLATA3 may modulate chromatin accessibility, thereby facilitating access of other transcriptional regulators to their target genes.

Conclusions

Our findings indicate that different homeotic factors regulate partly overlapping, yet also distinctive sets of target genes in a partly stage-specific fashion. By combining the information from DNA-binding and gene expression data, we are able to propose models of stage-specific regulatory interactions, thereby addressing dynamics of regulatory networks throughout flower development. Furthermore, MADS-domain TFs may regulate gene expression by alternative strategies, one of which is modulation of chromatin accessibility.

Identification of Arabidopsis knockout lines for genes of interest

Determining gene function through reverse genetics has been an important experimental approach in the field of flower development. The method largely relies on the availability of knockout lines for the gene of interest. Insertional mutagenesis can be performed using either T-DNA or transposable elements, but the former has been more frequently employed in Arabidopsis. A primary concern for working with insertional mutant lines is whether the respective insertion results in a complete or rather a partial loss of gene function. The effect of the insertion largely depends on its position with respect to the structure of the gene. In order to quickly identify and obtain knockout lines for genes of interest in Arabidopsis, more than 325,000 mapped insertion lines have been catalogued on indexed libraries and made publicly available to researchers. Online accessible databases provide information regarding the site of insertion, whether a mutant line is available in a homozygous or hemizygous state, and outline technical aspects for plant identification, such as primer design tools used for genotyping. In this chapter, we describe the procedure for isolating knockout lines for genes of interest in Arabidopsis.

The grapevine guard cell-related VvMYB60 transcription factor is involved in the regulation of stomatal activity and is differentially expressed in response to ABA and osmotic stress

BACKGROUND

Under drought, plants accumulate the signaling hormone abscisic acid (ABA), which induces the rapid closure of stomatal pores to prevent water loss. This event is trigged by a series of signals produced inside guard cells which finally reduce their turgor. Many of these events are tightly regulated at the transcriptional level, including the control exerted by MYB proteins. In a previous study, while identifying the grapevine R2R3 MYB family, two closely related genes, VvMYB30 and VvMYB60 were found with high similarity to AtMYB60, an Arabidopsis guard cell-related drought responsive gene.

RESULTS

Promoter-GUS transcriptional fusion assays showed that expression of VvMYB60 was restricted to stomatal guard cells and was attenuated in response to ABA. Unlike VvMYB30, VvMYB60 was able to complement the loss-of-function atmyb60-1 mutant, indicating that VvMYB60 is the only true ortholog of AtMYB60 in the grape genome. In addition, VvMYB60 was differentially regulated during development of grape organs and in response to ABA and drought-related stress conditions.

CONCLUSIONS

These results show that VvMYB60 modulates physiological responses in guard cells, leading to the possibility of engineering stomatal conductance in grapevine, reducing water loss and helping this species to tolerate drought under extreme climatic conditions.

Effect of pollination and fertilization on the expression of genes related to floral transition, hormone synthesis and berry development in grapevine

In the present work, the effect of assisted fertilization on anatomical, morphological and gene expression changes occurring in carpels and during early stages of berry development in Vitis vinifera were studied. Inflorescences were emasculated before capfall, immediately manually pollinated (EP) and fruit development was compared to emasculated but non-pollinated (ENP) and self-pollinated inflorescences (NESP). The diameter of berries derived from pollinated flowers (EP and NESP) was significantly higher than from non-pollinated flowers (ENP) at 21 days after emasculation/pollination (DAE), and a rapid increase in the size of the inner mesocarp, together with the presence of an embryo-like structure, were observed. The expression of gibberellin oxidases (GA20ox and GA2ox), anthranilate synthase (related to auxin synthesis) and cytokinin synthase coding genes was studied to assess the relationship between hormone synthesis and early berry development, while flower patterning genes were analyzed to describe floral transition. Significant expression changes were found for hormone-related genes, suggesting that their expression at early stages of berry development (13 DAE) is related to cell division and differentiation of mesocarp tissue at a later stage (21 DAE). Expression of hormone-related genes also correlates with the expression of VvHB13, a gene related to mesocarp expansion, and with an increased repression of floral patterning genes (PISTILLATA and TM6), which may contribute to prevent floral transition inhibiting fruit growth before fertilization takes place.

Arabidopsis paves the way: genomic and network analyses in crops

Arabidopsis genomic and network analyses have facilitated crop research towards the understanding of many biological processes of fundamental importance for agriculture. Genes that were identified through genomic analyses in Arabidopsis have been used to manipulate crop traits such as pathogen resistance, yield, water-use efficiency, and drought tolerance, with the effects being tested in field conditions. The integration of diverse Arabidopsis genome-wide datasets in probabilistic functional networks has been demonstrated as a feasible strategy to associate novel genes with traits of interest, and novel genomic methods continue to be developed. The combination of genome-wide location studies, using ChIP-Seq, with gene expression profiling data is affording a genome-wide view of regulatory networks previously delineated through genetic and molecular analyses, leading to the identification of novel components and of new connections within these networks.

Isolation of WDR and bHLH genes related to flavonoid synthesis in grapevine (Vitis vinifera L.)

Anthocyanins and tannins are two of the most abundant flavonoids found in grapevine, and their synthesis is derived from the phenylpropanoid pathway. As described for model species such as Arabidopsis thaliana, maize and petunia, the end-point branches of this pathway are tightly regulated by the combinatorial interaction of three families of regulatory factors; MYB, bHLH (also known as MYC) and WDR proteins. Among these, only MYB genes have been previously identified in grapes. Here, we report the isolation of the first members from the WDR and bHLH families found in Vitis vinifera, named WDR1, WDR2 and MYCA1. WDR1 contributed positively to the accumulation of anthocyanins when it was overexpressed in A. thaliana, although it was not possible to determine the function of WDR2 by ectopic expression. The sub-cellular localizations of WDR1 and MYCA1 were observed by means of GFP-fusion proteins, indicating both cytoplasm and nuclear localization, in contrast to the localization of a MYB factor exclusively in the nucleus. The expression patterns of these genes were quantified in coloured reproductive organs throughout development, and correlated with anthocyanin accumulation and the expression profiles of the flavonoid-related MYBA1-2, UFGT, and ANR genes. In vitro grapevine plantlets grown under high salt concentrations showed a cultivar-dependent response for anthocyanin accumulation, which correlated with the expression of MYBA1-2, MYCA1 and WDR1 genes. These results suggest that MYCA1 may regulate ANR and UFGT and that this last control is easier to distinguish whenever MYBA genes are absent or in low abundance. Future studies should address the specific interactions of these proteins and their quantitative contribution to flavonoid synthesis in grape berries.

Isolation of WDR and bHLH genes related to flavonoid synthesis in grapevine (Vitis vinifera L.)

Anthocyanins and tannins are two of the most abundant flavonoids found in grapevine, and their synthesis is derived from the phenylpropanoid pathway. As described for model species such as Arabidopsis thaliana, maize and petunia, the end-point branches of this pathway are tightly regulated by the combinatorial interaction of three families of regulatory factors; MYB, bHLH (also known as MYC) and WDR proteins. Among these, only MYB genes have been previously identified in grapes. Here, we report the isolation of the first members from the WDR and bHLH families found in Vitis vinifera, named WDR1, WDR2 and MYCA1. WDR1 contributed positively to the accumulation of anthocyanins when it was overexpressed in A. thaliana, although it was not possible to determine the function of WDR2 by ectopic expression. The sub-cellular localizations of WDR1 and MYCA1 were observed by means of GFP-fusion proteins, indicating both cytoplasm and nuclear localization, in contrast to the localization of a MYB factor exclusively in the nucleus. The expression patterns of these genes were quantified in coloured reproductive organs throughout development, and correlated with anthocyanin accumulation and the expression profiles of the flavonoid-related MYBA1-2, UFGT, and ANR genes. In vitro grapevine plantlets grown under high salt concentrations showed a cultivar-dependent response for anthocyanin accumulation, which correlated with the expression of MYBA1-2, MYCA1 and WDR1 genes. These results suggest that MYCA1 may regulate ANR and UFGT and that this last control is easier to distinguish whenever MYBA genes are absent or in low abundance. Future studies should address the specific interactions of these proteins and their quantitative contribution to flavonoid synthesis in grape berries.

Post-veraison sunlight exposure induces MYB-mediated transcriptional regulation of anthocyanin and flavonol synthesis in berry skins of Vitis vinifera

Anthocyanins, flavan-3-ols, and flavonols are the three major classes of flavonoid compounds found in grape berry tissues. Several viticultural practices increase flavonoid content in the fruit, but the underlying genetic mechanisms responsible for these changes have not been completely deciphered. The impact of post-veraison sunlight exposure on anthocyanin and flavonol accumulation in grape berry skin and its relation to the expression of different transcriptional regulators known to be involved in flavonoid synthesis was studied. Treatments consisting of removing or moving aside the basal leaves which shade berry clusters were applied. Shading did not affect sugar accumulation or gene expression of HEXOSE TRANSPORTER 1, although in the leaf removal treatment, these events were retarded during the first weeks of ripening. Flavonols were the most drastically reduced flavonoids following shading and leaf removal treatments, related to the reduced expression of FLAVONOL SYNTHASE 4 and its putative transcriptional regulator MYB12. Anthocyanin accumulation and the expression of CHS2, LDOX, OMT, UFGT, MYBA1, and MYB5a genes were also affected. Other regulatory genes were less affected or not affected at all by these treatments. Non-transcriptional control mechanisms for flavonoid synthesis are also suggested, especially during the initial stages of ripening. Although berries from the leaf removal treatment received more light than shaded fruits, malvidin-3-glucoside and total flavonol content was reduced compared with the treatment without leaf removal. This work reveals that flavonol-related gene expression responds rapidly to field changes in light levels, as shown by the treatment in which shaded fruits were exposed to light in the late stages of ripening. Taken together, this study establishes MYB-specific responsiveness for the effect of sun exposure and sugar transport on flavonoid synthesis.

Genetic and histological studies on the delayed systemic movement of Tobacco Mosaic Virus in Arabidopsis thaliana

BACKGROUND

Viral infections and their spread throughout a plant require numerous interactions between the host and the virus. While new functions of viral proteins involved in these processes have been revealed, current knowledge of host factors involved in the spread of a viral infection is still insufficient. In Arabidopsis thaliana, different ecotypes present varying susceptibilities to Tobacco mosaic virus strain U1 (TMV-U1). The rate of TMV-U1 systemic movement is delayed in ecotype Col-0 when compared with other 13 ecotypes.We followed viral movement through vascular tissue in Col-0 plants by electronic microscopy studies. In addition, the delay in systemic movement of TMV-U1 was genetically studied.

RESULTS

TMV-U1 reaches apical leaves only after 18 days post rosette inoculation (dpi) in Col-0, whereas it is detected at 9 dpi in the Uk-4 ecotype. Genetic crosses between Col-0 and Uk-4 ecotypes, followed by analysis of viral movement in F1 and F2 populations, revealed that this delayed movement correlates with a recessive, monogenic and nuclear locus. The use of selected polymorphic markers showed that this locus, denoted DSTM1 (Delayed Systemic Tobamovirus Movement 1), is positioned on the large arm of chromosome II. Electron microscopy studies following the virion's route in stems of Col-0 infected plants showed the presence of curved structures, instead of the typical rigid rods of TMV-U1. This was not observed in the case of TMV-U1 infection in Uk-4, where the observed virions have the typical rigid rod morphology.

CONCLUSION

The presence of defectively assembled virions observed by electron microscopy in vascular tissue of Col-0 infected plants correlates with a recessive delayed systemic movement trait of TMV-U1 in this ecotype.

Analysis of the grape MYB R2R3 subfamily reveals expanded wine quality-related clades and conserved gene structure organization across Vitis and Arabidopsis genomes

BACKGROUND

The MYB superfamily constitutes the most abundant group of transcription factors described in plants. Members control processes such as epidermal cell differentiation, stomatal aperture, flavonoid synthesis, cold and drought tolerance and pathogen resistance. No genome-wide characterization of this family has been conducted in a woody species such as grapevine. In addition, previous analysis of the recently released grape genome sequence suggested expansion events of several gene families involved in wine quality.

RESULTS

We describe and classify 108 members of the grape R2R3 MYB gene subfamily in terms of their genomic gene structures and similarity to their putative Arabidopsis thaliana orthologues. Seven gene models were derived and analyzed in terms of gene expression and their DNA binding domain structures. Despite low overall sequence homology in the C-terminus of all proteins, even in those with similar functions across Arabidopsis and Vitis, highly conserved motif sequences and exon lengths were found. The grape epidermal cell fate clade is expanded when compared with the Arabidopsis and rice MYB subfamilies. Two anthocyanin MYBA related clusters were identified in chromosomes 2 and 14, one of which includes the previously described grape colour locus. Tannin related loci were also detected with eight candidate homologues in chromosomes 4, 9 and 11.

CONCLUSION

This genome wide transcription factor analysis in Vitis suggests that clade-specific grape R2R3 MYB genes are expanded while other MYB genes could be well conserved compared to Arabidopsis. MYB gene abundance, homology and orientation within particular loci also suggests that expanded MYB clades conferring quality attributes of grapes and wines, such as colour and astringency, could possess redundant, overlapping and cooperative functions.

The N-homologue LRR domain adopts a folding which explains the TMV-Cg-induced HR-like response in sensitive tobacco plants

Following leaf infection with the tobacco mosaic virus (TMV), Nicotiana species that carry the disease resistance N gene develop a hypersensitive response (HR) that blocks the systemic movement of the virus. TMV-sensitive tobacco plants that lack the N gene develop classical disease symptoms following infection with most of the tobamoviruses. However, upon infection with TMV-Cg, these plants display a HR-like response that is unable to limit viral spread. We previously identified the NH gene in sensitive plants; this gene is homologous to the resistance N gene and both belong to the TIR/NBS/LRR family. Isolation and analysis of the NH transcript enabled the prediction of the amino acid sequence in which we detected a leucine-rich repeat domain, proposed to be involved in pathogen recognition. This domain is found in four of five classes of pathogen resistant proteins, in which sequence and structural changes may generate different specificities. In order to study the possible functional role of the LRR domain in the HR-like response, we developed a comparative three-dimensional model for the NH and N gene products, by means of functional and structural domains recognition, secondary structure prediction, domain assignment through profile Hidden Markov Models (HMM) and molecular dynamics (MD) simulations. Based on our results we postulate that the NH protein could adopt a LRR fold with a functional role in the HR-like response. Our two reliable LRR three-dimensional models (N-LRR, NH-LRR) can be used as structural frameworks for future experiments in which the structure-function relationships regarding the protein-protein interaction process may be revealed. Evolutionary aspects of the N and NH genes in Nicotiana species are also discussed.

Isolation of the three grape sub-lineages of B-class MADS-box TM6, PISTILLATA and APETALA3 genes which are differentially expressed during flower and fruit development

The B class of MADS-box floral homeotic genes specifies petal and stamen identity in angiosperms. While this group is one of the most studied in herbaceous plant species, it has remained largely uncharacterized in woody species such as grapevine. Although the B class PI/GLO and AP3/DEF clades have been extensively characterized in model species, the role of the TM6 subgroup within the AP3 clade is not completely understood, since it is absent in Arabidopsis thaliana. In this study, the coding regions of VvTM6 and VvAP3 and the genomic sequence of VvPI, were cloned. VvPI and AtPI were confirmed to be functional homologues by means of complementation of the pi Arabidopsis mutant. Expression analysis revealed that VvPI and VvAP3 transcripts are restricted almost exclusively to inflorescences, although VvPI was detected at low levels in leaves and roots. VvTM6 expresses throughout the plant, with higher levels in flowers and berries. A detailed chronological study of grape flower progression by light microscopy and temporal expression analysis throughout early and late developmental stages, revealed that VvPI expression increases during pollen maturation and decreases between the events of pollination and fertilization, before the cap fall. On the other hand, VvTM6 is expressed in the last stage of anther development. Specific expression of VvAP3 and VvPI was detected in petals and stamens within the flower, while VvTM6 was also expressed in carpels. Moreover, this work provides the first evidence for expression of a TM6-like gene throughout fruit growth and ripening. Even if these genes belong to the same genetic class they could act in different periods and/or tissues during reproductive organ development.

Identification and characterization of a novel tobacco mosaic virus resistance N gene homologue in Nicotiana tabacum plants

Nicotiana tabacum cv. Xanthi nn plants are susceptible to infection by most tobamoviruses (TMV). However, such plants display a partial hypersensitive resistance response (HR-like response) to TMV-Cg. The genetic mechanism of the HR-like response has yet not been determined, but it may involve a gene with a function similar to that of a resistance gene, responsible for HR in resistant plants. We have cloned a gene homologous to the resistance N gene, named NH, from Nicotiana Xanthi nn plants. The coding region of NH is 5.028 base pairs (bp) long and has 82.6% nucleotide identity with the N gene. In contrast to the N gene, the NH gene lacks intron 4 and does not have sites for alternative splicing of intron 3. Analysis of its sequence revealed that NH belongs to the TIR / NSB / LRR gene class. We were able to detect stable levels of NH-transcript in Nicotiana Xanthi nn plants from 0 to 18 h post-inoculation (hpi) with TMV-Cg. Transcript levels increased slightly at 24 hpi and dropped below basal values at 48 hpi. The NH transcript was also detected in a range of resistant Nicotiana plants (N. tabacum Xanthi NN, N. glutinosa, N. glauca and N. rustica) suggesting that NH is a homologue of the N gene, rather than an allele. We have cloned and characterised the NH gene (GenBank acc. no. bankit598573 AY535010) from nn susceptible plants and postulate that this gene might be involved in the HR-like response seen in these plants.