Establishing a reproducible approach to study cellular functions of plant cells with 3D bioprinting.

Capturing cell-to-cell signals in a three-dimensional (3D) environment is key to studying cellular functions. A major challenge in the current culturing methods is the lack of accurately capturing multicellular 3D environments. In this study, we established a framework for 3D bioprinting plant cells to study cell viability, cell division, and cell identity. We established long-term cell viability for bioprinted Arabidopsis and soybean cells. To analyze the generated large image datasets, we developed a high-throughput image analysis pipeline. Furthermore, we showed the cell cycle reentry of bioprinted cells for which the timing coincides with the induction of core cell cycle genes and regeneration-related genes, ultimately leading to microcallus formation. Last, the identity of bioprinted Arabidopsis root cells expressing endodermal markers was maintained for longer periods. The framework established here paves the way for a general use of 3D bioprinting for studying cellular reprogramming and cell cycle reentry toward tissue regeneration.

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.

tuxnet: a simple interface to process RNA sequencing data and infer gene regulatory networks.

Predicting gene regulatory networks (GRNs) from expression profiles is a common approach for identifying important biological regulators. Despite the increased use of inference methods, existing computational approaches often do not integrate RNA-sequencing data analysis, are not automated or are restricted to users with bioinformatics backgrounds. To address these limitations, we developed tuxnet, a user-friendly platform that can process raw RNA-sequencing data from any organism with an existing reference genome using a modified tuxedo pipeline (hisat 2 + cufflinks package) and infer GRNs from these processed data. tuxnet is implemented as a graphical user interface and can mine gene regulations, either by applying a dynamic Bayesian network (DBN) inference algorithm, genist, or a regression tree-based pipeline, rtp-star. We obtained time-course expression data of a PERIANTHIA (PAN) inducible line and inferred a GRN using genist to illustrate the use of tuxnet while gaining insight into the regulations downstream of the Arabidopsis root stem cell regulator PAN. Using rtp-star, we inferred the network of ATHB13, a downstream gene of PAN, for which we obtained wild-type and mutant expression profiles. Additionally, we generated two networks using temporal data from developmental leaf data and spatial data from root cell-type data to highlight the use of tuxnet to form new testable hypotheses from previously explored data. Our case studies feature the versatility of tuxnet when using different types of gene expression data to infer networks and its accessibility as a pipeline for non-bioinformaticians to analyze transcriptome data, predict causal regulations, assess network topology and identify key regulators.