Identification of genes differentially expressed during adventitious shoot induction in Pinus pinea cotyledons by subtractive hybridization and quantitative PCR

Stone Pine (Pinus pinea L.) High-Added-Value Genetics: An Overview

Stone pine (Pinus pinea L.) has received limited attention in terms of genetic research. However, genomic techniques hold promise for decoding the stone pine genome and contributing to developing a more resilient bioeconomy. Retrotransposon and specific genetic markers are effective tools for determining population-specific genomic diversity. Studies on the transcriptome and proteome have identified differentially expressed genes PAS1, CLV1, ATAF1, and ACBF involved in shoot bud formation. The stone pine proteome shows variation among populations and shows the industrial potential of the enzyme pinosylvin. Microsatellite studies have revealed low levels of polymorphism and a unique genetic diversity in stone pine, which may contribute to its environmental adaptation. Transcriptomic and proteomic analyses uncover the genetic and molecular responses of stone pine to fungal infections and nematode infestations, elucidating the defense activation, gene regulation, and the potential role of terpenes in pathogen resistance. Transcriptomics associated with carbohydrate metabolism, dehydrins, and transcription factors show promise as targets for improving stone pine's drought stress response and water retention capabilities. Stone pine presents itself as an important model tree for studying climate change adaptation due to its characteristics. While knowledge gaps exist, stone pine's genetic resources hold significant potential, and ongoing advancements in techniques offer prospects for future exploration.

In Vitro Plant Regeneration in Conifers: The Role of WOX and KNOX Gene Families

Conifers are a group of woody plants with an enormous economic and ecological importance. Breeding programs are necessary to select superior varieties for planting, but they have many limitations due to the biological characteristics of conifers. Somatic embryogenesis (SE) and de novo organogenesis (DNO) from in vitro cultured tissues are two ways of plant mass propagation that help to overcome this problem. Although both processes are difficult to achieve in conifers, they offer advantages like a great efficiency, the possibilities to cryopreserve the embryogenic lines, and the ability of multiplying adult trees (the main bottleneck in conifer cloning) through DNO. Moreover, SE and DNO represent appropriate experimental systems to study the molecular bases of developmental processes in conifers such as embryogenesis and shoot apical meristem (SAM) establishment. Some of the key genes regulating these processes belong to the WOX and KNOX homeobox gene families, whose function has been widely described in Arabidopsis thaliana. The sequences and roles of these genes in conifers are similar to those found in angiosperms, but some particularities exist, like the presence of WOXX, a gene that putatively participates in the establishment of SAM in somatic embryos and plantlets of Pinus pinaster.

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.

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.

Hormonal and gene dynamics in de novo shoot meristem formation during adventitious caulogenesis in cotyledons of Pinus pinea

KEY MESSAGE

Several members of WOX and KNOX gene families and several plant growth regulators, basically cytokinins and auxins, play a key role during adventitious caulogenesis in the conifer Pinus pinea. Similar to Arabidopsis thaliana, Pinus pinea shoot organogenesis is a multistep process. However, there are key differences between both species, which may alter the underlying physiological and genetic programs. It is unknown if the genic expression models during angiosperm development may be applicable to conifers. In this work, an analysis of the endogenous content of different plant growth regulators and the expression of genes putatively involved in adventitious caulogenesis in P. pinea cotyledons was conducted. A multivariate analysis of both datasets was also realized through partial least squares regression and principal component analysis to obtain an integral vision of the mechanisms involved in caulogenesis in P. pinea. Analyses show that cotyledons cultured in the presence of benzyladenine during long times (2-6 days) cluster separately from the rest of the samples, suggesting that the benzyladenine increase observed during the first hours of culture is sufficient to trigger the caulogenic response through the activation of specific developmental programs. In particular, the most relevant factors involved in this process are the cytokinins trans-zeatin, dihydrozeatin, trans-zeatin riboside and isopentenyl adenosine; the auxin indoleacetic acid; and the genes PpWUS, PpWOX5, PpKN2, PpKN3 and PipiRR1. WUS is functional in pines and has an important role in caulogenesis. Interestingly, WOX5 also seems to participate in the process, although its specific role has not been determined.

Comparative transcriptional and metabolic responses of Pinus pinea to a native and a non-native Heterobasidion species

Heterobasidion irregulare is a causal agent of root and butt-rot disease in conifers, and is native to North America. In 1944 it was introduced in central Italy in a Pinus pinea stand, where it shares the same niche with the native species Heterobasidion annosum. The introduction of a non-native pathogen may have significant negative effects on a naïve host tree and the ecosystem in which it resides, requiring a better understanding of the system. We compared the spatio-temporal phenotypic, transcriptional and metabolic host responses to inoculation with the two Heterobasidion species in a large experiment with P. pinea seedlings. Differences in length of lesions at the inoculation site (IS), expression of host genes involved in lignin pathway and in cell rescue and defence, and analysis of terpenes at both IS and 12 cm above the IS (distal site, DS), were assessed at 3, 14 and 35 days post inoculation (dpi). Results clearly showed that both species elicit similar physiological and biochemical responses in P. pinea seedlings. The analysis of host transcripts and total terpenes showed differences between inoculation sites and between pathogen and mock inoculated plants. Both pathogen and mock inoculations induced antimicrobial peptide and phenylalanine ammonia-lyase overexpression at IS beginning at 3 dpi; while at DS all the analysed genes, except for peroxidase, were overexpressed at 14 dpi. A significantly higher accumulation of terpenoids was observed at 14 dpi at IS, and at 35 dpi at DS. The terpene blend at IS showed significant variation among treatments and sampling times, while no significant differences were ever observed in DS tissues. Based on our results, H. irregulare does not seem to have competitive advantages over the native species H. annosum in terms of pathogenicity towards P. pinea trees; this may explain why the non-native species has not widely spread over the 73 years since its putative year of introduction into central Italy.

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.

Selection of Reference Genes for Real-Time Quantitative PCR in Pinus massoniana Post Nematode Inoculation

Pinus massoniaia Lamb has gained more and more attention as the most important tree species for timber and forestation in South China. Gene expression studies are of great importance to identify new and elite cultivars. Real-time quantitative PCR, a highly sensitive and specific method, is commonly used in the analysis of gene expression. The appropriate reference genes must be employed to normalize the calculation program for ascertaining repeatable and significant results. Herein, eleven housekeeping genes were evaluated during different stages of P. massoniana post nematode inoculation in this study. Three statistical approaches such as geNorm, NormFinder and BestKeeper were selected to analyze the stability of candidate genes. The results indicated that U2af and β-TUB were the most stable reference genes. These two genes could be used for the normalization in most of the experiments of P. massoniana, while Histone and AK were the least stable ones. In addition, EF expressed at the lowest average Ct value was the most abundant candidate gene. As an important gene associated with defense mechanisms, ABC transporter was analyzed by qRT-PCR, and the results were used to confirm the reliability of two genes. The selected reference genes in the present study will be conducive to future gene expression normalized by qRT-PCR in P. massoniana.

High-density SNP assay development for genetic analysis in maritime pine (Pinus pinaster)

Maritime pine provides essential ecosystem services in the south-western Mediterranean basin, where it covers around 4 million ha. Its scattered distribution over a range of environmental conditions makes it an ideal forest tree species for studies of local adaptation and evolutionary responses to climatic change. Highly multiplexed single nucleotide polymorphism (SNP) genotyping arrays are increasingly used to study genetic variation in living organisms and for practical applications in plant and animal breeding and genetic resource conservation. We developed a 9k Illumina Infinium SNP array and genotyped maritime pine trees from (i) a three-generation inbred (F2) pedigree, (ii) the French breeding population and (iii) natural populations from Portugal and the French Atlantic coast. A large proportion of the exploitable SNPs (2052/8410, i.e. 24.4%) segregated in the mapping population and could be mapped, providing the densest ever gene-based linkage map for this species. Based on 5016 SNPs, natural and breeding populations from the French gene pool exhibited similar level of genetic diversity. Population genetics and structure analyses based on 3981 SNP markers common to the Portuguese and French gene pools revealed high levels of differentiation, leading to the identification of a set of highly differentiated SNPs that could be used for seed provenance certification. Finally, we discuss how the validated SNPs could facilitate the identification of ecologically and economically relevant genes in this species, improving our understanding of the demography and selective forces shaping its natural genetic diversity, and providing support for new breeding strategies.

Normalizing gene expression by quantitative PCR during somatic embryogenesis in two representative conifer species: Pinus pinaster and Picea abies

Suitable internal control genes to normalize qPCR data from different stages of embryo development and germination were identified in two representative conifer species. Clonal propagation by somatic embryogenesis has a great application potentiality in conifers. Quantitative PCR (qPCR) is widely used for gene expression analysis during somatic embryogenesis and embryo germination. No single reference gene is universal, so a systematic characterization of endogenous genes for concrete conditions is fundamental for accuracy. We identified suitable internal control genes to normalize qPCR data obtained at different steps of somatic embryogenesis (embryonal mass proliferation, embryo maturation and germination) in two representative conifer species, Pinus pinaster and Picea abies. Candidate genes included endogenous genes commonly used in conifers, genes previously tested in model plants, and genes with a lower variation of the expression along embryo development according to genome-wide transcript profiling studies. Three different algorithms were used to evaluate expression stability. The geometric average of the expression values of elongation factor-1α, α-tubulin and histone 3 in P. pinaster, and elongation factor-1α, α-tubulin, adenosine kinase and CAC in P. abies were adequate for expression studies throughout somatic embryogenesis. However, improved accuracy was achieved when using other gene combinations in experiments with samples at a single developmental stage. The importance of studies selecting reference genes to use in different tissues or developmental stages within one or close species, and the instability of commonly used reference genes, is highlighted.

Characterization of a type-A response regulator differentially expressed during adventitious caulogenesis in Pinus pinaster

The molecular cloning and characterization of PipsRR1, a type-A response regulator in Pinus pinaster, is reported here. Type-A response regulators mediate downstream responses to cytokinin and act as negative feedback regulators of the signal transduction pathway. Some type-A response regulators in Arabidopsis have been related to de novo meristem formation. However, little information exists in Pinus spp. The PipsRR1 gene contains 5 exons, as do all type-A response regulators in Arabidopsis, and the deduced protein contains a receiver domain with the conserved DDK residues and a short C terminal extension. Expression analysis showed that the PipsRR1 gene is differentially expressed during the first phases of adventitious caulogenesis induced by benzyladenine in P. pinaster cotyledons, suggesting that PipsRR1 plays a role in caulogenesis in conifers. Additionally, a binary vector carrying the PipsRR1 promoter driving GFP:GUS expression was constructed to analyze the promoter activity in P. pinaster somatic embryos. The results of genetic transformation showed GUS activity during somatic embryo mass proliferation and embryo maturation.

EuroPineDB: a high-coverage web database for maritime pine transcriptome

Background

Pinus pinaster is an economically and ecologically important species that is becoming a woody gymnosperm model. Its enormous genome size makes whole-genome sequencing approaches are hard to apply. Therefore, the expressed portion of the genome has to be characterised and the results and annotations have to be stored in dedicated databases.

Description

EuroPineDB is the largest sequence collection available for a single pine species, Pinus pinaster (maritime pine), since it comprises 951 641 raw sequence reads obtained from non-normalised cDNA libraries and high-throughput sequencing from adult (xylem, phloem, roots, stem, needles, cones, strobili) and embryonic (germinated embryos, buds, callus) maritime pine tissues. Using open-source tools, sequences were optimally pre-processed, assembled, and extensively annotated (GO, EC and KEGG terms, descriptions, SNPs, SSRs, ORFs and InterPro codes). As a result, a 10.5× P. pinaster genome was covered and assembled in 55 322 UniGenes. A total of 32 919 (59.5%) of P. pinaster UniGenes were annotated with at least one description, revealing at least 18 466 different genes. The complete database, which is designed to be scalable, maintainable, and expandable, is freely available at: http://www.scbi.uma.es/pindb/. It can be retrieved by gene libraries, pine species, annotations, UniGenes and microarrays (i.e., the sequences are distributed in two-colour microarrays; this is the only conifer database that provides this information) and will be periodically updated. Small assemblies can be viewed using a dedicated visualisation tool that connects them with SNPs. Any sequence or annotation set shown on-screen can be downloaded. Retrieval mechanisms for sequences and gene annotations are provided.

Conclusions

The EuroPineDB with its integrated information can be used to reveal new knowledge, offers an easy-to-use collection of information to directly support experimental work (including microarray hybridisation), and provides deeper knowledge on the maritime pine transcriptome.

Expression profiling of a complex thaumatin-like protein family in western white pine

The protein content in the plant apoplast is believed to change dramatically as a result of host defense response upon infection with various pathogens. In this study, six novel thaumatin-like proteins (TLPs) were identified in western white pine (Pinus monticola) needle apoplast by a proteomic strategy using two-dimensional protein electrophoresis followed by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. Sequent cDNA cloning found that ten P. monticola TLP genes (PmTLP-L1 to -L6 and -S1 to -S4) were expressed in various tissues. Phylogenetic analysis demonstrated that these PmTLP genes belong to a large, complex, and highly diverse plant TLP family. Quantitative reverse transcription-polymerase chain reaction (RT-qPCR) using gene-specific primer pairs showed that each PmTLP gene exhibited a characteristic pattern of mRNA expression based on their unique organ distribution, seasonal regulation, and response to abiotic and biotic stresses. A time-course analysis at the early stages of infection by white pine blister rust pathogen Cronartium ribicola revealed that a coordinated upregulation of multiple PmTLP genes was involved in P. monticola major gene (Cr2) resistance. The structural and expressional differentiations suggest that the PmTLP family may contribute to host defense as well as other mechanism.

Caulogenic induction in cotyledons of stone pine (Pinus pinea): relationship between organogenic response and benzyladenine trends in selected families

Adventitious bud formation in stone pine cotyledons cultured in the presence of benzyladenine (BA) has been proposed as a model for the study of in vitro shoot organogenesis in conifers. This is because of its advantageous characteristics including the requirement of only one plant growth regulator (BA), the synchronous fashion of its induction, and the homogeneity and low degree of differentiation of cotyledons. Although optimal culture conditions have been developed and are currently in use, we still lack data for BA dynamics in cotyledons cultured under these conditions, and the morphological description of the early induction stages has not, until now, been approached from a histological perspective. Consequently, this is the focus of the present report. Additionally, we examined uptake and metabolism of BA in cotyledons from two selected families, previously characterized by, and selected for, the difference in the magnitude of their organogenic response. Media transfer experiments established that cotyledons should be in contact with 44.4 microM BA for at least 6h to obtain any caulogenic response (minimum shoot-induction period). Histological observations, carried out here for the first time in this species, determined that meristemoid structures had already begun to appear in explants within 12 h of culture. Moreover, results from the BA uptake and metabolism experiments indicated that the point at which explants reached the maximum concentration of active forms of BA (276.60 microM at 6 h) and the onset of the determination phase of shoot organogenesis were directly related. A direct relationship was also observed between the intensity of the caulogenic response in cotyledons from families 36 and 61 and the endogenous concentration of BA and its riboside at the start of the induction phase. Hence, family 36, characterized by its higher bud production, reached concentrations of 251.56 microM, while family 61, selected for its low bud-producing trait, only attained 175.80 microM. Finally, a correlation was observed between 6-benzylamino-9-[O-glucopyranosyl-(1-->3)ribofuranosyl]-purine values and the magnitude of the shoot organogenesis response.