Next Article in Journal
Components of the Complement Cascade Differ in Polycystic Ovary Syndrome
Previous Article in Journal
Transcriptome Analysis Reveals Genes and Pathways Associated with Salt Tolerance during Seed Germination in Suaeda liaotungensis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 20
10.3390/ijms232012218
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genome-Wide Identification of C2H2 ZFPs and Functional Analysis of BRZAT12 under Low-Temperature Stress in Winter Rapeseed (Brassica rapa)

by
Li Ma
1,
Jia Xu
1,2,
Xiaolei Tao
1,2,
Junyan Wu
1,2,
Wangtian Wang
1,
Yuanyuan Pu
1,2,
Gang Yang
1,2,
Yan Fang
1,
Lijun Liu
1,
Xuecai Li
1,2,* and
Wancang Sun
1,2,*
1
State Key Laboratory of Aridland Crop Science, Gansu Agricultural University, Lanzhou 730070, China
2
College of Agronomy, Gansu Agricultural University, Lanzhou 730070, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(20), 12218; https://doi.org/10.3390/ijms232012218
Submission received: 19 September 2022 / Revised: 7 October 2022 / Accepted: 11 October 2022 / Published: 13 October 2022
(This article belongs to the Topic Temperature Stress and Responses in Plants)
Browse Figures
Review Reports Versions Notes

Abstract

:
Zinc-finger protein (ZFP) transcription factors are among the largest families of transcription factors in plants. They participate in various biological processes such as apoptosis, autophagy, and stemness maintenance and play important roles in regulating plant growth and development and the response to stress. To elucidate the functions of ZFP genes in the low-temperature response of winter (Brassica rapa L.) B. rapa, this study identified 141 members of the C2H2 ZFP gene family from B. rapa, which are heterogeneously distributed on 10 chromosomes and have multiple cis-acting elements related to hormone regulation and abiotic stress of adversity. Most of the genes in this family contain only one CDS, and genes distributed in the same evolutionary branch share mostly the same motifs and are highly conserved in the evolution of cruciferous species. The genes were significantly upregulated in the roots and growth cones of ‘Longyou-7’, indicating that they play a role in the stress-response process of winter B. rapa. The expression level of the Bra002528 gene was higher in the strongly cold-resistant varieties than in the weakly cold-resistant varieties after low-temperature stress. The survival rate and BrZAT12 gene expression of trans-BrZAT12 Arabidopsis thaliana (Arabidopsis) were significantly higher than those of the wild-type plants at low temperature, and the enzyme activities in vivo were higher than those of the wild-type plants, indicating that the BrZAT12 gene could improve the cold resistance of winter B. rapa. BrZAT12 expression and superoxide dismutase and ascorbate peroxidase enzyme activities were upregulated in winter B. rapa after exogenous ABA treatment. BrZAT12 expression and enzyme activities decreased after the PD98059 treatment, and BrZAT12 expression and enzyme activities were higher than in the PD98059 treatment but lower than in the control after both treatments together. It is speculated that BrZAT12 plays a role in the ABA signaling process in which MAPKK is involved. This study provides a theoretical basis for the resolution of cold-resistance mechanisms in strong winter B. rapa.

1. Introduction

Winter rapeseed (B. rapa) is one of the most important oil crops in the arid and cold regions of northern China. At present, the northward migration of winter rapeseed has been successfully carried out in the north, and the cultivation area of winter rapeseed has gradually changed because of its strong cold resistance and its ability to winter smoothly under the severe cold environment in the north [1,2]. The fact that winter rapeseed can successfully winter under low temperatures and shows good cold resistance indicates that it is rich in cold-resistance genes and is thus an important material for studying cold resistance.
During plant growth and development, abiotic stresses such as low temperature, high salinity, and drought affect plant growth and crop yield [3,4]. Plants have evolved protective mechanisms to cope with stresses. These involve morphological, physiological, and molecular adaptations over a long evolutionary period [5,6]. Transcription factors, also known as trans-acting factors, are an important regulator of the environment in plant cells. They can bind with promoters to activate gene transcription, participate in plant growth and development regulation, participate in signal transduction, and regulate gene expression under stress [7,8]. Therefore, transcription factors related to plant-stress responses have been the focus of much research. The transcription factors that have been identified in plants that play a role when plants encounter adverse environments include NAC, ZFP, bHLH, bZIP, and WRKY [9,10].
ZFPs have been found to exist in a variety of plants and have various biological functions, such as regulating gene expression, DNA recognition, and participating in plant growth and development. A large number of ZFPs have been discovered, among the most popular ones are C2H2, which mostly relate to the stress response and plant growth and development regulation [11,12]. ZFPs are one of the most widely studied transcription factor families. According to their structure and characteristics, ZFPs can be divided into several subfamilies, including C2H2, CCCH, C3HC4, C2HC, C2HC5, C4, C4HC3, C6, and C8, where C stands for Cys (cysteine) and H for His (histidine) [3,11]. Different types of ZFPs have different spatial configurations, resulting in different functions of the proteins. Among them, the C2H2 ZFPs have the most eukaryotic members. Each ZFP contains two cysteine (Cys2) and two histidine (His2) residues, which bind stably to Zn2+ and specifically to a domain in the gene promoter regulated by the C2H2 ZFPs [11,13,14].
According to reports, members of the C2H2 ZFP gene family have been identified in a variety of plants, such as Arabidopsis, potato, sorghum, tomato, and poplar, and it has been indicated that C2H2 ZFPs in plants play an important regulatory role in abiotic stress responses [15,16,17,18,19]. The constitutive expression of ZAT12 in Arabidopsis caused a small but reproducible increase in freezing tolerance, suggesting a role for the ZAT12 regulator in cold acclimation. In addition, ZAT12 downregulated CBF gene expression, suggesting a role for ZAT12 in a negative regulatory loop that inhibits the expression of the CBF cold-response pathway [20]. It was shown that the BcZAT12-transformed tomato lines ZT1, ZT2, and ZT6 were better adapted to drought stress by accumulating the osmotic protein proline and increasing the antioxidant response triggered by the ZAT12 gene. Thus, the ZAT12-gene-transformed tomato cv. H-86 lines will prove useful for producing high-yield tomato crops in areas subjected to severe drought stress [21]. The overexpression of soybean GmZF1 significantly improved the proline and soluble-sugar content and reduced the malondialdehyde (MDA) content in transgenic lines exposed to cold stress, suggesting that transgenic Arabidopsis carrying the GmZF1 gene has an adaptive mechanism to cold stress. The overexpression of GmZF1 also increased the expression of the cold-regulated cor6.6 gene, probably by recognizing the protein–DNA binding site, indicating that GmZF1 from soybean can improve cold-stress tolerance in Arabidopsis by regulating the expression of cold-regulated genes in transgenic Arabidopsis [22]. Based on stress-induction analysis, mutation analysis, complementation analysis, and ectopic-expression analysis, many C2H2 ZFPs involved in abiotic stress signaling pathways were identified [11]. Phytohormones are responsible for resistance to abiotic stresses and are involved in the process of the response to various stresses through C2H2 ZFPs, especially ABA [23,24].
In this study, we analyzed the gene chromosome distribution, gene evolution, motifs, cis-acting elements, and covariance of the C2H2 ZFP gene family in winter rapeseed (B. rapa) based on genomic and transcriptomic data and investigated the expression pattern of C2H2 ZFP genes under abiotic-adversity stress. The BrZAT12 gene was identified in response to cold stress, and the gene was functionally validated. This study provides a theoretical basis for the analysis of the mechanism by which the C2H2 ZFP regulates cold resistance in winter rapeseed as well as for the cultivation of strong cold-resistant quality germplasm.

2. Results

2.1. Chromosome Distribution of the C2H2 ZFP Gene Family

The C2H2 ZFP gene family in B. rapa was identified by TBtools. A total of 141 members of the C2H2 ZFP gene family were identified. According to the visual analysis results of the chromosome distribution obtained by TBtools software (Figure 1, Supplementary Table S2), it was found that the members of the gene family were unevenly distributed on 10 chromosomes. Among them, most were distributed on chromosome A03, containing 24 members of the gene family, followed by A02 (18), A09 (18), A06 (17), A07 (15), A01 (12), A10 (12), A04 (10), and A05 (9), and the least were distributed on chromosome A08, including only six gene-family members.

2.2. Phylogenetic Analysis of the C2H2 ZFP Gene Family

Furthermore, the phylogeny of the C2H2 family members of B. rapa was studied. The phylogenetic relationships showed that the C2H2 ZFP family of B. rapa, B. napus, and Arabidopsis could be divided into 16 subfamilies (I–XVI) (Figure 2). The C2H2 ZFP family of B. rapa could be divided into 12 subfamilies (I–XII) (Figure 3A). Among them, subfamily IV contained, at most, 78 genes and the subfamily XI contained at least four genes. The genetic distances of subfamilies I, II, and III and other subfamilies were distant, while those of subfamilies XIV, XV, and XVI were relatively close. The C2H2 ZFPs of B. rapa, B. napus, and Arabidopsis were unevenly distributed in each subfamily. The phylogenetic homology of C2H2 transcription factors in B. rapa, B. napus, and Arabidopsis indicates that they may play the same or similar roles in some biological functions.

2.3. Gene Structure of the C2H2 ZFP Gene Family

Using MEME for conservative motif analysis, it was found that most of the genes distributed in the same evolutionary branch had the same motif (Figure 3B, Supplementary Table S3). Gene structure analysis, as shown in Figure 3C, and Supplementary Table S2, showed that most of the genes in this family have only one CDS. A number of stress-related response elements were identified, such as low temperature response elements, anaerobic response elements, and drought response elements. In addition, there were hormone-related response elements, such as IAA, SA, ABA, GA, and MeJA, indicating that the gene family has a variety of cis-acting elements and diverse functions (Figure 3D). The C2H2 family genes in B. rapa have a variety of functional elements related to hormonal and abiotic stress, suggesting that C2H2 family genes may regulate plant growth and development and a variety of stress responses through different hormones.

2.4. Genomic Collinearity Analysis of the C2H2 ZFP Gene Family

Intraspecific collinearity analysis of the C2H2 ZFP genes in B. rapa showed that 126 pairs of genes with collinear relationships, such as Bra000245 and Bra016900, Bra000245, Bra004660, Bra004660 and Bra003317, and Bra004660 and Bra007378, were distributed on 10 chromosomes (Figure 4A). This suggests that the C2H2 ZFP gene is highly conserved and under stabilizing selection, and may be produced by gene replication. It was found that there was a collinear relationship among the C2H2 ZFP genes in B. rapa, Arabidopsis, and B. napus, indicating that the gene family was highly conserved in the process of evolution (Figure 4B). The Ka/Ks of B. rapa and B. napus ranged from 0.044 to 0.755, with an average of 0.331. The Ka/Ks of B. rapa and Arabidopsis ranged from 0.072 to 0.54, with an average of 0.262. It is suggested that the C2H2 ZFP gene is relatively conserved in cruciferous species (Figure 4C, Supplementary Table S4).

2.5. Analysis of the Gene Expression Characteristics of Members of the Gene Family

The expression of 53 and 52 genes of the cold-resistant ‘Longyou-7’ was upregulated at 3 h and 24 h, respectively (Figure 5, Supplementary Table S5). Among them, 36 genes were upregulated at 3 h and 24 h under low-temperature stress, and 22 genes, including Bra000245, Bra001752, and Bra002528, were higher at 24 h than at 3 h. Additionally, the expression levels of 11 genes such as Bra000245, Bra002528, and Bra004312 in the highly cold-resistant ‘Longyou-7’ were higher than those in the weaky cold-resistant ‘Lenox’. In addition, there were 25 downregulated genes in ‘Longyou-7’ at 3 h and 24 h, among which 10 genes, such as Bra034396, Bra035609, and Bra024706, had lower expression levels than in ‘Lenox’. After 6 h of freezing stress, there were 16 upregulated genes, 62 downregulated genes, and 63 non-expressed genes in ‘Longyou-7’ and ‘Tianyou-2’. There were 21 upregulated genes, 51 downregulated genes, and 69 non-expressed genes in ‘Tianyou-2’. There were 22 genes expressed differently between ‘Longyou-7’ and ‘Tianyou-2’. Among them, seven genes had higher expression levels in ‘Longyou-7’ than in ‘Tianyou-2’, and 15 genes had higher expression levels in ‘Tianyou-2’ than in ‘Longyou-7’.

2.6. Analysis of the Expression Pattern of Gene Family Members

Ten genes differentially expressed under low-temperature stress were screened from the transcriptome results of two varieties differing in cold resistance. The expression patterns of these genes under cold, drought, high temperature, and salt stresses were further analyzed. The untreated rapeseed was used as a control (CK) for four abiotic stress treatments. The expression induced by Bra002528 and Bra004312 was more obvious under cold stress, in which the expression of Bra002528 in the growth cone of Bra002528 ‘Longyou-7’ reached 17 times that of CK and 5.5 times that of the ‘Lenox’ growth cone at 24 h. The expression of Bra002528 in the roots of ‘Longyou-7’ at 24 h was 24.7 times that of CK and 2.8 times that of ‘Lenox’, and the expression of the Bra004312 gene in ‘Longyou-7’ was 3 times higher than that in ‘Lenox’ at 24 h. The expression level of Bra004312 was 4.6 times higher than that in ‘Lenox’ at 24 h in the growth cone (Figure 6A). At 24 h of high-temperature stress treatment, the expression of Bra002528, Bra004312, and Bra020284 genes in ‘Longyou-7’ leaves was 2.8, 2.5, and 3 times higher than that in ‘Lenox’ leaves. Expression in stems was 1.2, 2.5, and 3.1 times higher than that in ‘Lenox’. Expression in roots was 4.2, 4.2, and 3 times higher than in ‘Lenox’ (Figure 6B).
Under drought stress, the expression of Bra002528 and Bra006691 was more significantly induced by drought. The expression of Bra002528 and Bra006691 in the leaves of ‘Longyou-7’ at 24 h was 7.9 and 14.4 times that of CK and 3.3 and 8.2 times that of ‘Lenox’, respectively, while in the growth cone, the expression of Bra002528 and Bra006691 in ‘Longyou-7’ at 24 h was 20.4 and 24.8 times that of CK and 5.5 and 13.4 times that of ‘Lenox’, respectively. The expression levels of these genes in the root of ‘Longyou-7’ at 24 h was 25.8 and 10.5 times that of CK, respectively, and 5.6 and 5.7 times that of ‘Lenox’, respectively (Figure 6C). After 24 h of salt stress, the expression of Bra004660 in the ‘Longyou-7’ growth cone was 26.5 times higher than that in CK, and 6.6 times higher than that in ‘Lenox’. The expression of Bra020284 in the ‘Longyou-7’ growth cone was 27.1 times higher than that in CK, and 6.9 times higher than that in ‘Lenox’ (Figure 6D).
With these results we were surprised to find that Bra002528, Bra004312, Bra004312, Bra020284, Bra006691, and Bra004660 genes were expressed at higher levels in ‘Longyou-7’ than ‘Lenox’ after stress, and that the expression was higher in growth cones and roots, followed by leaves and stems. In addition, we identified Bra002528, Bra020284, and Bra004312 genes in response to cold- and high-temperature stress, Bra006691 gene in response to drought stress, and Bra004660 gene in response to salt stress. In addition, individual genes were found to respond to multiple abiotic stresses, e.g., Bra002528 gene responded to cold, high-temperature and drought stress; Bra004312 responded to cold and high-temperature stress; Bra020284 responded to high-temperature and salt stress.

2.7. Screening and Identification of Transgenic Arabidopsis with the BrZAT12 Gene

After infesting Arabidopsis, the plants were screened using Basta for three generations, and the T3 generation of Arabidopsis plants exhibiting normal growth was transplanted into pots for the next experiment (Supplementary Figure S1A–C). The PCR amplification results showed that wild-type Arabidopsis plants had no amplified bands, and 12 of the randomly selected transgenic plants from Z1 to Z13 showed a 483 bp band after agarose gel electrophoresis, while the Z3 plants failed to show a band. The expression of the BrZAT12 gene in transgenic Arabidopsis was higher than that in wild-type Arabidopsis and was significantly different from that in wild-type Arabidopsis (Supplementary Figure S1D). The gene expression of Z5 plants was the highest, which was 4.65 times higher than that of wild-type Arabidopsis. The results showed that all the six plants that over-expressed BrZAT12, and Z5 had the highest transformation efficiency (Supplementary Figure S1E).

2.8. Phenotypic and Expression Analysis of BrZAT12-Transformed Plants under Low-Temperature Stress

It was found that the leaves of the wild-type and transgenic Arabidopsis plants with the BrZAT12 gene exhibited no significant change after 3 h of treatment, and the survival rate was similar. After 6 h of treatment, the leaves of the plants began to wilt, the survival rate of wild-type Arabidopsis decreased by about 83.3%, and the survival rate of transgenic plants with the BrZAT12 gene was 91.7% (Figure 7A,B). After 12 h of treatment, the wild-type plants froze, wilted, and died after being returned to normal temperature for a week, and the survival rate was just 16.7%, while the survival rate of the transgenic Arabidopsis was 62.5%, which was significantly different from that of the wild-type. The survival rate of wild-type Arabidopsis treated for 24 h was 4.17%, and most of the transgenic plants with the BrZAT12 gene returned to normal growth at room temperature, and the survival rate reached 41.7%, which was significantly different from that of the wild-type (Figure 7C). The expression of the BrZAT12 gene in transgenic Arabidopsis plants reached the highest level after 12 h of treatment, which was 18.8 times higher than that in the wild-type Arabidopsis plants treated for 12 h and was significantly different from that in the wild-type Arabidopsis plants treated for 12 h (Figure 7D). It is suggested that the BrZAT12 gene can improve the cold resistance of plants.

2.9. Physiological-Activity Analysis of the BrZAT12 Gene in Transgenic Plants at Low Temperature

After 24 h of low-temperature treatment, the SOD enzyme activity of the transgenic Arabidopsis was 3.17 times higher than that of the control and 1.90 times higher than that of the wild-type plants treated for 24 h. The activity of SOD in the wild-type plants reached the maximum at 24 h, which was 1.54 times higher than that of the untreated plants (Figure 7E). Regarding the POD enzyme activity, the enzyme activity of the transgenic Arabidopsis plants reached the highest at 24 h of treatment, which was 2.54 times higher than that of the untreated Arabidopsis and 1.42 times that of the wild-type Arabidopsis. There was a significant difference in POD activity between the transgenic Arabidopsis and wild-type Arabidopsis at 24 h (Figure 7C). The enzyme activity of the wild-type Arabidopsis was 1.91 times higher than that of the untreated Arabidopsis for 24 h. The CAT activity of the transgenic Arabidopsis was 2.13 times higher than that of the untreated Arabidopsis and 1.38 times that of the wild-type Arabidopsis treated for 24 h. The CAT activity of the wild-type Arabidopsis increased by 1.61 times as much as that of the untreated Arabidopsis at 24 h, and the difference was significant (Figure 7G). The activities of the three enzymes in the transgenic plants were higher than those in the wild-type plants. The results showed that the reactive oxygen species (ROS) scavenging ability of the transgenic plants with the BrZAT12 gene was stronger than that of the wild-type plants.

2.10. The BrZAT12 Gene Is Involved in the Regulation of ABA and MAPK Signaling Pathways

The results showed that the ascorbate peroxidase (APX) activity and BrZAT12 gene expression of ‘Longyou-7’ and ‘Lenox’ were significantly lower than those of CK after treatment with PD98059 at room temperature. The SOD and APX activities of ‘Longyou-7’ and ‘Lenox’ were 1.46 and 1.22 times higher than those of CK, respectively, while the SOD and APX activities of ‘Lenox’ were 1.18 and 1.10 times higher than those of CK, respectively. The expression of BrZAT12 in ‘Longyou-7’ was significantly higher than that in ‘Lenox’. The SOD and APX activities of ‘Longyou-7’ under low-temperature treatment were 1.22 and 1.21 times higher than those of CK, respectively, and those of ‘Lenox’ were 1.08 and 1.10 times higher than those of CK, respectively. The gene expression and enzyme activity of BrZAT12 after co-treatment were higher than those of the PD98059 treatment, but were lower than those of the control. Compared with the ambient treatment, the expression changes in ‘Longyou-7’ under low-temperature treatment were relatively small, indicating that ‘Longyou-7’ was less affected by the inhibitor PD98059 and exogenous ABA under low-temperature stress (Figure 8). It is speculated that the BrZAT12 gene may play a role in the involvement of MAPKK in ABA signaling processes.

3. Discussion

ZFP genes are abundant in plants, and they play an important role in plant growth and development and in regulating the plant response to various biotic and abiotic stresses [25,26]. Ninety-nine C2H2 ZFP genes were identified in tomato and were classified into four groups, and 32 C2H2 genes were identified in response to multiple stresses. Pathway analysis identified five gluconeogenic and sphingolipid-related pathways and endocytosis pathways as being enriched [18]. In this study, we identified, expressed, cloned, and analyzed the C2H2 ZFP gene family in B. rapa. A total of 141 members of the C2H2 gene family were identified, which were distributed on 10 chromosomes. Most genes in the same evolutionary branch contained the same conserved motif, which was similar to that of B. rapa. Collinear relationship analysis showed that the gene was highly conserved in the process of evolution. In addition, 77 and 246 members of the ZFP family were identified in Arabidopsis and B. napus, respectively (Figure 2). In view of the fact that genome duplication occurred about 4 million years after the Brassica genome split from Arabidopsis ancestors, the doubling of these ZFPs in diploid Brassica plants reflects genome duplication [27]. About 7500 years ago, the genomes of B. oleracea and B. rapa also fused to form B. napus, which led to a doubling of the number of ZFPs [27]. It is generally believed that this increase in gene expression will cause an imbalance in gene expression and, over time, lead to a decrease in excessive gene expression.
An earlier study identified a total of 146 Q-type C2H2-ZFPs, including 37 in B. oleracea, 35 in B. rapa, and 74 in B. napus. The sequence similarity and alignment of these genes on the chromosome remained essentially unchanged in B. napus compared to those inherited in B. oleracea and B. rapa. By contrast, protein sequences differed more between homologous genes in B. oleracea and B. rapa and were organized more differently on the chromosomes. Overall, these 146 proteins are highly conserved, especially in known motifs [28]. A search of the cabbage genome identified 25 proteins with 2-Q zinc-finger structures. Furthermore, a comparative map position analysis of Arabidopsis and B. oleracea genes revealed that Arabidopsis inherited all 24 blocks, but the final karyotype was reorganized and reduced to five chromosomes. These homozygous blocks are found in modern Brassica species such as B. rapa and B. oleracea. These results further support the idea that the B. oleracea genome originated from ancient polyploids [29,30,31].
Under different stress treatments, the Bra002528, Bra004312, Bra020284, Bra006691, and Bra004660 genes were significantly expressed in the root and growth cone of ‘Longyou-7’, indicating that they play a certain role in the stress response of B. rapa. Among them, the Bra002528 gene was significantly expressed in the growth cone and root of the cold-resistant B. rapa ‘Longyou-7’ and was significantly expressed in different stress environments, indicating that the Bra002528 (BrZAT12) gene plays an important role in regulating the response of B. rapa to stress [32,33]. After the TaZNF gene was transferred into Arabidopsis, it was found that the chlorophyll content, proline content, and stomatal reduction of the transgenic Arabidopsis were higher than those of the wild-type Arabidopsis under salt stress, and the salt resistance of the transgenic plants was enhanced [34]. C2H2 ZFPs enhance cold resistance by directly regulating cold-related downstream genes in plants. Transgenic Arabidopsis and tobacco (Nicotiana tabacum) overexpressing SCOF-1 showed increased expression of cold-regulated (COR) genes and enhanced cold tolerance. The SCOF-1 transgenic plants recovered from cold stress faster than the control, and the T2 generation SCOF-1 transgenic Arabidopsis still expressed more COR genes, such as COR15a, COR47, and RD29B, than the wild type at normal growth temperatures. SCOF-1 significantly enhanced the binding activity of SGBF-1 to the ABRE sequence in vitro, thereby promoting COR gene expression and enhancing cold tolerance. The interaction of these two proteins is required for the role of ABRE in COR gene expression to enhance cold tolerance [3,35,36]. In addition, SCOF-1 enhanced cold tolerance in transgenic sweet potato and transgenic potato [37,38].
It was shown that the C2H2 ZFPs could enhance plant cold resistance by increasing the level of osmotic substances. ZFP182 significantly enhanced cold tolerance in rice overexpression lines by increasing the expression of OsP5CS and OsLEA3 and the accumulation of osmoprotectants [39]. GmZF1 enhanced cold tolerance in transgenic Arabidopsis by increasing proline and soluble sugar contents and reducing membrane lipid peroxidation under low-temperature stress [22]. In this study, after low-temperature treatment of the transgenic BrZAT12 plants, it was found that there was less wilting of, and damage to, the transgenic Arabidopsis leaves compared to the wild-type Arabidopsis, and the survival rate of the transgenic Arabidopsis was significantly higher than that of the wild-type Arabidopsis, indicating that the BrZAT12 gene did improve the cold resistance of the plants. The enzyme activities of both the wild-type and BrZAT12 transgenic plants increased after low-temperature treatment, and the increasing trend of enzyme activities in the BrZAT12 transgenic Arabidopsis plants was higher than that in the wild-type plants. The results indicated that the BrZAT12 gene was associated with cold resistance, and the transgenic Arabidopsis with the BrZAT12 gene demonstrated improved cold resistance [40].
C2H2 ZFPs play an important role in plant tolerance to low temperatures through ABA-dependent and ABA-independent pathways [13,24]. In addition to the ABA signaling pathway, C2H2 ZFPs enhance abiotic stress resistance through the MAPK signaling pathway [7,11]. Thus, C2H2 ZFPs regulate abiotic stress responses through the ABA signaling pathway and MAPK signaling pathway, constituting an interactive and complex regulatory network [20,41]. It was reported that the overexpression of ZFP245 improved tolerance to drought and cold damage in rice. In addition, the overexpression of ZFP245 promoted rice resistance to H2O2 and increased the content of osmoregulatory substances. This indicates that ZFP245 can enhance tolerance to abiotic and oxidative stresses through proline biosynthesis and scavenging of ROS. Importantly, the overexpression of ZFP245 enhanced plant sensitivity to exogenous ABA, suggesting that ZFP245 is likely involved in abiotic stress responses through an ABA-dependent signaling pathway [42,43]. The overexpression of ZFP36 increased the expression level and activity of SOD and APX, which in turn enhanced tolerance to water and oxidative stresses. Meanwhile, ZFP36 regulates the expression level of NADPH oxidase and H2O2 production in the ABA signaling pathway [44,45].
ZAT12 is involved in abiotic-stress regulation through ABA non-dependent signaling pathways and the MAPK signaling pathway as an active repressor [46,47,48]. EIN3, which is phosphorylated by MPK3/6, can directly upregulate ZAT12 expression, indicating that ZAT12 may be involved in the MAPK signaling pathway [49,50,51]. In addition, the expression of ZAT12 and ZAT7 preceded the increased expression of APX1 in response to oxidative stress. The expression of ZAT12 and ZAT7 in the loss-of-function mutant apx1 increased with increasing levels of H2O2 [52,53]. Another study showed that CBF is also regulated by ZAT12 and acts upstream of ZAT10 and AZF1/2/3 through the MAPK signaling pathway. Thus, different ZAT proteins may be highly related and involved in responses to abiotic stresses in a similarly coordinated and cooperative manner [20,54,55]. The results showed that the activities of SOD and APX and the expression of the BrZAT12 gene in plants increased after exogenous ABA treatment, and the increased trend in enzyme activity of ‘Longyou-7’ was greater than that of ‘Lenox.’ However, the inhibitor treatment of MAPKK could decrease the enzyme activity and expression. The sensitivity of plants to ABA and inhibitors decreased at low temperature. After the combination of exogenous ABA and inhibitors, the enzyme activity and expression were still lower than those in the control group. It is speculated that the BrZAT12 gene may play a role in the ABA signaling process in which MAPKK is involved.

4. Materials and Methods

4.1. Plant Growth and Treatments

Two B. rapa varieties ‘Longyou-7’ (ultra-cold resistant, AA, 2n = 20), ‘Tianyou-2’ (cold resistant, AA, 2n = 20), and ‘Lenox’ (weakly-cold resistant, AA, 2n = 20) with differential cold resistance were selected. ‘Longyou-7’ and ‘Lenox’ seeds of uniform size and healthy fullness were selected, sterilized using 5% NaClO solution, rinsed with distilled water after sterilization, and placed evenly on filter paper for germination. After the seeds germinated, they were sown in pots. Four seeds were sown in each pot, and the pots were incubated in a light incubator (14 h light, 25 °C; 10 h dark, 20 °C) until the seven-leaf stage. The untreated rapeseed was used as a control (CK) for four abiotic stress treatments: low-temperature treatment (in a 4 °C incubator), high-temperature treatment (in a 40 °C incubator), drought treatment (simulating drought stress using 18% polyethylene glycol solution), and high salt treatment (simulating salt stress using 180 mmol/L NaCl solution). Four parts, namely, the leaves, roots, stems, and growth cones, were used for subsequent experiments at 0, 1, 3, and 24 h after various treatments. Another set of experiments was conducted to study the involvement of the BrZAT12 gene in the regulation of MAPK and ABA signaling. The seeds of ‘Longyou-7’ and ‘Lenox’ were cultured in MS medium for 1 week and then transferred to four groups of different treatment media for 10 days as follows: control group, MS medium; ABA treatment group, MS medium + ABA (100 μM); MAPKK inhibitor treatment group, MS medium + PD98059 (MAPKK kinase specificity inhibitor) (50 μM); ABA + MAPKK inhibitor treatment group, and MS medium + PD98059 (50 μM) + ABA (100 μM). The above samples were snap frozen in liquid nitrogen and subsequently stored in a refrigerator at −80 °C, followed by RNA extraction and reverse transcription, with three biological replicates set up for each treatment.

4.2. Identification of Members of the C2H2 ZFP Gene Family

In this study, B. rapa genome and gene structure annotation information with high homology were downloaded from the Brassica Database (http://brassicadb.org/brad/index.php/ (accessed on 13 March 2022)). The C2H2 ZFP domain (PF13912) HMM was downloaded from the Pfam (http://pfam.xfam.org/family/PF13912 (accessed on 13 March 2022)) [56,57]. Blast was used to screen the members of the C2H2 ZFP family in B. rapa. The CDD online tool (http://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi (accessed on 13 March 2022)) was used for further screening [28,58]. TBtools software (https://github.com/CJ-Chen/TBtools (accessed on 14 March 2022)) was used to determine the position of the C2H2 ZFP gene family members on chromosomes in B. rapa, and visual analysis was carried out [59].

4.3. Analysis of Members of the C2H2 ZFP Gene Family

Conservative motif analysis was carried out on the MEME (http://meme-suite.org/tools/meme (accessed on 17 March 2022)), and the number of motifs identified was set to 10 [60]. The upstream 2000 bp sequence of the gene initiation codon was selected by TBtools and submitted to PlantCare (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/ (accessed on 17 March 2022)) for cis-acting element analysis [61]. The evolutionary tree, motifs, gene structures, and cis-acting elements of the C2H2 gene family members of B. rapa were combined and visually analyzed using TBtools software. The C2H2 ZFP sequences of Brassica napus and Arabidopsis were downloaded from the PlantTFDB database (http://planttfdb.cbi.pku.edu.cn/ (accessed on 13 March 2022)) [62,63]. A phylogenetic tree was constructed using MEGA7.0 software (neighbor-joining, NJ) with the bootstrap value set to 1000, and the tree was edited using the Evolview online website (http://www.evolgenius.info/evolview/ (accessed on 19 March 2022)) [64,65].

4.4. Collinear Analysis of Gene Family Members

The collinearity of members of the C2H2 gene family in B. rapa was analyzed using the Amazing Super Circos tool of TBtools. The collinear relationship between B. rapa, B. napus, and Arabidopsis was visually analyzed by TBtools software. The Ka and Ks replacement rates of B. rapa, B. napus, and Arabidopsis were calculated using TBtools’ Calculator tool [59]. MapChart software was used to map the chromosomal location of the C2H2 gene family [66].

4.5. RNA Isolation, Reverse Transcription, qRT-PCR, and Transcriptome Expression Analysis

RNA Extraction Kits (TIANGEN Biotech Co., Ltd., DP432, Beijing, China) were used to extract total RNA from ‘Longyou-7’ and ‘Lenox’ and remove genomic DNA contamination. The cDNA was obtained using the PrimeScript™ RT Master Mix (Takara, RR036A, Beijing, China) according to the instructions. The cDNA was checked for quality and concentration using an ultra-trace UV-visible spectrophotometer (NanoVueTM Plus, Wilmington, DE, USA) and kept on standby. Ten differentially expressed genes in the transcriptome data were selected, and specific primers were designed for qRT-PCR (ABI QuantStudio 5, Thermo Fisher Scientific, Shanghai, China) under low-temperature, high-temperature, salt, and drought stress (Supplementary Table S1) according to the instructions of the TB Green® Premix Ex Taq™ II (Takara, RR820A, Beijing, China). The reaction process was as follows: 5 min at 95 °C, followed by 40 cycles of 5 s at 95 °C and 30 s at 60 °C, followed by 65–95 °C melting curve detection. The qRT-PCR efficiency of the genes was obtained by analyzing the standard curve of the cDNA gradient dilution. The gene fragment encoding B. rapa β-actin RNA was used as an internal control to normalize the amount of template cDNA. The relative expression values of each gene were calculated using the comparative 2−ΔΔCT method [67]. The RNA-seq libraries (SRP179662) and (SRP211768) were selected for gene-expression analysis [68,69]. Three biological replicates were set up for each experimental treatment.

4.6. Preliminary Verification of BrZAT12 Function

The BrZAT12 (Bra002528) gene of B. rapa was cloned by PCR based on specific primers (BrZFP1-F: ATGGTTGCTATTTCAGAGATCAAGTCGACG; BrZFP1-R: ACAAACAGGCCTTCCAAGTTC). With reference to Gateway technology, PCR amplification was carried out using primers (BrZAT12-vec-F: AAAAAAGCAGGCTTCATGGTTGCTATTTCAGAGATCAAGTCGACG; BrZAT12-vec-R: AGAAAGCTGGGTCACAAACAGGCCTTCCAAGTTC), and the recovered products were ligated with BP Clonase enzyme (Invitrogen, Carlsbad, CA, USA) and the pDONR vector. After transforming into coliform bacteria (DH5α), the positive clones were screened by smear sequencing. The positive clones were selected to interact with LR Clonase enzyme (Invitrogen, Carlsbad, CA, USA) to connect to the overexpression vector pEarly-Gate101, following the instructions for Agrobacterium tumefaciens GV3101 (Anyu Biotechnology Co., Ltd., Shanghai, China). After the bacteria were picked out from the coated plate, colony PCR was carried out and the positive bacterial solution was selected and preserved.
Transgenic Arabidopsis was obtained by the floral-dip method. The T3 seeds were obtained by screening by spraying with 0.01% Basta (10%) solution. The T3 seeds of transgenic Arabidopsis lines and wild-type Arabidopsis were sown in pots according to the experimental requirements. After low-temperature treatment in the incubator at −4 °C for 0, 3, 6, 12, and 24 h, one part was subjected to RNA extraction and physiological index determination, and the other part resumed normal growth, and the plant phenotype and survival rate were observed and counted after 1 week. Primers (BrZAT12-pcr-F: GCGAACTGTCTGATGCTC; BrZAT12-pcr-R: TTAGGCTTCTTGTGGCTC) were used to perform qRT-PCR on Arabidopsis. qUBC (F: ACAGCGAGAGAAAGTAGCAGA; qUBC-R: TTGATAAGAGCGGTCCATTTGAA) was used as an internal control for Arabidopsis. The CAT activity was calculated by measuring the decrease in absorbance at 240 nm with H2O2. POD activity and SOD was estimated according to the method developed by Gill [70,71]. Three biological replicates were set up for each of the above experimental treatments.

5. Conclusions

In this study, the C2H2 ZFP gene family was found to be distributed on 10 chromosomes of B. rapa. The genes were highly conserved during evolution. Quantitative analysis revealed that the Bra002528 gene was significantly expressed in the growth cones and roots of the cold-resistant rapeseed ‘Longyou-7’ as well as in different stress environments, indicating that this gene plays an important role in plant growth and development and the stress response. The survival rate of the trans-BrZAT12 gene in Arabidopsis was higher than that in the wild-type at low temperature, and the degree of leaf damage was lower than that in the wild-type. The activities of SOD, POD, and CAT were higher in the transgenic Arabidopsis than in the wild-type Arabidopsis. The expression of the BrZAT12 gene was induced by low temperature and played a role in the ABA signaling process involved in MAPKK, and it could enhance the cold resistance of plants. These results provide a theoretical basis for analyzing the cold-resistance mechanism of winter B. rapa and for increasing the cold resistance of northern winter rapeseed.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232012218/s1.

Author Contributions

Data curation, L.M., J.X., X.T., Y.P., Y.F. and L.L.; funding acquisition, L.M. and J.W.; investigation, W.W. and G.Y.; writing—original draft, L.M. and J.X.; writing—review and editing, L.M., X.L. and W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Youth Science and Technology Fund Program of Gansu Provincial Department of Science and Technology (21JR7RA835) and the China Agriculture Research System of MOF and MARA (CARS-12).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article and Supplementary Materials.

Acknowledgments

We would like to acknowledge all researchers in our laboratory for their help.

Conflicts of Interest

The authors have no conflict of interest to declare in relation to this article.

Abbreviations

B. rapaBrassica rapa
ArabidopsisArabidopsis thaliana
C2H2 ZFPC2H2 zinc-finger protein
B. napusBrassica napus
BrZAT12Brassica rapa zinc finger protein ZAT12-1
qRT-PCRReal-time fluorescence quantitative PCR
ABAAbscisic acid
HMMHidden Markov model
KsSynonymous substitution rate
KaNon-synonymous substitution rate
PODPeroxidase activity
SODSuperoxide dismutase activity
CATCatalase activity

References

  1. Wu, J.; Xu, X.; Liu, L.; Ma, L.; Pu, Y.; Wang, W.; Hua, X.; Song, J.; Liu, K.; Lu, G.; et al. Assembly of a Winter Turnip Rape (Brassica rapa L.) to Explore the Genetic Basis of Cold Tolerance. Front. Plant Sci. 2022, 13, 936958. [Google Scholar] [CrossRef] [PubMed]
  2. Ma, L.; Qi, W.; Bai, J.; Li, H.; Fang, Y.; Xu, J.; Xu, Y.; Zeng, X.; Pu, Y.; Wang, W.; et al. Genome-Wide Identification and Analysis of the Ascorbate Peroxidase (APX) Gene Family of Winter Rapeseed (Brassica Rapa L.) Under Abiotic Stress. Front. Genet. 2022, 12, 753624. [Google Scholar] [CrossRef]
  3. Han, G.; Lu, C.; Guo, J.; Qiao, Z.; Sui, N.; Qiu, N.; Wang, B. C2H2 Zinc Finger Proteins: Master Regulators of Abiotic Stress Responses in Plants. Front. Plant Sci. 2020, 11, 115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ding, Y.; Shi, Y.; Yang, S. Advances and Challenges in Uncovering Cold Tolerance Regulatory Mechanisms in Plants. New Phytol. 2019, 222, 1690–1704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhang, H.; Zhu, J.; Gong, Z.; Zhu, J.K. Abiotic Stress Responses in Plants. Nat. Rev. Genet. 2022, 23, 104–119. [Google Scholar] [CrossRef]
  6. Hassan, M.A.; Xiang, C.; Farooq, M.; Muhammad, N.; Yan, Z.; Hui, X.; Yuanyuan, K.; Bruno, A.K.; Lele, Z.; Jincai, L. Cold Stress in Wheat: Plant Acclimation Responses and Management Strategies. Front. Plant Sci. 2021, 12, 676884. [Google Scholar] [CrossRef] [PubMed]
  7. Chen, X.; Ding, Y.; Yang, Y.; Song, C.; Wang, B.; Yang, S.; Guo, Y.; Gong, Z. Protein Kinases in Plant Responses to Drought, Salt, and Cold Stress. J. Integr. Plant Biol. 2021, 63, 53–78. [Google Scholar] [CrossRef]
  8. Hwang, I.; Jung, H.J.; Park, J.I.; Yang, T.J.; Nou, I.S. Transcriptome Analysis of Newly Classified BZIP Transcription Factors of Brassica Rapa in Cold Stress Response. Genomics 2014, 104, 194–202. [Google Scholar] [CrossRef] [PubMed]
  9. Tian, Z.D.; Zhang, Y.; Liu, J.; Xie, C.H. Novel Potato C2H2-Type Zinc Finger Protein Gene, StZFP1, Which Responds to Biotic and Abiotic Stress, Plays a Role in Salt Tolerance. Plant Biol. 2010, 12, 689–697. [Google Scholar] [CrossRef]
  10. Xie, T.; Chen, C.; Li, C.; Liu, J.; Liu, C.; He, Y. Genome-Wide Investigation of WRKY Gene Family in Pineapple: Evolution and Expression Profiles during Development and Stress. BMC Genom. 2018, 19, 490. [Google Scholar] [CrossRef]
  11. Liu, Y.; Khan, A.R.; Gan, Y. C2H2 Zinc Finger Proteins Response to Abiotic Stress in Plants. Int. J. Mol. Sci. 2022, 23, 2730. [Google Scholar] [CrossRef]
  12. Selvaraj, M.G.; Jan, A.; Ishizaki, T.; Valencia, M.; Dedicova, B.; Maruyama, K.; Ogata, T.; Todaka, D.; Yamaguchi-Shinozaki, K.; Nakashima, K.; et al. Expression of the CCCH-Tandem Zinc Finger Protein Gene OsTZF5 under a Stress-Inducible Promoter Mitigates the Effect of Drought Stress on Rice Grain Yield under Field Conditions. Plant Biotechnol. J. 2020, 18, 1711–1721. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wang, K.; Ding, Y.; Cai, C.; Chen, Z.; Zhu, C. The Role of C2H2 Zinc Finger Proteins in Plant Responses to Abiotic Stresses. Physiol. Plant. 2019, 165, 690–700. [Google Scholar] [CrossRef] [PubMed]
  14. Wei, K.; Pan, S.; Li, Y. Functional Characterization of Maize C2H2 Zinc-Finger Gene Family. Plant Mol. Biol. Report. 2016, 34, 761–776. [Google Scholar] [CrossRef]
  15. Englbrecht, C.C.; Schoof, H.; Böhm, S. Conservation, Diversification and Expansion of C2H2 Zinc Finger Proteins in the Arabidopsis Thaliana Genome. BMC Genom. 2004, 5, 39. [Google Scholar] [CrossRef] [Green Version]
  16. Liu, Z.; Coulter, J.A.; Li, Y.; Zhang, X.; Meng, J.; Zhang, J.; Liu, Y. Genome-Wide Identification and Analysis of the Q-Type C2H2 Gene Family in Potato (Solanum tuberosum L.). Int. J. Biol. Macromol. 2020, 153, 327–340. [Google Scholar] [CrossRef] [PubMed]
  17. Cui, H.; Chen, J.; Liu, M.; Zhang, H.; Zhang, S.; Liu, D.; Chen, S. Genome-Wide Analysis of C2H2 Zinc Finger Gene Family and Its Response to Cold and Drought Stress in Sorghum [Sorghum Bicolor (L.) Moench]. Int. J. Mol. Sci. 2022, 23, 5571. [Google Scholar] [CrossRef] [PubMed]
  18. Zhao, T.; Wu, T.; Zhang, J.; Wang, Z.; Pei, T.; Yang, H.; Li, J.; Xu, X. Genome-Wide Analyses of the Genetic Screening of C2H2-Type Zinc Finger Transcription Factors and Abiotic and Biotic Stress Responses in Tomato (Solanum Lycopersicum) Based on RNA-Seq Data. Front. Genet. 2020, 11, 540. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, Q.; Wang, Z.; Xu, X.; Zhang, H.; Li, C. Genome-Wide Analysis of C2H2 Zinc-Finger Family Transcription Factors and Their Responses to Abiotic Stresses in Poplar (Populus Trichocarpa). PLoS ONE 2015, 10, e0134753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Vogel, J.T.; Zarka, D.G.; Van Buskirk, H.A.; Fowler, S.G.; Thomashow, M.F. Roles of the CBF2 and ZAT12 Transcription Factors in Configuring the Low Temperature Transcriptome of Arabidopsis. Plant J. 2005, 41, 195–211. [Google Scholar] [CrossRef] [PubMed]
  21. Chandra Rai, A.; Singh, M.; Shah, K.; Rai, A.C.; Singh, I.; Singh, M.; Shah, K.; Chandra Rai, A.; Singh, M.; Shah, K. Effect of Water Withdrawal on Formation of Free Radical, Proline Accumulation and Activities of Antioxidant Enzymes in ZAT12-Transformed Transgenic Tomato Plants. Plant Physiol. Biochem. 2012, 61, 108–114. [Google Scholar] [CrossRef] [PubMed]
  22. Yu, G.H.; Jiang, L.L.; Ma, X.F.; Xu, Z.S.; Liu, M.M.; Shan, S.G.; Cheng, X.G. A Soybean C2H2-Type Zinc Finger Gene GmZF1 Enhanced Cold Tolerance in Transgenic Arabidopsis. PLoS ONE 2014, 9, e109399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chong, L.; Guo, P.; Zhu, Y. Mediator Complex: A Pivotal Regulator of Aba Signaling Pathway and Abiotic Stress Response in Plants. Int. J. Mol. Sci. 2020, 21, 7755. [Google Scholar] [CrossRef] [PubMed]
  24. Ku, Y.S.; Sintaha, M.; Cheung, M.Y.; Lam, H.M. Plant Hormone Signaling Crosstalks between Biotic and Abiotic Stress Responses. Int. J. Mol. Sci. 2018, 19, 3206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Faraji, S.; Rasouli, S.H.; Kazemitabar, S.K. Genome-Wide Exploration of C2H2 Zinc Finger Family in Durum Wheat (Triticum Turgidum Ssp. Durum): Insights into the Roles in Biological Processes Especially Stress Response. BioMetals 2018, 31, 1019–1042. [Google Scholar] [CrossRef]
  26. Fadda, A.; Syed, N.; Mackeh, R.; Papadopoulou, A.; Suzuki, S.; Jithesh, P.V.; Kino, T. Genome-Wide Regulatory Roles of the C2H2-Type Zinc Finger Protein ZNF764 on the Glucocorticoid Receptor. Sci. Rep. 2017, 7, 41598. [Google Scholar] [CrossRef] [Green Version]
  27. Chalhoub, B.; Denoeud, F.; Liu, S.; Parkin, I.A.P.; 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]
  28. Lawrence, S.D.; Novak, N.G. Comparative Analysis of the Genetic Variability within the Q-Type C2H2 Zinc-Finger Transcription Factors in the Economically Important Cabbage, Canola and Chinese Cabbage Genomes. Hereditas 2018, 155, 29. [Google Scholar] [CrossRef] [Green Version]
  29. Lawrence, S.D.; Novak, N.G. Mining the Brassica Oleracea Genome for Q-Type C2H2 Zinc Finger Transcription Factor Genes. Plant Mol. Biol. Report. 2015, 33, 1611–1617. [Google Scholar] [CrossRef]
  30. Schranz, M.E.; Lysak, M.A.; Mitchell-Olds, T. The ABC’s of Comparative Genomics in the Brassicaceae: Building Blocks of Crucifer Genomes. Trends Plant Sci. 2006, 11, 535–542. [Google Scholar] [CrossRef]
  31. Liu, S.; Liu, Y.; Yang, X.; Tong, C.; Edwards, D.; Parkin, I.A.P.; Zhao, M.; Ma, J.; Yu, J.; Huang, S.; et al. The Brassica Oleracea Genome Reveals the Asymmetrical Evolution of Polyploid Genomes. Nat. Commun. 2014, 5, 3930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Xu, J.; Zhao, Y.; Hou, X.; Tao, X.; Sun, B.; Zeng, R.; Zhu, M.; Sun, W.; Wu, J.; Liu, L.; et al. Cloning and Functional Analysis of Zinc Finger Protein Gene BrZFP1 in Brassica Campestris. J. Agric. Biotechnol. 2022, 30, 236–248. [Google Scholar] [CrossRef]
  33. Davletova, S.; Schlauch, K.; Coutu, J.; Mittler, R. The Zinc-Finger Protein Zat12 Plays a Central Role in Reactive Oxygen and Abiotic Stress Signaling in Arabidopsis. Plant Physiol. 2005, 139, 847–856. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Ma, X.; Liang, W.; Gu, P.; Huang, Z. Salt Tolerance Function of the Novel C2H2-Type Zinc Finger Protein TaZNF in Wheat. Plant Physiol. Biochem. 2016, 106, 129–140. [Google Scholar] [CrossRef] [PubMed]
  35. Han, Y.C.; Fu, C.C. Cold-Inducible MaC2H2s Are Associated with Cold Stress Response of Banana Fruit via Regulating MaICE1. Plant Cell Rep. 2019, 38, 673–680. [Google Scholar] [CrossRef]
  36. Kim, J.C.; Lee, S.H.; Cheong, Y.H.; Yoo, C.M.; Lee, S.I.; Chun, H.J.; Yun, D.J.; Hong, J.C.; Lee, S.Y.; Lim, C.O.; et al. A Novel Cold-Inducible Zinc Finger Protein from Soybean, SCOF-1, Enhances Cold Tolerance in Transgenic Plants. Plant J. 2001, 25, 247–259. [Google Scholar] [CrossRef]
  37. Kim, Y.H.; Kim, M.D.; Park, S.C.; Yang, K.S.; Jeong, J.C.; Lee, H.S.; Kwak, S.S. SCOF-1-Expressing Transgenic Sweetpotato Plants Show Enhanced Tolerance to Low-Temperature Stress. Plant Physiol. Biochem. 2011, 49, 1436–1441. [Google Scholar] [CrossRef]
  38. Kim, Y.H.; Kim, M.D.; Park, S.C.; Jeong, J.C.; Kwak, S.S.; Lee, H.S. Transgenic Potato Plants Expressing the Cold-Inducible Transcription Factor SCOF-1 Display Enhanced Tolerance to Freezing Stress. Plant Breed. 2016, 135, 513–518. [Google Scholar] [CrossRef]
  39. Huang, J.; Sun, S.; Xu, D.; Lan, H.; Sun, H.; Wang, Z.; Bao, Y.; Wang, J.; Tang, H.; Zhang, H. A TFIIIA-Type Zinc Finger Protein Confers Multiple Abiotic Stress Tolerances in Transgenic Rice (Oryza sativa L.). Plant Mol. Biol. 2012, 80, 337–350. [Google Scholar] [CrossRef]
  40. Yilmaz, S.C.; Mittler, R. The Zinc Finger Network of Plants. Cell. Mol. Life Sci. 2008, 65, 1150–1160. [Google Scholar] [CrossRef]
  41. Kiełbowicz, M.A. Involvement of Plant C2H2-Type Zinc Finger Transcription Factors in Stress Responses. Plant Sci. 2012, 185–186, 78–85. [Google Scholar] [CrossRef]
  42. Huang, J.; Sun, S.J.; Xu, D.Q.; Yang, X.; Bao, Y.M.; Wang, Z.F.; Tang, H.J.; Zhang, H. Increased Tolerance of Rice to Cold, Drought and Oxidative Stresses Mediated by the Overexpression of a Gene That Encodes the Zinc Finger Protein ZFP245. Biochem. Biophys. Res. Commun. 2009, 389, 556–561. [Google Scholar] [CrossRef]
  43. Huang, J.; Wang, J.F.; Wang, Q.H.; Zhang, H.S. Identification of a Rice Zinc Finger Protein Whose Expression Is Transiently Induced by Drought, Cold but Not by Salinity and Abscisic Acid. DNA Seq.-J. DNA Seq. Mapp. 2005, 16, 130–136. [Google Scholar] [CrossRef]
  44. Zhang, H.; Liu, Y.; Wen, F.; Yao, D.; Wang, L.; Guo, J.; Ni, L.; Zhang, A.; Tan, M.; Jiang, M. A Novel Rice C2H2-Type Zinc Finger Protein, ZFP36, Is a Key Player Involved in Abscisic Acid-Induced Antioxidant Defence and Oxidative Stress Tolerance in Rice. J. Exp. Bot. 2014, 65, 5795–5809. [Google Scholar] [CrossRef]
  45. Huang, L.; Zhang, M.; Jia, J.; Zhao, X.; Huang, X.; Ji, E.; Ni, L.; Jiang, M. An Atypical Late Embryogenesis Abundant Protein OsLEA5 Plays a Positive Role in ABA-Induced Antioxidant Defense in Oryza sativa L. Plant Cell Physiol. 2018, 59, 916–929. [Google Scholar] [CrossRef]
  46. Brumbarova, T.; Le, C.T.T.; Ivanov, R.; Bauer, P. Regulation of ZAT12 Protein Stability: The Role of Hydrogen Peroxide. Plant Signal. Behav. 2016, 11, e1137408. [Google Scholar] [CrossRef] [Green Version]
  47. Ben Daniel, B.H.; Cattan, E.; Wachtel, C.; Avrahami, D.; Glick, Y.; Malichy, A.; Gerber, D.; Miller, G. Identification of Novel Transcriptional Regulators of Zat12 Using Comprehensive Yeast One-Hybrid Screens. Physiol. Plant. 2016, 157, 422–441. [Google Scholar] [CrossRef]
  48. Le, C.T.T.; Brumbarova, T.; Ivanov, R.; Stoof, C.; Weber, E.; Mohrbacher, J.; Fink-Straube, C.; Bauer, P. Zinc Finger of Arabidopsis Thaliana12 (Zat12) Interacts with Fer-like Iron Deficiency-Induced Transcription Factor (FIT) Linking Iron Deficiency and Oxidative Stress Responses. Plant Physiol. 2016, 170, 540–557. [Google Scholar] [CrossRef]
  49. Lingam, S.; Mohrbacher, J.; Brumbarova, T.; Potuschak, T.; Fink-Straube, C.; Blondet, E.; Genschik, P.; Bauer, P. Interaction between the BHLH Transcription Factor FIT and ETHYLENE INSENSITIVE3/ETHYLENE INSENSITIVE3-LIKE1 Reveals Molecular Linkage between the Regulation of Iron Acquisition and Ethylene Signaling in Arabidopsis. Plant Cell 2011, 23, 1815–1829. [Google Scholar] [CrossRef] [Green Version]
  50. Yoo, S.D.; Cho, Y.H.; Tena, G.; Xiong, Y.; Sheen, J. Dual Control of Nuclear EIN3 by Bifurcate MAPK Cascades in C2H4 Signalling. Nature 2008, 451, 789–795. [Google Scholar] [CrossRef]
  51. Peng, J.; Li, Z.; Wen, X.; Li, W.; Shi, H.; Yang, L.; Zhu, H.; Guo, H. Salt-Induced Stabilization of EIN3/EIL1 Confers Salinity Tolerance by Deterring ROS Accumulation in Arabidopsis. PLoS Genet. 2014, 10, e1004664. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Rizhsky, L.; Davletova, S.; Liang, H.; Mittler, R. The Zinc Finger Protein Zat12 Is Required for Cytosolic Ascorbate Peroxidase 1 Expression during Oxidative Stress in Arabidopsis. J. Biol. Chem. 2004, 279, 11736–11743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ciftci-Yilmaz, S.; Morsy, M.R.; Song, L.; Coutu, A.; Krizek, B.A.; Lewis, M.W.; Warren, D.; Cushman, J.; Connolly, E.L.; Mittler, R. The EAR-Motif of the Cys2/His2-Type Zinc Finger Protein Zat7 Plays a Key Role in the Defense Response of Arabidopsis to Salinity Stress. J. Biol. Chem. 2007, 282, 9260–9268. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Liu, X.M.; Nguyen, X.C.; Kim, K.E.; Han, H.J.; Yoo, J.; Lee, K.; Kim, M.C.; Yun, D.J.; Chung, W.S. Phosphorylation of the Zinc Finger Transcriptional Regulator ZAT6 by MPK6 Regulates Arabidopsis Seed Germination under Salt and Osmotic Stress. Biochem. Biophys. Res. Commun. 2013, 430, 1054–1059. [Google Scholar] [CrossRef]
  55. Nguyen, X.C.; Kim, S.H.; Lee, K.; Kim, K.E.; Liu, X.M.; Han, H.J.; Hoang, M.H.T.; Lee, S.W.; Hong, J.C.; Moon, Y.H.; et al. Identification of a C2H2-Type Zinc Finger Transcription Factor (ZAT10) from Arabidopsis as a Substrate of MAP Kinase. Plant Cell Rep. 2012, 31, 737–745. [Google Scholar] [CrossRef]
  56. Finn, R.D.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Mistry, J.; Mitchell, A.L.; Potter, S.C.; Punta, M.; Qureshi, M.; Sangrador-Vegas, A.; et al. The Pfam Protein Families Database: Towards a More Sustainable Future. Nucleic Acids Res. 2016, 44, D279–D285. [Google Scholar] [CrossRef]
  57. Mistry, J.; Finn, R.D.; Eddy, S.R.; Bateman, A.; Punta, M. Challenges in Homology Search: HMMER3 and Convergent Evolution of Coiled-Coil Regions. Nucleic Acids Res. 2013, 41, e121. [Google Scholar] [CrossRef] [Green Version]
  58. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional Classification of Proteins via Subfamily Domain Architectures. Nucleic Acids Res. 2017, 45, D200–D203. [Google Scholar] [CrossRef] [Green Version]
  59. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
  60. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.; Li, W.W.; Noble, W.S. MEME Suite: Tools for Motif Discovery and Searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  61. Lescot, M.; Déhais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van De Peer, Y.; Rouzé, P.; Rombauts, S. PlantCARE, a Database of Plant Cis-Acting Regulatory Elements and a Portal to Tools for in Silico Analysis of Promoter Sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef] [PubMed]
  62. Luo, X.; Tu, X.; Ding, Y.; Gao, G.; Deng, M. Expectation Pooling: An Effective and Interpretable Pooling Method for Predicting DNA-Protein Binding. Bioinformatics 2020, 36, 1405–1412. [Google Scholar] [CrossRef] [PubMed]
  63. Jin, J.; Zhang, H.; Kong, L.; Gao, G.; Luo, J. PlantTFDB 3.0: A Portal for the Functional and Evolutionary Study of Plant Transcription Factors. Nucleic Acids Res. 2014, 42, 1182–1187. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Subramanian, B.; Gao, S.; Lercher, M.J.; Hu, S.; Chen, W.H. Evolview v3: A Webserver for Visualization, Annotation, and Management of Phylogenetic Trees. Nucleic Acids Res. 2019, 47, W270–W275. [Google Scholar] [CrossRef] [PubMed]
  66. Voorrips, R.E. Mapchart: Software for the Graphical Presentation of Linkage Maps and QTLs. J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2-ΔΔCT Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  68. Ma, L.; Coulter, J.A.; Liu, L.; Zhao, Y.; Chang, Y.; Pu, Y.; Zeng, X.; Xu, Y.; Wu, J.; Fang, Y.; et al. Transcriptome Analysis Reveals Key Cold-Stress-Responsive Genes in Winter Rapeseed (Brassica Rapa L.). Int. J. Mol. Sci. 2019, 20, 1071. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Fang, Y.; Coulter, J.A.; Wu, J.; Liu, L.; Li, X.; Dong, Y.; Ma, L.; Pu, Y.; Sun, B.; Niu, Z.; et al. Identification of Differentially Expressed Genes Involved in Amino Acid and Lipid Accumulation of Winter Turnip Rape (Brassica Rapa L.) in Response to Cold Stress. PLoS ONE 2021, 16, e0245494. [Google Scholar] [CrossRef] [PubMed]
  70. Quiroga, M.; Guerrero, C.; Botella, M.A.; Barceló, A.; Amaya, I.; Medina, M.I.; Alonso, F.J.; De Forchetti, S.M.; Tigier, H.; Valpuesta, V. A Tomato Peroxidase Involved in the Synthesis of Lignin and Suberin. Plant Physiol. 2000, 122, 1119–1127. [Google Scholar] [CrossRef] [PubMed]
  71. Gill, T.; Sreenivasulu, Y.; Kumar, S.; Singh Ahuja, P. Over-Expression of Superoxide Dismutase Exhibits Lignification of Vascular Structures in Arabidopsis Thaliana. J. Plant Physiol. 2010, 167, 757–760. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 12218 g001 550
Figure 1. Chromosome distribution of C2H2 ZFP family members in B. rapa. The 141 members are distributed over 10 chromosomes.
Figure 1. Chromosome distribution of C2H2 ZFP family members in B. rapa. The 141 members are distributed over 10 chromosomes.
Ijms 23 12218 g001
Ijms 23 12218 g002 550
Figure 2. Phylogenetic tree of the C2H2 ZFP family. The 16 subfamilies are colored differently. Red stars represent B. rapa, blue circles represent B. napus, and purple triangles represent Arabidopsis. The abbreviations indicate the ZFP gene family known in Arabidopsis.
Figure 2. Phylogenetic tree of the C2H2 ZFP family. The 16 subfamilies are colored differently. Red stars represent B. rapa, blue circles represent B. napus, and purple triangles represent Arabidopsis. The abbreviations indicate the ZFP gene family known in Arabidopsis.
Ijms 23 12218 g002
Ijms 23 12218 g003 550
Figure 3. Analysis of the gene structure of the C2H2 ZFP gene family in B. rapa. (A) The 141 members are assigned to 12 sub-clusters. (B) Each colored box represents a motif, and the gray line represents a non-conserved sequence. (C) Exons and introns of each subgroup are indicated by green boxes and grey lines, respectively. (D) The specific colored boxes indicate the different cis-acting elements.
Figure 3. Analysis of the gene structure of the C2H2 ZFP gene family in B. rapa. (A) The 141 members are assigned to 12 sub-clusters. (B) Each colored box represents a motif, and the gray line represents a non-conserved sequence. (C) Exons and introns of each subgroup are indicated by green boxes and grey lines, respectively. (D) The specific colored boxes indicate the different cis-acting elements.
Ijms 23 12218 g003
Ijms 23 12218 g004 550
Figure 4. C2H2 ZFP gene family collinearity analysis. (A) The gene names are marked outwards on the corresponding chromosome. Duplicate gene pairs are highlighted with red lines. (B) The red, green, and orange indicate the chromosomes of Arabidopsis, B. rapa, and B. napus, respectively. Duplicate gene pairs are highlighted with blue connecting lines. (C) Ka, Ks, and Ka/Ks values of homologous gene pairs among species.
Figure 4. C2H2 ZFP gene family collinearity analysis. (A) The gene names are marked outwards on the corresponding chromosome. Duplicate gene pairs are highlighted with red lines. (B) The red, green, and orange indicate the chromosomes of Arabidopsis, B. rapa, and B. napus, respectively. Duplicate gene pairs are highlighted with blue connecting lines. (C) Ka, Ks, and Ka/Ks values of homologous gene pairs among species.
Ijms 23 12218 g004
Ijms 23 12218 g005 550
Figure 5. Expression of C2H2 ZFP genes in different varieties under cold and freezing stress. Heat maps were expressed using log2 values for each gene. The color scale represents the relative expression levels from low (blue) to high (red).
Figure 5. Expression of C2H2 ZFP genes in different varieties under cold and freezing stress. Heat maps were expressed using log2 values for each gene. The color scale represents the relative expression levels from low (blue) to high (red).
Ijms 23 12218 g005
Ijms 23 12218 g006 550
Figure 6. Heatmap of the expression of C2H2 ZFP genes in B. rapa under abiotic stress. (A–D) represent cold, high temperature, drought, and salt stress, respectively. Y, Z, J, and G in the column labels represent leaf, growth cone, stem, and root tissues of rapeseed, respectively.
Figure 6. Heatmap of the expression of C2H2 ZFP genes in B. rapa under abiotic stress. (A–D) represent cold, high temperature, drought, and salt stress, respectively. Y, Z, J, and G in the column labels represent leaf, growth cone, stem, and root tissues of rapeseed, respectively.
Ijms 23 12218 g006
Ijms 23 12218 g007 550
Figure 7. Phenotype, expression, and physiological activity of BrZAT12 transgenic Arabidopsis after low-temperature treatment. (A,B) Phenotype of BrZAT12 transgenic Arabidopsis after low-temperature treatment. (C) Survival rate of Arabidopsis plants after low-temperature treatment. (D) Quantitative RT-PCR test of BrZAT12 transgenic Arabidopsis at low temperature. (E–G) SOD, POD, and CAT activities, respectively. WT: wild-type Arabidopsis, OE: BrZAT12 transgenic Arabidopsis. a–f represent significant differences between treatments at p ≤ 0.05.
Figure 7. Phenotype, expression, and physiological activity of BrZAT12 transgenic Arabidopsis after low-temperature treatment. (A,B) Phenotype of BrZAT12 transgenic Arabidopsis after low-temperature treatment. (C) Survival rate of Arabidopsis plants after low-temperature treatment. (D) Quantitative RT-PCR test of BrZAT12 transgenic Arabidopsis at low temperature. (E–G) SOD, POD, and CAT activities, respectively. WT: wild-type Arabidopsis, OE: BrZAT12 transgenic Arabidopsis. a–f represent significant differences between treatments at p ≤ 0.05.
Ijms 23 12218 g007
Ijms 23 12218 g008 550
Figure 8. Expression of BrZAT12 and activities of SOD and APX in B. rapa under different treatments. (A–C) Room-temperature treatment and (D–F) low-temperature treatment. a–e represent significant differences between treatments at p ≤ 0.05.
Figure 8. Expression of BrZAT12 and activities of SOD and APX in B. rapa under different treatments. (A–C) Room-temperature treatment and (D–F) low-temperature treatment. a–e represent significant differences between treatments at p ≤ 0.05.
Ijms 23 12218 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, L.; Xu, J.; Tao, X.; Wu, J.; Wang, W.; Pu, Y.; Yang, G.; Fang, Y.; Liu, L.; Li, X.; et al. Genome-Wide Identification of C2H2 ZFPs and Functional Analysis of BRZAT12 under Low-Temperature Stress in Winter Rapeseed (Brassica rapa). Int. J. Mol. Sci. 2022, 23, 12218. https://doi.org/10.3390/ijms232012218

AMA Style

Ma L, Xu J, Tao X, Wu J, Wang W, Pu Y, Yang G, Fang Y, Liu L, Li X, et al. Genome-Wide Identification of C2H2 ZFPs and Functional Analysis of BRZAT12 under Low-Temperature Stress in Winter Rapeseed (Brassica rapa). International Journal of Molecular Sciences. 2022; 23(20):12218. https://doi.org/10.3390/ijms232012218

Chicago/Turabian Style

Ma, Li, Jia Xu, Xiaolei Tao, Junyan Wu, Wangtian Wang, Yuanyuan Pu, Gang Yang, Yan Fang, Lijun Liu, Xuecai Li, and et al. 2022. "Genome-Wide Identification of C2H2 ZFPs and Functional Analysis of BRZAT12 under Low-Temperature Stress in Winter Rapeseed (Brassica rapa)" International Journal of Molecular Sciences 23, no. 20: 12218. https://doi.org/10.3390/ijms232012218

APA Style

Ma, L., Xu, J., Tao, X., Wu, J., Wang, W., Pu, Y., Yang, G., Fang, Y., Liu, L., Li, X., & Sun, W. (2022). Genome-Wide Identification of C2H2 ZFPs and Functional Analysis of BRZAT12 under Low-Temperature Stress in Winter Rapeseed (Brassica rapa). International Journal of Molecular Sciences, 23(20), 12218. https://doi.org/10.3390/ijms232012218

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