Multi-omics atlas of combinatorial abiotic stress responses in wheat.

Field-grown crops rarely experience growth conditions in which yield can be maximized. Environmental stresses occur in combination, with advancements in crop tolerance further complicated by its polygenic nature. Strategic targeting of causal genes is required to meet future crop production needs. Here, we employed a systems biology approach in wheat (Triticum aestivum L.) to investigate physio-metabolic adjustments and transcriptome reprogramming involved in acclimations to heat, drought, salinity and all combinations therein. A significant shift in magnitude and complexity of plant response was evident across stress scenarios based on the agronomic losses, increased proline concentrations and 8.7-fold increase in unique differentially expressed transcripts (DETs) observed under the triple stress condition. Transcriptome data from all stress treatments were assembled into an online, open access eFP browser for visualizing gene expression during abiotic stress. Weighted gene co-expression network analysis revealed 152 hub genes of which 32% contained the ethylene-responsive element binding factor-associated amphiphilic repression (EAR) transcriptional repression motif. Cross-referencing against the 31 DETs common to all stress treatments isolated TaWRKY33 as a leading candidate for greater plant tolerance to combinatorial stresses. Integration of our findings with available literature on gene functional characterization allowed us to further suggest flexible gene combinations for future adaptive gene stacking in wheat. Our approach demonstrates the strength of robust multi-omics-based data resources for gene discovery in complex environmental conditions. Accessibility of such datasets will promote cross-validation of candidate genes across studies and aid in accelerating causal gene validation for crop resiliency.

Spatiotemporal transcriptomics and metabolic profiling provide insights into gene regulatory networks during lentil seed development.

Lentil (Lens culinaris Medik.) is a nutritious legume with seeds rich in protein, minerals and an array of diverse specialized metabolites. The formation of a seed requires regulation and tight coordination of developmental programs to form the embryo, endosperm and seed coat compartments, which determines the structure and composition of mature seed and thus its end-use quality. Understanding the molecular and cellular events and metabolic processes of seed development is essential for improving lentil yield and seed nutritional value. However, such information remains largely unknown, especially at the seed compartment level. In this study, we generated high-resolution spatiotemporal gene expression profiles in lentil embryo, seed coat and whole seeds from fertilization through maturation. Apart from anatomic differences between the embryo and seed coat, comparative transcriptomics and weighted gene co-expression network analysis revealed embryo- and seed coat-specific genes and gene modules predominant in specific tissues and stages, which highlights distinct genetic programming. Furthermore, we investigated the dynamic profiles of flavonoid, isoflavone, phytic acid and saponin in seed compartments across seed development. Coupled with transcriptome data, we identified sets of candidate genes involved in the biosynthesis of these metabolites. The global view of the transcriptional and metabolic changes of lentil seed tissues throughout development provides a valuable resource for dissecting the genetic control of secondary metabolism and development of molecular tools for improving seed nutritional quality.

Genomic selection for agronomic traits in a winter wheat breeding program.

KEY MESSAGE: rAMP-seq based genomic selection for agronomic traits has been shown to be a useful tool for winter wheat breeding programs by increasing the rate of genetic gain. Genomic selection (GS) is an effective strategy to employ in a breeding program that focuses on optimizing quantitative traits, which results in the ability for breeders to select the best genotypes. GS was incorporated into a breeding program to determine the potential for implementation on an annual basis, with emphasis on selecting optimal parents and decreasing the time and costs associated with phenotyping large numbers of genotypes. The design options for applying repeat amplification sequencing (rAMP-seq) in bread wheat were explored, and a low-cost single primer pair strategy was implemented. A total of 1870 winter wheat genotypes were phenotyped and genotyped using rAMP-seq. The optimization of training to testing population size showed that the 70:30 ratio provided the most consistent prediction accuracy. Three GS models were tested, rrBLUP, RKHS and feed-forward neural networks using the University of Guelph Winter Wheat Breeding Program (UGWWBP) and Elite-UGWWBP populations. The models performed equally well for both populations and did not differ in prediction accuracy (r) for most agronomic traits, with the exception of yield, where RKHS performed the best with an r = 0.34 and 0.39 for each population, respectively. The ability to operate a breeding program where multiple selection strategies, including GS, are utilized will lead to higher efficiency in the program and ultimately lead to a higher rate of genetic gain.

Evolutionary divergence in embryo and seed coat development of U's Triangle Brassica species illustrated by a spatiotemporal transcriptome atlas.

The economically valuable Brassica species include the six related members of U's Triangle. Despite the agronomic and economic importance of these Brassicas, the impacts of evolution and relatively recent domestication events on the genetic landscape of seed development have not been comprehensively examined in these species. Here we present a 3D transcriptome atlas for the six species of U's Triangle, producing a unique resource that captures gene expression data for the major subcompartments of the seed, from the unfertilized ovule to the mature embryo and seed coat. This comprehensive dataset for seed development in tetraploid and ancestral diploid Brassicas provides new insights into evolutionary divergence and expression bias at the gene and subgenome levels during the domestication of these valued crop species. Comparisons of gene expression associated with regulatory networks and metabolic pathways operating in the embryo and seed coat during seed development reveal differences in storage reserve accumulation and fatty acid metabolism among the six Brassica species. This study illustrates the genetic underpinnings of seed traits and the selective pressures placed on seed production, providing an immense resource for continued investigation of Brassica polyploid biology, genomics and evolution.

Wheat transcriptome profiling reveals abscisic and gibberellic acid treatments regulate early-stage phytohormone defense signaling, cell wall fortification, and metabolic switches following Fusarium graminearum-challenge.

BACKGROUND: Treatment of wheat with the phytohormones abscisic acid (ABA) and gibberellic acid (GA) has been shown to affect Fusarium head blight (FHB) disease severity. However, the molecular mechanisms underlying the elicited phenotypes remain unclear. Toward addressing this gap in our knowledge, global transcriptomic profiling was applied to the FHB-susceptible wheat cultivar 'Fielder' to map the regulatory responses effected upon treatment with ABA, an ABA receptor antagonist (AS6), or GA in the presence or absence of Fusarium graminearum (Fg) challenge. RESULTS: Spike treatments resulted in a total of 30,876 differentially expressed genes (DEGs) identified in 'Fielder' (26,004) and the Fg (4872) pathogen. Topology overlap and correlation analyses defined 9689 wheat DEGs as Fg-related across the treatments. Further enrichment analyses demonstrated that these included expression changes within 'Fielder' defense responses, cell structural metabolism, molecular transport, and membrane/lipid metabolism. Dysregulation of ABA and GA crosstalk arising from repression of 'Fielder' FUS3 was noted. As well, expression of a putative Fg ABA-biosynthetic cytochrome P450 was detected. The co-applied condition of Fg + ABA elicited further up-regulation of phytohormone biosynthesis, as well as SA and ET signaling pathways and cell wall/polyphenolic metabolism. In contrast, co-applied Fg + GA mainly suppressed phytohormone biosynthesis and signaling, while modulating primary and secondary metabolism and flowering. Unexpectedly, co-applied Fg + AS6 did not affect ABA biosynthesis or signaling, but rather elicited antagonistic responses tied to stress, phytohormone transport, and FHB disease-related genes. CONCLUSIONS: Observed exacerbation (misregulation) of classical defense mechanisms and cell wall fortifications upon ABA treatment are consistent with its ability to promote FHB severity and its proposed role as a fungal effector. In contrast, GA was found to modulate primary and secondary metabolism, suggesting a general metabolic shift underlying its reduction in FHB severity. While AS6 did not antagonize traditional ABA pathways, its impact on host defense and Fg responses imply potential for future investigation. Overall, by comparing these findings to those previously reported for four additional plant genotypes, an additive model of the wheat-Fg interaction is proposed in the context of phytohormone responses.

Ectopic expression of cDNAs from larkspur (Consolida ajacis) for increased synthesis of gondoic acid (cis-11 eicosenoic acid) and its positional redistribution in seed triacylglycerol of Camelina sativa.

MAIN CONCLUSION: A ß-ketoacyl-ACP-synthase II (KAS2) like enzyme and a lysophosphatidic acid acyltransferase (LPAT2) from Consolida ajacis catalyze gondoic acid biosynthesis and incorporation into the sn-2 position of seed TAG in engineered Camelina sativa. Gondoic acid (cis-11 eicosenoic acid, 20:1∆11) is the predominant very-long-chain fatty acid (VLCFA) in camelina (Camelina sativa) seed oil accounting for 12-15% of total triacylglycerol fatty acids. To explore the feasibility of engineering increased levels of this fatty acid in camelina seed, oils from a range of plant species were analyzed to identify those producing 20-Carbon (C20) fatty acids as the only VLCFAs in their seed oil. Seeds of Consolida and Delphinium species (Ranunculaceae) were found to contain moderate levels (0.2% to 25.5%) of C20 fatty acids without accompanying longer chain fatty acids. The C20 fatty acids were abundant in both sn-2 and sn-1/3 positions of seed TAG in Consolida, but were largely absent from the sn-2 position in Delphinium seed TAG. Through generation of a developing seed transcriptome, sequences were identified and cDNAs amplified from Consolida ajacis encoding a ß-ketoacyl-ACP-synthase II like protein (CaKAS2B) that lacked a predicted chloroplast transit peptide, and two homologues of Arabidopsis thaliana lysophosphatidic acid acyltransferase 2 (CaLPAT2a and CaLPAT2b). Expression of CaKAS2B in conventional (WT) camelina and a line previously engineered for high seed oleic acid content (HO) resulted in increased seed VLCFA content. Total VLCFA levels were raised from 24 to 35% and from 7 to 23% in T3 seed from representative transformants in the WT and HO backgrounds, respectively. Gondoic acid was the predominant VLCFA in transformed HO lines with low endogenous cytoplasmic fatty acid elongation activity, suggesting limited capacity of CaKAS2B to elongate beyond C20. Expression in camelina of CaLPAT2b resulted in significantly increased C20-VLCFA esterification at the sn-2 position of seed TAG with VLCFA levels of 33.8% in this position in one transformed line compared to 0.3% at sn-2 in the corresponding control line. Only small changes in total seed VLCFA content were observed in transformed lines implying that increased VLCFA esterification capacity in camelina results in positional redistribution of VLCFAs but does not significantly enhance flux through the fatty acid elongation pathway. The full potential of CaKAS2B and CaLPAT2a for the engineering of high gondoic acid levels in camelina remains to be determined. Seed fatty acid composition of Consolida and Delphinium also provides information that may be of value in the systematics of the Ranunculaceae.

Elucidating the biochemical basis of trans-16:1 fatty acid change in leaves during cold acclimation in wheat.

In plant cells, phosphatidylglycerol (PG) in the chloroplast has a characteristic trans-∆3-hexadecenoic acid (t16:1) at the sn-2 position. The t16:1 content in wheat leaf tissues decreases during cold treatment, but the significance of this fatty acid compositional change and the underlying biochemical mechanism remains poorly understood. Using a large collection of wheat cultivars displaying a varying capacity of freezing tolerance, we show for the first time under field conditions that this low temperature induced t16:1 change is associated with winter hardiness. To explore the metabolic mechanism responsible for the reduction of t16:1, we performed detailed lipid analysis and comparative transcriptome study with four selected wheat lines under cold acclimation. Our results show that wheat leaf tissues experience a gradual decrease in chloroplast lipid pathway activity during cold acclimation and that the decline in chloroplast lipid synthesis manifests itself in the decrease of t16:1 in PG. Comparative RNA-seq analyses with leaf tissues further reveal concerted transcriptome shifts indicating a rebalancing of chloroplast and cytosolic lipid synthesis during cold acclimation. Our study, thus, provides mechanistic understanding on chloroplast lipid adjustments as a "molecular ideotype" and the t16:1 change as a specific metabolite marker for screening freezing tolerance in wheat.

WheatCRISPR: a web-based guide RNA design tool for CRISPR/Cas9-mediated genome editing in wheat.

BACKGROUND: CRISPR/Cas9 gene editing has become a revolutionary technique for crop improvement as it can facilitate fast and efficient genetic changes without the retention of transgene components in the final plant line. Lack of robust bioinformatics tools to facilitate the design of highly specific functional guide RNAs (gRNAs) and prediction of off-target sites in wheat is currently an obstacle to effective application of CRISPR technology to wheat improvement. DESCRIPTION: We have developed a web-based bioinformatics tool to design specific gRNAs for genome editing and transcriptional regulation of gene expression in wheat. A collaborative study between the Broad Institute and Microsoft Research used large-scale empirical evidence to devise algorithms (Doech et al., 2016, Nature Biotechnology 34, 184-191) for predicting the on-target activity and off-target potential of CRISPR/SpCas9 (Streptococcus pyogenes Cas9). We applied these prediction models to determine on-target specificity and potential off-target activity for individual gRNAs targeting specific loci in the wheat genome. The genome-wide gRNA mappings and the corresponding Doench scores predictive of the on-target and off-target activities were used to create a gRNA database which was used as a data source for the web application termed WheatCRISPR. CONCLUSION: The WheatCRISPR tool allows researchers to browse all possible gRNAs targeting a gene or sequence of interest and select effective gRNAs based on their predicted high on-target and low off-target activity scores, as well as other characteristics such as position within the targeted gene. It is publicly available at https://crispr.bioinfo.nrc.ca/WheatCrispr/ .

Exogenous Abscisic Acid and Gibberellic Acid Elicit Opposing Effects on Fusarium graminearum Infection in Wheat.

Although the roles of salicylate (SA) and jasmonic acid (JA) have been well-characterized in Fusarium head blight (FHB)-infected cereals, the roles of other phytohormones remain more ambiguous. Here, the association between an array of phytohormones and FHB pathogenesis in wheat is investigated. Comprehensive profiling of endogenous hormones demonstrated altered cytokinin, gibberellic acid (GA), and JA metabolism in a FHB-resistant cultivar, whereas challenge by Fusarium graminearum increased abscisic acid (ABA), JA, and SA in both FHB-susceptible and -resistant cultivars. Subsequent investigation of ABA or GA coapplication with fungal challenge increased and decreased FHB spread, respectively. These phytohormones-induced effects may be attributed to alteration of the F. graminearum transcriptome because ABA promoted expression of early-infection genes, including hydrolases and cytoskeletal reorganization genes, while GA suppressed nitrogen metabolic gene expression. Neither ABA nor GA elicited significant effects on F. graminearum fungal growth or sporulation in axenic conditions, nor do these phytohormones affect trichothecene gene expression, deoxynivalenol mycotoxin accumulation, or SA/JA biosynthesis in F. graminearum-challenged wheat spikes. Finally, the combined application of GA and paclobutrazol, a Fusarium fungicide, provided additive effects on reducing FHB severity, highlighting the potential for combining fungicidal agents with select phytohormone-related treatments for management of FHB infection in wheat.

Understanding the biochemical basis of temperature-induced lipid pathway adjustments in plants.

Glycerolipid biosynthesis in plants proceeds through two major pathways compartmentalized in the chloroplast and the endoplasmic reticulum (ER). The involvement of glycerolipid pathway interactions in modulating membrane desaturation under temperature stress has been suggested but not fully explored. We profiled glycerolipid changes as well as transcript dynamics under suboptimal temperature conditions in three plant species that are distinctively different in the mode of lipid pathway interactions. In Arabidopsis thaliana, a 16:3 plant, the chloroplast pathway is upregulated in response to low temperature, whereas high temperature promotes the eukaryotic pathway. Operating under a similar mechanistic framework, Atriplex lentiformis at high temperature drastically increases the contribution of the eukaryotic pathway and correspondingly suppresses the prokaryotic pathway, resulting in the switch of lipid profile from 16:3 to 18:3. In wheat (Triticum aestivum), an 18:3 plant, low temperature also influences the channeling of glycerolipids from the ER to chloroplast. Evidence of differential trafficking of diacylglycerol moieties from the ER to chloroplast was uncovered in three plant species as another layer of metabolic adaptation under temperature stress. We propose a model that highlights the predominance and prevalence of lipid pathway interactions in temperature-induced lipid compositional changes.

Transcriptome Profiling Identifies Candidate Genes Associated with the Accumulation of Distinct Sulfur γ-Glutamyl Dipeptides in Phaseolus vulgaris and Vigna mungo Seeds.

Common bean (Phaseolus vulgaris) and black gram (Vigna mungo) accumulate γ-Glutamyl-S-methylcysteine and γ-Glutamyl-methionine in seed, respectively. Transcripts were profiled by 454 pyrosequencing data at a similar developmental stage coinciding with the beginning of the accumulation of these metabolites. Expressed sequence tags were assembled into Unigenes, which were assigned to specific genes in the early release chromosomal assembly of the P. vulgaris genome. Genes involved in multiple sulfur metabolic processes were expressed in both species. Expression of Sultr3 members was predominant in P. vulgaris, whereas expression of Sultr5 members predominated in V. mungo. Expression of the cytosolic SERAT1;1 and -1;2 was approximately fourfold higher in P. vulgaris while expression of the plastidic SERAT2;1 was twofold higher in V. mungo. Among BSAS family members, BSAS4;1, encoding a cytosolic cysteine desulfhydrase, and BSAS1;1, encoding a cytosolic O-acetylserine sulphydrylase were most highly expressed in both species. This was followed by BSAS3;1 encoding a plastidic ß-cyanoalanine synthase which was more highly expressed by 10-fold in P. vulgaris. The data identify BSAS3;1 as a candidate enzyme for the biosynthesis of S-methylcysteine through the use of methanethiol as substrate instead of cyanide. Expression of GLC1 would provide a complete sequence leading to the biosynthesis of γ-Glutamyl-S-methylcysteine in plastids. The detection of S-methylhomoglutathione in P. vulgaris suggested that homoglutathione synthetase may accept, to some extent, γ-Glutamyl-S-methylcysteine as substrate, which might lead to the formation of S-methylated phytochelatins. In conclusion, 454 sequencing was effective at revealing differences in the expression of sulfur metabolic genes, providing information on candidate genes for the biosynthesis of distinct sulfur amino acid γ-Glutamyl dipeptides between P. vulgaris and V. mungo.

Fatty acid composition of developing sea buckthorn (Hippophae rhamnoides L.) berry and the transcriptome of the mature seed.

BACKGROUND: Sea buckthorn (Hippophae rhamnoides L.) is a hardy, fruit-producing plant known historically for its medicinal and nutraceutical properties. The most recognized product of sea buckthorn is its fruit oil, composed of seed oil that is rich in essential fatty acids, linoleic (18:2 ω-6) and α-linolenic (18:3 ω-3) acids, and pulp oil that contains high levels of monounsaturated palmitoleic acid (16:1 ω-7). Sea buckthorn is fast gaining popularity as a source of functional food and nutraceuticals, but currently has few genomic resources; therefore, we explored the fatty acid composition of Canadian-grown cultivars (ssp. mongolica) and the sea buckthorn seed transcriptome using the 454 GS FLX sequencing technology. RESULTS: GC-MS profiling of fatty acids in seeds and pulp of berries indicated that the seed oil contained linoleic and α-linolenic acids at 33-36% and 30-36%, respectively, while the pulp oil contained palmitoleic acid at 32-42%. 454 sequencing of sea buckthorn cDNA collections from mature seeds yielded 500,392 sequence reads, which identified 89,141 putative unigenes represented by 37,482 contigs and 51,659 singletons. Functional annotation by Gene Ontology and computational prediction of metabolic pathways indicated that primary metabolism (protein>nucleic acid>carbohydrate>lipid) and fatty acid and lipid biosynthesis pathways were highly represented categories. Sea buckthorn sequences related to fatty acid biosynthesis genes in Arabidopsis were identified, and a subset of these was examined for transcript expression at four developing stages of the berry. CONCLUSION: This study provides the first comprehensive genomic resources represented by expressed sequences for sea buckthorn, and demonstrates that the seed oil of Canadian-grown sea buckthorn cultivars contains high levels of linoleic acid and α-linolenic acid in a close to 1:1 ratio, which is beneficial for human health. These data provide the foundation for further studies on sea buckthorn oil, the enzymes involved in its biosynthesis, and the genes involved in the general hardiness of sea buckthorn against environmental conditions.

Integration of deep transcript and targeted metabolite profiles for eight cultivars of opium poppy.

Recent advances in DNA sequencing technology and analytical mass spectrometry are providing unprecedented opportunities to develop the functional genomics resources required to investigate complex biological processes in non-model plants. Opium poppy produces a wide variety of benzylisoquinoline alkaloids (BIAs), including the pharmaceutical compounds codeine, morphine, noscapine and papaverine. A functional genomics platform to identify novel BIA biosynthetic and regulatory genes in opium poppy has been established based on the differential metabolite profile of eight selected cultivars. Stem cDNA libraries from each of the eight opium poppy cultivars were subjected to 454 pyrosequencing and searchable expressed sequence tag databases were created from the assembled reads. These deep and integrated metabolite and transcript databases provide a nearly complete representation of the genetic and metabolic variances responsible for the differential occurrence of specific BIAs in each cultivar as demonstrated using the biochemically well characterized pathway from tyrosine to morphine. Similar correlations between the occurrence of specific transcripts and alkaloids effectively reveals candidate genes encoding uncharacterized biosynthetic enzymes as shown using cytochromes P450 potentially involved in the formation of papaverine and noscapine.

Gene expression analysis of flax seed development.

BACKGROUND: Flax, Linum usitatissimum L., is an important crop whose seed oil and stem fiber have multiple industrial applications. Flax seeds are also well-known for their nutritional attributes, viz., omega-3 fatty acids in the oil and lignans and mucilage from the seed coat. In spite of the importance of this crop, there are few molecular resources that can be utilized toward improving seed traits. Here, we describe flax embryo and seed development and generation of comprehensive genomic resources for the flax seed. RESULTS: We describe a large-scale generation and analysis of expressed sequences in various tissues. Collectively, the 13 libraries we have used provide a broad representation of genes active in developing embryos (globular, heart, torpedo, cotyledon and mature stages) seed coats (globular and torpedo stages) and endosperm (pooled globular to torpedo stages) and genes expressed in flowers, etiolated seedlings, leaves, and stem tissue. A total of 261,272 expressed sequence tags (EST) (GenBank accessions LIBEST_026995 to LIBEST_027011) were generated. These EST libraries included transcription factor genes that are typically expressed at low levels, indicating that the depth is adequate for in silico expression analysis. Assembly of the ESTs resulted in 30,640 unigenes and 82% of these could be identified on the basis of homology to known and hypothetical genes from other plants. When compared with fully sequenced plant genomes, the flax unigenes resembled poplar and castor bean more than grape, sorghum, rice or Arabidopsis. Nearly one-fifth of these (5,152) had no homologs in sequences reported for any organism, suggesting that this category represents genes that are likely unique to flax. Digital analyses revealed gene expression dynamics for the biosynthesis of a number of important seed constituents during seed development. CONCLUSIONS: We have developed a foundational database of expressed sequences and collection of plasmid clones that comprise even low-expressed genes such as those encoding transcription factors. This has allowed us to delineate the spatio-temporal aspects of gene expression underlying the biosynthesis of a number of important seed constituents in flax. Flax belongs to a taxonomic group of diverse plants and the large sequence database will allow for evolutionary studies as well.

Integration of deep transcriptome and proteome analyses reveals the components of alkaloid metabolism in opium poppy cell cultures.

BACKGROUND: Papaver somniferum (opium poppy) is the source for several pharmaceutical benzylisoquinoline alkaloids including morphine, the codeine and sanguinarine. In response to treatment with a fungal elicitor, the biosynthesis and accumulation of sanguinarine is induced along with other plant defense responses in opium poppy cell cultures. The transcriptional induction of alkaloid metabolism in cultured cells provides an opportunity to identify components of this process via the integration of deep transcriptome and proteome databases generated using next-generation technologies. RESULTS: A cDNA library was prepared for opium poppy cell cultures treated with a fungal elicitor for 10 h. Using 454 GS-FLX Titanium pyrosequencing, 427,369 expressed sequence tags (ESTs) with an average length of 462 bp were generated. Assembly of these sequences yielded 93,723 unigenes, of which 23,753 were assigned Gene Ontology annotations. Transcripts encoding all known sanguinarine biosynthetic enzymes were identified in the EST database, 5 of which were represented among the 50 most abundant transcripts. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) of total protein extracts from cell cultures treated with a fungal elicitor for 50 h facilitated the identification of 1,004 proteins. Proteins were fractionated by one-dimensional SDS-PAGE and digested with trypsin prior to LC-MS/MS analysis. Query of an opium poppy-specific EST database substantially enhanced peptide identification. Eight out of 10 known sanguinarine biosynthetic enzymes and many relevant primary metabolic enzymes were represented in the peptide database. CONCLUSIONS: The integration of deep transcriptome and proteome analyses provides an effective platform to catalogue the components of secondary metabolism, and to identify genes encoding uncharacterized enzymes. The establishment of corresponding transcript and protein databases generated by next-generation technologies in a system with a well-defined metabolite profile facilitates an improved linkage between genes, enzymes, and pathway components. The proteome database represents the most relevant alkaloid-producing enzymes, compared with the much deeper and more complete transcriptome library. The transcript database contained full-length mRNAs encoding most alkaloid biosynthetic enzymes, which is a key requirement for the functional characterization of novel gene candidates.

Exploiting the wild crucifer Thlaspi arvense to identify conserved and novel genes expressed during a plant's response to cold stress.

Thlaspi arvense, a wild species from the Brassicaceae family, was shown to have a higher level of freezing tolerance than either of its close relatives, the model plant Arabidopsis thaliana or the crop Brassica napus (canola). Over 600 clones were sequenced from a subtractive cDNA library generated from cold treated T. arvense tissue, establishing that T. arvense shared significant sequence identity with both A. thaliana and B. napus (90-92%). In light of the strong sequence similarity between T. arvense and A. thaliana and to exploit the available genomics resources for Arabidopsis, the efficacy of using long 70 mer oligonucleotide whole genome Arabidopsis microarrays was tested for T. arvense. Gene expression in T. arvense leaf tissue during the very early stages of cold acclimation (or cold stress) was assayed at three time points and compared to an untreated control. This analysis highlights some of the difficulties and benefits of using cross-species microarray analysis. The data suggested that T. arvense responds in a similar fashion to cold stress as the model plant A. thaliana. However, for a number of genes quantitative differences in the level and timing of expression were identified. One of the most notable differences suggested that sulphur assimilation leading to the increased production of the methyl donor S-adenosyl-methionine was playing a role in the response of T. arvense to cold stress.

Maximizing the efficacy of SAGE analysis identifies novel transcripts in Arabidopsis.

The efficacy of using Serial Analysis of Gene Expression (SAGE) to analyze the transcriptome of the model dicotyledonous plant Arabidopsis was assessed. We describe an iterative tag-to-gene matching process that exploits the availability of the whole genome sequence of Arabidopsis. The expression patterns of 98% of the annotated Arabidopsis genes could theoretically be evaluated through SAGE and using an iterative matching process 79% could be identified by a tag found at a unique site in the genome. A total of 145,170 reliable experimental tags from two Arabidopsis leaf tissue SAGE libraries were analyzed, of which 29,632 were distinct. The majority (93%) of the 12,988 experimental tags observed greater than once could be matched within the Arabidopsis genome. However, only 78% were matched to a single locus within the genome, reflecting the complexities associated with working in a highly duplicated genome. In addition to a comprehensive assessment of gene expression in Arabidopsis leaf tissue, we describe evidence of transcription from pseudo-genes as well as evidence of alternative mRNA processing and anti-sense transcription. This collection of experimental SAGE tags could be exploited to assist in the on-going annotation of the Arabidopsis genome.