# De Novo Transcriptome Assembly, Gene Annotations, and Characterization of Functional Profiling Reveal Key Genes for Lead Alleviation in the Pb Hyperaccumulator Greek Mustard (Hirschfeldia incana L.)

Said El Hasnaoui, Mouna Fahr, Mohamed Zouine, Abdelaziz Smouni

► To cite this version:

Said El Hasnaoui, Mouna Fahr, Mohamed Zouine, Abdelaziz Smouni. De Novo Transcriptome Assembly, Gene Annotations, and Characterization of Functional Profiling Reveal Key Genes for Lead Alleviation in the Pb Hyperaccumulator Greek Mustard (Hirschfeldia incana L.). Current Issues in Molecular Biology, 2022, 44 (10), pp.4658-4675. 10.3390/cimb44100318. hal-04721322

HAL Id: hal-04721322

https://ut3-toulouseinp.hal.science/hal-04721322v1

Submitted on 4 Oct 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# De Novo Transcriptome Assembly, Gene Annotations, and Characterization of Functional Profiling Reveal Key Genes for Lead Alleviation in the Pb Hyperaccumulator Greek Mustard (*Hirschfeldia incana* L.)

Said El Hasnaoui <fnref n="1" />,<fnref n="2" /> ID, Mouna Fahr <fnref n="1" />,<fnref n="2" /> ID, Mohamed Zouine <fnref n="3" />,* and Abdelaziz Smouni <fnref n="1" />,<fnref n="2" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Hasnaoui, S.E.; Fahr, M.; Zouine, M.; Smouni, A. De Novo Transcriptome Assembly, Gene Annotations, and Characterization of Functional Profiling Reveal Key Genes for Lead Alleviation in the Pb Hyperaccumulator Greek Mustard (*Hirschfeldia incana* L.). *Curr. Issues Mol. Biol.* **2022**, *44*, 4658–4675.
https://doi.org/10.3390/cimb44100318

Academic Editors: Quan Zou and Ran Su

Received: 29 August 2022

Accepted: 27 September 2022

Published: 4 October 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Lead (Pb) contamination is a widespread environmental problem due to its toxicity to living organisms. *Hirschfeldia incana* L, a member of the Brassicaceae family, commonly found in the Mediterranean regions, is characterized by its ability to tolerate and accumulate Pb in soils and hydroponic cultures. This plant has been reported as an excellent model to assess the response of plants to Pb. However, the lack of genomic data for *H. incana* hinders research at the molecular level. In the present study, we carried out RNA deep transcriptome sequencing (RNA-seq) of *H. incana* under two conditions, control without Pb(NO₃)₂ and treatment with 100 µM of Pb(NO₃)₂ for 15 days. A total of 797.83 million reads were generated using Illumina sequencing technology. We assembled 77,491 transcript sequences with an average length of 959 bp and N50 of 1330 bp. Sequence similarity analyses and annotation of these transcripts were performed against the *Arabidopsis thaliana* nr protein database, Gene Ontology (GO), and KEGG databases. As a result, 13,046 GO terms and 138 KEGG maps were created. Under Pb stress, 577 and 270 genes were differentially expressed in roots and aboveground parts, respectively. Detailed elucidation of regulation of metal transporters, transcription factors (TFs), and plant hormone genes described the role of actors that allow the plant to fine-tune Pb stress responses. Our study revealed that several genes related to jasmonic acid biosynthesis and alpha-linoleic acid were upregulated, suggesting these components' implication in *Hirschfeldia incana* L responses to Pb stress. This study provides data for further genomic analyses of the biological and molecular mechanisms leading to Pb tolerance and accumulation in *Hirschfeldia incana* L.

**Keywords:** *Hirschfeldia incana* L.; Pb stress; RNA-seq; de novo assembly; transcriptome; jasmonic acid

## 1. Introduction

Lead (Pb) is known as one of the prominent heavy metals that occurs in the environment both from natural sources and mainly from anthropogenic activities [1]. Pb is known to be toxic to plants, animals, and microorganisms [1]. As one of the non-essential elements, Pb has no biological function and is toxic to living organisms even at low concentrations [2]. In plants, Pb causes various morphological, physiological, and biochemical effects [3]. However, several plant species have developed various tolerance strategies that are used to mitigate Pb toxicity [2]. These plant species can reduce the levels of Pb toxicity by the induction of both enzymatic and non-enzymatic pathways [4]. In *Brassica juncea*, it was

<footnote n="1">¹ Laboratoire de Biotechnologie et Physiologie Végétales, Centre de Biotechnologie Végétale et Microbienne Biodiversité et Environnement, Faculté des Sciences, Université Mohammed V de Rabat, Rabat 10000, Morocco</footnote>

<footnote n="2">² Laboratoire Mixte International Activité Minière Responsable “LMI-AMIR”, IRD/UM5/INAU, Rabat 10000, Morocco</footnote>

<footnote n="3">³ Laboratoire de Recherche en Sciences Végétales, Université de Toulouse, CNRS, UPS, Toulouse INP, 31320 Auzeville-Tolosane, France</footnote>

<footnote>* Correspondence: mohamed.zouine@toulouse-inp.fr (M.Z.); a.smouni@um5r.ac.ma (A.S.); Tel.: +212-6-61-77-16-62 (A.S.)</footnote>
reported that Pb toxicity increases SOD (superoxide dismutase) and APX (ascorbate peroxidase) levels with a decrease in CAT (Catalase) activity [5]. Hyper-tolerant species such as *Festuca ovina*, *Silene vulgaris*, *Noccaea caerulescens*, and *Peganum harmala* L. absorb and accumulate high levels of Pb in their root tissues and restrict Pb translocation to the aboveground parts [6,7]. Interestingly, plants with the capacity to hyper-accumulate Pb can accumulate concentrations higher than 1000 mg Kg⁻¹ of Pb in their aerial parts without showing toxicity symptoms [8]. *Hirschfeldia incana* L. commonly known as shortpod mustard is a perennial plant belonging to the Brassicaceae family currently found in the Mediterranean area and known for its capacity to grow and survive on highly heavy-metal-contaminated soils [9,10]. This plant has been identified and characterized as a Pb hyper-accumulator and a good model to evaluate the response of plants to Pb at physiological and molecular levels [8,10,11]. In addition, *H. incana* accumulates higher than 3 percent dry weight of Pb in its aerial parts when cultivated in a hydroponic medium containing 100 µM of Pb(NO₃)₂ [8,10]. Several plants have developed and formed specific physiological and molecular mechanisms to cope with Pb stress successfully. To date, many scientific reports have been published about the effects of Pb on plants and their tolerance mechanisms [2], but limited data are available on molecular mechanisms governing Pb uptake, translocation, and detoxification in plants [11]. Some molecular actors in Pb homeostatic processes have been identified including some transporters (*NtCBP4*, *AtCNGC1*, *AtATM3*, *AtPDR12*, *AtPDR8* ...) [12–17] and chelators (MTs, PCs, and GSH) [18]. Using heterologous microarray technology, several specific genes were identified as expressed under Pb exposure in *H. incana* [19]. Among these genes, *HiMT2b*, a MTs gene, was involved directly in Pb accumulation in aerial parts. Several genes implicated in ABA biosynthesis were reported to be overexpressed in roots and shoots of *H. incana* [19]. Moreover, research on *H. incana* is limited, and molecular mechanisms behind its Pb tolerance and accumulation are unclear due to the lack of knowledge on the genome of this plant. De novo assembly of transcriptome makes the study of RNAs feasible, even without a reference genome, by providing a large search transcriptome, characterizing and identifying genes that are expressed in roots and aboveground parts [20]. In this study, we first realized de novo assembly of RNA-seq transcriptome and gene annotation of *Hirschfeldia incana* L. transcriptome for the first time. We aimed to (i) explore the transcriptome of *H. incana* L. and (ii) find and characterize transcriptional changes that occur in this species in response to Pb stress.

## 2. Materials and Methods

### 2.1. Plant Material

Seeds of *Hirschfeldia incana* were collected from Oued El Heimer Pb mining site (34°26′88″ N–1°54′03″ W) which is located in the south of Oujda city in the eastern region of Morocco [21]. Pb highly contaminates the soil of this area with concentrations higher than 18,626 mg kg⁻¹ DW [8,11,22]. Two-week-old *H. incana* seedlings were cultivated in hydroponic Broughton and Dilworth (BD) medium [22] for 15 days. After acclimation, *H. incana* plants were transferred to BD medium (without phosphate) and added with 100 µM of Pb(NO₃)₂. This concentration used was based on the physiological and molecular responses of *H. incana* under Pb stress [8,10,20]. Control plants were cultivated in free phosphate BD medium without Pb. Plants (3 replicates per treatment) were harvested after 15 days of treatment for further experiments (Figure 1).

### 2.2. mRNA Extraction

mRNA was extracted from the roots and aboveground parts of *H. incana* plants. mRNA extraction was carried out by using the ReliaPrep™ RNA Tissue miniprep System (Promega, USA) kit. A post-treatment with an Invitrogen DNA free™ kit (Thermo Fisher Scientific, Waltham, MA, USA) eliminates DNA contamination. The concentrations of mRNA were assayed using Nanodrop 2000 (Thermo Scientific, Waltham, MA, USA). The quality of mRNA samples was assayed using an Agilent 2100 Bioanalyzer (Agilent
Technologies Canada, Inc., Mississauga, ON, Canada), and only mRNA samples with RIN ≥7 were used for further analyses (Figure 1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Experimental design for RNAseq and pipeline.</figcaption></figure>

## 2.3. cDNA Construction and Sequencing

In total, 12 mRNA samples were submitted to the Genomic platform INRAE center of Toulouse Occitanie, France, for further analyses, namely mRNA quality checking, cDNA preparation, and PE mRNA sequencing. The cDNA library was constructed and then sequenced using Illumina Inc to obtain a length of approximately 150 bp (Figure 1).

## 2.4. De Novo Assembly of Transcriptome

Quality control of the sequences generated by Illumina was assayed using the Fastqc v0.11.9 tool. Reads were used for de novo assembly using Trinity tool v2.14.0, with the default parameters such as K-mers equal to 25 and a minimum length equal to 200 bp [23]. After that, a second assembly was carried out using CD-HIT (http://weizhong-cluster.ucsd.edu/cd-hit/with 95% identity) (accessed on 10 June 2022) [24]. The quality of the secondary assembly was then evaluated using BUSCO tool version 5.4.3 [25]. In addition, the reads were mapped again to the transcripts generated by CD-Hit using the Bowtie2 tool
v2.4.5 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml) (accessed on 10 June 2022) (Figure 1). The transcripts generated in the study are available online: http://genoweb.toulouse.inra.fr/~mzouine/Hi_Incana/ (accessed on 10 June 2022).

## 2.5. Annotation and Function Profiling

Annotation of assembled transcripts was performed with BLASTX (E-Value $1e^{-10}$) on a local server using protein sequences of *Arabidopsis thaliana* downloaded from the ensemblgenomes.org database (http://ftp.ensemblgenomes.org/pub/plants/release-52/fasta/arabidopsis_thaliana/pep/Arabidopsis_thaliana.TAIR10.pep.all.fa.gz) (accessed on 10 June 2022). The BLASTX results (file represents each transcript ID with their TAIR ID) were combined with the Gene Ontology (GO) annotation file of *Arabidopsis thaliana* (https://www.arabidopsis.org/download_files/GO_and_PO_Annotations/Gene_Ontology_Annotations/geneassociation.tair.gz) (accessed on 10 June 2022) using Excel software 2021 for retrieving GO terms for each transcript (Figure 1). With these GO terms, all transcripts identified are classified into biological process (BP), molecular function (MF), and cellular component (CC) domains. Then, the ShinyGO v0.76.1 tool (http://bioinformatics.sdstate.edu/go/) (accessed on 10 June 2022) [26] and REVIGO tool (http://revigo.irb.hr/) (accessed on 10 June 2022) [27] were used to classify GO terms and to graphically represent the transcript functions.

## 2.6. Transcript Estimation and DEG Analysis

For the quantification of gene expression, we mapped the reads of each replicate against the transcriptome generated by CD-HIT using the galaxy server (https://usegalaxy.org) (accessed on 10 June 2022). In this analysis, the RSEM v1.3.3 algorithm [23] was used for transcripts estimation, then bowtie2 v2.4.5 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml) (accessed on 10 June 2022) was used as a tool for mapping the reads onto the transcriptome.

DEGs were inferred based on the normalized counts with R software by using the DESeq2 v3.15 package (https://bioconductor.org/packages/release/bioc/html/DESeq2.html) (accessed on 10 June 2022), and results were extracted using the restrictive method "Bonferroni" with a *p*-adj value <0.05 [28]. After that, the Venn-Diagram v1.7.3 algorithm was used to determine the numbers of DEGs shared between control and Pb-treated plants in roots and aboveground parts [29]. Then, the Gene Ontology and Enrichment Analysis tool available online (http://bioinformatics.sdstate.edu/go/) (accessed on 10 June 2022) and REVIGO tool [27] were used for functional annotation of expressed genes identified in roots and aboveground parts and to determine overrepresented GO terms across BP, MF, and MF domains. A *p*-value < 0.05 was used. The homologous genes found using BLASTX were used for pathway annotation through the KEGG database v103.1 (*Arabidopsis thaliana* used as a reference) (Figure 1). Roots and aboveground parts were harvested from six-week-old plants cultivated on Pb at 100 µM and the control. mRNAs were extracted (three replicates per condition) and sequenced using Illumina. Workflow of transcriptome analyses was performed. Whole transcriptome data were obtained and assembled. The assembled transcriptome was directly mapped onto the *A. thaliana* reference genome and then annotated and DEGs were measured to identify up- and downregulated genes under Pb exposure.

# 3. Results and Discussion

## 3.1. mRNA-seq and De Novo Transcriptome Assembly

mRNA samples extracted from roots and aboveground parts of *Hirschfeldia incana* plants were of high-quality RNA integrity number (superior to 7). The sequencing of mRNA generated millions of sequences with a length equal to 150 bp. An evaluation of the quality of all reads was performed using Fastqc. The results indicate that all reads (about 797.83 million reads) were of high quality with a phred-score ≥30. The percentage of GC was 44–47% for all the reads (Table S1). First, assembly analysis using the Trinity tool [23] generated 216,315 (211 Mb) transcripts with an average length equal to 979.57 bp and an
N50 = 1304 bp. The initial transcriptome assembled reported about 44.47 percent of GC content. The second transcriptome assembly using the CD-HIT tool [24], with 95 percent of identity, produced 77,491 (77 Mb) transcripts (Table 1). The number of transcripts found in this study was higher than that obtained in the transcriptome assembly in *Brassica juncea* with 53,669 transcripts [30] and lower than that mentioned in the reported transcriptome assembly in *Brassica napus* with 161,537 transcripts [31] (Table 1). The secondary transcripts of *Hirschfeldia incana* were characterized by N50 of 1330 bp and an average length of 959.83 bp. These findings were higher than those obtained in two species of the Brassicaceae family, namely *Brassica napus* (N50 = 1093 pb, average length = 693 pb) [31] and *Brassica juncea* (N50 = 1282 pb, average length = 953 pb) [30] but lower than those of *Brassica nigra* (N50 = 1482 pb, average length = 1173 pb) [32]. The secondary assembly of the *H. incana* transcriptome revealed 43.83% of GC content (Table 1). This percentage was slightly lower than GC content reported in *B. juncea* (44.38%) [30], and higher than GC content reported in *Brassica nigra* (37%) [32]. The reads were mapped again to the secondary assembly using the bowtie2 aligner tool to evaluate the consistency of the transcriptomes generated by CD-Hit. All reads were successfully mapped to the secondary assembly with a percentage higher than 93.5 percent (Table S1). The quality of the transcriptome assembly was verified using the BUSCO score algorithm. As shown in Table 1, a high percentage of BUSCO score was obtained; about 94.2% were complete, 2.8% were fragmented, and 3.0% were missing (Table S1). The percentage of the BUSCO reported in the present study was slightly lower than that reported in *Brassica nigra* [32].

<table><caption>**Table 1.** Transcriptome comparison of four Brassicaceae species including *H. incana*, *B. napus* [31], *B. juncea* [30], and *B. nigra* [32].</caption><thead><tr><th></th><th>H. incana</th><th>B. napus</th><th>B. nigra</th><th>B. juncea</th></tr></thead><tbody><tr><td colspan="5"><strong>Initial Transcripts (TRINITY)</strong></td></tr><tr><td>Total transcripts</td><td>216,315</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Total assembled bases</td><td>211,894,927</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Average length (pb)</td><td>979.57</td><td>-</td><td>-</td><td>-</td></tr><tr><td>N50</td><td>1304</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Percent GC</td><td>44.47</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="5"><strong>Secondary transcripts (CD-HIT)</strong></td></tr><tr><td>Total transcripts</td><td>77,491</td><td>161,537</td><td>NA</td><td>53,669</td></tr><tr><td>Total assembled bases</td><td>74,378,239</td><td>111,953,629</td><td>NA</td><td>51,151,545</td></tr><tr><td>Average length (pb)</td><td>959.83</td><td>693</td><td>1173</td><td>953</td></tr><tr><td>N50</td><td>1330</td><td>1093</td><td>1482</td><td>1282</td></tr><tr><td>Percent GC</td><td>43.83</td><td>NA</td><td>37</td><td>44.38</td></tr><tr><td>BUSCO %</td><td>94.2%</td><td>NA</td><td>97%</td><td>NA</td></tr><tr><td><strong>Number of genes</strong></td><td><strong>50,707</strong></td><td><strong>NA</strong></td><td><strong>56,331</strong></td><td><strong>NA</strong></td></tr></tbody></table>

The annotation of 77,491 secondary transcripts was performed using BLASTX on the protein database of *Arabidopsis thaliana*. About 72,306 annotated transcripts were obtained, with 11,862 duplicate transcripts removed and 60,444 unique transcripts (Table S3). GO annotation of *H. incana* transcripts was performed using the gene ontology annotation file of *Arabidopsis thaliana*. Afterward, GO annotations were classified into three essential domains, BP, CC, and MF (Figure 2; Table S3). Out of 60,444 annotated transcripts, 60,395 transcripts possess GO. About 13,046 own unique GO terms. In total, 7048 transcripts belong to the BP domain, with 1170 BP functions highly overrepresented (*p* < 0.05) (Table S3). The majority of them belong to the cellular process "GO:0009987", metabolic processes "GO:0008152", response to stimulus "GO:0050896", and response to stress "GO:0006950" (Figure 2; Table S3). Similar findings were reported in *Brassica napus* [31]. For MF, a total
of 4535 terms were found (Figure 2; Table S3). The most overrepresented activities were binding “GO:0005488” followed by catalytic activity “GO:0003824” (Figure 2; Table S3). The same activities were shown most represented in the MF category of *Brassica napus* [31]. However, 1463 GOs were found in the CC, most of them located in organelle “GO:0043226”, cytoplasm “GO:0005737”, and membrane “GO:0016020” (Figure 2; Table S3). Similar results were reported in *Brassica juncea* [30]. In addition, pathway analyses were performed using the KEGG database (*Arabidopsis thaliana* as the reference) to identify the pathways and understand the possible gene interactions. In total, 13,046 transcripts were assigned to 138 pathways. Among these pathways, metabolic pathways, biosynthesis of secondary metabolites, and plant hormone signal transduction were the top nine pathways represented by transcripts (Table S3). Our results indicate that the transcriptome generated in this study is successfully assembled and characterized by high quality, with relatively complete annotation information. The transcriptome assembly obtained provides helpful reference data and can be used to conduct subsequent functional genomics research in *Hirschfeldia incana*.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** The GO classification of *H. incana* transcripts in response to Pb stress.</figcaption></figure>
## 3.2. Differential Expressed Genes (DEGs) under Pb Exposure

For *H. incana*, the Venn diagram results indicate that most transcripts were shared between treated (Pb) and control (T) conditions. About 38,721 and 41,772 transcripts were shared between the two conditions in aboveground parts and roots, respectively (Figure 3A,B). Using the empirical criterion of ≥2-fold change for upregulated and ≤-2-fold change for downregulated and restrictive method "Bonferroni" with *p*-value < 0.05, a total of 577 and 270 genes were differentially expressed (DE) under Pb exposure in roots and aboveground parts, respectively (Figure 3C, Table S5). For roots, among 577 DEGs, 508 transcripts were upregulated, whereas 69 were downregulated under Pb stress (Figure 3C, Table S5). Without using the restrictive method, several studies detected a high number of genes in response to heavy-metal stress. For example, a total of 4614 DEGs were detected in roots of *Raphanus sativus* L., with 2154 upregulated genes and 2460 downregulated genes under Pb stress [33]. Another study reported 4682 genes expressed and 3599 inhibited in *B. juncea* root under Cd stress [34]. In our study, among the 270 genes, 244 were over-expressed, whereas 26 genes were inhibited in aboveground parts of *H. incana* in response to Pb stress (Figure 3C, Table S4). Other authors reported that under Pb stress, about 1641 over-expressed and 2884 inhibited genes in shoots of Pb hyper-tolerant species, namely *Fagopyrum tataricum* [3]. A study on *B. juncea* revealed 2021 upregulated and 992 downregulated genes in the shoots in response to Cd stress [34].

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Venn diagram and DEGs in roots and aboveground parts of *H. incana*: (A) Venn diagram for control (T) and treated (Pb) in aboveground parts; (B) Venn diagram for control (T) and treated (Pb) in roots; and (C) number of DEGs in different tissues.</figcaption></figure>

The GO enrichment analysis was performed and DEGs were annotated and classified into BP, CC, and MF domains (Table S5). For *H. incana* roots, the GO enrichment analysis of upregulated genes identified 432 CC, 221 MF, and 96 CC terms (Figure 4A, Table S5). Response to stress "GO:0006950", biosynthetic process "GO:0009058", and developmental process "GO:0032502" are the most represented GO terms in the BP domain (Figure 4A).

For the MF domain, catalytic activity "GO:0003824" and protein binding "GO:0005515" were the most represented (Figure 4A). In the CC domain, cytoplasm "GO:0005737", membrane "GO:0016020", plastid "GO:0009536", and cell wall "GO:0009505" were over-represented. Under heavy-metal stress, including Pb, the cell wall is the first barrier that protects the cell against Pb entry into the cytoplasm [35]. The cell wall contains several component such as polysaccharides and proteins that may bind Pb ions [36]. Such bind-
ing severely limits the Pb transport into the cell, thereby allowing cellular metabolism maintenance [37]. Our results showed that the CC group of genes is particularly enriched by Pb stress in roots of *H. incana* and contains genes expressed in the cell walls and cell membranes. Additionally, genes associated with the biosynthesis of polysaccharides such as cellulose were shown to be significantly upregulated (LogFC = 3.31) in response to Pb stress, which may result in considerable cell wall thickening [3].

A

<figure><img src="image_8.png" /></figure>

B

<figure><img src="image_9.png" /><figcaption>**Figure 4.** Top 12 GO classification of upregulated genes in roots (A) and aboveground parts (B) of *H. incana* under Pb exposure.</figcaption></figure>

Transcripts identified in *H. incana* roots exposed to Pb were related to the catalytic activity, protein binding, and nitrogen compound metabolic process. Similarly, in *R. sativus* L. roots, different BP and MF were reported under Pb stress [33]. These overrepresented GO terms were also related to catalytic activity and explained by the adverse effects of heavy metals in plants resulting in the inhibition of biomass production and plant growth [38].

For downregulated transcripts, the most represented GO, identified in roots of *H. incana* in response to Pb stress, are related to carbohydrate metabolic process "GO:0005975", defense response "GO:0006952", and proteolysis "GO:0006508", in BP (Table S5). Protein binding and DNA binding are most represented in the MF domain (Table S4). For CC, the overrepresented downregulated transcripts were located in the intracellular membrane-bounded organelle, cytoplasm, and nucleus (Table S5). In *R. sativus* roots, cell-wall organization and catalytic activity were the most represented in BP and MF domains, respectively, in downregulated genes under Pb stress [33]. The KEGG pathway enrichment analysis of upregulated transcripts in *H. incana* L. roots identified 67 enriched pathways in response to Pb stress (Table S5), 10 of which were overrepresented (Figure 5A, Table S5). Metabolic, biosynthesis of secondary metabolites, and carbon metabolism were the most significantly enriched pathways responding to Pb stress. For downregulated transcripts, 17 pathways were reduced in response to Pb stress (Figure 5A, Table S5), suggesting that pathways could differ between species in response to Pb stress [33].

In the aboveground parts, GO enrichment analysis of upregulated genes identified 382 BP terms, 245 MF terms, and 64 CC (Figure 4B, Table S5). For the upregulated transcripts, developmental process "GO:0032502", defense response "GO:0006952", and small
molecule metabolic process "GO:0044281" were the most represented in the BP domain (Figure 4B). In the MF domain, protein binding, cation binding, hydrolase activity, and oxidoreductase activity and those related to the CC domain include intracellular anatomical structure, cytoplasm "GO:0005737", membrane "GO:0016020", cell wall, chloroplast "GO:0009507", and Golgi apparatus GO:0005794 (Figure 4B; Table S5) were the most represented. These activities and functions were reported in DEGs in *F. tataricum* shoots in response to Pb exposure [3]. For downregulated genes detected in *H. incana* aboveground parts in response to Pb stress, 42 GO overrepresented in the BP domain include root morphogenesis "GO:0010015", cellular protein localization "GO:0034613", and response to light stimulus "GO:0009416" (Table S5). For MF, 25 GO including protein binding "GO:0005515", and nucleic acid binding "GO:0003676" (Table S5), and 24 GO related to CC include cytoplasm "GO:0005737", and nucleus "GO:0005634" (Table S5). However, 100% of the 10 most up- and down-expressed transcripts in roots and aboveground parts are associated with a GO term (Table 2). The KEGG analysis showed 57 different pathways affected by Pb stress in aboveground parts. The DEGs analysis revealed that the most upregulated transcripts were associated with different pathways, including metabolic, glucosinolate biosynthesis, and alpha-linolenic acid (Figure 5B; Table S5). Different pathways are identified in shoots of *F. tataricum* under Pb exposure [3]. For the downregulated genes, many pathways, such as plant autophagy and metabolic are overrepresented in aboveground parts of *H. incana* L. under Pb stress (Figure 5B; Table S5).

A

<figure><img src="image_10.png" /></figure>

B

<figure><img src="image_11.png" /><figcaption>**Figure 5.** Top pathways of DEGs in roots (A) and aboveground parts (B) of *H. incana* in response to Pb exposure.</figcaption></figure>

Previous studies on *H. incana* response to Pb reported that this perennial plant tolerates and accumulates a high level of Pb in its aboveground parts [8,10,20]. These results suggest that *H. incana* has developed adaptive mechanisms to enhance the capacity to survive, tolerate, and accumulate Pb and are certainly implicated in the activation of many essential transporters. A multitude of transporter families, including OCT, CNGC, and ABC which have been characterized as playing a crucial role in heavy-metal uptake, transport, distribution, and plant tolerance [39], were identified as upregulated in the present study. γ-aminobutyric acid (GABA) transporter was upregulated in *H. incana* aboveground parts in response to Pb stress. GABA is implicated in several physiological activities and contributes to the tolerance of plants to different stresses, including heavy
metals [40]. Previous studies showed that GABA was upregulated during Zn and Cr stress in *Nicotiana tabaccum* L. [41] and *Brassica juncea* [42]. AVP1, a proton pump located in the vacuolar membrane that enhances cadmium tolerance and accumulation in tobacco plants [43], was shown upregulated in aboveground parts of *H. incana* under Pb stress (Table S5). The ATP-Mg/Pi transporter (APC3: mitochondrial substrate carrier family protein) was shown upregulated in aboveground parts of *H. incana* during Pb stress (Table S5). AtAPC1 and AtAPC3 are known to transport phosphate compounds such as AMP, ADP, and ATP [44]. Abiotic stress, such as heavy-metal stress, causes increased energy demand [45], which requires high consumption of ATP and ADP [45]. The high energy consumption may be related to the expression of the APC3 gene under Pb stress.

<table><caption>**Table 2.** Example of the most upregulated DEGs (A) and downregulated DEGs (B) in roots and aboveground parts of *H. incana* in response to Pb stress. Bold numbers indicate the standard error for the Log foldchange.</caption><thead><tr><th colspan="6">(A)</th></tr><tr><th>Transcript ID</th><th>Transcript Name</th><th>Log2FC</th><th>p-Value</th><th>N° of GO Terms</th><th></th></tr></thead><tbody><tr><td colspan="6">Up-Regulated Transcripts in aboveground Parts of H. incana</td></tr><tr><td>HiInc_DN6940_c0_g2</td><td>Noc2p family</td><td>10.56</td><td>1.94</td><td>5.43 × 10<sup>-8</sup></td><td>8</td></tr><tr><td>HiInc_DN64666_c0_g1</td><td>RNA-directed DNA polymerase (Reverse transcriptase)-related family protein</td><td>8.55</td><td>1.47</td><td>6.97 × 10<sup>-9</sup></td><td>3</td></tr><tr><td>HiInc_DN54630_c0_g1</td><td>Translation initiation factor IF3-4</td><td>7.31</td><td>1.40</td><td>2.02 × 10<sup>-7</sup></td><td>8</td></tr><tr><td>HiInc_DN38797_c0_g1</td><td>Peroxidase</td><td>6.07</td><td>1.84</td><td>1.57 × 10<sup>-11</sup></td><td>4</td></tr><tr><td>HiInc_DN8233_c0_g1</td><td>Inorganic pyrophosphatase 2</td><td>5.65</td><td>0.86</td><td>6.74 × 10<sup>-11</sup></td><td>6</td></tr><tr><td>HiInc_DN250_c0_g1</td><td>Syringolide-induced protein</td><td>5.28</td><td>0.88</td><td>2.89 × 10<sup>-9</sup></td><td>7</td></tr><tr><td>HiInc_DN35271_c0_g1</td><td>Beta-glucosidase 32</td><td>5.13</td><td>1.64</td><td>3.11 × 10<sup>-12</sup></td><td>6</td></tr><tr><td>HiInc_DN5802_c0_g2</td><td>Monogalactosyldiacylglycerol synthase 2</td><td>5.10</td><td>0.57</td><td>4.03 × 10<sup>-19</sup></td><td>9</td></tr><tr><td>HiInc_DN12789_c0_g1</td><td>REF/SRPP-like protein</td><td>5.02</td><td>0.69</td><td>3.25 × 10<sup>-13</sup></td><td>7</td></tr><tr><td>HiInc_DN1096_c0_g4</td><td>Salicylate/benzoate carboxyl methyltransferase</td><td>4.89</td><td>0.92</td><td>1.36 × 10<sup>-7</sup></td><td>5</td></tr><tr><td colspan="6">Up-regulated transcripts in roots of H. incana</td></tr><tr><td>HiInc_DN6451_c0_g2</td><td>Endochitinase</td><td>11.68</td><td>1.75</td><td>2.60 × 10<sup>-11</sup></td><td>7</td></tr><tr><td>HiInc_DN64488_c0_g1</td><td>ATPase 9</td><td>10.29</td><td>1.96</td><td>1.52 × 10<sup>-7</sup></td><td>4</td></tr><tr><td>HiInc_DN5848_c0_g1</td><td>Accelerated cell death 11</td><td>9.89</td><td>1.96</td><td>5.01 × 10<sup>-7</sup></td><td>7</td></tr><tr><td>HiInc_DN27683_c0_g1</td><td>Putative endonuclease or glycosyl hydrolase</td><td>9.74</td><td>1.32</td><td>2.38 × 10<sup>-13</sup></td><td>2</td></tr><tr><td>HiInc_DN6940_c0_g1</td><td>Noc2p family</td><td>9.47</td><td>1.49</td><td>2.56 × 10<sup>-10</sup></td><td>4</td></tr><tr><td>HiInc_DN10792_c0_g1</td><td>Nucleic acid-binding</td><td>8.82</td><td>1.48</td><td>3.06 × 10<sup>-9</sup></td><td>3</td></tr><tr><td>HiInc_DN62699_c2_g1</td><td>BED zinc finger and hAT dimerization domain-containing protein DAYSLEEPER</td><td>8.70</td><td>1.63</td><td>1.03 × 10<sup>-7</sup></td><td>1</td></tr><tr><td>HiInc_DN63847_c0_g1</td><td>Alpha carbonic anhydrase 2</td><td>8.36</td><td>1.29</td><td>1.20 × 10<sup>-10</sup></td><td>4</td></tr><tr><td>HiInc_DN13103_c0_g1</td><td>60S ribosomal protein L10a-2</td><td>7.81</td><td>1.33</td><td>5.24 × 10<sup>-9</sup></td><td>2</td></tr><tr><td>HiInc_DN8575_c1_g1</td><td>Disease resistance protein (TIR-NBS class)</td><td>7.77</td><td>1.39</td><td>2.71 × 10<sup>-8</sup></td><td>13</td></tr><tr><td colspan="6">(B)</td></tr><tr><th>Transcript ID</th><th>Transcript name</th><th>Log2FC</th><th>p-value</th><th>N° of GO terms</th><th></th></tr><tr><td colspan="6">Down-regulated transcripts in aboveground parts of H. incana</td></tr><tr><td>HiInc_DN38456_c0_g1</td><td>transmembrane protein, putative (DGR2)</td><td>-2.10</td><td>0.25</td><td>1.20 × 10<sup>-16</sup></td><td>1</td></tr><tr><td>HiInc_DN19265_c0_g1</td><td>Chaperone DnaJ-domain superfamily protein</td><td>-2.17</td><td>0.43</td><td>7.06 × 10<sup>-7</sup></td><td>5</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Transcript ID</th><th>Transcript name</th><th>Log2FC</th><th>p-value</th><th>N° of GO terms</th></tr></thead><tbody><tr><td>HiInc_DN12488_c0_g1</td><td>Thioredoxin-like protein CXXS1</td><td>-2.19</td><td>0.34</td><td>3.79 × 10<sup>-6</sup></td><td>2</td></tr><tr><td>HiInc_DN54188_c0_g1</td><td>Tetratricopeptide repeat (TPR)-like</td><td>-2.24</td><td>0.45</td><td>1.02 × 10<sup>-6</sup></td><td>2</td></tr><tr><td>HiInc_DN9080_c0_g1</td><td>Tubulin alpha chain</td><td>-2.29</td><td>0.39</td><td>4.30 × 10<sup>-9</sup></td><td>15</td></tr><tr><td>HiInc_DN33913_c0_g1</td><td>Probable aquaporin PIP2-8</td><td>-2.37</td><td>0.42</td><td>2.20 × 10<sup>-8</sup></td><td>3</td></tr><tr><td>HiInc_DN7675_c0_g1</td><td>Gamma vacuolar processing enzyme</td><td>-2.51</td><td>0.35</td><td>2.21 × 10<sup>-12</sup></td><td>12</td></tr><tr><td>HiInc_DN10024_c0_g1</td><td>RNA-binding (RRM/RBD/RNP motifs)</td><td>-2.75</td><td>0.55</td><td>8.13 × 10<sup>-7</sup></td><td>6</td></tr><tr><td>HiInc_DN18907_c0_g1</td><td>RmlC-like cupins superfamily protein</td><td>-2.79</td><td>0.52</td><td>7.56 × 10<sup>-8</sup></td><td>2</td></tr><tr><td colspan="6"><strong>Down-regulated transcripts in roots of H. incana</strong></td></tr><tr><td>HiInc_DN20249_c0_g1</td><td>Ribosomal protein S3Ae</td><td>-3.27</td><td>0.53</td><td>8.09 × 10<sup>-10</sup></td><td>8</td></tr><tr><td>HiInc_DN1121_c0_g2</td><td>RmlC-like cupins superfamily protein</td><td>-2.02</td><td>0.22</td><td>8.16 × 10<sup>-19</sup></td><td>3</td></tr><tr><td>HiInc_DN2301_c0_g1</td><td>P-loop containing nucleoside triphosphate hydrolases superfamily protein</td><td>-2.06</td><td>0.33</td><td>5.82 × 10<sup>-10</sup></td><td>4</td></tr><tr><td>HiInc_DN315_c5_g1</td><td>Receptor-like protein 30</td><td>-2.07</td><td>0.39</td><td>1.50 × 10<sup>-7</sup></td><td>6</td></tr><tr><td>HiInc_DN6047_c0_g1</td><td>Protein kinase superfamily protein</td><td>-2.10</td><td>0.39</td><td>1.22 × 10<sup>-7</sup></td><td>2</td></tr><tr><td>HiInc_DN10269_c1_g2</td><td>UDP-glycosyltransferase 76B1</td><td>-2.42</td><td>0.49</td><td>1.07 × 10<sup>-6</sup></td><td>8</td></tr><tr><td>HiInc_DN38838_c0_g1</td><td>60S acidic ribosomal protein P2-3</td><td>-2.45</td><td>0.33</td><td>2.33 × 10<sup>-13</sup></td><td>2</td></tr><tr><td>HiInc_DN16964_c0_g1</td><td>Probable LRR receptor-like serine/threonine-protein kinase</td><td>-2.48</td><td>0.34</td><td>8.71 × 10<sup>-10</sup></td><td>2</td></tr><tr><td>HiInc_DN22651_c0_g3</td><td>Glutathione S-transferase F3</td><td>-2.51</td><td>0.26</td><td>7.18 × 10<sup>-7</sup></td><td>5</td></tr><tr><td>HiInc_DN24519_c1_g1</td><td>Probable LRR receptor-like serine/threonine-protein kinase</td><td>-2.51</td><td>0.50</td><td>5.27 × 10<sup>-7</sup></td><td>7</td></tr><tr><td>HiInc_DN13802_c0_g1</td><td>Tricyclene synthase</td><td>-2.64</td><td>0.26</td><td>5.76 × 10<sup>-24</sup></td><td>8</td></tr></tbody></table>

Several studies reported that heavy metals affect the absorption of Mg and the photosynthesis process. For example, Cr reduces the absorption of Mg resulting in lower chlorophyll content [46]. Pb causes an increase in the activity of some enzymes, specifically chlorophyllase. Additionally, Pb could modify and change the activity of photosynthesis [47]. In the present study, the treatment with 100 µM of Pb decreases the expression level of the Mg ion transporter and can directly reduce the photosynthesis activity under Pb stress. In this study, several magnesium ion transporters of the *MRS2* family were shown regulated by Pb. *MRS2-3*, *MRS2-7*, and *AtMRS2-11*, essential in transporting Mg in the chloroplast membrane system [48], were shown downregulated in aboveground parts of *H. incana* under Pb exposure

In *H. incana* roots, many transporters are significantly upregulated in response to Pb, such as *CNGC5*. *CNGC* family transporters are channels whereby metal ions can enter the cells and are involved in tolerance to heavy-metal stress [49]. In *Arabidopsis thaliana*, many members of the *CNGC* family such as *AtCNGC11* and *AtCNGC15* have been recorded in plant tolerance and uptake of Pb²⁺ and Cd²⁺ [50]. In the current study, the *HAK5* gene involved in K⁺ uptake [51] is shown upregulated in the root of *H. incana* in response to Pb stress. Metal toxicity induces oxidative stress, leading to ROS production and modification of several physiological activities, including degradation of several enzymes, amino acids, and proteins [52]. Potassium activates the antioxidant defense in plants and therefore increases oxidative stress tolerance by minimizing NADPH's activity and controlling photosynthetic activity, which helps to limit ROS production [53]. In the present study, upregulation of the *HAK5* gene in *H. incana* roots suggests the implication of potassium transporters in the reduction of the damage caused by Pb. To sum up, several genes coding for *ABC*, *OCT*, *MRS*, and *CNGC* transporter families were upregulated and may play an
essential role in *H. incana* Pb tolerance. Further functional analyses are necessary to validate their specific implication in response to Pb exposure.

Several genes involved in regulation, biosynthesis, and transport of glucosinolates such as SOT18 (sulfotransferase), GTR1 (major facilitator superfamily protein), APK2 (adenylyl-sulfate kinase 2), JAL23, BASS5, BGLU32 (beta-glucosidase 32), and IGMT5 in aboveground parts and SDI1 (sulphur deficiency-induced 1) in roots were upregulated in *H. incana* in response to Pb stress. Our finding is in accordance with several reports mentioning that under Zn, Cd, and Pb the biosynthesis of glucosinolates was higher in other *Brassica* species [54,55]. These authors suggest that glucosinolates play an important role in heavy-metal tolerance in different *Brassica* plant species [54,55].

Under Pb exposure, several transcription factors (TFs) were shown upregulated in aboveground parts such as WRKY40, MYB95, and MYB15 and in roots of *H. incana* such as ERF021 and MYB108. Transcription factors are implicated in the control of different transcriptional activities under heavy-metal stress [56,57]. Several studies reported an over-expression of MYB and WARKY families in different plants species in response to heavy-metal exposure [57,58]. MYB4 plays an essential role in controlling gene expres-sion in different pathways, specifically the phenylpropanoid pathway. In the present study, C4H (cinnamate-4-hydroxylase) and CCR2 (Cinnamoyl-CoA reductase 2), impli-cated in phenylpropanoid biosynthetic activities, were upregulated in aboveground parts and roots of *H. incana*, respectively, under Pb exposure. Some authors reported that the phenylpropanoid pathway is activated under heavy-metal exposure and involved in the protection of plants against ROS [59]. This type of TF is implicated directly in the control of glucosinolate production that is increased in *N. caerulescens* shoots after exposure to Cd [60]. In addition, several genes implicated in jasmonates (JAs) biosynthesis are shown to be significantly upregulated in aboveground parts of *H. incana* in response to Pb stress. Various TFs increase the concentrations of JAs and jasmonyl isoleucine by controlling the activation of some essential genes implicated directly in JAs biosynthesis such as LOX, AOS, and AOC in response to heavy-metal stress [61]. JAs are also essential for the induction of ERF1 in response to HMs stresses [62]. ERFs were reported as playing an essential role in the response of plants to Cd and Pb stress by regulating several genes' expression [33].

### 3.3. Jasmonates (JAs) and Plant Response to Pb—A Case Study

Jasmonates (JAs) are signaling molecules involved in plant development and plant responses to biotic and abiotic stresses [63]. In the present study, whole transcriptomic analysis of *Hirschfeldia incana* showed that under Pb stress several upregulated genes in roots and aboveground parts are involved in JAs biosynthesis (Figure 6). Seven genes were characterized as DEGs in aboveground parts, namely LOX2, in chloroplast, AOS in thylakoid and plastid, AOC2 and AOC4 in plastid, ACX1 and KAT5 in peroxisome, and the JAZ10 gene located in the nucleus (Table S4). In *H. incana* roots, the JAZ10 gene was shown upregulated under Pb exposure. Several studies reported that JAs synthesis and pathways are affected by HM stress. For example, Cd activates the expression of genes involved in JAs synthesis, thus, increasing the concentrations of JAs in the root tissues of *Arabidopsis thaliana* [64]. The authors mentioned that AOS and AOC genes were upregulated in *Arabidopsis* under Pb exposure [65]. A study by [66] reported that JAs are synthesized through several steps from lipoxygenases, followed by oxygenation of α-linolenic acid. According to [67,68], the JA biosynthesis is controlled by two crucial reduction and oxidation steps (Figure 6). In this study, DEGs analysis reveals that the JAs pathway is activated by Pb stress. This finding suggests that JAs play an essential role in *H. incana* response to Pb stress.

Alpha-linolenic acid is an essential element for the JA biosynthesis [66]. In our study, the FAD7 (fatty acid desaturase 7) gene was significantly upregulated with logFC = 2.40 in response to Pb stress (Table S5; Figure 6). The FAD gene is implicated directly in the biosynthesis of unsaturated fatty acid such as α-linolenic acid content and adversely affects JAs in response to heavy-metal stress [69]. Previous studies mentioned that fatty acid
desaturase genes are essential for stabilizing the plant membrane under abiotic stress, including heavy metals, by producing high amounts of fatty acids [70,71]. Fatty acids are essential components of cellular membranes and can play an important role in structural barriers for minimizing the effects of Pb in plants [72]. They contribute to inducible stress tolerance through controlling membrane fluidity [73]. Alpha-linoleic acid can also be used as a substrate for JA biosynthesis. Alpha-linolenic acid and JAs then emerge as a strategy to withstand Pb toxicity in *H. incana*.

<figure><img src="image_12.png" /><figcaption>**Figure 6.** Expression level of genes implicated in the JAs biosynthesis pathways in *H. incana* under Pb exposure. The red color represents the genes selected using the restrictive method "Bonferroni", and the light red represent the genes selected without using the restrictive method.</figcaption></figure>

**4. Conclusions**

This work presents the first evidence on the RNA-seq, de novo assembly, and functional annotation of the Greek mustard (*Hirschfeldia incana* L.) transcriptome, and the differential expression genes (DEGs) in roots and aboveground parts in response to Pb stress. The de novo assembly produced 77,491 transcripts, creating a reference transcriptome.
tome with an N50 of 1330 bp and 43.83 percent of GC. We report more than 60,444 annotated transcripts in *H. incana*. In total, 577 and 270 DEGs were identified in response to Pb stress in *H. incana* roots and aboveground parts, respectively. The GO term enrichment revealed that in *H. incana*, most identified genes involved in biosynthetic processes were significantly upregulated under Pb stress. Whereas genes involved in carbohydrate metabolic process, and defense response were significantly downregulated in roots. In aboveground parts, the majority of genes significantly upregulated are involved in lipid metabolic process, defense response, and response to jasmonic acid. In addition, several genes implicated in jasmonates biosynthesis were significantly upregulated under Pb exposure, suggesting that *H. incana* uses the JAs pathway as a strategy to alleviate Pb stress. The dataset generated in this study also contributes to the molecular resources of *Hirschfeldia incana*.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cimb44100318/s1, Table S1: RNA-Seq data summary; Table S2: Results of BlastX; Table S3: Annotations of total transcripts; Table S4: Gene expressions sorted by normal method; Table S5: Gene expressions sorted by restricted method; Table S6: Pathways identified using KEGG.

**Author Contributions:** Methodology, investigation, formal analysis, validation, visualization, writing—original draft, S.E.H.; funding acquisition, conceptualization, supervision, writing—original draft, M.F.; validation, visualization, M.Z.; project administration, funding acquisition, conceptualization, supervision, writing—original draft, A.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** Published with the support of "IRD: JEAI—PhytoMet" 2018–2020; IRD Joint Laboratory "LMI AMIR" 2020–2025"; Programme Prioritaire de recherche "PPR2-BIOMIVER-1108/16"-Ministère de l'Enseignement Supérieur, de la Recherche Scientifique et de la Formation des Cadres–Morocco and the COST Action "RoxyCOST (CA18210)" 2019–2023.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The transcripts generated in the study are available online: http://genoweb.toulouse.inra.fr/~mzouine/Hi_Incana/.

**Acknowledgments:** The authors would like to thank "IRD: JEAI—PhytoMet" 2018–2020; IRD Joint Laboratory "LMI AMIR" 2020–2025"; Programme Prioritaire de recherche "PPR2-BIOMIVER"-Ministère de l'Enseignement Supérieur, de la Recherche Scientifique et de la Formation des Cadres–Morocco and the COST Action "RoxyCOST (CA18210)" 2019–2023 for their support. We thank reviewers and editors for their time and fruitful comments.

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

**Abbreviations**

<table><tr><td>ACX1</td><td>Acyl-CoA oxidase 1;</td></tr><tr><td>AOC4</td><td>allene oxide cyclase 4;</td></tr><tr><td>AOS</td><td>allene oxide synthase;</td></tr><tr><td>BUSCO</td><td>Benchmarking Universal Single-Copy Orthologs;</td></tr><tr><td>CD-HIT</td><td>Cluster Database at High Identity with Tolerance;</td></tr><tr><td>DEGs</td><td>differentially gene expressions;</td></tr><tr><td>DE</td><td>differential expressed;</td></tr><tr><td>GO</td><td>Gene Ontology;</td></tr><tr><td>GTH</td><td>glutathione;</td></tr><tr><td>JAZ10</td><td>JASMONATE ZIM-domain protein 10;</td></tr><tr><td>JAs</td><td>jasmonates;</td></tr><tr><td>KAT5</td><td>3-keto-acyl-CoA thiolase 2;</td></tr><tr><td>KEGG</td><td>Kyoto Encyclopedia of Genes and Genomes;</td></tr><tr><td>LOX2</td><td>lipoxygenase;</td></tr></table>
<table><tr><td>HMs</td><td>heavy metals;</td></tr><tr><td>MTs</td><td>metallothioneins;</td></tr><tr><td>PCs</td><td>Phytochelatin;</td></tr><tr><td>REVIGO</td><td>reduce and visualize gene ontology;</td></tr><tr><td>RIN</td><td>RNA Integrity Number;</td></tr><tr><td>ROS</td><td>reactive oxygen species;</td></tr><tr><td>RSEM</td><td>RNA-Seq by Expectation-Maximization;</td></tr><tr><td>TFs</td><td>transcription factors;</td></tr><tr><td>13(S)-HPOT</td><td>(9Z,11E,15Z)-(13S)-13-Hydroperoxyoctadeca-9,11,15-trienoic acid;</td></tr><tr><td>12,13(S)-EOT</td><td>(9Z,15Z)-(13S)-12,13-Epoxyoctadeca-9,11,15-trienoic acid;</td></tr><tr><td>12-OPDA</td><td>(15Z)-12-Oxophyto-10,15-dienoic acid;</td></tr><tr><td>OPC8</td><td>8-[(1R,2R)-3-Oxo-2-(Z)-pent-2-enylcyclopentyl] octanoate</td></tr></table>

## References

1. Tangahu, B.V.; Abdullah, S.R.S.; Basri, H.; Idris, M.; Anuar, N.; Mukhlisin, M. A review on heavy metals (As, Pb, and Hg) uptake by plants through phytoremediation. *Int. J. Chem. Eng.* **2011**, *2011*, 939161. [CrossRef]

2. Fahr, M.; Laplaze, L.; Bendaou, N.; Hocher, V.; El Mzibri, M.; Bogusz, D.; Smouni, A. Effect of lead on root growth. *Front. Plant Sci.* **2013**, *4*, 175. [CrossRef] [PubMed]

3. Wang, L.; Zheng, B.; Yuan, Y.; Xu, Q.; Chen, P. Transcriptome profiling of *Fagopyrum tataricum* leaves in response to lead stress. *BMC Plant Biol.* **2020**, *20*, 54. [CrossRef] [PubMed]

4. Mourato, M.; Reis, R.; Martins, L.L. Characterization of plant antioxidative system in response to abiotic stresses: A focus on heavy metal toxicity. *Adv. Sel. Plant Physiol. Asp.* **2012**, *12*, 23–44.

5. John, R.; Ahmad, P.; Gadgil, K.; Sharma, S. Cadmium and lead-induced changes in lipid peroxidation, antioxidative enzymes and metal accumulation in *Brassica juncea* L. at three different growth stages. *Arch. Agron. Soil Sci.* **2009**, *55*, 395–405. [CrossRef]

6. Mahdavian, K.; Ghaderian, S.M.; Schat, H. Pb accumulation, Pb tolerance, antioxidants, thiols, and organic acids in metallicolous and non-metallicolous *Peganum harmala* L. under Pb exposure. *Environ. Exp. Bot.* **2016**, *126*, 21–31. [CrossRef]

7. Mohtadi, A.; Ghaderian, S.M.; Schat, H. A comparison of lead accumulation and tolerance among heavy metal hyperaccumulating and non-hyperaccumulating metallophytes. *Plant Soil* **2012**, *352*, 267–276. [CrossRef]

8. Auguy, F.; Fahr, M.; Moulin, P.; Brugel, A.; Laplaze, L.; Mzibri, M.E.; Filali-Maltouf, A.; Doumas, P.; Smouni, A. Lead tolerance and accumulation in *Hirschfeldia incana*, a Mediterranean Brassicaceae from metalliferous mine spoils. *PLoS ONE* **2013**, *8*, e61932. [CrossRef]

9. Siemens, J.; González, M.C.; Wolf, S.; Hofmann, C.; Greiner, S.; Du, Y.; Rausch, T.; Roitsch, T.; Ludwig-Müller, J. Extracellular invertase is involved in the regulation of clubroot disease in *Arabidopsis thaliana*. *Mol. Plant Pathol.* **2011**, *12*, 247–262. [CrossRef]

10. Fahr, M.; Laplaze, L.; El Mzibri, M.; Doumas, P.; Bendaou, N.; Hocher, V.; Bogusz, D.; Smouni, A. Assessment of lead tolerance and accumulation in metallicolous and non-metallicolous populations of *Hirschfeldia incana*. *Environ. Exp. Bot.* **2015**, *109*, 186–192. [CrossRef]

11. Pourrut, B.; Shahid, M.; Douay, F.; Dumat, C.; Pinelli, E. Molecular mechanisms involved in lead uptake, toxicity and detoxification in higher plants. In *Heavy Metal Stress in Plants*; Springer: Berlin/Heidelberg, Germany, 2013; pp. 121–147.

12. Arazi, T.; Sunkar, R.; Kaplan, B.; Fromm, H. A tobacco plasma membrane calmodulin-binding transporter confers $Ni^{2+}$ tolerance and $Pb^{2+}$ hypersensitivity in transgenic plants. *Plant J.* **1999**, *20*, 171–182. [CrossRef] [PubMed]

13. Gravot, A.; Lieutaud, A.; Verret, F.; Auroy, P.; Vavasseur, A.; Richaud, P. *AtHMA3*, a plant *P1B-ATPase*, functions as a Cd/Pb transporter in yeast. *FEBS Lett.* **2004**, *561*, 22–28. [CrossRef]

14. Morel, M.; Crouzet, J.; Gravot, A.; Auroy, P.; Leonhardt, N.; Vavasseur, A.; Richaud, P. *AtHMA3*, a *P1B-ATPase* allowing Cd/Zn/co/Pb vacuolar storage in *Arabidopsis*. *Plant Physiol.* **2009**, *149*, 894–904. [CrossRef] [PubMed]

15. Lee, M.; Lee, K.; Lee, J.; Noh, E.W.; Lee, Y. *AtPDR12* contributes to lead resistance in *Arabidopsis*. *Plant Physiol.* **2005**, *138*, 827–836. [CrossRef] [PubMed]

16. Kim, D.-Y.; Bovet, L.; Kushnir, S.; Noh, E.W.; Martinoia, E.; Lee, Y. *AtATM3* is involved in heavy metal resistance in *Arabidopsis*. *Plant Physiol.* **2006**, *140*, 922–932. [CrossRef]

17. Kim, D.Y.; Bovet, L.; Maeshima, M.; Martinoia, E.; Lee, Y. The ABC transporter *AtPDR8* is a cadmium extrusion pump conferring heavy metal resistance. *Plant J.* **2007**, *50*, 207–218. [CrossRef]

18. Cobbett, C.; Goldsbrough, P. Phytochelatins and metallothioneins: Roles in heavy metal detoxification and homeostasis. *Annu. Rev. Plant Biol.* **2002**, *53*, 159–182. [CrossRef]

19. Auguy, F.; Fahr, M.; Moulin, P.; El Mzibri, M.; Smouni, A.; Filali-Maltouf, A.; Béna, G.; Doumas, P. Transcriptome changes in *Hirschfeldia incana* in response to lead exposure. *Front. Plant Sci.* **2016**, *6*, 1231. [CrossRef]

20. Evangelistella, C.; Valentini, A.; Ludovisi, R.; Firrincieli, A.; Fabbrini, F.; Scalabrin, S.; Cattonaro, F.; Morgante, M.; Mugnozza, G.S.; Keurentjes, J.J. De novo assembly, functional annotation, and analysis of the giant reed (*Arundo donax* L.) leaf transcriptome provide tools for the development of a biofuel feedstock. *Biotechnol. Biofuels* **2017**, *10*, 138. [CrossRef]
21. Hasnaoui, S.E.; Fahr, M.; Keller, C.; Levard, C.; Angeletti, B.; Chaurand, P.; Triqui, Z.E.A.; Guedira, A.; Rhazi, L.; Colin, F. Screening of native plants growing on a Pb/Zn mining area in eastern Morocco: Perspectives for phytoremediation. *Plants* **2020**, *9*, 1458. [CrossRef]

22. Broughton, W.; Dilworth, M. Control of leghaemoglobin synthesis in snake beans. *Biochem. J.* **1971**, *125*, 1075–1080. [CrossRef] [PubMed]

23. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. *Nat. Protoc.* **2013**, *8*, 1494–1512. [CrossRef] [PubMed]

24. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. *Bioinformatics* **2012**, *28*, 3150–3152. [CrossRef] [PubMed]

25. Waterhouse, R.M.; Seppey, M.; Simão, F.A.; Manni, M.; Ioannidis, P.; Klioutchnikov, G.; Kriventseva, E.V.; Zdobnov, E.M. BUSCO applications from quality assessments to gene prediction and phylogenomics. *Mol. Biol. Evol.* **2018**, *35*, 543–548. [CrossRef] [PubMed]

26. Ge, S.X.; Jung, D.; Yao, R. ShinyGO: A graphical gene-set enrichment tool for animals and plants. *Bioinformatics* **2020**, *36*, 2628–2629. [CrossRef] [PubMed]

27. Supek, F.; Bošnjak, M.; Škunca, N.; Šmuc, T. REVIGO summarizes and visualizes long lists of gene ontology terms. *PLoS ONE* **2011**, *6*, e21800. [CrossRef]

28. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biol.* **2014**, *15*, 550. [CrossRef]

29. Chen, H.; Boutros, P.C. VennDiagram: A package for the generation of highly-customizable Venn and Euler diagrams in R. *BMC Bioinform.* **2011**, *12*, 35. [CrossRef]

30. Sharma, R.; Mishra, M.; Gupta, B.; Parsania, C.; Singla-Pareek, S.L.; Pareek, A. De novo assembly and characterization of stress transcriptome in a salinity-tolerant variety CS52 of *Brassica juncea*. *PLoS ONE* **2015**, *10*, e0126783. [CrossRef]

31. Yong, H.-Y.; Zou, Z.; Kok, E.-P.; Kwan, B.-H.; Chow, K.; Nasu, S.; Nanzyo, M.; Kitashiba, H.; Nishio, T. Comparative transcriptome analysis of leaves and roots in response to sudden increase in salinity in *Brassica napus* by RNA-seq. *BioMed Res. Int.* **2014**, *2014*, 395–467. [CrossRef]

32. Seol, Y.-J.; Kim, K.; Kang, S.-H.; Perumal, S.; Lee, J.; Kim, C.-K. The complete chloroplast genome of two *Brassica* species, *Brassica nigra* and *B. oleracea*. *Mitochondrial DNA Part A* **2017**, *28*, 167–168. [CrossRef] [PubMed]

33. Wang, Y.; Xu, L.; Chen, Y.; Shen, H.; Gong, Y.; Limera, C.; Liu, L. Transcriptome profiling of radish (*Raphanus sativus* L.) root and identification of genes involved in response to lead (Pb) stress with next generation sequencing. *PLoS ONE* **2013**, *8*, e66539. [CrossRef] [PubMed]

34. Thakur, S.; Choudhary, S.; Bhardwaj, P. Comparative transcriptome profiling under cadmium stress reveals the uptake and tolerance mechanism in *Brassica juncea*. *J. Plant Growth Regul.* **2019**, *38*, 1141–1152. [CrossRef]

35. Di Toppi, L.S.; Gabbrielli, R. Response to cadmium in higher plants. *Environ. Exp. Bot.* **1999**, *41*, 105–130. [CrossRef]

36. Allan, D.L.; Jarrell, W.M. Proton and copper adsorption to maize and soybean root cell walls. *Plant Physiol.* **1989**, *89*, 823–832. [CrossRef]

37. Peng, H.-Y.; Yang, X.-E.; Tian, S.-K. Accumulation and ultrastructural distribution of copper in *Elsholtzia splendens*. *J. Zhejiang Univ. Sci. B* **2005**, *6*, 311. [CrossRef]

38. Wang, P.; Zhang, S.; Wang, C.; Lu, J. Effects of Pb on the oxidative stress and antioxidant response in a Pb bioaccumulator plant *Vallisneria natans*. *Ecotoxicol. Environ. Saf.* **2012**, *78*, 28–34. [CrossRef]

39. Satoh-Nagasawa, N.; Mori, M.; Nakazawa, N.; Kawamoto, T.; Nagato, Y.; Sakurai, K.; Takahashi, H.; Watanabe, A.; Akagi, H. Mutations in rice (*Oryza sativa*) heavy metal ATPase 2 (OsHMA2) restrict the translocation of zinc and cadmium. *Plant Cell Physiol.* **2012**, *53*, 213–224. [CrossRef]

40. Hasan, M.; Alabdallah, N.M.; Alharbi, B.M.; Waseem, M.; Yao, G.; Liu, X.-D.; El-Gawad, A.; Hany, G.; El-Yazied, A.A.; Ibrahim, M.F. GABA: A Key Player in Drought Stress Resistance in Plants. *Int. J. Mol. Sci.* **2021**, *22*, 10136. [CrossRef]

41. Daş, Z.A.; Dimlioğlu, G.; Bor, M.; Özdemir, F. Zinc induced activation of GABA-shunt in tobacco (*Nicotiana tabaccum* L.). *Environ. Exp. Bot.* **2016**, *122*, 78–84. [CrossRef]

42. Mahmud, J.A.; Hasanuzzaman, M.; Nahar, K.; Rahman, A.; Hossain, M.; Fujita, M. γ-aminobutyric acid (GABA) confers chromium stress tolerance in *Brassica juncea* L. by modulating the antioxidant defense and glyoxalase systems. *Ecotoxicology* **2017**, *26*, 675–690. [CrossRef] [PubMed]

43. Khoudi, H.; Maatar, Y.; Gouiaa, S.; Masmoudi, K. Transgenic tobacco plants expressing ectopically wheat H⁺-pyrophosphatase (H⁺-PPase) gene TaVP1 show enhanced accumulation and tolerance to cadmium. *J. Plant Physiol.* **2012**, *169*, 98–103. [CrossRef] [PubMed]

44. Monné, M.; Miniero, D.V.; Obata, T.; Daddabbo, L.; Palmieri, L.; Vozza, A.; Nicolardi, M.C.; Fernie, A.R.; Palmieri, F. Functional characterization and organ distribution of three mitochondrial ATP–Mg/Pi carriers in *Arabidopsis thaliana*. *Biochim. Biophys. Acta BBA Bioenerg.* **2015**, *1847*, 1220–1230. [CrossRef] [PubMed]

45. Block, M.D.; Verduyn, C.; Brouwer, D.D.; Cornelissen, M. Poly (ADP-ribose) polymerase in plants affects energy homeostasis, cell death and stress tolerance. *Plant J.* **2005**, *41*, 95–106. [CrossRef]
46. Singh, H.P.; Mahajan, P.; Kaur, S.; Batish, D.R.; Kohli, R.K. Chromium toxicity and tolerance in plants. *Environ. Chem. Lett.* **2013**, *11*, 229–254. [CrossRef]

47. Sharma, A.; Kumar, V.; Shahzad, B.; Ramakrishnan, M.; Sidhu, G.P.S.; Bali, A.S.; Handa, N.; Kapoor, D.; Yadav, P.; Khanna, K. Photosynthetic response of plants under different abiotic stresses: A review. *J. Plant Growth Regul.* **2020**, *39*, 509–531. [CrossRef]

48. Drummond, R.; Tutone, A.; Li, Y.-C.; Gardner, R. A putative magnesium transporter AtMRS2-11 is localized to the plant chloroplast envelope membrane system. *Plant Sci.* **2006**, *170*, 78–89. [CrossRef]

49. Jha, S.K.; Sharma, M.; Pandey, G.K. Role of cyclic nucleotide gated channels in stress management in plants. *Curr. Genom.* **2016**, *17*, 315–329. [CrossRef]

50. Moon, J.Y.; Belloeil, C.; Ianna, M.L.; Shin, R. *Arabidopsis* CNGC family members contribute to heavy metal ion uptake in plants. *Int. J. Mol. Sci.* **2019**, *20*, 413. [CrossRef]

51. Nieves-Cordones, M.; Alemán, F.; Martínez, V.; Rubio, F. K⁺ uptake in plant roots. The systems involved, their regulation and parallels in other organisms. *J. Plant Physiol.* **2014**, *171*, 688–695. [CrossRef]

52. Kohli, S.K.; Handa, N.; Gautam, V.; Bali, S.; Sharma, A.; Khanna, K.; Arora, S.; Thukral, A.K.; Ohri, P.; Karpets, Y.V. ROS Signaling in Plants Under Heavy Metal Stress. In *Reactive Oxygen Species and Antioxidant Systems in Plants: Role and Regulation under Abiotic Stress*; Khan, M., Khan, N., Eds.; Springer: Singapore, 2017. [CrossRef]

53. Hasanuzzaman, M.; Bhuyan, M.; Nahar, K.; Hossain, M.; Mahmud, J.A.; Hossen, M.; Masud, A.A.C.; Fujita, M. Potassium: A vital regulator of plant responses and tolerance to abiotic stresses. *Agronomy* **2018**, *8*, 31. [CrossRef]

54. Akram, W.; Khan, W.U.; Shah, A.A.; Yasin, N.A.; Li, G. Liquiritoside Alleviated Pb Induced Stress in *Brassica rapa* subsp. *Parachinensis*: Modulations in Glucosinolate Content and Some Physiochemical Attributes. *Front. Plant Sci.* **2021**, *12*, 1799. [CrossRef] [PubMed]

55. Chowdhury, P. Glucosinolates and its role in mitigating abiotic and biotic stress in Brassicaceae. In *Plant Stress Physiology—Perspectives in Agriculture*; IntechOpen: London, UK, 2022.

56. Wu, H.; Chen, C.; Du, J.; Liu, H.; Cui, Y.; Zhang, Y.; He, Y.; Wang, Y.; Chu, C.; Feng, Z. Co-overexpression FIT with AtbHLH38 or AtbHLH39 in *Arabidopsis*-enhanced cadmium tolerance via increased cadmium sequestration in roots and improved iron homeostasis of shoots. *Plant Physiol.* **2012**, *158*, 790–800. [CrossRef] [PubMed]

57. Manzoor, Z.; Hassan, Z.; Ul-Allah, S.; Khan, A.A.; Sattar, A.; Shahzad, U.; Amin, H.; Hussain, M. Transcription factors involved in plant responses to heavy metal stress adaptation. In *Plant Perspectives to Global Climate Changes*; Elsevier: Amsterdam, The Netherlands, 2022; pp. 221–231.

58. Dutta, S.; Mitra, M.; Agarwal, P.; Mahapatra, K.; De, S.; Sett, U.; Roy, S. Oxidative and genotoxic damages in plants in response to heavy metal stress and maintenance of genome stability. *Plant Signal. Behav.* **2018**, *13*, e1460048. [CrossRef] [PubMed]

59. Sharma, A.; Shahzad, B.; Rehman, A.; Bhardwaj, R.; Landi, M.; Zheng, B. Response of phenylpropanoid pathway and the role of polyphenols in plants under abiotic stress. *Molecules* **2019**, *24*, 2452. [CrossRef] [PubMed]

60. Van de Mortel, J.E.; Schat, H.; Moerland, P.D.; Van Themaat, E.V.L.; Van Der Ent, S.; Blankestijn, H.; Ghandilyan, A.; Tsiatsiani, S.; Aarts, M.G. Expression differences for genes involved in lignin, glutathione and sulphate metabolism in response to cadmium in *Arabidopsis thaliana* and the related Zn/Cd-hyperacumulator *Thlaspi caerulescens*. *Plant Cell Environ.* **2008**, *31*, 301–324. [CrossRef]

61. Skibbe, M.; Qu, N.; Galis, I.; Baldwin, I.T. Induced plant defenses in the natural environment: *Nicotiana attenuata* WRKY3 and WRKY6 coordinate responses to herbivory. *Plant Cell* **2008**, *20*, 1984–2000. [CrossRef]

62. Raza, A.; Charagh, S.; Zahid, Z.; Mubarik, M.S.; Javed, R.; Siddiqui, M.H.; Hasanuzzaman, M. Jasmonic acid: A key frontier in conferring abiotic stress tolerance in plants. *Plant Cell Rep.* **2021**, *40*, 1513–1541. [CrossRef]

63. Chen, X.; Jiang, W.; Tong, T.; Chen, G.; Zeng, F.; Jang, S.; Gao, W.; Li, Z.; Mak, M.; Deng, F. Molecular interaction and evolution of jasmonate signaling with transport and detoxification of heavy metals and metalloids in plants. *Front. Plant Sci.* **2021**, *12*, 625. [CrossRef]

64. Lei, G.J.; Sun, L.; Sun, Y.; Zhu, X.F.; Li, G.X.; Zheng, S.J. Jasmonic acid alleviates cadmium toxicity in *Arabidopsis* via suppression of cadmium uptake and translocation. *J. Integr. Plant Biol.* **2020**, *62*, 218–227. [CrossRef]

65. Liu, T.; Liu, S.; Guan, H.; Ma, L.; Chen, Z.; Gu, H.; Qu, L.-J. Transcriptional profiling of *Arabidopsis* seedlings in response to heavy metal lead (Pb). *Environ. Exp. Bot.* **2009**, *67*, 377–386. [CrossRef]

66. Ahmad, P.; Raja, V.; Ashraf, M.; Wijaya, L.; Bajguz, A.; Alyemeni, M.N. Jasmonic acid (JA) and gibberellic acid (GA3) mitigated Cd-toxicity in chickpea plants through restricted cd uptake and oxidative stress management. *Sci. Rep.* **2021**, *11*, 19768. [CrossRef] [PubMed]

67. Acosta, I.F.; Farmer, E.E. Jasmonates. In *The Arabidopsis Book*; The American Society of Plant Biologists: Rockville, MD, USA, 2010; pp. 8–129.

68. Wasternack, C.; Hause, B. Jasmonates: Biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany. *Ann. Bot.* **2013**, *111*, 1021–1058. [CrossRef] [PubMed]

69. Li, C.; Liu, G.; Xu, C.; Lee, G.I.; Bauer, P.; Ling, H.-Q.; Ganal, M.W.; Howe, G.A. The tomato suppressor of prosystemin-mediated responses 2 gene encodes a fatty acid desaturase required for the biosynthesis of jasmonic acid and the production of a systemic wound signal for defense gene expression. *Plant Cell* **2003**, *15*, 1646–1661. [CrossRef]

70. Li, Q.; Lei, S.; Du, K.; Li, L.; Pang, X.; Wang, Z.; Wei, M.; Fu, S.; Hu, L.; Xu, L. RNA-seq based transcriptomic analysis uncovers α-linolenic acid and jasmonic acid biosynthesis pathways respond to cold acclimation in *Camellia japonica*. *Sci. Rep.* **2016**, *6*, 36463. [CrossRef]
71. Upchurch, R.G. Fatty acid unsaturation, mobilization, and regulation in the response of plants to stress. *Biotechnol. Lett.* **2008**, *30*, 967–977. [CrossRef]

72. Beisson, F.; Li, Y.; Bonaventure, G.; Pollard, M.; Ohlrogge, J.B. The acyltransferase GPAT5 is required for the synthesis of suberin in seed coat and root of *Arabidopsis*. *Plant Cell* **2007**, *19*, 351–368. [CrossRef]

73. Kachroo, P.; Kachroo, A.; Lapchyk, L.; Hildebrand, D.; Klessig, D.F. Restoration of defective cross talk in *ssi2* mutants: Role of salicylic acid, jasmonic acid, and fatty acids in *SSI2*-mediated signaling. *Mol. Plant Microbe Interact.* **2003**, *16*, 1022–1029. [CrossRef]