Transcriptional Roadmap to Seasonal Variation in Wood Formation of Norway Spruce

Leaf gene expression trajectories during the growing season are consistent between sites and years in American beech

Transcriptomics provides a versatile tool for ecological monitoring. Here, through genome-guided profiling of transcripts mapping to 33 042 gene models, expression differences can be discerned among multi-year and seasonal leaf samples collected from American beech trees at two latitudinally separated sites. Despite a bottleneck due to post-Columbian deforestation, the single nucleotide polymorphism-based population genetic background analysis has yielded sufficient variation to account for differences between populations and among individuals. Our expression analyses during spring-summer and summer-autumn transitions for two consecutive years involved 4197 differentially expressed protein coding genes. Using Populus orthologues we reconstructed a protein-protein interactome representing leaf physiological states of trees during the seasonal transitions. Gene set enrichment analysis revealed gene ontology terms that highlight molecular functions and biological processes possibly influenced by abiotic forcings such as recovery from drought and response to excess precipitation. Further, based on 324 co-regulated transcripts, we focused on a subset of GO terms that could be putatively attributed to late spring phenological shifts. Our conservative results indicate that extended transcriptome-based monitoring of forests can capture diverse ranges of responses including air quality, chronic disease, as well as herbivore outbreaks that require activation and/or downregulation of genes collectively tuning reaction norms maintaining the survival of long living trees such as the American beech.

Tree post-drought recovery: scenarios, regulatory mechanisms and ways to improve

Efficient post-drought recovery of growth and assimilation enables a plant to return to its undisturbed state and functioning. Unlike annual plants, trees suffer not only from the current drought, but also from cumulative impacts of consecutive water stresses which cause adverse legacy effects on survival and performance. This review provides an integrated assessment of ecological, physiological and molecular evidence on the recovery of growth and photosynthesis in trees, with a view to informing the breeding of trees with a better ability to recover from water stress. Suppression of recovery processes can result not only from stress damage but also from a controlled downshift of recovery as part of tree acclimation to water-limited conditions. In the latter case, recovery processes could potentially be activated by turning off the controlling mechanisms, but several obstacles make this unlikely. Tree phenology, and specifically photoperiodic constraints, can limit post-drought recovery of growth and photosynthesis, and targeting these constraints may represent a promising way to breed trees with an enhanced ability to recover post-drought. The mechanisms of photoperiod-dependent regulation of shoot, secondary and root growth and of assimilation processes are reviewed. Finally, the limitations and trade-offs of altering the photoperiodic regulation of growth and assimilation processes are discussed.

Life under the snow: A year-round transcriptome analysis of Antarctic mosses in natural habitats provides insight into the molecular adaptation of plants under extreme environment

Mosses are vital components of ecosystems, exhibiting remarkable adaptability across diverse habitats from deserts to polar ice caps. Sanionia uncinata (Hedw.) Loeske, a dominant Antarctic moss survives extreme environmental condition through perennial lifecycles involving growth and dormancy alternation. This study explores genetic controls and molecular mechanisms enabling S. uncinata to cope with seasonality of the Antarctic environment. We analysed the seasonal transcriptome dynamics of S. uncinata collected monthly from February 2015 to January 2016 in King George Island, Antarctica. Findings indicate that genes involved in plant growth were predominantly upregulated in Antarctic summer, while those associated with protein synthesis and cell cycle showed marked expression during the winter-to-summer transition. Genes implicated in cellular stress and abscisic acid signalling were highly expressed in winter. Further, validation included a comparison of the Antarctic field transcriptome data with controlled environment simulation of Antarctic summer and winter temperatures, which revealed consistent gene expression patterns in both datasets. This proposes a seasonal gene regulatory model of S. uncinate to understand moss adaptation to extreme environments. Additionally, this data set is a valuable resource for predicting genetic responses to climatic fluctuations, enhancing our knowledge of Antarctic flora's resilience to global climate change.

Seasonal gene expression signatures of delayed fertilization in Fagaceae

In the family Fagaceae, fertilization is delayed by several weeks to 1 year after pollination, leading to 1- or 2-year fruiting species depending on whether fruiting occurs in the same or the next year after flowering. To investigate physiological responses underlying the regulation of delayed fertilization, we monitored seasonal changes in genome-wide gene expression in tissues including leaves and buds over 2 years under natural conditions in one- (Quercus glauca) and 2-year fruiting species (Lithocarpus edulis). Genes associated with metabolic changes in response to winter cold, photosynthesis and cell proliferation, which are essential for survival and growth, showed highly conserved seasonal expression profiles between species. However, seasonal expression profiles diverged between species in genes associated with pollination, an important process contributing to the origin and maintenance of the reproductive barrier between plant species. By comparing seasonal progression of ovule development and gene expression in pistillate flowers, we revealed that ovules started developing after winter in the 2-year fruiting species, which could be linked to the activation of genes involved in fertilization and female gametophyte development after winter. These findings suggest that the 2-year fruiting species may have evolved a requirement of winter cold to prevent fertilization before winter and facilitate fertilization and embryo development in the following spring when temperature rises. This study offers new possibilities to explore the evolution of reproductive strategies in Fagaceae.

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.

Exploring the Seasonal Dynamics and Molecular Mechanism of Wood Formation in Gymnosperm Trees

Forests, comprising 31% of the Earth's surface, play pivotal roles in regulating the carbon, water, and energy cycles. Despite being far less diverse than angiosperms, gymnosperms account for over 50% of the global woody biomass production. To sustain growth and development, gymnosperms have evolved the capacity to sense and respond to cyclical environmental signals, such as changes in photoperiod and seasonal temperature, which initiate growth (spring and summer) and dormancy (fall and winter). Cambium, the lateral meristem responsible for wood formation, is reactivated through a complex interplay among hormonal, genetic, and epigenetic factors. Temperature signals perceived in early spring induce the synthesis of several phytohormones, including auxins, cytokinins, and gibberellins, which in turn reactivate cambium cells. Additionally, microRNA-mediated genetic and epigenetic pathways modulate cambial function. As a result, the cambium becomes active during the summer, resulting in active secondary xylem (i.e., wood) production, and starts to become inactive in autumn. This review summarizes and discusses recent findings regarding the climatic, hormonal, genetic, and epigenetic regulation of wood formation in gymnosperm trees (i.e., conifers) in response to seasonal changes.

Leaf-branch vulnerability segmentation occurs all year round for three temperate evergreen tree species

Vulnerability segmentation (VS) and Hydraulic segmentation (HS) hypotheses propose higher hydraulic resistance and vulnerability to embolism in leaves than in branches, respectively. The VS and HS are suggested as an acclimation strategy of trees to drought stress, but whether they occur during freezing stress has rarely been explored. We measured the leaf and branch hydraulic traits of three temperate evergreen tree species [Picea koraiensis (Korean spruce), Pinus koraiensis (Korean pine), and Pinus sylvestris var. mongolica (Mongolian pine)] during four seasons (winter, spring, summer, and autumn) across the year. We assessed the applicability of VS and HS all year round, particularly in winter. The water potential at which leaf hydraulic conductance lost 50% (P_50L), was more negative in winter than in summer, while higher leaf mass per area was obtained in winter. These results suggest that these species invest more carbon into leaf (including hydraulic systems) to acclimate to winter frost drought. Leaf and branch hydraulic conductance (K_mL and K_mB) were lower, and the percentage loss of branch hydraulic conductance (PLC_B) was higher in spring than in autumn. These results were probably because of more freeze-thaw cycles in spring (69 cycles) than in autumn (37 cycles). The water potential at which branch hydraulic conductance lost 50%, P_50B, was more negative than P_50L across the year. The values of VS (P_50L minus P_50B) were positive, i.e. leaf was more vulnerable than the branch in all species and across seasons, with higher values occurring in spring or autumn. However, K_mL positively correlated with K_mB, suggesting hydraulic coordination between leaf and branch, but did not support HS. Our findings indicate that leaf-branch vulnerability segmentation can occur all year round, including freezing stress, to protect branches from hydraulic failure in temperate evergreen conifers.

Plant-Programmed Cell Death-Associated Genes Participation in Pinus sylvestris L. Trunk Tissue Formation

Molecular genetic markers of various PCD (programmed cell death) variants during xylo- and phloemogenesis have been identified for the first time in Scots pine under lingonberry pine forest conditions in Northwest Russia (middle taiga subzone). PCD is a genetically determined process. Gene profiles of serine and cysteine proteases (endopeptidases), endonucleases, and metacaspases families are often considered markers of the final xylogenesis stage. In the present study, we examined the gene expression profiles of the BFN (bifunctional endonuclease) family-BFN, BFN1, BFN2, BFN3, and peptidase (cysteine endopeptidase, CEP and metacaspase, MC5) in the radial row, in addition to the vascular phloem and cambium (F1), differentiating xylem (F2), sapwood (SW), and transition zone during the active cambial growth period of uneven-aged pine trees (25-, 63- and 164-cambial age (c.a.) years old). We have shown that the expression patterns of the PCD-related genes did not depend on the cambial age but were largely determined by plant tissue type. In the radial row F1-F2-SW, we studied the activities of enzymes, including sucrose in metabolism (sucrose synthase, three forms of invertase); antioxidant system (AOS) enzymes (superoxide dismutase, catalase); and peroxidase andpolyphenol oxidase, which belonged to AOS enzymes and were involved in the synthesis of phenolic components of cell walls. The activity of the enzymes indicated that the trunk tissues of pine trees had varying metabolic status. Molecular genetic PCD regulation mechanisms during xylem vascular and mechanical element formation and parenchyma cells' PCD during the formation of Scots pine heartwood were discussed.

Regulation of PaRBOH1-mediated ROS production in Norway spruce by Ca2+ binding and phosphorylation

Plant respiratory burst oxidase homologs (RBOHs) are plasma membrane-localized NADPH oxidases that generate superoxide anion radicals, which then dismutate to H₂O₂, into the apoplast using cytoplasmic NADPH as an electron donor. PaRBOH1 is the most highly expressed RBOH gene in developing xylem as well as in a lignin-forming cell culture of Norway spruce (Picea abies L. Karst.). Since no previous information about regulation of gymnosperm RBOHs exist, our aim was to resolve how PaRBOH1 is regulated with a focus on phosphorylation. The N-terminal part of PaRBOH1 was found to contain several putative phosphorylation sites and a four-times repeated motif with similarities to the Botrytis-induced kinase 1 target site in Arabidopsis AtRBOHD. Phosphorylation was indicated for six of the sites in in vitro kinase assays using 15 amino-acid-long peptides for each of the predicted phosphotarget site in the presence of protein extracts of developing xylem. Serine and threonine residues showing positive response in the peptide assays were individually mutated to alanine (kinase-inactive) or to aspartate (phosphomimic), and the wild type PaRBOH1 and the mutated constructs transfected to human kidney embryogenic (HEK293T) cells with a low endogenous level of extracellular ROS production. ROS-producing assays with HEK cells showed that Ca²⁺ and phosphorylation synergistically activate the enzyme and identified several serine and threonine residues that are likely to be phosphorylated including a novel phosphorylation site not characterized in other plant species. These were further investigated with a phosphoproteomic study. Results of Norway spruce, the first gymnosperm species studied in relation to RBOH regulation, show that regulation of RBOH activity is conserved among seed plants.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

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.

Full-length transcriptomic identification of R2R3-MYB family genes related to secondary cell wall development in Cunninghamia lanceolata (Chinese fir)

BACKGROUND

R2R3-MYB is a class of transcription factor crucial in regulating secondary cell wall development during wood formation. The regulation of wood formation in gymnosperm has been understudied due to its large genome size. Using Single-Molecule Real-Time sequencing, we obtained full-length transcriptomic libraries from the developmental stem of Cunninghamia lanceolata, a perennial conifer known as Chinese fir. The R2R3-MYB of C. lanceolata (hereafter named as ClMYB) associated with secondary wall development were identified based on phylogenetic analysis, expression studies and functional study on transgenic line.

RESULTS

The evolutionary relationship of 52 ClMYBs with those from Arabidopsis thaliana, Eucalyptus grandis, Populus trichocarpa, Oryza sativa, two gymnosperm species, Pinus taeda, and Picea glauca were established by neighbour-joining phylogenetic analysis. A large number of ClMYBs resided in the woody-expanded subgroups that predominated with the members from woody dicots. In contrast, the woody-preferential subgroup strictly carrying the members of woody dicots contained only one candidate. The results suggest that the woody-expanded subgroup emerges before the gymnosperm/angiosperm split, while most of the woody-preferential subgroups are likely lineage-specific to woody dicots. Nine candidates shared the same subgroups with the A. thaliana orthologs, with known function in regulating secondary wall development. Gene expression analysis inferred that ClMYB1/2/3/4/5/26/27/49/51 might participate in secondary wall development, among which ClMYB1/2/5/26/27/49 were significantly upregulated in the highly lignified compression wood region, reinforcing their regulatory role associated with secondary wall development. ClMYB1 was experimentally proven a transcriptional activator that localised in the nucleus. The overexpression of ClMYB1 in Nicotiana benthamiana resulted in an increased lignin deposition in the stems. The members of subgroup S4, ClMYB3/4/5 shared the ERF-associated amphiphilic repression motif with AtMYB4, which is known to repress the metabolism of phenylpropanoid derived compounds. They also carried a core motif specific to gymnosperm lineage, suggesting divergence of the regulatory process compared to the angiosperms.

CONCLUSIONS

This work will enrich the collection of full-length gymnosperm-specific R2R3-MYBs related to stem development and contribute to understanding their evolutionary relationship with angiosperm species.

An atlas of the Norway spruce needle seasonal transcriptome

Boreal conifers possess a tremendous ability to survive and remain evergreen during harsh winter conditions and resume growth during summer. This is enabled by coordinated regulation of major cellular functions at the level of gene expression, metabolism, and physiology. Here we present a comprehensive characterization of the annual changes in the global transcriptome of Norway spruce (Picea abies) needles as a resource to understand needle development and acclimation processes throughout the year. In young, growing needles (May 15 until June 30), cell walls, organelles, etc., were formed, and this developmental program heavily influenced the transcriptome, explained by over-represented Gene Ontology (GO) categories. Later changes in gene expression were smaller but four phases were recognized: summer (July-August), autumn (September-October), winter (November-February), and spring (March-April), where over-represented GO categories demonstrated how the needles acclimated to the various seasons. Changes in the seasonal global transcriptome profile were accompanied by differential expression of members of the major transcription factor families. We present a tentative model of how cellular activities are regulated over the year in needles of Norway spruce, which demonstrates the value of mining this dataset, accessible in ConGenIE together with advanced visualization tools.

Concerted control of the LaRAV1-LaCDKB1;3 module by temperature during dormancy release and reactivation of larch

Dormancy release and reactivation of temperate-zone trees involve the temperature-modulated expression of cell-cycle genes. However, information on the detailed regulatory mechanism is limited. Here, we compared the transcriptomes of the stems of active and dormant larch trees, emphasizing the expression patterns of cell-cycle genes and transcription factors and assessed their relationships and responses to temperatures. Twelve cell-cycle genes and 31 transcription factors were strongly expressed in the active stage. Promoter analysis suggested that these 12 genes might be regulated by transcription factors from 10 families. Altogether, 73 cases of regulation between 16 transcription factors and 12 cell-cycle genes were predicted, while the regulatory interactions between LaMYB20 and LaCYCB1;1, and LaRAV1 and LaCDKB1;3 were confirmed by yeast one-hybrid and dual-luciferase assays. Last, we found that LaRAV1 and LaCDKB1;3 had almost the same expression patterns during dormancy release and reactivation induced naturally or artificially by temperature, indicating that the LaRAV1-LaCDKB1;3 module functions in the temperature-modulated dormancy release and reactivation of larch trees. These results provide new insights into the link between temperature and cell-cycle gene expression, helping to understand the temperature control of tree growth and development in the context of climate change.

Investigation Into Different Wood Formation Mechanisms Between Angiosperm and Gymnosperm Tree Species at the Transcriptional and Post-transcriptional Level

Enormous distinctions of the stem structure and cell types between gymnosperms and angiosperms tree species are expected to cause quite different wood physical and mechanical attributes, however, the molecular mechanisms underlying the differing wood morphology are still unclear. In this study, we compared the transcriptomes obtained by RNA-Seq between Populus alba × P. glandulosa clone 84K, and Larix kaempferi (Lamb.) Carr trees. Available genome resource served as reference for P. alba × P. glandulosa and the Iso-Seq results of a three-tissues mixture (xylem, phloem, and leaf) were used as the reference for L. kaempferi to compare the xylem-specifically expressed genes and their alternative splicing model. Through screening, we obtained 13,907 xylem-specifically expressed genes (5,954 up-regulated, 7,953 down-regulated) in the xylem of P. alba × P. glandulosa, and 2,596 xylem-specifically expressed genes (1,648 up-regulated, 948 down-regulated) in the xylem of L. kaempferi. From the GO and KEGG analyses, some genes associated with two wood formation-related pathways, namely those for phenylpropanoid biosynthesis, and starch and sucrose metabolism, were successfully screened. Then the distributions and gene expression models between P. alba × P. glandulosa and L. kaempferi in those pathways were compared, which suggested differential wood formation processes between the angiosperm and gymnosperm trees. Furthermore, a Weight Gene Co-expression Network Analysis (WGCNA) for total xylem-specifically expressed genes in two species was conducted, from which wood formation-related modules were selected to build a co-expression network for the two tree species. The genes within this co-expression network showed different co-expression relationships between the angiosperm and gymnosperm woody species. Comparing the alternative splicing events for wood formation-related genes suggests a different post-transcriptional regulation process exists between the angiosperm and gymnosperm trees. Our research thus provides the foundation for the in-depth investigation of different wood formation mechanisms of angiosperm and gymnosperm species.

Synergies and Entanglement in Secondary Cell Wall Development and Abiotic Stress Response in Trees

A major challenge for sustainable food, fuel, and fiber production is simultaneous genetic improvement of yield, biomass quality, and resilience to episodic environmental stress and climate change. For Populus and other forest trees, quality traits involve alterations in the secondary cell wall (SCW) of wood for traditional uses, as well as for a growing diversity of biofuels and bioproducts. Alterations in wood properties that are desirable for specific end uses can have negative effects on growth and stress tolerance. Understanding of the diverse roles of SCW genes is necessary for the genetic improvement of fast-growing, short-rotation trees that face perennial challenges in their growth and development. Here, we review recent progress into the synergies and antagonisms of SCW development and abiotic stress responses, particularly, the roles of transcription factors, SCW biogenesis genes, and paralog evolution.

Transcriptome and microRNA Sequencing Identified miRNAs and Target Genes in Different Developmental Stages of the Vascular Cambium in Cryptomeria fortunei Hooibrenk

Cryptomeria fortunei Hooibrenk is an important fast-growing coniferous timber species that is widely used in landscaping. Recently, research on timber quality has gained substantial attention in the field of tree breeding. Wood is the secondary xylem formed by the continuous inward division and differentiation of the vascular cambium; therefore, the development of the vascular cambium is particularly important for wood quality. In this study, we analyzed the transcriptomes of the cambial zone in C. fortunei during different developmental stages using Illumina HiSeq sequencing, focusing on general transcriptome and microRNA (miRNA) data. We performed functional annotation of the differentially expressed genes (DEGs) in the different stages identified by transcriptome sequencing and generated 15 miRNA libraries yielding 4.73 Gb of clean reads. The most common length of the filtered miRNAs was 21nt, accounting for 33.1% of the total filtered reads. We annotated a total of 32 known miRNA families. Some miRNAs played roles in hormone signal transduction (miR159, miR160, and miR166), growth and development (miR166 and miR396), and the coercion response (miR394 and miR395), and degradome sequencing showed potential cleavage sites between miRNAs and target genes. Differential expression of miRNAs and target genes and functional validation of the obtained transcriptome and miRNA data provide a theoretical basis for further elucidating the molecular mechanisms of cellular growth and differentiation, as well as wood formation in the vascular cambium, which will help improve the wood quality of C. fortunei.

Hunting monolignol transporters: membrane proteomics and biochemical transport assays with membrane vesicles of Norway spruce

Both the mechanisms of monolignol transport and the transported form of monolignols in developing xylem of trees are unknown. We tested the hypothesis of an active, plasma membrane-localized transport of monolignol monomers, dimers, and/or glucosidic forms with membrane vesicles prepared from developing xylem and lignin-forming tissue-cultured cells of Norway spruce (Picea abies L. Karst.), as well as from control materials, comprising non-lignifying Norway spruce phloem and tobacco (Nicotiana tabacum L.) BY-2 cells. Xylem and BY-2 vesicles transported both coniferin and p-coumaryl alcohol glucoside, but inhibitor assays suggested that this transport was through the tonoplast. Membrane vesicles prepared from lignin-forming spruce cells showed coniferin transport, but the Km value for coniferin was much higher than those of xylem and BY-2 cells. Liquid chromatography-mass spectrometry analysis of membrane proteins isolated from spruce developing xylem, phloem, and lignin-forming cultured cells revealed multiple transporters. These were compared with a transporter gene set obtained by a correlation analysis with a selected set of spruce monolignol biosynthesis genes. Biochemical membrane vesicle assays showed no support for ABC-transporter-mediated monolignol transport but point to a role for secondary active transporters (such as MFS or MATE transporters). In contrast, proteomic and co-expression analyses suggested a role for ABC transporters and MFS transporters.

Integrative analysis of wood biomass and developing xylem transcriptome provide insights into mechanisms of lignin biosynthesis in wood formation of Pinus massoniana

Lignin is an important renewable energy source as an excellent new battery fuel and ideal substitutes for the petrochemical industry. However, the molecular mechanism underlying lignin biosynthesis in wood formation of P. massoniana remains unexplored. Thus, an integrative analysis of wood biomass and the developing xylem transcriptome was performed to identify genes involved in lignin biosynthesis. A total of 1624 differentially expressed genes (DEGs) were identified, consisting of 797 upregulated and 827 downregulated genes (MaxG vs MinG). Additionally, 122 candidate genes and 17 DEGs were successfully annotated to the lignin biosynthesis pathway. All upregulated MYB and NAC genes were regulators of secondary cell wall formation. Moreover, the qRT-PCR analyses shown that 9 lignin biosynthesis-related genes and 7 transcription factor-encoding genes were upregulated (MaxG vs MinG), which indicated that the downregulation of lignin biosynthesis-related genes might be the possible causes of growth retardation and dwarf phenotype in some P. massoniana individuals. The identification of lignin biosynthesis-related genes can provide valuable genetic basis and resource for further researches on molecular mechanisms of lignin biosynthesis and contribute to the future investigations of bioengineering and synthetic biology to regulate lignin content in wood formation for the pulp and wood utilization industry.

Transcriptomic Profiling of Cryptomeria fortunei Hooibrenk Vascular Cambium Identifies Candidate Genes Involved in Phenylpropanoid Metabolism

Cryptomeria fortunei Hooibrenk (Chinese cedar) is a coniferous tree from southern China that has an important function in landscaping and timber production. Lignin is one of the key components of secondary cell walls, which have a crucial role in conducting water and providing mechanical support for the upward growth of plants. It is mainly biosynthesized via the phenylpropanoid metabolic pathway, of which the molecular mechanism remains so far unresolved in C. fortunei. In order to obtain further insight into this pathway, we performed transcriptome sequencing of the C. fortunei cambial zone at 5 successive growth stages. We generated 78,673 unigenes from transcriptome data, of which 45,214 (57.47%) were successfully annotated in the non-redundant protein database (NR). A total of 8975 unigenes were identified to be significantly differentially expressed between Sample_B and Sample_A after analyzing their expression profiles. Of the differentially expressed genes (DEGs), 6817 (75.96%) and 2158 (24.04%) were up- and down-regulated, respectively. 83 DEGs were involved in phenylpropanoid metabolism, 37 DEGs that encoded v-Myb avian myeloblastosis viral oncogene homolog (MYB) transcription factor (TF), and many candidates that encoded lignin synthesizing enzymes. These findings contribute to understanding the expression pattern of C. fortunei cambial zone transcriptome. Furthermore, our results provide additional insight towards understanding the molecular mechanisms of wood formation in C. fortunei.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

Ray Parenchymal Cells Contribute to Lignification of Tracheids in Developing Xylem of Norway Spruce

A comparative transcriptomic study and a single-cell metabolome analysis were combined to determine whether parenchymal ray cells contribute to the biosynthesis of monolignols in the lignifying xylem of Norway spruce (Picea abies). Ray parenchymal cells may function in the lignification of upright tracheids by supplying monolignols. To test this hypothesis, parenchymal ray cells and upright tracheids were dissected with laser-capture microdissection from tangential cryosections of developing xylem of spruce trees. The transcriptome analysis revealed that among the genes involved in processes typical for vascular tissues, genes encoding cell wall biogenesis-related enzymes were highly expressed in both developing tracheids and ray cells. Interestingly, most of the shikimate and monolignol biosynthesis pathway-related genes were equally expressed in both cell types. Nonetheless, 1,073 differentially expressed genes were detected between developing ray cells and tracheids, among which a set of genes expressed only in ray cells was identified. In situ single cell metabolomics of semi-intact plants by picoliter pressure probe-electrospray ionization-mass spectrometry detected monolignols and their glycoconjugates in both cell types, indicating that the biosynthetic route for monolignols is active in both upright tracheids and parenchymal ray cells. The data strongly support the hypothesis that in developing xylem, ray cells produce monolignols that contribute to lignification of tracheid cell walls.

How to organize bidirectional tissue production?

The cambium is a plant-borne stem cell system producing wood and bast, two distinct types of vascular tissues, in strictly opposite directions. Thereby, the cambium contributes substantially to terrestrial biomass accumulation and represents the basis for the formation of large plant bodies. Although the bidirectional mode of tissue production by a common stem cell pool holds interesting implications for developmental biology, functional domains of the cambium, and their interaction remained poorly defined for decades. Here, we summarize recent findings on domain organization of the cambium and discuss potential mechanisms important for its bipartite organization. By highlighting the conceptual implication for stem cell biology, we integrate our understanding of cambium regulation into a larger context.

Cloning and Functional Analysis of Lignin Biosynthesis Genes Cf4CL and CfCCoAOMT in Cryptomeria fortunei

Cryptomeria fortunei, also known as the Chinese cedar, is an important timber species in southern China. The primary component of its woody tissues is lignin, mainly present in secondary cell walls. Therefore, continuous lignin synthesis is crucial for wood formation. In this study, we aimed to discover key genes involved in lignin synthesis expressed in the vascular cambium of C. fortunei. Through transcriptome sequencing, we detected expression of two genes, 4CL and CCoAOMT, known to be homologous to enzymes involved in the lignin synthesis pathway. We studied the function of these genes through bioinformatics analysis, cloning, vascular cambium expression analysis, and transgenic cross-species functional validation studies. Our results show that Cf4CL and CfCCoAOMT do indeed function in the pathway of lignin synthesis and likely perform this function in C. fortunei. They are prime candidates for future (gene-editing) studies aimed at optimizing C. fortunei wood production.

Finding New Cell Wall Regulatory Genes in Populus trichocarpa Using Multiple Lines of Evidence

Understanding the regulatory network controlling cell wall biosynthesis is of great interest in Populus trichocarpa, both because of its status as a model woody perennial and its importance for lignocellulosic products. We searched for genes with putatively unknown roles in regulating cell wall biosynthesis using an extended network-based Lines of Evidence (LOE) pipeline to combine multiple omics data sets in P. trichocarpa, including gene coexpression, gene comethylation, population level pairwise SNP correlations, and two distinct SNP-metabolite Genome Wide Association Study (GWAS) layers. By incorporating validation, ranking, and filtering approaches we produced a list of nine high priority gene candidates for involvement in the regulation of cell wall biosynthesis. We subsequently performed a detailed investigation of candidate gene GROWTH-REGULATING FACTOR 9 (PtGRF9). To investigate the role of PtGRF9 in regulating cell wall biosynthesis, we assessed the genome-wide connections of PtGRF9 and a paralog across data layers with functional enrichment analyses, predictive transcription factor binding site analysis, and an independent comparison to eQTN data. Our findings indicate that PtGRF9 likely affects the cell wall by directly repressing genes involved in cell wall biosynthesis, such as PtCCoAOMT and PtMYB.41, and indirectly by regulating homeobox genes. Furthermore, evidence suggests that PtGRF9 paralogs may act as transcriptional co-regulators that direct the global energy usage of the plant. Using our extended pipeline, we show multiple lines of evidence implicating the involvement of these genes in cell wall regulatory functions and demonstrate the value of this method for prioritizing candidate genes for experimental validation.

Domain swap between two type-II metacaspases defines key elements for their biochemical properties

Type-II metacaspases are conserved cysteine proteases found in eukaryotes with oxygenic photosynthesis, including green plants and some algae, such as Chlamydomonas and Volvox. Genetic and biochemical studies showed that some members in this protease family could be involved in oxidative stress-induced cell death in higher plants, but their regulatory mechanisms remain unclear. Biochemically, two distinct classes of type-II metacaspases are exemplified by AtMC4 and AtMC9 from Arabidopsis, with AtMC4 activation dependent on calcium under neutral pH, whereas AtMC9 is active only under mildly acidic pH, regardless of the availability of calcium. Here, we constructed all six possible combinations between the p20, linker, and p10 domains from AtMC4 and AtMC9. Our results show that calcium stimulation of AtMC4 activity is associated with essential amino acids located in its p20 domain. In contrast, the acidic pH optimum trait is lost from AtMC9 if one or two of its domains are replaced by that from AtMC4, suggesting that multiple interactions between domains in AtMC9 may be responsible for this property. Consistent with this, we found conserved 'signature' residues in each of the three domains that distinguish AtMC4- and AtMC9-like classes of metacaspases. Tracing the origin of the AtMC9 class, we found evidence for its appearance between lycophytes and gymnosperms, coincident with the evolution of more complex root archetypes in terrestrial plants. Our work suggests that the distinctive properties of the AtMC9-like protease could be associated with special cellular physiology in the roots of gymnosperms and angiosperms that are distinct from lycophytes.

Co-expression network of transcription factors reveal ethylene-responsive element-binding factor as key regulator of wood phenotype in Eucalyptus tereticornis

Suitability of wood biomass for pulp production is dependent on the cellular architecture and composition of secondary cell wall. Presently, systems genetics approach is being employed to understand the molecular basis of trait variation and co-expression network analysis has enabled holistic understanding of complex trait such as secondary development. Transcription factors (TFs) are reported as key regulators of meristematic growth and wood formation. The hierarchical TF network is a multi-layered system which interacts with downstream structural genes involved in biosynthesis of cellulose, hemicelluloses and lignin. Several TFs have been associated with wood formation in tree species such as Populus, Eucalyptus, Picea and Pinus. However, TF-specific co-expression networks to understand the interaction between these regulators are not reported. In the present study, co-expression network was developed for TFs expressed during wood formation in Eucalyptus tereticornis and ethylene-responsive element-binding factor, EtERF2, was identified as the major hub transcript which co-expressed with other secondary cell wall biogenesis-specific TFs such as EtSND2, EtVND1, EtVND4, EtVND6, EtMYB70, EtGRAS and EtSCL8. This study reveals a probable role of ethylene in determining natural variation in wood properties in Eucalyptus species. Understanding this transcriptional regulation underpinning the complex bio-processing trait of wood biomass will complement the Eucalyptus breeding program through selection of industrially suitable phenotypes by marker-assisted selection.