Meta-Analysis of Common and Differential Transcriptomic Responses to Biotic and Abiotic Stresses in Arabidopsis thaliana

Environmental stresses adversely affect crop growth and yield, resulting in major losses to plants. These stresses occur simultaneously in nature, and we therefore conducted a meta-analysis in this study to identify differential and shared genes, pathways, and transcriptomic mechanisms involved in Arabidopsis response to biotic and abiotic stresses. The results showed a total of 436/21 significant up-/downregulated differentially expressed genes (DEGs) in response to biotic stresses, while 476 and 71 significant DEGs were respectively up- and downregulated in response to abiotic stresses in Arabidopsis thaliana. In addition, 21 DEGs (2.09%) were commonly regulated in response to biotic and abiotic stresses. Except for WRKY45 and ATXTH22, which were respectively up-/down- and down-/upregulated in response to biotic and abiotic stresses, other common DEGs were upregulated in response to all biotic and abiotic treatments. Moreover, the transcription factors (TFs) bHLH, MYB, and WRKY were the common TFs in response to biotic and abiotic stresses. In addition, ath-miR414 and ath-miR5658 were identified to be commonly expressed in response to both biotic and abiotic stresses. The identified common genes and pathways during biotic and abiotic stresses may provide potential candidate targets for the development of stress resistance breeding programs and for the genetic manipulation of crop plants.


Introduction
Plants are continuously exposed to a wide range of adverse challenges under natural environmental conditions, which adversely affect plant growth, development, crop quality, and quantity. They are faced with various biotic stresses, including bacteria, fungi, viruses, herbivorous insects, and nematodes [1]. Moreover, abiotic stresses, such as salinity, drought, and extreme temperature, affect many physiological processes in plants [2]. The regulation of proper responses is crucial to reduce the damage caused by these stresses [3]. Various stresses initiate highly complex responses that lead to specific modifications at the transcriptome, cellular, and physiological levels. The plant stress response to adverse conditions includes complex interactions among their respective signaling pathways [4,5]. Environmental factors, such as drought, extreme temperature, and salinity, potentially alter the outcome of plant-pathogen interactions by modulating plant responses to the Our results therefore highlighted that the most common DEGs in both biotic and abiotic stresses were upregulated. We conclude that a number of common pathways are induced by biotic and abiotic stresses. The common genes may play a role in the resistance against biotic and abiotic stresses and may be considered potential candidates in plant breeding programs. The DEGs including AT2G44890 (cytochrome P450) and AT1G09935 (phosphoglycerate mutase family protein) were the most upregulated ones in response to biotic stresses, while AT5G59720 (heat shock protein 18.2) and AT1G72660 (p-loop-containing nucleoside triphosphate hydrolases superfamily protein) were the highest upregulated DEGs in response to abiotic stresses. These genes showed the most upregulated expression in response to biotic and abiotic stresses. Previous studies have shown that cytochrome P450s play a key role in the biosynthesis and regulation of antioxidants, fatty acids, cell wall components, biopolymers, defense compounds, and phytohormones in plants [23,24]. These mediate plant defense responses and protect plants from biotic stresses including insects and diseases [25,26]. The expression of a phosphoglycerate mutase homologue in Arabidopsis is also reportedly regulated by hormones and induced by cyst and root-knot nematodes [27]. Heat shock protein 18.2 was highly expressed during abiotic stresses, which is in agreement with another study, and was found to be upregulated in response to heat, cold, and drought stresses [28,29]. Another highly expressed gene encoding a p-loop-containing nucleoside triphosphate hydrolase superfamily protein has been shown to play a role in salt tolerance in transgenic A. thaliana during abiotic stresses [30].

Gene Ontology
GO enrichment analysis was performed to identify the function of the differential and common up-/downregulated DEGs in response to biotic and abiotic stresses. A total Our results therefore highlighted that the most common DEGs in both biotic and abiotic stresses were upregulated. We conclude that a number of common pathways are induced by biotic and abiotic stresses. The common genes may play a role in the resistance against biotic and abiotic stresses and may be considered potential candidates in plant breeding programs. The DEGs including AT2G44890 (cytochrome P450) and AT1G09935 (phosphoglycerate mutase family protein) were the most upregulated ones in response to biotic stresses, while AT5G59720 (heat shock protein 18.2) and AT1G72660 (p-loop-containing nucleoside triphosphate hydrolases superfamily protein) were the highest upregulated DEGs in response to abiotic stresses. These genes showed the most upregulated expression in response to biotic and abiotic stresses. Previous studies have shown that cytochrome P450s play a key role in the biosynthesis and regulation of antioxidants, fatty acids, cell wall components, biopolymers, defense compounds, and phytohormones in plants [23,24]. These mediate plant defense responses and protect plants from biotic stresses including insects and diseases [25,26]. The expression of a phosphoglycerate mutase homologue in Arabidopsis is also reportedly regulated by hormones and induced by cyst and root-knot nematodes [27]. Heat shock protein 18.2 was highly expressed during abiotic stresses, which is in agreement with another study, and was found to be upregulated in response to heat, cold, and drought stresses [28,29]. Another highly expressed gene encoding a p-loop-containing nucleoside triphosphate hydrolase superfamily protein has been shown to play a role in salt tolerance in transgenic A. thaliana during abiotic stresses [30].

Gene Ontology
GO enrichment analysis was performed to identify the function of the differential and common up-/downregulated DEGs in response to biotic and abiotic stresses. A total of 457 and 547 DEGs in the samples exposed to biotic and abiotic stresses were, respectively. analyzed to conduct GO analysis. GO analysis grouped the upregulated DEGs into the categories of biological process, cellular component, and molecular function. The highest frequencies of biological process-related DEGs in the samples under biotic stress belonged to the response to stimulus (26.48%), biological regulation (17.58%), and response to stress (16.44%) categories. Concurrently, cellular (37.85%) and metabolic (34.66%) processes showed the highest number of biological process-related DEGs in response to abiotic stress. In the cellular component category, most of the upregulated genes under biotic stress were related to the cell part (48.4%) and cell (48.4%) categories, while cell part (50.2%) and cell (50.2%) comprised the most significantly upregulated genes of the cellular component during abiotic stress. The large functional group of DEGs related to cellular processes highlights a rearrangement in plant metabolism in response to stress adaptation [31]. Furthermore, catalytic activity (38.36%) and binding (37.21%) revealed the most enriched numbers of molecular function group in the samples under biotic stress. On the other hand, DNA binding (11.55%), transcription regulator (10.76%), and transcription factor (10.36%) activities showed the highest number of molecular function-related DEGs in the samples under abiotic stress ( Figure 2). of 457 and 547 DEGs in the samples exposed to biotic and abiotic stresses were, respectively. analyzed to conduct GO analysis. GO analysis grouped the upregulated DEGs into the categories of biological process, cellular component, and molecular function. The highest frequencies of biological process-related DEGs in the samples under biotic stress belonged to the response to stimulus (26.48%), biological regulation (17.58%), and response to stress (16.44%) categories. Concurrently, cellular (37.85%) and metabolic (34.66%) processes showed the highest number of biological process-related DEGs in response to abiotic stress. In the cellular component category, most of the upregulated genes under biotic stress were related to the cell part (48.4%) and cell (48.4%) categories, while cell part (50.2%) and cell (50.2%) comprised the most significantly upregulated genes of the cellular component during abiotic stress. The large functional group of DEGs related to cellular processes highlights a rearrangement in plant metabolism in response to stress adaptation [31]. Furthermore, catalytic activity (38.36%) and binding (37.21%) revealed the most enriched numbers of molecular function group in the samples under biotic stress. On the other hand, DNA binding (11.55%), transcription regulator (10.76%), and transcription factor (10.36%) activities showed the highest number of molecular function-related DEGs in the samples under abiotic stress ( Figure 2). The significant upregulated DEGs involved in biological processes (response to stimulus (GO:0050896), biological regulation (GO:0065007), and response to stress (GO:0006950)) during biotic stress treatment were assigned to the following functions: defense response to pathogens (viruses, fungi, oomycetes, and bacteria), response to stress, protein kinase activity, cellulase and beta-glucosidase activity, MAP kinase activity, innate immune response, camalexin biosynthetic process, protein phosphorylation, calmod- The significant upregulated DEGs involved in biological processes (response to stimulus (GO:0050896), biological regulation (GO:0065007), and response to stress (GO:0006950)) during biotic stress treatment were assigned to the following functions: defense response to pathogens (viruses, fungi, oomycetes, and bacteria), response to stress, protein kinase activity, cellulase and beta-glucosidase activity, MAP kinase activity, innate immune response, camalexin biosynthetic process, protein phosphorylation, calmodulin binding, induced systemic and systemic acquired resistance, signal transduction, circadian rhythm, leaf senescence, DNA-binding transcription factor activity, growth regulation, gene expression regulation, response to oxidative stress, protein ubiquitination, response to hormones, and programmed cell death. On the other hand, significant upregulated DEGs were involved in the following biological processes: cellular (GO:0009987) and metabolic (GO:0008152) processes in the abiotic stress treatment assigned to functions, including proteolysis involved in the cellular protein catabolic process, DNA-binding transcription factor activity, response to salt stress, response to desiccation, response to hormones, response to heat and cold, photosynthesis, protein ubiquitination, protein phosphorylation, heat and cold acclimation, response to water deprivation, regulation of cellular response to alkaline pH, protein binding, response to oxidative stress, chaperone binding, cellular respiration, response to osmotic stress, drought recovery, response to freezing and salt stress, protein stabilization, heat shock protein binding, proline biosynthetic process, abscisic acid biosynthetic process, signal transduction, malate dehydrogenase (decarboxylating) (NADP+) activity, lipid metabolic process, protein kinase activity, regulation of cellular protein metabolic process, lipid oxidation, circadian rhythm, and protein folding. In the cellular component, the most enriched GO in terms of the cell part and cell during biotic and abiotic stresses were signal transduction, cellular response to oxidative stress, cellular oxidant detoxification, respiratory burst, cell redox homeostasis, response to oxidative stress and reactive oxygen species, protein phosphorylation, defense response, kinase activity, cellular response to freezing, cell wall modification, pectin esterase inhibitor activity, chitinase activity, PAMPtriggered immunity signaling pathway, glutathione metabolic process, programmed cell death, lipid and carbohydrate metabolic process, pattern recognition receptor signaling pathway, proteolysis involved in cellular protein catabolic process, regulation of reactive oxygen species metabolic process, chaperone-mediated protein folding, cytochrome-c oxidase activity, mitochondrial electron transport, electron transfer activity, regulation of cellular response to alkaline pH, proline biosynthetic process, ATP synthesis-coupled electron transport, and gene silencing by RNA-directed DNA methylation. In addition, the significant upregulated DEGs involved in molecular functions, including catalytic activity (GO:0003824) and binding (GO:0005488) in the biotic stress treatment, contained activities such as a defense response to pathogens, response to stress, camalexin biosynthetic process, induced systemic resistance and systemic acquired resistance, protein binding, response to hormones, hormone-mediated signaling pathway, response to oxidative stress, protein kinase activity, protein phosphorylation, protein ubiquitination, cellulase activity, signal transduction, regulation of growth, MAP kinase activity, RNA binding, transcription regulatory region sequence-specific DNA binding, DNA-binding transcription factor activity, regulation of gene expression, RNA processing, respiratory burst, programmed cell death, production of siRNA involved in RNA interference, methylation, metalloendopeptidase activity, gene silencing by miRNA, ribonuclease III activity, innate immune response, signal transduction, glutathione transferase activity, and calmodulin binding. The most enriched DEGs for molecular functions, including DNA binding (GO:0003677), transcription regulator (GO:0030528), and transcription factor (GO:0003700) activities in the abiotic stress treatment contained functions such as response to stress, response to oxidative stress, protein ubiquitination, transcription regulatory region sequence-specific DNA binding, DNA-binding transcription factor activity, hyperosmotic salinity response, regulation of transcription and gene expression, response to hormones, response to water deprivation, heat acclimation, RNA secondary structure unwinding, RNA binding, mRNA export from the nucleus, alternative mRNA splicing, transcription coactivator activity, and the mitochondria-nucleus signaling pathway. Our findings are in agreement with other studies which showed that the large functional groups of DEGs were related to metabolic processes, regulatory function, and response to stimulus [16,31].

KEGG Analysis of DEGs
A total of five significant pathways, i.e., phenylpropanoid biosynthesis, plant-pathogen interaction, glutathione metabolism, amino sugar and nucleotide sugar metabolism, and tropane, piperidine, and pyridine alkaloid biosynthesis, were enriched with upregulated DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis during biotic stresses, whereas four significant enriched pathways in response to abiotic stresses were protein processing in the endoplasmic reticulum (ER), oxidative phosphorylation, plant hormone signal transduction, and carotenoid biosynthesis (Table 1). Among the pathways Plants 2022, 11, 502 6 of 15 characteristic of biotic stresses, phenylpropanoid and glutathione have been shown to mediate resistance to plant pathogens. Phenylpropanoids play a role in plant responses to biotic stresses and disease resistance [32,33]. An increase in phenylalanine content affects the biosynthesis of phenylpropanoid compounds, leading to increased resistance to oxidative stress and fungal infection [33]. Moreover, glutathione plays a crucial role in plant responses during biotic stresses and plant resistance to pathogens [34]. DEGs annotated to protein processing in the ER, oxidative phosphorylation, and plant hormone signal transduction were upregulated in response to abiotic stresses. The ER regulates stress responses in plants and is related to abiotic stress tolerance [35,36], mitigating stress damage and providing stress tolerance to plants [37]. ER stress signaling is involved in the response to abiotic stresses, including heat and osmotic stresses [38]. Oxidative phosphorylation produces adenosine triphosphate (ATP) required for providing the necessary energy for plant survival and development [39,40]. Cell death has been suggested to be controlled through various pathways, e.g., the perturbation of oxidative phosphorylation [41]. Phytohormones regulate a number of plant cellular processes, such as abiotic stress tolerance and plant defense response, and coordinate signal transduction pathways during an abiotic stress response. They are thereby considered potential targets for producing abiotic stress-tolerant plants [42].

Differential Stress-Responsive Groups
Here, we functionally grouped the highest upregulated genes in response to biotic and abiotic stresses. Our results revealed that the stress response types varied in both biotic and abiotic stresses, highlighting different stress-responsive pathways involved in plantstress interactions. Heat shock proteins (HSPs) and heat shock transcription factors (HSFs) were among the most upregulated responsive groups under abiotic stress. The role of HSPs has previously been reported in the development of tolerance against various abiotic stresses [43], which acts through the regulation of protein folding and unfolding, subcellular localization of proteins, signal transduction, transcriptional activation and degradation of unfolded proteins, resulting in protein quality control and cellular homeostasis [43][44][45]. HSFs are also considered one of the key components of the signal transduction cascade causing the regulation and expression of genes involved in abiotic stress responses [46]. On the other hand, a number of specific stress-responsive groups were revealed during biotic stresses. Calcium (Ca) and calmodulin (CaM)-related responses were found to be upregulated in response to biotic stresses. Ca and CaM are crucial messengers in stress signaling by eliciting responses to stresses [47]. In addition, glutathione S-transferases (GSTs)-related DEGs were upregulated during biotic stresses. Our findings are in agreement with other studies that showed GST genes to be specifically induced in response to plant pathogens, which may contribute to disease resistance [48]. Some GSTs are remarkably induced in the early stages of bacterial, viral, and fungal infections. GSTs detoxify toxic substances, attenuate oxidative stress, and participate in hormone transport [48]. Furthermore, we identified induced DEGs of pathogenesis-related (PR) proteins that were responsive to phytopathogens. PR proteins play key roles in plant innate immune systems and cause enhanced resistance against biotic stresses [49] via inducing plant cell wall reinforcement, suppressing pathogen enzymes, degrading PR proteins, and inducing defense mechanisms and their direct antimicrobial activity against pathogens [50]. In addition, our findings showed that a number of membrane-and cell wall-related DEGs were induced under biotic stresses. Induced Ca-related DEGs may be responsible for the biochemical modification of the membrane and the cell wall. Ca plays a key role in maintaining the structural rigidity of cell walls and membranes [51]. Pattern recognition receptors perceive pathogen-associated molecular patterns in the plasma membrane and initiate pattern-triggered immunity (PTI) during biotic stresses [11,52,53]. PTI can be triggered by host-derived molecules from the cell wall [54]. Cell wall proteins play crucial roles in plants during stress [55]. Plasma membrane-localized receptor-like kinases and receptor-like proteins such as pattern recognition receptors play a critical role in adapting to biotic stresses [56]. In addition, plants sense the pathogens and modify their cell wall components to induce a defense response [57].

Identification of DEGs Encoding Transcription Factors in Response to Biotic/Abiotic Stresses
In the present study, a total of 11 and 15 TF families were significantly upregulated in response to biotic and abiotic stresses, respectively (Table S1). WRKY family showed the highest TF number (15) in response to biotic stresses, whereas AP2/EREBP family included the most presented TF (15) during abiotic stresses ( Table 2). WRKY TFs are a large family of regulatory proteins and are involved in plant defense responses to biotic and abiotic stresses [58][59][60]. Bacterial, viral, and fungal infections can activate WRKY TF genes, which subsequently mediate resistance to plant pathogens [61][62][63][64][65][66]. On the other hand, AP2/EREBP, a large family of TF genes in plants [67], is involved in various abiotic stress responses, including genes related to drought [68], low temperature [69], and high salt concentrations [68]. In addition, TFs, i.e., bHLH, C2H2, C3H, CCAAT-HAP2, HSF, MYB, NAC, and WRKY were the common TFs in response to biotic and abiotic stresses. On the other hand, C2C2-CO-like, C2C2-Dof, G2-like, GRAS, Homeobox, and RAV were specifically detected in response to abiotic stresses, while TFs including REM and Trihelix were only identified against biotic stresses.

Protein-Protein Interaction
In this study, we constructed a PPI network and identified the most significant hub genes. The network showed 259/333 nodes and 1283/2189 edges during biotic and abiotic stresses, respectively. The most significant hubs in the biotic stress network were transmembrane protein (AT2G18690), calmodulin-binding protein (CBP60G), copper transport protein family (AT5G52760), hypothetical protein (XBAT34), and mitogen-activated protein kinase 11 (MPK11), while a mitochondrial ribosomal protein (ATMg00080), ribosomal protein L2 (RPL2), ATP synthase CF0 A subunit (ATPI), NADH dehydrogenase subunit (NDHI), and NADH dehydrogenase subunit 1 (NDHA) represented the most significant nodes in the abiotic stress network (Figure 3). CaM-binding proteins are involved in biotic and abiotic stress responses. Among them, CBP60g regulates SA biosynthesis induced by microbe-associated molecular patterns. The overexpression of CBP60g in Arabidopsis increased SA accumulation, tolerance to drought stress, expression of the defense genes, and resistance to Pseudomonas syringae [70]. Mitogen-activated protein kinases (MAPK) cause cellular signal transduction during stress responses. Treatment with flg22, a pathogen-associated molecular pattern, was found to increase MPK11 expression in A. thaliana [71]. Copper (Cu) acts as a redox-active cofactor in a number of plant proteins and is required for a wide number of critical biological processes including mitochondrial respiration, cell wall metabolism, oxidative stress protection, response to pathogens, and ethylene perception [72][73][74]. Cu is suggested CaM-binding proteins are involved in biotic and abiotic stress responses. Among them, CBP60g regulates SA biosynthesis induced by microbe-associated molecular patterns. The overexpression of CBP60g in Arabidopsis increased SA accumulation, tolerance to drought stress, expression of the defense genes, and resistance to Pseudomonas syringae [70]. Mitogen-activated protein kinases (MAPK) cause cellular signal transduction during stress responses. Treatment with flg22, a pathogen-associated molecular pattern, was found to increase MPK11 expression in A. thaliana [71]. Copper (Cu) acts as a redox-active cofactor in a number of plant proteins and is required for a wide number of critical biological processes including mitochondrial respiration, cell wall metabolism, oxidative stress protection, response to pathogens, and ethylene perception [72][73][74]. Cu is suggested to play a role in antimicrobial activity through the alteration and inhibition of protein activity, the denaturation of nucleic acids, and the modulation of plasma membrane permeabilization [75]. Xanthomonas oryzae has been shown to overcome rice defenses by regulating host Cu redistribution [74].
Mitochondrial F(1)F(0)-ATPase, an important enzyme in plant metabolism, provides ATP for cells. The overexpression of a mitochondrial ATP synthase small subunit gene reportedly confers resistance to several abiotic stresses (salts, drought, oxidative and cold stresses) in A. thaliana, suggesting F(1)F(0)-ATPase to play a role in stress tolerance [76]. ATP synthase subunits (e.g., ATPase F0) and NADH dehydrogenase were also found to change in plants under drought stress [77,78], and these changes may result in an increase of energy production in response to water deficit to inhibit stress damage [79]. NADH dehydrogenase enhances tolerance to environmental stresses, including drought and elevated light, and is thought to prevent cellular damage [80].

Identification of miRNAs Targeting Downregulated Genes
miRNAs are an extensive class of non-coding RNAs that play a key role in the regulation of various biological processes through post-transcriptional regulation of gene expression. In this study, we identified 84 miRNAs belonging to 33 miRNA families (Table S2), a high proportion of which are involved in various biological processes.
The analysis of miRNAs during biotic and abiotic stresses exhibited 5 and 25 miRNA families, respectively, showing their role in the downregulation of A. thaliana genes in response to stresses.
ath-miR172 was predicted to have the highest frequency among downregulated gene targets, with eight targets during biotic stresses, whereas ath-miR858 and ath-miR5665 were identified with the highest targets (eight and six targets) during abiotic stresses, significantly showing the highest number of miRNAs targeting downregulated genes. In addition, ath-miR414 and ath-miR5658 were identified to be commonly expressed in response to both biotic and abiotic stresses. Similarly, miR172 was found to significantly downregulate a number of gene targets during viral infections in A. thaliana [81]. miR172 also reportedly regulates genes during biotic and abiotic stresses [82,83]. MiR858 has been shown to target regulatory factors, resulting in the regulation of a number of plant functions, including hormonal and stress responses [84]. In addition, the role of miR858 has been established in drought stress response [85], and this miRNA was upregulated in response to salt stress [86]. Moreover, two common microRNAs (miR414, miR5658) detected in response to biotic and abiotic stresses were the key regulators of salt-responsive genes [87]. MiR5658 plays prominent roles in the regulation of processes involved in responses to biotic and abiotic stresses [88][89][90]. Our results may shed light on miRNA functions in plant responses to biotic and abiotic stresses.

Data Collection and Preprocessing
The raw expression data (RNA-seq) for biotic and abiotic responses were obtained from the ArrayExpress of the European Bioinformatics Institute (http://www.ebi.ac.uk/ arrayexpress; 15 September 2020). The data sets were filtered for Arabidopsis thaliana, to only include RNA-seq data. The search yielded 12 entries in 11 papers, as detailed in Table 3. After trimming the adaptor sequences and low-quality reads based on their Phred quality scores [91], we mapped the checked reads onto the Arabidopsis reference genome. Expression profiling analysis was carried out using CLC Genomic Workbench version 10 (CLC Bio, Qiagen). The raw expression data pertaining to each data set were normalized with counts per million (CPM). The workflow is presented in Figure 4.

Meta-Analysis of an Expression Dataset
We conducted a meta-analysis independently on an integrated dataset pertaining to DEGs in biotic/abiotic stresses. Each dataset was grouped into a stress class and a healthy class according to the type of stress involved. Before the meta-analysis, the SVA R package was applied to correct for the batch effect, according to the empirical Bayes method [103]. Fisher's method was used for detecting DEGs (FDR < 0.01) [104]. The log ratio of means (ROM) was used to measure gene expression level [105]. ROM was calculated using the following equation: where y gn , r gr , and r gs represent ROM and mean expression levels of the stress and healthy classes for each gene in the dataset, respectively. Data were preprocessed and analyzed using Bioconductor packages (http://www.bioconductor.org; 15 September 2020), including MetaMA.

Meta-Analysis of an Expression Dataset
We conducted a meta-analysis independently on an integrated dataset pertaining to DEGs in biotic/abiotic stresses. Each dataset was grouped into a stress class and a healthy class according to the type of stress involved. Before the meta-analysis, the SVA R package was applied to correct for the batch effect, according to the empirical Bayes method [103]. Fisher's method was used for detecting DEGs (FDR < 0.01) [104]. The log ratio of means (ROM) was used to measure gene expression level [105]. ROM was calculated using the following equation: where ygn, rgr, and rgs represent ROM and mean expression levels of the stress and healthy classes for each gene in the dataset, respectively. Data were preprocessed and analyzed using Bioconductor packages (http://www.bioconductor.org; 15 September 2020), including MetaMA.

Gene Enrichment Analysis and Functional Analysis
The genes were analyzed after being selected by the meta-analysis. In the first instance, the co-occurrence of DEG with the resulting list was represented with Venn's diagrams using Venny 2.0 (http://bioinfogp.cnb.csic.es/tools/venny/; 15 September 2020). Enrichment analysis of Gene Ontology (GO) was performed on significant DEGs, which had been obtained from the meta-analysis. The enrichment analysis was conducted using the AgriGO platform [106]. Information on gene ontology was extracted based on GO terms for biological processes, cellular components, and molecular functions when a significant threshold of FDR < 0.05 was obtained. The pathway analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) was used to elucidate the enriched pathways of the DEGs (http://david.abcc.ncifcrf.gov/; 15 September 2020). To identify transcription factors among the DEGs, a list of Arabidopsis transcription factors was obtained from the AGRIS database (https://agris-knowledgebase.org/AtTFDB/; 15 September 2020).

Protein-Protein Interactions and Network Construction
The network analysis of protein-protein interactions (PPI) was performed to uncover any plausible interactions among proteins with significantly different DEGs. The STRING database was employed to enable the PPI network analysis [107], and Cytoscape software was used for visualizing the interaction networks.

Prediction of Potential miRNAs
The identification of miRNAs associated with biotic and abiotic stresses is necessary for understanding how plants respond defensively to stresses. Identifying potential miR-NAs is possible with the psRNATarget server (http://plantgrn.noble.org/psRNATarget/; 15 September 2020), wherein we set the parameters to default, except for maximum expectation, which was set to 3 in this study.

Conclusions
This study provides valuable information on A. thaliana in terms of differential and common host responses against biotic and abiotic stresses. As plants are simultaneously exposed to both biotic and abiotic stresses in the field and in natural ecosystems, the information generated here will help us identify common and specific genes and pathways during biotic and abiotic stresses. The common and specific genes, as potential candidate targets, may provide new insights for developing stress-resistant breeding and the genetic manipulation of crop plants via integrating biotic and abiotic resistance traits. Defenserelated pathways will be determined based on this, and they may simultaneously develop plant resistance to biotic and abiotic stresses and improve crop productivity. Further studies are required to elucidate the molecular mechanisms and validate the function of responsive DEGs, TFs, and miRNAs that regulate plant responses to biotic and abiotic stresses.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/plants11040502/s1, Table S1: List of the identified DEGs by the meta-analysis, Table S2: Gene ontology enrichment analysis of DEGs under abiotic and biotic stresses, Table S3: Identified transcription factors in response to biotic and abiotic stresses, Table S4: miRNAs associated with the DEGs were discovered using the computational algorithm psRNATarget server.