留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建

李妙珍 李奕勋 陈静 张磊 廖承红 韩谦

李妙珍,李奕勋,陈静,等. 伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建[J]. 热带生物学报,2021, 12(3):347−355. DOI:10.15886/j.cnki.rdswxb.2021.03.011 doi:  10.15886/j.cnki.rdswxb.2021.03.011
引用本文: 李妙珍,李奕勋,陈静,等. 伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建[J]. 热带生物学报,2021, 12(3):347−355. DOI:10.15886/j.cnki.rdswxb.2021.03.011 doi:  10.15886/j.cnki.rdswxb.2021.03.011
LI Miaozhen, LI Yixun, CHEN Jing, ZHANG Lei, LIAO Chenghong, HAN Qian. Bioinformatics Analysis and Construction of Spatiotemporal Expression Profile of the Serotonin Receptor Protein Family in Aedes aegypti[J]. Journal of Tropical Biology, 2021, 12(3): 347-355. doi: 10.15886/j.cnki.rdswxb.2021.03.011
Citation: LI Miaozhen, LI Yixun, CHEN Jing, ZHANG Lei, LIAO Chenghong, HAN Qian. Bioinformatics Analysis and Construction of Spatiotemporal Expression Profile of the Serotonin Receptor Protein Family in Aedes aegypti[J]. Journal of Tropical Biology, 2021, 12(3): 347-355. doi: 10.15886/j.cnki.rdswxb.2021.03.011

伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建

doi: 10.15886/j.cnki.rdswxb.2021.03.011
基金项目: 国家自然科学基金(31960703,31860702)
详细信息
    第一作者:

    李妙珍(1995−),女. 海南大学生命科学与药学院2018级硕士研究生. E-mail:limiaozhen4956@163.com

    通信作者:

    廖承红(1977−),女,教授. 研究方向:病媒生物学. E-mail:liaochh@hainanu.edu.cn

    韩谦(1963−),男,教授. 研究方向:病媒生物学. Email:qianhan@hainanu.edu.cn

  • 中图分类号: R 384.1

Bioinformatics Analysis and Construction of Spatiotemporal Expression Profile of the Serotonin Receptor Protein Family in Aedes aegypti

  • 摘要: 5-羟色胺(5-HT)可通过多个受体介导来调控昆虫的生长发育和运动等生理活动,也是新型杀虫剂潜在的靶点。本研究对埃及伊蚊的5-羟色胺受体家族进行系统分析,即利用生物信息学方法、进化分析5-HT受体在物种间的同源性,构建系统进化树;使用同源建模的方法以人类5-羟色胺受体结构为模板,构建埃及伊蚊(Aedes aegypti)5-HT受体的三维结构模型;通过实时荧光定量PCR方法研究受体在埃及伊蚊幼虫期、蛹期和成蚊,未吸血雌蚊的头、胸、腹及附肢的表达情况。从NCBI和VectorBase数据库中得到埃及伊蚊的6个5-HT受体基因,生物信息学分析发现6个受体结构均具有7次跨膜区,其中在受体蛋白的二级结构中α-螺旋占有较高的比重。进化分析结果表明,埃及伊蚊5-HT受体基因分别位于3个不同的分支,各分支均与黑腹果蝇5-HT相应受体具有较高的同源性。表达谱构建结果表明,相比各龄幼虫和雌蚊,5-HT受体均在雄蚊中的表达量最高,且在雌蚊各组织中呈现明显的差异。5-HTR1B、5-HTR2B、5-HTR7B在头、胸、腹和附肢4种组织中均有表达;5-HTR1A在头、胸、腹中均有表达;5-HTR1D在胸腔中高表达,但在头部和附肢中不表达;5-HTR7A在附肢中表达丰度最高,而在头部无表达;该结果预示不同受体介导5-HT的功能可能不同。本研究明确了埃及伊蚊5-HT受体家族蛋白的序列特征、进化关系及组织表达特征,为埃及伊蚊5-HT受体功能的进一步研究奠定了基础。
  • 图  2  埃及伊蚊5-HTR家族保守结构域序列比对

    使用Clustal Omega对氨基酸进行多序列比对。“*”为100%同源,“:”为90%同源,“.”为80%同源。

    Fig.  2  Alignment of conserved domains of 5-HT proteins family in A.aegypti

    Multiple sequence alignments of amino acids were performed using Clustal Omega. ‘*’ indicates 100% homology; ‘:’ indicates 90% homology; ‘.’ indicates 80% homology.

    图  3  埃及伊蚊5-HT受体蛋白系统进化树(Boostrap值为1000)

    Drosophila melanogaster;黑腹果蝇;Rhagoletispomonella:苹果实蝇;Gryllusbimaculatus:双斑蟋;Rhipicephalus microplus:微小扇头蜱;Nephotettixcincticeps:黑尾叶蝉;Antheraeapernyi:柞蚕;Triboliumcastaneum:赤拟谷盗;Anoplophoraglabripennis:光肩星天牛;Culex quinquefasciatus:致倦库蚊;Melanaphis sacchari:高粱蚜;Anopheles albimanus:淡色按蚊;Aedes aegypti :埃及伊蚊。

    Fig.  3  Phylogenetic tree of 5-HT receptor proteins in A. aegypti

    图  4  埃及伊蚊5-HT受体不同发育时期表达谱

    A~F:分别代表5-HT受体6个亚基在不同发育时期的相对表达水平(*P<0.05,**P<0.01),****P<0.0001)。

    Fig.  4  Expression profile of 5-HT receptor in A. aegypti at different developmental stages

    A-F represent the relative expression levels of the 6 subunits of 5-HT receptor at different developmental stages (*P<0.05, **P<0.01),****P<0.0001).

    图  5  埃及伊蚊5-HT 受体不同组织表达谱

    A~F:分别代表5-HT受体6个亚基的在雌性成蚊不同组织中的相对表达水平 (*P<0.05, **P<0.01,***P<0.001)。

    Fig.  5  Expression profile of 5-HT receptor in A. aegypti at different tissues

    A-F: Relative expression levels of the 6 subunits of the 5-HT receptors n different tissues of female adult mosquitoes (*P<0.05, **P<0.01,***P<0.001).

    表  1  荧光定量PCR引物

    Table  1  Primers for fluorescence quantitative PCR

    基因名称 Gene引物序列 Sequence of primer(5′-3′)
    上游引物 F下游引物 R
    5-HTR1B TACTGGCGACGGTTATAGG CGACAATATCAGGTGGTTGG
    5-HTR1A CCCACGAGCAATAATACCC CCTTGCGTTTGGAGACC
    5-HTR1D GTTCTACAGTAAATCGTCCTACC CCCGTAGCAATACATCAGC
    5-HTR2B CGAATCTAATTGCTCCTCTCC CTGGACACTTCTGCTATTCC
    5-HTR7A ACGAAACCGCTAGAATACG GTGGCAGCGATATACAGG
    5-HTR7B TTCCTAGAGTACCTCCTAACC GAGTAGAGAGCTTGTTGACG
    RPS17 AAGAAGTGGCCATCATTCCA GGTCTCCGGGTCGACTTC
    下载: 导出CSV

    表  2  埃及伊蚊5-HT受体蛋白二级结构的主要元件比例

    Table  2  The proportion of major secondary protein structure of 5-HT receptor protein in A. aegypti.

    受体名称
    Receptor
    α-螺旋/%
    α-helix (%)
    延伸链/%
    Extended
    strand (%)
    β-转角/%
    Beta
    turn (%)
    无规卷曲/%
    Random
    coil (%)
    Aae5-HTR1B40.3812.741.9244.95
    Aae5-HTR1A32.0715.523.2849.14
    Aae5-HTR2B37.4119.174.9636.46
    Aae5-HTR1D28.09621.736.7143.46
    Aae5-HTR7A47.6713.562.9735.81
    Aae5-HTR7B47.6417.673.2331.47
    下载: 导出CSV
  • [1] RUT V,. HELEEN V, JOZEF VANDEN B Serotonin, serotonin receptors and their actions in insects [J]. Neurotransmitter, 2015, 2(1): e314.
    [2] SAWIN E R, RANGANATHAN R, HORVITZ H R. Caenorhabditis elegans locomotory rate is modulated by the environment through a dopaminergic pathway and by experience through a serotonergic pathway [J]. Neuron, 2000, 26(3): 619 − 631. doi:  10.1016/S0896-6273(00)81199-X
    [3] PRASAD P, OGAWA S, PARHAR IS. Role of serotonin in fish reproduction [J]. Front Neurosci, 2015(9): 195.
    [4] LIU Y, LUO J, CARLSSON M A, et al. Serotonin and insulin-like peptides modulate leucokinin-producing neurons that affect feeding and water homeostasis in Drosophila [J]. J Comp Neurol, 2015, 523(12): 1840 − 1863. doi:  10.1002/cne.23768
    [5] DONALDSON Z R, NAUTIYAL K M, AHMARI S E, et al. Genetic approaches for understanding the role of serotonin receptors in mood and behavior [J]. Curr Opin Neurobiol, 2013, 23(3): 399 − 406. doi:  10.1016/j.conb.2013.01.011
    [6] SOOD S, MORRISON JL, LIU H, et al. Role of endogenous serotonin in modulating genioglossus muscle activity in awake and sleeping rats [J]. Am J Respir Crit Care Med, 2005, 172(10): 1338 − 1347. doi:  10.1164/rccm.200502-258OC
    [7] HALFORD J C, HARROLD J A, LAWTON C L, et al. Serotonin (5-HT) drugs: effects on appetite expression and use for the treatment of obesity [J]. Curr Drug Targets, 2005, 6(2): 201 − 213. doi:  10.2174/1389450053174550
    [8] QI Y X, HUANG J, LI M Q, et al. Serotonin modulates insect hemocyte phagocytosis via two different serotonin receptors [J]. Elife, 2016(5): e12241.
    [9] FRENCH A S, SIMCOCK K L, ROLKE D, et al. The role of serotonin in feeding and gut contractions in the honeybee [J]. J Insect Physiol, 2014, 61: 8 − 15. doi:  10.1016/j.jinsphys.2013.12.005
    [10] LING L, RAIKHEL A S. Serotonin signaling regulates insulin-like peptides for growth, reproduction, and metabolism in the disease vector Aedes aegypti [J]. Proc Natl Acad Sci U S A, 2018, 115(42): E9822 − E9831. doi:  10.1073/pnas.1808243115
    [11] CAI M, LI Z, FAN F, et al. Design and synthesis of novel insecticides based on the serotonergic ligand 1-[(4-Aminophenyl)ethyl]-4-[3-(trifluoromethyl)phenyl]piperazine (PAPP) [J]. Journal of Agricultural and Food Chemistry, 2010, 58(5): 2624 − 2629. doi:  10.1021/jf902640u
    [12] SHARAN S, HILL C A. Advances in agrochemicals: ion channels and G Protein-Coupled Receptors (GPCRs) as targets for pest control[M]. America: American Chemical Society, 2017: 55 − 84.
    [13] LAN J, WANG Z, CHEN Z, et al. Identification of the Aedes aegyptin AChR gene family and molecular target of Spinosad [J]. Pest Management Science, 2021, 77(4): 1633 − 1641. doi:  10.1002/ps.6183
    [14] BHATT S, GETHING P W, BRADY O J, et al. The global distribution and burden of dengue [J]. Nature, 2013, 496(7446): 504 − 507. doi:  10.1038/nature12060
    [15] YIN W, ZHOU X E, YANG D, et al. Crystal structure of the human 5-HT1B serotonin receptor bound to an inverse agonist [J]. Cell Discovery, 2018, 4(1): 12. doi:  10.1038/s41421-018-0009-2
    [16] 王子贺, 兰坚强, 廖承红, 等. 天维菌素处理埃及伊蚊幼虫的转录组学分析[J]. 热带生物学报, 2020, 11(3): 281 − 287.
    [17] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2 −ΔΔ C T method [J]. Methods, 2001, 25(4): 402 − 408. doi:  10.1006/meth.2001.1262
    [18] ALMAGRO ARMENTEROS J J, TSIRIGOS K D, SØNDERBY C K, et al. Signal P 5.0 improves signal peptide predictions using deep neural networks [J]. Nature Biotechnology, 2019(37): 420 − 423.
    [19] BERMAN H M, BATTISTUZ T F, BHAT T N, et al. The protein data Bank [J]. Nucleic Acids Research, 2000, 28(1): 235 − 242. doi:  10.1093/nar/28.1.235
    [20] BERMAN H, HENRICK K F, NAKAMURA H, et al. The worldwide Protein Data Bank (wwPDB): ensuring a single, uniform archive of PDB data [J]. Nucleic Acids Research, 2007, 35(1): D301 − D303.
    [21] LASKOWSKI R A MMW, MOSS D S, THORNTON J M. PROCHECK - a program to check the stereochemical quality of protein structures [J]. J App Cryst, 1993, 26: 283 − 291. doi:  10.1107/S0021889892009944
    [22] MORRIS A L, MACARTHUR M W, HUTCHINSON E G, at al. Stereochemical quality of protein structure coordinates [J]. Proteins, 1992, 12(4): 345 − 364. doi:  10.1002/prot.340120407
    [23] WATERHOUSE A, BERTONI M, BIENERT S, et al. SWISS-MODEL: homology modelling of protein structures and complexes [J]. Nucleic Acids Research, 2018, 46(W1): W296 − W303. doi:  10.1093/nar/gky427
    [24] BIENERT S, WATERHOUSE A, DE BEER TJAART A P, et al. The SWISS-MODEL Repository—new features and functionality [J]. Nucleic Acids Research, 2017, 45(D1): D313 − D319. doi:  10.1093/nar/gkw1132
    [25] STUDER G, REMPFER C, WATERHOUSE A M, et al. QMEANDisCo—distance constraints applied on model quality estimation [J]. Bioinformatics, 2020, 36(6): 1765 − 1771. doi:  10.1093/bioinformatics/btz828
    [26] BENKERT P, BIASINI M, SCHWEDE T. Toward the estimation of the absolute quality of individual protein structure models [J]. Bioinformatics, 2011, 27(3): 343 − 350. doi:  10.1093/bioinformatics/btq662
    [27] BERTONI M, KIEFER F, BIASINI M, et al. Modeling protein quaternary structure of homo- and hetero-oligomers beyond binary interactions by homology [J]. Scientific Reports, 2017, 7(1): 10480. doi:  10.1038/s41598-017-09654-8
    [28] MADEIRA F, PARK YM, LEE J, et al. The EMBL-EBI search and sequence analysis tools APIs in 2019 [J]. Nucleic Acids Research, 2019, 47(W1): W636 − W641. doi:  10.1093/nar/gkz268
    [29] 陈世骧. 昆虫的变态类型与分类体系[J]. 昆虫学报, 1962(1): 1 − 15.
    [30] 李海银, 李艳, 陈曦, 等. 家蚕5-羟色胺受体基因的克隆及表达分析[J]. 中国农业科学, 2015, 48(5): 987 − 1001. doi:  10.3864/j.issn.0578-1752.2015.05.17
    [31] BORYCZ J, BORYCZ JA, KUBOW A, et al. Drosophila ABC transporter mutants white, brown and scarlet have altered contents and distribution of biogenic amines in the brain [J]. J Exp Biol, 2008, 211(Pt 21): 3454 − 3466.
    [32] JOHNSON O, BECNEL J, NICHOLS C D. Serotonin receptor activity is necessary for olfactory learning and memory in Drosophila melanogaster [J]. Neuroscience, 2011(192): 372 − 381.
    [33] HUSER A, ESCHMENT M, GÜLLÜ N, et al. Anatomy and behavioral function of serotonin receptors in Drosophila melanogaster larvae [J]. PLOS ONE, 2017, 12(8): e0181865. doi:  10.1371/journal.pone.0181865
    [34] MATTHEWS B J, YOUNGER M A, VOSSHALL L B. The ion channel ppk301 controls freshwater egg-laying in the mosquito Aedes aegypti [J]. Elife, 2019(8): e43963.
  • [1] 黄兹宝, 陈咪咪, 张勇, 张小坡, 董琳, 张彩云.  无根藤中阿朴菲类生物碱的化学成分分析 . 热带生物学报, 2024, 15(2): 251-258. doi: 10.15886/j.cnki.rdswxb.20230043
    [2] 李佳雪, 丁一, 王猛, 李涛, 郭攀阳, 刘成立, 韦双双, 黄家权, 李洪立, 胡文斌, 汤华.  火龙果HubHLH基因家族的全基因组分析及其对冬季补光诱导开花的表达响应 . 热带生物学报, 2024, 15(2): 198-209. doi: 10.15886/j.cnki.rdswxb.20220108
    [3] 董晨, 郑雪文, 王弋, 全振炫, 李伟才.  荔枝GRF基因家族的全基因组鉴定及表达分析 . 热带生物学报, 2024, 15(2): 190-197. doi: 10.15886/j.cnki.rdswxb.20230099
    [4] 金海峰, 王朝政, 侯清芳, 咸利民, 张华剑, 吴少英.  普通大蓟马MuRhodopsin基因的全长克隆及生物信息学分析 . 热带生物学报, 2023, 14(6): 651-659. doi: 10.15886/j.cnki.rdswxb.20220105
    [5] 翁可欣, 张明亮, 李力.  微生物合成5-羟基色氨酸的研究进展 . 热带生物学报, 2023, 14(1): 42-49. doi: 10.15886/j.cnki.rdswxb.2023.01.016
    [6] 全飞, 兰国玉, 魏亚情, 李明美, 孙树晴, 杜昊楠.  基于宏基因组测序的橡胶林病原微生物的分析 . 热带生物学报, 2022, 13(1): 27-35. doi: 10.15886/j.cnki.rdswxb.2022.01.005
    [7] 符丹凤, 刘洪涛, 杨明秋, 何玉贵.  斑节对虾PmML1PmML2基因的鉴定及表达分析 . 热带生物学报, 2022, 13(5): 440-450. doi: 10.15886/j.cnki.rdswxb.2022.05.003
    [8] 刘应帅, 余瑞, 郑彬彬, 刘嘉慧, 宋奇, 陈荣昊, 严哲.  海南岛森林植被NEP季节性时空变化规律及气候驱动因素分析 . 热带生物学报, 2022, 13(2): 166-176. doi: 10.15886/j.cnki.rdswxb.2022.02.008
    [9] 李光伟, 邢峰华, 敖杰, 毛志远.  海南岛30年空中云水资源时空特征分析 . 热带生物学报, 2022, 13(4): 331-338. doi: 10.15886/j.cnki.rdswxb.2022.04.003
    [10] 代霖欣, 周小磊, 杨术, 侯丹清, 孙成波.  高位处理水精养虾池水体浮游生物的时空变化 . 热带生物学报, 2022, 13(5): 429-439. doi: 10.15886/j.cnki.rdswxb.2022.05.002
    [11] 符传博, 佟金鹤, 徐文帅, 刘丽君.  海南岛臭氧污染时空特征及其成因分析 . 热带生物学报, 2022, 13(4): 404-409. doi: 10.15886/j.cnki.rdswxb.2022.04.012
    [12] 武映宏, 黄东益, 王永, 张大鹏, 刘蕊, 孙熹微, 肖勇.  椰子ALDH基因家族的鉴定及生物信息学分析 . 热带生物学报, 2022, 13(6): 541-549. doi: 10.15886/j.cnki.rdswxb.2022.06.001
    [13] 张忠辉, 邓渊, 强奇, 陈日东, 杨君, 刘贤青, 王守创, 袁潜华, 罗杰, 罗越华.  基于高效液相色谱-质谱联用技术解析山栏稻营养品质的代谢物谱 . 热带生物学报, 2021, 12(4): 419-427. doi: 10.15886/j.cnki.rdswxb.2021.04.003
    [14] 费小雯, 张阳, 李亚军, 邓晓东.  伊蚊3HKT基因RNAi载体构建及口服对伊蚊的致死作用 . 热带生物学报, 2021, 12(3): 356-362. doi: 10.15886/j.cnki.rdswxb.2021.03.012
    [15] 董晨, 魏永赞, 王弋, 郑雪文, 李伟才.  转录组荔枝Dof基因家族的鉴定及其表达 . 热带生物学报, 2021, 12(1): 7-14. doi: 10.15886/j.cnki.rdswxb.2021.01.002
    [16] 胡发广, 毕晓菲, 付兴飞, 李贵平, 李亚男, 吕玉兰, 黄家雄.  天牛信息素的研究进展 . 热带生物学报, 2020, 11(4): 517-522. doi: 10.15886/j.cnki.rdswxb.2020.04.015
    [17] 闫冰玉, 巩笑笑, 谭玉荣, 王丹, 高璇, 张恒, 李双江, 王鹏, 潘英文, 刘进平.  文心兰OnRR10基因的克隆及表达分析 . 热带生物学报, 2020, 11(3): 288-295. doi: 10.15886/j.cnki.rdswxb.2020.03.005
    [18] 宋娜, 张刘宁颖, 曹敏, 郭文雅, 洪雨慧, 吴金山, 陈银华, 于晓惠.  木薯ERF转录因子调控的靶基因筛选与表达分析 . 热带生物学报, 2020, 11(2): 177-189. doi: 10.15886/j.cnki.rdswxb.2020.02.008
    [19] 白雪杨, 陈秀珍, 黄天帆, 江行玉, 周扬.  盐胁迫下拟南芥SCAMP基因克隆和表达的生物信息学分析 . 热带生物学报, 2020, 11(2): 138-144. doi: 10.15886/j.cnki.rdswxb.2020.02.003
    [20] 王子贺, 兰坚强, 廖承红, 韩谦.  天维菌素处理埃及伊蚊幼虫的转录组学分析 . 热带生物学报, 2020, 11(3): 281-287. doi: 10.15886/j.cnki.rdswxb.2020.03.004
  • 加载中
图(5) / 表 (2)
计量
  • 文章访问数:  508
  • HTML全文浏览量:  142
  • PDF下载量:  34
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-01-07
  • 修回日期:  2021-05-19
  • 网络出版日期:  2021-10-08
  • 刊出日期:  2021-10-11

伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建

doi: 10.15886/j.cnki.rdswxb.2021.03.011
    基金项目:  国家自然科学基金(31960703,31860702)
    作者简介:

    李妙珍(1995−),女. 海南大学生命科学与药学院2018级硕士研究生. E-mail:limiaozhen4956@163.com

    通讯作者: 廖承红(1977−),女,教授. 研究方向:病媒生物学. E-mail:liaochh@hainanu.edu.cn韩谦(1963−),男,教授. 研究方向:病媒生物学. Email:qianhan@hainanu.edu.cn
  • 中图分类号: R 384.1

摘要: 5-羟色胺(5-HT)可通过多个受体介导来调控昆虫的生长发育和运动等生理活动,也是新型杀虫剂潜在的靶点。本研究对埃及伊蚊的5-羟色胺受体家族进行系统分析,即利用生物信息学方法、进化分析5-HT受体在物种间的同源性,构建系统进化树;使用同源建模的方法以人类5-羟色胺受体结构为模板,构建埃及伊蚊(Aedes aegypti)5-HT受体的三维结构模型;通过实时荧光定量PCR方法研究受体在埃及伊蚊幼虫期、蛹期和成蚊,未吸血雌蚊的头、胸、腹及附肢的表达情况。从NCBI和VectorBase数据库中得到埃及伊蚊的6个5-HT受体基因,生物信息学分析发现6个受体结构均具有7次跨膜区,其中在受体蛋白的二级结构中α-螺旋占有较高的比重。进化分析结果表明,埃及伊蚊5-HT受体基因分别位于3个不同的分支,各分支均与黑腹果蝇5-HT相应受体具有较高的同源性。表达谱构建结果表明,相比各龄幼虫和雌蚊,5-HT受体均在雄蚊中的表达量最高,且在雌蚊各组织中呈现明显的差异。5-HTR1B、5-HTR2B、5-HTR7B在头、胸、腹和附肢4种组织中均有表达;5-HTR1A在头、胸、腹中均有表达;5-HTR1D在胸腔中高表达,但在头部和附肢中不表达;5-HTR7A在附肢中表达丰度最高,而在头部无表达;该结果预示不同受体介导5-HT的功能可能不同。本研究明确了埃及伊蚊5-HT受体家族蛋白的序列特征、进化关系及组织表达特征,为埃及伊蚊5-HT受体功能的进一步研究奠定了基础。

English Abstract

李妙珍,李奕勋,陈静,等. 伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建[J]. 热带生物学报,2021, 12(3):347−355. DOI:10.15886/j.cnki.rdswxb.2021.03.011 doi:  10.15886/j.cnki.rdswxb.2021.03.011
引用本文: 李妙珍,李奕勋,陈静,等. 伊蚊5-HT受体家族生物信息学分析和时空表达谱的构建[J]. 热带生物学报,2021, 12(3):347−355. DOI:10.15886/j.cnki.rdswxb.2021.03.011 doi:  10.15886/j.cnki.rdswxb.2021.03.011
LI Miaozhen, LI Yixun, CHEN Jing, ZHANG Lei, LIAO Chenghong, HAN Qian. Bioinformatics Analysis and Construction of Spatiotemporal Expression Profile of the Serotonin Receptor Protein Family in Aedes aegypti[J]. Journal of Tropical Biology, 2021, 12(3): 347-355. doi: 10.15886/j.cnki.rdswxb.2021.03.011
Citation: LI Miaozhen, LI Yixun, CHEN Jing, ZHANG Lei, LIAO Chenghong, HAN Qian. Bioinformatics Analysis and Construction of Spatiotemporal Expression Profile of the Serotonin Receptor Protein Family in Aedes aegypti[J]. Journal of Tropical Biology, 2021, 12(3): 347-355. doi: 10.15886/j.cnki.rdswxb.2021.03.011
  • 5-羟色胺(5-hydroxytryptamine,5-HT)又称血清素,是广泛分布于脊椎动物和无脊椎动物的中枢神经系统和外周神经组织中的一种重要的神经递质[1]。5-HT通过与5-HT受体(5-hydroxytryptamine receptor,5-HTR)结合调控许多重要的生理和行为过程,如运动[2]、繁殖[3]、性活动、进食[4]、情绪和感情[5]、睡眠[6]、食欲等[7]。5-HT受体分型复杂,至少有7大类14种亚型在人类基因组中被发现,除5-HT3受体是配体门控离子通道外,其余5-HT受体都是G蛋白偶联受体。根据保守的氨基酸序列和第二信使激活系统,昆虫的5-HT受体可分为5-HTR1、5-HTR2和5-HTR73种[1]

    5-HT和5-HTR在昆虫的免疫系统中发挥关键的调节作用,研究发现菜青虫中的(Pieris rapae5-HTR1B5-HTR2B能介导5-HT引起的血细胞吞噬作用[8],该研究在分子水平上证实了昆虫也与哺乳动物一样,可以通过5-HT来调节神经系统与免疫系统的相互作用,为利用5-HT类药物特异性干扰昆虫免疫系统,增强生物杀虫剂效力提供了可能性。研究还发现通过激活神经元释放出的大量的5-HT,会抑制黑腹果蝇(Drosophila melanogaster)的进食和交配行为[4]。同样,在蜜蜂(Apis mellifera)大脑中注射5-HT,其进食行为受到抑制,但是在肠道中注射5-HT则会刺激肠道肌肉的收缩[9]。在埃及伊蚊(Aedes aegypti)中,使用CRISPR-Cas9技术敲除5-HTR2B会导致埃及伊蚊的生长抑制,并且减少成蚊的脂质积累[10]。目前,研究人员还对5-HTR1A的激动剂衍生物进行设计、合成和生测,发现其能够影响昆虫的生长和幼虫的活性[11],此结果更加表明了5-HT受体可作为有害昆虫的药物靶标。随着全球杀虫剂的滥用,传统小分子杀虫剂的耐药性成为控制节肢动物传播传染病的一个难题,需要具有新的作用方式和对抗性害虫种群活性的杀虫剂化学制剂作为替代。破坏节肢动物G蛋白偶联受体(GPCR)靶标的小分子杀虫剂逐渐成为人们关注的重点,将用于开发下一代杀虫剂[12-13]

    埃及伊蚊(Aedes aegypti)是基孔肯尼雅热、黄热病和登革热等急性传染病的主要传播媒介,每年会造成数亿人感染病毒[14],严重威胁人类健康。已有研究报道,埃及伊蚊吸血行为与5-HT的浓度有关[10]。目前,有关埃及伊蚊5-HT受体家族的功能与结构尚未分析透彻。本研究首先对文献已报道的昆虫5-HT受体蛋白进行生物信息学分析,然后采用多序列比对和系统进化树构建法对5-HT受体基因进行进化分析,最后对每个受体进行不同发育时期的表达谱分析,为埃及伊蚊5-HT受体蛋白的功能研究提供基础,进而为新型杀虫剂寻找靶标。

    • 埃及伊蚊(Rockfeller株)来自军事科学院军事医学研究院微生物流行病研究所。

    • 通过比对NCBI和VectorBase 2个数据库序列,获取埃及伊蚊5-HT受体序列,并下载5-HT受体蛋白的氨基酸序列。利用TMHMM Server v2.0(http://www.cbs.dtu.dk/services/ TMHMM-2.0/)程序在线对蛋白序列进行跨膜结构预测;利用在线工具SOPMA(https://npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl?page=npsa_sopma.html)对埃及伊蚊5-HT受体蛋白的二级结构进行预测;使用Signal P 5.0在线软件预测埃及伊蚊5-HT受体的信号肽;使用在线软件Clustal Omega对不同昆虫的5-HT受体氨基酸序列进行多序列比对,用邻接法(Neighbor-Joining)对氨基酸序列构建进化树。Bootstrap设为1000,其余参数采用默认设置,分析不同昆虫的5-HT受体的进化关系。

      比较Protein Data Base数据库中已知5-HT受体蛋白的晶体结构,采用GeneDoc计算数据库序列相似性分析,发现埃及伊蚊5-HTR7B与已知人类5-HTR1B的序列相似性最高,为38.61%。利用SWISS-MODEL(https://swissmodel.expasy.org/interactive),以Hsap5-HTR1B蛋白的晶体结构[15](PDB:V)为模板,构建埃及伊蚊5-HTR7B三级同源模型。

    • (1)埃及伊蚊总RNA的提取。分别收集50只1龄幼虫,50只2龄幼虫,30只3龄幼虫,10只4龄幼虫,4只白色雌蛹,4只白色雄蛹,4只黑色雌蛹、4只黑色雄蛹,2只雄蚊,2只未吸血的雌蚊,实验重复3次;解剖20只未吸血雌蚊,分别收集其头、胸、腹、附肢4个部位,实验重复3次。采用Trizol法[16]提取埃及伊蚊不同样品的总RNA,用1%的琼脂糖凝胶电泳检测RNA的完整性,用微量核酸检测仪检测RNA浓度和纯度。(2)cDNA模板的合成。使用TAKARA逆转录试剂盒合成cDNA:gDNA Eraser Buffer 2.0 μL, gDNA Eraser 1.0 μL,Total RNA 1.0 μL,RNase Free dH2O 6.0 μL,混匀后,42 ℃反应2 min,转移至冰上,再加入PrimeScript RT Enzyme Mix 10.0 μL,RT Primer Mix 1.0 μL,5×PrimeScript Buffer 2 (for Real Time) 4.0 μL,RNase Free dH2O 4.0 μL,混匀,37 ℃反应15 min,85 ℃反应5 s,−20 ℃保存备用。(3)荧光定量PCR检测不同发育时期和不同组织中的表达量。采用荧光定量试剂盒(ROCHE)对各个基因进行qRT-PCR分析。反应体系为:实时荧光定量PCR预混液5 μL,上、下游引物各0.5 μL(10 μL·mL−1),模板cDNA 1 μL,ddH2O 3 μL。使用Roche公司的LightCycler96进行实时扩增检测。反应采用两步法PCR扩增标准程序:95 ℃ 20 s,95 ℃ 5 s,58 ℃ 30 s,共40个循环。以RPS17为内参基因[16],每个样品设置3次重复,实验所用引物见表1

      表 1  荧光定量PCR引物

      Table 1.  Primers for fluorescence quantitative PCR

      基因名称 Gene引物序列 Sequence of primer(5′-3′)
      上游引物 F下游引物 R
      5-HTR1B TACTGGCGACGGTTATAGG CGACAATATCAGGTGGTTGG
      5-HTR1A CCCACGAGCAATAATACCC CCTTGCGTTTGGAGACC
      5-HTR1D GTTCTACAGTAAATCGTCCTACC CCCGTAGCAATACATCAGC
      5-HTR2B CGAATCTAATTGCTCCTCTCC CTGGACACTTCTGCTATTCC
      5-HTR7A ACGAAACCGCTAGAATACG GTGGCAGCGATATACAGG
      5-HTR7B TTCCTAGAGTACCTCCTAACC GAGTAGAGAGCTTGTTGACG
      RPS17 AAGAAGTGGCCATCATTCCA GGTCTCCGGGTCGACTTC
    • 根据定量检测结果,以2−△△Ct[17]计算样品间基因表达量的差异倍数,用Graphpad Prism 6.0软件T检验法进行数据的统计分析,用LSD法进行多重比较检验,显著性检验水平P < 0.05。

    • 根据数据库数据分析,发现埃及伊蚊有5-HTR1A、5-HTR1B、5-HTR1D、5-HTR2B、5-HTR7A、5-HTR7B共3个5-HT受体家族6个亚基,TMHMM数据跨膜螺旋结构分析结果表明,它们均有7个跨膜区,属于典型的G蛋白偶联受体;利用SignalP 5.0在线软件预测埃及伊蚊5-HT受体蛋白的氨基酸序列[18],发现除5-HTR7B有可能存在SP信号肽外,其他亚基均不存在信号肽。表2结果表明,埃及伊蚊5-HT受体蛋白主要由α螺旋、延伸链以及无规则卷曲构成,并且还有少量的β-转角,其中α螺旋的比例非常高,约占二级结构组成原件的28.09%~47.64%,5-HTR7Bα螺旋比例最高,为 47.64%;5-HTR1Dα螺旋比例最低。每种受体蛋白都含有无规则卷曲,其中5-HTR1A的无规则卷曲比例为49.14%,是6种受体中比例最高的;5-HTR7B的无规则卷曲比例最少。在埃及伊蚊5-HT受体家族的二级结构中,β-转角的含量相对较少,其中5-HTR1Bβ-转角比例最少,5-HTR1Dβ-转角比例最高,为6.71%。

      表 2  埃及伊蚊5-HT受体蛋白二级结构的主要元件比例

      Table 2.  The proportion of major secondary protein structure of 5-HT receptor protein in A. aegypti.

      受体名称
      Receptor
      α-螺旋/%
      α-helix (%)
      延伸链/%
      Extended
      strand (%)
      β-转角/%
      Beta
      turn (%)
      无规卷曲/%
      Random
      coil (%)
      Aae5-HTR1B40.3812.741.9244.95
      Aae5-HTR1A32.0715.523.2849.14
      Aae5-HTR2B37.4119.174.9636.46
      Aae5-HTR1D28.09621.736.7143.46
      Aae5-HTR7A47.6713.562.9735.81
      Aae5-HTR7B47.6417.673.2331.47

      用Proteindatabase中已知Hsap5-HTR1B图1-A)的晶体结构做为模板,根据SWISS-MODEL同源建模软件构建埃及伊蚊5-HTR7B图1-B)的三级结构。其GMQE值为0.4,QMEAN值为−3.93,均在合理范围之内[19-20]。拉氏图是用来描述蛋白质结构中氨基酸残基二面角ψφ是否在合理区域的一种可视化方法,可以反映出蛋白质构象是否合理[21]。结果表明,埃及伊蚊5-HTR7B所有氨基酸都存在在绿色和浅绿色的允许区域,无氨基酸在白色不允许区域(图1-C),证明这些氨基酸残基具有较好的构象空间[22]。结合GMQE值、QMEAN值和拉氏图可知,所建模型合理[23-27]。同源建模结果显示,Aae5-HTR7B具有7个α螺旋,是典型的7次跨膜G蛋白偶联受体。与Hsap5-HTR1B相比,在α4~α7间的loop差别较大。应用Clustal Omega 软件[28]对埃及伊蚊5-HT受体基因家族成员完整ORF进行同源性和保守域比对,结果表明,埃及伊蚊5-HT受体基因家族之间相似性较低,其保守域的同源性约为40%(图2),主要集中在跨膜区。

      图  1  埃及伊蚊5-HT受体蛋白的三级结构建模

      Figure 1.  Tertiary structure modeling of 5-HT receptor protein in A. aegypti

      图  2  埃及伊蚊5-HTR家族保守结构域序列比对

      Figure 2.  Alignment of conserved domains of 5-HT proteins family in A.aegypti

    • 利用MEGA7的邻接法对6个受体基因构建系统进化树(图3)。高粱蚜虫(Melanaphis sacchari)作为半翅目昆虫[29],与埃及伊蚊亲缘关系较远,因此,可以其作为外群构建系统进化树。结果表明,该进化树分为3个分支,Aae5-HTR1D与淡色按蚊的5-HTR1D在同一分支;Aae5-HTR1A与黑腹果蝇5-HTR1A为同一分支,其和致倦库蚊5-HTR2同源性较高;Aae 5-HTR2B从黑腹果蝇5-HTR2B分支上演化而来,并且与Aae5-HTR1D更接近;Aae5-HTR7A5-HTR7B受体与黑腹果蝇5-HTR7聚类为一支。

      图  3  埃及伊蚊5-HT受体蛋白系统进化树(Boostrap值为1000)

      Figure 3.  Phylogenetic tree of 5-HT receptor proteins in A. aegypti

    • 5-HT受体基因在各发育时期均有表达,如图4-A4-E结果所示,5-HTR1A5-HTR1B5-HTR1D5-HTR2B5-HTR7A 5个受体基因在1龄到4龄幼虫中表达水平较低,到蛹期有显著增加(P< 0.05),而5-HTR7B基因在1 ~4龄阶段,表达量与蛹期水平相近,随幼虫发育到成虫,表达量逐渐增高(图4-F)。结果同时表明,显示6个受体基因的表达在雌雄成体中具有明显的性别差异,所有的雄性成体的表达量均高于雌性成体(图5-A5-F),这与家蚕的表达模式不同[30]

      图  4  埃及伊蚊5-HT受体不同发育时期表达谱

      Figure 4.  Expression profile of 5-HT receptor in A. aegypti at different developmental stages

      图  5  埃及伊蚊5-HT 受体不同组织表达谱

      Figure 5.  Expression profile of 5-HT receptor in A. aegypti at different tissues

    • 6个亚基在不同组织中的表达模式差异较大。Aae5-HTR1A 在蚊虫头部的表达量最高,其次是胸腔,在腹部也有较少的表达量,但其在腿部并无表达(图5-A);Aae 5-HTR1B在头部有较高的表达丰度,但是在胸腔、腹部和附肢中表达量较少(图5-B);Aae5-HTR1D在各个组织中的表达量相比于其他受体是最低的,并且它是唯一只在胸腔和腹部表达的5-HT受体,其表达丰度最高的组织是胸腔,在腹部也有较少的表达量(图5-C);Aae5-HTR2B在头、胸、腹和附肢中均存在表达丰度,其中头部的表达量最高,腹部的表达量最少(图5-D);Aae5-HTR7A是唯一在附肢中具有高表达的受体,但其在头部中没有表达丰度,在腹部中存在较少的表达(图5-E);Aae5-HTR7BAae5-HTR2BD 的组织表达丰度相似。

    • 本研究通过生物信息学方法分析了埃及伊蚊5-HT受体的结构,通过同源建模确定5-HT受体具有7个跨膜区,是典型的G蛋白偶联受体,并且跨膜区结构保守。该结果与果蝇、家蚕等昆虫的受体家族结构一致[30-31]

      不同的5-HT受体在昆虫体内分布不同,与5-HT结合在昆虫体内发挥的生理功能也不同[1]。有研究表明,Dm5-HTR1A、Dm5-HTR2ADm5-HTR7在头部表达,主要参与果蝇的短期学习和记忆[32];黑腹果蝇Dm5HT1A受体、Dm5HT1B受体和Dm5HT7受体在幼虫嗅觉选择行为及食欲和厌恶联想嗅觉学习和记忆中具调节性[33]。当黑腹果蝇5-HT受体缺失或表达下调的情况下,果蝇由于血细胞吞噬能力的降低从而在被病原菌感染后更容易死亡。本研究发现Aae5-HTR1AAae5-HTR1B、Aae5-HTR1DAae5-HTR7B在头部特异性高表达,而蚊虫头部具有很多气味性受体,因此,这4个受体可能与蚊虫的吸血、交配、寻找宿主等行为有关,但后续仍需对此进行验证。有研究报道Aae5-HTR2B基因被敲除后埃及伊蚊的生长和繁殖受到影响[10]。本研究还发现Aae5-HTR1D在头部和四肢中并无表达,但是在胸腔中具有显著性高表达,暗示其可能与胸腹部消化功能相关。而Aae5-HTR7A在头部没有表达,是唯一在四肢中特异性高表达的5-HT受体,蚊虫附肢中存在丰富的感受器,能感受外界环境的气味以及环境的温湿度/盐度,因此,推测该受体可能参与外界温湿度或者气味的感知[34]。通过构建埃及伊蚊5-HT受体基因家族时空表达谱还发现,埃及伊蚊5-HT受体在蚊虫的各个时期都有表达,说明5-HT受体在埃及伊蚊的整个生活史中均发挥作用,可能是蚊虫基本生理活动必需的。其中,Aae5-HTR1D、Aae5-HTR1B、Aae5-HTR7A在幼虫阶段的表达较低,但是在蛹期和成蚊期表达量高,暗示这些受体可能主要与蚊虫的幼虫后期发育及成体的生理功能有关。

      目前研究认为,在昆虫体内5-HT作为重要的神经递质,对昆虫的生长发育以及行为活动等方面都具有许多重要作用,而其发挥作用主要通过分布在体内不同位置的受体介导。因此,通过对埃及伊蚊5-HT受体家族的序列分析和表达谱构建,进而进一步研究其功能,将对探索控制蚊虫及传播病原寻找新型药物提供新的思路。

参考文献 (34)

目录

    /

    返回文章
    返回