The gene expression landscape of pine seedling tissues

Properties and Functional Analysis of Two Chorismate Mutases from Maritime Pine

Through the shikimate pathway, a massive metabolic flux connects the central carbon metabolism with the synthesis of chorismate, the common precursor of the aromatic amino acids phenylalanine, tyrosine, and tryptophan, as well as other compounds, including salicylate or folate. The alternative metabolic channeling of chorismate involves a key branch-point, finely regulated by aromatic amino acid levels. Chorismate mutase catalyzes the conversion of chorismate to prephenate, a precursor of phenylalanine and tyrosine and thus a vast repertoire of fundamental derived compounds, such as flavonoids or lignin. The regulation of this enzyme has been addressed in several plant species, but no study has included conifers or other gymnosperms, despite the importance of the phenolic metabolism for these plants in processes such as lignification and wood formation. Here, we show that maritime pine (Pinus pinaster Aiton) has two genes that encode for chorismate mutase, PpCM1 and PpCM2. Our investigations reveal that these genes encode plastidial isoenzymes displaying activities enhanced by tryptophan and repressed by phenylalanine and tyrosine. Using phylogenetic studies, we have provided new insights into the possible evolutionary origin of the cytosolic chorismate mutases in angiosperms involved in the synthesis of phenylalanine outside the plastid. Studies based on different platforms of gene expression and co-expression analysis have allowed us to propose that PpCM2 plays a central role in the phenylalanine synthesis pathway associated with lignification.

Under Stress: Searching for Genes Involved in the Response of Abies pinsapo Boiss to Climate Change

Currently, Mediterranean forests are experiencing the deleterious effects of global warming, which mainly include increased temperatures and decreased precipitation in the region. Relict Abies pinsapo fir forests, endemic in the southern Iberian Peninsula, are especially sensitive to these recent environmental disturbances, and identifying the genes involved in the response of this endangered tree species to climate-driven stresses is of paramount importance for mitigating their effects. Genomic resources for A. pinsapo allow for the analysis of candidate genes reacting to warming and aridity in their natural habitats. Several members of the complex gene families encoding late embryogenesis abundant proteins (LEAs) and heat shock proteins (HSPs) have been found to exhibit differential expression patterns between wet and dry seasons when samples from distinct geographical locations and dissimilar exposures to the effects of climate change were analyzed. The observed changes were more perceptible in the roots of trees, particularly in declining forests distributed at lower altitudes in the more vulnerable mountains. These findings align with previous studies and lay the groundwork for further research on the molecular level. Molecular and genomic approaches offer valuable insights for mitigating climate stress and safeguarding this endangered conifer.

Pine has two glutamine synthetase paralogs, GS1b.1 and GS1b.2, exhibiting distinct biochemical properties

The enzyme glutamine synthetase (EC 6.3.1.2) is mainly responsible for the incorporation of inorganic nitrogen into organic molecules in plants. In the present work, a pine (Pinus pinaster) GS1 (PpGS1b.2) gene was identified, showing a high sequence identity with the GS1b.1 gene previously characterized in conifers. Phylogenetic analysis revealed that the presence of PpGS1b.2 is restricted to the genera Pinus and Picea and is not found in other conifers. Gene expression data suggest a putative role of PpGS1b.2 in plant development, similar to other GS1b genes from angiosperms, suggesting evolutionary convergence. The characterization of GS1b.1 and GS1b.2 at the structural, physicochemical, and kinetic levels has shown differences even though they have high sequence homology. GS1b.2 had a lower optimum pH (6 vs. 6.5) and was less thermally stable than GS1b.1. GS1b.2 exhibited positive cooperativity for glutamate and substrate inhibition for ammonium. However, GS1b.1 exhibited substrate inhibition behavior for glutamate and ATP. Alterations in the kinetic characteristics produced by site-directed mutagenesis carried out in this work strongly suggest an implication of amino acids at positions 264 and 267 in the active center of pine GS1b.1 and GS1b.2 being involved in affinity toward ammonium. Therefore, the amino acid differences between GS1b.1 and GS1b.2 would support the functioning of both enzymes to meet distinct plant needs.

SNP Detection in Pinus pinaster Transcriptome and Association with Resistance to Pinewood Nematode

Pinewood nematode (PWN, Bursaphelenchus xylophilus) is the causal agent of pine wilt disease (PWD), which severely affects Pinus pinaster stands in southwestern Europe. Despite the high susceptibility of P. pinaster, individuals of selected half-sib families have shown genetic variability in survival after PWN inoculation, indicating that breeding for resistance can be a valuable strategy to control PWD. In this work, RNA-seq data from susceptible and resistant plants inoculated with PWN were used for SNP discovery and analysis. A total of 186,506 SNPs were identified, of which 31 were highly differentiated between resistant and susceptible plants, including SNPs in genes involved in cell wall lignification, a process previously linked to PWN resistance. Fifteen of these SNPs were selected for validation through Sanger sequencing and 14 were validated. To evaluate SNP-phenotype associations, 40 half-sib plants were genotyped for six validated SNPs. Associations with phenotype after PWN inoculation were found for two SNPs in two different genes (MEE12 and PCMP-E91), as well as two haplotypes of HIPP41, although significance was not maintained following Bonferroni correction. SNPs here detected may be useful for the development of molecular markers for PWD resistance and should be further investigated in future association studies.

Identification of Metabolic Pathways Differentially Regulated in Somatic and Zygotic Embryos of Maritime Pine

Embryogenesis is a complex phase of conifer development involving hundreds of genes, and a proper understanding of this process is critical not only to produce embryos with different applied purposes but also for comparative studies with angiosperms. A global view of transcriptome dynamics during pine somatic and zygotic embryogenesis is currently missing. Here, we present a genome-wide transcriptome analysis of somatic and zygotic embryos at three developmental stages to identify conserved biological processes and gene functions during late embryogenesis. Most of the differences became more significant as the developmental process progressed from early to cotyledonary stages, and a higher number of genes were differentially expressed in somatic than in zygotic embryos. Metabolic pathways substantially affected included those involved in amino acid biosynthesis and utilization, and this difference was already observable at early developmental stages. Overall, this effect was found to be independent of the line (genotype) used to produce the somatic embryos. Additionally, transcription factors differentially expressed in somatic versus zygotic embryos were analyzed. Some potential hub regulatory genes were identified that can provide clues as to what transcription factors are controlling the process and to how the observed differences between somatic and zygotic embryogenesis in conifers could be regulated.

miR160 Interacts in vivo With Pinus pinaster AUXIN RESPONSE FACTOR 18 Target Site and Negatively Regulates Its Expression During Conifer Somatic Embryo Development

MicroRNAs (miRNAs) are key regulators of several plant developmental processes including embryogenesis. Most miRNA families are conserved across major groups of plant species, but their regulatory roles have been studied mainly in model species like Arabidopsis and other angiosperms. In gymnosperms, miRNA-dependent regulation has been less studied since functional approaches in these species are often difficult to establish. Given the fundamental roles of auxin signaling in somatic embryogenesis (SE) induction and embryo development, we investigated a previously predicted interaction between miR160 and a putative target encoding AUXIN RESPONSE FACTOR 18 in Pinus pinaster (PpARF18) embryonic tissues. Phylogenetic analysis of AUXIN RESPONSE FACTOR 18 (ARF18) from Pinus pinaster and Picea abies, used here as a model system of conifer embryogenesis, showed their close relatedness to AUXIN RESPONSE FACTOR (ARF) genes known to be targeted by miR160 in other species, including Arabidopsis ARF10 and ARF16. By using a luciferase (LUC) reporter system for miRNA activity in Arabidopsis protoplasts, we have confirmed that P. pinaster miR160 (ppi-miR160) interacts in vivo with PpARF18 target site. When the primary miR160 from P. pinaster was overexpressed in protoplasts under non-limiting levels of ARGONAUTE1, a significant increase of miR160 target cleavage activity was observed. In contrast, co-expression of the primary miRNA and the target mimic MIM160 led to a decrease of miR160 activity. Our results further support that this interaction is functional during consecutive stages of SE in the conifer model P. abies. Expression analyses conducted in five stages of development, from proembryogenic masses (PEMs) to the mature embryo, show that conifer ARF18 is negatively regulated by miR160 toward the fully developed mature embryo when miR160 reached its highest expression level. This study reports the first in vivo validation of a predicted target site of a conifer miRNA supporting the conservation of miR160 interaction with ARF targets in gymnosperms. The approach used here should be useful for future characterization of miRNA functions in conifer embryogenesis.

Ammonium regulates the development of pine roots through hormonal crosstalk and differential expression of transcription factors in the apex

Ammonium is a prominent source of inorganic nitrogen for plant nutrition, but excessive amounts can be toxic for many species. However, most conifers are tolerant to ammonium, a relevant physiological feature of this ancient evolutionary lineage. For a better understanding of the molecular basis of this trait, ammonium-induced changes in the transcriptome of maritime pine (Pinus pinaster Ait.) root apex have been determined by laser capture microdissection and RNA sequencing. Ammonium promoted changes in the transcriptional profiles of multiple transcription factors, such as SHORT-ROOT, and phytohormone-related transcripts, such as ACO, involved in the development of the root meristem. Nano-PALDI-MSI and transcriptomic analyses showed that the distributions of IAA and CKs were altered in the root apex in response to ammonium nutrition. Taken together, the data suggest that this early response is involved in the increased lateral root branching and principal root growth, which characterize the long-term response to ammonium supply in pine. All these results suggest that ammonium induces changes in the root system architecture through the IAA-CK-ET phytohormone crosstalk and transcriptional regulation.

The amino acid permease PpAAP1 mediates arginine transport in maritime pine

Forest trees have access to diverse nitrogenous compounds in the soil such as ammonium, nitrate and amino acids. Recent progress has been made in the identification and characterization of ammonium and nitrate transporters. However, much more limited is our current knowledge of amino acid transport systems despite their relevance to fully understanding nitrogen nutrition in trees. In the present study, we have identified 10 genes encoding putative amino acid permeases of the AAP family in maritime pine (Pinus pinaster Ait.). Four members of this family, PpAAP1, PpAAP2, PpAAP3 and PpAAP4 were phylogenetically related to AtAAP5, involved in arginine transport in Arabidopsis thaliana. One of these genes, PpAAP1, exhibited enhanced expression levels in maritime pine roots when arginine was externally supplied. PpAAP1 was functionally characterized by complementation of a yeast mutant strain defective in the transport of arginine, allowing yeast to take up [14C]-arginine with high affinity. Furthermore, PpAAP1 was able to restore the severely affected root uptake of arginine displayed by AtAAP5 T-DNA mutants. Uptake rates of 15N-labelled arginine were significantly higher in PpAAP1-overexpressing plants when compared to wild-type and AtAAP5 mutant plants. Taken together, our results indicate that PpAAP1 is a high-affinity arginine transporter in maritime pine.

Atlas of tissue-specific and tissue-preferential gene expression in ecologically and economically significant conifer Pinus sylvestris

Despite their ecological and economical importance, conifers genomic resources are limited, mainly due to the large size and complexity of their genomes. Additionally, the available genomic resources lack complete structural and functional annotation. Transcriptomic resources have been commonly used to compensate for these deficiencies, though for most conifer species they are limited to a small number of tissues, or capture only a fraction of the genes present in the genome. Here we provide an atlas of gene expression patterns for conifer Pinus sylvestris across five tissues: embryo, megagametophyte, needle, phloem and vegetative bud. We used a wide range of tissues and focused our analyses on the expression profiles of genes at tissue level. We provide comprehensive information of the per-tissue normalized expression level, indication of tissue preferential upregulation and tissue-specificity of expression. We identified a total of 48,001 tissue preferentially upregulated and tissue specifically expressed genes, of which 28% have annotation in the Swiss-Prot database. Even though most of the putative genes identified do not have functional information in current biological databases, the tissue-specific patterns discovered provide valuable information about their potential functions for further studies, as for example in the areas of plant physiology, population genetics and genomics in general. As we provide information on tissue specificity at both diploid and haploid life stages, our data will also contribute to the understanding of evolutionary rates of different tissue types and ploidy levels.

Chia (Salvia hispanica) Gene Expression Atlas Elucidates Dynamic Spatio-Temporal Changes Associated With Plant Growth and Development

Chia (Salvia hispanica L.), now a popular superfood and a pseudocereal, is one of the richest sources of dietary nutrients such as protein, fiber, and polyunsaturated fatty acids (PUFAs). At present, the genomic and genetic information available in the public domain for this crop are scanty, which hinders an understanding of its growth and development and genetic improvement. We report an RNA-sequencing (RNA-Seq)-based comprehensive transcriptome atlas of Chia sampled from 13 tissue types covering vegetative and reproductive growth stages. We used ~355 million high-quality reads of total ~394 million raw reads from transcriptome sequencing to generate de novo reference transcriptome assembly and the tissue-specific transcript assemblies. After the quality assessment of the merged assemblies and implementing redundancy reduction methods, 82,663 reference transcripts were identified. About 65,587 of 82,663 transcripts were translated into 99,307 peptides, and we were successful in assigning InterPro annotations to 45,209 peptides and gene ontology (GO) terms to 32,638 peptides. The assembled transcriptome is estimated to have the complete sequence information for ~86% of the genes found in the Chia genome. Furthermore, the analysis of 53,200 differentially expressed transcripts (DETs) revealed their distinct expression patterns in Chia's vegetative and reproductive tissues; tissue-specific networks and developmental stage-specific networks of transcription factors (TFs); and the regulation of the expression of enzyme-coding genes associated with important metabolic pathways. In addition, we identified 2,411 simple sequence repeats (SSRs) as potential genetic markers from the transcripts. Overall, this study provides a comprehensive transcriptome atlas, and SSRs, contributing to building essential genomic resources to support basic research, genome annotation, functional genomics, and molecular breeding of Chia.

Insights Into the Mechanisms Implicated in Pinus pinaster Resistance to Pinewood Nematode

Pine wilt disease (PWD), caused by the plant-parasitic nematode Bursaphelenchus xylophilus, has become a severe environmental problem in the Iberian Peninsula with devastating effects in Pinus pinaster forests. Despite the high levels of this species' susceptibility, previous studies reported heritable resistance in P. pinaster trees. Understanding the basis of this resistance can be of extreme relevance for future programs aiming at reducing the disease impact on P. pinaster forests. In this study, we highlighted the mechanisms possibly involved in P. pinaster resistance to PWD, by comparing the transcriptional changes between resistant and susceptible plants after infection. Our analysis revealed a higher number of differentially expressed genes (DEGs) in resistant plants (1,916) when compared with susceptible plants (1,226). Resistance to PWN is mediated by the induction of the jasmonic acid (JA) defense pathway, secondary metabolism pathways, lignin synthesis, oxidative stress response genes, and resistance genes. Quantification of the acetyl bromide-soluble lignin confirmed a significant increase of cell wall lignification of stem tissues around the inoculation zone in resistant plants. In addition to less lignified cell walls, susceptibility to the pine wood nematode seems associated with the activation of the salicylic acid (SA) defense pathway at 72 hpi, as revealed by the higher SA levels in the tissues of susceptible plants. Cell wall reinforcement and hormone signaling mechanisms seem therefore essential for a resistance response.

Construction of a reference transcriptome for the analysis of male sterility in sugi (Cryptomeria japonica D. Don) focusing on MALE STERILITY 1 (MS1)

Sugi (Cryptomeria japonica D. Don) is an important conifer used for afforestation in Japan. As the genome of this species is 11 Gbps, it is too large to assemble within a short timeframe. Transcriptomics is one approach that can address this deficiency. Here we designed a workflow consisting of three stages to de novo assemble transcriptome using Oases and Trinity. The three transcriptomic stage used were independent assembly, automatic and semi-manual integration, and refinement by filtering out potential contamination. We identified a set of 49,795 cDNA and an equal number of translated proteins. According to the benchmark set by BUSCO, 87.01% of cDNAs identified were complete genes, and 78.47% were complete and single-copy genes. Compared to other full-length cDNA resources collected by Sanger and PacBio sequencers, the extent of the coverage in our dataset was the highest, indicating that these data can be safely used for further studies. When two tissue-specific libraries were compared, there were significant expression differences between male strobili and leaf and bark sets. Moreover, subtle expression difference between male-fertile and sterile libraries were detected. Orthologous genes from other model plants and conifer species were identified. We demonstrated that our transcriptome assembly output (CJ3006NRE) can serve as a reference transcriptome for future functional genomics and evolutionary biology studies.

Systematic investigation and expression profiles of the GbR2R3-MYB transcription factor family in ginkgo (Ginkgo biloba L.)

As one of the largest families of transcription factors, the R2R3-MYB family plays a significant role in plant growth, development, and response to hormone and environmental stress. To explore its evolutionary mechanism and potential function in Ginkgo biloba, a gymnosperm of great economic and ecological value, we presented a comprehensive analysis of the R2R3-MYB genes in ginkgo. Sixty-nine GbR2R3-MYB genes were identified and these genes could be classified into 33 groups based on the characteristics of the amino acid sequence of the R2R3-MYB domain and gene structure. Syntenic analyses indicated that few tandem and segmental duplications possibly resulted in the contraction of the GbR2R3-MYB gene family. Based on the transcriptome data, expression profiles of eight different tissues and different developmental stages of leaf and kernel showed that GbR2R3-MYB genes had distinct temporal and spatial expression characteristics. Specific expression patterns of the sixteen GbR2R3-MYB genes were also identified in response to different abiotic stresses and hormonal exposures. Further investigation revealed that GbR2R3-MYB19 was located in the nucleus and possessed transcriptional activity, implying its potential roles in the regulation of multiple biological processes. Our findings provide a robust basis for future comprehensive evolutionary and functional analyses of GbR2R3-MYB genes in ginkgo.

Characterization of the KNOTTED1-LIKE HOMEOBOX (KNOX) gene family in Pinus pinaster Ait

KNOTTED1-LIKE HOMEOBOX (KNOX) genes are a family of plant-specific homeobox transcription factors with important roles in plant development that have been classified into two subfamilies with differential expression domains and functions. Studies in angiosperms have shown that class I members are related to the maintenance of meristem homeostasis and leaf development, whereas class II members promote differentiation of tissues and organs. However, little is known about its diversification and function in gymnosperms. By combining PCR-based detection and transcriptome data analysis, we identified four class I and two class II KNOX genes in Pinus pinaster. Expression analyses showed that class I members were mainly expressed in meristematic regions and differentiating tissues, with practically no expression in lateral organs, whereas expression of class II members was restricted to lateral organs. Furthermore, overexpression of P. pinaster KNOX genes in Arabidopsis thaliana caused similar phenotypic effects to those described for their angiosperms counterparts. This is the first time to our knowledge that functional analyses of class II members are reported in a conifer species. These results suggest a high conservation of the KNOX gene family throughout seed plants, as the functional differentiation of both subfamilies observed in angiosperms might be partially conserved in gymnosperms.

Enzymes Involved in the Biosynthesis of Arginine from Ornithine in Maritime Pine (Pinus pinaster Ait.)

The amino acids arginine and ornithine are the precursors of a wide range of nitrogenous compounds in all living organisms. The metabolic conversion of ornithine into arginine is catalyzed by the sequential activities of the enzymes ornithine transcarbamylase (OTC), argininosuccinate synthetase (ASSY) and argininosuccinate lyase (ASL). Because of their roles in the urea cycle, these enzymes have been purified and extensively studied in a variety of animal models. However, the available information about their molecular characteristics, kinetic and regulatory properties is relatively limited in plants. In conifers, arginine plays a crucial role as a main constituent of N-rich storage proteins in seeds and serves as the main source of nitrogen for the germinating embryo. In this work, recombinant PpOTC, PpASSY and PpASL enzymes from maritime pine (Pinus pinaster Ait.) were produced in Escherichia coli to enable study of their molecular and kinetics properties. The results reported here provide a molecular basis for the regulation of arginine and ornithine metabolism at the enzymatic level, suggesting that the reaction catalyzed by OTC is a regulatory target in the homeostasis of ornithine pools that can be either used for the biosynthesis of arginine in plastids or other nitrogenous compounds in the cytosol.

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Structural and Functional Characteristics of Two Molecular Variants of the Nitrogen Sensor PII in Maritime Pine

High levels of nitrogen are stored as arginine during the last stages of seed formation in maritime pine (Pinus pinaster Aiton). The protein sensor PII regulates the feedback inhibition of arginine biosynthesis through interaction with the key enzyme N-acetylglutamate kinase (NAGK). In this study, the structural and functional characteristics of PII have been investigated in maritime pine to get insights into the regulation of arginine metabolism. Two different forms of PII have been identified, PpPIIa and PpPIIb, which differ in their amino acid sequence and most likely correspond to splicing variants of a single gene in the pine genome. Two PII variants are also present in other pine species but not in other conifers such as spruces. PpPIIa and PpPIIb are trimeric proteins for which structural modeling predicts similar tridimensional protein core structures. Both are located in the chloroplast, where the PII-target enzyme PpNAGK is also found. PpPIIa, PpPIIb, and PpNAGK have been recombinantly produced to investigate the formation of NAGK-PII complexes. The interaction of PpPIIa/PpPIIb and PpNAGK may be enhanced by glutamine and contribute to relieve the feedback inhibition of PpNAGK by arginine. Expression analysis of PpPII genes revealed that PpIIa transcripts were predominant during embryogenesis and germination. The potential roles of PpPIIa and PpPIIb in the regulation of arginine metabolism of maritime pine are discussed.

Gene Expression Maps in Plants: Current State and Prospects

For many years, progress in the identification of gene functions has been based on classical genetic approaches. However, considerable recent omics developments have brought to the fore indirect but high-resolution methods of gene function identification such as transcriptomics, proteomics, and metabolomics. A transcriptome map is a powerful source of functional information and the result of the genome-wide expression analysis of a broad sampling of tissues and/or organs from different developmental stages and/or environmental conditions. In plant science, the application of transcriptome maps extends from the inference of gene regulatory networks to evolutionary studies. However, only some of these data have been integrated into databases, thus enabling analyses to be conducted without raw data; without this integration, extensive data preprocessing is required, which limits data usability. In this review, we summarize the state of plant transcriptome maps, analyze the problems associated with the combined analysis of large-scale data from various studies, and outline possible solutions to these problems.

An update to database TraVA: organ-specific cold stress response in Arabidopsis thaliana

BACKGROUND

Transcriptome map is a powerful tool for a variety of biological studies; transcriptome maps that include different organs, tissues, cells and stages of development are currently available for at least 30 plants. Some of them include samples treated by environmental or biotic stresses. However, most studies explore only limited set of organs and developmental stages (leaves or seedlings). In order to provide broader view of organ-specific strategies of cold stress response we studied expression changes that follow exposure to cold (+ 4 °C) in different aerial parts of plant: cotyledons, hypocotyl, leaves, young flowers, mature flowers and seeds using RNA-seq.

RESULTS

The results on differential expression in leaves are congruent with current knowledge on stress response pathways, in particular, the role of CBF genes. In other organs, both essence and dynamics of gene expression changes are different. We show the involvement of genes that are confined to narrow expression patterns in non-stress conditions into stress response. In particular, the genes that control cell wall modification in pollen, are activated in leaves. In seeds, predominant pattern is the change of lipid metabolism.

CONCLUSIONS

Stress response is highly organ-specific; different pathways are involved in this process in each type of organs. The results were integrated with previously published transcriptome map of Arabidopsis thaliana and used for an update of a public database TraVa: http://travadb.org/browse/Species=AthStress .

PpNAC1, a main regulator of phenylalanine biosynthesis and utilization in maritime pine

The transcriptional regulation of phenylalanine metabolism is particularly important in conifers, long-lived species that use large amounts of carbon in wood. Here, we show that the Pinus pinaster transcription factor, PpNAC1, is a main regulator of phenylalanine biosynthesis and utilization. A phylogenetic analysis classified PpNAC1 in the NST proteins group and was selected for functional characterization. PpNAC1 is predominantly expressed in the secondary xylem and compression wood of adult trees. Silencing of PpNAC1 in P. pinaster results in the alteration of stem vascular radial patterning and the down-regulation of several genes associated with cell wall biogenesis and secondary metabolism. Furthermore, transactivation and EMSA analyses showed that PpNAC1 is able to activate its own expression and PpMyb4 promoter, while PpMyb4 is able to activate PpMyb8, a transcriptional regulator of phenylalanine and lignin biosynthesis in maritime pine. Together, these results suggest that PpNAC1 is a functional ortholog of the ArabidopsisSND1 and NST1 genes and support the idea that key regulators governing secondary cell wall formation could be conserved between gymnosperms and angiosperms. Understanding the molecular switches controlling wood formation is of paramount importance for fundamental tree biology and paves the way for applications in conifer biotechnology.

Analysis of the WUSCHEL-RELATED HOMEOBOX gene family in Pinus pinaster: New insights into the gene family evolution

WUSCHEL-RELATED HOMEOBOX (WOX) genes are key players controlling stem cells in plants and can be divided into three clades according to the time of their appearance during plant evolution. Our knowledge of stem cell function in vascular plants other than angiosperms is limited, they separated from gymnosperms ca 300 million years ago and their patterning during embryogenesis differs significantly. For this reason, we have used the model gymnosperm Pinus pinaster to identify WOX genes and perform a thorough analysis of their gene expression patterns. Using transcriptomic data from a comprehensive range of tissues and stages of development we have shown three major outcomes: that the P. pinaster genome encodes at least fourteen members of the WOX family spanning all the major clades, that the genome of gymnosperms contains a WOX gene with no homologues in angiosperms representing a transitional stage between intermediate- and WUS-clade proteins, and that we can detect discrete WUS and WOX5 transcripts for the first time in a gymnosperm.