Coexpression-based clustering of Arabidopsis root genes predicts functional modules in early phosphate deficiency signaling

Transcriptomic Response of White Lupin Roots to Short-Term Sucrose Treatment

White lupin (Lupinus albus) has become a model plant for understanding plant adaptations to phosphorus (P) and iron (Fe) deficiency, two major limiting factors for plant productivity. In response to both nutrient deficiencies, white lupin forms cluster roots, bottle-brush-like root structures that aid in P and Fe acquisition from soil. While the cluster root function is well-studied, not much is known about the signaling pathways involved in sensing and responding to a P and Fe deficiency. Sucrose has been identified as a long-distance signal sent in increased concentrations from shoot to root in response to both a P and Fe deficiency. Thus, sucrose plays a dual role both as a signal and as a major source of energy for the root. To unravel the responses to sucrose as a signal, we performed an Illumina paired-end cDNA sequencing of white lupin roots treated with sucrose for 20, 40 or 80 min, compared to untreated controls (0 min). We identified 634 up-regulated and 956 down-regulated genes in response to sucrose. Twenty minutes of sucrose treatment showed the most responses, with the ethylene-activated signaling pathway as the most enriched Gene Ontology (GO) category. The number of up-regulated genes decreased at 40 min and 80 min, and protein dephosphorylation became the most enriched category. Taken together, our findings indicate active responses to sucrose as a signal at 20 min after a sucrose addition, but fewer responses and a potential resetting of signal transduction pathways by the dephosphorylation of proteins at 40 and 80 min.

Comparative transcriptomics pinpoints conserved and specific transcriptional responses to phosphorus limitation

Translating biological knowledge from Arabidopsis to crop species is important to advance agriculture and secure food production in the face of dwindling fertilizer resources and biotic and abiotic stresses. However, it is often not trivial to identify functional homologs (orthologs) of Arabidopsis genes in crops. Combining sequence and expression data can improve the correct prediction of orthologs. Here, we conducted a large-scale RNA sequencing based transcriptomics study of Arabidopsis, Medicago, Brachypodium, and Setaria grown side-by-side in phosphorus (P)-sufficient and P-limited conditions to generate comparable transcriptomics datasets. Comparison of top 200 P-limitation-induced genes in Arabidopsis revealed that ~80% of these genes have identifiable close homologs in the other three species but only ~50% retain their P-limitation response in the legume and grasses. Most of the hallmark genes of the P-starvation response were found conserved in all four species. This study reveals many known, novel, unannotated, conserved, and species-specific forms of regulation of the transcriptional P-starvation response. Identification and experimental verification of expressologs by independent RT-qPCR for P-limitation marker genes in Prunus showed the usefulness of comparative transcriptomics in pinpointing the functional orthologs in diverse crop species. This study provides an unprecedented resource for functional genomics and translational research to create P-efficient crops.

GLABRA2 transcription factor integrates arsenic tolerance with epidermal cell fate determination

Arsenic poses a global threat to living organisms, compromising crop security and yield. Limited understanding of the transcriptional network integrating arsenic-tolerance mechanisms with plant developmental responses hinders the development of strategies against this toxic metalloid. Here, we conducted a high-throughput yeast one-hybrid assay using as baits the promoter region from the arsenic-inducible genes ARQ1 and ASK18 from Arabidopsis thaliana, coupled with a transcriptomic analysis, to uncover novel transcriptional regulators of the arsenic response. We identified the GLABRA2 (GL2) transcription factor as a novel regulator of arsenic tolerance, revealing a wider regulatory role beyond its established function as a repressor of root hair formation. Furthermore, we found that ANTHOCYANINLESS2 (ANL2), a GL2 subfamily member, acts redundantly with this transcription factor in the regulation of arsenic signaling. Both transcription factors act as repressors of arsenic response. gl2 and anl2 mutants exhibit enhanced tolerance and reduced arsenic accumulation. Transcriptional analysis in the gl2 mutant unveils potential regulators of arsenic tolerance. These findings highlight GL2 and ANL2 as novel integrators of the arsenic response with developmental outcomes, offering insights for developing safer crops with reduced arsenic content and increased tolerance to this hazardous metalloid.

Arabidopsis PROTODERMAL FACTOR2 binds lysophosphatidylcholines and transcriptionally regulates phospholipid metabolism

Plant homeodomain leucine zipper IV (HD-Zip IV) transcription factors (TFs) contain an evolutionarily conserved steroidogenic acute regulatory protein (StAR)-related lipid transfer (START) domain. While the START domain is required for TF activity, its presumed role as a lipid sensor is not clear. Here we used tandem affinity purification from Arabidopsis cell cultures to demonstrate that PROTODERMAL FACTOR2 (PDF2), a representative member that controls epidermal differentiation, recruits lysophosphatidylcholines (LysoPCs) in a START-dependent manner. Microscale thermophoresis assays confirmed that a missense mutation in a predicted ligand contact site reduces lysophospholipid binding. We additionally found that PDF2 acts as a transcriptional regulator of phospholipid- and phosphate (Pi) starvation-related genes and binds to a palindromic octamer with consensus to a Pi response element. Phospholipid homeostasis and elongation growth were altered in pdf2 mutants according to Pi availability. Cycloheximide chase experiments revealed a role for START in maintaining protein levels, and Pi starvation resulted in enhanced protein destabilization, suggesting a mechanism by which lipid binding controls TF activity. We propose that the START domain serves as a molecular sensor for membrane phospholipid status in the epidermis. Our data provide insights toward understanding how the lipid metabolome integrates Pi availability with gene expression.

Transcriptional and Metabolomic Analyses Reveal That GmESR1 Increases Soybean Seed Protein Content Through the Phenylpropanoid Biosynthesis Pathway

Soybeans are an economically vital food crop, which is employed as a key source of oil and plant protein globally. This study identified an EREBP-type transcription factor, GmESR1 (Enhance of Shot Regeneration). GmESR1 overexpression has been observed to significantly increase seed protein content. Furthermore, the molecular mechanism by which GmESR1 affects protein accumulation through transcriptome and metabolomics was also identified. The transcriptomic and metabolomic analyses identified 95 differentially expressed genes and 83 differentially abundant metabolites during the seed mid-maturity stage. Co-analysis strategies revealed that GmESR1 overexpression inhibited the biosynthesis of lignin, cellulose, hemicellulose, and pectin via the phenylpropane biosynthetic pathway, thereby redistributing biomass within cells. The key genes and metabolites impacted by this biochemical process included Gm4CL-like, GmCCR, Syringin, and Coniferin. Moreover, it was also found that GmESR1 binds to (AATATTATCATTAAGTACGGAC) during seed development and inhibits the transcription of GmCCR. GmESR1 overexpression also enhanced sucrose transporter gene expression during seed development and increased the sucrose transport rate. These results offer new insight into the molecular mechanisms whereby GmESR1 increases protein levels within soybean seeds, guiding future molecular-assisted breeding efforts aimed at establishing high-protein soybean varieties.

The centrality of redox regulation and sensing of reactive oxygen species in abiotic and biotic stress acclimatization

During land plant evolution, the number of genes encoding for components of the thiol redox regulatory network and the generator systems of reactive oxygen species (ROS) expanded, tentatively indicating that they have a role in tailored environmental acclimatization. This hypothesis has been validated both experimentally and theoretically during the last few decades. Recent developments of dynamic redox-sensitive GFP (roGFP)-based in vivo sensors for H2O2 and the redox potential of the glutathione pool have paved the way for dissecting the kinetics changes that occur in these crucial parameters in response to environmental stressors. The versatile cellular redox sensory and response regulatory system monitors alterations in redox metabolism and controls the activity of redox target proteins, and thereby affects most, if not all, cellular processes ranging from transcription to translation and metabolism. This review uses examples to describe the role of the redox- and ROS-dependent regulatory network in realising the appropriate responses to diverse environmental stresses. The selected case studies concern different environmental challenges, namely excess excitation energy, the heavy metal cadmium and the metalloid arsenic, nitrogen or phosphate shortages as examples for nutrient deficiency, wounding, and nematode infestation. Each challenge affects the redox-regulatory and ROS network, but our present state of knowledge also points toward pressing questions that remain open in relation to the translation of redox regulation to environmental acclimatization.

MicroRNA1432 regulates rice drought stress tolerance by targeting the CALMODULIN-LIKE2 gene

Due to climate change, drought has become a major threat to rice (Oryza sativa L.) growth and yield worldwide. Understanding the genetic basis of drought tolerance in rice is therefore of great importance. Here, we identified a microRNA, miR1432, which regulates rice drought tolerance by targeting the CALMODULIN-LIKE2 (OsCaML2) gene. Mutation of MIR1432 or suppression of miR1432 expression significantly impaired seed germination and seedling growth under drought-stress conditions. Molecular analysis demonstrated that miR1432 affected rice drought tolerance by directly targeting OsCaML2, which encodes an EF-hand chiral calcium-binding protein. Overexpression of a miR1432-resistant form of OsCaML2 (OEmCaML2) phenocopied the mir1432 mutant and miR1432 suppression plants. Furthermore, the suppression of miR1432 severely affected the expression of genes involved in responses to stimulation, metabolism and signal transduction, especially the mitogen-activated protein kinase (MAPK) pathway and hormone transduction pathway in rice under drought stress. Thus, our findings show that the miR1432-OsCaML2 module plays an important role in the regulation of rice drought tolerance, suggesting its potential utilization in developing molecular breeding strategies that improve crop drought tolerance.

Phosphate starvation regulates cellulose synthesis to modify root growth

In the model plant Arabidopsis (Arabidopsis thaliana), the absence of the essential macro-nutrient phosphate reduces primary root growth through decreased cell division and elongation, requiring alterations to the polysaccharide-rich cell wall surrounding the cells. Despite its importance, the regulation of cell wall synthesis in response to low phosphate levels is not well understood. In this study, we show that plants increase cellulose synthesis in roots under limiting phosphate conditions, which leads to changes in the thickness and structure of the cell wall. These changes contribute to the reduced growth of primary roots in low-phosphate conditions. Furthermore, we found that the cellulose synthase complex (CSC) activity at the plasma membrane increases during phosphate deficiency. Moreover, we show that this increase in the activity of the CSC is likely due to alterations in the phosphorylation status of cellulose synthases in low-phosphate conditions. Specifically, phosphorylation of CELLULOSE SYNTHASE 1 (CESA1) at the S688 site decreases in low-phosphate conditions. Phosphomimic versions of CESA1 with an S688E mutation showed significantly reduced cellulose induction and primary root length changes in low-phosphate conditions. Protein structure modeling suggests that the phosphorylation status of S688 in CESA1 could play a role in stabilizing and activating the CSC. This mechanistic understanding of root growth regulation under limiting phosphate conditions provides potential strategies for changing root responses to soil phosphate content.

Alkalinity modulates a unique suite of genes to recalibrate growth and pH homeostasis

Alkaline soils pose a conglomerate of constraints to plants, restricting the growth and fitness of non-adapted species in habitats with low active proton concentrations. To thrive under such conditions, plants have to compensate for a potential increase in cytosolic pH and restricted softening of the cell wall to invigorate cell elongation in a proton-depleted environment. To discern mechanisms that aid in the adaptation to external pH, we grew plants on media with pH values ranging from 5.5 to 8.5. Growth was severely restricted above pH 6.5 and associated with decreasing chlorophyll levels at alkaline pH. Bicarbonate treatment worsened plant performance, suggesting effects that differ from those exerted by pH as such. Transcriptional profiling of roots subjected to short-term transfer from optimal (pH 5.5) to alkaline (pH 7.5) media unveiled a large set of differentially expressed genes that were partially congruent with genes affected by low pH, bicarbonate, and nitrate, but showed only a very small overlap with genes responsive to the availability of iron. Further analysis of selected genes disclosed pronounced responsiveness of their expression over a wide range of external pH values. Alkalinity altered the expression of various proton/anion co-transporters, possibly to recalibrate cellular proton homeostasis. Co-expression analysis of pH-responsive genes identified a module of genes encoding proteins with putative functions in the regulation of root growth, which appears to be conserved in plants subjected to low pH or bicarbonate. Our analysis provides an inventory of pH-sensitive genes and allows comprehensive insights into processes that are orchestrated by external pH.

Genome-wide identification of the WRKY gene family in Camellia oleifera and expression analysis under phosphorus deficiency

Camellia oleifera Abel. is an economically important woody edible-oil species that is mainly cultivated in hilly areas of South China. The phosphorus (P) deficiency in the acidic soils poses severe challenges for the growth and productivity of C. oleifera. WRKY transcription factors (TFs) have been proven to play important roles in biological processes and plant responses to various biotic/abiotic stresses, including P deficiency tolerance. In this study, 89 WRKY proteins with conserved domain were identified from the C. oleifera diploid genome and divided into three groups, with group II further classified into five subgroups based on the phylogenetic relationships. WRKY variants and mutations were detected in the gene structure and conserved motifs of CoWRKYs. Segmental duplication events were considered as the primary driver in the expanding process of WRKY gene family in C. oleifera. Based on transcriptomic analysis of two C. oleifera varieties characterized with different P deficiency tolerances, 32 CoWRKY genes exhibited divergent expression patterns in response to P deficiency stress. qRT-PCR analysis demonstrated that CoWRKY11, -14, -20, -29 and -56 had higher positive impact on P-efficient CL40 variety compared with P-inefficient CL3 variety. Similar expression trends of these CoWRKY genes were further observed under P deficiency with longer treatment period of 120d. The result indicated the expression sensitivity of CoWRKYs on the P-efficient variety and the C. oleifera cultivar specificity on the P deficiency tolerance. Tissue expression difference showed CoWRKYs may play a crucial role in the transportation and recycling P in leaves by affecting diverse metabolic pathways. The available evidences in the study conclusively shed light on the evolution of the CoWRKY genes in C. oleifera genome and provided a valuable resource for further investigation of functional characterization of WRKY genes involved to enhance the P deficiency tolerance in C. oleifera.

Phosphate-deprivation and damage signalling by extracellular ATP

Phosphate deprivation compromises plant productivity and modulates immunity. DAMP signalling by extracellular ATP (eATP) could be compromised under phosphate deprivation by the lowered production of cytosolic ATP and the need to salvage eATP as a nutritional phosphate source. Phosphate-starved roots of Arabidopsis can still sense eATP, indicating robustness in receptor function. However, the resultant cytosolic free Ca²⁺ signature is impaired, indicating modulation of downstream components. This perspective on DAMP signalling by extracellular ATP (eATP) addresses the salvage of eATP under phosphate deprivation and its promotion of immunity, how Ca²⁺ signals are generated and how the Ca²⁺ signalling pathway could be overcome to allow beneficial fungal root colonization to fulfill phosphate demands. Safe passage for an endophytic fungus allowing root colonization could be achieved by its down-regulation of the Ca²⁺ channels that act downstream of the eATP receptors and by also preventing ROS accumulation, thus further impairing DAMP signalling.

Combined Approach of GWAS and Phylogenetic Analyses to Identify New Candidate Genes That Participate in Arabidopsis thaliana Primary Root Development Using Cellular Measurements and Primary Root Length

Genome-wide association studies (GWAS) have allowed the identification of different loci associated with primary root (PR) growth, and Arabidopsis is an excellent model for these studies. The PR length is controlled by cell proliferation, elongation, and differentiation; however, the specific contribution of proliferation and differentiation in the control of PR growth is still poorly studied. To this end, we analyzed 124 accessions and used a GWAS approach to identify potential causal genomic regions related to four traits: PR length, growth rate, cell proliferation and cell differentiation. Twenty-three genes and five statistically significant SNPs were identified. The SNP with the highest score mapped to the fifth exon of NAC048 and this change makes a missense variant in only 33.3% of the accessions with a large PR, compared with the accessions with a short PR length. Moreover, we detected five more SNPs in this gene and in NAC3 that allow us to discover closely related accessions according to the phylogenetic tree analysis. We also found that the association between genetic variants among the 18 genes with the highest scores in our GWAS and the phenotypic classes into which we divided our accessions are not straightforward and likely follow historical patterns.

Noncoding RNAs responsive to nitric oxide and their protein-coding gene targets shed light on root hair formation in Arabidopsis thaliana

An overview of the total Arabidopsis thaliana transcriptome, described previously by our research group, pointed some noncoding RNA (ncRNA) as participants in the restoration of hair-root phenotype in A. thaliana rhd6 mutants, leading us to a deeper investigation. A transcriptional gene expression profiling of seedling roots was performed aiming to identify ncRNA responsive to nitric oxide (GSNO) and auxin (IAA), and their involvement in root hair formation in the rhd6 null mutant. We identified 3,631 ncRNAs, including new ones, in A. thaliana and differential expression (DE) analysis between the following: 1) GSNO-treated rhd6 vs. untreated rhd6, 2) IAA-treated rhd6 vs. untreated rhd6, 3) GSNO-treated rhd6 vs. IAA-treated rhd6, and 4) WS-2 vs. untreated rhd6 detected the greatest number of DE genes in GSNO-treated rhd6. We detected hundreds of in silico interactions among ncRNA and protein-coding genes (PCGs), highlighting MIR5658 and MIR171 precursors highly upregulated in GSNO-treated rhd6 and wild type, respectively. Those ncRNA interact with many DE PCGs involved in hormone signaling, cell wall development, transcription factors, and root hair formation, becoming candidate genes in cell wall modulation and restoration of root hair phenotype by GSNO treatment. Our data shed light on how GSNO modulates ncRNA and their PCG targets in A. thaliana root hair formation.

Phosphate transporter PHT1;1 is a key determinant of phosphorus acquisition in Arabidopsis natural accessions

Phosphorus (P) is a mineral nutrient essential for plant growth and development, but most P in the soil is unavailable for plants. To understand the genetic basis of P acquisition regulation, we performed genome-wide association studies (GWASs) on a diversity panel of Arabidopsis (Arabidopsis thaliana). Two primary determinants of P acquisition were considered, namely, phosphate (Pi)-uptake activity and PHOSPHATE TRANSPORTER 1 (PHT1) protein abundance. Association mapping revealed a shared significant peak on chromosome 5 (Chr5) where the PHT1;1/2/3 genes reside, suggesting a connection between the regulation of Pi-uptake activity and PHT1 protein abundance. Genes encoding transcription factors, kinases, and a metalloprotease associated with both traits were also identified. Conditional GWAS followed by statistical analysis of genotype-dependent PHT1;1 expression and transcriptional activity assays revealed an epistatic interaction between PHT1;1 and MYB DOMAIN PROTEIN 52 (MYB52) on Chr1. Further, analyses of F1 hybrids generated by crossing two subgroups of natural accessions carrying specific PHT1;1- and MYB52-associated single nucleotide polymorphisms (SNPs) revealed strong effects of these variants on PHT1;1 expression and Pi uptake activity. Notably, the soil P contents in Arabidopsis habitats coincided with PHT1;1 haplotype, emphasizing how fine-tuned P acquisition activity through natural variants allows environmental adaptation. This study sheds light on the complex regulation of P acquisition and offers a framework to systematically assess the effectiveness of GWAS approaches in the study of quantitative traits.

Spatial origin of the extracellular ATP-induced cytosolic calcium signature in Arabidopsis thaliana roots: wave formation and variation with phosphate nutrition

Extracellular ATP (eATP) increases cytosolic free calcium ([Ca²⁺ ]_cyt ) as a specific second messenger 'signature' through the plasma membrane DORN1/P2K1 receptor. Previous studies revealed a biphasic signature in Arabidopsis thaliana roots that is altered by inorganic phosphate (Pi) deprivation. The relationship between the two phases of the signature and possible wave formation have been tested as a function of Pi nutrition. The bioluminescent aequorin and intensiometric GCaMP3 reporters were used to resolve the spatial origin of the eATP [Ca²⁺ ]_cyt signature in Arabidopsis root tips. Application of eATP only to the root apex allowed [Ca²⁺ ]_cyt wave resolution without the confounding effects of eATP delivery by superfusion. The first apical millimetre of the root generates the first [Ca²⁺ ]_cyt increase by eATP, regardless of nutritional status. The second increase occurs sub-apically in the root hair zone, has some autonomy and is significantly reduced in Pi-starved roots. A significant component of the Pi-replete signature does not require DORN1/P2K1, but Pi-starved roots appear to have an absolute requirement for that receptor. Application of eATP specifically to the root apex provides evidence for cell-to-cell propagation of a [Ca²⁺ ]_cyt wave that diminishes sub-apically. The apex maintains a robust [Ca²⁺ ]_cyt increase (even under Pi starvation) that is the basis of a propagative wave, with implications for the ability of the root's eATP signalling systems to signal systemically. Partial autonomy of the sub-apical region may be relevant to the perception of eATP from microbes. eATP-induced [Ca²⁺ ]_cyt increase may not have always have an obligate requirement for DORN1/P2K1.

Arabinogalactan Proteins: Focus on the Role in Cellulose Synthesis and Deposition during Plant Cell Wall Biogenesis

Arabinogalactan proteins (AGPs) belong to a family of glycoproteins that are widely present in plants. AGPs are mostly composed of a protein backbone decorated with complex carbohydrate side chains and are usually anchored to the plasma membrane or secreted extracellularly. A trickle of compelling biochemical and genetic evidence has demonstrated that AGPs make exciting candidates for a multitude of vital activities related to plant growth and development. However, because of the diversity of AGPs, functional redundancy of AGP family members, and blunt-force research tools, the precise functions of AGPs and their mechanisms of action remain elusive. In this review, we put together the current knowledge about the characteristics, classification, and identification of AGPs and make a summary of the biological functions of AGPs in multiple phases of plant reproduction and developmental processes. In addition, we especially discuss deeply the potential mechanisms for AGP action in different biological processes via their impacts on cellulose synthesis and deposition based on previous studies. Particularly, five hypothetical models that may explain the AGP involvement in cellulose synthesis and deposition during plant cell wall biogenesis are proposed. AGPs open a new avenue for understanding cellulose synthesis and deposition in plants.

Hypoxia-Induced Aquaporins and Regulation of Redox Homeostasis by a Trans-Plasma Membrane Electron Transport System in Maize Roots

In plants, flooding-induced oxygen deficiency causes severe stress, leading to growth reduction and yield loss. It is therefore important to understand the molecular mechanisms for adaptation to hypoxia. Aquaporins at the plasma membrane play a crucial role in water uptake. However, their role during hypoxia and membrane redox changes is still not fully understood. The influence of 24 h hypoxia induction on hydroponically grown maize (Zea mays L.) was investigated using an oil-based setup. Analyses of physiological parameters revealed typical flooding symptoms such as increased ethylene and H₂O₂ levels, an increased alcohol dehydrogenase activity, and an increased redox activity at the plasma membrane along with decreased oxygen of the medium. Transcriptomic analysis and shotgun proteomics of plasma membranes and soluble fractions were performed to determine alterations in maize roots. RNA-sequencing data confirmed the upregulation of genes involved in anaerobic metabolism, biosynthesis of the phytohormone ethylene, and its receptors. Transcripts of several antioxidative systems and other oxidoreductases were regulated. Mass spectrometry analysis of the plasma membrane proteome revealed alterations in redox systems and an increased abundance of aquaporins. Here, we discuss the importance of plasma membrane aquaporins and redox systems in hypoxia stress response, including the regulation of plant growth and redox homeostasis.

Ethylene Response Factor109 Attunes Immunity, Photosynthesis, and Iron Homeostasis in Arabidopsis Leaves

Iron (Fe) is an essential micronutrient element for all organisms including plants. Chlorosis of young leaves is a common symptom of Fe deficiency, reducing the efficiency of photosynthesis, and, ultimately, crop yield. Previous research revealed strong responsiveness of the putative key transcription factor ERF109 to the Fe regime. To elucidate the possible role of ERF109 in leaf Fe homeostasis and photosynthesis, we subjected Arabidopsis thaliana erf109 knockout lines and Col-0 wild-type plants to transcriptome profiling via RNA-seq. The transcriptome profile of Fe-sufficient erf109 leaves showed a 71% overlap with Fe-deficient Col-0 plants. On the other hand, genes that were differentially expressed between Fe-deficient and Fe-sufficient Col-0 plants remained unchanged in erf109 plants under conditions of Fe deficiency. Mutations in ERF109 increased the expression of the clade Ib bHLH proteins bHLH38, bHLH39, bHLH101, the nicotianamine synthase NAS4, and the Fe storage gene FER1. Moreover, mutations in ERF109 led to significant down-regulation of defense genes, including CML37, WRKY40, ERF13, and EXO70B2. Leaves of erf109 exhibited increased Fe levels under both Fe-sufficient and Fe-deficient conditions. Reduced Fv/Fm and Soil Plant Analysis Development (SPAD) values in erf109 lines under Fe deficiency indicate curtailed ability of photosynthesis relative to the wild-type. Our findings suggest that ERF109 is a negative regulator of the leaf response to Fe deficiency. It further appears that the function of ERF109 in the Fe response is critical for regulating pathogen defense and photosynthetic efficiency. Taken together, our study reveals a novel function of ERF109 and provides a systematic perspective on the intertwining of the immunity regulatory network and cellular Fe homeostasis.

An improved repertoire of splicing variants and their potential roles in Arabidopsis photomorphogenic development

Background

Light switches on the photomorphogenic development of young plant seedlings, allowing young seedlings to acquire photosynthetic capacities and gain survival fitness. Light regulates gene expression at all levels of the central dogma, including alternative splicing (AS) during the photomorphogenic development. However, accurate determination of full-length (FL) splicing variants has been greatly hampered by short-read RNA sequencing technologies.

Result

In this study, we adopt PacBio isoform sequencing (Iso-seq) to overcome the limitation of the short-read RNA-seq technologies. Normalized cDNA libraries used for Iso-seq allows for comprehensive and effective identification of FL AS variants. Our analyses reveal more than 30,000 splicing variant models from approximately 16,500 gene loci and additionally identify approximately 700 previously unannotated genes. Among the variants, approximately 12,000 represent new gene models. Intron retention (IR) is the most frequently observed form of variants, and many IR-containing AS variants show evidence of engagement in translation. Our study reveals the formation of heterodimers of transcription factors composed of annotated and IR-containing AS variants. Moreover, transgenic plants overexpressing the IR forms of two B-BOX DOMAIN PROTEINs exhibits light-hypersensitive phenotypes, suggesting their regulatory roles in modulating optimal light responses.

Conclusions

This study provides an accurate and comprehensive portrait of full-length transcript isoforms and experimentally confirms the presence of de novo synthesized AS variants that impose regulatory functions in photomorphogenic development in Arabidopsis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02620-2.

Comprehensive In Silico Analysis and Transcriptional Profiles Highlight the Importance of Mitochondrial Dicarboxylate Carriers (DICs) on Hypoxia Response in Both Arabidopsis thaliana and Eucalyptus grandis

Plant dicarboxylate carriers (DICs) transport a wide range of dicarboxylates across the mitochondrial inner membrane. The Arabidopsis thalianaDIC family is composed of three genes (AtDIC1, 2 and 3), whereas two genes (EgDIC1 and EgDIC2) have been retrieved in Eucalyptus grandis. Here, by combining in silico and in planta analyses, we provide evidence that DICs are partially redundant, important in plant adaptation to environmental stresses and part of a low-oxygen response in both species. AtDIC1 and AtDIC2 are present in most plant species and have very similar gene structure, developmental expression patterns and absolute expression across natural Arabidopsis accessions. In contrast, AtDIC3 seems to be an early genome acquisition found in Brassicaceae and shows relatively low (or no) expression across these accessions. In silico analysis revealed that both AtDICs and EgDICs are highly responsive to stresses, especially to cold and submergence, while their promoters are enriched for stress-responsive transcription factors binding sites. The expression of AtDIC1 and AtDIC2 is highly correlated across natural accessions and in response to stresses, while no correlation was found for AtDIC3. Gene ontology enrichment analysis suggests a role for AtDIC1 and AtDIC2 in response to hypoxia, and for AtDIC3 in phosphate starvation. Accordingly, the investigated genes are induced by submergence stress in A. thaliana and E. grandis while AtDIC2 overexpression improved seedling survival to submergence. Interestingly, the induction of AtDIC1 and AtDIC2 is abrogated in the erfVII mutant that is devoid of plant oxygen sensing, suggesting that these genes are part of a conserved hypoxia response in Arabidopsis.

The versatility of plant organic acid metabolism in leaves is underpinned by mitochondrial malate-citrate exchange

Malate and citrate underpin the characteristic flexibility of central plant metabolism by linking mitochondrial respiratory metabolism with cytosolic biosynthetic pathways. However, the identity of mitochondrial carrier proteins that influence both processes has remained elusive. Here we show by a systems approach that DICARBOXYLATE CARRIER 2 (DIC2) facilitates mitochondrial malate-citrate exchange in vivo in Arabidopsis thaliana. DIC2 knockout (dic2-1) retards growth of vegetative tissues. In vitro and in organello analyses demonstrate that DIC2 preferentially imports malate against citrate export, which is consistent with altered malate and citrate utilization in response to prolonged darkness of dic2-1 plants or a sudden shift to darkness of dic2-1 leaves. Furthermore, isotopic glucose tracing reveals a reduced flux towards citrate in dic2-1, which results in a metabolic diversion towards amino acid synthesis. These observations reveal the physiological function of DIC2 in mediating the flow of malate and citrate between the mitochondrial matrix and other cell compartments.

Transcriptome provides potential insights into how calcium affects the formation of stone cell in Pyrus

Background

The content of stone cells in pears has a great influence on taste. Stone cells are formed by the accumulation of lignin. The treatment of exogenous calcium can affect the lignin synthesis, but this Ca-mediated mechanism is still unclear. In this study, the author performed a comparative transcriptomic analysis of callus of pears (Pyrus x bretschneideri) treated with calcium nitrate Ca (NO₃)₂ to investigate the role of calcium in lignin synthesis.

Results

There were 2889 differentially expressed genes (DEGs) detected between the Control and Ca (NO₃)₂ treatment in total. Among these 2889 DEGs, not only a large number of genes related to Ca single were found, but also many genes were enriched in secondary metabolic pathway, especially in lignin synthesis. Most of them were up-regulated during the development of callus after Ca (NO₃)₂ treatment. In order to further explore how calcium nitrate treatment affects lignin synthesis, the author screened genes associated with transduction of calcium signal in DEGs, and finally found CAM, CML, CDPK, CBL and CIPK. Then the author identified the PbCML3 in pears and conducted relevant experiments finding the overexpression of PbCML3 would increase the content of pear stone cells, providing potential insights into how Ca treatment enhances the stone cell in pears.

Conclusions

Our deep analysis reveals the effects of exogenous calcium on calcium signal and lignin biosynthesis pathway. The function of PbCML3 on stone cells formation was verified in pear.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08161-5.

IRONMAN tunes responses to iron deficiency in concert with environmental pH

Iron (Fe) is an essential mineral element that governs the composition of natural plant communities and limits crop yield in agricultural ecosystems due to its extremely low availability in most soils, particularly at alkaline pH. To extract sufficient Fe from the soil under such conditions, some plants, including Arabidopsis (Arabidopsis thaliana), secrete Fe-mobilizing phenylpropanoids, which mobilize sparingly soluble Fe hydroxides by reduction and chelation. We show here that ectopic expression of the peptides IRONMAN (IMA1) and IMA2 improves growth on calcareous soil by inducing biosynthesis and secretion of the catecholic coumarin 7,8-dihydroxy-6-methoxycoumarin (fraxetin) via increased expression of MYB72 and SCOPOLETIN 8-HYDROXYLASE, a response that is strictly dependent on elevated environmental pH (pHe). By contrast, transcription of the cytochrome P450 family protein CYP82C4, catalyzing the subsequent hydroxylation of fraxetin to sideretin, which forms less stable complexes with iron, was strongly repressed under such conditions. We concluded that IMA peptides regulate processes supporting Fe uptake at both acidic and elevated pH by controlling gene expression upstream of or in concert with a putative pHe signal, adapting the plant to prevailing edaphic conditions. This regulatory pattern confers tolerance to calcareous soils by extending the pH range in which Fe can be efficiently absorbed from the soil. Our results further suggest that pHe calibrates the activities of components of the Fe deficiency response, accentuating processes that are most efficient under the prevailing conditions. Altering the expression of IMA peptides provides a route for generating plants adapted to calcareous soils.

Genome-wide analysis reveals widespread roles for RcREM genes in floral organ development in Rosa chinensis

Members of the REM (Reproductive Meristem) gene family are expressed primarily in reproductive meristems and floral organs. However, their evolution and their functional profiles in flower development remain poorly understood. Here, we performed genome-wide identification and evolutionary analysis of the REM gene family in Rosaceae. This family has been greatly expanded in rose (Rosa chinensis) compared to other species, primarily through tandem duplication. Expression analysis revealed that most RcREM genes are specifically expressed in reproductive organs and that their specific expression patterns are dramatically altered in rose plants with mutations affecting floral organs. Protein-protein interaction analysis indicated that RcREM14 interact with RcAP1 (one of the homology of A class genes in ABCDE model), highlighting the roles of RcREM genes in floral organ identity. Finally, co-expression network analysis indicated that RcREM genes are co-expressed with a high proportion of key genes that regulate flowering time, floral organ development, and cell proliferation and expansion in R. chinensis.

Comparative transcriptomic analysis uncovers conserved pathways involved in adventitious root formation in poplar

Cutting propagation is widely used in establishing poplar plantations, and this approach requires efficient adventitious root (AR) forming capacities. Although poplar species are considered to form roots easily, interspecific variations in AR formation are still observed. To better understand the gene regulatory network underlying the conserved modified pathways that are essential for AR formation in poplar species, comparative transcriptomic approaches were applied to identify the conserved common genes that were differentially expressed during the AR formation processes in two poplar species (Populus × euramericana and P. simonii) in woody plant medium (WPM). A total of 2146 genes were identified as conserved genes that shared similar gene expression profiles in at least one comparison. These conserved genes were enriched in diverse hormone signaling pathways, as well as the mitogen-associated protein kinase (MAPK) signaling pathway, suggesting an important role for signaling transduction in coordinating external stimuli and endogenous physiological status during AR regulation in poplar. Furthermore, the co-expression network analysis of conserved genes allowed identification of several co-expressed modules (CM) that are co-expressed with distinct biological functions, for instance, CM1 was enriched in defense response and hormone signaling, CM2 and CM3 were overrepresented in defense response-related pathways and for cell cycle, respectively. These results suggest that the AR formation processes in poplar were finely tuned at the transcriptomic level by integrating multiple biological processes essential for AR formation. Our results suggest conserved machinery for AR formation in poplar and generated informative gene co-expression networks that describe the basis of AR formation in these species.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01054-7.

Single-cell transcriptome atlas and chromatin accessibility landscape reveal differentiation trajectories in the rice root

Root development relies on the establishment of meristematic tissues that give rise to distinct cell types that differentiate across defined temporal and spatial gradients. Dissection of the developmental trajectories and the transcriptional networks that underlie them could aid understanding of the function of the root apical meristem in both dicots and monocots. Here, we present a single-cell RNA (scRNA) sequencing and chromatin accessibility survey of rice radicles. By temporal profiling of individual root tip cells we reconstruct continuous developmental trajectories of epidermal cells and ground tissues, and elucidate regulatory networks underlying cell fate determination in these cell lineages. We further identify characteristic processes, transcriptome profiles, and marker genes for these cell types and reveal conserved and divergent root developmental pathways between dicots and monocots. Finally, we demonstrate the potential of the platform for functional genetic studies by using spatiotemporal modeling to identify a rice root meristematic mutant from a cell-specific gene cohort.

Metals and methylotrophy: Via global gene expression studies

A number of minerals, such as copper, cobalt, and rare earth elements (REE), are essential modulators of microbial one-carbon metabolism. This chapter provides an overview of the gene expression study design and analysis protocols for uncovering REE-induced changes in methylotrophic bacteria. By interrogating relationships and differences in total gene expression induced by mineral micronutrients, a deeper understanding of gene regulation at a systems scale can be gained. With careful design and execution of RNA-sequencing experiments, thorough processing and assessment of read quality can be utilized to assess and adjust for possible biases. By ensuring only quality data are utilized in downstream processes, differential gene expression, overrepresented analyses, and gene-set enrichment analyses provide reliable and reproducible representation of pathways and functions which are being affected by changes in environmental conditions.

Arabinogalactan Proteins in Plant Roots - An Update on Possible Functions

Responsiveness to environmental conditions and developmental plasticity of root systems are crucial determinants of plant fitness. These processes are interconnected at a cellular level with cell wall properties and cell surface signaling, which involve arabinogalactan proteins (AGPs) as essential components. AGPs are cell-wall localized glycoproteins, often GPI-anchored, which participate in root functions at many levels. They are involved in cell expansion and differentiation, regulation of root growth, interactions with other organisms, and environmental response. Due to the complexity of cell wall functional and regulatory networks, and despite the large amount of experimental data, the exact molecular mechanisms of AGP-action are still largely unknown. This dynamically evolving field of root biology is summarized in the present review.

Molecular basis for neofunctionalization of duplicated E3 ubiquitin ligases underlying adaptation to drought tolerance in Arabidopsis thaliana

Multigene families in plants expanded from ancestral genes via gene duplication mechanisms constitute a significant fraction of the coding genome. Although most duplicated genes are lost over time, many are retained in the genome. Clusters of tandemly arrayed genes are commonly found in the plant genome where they can promote expansion of gene families. In the present study, promoter fusion to the GUS reporter gene was used to examine the promoter architecture of duplicated E3 ligase genes that are part of group C in the Arabidopsis thaliana ATL family. Acquisition of gene expression by AtATL78, possibly generated from defective AtATL81 expression, is described. AtATL78 expression was purportedly enhanced by insertion of a TATA box within the core promoter region after a short tandem duplication that occurred during evolution of Brassicaceae lineages. This gene is associated with an adaptation to drought tolerance of A. thaliana. These findings also suggest duplicated genes could serve as a reservoir of tacit genetic information, and expression of these duplicated genes is activated upon acquisition of core promoter sequences. Remarkably, drought transcriptome profiling in response to rehydration suggests that ATL78-dependent gene expression predominantly affects genes with root-specific activities.

Phosphate Deprivation Can Impair Mechano-Stimulated Cytosolic Free Calcium Elevation in Arabidopsis Roots

The root tip responds to mechanical stimulation with a transient increase in cytosolic free calcium as a possible second messenger. Although the root tip will grow through a heterogeneous soil nutrient supply, little is known of the consequence of nutrient deprivation for such signalling. Here, the effect of inorganic phosphate deprivation on the root’s mechano-stimulated cytosolic free calcium increase is investigated. Arabidopsisthaliana (cytosolically expressing aequorin as a bioluminescent free calcium reporter) is grown in zero or full phosphate conditions, then roots or root tips are mechanically stimulated. Plants also are grown vertically on a solid medium so their root skewing angle (deviation from vertical) can be determined as an output of mechanical stimulation. Phosphate starvation results in significantly impaired cytosolic free calcium elevation in both root tips and whole excised roots. Phosphate-starved roots sustain a significantly lower root skewing angle than phosphate-replete roots. These results suggest that phosphate starvation causes a dampening of the root mechano-signalling system that could have consequences for growth in hardened, compacted soils.

A glutathione-dependent control of the indole butyric acid pathway supports Arabidopsis root system adaptation to phosphate deprivation

Root system architecture results from a highly plastic developmental process to adapt to environmental conditions. In particular, the development of lateral roots and root hair growth are constantly optimized to the rhizosphere properties, including biotic and abiotic constraints. The development of the root system is tightly controlled by auxin, the driving morphogenic hormone in plants. Glutathione, a major thiol redox regulator, is also critical for root development but its interplay with auxin is scarcely understood. Previous work showed that glutathione deficiency does not alter root responses to indole acetic acid (IAA), the main active auxin in plants. Because indole butyric acid (IBA), another endogenous auxinic compound, is an important source of IAA for the control of root development, we investigated the crosstalk between glutathione and IBA during root development. We show that glutathione deficiency alters lateral roots and root hair responses to exogenous IBA but not IAA. Detailed genetic analyses suggest that glutathione regulates IBA homeostasis or conversion to IAA in the root cap. Finally, we show that both glutathione and IBA are required to trigger the root hair response to phosphate deprivation, suggesting an important role for this glutathione-dependent regulation of the auxin pathway in plant developmental adaptation to its environment.

Scion genotypes exert long distance control over rootstock transcriptome responses to low phosphate in grafted grapevine

BACKGROUND

Grafting is widely used in horticulture and rootstocks are known to modify scion growth and adaptation to soil conditions. However, the role of scion genotype in regulating rootstock development and functioning has remained largely unexplored. In this study, reciprocal grafts of two grapevine genotypes were produced as well as the corresponding homo-graft controls. These plants were subjected to a low phosphate (LP) treatment and transcriptome profiling by RNA sequencing was done on root samples collected 27 h after the onset of the LP treatment.

RESULTS

A set of transcripts responsive to the LP treatment in all scion/rootstock combinations was identified. Gene expression patterns associated with genetic variation in response to LP were identified by comparing the response of the two homo-grafts. In addition, the scion was shown to modify root transcriptome responses to LP in a rootstock dependent manner. A weighted gene co-expression network analysis identified modules of correlated genes; the analysis of the association of these modules with the phosphate treatment, and the scion and rootstock genotype identified potential hub genes.

CONCLUSIONS

This study provides insights into the response of grafted grapevine to phosphate supply and identifies potential shoot-to-root signals that could vary between different grapevine genotypes.

Genome-Wide Dissection of the CRF Gene Family in Brassica napus Indicates that BnaCRF8s Specifically Regulate Root Architecture and Phosphate Homeostasis against Phosphate Fluctuation in Plants

Phosphorus (P) is an essential macronutrient required for plant growth and development. The involvement of cytokinin response factors (CRFs) in phosphate (Pi) homeostasis and lateral root (LR) initiation in Arabidopsis has been revealed. However, little is known in oil crops. Here, we performed genome-wide dissection of the CRF family in Brassica napus to identify 44 members, which were evolutionally classified into 6 subgroups. Among them, four BnaCRF8 genes were strongly upregulated by P deprivation, and were selected to be further investigated. Time course qRT-PCR analyses showed that four BnaCRF8 genes were enhanced dramatically after 12 h of P stress. Analyses of the subcellular localization in tobacco leaves indicated that BnaA7.CRF8 and BnaC2.CRF8 were localized in the nucleus. The expression of BnaCRF8 genes had constant negative effects on primary root growth and LR initiation and growth, and it reduced Pi acquisition and plant growth in Arabidopsis. Moreover, the expression of Pi homeostasis-related genes was modulated in BnaA7.CRF8 overexpression plants. These results suggest that BnaCRF8 genes might negatively regulate root architecture and plant growth through transcriptional modification of Pi homeostasis-related components. Overall, this study suggests that upregulation of BnaCRF8 genes might be a smart adaptive strategy to cope with continuous Pi deficiency in the environment.

PEATmoss (Physcomitrella Expression Atlas Tool): a unified gene expression atlas for the model plant Physcomitrella patens

Physcomitrella patens is a bryophyte model plant that is often used to study plant evolution and development. Its resources are of great importance for comparative genomics and evo-devo approaches. However, expression data from Physcomitrella patens were so far generated using different gene annotation versions and three different platforms: CombiMatrix and NimbleGen expression microarrays and RNA sequencing. The currently available P. patens expression data are distributed across three tools with different visualization methods to access the data. Here, we introduce an interactive expression atlas, Physcomitrella Expression Atlas Tool (PEATmoss), that unifies publicly available expression data for P. patens and provides multiple visualization methods to query the data in a single web-based tool. Moreover, PEATmoss includes 35 expression experiments not previously available in any other expression atlas. To facilitate gene expression queries across different gene annotation versions, and to access P. patens annotations and related resources, a lookup database and web tool linked to PEATmoss was implemented. PEATmoss can be accessed at https://peatmoss.online.uni-marburg.de.

Low-Phosphate Chromatin Dynamics Predict a Cell Wall Remodeling Network in Rice Shoots

Phosphorus (P) is an essential plant macronutrient vital to fundamental metabolic processes. Plant-available P is low in most soils, making it a frequent limiter of growth. Declining P reserves for fertilizer production exacerbates this agricultural challenge. Plants modulate complex responses to fluctuating P levels via global transcriptional regulatory networks. Although chromatin structure plays a substantial role in controlling gene expression, the chromatin dynamics involved in regulating P homeostasis have not been determined. Here we define distinct chromatin states across the rice (Oryza sativa) genome by integrating multiple chromatin marks, including the H2A.Z histone variant, H3K4me3 modification, and nucleosome positioning. In response to P starvation, 40% of all protein-coding genes exhibit a transition from one chromatin state to another at their transcription start site. Several of these transitions are enriched in subsets of genes differentially expressed under P deficiency. The most prominent subset supports the presence of a coordinated signaling network that targets cell wall structure and is regulated in part via a decrease of H3K4me3 at transcription start sites. The P starvation-induced chromatin dynamics and correlated genes identified here will aid in enhancing P use efficiency in crop plants, benefitting global agriculture.

The polyadenylation factor FIP1 is important for plant development and root responses to abiotic stresses

Root development and its response to environmental changes is crucial for whole plant adaptation. These responses include changes in transcript levels. Here, we show that the alternative polyadenylation (APA) of mRNA is important for root development and responses. Mutations in FIP1, a component of polyadenylation machinery, affects plant development, cell division and elongation, and response to different abiotic stresses. Salt treatment increases the amount of poly(A) site usage within the coding region and 5' untranslated regions (5'-UTRs), and the lack of FIP1 activity reduces the poly(A) site usage within these non-canonical sites. Gene ontology analyses of transcripts displaying APA in response to salt show an enrichment in ABA signaling, and in the response to stresses such as salt or cadmium (Cd), among others. Root growth assays show that fip1-2 is more tolerant to salt but is hypersensitive to ABA or Cd. Our data indicate that FIP1-mediated alternative polyadenylation is important for plant development and stress responses.

The Nutrient Response Transcriptional Regulome of Arabidopsis

Plants respond actively to changes in their environment. Variations in nutrient availability elicit substantial transcriptional reprogramming, and we aimed to systematically describe these adjustments and identify the regulators responsible. Using gene coexpression analysis based on 13 different nutrient availability anomalies, we defined and analyzed nutrient stress response signatures. We identified known regulators and could predict functions in nutrient responses for transcriptional regulators previously associated with other processes, thus linking development and environmental interaction. Three of the identified transcriptional regulators, PIF4, HY5, and NF-Y, known from their role in light signaling, targeted a substantial part of the network and may participate in remodeling the global Arabidopsis transcriptome in response to variations of nutrient availability. We present gene coexpression and transcriptional regulation networks, which can be used as tools to further explore regulatory events and dependencies even by users with basic informatics skills.

•

Gene coexpression analysis is a powerful tool for elucidating nutrient stress

•

Nutrient stress elicits unique signatures of modular transcriptional response

•

Master transcriptional regulators coordinate plant growth and nutrient utilization

•

Analysis suggests PIF4, HY5, and the NF-Y to be master regulators

Phosphate Starvation Alters Abiotic-Stress-Induced Cytosolic Free Calcium Increases in Roots

Phosphate (Pi) deficiency strongly limits plant growth, and plant roots foraging the soil for nutrients need to adapt to optimize Pi uptake. Ca²⁺ is known to signal in root development and adaptation but has to be tightly controlled, as it is highly toxic to Pi metabolism. Under Pi starvation and the resulting decreased cellular Pi pool, the use of cytosolic free Ca²⁺ ([Ca²⁺]_cyt) as a signal transducer may therefore have to be altered. Employing aequorin-expressing Arabidopsis (Arabidopsis thaliana), we show that Pi starvation, but not nitrogen starvation, strongly dampens the [Ca²⁺]_cyt increases evoked by mechanical, salt, osmotic, and oxidative stress as well as by extracellular nucleotides. The altered root [Ca²⁺]_cyt response to extracellular ATP manifests during seedling development under chronic Pi deprivation but can be reversed by Pi resupply. Employing ratiometric imaging, we delineate that Pi-starved roots have a normal response to extracellular ATP at the apex but show a strongly dampened [Ca²⁺]_cyt response in distal parts of the root tip, correlating with high reactive oxygen species levels induced by Pi starvation. Excluding iron, as well as Pi, rescues this altered [Ca²⁺]_cyt response and restores reactive oxygen species levels to those seen under nutrient-replete conditions. These results indicate that, while Pi availability does not seem to be signaled through [Ca²⁺]_cyt, Pi starvation strongly affects stress-induced [Ca²⁺]_cyt signatures. These data reveal how plants can integrate nutritional and environmental cues, adding another layer of complexity to the use of Ca²⁺ as a signal transducer.

Glycosylphosphatidylinositol-Anchored Proteins in Arabidopsis and One of Their Common Roles in Signaling Transduction

Diverse proteins are found modified with glycosylphosphatidylinositol (GPI) at their carboxyl terminus in eukaryotes, which allows them to associate with membrane lipid bilayers and anchor on the external surface of the plasma membrane. GPI-anchored proteins (GPI-APs) play crucial roles in various processes, and more and more GPI-APs have been identified and studied. In this review, previous genomic and proteomic predictions of GPI-APs in Arabidopsis have been updated, which reveal their high abundance and complexity. From studies of individual GPI-APs in Arabidopsis, certain GPI-APs have been found associated with partner receptor-like kinases (RLKs), targeting RLKs to their subcellular localization and helping to recognize extracellular signaling polypeptide ligands. Interestingly, the association might also be involved in ligand selection. The analyses suggest that GPI-APs are essential and widely involved in signal transduction through association with RLKs.

Tissue-specific transcriptomic profiling of Plantago major provides insights for the involvement of vasculature in phosphate deficiency responses

The vasculature of higher plants is important with transport of both nutrient and information molecules. To understand the correspondence of this tissue in molecular responses under phosphate (Pi) deficiency, Plantago major, a model plant for vasculature biology study, was chosen in our analysis. After RNA-Seq and de novo transcriptome assembly of 24 libraries prepared from the vasculature of P. major, 37,309 unigenes with a mean length of 1571 base pairs were obtained. Upon 24 h of Pi deficiency, 237 genes were shown to be differentially expressed in the vasculature of P. major. Among these genes, only 27 have been previously identified to be specifically expressed in the vasculature tissues in other plant species. Temporal expression of several marker genes associated with Pi deficiency showed that the time period of first 24 h is at the beginning stage of more dynamic expression patterns. In this study, we found several physiological processes, e.g., "phosphate metabolism and remobilization", "sucrose metabolism, loading and synthesis", "plant hormone metabolism and signal transduction", "transcription factors", and "metabolism of other minerals", were mainly involved in early responses to Pi deficiency in the vasculature. A number of vasculature genes with promising roles in Pi deficiency adaptation have been identified and deserve further functional characterization. This study clearly demonstrated that plant vasculature is actively involved in Pi deficiency responses and understanding of this process may help to create plants proficient to offset Pi deficiency.

Molecular Background of Pi Deficiency-Induced Root Hair Growth in Brassica carinata - A Fasciclin-Like Arabinogalactan Protein Is Involved

Formation of longer root hairs under limiting phosphate (P) conditions can increase the inorganic P (Pi) uptake. Here, regulatory candidate genes for Pi deficiency-induced root hair growth were identified by comparison of massive analysis of cDNA ends (MACE) provided expression profiles of two Brassica carinata cultivars (cv.) differing in their root hair response to Pi deficiency: cv. Bale develops longer root hairs under Pi deficiency, but not cv. Bacho. A split-root experiment was conducted for the differentiation between locally and systemically regulated genes. Furthermore, plants were exposed to nitrogen and potassium deficiency to identify P-specific reacting genes. The latter were knocked out by CRISPR/Cas9 and the effect on the root hair length was determined. About 500 genes were differentially expressed under Pi deficiency in cv. Bale, while these genes did not respond to the low P supply in cv. Bacho. Thirty-three candidate genes with a potential regulatory role were selected and the transcriptional regulation of 30 genes was confirmed by quantitative PCR. Only five candidate genes seemed to be either exclusively regulated locally (two) or systemically (three), whereas 25 genes seemed to be involved in both local and systemic signaling pathways. Potassium deficiency affected neither the root hair length nor the expression of the 30 candidate genes. By contrast, both P and nitrogen deficiency increased the root hair length, and both affected the transcript levels in 26 cases. However, four genes reacted specifically to Pi starvation. These genes and, additionally, INORGANIC PHOSPHATE TRANSPORTER 1 (BcPHT1) were targeted by CRISPR/Cas9. However, even if the transcript levels of five of these genes were clearly decreased, FASCICLIN-LIKE ARABINOGALACTAN PROTEIN 1 (BcFLA1) was the only gene whose downregulation reduced the root hair length in transgenic hairy roots under Pi-deficient conditions. To the best of our knowledge, this is the first study describing a fasciclin-like arabinogalactan protein with a predicted role in the Pi deficiency-induced root hair elongation.

Feeding the Walls: How Does Nutrient Availability Regulate Cell Wall Composition?

Nutrients are critical for plants to grow and develop, and nutrient depletion severely affects crop yield. In order to optimize nutrient acquisition, plants adapt their growth and root architecture. Changes in growth are determined by modifications in the cell walls surrounding every plant cell. The plant cell wall, which is largely composed of complex polysaccharides, is essential for plants to attain their shape and to protect cells against the environment. Within the cell wall, cellulose strands form microfibrils that act as a framework for other wall components, including hemicelluloses, pectins, proteins, and, in some cases, callose, lignin, and suberin. Cell wall composition varies, depending on cell and tissue type. It is governed by synthesis, deposition and remodeling of wall components, and determines the physical and structural properties of the cell wall. How nutrient status affects cell wall synthesis and organization, and thus plant growth and morphology, remains poorly understood. In this review, we aim to summarize and synthesize research on the adaptation of root cell walls in response to nutrient availability and the potential role of cell walls in nutrient sensing.

iTRAQ-based analysis of the Arabidopsis proteome reveals insights into the potential mechanisms of anthocyanin accumulation regulation in response to phosphate deficiency

Phosphate (Pi) deficiency significantly limits plant growth in natural and agricultural systems. Accumulation of anthocyanins in shoots is a common response of Arabidopsis thaliana to Pi deficiency. To elucidate the mechanisms underlying Pi deficiency-induced anthocyanin accumulation, we employed a proteomic approach based on isobaric tags for relative and absolute quantification (iTRAQ) to investigate protein expression profiles of Arabidopsis thaliana seedlings subjected to Pi deficiency for 7 days. In total, 5,106 proteins were identified, of which 156 displayed significant changes in abundance upon Pi deficiency. Bioinformatics analysis indicated that flavonoid biosynthesis was the most significantly elevated metabolic process under Pi deficiency. We further examined the potential role of the flavonoid biosynthetic pathway using a dihydroflavonol 4-reductase (DFR) mutant (tt3) and quantitative RT-PCR (qRT-PCR) analysis, and found that the tt3 mutant was deprived of transcriptional up-regulation of three genes related to anthocyanin biosynthesis, modification and transport under Pi deficiency. These results showed that Pi deficiency probably enhances the anthocyanin accumulation by promoting the flavonoid biosynthesis. The exact functions of these proteins remain to be examined. Nevertheless, our study increases the understanding of the mechanisms implicated in the anthocyanin accumulation induced by Pi deficiency and adaptive responses of plants to Pi starvation.

Towards Understanding Plant Calcium Signaling through Calmodulin-Like Proteins: A Biochemical and Structural Perspective

Ca2+ ions play a key role in a wide variety of environmental responses and developmental processes in plants, and several protein families with Ca2+-binding domains have evolved to meet these needs, including calmodulin (CaM) and calmodulin-like proteins (CMLs). These proteins have no catalytic activity, but rather act as sensor relays that regulate downstream targets. While CaM is well-studied, CMLs remain poorly characterized at both the structural and functional levels, even if they are the largest class of Ca2+ sensors in plants. The major structural theme in CMLs consists of EF-hands, and variations in these domains are predicted to significantly contribute to the functional versatility of CMLs. Herein, we focus on recent advances in understanding the features of CMLs from biochemical and structural points of view. The analysis of the metal binding and structural properties of CMLs can provide valuable insight into how such a vast array of CML proteins can coexist, with no apparent functional redundancy, and how these proteins contribute to cellular signaling while maintaining properties that are distinct from CaM and other Ca2+ sensors. An overview of the principal techniques used to study the biochemical properties of these interesting Ca2+ sensors is also presented.

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.

The WUSCHEL-related homeobox 5a (PtoWOX5a) is involved in adventitious root development in poplar

Adventitious rooting is an essential step in vegetative propagation. Currently, the mechanism that regulates adventitious root (AR) development in woody plants is poorly understood. This work demonstrates that Populus tomentosa WUSCHEL-related homeobox 5a (PtoWOX5a) transcription factor is involved in AR development in poplar. PtoWOX5a was specifically expressed in the AR tip and lateral root tip during AR and lateral root regeneration from the stem segment. Phenotypic complementation experiments indicated that the PtoWOX5a can functionally complement AtWOX5 in quiescent center (QC) cells. Overexpression of PtoWOX5a introduces significant developmental phenotypes in roots and leaves, such as increased AR number, decreased AR length, swollen AR tip and lateral root tip, and decreased leaf number and area. The conserved mechanism of D-type cyclins (CYCD) repression mediated by WOX5 was confirmed in poplar. The co-expression network of PtWOX5a was constructed, which provided clues to reveal the molecular mechanism of PtoWOX5a in AR development in poplar. Taken together, our results suggest that the PtoWOX5a is involved in AR development though cooperating with a series of functional genes.

Melonet-DB, a Grand RNA-Seq Gene Expression Atlas in Melon (Cucumis melo L.)

Melon (Cucumis melo L.) is an important Cucurbitaceae crop produced worldwide, exhibiting wide genetic variations and comprising both climacteric and non-climacteric fruit types. The muskmelon cultivar "'Earl's favorite Harukei-3 (Harukei-3)"' known for its sweetness and rich aroma is used for breeding of high-grade muskmelon in Japan. We conducted RNA sequencing (RNA-seq) transcriptome studies in 30 different tissues of the 'Harukei-3' melon. These included root, stems, leaves, flowers, regenerating callus and ovaries, in addition to the flesh and peel sampled at seven stages of fruit development. The expression patterns of 20,752 genes were determined with fragments per kilobase of transcript per million fragments sequenced (FPKM) >1 in at least one tissue. Principal component analysis distinguished 30 melon tissues based on the global gene expression profile and, further, the weighted gene correlation network analysis classified melon genes into 45 distinct coexpression groups. Some coexpression groups exhibited tissue-specific gene expression. Furthermore, we developed and published web application tools designated "'Gene expression map viewer"' and "'Coexpression viewer"' on our website Melonet-DB (http://melonet-db.agbi.tsukuba.ac.jp/) to promote functional genomics research in melon. By using both tools, we analyzed melon homologs of tomato fruit ripening regulators such as E8, RIPENING-INHIBITOR (RIN) and NON-RIPENING (NOR). The "'Coexpression viewer"' clearly distinguished fruit ripening-associated melon RIN/NOR/CNR homologs from those expressed in other tissues. In addition, several other MADS-box, NAM/ATAF/CUC (NAC) and homeobox transcription factor genes were identified as fruit ripening-associated genes. Our tools provide useful information for research not only on melon but also on other fleshy fruit plants.

Nitrate induction of root hair density is mediated by TGA1/TGA4 and CPC transcription factors in Arabidopsis thaliana

Root hairs are specialized cells that are important for nutrient uptake. It is well established that nutrients such as phosphate have a great influence on root hair development in many plant species. Here we investigated the role of nitrate on root hair development at a physiological and molecular level. We showed that nitrate increases root hair density in Arabidopsis thaliana. We found that two different root hair defective mutants have significantly less nitrate than wild-type plants, suggesting that in A. thaliana root hairs have an important role in the capacity to acquire nitrate. Nitrate reductase-null mutants exhibited nitrate-dependent root hair phenotypes comparable with wild-type plants, indicating that nitrate is the signal that leads to increased formation of root hairs. We examined the role of two key regulators of root hair cell fate, CPC and WER, in response to nitrate treatments. Phenotypic analyses of these mutants showed that CPC is essential for nitrate-induced responses of root hair development. Moreover, we showed that NRT1.1 and TGA1/TGA4 are required for pathways that induce root hair development by suppression of longitudinal elongation of trichoblast cells in response to nitrate treatments. Our results prompted a model where nitrate signaling via TGA1/TGA4 directly regulates the CPC root hair cell fate specification gene to increase formation of root hairs in A. thaliana.

Comparative transcriptome analysis of nodules of two Mesorhizobium-chickpea associations with differential symbiotic efficiency under phosphate deficiency

Phosphate (Pi) deficiency is known to be a major limitation for symbiotic nitrogen fixation (SNF), and hence legume crop productivity globally. However, very little information is available on the adaptive mechanisms, particularly in the important legume crop chickpea (Cicer arietinum L.), which enable nodules to respond to low-Pi availability. Thus, to elucidate these mechanisms in chickpea nodules at molecular level, we used an RNA sequencing approach to investigate transcriptomes of the nodules in Mesorhizobium mediterraneum SWRI9-(MmSWRI9)-chickpea and M. ciceri CP-31-(McCP-31)-chickpea associations under Pi-sufficient and Pi-deficient conditions, of which the McCP-31-chickpea association has a better SNF capacity than the MmSWRI9-chickpea association during Pi starvation. Our investigation revealed that more genes showed altered expression patterns in MmSWRI9-induced nodules than in McCP-31-induced nodules (540 vs. 225) under Pi deficiency, suggesting that the Pi-starvation-more-sensitive MmSWRI9-induced nodules required expression change in a larger number of genes to cope with low-Pi stress than the Pi-starvation-less-sensitive McCP-31-induced nodules. The functional classification of differentially expressed genes (DEGs) was examined to gain an understanding of how chickpea nodules respond to Pi starvation, caused by soil Pi deficiency. As a result, more DEGs involved in nodulation, detoxification, nutrient/ion transport, transcriptional factors, key metabolic pathways, Pi remobilization and signalling were found in Pi-starved MmSWRI9-induced nodules than in Pi-starved McCP-31-induced nodules. Our findings have enabled the identification of molecular processes that play important roles in the acclimation of nodules to Pi deficiency, ultimately leading to the development of Pi-efficient chickpea symbiotic associations suitable for Pi-deficient soils.