Abstract
Cyclophilins (CYPs) are a group of ubiquitous proteins characterized by their ability to bind to the immunosuppressive drug cyclosporin A. The CYP family occurs in a wide range of organisms and contains a conserved peptidyl-prolyl cis/trans isomerase domain. In addition to fulfilling a basic role in protein folding, CYPs may also play diverse important roles, e.g. in protein degradation, mRNA processing, development, and stress responses. We performed a genome-wide database survey and identified a total of 94 CYP genes encoding 91 distinct proteins. Sequence alignment analysis of the putative BnCYP cyclophilin-like domains revealed highly conserved motifs. By using RNA-Seq, we could verify the presence of 77 BnCYP genes under control conditions. To identify phloem-specific BnCYP proteins in a complementary approach, we used LC-MS/MS to determine protein abundances in leaf and phloem extracts. We detected 26 BnCYPs in total with 12 being unique to phloem sap. Our analysis provides the basis for future studies concentrating on the functional characterization of individual members of this gene family in a plant of dual importance: as a crop and a model system for polyploidization and long-distance signalling.
Similar content being viewed by others
Introduction
As one of the most important crops for nutritional oil, fodder, biodiesel, chemical and pharmaceutical products, Brassica napus (oilseed rape or canola) is widespread in agriculture, especially in the European Union, China, and Canada. Beside its essential agricultural significance, B. napus has also become a model plant for studying long-distance signalling1,2,3,4,5 and evolutionary consequences of polyploidization6, 7. Brassica napus is used as a model plant for studying long-distance signalling, because methods for the collection of sufficient quantities of pure xylem and phloem sap are well established2, 4. Compared to crops like wheat, soybean or rice, the domestication of B. napus was more recent. It is assumed that chromosome doubling occurred after spontaneous hybridization between Brassica rapa (Asian cabbage or turnip rape, 2n = 2 × 10 = 20, genome AA) and Brassica oleracea (Mediterranean cabbage, 2n = 2 × 9 = 18, genome CC)8. The profitable cultivation of B. napus (genome AACC, 2n = 38) has high importance for the economy9.
A group of proteins involved in diverse fundamental cellular functions in many different organisms is called immunophilins, originally discovered as receptors for immunosuppressive drugs in mammals. The family of immunophilins consists of two major groups, FK506-binding proteins (FKBPs)10, 11 and cyclophilins (CYPs)12, 13. Despite the lack of sequence and structure similarity, FKBPs and CYPs each possess a conserved domain responsible for their common peptidyl-prolyl cis/trans isomerase (PPIase) activity, catalyzing the rate-limiting rotation of X-proline peptide bonds from a cis to a trans conformation14. These domains are called the FKBP and cyclophilin-like domain (CLD), respectively. An additional group of proteins exhibiting PPIase activity, the parvulins, cannot be classified as immunophilins in a strict sense, since they do not bind to any known immunosuppressant molecule15. The drugs of CYPs and FKBPs are cyclosporin A (CsA) and FK506/rapamycin, respectively, that bind to the catalytic pocket of the PPIase domain16, thereby inhibiting its activity and forming a high-affinity binding site for the interaction with calcineurin17. However, these drugs do not occur naturally in cells, therefore the consequences of drug treatment have clinical but no physiological relevance18.
Human CYPs were first characterized by their ability to bind the drug cyclosporin A in 198412. The first plant CYPs were described in 1990, where CYP cDNA sequences were identified in tomato (Lycopersicon esculentum), maize (Zea mays), and oilseed rape (B. napus)19. The ubiquitous CYPs show highly conserved structural features20 and are involved in several fundamental cellular functions, including protein folding21,22,23, protein trafficking24, signalling25,26,27,28,29, pathogen response30, apoptosis31, RNA-binding, regulation of gene expression or transcription32,33,34,35,36, and plant stress responses37,38,39,40.
Interestingly, CYPs were found to be a prominent and abundant class of proteins in the phloem long-distance system of higher plants including castor bean and oilseed rape2, 41. Phloem sieve elements lose their ability for transcription and translation during their maturation into transport tubes42. However, phloem exudate contains a complex set of proteins, some of which have been implicated with long-distance signalling43. Such signalling proteins are imported from the adjacent companion cells through plasmodesmata. It was proposed that phloem CYPs might act as molecular chaperones in this process, potentially involved in refolding of non-cell-autonomous proteins after entry into the translocation stream2, 41.
The growing number of sequenced genomes allowed the identification of whole sets of CYPs in various organisms by sequence comparisons. 19 CYPs were detected in the human genome, whereas Saccharomyces cerevisiae possesses 8, Schizosaccharomyces pombe 9, Caenorhabditis elegans 17, Drosophila melanogaster 1444, 45, and the fungus Leptosphaeria maculans 1246 CYPs. Compared to other organisms, plants have a higher number of CYPs with 29 encoded in A. thaliana 47,48,49, and 27 in rice37. To date, soybean (Glycine max) is reported to have the largest set of CYPs with 62 members50.
Based on the recently sequenced B. napus genome51, the major aim of the present study was the identification and classification of CYP-like proteins in this economically important species. A total number of 94 genes belonging to the CYP gene family (resulting in 91 different CYP proteins) could be identified in the sequenced B. napus cultivar ‘Darmor-bzh’. By applying transcriptome analysis of the B. napus cultivar ‘Drakkar’ we could confirm the expression of 77 BnCYP genes in leaves. To identify CYP proteins specifically occurring in phloem sap, we performed complementary protein analyses by LC-MS/MS on leaf and phloem extracts and found 26 different BnCYP proteins in total, 12 only present in phloem sap.
Results and Discussion
Identification of cyclophilins in the Brassica napus genome
Putative CYPs of B. napus containing full length or partial CYP-like domains were identified by BLASTp of A. thaliana CYPs47, 48. The 94 B. napus (cultivar ‘Darmor-bzh’) gene sequences determined by this approach resulted in 91 distinct proteins which were subjected to additional analysis like sequence alignments with the respective Arabidopsis homologs and the verification of CLDs. The analysis showed that B. napus contains the largest CYP gene family known to date with 94 genes followed by soybean (Glycine max) with 62 members50.
As proposed by Romano et al.47 and He et al.48 we used the nomenclature BnCYP (Brassica napus cyclophilin) followed by the molecular weight and a consecutive number for genes encoding proteins with similar molecular weight. We used the molecular weight of the immature proteins as the basis, because the prediction of potential cleavable signal peptides by different prediction tools was not unambiguous. Table 1 summarizes the information about all identified BnCYPs.
The cyclophilin-like domain
Since CYPs are characterized by the highly conserved CLD, its occurrence in the identified potential CYP sequences was verified. Most of the BnCYPs contained full length CLDs, whereas CYPs with a molecular weight below 17 kDa have partial CLDs missing essential residues or whole secondary structure parts and thereby might be lacking PPIase activity.
As a consequence of the CLD analysis, two proteins annotated as CYPs in the first place, BnaC05g10020D (30 kDa) and BnaC04g45810D (54 kDa), were excluded from subsequent studies. In addition, for one of the putative low molecular weight CYPs, BnCYP7-1, no CLD was predicted. Nevertheless, it was retained, because a BLASTp search and sequence alignment suggested a partial CLD.
Figure 1 depicts the conserved sequence of all identified B. napus CLDs aligned with the secondary structure of human cyclophilin A (hCYPA) as a reference (for a detailed view of the multiple alignment see Supplementary Fig. S1). The crystal structure of hCYPA shows eight β-strands forming a β-barrel structure, plus two additional α-helices located at the top and bottom52. Human CYPA is often referred to as the “archetypal” CYP. Residues important for CYP function as well as residues promoting the secondary structure are conserved, whereas gaps in the conserved sequence represent insertions in individual members of this protein family. In previous studies it has been demonstrated that the highly conserved amino acid W121 of hCYPA (Fig. 1, W225) is required for CsA binding and interacts directly with the phosphatase calcineurin, but it is not essential for PPIase activity53, 54. This amino acid was not present in 57 out of the 91 BnCYP proteins. Further highly conserved amino acids fundamental for the PPIase activity of hCYPA are R55 (Fig. 1, R87), F60 (Fig. 1, F102;), and H126 (Fig. 1, H232)53. Three single mutations of these amino acids reduced the original activity of the human wild-type isomerase to less than 1%. However, the mutant proteins were still able to bind CsA. Since some of the BnCYPs did not possess all of these three amino acids, their PPIase activity needs to be experimentally confirmed. The motif VXGXV of the hydrophobic region around β-sheet-VII is reported to be conserved among all AtCYP proteins47. This motif was also present in 62 BnCYP protein sequences (Fig. 1, V234, G236, V238).
Exceptional insertions of amino acids in plant CYPs between β-sheet-I and β-sheet-II (β-I/β-II), β-sheet-IV and β-sheet-V (β-IV/β-V), and β-sheet-VI and α-helix-II (β-VI/α-II)20 have also been observed in this study. In more detail, 3 amino acids are inserted between β-I/β-II (Fig. 1, residues 25–27) in 16 BnCYPs. Yet, 45 BnCYPs do not show this insertion. Interestingly, we also observed deletions in 30 BnCYPs in this region. Only a minority of BnCYPs contains an insertion between β-IV/β-V (Fig. 1, residues 133–152) with 1 to 19 amino acids in length, but also here deletions occurred in some cases. Insertions between β-VI/α-II (Fig. 1, residues 216–221) are only present in 5 BnCYPs with medium and higher molecular masses. An additional insertion between α-helix-I and β-sheet-III (α-I/β-III junction) has been described by Romano et al.47. Here, 8 to 11 amino acids are inserted in several AtCYPs and 3 to 4 amino acids in chloroplast variants. For B. napus a 7 amino acid insertion occurs in several CYPs. For example, except for BnCYP18-1, all 18 and 19 kDa BnCYPs show this insertion.
Domain architecture of cyclophilins and homology
Depending on the domains present, CYPs are classified as single- (SD) and multi-domain (MD) forms48. SD CYPs contain only a CLD which was the case for most of the BnCYPs with a rather low molecular weight (18–30 kDa). While there are mostly two BnCYP homologs with high sequence identity (>80%) corresponding to each AtCYP, caused by polyploidization, there are only a few exceptions with only one corresponding BnCYP homolog. One example is BnCYP18-1 which is the only B. napus homolog of AtCYP18-2 with a full length CLD (with 95% sequence identity, Table 2). It has been shown that the Arabidopsis CYP AtCYP18-2 is recruited by AtSKIP to the nucleus to regulate pre-mRNA splicing55. Because of the high homology, the oilseed rape variant might fulfil a similar role in the nucleus of B. napus cells.
For 12 B. napus CYPs additional domains are predicted which promote diverse capabilities like interaction (protein-protein, protein-DNA, protein-RNA) and modification (ubiquitination) (Fig. 2). These seem to be conserved among Brassicaceae, since similar domain structures exist in the corresponding A. thaliana homologs.
The MD CYPs BnCYP40-1 and BnCYP40-2 are characterized by a tetratricopeptide repeat (TPR) motif with two TPRs, each 34 amino acids long. Such motifs mediate protein-protein interactions and can thereby assist in the assembly of multi-protein complexes. BnCYP40-1 and BnCYP40-2 show a high degree of sequence identity (92 and 93%, respectively) to AtCYP40 (Table 2). A cytoplasmic interaction partner of AtCYP40 is HSP90, which mediates its recruitment to an intermediate RISC complex by its TPR motifs56, 57. Thereby, AtCYP40 forms, together with HSP90, a complex with AGO1 and a small RNA duplex before the mature RISC complex consisting of AGO1 and a siRNA or miRNA strand is formed56. The mechanism of CYP40-HSP90 binding is conserved between different species. Besides plants, human HSP90 is also known to bind to hCYP40, the human AtCYP40 homolog58. Thus, a similar function for the highly identical CYPs BnCYP40-1 and BnCYP40-2 can be assumed.
Other CYPs showing classical protein-protein interaction domains are BnCYP70-1 and BnCYP70-2, both containing four tryptophan-aspartic acid (WD40) repeats. These are short structural motifs typically forming a four stranded anti-parallel β-sheet. Multiple copies build a circular β-propeller structure promoting protein-protein interactions. Both show 91% sequence identity to their homolog in Arabidopsis, AtCYP71 (Table 2), which is located in the nucleus and functions in the regulation of gene expression and organogenesis. The WD40 domain enables AtCYP71 to associate with histone H3 affecting its methylation35. Data from yeast PPIases provide evidence that they are responsible for histone modification by utilizing the peptidyl-prolyl bond isomerisation as a molecular switch for transcription regulation59. Furthermore, AtCYP71 was shown to interact with FAS1 (a subunit of Chromatin Assembly Factor-1) and LHP1 (a heterochromatin protein)60, suggesting that AtCYP71 is involved in histone modification and chromatin assembly.
BnCYP47-2 and BnCYP47-3 both contain a putative PsbQ-like domain and have only small differences in their N-terminal amino acid sequences. Their closest homolog is AtCYP38 with 86 and 84% sequence identity (Table 2), which also contains a PsbQ-like domain. This domain is typical for proteins localized in the chloroplast, but its function is mostly unknown. AtCYP38 is the only Arabidopsis CYP for which a crystal structure is available. This revealed the additional PsbQ-like domain61. AtCYP38 is experimentally assigned to the chloroplast thylakoid lumen, does not show any PPIase activity62, 63, and plays a critical role in the assembly and maintenance of photosystem II61, 64, 65. It is suggested to interact with the E-loop of chlorophyll protein47 (CP47), a component of the photosystem II (PSII) complex61, via its CLD. Furthermore it is proposed to be responsible for proper folding and insertion of D1 and CP43, both components of PSII65.
There are also two BnCYPs showing classical domains for protein modification. BnCYP65-1 and BnCYP65-2 both contain a Zinc finger U-box motif. Therefore, they might be involved in the ubiquitin degradation pathway.
Four BnCYPs contain sequences characteristic for RNA-interacting proteins. BnCYP52 and BnCYP55 both possess a coiled coil domain. This domain represents a structural motif with coiled α-helices. Proteins containing such a domain can, for example, be transcription factors involved in the regulation of gene expression. The respective A. thaliana homolog is AtCYP57, which shows, besides a predicted coiled coil domain, also an S/K-R/E-rich region47. Sequence comparisons of BnCYP52, BnCYP55 and AtCYP57 revealed a high sequence identity in the AtCYP57 S/K-R/E-rich region (data not shown). This region is supposed to mediate the interaction with ribonucleoproteins66.
BnCYP112 contains an RNA recognition motif (RRM), a Zinc finger motif (CCHC-type), a positively charged region (arginine, lysine, glutamate) and an actin-like ATPase domain. Besides AtCYP59 in A. thaliana, homologs to BnCYP112 can be found in different species like P. tetraurelia, S. pombe, C. elegans, D. melanogaster, and H. sapiens 32. This group of proteins is called CRIPs, for cyclophilin-RNA interacting proteins32. Characteristic for this group of proteins is an RRM in addition to the CLD. Nevertheless, CRIPs differ in their C-terminal region. AtCYP59 and BnCYP112 are the only members containing a Zinc finger motif. AtCYP59 interacts with the C-terminal domain of the largest subunit of RNA polymerase II and thereby influences its phosphorylation state. Furthermore, binding of an RNA transcript decreases PPIase activity of AtCYP59. This might modulate RNA polymerase II activity36. Therefore, it is suggested that AtCYP59 connects pre-mRNA processing and transcription36, 67. Moreover, BnCYP112 is the only homolog of all CRIPs showing an additional actin-like ATPase domain.
BnCYP146 is the largest CYP found in oilseed rape. It contains a transmembrane domain and a Fip1 motif. In yeast, the Fip1 protein ensures the polyadenylation of mRNAs by interacting with poly(A) polymerase. There is no A. thaliana homolog showing the same domain structure as BnCYP146, but its CLD is closely related to AtCYP21-4 (Table 2).
Subcellular localization and phylogenetic relationships
The phylogenetic analysis of the full length amino acid sequences revealed clusters of proteins with high sequence similarity (Fig. 3). Often, these groups possess the same additional domains, consistent subcellular localization, and the same A. thaliana CYP homolog.
As previously described for other organisms (e.g. A. thaliana homologs, see Table 2), B. napus CYPs are targeted to all intercellular organelles. At least 50 BnCYPs are predicted to be located in the cytosol. This result is comparable to the subcellular distribution of the A. thaliana CYPs where 14 of 29 are either predicted or experimentally proven to be located in the cytosol (Table 2). Furthermore, 14 BnCYPs are predicted to be located in chloroplasts, 7 in mitochondria, 13 in the nucleus, and 7 are secreted (summarized in Fig. 3, for a detailed prediction see Supplementary Table S1). In addition, some BnCYPs contain nuclear localization signals (NLS). Proteins predicted to be cytosolic or nuclear which possess NLS might relocate between both compartments.
Genomic distribution of cyclophilins
CYP genes are distributed on all 19 chromosomes of the B. napus genome (Fig. 4a). There is no chromosome that is not encoding any CYP. Chromosomes A07 and C06 contain only one CYP each. Chromosome A01 has the largest number with 11 CYP genes. All other chromosomes contain between 2 and 10 CYP genes (Fig. 4b). The 4 CYP genes BnaAnng41240D (BnCYP12-2), BnaAnng15590D (BnCYP21-3), BnaAnng17350D (BnCYP21-6), and BnaAnng12550D (BnCYP62) were mapped to the A chromosomes, but without a detailed information regarding the exact chromosomal location. Likewise, the 4 CYP genes BnaCnng42430D (BnCYP8), BnaCnng08980D (BnCYP21-1), BnaCnng32070D (BnCYP22-1), and BnaCnng32180D (BnCYP26-5) were mapped to the C chromosomes, but again without a detailed information regarding the exact chromosomal location. The A chromosomes contain 42 CYP genes and 4 not fully assigned genes, and the C chromosomes contain 44 CYP genes and 4 not fully assigned genes. Thus, both chromosome sets (A and C), which originate from the B. rapa and B. oleracea genome, contain a similar number of CYP genes.
Interestingly, there are three CYP proteins that are each encoded by two of the six following genes BnaA05g08140D and BnaC04g09170D (BnCYP18-1), BnaA08g16920D and BnaC03g60160D (BnCYP18-4), or BnaA10g29520D and BnaC09g34020D (BnCYP21-9). Thus, one copy of these genes occurs on each of the chromosomal sets. Alignments of these CYP genes from the A and respective C chromosome showed slightly different nucleotide sequences. Nevertheless, these encode identical amino acid sequences (data not shown).
Due to the two progenitor chromosomal sets, several CYP genes encode proteins with high sequence homologies. Most often these occur as pairs, each gene originating from one progenitor chromosomal set. They are highlighted with identical colours on the A and C chromosomes in Fig. 4 and are from now on called isoforms. Some of them are located on the respective A and C chromosome, others are spread on different chromosomes, probably due to chromosomal rearrangements of the B. napus genome (as recently described by Cheng et al.68).
In summary, the allopolyploidy of oilseed rape results in two isoforms of many proteins. These show high sequence identities on the nucleotide and amino acid levels and originate from either the A or C genome. Due to their similarity, they cluster in the phylogenetic tree and often the same localization is predicted (Fig. 3). They might either possess the ability to replace each other or might be specialized for certain tasks and interaction partners.
One of the biggest groups of isoforms among the BnCYPs is the group of 18 kDa proteins. These contain a full-length CLD with a high degree of conservation, a cytosolic localization, and can be further subdivided by their homology. The phylogenetic tree reveals isoform pairs: BnCYP18-2 and BnCYP18-3, BnCYP18-5 and BnCYP18-8, and BnCYP18-6 and BnCYP18-7 (Fig. 3). Each pair has a representative gene on chromosome A and C, respectively (Fig. 4a). Interestingly, BnCYP18-4 which is closely related to the pair BnCYP18-6/18-7 is encoded by two genes located on C03 and A08 that might derive from a gene duplication event. The A. thaliana homologs of the described 18 kDa isoforms are AtCYP18-3, AtCYP18-4, and AtCYP19-1 which are shown to be important for various processes from plant growth to pathogen defence (summarized in Table 2). In contrast to the other 18 kDa BnCYPs, BnCYP18-1 is the only 18 kDa CYP with introns and without the prominent insertion between α-helix-I and β-sheet-III. Moreover, its localization is predicted to be nuclear instead of cytosolic (Fig. 3).
Expression analysis of the predicted BnCYPs on transcript and protein levels
To investigate the mRNA expression of the bioinformatically predicted CYP genes, RNA-Seq of RNA isolated from leaves was performed. This approach revealed 77 expressed BnCYPs under the applied conditions (summarized in Table 3). The expression pattern indicates a large variance in the abundance of BnCYP transcripts (Fig. 5, for raw data see Supplementary Table S2). The results thus show that most of the predicted CYP genes are indeed transcribed. Transcripts for CYP genes not detected in this study could be expressed in other tissues, specific cell types, or under different experimental conditions.
In the analysed leaf sample, the 18 kDa and 27 kDa family members belong to the strongest expressed ones, with BnCYP18-5 showing the highest read count. Interestingly, some isoform pairs had differing read counts, e.g. BnCYP18-2/BnCYP18-3, BnCYP18-4/BnCYP18-6/BnCYP18-7, BnCYP27-3/BnCYP27-5.
Since CYPs are a prominent protein class in the phloem long-distance transport system of higher plants including oilseed rape with a potential function as a molecular chaperone important for protein long-distance transport2, 41, we performed LC-MS/MS analysis of phloem protein extract and compared the CYP protein profile to that of leaf extract. Table 3 summarizes the BnCYPs identified from leaf and phloem protein extracts. In total, 26 BnCYPs could be detected at the protein level under the applied conditions with six being unique to leaves and 12 to phloem exudate. Since sieve elements do not have a functional transcription machinery, RNA-Seq was not performed with phloem samples. However, phloem samples contain a specific set of mobile RNAs that have in part been implicated with long-distance signalling43, but this was not the subject of the present study.
Most of the identified BnCYP proteins are predicted to be localized in the cytosol, but some potentially chloroplastic, mitochondrial and nuclear BnCYPs were found as well. Besides several SD CYPs, also one MD CYP, the putative chloroplastic protein BnCYP47-2, could be identified in leaf extract. The analysis of phloem sap revealed a low molecular weight CYP, BnCYP13, which is suggested to contain only a partial CLD. Figure 6 shows the abundance of the identified CYP proteins in the two examined organs (more details in Supplementary Tables S3, S4, S5, S6). Whereas BnCYP18-5, BnCYP24-1 and BnCYP27-4 were the most abundant CYP proteins in leaf extract, only BnCYP18-5 was the most abundant one in phloem sap with a considerably higher relative abundance observed than for the other identified BnCYP proteins from this compartment. Interestingly, the most abundant CYP protein in both compartments is BnCYP18-5. The finding of so many CYPs in the phloem transport system suggests an essential function of this protein family in this compartment, probably in protein transport or long-distance signalling.
Conclusions
Cyclophilins are ubiquitous proteins that constitute a multigene family in higher organisms. Exceptionally high numbers of CYPs have been found in plants, underlining their essential importance in many essential physiological processes. However, the physiological roles of most CYPs in plants are not well understood.
The present study applied bioinformatic tools for a genome-wide identification of CYPs in the important oil crop B. napus. Sequence similarity searches with known Arabidopsis CYPs, sequence alignments and CLD prediction identified a surprisingly high number of 94 CYP-coding genes. Therefore, B. napus contains the highest number of CYPs known so far. As in other plants, CYPs are predicted to be localized in all compartments, most of them being probably cytosolic. Most BnCYPs are single-domain proteins.
Transcriptome analysis confirmed the expression of 77 distinct CYPs in the B. napus cultivar ‘Drakkar’ from leaf material under normal growth conditions. The occurrence of 26 BnCYP proteins was confirmed by LC-MS/MS analysis. It is likely that additional CYPs are expressed on transcript and protein levels in different cell types, plant parts, or under different environmental conditions. Interestingly, 12 BnCYPs were exclusively found in phloem samples and not in leaf extract supporting a fundamental and specific role in this specialized compartment.
Future studies must now focus on the functional characterization of the high number of CYPs in B. napus in order to better understand the diverse roles of CYPs in oilseed rape and in plant biology in more general. In this regard, elucidating the role(s) of the phloem CYPs in protein refolding and long-distance signalling will be of special interest.
Materials and Methods
Sequence analysis
Brassica napus CYPs were identified by BLASTp of the Arabidopsis thaliana CYP18-1 (At1g01940) and CYP19-1 (At2g16600) against the B. napus genome sequence database51 deposited at NCBI (http://www.ncbi.nlm.nih.gov/)69. BLASTp searches with the remaining AtCYPs revealed the same BnCYPs as already identified. Amino acid and cDNA sequences were obtained by the European Nucleotide Archive (http://www.ebi.ac.uk/ena/data/view/PRJEB5043). To identify A. thaliana homologs for the individual BnCYPs, the amino acid sequences of BnCYPs were used as queries for a BLAST search on UniProtKB (http://www.uniprot.org/)70.
All identified proteins were analyzed for the presence of a CLD and potential additional domains with InterPro (http://www.ebi.ac.uk/interpro/)71 and drawn by CorelDRAW. The theoretical isoelectric point was determined by the ProteinProspector Tool MS-Digest (http://prospector.ucsf.edu/prospector/cgi-bin/msform.cgi?form=msdigest), subcellular localization was predicted by LocTree3 (https://rostlab.org/services/loctree3/)72 and nuclear localization signals (NLS) by using the NLS mapper (http://nls-mapper.iab.keio.ac.jp/cgi-bin/NLS_Mapper_form.cgi)73.
Chromosome mapping of the CYPs was performed using the Brassica napus Genome Browser (http://www.genoscope.cns.fr/brassicanapus/cgi-bin/gbrowse/colza/) available by Genoscope - Centre National de Séquençage and redrawn by CorelDRAW.
Protein sequence alignment and phylogenetic analysis
Sequences were aligned using ClustalOmega (http://www.ebi.ac.uk/Tools/msa/clustalo/)74 and displayed with Jalview 2.9.0b275. The secondary structure annotation is based on the structure of the crystallized human CYPA (4N1M.pdb).
The phylogenetic tree was calculated with ClustalW2 Phylogeny (http://www.ebi.ac.uk/Tools/phylogeny/clustalw2_phylogeny/) by using the multiple alignment from ClustalOmega, and subsequently processed with iTOL (http://itol.embl.de/)76, 77.
Plant material and growth conditions
Brassica napus cultivar ‘Drakkar’ plants were grown in 19 cm pots on soil (LAT-Terra Standard P, Industrie-Erdenwerk Archut, Germany) in a glasshouse under controlled conditions with 70% humidity and a 16 h/8 h light/dark (day/night) and 22 °C/18 °C (day/night) cycle. Plants were watered once per day and fertilized with 2 g/l Osmocote Exact Standard High K (Scotts, the Netherlands).
Expression profile of cyclophilins
Transcriptome data were generated by GAMAVIR, a tri-national research activity aiming at characterizing plant:virus interactions in rapeseed (ANR-13-KBBE-0005). PolyA RNA isolated from leaf disks of 6 weeks-old Brassica napus cultivar ‘Drakkar’ plants (n = 3) was sequenced using Illumina technology. Clean RNAseq reads were aligned against the re-sequenced Brassica napus cultivar ‘Drakkar’ genome using Tophat2 v2.0.1378 and counted with Samtools v1.179.
Protein extraction for proteomics
Phloem sap was collected as described previously by Giavalisco et al.2 at the inflorescence of oilseed rape plants before flowering. Phloem sap samples were collected on ice four times (n = 4), frozen in liquid nitrogen and stored at −80 °C until further processing. Proteins were precipitated in 4 volumes of 90% (v/v) acetone, 10% (v/v) methanol, 10 mM DTT over night at −20 °C. The precipitates were pelleted at 14,000 × g at 4 °C for 15 min, washed twice with 100% acetone and air-dried.
Leaf material was harvested from three different plants (n = 3) and for each sample 100 mg material was grinded in liquid nitrogen. Proteins were extracted with 800 µl extraction buffer (50 mM MOPS/pH 7.5, 5% glycerol, 0.55% PVPP, 0.5% Nonidet P-40, 5 mM L-ascorbic acid, 5 mM DTT, 1x protease inhibitor (cOmplete Protease Inhibitor Cocktail, Roche), 1x phosphatase inhibitor (PhosStop, Phosphatase Inhibitor Cocktail, Roche). Centrifuging at 14,000 × g at 4 °C for 15 min allowed the separation of soluble proteins from insoluble material. Proteins were precipitated as described above by acetone/methanol/DTT.
Analysis of phloem sap and leaf proteins by LC-MS/MS
For the analysis of phloem proteins the extracted and precipitated phloem sap proteins were resuspended in extraction buffer (6 M urea, 2 M thiourea, 15 mM DTT, 2% CHAPS). Here, each protein pellet was dissolved in exactly the same volume deployed for the initial phloem sap precipitation. Once the proteins were in solution, samples were sonicated for 10 min in a sonication bath, followed by 30 min incubation on an orbital shaker (100 rpm) at room temperature. Solubilised proteins were centrifuged at 10,000 × g for 5 min and the protein concentration was determined from the collected supernatant.
For the analysis of leaf tissue proteins, the precipitated proteins were resuspended in sufficient protein extraction buffer to provide a final protein concentration of 2 µg/µl.
25 µg of phloem sap and 50 µg of leaf tissue protein extract were then digested in solution using a Trypsin/Lys-C mixture (Mass Spec Grade, Promega, Madison, WI, USA) according to the instruction manual. After the digestion, the samples were desalted using C18-stage tips as described in Rappsilber et al.80.
After the elution of the digested and desalted peptides from C18-stage tips, the samples were concentrated to near dryness in a SpeedVac and the peptide mixtures were reconstituted in 30 µl resuspension buffer (5% acetonitrile, 0.1% formic acid in water). 5 µl of this peptide mix was injected onto an Acclaim PepMap RSLC HPLC column (75 µm × 15 cm, Thermo Scientific) connected to the EASY-nLC 1000 system (Thermo Scientific). The eluting peptides were then analyzed on a Q Exactive Plus (Thermo Scientific, Bremen, Germany) high-resolution mass spectrometer.
The peptides were separated using a binary buffer system of 0.1% formic acid in water (Buffer A) and 60% acetonitrile containing 0.1% formic (Buffer B). The flow rate was adjusted to 300 nl/min. Peptides were eluted with on a linear gradient of 0–40% buffer B for 50 min followed by a linear gradient between 40–80% buffer B for additional 30 min. Peptides were analyzed in the mass spectrometer using one full scan (300–1600 m/z, R = 70,000 at 200 m/z), followed by up to fifteen data-dependent MS/MS scans (Top 15 approach) with higher-energy collisional dissociation (HCD) at a resolution of 17,500 at 200 m/z. Dynamic exclusion was set to 15 s.
Raw data were processed using the Progenesis QI for proteomics (Progenesis QI for Proteomics Version 3.0, Nonlinear Dynamics, Newcastle, UK) software and the protein sequences of all identified CYPs from Brassica napus. Protein identifications were filtered with a false discovery rate better than 1%, at least two peptides, one unique peptide and a score of 50.
References
Buhtz, A., Pieritz, J., Springer, F. & Kehr, J. Phloem small RNAs, nutrient stress responses, and systemic mobility. BMC Plant Biology 10, 64, doi:10.1186/1471-2229-10-64 (2010).
Giavalisco, P., Kapitza, K., Kolasa, A., Buhtz, A. & Kehr, J. Towards the proteome of Brassica napus phloem sap. Proteomics 6, 896–909, doi:10.1002/pmic.200500155 (2006).
Pant, B. D., Buhtz, A., Kehr, J. & Scheible, W. R. MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. The Plant Journal 53, 731–738, doi:10.1111/tpj.2008.53.issue-5 (2008).
Kehr, J., Buhtz, A. & Giavalisco, P. Analysis of xylem sap proteins from Brassica napus. BMC Plant Biology 5, 11, doi:10.1186/1471-2229-1185-1111 (2005).
Ostendorp, A., Pahlow, S., Deke, J., Thieß, M. & Kehr, J. Protocol: optimisation of a grafting protocol for oilseed rape (Brassica napus) for studying long-distance signalling. Plant Methods 12, 1–8, doi:10.1186/s13007-016-0122-x (2016).
Kagale, S. et al. Polyploid evolution of the Brassicaceae during the Cenozoic era. Plant Cell 26, 2777–2791, doi:10.1105/tpc.114.126391 (2014).
Schmutzer, T. et al. Species-wide genome sequence and nucleotide polymorphisms from the model allopolyploid plant Brassica napus. Scientific data 2, 150072, doi:10.1038/sdata.2015.72 (2015).
Iñiguez-Luy, F. L. F., M. L. In Genetics and Genomics of the Brassicaceae (ed I.; Schmidt Bancroft, R.) 291–322 (New York, NY: Springer, 2011).
Friedt, W. & Snowdon, R. In Oil Crops Vol. 4 Handbook of Plant Breeding (eds Johann Vollmann & Istvan Rajcan) Ch. 4, 91–126 (Springer New York, 2010).
Harding, M. W., Galat, A., Uehling, D. E. & Schreiber, S. L. A receptor for the immunosuppressant FK506 is a cis-trans peptidyl-prolyl isomerase. Nature 341, 758–760, doi:10.1038/341758a0 (1989).
Bierer, B. E. et al. Two distinct signal transmission pathways in T lymphocytes are inhibited by complexes formed between an immunophilin and either FK506 or rapamycin. Proc Natl Acad Sci USA 87, 9231–9235, doi:10.1073/pnas.87.23.9231 (1990).
Handschumacher, R. E., Harding, M. W., Rice, J., Drugge, R. J. & Speicher, D. W. Cyclophilin: a specific cytosolic binding protein for cyclosporin A. Science 226, 544–547, doi:10.1126/science.6238408 (1984).
Harding, M. W., Handschumacher, R. E. & Speicher, D. W. Isolation and amino acid sequence of cyclophilin. J Biol Chem 261, 8547–8555 (1986).
Kiefhaber, T., Quaas, R., Hahn, U. & Schmid, F. X. Folding of ribonuclease T1. 2. Kinetic models for the folding and unfolding reactions. Biochemistry 29, 3061–3070, doi:10.1021/bi00464a024 (1990).
Rahfeld, J. U. et al. Confirmation of the existence of a third family among peptidyl-prolyl cis/trans isomerases. Amino acid sequence and recombinant production of parvulin. FEBS Lett 352, 180–184, doi:10.1016/0014-5793(94)00932-5 (1994).
Galat, A. Peptidylprolyl cis/trans isomerases (immunophilins): biological diversity–targets–functions. Current topics in medicinal chemistry 3, 1315–1347, doi:10.2174/1568026033451862 (2003).
Liu, J. et al. Calcineurin is a common target of cyclophilin-cyclosporin A and FKBP-FK506 complexes. Cell 66, 807–815, doi:10.1016/0092-8674(91)90124-H (1991).
Barik, S. Immunophilins: for the love of proteins. Cell Mol Life Sci 63, 2889–2900, doi:10.1007/s00018-006-6215-3 (2006).
Gasser, C. S., Gunning, D. A., Budelier, K. A. & Brown, S. M. Structure and expression of cytosolic cyclophilin/peptidyl-prolyl cis-trans isomerase of higher plants and production of active tomato cyclophilin in Escherichia coli. Proc Natl Acad Sci USA 87, 9519–9523, doi:10.1073/pnas.87.24.9519 (1990).
Galat, A. Variations of sequences and amino acid compositions of proteins that sustain their biological functions: An analysis of the cyclophilin family of proteins. Arch Biochem Biophys 371, 149–162, doi:10.1006/abbi.1999.1434 (1999).
Klappa, P., Freedman, R. B. & Zimmermann, R. Protein disulphide isomerase and a lumenal cyclophilin-type peptidyl prolyl cis-trans isomerase are in transient contact with secretory proteins during late stages of translocation. Eur J Biochem 232, 755–764, doi:10.1111/ejb.1995.232.issue-3 (1995).
Kern, G., Kern, D., Schmid, F. X. & Fischer, G. A kinetic analysis of the folding of human carbonic anhydrase II and its catalysis by cyclophilin. J Biol Chem 270, 740–745, doi:10.1074/jbc.270.2.740 (1995).
Zander, K. et al. Cyclophilin A interacts with HIV-1 Vpr and is required for its functional expression. J Biol Chem 278, 43202–43213, doi:10.1074/jbc.M305414200 (2003).
Baker, E. K., Colley, N. J. & Zuker, C. S. The cyclophilin homolog NinaA functions as a chaperone, forming a stable complex in vivo with its protein target rhodopsin. EMBO J 13, 4886–4895 (1994).
Allain, F., Denys, A. & Spik, G. Characterization of surface binding sites for cyclophilin B on a human tumor T-cell line. J Biol Chem 269, 16537–16540 (1994).
Schiene-Fischer, C. & Yu, C. Receptor accessory folding helper enzymes: the functional role of peptidyl prolyl cis/trans isomerases. FEBS Lett 495, 1–6, doi:10.1016/S0014-5793(01)02326-2 (2001).
Brazin, K. N., Mallis, R. J., Fulton, D. B. & Andreotti, A. H. Regulation of the tyrosine kinase Itk by the peptidyl-prolyl isomerase cyclophilin A. Proc Natl Acad Sci USA 99, 1899–1904, doi:10.1073/pnas.042529199 (2002).
Bukrinsky, M. I. Cyclophilins: unexpected messengers in intercellular communications. Trends Immunol 23, 323–325, doi:10.1016/S1471-4906(02)02237-8 (2002).
Jing, H. et al. Peptidyl-prolyl isomerization targets rice Aux/IAAs for proteasomal degradation during auxin signalling. Nat Commun 6, 7395, doi:10.1038/ncomms8395 (2015).
Pogorelko, G. V. et al. Characterization of three Arabidopsis thaliana immunophilin genes involved in the plant defense response against Pseudomonas syringae. Gene 538, 12–22, doi:10.1016/j.gene.2014.01.029 (2014).
Lin, D. T. & Lechleiter, J. D. Mitochondrial targeted cyclophilin D protects cells from cell death by peptidyl prolyl isomerization. J Biol Chem 277, 31134–31141, doi:10.1074/jbc.M112035200 (2002).
Krzywicka, A. et al. KIN241: a gene involved in cell morphogenesis in Paramecium tetraurelia reveals a novel protein family of cyclophilin-RNA interacting proteins (CRIPs) conserved from fission yeast to man. Mol Microbiol 42, 257–267, doi:10.1046/j.1365-2958.2001.02634.x (2001).
Anderson, M. et al. A new family of cyclophilins with an RNA recognition motif that interact with members of the trx/MLL protein family in Drosophila and human cells. Dev Genes Evol 212, 107–113, doi:10.1007/s00427-002-0213-8 (2002).
Dubourg, B. et al. The human nuclear SRcyp is a cell cycle-regulated cyclophilin. J Biol Chem 279, 22322–22330, doi:10.1074/jbc.M400736200 (2004).
Li, H. et al. A WD40 domain cyclophilin interacts with histone H3 and functions in gene repression and organogenesis in Arabidopsis. Plant Cell 19, 2403–2416, doi:10.1105/tpc.107.053579 (2007).
Bannikova, O. et al. Identification of RNA targets for the nuclear multidomain cyclophilin atCyp59 and their effect on PPIase activity. Nucleic Acids Res 41, 1783–1796, doi:10.1093/nar/gks1252 (2013).
Ahn, J. C. et al. Classification of rice (Oryza sativa L. Japonica nipponbare) immunophilins (FKBPs, CYPs) and expression patterns under water stress. BMC Plant Biol 10, 253, doi:10.1186/1471-2229-10-253 (2010).
Dominguez-Solis, J. R. et al. A cyclophilin links redox and light signals to cysteine biosynthesis and stress responses in chloroplasts. Proc Natl Acad Sci USA 105, 16386–16391, doi:10.1073/pnas.0808204105 (2008).
Ruan, S. L. et al. Proteomic identification of OsCYP2, a rice cyclophilin that confers salt tolerance in rice (Oryza sativa L.) seedlings when overexpressed. BMC Plant Biol 11, 34, doi:10.1186/1471-2229-11-34 (2011).
Park, S. W. et al. Cyclophilin 20-3 relays a 12-oxo-phytodienoic acid signal during stress responsive regulation of cellular redox homeostasis. Proc Natl Acad Sci USA 110, 9559–9564, doi:10.1073/pnas.1218872110 (2013).
Gottschalk, M. et al. Ricinus communis cyclophilin: functional characterisation of a sieve tube protein involved in protein folding. Planta 228, 687–700, doi:10.1007/s00425-008-0771-8 (2008).
Esau, K. The phloem. Vol. 5 (Gebr. Borntraeger, 1969).
Kehr, J. Long-distance transport of macromolecules through the phloem. F1000 Biol Rep 1, 31, doi:10.3410/B1-31 (2009).
Galat, A. Function-dependent clustering of orthologues and paralogues of cyclophilins. Proteins 56, 808–820, doi:10.1002/prot.20156 (2004).
Pemberton, T. J. & Kay, J. E. The cyclophilin repertoire of the fission yeast Schizosaccharomyces pombe. Yeast 22, 927–945, doi:10.1002/yea.1288 (2005).
Singh, K., Zouhar, M., Mazakova, J. & Rysanek, P. Genome wide identification of the immunophilin gene family in Leptosphaeria maculans: a causal agent of Blackleg disease in Oilseed Rape (Brassica napus). OMICS 18, 645–657, doi:10.1089/omi.2014.0081 (2014).
Romano, P. G., Horton, P. & Gray, J. E. The Arabidopsis cyclophilin gene family. Plant Physiol 134, 1268–1282, doi:10.1104/pp.103.022160 (2004).
He, Z., Li, L. & Luan, S. Immunophilins and parvulins. Superfamily of peptidyl prolyl isomerases in Arabidopsis. Plant Physiol 134, 1248–1267, doi:10.1104/pp.103.031005 (2004).
Vasudevan, D. et al. Plant immunophilins: a review of their structure-function relationship. Biochim Biophys Acta 1850, 2145–58, doi:10.1016/j.bbagen.2014.12.017 (2014).
Mainali, H. R., Chapman, P. & Dhaubhadel, S. Genome-wide analysis of Cyclophilin gene family in soybean (Glycine max). BMC Plant Biol 14, 282, doi:10.1186/s12870-014-0282-7 (2014).
Chalhoub, B. et al. Plant genetics. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 345, 950–953, doi:10.1126/science.1253435 (2014).
Ke, H. M., Zydowsky, L. D., Liu, J. & Walsh, C. T. Crystal structure of recombinant human T-cell cyclophilin A at 2.5 A resolution. Proc Natl Acad Sci USA 88, 9483–9487, doi:10.1073/pnas.88.21.9483 (1991).
Zydowsky, L. D. et al. Active site mutants of human cyclophilin A separate peptidyl-prolyl isomerase activity from cyclosporin A binding and calcineurin inhibition. Protein Sci 1, 1092–1099, doi:10.1002/pro.5560010903 (1992).
Liu, J., Chen, C. M. & Walsh, C. T. Human and Escherichia coli cyclophilins: sensitivity to inhibition by the immunosuppressant cyclosporin A correlates with a specific tryptophan residue. Biochemistry 30, 2306–2310, doi:10.1021/bi00223a003 (1991).
Lee, S. S. et al. Rice cyclophilin OsCYP18-2 is translocated to the nucleus by an interaction with SKIP and enhances drought tolerance in rice and Arabidopsis. Plant Cell Environ 38, 2071–87, doi:10.1111/pce.12531 (2015).
Iki, T., Yoshikawa, M., Meshi, T. & Ishikawa, M. Cyclophilin 40 facilitates HSP90-mediated RISC assembly in plants. EMBO J 31, 267–278, doi:10.1038/emboj.2011.395 (2012).
Earley, K. W. & Poethig, R. S. Binding of the cyclophilin 40 ortholog SQUINT to Hsp90 protein is required for SQUINT function in Arabidopsis. J Biol Chem 286, 38184–38189, doi:10.1074/jbc.M111.290130 (2011).
Ratajczak, T. & Carrello, A. Cyclophilin 40 (CyP-40), mapping of its hsp90 binding domain and evidence that FKBP52 competes with CyP-40 for hsp90 binding. J Biol Chem 271, 2961–2965, doi:10.1074/jbc.271.6.2961 (1996).
Nelson, C. J., Santos-Rosa, H. & Kouzarides, T. Proline isomerization of histone H3 regulates lysine methylation and gene expression. Cell 126, 905–916, doi:10.1016/j.cell.2006.07.026 (2006).
Li, H. & Luan, S. The cyclophilin AtCYP71 interacts with CAF-1 and LHP1 and functions in multiple chromatin remodeling processes. Mol Plant 4, 748–758, doi:10.1093/mp/ssr036 (2011).
Vasudevan, D., Fu, A., Luan, S. & Swaminathan, K. Crystal structure of Arabidopsis cyclophilin38 reveals a previously uncharacterized immunophilin fold and a possible autoinhibitory mechanism. Plant Cell 24, 2666–2674, doi:10.1105/tpc.111.093781 (2012).
Shapiguzov, A., Edvardsson, A. & Vener, A. V. Profound redox sensitivity of peptidyl-prolyl isomerase activity in Arabidopsis thylakoid lumen. FEBS Lett 580, 3671–3676, doi:10.1016/j.febslet.2006.05.054 (2006).
Edvardsson, A., Eshaghi, S., Vener, A. V. & Andersson, B. The major peptidyl-prolyl isomerase activity in thylakoid lumen of plant chloroplasts belongs to a novel cyclophilin TLP20. FEBS Lett 542, 137–141, doi:10.1016/S0014-5793(03)00366-1 (2003).
Fu, A. et al. A chloroplast cyclophilin functions in the assembly and maintenance of photosystem II in Arabidopsis thaliana. Proc Natl Acad Sci USA 104, 15947–15952, doi:10.1073/pnas.0707851104 (2007).
Sirpio, S. et al. AtCYP38 ensures early biogenesis, correct assembly and sustenance of photosystem II. Plant J 55, 639–651, doi:10.1111/j.1365-313X.2008.03532.x (2008).
Weighardt, F. et al. A novel hnRNP protein (HAP/SAF-B) enters a subset of hnRNP complexes and relocates in nuclear granules in response to heat shock. J Cell Sci 112 (Pt 10), 1465–1476 (1999).
Gullerova, M., Barta, A. & Lorkovic, Z. J. AtCyp59 is a multidomain cyclophilin from Arabidopsis thaliana that interacts with SR proteins and the C-terminal domain of the RNA polymerase II. RNA 12, 631–643, doi:10.1261/rna.2226106 (2006).
Cheng, F., Wu, J. & Wang, X. Genome triplication drove the diversification of Brassica plants. Horticulture research 1, 14024, doi:10.1038/hortres.2014.24 (2014).
Coordinators, N. R. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res 44, D7–19, doi:10.1093/nar/gkv1290 (2016).
UniProt, C. UniProt: a hub for protein information. Nucleic Acids Res 43, D204–212, doi:10.1093/nar/gku989 (2015).
Mitchell, A. et al. The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Res 43, D213–221, doi:10.1093/nar/gku1243 (2015).
Goldberg, T. et al. LocTree3 prediction of localization. Nucleic Acids Res 42, W350–355, doi:10.1093/nar/gku396 (2014).
Kosugi, S., Hasebe, M., Tomita, M. & Yanagawa, H. Systematic identification of cell cycle-dependent yeast nucleocytoplasmic shuttling proteins by prediction of composite motifs. Proc Natl Acad Sci USA 106, 10171–10176, doi:10.1073/pnas.0900604106 (2009).
Sievers, F. et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7, 539–539, doi:10.1038/msb.2011.75 (2011).
Waterhouse, A. M., Procter, J. B., Martin, D. M., Clamp, M. & Barton, G. J. Jalview Version 2–a multiple sequence alignment editor and analysis workbench. Bioinformatics 25, 1189–1191, doi:10.1093/bioinformatics/btp033 (2009).
Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics 23, 127–128, doi:10.1093/bioinformatics/btl529 (2007).
Letunic, I. & Bork, P. Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res 39, W475–478, doi:10.1093/nar/gkr201 (2011).
Kim, D. et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biology 14, R36, doi:10.1186/gb-2013-14-4-r36 (2013).
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079, doi:10.1093/bioinformatics/btp352 (2009).
Rappsilber, J., Mann, M. & Ishihama, Y. Protocol for micro-purification, enrichment, pre-fractionation and storage of peptides for proteomics using StageTips. Nat Protoc 2, 1896–1906, doi:10.1038/nprot.2007.261 (2007).
Sakuma, Y. et al. Dual function of an Arabidopsis transcription factor DREB2A in water-stress-responsive and heat-stress-responsive gene expression. Proc Natl Acad Sci USA 103, 18822–18827, doi:10.1073/pnas.0605639103 (2006).
Trupkin, S. A., Mora-Garcia, S. & Casal, J. J. The cyclophilin ROC1 links phytochrome and cryptochrome to brassinosteroid sensitivity. Plant J 71, 712–723, doi:10.1111/j.1365-313X.2012.05013.x (2012).
Li, M. et al. Proline isomerization of the immune receptor-interacting protein RIN4 by a cyclophilin inhibits effector-triggered immunity in Arabidopsis. Cell Host Microbe 16, 473–483, doi:10.1016/j.chom.2014.09.007 (2014).
Kovalev, N. & Nagy, P. D. Cyclophilin A binds to the viral RNA and replication proteins, resulting in inhibition of tombusviral replicase assembly. J Virol 87, 13330–13342, doi:10.1128/JVI.02101-13 (2013).
Deng, W. et al. Agrobacterium VirD2 protein interacts with plant host cyclophilins. Proc Natl Acad Sci USA 95, 7040–7045, doi:10.1073/pnas.95.12.7040 (1998).
Ma, X., Song, L., Yang, Y. & Liu, D. A gain-of-function mutation in the ROC1 gene alters plant architecture in Arabidopsis. New Phytol 197, 751–762, doi:10.1111/nph.12056 (2013).
Coaker, G., Falick, A. & Staskawicz, B. Activation of a phytopathogenic bacterial effector protein by a eukaryotic cyclophilin. Science 308, 548–550, doi:10.1126/science.1108633 (2005).
Coaker, G., Zhu, G., Ding, Z., Van Doren, S. R. & Staskawicz, B. Eukaryotic cyclophilin as a molecular switch for effector activation. Mol Microbiol 61, 1485–1496, doi:10.1111/j.1365-2958.2006.05335.x (2006).
Aumuller, T., Jahreis, G., Fischer, G. & Schiene-Fischer, C. Role of prolyl cis/trans isomers in cyclophilin-assisted Pseudomonas syringae AvrRpt2 protease activation. Biochemistry 49, 1042–1052, doi:10.1021/bi901813e (2010).
Stangeland, B. et al. Molecular analysis of Arabidopsis endosperm and embryo promoter trap lines: reporter-gene expression can result from T-DNA insertions in antisense orientation, in introns and in intergenic regions, in addition to sense insertion at the 5′ end of genes. J Exp Bot 56, 2495–2505, doi:10.1093/jxb/eri242 (2005).
Kaur, G. et al. Characterization of Peptidyl-Prolyl Cis-Trans Isomerase- and Calmodulin-Binding Activity of a Cytosolic Arabidopsis thaliana Cyclophilin AtCyp19-3. PloS one 10, e0136692, doi:10.1371/journal.pone.0136692 (2015).
Saito, T. et al. Expression of a gene for cyclophilin which contains an amino-terminal endoplasmic reticulum-targeting signal. Plant Cell Physiol 40, 77–87, doi:10.1093/oxfordjournals.pcp.a029477 (1999).
Grebe, M. et al. A conserved domain of the arabidopsis GNOM protein mediates subunit interaction and cyclophilin 5 binding. Plant Cell 12, 343–356, doi:10.1105/tpc.12.3.343 (2000).
Kamauchi, S., Nakatani, H., Nakano, C. & Urade, R. Gene expression in response to endoplasmic reticulum stress in Arabidopsis thaliana. FEBS J 272, 3461–3476, doi:10.1111/j.1742-4658.2005.04770.x (2005).
Jackson, K. & Soll, D. Mutations in a new Arabidopsis cyclophilin disrupt its interaction with protein phosphatase 2A. Mol Gen Genet 262, 830–838, doi:10.1007/s004380051147 (1999).
Schubert, M. et al. Proteome map of the chloroplast lumen of Arabidopsis thaliana. J Biol Chem 277, 8354–8365, doi:10.1074/jbc.M108575200 (2002).
Ingelsson, B., Shapiguzov, A., Kieselbach, T. & Vener, A. V. Peptidyl-prolyl isomerase activity in chloroplast thylakoid lumen is a dispensable function of immunophilins in Arabidopsis thaliana. Plant Cell Physiol 50, 1801–1814, doi:10.1093/pcp/pcp122 (2009).
Zhang, Y. et al. The cyclophilin CYP20-2 modulates the conformation of BRASSINAZOLE-RESISTANT1, which binds the promoter of FLOWERING LOCUS D to regulate flowering in Arabidopsis. Plant Cell 25, 2504–2521, doi:10.1105/tpc.113.110296 (2013).
Sirpio, S., Holmstrom, M., Battchikova, N. & Aro, E. M. AtCYP20-2 is an auxiliary protein of the chloroplast NAD(P)H dehydrogenase complex. FEBS Lett 583, 2355–2358, doi:10.1016/j.febslet.2009.06.031 (2009).
Edvardsson, A., Shapiguzov, A., Petersson, U. A., Schroder, W. P. & Vener, A. V. Immunophilin AtFKBP13 sustains all peptidyl-prolyl isomerase activity in the thylakoid lumen from Arabidopsis thaliana deficient in AtCYP20-2. Biochemistry 46, 9432–9442, doi:10.1021/bi700426q (2007).
Lippuner, V. et al. Cloning and characterization of chloroplast and cytosolic forms of cyclophilin from Arabidopsis thaliana. J Biol Chem 269, 7863–7868 (1994).
Motohashi, K., Koyama, F., Nakanishi, Y., Ueoka-Nakanishi, H. & Hisabori, T. Chloroplast cyclophilin is a target protein of thioredoxin. J Biol Chem 278, 31848–31852, doi:10.1074/jbc.M304258200 (2003).
Laxa, M., Konig, J., Dietz, K. J. & Kandlbinder, A. Role of the cysteine residues in Arabidopsis thaliana cyclophilin CYP20-3 in peptidyl-prolyl cis-trans isomerase and redox-related functions. Biochem J 401, 287–297, doi:10.1042/BJ20061092 (2007).
Liebthal, M. et al. Redox-Dependent Conformational Dynamics of Decameric 2-Cysteine Peroxiredoxin and its Interaction with Cyclophilin 20-3. Plant and Cell Physiology 57, 1415–1425, doi:10.1093/pcp/pcw031 (2016).
Kim, J. H., Nguyen, N. H., Nguyen, N. T., Hong, S. W. & Lee, H. Loss of all three calreticulins, CRT1, CRT2 and CRT3, causes enhanced sensitivity to water stress in Arabidopsis. Plant Cell Rep 32, 1843–1853, doi:10.1007/s00299-013-1497-z (2013).
Tomasic Paic, A. & Fulgosi, H. Chloroplast immunophilins. Protoplasma 253, 249–58, doi:10.1007/s00709-015-0828-z (2015).
Smith, M. R. et al. Cyclophilin 40 is required for microRNA activity in Arabidopsis. Proc Natl Acad Sci USA 106, 5424–5429, doi:10.1073/pnas.0812729106 (2009).
Acknowledgements
We are grateful to S. Graindorge and D. Pflieger at the IBMP BioImage and Bioinformatics platform for providing access to the bioinformatic resources and the analysis of RNAseq data. We would like to acknowledge the financial contribution to the research activities by a Career Integration Grant (CIG; PCIG14-GA-2013-63 0734) by the European Commission within the 7th framework program, the grant LFF-GK06 ‘DELIGRAH’ (Landesforschungsförderung Hamburg), and a DFG grant (DFG KE 856_6-1) awarded to J.K., and funding by the Agence National de la Recherche (ANR-13-KBBE-0005-01) awarded to M.H.
Author information
Authors and Affiliations
Contributions
P.H. and M.T. performed the original database searches and contributed to the design of the project. P.H. performed bioinformatic predictions, phloem sap sampling and protein extraction. K.A. performed RNA extraction, K.A. and M.H. contributed to RNA-Seq data analysis. K.B. and P.G. performed LC-MS/MS. P.H. and P.G. analyzed the protein data. J.K. and M.H. participated in data analysis and discussion. P.H. and M.T. drafted the initial manuscript, J.K. conceived of and supervised the study and complemented the writing. All authors reviewed the manuscript and approved the final article.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hanhart, P., Thieß, M., Amari, K. et al. Bioinformatic and expression analysis of the Brassica napus L. cyclophilins. Sci Rep 7, 1514 (2017). https://doi.org/10.1038/s41598-017-01596-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-01596-5
This article is cited by
-
Arabidopsis cyclophilins direct intracellular transport of mobile mRNA via organelle hitchhiking
Nature Plants (2024)
-
Characterization and regulation of salt upregulated cyclophilin from a halotolerant strain of Penicillium oxalicum
Scientific Reports (2023)
-
Genome-wide identification and stress response analysis of cyclophilin gene family in apple (Malus × domestica)
BMC Genomics (2022)
-
Genome-wide characterization of peptidyl-prolyl cis–trans isomerases in Penicillium and their regulation by salt stress in a halotolerant P. oxalicum
Scientific Reports (2021)
-
Enzyme activity and structural features of three single-domain phloem cyclophilins from Brassica napus
Scientific Reports (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.