文章信息
- 于秋寒, 平洪领, 闫允君, 史会来, 吕振明, 杨静文. 2018.
- YU Qiu-han, PING Hong-ling, YAN Yun-jun, SHI Hui-lai, LÜ Zhen-ming, YANG Jing-wen. 2018.
- 刺参GPR54-like受体基因克隆及特征分析
- Molecular cloning and characterization of GPR54-like receptor in the sea cucumber Apostichopus japonicas
- 海洋科学, 42(3): 92-100
- Marine Sciences, 42(3): 92-100.
- http://dx.doi.org/10.11759/hykx20171220001
-
文章历史
- 收稿日期:2017-12-02
- 修回日期:2018-01-14
2. 浙江省海洋水产研究所, 浙江 舟山 316022
2. Marine Fisheries Research Institute of Zhejiang Province, Zhoushan 316022, China
GPR54(G-protein Coupled Receptor 54)是神经肽Kisspeptin的专一受体, 通过配体Kisspeptin激活后介导细胞信号转导进而行使一系列生理功能。Kisspeptin属于RF-酰胺肽家族, 其前体肽由145个氨基酸编码组成, 通过不同位点水解酶切, 可形成片段长度不同的酰胺化短肽, 常见活性肽为: Kisspeptin-54 (Kp54)、Kisspeptin-14(Kp14)、Kisspeptin-13 (Kp13)和Kisspeptin-10(Kp10)等。Lee等[1]1996年首次克隆鉴定了Kisspeptin编码基因KISS-1, 并认为Kisspeptin具有明显抑制肿瘤转移的生理功能; 为进一步研究其功能的调控机理, Lee等[2]开始关注其受体, 并从鼠脑中克隆得到受体GPR54;此后, 人们又鉴定了具有酰胺化末端的Kisspeptin神经肽是受体GPR54的内源性配体[3-5]。在后续的Kisspeptin/GPR54系统生理功能及临床医学研究中, 该系统在生殖发育调控及青春期启动发育等生理过程中的关键作用得到了进一步揭示[6]。
近年来, Kisspeptin/GPR54信号通路的生殖调控研究亦受到了人们的广泛关注并已经建立了完整的理论体系, 目前认为在哺乳动物中其主要作用过程主要有两条: (1)配体Kisspeptin通过激活下丘脑GnRH神经元表面GPR54受体, 介导一系列细胞信号转导, 刺激分泌释放GnRH, 进而通过HPG轴的调控生殖发育启动[7-9]; (2)刺激垂体GTH神经元, 介导细胞信号转导, 调控黄体生成素(LH)和促卵泡素(FSH)分泌, 进而影响生殖行为[10-11]。同时, 研究表明, 不同动物中GPR54所介导的信号转导通路相对保守, 主要通过与Gαq/11及Gαs蛋白偶联两条信号途径[12-13]。然而, 由于进化距离等原因, 有关Kisspeptin/GPR54系统的研究主要集中在脊椎动物中, 低等动物中仅有文昌鱼的Kisspeptin/GPR54系统有报道[14], 其他无脊椎动物中尚未证明该系统的存在, 更无相关生理调控功能的报道。
刺参(Apostichopus japonicus)属棘皮动物门(Echinodermata)、海参纲(Holothuroidea)、楯手目(Aspidochirotida)、刺参科(Stichopodidae)、仿刺参属(Apostichopus), 属于后口动物, 在无脊椎动物中进化地位高[15]。因此相比其他更为低等的无脊椎动物, 刺参中鉴定研究Kisspeptin/GPR54系统较为可行。中国的刺参主要分布于辽宁、山东、河北等北方沿海, 是中国重要的海水养殖经济动物[16]。由于刺参具有再生、夏眠、排脏、自溶等较为独特的生理现象, 其相关基础研究得到了人们的关注, 刺参全基因组测序的完成[17], 为进一步深入探究刺参生物学提供了理论基础。本文利用分子克隆、生物信息学分析、表达载体构建、体外表达等技术, 获得了刺参GPR54-like (Aj GPR54L)编码基因序列全长, 分析了GFP荧光蛋白融合表达AjGPR54L的亚细胞定位以及AjGPR54L在不同组织的转录表达差异, 初步探讨了GPR54受体结构与功能特征, 以期为深入揭示GPR54受体在刺参中的功能活性及调控机制奠定基础。
1 材料与方法 1.1 材料本文所用刺参于2016年4月取自山东青岛胶南刺参养殖场, 采捕刺参时池内水温为12.35℃, 体质量为79.7 g±4.7 g; 低温运至实验室, 暂养驯化一周, 池内水温控制在16℃±0.5℃, 盐度为31.0±0.5, 24 h不间断充气, 每天按时定量投喂人工配合饲料; 活体解剖获取刺参组织(呼吸树、消化道、体壁、肌肉、神经环、雌性性腺), 液氮冷冻保存。
1.2 总RNA提取将上述组织样品从液氮中取出, 液氮研磨后取30 mg样品, 参照Invitrogen TRIzol Reagent说明提取获得的总RNA样品用琼脂糖电泳检查RNA完整性, NanoDrop 2000超微量分光光度计检测RNA浓度, RNA置于-80℃冰箱保存备用。
1.3 刺参GPR54L全长cDNA扩增以提取的RNA为模版, Oligo-dT为逆转录引物, 使用Clontech公司SMARTScribe Reverse Transcriptase 42℃ 1h合成cDNA第一链, 根据转录组基因库中获取的GPR54核心序列信息, 通过序列分析, 使用primer 5.0引物设计软件设计5′RACE引物和3′RACE引物, 与试剂盒中的引物结合通过巢式PCR扩增GPR54基因全长序列的5′端和3′端。扩增后的PCR产物经琼脂糖电泳后, 用DNA回收试剂盒(Axygen)纯化回收后连接到pMD18-T(TaKaRa)载体上, 转化至E.coli DH5α感受态细胞, 经蓝白斑筛选, 挑选阳性克隆, 于LB培养基37℃摇床培养8~12 h后提取质粒, 送生工生物工程(上海)有限公司进行双向测序。所用引物序列见表 1。
引物 | 序列(5′-3′) | 功能 |
5′-outer | TGATAGTGGCCAAGTATCGTTC | AjGPR54L-5′ RACE |
5′-inner | GTCAGGTCCGAACACGCCAA | |
3′-outer | TACCGTCGTTCCACTAGTATTT | AjGPR54L-3′ RACE |
3′-inner | CGCGGATCCTCCACTAGTGATTTCACTATAGG | |
AjGPR54L-vec-F | CCCAAGCTTATGGCTGACAACCTG | AjGPR54L-EGFP 载体构建 |
AjGPR54L-vec-R | CGGGATCCTACAACTGAGAGCTAAGT | |
AjGPR54L-q-F | TACAGAGGTATTTAGCGGATTGACG | AjGPR54L-qRT-PCR |
AjGPR54L-q-R | GGTCTTCATCTTCCTACTGCCATT | |
ACTB-F | AAGGTTATGCTCTTCCTCACGC | 内参基因 |
ACTB-R | GATGTCACGGACGATTTCACG | |
TUBB-F | CACCACGTGGACTCAAAATG | 内参基因 |
TUBB-R | GAAAGCCTTACGACGGAACA |
为了更准确地了解刺参相关受体的序列特征和系统进化关系, 笔者对基因序列以及氨基酸序列进行了分析(表 2), 并将获得的数据进行整理比对。
分析内容功能与作用 | 网络平台或软件名称 |
预测等电点和相对分子质量 | Compute pI/Mw (http://web.expasy.org/compute_pi/) |
查找糖基化位点(N-X-S/T) | NetNGlyc 1.0 Server (http://www.cbs.dtu.dk/services/NetNGlyc/) |
查找磷酸化位点 | NetPhos 2.0 Server (http://www.cbs.dtu.dk/services/NetPhos/) |
预测二硫键位置 | ScanProsite tool (http://prosite.expasy.org/scanprosite/) |
寻找信号肽 | SignalP 4.1 Server (http://www.cbs.dtu.dk/services/SignalP/) |
序列同源性比对分析 | NCBI进行Blast(http://blast.ncbi.nlm.nih.gov/Blast.cgi) |
预测蛋白跨膜结构域 | TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM/) |
预测蛋白功能结构域 | InterPro (http://www.ebi.ac.uk/interpro/scan.html) |
预测蛋白三级结构 | SWISS-MODEL (http://swissmodel.expasy.org/) |
序列多重比较 | BioEdit软件 |
构建系统进化树 | MEGA5.0软件 |
预测氨基酸序列 | DNAMAN8.0 |
基于AjGPR54L基因序列特征和表达载体特征, 使用primer5.0在线软件设计正向引物(AjGPR54L- vec-F)和反向引物(AjGPR54L-vec-R)扩增AjGPR54L CDS。使用限制性内切酶HindШ和BamHI双酶切扩增产物及空白质粒后, 将PCR扩增产物连接到pEGFP-N1表达载体上, 经转化扩增提取质粒、双酶切电泳检测为阳性后, 测序验证载体准确性, 并将构建的载体命名为AjGPR54L-EGFP。
1.6 转染及共聚焦显微镜观察HEK293细胞采用DMEM高糖完全培养基培养(DMEM高糖培养基+10%胎牛血清+1%双抗)。细胞培养在37℃、5% CO2、湿度衡定的恒温培养箱中。使用X-treme GENE HP DNA Transfection Reagent (Roche)按其说明将AjGPR54L-EGFP转染至HEK293细胞, 12~16 h后, 将瞬时转染AjGPR54L-EGFP的HEK293细胞经胰蛋白酶处理后接种于载玻片上。24 h后, 膜荧光探针DiI染膜, 细胞核染料DAPI染核, 多聚甲醛固定后封片, 置于激光扫描共聚焦显微镜(Leica TCS SP5II)下检测荧光信号分布情况。
1.7 实时荧光定量PCR基于AjGPR54L ORF序列设计特异性引物(表 1)。以1.3中制备的cDNA为模板, 选择刺参β-actin (ACTB)和β-tubulin(TUBB)作为内参基因[18], 使用SYBR PrimeScriptTM逆转录试剂盒在ABI 7500 Real Time PCR仪上进行定量扩增, 实验重复3次。标准曲线和Ct值等由ABI7500 Real Time PCR仪自带的ABI7500 Real time PCR System Detection软件自动完成。使用2-△△Ct方法测定不同组织GPR54基因的表达量水平[19-20]。用SPSS 17.0软件进行单因素方差分析(one way ANOVA), Duncan法多重比较检验组间差异的显著性, P < 0.05为差异显著。使用软件Origin 8.1绘制柱形图。
2 结果 2.1 AjGPR54L cDNA全长及预测氨基酸序列克隆获得的AjGPR54L序列全长为1454 bp, 包含993 bp的ORF(编码330个氨基酸)、191 bp的5′UTR和270 bp的3′UTR。预测蛋白质的相对分子质量为37.83 ku、等电点(pI)为9.81。通过预测发现AjGPR54L的氨基酸序列包含: 1个N-连接糖基化位点, 5个ASN-XAA-Ser/Thr序列, 1个细胞外的N端结构域(N9), 34个磷酸化位点, 包括13个丝氨酸(S27、S29、S98、S123、S 227、S230、S263、S268、S272、S294、S307、S322、S328)、16个苏氨酸(T4、T18、T33、T58、T65、T89、T94、T126、T137、T176、T182、T211、T222、T234、T271、T304)和5个酪氨酸(Y7、Y125、Y167、Y301、Y303)位点, 典型的7次跨膜结构, 3个胞外环(EC)和3个胞内环(IC)见图 1和图 2。另外, AjGPR54L具有GPCR视紫红质蛋白家族的典型特征, 在IC3区有1个保守的半胱氨酸, 1个典型的7tm_1结构域, 因此可判断AjGPR54L属于GPCR视紫红质蛋白家族成员。
通过AjGPR54L氨基酸序列与文昌鱼、斑马鱼、小鼠、人等已鉴定活性功能的GPR54氨基酸序列的多序列比对分析发现, AjGPR54L与斑马鱼GPR54a同源性最高, 为30.0%, 而且与胞外N端和C端相比, 7次跨膜结构域相对比较保守(图 2)。利用Predict Protein网络平台预测了AjGPR54L的二级结构及其蛋白结合位点, 同时用SWISS-MODEL结构预测软件, 以Protein Data Bank中的5gli.1.A片段为模板, 加入人工矫正, 预测获得AjGPR54L的三维结构图(图 3)。
2.2 AjGPR54L序列系统进化树为了探究AjGPR54L和其他物种GPR54之间的进化关系, 本研究从Genbank数据库中选取21条不同物种的GPR54氨基酸序列(包括预测信息), 用MEGA 5.0软件以最大似然法构建系统进化树(图 4)。通过BLAST同源搜索, 发现根据cDNA序列预测的刺参AjGPR54L首先与预测的紫球海胆GPR54聚为一支, 再与文昌鱼GPR54L-1聚为一支, 并与脊椎动物GPR54s相距较远。
2.3 AjGPR54L-EGFP在HEK293细胞中的亚细胞定位为分析AjGPR54L亚细胞定位, 本文将构建的AjGPR54L-EGFP的融合表达质粒瞬时转染至HEK293细胞, 24 h后, DiI染膜, DAPI染核, 多聚甲醛固定后封片, 置于激光共聚焦扫描显微镜观察。观察结果显示:细胞成功转染AjGPR54L-EGFP, 可见明显的绿色荧光信号, 且分布在细胞膜上, 与DiI红色荧光重叠(图 5)。这一结果表明HEK293细胞中表达的AjGPR54L受体蛋白定位于细胞膜, AjGPR54L- EGFP蛋白中C-末端EGFP标记不影响AjGPR54L表达。
2.4 AjGPR54L的组织表达用实时荧光定量PCR方法, 以刺参ACTB和TUBB为内参基因, 对基因表达量进行校正, 检测成体刺参GPR54基因在刺参呼吸树、消化道、体壁、肌肉、神经环、雌性性腺中的定量表达, 结果(图 6)表明: AjGPR54L在刺参不同组织中均有表达但表达量存在差异, 在神经环中表达量显著高于其他组织; 在呼吸树、消化道、体壁、肌肉及雌性性腺中表达量差异不显著。
3 讨论本文克隆得到的刺参GPR54 cDNA序列编码一个330个氨基酸的成熟蛋白。对所编码的氨基酸结构域进行分析表明, AjGPR54L具有7个典型的GPCR家族跨膜结构, 属于视紫红质蛋白家族成员。另外, 预测其含有1个潜在的N-连接糖基化位点、5个Asn-Xaa-Ser/Thr序列子和34个磷酸化位点, 这些位点可能参与调节蛋白质运输、定位、功能及信号转导等[21]。进一步的多序列比对分析表明:与其他物种的GPR54蛋白序列相比, AjGPR54L蛋白序列具有一定保守性, 尤其是在跨膜区。但与不同物种GPR54的N端和C端区域相似性较低, 可能与不同物种中对不同内源配体的结合特性不同及胞内信号转导相关因子结构差异相关。
GPCRs是众多生物调控及基础医学研究的重要研究对象, 目前研究已建立了许多GPCR超家族受体蛋白的三维结构、配体结合过程以及胞内信号转导的介导过程等研究模型[22-23]。GPR54属于GPCR超家族受体蛋白, 通过结合激动剂kisspeptin, 进一步激活胞内偶联的G蛋白并介导信号转导[12]。本研究利用生物信息学方法建立了AjGPR54L三维结构模型并分析了它的二级结构, 预测了其蛋白结合位点, 这些结合位点可能对GPR54与Kisspeptin配体结合及胞内信号因子相互作用相关。
G蛋白偶联受体位于细胞膜上, 具有7次跨膜结构, 通过感受细胞外的分子以激活胞内的信号转导通路。AjGPR54L的细胞膜定位对其功能活动至关重要, 本文发现其主要在细胞膜表达, 这表明所表达的AjGPR54L-EGFP细胞膜定位特性与7次跨膜的蛋白质三维结构特征一致, 且证明其C-末端融合表达的EGFP不影响AjGPR54L的表达及亚细胞定位。从而为进一步查明该受体与其配体之间的相互作用、以及阐明信号转导与调控机制提供了基础。
在脊椎动物中, Kisspeptin/GPR54系统主要行使生殖调控等生理功能[6, 24-25]。人体的青春发育启动直接受到Kisspeptin/GPR54系统的HPG轴的调控影响[26], 而GPR54受体基因表达主要分布于多个神经组织, 以接受Kisspeptin神经元分泌调控。但研究亦发现GPR54受体基因在胰脏等外周组织中亦有较高的表达分布[3], 说明Kisspeptin/GPR54系统亦可能参与了其他生理功能的调控, 这种GPR54受体基因的表达分布模式在鼠等其他哺乳动物中亦有类似发现[27-28]。本文通过qRT-PCR分析, 表明AjGPR54L在刺参多组织中广泛表达并存在明显差异, 神经环中表达显著高于其他组织。神经环是刺参的重要神经器官, AjGPR54L在神经环中的高水平表达, 与已报道的脊椎动物中GPR54的表达有相似之处, 但该基因是否参与某些神经元调控并通过神经内分泌系统调节某些生理过程, 尚有待进一步研究。
[1] |
Lee J H, Miele M E, Hicks D J, et al. KiSS-1, a novel human malignant melanoma metastasis-suppressor gene[J]. J Natl Cancer Inst, 1996, 88(23): 1731-1737. DOI:10.1093/jnci/88.23.1731 |
[2] |
Lee D K, Nguyen T, O'Neill G P, et al. Discovery of a receptor related to the galanin receptors[J]. FEBS Lett, 1999, 446: 103-107. DOI:10.1016/S0014-5793(99)00009-5 |
[3] |
Kotani M, Detheux M, Vandenbogaerde A, et al. The Metastasis Suppressor Gene KiSS-1 Encodes Kisspeptins, the Natural Ligands of the Orphan G Protein-coupled Receptor GPR54[J]. Journal of Biological Chemistry, 2001, 276(37): 34631-34636. DOI:10.1074/jbc.M104847200 |
[4] |
Muir A I, Chamberlain L, Elshourbagy N A, et al. AXOR12, a Novel Human G Protein-coupled Receptor, Activated by the Peptide KiSS-1*[J]. Journal of Biological Chemistry, 2001, 276(31): 28969-28975. DOI:10.1074/jbc.M102743200 |
[5] |
Ohtaki T, Shintani Y, Honda S, et al. Metastasis suppressor gene KiSS-1 encodes peptide ligand of a G-protein-coupled receptor[J]. Nature, 2001, 411(6837): 613. DOI:10.1038/35079135 |
[6] |
De R N, Genin E, Carel J C, et al. Hypogonadotropic hypogonadism due to loss of function of the KiSS1-derived peptide receptor GPR54[J]. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(19): 10972-10976. DOI:10.1073/pnas.1834399100 |
[7] |
Thompson E L, Murphy K G. Modulation of the hypothalamic-pituitary-gonadal axis by selective ligands of the KISS1R[J]. Current Opinion in Investigational Drugs, 2010, 11(4): 432-439. |
[8] |
Javed Z, Qamar U, Sathyapalan T. The role of kisspeptin signalling in the hypothalamic-pituitary-gonadal axis-current perspective[J]. Endokrynologia Polska, 2015, 66(6): 534. DOI:10.5603/EP.2015.0066 |
[9] |
Gottsch M L, Cunningham M J, Smith J T, et al. A role for kisspeptins in the regulation of gonadotropin secretion in the mouse[J]. Endocrinology, 2004, 145(9): 4073-4077. DOI:10.1210/en.2004-0431 |
[10] |
Saito H, Sawada T, Yaegashi T, et al. Kisspeptin-10 stimulates the release of luteinizing hormone and testosterone in pre-and post-pubertal male goats[J]. Animal Science Journal, 2012, 83(6): 487-492. DOI:10.1111/j.1740-0929.2011.00978.x |
[11] |
Gopurappilly R, Ogawa S, Parhar I S. Functional Significance of GnRH and Kisspeptin, and Their Cognate Receptors in Teleost Reproduction[J]. Frontiers in Endocrinology, 2013, 4(24): 24. |
[12] |
Castaño J P, Martínez-Fuentes A J, Gutiérrez-Pascual E, et al. Intracellular signaling pathways activated by kisspeptins through GPR54:do multiple signals underlie function diversity[J]. Peptides, 2009, 30(1): 10-15. DOI:10.1016/j.peptides.2008.07.025 |
[13] |
Pampillo M, Camuso N, Taylor J E, et al. Regulation of GPR54 signaling by GRK2 and beta-arrestin[J]. Molecular Endocrinology, 2009, 23(12): 2060. DOI:10.1210/me.2009-0013 |
[14] |
Wang P, Wang M, Ji G, et al. Demonstration of a Functional Kisspeptin/Kisspeptin Receptor System in Amphioxus With Implications for Origin of Neuroendocrine Regulation[J]. Endocrinology, 2017, 158(5): 1461-1473. DOI:10.1210/en.2016-1848 |
[15] |
廖玉麟. 中国动物志, 棘皮动物门, 海参纲[M]. 北京: 科学出版社, 1997. Liao Yulin. Zoology, Echinoderma, Holothuroidea[M]. Beijing: Science Press, 1997. |
[16] |
常亚青, 隋锡林, 李俊. 刺参增养殖业现状、存在问题与展望[J]. 水产科学, 2006, 254(4): 198-201. Chang Yaqing, Sui Xilin, Li Jun. Current situation, existing problems and prospects of Apostichopus japonicus breeding[J]. Aquatic Science, 2006, 254(4): 198-201. |
[17] |
Zhang X, Sun L, Yuan J, et al. The sea cucumber genome provides insights into morphological evolution and visceral regeneration[J]. Plos Biology, 2017, 15(10): e2003790. DOI:10.1371/journal.pbio.2003790 |
[18] |
Wang T, Yang Z, Zhou N, et al. Identification and functional characterisation of 5-HT4 receptor in sea cucumber Apostichopus japonicus(Selenka)[J]. Scientific Reports, 2017, 7: 40247. DOI:10.1038/srep40247 |
[19] |
Livak K J, Schmittgen T D, 2001. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(Delta DeltaC)method[J]. Methods, 2001, 25(4): 402-408. DOI:10.1006/meth.2001.1262 |
[20] |
Breukelen F V, Krumschnabel G, Podrabsky J E. Vertebrate cell death in energy-limited conditions and how to avoid it:what we might learn from mammalian hibernators and other stress-tolerant vertebrates[J]. Apoptosis An International Journal on Programmed Cell Death, 2010, 15(3): 386. DOI:10.1007/s10495-010-0467-y |
[21] |
Arnold J N, Wormald M R, Sim R B, et al. The impact of glycosylation on the biological function and structure of human immunoglobulins[J]. Annual Review of Immunology, 2007, 25(25): 21. |
[22] |
Kristiansen K. Molecular mechanisms of ligand binding, signaling, and regulation within the superfamily of G-protein-coupled receptors:molecular modeling and mutagenesis approaches to receptor structure and function[J]. Pharmacology & Therapeutics, 2004, 103(1): 21-80. |
[23] |
Katritch V, Cherezov V, Stevens R C. Structure-function of the G protein-coupled receptor superfamily[J]. Annual Review of Pharmacology and Toxicology, 2013, 53: 531-536. DOI:10.1146/annurev-pharmtox-032112-135923 |
[24] |
Saito H, Sawada T, Yaegashi T, et al. Kisspeptin-10 stimulates the release of luteinizing hormone and testosterone in pre-and post-pubertal male goats[J]. Animal Science Journal, 2012, 83(6): 487-492. DOI:10.1111/j.1740-0929.2011.00978.x |
[25] |
Funes S, Hedrick J A, Vassileva G, et al. The KiSS-1 receptor GPR54 is essential for the development of the murine reproductive system[J]. Biochemical & Biophysical Research Communications, 2003, 312(4): 1357-1363. |
[26] |
Roux N D, Genin E, Carel J, et al. Hypogonadotropic hypogonadism due to loss of function of the KiSS1-derived peptide receptor GPR54[J]. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(19): 10972-10976. DOI:10.1073/pnas.1834399100 |
[27] |
Hauge-Evans A C, Richardson C C, Milne H M, et al. A role for kisspeptin in islet function[J]. Diabetologia, 2006, 49(9): 2131-2135. DOI:10.1007/s00125-006-0343-z |
[28] |
Dudek M, Kołodziejski P A, Pruszyńska-Oszmałek E, et al. Effects of high-fat diet-induced obesity and diabetes on Kiss1 and GPR54 expression in the hypothalamic-pituitary-gonadal (HPG) axis and peripheral organs (fat, pancreas and liver) in male rats[J]. Neuropeptides, 2016, 56: 41-49. DOI:10.1016/j.npep.2016.01.005 |