Comparative analysis of early embryonic sunflower cDNA libraries

Molecular aspects of somatic-to-embryogenic transition in plants

Somatic embryogenesis (SE) is a model system for understanding the physiological, biochemical, and molecular biological events occurring during plant embryo development. Plant somatic cells have the ability to undergo sustained divisions and give rise to an entire organism. This remarkable feature is called plant cell totipotency. SE is a notable illustration of plant totipotency and involves reprogramming of development in somatic cells toward the embryogenic pathway. Plant growth regularities, especially auxins, are key components as their exogenous application recapitulates the embryogenic potential of the mitotically quiescent somatic cells. It has been observed that there are genetic and also physiological factors that trigger in vitro embryogenesis in various types of plant somatic cells. Analysis of the proteome and transcriptome has led to the identification and characterization of certain genes involved in SE. Most of these genes, however, are upregulated only in the late developmental stages, suggesting that they do not play a direct role in the vegetative-to-embryogenic transition. However, the molecular bases of those triggering factors and the genetic and biochemical mechanisms leading to in vitro embryogenesis are still unknown. Here, we describe the plant factors that participate in the vegetative-to-embryogenic transition and discuss their possible roles in this process.

SSRs and INDELs mined from the sunflower EST database: abundance, polymorphisms, and cross-taxa utility

Simple sequence repeats (SSRs) are abundant and frequently highly polymorphic in transcribed sequences and widely targeted for marker development in eukaryotes. Sunflower (Helianthus annuus) transcript assemblies were built and mined to identify SSRs and insertions-deletions (INDELs) for marker development, comparative mapping, and other genomics applications in sunflower. We describe the spectrum and frequency of SSRs identified in the sunflower EST database, a catalog of 16,643 EST-SSRs, a collection of 484 EST-SSR and 43 EST-INDEL markers developed from common sunflower ESTs, polymorphisms of the markers among the parents of several intraspecific and interspecific mapping populations, and the transferability of the markers to closely and distantly related species in the Compositae. Of 17,904 unigenes in the transcript assembly, 1,956 (10.9%) harbored one or more SSRs with repeat counts of n > or = 5. EST-SSR markers were 1.6-fold more polymorphic among exotic than elite genotypes and 0.7-fold less polymorphic than non-genic SSR markers. Of 466 EST-SSR or INDEL markers screened for cross-species amplification and polymorphisms, 413 (88.6%) amplified alleles from one or more wild species (H. argophyllus, H. tuberosus, H. anomalus, H. paradoxus, and H. deserticola), whereas 69 (14.8%) amplified alleles from safflower (Carthamus tinctorius) and 67 (14.4%) amplified alleles from lettuce (Lactuca sativa); hence, only a fraction were transferable to distantly related genera in the Compositae, whereas most were transferable to wild relatives of H. annuus. Several thousand additional SSRs were identified in the EST database and supply a wealth of templates for EST-SSR marker development in sunflower.

SSRs and INDELs mined from the sunflower EST database: abundance, polymorphisms, and cross-taxa utility

Simple sequence repeats (SSRs) are abundant and frequently highly polymorphic in transcribed sequences and widely targeted for marker development in eukaryotes. Sunflower (Helianthus annuus) transcript assemblies were built and mined to identify SSRs and insertions-deletions (INDELs) for marker development, comparative mapping, and other genomics applications in sunflower. We describe the spectrum and frequency of SSRs identified in the sunflower EST database, a catalog of 16,643 EST-SSRs, a collection of 484 EST-SSR and 43 EST-INDEL markers developed from common sunflower ESTs, polymorphisms of the markers among the parents of several intraspecific and interspecific mapping populations, and the transferability of the markers to closely and distantly related species in the Compositae. Of 17,904 unigenes in the transcript assembly, 1,956 (10.9%) harbored one or more SSRs with repeat counts of n > or = 5. EST-SSR markers were 1.6-fold more polymorphic among exotic than elite genotypes and 0.7-fold less polymorphic than non-genic SSR markers. Of 466 EST-SSR or INDEL markers screened for cross-species amplification and polymorphisms, 413 (88.6%) amplified alleles from one or more wild species (H. argophyllus, H. tuberosus, H. anomalus, H. paradoxus, and H. deserticola), whereas 69 (14.8%) amplified alleles from safflower (Carthamus tinctorius) and 67 (14.4%) amplified alleles from lettuce (Lactuca sativa); hence, only a fraction were transferable to distantly related genera in the Compositae, whereas most were transferable to wild relatives of H. annuus. Several thousand additional SSRs were identified in the EST database and supply a wealth of templates for EST-SSR marker development in sunflower.

A comprehensive analysis of the combined effects of high light and high temperature stresses on gene expression in sunflower

BACKGROUND AND AIMS

Although high light (HL) and high temperature (HT) stresses have been extensively investigated, a global analysis of their combined effects on the transcriptome of any plant species has not yet been described. Sunflower is an agronomically important oil crop frequently subjected to these stress factors. Because results in model plants may not always translate well to crop plants, responses of sunflower (Helianthus annuus) to HL, HT and a combination of both stresses were analysed by profiling gene expression in leaves and immature seeds.

METHODS

Plants were grown in HL (600 microE m(-2) s(-1)), HT (35 degrees C) and a combination of HL and HT (HL + HT), and gene expression in leaves and immature seeds was profiled using cDNA microarrays containing more than 8000 putative unigenes.

KEY RESULTS

Using two-way analysis of variance, 105, 55 and 129 cDNA clones were identified showing significant changes in steady-state transcript levels, across the two tissues, in response to HL, HT and HL + HT, respectively. A significant number of these transcripts were found to be specific to each stress. Comparing gene expression profiles between leaves and immature seeds revealed that 89, 113 and 186 cDNA clones can be considered as differentially expressed in response to HL, HT and HL + HT, respectively. More than half of the cDNA clones showing significant differences between embryo and leaf tissues in response to HL + HT were specific to this stress. Significant differences between leaves and seeds shared by all three stress treatments were observed for only eight genes.

CONCLUSIONS

Taken together, these results indicate that vegetative and reproductive tissues employ different transcriptome responses to these stress treatments. Careful examination of the putative functions of these genes revealed novel and specific responses. The potential roles of many of the differentially expressed genes in stress tolerance are mentioned and discussed.

Analysis of expression profile of selected genes expressed during auxin-induced somatic embryogenesis in leaf base system of wheat (Triticum aestivum) and their possible interactions

Somatic embryogenesis is a notable illustration of plant totipotency and involves reprogramming of development in somatic cells toward the embryogenic pathway. Auxins are key components as their exogenous application recuperates the embryogenic potential of the mitotically quiescent somatic cells. In order to unravel the molecular basis of somatic embryogenesis, cDNA library was made from the regeneration proficient wheat leaf base segments treated with auxin. In total, 1440 clones were sequenced and among these 1,196 good quality sequences were assembled into 270 contigs and 425 were singletons. By reverse northern analysis, a total of 57 clones were found to be upregulated during somatic embryogenesis, 64 during 2,4-D treatment, and 170 were common to 2,4-D treatment and somatic embryogenesis. A substantial number of genes involved in hormone response, signal transduction cascades, defense, anti-oxidation, programmed cell death/senescence and cell division were identified and characterized partially. Analysis of data of select genes suggests that the induction phase of somatic embryogenesis is accompanied by the expression of genes that may also be involved in zygotic embryogenesis. The developmental reprogramming process may in fact involve multiple cellular pathways and unfolding of as yet unknown molecular events. Thus, an interaction network draft using bioinformatics and system biology strategy was constructed. The outcome of a systematic and comprehensive analysis of somatic embryogenesis associated interactome in a monocot leaf base system is presented.

Carbon conversion efficiency and central metabolic fluxes in developing sunflower (Helianthus annuus L.) embryos

The efficiency with which developing sunflower embryos convert substrates into seed storage reserves was determined by labeling embryos with [U-(14)C6]glucose or [U-(14)C5]glutamine and measuring their conversion to CO2, oil, protein and other biomass compounds. The average carbon conversion efficiency was 50%, which contrasts with a value of over 80% previously observed in Brassica napus embryos (Goffman et al., 2005), in which light and the RuBisCO bypass pathway allow more efficient conversion of hexose to oil. Labeling levels after incubating sunflower embryos with [U-(14)C4]malate indicated that some carbon from malate enters the plastidic compartment and contributes to oil synthesis. To test this and to map the underlying pattern of metabolic fluxes, separate experiments were carried out in which embryos were labeled to isotopic steady state using [1-(13)C1]glucose, [2-(13)C1]glucose, or [U-(13)C5]glutamine. The resultant labeling in sugars, starch, fatty acids and amino acids was analyzed by NMR and GC-MS. The fluxes through intermediary metabolism were then quantified by computer-aided modeling. The resulting flux map accounted well for the labeling data, was in good agreement with the observed carbon efficiency, and was further validated by testing for agreement with gas exchange measurements. The map shows that the influx of malate into oil is low and that flux through futile cycles (wasting ATP) is low, which contrasts with the high rates previously determined for growing root tips and heterotrophic cell cultures.

Transcriptional profiles of primary metabolism and signal transduction-related genes in response to water stress in field-grown sunflower genotypes using a thematic cDNA microarray

A sunflower cDNA microarray containing about 800 clones covering major metabolic and signal transduction pathways was used to study gene expression profiles in leaves and embryos of drought-tolerant and -sensitive genotypes subjected to water-deficit stress under field conditions. Using two-step ANOVA normalization and analysis models, we identified 409 differentially expressed genes among genotypes, water treatment and organs. The majority of the cDNA clones differentially expressed under water stress was found to display opposite gene expression profiles in drought-tolerant genotype compared to drought-sensitive genotype. These dissimilarities suggest that the difference between tolerant and non-tolerant plants seems to be associated with changes in qualitative but not quantitative mRNA expression. Comparing leaves and embryos, 82 cDNA clones showing organ-specific variation in gene expression levels were identified in response to water stress across genotypes. Genes related to amino acids and carbohydrates metabolisms, and signal transduction were induced in embryos and repressed in leaves; suggesting that vegetative and reproductive organs respond differentially to water stress. Adaptive mechanisms controlling water deficit tolerance are proposed and discussed.

Transcript profiling and identification of molecular markers for early microspore embryogenesis in Brassica napus

Isolated microspores of Brassica napus are developmentally programmed to form gametes; however, microspores can be reprogrammed through stress treatments to undergo appropriate divisions and form embryos. We are interested in the identification and isolation of factors and genes associated with the induction and establishment of embryogenesis in isolated microspores. Standard and normalized cDNA libraries, as well as subtractive cDNA libraries, were constructed from freshly isolated microspores (0 h) and microspores cultured for 3, 5, or 7 d under embryogenesis-inducing conditions. Library comparison tools were used to identify shifts in metabolism across this time course. Detailed expressed sequence tag analyses of 3 and 5 d cultures indicate that most sequences are related to pollen-specific genes. However, semiquantitative and real-time reverse transcription-polymerase chain reaction analyses at the initial stages of embryo induction also reveal expression of embryogenesis-related genes such as BABYBOOM1, LEAFY COTYLEDON1 (LEC1), and LEC2 as early as 2 to 3 d of microspore culture. Sequencing results suggest that embryogenesis is clearly established in a subset of the microspores by 7 d of culture and that this time point is optimal for isolation of embryo-specific expressed sequence tags such as ABSCISIC ACID INSENSITIVE3, ATS1, LEC1, LEC2, and FUSCA3. Following extensive polymerase chain reaction-based expression profiling, 16 genes were identified as unequivocal molecular markers for microspore embryogenesis in B. napus. These molecular marker genes also show expression during zygotic embryogenesis, underscoring the common developmental pathways that function in zygotic and gametic embryogenesis. The quantitative expression values of several of these molecular marker genes are shown to be predictive of embryogenic potential in B. napus cultivars (e.g. 'Topas' DH4079, 'Allons,' 'Westar,' 'Garrison').

Heterosis associated gene expression in maize embryos 6 days after fertilization exhibits additive, dominant and overdominant pattern

Heterosis is important for conventional plant breeding and is intensively used to increase the productivity of crop plants. Genetic processes shortly after fertilization might be of particular importance with respect to heterosis, because coordination of the diverse genomes establishes a basis for future performance of the sporophyte. Here we demonstrate a strong crossbreeding advantage of hybrid maize embryos as early as 6 days after fertilization in a modern maize hybrid and provide the first embryo specific analysis of associated gene expression pattern at this early stage of development. We identified differentially expressed genes between hybrid embryos and the parental genotypes by a combined approach of suppression subtractive hybridization and differential screening by microarray hybridizations. Association of heterosis in embryos with genes related to signal transduction and other regulatory processes was implied by the enrichment of these functional classes among the identified gene set. Quantitative RT-PCR analysis validated the expression pattern of 7 of 12 genes analysed and revealed predominantly additive, but also dominant and overdominant expression patterns in hybrid embryos. These patterns indicate that gene regulatory interactions among parental alleles act at this early developmental stage and the genes identified provide entry points for the exploration of gene regulatory networks associated with the specification of the phenomenon heterosis in the plant life cycle.

Primary metabolic pathways and signal transduction in sunflower (Helianthus annuus L.): comparison of transcriptional profiling in leaves and immature embryos using cDNA microarrays

The early stage of embryo development is a critical step in plant production. To identify genes with potential roles in the early sunflower seed development, a cDNA microarray approach was employed. We developed a thematic cDNA microarray containing clones representing high sequence similarities with known or predicted Arabidopsis genes implicated in different metabolic and signal transduction pathways. This 800-element cDNA array was used to compare the expression patterns in leaves and immature embryos (2 mm and 6 mm). Statistical analysis, using two-step ANOVA, revealed that 143 cDNA clones can be considered as differentially expressed. Of these, 62 clones were found to be up-regulated in leaves, 81 in embryos whereas only seven clones displayed increased level of mRNA in the 6 mm embryos when compared with 2 mm embryos. The differentially expressed clones are distributed among many metabolic and signal transduction pathways. For example, genes related to fatty acid metabolism and amino acid biosynthesis exhibited preferential expression patterns in immature embryos. Also, clones potentially encoding enzymes involved in the metabolism of ascorbate and aldarate, pyruvate, propanoate and inositol, and citrate cycle were found to be up-regulated in embryos. In contrast, cDNA clones putatively involved in energy metabolism were more abundant in leaves than embryos. Clones encoding potential signal transduction components including receptors, protein kinases, protein phosphatases, and transcription factors were also identified, with preferential expression profiles in immature embryos. The expression patterns derived from this study provide initial characterization of metabolic pathways and signalling transduction networks occurring in the early stage of sunflower seed development.

A cDNA microarray approach to decipher sunflower (Helianthus annuus) responses to the necrotrophic fungus Phoma macdonaldii

To identify the genes involved in the partial resistance of sunflower (Helianthus annuus) to the necrotrophic fungus Phoma macdonaldii, we developed a 1000-element cDNA microarray containing carefully chosen genes putatively involved in primary metabolic pathways, signal transduction and biotic stress responses. A two-pass general linear model was used to normalize the data and then to detect differentially expressed genes. This method allowed us to identify 38 genes differentially expressed among genotypes, treatments and times, mainly belonging to plant defense, signaling pathways and amino acid metabolism. Based on a set of genes whose differential expression was highly significant, we propose a model in which negative regulation of a dual-specificity MAPK phosphatase could be implicated in sunflower defense mechanisms against the pathogen. The resulting activation of the MAP kinase cascade could subsequently trigger defense responses (e.g. thaumatin biosynthesis and phenylalanine ammonia lyase activation), under the control of transcription factors belonging to MYB and WRKY families. Concurrently, the activation of protein phosphatase 2A (PP2A), which is implicated in cell death inhibition, could limit pathogen development. The results reported here provide a valuable first step towards the understanding and analysis of the P. macdonaldii-sunflower interaction.

Transcriptional profiling of sunflower plants growing under low temperatures reveals an extensive down-regulation of gene expression associated with chilling sensitivity

Being able to sow early to maximize the growing season and to escape drought stress has increased the importance of low-temperature tolerance in sunflower. Yet knowledge about the molecular basis of sunflower response to low temperature is still lacking. To address this issue, nylon microarrays containing >8000 putative unigenes were developed and used. Early- and late-flowering genotypes were sown at 15 degrees C and grown until the two-leaf stage when they were subjected to 7 degrees C until the four-leaf stage. The transcriptional profiles of low temperature-grown plants (15 degrees C and 7 degrees C) were compared with those grown under standard conditions (25 degrees C). Two-step ANOVA normalization and analysis models were used to identify the differentially expressed genes. A total of 108 cDNA clones having a P-value <10(-3) were found to be differentially expressed between the low temperature-grown plants (15 degrees C and 7 degrees C) and their corresponding two-leaf- and four-leaf-stage controls across the two genotypes. About 90% of these genes were down-regulated. This includes genes potentially involved in the metabolism of carbohydrate and energy, protein synthesis, signal transduction, and transport function. Comparing gene expression profiles at 15 degrees C and 7 degrees C revealed that only four genes can be considered as differentially expressed, in both genotypes, suggesting that similar genetic programmes underlie the response of sunflower plants to these temperature regimes. The analysis also revealed that early- and late-flowering genotypes respond similarly to low-temperature tolerance as justified by the low number of genes showing a significant genotype x treatment interaction effect. It seems likely that the down-regulation and/or non-induction of genes having a critical role in low-temperature tolerance may be responsible for the sensitivity of sunflower plants to low-temperature tolerance. The results reported provide an initial characterization of the transcriptome activity of sunflower, as a chilling-sensitive plant under suboptimal temperatures, and could be of importance to reveal the potential differences between chilling-sensitive and chilling-tolerant species.

A 6374 unigene set corresponding to low abundance transcripts expressed following fertilization in Solanum chacoense Bitt, and characterization of 30 receptor-like kinases

In order to characterize regulatory genes that are expressed in ovule tissues after fertilization we have undertaken an EST sequencing project in Solanum chacoense, a self-incompatible wild potato species. Two cDNA libraries made from ovule tissues covering embryo development from zygote to late torpedo-stage were constructed and plated at high density on nylon membranes. To decrease EST redundancy and enrich for transcripts corresponding to weakly expressed genes a self-probe subtraction method was used to select the colonies harboring the genes to be sequenced. 7741 good sequences were obtained and, from these, 6374 unigenes were isolated. Thus, the self-probe subtraction resulted in a strong enrichment in singletons, a decrease in the number of clones per contigs, and concomitantly, an enrichment in the total number of unigenes obtained (82%). To gain insights into signal transduction events occurring during embryo development all the receptor-like kinases (or protein receptor kinases) were analyzed by quantitative real-time RT-PCR. Interestingly, 28 out of the 30 RLK isolated were predominantly expressed in ovary tissues or young developing fruits, and 23 were transcriptionaly induced following fertilization. Thus, the self-probe subtraction did not select for genes weakly expressed in the target tissue while being highly expressed elsewhere in the plant. Of the receptor-like kinases (RLK) genes isolated, the leucine-rich repeat (LRR) family of RLK was by far the most represented with 25 members covering 11 LRR classes.