Genome-wide identification and comparative evolutionary analysis of the Dof transcription factor family in physic nut and castor bean

DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family involved in plant growth, development and stress responses. This study presents a genome-wide comparison of Dof family genes in physic nut (Jatropha curcas) and castor bean (Ricinus communis), two Euphorbiaceae plants that have not experienced any recent whole-genome duplication. A total of 25 or 24 Dof genes were identified from physic nut and castor genomes, respectively, where JcDof genes are distributed across nine out of 11 chromosomes. Phylogenetic analysis assigned these genes into nine groups representing four subfamilies, and 24 orthologous groups were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs. Conserved microsynteny was observed between physic nut and castor Dof-coding scaffolds, which allowed anchoring of 23 RcDof genes to nine physic nut chromosomes. In contrast to how no recent duplicate was present in castor, two tandem duplications and one gene loss were found in the Dof gene family of physic nut. Global transcriptome profiling revealed diverse patterns of Jc/RcDof genes over various tissues, and key Dof genes involved in flower development and stress response were also identified in physic nut. These findings provide valuable information for further studies of Dof genes in physic nut and castor.


INTRODUCTION
DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family, which is defined by the presence of the highly conserved Dof domain at the N-terminus (Yanagisawa, 2002).The Dof domain, which includes approximately 52 amino acids with a C 2 C 2 zinc finger structure, mediates DNA binding as well as protein-protein interactions (Yanagisawa & Schmidt, 1999;Rueda-Romero et al., 2012).
genes, and semi-quantitative RT-PCR analysis revealed that 18 RcDof genes are responsive to abscisic acid (ABA) and/or gibberellic acid (Jin, Chandrasekaran & Liu, 2014).However, the homology search performed in these two studies was based on proteomes derived from automatic genome annotation (Jin, Chandrasekaran & Liu, 2014;Wang et al., 2018), leading to omission of several family members as well as incorrect gene expression annotation.Moreover, species-specific evolution of the family in these two plants is yet to be resolved.
In this study, we would like to present a genome-wide identification of the complete set of Dof family genes in physic nut and castor by taking advantage of the available genome as well as various transcriptome datasets.After an expert revision of their gene structures, sequence characteristics, chromosome locations, evolutionary relationships, and expression profiles were further investigated and comprehensively compared.These findings provide valuable information for further studies of Dof genes in physic nut and castor.

Identification and manual curation of Dof family genes
The hidden Markov model (HMM)-based profile of the Dof domain (PF02701) was obtained from the Pfam database (http://pfam.sanger.ac.uk/), and was used as a query to perform the HMMER search (Version 3.0, http://hmmer.janelia.org/).The resulted protein sequences were further adopted for the TBLASTN search (Altschul et al., 1997), and the positive genomic sequences with an E-value of less than 1e-5 were predicted as described before (Zou et al., 2015a).Predicted gene models were further validated with available cDNAs, ESTs, and RNA-seq reads.The presence of the Dof domain was confirmed using SMART (http://smart.embl-heidelberg.de/).Homology search for nucleotides or ESTs and expression annotation using RNA-seq data were performed as described before (Zou et al., 2015b).

Sequence alignment and phylogenetic analysis
Multiple sequence alignments were performed using MUSCLE (Edgar, 2004) and phylogenetic trees were constructed using MEGA (version 6.0) (Tamura et al., 2013) with the maximum likelihood method (bootstrap: 1,000 replicates).The alignment of Dof domains was displayed using Boxshade (http://www.ch.embnet.org/software/BOX_form. html), and gene structures were displayed using GSDS (Hu et al., 2015) via aligning the cDNA to the gene sequence.For synteny analysis, the all-to-all BLASTP method was used to identify duplicate pairs as previously described (Zou et al., 2018), and syntenic blocks were inferred using MCScanX (Wang et al., 2012).Considering the relatively fragmented status of the current castor genome assembly, MicroSyn (Cai et al., 2011) was also used to detect microsynteny.Best reciprocal hit (BRH)-based BLAST analysis described before (Zou et al., 2017a;Zou, Xie & Yang, 2017b) was also used to define orthologs across different species, where the E-value was set to 1e-10.Systematic names were assigned based on the phylogenetic analysis for castor, or chromosome distribution for physic nut.

Protein properties and conserved motif analysis
Theoretical molecular weight (MW), isoelectric point (pI), and grand average of hydropathicity (GRAVY) were determined using ProtParam (http://web.expasy.org/protparam/).Protein subcellular localization was predicted using CELLO (http://cello.life.nctu.edu.tw/).Conserved motifs in different Dof proteins were analyzed using MEME (http://meme-suite.org/tools/meme),and the optimized parameters were: any number of repetitions; maximum number of motifs, 15; and, the optimum width of each motif, between six and 100 residues.The MAST program (http://meme-suite.org/tools/mast) was also used to search detected motifs in protein databases.

Characterization of 25 Dof family genes in physic nut
As shown in Table 1, the search of the physic nut genome resulted in 25 Dof-coding loci from 16 out of the 6,023 scaffolds.Compared with the original genome annotation (Wu et al., 2015), one more locus was identified from scaffold241 (denoted JcDof6.1,see File S1).Interestingly, this gene was also identified by Wang et al. (2018), however, JcDof7.4,another member encoded by scaffold341, was not reported in their study.

Characterization of 24 Dof family genes in castor
A total of 24 Dof-coding loci were identified from 21 castor scaffolds (Table 1).
As observed in physic nut, most of these castor scaffolds contain a single Dof gene, whereas scaffold30138, scaffold30190, and scaffold29676 harbor two instead.The expression of all these genes was supported by available ESTs and/or RNA-seq reads.Compared with the automatic genome annotation (Chan et al., 2010), nine predicted gene models were manually optimized based on read alignment (see Files S6-S14).Conserved microsynteny was observed between castor and physic nut Dof-coding scaffolds.Although the current castor genome is fragmented in 25,763 scaffolds, 23 out of the 24 RcDof genes could be anchored to nine physic nut chromosomes based on synteny analysis (Fig. 1).
To facilitate the transfer of functional information obtained in model plants, orthologs of Jc/RcDofs were further identified in Arabidopsis and rice, two well-studied species of dicotyledons or monocotyledons respectively.The detailed orthologous relationships are shown in Table 2, whereas 24 orthologous groups (OGs) were also proposed based on the definition of a certain OG with at least one member found in at least two out of four examined species.Genes whose functions have been characterized were also summarized in each OG that they belong.Nine phylogenetic groups identified above were shown to contain one to five OGs.Consistent with the phylogenetic analysis, species-specific gene loss was observed and several OGs were only found in physic nut and castor (i.e., OG-1b and OG-2c) (Table 2).

Exon-intron structures, sequence features, and conserved motifs
The exon-intron structure was analyzed based on optimized gene models and results showed that these genes contain no or a single intron (Table 1).Approximately 48.00% of JcDof genes and 43.48% of RcDof genes are intronless.Similar exon-intron structure was also observed in rice and Arabidopsis, though two-intron-containing genes (i.e., AtDof3.6 and AtDof3.7,accounting for 5.6%) were also found in Arabidopsis (Table S1).Generally, the exon-intron structure is highly conserved within a certain phylogenetic group especially within a certain OG: Groups 2-5 usually contain no intron, whereas Groups 1 and 6-9 feature one intron with the exception of OG-7c.
Deduced JcDof and RcDof proteins consist of 160-518 or 162-506 amino acid residues, respectively.The MW ranges from 18.22 to 55.72 kDa for JcDofs, or from 18.46 to 54.98 kDa for RcDofs.The predicted pI value ranges from 4.65 to 9.44 for JcDofs, or from 4.48 to 9.51 for RcDofs.Nevertheless, the average protein length, MW and pI value are similar between two species, that is, 323 residues, 34.92 kDa or 8.00 in physic nut; 334 residues, 36.23 kDa or 7.78 in castor.The GRAVY value of all Dof proteins is less than 0, indicating their hydrophilic feature.All these proteins were shown to harbor a single  Dof domain with the Cx 2 Cx 21 Cx 2 C zinc finger pattern, and were predicted to localize to the nucleus (Fig. 3).Conserved motifs were also analyzed using MEME and results were shown in Fig. 2 and Fig. S1.Among 15 motifs identified, the broadly distributed Motif 1 is characterized as the Dof domain, in contrast, little is known about other motifs: Motif 8 is present in all members of Group 1, and several members of Groups 5 and 8; Motif 10 is widely present in Group 2, and several members of Groups 3 and 4; Motif 11 is present in all members of Group 1, and several members of Group 2; Motif 15 is present in all members of Group 9, and several members of Group 7; Motifs 2, 4, 5, 7, 9, 12 is limited to Group 1; Motifs 3 and 13 are limited to Group 2; Motif 6 is limited to Group 9; and, Motif 14 is limited to OG-1d of Group 1 (Fig. 2).
Gain or loss of certain motifs was observed within a certain group even between different OGs in the same group.Group 1 usually contains Motifs 8, 11, 1, 9, 4, 5, 12, 2, and 7 in orders; however, OG-1a only harbors Motifs 8 and 1, suggesting their early divergence; moreover, OG-1d harbors Motif 14 instead of Motifs 5 and 12 found in other OGs.OG-2b of Group 2 usually contains Motifs 11, 1, 13, 3, 10 in orders, in contrast, OG-2c only harbors Motif 1 and OG-2a has Motifs 1, 3, and 10, supporting this group can be further divided into three subgroups.Interestingly, JcDof7.3 has also lost Motif 10.OG-4b of Group 4 contains Motifs 1 and 10, whereas OG-4a only harbors Motif 1, supporting this group can be subdivided into two subgroups.OG-5b of Group 5 contains Motifs 1 and 8, whereas OG-5a only harbors Motif 1, supporting this group can also be subdivided into two subgroups.Group 8 usually features a single motif; however, OG-8c gained Motif 8 (Fig. 2).

Expression profiles of JcDof genes
Despite the expression of all JcDof genes, transcriptional profiling showed that JcDof6.1 was lowly expressed in all tissues examined.Based on the FPKM value, the total transcripts of JcDof genes were shown to be most abundant in flower bud, moderate in expanded leaf and leafage, and relatively low in seed and root.Most genes were expressed in all five examined tissues, though their transcript levels are variable.According to the expression pattern, these genes can be divided into five main clusters: Cluster I is preferentially expressed in root, where JcDof2.2 and JcDof6.3 were also abundant in leafage and flower bud; Cluster II is predominantly expressed in seed as well as high abundance in leafage; Cluster III is mostly expressed in leafage, where JcDof8.2 and JcDof9.1 were also abundant in expanded leaf; Cluster IV is typically expressed in expanded leaf, where JcDof3.4 and JcDof5.2 were also abundant in root or leafage respectively; Clusters V and VI are mostly expressed in flower bud, where Cluster V members were also highly expressed in seed (Fig. 4).Several key family members were also identified in a certain tissue: JcDof7.2represents the most expressed gene in both flower bud and root; JcDof5.3represents the most expressed gene in both expanded leaf and seed; and, JcDof9.1 represents the most expressed gene in leafage.Low seed yield, which was caused by the lower number of female flowers, is a bottleneck for commercial-scale exploration of physic nut.Transcriptome profile analysis of five typical stages during floral sex differentiation revealed thousands of DEGs involved in ovule and stamen development (Hui et al., 2017).With a focus on Dof family genes, we found that the total transcripts were most abundant in IND, followed by PID2, STD1, STD2, and PID1 (Fig. 5).Compared with IND, eight genes were significantly downregulated in PID1; three genes were significantly regulated in PID2, including one upregulated and two downregulated; 10 genes were significantly downregulated in STD1; 11 genes were significantly regulated in STD2, including one upregulated and 10 downregulated.Compared with PID1, nine genes were significantly regulated in PID2, including one downregulated and eight upregulated.Compared with STD1, six genes were significantly regulated in STD2, including two upregulated and four downregulated (Fig. 5; Table S2).More remarkable, the expression pattern of JcDof9.3 suggests its key role in the later stage of both ovule and stamen development.
Anthracnose, caused by C. gloeosporioides, is an important disease causing damages on physic nut leaves, stems, and fruits (Ferreira-Pinto, Silva & Santos, 2011).To analyze the response of JcDof genes to this pathogen, two different genotypes (i.e., the susceptible RJ127 and the resistant 9-1) were artificially infected with pure isolate of C. gloeosporioides, and the induced leaves were collected at 2, 24, 96, and 144 h after infection and pooled for Illumina RNA-seq.Differential expression analysis showed that three JcDof genes were significantly regulated, that is, JcDof3.3,JcDof7.2, and JcDof9.3:JcDof3.3 was upregulated in both cultivars, whereas JcDof7.2 and JcDof9.3exhibit an opposite trend (Fig. 6), implying something different regulation mechanism of these two cultivars.

Tissue-specific expression profiles of RcDof genes
Similar to its ortholog in physic nut (i.e., JcDof6.1), the transcript level of RcDof14 was extremely low in all examined tissues, that is, leaf, male flower, endosperm I/III, endosperm V/VI, developing seed, and germinating seed.The total transcripts were most abundant in leaf (Class I); followed by male flower and developing seed (Class II, accounting for 79.3-79.4% of Class I); moderate in endosperm I/III and germinating seed (Class III, accounting for 38.2-56.0% of Class I); and, relatively low in endosperm V/VI (Class IV, accounting for 13.7% of Class I).Based on their expression patterns, these genes were grouped into seven main clusters: Cluster I is mostly expressed in germinating seed, including RcDof04 and RcDof12, which were also abundant in male flower or endosperm I/III respectively; Cluster II is preferentially expressed in developing seed; Cluster III is predominantly expressed in endosperm I/III, where RcDof19, RcDof06, and RcDof18 were also abundant in developing seed; Cluster IV is typically expressed in male flower; Cluster V is preferentially expressed in leaf and male flower; Cluster VI is preferentially expressed in leaf and developing seed; and, Cluster VII is preferentially expressed in leaf.As for a specific tissue, several key isoforms were also identified: RcDof09 represents the most and second most expressed gene in leaf or male flower, respectively; RcDof13 represents the most expressed gene in male flower; RcDof19 represents the most expressed gene in both developing seed and endosperm I/III; RcDof12 represents the most expressed gene in both germinating seed and endosperm V/VI; RcDof18 represents the second most expressed genes in germinating seed, endosperm I/III and endosperm V/VI; RcDof18 represents the second most expressed genes in developing seed (Fig. 7).Gupta et al., 2015).By contrast, research in non-model plants is still limited.In physic nut and castor, two groups have attempted to identify Dof genes from the genomic level, which resulted in 24 or 21 family members, respectively (Jin, Chandrasekaran & Liu, 2014;Wang et al., 2018).However, due to the proteomes used in their studies are not well annotated, identification of the gene family seems to be far from complete.In the present study, the original genome sequences and various transcriptome data were also added in homology search.As expected, relatively more Dof genes, that is, 25 or 24, were identified from physic nut and castor genomes, respectively (Table 1).It is worth noting that, besides one gene (JcDof6.1)escaping from the original genome annotation, the gene models of 13 other genes were manually optimized based on available mRNA, that is, four in physic nut and nine in castor (Files S2-S14).Similar results were also observed when mining other gene families, for example, aquaporin, WRKY, papain-like cysteine protease, and respiratory burst oxidase homolog (Zou et al., 2015b(Zou et al., , 2016a(Zou et al., , 2016b(Zou et al., , 2018;;Zou, Yang & Zhang, 2019a), suggesting a manual curation of computationally predicted gene structures is highly needed in these two species.As for other non-model plants whose genomes were recently released, a similar expert revision is highly recommended.
Based on phylogenetic analysis, Dof genes in Arabidopsis and rice were previously divided into four subfamilies and 10 groups named A, B1, B2, C1, C2.1, C2.2, C3, D1, D2, and D3 (Lijavetzky, Carbonero & Vicente-Carbajosa, 2003).In that study, sequences of the conserved Dof domain were used, however, this strategy usually cannot resolve the evolutionary relationships of ancient and recent duplicates.Thereby, the full-length proteins were adopted in our study.Moreover, orthologs across different species including model plants Arabidopsis and rice were identified using BRH-based BLAST analysis (as discussed below), where the orthologous relationships between physic nut and castor were also confirmed by synteny analysis.Our results showed that, except for C3 and D3 that is specific to Arabidopsis or rice, respectively, most groups can be found in both physic nut and castor (Fig. S2).Moreover, B2 is obviously split into two clusters, which is also supported by exon-intron structures and conserved motifs (Fig. 2).Thereby, we updated the subclassification as Group 1 (D1), 2 (D2), 3 (B2), 4 (A), 5 (C2.2), 6 (B2), 7 (C1), 8 (C2.1), and 9 (B1).Among them, Group 1 represents the largest group not only in physic nut and castor but also in Arabidopsis and rice (Fig. S2).
Species-specific evolution of the Dof gene family in physic nut and castor was further investigated.Consistent with their taxonomic classification standing within the same plant family, a significant level of syntenic relation and nearly one-to-one orthologous relationship were observed between physic nut and castor Dof genes.The 25 JcDof genes are distributed across nine physic nut chromosomes, and 23 out of 24 RcDof genes could be anchored to these chromosomes based on synteny analysis.To infer putative functions of Jc/RcDof genes, 24 OGs were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs.Compared with castor, two tandem duplications (JcDof7.2/JcDof7.3/JcDof7.4)and one gene loss (the ortholog of RcDof20) were found in physic nut (Table 2).The evolution pattern is something different from the WRKY transcription factor gene family, where one gene loss, one intron loss and one proximal duplication were found in castor in comparison to physic nut (Zou et al., 2016b).Interestingly, several OGs have no members in Arabidopsis and/or rice, suggesting the occurrence of species-specific gene loss (Table 2).
Putative roles of Jc/RcDof genes can be inferred from their expression profiles and their orthologs in Arabidopsis and/or rice, whose functions have been characterized (Figs.4-7; Table 2).Among nine phylogenetic groups identified, Group 1 includes the most of five OGs, where OG-1a and OG-1e are widely distributed.OG-1d, which is absent from Arabidopsis, is widely present in flowering plants including papaya (Carica papaya, i.e., evm.TU.supercontig_161.44), a species standing very close to Arabidopsis, supporting species-specific gene loss.OG-1c, which is absent from rice, may be dicot-specific, since it can be traced back to Aquilegia coerulea, a member of the basal-most eudicot clade (Sharma et al., 2014;Zou et al., 2018).OG-1b, which is absent from both Arabidopsis and rice, may appear in the last common ancestor of core eudicots along with the c event, since it is present in eudicots but not A. coerulea.Thus far, four members of OG-1a (i.e., AtDof5.5,AtDof3.3,OsDof3.1, and JcDof5.2) and two members of OG-1c (i.e., AtDof5.2 and JcDof9.1)were shown to be involved in photoperiodic control of flowering (Fornara et al., 2009;Li et al., 2009;Yang et al., 2011), where AtDof3.3 is also associated with the tolerance to drought, cold and osmotic stresses (Corrales et al., 2017).AtDof1.5, a member of OG-1e, is involved in phytochrome signaling (Park et al., 2003), and RcDof05 was shown to be upregulated by ABA in leaf (Jin, Chandrasekaran & Liu, 2014), implying their similar functions.
Group 2, which includes three OGs, may be dicot-specific, since they are widely present in eudicots but not monocots.OG-2b and OG-2c can also been found in A. coerulea, that is, Aqcoe5G171900 or Aqcoe4G144000, respectively, suggesting their early diversification.The absence of OG-2c from Arabidopsis may be due to species-specific gene loss, since it is also present in papaya (i.e., evm.TU.supercontig_1117.2).By contrast, OG-2a may be generated in core eudicots along with the c event.Compared with castor, OG-2b has been highly expanded via tandem duplication in physic nut (Table 2).Thus far, two members of OG-2b, that is, AtDof3.4 and AtDof5.8,have been characterized, which were shown to participate in regulation of cell size and number, vascular development, stomatal functioning and morphogenesis (Konishi & Yanagisawa, 2007;Skirycz et al., 2008;Negi et al., 2013).
Group 3 contains two OGs that can be traced back to Amborella trichopoda, the sister lineage to all other flowering plants (Amborella Genome Project, 2013), though species-specific loss of OG-3b was observed in Arabidopsis (Table 2).The sole gene characterized in this group, that is, AtDof4.7 (a member of OG-3a), is associated with floral organ abscission deficiency (Wei et al., 2010).
Group 4 contains two OGs that are widely present in flowering plants including A. trichopoda.However, species-specific loss of OG-4a was observed in rice as well as a high number of other monocots.For two genes characterized in this group, AtDof1.7 (a member of OG-4a) is associated with carbon skeleton production, nitrogen assimilation and reduction of glucose level (Yanagisawa et al., 2004), whereas AtDof5.4(a member of OG-4b) is involved in regulation of cell cycle progression, cell growth, root growth and callus formation (Xu et al., 2016;Ramirez-Parra et al., 2017;Rymen et al., 2017).
Group 7 contains three OGs that can be traced back to A. trichopoda.OG-7b, which is absent from rice, is widely present in monocots as well as A. trichopoda (i.e., evm_27.TU.AmTr_v1.0_scaffold00080.30).OG-7c, which is absent from Arabidopsis, is widely present in dicots as well as papaya (i.e., evm.TU.supercontig_464.1),suggesting species-specific gene loss.Among two genes characterized in this group, AtDof5.6 (a member of OG-7a) is associated with interfascicular cambium formation and vascular tissue development (Guo et al., 2009), whereas AtDof3.2(a member of OG-7b) is involved in regulation of germination in non-after-ripened seeds (Rueda-Romero et al., 2012).
Group 8 includes three OGs: OG-8c is widely present in flowering plants including A. trichopoda, whereas OG-8a and OG-8b may be dicot-specific that appeared in core eudicots along with the c event.AtDof3.7 and AtDof2.5, two members of OG-8a, were shown to control seed germination via a maternal switch (Papi et al., 2000(Papi et al., , 2002;;Gabriele et al., 2010).AtDof4.1, a member of OG-8b, is associated with selective intercellular trafficking (Chen et al., 2013).OsDof4.1, a member of OG-8c, is involved in ammonium transport and nitrogen distribution (Wu et al., 2017).
By contrast, little information is available for another two groups, that is, Group 5 (two OGs) and Group 6 (a single OG), though they are widely distributed (Table 2).

CONCLUSION
Our paper presents the genome-wide comparative evolutionary analysis of Dof family genes in physic nut and castor, two economically important Euphorbiaceae plants that have not experienced any recent WGD.The family number of 25 JcDof or 24 RcDof genes is relatively smaller than that in other terrestrial plants reported thus far.Nevertheless, the family is highly diverse and nine groups representing previously defined four subfamilies were found.Moreover, among 24 OGs identified, several OGs were absent from Arabidopsis and/or rice, suggesting species-specific evolution and their specific functions in these two and/or other species.Conserved microsynteny and nearly one-to-one orthologous relationships were observed between JcDof and RcDof genes, though two tandem duplications and one gene loss were found in the physic nut Dof gene family.Furthermore, their expression profiles and putative roles were also investigated and discussed with functional information available in Arabidopsis and/or rice.These findings will not only facilitate further studies of Dof genes in physic nut and castor, but also provide a useful reference to analyze the Dof gene family in other species.

Figure 1
Figure 1 Chromosomal locations of 25 JcDof genes and their collinear genes in castor.Shown are nine chromosomes encoding Dof genes and the chromosome serial number is indicated at the top of each chromosome.Twenty-three RcDof genes shown just behind their collinear genes in physic nut were marked in orange.Full-size  DOI: 10.7717/peerj.6354/fig-1

Figure 2
Figure 2 Structural and phylogenetic analysis of physic nut and castor Dof genes.(A) The unrooted phylogenetic tree resulting from full-length Dof proteins with MEGA6, where the distance scale denotes the number of amino acid substitutions per site.(B) The graphic representation of exonintron structures displayed using GSDS.(C) The distribution of conserved motifs among Dof proteins, where different motifs are represented by different color blocks as indicated at the bottom of the figure and the same color block in different proteins indicates a certain motif.Full-size  DOI: 10.7717/peerj.6354/fig-2

Figure 3
Figure 3 Dof domain sequence alignment of physic nut and castor Dof proteins.The four cysteine residues putatively responsible of the zinc-finger structure are indicated.Identical and similar amino acids are highlighted in black or dark gray, respectively.Full-size  DOI: 10.7717/peerj.6354/fig-3

Figure 4
Figure 4 Tissue-specific expression profiles of JcDof genes.Color scale represents FPKM normalized log 10 transformed counts where green indicates low expression and red indicates high expression.Full-size  DOI: 10.7717/peerj.6354/fig-4

Figure 5
Figure 5 Expression profiles of JcDof genes over five typical flower stages.Color scale represents FPKM normalized log 10 transformed counts where green indicates low expression and red indicates high expression.Full-size  DOI: 10.7717/peerj.6354/fig-5

Figure 7
Figure 7 Tissue-specific expression profiles of RcDof genes.Color scale represents FPKM normalized log 10 transformed counts where green indicates low expression and red indicates high expression.Full-size  DOI: 10.7717/peerj.6354/fig-7

Table 1
Dof family genes identified in this study.