Apparent homology of expressed genes from wood-forming tissues of loblolly pine (Pinus taeda L.) with Arabidopsis thaliana

Haplotype-resolved genome assembly and resequencing analysis provide insights into genome evolution and allelic imbalance in Pinus densiflora

Haplotype-level allelic characterization facilitates research on the functional, evolutionary and breeding-related features of extremely large and complex plant genomes. We report a 21.7-Gb chromosome-level haplotype-resolved assembly in Pinus densiflora. We found genome rearrangements involving translocations and inversions between chromosomes 1 and 3 of Pinus species and a proliferation of specific long terminal repeat (LTR) retrotransposons (LTR-RTs) in P. densiflora. Evolutionary analyses illustrated that tandem and LTR-RT-mediated duplications led to an increment of transcription factor (TF) genes in P. densiflora. The haplotype sequence comparison showed allelic imbalances, including presence-absence variations of genes (PAV genes) and their functional contributions to flowering and abiotic stress-related traits in P. densiflora. Allele-aware resequencing analysis revealed PAV gene diversity across P. densiflora accessions. Our study provides insights into key mechanisms underlying the evolution of genome structure, LTR-RTs and TFs within the Pinus lineage as well as allelic imbalances and diversity across P. densiflora.

Seasonal Developing Xylem Transcriptome Analysis of Pinus densiflora Unveils Novel Insights for Compression Wood Formation

Wood is the most important renewable resource not only for numerous practical utilizations but also for mitigating the global climate crisis by sequestering atmospheric carbon dioxide. The compressed wood (CW) of gymnosperms, such as conifers, plays a pivotal role in determining the structure of the tree through the reorientation of stems displaced by environmental forces and is characterized by a high content of lignin. Despite extensive studies on many genes involved in wood formation, the molecular mechanisms underlying seasonal and, particularly, CW formation remain unclear. This study examined the seasonal dynamics of two wood tissue types in Pinus densiflora: CW and opposite wood (OW). RNA sequencing of developing xylem for two consecutive years revealed comprehensive transcriptome changes and unique differences in CW and OW across seasons. During growth periods, such as spring and summer, we identified 2255 transcripts with differential expression in CW, with an upregulation in lignin biosynthesis genes and significant downregulation in stress response genes. Notably, among the laccases critical for monolignol polymerization, PdeLAC17 was found to be specifically expressed in CW, suggesting its vital role in CW formation. PdeERF4, an ERF transcription factor preferentially expressed in CW, seems to regulate PdeLAC17 activity. This research provides an initial insight into the transcriptional regulation of seasonal CW development in P. densiflora, forming a foundation for future studies to enhance our comprehension of wood formation in gymnosperms.

Decoding exceptional plant traits by comparative single-cell genomics

Some plants have acquired traits of remarkable adaptive value to thrive under stress. Transferring these traits to crops could improve agriculture, but uncovering the toolkit required has remained largely elusive. We propose that single-cell genomics offers a framework to compare species with contrasting developmental traits and to identify the regulators of evolutionary innovations.

Genome-Wide Analysis and Expression Profiling of Glutathione Reductase Gene Family in Oat (Avena sativa) Indicate Their Responses to Abiotic Stress during Seed Imbibition

Abiotic stress disturbs plant cellular redox homeostasis, inhibiting seed germination and plant growth. This is a crucial limitation to crop yield. Glutathione reductase (GR) is an important component of the ascorbate-glutathione (AsA-GSH) cycle which is involved in multiple plant metabolic processes. In the present study, GRs in A. sativa (AsGRs) were selected to explore their molecular characterization, phylogenetic relationship, and RNA expression changes during seed imbibition under abiotic stress. Seven AsGR genes were identified and mapped on six chromosomes of A, C, and D subgenomes. Phylogenetic analysis and subcellular localization of AsGR proteins divided them into two sub-families, AsGR1 and AsGR2, which were predicted to be mainly located in cytoplasm, mitochondrion, and chloroplast. Cis-elements relevant to stress and hormone responses are distributed in promoter regions of AsGRs. Tissue-specific expression profiling showed that AsGR1 genes were highly expressed in roots, leaves, and seeds, while AsGR2 genes were highly expressed in leaves and seeds. Both AsGR1 and AsGR2 genes showed a decreasing-increasing expression trend during seed germination under non-stress conditions. In addition, their responses to drought, salt, cold, copper, H₂O₂, and ageing treatments were quite different during seed imbibition. Among the seven AsGR genes, AsGR1-A, AsGR1-C, AsGR2-A, and AsGR2-D responded more significantly, especially under drought, ageing, and H₂O₂ stress. This study has laid the ground for the functional characterization of GR and the improvement of oat stress tolerance and seed vigor.

Wood Formation under Changing Environment: Omics Approaches to Elucidate the Mechanisms Driving the Early-to-Latewood Transition in Conifers

The global change scenarios highlight the urgency of clarifying the mechanisms driving the determination of wood traits in forest trees. Coniferous xylem is characterized by the alternation between earlywood (EW) and latewood (LW), on which proportions the wood density depend, one of the most important mechanical xylem qualities. However, the molecular mechanisms triggering the transition between the production of cells with the typical features of EW to the LW are still far from being completely elucidated. The increasing availability of omics resources for conifers, e.g., genomes and transcriptomes, would lay the basis for the comprehension of wood formation dynamics, boosting both breeding and gene-editing approaches. This review is intended to introduce the importance of wood formation dynamics and xylem traits of conifers in a changing environment. Then, an up-to-date overview of the omics resources available for conifers was reported, focusing on both genomes and transcriptomes. Later, an analysis of wood formation studies using omics approaches was conducted, with the aim of elucidating the main metabolic pathways involved in EW and LW determination. Finally, the future perspectives and the urgent needs on this research topic were highlighted.

Functional and comparative genomics reveals conserved noncoding sequences in the nitrogen-fixing clade

Nitrogen is one of the most inaccessible plant nutrients, but certain species have overcome this limitation by establishing symbiotic interactions with nitrogen-fixing bacteria in the root nodule. This root-nodule symbiosis (RNS) is restricted to species within a single clade of angiosperms, suggesting a critical, but undetermined, evolutionary event at the base of this clade. To identify putative regulatory sequences implicated in the evolution of RNS, we evaluated the genomes of 25 species capable of nodulation and identified 3091 conserved noncoding sequences (CNS) in the nitrogen-fixing clade (NFC). We show that the chromatin accessibility of 452 CNS correlates significantly with the regulation of genes responding to lipochitooligosaccharides in Medicago truncatula. These included 38 CNS in proximity to 19 known genes involved in RNS. Five such regions are upstream of MtCRE1, Cytokinin Response Element 1, required to activate a suite of downstream transcription factors necessary for nodulation in M. truncatula. Genetic complementation of an Mtcre1 mutant showed a significant decrease of nodulation in the absence of the five CNS, when they are driving the expression of a functional copy of MtCRE1. CNS identified in the NFC may harbor elements required for the regulation of genes controlling RNS in M. truncatula.

Seasonal changes in cambium activity from active to dormant stage affect the formation of secondary xylem in Pinus tabulaeformis Carr

Understanding the changing patterns of vascular cambium during seasonal cycles is crucial to reveal the mechanisms that control cambium activity and wood formation, but this area has been underexplored, especially in conifers. Here, we quantified the changing cellular morphology patterns of cambial zones during the active, transition and dormant stages. With the help of toluidine blue and periodic acid-Schiff staining to visualize cell walls and identify their constituents, we observed decreasing cambial cell layers, thickening of newly formed xylem cell walls and increased polysaccharide granules in phloem from June to the following March over the course of our collecting period. Pectin immunofluorescence showed that dormant-stage cambium can produce highly abundant de-esterified homogalacturonan and (1-4)-β-d-galactan epitopes, whereas active cambium can strong accumulate high methylesterified homogalacturonan. Calcofluor white staining and confocal Raman spectroscopy analysis revealed regular changes in the chemical composition of cell walls, such as relative lower cellulose deposition in transition stage in vascular cambium, and higher lignin accumulation was found in dormant stage in secondary xylem. Moreover, real-time quantitative polymerase chain reaction analysis suggested that various IAA (Aux/IAA protein), CesA, CslA and HDZ genes, as well as NAC, PME3 and PME4, may be involved in cambium activities and secondary xylem formation. Taken together, these findings provide new information about cambium activity and cell differentiation in the formation, structure and chemistry in conifers during the active-dormant transition.

Transcriptomic and epigenomic remodeling occurs during vascular cambium periodicity in Populus tomentosa

Trees in temperate regions exhibit evident seasonal patterns, which play vital roles in their growth and development. The activity of cambial stem cells is the basis for regulating the quantity and quality of wood, which has received considerable attention. However, the underlying mechanisms of these processes have not been fully elucidated. Here we performed a comprehensive analysis of morphological observations, transcriptome profiles, the DNA methylome, and miRNAs of the cambium in Populus tomentosa during the transition from dormancy to activation. Anatomical analysis showed that the active cambial zone exhibited a significant increase in the width and number of cell layers compared with those of the dormant and reactivating cambium. Furthermore, we found that differentially expressed genes associated with vascular development were mainly involved in plant hormone signal transduction, cell division and expansion, and cell wall biosynthesis. In addition, we identified 235 known miRNAs and 125 novel miRNAs. Differentially expressed miRNAs and target genes showed stronger negative correlations than other miRNA/target pairs. Moreover, global methylation and transcription analysis revealed that CG gene body methylation was positively correlated with gene expression, whereas CHG exhibited the opposite trend in the downstream region. Most importantly, we observed that the number of CHH differentially methylated region (DMR) changes was the greatest during cambium periodicity. Intriguingly, the genes with hypomethylated CHH DMRs in the promoter were involved in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interactions during vascular cambium development. These findings improve our systems-level understanding of the epigenomic diversity that exists in the annual growth cycle of trees.

Exploring the loblolly pine (Pinus taeda L.) genome by BAC sequencing and Cot analysis

Loblolly pine (LP; Pinus taeda L.) is an economically and ecologically important tree in the southeastern U.S. To advance understanding of the loblolly pine (LP; Pinus taeda L.) genome, we sequenced and analyzed 100 BAC clones and performed a Cot analysis. The Cot analysis indicates that the genome is composed of 57, 24, and 10% highly-repetitive, moderately-repetitive, and single/low-copy sequences, respectively (the remaining 9% of the genome is a combination of fold back and damaged DNA). Although single/low-copy DNA only accounts for 10% of the LP genome, the amount of single/low-copy DNA in LP is still 14 times the size of the Arabidopsis genome. Since gene numbers in LP are similar to those in Arabidopsis, much of the single/low-copy DNA of LP would appear to be composed of DNA that is both gene- and repeat-poor. Macroarrays prepared from a LP bacterial artificial chromosome (BAC) library were hybridized with probes designed from cell wall synthesis/wood development cDNAs, and 50 of the "targeted" clones were selected for further analysis. An additional 25 clones were selected because they contained few repeats, while 25 more clones were selected at random. The 100 BAC clones were Sanger sequenced and assembled. Of the targeted BACs, 80% contained all or part of the cDNA used to target them. One targeted BAC was found to contain fungal DNA and was eliminated from further analysis. Combinations of similarity-based and ab initio gene prediction approaches were utilized to identify and characterize potential coding regions in the 99 BACs containing LP DNA. From this analysis, we identified 154 gene models (GMs) representing both putative protein-coding genes and likely pseudogenes. Ten of the GMs (all of which were specifically targeted) had enough support to be classified as intact genes. Interestingly, the 154 GMs had statistically indistinguishable (α = 0.05) distributions in the targeted and random BAC clones (15.18 and 12.61 GM/Mb, respectively), whereas the low-repeat BACs contained significantly fewer GMs (7.08 GM/Mb). However, when GM length was considered, the targeted BACs had a significantly greater percentage of their length in GMs (3.26%) when compared to random (1.63%) and low-repeat (0.62%) BACs. The results of our study provide insight into LP evolution and inform ongoing efforts to produce a reference genome sequence for LP, while characterization of genes involved in cell wall production highlights carbon metabolism pathways that can be leveraged for increasing wood production.

Transcriptional Roadmap to Seasonal Variation in Wood Formation of Norway Spruce

Seasonal cues influence several aspects of the secondary growth of tree stems, including cambial activity, wood chemistry, and transition to latewood formation. We investigated seasonal changes in cambial activity, secondary cell wall formation, and tracheid cell death in woody tissues of Norway spruce (Picea abies) throughout one seasonal cycle. RNA sequencing was performed simultaneously in both the xylem and cambium/phloem tissues of the stem. Principal component analysis revealed gradual shifts in the transcriptomes that followed a chronological order throughout the season. A notable remodeling of the transcriptome was observed in the winter, with many genes having maximal expression during the coldest months of the year. A highly coexpressed set of monolignol biosynthesis genes showed high expression during the period of secondary cell wall formation as well as a second peak in midwinter. This midwinter peak in expression did not trigger lignin deposition, as determined by pyrolysis-gas chromatography/mass spectrometry. Coexpression consensus network analyses suggested the involvement of transcription factors belonging to the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES and MYELOBLASTOSIS-HOMEOBOX families in the seasonal control of secondary cell wall formation of tracheids. Interestingly, the lifetime of the latewood tracheids stretched beyond the winter dormancy period, correlating with a lack of cell death-related gene expression. Our transcriptomic analyses combined with phylogenetic and microscopic analyses also identified the cellulose and lignin biosynthetic genes and putative regulators for latewood formation and tracheid cell death in Norway spruce, providing a toolbox for further physiological and functional assays of these important phase transitions.

Why we need more non-seed plant models

Contents 355 I. 355 II. 356 III. 356 IV. 357 V. 358 VI. 359 359 References 359 SUMMARY: Out of a hundred sequenced and published land plant genomes, four are not of flowering plants. This severely skewed taxonomic sampling hinders our comprehension of land plant evolution at large. Moreover, most genetically accessible model species are flowering plants as well. If we are to gain a deeper understanding of how plants evolved and still evolve, and which of their developmental patterns are ancestral or derived, we need to study a more diverse set of plants. Here, I thus argue that we need to sequence genomes of so far neglected lineages, and that we need to develop more non-seed plant model species.

Cell Wall Ultrastructure of Stem Wood, Roots, and Needles of a Conifer Varies in Response to Moisture Availability

The composition, integrity, and architecture of the macromolecular matrix of cell walls, collectively referred to as cell wall ultrastructure, exhibits variation across species and organs and among cell types within organs. Indirect approaches have suggested that modifications to cell wall ultrastructure occur in response to abiotic stress; however, modifications have not been directly observed. Glycome profiling was used to study cell wall ultrastructure by examining variation in composition and extractability of non-cellulosic glycans in cell walls of stem wood, roots, and needles of loblolly pine saplings exposed to high and low soil moisture. Soil moisture influenced physiological processes and the overall composition and extractability of cell wall components differed as a function of soil moisture treatments. The strongest response of cell wall ultrastructure to soil moisture was increased extractability of pectic backbone epitopes in the low soil moisture treatment. The higher abundance of these pectic backbone epitopes in the oxalate extract indicate that the loosening of cell wall pectic components could be associated with the release of pectic signals as a stress response. The increased extractability of pectic backbone epitopes in response to low soil moisture availability was more pronounced in stem wood than in roots or needles. Additional responses to low soil moisture availability were observed in lignin-associated carbohydrates released in chlorite extracts of stem wood, including an increased abundance of pectic arabinogalactan epitopes. Overall, these results indicate that cell walls of loblolly pine organs undergo changes in their ultrastructural composition and extractability as a response to soil moisture availability and that cell walls of the stem wood are more responsive to low soil moisture availability compared to cell walls of roots and needles. To our knowledge, this is the first direct evidence, delineated by glycomic analyses, that abiotic stress affects cell wall ultrastructure. This study is also unique in that glycome profiling of pine needles has never before been reported.

Conifer genomics and adaptation: at the crossroads of genetic diversity and genome function

Conifers have been understudied at the genomic level despite their worldwide ecological and economic importance but the situation is rapidly changing with the development of next generation sequencing (NGS) technologies. With NGS, genomics research has simultaneously gained in speed, magnitude and scope. In just a few years, genomes of 20-24 gigabases have been sequenced for several conifers, with several others expected in the near future. Biological insights have resulted from recent sequencing initiatives as well as genetic mapping, gene expression profiling and gene discovery research over nearly two decades. We review the knowledge arising from conifer genomics research emphasizing genome evolution and the genomic basis of adaptation, and outline emerging questions and knowledge gaps. We discuss future directions in three areas with potential inputs from NGS technologies: the evolutionary impacts of adaptation in conifers based on the adaptation-by-speciation model; the contributions of genetic variability of gene expression in adaptation; and the development of a broader understanding of genetic diversity and its impacts on genome function. These research directions promise to sustain research aimed at addressing the emerging challenges of adaptation that face conifer trees.

Two iron-regulated transporter (IRT) genes showed differential expression in poplar trees under iron or zinc deficiency

Two iron-regulated transporter (IRT) genes were cloned from the iron chlorosis resistant (PtG) and susceptible (PtY) Populus tremula 'Erecta' lines. Nucleotide sequence analysis showed no significant difference between PtG and PtY. The predicted proteins contain a conserved ZIP domain with 8 transmembrane (TM) regions. A ZIP signature sequence was found in the fourth TM domain. Phylogenetic analysis revealed that PtIRT1 was clustered with tomato and tobacco IRT genes that are highly responsible to iron deficiency. The PtIRT3 gene was clustered with the AtIRT3 gene that was related to zinc and iron transport in plants. Tissue specific expression indicated that PtIRT1 only expressed in the root, while PtIRT3 constitutively expressed in all tested tissues. Under iron deficiency, the expression of PtIRT1 was dramatically increased and a significantly higher transcript level was detected in PtG than in PtY. Iron deficiency also enhanced the expression of PtIRT3 in PtG. On the other hand, zinc deficiency down-regulated the expression of PtIRT1 and PtIRT3 in both PtG and PtY. Zinc accumulated significantly under iron-deficient conditions, whereas the zinc deficiency showed no significant effect on iron accumulation. A yeast complementation test revealed that the PtIRT1 and PtIRT3 genes could restore the iron uptake ability under the iron uptake-deficiency condition. The results will help understand the mechanisms of iron deficiency response in poplar trees and other woody species.

Comparative interrogation of the developing xylem transcriptomes of two wood-forming species: Populus trichocarpa and Eucalyptus grandis

Wood formation is a complex developmental process governed by genetic and environmental stimuli. Populus and Eucalyptus are fast-growing, high-yielding tree genera that represent ecologically and economically important species suitable for generating significant lignocellulosic biomass. Comparative analysis of the developing xylem and leaf transcriptomes of Populus trichocarpa and Eucalyptus grandis together with phylogenetic analyses identified clusters of homologous genes preferentially expressed during xylem formation in both species. A conserved set of 336 single gene pairs showed highly similar xylem preferential expression patterns, as well as evidence of high functional constraint. Individual members of multi-gene orthologous clusters known to be involved in secondary cell wall biosynthesis also showed conserved xylem expression profiles. However, species-specific expression as well as opposite (xylem versus leaf) expression patterns observed for a subset of genes suggest subtle differences in the transcriptional regulation important for xylem development in each species. Using sequence similarity and gene expression status, we identified functional homologs likely to be involved in xylem developmental and biosynthetic processes in Populus and Eucalyptus. Our study suggests that, while genes involved in secondary cell wall biosynthesis show high levels of gene expression conservation, differential regulation of some xylem development genes may give rise to unique xylem properties.

Development of genetic markers in Eucalyptus species by target enrichment and exome sequencing

The advent of next-generation sequencing has facilitated large-scale discovery, validation and assessment of genetic markers for high density genotyping. The present study was undertaken to identify markers in genes supposedly related to wood property traits in three Eucalyptus species. Ninety four genes involved in xylogenesis were selected for hybridization probe based nuclear genomic DNA target enrichment and exome sequencing. Genomic DNA was isolated from the leaf tissues and used for on-array probe hybridization followed by Illumina sequencing. The raw sequence reads were trimmed and high-quality reads were mapped to the E. grandis reference sequence and the presence of single nucleotide variants (SNVs) and insertions/ deletions (InDels) were identified across the three species. The average read coverage was 216X and a total of 2294 SNVs and 479 InDels were discovered in E. camaldulensis, 2383 SNVs and 518 InDels in E. tereticornis, and 1228 SNVs and 409 InDels in E. grandis. Additionally, SNV calling and InDel detection were conducted in pair-wise comparisons of E. tereticornis vs. E. grandis, E. camaldulensis vs. E. tereticornis and E. camaldulensis vs. E. grandis. This study presents an efficient and high throughput method on development of genetic markers for family– based QTL and association analysis in Eucalyptus.

Unique features of the loblolly pine (Pinus taeda L.) megagenome revealed through sequence annotation

The largest genus in the conifer family Pinaceae is Pinus, with over 100 species. The size and complexity of their genomes (∼20-40 Gb, 2n = 24) have delayed the arrival of a well-annotated reference sequence. In this study, we present the annotation of the first whole-genome shotgun assembly of loblolly pine (Pinus taeda L.), which comprises 20.1 Gb of sequence. The MAKER-P annotation pipeline combined evidence-based alignments and ab initio predictions to generate 50,172 gene models, of which 15,653 are classified as high confidence. Clustering these gene models with 13 other plant species resulted in 20,646 gene families, of which 1554 are predicted to be unique to conifers. Among the conifer gene families, 159 are composed exclusively of loblolly pine members. The gene models for loblolly pine have the highest median and mean intron lengths of 24 fully sequenced plant genomes. Conifer genomes are full of repetitive DNA, with the most significant contributions from long-terminal-repeat retrotransposons. In depth analysis of the tandem and interspersed repetitive content yielded a combined estimate of 82%.

Generation, functional annotation and comparative analysis of black spruce (Picea mariana) ESTs: an important conifer genomic resource

Background

EST (expressed sequence tag) sequences and their annotation provide a highly valuable resource for gene discovery, genome sequence annotation, and other genomics studies that can be applied in genetics, breeding and conservation programs for non-model organisms. Conifers are long-lived plants that are ecologically and economically important globally, and have a large genome size. Black spruce (Picea mariana), is a transcontinental species of the North American boreal and temperate forests. However, there are limited transcriptomic and genomic resources for this species. The primary objective of our study was to develop a black spruce transcriptomic resource to facilitate on-going functional genomics projects related to growth and adaptation to climate change.

Results

We conducted bidirectional sequencing of cDNA clones from a standard cDNA library constructed from black spruce needle tissues. We obtained 4,594 high quality (2,455 5' end and 2,139 3' end) sequence reads, with an average read-length of 532 bp. Clustering and assembly of ESTs resulted in 2,731 unique sequences, consisting of 2,234 singletons and 497 contigs. Approximately two-thirds (63%) of unique sequences were functionally annotated. Genes involved in 36 molecular functions and 90 biological processes were discovered, including 24 putative transcription factors and 232 genes involved in photosynthesis. Most abundantly expressed transcripts were associated with photosynthesis, growth factors, stress and disease response, and transcription factors. A total of 216 full-length genes were identified. About 18% (493) of the transcripts were novel, representing an important addition to the Genbank EST database (dbEST). Fifty-seven di-, tri-, tetra- and penta-nucleotide simple sequence repeats were identified.

Conclusions

We have developed the first high quality EST resource for black spruce and identified 493 novel transcripts, which may be species-specific related to life history and ecological traits. We have also identified full-length genes and microsatellite-containing ESTs. Based on EST sequence similarities, black spruce showed close evolutionary relationships with congeneric Picea glauca and Picea sitchensis compared to other Pinaceae members and angiosperms. The EST sequences reported here provide an important resource for genome annotation, functional and comparative genomics, molecular breeding, conservation and management studies and applications in black spruce and related conifer species.

The Arabidopsis wood model-the case for the inflorescence stem

Arabidopsis thaliana has successfully served as a model to discover genes and proteins that have roles in a wide range of plant traits, including wood-related traits, such as lignin, cellulose and hemicellulose biosynthesis, secondary growth regulation, and secondary cell wall synthesis. Both the radially thickened hypocotyl and the inflorescence stem (flower stalk) have been studied. In this review, we address lingering doubts regarding the utility of Arabidopsis as a model for wood development by highlighting studies that provide new biochemical and biophysical evidence that extend support for the Arabidopsis inflorescence stem as a model for wood development beyond what is currently thought. We describe different aspects of Arabidopsis that make it a highly versatile tool for the study of wood development. One would likely utilise the radially thickened hypocotyl because of its more fully developed vascular cambium for traits related specifically to secondary (i.e. cambial) growth. It is more productive to utilise the inflorescence stem for wood-like biophysical traits. Accession variation has been underexploited as a powerful method to discover genes governing wood-like traits. We discuss recent findings that survey the accession variation in Arabidopsis for biochemical and biophysical properties of various wood traits, such as microfibril angle, tensile strength and cellulose/hemicellulose content. Furthermore we discuss how larger-scale studies of this nature using plants grown in long days (as opposed to the current short-day paradigm) could accelerate gene discovery and our understanding of cell wall and wood development. We highlight some relatively unexplored areas of research relating to the secondary cell wall composition, architecture and biophysical properties of the inflorescence stem, and how these traits are relevant to wood formation. The Arabidopsis inflorescence stem has other characteristics, expressed genes and traits held in common with woody species that have not been widely characterised or discussed to date. We discuss how this conservation may indicate the more general potential for "true" woodiness in herbaceous species, in the context of so-called secondary woodiness.

The regulation of cambial activity in Chinese fir (Cunninghamia lanceolata) involves extensive transcriptome remodeling

Chinese fir (Cunninghamia lanceolata), a commercially important tree for the timber and pulp industry, is widely distributed in southern China and northern Vietnam, but its large and complex genome has hindered the development of genomic resources. Few efforts have focused on analysis of the modulation of transcriptional networks in vascular cambium during the transition from active growth to dormancy in conifers. Here, we used Illumina sequencing to analyze the global transcriptome alterations at the different stages of vascular cambium development in Chinese fir. By analyzing dynamic changes in the transcriptome of vascular cambium based on our RNA sequencing (RNA-Seq) data at the dormant, reactivating and active stages, many potentially interesting genes were identified that encoded putative regulators of cambial activity, cell division, cell expansion and cell wall biosynthesis and modification. In particular, the genes involved in transcriptional regulation and hormone signaling were highlighted to reveal their biological importance in the cambium development and wood formation. Our results reveal the dynamics of transcriptional networks and identify potential key components in the regulation of vascular cambium development in Chinese fir, which will contribute to the in-depth study of cambial differentiation and wood-forming candidate genes in conifers.

Whole-exome targeted sequencing of the uncharacterized pine genome

The large genome size of many species hinders the development and application of genomic tools to study them. For instance, loblolly pine (Pinus taeda L.), an ecologically and economically important conifer, has a large and yet uncharacterized genome of 21.7 Gbp. To characterize the pine genome, we performed exome capture and sequencing of 14 729 genes derived from an assembly of expressed sequence tags. Efficiency of sequence capture was evaluated and shown to be similar across samples with increasing levels of complexity, including haploid cDNA, haploid genomic DNA and diploid genomic DNA. However, this efficiency was severely reduced for probes that overlapped multiple exons, presumably because intron sequences hindered probe:exon hybridizations. Such regions could not be entirely avoided during probe design, because of the lack of a reference sequence. To improve the throughput and reduce the cost of sequence capture, a method to multiplex the analysis of up to eight samples was developed. Sequence data showed that multiplexed capture was reproducible among 24 haploid samples, and can be applied for high-throughput analysis of targeted genes in large populations. Captured sequences were de novo assembled, resulting in 11 396 expanded and annotated gene models, significantly improving the knowledge about the pine gene space. Interspecific capture was also evaluated with over 98% of all probes designed from P. taeda that were efficient in sequence capture, were also suitable for analysis of the related species Pinus elliottii Engelm.

Open access to tree genomes: the path to a better forest

An open-access culture and a well-developed comparative-genomics infrastructure must be developed in forest trees to derive the full potential of genome sequencing in this diverse group of plants that are the dominant species in much of the earth's terrestrial ecosystems.

The genomic architecture and association genetics of adaptive characters using a candidate SNP approach in boreal black spruce

BACKGROUND

The genomic architecture of adaptive traits remains poorly understood in non-model plants. Various approaches can be used to bridge this gap, including the mapping of quantitative trait loci (QTL) in pedigrees, and genetic association studies in non-structured populations. Here we present results on the genomic architecture of adaptive traits in black spruce, which is a widely distributed conifer of the North American boreal forest. As an alternative to the usual candidate gene approach, a candidate SNP approach was developed for association testing.

RESULTS

A genetic map containing 231 gene loci was used to identify QTL that were related to budset timing and to tree height assessed over multiple years and sites. Twenty-two unique genomic regions were identified, including 20 that were related to budset timing and 6 that were related to tree height. From results of outlier detection and bulk segregant analysis for adaptive traits using DNA pool sequencing of 434 genes, 52 candidate SNPs were identified and subsequently tested in genetic association studies for budset timing and tree height assessed over multiple years and sites. A total of 34 (65%) SNPs were significantly associated with budset timing, or tree height, or both. Although the percentages of explained variance (PVE) by individual SNPs were small, several significant SNPs were shared between sites and among years.

CONCLUSIONS

The sharing of genomic regions and significant SNPs between budset timing and tree height indicates pleiotropic effects. Significant QTLs and SNPs differed quite greatly among years, suggesting that different sets of genes for the same characters are involved at different stages in the tree's life history. The functional diversity of genes carrying significant SNPs and low observed PVE further indicated that a large number of polymorphisms are involved in adaptive genetic variation. Accordingly, for undomesticated species such as black spruce with natural populations of large effective size and low linkage disequilibrium, efficient marker systems that are predictive of adaptation should require the survey of large numbers of SNPs. Candidate SNP approaches like the one developed in the present study could contribute to reducing these numbers.

Transcriptome characterisation of Pinus tabuliformis and evolution of genes in the Pinus phylogeny

BACKGROUND

The Chinese pine (Pinus tabuliformis) is an indigenous conifer species in northern China but is relatively underdeveloped as a genomic resource; thus, limiting gene discovery and breeding. Large-scale transcriptome data were obtained using a next-generation sequencing platform to compensate for the lack of P. tabuliformis genomic information.

RESULTS

The increasing amount of transcriptome data on Pinus provides an excellent resource for multi-gene phylogenetic analysis and studies on how conserved genes and functions are maintained in the face of species divergence. The first P. tabuliformis transcriptome from a normalised cDNA library of multiple tissues and individuals was sequenced in a full 454 GS-FLX run, producing 911,302 sequencing reads. The high quality overlapping expressed sequence tags (ESTs) were assembled into 46,584 putative transcripts, and more than 700 SSRs and 92,000 SNPs/InDels were characterised. Comparative analysis of the transcriptome of six conifer species yielded 191 orthologues, from which we inferred a phylogenetic tree, evolutionary patterns and calculated rates of gene diversion. We also identified 938 fast evolving sequences that may be useful for identifying genes that perhaps evolved in response to positive selection and might be responsible for speciation in the Pinus lineage.

CONCLUSIONS

A large collection of high-quality ESTs was obtained, de novo assembled and characterised, which represents a dramatic expansion of the current transcript catalogues of P. tabuliformis and which will gradually be applied in breeding programs of P. tabuliformis. Furthermore, these data will facilitate future studies of the comparative genomics of P. tabuliformis and other related species.

Sequencing of the needle transcriptome from Norway spruce (Picea abies Karst L.) reveals lower substitution rates, but similar selective constraints in gymnosperms and angiosperms

BACKGROUND

A detailed knowledge about spatial and temporal gene expression is important for understanding both the function of genes and their evolution. For the vast majority of species, transcriptomes are still largely uncharacterized and even in those where substantial information is available it is often in the form of partially sequenced transcriptomes. With the development of next generation sequencing, a single experiment can now simultaneously identify the transcribed part of a species genome and estimate levels of gene expression.

RESULTS

mRNA from actively growing needles of Norway spruce (Picea abies) was sequenced using next generation sequencing technology. In total, close to 70 million fragments with a length of 76 bp were sequenced resulting in 5 Gbp of raw data. A de novo assembly of these reads, together with publicly available expressed sequence tag (EST) data from Norway spruce, was used to create a reference transcriptome. Of the 38,419 PUTs (putative unique transcripts) longer than 150 bp in this reference assembly, 83.5% show similarity to ESTs from other spruce species and of the remaining PUTs, 3,704 show similarity to protein sequences from other plant species, leaving 4,167 PUTs with limited similarity to currently available plant proteins. By predicting coding frames and comparing not only the Norway spruce PUTs, but also PUTs from the close relatives Picea glauca and Picea sitchensis to both Pinus taeda and Taxus mairei, we obtained estimates of synonymous and non-synonymous divergence among conifer species. In addition, we detected close to 15,000 SNPs of high quality and estimated gene expression differences between samples collected under dark and light conditions.

CONCLUSIONS

Our study yielded a large number of single nucleotide polymorphisms as well as estimates of gene expression on transcriptome scale. In agreement with a recent study we find that the synonymous substitution rate per year (0.6 × 10-09 and 1.1 × 10-09) is an order of magnitude smaller than values reported for angiosperm herbs. However, if one takes generation time into account, most of this difference disappears. The estimates of the dN/dS ratio (non-synonymous over synonymous divergence) reported here are in general much lower than 1 and only a few genes showed a ratio larger than 1.

Transcriptional analysis of differentially expressed genes in response to stem inclination in young seedlings of pine

The gravitropic response in trees is a widely studied phenomenon, however understanding of the molecular mechanism involved remains unclear. The purpose of this work was to identify differentially expressed genes in response to inclination using a comparative approach for two conifer species. Young seedlings were subjected to inclination and samples were collected at four different times points. First, suppression subtractive hybridisation (SSH) was used to identify differentially regulated genes in radiata pine (Pinus radiata D. Don). cDNA libraries were constructed from the upper and lower part of inclined stems in a time course experiment, ranging from 2.5 h to 1 month. From a total of 3092 sequences obtained, 2203 elements were assembled, displaying homology to a public database. A total of 942 unigene elements were identified using bioinformatic tools after redundancy analysis. Of these, 614 corresponded to known function genes and 328 to unknown function genes, including hypothetical proteins. Comparative analysis between radiata pine and maritime pine (Pinus pinaster Ait.) was performed to validate the differential expression of relevant candidate genes using qPCR. Selected genes were involved in several functional categories: hormone regulation, phenylpropanoid pathway and signal transduction. This comparative approach for the two conifer species helped determine the molecular gene pattern generated by inclination, providing a set of Pinus gene signatures that may be involved in the gravitropic stress response. These genes may also represent relevant candidate genes involved in the gravitropic response and potentially in wood formation.

MicroRNAs in trees

MicroRNAs (miRNAs) are 20-24 nucleotide long molecules processed from a specific class of RNA polymerase II transcripts that mainly regulate the stability of mRNAs containing a complementary sequence by targeted degradation in plants. Many features of tree biology are regulated by miRNAs affecting development, metabolism, adaptation and evolution. MiRNAs may be modified and harnessed for controlled suppression of specific genes to learn about gene function, or for practical applications through genetic engineering. Modified (artificial) miRNAs act as dominant suppressors and are particularly useful in tree genetics because they bypass the generations of inbreeding needed for fixation of recessive mutations. The purpose of this review is to summarize the current status of information on miRNAs in trees and to guide future studies on the role of miRNAs in the biology of woody perennials and to illustrate their utility in directed genetic modification of trees.

Transcriptome profiling in conifers and the PiceaGenExpress database show patterns of diversification within gene families and interspecific conservation in vascular gene expression

Background

Conifers have very large genomes (13 to 30 Gigabases) that are mostly uncharacterized although extensive cDNA resources have recently become available. This report presents a global overview of transcriptome variation in a conifer tree and documents conservation and diversity of gene expression patterns among major vegetative tissues.

Results

An oligonucleotide microarray was developed from Picea glauca and P. sitchensis cDNA datasets. It represents 23,853 unique genes and was shown to be suitable for transcriptome profiling in several species. A comparison of secondary xylem and phelloderm tissues showed that preferential expression in these vascular tissues was highly conserved among Picea spp. RNA-Sequencing strongly confirmed tissue preferential expression and provided a robust validation of the microarray design. A small database of transcription profiles called PiceaGenExpress was developed from over 150 hybridizations spanning eight major tissue types. In total, transcripts were detected for 92% of the genes on the microarray, in at least one tissue. Non-annotated genes were predominantly expressed at low levels in fewer tissues than genes of known or predicted function. Diversity of expression within gene families may be rapidly assessed from PiceaGenExpress. In conifer trees, dehydrins and late embryogenesis abundant (LEA) osmotic regulation proteins occur in large gene families compared to angiosperms. Strong contrasts and low diversity was observed in the dehydrin family, while diverse patterns suggested a greater degree of diversification among LEAs.

Conclusion

Together, the oligonucleotide microarray and the PiceaGenExpress database represent the first resource of this kind for gymnosperm plants. The spruce transcriptome analysis reported here is expected to accelerate genetic studies in the large and important group comprised of conifer trees.

Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis

MicroRNAs (miRNAs) are emerging as essential regulators of biological processes. Somatic embryogenesis is one of the most important techniques for gymnosperm-breeding programs, but there is little understanding of its underlying mechanism. To investigate the roles of miRNAs during somatic embryogenesis in larch, we constructed a small RNA library from somatic embryos. High-throughput sequencing of the library identified 83 conserved miRNAs from 35 families, 16 novel miRNAs, and 14 plausible miRNA candidates, with a high proportion specific to larch or gymnosperms. qRT-PCR analysis demonstrated that both the conserved and novel or candidate miRNAs were expressed in larch. Several miRNA precursor sequences were obtained via RACE. We predicted 110 target genes using bioinformatics, and validated 9 of them by 5' RACE. 11 conserved miRNA families including 17 miRNAs with critical functions in plant development and six target mRNAs were detected by qRT-PCR in the larch SE. Stage-specific expression of miRNAs and their targets indicate their possible modulation on SE of larch: miR171a/b might exert function on PEMs, while miR171c acts in the induction process of larch SE; miR397 and miR398 mainly involved in modulation of PEM propagation and transition to single embryo; miR162 and miR168 exert their regulatory function during total SE process, especially during stages 5-8; miR156, miR159, miR160, miR166, miR167, and miR390 might play regulatory roles during cotyledonary embryo development. These findings indicate that larch and possibly other gymnosperms have complex mechanisms of gene regulation involving specific and common miRNAs operating post-transcriptionally during embryogenesis.

Deciphering the enigma of lignification: precursor transport, oxidation, and the topochemistry of lignin assembly

Plant lignification is a tightly regulated complex cellular process that occurs via three sequential steps: the synthesis of monolignols within the cytosol; the transport of monomeric precursors across plasma membrane; and the oxidative polymerization of monolignols to form lignin macromolecules within the cell wall. Although we have a reasonable understanding of monolignol biosynthesis, many aspects of lignin assembly remain elusive. These include the precursors' transport and oxidation, and the initiation of lignin polymerization. This review describes our current knowledge of the molecular mechanisms underlying monolignol transport and oxidation, discusses the intriguing yet least-understood aspects of lignin assembly, and highlights the technologies potentially aiding in clarifying the enigma of plant lignification.

A white spruce gene catalog for conifer genome analyses

Several angiosperm plant genomes, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), poplar (Populus trichocarpa), and grapevine (Vitis vinifera), have been sequenced, but the lack of reference genomes in gymnosperm phyla reduces our understanding of plant evolution and restricts the potential impacts of genomics research. A gene catalog was developed for the conifer tree Picea glauca (white spruce) through large-scale expressed sequence tag sequencing and full-length cDNA sequencing to facilitate genome characterizations, comparative genomics, and gene mapping. The resource incorporates new and publicly available sequences into 27,720 cDNA clusters, 23,589 of which are represented by full-length insert cDNAs. Expressed sequence tags, mate-pair cDNA clone analysis, and custom sequencing were integrated through an iterative process to improve the accuracy of clustering outcomes. The entire catalog spans 30 Mb of unique transcribed sequence. We estimated that the P. glauca nuclear genome contains up to 32,520 transcribed genes owing to incomplete, partially sequenced, and unsampled transcripts and that its transcriptome could span up to 47 Mb. These estimates are in the same range as the Arabidopsis and rice transcriptomes. Next-generation methods confirmed and enhanced the catalog by providing deeper coverage for rare transcripts, by extending many incomplete clusters, and by augmenting the overall transcriptome coverage to 38 Mb of unique sequence. Genomic sample sequencing at 8.5% of the 19.8-Gb P. glauca genome identified 1,495 clusters representing highly repeated sequences among the cDNA clusters. With a conifer transcriptome in full view, functional and protein domain annotations clearly highlighted the divergences between conifers and angiosperms, likely reflecting their respective evolutionary paths.

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.

Quantitative evolution of transposable and satellite DNA sequences in Picea species

Clones containing tandemly arranged repeats belonging to two distinct sequence families, (i) PAG004P22F (2F) and PAG004E03C (3C) or (ii) Ty3/gypsy- (8R; PAG004B08R) and Ty1/copia-like sequences (9R; PAG007F19R), were selected from a randomly sheared total genomic DNA library of Picea abies . The inserts were used as probes in dot-blot hybridizations to genomic DNA of P. abies, Picea orientalis , Picea pungens , and Picea pungens var. glauca. All these entities are diploid and share the same chromosome number (2n = 24), but the genome sizes differ largely. The redundancy (copy number per 1C DNA) of sequences related to each probe varied greatly between the genomes. No significant correlation was found between the genome size and the copy number of sequences in any family. The quantitative ratios varied greatly (in each genome) between the two families of satellite DNA, between the sequences that represented copia or gypsy retrotransposons, and between tandemly arranged sequences and retroelements as a whole, suggesting that there is no common factor that controls the quantitative evolution of repeats belonging to different sequence families during speciation in Picea.

Forest tree genomics: growing resources and applications

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment - enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

Sequestration and transport of lignin monomeric precursors

Lignin is the second most abundant terrestrial biopolymer after cellulose. It is essential for the viability of vascular plants. Lignin precursors, the monolignols, are synthesized within the cytosol of the cell. Thereafter, these monomeric precursors are exported into the cell wall, where they are polymerized and integrated into the wall matrix. Accordingly, transport of monolignols across cell membranes is a critical step affecting deposition of lignin in the secondarily thickened cell wall. While the biosynthesis of monolignols is relatively well understood, our knowledge of sequestration and transport of these monomers is sketchy. In this article, we review different hypotheses on monolignol transport and summarize the recent progresses toward the understanding of the molecular mechanisms underlying monolignol sequestration and transport across membranes. Deciphering molecular mechanisms for lignin precursor transport will support a better biotechnological solution to manipulate plant lignification for more efficient agricultural and industrial applications of cell wall biomass.

Subgroup 4 R2R3-MYBs in conifer trees: gene family expansion and contribution to the isoprenoid- and flavonoid-oriented responses

Transcription factors play a fundamental role in plants by orchestrating temporal and spatial gene expression in response to environmental stimuli. Several R2R3-MYB genes of the Arabidopsis subgroup 4 (Sg4) share a C-terminal EAR motif signature recently linked to stress response in angiosperm plants. It is reported here that nearly all Sg4 MYB genes in the conifer trees Picea glauca (white spruce) and Pinus taeda (loblolly pine) form a monophyletic clade (Sg4C) that expanded following the split of gymnosperm and angiosperm lineages. Deeper sequencing in P. glauca identified 10 distinct Sg4C sequences, indicating over-representation of Sg4 sequences compared with angiosperms such as Arabidopsis, Oryza, Vitis, and Populus. The Sg4C MYBs share the EAR motif core. Many of them had stress-responsive transcript profiles after wounding, jasmonic acid (JA) treatment, or exposure to cold in P. glauca and P. taeda, with MYB14 transcripts accumulating most strongly and rapidly. Functional characterization was initiated by expressing the P. taeda MYB14 (PtMYB14) gene in transgenic P. glauca plantlets with a tissue-preferential promoter (cinnamyl alcohol dehydrogenase) and a ubiquitous gene promoter (ubiquitin). Histological, metabolite, and transcript (microarray and targeted quantitative real-time PCR) analyses of PtMYB14 transgenics, coupled with mechanical wounding and JA application experiments on wild-type plantlets, allowed identification of PtMYB14 as a putative regulator of an isoprenoid-oriented response that leads to the accumulation of sesquiterpene in conifers. Data further suggested that PtMYB14 may contribute to a broad defence response implicating flavonoids. This study also addresses the potential involvement of closely related Sg4C sequences in stress responses and plant evolution.

Comparative genomics reveals conservative evolution of the xylem transcriptome in vascular plants

BACKGROUND

Wood is a valuable natural resource and a major carbon sink. Wood formation is an important developmental process in vascular plants which played a crucial role in plant evolution. Although genes involved in xylem formation have been investigated, the molecular mechanisms of xylem evolution are not well understood. We use comparative genomics to examine evolution of the xylem transcriptome to gain insights into xylem evolution.

RESULTS

The xylem transcriptome is highly conserved in conifers, but considerably divergent in angiosperms. The functional domains of genes in the xylem transcriptome are moderately to highly conserved in vascular plants, suggesting the existence of a common ancestral xylem transcriptome. Compared to the total transcriptome derived from a range of tissues, the xylem transcriptome is relatively conserved in vascular plants. Of the xylem transcriptome, cell wall genes, ancestral xylem genes, known proteins and transcription factors are relatively more conserved in vascular plants. A total of 527 putative xylem orthologs were identified, which are unevenly distributed across the Arabidopsis chromosomes with eight hot spots observed. Phylogenetic analysis revealed that evolution of the xylem transcriptome has paralleled plant evolution. We also identified 274 conifer-specific xylem unigenes, all of which are of unknown function. These xylem orthologs and conifer-specific unigenes are likely to have played a crucial role in xylem evolution.

CONCLUSIONS

Conifers have highly conserved xylem transcriptomes, while angiosperm xylem transcriptomes are relatively diversified. Vascular plants share a common ancestral xylem transcriptome. The xylem transcriptomes of vascular plants are more conserved than the total transcriptomes. Evolution of the xylem transcriptome has largely followed the trend of plant evolution.

Evolution of development of vascular cambia and secondary growth

Secondary growth from vascular cambia results in radial, woody growth of stems. The innovation of secondary vascular development during plant evolution allowed the production of novel plant forms ranging from massive forest trees to flexible, woody lianas. We present examples of the extensive phylogenetic variation in secondary vascular growth and discuss current knowledge of genes that regulate the development of vascular cambia and woody tissues. From these foundations, we propose strategies for genomics-based research in the evolution of development, which is a next logical step in the study of secondary growth.

A microarray-based method for the parallel analysis of genotypes and expression profiles of wood-forming tissues in Eucalyptus grandis

BACKGROUND

Fast-growing Eucalyptus grandis trees are one of the most efficient producers of wood in South Africa. The most serious problem affecting the quality and yield of solid wood products is the occurrence of end splitting in logs. Selection of E. grandis planting stock that exhibit preferred wood qualities is thus a priority of the South African forestry industry. We used microarray-based DNA-amplified fragment length polymorphism (AFLP) analysis in combination with expression profiling to develop fingerprints and profile gene expression of wood-forming tissue of seven different E. grandis trees.

RESULTS

A 1578-probe cDNA microarray was constructed by arraying 768 cDNA-AFLP clones and 810 cDNA library clones from seven individual E. grandis trees onto silanised slides. The results revealed that 32% of the spotted fragments showed distinct expression patterns (with a fold change of at least 1.4 or -1.4 and a p value of 0.01) could be grouped into clusters representing co-expressed genes. Evaluation of the binary distribution of cDNA-AFLP fragments on the array showed that the individual genotypes could be discriminated.

CONCLUSION

A simple, yet general method was developed for genotyping and expression profiling of wood-forming tissue of E. grandis trees differing in their splitting characteristics and in their lignin contents. Evaluation of gene expression profiles and the binary distribution of cDNA-AFLP fragments on the chip suggest that the prototype chip developed could be useful for transcript profiling and for the identification of Eucalyptus trees with preferred wood quality traits in commercial breeding programmes.

Sequence analysis and functional characterization of the promoter of the Picea glauca Cinnamyl Alcohol Dehydrogenase gene in transgenic white spruce plants

The enzyme Cinnamyl Alcohol Dehydrogenase (CAD) catalyses the last step of lignin monomer synthesis, and is considered as a molecular marker of cell wall lignification in different plants species. Here, we report the isolation and analysis of 5' flanking genomic DNA regions upstream to the CAD gene, from two conifers, i.e. white spruce (Picea glauca (Moench) Voss) and loblolly pine (Pinus taeda L.). Sequence comparisons with available CAD gene promoters from angiosperms highlighted the conservation of cis-elements matching MYB, WRKY and bHLH binding sites. Functional characterization of the P. glauca CAD promoter used P. glauca seedlings stably transformed with a DNA fragment of 1,163 base pairs (PgCAD) fused to the beta-glucuronidase (GUS) gene. Histochemical observations of different vegetative organs of the transgenic trees showed that this sequence was sufficient to drive GUS expression in lignifying tissues, and more specifically in differentiating xylem cells. Quantitative RT-PCR experiments also indicated that the native CAD gene was preferentially expressed in differentiating xylem both in stems and roots. In addition, GUS expression driven by the PgCAD promoter was wound-inducible which was consistent with the accumulation of CAD mRNA in response to jasmonate application and mechanical wounding. The spruce CAD promoter represents a valuable tool for research and biotechnology applications related to xylem and wood.

Generation and analysis of expressed sequence tags from six developing xylem libraries in Pinus radiata D. Don

BACKGROUND

Wood is a major renewable natural resource for the timber, fibre and bioenergy industry. Pinus radiata D. Don is the most important commercial plantation tree species in Australia and several other countries; however, genomic resources for this species are very limited in public databases. Our primary objective was to sequence a large number of expressed sequence tags (ESTs) from genes involved in wood formation in radiata pine.

RESULTS

Six developing xylem cDNA libraries were constructed from earlywood and latewood tissues sampled at juvenile (7 yrs), transition (11 yrs) and mature (30 yrs) ages, respectively. These xylem tissues represent six typical development stages in a rotation period of radiata pine. A total of 6,389 high quality ESTs were collected from 5,952 cDNA clones. Assembly of 5,952 ESTs from 5' end sequences generated 3,304 unigenes including 952 contigs and 2,352 singletons. About 97.0% of the 5,952 ESTs and 96.1% of the unigenes have matches in the UniProt and TIGR databases. Of the 3,174 unigenes with matches, 42.9% were not assigned GO (Gene Ontology) terms and their functions are unknown or unclassified. More than half (52.1%) of the 5,952 ESTs have matches in the Pfam database and represent 772 known protein families. About 18.0% of the 5,952 ESTs matched cell wall related genes in the MAIZEWALL database, representing all 18 categories, 91 of all 174 families and possibly 557 genes. Fifteen cell wall-related genes are ranked in the 30 most abundant genes, including CesA, tubulin, AGP, SAMS, actin, laccase, CCoAMT, MetE, phytocyanin, pectate lyase, cellulase, SuSy, expansin, chitinase and UDP-glucose dehydrogenase. Based on the PlantTFDB database 41 of the 64 transcription factor families in the poplar genome were identified as being involved in radiata pine wood formation. Comparative analysis of GO term abundance revealed a distinct transcriptome in juvenile earlywood formation compared to other stages of wood development.

CONCLUSION

The first large scale genomic resource in radiata pine was generated from six developing xylem cDNA libraries. Cell wall-related genes and transcription factors were identified. Juvenile earlywood has a distinct transcriptome, which is likely to contribute to the undesirable properties of juvenile wood in radiata pine. The publicly available resource of radiata pine will also be valuable for gene function studies and comparative genomics in forest trees.

A conifer genomics resource of 200,000 spruce (Picea spp.) ESTs and 6,464 high-quality, sequence-finished full-length cDNAs for Sitka spruce (Picea sitchensis)

Background

Members of the pine family (Pinaceae), especially species of spruce (Picea spp.) and pine (Pinus spp.), dominate many of the world's temperate and boreal forests. These conifer forests are of critical importance for global ecosystem stability and biodiversity. They also provide the majority of the world's wood and fiber supply and serve as a renewable resource for other industrial biomaterials. In contrast to angiosperms, functional and comparative genomics research on conifers, or other gymnosperms, is limited by the lack of a relevant reference genome sequence. Sequence-finished full-length (FL)cDNAs and large collections of expressed sequence tags (ESTs) are essential for gene discovery, functional genomics, and for future efforts of conifer genome annotation.

Results

As part of a conifer genomics program to characterize defense against insects and adaptation to local environments, and to discover genes for the production of biomaterials, we developed 20 standard, normalized or full-length enriched cDNA libraries from Sitka spruce (P. sitchensis), white spruce (P. glauca), and interior spruce (P. glauca-engelmannii complex). We sequenced and analyzed 206,875 3'- or 5'-end ESTs from these libraries, and developed a resource of 6,464 high-quality sequence-finished FLcDNAs from Sitka spruce. Clustering and assembly of 147,146 3'-end ESTs resulted in 19,941 contigs and 26,804 singletons, representing 46,745 putative unique transcripts (PUTs). The 6,464 FLcDNAs were all obtained from a single Sitka spruce genotype and represent 5,718 PUTs.

Conclusion

This paper provides detailed annotation and quality assessment of a large EST and FLcDNA resource for spruce. The 6,464 Sitka spruce FLcDNAs represent the third largest sequence-verified FLcDNA resource for any plant species, behind only rice (Oryza sativa) and Arabidopsis (Arabidopsis thaliana), and the only substantial FLcDNA resource for a gymnosperm. Our emphasis on capturing FLcDNAs and ESTs from cDNA libraries representing herbivore-, wound- or elicitor-treated induced spruce tissues, along with incorporating normalization to capture rare transcripts, resulted in a rich resource for functional genomics and proteomics studies. Sequence comparisons against five plant genomes and the non-redundant GenBank protein database revealed that a substantial number of spruce transcripts have no obvious similarity to known angiosperm gene sequences. Opportunities for future applications of the sequence and clone resources for comparative and functional genomics are discussed.

Characterization of expressed sequence tags from a full-length enriched cDNA library of Cryptomeria japonica male strobili

Background

Cryptomeria japonica D. Don is one of the most commercially important conifers in Japan. However, the allergic disease caused by its pollen is a severe public health problem in Japan. Since large-scale analysis of expressed sequence tags (ESTs) in the male strobili of C. japonica should help us to clarify the overall expression of genes during the process of pollen development, we constructed a full-length enriched cDNA library that was derived from male strobili at various developmental stages.

Results

We obtained 36,011 expressed sequence tags (ESTs) from either one or both ends of 19,437 clones derived from the cDNA library of C. japonica male strobili at various developmental stages. The 19,437 cDNA clones corresponded to 10,463 transcripts. Approximately 80% of the transcripts resembled ESTs from Pinus and Picea, while approximately 75% had homologs in Arabidopsis. An analysis of homologies between ESTs from C. japonica male strobili and known pollen allergens in the Allergome Database revealed that products of 180 transcripts exhibited significant homology. Approximately 2% of the transcripts appeared to encode transcription factors. We identified twelve genes for MADS-box proteins among these transcription factors. The twelve MADS-box genes were classified as DEF/GLO/GGM13-, AG-, AGL6-, TM3- and TM8-like MIKC^C genes and type I MADS-box genes.

Conclusion

Our full-length enriched cDNA library derived from C. japonica male strobili provides information on expression of genes during the development of male reproductive organs. We provided potential allergens in C. japonica. We also provided new information about transcription factors including MADS-box genes expressed in male strobili of C. japonica. Large-scale gene discovery using full-length cDNAs is a valuable tool for studies of gymnosperm species.

Tracking monolignols during wood development in lodgepole pine

Secondary xylem (wood) formation in gymnosperms requires that the tracheid protoplasts first build an elaborate secondary cell wall from an array of polysaccharides and then reinforce it with lignin, an amorphous, three-dimensional product of the random radical coupling of monolignols. The objective of this study was to track the spatial distribution of monolignols during development as they move from symplasm to apoplasm. This was done by feeding [(3)H]phenylalanine ([(3)H]Phe) to dissected cambium/developing wood from lodgepole pine (Pinus contorta var latifolia) seedlings, allowing uptake and metabolism, then rapidly freezing the cells and performing autoradiography to detect the locations of the monolignols responsible for lignification. Parallel experiments showed that radioactivity was incorporated into polymeric lignin and a methanol-soluble pool that was characterized by high-performance liquid chromatography. [(3)H]Phe was incorporated into expected lignin precursors, such as coniferyl alcohol and p-coumaryl alcohol, as well as pinoresinol. Coniferin, the glucoside of coniferyl alcohol, was detected by high-performance liquid chromatography but was not radioactively labeled. With light microscopy, radiolabeled phenylpropanoids were detected in the rays as well as the tracheids, with the two cell types showing differential sensitivity to inhibitors of protein translation and phenylpropanoid metabolism. Secondary cell walls of developing tracheids were heavily labeled when incubated with [(3)H]Phe. Inside the cell, cytoplasm was most strongly labeled followed by Golgi and low-vacuole label. Inhibitor studies suggest that the Golgi signal could be attributed to protein, rather than phenylpropanoid, origins. These data, produced with the best microscopy tools that are available today, support a model in which unknown membrane transporters, rather than Golgi vesicles, export monolignols.

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce (Picea glauca) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations (F(ST)) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates (H(E) = 0.270 and F(ST) = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F(ST) at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.

Identification and analysis of differentially expressed genes in differentiating xylem of Chinese fir (Cunninghamia lanceolata) by suppression subtractive hybridization

Wood is an important raw material for global industries with rapidly increasing demand. To isolate the genes differentially expressed during xylogenesis of Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), we used a novel system. Forward and reverse subtracted cDNA libraries were constructed using the suppression subtractive hybridization method; for the forward library we used cDNA from the mutant Dugansha as the tester and cDNA from the wild-type clone Jurong 0 as the driver, and for the reverse library we used Jurong 0 cDNA as the tester and Dugansha cDNA as the driver. Transcriptional profiling was performed using a macroarray with 4 digoxigenin-labeled probes. We obtained 618 and 409 clones from the forward and the reverse subtracted library, respectively. A total of 405 unique expressed sequence tags (ESTs) were obtained. Forty percent of the ESTs exhibited homologies with proteins of known function and fell into 4 major classes: metabolism, cell wall biogenesis and remodeling, signal transduction, and stress. Real-time PCR was performed to confirm the results. The expression levels of 11 selected ESTs were consistent with both macroarray and real-time PCR results. The systematic analysis of genes involved in wood formation in Chinese fir provides valuable insights into the molecular mechanisms involved in xylem differentiation and is an important resource for forest research that can be directed toward understanding the genetic control of wood formation and future endeavors to modify wood and fiber properties for industrial use.