AspWood: High-Spatial-Resolution Transcriptome Profiles Reveal Uncharacterized Modularity of Wood Formation in Populus tremula

Time-sequential detection of quantitative trait loci and candidate genes underlying the dynamic growth of Salix suchowensis

Elucidating the genetic factors underlying long-term biological processes remains challenging since the relevant genes and their effects may vary across different developmental stages. In this study, we carried out a large-scale field trial of the progeny of an F1 full-sib pedigree of Salix suchowensis and measured plant height and ground diameter periodically over a time course of 240 days. With the obtained data, we characterized plant growth rhythms and performed time-sequential analyses of quantitative trait loci underlying the dynamic growth of the plants. The dynamic mapping of quantitative trait loci revealed that stem height and ground diameter were under the control of four quantitative trait loci, and the effects of these quantitative trait loci varied greatly throughout the growth process, in which two quantitative trait loci were found to exert a pleiotropic effect determining the correlation between stem height and ground diameter. The analysis of candidate genes in the target genetic intervals showed that the pleiotropic effect of the two quantitative trait loci arises from the colocalization of genes with independent effects on stem height and ground diameter. Further examination of the expression patterns of the candidate genes indicated that height and circumference growth involve different activities of leaf and cambium tissues. This study provides unprecedented information to help us understand the dynamic growth of plants and presents an applicable strategy for elucidating the genetic mechanism underlying a long-term biological process by using plant growth as an example.

Wood Formation Modeling - A Research Review and Future Perspectives

Wood formation has received considerable attention across various research fields as a key process to model. Historical and contemporary models of wood formation from various disciplines have encapsulated hypotheses such as the influence of external (e.g., climatic) or internal (e.g., hormonal) factors on the successive stages of wood cell differentiation. This review covers 17 wood formation models from three different disciplines, the earliest from 1968 and the latest from 2020. The described processes, as well as their external and internal drivers and their level of complexity, are discussed. This work is the first systematic cataloging, characterization, and process-focused review of wood formation models. Remaining open questions concerning wood formation processes are identified, and relate to: (1) the extent of hormonal influence on the final tree ring structure; (2) the mechanism underlying the transition from earlywood to latewood in extratropical regions; and (3) the extent to which carbon plays a role as "active" driver or "passive" substrate for growth. We conclude by arguing that wood formation models remain to be fully exploited, with the potential to contribute to studies concerning individual tree carbon sequestration-storage dynamics and regional to global carbon sequestration dynamics in terrestrial vegetation models.

Mobile forms of carbon in trees: metabolism and transport

Plants constitute 80% of the biomass on earth, and almost two-thirds of this biomass is found in wood. Wood formation is a carbon (C)-demanding process and relies on C transport from photosynthetic tissues. Thus, understanding the transport process is of major interest for understanding terrestrial biomass formation. Here, we review the molecules and mechanisms used to transport and allocate C in trees. Sucrose is the major form in which C is transported in plants, and it is found in the phloem sap of all tree species investigated so far. However, in several tree species, sucrose is accompanied by other molecules, notably polyols and the raffinose family of oligosaccharides. We describe the molecules that constitute each of these transport groups, and their distribution across different tree species. Furthermore, we detail the metabolic reactions for their synthesis, the mechanisms by which trees load and unload these compounds in and out of the vascular system, and how they are radially transported in the trunk and finally catabolized during wood formation. We also address a particular C recirculation process between phloem and xylem that occurs in trees during the annual cycle of growth and dormancy. A search of possible evolutionary drivers behind the diversity of C-carrying molecules in trees reveals no consistent differences in C transport mechanisms between angiosperm and gymnosperm trees. Furthermore, the distribution of C forms across species suggests that climate-related environmental factors will not explain the diversity of C transport forms. However, the consideration of C-transport mechanisms in relation to tree-rhizosphere coevolution deserves further attention. To conclude the review, we identify possible future lines of research in this field.

Characterization and Interaction Analysis of the Secondary Cell Wall Synthesis-Related Transcription Factor PmMYB7 in Pinus massoniana Lamb

In vascular plants, the importance of R2R3-myeloblastosis (R2R3-MYB) transcription factors (TFs) in the formation of secondary cell walls (SCWs) has long been a controversial topic due to the lack of empirical evidence of an association between TFs and downstream target genes. Here, we found that the transcription factor PmMYB7, which belongs to the R2R3-MYB subfamily, is involved in lignin biosynthesis in Pinus massoniana. PmMYB7 was highly expressed in lignified tissues and upon abiotic stress. As a bait carrier, the PmMYB7 protein had no toxicity or autoactivation in the nucleus. Forty-seven proteins were screened from the P. massoniana yeast library. These proteins were predicted to be mainly involved in resistance, abiotic stress, cell wall biosynthesis, and cell development. We found that the PmMYB7 protein interacted with caffeoyl CoA 3-O-methyltransferase-2 (PmCCoAOMT2)—which is involved in lignin biosynthesis—but not with beta-1, 2-xylosyltransferase (PmXYXT1) yeast two-hybrid (Y2H) studies. Our in vivo coimmunoprecipitation (Co-IP) assay further showed that the PmMYB7 and PmCCoAOMT2 proteins could interact. Therefore, we concluded that PmMYB7 is an upstream TF that can interact with PmCCoAOMT2 in plant cells. These findings lay a foundation for further research on the function of PmMYB7, lignin biosynthesis and molecular breeding in P. massoniana.

pHBMT1, a BAHD-family monolignol acyltransferase, mediates lignin acylation in poplar

Poplar (Populus) lignin is naturally acylated with p-hydroxybenzoate ester moieties. However, the enzyme(s) involved in the biosynthesis of the monolignol-p-hydroxybenzoates have remained largely unknown. Here, we performed an in vitro screen of the Populus trichocarpa BAHD acyltransferase superfamily (116 genes) using a wheatgerm cell-free translation system and found five enzymes capable of producing monolignol-p-hydroxybenzoates. We then compared the transcript abundance of the five corresponding genes with p-hydroxybenzoate concentrations using naturally occurring unrelated genotypes of P. trichocarpa and revealed a positive correlation between the expression of p-hydroxybenzoyl-CoA monolig-nol transferase (pHBMT1, Potri.001G448000) and p-hydroxybenzoate levels. To test whether pHBMT1 is responsible for the biosynthesis of monolignol-p-hydroxybenzoates, we overexpressed pHBMT1 in hybrid poplar (Populus alba × P. grandidentata) (35S::pHBMT1 and C4H::pHBMT1). Using three complementary analytical methods, we showed that there was an increase in soluble monolignol-p-hydroxybenzoates and cell-wall-bound monolignol-p-hydroxybenzoates in the poplar transgenics. As these pendent groups are ester-linked, saponification releases p-hydroxybenzoate, a precursor to parabens that are used in pharmaceuticals and cosmetics. This identified gene could therefore be used to engineer lignocellulosic biomass with increased value for emerging biorefinery strategies.

Laying it on thick: a study in secondary growth

The development of secondary vascular tissue enhances the transport capacity and mechanical strength of plant bodies, while contributing a huge proportion of the world's biomass in the form of wood. Cell divisions in the cambium, which constitutes the vascular meristem, provide progenitors from which conductive xylem and phloem are derived. The cambium is a somewhat unusual stem cell population in two respects, making it an interesting subject for developmental research. Firstly, it arises post-germination, and thus represents a model for understanding stem cell initiation beyond embryogenesis. Secondly, xylem and phloem differentiate on opposing sides of cambial stem cells, making them bifacial in nature. Recent discoveries in Arabidopsis thaliana have provided insight into the molecular mechanisms that regulate the initiation, patterning, and maintenance of the cambium. In this review, the roles of intercellular signalling via mobile transcription factors, peptide-receptor modules, and phytohormones are described. Crosstalk between these regulatory pathways is becoming increasingly apparent, yet the underlying mechanisms are not fully understood. Future study of the interaction between multiple independently identified regulators, as well as the functions of their orthologues in trees, will deepen our understanding of radial growth in plants.

A High-Throughput Screening System for Populus Wood-Associated Transcription Factors and Its Application to Lignin Regulation

Wood formation of trees is a complex and costly developmental process, whose regulatory network is involved in the protein-protein and protein-DNA interactions. To detect such interactions in wood development, we developed a high-throughput screening system with 517 Gal4-AD-wood-associated transcription factors (TFs) library from Populus alba × P. glandulosa cv "84K." This system can be used for screening the upstream regulators and interacting proteins of targets by mating-based yeast-one hybrid (Y1H) and yeast-two-hybrid (Y2H) method, respectively. Multiple regulatory modules of lignin biosynthesis were identified based on this Populus system. Five TFs interacted with the 500-bp promoter fragment of PHENYLALANINE AMMONIA-LYASE 2 (PAL2), the first rate-limiting enzyme gene in the lignin biosynthesis pathway, and 10 TFs interacted with PaMYB4/LTF1, a key regulator of lignin biosynthesis. Some of these interactions were further validated by EMSA and BiFC assays. The TF-PaPAL2 promoter interaction and TF-PaMYB4 interaction revealed a complex mechanism governing the regulation of lignin synthesis in wood cells. Our high-throughput Y1H/Y2H screening system may be an efficient tool for studying regulatory network of wood formation in tree species.

The XTH Gene Family in Schima superba: Genome-Wide Identification, Expression Profiles, and Functional Interaction Network Analysis

Xyloglucan endotransglucosylase/hydrolase (XTH), belonging to glycoside hydrolase family 16, is one of the key enzymes in plant cell wall remodeling. Schima superba is an important timber and fireproof tree species in southern China. However, little is known about XTHs in S. superba. In the present study, a total of 34 SsuXTHs were obtained, which were classified into three subfamilies based on the phylogenetic relationship and unevenly distributed on 18 chromosomes. Furthermore, the intron-exon structure and conserved motif composition of them supported the classification and the members belonging to the same subfamily shared similar gene structures. Segmental and tandem duplication events did not lead to SsuXTH gene family expansion, and strong purifying selection pressures during evolution led to similar structure and function of SsuXTH gene family. The interaction network and cis-acting regulatory elements analysis revealed the SsuXTH expression might be regulated by multiple hormones, abiotic stresses and transcription factors. Finally, expression profiles and GO enrichment analysis showed most of the tandem repeat genes were mainly expressed in the phloem and xylem and they mainly participated in glycoside metabolic processes through the transfer and hydrolysis of xyloglucan in the cell wall and then regulated fiber elongation.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BACKGROUND

Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.

RESULTS

We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.

CONCLUSIONS

qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Transcriptional landscape of highly lignified poplar stems at single-cell resolution

BACKGROUND

Plant secondary growth depends on the activity of the vascular cambium, which produces xylem and phloem. Wood derived from xylem is the most abundant form of biomass globally and has played key socio-economic and subsistence roles throughout human history. However, despite intensive study of vascular development, the full diversity of cell types and the gene networks engaged are still poorly understood.

RESULTS

Here, we have applied an optimized protoplast isolation protocol and RNA sequencing to characterize the high-resolution single-cell transcriptional landscape of highly lignified poplar stems. We identify 20 putative cell clusters with a series of novel cluster-specific marker genes and find that these cells are highly heterogeneous based on the transcriptome. Analysis of these marker genes' expression dynamics enables reconstruction of the cell differentiation trajectories involved in phloem and xylem development. We find that different cell clusters exhibit distinct patterns of phytohormone responses and emphasize the use of our data to predict potential gene redundancy and identify candidate genes related to vascular development in trees.

CONCLUSIONS

These findings establish the transcriptional landscape of major cell types of poplar stems at single-cell resolution and provide a valuable resource for investigating basic principles of vascular cell specification and differentiation in trees.

Integrated Transcriptomic and Proteomic Analysis in the Roadmap of the Xylem Development Stage in Populus tomentosa

Xylem development plays an important role in the wood formation of plants. In this study, we found that xylem development was a rapid thickening process characterized by initially rapid increases in the number of tracheary elements and fiber cells and the thickness of the secondary walls that later plateaued. Transcriptome analysis showed that the xylan and lignin biosynthetic pathways, which are involved in the early rapid thickening of the xylem, were mainly upregulated in the second month. The expression of a total of 124 transcription factors (TFs), including 28 NAC TFs and 31 MYB TFs, peaked in 2- and 3-month-old plants compared with 1-month-old plants. Based on previous studies and the key cis-acting elements secondary wall NAC-binding elements, secondary wall MYB-responsive elements, W-box and TGTG[T/G/C], 10 TFs related to xylem development, 50 TFs with unknown function, 98 cell wall biosynthetic genes, and 47 programmed cell death (PCD) genes were used to construct a four-layer transcriptional regulatory network (TRN) with poplar NAC domain TFs to characterize the transcriptional regulation of cell wall biosynthesis and PCD in Populus tomentosa. The proteome revealed that post-transcriptional modification may be widely involved in lignification development. Overall, our results revealed that xylem development is a rapid thickening process in P. tomentosa, and expression patterns varied temporally from cell division to cell death.

Single-cell RNA sequencing reveals a high-resolution cell atlas of xylem in Populus

High-throughput single-cell RNA sequencing (scRNA-seq) has advantages over traditional RNA-seq to explore spatiotemporal information on gene dynamic expressions in heterogenous tissues. We performed Drop-seq, a method for the dropwise sequestration of single cells for sequencing, on protoplasts from the differentiating xylem of Populus alba × Populus glandulosa. The scRNA-seq profiled 9,798 cells, which were grouped into 12 clusters. Through characterization of differentially expressed genes in each cluster and RNA in situ hybridizations, we identified vessel cells, fiber cells, ray parenchyma cells and xylem precursor cells. Diffusion pseudotime analyses revealed the differentiating trajectory of vessels, fiber cells and ray parenchyma cells and indicated a different differentiation process between vessels and fiber cells, and a similar differentiation process between fiber cells and ray parenchyma cells. We identified marker genes for each cell type (cluster) and key candidate regulators during developmental stages of xylem cell differentiation. Our study generates a high-resolution expression atlas of wood formation at the single cell level and provides valuable information on wood formation.

The transcription factor PagLBD3 contributes to the regulation of secondary growth in Populus

LATERAL ORGAN BOUNDARIES DOMAIN (LBD) genes encode plant-specific transcription factors that participate in regulating various developmental processes. In this study, we genetically characterized PagLBD3 encoding an important regulator of secondary growth in poplar (Populus alba × Populus glandulosa). Overexpression of PagLBD3 increased stem secondary growth in Populus with a significantly higher rate of cambial cell differentiation into phloem, while dominant repression of PagLBD3 significantly decreased the rate of cambial cell differentiation into phloem. Furthermore, we identified 1756 PagLBD3 genome-wide putative direct target genes (DTGs) through RNA sequencing (RNA-seq)-coupled DNA affinity purification followed by sequencing (DAP-seq) assays. Gene Ontology analysis revealed that genes regulated by PagLBD3 were enriched in biological pathways regulating meristem development, xylem development, and auxin transport. Several central regulator genes for vascular development, including PHLOEM INTERCALATED WITH XYLEM (PXY), WUSCHEL RELATED HOMEOBOX4 (WOX4), Secondary Wall-Associated NAC Domain 1s (SND1-B2), and Vascular-Related NAC-Domain 6s (VND6-B1), were identified as PagLBD3 DTGs. Together, our results indicate that PagLBD3 and its DTGs form a complex transcriptional network to modulate cambium activity and phloem/xylem differentiation.

Between Stress and Response: Function and Localization of Mechanosensitive Ca2+ Channels in Herbaceous and Perennial Plants

Over the past three decades, how plants sense and respond to mechanical stress has become a flourishing field of research. The pivotal role of mechanosensing in organogenesis and acclimation was demonstrated in various plants, and links are emerging between gene regulatory networks and physical forces exerted on tissues. However, how plant cells convert physical signals into chemical signals remains unclear. Numerous studies have focused on the role played by mechanosensitive (MS) calcium ion channels MCA, Piezo and OSCA. To complement these data, we combined data mining and visualization approaches to compare the tissue-specific expression of these genes, taking advantage of recent single-cell RNA-sequencing data obtained in the root apex and the stem of Arabidopsis and the Populus stem. These analyses raise questions about the relationships between the localization of MS channels and the localization of stress and responses. Such tissue-specific expression studies could help to elucidate the functions of MS channels. Finally, we stress the need for a better understanding of such mechanisms in trees, which are facing mechanical challenges of much higher magnitudes and over much longer time scales than herbaceous plants, and we mention practical applications of plant responsiveness to mechanical stress in agriculture and forestry.

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.

Elongation of wood fibers combines features of diffuse and tip growth

Xylem fibers are highly elongated cells that are key constituents of wood, play major physiological roles in plants, comprise an important terrestrial carbon reservoir, and thus have enormous ecological and economic importance. As they develop, from fusiform initials, their bodies remain the same length while their tips elongate and intrude into intercellular spaces. To elucidate mechanisms of tip elongation, we studied the cell wall along the length of isolated, elongating aspen xylem fibers and used computer simulations to predict the forces driving the intercellular space formation required for their growth. We found pectin matrix epitopes (JIM5, LM7) concentrated at the tips where cellulose microfibrils have transverse orientation, and xyloglucan epitopes (CCRC-M89, CCRC-M58) in fiber bodies where microfibrils are disordered. These features are accompanied by changes in cell wall thickness, indicating that while the cell wall elongates strictly at the tips, it is deposited all over fibers. Computer modeling revealed that the intercellular space formation needed for intrusive growth may only require targeted release of cell adhesion, which allows turgor pressure in neighboring fiber cells to 'round' the cells creating spaces. These characteristics show that xylem fibers' elongation involves a distinct mechanism that combines features of both diffuse and tip growth.

Wood Formation under Severe Drought Invokes Adjustment of the Hormonal and Transcriptional Landscape in Poplar

Drought is a severe environmental stress that exerts negative effects on plant growth. In trees, drought leads to reduced secondary growth and altered wood anatomy. The mechanisms underlying wood stress adaptation are not well understood. Here, we investigated the physiological, anatomical, hormonal, and transcriptional responses of poplar to strong drought. Drought-stressed xylem was characterized by higher vessel frequencies, smaller vessel lumina, and thicker secondary fiber cell walls. These changes were accompanied by strong increases in abscisic acid (ABA) and antagonistic changes in salicylic acid in wood. Transcriptional evidence supported ABA biosynthesis and signaling in wood. Since ABA signaling activates the fiber-thickening factor NST1, we expected upregulation of the secondary cell wall (SCW) cascade under stress. By contrast, transcription factors and biosynthesis genes for SCW formation were down-regulated, whereas a small set of cellulose synthase-like genes and a huge array of genes involved in cell wall modification were up-regulated in drought-stressed wood. Therefore, we suggest that ABA signaling monitors normal SCW biosynthesis and that drought causes a switch from normal to "stress wood" formation recruiting a dedicated set of genes for cell wall biosynthesis and remodeling. This proposition implies that drought-induced changes in cell wall properties underlie regulatory mechanisms distinct from those of normal wood.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

PopulusPtERF85 Balances Xylem Cell Expansion and Secondary Cell Wall Formation in Hybrid Aspen

Secondary growth relies on precise and specialized transcriptional networks that determine cell division, differentiation, and maturation of xylem cells. We identified a novel role for the ethylene-induced Populus Ethylene Response Factor PtERF85 (Potri.015G023200) in balancing xylem cell expansion and secondary cell wall (SCW) formation in hybrid aspen (Populus tremula x tremuloides). Expression of PtERF85 is high in phloem and cambium cells and during the expansion of xylem cells, while it is low in maturing xylem tissue. Extending PtERF85 expression into SCW forming zones of woody tissues through ectopic expression reduced wood density and SCW thickness of xylem fibers but increased fiber diameter. Xylem transcriptomes from the transgenic trees revealed transcriptional induction of genes involved in cell expansion, translation, and growth. The expression of genes associated with plant vascular development and the biosynthesis of SCW chemical components such as xylan and lignin, was down-regulated in the transgenic trees. Our results suggest that PtERF85 activates genes related to xylem cell expansion, while preventing transcriptional activation of genes related to SCW formation. The importance of precise spatial expression of PtERF85 during wood development together with the observed phenotypes in response to ectopic PtERF85 expression suggests that PtERF85 contributes to the transition of fiber cells from elongation to secondary cell wall deposition.

Efficient method for isolation of high-quality RNA from Psidium guajava L. tissues

Acquiring high-quality RNA in sufficient amounts is crucial in plant molecular biology and genetic studies. Several methods for RNA extraction from plants are available in the literature, mainly due to the great biochemical diversity present in each species and tissue, which can complicate or prevent the extraction. Psidium guajava (Myrtaceae family) is a perennial fruit tree of medicinal and economic value; nevertheless, only a few molecular studies are available for the species. One reason is the difficulty in obtaining RNA due to the content of the samples, which are rich in polyphenols, polysaccharides, and secondary metabolites. Furthermore, there are few studies available for the isolation of RNA from guava or Psidium samples, which hampers advances in the study of the genus. Here, quality and yields of RNA isolates were compared using six extraction protocols: two protocols based on the application of cetyltrimethylammonium bromide (CTAB) lysis buffer, one protocol which uses the TRIzol reagent, one which applies guanidine thiocyanate lysis buffer followed by organic phase extraction, and two commercial kits (PureLink RNA Mini Kit and RNeasy Plant Mini Kit). The CTAB-based method provided the highest RNA yields and quality for five different tissues (flower bud, immature leaf, young leaf, mature leaf, and root), genotypes, and stress conditions. For the most efficient protocol, the average yield of RNA from guava leaves was 203.06 μg/g of tissue, and the A₂₆₀/A₂₈₀ and A₂₆₀/A₂₃₀ ratios were 2.1 and 2.2, respectively. RT-qPCR analysis demonstrated that the purity of the samples was sufficient for molecular biology experiments. CTAB-based methods for RNA isolation were found to be the most efficient, providing the highest RNA yields and quality for tissues from P. guajava. Additionally, they were compatible for downstream RNA-based applications, besides being simple and cost-effective.

Knockdown of miR393 Promotes the Growth and Biomass Production in Poplar

Short tandem target mimic (STTM), which is composed of two short sequences mimicking small RNA target sites, separated by a linker of optimal size, can block the functions of all members in a miRNA family. microRNA393 (miR393), which is one of the conserved miRNA families in plants, can regulate plant root growth, leaf development, plant architecture, and stress resistance. In order to verify the role of miR393 in the secondary growth of trees, we created its STTM transgenic poplar lines (STTM393). The expression of miR393 in STTM393 lines was reduced by over 10 times compared with the control plants. STTM393 lines showed promoted growth with about 20% higher, 15% thicker, and 2-4 more internodes than the control plants after 3 months of growth. The cross-section of the stems showed that STTM393 lines had wider phloem, xylem, and more cambium cell layers than control plants, and the lignin content in STTM393 lines was also higher as revealed by staining and chemical determination. Based on the transcriptome analysis, the genes related to the auxin signaling pathway, cell cyclin, cell expansion, and lignin synthesis had higher expression in STTM393 lines than that in control plants. The higher expression levels of FBL family members suggested that the auxin signaling pathway was strengthened in STTM393 lines to promote plant growth. Therefore, the knockdown of miR393 using the STTM approach provides a way to improve poplar growth and biomass production.

OMICs, Epigenetics, and Genome Editing Techniques for Food and Nutritional Security

The incredible success of crop breeding and agricultural innovation in the last century greatly contributed to the Green Revolution, which significantly increased yields and ensures food security, despite the population explosion. However, new challenges such as rapid climate change, deteriorating soil, and the accumulation of pollutants require much faster responses and more effective solutions that cannot be achieved through traditional breeding. Further prospects for increasing the efficiency of agriculture are undoubtedly associated with the inclusion in the breeding strategy of new knowledge obtained using high-throughput technologies and new tools in the future to ensure the design of new plant genomes and predict the desired phenotype. This article provides an overview of the current state of research in these areas, as well as the study of soil and plant microbiomes, and the prospective use of their potential in a new field of microbiome engineering. In terms of genomic and phenomic predictions, we also propose an integrated approach that combines high-density genotyping and high-throughput phenotyping techniques, which can improve the prediction accuracy of quantitative traits in crop species.

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.

High-spatiotemporal-resolution transcriptomes provide insights into fruit development and ripening in Citrus sinensis

Citrus fruit has a unique structure with soft leathery peel and pulp containing vascular bundles and several segments with many juice sacs. The function and morphology of each fruit tissue are different. Therefore, analysis at the organ-wide or mixed-tissue level inevitably obscures many tissue-specific phenomena. High-throughput RNA sequencing was used to profile Citrus sinensis fruit development based on four fruit tissue types and six development stages from young fruits to ripe fruits. Using a coexpression network analysis, modules of coexpressed genes and hub genes of tissue-specific networks were identified. Of particular, importance is the discovery of the regulatory network of phytohormones during citrus fruit development and ripening. A model was proposed to illustrate how ABF2 mediates the ABA signalling involved in sucrose transport, chlorophyll degradation, auxin homoeostasis, carotenoid and ABA biosynthesis, and cell wall metabolism during citrus fruit development. Moreover, we depicted the detailed spatiotemporal expression patterns of the genes involved in sucrose and citric acid metabolism in citrus fruit and identified several key genes that may play crucial roles in sucrose and citric acid accumulation in the juice sac, such as SWEET15 and CsPH8. The high spatial and temporal resolution of our data provides important insights into the molecular networks underlying citrus fruit development and ripening.

BEL1-like Homeodomain Protein BLH6a Is a Negative Regulator of CAl5H2 in Sinapyl Alcohol Monolignol Biosynthesis in Poplar

Lignin is one of the major components of xylem cell walls in tree stems. The lignin in the wood of most flowering plants (dicotyledonous angiosperms) is typically polymerized from three monolignol precursors, coniferyl alcohol, sinapyl alcohol, and p-coumaroyl alcohol, resulting in guaiacyl (G), syringyl (S), and hydroxyphenyl (H) subunits, respectively. In this study, we focus on the transcriptional regulation of a coniferaldehyde 5-hydroxylase (CAld5H2) gene, which encodes a key enzyme for sinapyl alcohol biosynthesis. We carried out a yeast one-hybrid (Y1H) screen to identify candidate upstream transcription factors (TFs) regulating CAld5H2. We obtained 12 upstream TFs as potential regulators of CAld5H2. One of these TF genes, BLH6a, encodes a BEL1-like homeodomain (BLH) protein and negatively regulated the CAld5H2 promoter activity. The direct regulation of CAld5H2 promoter by BLH6a was supported by chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) and dominant repression of BLH6a in transgenic plants. Luciferase complementation imaging analyses showed extensive protein-protein interactions among these 12 TFs. We propose that BLH6a is a negative regulator of CAld5H2, which acts through combinatorial regulation of multiple TFs for sinapyl alcohol (S monolignol) biosynthesis in poplar.

Xylem Parenchyma-Role and Relevance in Wood Functioning in Trees

Woody plants are characterised by a highly complex vascular system, wherein the secondary xylem (wood) is responsible for the axial transport of water and various substances. Previous studies have focused on the dead conductive elements in this heterogeneous tissue. However, the living xylem parenchyma cells, which constitute a significant functional fraction of the wood tissue, have been strongly neglected in studies on tree biology. Although there has recently been increased research interest in xylem parenchyma cells, the mechanisms that operate in these cells are poorly understood. Therefore, the present review focuses on selected roles of xylem parenchyma and its relevance in wood functioning. In addition, to elucidate the importance of xylem parenchyma, we have compiled evidence supporting the hypothesis on the significance of parenchyma cells in tree functioning and identified the key unaddressed questions in the field.

A Small Guanosine Triphosphate Binding Protein PagRabE1b Promotes Xylem Development in Poplar

Rab GTPases are the subfamily of the small guanosine triphosphate (GTP)-binding proteins which participated in the regulation of various biological processes. Recent studies have found that plant Rabs play some specific functions. However, the functions of Rabs in xylem development in trees remain unclear. In this study, functional identification of PagRabE1b in Populus was performed. Quantitative reverse transcription PCR (qRT-PCR) results showed that PagRabE1b was highly accumulated in stems, especially in phloem and xylem tissues. Overexpression of PagRabE1b in poplar enhanced programmed cell death (PCD) and increased the growth rate and the secondary cell wall (SCW) thickness. Quantitative analysis of monosaccharide content showed that various monosaccharides were significantly increased in secondary xylem tissues of the overexpressed lines. Flow cytometry analysis revealed that the number of apoptotic cells in PagRabE1b-OE lines is more than a wild type (WT), which indicated that PagRabE1b may play an important role in PCD. Further studies showed that overexpression of PagRabE1b increased the expression level of genes involved in SCW biosynthesis, PCD, and autophagy. Collectively, the results suggest that PagRabE1b plays a positive role in promoting the xylem development of poplar.

The BOP-type co-transcriptional regulator NODULE ROOT1 promotes stem secondary growth of the tropical Cannabaceae tree Parasponia andersonii

Tree stems undergo a massive secondary growth in which secondary xylem and phloem tissues arise from the vascular cambium. Vascular cambium activity is driven by endogenous developmental signalling cues and environmental stimuli. Current knowledge regarding the genetic regulation of cambium activity and secondary growth is still far from complete. The tropical Cannabaceae tree Parasponia andersonii is a non-legume research model of nitrogen-fixing root nodulation. Parasponia andersonii can be transformed efficiently, making it amenable for CRISPR-Cas9-mediated reverse genetics. We considered whether P. andersonii also could be used as a complementary research system to investigate tree-related traits, including secondary growth. We established a developmental map of stem secondary growth in P. andersonii plantlets. Subsequently, we showed that the expression of the co-transcriptional regulator PanNODULE ROOT1 (PanNOOT1) is essential for controlling this process. PanNOOT1 is orthologous to Arabidopsis thaliana BLADE-ON-PETIOLE1 (AtBOP1) and AtBOP2, which are involved in the meristem-to-organ-boundary maintenance. Moreover, in species forming nitrogen-fixing root nodules, NOOT1 is known to function as a key nodule identity gene. Parasponia andersonii CRISPR-Cas9 loss-of-function Pannoot1 mutants are altered in the development of the xylem and phloem tissues without apparent disturbance of the cambium organization and size. Transcriptomic analysis showed that the expression of key secondary growth-related genes is significantly down-regulated in Pannoot1 mutants. This allows us to conclude that PanNOOT1 positively contributes to the regulation of stem secondary growth. Our work also demonstrates that P. andersonii can serve as a tree research system.

Synonymous mutation in Growth Regulating Factor 15 of miR396a target sites enhances photosynthetic efficiency and heat tolerance in poplar

Heat stress damages plant tissues and induces multiple adaptive responses. Complex and spatiotemporally specific interactions among transcription factors (TFs), microRNAs (miRNAs), and their targets play crucial roles in regulating stress responses. To explore these interactions and to identify regulatory networks in perennial woody plants subjected to heat stress, we integrated time-course RNA-seq, small RNA-seq, degradome sequencing, weighted gene correlation network analysis, and multi-gene association approaches in poplar. Results from Populus trichocarpa enabled us to construct a three-layer, highly interwoven regulatory network involving 15 TFs, 45 miRNAs, and 77 photosynthetic genes. Candidate gene association studies in a population of P. tomentosa identified 114 significant associations and 696 epistatic SNP–SNP pairs that were linked to 29 photosynthetic and growth traits (P<0.0001, q<0.05). We also identified miR396a and its target, Growth-Regulating Factor 15 (GRF15) as an important regulatory module in the heat-stress response. Transgenic plants of hybrid poplar (P. alba × P. glandulosa) overexpressing a GRF15 mRNA lacking the miR396a target sites exhibited enhanced heat tolerance and photosynthetic efficiency compared to wild-type plants. Together, our observations demonstrate that GRF15 plays a crucial role in responding to heat stress, and they highlight the power of this new, multifaceted approach for identifying regulatory nodes in plants.

Expression and localization of SWEETs in Populus and the effect of SWEET7 overexpression in secondary growth

In trees, wood formation needs carbon import from the photosynthetic source tissues. Sugar transporters play important roles in carbohydrate transport into wood-forming cells. Sugars will eventually be exported transporters (SWEETs) play essential roles in many physiological processes. However, the roles of this family in the growth and development of woody plants have not been systematically investigated. In this study, 27 SWEET genes were identified in the Populus trichocarpa genome. These SWEET genes were classified into four clades based on their phylogenetic relationships, gene structures, conserved motifs and chromosomal locations. Representative SWEET members from each clade were selected for further studies. The PagSWEETs of Populus alba × Populus glandulosa were localized to plasma membrane, vacuolar, endoplasmic reticulum or Golgi. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that PagSWEETs have distinct expression patterns in various tissues, and PagSWEET5, 7, 10b, 10c, 15b, 17a and 17c exhibited high expression levels in stems. PagSWEET7 is localized to the cytoplasmic membrane and specifically expressed in the phloem as detected by histochemical GUS ($\beta $ - glucuronidase) assays. Xylem production and xylem sugar content were greater in developing wood of PagSWEET7 overexpression than wild-type lines. Collectively, these results provide valuable information for further investigating functions of PagSWEET genes, and identify PagSWEET7 as a candidate gene for using biotechnology to modify the wood formation in poplar.

Cytokinin signaling localized in phloem noncell-autonomously regulates cambial activity during secondary growth of Populus stems

The regulation of cytokinin on secondary vascular development has been uncovered by modulating cytokinin content. However, it remains unclear how cytokinin enriched in developing secondary phloem regulates cambium activity in poplar. Here, we visualized the gradient distribution of cytokinin with a peak in the secondary phloem of poplar stem via immunohistochemical imaging, and determined the role of phloem-located cytokinin signaling during wood formation. We generated transgenic poplar harboring cytokinin oxidase/dehydrogenase (CKX)2, a gene encoding a cytokinin degrading enzyme, driven by the phloem-specific CLE41b promoter, indicating that the disruption of the cytokinin gradient pattern restricts the cambial activity. The RNA interference-based knockdown of the histidine kinase (HK) genes encoding cytokinin receptors specifically in secondary phloem significantly compromised the division activity of cambial cells, whereas the phloem-specific expression of a type-B response regulator (RR) transcription factor stimulated cambial proliferation, providing evidence for the noncell-autonomous regulation of local cytokinin signaling on the cambial activity. Moreover, the cambium-specific knockdown of HKs also led to restricted cambial activity, and the defects were aggravated by the reduced cytokinin accumulation. Our results showed that local cytokinin signaling in secondary phloem regulates cambial activity noncell-autonomously, and coordinately with its local signaling in cambium.

Tissue-specific transcriptome profiling of the Arabidopsis inflorescence stem reveals local cellular signatures

Genome-wide gene expression maps with a high spatial resolution have substantially accelerated plant molecular science. However, the number of characterized tissues and growth stages is still small due to the limited accessibility of most tissues for protoplast isolation. Here, we provide gene expression profiles of the mature inflorescence stem of Arabidopsis thaliana covering a comprehensive set of distinct tissues. By combining fluorescence-activated nucleus sorting and laser-capture microdissection with next-generation RNA sequencing, we characterized the transcriptomes of xylem vessels, fibers, the proximal and distal cambium, phloem, phloem cap, pith, starch sheath, and epidermis cells. Our analyses classified more than 15,000 genes as being differentially expressed among different stem tissues and revealed known and novel tissue-specific cellular signatures. By determining overrepresented transcription factor binding regions in the promoters of differentially expressed genes, we identified candidate tissue-specific transcriptional regulators. Our datasets predict the expression profiles of an exceptional number of genes and allow hypotheses to be generated about the spatial organization of physiological processes. Moreover, we demonstrate that information about gene expression in a broad range of mature plant tissues can be established at high spatial resolution by nuclear mRNA profiling. Tissue-specific gene expression values can be accessed online at https://arabidopsis-stem.cos.uni-heidelberg.de/.

Transcriptomic Analysis of Seasonal Gene Expression and Regulation during Xylem Development in “Shanxin” Hybrid Poplar (Populus davidiana × Populus bolleana)

Xylem development is a key process for wood formation in woody plants. To study the molecular regulatory mechanisms related to xylem development in hybrid poplar P. davidiana × P. bolleana, transcriptome analyses were conducted on developing xylem at six different growth stages within a single growing season. Xylem development and differentially expressed genes in the six time points were selected for a regulatory analysis. Xylem development was observed in stem sections at different growth stages, which showed that xylem development extended from the middle of April to early August and included cell expansion and secondary cell wall biosynthesis. An RNA-seq analysis of six samples with three replicates was performed. After transcriptome assembly and annotation, the differentially expressed genes (DEGs) were identified, and a Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and expression analysis of the DEGs were performed on each sample. On average, we obtained >20 million clean reads per sample, which were assembled into 84,733 nonredundant transcripts, of which there were 17,603 unigenes with lengths >1 kb. There were 14,890 genes that were differentially expressed among the six stages. The upregulated DEGs were enriched in GO terms related to cell wall biosynthesis between S1 vs. S2 or S3 vs. S4 and, in GO terms, related to phytohormones in the S1 vs. S2 or S4 vs. S5 comparisons. The downregulated DEGs were enriched in GO terms related to cell wall biosynthesis between S4 vs. S5 or S5 vs. S6 and, in GO terms, related to hormones between S1 vs. S2 or S2 vs. S3. The KEGG pathways in the DEGs related to “phenylpropanoid biosynthesis”, “plant hormone signal transduction” and “starch and sucrose metabolism” were significantly enriched among the different stages. The DEGs related to cell expansion, polysaccharide metabolism and synthesis, lignin synthesis, transcription factors and hormones were identified. The identification of genes involved in the regulation of xylem development will increase our understanding of the molecular regulation of wood formation in trees and, also, offers potential targets for genetic manipulation to improve the properties of wood.

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.

Genome-Wide Identification of Populus Malectin/Malectin-Like Domain-Containing Proteins and Expression Analyses Reveal Novel Candidates for Signaling and Regulation of Wood Development

Malectin domain (MD) is a ligand-binding protein motif of pro- and eukaryotes. It is particularly abundant in Viridiplantae, where it occurs as either a single (MD, PF11721) or tandemly duplicated domain (PF12819) called malectin-like domain (MLD). In herbaceous plants, MD- or MLD-containing proteins (MD proteins) are known to regulate development, reproduction, and resistance to various stresses. However, their functions in woody plants have not yet been studied. To unravel their potential role in wood development, we carried out genome-wide identification of MD proteins in the model tree species black cottonwood (Populus trichocarpa), and analyzed their expression and co-expression networks. P. trichocarpa had 146 MD genes assigned to 14 different clades, two of which were specific to the genus Populus. 87% of these genes were located on chromosomes, the rest being associated with scaffolds. Based on their protein domain organization, and in agreement with the exon-intron structures, the MD genes identified here could be classified into five superclades having the following domains: leucine-rich repeat (LRR)-MD-protein kinase (PK), MLD-LRR-PK, MLD-PK (CrRLK1L), MLD-LRR, and MD-Kinesin. Whereas the majority of MD genes were highly expressed in leaves, particularly under stress conditions, eighteen showed a peak of expression during secondary wall formation in the xylem and their co-expression networks suggested signaling functions in cell wall integrity, pathogen-associated molecular patterns, calcium, ROS, and hormone pathways. Thus, P. trichocarpa MD genes having different domain organizations comprise many genes with putative foliar defense functions, some of which could be specific to Populus and related species, as well as genes with potential involvement in signaling pathways in other tissues including developing wood.

Identifying transcription factors that reduce wood recalcitrance and improve enzymatic degradation of xylem cell wall in Populus

Developing an efficient deconstruction step of woody biomass for biorefinery has been drawing considerable attention since its xylem cell walls display highly recalcitrance nature. Here, we explored transcriptional factors (TFs) that reduce wood recalcitrance and improve saccharification efficiency in Populus species. First, 33 TF genes up-regulated during poplar wood formation were selected as potential regulators of xylem cell wall structure. The transgenic hybrid aspens (Populus tremula × Populus tremuloides) overexpressing each selected TF gene were screened for in vitro enzymatic saccharification. Of these, four transgenic seedlings overexpressing previously uncharacterized TF genes increased total glucan hydrolysis on average compared to control. The best performing lines overexpressing Pt × tERF123 and Pt × tZHD14 were further grown to form mature xylem in the greenhouse. Notably, the xylem cell walls exhibited significantly increased total xylan hydrolysis as well as initial hydrolysis rates of glucan. The increased saccharification of Pt × tERF123-overexpressing lines could reflect the improved balance of cell wall components, i.e., high cellulose and low xylan and lignin content, which could be caused by upregulation of cellulose synthase genes upon the expression of Pt × tERF123. Overall, we successfully identified Pt × tERF123 and Pt × tZHD14 as effective targets for reducing cell wall recalcitrance and improving the enzymatic degradation of woody plant biomass.

The Spatio-Temporal Distribution of Cell Wall-Associated Glycoproteins During Wood Formation in Populus

Plant cell wall associated hydroxyproline-rich glycoproteins (HRGPs) are involved in several aspects of plant growth and development, including wood formation in trees. HRGPs such as arabinogalactan-proteins (AGPs), extensins (EXTs), and proline rich proteins (PRPs) are important for the development and architecture of plant cell walls. Analysis of publicly available gene expression data revealed that many HRGP encoding genes show tight spatio-temporal expression patterns in the developing wood of Populus that are indicative of specific functions during wood formation. Similar results were obtained for the expression of glycosyl transferases putatively involved in HRGP glycosylation. In situ immunolabelling of transverse wood sections using AGP and EXT antibodies revealed the cell type specificity of different epitopes. In mature wood AGP epitopes were located in xylem ray cell walls, whereas EXT epitopes were specifically observed between neighboring xylem vessels, and on the ray cell side of the vessel walls, likely in association with pits. Molecular mass and glycan analysis of AGPs and EXTs in phloem/cambium, developing xylem, and mature xylem revealed clear differences in glycan structures and size between the tissues. Separation of AGPs by agarose gel electrophoresis and staining with β-D-glucosyl Yariv confirmed the presence of different AGP populations in phloem/cambium and xylem. These results reveal the diverse changes in HRGP-related processes that occur during wood formation at the gene expression and HRGP glycan biosynthesis levels, and relate HRGPs and glycosylation processes to the developmental processes of wood formation.

A metabolite roadmap of the wood-forming tissue in Populus tremula

Wood, or secondary xylem, is the product of xylogenesis, a developmental process that begins with the proliferation of cambial derivatives and ends with mature xylem fibers and vessels with lignified secondary cell walls. Fully mature xylem has undergone a series of cellular processes, including cell division, cell expansion, secondary wall formation, lignification and programmed cell death. A complex network of interactions between transcriptional regulators and signal transduction pathways controls wood formation. However, the role of metabolites during this developmental process has not been comprehensively characterized. To evaluate the role of metabolites during wood formation, we performed a high spatial resolution metabolomics study of the wood-forming zone of Populus tremula, including laser dissected aspen ray and fiber cells. We show that metabolites show specific patterns within the wood-forming zone, following the differentiation process from cell division to cell death. The data from profiled laser dissected aspen ray and fiber cells suggests that these two cell types host distinctly different metabolic processes. Furthermore, by integrating previously published transcriptomic and proteomic profiles generated from the same trees, we provide an integrative picture of molecular processes, for example, deamination of phenylalanine during lignification is of critical importance for nitrogen metabolism during wood formation.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

ACAULIS5 Is Required for Cytokinin Accumulation and Function During Secondary Growth of Populus Trees

In the primary root and young hypocotyl of Arabidopsis, ACAULIS5 promotes translation of SUPPRESSOR OF ACAULIS51 (SAC51) and thereby inhibits cytokinin biosynthesis and vascular cell division. In this study, the relationships between ACAULIS5, SAC51 and cytokinin biosynthesis were investigated during secondary growth of Populus stems. Overexpression of ACAULIS5 from the constitutive 35S promoter in hybrid aspen (Populus tremula × Populus tremuloides) trees suppressed the expression level of ACAULIS5, which resulted in low levels of the physiologically active cytokinin bases as well as their direct riboside precursors in the transgenic lines. Low ACAULIS5 expression and low cytokinin levels of the transgenic trees coincided with low cambial activity of the stem. ACAULIS5 therefore, contrary to its function in young seedlings in Arabidopsis, stimulates cytokinin accumulation and cambial activity during secondary growth of the stem. This function is not derived from maturing secondary xylem tissues as transgenic suppression of ACAULIS5 levels in these tissues did not influence secondary growth. Interestingly, evidence was obtained for increased activity of the anticlinal division of the cambial initials under conditions of low ACAULIS5 expression and low cytokinin accumulation. We propose that ACAULIS5 integrates auxin and cytokinin signaling to promote extensive secondary growth of tree stems.

MYB Transcription Factor161 Mediates Feedback Regulation of Secondary wall-associated NAC-Domain1 Family Genes for Wood Formation

Wood formation is a complex process that involves cell differentiation, cell expansion, secondary wall deposition, and programmed cell death. We constructed a four-layer wood formation transcriptional regulatory network (TRN) in Populus trichocarpa (black cottonwood) that has four Secondary wall-associated NAC-Domain1 (PtrSND1) transcription factor (TF) family members as the top-layer regulators. We characterized the function of a MYB (PtrMYB161) TF in this PtrSND1-TRN, using transgenic P trichocarpa cells and whole plants. PtrMYB161 is a third-layer regulator that directly transactivates five wood formation genes. Overexpression of PtrMYB161 in P. trichocarpa (OE-PtrMYB161) led to reduced wood, altered cell type proportions, and inhibited growth. Integrative analysis of wood cell-based chromatin-binding assays with OE-PtrMYB161 transcriptomics revealed a feedback regulation system in the PtrSND1-TRN, where PtrMYB161 represses all four top-layer regulators and one second-layer regulator, PtrMYB021, possibly affecting many downstream TFs in, and likely beyond, the TRN, to generate the observed phenotypic changes. Our data also suggested that the PtrMYB161's repressor function operates through interaction of the base PtrMYB161 target-binding system with gene-silencing cofactors. PtrMYB161 protein does not contain any known negative regulatory domains. CRISPR-based mutants of PtrMYB161 in P. trichocarpa exhibited phenotypes similar to the wild type, suggesting that PtrMYB161's activator functions are redundant among many TFs. Our work demonstrated that PtrMYB161 binds to multiple sets of target genes, a feature that allows it to function as an activator as well as a repressor. The balance of the two functions may be important to the establishment of regulatory homeostasis for normal growth and development.

MiR319a-targeted PtoTCP20 regulates secondary growth via interactions with PtoWOX4 and PtoWND6 in Populus tomentosa

Secondary growth is a key characteristic of trees, which requires the coordination of multiple regulatory mechanisms including transcriptional regulators and microRNAs (miRNAs). However, the roles of microRNAs in the regulation of secondary growth need to be explored in depth. Here, the role of miR319a and its target, PtoTCP20, in the secondary growth of Populus tomentosa stem was investigated using genetic and molecular analyses. The expression level of miR319a gradually decreased from primary to secondary growth in P. tomentosa, while that of PtoTCP20 gradually increased. MiR319a overexpression in seedlings resulted in delayed secondary growth and decreased xylem production, while miR319a knockdown and PtoTCP20 overexpression promoted secondary growth and increased xylem production. Further analysis showed that PtoTCP20 interacted with PtoWOX4a and activated PtoWND6 transcription in vitro and in vivo. Our data show that PtoTCP20 controls vascular cambium proliferation by binding to PtoWOX4a, and promotes secondary xylem differentiation by activating PtoWND6 transcription, thereby regulating secondary growth in P. tomentosa. Our findings provide insight into the molecular mechanisms underlying secondary growth in trees.

Leaf shape in Populus tremula is a complex, omnigenic trait

Leaf shape is a defining feature of how we recognize and classify plant species. Although there is extensive variation in leaf shape within many species, few studies have disentangled the underlying genetic architecture. We characterized the genetic architecture of leaf shape variation in Eurasian aspen (Populus tremula L.) by performing genome-wide association study (GWAS) for physiognomy traits. To ascertain the roles of identified GWAS candidate genes within the leaf development transcriptional program, we generated RNA-Seq data that we used to perform gene co-expression network analyses from a developmental series, which is publicly available within the PlantGenIE resource. We additionally used existing gene expression measurements across the population to analyze GWAS candidate genes in the context of a population-wide co-expression network and to identify genes that were differentially expressed between groups of individuals with contrasting leaf shapes. These data were integrated with expression GWAS (eQTL) results to define a set of candidate genes associated with leaf shape variation. Our results identified no clear adaptive link to leaf shape variation and indicate that leaf shape traits are genetically complex, likely determined by numerous small-effect variations in gene expression. Genes associated with shape variation were peripheral within the population-wide co-expression network, were not highly connected within the leaf development co-expression network, and exhibited signatures of relaxed selection. As such, our results are consistent with the omnigenic model.

Genetic control of tracheid properties in Norway spruce wood

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

Characterization of the UDP-glycosyltransferase UGT72 Family in Poplar and Identification of Genes Involved in the Glycosylation of Monolignols

Monolignols are the building blocks for lignin polymerization in the apoplastic domain. Monolignol biosynthesis, transport, storage, glycosylation, and deglycosylation are the main biological processes partaking in their homeostasis. In Arabidopsis thaliana, members of the uridine diphosphate-dependent glucosyltransferases UGT72E and UGT72B subfamilies have been demonstrated to glycosylate monolignols. Here, the poplar UGT72 family, which is clustered into four groups, was characterized: Group 1 UGT72AZ1 and UGT72AZ2, homologs of Arabidopsis UGT72E1-3, as well as group 4 UGT72B37 and UGT72B39, homologs of Arabidopsis UGT72B1-3, glycosylate monolignols. In addition, promoter-GUS analyses indicated that poplar UGT72 members are expressed within vascular tissues. At the subcellular level, poplar UGT72s belonging to group 1 and group 4 were found to be associated with the nucleus and the endoplasmic reticulum. However, UGT72A2, belonging to group 2, was localized in bodies associated with chloroplasts, as well as possibly in chloroplasts. These results show a partial conservation of substrate recognition between Arabidopsis and poplar homologs, as well as divergent functions between different groups of the UGT72 family, for which the substrates remain unknown.

Network aggregation improves gene function prediction of grapevine gene co-expression networks

Aggregation across multiple networks highlights robust co-expression interactions and improves the functional connectivity of grapevine gene co-expression networks. In recent years, the rapid accumulation of transcriptome datasets from diverse experimental conditions has enabled the widespread use of gene co-expression network (GCN) analysis in plants. In grapevine, GCN analysis has shown great promise for gene function prediction, however, measurable progress is currently lacking. Using accumulated microarray datasets from the grapevine whole-genome array (33 experiments, 1359 samples), we explored how meta-analysis through aggregation influences the functional connectivity (performance) of derived networks using guilt-by-association neighbor voting. Two annotation schemes, i.e. MapMan BIN and Pfam, at two sparsity thresholds, i.e. top 100 (stringent) and 300 (relaxed) ranked genes were evaluated. We observed that aggregating across multiple networks improves performance dramatically, with the aggregate outperforming the majority of functional terms across individual networks. Network sparsity and size (i.e. the number of samples and aggregates) were key factors influencing performance while the choice of annotation scheme had little. Systematic comparison with various state-of-the-art microarray and RNA-seq networks was also performed, however, none outperformed the aggregate microarray network despite having good predictive performance. Repeating these series of tests using a functional enrichment-based performance metric also showed remarkably consistent findings with guilt-by-association neighbor voting. To demonstrate its functionality, we explore the function and transcriptional regulation of grapevine EXPANSIN genes. We envisage that network aggregation will offer new and unique opportunities for gene function prediction in future grapevine functional genomics studies. To this end, we make the aggregate networks and associated metadata publicly available at VTC-Agg (https://sites.google.com/view/vtc-agg).

Overexpression of a developing xylem cDNA library in transgenic poplar generates high mutation rate specific to wood formation

We investigated feasibility of the Full-length complementary DNA OvereXpression (FOX) system as a mutagenesis approach in poplar, using developing xylem tissue. The main goal was to assess the overall mutation rate and if the system will increase instances of mutants affected in traits linked to the xylem tissue. Indeed, we found a high mutation rate of 17.7%, whereas 80% of all mutants were significantly affected in cellulose, lignin and/or hemicellulose. Cell wall biosynthesis is a major process occurring during xylem development. Enrichment of mutants affected in cell wall composition suggests that the tissue source for the FOX library influenced the occurrence of mutants affected in a trait linked to this tissue. Additionally, we found that FLcDNAs from mutants affected in cell wall composition were homologous to genes known to be involved in cell wall biosynthesis and most recovered FLcDNAs corresponded to genes whose native expression was highest in xylem. We characterized in detail a mutant line with increased diameter. The phenotype was caused by a poplar homolog of LONELY GUY 1 (LOG1), which encodes an enzyme in cytokinin biosynthesis and significantly increased xylem proliferation. The causative role of LOG1 in the observed phenotype was further reaffirmed by elevated cytokinin concentration in the mutant and recapitulation overexpression experiment wherein multiple independent lines phenocopied the original FOX mutant. Our experiments show that the FOX approach can be efficiently used for gene discovery and molecular interrogation of traits specific to woody perennial growth and development.