Classification and expression analysis of cucumber (Cucumis sativus L.) callose synthase (CalS) family genes and localization of CsCalS4, a protein involved in pollen development

Abstract Callose plays an important role in plant development and response to multiple biotic and abiotic stresses. Callose synthase (CalS; EC 2.4.1.34), the main enzyme in callose synthesis, has been studied in a number of crops but not in cucumber (Cucumis sativus L.). In this study, 10 genes encoding putative callose synthase were identified and classified into four clades. Expression levels of the 10 CsCalS genes, which varied considerably, were detected in most analyzed organs. Transcription of CsCalS1 and CsCalS8 was obviously induced by downy or powdery mildews, which are pathogenic fungi of cucumber. CsCalS4, which was highly expressed in male flowers, was chosen for further study. Expression analysis revealed that CsCalS4 was mainly expressed in anther, and highly expressed at the stage of anther meiosis (St10), and uninuclear pollen appearance stage (St11). Histochemical and subcellular localization analyses indicated that CsCalS4 was mainly present in pollen grains and primarily localized in the plasma membrane. The results of this study provided a basic understanding of CsCalS genes and suggest the involvement of CsCalS4 in pollen development.

Callose synthase is involved in various biological processes associated with plant growth, development and stress responses. First, CalS is necessary for the normal structure and normal functioning of the material transport system. As mentioned above, 12 AtCalS genes have been identified and characterized in Arabidopsis [9]. AtCalS7 (AtGSL7) participates in sieve plate pore formation and directly affects the transport function of phloem [13,14]. In addition, almost all cells are connected by plasmodesmata, which controls the transport of nutrients, signalling molecules and so on [1]. Compared with the wild type, the Arabidopsis mutant AtCalS10 (gsl8/chorus) showed reduced callose deposition in the plasmodesmata and cell plate, indicating that AtCalS10 (GSL8) was required for cytokinesis and cell patterning [15][16][17]. Second, CalS participates in pollen development. Mutational and correlative analyses have suggested or proven that AtCalS5, AtCalS9, AtCalS10, AtCalS11 and AtCalS12 encode isoforms with important functions in pollen development [18][19][20][21]. Although all five of these genes participate in pollen development, their functions in this process are different. AtCalS5 plays a major role in callose wall synthesis during pollen formation and participates in callose barrier and pollen tube wall formation during pollen tube elongation [19,22]. AtCalS9 is required for symmetric and asymmetric mitosis during gamete formation [20,21,23]. AtCals10 is directly involved in the initial positioning of cell plates in asymmetric mitosis during pollen microgamete formation and directly affects stomatal formation in plant leaves [17,20]. AtCals11 and AtCals12 play a crucial role in callose plate synthesis during microspore tetrad formation [18]. The third, CalS plays an important role in the response to stress. Callose deposition occurs within minutes of damage, as well as during physiological or biotic stress. CalS12 (GSL5) reportedly accounts for callose synthesis in sporophytic tissues during injury or pathogen invasion [24][25][26]. The gsl5 mutant was unable to synthesize callose at papillae and made plants more resistant to pathogens. It suggested that callose deficiency in gsl5 mutants may enhance SA signalling, thereby enhancing resistance to pathogens [24].
Cucumber is one of the world's most widely cultivated and consumed vegetable crops. nevertheless, no reports have appeared concerning CalS homologs in cucurbitaceous crops. In this study, we identified and isolated CalS family members, analyzed their expression levels and predicted their functions, with a special focus on CsCalS4. our findings could aid future studies on cucumber growth, development and stress resistance.

Plant material
Seeds of wild-type cucumber (Cucumis sativus 'Xintaimici') were sown in 20 × 20 cm plastic pots containing vermiculite and then incubated at 25/18 °C (day/night) under a 16-h/8-h light/dark regime in a phytotron. The length of floral bud served as a criterion to define the developmental stages of cucumber male flowers, as described by Sun et al. [27].
Arabidopsis ecotype Col-0 was grown in potting soil in a growth chamber under short-day (8-h light/16-h dark) or long-day (16-h light/8-h dark) conditions at 22 °C and 55% relative humidity. Arabidopsis transformation was performed with Agrobacterium tumefaciens strain GV3101 [28] using the floral dip method [29].

Identification of CalS sequences
CalS family genes of cucumber were identified by BLASTP searching against the Arabidopsis genome (https://www.arabidopsis.org). Sequences of the different CsCalS genes were obtained from the Cucurbit Genomics databases (http://cucurbitgenomics.org/). Synteny analyses between cucumber and Arabidopsis CalS sequences were based on identity percentages returned using BLASTn (https://blast.ncbi.nlm.nih.gov/). The resulting 10 putative cucumber CalS genes were named CsCalS1-10.
The complete deduced amino acid sequences of CalS genes from cucumber and Arabidopsis were used to construct a phylogenetic tree. The sequences were subjected to multiple alignment in CLUSTAL W and analyzed by the neighbor-joining method in MEGA v6.0 [30].

Gene structure and transmembrane domain (TMD) analysis
Information on CsCalS genes, including locus tags and accession numbers, were acquired from the Cucurbit Genomics databases and nCBI (national Center for Biotechnology Information). The exon-intron structure was analyzed using the GSDS2.0 website (http://gsds. cbi.pku.edu.cn/). Protein sequence analysis was performed using the SMART website (http://smart. embl-heidelberg.de/), and the biochemical properties were predicted based on the amino acid sequence with the ProtParam online software (http://web.expasy. org/protparam/). Functional domain analysis was carried out using the Prosite database (https://prosite. expasy.org/). TMDs were analyzed using TMHMM (http://www.cbs.dtu.dk/services/TMHMM/).

Preparation and inoculation of fungal pathogen spore suspension from cucumber downy and powdery mildews
Fresh sporangia were collected and brushed into Petri dishes filled with sterile distilled water. After three rounds of filtration, a spore suspension was prepared. Suspended droplets were placed on a haemocytometer, and the concentration of spores in the suspension was determined under a 40× microscope. The concentrations of downy mildew and powdery mildew spores were 5 × 10 5 spores/mL and 1 × 10 5 spores/mL, respectively.
Six cucumber seedlings (at the six-leaf stage) with the same growth status were chosen as test plants, with four leaves per plant selected for inoculation and numbering. The front and back sides of each leaf were spray-inoculated with 2 mL of cucumber downy mildew/powdery mildew spore suspension at 17:00. Water was used as a control. The inoculated cucumber leaves were quickly covered with a self-sealing bag. Samples of treated and control plants were taken 6 and 24 h after inoculation. All samples were rinsed and immediately frozen in liquid nitrogen before storage at −80 °C.

RNA extraction and cDNA synthesis
Total RnA was extracted from the cucumber samples using a Column Plant RnAout kit (Tiandz, Beijing, China; http://www.tiandz.com). The extracted total RnA was treated with DnAase to remove residual DnA. Approximately 3 μg of purified RnA was used as a template for single-stranded cDnA synthesis using PowerScript reverse transcriptase (Invitrogen, Carlsbad, CA, USA).

Real-time quantitative polymerase chain reaction (qRT-PCR)
qRT-PCR was performed on an ABI 7500 system (Bio-Rad, United States) using the SYBR Green detection protocol (Takara). TUA was used as an internal control. Primers used for the qRT-PCR analysis are given in Supplemental Table S1. The mean expression level of relevant genes was calculated using the 2 −ΔΔCt method [31].

Construction of a β-glucuronidase (GUS) vector and GUS staining
To generate a CsCalS4 promoter::GUS construct, a 1.8-kb fragment of the CsCalS4 promoter (−1860 to −1 bp from the ATG start codon) was amplified from genomic DnA using primers PCsCalS4F and PCsCalS4F (Supplemental Table S1) and then inserted between the PstI and BamHI sites of a pCAM-BIA1391 vector.
For GUS staining, the transgenic Arabidopsis samples were incubated in GUS staining solution overnight at 37 °C using the method described by Jefferson et al. [32]. After staining, the samples were washed with 80% ethanol and observed.

Localization of CsCalS4::GFP in epidermal cells of onion (Allium cepa)
To study the subcellular localization of the CsCalS4 protein, the 1254-bp coding sequence of CsCalS4 (4504-5757 bp from the ATG start codon, including eight C-terminal TMDs) was amplified from cDnA using primers GCalS4F and GCalS4R (Supplemental Table S1). The resulting PCR product was then fused to the upstream region of green fluorescent protein (GFP) between the SmaI (5′-end)/BamHI (3′-end) sites of a pEZS-nL vector (http://deepgreen.stanford.edu).
Using a Bio-Rad PDS-1000/He particle delivery system, strips of onion bulb epidermis were bombarded with gold particles containing the CsCalS4::GFP fusion plasmid. The epidermal strips were placed on Murashige-Skoog medium and incubated in darkness at 22 °C for 24 h. GFP fluorescence was then observed under a confocal laser-scanning microscope (nikon C1, Tokyo, Japan) at a 488-nm excitation wavelength.

Identification and phylogenetic analysis of the cucumber CalS family
Ten candidate CsCalS genes were identified and named CsCalS1-10 ( Table 1). The CsCalS sequences were based on BLAST searches against nCBI and Cucurbit Genomics databases (http://cucurbitgenomics.org/). The deduced CsCalS polypeptides ranged from 1767 to 1951 amino acid residues. Eight of the deduced cucumber polypeptides (CsCalS1-8) had an average size of 1930 amino acid residues, whereas two polypeptides (CsCalS9-10) averaged 1768 residues, which was approximately 160 amino acids shorter than the first eight polypeptides (Table 1). Their isoelectric points (pI) were from 8.36 to 9.12, indicating that they were alkaline proteins. According to gene structural modelling, most CsCalS genes have 40-50 exons; the exceptions are CsCalS9 and CsCalS10, with only one exon each (Figure 1(A)). To analyze the structural characteristics of CsCalS genes, we examined their conserved regions and aligned their putative protein sequences. According to our results, CsCalS proteins have pfam14288 (FKS1_dom1;1,3-beta-glucan synthase subunit FKS1, domain-1) and pfam02364 (Glucan_synthase; 1,3-beta-glucan synthase component) domains. In addition, we found that CsCalS1-6, but not CsCalS7-10, have a pfam04652 (Vta1) domain.
All CsCalSs appeared to be membrane proteins with 13-18 TMDs (Figure 1(B)). To confirm the presence of the TMDs, the amino acid sequences were submitted to the TMHMM v.2.07 server. All CsCalS proteins contain multiple TMDs clustered into two regions (n-terminal and C-terminal), leaving a large hydrophilic central loop. The n-terminal TMDs are outside the membrane in CsCalS1, CsCalS2, CsCalS5 and CsCalS9, and inside the membrane in CsCalS3, CsCalS4, CsCalS6, CsCalS7, CsCalS8 and CsCalS10. The central loop contains a putative catalytic domain that can be further subdivided into two domains: a UDP-glucose binding domain and a glycosyltransferase domain (Table 1 and Figure 1(B)).
A phylogenetic tree based on the 10 CsCalS sequences and 12 AtCalS sequences was constructed using MEGA v6.0. In the tree shown in Figure 2, the cucumber sequences cluster into four main subfamilies. The first subfamily contains four cucumber sequences (CsCalS3, CsCalS4, CsCalS5 and CsCalS6) and four sequences from Arabidopsis (AtCalS3, AtCalS4, AtCalS5 and AtCalS52). The second subfamily contains two sequences each from cucumber (CsCalS9 and CsCalS10) and Arabidopsis (AtCalS11 and AtCalS12), and the third subfamily comprises two cucumber (CsCalS7 and CsCalS8) and two Arabidopsis (AtCalS9 and AtCalS10) sequences as well. The fourth subfamily includes two sequences from cucumber (CsCalS1 and CsCalS2) and four from Arabidopsis (AtCalS1, AtCalS2, AtCalS7 and AtCalS8).

Expression patterns of CsCalS genes in cucumber
The CsCalS expression patterns in various tissues and organs, including roots, stems, leaves, male flowers, female flowers and fruits, were analyzed by qRT-PCR ( Figure 3). Although the 10 CsCalS genes were expressed in most selected organs, their expression levels were similar or varied considerably. For example, CsCalS5, CsCalS6, CsCalS7 and CsCalS8 have similar expression characteristics: they were mainly expressed in stems, male and female flowers. The highest expression of CsCalS5, CsCalS6 and CsCalS7 was found in female flowers, whilst CsCalS8 was found in male flowers. CsCalS1 and CsCalS10 were more highly expressed in the stem and female flowers than in other tissues. CsCalS9 was highly expressed in male and female flowers. notably, CsCalS4 was mainly expressed in male flowers, where the expression level was nearly 20 times higher than in other parts, whereas there was hardly expression of CsCalS3 in male flowers. CsCalS2 was detected at similar expression levels in all selected organs, thus suggesting a role as a housekeeping gene. Genes with different expression patterns may have different functions. These results indicate that CsCalS genes are involved in various physiological and developmental processes.

Expression patterns of CsCalS genes in response to biotic stress
Downy mildew and powdery mildew are the two main diseases severely affecting cucumber production. Callose plays an essential role in the plant defense responses to invading pathogens. To determine whether CsCalS is involved in response to pathogen stress, wild-type cucumber plants were treated with cucumber downy mildew and powdery mildew pathomycetes Pseudoperonospora cubensis and Sphaerotheca fuliginea, respectively. Water was used as a control. CsCalS gene expressions were then examined in the leaves (Figure 4).
Under the control condition, the change in genes expression with time was not obvious (Figure 4(A)). Infection with cucumber downy mildew modulated by Pseudoperonospora cubensis caused a decrease in the transcript levels of all CsCalS genes 6 h following incubation. Except for the expression levels of CsCalS1 and CsCalS7, which showed some degree of increase or recovery, the gene expression levels were still decreased at 24 h. noteworthily, the expression level of CsCalS1 was dramatically increased at 24 h and was nearly two times that of the control (Figure 4(B)). Cucumber powdery mildew modulated by Sphaerotheca fuliginea generally caused the transcript levels of CsCalS genes to increase compared with the control as the infection time increased; the exceptions were CsCalS1, CsCalS4, CsCalS9 and CsCalS10, whose expression levels were changed slightly. In addition, the transcript levels of CsCalS3 and CsCalS8 were approximately five to six times higher than that of the control after 6 h of incubation, and CsCalS8 expression continued to increase to a level nearly 10 times higher than the control after 24 h of incubation (Figure 4(C)). These results indicate that CsCalS1 and CsCalS8 transcription was induced by cucumber downy mildew and powdery mildew disease, respectively.

Amino acid sequence structural analysis and subcellular localization of CsCalS4
According to our analysis of CsCalS gene expression patterns in cucumber, CsCalS4 is mainly expressed in  male flowers. We therefore further focussed on the expression characteristics and localization of CsCalS4. The amino acid sequence analyses indicated that CsCalS4 possesses an FKS1 domain, Vta-like domains and a large hydrophilic domain thought to provide glycosyl transferase catalytic activity at the C-terminus ( Figure 5(A)). The predicted topology of CsCalS4 includes 14 TMDs: six at the n-terminus and eight at the C-terminus. Structurally, CsCalS4 consists of a 5′ tail comprising 482 amino acids and an inner loop 747 amino acids in length ( Figure 5(B)). Comparison of the CsCalS4 polypeptide and AtCalS5 revealed that both sequences share highly conserved regions in their inner loops (Supplemental Figure S1).
In order to detect the subcellular localization of CsCalS4, the 1254-bp CsCalS4 oRF was fused with the GFP reporter gene to produce the CsCalS4-GFP fusion protein. Transient expression experiments in onion epidermal cells showed that CsCalS4-GFP was mainly located in the plasma membrane, whilst the transformed onion epidermal cells containing only GFP showed fluorescence throughout the whole structure ( Figure 5(C)). These results indicated CsCalS4 was a plasma membrane protein.

CsCalS4 expression during flower development and histochemical localization
To examine the spatial and temporal expression of CsCalS4 during flower development, we further examined the transcript levels of CsCalS4 in petals, sepals and anthers of male cucumber flowers. We observed that CsCalS4 was mainly expressed in anthers and only rarely in sepals and petals (Figure 6(A)). Cucumber anther development can be divided into 14 stages on the basis of floret length [27]. According to our results, CsCalS4 was weakly expressed in the stage of locules formation (St8) and increased during microsporocytes formation stage (St9). It was highly expressed at the stage of anther meiosis (St10) and uninuclear pollen appearance stages (St11). The transcript levels decreased at mature pollen formation and pollen release stages (St12 and St13) (Figure 6(B)). These results suggest that CsCalS4 participates in microspore formation and pollen development.
To confirm the CsCalS4 expression patterns and obtain detailed tissue-level expression profiles during plant development, we generated transgenic Arabidopsis plants expressing the GUS gene under the control of the CsCalS4 promoter. GUS activity was primarily limited to the anthers, and observation of anthers under magnification revealed that the staining was restricted to the pollen grains. GUS staining was also observed in the stigmas of pollinated flowers but not in unpollinated ones (Figure 6(C)). These observations indicate that CsCalS4 is mainly expressed in anther pollen grains.

Discussion
The CalS gene family is widely distributed in plants and participates in many significant physiological processes, but members of this family have only been identified and characterized at the genome level in a few species. In this study, we explored 10 CsCalS genes in cucumber. These values are consistent with the 2000-amino acid size of most proteins encoded by CalS genes, which are larger than those of most plant genes. The 10 genes are distributed on the same or different chromosomes, and the number of exons and introns varies greatly amongst them (Table 1 and Figure 1(A)). As reviewed by Chen and Kim [1], CalS proteins contain multiple TMDS, which cluster into two regions (n-terminal and C-terminal), leaving a large hydrophilic central loop facing the cytoplasm ( Figure  1(B)). The loop contains a putative catalytic domain, which is further subdivided into two domains: a UDP-glucose binding domain and a glycosyltransferase domain [33].
our phylogenetic analysis of Arabidopsis and cucumber CalS proteins divided the family into four main subfamilies ( Figure 2). As revealed by previous functional characterizations of CalS genes, members of the same CalS subfamily in different plant species may have similar gene functions, whereas members belonging to different subfamilies exhibit partially redundant roles. Expression analysis showed that each CsCalS is delicately regulated and displays a unique expression pattern (Figure 3). According to this study, the gene encoding CsCalS4, whose deduced polypeptide sequence is nearly 65% identical to that of the orthologous AtCalS5, and mRnA mainly expressed in male flowers, is likely involved in pollen development. CsCalS2, similar to AtCalS7 [13], has two unique TMDs in the n-terminal 550-amino acid region (Figure 1(B)). These two unique TMDs may play a role in targeting the subcellular location of CalS in the phloem sieve plate [34]. The expression of CsCalS2 was detected at similar expression levels in all organs (Figure 3). This suggests that CsCalS2 may be a housekeeping gene, or may play an important role in phloem development and transport in cucumber. It is worth noting that many CsCalSs were expressed highly in the stem, male and female flowers. Previous studies have shown that callose synthase activity was regulated in complex  ways [35]. As reviewed by Wu et al. [36], CalS complexes including UDP-glucose transferase 1 (UGT1), Rho-like GTPase (Rop) and sucrose synthase (SUS) are involved in activity regulation or targeting [4,10,15]. SUS catalytic product UDPG was used as the substrate for callose synthesis. our previous studies have found that CsSUS is mainly located in the vascular tissues involved in stress tolerance, flower and fruit development [37,38]. As protein complexes, CsCalS may also participate in the process of stress tolerance, flower and fruit development.
Callose plays an important role in plant defense response. In plants invaded by a pathogenic fungus, callose may be deposited between the plasma membrane and cell wall or around stomata or may affect the plasmodesmata transport following pathogen attack [39,40]. In Arabidopsis, AtCalS1, AtCalS5, AtCalS9, AtCalS10 and AtCalS12 participate in the defense response to pathogen infection [26]. AtCalS12 is closely associated with callose synthesis after pathogen invasion; the expression of AtCalS12 is induced by downy mildew (Hyaloperonospora), and its mutants exhibited reduced callose synthesis at the site of wounding and pathogen invasion [26,41,42]. As another example of the important role of callose synthase in the resistance to fungal pathogens, increased callose deposition can enhance the resistance of plants to powdery mildew (Golovinomyces cichoracearum) [43]. Downy mildew and powder y mildew respec tively caused by Pseudoperonospora cubensis and Sphaerotheca fuliginea are the main destructive cucumber diseases in agricultural areas with warm and humid climate conditions. In the present study, CsCalS1 was up-regulated by P. cubensis, and CsCalS3 and CsCalS8 were up-regulated by S. fuliginea. These results suggest that these CsCalS genes may be involved in defense responses against downy mildew and powdery mildew disease attack. These findings will help better understand the role of CsCalS in fungal disease resistance and in functional analysis of the CsCalS gene family.
Callose participates in the process of plant sexual reproduction. Callose not only acts as mechanical isolation and chemical barrier to the sex cells (such as macrospore mother cells and microspore mother cells) that are about to divide, but also is necessary for the formation of the outer wall of pollen and the maintenance of pollen fertility and vitality during pollen development [22,44]. Callose walls formed during microspore formation can be divided into two types: the peripheral wall surrounding the four microspores and the interstitial wall separating the four microspores in the tetrad stage [45]. In this study, we found CsCalS4, which is a plasma membrane protein ( Figure  5(C)), mRnA mainly expressed in anthers (Figure 6(A)). Further investigation revealed that CsCalS4 mRnA was expressed weakly in the stage of locules formation (St8), increased in the microsporocytes formation stage (St9), and had a high expression level at the stage of anther meiosis (St10) and uninuclear pollen appearance stage (St11) (Figure 6(B)), which is similar to that of rice OsGSL5 [12,19]. The knockout mutants (gsl5) and T-DnA insert mutants (GSL5-Ri) of rice failed to form normal pollen walls, resulting in decreased pollen viability and fertility [12] CsCalS4 shared high sequence similarity with AtCalS5 (Supplemental Figure  S1), which is also observed mainly in pollen grains. However, two studies have yielded conflicting data on the function of the AtCalS5 gene during microgametophyte development and pollen viability. In particular, Dong et al. [22] identified two T-DnA deletion mutants of cals5 and found that the cals5 mutants showed male sterility and lacked normal callose walls, which affected the exosporal pattern of microspores. In contrast, three additional T-DnA alleles of cals5 were isolated from Arabidopsis, which all altered the pattern of envelope deposition but did not affect pollen fertility. The calS5-3 mutation resulted in loss of callose walls and plugs in pollen tubes, whilst these pollen tubes were still able to fertilize, suggesting that a structured outer wall is not necessary for pollen development and fertilization ability [1,19]. Taken together, these two studies thus underscore the importance of generating and analyzing the allele family under investigation. nevertheless, plants overexpressing AtCalS5 also have microspore callose wall defects, a rough pollen outer wall structure, wrinkled anthers and partial male sterility [23]. Consequently, normal pollen development may require the expression pattern and expression level of AtCalS5 to be normal. Therefore, whether CsCalS4 was involved in the development of male flowers required further detailed analysis.
Different members of the gene family may have different functions, may participate in different stages of plant development, or may have similar functions, resulting in functional redundancy. In Arabidopsis, various AtCalS genes are also expressed in pollen grains and have been reported to affect pollen development. The interstitial wall biosynthesis depends on the combined function of AtCalS11 and AtCalS12 [18]. AtCalS9 and AtCalS10 play a role in the microspore entry into mitosis. Impaired AtCalS9 function resulted in dissymmetry of microspore division, irregular deposition of callose, and phagocytosis of germ cells by the cytoplasm of vegetative cells [16,20,21]. AtCalS10 is necessary for microspore asymmetric mitosis [19,20]. In addition, inserted mutant plants of AtCalS7 also produced abnormal pollen grains and aborted embryos [13]. In this study, except CsCalS4, all CsCalS7, CsCalS78 and CsCalS9 were widely expressed in male flowers ( Figure 3). We cannot rule out that these genes also play an important role in pollen development.

Conclusions
In the present study, we performed a genome-wide analysis of the characteristics, transmembrane domain, phylogenetic relationships, expression levels of 10 CsCalS genes in different tissues of cucumber plants and under biotic stress and the location of CsCalS4. The results indicated that CsCalS2, similar to AtCalS7, may play an important role in phloem development and transport in cucumber. CsCalS1 and CsCalS8 were up-regulated by Pseudoperonospora cubensis and Sphaerotheca fuliginea, respectively, suggested that these CsCalS genes may be involved in the defense responses against downy mildew and powdery mildew disease attack. CsCalS4, encodes a plasma membrane protein, expressed mainly in pollen grains at St10 and St11 of cucumber male flowers. This indicated that CsCalS4 participated in the crucial pollen formation and maturation stages in cucumber male flowers. In addition, CsCalS7, CsCalS8 and CsCalS9 were also expressed in male flowers. They may also perform unique or redundant functions during male flower development. our research results provide new and effective information for further understanding of the function of the CsCalS gene family in cucumber development and response to multiple biotic and abiotic stresses.

Data sharing
All data that support the findings reported in this study are available from the corresponding author upon reasonable request.