- Research article
- Open Access
Genome-wide identification of cyclophilin genes in Gossypium hirsutum and functional characterization of a CYP with antifungal activity against Verticillium dahliae
BMC Plant Biology volume 19, Article number: 272 (2019)
Cyclophilins (CYPs), belonging to the peptidyl prolyl cis/trans isomerase (PPIase) superfamily, play important roles during plant responses to biotic and abiotic stresses.
Here, a total of 79 CYPs were identified in the genome of Gossypium hirsutum. Of which, 65 GhCYPs only contained one cyclophilin type PPIase domain, others 14 GhCYPs contain additional domains. A number of cis-acting elements related to phytohormone signaling were predicated in the upstream of GhCYPs ORF. The expression analysis revealed that GhCYPs were induced in response to cold, hot, salt, PEG and Verticillium dahliae infection. In addition, the functional importance of GhCYP-3 in Verticillium wilt resistance was also presented in this study. GhCYP-3 showed both cytoplasmic and nuclear localization. Overexpression of GhCYP-3 in Arabidopsis significantly improved Verticillium wilt resistance of the plants. Recombinant GhCYP-3 displayed PPIase activity and evident inhibitory effects on V. dahliae in vitro. Moreover, the extracts from GhCYP-3 transgenic Arabidopsis displayed significantly inhibit activity to conidia germinating and hyphal growth of V. dahliae.
Our study identified the family members of cotton CYP genes using bioinformatics tools. Differential expression patterns of GhCYPs under various abiotic stress and V. dahliae infection conditions provide a comprehensive understanding of the biological functions of candidate genes. Moreover, GhCYP-3 involved in the resistance of cotton to V. dahliae infection presumably through antifungal activity.
The allotetraploid Gossypium hirsutum L., the most important fiber crop, is cultivated worldwide because of its high yield . Environmental stresses, such as cold, drought, heat, salinity, various pests and pathogens, threaten cotton growth, yield and fiber quality. For example, an area of about 300 million hectares of cotton is subject to Verticillium wilt and the economic loss is more than RMB 12 billion every year in China. Verticillium wilt, the most serious disease to influence cotton production in China, caused by soil-borne fungus Verticillium dahliae. At present, few resistant varieties have been cultivated by traditional cross-breeding in G. hirsutum, which contributes 95% of the total cotton yield in the world . Genetic engineering using plant resistance genes is becoming an alternative to improve cotton resistance to V. dahliae.
Cyclophilins (CYPs), possessing peptidyl prolyl cis/trans isomerase (PPIase) activity, are classified in the immunophilin family of proteins . CYPs play important roles in various biological processes, including transcription regulation , protein folding , signal transduction  and ROS (reactive oxygen species) regulation . With the availability of whole genome sequencing, the identification and characterization of plant CYPs are carried out mainly in Arabidopsis thaliana (29 AtCYPs) , Oryza sativa (27 OsCYPs)  and Glycine max (62 GmCYPs) . The majority of studies reveal the involvement of plant CYPs mostly in different types of abiotic stress. For example, Arabidopsis CyPs showed evidence of response to wounding . Rice OsCYP19-4 showed over 10-fold up-regulation in response to cold. Overexpressing of OsCYP19-4 could enhance rice plants cold-resistance with significantly increased tiller and spike numbers, and consequently enhanced grain weight . Transgenic plants overexpressing OsCYP21-4 exhibited increased tolerance to salinity and hydrogen peroxide treatment . Ectopic expression of pigeon pea (Cajanus cajan L.) CYP, CcCYP, in Arabidopsis exhibited high-level tolerance against drought, salinity and extreme temperatures . Against biotic stress, especially against pathogen infection, only several plant CYPs have been studied in plant-pathogen system. Pepper cyclophilin (CACYP1) gene expression increased in response to Xanthomonas campestris pv. vesicatoria and Colletotrichum gloeosporioides . Fungal infection with Fusariumsolani f. speumartii increased the level of Solanum tuberosusm CyP gene StCyP mRNA in tubers . The expression of V. vinifera VviCyP was highly induced by Plasmopara viticola . In cotton, a cyclophilin-like gene GhCyp1 was cloned from G. hirsutum. Overexpression of GhCyp1 in transgenic tobacco plants conferred higher tolerance to salt stress and Pseudomonas syringae pv. tabaci infection compared with control plants .
In 2015, the genome of G. hirsutum L. acc. Texas Marker-1 (TM-1) was sequenced, more than 70, 000 protein-coding genes were predicted (NAU version 1.1) [19, 20]. Recently, an improved de novo–assembled genome for G. hirsutum L. acc. TM-1 were generated (NAU version 2.1) . The genome-sequencing project facilitates the survey of all CYP genes in cotton. In the present study, the CYP gene family members in G. hirsutum and their expression patterns under various abiotic stresses and on V. dahliae infection were systematically investigated. Furthermore, the function of GhCYP-3 was analyzed to reveal its role in cotton resistance to V. dahliae infection. Our study will enlighten the novel insights into the function of CYP genes in plant against multivariate stress responses in the future and provide more candidate genes for resistance breeding in cotton.
Up to 79 CYPs were identified in the genome of G. hirsutum TM-1
A local BLASTP search was performed with the Arabidopsis CYP proteins as query, which resulted in 79 CYP candidates from G. hirsutum NAU version 1.1, 74 CYPs from G. hirsutum JGI version 1.0 and 78 CYPs from G. hirsutum NAU version 2.1 (Table 1). These candidates were submitted to Pfam to confirm the existence of cyclophilin type PPIase domain (CLD, PF00160) and named GhCYP-1 to GhCYP-79. The characteristics of the individual CYP, including CDS length, protein length, molecular weight, and isoelectric point (pI) were presented in Table 1. The protein length varied from 69 amino acid (aa) residues (GhCYP-70) to 801 aa (GhCYP-55). The molecular weight ranged from 7.5 kDa (GhCYP-70) to 90.5 kDa (GhCYP-55), and the pI values ranged from 4.6 to 12.0. Most of the GhCYPs were expected to be in the cytoplasm. Also, some GhCYPs exhibited chloroplast, mitochondrial, nuclear and extracellular localization. Of the 79 GhCYPs, 65 only contained one CLD domain, but the remaining 14 GhCYPs contain additional domains, including tetratricopeptide-like repeats (TPR, PF00515, PF07719, PF13181, PF13414), Zinc finger (zf-CCHC, PF00098), RNA recognition motif (RRM, PF00076) and WD40 (PF00400) (Fig. 1).
Cis-elements potentially related to hormonal signal for GhCYPs
Here we surveyed the presence of cis-elements potentially related to the hormonal signal, in the -2 kb 5′ flanking region upstream to the start codon of these GhCYPs. In total six types of hormones related cis-elements in the promoters were predicted (Fig. 2). Of these GhCYPs, 66 GhCYPs had ethylene (ET) responsive element (ERE), 35 GhCYPs contained salicylic acid (SA) responsive element (TCA-element), 47 GhCYPs harbored abscisic acid (ABA) responsive element (ABRE), 38 GhCYPs possessed gibberellin (GA) responsive element (P-box; TATC-box), 47 GhCYPs contained methyl jasmonate (MeJA) responsive element (CGTCA-motif), and 26 possessed auxin responsive element (TGA-box; AuxRR-core). In total, 217 cis-elements related to ET, 150 cis-elements related to MeJA, 132 cis-elements related to ABA, 57 cis-elements related to GA, 43 cis-elements related to SA and 34 cis-elements related to auxin were identified in all GhCYPs (Fig. 2). The enrichment of hormone-responsive cis-elements in the upstream of these GhCYPs suggests that they are likely to be involved in plant responses to various hormone signal pathways.
Expression patterns of GhCYPs under various abiotic stresses
Expression profiles of GhCYPs were examined in roots of cotton plants under four different abiotic stress conditions using high-throughput RNA-seq data (Additional file 2: Table S2). The transcripts with low Fragments Per Kilobase of exon per Million fragments mapped (FPKM) were probably false assembly. Therefore, in this study GhCYPs with FPKM >10 and present in at least two samples were identified as potentially expressed transcripts. We focused on the significantly differentially expressed genes (fold change [FC] > 2 or FC < 0.5) in the various stresses. GhCYPs that showed significant differential expression at least one treatment time were listed in Fig. 3. In the cold treatment, 14 GhCYPs were commonly down-regulated significantly apart from GhCYP-49 at 1 hct (hours after cold treated) and 3 hct and GhCYP-52 at 6 hct, which showed up-regulated significantly. At 1 hht (hours after hot treated), 22 of 26 GhCYPs are response to hot treatment. Of these, only the expression level of GhCYP-49 and GhCYP-61 was up-regulated significantly. At 3 hht, 6 hht and 12 hht, only 12 GhCYPs exhibited differential expression pattern.. Of these, up to 11 genes were up-regulated. Most GhCYP genes were found to be down-regulated under the condition of salt treatment apart from GhCYP-3, GhCYP-24, GhCYP-38, GhCYP-42, GhCYP-49, GhCYP-57 and GhCYP-61. After PEG (polyethylene glycol) treatment, only 8 GhCYPs displayed differential expression pattern, of which, GhCYP-3, GhCYP-42, GhCYP-49 and GhCYP-73 were up-regulated, GhCYP-41, GhCYP-47 and GhCYP-75 were down-regulated. Only the expression of GhCYP-38 showed be down-regulated at 3 hpt (hours after PEG treated) and 6 hpt, and then be up-regulated at 12 hpt. These expression patterns suggest that CYPs undertake multiple functions to help the cotton counter various complex environmental challenges.
Expression profiles of GhCYPs under the stress of V. dahliae
To gain a better understanding of the roles of CYP family genes in the cotton resistance against V. dahliae, the expression profiles were obtained by RNA-Seq. A total of 30 GhCYP genes showed differential expression at least one time point (hpi, hours post inoculation), including 19 upregulated genes and 11 downregulated genes (Fig. 4). Notably, levels of GhCYP-10, GhCYP-22, GhCYP-48 and GhCYP-59 were up-regulated at all five treatment time points. Additionally, the expression of GhCYP-11, GhCYP-27, GhCYP-37, GhCYP-45 and GhCYP-74 was down-regulated in most time points. These results revealed that GhCYPs were associated with the interaction between cotton and V. dahliae.
GhCYP-3 contains conserved amino acid residues and has PPIase activity
GhCYP-3, cloned from G. hirsutum cv. JM20, contains a single cyclophilin domain, 173 amino acid residues with a calculated molecular mass of 18.2 kDa and a pI of 8.34. Alignment with previously characterized Arabidopsis and human CYP (AtCYP19-1 and hCypA) revealed that GhCYP-3 contains seven conserved amino acid residues that critically affect PPIase activity (Fig. 5a). GhCYP-3 was found to be expressed in root, stem and leaf of unchallenged cotton plants with V. dalihae (Fig. 5b). Translation fusion of GhCYP-3 with GFP was constructed under the control of 35S promoter, and then transiently expressed in onion epidermal cells. Fluorescent imaging of the GhCYP-3-GFP bombarded onion epidermal cells showed both cytoplasm and nuclear localization (Fig. 5c).
The GhCYP-3 ORF was cloned into pET-32a at SacI and BglII sites, thus the recombinant plasmid was constructed (Fig. 6a). SDS-PAGE and western blot analysis revealed that the fusion TrxA-6×His-S-tag-6×His (THS, 20.4 kDa; empty vector as control) and TrxA-6×His-S-tag-GhCYP3 (THS-CYP, 34.6 kDa) proteins were highly expressed in the E. coli BL21(DE3) at 37°C with 1 mM IPTG for 4 h (Fig. 6b). The recombinant THS-CYP protein was purificated using the nickel-affinity. Further, we measured the PPIase activity of the recombinant protein by a coupled assay using synthetic peptide synthetic peptide succinyl- Ala-Ala-Pro-Phe-p-nitroanilide. The isomerization of the peptide substrate was observed in the presence of recombinant GhCYP-3, showing OD values higher than that of the blank, a negative control (Fig. 6c). These results indicated that GhCYP-3 has PPIase activity in vitro.
GhCYP-3 expression is upregulated in response to V. dahliae invasion
To analyze the expression pattern of GhCYP-3 in three cotton cultivars with different degrees of resistance to V. dahliae, the roots of two-week-old cotton seedlings were sampled at 0, 4, 8, 12, 24 and 48 hpi. The qRT-PCR results revealed that the transcriptional levels of GhCYP-3, compared with 0 hpi, were significantly upregulated in two-resistant cultivars Pima90-53 and JM-20 at 4 hpi. However, at the same inoculation time, GhCYP-3 in susceptible Han-208 was strongly downregulated and reached the highest expression level until 24 hpi (Fig 7). These results suggested that GhCYP-3 is upregulated and earlier involved in the cotton interaction with V. dahliae in resistant than in susceptible cultivars.
Overexpression of GhCYP-3 in Arabidopsis improves plant Verticillium wilt resistance
To further evaluate whether the GhCYP-3 functions in plant resistance to Verticillium wilt, we designed primers to amplify the ORF of GhCYP-3 into the plant expression vector pBI121 (Fig. 8a), which will make GhCYP-3 overexpression in transgenic Arabidopsis plants. By selection on kanamycin-containing medium (Fig. 8b), ten potential T1 transgenic lines were generated. After PCR-verification of gene insertion (Fig. 8c), two lines of the T3 generation displaying GhCYP-3 overexpression (Fig. 8d) were selected for analysis of PPIase activity and Verticillium wilt resistance.
Under unstressed conditions (treated with H2O), no significant PPIase activity was enhanced in two GhCYP3-expressing transgenic lines compared to the WT (wild type). However, these two transgenic lines, infected with V. dahliae for 3 d, showed increased PPIase activity of 0.51±0.03 (L1) and 0.45±0.03 nmol·s-1·mg-1 protein (L2) compared to 0.35±0.03 nmol·s-1·mg-1 protein observed in the WT (Fig. 8e).
In response to V. dahliae infection, two independent overexpressing lines displayed less wilting and smaller degree of leaf etiolation in comparison with WT plants at 15 dpi (Fig. 9a). Disease evaluation further indicated that the disease indices of transgenic lines were significantly lower than that of the WT (Fig. 9b). Furthermore, we examined the colonization in stems of infected WT and transgenic plants by isolation and cultivation of V. dahliae on PDA (potato dextrose agar) for 8 days. As a result, less fungal colonies came out from transgenic plants comparing to WT (Fig. 9c and d). These fungi recovery assay indicated that GhCYP-3 has a specialized effect of growth inhibition on V. dahliae.
GhCYP-3 exhibited obvious inhibitory effects on V. dahliae
We further assessed the inhibitory effects of GhCYP-3 on V. dahliae. As showed in Fig. 10a, distinct inhibition zones formed around the disc containing recombinant GhCYP-3. Additionally, V. dahliae spores were inoculated with plant extracts from WT and transgenic lines and then be spread on PDA plates. After 48 h, extracts from all the plants significantly reduced the number of colonies compared to the H2O (as control). However, the number of fungal colonies on plates contain transgenes extracts was significantly less than the WT (Fig. 10b and c). These results indicated that GhCYP-3 can efficiently inhibit conidia germinating and hyphae growth of V. dahliae in vitro.
CYP genes family has been systematically analyzed in several plants, such as Arabidopsis , rice  and soybean . Recently, 78 CYP genes were also identified in G. hirsutum L. acc. TM-1 from Illumina paired-end genomic sequencing (NAU version 1.1) by Chen et al. . Now, NAU version 2.1, highly accurate reference grade genome assemblies and annotations for G. hirsutum, was generated . A fairly large number of gaps and erroneous assemblies were successfully filled and corrected in new NAU version 2.1 . Therefore, in the present study we identified a total of 79 CYP genes by integrating NAU version 1.1, JGI database and NAU version 2.1 (Table 1). Compared to the previous report from Chen et al., the supplementary three CYPs (GhCYP-49, GhCYP-60 and GhCYP-70) were identified in this study. Besides, two CYPs (Gh_Sca140771G01 and Gh_A12G1281), undiscoverable in JGI database and NAU version 2.1, were removed. It thus made the candidate CYP genes in G. hirsutum more reliable.
All 79 CYPs in G. hirsutum have conservative PPIase domain (CLD). Of which, 14 GhCYPs are multi-domain proteins (Fig. 1). In addition to the CLD, GhCYP-14 and GhCYP-51 contain RRM and zf-CCHC. The RRM domain was found in proteins involved in RNA processing where it mediates binding to various RNAs to execute both housekeeping functions and regulatory mechanisms . Proteins that contain zinc fingers typically interact with DNA and RNA, and serve primarily to alter the binding specificity of a particular protein. In addition, GhCYP-14 and GhCYP-51 were predicted to be localized in the nuclear (Table 1). Therefore, these two CYPs possibly mediate ribosomal association while the CLD catalyzes peptidyl prolyl cis-trans isomerization of nascent polypeptides. GhCYP-19 and GhCYP-56 contain WD40 domain that generally serve as a rigid scaffold for protein interactions . Other 10 multi-domain GhCYPs contain TPR domain, which mainly act as interactive scaffolds in the formation of protein complexes and regulators of RNA metabolism involved in the immune response [26, 27].
Regulation of gene expression via specific cis-regulatory elements in the promoter regions has evolved as a major adaptive mechanism to respond to environmental stress in plants . Phytohormones are critical to the regulation of plant development and defense . Thus, the analysis of the putative cis-regulatory elements relating to hormone helps to advance our understanding of GhCYPs involving stress tolerance in cotton. Six hormones (ET, MeJA, ABA, GA, SA and Auxin) responsive regulatory elements were detected in the potential promoter regions of GhCYPs (Fig. 2), indicating that GhCYPs involve different hormone-mediated signaling pathways.
In 2010, we isolated 203 ESTs from a cDNA library using suppression subtractive hybridization (SSH) with a resistant upland cotton cultivar Jimian20 induced with V. dahliae. Of which, an EST encoded a partial polypeptide with homology to CYP . This is the first report that CYP is involved in the interaction between cotton and phytopathogen. In 2011, GhCyp1 was cloned from G. hirsutum cv. Zhongmian 35. Overexpression of GhCyp1 in transgenic tobacco plants conferred higher tolerance to salt stress and P. syringae pv. tabaci infection compared with control plants . Here the RNA-seq expression analysis of the cyclophilin gene families for G. hirsutum revealed that many of the genes potentially play important roles in various stress response (Figs. 3 and 4). So far as we know, almost all of the reported CYP genes were involved in plant stress reactions by their up-regulation expression, such as OsCYP19-4 and OsCYP21-4 from rice against cold  and salt , CcCYP from pigeon pea against drought, salinity and extreme temperatures , CyPs from Arabidopsis against wounding , and several cyclophilin genes in response to phytopathogen infection [15,16,17]. Therefore, significantly higher expression of CYP genes in cotton, such as GhCYP-49 and GhCYP-52 to cold, GhCYP-3, GhCYP-24, GhCYP-42, GhCYP-49 GhCYP-57 and GhCYP-61 to salt, GhCYP-3, GhCYP-24, GhCYP-27, GhCYP-38, GhCYP-42, GhCYP-49 and GhCYP-61 to hot, GhCYP-3, GhCYP-42, GhCYP-49 and GhCYP-73 to PEG, and GhCYP-3, GhCYP-5, GhCYP-10, GhCYP-17, GhCYP-18, GhCYP-20, GhCYP-22, GhCYP-23, GhCYP-28, GhCYP-30, GhCYP-33, GhCYP-46, GhCYP-48, GhCYP-51, GhCYP-54, GhCYP-59, GhCYP-62, GhCYP-64 and GhCYP-67 to V. dahliae suggests likely functional importance under these stress conditions. However, further functional studies are required to unravel the precise role of these candidates during cotton response to biotic and abiotic stress conditions.
GhCYP-3 contained the above-mentioned EST sequence related to V. dahliae-infection and showed 98% similarity in amino acid sequences with GhCyp1. Thus, GhCYP-3 potentially plays an important role in regulating cotton immune response. The results of expression (Fig. 7) and overexpression analysis (Fig. 9) further supposed that GhCYP-3 was involved in the cotton defence to V. dahliae. Plant CYPs could locate in multiple cell organelles, such as ER , chloroplast , Golgi , cytoplasm and nucleus [32, 33]. Interestingly, the GhCYP3-GFP fusion protein was localized in the cell cytoplasm as well as the cell nucleus in onion epidermal cells (Fig. 5c). Nuclear localization was proposed to play a possible role in the regulation of gene expression . GmCYP1 was as a “helper” that activates the enzymatic activity of a Phytophthora sojae RXLR effector Avr3b in a PPIase activity-dependent manner . Furthermore, GmCYP1 was demonstrated to interact with the isoflavonoid regulators GmMYB176 and 14-3-3 protein, suggesting that it participates in isoflavonoids metabolism and plays role in defense . AtCYP57 was proved to be involved in the A. thaliana response to P. syringae infection by influencing callose accumulation and PAD4 (peptidyl arginine deiminase type 4) expression, which interacts with EDS1 (enhanced disease susceptibility 1) to provide basal immune response in plants. Nucleus location makes authors inferred that AtCYP57 could directly regulate the translation of defence genes . Therefore, we inferred that GhCYP-3 play resistance function, may like AtCYP57 and GmCYP1, in the cell nucleus by directly interacting with some transcription factor to regulate the translation of defence genes, which lead to the production of antimicrobial metabolites.
Alternatively, GhCYP-3 was also located in the cytoplasm, indicating that it needs to play some extra roles, presumably including direct antifungal activity. Antifungal activities of CYPs have been reported from some plants including ginseng , Chinese cabbage , chickpea  and black-eyed pea . In our study, recombinant GhCYP-3 displayed evident inhibitory effects on V. dahliae on the plate (Fig. 10a). In a state of nature, the concentration of CYP in plant is hardly so high in the plate. Otherwise, there was no definite evidence to suggest that CYP is involved in plant resistance to pathogens by direct antifungal activity in vivo at present. However, we do not rule it out, because the extracts from GhCYP-3 transgenic Arabidopsis displayed significantly inhibit activity to conidia germinating and hyphae growth of V. dahliae (Fig. 10b). Many characterized antifungal proteins active on the fungal cell wall, plasma membrane and or intracellular targets . For example, Buforin 2 from the stomach tissue of Bufo bufo gargarizans can translocate through the plasma membrane and exhibits antifungal activity upon interaction with fungal DNA and RNA . During the infection process, V. dahliae forms penetration peg and specialized fungus-host interface to secret secretory effector proteins . This is probably a chance for antifungal proteins including GhCYP-3 produced by host plant encountering and entering into the cytoplasm of V. dahliae. Additionally, we identified 10 putative CYPs (VdCYPs) in V. dahliae strain VdLs.17 genome with very high sequence similarity with GhCYP-3 (Additional file 3: Table S3). These mean that GhCYP-3 may replace VdCYPs and carry out the same biological function in the cytoplasm of V. dahliae. Thus, we inferred that the antifungal activity of GhCYP-3 has been shown to be due to the effect on the normal function of VdCYPs, which is essential for the development of V. dahliae. Nevertheless, further research is needed to confirm this bold deduction.
This is the first systematic analysis of CYP family genes in cotton aiming to help clarifying the gene sequence characteristics and expression patterns. The putative cis-regulatory elements predication and expression divergence suggested that the GhCYP genes are involved in multiple phytohormone regulation pathways and responses to various abiotic stress and V. dahliae infection. These results will provide potential clues for the selection of candidate genes for further in-depth study on the functional characterization. Furthermore, GhCYP-3 showed both cytoplasmic and nuclear localization. Heterologous overexpression of GhCYP-3 in Arabidopsis significantly improved Verticillium wilt resistance of the plants. Recombinant GhCYP-3 and the extracts from GhCYP-3 transgenic Arabidopsis displayed significantly inhibit activity to V. dahliae. These results indicated that GhCYP-3 was associated with the resistance of cotton to V. dahliae infection presumably through antifungal activity, and it will offer an important candidate gene for Verticillium wilt tolerance in cotton molecular breeding.
Identification and characterization of CYP family genes in G. hirsutum
Genome assemblies of G. hirsutum TM-1 from Nanjing Agricultural University (NAU version 1.1 and version 2.1) and JGI (version 1.0) were retrieved from the CottonFGD website (https://cottonfgd.org/). Arabidopsis CYPs, accessed from TAIR website (https://www.arabidopsis.org/) were used as query to identify putative CYPs in G. hirsutum genomes by local BLAST using BioEdit software (Ibis Biosciences, Carlsbad, CA, USA). The identified GhCYPs were further verified through the Pfam database (http://pfam.xfam.org). The molecular weight (MW) and isoelectric point (pI) of each protein were calculated using ExPASy program (http://www.expasy.org/). The signal peptide was predicted with the program of SignalP 4.1. The amino acid sequences were aligned with DNAMAN software (Vers. 7; Lynnon Corporation, Quebec, Canada), using default parameters. CELLO v2.5 (http://cello.life.nctu.edu.tw/) was used for the subcellular localization prediction of GhCYPs. The putative cis-acting elements in the promoter regions were predicted using NAU version 2.1 database with the Plant CARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).
Genome-wide expression analysis of GhCYPs
A genome-wide expression analysis of the cotton CYP genes in various abiotic stresses and V. dahliae infection was performed using high-through RNA sequence data, which was downloaded from NCBI databases (SRP044705) and extracted from our RNA-seq data . Genes with FPKM ≥ 10 were used for further expression analysis.
Plant materials and V. dahliae strains
The cotton seeds G. hirsutum cv. Ji Mian 20 (JM20), Han208, CCRI8 and G. barbadense cv. Pima90-53 were preserved at the North China Key Laboratory for Crop Germplasm Resources of Education Ministry, Hebei Agricultural University, Baoding, China. Cotton seedlings were grown in commercial sterilized soil at 28°C /25°C (day/night) temperatures with a 16-h-light/8-h-dark regime. A. thaliana was grown in pots containing vermiculite soil with temperature at 23°C day and 20°C night, under a 16/8 h photoperiod. V. dahliae strain Linxi2-1 was isolated from a symptomatic upland cotton plants growing in agricultural fields in Linxi county, Hebei Province, China . These highly aggressive defoliating V. dahliae strains were maintained on PDA. The conidial suspension was prepared according to previous description  and adjusted to 107 spores per milliliter and 106 spores per milliliter with distilled water was used to the inoculation of cotton and Arabidopsis, respectively. The plant was infected with V. dahliae using soil drench method . 10 mL of the conidial suspension was directly injected with a needle without piercing into the bottom of each pot. Seedlings received sterile water in the same manner were used as control.
Gene cloning and subcellular localization
Total RNA was extracted from leaf tissues of JM20 with an RNA plant plus reagent (TIANGEN Biotech, China). First-strand cDNA was synthesized from an aliquot of 1 μg of total RNA with a PrimeScript™ RT Reagent Kit and gDNA Eraser (TaKaRa, China). GhCYP-3 was amplified with primers CYP-F1 and CYP-R1 (Additional file 1: Table S1), designed based on the sequences of Gh_A01G1361 (G. hirsutum L. acc. TM-1) . A GhCYP-3-GFP fusion construct under the control of the 35S promoter was generated by cloning the ORF into the SalI and BamHI sites of the binary vector pCamE. The vector expressing GFP alone served as control. Protein subcellular localization in onion (Allium cepa) epidermal cells was determined according to the protocol of Yang (2015) .
qRT-PCR and semi-quantitative PCR
Total RNA and cDNA were prepared by the method described above. The qRT-PCR was performed on a CFX96 Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) using the qPCR kit for SYBR Green (TaKaRa, China). qPCR conditions consisted of one cycle of 3 min at 95 °C, followed by 40 cycles of 15 s at 95 °C, 15 s at 58 °C, and 20 s at 72 °C. The cotton ubiquitin 14 (UBQ14) gene served as an internal standard . Fold-changes in expression were calculated via the 2−ΔCt method. Semi-quantitative RT-PCR was done with an Applied Biosystems® 2720 Thermal Cycler with Arabidopsis AtACT2 (AT3G18780) as the internal standard . The resulting products were resolved on a 1.5% agarose gel. All primers are listed in Additional file 1: Table S1 Three biological and three technical replicates were analyzed for all quantitative experiments.
Generation and evaluation resistance of transgenic Arabidopsis
GhCYP-3 was cloned into pBI121 vector at the XbaI and SacI sites by PCR with primers CYP-X and CYP-S (Additional file 1: Table S1). The chimeric construct was introduced into Agrobacterium tumefaciens strain GV3101 for Arabidopsis transformation using the floral dip method . Putative transformants were selected on MS (Murashige and Skoog) medium containing 50 mg L-1 kanamycin, and then be further verified by PCR for gene insertion and semi-quantitative RT-PCR for gene expression. Independent T1 transgenic lines were used to produce the T3 generations, which were randomly chosen as representative lines and subjected to analysis. Disease severity of Arabidopsis plants caused by V. dahliae was assessed according to symptoms manifested on the leaves. The disease index (DI) was calculated as previously described . Fungal recovery assay was also used to evaluate the resistance of the plant through assessment of Verticillium colonization recovered from stem sections according to the method of Fradin (2009) . In each treatment, 18 individual plants were used and all the experiments were repeated thrice.
Purification and PPIase activity assay of recombinant protein
The GhCYP-3 ORF was cloned into expression vector pET-32a (+) (Novagen, Darmstadt, Germany) with forward primer CYP-Bg-F and reverse primer CYP-Sa-R (Additional file 1: Table S1), which will introduce the BglII and SacI site into the 5′ and 3′ end of the ORF, respectively. To induce the expression of GhCYP-3 protein with a His-tag in E. coli BL21(DE3) (TransGen Biotech, Beijing, China), a final concentration of 1.0 mmol·L−1 isopropyl-β-D-thiogalactopyranoside (IPTG) was added into the culture when the OD600 value reached 0.4–0.6, and the culture was allowed to continue growing for 4 – 6 h before harvesting. The proteins were separated on a SDS–PAGE gel and detected by western blot using anti His-Tag mouse monoclonal antibody (1:5000; CW Biotech, Beijing, China). The recombinant protein was purified using a 6×His-Tagged Protein Purification Kit (CW Biotech, Beijing, China). The PPIase activity of the recombinant protein was assayed in vitro using the tetrapeptide substrate Suc-AAPF-pNA (N-succinyl-Ala-Leu-Pro-Phe-p-nitroanilide; Sigma-Aldrich, Ontario, Canada) in a Shimadzu UV-2450 spectrophotometer (Shimadzu, Kyoto, Japan) as described by Yoon (2016) . Three biological and two technical replicates were analyzed for all measurement.
Assay of antifungal activity
The antifungal activity of the recombinant protein was tested against V. dahliae using filter paper disc diffusion method. V. dahliae spores were uniformly spread on the PDA medium plates and then cultured for 60 h at 28°C. Sterilized blank paper discs of 6 mm diameters, impregnated with tested protein, were placed on the surface of PDA medium previously spread with V. dahliae. The plates were inoculated at 25°C for 10 d. The antifungal activity of extracts from Arabidopsis plants transformed with GhCYP-3 was performed as described previously  with the following modifications. Conidial suspension adjusted to a density of 105 conidia ml-1, and germinated on PDA overnight at 25 °C prior to assay. Total homogenates (5 g) from ten Arabidopsis plants were prepared by directly grinding plant leaves into a fine powder in liquid nitrogen with no buffer added. Subsequently, extracts were collected by centrifugation at 10000 g for 10 min at 25 °C. Conidial suspensions (25 μl) were mixed with 225 μl of extract, and incubated for 1 h at 25°C. The mixture (50 μl) were spread onto PDA plates and incubated at 25 °C for 48 h and fungal colonies enumerated. All experiments were repeated three times independently.
All experiments were performed at least three times for each determination. Statistical analysis was performed using GraphPad Prism® 6 software (Graph Pad, San Diego, CA, USA). Unless otherwise indicated, data were evaluated by using analysis of variance (ANOVA) followed by Dunnett's multiple comparisons test.
Availability of data and materials
The data generated or analyzed during the current study are included in this published article and its supplemental data files and available from the corresponding author on reasonable request.
Abscisic acid responsive element
Cyclophilin type PPIase domain
Enhanced disease susceptibility 1
Ethylene responsive element
Fragments per kilobase of exon per million fragments mapped
Hours after cold treated
Hours after hot treated
Hours post inoculation
Hours after PEG treated
Murashige and Skoog
- PAD4 :
Peptidyl arginine deiminase type 4
Potato dextrose agar
Peptidyl prolyl cis/trans isomerase
Reactive oxygen species
RNA recognition motif
Chen ZJ, Scheffler BE, Dennis E, Triplett BA, Zhang T, Guo W, Chen X, Stelly DM, Rabinowicz PD, Town CD, et al. Toward sequencing cotton (Gossypium) genomes. Plant Physiol. 2007;145(4):1303–10.
Zhang J, Flynn R, Baral JB, Bajaj S, Hughs SE, Percy RGJE. Germplasm evaluation and transfer of Verticillium wilt resistance from Pima (Gossypium barbadense) to Upland cotton (G. hirsutum). Euphytica. 2012;187(2):147–60.
Romano P, Gray J, Horton P, Luan S. Plant immunophilins: functional versatility beyond protein maturation. New Phytol. 2005;166(3):753–69.
Olga B, Marek Z, Yamile M, Tatsiana S, Maria K, Andrea BJNAR. Identification of RNA targets for the nuclear multidomain cyclophilin atCyp59 and their effect on PPIase activity. Nucleic Acids Res. 2013;41(3):1783–96.
Kumari S, Roy S, Singh P, Singla-Pareek SL, Pareek A. Cyclophilins: proteins in search of function. Plant Signaling & Behavior. 2013;8(1):e22734.
Kristine NB, Mallis RJ, Fulton DB, Andreotti AH. Regulation of the tyrosine kinase Itk by the peptidyl-prolyl isomerase cyclophilin A. Proc Natl Acad Sci U S A. 2002;99(4):1899–904.
Hong F, Lee J, Song JW, Lee SJ, Ahn H, Cho JJ, Ha J, Kim SS. Cyclosporin A blocks muscle differentiation by inducing oxidative stress and inhibiting the peptidyl-prolyl-cis-trans isomerase activity of cyclophilin A: cyclophilin A protects myoblasts from cyclosporin A-induced cytotoxicity. FASEB J. 2002;16(16):1633–5.
Romano PG, Horton P, Gray JE. The Arabidopsis cyclophilin gene family. Plant Physiol. 2004;134(4):1268–82.
Juncheul A, Daewon K, You YN, Minsook S, Jeongmee P, Hyunsik H, Beomgi K, Sheng L, Hongseog P, Hyesun C. Classification of rice (Oryza sativa L. japonica Nipponbare) immunophilins (FKBPs, CYPs) and expression patterns under water stress. BMC Plant Biol. 2010;10(1):1–22.
Mainali HR, Chapman P, Dhaubhadel S. Genome-wide analysis of cyclophilin gene family in soybean ( Glycine max ). BMC Plant Biol. 2014;14(1):1–11.
Chou IT, Gasser CS. Characterization of the cyclophilin gene family of Arabidopsis thaliana and phylogenetic analysis of known cyclophilin proteins. Plant Mol Biol. 1997;35(6):873.
Yoon DH, Sang SL, Park HJ, Lyu JI, Chong WS, Liu JR, Kim BG, Ahn JC, Cho HS. Overexpression of OsCYP19-4 increases tolerance to cold stress and enhances grain yield in rice (Oryza sativa). J Exp Bot. 2016;67(1):69–82.
Lee SS, Park HJ, Jung WY, Lee A, Yoon DH, You YN, Kim HS, Kim BG, Ahn JC, Cho HS. OsCYP21-4, a novel Golgi-resident cyclophilin, increases oxidative stress tolerance in rice. Front Plant Sci. 2015;6:797.
Sekhar K, Priyanka B, Reddy VD, Rao KV. Isolation and characterization of a pigeonpea cyclophilin (CcCYP) gene, and its over-expression in Arabidopsis confers multiple abiotic stress tolerance. Plant Cell Environ. 2010;33(8):1324–38.
Kong HY, Lee SC, Hwang BK. Expression of pepper cyclophilin gene is differentially regulated during the pathogen infection and abiotic stress conditions. Physiol Mol Plant Pathol. 2001;59(4):189–99.
Godoy AV, Lazzaro AS, Casalongué CA, San Segundo B. Expression of a Solanum tuberosum cyclophilin gene is regulated by fungal infection and abiotic stress conditions. Plant Sci. 2000;152(2):123–34.
Figueiredo A, Monteiro F, Pais MS, Rex M, Töpfer R, Zyprian E. Cyclophilin: A downy mildew resistance gene candidate in grapevine. Acta Horticulturae. 2014;1046(1046):371–8.
Zhu C, Wang Y, Li Y, Bhatti KH, Tian Y, Wu J. Overexpression of a cotton cyclophilin gene (GhCyp1) in transgenic tobacco plants confers dual tolerance to salt stress and Pseudomonas syringae pv. tabaci infection. Plant Physiol Biochem. 2011;49(11):1264–71.
Li FG, Fan GY, Lu CR, Xiao GH, Zou CS, Kohel RJ, Ma ZY, Shang HH, Ma XF, Wu JY, et al. Genome sequence of cultivated Upland cotton (Gossypium hirsutum TM-1) provides insights into genome evolution. Nat Biotechnol. 2015;33(5):524–30.
Zhang TZ, Hu Y, Jiang WK, Fang L, Guan XY, Chen JD, Zhang JB, Saski CA, Scheffler BE, Stelly DM, et al. Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat Biotechnol. 2015;33(5):531–7.
Hu Y, Chen J, Fang L, Zhang Z, Ma W, Niu Y, Ju L, Deng J, Zhao T, Lian J, et al. Gossypium barbadense and Gossypium hirsutum genomes provide insights into the origin and evolution of allotetraploid cotton. Nat Genet. 2019;51(4):739–48.
Trivedi DK, Yadav S, Vaid N, Tuteja N. Genome wide analysis of cyclophilin gene family from rice and Arabidopsis and its comparison with yeast. Plant Signal Behav. 2012;7(12):1653.
Chen Q, Chen Q-J, Sun G-Q, Zheng K, Yao Z-P, Han Y-H, Wang L-P, Duan Y-J, Yu D-Q, Qu Y-Y. Genome-wide identification of cyclophilin gene family in cotton and expression analysis of the fibre development in Gossypium barbadense. Int J Mol Sci. 2019;20(2):349.
Albà MM, Pagès M. Plant proteins containing the RNA-recognition motif. Trends Plant Sci. 1998;3(97):15–21.
Smith TF, Gaitatzes C, Saxena K, Neer EJ. The WD repeat: a common architecture for diverse functions. Trends Biochem Sci. 1999;24(5):181–5.
Sane AP, Stein B, Westhoff P. The nuclear gene HCF107 encodes a membrane-associated R-TPR (RNA tetratricopeptide repeat)-containing protein involved in expression of the plastidial psbH gene in Arabidopsis. Plant J. 2005;42(5):720–30.
Chen L, Hamada S, Fujiwara M, Zhu T, Thao NP, Wong HL, Krishna P, Ueda T, Kaku H, Shibuya N, et al. The Hop/Sti1-Hsp90 chaperone complex facilitates the maturation and transport of a PAMP receptor in rice innate immunity. Cell Host Microbe. 2010;7(3):185–96.
Dirk W, Roman B, Joachim S. The regulatory code for transcriptional response diversity and its relation to genome structural properties in A. thaliana. Plos Genetics. 2007;3(2):e11.
Wang C, Liu Y, Li S-S, Han G-Z. Insights into the origin and evolution of plant hormone signaling machinery. Plant Physiol. 2015;167(3):872–86.
Zhang CY, Wang XF, Zhang GY, Lq W, Chi JN, Li ZK, Ma ZY. ESTs analysis of suppression subtractive hybridization library from upland cotton resistant cultivar infected by Verticillium dahliae. Cotton Sci. 2010;22(01):17–22.
Saito T, Niwa Y, Ashida H, Tanaka K, Kawamukai M, Matsuda H, Nakagawa T. Expression of a gene for cyclophilin which contains an amino-terminal endoplasmic reticulum-targeting signal. Plant Cell Physiol. 1999;40(1):77–87.
Pogorelko GV, Mokryakova M, Fursova OV, Abdeeva I, Piruzian ES, Bruskin SA. Characterization of three Arabidopsis thaliana immunophilin genes involved in the plant defense response against Pseudomonas syringae. Gene. 2014;538(1):12–22.
Mainali HR, Vadivel AK, Li X, Gijzen M, Dhaubhadel S. Soybean cyclophilin GmCYP1 interacts with an isoflavonoid regulator GmMYB176. Sci Rep. 2017;7:39550.
Kong G, Zhao Y, Jing M, Huang J, Yang J, Xia Y, Kong L, Ye W, Xiong Q, Qiao Y. The activation of phytophthora effector Avr3b by plant cyclophilin is required for the nudix hydrolase activity of Avr3b. PLoS Pathog. 2015;11(8):e1005139.
Zhang H, Wang J, Li S, Wang S, Liu M, Wang W, Zhao Y. Molecular cloning, expression, purification and functional characterization of an antifungal cyclophilin protein from Panax ginseng. Biomedical Reports. 2017;7(6):527–31.
Lee JR, Park S-C, Kim J-Y, Lee SS, Park Y, Cheong G-W, Hahm K-S, Lee SY. Molecular and functional characterization of a cyclophilin with antifungal activity from Chinese cabbage. Biochem Biophys Res Commun. 2007;353(3):672–8.
Ye XY, Ng TB. Isolation of a new cyclophilin-like protein from chickpeas with mitogenic, antifungal and anti-HIV-1 reverse transcriptase activities. Life Sci. 2002;70(10):1129–38.
Ye XY, Ng TB. Isolation of Unguilin, a cyclophilin-like protein with anti-mitogenic, antiviral, and antifungal activities, from Black-Eyed Pea. J Protein Chem. 2001;20(5):353–9.
Theis T, Stahl U. Antifungal proteins: targets, mechanisms and prospective applications. Cell Mol Life Sci. 2004;61(4):437–55.
Chan BP, Kim HS, Sun CK. Mechanism of action of the antimicrobial peptide Buforin II: Buforin II kills microorganisms by penetrating the cell membrane and inhibiting cellular functions. Biochem Biophys Res Commun. 1998;244(1):253–7.
Zhou TT, Zhao YL, Guo HS. Secretory proteins are delivered to the septin-organized penetration interface during root infection by Verticillium dahliae. PLoS Pathog. 2017;13(3):e1006275.
Zhang Y, Wang X, Rong W, Yang J, Li Z, Wu L, Zhang G, Ma Z. Histochemical analyses reveal that stronger intrinsic defenses in Gossypium barbadense than in G. hirsutum are associated with resistance to Verticillium dahliae. Mol Plant Microbe Interact. 2017;30(12):984–96.
Wang GN, Zhao GY, Yue XW, Li ZK, Zhang Y, Zhang GY, Wu LQ, Wang XF, Ma ZY. Pathogenicity and ISSR genetic differentiation of Verticillium dahliae isolates from cotton growing areas of Hebei Province. Cotton Science. 2012;24(4):348–57.
Yang J, Ji L, Wang X, Zhang Y, Wu L, Yang Y, Ma Z. Overexpression of 3-deoxy-7-phosphoheptulonate synthase gene from Gossypium hirsutum enhances Arabidopsis resistance to Verticillium wilt. Plant Cell Rep. 2015;34(8):1429–41.
Artico S, Nardeli SM, Brilhante O, Grossi-de-Sa MF, Alves-Ferreira M. Identification and evaluation of new reference genes in Gossypium hirsutum for accurate normalization of real-time quantitative RT-PCR data. BMC Plant Biol. 2010;10:49.
Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible W-R. Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005;139(1):5–17.
Clough SJ, Bent AF. Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. Plant J. 1998;16(6):735–43.
Zhang B, Yang Y, Chen T, Yu W, Liu T, Li H, Fan X, Ren Y, Shen D, Liu L, et al. Island cotton Gbve1 gene encoding a receptor-like protein confers resistance to both defoliating and non-defoliating isolates of Verticillium dahliae. PLoS One. 2012;7(12):e51091.
Fradin EF, Zhang Z, Juarez Ayala JC, Castroverde CD, Nazar RN, Robb J, Liu CM, Thomma BP. Genetic dissection of Verticillium wilt resistance mediated by tomato Ve1. Plant Physiol. 2009;150(1):320–32.
Cary JW, Jaynes JM, Cleveland TRK. Transgenic expression of a gene encoding a synthetic antimicrobial peptide results in inhibition of fungal growth in vitro and in planta. Plant Sci. 2000;154(2):171–81.
All authors are grateful to the laboratory members for help, advice and discussion.
This work was supported by grants from the National Key Research and Development Program of China (2016YFD0101006), the Top Talent Fund of Hebei Province and the Fund of Research Group Construction for Crop Science in Hebei Agricultural University (TD2016C328). None of the funding bodies were involved with the research other than providing funds to support the work.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Yang, J., Wang, G., Ke, H. et al. Genome-wide identification of cyclophilin genes in Gossypium hirsutum and functional characterization of a CYP with antifungal activity against Verticillium dahliae. BMC Plant Biol 19, 272 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12870-019-1848-1
- Verticillium dahliae
- Antifungal activity