妃子笑荔枝果肉材料取自中国热带农业科学院南亚热带作物研究所荔枝种植园。分别取妃子笑果肉发育的3个时期（花后35 ，56，70 d）的样品，3个生物学重复，迅速置于液氮中，用于后续实验。

转录组测序委托广州基迪奥生物公司，利用Illumina HiSeqTM4000平台进行测序，对原始测序数据进行过滤、去接头和低质量值数据后，用reads组装软件Trinity 软件进行序列组装并获得Unigene 序列，共拼接组装获得61648个Unigenes。从植物转录因子数据库下载Dof转录因子序列作为探针，对本地妃子笑荔枝果肉转录组数据进行搜索，筛选出LcDof并利用SMART和Pfam 验证，删除不完整的基因。

分别采用ProtParam、Plant-mPlc Sever、MEME和SMART在线软件对Dof基因家族的相应基本理化性质、亚细胞定位、保守基序、基序功能注释进行生物信息学分析。MEME参数设置：基序最大值设为15，基序长度范围为6～100个氨基酸，其他参数为默认值。将荔枝、拟南芥和水稻中的Dof氨基酸序列通过Clustal X 1.83进行多重序列比对，进一步通过MEGA6.0邻接法（neighbor joining，NJ）将比对结果构建系统发育树。进化树参数设置：P-距离（P-distance）模型，成对删除（pairwise deletion）空位选项，Bootstrap method取值1 000。

Dof基因家族的表达数据基于课题组前期构建的“妃子笑”果肉发育RNA-seq数据库中的表达数据的rpkm值，利用Heml热图软件将Dof基因家族的表达数据均一化处理后制作热图。

以拟南芥和水稻的Dof基因家族序列为探针，采用本地blast方法对荔枝果实发育RNA-seq数据库进行比对检索，同时利用Dof为关键词在RNA-seq数据库直接搜索；将得到的数据结果进一步进行整合分析，去除冗余序列，得到20个Dof蛋白序列。随后采用SMART和Pfam在线软件对Dof蛋白序列进行保守结构域分析，最终获得19个Dof蛋白序列（表1）。Dof基因家族成员编号按照Dof序列在RNA-seq数据库中的Unigene ID号出现顺序编号（表1）。LcDof蛋白序列长度范围为157～497 aa，平均长度为330.68 aa，LcDof18蛋白序列最短，LcDof15蛋白序列最长；LcDof蛋白序列的分子量（MW）为17.70～54.35 kDa，平均分子量为35.95 kDa。LcDof等电点（pI）在4.49～9.42，19个Dof蛋白中有6个等电点小于7显酸性；13个等电点大于7显碱性；LcDof蛋白平均等电点大于7，表明LcDof为弱碱性，在碱性的亚细胞环境中发挥作用。分析LcDof家族不稳定指数发现，LcDof4/9/11/14蛋白不稳定指数<40，为稳定蛋白，其余均为不稳定蛋白。分析LcDof家族脂肪族氨基酸指数发现，Dof家族的脂肪族氨基酸指数分布在46.2～66.32范围，脂肪族氨基酸指数与蛋白的热稳定性相关，表明Dof家族蛋白质间的热稳定性存在差异。LcDof家族蛋白质疏水性指数（Grand average of hydropathicity，GRAVY）均<0，表明LcDof家族蛋白均为亲水性蛋白。亚细胞定位预测结果表明，LcDof蛋白均定位于细胞核，这与转录因子定位于细胞核结果相符。

为了进一步了解 LcDof蛋白结构特征，利用MEME在线软件分析19个LcDof蛋白的保守基序，不同保守结构域在LcDof的位置如图1所示。LcDof基因家族中含有15个保守基序，并将得到的15个保守基序进一步进行功能注释（表2）的结果表明：15个基序中基序1出现在所有的蛋白质中，为N−末端高度保守的锌指结构域（zinc-finger Dof domain，zf-Dof），基序8、15为低密度复杂区，其余的12个基序没有对应的注释，功能未知。虽然荔枝Dof基因家族成员均含有基序1，但是LcDof基因家族成员之间包含的保守基序数目及种类存在一定的差异，其中LcDof2、LcDof15和LcDof16含有的基序数目最多，均含有11个保守基序；LcDof13含有8个基序；LcDof3、LcDof6和LcDof11均含2个基序（motif 1和motif12）；LcDof8、LcDof9和LcDof17均含2个基序（motif 1和motif14）。LcDof4仅含有1个基序。在进化树中关系较近的LcDof成员间有类似的保守基序，如GroupIV中的LcDof2、LcDof15和LcDof16。LcDof中保守结构域组成相似的成员可能具有相近的基因功能。

图  1  LcDof基因家族成员保守基序分析

Figure 1.  Distribution of conserved motifs in the LcDof gene family members

通过对荔枝Dof基因家族19个家族成员进行多重序列比对，抽取保守结构域进行观察分析（图2），发现19个LcDof蛋白结构域高度保守，均包含CX₂CX₂₁CX₂C保守基序，构成了C2-C2型单锌指结构（Zinc-finger）。为进一步了解Dof基因家族在荔枝中的进化关系和生物学功能，分别以荔枝19个Dof（LcDof）、拟南芥中36个Dof（AtDof）和水稻中的30个Dof（OsDof）蛋白序列构建系统进化树 (图3)。聚类结果表明：19个荔枝Dof家族成员分别聚为4个亚家族（Group I-Group IV）个，其中GroupIV中含有LcDof基因家族成员数最多，有7个LcDof，分别为LcDof2、LcDof4、LcDof9、LcDof13、LcDof15、LcDof16、LcDof18，占基因家族总数的36.84%。其次为Group I，含有LcDof基因家族成员数最多，有6个LcDof，分别为LcDof1、LcDof3、LcDof7、LcDof11、LcDof12和LcDof15，占基因家族总数的31.58%。第三为Group II，含有5个LcDof基因家族成员，占基因家族总数的26.32%，分别为LcDof6、LcDof10、LcDof14、LcDof17、LcDof19。Group III含有LcDof基因家族成员最少，仅有1个LcDof8，占基因家族总数的5.26%。荔枝Dof家族成员与拟南芥Dof家族成员在进化上亲缘关系较近，而与水稻Dof家族成员关系较远，其中AtDof2.1和LcDof7，AtDof1.4和LcDof1，AtDof5.4和LcDof6，AtDof1.2和LcDof3为直系同源基因，推测以上荔枝中Dof与拟南芥的Dof在生物学功能上相似。LcDof5和LcDof12，LcDof10和LcDof19，LcDof4和LcDof9，LcDof2和LcDof15为旁系同源基因，推测荔枝Dof转录因子经历了基因复制事件，有可能存在功能的冗余。

图  2  LcDof基因家族成员保守结构域

Figure 2.  Conserved domain of LcDof gene family members

图  3  荔枝Dof、拟南芥Dof和水稻Dof的邻接法系统发生树

Figure 3.  Neighbor-joining phylogenetic tree of Dof in litchi, Arabidopsis and rice

为研究荔枝Dof基因家族在荔枝果实不同发育时期的表达情况，利用妃子笑荔枝（果肉不同发育时期）的RNA-Seq 转录组数据库，找到候选的19个Dof基因对应转录本的RPKM值，然后用Heml热图软件对LcDofs的RPKM值取对数值转换制作聚类热图（图4）。图4显示，在果肉发育的不同时期均检测到19 个LcDof基因的表达，但表达丰度不同；其中LcDof7、LcDof9、LcDof12、LcDof15在果肉不同发育时期表达量均较强，而LcDof3、LcDof10、LcDof16、LcDof17、LcDof19表达量较低。Group Ⅰ中LcDof3、LcDof5、LcDof12表达规律相似，Group Ⅱ中LcDof6、LcDof10、LcDof17表达规律相似，Group Ⅳ中LcDo2、LcDof9、LcDof13和LcDof15表达规律相似，推测以上基因存在相近的基因功能。

图  4  荔枝果肉不同发育阶段Dof家族基因的表达热图

Figure 4.  Heat map for expression of LcDof genes in the pulp of litchi fruit at different fruit development stages

植物转录因子以多基因家族的形式存在，在植物生长发育及其基因表达调控过程中发挥重要的作用。随着高通量测序技术的发展，目前越来越多的物种基因组公布，在组学水平上对基因家族进行鉴定和生物信息学分析，已成为基因组学的研究热点。Dof转录因子作为植物中特有的一类转录因子，在植物生长发育和调控基因表达过程中发挥着重要的作用。目前关于植物Dof基因家族在包括模式植物拟南芥和水稻在内的多个物种中进行了鉴定和分析^[5-19]。本研究参考拟南芥和水稻Dof家族基因的信息，基于荔枝果实发育转录组数据库，利用生物信息学方法，共鉴定了19个Dof基因家族成员。通过对LcDof基因家族中各蛋白的基本理化性质、保守基序、保守结构域、系统进化、基因表达进行了生物学分析。分析表明，LcDof转录因子编码的蛋白质为亲水性蛋白，定位于细胞核，这与马铃薯^[11]、甜瓜^[20]中的报道一致。通过系统进化分析，对荔枝Dof基因家族进行了详细的分类和命名。多重序列比对后利用Dof序列全长构建系统发育树，将LcDof分为4个不同的亚家族，该结果与在拟南芥、水稻、马铃薯、香蕉等物种上面的研究结果一致^[5,11,15]，拟南芥中Dof基因家族的37个成员、水稻的30个成员、马铃薯35个成员、香蕉中74个成员分为4个亚家族。进化树中聚类在1个亚家族的Dof含有的基序数目大致相似，比如亚家族4中的LcDof2、LcDof15和LcDof16，含有相同的基序，推测以上3个Dof具有相近的生物学功能。进化树中聚类在1个亚家族的Dof在荔枝果实发育中基因表达规律大致相似，例如亚家族1中LcDof3、LcDof5和LcDof12，推测以上3个Dof具有相近的基因功能。通过对荔枝Dof基因家族蛋白保守结构域分析可以看出，荔枝Dof蛋白的N−末端的DNA结合域具有高度的保守性，包含典型的保守结构域C2-C2的C2-C2型单锌指结构（Zinc-finger），且C−末端转录调控结构域复杂多变，这与前人的研究结果一致^[5]。本研究同时对荔枝果实发育转录组中Dof蛋白的基因表达特征进行了全面的分析，其中进化关系近的LcDof成员表达规律相似，这与香蕉中的MaDof的研究结果类似^[15]。