De novo assembly of maritime pine transcriptome: implications for forest breeding and biotechnology

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.

Conifer Biotechnology: An Overview

The peculiar characteristics of conifers determine the difficulty of their study and their great importance from various points of view. However, their study faces numerous important scientific, methodological, cultural, economic, social, and legal challenges. This paper presents an approach to several of those challenges and proposes a multidisciplinary scientific perspective that leads to a holistic understanding of conifers from the perspective of the latest technical, computer, and scientific advances. This review highlights the deep connection that all scientific contributions to conifers can have in each other as fully interrelated communicating vessels.

Comparative transcriptome analysis reveals the resistance regulation mechanism and fungicidal activity of the fungicide phenamacril in Fusarium oxysporum

Fusarium oxysporum (Fo) is an important species complex of soil-borne pathogenic fungi that cause vascular wilt diseases of agricultural crops and some opportunistic diseases of humans. The fungicide phenamacril has been extensively reported to have antifungal activity against Fusarium graminearum and Fusarium fujikuroi. In this study, we found that the amino acid substitutions (V151A and S418T) in Type I myosin FoMyo5 cause natural low resistance to phenamacril in the plant pathogenic Fo isolates. Therefore, we compared the transcriptomes of two phenamacril-resistant Fo isolates FoII5, Fo1st and one phenamacril-sensitive isolate Fo3_a after 1 μg/mL phenamacril treatment. Among the 2728 differentially expressed genes (DEGs), 14 DEGs involved in oxidation–reduction processes and MFS transporters, were significantly up-regulated in phenamacril-resistant isolates. On the other hand, 14 DEGs involved in ATP-dependent RNA helicase and ribosomal biogenesis related proteins, showed significantly down-regulated expression in both phenamacril-resistant and -sensitive isolates. These results indicated that phenamacril not only seriously affected the cytoskeletal protein binding and ATPase activity of sensitive isolate, but also suppressed ribosome biogenesis in all the isolates. Hence, this study helps us better understand resistance regulation mechanism and fungicidal activity of phenamacril and provide reference for the development of new fungicides to control Fo.

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.

Transcriptional Regulation of Pine Male and Female Cone Initiation and Development: Key Players Identified Through Comparative Transcriptomics

With long reproductive timescales, large complex genomes, and a lack of reliable reference genomes, understanding gene function in conifers is extremely challenging. Consequently, our understanding of which genetic factors influence the development of reproductive structures (cones) in monoecious conifers remains limited. Genes with inferred roles in conifer reproduction have mostly been identified through homology and phylogenetic reconstruction with their angiosperm counterparts. We used RNA-sequencing to generate transcriptomes of the early morphological stages of cone development in the conifer species Pinus densiflora and used these to gain a deeper insight into the transcriptional changes during male and female cone development. Paired-end Illumina sequencing was used to generate transcriptomes from non-reproductive tissue and male and female cones at four time points with a total of 382.82 Gbp of data generated. After assembly and stringent filtering, a total of 37,164 transcripts were retrieved, of which a third were functionally annotated using the Mercator plant pipeline. Differentially expressed gene (DEG) analysis resulted in the identification of 172,092 DEGs in the nine tissue types. This, alongside GO gene enrichment analyses, pinpointed transcripts putatively involved in conifer reproductive structure development, including co-orthologs of several angiosperm flowering genes and several that have not been previously reported in conifers. This study provides a comprehensive transcriptome resource for male and early female cone development in the gymnosperm species Pinus densiflora. Characterisation of this resource has allowed the identification of potential key players and thus provides valuable insights into the molecular regulation of reproductive structure development in monoecious conifers.

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.

Constitutive Overexpression of a Conifer WOX2 Homolog Affects Somatic Embryo Development in Pinus pinaster and Promotes Somatic Embryogenesis and Organogenesis in Arabidopsis Seedlings

Although full sequence data of several embryogenesis-related genes are available in conifers, their functions are still poorly understood. In this study, we focused on the transcription factor WUSCHEL-related HOMEOBOX 2 (WOX2), which is involved in determination of the apical domain during early embryogenesis, and is required for initiation of the stem cell program in the embryogenic shoot meristem of Arabidopsis. We studied the effects of constitutive overexpression of Pinus pinaster WOX2 (PpWOX2) by Agrobacterium-mediated transformation of P. pinaster somatic embryos and Arabidopsis seedlings. Overexpression of PpWOX2 during proliferation and maturation of somatic embryos of P. pinaster led to alterations in the quantity and quality of cotyledonary embryos. In addition, transgenic somatic seedlings of P. pinaster showed non-embryogenic callus formation in the region of roots and subsequently inhibited root growth. Overexpression of PpWOX2 in Arabidopsis promoted somatic embryogenesis and organogenesis in a part of the transgenic seedlings of the first and second generations. A concomitant increased expression of endogenous embryogenesis-related genes such as AtLEC1 was detected in transgenic plants of the first generation. Various plant phenotypes observed from single overexpressing transgenic lines of the second generation suggest some significant interactions between PpWOX2 and AtWOX2. As an explanation, functional redundancy in the WOX family is suggested for seed plants. Our results demonstrate that the constitutive high expression of PpWOX2 in Arabidopsis and P. pinaster affected embryogenesis-related traits. These findings further support some evolutionary conserved roles of this gene in embryo development of seed plants and have practical implications toward somatic embryogenesis induction in conifers.

Polygenic adaptation and negative selection across traits, years and environments in a long-lived plant species (Pinus pinaster Ait., Pinaceae)

A decade of genetic association studies in multiple organisms suggests that most complex traits are polygenic, i.e., they have a genetic architecture determined by numerous loci each with small effect-size. Thus, determining the degree of polygenicity and its variation across traits, environments and time is crucial to understand the genetic basis of phenotypic variation. We applied multilocus approaches to estimate the degree of polygenicity of fitness-related traits in a long-lived plant (Pinus pinaster Ait., maritime pine) and to analyze this variation across environments and years. We evaluated five categories of fitness-related traits (survival, height, phenology, functional, and biotic-stress response traits) in a clonal common-garden network, planted in contrasted environments (over 12,500 trees). Most of the analyzed traits showed evidence of local adaptation based on Q_st -F_st comparisons. We further observed a remarkably stable degree of polygenicity, averaging 6% (range of 0-27%), across traits, environments and years. We detected evidence of negative selection, which could explain, at least partially, the high degree of polygenicity. Because polygenic adaptation can occur rapidly, our results suggest that current predictions on the capacity of natural forest tree populations to adapt to new environments should be revised, especially in the current context of climate change.

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.

Transcriptome analysis identifies CsNRT genes involved in nitrogen uptake in tea plants, with a major role of CsNRT2.4

Tea trees have a high demand for nitrogen (N) fertilizer to improve the yield and quality of tea. In this research, transcriptome analysis revealed the effect of N starvation and resupply upon N uptake in tea plants. We identified 4098 differentially expressed genes (DEGs) that were significantly enriched in amino acid and N metabolism and were extensively mapped to the tea genome. The CsNRT gene family plays vital roles in the nitrogen uptake of tea plants. The full CDS sequences of CsNRT1.1, CsNRT1.2, CsNRT1.5, CsNRT1.7, CsNRT2.4, CsNRT2.5, CsNRT3.1 and CsNRT3.2 were cloned. One-year-old cutting seedlings of Zhongcha302 (ZC302) were selected for hydroponic culture and were used for gene expression analysis. The seedlings were resupplied with 0.2 and 2 mM N after N starvation. The results of the gene expression under different N treatments and in various tissues indicated that the expression of CsNRT2.4 was highly expressed in tea roots and was greatly induced by N. Overexpressed CsNRT2.4 in transgenic Arabidopsis thaliana increased the root lengths and fresh weights and improved the NO₃^- uptake rate in the Arabidopsis roots at a low NO₃^- level. Thus, we inferred that CsNRT2.4 was a key gene for N uptake in tea plant roots. This study provides new insights into the molecular mechanisms of tea plant responses to N resupply and reveals hub genes for improving nitrogen usage efficiency (NUE) in tea plants.

Proteometabolomic characterization of apical bud maturation in Pinus pinaster

Bud maturation is a physiological process that implies a set of morphophysiological changes that lead to the transition of growth patterns from young to mature. This transition defines tree growth and architecture, and in consequence traits such as biomass production and wood quality. In Pinus pinaster Aiton, a conifer of great timber value, bud maturation is closely related to polycyclism (multiple growth periods per year). This process causes a lack of apical dominance, and consequently increased branching that reduces its timber quality and value. However, despite its importance, little is known about bud maturation. In this work, proteomics and metabolomics were employed to study apical and basal sections of young and mature buds in P. pinaster. Proteins and metabolites in samples were described and quantified using (n)UPLC-LTQ-Orbitrap. The datasets were analyzed employing an integrative statistical approach, which allowed the determination of the interactions between proteins and metabolites and the different bud sections and ages. Specific dynamics of proteins and metabolites such as histones H3 and H4, ribosomal proteins L15 and L12, chaperonin TCP1, 14-3-3 protein gamma, gibberellins A1, A3 and A8, strigolactones and abscisic acid, involved in epigenetic regulation, proteome remodeling, hormonal signaling and abiotic stress pathways showed their potential role during bud maturation. Candidates and pathways were validated employing interaction databases and targeted transcriptomics. These results increase our understanding of the molecular processes behind bud maturation, a key step towards improving timber production and natural pine forests management in a future scenario of climate change. However, further studies are necessary using different P. pinaster populations that show contrasting wood quality and stress tolerance in order to generalize the results.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

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.

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.

Phylogenomics indicates the "living fossil" Isoetes diversified in the Cenozoic

The fossil record provides an invaluable insight into the temporal origins of extant lineages of organisms. However, establishing the relationships between fossils and extant lineages can be difficult in groups with low rates of morphological change over time. Molecular dating can potentially circumvent this issue by allowing distant fossils to act as calibration points, but rate variation across large evolutionary scales can bias such analyses. In this study, we apply multiple dating methods to genome-wide datasets to infer the origin of extant species of Isoetes, a group of mostly aquatic and semi-aquatic isoetalean lycopsids, which closely resemble fossil forms dating back to the Triassic. Rate variation observed in chloroplast genomes hampers accurate dating, but genome-wide nuclear markers place the origin of extant diversity within this group in the mid-Paleogene, 45-60 million years ago. Our genomic analyses coupled with a careful evaluation of the fossil record indicate that despite resembling forms from the Triassic, extant Isoetes species do not represent the remnants of an ancient and widespread group, but instead have spread around the globe in the relatively recent past.

The transcriptome of Pinus pinaster under Fusarium circinatum challenge

BACKGROUND

Fusarium circinatum, the causal agent of pitch canker disease, poses a serious threat to several Pinus species affecting plantations and nurseries. Although Pinus pinaster has shown moderate resistance to F. circinatum, the molecular mechanisms of defense in this host are still unknown. Phytohormones produced by the plant and by the pathogen are known to play a crucial role in determining the outcome of plant-pathogen interactions. Therefore, the aim of this study was to determine the role of phytohormones in F. circinatum virulence, that compromise host resistance.

RESULTS

A high quality P. pinaster de novo transcriptome assembly was generated, represented by 24,375 sequences from which 17,593 were full length genes, and utilized to determine the expression profiles of both organisms during the infection process at 3, 5 and 10 days post-inoculation using a dual RNA-sequencing approach. The moderate resistance shown by Pinus pinaster at the early time points may be explained by the expression profiles pertaining to early recognition of the pathogen, the induction of pathogenesis-related proteins and the activation of complex phytohormone signaling pathways that involves crosstalk between salicylic acid, jasmonic acid, ethylene and possibly auxins. Moreover, the expression of F. circinatum genes related to hormone biosynthesis suggests manipulation of the host phytohormone balance to its own benefit.

CONCLUSIONS

We hypothesize three key steps of host manipulation: perturbing ethylene homeostasis by fungal expression of genes related to ethylene biosynthesis, blocking jasmonic acid signaling by coronatine insensitive 1 (COI1) suppression, and preventing salicylic acid biosynthesis from the chorismate pathway by the synthesis of isochorismatase family hydrolase (ICSH) genes. These results warrant further testing in F. circinatum mutants to confirm the mechanism behind perturbing host phytohormone homeostasis.

Connections Between Amino Acid Metabolisms in Plants: Lysine as an Example

Extensive efforts have been made to fortify essential amino acids and boost nutrition in plants, but unintended effects on growth and physiology are also observed. Understanding how different amino acid metabolisms are connected with other biological pathways is therefore important. In addition to protein synthesis, amino acid metabolism is also tightly linked to energy and carbohydrate metabolism, the carbon-nitrogen budget, hormone and secondary metabolism, stress responses, and so on. Here, we update the currently available information on the connections between amino acid metabolisms, which tend to be overlooked in higher plants. Particular emphasis was placed on the connections between lysine metabolism and other pathways, such as tryptophan metabolism, the tricarboxylic acid cycle, abiotic and biotic stress responses, starch metabolism, and the unfolded protein response. Interestingly, regulation of lysine metabolism was found to differ between plant species, as is the case between dicots and monocots. Determining the metabolic connection between amino acid metabolisms will help improve our understanding of the metabolic flux, supporting studies on crop nutrition.

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.

Characterization and Application of EST-SSR Markers Developed From the Transcriptome of Amentotaxus argotaenia (Taxaceae), a Relict Vulnerable Conifer

Amentotaxus argotaenia (Taxaceae) is a vulnerable coniferous species with preference for shade and moist environment. Accurate estimation of genetic variation is crucial for its conservation, especially in the context of global warming. In this study, we acquired a transcriptome from A. argotaenia leaves using Illumina sequencing and de novo assembled 62,896 unigenes, of which 5510 EST-SSRs were detected. Twenty-two polymorphic EST-SSRs were successfully developed and further used to investigate genetic variation, linkage disequilibrium, and bottleneck signatures of A. argotaenia. The results showed that A. argotaenia had moderate genetic variation and high genetic differentiation, which may provide raw material to protect against climatic changes and accelerate local adaptation, respectively. No bottlenecks were found to occur in A. argotaenia. Our study not only showed that these EST markers are very effective in population genetic analysis but also lay a solid foundation for further investigating adaptive evolution and conservation strategies of A. argotaenia.

Utilization of Tissue Ploidy Level Variation in de Novo Transcriptome Assembly of Pinus sylvestris

Compared to angiosperms, gymnosperms lag behind in the availability of assembled and annotated genomes. Most genomic analyses in gymnosperms, especially conifer tree species, rely on the use of de novo assembled transcriptomes. However, the level of allelic redundancy and transcript fragmentation in these assembled transcriptomes, and their effect on downstream applications have not been fully investigated. Here, we assessed three assembly strategies for short-reads data, including the utility of haploid megagametophyte tissue during de novo assembly as single-allele guides, for six individuals and five different tissues in Pinus sylvestris. We then contrasted haploid and diploid tissue genotype calls obtained from the assembled transcriptomes to evaluate the extent of paralog mapping. The use of the haploid tissue during assembly increased its completeness without reducing the number of assembled transcripts. Our results suggest that current strategies that rely on available genomic resources as guidance to minimize allelic redundancy are less effective than the application of strategies that cluster redundant assembled transcripts. The strategy yielding the lowest levels of allelic redundancy among the assembled transcriptomes assessed here was the generation of SuperTranscripts with Lace followed by CD-HIT clustering. However, we still observed some levels of heterozygosity (multiple gene fragments per transcript reflecting allelic redundancy) in this assembled transcriptome on the haploid tissue, indicating that further filtering is required before using these assemblies for downstream applications. We discuss the influence of allelic redundancy when these reference transcriptomes are used to select regions for probe design of exome capture baits and for estimation of population genetic diversity.

A high-density exome capture genotype-by-sequencing panel for forestry breeding in Pinus radiata

Development of genome-wide resources for application in genomic selection or genome-wide association studies, in the absence of full reference genomes, present a challenge to the forestry industry, where longer breeding cycles could benefit from the accelerated selection possible through marker-based breeding value predictions. In particular, large conifer megagenomes require a strategy to reduce complexity, whilst ensuring genome-wide coverage is achieved. Using a transcriptome-based reference template, we have successfully developed a high density exome capture genotype-by-sequencing panel for radiata pine (Pinus radiata D.Don), capable of capturing in excess of 80,000 single nucleotide polymorphism (SNP) markers with a minor allele frequency above 0.03 in the population tested. This represents approximately 29,000 gene models from a core set of 48,914 probes. A set of 704 SNP markers capable of pedigree reconstruction and differentiating individual genotypes were tested within two full-sib mapping populations. While as few as 70 markers could reconstruct parentage in almost all cases, the impact of missing genotypes was noticeable in several offspring. Therefore, 60 sets of 110 randomly selected SNP markers were compared for both parentage reconstruction and clone differentiation. The performance in parentage reconstruction showed little variation over 60 iterations. However, there was notable variation in discriminatory power between closely related individuals, indicating a higher density SNP marker panel may be required to elucidate hidden relationships in complex pedigrees.

Small RNA profiling in Pinus pinaster reveals the transcriptome of developing seeds and highlights differences between zygotic and somatic embryos

Regulation of seed development by small non-coding RNAs (sRNAs) is an important mechanism controlling a crucial phase of the life cycle of seed plants. In this work, sRNAs from seed tissues (zygotic embryos and megagametophytes) and from somatic embryos of Pinus pinaster were analysed to identify putative regulators of seed/embryo development in conifers. In total, sixteen sRNA libraries covering several developmental stages were sequenced. We show that embryos and megagametophytes express a large population of 21-nt sRNAs and that substantial amounts of 24-nt sRNAs were also detected, especially in somatic embryos. A total of 215 conserved miRNAs, one third of which are conifer-specific, and 212 high-confidence novel miRNAs were annotated. MIR159, MIR171 and MIR394 families were found in embryos, but were greatly reduced in megagametophytes. Other families, like MIR397 and MIR408, predominated in somatic embryos and megagametophytes, suggesting their expression in somatic embryos is associated with in vitro conditions. Analysis of the predicted miRNA targets suggests that miRNA functions are relevant in several processes including transporter activity at the cotyledon-forming stage, and sulfur metabolism across several developmental stages. An important resource for studying conifer embryogenesis is made available here, which may also provide insightful clues for improving clonal propagation via somatic embryogenesis.

Completion of the cytosolic post-chorismate phenylalanine biosynthetic pathway in plants

In addition to being a vital component of proteins, phenylalanine is also a precursor of numerous aromatic primary and secondary metabolites with broad physiological functions. In plants phenylalanine is synthesized predominantly via the arogenate pathway in plastids. Here, we describe the structure, molecular players and subcellular localization of a microbial-like phenylpyruvate pathway for phenylalanine biosynthesis in plants. Using a reverse genetic approach and metabolic flux analysis, we provide evidence that the cytosolic chorismate mutase is responsible for directing carbon flux towards cytosolic phenylalanine production via the phenylpyruvate pathway. We also show that an alternative transcription start site of a known plastidial enzyme produces a functional cytosolic prephenate dehydratase that catalyzes the conversion of prephenate to phenylpyruvate, the intermediate step between chorismate mutase and phenylpyruvate aminotransferase. Thus, our results complete elucidation of phenylalanine biosynthesis via phenylpyruvate in plants, showing that this pathway splits from the known plastidial arogenate pathway at chorismate, instead of prephenate as previously thought, and the complete pathway is localized in the cytosol.

Completion of the cytosolic post-chorismate phenylalanine biosynthetic pathway in plants

In addition to being a vital component of proteins, phenylalanine is also a precursor of numerous aromatic primary and secondary metabolites with broad physiological functions. In plants phenylalanine is synthesized predominantly via the arogenate pathway in plastids. Here, we describe the structure, molecular players and subcellular localization of a microbial-like phenylpyruvate pathway for phenylalanine biosynthesis in plants. Using a reverse genetic approach and metabolic flux analysis, we provide evidence that the cytosolic chorismate mutase is responsible for directing carbon flux towards cytosolic phenylalanine production via the phenylpyruvate pathway. We also show that an alternative transcription start site of a known plastidial enzyme produces a functional cytosolic prephenate dehydratase that catalyzes the conversion of prephenate to phenylpyruvate, the intermediate step between chorismate mutase and phenylpyruvate aminotransferase. Thus, our results complete elucidation of phenylalanine biosynthesis via phenylpyruvate in plants, showing that this pathway splits from the known plastidial arogenate pathway at chorismate, instead of prephenate as previously thought, and the complete pathway is localized in the cytosol.

De novo transcriptome sequencing and SSR markers development for Cedrela balansae C.DC., a native tree species of northwest Argentina

The endangered Cedrela balansae C.DC. (Meliaceae) is a high-value timber species with great potential for forest plantations that inhabits the tropical forests in Northwestern Argentina.Research on this species is scarce because of the limited genetic and genomic information available. Here, we explored the transcriptome of C. balansae using 454 GS FLX Titanium next-generation sequencing (NGS) technology. Following de novo assembling, we identified 27,111 non-redundant unigenes longer than 200 bp, and considered these transcripts for further downstream analysis. The functional annotation was performed searching the 27,111 unigenes against the NR-Protein and the Interproscan databases. This analysis revealed 26,977 genes with homology in at least one of the Database analyzed. Furthermore, 7,774 unigenes in 142 different active biological pathways in C. balansae were identified with the KEGG database. Moreover, after in silico analyses, we detected 2,663 simple sequence repeats (SSRs) markers. A subset of 70 SSRs related to important "stress tolerance" traits based on functional annotation evidence, were selected for wet PCR-validation in C. balansae and other Cedrela species inhabiting in northwest and northeast of Argentina (C. fissilis, C. saltensis and C. angustifolia). Successful transferability was between 77% and 93% and thanks to this study, 32 polymorphic functional SSRs for all analyzed Cedrela species are now available. The gene catalog and molecular markers obtained here represent a starting point for further research, which will assist genetic breeding programs in the Cedrela genus and will contribute to identifying key populations for its preservation.

Transcriptome Comparison Reveals Distinct Selection Patterns in Domesticated and Wild Agave Species, the Important CAM Plants

Agave species are an important family of crassulacean acid metabolism (CAM) plants with remarkable tolerance to heat and drought stresses (Agave deserti) in arid regions and multiple agricultural applications, such as spirit (Agave tequilana) and fiber (Agave sisalana) production. The agave genomes are commonly too large to sequence, which has significantly restricted our understanding to the molecular basis of stress tolerance and economic traits in agaves. In this study, we collected three transcriptome databases for comparison to reveal the phylogenic relationships and evolution patterns of the three agave species. The results indicated the close but distinctly domesticated relations between A. tequilana and A. sisalana. Natural abiotic and biotic selections are very important factors that have contributed to distinct economic traits in agave domestication together with artificial selection. Besides, a series of candidate unigenes regulating fructan, fiber, and stress response-related traits were identified in A. tequilana, A. sisalana, and A. deserti, respectively. This study represents the first transcriptome comparison within domesticated and wild agaves, which would serve as a guidance for further studies on agave evolution, environmental adaptation, and improvement of economically important traits.

Identification and expression profiling of genes governing lignin biosynthesis in Casuarina equisetifolia L

Casuarina equisetifolia L. is an important multi-purpose, fast growing and widely planted tree species native to tropical and subtropical coastlines of Australia, Southeast Asia, Malaysia, Melanesia, Polynesia and New Caledonia. It is a nitrogen-fixing tree mainly used for charcoal making, construction poles, landscaping, timber, pulp, firewood, windbreaks, shelterbelts, soil erosion and sand dune stabilization. Casuarina wood is presently used for paper and pulp production. Raw material with reduced lignin is highly preferred to increase the pulp yield. Hence, understanding the molecular regulation of wood formation in this tree species is vital for selecting industrially suitable phenotypes for breeding programs. The lignin biosynthetic pathway has been extensively studied in tree species like Eucalypts, poplars, pines, Picea, Betula and Acacia sp. However, studies on wood formation at molecular level is presently lacking in casuarinas. Hence, in the present study, the transcriptome of the developing secondary tissues of 15 years old Casuarina equiseitfolia subsp. equisetifolia was sequenced, de novo assembled, annotated and mapped to functional pathways. Transcriptome sequencing generated a total of 26,985 transcripts mapped to 31 pathways. Mining of the annotated data identified nine genes involved in lignin biosynthesis pathway and relative expression of the transcripts in four tissues including scale-like leaves, needle-like brachlets, wood and root were documented. The expression of CeCCR1 and CeF5H were found to be significantly high in wood tissues, while maximum expression of CeHCT was documented in stem. Additionally, CeTUBA and CeH2A were identified as the most stable reference transcript for normalization of qRT-PCR data in C. equisetifolia. The present study is the first wood genomic resource in C. equisetifolia, which will be valuable for functional genomics research in this genus.

Approaches to variant discovery for conifer transcriptome sequencing

There is a wide diversity of bioinformatic tools available for the assembly of next generation sequence and subsequence variant calling to identify genetic markers at scale. Integration of genomics tools such as genomic selection, association studies, pedigree analysis and analysis of genetic diversity, into operational breeding is a goal for New Zealand’s most widely planted exotic tree species, Pinus radiata. In the absence of full reference genomes for large megagenomes such as in conifers, RNA sequencing in a range of genotypes and tissue types, offers a rich source of genetic markers for downstream application. We compared nine different assembler and variant calling software combinations in a single transcriptomic library and found that Single Nucleotide Polymorphism (SNPs) discovery could vary by as much as an order of magnitude (8,061 SNPs up to 86,815 SNPs). The assembler with the best realignment of the packages trialled, Trinity, in combination with several variant callers was then applied to a much larger multi-genotype, multi-tissue transcriptome and identified 683,135 in silico SNPs across a predicted 449,951 exons when mapped to the Pinus taeda ver 1.01e reference.

Application of EST-SSR markers developed from the transcriptome of Torreya grandis (Taxaceae), a threatened nut-yielding conifer tree

Torreya grandis (Taxaceae) is an ancient conifer species endemic to southeast China. Because of its nutrient-rich and delicious seeds, this species has been utilized for centuries by the Chinese. However, transcriptome data and transcriptome-derived microsatellite markers for population genetics studies are still insufficient for understanding of this species’ genetic basis. In this study, a transcriptome from T. grandis leaves was generated using Illumina sequencing. A total of 69,920 unigenes were generated after de novo assembly, and annotated by searching against seven protein databases. In addition, 2,065 expressed sequence tag–simple sequence repeats (EST-SSRs) were detected, with the distribution frequency of 2.75% of total unigenes and average number of 0.03 SSRs per unigene. Among these EST-SSRs, 1,339 primer pairs were successfully designed, and 106 primer pairs were randomly selected for the development of potential molecular markers. Among them, 11 EST-SSR markers revealed a moderate level of genetic diversity, and were used to investigate the population structure of T. grandis. Two different genetic groups within this species were revealed using these EST-SSR markers, indicating that these markers developed in this study can be effectively applied to the population genetic analysis of T. grandis.

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.

Transcriptome analysis of Pinus halepensis under drought stress and during recovery

Forest trees use various strategies to cope with drought stress and these strategies involve complex molecular mechanisms. Pinus halepensis Miller (Aleppo pine) is found throughout the Mediterranean basin and is one of the most drought-tolerant pine species. In order to decipher the molecular mechanisms that P. halepensis uses to withstand drought, we performed large-scale physiological and transcriptome analyses. We selected a mature tree from a semi-arid area with suboptimal growth conditions for clonal propagation through cuttings. We then used a high-throughput experimental system to continuously monitor whole-plant transpiration rates, stomatal conductance and the vapor pressure deficit. The transcriptomes of plants were examined at six physiological stages: pre-stomatal response, partial stomatal closure, minimum transpiration, post-irrigation, partial recovery and full recovery. At each stage, data from plants exposed to the drought treatment were compared with data collected from well-irrigated control plants. A drought-stressed P. halepensis transcriptome was created using paired-end RNA-seq. In total, ~6000 differentially expressed, non-redundant transcripts were identified between drought-treated and control trees. Cluster analysis has revealed stress-induced down-regulation of transcripts related to photosynthesis, reactive oxygen species (ROS)-scavenging through the ascorbic acid (AsA)-glutathione cycle, fatty acid and cell wall biosynthesis, stomatal activity, and the biosynthesis of flavonoids and terpenoids. Up-regulated processes included chlorophyll degradation, ROS-scavenging through AsA-independent thiol-mediated pathways, abscisic acid response and accumulation of heat shock proteins, thaumatin and exordium. Recovery from drought induced strong transcription of retrotransposons, especially the retrovirus-related transposon Tnt1-94. The drought-related transcriptome illustrates this species' dynamic response to drought and recovery and unravels novel mechanisms.

The role of arginine metabolic pathway during embryogenesis and germination in maritime pine (Pinus pinaster Ait.)

Vegetative propagation through somatic embryogenesis is critical in conifer biotechnology towards multivarietal forestry that uses elite varieties to cope with environmental and socio-economic issues. An important and still sub-optimal process during in vitro maturation of somatic embryos (SE) is the biosynthesis and deposition of storage proteins, which are rich in amino acids with high nitrogen (N) content, such as arginine. Mobilization of these N-rich proteins is essential for the germination and production of vigorous somatic seedlings. Somatic embryos accumulate lower levels of N reserves than zygotic embryos (ZE) at a similar stage of development. To understand the molecular basis for this difference, the arginine metabolic pathway has been characterized in maritime pine (Pinus pinaster Ait.). The genes involved in arginine metabolism have been identified and GFP-fusion constructs were used to locate the enzymes in different cellular compartments and clarify their metabolic roles during embryogenesis and germination. Analysis of gene expression during somatic embryo maturation revealed high levels of transcripts for genes involved in the biosynthesis and metabolic utilization of arginine. By contrast, enhanced expression levels were only observed during the last stages of maturation and germination of ZE, consistent with the adequate accumulation and mobilization of protein reserves. These results suggest that arginine metabolism is unbalanced in SE (simultaneous biosynthesis and degradation of arginine) and could explain the lower accumulation of storage proteins observed during the late stages of somatic embryogenesis.

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.

Glutamate synthases from conifers: gene structure and phylogenetic studies

BACKGROUND

Plants synthesize glutamate from ammonium by the combined activity of the enzymes glutamine synthetase (GS) and glutamate synthase (GOGAT) through the glutamate synthase cycle. In plants, there are two forms of glutamate synthases that differ in their electron donors, NADH-GOGAT (EC 1.4.1.14) and Fd-GOGAT (EC 1.4.7.1), which have differential roles either in primary ammonia assimilation or in the reassimilation of ammonium from different catabolic processes. Glutamate synthases are complex iron-sulfur flavoproteins containing functional domains involved in the control and coordination of their catalytic activities in annual plants. In conifers, partial cDNA sequences for GOGATs have been isolated and used for gene expression studies. However, knowledge of the gene structure and of phylogenetic relationships with other plant enzymes is quite scant.

RESULTS

Technological advances in conifer megagenomes sequencing have made it possible to obtain full-length cDNA sequences encoding Fd- and NADH-GOGAT from maritime pine, as well as BAC clones containing sequences for NADH-GOGAT and Fd-GOGAT genes. In the current study, we studied the genomic organization of pine GOGAT genes, the size of their exons/introns, copy numbers in the pine genome and relationships with other plant genes. Phylogenetic analysis was performed, and the degree of preservation and dissimilarity of key domains for the catalytic activities of these enzymes in different taxa were determined.

CONCLUSIONS

Fd- and NADH-GOGAT are encoded by single-copy genes in the maritime pine genome. The Fd-GOGAT gene is extremely large spanning more than 330 kb and the presence of very long introns highlights the important contribution of LTR retrotransposons to the gene size in conifers. In contrast, the structure of the NADH-GOGAT gene is similar to the orthologous genes in angiosperms. Our phylogenetic analysis indicates that these two genes had different origins during plant evolution. The results provide new insights into the structure and molecular evolution of these essential genes.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

The gene expression landscape of pine seedling tissues

Conifers dominate vast regions of the Northern hemisphere. They are the main source of raw materials for timber industry as well as a wide range of biomaterials. Despite their inherent difficulties as experimental models for classical plant biology research, the technological advances in genomics research are enabling fundamental studies on these plants. The use of laser capture microdissection followed by transcriptomic analysis is a powerful tool for unravelling the molecular and functional organization of conifer tissues and specialized cells. In the present work, 14 different tissues from 1-month-old maritime pine (Pinus pinaster) seedlings have been isolated and their transcriptomes analysed. The results increased the sequence information and number of full-length transcripts from a previous reference transcriptome and added 39 841 new transcripts. In total, 2376 transcripts were ubiquitously expressed in all of the examined tissues. These transcripts could be considered the core 'housekeeping genes' in pine. The genes have been clustered in function to their expression profiles. This analysis reduced the number of profiles to 38, most of these defined by their expression in a unique tissue that is much higher than in the other tissues. The expression and localization data are accessible at ConGenIE.org (http://v22.popgenie.org/microdisection/). This study presents an overview of the gene expression distribution in different pine tissues, specifically highlighting the relationships between tissue gene expression and function. This transcriptome atlas is a valuable resource for functional genomics research in conifers.

Molecular fundamentals of nitrogen uptake and transport in trees

Nitrogen (N) is frequently a limiting factor for tree growth and development. Because N availability is extremely low in forest soils, trees have evolved mechanisms to acquire and transport this essential nutrient along with biotic interactions to guarantee its strict economy. Here we review recent advances in the molecular basis of tree N nutrition. The molecular characteristics, regulation, and biological significance of membrane proteins involved in the uptake and transport of N are addressed. The regulation of N uptake and transport in mycorrhized roots and transcriptome-wide studies of N nutrition are also outlined. Finally, several areas of future research are suggested.

Identification of Norway Spruce MYB-bHLH-WDR Transcription Factor Complex Members Linked to Regulation of the Flavonoid Pathway

Transcription factors (TFs) forming MYB-bHLH-WDR complexes are known to regulate the biosynthesis of specialized metabolites in angiosperms through an intricate network. These specialized metabolites participate in a wide range of biological processes including plant growth, development, reproduction as well as in plant immunity. Studying the regulation of their biosynthesis is thus essential. While MYB (TFs) have been previously shown to control specialized metabolism (SM) in gymnosperms, the identity of their partners, in particular bHLH or WDR members, has not yet been revealed. To gain knowledge about MYB-bHLH-WDR transcription factor complexes in gymnosperms and their regulation of SW, we identified two bHLH homologs of AtTT8, six homologs of the MYB transcription factor AtTT2 and one WDR ortholog of AtTTG1 in Norway spruce. We investigated the expression levels of these genes in diverse tissues and upon treatments with various stimuli including methyl-salicylate, methyl-jasmonate, wounding or fungal inoculation. In addition, we also identified protein-protein interactions among different homologs of MYB, bHLH and WDR. Finally, we generated transgenic spruce cell lines overexpressing four of the Norway spruce AtTT2 homologs and observed differential regulation of genes in the flavonoid pathway and flavonoid contents.

Characterization of Three L-Asparaginases from Maritime Pine (Pinus pinaster Ait.)

Asparaginases (ASPG, EC 3.5.1.1) catalyze the hydrolysis of the amide group of L-asparagine producing L-aspartate and ammonium. Three ASPG, PpASPG1, PpASPG2, and PpASPG3, have been identified in the transcriptome of maritime pine (Pinus pinaster Ait.) that were transiently expressed in Nicotiana benthamiana by agroinfection. The three recombinant proteins were processed in planta to active enzymes and it was found that all mature forms exhibited double activity asparaginase/isoaspartyl dipeptidase but only PpASPG1 was able to catalyze efficiently L-asparagine hydrolysis. PpASPG1 contains a variable region of 77 amino acids that is critical for proteolytic processing of the precursor and is retained in the mature enzyme. Furthermore, the functional analysis of deletion mutants demonstrated that this protein fragment is required for specific recognition of the substrate and favors enzyme stability. Potassium has a limited effect on the activation of maritime pine ASPG what is consistent with the lack of a critical residue essential for interaction of cation. Taken together, the results presented here highlight the specific features of ASPG from conifers when compared to the enzymes from angiosperms.

Integrated analysis of 454 and Illumina transcriptomic sequencing characterizes carbon flux and energy source for fatty acid synthesis in developing Lindera glauca fruits for woody biodiesel

BACKGROUND

Lindera glauca fruit with high quality and quantity of oil has emerged as a novel potential source of biodiesel in China, but the molecular regulatory mechanism of carbon flux and energy source for oil biosynthesis in developing fruits is still unknown. To better develop fruit oils of L. glauca as woody biodiesel, a combination of two different sequencing platforms (454 and Illumina) and qRT-PCR analysis was used to define a minimal reference transcriptome of developing L. glauca fruits, and to construct carbon and energy metabolic model for regulation of carbon partitioning and energy supply for FA biosynthesis and oil accumulation.

RESULTS

We first analyzed the dynamic patterns of growth tendency, oil content, FA compositions, biodiesel properties, and the contents of ATP and pyridine nucleotide of L. glauca fruits from seven different developing stages. Comprehensive characterization of transcriptome of the developing L. glauca fruit was performed using a combination of two different next-generation sequencing platforms, of which three representative fruit samples (50, 125, and 150 DAF) and one mixed sample from seven developing stages were selected for Illumina and 454 sequencing, respectively. The unigenes separately obtained from long and short reads (201, and 259, respectively, in total) were reconciled using TGICL software, resulting in a total of 60,031 unigenes (mean length = 1061.95 bp) to describe a transcriptome for developing L. glauca fruits. Notably, 198 genes were annotated for photosynthesis, sucrose cleavage, carbon allocation, metabolite transport, acetyl-CoA formation, oil synthesis, and energy metabolism, among which some specific transporters, transcription factors, and enzymes were identified to be implicated in carbon partitioning and energy source for oil synthesis by an integrated analysis of transcriptomic sequencing and qRT-PCR. Importantly, the carbon and energy metabolic model was well established for oil biosynthesis of developing L. glauca fruits, which could help to reveal the molecular regulatory mechanism of the increased oil production in developing fruits.

CONCLUSIONS

This study presents for the first time the application of an integrated two different sequencing analyses (Illumina and 454) and qRT-PCR detection to define a minimal reference transcriptome for developing L. glauca fruits, and to elucidate the molecular regulatory mechanism of carbon flux control and energy provision for oil synthesis. Our results will provide a valuable resource for future fundamental and applied research on the woody biodiesel plants.

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

Background

Increasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.

Results

The populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.

Conclusions

This study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping

Background

Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola.

Results

Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed.

Conclusions

We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3079-2) contains supplementary material, which is available to authorized users.

Performance of genomic prediction within and across generations in maritime pine

BACKGROUND

Genomic selection (GS) is a promising approach for decreasing breeding cycle length in forest trees. Assessment of progeny performance and of the prediction accuracy of GS models over generations is therefore a key issue.

RESULTS

A reference population of maritime pine (Pinus pinaster) with an estimated effective inbreeding population size (status number) of 25 was first selected with simulated data. This reference population (n = 818) covered three generations (G0, G1 and G2) and was genotyped with 4436 single-nucleotide polymorphism (SNP) markers. We evaluated the effects on prediction accuracy of both the relatedness between the calibration and validation sets and validation on the basis of progeny performance. Pedigree-based (best linear unbiased prediction, ABLUP) and marker-based (genomic BLUP and Bayesian LASSO) models were used to predict breeding values for three different traits: circumference, height and stem straightness. On average, the ABLUP model outperformed genomic prediction models, with a maximum difference in prediction accuracies of 0.12, depending on the trait and the validation method. A mean difference in prediction accuracy of 0.17 was found between validation methods differing in terms of relatedness. Including the progenitors in the calibration set reduced this difference in prediction accuracy to 0.03. When only genotypes from the G0 and G1 generations were used in the calibration set and genotypes from G2 were used in the validation set (progeny validation), prediction accuracies ranged from 0.70 to 0.85.

CONCLUSIONS

This study suggests that the training of prediction models on parental populations can predict the genetic merit of the progeny with high accuracy: an encouraging result for the implementation of GS in the maritime pine breeding program.

Performance of genomic prediction within and across generations in maritime pine

BACKGROUND

Genomic selection (GS) is a promising approach for decreasing breeding cycle length in forest trees. Assessment of progeny performance and of the prediction accuracy of GS models over generations is therefore a key issue.

RESULTS

A reference population of maritime pine (Pinus pinaster) with an estimated effective inbreeding population size (status number) of 25 was first selected with simulated data. This reference population (n = 818) covered three generations (G0, G1 and G2) and was genotyped with 4436 single-nucleotide polymorphism (SNP) markers. We evaluated the effects on prediction accuracy of both the relatedness between the calibration and validation sets and validation on the basis of progeny performance. Pedigree-based (best linear unbiased prediction, ABLUP) and marker-based (genomic BLUP and Bayesian LASSO) models were used to predict breeding values for three different traits: circumference, height and stem straightness. On average, the ABLUP model outperformed genomic prediction models, with a maximum difference in prediction accuracies of 0.12, depending on the trait and the validation method. A mean difference in prediction accuracy of 0.17 was found between validation methods differing in terms of relatedness. Including the progenitors in the calibration set reduced this difference in prediction accuracy to 0.03. When only genotypes from the G0 and G1 generations were used in the calibration set and genotypes from G2 were used in the validation set (progeny validation), prediction accuracies ranged from 0.70 to 0.85.

CONCLUSIONS

This study suggests that the training of prediction models on parental populations can predict the genetic merit of the progeny with high accuracy: an encouraging result for the implementation of GS in the maritime pine breeding program.

Deciphering the molecular basis of ammonium uptake and transport in maritime pine

Ammonium is the predominant form of inorganic nitrogen in the soil of coniferous forests. Despite the ecological and economic importance of conifers, the molecular basis of ammonium uptake and transport in this group of gymnosperms is largely unknown. In this study, we describe the functional characterization of members of the AMT gene family in Pinus pinaster: PpAMT1.1, PpAMT1.2 and PpAMT1.3 (subfamily 1) and PpAMT2.1 and PpAMT2.3 (subfamily 2). Our phylogenetic analysis indicates that in conifers, all members of the AMT1 subfamily evolved from a common ancestor that is evolutionarily related to the ancient PpAMT1.2 gene. Individual AMT genes are developmentally and nutritionally regulated, and their transcripts are specifically distributed in different organs. PpAMT1.3 was predominantly expressed in the roots, particularly during N starvation and mycorrhizal interaction, whereas PpAMT2.3 was preferentially expressed in lateral roots. Immunolocalization studies of roots with varied nitrogen availability revealed that PpAMT1 and PpAMT2 proteins play complementary roles in the uptake of external ammonium. Heterologous expression in yeast and Xenopus oocytes revealed that the AMT genes encode functional transporters with different kinetics and with different capacities for ammonium transport. Our results provide new insights on how nitrogen is acquired and transported in conifers.

Identification of a small protein domain present in all plant lineages that confers high prephenate dehydratase activity

l-Phenylalanine serves as a building block for the biosynthesis of proteins, but also as a precursor for a wide range of plant-derived compounds essential for plants and animals. Plants can synthesize Phe within the plastids using arogenate as a precursor; however, an alternative pathway using phenylpyruvate as an intermediate, described for most microorganisms, has recently been proposed. The functionality of this pathway requires the existence of enzymes with prephenate dehydratase (PDT) activity (EC 4.2.1.51) in plants. Using phylogenetic studies, functional complementation assays in yeast and biochemical analysis, we have identified the enzymes displaying PDT activity in Pinus pinaster. Through sequence alignment comparisons and site-directed mutagenesis we have identified a 22-amino acid region conferring PDT activity (PAC domain) and a single Ala314 residue critical to trigger this activity. Our results demonstrate that all plant clades include PAC domain-containing ADTs, suggesting that the PDT activity, and thus the ability to synthesize Phe using phenylpyruvate as an intermediate, has been preserved throughout the evolution of plants. Moreover, this pathway together with the arogenate pathway gives plants a broad and versatile capacity to synthesize Phe and its derived compounds. PAC domain-containing enzymes are also present in green and red algae, and glaucophytes, the three emerging clades following the primary endosymbiont event resulting in the acquisition of plastids in eukaryotes. The evolutionary prokaryotic origin of this domain is discussed.

Dataset of UV induced changes in nuclear proteome obtained by GeLC-Orbitrap/MS in Pinus radiata needles

Although responses to UV stress have been characterised at system and cellular levels, the dynamics of the nuclear proteome triggered in this situation are still unknown, despite its essential role in regulating gene expression and in last term plant physiology. To fill this gap, we characterised the variations in the nuclear proteome after 2 h and 16 h (8 h/day) of UV irradiation by using state-of-the-art mass spectrometry-based shotgun proteomics methods combined with novel bioinformatics workflows that were employed in the manuscript entitled "The variations in the nuclear proteome reveal new transcription factors and mechanisms involved in UV stress response in Pinus radiata" (Pascual et al., 2016) [1]. We employed in-gel digestion followed by a 120 min gradient prior to MS analysis. Data was processed following two approaches: a database dependent employing the SEQUEST algorithm and custom databases, and a database independent by mass accuracy precursor alignment (MAPA). 388 proteins were identified by SEQUEST search and 9094 m/z were quantified by MAPA. Significant m/z were de novo sequenced using the Novor algorithm. We present here the complete datasets and the analysis workflow.