Plant cell walls constitute complex polysaccharidic/proteinaceous networks whose biosynthesis and dynamics implicate several cell compartments. The synthesis and remodeling of homogalacturonan pectins involve Golgi-localized methylation/acetylation and subsequent cell wall-localized demethylation/deacetylation. So far, TRICHOME BIREFRINGENCE-LIKE (TBL) family members have been described as Golgi-localized acetyltransferases targeting diverse hemicelluloses or pectins. Using seed mucilage secretory cells (MSCs) from Arabidopsis thaliana, we demonstrate the atypical localization of TBL38 restricted to a cell wall microdomain. A tbl38 mutant displays an intriguing homogalacturonan immunological phenotype in this cell wall microdomain and in an MSC surface-enriched abrasion powder. Mass spectrometry oligosaccharide profiling of this fraction reveals an increased homogalacturonan acetylation phenotype. Finally, TBL38 displays pectin acetylesterase activity in vitro. These results indicate that TBL38 is an atypical cell wall-localized TBL that displays a homogalacturonan acetylesterase activity rather than a Golgi-localized acetyltransferase activity as observed in previously studied TBLs. TBL38 function during seed development is discussed.

Plant cell walls are made of complex polysaccharidic/proteinaceous network whose biosynthesis and dynamics implicate several cell compartments and impact plant development. The synthesis and remodeling of homogalacturonan pectins is associated with multiple developmental processes ranging from growth to response to biotic/abiotic stress. It encompasses Golgi-localized methylation and acetylation and subsequent demethylation and deacetylation in the cell wall. In the last decade, our comprehension of plant polysaccharides acetylation has increased significantly thanks to the study of the TRICHOME BIREFRINGENCE-LIKE (TBL) protein family. TBLs are mostly described as Golgi-localized acetyltransferases specifically targeting diverse hemicelluloses or pectins. Various tbl mutants showed altered wall mechanical properties and dynamics. Here, we study TBL38 that is co-expressed with PECTIN METHYLESTERASE INHIBITOR6 (PMEI6) and PEROXIDASE 36 (PRX36) during the development of Arabidopsis seed mucilage secretory cells (MSCs). We demonstrate the atypical TBL38 cell wall localization restricted to the PMEI6/PRX36 MSC cell wall microdomain. A tbl38 mutant displays an intriguing homogalacturonan immunological phenotype in this cell wall microdomain and in a MSC surface-enriched abrasion powder. This fraction was further characterized by mass spectrometry oligosaccharide profiling revealing an increased homogalacturonan acetylation phenotype. Finally, a recombinant TBL38 is shown to display pectin acetylesterase activity in vitro. These results indicate that TBL38 is an atypical cell wall-localized TBL that displays a homogalacturonan acetylesterase activity rather than a Golgi-localized acetyltransferase activity as observed in previously studied TBLs. TBL38 function during seed development is discussed.

Root Hairs (RHs) growth is influenced by endogenous and by external environmental signals that coordinately regulate its final cell size. We have recently determined that RH growth was unexpectedly boosted when Arabidopsis thaliana seedlings are cultivated at low temperatures. It was proposed that RH growth plasticity in response to low temperature was linked to a reduced nutrient availability in the media. Here, we explore the molecular basis of this RH growth response by using a Genome Wide Association Study (GWAS) approach using Arabidopsis thaliana natural accessions. We identify the poorly characterized PEROXIDASE 62 (PRX62) and a related protein PRX69 as key proteins under moderate low temperature stress. Strikingly, a cell wall protein extensin (EXT) reporter reveals the effect of peroxidase activity on EXT cell wall association at 10 °C in the RH apical zone. Collectively, our results indicate that PRX62, and to a lesser extent PRX69, are key apoplastic PRXs that modulate ROS-homeostasis and cell wall EXT-insolubilization linked to RH elongation at low temperature.

Root hair cells are important sensors of soil conditions. They grow towards and absorb water-soluble nutrients. This fast and oscillatory growth is mediated by continuous remodeling of the cell wall. Root hair cell walls contain polysaccharides and hydroxyproline-rich glycoproteins, including extensins (EXTs). Class-III peroxidases (PRXs) are secreted into the apoplastic space and are thought to trigger either cell wall loosening or polymerization of cell wall components, such as Tyr-mediated assembly of EXT networks (EXT-PRXs). The precise role of these EXT-PRXs is unknown. Using genetic, biochemical, and modeling approaches, we identified and characterized three root-hair-specific putative EXT-PRXs, PRX01, PRX44, and PRX73. prx01,44,73 triple mutation and PRX44 and PRX73 overexpression had opposite effects on root hair growth, peroxidase activity, and ROS production, with a clear impact on cell wall thickness. We use an EXT fluorescent reporter with contrasting levels of cell wall insolubilization in prx01,44,73 and PRX44-overexpressing background plants. In this study, we propose that PRX01, PRX44, and PRX73 control EXT-mediated cell wall properties during polar expansion of root hair cells.

Plant cell walls display cellular and subcellular specificities. At the subcellular level, wall regional territories with specific compositions are necessary for mac-roscopic developmental processes. These regional specificities were named dif-ferently throughout the years, and are unified here under the term 'cell-wall microdomains' that define the local composition and organization of wall poly-mers underlying territories of wall loosening and/or softening or stiffening. We re-view the occurrence and developmental role of wall microdomains in different cell types. We primarily focus on the contribution of two categories of wall -remodeling molecular actors: fine-tuning of homogalacturonan (HG; pectin) demethylesterification patterns and two classes of oxidoreductases [class III peroxidases (CIII PRXs) and laccases (LACs)], but we also highlight two different molecular scaffolds recently identified for positioning specific CIII PRXs.

High-throughput data generated by new biotechnologies require specific and adapted statistical treatment in order to be efficiently used in biological studies. In this article, we propose a powerful framework to manage and analyse multi-omics heterogeneous data to carry out an integrative analysis. We have illustrated this using the mixOmics package for R software as it specifically addresses data integration issues. Our work also aims at applying the most recent functionalities of mixOmics to real datasets. Although multi-block integrative methodologies exist, we hope to encourage a more widespread use of such approaches in an operational framework by biologists. We have used natural populations of the model plant Arabidopsis thaliana in this work, but the framework proposed is not limited to this plant and can be deployed whatever the organisms of interest and the biological question may be. Four omics datasets (phenomics, metabolomics, cell wall proteomics and transcriptomics) were collected, analysed and integrated to study the cell wall plasticity of plants exposed to sub-optimal temperature growth conditions. The methodologies presented here start from basic univariate statistics leading to multi-block integration analysis. We have also highlighted the fact that each method, either unsupervised or supervised, is associated with one biological issue. Using this powerful framework enabled us to arrive at novel conclusions on the biological system, which would not have been possible using standard statistical approaches.

Root Hairs (RHs) growth is highly influenced by endogenous as well as by external environmental signals that coordinately regulate its final cell size. RHs actively expand the root surface responsible for nutrient uptake and water absorption. We have recently determined that RH growth was unexpectedly boosted when Arabidopsis thaliana seedlings are cultivated at low temperatures. It was proposed that RH growth plasticity in response to low temperature was linked to a reduced nutrient availability in the media. Here, we explored the molecular basis of this strong RH growth response by using the Genome Wide Association Studies (GWAS) approach on Arabidopsis thaliana natural accessions. We identified the poorly characterized PEROXIDASE 62 (PRX62) as a key protein triggering this conditional growth under a moderate low-temperature stress. In addition, we identified the related protein PRX69 as an important factor in this developmental process. The prx62 prx69 double mutant and the PRX62 and PRX69 over-expressing lines showed contrasting RH phenotypes, peroxidase activities and cyt/apoReactive Oxygen Species (ROS) levels. Strikingly, a cell wall protein extensin (EXT) reporter revealed the effect of peroxidase activity on the EXT cell wall association at 10°C in the RH apical zone. EXT cell wall insolubilization was enhanced at 10°C, which was completely abolished under the PRX inhibitor salicylhydroxamic acid (SHAM) treatment. Finally, we demonstrated that the Root Hair defective 6-like 4 (RSL4) bHLH family transcription factor directly controls the expression of PRX69. Collectively, our results indicate that both PRX62 and PRX69 are key apoplastic PRXs that modulate ROS-homeostasis and cell wall EXT-insolubilization linked to RH elongation at low-temperature.

### Competing Interest Statement The authors have declared no competing interest.

The Class III peroxidases (CIII Prxs) belong to a plant-specific multigene family. Thanks to their double catalytic cycle they can oxidize compounds or release reactive oxygen species (ROS). They are either involved in different cell wall stiffening processes such as lignification and suberization, in cell wall loosening or defense mechanisms. Germination is an important developmental stage requiring specific peroxidase activity. However, little is known about which isoforms are involved. Five CIII Prx encoding genes: AtPrx04, AtPrx16, AtPrx62, AtPrx69, and AtPrx71 were identified from published microarray data mining. Delayed or induced testa and endosperm rupture were observed for the corresponding CIII Prx mutant lines indicating either a gene-specific inducing or repressing role during germination, respectively. Via in situ hybridization AtPrx16, AtPrx62, AtPrx69 and AtPrx71 transcripts were exclusively localized to the micropylar endosperm facing the radicle, and transcriptomic data analysis enabled positioning the five CIII Prxs in a co-expression network enriched in germination, cell wall, cell wall proteins and xyloglucan hits. Evidence were produced showing that the five CIII Prxs were cell wall-targeted proteins and that the micropylar endosperm displayed a complex cell wall domain topochemistry. Finally, we drew a spatio-temporal model highlighting the fine sequential gene expression and the possible involvement of micropylar endosperm cell wall domains to explain the non-redundant cell wall stiffening and loosening functions of the CIII Prxs in a single cell type. We also highlighted the necessity of a peroxidase homeostasis to accurately control the micropylar endosperm cell wall dynamics during Arabidopsis germination events.

The non-animal peroxidases belong to a superfamily of oxidoreductases that reduce hydrogen peroxide and oxidize numerous substrates. Since their initial characterization in 1992, a number of studies have provided an understanding of the origin and evolution of this protein family. Here, we report a comprehensive evolutionary analysis of non-animal peroxidases using integrated in silico and biochemical approaches. Thanks to the availability of numerous genomic sequences from more than 2500 species belonging to 14 kingdoms together with expert and comprehensive annotation of peroxidase sequences that have been centralized in a dedicated database, we have been able to use phylogenetic reconstructions to increase our understanding of the evolutionary processes underlying the diversification of non-animal peroxidases. We analysed the distribution of all non-animal peroxidases in more than 200 eukaryotic organisms in silico. First, we show that the presence or absence of non-animal peroxidases correlates with the presence or absence of certain organelles or with specific biological processes. Examination of almost 2000 organisms determined that ascorbate peroxidases (APxs) and cytochrome c peroxidases (CcPs) are present in those containing chloroplasts and mitochondria, respectively. Plants, which contain both organelles, are an exception and contain only APxs without CcP. Class II peroxidases (CII Prxs) are only found in fungi with wood-decay and plant-degradation abilities. Class III peroxidases (CIII Prxs) are only found in streptophyte algae and land plants, and have been subjected to large family expansion. Biochemical activities of APx, CcP, and CIII Prx assessed using protein extracts from 30 different eukaryotic organisms support the distribution of the sequences resulting from our in silico analysis. The biochemical results confirmed both the presence and classification of the non-animal peroxidase encoding sequences.

AbstractRoot hair cells are important sensors of soil conditions. Expanding several hundred times their original size, root hairs grow towards and absorb water-soluble nutrients. This rapid growth is oscillatory and is mediated by continuous remodelling of the cell wall. Root hair cell walls contain polysaccharides and hydroxyproline-rich glycoproteins including extensins (EXTs).Class-III peroxidases (PRXs) are secreted into the apoplastic space and are thought to trigger either cell wall loosening, mediated by oxygen radical species, or polymerization of cell wall components, including the Tyr-mediated assembly of EXT networks (EXT-PRXs). The precise role of these EXT-PRXs is unknown.Using genetic, biochemical, and modeling approaches, we identified and characterized three root hair-specific putative EXT-PRXs, PRX01, PRX44, and PRX73. The triple mutantprx01,44,73and the PRX44 and PRX73 overexpressors had opposite phenotypes with respect to root hair growth, peroxidase activity and ROS production with a clear impact on cell wall thickness.Modeling and docking calculations suggested that these three putative EXT-PRXs may interact with non-O-glycosylated sections of EXT peptides that reduce the Tyr-to-Tyr intra-chain distances in EXT aggregates and thereby may enhance Tyr crosslinking. These results suggest that these three putative EXT-PRXs control cell wall properties during the polar expansion of root hair cells.

Exogenous low pH stress causes cell death in root cells, limiting root development, and agricultural production. Different lines of evidence suggested a relationship with cell wall (CW) remodeling players. We investigated whether class III peroxidase (CIII Prx) total activity, CIII Prx candidate gene expression, and reactive oxygen species (ROS) could modify CW structure during low pH-induced cell death in Arabidopsis thaliana roots. Wild-type roots displayed a good spatio-temporal correlation between the low pH-induced cell death and total CIII Prx activity in the early elongation (EZs), transition (TZs), and meristematic (MZs) zones. In situ mRNA hybridization showed that AtPrx62 transcripts accumulated only in roots treated at pH 4.6 in the same zones where cell death was induced. Furthermore, roots of the atprx62-1 knockout mutant showed decreased cell mortality under low pH compared to wild-type roots. Among the ROS, there was a drastic decrease in O2●− levels in the MZs of wild-type and atprx62-1 roots upon low pH stress. Together, our data demonstrate that AtPrx62 expression is induced by low pH and that the produced protein could positively regulate cell death. Whether the decrease in O2●− level is related to cell death induced upon low pH treatment remains to be elucidated.

In the global warming context, plant adaptation occurs, but the underlying molecular mechanisms are poorly described. Studying natural variation of the model plant Arabidopsisthaliana adapted to various environments along an altitudinal gradient should contribute to the identification of new traits related to adaptation to contrasted growth conditions. The study was focused on the cell wall (CW) which plays major roles in the response to environmental changes. Rosettes and floral stems of four newly-described populations collected at different altitudinal levels in the Pyrenees Mountains were studied in laboratory conditions at two growth temperatures (22 vs. 15 °C) and compared to the well-described Col ecotype. Multi-omic analyses combining phenomics, metabolomics, CW proteomics, and transcriptomics were carried out to perform an integrative study to understand the mechanisms of plant adaptation to contrasted growth temperature. Different developmental responses of rosettes and floral stems were observed, especially at the CW level. In addition, specific population responses are shown in relation with their environment and their genetics. Candidate genes or proteins playing roles in the CW dynamics were identified and will deserve functional validation. Using a powerful framework of data integration has led to conclusions that could not have been reached using standard statistical approaches.

Natural variations help in identifying genetic mechanisms of morphologically and developmentally complex traits. Mountainous habitats provide an altitudinal gradient where one species encounters different abiotic conditions. We report the study of 341 individuals of Arabidopsis thaliana derived from 30 natural populations not belonging to the 1001 genomes, collected at increasing altitudes, between 200 and 1800 m in the Pyrenees. Class III peroxidases and ribosomal RNA sequences were used as markers to determine the putative genetic relationships among these populations along their altitudinal gradient. Using Bayesian-based statistics and phylogenetic analyses, these Pyrenean populations appear with significant divergence from the other regional accessions from 1001 genome (i.e., from north Spain or south France). Individuals of these populations exhibited varying phenotypic changes, when grown at sub-optimal temperature (22 vs. 15°C). These phenotypic variations under controlled conditions reflected intraspecific morphological variations. This study could bring new information regarding the west European population structure of A. thaliana and its phenotypic variations at different temperatures. The integrative analysis combining genetic, phenotypic variation and environmental datasets is used to analyze the acclimation of population in response to temperature changes. Regarding their geographical proximity and environmental diversity, these populations represent a tool of choice for studying plant response to temperature variation.HIGHLIGHTS-Studying the natural diversity of A. thaliana in the Pyrenees mountains helps to understand European population structure and to evaluate the phenotypic trait variation in response to climate change.

(m)RNA spatiotemporal pattern of distribution is of key importance to decipher gene function. In this post-genomic era, numerous transcriptomic studies are made publicly available, sometimes reaching a tissular resolution and even more rarely the cellular level. This "one tissue-numerous genes" information can be completed by the reverse "one gene-numerous tissues" picture through traditional RNA in situ hybridization (ISH). Here, we present a method including (1) principles of transcriptomic data mining to be performed prior and following ISH and (2) a detailed step-by-step medium-throughput ISH protocol performed on serial sections from tissue microarrays. In a recent work, we implemented this method for 39 selected genes studied by medium-throughput ISH complementing an existing tissue-specific transcriptomic dataset focused on the model plant Arabidopsis seed development kinetics (Francoz et al., Scientific Reports 6:24644, 2016). This full integration of ISH and transcriptomics demonstrated the complementarity of both techniques in terms of tissue/cell specificity, signal sensitivity, gene specificity, and spatiotemporal resolution.

This article provides experimental data describing the RNA and the cell wall protein profiles of rosettes and flower stems of five Arabidopsis thaliana ecotypes. Four newly-described Pyrenees ecotypes [1] are analyzed in addition to the well-described and sequenced Columbia (Col) ecotype of A. thaliana. All five ecotypes have been grown at two different temperatures, 22 °C and 15 °C. We provide transcriptomics and cell wall proteomics data regarding (i) rosettes at the bolting stage, and (ii) floral stems at the first flower stage. These data are a valuable resource to study the adaptation of A. thaliana ecotypes to sub-optimal temperature growth conditions.

Plant cell walls are made of polysaccharidic-proteinaceous complex matrices. Molecular interactions governing their organization remain understudied. We take advantage of the highly dynamic cell walls of Arabidopsis seed mucilage secretory cells to propose a hierarchical multi-molecular interaction model within a cell wall domain. We show that the PECTINMETHYLESTERASE INHIBITOR6 activity creates a partially demethylesterified pectin pattern acting as a platform allowing positioning of PEROXIDASE36 in a remote primary cell wall domain during early development. This allows triggering the loosening of this domain during later development, in turn leading to proper physiological function upon mature seed imbibition and germination. We anticipate that this pioneer example of molecular scaffold within a cell wall domain is more widespread through other combinations of the individual molecular players all belonging to large multigenic families. These results highlight the role of cell wall polysaccharide-protein interactions in the organization of cell wall domains.

With the global temperature change, plant adaptations are predicted, but little is known about the molecular mechanisms underlying them. Arabidopsis thaliana is a model plant adapted to various environmental conditions, in particular able to develop along an altitudinal gradient. Two ecotypes, Columbia (Col) growing at low altitude, and Shahdara (Sha) growing at 3400m, have been studied at optimal and sub-optimal growth temperature (22°C vs 15°C). Macro- and micro-phenotyping, cell wall monosaccharides analyses, cell wall proteomics, and transcriptomics have been performed in order to accomplish an integrative analysis. The analysis has been focused on cell walls (CWs) which are assumed to play roles in response to environmental changes. At 15°C, both ecotypes presented characteristic morphological traits of low temperature growth acclimation such as reduced rosette diameter, increased number of leaves, modifications of their CW composition and cuticle reinforcement. Altogether, the integrative analysis has allowed identifying several candidate genes/proteins possibly involved in the cell wall modifications observed during the temperature acclimation response.

Root hair polar growth is endogenously controlled by auxin and sustained by oscillating levels of reactive oxygen species (ROS). These cells extend several hundred-fold their original size toward signals important for plant survival. Although their final cell size is of fundamental importance, the molecular mechanisms that control it remain largely unknown. Here we show that ROS production is controlled by the transcription factor RSL4, which in turn is transcriptionally regulated by auxin through several auxin response factors (ARFs). In this manner, auxin controls ROS-mediated polar growth by activating RSL4, which then up-regulates the expression of genes encoding NADPH oxidases (also known as RESPIRATORY BURST OXIDASE HOMOLOG proteins) and class III peroxidases, which catalyze ROS production. Chemical or genetic interference with ROS balance or peroxidase activity affects root hair final cell size. Overall, our findings establish a molecular link between auxin and ROS-mediated polar root hair growth.