Genomic Survey, Characterization, and Expression Profile Analysis of the SBP Genes in Pineapple (Ananas comosus L.)

Gene expression is regulated by transcription factors, which play many significant developmental processes. SQUAMOSA promoter-binding proteins (SBP) perform a variety of regulatory functions in leaf, flower, and fruit development, plant architecture, and sporogenesis. 16 SBP genes were identified in pineapple and were divided into four groups on basis of phylogenetic analysis. Five paralogs in pineapple for SBP genes were identified with Ka/Ks ratio varied from 0.20 for AcSBP14 and AcSBP15 to 0.36 for AcSBP6 and AcSBP16, respectively. 16 SBP genes were located on 12 chromosomes out of 25 pineapple chromosomes with highly conserved protein sequence structures. The isoionic points of SBP ranged from 6.05 to 9.57, while molecular weight varied from 22.7 to 121.9 kD. Expression profiles of SBP genes revealed that AcSBP7 and AcSBP15 (leaf), AcSBP13, AcSBP12, AcSBP8, AcSBP16, AcSBP9, and AcSBP11 (sepal), AcSBP6, AcSBP4, and AcSBP10 (stamen), AcSBP14, AcSBP1, and AcSBP5 (fruit) while the rest of genes showed low expression in studied tissues. Four genes, that is, AcSBP11, AcSBP6, AcSBP4, and AcSBP12, were highly expressed at 4°C, while AcSBP16 were upregulated at 45°C. RNA-Seq was validated through qRT-PCR for some genes. Salt stress-induced expression of two genes, that is, AcSBP7 and AcSBP14, while in drought stress, AcSBP12 and AcSBP15 were highly expressed. Our study lays a foundation for further gene function and expression studies of SBP genes in pineapple.


Introduction
Pineapple (Ananas comosus) is a tropical plant with edible multiple fruits, can be consumed fresh, cooked, juiced, or preserved and is the most economically significant plant in the Bromeliaceae family [1,2]. Pineapple fruits are important source of vitamins such as vitamin A and vitamin B1. In addition, it also contains a protein-digesting enzyme bromelain, a mixture of proteolytic enzymes with potential that may be used as a therapeutic agent for a variety of clinical disorders in the future [3].
Gene expression is regulated by transcription factors (TFs). Transcription factors are proteins of sequencespecific DNA binding with capability of activating and/or repressing transcription (Luscombe et al. 2000). In plants including the pineapple, transcription factors play essential roles in the regulation networks of many significant developmental processes. Genes in the SQUAMOSA promoterbinding protein (SBP-box) gene family encode SBP proteins, which perform a variety of regulatory functions that are involved in leaf development, vegetative phase change, flower and fruit development, plant architecture, sporogenesis, gibberellic acid signaling, and toxin response [4]. The SBP domain is a highly conserved DNA-binding domain and contains 76 amino acids and two zinc-binding sites, Cys-Cys-His-Cys and Cys-Cys-Cys-His [5,6]. AmSBP1 and AmSBP2 were identified as the first SBP-box genes in Antirrhinum majus discovered on the basis of their ability to bind with the promoter sequence of the floral meristem [5]. The first identified SBP-box gene in Arabidopsis was SPL3 (SQUAMOSA promoter-binding protein-like 3), which has an essential role in the floral transition [7].
Later, SBP-box genes were identified, isolated, and characterized in many plants, including Arabidopsis [6], silver birch [8], Salvia miltiorrhiza [9], rice [10], maize [11], tomato [12], grape [13], and Gossypium hirsutum [14]. Furthermore, recent studies showed that SBP-box genes show response to biotic and abiotic stresses and are also involved in signaling pathways in many species, for instance, AtSPL14 displaying sensitivity to the programmed cell death, inducing fungal toxin fumonisin B1 [15]. In grapes, VpSBP5 has the ability to show resistance to Erysiphenecator by activating the SA-induced systemic acquired resistance pathway and MeJA-induced wound signaling pathway [13]. Recently, it is considered that SBP-box genes are involved in the defense mechanism of Phytophthora blight in pepper.
The growth and development of pineapple can be affected by stress conditions such as soil type, temperature, salinity, drought, cold, and water. Importantly, gene expression is essential also for organism growth and development of the plant which could be deregulated under the stress conditions. Abiotic stresses adversely affect the growth and productivity of plants. Generally, plants are exposed to abiotic stresses like drought, salinity, low and high temperature, salt, and chemical stresses. Genetic stability and productivity of crops are affected by all these stresses. More than 50% of crop yield is reduced by abiotic stresses worldwide [16]. On the basis of genome sequence analysis, we identified 16 SBP genes in pineapple and divided it into seven groups to show the phylogenetic relationship with rice and Arabidopsis. Also, this study helps to identify the gene structure of SBP gene family, phylogenetic study, protein motifs, gene location on chromosomes, and RNA-Seq analysis for different tissues. Under four abiotic stresses (cold, drought, salt, and heat), the gene expression profiles were assessed. This study will help to identify SBP genes which are showing responses to different stress conditions and can be utilized to create abiotic stress resistant varieties through breeding programs to address the issue of yield loss due to abiotic stress. The effect of SBP proteins on gene expression is well characterized in other plants under stress conditions but is yet to be characterized in pineapple. Therefore, the present study was conducted to characterize SBP genes and the impact of SPB proteins on the gene expression in pineapple under abiotic stress conditions.

Genome-Wide Identification of SQUAMOSA (SBP) Gene
Family Homologue Sequences. To identify the SBP gene family in pineapple (Ananas comosus), the Hidden Markov Model (HMM) profile of the SBP domain (PF03110) was downloaded from the Pfam database. The SBP protein sequences of pineapple, Arabidopsis thaliana, and rice were obtained from the Phytozome, UniPort, and RGAP databases, respectively [17]. For the identification of SBP transcription factor-coding genes of Ananas comosus, the Hidden Markov Model (HMM) profile of the SBP domain was used to query pineapple (Ananas comosus) genome database as well as nonredundant protein database with the help of the BLAST program. Primers were designed to amplify the sequences from PrimerQuest tool software.

Characteristics of SBP Genes in Pineapple.
To identify the characteristics of SBP genes in pineapple, ExPASy server was used to calculate the isoionic point (IP) and molecular weights of SBP genes, while the number of amino acids, introns, exons, and open reading frame (ORF) was calculated manually.

Phylogenetic Analysis.
Multiple sequence alignments were carried out based on amino acid sequences of SBP gene family using MUSCLE software. To study the evolutionary relationship of SBP gene family in pineapple with Arabidopsis thaliana and rice, a phylogenetic tree was constructed based on sequence alignments using RAxML and the maximum likelihood method with the following parameters: JTT model, pair wise gap deletion, and 1000 bootstraps [18]. Furthermore, for tree construction maximum likelihood, minimal evolution and PhyML methods were also applied to validate the results of the neighbour-joining (NJ) method. To find out phylogenetic relationship, Ka/Ks values of SBP genes were also conducted.

Calculation of Ka/Ks
Values. Pairwise alignment of SBP gene encoding sequences of the orthologous and paralogous pairs was carried out using bio-pipeline-master software to calculate Ka/Ks value (https://github.com/tanghaibao/biopipeline/tree/master/synonymous_calculation).

SBP Gene Location on Chromosomes.
To derive the location of SBP genes on chromosomes, genes were mapped on chromosomes using Mapinspect software according to the start site (first nucleotide) of gene in pineapple. Start points and stop points of SBP genes were calculated from pineapple genome.
2.6. Gene Structure and Conserved Motif Analysis. The structural organizations of exon-intron of SBP gene family in pineapple were determined using the Gene Structure Display Server. Furthermore, by computer alignments and manual analysis, the positions of exons, SBP domains, and CDS were identified [19]. Subsequently, MEME program was employed to identify the motifs of SBP protein sequences in pineapple genome, with the motif length set to 6-100 and motif sites to 2-120. The maximum number of motifs was set to 10. The distribution of one single motif was "any number of repetitions," and the other parameter was "search-given strand only." 2.7. Plant Material and Growth Conditions. In this study, pineapple (Ananas comosus) variety MD2 was used, which are provided by the Haixia Institute of Science and Technology, Center for Genomics and Biotechnology, Fujian Agriculture and Forestry University, Fujian, China. Pineapple plants were grown in medium plastic pots containing soil and were placed in a greenhouse at 30°C with light availability of 60-70 mMol photons m −2 s −1 , under 70% humidity, with 16 hours light or an 8-hour dark photoperiod.
2.8. RNA-Seq for Specific Tissues in Pineapple. Four stages of sepal, three stages of petal, five stages of stamen, and seven stages of ovule samples were collected from pineapple MD2 variety. Samples were stored in liquid nitrogen, and RNA was extracted through Plant RNeasy Mini Kit (Qiagen, Hilden, Germany) following manufacturer protocol. The cDNA libraries were constructed using the NEBNext Ultra™ RNA Library Prep Kit for Illumina (NEB) following standard protocols. Briefly, the mRNA was fragmented and the fragmented mRNA was converted into a double-stranded cDNA. The cDNA ends were polished and ligated to NEBNext adapters. The adaptor-ligated DNA was amplified by 15 cycles of PCR followed by purification to obtain the final library for sequencing. The DNA yield and fragment insert size distribution of the library were determined on the Agilent Bioanalyzer. The libraries were sequenced on the HiSeq 2500 sequencing instrument using 100 bp pair-end protocols at the Center for Genomics and Biotechnology in Fujian Agriculture and Forestry University [20]. RNA-Seq data from the following tissues, that is, flower, leaf, root, and six stages of fruit, was extracted from a previous study [21]. A heat expression map was developed for sixteen SBP genes to show the expression pattern in different tissues in pineapple.
Additionally, the accuracy of RNA-Seq data was confirmed through qRT-PCR. Total RNA was extracted from specific tissues (leaf, root, flower, stamen, and fruit) and qRT-PCR was run as we do for abiotic stress study (shown below).
2.9. Abiotic Stress Treatment. Fully grown pineapple plants (3 biological replicates) were exposed to the four different types of abiotic stresses (cold, heat, drought, and salt). With cold stress at 4°C [22] and heat stress at 45°C, 400 mMol/mL mannitol and 400 mMol/mL NaCl were applied for drought stress and salt stress, respectively, and samples were collected at 24-hour and 48-hour time points, respectively. Samples were stored in liquid nitrogen at −80°C, RNA was extracted using RNA extraction kit (Omega Bio-Tek), cDNA was extracted using cDNA extraction kit (TransGen), and qPCR was performed using qPCR kit (TransGen), while the cis-DNA-acting elements of the pineapple SBP gene promoter were carried out to test the function of SBP genes, especially against the abiotic response [23].

Quantitative Real-Time PCR.
To investigate expression profiles of sixteen SBP genes in pineapple under abiotic stresses, pineapple (MD2) plants were exposed to four different treatments, cold, heat, drought, and mannitol for 24 h and 48 h and changes in the gene expression were analyzed using qRT-PCR. Quantitative real-time PCR was performed to determine the relative transcript levels of selected SBP family genes in pineapple according to the manufacturer's instructions for the Bio-Rad real-time PCR system (Foster City, CA, USA). PrimerQuest tool was used to design primers for SBP gene family (http://www.idtdna.com/PrimerQuest/ Home/Index?Display=AdvancedParams). The actin gene (housekeeping gene) from pineapple was used as a reference. Total RNA was extracted using RNA extraction kit (Omega Bio-Tek) following manufacturer's protocol, and cDNA was synthesized from RNA by using cDNA extraction kit (Trans-Gen) with gDNA remover. Quantitative real-time PCR (qRT-PCR) was performed using qRT-PCR kit (TransGen). Real-time PCR system in a total volume of 20 μL with the following steps: denaturation at 95°C for 30 s, annealing at 95°C for 5 s, and extension at 60°C for 30 s followed by 39 cycles. SYBR green fluorescence was measured at each cycle to monitor the amplification of target SBP family genes. The Ct (threshold cycle) values were measured at each cycle, Ct (threshold cycle) is defined as the real-time PCR cycle at which a statistically significant increase in reporter fluorescence is first detected. Three biological replicates were carried out. The relative transcript levels of the analyzed SBP family genes were normalized to the transcript levels of AcoActin gene. Livak method (−ΔΔCT) was used to calculate the expression level of SBP family genes in pineapple [24].

Identification of SQUAMOSA (SBP) Gene Family in
Pineapple. To identify the SBP gene family in pineapple (Ananas comosus), the Hidden Markov Model (HMM) profile of the SBP domain (PF03110) was downloaded from the Pfam database. For the identification of SBP transcription factor-coding genes of Ananas comosus, the HMM profile of the SBP domain was used to query pineapple (Ananas comosus) genome database as well as nonredundant protein database with the help of the BLAST program. Sixteen SBP gene sequences were identified from the HMMER database in pineapple genome.

Characterization of SBP Genes in Pineapple.
ExPASy server was used to calculate the isoionic point (IP) and molecular weights of SBP genes. The isoelectronic point or isoionic point (IP) is the pH, at which the amino acid does not migrate in an electric field. The lowest IP value was 6.05 for AcSBP5, while the highest IP value was 9.57 for AcSBP16. The molecular weight ranges from 22.7 to 121.9 kD for AcSBP8 and AcSBP2. SBP gene, that is, AcSBP2, has the highest number of amino acids (1113) and ORF (3339), followed by SBP gene, that is, AcSBP8, having the lowest number of amino acids (210) and ORF (630) ( Table 1).

Phylogenetic Analysis of SBP Gene Family in Pineapple.
Multiple sequence alignments were carried out based on full-length amino acid sequences of SBP gene family using MUSCLE software. To study the evolutionary relationship of SBP gene family in pineapple with Arabidopsis thaliana and rice, an unrooted phylogenetic tree was constructed using RAxML and the maximum likelihood method. The bootstrap test was performed with 1000 iterations based on the alignment of the SBP domain sequences. The phylogenetic analysis showed sixteen (16) SBP genes in pineapple (Ananas comosus) genome, seventeen (17) SBP genes in Arabidopsis thaliana, nineteen (19) SBP genes in Oryza sativa (Figure 1(a)). An individual phylogenetic tree was also constructed for 16 pineapple SBP gene family and divided it into 7 groups (A, B, C, D, E, F, and G) through RAxML and the maximum likelihood method [25,26] (Figure 1(b)).

Calculation of Ka/Ks
Values. The divergence and evolutionary relationship were calculated by Ka, Ks, and Ka/Ks ratios between SBP genes. We identified 5 paralogs in pineapple for SBP genes using bio-pipeline-master (https://github. com/tanghaibao/bio-pipeline/tree/master/synonymous_cal culation). Ka/Ks ratio varied from 0.20 for AcSBP14 and AcSBP15 to 0.36 for AcSBP6 and AcSBP16 ( Table 2). The trend of Ka/Ks ratio in the coding sequences of each gene pairs indicates that the SBP domains have undergone strong purifying selection (Ka/Ks < < 1) during the process of evolution. In general, Ka/Ks ratio less than 1, equal to 1, and greater than 1 means negative or stabilizing selection, neutral selection, and positive selection, respectively [27].

SBP Gene Location on Chromosomes.
To derive the SBP gene location on chromosomes, Mapinspect software was used to locate the SBP genes on chromosomes. It was found that sixteen SBP-box genes were located on twelve out of twenty-five chromosomes in pineapple genome: chromosomes 1, 2, 3, 4, 5, 6, 7, 12, 17, 19, 20, and 23. Two of the SBP genes were located on chromosome 20, two on chromosome 4, two on chromosome 3, and two on chromosome 19 ( Figure 2).

Gene Structure and Conserved Motif Analysis of SBP
Gene Family in Pineapple. Gene Structure Display Server was employed to analyze the structural organization of intron-exon of SBP gene family in pineapple. Gene structure analysis showed that most of the genes have the same number of introns. SBP gene family has 2 introns and 3 exons for most of the genes. Five SBP genes have the highest number of introns and exons, that is, AcSBP2, AcSBP1, AcSBP5, AcSBP14, and AcSBP15, indicating structural similarities (Figure 3(a)) ( Table 1). Ten protein motifs were identified through MEME program (Figure 3(b)).

Expression Profiles of SBP Genes in Different Pineapple
Tissues. To study the expression of SBP gene family in pineapple, we used these RNA-Seq data and further performed RNA-Seq analysis for four stages of sepal, three stages of petal, five stages of stamen, and seven stages of ovule [28].
RNA-Seq expression profiles for SBP gene family in pineapple revealed that there was only one gene AcSBP15 showed low expression in flower and leaf tissues, indicating SBP gene family having a little role in flower and leaf tissue development. In root tissues, two SBP genes, that is, AcSBP7 and AcSBP15, were highly expressed. Sepal is considered as the first whorl of flower, usually green in color. It provides safety and nutrients in initial stages of flower development. At sepal developmental stage 1, two genes, that is, AcSBP16 and AcSBP9, were highly expressed, while two genes, that is, AcSBP11 and AcSBP13, were moderately expressed at stage 2. High expression was shown by only one gene, that is, AcSBP11 at stage 2. On the other hand, at all four stages of sepal development, only one gene, that is, AcSBP12, was moderately expressed. Petals are collectively called as corolla. A T S P L 9 A T S P L 1 5 There was no specific SBP gene expression at the three different stages of petal development, suggesting that SBP genes may not involve in petal formation. Stamen is the male reproductive organ of flower producing pollen grains. Collectively, stamens form androecium. RNA-Seq analysis was carried out for the 5 developmental stages of stamen revealed; two genes AcSBP4 and AcSBP6 were highly expressed at the first developmental stage of stamen. One gene, that is, AcSBP10, was moderately expressed at stage 1 and stage 2.
Moderate expression of SBP gene, that is, AcSBP2, was observed at stage 3 and 4 of stamen development. Pineapple fruit is of great economical and nutritional value as being consumed as fresh, cooked, juiced, and preserved form. RNA-Seq analysis performed for 6 stages of fruit development indicated 3 genes, that is, AcSBP14, AcSBP1, and AcSBP5, were moderately expressed at stage 1, followed by gene AcSBP5 that also showed moderate expression at stage 2, while the rest of the SBP genes were not expressed at the four stages of fruit development. So, for the improvement of the pineapple crop, further studies can be done related to genes that are highly expressed in SBP genes showing high RNA-Seq expression regarding fruit development. Furthermore, ovule, plant structure that develops into a seed when fertilized. Ovule development was divided into 7 stages. RNA-Seq analysis showed no specific expression at all stages of ovule development ( Figure 4). The reads and depths obtained during RNA-Seq for each tissue are listed in additional data. The highest reads were observed for fruit developmental stage 6 (10.69 G) with the highest depths (20.32) while the lowest reads were reported for fruit developmental stage 1 (2.57 G) with sequence depths (4.88) (Supplementary Table 3s available online at https://doi.org/ 10.1155/2017/1032846). RNA-Seq data for four genes, that is, AcSBP2, AcSBP5, AcSBP7, and AcSBP15, were selected, and their RNA-Seq expression was confirmed through qRT-PCR in five different tissues (leaf, root, flower, stamen, and fruit). The relative expression of qRT-PCR for these genes showed consistent results with the RNA-Seq expression. The data of relative expression of these genes and FPKM data of RNA-Seq were shown using a heat map ( Figure 5).

Expression Pattern of SBP Genes under Cold Stress.
Pineapple is cold sensitive and cannot withstand cold and mostly all pineapple varieties are injured after exposure to 4°C for 24 hours. After exposure to cold stress at 4°C, four genes, that is, AcSBP11, AcSBP6, AcSBP4, and AcSBP12, were highly expressed, while 3 genes, that is, AcSBP16, AcSBP13, and AcSBP5 were moderately expressed followed by low expression of two genes, that is, AcSBP2 and AcSBP8 after 48 h of cold treatment. On the other hand, only one gene AcSBP3 showed high expression on 24 h as compared to 48 h treatment. The result indicates that these genes can be utilized to develop cold-resistant varieties in pineapple ( Figure 6).

Expression Pattern of SBP Genes under Heat Stress.
In this study, heat stress was applied at 45°C for 24 h and 48 h intervals, showing that after exposure to heat stress, gene expression of two genes, that is, AcSBP16 and AcSBP1, were noted to be increased at 45°C after 48 h. Moreover, the rest of all SBP genes were moderately expressed after 48 h of treatment except one gene AcSBP14, which exhibited high expression after 24 h as compared to after 48 h. Therefore, there is no significant difference in SPB family gene expression when the plants are exposed to heat stress. That reveals that SBP family genes are not significantly responsive to heat stress ( Figure 6).

Expression Pattern of SBP Genes under Salt Stress.
To evaluate the gene expression of SBP gene family under saline conditions, pineapple plants were exposed to salt stress as 400 mMol/mL NaCl for 24 h and 48 h. Results showed that two genes, that is, AcSBP7 and AcSBP14, were highly expressed, two genes, that is, AcSBP5 and AcSBP10, exhibited very low gene expression on 48 h. After 48 h of treatment, the expression pattern of four genes, that is, AcSBP2, AcSBP9, AcSBP1, and AcSBP15, was moderate while only one gene AcSBP8 was highly expressed on 24 h but reduced expression was observed on 48 h. This dynamic response of SBP genes to salt stress predicts an important role of SBP genes in response to salt stress ( Figure 6).

Expression Pattern of SBP Genes under Drought Stress.
Pineapple plants were subjected to drought stress (mannitol in concentration of 400 mMol/mL). For drought stress, two genes, that is, AcSBP12 and AcSBP15, exhibited consistent high gene expression as compared to other SBP genes which were moderately expressed followed by very low expression of one gene AcSBP16 after 48 h of treatment. On the basis of these findings, an innovative role of SBP genes can be predicted against drought stress ( Figure 6). The cis-DNA-acting elements of the pineapple SBP gene promoter were carried out to test the function of SBP genes, especially the abiotic response [23] (supplementary Table 1s).

Discussion
Pineapple is the most economically important plant in the Bromeliaceae family [1,2]. SQUAMOSA promoter-binding protein (SBP-box) gene family encodes SBP proteins which play a vital role in plant development, signaling, and defense mechanisms. Therefore, it is important to study the characterization, analysis, and gene expression profiles of SBP gene family through genome-wide data in response to abiotic stresses in pineapple. This study was focused on 16 SBP genes in pineapple genome.  Arabidopsis thaliana and rice, showing high similarity and conservation (Figure 1(a)). SBP gene family could be divided into seven groups on the basis of structural similarities (Figure 1(b)). SBP gene family could be divided into six groups in pepper on the basis of phylogenetic analysis [29]. SBP gene family was divided into four groups in melon [30]. The structural organization of intron-exon showed high variation ranged from 1 to 11 in one gene, similar as in Arabidopsis and rice genome (Table 1). Maximum number of introns reported in pepper ranged from 0 to 11 [29] and in melon, ranged from 0 to 10 [30].
To further highlight the structure of SBP gene family, this study showed ten protein motifs, suggesting that SBP gene family was highly conserved in pineapple as the same as in rice and Arabidopsis (Figure 3(c)). High conservation among protein motifs in SBP gene family was also observed in melon, pepper, and apple [4,29,30]. The resemblance and structural variation among these motifs can be further studied to provide novel information. The results showed that SBP genes were located on 12 chromosomes in pineapple genome while chromosome 4, 3, 19, and 20 contains more than one SBP genes ( Figure 2). Ka, Ks, and Ka/Ks ratios between SBP genes showed divergence and evolutionary relationship among paralogs (Table 2). In general, Ka/Ks ratio less than 1, equal to 1, and greater than 1 means negative or stabilizing selection, neutral selection, and positive selection, respectively.

RNA-Seq Analysis in Different
Tissues. RNA-Seq is used as a powerful tool to study gene expression profiles. In this study, RNA-Seq was carried out to analyze the expression pattern of sixteen SBP genes in pineapple. We used RNA-Seq data of flower, leaf, root, and six stages of fruit provided by [21] and performed RNA-Seq analysis for further tissues, that is, four stages of sepal, three stages of petal, five stages of stamen, and seven stages of ovule (Figure 4). RNA-Seq analysis revealed that only one SBP gene showed low expression in flower and leaf tissues. High expression was shown by two SBP genes in root, indicating a vital role of SBP genes in the development of root. At developmental stages of sepal S1 and S2, high and moderate expressions of SBP genes were observed. No specific expression was noted at the three stages of development in petal, showing no role of SBP gene family in petal formation. Two SBP genes AcSBP4 and AcSBP6 were highly expressed at stamen S1 stage followed by moderate expression of only one SBP gene at stamen S1 and S2 stages while at stamen S3 and S4 stage, only one gene was moderately expressed. These SBP genes can be used to improve flower reproductive development in pineapple. Four SBP genes, that is, AcSBP14, AcSBP1, AcSBP5, and AcSBP5, showed high expression at the S1 and S2 stages of fruit development while none of the genes were expressed at other four stages of fruit development. So, these SBP genes can be used to improve the pineapple varieties to have more yields because pineapple fruit is of great importance containing bromelain, a digestive enzyme. RNA-Seq analysis showed that SBP genes were not specifically expressed in the developmental stages of ovule. Therefore, it is possible to suggest that SBP gene family have no specific role in the development of ovule (Figure 4). The above finding was supported up to some extent in pepper where there one SBP gene AcSBP02 had expression in root, leaf, green fruit, and mature fruit but no expression was observed in flower in pepper [29]. Our RNA-Seq was accurate as we confirmed it through qRT-PCR. RNA-Seq expression of different genes in different tissues was confirmed through qRT-PCR; the results were consistent with the RNA-Seq data.

Response of SBP Genes to Abiotic Stress.
Recent studies suggested that SBP-box genes are involved in defense mechanism and signal transduction, showing response to biotic and abiotic stress. AtSPL14 has been involved in programmed cell death, inducing fungal toxin fumonisin B1 [15]. AtSPL genes are coexpressed with two TFs, TGA1, and WRKY65 which can be induced by pathogens and regulate the expression of genes which are considered to be involved in response to stress, such as pathogenesis-related 1 (PR-1) protein and glutathione s-transferase 6 (GST6) [31].

Conclusions
A comprehensive analysis of SBP gene family was carried out in pineapple. A total of sixteen SBP genes were identified in pineapple genome. Furthermore, these 16 SBP genes were classified according to phylogenetic tree, gene structure, conserved amino acid sequences, protein motifs, highthroughput sequencing, and RNA-Seq expression profiles for specific tissue at different developmental stages in pineapple genome. RNA-Seq was also validated through qRT-PCR. SBP genes were characterized by different traits, that is, molecular weight, isoelectric point, ORF, number of intronexon, and number of amino acids. Also, this study accompanied by the effect of SBP genes against various abiotic stresses (cold, heat, salt, and drought) and expression profiles were analyzed by qRT-PCR. cis-DNA-acting elements were also identified. The results obtained from this investigation are useful information to highlight the molecular background of SBP gene family in pineapple. The results obtained were according to our expectations and the information published related to SBP gene family studies in other species but we consider that SBP gene family response to abiotic (1) AcSBP3 (2) AcSBP11 (3) AcSBP4 (4) AcSBP12 (5) AcSBP13 (6) AcSBP7 (7) AcSBP10 (8) AcSBP6 (9) AcSBP16 (10) AcSBP9 (11) AcSBP8 (12) AcSBP2 (13) AcSBP1      qRT-PCR was used to analyze the relative expression level of each SBP gene. The expression level of pineapple actin was used as the internal control to standardize RNA samples for each reaction, and the expression at 0 h was taken as 1 (data not shown). Expression level can be understood using the given scale.
stress is a different parameter that gives a unique understanding of SBP gene family and their function in different tissues. Our study has the efficacy for understanding how SBP genes can be utilized for its gene function and pineapple crop improvement.