Physiological and proteomics analyses of Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.) responses to Phytophthora cinnamomi

Proteomic and Metabolomic Analysis of the Quercus ilex-Phytophthora cinnamomi Pathosystem Reveals a Population-Specific Response, Independent of Co-Occurrence of Drought

Holm oak (Quercus ilex) is considered to be one of the major structural elements of Mediterranean forests and the agrosilvopastoral Spanish "dehesa", making it an outstanding example of ecological and socioeconomic sustainability in forest ecosystems. The exotic Phytophthora cinnamomi is one of the most aggressive pathogens of woody species and, together with drought, is considered to be one of the main drivers of holm oak decline. The effect of and response to P. cinnamomi inoculation were studied in the offspring of mother trees from two Andalusian populations, Cordoba and Huelva. At the two locations, acorns collected from both symptomatic (damaged) and asymptomatic (apparently healthy) trees were sampled. Damage symptoms, mortality, and chlorophyll fluorescence were evaluated in seedlings inoculated under humid and drought conditions. The effect and response depended on the population and were more apparent in Huelva than in Cordoba. An integrated proteomic and metabolomic analysis revealed the involvement of different metabolic pathways in response to the pathogen in both populations, including amino acid metabolism pathways in Huelva, and terpenoid and flavonoid biosynthesis in Cordoba. However, no differential response was observed between seedlings inoculated under humid and drought conditions. A protective mechanism of the photosynthetic apparatus was activated in response to defective photosynthetic activity in inoculated plants, which seemed to be more efficient in the Cordoba population. In addition, enzymes and metabolites of the phenylpropanoid and flavonoid biosynthesis pathways may have conferred higher resistance in the Cordoba population. Some enzymes are proposed as markers of resilience, among which glyoxalase I, glutathione reductase, thioredoxin reductase, and cinnamyl alcohol dehydrogenase are candidates.

Silicon modifies leaf nutriome and improves growth of oak seedlings exposed to phosphorus deficiency and Phytophthora plurivora infection

Beneficial effects of silicon (Si) on plants have primarily been studied in crop species under single stress. Moreover, nutrient acquisition-based responses to combination of biotic and abiotic stresses (a common situation in natural habitats) have rarely been reported, in particular in conjunction with soil amendments with Si. Pedunculate oak (Quercus robur L.), one of the ecologically and economically most important tree species in Europe, is facing a severe decline due to combined stresses, but also problems in assisted regeneration in nurseries. Here, we studied the effect of Si supply on the leaf nutriome, root traits and overall growth of 12-weeks-old oak seedlings exposed to abiotic stress [low phosphorus (P) supply], biotic stress (Phytophthora plurivora root infection), and their combination. The application of Si had the strongest ameliorative effect on growth, root health and root phenome under the most severe stress conditions (i.e., combination of P deficiency and P. plurivora root infection), where it differentially affected the uptake and leaf accumulation in 11 out of 13 analysed nutrients. Silicon supply tended to reverse the pattern of change of some, but not all, leaf nutrients affected by stresses: P, boron (B) and magnesium (Mg) under P deficiency, and P, B and sulphur (S) under pathogen attack, but also nickel (Ni) and molybdenum (Mo) under all three stresses. Surprisingly, Si affected some nutrients that were not changed by a particular stress itself and decreased leaf Mg levels under all the stresses. On the other hand, pathogen attack increased leaf accumulation of Si. This exploratory work presents the complexity of nutrient crosstalk under three stresses, and opens more questions about genetic networks that control plant physiological responses. Practically, we show a potential of Si application to improve P status and root health in oak seedlings, particularly in nurseries.

Proteomics research in forest trees: A 2012-2022 update

This review is a compilation of proteomic studies on forest tree species published in the last decade (2012-2022), mostly focused on the most investigated species, including Eucalyptus, Pinus, and Quercus. Improvements in equipment, platforms, and methods in addition to the increasing availability of genomic data have favored the biological knowledge of these species at the molecular, organismal, and community levels. Integration of proteomics with physiological, biochemical and other large-scale omics in the direction of the Systems Biology, will provide a comprehensive understanding of different biological processes, from growth and development to responses to biotic and abiotic stresses. As main issue we envisage that proteomics in long-living plants will thrive light on the plant responses and resilience to global climate change, contributing to climate mitigation strategies and molecular breeding programs. Proteomics not only will provide a molecular knowledge of the mechanisms of resilience to either biotic or abiotic stresses, but also will allow the identification on key gene products and its interaction. Proteomics research has also a translational character being applied to the characterization of the variability and biodiversity, as well as to wood and non-wood derived products, traceability, allergen and bioactive peptides identification, among others. Even thought, the full potential of proteomics is far from being fully exploited in forest tree research, with PTMs and interactomics being reserved to plant model systems. The most outstanding achievements in forest tree proteomics in the last decade as well as prospects are discussed.

Multiomics Molecular Research into the Recalcitrant and Orphan Quercus ilex Tree Species: Why, What for, and How

The holm oak (Quercus ilex L.) is the dominant tree species of the Mediterranean forest and the Spanish agrosilvopastoral ecosystem, “dehesa.” It has been, since the prehistoric period, an important part of the Iberian population from a social, cultural, and religious point of view, providing an ample variety of goods and services, and forming the basis of the economy in rural areas. Currently, there is renewed interest in its use for dietary diversification and sustainable food production. It is part of cultural richness, both economically (tangible) and environmentally (intangible), and must be preserved for future generations. However, a worrisome degradation of the species and associated ecosystems is occurring, observed in an increase in tree decline and mortality, which requires urgent action. Breeding programs based on the selection of elite genotypes by molecular markers is the only plausible biotechnological approach. To this end, the authors’ group started, in 2004, a research line aimed at characterizing the molecular biology of Q. ilex. It has been a challenging task due to its biological characteristics (long life cycle, allogamous, high phenotypic variability) and recalcitrant nature. The biology of this species has been characterized following the central dogma of molecular biology using the omics cascade. Molecular responses to biotic and abiotic stresses, as well as seed maturation and germination, are the two main objectives of our research. The contributions of the group to the knowledge of the species at the level of DNA-based markers, genomics, epigenomics, transcriptomics, proteomics, and metabolomics are discussed here. Moreover, data are compared with those reported for Quercus spp. All omics data generated, and the genome of Q. ilex available, will be integrated with morphological and physiological data in the systems biology direction. Thus, we will propose possible molecular markers related to resilient and productive genotypes to be used in reforestation programs. In addition, possible markers related to the nutritional value of acorn and derivate products, as well as bioactive compounds (peptides and phenolics) and allergens, will be suggested. Subsequently, the selected molecular markers will be validated by both genome-wide association and functional genomic analyses.

Protein Markers for the Identification of Cork Oak Plants Infected with Phytophthora cinnamomi by Applying an (α, β)-k-Feature Set Approach

Cork oak decline in Mediterranean forests is a complex phenomenon, observed with remarkable frequency in the southern part of the Iberian Peninsula, causing the weakening and death of these woody plants. The defoliation of the canopy, the presence of dry peripheral branches, and exudations on the trunk are visible symptoms used for the prognosis of decline, complemented by the presence of Phytophthora cinnamomi identified in the rhizosphere of the trees and adjacent soils. Recently, a large proteomic dataset obtained from the leaves of cork oak plants inoculated and non-inoculated with P. cinnamomi has become available. We explored it to search for an optimal set of proteins, markers of the biological pattern of interaction with the oomycete. Thus, using published data from the cork oak leaf proteome, we mathematically modelled the problem as an α, β-k-Feature Set Problem to select molecular markers. A set of proteins (features) that represent dominant effects on the host metabolism resulting from pathogen action on roots was found. These results contribute to an early diagnosis of biochemical changes occurring in cork oak associated with P. cinnamomi infection. We hypothesize that these markers may be decisive in identifying trees that go into decline due to interactions with the pathogen, assisting the management of cork oak forest ecosystems.

Metabolomic and Physiological Changes in Fagus sylvatica Seedlings Infected with Phytophthora plurivora and the A1 and A2 Mating Types of P. ×cambivora

Phytophthora infections are followed by histological alterations, physiological and metabolomic adjustments in the host but very few studies contemplate these changes simultaneously. Fagus sylvatica seedlings were inoculated with A1 and A2 mating types of the heterothallic P. ×cambivora and with the homothallic P. plurivora to identify plant physiological and metabolomic changes accompanying microscope observations of the colonization process one, two and three weeks after inoculation. Phytophthora plurivora-infected plants died at a faster pace than those inoculated with P. ×cambivora and showed higher mortality than P. ×cambivora A1-infected plants. Phytophthora ×cambivora A1 and A2 caused similar progression and total rate of mortality. Most differences in the physiological parameters between inoculated and non-inoculated plants were detected two weeks after inoculation. Alterations in primary and secondary metabolites in roots and leaves were demonstrated for all the inoculated plants two and three weeks after inoculation. The results indicate that P. plurivora is more aggressive to Fagus sylvatica seedlings than both mating types of P. ×cambivora while P. ×cambivora A1 showed a slower infection mode than P. ×cambivora A2 and led to minor plant metabolomic adjustments.

Micropropagation, Characterization, and Conservation of Phytophthora cinnamomi-Tolerant Holm Oak Mature Trees

Holm oak populations have deteriorated drastically due to oak decline syndrome. The first objective of the present study was to investigate the use of axillary budding and somatic embryogenesis (SE) to propagate asymptomatic holm oak genotypes identified in disease hotspots in Spain. Axillary budding was achieved in two out of six tolerant genotypes from the south-western region and in two out of four genotypes from the Mediterranean region. Rooting of shoots cultured on medium supplemented with 3 mg L−1 of indole-3-acetic acid plus 0.1 mg L−1 α-naphthalene acetic acid was achieved, with rates ranging from 8 to 36%. Shoot cultures remained viable after cold storage for 9–12 months; this procedure is therefore suitable for medium-term conservation of holm oak germplasm. SE was induced in two out of the three genotypes tested, by using nodes and shoot tips cultured in medium without plant growth regulators. In vitro cloned progenies of the tolerant genotypes PL-T2 and VA5 inhibited growth of Phytophthora cinnamomi mycelia when exposed to the oomycete in vitro. Significant differences in total phenol contents and in the expression profiles of genes regulating phenylpropanoid biosynthesis were observed between in vitro cultured shoots derived from tolerant trees and cultures established from control genotypes.

Effect and Response of Quercus ilex subsp. ballota [Desf.] Samp. Seedlings From Three Contrasting Andalusian Populations to Individual and Combined Phytophthora cinnamomi and Drought Stresses

Quercus ilex L. is the dominant species in the Mediterranean forest and agrosilvopastoral ecosystem "dehesa." Currently, this forest species is threatened by natural and anthropogenic agents, especially by the decline syndrome, which is caused by Phytophthora cinnamomi and drought periods. Although the morphological and physiological responses of Q. ilex to combined stress (P. cinnamomi and drought) have been examined already, little is known at the molecular level. In this study, we studied the effect and response of 8-month seedlings from three contrasting Andalusian populations (Seville [Se], Granada [Gr], and Almeria [Al]) to the individual and combined stresses of P. cinnamomi and drought from morphological, physiological, biochemical, and proteomics data. Whereas, seedling damage (leaf chlorosis and necrosis) and mortality were greater under the combined stresses in the three populations, the effect of each individual stress was population-dependent. Resilient individuals were found in all the populations at different percentages. The decrease in leaf chlorophyll fluorescence, photosynthetic activity, and stomatal conductance observed in undamaged seedlings was greater in the presence of both stresses, the three populations responding similarly to drought and P. cinnamomi. Biochemical and proteomic analyses of undamaged seedlings from the two most markedly contrasting populations (Se and Al) revealed the absence of significant differences in the contents in photosynthetic pigments, amino acids, and phenolics among treatments. The Se and Al populations exhibited changes in protein profile in response to the different treatments, with 83 variable proteins in the former population and 223 in the latter. Variable proteins belonged to 16 different functional groups, the best represented among which were protein folding, sorting and degradation, carbohydrate, amino acid, and secondary metabolism, photosynthesis, and ROS scavenging. While photosynthetic proteins were mainly downaccumulated, those of stress-responsive were upaccumulated. Although no treatment-specific response was observed in any functional group, differences in abundance were especially marked under the combined stresses. The following variable proteins are proposed as putative markers for resilience in Q. ilex, namely, aldehyde dehydrogenase, glucose-6-phosphate isomerase, 50S ribosomal protein L5, and α-1,4-glucan-protein synthase [UDP-forming].

Proteomics and plant biology: contributions to date and a look towards the next decade

INTRODUCTION

This review presents the view of the author, that is opinionable and even speculative, on the field of proteomics, its application to plant biology knowledge, and translation to biotechnology. Written in a more academic than scientific style, it is based on past original and review articles by the author´s group, and those published by leading scientists in the last two years.

AREAS COVERED

Starting with a general definition and references to historical milestones, it covers sections devoted to the different platforms employed, the plant biology discourse in the protein language, challenges and future prospects, ending with the author opinion.

EXPERT OPINION

In 25 years, five proteomics platform generations have appeared. We are now moving from proteomics to Systems Biology. While feasible with model organisms, proteomics of orphan species remains challenging. Proteomics, even in its simplest approach, sheds light on plant biological processes, central dogma, and molecular bases of phenotypes of interest, and it can be translated to areas such as food traceability and allergen detection. Proteomics should be validated and optimized to each experimental system, objectives, and hypothesis. It has limitations, artifacts, and biases. We should not blindly accept proteomics data and just create a list of proteins, networks, and avoid speculative biological interpretations. From the hundred to thousand proteins identified and quantified, it is important to obtain a focus and validate some of them, otherwise it is merely. We are starting to have the protein pieces, so let, from now, build the proteomics and biological puzzle.

Proteomics Data Analysis for the Identification of Proteins and Derived Proteotypic Peptides of Potential Use as Putative Drought Tolerance Markers for Quercus ilex

Drought is one of the main causes of mortality in holm oak (Quercus ilex) seedlings used in reforestation programs. Although this species shows high adaptability to the extreme climate conditions prevailing in Southern Spain, its intrinsic genetic variability may play a role in the differential response of some populations and individuals. The aim of this work was to identify proteins and derived proteotypic peptides potentially useful as putative markers for drought tolerance in holm oak by using a targeted post-acquisition proteomics approach. For this purpose, we used a set of proteins identified by shotgun (LC-MSMS) analysis in a drought experiment on Q. ilex seedlings from four different provenances (viz. the Andalusian provinces Granada, Huelva, Cadiz and Seville). A double strategy involving the quantification of proteins and target peptides by shotgun analysis and post-acquisition data analysis based on proteotypic peptides was used. To this end, an initial list of proteotypic peptides from proteins highly represented under drought conditions was compiled that was used in combination with the raw files from the shotgun experiment to quantify the relative abundance of the fragment’s ion peaks with the software Skyline. The most abundant peptides under drought conditions in at least two populations were selected as putative markers of drought tolerance. A total of 30 proteins and 46 derived peptides belonging to the redox, stress-related, synthesis,-folding and degradation, and primary and secondary metabolism functional groups were thus identified. Two proteins (viz., subtilisin and chaperone GrpE protein) were found at increased levels in three populations, which make them especially interesting for validation drought tolerance markers in subsequent experiments.

Molecular Research on Stress Responses in Quercus spp.: From Classical Biochemistry to Systems Biology through Omics Analysis

The genus Quercus (oak), family Fagaceae, comprises around 500 species, being one of the most important and dominant woody angiosperms in the Northern Hemisphere. Nowadays, it is threatened by environmental cues, which are either of biotic or abiotic origin. This causes tree decline, dieback, and deforestation, which can worsen in a climate change scenario. In the 21st century, biotechnology should take a pivotal role in facing this problem and proposing sustainable management and conservation strategies for forests. As a non-domesticated, long-lived species, the only plausible approach for tree breeding is exploiting the natural diversity present in this species and the selection of elite, more resilient genotypes, based on molecular markers. In this direction, it is important to investigate the molecular mechanisms of the tolerance or resistance to stresses, and the identification of genes, gene products, and metabolites related to this phenotype. This research is being performed by using classical biochemistry or the most recent omics (genomics, epigenomics, transcriptomics, proteomics, and metabolomics) approaches, which should be integrated with other physiological and morphological techniques in the Systems Biology direction. This review is focused on the current state-of-the-art of such approaches for describing and integrating the latest knowledge on biotic and abiotic stress responses in Quercus spp., with special reference to Quercus ilex, the system on which the authors have been working for the last 15 years. While biotic stress factors mainly include fungi and insects such as Phytophthora cinnamomi, Cerambyx welensii, and Operophtera brumata, abiotic stress factors include salinity, drought, waterlogging, soil pollutants, cold, heat, carbon dioxide, ozone, and ultraviolet radiation. The review is structured following the Central Dogma of Molecular Biology and the omic cascade, from DNA (genomics, epigenomics, and DNA-based markers) to metabolites (metabolomics), through mRNA (transcriptomics) and proteins (proteomics). An integrated view of the different approaches, challenges, and future directions is critically discussed.

Mass spectrometry-based forest tree metabolomics

Research in forest tree species has advanced slowly when compared with other agricultural crops and model organisms, mainly due to the long-life cycles, large genome sizes, and lack of genomic tools. Additionally, trees are complex matrices, and the presence of interferents (e.g., oleoresins and cellulose) challenges the analysis of tree tissues with mass spectrometry (MS)-based analytical platforms. In this review, advances in MS-based forest tree metabolomics are discussed. Given their economic and ecological significance, particular focus is given to Pinus, Quercus, and Eucalyptus forest tree species to better understand their metabolite responses to abiotic and biotic stresses in the current climate change scenario. Furthermore, MS-based metabolomics technologies produce large and complex datasets that require expertize to adequately manage, process, analyze, and store the data in dedicated repositories. To ensure that the full potential of forest tree metabolomics data are translated into new knowledge, these data should comply with the FAIR principles (i.e., Findable, Accessible, Interoperable, and Re-usable). It is essential that adequate standards are implemented to annotate metadata from forest tree metabolomics studies as is already required by many science and governmental agencies and some major scientific publishers. © 2019 John Wiley & Sons Ltd. Mass Spec Rev 40:126-157, 2021.

Disclosing proteins in the leaves of cork oak plants associated with the immune response to Phytophthora cinnamomi inoculation in the roots: A long-term proteomics approach

The pathological interaction between oak trees and Phytophthora cinnamomi has implications in the cork oak decline observed over the last decades in the Iberian Peninsula. During host colonization, the phytopathogen secretes effector molecules like elicitins to increase disease effectiveness. The objective of this study was to unravel the proteome changes associated with the cork oak immune response triggered by P. cinnamomi inoculation in a long-term assay, through SWATH-MS quantitative proteomics performed in the oak leaves. Using the Arabidopis proteome database as a reference, 424 proteins were confidently quantified in cork oak leaves, of which 80 proteins showed a p-value below 0.05 or a fold-change greater than 2 or less than 0.5 in their levels between inoculated and control samples being considered as altered. The inoculation of cork oak roots with P. cinnamomi increased the levels of proteins associated with protein-DNA complex assembly, lipid oxidation, response to endoplasmic reticulum stress, and pyridine-containing compound metabolic process in the leaves. In opposition, several proteins associated with cellular metabolic compound salvage and monosaccharide catabolic process had significantly decreased abundances. The most significant abundance variations were observed for the Ribulose 1,5-Bisphosphate Carboxylase small subunit (RBCS1A), Heat Shock protein 90–1 (Hsp90-1), Lipoxygenase 2 (LOX2) and Histone superfamily protein H3.3 (A8MRLO/At4G40030) revealing a pertinent role for these proteins in the host-pathogen interaction mechanism. This work represents the first SWATH-MS analysis performed in cork oak plants inoculated with P. cinnamomi and highlights host proteins that have a relevant action in the homeostatic states that emerge from the interaction between the oomycete and the host in the long term and in a distal organ.

The Impacts of Drought Stress and Phytophthora cinnamomi Infection on Short-Term Water Relations in Two Year-Old Eucalyptus obliqua

The effects of drought stress, Phytophthora cinnamomi infection and their interaction on water relations and growth were examined for 28 days on two year-old potted trees of Eucalyptus obliqua (L’Hér.). There were significant effects of drought stress on plant photosynthesis, stomatal conductance, biomass accumulation, plant water potential at turgor loss point and the bulk modulus of elasticity. E. obliqua was successfully infected but the trees showed only mild symptoms. Infection with P. cinnamomi led to a significant reduction in the root biomass and root-to-shoot ratio in well-watered and droughted plants but did not impact water relations. There was no observable cumulative effect of drought and P. cinnamomi infection. There are multiple potential reasons why P. cinnamomi infection did not lead to drought-like symptoms in E. obliqua, including short experimental duration, delayed infection symptoms, potential resistance of E. obliqua and a possible lower aggressiveness of the P. cinnamomi strain. Hence, our results indicate that P. cinnamomi infection will not always lead to immediate short-term symptoms, and that plants that are mildly symptomatic respond very similar to drought stress compared to non-infected trees.

Holm Oak Somatic Embryogenesis: Current Status and Future Perspectives

Quercus ilex (holm oak) is one of the most representative trees in the Mediterranean basin, but now the sustainability of its ecosystems is at serious risk due to the lack of natural regeneration and to the presence of a severe disease called oak decline that has caused the death of thousands of trees. The application of biotechnological tools, such as somatic embryogenesis, allows programs of genetic improvement of the species to be speeded up and helps in the conservation of its ecosystems. Somatic embryogenesis is currently considered one of the main biotechnological techniques that has demonstrated significant benefits when has applied to forest tree species, providing significant advantages such as mass propagation, genetic transformation, application of synthetic seed, and cryopreservation of elite genotypes. In this report, the state of the art of somatic embryogenesis in holm oak is reviewed. Factors affecting the induction (plant donor age, type of explant, or plant growth regulators) and maintenance and proliferation of the embryogenic cultures are summarized. Advances on the conversion of somatic embryos into plants and on the acclimatization of these plantlets, as well as the results obtained on the application of the genetic transformation and the cryopreservation procedures to holm oak embryogenic cultures, are also presented.

Proteomics, Holm Oak (Quercus ilex L.) and Other Recalcitrant and Orphan Forest Tree Species: How do They See Each Other?

Proteomics has had a big impact on plant biology, considered as a valuable tool for several forest species, such as Quercus, Pines, Poplars, and Eucalyptus. This review assesses the potential and limitations of the proteomics approaches and is focused on Quercus ilex as a model species and other forest tree species. Proteomics has been used with Q. ilex since 2003 with the main aim of examining natural variability, developmental processes, and responses to biotic and abiotic stresses as in other species of the genus Quercus or Pinus. As with the progress in techniques in proteomics in other plant species, the research in Q. ilex moved from 2-DE based strategy to the latest gel-free shotgun workflows. Experimental design, protein extraction, mass spectrometric analysis, confidence levels of qualitative and quantitative proteomics data, and their interpretation are a true challenge with relation to forest tree species due to their extreme orphan and recalcitrant (non-orthodox) nature. Implementing a systems biology approach, it is time to validate proteomics data using complementary techniques and integrate it with the -omics and classical approaches. The full potential of the protein field in plant research is quite far from being entirely exploited. However, despite the methodological limitations present in proteomics, there is no doubt that this discipline has contributed to deeper knowledge of plant biology and, currently, is increasingly employed for translational purposes.

Regulation by biotic stress of tannins biosynthesis in Quercus ilex: Crosstalk between defoliation and Phytophthora cinnamomi infection

Sustainability of the Mediterranean forest is threatened by oak decline, a disease of holm oak and other Quercus species that is initiated by infection with the oomycete Phytophthora cinnamomi. Focusing on the role of tannins in the chemical defense of plants, this work investigated whether tannins content in Quercus ilex is regulated by biotic stress. Screening of published genomes allowed the identification of Quercus sequences encoding enzymes for early steps of the biosynthesis of phenolic compounds, like hydrolysable tannins and condensed tannins (CT) among others, plus genes involved in the late steps of CT biosynthesis. Four days after treatment of Q. ilex seedlings by mechanical defoliation, P. cinnamomi infection and both stressors simultaneously, mRNA concentrations for tannins biosynthesis enzymes were measured in leaves. Among the transcript amount for shikimate dehydrogenase (SDH, EC 1.1.1.25), anthocyanidin reductase (EC 1.3.1.77), anthocyanidin synthase (EC 1.14.11.19) and leucoanthocyanidine reductase (EC 1.17.1.3), defoliation induced gene expression for SDH2 isoenzyme. About 4 days after infection of roots by P. cinnamomi, this up-regulation was canceled and SDH enzyme activity decreased. Furthermore, during this late stage of biotrophic interaction the pathogen switched off the correlation engaged by defoliation between the expression of SDH1 and SDH2 encoding genes and chemical defenses corresponding to total tannins, which were down-regulated. Thus, tannins biosynthesis in seedlings of Q. ilex is induced after mechanical defoliation whereas infection by the pathogen interferes with this regulation, potentially increasing the susceptibility of plants to herbivory and aggravating the impact of biotic stress.

Ion Torrent and lllumina, two complementary RNA-seq platforms for constructing the holm oak (Quercus ilex) transcriptome

Transcriptome analysis is widely used in plant biology research to explore gene expression across a large variety of biological contexts such as those related to environmental stress and plant-pathogen interaction. Currently, next generation sequencing platforms are used to obtain a high amount of raw data to build the transcriptome of any plant. Here, we compare Illumina and Ion Torrent sequencing platforms for the construction and analysis of the holm oak (Quercus ilex) transcriptome. Genomic analysis of this forest tree species is a major challenge considering its recalcitrant character and the absence of previous molecular studies. In this study, Quercus ilex raw sequencing reads were obtained from Illumina and Ion Torrent and assembled by three different algorithms, MIRA, RAY and TRINITY. A hybrid transcriptome combining both sequencing technologies was also obtained in this study. The RAY-hybrid assembly generated the most complete transcriptome (1,116 complete sequences of which 1,085 were single copy) with a E90N50 of 1,122 bp. The MIRA-Illumina and TRINITY-Ion Torrent assemblies annotated the highest number of total transcripts (62,628 and 74,058 respectively). MIRA-Ion Torrent showed the highest number of shared sequences (84.8%) with the oak transcriptome. All the assembled transcripts from the hybrid transcriptome were annotated with gene ontology grouping them in terms of biological processes, molecular functions and cellular components. In addition, an in silico proteomic analysis was carried out using the translated assemblies as databases. Those from Ion Torrent showed more proteins compared to the Illumina and hybrid assemblies. This new generated transcriptome represents a valuable tool to conduct differential gene expression studies in response to biotic and abiotic stresses and to assist and validate the ongoing Q. ilex whole genome sequencing.

Toward characterizing germination and early growth in the non-orthodox forest tree species Quercus ilex through complementary gel and gel-free proteomic analysis of embryo and seedlings

By using two complementary proteomics, gel-based and gel-free (shotgun) approaches, the protein profiles of the non-orthodox forest tree species Quercus ilex seeds during germination and early seedling growth have been compared. Proteins were extracted from embryo axis, radicle and shoot tissues at different developmental stages. Proteins were subjected to one- and two-dimensional gel electrophoresis. A multivariate analysis (PCA) revealed that SDS-PAGE clearly separated germination (0-24 h post-imbibition), postgermination (72-216 h post-imbibition) and early seedling growth stages (2 weeks post-imbibition). Image analysis of the two-dimensional gels revealed a total of 732 spots, 103 of which were significantly variable among developmental stages. After MALDI-TOF/TOF MS analysis, 90 spots were identified, belonging to six main functional categories: carbohydrate, amino acids, energy, and protein metabolism, biosynthesis of secondary metabolites, and redox processes. The gel-based approach disclosed important metabolic changes that occurred in the holm oak seed after the germination. However, few proteins were significantly altered during the germination period (from 0 h to 24 h post imbibition) and, because of that, a further shotgun analysis was therefore used to analyse changes in the protein profile during seed germination. Up to 1250 proteins could be confidently identified, with 153 being variable. They belonged to the main functional categories of carbohydrate, amino acids and secondary metabolism, protein degradation, and responses to abiotic stress. The accumulation of proteases and amino acids metabolism proteins in mature seeds can be reflecting the production of energy from the mobilization of storage proteins to start germination. These results, therefore, corroborate the hypothesis that the mature non-orthodox seeds of Q. ilex have all the machinery necessary for rapidly resuming metabolic activities and starting the germination process, in contrast to that occurs in orthodox seeds, which metabolic activity ceases in mature dry seeds. The use of a genus-specific database combined with the public Viridiplantae database improved the quality and quantity of protein identification in this orphan species. In addition, both proteomics approaches (gel-based and shotgun) were complementary, with shotgun increasing by over two-fold the coverage of the proteome analysed. Both approaches provided similar results and supported the same conclusions on the metabolic switch experienced by the seed upon germination. SIGNIFICANCE: The optimal seed germination is a prerequisite for successful seedling establishment and plant vigour, being of great relevance in the case of crops and commercial woody plants. By using a complementary gel-based and gel-free proteomic strategy we have study the protein profiles of the non-orthodox forest tree species Quercus ilex seeds during germination and early seedling growth. The contribution of this work is of great importance, due to the complemented proteomic approaches giving similar clues to the metabolic state of the mature Q. ilex seed before the germination starts, and the metabolic switch experienced by the imbibed acorn until the seedling is established.

Holm oak proteomic response to water limitation at seedling establishment stage reveals specific changes in different plant parts as well as interaction between roots and cotyledons

Quercus ilex is a dominant tree species in the Mediterranean region with double economic and ecological importance and increasing use in reforestation. Seedling establishment is extremely vulnerable to environmental stresses, particularly drought. A time course study on physiological and proteomic response of holm oak to water limitation stress and recovery during early heterotrophic growth is reported. Applied stress led to diminution in plant water content and root growth, oxidative stress in roots and some alterations in the anti-oxidative protection. Plant parts differed substantially in soluble sugar and free phenolic content, and in their changes during stress and recovery. Proteomic response in holm oak roots and cotyledons was estimated using combined 1-DE/2-DE approach and protein identification by MALDI TOF-TOF PMF and MS/MS. A total of 127 differentially abundant protein species (DAPs) were identified. DAPs related to starch metabolism, lipid to sugar conversion, reserve proteins and their mobilization were typical for cotyledons. DAPs in roots were involved in sugar utilization, secondary metabolism and defense, including pathogenesis related proteins from PR-5 and PR-10 families. Results emphasize specific proteome signatures of separate plant parts as well as importance of sink-source interaction between root and cotyledon in the time course of stress and in recovery.

Differences in the Response to Acute Drought and Phytophthora cinnamomi Rands Infection in Quercus ilex L. Seedlings

The sustainability of “dehesas” is threatened by the Holm oak decline. It is thought that the effects of root rot on plant physiology vary depending on external stress factors. Plant growth and biomass allocation are useful tools to characterize differences in the response to drought and infection. The study of physiological responses together with growth patterns will clarify how and to what extent root rot is able to damage the plant. A fully factorial experiment, including drought and Phytophtora cinnamomi Rands infection as factors, was carried out with Quercus ilex L. seedlings. Photosynthesis, biomass allocation and root traits were assessed. Photosynthetic variables responded differently to drought and infection over time. The root mass fraction showed a significant reduction due to infection. P. cinnamomi root rot altered the growth patterns. Plants could not recover from the physiological effects of infection only when the root rot coincided with water stress. Without additional stressors, the strategy of our seedlings in the face of root rot was to reduce the biomass increment and reallocate resources. Underlying mechanisms involved in plant-pathogen interactions should be considered in the study of holm oak decline, beyond the consideration of water stress as the primary cause of tree mortality.

A Multi-Omics Analysis Pipeline for the Metabolic Pathway Reconstruction in the Orphan Species Quercus ilex

Holm oak (Quercus ilex) is the most important and representative species of the Mediterranean forest and of the Spanish agrosilvo-pastoral "dehesa" ecosystem. Despite its environmental and economic interest, Holm oak is an orphan species whose biology is very little known, especially at the molecular level. In order to increase the knowledge on the chemical composition and metabolism of this tree species, the employment of a holistic and multi-omics approach, in the Systems Biology direction would be necessary. However, for orphan and recalcitrant plant species, specific analytical and bioinformatics tools have to be developed in order to obtain adequate quality and data-density before to coping with the study of its biology. By using a plant sample consisting of a pool generated by mixing equal amounts of homogenized tissue from acorn embryo, leaves, and roots, protocols for transcriptome (NGS-Illumina), proteome (shotgun LC-MS/MS), and metabolome (GC-MS) studies have been optimized. These analyses resulted in the identification of around 62629 transcripts, 2380 protein species, and 62 metabolites. Data are compared with those reported for model plant species, whose genome has been sequenced and is well annotated, including Arabidopsis, japonica rice, poplar, and eucalyptus. RNA and protein sequencing favored each other, increasing the number and confidence of the proteins identified and correcting erroneous RNA sequences. The integration of the large amount of data reported using bioinformatics tools allows the Holm oak metabolic network to be partially reconstructed: from the 127 metabolic pathways reported in KEGG pathway database, 123 metabolic pathways can be visualized when using the described methodology. They included: carbohydrate and energy metabolism, amino acid metabolism, lipid metabolism, nucleotide metabolism, and biosynthesis of secondary metabolites. The TCA cycle was the pathway most represented with 5 out of 10 metabolites, 6 out of 8 protein enzymes, and 8 out of 8 enzyme transcripts. On the other hand, gaps, missed pathways, included metabolism of terpenoids and polyketides and lipid metabolism. The multi-omics resource generated in this work will set the basis for ongoing and future studies, bringing the Holm oak closer to model species, to obtain a better understanding of the molecular mechanisms underlying phenotypes of interest (productive, tolerant to environmental cues, nutraceutical value) and to select elite genotypes to be used in restoration and reforestation programs, especially in a future climate change scenario.

Phytophthora cinnamomi

Phytophthora cinnamomi is one of the most devastating plant pathogens in the world. It infects close to 5000 species of plants, including many of importance in agriculture, forestry and horticulture. The inadvertent introduction of P. cinnamomi into natural ecosystems, including a number of recognized Global Biodiversity Hotspots, has had disastrous consequences for the environment and the biodiversity of flora and fauna. The genus Phytophthora belongs to the Class Oomycetes, a group of fungus-like organisms that initiate plant disease through the production of motile zoospores. Disease control is difficult in agricultural and forestry situations and even more challenging in natural ecosystems as a result of the scale of the problem and the limited range of effective chemical inhibitors. The development of sustainable control measures for the future management of P. cinnamomi requires a comprehensive understanding of the cellular and molecular basis of pathogen development and pathogenicity. The application of next-generation sequencing technologies to generate genomic and transcriptomic data promises to underpin a new era in P. cinnamomi research and discovery. The aim of this review is to integrate bioinformatic analyses of P. cinnamomi sequence data with current knowledge of the cellular and molecular basis of P. cinnamomi growth, development and plant infection. The goal is to provide a framework for future research by highlighting potential pathogenicity genes, shedding light on their possible functions and identifying suitable targets for future control measures.

TAXONOMY

Phytophthora cinnamomi Rands; Kingdom Chromista; Phylum Oomycota or Pseudofungi; Class Oomycetes; Order Peronosporales; Family Peronosporaceae; genus Phytophthora.

HOST RANGE

Infects about 5000 species of plants, including 4000 Australian native species. Host plants important for agriculture and forestry include avocado, chestnut, macadamia, oak, peach and pineapple.

DISEASE SYMPTOMS

A root pathogen which causes rotting of fine and fibrous roots, but which can also cause stem cankers. Root damage may inhibit water movement from roots to shoots, leading to dieback of young shoots. USEFUL WEBSITES: http://fungidb.org/fungidb/; http://genome.jgi.doe.gov/Phyci1/Phyci1.home.html; http://www.ncbi.nlm.nih.gov/assembly/GCA_001314365.1; http://www.ncbi.nlm.nih.gov/assembly/GCA_001314505.1.

Characterization of the orthodox Pinus occidentalis seed and pollen proteomes by using complementary gel-based and gel-free approaches

This work presents an analysis of Pinus occidentalis pollen and seed proteomes, in which both gel-based and gel-free approaches have been used. Proteins were extracted from P. occidentalis seeds and pollen by using the TCA/acetone/phenol precipitation protocol, and protein extracts were subjected to 1- and 2-DE coupled to MALDI-TOF-TOF as well as to shotgun (nLC-LTQ-Orbitrap) analysis. All bands (1-DE) and the most abundant spots (2-DE) were excised, trypsin digested and the resulting peptides analyzed by MALDI TOF/TOF. In order to increase the proteome coverage, a gel free approach was used. Proteins were identified from mass spectra by using three different databases, including UniProtKB, NCBI and a Pinus spp. custom database [2]. The gel-based approach resulted in 42 (seeds) and 94 (pollen) protein identifications, while the shotgun approach permitted the identification of 187 (seed) and 960 (pollen) proteins. Proteins were classified based on their corresponding functional categories. In seeds, storage proteins were the most abundant ones, and some allergens and proteases were also identified. In pollen proteins related to general metabolism were the most predominant. Data are compared and discussed from a methodological and biological point of view, taking into account the particularities of the seed and pollen organs.

BIOLOGICAL SIGNIFICANCE

In this work we characterized P. occidentalis proteome with seeds and pollen samples implementing two complementary approaches for the analysis. We found a high content of storage protein, stress response and metabolism related proteins in the seed proteome. Similarly, in the pollen proteome we found predominant groups of proteins related to metabolism and stress response.

2-DE proteomics analysis of drought treated seedlings of Quercus ilex supports a root active strategy for metabolic adaptation in response to water shortage

Holm oak is a dominant tree in the western Mediterranean region. Despite being well adapted to dry hot climate, drought is the main cause of mortality post-transplanting in reforestation programs. An active response to drought is critical for tree establishment and survival. Applying a gel-based proteomic approach, dynamic changes in root proteins of drought treated Quercus ilex subsp. Ballota [Desf.] Samp. seedlings were followed. Water stress was applied on 20 day-old holm oak plantlets by water limitation for a period of 10 and 20 days, each followed by 10 days of recovery. Stress was monitored by changes in water status, plant growth, and electrolyte leakage. Contrary to leaves, holm oak roots responded readily to water shortage at physiological level by growth inhibition, changes in water status and membrane stability. Root proteins were extracted using trichloroacetate/acetone/phenol protocol and separated by two-dimensional electrophoresis. Coomassie colloidal stained gel images were analyzed and spot intensity data subjected to multivariate statistical analysis. Selected consistent spots in three biological replicas, presenting significant changes under stress, were subjected to MALDI-TOF mass spectrometry (peptide mass fingerprinting and MS/MS). For protein identification, combined search was performed with MASCOT search engine over NCBInr Viridiplantae and Uniprot databases. Data are available via ProteomeXchange with identifier PXD002484. Identified proteins were classified into functional groups: metabolism, protein biosynthesis and proteolysis, defense against biotic stress, cellular protection against abiotic stress, intracellular transport. Several enzymes of the carbohydrate metabolism decreased in abundance in roots under drought stress while some related to ATP synthesis and secondary metabolism increased. Results point at active metabolic adjustment and mobilization of the defense system in roots to actively counteract drought stress.

Multiplex staining of 2-DE gels for an initial phosphoproteome analysis of germinating seeds and early grown seedlings from a non-orthodox specie: Quercus ilex L. subsp. ballota [Desf.] Samp

As a preliminary step in the phosphoproteome analysis of germinating seeds (0 and 24 h after seed imbibition) and early grown seedlings (216 h after seed imbibition) from a non-orthodox sp. Quercus ilex, a multiplex (SYPRO-Ruby and Pro-Q DPS) staining of high-resolution 2-DE gels was used. By using this protocol it was possible to detect changes in protein-abundance and/or phosphorylation status. This simple approach could be a good complementary alternative to the enrichment protocols used in the search for phosphoprotein candidates. While 482 spots were visualized with SYPRO-Ruby, 222 were with Pro-Q DPS. Statistically significant differences in spot intensity were observed among samples, these corresponding to 85 SYPRO-Ruby-, 20 Pro-Q-DPS-, and 35 SYPRO-Ruby and Pro-Q-DPS-stained spots. Fifty-five phosphoprotein candidates showing qualitative or quantitative differences between samples were subjected to MALDI-TOF-TOF MS analysis, with 20 of them being identified. Identified proteins belonged to five different functional categories, namely: carbohydrate and amino acid metabolism, defense, protein folding, and oxidation-reduction processes. With the exception of a putative cyclase, the other 19 proteins had at least one orthologous phosphoprotein in Arabidopsis thaliana, Medicago truncatula, N. tabacum, and Glycine max. Out of the 20 identified, seven showed differences in intensity in Pro-Q-DPS but not in SYPRO-Ruby-stained gels, including enzymes of the glycolysis and amino acid metabolism. This bears out that theory the regulation of these enzymes occurs at the post-translational level by phosphorylation with no changes at the transcriptional or translational level. This is different from the mechanism reported in orthodox seeds, in which concomitant changes in abundance and phosphorylation status have been observed for these enzymes.

Comparative proteomics of Euphorbia kansui Liou milky sap at two different developmental stages

Euphorbia kansui Liou is a unique traditional Chinese medicinal herb. Its milky sap proteins play important roles in laticifer development, synthesis and transport of its biologically active substances. A proteomic approach was applied to analyze the E. kansui latex proteins related to laticifer development and secondary metabolite synthesis by using sodium dodecyl sulfate-polyacrylamide gel electrophoresis and two-dimensional polyacrylamide gel electrophoresis. A total of 125 milky sap proteins associated with development of laticifers, disease and defense, and general metabolism were identified, and 19 differentially expressed proteins at two different developmental stages of laticifers were successfully detected. Peroxidase, cytochrome P450 mono-oxygenase superfamily, lipoxygenase, and multidrug resistance protein ABC transporter family may be involved in laticifer development, secondary metabolite synthesis and transport, and plant physiology.

Improving the quality of protein identification in non-model species. Characterization of Quercus ilex seed and Pinus radiata needle proteomes by using SEQUEST and custom databases

Nowadays the most used pipeline for protein identification consists in the comparison of the MS/MS spectra to reference databases. Search algorithms compare obtained spectra to an in silico digestion of a sequence database to find exact matches. In this context, the database has a paramount importance and will determine in a great deal the number of identifications and its quality, being this especially relevant for non-model plant species. Using a single Viridiplantae database (NCBI, UniProt) and TAIR is not the best choice for non-model species since they are underrepresented in databases resulting in poor identification rates. We demonstrate how it is possible to improve the rate and quality of identifications in two orphan species, Quercus ilex and Pinus radiata, by using SEQUEST and a combination of public (Viridiplantae NCBI, UniProt) and a custom-built specific database which contained 593,294 and 455,096 peptide sequences (Quercus and Pinus, respectively). These databases were built after gathering and processing (trimming, contiging, 6-frame translation) publicly available RNA sequences, mostly ESTs and NGS reads. A total of 149 and 1533 proteins were identified from Quercus seeds and Pinus needles, representing a 3.1- or 1.5-fold increase in the number of protein identifications and scores compared to the use of a single database. Since this approach greatly improves the identification rate, and is not significantly more complicated or time consuming than other approaches, we recommend its routine use when working with non-model species.

BIOLOGICAL SIGNIFICANCE

In this work we demonstrate how the construction of a custom database (DB) gathering all available RNA sequences and its use in combination with Viridiplantae public DBs (NCBI, UniProt) significantly improve protein identification when working with non-model species. Protein identification rate and quality is higher to those obtained in routine procedures based on using only one database (commonly Viridiplantae from NCBI), as we demonstrated analyzing Quercus seeds and Pine needles. The proposed approach based on the building of a custom database is not difficult or time consuming, so we recommend its routine use when working with non-model species. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Non-model organisms, a species endangered by proteogenomics

Previously, large-scale proteomics was possible only for organisms whose genomes were sequenced, meaning the most common model organisms. The use of next-generation sequencers is now changing the deal. With "proteogenomics", the use of experimental proteomics data to refine genome annotations, a higher integration of omics data is gaining ground. By extension, combining genomic and proteomic data is becoming routine in many research projects. "Proteogenomic"-flavored approaches are currently expanding, enabling the molecular studies of non-model organisms at an unprecedented depth. Today draft genomes can be obtained using next-generation sequencers in a rather straightforward way and at a reasonable cost for any organism. Unfinished genome sequences can be used to interpret tandem mass spectrometry proteomics data without the need for time-consuming genome annotation, and the use of RNA-seq to establish nucleotide sequences that are directly translated into protein sequences appears promising. There are, however, certain drawbacks that deserve further attention for RNA-seq to become more efficient. Here, we discuss the opportunities of working with non-model organisms, the proteomic methods that have been used until now, and the dramatic improvements proffered by proteogenomics. These put the distinction between model and non-model organisms in great danger, at least in terms of proteomics!

BIOLOGICAL SIGNIFICANCE

Model organisms have been crucial for in-depth analysis of cellular and molecular processes of life. Focusing the efforts of thousands of researchers on the Escherichia coli bacterium, Saccharomyces cerevisiae yeast, Arabidopsis thaliana plant, Danio rerio fish and other models for which genetic manipulation was possible was certainly worthwhile in terms of fundamental and invaluable biological insights. Until recently, proteomics of non-model organisms was limited to tedious, homology-based techniques, but today draft genomes or RNA-seq data can be straightforwardly obtained using next-generation sequencers, allowing the establishment of a draft protein database for any organism. Thus, proteogenomics opens new perspectives for molecular studies of non-model organisms, although they are still difficult experimental organisms. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Physiological and proteomic analyses of drought stress response in Holm oak provenances

Responses to drought stress by water withholding have been studied in 1 year old Holm oak (Quercus ilex subsp. ballota [Desf.] Samp.) seedlings from seven provenances from Andalusia (southern Spain). Several physiological parameters, including predawn xylem water potentials and relative water content in soil, roots, and leaves as well as maximum quantum efficiency and yield of PSII were evaluated for 28 days in both irrigated and nonirrigated seedlings. The leaf proteome map of the two provenances that show the extreme responses (Seville, GSE, is the most susceptible, while Almerı́a, SSA, is the least susceptible) was obtained. Statistically significant variable spots among provenances and treatments were subjected to MALDI-TOF/TOF-MS/MS analysis for protein identification. In response to drought stress, ~12.4% of the reproducible spots varied significantly depending on the treatment and the population. These variable proteins were mainly chloroplastic and belonged to the metabolism and defense/stress functional categories. The 2-DE protein profile of nonirrigated seedlings was similar in both provenances. Physiological and proteomics data were generally in good agreement. The general trend was a decrease in protein abundance upon water withholding in both provenances, mainly in those involved in ATP synthesis and photosynthesis. This decrease, moreover, was most marked in the most susceptible population compared with the less susceptible one.