Sara Zenoni

Characterization of Transcriptional Complexity during Berry Development in Vitis vinifera Using RNA-Seq

The development of massively parallel sequencing technologies enables the sequencing of total cDNA (RNA-Seq) to derive accurate measure of individual gene expression, differential splicing activity, and to discover novel regions of transcription, dramatically changing the way that the functional complexity of transcriptomes can be studied. Here we report on the first use of RNA-Seq to gain insight into the wide range of transcriptional responses that are associated with berry development in Vitis vinifera ‘Corvina’. More than 59 million sequence reads, 36 to 44 bp in length, were generated from three developmental stages: post setting, véraison, and ripening. The sequence reads were aligned onto the 8.4-fold draft sequence of the Pinot Noir 40024 genome and then analyzed to measure gene expression levels, to detect alternative splicing events, and expressed single nucleotide polymorphisms. We detected 17,324 genes expressed during berry development, 6,695 of which were expressed in a stage-specific manner, suggesting differences in expression for genes in numerous functional categories and a significant transcriptional complexity. This exhaustive overview of gene expression dynamics demonstrates the utility of RNA-Seq for identifying single nucleotide polymorphisms and splice variants and for describing how plant transcriptomes change during development.

The grapevine expression atlas reveals a deep transcriptome shift driving the entire plant into a maturation program.

The authors developed a comprehensive transcriptome atlas in grapevine by comparing the genes expressed in 54 diverse samples accounting for ∼91% of all known grapevine genes. Using a panel of different statistical techniques, they found that the whole plant undergoes transcriptomic reprogramming, driving it towards maturity. We developed a genome-wide transcriptomic atlas of grapevine (Vitis vinifera) based on 54 samples representing green and woody tissues and organs at different developmental stages as well as specialized tissues such as pollen and senescent leaves. Together, these samples expressed ∼91% of the predicted grapevine genes. Pollen and senescent leaves had unique transcriptomes reflecting their specialized functions and physiological status. However, microarray and RNA-seq analysis grouped all the other samples into two major classes based on maturity rather than organ identity, namely, the vegetative/green and mature/woody categories. This division represents a fundamental transcriptomic reprogramming during the maturation process and was highlighted by three statistical approaches identifying the transcriptional relationships among samples (correlation analysis), putative biomarkers (O2PLS-DA approach), and sets of strongly and consistently expressed genes that define groups (topics) of similar samples (biclustering analysis). Gene coexpression analysis indicated that the mature/woody developmental program results from the reiterative coactivation of pathways that are largely inactive in vegetative/green tissues, often involving the coregulation of clusters of neighboring genes and global regulation based on codon preference. This global transcriptomic reprogramming during maturation has not been observed in herbaceous annual species and may be a defining characteristic of perennial woody plants.

Genomic and transcriptomic analysis of the AP2/ERF superfamily in Vitis vinifera

BackgroundThe AP2/ERF protein family contains transcription factors that play a crucial role in plant growth and development and in response to biotic and abiotic stress conditions in plants. Grapevine (Vitis vinifera) is the only woody crop whose genome has been fully sequenced. So far, no detailed expression profile of AP2/ERF-like genes is available for grapevine.ResultsAn exhaustive search for AP2/ERF genes was carried out on the Vitis vinifera genome and their expression profile was analyzed by Real-Time quantitative PCR (qRT-PCR) in different vegetative and reproductive tissues and under two different ripening stages.One hundred and forty nine sequences, containing at least one ERF domain, were identified. Specific clusters within the AP2 and ERF families showed conserved expression patterns reminiscent of other species and grapevine specific trends related to berry ripening. Moreover, putative targets of group IX ERFs were identified by co-expression and protein similarity comparisons.ConclusionsThe grapevine genome contains an amount of AP2/ERF genes comparable to that of other dicot species analyzed so far. We observed an increase in the size of specific groups within the ERF family, probably due to recent duplication events. Expression analyses in different aerial tissues display common features previously described in other plant systems and introduce possible new roles for members of some ERF groups during fruit ripening. The presented analysis of AP2/ERF genes in grapevine provides the bases for studying the molecular regulation of berry development and the ripening process.

Identification of Putative Stage-Specific Grapevine Berry Biomarkers and Omics Data Integration into Networks

The analysis of grapevine (Vitis vinifera) berries at the transcriptomic, proteomic, and metabolomic levels can provide great insight into the molecular events underlying berry development and postharvest drying (withering). However, the large and very different data sets produced by such investigations are difficult to integrate. Here, we report the identification of putative stage-specific biomarkers for berry development and withering and, to our knowledge, the first integrated systems-level study of these processes. Transcriptomic, proteomic, and metabolomic data were integrated using two different strategies, one hypothesis free and the other hypothesis driven. A multistep hypothesis-free approach was applied to data from four developmental stages and three withering intervals, with integration achieved using a hierarchical clustering strategy based on the multivariate bidirectional orthogonal projections to latent structures technique. This identified stage-specific functional networks of linked transcripts, proteins, and metabolites, providing important insights into the key molecular processes that determine the quality characteristics of wine. The hypothesis-driven approach was used to integrate data from three withering intervals, starting with subdata sets of transcripts, proteins, and metabolites. We identified transcripts and proteins that were modulated during withering as well as specific classes of metabolites that accumulated at the same time and used these to select subdata sets of variables. The multivariate bidirectional orthogonal projections to latent structures technique was then used to integrate the subdata sets, identifying variables representing selected molecular processes that take place specifically during berry withering. The impact of this holistic approach on our knowledge of grapevine berry development and withering is discussed.

Genome-wide analysis of the grapevine stilbene synthase multigenic family: genomic organization and expression profiles upon biotic and abiotic stresses

BackgroundPlant stilbenes are a small group of phenylpropanoids, which have been detected in at least 72 unrelated plant species and accumulate in response to biotic and abiotic stresses such as infection, wounding, UV-C exposure and treatment with chemicals. Stilbenes are formed via the phenylalanine/polymalonate-route, the last step of which is catalyzed by the enzyme stilbene synthase (STS), a type III polyketide synthase (PKS). Stilbene synthases are closely related to chalcone synthases (CHS), the key enzymes of the flavonoid pathway, as illustrated by the fact that both enzymes share the same substrates. To date, STSs have been cloned from peanut, pine, sorghum and grapevine, the only stilbene-producing fruiting-plant for which the entire genome has been sequenced. Apart from sorghum, STS genes appear to exist as a family of closely related genes in these other plant species.ResultsIn this study a complete characterization of the STS multigenic family in grapevine has been performed, commencing with the identification, annotation and phylogenetic analysis of all members and integration of this information with a comprehensive set of gene expression analyses including healthy tissues at differential developmental stages and in leaves exposed to both biotic (downy mildew infection) and abiotic (wounding and UV-C exposure) stresses. At least thirty-three full length sequences encoding VvSTS genes were identified, which, based on predicted amino acid sequences, cluster in 3 principal groups designated A, B and C. The majority of VvSTS genes cluster in groups B and C and are located on chr16 whereas the few gene family members in group A are found on chr10. Microarray and mRNA-seq expression analyses revealed different patterns of transcript accumulation between the different groups of VvSTS family members and between VvSTSs and VvCHSs. Indeed, under certain conditions the transcriptional response of VvSTS and VvCHS genes appears to be diametrically opposed suggesting that flow of carbon between these two competing metabolic pathways is tightly regulated at the transcriptional level.ConclusionsThis study represents an overview of the expression pattern of each member of the STS gene family in grapevine under both constitutive and stress-induced conditions. The results strongly indicate the existence of a transcriptional subfunctionalization amongst VvSTSs and provide the foundation for further functional investigations about the role and evolution of this large gene family. Moreover, it represents the first study to clearly show the differential regulation of VvCHS and VvSTS genes, suggesting the involvement of transcription factors (TFs) in both the activation and repression of these genes.

Downregulation of the Petunia hybrida α-Expansin Gene PhEXP1 Reduces the Amount of Crystalline Cellulose in Cell Walls and Leads to Phenotypic Changes in Petal Limbs

The expansins comprise a family of proteins that appear to be involved in the disruption of the noncovalent bonds between cellulose microfibrils and cross-linking glycans, thereby promoting wall creep. To understand better the expansion process in Petunia hybrida (petunia) flowers, we isolated a cDNA corresponding to the PhEXP1 α-expansin gene of P. hybrida. Evaluation of the tissue specificity and temporal expression pattern demonstrated that PhEXP1 is preferentially expressed in petal limbs during development. To determine the function of PhEXP1, we used a transgenic antisense approach, which was found to cause a decrease in petal limb size, a reduction in the epidermal cell area, and alterations in cell wall morphology and composition. The diminished cell wall thickness accompanied by a reduction in crystalline cellulose indicates that the activity of PhEXP1 is associated with cellulose metabolism. Our results suggest that expansins play a role in the assembly of the cell wall by affecting either cellulose synthesis or deposition.

The plasticity of the grapevine berry transcriptome

BackgroundPhenotypic plasticity refers to the range of phenotypes a single genotype can express as a function of its environment. These phenotypic variations are attributable to the effect of the environment on the expression and function of genes influencing plastic traits. We investigated phenotypic plasticity in grapevine by comparing the berry transcriptome in a single clone of the vegetatively-propagated common grapevine species Vitis vinifera cultivar Corvina through 3 consecutive growth years cultivated in 11 different vineyards in the Verona area of Italy.ResultsMost of the berry transcriptome clustered by year of growth rather than common environmental conditions or viticulture practices, and transcripts related to secondary metabolism showed high sensitivity towards different climates, as confirmed also by metabolomic data obtained from the same samples. When analyzed in 11 vineyards during 1 growth year, the environmentally-sensitive berry transcriptome comprised 5% of protein-coding genes and 18% of the transcripts modulated during berry development. Plastic genes were particularly enriched in ontology categories such as transcription factors, translation, transport, and secondary metabolism. Specific plastic transcripts were associated with groups of vineyards sharing common viticulture practices or environmental conditions, and plastic transcriptome reprogramming was more intense in the year characterized by extreme weather conditions. We also identified a set of genes that lacked plasticity, showing either constitutive expression or similar modulation in all berries.ConclusionsOur data reveal candidate genes potentially responsible for the phenotypic plasticity of grapevine and provide the first step towards the characterization of grapevine transcriptome plasticity under different agricultural systems.

The Phenylpropanoid Pathway Is Controlled at Different Branches by a Set of R2R3-MYB C2 Repressors in Grapevine

A set of transcriptional repressors negatively regulates the expression of genes involved in different branches of the phenylpropanoid pathway. 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.

De novo transcriptome characterization of Vitis vinifera cv. Corvina unveils varietal diversity.

BackgroundPlants such as grapevine (Vitis spp.) display significant inter-cultivar genetic and phenotypic variation. The genetic components underlying phenotypic diversity in grapevine must be understood in order to disentangle genetic and environmental factors.ResultsWe have shown that cDNA sequencing by RNA-seq is a robust approach for the characterization of varietal diversity between a local grapevine cultivar (Corvina) and the PN40024 reference genome. We detected 15,161 known genes including 9463 with novel splice isoforms, and identified 2321 potentially novel protein-coding genes in non-annotated or unassembled regions of the reference genome. We also discovered 180 apparent private genes in the Corvina genome which were missing from the reference genome.ConclusionsThe de novo assembly approach allowed a substantial amount of the Corvina transcriptome to be reconstructed, improving known gene annotations by robustly defining gene structures, annotating splice isoforms and detecting genes without annotations. The private genes we discovered are likely to be nonessential but could influence certain cultivar-specific characteristics. Therefore, the application of de novo transcriptome assembly should not be restricted to species lacking a reference genome because it can also improve existing reference genome annotations and identify novel, cultivar-specific genes.

Genome and transcriptome analysis of the grapevine (Vitis vinifera L.) WRKY gene family

The plant WRKY gene family represents an ancient and complex class of zinc-finger transcription factors (TFs) that are involved in the regulation of various physiological processes, such as development and senescence, and in plant response to many biotic and abiotic stresses. Despite the growing number of studies on the genomic organisation of WRKY gene family in different species, little information is available about this family in grapevine (Vitis vinifera L.). In the present study, a total number of 59 putative grapevine WRKY transcription factors (VvWRKYs) were identified based on the analysis of various genomic and proteomic grapevine databases. According to their structural and phylogentic features, the identified grapevine WRKY transcription factors were classified into three main groups. In order to shed light into their regulatory roles in growth and development as well as in response to biotic and abiotic stress in grapevine, the VvWRKYs expression profiles were examined in publicly available microarray data. Bioinformatics analysis of these data revealed distinct temporal and spatial expression patterns of VvWRKYs in various tissues, organs and developmental stages, as well as in response to biotic and abiotic stresses. To also extend our analysis to situations not covered by the arrays and to validate our results, the expression profiles of selected VvWRKYs in response to drought stress, Erysiphe necator (powdery mildew) infection, and hormone treatments (salicilic acid and ethylene), were investigated by quantitative real-time reverse transcription PCR (qRT-PCR). The present study provides a foundation for further comparative genomics and functional studies of this important class of transcriptional regulators in grapevine.