Richard L. Tillett

Transcriptomic and metabolite analyses of Cabernet Sauvignon grape berry development

BackgroundGrape berry development is a dynamic process that involves a complex series of molecular genetic and biochemical changes divided into three major phases. During initial berry growth (Phase I), berry size increases along a sigmoidal growth curve due to cell division and subsequent cell expansion, and organic acids (mainly malate and tartrate), tannins, and hydroxycinnamates accumulate to peak levels. The second major phase (Phase II) is defined as a lag phase in which cell expansion ceases and sugars begin to accumulate. Véraison (the onset of ripening) marks the beginning of the third major phase (Phase III) in which berries undergo a second period of sigmoidal growth due to additional mesocarp cell expansion, accumulation of anthocyanin pigments for berry color, accumulation of volatile compounds for aroma, softening, peak accumulation of sugars (mainly glucose and fructose), and a decline in organic acid accumulation. In order to understand the transcriptional network responsible for controlling berry development, mRNA expression profiling was conducted on berries of V. vinifera Cabernet Sauvignon using the Affymetrix GeneChip®Vitis oligonucleotide microarray ver. 1.0 spanning seven stages of berry development from small pea size berries (E-L stages 31 to 33 as defined by the modified E-L system), through véraison (E-L stages 34 and 35), to mature berries (E-L stages 36 and 38). Selected metabolites were profiled in parallel with mRNA expression profiling to understand the effect of transcriptional regulatory processes on specific metabolite production that ultimately influence the organoleptic properties of wine.ResultsOver the course of berry development whole fruit tissues were found to express an average of 74.5% of probes represented on the Vitis microarray, which has 14,470 Unigenes. Approximately 60% of the expressed transcripts were differentially expressed between at least two out of the seven stages of berry development (28% of transcripts, 4,151 Unigenes, had pronounced (≥2 fold) differences in mRNA expression) illustrating the dynamic nature of the developmental process. The subset of 4,151 Unigenes was split into twenty well-correlated expression profiles. Expression profile patterns included those with declining or increasing mRNA expression over the course of berry development as well as transient peak or trough patterns across various developmental stages as defined by the modified E-L system. These detailed surveys revealed the expression patterns for genes that play key functional roles in phytohormone biosynthesis and response, calcium sequestration, transport and signaling, cell wall metabolism mediating expansion, ripening, and softening, flavonoid metabolism and transport, organic and amino acid metabolism, hexose sugar and triose phosphate metabolism and transport, starch metabolism, photosynthesis, circadian cycles and pathogen resistance. In particular, mRNA expression patterns of transcription factors, abscisic acid (ABA) biosynthesis, and calcium signaling genes identified candidate factors likely to participate in the progression of key developmental events such as véraison and potential candidate genes associated with such processes as auxin partitioning within berry cells, aroma compound production, and pathway regulation and sequestration of flavonoid compounds. Finally, analysis of sugar metabolism gene expression patterns indicated the existence of an alternative pathway for glucose and triose phosphate production that is invoked from véraison to mature berries.ConclusionThese results reveal the first high-resolution picture of the transcriptome dynamics that occur during seven stages of grape berry development. This work also establishes an extensive catalog of gene expression patterns for future investigations aimed at the dissection of the transcriptional regulatory hierarchies that govern berry development in a widely grown cultivar of wine grape. More importantly, this analysis identified a set of previously unknown genes potentially involved in critical steps associated with fruit development that can now be subjected to functional testing.

Tissue-specific mRNA expression profiling in grape berry tissues

BackgroundBerries of grape (Vitis vinifera) contain three major tissue types (skin, pulp and seed) all of which contribute to the aroma, color, and flavor characters of wine. The pericarp, which is composed of the exocarp (skin) and mesocarp (pulp), not only functions to protect and feed the developing seed, but also to assist in the dispersal of the mature seed by avian and mammalian vectors. The skin provides volatile and nonvolatile aroma and color compounds, the pulp contributes organic acids and sugars, and the seeds provide condensed tannins, all of which are important to the formation of organoleptic characteristics of wine. In order to understand the transcriptional network responsible for controlling tissue-specific mRNA expression patterns, mRNA expression profiling was conducted on each tissue of mature berries of V. vinifera Cabernet Sauvignon using the Affymetrix GeneChip®Vitis oligonucleotide microarray ver. 1.0. In order to monitor the influence of water-deficit stress on tissue-specific expression patterns, mRNA expression profiles were also compared from mature berries harvested from vines subjected to well-watered or water-deficit conditions.ResultsOverall, berry tissues were found to express approximately 76% of genes represented on the Vitis microarray. Approximately 60% of these genes exhibited significant differential expression in one or more of the three major tissue types with more than 28% of genes showing pronounced (2-fold or greater) differences in mRNA expression. The largest difference in tissue-specific expression was observed between the seed and pulp/skin. Exocarp tissue, which is involved in pathogen defense and pigment production, showed higher mRNA abundance relative to other berry tissues for genes involved with flavonoid biosynthesis, pathogen resistance, and cell wall modification. Mesocarp tissue, which is considered a nutritive tissue, exhibited a higher mRNA abundance of genes involved in cell wall function and transport processes. Seeds, which supply essential resources for embryo development, showed higher mRNA abundance of genes encoding phenylpropanoid biosynthetic enzymes, seed storage proteins, and late embryogenesis abundant proteins. Water-deficit stress affected the mRNA abundance of 13% of the genes with differential expression patterns occurring mainly in the pulp and skin. In pulp and seed tissues transcript abundance in most functional categories declined in water-deficit stressed vines relative to well-watered vines with transcripts for storage proteins and novel (no-hit) functional assignments being over represented. In the skin of berries from water-deficit stressed vines, however, transcripts from several functional categories including general phenypropanoid and ethylene metabolism, pathogenesis-related responses, energy, and interaction with the environment were significantly over-represented.ConclusionThese results revealed novel insights into the tissue-specific expression mRNA expression patterns of an extensive repertoire of genes expressed in berry tissues. This work also establishes an extensive catalogue of gene expression patterns for future investigations aimed at the dissection of the transcriptional regulatory hierarchies that govern tissue-specific expression patterns associated with tissue differentiation within berries. These results also confirmed that water-deficit stress has a profound effect on mRNA expression patterns particularly associated with the biosynthesis of aroma and color metabolites within skin and pulp tissues that ultimately impact wine quality.

Transcriptomic analysis of the late stages of grapevine (Vitis vinifera cv. Cabernet Sauvignon) berry ripening reveals significant induction of ethylene signaling and flavor pathways in the skin

BackgroundGrapevine berry, a nonclimacteric fruit, has three developmental stages; the last one is when berry color and sugar increase. Flavors derived from terpenoid and fatty acid metabolism develop at the very end of this ripening stage. The transcriptomic response of pulp and skin of Cabernet Sauvignon berries in the late stages of ripening between 22 and 37 °Brix was assessed using whole-genome micorarrays.ResultsThe transcript abundance of approximately 18,000 genes changed with °Brix and tissue type. There were a large number of changes in many gene ontology (GO) categories involving metabolism, signaling and abiotic stress. GO categories reflecting tissue differences were overrepresented in photosynthesis, isoprenoid metabolism and pigment biosynthesis. Detailed analysis of the interaction of the skin and pulp with °Brix revealed that there were statistically significantly higher abundances of transcripts changing with °Brix in the skin that were involved in ethylene signaling, isoprenoid and fatty acid metabolism. Many transcripts were peaking around known optimal fruit stages for flavor production. The transcript abundance of approximately two-thirds of the AP2/ERF superfamily of transcription factors changed during these developmental stages. The transcript abundance of a unique clade of ERF6-type transcription factors had the largest changes in the skin and clustered with genes involved in ethylene, senescence, and fruit flavor production including ACC oxidase, terpene synthases, and lipoxygenases. The transcript abundance of important transcription factors involved in fruit ripening was also higher in the skin.ConclusionsA detailed analysis of the transcriptome dynamics during late stages of ripening of grapevine berries revealed that these berries went through massive transcriptional changes in gene ontology categories involving chemical signaling and metabolism in both the pulp and skin, particularly in the skin. Changes in the transcript abundance of genes involved in the ethylene signaling pathway of this nonclimacteric fruit were statistically significant in the late stages of ripening when the production of transcripts for important flavor and aroma compounds were at their highest. Ethylene transcription factors known to play a role in leaf senescence also appear to play a role in fruit senescence. Ethylene may play a bigger role than previously thought in this non-climacteric fruit.

Identification of tissue-specific, abiotic stress-responsive gene expression patterns in wine grape (Vitis vinifera L.) based on curation and mining of large-scale EST data sets

BackgroundAbiotic stresses, such as water deficit and soil salinity, result in changes in physiology, nutrient use, and vegetative growth in vines, and ultimately, yield and flavor in berries of wine grape, Vitis vinifera L. Large-scale expressed sequence tags (ESTs) were generated, curated, and analyzed to identify major genetic determinants responsible for stress-adaptive responses. Although roots serve as the first site of perception and/or injury for many types of abiotic stress, EST sequencing in root tissues of wine grape exposed to abiotic stresses has been extremely limited to date. To overcome this limitation, large-scale EST sequencing was conducted from root tissues exposed to multiple abiotic stresses.ResultsA total of 62,236 expressed sequence tags (ESTs) were generated from leaf, berry, and root tissues from vines subjected to abiotic stresses and compared with 32,286 ESTs sequenced from 20 public cDNA libraries. Curation to correct annotation errors, clustering and assembly of the berry and leaf ESTs with currently available V. vinifera full-length transcripts and ESTs yielded a total of 13,278 unique sequences, with 2302 singletons and 10,976 mapped to V. vinifera gene models. Of these, 739 transcripts were found to have significant differential expression in stressed leaves and berries including 250 genes not described previously as being abiotic stress responsive. In a second analysis of 16,452 ESTs from a normalized root cDNA library derived from roots exposed to multiple, short-term, abiotic stresses, 135 genes with root-enriched expression patterns were identified on the basis of their relative EST abundance in roots relative to other tissues.ConclusionsThe large-scale analysis of relative EST frequency counts among a diverse collection of 23 different cDNA libraries from leaf, berry, and root tissues of wine grape exposed to a variety of abiotic stress conditions revealed distinct, tissue-specific expression patterns, previously unrecognized stress-induced genes, and many novel genes with root-enriched mRNA expression for improving our understanding of root biology and manipulation of rootstock traits in wine grape. mRNA abundance estimates based on EST library-enriched expression patterns showed only modest correlations between microarray and quantitative, real-time reverse transcription-polymerase chain reaction (qRT-PCR) methods highlighting the need for deep-sequencing expression profiling methods.

Expressed sequence tag (EST) profiling in hyper saline shocked Dunaliella salina reveals high expression of protein synthetic apparatus components

The unicellular, halotolerant, green alga, Dunaliella salina (Chlorophyceae) has the unique ability to adapt and grow in a wide range of salt conditions from about 0.05 to 5.5M. To better understand the molecular basis of its salinity tolerance, a complementary DNA (cDNA) library was constructed from D. salina cells adapted to 2.5M NaCl, salt-shocked at 3.4M NaCl for 5h, and used to generate an expressed sequence tag (EST) database. ESTs were obtained for 2831 clones representing 1401 unique transcripts. Putative functions were assigned to 1901 (67.2%) ESTs after comparison with protein databases. An additional 154 (5.4%) ESTs had significant similarity to known sequences whose functions are unclear and 776 (27.4%) had no similarity to known sequences. For those D. salina ESTs for which functional assignments could be made, the largest functional categories included protein synthesis (35.7%), energy (photosynthesis) (21.4%), primary metabolism (13.8%) and protein fate (6.8%). Within the protein synthesis category, the vast majority of ESTs (80.3%) encoded ribosomal proteins representing about 95% of the approximately 82 subunits of the cytosolic ribosome indicating that D. salina invests substantial resources in the production and maintenance of protein synthesis. The increased mRNA expression upon salinity shock was verified for a small set of selected genes by real-time, quantitative reverse-transcription-polymerase chain reaction (qRT-PCR). This EST collection also provided important new insights into the genetic underpinnings for the biosynthesis and utilization of glycerol and other osmoprotectants, the carotenoid biosynthetic pathway, reactive oxygen-scavenging enzymes, and molecular chaperones (heat shock proteins) not described previously for D. salina. EST discovery also revealed the existence of RNA interference and signaling pathways associated with osmotic stress adaptation. The unknown ESTs described here provide a rich resource for the identification of novel genes associated with the mechanistic basis of salinity stress tolerance and other stress-adaptive traits.

Sporobolus stapfianus: Insights into desiccation tolerance in the resurrection grasses from linking transcriptomics to metabolomics

BackgroundUnderstanding the response of resurrection angiosperms to dehydration and rehydration is critical for deciphering the mechanisms of how plants cope with the rigors of water loss from their vegetative tissues. We have focused our studies on the C4 resurrection grass, Sporobolus stapfianus Gandoger, as a member of a group of important forage grasses.MethodsWe have combined non-targeted metabolomics with transcriptomics, via a NimbleGen array platform, to develop an understanding of how gene expression and metabolite profiles can be linked to generate a more detailed mechanistic appreciation of the cellular response to both desiccation and rehydration.ResultsThe rehydration transcriptome and metabolome are primarily geared towards the rapid return of photosynthesis, energy metabolism, protein turnover, and protein synthesis during the rehydration phase. However, there are some metabolites associated with ROS protection that remain elevated during rehydration, most notably the tocopherols. The analysis of the dehydration transcriptome reveals a strong concordance between transcript abundance and the associated metabolite abundance reported earlier, but only in responses that are directly related to cellular protection during dehydration: carbohydrate metabolism and redox homeostasis. The transcriptome response also provides strong support for the involvement of cellular protection processes as exemplified by the increases in the abundance of transcripts encoding late embryogenesis abundant (LEA) proteins, anti-oxidant enzymes, early light-induced proteins (ELIP) proteins, and cell-wall modification enzymes. There is little concordance between transcript and metabolite abundance for processes such as amino acid metabolism that do not appear to contribute directly to cellular protection, but are nonetheless important for the desiccation tolerant phenotype of S. stapfianus.ConclusionsThe transcriptomes of both dehydration and rehydration offer insight into the complexity of the regulation of responses to these processes that involve complex signaling pathways and associated transcription factors. ABA appears to be important in the control of gene expression in both the latter stages of the dehydration and the early stages of rehydration. These findings add to the growing body of information detailing how plants tolerate and survive the severe cellular perturbations of dehydration, desiccation, and rehydration.

Comparative transcriptomics of mountain pine beetle pheromone-biosynthetic tissues and functional analysis of CYP6DE3

BackgroundThe mountain pine beetle (MPB, Dendroctonus ponderosae Hopkins) is a highly destructive pest of pine forests in western North America. During flight to a new host tree and initiation of feeding, mountain pine beetles release aggregation pheromones. The biosynthetic pathways of these pheromones are sex-specific and localized in the midgut and fat body, but the enzymes involved have not all been identified or characterized.ResultsWe used a comparative RNA-Seq analysis between fed and unfed male and female MPB midguts and fat bodies to identify candidate genes involved in pheromone biosynthesis. The 13,407 potentially unique transcripts showed clear separation based on feeding state and gender. Gene co-expression network construction and examination using petal identified gene groups that were tightly connected. This, as well as other co-expression and gene ontology analyses, identified all four known pheromone biosynthetic genes, confirmed the tentative identification of four others from a previous study, and suggested nine novel candidates. One cytochrome P450 monooxygenase, CYP6DE3, identified as a possible exo-brevicomin-biosynthetic enzyme in this study, was functionally characterized and likely is involved in resin detoxification rather than pheromone biosynthesis.ConclusionsOur analysis supported previously characterized pheromone-biosynthetic genes involved in exo-brevicomin and frontalin biosynthesis and identified a number of candidate cytochrome P450 monooxygenases and a putative cyclase for further studies. Functional analyses of CYP6DE3 suggest its role in resin detoxification and underscore the limitation of using high-throughput data to tentatively identify candidate genes. Further functional analyses of candidate genes found in this study should lead to the full characterization of MPB pheromone biosynthetic pathways and the identification of molecular targets for possible pest management strategies.

Previously Unidentified Single Nucleotide Polymorphisms in HIV/AIDS Cases Associate with Clinical Parameters and Disease Progression

The genetic background of an individual plays an important role in the progression of HIV infection to AIDS. Identifying previously unknown or uncharacterized single nucleotide polymorphisms (SNPs) that associate with disease progression may reveal important therapeutic targets and provide a greater understanding of disease pathogenesis. In the present study, we employed ultra-high multiplex PCR on an Ion Torrent next-generation sequencing platform to sequence 23 innate immune genes from 94 individuals with HIV/AIDS. This data was used to identify potential associations of SNPs with clinical parameters and disease progression. SNPs that associated with an increased viral load were identified in the genes for the interleukin 15 receptor (IL15RA), toll-like receptor 7 (TLR7), tripartite motif-containing protein 5 (TRIM5), and two killer-cell immunoglobulin-like receptors (KIR2DL1 and KIR2DL3). Additionally, SNPs that associated with progression from HIV infection to AIDS were identified in two 2′-5′-oligoadenylate synthetase genes (OAS2 and OAS3). In contrast, other SNPs identified in OAS2 and OAS3 genes, as well as in the TRIM5 and KIR2DS4 genes, were associated with a slower progression of disease. Taken together, our data demonstrates the utility of ultra-high multiplex PCR in identifying polymorphisms of potential clinical significance and further,identifies SNPs that may play a role in HIV pathogenesis.

The common transcriptional subnetworks of the grape berry skin in the late stages of ripening

BackgroundWine grapes are important economically in many countries around the world. Defining the optimum time for grape harvest is a major challenge to the grower and winemaker. Berry skins are an important source of flavor, color and other quality traits in the ripening stage. Senescent-like processes such as chloroplast disorganization and cell death characterize the late ripening stage.ResultsTo better understand the molecular and physiological processes involved in the late stages of berry ripening, RNA-seq analysis of the skins of seven wine grape cultivars (Cabernet Franc, Cabernet Sauvignon, Merlot, Pinot Noir, Chardonnay, Sauvignon Blanc and Semillon) was performed. RNA-seq analysis identified approximately 2000 common differentially expressed genes for all seven cultivars across four different berry sugar levels (20 to 26 °Brix). Network analyses, both a posteriori (standard) and a priori (gene co-expression network analysis), were used to elucidate transcriptional subnetworks and hub genes associated with traits in the berry skins of the late stages of berry ripening. These independent approaches revealed genes involved in photosynthesis, catabolism, and nucleotide metabolism. The transcript abundance of most photosynthetic genes declined with increasing sugar levels in the berries. The transcript abundance of other processes increased such as nucleic acid metabolism, chromosome organization and lipid catabolism. Weighted gene co-expression network analysis (WGCNA) identified 64 gene modules that were organized into 12 subnetworks of three modules or more and six higher order gene subnetworks. Some gene subnetworks were highly correlated with sugar levels and some subnetworks were highly enriched in the chloroplast and nucleus. The petal R package was utilized independently to construct a true small-world and scale-free complex gene co-expression network model. A subnetwork of 216 genes with the highest connectivity was elucidated, consistent with the module results from WGCNA. Hub genes in these subnetworks were identified including numerous members of the core circadian clock, RNA splicing, proteolysis and chromosome organization. An integrated model was constructed linking light sensing with alternative splicing, chromosome remodeling and the circadian clock.ConclusionsA common set of differentially expressed genes and gene subnetworks from seven different cultivars were examined in the skin of the late stages of grapevine berry ripening. A densely connected gene subnetwork was elucidated involving a complex interaction of berry senescent processes (autophagy), catabolism, the circadian clock, RNA splicing, proteolysis and epigenetic regulation. Hypotheses were induced from these data sets involving sugar accumulation, light, autophagy, epigenetic regulation, and fruit development. This work provides a better understanding of berry development and the transcriptional processes involved in the late stages of ripening.

A comparison of heat-stress transcriptome changes between wild-type Arabidopsis pollen and a heat-sensitive mutant harboring a knockout of cyclic nucleotide-gated cation channel 16 (cngc16)

BackgroundIn flowering plants, the male gametophyte (pollen) is one of the most vulnerable cells to temperature stress. In Arabidopsis thaliana, a pollen-specific CyclicNucleotide-Gated cationChannel 16 (cngc16), is required for plant reproduction under temperature-stress conditions. Plants harboring a cncg16 knockout are nearly sterile under conditions of hot days and cold nights. To understand the underlying cause, RNA-Seq was used to compare the pollen transcriptomes of wild type (WT) and cngc16 under normal and heat stress (HS) conditions.ResultsHere we show that a heat-stress response (HSR) in WT pollen resulted in 2102 statistically significant transcriptome changes (≥ 2-fold changes with adjusted p-value ≤0.01), representing approximately 15% of 14,226 quantified transcripts. Of these changes, 89 corresponded to transcription factors, with 27 showing a preferential expression in pollen over seedling tissues. In contrast to WT, cngc16 pollen showed 1.9-fold more HS-dependent changes (3936 total, with 2776 differences between WT and cngc16). In a quantitative direct comparison between WT and cngc16 transcriptomes, the number of statistically significant differences increased from 21 pre-existing differences under normal conditions to 192 differences under HS. Of the 20 HS-dependent changes in WT that were most different in cngc16, half corresponded to genes encoding proteins predicted to impact cell wall features or membrane dynamics.ConclusionsResults here define an extensive HS-dependent reprogramming of approximately 15% of the WT pollen transcriptome, and identify at least 27 transcription factor changes that could provide unique contributions to a pollen HSR. The number of statistically significant transcriptome differences between WT and cngc16 increased by more than 9-fold under HS, with most of the largest magnitude changes having the potential to specifically impact cell walls or membrane dynamics, and thereby potentiate cngc16 pollen to be hypersensitive to HS. However, HS-hypersensitivity could also be caused by the extensive number of differences throughout the transcriptome having a cumulative effect on multiple cellular pathways required for tip growth and fertilization. Regardless, results here support a model in which a functional HS-dependent reprogramming of the pollen transcriptome requires a specific calcium-permeable Cyclic Nucleotide-Gated cation Channel, CNGC16.