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.

Comparative proteomics of Pinus-Fusarium circinatum interactions reveal metabolic clues to biotic stress resistance.

Fusarium circinatum, causing pine pitch canker (PPC), affects conifers productivity and health worldwide. Selection and breeding for resistance arises as the most promising approach to fight PPC. Therefore, it is crucial to explore the response of hosts with varying levels of susceptibility to PPC to unveil the genes/pathways behind these phenotypes. We evaluated the dynamics of the needle proteome of a susceptible (Pinus radiata) and a relatively resistant (Pinus pinea) species upon F. circinatum inoculation by GeLC-MS/MS. Integration with physiological data and validation of key genes by qPCR allowed to identify core pathways regulating these contrasting responses. In P. radiata, the pathogen may target both the secondary metabolism to negatively regulate immune response and chloroplast redox proteins to increase energy-producing pathways for amino acid production in its favour. In contrast, chloroplast redox regulation may assure redox homeostasis in P. pinea, as well as nonenzymatic antioxidants. The presence of membrane trafficking-related proteins exclusively in P. pinea likely explains its defence response against F. circinatum. A crosstalk between abscisic acid and epigenetic regulation of gene expression is also proposed in PPC response. These results are useful to support breeding programs aiming to achieve PPC resistance.

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.

Quantitative Analysis of Target Peptides Related to Resistance Against Ascochyta Blight (Peyronellaea pinodes) in Pea.

Peyronellaea pinodes causes Ascochyta blight, one of the major diseases in pea worldwide. Cultivated pea plants have a low resistance to this disease. Although quantitative trait loci (QTLs) involved in the resistance to Ascochyta blight have been identified, the specific genes associated with these QTLs remain unknown, which makes marker-assisted selection difficult. Complex traits alter proteins and their abundance. Quantitative estimation of proteins in pea might therefore be useful in selecting potential markers for breeding. In this work, we developed a strategy using a combination of shotgun proteomics (viz., high performance liquid chromatography-mass spectrometry data-dependent acquisition) and data-independent acquisition (DIA) analysis, to identify putative protein markers associated with resistance to Ascochyta blight and explored its use for breeding selection. For this purpose, an initial list of target peptides based on proteins closely related to resistance to P. pinodes was compiled by using two genotypes with contrasting responses to the disease. Then, targeted data analysis (viz., shotgun proteomics-DIA) was used for constitutive quantification of the target peptides in a representative number of the recombinant inbred line population segregated for resistance as derived from a cross between the two genotypes. Finally, a peptide panel of potential markers for resistance to P. pinodes was built. The results thus obtained are discussed and compared with those of previous gene expression studies using the same parental pea genotypes responding to the pathogen. Also, a molecular defense mechanism against Ascochyta blight in pea is proposed. To the authors' knowledge, this is the first time a targeted proteomics approach based on data analysis has been used to identify peptides associated with resistance to this disease.

Quercus ilex pollen allergen, Que i 1, responsible for pollen food allergy syndrome caused by fruits in Spanish allergic patients.

BACKGROUND: Pollen food allergy syndrome (PFAS) related to PR10 from vegetables is common in northern Europe, whereas in Mediterranean countries PFAS has been preferentially associated with profilins. However, there are pollen-allergic patients reactive to Bet v 1 in birch-free regions. Since it cannot be the primary sensitizer, there has to be another culprit. Quercus ilex is a good candidate as it belongs to the order Fagales. This order includes trees with highly sensitizing pollen such as alder, hazel, hornbeam, oak and chestnut because of the presence of PR10 allergens. PR10 allergens have indeed been described in other Quercus species. OBJECTIVE: Our goals were to determine the rate of sensitization to Q. ilex in central Spain and the associated frequency of PFAS; secondly to identify and clone the Q. ilex allergen PR10. METHODS: We included 224 allergic patients with respiratory symptoms to estimate the rate of sensitization. A skin prick test (SPT) and ImmunoCAP were performed. A total of 38 Q. ilex-sensitized patients were tested using Western blotting to determine the rate of Que i 1. Peptides from Que i 1 were analysed by MALDI-TOF/TOF and Orbitrap LC-MSMS. The Que i 1 sequence was first obtained from the Holm oak transcriptome then cloned and expressed in bacteria. RESULTS: 59.8% of pollen-allergic children were sensitized to Q. ilex. We described and cloned the Q. ilex PR10, Que i 1, which has a sensitization rate of 60.5% and was recognized by 65.4% patients reporting PFAS. CONCLUSION AND CLINICAL RELEVANCE: Sensitization to Q. ilex pollen has increased significantly since 1995. This sensitization could be important, as the presence of PFAS in this population is higher than in patients not sensitized to Q. ilex. The first Q. ilex allergen has been described and is related to PFAS in Spanish patients sensitized to PR10 but not exposed to birch pollen.

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.

Drought and Recovery: Independently Regulated Processes Highlighting the Importance of Protein Turnover Dynamics and Translational Regulation in Medicago truncatula.

Climate change in conjunction with population growth necessitates a systems biology approach to characterize plant drought acclimation as well as a more thorough understanding of the molecular mechanisms of stress recovery. Plants are exposed to a continuously changing environment. Extremes such as several weeks of drought are followed by rain. This requires a molecular plasticity of the plant enabling drought acclimation and the necessity of deacclimation processes for recovery and continuous growth.During drought stress and subsequent recovery, the metabolome and proteome are regulated through a sequence of molecular processes including synthesis and degradation and molecular interaction networks are part of this regulatory process. In order to study this complex regulatory network, a comprehensive analysis is presented for the first time, investigating protein turnover and regulatory classes of proteins and metabolites during a stress recovery scenario in the model legume Medicago truncatula The data give novel insights into the molecular capacity and differential processes required for acclimation and deacclimation of severe drought stressed plants.Functional cluster and network analyses unraveled independent regulatory mechanisms for stress and recovery with different dynamic phases that during the course of recovery define the plants deacclimation from stress. The combination of relative abundance levels and turnover analysis revealed an early transition phase that seems key for recovery initiation through water resupply and is independent from renutrition. Thus, a first indication for a metabolite and protein-based load capacity was observed necessary for the recovery from drought, an important but thus far ignored possible feature toward tolerance. The data indicate that apart from the plants molecular stress response mechanisms, plasticity may be related to the nutritional status of the plant prior to stress initiation. A new perspective and possible new targets as well as metabolic mechanisms for future plant-bioengineering toward enhanced drought stress tolerance are presented.

Label-free quantitative proteomic analysis of tolerance to drought in Pisum sativum.

Abiotic stresses caused by adverse environmental conditions are responsible for heavy economic losses on pea crop, being drought one of the most important abiotic constraints. Development of pea cultivars well adapted to dry conditions has been one of the major tasks in breeding programs. The increasing food requirements drive the necessity to broaden the molecular basis of tolerance to drought to develop pea cultivars well adapted to dry conditions. We have used a shotgun proteomic approach (nLC-MSMS) to study the tolerance to drought in three pea genotypes that were selected based on differences in the level of water deficit tolerance. Multivariate statistical analysis of data unraveled 367 significant differences of 700 identified when genotypes and/or treatment were compared. More than half of the significantly changed proteins belong to primary metabolism and protein regulation categories. We propose different mechanisms to cope drought in the genotypes studied. Maintenance of the primary metabolism and protein protection seems a strategy for drought tolerance. On the other hand susceptibility might be related to maintenance of the homeostatic equilibrium, a very energy consuming process. Data are available via ProteomeXchange with identifier PXD004587.

Heat-Treatment-Responsive Proteins in Different Developmental Stages of Tomato Pollen Detected by Targeted Mass Accuracy Precursor Alignment (tMAPA).

Recently, we have developed a quantitative shotgun proteomics strategy called mass accuracy precursor alignment (MAPA). The MAPA algorithm uses high mass accuracy to bin mass-to-charge (m/z) ratios of precursor ions from LC-MS analyses, determines their intensities, and extracts a quantitative sample versus m/z ratio data alignment matrix from a multitude of samples. Here, we introduce a novel feature of this algorithm that allows the extraction and alignment of proteotypic peptide precursor ions or any other target peptide from complex shotgun proteomics data for accurate quantification of unique proteins. This strategy circumvents the problem of confusing the quantification of proteins due to indistinguishable protein isoforms by a typical shotgun proteomics approach. We applied this strategy to a comparison of control and heat-treated tomato pollen grains at two developmental stages, post-meiotic and mature. Pollen is a temperature-sensitive tissue involved in the reproductive cycle of plants and plays a major role in fruit setting and yield. By LC-MS-based shotgun proteomics, we identified more than 2000 proteins in total for all different tissues. By applying the targeted MAPA data-processing strategy, 51 unique proteins were identified as heat-treatment-responsive protein candidates. The potential function of the identified candidates in a specific developmental stage is discussed.

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.

A first draft genome of holm oak (Quercus ilex subsp. ballota), the most representative species of the Mediterranean forest and the Spanish agrosylvopastoral ecosystem "dehesa".

The holm oak (Quercus ilex subsp. ballota) is the most representative species of the Mediterranean Basin and the agrosylvopastoral Spanish "dehesa" ecosystem. Being part of our life, culture, and subsistence since ancient times, it has significant environmental and economic importance. More recently, there has been a renewed interest in using the Q. ilex acorn as a functional food due to its nutritional and nutraceutical properties. However, the holm oak and its related ecosystems are threatened by different factors, with oak decline syndrome and climate change being the most worrying in the short and medium term. Breeding programs informed by the selection of elite genotypes seem to be the most plausible biotechnological solution to rescue populations under threat. To achieve this and other downstream analyses, we need a high-quality and well-annotated Q. ilex reference genome. Here, we introduce the first draft genome assembly of Q. ilex using long-read sequencing (PacBio). The assembled nuclear haploid genome had 530 contigs totaling 842.2 Mbp (N50 = 3.3 Mbp), of which 448.7 Mb (53%) were repetitive sequences. We annotated 39,443 protein-coding genes of which 94.80% were complete and single-copy genes. Phylogenetic analyses showed no evidence of a recent whole-genome duplication, and high synteny of the 12 chromosomes between Q. ilex and Quercus lobata and between Q. ilex and Quercus robur. The chloroplast genome size was 142.3 Kbp with 149 protein-coding genes successfully annotated. This first draft should allow for the validation of omics data as well as the identification and functional annotation of genes related to phenotypes of interest such as those associated with resilience against oak decline syndrome and climate change and higher acorn productivity and nutraceutical value.

Proteomics analysis of date palm leaves affected at three characteristic stages of brittle leaf disease.

Proteomics analysis has been performed in leaf tissue from field date palm trees showing the brittle leaf disease (BLD) or maladie des feuilles cassantes, the main causal agent of the date palm decline in south Tunisia. To study the evolution of the disease, proteins from healthy and affected leaves taken at three disease stages (S1, S2 and S3) were trichloroacetic acid acetone extracted and subjected to two-dimensional gel electrophoresis (5-8 pH range). Statistical analysis showed that the protein abundance profile is different enough to differentiate the affected leaves from the healthy ones. Fifty-eight variable spots were successfully identified by matrix-assisted laser desorption ionization time of flight, 60 % of which corresponded to chloroplastic ones being involved in the photosynthesis electronic chain and ATP synthesis, metabolic pathways implicated in the balance of the energy, and proteases. Changes in the proteome start at early disease stage (S1), and are greatest at S2. In addition to the degradation of the ribulose-1.5-bisphosphate carboxylase oxygenase in affected leaflets, proteins belonging to the photosynthesis electronic chain and ATP synthesis decreased following the disease, reinforcing the relationship between BLD and manganese deficiency. The manganese-stabilizing proteins 33 kDa, identified in the present work, can be considered as protein biomarkers of the disease, especially at early disease step.

Multiomic Data Integration in the Analysis of Drought-Responsive Mechanisms in Quercus ilex Seedlings.

The integrated analysis of different omic layers can provide new knowledge not provided by their individual analysis. This approach is also necessary to validate data and reveal post-transcriptional and post-translational mechanisms of gene expression regulation. In this work, we validated the possibility of applying this approach to non-model species such as Quercus ilex. Transcriptomics, proteomics, and metabolomics from Q. ilex seedlings subjected to drought-like conditions under the typical summer conditions in southern Spain were integrated using a non-targeted approach. Two integrative approaches, PCA and DIABLO, were used and compared. Both approaches seek to reduce dimensionality, preserving the maximum information. DIABLO also allows one to infer interconnections between the different omic layers. For easy visualization and analysis, these interconnections were analyzed using functional and statistical networks. We were able to validate results obtained by analyzing the omic layers separately. We identified the importance of protein homeostasis with numerous protease and chaperones in the networks. We also discovered new key processes, such as transcriptional control, and identified the key function of transcription factors, such as DREB2A, WRKY65, and CONSTANS, in the early response to drought.

Identification of Proteases and Protease Inhibitors in Seeds of the Recalcitrant Forest Tree Species Quercus ilex.

Proteases and protease inhibitors have been identified in the recalcitrant species Quercus ilex using in silico and wet methods, with focus on those present in seeds during germination. In silico analyses showed that the Q. ilex transcriptome database contained 2,240 and 97 transcripts annotated as proteases and protease inhibitors, respectively. They belonged to the different families according to MEROPS, being the serine and metallo ones the most represented. The data were compared with those previously reported for other Quercus species, including Q. suber, Q. lobata, and Q. robur. Changes in proteases and protease inhibitors alongside seed germination in cotyledon and embryo axis tissues were assessed using proteomics and in vitro and in gel activity assays. Shotgun (LC-MSMS) analysis of embryo axes and cotyledons in nonviable (NV), mature (T1) and germinated (T3) seeds allowed the identification of 177 proteases and 12 protease inhibitors, mostly represented by serine and metallo types. Total protease activity, as determined by in vitro assays using azocasein as substrate, was higher in cotyledons than in embryo axes. There were not differences in activity among cotyledon samples, while embryo axis peaked at germinated T4 stage. Gel assays revealed the presence of protease activities in at least 10 resolved bands, in the Mr range of 60-260 kDa, being some of them common to cotyledons and embryo axes in either nonviable, mature, and germinated seeds. Bands showing quantitative or qualitative changes upon germination were observed in embryo axes but not in cotyledons at Mr values of 60-140 kDa. Proteomics shotgun analysis of the 10 bands with protease activity supported the results obtained in the overall proteome analysis, with 227 proteases and 3 protease inhibitors identified mostly represented by the serine, cysteine, and metallo families. The combined use of shotgun proteomics and protease activity measurements allowed the identification of tissue-specific (e.g., cysteine protease inhibitors in embryo axes of mature acorns) and stage-specific proteins (e.g., those associated with mobilization of storage proteins accumulated in T3 stage). Those proteins showing differences between nonviable and viable seeds could be related to viability, and those variables between mature and germinated could be associated with the germination process. These differences are observed mostly in embryo axes but not in cotyledons. Among them, those implicated in mobilization of reserve proteins, such as the cathepsin H cysteine protease and Clp proteases, and also the large number of subunits of the CNS and 26S proteasome complex differentially identified in embryos of the several stages suggests that protein degradation via CNS/26S plays a major role early in germination. Conversely, aspartic proteases such as nepenthesins were exclusively identified in NV seeds, so their presence could be used as indicator of nonviability.

Proteomic Changes in Paspalum fasciculatum Leaves Exposed to Cd Stress.

(1) Background: Cadmium is a toxic heavy metal that is widely distributed in water, soil, and air. It is present in agrochemicals, wastewater, battery waste, and volcanic eruptions. Thus, it can be absorbed by plants and enter the trophic chain. P. fasciculatum is a plant with phytoremediation capacity that can tolerate Cd stress, but changes in its proteome related to this tolerance have not yet been identified. (2) Methods: We conducted a quantitative analysis of the proteins present in P. fasciculatum leaves cultivated under greenhouse conditions in mining soils doped with 0 mg kg-1 (control), 30 mg kg-1, or 50 mg kg-1. This was carried out using the label-free shotgun proteomics technique. In this way, we determined the changes in the proteomes of the leaves of these plants, which allowed us to propose some tolerance mechanisms involved in the response to Cd stress. (3) Results: In total, 329 variable proteins were identified between treatments, which were classified into those associated with carbohydrate and energy metabolism; photosynthesis; structure, transport, and metabolism of proteins; antioxidant stress and defense; RNA and DNA processing; and signal transduction. (4) Conclusions: Based on changes in the differences in the leaf protein profiles between treatments, we hypothesize that some proteins associated with signal transduction (Ras-related protein RABA1e), HSPs (heat shock cognate 70 kDa protein 2), growth (actin-7), and cellular development (actin-1) are part of the tolerance response to Cd stress.

Application and optimization of label-free shotgun approaches in the study of Quercus ilex.

Advances in proteomic equipment, algorithms and wet protocols are being increasingly reported. Each step in the experimental workflow must be adapted and optimized to the target experimental system and objectives. The influence of the amount of peptides loaded onto the column in shotgun platforms has rarely been considered to date even though it dictates the confidence with which proteins can be identified and quantified. An experiment using variable dilutions of protein equivalent mixtures of root, leaf and seed tissue extracts of Quercus ilex was performed by subjecting BSA protein equivalent amounts of 1-100 µg to SDS-PAGE, the resulting bands being trypsin digested and peptides (10-1000 ng protein equivalents) loaded onto an LC column. Mass spectra were used to identify proteins against the in-house Q. ilex transcriptome database. Determinations included SEQUEST quantification (average of the three most abundant distinct peptides for each protein) and proteotypic peptides. The number of proteins identified was found to depend on peptide load and to peak at 2054 with 600 ng. Smaller loads led to linearly decreasing identifications from 1859 with 400 ng to 495 with 10 ng. Both quantification strategies provided similar results. The linear dynamic range was from 100 to 600 ng.

GeLC-Orbitrap/MS and 2-DE-MALDI-TOF/TOF comparative proteomics analysis of seed cotyledons from the non-orthodox Quercus ilex tree species.

Gel electrophoresis-based and shotgun approaches are the most employed proteomic platforms in plant biology research, with the latter replacing the former in the last years. We have compared 2-DE-MALDI-TOF/TOF and GeLC-Orbitrap/MS analyses using the same protein extracts from Quercus ilex cotyledons at different development stages. The results obtained (ProteomeXchange available data, PXD020603) showed that both platforms were complementary, showing common and specific proteins identified in each case, but leading to similar biological conclusions. Protein analysis identified 562 spots in gel-based (292 variables) and 2409 proteins in shotgun (560 variables), that were detected with both platforms and represent common key pathways related to maturation and germination. The main differences concern hormone metabolism, storage and late embryogenesis abundant proteins. Deeper proteome coverage was obtained with the shotgun approach, with a greater number of metabolic pathways represented, as gibberellin biosynthesis, not observed in the gel-based analysis. Nevertheless, several storage proteins, highly abundant in cotyledons and well represented in gel-based platform were not identified using the shotgun platform. These results support that when analyzing any plant biological process, the use of both platforms is complementary rather than redundant, that favors an in-depth proteomic analysis and a more confident biological interpretation of the data obtained.

Changes in the transcript and protein profiles of Quercus ilex seedlings in response to drought stress.

Quercus ilex is the dominant tree species in natural forest ecosystems across the Mediterranean Basin and in the agrosilvopastoral system dehesa, which has a high ecological and economical significance. As in other forestry species, survival in Q. ilex is threatened by long periods of drought. This paper reports the transcriptome and proteome profiles of 6-month-old seedlings subjected to severe drought conditions. Drought was imposed by water withholding in seedlings grown in perlite for 28 days. Seedling leaves were collected when leaf fluorescence had decreased by 20% and 45% relative to well-watered seedlings. The transcriptome and proteome were analyzed by using Illumina and shotgun platforms. The quality and confidence of the mRNA and protein identifications and quantifications were assessed, obtaining 25,169 transcripts and 3312 proteins. Variable transcripts and proteins were analyzed by Venn diagram, Pearson's correlation, GO enrichment, KEGG pathways, multivariate analysis and interaction networks. Despite the poor correlation between mRNA and protein, both platforms gave a complementary view of the changes in the abundance of several gene products under drought conditions and indicated that gene expression regulation and translation to phenotype is quite complex and gene-specific. As a general tendency, while transcripts and proteins of the metabolism were down-accumulated, those of stress related were up-accumulated. Out of the variable dataset, four gene products (viz., FtSH6, CLPB1, CLPB3, and HSP22) were up-accumulated at both omics levels at the two surveyed times, being the first work where they are described in drought response in forest species. These chaperones and proteases could be considered as potential drought tolerance markers to be used in the selection of elite, resilient genotypes, and in breeding programs.

Combining P and Zn fertilization to enhance yield and grain quality in maize grown on Mediterranean soils.

The main aim of this study was to elucidate the effect of individual and joint fertilization with P and Zn on maize plants grown on typical Mediterranean soils with a limited Zn availability. For this purpose, we examined the effects of P and Zn fertilization individually and in combination on growth, yield and grain protein content in maize grown in pots filled with three different Mediterranean soils (LCV, FER and INM). Phosphorus and Zn translocation to grain was impaired, and aboveground dry matter and yield at harvest reduced by 8-85% (LCV and FER), in plants treated with Zn or P alone relative to unfertilized (control) plants. In contrast, joint fertilization with P and Zn enhanced translocation of these nutrients to grain and significantly increased aboveground dry matter (30% in LCV, 50% in FER and 250% in INM) and grain Zn availability in comparison with control plants. Also, joint application of both nutrients significantly increased grain P (LCV) and Zn (LCV and FER) use efficiency relative P and Zn, respectively, alone. Yield was increased between 31% in LCV and 121% in FER relative to control plants, albeit not significantly. Fertilization with P or Zn significantly influenced the abundance of specific proteins affecting grain quality (viz., storage, lys-rich and cell wall proteins), which were more abundant in mature grains from plants fertilized with Zn alone and, to a lesser extent, P + Zn. Sustainable strategies in agriculture should consider P-Zn interactions in maize grown on soils with a limited availability of Zn, where Zn fertilization is crucial to ensure grain quality.

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].