基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析

热带作物学报2022, 43(5): 882 892 Chinese Journal of Tropical Crops
收稿日期 2021-10-27；修回日期 2021-12-31
基金项目 海南省自然科学基金项目（No. 321MS077）；广东省重点研发计划项目（No. 2018B020202011）；国家荔枝龙眼
产业技术体系项目（No. CARS-32-20）。

作者简介 董 晨（1981—），女，硕士，副研究员，研究方向：热带果树生理生化及分子生物学。

*通信作者（Corresponding
author ）：李伟才（LI Weicai ），E-mail ：***************。

基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析
董 晨，李金枝，郑雪文，王 弋，李伟才*
中国热带农业科学院南亚热带作物研究所/农业农村部热带果树生物学重点实验室，广东湛江 524091
摘 要：泛素结合酶在泛素化级联反应中居于中心位置，是连接泛素活化酶和泛素连接酶的桥梁，在泛素化系统中发挥着关键作用。

本研究基于荔枝转录组数据库，利用生物信息学方法对转录组数据库泛素结合酶UBC 基因家族进行鉴定，最终得到28个LcUBC 基因家族成员。

通过荧光定量PCR 技术对LcUBC 基因家族在‘妃子笑’荔枝不同组织中进行时空表达分析，并对花穗对照和烯效唑处理花穗发育的不同时期进行表达分析。

结果表明：LcUBC 基因家族成员对应所编码的氨基酸数分布在85~1146 aa 之间，等电点大小在4.03~9.61之间，5个LcUBC 蛋白为稳定蛋白，其余均为不稳定蛋白。

2个LcUBC 蛋白定位于细胞质，2个LcUBC 蛋白定位于细胞质和细胞核，1个LcUBC 蛋白定位于内质网，其余蛋白均定位于细胞核。

系统进化分析结果表明，LcUBC 蛋白分为10个亚家族；28个LcUBC 基因家族成员在不同组织的表达量存在差异；烯效唑处理‘妃子笑’荔枝花穗与对照相比，烯效唑处理21 d 后，LcUBC 基因家族成员的多个基因表达量较高。

这是在转录组水平分析LcUBC 基因家族成员，本研究结果对今后该基因家族的分类、克隆和功能研究提供了参考依据。

关键词：荔枝；泛素结合酶；基因家族；表达分析 中图分类号：S667.1 文献标识码：A
Transcriptome-wide Identification and Analysis of UBC Gene Family
in Litchi chinensis Sonn.
DONG Chen, LI Jinzhi, ZHENG Xuewen, WANG Yi, LI Weicai *
Institute of South Subtropical Corp Research, Chinese Academy of Tropical Agricultural Science / Key Laboratory of Tropical Fruit
Biology, Ministry of Agriculture and Rural Affairs, Zhanjiang, Guangdong 524091, China
Abstract: Ubiquitin-conjugating enzymes play a central role in the ubiquitination cascade. They are the bridge between ubiquitin activating enzymes and ubiquitin ligases, and play a key role in the ubiquitination system. Through the local BLAST search and keyword search of transcriptome database, 32 ubiquitin-conjugating enzyme gene family members, including 5 UEV and 28 UBC genes, were obtained through smart screening of conserved domains. Only LcUBC1– LcUBC28 was studied in the article. Based on the litchi transcriptome database, the ubiquitin-conjugating enzyme gene family was identified by bioinformatics methods, and the temporal and spatial expression of ubiquitin-conjugating en-zyme gene family in different tissues and organs of Feizixiao litchi was analyzed by the fluorescence quantitative PCR technology. Meanwhile, the expression of the control and uniconazole treatment at different stages of flower develop-ment was analyzed. The results showed that the amino acids were 85–1146 aa. The isoelectric point ranged from 4.03 to 9.61. Five LcUBC proteins were stable, and the others LcUBC proteins were unstable. Two LcUBC proteins were lo-cated in cytoplasm, two LcUBC proteins in cytoplasm and nucleus, one LcUBC protein in endoplasmic reticulum and the rest LcUBC proteins in nucleus. There were great differences in the length of the protein and the changes in the characteristics of the protein, indicating that the proteins of the ubiquitin protease gene family had different characteris-tics. In the multiple sequence alignment, all 28 LcUBC genes contained an evolutionarily highly conserved amino acid
All Rights Reserved.
第5期董晨等：基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析 883
residue cysteine active site, and the C-terminal tryptophan (W) located behind the cysteine active site. The 28 LcUBC
proteins were also highly conserved, and the “HPNI” conserved motif was located 7–10 amino acids before the active
site of cysteine. The LcUBC gene family was predicted to contain 10 conserved motifs. There were certain differences
in the number and types of conserved motifs contained in the members of the LcUBC gene family. Among them, the
conserved motifs 1 and 2 were very important conserved motifs in the LcUBC domain. Phylogenetic analysis showed
that the litchi UBC protein was divided into 10 subfamilies, namely UBC1, UBC2/11, UBC3/7, UBC4/5, UBC6, UBC8,
UBC10/12, UBC13, UBC14, UBC15. The 28 LcUBC genes were expressed in at least one tissue, and the expression
levels of different genes in different tissues were different. Among them, the expression of LcUBC4, LcUBC5, and
LcUBC18 in seeds were significantly higher than that of other tissues, and the expression of LcUBC7 in male flowers
was significantly higher than that of other tissues. The expression level of female flowers of LcUBC14 and LcUBC15
was significantly higher than that of other tissues, and LcUBC10 was highly expressed in pericarp and old leaves.
Compared with the control, the expression levels of LcUBC genes in different developmental stages of Feizixiao litchi
flower spikes treated with uniconazole were higher at 21 days after the treatment with Uniconazole. This is the first time
that the litchi UBC gene family was analyzed at the transcriptome level. The results of the study would provide some
referrences for future studies on the classification, cloning and function of this gene family.
Keywords: Litchi chinensis Sonn.; ubiquitin-conjugating enzyme; genes family; expression analysis
DOI: 10.3969/j.issn.1000-2561.2022.05.002
荔枝（Litchi chinensis Sonn.）是热带、亚热带发展中国家的最重要的作物之一。

目前荔枝产业已成为中国华南热作产业的支柱产业之一，对中国南方农村经济的发展和热区农民收入水平的提高发挥了举足轻重的作用。

2018年是我国有史以来最丰产的一年，据不完全统计，全国荔枝种植面积约为54.2万hm2，产量约为301万t，种植面积排在苹果、柑橘、梨、葡萄和桃之后，属于大宗果树，但从总产量来看却还算“小水果”。

2018年荔枝结果面积按照总种植面积的80%计算，平均单产不足7.0 t/hm2[1]。

‘妃子笑’具有花穗长、花量大、坐果率低的特点，坐果率低是其生产中亟待解决的一个“瓶颈”问题，调控坐果的分子机理是当今果实发育领域研究的一个热点问题。

泛素-蛋白酶体途径（ubiquitin-proteasome pathway，UPP）广泛存在于真核生物中，并且在维持细胞功能、细胞周期运转、抵御环境胁迫、胚胎发育、激素响应和衰老等方面发挥着重要的作用[2-5]。

在泛素化过程中，关键酶主要包括泛素活化酶（E1）、泛素结合酶（E2）和泛素-蛋白连接酶（E3）。

编码E2及E2-like的基因已报道37个，根据结构和序列的相似性可分为12个亚家族或14个组，E2-like基因8个[6]。

E2蛋白的典型特征是包含一个长约140~200 aa的保守催化结构域——UBC结构域，内含有1个高度保守的半胱氨酸活性位点，另外还存在一类UEV （ubiqutin E2 variants）蛋白，该蛋白家族在序列和结构上与泛素结合酶相似，但是缺少半胱氨酸催化位点，其功能也与典型的泛素结合酶蛋白有所不同。

E2多以基因家族的形式存在，参与许多植物生理活动[7]。

泛素结合酶在泛素化级联反应中居于中心位置，是连接泛素活化酶和泛素连接酶的桥梁，在泛素化系统中发挥关键的作用。

目前，对E2基因的研究多集中在对DNA的修复进程。

如在核酸修复途径中，UBC35/UBC36和 Mms2p/Uev1a蛋白形成复合物，在赖氨酸位点催化形成泛素链，来引发核酸修复的信号转导，完成对受损核酸的修复[8]。

UBC13协同E3连接酶RNF 8相互作用形成复合物，DNA损伤信号通过依赖泛素化的信号途径来传递至染色质，引发染色质蛋白的磷酸化，促使RNF 8-UBC 13复合物定位到受损伤的核酸部位，产生泛素链后触发DNA修复信号通路，使RAP 80蛋白和整个Brca1A复合物定位到受损伤的DNA部位，对损伤的DNA进行修复[9]。

另外还有研究表明RNF 8-UBC 13复合物可独立的行使底物蛋白的泛素化，并促进DNA损伤的调节子53BP1蛋白的积累[10]。

泛素结合酶在植物衰老和果实发育过程中起重要作用。

PICTON等[11]首次从番茄中克隆到泛素结合酶基因，研究发现泛素结合酶基因在番茄叶片衰老和果实成熟过程中表达。

香蕉采后早期抑制差减文库中筛选得到香蕉中的泛素结合酶基因MaUCEI，该基因介导细胞内小的调节蛋白的降解，与酵母中的UBC 4/UBC 5结合，推测该基
All Rights Reserved.
884 热带作物学报第43卷
因可能通过影响果实成熟过程中淀粉酶的活性或含量而参与对淀粉降解的调控，在果实成熟衰老的分解代谢中发挥作用[2]。

2013年日本学者研究表明，大豆近等基因系成熟基因座等位基因E2/E3对基因型的修饰可以提供新品种，与对照相比，新品种花期较晚且产量高，可以适应不同的纬度环境，并可保留原来的种子质量[12]。

WANG等[13]研究发现，番茄中2个E2基因（SlUBC 32和SlUBC 41）参与调控番茄果实成熟，但具体分子调控机制尚不清楚。

DONG等[14]研究香蕉花后果实不同发育阶段E2基因家族的表达情况，结果发现E2基因家族不同成员在果实发育不同阶段行使功能，推测这些基因在果实发育过程中发挥关键作用。

陈曙等[15]研究低磷处理下玉米幼苗不同组织中E2基因的表达情况，结果发现ZmUBC 17基因具有组织表达特异性，在叶片中大量表达，推测ZmUBC 17基因是通过调控玉米对磷元素吸收和转运间接影响光合作用效率。

通过对拟南芥基因组数据库搜索显示，E2基因以多拷贝形式存在，造成E2蛋白功能冗余，这可能是使E2基因的研究进展不大的主要原因[7]。

近年来，随着基因组测序技术的不断提高，越来越多的植物完成全基因组的测序工作，越来越多的植物物种中泛素结合酶基因家族在基因组学上得到鉴定，如玉米[16]、可可[17]、水稻[18]、葡萄[19]、普通油茶[20]和香蕉[14]等。

目前，关于荔枝基因家族的研究有XTH[21]、ARF[22-23]基因家族，关于荔枝泛素结合酶基因家族仅仅研究了LcUBC12基因[24]。

本研究基于荔枝转录组数据库，利用生物信息学方法对转录组数据库UBC基因家族进行鉴定，并通过荧光定量PCR（qPCR）技术对LcUBC 基因家族在‘妃子笑’荔枝不同组织中进行时空表达分析，并对花穗对照和烯效唑处理花穗发育的不同时期进行表达分析，挖掘关键的LcUBC基因家族成员，为今后LcUBC基因功能的研究奠定基础。

本研究结果将为从蛋白质降解方面深化荔枝坐果机理的研究和生产上调控荔枝坐果技术措施提供理论依据。

1 材料与方法
1.1 材料
供试品种为‘妃子笑’荔枝，取自广东省湛江市麻章区湖光镇湖秀路一号南亚热带作物研究所荔枝栽培示范园，荔枝UBC蛋白序列来源于本课题组前期构建的‘妃子笑’花穗发育转录组数据库（GenBank accession SRP092890）。

1.2 方法
1.2.1 LcUBC基因家族的鉴定与分析利用模式植物拟南芥和水稻UBC序列作为探针序列，运用本地Blast软件对荔枝测序转录组数据库中的unigene进行搜索；同时，利用关键词“Ubiquitin- conjugating enzymes”和“UBC”直接检索，得到所有可能的UBC序列后，对检索结果进行整合分析，去除重复序列后，利用Pfam在线软件对序列保守结构域进行鉴定，去除不含LcUBC基因家族保守结构域的序列。

1.2.2 LcUBC蛋白的基本理化性质、序列比对和进化分析运用ProtParam在线软件（http://web. /protparam/）分析预测LcUBC蛋白序列的分子量、等电点、不稳定系数、脂肪指数和疏水性等基本的理化性质；采用Plant-mPLoc Server 在线软件（/bioinf/ plant-multi/#）分析LcUBC蛋白的亚细胞定位；运用MEME 4.12.0在线软件（/ meme/）分析LcUBC家族蛋白的保守基序，参数设置：基序的最大数目设置为10, 基序长度设为6~150个氨基酸，其他参数为默认值。

通过MEGA 6.0的MUSCLE程序比对氨基酸多重序列，参数为程序默认参数。

利用多重序列比对产生的LcUBC保守结构域145个氨基酸残基后续进行NJ（Neighbor-Joining）分析。

采用MEGA 6.0软件构建进化树，Bootstrap值设为1000次重复，以获得更为可靠的分支聚类，其他参数为默认参数。

1.2.3 LcUBC基因家族在花穗调控下的表达特征分析选择树龄相同、长势相近的‘妃子笑’荔枝树，分别选取荔枝树的不同组织如嫩叶、根、茎、老叶、雌花、雄花、果肉、果皮、种子，液氮速冻后存放于–80℃冰箱中，用于RNA的提取。

花穗处理：选择树龄相同、长势相近的‘妃子笑’荔枝树，当花穗长度抽生至约18 cm时，进行如下操作：（1）对照组（CK），不做任何处理；（2）烯效唑（Un）处理，用50 mg/L烯效唑（品牌：四川国光，有效成分含量5%）喷施，喷至叶面滴水；分别取对照及烯效唑处理后0、7、14、21、28、35、42 d的花穗，液氮速冻后存放于–80℃冰箱中，用于RNA的提取。

选用前期课题组筛选的LcActin为内参，根据转录组中注释的LcUBC的核苷酸序列，利用在线
All Rights Reserved.
第5期董晨等：基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析 885
软件Primer 3 plus设计荧光定量引物（表1），通过Blast分析引物的特异性，引物序列委托广州艾基生物技术有限公司合成。

将处理好的材料按照RNA提取标准进行磨样，采用华越洋生物科技有限公司RAN提取试剂盒QK型，分别对根、茎、嫩叶、老叶、雌花、雄花、果肉、果皮、种子和花穗的总RNA进行提取，提取1 μL检测RNA的质量和浓度，样品符合要求后利用M-MLV逆转录酶（TaKaRa）反转录合成cDNA。

qPCR分析采用SYBR Premix Ex Taq™（TaKaRa）试剂盒，Roche480Ⅱ定量PCR系统（瑞士），采用384孔板，qPCR反应体系为10 μL，其中cDNA 1 μL（相当于25 ng总RNA），10 µmol/L基因特异性引物1 µL，2×SYBR Premix ExTaq 5 μL，最后用ddH2O 补足10 μL。

qPCR反应条件：50℃预变性2 min，95℃预变性2 min；95℃变性15 s，56℃退火15 s，72℃延伸30 s，共55个循环。

每个样品3个重复。

采用2–∆∆C T法计算基因相对表达量。

表1 LcUBC基因家族引物序列
Tab. 1 Primer sequences of LcUBC gene family
基因Gene ID 正向引物（5'–3'）Forward primer (5'–3') 反向引物（5'–3'）Reverse primer (5'–3')
LcUBC1 Unigene0001844 TCGTACAGTTGGAGGCTTGT TGCAAGGAAGAGGCTGATGA LcUBC2 Unigene0001927 TTCGAAATGGTCAGTGCTGG TTGGCAAGCGACGATAATGG LcUBC3 Unigene0003053 AGCTAGCTGTGTGACCATCA ACCTCAGCCCTGGAAAGAAA LcUBC4 Unigene0003214 GAAAGGTGATGCTGGTGGTG GCCCTCAGAAGTCCAATCCT LcUBC5 Unigene0010286 TCTCAAGCAAGGGAATCGGT ATGAGGTCCCATGTGATGCA LcUBC6 Unigene0014473 CAGGAGGTCAAGGAGATGCA ACGGGTACTCAGCAGGTAAC LcUBC7 Unigene0014483 CATTTGCCTCCTGACCCATG GCTACTGCCTCATCACTGGA LcUBC8 Unigene0014555 GGGAGGGCTGTTGGGATAAT ATTCATGGCAGCGTCTCAAG LcUBC9 Unigene0015820 GGTTGTTTCGTCCAGCACAT TGGGTTTGGATCTGTGAGCA LcUBC10 Unigene0016542 ATTGACCGAGGGCTCTGAAA ACTCTTCCCATTTGCCCAGA LcUBC11 Unigene0029596 CACCTGCTTCCATCCCAATG TGTGCTGCTTGGCTGTTTAG LcUBC12 Unigene0032422 ATCCTCCTACCTCTTGCAGC CCTTTGTCCTGAATGCGACC LcUBC13 Unigene0032423 TAGTCCTTACGCTGGTGGTG ACCAATGGATCATCGGGGTT LcUBC14 Unigene0033691 GTGATCGAAGTCAACGGAGC CTAACAGTCATGGCAGGGGA LcUBC15 Unigene0034885 TGCTTCGGTTTCTCAGAGGT TGGATGGAAGCACATGGTCT LcUBC16 Unigene0035356 ACGAAGGAGGTGTCTGGAAA TACCAAACATCGGGCTCCAT LcUBC17 Unigene0036703 CCTCTCTTTTCCAATGCCCG ACATGTCATCAGCCACAGGA LcUBC18 Unigene0038283 AGCGGACTTCTTCTCCCAAA TGAGCCTTTTGAACACGTGG LcUBC19 Unigene0039504 CGAAGCAACAAACTGGACCA GTGAAACCATCCCATCTGCC LcUBC20 Unigene0039859 AGGAACTCACAGATCTCGGC GGCTTGAAGGGGTAGTCTGT LcUBC21 Unigene0039860 CTCAAGCGCATCAACAAGGA AGTCTGTGGGGAAGTGGATG LcUBC22 Unigene0042419 AAATTCTGCGATCCTGTGCC AGATTCCCTCCTCCGCTTTC LcUBC23 Unigene0042906 TGTCTGCACCTGAACCTGAT CACCACCGACAGCTTCTAGA LcUBC24 Unigene0043251 CTCCGAGAAATGGCCAACAG AATCGAACCTTTGGAGCAGC LcUBC25 Unigene0044469 GGAGCTAGTGGTGATGGTGA CATGCCAGGATTTACACCGG LcUBC26 Unigene0050724 GGAATTGCAGAAATGGCCGA GCAACAATAAGCGTGCAGTC LcUBC27 Unigene0051029 AGAGTTCTGTGACCAAGCGA CAGAAGAAGTTGTGGCCGAA LcUBC28 Unigene0052502 TCTGTGTCATGAACCGTCCA ATGTTGTTGCACGTCCTTCC
2 结果与分析
2.1 LcUBC基因家族的鉴定
首先对转录组数据库的本地Blast进行比对检索和关键词搜素，再通过SMART对保守结构域进行筛选，最终得到28条LcUBC蛋白序列。

为了便于研究，按照转录组中的基因ID大小前后顺序统一为LcUBC蛋白编号为LcUBC 1~LcUBC 28（表2），对应所编码的氨基酸数分布在85~
All Rights Reserved.
886 热带作物学报第43卷
1149 aa之间；预测的理论分子量大小范围在9.668~102.307 kDa之间。

等电点大小范围在4.03~9.61之间，其中 18个LcUBC蛋白的等电点小于7，显酸性；10个LcUBC蛋白等电点大于7，显碱性；LcUBC蛋白平均等电点小于7，表明UBC蛋白为弱酸性，在酸性的亚细胞环境中发挥作用。

不稳定指数分析结果发现，5个LcUBC蛋白（LcUBC 4、LcUBC 9、LcUBC 18、LcUBC 23和LcUBC 25）的不稳定指数小于40，为稳定蛋白，其余均为不稳定蛋白。

平均亲水系数结果表明，除了LcUBC 4蛋白疏水，其余均为亲水蛋白。

亚细胞定位分析结果表明，LcUBC 7和LcUBC 16蛋白定位于细胞质，LcUBC 1和LcUBC 24蛋白定位于细胞质和细胞核，LcUBC 6蛋白定位于内质网，其余蛋白均定位于细胞核。

综合上述分析结果，无论从氨基酸序列的长度还是蛋白的特性变化分析，LcUBC蛋白都有很大差异，这表明LcUBC基因家族蛋白具有不同特性。

表2 LcUBC基因家族信息Tab. 2 LcUBC gene family information
基因Gene
转录组ID
Transcriptome ID
UBC结构域
UBC
domain
蛋白长度
Protein
length/aa
分子量
Molecular
weight/Da
等电点
pI
不稳定指数
Instability
index
脂肪系数
Aliphatic
indexr
平均亲水系数
Average hydro-
philic coefficient
亚细胞定位
Subcellular
localization
LcUBC1 Unigene0001844 11-154 156 17
792.11 5.8369.97 70.06 –0.568 Cytoplasm,
Nucleus LcUBC2 Unigene0001927 4-147 148 16
492.067.7245.84 77.16 –0.272 Nucleus LcUBC3 Unigene0003053 1-68 104 11
765.79 4.0342.97 69.42 –0.830 Nucleus LcUBC4 Unigene0003214 5-170 260 29
379.268.2933.19 84.69 0.002 Nucleus LcUBC5 Unigene0010286 1-153 282 31
162.01 5.6448.65 76.77 –0.651 Nucleus LcUBC6 Unigene0014473 14-174 276 31
084.15 5.2557.42 72.83 –0.629 Endoplasmic
reticulum LcUBC7 Unigene0014483 4-145 145 16
772.25 4.8647.28 69.79 –0.401 Cytoplasm LcUBC8 Unigene0014555 6-163 165 18
454.92 5.2047.60 77.94 –0.353 Nucleus LcUBC9 Unigene0015820 3-84 85 9
668.187.9937.45 91.76 –0.165 Nucleus LcUBC10 Unigene0016542 676-833 916 102
307.90 4.7940.20 77.14 –0.336 Nucleus LcUBC11 Unigene0029596 36-178 178 19
701.27 5.0262.66 75.62 –0.363 Nucleus LcUBC12 Unigene0032422 4-147 148 16
548.137.7249.03 75.81 –0.312 Nucleus LcUBC13 Unigene0032423 4-147 148 16
575.157.7248.49 75.81 –0.330 Nucleus LcUBC14 Unigene0033691 18-161 161 18
341.888.3568.35 76.89 –0.455 Nucleus LcUBC15 Unigene0034885 38-181 185 20
677.37 5.5856.17 66.92 –0.470 Nucleus LcUBC16 Unigene0035356 3-96 97 11
042.66 5.0249.56 76.19 –0.327 Cytoplasm LcUBC17 Unigene0036703 4-147 148 16
498.077.6844.66 80.41 –0.267 Nucleus LcUBC18 Unigene0038283 61-218 351 39
250.809.4937.79 80.83 –0.326 Nucleus LcUBC19 Unigene0039504 895-1055 1146 126
438.45 4.5740.24 75.75 –0.570 Nucleus LcUBC20 Unigene0039859 4-147 147 16
495.81 6.4041.47 79.66 –0.324 Nucleus LcUBC21 Unigene0039860 4-115 115 12
662.53 6.2344.45 83.13 –0.131 Nucleus LcUBC22 Unigene0042419 299-484 581 64
792.80 5.3046.48 77.30 –0.350 Nucleus LcUBC23 Unigene0042906 5-150 193 21
076.06 4.7237.71 89.90 –0.129 Nucleus LcUBC24 Unigene0043251 8-151 153 17
205.85 6.7449.03 92.55 –0.319 Cytoplasm,
Nucleus LcUBC25 Unigene0044469 2-81 195 20
831.029.6138.30 89.69 –0.268 Nucleus LcUBC26 Unigene0050724 18-130 130 14
519.60 6.1658.43 89.23 –0.135 Nucleus LcUBC27 Unigene0051029 30-178 209 24
455.55 4.4450.54 81.58 –0.330 Nucleus LcUBC28 Unigene0052502 8-162 238 27
166.319.0759.69 66.30 –0.511 Nucleus All Rights Reserved.
第5期董晨等：基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析 887
2.2 LcUBC基因家族蛋白保守基序分析
通过多重序列比对LcUBC基因家族的氨基酸序列全长，比较不同LcUBC家族成员间UBC 保守结构域的特点（图1），多重序列比对中，28个LcUBC基因均含有进化上高度保守的氨基酸残基半胱氨酸的活性位点，位于半胱氨酸活性位点后的C-端色氨酸（W）在28个LcUBC蛋白中也高度保守，“HPNI”保守基序位于半胱氨酸活性位点之前的7~10位氨基酸。

图1 LcUBC基因家族成员蛋白保守区的
多重序列比对分析
Fig. 1 Multiple sequence alignment analysis of LcUBC gene family members conserved regions
通过在线软件MEME分析LcUBC蛋白中的保守基序，结果见图2。

预测LcUBC基因家族中含有10个保守基序，LcUBC基因家族成员之间所包含的保守基序数目及种类存在一定的差异。

在LcUBC基因家族中有25个成员含Motif 1基序，有15个成员含Motif 2基序。

通过SMART 分析Motif 1、Motif 2为LcUBC蛋白结构域中非常重要的保守基序。

为进一步了解LcUBC基因家族成员之间的进化关系，利用模式物种拟南芥中的37个AtUBC 和酵母中的13个ScUBC同荔枝中的28个LcUBC 构建系统进化树（图3），结果表明，荔枝中的28个LcUBC分为10个亚家族，分别为UBC1、UBC2/11、UBC3/7、UBC4/5、UBC6、UBC8、UBC10/12、UBC13、UBC14、UBC15，其中UBC14和UBC15在酵母基因组中缺失，推测UBC14和UBC15是植物中特有的或是酵母进化过程中丢失。

2.3 LcUBC基因家族时空表达分析
UBC参与DNA的损失修复、生长、发育及胁迫响应等方面，为了研究LcUBC基因家族在不同发育过程中的功能，分析9个不同组织（如嫩叶、嫩茎、果皮、种子、根、雄花、老叶、雌花、果肉）中LcUBC基因的表达情况（图4）。

在植物正常的发育过程中，不是所有的UBC基因在不同组织部位都表达，28个LcUBC基因家族成员中至少在某一个组织部位表达（图5），推测这些基因有可能是在生物胁迫或非生物胁迫条件下表达。

其中，LcUBC4、LcUBC5、LcUBC18在种子中的表达显著高于其他组织，LcUBC7在雄花中的表达上显著较高于其他组织，LcUBC14、LcUBC15在雌花中的表达量显著高于其他组织，LcUBC10在果皮和老叶中高表达。

2.4 LcUBC基因家族在花穗调控表达模式
利用qPCR分析LcUBC基因家族在荔枝花穗调控的表达量（图5），结果表明，与对照相比，LcUBC基因家族在烯效唑处理后表达量存在显著变化，经烯效唑处理后，LcUBC27在7、14、21、28、35 d表达量上调，LcUBC2、LcUBC7、LcUBC12、LcUBC13、LcUBC15、LcUBC16、LcUBC17、LcUBC28在21 d表达量最高，LcUBC 20、LcUBC21在35 d表达量最高，LcUBC16表达呈现先上升后下降的趋势，在21 d表达量最高。

LcUBC12、LcUBC 23、LcUBC24、LcUBC25与对照相比表达有变化，但差异不显著。

3 讨论
本研究参考拟南芥和水稻UBC家族基因的信息，基于烯效唑处理荔枝花穗转录组数据，利用生物信息学方法全面分析鉴定LcUBC基因家族。

结果共鉴定出28个LcUBC基因家族成员，并对其基本理化性质、亚细胞定位进行了分析。

通过对LcUBC基因家族蛋白保守结构域分析可看出，LcUBC蛋白的结构域具有高度的保守性，包含典型的保守结构域UBC domain，且含有半胱氨酸活性位点，这与香蕉中UBC蛋白的分析结果一致[14]。

All Rights Reserved.
888 热带作物学报第43卷All Rights Reserved.
A：LcUBC保守基序位置；B：保守基序序列。

A: Conserved motif location of LcUBC; B: Conserved motif consensus.
图2 LcUBC基因家族成员保守基序分析
Fig. 2 Distribution of conserved motifs in LcUBC gene family members
第5期董晨等：基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析 889
All Rights Reserved.
图3荔枝、拟南芥和酵母中UBC的邻接法系统进化树
Fig. 3 Neighbor-joining phylogenetic tree of UBC in Litchi, Arabidopsis and Yeast
图4 LcUBC基因家族时空表达分析
Fig. 4 Spatiotemporal expression analysis of LcUBC gene family
890 热带作物学报第43卷All Rights Reserved.
*表示差异显著（P<0.05）。

* indicates significant difference (P<0.05).
图5烯效唑处理花穗发育不同时期LcUBC基因家族表达分析
Fig. 5 Analysis of LcUBC gene family expression at different stages of flower panicle development in uniconazole treatment
第5期董晨等：基于转录组的荔枝泛素结合酶基因家族的鉴定及表达分析 891
通过转录组及qPCR的结果分析，在荔枝的时空表达和烯效唑处理后的花穗调控中，大多数LcUBC基因家族的表达水平发生了变化，推测这些基因很可能通过参与激素信号转导途径或调控果实相关转录因子的水平，进而影响荔枝果实成熟过程中的各种生理变化。

LcUBC 4、LcUBC 5、LcUBC8、LcUBC18、LcUBC19在种子中表达量较高，推测这5个LcUBC家族基因在果实种子发育过程中发挥着重要的作用；LcUBC 6、LcUBC14、LcUBC15、LcUBC23、LcUBC24在雌花中表达量较高，LcUBC 7、LcUBC 26在雄花中表达量较高，推测这7个LcUBC家族基因在果实花发育过程中发挥着重要的作用；LcUBC 10、LcUBC12、LcUBC18、LcUBC27在老叶中表达量较高，而LcUBC 13、LcUBC17、LcUBC20、LcUBC21在果皮中表达量较高。

多个LcUBC家族基因的表达水平在烯效唑处理后第21天的表达量较高，这表明LcUBC基因家族在该时期调控荔枝花穗发育过程中发挥着关键作用。

该研究结果为进一步研究LcUBC基因家族的生物学功能奠定了坚实的基础。

参考文献
[1] 李建国, 王惠聪, 周碧燕, 赵明磊, 李彩琴, 夏瑞, 黄旭
明. 荔枝花果发育生理和分子生物学研究进展[J]. 华南农
业大学学报, 2019, 40(5): 119-127.
LI J G, WANG H C, ZHOU B Y, ZHAO M L, LI C Q, XIA
R, HUANG X M. Research advances in physiology and mo-
lecular biology of flower and fruit development in litchi[J].
Journal of South China Agricultural University, 2019, 40(5): 119-127. (in Chinese)
[2] 董发才, 宋纯鹏. 植物细胞中的泛素及其生理功能[J]. 植
物生理学报, 1999, 35(1): 54-59.
DONG F C, SONG C P. The ubiquitin and its physilolgical functions in plants[J]. Plant Physiology Journal, 1999, 35(1):
54-59. (in Chinese)
[3] 覃碧. 泛素蛋白酶体途径及其对植物激素信号转导的
调控[J]. 热带农业科学, 2013, 33(4): 39-45.
QI B. Ubiquitin-proteasome pathway and its regulation of plant hormone signaling[J]. Chinese Journal of Tropical Ag-
riculture, 2013, 33(4): 39-45. (in Chinese)
[4] 王园, 王甲水, 谢学立, 雷晓明, 金志强. 香蕉泛素结
合酶基因与果实成熟关系的研究[J]. 园艺学报, 2010, 37(5): 705-712.
WANG Y, WANG J S, XIE X L, LEI X M, JIN Z Q, Studies
of the relationship between MaUCE1 and banana fruit rip-
ening[J]. Acta Horticulturae Sinica, 2010, 37(5): 705-712.
(in Chinese)
[5] 盛仙永, 胡正海, 林金星. 泛素/蛋白酶体途径及其在高等
植物有性生殖中的作用[J]. 西北植物学报, 2004, 24(8):
1527-1536.
SHENG X Y, HU Z H．LIN J X. Ubiquitin/proteasome
pathway and it’s role in high plant sexual reproduction[J].
Acta Botanica Boreali-Occidentalia Sinica, 2004, 24(8):
1527-1536. (in Chinese)
[6] BACHMAIR A, NOVATCHKOVA M, POTUSCHAK T,
EISENHABER F. Ubiquitylation in plants: a post-genomic
look at a post-translational modification[J]. Trends in Plant
Science, 2001, 6(10): 463-470.
[7] 王金利, 史胜青, 贾利强, 江泽平. 植物泛素结合酶E2功
能研究进展[J]. 生物技术通报, 2010(4): 7-10.
WANG J L, SHI S Q, JIA L Q, JIANG Z P. Progress on
functions of ubiquitin-conjugating enzyme (E2) in plants[J].
Biotechnology Bulletin, 2010(4): 7-10. (in Chinese)
[8] HOEGE C, PFANDER B, MOLDOVAN G,
PYROWOLAKIS G, JENTSCH S. RAD6-dependent DNA
repair is linked to modification of PCNA by ubiquitin and
SUMO[J]. Nature, 2002, 419(6903): 135-141.
[9] BIN W, STEPHEN J E. Ubc13/Rnf8 ubiquitin ligases con-
trol foci formation of the Rap80/Abraxas/Brca1/Brcc36
complex in response to DNA damage[J]. Proceedings of the
National Academy of Sciences of the United States of
America, 2007, 104(52): 20759-20763.
[10] HUEN M, HUANG J, YUAN J S, YAMAMOTO M,
AKIRA S Z, ASHLEY C, XIAO W, CHEN J J. Noncanoni-
cal E2 variant-independent function of UBC13 in promoting
checkpoint protein assembly[J]. Molecular & Cellular Biol-
ogy, 2008, 28(19): 6104-6112.
[11] PICTON S, GRAY J E, LOWE A, BARTON S L,
GRIERSON D. Sequence of a cloned tomato ubiquitin con-
jugating enzyme[J].
Plant Physiology,
1993, 103(4): 1471-1472.
[12] YAMADA T, HAJIKA M, YAMADA N, HIRATA K,
OKABE A, OKI N, TAKAHASHI K, SEKI K, OKANO K,
FUJITA Y, KAGA A, SHIMIZU T, SAYAMA T,
ISHIMOTO M. Effects on flowering and seed yield of dom-
inant alleles at maturity loci E2 and E3 in a Japanese cultivar,
Enrei[J]. Breeding Science, 2012, 61(5): 653-660.
[13] WANG Y Y, WANG W H, CAI J H, ZHANG Y R, QIN G
Z, TIAN S P. Tomato nuclear proteome reveals the involve-
ment of specific E2 ubiquitin-conjugating enzymes in fruit
ripening[J]. Genome Biology, 2014, 15(12): 548.
[14] DONG C, HU H G, JUE D W, ZHAO Q F, CHEN H L, XIE
All Rights Reserved.
892 热带作物学报第43卷
J H, JIA L Q. The banana E2 gene family: Genomic identi-
fication, characterization, expression profiling analysis[J].
Plant Science, 2016, 245: 11-24.
[15] 陈曙, 赵秋芳, 陈宏良, 贾利强. 玉米泛素结合酶基因
家族的生物信息学及表达分析[J]. 南方农业学报, 2018, 49(8): 1476-1483.
CHEN S, ZHAO Q F, CHEN H L, JIA L Q. Bioinformatics
and expression analysis of ubiquitin conjugating enzyme gene family in maize[J]. Journal of Southern Agriculture, 2018, 49(8): 1476-1483. (in Chinese)
[16] JUE D W, SANG X L, LU S Q, DONG C, ZHAO Q F,
CHEN H L, JIA L Q. Genome-wide identification, phyloge-
netic and expression analyses of the ubiquitin-conjugating enzyme gene family in maize[J]. PLoS One, 2015, 10(11): e143488.
[17] 赵秋芳, 董晨, 决登伟, 陈宏良, 贾利强. 可可泛素结
合酶基因家族的生物信息学初步分析[J]. 热带作物学报,
2016, 37(9): 1732-1740.
ZHAO Q F, DONG C, JUE D W, CHEN H L, JIA L Q. Bi-
oinformatics preliminary analysis of ubiquitin conjugating enzyme gene family in cocoa[J]. Chineese Journal of Tropi-
cal Crops, 2016, 37(9): 1732-1740. (in Chinese)
[18] 刘鑫, 张恒, 阚虎飞, 周立帅, 黄昊, 宋林林, 翟焕
趁, 张君, 鲁国东. 水稻泛素结合酶基因家族的生物信
息学与表达分析[J]. 中国水稻科学, 2016, 30(3): 223-231.
LIU X, ZHANG H, KAN H F, ZHOU L S, HUANG H, SONG L L, ZHAI H C, ZHANG J, LU G D. Bioinformatic and expression analysis of rice ubiquitin-conjugating enzyme
gene family[J]. Chinese Journal of Rice Science, 2016, 30(3): 223-231. (in Chinese)
[19] 高营营. 葡萄泛素结合酶E2基因家族全基因组鉴定及功
能结构分析[D]. 武汉: 中国科学院大学(中国科学院武汉
植物园), 2018.
GAO Y Y. Genome-wide identification and func-
tion-structure analysis of ubiquitin-conjugating E2 genes family in grape[D]. Wuhan: University of Chinese Academy
of Sciences (Wuhan Botanical Garden, Chinese Academy of
Sciences), 2018. (in Chinese) [20] 林萍, 姚小华, 曹永庆, 龙伟, 王开良, 滕建华. 普通
油茶泛素结合酶UBE2-J2的cDNA序列及蛋白质结构分
析[J]. 林业科学研究, 2013, 26(6): 744-751.
LIN P, YAO X H, CAO Y Q, LONG W, WANG K L, TENG J H. Characterization of a novel ubiquitin-conjugating
enzyme from Camellia oleifera[J]. Forest Research, 2013, 26(6): 744-751. (in Chinese)
[21] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才. 基于转录组
的荔枝XTH家族基因的鉴定及分析[J]. 分子植物育种,
2019, 17(12): 3865-3873.
DONG C, WEI Y Z, WANG Y, ZHENG X W, LI W C.
Transcriptome-wide identification and analysis of xyloglu-
can endotransglucosylase/hydrolase (XTH) family gene in Litchi chinensis Sonn.[J]. Molecular Plant Breeding, 2019, 17(12): 3865-3873. (in Chinese)
[22] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才. 基于转录组
的荔枝ARF基因家族的鉴定及表达分析[J]. 热带作物学
报, 2017, 38(8): 1485-1491.
DONG C, WEI Y Z, WANG Y, ZHENG X W, LI W C.
Transcriptome-wide identification and expression profiling of the auxin response factor (ARF) gene family in Litchi
chinensis Sonn.[J]. Chinese Journal of Tropical Crops, 2017,
38(8): 1485-1491. (in Chinese)
[23] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才. 不同处理方
式的妃子笑荔枝花穗ARF基因家族的表达分析[J]. 热带
作物学报, 2018, 39(6): 1122-1127.
DONG C, WEI Y Z, WANG Y, ZHENG X W, LI W C. Ex-
pression analysis of auxin response factor (ARF) gene family
in Litchi chinesis Sonn.[J]. Chinese Journal of Tropical Crops, 2018, 39(6): 1122-1127. (in Chinese)
[24] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才, 石胜友. 荔枝
泛素结合酶基因(LcUBC12)的生物信息学及表达特性分
析[J]. 南方农业学报, 2020, 51(5): 1091-1097.
DONG C, WEI Y Z, WANG Y, ZHENG X W, LI W C, SHI
S Y. Bioinformatics and expression characteristics analysis of Litchi chinesis Sonn. ubiquitin-conjugating enzyme gene (LcUBC12)[J]. Journal of Southern Agriculture, 2020, 51(5):
1091-1097. (in Chinese)
All Rights Reserved.。