哈维氏弧菌与溶藻弧菌共有核酸适配体的筛选及鉴定

    刘慧敏 郝聚敏 许净 鄢庆枇 高晗 郑江

    

    

    

    摘?要?共有位点是指不同微生物上结构相同或相似的位点，有些共有位点是同种或同属微生物所共有的，可以作为微生物种属鉴定的标记位点。因此，筛选这些共有位点的共有核酸适配体，不仅有助于对这些微生物进行总量分析，对于相关微生物的种属鉴定也具有重要意义。本研究以哈维氏弧菌（Vibrio harveyi）和溶藻弧菌（Vibrio alginolyticus）为靶目标，筛选这两种弧菌共有位点的共有核酸适配体。经过8轮筛选，适配体富集库的亲和力提高了24.1倍，同源性分析表明，适配体富集库呈明显的收敛性，表现了较好的筛选效率。对2个核酸适配体C14和C22的亲和性和特异性进行了研究，结果表明，这两个核酸适配体对哈维氏弧菌和溶藻弧菌都有较好的亲和性和特异性，但对哈维氏弧菌的亲和力要显著高于溶藻弧菌，而对非弧菌属的嗜水气单胞菌（Aeromonas hydrophila）和迟钝爱德华氏菌（Edwardsiella tarda）的亲和力较弱。C14对哈维氏弧菌和溶藻弧菌的亲和常数（Kd）分别为55.76和82.88 nmol/L，C22对哈维氏弧菌和溶藻弧菌的亲和常数分别为33.97和43.95 nmol/L。二级结构分析显示，两个核酸适配体都以茎环结构为主，其中C22的3′端结构比C14更复杂。荧光显微镜下，与核酸适配体孵育的哈维氏弧菌、溶藻弧菌、嗜水气单胞菌、迟钝爱德华氏菌和大肠杆菌（Escherichia coli）中，只有哈维氏弧菌和溶藻弧菌呈现出明显亮黄色的荧光，其它3种菌则无显著荧光，进一步证明了两个核酸适配体对这两种弧菌有较好的亲和性和特异性。本研究为这两种弧菌的识别检测及水产病原弧菌的种属鉴定奠定了基础。

    关键词?哈维氏弧菌;溶藻弧菌;共有核酸适配体;筛选;荧光显微镜

    1?引 言

    哈维氏弧菌（Vibrio harveryi）和溶藻弧菌（Vibrio alginolyticus）等水产病原弧菌是养殖环境中的重要病原菌，可感染石斑鱼、银鲷、军曹鱼、鲈鱼、对虾、蟹类等多种水产养殖品种，给养殖业造成巨大的损失[1～4]。另外，它们还可通过食物和水传播，进而感染人类，严重威胁公共安全[5，6]。因此，对这些弧菌进行快速的检测和鉴定十分必要。但由于弧菌的種类多、变异快、基因相似度较高，传统的微生物培养法和16SrDNA等检测方法对弧菌的检测鉴定效果并不理想[7，8]。

    核酸适配体是经体外筛选得到的、长度>10 nt的寡核苷酸序列，可以高效、特异地与靶目标结合，在病原微生物的识别检测、食品安全和靶向治疗等领域得到广泛的关注[9～12]。然而，目前筛选出的病原微生物的核酸适配体大多识别的是微生物上的特异性位点。

    微生物的表面存在有多种位点，有些位点是某种微生物所独有的特异性位点，有些位点则是一些微生物上结构相同或相似的位点，即共有位点。共有位点中有些是同种或同属微生物所共有的，具有种属特异性，可以作为微生物种属鉴定的标记位点。对这些共有位点进行识别，筛选其共有核酸适配体，不仅有助于对这些微生物进行总量分析，对于相关微生物的种属鉴定也具有重要意义。对于不同病原微生物中共有位点的识别及其核酸适配体的筛选，目前还未见报道。本研究以溶藻弧菌和哈维氏弧菌为共同的靶目标，筛选能特异性识别这两种弧菌共有位点的核酸适配体，为后续这两种微生物的识别检测及其种属的识别鉴定奠定了基础。

    2?实验部分

    2.1?仪器与试剂

    BIO-TEK全自动酶标仪（美国Synergy HT公司）;MG96G PCR仪（杭州朗基科学仪器有限公司）;BioDoc-It凝胶成像一体机（美国UVP公司）;TGL-16C高速台式离心机（上海安亭科学仪器厂）;JM-250电泳仪和电泳槽（大连竟迈生物科技有限公司）;SW-CJ超净工作台（苏州泰安空气技术有限公司）;HI98127 pH计（意大利Hanna公司）;DM5000B荧光显微镜（德国莱卡公司）。

    随机ssDNA文库： 5′-TCA GTC GCT TCG CCG TCT CCT TC-N35-GCA CAA GAG GGA GAC CCC AGA GGG-3′ （两端为与引物结合的固定序列，中间为35nt的随机序列）;引物P1： 5′-TCA GTC GCT TCG CCG TCT CCT TC-3′，引物P2： 5′-CCC TCT GGG GTC TCC CTC TTG TGC-3′，引物P3： 5′端标记有地高辛的引物P1;5′端标记有FAM荧光基团的核酸适配体C14： FAM-5′-TCA GTC GCT TCG CCG TCT CCT TC GGG GGC GCG GTG AGG GGC TGC ACA AGA GGG AGG CAC AAG AGG GAG GCA CAA GAG GGA GAG CCC AGA GGG-3′;5′端标记有FAM荧光基团的核酸适配体C22： FAM-5′-TCA GTC GCT TCG CCG TCT CCT TC TGC AGG GCC AGA ACA GGG GGA AGG CAC AAG AGG GAG CAC AAG AGG GAG GCA CAA GAG GGA GAG CCC AGA GGG-3′。以上文库、引物及核酸适配体序列均由上海生工生物公司合成和标记。Dreamtaq DNA聚合酶和缓冲液购自Fermentas公司;dNTP购自Generay Biotech公司;DNA marker购自天根生化科技有限公司;琼脂糖购自厦门太阳马生物工程有限公司;辣根过氧化物酶标记的兔抗地高辛抗体购自北京博奥森生物技术公司。实验用水均为超纯水。

    20×结合缓冲液： ?5.844 g NaCl、 3.725 g KCl、6.06 g三羟甲基氨基甲烷（Tris-HCl）和2.033 g MgCl2·6H2O用超纯水溶解并定容至100 mL，调pH至7.4，使用时用超纯水稀释为2×和1×缓冲液。番红染色液： 番红O（safranine）2.5 g，溶于100 mL 95%乙醇中，配成番红原液，于密闭的棕色瓶中避光保存。使用时稀释5倍，配成番红染色液。

    2.2?微生物及培养基

    目标菌： 哈维氏弧菌（V. harveryi）、溶藻弧菌（V. alginolyticus）;非目标菌： 嗜水气单胞菌（Aeromonas hydrophila）、迟钝爱德华氏菌（Edwardsiella tarda）、 大肠杆菌（Escherichia coli）。大肠杆菌采用LB培养基进行培养，其它菌采用胰蛋白胨大豆肉汤培养基（TSB）进行培养。所用微生物由集美大学病原微生物实验室提供。

    2.3?共有核酸适配体的筛选

    2.3.1?ssDNA文库与目标菌哈维氏弧菌结合?筛选流程如图1所示，具体过程如下： 分别取哈维氏弧菌菌液（含菌量为4×108个），6000 r/min离心5 min，弃上清液，1×结合缓冲液洗涤1次，2×结合缓冲液100 μL悬浮菌沉淀;然后取10 μmol/L 随机ssDNA文库30 μL，用2×结合缓冲液稀释至100 μL，95℃变性5 min，冰浴10 min后，加入到上述菌悬液中，30℃，100 r/min摇振结合一定时间： 第1和第2轮的筛选分别结合2 h，第3轮筛选结合1 h，第4轮以后的筛选结合0.5 h。完成结合后的菌悬液于6000 r/min离心5 min，弃上清液，用1×结合缓冲液洗沉淀1～3次（其中，第1和第2轮筛选洗涤1次，第3轮筛选洗涤2次，第4轮以后的筛选洗涤3次），洗去未与目标菌哈维氏弧菌结合的或结合较弱的ssDNA。再在沉淀中加入1×结合缓冲液100 μL重悬，96℃加热5 min，使与目标菌结合的ssDNA变性，从而与目标菌分离，15000 r/min离心10 min，取上清液，得到可与目标菌哈维氏弧菌结合的ssDNA。

    2.3.2?与目标菌溶藻弧菌结合

    将分离到的可与哈维氏弧菌结合的ssDNA（约100 μL），按照2.3.1节的方法与溶藻弧菌菌液结合和分离，从而得到能同时结合哈维氏弧菌和溶藻弧菌的ssDNA。

    2.3.3?不对称PCR扩增?25 μL PCR反应体系： 1 μL 10 μmol/L引物P1，1 μL 0.2 μmol/L引物P2，0.4 μL 10 mmol/L dNTP，2 μL 10×PCR Buffer，1.2 μL 25mmol/L MgCl2，2 μL模板（即经过2.3.1和2.3.2节后分离得到的ssDNA），0.2 μL 5 U/μL DNA聚合酶，17.2 μL ddH2O。PCR循环参数： 94℃预变性4 min，94℃变性30 s，58℃退火30 s，72℃延伸20 s;?40个循环，最后于72℃延伸7 min。

    2.3.4?循环重复筛选?取不对称PCR产物100 μL作为次级文库，依次与目标菌哈维氏弧菌和溶藻弧菌（数量各为4×108个）结合，重复上述2.3.1至2.3.3节中的步骤，共进行8个循环的筛选。

    2.3.5?反筛?为提高核酸适配体的特异性，分别在第4轮用非目标菌嗜水气单胞菌进行了一次反筛，在第6轮用迟钝爱德华氏菌进行了一次反筛，以去除能与嗜水气单胞菌和迟钝爱德华氏菌结合的ssDNA。具体方法如下： 取非目标菌菌液（含菌量为4×108个），按2.3.1节的方法与次级文库100 μL结合0.5 h，然后以6000 r/min离心5 min，取上清液，再按2.3.1到2.3.3节中步骤进行正常筛选。

    2.4?克隆测序及结构分析

    以最后一轮（即第8轮）筛选产物作为模板，进行常规PCR扩增。25 μL PCR体系： 0.5 μL 10 μmol/L引物P1，0.5 μL 10 μmol/L引物P2，0.5 μL 10 mmol/L dNTP，2.5 μL 10×PCR buffer，1.5 μL 25 mmol/L MgCl2，2 μL模板，0.25 μL 5 U/μL DNA聚合酶，17.25 μL ddH2O;热力学循环参数： 94℃预变性4 min，94℃变性30 s，58℃退火30 s，72℃延伸20 s;扩增35个循环;最后72℃延伸7 min。PCR产物由上海美吉生物医药科技有限公司进行克隆测序。对测序获得的一级结构序列用DNAMAN软件进行同源性分析，并构建同源树，二级结构用RNAstructure4.6软件按能量最低原则进行模拟分析。

    2.5?地高辛-抗体-过氧化物酶法测定核酸适配体的亲和力

    参考文献[13]方法，利用地高辛-抗地高辛抗体-过氧化物酶法测定核酸适配体的亲和力。其原理为： 有地高辛标记的核酸适配体与菌结合后形成地高辛标记的复合物，用辣根过氧化物酶标记的抗地高辛的抗体识别和结合菌复合物上的地高辛，形成辣根过氧化物酶-抗體-地高辛-核酸适配体-菌的复合物，最后通过辣根过氧化物酶的显色反应产生的OD值表征与菌结合的适配体的亲和力，显色越深，OD值越高，亲和力越高。其中，筛选过程中筛选产物的亲和力测定是通过PCR将地高辛标记到核酸适配体上，而纯核酸适配体的亲和力测定中，则是直接采用已标记有地高辛的核酸适配体，无需PCR扩增。具体方法如下所述。

    2.5.1?筛选产物的亲和力测定?筛选产物是筛选过程中获得的含有多种适配体的混合物，对其进行亲和力测定可以跟踪反映筛选进行的效果。分别取哈维氏弧菌和溶藻弧菌菌液（各4×108个），1×结合缓冲液洗涤1次，弃上清液;2×结合缓冲液100 μL悬浮菌沉淀后，将两种菌悬液混合作为靶目标菌液。同时用引物P2和P3（即带地高辛标记的P1引物）对各轮筛选产物进行不对称PCR扩增（不对称PCR参数同上，用引物P3代替引物P1），得到各轮筛选产物的不对称PCR产物。取这些带有地高辛标记的不对称PCR产物各50 μL于离心管中作为实验组，同时用2×结合缓冲液50 μL代替不对称PCR产物作为对照组，于95℃变性5 min，然后冰浴10 min后，加到上述100 μL的目标菌悬液中，30℃下，100 r/min摇床结合30 min;?6000 r/min离心5 min，弃上清液，1×结合缓冲液洗涤1次，然后加入100 μL用TBS稀释过的辣根过氧化物酶标记兔抗地高辛抗体（稀释的体积比为抗体∶TBS=1∶1000），充分混合后反应10 min;离心，弃上清液，1×结合缓冲液洗涤3次，加入200 μL TMB显色液避光显色10 min后，用200 μL 2 mol/L H2SO4终止反应，用酶标仪于450 nm测定OD值，获得相应的亲和力： 亲和力= OD450（实验组）-OD450（对照组）。

    2.5.2?核酸适配体的亲和力测定?分别取标有地高辛的核酸适配体各50 μL（0.2 μmol/L）作为实验组，用2×结合缓冲液50 μL代替适配体作为对照组，按2.5.1节的方法与菌液结合、洗涤后进行显色反应，用酶标仪测450 nm的OD值，获得相应适配体的亲和力： 亲和力= OD450（实验组）-OD450（对照组）。

    2.6?亲和常数的测定

    将5′端标记地高辛的核酸适配体稀释到不同浓度，浓度范围为0～140 nmol/L，然后分别与等量的哈维氏弧菌（4×108个）结合，按2.5.2节方法测定相应浓度下此适配体的亲和力。以适配体浓度为横坐标，亲和力为纵坐标，用Origin 8.0软件的Hyperbola函数对数据进行拟合分析，获得相应适配体与哈维氏弧菌的结合曲线、亲和常数Kd值及其拟合系数。按同样的方法可得到适配体与溶藻弧菌的结合曲线、亲和常数和拟合系数。

    2.7?荧光显微镜观察验证适配体的亲和性和特异性

    用2×结合缓冲溶液将浓度为10 μmol/L的5′端标记有荧光基团FAM的核酸适配体稀释3000倍，然后于95℃变性5 min，冰浴10 min。各取100 μL浓度为108 cfu/mL菌液于6个离心管中，在3管中加入100 μL上述处理过的核酸适配体作为实验组，另外3管加入100 μL 2×缓冲溶液作为空白对照组，然后在30℃、100 r/min摇床结合30 min。然后6000 r/min离心5 min，弃上清液，2×结合缓冲溶液洗涤2次后，用100 μL 2×结合缓冲溶液重悬。取重悬后的混合菌液1 μL涂于玻片上，酒精灯加热固定，再用番红染液染色1 min后，用超纯水洗净，晾干后，置于荧光显微镜下观察对照组和实验组在明场（Bright field）和荧光（Fluorescence）下的图片。

    2.8?统计分析

    实验数据采用Excel软件中的t-检验函数进行组间差异分析，p<0.05为有显著差异，p<0.01为有极显著差异。

    3?结果与讨论

    3.1?核酸适配体的筛选

    筛选过程中亲和力随着筛选轮数的增加而上升，从第1轮的0.0150上升到第8轮的0.3765，上升了24.1倍（图2），说明随着筛选进行，结合到目标菌上的核酸适配体逐步增加;最后的第7轮和第8轮的亲和力已无显著差异（p>0.05），说明结合到两个目标菌哈维氏弧菌和溶藻弧菌上的核酸适配体数量已基本饱和。对第8轮的筛选产物进行PCR扩增和测序，并利用DNAMAN软件对测序获得的序列进行同源性分析，获得同源树（图3），同源树节点处的百分数表示该节点左边所涉及的核酸适配体之间的同源性比率。由图3可见，在第8轮的筛选产物中，有相当多的核酸适配体，其同源性已达到90%以上，有些甚至达到了100%，而最初筛选所用的随机文库中，其随机ssDNA的同源性理论上是接近零的，这说明随着筛选进行，筛选文库中ssDNA的同源性在快速增加，适配体序列在一级结构上呈现出明显的收敛性。这也说明，经过8轮筛选，初始随机文库中的ssDNA呈现出明显定向进化趋势，文库中与目标菌有较好亲和力的核酸适配体得到了较好的富集。

    有关病原菌核酸适配体的筛选已有文献报道，但不同的研究针对不同的靶細菌，筛选效率差异较大。Marton等[14]以大肠杆菌为靶细菌，经过12轮筛选，获得的富集库的亲和力比初始文库大约提高了6倍，Yang等[15]以沙门氏菌为靶目标，经过17轮筛选，富集库对靶细菌的亲和力提高了9倍以上，Yu等[16]以溶藻弧菌为靶目标，经过7～9轮筛选，其富集库的亲和力提高了约4倍。本研究则以溶藻弧菌和哈维氏弧菌作为共同靶细菌，经过7～8轮筛选，获得了亲和力提高了24.1倍的富集库，显示出较高的筛选效率，筛选效率的提高可能与反筛和逐步增加的洗涤次数有关。在第4轮和第6轮分别增加了对嗜水气单胞菌（A. hydrophila）和迟钝爱德华氏菌（E. tarda）的反筛，较好去除了一些非特异性结合的ssDNA。同时，在第1和第2轮分别洗涤1次，第3轮洗涤2次，第4～8轮都分别洗涤3次，有效去除了一些亲和力较弱的核酸适配体，从而使最终富集文库中适配体的亲和力得到较大提升。

    3.2?共有核酸适配体的亲和性和特异性

    图3的同源树是根据适配体序列的同源性大小，采用DNAMAN软件通过同源性分析得到的，其节点处的百分数表示所涉及的核酸适配体之间的同源性比率。从图3可见，同源性较高、适配体数量较多的主要有两个群体，一个是同源性在87%以上、由14个核酸适配体组成的群体（C1、C17、C16、C39、C22、C34、C5、C8、C9、C42、C32、C23、C35、C24），另一个是同源性在92%以上、由13个适配体组成的群体（C14、C25、C48、C7、C37、C28、C29、C15、C27、C30、C6、C31、C38），其它适配体由于数量较少或同源性不够高，缺乏代表性。因此，从有代表性的两个群体中各选择出一个核酸适配体，即适配体C14和C22，进行亲和性和特异性考察、亲和常数的测定及其二级结构模拟。图4显示了这两个共有核酸适配体对目标菌的亲和性和特异性，可见两个核酸适配体对目标菌哈维氏弧菌和溶藻弧菌的亲和力都显著高于非目标菌嗜水气单胞菌和迟钝爱德华氏菌（p<0.01），对哈维氏弧菌的亲和力也显著高于溶藻弧菌（p<0.01），说明这两个核酸适配体不仅能区分目标菌和非目标菌，还能区分同为弧菌属的哈维氏弧菌和溶藻弧菌，显示出较好的亲和性和特异性。

    3.3?共有核酸适配体的亲和常数及其二级结构

    两个核酸适配体C14和C22对目标菌哈维氏弧菌和溶藻弧菌的亲和常数Kd及其拟合曲线如图5所示，两个核酸适配体对哈维氏弧菌的Kd都比溶藻弧菌低，说明两个适配体对哈维氏弧菌的亲和力要更高，这与上述亲和性和特异性的研究结果是一致的。

    目前报道的核酸适配体的Kd多在pmol/L～nmol/L水平[17]。水产病原弧菌中，哈维氏弧菌核酸适配体的Kd在10～50 nmol/L之间[18]，副溶血弧菌核酸适配体的Kd在10～30 nmol/L之间[19]，溶藻弧菌适配体的Kd则在10～100 nmol/L之间[13，16，20]，而本研究得到的这两个共有适配体，其对哈维氏弧菌和溶藻弧菌的Kd也都是介于10～100 nmol/L，与上述研究结果相近。

    软件模拟出的两个核酸适配体（含两端固定序列）的二级结构如图6所示，以茎环和口袋结构为主，茎以G-C配对为主。两个适配体都含有6个封闭环，C14的封闭环较小，C22的封闭环较大，而且C22的3′端结构要比C14更复杂些。越复杂的结构，其Kd通常会较小，亲和力通常也较高，这可能是因为较复杂的结构可能有较多的位点能与靶目标结合，从而能够提高其对靶目标的亲和力，这与本研究组前期的研究结果相近[18]。

    3.4?荧光显微镜鉴定

    如图7所示，目标菌哈维氏弧菌、溶藻弧菌和非目标菌嗜水气单胞菌、迟钝爱德华氏菌、大肠杆菌自身均无明显的荧光（如图7的Vh，Va，Ah，Et，Ec），加入FAM荧光标记的适配体C14或C22后，目标菌哈维氏弧菌和溶藻弧菌均呈显明显的黄绿色荧光（如图7的Vh+C14，Vh+C22，Va+C14，Va+C22），但非目标菌嗜水气单胞菌、迟钝爱德华氏菌和大肠杆菌则未呈现出明显的荧光（如图7的Ah+C14，Ah+C22，Et+C14，Et+C22，Ec+C14，Ec+C22），说明这两个核酸适配体都能与哈维氏弧菌、溶藻弧菌特异结合，使这两种目标菌呈现出荧光，但对非目标菌嗜水气单胞菌、迟钝爱德华氏菌和大肠杆菌则基本不结合，进一步证明这两个核酸适配体对哈维氏弧菌和溶藻弧菌有较好的亲和性和特异性。

    4?结 论

    经过8轮筛选，获得了能识别哈维氏弧菌和溶藻弧菌的共有核酸适配体富集文库，富集文库的亲和力提高了24倍以上。研究了两个核酸适配体C14和C22的亲和性和特异性，结果表明，两个核酸适配体对哈维氏弧菌和溶藻弧菌都有较好的亲和性和特异性，但对哈维氏弧菌的亲和力要显著高于溶藻弧菌。Kd的测定结果也表明，两个核酸适配体对哈维氏弧菌的Kd都比溶藻弧菌低，说明两个适配体对哈维氏弧菌的亲和力要更高。最后，用荧光显微镜进一步证明了这两个核酸适配体对哈维氏弧菌和溶藻弧菌都有较好的亲和性和特异性，可用于这两种弧菌的识别鉴定。利用核酸适配体进行微生物的种属鉴定具有快速、准确、便捷的特点，但是由于核酸适配体是通过空间结构来识别微生物的，因此当微生物的形态结构出现较大变化，影响到相应识别位点的结构时，对核酸适配体的识别鉴定效果也会产生一定影响。

    References

    1?Rameshkumar P，Nazar A K A，Pradeep M A，Kalidas C，Jayakumar R ，Tamilmani G，Sakthivel M，Samal A K，Sirajudeen S，Venkatesan V，Nazeera B M. Lett. Appl. Microbiol.，2017，65（5）： 423-430

    2?Austin B，Austin D A. Bacterial Fish Pathogens. Heidelberg： Springer，2012： 365-396

    3?Selvin J，Lipton A P. Dis. Aquat.Org.，2003，57（1-2）： 147-150

    4?Aguirre-Guzmán G，Ruíz H M，Ascencio F. Aquacult. Res.，2004，35（15）： 1395-1404

    5?Schmidt U，Chmel H，Cobbs C. J. Clin. Microbiol.，1979，10（5）： 666-668

    6?Austin B. Vet. Microbiol.，2010，140（3-4）： 310-317

    7?Wei S，Zhao H，Xian Y Y，Hussain M A，Wu X Y. Diagn. Microbiol. Infect. Dis.，2014，79（2）： 115-118

    8?Sawabe T，Kita T K，Thompson F L. J. Bacteriol.，2007，189 （21）： 7932-7936

    9?Soontornworajit B，Wang Y. Anal. Bioanal. Chem.，2011，399（4）： 1591-1599

    10?Farokhzad O C，Cheng J J，Teply B A，Sherifi I，Jon S，Kantoff P W ，Richie J P，Langer R. Proc. Natl. Acad. Sci. USA，2006，103（16）： 6315-6320

    11?Yazdian-Robati R，Arab A，Ramezani M，Abnous K，Taghdisi S M. Int. J. Pharm.，2017，529（1-2）： 44-54

    12?YANG Ge，ZHU Chao，LIU Xiao-Hui，WANG Yong，QU Feng. Chinese J. Anal. Chem.，2018，46（10）： 1595-1603

    杨 歌，朱 超，刘晓慧，王 勇，屈 锋. 分析化学，2018，46（10）： 1595-1603

    13?Tang X M，Zheng J，Yan Q P，Li Z B，Li Y B. Biotechnol. Lett.，2013，35（6）： 909-914

    14?Marton S，Cleto F，Krieger M A，Cardoso J. PLoS One，2016，11（4）： e0153637

    15?Yang M，Peng Z H，Ning Y，Chen Y Z，Zhou Q，Deng L. Sensors，2013，13： 6865-6881

    16?Yu Q，Liu M Z，Su H F，Xiao H H，Wu S T，Qin X L，Li S Q，Mi H Z，Lu Z J，Shi D Q，Li P F. J. Fish Dis.，2019，42（6）： 851-858

    17?Harleen K. BBA-Gen. Subjects，2018，1862： 2323-2329

    18?ZHENG Jiang，HAO Ju-Min，SONG Lin-Sheng，LIU Rui. Progress in Biochemistry and Biophysics，2014，41（7）： 704-711

    郑 江，郝聚敏，宋林生，刘 瑞. 生物化学与生物物理進展，2014，41（7）： 704-711

    19?Duan N，Wu S J，Chen X J，Huang Y K，Wang Z P. J. Agric. Food Chem.，?2012，60（16）： 4034-4038

    20?Zheng J，Tang X M，Wu R X，Yan Q P，Tang H，Luo J W，Niu S F，Qu Y K，Sun L W. LWT-Food Sci. Technol.，2015，64（2）： 1138-1142