Allele specific expression of Dof genes responding to hormones and abiotic stresses in sugarcane
Authors:
Mingxing Cai aff001; Jishan Lin aff001; Zeyun Li aff001; Zhicong Lin aff002; Yaying Ma aff001; Yibin Wang aff001; Ray Ming aff001
Authors place of work:
College of Life Sciences, Center for Genomics and Biotechnology, Fujian Provincial Key Laboratory of Haixia Applied Plant Systems Biology, Fujian Agriculture and Forestry University, Fuzhou, Fujian, China
aff001; College of Crop Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian, China
aff002; Department of Plant Biology, University of Illinois at Urbana-Champaign, Urbana, IL, United States of America
aff003
Published in the journal:
PLoS ONE 15(1)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0227716
Summary
Dof transcription factors plant-specific and associates with growth and development in plants. We conducted comprehensive and systematic analyses of Dof transcription factors in sugarcane, and identified 29 SsDof transcription factors in sugarcane genome. Those SsDof genes were divided into five groups, with similar gene structures and conserved motifs within the same groups. Segmental duplications are predominant in the evolution of Dof in sugarcane. Cis-element analysis suggested that the functions of SsDofs were involved in growth and development, hormones and abiotic stresses responses in sugarcane. Expression patterns indicated that SsDof7, SsDof23 and SsDof24 had a comparatively high expression in all detected tissues, indicating these genes are crucial in sugarcane growth and development. Moreover, we examined the transcription levels of SsDofs under four plant hormone treatments, SsDof7-3 and SsDof7-4 were down-regulated after ABA treatment, while SsDof7-1 and SsDof7-2 were induced after the same treatment, indicating different alleles may play different roles in response to plant hormones. We also analyzed SsDofs’ expression profiling under four abiotic stresses, SsDof5 and SsDof28 significantly responded to these four stresses, indicating they are associate with abiotic stresses responses. Collectively, our results yielded allele specific expression of Dof genes responding to hormones and abiotic stresses in sugarcane, and their cis-elements could be crucial for sugarcane improvement.
Keywords:
Gene expression – Plant genomics – Sequence motif analysis – Transcription factors – Seedlings – Plant hormones – Arabidopsis thaliana – Sugarcane
Introduction
Dof (DNA-binding with one finger) transcription factors (TFs) are associated with growth and development in plants. A typical DNA-binding domain (C2/C2) exists in all Dof transcription factors and the C2/C2 domain is composed of about 52 amino acids. The C2/C2 domain contains a single zinc finger, which is beneficial for combining the 5 ′-(T/A)AAAG-3 ′ sequence with a conversed target DNA sequence [1]. The C-terminal of Dof transcription factors play important roles in transcription regulation, including interaction with diverse regulatory proteins [2].
The functions of Dof TFs have been identified in many plants. AtDAG1, a Dof transcription factor, was identified to be involved in light-quality response in Arabidopsis. In maize, Dof1 and Dof2 were identified to promote regulation of carbohydrate metabolism [3]. In potato, researches have been confirmed that Dof transcription factors StCDF1 was involved in the regulation of tuber development through restraining the expression of CO1/2 in potato [4]. In rice, over-expression of OsDof12 promoted early flowering [5]. In tomato, over-expression of Dof transcription factors SICDF3 promoted late flowering in transgenic Arabidopsis plants [6]. In Jatropha curcas, JcDof3 was regulated by circadian clock and identified to regulate flowering time [7]. Dof transcription factors were also identified to play important roles in plant hormonal signaling. In barley, HvDof19 was reported to repress the hydrolase gene when the barley aleurone was germinating [8]. OsDof3 was associate with gibberellin-related expression during germination in rice [9,10]. In Arabidopsis, Dof transcription factor OBP1 could regulate gene expression when responding to plant hormones, such as salicylic acid and auxin [11]. In addition, previous studies showed that Dof TFs are involved in abiotic stresses responses. In Arabidopsis, OBP1 was identified to play important roles in regulating the gene expressions responding to the signals of oxidative stresses [11]. In tomato, over-expression of Dof transcription factors SlCDF1 and SlCDF3 could influence the salt and drought responses of transgenic plants in Arabidopsis [6].
Dof transcription factors in different plant species have been studied in past years, such as Arabidopsis [12], rice [10], cucumber [13] and soybean [14]. However, information about Dof genes is lacking in sugarcane (Saccharum spp., Poaceae). Sugarcane is a major crop in producing biofuel and sugar, accounting for about 40% of ethanol production and 80% of sugar production all over the world [15]. The sugarcane (Saccharum spontaneum) genome was sequenced and genomic resources are available for detailed analysis of target genes [16]. We performed a comprehensive and systematic analysis to investigate the Dof genes in sugarcane genome and 29 SsDof genes were identified in sugarcane. These transcription factors were thoroughly analyzed on sequence phylogeny, exon and intron structure, motif patterns, chromosome location, duplication events and cis-element analysis. We examined the expression profiling of SsDofs in various developmental stages and tissues in sugarcane. We also analyzed the transcription levels of SsDofs under different treatments of abiotic stresses and plant hormones.
Materials and methods
Plant material and treatments
We used SES208 (Saccharum spontaneum, 2n = 8x = 64) as plant materials in our study. And these sugarcane plants grew in the green house at Fujian Agriculture and Forestry University.
For analyzing transcription levels of SsDof genes in different tissues and stages: root samples were obtained from root in seedling stage (45 days old), including the top of root (below the root hair, root-t), the middle of root (root-m) and the base of the root (root-b). Stem and leaf samples were from 9 months old premature internode (pre-m-stem3, pre-m-stem6 and pre-m-stem9), 12 months old mature internode (m-stem3, m-stem6 and m-stem9) and leaf (leaf-b, leaf-m and leaf-u) as previously described [17–19].
For analyzing transcription levels of sugarcane Dof genes under four plant hormones: the whole sugarcane seedlings (45 days old) were subjected to four plant hormones (ABA, GA, Auxin and Ethylene, purchased from Solarbio company), the leaves were collected at 24 h after treatments.
For analyzing transcription levels of sugarcane Dof genes in seedling stage and under four treatments by RT-qPCR: tissue samples were obtained from root, stem and leaf in seedling stage (45 days old). As for the cold and heat applications, the sugarcane seedlings were grown at 4°C and 38°C (artificial climate chamber from Yiheng company) for 4, 8, 12 and 24 h, respectively. In addition, the whole seedlings were performed with 15% PEG6000 (purchased from Takara company) and 100 mM NaCl (purchased from Takara company) for 4, 8, 12 and 24 h respectively.
Identification of Dof genes in Saccharum spontaneum
We obtained the sequences of Dof genes in Arabidopsis thaliana and Oryza sativa from Arabidopsis genome (http://www.arabidopsis.org/) and rice genome (http://rice.plantbiology.msu.edu/). Then we performed BLASTN to identify all Dof homolog hits in Saccharum spontaneum genome. We collected all non-redundant hits whose values were less than 1E-5. And we used the PFAM program (http://pfam.sanger.ac.uk/) and SMART program (http://smart.embl-heidelberg.de/) to further confirm the existence of Dof domain (PF002701).Then we used the GENSCAN program (http://genes.mit.edu/GENSCAN.html) to verify the sequences identified [20]. We used the ExPASy program (https://web.expasy.org/protparam/) to check the molecular weights (MW) and isoelectric points (PI) of all sequences.
Sequence analysis
We performed the ClustalW to investigate multiple sequence alignments of SsDof protein sequences. We checked the distribution of amino-acids of SsDof domains with WebLogo program (http://weblogo.berkeley.edu/logo.cgi). By performing GSDS program (http://gsds.cbi.pku.edu.cn) [21], we investigated exon and intron compositions of SsDof genes. We checked conserved motifs composition of sugarcane Dof proteins by MEME program (http://meme.nbcr.net/meme/intro.html) [22].
Phylogenetic analysis of SsDof genes in sugarcane
Based on multiple sequence alignments of SsDof and AtDof proteins and all sugarcane Dof genes could be divided into various groups. We performed phylogenetic analysis with MEGA5.0. Sequence of Dof proteins from Arabidopsis and sorghum were obtained from literature [23]. The phylogenetic tree image was enhanced by the Evolview online program (http://www.evolgenius.info/evolview).
Chromosomal distribution and gene duplication
The genomic and CDS sequences of SsDof genes were obtained from Saccharum spontaneum genome. We checked the gene duplications of SsDof genes by BLAST search in the genome. The chromosomal distribution of SsDof genes was generated by Circos software (http://circos.ca/).
Ka/Ks values of the sugarcane Dof genes
We investigated the nonsynonymous substitution rate (Ka) and synonymous substitution rate (Ks) with KaKs_Calculator v2.0 [24,25]. We calculated the divergence time of SsDof genes with the formula T = Ks/ (2 ×6.1 ×10−9) ×10−6 Mya [26].
Cis-element analysis
We extracted the 1.5kb upstream sequence of SsDof genes promoter. With the PlantPAN program [27] (http://plantpan.itps.ncku.edu.tw/) and PlantCARE program [28] (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/), we investigated the cis-elements of SsDof genes and collected the cis-element about growth and development, abiotic stresses and hormones responses in plant. The heatmap of cis-elements of SsDof genes was performed by TBtools software [29].
Expression profiling of sugarcane Dof genes by RNA-seq
RNA-seq was carried out using Illumina NovaSeq. We use the S. spontaneum AP85-441 genome as the reference genome to align the reads of SES208. Using Trinity software (https://github.com/trinityrnaseq/trinityrnaseq/wiki), we counted mappable reads from SES208 and normalized the FPKM values of each sample. Then, we used the TBtools software to generate the heatmap. RPKM value of SsDofs used in this study were shown in S5 Table.
Expression levels of SsDof genes based on qRT-PCR
We isolated RNA of sugarcane sample using Trizol [30] (purchased from Solarbio company). The Roche Lightcyler® 480 instrument was used to perform the quantitative RT-PCR. We selected the GAPDH (glyceraldehyde-3-phosphate dehydrogenase) gene as the internal standards for normalization [31], and each treatment was carried out with three replications. The expression levels of SsDof genes were calculated by the 2-ΔΔCt methods [32]. The primers of SsDofs performed were shown in S6 Table.
Results
Identification of SsDof genes in sugarcane
29 SsDof genes were identified in the sugarcane S. spontaneum AP85-441 genome and these SsDof genes were named as SsDof1-SsDof29. The alleles, tandem duplicates and paralogs of each SsDof are named by “-1” to “-7” with gene name (Table 1). Among these 29 SsDofs, four SsDofs have four alleles (SsDof5, SsDof6, SsDof7 and SsDof13), ten SsDofs have three alleles (SsDof1, SsDof3, SsDof11, SsDof12, SsDof20, SsDof22, SsDof25, SsDof26, SsDof27 and SsDof28), ten SsDofs have two alleles (SsDof2, SsDof4, SsDof9, SsDof14, SsDof15, SsDof16, SsDof18, SsDof19, SsDof23 and SsDof24), five SsDofs have only one alleles (SsDof8, SsDof10, SsDof17, SsDof21, SsDof29). In addition, ten SsDofs have one paralog (SsDof4, SsDof6, SsDof18, SsDof20, SsDof22, SsDof24, SsDof25, SsDof26, SsDof27 and SsDof28), SsDof8 and SsDof13 have two paralogs, SsDof1 have four paralogs. SsDof1, SsDof8 and SsDof20 have one tandem duplicate respectively (S1 Table).
The Open Reading Frame length of SsDofs ranged from 504 bp (SsDof15-2) to 2337 bp (SsDof13-4) (Tables 2 and 3). The encoding peptides of SsDofs ranged 167 to 778 amino acids. The molecular weight (Mw) of SsDofs ranged from 17096.27 Da to 86535.27 Da. The theoretical PI values of SsDofs varied from 4.74 (SsDof17) to 11.58 (SsDof26-2).
To explore the distribution of the homologous sequences at each position, we performed the multiple alignment analysis with SsDofs’ amino acid sequences. It was indicated that all SsDofs possess a representative DNA binding domain of 52 amino acids that included a single C2/C2 zinc finger structure. (Fig 1).
Multiple alignments of all sugarcane Dof domains were performed to plot the sequence logos with ClustalW. The information content of Dof domains at each position was indicated with the bits score. The four conserved cysteine residues of SsDof domains were indicated with the asterisks.
Phylogenetic relationships of Dof genes in sugarcane, sorghum and Arabidopsis
The amino acid sequences of all SsDofs with 36 AtDofs [33] and 28 SbDofs [34] were used to construct an unrooted phylogenetic tree (Fig 2 and S2 Table). Similar to earlier reports of AtDofs, the Dof proteins of three plants would be divided to five groups (group A, B, C, D and E). Group E contains the most Dof genes (53), accounting for 34.6%. Group A, B, C and D contain 26, 18, 32 and 24 Dof genes, respectively. Additionally, five SsDof genes belong to Group A (SsDof1, SsDof4, SsDof8, SsDof25, SsDof27); three SsDof genes belong to Group B (SsDof2, SsDof21, SsDof23); five SsDof genes belong to Group C (SsDof17, SsDof19, SsDof20, SsDof22, SsDof29); six SsDof genes belong to Group C (SsDof3, SsDof9, SsDof15, SsDof18, SsDof24, SsDof26); eleven SsDof genes belong to Group E (SsDof4, SsDof5, SsDof6, SsDof7, SsDof10, SsDof11, SsDof12, SsDof13, SsDof14, SsDof16 and SsDof28). Based on the phylogenetic tree, five pairs of putative orthologs from Saccharum spontaneum and Sorghum bicolor were also identified, such as SsDof29/SbDof25, SsDof17/SbDof11, SsDof15-1/SbDof13, SsDof11-2/SbDof8 and SsDof10/SbDof6.
All Dof genes were divided into five groups (A, B, C, D and E) with different color arcs. The blue squares, green circles and red asterisks represent Dof genes from Arabidopsis, sorghum, and sugarcane, respectively.
Motif composition and gene structure of sugarcane Dof gene family
We performed the MEME program to investigate the motif patterns of SsDof proteins. And 25 motifs were checked in SsDofs protein sequences (Fig 3A and 3B). Similar to the results in Arabidopsis [35], soybean [14], cucumber [13] and tomato [36], our results suggested that SsDof genes were highly conserved in sugarcane. The motif1 was the conserved Dof domain and distributed in each SsDof proteins. In addition, the motif patterns of SsDof proteins have similar compositions within the same group. For instance, in group Ⅰ, 10 motifs (1, 4, 7, 8, 10, 11, 15, 16, 21, 25) were the conserved motifs. There were 12 conserved motifs (1, 7, 8, 11, 12, 13, 14, 15, 16, 18, 21, 24) in group Ⅱ. And group Ⅲ contained the most numbers of motifs, including 17 conserved motifs, while group Ⅳ had only one conserved motif (Dof domain). These results indicated that there would be some similar functions of SsDof genes within the same group.
(A) The phylogenetic tree was based on multiple sequence alignments of SsDof proteins using MEGA 5 and divided into four groups (Ⅰ, Ⅱ, Ⅲ and Ⅳ). (B) The conserved motif of SsDof proteins. The motifs are shown in different color boxes with numbers 1–25. (C) The gene structures of SsDof genes. The green boxes indicate CDS; the yellow boxes indicate untranslated 5′- and 3′-regions; the black lines indicate introns.
To investigate the evolution of SsDof genes in sugarcane, we examined the gene structure of SsDof genes. As depicted in Fig 3A and 3C, the number of introns of SsDofs was no more than 5. Thirty-nine (43.8%) alleles and paralogs were intronless, whereas thirty-seven (41.6%) alleles and paralogs contained one intron. In addition, some SsDofs groups showed similar gene structure compositions. For instance, SsDofs in group Ⅲ had the most numbers of introns, in which SsDof13-4 had five introns. SsDofs in group Ⅳ were intronless except SsDof23 including one intron. In groupⅠ, the number of introns of SsDofs various from 0 to 3.
Chromosomal location and duplication of sugarcane Dof genes
SsDofs were unevenly distributed in 27 of the 32 chromosomes of S. spontaneum AP85-441 except chromosome 6C, 6D, 8A, 8C and 8D (Fig 4). Chromosome 1A and 1D contained eight SsDofs followed by seven SsDofs in chromosomes 1B and 1C. There was only one SsDofs in chromosomes 4A, 4D, 5A, 5B, 5C, 7D and 8B. There was no correlation between the number of SsDof genes and the length of sugarcane chromosomes.
Green lines indicate segmental duplication gene pairs, while the red lines show tandem duplication gene pairs. The chromosome numbers are shown in the center of sugarcane chromosomes.
Transposition events, tandem and segmental duplications are the primary reasons of gene family expansions [37]. Tandem duplication events happen when two or more genes duplicate within 200kb chromosome region [38], while segmental duplication events mean gene duplications happened in different chromosomes [39]. In this study, 49 pairs of duplicated genes were identified (S3 Table). Among these duplicated SsDof gene pairs, three gene pairs are tandem duplications (SsDof1-1/SsDof1-2, SsDof8-1/SsDof8-3, SsDof20-2/SsDof20-3), and the other forty-six gene pairs belong to segmental duplications.
The Ks, Ka and Ka/Ks ratio were calculated to investigate the divergence time of the duplication blocks. The duplications of SsDofs in S. spontaneum AP85-441 occurred approximately 0.21 Mya (million years ago) to 15.60 Mya with an average of 1.97 Mya (Table 4). SsDof1, SsDof2, SsDof3, SsDof11, SsDof14, SsDof20, SsDof23, SsDof27 and SsDof28 had undergone purifying selection because their Ka/Ks ratio were lower than 1, whereas SsDof8 had undergone positive selection as its Ka/Ks ratio was higher than 1. These results indicate that different SsDofs were under different selective constraints relating to their functions.
Cis-elements analysis of SsDof genes in sugarcane
We checked the cis-elements of SsDof genes and collected the cis-elements for growth and development, plant hormones and abiotic stresses responses in plants (Fig 5). For plant growth and development, the most frequent cis-elements identified were G-box and Sp1elements, which are related to light responses. The ABRE elements and TGACG motifs and CGTCA motifs were the most frequent elements for plant hormones-related cis-elements. For abiotic stress responses, ARE element included the most numbers of elements (S4 Table). Additionally, the promoter of SsDof13 contained most cis-elements of MYB binding site motifs. The promoter of SsDof20 contained most ABA responsive cis-elements, whereas SsDof3 contained most MeJA responsive elements and SsDof13 contained most anaerobic induction elements (Fig 5).
The numbers of cis-elements of SsDof genes for growth and development, phytohormone response, and stress response are indicated with different color boxes.
The cis-elements of SsDofs’ alleles distributed differently (Fig 5). The number of cis-elements for plant hormone responses distributed differently in SsDofs’ alleles. For example, the numbers of MeJA-responsive elements were different between alleles of SsDof3, SsDof4, SsDof5, SsDof6, SsDof9, SsDof13 and SsDof25. The numbers of gibberellin response elements were detected differently in alleles of SsDof3, SsDof4, SsDof5 and SsDof26. And the numbers of auxin responsive elements were also detected differently between the alleles of SsDof2, SsDof3, SsDof15 and SsDof28. In addition, the number of cis-elements for plant growth and development distributed differently in some SsDofs’ alleles, such as the numbers of light responses elements between alleles of SsDof1, SsDof3, SsDof5, SsDof13, SsDof18, SsDof20 and SsDof27, the numbers of the MYB binding site elements between alleles of SsDof3, SsDof4, SsDof5, SsDof13 and SsDof22.
Expression profiling of SsDof genes
To investigate the expression profiling of SsDofs, we examined the transcription levels of SsDofs in different tissues and stages, including root in seedling stage, stem in premature and mature stage, and leaf in mature stage (Fig 6). Among these SsDofs, SsDof1-2, SsDof26-2 and SsDof13-5 was not expressed in all samples, which may have special temporal expression patterns not examined in our libraries. And forty SsDofs (44.9%) were expressed in all samples. SsDof7, SsDof23 and SsDof24 had a high expression in all detected tissues. The expression of SsDof1-1, SsDof3-2, SsDof3-3, SsDof4-2, SsDof4-3, SsDof11-1, SsDof26-1, SsDof26-3 and SsDof26-4 were only detected in leaves, indicating that they may be involved in leaf development. Additionally, SsDof4-1, SsDof11-2, SsDof11-3, SsDof22-2 and SsDof28-3 only expressed in roots and leaves, indicating that these SsDofs may be associated with leaf and root development.
The root samples were obtained from root in seedling stage (45 days old), including the top of root (below the root hair, root-t), the middle of root (root-m) and the base of the root (root-b). The stem and leaf samples were from 9 months old premature internode (pre-m-stem3, pre-m-stem6 and pre-m-stem9), 12 months old mature internode (m-stem3, m-stem6 and m-stem9) and leaf (leaf-b, leaf-m and leaf-u).
Expression profiling of SsDofs’ alleles displayed differently. Some alleles of SsDof genes displayed similar expression profiling, such as the alleles of SsDof3, SsDof6, SsDof17, SsDof24 and SsDof26. However, the expression patterns were different for many SsDofs’ alleles. For example, SsDof1-3, SsDof1-4, SsDof1-6 and SsDof1-7 showed comparatively higher levels of expression in root, while SsDof1-1 and SsDof1-2 showed low expressions. SsDof7-3 had a high expression in all detected samples, while SsDof7-1, SsDof7-2 and SsDof7-4 had comparatively lower levels of expression in all detected tissues (Fig 6).
In order to verify the transcriptome data, we carried out the quantitative real-time PCR experiments. SsDof10, SsDof20 and SsDof23 showed comparatively higher levels of expression in root, while SsDof3, SsDof4, SsDof5, SsDof13, SsDof18, SsDof22, SsDof24 and SsDof28 showed relatively higher levels of expression in leaf (Fig 7). All the 12 SsDofs showed a very low level of expression in stem except SsDof17 and SsDof23. Our results were identical to the expression profiling of SsDof genes detected by RNA-Seq.
The tissue samples were obtained from root, stem and leaf in seedling stage (45 days old). The standard deviation was indicated with the vertical bars.
Expression profiling of SsDof genes responding to plant hormones
To investigate the expression profiling of SsDof genes responding to plant hormones, we examined their transcription levels under four plant hormones treatments (ABA, GA, Auxin and Ethylene). As shown in Fig 8, SsDof7 and SsDof10 were up-regulated under ABA treatment, but SsDof18 were down-regulated. After GA treatment, SsDof13 and SsDof24 were up-regulated, but SsDof18 were down-regulated. The expression of SsDof10, SsDof13 and SsDof24 increased after IAA treatment, but the transcription levels of SsDof8 and SsDof18 reduced. Under ET treatment, SsDof7 were up-regulated whereas SsDof9 and SsDof18 were down-regulated. Interestingly, after four plant hormones treatment, SsDof10 and SsDof13 were up-regulated whereas SsDof18 were down-regulated.
The whole sugarcane seedlings (45 days old) were subjected to four plant hormones (ABA, GA, Auxin and Ethylene) and the leaves were collected at 24 h after treatments. The heatmap was generated by TBtools with the FPKM values of each tissue sample.
In addition, some alleles of SsDof genes displayed similar expression profiling, such as SsDof2-1, SsDof2-2, SsDof3-1, SsDof3-2, SsDof3-3, SsDof11-1, SsDof11-2 and SsDof11-3. However, some of SsDofs’ alleles showed opposite expression pattern. For example, SsDof7-1 and SsDof7-2 were up-regulated after ABA treatment, while SsDof7-3 and SsDof7-4 were down-regulated.
Expression profiling of SsDof genes responding to abiotic stresses
12 SsDof members were selected from 29 sugarcane SsDof genes to investigate the expression profiling under various abiotic stresses. Then we conducted qRT-PCR experiments to observe their expression patterns after four treatments (4°C, 38°C, NaCl, PEG). As shown in Figs 9 and 10, SsDof5 and SsDof28 were obviously responding to all four treatments. The expression of SsDof5, SsDof10, SsDof18 and SsDof28 increased after these four treatments. All 12 SsDof genes were induced after cold treatment whereas SsDof3, SsDof4, SsDof5, SsDof17 and SsDof28 were induced after heat treatment. The transcription levels of SsDof4 and SsDof17 decreased after salt treatment. After different treatments, some SsDof genes showed opposite expression patterns. For example, SsDof17 was obviously up-regulated after cold and heat treatment whereas was down-regulated by salt treatment.
The tissue samples were obtained from leave in seeding stage (45 days old). The standard deviation was indicated with the vertical bars. (A). The sugarcane seedlings were performed with 4°C for 4, 8, 12 and 24 h respectively. (B). The sugarcane seedlings were performed with 38°C for 4, 8, 12 and 24 h respectively.
The tissue samples were obtained from leaf in seeding stage (45 days old). The standard deviation was indicated with the vertical bars. (A). The sugarcane seedlings were performed with 100 mM NaCl for 4, 8, 12 and 24 h respectively. (B). The sugarcane seedlings were performed with 15% PEG6000 for 4, 8, 12 and 24 h respectively.
Discussion
Gene expression profiles provide valuable clues for gene function. In Arabidopsis, AtDof5.8 was involved in processes of vascular development [40,41]. SsDof23 is orthologous to AtDof5.8 and had a high expression in roots and stem, indicating that SsDof23 may contain some similar functions in the development of sugarcane vascular tissues. Moreover, SsDof1, SsDof8, SsDof25 and SsDof27 are orthologous to AtDof5.7, which had been confirmed to control the differentiation of guard cells by controlling the transcription levels of genes [42]. Interestingly, the expression profiling of most SsDofs’ alleles displayed differently, such as the alleles of SsDof1, SsDof1-4 had a high expression in roots, whereas SsDof1-1 and SsDof1-5 had a low expression in roots. And the promoter regions of SsDof1-4, SsDof1-1 and SsDof1-5 contained different numbers of cis-elements for plant growth and development. These results indicated that allele specific expression of SsDof genes may be associated with cis-elements for plant growth and development.
There are previous studies about Dof genes in response to plant hormones. In potato, StDof genes showed either ABA-independent or ABA-dependent expression profiling [43]. In tobacco, NtBBF1was reported to facilitate the auxin-inducible gene expression [44]. In our study, SsDof10, SsDof13-1, SsDof13-2, SsDof13-3, SsDof13-6 and SsDof24-1 were up-regulated under four plant hormone treatments, whereas SsDof18-2 and SsDof18-3 were down-regulated, indicating these SsDofs may play important roles in response to phytohormones. Interestingly, SsDof7-1 and SsDof7-2 were up-regulated after ABA treatment, while SsDof7-3 and SsDof7-4 were repressed. Meanwhile, SsDof7-3 and SsDof7-4 had more abscisic acid responsive elements than SsDof7-1 or SsDof7-2 in their promoter regions. Our results suggested that allele specific expression of SsDof genes responding to hormones may be associated with cis-elements for plant hormones.
Cis-elements play critical roles in regulating phytohormones and abiotic stresses responses in plants [45,46]. The most cis-elements we have identified are those associated with light responsive, indicating light signals may play critical roles in transcriptional regulation of SsDofs in S. spontaneum AP85-441. Moreover, we also identified numbers of cis-elements about plant hormones and abiotic stresses in promoter regions of SsDofs. Meanwhile, most SsDofs were responsive to phytohormones and abiotic stresses detected by our data. These results suggested that SsDof genes may be involved in responding to phytohormones and abiotic stresses.
Dof genes have been reported to be associate with abiotic stresses responses. In Arabidopsis, the transcription levels of AtDof1.1 was up-regulated for three times under MeJA treatment, damaging the plant tissues [47]. In Chinese cabbage, many BraDof genes were induced obviously after cold, heat, salt and drought stresses. In tomato, SlCDF1-5 was obviously up-regulated after osmotic, cold, heat and salt treatments. Similar to previous researches, many SsDof genes were induced or repressed under cold, heat, salt and drought stresses, indicating that SsDof genes may be involved in responding to abiotic stresses. Interestingly, those SsDof genes induced were always detected about 4 hours after abiotic stresses treatments, indicating SsDofs’ expression increased immediately under cold and heat stresses. Under diverse treatments, some SsDof genes presented reverse expression patterns. For instance, SsDof17 was significantly induced by cold and heat treatment, whereas was repressed by salt treatment. Our study demonstrated that SsDof genes may play important roles in responding to various abiotic stresses in sugarcane.
Gene, genome, and segmental duplications are reported to be associate with genetic novelty [48–50]. The sugarcane genome was identified to undergo two WGD events after divergence from its closest relative and detailed analysis of the genome showed duplications in other gene families [51–56]. The duplications of SsDofs in sugarcane originated from approximately 0.21 Mya to 15.60 Mya, which indicated the duplications of SsDofs in sugarcane took place prior and after the divergence of sugarcane and sorghum. Moreover, we identified forty-nine pairs of duplicated SsDof gene pairs, including forty-six pairs of segmentally duplicated genes and three pairs of tandemly duplicated genes. This result suggested that segmental duplications are predominant in the evolution of SsDof in sugarcane.
Conclusions
We performed a comprehensive and systematic analysis to investigate the Dof genes in sugarcane genome and 29 SsDof genes were identified. Those SsDof genes were divided into five groups, with similar gene structures and motif patterns in the same group. Forty-nine pairs of duplicated SsDof genes were identified in sugarcane chromosomes. The duplications of SsDof genes originated from approximately 0.21 Mya to 15.60 Mya. Cis-element analysis suggested that the functions of SsDofs were involved in growth and development, hormone and abiotic stress responses in sugarcane. Expression patterns indicated that SsDof genes are crucial in sugarcane growth and development. The transcription levels of SsDofs under plant hormone treatments indicated that different alleles may play different roles in response to plant hormones. SsDofs’ expression profiling under four abiotic stresses indicated that they are involved in abiotic stress responses in sugarcane. This work provides a foundation for further functional analysis of SsDof genes in sugarcane.
Supporting information
S1 Table [xlsx]
List of genes identified in this study.
S2 Table [docx]
The full length Dof protein sequences in (..) and () used in phylogenetic tree construction.
S3 Table [xlsx]
Duplicated genes in .
S4 Table [docx]
Total -elements of genes.
S5 Table [xlsx]
RNA-seq data of genes.
S6 Table [docx]
The primers of genes in this study.
Zdroje
1. Yanagisawa S, Schmidt RJ (1999) Diversity and similarity among recognition sequences of Dof transcription factors. Plant J 17: 209–214. doi: 10.1046/j.1365-313x.1999.00363.x 10074718
2. Ma J, Li MY, Wang F, Tang J, Xiong AS (2015) Genome-wide analysis of Dof family transcription factors and their responses to abiotic stresses in Chinese cabbage. BMC Genomics 16: 33. doi: 10.1186/s12864-015-1242-9 25636232
3. Yanagisawa S (2000) Dof1 and Dof2 transcription factors are associated with expression of multiple genes involved in carbon metabolism in maize. Plant J 21: 281–288. doi: 10.1046/j.1365-313x.2000.00685.x 10758479
4. Kloosterman B, Abelenda JA, Gomez Mdel M, Oortwijn M, de Boer JM, Kowitwanich K, et al. (2013) Naturally occurring allele diversity allows potato cultivation in northern latitudes. Nature 495: 246–250. doi: 10.1038/nature11912 23467094
5. Iwamoto M, Higo K, Takano M (2009) Circadian clock- and phytochrome-regulated Dof-like gene, Rdd1, is associated with grain size in rice. Plant Cell Environ 32: 592–603. doi: 10.1111/j.1365-3040.2009.01954.x 19210638
6. Corrales AR, Nebauer SG, Carrillo L, Fernandez-Nohales P, Marques J, Renau-Morata B, et al. (2014) Characterization of tomato Cycling Dof Factors reveals conserved and new functions in the control of flowering time and abiotic stress responses. J Exp Bot 65: 995–1012. doi: 10.1093/jxb/ert451 24399177
7. Yang J, Yang MF, Zhang WP, Chen F, Shen SH (2011) A putative flowering-time-related Dof transcription factor gene, JcDof3, is controlled by the circadian clock in Jatropha curcas. Plant Sci 181: 667–674. doi: 10.1016/j.plantsci.2011.05.003 21958709
8. Moreno-Risueno MA, Diaz I, Carrillo L, Fuentes R, Carbonero P (2007) The HvDOF19 transcription factor mediates the abscisic aciddependent repression of hydrolase genes in germinating barley aleurone. Plant Journal 51: 352–365. doi: 10.1111/j.1365-313X.2007.03146.x 17565581
9. Washio K (2001) Identification of Dof proteins with implication in the gibberellin-regulated expression of a peptidase gene following the germination of rice grains. Biochim Biophys Acta 1520: 54–62. doi: 10.1016/s0167-4781(01)00251-2 11470159
10. Washio K (2003) Functional dissections between GAMYB and Dof transcription factors suggest a role for protein-protein associations in the gibberellin-mediated expression of the RAmy1A gene in the rice aleurone. Plant Physiol 133: 850–863. doi: 10.1104/pp.103.027334 14500792
11. Chen W, Chao G, Singh KB (1996) The promoter of a H2O2-inducible, Arabidopsis glutathione S-transferase gene contains closely linked OBF- and OBP1-binding sites. Plant J 10: 955–966. doi: 10.1046/j.1365-313x.1996.10060955.x 9011080
12. Papi M, Sabatini S, Bouchez D, Camilleri C, Costantino P, Vittorioso P (2000) Identification and disruption of an Arabidopsis zinc finger gene controlling seed germination. Genes Dev 14: 28–33. 10640273
13. Wen CL, Cheng Q, Zhao L, Mao A, Yang J, Yu S, et al. (2016) Identification and characterisation of Dof transcription factors in the cucumber genome. Sci Rep 6: 23072. doi: 10.1038/srep23072 26979661
14. Guo Y, Qiu LJ (2013) Genome-Wide Analysis of the Dof Transcription Factor Gene Family Reveals Soybean-Specific Duplicable and Functional Characteristics (Retracted article. See vol. 11, e0167019, 2016). Plos One 8.
15. Zhang J, Zhou M, Walsh J, Lin Z, Ming R (2013) Sugarcane Genetics and Genomics: John Wiley & Sons Ltd. doi: 10.1186/1471-2164-14-314
16. Zhang J, Zhang X, Tang H, Zhang Q, Hua X, Ma X, et al. (2018) Allele-defined genome of the autopolyploid sugarcane Saccharum spontaneum L. Nat Genet 50: 1565–1573. doi: 10.1038/s41588-018-0237-2 30297971
17. Chen Y, Zhang Q, Hu W, Zhang X, Wang L, Hua X, et al. (2017) Evolution and expression of the fructokinase gene family in Saccharum. BMC Genomics 18: 197. doi: 10.1186/s12864-017-3535-7 28222695
18. Hu W, Hua X, Zhang Q, Wang J, Shen Q, Zhang X, et al. (2018) New insights into the evolution and functional divergence of the SWEET family in Saccharum based on comparative genomics. BMC Plant Biol 18: 270. doi: 10.1186/s12870-018-1495-y 30404601
19. Zhang Q, Hu W, Zhu F, Wang L, Yu Q, Ming R, et al. (2016) Structure, phylogeny, allelic haplotypes and expression of sucrose transporter gene families in Saccharum. BMC Genomics 17: 88. doi: 10.1186/s12864-016-2419-6 26830680
20. Stormo GD (2000) Gene-finding approaches for eukaryotes. Genome Res 10: 394–397. doi: 10.1101/gr.10.4.394 10779479
21. Guo AY, Zhu QH, Chen X, Luo JC (2007) [GSDS: a gene structure display server]. Yi Chuan 29: 1023–1026.
22. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. (2009) MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res 37: W202–208. doi: 10.1093/nar/gkp335 19458158
23. Wu Z, Cheng J, Cui J, Xu X, Liang G, Luo X, et al. (2016) Genome-Wide Identification and Expression Profile of Dof Transcription Factor Gene Family in Pepper (Capsicum annuum L.). Front Plant Sci 7: 574. doi: 10.3389/fpls.2016.00574 27200047
24. Wang D, Zhang Y, Zhang Z, Zhu J, Yu J (2010) KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics 8: 77–80. doi: 10.1016/S1672-0229(10)60008-3 20451164
25. Posada D (2003) Using MODELTEST and PAUP* to select a model of nucleotide substitution. Curr Protoc Bioinformatics Chapter 6: Unit 6 5.
26. Lynch M, Conery JS (2000) The evolutionary fate and consequences of duplicate genes. Science 290: 1151–1155. doi: 10.1126/science.290.5494.1151 11073452
27. Chow CN, Lee TY, Hung YC, Li GZ, Tseng KC, Liu YH, et al. (2019) PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants. Nucleic Acids Res 47: D1155–D1163. doi: 10.1093/nar/gky1081 30395277
28. Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Van de Peer Y, et al. (2002) PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res 30: 325–327. doi: 10.1093/nar/30.1.325 11752327
29. Chen C, Chen H, He Y, Xia R (2018) TBtools, a <span class = "underline">Toolkit for <span class = "underline">Biologists integrating various biological data handling <span class = "underline">tools with a user-friendly interface. bioRxiv.
30. Ma J, He Y, Wu C, Liu H, Hu Z, Sun G Cloning and Molecular Characterization of a SERK Gene Transcriptionally Induced During Somatic Embryogenesis inAnanas comosus cv. Shenwan. Plant Molecular Biology Reporter 30: 195–203.
31. Ling H, Wu Q, Guo J, Xu L, Que Y (2014) Comprehensive selection of reference genes for gene expression normalization in sugarcane by real time quantitative rt-PCR. PLoS One 9: e97469. doi: 10.1371/journal.pone.0097469 24823940
32. Pfaffl MW (2001) A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res 29: e45. doi: 10.1093/nar/29.9.e45 11328886
33. Yanagisawa S (2002) The Dof family of plant transcription factors. Trends Plant Sci 7: 555–560. doi: 10.1016/s1360-1385(02)02362-2 12475498
34. Kushwaha H, Gupta S, Singh VK, Rastogi S, Yadav D (2011) Genome wide identification of Dof transcription factor gene family in sorghum and its comparative phylogenetic analysis with rice and Arabidopsis. Mol Biol Rep 38: 5037–5053. doi: 10.1007/s11033-010-0650-9 21161392
35. Lijavetzky D, Carbonero P, Vicente-Carbajosa J (2003) Genome-wide comparative phylogenetic analysis of the rice and Arabidopsis Dof gene families. BMC Evol Biol 3: 17. doi: 10.1186/1471-2148-3-17 12877745
36. Cai X, Zhang Y, Zhang C, Zhang T, Hu T, Ye J, et al. (2013) Genome-wide analysis of plant-specific Dof transcription factor family in tomato. J Integr Plant Biol 55: 552–566. doi: 10.1111/jipb.12043 23462305
37. Cannon SB, Mitra A, Baumgarten A, Young ND, May G (2004) The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol 4: 10. doi: 10.1186/1471-2229-4-10 15171794
38. Holub EB (2001) The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet 2: 516–527. doi: 10.1038/35080508 11433358
39. Liu Y, Jiang H, Chen W, Qian Y, Ma Q, Cheng B, et al. Genome-wide analysis of the auxin response factor (ARF) gene family in maize (Zea mays). Plant Growth Regulation 63: 225–234.
40. Konishi M, Yanagisawa S (2007) Sequential activation of two Dof transcription factor gene promoters during vascular development in Arabidopsis thaliana. Plant Physiol Biochem 45: 623–629. doi: 10.1016/j.plaphy.2007.05.001 17583520
41. Gardiner J, Sherr I, Scarpella E (2010) Expression of DOF genes identifies early stages of vascular development in Arabidopsis leaves. Int J Dev Biol 54: 1389–1396. doi: 10.1387/ijdb.093006jg 20563990
42. Negi J, Moriwaki K, Konishi M, Yokoyama R, Nakano T, Kusumi K, et al. (2013) A Dof transcription factor, SCAP1, is essential for the development of functional stomata in Arabidopsis. Curr Biol 23: 479–484. doi: 10.1016/j.cub.2013.02.001 23453954
43. Venkatesh J, Park SW (2015) Genome-wide analysis and expression profiling of DNA-binding with one zinc finger (Dof) transcription factor family in potato. Plant Physiol Biochem 94: 73–85. doi: 10.1016/j.plaphy.2015.05.010 26046625
44. Baumann K, De Paolis A, Costantino P, Gualberti G (1999) The DNA binding site of the Dof protein NtBBF1 is essential for tissue-specific and auxin-regulated expression of the rolB oncogene in plants. Plant Cell 11: 323–334. doi: 10.1105/tpc.11.3.323 10072394
45. Le DT, Nishiyama R, Watanabe Y, Vankova R, Tanaka M, Seki M, et al. (2012) Identification and expression analysis of cytokinin metabolic genes in soybean under normal and drought conditions in relation to cytokinin levels. PLoS One 7: e42411. doi: 10.1371/journal.pone.0042411 22900018
46. Walther D, Brunnemann R, Selbig J (2007) The regulatory code for transcriptional response diversity and its relation to genome structural properties in A. thaliana. PLoS Genet 3: e11. doi: 10.1371/journal.pgen.0030011 17291162
47. Skirycz A, Reichelt M, Burow M, Birkemeyer C, Rolcik J, Kopka J, et al. (2006) DOF transcription factor AtDof1.1 (OBP2) is part of a regulatory network controlling glucosinolate biosynthesis in Arabidopsis. Plant J 47: 10–24. doi: 10.1111/j.1365-313X.2006.02767.x 16740150
48. Taylor JS, Raes J (2004) Duplication and divergence: the evolution of new genes and old ideas. Annu Rev Genet 38: 615–643. doi: 10.1146/annurev.genet.38.072902.092831 15568988
49. Moore RC, Purugganan MD (2005) The evolutionary dynamics of plant duplicate genes. Curr Opin Plant Biol 8: 122–128. doi: 10.1016/j.pbi.2004.12.001 15752990
50. Semon M, Wolfe KH (2007) Consequences of genome duplication. Curr Opin Genet Dev 17: 505–512. doi: 10.1016/j.gde.2007.09.007 18006297
51. Wang J, Li Y, Zhu F, Ming R, Chen L-Q Genome-Wide Analysis of Nitrate Transporter (NRT/NPF) Family in Sugarcane Saccharum spontaneum L. Tropical Plant Biology.
52. Nelson DR (2019) Cytochrome P450s in the sugarcane Saccharum spontaneum. Tropical Plant Biology 12: 150–157.
53. Ma PP, Yuan Y, Shen QC, Jiang Q, Hua XT, Zhang Q, et al. (2019) Evolution and Expression Analysis of Starch Synthase Gene Families in Saccharum spontaneum. Tropical Plant Biology 12: 158–173.
54. Shi Y, Xu H, Shen Q, Lin J, Zhang J (2019) Comparative Analysis of SUS Gene Family between Saccharum officinarum and Saccharum spontaneum. Tropical Plant Biology.
55. Lin JS, Zhu MT, Cai MX, Zhang WP, Fatima M, Jia HF, et al. (2019) Identification and Expression Analysis of TCP Genes in Saccharum spontaneum L. Tropical Plant Biology 12: 206–218.
56. Zhang WP, Lin JS, Dong F, Ma Q, Wu SL, Ma XY, et al. (2019) Genomic and Allelic Analyses of Laccase Genes in Sugarcane (Saccharum spontaneum L.). Tropical Plant Biology 12: 219–229.
Článek vyšel v časopise
PLOS One
2020 Číslo 1
- Zvoní budík? Přispěte si ještě chvilku, bude vám to myslet rychleji!
- Metamizol jako analgetikum první volby: kdy, pro koho, jak a proč?
- Není statin jako statin aneb praktický přehled rozdílů jednotlivých molekul
- Vědci vytvořili první funkční mozkovou tkáň pomocí 3D tisku
- Umíme správně řešit osteporotické zlomeniny?
Nejčtenější v tomto čísle
- Severity of misophonia symptoms is associated with worse cognitive control when exposed to misophonia trigger sounds
- Chemical analysis of snus products from the United States and northern Europe
- Calcium dobesilate reduces VEGF signaling by interfering with heparan sulfate binding site and protects from vascular complications in diabetic mice
- Effect of Lactobacillus acidophilus D2/CSL (CECT 4529) supplementation in drinking water on chicken crop and caeca microbiome