The catalase gene family in cucumber: genome-wide identification and organization

Abstract

Catalase (CAT) is a common antioxidant enzyme in almost all living organisms. Currently, detailed reports on cucumber (Cucumis sativus L.) CAT (CsCAT) genes and tissue expression profiling are limited. In the present study, four candidate CsCAT genes were identified in cucumber. Phylogenetic analysis indicated that CsCAT1-CsCAT3 are closely related to Arabidopsis AtCAT1-AtCAT3, but no obvious counterpart was observed for CsCAT4. Intron/exon structure analysis revealed that only one of the 15 positions was completely conserved. Motif analysis showed that, unlike the CAT genes of other species, none of CsCAT genes contained all 10 motifs. Expression data showed that transcripts of all of the CsCAT genes, except CsCAT4, were detected in five tissues. Moreover, their transcription levels displayed differences under different stress treatments.

Keywords
Cucumis sativus L.; catalase; phylogenetic analysis; gene family; motif

Introduction

Catalase (CAT, EC 1.11.1.6) is a common antioxidant enzyme found in nearly all living organisms. It consists of four ferriprotoporphyrin groups per molecule. The typical catalase reaction involves the decomposition of two molecules of hydrogen peroxide (H₂O₂) to water (H₂O) and oxygen (O₂) (Scandalios, 1987). Catalase exists preferentially in peroxisomes. It is also detected in cytosol, mitochondria, and chloroplasts (Mullen et al., 1997). Plant catalase is usually encoded by a small gene family. Arabidopsis, tobacco, maize, and pumpkins were each found to contain three members of this family, and two each were identified in cottonseed and Hordeum vulgare (Kendall et al., 1983; Smith et al., 1984; Willekens et al., 1995; Frugoli et al., 1996; Guan and Scandalios, 1996; Esaka et al., 1997; Mullen et al., 1997; Iwamoto et al., 2000).

Catalase plays a critical role in plant development, defense, and senescence. Catalase scavenges H₂O₂ generated during mitochondrial electron transport, the β-oxidation of fatty acids, and most importantly photorespiratory oxidation (Yang and Poovaiah, 2002). Catalase is required when light-dependent plants produce H₂O₂ via photorespiration in the peroxisomes. Mutants that lack catalase activity are inviable under conditions in which photorespiration occurs or would occur. Catalase is essential for breaking down the concomitant generation of H₂O₂ produced during the first step of β-oxidation during the germination of oilseed plants. Transgenic tobacco overexpressing the E. coli katE gene is tolerant to high irradiance under drought conditions, though wild plants suffer severe photosynthesis-induced damage under the same conditions (Shikanai et al., 1998). Arabidopsis catalase 2 knock-out mutants (cat2) contain greater amounts of H₂O₂, and this is associated with the spread of necrotic lesions (Queval et al., 2007). In tobacco CAT1 antisense lines, catalase activity is severely reduced and necrotic lesions develop on some of the lower leaves as the H₂O₂ level increases (Takahashi et al., 1997). Exogenous application of sweet potato catalase SPCAT1 fusion protein delays or alleviates ethephon-mediated leaf senescence and H₂O₂ elevation, which suggests that SPCAT1 may play a physiological role in H₂O₂ homeostasis in leaves, as indicated by developmental cues and environmental stimuli (Chen et al., 2012).

The transcription level of various plant catalases is regulated both temporally and spatially and responds differentially to developmental and environmental stimuli (Guan and Scandalios, 1996; Zimmermann et al., 2006; Du et al., 2008). Tobacco CAT1 and CAT2 mRNA transcripts are detected in non-senescent leaves, but CAT3 is detected in both non-senescent and senescing leaves (Niewiadomska et al., 2009). Maize CAT1 and CAT3 are expressed in kernels throughout their development, but CAT2 is detected only during the later stages of kernel development (Acevedo and Scandalios, 1990; Acevedo et al., 1991). In Arabidopsis, CAT1 and CAT2 are mainly expressed in leaves and siliques, whereas CAT3 is mainly expressed in stem and root. The expression of CAT2 and CAT3 has been found to be controlled by circadian rhythms. CAT2 has been shown to be activated by cold and drought stresses, but CAT3 is mainly activated by abscisic acid, oxidative treatments, and senescence (Du et al., 2008; Hackenberg et al., 2013; Li et al., 2013; Zou et al., 2015). In hot pepper, different organ-specific expression patterns for CaCAT1-CaCAT3 are related to circadian rhythms and stress treatments (Lee and An, 2005). In Scots pine, CAT is involved in embryogenesis and cell death processes (Vuosku et al., 2015).

Cucumbers are an economically and nutritionally important vegetable crop cultivated worldwide. It belongs to the Cucurbitaceae family. The CAT genes exist as a small family in various plant species, but detailed reports on cucumber CAT genes and tissue- and stress-specific expression profiles are limited. The recent sequencing of the whole cucumber genome and deep sequencing have made genome-wide analysis of CAT genes in cucumbers possible (Huang et al., 2009). In the present study, a genome-wide gene identification, phylogenetic analysis, and expression analysis of CAT genes in cucumber and a comparative analysis of cucumber CAT genes with those of Arabidopsis, rice, and poplar were conducted. The present findings may serve as an important reference for functional studies of cucumber CAT genes.

Materials and Methods

Database search for cucumber CAT genes

A cucumber catalase sequence (GenBank number GU248529) was used as a query sequence for TBLASTN searches (Altschul et al., 1990) of CAT genes encoded in the cucumber genome. The cucumber genome sequence from the Cucumber Genome Initiative (CuGI) obtained and published by the Institute of Vegetables and Flowers of the Chinese Academy of Agricultural Sciences (IVF-CAAS) was used. Default parameters in the TBLASTN searches were wordsize two and extension 11. Redundant sequences with the same scaffold or chromosome location were removed from the data set.

To further confirm the existence of these hypothetical CAT genes, the cDNA sequences were first conceptually translated into amino-acid sequences and then searched for the catalase domain using the Simple Modular Architecture Research Tool (SMART) (Letunic et al., 2004).

Tree building

Multiple sequence alignments were performed on the CAT protein sequences using Clustal X (http://www.clustal.org/) with the default parameters (Larkin et al., 2007), and the alignments were then manually adjusted. A phylogenetic tree was constructed with the aligned CAT protein sequences using MEGA4 (http://www.megasoftware.net/mega4/mega.html) (Tamura et al., 2007) and the neighbor-joining (NJ) method with the following parameters: Poisson correction, pairwise deletion, and bootstrap (1,000 replicates). The constructed tree file was visualized using TreeView 1.6.6 (http://taxonomy.zoology.gla.ac.uk/rod/treeview.html) (Page, 1996).

Intron/exon structure, genome distribution, and segmental duplication

The DNA and cDNA sequences corresponding to each predicted gene from the cucumber genome and annotation database CuGI were downloaded, and the intron distribution pattern and splicing phase were then analyzed using the web-based bioinformatics tool Gene Structure Display Server (GSDS; http://gsds.cbi.pku.edu.cn/). To obtain information on the location of cucumber CAT gene, a map of the distribution of CsCAT genes throughout the cucumber genome was drawn using the MapInspect tool. To detect segmental duplication events, the 100 kb DNA segments flanking each CsCAT gene were analyzed. Regions in different linkage groups that contained ≥ 6 homologous pairs with < 25 nonhomologous gene interventions were defined as duplicated segments. A gene pair was considered tandemly duplicated when the genes were separated by < 5 intervening genes and shared ≥ 40% amino acid sequence similarity. BioEdit 5.0.6 (http://www.mbio.ncsu.edu/BioEdit/bioedit.html) (Hall, 1999) was used to analyze the CsCAT homologs for similarity on the phylogenetic tree.

Conserved motif prediction

Cucumber CAT proteins were searched for conserved motifs using the Multiple Em for Motif Elicitation to find similar sequences shared by these genes (MEME; http://meme.nbcr.net/meme/tools/meme) (Bailey and Elkan, 1994).

Reverse Transcription (RT)-PCR analysis of cucumber CAT genes

PCR primers were designed to avoid conserved regions. The primer sequences are shown in detail in Table S1. Seeds of the ‘Chinese long' 9930 inbred line, commonly used in modern cucumber breeding (Huang et al., 2009), were germinated and grown in trays containing soil mixture (peat:sand:pumice, 1:1:1, v/v/v). Plants were adequately watered and grown at day/night temperature cycles of 24/18°C with a 16 h photoperiod. For the salt, abscisic acid (ABA) and H₂O₂ treatments, seedlings were grown in liquid Murashige and Skoog (MS) medium containing 200 mM NaCl, 100 μM ABA and 10 mM H₂O₂, respectively. For the cold treatment, seedlings in the growth chamber were transferred to 4 °C under light conditions. The drought treated seedlings were desiccated.

After treatment for 0, 1, 3, and 6 h, whole seedlings were frozen in liquid nitrogen. Total RNAs of the roots, stems, leaves, flowers, and fruit at the 20 main-stem nodes stages collected under cold, salt, and drought treatment conditions were isolated using the TRIzol Reagent (Tiangen Biotech Co., Ltd. Beijing, China). DNase-treated RNA samples (0.5 μg) were reverse-transcribed using M-MLV reverse transcriptase (Invitogen). The reverse transcription reactions were performed at 42 °C for 1 h using 2 μM oligo-dT18 primer. Two microliters of the first strand cDNAs were used as templates for PCR amplification with a pair of gene-specific primers (Table S1). Samples were denatured for 5 min at 94 °C and then run for about 30 cycles of 30 s at 94 °C, 30 s at 55 °C and 45 s at 72 °C with a final extension of 5 min at 72 °C. For an accurate comparison and quantification of the transcript levels, the exponential phase of PCR amplification was determined by establishing the number of PCR cycles where the products exhibited an exponential phase: 25 cycles for actin PCR products and 30 cycles for CAT PCR products. The PCR products were separated by 1.5% agarose gels containing ethidium bromide and photographed under UV light. The results were confirmed using three independent biological replicates. The cucumber ACTIN cDNA fragment (161 bp) was used as an internal standard for normalizing cDNA concentration variations.

Results

Identification of four CsCAT genes

To identify the full complement of CAT genes in cucumber genomes, a cucumber catalase sequence (GenBank accession number GU248529) was used as a BLAST query sequence, and four candidate sequences were identified. The same four sequences were also obtained from Cucumber Genome Initiative (CuGI) using the HMM program based on multiple sequence alignment results of Arabidopsis CAT domain sequences. To further verify the reliability of these candidate sequences, SMART analysis was performed and all four sequences showed a typical CAT domain. The four cucumber CAT genes (CsCAT; Table 1) were subjected to further analysis. The number designation was based on the order of multiple sequence alignments. For study purposes, each was provisionally distinguished by a generic name, viz., CsCAT1-CsCAT4.

Table 1
The CsCAT genes identified in this study. The number designations are based on the order of the multiple sequence alignments.

Serial No.	Gene Name	Gene ID	Chromosamal location	CuGI(5'-3')	Length (a.a.)	Group name
1	CsCAT1	Csa014901	4	14579559-14582377	536	III
2	CsCAT2	Csa014900	4	14584724-14588175	458	I
3	CsCAT3	Csa013194	6	14161472-14164768	483	II
4	CsCAT4	Csa022770	Scaffold000416	9221-11455	744	II

Phylogenetic analysis of the CsCAT genes

To investigate the evolutionary relationships of cucumber CsCATs to those in other species, an unrooted phylogenetic tree using bootstrap analysis (1,000 replicates) was built from alignments of the complete protein sequences of four cucumber CsCATs, three rice OsCATs, three Arabidopsis AtCATs, and three poplar PtCATs (Figure 1A). The results showed that, like the findings observed in other species (Mhamdi et al., 2010a), these 13 CATs could be divided into three classes (I–III), containing five, four, and four genes, respectively (Figure 1A). All three poplar CAT genes (PtCAT1-PtCAT3) existed in the form of an exclusive cluster belonging to class I. Among the four cucumber CsCATs, CsCAT1 was grouped in class III, CsCAT2 in class I, and CsCAT3 and CsCAT4 in class II. Compared to genes from rice and poplar, it seemed that CsCAT1, CsCAT2 and CsCAT3 were more related to AtCAT1, AtCAT2, and AtCAT3 in Arabidopsis, respectively. CsCAT4 was at the basal position of the tree and far distant from the other clades.

Figure 1
Plant catalase genes. (A) Neighbor-joining tree for CAT sequences of cucumber and other species. (B) Distribution of conserved motifs in four cucumber CsCATs, three rice OsCATs, three Arabidopsis AtCATs, and three poplar PtCAT proteins identified using the MEME search tool. Each motif is represented by a colored box and a number. The order of the motifs corresponds to the position of motifs in individual protein sequences. Details regarding the motifs are given in Table S3.

Sequence analysis and multiple sequence alignment

The open reading frame length of the four CsCAT genes varied from 1377 (CsCAT2) to 2235 bp (CsCAT4), encoding polypeptides of 458–744 amino acids, with a predicted molecular weight range of 53.49–82.90 kDa. Multiple alignment demonstrated that all CsCAT proteins contained a highly conserved region of about 380 amino acid residues in their N-terminal portion and 60 amino acid residues in their C-terminal portion corresponding to the catalase and catalase-rel domains, respectively. The homologies ranged from 40% to 76%, with an average homology of only 57.33%, suggesting a high sequence diversity between cucumber CsCAT genes. Detailed information is given in Table S2. Multiple sequence alignment showed three conserved catalytic amino acid residues for CAT enzyme (His-71, Asn-144, and Tyr-354 in the crab sequence) to be completely conserved in all CsCATs except CsCAT2 (Figure S1). The catalase proximal active site signature FDRERIPERVVHAKGAGA (residues 60–77), was also found to be highly conserved. Another sequence, RLFSYNDTH (residues 350–358), representing the proximal heme-ligand signature was only conserved in CsCAT1-CsCAT3 (Figure S1).

The plant CAT enzyme activity is known to be mainly localized in the peroxisome, which was targeted by the particular peroxisomal targeting signal (PTS). The consensus amino acid sequence for putative peroxisomal targeting, Ser/Glu/Cys-Lys/Arg/His-Leu, was located at nine amino acids from the carboxy terminus (Subramani, 1998). The current alignment analysis results showed that the classical peroxisomal targeting sequence was only observable at the C-terminus of CsCAT2. Another internal consensus tripeptide PTS1-like motif (QKL/I/V) was reported by Kamigaki (Kamigaki et al., 2003). It is located upstream of the Ser/Glu/Cys-Lys/Arg/His-Leu motif and was also found at the same sites as CsCAT2 and CsCAT3 (Figure S1). The PTS1-like motif sequence can be attributed to the efficient importation of catalase into peroxisomes.

Structure and evolution of the CsCAT genes

The pattern of intron positioning can provide some clues regarding evolutionary history. To investigate intron number and positions, a comparison of the full-length cDNA sequences with the corresponding genomic DNA sequences was performed. A previous study revealed that the rice and Arabidopsis CAT genes contained six to seven introns (Mhamdi et al., 2010b). A similar situation was observed in cucumber CsCAT1–CsCAT3, which had five to eight introns. CsCAT4 was an exception to this rule, showing no introns (Figure 2A). Further analysis showed a total of 15 intron positions in CsCAT1–CsCAT3. Among them, only one position (1) was completely conserved in three genes (indicated with a solid triangle). Of the remaining 14 positions, four positions (2, 10, 12, 14) were highly conserved between CsCAT1 and CsCAT2. However, no position was conserved between CsCAT1 and CsCAT3 or between CsCAT2 and CsCAT3. The intron phase analysis demonstrated that all the positions showed phase 0 (splicing occurring after the third nucleotide of the codon), except positions 12 and 13 with phase 2 (splicing occurring after the second nucleotide of the codon) and position two with phase 1.

Figure 2
Sequence analysis of plant catalase genes. (A) Nucleotide sequences surrounding splice sites in three cucumber CsCAT genes and the consensus sequence for the proto-splice site. Amino acid residues and codons (solid triangles) that are split by introns are also shown. (B) Amino acid sequences flanking the intron splice sites (underlined) in three cucumber CsCATs, two rice OsCATs, and three Arabidopsis AtCATs are shown. The amino acid Thr that is split by intron, which is widely conserved in the CAT family, is underlined.

In Arabidopsis, the CAT genes had one intron that splits the Thr in the consensus sequence of SS(M/L)T(V/I)G(P/T/E)RGP (Frugoli et al., 1996). The current analysis showed that CsCAT1 and CsCAT2 contained one intron at the conserved Thr site, as expected, but no such situation was observed in CsCAT3 (Figure 2B).

CsCAT conserved motif prediction

To investigate the characterized regions of cucumber CAT proteins, the online MEME motif search tool was used to analyze the distribution of the motifs in CsCATs in rice, Arabidopsis, and poplar. A total of 10 motifs, here called motifs 1–10, were identified (Figure 1B; Table S3). Among them, motifs 1–9 represented the catalase domain and motif 10 the catalase-rel region. It is worthy to note that the 10 motifs were all observed in rice, Arabidopsis, and poplar CATs, but no CsCAT gene contained all 10 motifs. For example, motif 1 was absent from CsCAT1, CsCAT3, and CsCAT4, motif 2 was absent from CsCAT2, and motifs 8–10 were absent from CsCAT4.

Detection of 4 CsCAT genes in different tissues and stress conditions

RT-PCR analysis was performed on RNA from root, stem, leaves, flowers, and fruit to assess the expression patterns of the four CsCAT genes (Figure 3A). The RT-PCR results showed that CsCAT1–CsCAT3 are expressed in all tissues investigated, but no CsCAT4 expression was detected in any of these. CsCAT1 had expression patterns similar to CsCAT2, with a higher transcript signal in root, leaves, and fruit. CsCAT3 displayed noticeably high expression signals in stem, leaves, flower, and fruit. Catalase is known to play a vital role in response to a subset of biotic and abiotic stress reactions.

Figure 3
RT-PCR analysis of cucumber CsCAT genes, (A) in different tissues and, (B) under salt, cold, drought, H₂O_2-, and ABA stress conditions for 0, 1, 3, and 6 h. RT-PCR was performed using primers specific to the CsCAT genes. PCR products were run on 1.5% agarose gels. CsACTIN primers that generated a 161 bp product, which was used as the internal standards for each CsCAT gene. Sources of the samples are as follows: (R) root, (S) stem, (L) leaves, (F) flowers, and (Fr) fruit.

To investigate the response of CsCAT genes to various abiotic stresses, RT-PCR was performed on 7-day-old seedlings treated under five different sets of abiotic stress conditions (salt, cold, drought, H₂O₂, and ABA) (Figure 3B). The data indicated that CsCAT1–CsCAT3 expression levels were reduced or increased relative to controls in at least one of the stress conditions examined. The level of transcription of CsCAT1 was increased in ABA treatment and that of CsCAT2 was enhanced by ABA and in the sample 1 for H₂O₂, but levels of CsCAT1 and CsCAT2 were moderately diminished for salt and drought. CsCAT3 transcript levels were slightly increased under salt, drought, and ABA stress conditions. Its transcription levels were enhanced under salt and drought conditions but reduced under ABA stress conditions.

Discussion

Generally, plant catalases comprise a small gene family only. There are three catalases each in Arabidopsis thaliana, Nicotiana tabacum, and Zea mays and two each in Hordeum vulgare and cottonseed (Kendall et al., 1983; Smith et al., 1984; Willekens et al., 1995; Frugoli et al., 1996; Guan and Scandalios, 1996; Esaka et al., 1997; Mullen et al., 1997; Iwamoto et al., 2000). The current study showed that cucumber contains four CAT genes (CsCAT1-CsCAT4), i.e. at least one more than those of the above mentioned species. The genomic distribution analysis showed that CsCAT1–CsCAT3 are located on chromosome scaffolds, whereas CsCAT4 was found on an unassembled sequence, scaffold000416. Hence, if CsCAT4 is removed from the dataset, the number of CAT genes in cucumber genome would drop to three, which makes it comparable to those in other reported species.

Sequence analysis showed that the catalytic amino acids, proximal heme-ligand signature sequence and catalase proximal active site signature were highly conserved among CsCAT1–CsCAT3, suggesting that these genes are likely to have functions similar to those of catalases in other species. The predicted peroxisome targeting signal analysis showed that CsCAT2 and CsCAT3 contain the putative peroxisomal targeting signals S/E/C-K/R/H-L, QKL/I/V, or both, indicating that CsCAT2 and CsCAT3 could be peroxisomal catalases. Nonetheless, further experimental analysis is needed for confirmation.

The phylogenetic tree analysis indicated that cucumber CsCAT1, CsCAT2, and CsCAT3 are more closely related to AtCAT1, AtCAT2, and AtCAT3 in Arabidopsis, respectively, than to other counterparts in other species. AtCAT1 and AtCAT2 were found to be mainly expressed in leaves and siliques, but AtCAT3 was preferentially expressed in stems and roots (Du et al., 2008).

The current expression analysis showed that CsCAT1 and CsCAT2 are represented by transcripts in photosynthetic tissue, such as leaves, but also in root and fruit. CsCAT3 showed transcripts in flowers and fruit, and also in stem and leaves. Based on the expression pattern comparison, CsCATs exhibited an expression pattern similar to their Arabidopsis counterparts. This suggests that the three CsCAT genes might play similar roles in cucumber development.

Under different stress treatment conditions, CsCAT1 was only activated by abscisic acid treatment, whereas the abundance of its Arabidopsis counterpart, the CAT1 transcript, has been shown significantly increased by cold, drought, abscisic acid, and oxidative treatments (Du et al., 2008). CsCAT2 transcription was reduced in response to drought, whereas Arabidopsis CAT2 transcription was elevated (Du et al., 2008). Similarly, CsCAT3 transcription was enhanced under salt and drought stress conditions and reduced under ABA stress conditions, while Arabidopsis CAT3 was increased by abscisic acid and oxidative treatments (Du et al., 2008). Based on the analysis of different environment stimuli, obvious differences were observed between cucumber CsCAT and their Arabidopsis counterparts. This leads to conclude that the activation of cucumber CsCATs may differ in response to different abiotic and biotic stresses.

Intron analysis demonstrated that CsCAT1 and CsCAT2 each contained one conserved intron and Thr in the consensus sequence of SS(M/L)T(V/I)G(P/T/E)RGP as expected, but this was not observed in CsCAT3, suggesting that CsCAT3 might have evolved more recently. Motif analysis showed that none of the CsCAT genes contained all of the expected motifs. This makes them different from the CAT genes in other species, and suggests that some motifs might have been lost from cucumber after the divergence of monocots and dicots, or may have evolved solely in cucumbers after the divergence.

In summary, we performed extensive analyses of the four cucumber CAT genes and compared them to three rice OsCATs, three Arabidopsis AtCATs, and three poplar PtCATs. The 13 CAT genes clustered into three classes (I–III), which were in general agreement with reported results (Mhamdi et al., 2010a). The protein sequence analysis showed that CsCAT2 and CsCAT3 might be peroxisomal catalases. The expression analysis of the four cucumber CsCAT genes in root, stem, leaves, flowers, and fruit showed that all genes were expressed in at least one tissue. Furthermore, their expression patterns displayed differences when exposed to stress conditions (cold, salt, drought, H₂O₂, and ABA treatment). The comprehensive data collected here may be useful for future analysis of the biological functions of CAT family genes in cucumber growth, development, and responses to different stresses.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (31460522), the Natural Science Foundation of Jiangxi, China (20122BAB214017).

References

• Acevedo A and Scandalios JG (1990) Expression of the catalase and superoxide dismutase genes in mature pollen in maize. Theor Appl Genet 80:705-711.
• Acevedo A, Williamson JD and Scandalios JG (1991) Photoregulation of the Cat2 and Cat3 catalase genes in pigmented and pigment-deficient maize: The circadian regulation of Cat3 is superimposed on its quasi-constitutive expression in maize leaves. Genetics 127:601-607.
• Altschul S, Gish W, Miller W, Myers EW and Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215:403-410.
• Bailey TL and Elkan C (1994) Fitting a mixture model by expectation maximization to discover motifs in biopolymers. Proc Int Conf Intell Syst Mol Biol 2:28-36.
• Chen HJ, Wu SD, Huang GJ, Shen C, Afiyanti M, Li WJ and Lin YH (2012) Expression of a cloned sweet potato catalase SPCAT1 alleviates ethephon-mediated leaf senescence and H(2)O(2) elevation. J Plant Physiol 169:86-97.
• Du YY, Wang PC, Chen J and Song CP (2008) Comprehensive functional analysis of the catalase gene family in Arabidopsis thaliana J Integr Plant Biol 50:1318-1326.
• Esaka M, Yamada N, Kitabayashi M, Setoguchi Y, Tsugeki R, Kondo M and Nishimura M (1997) cDNA cloning and differential gene expression of three catalases in pumpkin. Plant Mol Biol 33:141-155.
• Frugoli JA, Zhong HH, Nuccio ML, McCourt P, McPeek MA, Thomas TL and McClung CR (1996) Catalase is encoded by a multigene family in Arabidopsis thaliana (L.) Heynh. Plant Physiol 112:327-336.
• Guan L and Scandalios JG (1996) Molecular evolution of maize catalases and their relationship to other eukaryotic and prokaryotic catalases. J Mol Evol 42:570-579.
• Hackenberg T, Juul T, Auzina A, Gwizdz S, Malolepszy A, Van Der Kelen K, Dam S, Bressendorff S, Lorentzen A, Roepstorff P, et al. (2013) Catalase and NO CATALASE ACTIVITY1 promote autophagy-dependent cell death in Arabidopsis. Plant Cell 25:4616-4626.
• Hall TA (1999) BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Seri 41:95-98.
• Huang S, Li R, Zhang Z, Li L, Gu X, Fan W, Lucas WJ, Wang X, Xie B, Ni P, et al. (2009) The genome of the cucumber, Cucumis sativus L. Nat Genet 41:1275-1281.
• Iwamoto M, Higo H and Higo K (2000) Differential diurnal expression of rice catalase genes: The 5'-flanking region of CatA is not sufficient for circadian control. Plant Sci 151:39-46.
• Kamigaki A, Mano S, Terauchi K, Nishi Y, Tachibe-Kinoshita Y, Nito K, Kondo M, Hayashi M, Nishimura M and Esaka M (2003) Identification of peroxisomal targeting signal of pumpkin catalase and the binding analysis with PTS1 receptor. Plant J 33:161-175.
• Kendall AC, Keys AJ, Turner JC, Lea PJ and Miflin BJ (1983) The isolation and characterisation of a catalase-deficient mutant of barley (Hordeum vulgare L.). Planta 159:505-511.
• Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al. (2007) Clustal W and Clustal X version 2.0. Bioinformatics 23:2947-2948.
• Lee SH and An CS (2005) Differential expression of three catalase genes in hot pepper (Capsicum annuum L.). Mol Cells 20:247-255.
• Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP and Bork P (2004) SMART 4.0: Towards genomic data integration. Nucleic Acids Res. 32:D142-D144.
• Li Y, Chen L, Mu J and Zuo J (2013) Lesion Simulating Disease1 interacts with catalases to regulate hypersensitive cell death in Arabidopsis. Plant Physiol 163:1059-1070.
• Mhamdi A, Queval G, Chaouch S, Vanderauwera S, Van Breusegem F and Noctor G (2010a) Catalase function in plants: A focus on Arabidopsis mutants as stress-mimic models. J Exp Bot 61:4197-4220.
• Mhamdi A, Hager J, Chaouch S, Queval G, Han Y, Taconnat L, Saindrenan P, Gouia H, Issakidis-Bourguet E, Renou JP, et al. (2010b) Arabidopsis Glutathione Reductase1 plays a crucial role in leaf responses to intracellular hydrogen peroxide and in ensuring appropriate gene expression through both salicylic acid and jasmonic acid signaling pathways. Plant Physiol 153:1144-1160.
• Mullen RT, Lee MS and Trelease RN (1997) Identification of the peroxisomal targeting signal for cottonseed catalase. Plant J 12:313-322.
• Niewiadomska E, Polzien L, Desel C, Rozpadek P, Miszalski Z and Krupinska, K (2009) Spatial patterns of senescence and development-dependent distribution of reactive oxygen species in tobacco (Nicotiana tabacum) leaves. J Plant Physiol 166:1057-1068.
• Page RD (1996) TreeView: An application to display phylogenetic trees on personal computers. Comput Appl Biosci 12:357-358.
• Queval G, Issakidis-Bourguet E, Hoeberichts FA, Vandorpe M, Gakiere B, Vanacker H, Miginiac-Maslow M, Van Breusegem F and Noctor G (2007) Conditional oxidative stress responses in the Arabidopsis photorespiratory mutant cat2 demonstrate that redox state is a key modulator of daylength-dependent gene expression, and define photoperiod as a crucial factor in the regulation of H2O2-induced cell death. Plant J 52:640-657.
• Scandalios JG (1987) The antioxidant enzyme genes Cat and Sod of maize: Regulation, functional significance, and molecular biology. Isozymes Curr Top Biol Med Res 14:19-44.
• Shikanai T, Takeda T, Yamauchi H, Sano S, Tomizawa KI, Yokota A and Shigeoka S (1998) Inhibition of ascorbate peroxidase under oxidative stress in tobacco having bacterial catalase in chloroplasts. FEBS Lett 428:47-51.
• Smith IK, Kendall AC, Keys AJ, Turner JC and Lea PJ (1984) Increased levels of glutathione in a catalase-deficient mutant of barley (Hordeum vulgare L.). Plant Sci Lett 37:29-33.
• Subramani S (1998) Components involved in peroxisome import, biogenesis, proliferation, turnover, and movement. Physiol Rev 78:171-188.
• Takahashi H, Chen Z, Du H, Liu Y and Klessig DF (1997) Development of necrosis and activation of disease resistance in transgenic tobacco plants with severely reduced catalase levels. Plant J 11:993-1005.
• Tamura K, Dudley J, Nei M and Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24:1596-1599.
• Vuosku J, Sutela S, Kestila J, Jokela A, Sarjala T and Haggman H (2015) Expression of catalase and retinoblastoma-related protein genes associates with cell death processes in Scots pine zygotic embryogenesis. BMC Plant Biol 15:e88.
• Willekens H, Inzé D, Van Montagu M and Van Camp W (1995) Catalases in plants. Mol Breed 1:207-228.
• Yang T and Poovaiah BW (2002) Hydrogen peroxide homeostasis: Activation of plant catalase by calcium/calmodulin. Proc Natl Acad Sci USA 99:4097-4102.
• Zimmermann P, Heinlein C, Orendi G and Zentgraf U (2006) Senescence-specific regulation of catalases in Arabidopsis thaliana (L.) Heynh. Plant Cell Environ 29:1049-1060.
• Zou JJ, Li XD, Ratnasekera D, Wang C, Liu WX, Song LF, Zhang WZ and Wu WH (2015) Arabidopsis calcium- dependent protein kinase8 and catalase3 function in abscisic acid-mediated signaling and H2O2 homeostasis in stomatal guard cells under drought stress. Plant Cell 27:1445-1460.

Supplementary material

The following online material is available for this article:

Figure S1 - Multiple alignment of the deduced catalase amino acid sequences.

Table S1 - Primers used in RT-PCR.

Table S2 - Sequence alignment between cucumber four CsCAT genes.

Table S3 - Summary of the conserved motifs (CMs) within the CsCAT family.

• Associate Editor: Adriana S. Hemerly

Publication Dates

• Publication in this collection
  25 July 2016
• Date of issue
  Jul-Sep 2016

History

• Received
  21 July 2015
• Accepted
  05 Jan 2016