Identification and analysis of safener-inducible expressed sequence tags in Populus using a cDNA microarray

Transcriptome Profiling and Genome-Wide Association Studies Reveal GSTs and Other Defense Genes Involved in Multiple Signaling Pathways Induced by Herbicide Safener in Grain Sorghum

Herbicide safeners protect cereal crops from herbicide injury by inducing genes and proteins involved in detoxification reactions, such as glutathione S-transferases (GSTs) and cytochrome P450s (P450s). Only a few studies have characterized gene or protein expression profiles for investigating plant responses to safener treatment in cereal crops, and most transcriptome analyses in response to safener treatments have been conducted in dicot model species that are not protected by safener from herbicide injury. In this study, three different approaches were utilized in grain sorghum (Sorghum bicolor (L.) Moench) to investigate mechanisms involved in safener-regulated signaling pathways. An initial transcriptome analysis was performed to examine global gene expression in etiolated shoot tissues of hybrid grain sorghum following treatment with the sorghum safener, fluxofenim. Most upregulated transcripts encoded detoxification enzymes, including P450s, GSTs, and UDP-dependent glucosyltransferases (UGTs). Interestingly, several of these upregulated transcripts are similar to genes involved with the biosynthesis and recycling/catabolism of dhurrin, an important chemical defense compound, in these seedling tissues. Secondly, 761 diverse sorghum inbred lines were evaluated in a genome-wide association study (GWAS) to determine key molecular-genetic factors governing safener-mediated signaling mechanisms and/or herbicide detoxification. GWAS revealed a significant single nucleotide polymorphism (SNP) associated with safener-induced response on chromosome 9, located within a phi-class SbGST gene and about 15-kb from a different phi-class SbGST. Lastly, the expression of these two candidate SbGSTs was quantified in etiolated shoot tissues of sorghum inbred BTx623 in response to fluxofenim treatment. SbGSTF1 and SbGSTF2 transcripts increased within 12-hr after fluxofenim treatment but the level of safener-induced expression differed between the two genes. In addition to identifying specific GSTs potentially involved in the safener-mediated detoxification pathway, this research elucidates a new direction for studying both constitutive and inducible mechanisms for chemical defense in cereal crop seedlings.

Genetic manipulation of putrescine biosynthesis reprograms the cellular transcriptome and the metabolome

BACKGROUND

With the increasing interest in metabolic engineering of plants using genetic manipulation and gene editing technologies to enhance growth, nutritional value and environmental adaptation, a major concern is the potential of undesirable broad and distant effects of manipulating the target gene or metabolic step in the resulting plant. A comprehensive transcriptomic and metabolomic analysis of the product may shed some useful light in this regard. The present study used these two techniques with plant cell cultures to analyze the effects of genetic manipulation of a single step in the biosynthesis of polyamines because of their well-known roles in plant growth, development and stress responses.

RESULTS

The transcriptomes and metabolomes of a control and a high putrescine (HP) producing cell line of poplar (Populus nigra x maximowiczii) were compared using microarrays and GC/MS. The HP cells expressed an ornithine decarboxylase transgene and accumulated several-fold higher concentrations of putrescine, with only small changes in spermidine and spermine. The results show that up-regulation of a single step in the polyamine biosynthetic pathway (i.e. ornithine → putrescine) altered the expression of a broad spectrum of genes; many of which were involved in transcription, translation, membrane transport, osmoregulation, shock/stress/wounding, and cell wall metabolism. More than half of the 200 detected metabolites were significantly altered (p ≤ 0.05) in the HP cells irrespective of sampling date. The most noteworthy differences were in organic acids, carbohydrates and nitrogen-containing metabolites.

CONCLUSIONS

The results provide valuable information about the role of polyamines in regulating nitrogen and carbon use pathways in cell cultures of high putrescine producing transgenic cells of poplar vs. their low putrescine counterparts. The results underscore the complexity of cellular responses to genetic perturbation of a single metabolic step related to nitrogen metabolism in plants. Combined with recent studies from our lab, where we showed that higher putrescine production caused an increased flux of glutamate into ornithine concurrent with enhancement in glutamate production via additional nitrogen and carbon assimilation, the results from this study provide guidance in designing transgenic plants with increased nitrogen use efficiency, especially in plants intended for non-food/feed applications (e.g. increased biomass production for biofuels).

Transgenic poplar trees expressing yeast cadmium factor 1 exhibit the characteristics necessary for the phytoremediation of mine tailing soil

Genetic engineering of plants for phytoremediation is thought to be possible based on results using model plants expressing genes involved in heavy metal resistance, which improve the plant's tolerance of heavy metals and accumulation capacity. The next step of progress in this technology requires the genetic engineering of plants that produce large amounts of biomass and the testing of these transgenic plants in contaminated soils. Thus, we transformed a sterile line of poplar Populus alba X P. tremula var. glandulosa with a heavy metal resistance gene, ScYCF1 (yeast cadmium factor 1), which encodes a transporter that sequesters toxic metal(loid)s into the vacuoles of budding yeast, and tested these transgenic plants in soil taken from a closed mine site contaminated with multiple toxic metal(loid)s under greenhouse and field conditions. The YCF1-expressing transgenic poplar plants exhibited enhanced growth, reduced toxicity symptoms, and increased Cd content in the aerial tissue compared to the non-transgenic plants. Furthermore, the plants accumulated increased amounts of Cd, Zn, and Pb in the root, because they could establish an extensive root system in mine tailing soil. These results suggest that the generation of YCF1-expressing transgenic poplar represents the first step towards producing plants for phytoremediation. The YCF1-expressing poplar may be useful for phytostabilization and phytoattenuation, especially in highly contaminated regions, where wild-type plants cannot survive.

Microarray and suppression subtractive hybridization analyses of gene expression in hybrid poplar (Populus alba × Populus tremula var. glandulosa) cell suspension cultures after exposure to NaCl

The gene expression profiles of hybrid poplar (Populus alba × Populus tremula var. glandulosa) cells in suspension culture after exposure to salinity (NaCl) induced stress were examined by constructing two suppression subtractive hybridization (SSH) libraries. cDNA from non-treated cells was used as a driver and cDNA samples from cell suspension cultures exposed to 150 mM NaCl for 2 or 10 h were used as testers. Randomly selected clones from each SSH library were sequenced and 727 high-quality expressed sequence tags (ESTs) were obtained and analyzed. Four novel ESTs were identified. Between the two libraries, 542 unique SSH clones were selected for placement on a cDNA microarray. In total, 18 differentially expressed genes were identified with 4 and 12 genes being significantly differentially expressed 2 and 10 h after the treatment, respectively. Genes related to metabolism and protein synthesis and several genes whose protein products are implicated in salt or other abiotic stress-related responses were expressed in the salt-stressed cells.

Safeners recruit multiple signalling pathways for the orchestrated induction of the cellular xenobiotic detoxification machinery in Arabidopsis

Safeners enhance herbicide tolerance in crop plants but not in target weeds, thus improving herbicide selectivity. The safeners isoxadifen-ethyl and mefenpyr-diethyl protect cereal crops from sulfonyl urea herbicides in postemergence application. The two safeners were shown here to induce the cellular xenobiotic detoxification machinery in Arabidopsis thaliana when applied to leaves in a way mimicking field application. Gene expression profiling revealed the induction of 446 genes potentially involved in the detoxification process. Transgenic Arabidopsis plants expressing a reporter gene under control of a safener-responsive maize promoter were used as a model system to study the safener signalling pathway. Reporter gene analysis in the tga2/3/5/6, sid2-2 and npr1 mutants as compared with the wild-type background showed that safener inducibility required TGA transcription factors and salicylic acid (SA) in a NON-EXPRESSOR of PR-1 (NPR1)-independent pathway converging on two as-1 promoter elements. For the majority of the safener-responsive Arabidopsis genes, a similar dependence on TGA transcription factors and/or SA was shown by gene expression profiling in wild-type plants as compared with the tga2/3/5/6 and sid2-2 mutants. Thirty-eight percent of the genes, however, were induced by safeners in a TGA/SA-independent manner. These genes are likely to be controlled by WRKY transcription factors and cognate W-boxes in their promoters.

Xenobiotic responsiveness of Arabidopsis thaliana to a chemical series derived from a herbicide safener

Plants respond to synthetic chemicals by eliciting a xenobiotic response (XR) that enhances the expression of detoxifying enzymes such as glutathione transferases (GSTs). In agrochemistry, the ability of safeners to induce an XR is used to increase herbicide detoxification in cereal crops. Based on the responsiveness of the model plant Arabidopsis thaliana to the rice safener fenclorim (4,6-dichloro-2-phenylpyrimidine), a series of related derivatives was prepared and tested for the ability to induce GSTs in cell suspension cultures. The XR in Arabidopsis could be divided into rapid and slow types depending on subtle variations in the reactivity (electrophilicity) and chemical structure of the derivatives. In a comparative microarray study, Arabidopsis cultures were treated with closely related compounds that elicited rapid (fenclorim) and slow (4-chloro-6-methyl-2-phenylpyrimidine) XRs. Both chemicals induced major changes in gene expression, including a coordinated suppression in cell wall biosynthesis and an up-regulation in detoxification pathways, whereas only fenclorim selectively induced sulfur and phenolic metabolism. These transcriptome studies suggested several linkages between the XR and oxidative and oxylipin signaling. Confirming links with abiotic stress signaling, suppression of glutathione content enhanced GST induction by fenclorim, whereas fatty acid desaturase mutants, which were unable to synthesize oxylipins, showed an attenuated XR. Examining the significance of these studies to agrochemistry, only those fenclorim derivatives that elicited a rapid XR proved effective in increasing herbicide tolerance (safening) in rice.

Living with high putrescine: expression of ornithine and arginine biosynthetic pathway genes in high and low putrescine producing poplar cells

Arginine (Arg) and ornithine (Orn), both derived from glutamate (Glu), are the primary substrates for polyamine (PA) biosynthesis, and also play important roles as substrates and intermediates of overall N metabolism in plants. Their cellular homeostasis is subject to multiple levels of regulation. Using reverse transcription quantitative PCR (RT-qPCR), we studied changes in the expression of all genes of the Orn/Arg biosynthetic pathway in response to up-regulation [via transgenic expression of mouse Orn decarboxylase (mODC)] of PA biosynthesis in poplar (Populus nigra × maximowiczii) cells grown in culture. Cloning and sequencing of poplar genes involved in the Orn/Arg biosynthetic pathway showed that they have high homology with similar genes in other plants. The expression of the genes of Orn, Arg and PA biosynthetic pathway fell into two hierarchical clusters; expression of one did not change in response to high putrescine, while members of the other cluster showed a shift in expression pattern during the 7-day culture cycle. Gene expression of branch point enzymes (N-acetyl-Glu synthase, Orn aminotransferase, Arg decarboxylase, and spermidine synthase) in the sub-pathways, constituted a separate cluster from those involved in intermediary reactions of the pathway (N-acetyl-Glu kinase, N-acetyl-Glu-5-P reductase, N-acetyl-Orn aminotransferase, N (2)-acetylOrn:N-acetyl-Glu acetyltransferase, N (2)-acetyl-Orn deacetylase, Orn transcarbamylase, argininosuccinate synthase, carbamoylphosphate synthetase, argininosuccinate lyase, S-adenosylmethionine decarboxylase, spermine synthase). We postulate that expression of all genes of the Glu-Orn-Arg pathway is constitutively coordinated and is not influenced by the increase in flux rate through this pathway in response to increased utilization of Orn by mODC; thus the pathway involves mostly biochemical regulation rather than changes in gene expression. We further suggest that Orn itself plays a major role in the regulation of this pathway.

Analysis of 4,664 high-quality sequence-finished poplar full-length cDNA clones and their utility for the discovery of genes responding to insect feeding

Background

The genus Populus includes poplars, aspens and cottonwoods, which will be collectively referred to as poplars hereafter unless otherwise specified. Poplars are the dominant tree species in many forest ecosystems in the Northern Hemisphere and are of substantial economic value in plantation forestry. Poplar has been established as a model system for genomics studies of growth, development, and adaptation of woody perennial plants including secondary xylem formation, dormancy, adaptation to local environments, and biotic interactions.

Results

As part of the poplar genome sequencing project and the development of genomic resources for poplar, we have generated a full-length (FL)-cDNA collection using the biotinylated CAP trapper method. We constructed four FLcDNA libraries using RNA from xylem, phloem and cambium, and green shoot tips and leaves from the P. trichocarpa Nisqually-1 genotype, as well as insect-attacked leaves of the P. trichocarpa × P. deltoides hybrid. Following careful selection of candidate cDNA clones, we used a combined strategy of paired end reads and primer walking to generate a set of 4,664 high-accuracy, sequence-verified FLcDNAs, which clustered into 3,990 putative unique genes. Mapping FLcDNAs to the poplar genome sequence combined with BLAST comparisons to previously predicted protein coding sequences in the poplar genome identified 39 FLcDNAs that likely localize to gaps in the current genome sequence assembly. Another 173 FLcDNAs mapped to the genome sequence but were not included among the previously predicted genes in the poplar genome. Comparative sequence analysis against Arabidopsis thaliana and other species in the non-redundant database of GenBank revealed that 11.5% of the poplar FLcDNAs display no significant sequence similarity to other plant proteins. By mapping the poplar FLcDNAs against transcriptome data previously obtained with a 15.5 K cDNA microarray, we identified 153 FLcDNA clones for genes that were differentially expressed in poplar leaves attacked by forest tent caterpillars.

Conclusion

This study has generated a high-quality FLcDNA resource for poplar and the third largest FLcDNA collection published to date for any plant species. We successfully used the FLcDNA sequences to reassess gene prediction in the poplar genome sequence, perform comparative sequence annotation, and identify differentially expressed transcripts associated with defense against insects. The FLcDNA sequences will be essential to the ongoing curation and annotation of the poplar genome, in particular for targeting gaps in the current genome assembly and further improvement of gene predictions. The physical FLcDNA clones will serve as useful reagents for functional genomics research in areas such as analysis of gene functions in defense against insects and perennial growth. Sequences from this study have been deposited in NCBI GenBank under the accession numbers EF144175 to EF148838.

Functional annotation of 19,841 Populus nigra full-length enriched cDNA clones

Background

Populus is one of favorable model plants because of its small genome. Structural genomics of Populus has reached a breakpoint as nucleotides of the entire genome have been determined. Reaching the post genome era, functional genomics of Populus is getting more important for well-comprehended plant science. Development of bioresorce serving functional genomics is making rapid progress. Huge efforts have achieved deposits of expressed sequence tags (ESTs) in various plant species consequently accelerating functional analysis of genes. ESTs from full-length cDNA clones are especially powerful for accurate molecular annotation. We promoted collection and annotation of the ESTs from Populus full-length enriched cDNA clones as part of functional genomics of tree species.

Results

We have been collecting the full-length enriched cDNA of the female poplar (Populus nigra var. italica) for years. By sequencing P. nigra full-length (PnFL) cDNA libraries, we generated about 116,000 5'-end or 3'-end ESTs corresponding to 19,841 nonredundant PnFL clones. Population of PnFL cDNA clones represents 44% of the predicted genes in the Populus genome.

Conclusion

Our resource of P. nigra full-length enriched clones is expected to provide valuable tools to gain further insight into genome annotation and functional genomics in Populus.

Differentially expressed cDNAs at the early stage of banana ripening identified by suppression subtractive hybridization and cDNA microarray

The banana (Musa acuminate L. AAA group) fruit undergoes a postharvest ripening process, which plays an important role in improving the quality and extending the shelf life of bananas. To manipulate postharvest banana ripening, a better understanding of the mechanism of postharvest ripening is necessary. The isolation of mRNA transcripts encoding proteins associated with the ripening process is a powerful tool for this purpose. To isolate differentially expressed genes at the early stage of postharvest banana ripening, a forward suppression subtractive hybridization (SSH) cDNA library was constructed. SSH was performed with cDNA from banana fruit on the day of harvest as the "driver" and cDNA from banana fruit 2 days postharvest (DPH) as the "tester." A total of 289 clones in the SSH library were sequenced. BLASTX results revealed that 191 cDNAs had significant sequence homologies with known sequences in the NCBI database. Of the 191 cDNAs, 138 were singletons, and 53 belonged to divergent clusters containing 2-8 sequences. The identified cDNAs encoded proteins involved in cellular processes such as: metabolism; protein destination and storage; protein synthesis; signal transduction; transport and intracellular traffic; cell structure, growth, and division; transcription and post-transcription; and disease and defense. To characterize differentially expressed cDNAs in the SSH library, cDNA microarray analysis was conducted. A total of 26 cDNAs in the 2-DPH banana fruit were found to be up-regulated and these results were confirmed by using reverse transcriptase-polymerase chain reaction (RT-PCR). The information generated in this study provides new clues to aid in the understanding of banana ripening.

Identification of regulatory pathways involved in the reacquisition of root growth after salt stress in Medicago truncatula

Root growth and function are determined by the action of environmental stresses through specific genes that adapt root development to these restrictive conditions. We have defined in vitro conditions affecting the growth and recovery of Medicago truncatula roots after a salt stress. A dedicated macroarray containing 384 genes, based on a large-scale subtractive hybridization approach, was constructed and used to analyze gene expression during salt stress and recovery of root growth from this stress. Several potential regulatory genes were identified as being linked to this recovery process: a novel RNA-binding protein, a small G-protein homologous to ROP9, a receptor-like kinase, two TF IIIA-like and an AP2-like transcription factors (TF), MtZpt2-1, MtZpt2-2 and MtAp2, and a histidine kinase associated with cytokinin transduction pathways. The two ZPT2-type TFs were also rapidly induced by cold stress in roots. By analyzing transgenic M. truncatula plants showing reduced expression levels of both TFs and affected in their capacity to recover root growth after a salt stress, we identified potential target genes that were either activated or repressed in these plants. Overexpression of MtZpt2-1 in roots conferred salt tolerance and affected the expression of three putative targets in the predicted manner: a cold-regulated A (CORA) homolog, a flower-promoting factor (FPF1) homolog and an auxin-induced proline-rich protein (PRP) gene. Hence, regulatory networks depending on TFIIIA-like transcription factors are involved in the control of root adaptation to salt stress.

Safeners coordinately induce the expression of multiple proteins and MRP transcripts involved in herbicide metabolism and detoxification inTriticum tauschii seedling tissues

Chemicals called safeners protect cereal crops from herbicide toxicity. Proteomic methods (2-D PAGE and LC-MS/MS) were utilized to identify safener- and/or herbicide-regulated proteins in three tissues (root, leaf, and coleoptile) of Triticum tauschii seedlings to better understand a safener's mechanism of action. Growth experiments showed that the safener cloquintocet-mexyl protected seedlings from injury by the herbicide dimethenamid. In total, 29 safener-induced and 10 herbicide-regulated proteins were identified by LC-MS/MS. These proteins were classified into two major categories based on their expression patterns, and were further classified into several functional groups. Surprisingly, mutually exclusive sets of proteins were identified following herbicide or safener treatment, suggesting that different signaling pathways may be recruited. Safener-responsive proteins, mostly involved in xenobiotic detoxification, also included several new proteins that had not been previously identified as safener-responsive, whereas herbicide-regulated proteins belonged to several classes involved in general stress responses. Quantitative RT-PCR revealed that multidrug resistance-associated protein (MRP) transcripts were highly induced by safeners and two MRP genes were differentially expressed. Our results indicate that safeners protect T. tauschii seedlings from herbicide toxicity by coordinately inducing proteins involved in an entire herbicide detoxification pathway mainly in the coleoptile and root, thereby protecting new leaves from herbicide injury.

UPSC-BASE--Populus transcriptomics online

The increasing accessibility and use of microarrays in transcriptomics has accentuated the need for purpose-designed storage and analysis tools. Here we present UPSC-BASE, a database for analysis and storage of Populus DNA microarray data. A microarray analysis pipeline has also been established to allow consistent and efficient analysis (from small to large scale) of samples in various experimental designs. A range of optimized experimental protocols is provided for each step in generating the data. Within UPSC-BASE, researchers can perform standard and advanced microarray analysis procedures in a user-friendly environment. Background corrections, normalizations, quality-control tools, visualizations, hypothesis tests and export tools are provided without requirements for expert-level knowledge. Although the database has been developed primarily for handling Populus DNA microarrays, most of the tools are generic and can be used for other types of microarray. UPSC-BASE is also a repository of Populus microarray information, providing data from 21 experiments on a total of 407 microarray hybridizations in the public domain of the database. There are also an additional 10 experiments containing 347 hybridizations, where the automatically analysed data are searchable.

Genomics of hybrid poplar (Populus trichocarpa× deltoides) interacting with forest tent caterpillars (Malacosoma disstria): normalized and full-length cDNA libraries, expressed sequence tags, and a cDNA microarray for the study of insect-induced defences

As part of a genomics strategy to characterize inducible defences against insect herbivory in poplar, we developed a comprehensive suite of functional genomics resources including cDNA libraries, expressed sequence tags (ESTs) and a cDNA microarray platform. These resources are designed to complement the existing poplar genome sequence and poplar (Populus spp.) ESTs by focusing on herbivore- and elicitor-treated tissues and incorporating normalization methods to capture rare transcripts. From a set of 15 standard, normalized or full-length cDNA libraries, we generated 139,007 3'- or 5'-end sequenced ESTs, representing more than one-third of the c. 385,000 publicly available Populus ESTs. Clustering and assembly of 107,519 3'-end ESTs resulted in 14,451 contigs and 20,560 singletons, altogether representing 35,011 putative unique transcripts, or potentially more than three-quarters of the predicted c. 45,000 genes in the poplar genome. Using this EST resource, we developed a cDNA microarray containing 15,496 unique genes, which was utilized to monitor gene expression in poplar leaves in response to herbivory by forest tent caterpillars (Malacosoma disstria). After 24 h of feeding, 1191 genes were classified as up-regulated, compared to only 537 down-regulated. Functional classification of this induced gene set revealed genes with roles in plant defence (e.g. endochitinases, Kunitz protease inhibitors), octadecanoid and ethylene signalling (e.g. lipoxygenase, allene oxide synthase, 1-aminocyclopropane-1-carboxylate oxidase), transport (e.g. ABC proteins, calreticulin), secondary metabolism [e.g. polyphenol oxidase, isoflavone reductase, (-)-germacrene D synthase] and transcriptional regulation [e.g. leucine-rich repeat transmembrane kinase, several transcription factor classes (zinc finger C3H type, AP2/EREBP, WRKY, bHLH)]. This study provides the first genome-scale approach to characterize insect-induced defences in a woody perennial providing a solid platform for functional investigation of plant-insect interactions in poplar.

Transcriptome analysis of ozone-responsive genes in leaves of European beech (Fagus sylvatica L.)

Suppression subtractive hybridization (SSH) was performed to isolate cDNAs representing genes that are differentially expressed in leaves of Fagus sylvatica upon ozone exposure. 1248 expressed sequence tags (ESTs) were obtained from 2 subtractive libraries containing early and late ozone-responsive genes. Sequences of 1139 clones (91 %) matched the EBI/NCBI database entries. For 578 clones, no putative function could be assigned. Most abundant transcripts were O-methyltransferases, representing 7 % of all sequenced clones. ESTs were organized into 12 functional categories according to the MIPS database. Among them, 12 % (early)/15 % (late) were associated with disease and defence, 19/11 % with cell structure, 4/10 % with signal transduction, and 9/6 % with transcription. The expression pattern of selected ESTs (ribulose-1,5-bisphosphate carboxylase/oxygenase small subunit [rbcS], WRKY-type transcription factor, ultraviolet-B-repressible protein, aquaporine, glutathione S-transferase, catalase, caffeic acid O-methyltransferase, and pathogenesis-related protein 1 [PR1]) was analysed by quantitative real-time RT-PCR (qRT-PCR) which confirmed changed transcript levels upon ozone treatment of European beech saplings. The ESTs characterized will contribute to a better understanding of forest tree genomics and also to a comparison of ozone-responsive genes in woody and herbaceous plants.

Generation of subspecies level-specific microbial diagnostic microarrays using genes amplified from subtractive suppression hybridization as microarray probes

The generation of microarray probes with specificity below the species level is an ongoing challenge, not least because the high-throughput detection of microorganisms would be an efficient means of identifying environmentally relevant microbes. Here, we describe how suppression subtractive hybridization (SSH) can be applied to the production of microarray probes that are useful for microbial differentiation at the subspecies level. SSH was used to initially isolate unique genomic sequences of nine Salmonella strains, and these were validated in quadruplicate by microarray analysis. The results obtained indicate that a large group of genes subtracted by SSH could serve together, as one probe, for detecting a microbial subspecies. Similarly, the whole microbial genome (not subjected to SSH) can be used as a species-specific probe. The detailed methods described herein could be used and adapted for the estimation of any cultivable bacteria from different environments.