{'fraction_non_alphanumeric': 0.03420878975847961, 'fraction_numerical': 0.01808103471030751, 'mean_word_length': 5.366324987397076, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 1, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Insights from a Multi-Omics Integration (MOI) Study in Oil Palm (Elaeis guineensis Jacq.) Response to Abiotic Stresses: Part Two—Drought

Drought and salinity are two of the most severe abiotic stresses affecting agriculture worldwide and bear some similarities regarding the responses of plants to them. The first is also known as osmotic stress and shows similarities mainly with the osmotic effect, the first phase of salinity stress. Multi-Omics Integration (MOI) offers a new opportunity for the non-trivial challenge of unraveling the mechanisms behind multigenic traits, such as drought and salinity resistance. The current study carried out a comprehensive, large-scale, single-omics analysis (SOA) and MOI studies on the leaves of young oil palm plants submitted to water deprivation. After performing SOA, 1955 DE enzymes from transcriptomics analysis, 131 DE enzymes from proteomics analysis, and 269 DE metabolites underwent MOI analysis, revealing several pathways affected by this stress, with at least one DE molecule in all three omics platforms used. Moreover, the similarities and dissimilarities in the molecular response of those plants to those two abiotic stresses underwent mapping. Cysteine and methionine metabolism (map00270) was the most affected pathway in all scenarios evaluated. The correlation analysis revealed that 91.55% of those enzymes expressed under both stresses had similar qualitative profiles, corroborating the already known fact that plant responses to drought and salinity show several similarities. At last, the results shed light on some candidate genes for engineering crop species resilient to both abiotic stresses.


Introduction
The palm oil industry established itself in places around the world where the oil palm (Elaeis guineensis Jacq.), the number one source of vegetable oil consumed today, exists in its natural environment of high rainfall throughout the year. Those places are in an area of 5 • latitude from the Equator line, and as one moves further north or south in the tropical belt, one experiences long periods of drought throughout the year when oil palm does not receive the physiological water level required to maintain productivity. Consequently, the oil palm grower must embrace artificial irrigation, together with the production cost linked to it, not to mention the risk of facing problems of soil salinization.
A recent report by the UN Convention to Combat Desertification, entitled Drought in Numbers 2022-Restoration to Readiness and Resilience, brings to light the dimension of the damages caused by this abiotic stress on the planet and shows to some extent what is needed to overcome this problem and sustainably ensure food security [1]. Drought stress is one of the most severe environmental stresses affecting agriculture worldwide. The efficient and effective use of genetics and biotechnology is among the most proactive interventions necessary to mitigate the effects of water deprivation on plant productivity [2]. Plant response to drought is a complex biological system that needs investigation via an integrative approach that must involve both drought-tolerant and intolerant genotypes. Resistance to drought in plants results from the dynamic interaction of various parts of this system that need deep understanding in order to enable the design of approaches for the effective vertical or horizontal transfer of such traits to drought intolerant genotypes [2,3]. A Multi-Omics Integration (MOI) study [4][5][6], made viable due to the recent massive generation of genomics, transcriptomics, proteomics, and metabolomics datasets, is helping the non-trivial challenge of unraveling the mechanisms behind this multigenic trait, and mapping not only resistance genes but also their hot spots for precision genome editing [7].
Plant responses to drought and salinity stresses show several similarities, especially in the first phase of the latter, the osmotic effect [8][9][10]. Distinct studies showed the existence of genes able to promote resistance to more than one abiotic stress at once. For instance, Brini et al. [11] showed that the transgenic Arabidopsis thaliana plants overexpressing two distinct genes from wheat were much more resistant to both stresses than the isogenic wild-type plants. The same was true for Feng et al. [12], who overexpressed a gene from pepper in A. thaliana and showed that this gene positively regulates heat, salt, and drought resistance in pepper. In another example, Zhang et al. [13] overexpressed three genes from sweet potato in sweet potato and obtained genetically modified plants tolerant to salt and drought stresses.
Jha et al. [14], Zargar et al. [15], Muthuramalingam et al. [16], and Bittencourt et al. [17] are a few examples of recent reports dealing with the use of MOI approaches toward understanding plant responses to these abiotic stresses. The current study is a step forward in our research activities on characterizing the biochemical and molecular response of oil palm to abiotic stress. After carrying out a comprehensive, large-scale, single-omics analysis (SOA) and an MOI analysis of the metabolome, transcriptome, and proteome profiles on the leaves of young oil palm plants submitted to salinity stress [17], we here report a similar study characterizing the response of this oilseed crop to water deprivation. Furthermore, the differentially expressed (DE) molecules responsive to both drought (from this study) and salt (from [17]) underwent additional MOI analysis and correlation study to identify the commonalities and differences in the response of such plants to both stresses.
This group of 8421 DE protein sequences was submitted to analysis in the BlastKOALA annotation tool for the K number assignment, revealing 4133 DE proteins with K number and 1955 with Enzyme Commission (EC) numbers. This set of 1955 DE enzymes underwent gene ontology analyses to classify them accordingly to biological process (BP), molecular function (MF), and cellular component (CC).
The BP subgroups with more representatives were protein phosphorylation, followed by carbohydrate metabolic process, phosphorylation, proteolysis, and cellular oxidant detoxification. For MF, the most populated subgroups were ATP binding, metal ion binding, iron ion binding, and heme binding. Finally, for the CC subgroup, the integral component of the membrane came in first, followed by the cytoplasm, nucleus, and cytosol components ( Figure 1a). Furthermore, the DE enzymes underwent classification according to the EC number, a numerical classification scheme for enzymes based on the chemical reaction. The first classification level involves a general type of enzyme-catalyzed reaction that ranges from one to six, and the oxidoreductases (EC 1), transferases (EC 2), and hydrolases (EC 3) classes dominated. Enzymes acting on paired donors, with the incorporation or reduction in molecular oxygen (EC 1.14) was the most populated group among the subclasses of oxidoreductases, with 86 hits. The most representative subclass of transferases included those with a function of transferring phosphorus-containing groups, glycosyltransferases, and acyltransferases, with 218, 123, and 86 hits, respectively. Finally, the hydrolases had the following subclasses standing out with over 50 hits: compounds involved and acting on ester bonds, acting on peptide bonds (peptidases), glycosylases, and acting on acid anhydrides (Figure 1b). The BP subgroups with more representatives were protein phosphorylation, followed by carbohydrate metabolic process, phosphorylation, proteolysis, and cellular oxidant detoxification. For MF, the most populated subgroups were ATP binding, metal ion binding, iron ion binding, and heme binding. Finally, for the CC subgroup, the integral component of the membrane came in first, followed by the cytoplasm, nucleus, and cytosol components (Figure 1a). Furthermore, the DE enzymes underwent classification according to the EC number, a numerical classification scheme for enzymes based on the chemical reaction. The first classification level involves a general type of enzyme-catalyzed reaction that ranges from one to six, and the oxidoreductases (EC 1), transferases (EC 2), and hydrolases (EC 3) classes dominated. Enzymes acting on paired donors, with the incorporation or reduction in molecular oxygen (EC 1.14) was the most populated group among the subclasses of oxidoreductases, with 86 hits. The most representative subclass of transferases included those with a function of transferring phosphorus-containing groups, glycosyltransferases,

Oil Palm Proteome under Drought Stress
It was possible to identify with high confidence (FDR ≤ 0.01) a total of 3659 and 3824 peptides in control and stressed samples, respectively, which infers up to 1859 protein entries from E. guineensis proteome (Uniprot) in both conditions (Table 1).
Seven hundred and seventy-nine proteins were inferred from more than four peptides and five hundred and ten had at least one proteotypic peptide observation. All peptides and proteins identified and the list of 1085 proteins according to the maximum parsimony  Tables S5 and S6). Fifty-six proteins showed statistically significant differences in their abundance between stressed and control samples; among them, there were thirty-five proteins upregulated and twenty-one downregulated (Supplementary Table S7). Our differential abundance analysis considered proteins identified at least in two replicates in each condition. A group of 284 proteins that attended to the statistical criteria of PatternLab V software underwent functional annotation and MOI analyses, being 76 found exclusively in the control plants, 152 found only in the drought-stressed plants, and 56 showed statistically significant differences in their abundance between the two treatments used. Those 284 DE proteins underwent analysis in the BlastKOALA annotation tool, revealing 215 DE proteins with K numbers and 131 with EC numbers.
This set of 131 DE enzymes underwent gene ontology analyses to classify them accordingly to biological process (BP), molecular function (MF), and cellular component (CC). The BP top subgroups were carbohydrate metabolic processes, followed by cellular oxidant detoxification and glycolytic process, with 13, 10, and 10 hits, respectively. For MF, the proteins were mainly from subgroups of ATP binding, metal ion binding, magnesium ion binding, and pyridoxal phosphate binding, all with over ten hits. Finally, the CC of the cytoplasm came in first, followed by the cytosol, integral component of membrane, and mitochondrion, all with over ten hits (Figure 2a).
The prevalent chemical reactions according to EC were Oxireductases (EC 1), transferases (EC 2), and hydrolases (EC 3) classes. In the subclasses of oxidoreductases, the top group was acting on a peroxide as an acceptor, followed by acting on the CH-OH group of donors and acting on the aldehyde or oxo group of donors. The most representative subclass of transferases class was the one with transferring phosphorus-containing groups ( Figure 2b).
The results from transcriptome and proteome single-omics analysis (SOA), when comparing drought-stressed young oil palm plants against the control ones, show that the ATP binding came as the number one BP affected by the stress, followed by the metal ion binding process. Meanwhile, transferring phosphorus-containing groups, glycosyltransferases, and acyltransferases were subclasses of prevalent chemical reactions in both SOA (Figures 1 and 2).

Oil Palm Metabolome under Drought Stress
The statistical analysis performed on the MetaboAnalyst has returned 3363 and 2538 peaks in the positive and negative polar fractions, respectively. Ninety-six were up-and seventy-one downregulated in the polar-positive. Two hundred ninety-one were up-and nine downregulated, in the polar-negative.
The results from transcriptome and proteome single-omics analysis (SOA), when comparing drought-stressed young oil palm plants against the control ones, show that the ATP binding came as the number one BP affected by the stress, followed by the metal ion binding process. Meanwhile, transferring phosphorus-containing groups, glycosyltransferases, and acyltransferases were subclasses of prevalent chemical reactions in both SOA (Figures 1 and 2).

Oil Palm Metabolome under Drought Stress
The statistical analysis performed on the MetaboAnalyst has returned 3363 and 2538 peaks in the positive and negative polar fractions, respectively. Ninety-six were up-and

Integrating Oil Palm Transcriptome, Proteome and Metabolome under Drought Stress
The approach used to integrate the datasets from the three omics platforms was a pathway-based mapping approach, similar to that previously used by Bittencourt et al. [17] ( Figure 3).
All enzymes (from SOA analysis-transcriptomics or proteomics) and metabolites (from metabolomics SOA analysis), selected as differentially expressed in the leaves of young oil palm plants (stressed/control), underwent MOI analysis.
When applying the Omics Fusion platform to perform the MOI analysis, a group of 56 pathways appeared as affected by drought stress, with at least one molecule differentially expressed in each of the three omics platforms used. Those pathways with ≥20 unique molecules differentially expressed are shown in Table 2. The Cysteine and

Commonalities and Dissimilarity in the Transcriptome, Proteome, and Metabolome Profiles of Young Oil Palm Plants Separately Submitted to Drought and Salinity Stresses
To generate the data regarding the commonalities and dissimilarity in the response of young oil palm plants separately submitted to drought and salinity stresses, all datasets generated in this present study underwent a comparison analysis with the equivalent datasets that had been previously generated and used in Bittencourt et al. [17].
In the transcriptomics analysis, there is a positive correlation between the behavior of the 554 enzymes differentially expressed under drought and salinity stresses (Supplementary Table S9) when 91.15% of the enzymes had the same qualitative profile (upregulated or downregulated), independently of the stress condition ( Figure 5). This result is a strong amount of evidence that there are several commonalities regarding gene expressionmolecular dynamics-in young oil palm plants submitted separately to these two abiotic stresses. It is also a valuable source of information to further study the osmotic effect. However, 8 Table S8). The tyrosine metabolism (map0 phenylpropanoid biosynthesis (map00940); arginine biosynthesis (map0 isoquinoline alkaloid biosynthesis (map00950); pyrimidine metabolism (map0 arginine and proline metabolism (map00330); betalain biosynthesis (map0 monobactam biosynthesis (map00261); valine, leucine and isoleucine biosyn (map00290); pentose phosphate pathway (map00030); aminoacyl-tRNA biosyn (map00970); flavonoid biosynthesis (map00941); histidine metabolism (map0 stilbenoid, diarylheptanoid and gingerol biosynthesis (map00945); riboflavin metab (map00740); phenylalanine, tyrosine and tryptophan biosynthesis (map00400); an carbon pool by folate (map00670) pathways were the most affected pathways with p ≤ 0.05; in this order, from top to bottom.

Integrating Oil Palm Transcriptome, Proteome and Metabolome under Drought Stress
The approach used to integrate the datasets from the three omics platforms pathway-based mapping approach, similar to that previously used by Bittencour [17] (Figure 3). All enzymes (from SOA analysis-transcriptomics or proteomics) and metab In the case of the proteomics analysis, a group of enzymes were also differentially expressed under both abiotic stresses ( Table 3). Most of them had the same qualitative profile independently of the stress condition and a minority had different ones. For instance, enzymes A0A6J0PP21 and A0A6I9S451, two isoforms coded by the same gene (LOC105053549) in the oil palm genome and described as L-ascorbate peroxidases in the KEGG database, were upregulated under drought stress and downregulated almost to no detection under salt stress.
In the case of the metabolomics analysis, a group of 13 metabolites was also differentially expressed in both scenarios. Among them, eight had the same qualitative profile, independently of the stress and five had different ones (Figure 7). The metabolite Melibiitol (C05399) was the one that showed the highest quantitative difference in expression, with a Log 2 (FC) of −4.698 (or 3.85% of the initial expression level in the control treatment) under salinity stress and of 1.53 (or 190% on the top of the initial expression level) under drought stress.
At last, when applying the Omics Fusion platform to perform the MOI analysis using only the proteins and metabolites DE in both scenarios, Cysteine and methionine metabolism came first in the list of most affected pathways, with 15 full-length transcripts coding for enzymes, two enzymes (from the proteomics analysis), and two metabolites ( Table 4).

Commonalities and Dissimilarity in the Transcriptome, Proteome, and Metabolome Profiles of Young Oil Palm Plants Separately Submitted to Drought and Salinity Stresses
To generate the data regarding the commonalities and dissimilarity in the response of young oil palm plants separately submitted to drought and salinity stresses, all datasets result is a strong amount of evidence that there are several commonalities regarding gene expression-molecular dynamics-in young oil palm plants submitted separately to these two abiotic stresses. It is also a valuable source of information to further study the osmotic effect. However, 8.85% of those enzymes had different qualitative profiles (Figure 5), where 29 enzymes were downregulated under salt stress and upregulated under drought (Figure 6a), and 20 had the opposite behavior (Figure 6b).     In the case of the proteomics analysis, a group of enzymes were also differentially expressed under both abiotic stresses ( Table 3). Most of them had the same qualitative profile independently of the stress condition and a minority had different ones. For instance, enzymes A0A6J0PP21 and A0A6I9S451, two isoforms coded by the same gene (LOC105053549) in the oil palm genome and described as L-ascorbate peroxidases in the KEGG database, were upregulated under drought stress and downregulated almost to no detection under salt stress.   Figure 7. Heatmap of the 13 differentially expressed metabolites that were found in the two scenarios studied-drought and salinity stresses.

Discussion
Plant response to water deprivation shows several similarities with the first phase of salinity stress, the osmotic effect [8][9][10]. The current study is a step further in a study in our lab that applies different omics platforms to characterize the morphophysiological and molecular responses of young oil palms to two abiotic stresses-salinity and drought [17,[20][21][22][23].
The young oil palm plants used in the study by Bittencourt et al. [17] and the ones used in the present study are clones of a specific plant-genotype AM33, a Deli x Ghana from ASD Costa Rica. Moreover, they were of almost similar age, cultivated under the same main environmental variables-as they were submitted to stress almost at the same time in the same greenhouse, and the stresses lasted 12 days in the case of salinity and 14 days in the case of drought, or almost the same time. Consequently, one can say that it is a unique opportunity to use all omics datasets-transcriptomics, proteomics, and metabolomics-from either one of the studies to map the similarities and dissimilarities in the response of those plants to those abiotic stresses. It is important to state that, in the case of salinity stress, plants were under field conditions throughout the entire experiment to ensure that the osmotic effect was due to the excess of salt in the soil and not otherwise [20].
In the present study, single-omics analyses showed that the most affected biological process subcategories in the transcriptome were protein phosphorylation, carbohydrate metabolic, and phosphorylation processes; meanwhile, in the proteome, the most affected ones were carbohydrate metabolism, cellular oxidant detoxification, and glycolytic processes. The transcriptomics and proteomics results, whether under salinity [17] or drought stress (present study), show that the most affected biological process subcategories in the transcriptome and proteome profiles are very much the same in the apical leaf of young oil palm plants.
Stress mitigation results from the action of several distinct proteins, allowing plants to survive under unfavorable conditions, including water deprivation. Post-translation modifications, such as protein phosphorylation, promote changes in those proteins that already existed in the cells before the onset of the stress and play roles important in plant resistance to abiotic stresses [24,25]. Carbohydrate metabolism also plays a role in plant stress perception and signal transduction, besides being a substrate for energy production and a mediator of osmotic regulation and carbohydrate distribution [26].
The number one molecular function affected by drought or salinity stress in the apical leaf of young oil palm plants, revealed by both the transcriptome and proteome profiles obtained here and by Bittencourt et al. [17], was ATP binding. ATP binding proteins (ABPs) have a binding site that allows adenosine-5'-triphosphate (ATP) molecules to interact, and they have roles important in membrane transport and the regulation of various metabolic processes. The ATP-binding cassette (ABC) transporters belong to a large protein family and are present in the membranes of chloroplasts, mitochondria, peroxisomes, and vacuoles. They participate in several biological processes of regulation in plants, including resistance to drought and salinity stresses, and are considered reliable targets for developing stresstolerant crop plants in the future [27]. Consequently, it was no surprise that the number one cellular component affected by the two abiotic stresses in question, either under drought or salinity stress, was an integral component of a membrane in the transcriptome profile and cytoplasm in the proteome profile.
The oxidoreductases, transferases, and hydrolases classes of enzymes had the highest number of positive hits, similar to what we saw in Bittencourt et al. [17]. Within those classes, one can find several subclasses equally affected by drought and salinity stress. Several successful cases of resistance to abiotic stresses achieved by the overexpression of genes coding for oxidoreductases, transferases, or hydrolases enzymes are available in the scientific literature [28][29][30], showing that here also resides reliable targets to generate stress-tolerant crop plants.
The integration of transcriptomics, proteomics, and metabolomics datasets to gain further insights into the mechanism behind the response of young oil palm plants to drought stress used a much larger number of differentially expressed molecules compared to the MOI study carried out previously by Bittencourt et al. [17], when a total of 510 enzymes and 19 metabolites, differentially expressed under salinity stress, were used. Here, 1955 DE enzymes from transcriptomics analysis, 131 DE enzymes from proteomics analysis, and 269 differentially expressed metabolites (DEMs) underwent MOI analysis. The consequent result was a significant increase in the number of unique DE molecules in the most affected pathways.
The present MOI study showed that the Cysteine and methionine metabolism was the most affected pathway in the apical leaf of young oil palm plants under drought stress, with 47 unique DE molecules. In Bittencourt et al. [17], this pathway was the number one in the rank of affected ones, tied together with Glycolysis/Gluconeogenesis pathway (map00010), both with 20 unique DE molecules. The latter turned out to be the seventh most affected pathway in the list of the most affected ones in young oil palm plants under drought stress, with 30 unique DE molecules. Glycine, serine, and threonine metabolism (map00260) and Amino sugar and nucleotide sugar metabolism (map00520) were present in both lists of most affected pathways, with 34 and 35 unique DE molecules, respectively, in plants under drought stress, and 14 and 12 DE molecules, respectively, under salinity stress [17]. At last, in the present study, when performing the MOI analysis only with the molecules (proteins and metabolites) differentially expressed in both scenarios-drought and salinity stresses-the Cysteine and methionine metabolism pathway again turned out to be the most affected.
Even though this study identified other pathways, further discussion will concentrate only on cysteine and methionine metabolism. Cysteine and methionine are sulfurcontaining amino acids. Besides being the metabolic sulfide donor for all cellular components containing reduced sulfur, cysteine is also a precursor of methionine, glutathione, phytochelatins, iron-sulfur clusters, vitamin cofactors, and multiple secondary metabolites, and the incorporation of sulfide into cysteine takes place in chloroplasts, mitochondria, and the cytoplasm [31]. Our results show that L-serine O-acetyltransferase, the first enzyme in cysteine biosynthesis, which catalyzes the synthesis of O-acetylserine (OAS) from acetyl-CoA (AcCoA) and L-serine (L-Ser), is upregulated in the apical leaf of young oil palm plants under saline (Log 2 (FC) = 3.4) as well as drought stress (Log 2 (FC) = 4.7). The same is true for the cysteine synthase, which catalyzes the synthesis of L-cysteine and acetate from OAS and hydrogen sulfide and showed a~2-fold and~6-fold increase, respectively, under saline and drought stresses.
Several enzymes from the omnipresent methionine salvage pathway (MSP), which uses dioxygen to regenerate methionine, associated with a ratchet-like step that prevents methionine back degradation [32], were differentially expressed in the apical leaf of young oil palm plants under saline and drought stresses. Enzymes from this pathway, such as S-adenosylmethionine synthetase, S-adenosylmethionine decarboxylase, and tyrosine aminotransferase, were positively regulated under both abiotic stresses. Meanwhile, 5-methylthioribose kinase and L-glutamine-4-(methylsulfanyl)-2-oxobutanoate aminotransferase were both downregulated. Resistance to salt and oxidative stress was achieved in Arabidopsis plants by the overexpression of a S-adenosylmethionine synthetase gene from salt-resistance sugar beet [33]. Barley stripe mosaic virus-based virus-induced gene silencing (BSMV-VIGS) of a S-adenosylmethionine synthetase gene in the Tibetan wild barley genotypes XZ5 severely compromised its resistance to both drought and salinity stress [34].
Besides mapping the commonalities in the molecular response of young oil palm plants to drought and salinity stresses, the present study also reported the molecules that displayed the opposite behavior when submitted separately to those two stresses. Stachyose synthetase (XP_010931766.1) is the enzyme that showed to most opposite behavior, downregulated under salt stress (Log 2 (FC) = −5.5) and upregulated in drought (Log 2 (FC) = 4.9). This enzyme catalyzes the reaction that produces myo-inositol plus stachyose from alpha-D-Galactosyl-(1->3)-1D-myo-inositol and Raffinose and was regulated positively in the drought-tolerant conifer Larix olgensis A. Henry under PEG-simulated drought stress [35]. According to Okemo et al. [36], the shoots of Tripogon loliiformis-an Australian resurrection plant-accumulated stachyose synthase transcripts and stachyose during drought stress. Stachyose triggers the formation of the hallmarks of plant apoptotic-like cell death in the desiccation-sensitive Nicotiana benthamiana but not the resilient T. loliiformis. Because of that, the stachyose synthase gene from this species is a potential candidate for improving stress resistance.
The present study, and the one by Bittencourt et al. [17], used transcriptome datasets from the apical leaf of young oil palm plants already showing extensive necrosis [20,23]. The reason stachyose synthase gene in E. guineensis increased approximately 32-fold under drought but decreased by almost 100% under saline stress, and considering that the transcriptome profiles used came from leaves with extensive necrosis, this is a question that needs further study to be answered.
The enzyme that showed to second most opposite behavior, downregulated under salt stress (Log 2 (FC) = −2.3) and upregulated in drought (Log 2 (FC) = 3.1), is a Topoisomerase II (XP_010931766.1). John et al. [37] over-expressed a topoisomerase II gene from tobacco in Nicotiana tabaccum and observed an increase in the resistance to salt stress.
The enzyme (XP_010934998.1) that showed the most opposite behavior, upregulated under salt stress (Log 2 (FC) = 1.9) and downregulated in drought (Log 2 (FC) = −3.1), is a member of the Cytochrome P450, E-class, group IV (IPR002403), protein family. Cytochrome P450s (CYPs) act as versatile catalysts and play a crucial role in the biosynthesis of secondary metabolites, antioxidants, and phytohormones in higher plants. The expression of several CYP genes is regulated in response to environmental stresses, with a role important in the crosstalk between abiotic and biotic stress responses, which make them candidate genes for engineering crop species resilient to biotic and abiotic stresses [38].

Plant Material, Growth Conditions, Experimental Design, and Drought Stress
The oil palm plants used in this study are clones of the ones used in the Bittencourt et al. [17] study. All plants-from both studies-came from the same embryogenic calluses. The young oil palm plants used in both studies were clones regenerated out of embryogenic calluses obtained from the leaves of an adult plant-genotype AM33, a Deli x Ghana from ASD Costa Rica; and were subjected to treatments when they were in the growth stage known as "bifid" saplings. Before starting the experiments, plants were standardized according to their developmental stage, size, and the number of leaves.
The experiment consisted of two water availability levels (field capacity-control and water deprivation-stressed), with four replicates in a completely randomized design. For the omics (transcriptomics, metabolomics, and proteomics) analysis described in the present study, we collected the apical leaves from control and stressed plants 14 days after imposing the treatments (DAT). See Rodrigues-Neto et al. [21] and Salgado et al. [23] for additional details regarding growth conditions, experimental design, drought stress, and the morphophysiological response of the oil palm plants.

Single-Omics Data Analysis, Functional Annotation, and Integratomics Analysis
Samples harvested from control and stressed plants were immediately immersed in liquid nitrogen and stored at −80 • C until RNA/Protein/Metabolite extraction; three replicates (plants) per treatment for RNA and protein extraction and four for the extraction of metabolites.
See Salgado et al. [22] for additional details on total RNA extraction, quality analysis, library preparation, and sequencing. The RNA-Seq raw sequence data are in the Sequence Read Archive (SRA) database of the National Center for Biotechnology Information (NCBI) under the BioProject number PRJNA573093. The transcriptomics analysis was performed using the OmicsBox platform-version 2.0.36 [39], as previously described [22]. The oil palm genome [18,19]-downloaded from NCBI (BioProject PRJNA268357; BioSample SAMN02981535) in September 2021-was the reference genome for RNA-Seq data alignment. The pairwise differential expression analysis between experimental conditions (Stressed Plants X Control) was performed as described [17].
The protein extraction protocol used here was the same used before [17]. The same happened for protein preparation and LC-MS/MS analysis at the GenOne company (Rio de Janeiro, RJ, Brazil). Data were acquired in biological triplicates using the the Xcalibur software (version 2.0.7) (Thermo Scientific, Waltham, MA, USA). Raw files were processed using the PatternLab for Proteomics V software [40], also as described by [17].
The metabolite extraction protocol used here was the same used by [17]. The Data-Analysis 4.2 software (Bruker Daltonics, Bremen, Germany) was the first used to analyze the raw data from UHPLC-MS as mzMXL files, and the pre-processing of data, peak detection, alignment of retention times and the statistical analysis at the pre-processing stage was similarly used by [17]. Then, a dataset was created from control and stressed plants. The pre-processed data (csv file) underwent analysis in the Statistical Analysis, MS Peaks to Pathway, and pathway analysis modules of the MetaboAnalyst 5.0 [41,42], exactly as previously described [17]. The differentially expressed peaks (DEPs) selected were those passing the criteria of false rate discovery (FDR) ≤ 0.05 and Log 2 (fold change [FC]) = 1.
The differentially expressed (DE) full-length ORFs and proteins identified in the OmicsBox and PatternLab V analysis underwent KOALA annotation before MOI analysis (Figure 1). The approach used to integrate the three omics was pathway-based mapping, performed using the Omics Fusion platform [43], as shown in Figure 1 and previously described [17].

Correlation Analysis of Differentially Expressed Molecules under Two Distinct Scenarios-Salinity and Drought Stresses
To perform a correlation analysis of differentially expressed (DE) full-length RNAs, proteins, and metabolites under two distinct scenarios, we used sets of data with the DE molecules from this present study and re-used the respective ones from [17]. First, to check the data distribution, we used the Data Overview module of Omics Fusion [43], the web platform for integrative analysis of Omics data, and then the Scatter Plot one for the correlation analysis between the sets of data-a pairwise combination of the different molecules and scenarios evaluated. The input data was the Log 2 (FC) from the DE molecules obtained from the single-omics analysis.

Conclusions
The SOA and MOI studies here reported generated new insights on the early response of young oil palm plants to drought stress, pointing out genes, proteins, metabolites, and pathways directly affected by this stress.
The MOI study performed using the molecules differentially expressed under drought (from this study) and salinity stress (from [14]) gave a unique opportunity to map the similarities and dissimilarities in the molecular response of those plants to those abiotic stresses of very high importance.
All MOI studies showed that cysteine and methionine metabolism was the most affected pathway, independent of the abiotic stress applied. This was also the case when only those molecules differently expressed under both abiotic stresses were considered.
The correlation analysis revealed that only 8.85% of those enzymes which were differently expressed under both stresses had different qualitative profiles, corroborating the already known fact that plant responses to drought and salinity show several similarities, with emphasis on the osmotic effect.
Finally, besides allowing us to map the pathways most affected by two abiotic stresses in the leaves of young oil palm plants, this study also generated a list of genes, proteins, and metabolites responsive to those stresses. This list can be used as valuable information to help those interested in the prospection and validation of candidate genes to promote resistance to both abiotic stresses once overexpressed in oil palm or other plant species.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/plants11202786/s1, Table S1: Proteins from the oil palm genome differentially expressed (FDR ≤ 0.05) in the leaves of young oil palm plants under drought stress. FC-Fold Change, CPM-Count Per Million, FDR-False Discovery Rate, Kegg Orthology number, and E.C. number; Table S2: List of all peptides confidently identified in oil palm conditions; Table S3: List of all proteins confidently identified in oil palm conditions; Table S4: All proteins confidently identified in oil palm conditions with maximum parsimony criterion; Table S5: List of the proteins detected only in control biological condition; Table S6: List of the proteins detected only in stressed biological condition; Table S7: List of all differentially expressed proteins detected in both biological conditions (Stressed and Control) with statistical significance (FDR ≤ 0.05); Table S8: List of metabolites identified in the leaves of young oil palm plants submitted to drought stress via metabolomics analysis, after submitting the differentially expressed (DE) peaks to the pathway topology analysis module in MetaboAnalyst 5.0; Table S9: group of 554 enzymes differentially expressed under both abiotic stresses-drought and salinity. Funding: This work was funded by the Brazilian Innovation Agency (Financiadora de Estudos e Projetos-FINEP-DendêPalm Project-01. 13.0315.00) and the Coordination for the Improvement of Higher Education Personnel-CAPES (Scholarships for C.B.B., T.L.C.d.S., Í.d.O.B. and L.R.V.). The funding bodies covered all student and project costs but were not involved in the design, collection, analysis and interpretation of the data or the preparation of the manuscript.

Data Availability Statement:
The datasets used and/or analyzed in the current study are available from the corresponding author on reasonable request.

Conflicts of Interest:
The authors declare no conflict of interest.