Discriminative gene co-expression network analysis uncovers novel modules involved in the formation of phosphate deficiency-induced root hairs in Arabidopsis

Silveira S, Brommonschenkel SH, Fontes EPB. RLPredictiOme, a Machine Learning-Derived Method for High-Throughput Prediction of Plant Receptor-like Proteins, Reveals Novel Classes of Transmembrane Receptors

Cell surface receptors play essential roles in perceiving and processing external and internal signals at the cell surface of plants and animals. The receptor-like protein kinases (RLK) and receptor-like proteins (RLPs), two major classes of proteins with membrane receptor configuration, play a crucial role in plant development and disease defense. Although RLPs and RLKs share a similar single-pass transmembrane configuration, RLPs harbor short divergent C-terminal regions instead of the conserved kinase domain of RLKs. This RLP receptor structural design precludes sequence comparison algorithms from being used for high-throughput predictions of the RLP family in plant genomes, as has been extensively performed for RLK superfamily predictions. Here, we developed the RLPredictiOme, implemented with machine learning models in combination with Bayesian inference, capable of predicting RLP subfamilies in plant genomes. The ML models were simultaneously trained using six types of features, along with three stages to distinguish RLPs from non-RLPs (NRLPs), RLPs from RLKs, and classify new subfamilies of RLPs in plants. The ML models achieved high accuracy, precision, sensitivity, and specificity for predicting RLPs with relatively high probability ranging from 0.79 to 0.99. The prediction of the method was assessed with three datasets, two of which contained leucine-rich repeats (LRR)-RLPs from Arabidopsis and rice, and the last one consisted of the complete set of previously described Arabidopsis RLPs. In these validation tests, more than 90% of known RLPs were correctly predicted via RLPredictiOme. In addition to predicting previously characterized RLPs, RLPredictiOme uncovered new RLP subfamilies in the Arabidopsis genome. These include probable lipid transfer (PLT)-RLP, plastocyanin-like-RLP, ring finger-RLP, glycosyl-hydrolase-RLP, and glycerophosphoryldiester phosphodiesterase (GDPD, GDPDL)-RLP subfamilies, yet to be characterized. Compared to the only Arabidopsis GDPDL-RLK, molecular evolution studies confirmed that the ectodomain of GDPDL-RLPs might have undergone a purifying selection with a predominance of synonymous substitutions. Expression analyses revealed that predicted GDPGL-RLPs display a basal expression level and respond to developmental and biotic signals. The results of these biological assays indicate that these subfamily members have maintained functional domains during evolution and may play relevant roles in development and plant defense. Therefore, RLPredictiOme provides a framework for genome-wide surveys of the RLP superfamily as a foundation to rationalize functional studies of surface receptors and their relationships with different biological processes.

The Root Hair Development of Pectin Polygalacturonase PGX2 Activation Tagging Line in Response to Phosphate Deficiency

Pectin, cellulose, and hemicellulose constitute the primary cell wall in eudicots and function in multiple developmental processes in plants. Root hairs are outgrowths of specialized epidermal cells that absorb water and nutrients from the soil. Cell wall architecture influences root hair development, but how cell wall remodeling might enable enhanced root hair formation in response to phosphate (P) deficiency remains relatively unclear. Here, we found that POLYGALACTURONASE INVOLVED IN EXPANSION 2 (PGX2) functions in conditional root hair development. Under low P conditions, a PGX2 activation tagged line (PGX2^AT ) displays bubble-like root hairs and abnormal callose deposition and superoxide accumulation in roots. We found that the polar localization and trafficking of PIN2 are altered in PGX2^AT roots in response to P deficiency. We also found that actin filaments were less compact but more stable in PGX2^AT root hair cells and that actin filament skewness in PGX2^AT root hairs was recovered by treatment with 1-N-naphthylphthalamic acid (NPA), an auxin transport inhibitor. These results demonstrate that activation tagging of PGX2 affects cell wall remodeling, auxin signaling, and actin microfilament orientation, which may cooperatively regulate root hair development in response to P starvation.

Tackling Plant Phosphate Starvation by the Roots

Plant responses to phosphate deprivation encompass a wide range of strategies, varying from altering root system architecture, entering symbiotic interactions to excreting root exudates for phosphorous release, and recycling of internal phosphate. These processes are tightly controlled by a complex network of proteins that are specifically upregulated upon phosphate starvation. Although the different effects of phosphate starvation have been intensely studied, the full extent of its contribution to altered root system architecture remains unclear. In this review, we focus on the effect of phosphate starvation on the developmental processes that shape the plant root system and their underlying molecular pathways.

Sub-epidermal Expression of ENHANCER OF TRIPTYCHON AND CAPRICE1 and Its Role in Root Hair Formation Upon Pi Starvation

Root hair patterning is best studied in Arabidopsis thaliana. A pattern of root hair and non-root hair files is governed by a gene-regulatory network of activators and inhibitors. Under phosphate starvation conditions, extra root hairs are formed in non-root hair positions. This raises the question, whether and how this environmental stimulus is mediated by the known root hair gene network. In this study, we provide genetic and molecular data on the role of ETC1 in the phosphate starvation induced ectopic root hair formation. We show that the expression in the epidermis is irregular and reduced and that a new expression domain is induced in the sub-epidermis. By expressing ETC1 in the sub-epidermis, we show that this is sufficient to induce extra root hair formation in N-files. This suggests that the phosphate induced expressional switch from epidermal to epidermal plus sub-epidermal expression of ETC1 is one environmental input to the underlying patterning network.

Transcriptome analysis for identifying possible gene regulations during maize root emergence and formation at the initial growth stage

The root plays an important role during plant development and growth, i.e., the plant body maintenance, nutrient storage, absorption of water, oxygen and nutrient from the soil, and storage of water and carbohydrates, etc. The objective of this study was attempted to determine root-specific genes at the initial developmental stages of maize by using network-based transcriptome analysis. The raw data obtained using RNA-seq were filtered for quality control of the reads with the FASTQC tool, and the filtered reads were pre-proceed using the TRIMMOMATIC tool. The enriched BINs of the DEGs were detected using PageMan analysis with the ORA_FISHER statistical test, and genes were assigned to metabolic pathways by using the MapMan tool, which was also used for detecting transcription factors (TFs). For reconstruction of the co-expression network, we used the algorithm for the reconstruction of accurate cellular networks (ARACNE) in the R package, and then the reconstructed co-expression network was visualized using the Cytoscape tool. RNA-seq. was performed using maize shoots and roots at different developmental stages of root emergence (6-10 days after planting, VE) and 1 week after plant emergence (V2). A total of 1286 differentially expressed genes (DEGs) were detected in both tissues. Many DEGs involved in metabolic pathways exhibited altered mRNA levels between VE and V2. In addition, we observed gene expression changes for 113 transcription factors and found five enriched cis-regulatory elements in the 1-kb upstream regions of both DEGs. The network-based transcriptome analysis showed two modules as co-expressed gene clusters differentially expressed between the shoots and roots during plant development. The DEGs of one module exhibited gene expressional coherence in the maize root tips, suggesting that their functional relationships are associated with the initial developmental stage of the maize root. Finally, we confirmed reliable mRNA levels of the hub genes in the potential sub-network related to initial root development at the different developmental stages of VE, V2, and 2 weeks after plant emergence.

The Deubiquitinase OTU5 Regulates Root Responses to Phosphate Starvation

Phosphorus, taken up by plants as inorganic phosphate (Pi), is an essential but often growth-limiting mineral nutrient for plants. As part of an orchestrated response to improve its acquisition, insufficient Pi supply triggers alterations in root architecture and epidermal cell morphogenesis. Arabidopsis (Arabidopsis thaliana) mutants defective in the expression of the OVARIAN TUMOR DOMAIN-CONTAINING DEUBIQUITINATING ENZYME5 (OTU5) exhibited a constitutive Pi deficiency root phenotype, comprising the formation of long and dense root hairs and attenuated primary root growth. Quantitative protein profiling of otu5 and wild-type roots using the isobaric tag for relative and absolute quantification methodology revealed genotype- and Pi-dependent alterations in protein profiles. In otu5 plants, Pi starvation caused a short-root-hair phenotype and decreased abundance of a suite of Pi-responsive root hair-related proteins. Mutant plants also showed the accumulation of proteins involved in chromatin remodeling and altered distribution of reactive oxygen species along the root, which may be causative for the alterations in root hair morphogenesis. The root hair phenotype of otu5 was synergistic to that of actin-related protein6 (arp6), harboring a mutation in the SWR1 chromatin-remodeling complex. Genetic analysis of otu5/arp6 double mutants suggests independent but functionally related roles of the two proteins in chromatin organization. The root hair phenotype of otu5 is not caused by a general up-regulation of the Pi starvation response, indicating that OTU5 acts downstream of or interacts with Pi signaling. It is concluded that OTU5 is involved in the interpretation of environmental information, probably by altering chromatin organization and maintaining redox homeostasis.

Molecular networks orchestrating plant cell growth

Plant cell growth can broadly be categorized into either diffuse or tip growth. Here we compare gene regulatory networks (GRNs) controlling growth of hypocotyls and root hairs as examples for diffuse and tip growth, respectively. Accumulating evidence shows that GRNs in both cell types are multi-layered in structure and fine-tuned by transcriptional and post-translational mechanisms. We discuss how these GRNs regulate the expression of proteins controlling cell wall remodeling or other growth regulatory processes. Finally, we highlight how specific regulators within GRNs adjust plant cell growth in response to variable environmental conditions.