Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML.

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

Advanced Cataloging of Lysine-63 Polyubiquitin Networks by Genomic, Interactome, and Sensor-Based Proteomic Analyses.

The lack of resolution when studying the many different ubiquitin chain types found in eukaryotic cells has been a major hurdle to our understanding of their specific roles. We currently have very little insight into the cellular and physiological functions of Lys-63 (K63)-linked ubiquitin chains, although they are the second most abundant forms of ubiquitin in plant cells. To overcome this problem, we developed several large-scale approaches to characterize (1) the E2-E3 ubiquitination machinery driving K63-linked ubiquitin chain formation and (2) K63 polyubiquitination targets to provide a comprehensive picture of K63 polyubiquitin networks in Arabidopsis (Arabidopsis thaliana). Our work identified the ubiquitin-conjugating enzymes (E2s) UBC35/36 as the major drivers of K63 polyubiquitin chain formation and highlights the major role of these proteins in plant growth and development. Interactome approaches allowed us to identify many proteins that interact with the K63 polyubiquitination-dedicated E2s UBC35/36 and their cognate E2 variants, including more than a dozen E3 ligases and their putative targets. In parallel, we improved the in vivo detection of proteins decorated with K63-linked ubiquitin chains by sensor-based proteomics, yielding important insights into the roles of K63 polyubiquitination in plant cells. This work strongly increases our understanding of K63 polyubiquitination networks and functions in plants.

Plant Cell Wall Proteomes: The Core of Conserved Protein Families and the Case of Non-Canonical Proteins.

Plant cell wall proteins (CWPs) play critical roles during plant development and in response to stresses. Proteomics has revealed their great diversity. With nearly 1000 identified CWPs, the Arabidopsis thaliana cell wall proteome is the best described to date and it covers the main plant organs and cell suspension cultures. Other monocot and dicot plants have been studied as well as bryophytes, such as Physcomitrella patens and Marchantia polymorpha. Although these proteomes were obtained using various flowcharts, they can be searched for the presence of members of a given protein family. Thereby, a core cell wall proteome which does not pretend to be exhaustive, yet could be defined. It comprises: (i) glycoside hydrolases and pectin methyl esterases, (ii) class III peroxidases, (iii) Asp, Ser and Cys proteases, (iv) non-specific lipid transfer proteins, (v) fasciclin arabinogalactan proteins, (vi) purple acid phosphatases and (vii) thaumatins. All the conserved CWP families could represent a set of house-keeping CWPs critical for either the maintenance of the basic cell wall functions, allowing immediate response to environmental stresses or both. Besides, the presence of non-canonical proteins devoid of a predicted signal peptide in cell wall proteomes is discussed in relation to the possible existence of alternative secretion pathways.

A Tandem Mass Tags (TMTs) labeling approach highlights differences between the shoot proteome of two Arabidopsis thaliana ecotypes, Col-0 and Ws.

Arabidopsis has become a powerful model to study morphogenesis, plant growth, development but also plant response to environmental conditions. Over 1000 Arabidopsis genomes are available and show natural genetic variations. Among them, the main reference accessions Wassilewskija (Ws) and Columbia (Col-0), originally growing at contrasted altitudes and temperatures, are widely studied, but data contributing to their molecular phenotyping are still scarce. A global quantitative proteomics approach using isobaric stable isotope labeling (Tandem Mass Tags, TMT) was performed on Ws and Col-0. Plants have been hydroponically grown at 16 h/8 h (light/dark cycle) at 23°C day/19°C night for three weeks. A TMT labeling of the proteins extracted from their shoots has been performed and showed a differential pattern of protein abundance between them. These results have allowed identifying several proteins families possibly involved in the differential responses observed for Ws and Col-0 during plant development and upon environmental changes. In particular, Ws and Col-0 mainly differ in photosynthesis, cell wall-related proteins, plant defense/stress, ROS scavenging enzymes/redox homeostasis and DNA/RNA binding/transcription/translation/protein folding.

The Same against Many: AtCML8, a Ca2+ Sensor Acting as a Positive Regulator of Defense Responses against Several Plant Pathogens.

Calcium signals are crucial for the activation and coordination of signaling cascades leading to the establishment of plant defense mechanisms. Here, we studied the contribution of CML8, an Arabidopsis calmodulin-like protein in response to Ralstonia solanacearum and to pathogens with different lifestyles, such as Xanthomonas campestris pv. campestris and Phytophtora capsici. We used pathogenic infection assays, gene expression, RNA-seq approaches, and comparative analysis of public data on CML8 knockdown and overexpressing Arabidopsis lines to demonstrate that CML8 contributes to defense mechanisms against pathogenic bacteria and oomycetes. CML8 gene expression is finely regulated at the root level and manipulated during infection with Ralstonia, and CML8 overexpression confers better plant tolerance. To understand the processes controlled by CML8, genes differentially expressed at the root level in the first hours of infection have been identified. Overexpression of CML8 also confers better tolerance against Xanthomonas and Phytophtora, and most of the genes differentially expressed in response to Ralstonia are differentially expressed in these different pathosystems. Collectively, CML8 acts as a positive regulator against Ralstonia solanaceraum and against other vascular or root pathogens, suggesting that CML8 is a multifunctional protein that regulates common downstream processes involved in the defense response of plants to several pathogens.

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.

Identification of new signalling peptides through a genome-wide survey of 250 fungal secretomes.

BACKGROUND: Many small peptides regulate eukaryotic cell biology. In fungi, some of these peptides are produced after KEX2 protease activity on proteins displaying repetitions of identical or nearly identical motifs. Following this endoprotease activity, peptides are released in the extracellular space. This type of protein maturation is involved in the production of the α-type sexual pheromone in Ascomycota. In other cases, this processing allows the production of secreted peptides regulating fungal cell wall structure or acting as mycotoxins. In this work, we report for the first time a genome-wide search of KEX2-processed repeat proteins that we call KEPs. We screened the secreted proteins of 250 fungal species to compare their KEP repertoires with regard to their lifestyle, morphology or lineage. RESULTS: Our analysis points out that nearly all fungi display putative KEPs, suggesting an ancestral origin common to all opisthokonts. As expected, our pipeline identifies mycotoxins but also α-type sexual pheromones in Ascomycota that have not been explored so far, and unravels KEP-derived secreted peptides of unknown functions. Some species display an expansion of this class of proteins. Interestingly, we identified conserved KEPs in pathogenic fungi, suggesting a role in virulence. We also identified KEPs in Basidiomycota with striking similarities to Ascomycota α-type sexual pheromones, suggesting they may also play alternative roles in unknown signalling processes. CONCLUSIONS: We identified putative, new, unexpected secreted peptides that fall into different functional categories: mycotoxins, hormones, sexual pheromones, or effectors that promote colonization during host-microbe interactions. This wide survey will open new avenues in the field of small-secreted peptides in fungi that are critical regulators of their intimate biology and modulators of their interaction with the environment.

A systems biology view of wood formation in Eucalyptus grandis trees submitted to different potassium and water regimes.

Wood production in fast-growing Eucalyptus grandis trees is highly dependent on both potassium (K) fertilization and water availability but the molecular processes underlying wood formation in response to the combined effects of these two limiting factors remain unknown. E. grandis trees were submitted to four combinations of K-fertilization and water supply. Weighted gene co-expression network analysis and MixOmics-based co-regulation networks were used to integrate xylem transcriptome, metabolome and complex wood traits. Functional characterization of a candidate gene was performed in transgenic E. grandis hairy roots. This integrated network-based approach enabled us to identify meaningful biological processes and regulators impacted by K-fertilization and/or water limitation. It revealed that modules of co-regulated genes and metabolites strongly correlated to wood complex traits are in the heart of a complex trade-off between biomass production and stress responses. Nested in these modules, potential new cell-wall regulators were identified, as further confirmed by the functional characterization of EgMYB137. These findings provide new insights into the regulatory mechanisms of wood formation under stressful conditions, pointing out both known and new regulators co-opted by K-fertilization and/or water limitation that may potentially promote adaptive wood traits.

Comparative genomics of Rhizophagus irregularis, R. cerebriforme, R. diaphanus and Gigaspora rosea highlights specific genetic features in Glomeromycotina.

Glomeromycotina is a lineage of early diverging fungi that establish arbuscular mycorrhizal (AM) symbiosis with land plants. Despite their major ecological role, the genetic basis of their obligate mutualism remains largely unknown, hindering our understanding of their evolution and biology. We compared the genomes of Glomerales (Rhizophagus irregularis, Rhizophagus diaphanus, Rhizophagus cerebriforme) and Diversisporales (Gigaspora rosea) species, together with those of saprotrophic Mucoromycota, to identify gene families and processes associated with these lineages and to understand the molecular underpinning of their symbiotic lifestyle. Genomic features in Glomeromycotina appear to be very similar with a very high content in transposons and protein-coding genes, extensive duplications of protein kinase genes, and loss of genes coding for lignocellulose degradation, thiamin biosynthesis and cytosolic fatty acid synthase. Most symbiosis-related genes in R. irregularis and G. rosea are specific to Glomeromycotina. We also confirmed that the present species have a homokaryotic genome organisation. The high interspecific diversity of Glomeromycotina gene repertoires, affecting all known protein domains, as well as symbiosis-related orphan genes, may explain the known adaptation of Glomeromycotina to a wide range of environmental settings. Our findings contribute to an increasingly detailed portrait of genomic features defining the biology of AM fungi.

Genomics analysis of Aphanomyces spp. identifies a new class of oomycete effector associated with host adaptation.

BACKGROUND: Oomycetes are a group of filamentous eukaryotic microorganisms that have colonized all terrestrial and oceanic ecosystems, and they include prominent plant pathogens. The Aphanomyces genus is unique in its ability to infect both plant and animal species, and as such exemplifies oomycete versatility in adapting to different hosts and environments. Dissecting the underpinnings of oomycete diversity provides insights into their specificity and pathogenic mechanisms. RESULTS: By carrying out genomic analyses of the plant pathogen A. euteiches and the crustacean pathogen A. astaci, we show that host specialization is correlated with specialized secretomes that are adapted to the deconstruction of the plant cell wall in A. euteiches and protein degradation in A. astaci. The A. euteiches genome is characterized by a large repertoire of small secreted protein (SSP)-encoding genes that are highly induced during plant infection, and are not detected in other oomycetes. Functional analysis revealed an SSP from A. euteiches containing a predicted nuclear-localization signal which shuttles to the plant nucleus and increases plant susceptibility to infection. CONCLUSION: Collectively, our results show that Aphanomyces host adaptation is associated with evolution of specialized secretomes and identify SSPs as a new class of putative oomycete effectors.

Plant extracellular vesicles are incorporated by a fungal pathogen and inhibit its growth.

Extracellular vesicles (EV) are membrane particles released by cells into their environment and are considered to be key players in intercellular communication. EV are produced by all domains of life but limited knowledge about EV in plants is available, although their implication in plant defense has been suggested. We have characterized sunflower EV and tested whether they could interact with fungal cells. EV were isolated from extracellular fluids of seedlings and characterized by transmission electron microscopy and proteomic analysis. These nanovesicles appeared to be enriched in cell wall remodeling enzymes and defense proteins. Membrane-labeled EV were prepared and their uptake by the phytopathogenic fungus Sclerotinia sclerotiorum was verified. Functional tests further evaluated the ability of EV to affect fungal growth. Spores treated with plant EV showed growth inhibition, morphological changes, and cell death. Conclusive evidence on the existence of plant EV is presented and we demonstrate their ability to interact with and kill fungal cells. Our results introduce the concept of cell-to-cell communication through EV in plants.

Special trends in CBF and DREB2 groups in Eucalyptus gunnii vs Eucalyptus grandis suggest that CBF are master players in the trade-off between growth and stress resistance.

Annotation of the Eucalyptus grandis genome showed a large amplification of the dehydration-responsive element binding 1/C-repeat binding factor (DREB1/CBF) group without recent DREB2 gene duplication compared with other plant species. The present annotation of the CBF and DREB2 genes from a draft of the Eucalyptus gunnii genome sequence reveals at least one additional CBF copy in the E. gunnii genome compared with E. grandis, suggesting that this group is still evolving, unlike the DREB2 group. This study aims to investigate the redundancy/neo- or sub-functionalization of the duplicates and the relative involvement of the two groups in abiotic stress responses in both E. grandis and E. gunnii (lower growth but higher cold resistance). A comprehensive transcriptional analysis using high-throughput quantitative real-time polymerase chain reaction (qRT-PCR) was performed on leaves, stems and roots from the two Eucalyptus species after cold, heat or drought treatment. A large CBF cluster accounted for most of the cold response in all the organs, whereas heat and drought responses mainly involved a small CBF cluster and the DREB2 genes. In addition, CBF putative target genes, known to be involved in plant tolerance and development, were found to be cold-regulated. The higher transcript amounts of both the CBF and target genes in the cold tolerant E. gunnii contrasted with the higher CBF induction rates in the fast growing E. grandis. Altogether, the present results, in agreement with previous data about Eucalyptus transgenic lines over-expressing CBF, suggest that these factors, which promote both stress protection and growth limitation, participate in the trade-off between growth and resistance in this woody species.

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.

Genomic Signature of Selective Sweeps Illuminates Adaptation of Medicago truncatula to Root-Associated Microorganisms.

Medicago truncatula is a model legume species used to investigate plant-microorganism interactions, notably root symbioses. Massive population genomic and transcriptomic data now available for this species open the way for a comprehensive investigation of genomic variations associated with adaptation of M. truncatula to its environment. Here we performed a fine-scale genome scan of selective sweep signatures in M. truncatula using more than 15 million single nucleotide polymorphisms identified on 283 accessions from two populations (Circum and Far West), and exploited annotation and published transcriptomic data to identify biological processes associated with molecular adaptation. We identified 58 swept genomic regions with a 15 kb average length and comprising 3.3 gene models on average. The unimodal sweep state probability distribution in these regions enabled us to focus on the best single candidate gene per region. We detected two unambiguous species-wide selective sweeps, one of which appears to underlie morphological adaptation. Population genomic analyses of the remaining 56 sweep signatures indicate that sweeps identified in the Far West population are less population-specific and probably more ancient than those identified in the Circum population. Functional annotation revealed a predominance of immunity-related adaptations in the Circum population. Transcriptomic data from accessions of the Far West population allowed inference of four clusters of coregulated genes putatively involved in the adaptive control of symbiotic carbon flow and nodule senescence, as well as in other root adaptations upon infection with soil microorganisms. We demonstrate that molecular adaptations in M. truncatula were primarily triggered by selective pressures from root-associated microorganisms.

Cell wall proteome of sugarcane stems: comparison of a destructive and a non-destructive extraction method showed differences in glycoside hydrolases and peroxidases.

BACKGROUND: Sugarcane has been used as the main crop for ethanol production for more than 40 years in Brazil. Recently, the production of bioethanol from bagasse and straw, also called second generation (2G) ethanol, became a reality with the first commercial plants started in the USA and Brazil. However, the industrial processes still need to be improved to generate a low cost fuel. One possibility is the remodeling of cell walls, by means of genetic improvement or transgenesis, in order to make the bagasse more accessible to hydrolytic enzymes. We aimed at characterizing the cell wall proteome of young sugarcane culms, to identify proteins involved in cell wall biogenesis. Proteins were extracted from the cell walls of 2-month-old culms using two protocols, non-destructive by vacuum infiltration vs destructive. The proteins were identified by mass spectrometry and bioinformatics. RESULTS: A predicted signal peptide was found in 84 different proteins, called cell wall proteins (CWPs). As expected, the non-destructive method showed a lower percentage of proteins predicted to be intracellular than the destructive one (33% vs 44%). About 19% of CWPs were identified with both methods, whilst the infiltration protocol could lead to the identification of 75% more CWPs. In both cases, the most populated protein functional classes were those of proteins related to lipid metabolism and oxido-reductases. Curiously, a single glycoside hydrolase (GH) was identified using the non-destructive method whereas 10 GHs were found with the destructive one. Quantitative data analysis allowed the identification of the most abundant proteins. CONCLUSIONS: The results highlighted the importance of using different protocols to extract proteins from cell walls to expand the coverage of the cell wall proteome. Ten GHs were indicated as possible targets for further studies in order to obtain cell walls less recalcitrant to deconstruction. Therefore, this work contributed to two goals: enlarge the coverage of the sugarcane cell wall proteome, and provide target proteins that could be used in future research to facilitate 2G ethanol production.

Genome of an arbuscular mycorrhizal fungus provides insight into the oldest plant symbiosis.

The mutualistic symbiosis involving Glomeromycota, a distinctive phylum of early diverging Fungi, is widely hypothesized to have promoted the evolution of land plants during the middle Paleozoic. These arbuscular mycorrhizal fungi (AMF) perform vital functions in the phosphorus cycle that are fundamental to sustainable crop plant productivity. The unusual biological features of AMF have long fascinated evolutionary biologists. The coenocytic hyphae host a community of hundreds of nuclei and reproduce clonally through large multinucleated spores. It has been suggested that the AMF maintain a stable assemblage of several different genomes during the life cycle, but this genomic organization has been questioned. Here we introduce the 153-Mb haploid genome of Rhizophagus irregularis and its repertoire of 28,232 genes. The observed low level of genome polymorphism (0.43 SNP per kb) is not consistent with the occurrence of multiple, highly diverged genomes. The expansion of mating-related genes suggests the existence of cryptic sex-related processes. A comparison of gene categories confirms that R. irregularis is close to the Mucoromycotina. The AMF obligate biotrophy is not explained by genome erosion or any related loss of metabolic complexity in central metabolism, but is marked by a lack of genes encoding plant cell wall-degrading enzymes and of genes involved in toxin and thiamine synthesis. A battery of mycorrhiza-induced secreted proteins is expressed in symbiotic tissues. The present comprehensive repertoire of R. irregularis genes provides a basis for future research on symbiosis-related mechanisms in Glomeromycota.

Comprehensive genome-wide analysis of the Aux/IAA gene family in Eucalyptus: evidence for the role of EgrIAA4 in wood formation.

Auxin plays a pivotal role in various plant growth and development processes, including vascular differentiation. The modulation of auxin responsiveness through the auxin perception and signaling machinery is believed to be a major regulatory mechanism controlling cambium activity and wood formation. To gain more insights into the roles of key Aux/IAA gene regulators of the auxin response in these processes, we identified and characterized members of the Aux/IAA family in the genome of Eucalyptus grandis, a tree of worldwide economic importance. We found that the gene family in Eucalyptus is slightly smaller than that in Populus and Arabidopsis, but all phylogenetic groups are represented. High-throughput expression profiling of different organs and tissues highlighted several Aux/IAA genes expressed in vascular cambium and/or developing xylem, some showing differential expression in response to developmental (juvenile vs. mature) and/or to environmental (tension stress) cues. Based on the expression profiles, we selected a promising candidate gene, EgrIAA4, for functional characterization. We showed that EgrIAA4 protein is localized in the nucleus and functions as an auxin-responsive repressor. Overexpressing a stabilized version of EgrIAA4 in Arabidopsis dramatically impeded plant growth and fertility and induced auxin-insensitive phenotypes such as inhibition of primary root elongation, lateral root emergence and agravitropism. Interestingly, the lignified secondary walls of the interfascicular fibers appeared very late, whereas those of the xylary fibers were virtually undetectable, suggesting that EgrIAA4 may play crucial roles in fiber development and secondary cell wall deposition.