Evolutionary Analysis of Six Gene Families Part of the Reactive Oxygen Species (ROS) Gene Network in Three Brassicaceae Species.

Climate change is expected to intensify the occurrence of abiotic stress in plants, such as hypoxia and salt stresses, leading to the production of reactive oxygen species (ROS), which need to be effectively managed by various oxido-reductases encoded by the so-called ROS gene network. Here, we studied six oxido-reductases families in three Brassicaceae species, Arabidopsis thaliana as well as Nasturtium officinale and Eutrema salsugineum, which are adapted to hypoxia and salt stress, respectively. Using available and new genomic data, we performed a phylogenomic analysis and compared RNA-seq data to study genomic and transcriptomic adaptations. This comprehensive approach allowed for the gaining of insights into the impact of the adaptation to saline or hypoxia conditions on genome organization (gene gains and losses) and transcriptional regulation. Notably, the comparison of the N. officinale and E. salsugineum genomes to that of A. thaliana highlighted changes in the distribution of ohnologs and homologs, particularly affecting class III peroxidase genes (CIII Prxs). These changes were specific to each gene, to gene families subjected to duplication events and to each species, suggesting distinct evolutionary responses. The analysis of transcriptomic data has allowed for the identification of genes related to stress responses in A. thaliana, and, conversely, to adaptation in N. officinale and E. salsugineum.

A powerful framework for an integrative study with heterogeneous omics data: from univariate statistics to multi-block analysis.

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.

Class III Peroxidases in Response to Multiple Abiotic Stresses in Arabidopsis thaliana Pyrenean Populations.

Class III peroxidases constitute a plant-specific multigene family, where 73 genes have been identified in Arabidopsis thaliana. These genes are members of the reactive oxygen species (ROS) regulatory network in the whole plant, but more importantly, at the root level. In response to abiotic stresses such as cold, heat, and salinity, their expression is significantly modified. To learn more about their transcriptional regulation, an integrative phenotypic, genomic, and transcriptomic study was executed on the roots of A. thaliana Pyrenean populations. Initially, the root phenotyping highlighted 3 Pyrenean populations to be tolerant to cold (Eaux), heat (Herr), and salt (Grip) stresses. Then, the RNA-seq analyses on these three populations, in addition to Col-0, displayed variations in CIII Prxs expression under stressful treatments and between different genotypes. Consequently, several CIII Prxs were particularly upregulated in the tolerant populations, suggesting novel and specific roles of these genes in plant tolerance against abiotic stresses.

Class III Peroxidases PRX01, PRX44, and PRX73 Control Root Hair Growth in Arabidopsis thaliana.

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.

Primary transcripts of microRNAs encode regulatory peptides.

MicroRNAs (miRNAs) are small regulatory RNA molecules that inhibit the expression of specific target genes by binding to and cleaving their messenger RNAs or otherwise inhibiting their translation into proteins. miRNAs are transcribed as much larger primary transcripts (pri-miRNAs), the function of which is not fully understood. Here we show that plant pri-miRNAs contain short open reading frame sequences that encode regulatory peptides. The pri-miR171b of Medicago truncatula and the pri-miR165a of Arabidopsis thaliana produce peptides, which we term miPEP171b and miPEP165a, respectively, that enhance the accumulation of their corresponding mature miRNAs, resulting in downregulation of target genes involved in root development. The mechanism of miRNA-encoded peptide (miPEP) action involves increasing transcription of the pri-miRNA. Five other pri-miRNAs of A. thaliana and M. truncatula encode active miPEPs, suggesting that miPEPs are widespread throughout the plant kingdom. Synthetic miPEP171b and miPEP165a peptides applied to plants specifically trigger the accumulation of miR171b and miR165a, leading to reduction of lateral root development and stimulation of main root growth, respectively, suggesting that miPEPs might have agronomical applications.

Effects of Dielectric Barrier Ambient Air Plasma on Two Brassicaceae Seeds: Arabidopsis thaliana and Camelina sativa.

We investigated low-temperature plasma effects on two Brassicaceae seeds (A. thaliana and C. sativa) using dielectric barrier discharge in air. Comparisons of plasma treatments on seeds showed distinct responses on germination rate and speed. Optimal treatment time giving optimal germination is 15 min for A. thaliana with 85% increase compared to control after 48 h of germination and 1 min for C. sativa with 75% increase compared to control after 32 h of germination. Such germination increases are associated with morphological changes shown by SEM of seed surface. For better understanding at the biochemical level, seed surfaces were analyzed using gas chromatography-mass spectrometry which underlined changes of lipidic composition. For both treated seeds, there is a decrease of saturated (palmitic and stearic) fatty acids while treated C. sativa showed a decrease of unsaturated (oleic and linoleic) acids and treated A. thaliana an increase of unsaturated ones. Such lipid changes, specifically a decrease of hydrophobic saturated fatty acids, are coherent with the other analyses (SEM, water uptake and contact angle). Moreover, an increase in A. thaliana of unsaturated acids (very reactive) probably neutralizes plasma RONS effects thus needing longer plasma exposure time (15 min) to reach optimal germination. For C. sativa, 1 min is enough because unsaturated linoleic acid becomes lower in treated C. sativa (1.2 × 107) compared to treated A. thaliana (3.7 × 107).

Seed mucilage evolution: Diverse molecular mechanisms generate versatile ecological functions for particular environments.

Plant myxodiasporous species have the ability to release a polysaccharidic mucilage upon imbibition of the seed (myxospermy) or the fruit (myxocarpy). This is a widespread capacity in angiosperms providing multiple ecological functions including higher germination efficiency under environmental stresses. It is unclear whether myxodiaspory has one or multiple evolutionary origins and why it was supposedly lost in several species. Here, we summarize recent advances on three main aspects of myxodiaspory. (a) It represents a combination of highly diverse traits at different levels of observation, ranging from the dual tissular origin of mucilage secretory cells to diverse mucilage polysaccharidic composition and ultrastructural organization. (b) An asymmetrical selection pressure is exerted on myxospermy-related genes that were first identified in Arabidopsis thaliana. The A. thaliana and the flax intra-species mucilage variants show that myxospermy is a fast-evolving trait due to high polymorphism in a few genes directly acting on mucilage establishment. In A. thaliana, these actors are downstream of a master regulatory complex and an original phylogenetic overview provided here illustrates that this complex has sequentially evolved after the common ancestor of seed plants and was fully established in the common ancestor of the rosid clade. (c) Newly identified myxodiaspory ecological functions indicate new perspectives such as soil microorganism control and plant establishment support.

Global analysis of non-animal peroxidases provides insights into the evolution of this gene family in the green lineage.

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.

The genome of Eucalyptus grandis.

Eucalypts are the world's most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

Molecular link between auxin and ROS-mediated polar growth.

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.

Plant Cell Wall Proteins and Development.

Plant cell walls surround cells and provide both external protection and a means of cell-to-cell communication [...].

The Class III Peroxidase Encoding Gene AtPrx62 Positively and Spatiotemporally Regulates the Low pH-Induced Cell Death in Arabidopsis thaliana Roots.

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.

The Cell Wall PAC (Proline-Rich, Arabinogalactan Proteins, Conserved Cysteines) Domain-Proteins Are Conserved in the Green Lineage.

Plant cell wall proteins play major roles during plant development and in response to environmental cues. A bioinformatic search for functional domains has allowed identifying the PAC domain (Proline-rich, Arabinogalactan proteins, conserved Cysteines) in several proteins (PDPs) identified in cell wall proteomes. This domain is assumed to interact with pectic polysaccharides and O-glycans and to contribute to non-covalent molecular scaffolds facilitating the remodeling of polysaccharidic networks during rapid cell expansion. In this work, the characteristics of the PAC domain are described in detail, including six conserved Cys residues, their spacing, and the predicted secondary structures. Modeling has been performed based on the crystal structure of a Plantago lanceolata PAC domain. The presence of ß-sheets is assumed to ensure the correct folding of the PAC domain as a ß-barrel with loop regions. We show that PDPs are present in early divergent organisms from the green lineage and in all land plants. PAC domains are associated with other types of domains: Histidine-rich, extensin, Proline-rich, or yet uncharacterized. The earliest divergent organisms having PDPs are Bryophytes. Like the complexity of the cell walls, the number and complexity of PDPs steadily increase during the evolution of the green lineage. The association of PAC domains with other domains suggests a neo-functionalization and different types of interactions with cell wall polymers.

Algal ancestor of land plants was preadapted for symbiosis.

Colonization of land by plants was a major transition on Earth, but the developmental and genetic innovations required for this transition remain unknown. Physiological studies and the fossil record strongly suggest that the ability of the first land plants to form symbiotic associations with beneficial fungi was one of these critical innovations. In angiosperms, genes required for the perception and transduction of diffusible fungal signals for root colonization and for nutrient exchange have been characterized. However, the origin of these genes and their potential correlation with land colonization remain elusive. A comprehensive phylogenetic analysis of 259 transcriptomes and 10 green algal and basal land plant genomes, coupled with the characterization of the evolutionary path leading to the appearance of a key regulator, a calcium- and calmodulin-dependent protein kinase, showed that the symbiotic signaling pathway predated the first land plants. In contrast, downstream genes required for root colonization and their specific expression pattern probably appeared subsequent to the colonization of land. We conclude that the most recent common ancestor of extant land plants and green algae was preadapted for symbiotic associations. Subsequent improvement of this precursor stage in early land plants through rounds of gene duplication led to the acquisition of additional pathways and the ability to form a fully functional arbuscular mycorrhizal symbiosis.

Cell wall proteome analysis of Arabidopsis thaliana mature stems.

Plant stems carry flowers necessary for species propagation and need to be adapted to mechanical disturbance and environmental factors. The stem cell walls are different from other organs and can modify their rigidity or viscoelastic properties for the integrity and the robustness required to withstand mechanical impacts and environmental stresses. Plant cell wall is composed of complex polysaccharide networks also containing cell wall proteins (CWPs) crucial to perceive and limit the environmental effects. The CWPs are fundamental players in cell wall remodeling processes, and today, only 86 have been identified from the mature stems of the model plant Arabidopsis thaliana. With a destructive method, this study has enlarged its coverage to 302 CWPs. This new proteome is mainly composed of 27.5% proteins acting on polysaccharides, 16% proteases, 11.6% oxido-reductases, 11% possibly related to lipid metabolism and 11% of proteins with interacting domains with proteins or polysaccharides. Compared to stem cell wall proteomes already available (Brachypodium distachyon, Sacharum officinarum, Linum usitatissimum, Medicago sativa), that of A. thaliana stems has a higher proportion of proteins acting on polysaccharides and of proteases, but a lower proportion of oxido-reductases.

Automatic multigenic family annotation: risks and solutions.

A major challenge facing bioinformatics today is the efficient annotation of the exponential flow of genomic data. This has led to an increasing dependence on automatic annotation procedures, despite the relatively high error rates of these programs, particularly for multigenic families. We discuss here the errors and biases introduced by automatic genome annotations, focusing on issues with structural annotations of gene families, and suggest ways to overcome these limitations.

The class III peroxidase PRX17 is a direct target of the MADS-box transcription factor AGAMOUS-LIKE15 (AGL15) and participates in lignified tissue formation.

Several physiological functions have been attributed to class III peroxidases (PRXs) in plants, but the in planta role of most members of this family still remains undetermined. Here, we report the first functional characterization of PRX17 (At2g22420), one of the 73 members of this family in Arabidopsis thaliana. Localization of PRX17 was examined by transient expression in Nicotiana benthamiana. Loss- and gain-of-function mutants in A. thaliana were studied. Regulation at the gene and protein levels was analyzed using ß-glucuronidase (GUS) activity, quantitative reverse transcriptase (qRT)-PCR, zymography, and chromatin immunoprecipitation. Phenotypes were characterized including lignin and xylan contents. PRX17 was expressed in various tissues, including vascular tissues, and PRX17 was localized to the cell wall. In prx17, the lignin content was reduced in the stem and siliques and bolting was delayed, while the opposite phenotype was observed in 35S:PRX17 plants, together with a significant increase of lignin and xylan immunofluorescence signal. Finally, we demonstrated that the transcription factor AGAMOUS-LIKE15 (AGL15) binds to the PRX17 promoter and regulates PRX17 expression level. This converging set of structural, transcriptomic and physiological data suggests that PRX17, under the control of AGL15, contributes to developmental programs by playing an essential role in regulating age-dependent lignified tissue formation, including changes in cell wall properties.

Arabidopsis thaliana root cell wall proteomics: Increasing the proteome coverage using a combinatorial peptide ligand library and description of unexpected Hyp in peroxidase amino acid sequences.

Plant cell walls (CWs) contain a large proportion of polysaccharides (90-95% of CW mass) and proteins (5-10%) that play major roles in CW plasticity during development and in response to environmental cues. Here, we present CW proteomics data of Arabidopsis thaliana roots. Plants were cultivated in hydroponic conditions. CW protein (CWP) extracts were prepared and analyzed in two different ways in order to enlarge the coverage of the root CW proteome: proteins were analyzed either directly or following an affinity chromatography on a combinatorial peptide ligand library (CPLL) to reduce the concentration dynamic range. Proteins were identified by LC-MS/MS and bioinformatics. Altogether, 424 proteins having predicted signal peptides have been identified (CWPs). CPLL permitted to identify low-abundant CWPs never described before, thus enlarging the coverage of the root CW proteome. The number of oxidoreductases is particularly high and includes a large collection of class III peroxidases (CIII Prxs; 38 out of the 73 A. thaliana CIII Prxs). For the first time, hydroxyproline residues were localized at conserved positions in CIII Prx amino acid sequences.

An enlarged cell wall proteome of Arabidopsis thaliana rosettes.

Plant cells are surrounded by cell walls playing many roles during development and in response to environmental constraints. Cell walls are mainly composed of polysaccharides (cellulose, hemicelluloses and pectins), but they also contain proteins which are critical players in cell wall remodeling processes. Today, the cell wall proteome of Arabidopsis thaliana, a major dicot model plant, comprises more than 700 proteins predicted to be secreted (cell wall proteins-CWPs) identified in different organs or in cell suspension cultures. However, the cell wall proteome of rosettes is poorly represented with only 148 CWPs identified after extraction by vacuum infiltration. This new study allows enlarging its coverage. A destructive method starting with the purification of cell walls has been performed and two experiments have been compared. They differ by the presence/absence of protein separation by a short 1D-electrophoresis run prior to tryptic digestion and different gradient programs for peptide separation before mass spectrometry analysis. Altogether, the rosette cell wall proteome has been significantly enlarged to 361 CWPs, among which 213 newly identified in rosettes and 57 newly described. The identified CWPs fall in four major functional classes: 26.1% proteins acting on polysaccharides, 11.1% oxido-reductases, 14.7% proteases and 11.7% proteins possibly related to lipid metabolism.

The Arabidopsis Class III Peroxidase AtPRX71 Negatively Regulates Growth under Physiological Conditions and in Response to Cell Wall Damage.

The structure of the cell wall has a major impact on plant growth and development, and alteration of cell wall structural components is often detrimental to biomass production. However, the molecular mechanisms responsible for these negative effects are largely unknown. Arabidopsis (Arabidopsis thaliana) plants with altered pectin composition because of either the expression of the Aspergillus niger polygalacturonase II (AnPGII; 35S:AnPGII plants) or a mutation in the QUASIMODO2 (QUA2) gene that encodes a putative pectin methyltransferase (qua2-1 plants), display severe growth defects. Here, we show that expression of Arabidopsis PEROXIDASE71 (AtPRX71), encoding a class III peroxidase, strongly increases in 35S:AnPGII and qua2-1 plants as well as in response to treatments with the cellulose synthase inhibitor isoxaben, which also impairs cell wall integrity. Analysis of atprx71 loss-of-function mutants and plants overexpressing AtPRX71 indicates that this gene negatively influences Arabidopsis growth at different stages of development, likely limiting cell expansion. The atprx71-1 mutation partially suppresses the dwarf phenotype of qua2-1, suggesting that AtPRX71 contributes to the growth defects observed in plants undergoing cell wall damage. Furthermore, AtPRX71 seems to promote the production of reactive oxygen species in qua2-1 plants as well as plants treated with isoxaben. We propose that AtPRX71 contributes to strengthen cell walls, therefore restricting cell expansion, during normal growth and in response to cell wall damage.