Stem-cell-ubiquitous genes spatiotemporally coordinate division through regulation of stem-cell-specific gene networks

Deciphering the molecular logic of WOX5 function in the root stem cell organizer

Plant and animal stem cells receive signals from their surrounding cells to stay undifferentiated. In the Arabidopsis root, the quiescent center (QC) acts as a stem cell organizer, signaling to the neighboring stem cells. WOX5 is a central transcription factor regulating QC function. However, due to the scarcity of QC cells, WOX5 functions in the QC are largely unexplored at a genomic scale. Here, we unveil the transcriptional and epigenetic landscapes of the QC and the role of WOX5 within them. We find that WOX5 functions both as a transcriptional repressor and activator, affecting histone modifications and chromatin accessibility. Our data expand on known WOX5 functions, such as the regulation of differentiation, cell division, and auxin biosynthesis. We also uncover unexpected WOX5-regulated pathways involved in nitrate transport and the regulation of basal expression levels of genes associated with mature root tissues. These data suggest a role for QC cells as reserve stem cells and primed cells for prospective progenitor fates. Taken together, these findings offer insights into the role of WOX5 at the QC and provide a basis for further analyses to advance our understanding of the nature of plant stem cell organizers.

Machine Learning Inference of Gene Regulatory Networks in Developing Mimulus Seeds

The angiosperm seed represents a critical evolutionary breakthrough that has been shown to propel the reproductive success and radiation of flowering plants. Seeds promote the rapid diversification of angiosperms by establishing postzygotic reproductive barriers, such as hybrid seed inviability. While prezygotic barriers to reproduction tend to be transient, postzygotic barriers are often permanent and therefore can play a pivotal role in facilitating speciation. This property of the angiosperm seed is exemplified in the Mimulus genus. In order to further the understanding of the gene regulatory mechanisms important in the Mimulus seed, we performed gene regulatory network (GRN) inference analysis by using time-series RNA-seq data from developing hybrid seeds from a viable cross between Mimulus guttatus and Mimulus pardalis. GRN inference has the capacity to identify active regulatory mechanisms in a sample and highlight genes of potential biological importance. In our case, GRN inference also provided the opportunity to uncover active regulatory relationships and generate a reference set of putative gene regulations. We deployed two GRN inference algorithms-RTP-STAR and KBoost-on three different subsets of our transcriptomic dataset. While the two algorithms yielded GRNs with different regulations and topologies when working with the same data subset, there was still significant overlap in the specific gene regulations they inferred, and they both identified potential novel regulatory mechanisms that warrant further investigation.

Rewiring gene circuitry for plant improvement

Aspirations for high crop growth and yield, nutritional quality and bioproduction of materials are challenged by climate change and limited adoption of new technologies. Here, we review recent advances in approaches to profile and model gene regulatory activity over developmental and response time in specific cells, which have revealed the basis of variation in plant phenotypes: both redeployment of key regulators to new contexts and their repurposing to control different slates of genes. New synthetic biology tools allow tunable, spatiotemporal regulation of transgenes, while recent gene-editing technologies enable manipulation of the regulation of native genes. Ultimately, understanding how gene circuitry is wired to control form and function across varied plant species, combined with advanced technology to rewire that circuitry, will unlock solutions to our greatest challenges in agriculture, energy and the environment.

Changes in cell wall composition due to a pectin biosynthesis enzyme GAUT10 impact root growth

Arabidopsis (Arabidopsis thaliana) root development is regulated by multiple dynamic growth cues that require central metabolism pathways such as β-oxidation and auxin. Loss of the pectin biosynthesizing enzyme GALACTURONOSYLTRANSFERASE 10 (GAUT10) leads to a short-root phenotype under sucrose-limited conditions. The present study focused on determining the specific contributions of GAUT10 to pectin composition in primary roots and the underlying defects associated with gaut10 roots. Using live-cell microscopy, we determined reduced root growth in gaut10 is due to a reduction in both root apical meristem size and epidermal cell elongation. In addition, GAUT10 was required for normal pectin and hemicellulose composition in primary Arabidopsis roots. Specifically, loss of GAUT10 led to a reduction in galacturonic acid and xylose in root cell walls and altered the presence of rhamnogalacturonan-I (RG-I) and homogalacturonan (HG) polymers in the root. Transcriptomic analysis of gaut10 roots compared to wild type uncovered hundreds of genes differentially expressed in the mutant, including genes related to auxin metabolism and peroxisome function. Consistent with these results, both auxin signaling and metabolism were modified in gaut10 roots. The sucrose-dependent short-root phenotype in gaut10 was linked to β-oxidation based on hypersensitivity to indole-3-butyric acid (IBA) and an epistatic interaction with TRANSPORTER OF IBA1 (TOB1). Altogether, these data support a growing body of evidence suggesting that pectin composition may influence auxin pathways and peroxisome activity.

Network Biology Analyses and Dynamic Modeling of Gene Regulatory Networks under Drought Stress Reveal Major Transcriptional Regulators in Arabidopsis

Drought is one of the most serious abiotic stressors in the environment, restricting agricultural production by reducing plant growth, development, and productivity. To investigate such a complex and multifaceted stressor and its effects on plants, a systems biology-based approach is necessitated, entailing the generation of co-expression networks, identification of high-priority transcription factors (TFs), dynamic mathematical modeling, and computational simulations. Here, we studied a high-resolution drought transcriptome of Arabidopsis. We identified distinct temporal transcriptional signatures and demonstrated the involvement of specific biological pathways. Generation of a large-scale co-expression network followed by network centrality analyses identified 117 TFs that possess critical properties of hubs, bottlenecks, and high clustering coefficient nodes. Dynamic transcriptional regulatory modeling of integrated TF targets and transcriptome datasets uncovered major transcriptional events during the course of drought stress. Mathematical transcriptional simulations allowed us to ascertain the activation status of major TFs, as well as the transcriptional intensity and amplitude of their target genes. Finally, we validated our predictions by providing experimental evidence of gene expression under drought stress for a set of four TFs and their major target genes using qRT-PCR. Taken together, we provided a systems-level perspective on the dynamic transcriptional regulation during drought stress in Arabidopsis and uncovered numerous novel TFs that could potentially be used in future genetic crop engineering programs.

Cell Cycle Dynamics during Stomatal Development: Window of MUTE Action and Ramification of Its Loss-of-Function on an Uncommitted Precursor

Plants develop in the absence of cell migration. As such, cell division and differentiation need to be coordinated for functional tissue formation. Cellular valves on the plant epidermis, stomata, are generated through a stereotypical sequence of cell division and differentiation events. In Arabidopsis, three master regulatory transcription factors, SPEECHLESS (SPCH), MUTE and FAMA, sequentially drive initiation, proliferation and differentiation of stomata. Among them, MUTE switches the cell cycle mode from proliferative asymmetric division to terminal symmetric division and orchestrates the execution of the single symmetric division event. However, it remains unclear to what extent MUTE regulates the expression of cell cycle genes through the symmetric division and whether MUTE accumulation itself is gated by the cell cycle. Here, we show that MUTE directly upregulates the expression of cell cycle components throughout the terminal cell cycle phases of a stomatal precursor, not only core cell cycle engines but also check-point regulators. Time-lapse live imaging using the multicolor Plant Cell Cycle Indicator revealed that MUTE accumulates up to the early G2 phase, whereas its successor and direct target, FAMA, accumulate at late G2 through terminal mitosis. In the absence of MUTE, meristemoids fail to differentiate and their G1 phase elongates as they reiterate asymmetric divisions. Together, our work provides the framework of cell cycle and master regulatory transcription factors to coordinate a single symmetric cell division and suggests a mechanism for the eventual cell cycle arrest of an uncommitted stem-cell-like precursor at the G1 phase.

Charting plant gene functions in the multi-omics and single-cell era

Despite the increased access to high-quality plant genome sequences, the set of genes with a known function remains far from complete. With the advent of novel bulk and single-cell omics profiling methods, we are entering a new era where advanced and highly integrative functional annotation strategies are being developed to elucidate the functions of all plant genes. Here, we review different multi-omics approaches to improve functional and regulatory gene characterization and highlight the power of machine learning and network biology to fully exploit the complementary information embedded in different omics layers. Finally, we discuss the potential of emerging single-cell methods and algorithms to further increase the resolution, allowing generation of functional insights about plant biology.

Biomaterials based on hyaluronic acid, collagen and peptides for three-dimensional cell culture and their application in stem cell differentiation

In recent decades, three-dimensional (3D) cell culture technologies have been developed rapidly in the field of tissue engineering and regeneration, and have shown unique advantages and great prospects in the differentiation of stem cells. Herein, the article reviews the progress and advantages of 3D cell culture technologies in the field of stem cell differentiation. Firstly, 3D cell culture technologies are divided into two main categories: scaffoldless and scaffolds. Secondly, the effects of hydrogels scaffolds and porous scaffolds on stem cell differentiation in the scaffold category were mainly reviewed. Among them, hydrogels scaffolds are divided into natural hydrogels and synthetic hydrogels. Natural materials include polysaccharides, proteins, and their derivatives, focusing on hyaluronic acid, collagen and polypeptides. Synthetic materials mainly include polyethylene glycol (PEG), polyacrylic acid (PAA), polyvinyl alcohol (PVA), etc. In addition, since the preparation techniques have a large impact on the properties of porous scaffolds, several techniques for preparing porous scaffolds based on different macromolecular materials are reviewed. Finally, the future prospects and challenges of 3D cell culture in the field of stem cell differentiation are reviewed. This review will provide a useful guideline for the selection of materials and techniques for 3D cell culture in stem cell differentiation.

A Practical Guide to Inferring Multi-Omics Networks in Plant Systems

The inference of gene regulatory networks can reveal molecular connections underlying biological processes and improve our understanding of complex biological phenomena in plants. Many previous network studies have inferred networks using only one type of omics data, such as transcriptomics. However, given more recent work applying multi-omics integration in plant biology, such as combining (phospho)proteomics with transcriptomics, it may be advantageous to integrate multiple omics data types into a comprehensive network prediction. Here, we describe a state-of-the-art approach for integrating multi-omics data with gene regulatory network inference to describe signaling pathways and uncover novel regulators. We detail how to download and process transcriptomics and (phospho)proteomics data for network inference, using an example dataset from the plant hormone signaling field. We provide a step-by-step protocol for inference, visualization, and analysis of an integrative multi-omics network using currently available methods. This chapter serves as an accessible guide for novice and intermediate bioinformaticians to analyze their own datasets and reanalyze published work.

POPEYE intercellular localization mediates cell-specific iron deficiency responses

Plants must tightly regulate iron (Fe) sensing, acquisition, transport, mobilization, and storage to ensure sufficient levels of this essential micronutrient. POPEYE (PYE) is an iron responsive transcription factor that positively regulates the iron deficiency response, while also repressing genes essential for maintaining iron homeostasis. However, little is known about how PYE plays such contradictory roles. Under iron-deficient conditions, pPYE:GFP accumulates in the root pericycle while pPYE:PYE-GFP is localized to the nucleus in all Arabidopsis (Arabidopsis thaliana) root cells, suggesting that PYE may have cell-specific dynamics and functions. Using scanning fluorescence correlation spectroscopy and cell-specific promoters, we found that PYE-GFP moves between different cells and that the tendency for movement corresponds with transcript abundance. While localization to the cortex, endodermis, and vasculature is required to manage changes in iron availability, vasculature and endodermis localization of PYE-GFP protein exacerbated pye-1 defects and elicited a host of transcriptional changes that are detrimental to iron mobilization. Our findings indicate that PYE acts as a positive regulator of iron deficiency response by regulating iron bioavailability differentially across cells, which may trigger iron uptake from the surrounding rhizosphere and impact root energy metabolism.

Temporal and spatial auxin responsive networks in maize primary roots

Auxin is a key regulator of root morphogenesis across angiosperms. To better understand auxin-regulated networks underlying maize root development, we have characterized auxin-responsive transcription across two time points (30 and 120 min) and four regions of the primary root: the meristematic zone, elongation zone, cortex and stele. Hundreds of auxin-regulated genes involved in diverse biological processes were quantified in these different root regions. In general, most auxin-regulated genes are region unique and are predominantly observed in differentiated tissues compared with the root meristem. Auxin gene regulatory networks were reconstructed with these data to identify key transcription factors that may underlie auxin responses in maize roots. Additionally, Auxin-Response Factor subnetworks were generated to identify target genes that exhibit tissue or temporal specificity in response to auxin. These networks describe novel molecular connections underlying maize root development and provide a foundation for functional genomic studies in a key crop.

Integrative omics approaches for biosynthetic pathway discovery in plants

Covering: up to 2022With the emergence of large amounts of omics data, computational approaches for the identification of plant natural product biosynthetic pathways and their genetic regulation have become increasingly important. While genomes provide clues regarding functional associations between genes based on gene clustering, metabolome mining provides a foundational technology to chart natural product structural diversity in plants, and transcriptomics has been successfully used to identify new members of their biosynthetic pathways based on coexpression. Thus far, most approaches utilizing transcriptomics and metabolomics have been targeted towards specific pathways and use one type of omics data at a time. Recent technological advances now provide new opportunities for integration of multiple omics types and untargeted pathway discovery. Here, we review advances in plant biosynthetic pathway discovery using genomics, transcriptomics, and metabolomics, as well as recent efforts towards omics integration. We highlight how transcriptomics and metabolomics provide complementary information to link genes to metabolites, by associating temporal and spatial gene expression levels with metabolite abundance levels across samples, and by matching mass-spectral features to enzyme families. Furthermore, we suggest that elucidation of gene regulatory networks using time-series data may prove useful for efforts to unwire the complexities of biosynthetic pathway components based on regulatory interactions and events.

BRAVO self-confined expression through WOX5 in the Arabidopsis root stem-cell niche

In animals and plants, stem-cell niches are local microenvironments that are tightly regulated to preserve their unique identity while communicating with adjacent cells that will give rise to specialized cell types. In the primary root of Arabidopsis thaliana, two transcription factors, BRAVO and WOX5, among others, are expressed in the stem-cell niche. Intriguingly, BRAVO, a repressor of quiescent center divisions, confines its own gene expression to the stem-cell niche, as evidenced in a bravo mutant background. Here, we propose through mathematical modeling that BRAVO confines its own expression domain to the stem-cell niche by attenuating a WOX5-dependent diffusible activator of BRAVO. This negative feedback drives WOX5 activity to be spatially restricted as well. The results show that WOX5 diffusion and sequestration by binding to BRAVO are sufficient to drive the experimentally observed confined BRAVO expression at the stem-cell niche. We propose that the attenuation of a diffusible activator can be a general mechanism acting at other stem-cell niches to spatially confine genetic activity to a small region while maintaining signaling within them and with the surrounding cells.

Gene regulatory networks shape developmental plasticity of root cell types under water extremes in rice

Understanding how roots modulate development under varied irrigation or rainfall is crucial for development of climate-resilient crops. We established a toolbox of tagged rice lines to profile translating mRNAs and chromatin accessibility within specific cell populations. We used these to study roots in a range of environments: plates in the lab, controlled greenhouse stress and recovery conditions, and outdoors in a paddy. Integration of chromatin and mRNA data resolves regulatory networks of the following: cycle genes in proliferating cells that attenuate DNA synthesis under submergence; genes involved in auxin signaling, the circadian clock, and small RNA regulation in ground tissue; and suberin biosynthesis, iron transporters, and nitrogen assimilation in endodermal/exodermal cells modulated with water availability. By applying a systems approach, we identify known and candidate driver transcription factors of water-deficit responses and xylem development plasticity. Collectively, this resource will facilitate genetic improvements in root systems for optimal climate resilience.

RNA-Seq and Gene Regulatory Network Analyses Uncover Candidate Genes in the Early Defense to Two Hemibiotrophic Colletorichum spp. in Strawberry

Two hemibiotrophic pathogens, Colletotrichum acutatum (Ca) and C. gloeosporioides (Cg), cause anthracnose fruit rot and anthracnose crown rot in strawberry (Fragaria × ananassa Duchesne), respectively. Both Ca and Cg can initially infect through a brief biotrophic phase, which is associated with the production of intracellular primary hyphae that can infect host cells without causing cell death and establishing hemibiotrophic infection (HBI) or quiescent (latent infections) in leaf tissues. The Ca and Cg HBI in nurseries and subsequent distribution of asymptomatic infected transplants to fruit production fields is the major source of anthracnose epidemics in North Carolina. In the absence of complete resistance, strawberry varieties with good fruit quality showing rate-reducing resistance have frequently been used as a source of resistance to Ca and Cg. However, the molecular mechanisms underlying the rate-reducing resistance or susceptibility to Ca and Cg are still unknown. We performed comparative transcriptome analyses to examine how rate-reducing resistant genotype NCS 10-147 and susceptible genotype ‘Chandler’ respond to Ca and Cg and identify molecular events between 0 and 48 h after the pathogen-inoculated and mock-inoculated leaf tissues. Although plant response to both Ca and Cg at the same timepoint was not similar, more genes in the resistant interaction were upregulated at 24 hpi with Ca compared with those at 48 hpi. In contrast, a few genes were upregulated in the resistant interaction at 48 hpi with Cg. Resistance response to both Ca and Cg was associated with upregulation of MLP-like protein 44, LRR receptor-like serine/threonine-protein kinase, and auxin signaling pathway, whereas susceptibility was linked to modulation of the phenylpropanoid pathway. Gene regulatory network inference analysis revealed candidate transcription factors (TFs) such as GATA5 and MYB-10, and their downstream targets were upregulated in resistant interactions. Our results provide valuable insights into transcriptional changes during resistant and susceptible interactions, which can further facilitate assessing candidate genes necessary for resistance to two hemibiotrophic Colletotrichum spp. in strawberry.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

A single-cell Arabidopsis root atlas reveals developmental trajectories in wild-type and cell identity mutants

In all multicellular organisms, transcriptional networks orchestrate organ development. The Arabidopsis root, with its simple structure and indeterminate growth, is an ideal model for investigating the spatiotemporal transcriptional signatures underlying developmental trajectories. To map gene expression dynamics across root cell types and developmental time, we built a comprehensive, organ-scale atlas at single-cell resolution. In addition to estimating developmental progressions in pseudotime, we employed the mathematical concept of optimal transport to infer developmental trajectories and identify their underlying regulators. To demonstrate the utility of the atlas to interpret new datasets, we profiled mutants for two key transcriptional regulators at single-cell resolution, shortroot and scarecrow. We report transcriptomic and in vivo evidence for tissue trans-differentiation underlying a mixed cell identity phenotype in scarecrow. Our results support the atlas as a rich community resource for unraveling the transcriptional programs that specify and maintain cell identity to regulate spatiotemporal organ development.

Analyzing immune response to engineered hydrogels by hierarchical clustering of inflammatory cell subsets

Understanding the immune response to hydrogel implantation is critical for the design of immunomodulatory biomaterials. To study the progression of inflammation around poly(ethylene glycol) hydrogels presenting Arg-Gly-Asp (RGD) peptides and vascular endothelial growth factor, we used temporal analysis of high-dimensional flow cytometry data paired with intravital imaging, immunohistochemistry, and multiplexed proteomic profiling. RGD-presenting hydrogels created a reparative microenvironment promoting CD206+ cellular infiltration and revascularization in wounded dorsal skin tissue. Unbiased clustering algorithms (SPADE) revealed significant phenotypic transition shifts as a function of the cell-adhesion hydrogel properties. SPADE identified an intermediate macrophage subset functionally regulating in vivo cytokine secretion that was preferentially recruited for RGD-presenting hydrogels, whereas dendritic cell subsets were preferentially recruited to RDG-presenting hydrogels. Last, RGD-presenting hydrogels controlled macrophage functional cytokine secretion to direct polarization and vascularization. Our studies show that unbiased clustering of single-cell data provides unbiased insights into the underlying immune response to engineered materials.

Comparative Analysis of Seed Transcriptome and Coexpression Analysis Reveal Candidate Genes for Enhancing Seed Size/Weight in Brassica juncea

Seed size/weight is a multigenic trait that is governed by complex transcriptional regulatory pathways. An understanding of the genetic basis of seed size is of great interest in the improvement of seed yield and quality in oilseed crops. A global transcriptome analysis was performed at the initial stages of seed development in two lines of Brassica juncea, small-seeded EH-2 and large-seeded PJ. The anatomical analyses revealed significant differences in cell number and cell size in the outer layer of the seed coat between EH-2 and PJ. Pairwise comparisons at each developmental stage identified 5,974 differentially expressed genes (DEGs) between the two lines, of which 954 genes belong to different families of transcription factors. Two modules were found to be significantly correlated with an increased seed size using weighted gene coexpression network analysis. The DEG and coexpression datasets were integrated with the thousand seed weight (Tsw) quantitative trait loci (QTL) mapped earlier in the EPJ (EH-2 × PJ) doubled haploid (DH) population, which identified forty potential key components controlling seed size. The candidate genes included genes regulating the cell cycle, cell wall biogenesis/modification, solute/sugar transport, and hormone signaling. The results provide a valuable resource to widen the current understanding of regulatory mechanisms underlying seed size in B. juncea.

Gene regulatory networks for compatible versus incompatible grafts identify a role for SlWOX4 during junction formation

Grafting has been adopted for a wide range of crops to enhance productivity and resilience; for example, grafting of Solanaceous crops couples disease-resistant rootstocks with scions that produce high-quality fruit. However, incompatibility severely limits the application of grafting and graft incompatibility remains poorly understood. In grafts, immediate incompatibility results in rapid death, but delayed incompatibility can take months or even years to manifest, creating a significant economic burden for perennial crop production. To gain insight into the genetic mechanisms underlying this phenomenon, we developed a model system using heterografting of tomato (Solanum lycopersicum) and pepper (Capsicum annuum). These grafted plants express signs of anatomical junction failure within the first week of grafting. By generating a detailed timeline for junction formation, we were able to pinpoint the cellular basis for this delayed incompatibility. Furthermore, we inferred gene regulatory networks for compatible self-grafts and incompatible heterografts based on these key anatomical events, which predict core regulators for grafting. Finally, we examined the role of vascular development in graft formation and uncovered SlWOX4 as a potential regulator of graft compatibility. Following this predicted regulator up with functional analysis, we show that Slwox4 homografts fail to form xylem bridges across the junction, demonstrating that indeed, SlWOX4 is essential for vascular reconnection during grafting, and may function as an early indicator of graft failure.

Gene regulatory networks for compatible versus incompatible grafts identify a role for SlWOX4 during junction formation

Grafting has been adopted for a wide range of crops to enhance productivity and resilience; for example, grafting of Solanaceous crops couples disease-resistant rootstocks with scions that produce high-quality fruit. However, incompatibility severely limits the application of grafting and graft incompatibility remains poorly understood. In grafts, immediate incompatibility results in rapid death, but delayed incompatibility can take months or even years to manifest, creating a significant economic burden for perennial crop production. To gain insight into the genetic mechanisms underlying this phenomenon, we developed a model system using heterografting of tomato (Solanum lycopersicum) and pepper (Capsicum annuum). These grafted plants express signs of anatomical junction failure within the first week of grafting. By generating a detailed timeline for junction formation, we were able to pinpoint the cellular basis for this delayed incompatibility. Furthermore, we inferred gene regulatory networks for compatible self-grafts and incompatible heterografts based on these key anatomical events, which predict core regulators for grafting. Finally, we examined the role of vascular development in graft formation and uncovered SlWOX4 as a potential regulator of graft compatibility. Following this predicted regulator up with functional analysis, we show that Slwox4 homografts fail to form xylem bridges across the junction, demonstrating that indeed, SlWOX4 is essential for vascular reconnection during grafting, and may function as an early indicator of graft failure.

Future Challenges in Plant Systems Biology

Plant systems biology is currently facing several important challenges, whose nature depend on the considered frame of reference and associated scale. This review covers some of the issues associated respectively with the molecular, tissue, and whole-plant scales, as well as discusses the potential for latest advances in synthetic biology and machine-learning methods to be of use in the future of plant systems biology.

Regulation of stem cell fate and function by using bioactive materials with nanoarchitectonics for regenerative medicine

Nanoarchitectonics has emerged as a post-nanotechnology concept. As one of the applications of nanoarchitectonics, this review paper discusses the control of stem cell fate and function as an important issue. For hybrid nanoarchitectonics involving living cells, it is crucial to understand how biomaterials and their nanoarchitected structures regulate behaviours and fates of stem cells. In this review, biomaterials for the regulation of stem cell fate are firstly discussed. Besides multipotent differentiation, immunomodulation is an important biological function of mesenchymal stem cells (MSCs). MSCs can modulate immune cells to treat multiple immune- and inflammation-mediated diseases. The following sections summarize the recent advances of the regulation of the immunomodulatory functions of MSCs by biophysical signals. In the third part, we discussed how biomaterials direct the self-organization of pluripotent stem cells for organoid. Bioactive materials are constructed which mimic the biophysical cues of in vivo microenvironment such as elasticity, viscoelasticity, biodegradation, fluidity, topography, cell geometry, and etc. Stem cells interpret these biophysical cues by different cytoskeletal forces. The different cytoskeletal forces lead to substantial transcription and protein expression, which affect stem cell fate and function. Regulations of stem cells could not be utilized only for tissue repair and regenerative medicine but also potentially for production of advanced materials systems. Materials nanoarchitectonics with integration of stem cells and related biological substances would have high impacts in science and technology of advanced materials.

Integrated omics networks reveal the temporal signaling events of brassinosteroid response in Arabidopsis

Brassinosteroids (BRs) are plant steroid hormones that regulate cell division and stress response. Here we use a systems biology approach to integrate multi-omic datasets and unravel the molecular signaling events of BR response in Arabidopsis. We profile the levels of 26,669 transcripts, 9,533 protein groups, and 26,617 phosphorylation sites from Arabidopsis seedlings treated with brassinolide (BL) for six different lengths of time. We then construct a network inference pipeline called Spatiotemporal Clustering and Inference of Omics Networks (SC-ION) to integrate these data. We use our network predictions to identify putative phosphorylation sites on BES1 and experimentally validate their importance. Additionally, we identify BRONTOSAURUS (BRON) as a transcription factor that regulates cell division, and we show that BRON expression is modulated by BR-responsive kinases and transcription factors. This work demonstrates the power of integrative network analysis applied to multi-omic data and provides fundamental insights into the molecular signaling events occurring during BR response.

The DREAM complex represses growth in response to DNA damage in Arabidopsis

The DNA of all organisms is constantly damaged by physiological processes and environmental conditions. Upon persistent damage, plant growth and cell proliferation are reduced. Based on previous findings that RBR1, the only Arabidopsis homolog of the mammalian tumor suppressor gene retinoblastoma, plays a key role in the DNA damage response in plants, we unravel here the network of RBR1 interactors under DNA stress conditions. This led to the identification of homologs of every DREAM component in Arabidopsis, including previously not recognized homologs of LIN52. Interestingly, we also discovered NAC044, a mediator of DNA damage response in plants and close homolog of the major DNA damage regulator SOG1, to directly interact with RBR1 and the DREAM component LIN37B. Consistently, not only mutants in NAC044 but also the double mutant of the two LIN37 homologs and mutants for the DREAM component E2FB showed reduced sensitivities to DNA-damaging conditions. Our work indicates the existence of multiple DREAM complexes that work in conjunction with NAC044 to mediate growth arrest after DNA damage.

Spatiotemporal Gene Expression Profiling and Network Inference: A Roadmap for Analysis, Visualization, and Key Gene Identification

Gene expression data analysis and the prediction of causal relationships within gene regulatory networks (GRNs) have guided the identification of key regulatory factors and unraveled the dynamic properties of biological systems. However, drawing accurate and unbiased conclusions requires a comprehensive understanding of relevant tools, computational methods, and their workflows. The topics covered in this chapter encompass the entire workflow for GRN inference including: (1) experimental design; (2) RNA sequencing data processing; (3) differentially expressed gene (DEG) selection; (4) clustering prior to inference; (5) network inference techniques; and (6) network visualization and analysis. Moreover, this chapter aims to present a workflow feasible and accessible for plant biologists without a bioinformatics or computer science background. To address this need, TuxNet, a user-friendly graphical user interface that integrates RNA sequencing data analysis with GRN inference, is chosen for the purpose of providing a detailed tutorial.

Surface acoustic wave (SAW) techniques in tissue engineering

Recently, the introduction of surface acoustic wave (SAW) technique for microfluidics has drawn a lot of attention. The pattern and mutual communication in cell layers, tissues, and organs play a critical role in tissue homeostasis and regeneration and may contribute to disease occurrence and progression. Tissue engineering aims to repair and regenerate damaged organs, depending on biomimetic scaffolds and advanced fabrication technology. However, traditional bioengineering synthesis approaches are time-consuming, heterogeneous, and unmanageable. It is hard to pattern cells in scaffolds effectively with no impact on cell viability and function. Here, we summarize a biocompatible, easily available, label-free, and non-invasive tool, surface acoustic wave (SAW) technique, which is getting a lot of attention in tissue engineering. SAW technique can realize accurate sorting, manipulation, and cells' pattern and rapid formation of spheroids. By integrating several SAW devices onto lab-on-a-chip platforms, tissue engineering lab-on-a-chip system was proposed. To the best of our knowledge, this is the first report to summarize the application of this novel technique in the field of tissue engineering.

Reconstruction of lateral root formation through single-cell RNA sequencing reveals order of tissue initiation

Postembryonic organogenesis is critical for plant development. Underground, lateral roots (LRs) form the bulk of mature root systems, yet the ontogeny of the LR primordium (LRP) is not clear. In this study, we performed the single-cell RNA sequencing through the first four stages of LR formation in Arabidopsis. Our analysis led to a model in which a single group of precursor cells, with a cell identity different from their pericycle origins, rapidly reprograms and splits into a mixed ground tissue/stem cell niche fate and a vascular precursor fate. The ground tissue and stem cell niche fates soon separate and a subset of more specialized vascular cells form sucrose transporting phloem cells that appear to connect to the primary root. We did not detect cells resembling epidermis or root cap, suggesting that outer tissues may form later, preceding LR emergence. At this stage, some remaining initial precursor cells form the primordium flanks, while the rest create a reservoir of pluripotent cells that are able to replace the LR if damaged. Laser ablation of the central and lateral LRP regions showed that remaining cells restart the sequence of tissue initiation to form a LR. Collectively, our study reveals an ontological hierarchy for LR formation with an early and sequential split of main root tissues and stem cells.

Metaphloem development in the Arabidopsis root tip

The phloem transport network is a major evolutionary innovation that enabled plants to dominate terrestrial ecosystems. In the growth apices, the meristems, apical stem cells continuously produce early 'protophloem'. This is easily observed in Arabidopsis root meristems, in which the differentiation of individual protophloem sieve element precursors into interconnected conducting sieve tubes is laid out in a spatio-temporal gradient. The mature protophloem eventually collapses as the neighboring metaphloem takes over its function further distal from the stem cell niche. Compared with protophloem, metaphloem ontogenesis is poorly characterized, primarily because its visualization is challenging. Here, we describe the improved TetSee protocol to investigate metaphloem development in Arabidopsis root tips in combination with a set of molecular markers. We found that mature metaphloem sieve elements are only observed in the late post-meristematic root, although their specification is initiated as soon as protophloem sieve elements enucleate. Moreover, unlike protophloem sieve elements, metaphloem sieve elements only differentiate once they have fully elongated. Finally, our results suggest that metaphloem differentiation is not directly controlled by protophloem-derived cues but rather follows a distinct, robust developmental trajectory.

Transcriptome profiling reveals the spatial-temporal dynamics of gene expression essential for soybean seed development

BACKGROUND

Seeds are the economic basis of oilseed crops, especially soybeans, the most widely cultivated oilseed crop worldwide. Seed development is accompanied by a multitude of diverse cellular processes, and revealing the underlying regulatory activities is critical for seed improvement.

RESULTS

In this study, we profiled the transcriptomes of developing seeds at 20, 25, 30, and 40 days after flowering (DAF), as these stages represent critical time points of seed development from early to full development. We identified a set of highly abundant genes and highlighted the importance of these genes in supporting nutrient accumulation and transcriptional regulation for seed development. We identified 8925 differentially expressed genes (DEGs) that exhibited temporal expression patterns over the course and expression specificities in distinct tissues, including seeds and nonseed tissues (roots, stems, and leaves). Genes specific to nonseed tissues might have tissue-associated roles, with relatively low transcript abundance in developing seeds, suggesting their spatially supportive roles in seed development. Coexpression network analysis identified several underexplored genes in soybeans that bridge tissue-specific gene modules.

CONCLUSIONS

Our study provides a global view of gene activities and biological processes critical for seed formation in soybeans and prioritizes a set of genes for further study. The results of this study help to elucidate the mechanism controlling seed development and storage reserves.

Precise transcriptional control of cellular quiescence by BRAVO/WOX5 complex in Arabidopsis roots

Understanding stem cell regulatory circuits is the next challenge in plant biology, as these cells are essential for tissue growth and organ regeneration in response to stress. In the Arabidopsis primary root apex, stem cell-specific transcription factors BRAVO and WOX5 co-localize in the quiescent centre (QC) cells, where they commonly repress cell division so that these cells can act as a reservoir to replenish surrounding stem cells, yet their molecular connection remains unknown. Genetic and biochemical analysis indicates that BRAVO and WOX5 form a transcription factor complex that modulates gene expression in the QC cells to preserve overall root growth and architecture. Furthermore, by using mathematical modelling we establish that BRAVO uses the WOX5/BRAVO complex to promote WOX5 activity in the stem cells. Our results unveil the importance of transcriptional regulatory circuits in plant stem cell development.

Unraveling Root Development Through Single-Cell Omics and Reconstruction of Gene Regulatory Networks

Over the last decades, research on postembryonic root development has been facilitated by "omics" technologies. Among these technologies, microarrays first, and RNA sequencing (RNA-seq) later, have provided transcriptional information on the underlying molecular processes establishing the basis of System Biology studies in roots. Cell fate specification and development have been widely studied in the primary root, which involved the identification of many cell type transcriptomes and the reconstruction of gene regulatory networks (GRN). The study of lateral root (LR) development has not been an exception. However, the molecular mechanisms regulating cell fate specification during LR formation remain largely unexplored. Recently, single-cell RNA-seq (scRNA-seq) studies have addressed the specification of tissues from stem cells in the primary root. scRNA-seq studies are anticipated to be a useful approach to decipher cell fate specification and patterning during LR formation. In this review, we address the different scRNA-seq strategies used both in plants and animals and how we could take advantage of scRNA-seq to unravel new regulatory mechanisms and reconstruct GRN. In addition, we discuss how to integrate scRNA-seq results with previous RNA-seq datasets and GRN. We also address relevant findings obtained through single-cell based studies and how LR developmental studies could be facilitated by scRNA-seq approaches and subsequent GRN inference. The use of single-cell approaches to investigate LR formation could help to decipher fundamental biological mechanisms such as cell memory, synchronization, polarization, or pluripotency.

A hybrid model connecting regulatory interactions with stem cell divisions in the root

Stem cells give rise to the entirety of cells within an organ. Maintaining stem cell identity and coordinately regulating stem cell divisions is crucial for proper development. In plants, mobile proteins, such as WUSCHEL-RELATED HOMEOBOX 5 (WOX5) and SHORTROOT (SHR), regulate divisions in the root stem cell niche. However, how these proteins coordinately function to establish systemic behaviour is not well understood. We propose a non-cell autonomous role for WOX5 in the cortex endodermis initial (CEI) and identify a regulator, ANGUSTIFOLIA (AN3)/GRF-INTERACTING FACTOR 1, that coordinates CEI divisions. Here, we show with a multi-scale hybrid model integrating ordinary differential equations (ODEs) and agent-based modeling that quiescent center (QC) and CEI divisions have different dynamics. Specifically, by combining continuous models to describe regulatory networks and agent-based rules, we model systemic behaviour, which led us to predict cell-type-specific expression dynamics of SHR, SCARECROW, WOX5, AN3 and CYCLIND6;1, and experimentally validate CEI cell divisions. Conclusively, our results show an interdependency between CEI and QC divisions.

BAM1/2 receptor kinase signaling drives CLE peptide-mediated formative cell divisions in Arabidopsis roots

Proper elaboration of the plant body plan requires that cell division patterns are coordinated during development in complex tissues. Activation of cell cycle machinery is critical for this process, but it is not clear how or if this links to cell-to-cell communication networks that are important during development. Here we show that key cell divisions that generate the plant root are controlled by cell-to-cell signaling peptides which act through plant-specific receptor kinases to control expression of a specific cyclinD cell cycle regulatory gene. We show that cyclinD gene expression depends on both receptor signaling and the SHORT-ROOT transcription factor to ensure timely and robust cell division patterns.

Cell division is often regulated by extracellular signaling networks to ensure correct patterning during development. In Arabidopsis, the SHORT-ROOT (SHR)/SCARECROW (SCR) transcription factor dimer activates CYCLIND6;1 (CYCD6;1) to drive formative divisions during root ground tissue development. Here, we show plasma-membrane-localized BARELY ANY MERISTEM1/2 (BAM1/2) family receptor kinases are required for SHR-dependent formative divisions and CYCD6;1 expression, but not SHR-dependent ground tissue specification. Root-enriched CLE ligands bind the BAM1 extracellular domain and are necessary and sufficient to activate SHR-mediated divisions and CYCD6;1 expression. Correspondingly, BAM-CLE signaling contributes to the restriction of formative divisions to the distal root region. Additionally, genetic analysis reveals that BAM-CLE and SHR converge to regulate additional cell divisions outside of the ground tissues. Our work identifies an extracellular signaling pathway regulating formative root divisions and provides a framework to explore this pathway in patterning and evolution.

The Wheat GENIE3 Network Provides Biologically-Relevant Information in Polyploid Wheat

Gene regulatory networks are powerful tools which facilitate hypothesis generation and candidate gene discovery. However, the extent to which the network predictions are biologically relevant is often unclear. Recently a GENIE3 network which predicted targets of wheat transcription factors was produced. Here we used an independent RNA-Seq dataset to test the predictions of the wheat GENIE3 network for the senescence-regulating transcription factor NAM-A1 (TraesCS6A02G108300). We re-analyzed the RNA-Seq data against the RefSeqv1.0 genome and identified a set of differentially expressed genes (DEGs) between the wild-type and nam-a1 mutant which recapitulated the known role of NAM-A1 in senescence and nutrient remobilisation. We found that the GENIE3-predicted target genes of NAM-A1 overlap significantly with the DEGs, more than would be expected by chance. Based on high levels of overlap between GENIE3-predicted target genes and the DEGs, we identified candidate senescence regulators. We then explored genome-wide trends in the network related to polyploidy and found that only homeologous transcription factors are likely to share predicted targets in common. However, homeologs which vary in expression levels across tissues are less likely to share predicted targets than those that do not, suggesting that they may be more likely to act in distinct pathways. This work demonstrates that the wheat GENIE3 network can provide biologically-relevant predictions of transcription factor targets, which can be used for candidate gene prediction and for global analyses of transcription factor function. The GENIE3 network has now been integrated into the KnetMiner web application, facilitating its use in future studies.

Protein complex stoichiometry and expression dynamics of transcription factors modulate stem cell division

Stem cells divide and differentiate to form all of the specialized cell types in a multicellular organism. In the Arabidopsis root, stem cells are maintained in an undifferentiated state by a less mitotically active population of cells called the quiescent center (QC). Determining how the QC regulates the surrounding stem cell initials, or what makes the QC fundamentally different from the actively dividing initials, is important for understanding how stem cell divisions are maintained. Here we gained insight into the differences between the QC and the cortex endodermis initials (CEI) by studying the mobile transcription factor SHORTROOT (SHR) and its binding partner SCARECROW (SCR). We constructed an ordinary differential equation model of SHR and SCR in the QC and CEI which incorporated the stoichiometry of the SHR-SCR complex as well as upstream transcriptional regulation of SHR and SCR. Our model prediction, coupled with experimental validation, showed that high levels of the SHR-SCR complex are associated with more CEI division but less QC division. Furthermore, our model prediction allowed us to propose the putative upstream SHR regulators SEUSS and WUSCHEL-RELATED HOMEOBOX 5 and to experimentally validate their roles in QC and CEI division. In addition, our model established the timing of QC and CEI division and suggests that SHR repression of QC division depends on formation of the SHR homodimer. Thus, our results support that SHR-SCR protein complex stoichiometry and regulation of SHR transcription modulate the division timing of two different specialized cell types in the root stem cell niche.

Integrated multi-omics framework of the plant response to jasmonic acid

Understanding the systems-level actions of transcriptional responses to hormones provides insight into how the genome is reprogrammed in response to environmental stimuli. Here, we investigated the signalling pathway of the hormone jasmonic acid (JA), which controls a plethora of critically important processes in plants and is orchestrated by the transcription factor MYC2 and its closest relatives in Arabidopsis thaliana. We generated an integrated framework of the response to JA, which spans from the activity of master and secondary regulatory transcription factors, through gene expression outputs and alternative splicing, to protein abundance changes, protein phosphorylation and chromatin remodelling. We integrated time-series transcriptome analysis with (phospho)proteomic data to reconstruct gene regulatory network models. These enabled us to predict previously unknown points of crosstalk of JA to other signalling pathways and to identify new components of the JA regulatory mechanism, which we validated through targeted mutant analysis. These results provide a comprehensive understanding of how a plant hormone remodels cellular functions and plant behaviour, the general principles of which provide a framework for analyses of cross-regulation between other hormone and stress signalling pathways.

A Reservoir of Pluripotent Phloem Cells Safeguards the Linear Developmental Trajectory of Protophloem Sieve Elements

Plant cells can change their identity based on positional information, a mechanism that confers developmental plasticity to plants. This ability, common to distinct multicellular organisms, is particularly relevant for plant phloem cells. Protophloem sieve elements (PSEs), one type of phloem conductive cells, act as the main organizers of the phloem pole, which comprises four distinct cell files organized in a conserved pattern. Here, we report how Arabidopsis roots generate a reservoir of meristematic phloem cells competent to swap their cell identities. Although PSE misspecification induces cell identity hybridism, the activity of RECEPTOR LIKE PROTEIN KINASE 2 (RPK2) by perceiving CLE45 peptide contributes to restrict PSE identity to the PSE position. By maintaining a spatiotemporal window when PSE and PSE-adjacent cells' identities are interchangeable, CLE45 signaling endows phloem cells with the competence to re-pattern a functional phloem pole when protophloem fails to form.