Next Article in Journal
FcγRs and Their Relevance for the Activity of Anti-CD40 Antibodies
Next Article in Special Issue
The Citrus Laccase Gene CsLAC18 Contributes to Cold Tolerance
Previous Article in Journal
Effects of Sarcoptes scabiei Translationally Controlled Tumor Protein (TCTP) on Histamine Release and Degranulation of KU812 Cells
Previous Article in Special Issue
Overexpression of MxbHLH18 Increased Iron and High Salinity Stress Tolerance in Arabidopsis thaliana
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification and Variation Analysis of JAZ Family Reveals BnaJAZ8.C03 Involved in the Resistance to Plasmodiophora brassicae in Brassica napus

Key Laboratory of Biology and Genetic Improvement of Oil Crops, Ministry of Agriculture and Rural Affairs, Oil Crop Research Institute, Chinese Academy of Agricultural Sciences, Wuhan 430062, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(21), 12862; https://doi.org/10.3390/ijms232112862
Submission received: 22 September 2022 / Revised: 10 October 2022 / Accepted: 20 October 2022 / Published: 25 October 2022
(This article belongs to the Special Issue Function and Mechanism Analysis of Plant Stress Resistance Genes)

Abstract

:
Clubroot caused by Plasmodiophora brassicae led to a significant decrease in the yield and quality of Brassica napus, one of the most important oil crops in the world. JAZ proteins are an essential repressor of jasmonates (JAs) signaling cascades, which have been reported to regulate the resistance to P. brassicae in B. napus. In this study, we identified 51, 25 and 26 JAZ proteins in B. napus, B. rapa and B. oleracea, respectively. Phylogenetic analysis displayed that the notedJAZ proteins were divided into six groups. The JAZ proteins clustered in the same group shared a similar motif composition and distribution order. The 51 BnaJAZs were not evenly assigned on seventeen chromosomes in B. napus, except for A04 and C07. The BnaJAZs of the AtJAZ7/AtJAZ8 group presented themselves to be significantly up-regulated after inoculation by P. brassicae. Variation analysis in a population with a specific resistance performance in P. brassicae displayed a 64 bp translocation in BnaC03T0663300ZS (BnaJAZ8.C03, homologous to AtJAZ8) with an 8% reduction in the disease index on average. Through protein–protein interaction analysis, 65 genes were identified that might be involved in JAZ8 regulation of resistance to P. brassicae in B. napus, which provided new clues for understanding the resistance mechanism to P. brassicae.

1. Introduction

Jasmonates (JAs) are one of the important defense-signaling hormones [1], which are derived from α-linoleic acid released by the oxidation of membrane lipids [2]. The main bioactive form of JAs, JA-Ile, is perceived to ultimately trigger a set of transcriptional reprogramming and produce specialized metabolites or proteins in response to abiotic/biotic stresses [3,4]. The core JAs signaling pathway constitutes several functional modules, including a repertoire of COI1–JAZ (Coronatine Insensitive1–Jasmonate-Zim Domain) co-receptors [5]. COI1 is identified to encode an F-box protein [6], which provides substrate specificity for the Skp–Cullin–F-box-type E3 ubiquitin ligase complex (SCFCOI1) [7]. JAZ proteins target SCFCOI1 [8], which plays a central role within the signaling cascades triggered by JAs. In plant cells with a low concentration of JAs, the expression of JA early response genes is inhibited by JAZ proteins binding to transcription factors (TFs) or an adaptor protein NINJA [9]. When external stimuli stimulate plants, the content of JAs in plants increases, accelerating the binding of SCFCOI1 with JAZ protein and JAZ ubiquitination. As a result, the JAZ proteins are degraded by 26S proteasome, releasing TFs inhibited by JAZ proteins and finally activating the expression of JA early response genes [10].
JAZ proteins are the plant-specific zinc finger proteins belonging to the superfamily of TIFY [11]. The JAZ proteins contain a highly conserved TIF[F/Y]XG motif within the ZIM domain (also known as TIFY domain) near its N-terminal, which mediates the interaction of proteins [12,13,14]. JAZ proteins also contain a conserved Jas domain (also known as CCT_2 domain) composed of 12–29 amino acids (SLX2FX2KRX2RX5PY) near its C-terminal [15], a variant of CCT motif. The Jas domain plays a dual role in JAs signaling pathway, mediates the interaction between JAZ and its downstream TFs, and promotes the combination with COI1 [16]. The JAZ proteins lacking a Jas sequence showed JA-dominant insensitivity [17]. The JAZ family has large numbers of genes. Previously, the 15, 14 and 16 JAZ genes were found in rice [18], wheat [19] and maize [20], respectively. In Arabidopsis thaliana, 13 JAZ genes have been reported. Among them, AtJAZ1-12 encodes typical TIFY JAZ proteins, which conserve the ZIM and TIFY domains. AtJAZ13 encodes a non-TIFY JAZ protein that contains divergent ZIM, Jas and EAR motifs. Compared with other JAZ proteins, AtJAZ13 has a ser-rich C-terminal tail for phosphorylation [21]. The identified JAZ genes have been proven to participate in plant development and defense reactions. For example, AtJAZ2 prevents pathogen penetration by triggering stomata closure [22]. The overexpression JAZ2 from Glycine soja could enhance tolerance to salt and alkali stresses [23]. In addition, the overexpression JAZ4 from grapes could improve the resistance to powdery mildew [24]. The expression of JAZ1 and JAZ5 could be affected by WRKY57, which compromised the resistance to Botrytis cinerea in Arabidopsis [25]. In addition, OsJAZ9 has been reported to regulate salt stress or potassium deficiency tolerance in rice [26]. AtJAZ10 was induced by wounding and negatively regulated the disease symptom development during Pseudomonas syringae pathogenesis [27].
Brassica napus (AACC) is one of the most important oil crops in China and the world, whose safety production is very important for supplying edible vegetable oil. However, clubroot, caused by Plasmodiophora brassicae, has rapidly spread worldwide, resulting in significant decreases in yield and quality in B. napus that result in significant economic losses yearly [28].
Breeding varieties with durable and stable resistance is the most economical and effective method to prevent clubroot. Moreover, identifying genes and regulatory pathways is the key to the genetic improvement of resistance to P. brassicae. In a previous study, we found that JAZ proteins, which are the repressors of JAs signaling pathway, are involved in the regulation of resistance to P. brassicae through compared transcriptome analysis in B. napus [29]. Even so, there are many members of the JAZ family in B. napus. It is unknown which JAZ proteins are involved in the regulation of resistance to P. brassicae in B. napus exactly. We focused on the JAZ genes in this study and performed genome-wide identification in B. napus using bioinformatics analysis. Furthermore, we analyzed the chromosome location and structure characteristics. Then, phylogenetic analysis of JAZs among B. rapa, B. oleracea and B. napus was conducted. In addition, expression levels of BnaJAZs in response to P. brassicae and variations of differential expressional BnaJAZs in the population were analyzed. The results of this study could provide critical clues for us to understand the mechanism of resistance to P. brassicae in B. napus.

2. Results

2.1. Genome-Wide Identification of JAZ in B. napus, B. rapa and B. oleracea

To identify all the JAZ proteins in B. napus, B. rapa and B. oleracea, we performed a Blastp and HMMER search simultaneously against the protein sequences of B. napus (ZS11), B. rapa (Chiifu V3.0) and B. oleracea (JZS V2.0). The merged candidate proteins were further analyzed using NCBI-CDD and SMART to ascertain the integrity of the conserved domains. Consequently, we identified 51, 25 and 26 JAZs in B. napus, B. rapa and B. oleracea, respectively. The amino acid (aa) residues of all the 102 JAZs ranged from 116 to 573, with an average length of 258 aa. The molecular weight of these JAZs varied from 13.1 to 62.8 kDa. In addition, the isoelectric point values were distributed from 4.96 to 10.07. The negative GRAVY scores indicated that all the JAZs were hydrophilic proteins. Subcellular localization showed that all the JAZs were predicted in the nuclei. It is worth noting that besides the nucleus, BnaC03T0060300ZS and BnaA03T0052500ZS were also found in the cell membrane or chloroplast (Table S2), which may affect multi-located proteins.

2.2. Phylogenetic Analysis of JAZ Family among Arabidopsis, B. napus, B. rapa, and B. oleracea

Based on the JAZ protein sequences from Arabidopsis, B. napus, B. rapa and B. oleracea, a neighbor-joining (NJ) phylogenetic tree was constructed using MEGA11 software (Figure 1). The phylogenetic analysis revealed that all the JAZs could be divided into six groups: group I (purple, homologous to AtJAZ7/AtJAZ8), group II (yellow, homologous to AtJAZ10), group III (red, homologous to AtJAZ3/AtJAZ4/AtJAZ9), group IV (blue, homologous to AtJAZ11/AtJAZ12), group V (green, homologous to AtJAZ5/AtJAZ6), and group VI (orange, homologous to AtJAZ1/AtJAZ2). The tree showed that group VI contained the maximal JAZ members (24 proteins), while group IV had the least members (seven proteins). We found that the number of JAZ in B. napus was equal to the total of B. oleracea and B. rapa in most groups. Nevertheless, in group V, the total number of JAZs in B. oleracea and B. rapa was more than that in B. napus. Furthermore, we found the numbers of JAZs on the A genome from B. rapa, and C genome from B. oleracea were equal to that on A-subgenome and C-subgenome from B. napus, respectively, in all the groups except for the group I, IV, and V. Furthermore, AtJAZ1, AtJAZ2, AtJAZ5, and AtJAZ9 had six copies in B. napus, and also had three copies in both B. rapa and B. oleracea. Similarly, AtJAZ3 had four copies in B. napus and two copies in B. rapa and B. oleracea.
These results further confirmed the genome duplication events in B. napus speciation. It was worth noting that AtJAZ10 had five copies in both B. rapa and B. oleracea, and 10 copies in rapeseed. However, AtJAZ4 did not have a homologue in B. rapa or B. oleracea, which meant that AtJAZ4 was the unique gene of Arabidopsis. Moreover, AtJAZ7 only had one homologue in B. oleracea and no homologue in B. rapa. Correspondingly, AtJAZ11 only had one homologue in B. rapa and no homologue in B. oleracea. It was speculated that these two genes might be lost when B. rapa or B. oleracea speciation. In addition, AtJAZ12 only had one copy in both B. rapa and B. oleracea, which implied that AtJAZ12 was not expanded. The above results implied that the different JAZ genes/clades might be divisional during the formation and evolution of Brassica species.

2.3. Conserved Motifs and Gene Structure Analysis Displayed That the JAZ Proteins Clustered in the Same Group Shared Similar Motif Composition and Distribution Order

The 102 JAZ proteins were typical and contained both the Jas domain and ZIM domain. Based on that, BnaA06T0133300ZS and BnaA09T0610300ZS had Abhydrolase_3 domains at their C-terminal as well. Conserved motifs analysis showed that 10 conserved motifs were detected, and motifs 2 and 4 existed in all the JAZs. Therefore, it was speculated that sequences of motifs 2 and 4 encoded the conserved functional domain of JAZ protein. However, the distribution of motifs 2 and 4 was variable in different JAZs (Figure 2A). In addition, motif 6 was found only in a branch (homologous to AtJAZ9) of group III. Meanwhile, motif 9 was found only in a branch of group II. Notably, besides all the members of group VI, motif 7 was also presented in BolC08g014090.2J and BnaC08T0129100ZS belonging to group II. Similarly, besides all the members of group III, motif 10 was also found in a branch (homologous to AtJAZ6) of group V. On the whole, the JAZ proteins, clustered in the same group, shared similar motif composition and distribution order. For example, group VI, including the maximal JAZ members, had a conserved “motif 7-2-4-1” and a variable motif 3 or 8. The results showed that the ZIM domain is mostly located in motif 2, and the Jas domain is mostly located in motif 1 or 5. In more detail, the Jas domains of cluster I (group III, IV, V, and VI) were located in motif 1. Instead, the Jas domains of cluster II (group I and II) were located in motif 5. Based on that, it was speculated that the functional difference between clusters I and II might be caused by the difference between motifs 1 and 5. Gene structure displayed that the number of exons and introns of all the JAZ genes was varied from 2 to 9 and 1 to 8, respectively (Figure 2B). It could be found that group I had the least exons (2–3), and the exon numbers of group II had the widest range of variation (3–9).

2.4. Chromosomal Location and Expression Pattern of JAZ Genes in B. napus under Multiple Treatments

All the 51 BnaJAZs were assigned to seventeen chromosomes in B. napus except for A04 and C07, 24 of which were located in the A sub-genome and 27 in the C sub-genome (Figure 3). The results displayed that the distribution of BnaJAZs on each chromosome was uneven. The most BnaJAZs (five genes) were present on the chromosomes A02, A07, C02, C06 and C08. Based on the gene expression data from BnIR (http://yanglab.hzau.edu.cn/, accessed on 18 August 2022), it was displayed the expression change of JAZ genes in ZS11 under hormone treatment, including IAA, GA, ABA, JA and adversity treatments, including salt, drought and freezing. After IAA, GA, and ABA treatment, the expression level of all JAZs did not present significant changes in the roots and leaves compared with the control. After JA treatment, the homologous genes of AtJAZ1/2/6/9 in B. napus showed significantly up-regulated expression in both roots and leaves, while homologous genes of other AtJAZs showed certain differentiation on expression level. The homologous genes of AtJAZ7/8/10 were unchanged in leaves, but some were up-regulated in the roots. About seven BnaJAZs were unchanged in both the roots and leaves. Furthermore, we analyzed the expression level of all 51 BnaJAZs under abiotic stress. The results showed that the homologue of AtJAZ1 in B. napus showed a response to salt, drought and freeze damage stress, especially in roots. However, the homologue of AtJAZ9 was strongly up-regulated only in leaves after 3 h under drought stress (Figure 4).

2.5. Haplotype Analysis of SNPs and SVs Reveals BnaJAZ8.C03 Associated with the Resistance to P. brassicae

To explore the response of the JAZ family to P. brassicae, we conducted qRT-PCR to analyze the difference in JAZ gene expression between inoculated and un-inoculated by P. brassicae in resistant (“28669”) and susceptible (YJ8) line. The results showed that there were no patterns in the expression of BnaJAZs in groups II, III, V, and VI. Meanwhile, the expression differential fold of these genes between inoculation and un-inoculation was small, except for an obvious differential expressed gene (DEG), BnaA02T0200100ZS, homologous to AtJAZ6 in group V (Figure S1). It was up-regulated fivefold in “28669” compared with the control, while the change was less than 1.5-fold in YJ8. In addition, five BnaJAZs belonging to group I (homologous to AtJAZ7/8) showed obvious up-regulation in both resistant and susceptible lines. The degree of up-regulation in resistant was higher than in susceptible (Figure 5A). Likewise, the expression of BnaJAZs belonging to group IV (homologous to AtJAZ11/12) was up-regulated, but the degree of which in susceptible was higher than that in resistant (Figure 5B). Because the expression differential fold of BnaJAZs belonging to group IV between inoculation and un-inoculation was very small, we further performed the SNPs and SVs detection only in BnaJAZs of group I and BnaA02T0200100ZS using the whole genome re-sequencing data including 418 rapeseed accessions [30]. A total of 17 SNPs were detected in four BnaJAZs, and three SVs were detected in two BnaJAZs. Among them, the BnaA02T0200100ZS possessed the maximal SNPs (10), and the BnaC03T0663300ZS had the maximal SVs (2). In conjunction with this population’s resistance to P. brassicae, we further analyzed whether the SNPs or SVs of DEGs affected the resistance to P. brassicae. The results showed no significant differences in resistance to P. brassicae among different haplotypes divided by SNPs in all the above BnaJAZs. After mapping the re-sequence data by pair-end sequencing through an Illumina HiSeq 4000 to the reference genome ZS11, the gray bars represent that normal reads which paired mapped to the BnaC03T0663300ZS (BnaJAZ8.C03, homologous to AtJAZ8) region (Figure S2). While the peak green bars represent the abnormal reads, one end was not mapped to the BnaC03T0663300ZS region of chrC03, but to the chrA08 region (Figure 6A). Based on the pair-end reads mapping position of the BnaC03T0663300ZS region in the natural population containing 418 accessions, the accessions can be divided into two groups, one containing this translocation and the other not containing the translocation at the BnaC03T0663300ZS region (Figure 6A and Figure S2).
The phenotypes of resistant and susceptible lines are shown in Figure 6B. There were 368 accessions that did not contain this translocation, and 37 accessions contained this translocation in the population. Thirteen accessions have not acquired the phenotype. The average disease indexes of accessions with and without this translocation were 74.4%, and 68.5%, respectively. Compared with the accessions without this translocation, the average disease indexes of accessions with this translocation had an 8% reduction. The t-test showed that the difference reached a significant level (Figure 6C). The above results implied that the 64 bp translocation in BnaJAZ8.C03 could significantly change resistance, which may involve the resistance to P. brassicae in B. napus.

2.6. Interaction Network of Key JAZs Identified Their Functional Partner in B. napus

Based on the related proteins in Arabidopsis, we predicted the potential proteins interacting with JAZ8 using the STRING database. A total of 21 proteins were detected in the PPI network, including four JAZs, which proved that JAZ8 might interact with some JAZ proteins, such as JAZ1/5/10. In addition to some well-known important components of the JAs signaling pathway, such as COI, NINJA, MYC, bHLH and TPL, other proteins have also been identified, such as IAAs, WRKY57, TCP15, SKP1, TRO and so on (Figure 7). It has been reported that WRKY57 affected the expression of JA ZIM-Domain genes transcriptionally to compromise Botrytis cinerea resistance. The overexpression of AtWRKY57 could enhance the drought tolerance in Arabidopsis and rice [25,31]. WRKY57 have also been reported to participate in the suppression of abscisic acid at the early infection stage of Verticillium longisporum in B. napus [32].
Furthermore, SKP1, a core component of the SCF family of E3 ubiquitin ligases, was also detected, which provides specificity in binding to ubiquitin ligase substrate proteins. Studies showed that the SKP1-like protein could influence the resistance to powdery mildew fungus by controlling the abundance of the susceptibility factor RACB in barley [33]. Notably, four IAAs (IAA8/10/16/19) were predicted in this interaction network, suggesting that besides the JA signal transduction pathway, JAZ is also involved in IAA biosynthesis or signal transduction. Based on the homologous collinearity between B. napus and Arabidopsis, 65 genes were identified that likely took part in the JAZ8 regulation of clubroot resistance in B. napus in addition to BnaJAZs.

3. Discussion

JAZ proteins have been reported to participate in the regulation of stress responses in plants. Comprehensive studies of JAZ gene family have been carried out in maize [20], tea plants [34], cotton [35] and other plant species. Most of the time, JAZ, as a subfamily, was mentioned in TIFY gene family research on rice [18], wheat [36], grape [37], soybean [38], cucumber [39] and so on. In this study, we only focused on the JAZ family. There were 51, 25 and 26 JAZ genes identified in B. napus, B. rapa and B. oleracea, respectively. This result is inconsistent with previous reports, in which 52, 26, and 27 JAZ genes have been reported in B. napus, B. rapa and B. oleracea [40]. The 22 and 21 JAZ genes have also been reported in B. oleracea [41] and B. rapa [42]. We considered that this inconsistency might be due to the difference in the reference genome. On the other hand, the divergent conserved domain formed the non-typical TIFY-JAZ, which also could result in this difference.
In general, the result of the phylogenetic tree was similar to the previous study based on the AtJAZs [43]. Because B. napus was generated 5000 to 10,000 years ago by the natural hybridization of its two progenitor diploids, B. rapa and B. oleracea, B. napus was an ideal model species to study the evolution of allopolyploid. Moreover, B. rapa and B. oleracea were produced by extensive triplication of their ancestral species at the genomic level. The three noted species are believed to share a common ancestor with Arabidopsis. Therefore, one gene in Arabidopsis had an average of 4–6 copies in rapeseed. It could be that most orthologous gene pairs in B. rapa and B. oleracea remained as homologous pairs in B. napus from the tree [44]. In addition, the phylogenetic tree showed the triplication events of the Brassica species in the evolution of JAZ gene family to a certain extent. Some of the AtJAZs had three homologues in the Ar or Co genomes, and some AtJAZs experienced expansion or loss after triplication.
Genomic structural variations (SVs) refer to DNA genetic variation over 50 bp in length and are comprised of insertion/deletion (InDel), duplication, inversion and translocation. Numerous studies to identify genes associated with important traits were based on single nucleotide polymorphism (SNP) data, while the SVs were ignored. The accuracy of SVs could be effectively improved by integrating the results of four software SV detection. It laid a foundation for subsequent research. However, the accuracy of SVs needs a further increase compared with the SNPs by improving the sequencing technology and bioinformatics computing methods. SVs not only lead to the fusion of different genes and generate new genes but also play a decisive role in the genetic domestication of crops and the regulation of agronomic trait genes in maize [45], soybeans [46], tomatoes [47], peaches [48,49,50] and other crops.
In this study, a translocation of a key JAZ, BnaC03T0663300ZS, was detected to be associated with the resistance to P. brassicae in B. napus. This gene was homologous to the AtJAZ8, which has been reported to participate in the signaling response during Agrobacterium infection to Arabidopsis. The overexpression of AtJAZ8 in Arabidopsis could inhibit the activity of VirE3 as a plant transcriptional regulator, activating AtPR1 gene expression and repressing the expression of AtPDF1.2, which significantly decreased the numbers of tumors formed [51]. In addition, BnaA02T0200100ZS was screened through the differential expression analysis, although no SV was detected. BnaA02T0200100ZS was homologous to AtJAZ6, reported previously, and played a crucial role in the rhythmic susceptibility of Arabidopsis to B. cinerea [52].
Based on the above research, we identified 65 proteins that interact with JAZ8 directly or indirectly in B. napus through protein–protein interaction analysis and homologous colinearity. Some genes have been reported to participate in Jas’ signaling pathways; others have also been reported to be involved in regulating host resistance to other diseases. The above information could provide some clues for understanding the resistance mechanism of P. brassicae in B. napus.

4. Materials and Methods

4.1. Plant Materials, Sequencing and Phenotype Data

This study obtained the seeds of 418 rapeseed accessions from the National Mid-term Gene Bank of Oil Crops Research Institute, Chinese Academy of Agriculture Sciences (OCRI, CAAS). The re-sequencing data of 418 accessions were downloaded from our previous study [30]. The resistance to P. brassicae data of 418 accessions was collected by artificial inoculation in 2015 (China, Shenyang) [53]. Among the accessions, we identified a resistant accession, “28669”, and a susceptible accession “, YJ-8”, which came from the National Mid-term Gene Bank of OCRI, CAAS. The total RNA of the root tissue 60 h post-inoculation by P. brassicae of “28669” and “YJ-8” were extracted to conduct the qRT-PCR [29]. The B. napus, B. oleracea, and B. rapa materials for sequence analyses were cv. ZS11 (V1.0) [54], JZS (V2.0) [55], and Chiifu (V3.0) [56], respectively.

4.2. Identification of the JAZs in B. napus, B. rapa and B. oleracea

To identify the JAZ family members, the genomic sequence of the B. napus cv. ZS11 HZAU V1.0 [54], B. oleracea cv. JZS V2.0 [55] and B. rapa cv. Chiifu V3.0 [56] were downloaded from the website http://www.brassicadb.cn/#/Download/, accessed on 28 April 2022. Twelve Arabidopsis typical TIFY-JAZ protein sequences were obtained from TAIR (https://www.arabidopsis.org/, accessed on 28 April 2022), which were used as the queries for searching the JAZ proteins in B. napus, B. rapa and B. oleracea by BLASTP (E-value < 1 × 10−6) using TBtools [57]. Meanwhile, the Hidden Markov Model (HMM) profiles of the Jas domain (PF09425) and ZIM domain (PF06200) were downloaded from Pfam (http://pfam.xfam.org/, accessed on 23 May 2022), which were used for protein screening. The conserved domains of candidate JAZs were further confirmed using NCBI-CDD (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi, accessed on 24 May 2022) and SMART (http://smart.embl-heidelberg.de/, accessed on 24 May 2022). Furthermore, the molecular weight (Mw) and theoretical isoelectric point (pI) of certain JAZ proteins were analyzed in ExPASy (https://web.expasy.org/compute_pi/, accessed on 11 June 2022). Moreover, the subcellular locations of all the JAZ proteins were predicted using ProtComp9.0 from Softberry (http://www.softberry.com/, accessed on 11 June 2022).

4.3. Phylogenetic, Gene Structure and Conserved Motif Analyses

To explore the evolutionary characteristics of JAZs from Arabidopsis, B. napus, B. rapa and B. oleracea, multiple JAZ protein sequences were aligned by MUSLE [58]. Then, the alignment results were used for constructing an un-rooted neighbor-joining (NJ) phylogenetic tree using MEGA 11 software [59] with 1000 bootstrap replicates. Finally, we beautified the phylogenetic tree in iTOL (https://itol.embl.de/, accessed on 26 May 2022). The conserved motifs were searched in MEME (https://meme-suite.org/meme/tools/meme, accessed on 27 May 2022), and the detection parameter of conserved motifs was set to 10 for fear of missing some rare motifs that may have specific functions. TBtools software was used for drawing the exons, introns, conserved domain coverage and motif distribution.

4.4. Chromosomal Location and Expression Pattern Analysis of JAZ Genes in B. napus under Multiple Treatment

Based on the genomic sequence annotation file of B. napus cv. ZS11 obtained from the BnPIR database (http://cbi.hzau.edu.cn/bnapus/, accessed on 28 Aprial 2022), the chromosomal positions of all the BnaJAZs were drawn using TBtools. The expression data of BnaJAZs in leaf and root under hormone and adversity treatment were acquired from BnIR (http://yanglab.hzau.edu.cn/BnIR, accessed on 18 August 2022). The qRT-PCR results analysis was performed on LightCycler 480 SYBR Green I Mastermix, and a LightCycler 480II real-time PCR system (Roche, Switzerland). The transcript abundance was calculated from three biological and three technical replicates with Bna.Actin7 as an internal control. The fold-change was estimated using the 2−ΔΔCt [60]. The gene-specific primer sequences (Table S1) and heat maps were also designed and drawn with TBtools.

4.5. Structural Variation Detection, Resistance Difference Analysis and Interaction Network Construction of Target Genes

Based on the re-sequencing data [30], we carried out the population’s SNPs and structural variation (SVs) detection. Therein, the SVs were detected by using three tools, LUMPY 0.2.13 [61], Manta 1.6.0 [62], and smoove 0.2.8 (https://github.com/brentp/smoove, accessed on 20 August 2022). To improve the accuracy, the variants detected by at least two tools were kept using SURVIVOR 1.0.6 [63]. SPASS software detected the significance of the resistance differences among different haplotypes divided by SNPs or SVs in each target gene. The interaction networks of the target JAZs were generated using the STRING database (https://cn.string-db.org/, accessed on 24 August 2022).

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232112862/s1. Figure S1: Expression profiles of BnaJAZs belonging to orange, green, red and yellow groups divided by phylogenetic analysis in the roots of resistant and susceptible accessions after 60 h infected by P. brassicae; Figure S2: The differential performance of reads in BnaC03T0663300ZS (BnaJAZ8.C03) between resistant and susceptible accessions; Table S1: Gene-specific primers sequence of qRT-PCR; Table S2: Features of 102 JAZ proteins in B. rapa, B. oleracea and B. napus.

Author Contributions

L.L. and G.C. conceived the overall study. L.L. performed and analyzed the experiments and wrote the manuscript. G.J., W.G. and F.Q. assisted with the bioinformatics analyses and graph diagrams. H.L. assisted with conducting the experiments. G.J., G.C. and X.W. revised the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by grants from the National Natural Science Foundation of China (grant no. 32201767, 32001892, and 32122063), the China Agriculture Research System (CARS-12), and the Agricultural Science and the Technology Innovation Project of the Chinese Academy of Agricultural Sciences (CAAS-ASTIP, CAAS-ZDRW202105).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nguyen, T.H.; Goossens, A.; Lacchini, E. Jasmonate: A hormone of primary importance for plant metabolism. Curr. Opin. Plant Biol. 2022, 67, 102197. [Google Scholar] [CrossRef] [PubMed]
  2. Wasternack, C.; Feussner, I. The Oxylipin Pathways: Biochemistry and Function. Annu. Rev. Plant Biol. 2018, 69, 363–386. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, X.; Wang, D.D.; Fang, X.; Chen, X.Y.; Mao, Y.B. Plant Specialized Metabolism Regulated by Jasmonate Signaling. Plant Cell Physiol. 2019, 60, 2638–2647. [Google Scholar] [CrossRef] [PubMed]
  4. Wasternack, C.; Strnad, M. Jasmonates are signals in the biosynthesis of secondary metabolites—Pathways, transcription factors and applied aspects—A brief review. New Biotechnol. 2019, 48, 1–11. [Google Scholar] [CrossRef]
  5. Zhang, F.; Yao, J.; Ke, J.; Zhang, L.; Lam, V.Q.; Xin, X.F.; Zhou, X.E.; Chen, J.; Brunzelle, J.; Griffin, P.R. Structural basis of JAZ repression of MYC transcription factors in jasmonate signalling. Nature 2015, 525, 269–273. [Google Scholar] [CrossRef] [Green Version]
  6. Xie, D. COI1: An Arabidopsis Gene Required for Jasmonate-Regulated Defense and Fertility. Science 1998, 280, 1091–1094. [Google Scholar] [CrossRef]
  7. Pauwels, L.; Goossens, A. The JAZ proteins: A crucial interface in the jasmonate signaling cascade. Plant Cell 2011, 23, 3089–3100. [Google Scholar] [CrossRef] [Green Version]
  8. Thines, B.; Katsir, L.; Melotto, M.; Niu, Y.; Mandaokar, A.; Liu, G.; Nomura, K.; He, S.Y.; Howe, G.A.; Browse, J. JAZ repressor proteins are targets of the SCF(COI1) complex during jasmonate signalling. Nature 2007, 448, 661–665. [Google Scholar] [CrossRef]
  9. Chico, J.M.; Chini, A.; Fonseca, S.; Solano, R. JAZ repressors set the rhythm in jasmonate signaling. Curr. Opin. Plant Biol. 2008, 11, 486–494. [Google Scholar] [CrossRef]
  10. Yan, J.; Li, H.; Li, S.; Yao, R.; Deng, H.; Xie, Q.; Xie, D. The Arabidopsis F-box protein CORONATINE INSENSITIVE1 is stabilized by SCFCOI1 and degraded via the 26S proteasome pathway. Plant Cell 2013, 25, 486–498. [Google Scholar] [CrossRef]
  11. Derelle, E.; Ferraz, C.; Rombauts, S.; Rouze, P.; Worden, A.Z.; Robbens, S.; Partensky, F.; Degroeve, S.; Echeynie, S.; Cooke, R.; et al. Genome analysis of the smallest free-living eukaryote Ostreococcus tauri unveils many unique features. Proc. Natl. Acad. Sci. USA 2006, 103, 11647–11652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Chung, H.S.; Howe, G.A. A critical role for the TIFY motif in repression of jasmonate signaling by a stabilized splice variant of the JASMONATE ZIM-domain protein JAZ10 in Arabidopsis. Plant Cell 2009, 21, 131–145. [Google Scholar] [CrossRef] [Green Version]
  13. Pauwels, L.; Barbero, G.F.; Geerinck, J.; Tilleman, S.; Grunewald, W.; Perez, A.C.; Chico, J.M.; Bossche, R.V.; Sewell, J.; Gil, E.; et al. NINJA connects the co-repressor TOPLESS to jasmonate signalling. Nature 2010, 464, 788–791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Chini, A.; Fonseca, S.; Chico, J.M.; Fernandez-Calvo, P.; Solano, R. The ZIM domain mediates homo- and heteromeric interactions between Arabidopsis JAZ proteins. Plant J. 2009, 59, 77–87. [Google Scholar] [CrossRef] [PubMed]
  15. Yan, Y.; Stolz, S.; Chetelat, A.; Reymond, P.; Pagni, M.; Dubugnon, L.; Farmer, E.E. A downstream mediator in the growth repression limb of the jasmonate pathway. Plant Cell 2007, 19, 2470–2483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Melotto, M.; Mecey, C.; Niu, Y.; Chung, H.S.; Katsir, L.; Yao, J.; Zeng, W.; Thines, B.; Staswick, P.; Browse, J.; et al. A critical role of two positively charged amino acids in the Jas motif of Arabidopsis JAZ proteins in mediating coronatine- and jasmonoyl isoleucine-dependent interactions with the COI1 F-box protein. Plant J. 2008, 55, 979–988. [Google Scholar] [CrossRef] [Green Version]
  17. Yamada, S.; Kano, A.; Tamaoki, D.; Miyamoto, A.; Shishido, H.; Miyoshi, S.; Taniguchi, S.; Akimitsu, K.; Gomi, K. Involvement of OsJAZ8 in jasmonate-induced resistance to bacterial blight in rice. Plant Cell Physiol. 2012, 53, 2060–2072. [Google Scholar] [CrossRef] [Green Version]
  18. Ye, H.; Du, H.; Tang, N.; Li, X.; Xiong, L. Identification and expression profiling analysis of TIFY family genes involved in stress and phytohormone responses in rice. Plant Mol. Biol. 2009, 71, 291–305. [Google Scholar] [CrossRef]
  19. Wang, Y.K.; Qiao, L.Y.; Bai, J.F.; Duan, W.J.; Yuan, S.H. Genome-wide characterization of JASMONATE-ZIM DOMAIN transcription repressors in wheat (Triticum aestivum L.). BMC Genom. 2017, 18, 152. [Google Scholar] [CrossRef] [Green Version]
  20. Han, Y.; Luthe, D. Identification and evolution analysis of the JAZ gene family in maize. BMC Genom. 2021, 22, 256. [Google Scholar] [CrossRef]
  21. Thireault, C.; Shyu, C.; Yoshida, Y.; St Aubin, B.; Campos, M.L.; Howe, G.A. Repression of jasmonate signaling by a non-TIFY JAZ protein in Arabidopsis. Plant J. 2015, 82, 669–679. [Google Scholar] [CrossRef]
  22. Gimenez-Ibanez, S.; Boter, M.; Ortigosa, A.; Garcia-Casado, G.; Chini, A.; Lewsey, M.G.; Ecker, J.R.; Ntoukakis, V.; Solano, R. JAZ2 controls stomata dynamics during bacterial invasion. New Phytol. 2017, 213, 1378–1392. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Zhu, D.; Cai, H.; Luo, X.; Bai, X.; Deyholos, M.K.; Chen, Q.; Chen, C.; Ji, W.; Zhu, Y. Over-expression of a novel JAZ family gene from Glycine soja, increases salt and alkali stress tolerance. Biochem. Biophys. Res. Commun. 2012, 426, 273–279. [Google Scholar] [CrossRef] [PubMed]
  24. Zhang, G.; Yan, X.; Zhang, S.; Zhu, Y.; Zhang, X.; Qiao, H.; van Nocker, S.; Li, Z.; Wang, X. The jasmonate-ZIM domain gene VqJAZ4 from the Chinese wild grape Vitis quinquangularis improves resistance to powdery mildew in Arabidopsis thaliana. Plant Physiol. Biochem. 2019, 143, 329–339. [Google Scholar] [CrossRef]
  25. Jiang, Y.; Yu, D. The WRKY57 transcription factor affects the expression of Jasmonate ZIM-Domain genes transcriptionally to compromise Botrytis cinerea resistance. Plant Physiol. 2016, 171, 2771–2782. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Xiong, L.Z.; Zhang, T.; Wu, H.; Ye, H.Y.; Yao, R.F. OsJAZ9 acts as a transcriptional regulator in jasmonate signaling and modulates salt stress tolerance in rice. Plant Sci. 2015, 232, 1–12. [Google Scholar]
  27. Demianski, A.J.; Chung, K.M.; Kunkel, B.N. Analysis of Arabidopsis JAZ gene expression during Pseudomonas syringae pathogenesis. Mol. Plant Pathol. 2012, 13, 46–57. [Google Scholar] [CrossRef]
  28. Chai, A.L.; Xie, X.W.; Shi, Y.X.; Li, B.J. Research status of clubroot (Plasmodiophora brassicae) on cruciferous crops in China. Can. J. Plant Pathol. 2014, 36 (Suppl. S1), 142–153. [Google Scholar] [CrossRef]
  29. Li, L.; Long, Y.; Li, H.; Wu, X. Comparative transcriptome analysis reveals key pathways and hub genes in rapeseed during the early stage of Plasmodiophora brassicae infection. Front. Genet. 2019, 10, 1275. [Google Scholar] [CrossRef] [Green Version]
  30. Hu, J.; Chen, B.; Zhao, J.; Zhang, F.; Xie, T.; Xu, K.; Gao, G.; Yan, G.; Li, H.; Li, L.; et al. Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding. Nat. Genet. 2022, 54, 694–704. [Google Scholar] [CrossRef]
  31. Jiang, Y.; Qiu, Y.; Hu, Y.; Yu, D. Heterologous expression of AtWRKY57 confers drought tolerance in Oryza sativa. Front. Plant Sci. 2016, 7, 145. [Google Scholar] [CrossRef] [PubMed]
  32. Behrens, F.H.; Schenke, D.; Hossain, R.; Ye, W.; Schemmel, M.; Bergmann, T.; Hader, C.; Zhao, Y.; Ladewig, L.; Zhu, W.; et al. Suppression of abscisic acid biosynthesis at the early infection stage of Verticillium longisporum in oilseed rape (Brassica napus). Mol. Plant Pathol. 2019, 20, 1645–1661. [Google Scholar] [CrossRef] [PubMed]
  33. Reiner, T.; Hoefle, C.; Huckelhoven, R. A barley SKP1-like protein controls abundance of the susceptibility factor RACB and influences the interaction of barley with the barley powdery mildew fungus. Mol. Plant Pathol. 2016, 17, 184–195. [Google Scholar] [CrossRef] [PubMed]
  34. Shen, J.; Zou, Z.; Xing, H.; Duan, Y.; Zhu, X.; Ma, Y.; Wang, Y.; Fang, W. Genome-wide analysis reveals stress and hormone responsive patterns of JAZ family genes in Camellia Sinensis. Int. J. Mol. Sci. 2020, 21, 2433. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Li, W.; Xia, X.C.; Han, L.H.; Ni, P.; Yan, J.Q.; Tao, M.; Huang, G.Q.; Li, X.B. Genome-wide identification and characterization of JAZ gene family in upland cotton (Gossypium hirsutum). Sci. Rep. 2017, 7, 2788. [Google Scholar] [CrossRef] [Green Version]
  36. Xie, S.; Cui, L.; Lei, X.; Yang, G.; Li, J.; Nie, X.; Ji, W. The TIFY gene family in wheat and its progenitors: Genome-Wide iIdentification, evolution and expression analysis. Curr. Genom. 2019, 20, 371–388. [Google Scholar] [CrossRef]
  37. Zhang, Y.; Gao, M.; Singer, S.D.; Fei, Z.; Wang, H.; Wang, X. Genome-wide identification and analysis of the TIFY gene family in grape. PLoS ONE 2012, 7, e44465. [Google Scholar] [CrossRef] [Green Version]
  38. Zhu, D.; Bai, X.; Luo, X.; Chen, Q.; Cai, H.; Ji, W.; Zhu, Y. Identification of wild soybean (Glycine soja) TIFY family genes and their expression profiling analysis under bicarbonate stress. Plant Cell Rep. 2013, 32, 263–272. [Google Scholar] [CrossRef]
  39. Dai, Z.; Dong, S.; Miao, H.; Liu, X.; Han, J.; Li, C.; Gu, X.; Zhang, S. Genome-wide identification of TIFY genes and their response to various pathogen infections in cucumber (Cucumis sativus L.). Sci. Hortic. 2022, 295, 110814. [Google Scholar] [CrossRef]
  40. He, X.; Kang, Y.; Li, W.; Liu, W.; Xie, P.; Liao, L.; Huang, L.; Yao, M.; Qian, L.; Liu, Z.; et al. Genome-wide identification and functional analysis of the TIFY gene family in the response to multiple stresses in Brassica napus L. BMC Genom. 2020, 21, 736. [Google Scholar] [CrossRef]
  41. Liu, X.; Zhao, C.; Yang, L.; Zhang, Y.; Wang, Y.; Fang, Z.; Lv, H. Genome-wide identification, expression profile of the TIFY gene family in Brassica oleracea var. capitata, and their divergent response to various pathogen infections and phytohormone treatments. Genes 2020, 11, 127. [Google Scholar]
  42. Gopal, S.; Jong-In, P.; Abdul, K.M.; Ill-Sup, N. A genome-wide analysis reveals stress and hormone responsive patterns of TIFY family genes in Brassica rapa. Front. Plant Sci. 2016, 7, 936. [Google Scholar]
  43. Bai, Y.; Meng, Y.; Huang, D.; Qi, Y.; Chen, M. Origin and evolutionary analysis of the plant-specific TIFY transcription factor family. Genomics 2011, 98, 128–136. [Google Scholar] [CrossRef] [PubMed]
  44. Chalhoub, B.; Denoeud, F.; Liu, S.; Parkin, I.A.; Tang, H.; Wang, X.; Chiquet, J.; Belcram, H.; Tong, C.; Samans, B.; et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 2014, 345, 950–953. [Google Scholar] [CrossRef] [Green Version]
  45. Yang, N.; Liu, J.; Gao, Q.; Gui, S.; Chen, L.; Yang, L.; Huang, J.; Deng, T.; Luo, J.; He, L.; et al. Genome assembly of a tropical maize inbred line provides insights into structural variation and crop improvement. Nat. Genet. 2019, 51, 1052–1059. [Google Scholar] [CrossRef] [Green Version]
  46. Liu, Y.; Du, H.; Li, P.; Shen, Y.; Peng, H.; Liu, S.; Zhou, G.A.; Zhang, H.; Liu, Z.; Shi, M.; et al. Pan-genome of wild and cultivated soybeans. Cell 2020, 182, 162–176.e113. [Google Scholar] [CrossRef]
  47. Sun, H.J.; Uchii, S.; Watanabe, S.; Ezura, H. A highly efficient transformation protocol for Micro-Tom, a model cultivar for tomato functional genomics. Plant Cell Physiol. 2006, 47, 426–431. [Google Scholar] [CrossRef]
  48. Guan, J.; Xu, Y.; Yu, Y.; Fu, J.; Ren, F.; Guo, J.; Zhao, J.; Jiang, Q.; Wei, J.; Xie, H. Genome structure variation analyses of peach reveal population dynamics and a 1.67 Mb causal inversion for fruit shape. Genome Biol. 2021, 22, 13. [Google Scholar] [CrossRef]
  49. Guo, J.; Cao, K.; Deng, C.; Li, Y.; Zhu, G.; Fang, W.; Chen, C.; Wang, X.; Wu, J.; Guan, L.; et al. An integrated peach genome structural variation map uncovers genes associated with fruit traits. Genome Biol. 2020, 21, 258. [Google Scholar] [CrossRef]
  50. Zhou, H.; Ma, R.; Gao, L.; Zhang, J.; Zhang, A.; Zhang, X.; Ren, F.; Zhang, W.; Liao, L.; Yang, Q.; et al. A 1.7-Mb chromosomal inversion downstream of a PpOFP1 gene is responsible for flat fruit shape in peach. Plant Biotechnol. J. 2021, 19, 192–205. [Google Scholar] [CrossRef]
  51. Li, S.; Xu, B.; Niu, X.; Lu, X.; Cheng, J.; Zhou, M.; Hooykaas, P.J.J. JAZ8 interacts with VirE3 Attenuating Agrobacterium mediated root Tumorigenesis. Front. Plant Sci. 2021, 12, 685533. [Google Scholar] [CrossRef] [PubMed]
  52. Ingle, R.A.; Stoker, C.; Stone, W.; Adams, N.; Smith, R.; Grant, M.; Carre, I.; Roden, L.C.; Denby, K.J. Jasmonate signalling drives time-of-day differences in susceptibility of Arabidopsis to the fungal pathogen Botrytis cinerea. Plant J. 2015, 84, 937–948. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Li, L.; Luo, Y.; Chen, B.; Xu, K.; Zhang, F.; Li, H.; Huang, Q.; Xiao, X.; Zhang, T.; Hu, J.; et al. A Genome-wide association study reveals new loci for resistance to clubroot disease in Brassica napus. Front. Plant Sci. 2016, 7, 1483. [Google Scholar] [CrossRef] [PubMed]
  54. Song, J.M.; Liu, D.X.; Xie, W.Z.; Yang, Z.; Guo, L.; Liu, K.; Yang, Q.Y.; Chen, L.L. BnPIR: Brassica napus pan-genome information resource for 1689 accessions. Plant Biotechnol. J. 2020, 19, 412–414. [Google Scholar] [CrossRef] [PubMed]
  55. Cai, X.; Wu, J.; Liang, J.; Lin, R.; Zhang, K.; Cheng, F.; Wang, X. Improved Brassica oleracea JZS assembly reveals significant changing of LTR-RT dynamics in different morphotypes. Theor. Appl. Genet. 2020, 133, 13. [Google Scholar] [CrossRef]
  56. Zhang, L.; Cai, X.; Wu, J.; Liu, M.; Wang, X. Improved Brassica rapa reference genome by single-molecule sequencing and chromosome conformation capture technologies. Hortic. Res. 2018, 5, 11. [Google Scholar] [CrossRef] [Green Version]
  57. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An iIntegrative tToolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 9. [Google Scholar] [CrossRef]
  58. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  59. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef]
  60. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2 (-Delta Delta C(T)). Method 2001, 25, 402–408. [Google Scholar] [CrossRef]
  61. Layer, R.M.; Hall, I.M.; Quinlan, A.R. LUMPY: A probabilistic framework for structural variant discovery. Genome Biol. 2014, 15, R84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Chen, X.; Schulz-Trieglaff, O.; Shaw, R.; Barnes, B.; Schlesinger, F.; Källberg, M.; Cox, A.J.; Kruglyak, S.; Saunders, C.T. Manta: Rapid detection of structural variants and indels for germline and cancer sequencing applications. Bioinformatics 2016, 32, 1220–1222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Jeffares, D.C.; Jolly, C.; Hoti, M.; Speed, D.; Sedlazeck, F.J. Transient structural variations have strong effects on quantitative traits and reproductive isolation in fission yeast. Nat. Commun. 2017, 8, 14061. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Phylogenetic analysis of 114 JAZ proteins from Arabidopsis, B. napus, B. rapa, and B. oleracea. The JAZ proteins from B. napus (51), B. rapa (25), B. oleracea (26), and Arabidopsis (12) are marked with stars, hollow circles, solid circles, and triangles, respectively. The number of bootstrap values was based on 1000 iterations. The hollow and solid stars represent genes located in the A and C subgenomes in B. napus, respectively. The neighbor-joining (NJ) tree shows six groups marked by Roman numerals and different colors: group I (purple) homologous to AtJAZ7/AtJAZ8, group II (yellow) homologous to AtJAZ10, group III (red) homologous to AtJAZ3/AtJAZ4/AtJAZ9, group IV (blue) homologous to AtJAZ11/AtJAZ12, group V (green) homologous to AtJAZ5/AtJAZ6 and group VI (orange) homologous to AtJAZ1/AtJAZ2.
Figure 1. Phylogenetic analysis of 114 JAZ proteins from Arabidopsis, B. napus, B. rapa, and B. oleracea. The JAZ proteins from B. napus (51), B. rapa (25), B. oleracea (26), and Arabidopsis (12) are marked with stars, hollow circles, solid circles, and triangles, respectively. The number of bootstrap values was based on 1000 iterations. The hollow and solid stars represent genes located in the A and C subgenomes in B. napus, respectively. The neighbor-joining (NJ) tree shows six groups marked by Roman numerals and different colors: group I (purple) homologous to AtJAZ7/AtJAZ8, group II (yellow) homologous to AtJAZ10, group III (red) homologous to AtJAZ3/AtJAZ4/AtJAZ9, group IV (blue) homologous to AtJAZ11/AtJAZ12, group V (green) homologous to AtJAZ5/AtJAZ6 and group VI (orange) homologous to AtJAZ1/AtJAZ2.
Ijms 23 12862 g001
Figure 2. Phylogenetic analysis (A), Conserved motif (B) and gene structure (C) of 114 JAZ proteins from Arabidopsis, B. napus, B. rapa, and B. oleracea.
Figure 2. Phylogenetic analysis (A), Conserved motif (B) and gene structure (C) of 114 JAZ proteins from Arabidopsis, B. napus, B. rapa, and B. oleracea.
Ijms 23 12862 g002
Figure 3. Chromosomal mapping of 51 BnaJAZs. The length of the chromosome is represented in Mb. The chromosome number is presented on the bottom of each chromosome, and the names of BnaJAZs were at their right.
Figure 3. Chromosomal mapping of 51 BnaJAZs. The length of the chromosome is represented in Mb. The chromosome number is presented on the bottom of each chromosome, and the names of BnaJAZs were at their right.
Ijms 23 12862 g003
Figure 4. Expression patterns of 51 BnaJAZs under abiotic stress. The expression level of BnaJAZs in leaves and roots after 1 h and 3 h of treatment with various hormones and adversities. The homologous genes of BnaJAZs in Arabidopsis are displayed at the right of BnaJAZs. The type and time of treatments are shown at the bottom of the expressional heat map. The tissues of the samples were on the top of the expressional heat map.
Figure 4. Expression patterns of 51 BnaJAZs under abiotic stress. The expression level of BnaJAZs in leaves and roots after 1 h and 3 h of treatment with various hormones and adversities. The homologous genes of BnaJAZs in Arabidopsis are displayed at the right of BnaJAZs. The type and time of treatments are shown at the bottom of the expressional heat map. The tissues of the samples were on the top of the expressional heat map.
Ijms 23 12862 g004
Figure 5. Expression profiles of BnaJAZs belonging to group I (A) and group IV (B) after infection of P. brassicae. The columns present the relative expressional foldchange compared with the mock values in the root of resistant and susceptible accessions. The relative expression levels were analyzed by qRT-PCR. The characters on the upper panel were presented with the p-values of a t-test. * presented the significant difference at a 95% confidence interval.
Figure 5. Expression profiles of BnaJAZs belonging to group I (A) and group IV (B) after infection of P. brassicae. The columns present the relative expressional foldchange compared with the mock values in the root of resistant and susceptible accessions. The relative expression levels were analyzed by qRT-PCR. The characters on the upper panel were presented with the p-values of a t-test. * presented the significant difference at a 95% confidence interval.
Ijms 23 12862 g005
Figure 6. Structural variation and analysis of significance of the differences in resistance. The translocation was detected in BnaJAZ8.C03, which resulted in the differences in resistance to P. brassicae. The black opposite arrows represent the normal read pairs. The green opposite arrows represent the abnormal read pairs, one end of which was located on the BnaC03 chromosome, but the other end was derived from BnaA08 (A). The performance of resistant (“28669”) and susceptible (“YJ-8”) accessions after 42 days infected by P. brassicae (B). The significance of resistance difference between accessions with and without the translocation in population. Haplotype A represented the average disease index of accessions without the translocation (including 368 accessions), and the haplotype B represented the average disease index of accessions with translocation (including 37 accessions). The p-value of the t-test is followed by the asterisk (C).
Figure 6. Structural variation and analysis of significance of the differences in resistance. The translocation was detected in BnaJAZ8.C03, which resulted in the differences in resistance to P. brassicae. The black opposite arrows represent the normal read pairs. The green opposite arrows represent the abnormal read pairs, one end of which was located on the BnaC03 chromosome, but the other end was derived from BnaA08 (A). The performance of resistant (“28669”) and susceptible (“YJ-8”) accessions after 42 days infected by P. brassicae (B). The significance of resistance difference between accessions with and without the translocation in population. Haplotype A represented the average disease index of accessions without the translocation (including 368 accessions), and the haplotype B represented the average disease index of accessions with translocation (including 37 accessions). The p-value of the t-test is followed by the asterisk (C).
Ijms 23 12862 g006
Figure 7. Interaction network of target JAZ proteins. Network nodes represent proteins interactive with JAZ8 directly or indirectly in Arabidopsis. Edges present protein–protein associations. The red gene names at the side of the interactive proteins were the homologous genes of interactive proteins in B. napus.
Figure 7. Interaction network of target JAZ proteins. Network nodes represent proteins interactive with JAZ8 directly or indirectly in Arabidopsis. Edges present protein–protein associations. The red gene names at the side of the interactive proteins were the homologous genes of interactive proteins in B. napus.
Ijms 23 12862 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, L.; Ji, G.; Guan, W.; Qian, F.; Li, H.; Cai, G.; Wu, X. Genome-Wide Identification and Variation Analysis of JAZ Family Reveals BnaJAZ8.C03 Involved in the Resistance to Plasmodiophora brassicae in Brassica napus. Int. J. Mol. Sci. 2022, 23, 12862. https://doi.org/10.3390/ijms232112862

AMA Style

Li L, Ji G, Guan W, Qian F, Li H, Cai G, Wu X. Genome-Wide Identification and Variation Analysis of JAZ Family Reveals BnaJAZ8.C03 Involved in the Resistance to Plasmodiophora brassicae in Brassica napus. International Journal of Molecular Sciences. 2022; 23(21):12862. https://doi.org/10.3390/ijms232112862

Chicago/Turabian Style

Li, Lixia, Gaoxiang Ji, Wenjie Guan, Fang Qian, Hao Li, Guangqin Cai, and Xiaoming Wu. 2022. "Genome-Wide Identification and Variation Analysis of JAZ Family Reveals BnaJAZ8.C03 Involved in the Resistance to Plasmodiophora brassicae in Brassica napus" International Journal of Molecular Sciences 23, no. 21: 12862. https://doi.org/10.3390/ijms232112862

APA Style

Li, L., Ji, G., Guan, W., Qian, F., Li, H., Cai, G., & Wu, X. (2022). Genome-Wide Identification and Variation Analysis of JAZ Family Reveals BnaJAZ8.C03 Involved in the Resistance to Plasmodiophora brassicae in Brassica napus. International Journal of Molecular Sciences, 23(21), 12862. https://doi.org/10.3390/ijms232112862

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop