畜牧与兽医2020年第52卷第2期 于农淇,鄢胜飞,覃广胜,等.水牛CART基因的克隆测序及蛋白质结构与功能预测[J].畜牧与兽医,2020,52(2):1-6 YU N Q, YAN S F, QIN G S, et al. Cloning, sequencing and bioinformatics analysis of the CART gene in buffalo [J. Animal Husbandry Veterinary Medi- cine,2020,52(2):1-6. 水牛CART基因的克隆测序及蛋白质结构与功能预测 于农淇,鄢胜飞,覃广胜,尚江华,郑海英,杨春艳 (中国农业科学院广西水牛研究所/广西水牛遗传繁育重点实验室,广西南宁530001) 摘要:为探讨水牛可卡因-苯丙胺调节转录肽( cocaine- and amphetamine-regulated transcript,CART)基因功能,以水牛卵巢组织RNA转录的 cDNA为模板,参考 GenBank预测的水牛CART基因组序列设计引物,采用PCR和在线生信分析软件对CART基因进行克隆和蛋白质结构预测分 析。结果显示:水牛CART基因CDs区全序列长为351bp,编码116个氨基酸;对所得核酸序列和蛋白质结构分析预测显示,水牛CART基因与 美洲野牛、野牦牛、普通牛、羊、双峰驼、单峰驼、野猪、马、非洲野驴、人类及猕猴相似程度分别为100%、%8%、100%、98%、91%、9%」 78%、85%、85%、81%、90%;Me分析显示在多胃反刍动物间有较高同源性;蛋白质结构预测显示,CART蛋白含有N端信号肽,是分泌蛋 白,功能与神经、辅酶和离子结合相关。 关键词:CART;水牛;蛋白质功能预测 中图分类号:S823.83 文献标志码:A文章编号:0529-5130(2020)02-0001-06 Cloning, sequencing and bioinformatics analysis of the cart gene in buffalo YU Nongqi, YAN Shengfei, QIN Guangsheng, SHANG Jianghua, ZHENG Haiying, YANG Chunyan Guangxi Key Laboratory of Buffalo Genetics, Breeding and Reproduction/Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China) Abstract: In this study, the cocaine-and amphetamine-regulated transcript( CART) gene of buffalos was cloned and its protein structure was predicted. The cDNA transcribed by RNA was extracted from the ovaries of buffalos as templates, primers were designed according to the genomic sequence of CARt predicted by GenBank so as to amplify and sequence the CDS region of the Cart gene of theanimals. The results showed that the total length of the CDS region of the buffalo CART gene was 351 bp, encoding 116 amino acids. The analysis of nucleic acid sequence and protein structure indicated that the similarities between the buffalo Cart gene and that of American bison, wild yak, common cattle, sheep, bactrian camel and single hump camel, wild boar, horse and African wild donkey, human and macaque were 100%, 98% 100%6,98%,91%,91%, 78%,85%,85%,81% and 90%, respectively. The prediction of protein structure indicated that the CArT pro- tein contained N-terminal signal peptides confirming it to be a secretory protein. The present study laid a foundation for further research on the function of the Cart gene in buffalos Keywords: CART; buffalo; protein function prediction 可卡因-苯丙胺调节转录肽( cocaine- and an关2-4。水牛CART基因( Gene id:102413682,NW phetamine- regulated transcript,CART)广泛分布于中005785813)全长为2059bp,有3个CDS区,mRNA 枢神经系统和外周神经系统,被认为是参与食物摄入全长775b,CDS全长312bp。前人对水牛CART及 调控的主要无食欲肽。CART是一种关键的神经递其基因的研究不多,尤其缺少奶水牛上该基因结构、 质和激素,参与调节多种生物过程,包括食物摄入、表达蛋白分析的相关研究。本研究通过克隆奶水牛 维持体重、奖励和成瘾、应激反应、心理刺激效应和CART基因CDS区序列,分析研究其编码蛋白质的氨 内分泌功能,并被报道与家畜繁殖力、生产性能有基酸结构及蛋白定位、互作,丰富了CART基因数据 信息,为今后CART对水牛内分泌和生长繁殖作用的 收稿日期:2019-03-21;修回日期:2019-12-16 研究提供数据支撑。 基因项目:广西科技计划项目(桂科AB16380040);广西水牛研 所基本科研业务费项目(水牛基180503) 1材料与方法 第一作者:于农淇(1991-):女,壮族,硕士 令·通信作者:杨春艳,副研究员,研究方向为动物遗传育种与繁1.1材料 E-mail:86748227(@qg.como 试验所用水牛卵巢组织样本来自广西南宁屠宰
于农淇, 鄢胜飞, 覃广胜, 等. 水牛 CART 基因的克隆测序及蛋白质结构与功能预测 [J]. 畜牧与兽医, 2020, 52 (2): 1-6. YU N Q, YAN S F, QIN G S, et al. Cloning, sequencing and bioinformatics analysis of the CART gene in buffalo [J]. Animal Husbandry & Veterinary Medi⁃ cine, 2020, 52 (2): 1-6. 水牛 CART 基因的克隆测序及蛋白质结构与功能预测 于农淇, 鄢胜飞, 覃广胜, 尚江华, 郑海英, 杨春艳∗ (中国农业科学院广西水牛研究所/ 广西水牛遗传繁育重点实验室, 广西 南宁 530001) 摘要: 为探讨水牛可卡因-苯丙胺调节转录肽 ( cocaine-and amphetamine-regulated transcript, CART) 基因功能, 以水牛卵巢组织 RNA 转录的 cDNA 为模板, 参考 GenBank 预测的水牛 CART 基因组序列设计引物, 采用 PCR 和在线生信分析软件对 CART 基因进行克隆和蛋白质结构预测分 析。 结果显示: 水牛 CART 基因 CDS 区全序列长为 351 bp, 编码 116 个氨基酸; 对所得核酸序列和蛋白质结构分析预测显示, 水牛 CART 基因与 美洲野牛、 野牦牛、 普通牛、 羊、 双峰驼、 单峰驼、 野猪、 马、 非洲野驴、 人类及猕猴相似程度分别为 100%、 98%、 100%、 98%、 91%、 91%、 78%、 85%、 85%、 81%、 90%; Mega 分析显示在多胃反刍动物间有较高同源性; 蛋白质结构预测显示, CART 蛋白含有 N 端信号肽, 是分泌蛋 白, 功能与神经、 辅酶和离子结合相关。 关键词: CART; 水牛; 蛋白质功能预测 中图分类号: S823 83 文献标志码: A 文章编号: 0529-5130(2020)02-0001-06 Cloning, sequencing and bioinformatics analysis of the CART gene in buffalo YU Nongqi, YAN Shengfei, QIN Guangsheng, SHANG Jianghua, ZHENG Haiying, YANG Chunyan ∗ (Guangxi Key Laboratory of Buffalo Genetics, Breeding and Reproduction / Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China) Abstract: In this study, the cocaine-and amphetamine-regulated transcript (CART) gene of buffalos was cloned and its protein structure was predicted. The cDNA transcribed by RNA was extracted from the ovaries of buffalos as templates, primers were designed according to the genomic sequence of CART predicted by GenBank so as to amplify and sequence the CDS region of the CART gene of theanimals. The results showed that the total length of the CDS region of the buffalo CART gene was 351 bp, encoding 116 amino acids. The analysis of nucleic acid sequence and protein structure indicated that the similarities between the buffalo CART gene and that of American bison, wild yak, common cattle, sheep, bactrian camel and single hump camel, wild boar, horse and African wild donkey, human and macaque were 100%, 98%, 100%, 98%, 91%, 91%, 78%, 85%, 85%, 81% and 90%, respectively. The prediction of protein structure indicated that the CART pro⁃ tein contained N-terminal signal peptides confirming it to be a secretory protein. The present study laid a foundation for further research on the function of the CART gene in buffalos. Keywords: CART; buffalo; protein function prediction 可卡因-苯丙胺调节转录肽 ( cocaine - and am⁃ phetamine-regulated transcript, CART) 广泛分布于中 枢神经系统和外周神经系统, 被认为是参与食物摄入 调控的主要无食欲肽[1] 。 CART 是一种关键的神经递 质和激素, 参与调节多种生物过程, 包括食物摄入、 维持体重、 奖励和成瘾、 应激反应、 心理刺激效应和 内分泌功能, 并被报道与家畜繁殖力、 生产性能有 收稿日期: 2019-03-21 ; 修回日期: 2019-12-16 基因项目: 广西科技计划项目 (桂科 AB16380040); 广西水牛研 究所基本科研业务费项目 (水牛基 180503) 第一作者: 于农淇 (1991-): 女, 壮族, 硕士 ∗ 通信作者: 杨春艳, 副研究员, 研究方向为动物遗传育种与繁 殖, E-mail: 86748227@ qq com。 关[2-4] 。 水牛 CART 基因(Gene ID: 102413682,NW_ 005785813)全长为 2 059 bp, 有 3 个 CDS 区, mRNA 全长 775 bp, CDS 全长 312 bp。 前人对水牛 CART 及 其基因的研究不多, 尤其缺少奶水牛上该基因结构、 表达蛋白分析的相关研究。 本研究通过克隆奶水牛 CART 基因 CDS 区序列, 分析研究其编码蛋白质的氨 基酸结构及蛋白定位、 互作, 丰富了 CART 基因数据 信息, 为今后 CART 对水牛内分泌和生长繁殖作用的 研究提供数据支撑。 1 材料与方法 1 1 材料 试验所用水牛卵巢组织样本来自广西南宁屠宰 畜牧与兽医 2020 年 第 52 卷 第 2 期 ·1·
Animal Husbandry Veterinary Medicine 2020 Vol. 52 No. 2 场,卵巢离体后立刻放入液氮中保存,转移至实验室 TRIzol加入0.6mL异丙醇,混匀放置5~10min,4 后置于液氮罐中待用。 ℃12000r/min离心15min。每毫升 trizol加入1 PremiⅸxTq酶、DL2000 DNA marker、反转录试mL75%乙醇,温和震荡使沉淀悬浮,4℃7500r 剂盒、pMDl8-T载体、T4连接酶、DH5α感受态大min离心5min。弃去酒精,所得沉淀室温干燥5~10 肠杆菌购自 TaKara公司,DNA纯化回收试剂盒、质min后用60℃DEPC处理纯净水溶解,紫外分光光 粒小提试剂盒购自天根公司,TRol试剂、氯仿、无度计检测浓度及质量。 水乙醇等试剂购自天地扬公司。 依照反转录试剂盒操作说明操作反转录cDNA。 1.2方法 所得cDNA溶液置于-20℃保存。 1.2.1水牛卵巢组织RNA提取及反转录 1.2.2引物设计与合成 将水牛卵巢组织剪成小块放入研钵,加入液氮反 参考NCBI数据库( Gene id:102413682,NW 复、充分研磨。每50~100mg组织加入1 mL TRIzol,005785813)中预测的水牛CART基因mRNA序列 再用电动匀浆器充分匀浆1~2min,12000r/min离使用NCBI在线引物设计软件设计引物扩增其CDS区 心5min,弃沉淀,每毫升 trizol加入200μ氯仿,序列,引物由北京奥科鼎盛生物科技有限公司合成 手动混匀15s,室温放置10min。4℃1200r/min引物信息见表1。 离心15min,取上层水相转移至新的离心管,每毫升 表1水牛CART基因扩增引物序列 引物名称 引物序列(5′→3′) 退火温度/℃ 扩增产物长度/bp F. GCCTACAGACGGTTGACCC CART-CDS 573 R. TIGCTAAAGCCAGACTCCAGG 1.2.3PCR扩增及载体构建 SwisS-modEl(htTp: //swiss-model. expasy. org/)1es PCR反应体系25μL: Premix Taq12.5μL,软件分析蛋白质理化性质、二级和三级结构,利用 dH2O9.5μ,100ng/μ L CDNA模板1μL,10 ng/ PSORT II prediction在线软件(hts:// psort.hge.jn μL上下游引物各1μL,dH2O15.75μL。PCR反应fom2htm)预测CART蛋白的亚细胞定位, STRING 条件:98℃预变性3min;98℃变性30s,62℃退蛋白质互作数据库检索水牛CART蛋白质序列,在 火30s.72℃延伸90s,共35个循环;72℃延伸6http://kiharalab.org/web/results.phpjob_id=49337 min,4℃保存。PCR产物用2%琼脂糖凝胶电泳检网站上进行蛋白质功能预测。 测,用DNA纯化回收试剂盒切胶回收。将回收的 PCR产物与pMD18-T载体通过T4连接酶过夜连接, 2结果与分析 转化感受态大肠杆菌DH5α,37℃、100r/min培养12.1水牛CART基因的PCR扩增 h,涂于含有Amp、X-gl、IPTG的IB固体培养基 紫外分光光度计检测所得RNA纯度及浓度较好。 上,37℃培养12h后挑取阳性单克隆扩增培养。用以反转录所得ωDNA为模板,用所设计引物扩增结果 质粒小提试剂盒提取质粒后PCR检测鉴定,将检测如图1所示。由图可见,所设计引物PCR扩增产物 阳性的质粒送北京奧科鼎盛生物科技有限公司测序条带清晰且特异性良好,与预期大小相符,可进行下 鉴定 步试验。 1.2.4生物信息学分析 2.2水牛CART基因克隆及测序 CART基因序列比对及同源性使用NCBI在线分 扩增产物切胶回收后转化大肠杆菌,经挑选培养 析,使用MEGA7.0软件构建系统进化树。通过在线后送测序鉴定。菌液PCR产物约570bp大小(图2) 软件 TMHMM2.0(htp:/w.chs.du.dk/ services/为阳性克隆,将扩大繁殖所得菌液送测序公司测序。 TMHMM/)预测分析跨膜结构及信号肽,使用Prot测序结果与NCBI序列对比显示扩增获得序列与 param(https://web.expasy.org/cgi-bin/protparam/pot-CART基因CDS区域一致,与预期相同(图3、图 param)、 Prot scale(htts://web. expasy.og/mp/4)。所获得的水牛CART基因CDS序列全长为351 scores55205.txt)、 SOPMA(htp://npsa- pbil. ibep.f/bp,编码116个氨基酸。 cgi-bin/npsa_automat. pl? page= npsa_sopma. html)
场, 卵巢离体后立刻放入液氮中保存, 转移至实验室 后置于液氮罐中待用。 Premix Taq 酶、 DL2000 DNA marker、 反转录试 剂盒、 pMD18-T 载体、 T4 连接酶、 DH5α 感受态大 肠杆菌购自 TaKaRa 公司, DNA 纯化回收试剂盒、 质 粒小提试剂盒购自天根公司, TRIzol 试剂、 氯仿、 无 水乙醇等试剂购自天地扬公司。 1 2 方法 1 2 1 水牛卵巢组织 RNA 提取及反转录 将水牛卵巢组织剪成小块放入研钵, 加入液氮反 复、 充分研磨。 每 50~100 mg 组织加入 1 mL TRIzol, 再用电动匀浆器充分匀浆 1 ~ 2 min, 12 000 r/ min 离 心 5 min, 弃沉淀, 每毫升 TRIzol 加入 200 μL 氯仿, 手动混匀 15 s, 室温放置 10 min。 4 ℃ 12 000 r/ min 离心 15 min, 取上层水相转移至新的离心管, 每毫升 TRIzol 加入 0 6 mL 异丙醇, 混匀放置 5 ~ 10 min, 4 ℃ 12 000 r/ min 离心 15 min。 每毫升 TRIzol 加入 1 mL 75%乙醇, 温和震荡使沉淀悬浮, 4 ℃ 7 500 r/ min 离心 5 min。 弃去酒精, 所得沉淀室温干燥 5 ~ 10 min 后用 60 ℃ DEPC 处理纯净水溶解, 紫外分光光 度计检测浓度及质量。 依照反转录试剂盒操作说明操作反转录 cDNA。 所得 cDNA 溶液置于-20 ℃保存。 1 2 2 引物设计与合成 参考 NCBI 数据库 ( Gene ID: 102413682, NW _ 005785813) 中预测的水牛 CART 基因 mRNA 序列, 使用 NCBI 在线引物设计软件设计引物扩增其 CDS 区 序列, 引物由北京奥科鼎盛生物科技有限公司合成。 引物信息见表 1。 表 1 水牛 CART 基因扩增引物序列 引物名称 引物序列 (5′ →3′ ) 退火温度/ ℃ 扩增产物长度/ bp CART-CDS F: GCCTACAGACGGTTGACCC R: TTGCTAAAGCCAGACTCCAGG 60 573 1 2 3 PCR 扩增及载体构建 PCR 反 应 体 系 25 μL: Premix Taq 12 5 μL, ddH2O 9 5 μL, 100 ng / μL cDNA 模板 1 μL, 10 ng / μL 上下游引物各 1 μL, ddH2O 15 75 μL。 PCR 反应 条件: 98 ℃ 预变性 3 min; 98 ℃ 变性 30 s, 62 ℃ 退 火 30 s, 72 ℃延伸 90 s, 共 35 个循环; 72 ℃ 延伸 5 min, 4 ℃ 保存。 PCR 产物用 2%琼脂糖凝胶电泳检 测, 用 DNA 纯化回收试剂盒切胶回收。 将回收的 PCR 产物与 pMD18-T 载体通过 T4 连接酶过夜连接, 转化感受态大肠杆菌 DH5α, 37 ℃ 、 100 r/ min 培养 1 h, 涂于含有 Amp、 X-gal、 IPTG 的 LB 固体培养基 上, 37 ℃培养 12 h 后挑取阳性单克隆扩增培养。 用 质粒小提试剂盒提取质粒后 PCR 检测鉴定, 将检测 阳性的质粒送北京奥科鼎盛生物科技有限公司测序 鉴定。 1 2 4 生物信息学分析 CART 基因序列比对及同源性使用 NCBI 在线分 析, 使用 MEGA 7 0 软件构建系统进化树。 通过在线 软件 TMHMM 2 0 ( http: / / www. cbs. dtu. dk / services/ TMHMM/ )预测分析跨膜结构及信号肽, 使用 Prot param(https: / / web.expasy.org / cgi-bin / protparam / prot⁃ param )、 Prot scale ( https: / / web. expasy. org / tmp / scores155205 txt)、SOPMA( http: / / npsa - pbil. ibcp. fr/ cgi- bin / npsa _ automat. pl? page = npsa _ sopma. html)、 SWISS-MODEL( http: / / swiss-model. expasy. org / )在线 软件分析蛋白质理化性质、 二级和三级结构, 利用 PSORT II prediction 在线软件 ( https: / / psort. hgc. jp / form2 html)预测 CART 蛋白的亚细胞定位, STRING 蛋白质互作数据库检索水牛 CART 蛋白质序列, 在 http: / / kiharalab. org / web / results. php? job _ id = 49337 网站上进行蛋白质功能预测。 2 结果与分析 2 1 水牛 CART 基因的 PCR 扩增 紫外分光光度计检测所得 RNA 纯度及浓度较好。 以反转录所得 cDNA 为模板, 用所设计引物扩增结果 如图 1 所示。 由图可见, 所设计引物 PCR 扩增产物 条带清晰且特异性良好, 与预期大小相符, 可进行下 一步试验。 2 2 水牛 CART 基因克隆及测序 扩增产物切胶回收后转化大肠杆菌, 经挑选培养 后送测序鉴定。 菌液 PCR 产物约 570 bp 大小 (图 2) 为阳性克隆, 将扩大繁殖所得菌液送测序公司测序。 测序结果与 NCBI 序列对比显示扩增获得序列与 CART 基因 CDS 区域一致, 与预期相同 (图 3、 图 4)。 所获得的水牛 CART 基因 CDS 序列全长为 351 bp, 编码 116 个氨基酸。 ·2· Animal Husbandry & Veterinary Medicine 2020 Vol 52 No 2
畜牧与兽医2020年第52卷第2期 2000 750 570bp M.DL000 Marker;1~3.PCR产物 图1水牛CART基因PCR产物电泳分析 M.DL2000 Marker;neg-1,neg-2,neg-3.阴性克隆PCR产物 os-1.阳性克隆PCR产物 图2水牛CART基因阳性单克隆PCR产物电泳分析 CGAGTTTCAGCACCATGGAGAGCC CTGCTACCTTTGCTGGGCGCCTTAGCCCA GCCGAGCTCCAGCCCCGAGCCCTGGACATCTACTCCGCCGTGGAGGATGCCTCC G AL A 0 E D A E L 0 PR A L D I FAAGAAGCTCAAGAGTAAACGTATTCCAATTT) YE K KYG 0 v: C AG E0C AV R K GKL 5 C AATTCCTTCCTCCTGAAGTGCTTATGAGGCACCCACGCTTCTCCATAAGCCCCCTCCCCACT NS FL LK C 图3水牛CART基因克隆测序及预测的蛋白编码 起始位点ATG 终止位点TGA 342bp 426bp 858bp 390bp 5°-UTR 外显子1 内含子 外显子2 内含子 外显子3 图4水牛CART基因组序列结构示意图 2.3水牛CART基因同源性及系统进化树分析 与猕猴相似程度为91%,与马、驴相似程度为85% NCBI序列相似性分析(htps:// blast.ncbi.nlm.与人类、野猪相似程度分别为81%、78%。Mega7.0 nih. gov/ Blast.cgi)对所得水牛序列与 GenBank中与美分析结果见图5,由图可知水牛与黄牛、牦牛和美洲 洲野牛、野牦牛、普通牛、羊、双峰驼与单峰驼、野野牛,以及绵羊最早聚为一类,同源性较高,说明 猪、马和驴、人类与猕猴基因序列进行比对分析,使CART基因在多胃反刍动物中保守性较强。单峰驼与 用Mega7.0软件构建系统进化树,参考序列及相似双峰驼、野猪、马与驴各聚为一类,共3条分支,与 性分析结果见表2。所得水牛CART基因序列与美洲牛羊共聚为一大分支,人类与猕猴聚为一类,为另 野牛、普通牛相似程度为100%,与野牦牛、羊相似大分支。聚类分析可见马与驴分化最早,最晚为单峰 程度为哭8%,与双峰驼、单峰驼相似程度为%,驼与双峰驼
M. DL2000 Marker; 1~ 3 PCR 产物 图 1 水牛 CART 基因 PCR 产物电泳分析 M. DL2000 Marker; neg- 1, neg - 2, neg - 3 阴性克隆 PCR 产物; pos-1 阳性克隆 PCR 产物 图 2 水牛 CART 基因阳性单克隆 PCR 产物电泳分析 图 3 水牛 CART 基因克隆测序及预测的蛋白编码 图 4 水牛 CART 基因组序列结构示意图 2 3 水牛 CART 基因同源性及系统进化树分析 NCBI 序列相似性分析 ( https: / / blast. ncbi. nlm. nih.gov / Blast.cgi)对所得水牛序列与 GenBank 中与美 洲野牛、 野牦牛、 普通牛、 羊、 双峰驼与单峰驼、 野 猪、 马和驴、 人类与猕猴基因序列进行比对分析, 使 用 Mega 7 0 软件构建系统进化树, 参考序列及相似 性分析结果见表 2。 所得水牛 CART 基因序列与美洲 野牛、 普通牛相似程度为 100%, 与野牦牛、 羊相似 程度为 98%, 与双峰驼、 单峰驼相似程度为 91%, 与猕猴相似程度为 91%, 与马、 驴相似程度为 85%, 与人类、 野猪相似程度分别为 81%、 78%。 Mega 7 0 分析结果见图 5, 由图可知水牛与黄牛、 牦牛和美洲 野牛, 以及绵羊最早聚为一类, 同源性较高, 说明 CART 基因在多胃反刍动物中保守性较强。 单峰驼与 双峰驼、 野猪、 马与驴各聚为一类, 共 3 条分支, 与 牛羊共聚为一大分支, 人类与猕猴聚为一类, 为另一 大分支。 聚类分析可见马与驴分化最早, 最晚为单峰 驼与双峰驼。 畜牧与兽医 2020 年 第 52 卷 第 2 期 ·3·
Animal Husbandry Veterinary Medicine 2020 Vol. 52 No. 2 表2水牛CART基因与11个物种编码区相似性分析 E值/随机匹配可能性与水牛基因匹配一致性/% 参考序列 美洲野牛( Bison bison bison) 7.00E-107 XM010853569.1 野牦牛( Bos mutus) CP027088.1 双峰驼( Camelus bactrianus) 单峰驼( Camelus dromedarii 1.00E-73 XM_o10974976.1 野猪( Sus scrofa) 马( Equus caballus) XM023618226 驴( Equus asinus) XM014851944.1 人类( Homo sapiens) 猕猴( Macaca mulatta) 2.00E-57 NM001265877.1 agr Bison bison bison 2068(NW 011494901.1 g! os mutus2196N07005930741) Camelus bactrianus 2210(NW 011511773.1) Camelus dromedarius 2107(NW 011591039.1) Sus scrofa 2035(NC 010458. 4) 100usa10N004638051 5CART基因系统进化树 2.4水牛CART蛋白理化性质分析 经 Prot Paran在线软件预测显示,水牛CART基因 由116个氨基酸残基组成,分子式为CsHN1s5O3S8 分子量为12790.18,带负电荷的残基总数(Asp+ Gu)为15个,带正电残基总数(Arg+Lys)为17 个,理论p值为8.22,属于碱性蛋白。消光值在 4845~4470之间,不稳定性指数(Ⅱ)计算为 38.67,属于稳定蛋白。 Prot scale在线预测分析显 示,水牛CART蛋白的总平均亲水性为-0.006,属 于两性氨基酸(图6)。 位点 图6水牛CART基因蛋白疏水性分析
表 2 水牛 CART 基因与 11 个物种编码区相似性分析 物种 E 值/ 随机匹配可能性 与水牛基因匹配一致性/ % 参考序列 美洲野牛(Bison bison bison) 7 00E-107 100 XM_010853569.1 野牦牛(Bos mutus) 0 98 CP027088 1 普通牛(Bos taurus) 0 100 AY603972 1 羊(Ovis aries) 4 00E-99 98 XM_015101190 1 双峰驼(Camelus bactrianus) 7 00E-74 91 XM_010949313 1 单峰驼(Camelus dromedarius) 1 00E-73 91 XM_010974976 1 野猪(Sus scrofa) 0 78 EF581838 1 马(Equus caballus) 2 00E-54 85 XM_023618226 1 驴(Equus asinus) 1 00E-53 85 XM_014851944 1 人类(Homo sapiens) 3 00E-180 81 NG_015988 1 猕猴(Macaca mulatta) 2 00E-57 90 NM_001265877 1 图 5 CART 基因系统进化树 2 4 水牛 CART 蛋白理化性质分析 经 ProtParam 在线软件预测显示,水牛 CART 基因 由 116 个氨基酸残基组成,分子式为 C566H949N155O163 S8 , 分子量为 12 790 18, 带负电荷的残基总数 (Asp + Glu) 为 15 个, 带正电残基总数 (Arg +Lys) 为 17 个, 理论 pI 值为 8 22, 属于碱性蛋白。 消光值在 4 845~ 4 470 之 间, 不 稳 定 性 指 数 ( Ⅱ) 计 算 为 38 67, 属于稳定蛋白。 Prot scale 在线预测分析显 示, 水牛 CART 蛋白的总平均亲水性为- 0 006, 属 于两性氨基酸 (图 6)。 图 6 水牛 CART 基因蛋白疏水性分析 ·4· Animal Husbandry & Veterinary Medicine 2020 Vol 52 No 2
畜牧与兽医2020年第52卷第2期 2.5水牛CART蛋白信号肽及跨膜结构预测分析 2.6水牛CART蛋白二级和三级结构预测及亚细胞 TMHMM2.0预测结果显示,水牛CART蛋白非定位 跨膜蛋白,序列1~30位氨基酸处存在1段N端信号 使用 SOPMA分析显示,水牛CART蛋白二级结 肽,是分泌蛋白(图7) 构主要为α螺旋,无规则卷曲和延伸链,分别占 50.86%、37.93%、10.34%,属于α类蛋白质(图 8)。采用 SWISS-MODEL在线软件预测水牛CART蛋 白三级结构,结果显示该蛋白主要由α-螺旋和无规 则卷曲组成,与二级结构预测结果一致。利用 PSORTⅡ prediction在线软件预测CART蛋白的亚细 胞定位,结果显示该蛋白在在细胞外、线粒体、空 泡、细胞核的分布比例分别为667%、11.1% 跨膜区 膜内区 膜外区 l1.1%、11.1%。 图7水牛CART蛋白跨膜结构预测分析结果 长竖线.α螺旋;中长竖线.延伸链;短竖线.无规则卷曲 图8水牛CART蛋白二级结构预测 2.7蛋白质互作及功能预测 使用 STRING蛋白质互作数据库检索水牛CART 蛋白质序列,检索到该序列与牛可卡因-苯丙胺调解 转录肽前肽( cocaine- and amphetamine- regulated tran- AGRP script peptide concentrations, CARTPT)蛋白质一致, 互作网络结果如下图。由图可见并不能找到有确切证 CDI 据(论文和试验数据)证明与 CARTPT互作的蛋白 CARTPT 质,与之共表达的蛋白有白细胞分化抗原19( cluster of differentiation19,CD19)、神经肽Y( neuropeptide NerSEs Y,NPY)和促肾上腺皮质激素释放激素( corticotro- pin releasing hormone,CRH),其余蛋白质均为文本 关联(图9)。在htp:// kiharalab.og/web/ index.php 网站上进行蛋白质功能预测,预测结果显示,未能预 测出达到置信度的分子功能,找到35条低置信度生 化通路,PFP分数最高的5条为G蛋白偶联受体信 CARTPT:可卡因-苯丙胺调解转录肽前肽;CD19:白细胞分化 号通路、突触传递、细胞饥锇反应、神经冲动传递的抗原19;NPY:神经肽Y;CHH:促肾上腺皮质激素释放激素 正向调节和骨吸收负调节。亚细胞定位预测结果为细 POMC:阿黑皮素原;CCK:胆囊收缩素;LEP:瘦素;AGRP:刺鼠 基因相关蛋白;cR℃3:环状腺嘌呤核糖核苷酸响应蛋白3 胞间隙和分泌泡,与 PSORTⅡ prediction预测结果Nss9:.肿瘤坏死因子受体超家族成员9:MLN:间皮素深色线为 致 共表达,浅色线为基因相邻 图9水牛CART基因蛋白质互作数据库检索结果
2 5 水牛 CART 蛋白信号肽及跨膜结构预测分析 TMHMM 2 0 预测结果显示, 水牛 CART 蛋白非 跨膜蛋白, 序列 1~30 位氨基酸处存在 1 段 N 端信号 肽, 是分泌蛋白 (图 7)。 图 7 水牛 CART 蛋白跨膜结构预测分析结果 2 6 水牛 CART 蛋白二级和三级结构预测及亚细胞 定位 使用 SOPMA 分析显示, 水牛 CART 蛋白二级结 构主要为 α 螺旋, 无规则卷曲和延伸链, 分别占 50 86%、 37 93%、 10 34%, 属于 α 类蛋白质 ( 图 8)。 采用 SWISS-MODEL 在线软件预测水牛 CART 蛋 白三级结构, 结果显示该蛋白主要由 α-螺旋和无规 则 卷 曲 组 成, 与 二 级 结 构 预 测 结 果 一 致。 利 用 PSORT Ⅱ prediction 在线软件预测 CART 蛋白的亚细 胞定位, 结果显示该蛋白在在细胞外、 线粒体、 空 泡、 细 胞 核 的 分 布 比 例 分 别 为 66 7%、 11 1%、 11 1%、 11 1%。 长竖线. α 螺旋; 中长竖线. 延伸链; 短竖线. 无规则卷曲 图 8 水牛 CART 蛋白二级结构预测 2 7 蛋白质互作及功能预测 使用 STRING 蛋白质互作数据库检索水牛 CART 蛋白质序列, 检索到该序列与牛可卡因-苯丙胺调解 转录肽前肽 (cocaine-and amphetamine-regulated tran⁃ script peptide concentrations, CARTPT) 蛋白质一致, 互作网络结果如下图。 由图可见并不能找到有确切证 据 (论文和试验数据) 证明与 CARTPT 互作的蛋白 质, 与之共表达的蛋白有白细胞分化抗原 19 (cluster of differentiation 19, CD19)、 神经肽 Y ( neuropeptide Y, NPY) 和促肾上腺皮质激素释放激素 (corticotro⁃ pin releasing hormone, CRH), 其余蛋白质均为文本 关联 (图 9)。 在 http: / / kiharalab. org / web / index. php 网站上进行蛋白质功能预测, 预测结果显示, 未能预 测出达到置信度的分子功能, 找到 35 条低置信度生 化通路, PFP 分数最高的 5 条为 G 蛋白偶联受体信 号通路、 突触传递、 细胞饥饿反应、 神经冲动传递的 正向调节和骨吸收负调节。 亚细胞定位预测结果为细 胞间隙和分泌泡, 与 PSORT Ⅱ prediction 预测结果 一致。 CARTPT: 可卡因-苯丙胺调解转录肽前肽; CD19: 白细胞分化 抗原 19; NPY: 神 经 肽 Y; CRH: 促 肾 上 腺 皮 质 激 素 释 放 激 素; POMC: 阿黑皮素原; CCK: 胆囊收缩素; LEP: 瘦素; AGRP: 刺鼠 基因 相 关 蛋 白; CRTC3: 环 状 腺 嘌 呤 核 糖 核 苷 酸 响 应 蛋 白 3; TNFRSF9: 肿瘤坏死因子受体超家族成员 9; MSLN: 间皮素深色线为 共表达, 浅色线为基因相邻 图 9 水牛 CART 基因蛋白质互作数据库检索结果 畜牧与兽医 2020 年 第 52 卷 第 2 期 ·5·
Animal Husbandry Veterinary Medicine 2020 Vol. 52 No. 2 3讨论 损伤疾病治疗上具有潜在用途。 水牛具有耐热、耐粗饲和抗病性强等特点,是南 结论 方地区一种重要家畜。近年来随着国外优良品种的引 本试验成功克隆了水牛CART基因CDS区,CDS 进,水牛的经济价值进一步提高,但繁殖力低下问题区全序列长为312bp,编码116个氨基酸。对所得核 直未能解决,是限制水牛产业发展的重要因素。研酸序列和蛋白质结构分析预测显示,水牛CART基因 究和调査水牛繁殖相关基因,从分子层面对其特性进序列与美洲野牛、野牦牛、普通牛、羊、双峰驼与单 行解释,是提高水牛繁殖力的一个重要途径。CART峰驼、野猪、马和非洲野驴、人类与猕猴相似程度分 是一种类生长抑制素多肽,同时兼有神经递质功能。别为100%、%8%、100%、98%、91%、91%、78%、85%、 CART可在中枢神经系统、外周交感神经、感觉或肠85%、81%、90%,说明该基因在反刍动物间有较高同 神经元以及神经内分泌细胞(包括胰岛细胞、胃肠源性,CART蛋白含有N端信号肽,是一种分泌蛋 道黏膜层以及甲状腺和肾上腺髓质)中产生,在下白,功能与神经、辅酶和离子结合相关。 丘脑垂体功能的调节中起着重要的作用,一般被认为 能够抑制食物摄入、刺激能量消耗和调节下丘脑-垂参考文献 体轴,影响动物的生产性状及繁殖性能6。本试验 以广西本地水牛为研究对象,对水牛CART基因[] MAKOWSKA K, RYTEL L, LECH P,etal. Cocaine- and amphet cDNA进行了测序并进行了蛋白质结构及功能预测分 amine-regulated transcript( CART) peptide in the enteric nen 析,在 Kiharalab的预测结果中,PFP分值前二即为 system of the porcine esophagus [J.C R Biol, 2018, 341(6) 325-333 神经肽激素活性和催化活性,其余功能则与构成辅(2]LAUJ, HERZOG H. CART in the regulation of appetite and energ 酶、离子结合相关。但由于 Kiharalab是基于已有研 homeostasis [J]. Front Neurosci, 2014, 8: 313 究进行预测分析,而该基因和蛋白质在家畜上的研究[3]李鹏飞,毕锡麟,王错等.CART在不同发育阶段牛卵泡颗粒细 太少,预测结果并没有达到置信分值,仅能作为一个 胞中的表达和定位[J].中国农业科学,2016,49(12) 参考。 2389-2396 [4]赵栓平,徐磊,贾玉堂.大别山牛CART基因多态性与体尺性状 CART也参与调控颗粒细胞(GCs)卵泡刺激素 相关性研究[J].中国牛业科学,2017,43(5):14-17 信号转导,在与选择优势卵泡相关的雌二醇生产负调 [5] OPONOWICZ A, KOZEOWSKA A, GONKOWSKI S, et al.Chan- 控中具有潜在的特定阶段调节作用。据Ma等 tran- 报道,体外受精妇女的体重指数、GCs中CART script-containing neural structures in the human colon affected by RNA的表达量以及卵泡液中CART肽的水平均呈显 the neoplastic process [. Int J Mol Sci, 2018, 19(2):414. 著正相关,而CART有抑制GC芳香化酶mRNA的表 [6] MORTENSEN A H, CAMPER S A. Cocaine-and amphetamine reg- 达和类固醇生成,导致雌二醇合成降低,卵巢卵泡发 ulated transcript( CART) peptide is expressed in precursor cells and 生改变。Lw等报道了CART基因控制肉牛颗粒细 atotropes of the mouse pituitary gland [. PLos One,2016,11 胞雌二醇表达,说明CART可能与家畜繁殖性能相 (9):c0160068 关。本试验还对cDNA进行了同源性和MEGA分析,7)sENA,Bc0pA,J川 IENEZ-KRASS. Cocaine-ar 结果显示CART基因在反刍动物(牛、羊、水牛) amphetamine-regulated transcript regulation of follicle-stimulating hormone signal transduction in bovine granulosa cells [J].Endocri- 中保守性较强,据此也可以推测CART在这几种动物 nology,2007,148(9):4400-4410. 中功能相近。根据赵栓平等的报道,大别山牛[8]MAx, HAYES E, PRIZANT H,eta.lpin- induced CART( CART基因SNP突变与其体尺性状相关,可以大胆推 ine-and amphetamine-regulated transcript)is a novel intraovarian 测CART与水牛体尺性状、繁殖性能也存在关系。根 mediator of obesity ted infertility in females [ J] 据cDNA进行的蛋白质结构预测分析结果显示,水牛 Endocrinology,2016,157(3):1248-1257 [9 LV L, JIMENEZ-KRASSEL F, SEN A, et al. Evidence supporting CART蛋白属于两性蛋白,其上存在一段N端信号 肽,可能是一种分泌蛋白。在对水牛CART蛋白质功 (CARTPT) in control of granulosa cell estradiol production 能的预测中,预测到其具有与辅酶、离子等结合的功 associated with dominant follicle selection in cattle [J].Biol 能,预测到的相关通路多为神经递质和基础代谢,这 Reprod,2009,81(3):580-586. 说明CABT基因与蛋白质具有多种功能。王路娜[10王路娜,陈一冰,张均等可卡因-苯丙胶调节转录肽对缺血 再灌注模型小鼠皮质突触可塑性的影响[J].中国脑血管病杂 等报道了注射CART可降低缺血-再灌注(L/R) 017,14(3);127-132 小鼠皮质梗死区周围突触的减少,说明CART在神经
3 讨论 水牛具有耐热、 耐粗饲和抗病性强等特点, 是南 方地区一种重要家畜。 近年来随着国外优良品种的引 进, 水牛的经济价值进一步提高, 但繁殖力低下问题 一直未能解决, 是限制水牛产业发展的重要因素。 研 究和调查水牛繁殖相关基因, 从分子层面对其特性进 行解释, 是提高水牛繁殖力的一个重要途径。 CART 是一种类生长抑制素多肽, 同时兼有神经递质功能。 CART 可在中枢神经系统、 外周交感神经、 感觉或肠 神经元以及神经内分泌细胞 (包括胰岛细胞、 胃肠 道黏膜层以及甲状腺和肾上腺髓质) 中产生, 在下 丘脑垂体功能的调节中起着重要的作用, 一般被认为 能够抑制食物摄入、 刺激能量消耗和调节下丘脑-垂 体轴, 影响动物的生产性状及繁殖性能[5-6] 。 本试验 以广西本地水牛为研究对象, 对水牛 CART 基因 cDNA 进行了测序并进行了蛋白质结构及功能预测分 析, 在 Kiharalab 的预测结果中, PFP 分值前二即为 神经肽激素活性和催化活性, 其余功能则与构成辅 酶、 离子结合相关。 但由于 Kiharalab 是基于已有研 究进行预测分析, 而该基因和蛋白质在家畜上的研究 太少, 预测结果并没有达到置信分值, 仅能作为一个 参考。 CART 也参与调控颗粒细胞 (GCs) 卵泡刺激素 信号转导, 在与选择优势卵泡相关的雌二醇生产负调 控中具有潜在的特定阶段调节作用[7-8] 。 据 Ma 等[8] 报道, 体 外 受 精 妇 女 的 体 重 指 数、 GCs 中 CART mRNA 的表达量以及卵泡液中 CART 肽的水平均呈显 著正相关, 而 CART 有抑制 GC 芳香化酶 mRNA 的表 达和类固醇生成, 导致雌二醇合成降低, 卵巢卵泡发 生改变。 Lv 等[9] 报道了 CART 基因控制肉牛颗粒细 胞雌二醇表达, 说明 CART 可能与家畜繁殖性能相 关。 本试验还对 cDNA 进行了同源性和 MEGA 分析, 结果显示 CART 基因在反刍动物 ( 牛、 羊、 水牛) 中保守性较强, 据此也可以推测 CART 在这几种动物 中功能相近。 根据赵栓平等[4] 的报道, 大别山牛 CART 基因 SNP 突变与其体尺性状相关, 可以大胆推 测 CART 与水牛体尺性状、 繁殖性能也存在关系。 根 据 cDNA 进行的蛋白质结构预测分析结果显示, 水牛 CART 蛋白属于两性蛋白, 其上存在一段 N 端信号 肽, 可能是一种分泌蛋白。 在对水牛 CART 蛋白质功 能的预测中, 预测到其具有与辅酶、 离子等结合的功 能, 预测到的相关通路多为神经递质和基础代谢, 这 说明 CART 基因与蛋白质具有多种功能。 王路娜 等[10]报道了注射 CART 可降低缺血-再灌注 ( I/ R) 小鼠皮质梗死区周围突触的减少, 说明 CART 在神经 损伤疾病治疗上具有潜在用途。 4 结论 本试验成功克隆了水牛 CART 基因 CDS 区, CDS 区全序列长为 312 bp, 编码 116 个氨基酸。 对所得核 酸序列和蛋白质结构分析预测显示, 水牛 CART 基因 序列与美洲野牛、 野牦牛、 普通牛、 羊、 双峰驼与单 峰驼、 野猪、 马和非洲野驴、 人类与猕猴相似程度分 别为 100%、98%、100%、98%、91%、91%、78%、85%、 85%、81%、90%, 说明该基因在反刍动物间有较高同 源性, CART 蛋白含有 N 端信号肽, 是一种分泌蛋 白, 功能与神经、 辅酶和离子结合相关。 参考文献: [1] MAKOWSKA K, RYTEL L, LECH P, et al. Cocaine- and amphet⁃ amine-regulated transcript ( CART) peptide in the enteric nervous system of the porcine esophagus [ J]. C R Biol, 2018, 341 ( 6): 325-333. [2] LAU J, HERZOG H. CART in the regulation of appetite and energy homeostasis [J]. Front Neurosci, 2014, 8: 313. [3] 李鹏飞, 毕锡麟, 王锴等. CART 在不同发育阶段牛卵泡颗粒细 胞中的表达和定位 [ J]. 中国农业科学, 2016, 49 ( 12 ): 2389-2396. [4] 赵栓平, 徐磊, 贾玉堂. 大别山牛 CART 基因多态性与体尺性状 相关性研究 [J]. 中国牛业科学, 2017, 43 (5): 14-17. [5] OPONOWICZ A, KOZŁOWSKA A, GONKOWSKI S, et al. Chan⁃ ges in the distribution of cocaine- and amphetamine-regulated tran⁃ script-containing neural structures in the human colon affected by the neoplastic process [ J]. Int J Mol Sci, 2018, 19 ( 2): 414- 419. [6] MORTENSEN A H, CAMPER S A. Cocaine-and amphetamine reg⁃ ulated transcript (CART) peptide is expressed in precursor cells and somatotropes of the mouse pituitary gland [J]. PLoS One, 2016, 11 (9): e0160068. [7] SEN A, BETTEGOWDA A, JIMENEZ-KRASSEL F. Cocaine-and amphetamine- regulated transcript regulation of follicle - stimulating hormone signal transduction in bovine granulosa cells [J]. Endocri⁃ nology, 2007, 148 (9): 4400-4410. [8] MA X, HAYES E, PRIZANT H, et al. Leptin-induced CART (co⁃ caine-and amphetamine-regulated transcript) is a novel intraovarian mediator of obesity - related infertility in females [ J ]. Endocrinology, 2016, 157 (3): 1248-1257. [9] LV L, JIMENEZ-KRASSEL F, SEN A, et al. Evidence supporting a role for cocaine - and amphetamine - regulated transcript ( CARTPT ) in control of granulosa cell estradiol production associated with dominant follicle selection in cattle [ J ]. Biol Reprod, 2009, 81 (3): 580-586. [10] 王路娜, 陈一冰, 张均等. 可卡因-苯丙胺调节转录肽对缺血- 再灌注模型小鼠皮质突触可塑性的影响 [J]. 中国脑血管病杂 志, 2017, 14 (3): 127-132. ·6· Animal Husbandry & Veterinary Medicine 2020 Vol 52 No 2