Original Research ARTICLE
A Genome-Wide Analysis Reveals Stress and Hormone Responsive Patterns of TIFY Family Genes in Brassica rapa
- Department of Horticulture, Sunchon National University, Suncheon, South Korea
The TIFY family is a plant-specific group of proteins with a diversity of functions and includes four subfamilies, viz. ZML, TIFY, PPD, and JASMONATE ZIM-domain (JAZ) proteins. TIFY family members, particularly JAZ subfamily proteins, play roles in biological processes such as development and stress and hormone responses in Arabidopsis, rice, chickpea, and grape. However, there is no information about this family in any Brassica crop. This study identifies 36 TIFY genes in Brassica rapa, an economically important crop species in the Brassicaceae. An extensive in silico analysis of phylogenetic grouping, protein motif organization and intron-exon distribution confirmed that there are four subfamilies of BrTIFY proteins. Out of 36 BrTIFY genes, we identified 21 in the JAZ subfamily, seven in the TIFY subfamily, six in ZML and two in PPD. Extensive expression profiling of 21 BrTIFY JAZs in various tissues, especially in floral organs and at different flower growth stages revealed constitutive expression patterns, which suggest that BrTIFY JAZ genes are important during growth and development of B. rapa flowers. A protein interaction network analysis also pointed to association of these proteins with fertility and defense processes of B. rapa. Using a low temperature-treated whole-genome microarray data set, most of the JAZ genes were found to have variable transcript abundance between the contrasting inbred lines Chiifu and Kenshin of B. rapa. Subsequently, the expression of all 21 BrTIFY JAZs in response to cold stress was characterized in the same two lines via qPCR, demonstrating that nine genes were up-regulated. Importantly, the BrTIFY JAZs showed strong and differential expression upon JA treatment, pointing to their probable involvement in JA-mediated growth regulatory functions, especially during flower development and stress responses. Additionally, BrTIFY JAZs were induced in response to salt, drought, Fusarium, ABA, and SA treatments, and six genes (BrTIFY3a, 3b, 6a, 9a, 9b, and 9c) were identified to have co-responsive expression patterns. The extensive annotation and transcriptome profiling reported in this study will be useful for understanding the involvement of TIFY genes in stress resistance and different developmental functions, which ultimately provides the basis for functional characterization and exploitation of the candidate TIFY genes for genetic engineering of B. rapa.
The TIFY gene family is plant-specific and encodes putative transcription factors (TFs) characterized by a highly conserved motif (TIF[F/Y]XG) positioned within an approximately 28-amino-acid TIFY domain (Vanholme et al., 2007; Bai et al., 2011). ZIM (ZINC-FINGER EXPRESSED IN INFLORESCENCE MERISTEM) from Arabidopsis was the first identified member of TIFY gene family (Nishii et al., 2000). Based on phylogenetic and structural analyses, genes of this family can be divided into four groups. The TIFY subfamily proteins contain only the TIFY domain (Staswick, 2008), whereas ZIM and ZIM-like (ZML) proteins, together forming the ZML subfamily, also contain both a C2C2-GATA zinc-finger DNA-binding domain and a CCT domain (CONSTANS, CO-like, TOC1). The Jasmonate ZIM-domain (JAZ) subfamily proteins contain a conserved Jas motif of approximately 27 amino acids near the C-terminus in addition to TIFY domains. These Jas sequences have similarity to the N-terminal portion of the CCT domain, with the characteristic motif SLX2FX2KRX2RX5PY (Staswick, 2008; Chung et al., 2009). By contrast, PEAPOD (PPD) subfamily proteins possess a characteristic N-terminal PPD domain and a modified Jas motif that lacks the conserved PY (Proline-Tyrosine) at their C-terminal region (Chung et al., 2009).
To date, most of the functionally validated TIFY genes are from Arabidopsis. So far, 18 TIFY genes in Arabidopsis, 22 in rice, 22 in poplar, 19 in grape, 30 in apple, and 34 in Glycine soja (wild soybean) have been identified (Vanholme et al., 2007; Ye et al., 2009; Zhang et al., 2012; Zhu et al., 2013; Li et al., 2015). TIFY genes are involved in the regulation of diverse aspects of plant development and responses to abiotic stresses and phytohormone treatments. For instance, overexpression of ZIM/AtTIFY1 in Arabidopsis produced elongated petioles and hypocotyls (Shikata et al., 2004), while AtTIFY4a (PPD1) and AtTIFY4b (PPD2) have been reported to be involved in the synchronization of leaf growth (White, 2006). Among the TIFY family proteins, JAZ subfamily members are the best characterized and play pivotal roles in the jasmonic acid (JA) pathway (Chini et al., 2007; Thines et al., 2007; Yan et al., 2007). JA and its bioactive derivatives have been identified as prime factors in plant responses to both biotic (e.g., wounding, pathogen infection, and insect attack) and abiotic (e.g., drought, low temperature, salinity, ozone exposure) stresses (Devoto and Turner, 2005). AtTIFY10b/JAZ1 is an important player in JA signal transduction. To stop their repression of JA-response genes, JAZ1 is degraded in the presence of JA through the SCFCOI1-dependent 26S proteasome pathway (Thines et al., 2007). JAZ proteins in different studies have been identified as repressors of JA-regulated transcription in Arabidopsis. Specifically, in the absence of JA, JAZ proteins interact with and repress the basic helix-loop-helix (bHLH) TFs MYC2 and MYC3. These TFs bind directly to DNA sequences and regulate the expression of downstream JA-responsive genes (Chini et al., 2007; Thines et al., 2007). JAZ genes show differential induction in response to JA treatment, herbivory, wounding, Pseudomonas syringae infection, and environmental stresses such as drought, low temperature and salinity (Chung et al., 2008; Ye et al., 2009; Demianski et al., 2012). JAZ genes have also been induced in response to different nutritional stress like phosphorus and bicarbonate (Aparicio-Fabre et al., 2013; Zhu et al., 2013). A more recent report on the involvement of JAZ genes against nutritional deficiency revealed their specificity with different nutrient stimuli (macro and micro nutrients) in rice and chickpea, where majority of the JAZ genes were induced by potassium (K) deficiency (Singh et al., 2015).
Brassica is a diversified crop genus that provides oil, vegetables, condiments, dietary fiber, and vitamin C (Talalay and Fahey, 2001). Brassica rapa, one of the most important members of the Brassica genus, has several subspecies and is a widely used vegetable in Asia. Over the course of evolution, the Brassicaceae genome experienced two whole-genome duplication (WGD) events as Arabidopsis beta and Arabidopsis alpha (Bowers et al., 2003). WGD results in gene duplication that is typically followed by significant gene loss (Lee et al., 2013). Evolutionary data indicate that B. rapa ssp. pekinensis is closely related to Arabidopsis and that since its divergence from Arabidopsis 13 to 17 Mya (Million years ago) it has experienced a whole-genome triplication (Town et al., 2006; Wang et al., 2011). This lineage thus provides a model for studying gene fractionation (Tang and Lyons, 2012; Kim et al., 2013). B. rapa as a representative of the Brassica “A” genome has been used as the model species and therefore, was selected for genome sequencing [Brassica Genome Gateway (http://brassica.nbi.ac.uk/); The Korea Brassica rapa Genome Project (http://www.brassica-rapa.org/BRGP/index.jsp)]. This species has a relatively small genome [approximately 529 megabase pairs (Mbp)] compared to other Brassica species and has proved a useful model for studying polyploidy.
Though, TIFY family genes have been investigated in Arabidopsis, which is a close relative of B. rapa, we set out to study this gene family in the relatively complex and diverse genome of B. rapa. The recent completion of the B. rapa Chiifu-401-42 genome sequencing (Wang et al., 2011) enabled us to carry out genome-wide research on TIFY proteins in B. rapa. We identified 36 TIFY genes in B. rapa. The putative BrTIFY genes were classified into four subfamilies and their gene structures, chromosomal locations, and potential gene duplicates were subsequently analyzed. From a low temperature-treated microarray data set, responses of all BrTIFY genes were investigated, comparing their expression between two contrasting inbred lines Chiifu and Kenshin of B. rapa. We further analyzed the responses of the BrTIFY JAZ genes to cold, salt, drought, Fusarium oxysporum, ABA, JA, and SA via quantitative RT-PCR. In an organ expression study, all the BrTIFY JAZ genes showed constitutive expression during the flower development stages and these results highlight the need for further study into the role(s) of BrTIFY JAZ genes in the development of B. rapa flowers.
Materials and Methods
Identification of B. rapa TIFY Genes
A SWISSPROT search on the Brassica database (BRAD) was conducted using keyword “TIFY” and a blast search using characteristic TIFY domain sequences was also employed in the BRAD to predict B. rapa TIFY proteins (V 1.5.tar.gz; http://brassicadb.org/brad/; Cheng et al., 2011). CDS and Protein sequences of the retrieved B. rapa TIFY genes were obtained from the Brassica database (http://brassicadb.org/brad/; Cheng et al., 2011). “TIFY” domains were confirmed using the web tool “SMART” from EMBL (http://smart.embl.de/smart/set_mode.cgi?GENOMIC=1) and by using the Basic Local Alignment Search Tool (BLAST) (http://www.ncbi.nlm.nih.gov/BLAST/) a protein homology study was conducted to confirm the putative TIFY genes in B. rapa.
Sequence Analysis of B. rapa TIFY Genes
A thorough investigation on the primary structures (molecular weight and isoelectric point) of the B. rapa TIFY proteins was carried out using protParam (http://expasy.org/tools/protparam.html). “PIR,” an online web tool (http://pir.georgetown.edu/pirwww/search/multialn.shtml) was used to align the TIFY protein sequences. Number of introns and exons was predicted using the Gene Structure Display Server (GSDS) web tool aligning the CDS with the genomic sequences (Guo et al., 2007). The subcellular location of TIFY proteins in B. rapa was determined using ProtComp 9.0 from Softberry (http://linux1.softberry.com/berry.phtml) and Blast2GO software (http://www.blast2go.de). For the analysis of the B. rapa TIFY protein motifs MEME software (Multiple Em for Motif Elicitation, V4.9.0) was employed (Bailey et al., 2006). A search on the MEME suit was executed to identify distinctive motifs with the following parameters: (1) width of optimum motif ≥6 and ≤ 200; (2) maximum number of motifs to identify = 10. Important protein motifs were drawn manually. Online tool STRING 10 (http://string-db.org/) was employed to predict the TIFY protein association and interactions in B. rapa.
Phylogenetic Analysis of B. rapa TIFY Proteins
TIFY protein sequences from B. rapa, rice, and Arabidopsis were aligned using ClustalX (Thompson et al., 1997). A phylogenetic tree was generated with the help of MEGA6.06 software using the Neighbor-Joining (NJ) algorithm (Tamura et al., 2013). Significance of the nodes was analyzed through bootstrapping with 1000 replicates. Pairwise gap deletion mode was used to assume that the divergent domains might contribute to the topology of the NJ tree.
Chromosomal Locations and Gene Duplication of BrTIFYs
A BLAST search (http://www.ncbi.nlm.nih.gov/BLAST/) was employed using all the BrTIFY genes against each other and gene duplications were considered when both the similarity and query coverage percentage of the candidate genes were >80% (Kong et al., 2013). Chromosomal positions of the candidate BrTIFY genes along with their subgenome information were obtained from the Brassica database (http://brassicadb.org/brad/; Cheng et al., 2011). Gene positions in the chromosomes and duplication lines among genes on the chromosomes were drawn manually. We aligned the duplicated BrTIFY CDS sequences using ClustalW (http://www.genome.jp/tools/clustalw/) to predict the divergences period. We estimated the synonymous rate (Ks), non-synonymous rate (Ka), and evolutionary constriction (Ka/Ks) between the duplicated pairs of B. rapa TIFY genes using codeml in the PAML package (Goldman and Yang, 1994). The formula T = Ks/2R was employed to calculate the divergence time, where T refers to divergence time, Ks refers to the synonymous substitutions per site, and R, for the rate of divergence of plant's nuclear genes. R-value was considered as 1.5 × 10−8 synonymous substitutions per site per year in case of dicotyledonous plants (Koch et al., 2000).
Collection and Preparation of Plant Material
For plant samples, we cultivated B. rapa “SUN-3061” at the Department of Horticulture, Sunchon National University, Korea. Fresh roots, stems, leaves, and flower buds of mature plants were harvested for the organ study. Collected plant parts were frozen immediately in liquid nitrogen, and stored at −80°C for RNA isolation. Two inbred lines of B. rapa ssp. pekinensis, Chiifu and Kenshin having two contrasting geographic origins were used for the cold stress treatment at 4°C. Chiifu originated in temperate regions and Kenshin originated in subtropical and tropical regions (Lee et al., 2013). Plants were cultured in a semi-solid medium maintaining an aseptic condition for 2 weeks. To induce the cold stress, plant samples (Chiifu and Kenshin) were transferred to the incubator at 4°C. For rest of the abiotic and biotic stress treatments, only B. rapa Chiifu was used. To impose the drought/desiccation stress plants were kept on the Whatmann 3 mm filter sheets. For salt, abscisic acid (ABA), jasmonic acid (JA), and salicylic acid (SA) treatments, plants were shifted to rectangular petri dishes (72 × 72 × 100 mm) containing 200 mM NaCl, 100 μM ABA, JA, and SA, respectively (Yanhui et al., 2006; Saha et al., 2015). Sampling of all the abiotic stress treated plants were done at 0, 1, 4, 12, 24, and 48 h.
For biotic stress treatment, a fungal inoculation by Fusarium oxysporum f.sp. conglutinans in B. rapa Chiifu were done according to Ahmed et al. (2013). Fungus- and mock-treated plants were sampled at 0, 3, 6 h, 1, 3, 6, 9, and 12 d. Sampling were done from the local (4th) and systemic (5th) leaves. Inoculations were done three times, and infection was verified by monitoring disease lesions on plant leaves. Stress treated leaf samples were collected and frozen immediately in liquid nitrogen and stored at −80°C until RNA isolation.
Microarray Analysis of BrTIFY Expression under Low Temperature Treatment
Low temperature-treated microarray data for 36 BrTIFY genes were recovered from the data of Jung et al. (2014). For that data, two contrasting B. rapa inbred lines Chiifu and Kenshin were treated with different cold and freezing temperatures viz. 4, 0, −2, and −4°C for 2 h. The transcript abundance values of the 36 BrTIFY genes were collected and analyzed further to generate a heat map using Cluster 3.0 and tree view software (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm) with the centroid clustering option.
Transcript Level Analysis of BrTIFY Genes
An AMV one-step RT-PCR kit (Takara, Japan) was used for RT-PCR analysis. B. rapa candidate TIFY gene specific primers were used for RT-PCR expression analysis (Supplementary Table 1). For internal control Actin primers from B. rapa (FJ969844) were used. For PCR analysis 50 ng cDNA from the plant and flower organs was used as templates in master mixes comprised of 20 pmol each primer, 150 μM each dNTP, 1.2 U Taq polymerase, 1x Taq polymerase buffer, and double-distilled H2O in a total volume of 20 μL in 0.5-mL PCR tubes. The samples were kept in a thermal cycler for PCR reaction to the following conditions: pre-denaturing at 94°C for 5 min, followed by 30 cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, and extension at 72°C for 45 s, with a final extension for 5 min at 72°C.
The real-time quantitative PCR was conducted using 1 μL cDNA in a 20-μL reaction volume employing iTaqTM SYBR® Green Super-mix with ROX (California, USA). List of specific primers used for real-time PCR has given in Supplementary Table 1. For real-time PCR analysis, prepared samples were subjected to the following conditions: 10 min at 95°C, followed by 40 cycles at 95°C for 20 s, 58°C for 20 s, and 72°C for 25 s. The fluorescence was determined observing the last step of each cycle, and three replications were used per sample. Light cycler® 96 SW 1.1 (Roche, Germany) was used for amplification detection and data analysis.
Identification of TIFY Family Genes in the B. rapa Genome
To identify all putative TIFY genes in B. rapa, we performed a search in the Swissprot annotations of the Brassica database (BRAD) using key word “TIFY.” (http://brassicadb.org/brad/). A set of 36 candidate TIFY genes was identified (Supplementary Table 2) and their protein and coding sequences (CDS) were retrieved (http://brassicadb.org/brad/). Besides, a blast search against characteristic TIFY domain sequences in the BRAD revealed 35 B. rapa TIFY proteins. We crosschecked all the TIFY protein sequences of both the sources and 36 B. rapa TIFY proteins were primarily predicted for further in silico analysis. A domain search using EMBL (http://smart.embl.de/smart/set_mode.cgi?GENOMIC=1) (“Pfam” web server option was also employed at the same time) with the corresponding B. rapa candidate TIFY protein sequences also confirmed 35 of them to contain the characteristic “TIFY” domain. The remaining candidate TIFY protein, which apparently lacks a “TIFY” domain, showed significant similarity to a TIFY protein of Arabidopsis (AT2G34600; JASMONATE-ZIM-DOMAIN PROTEIN 7) at the amino acid level (Supplementary Table 3). Notably, this number of TIFY genes in B. rapa (36) genome is the highest among all the crops have been reported to date (Vanholme et al., 2007; Ye et al., 2009; Zhang et al., 2012; Zhu et al., 2013; Li et al., 2015; Figure 1). The putative TIFY genes in B. rapa were named (BrTIFY1a-BrTIFY11e) based on the numbering system used in phylogenetic tree of Arabidopsis and rice (Supplementary Table 2).
Sequence Analysis of BrTIFY Genes
The predicted sizes of the 36 BrTIFYs ranged from 112 to 363 amino acids (aa), and the predicted isoelectric point varied from 4.48 to 9.98 (Table 1). Further analysis of the protein domain organization showed that among 35 TIFY domain-containing proteins, 32 had both a TIFY domain and a Jas motif (Staswick, 2008) of SLX2FX2KRX2RX5PY. Only five proteins were found to bear a TIFY domain along with the CCT motif and C2C2-GATA zing finger domain that are characteristic of the ZML subfamily (Figure 2, Supplementary Figure 1). Two TIFY proteins containing the Jas motif were lacking the conserved proline (P) and tyrosine (Y) at the C-termini (Chung et al., 2009), and contained a PPD domain, which indicated that they were likely PPD subfamily proteins. Generally, a TIFY protein containing a TIFY domain along with a Jas motif has been considered as a member of the JAZ subfamily. Jas motifs are found to be located in the C-terminal region of the JAZ proteins which are essential for the interaction with LRR domain of COI1 (Katsir et al., 2008b). A highly conserved KR(K/N/R)(E/D)R sequence in the C-terminal end of the Jas motif has been observed which has been reported essential for nuclear localization signal of JAZ proteins (Withers et al., 2012; Supplementary Figure 2). More specifically, we also found conserved arginine (R) within the Jas motif (Figure 2), which is important for JAZ nuclear localization and interaction with MYC2 (Withers et al., 2012). Besides, seven proteins were found to contain the TIFY domain motif only. Accordingly, these seven proteins were supposed to be TIFY subfamily proteins.
Figure 2. (A) Showing schematic representation of intron-exon distribution and motifs identified in B. rapa TIFY proteins. Intron–exon organization patterns of 36 TIFY genes are shown, along with their intron splicing patterns. On the right, different motifs (motif 1–5) are indicated by different colors, and the names of all members are shown on the left side of the figure, along with their phylogenetic relatedness, and (B) showing two motif logos (“TIFY” and “Jas”) with their amino acid conservation patterns.
To examine the variability among TIFY domains, we collected all the TIFY domain sequences and aligned them, revealing characteristic TIFY domains in different combinations such as TI[F/S]YXG, T[L/I]S[F/Y]XG, TVSYXG, VIFYXG, and TIFF[G/R]G (Figure 3). Specifically, all the TIFY (except BrTIFY5c) and PPD subfamily proteins and 14 out of 21 JAZ proteins shared a common sequence “TIFYXG” in the TIFY domain. Notably, all the phenylalanine (F) positions of the “TIFY” motif in the ZML proteins were replaced by serine (S). In addition, they showed inconsistent “TIFY” motif conservation as T[L/I]S[F/Y]XG in their domain sequences.
Figure 3. TIFY domain alignment showing conserved “TIFY[X]G” motif among 36 B. rapa TIFY proteins. Four subfamily proteins have been marked on the left of the figure.
The DNA sequences of the 36 BrTIFY genes were determined based on the B. rapa whole-genome sequence. Exon-intron analysis of 36 BrTIFY genes revealed that most of the genes contained more than one exon (up to 9) with fairly consistent splicing patterns within different groups (Figure 2). Specifically, genes from the ZML subfamily bore 7–8 exons with the splicing pattern “0/2.” Two PPD subfamily genes possessed the most exons (8–9) and had “1/0/2” splicing patterns. The JAZ subfamily members were identified to have 4–5 exons with a splicing pattern “1/2” in most cases. Finally, the TIFY subfamily genes exhibited the highest variability; they contained 1 to 6 exons with splicing pattern “0/1/2.” We have studied the subcellular localization of the 36 BrTIFY proteins using Protcomp 9.0 from Softberry and Blast2GO software. Consistent with being transcription factor family proteins, all the BrTIFYs were predicted to be localized in the nucleus (Table 1).
Phylogenetic Analysis of B. rapa TIFY Proteins
A phylogenetic tree was generated using the sequence alignments of the 36 full-length TIFY proteins to analyze the evolutionary relationships among the BrTIFY proteins and 18 TIFY proteins from Arabidopsis as well as 20 from rice. The 74 TIFY proteins were grouped into ten clades (C1-10) (Figure 4). Among the 10 clades, C1 and C3 were formed with the TIFY proteins containing only the TIFY domain, and these two clades were named as TIFY I and TIFY II. The TIFY I group included five proteins (BrTIFY5a, 5b, 5c, 5d, and 5e) and the TIFYII group had two proteins (BrTIFY8a and 8b). In C2, ZML proteins from B. rapa, Arabidopsis and rice were gathered together. Six BrTIFY proteins (BrTIFY1a, 1b, 1c, 2a, 2b, and 2c) were grouped with corresponding proteins from the Arabidopsis and rice ZML subfamily. Surprisingly, BrTIFY1c, which does not have a CCT or a GATA-Zn finger domain, was placed in the ZML subfamily. Significant sequence similarity outside of the CCT and GATA-Zn finger domain might explain its presence in this subfamily. Two PPD proteins (BrTIFY4a and 4b) from B. rapa formed C5. These two proteins contain partial Jas (lack conserved PY) motif and additionally possess a PPD domain in the N-terminal region. Due to Jas motifs in their C-terminal regions, these two proteins had close relationships with C4, which also contained proteins with Jas motifs and was designated as JAZ I. Notably, 21 BrTIFY JAZ proteins were found in six clades (C4, C6, C7, C8, C9, and C10, which we named as JAZ I, JAZ II, JAZ III, JAZ IV, JAZ V, and JAZ VI), indicating that there is a broader phylogenetic relationship within this subset of proteins. Specifically, 16 BrTIFY (BrTIFY3a, 3b, 6a, 6b, 7a, 7b, 7b, 9a, 9b, 9c, 10a, 10b, 10c, 10d, 10e, and 10f) proteins were grouped within four JAZ subgroups (JAZ I, JAZ II, JAZ III, and JAZ VI) with proteins from Arabidopsis and rice. The other 5 BrTIFY JAZ proteins (BrTIFY11a, 11b, 11c, 11d, and 11e) were grouped in C9 (JAZ V) with related proteins from Arabidopsis. A rice-specific JAZ group was formed with six OsTIFY JAZs in C8, which is designated as JAZ IV. It should be noted that in the phylogenetic tree, BrTIFY proteins were more commonly present with Arabidopsis TIFY proteins in a sister clade rather than forming clades with rice TIFY proteins. This signifies that the majority of B. rapa TIFY proteins are more closely related to those of Arabidopsis than those of rice, which is consistent with the fact that both B. rapa and Arabidopsis are eudicots and diverged more recently from a common ancestor than from the lineage that leads to monocots. Some interesting evolutionary relationships were observed from the above phylogenetic classification. For example, there was an approximately 1:3 ratio in the number of Arabidopsis TIFY homologs to those in B. rapa. Comparative genomic analysis has confirmed that B. rapa has experienced genome triplication since its divergence from Arabidopsis (Song et al., 2014).
Figure 4. Phylogenetic relatedness and subgroup designations of TIFY proteins from B. rapa, Arabidopsis and rice. Maximum-likelihood tree showing the relationships among 36 TIFY proteins from B. rapa, 18 from Arabidopsis and 20 from rice. The numbers on the branches indicate bootstrap support values from 1000 replications.
Chromosomal Location of TIFY Genes and Their Expansion Patterns in B. rapa
We retrieved the physical locations of the 36 TIFY genes on the 10 chromosomes of B. rapa (Figure 5). The highest number of TIFY genes was found on chromosome 8 (6 genes). Chromosomes 1, 2, and 7 had 5 TIFY genes each and chromosome 4 contained only 1 TIFY gene. We also predicted the proportions of BrTIFY genes among three subgenomes of B. rapa (LF, MF1, and MF2). The most BrTIFY genes (15) were found in the LF subgenome, whereas 13 and 8 genes were found in the MF1 and MF2 subgenomes, respectively. Duplication analysis revealed that 22 out of 36 BrTIFY genes (61.1%) were present in two or more copies. This gene duplication likely occurred as a result of whole genome duplication (WGD) or segmental duplications (black dotted lines in Figure 5); no evidence of tandem duplication was observed in any of the chromosomes. Higher frequencies of segmental duplications generated many homologs of TIFY genes along all chromosomes of B. rapa (black dotted lines in Figure 5). Moreover, we estimated the Ka/Ks ratios as less than one (1) for all the duplicated TIFY gene pairs in B. rapa genome which indicate about their purifying selections. The divergence time for duplicated BrTIFY genes revealed that it took 2.5–13.05 MYA (Table 2). The divergence period of six TIFY genes (BrTIFY6a-BrTIFY6b, BrTIFY10a-BrTIFY10b, BrTIFY10a-BrTIFY10c, BrTIFY10b-BrTIFY10c, BrTIFY11a-BrTIFY11b, and BrTIFY11c-BrTIFY11e) ranged from 5.2 to 7.6 MYA, indicating that divergence of these genes took place during the Brassica triplication events (5~9 MYA) (Woodhouse et al., 2014). Five gene pairs (BrTIFY3a-BrTIFY3b, BrTIFY4a-BrTIFY4b, BrTIFY5a-BrTIFY5b, BrTIFY5d-BrTIFY5e, and BrTIFY7b-BrTIFY7c) ranged their divergence time from 11.14 to 13.05 MYA, which indicate that these duplication events occurred during the divergence of B. rapa from Arabidopsis (9.6–16.1 MYA) (Wang et al., 2011). Finally, seven pairs of duplicated BrTIFY genes (BrTIFY2b-BrTIFY2c, BrTIFY8a-BrTIFY8b, BrTIFY9a-BrTIFY9b, BrTIFY9a-BrTIFY9c, BrTIFY9b-BrTIFY9c, BrTIFY10d-BrTIFY10e, and BrTIFY10d-BrTIFY10f) were identified as they diverged recently (2.5–4.22 MYA).
Figure 5. Distribution of TIFY genes on the 10 chromosomes and three subgenomes of B. rapa. Different colored blocks on chromosomes are representing three subgenomes. TIFY genes are shown on the right of each chromosome. The positive (+) and negative (−) signs following each gene represent forward and reverse orientation of the respective gene, respectively, on the chromosome or subgenome. Genes are in duplicated segments of genome are joined by black dotted lines. Gene positions and chromosome size can be measured using the scale on the left of the figure.
Expression Patterns of TIFY Genes in B. rapa
Expression Profiles of BrTIFY JAZ Genes in Various Organs
To gain insight into the potential developmental roles of B. rapa TIFY genes, we selected the JAZ genes for further analysis, as these genes make up the largest portion of the TIFY genes and have been reported to be involved in growth regulatory functions related to organ, especially flower development (Chini et al., 2007; Chung et al., 2009; Song et al., 2011; Hori et al., 2014; Yuan and Zhang, 2015). First, we conducted semi-quantitative RT-PCR analysis using specific primers with equal amounts of cDNA templates prepared from the mRNA of roots, stems, leaves, and flower buds of healthy, unstressed B. rapa plants to study the organ-specific expression patterns of JAZ genes. The genes showed differential expression patterns, with only four genes (BrTIFY3a, 3b, 7c, and 10b) expressed in all four tested organs (Figure 6). All 21 JAZ subfamily genes were predominantly expressed in the flower buds of B. rapa, which indicates that genes from this subfamily might have important roles related to growth and development of B. rapa flowers.
Figure 6. RT-PCR expression of 21 JAZ subfamily genes in different tissues of B. rapa. Data from lane 1–14 showing expression in roots (R), stems (S), leaves (L), flower buds (Fb), sepal (SP), petal (Pt), stamen (St), pistil (Pi), and six flower growth stages (young to mature buds are marked S1–S6 at the top of the figure) of B. rapa. Name of genes are indicated on the left of the figure.
Subsequently, we investigated the expression of the BrTIFY JAZ genes in four floral tissues viz. sepal, petal, stamen, and pistil of B. rapa flowers. Other than BrTIFY11a, which had no expression in any of the floral tissues, all of these genes exhibited widespread expression in most of the floral organs. The omnipresent expression of these genes in the floral tissues of B. rapa further indicates that these genes might be important factors for flower development in B. rapa. For further verification, we observed the expression patterns of all 21 JAZ genes throughout the development of B. rapa flowers (young to mature bud stages; S1–S6; Figure 6). We found all of these genes to be expressed in different growth stages of flower buds in two patterns. Four genes (BrTIFY3a, 3b, 7c, and 10b) expressed ubiquitously with equal transcript abundance levels from young to the mature bud stage. Surprisingly, BrTIFY11a also showed widespread expression in all the flower bud development stages. Besides, we observed that the rest of the JAZ genes (16 genes) showed gradual increments of their transcript levels from young to the mature bud stage. These genes might play greater roles in the mature buds than in the young buds.
Microarray Expression Analysis of BrTIFY Genes in Response to Cold and Freezing Treatments
We studied microarray expression data of the 36 TIFY genes using a previously published data set, wherein two contrasting inbred lines of B. rapa, Chiifu and Kenshin, were treated with cold and freezing stresses (4, 0, −2, and −4°C) for 2 h (Jung et al., 2014; Figure 7). Chiifu and Kenshin have contrasting geographic origins; Chiifu originated in temperate regions and Kenshin originated in subtropical and tropical regions. Hence, they are expected to respond differently to cold and freezing temperatures. In fact, Kenshin experiences severe damage upon exposure to low temperatures, while Chiifu does not exhibit such damage (Dong et al., 2014). To investigate expression patterns, we constructed a heat map using the differential transcript values for the TIFY genes in response to cold and freezing treatments in the two contrasting inbred lines of B. rapa. We divided the expression patterns of the 36 TIFY genes into five clusters (Cl1 to Cl5). In Cl1, eight TIFY genes [two PPD (BrTIFY4a, 4b), two TIFY (BrTIFY8a, 8b), two ZML (BrTIFY1a, 1c), and two JAZ (BrTIFY10f, 11d)] showed higher expression in response to both cold and freezing temperature in Kenshin than in Chiifu. The eight TIFY genes (six JAZ, one TIFY and one ZML genes) in Cl2 showed higher expression during the normal temperature (22°C) in both lines. Notably, three Chiifu JAZ genes (BrTIFY3a, 3b, and 7b) in this cluster showed the highest expression in response to freezing temperature (−4°C). In Cl3, seven genes formed a group, four of which were from TIFY (BrTIFY5a, 5c, 5d, and 5e) and the other three (BrTIFY9b, 11c, and 11e) were from JAZ subfamily. Generally, genes of this cluster showed very low expression upon cold and freezing in the two lines. In Cl4, nine TIFY genes in B. rapa Chiifu exhibited the highest transcript abundance in response to cold and freezing temperatures. Surprisingly, among these nine genes, seven were from the JAZ subfamily (BrTIFY6a, 6b, 7a, 9c, 10b, 10c, and 11b). Finally, in Cl5, four TIFY genes (three JAZ; BrIFY11a, 11d, 11e, and one ZML; BrTIFY1b) showed higher transcript abundance upon cold and freezing in both of the B. rapa lines. As a whole, we found that most of the JAZ genes in Chiifu were more highly induced by cold and freezing than those in Kenshin. Thus, we speculate that the set of TIFY genes exhibiting higher expression in Chiifu upon cold and freezing treatments might play roles in the cold and freezing tolerance of B. rapa Chiifu. Conversely, the TIFY genes that showed higher expression in Kenshin upon cold and freezing temperature might be factors that make Kenshin cold and freezing susceptible.
Figure 7. Microarray expression represented as heat map and hierarchical clustering of 36 TIFY genes against different low temperatures such as control (C1 and K1), 4°C (C2 and K2), 0°C (C3 and K3), −2°C (C4 and K4), and −4°C (C5 and K5) in B. rapa. Here, “C” and “K” stand for “Chiifu” and “Kenshin” two inbred lines of B. rapa, respectively. Expression clusters are shown in the left (Cl1–Cl5) and gene name against each expression has mentioned at right side. Color bars with values at the base represent differential expression in microarray.
Expression Analysis of JAZ Genes under Abiotic Stress Conditions in B. rapa
Different abiotic stresses such as low temperature, drought and salinity adversely affect plant growth and development. In the microarray expression data, we found that most of the TIFY genes, specifically the JAZ genes showed comparatively higher expression under low temperature conditions. Therefore, to further elucidate and to validate the responses of the JAZ genes to cold stress we used two inbred lines of B. rapa, Chiifu and Kenshin for qPCR expression analysis (Figures 8A,B). Additionally, we examined responses of these genes against salt and drought treatment (Figures 8C,D). We studied all the expressions on a time course basis (0, 1, 4, 12, 24, and 48 h).
Figure 8. Showing RT-qPCR expression analyses of 20 B. rapa JAZ genes in response to (A,B) cold, (C) salt, and (D) drought stresses (0–48 h). The error bars represent the standard error of the means of three independent replicates. Values denoted by the same letter did not differ significantly at P < 0.05 according to Duncan's multiple range tests.
In Chiifu, most of the JAZ genes were up-regulated in response to cold treatment. For example, BrTIFY3a, 6a, 6b, 7b, 9b, 9c, 10c, 10f, 11a, 11b, and 11c exhibited several fold up-regulation compared to the control under cold treatment. The highest up-regulation was shown by BrTIFY11b (>20-fold). BrTIFY6a and 10c showed gradually increasing expression from 0 to 48 h, when they up-regulated by more than 8- and 15-fold, respectively. BrTIFY3a, 6a, 6b, 7b, 7c, 9a, 9b, 9c, 10c, 10e, 11a, 11b, 11c, and 11d showed down-regulation at 1 h just after being induced by cold treatment. After that, from 4 to 48 h they showed differential up-and down-regulation during the cold period. In comparison with the microarray data, we also found that these genes also had differential expression against the cold and freezing temperatures. In general, the up-regulated BrTIFY JAZ genes showed their highest transcript levels at 48 h and exhibited differential expression patterns at other time points. We also found down-regulation of BrTIFY3b, 10a, 10b, 10d, and 11e throughout the cold treatment period and their down-regulation might be related with the cold resistance of Chiifu (Figure 8A). By contrast, in Kenshin, most of the JAZ genes were actually down-regulated or up-regulated to lower levels in response to cold stress. For example, BrTIFY 9b and 9c in Kenshin were down-regulated gradually from 0 to 48 h during cold stress. The JAZ genes highly up-regulated in Chiifu (BrTIFY 6a, 10c, and 11b) also showed up-regulation in Kenshin, but their levels of up-regulation were only about half as high (Figure 8B). These results suggest that the highly up-regulated JAZ genes in Chiifu might be related to the cold resistance of this line and, by contrast, the very low activity or no activity of these genes in response to cold stress in Kenshin might be responsible for the cold susceptibility of this line at differential levels.
Under salt stress in Chiifu, we identified some JAZ candidates that were up-regulated. In general, most of the Chiifu JAZ genes were highly induced at 4 h of treatment. And, at 1 h, most of these genes were induced only at low levels or were down-regulated. The highest up-regulation was exhibited by BrTIFY9b, which showed >100-fold up-regulation at 4 h. Two other members of the JAZ II (BrTIFY9a and 9c) also showed up-regulation, by >40- and >25-fold, respectively, compared to the reference. At the same time, other members of JAZ subfamily viz. BrTIFY3a, 3b, 6a, 6b, 10d, 10e, 10f, 11c, and 11e showed considerable and differential response to salt stress. BrTIFY3a, 6a, 7b, 7c, 9a, 9b, 9c, 10c, 10e, and 11b showed their highest up-regulation at 4 h which has been down-regulated a bit afterwards and remained static till 48 h. Besides, BrTIFY6b, 10d, 10f, 11a, 11c, and 11e were highly induced at 24/48 h of salt treatment. We found only BrTIFY11d to be down-regulated completely after salt treatment (Figure 8C).
In the case of exposure to drought stress, JAZ genes in Chiifu showed variable transcript abundance at different time points. Notably, 18 JAZ genes out of 20 were up-regulated in response to drought stress. Most of the genes started to up-regulate differentially from 4 h to onward. The highest induction was recorded for BrTIFY3a, which showed >15-fold up-regulation compared to the control. BrTIFY9b, 9c, 10c, and 10e also exhibited notable up-regulation in response to drought stress. Among the drought responsive JAZ genes BrTIFY3a, 3b, 6b, 9a, 9b, 9c, 10a, 10b, 10c, 10f, and 11e showed highest up-regulation during late hours of treatment that is from 24 to 48 h. Only two JAZ genes, viz. BrTIFY7c and 11d, showed complete down-regulation in response to drought treatment (Figure 8D).
It was interesting that most of the JAZ genes first showed up-regulation during early hours of salt and drought treatment followed by slight down regulation and then again up-regulated. The JAZ genes highly induced in response to salt and drought treatments might be good candidates for developing resistance against those stresses in B. rapa.
Expression Patterns of JAZ Genes in Response to Hormone Treatments in B. rapa
Phytohormones play important roles allowing plants to respond against stresses (Fujita et al., 2006; Huang et al., 2008). JAZ subfamily genes have been demonstrated to function in the JA signaling pathway in Arabidopsis (Chini et al., 2007; Thines et al., 2007). Here, we studied the responses of JAZ genes to JA, ABA, and SA in B. rapa Chiifu. Following treatment with JA, all JAZ genes in B. rapa were induced highly, where the highest up-regulation was recorded for BrTIFY9b (>800-fold), while 9a and 9c showed >700 and >500-fold up-regulation, respectively. The other JAZ genes had relative transcript abundance from >20- to 100-fold up-regulation at 48 h. Over the rest of the time course, they showed variable expression. Notably, all the genes showed the highest up-regulation at 48 h (Figure 9A). In response to SA treatment, most of the JAZ genes were induced moderately, with the highest up-regulation shown for BrTIFY11d (>5-fold). Conversely, BrTIFY9a, 10a, 10d, and 10e were down-regulated completely after exogenous SA treatment. The rest of the genes showed up-regulation at least at one time point in response to SA, where the maximum fold change with respect to the control was >3 (Figure 9B). After exogenous ABA treatment in Chiifu, JAZ genes were expressed differentially. BrTIFY3a, 3b, 6a, 6b, 7b, 9a, 9b, 9c, 10a, 10c, 10d, 10e, 10f, 11a, 11c, and 11e showed drastic up-regulation and the highest induction for each gene was recorded at 4 h. Their expression during other time points (1, 12, 24, and 48 h) were very low. Notably, BrTIFY6a, 6b, 7b, 10a, 10d, 10e, and 11b were found to be down-regulated below the reference (0 h) level. The expression of BrTIFY9b was consistent and showed the highest up-regulation (>35-fold) upon ABA treatment among all of the JAZ genes. Two other copies of this gene (9a and 9c) also showed up-regulation after ABA treatment. Besides, only two genes BrTIFY7c and 10b were down-regulated completely after ABA treatment (Figure 9C).
Figure 9. Showing RT-qPCR expression analyses of 20 B. rapa JAZ genes in response to (A) JA, (B) SA, and (C) ABA treatment (0–48 h). The error bars represent the standard error of the means of three independent replicates. Values denoted by the same letter did not differ significantly at P < 0.05 according to Duncan's multiple range tests.
Expression Patterns of JAZ Genes upon Fusarium Infection in B. rapa
TIFY genes are involved in responses to different biotic elicitors in Arabidopsis and grape (Demianski et al., 2012; Zhang et al., 2012). Zhang et al. (2012) reported that TIFY genes in grape showed significant induction in response to fungal infection. In this study, we analyzed the expression of JAZ genes in response to F. oxysporum f.sp. conglutinans, which specifically attacks Brassica species and causes wilt and root rot diseases. After Fusarium infection of B. rapa Chiifu, we observed the responses of JAZ genes up to 12 days (Figure 10, Supplementary Figure 3). Among 21 JAZ genes, 16 showed several fold up-regulations during Fusarium infection compared to mock treatments. In general, most of the genes were strongly induced at early time points (3 and 6 h) and were down-regulated especially during the period from 1 to 9 d, followed by up-regulation at 12 d. BrTIFY9b showed the highest up-regulation, by >60-fold at 3 h compared to the control, which was 6-fold greater induction than that exhibited in the mock treatment (10-fold). The expression then decreased from 6 h to 9 d, followed by up-regulation again at 12 d by 3-fold compared to mock treatment. BrTIFY9c is a close homolog of BrTIFY9b and other members of the JAZ subfamily, viz. BrTIFY3a, 3b, 6a, 7b, 10c, 10d, and 10e, also showed similar patterns of expression after Fusarium infection. Other JAZ members such as BrTIFY10a, 10b, 11b, 11c, 11d, and 11e showed differential up-regulation throughout the infection period. Most of these genes had higher up-regulation by several fold during infection compared to mock treatment. Notably, this group of JAZ genes showed relatively higher expression at the later stage of infection (12 d).
Figure 10. Showing RT-qPCR expression analyses of 16 JAZ genes in B. rapa after infection with Fusarium oxysporum f.sp. conglutinans (0–12 d). The error bars represent the standard error of the means of three independent replicates. Values denoted by the same letter did not differ significantly at P < 0.05 according to Duncan's multiple range tests.
In a nutshell, from the stress treatment experiments, we identified several JAZ genes showing up-regulation in response to different abiotic, biotic and hormone treatments that might be good candidates for developing stress resistance in B. rapa. In particular, BrTIFY9a, 9b, 9c and several genes from the BrIFY10 and BrIFY11 subgroups were strongly induced by most of the stresses.
Analysis of TIFY Protein Interaction Networks in B. rapa
We conducted a network interaction analysis for all 36 B. rapa TIFY proteins to predict their interactions and associations based on Arabidopsis proteins using STRING software (Figure 11). Arabidopsis CORONATINE INSENSITIVE 1 (COI1), NOVEL INTERACTOR OF JAZ (NINJA) and bHLH transcription factor family members MYC2 and MYC3, predicted to function in the processes of JA signaling in which they were supposed to interact with all of the B. rapa TIFY proteins. COI1, NINJA, and MYC2/MYC3 proteins have been reported to take part in JA signaling and different JA-dependent functions viz. root growth, plant fertility, biotic and abiotic stress resistance (Xie et al., 1998; Chini et al., 2007; Thines et al., 2007; Browse, 2009). We identified constitutive expression of B. rapa JAZ proteins in different tissues and in different growth stages during flower development. Most of these JAZ proteins were also found to be induced strongly in response to JA, ABA, environmental stresses including cold, salt, and drought and necrotrophic pathogen Fusarium. We speculate that these proteins might interact with COI1, NINJA, and MYC2/MYC3 in a JA-dependent manner for their functions, especially against different stresses and during flower development in B. rapa.
Figure 11. Interaction network of 36 TIFY proteins identified in B. rapa with their predicted functional partners based on related proteins in Arabidopsis. Stronger associations are represented by thicker lines.
Members of the TIFY family have been demonstrated to be putative TFs with various stress responsiveness (Ye et al., 2009; Demianski et al., 2012; Zhang et al., 2012). Particularly, the JAZ subfamily genes are the best characterized as responsive to stresses to date (Li et al., 2015). However, there is no information regarding the expression and functions of this gene family in economically important crop species B. rapa, for which production is affected by various abiotic and biotic stresses. Hence, we conducted this genome-wide identification and characterization of B. rapa TIFY genes, including their tissue-specific expression in different organs and in response to different abiotic, biotic, and hormone treatments.
Specifically, we studied the expression patterns of JAZ subfamily genes in different tissues and found their widespread expression in flower buds (Figure 6). Subsequent analysis of these genes in different floral tissues and different flower growth stages also revealed their likely involvement in the development of B. rapa flowers. Members of the JAZ subfamily have been reported to play roles in flower induction (Kim et al., 2011), stamen development (Song et al., 2011), and defense responses against biotic (Shoji et al., 2008; Sun et al., 2011; Demianski et al., 2012; Oh et al., 2012) and abiotic (Ye et al., 2009; Seo et al., 2011; Ismail et al., 2012; Zhu et al., 2012) stresses. Specifically, these genes have been associated with JA responses, where jasmonates stimulate the Skp1-Cullin-F-box (SCF) E3 ubiquitin ligase complex SCFCOI1 -mediated degradation of JAZ proteins to derepress transcription factors like MYC2 and MYC3 (Yan et al., 2007; Melotto et al., 2008; Bai et al., 2011). In this mechanism, up-regulation of certain JAZ genes is considered to be a method of fine-tuning JA responses (Chung et al., 2009).
Rice JAZ genes are induced by abiotic stresses like low temperature, drought and salinity, and OsTIFY11a was found to improve stress tolerance (Ye et al., 2009). In grape, 11 TIFY genes were found responsive to osmotic and cold stresses (Zhang et al., 2012). In our study, we found JAZ subfamily genes that showed differential up-regulation against cold, salt and drought treatments in B. rapa. From low temperature-treated whole-genome microarray data and subsequent validation via qPCR, we identified some good candidate B. rapa JAZ genes (BrTIFY3a, 6a, 6b, 7b, 9b, 9c, 10c, 10f, 11a, and 11c) to be involved in response to low temperatures. This group of genes was also implicated in salinity and drought stress response, during which they showed stronger induction than under cold stress. More specifically, these genes had the highest up-regulation in response to salt stress among three abiotic stresses (in case of BrTIFY9a, 9b, and 9c; Figures 7, 8).
In response to hormone treatments, B. rapa JAZ genes were highly regulated by JA and ABA, but not by SA. Several studies have been demonstrated that JA treatment and environmental cues rapidly trigger JAZ gene expression, which might be responsible for moderating the response to JA (Chini et al., 2007; Thines et al., 2007; Katsir et al., 2008a). Functional evidence shows that, in response to different phytohormones, JAZ genes seem to have wider regulatory roles in various aspects of plant development than those performed by other subfamily genes of TIFY (Zhang et al., 2012). In this study, we found that all B. rapa JAZ genes were strongly induced and showed high up-regulation (by >15-fold to >800-fold compared to the reference), especially at the final time point (48 h) of JA treatment (Figure 9A). However, against SA, 15 genes showed up-regulation at lower levels (maximum >5 fold) and the rest of the genes were down-regulated and became inactive (Figure 9B). ABA plays key roles in helping the plant to adjust to adverse environmental conditions like drought, cold and salinity (Davies and Jones, 1991). For instance, 20 rice TIFY genes exhibited strong induction in response to ABA and most of these genes were also induced by abiotic stresses (Ye et al., 2009). In our analysis, most of the abiotic stress-responsive JAZ genes were also induced by exogenous ABA treatment. Only three genes (BrTIFY10a, 11a, and 11b) were seemingly not regulated by ABA (Figure 9C). It is possible that both ABA-dependent and ABA-independent signaling pathways are responsible for regulating the expression of B. rapa JAZ genes.
JAZ genes play roles against different biotic agents like insect herbivory and pathogen infection (Shoji et al., 2008; Sun et al., 2011; Demianski et al., 2012; Oh et al., 2012) and many of these resistance mechanisms are mediated by a variety of signaling molecules that include JA, SA, and ET (ethylene) (Bari and Jones, 2009; Pieterse et al., 2009). Jasmonates might play roles in plant resistance against necrotrophic pathogens like Alternaria brassicicola and Botrytis cinerea (Glazebrook, 2005). In this study, the majority of the JAZ genes (16 out of 21 genes) showed responsiveness to the necrotrophic pathogen F. oxysporum in the form of differential up-regulation of their transcript levels (Figure 10). Notably, several genes exhibited up-regulation during the early period (3 and 6 h) of infection (BrTIFY3a, 3b, 6a, 7b, 9b, 9c, 10c, 10d, and 10e). Another set of JAZ genes (BrTIFY10a, 10b, 11b, 11c, 11d, and 11e) showed differential up-regulation throughout the infection period (3 h–12 d), with the highest expression at 12 d. Thus, both groups of genes are considered to be responsive at differential levels to Fusarium infection in B. rapa. This phenomenon may be explained by the fact that expression of B. rapa JAZ genes could be induced by a necrotrophic pathogen like Fusarium through triggering activation of JA-dependent defense responses. In contrary, there is some previous evidence that JAZ genes are induced by biotrophic pathogens. For instance, Arabidopsis JAZ genes showed up-regulation in response to the biotrophic pathogen Pseudomonas syringae. Generally, though, the SA-dependent pathway is considered to be the regulating factor in activation of defense responses to biotrophic infection (Glazebrook, 2005). However, up-regulation of Arabidopsis JAZ following infection with P. syringae might be related to the fact that the virulence factor coronatine is an inducer of JA/ET signaling (Glazebrook, 2005; Demianski et al., 2012).
Furthermore, we may look for some intersecting points among the signaling pathways of JA, ABA, and SA. Expression of genes in response to stresses might converge or interact in some points of the signaling pathways. For example, ABA and JA have been reported to act synergistically during drought stress signaling. The overlapping signaling cascade of ABA and JA certainly needs a common player to interact. MYC2 has been identified to play roles for multiple hormone signaling pathways (Riemann et al., 2015). It has also been reported that TF MYC2 activates downstream genes like RD22/RD26 whose expression is induced by pathogens, drought, salt, ABA, JA treatment (Fujita et al., 2004, 2006; Harb et al., 2010). In our study, JAZ genes were induced by multiple treatments like JA, ABA, SA as well as by abiotic stresses including low temperature, drought, salinity, and pathogen infection. We may postulate that JAZ genes in B. rapa responded against different stresses following different hormone signaling pathways and TFs like MYC2/3/4 might be a common interacting point to regulate different downstream genes.
In conclusion, among the four TIFY subfamilies, JAZ proteins have been characterized extensively as key factors in JA signaling, whereas less is known about the other subfamilies. However, remarkable progress in the identification and characterization of TIFY family genes in different higher plants has been made using functional evidence regarding these genes in the model plant Arabidopsis. Our current comprehensive and systemic analysis of TIFY genes in the B. rapa genome identified seven TIFY, six ZML, two PPD and 21 JAZ genes. The JAZ family genes in B. rapa Chiifu were found to be induced by low temperature more highly than members of the other subfamilies. The JAZ genes had constitutive expression patterns, mainly in the reproductive parts and during flower development in B. rapa. We further identified potential involvement of most of the JAZ genes in tolerance to cold, salt and drought. B. rapa JAZ genes were highly induced by JA/ABA and moderately by SA. JAZ genes were also induced strongly in response to the necrotrophic fungal pathogen Fusarium. Taken as a whole, our study identified several B. rapa JAZ genes (BrTIFY3a, 3b, 6a, 9a, 9b, and 9c) that showed corresponsive expression against all the stresses and phytohormone treatments. This study thus supplements our understanding of this gene family in plants and provides the basis for future functional studies on TIFY genes in B. rapa.
IN, JP planned the work and provided suggestions to organize data to first author and also provided editorial comments. GS designed and conducted all experiments under the supervision of IN, JP. MK assisted GS with experimentation and data collection. All authors read and agreed the manuscript contents.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This research was supported by Golden Seed Project (Center for Horticultural Seed Development, No. 213003-04-3-SB110, 213003-04-3-SB420) Ministry of Agriculture, Food and Rural Affairs (MAFRA), Ministry of Oceans and Fisheries (MOF), Rural Development Administration (RDA) and Korea Forest Service (KFS).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00936
Ahmed, N. U., Park, J. I., Jung, H. J., Kang, K. K., Lim, Y. P., Hur, Y., et al. (2013). Molecular characterization of thaumatin family genes related to stresses in Brassica rapa. Sci. Hortic. 152, 26–34. doi: 10.1016/j.scienta.2013.01.007
Aparicio-Fabre, R., Guillén, G., Loredo, M., Arellano, J., Valdés-López, O., Ramírez, M., et al. (2013). Common bean (Phaseolus vulgaris L.) PvTIFY orchestrates global changes in transcript profile response to jasmonate and phosphorus deficiency. BMC plant Biol. 13:26. doi: 10.1186/1471-2229-13-26
Bai, Y., Meng, Y., Huang, D., Qi, Y., and Chen, M. (2011). Origin and evolutionary analysis of the plant-specific TIFY transcription factor family. Genomics 98, 128–136. doi: 10.1016/j.ygeno.2011.05.002
Bowers, J. E., Abbey, C., Anderson, S., Chang, C., Draye, X., Hoppe, A. H., et al. (2003). A high-density genetic recombination map of sequence-tagged sites for Sorghum, as a framework for comparative structural and evolutionary genomics of tropical grains and grasses. Genetics 165, 367–386. Available online at: http://www.genetics.org/content/165/1/367.full
Chini, A., Fonseca, S., Fernandez, G., Adie, B., Chico, J. M., Lorenzo, O., et al. (2007). The JAZ family of repressors is the missing link in jasmonate signalling. Nature 448, 666–671. doi: 10.1038/nature06006
Chung, H. S., Koo, A. J., Gao, X., Jayanty, S., Thines, B., Jones, A. D., et al. (2008). Regulation and function of Arabidopsis JASMONATE ZIM-domain genes in response to wounding and herbivory. Plant Physiol. 146, 952–964. doi: 10.1104/pp.107.115691
Demianski, A. J., Chung, K. M., and Kunkel, B. N. (2012). Analysis of Arabidopsis JAZ gene expression during Pseudomonas syringae pathogenesis. Mol. Plant Pathol. 13, 46–57. doi: 10.1111/j.1364-3703.2011.00727.x
Dong, X., Im, S. B., Lim, Y. P., Nou, I. S., and Hur, Y. (2014). Comparative transcriptome profiling of freezing stress responsiveness in two contrasting Chinese cabbage genotypes, Chiifu and Kenshin. Genes Genomics 36, 215–227. doi: 10.1007/s13258-013-0160-y
Fujita, M., Fujita, Y., Noutoshi, Y., Takahashi, F., Narusaka, Y., Yamaguchi-Shinozaki, K., et al. (2006). Crosstalk between abiotic and biotic stress responses: a current view from the points of convergence in the stress signaling networks. Curr. Opin. Plant Biol. 9, 436–442. doi: 10.1016/j.pbi.2006.05.014
Fujita, M., Fujita, Y., Maruyama, K., Seki, M., Hiratsu, K., Ohme-Takagi, M., et al. (2004). A dehydration induced NAC protein, RD26, is involved in a novel ABA- dependent stress-signaling pathway. Plant J. 39, 863–876. doi: 10.1111/j.1365-313X.2004.02171.x
Harb, A., Krishnan, A., Ambavaram, M. M., and Pereira, A. (2010). Molecular and physiological analysis of drought stress in Arabidopsis reveals early responses leading to acclimation in plant growth. Plant Physiol. 154, 1254–1271. doi: 10.1104/pp.110.161752
Hori, Y., Kurotani, K. I., Toda, Y., Hattori, T., and Takeda, S. (2014). Overexpression of the JAZ factors with mutated Jas domains causes pleiotropic defects in Rice spikelet development. Plant Signal. Behav. 9:e970414. doi: 10.4161/15592316.2014.970414
Huang, D., Wu, W., Abrams, S. R., and Cutler, A. J. (2008). The relationship of drought-related gene expression in Arabidopsis thaliana to hormonal and environmental factors. J. Exp. Bot. 59, 2991–3007. doi: 10.1093/jxb/ern155
Jung, H. J., Dong, X., Park, J. I., Thamilarasan, S. K., Lee, S. S., Kim, Y. K., et al. (2014). Genome-wide transcriptome analysis of two contrasting Brassica rapa doubled haploid lines under cold-stresses using Br135K oligomeric chip. PLoS ONE 9:e106069. doi: 10.1371/journal.pone.0106069
Katsir, L., Schilmiller, A. L., Staswick, P. E., He, S. Y., and Howe, G. A. (2008b). COI1 is a critical component of a receptor for jasmonate and the bacterial virulence factor coronatine. Proc. Natl. Acad. Sci. U.S.A. 105, 7100–7105. doi: 10.1073/pnas.0802332105
Kim, K. C., Han, J. A., Lee, J., Maeng, J., and Hur, Y. (2011). Gene encoding PnFL-2 with TIFY and CCT motifs may control floral induction in Pharbitis nil. Genes Genomics 33, 229–236. doi: 10.1007/s13258-010-0174-7
Kim, Y. B., Li, X., Kim, S. J., Kim, H. H., Lee, J., Kim, H., et al. (2013). MYB transcription factors regulate glucosinolate biosynthesis in different organs of Chinese cabbage (Brassica rapa ssp. pekinensis). Molecules 18, 8682–8695. doi: 10.3390/molecules18078682
Koch, M. A., Haubold, B., and Mitchell-Olds, T. (2000). Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). Mol. Biol. Evol. 17, 1483–1498. Available online at: http://mbe.oxfordjournals.org/content/17/10/1483.short
Kong, X., Lv, W., Jiang, S., Zhang, D., Cai, G., Pan, J., et al. (2013). Genome-wide identification and expression analysis of calcium-dependent protein kinase in maize. BMC genomics 14:433. doi: 10.1186/1471-2164-14-433
Li, X., Yin, X., Wang, H., Li, J., Guo, C., Gao, H., et al. (2015). Genome-wide identification and analysis of the apple (Malus × domestica Borkh.) TIFY gene family. Tree Genet. Genomes 11, 1–13. doi: 10.1007/s11295-014-0808-z
Melotto, M., Mecey, C., Niu, Y., Chung, H. S., Katsir, L., Yao, J., et al. (2008). 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. 55, 979–988. doi: 10.1111/j.1365-313X.2008.03566.x
Nishii, A., Takemura, M., Fujita, H., Shikata, M., Yokota, A., and Kohchi, T. (2000). Characterization of a novel gene encoding a putative single zinc-finger protein, ZIM, expressed during the reproductive phase in Arabidopsis thaliana. Biosci. Biotechnol. Biochem. 64, 1402–1409. doi: 10.1271/bbb.64.1402
Oh, Y., Baldwin, I. T., and Gális, I. (2012). NaJAZh regulates a subset of defense responses against herbivores and spontaneous leaf necrosis in Nicotiana attenuata plants. Plant Physiol. 159, 769–788. doi: 10.1104/pp.112.193771
Riemann, M., Dhakarey, R., Hazman, M., Miro, B., Kohli, A., and Nick, P. (2015). Exploring jasmonates in the hormonal network of drought and salinity responses. Front. Plant Sci. 6:1077. doi: 10.3389/fpls.2015.01077
Saha, G., Park, J. I., Jung, H. J., Ahmed, N. U., Kayum, A., Chung, M. Y., et al. (2015). Genome-wide identification and characterization of MADS-box family genes related to organ development and stress resistance in Brassica rapa. BMC Genomics 16:178. doi: 10.1186/s12864-015-1349-z
Seo, J. S., Joo, J., Kim, M. J., Kim, Y. K., Nahm, B. H., Song, S. I., et al. (2011). OsbHLH148, a basic helix-loop-helix protein, interacts with OsJAZ proteins in a jasmonate signaling pathway leading to drought tolerance in rice. Plant J. 65, 907–921. doi: 10.1111/j.1365-313X.2010.04477.x
Shikata, M., Matsuda, Y., Ando, K., Nishii, A., Takemura, M., Yokota, A., et al. (2004). Characterization of Arabidopsis ZIM, a member of a novel plant-specific GATA factor gene family. J. Exp. Bot. 55, 631–639. doi: 10.1093/jxb/erh078
Singh, A. P., Pandey, B. K., Deveshwar, P., Narnoliya, L., Parida, S. K., and Giri, J. (2015). JAZ repressors: potential involvement in nutrients deficiency response in rice and chickpea. Front. Plant. Sci. 6:975. doi: 10.3389/fpls.2015.00975
Song, S., Qi, T., Huang, H., Ren, Q., Wu, D., Chang, C., et al. (2011). The jasmonate-ZIM domain proteins interact with the R2R3-MYB transcription factors MYB21 and MYB24 to affect jasmonate-regulated stamen development in Arabidopsis. Plant Cell 23, 1000–1013. doi: 10.1105/tpc.111.083089
Song, X. M., Huang, Z. N., Duan, W. K., Ren, J., Liu, T. K., Li, Y., et al. (2014). Genome-wide analysis of the bHLH transcription factor family in Chinese cabbage (Brassica rapa ssp. pekinensis). Mol. Genet. Genomics 289, 77–91. doi: 10.1007/s00438-013-0791-3
Talalay, P., and Fahey, J. W. (2001). Phytochemicals from cruciferous plants protect against cancer by modulating carcinogen metabolism. J. Nutr. 131, 3027S–3033S. Available online at: http://jn.nutrition.org/content/131/11/3027S.short
Thines, B., Katsir, L., Melotto, M., Niu, Y., Mandaokar, A., Liu, G., et al. (2007). JAZ repressor proteins are targets of the SCFCOI1 complex during jasmonate signalling. Nature 448, 661–665. doi: 10.1038/nature05960
Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F., and Higgins, D. G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25, 4876–4882. doi: 10.1093/nar/25.24.4876
Town, C. D., Cheung, F., Maiti, R., Crabtree, J., Haas, B. J., Wortman, J. R., et al. (2006). Comparative genomics of Brassica oleracea and Arabidopsis thaliana reveal gene loss, fragmentation, and dispersal after polyploidy. Plant Cell 18, 1348–1359. doi: 10.1105/tpc.106.041665
Withers, J., Yao, J., Mecey, C., Howe, G. A., Melotto, M., and He, S. Y. (2012). Transcription factor-dependent nuclear localization of a transcriptional repressor in jasmonate hormone signaling. Proc. Natl. Acad. Sci. U.S.A. 109, 20148–20153. doi: 10.1073/pnas.1210054109
Woodhouse, M. R., Cheng, F., Pires, J. C., Lisch, D., Freeling, M., and Wang, X. (2014). Origin, inheritance, and gene regulatory consequences of genome dominance in polyploids. Proc. Natl. Acad. Sci. U. S. A. 111, 5283–5288. doi: 10.1073/pnas.1402475111
Xie, D. X., Feys, B. F., James, S., Nieto-Rostro, M., and Turner, J. G. (1998). COI1: an Arabidopsis gene required for jasmonate-regulated defense and fertility. Science 280, 1091–1094. doi: 10.1126/science.280.5366.1091
Yan, Y., Stolz, S., Chételat, A., Reymond, P., Pagni, M., Dubugnon, L., et al. (2007). A downstream mediator in the growth repression limb of the jasmonate pathway. Plant Cell 19, 2470–2483. doi: 10.1105/tpc.107.050708
Yanhui, C., Xiaoyuan, Y., Kun, H., Meihua, L., Jigang, L., Zhaofeng, G., et al. (2006). The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol. Biol. 60, 107–124. doi: 10.1007/s11103-005-2910-y
Ye, H., Du, H., Tang, N., Li, X., and Xiong, L. (2009). Identification and expression profiling analysis of TIFY family genes involved in stress and phytohormone responses in rice. Plant Mol. Biol. 71, 291–305. doi: 10.1007/s11103-009-9524-8
Zhang, Y., Gao, M., Singer, S. D., Fei, Z., Wang, H., and Wang, X. (2012). Genome-wide identification and analysis of the TIFY gene family in grape. PLoS ONE 2012:e44465. doi: 10.1371/journal.pone.0044465
Zhu, D., Bai, X., Luo, X., Chen, Q., Cai, H., Ji, W., et al. (2013). Identification of wild soybean (Glycine soja) TIFY family genes and their expression profiling analysis under bicarbonate stress. Plant Cell Rep. 32, 263–272. doi: 10.1007/s00299-012-1360-7
Zhu, D., Cai, H., Luo, X., Bai, X., Deyholos, M. K., Chen, Q., et al. (2012). Over-expression of a novel JAZ family gene from Glycine soja, increases salt and alkali stress tolerance. Biochem. Biophys. Res. Commun. 426, 273–279. doi: 10.1016/j.bbrc.2012.08.086
Keywords: TIFY, JAZ, abiotic and biotic stress, hormone, gene expression, Brassica rapa, Chiifu, Kenshin
Citation: Saha G, Park J-I, Kayum MA and Nou I-S (2016) A Genome-Wide Analysis Reveals Stress and Hormone Responsive Patterns of TIFY Family Genes in Brassica rapa. Front. Plant Sci. 7:936. doi: 10.3389/fpls.2016.00936
Received: 30 January 2016; Accepted: 13 June 2016;
Published: 28 June 2016.
Edited by:Keqiang Wu, National Taiwan University, Taiwan
Reviewed by:Junxian He, The Chinese University of Hong Kong, Hong Kong
Jitender Giri, National Institute of Plant Genome Research, India
Copyright © 2016 Saha, Park, Kayum and Nou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ill-Sup Nou, email@example.com