留言板

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

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

锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析

叶碧欢 李海波 叶卫邦 陈友吾

叶碧欢, 李海波, 叶卫邦, 陈友吾. 锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析[J]. 浙江农林大学学报, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
引用本文: 叶碧欢, 李海波, 叶卫邦, 陈友吾. 锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析[J]. 浙江农林大学学报, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
YE Bihuan, LI Haibo, YE Weibang, CHEN Youwu. Sequence analysis of COⅠ, COⅡ, Cytb and 28S gene fragments from Apriona swainsoni larvae[J]. Journal of Zhejiang A&F University, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
Citation: YE Bihuan, LI Haibo, YE Weibang, CHEN Youwu. Sequence analysis of COⅠ, COⅡ, Cytb and 28S gene fragments from Apriona swainsoni larvae[J]. Journal of Zhejiang A&F University, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015

锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析

doi: 10.11833/j.issn.2095-0756.2020.02.015
基金项目: 

浙江省科技计划项目 2016F30012

浙江省省属科研院所资助项目 2019F1068-1

详细信息
    作者简介: 叶碧欢, 助理研究员, 博士, 从事森林病虫害综合防治等研究。E-mail:393938188@qq.com
  • 中图分类号: S763.38

Sequence analysis of COⅠ, COⅡ, Cytb and 28S gene fragments from Apriona swainsoni larvae

  • 摘要:   目的  基于分子生物学手段,探讨可快速准确鉴别锈色粒肩天牛Apriona swainsoni幼虫与其他天牛幼虫的目标基因,为该保健昆虫的食药用安全提供保障。  方法  利用线粒体COⅠ、COⅡ、Cytb基因和28S细胞核基因的通用引物分别克隆获得4种基因片段,并进行同源序列检索、多重比对、序列信息分析及进化树构建。  结果  基因克隆最终获得锈色粒肩天牛4种基因序列的片段大小分别为817、545、434和1 088 bp。特殊位点分布结果显示:28S的保守位点占比最高,其后依次为CytbCOⅠ和COⅡ,变异率则正好相反。COⅠ、COⅡ和Cytb基因片段的碱基组成和替换信息特点均表现为A+T含量>G+C含量,颠换高于转换,28S基因片段则表现为A+T含量<G+C含量,颠换低于转换。序列信息分析结果表明:线粒体基因和核基因的序列结构存在差异,核基因较线粒体基因更为保守。4种基因序列的分子进化树均能将锈色粒肩天牛幼虫与桑天牛A.germari等其他天牛幼虫区分开。  结论  COⅠ、COⅡ、Cytb和28S基因均可作为鉴定该虫的分子生物学参考依据。
  • 图  1  锈色粒肩天牛及天牛科其他昆虫基于COⅠ、COⅡ、Cytb和28S基因序列的分子进化关系

    Figure  1  Molecular phylogenetic relationship of A. stoainsoni and other insects of Cerambycidae based on COⅠ, COⅡ, Cytb and 28S gene sequences

    表  1  本研究供试引物信息

    Table  1.   Information of tested primers

    扩增基因片段 引物序列(5'→3') 退火温度Tm/℃ 参考文献
    CO 2195:TTGATTTTTTGGTCATCCAGAAGT 54 [22, 29]
    3014:TCCAATGCACTAATCTGCCATATTA
    CO COSz-F:TGCTTCAAGATAGAGCCTCTCC 54 [21, 28]
    COSz-R:GGTTTGCTCCACAGATTTCAG
    Cytb CytbF:TATGTACTACCATGAGGACAAATAT 58 [25, 30]
    CytbR:ATTACACCTCCTAATTTATTAGGAAT
    28S D2-3665F:AGAGAGAGTTCAAGAGTACGTG 58 [22, 31]
    D5-4749R:GTTACACACTCCTTAGCGGA
    下载: 导出CSV

    表  2  COⅠ、COⅡ、Cytb和28S基因同源序列的特殊位点分布

    Table  2.   Distribution of specific sites of COⅠ, COⅡ, Cytb and 28S gene homologous sequences

    基因名称 保守位点 变异性位点 简约信息位点
    数量/个
    单突变位点数
    量/个
    片段大小/bp
    数量/个 百分比/% 数量/个 百分比/%
    CO 473 63.57 271 36.42 218 53 744
    CO 288 54.86 237 45.14 190 47 525
    Cytb 303 66.45 153 33.55 130 23 456
    28 480 73.28 166 25.34 92 74 655
    下载: 导出CSV

    表  3  COⅠ、COⅡ、Cytb和28S基因同源序列碱基组成特点

    Table  3.   Base composition characteristic of COⅠ, COⅡ, Cytb and 28S homologous sequences

    基因
    名称
    碱基含量/% A+T含量/% G+C含量/%
    T C A G
    CO 38.75 15.37 30.97 14.90 69.72 30.27
    CO 40.90 15.96 32.91 10.23 73.81 26.19
    Cytb 39.70 17.76 31.63 10.91 72.53 27.47
    28 21.30 28.78 17.42 32.50 38.72 61.28
    下载: 导出CSV

    表  4  COⅠ、COⅡ、Cytb和28S基因同源序列的碱基替换模式

    Table  4.   Base substitution mode of COⅠ, COⅡ, Cytb and 28S gene homologous sequences

    基因
    名称
    位点 转换 颠换 R 碱基替换频次
    TT TC TA TG CC CA CG AA AG GG
    CO 第1 27 46 0.57 80 23 39 1 5 6 0 89 4 0
    第2 14 6 2.46 64 10 3 1 34 1 0 62 5 67
    第3 3 3 1.12 104 2 1 0 54 0 1 48 1 36
    总计 44 55 0.80 248 35 43 2 93 7 1 199 10 103
    CO 第1 14 10 1.36 46 8 5 2 27 3 0 47 6 32
    第2 3 3 1.07 72 3 2 0 36 1 0 45 0 16
    第3 20 34 0.58 64 18 30 0 3 4 0 55 2 0
    总计 37 47 0.79 182 29 37 2 66 8 0 147 8 48
    Cytb 第1 18 27 0.68 50 17 22 0 7 5 0 48 2 1
    第2 10 7 1.32 40 7 5 1 23 1 1 42 3 29
    第3 2 1 2.13 63 2 1 0 34 0 0 35 0 16
    总计 30 35 0.86 153 26 28 1 64 6 1 125 5 46
    28 第1 7 6 1.24 37 4 2 2 48 1 2 32 3 73
    第2 12 8 1.58 44 10 3 2 53 1 2 36 3 51
    第3 11 6 1.92 30 8 1 1 59 2 1 29 3 67
    总计 30 20 1.50 111 22 6 5 160 4 5 97 9 191
    下载: 导出CSV
  • [1] 王艳平, 温俊宝, 秦中云.林业危险性有害生物锈色粒肩天牛[J].植物检疫, 2006, 20(3):168-170.

    WANG Yanping, WEN Junbao, QIN Zhongyun. The forestry harmful pest-Apriona swainsoni[J]. Plant Quarant, 2006, 20(3):168-170.
    [2] 王晓红, 杨忠岐, 唐艳龙, 等.锈色粒肩天牛及其防治技术研究进展[J].中国森林病虫, 2010, 29(4):24-27.

    WANG Xiaohong, YANG Zhongqi, TANG Yanlong, et al. Advances in research on the damage and control techniques of Apriona swainsoni[J]. For Pest Dis, 2010, 29(4):24-27.
    [3] 刘桂华, 唐燕平.锈色粒肩天牛与寄主树种关系的初步研究[J].林业科学, 2002, 38(3):106-113.

    LIU Guihua, TANG Yanping. The relationships between Apriona swainsoni and its host trees[J]. Sci Silv Sin, 2002, 38(3):106-113.
    [4] 卢希平, 刘玉, 尹淑艳, 等.锈色粒肩天牛生物学观察[J].中国森林病虫, 2005, 24(4):11-13.

    LU Xiping, LIU Yu, YIN Shuyan, et al. Obsetrvation on the bionomics of Apriona swainsoni[J]. For Pest Dis, 2005, 24(4):11-13.
    [5] 刘海军, 骆有庆, 温俊宝, 等.北京地区红脂大小蠹、美国白蛾和锈色粒肩天牛风险评价[J].北京林业大学学报, 2005, 27(2):81-87.

    LIU Haijun, LUO Youqing, WEN Junbao, et al. Pest risk assessment of Dendroctonus valens, Hyphantira cunea and Apriona swainsoni in Beijing area[J]. J Beijing For Univ, 2005, 27(2):81-87.
    [6] 孙丹萍, 徐睿, 梁铁旺, 等.锈色粒肩天牛生态学特性研究[J].西北林学院学报, 2016, 22(5):114-116.

    SUN Danping, XU Rui, LIANG Tiewang, et al. A study on ecological characteristics of Apriona swainsoni population[J]. J Northwest For Univ, 2016, 22(5):114-116.
    [7] 田志超, 徐华潮, 洪沛, 等.斗米虫人工培养的初步研究[J].浙江农业科学, 2012, 53(8):1159-1160.

    TIAN Zhichao, XU Huachao, HONG Pei, et al. Primary research on artificial cultivation of Apriona swainsoni[J]. J Zhejiang Agric Sci, 2012, 53(8):1159-1160.
    [8] 张胜碧.斗米虫对小鼠毒性反应的实验研究[J].海峡药学, 2013, 25(2):29-31.

    ZHANG Shengbi. Experimental study of the toxic effects of Apriona swainsoni on mice[J]. Strait Pharma J, 2013, 25(2):29-31.
    [9] 王文辉, 严员英, 刘晓华, 等.人工饲料对锈色粒肩天牛生长发育的影响研究[J].中国森林病虫, 2017, 36(6):26-28.

    WANG Wenhui, YAN Yuanying, LIU Xiaohua, et al. Effects of artificial diets on growth and development of Apriona swainsoni[J]. For Pest Dis, 2017, 36(6):26-28.
    [10] 叶苓, 林美熙, 黄泰康, 等.虫类本草[M].北京:中国医药科技出版社, 2002:162-163.
    [11] 胡晓峰.虫蛇药用巧治百病[M].沈阳:辽宁科学技术出版社, 2000:10-13.
    [12] 张凯, 嵇保中, 刘曙雯, 等.南京地区天牛幼虫分类研究[J].南京林业大学学报(自然科学版), 2012, 36(5):102-108.

    ZHANG Kai, JI Baozhong, LIU Shuwen, et al. Taxonomic study on the larvae of longicorn beetles in Nanjing[J]. J Nanjing For Univ Nat Sci Ed, 2012, 36(5):102-108.
    [13] HEBERT P D N, CYWINSKA A, BALL S L, et al. Biological identifications through DNA barcodes[J]. Proc Royal Soc London Ser B Biol Sci, 2003, 270:313-321.
    [14] 江亚杰, 汪中明, 张涛, 等.基于DNA条形码技术的储粮害虫碎片鉴定研究[J].中国粮油学报, 2017, 32(8):131-135.

    JIANG Yajie, WANG Zhongming, ZHANG Tao, et al. Molecular identification of stored grain pest fragments by DNA barcoding[J]. J Chin Cereals Oils Assoc, 2017, 32(8):131-135.
    [15] SCHWARZBACH A E, RICKLEFS R E. The use of molecular data in mangrove plant research[J]. Wetlands Ecol Manage, 2001, 9(3):205-211.
    [16] SCHINDEL D E, MILLER S E. DNA barcoding a useful tool for taxonomists[J]. Nature, 2005, 435(7038):17.
    [17] FOLMER O, BLACK M, HOEH W, et al. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit Ⅰ from diverse metazoan invertebrates[J]. Mol Mar Biol Biotechnol, 1994, 3:294-299.
    [18] SMITH M A, POYARKOV N A, HEBERT P D N. COI DNA barcoding amphibians:take the chance, meet the challenge[J]. Mol Ecol Resour, 2008, 8:235-236.
    [19] KIM K G, HONG M Y, KIM M J, et al. Complete mitochondrial genome sequence of the yellow-spotted long-horned beetle Psacothea hilaris (Coleoptera:Cerambycidae) and phylogenetic analysis among Coleopteran insects[J]. Mol Cells, 2009, 27(4):429-441.
    [20] 李枷霖.皮蠹科DNA条码分类技术研究[D].苏州: 苏州大学, 2015.

    LI Jialin. Classificatoin Technique Research of DNA Barcode in Dermestidae[D]. Suzhou: Soochow University, 2015.
    [21] CORRÊA A S, VINSON C C, BRAGA L S, et al. Ancient origin and recent range expansion of the maize weevil Sitophilus zeamais, and its genealogical relationship to the rice weevil S. oryzae[J]. Bull Entomol Res, 2017, 107(1):9-20.
    [22] ANGELINI D R, JOCKUSCH E L. Relationships among pest flour beetles of the genus Tribolium (Tenebrionidae) inferred from multiple molecular markers[J]. Mol Phylogenet Evol, 2008, 46(1):127-141.
    [23] KIM C G, ZHOU H Z, IMURA Y, et al. Pattern of morphological diversification in the Leptocarabus ground beetles (Coleoptera:Carabidae) as deduced from mitochondrial ND5 gene and nuclear 28S rDNA sequences[J]. Mol Biol Evol, 2000, 17(1):137-145.
    [24] SEVCIK J, DAVID K, RULIK B. A new species of Docosia Winnertz from Central Europe, with DNA barcoding based on four gene markers (Diptera, Mycetophilidae)[J]. Zookeys, 2016, 549(549):127-143.
    [25] 王琦琦, 陈丹丹, 唐光辉, 等.核桃举肢蛾线粒体COⅠ和Cytb基因片段序列分析[J].西北林学院学报, 2016, 31(1):145-152.

    WANG Qiqi, CHEN Dandan, TANG Guanghui, et al. Sequence analysis of the fragments of COⅠ and Cytb genes from Atrijuglans hetaohei[J]. J Northwest For Univ, 2016, 31(1):145-152.
    [26] BARKER S C. Distinguishing species and populations of rhipicephaline ticks with ITS2 ribosomal RNA[J]. J Parasitol, 1998, 84(5):887-892.
    [27] JI Yajie, ZHANG Dexing, HE Lijun. Evolutionary conservation and versatility of a new set of primers for amplifying the ribosomal internal transcribed spacer regions in insects and other invertebrates[J]. Mol Ecol Resour, 2003, 3(4):581-585.
    [28] 陈光辉, 李焱, 刘东, 等.玉米象Sitophilus zeamais DNA条形码研究[J].应用昆虫学报, 2018, 55(3):454-463.

    CHEN Guanghui, LI Yan, LIU Dong, et al. DNA barcoding of Sitophilus zeamais (Motschulsky)[J]. Chin J Appl Entomol, 2018, 55(3):454-463.
    [29] SIMON C, FRATI F, BECKENBACH A, et al. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers[J]. Ann Entomol Soc Am, 1994, 87(6):651-701.
    [30] 毛增辉, 郝家胜, 朱国萍, 等.菜粉蝶线粒体基因组的全序列测定和分析[J].昆虫学报, 2010, 53(11):1295-1304.

    MAO Zenghui, HAO Jiasheng, ZHU Guoping, et al. Sequencing and analysis of the complete mitochondrial genome of Pieris rapae Linnaeus (Lepidoptera:Pieridae)[J]. Acta Entomol Sin, 2010, 53(11):1295-1304.
    [31] HANCOCK J M, TAUTZ D, DOVER G A. Evolution of the secondary structures and compensatory mutations of the ribosomal RNAs of Drosophila melanogaster[J]. Mol Biol Evol, 1988, 5:393-414.
    [32] SCHIMITZ J, MORITZR F A. Sequence analysis of the D1 and D2 regions of 28S rDNA in the hornet (Vespa crabro) (Hymenoptera, Vespinae)[J]. Insect Mol Biol, 1994, 3(4):273-277.
    [33] SOTA T, BOCAK L, HAYASHI M. Molecular phylogeny and historical biogeography of the holarctic wetland leaf beetle of the genus Plateumaris[J]. Mol Phylogenet Evol, 2008, 46(1):183-192.
    [34] SASAKAWA K, KUBOTA K. Phylogeny of ground beetles subgenus Nialoe (s. lat.) Tanaka (Coleoptera:Carabidae; genus Pterostichus):a molecular phylogenetic approach:original article[J]. Entomol Sci, 2009, 12(3):308-313.
    [35] YUJIE L, YARU Z, ZHENG Y W, et al. DNA barcoding, species-specific PCR for the identification of three stored-product pest species of genus Palorus (Coleoptera:Tenebrionidae)[J]. J Stored Prod Res, 2018, 78:32-38.
    [36] BUENAVENTUR E, CÉSAR V C, WOLFF M, et al. DNA barcoding for identifying synanthropic flesh flies (Diptera, Sarcophagidae) of Colombia[J]. Acta Trop, 2018, 182:291-297.
    [37] ASHFAQ M, SABIR J S M, El-ANSARY H O, et al. Insect diversity in the Saharo-Arabian region:revealing a little-studied fauna by DNA barcoding[J]. PLoS One, 2018, 13(7):e0199965. doi:10.1371/journal.pone. 019996.
    [38] 杨晋英.五倍子蚜一新种的形态学及分子鉴定[D].太原: 山西大学, 2012.

    YANG Jinying. Morphological and Molecular Identification of a New Species of Rhus Gall Aphids[D]. Taiyuan: Shanxi University, 2012.
    [39] CHEN Rui, JIANG Liyun, QIAO Gexia. The effectiveness of three regions in mitochondrial genome for aphid DNA barcoding:a case in Lachininae[J]. PLoS One, 2012, 7(10):e46190. doi:10.1371/journal.pone. 0046190.
    [40] 张健.吉林省天牛科昆虫分类学研究[D].长春: 东北师范大学, 2011.

    ZHANG Jian. Taxonomic Study on the Cerambycidae in Jilin Province[D]. Changchun: Northeast Normal University, 2011.
    [41] 魏子涵.天牛核糖体基因18S rDNA和28S rDNA的分子特征研究[D].郑州: 河南农业大学, 2014.

    WEI Zihan. Moleculat Characteristics Reaserch of 18S Ribosomal DNA and 28S Ribosomal DNA of Longicorn Bettles (Coleptera: Cerambycidae)[D]. Zhengzhou: Henan Agricultural University.
    [42] 李成, 魏颖, 袁金钱, 等.旋毛虫28S rRNA基因片段的克隆及其在分类学上的应用[J].东北农业大学学报, 2011, 42(9):80-83.

    LI Cheng, WEI Ying, YUAN Jinqian, et al. Cloning and application of 28S rRNA gene fragment of Trichinella spiralis on taxonmy[J]. J Northeast Agric Univ, 2011, 42(9):80-83.
  • [1] 叶晓明, 钱宇汀, 叶雯, 沈黄莹, 曾燕如, 喻卫武, 戴文圣.  香榧绿藻的生物学特性及物种鉴定 . 浙江农林大学学报, 2019, 36(4): 629-637. doi: 10.11833/j.issn.2095-0756.2019.04.001
    [2] 刘鹏, 刘美娟, 刘恒圣, 金明霞.  不同人工饲料对斗米虫生长发育的影响 . 浙江农林大学学报, 2018, 35(6): 1128-1132. doi: 10.11833/j.issn.2095-0756.2018.06.017
    [3] 赵艺, 徐华潮, 马艳, 史黎央.  虫酰肼和灭幼脲对锈色粒肩天牛氧化酶和解毒酶活性的影响 . 浙江农林大学学报, 2018, 35(1): 174-177. doi: 10.11833/j.issn.2095-0756.2018.01.023
    [4] 卜宇飞, 祁骁杰, 温俊宝, 许志春.  2种天牛幼虫的4类声行为特征 . 浙江农林大学学报, 2017, 34(1): 50-55. doi: 10.11833/j.issn.2095-0756.2017.01.008
    [5] 朱宁, 张冬勇, 沈红霞, 胡琴, 樊建庭.  光肩星天牛产卵偏好性及寄主植物挥发物成分分析 . 浙江农林大学学报, 2017, 34(6): 1059-1064. doi: 10.11833/j.issn.2095-0756.2017.06.013
    [6] 张冬勇, 周书永, 朱宁, 樊建庭.  山核桃天牛危害调查及种类的鉴别方法 . 浙江农林大学学报, 2017, 34(1): 42-49. doi: 10.11833/j.issn.2095-0756.2017.01.007
    [7] 王紫薇, 徐华潮, 张娓娓, 王佩星.  光肩星天牛对寄主的选择及主要寄主挥发物的化学成分分析 . 浙江农林大学学报, 2016, 33(4): 558-563. doi: 10.11833/j.issn.2095-0756.2016.04.002
    [8] 杜和芬, 王佩星, 徐华潮, 张娓娓, 王紫薇.  光肩星天牛对山核桃挥发性组分的触角电位分析 . 浙江农林大学学报, 2016, 33(1): 166-171. doi: 10.11833/j.issn.2095-0756.2016.01.022
    [9] 庞景, 童再康, 黄华宏, 林二培, 刘琼瑶.  杉木纤维素合成酶基因CesA的克隆及表达分析 . 浙江农林大学学报, 2015, 32(1): 40-46. doi: 10.11833/j.issn.2095-0756.2015.01.006
    [10] 黄俊浩, 吴时英, 高磊, 杨露菁.  中国新记录种——香樟齿喙象的鉴别与为害 . 浙江农林大学学报, 2014, 31(5): 764-767. doi: 10.11833/j.issn.2095-0756.2014.05.016
    [11] 沈辰, 裘佳妮, 黄坚钦.  山核桃COP1 E3连接酶的全长克隆及表达分析 . 浙江农林大学学报, 2014, 31(6): 831-837. doi: 10.11833/j.issn.2095-0756.2014.06.002
    [12] 刘攀峰, 乌云塔娜, 杜兰英, 吴敏, 黄海燕, 杜红岩.  杜仲2-甲基-D-赤藓糖醇-2,4-环焦磷酸合酶基因全长cDNA克隆与序列分析 . 浙江农林大学学报, 2014, 31(3): 410-416. doi: 10.11833/j.issn.2095-0756.2014.03.013
    [13] 詹伟君, 罗优波, 张立斌.  松墨天牛优良球孢白僵菌菌株生物学及毒力研究 . 浙江农林大学学报, 2013, 30(1): 83-89. doi: 10.11833/j.issn.2095-0756.2013.01.012
    [14] 杨希宏, 黄有军, 陈芳芳, 黄坚钦.  山核桃FLOWERING LOCUS C同源基因鉴定与表达分析 . 浙江农林大学学报, 2013, 30(1): 1-8. doi: 10.11833/j.issn.2095-0756.2013.01.001
    [15] 阮佳佳, 袁进强, 杨仙玉.  日本蟾蜍皮肤claudin-4 cDNA的克隆与序列分析 . 浙江农林大学学报, 2013, 30(4): 561-566. doi: 10.11833/j.issn.2095-0756.2013.04.016
    [16] 宋敏国, 袁进强, 杨仙玉, 张姝芳, 诸葛慧, 徐跃.  日本蟾蜍皮肤胸腺素α原cDNA的克隆及序列分析 . 浙江农林大学学报, 2013, 30(3): 401-405. doi: 10.11833/j.issn.2095-0756.2013.03.016
    [17] 王国立, 安华明, 秦巧平, 李孟娇, 刘真真, 陈佳莹, 周倩, 张岚岚.  柑橘果实成熟特异基因CsPMEI/InvI的克隆与序列分析 . 浙江农林大学学报, 2013, 30(3): 336-342. doi: 10.11833/j.issn.2095-0756.2013.03.005
    [18] 王策, 秦静静, 甘红豪1, 李红, 罗志斌.  毛果杨全基因组磷酸根转运蛋白家族成员序列分析 . 浙江农林大学学报, 2012, 29(4): 516-526. doi: 10.11833/j.issn.2095-0756.2012.04.006
    [19] 王正加, 黄有军, 夏国华, 郑炳松, 金松恒, 黄坚钦.  山核桃APETALA1同源基因的克隆与序列分析 . 浙江农林大学学报, 2008, 25(4): 427-430.
    [20] 李兰英, 高岚, 温亚利, 黄水灵, .  森林病害的环境影响分析 . 浙江农林大学学报, 2006, 23(5): 491-496.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2020.02.015

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2020/2/303

图(1) / 表(4)
计量
  • 文章访问数:  1784
  • HTML全文浏览量:  512
  • PDF下载量:  31
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-10
  • 修回日期:  2019-08-23
  • 刊出日期:  2020-04-20

锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析

doi: 10.11833/j.issn.2095-0756.2020.02.015
    基金项目:

    浙江省科技计划项目 2016F30012

    浙江省省属科研院所资助项目 2019F1068-1

    作者简介:

    叶碧欢, 助理研究员, 博士, 从事森林病虫害综合防治等研究。E-mail:393938188@qq.com

  • 中图分类号: S763.38

摘要:   目的  基于分子生物学手段,探讨可快速准确鉴别锈色粒肩天牛Apriona swainsoni幼虫与其他天牛幼虫的目标基因,为该保健昆虫的食药用安全提供保障。  方法  利用线粒体COⅠ、COⅡ、Cytb基因和28S细胞核基因的通用引物分别克隆获得4种基因片段,并进行同源序列检索、多重比对、序列信息分析及进化树构建。  结果  基因克隆最终获得锈色粒肩天牛4种基因序列的片段大小分别为817、545、434和1 088 bp。特殊位点分布结果显示:28S的保守位点占比最高,其后依次为CytbCOⅠ和COⅡ,变异率则正好相反。COⅠ、COⅡ和Cytb基因片段的碱基组成和替换信息特点均表现为A+T含量>G+C含量,颠换高于转换,28S基因片段则表现为A+T含量<G+C含量,颠换低于转换。序列信息分析结果表明:线粒体基因和核基因的序列结构存在差异,核基因较线粒体基因更为保守。4种基因序列的分子进化树均能将锈色粒肩天牛幼虫与桑天牛A.germari等其他天牛幼虫区分开。  结论  COⅠ、COⅡ、Cytb和28S基因均可作为鉴定该虫的分子生物学参考依据。

English Abstract

叶碧欢, 李海波, 叶卫邦, 陈友吾. 锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析[J]. 浙江农林大学学报, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
引用本文: 叶碧欢, 李海波, 叶卫邦, 陈友吾. 锈色粒肩天牛幼虫的COⅠ、COⅡ、Cytb和28S基因片段序列分析[J]. 浙江农林大学学报, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
YE Bihuan, LI Haibo, YE Weibang, CHEN Youwu. Sequence analysis of COⅠ, COⅡ, Cytb and 28S gene fragments from Apriona swainsoni larvae[J]. Journal of Zhejiang A&F University, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
Citation: YE Bihuan, LI Haibo, YE Weibang, CHEN Youwu. Sequence analysis of COⅠ, COⅡ, Cytb and 28S gene fragments from Apriona swainsoni larvae[J]. Journal of Zhejiang A&F University, 2020, 37(2): 303-310. doi: 10.11833/j.issn.2095-0756.2020.02.015
  • 锈色粒肩天牛Apriona swainsoni隶属于鞘翅目Coleoptera天牛科Cerambycidae沟颈天牛亚科Lamiinae白条天牛族Batocerini粒肩天牛属Apriona[1-2],可危害国槐Sophora japonica、黄檀Dalbergia hupehana、三叉蕨Tectaria subtriphylla、悬铃木Plantanus acerifoli等多种植物[3-6],而其寄生于云实Caesalpinia decapetala树茎的幼虫,异名云实蛀虫、云实蠹虫、黄牛刺虫、黄寮刺虫,在民间一直都有食服的习惯,因“一斗米换一条虫”而得名“斗米虫”[7-9]。近年来,斗米虫因其良好的药食功效[7-8, 10-11]日益获得认同。斗米虫野外采集困难,人工规模化生产尚未成熟,其市场价格因此逐步走高。受经济利益驱使,目前市场上时常发生其他天牛幼虫伪充斗米虫的现象,导致误食中毒等事件,扰乱药食昆虫市场的秩序。天牛科昆虫分类主要依据其成虫形态特征,幼虫形态描述信息缺乏,幼虫分类存在一定难度[12]。在实际工作中也发现,斗米虫与其他天牛幼虫,例如桑天牛A. germari幼虫,在形态方面具有一定的相似性,仅靠幼虫形态特征很容易鉴别错误。此外,因加工炮制等处理导致幼虫形态被破坏,也提高了鉴别难度。分子生物学手段的应用在很大程度上弥补了物种形态学鉴定的不足[13-16]。线粒体细胞色素氧化酶C亚基Ⅰ(Cytochrome C oxidase subunitⅠ,COⅠ)[17-20]CO[21]、16S核糖体RNA(16S rRNA,16S[22]ND基因(NADH dehydrogenase subunit,ND[23]和细胞色素B(Cytochrome b,Cytb[24-25]等,以及细胞核28S核糖体RNA(28S rRNA,28S[23]和核糖体DNA的内转录间隔区(Internal transcribed spacer,ITS)[26-28]等多种基因已广泛应用于多种昆虫的分类鉴别。目前尚未发现锈色粒肩天牛幼虫线粒体基因或核基因序列鉴定的相关报道,因此,本研究对锈色粒肩天牛幼虫的3种线粒体编码基因(COⅠ、COⅡ和Cytb),以及1种核糖体大亚基编码基因(28S)共4种序列展开了研究。此次也对来自3个地区的桑天牛A. germari幼虫样品进行了测序,以补充桑天牛幼虫上述4种基因的序列信息,提高天牛幼虫进化树拓扑结构的准确性。本研究通过多种基因片段的序列分析,为锈色粒肩天牛幼虫的准确鉴定提供依据,也为该虫的食药用安全提供保障。

    • 用于测序的4份锈色粒肩天牛幼虫样品分别来自于浙江金华(S1)、江西上饶(S2)、浙江开化(S3)、浙江淳安(S4),3份桑天牛幼虫样品分别来自于浙江金华(G1)、江西上饶(G2)、浙江开化(G3),幼虫活体置于-4 ℃冰箱保存备用。

    • TIANamp Genomic DNA Kit,天根生物科技(北京)有限公司;DNA凝胶回收试剂盒(GE0101-50)、感受态细胞(Escherichia coli DH5α)、TS-GelRed核酸染料(10 000×水溶液),北京擎科新业生物技术有限公司;pGEM-T Easy载体,美国Promega公司;Regular Agarose,西班牙Biowest公司;50×TAE Tris-乙酸电泳缓冲液,福州Phygene生物科技有限公司。

    • DK-S26电热恒温水浴锅,上海森信实验仪器有限公司;FRESCO 21微量冷冻离心机、NanoDrop 2000超微量分光光度计,美国Thermo Scientific公司;Life ECO扩增仪,杭州博日科技有限公司;DYY-12C型电泳仪,北京六一生物科技有限公司;Bio-Rad凝胶成像系统,美国伯乐公司。

    • 所采用的COⅠ、COⅡ、Cytb和28S基因片段扩增引物信息详见表 1,供试引物委托北京擎科新业生物技术有限公司合成。

      表 1  本研究供试引物信息

      Table 1.  Information of tested primers

      扩增基因片段 引物序列(5'→3') 退火温度Tm/℃ 参考文献
      CO 2195:TTGATTTTTTGGTCATCCAGAAGT 54 [22, 29]
      3014:TCCAATGCACTAATCTGCCATATTA
      CO COSz-F:TGCTTCAAGATAGAGCCTCTCC 54 [21, 28]
      COSz-R:GGTTTGCTCCACAGATTTCAG
      Cytb CytbF:TATGTACTACCATGAGGACAAATAT 58 [25, 30]
      CytbR:ATTACACCTCCTAATTTATTAGGAAT
      28S D2-3665F:AGAGAGAGTTCAAGAGTACGTG 58 [22, 31]
      D5-4749R:GTTACACACTCCTTAGCGGA
    • 取待提取幼虫,用无菌水清洗数遍,滤纸吸干水分,加液氮彻底磨碎。基因组DNA的提取参照天根TIANamp Genomic DNA Kit试剂盒操作说明书。利用NanoDrop 2000超微量分光光度计对样品基因组DNA进行核酸纯度及浓度的检测,于-20 ℃冰箱储藏备用。

    • PCR扩增反应体系为(总体积20 μL):2×TSINGKE Master Mix 10 μL,10 umol·L-1引物对各1.5 uL,模板DNA 3 uL(20 mg·L-1),ddH2O补足至20 μL。扩增程序为:94 ℃预变性5 min后;94 ℃变性30 s,退火45 s(不同引物对的退火温度Tm详见表 1),72 ℃延伸2 min,共35个循环;最后于72 ℃补平7 min,终止温度为4 ℃。PCR扩增反应在LifeECO基因扩增仪(杭州博日)上进行。取PCR扩增产物3 μL,点样于质量分数为1.5%的琼脂糖凝胶上,TS-GelRed核酸染料染色,1×TAE缓冲液中电泳30~40 min(120 V),扩增图谱摄取采用Bio-Rad凝胶成像分析仪。PCR扩增产物经切胶回收、纯化、pGEM-Teasy载体连接、感受态细胞转化(E. coli DH5α)和蓝白斑筛选后,随机挑选5个阳性克隆由北京擎科新业生物技术有限公司进行双向测序。

    • 测序结果利用Contig Express软件进行拼接,得到目的基因片段在美国国家生物信息中心(NCBI)数据库中进行BLAST比对和相似性检索,以确定所获片段是否为目的基因并获得同源序列。利用MEGA 6.0软件将本研究获得的基因片段与GenBank检索获得的同源序列进行多重比对,两端对齐,进行序列信息分析,包括碱基组成、变异位点、简约信息位点等。同时运用邻接法(Neighbor-joining,NJ)构建进化树,模型为K2P(Kimura 2 parameter),Bootstrap置信值估算重复数设定为1 000次。

    • 表 1所示:引物成功扩增获得2种天牛幼虫的4种基因片段。经测序、拼接并去除上下游引物及两端不稳定序列,锈色粒肩天牛幼虫最终的COⅠ、COⅡ、Cytb和28S基因片段大小分别为817、545、434和1 088 bp,桑天牛幼虫基因片段大小则分别为783、545、435和1 096 bp。将各基因片段在GenBank进行BLAST,比对结果显示所获基因片段与其他天牛等昆虫的目标基因片段的相似性达80%以上,说明本研究成功获得2种天牛幼虫的4种不同基因片段。

    • 将锈色粒肩天牛幼虫的4种基因序列与桑天牛幼虫序列以及GenBank中比对获得的其他天牛同源序列进行多重比较并分析各基因序列信息(表 2~4)。其中,COⅠ基因同源序列来自天牛科7属16种,COⅡ基因同源序列来自天牛科12属15种,Cytb基因同源序列来自天牛科3属4种,28S基因同源序列来自天牛科7属7种。

      表 2  COⅠ、COⅡ、Cytb和28S基因同源序列的特殊位点分布

      Table 2.  Distribution of specific sites of COⅠ, COⅡ, Cytb and 28S gene homologous sequences

      基因名称 保守位点 变异性位点 简约信息位点
      数量/个
      单突变位点数
      量/个
      片段大小/bp
      数量/个 百分比/% 数量/个 百分比/%
      CO 473 63.57 271 36.42 218 53 744
      CO 288 54.86 237 45.14 190 47 525
      Cytb 303 66.45 153 33.55 130 23 456
      28 480 73.28 166 25.34 92 74 655

      表 3  COⅠ、COⅡ、Cytb和28S基因同源序列碱基组成特点

      Table 3.  Base composition characteristic of COⅠ, COⅡ, Cytb and 28S homologous sequences

      基因
      名称
      碱基含量/% A+T含量/% G+C含量/%
      T C A G
      CO 38.75 15.37 30.97 14.90 69.72 30.27
      CO 40.90 15.96 32.91 10.23 73.81 26.19
      Cytb 39.70 17.76 31.63 10.91 72.53 27.47
      28 21.30 28.78 17.42 32.50 38.72 61.28

      表 4  COⅠ、COⅡ、Cytb和28S基因同源序列的碱基替换模式

      Table 4.  Base substitution mode of COⅠ, COⅡ, Cytb and 28S gene homologous sequences

      基因
      名称
      位点 转换 颠换 R 碱基替换频次
      TT TC TA TG CC CA CG AA AG GG
      CO 第1 27 46 0.57 80 23 39 1 5 6 0 89 4 0
      第2 14 6 2.46 64 10 3 1 34 1 0 62 5 67
      第3 3 3 1.12 104 2 1 0 54 0 1 48 1 36
      总计 44 55 0.80 248 35 43 2 93 7 1 199 10 103
      CO 第1 14 10 1.36 46 8 5 2 27 3 0 47 6 32
      第2 3 3 1.07 72 3 2 0 36 1 0 45 0 16
      第3 20 34 0.58 64 18 30 0 3 4 0 55 2 0
      总计 37 47 0.79 182 29 37 2 66 8 0 147 8 48
      Cytb 第1 18 27 0.68 50 17 22 0 7 5 0 48 2 1
      第2 10 7 1.32 40 7 5 1 23 1 1 42 3 29
      第3 2 1 2.13 63 2 1 0 34 0 0 35 0 16
      总计 30 35 0.86 153 26 28 1 64 6 1 125 5 46
      28 第1 7 6 1.24 37 4 2 2 48 1 2 32 3 73
      第2 12 8 1.58 44 10 3 2 53 1 2 36 3 51
      第3 11 6 1.92 30 8 1 1 59 2 1 29 3 67
      总计 30 20 1.50 111 22 6 5 160 4 5 97 9 191

      所有序列多重比较,两端对齐去除冗余位点后,最终得到不同长度的各基因片段(456~744 bp),并对所有样品4种基因序列中含有保守位点、变异性位点、简约信息位点和单突变位点的数量及百分比进行了分析统计。表 2显示:28S的保守位点所占比例最高(73.28%),其后依次为Cytb(66.45%)、COⅠ(63.57%)和COⅡ(54.86%);4种基因的变异率则恰好相反,COⅡ序列的变异率最高(45.14%),其后依次为COⅠ(36.42%)和Cytb(33.55%),28S变异率最低(25.34%)。

      COⅠ、COⅡ、Cytb和28S各基因所比对的同源序列碱基组成特点如表 3所示。其中,COⅠ、COⅡ和Cytb基因的T碱基含量最高,其次是A碱基,A+T含量(69.72%、73.81%、72.53%)均明显高于G+C含量(30.27%、26.19%、27.47%),而28S基因的G碱基含量最高,其次是C碱基,A+T含量(38.72%)则明显低于G+C含量(61.28 %)。

      COⅠ、COⅡ、Cytb和28S各基因所比对的同源序列碱基替换信息见表 4。从表 4可知:COⅠ、COⅡ和Cytb基因序列的颠换均高于转换,转换与颠换比(R)平均值约0.8,颠换以T和A之间为主;而28S基因序列R平均值为1.5,即转换高于颠换,转换以T和C之间为主。

    • 分别基于COⅠ、COⅡ、Cytb和28S基因片段,运用邻接法构建锈色粒肩天牛与其他天牛的分子进化树(图 1A~D)。可知,基于COⅠ和28S基因序列的分子进化树中,同属的昆虫均聚集为一类,其中粒肩天牛属的锈色粒肩天牛或桑天牛同种不同样品亦分别聚成一小分支,这与传统的分类方法相一致(图 1AD)。基于COⅡ和Cytb基因序列的分子进化树中,锈色粒肩天牛和桑天牛各自不同样品也以较高的置信值分别聚成一支,但这2种同属的昆虫并未如COⅠ基因序列的进化树聚为一大类。COⅡ基因进化树中,锈色粒肩天牛与星斑天牛属Psacothea的黄星天牛Psacothea hilaris聚集,桑天牛则与白天牛属Olenecamptus的六星天牛Olenecamptus clarus聚在一起。Cytb基因进化树中,锈色粒肩天牛先与厚花天牛属Pachyta的双斑厚花天牛Pachyta bicuneata聚集,再与桑天牛聚类(图 1B~C)。

      图  1  锈色粒肩天牛及天牛科其他昆虫基于COⅠ、COⅡ、Cytb和28S基因序列的分子进化关系

      Figure 1.  Molecular phylogenetic relationship of A. stoainsoni and other insects of Cerambycidae based on COⅠ, COⅡ, Cytb and 28S gene sequences

    • 斗米虫等昆虫天然保健产品因其特有的功效,逐渐被大众熟知和接受,鉴于天牛幼虫形态学鉴别存在的不足之处,需要依靠分子生物学手段完善以提高天牛幼虫鉴别的准确性。相对于其他分子标记法,DNA测序在种及以上高级分类阶元的分析结果更为准确直接,可作为一种有用的“分类”辅助工具[15]。至今为止,国内外已将多种线粒体基因和核基因广泛运用于昆虫的系统发育、遗传进化和物种分类等研究中[32-37]。线粒体基因严格遵守母系遗传方式,进化速度适中,相对保守的同时又有足够的变异,特别适用于解决较低阶元的系统发育关系[24, 38-39],28S基因则是真核生物染色体上编码核糖体大亚基的基因,在生物体内含量较大,在进化过程中较线粒体基因保守,尤其在高级分类阶元系统发育关系研究中应用广泛[40-41]

      本研究的3种线粒体蛋白编码基因(COⅠ、COⅡ和Cytb)和1种核基因(28S)为昆虫分类的常用基因。锈色粒肩天牛的上述4种基因片段及其同源序列的序列信息分析,结果表明线粒体基因(COⅠ、COⅡ和Cytb)和核基因(28S)的组成结构等存在差异。3种线粒体基因同源序列的变异率比28S基因片段高,这与核基因较线粒体基因更为保守这一观点相符合[40]。此外,3种线粒体基因同源序列表现出A和T碱基偏倚,颠换高于转换,而28S同源序列则相反(G和C碱基偏倚且转换高于颠换)。此次所选择分析的锈色粒肩天牛28S基因及其同源序列的碱基偏倚性与其他研究报道的天牛科28S基因相同,其碱基组成特点也与28S rDNA相符合[40]

      不同物种,其分子分类适用的基因也存在差异,根据不同物种选择合适的基因十分关键。CHEN等[39]以大蚜亚科Lachninae为研究对象,评估了3种线粒体基因(COⅠ、COⅡ和Cytb)作为DNA条形码的有效性,发现COⅡ基因较其他两者更合适于不同大蚜的鉴定,Cytb因其序列多样性高,是蚜虫种群遗传学研究的有效标记。杨晋英[38]研究表明:COⅠ、COⅡ和Cytb基因适用于五倍子蚜分类,其分子聚类关系与形态学分类相一致。李成等[42]研究表明:28S非常适合于旋毛虫Trichinella spiralis的分类鉴定及其分子进化关系的研究。本研究显示:利用COⅠ、COⅡ、Cytb和28S基因的通用性引物均能有效地扩增锈色粒肩天牛和桑天牛幼虫的目的片段。基于这4种基因序列构建的进化树亦都能以较高的置信值将锈色粒肩天牛幼虫单独聚成一支,与其他种的天牛区分开来,表明上述4种基因均可作为快速鉴别锈色粒肩天牛幼虫与桑天牛等其他天牛幼虫的分子参考依据。

      总体而言,COⅠ、COⅡ、Cytb和28S这4种基因的进化树拓扑结构与传统分类结果基本相一致,但COⅡ和Cytb基因进化树其属级别的聚类稍有差别,其原因可能是本研究涉及的COⅡ和Cytb序列仅是该基因的部分序列,其所覆盖的信息还不够充足,可通过增加昆虫的目标基因序列数量或长度,以获得更为理想的系统进化树。

参考文献 (42)

目录

    /

    返回文章
    返回