Identification and Analysis of Expression Patterns of the Caleosin Genes in Hickory (Carya cathayensis Sarg.)

By inergency On Mar 28, 2024

1. Introduction

Lipids are essential for various cellular functions, including energy homeostasis, membrane remodeling, and lipid signaling [1]. In plant cells, lipids are primarily found in the form of triacylglycerol (TAG), which is synthesized in the endoplasmic reticulum. These TAGs form individual oil bodies on the endoplasmic reticulum membrane, which, via the acquisition of specific proteins, support the growth and expansion of oil bodies [2]. These oil bodies play a role in providing energy for seed germination and early seedling growth and development [3]. The oil body is one of the smallest organelles in plant cells, typically measuring between 0.5–2.5 μm in diameter; however, in some species, such as soybean (Glycine max), the diameter of the oil body can reach 5–7.5 μm [4]. In contrast to organelles such as mitochondria and chloroplasts, which are surrounded by a phospholipid bilayer, oil bodies are characterized by a phospholipid monolayer and a high concentration of triacylglycerol, along with embedded proteins.

Proteomics analyses have identified 100–150 proteins in oil bodies isolated from mammalian cells, primarily associated with metabolic pathways concerning lipid synthesis, degradation, membrane transport, and protein degradation [5]. In seed plants, the predominant surface mosaic proteins of oil bodies are oleosin (OLEO), caleosin (CLO), and steroleosin (SLO) [6]. Proteomics have unveiled low-abundance oil body proteins such as LD-ASSOCIATED PROTEIN (LDAP) and LDAP-INTERACTING PROTEIN (LDIP) [7,8]. Notably, oleosin is the most prevalent oil body surface protein, constituting 79% of the total proteins [9]. In seeds, oleosin proteins play a crucial role in maintaining the integrity and stability of oil bodies by preventing their aggregation. The efficient mobilization of TAGs in plant oil bodies necessitates the degradation of oleosin proteins. Therefore, ensuring the degradation of oleosin proteins is crucial for the effective mobilization of TAGs [10,11]. Four important oleosin genes have been identified in Brassica napus, and it was found that the overexpression of BnaOLE1, BnaOLE2, and BnaOLE4 led to an increase in the oil content of the transgenic lines in Arabidopsis thaliana. Moreover, the overexpression of BnaOLE2 and BnaOLE4 significantly increased the size of the oil bodies and the linoleic acid content, as well as the thousand-grain weight of Arabidopsis thaliana seeds [12].

Caleosin (CLO/PXGs) is another highly abundant oil-body protein that possesses a Ca²⁺ binding domain and is initially identified within the oil bodies of sesame seeds [13]. Caleosin plays important roles in various developmental stages, including plant seed germination, seedling tissue differentiation, leaf senescence, pollen maturation, and seed maturation [14]. Both oleosin and caleosin contribute to the structural stability of oil bodies [15], with caleosin capable of assuming this role in the absence of oleosin [16]. During seed germination, caleosin is involved in the degradation of storage lipids in oil bodies. Arabidopsis mutants (Atclo1-1, Atclo1-2) exhibited a reduction in the number of oil bodies present in vacuole at 48 h and 60 h post-germination compared to wild-type, resulting in a relative delay in the breakdown process of lipids [17]. Introduction of the caleosin gene from the oilseed plant Ricinus communis into tobacco leaves resulted in an 8.1% increase in oil content compared to wild-type tobacco leaves [18]. In addition, caleosins play a crucial role in plant response to adversity stress. For example, in rice (Oryza sativa), the OsClo5 gene enhances the cold sensitivity of plants by inhibiting jasmonic acid signaling and synthesis [19]. Additionally, the Atclo3 gene is strongly expressed under abiotic stresses such as drought, high salt, and ABA treatment, and it was further found that clo3 mutants were sensitive to drought and salt stresses [20]. Thus, these results support the involvement of caleosin in regulating oil body stability and lipid degradation in plant species, as well as playing a role in abiotic stresses on plants.

Hickory (Carya cathayensis) and pecan (Carya illinoensis) are part of the hickory genus in the pecan family. Carya cathayensis, known as Chinese hickory, is economically important for its high phenolic content and antioxidant capacity, the highest among major tree nuts [21]. Hickory has a high oil content of 70%, with over 90% unsaturated fatty acids, including oleic acid, linoleic acid, and linolenic acid [22,23]. The molecular mechanism underlying the high oil content of hickory kernels remains poorly understood. Previous research revealed an increase in the expression of oleosin, caleosin, and steroleosin-related genes during hickory seed kernels seed development [24]. Specifically, seven oleosin genes in the hickory oleosin gene family showed varying levels of up-regulated expression at different stages of embryo development, with Cca0520S0025 and Cca0791S0002 showing significant up-regulation during the rapid growth phase of hickory embryo oil accumulation, suggesting its potential involvement in the process of hickory oil accumulation [25]. However, the caleosin genes involved in the oil accumulation of hickory remain unclear. In this study, the commercially cultivated hickory nut species Carya cathayensis (representative of the East Asian hickory genus) was taken as the research object, with Carya illinoensi (representative of the eastern North American hickory genus) and Juglans regia (Walnut) as the reference species. Based on published genomic information, the caleosin genes were identified and analyzed using bioinformatics. To further screen candidate caleosin genes involved in embryo development of hickory in order to provide basic data for the molecular mechanism of high oil accumulation and withstand adverse conditions in hickory.

2. Results

2.1. Identification of the Carya cathayensi Caleosin Genes

The Arabidopsis caleosin genes were used as bait sequences [26], and local blasts were performed among the hickory, pecan, and walnut protein libraries. The results (Table 1) showed that four caleosin proteins were identified in the hickory, and three caleosin proteins were identified in both the pecan and walnut protein libraries.

The physicochemical properties analysis of the proteins in the hickory caleosin genes revealed a range of amino acids encoded (142aa to 239aa), molecular weights (16.90 kDa to 26.86 kDa), isoelectric points (5.88 to 9.26), and instability coefficients (47.74 to 38.15). All caleosin proteins were identified as hydrophilic proteins. Subcellular localization predictions indicated that CcaCLO2 and CcCLO3 were located in the nucleus, CcaCLO1 in the chloroplasts, and CcaCLO4 in the cytoplasm. However, CLO4 exhibited notable disparities in amino acid length compared to the other three caleosin proteins and lacked significant structural characteristics, thus justifying its omission from further analysis.

The Carya illinoinensis caleosin proteins exhibited a range of amino acids encoded (203aa to 239aa), protein molecular weight (23.03 kDa to 27.08 kDa), isoelectric point (5.89 to 9.37), instability index (37.88 to 47.72), and fatty acidity coefficient (69.91 to 78.41). Additionally, the walnut caleosin genes encoded (239aa to 201aa) amino acids, along with protein molecular weight of (22.80 kDa to 26.93 kDa), isoelectric point of (5.70 to 9.40), instability coefficient of (40.35 to 45.71), and Aliphatic index of (75.22 to 80.41). Furthermore, the caleosins of the Carya illinoensis and walnut are also hydrophilic proteins. Subcellular localization predictions for Carya illinoinensis and walnut caleosins indicated that CilCLO2, CilCLO3, JreCLO1, and JreCLO2 were found in the nucleus, while CilCLO3 and JreCLO3 were located in the chloroplast.

2.2. Phylogenetic Tree Analysis

The evolutionary relationships of caleosin proteins, including those from Carya cathayensis, Carya illinoinensis, Juglans regia, Arabidopsis thaliana, Triticum aestivum, Oryza sativa, and Glycine max, were analyzed using MEGA 7.0 software, revealing a total of 45 caleosin proteins (Figure 1). The caleosin proteins can be mainly grouped into three subfamilies. The results indicate that CcaCLO1 (Cca0849S0013) belongs to the group I, while CcaCLO2 (Cca0601S0075) and CcaCLO3 (Cca0899S0091) belong to the group II. CcaCLO1 (Cca0849S0013), CilCLO3 (Cil0992S0011), JreCLO3 (JreChr01G11961) are part of the group I, while members CcaCLO2 (Cca0601S0075), CcaCLO3 (Cca0899S0091), CilCLO1 (Cil0951S0139), CilCLO2 (Cil0001S0015), JreCLO1 (JreChr05G11280), JreCLO2 (JreChr06G12167) cluster together within the group II subfamily.

Furthermore, the results demonstrate that monocotyledons such as maize, wheat, and rice form a distinct cluster in group III, while dicotyledons, including Carya cathayensis, Carya illinoinensis, Juglans regia, and Glycine max, are grouped together in group II. Notably, group I comprises both monocotyledons and dicotyledons.

2.3. Identification of Conserved Structural Domains and Motifs of Caleosin Proteins

A structural domain resembling the oleosin “proline knot” was identified in sesame, situated within the central hydrophobic region of caleosin and contributing to oil body formation. [27]. The “proline knot” of caleosin can be represented as “PX3PSPX2P”. Conserved proline junction motifs were examined in Carya cathayensis, Carya illinoinensis, Juglans regia, and Arabidopsis thaliana. The N-terminal domain of caleosin contains a single Ca²⁺ bing motif, and the C-terminal domain contains one tyrosine kinase and three casein kinase II phosphorylation sites (Figure 2A).

Analysis of the conserved motifs revealed nine highly conserved motifs among hickory caleosin family members. CcaCLO1, CcaCLO2, and CcaCLO3 all possess motifs 1~3, 5, 7, 8 (Figure 2B). CilCLO3 from Carya illinoensis and JreCLO3 from walnut are grouped together with CcaCLO1, and all three contain motifs 1–5, 7–9. Additionally, CilCLO1 and JreCLO2 clustered in a single unit with CcaCLO2, and all three contain motifs 1–3, 5–8. CilCLO2 and JreCLO1 clustered in a single unit with CcaCLO3, and all three contain motifs 1–3, 5–8.

The gene structure analysis revealed that three hickory caleosin proteins exhibited a six-exon configuration. Similarly, all members of the Carya illinoensis caleosin family displayed a six-exon structure, while walnut caleosins were characterized by five exons.

2.4. Analysis of Promoter-Acting Elements of the Caleosin Genes

A total of 13 acting elements were identified in Carya cathayensis that were implicated in stress response, hormone signaling, and tissue expression. Specifically, these elements included regulatory elements involved in light responsiveness elements (13), regulatory elements essential for anaerobic induction (3), drought-induced MYB binding sites (2), acting elements associated with defense and stress response (1), and cis-acting elements involved in low-temperature response (1) (Figure 3). Analysis of cis-regulatory elements in promoters of caleosin genes of Carya cathayensis, Carya illinoinensis, and Juglans regiaIn terms of hormone signaling, there are a total of 10 cis-acting elements involved in response to abscisic acid (ABA) and 10 cis-acting regulatory elements involved in response to jasmonate (MeJA). Moreover, there are three cis-acting elements associated with the gibberellins (GA) response and two cis-acting elements involved in the response to salicylic acid (SA). Additionally, a comparative analysis of cis-acting elements in the promoters of the caleosins in Carya illinoinensis and Juglans regia revealed the presence of 15 acting elements, encompassing responses to stress, hormone signaling, and tissue expression. Interestingly, a cis-acting regulatory element involved in seed-specific regulation was only found in Carya cathayensis.

2.5. Structural Analysis of Caleosin Proteins in Carya cathayensis

The caleosin protein structures of Carya cathayensis were predicted using the online prediction software SOPMA (https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=/NPSA/npsa_sopma.html, accessed on 26 May 2023). The results indicate that hickory caleosin proteins are composed of four structural elements: α-helix, β-turns, extended chains, and irregular curls. The α-helix and irregular coils represent the highest percentages in CcaCLO1, at 40.59% and 43.10%, respectively. In the beta-turn structure, CcaCLO3 has the highest percentage, while CcaCLO2 exhibits the highest percentage of extended strands (Table 2). Furthermore, the analysis of transmembrane structures indicates that hickory caleosin proteins exhibit characteristics of transmembrane proteins.

The tertiary structure of the caleosin proteins in Carya cathayensis was investigated using the online prediction software SWISS-MODEL (https://swissmodel.expasy.org/, accessed on 26 May 2023). The modeling results demonstrated a high level of credibility, with GMQE scores exceeding 0.9 for all proteins. The hickory caleosin proteins exist in their active form as monomers, and analysis revealed the presence of a ‘proline knot’ in CcaCLO1, CcaCLO2, and CcaCLO3 proteins (Figure 4).

2.6. Analysis of the Expression Pattern of the Hickory Caleosin Genes during Embryo Development

The expression profile of the caleosin genes identified in the Carya cathayensis genome was evaluated for embryos at various stages of development. The findings from the qRT-PCR expression levels were consistent with the trends observed in transcriptome data (Figure 5). The investigation revealed that the gene expression of CcaCLO1 and CcaCLO3 increased as hickory embryos developed and matured, while CcaCLO2 decreased. Additionally, the transcriptome data indicated that the expression of CcaCLO1 was significantly higher than that of CcaCLO2 and CcaCLO3 during hickory embryo development. Furthermore, the qRT-PCR results demonstrated a significant up-regulation of CcaCLO1 starting from the S2 stage, which corresponds to the rapid accumulation stage of oil content [24]. This suggests that the CcaCLO1 gene in hickory is responsive to the embryonic developmental process.

2.7. Subcellular Localization of CcaCLO1 and CcaCLO2

We previously showed that CcaCLO1 is localized to the oil body using the Arabidopsis transient transformation system [28]. To further verify Carya cathayensis, we chose to use the callus of Carya cathayensis for transient transformation. Given the low efficiency of transient transformation of hickory callus, we only found CcaCLO1 localized to the oil body, which is in line with our previous results. Moreover, using the transient transformation system of Arabidopsis, we also found that CcaCLO2 was localized to the oil body. As a positive control for oil body localization, we utilized Nile red, an oil body-specific dye. Our results demonstrated that the green fluorescence emitted by CcaCLO1 and CcaCLO2 co-localized with the red fluorescence of Nile red, indicating that CcaCLO1 and CcaCLO2 are localized within the oil body (Figure 6).

2.8. Analysis of the Expression Pattern of the Hickory Caleosin Genes under Salt Stress

In the presence of 150 mM NaCl treatment, the expression levels of Cca0849S0013 (CcaCLO1), Cca0601S0075 (CcaCLO2), and Cca0899S0091 (CcaCLO3) exhibited distinct temporal patterns. Specifically, CcaCLO1 displayed a gradual increase in expression over time, reaching its peak at 12 h post-treatment. In contrast, CcaCLO2 demonstrated a fluctuating expression profile, with peak levels observed at 4 h followed by a subsequent decrease and subsequent increase. Similarly, CcaCLO3 exhibited an initial rise in expression levels, followed by a decline, with peak expression occurring at 4 h post-treatment (Figure 7).

3. Discussion

Caleosin genes have been identified in multiple plant species, such as Arabidopsis thaliana, Brassica napus, Glycine max, and Oryza sativa. The number of genes identified in each species varies, with 8, 22, 7, and 8 genes, respectively [26,29,30,31]. These findings indicate the presence of caleosin in diverse plant species, with significant variation in the number of caleosin family members. In the genome of hazelnut of Betulaceae (Corylus heterophylla), five caleosin genes were identified [32]. In this study, we used bioinformatics to identify three caleosin proteins in the genomes of hickory, pecan, and walnut. The molecular weights of these caleosin proteins ranged from 20–27 kD, and CcaCLO1, CilCLO3, and JreCLO3 belong to the H-caleosin protein group, in accordance with Shen [26]. Previous research suggested that caleosins evolved differently in Magnoliaceae, Sphagnums, and Bryophytes [33]. Forty-five species of magnoliophyte plants were analyzed in a phylogenetic study, revealing that caleosins have evolved differently in monocotyledons and dicotyledons. Furthermore, the caleosins in Carya cathayensis, Carya illinoinensis, and Juglans regia are more closely related to each other than to other species.

Higher plant caleosins are characterized by three distinct regions: a hydrophilic N-terminal region containing a single Ca²⁺-binding EF-hand, a central hydrophobic domain known as a “proline knot”, and a C-terminal hydrophilic region with four potential phosphorylation sites [27]. We also found three distinct regions: Carya cathayensis, Carya illinoinensis, and Juglans regia. The conserved structural features of caleosin suggest that its function is likely influenced by calcium binding and phosphorylation. While phosphorylation by casein kinase II has been observed in plants, tyrosine phosphorylation is considered less probable [34].

Caleosins are significant members of the proteins of the oil body, and they play a crucial role in its functions in oil accumulation and embryo development. Liu discovered that doubly protruding seeds (clo1 × clo2) had a significant reduction in oil content and seed weight in Arabidopsis. In addition, overexpression of CLO1 in seedlings and BY2 cells increased triacylglycerol content up to 73.6% [35]. The expression of CcaCLO1 (Cca0849S0013) significantly increased during embryo development, indicating a potential response to oil accumulation in hickory. Our previous preoteome results reveal the dynamic proportion of hickory caleosin (CcaCLO1) increased parallelly with embryo development [28]. Based on the transient transformation system of hickory and Arabidopsis, it is proved that CLO1 is located in the oil body. Thus, CcaCLO1 (Cca0849S0013) played a function during embryo maturation in hickory. The latest research showed that CALEOSIN 1 (CLO1), CALEOSIN 2 (CLO2), and CALEOSIN 3 (CLO3) interact with ATG8 proteins and possess putative ATG8-interacting motifs (AIMs) and then facilitate lipid droplet microautophagy in Arabidopsis seedlings [36]. However, the molecular mechanism of CcaCLO1 in regulating lipid accumulation lipolysis needs further analysis.

Caleosin has been demonstrated to exhibit a response to abiotic stresses, particularly salt stress. Kim observed a decrease in the expression of AtCLO4 in Arabidopsis following exposure to salt stress [37]. Furthermore, a study conducted by Charuchinda et al. revealed that subjecting TISTR 8580 cells from Chlorella vulgaris to 0.3 M NaCl resulted in a notable increase in triacylglycerol levels and mRNA expression levels under salt stress conditions [38]. Additionally, Xue’s research indicated that enhancing salt tolerance in Dendranthema morifolium was achieved via the overexpression of the caleosin gene BrRD20 [39]. Following the treatment of rice seedlings with NaCl, Jing et al. observed a significant increase in OsClo5 expression with prolonged treatment duration. Furthermore, overexpressing OsClo5 strains exhibited lower seed germination rate, plant stem length, and seedling weight compared to the wild type under 150 mM NaCl treatment. Conversely, OsClo5 silenced strains demonstrated higher seed germination rate, stem length, and seedling weight than the wild type [40]. In the present investigation, it was noted that the genes CcaCLO1 (Cca0849S0013) and CcaCLO2 (Cca0601S0075) exhibited an increased expression pattern in hickory callus tissue subjected to NaCl stress. These findings imply that the up-regulation of these genes in response to salt stress may be indicative of their functional significance.

Carya cathayensis, commonly referred to as Chinese hickory, is a member of the pecan genus within the Juglans family and is recognized as one of the most prolific woody oil-producing plants in China [24,41,42]. The primary focus of this research was to analyze the caleosin proteins present in hickory, with a specific emphasis on identifying the role of CcCLO1 in both oil accumulation and its response to salt stress. Via examination of the phylogenetic tree, it was observed that there are analogous proteins present in both pecan and walnut, namely Cil0992S0011 and JreChr01G11961, respectively. As a result, it can be inferred that the function of CcCLO1 closely resembles that of its homologous proteins in other species. However, due to the absence of a gene transformation system for Chinese hickory or pecan, functional validations were conducted on Arabidopsis [42,43]. However, as Arabidopsis is not an oil-producing plant, further research is needed to investigate the network mechanism of CcCLO1 in lipid accumulation, specifically in Chinese hickory or pecan. Additionally, ectopic expression of proteins in different plant species may result in mis-localization [44], necessitating enhancements to the transient transformation system for hickory. Currently, we are employing the procedures outlined in our previous study on bamboo transient transformation for the manipulation of hickory callus. However, the current efficiency levels are suboptimal and are undergoing refinement. Furthermore, future studies could involve screening for regulatory factors interacting with CcCLO1 and upstream regulatory factors using molecular biological methods to elucidate the functions of CcCLO1. Overall, elucidating the high oil content of hickory nuts lays a theoretical foundation.

Our team at Inergency is excited to announce our partnership with @Disasters Expo Europe, the leading event in disaster management.

Join us on the 15th &16th of May at the Messe Frankfurt to explore the latest solutions shaping disaster preparation, response and recovery.

To all our members, followers, and subscribers in the industry, This is YOUR opportunity to explore the intersection of sustainability and disaster management, connect with industry leaders, and stay at the forefront of emergency preparedness.

Don't miss out! Secure your complimentary tickets now and be part of the conversation driving innovation in disaster resilience. ➡️ https://lnkd.in/dNuzyQEh

And in case you missed it, here is our ultimate road trip playlist is the perfect mix of podcasts, and hidden gems that will keep you energized for the entire journey