Merit of integrating in situ transcriptomics and anatomical information for cell annotation and lineage construction in single-cell analyses of Populus.

Cell type annotation and lineage construction are two of the most critical tasks conducted in the analyses of single-cell RNA sequencing (scRNA-seq). Four recent scRNA-seq studies of differentiating xylem propose four models on differentiating xylem development in Populus. The differences are mostly caused by the use of different strategies for cell type annotation and subsequent lineage interpretation. Here, we emphasize the necessity of using in situ transcriptomes and anatomical information to construct the most plausible xylem development model.

An AGO10:miR165/6 module regulates meristem activity and xylem development in the Arabidopsis root.

The RNA-silencing effector ARGONAUTE10 influences cell fate in plant shoot and floral meristems. ARGONAUTE10 also accumulates in the root apical meristem (RAM), yet its function(s) therein remain elusive. Here, we show that ARGONAUTE10 is expressed in the root cell initials where it controls overall RAM activity and length. ARGONAUTE10 is also expressed in the stele, where post-transcriptional regulation confines it to the root tip's pro-vascular region. There, variations in ARGONAUTE10 levels modulate metaxylem-vs-protoxylem specification. Both ARGONAUTE10 functions entail its selective, high-affinity binding to mobile miR165/166 transcribed in the neighboring endodermis. ARGONAUTE10-bound miR165/166 is degraded, likely via SMALL-RNA-DEGRADING-NUCLEASES1/2, thus reducing miR165/166 ability to silence, via ARGONAUTE1, the transcripts of cell fate-influencing transcription factors. These include PHABULOSA (PHB), which controls meristem activity in the initials and xylem differentiation in the pro-vasculature. During early germination, PHB transcription increases while dynamic, spatially-restricted transcriptional and post-transcriptional mechanisms reduce and confine ARGONAUTE10 accumulation to the provascular cells surrounding the newly-forming xylem axis. Adequate miR165/166 concentrations are thereby channeled along the ARGONAUTE10-deficient yet ARGONAUTE1-proficient axis. Consequently, inversely-correlated miR165/166 and PHB gradients form preferentially along the axis despite ubiquitous PHB transcription and widespread miR165/166 delivery inside the whole vascular cylinder.

ARGONAUTE10 controls cell fate specification and formative cell divisions in the Arabidopsis root.

A key question in plant biology is how oriented cell divisions are integrated with patterning mechanisms to generate organs with adequate cell type allocation. In the root vasculature, a gradient of miRNA165/6 controls the abundance of HD-ZIP III transcription factors, which in turn control cell fate and spatially restrict vascular cell proliferation to specific cells. Here, we show that vascular development requires the presence of ARGONAUTE10, which is thought to sequester miRNA165/6 and protect HD-ZIP III transcripts from degradation. Our results suggest that the miR165/6-AGO10-HDZIP III module acts by buffering cytokinin responses and restricting xylem differentiation. Mutants of AGO10 show faster growth rates and strongly enhanced survival under severe drought conditions. However, this superior performance is offset by markedly increased variation and phenotypic plasticity in sub-optimal carbon supply conditions. Thus, AGO10 is required for the control of formative cell division and coordination of robust cell fate specification of the vasculature, while altering its expression provides a means to adjust phenotypic plasticity.

Blue light photoreceptor cryptochrome 1 promotes wood formation and anthocyanin biosynthesis in Populus.

Blue light photoreceptor cryptochrome 1 (CRY1) in herbaceous plants plays crucial roles in various developmental processes, including cotyledon expansion, hypocotyl elongation and anthocyanin biosynthesis. However, the function of CRY1 in perennial trees is unclear. In this study, we identified two ortholog genes of CRY1 (PagCRY1a and PagCRY1b) from Populus, which displayed high sequence similarity to Arabidopsis CRY1. Overexpression of PagCRY1 substantially inhibited plant growth and promoted secondary xylem development in Populus, while CRISPR/Cas9-mediated knockout of PagCRY1 enhanced plant growth and delayed secondary xylem development. Moreover, overexpression of PagCRY1 dramatically increased anthocyanin accumulation. The further analysis supported that PagCRY1 functions specifically in response to blue light. Taken together, our results demonstrated that modulating the expression of blue light photoreceptor CRY1 ortholog gene in Populus could significantly influence plant biomass production and the process of wood formation, laying a foundation for further investigating the light-regulated tree growth.

Transcription factor PagMYB31 positively regulates cambium activity and negatively regulates xylem development in poplar.

Wood formation involves consecutive developmental steps, including cell division of vascular cambium, xylem cell expansion, secondary cell wall (SCW) deposition, and programmed cell death. In this study, we identified PagMYB31 as a coordinator regulating these processes in Populus alba × Populus glandulosa and built a PagMYB31-mediated transcriptional regulatory network. PagMYB31 mutation caused fewer layers of cambial cells, larger fusiform initials, ray initials, vessels, fiber and ray cells, and enhanced xylem cell SCW thickening, showing that PagMYB31 positively regulates cambial cell proliferation and negatively regulates xylem cell expansion and SCW biosynthesis. PagMYB31 repressed xylem cell expansion and SCW thickening through directly inhibiting wall-modifying enzyme genes and the transcription factor genes that activate the whole SCW biosynthetic program, respectively. In cambium, PagMYB31 could promote cambial activity through TRACHEARY ELEMENT DIFFERENTIATION INHIBITORY FACTOR (TDIF)/PHLOEM INTERCALATED WITH XYLEM (PXY) signaling by directly regulating CLAVATA3/ESR-RELATED (CLE) genes, and it could also directly activate WUSCHEL HOMEOBOX RELATED4 (PagWOX4), forming a feedforward regulation. We also observed that PagMYB31 could either promote cell proliferation through the MYB31-MYB72-WOX4 module or inhibit cambial activity through the MYB31-MYB72-VASCULAR CAMBIUM-RELATED MADS2 (VCM2)/PIN-FORMED5 (PIN5) modules, suggesting its role in maintaining the homeostasis of vascular cambium. PagMYB31 could be a potential target to manipulate different developmental stages of wood formation.

Analyses of high spatial resolution datasets identify genes associated with multi-layered secondary cell wall thickening in Pinus bungeana.

BACKGROUND AND AIMS: Secondary cell wall (SCW) thickening is a major cellular developmental stage determining wood structure and properties. Although the molecular regulation of cell wall deposition during tracheary element differentiation has been well established in primary growth systems, less is known about the gene regulatory processes involved in the multi-layered SCW thickening of mature trees. METHODS: Using third-generation [long-read single-molecule real-time (SMRT)] and second-generation [short-read sequencing by synthesis (SBS)] sequencing methods, we established a Pinus bungeana transcriptome resource with comprehensive functional and structural annotation for the first time. Using these approaches, we generated high spatial resolution datasets for the vascular cambium, xylem expansion regions, early SCW thickening, late SCW thickening and mature xylem tissues of 71-year-old Pinus bungeana trees. KEY RESULTS: A total of 79 390 non-redundant transcripts, 31 808 long non-coding RNAs and 5147 transcription factors were annotated and quantified in different xylem tissues at all growth and differentiation stages. Furthermore, using this high spatial resolution dataset, we established a comprehensive transcriptomic profile and found that members of the NAC, WRKY, SUS, CESA and LAC gene families are major players in early SCW formation in tracheids, whereas members of the MYB and LBD transcription factor families are highly expressed during late SCW thickening. CONCLUSIONS: Our results provide new molecular insights into the regulation of multi-layered SCW thickening in conifers. The high spatial resolution datasets provided can serve as important gene resources for improving softwoods.

Ectopic assembly of an auxin efflux control machinery shifts developmental trajectories.

Polar auxin transport in the Arabidopsis (Arabidopsis thaliana) root tip maintains high auxin levels around the stem cell niche that gradually decrease in dividing cells but increase again once they transition toward differentiation. Protophloem differentiates earlier than other proximal tissues and employs a unique auxin "canalization" machinery that is thought to balance auxin efflux with retention. It consists of a proposed activator of PIN-FORMED (PIN) auxin efflux carriers, the cAMP-, cGMP- and Calcium-dependent (AGC) kinase PROTEIN KINASE ASSOCIATED WITH BRX (PAX); its inhibitor, BREVIS RADIX (BRX); and PHOSPHATIDYLINOSITOL-4-PHOSPHATE-5-KINASE (PIP5K) enzymes, which promote polar PAX and BRX localization. Because of a dynamic PAX-BRX-PIP5K interplay, the net cellular output of this machinery remains unclear. In this study, we deciphered the dosage-sensitive regulatory interactions among PAX, BRX, and PIP5K by their ectopic expression in developing xylem vessels. The data suggest that the dominant collective output of the PAX-BRX-PIP5K module is a localized reduction in PIN abundance. This requires PAX-stimulated clathrin-mediated PIN endocytosis upon site-specific phosphorylation, which distinguishes PAX from other AGC kinases. An ectopic assembly of the PAX-BRX-PIP5K module is sufficient to cause cellular auxin retention and affects root growth vigor by accelerating the trajectory of xylem vessel development. Our data thus provide direct evidence that local manipulation of auxin efflux alters the timing of cellular differentiation in the root.

Evaluation of reference genes and characterization of the MYBs in xylem radial change of Chinese fir stem.

The radial change (RC) of tree stem is the process of heartwood formation involved in complex molecular mechanism. Chinese fir (Cunninghamia lanceolata (Lamb.) Hook.), an evergreen species, is an important fast-growing timber tree in southern China. In this study, the top four stable genes (IDH, UBC2, RCA and H2B) were selected in RC tissues of 15 years old Chinese fir stem (RC15) and the genes (H2B, 18S, TIP41 and GAPDH) were selected in RC tissues of 30 years old Chinese fir stem (RC30). The stability of the reference genes is higher in RC30 than in RC15. Sixty-one MYB transcripts were obtained on the PacBio Sequel platform from woody tissues of one 30 years old Chinese fir stem. Based on the number of MYB DNA-binding domain and phylogenetic relationships, the ClMYB transcripts contained 21 transcripts of MYB-related proteins (1R-MYB), 39 transcripts of R2R3-MYB proteins (2R-MYB), one transcript of R1R2R3-MYB protein (3R-MYB) belonged to 18 function-annotated clades and two function-unknown clades. In RC woody tissues of 30 years old Chinese fir stem, ClMYB22 was the transcript with the greatest fold change detected by both RNA-seq and qRT-PCR. Reference genes selected in this study will be helpful for further verification of transcript abundance patterns during the heartwood formation of Chinese fir.

A multi-organ metabolic model of tomato predicts plant responses to nutritional and genetic perturbations.

Predicting and understanding plant responses to perturbations require integrating the interactions between nutritional sources, genes, cell metabolism, and physiology in the same model. This can be achieved using metabolic modeling calibrated by experimental data. In this study, we developed a multi-organ metabolic model of a tomato (Solanum lycopersicum) plant during vegetative growth, named Virtual Young TOmato Plant (VYTOP) that combines genome-scale metabolic models of leaf, stem and root and integrates experimental data acquired from metabolomics and high-throughput phenotyping of tomato plants. It is composed of 6,689 reactions and 6,326 metabolites. We validated VYTOP predictions on five independent use cases. The model correctly predicted that glutamine is the main organic nutrient of xylem sap. The model estimated quantitatively how stem photosynthetic contribution impacts exchanges between the different organs. The model was also able to predict how nitrogen limitation affects vegetative growth and the metabolic behavior of transgenic tomato lines with altered expression of core metabolic enzymes. The integration of different components, such as a metabolic model, physiological constraints, and experimental data, generates a powerful predictive tool to study plant behavior, which will be useful for several other applications, such as plant metabolic engineering or plant nutrition.

A vacuolar hexose transport is required for xylem development in the inflorescence stem.

In Angiosperms, the development of the vascular system is controlled by a complex network of transcription factors. However, how nutrient availability in the vascular cells affects their development remains to be addressed. At the cellular level, cytosolic sugar availability is regulated mainly by sugar exchanges at the tonoplast through active and/or facilitated transport. In Arabidopsis (Arabidopsis thaliana), among the genes encoding tonoplastic transporters, SUGAR WILL EVENTUALLY BE EXPORTED TRANSPORTER 16 (SWEET16) and SWEET17 expression has been previously detected in the vascular system. Here, using a reverse genetics approach, we propose that sugar exchanges at the tonoplast, regulated by SWEET16, are important for xylem cell division as revealed in particular by the decreased number of xylem cells in the swt16 mutant and the accumulation of SWEET16 at the procambium-xylem boundary. In addition, we demonstrate that transport of hexoses mediated by SWEET16 and/or SWEET17 is required to sustain the formation of the xylem secondary cell wall. This result is in line with a defect in the xylem cell wall composition as measured by Fourier-transformed infrared spectroscopy in the swt16swt17 double mutant and by upregulation of several genes involved in secondary cell wall synthesis. Our work therefore supports a model in which xylem development partially depends on the exchange of hexoses at the tonoplast of xylem-forming cells.

Limited plasticity of anatomical and hydraulic traits in aspen trees under elevated CO2 and seasonal drought.

The timing of abiotic stress elicitors on wood formation largely affects xylem traits that determine xylem efficiency and vulnerability. Nonetheless, seasonal variability of elevated CO2 (eCO2) effects on tree functioning under drought remains largely unknown. To address this knowledge gap, 1-year-old aspen (Populus tremula L.) trees were grown under ambient (±445 ppm) and elevated (±700 ppm) CO2 and exposed to an early (spring/summer 2019) or late (summer/autumn 2018) season drought event. Stomatal conductance and stem shrinkage were monitored in vivo as xylem water potential decreased. Additional trees were harvested for characterization of wood anatomical traits and to determine vulnerability and desorption curves via bench dehydration. The abundance of narrow vessels decreased under eCO2 only during the early season. At this time, xylem vulnerability to embolism formation and hydraulic capacitance during severe drought increased under eCO2. Contrastingly, stomatal closure was delayed during the late season, while hydraulic vulnerability and capacitance remained unaffected under eCO2. Independently of the CO2 treatment, elastic, and inelastic water pools depleted simultaneously after 50% of complete stomatal closure. Our results suggest that the effect of eCO2 on drought physiology and wood traits are small and variable during the growing season and question a sequential capacitive water release from elastic and inelastic pools as drought proceeds.

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.

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.

Exogenous Brassinosteroid Facilitates Xylem Development in Pinus massoniana Seedlings.

Brassinosteroids (BRs) are known to be essential regulators for wood formation in herbaceous plants and poplar, but their roles in secondary growth and xylem development are still not well-defined, especially in pines. Here, we treated Pinus massoniana seedlings with different concentrations of exogenous BRs, and assayed the effects on plant growth, xylem development, endogenous phytohormone contents and gene expression within stems. Application of exogenous BR resulted in improving development of xylem more than phloem, and promoting xylem development in a dosage-dependent manner in a certain concentration rage. Endogenous hormone determination showed that BR may interact with other phytohormones in regulating xylem development. RNA-seq analysis revealed that some conventional phenylpropanoid biosynthesis- or lignin synthesis-related genes were downregulated, but the lignin content was elevated, suggesting that new lignin synthesis pathways or other cell wall components should be activated by BR treatment in P. massoniana. The results presented here reveal the foundational role of BRs in regulating plant secondary growth, and provide the basis for understanding molecular mechanisms of xylem development in P. massoniana.

Abscisic acid signaling activates distinct VND transcription factors to promote xylem differentiation in Arabidopsis.

Plants display remarkable abilities to adjust growth and development to environmental conditions, such as the amount of available water. This developmental plasticity is apparent not only in root and shoot growth rates, but also in tissue patterning and cell morphology.1,2 We have previously shown that in response to limited water availability, Arabidopsis thaliana root displays changes in xylem morphology, mediated by the non-cell-autonomous action of abscisic acid, ABA.2 Here, we show, through analyses of ABA response reporters and tissue-specific suppression of ABA signaling, that xylem cells themselves act as primary signaling centers governing both xylem cell fate and xylem differentiation rate, revealing the cell-autonomous control of multiple aspects of xylem development by ABA. ABA rapidly activates the expression of genes encoding VASCULAR-RELATED NAC DOMAIN (VND) transcription factors. Molecular and genetic analyses revealed that the two ABA-mediated xylem developmental changes are regulated by distinct members of this transcription factor family, with VND2 and VND3 promoting differentiation rate of metaxylem cells, while VND7 promotes the conversion of metaxylem toward protoxylem morphology. This phenomenon shows how different aspects of developmental plasticity can be interlinked, yet genetically separable. Moreover, similarities in phenotypic and molecular responses to ABA in diverse species indicate evolutionary conservation of the ABA-xylem development regulatory network among eudicots. Hence, this study gives molecular insights into how environmental stress modifies plant vascular anatomy and has potential relevance for water use optimization and adaptation to drought conditions.

BdERECTA controls vasculature patterning and phloem-xylem organization in Brachypodium distachyon.

BACKGROUND: The vascular system of plants consists of two main tissue types, xylem and phloem. These tissues are organized into vascular bundles that are arranged into a complex network running through the plant that is essential for the viability of land plants. Despite their obvious importance, the genes involved in the organization of vascular tissues remain poorly understood in grasses. RESULTS: We studied in detail the vascular network in stems from the model grass Brachypodium distachyon (Brachypodium) and identified a large set of genes differentially expressed in vascular bundles versus parenchyma tissues. To decipher the underlying molecular mechanisms of vascularization in grasses, we conducted a forward genetic screen for abnormal vasculature. We identified a mutation that severely affected the organization of vascular tissues. This mutant displayed defects in anastomosis of the vascular network and uncommon amphivasal vascular bundles. The causal mutation is a premature stop codon in ERECTA, a LRR receptor-like serine/threonine-protein kinase. Mutations in this gene are pleiotropic indicating that it serves multiple roles during plant development. This mutant also displayed changes in cell wall composition, gene expression and hormone homeostasis. CONCLUSION: In summary, ERECTA has a pleiotropic role in Brachypodium. We propose a major role of ERECTA in vasculature anastomosis and vascular tissue organization in Brachypodium.

PagGRF12a interacts with PagGIF1b to regulate secondary xylem development through modulating PagXND1a expression in Populus alba × P. glandulosa.

Growth-regulating factors (GRFs) are important regulators of plant development and growth, but their possible roles in xylem development in woody plants remain unclear. Here, we report that Populus alba × Papulus glandulosa PagGRF12a negatively regulates xylem development in poplar. PagGRF12a is expressed in vascular tissues. Compared to non-transgenic control plants, transgenic poplar plants overexpressing PagGRF12a exhibited reduced xylem width and plants with repressed expression of PagGRF12a exhibited increased xylem width. Xylem NAC domain 1 (XND1) encodes a NAC domain transcription factor that regulates xylem development and transcriptional analyses revealed that PagXND1a is highly upregulated in PagGRF12a-overexpressing plants and downregulated in PagGRF12a-suppressed plants, indicating that PagGRF12a may regulate xylem development through PagXND1a. Transient transcriptional assays and chromatin immunoprecipitation-polymerase chain reaction assays confirmed that PagGRF12a directly upregulates PagXND1a. In addition, PagGRF12a interacts with the GRF-Interacting Factor (GIF) PagGIF1b, and this interaction enhances the effects of PagGRF12a on PagXND1a. Our results indicate that PagGRF12a inhibits xylem development by upregulating the expression of PagXND1a.

MYB Transcription Factors and Its Regulation in Secondary Cell Wall Formation and Lignin Biosynthesis during Xylem Development.

The secondary wall is the main part of wood and is composed of cellulose, xylan, lignin, and small amounts of structural proteins and enzymes. Lignin molecules can interact directly or indirectly with cellulose, xylan and other polysaccharide molecules in the cell wall, increasing the mechanical strength and hydrophobicity of plant cells and tissues and facilitating the long-distance transportation of water in plants. MYBs (v-myb avian myeloblastosis viral oncogene homolog) belong to one of the largest superfamilies of transcription factors, the members of which regulate secondary cell-wall formation by promoting/inhibiting the biosynthesis of lignin, cellulose, and xylan. Among them, MYB46 and MYB83, which comprise the second layer of the main switch of secondary cell-wall biosynthesis, coordinate upstream and downstream secondary wall synthesis-related transcription factors. In addition, MYB transcription factors other than MYB46/83, as well as noncoding RNAs, hormones, and other factors, interact with one another to regulate the biosynthesis of the secondary wall. Here, we discuss the biosynthesis of secondary wall, classification and functions of MYB transcription factors and their regulation of lignin polymerization and secondary cell-wall formation during wood formation.

Developmental and biophysical determinants of grass leaf size worldwide.

One of the most notable ecological trends-described more than 2,300  years ago by Theophrastus-is the association of small leaves with dry and cold climates, which has recently been recognized for eudicotyledonous plants at a global scale1-3. For eudicotyledons, this pattern has been attributed to the fact that small leaves have a thinner boundary layer that helps to avoid extreme leaf temperatures4 and their leaf development results in vein traits that improve water transport under cold or dry climates5,6. However, the global distribution of leaf size and its adaptive basis have not been tested in the grasses, which represent a diverse lineage that is distinct in leaf morphology and that contributes 33% of terrestrial primary productivity (including the bulk of crop production)7. Here we demonstrate that grasses have shorter and narrower leaves under colder and drier climates worldwide. We show that small grass leaves have thermal advantages and vein development that contrast with those of eudicotyledons, but that also explain the abundance of small leaves in cold and dry climates. The worldwide distribution of leaf size in grasses exemplifies how biophysical and developmental processes result in convergence across major lineages in adaptation to climate globally, and highlights the importance of leaf size and venation architecture for grass performance in past, present and future ecosystems.