趙彥花 區(qū)又君 溫久福 李加兒 周慧

摘要：【目的】篩選出可用于黃唇魚（Bahaba flavolabiata）遺傳資源評價及親緣鑒定的SSR分子標(biāo)記，為開展其保護(hù)遺傳學(xué)研究打下基礎(chǔ)?！痉椒ā坷肐llumina HiSeq 4000測序平臺對野生黃唇魚樣本進(jìn)行轉(zhuǎn)錄組測序，采用TGICL對轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行聚類去冗余以獲得Unigene，并通過生物信息學(xué)方法進(jìn)行功能注釋分類、代謝通路和SSR特征分析等?！窘Y(jié)果】黃唇魚混合組織的轉(zhuǎn)錄組測序共獲得13.43 Gb數(shù)據(jù)，經(jīng)組裝并去冗余后獲得65047條Unigenes。采用BLAST對組裝獲得的Unigenes進(jìn)行七大功能數(shù)據(jù)庫（NR、NT、GO、COG、KEGG、Swissprot和Interpro）注釋，結(jié)果共有55583條Unigenes被注釋（占85.45%）。通過MISA對黃唇魚的Unigene進(jìn)行SSR搜索，共檢測到29797個SSR位點(diǎn)分布于19664條Unigenes中。SSR位點(diǎn)分布類型及其特征分析結(jié)果顯示，最主要的重復(fù)類型為二核苷酸重復(fù)基元（11855個，占39.79%），其次是單核苷酸重復(fù)基元（9695個，占32.54%）和三核苷酸重復(fù)基元（6663個，22.36%），而四核苷酸、五核苷酸和六核苷酸重復(fù)基元的數(shù)量相對較少。二核苷酸重復(fù)基元類型中重復(fù)最多的基序是AC/GT（8355條），在三核苷酸重復(fù)基元類型中重復(fù)最多的基序是AGG/CCT（1866條）。從隨機(jī)選取的186對SSR引物中，最終篩選出47對擴(kuò)增穩(wěn)定性好且具有多態(tài)性的SSR引物，以二核苷酸和四核苷酸重復(fù)基元的擴(kuò)增結(jié)果居多?！窘Y(jié)論】通過轉(zhuǎn)錄組測序能有效篩選出具有多樣性的黃唇魚SSR分子標(biāo)記，可為黃唇魚的遺傳多樣性分析、親緣鑒定和遺傳圖譜構(gòu)建等后續(xù)研究提供基礎(chǔ)資料。

關(guān)鍵詞： 黃唇魚;SSR分子標(biāo)記;轉(zhuǎn)錄組測序;功能注釋;多態(tài)性

中圖分類號： S965.299? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼： A 文章編號：2095-1191（2019）09-2078-10

Abstract：【Objective】SSR markers， which could be used to evaluate genetic resources and identify genetic relatives of Bahaba flavolabiata， were selected to lay a foundation for the study of conservation genetics. 【Method】Illumina HiSeq 4000 sequencing platform was used to carry out transcriptome sequencing on wild B. flavolabiata. TGICL was used to cluster and remove redundancy from transcriptome data to obtain Unigene. Functional annotation classification， metabolic pathway and SSR characteristics were analyzed by bioinformatics. 【Result】Transcriptome sequencing of mixed tissues of B. flavolabiata obtained a total of 13.43 Gb data， and 65047 Unigenes were obtained after assembly and redundancy removal. Unigenes generated by BLAST were annotated in seven functional databases（NR， NT， GO， COG， KEGG， Swissprot and Interpro）， and a total of 55583 Unigenes were annotated（accounting for 85.45%）. Through SSR search on Unigenes of B. flavolabiata by MISA， a total of 29797 SSRs were detected distributed in 19664 Unigenes. SSR loci distribution type and its characteristic analysis results showed that， the major repeat types were dinucleotide repeat units （11855， accounting for 39.79%）， followed by single nucleotide repeat units（9695， accounting for 32.54%） and trinucleo-tide repeat units（6663， accounting for 22.36%）， while the number of tetrapotide， pentapotide and hexapotide repeat units was relatively small. The most repeated motif in dinucleotide repeats was AC/GT（8355）， while the most repeated motif in dinucleotide repeats was AGG/CCT（1866）. From the 186 pairs of SSR primers randomly selected， 47 pairs of SSR pri-mers with good amplification stability and polymorphism were finally selected， and the amplification results of dinucleotide and tetrauleotide repeat primitives were the majority. 【Conclusion】SSR molecular markers with diversity can be effectively screened out by transcriptome sequencing， which can provide basic data for the genetic diversity analysis， genetic identification and genetic mapping construction of B. flavolabiata.

Key words： Bahaba flavolabiata; SSR makers; transcriptome sequencing; functional annotation; polymorphism

0 引言

【研究意義】微衛(wèi)星（Microsatellite）廣泛分布于真核生物基因組中，是由1～6個核苷酸單位組成的簡單串聯(lián)重復(fù)序列，因此又稱為簡單序列重復(fù)（Simple sequence repeats，SSR）。SSR廣泛分布于基因組的不同位置，其高度多態(tài)性主要來源于串聯(lián)數(shù)目的不同（熊良偉等，2018）。SSR分子標(biāo)記具有分布數(shù)量豐富、多態(tài)性高、共顯遺傳、擴(kuò)增穩(wěn)定和易檢測等優(yōu)點(diǎn)（王忠華等，2008），在種質(zhì)遺傳多樣性分析、親緣關(guān)系鑒定、遺傳圖譜構(gòu)建及分子標(biāo)記輔助育種等領(lǐng)域已得到廣泛應(yīng)用（趙志英等，2018），且逐漸發(fā)展成為經(jīng)濟(jì)魚類遺傳結(jié)構(gòu)分析、親緣關(guān)系鑒定及遺傳圖譜構(gòu)建等方面的有效工具（孫效文等，2008;謝子強(qiáng)等，2013;何福玲等，2017），但SSR分子標(biāo)記具有種質(zhì)特異性，因此使用前必須提前進(jìn)行篩選開發(fā)?！厩叭搜芯窟M(jìn)展】目前，開發(fā)SSR分子標(biāo)記最常用的方法主要有：①從公共數(shù)據(jù)庫中查找SSR位點(diǎn);②近緣物種間引物轉(zhuǎn)移擴(kuò)增;③從基因組DNA中篩選SSR位點(diǎn)（王玲玲等，2017）。第①種方法經(jīng)濟(jì)方便，但局限于已開發(fā)出SSR位點(diǎn)的物種，且引物數(shù)量不會增加。第②種方法在不同物種間共用SSR引物時，其盲目性較大，具有一定的風(fēng)險性和局限性。對于無法從前兩種方法獲得SSR引物的物種，可從基因組DNA中篩選SSR位點(diǎn)，但需經(jīng)建庫、雜交、分離及測序等步驟，工作量大、效率低。采用傳統(tǒng)磁珠富集法開發(fā)SSR分子標(biāo)記效率低且成本高（Zane et al.，2002），而利用高通量測序技術(shù)開發(fā)SSR分子標(biāo)記能極大提高開發(fā)效率，并有效降低成本。轉(zhuǎn)錄組測序技術(shù)可快速獲得生物體內(nèi)較全面的轉(zhuǎn)錄本信息，尤其對于缺乏基因組信息的魚類，轉(zhuǎn)錄組測序技術(shù)是獲得大量分子標(biāo)記的有效手段（王明等，2015）。近年來，隨著轉(zhuǎn)錄組測序技術(shù)的快速發(fā)展，通過高通量測序技術(shù)測定轉(zhuǎn)錄組cDNA序列揭示特定細(xì)胞或組織中表達(dá)的全部基因已在許多物種上得到廣泛應(yīng)用，并開發(fā)出大量基于轉(zhuǎn)錄組的分子標(biāo)記，如夏威夷石斑魚（Epinephelusquernus）（Rivera et al.，2003）、金槍魚（Thunnus thunnus）（Clark et al.，2004）、紅鰭東方鲀（Fugu rubripes）（Furukawa et al.，2004）、大西洋鮭 （Salmo salar）（King et al.，2005）、大菱鲆（Scophthalmus maximus）（Pardo et al.，2005）、羅氏沼蝦（Macrobrachium rosenbergii）（孫成飛等，2015）、二長棘鯛（Parargyrops edita）（楊兵等，2015）和黃姑魚（Nibea albiflora）（龔詩琦等，2016）等?！颈狙芯壳腥朦c(diǎn)】黃唇魚（Bahaba flavolabiata）是我國特有種，俗稱金錢鮸、金錢鳘、大鷗和白花等，隸屬于硬骨魚綱鱸形目石首魚科黃唇魚屬（區(qū)又君等，2010，2011，2012，2013，2014），一直被視為上等補(bǔ)品，尤其是其鰾（俗稱魚膠）甚為珍貴，具有滋補(bǔ)肝腎之功效。由于生態(tài)環(huán)境惡化和過度捕撈等原因，導(dǎo)致黃唇魚自然資源急劇下降，于1988年被列為國家二級重點(diǎn)保護(hù)野生動物，2006年被世界自然保護(hù)聯(lián)盟（IUCN）紅色名錄列為極度瀕危物種，因而制約了黃唇魚生物學(xué)研究的開展。目前，基于SSR分子標(biāo)記的黃唇魚遺傳多樣性分析鮮見報道?！緮M解決的關(guān)鍵問題】利用Illumina HiSeq 4000測序平臺對野生黃唇魚樣本進(jìn)行轉(zhuǎn)錄組測序，采用TGICL對轉(zhuǎn)錄組數(shù)據(jù)進(jìn)行聚類去冗余以獲得Unigene，并通過生物信息學(xué)方法進(jìn)行功能注釋分類、代謝通路和SSR特征分析等，旨在篩選出可用于黃唇魚遺傳資源評價及親緣鑒定的SSR分子標(biāo)記，為開展其保護(hù)遺傳學(xué)研究打下基礎(chǔ)。

1 材料與方法

1. 1 試驗(yàn)材料

轉(zhuǎn)錄組測序樣本為珠江口海域誤捕的黃唇魚，共1尾。用于多態(tài)性引物篩選的樣本來自同一海域的野生黃唇魚，共6尾，剪取其鰭條后酒精保存?zhèn)溆谩?/p>

1. 2 黃唇魚轉(zhuǎn)錄組測序

委托深圳華大基因科技有限公司進(jìn)行黃唇魚各組織總RNA提取，采用Illumina HiSeq 4000進(jìn)行測序，獲得黃唇魚轉(zhuǎn)錄組的原始序列。原始序列經(jīng)數(shù)據(jù)過濾后得到高質(zhì)量的Clean reads，以Trinity對Clean reads進(jìn)行組裝，然后采用TGICL對轉(zhuǎn)錄本進(jìn)行聚類去冗余以獲得Unigene。利用MISA對Unigene的SSR位點(diǎn)進(jìn)行檢測，單核苷酸重復(fù)大于12次、二核苷酸重復(fù)大于6次、三核苷酸和四核苷酸重復(fù)大于5次、五核苷酸和六核苷酸重復(fù)次數(shù)大于4次的判斷為SSR位點(diǎn)。

1. 3 基因組DNA提取與檢測

取6尾黃唇魚樣本尾鰭各30.0 mg左右，使用海洋動物組織基因組DNA試劑盒[天根生化科技（北京）有限公司]提取基因組DNA。取6.0 μL左右的DNA樣品進(jìn)行1%瓊脂糖凝膠電泳檢測，并用NanoDrop ND-1000紫外分光光度儀檢測其濃度和質(zhì)量;獲得的DNA用雙蒸水稀釋成100 ng/μL工作液，-20 ℃保存?zhèn)溆谩?/p>

1. 4 SSR引物設(shè)計與驗(yàn)證

1. 4. 1 SSR引物設(shè)計 使用Primer Premier 3.0對含SSR位點(diǎn)的Unigenes進(jìn)行引物設(shè)計，共設(shè)計出26949對引物。引物設(shè)計原則：GC含量為40%～60%，引物長度為18～23 bp，退火溫度為55～65 ℃，目的片段長度為80～200 bp。從上述引物中隨機(jī)選取186對引物，委托廣州擎科生物技術(shù)有限公司進(jìn)行引物合成，篩選條件：二核苷酸重復(fù)基元在8次以上，三核苷酸和四核苷酸重復(fù)基元在5次以上，五核苷酸和六核苷酸重復(fù)基元在4次以上。

1. 4. 2 PCR擴(kuò)增 在Eppendorf普通梯度PCR儀中進(jìn)行PCR擴(kuò)增，反應(yīng)體系10.0 μL，包括100 ng/μL基因組DNA 1.0 μL，2×PCR Mix 5.0 μL，正、反向引物各0.4 μL，超純水3.2 μL。擴(kuò)增程序：94 ℃預(yù)變性3 min;94 ℃ 30 s，55～65 ℃ 30 s，72 ℃ 30 s，進(jìn)行35個循環(huán);最后72 ℃延伸10 min，4 ℃保存。最終確定各引物最適退火溫度為58 ℃，PCR產(chǎn)物以3%瓊脂糖凝膠電泳進(jìn)行檢測，篩選出121對可擴(kuò)增出目的片段的SSR引物。

1. 4. 3 SSR引物多態(tài)性檢測 將121對SSR引物擴(kuò)增的PCR產(chǎn)物送至廣州擎科生物技術(shù)有限公司進(jìn)行測序和多態(tài)性檢測。測序采用ABI 3730XL儀器檢測得到測序峰圖，然后以Chromas查看峰圖，找到出現(xiàn)多態(tài)性的引物擴(kuò)增結(jié)果，確定具有多態(tài)性的引物。多態(tài)性測序結(jié)果的判斷方法：查看峰圖找到重復(fù)序列位置，確認(rèn)出現(xiàn)套峰的位置是否位于重復(fù)序列區(qū)間。另外，結(jié)合兩個終止峰（A堿基）位置判斷重復(fù)數(shù)的差異數(shù)量，從而確定引物的擴(kuò)增結(jié)果是否存在多態(tài)性。

2 結(jié)果與分析

2. 1 黃唇魚轉(zhuǎn)錄組測序結(jié)果

黃唇魚混合組織的轉(zhuǎn)錄組測序共獲得13.43 Gb數(shù)據(jù)，經(jīng)組裝并去冗余后獲得65047條Unigenes，其總長度為69382127 bp、平均長度為1066 bp，N50（按Unigene長度從大到小排序后逐個累加至所有Unigenes總長度的50%時，最后一個累加的數(shù)值即為N50，用于衡量組裝的連續(xù)性，數(shù)值越大說明組裝效果越好）為2032 bp，GC含量為46.65%。Unigene長度分布情況如圖1所示，其中，最小片段長度為300 bp，數(shù)量最多，有17704條;片段長度大于3000 bp的Unigenes有5013條。

2. 2 Unigene功能注釋結(jié)果

采用BLAST對組裝獲得的Unigenes進(jìn)行七大功能數(shù)據(jù)庫（NR、NT、GO、COG、KEGG、Swissprot和Interpro）注釋，注釋結(jié)果詳見表1。其中，38086條Unigenes注釋到NR數(shù)據(jù)庫（58.55%），54929條Unigenes注釋到NT數(shù)據(jù)庫（84.45%），32384條Unigenes注釋到Swissprot數(shù)據(jù)庫（49.79%），12572條Unigenes注釋到COG數(shù)據(jù)庫（19.33%），31228條Unigenes注釋到KEGG數(shù)據(jù)庫（48.01%），3099條Unigenes注釋到GO數(shù)據(jù)庫（4.76%），28318條Unigenes注釋到Interpro數(shù)據(jù)庫（43.53%）。NR數(shù)據(jù)庫屬于非冗余蛋白序列數(shù)據(jù)庫，其數(shù)據(jù)來源于GenPept、Swissprot、PIR、PDF、PDB和NCBI RefSeq。根據(jù)NR注釋結(jié)果統(tǒng)計注釋基因的物種分布情況，結(jié)果（圖2）顯示，80.68%黃唇魚Unigenes能匹配到大黃魚（Larimichthys crocea）的相應(yīng)蛋白中，表明在所比對的魚類中，黃唇魚與大黃魚的同源性最高，二者親緣關(guān)系最近;其次為深裂眶鋸雀鯛（Stegastes partitus）（2.98%）、虹鱒（Oncorhynchus mykiss）（1.21%）和尼羅羅非魚（Oreochromis niloticus）（1.18%）;匹配到其他物種上的Unigege比例為13.95%。

2. 3 Unigene功能分類結(jié)果

2. 3. 1 COG分類 COG是直系同源家族蛋白數(shù)據(jù)庫，通過COG數(shù)據(jù)庫比對可進(jìn)行功能注釋、歸類及蛋白進(jìn)化分析。黃唇魚轉(zhuǎn)錄組中被COG功能注釋的Unigenes共12572條，而這些Unigenes又被分成25個類別（圖3）。其中，與功能預(yù)測類相關(guān)的Unigene數(shù)量最多（4680條，占37.23%）;其次是與轉(zhuǎn)錄（2015條，占16.03%），復(fù)制、重組和修復(fù)（1988條，占15.81%），翻譯、核糖體結(jié)構(gòu)和生物合成（1821條，占14.48%），信號傳導(dǎo)機(jī)制（1468條，占11.68%），翻譯后修飾、蛋白折疊和分子伴侶（1296條，10.31%）等基因工程類相關(guān)的Unigene;與細(xì)胞運(yùn)動類序列細(xì)胞周期調(diào)控、細(xì)胞分裂和染色體分區(qū)（1515條，12.05%）及代謝類序列碳水化合物轉(zhuǎn)運(yùn)和代謝（1419條，11.29%）相關(guān)的Unigene數(shù)量也較多;而與核結(jié)構(gòu)相關(guān)的Unigene數(shù)量最少，僅有6條。

2. 3. 2 GO分類 GO是一個國際標(biāo)準(zhǔn)的基因功能分類體系，利用GO數(shù)據(jù)庫進(jìn)行比對分析可將Unigenes分為生物工程、細(xì)胞組分和分子功能三大功能類別。由圖4可看出，生物工程過程可細(xì)分為25類，其中與細(xì)胞進(jìn)程相關(guān)的Unigene數(shù)量最多（1804條），其次是與單細(xì)胞生物進(jìn)程相關(guān)的Unigene（1535條），與細(xì)胞殺傷相關(guān)的Unigene數(shù)量最少（1條）;細(xì)胞組分可細(xì)分為18類，其中與細(xì)胞（1291條）、細(xì)胞部分（1271條）和細(xì)胞膜（1046條）相關(guān)的Unigene數(shù)量較多，而與細(xì)胞外基質(zhì)組分（1條）、病毒顆粒（4條）和病毒顆粒部分（4條）相關(guān)的Unigene數(shù)量最少;分子功能可細(xì)分為14類，其中與結(jié)合相關(guān)的Unigene數(shù)量最多（1628條），其次是與催化活性相關(guān)的Unigene（1020條），與化學(xué)排斥物活性和蛋白標(biāo)簽相關(guān)的Unigene數(shù)量最少，均為1條。

2. 3. 3 KEGG分類 KEGG是處理基因組、生物通路、疾病、藥物和化學(xué)物質(zhì)間聯(lián)系的集成數(shù)據(jù)庫。本研究根據(jù)功能不同，將KEGG數(shù)據(jù)庫劃分為細(xì)胞過程、環(huán)境信息處理、基因信息處理、人類疾病、新陳代謝和生物系統(tǒng)六大通徑，并細(xì)分為44類（圖5）。通過KEGG數(shù)據(jù)庫比對分析，得到注釋的Unigene共31228條。其中，與信號分子和相互作用相關(guān)的Unigene數(shù)量最多（5843條），其次是與癌癥在視圖相關(guān)的Unigene，有3288條，而與生物合成和次生代謝物相關(guān)的Unigene數(shù)量最少，僅19條。

2. 4 SSR特征分析結(jié)果

通過MISA對黃唇魚的Unigene進(jìn)行SSR搜索，共檢測到29797個SSR位點(diǎn)分布于19664條Unigenes中。其中，最主要的重復(fù)類型為二核苷酸重復(fù)基元（11855個，占39.79%），其次是單核苷酸重復(fù)基元（9695個，占32.54%）和三核苷酸重復(fù)基元（6663個，22.36%），而四核苷酸、五核苷酸和六核苷酸重復(fù)基元的數(shù)量相對較少（圖6）。二核苷酸和三核苷酸重復(fù)基元種類較豐富，分別有4和10種。在二核苷酸重復(fù)基元類型中重復(fù)最多的基序是AC/GT（8355條），在三核苷酸重復(fù)基元類型中重復(fù)最多的基序是AGG/CCT（1866條）。在6種核苷酸重復(fù)基元分布中，單核苷酸重復(fù)基元次數(shù)分布在12～132次，二核苷酸重復(fù)基元次數(shù)分布在6～74次，三核苷酸重復(fù)基元次數(shù)分布在5～49次，四核苷酸重復(fù)基元次數(shù)分布在5～33次，五核苷酸重復(fù)基元次數(shù)分布在4～13次，六核苷酸重復(fù)基元次數(shù)分布在4～12次，其中在6種核苷酸重復(fù)基元類型中均可見的重復(fù)次數(shù)為12次，而重復(fù)次數(shù)最多的為6次。

2. 5 SSR引物有效性驗(yàn)證和多態(tài)性檢測結(jié)果

在SSR引物預(yù)擴(kuò)增試驗(yàn)中，以合成的186對SSR引物對6尾黃唇魚基因組DNA進(jìn)行擴(kuò)增，PCR產(chǎn)物以3.0%瓊脂糖凝膠電泳進(jìn)行檢測，結(jié)果篩選出121對能擴(kuò)增出目的條帶的SSR引物，擴(kuò)增成功率達(dá)65.05%。將121對SSR引物的PCR擴(kuò)增產(chǎn)物送至廣州擎科生物技術(shù)有限公司進(jìn)行測序和多態(tài)性檢測，經(jīng)比對分析最終篩選出47對擴(kuò)增穩(wěn)定性好且具有多態(tài)性的SSR引物（表2），多態(tài)性比例為38.84%，其中以二核苷酸和四核苷酸重復(fù)基元的擴(kuò)增結(jié)果居多。

3 討論

黃唇魚是我國特有種，但近年來因環(huán)境污染及過度捕撈等原因，導(dǎo)致其資源量急劇下降，因此，亟需開展黃唇魚保護(hù)遺傳學(xué)的相關(guān)研究。本研究利用高通量測序技術(shù)對黃唇魚轉(zhuǎn)錄組遺傳資源進(jìn)行挖掘，共獲得13.43 Gb數(shù)據(jù)，經(jīng)組裝并去冗余后得到65047條Unigenes，其總長度為69382127 bp，平均長度為1066 bp，N50為2032 bp，GC含量為46.65%，表明黃唇魚轉(zhuǎn)錄組數(shù)據(jù)具有較高質(zhì)量，可為開展黃唇魚保護(hù)遺傳學(xué)研究提供支撐。在黃唇魚基因組信息缺乏的情況下，利用BLAST對組裝獲得的Unigenes進(jìn)行七大功能數(shù)據(jù)庫（NR、NT、GO、COG、KEGG、Swissprot和Interpro）注釋，共有55583條Unigenes被注釋（占85.45%），通過Unigene注釋可為黃唇魚功能基因的發(fā)掘提供參考序列。尤其是通過GO、COG和KEGG功能分類及代謝通路等分析，不僅可為黃唇魚相關(guān)基因克隆及功能分析等提供基礎(chǔ)數(shù)據(jù)，還能為其保護(hù)遺傳學(xué)研究提供大量可利用的基因資源。

轉(zhuǎn)錄組測序是開發(fā)SSR分子標(biāo)記的一種有效方法，目前已從橄欖蚌（Solenaia oleivor）（徐艷，2014）、施氏鱘（Acipenser schrenckii）（孔杰等，2015）、翹嘴鱖（Siniperca chuatsi）（袁文成，2015）、安吉小鯢（Hynobius amjiensis）（王宇等，2017）和中華蜜蜂（Apiscerana cerana）（熊翠玲等，2017）的轉(zhuǎn)錄組成功開發(fā)出大量SSR分子標(biāo)記。本研究通過對黃唇魚進(jìn)行轉(zhuǎn)錄組測序，檢測到29797個SSR位點(diǎn)分布于19664條Unigenes中。黃唇魚高通量測序結(jié)果顯示這些SSR位點(diǎn)類型呈多樣性，包含單核苷酸～六核苷酸重復(fù)基元，且不同重復(fù)類型的數(shù)量差異明顯，其中，最多的重復(fù)類型為二核苷酸重復(fù)基元，其次是單核苷酸和三核苷酸重復(fù)基元。在二核苷酸重復(fù)基元類型中重復(fù)最多的基序是AC/GT，占二核苷酸重復(fù)基元的70.48%，與在團(tuán)頭魴（Megalobrama amblycephala）（曾聰?shù)龋?013）和翹嘴鱖（S. chuatsi）（袁文成，2015）上的研究結(jié)果相同，但與許多植物以單核苷酸重復(fù)基元為主的結(jié)論不同，說明不同物種的堿基重復(fù)基序存在種屬差異。

本研究從隨機(jī)挑選的186對SSR引物中篩選出121對能擴(kuò)增出目的條帶的SSR引物，擴(kuò)增成功率達(dá)65.05%。進(jìn)一步進(jìn)行測序分析和多態(tài)性檢測，結(jié)果篩選出47對擴(kuò)增穩(wěn)定性好且具有多態(tài)性的SSR引物，多態(tài)性比例為38.84%。本研究中黃唇魚SSR多態(tài)性引物擴(kuò)增效率不高的原因可能是：（1）隨機(jī)篩選的SSR引物存在不均性;（2）某些序列拼接錯誤或引物設(shè)計質(zhì)量不高;（3）用于多態(tài)性驗(yàn)證的黃唇魚樣本數(shù)量較少導(dǎo)致群體等位基因數(shù)和遺傳多樣性相應(yīng)較低，但具體原因有待進(jìn)一步探究。

4 結(jié)論

通過轉(zhuǎn)錄組測序能有效篩選出具有多樣性的黃唇魚SSR分子標(biāo)記，可為黃唇魚遺傳多樣性分析、親緣鑒定和遺傳圖譜構(gòu)建等后續(xù)研究提供基礎(chǔ)資料。

參考文獻(xiàn)：

龔詩琦，王志勇，肖世俊，林愛強(qiáng)，謝仰杰. 2016. 黃姑魚轉(zhuǎn)錄組SSR的開發(fā)與驗(yàn)證[J]. 集美大學(xué)學(xué)報（自然科學(xué)版），21（4）：241-246. [Gong S Q，Wang Z Y，Xiao S J，Lin A Q，Xie Y J. 2016. Development and verification of SSR based on transcription of yellow drum，Nibea albiflora[J]. Journal of Jimei University（Natural Science），21（4）：241-246.]

何福玲，凌正寶，肖俊，郭忠寶，鐘歡，楊弘，陳文治，唐瞻楊，單丹，梁軍能，甘西，羅永巨. 2017. 兩個人工選擇奧利亞羅非魚群體系統(tǒng)發(fā)育及其遺傳多樣性分析[J]. 南方農(nóng)業(yè)學(xué)報，48（2）：341-349. [He F L，Ling Z B，Xiao J，Guo Z B，Zhong H，Yang H，Chen W Z，Tang Z Y，Shan D，Liang J N，Gan X，Luo Y J. 2017. Phylogeny and genetic diversity of two artificial selection Oreochromis aureus populations[J]. Journal of Southern Agriculture，48（2）：341-349.]

孔杰，石婧，向燕，關(guān)梅. 2015. 基于微衛(wèi)星標(biāo)記的養(yǎng)殖施氏鱘親魚群體遺傳多樣性[J]. 貴州農(nóng)業(yè)科學(xué)，43（12）：119-122. [Kong J，Shi J，Xiang Y，Guan M. 2015. Genetic diversity of Acipenser schrenckii from hatchery with microsatellite to genetic markers[J]. Guizhou Agricultural Scien-ces，43（12）：119-122.]

區(qū)又君，廖銳，李加兒，勾效偉. 2010. 野生黃唇魚Bahaba flavolabiata肌肉和鰾營養(yǎng)成分的分析與評價[J]. 海洋學(xué)報，32（3）：113-120. [Ou Y J，Liao R，Li J E，Gou X W. 2010. The analysis and evaluation of nutrition composition in muscle and air bladder of wild Bahaba flavolabiata[J]. Acta Oceanologica Sinica，32（3）：114-119.]

區(qū)又君，廖銳，李加兒，勾效偉. 2011. 黃唇魚的耳石形態(tài)和微結(jié)構(gòu)特征[J]. 廣東農(nóng)業(yè)科學(xué)，38（12）：123-124. [Ou Y J，Liao R，Li J E，Gou X W. 2011. The morphology of otolith and its characteristics of microstructure in Bahaba flavolabiata[J]. Guangdong Agricultural Science，38（12）：123-124.]

區(qū)又君，廖銳，李加兒，勾效偉. 2013. 黃唇魚矢耳石的生長特性及與魚體生長的關(guān)系[J]. 江蘇農(nóng)業(yè)科學(xué)，41（6）：189-192. [Ou Y J，Liao R，Li J E，Gou X W. 2013. The growth characteristics of sagittal otolith of Bahaba flavolabiata and its relationship with the growth of fish body[J]. Jiangsu Agricultural Science，41（6）：189-192.]

區(qū)又君，廖銳，李加兒，勾效偉. 2014. 采用耳石日輪研究珠江口黃唇魚幼魚的產(chǎn)卵期及生長[J]. 生態(tài)科學(xué)，33（2）：209-212. [Ou Y J，Liao R，Li J E，Gou X W. 2014. Stu-dies on the spawning period and growth of young Bahaba flavolabiata in estuary of Pearl River by otolith determination[J]. Ecological Science，33（2）：209-212.]

區(qū)又君，廖銳，李加兒，李劉冬，勾效偉. 2012. 4種石首魚耳石形態(tài)特征的比較[J]. 華南農(nóng)業(yè)大學(xué)學(xué)報，33（2）：203-210. [Ou Y J，Liao R，Li J E，Li L D，Gou X W. 2012. Comparison of morphological characteristics of otolith in four Sciaenid fishes[J]. Journal of South China Agricultural University，33（2）：203-210.]

孫成飛，葉星，董浚鍵，田園園，梁健輝. 2015. 羅氏沼蝦6個養(yǎng)殖群體遺傳多樣性的微衛(wèi)星分析[J]. 南方水產(chǎn)科學(xué)，11（2）：20-26. [Sun C F，Ye X，Dong J J，Tian Y Y，Liang J H. 2015. Genetic diversity analysis of six cultured populations of Macrobrachium rosenbergii using microsatellite markers[J]. South China Fisheries Science，11（2）：20-26.]

孫效文，張曉鋒，趙瑩瑩，張研，賈智英，常玉梅，魯翠云，梁利群. 2008. 水產(chǎn)生物微衛(wèi)星標(biāo)記技術(shù)研究進(jìn)展及其應(yīng)用[J]. 中國水產(chǎn)科學(xué)，15（4）：689-703. [Sun X W，Zhang X F，Zhao Y Y，Zhang Y，Jia Z Y，Chang Y M，Lu C Y，Liang L Q. 2008. Development and application of microsatellite markers in aquatic species[J]. Journal of Fishery Sciences of China，15（4）：689-703.]

王玲玲，陳東亮，黃叢林，邢震. 2017. SSR分子標(biāo)記技術(shù)在植物研究中的應(yīng)用[J]. 安徽農(nóng)業(yè)科學(xué)，45（36）：123-126. [Wang L L，Chen D L，Huang C L，Xing Z. 2017. Application of SSR molecular marker technique in plant research[J]. Journal of Anhui Agricultural Sciences，45（36）：123-126.]

王明，應(yīng)東山，王琴飛，李莉萍，黃建鋒，張如蓮. 2015. 基于454 GSFLX高通量測序的芒果微衛(wèi)星標(biāo)記特征分析[J]. 南方農(nóng)業(yè)學(xué)報，46（10）：1745-1749. [Wang M，Ying D S，Wang Q F，Li L P，Huang J F，Zhang R L. 2015. Deve-lopment and characteristics analysis of microsatellite mar-ker in mango（Mangifera indica L.） by using 454 GSFLX high-throughput sequencing[J]. Journal of Southern Agriculture，46（10）：1745-1749.]

王宇，李丹丹，邵晨. 2017. 安吉小鯢微衛(wèi)星標(biāo)記及其應(yīng)用：CN106868186A[P]. 2017-06-20. [Wang Y，Li D D，Shao C. 2017. Microsatellite marker of Hynobius amjiensis and its application：CN106868186A[P]. 2017-06-20.]

王忠華，董西征，錢國英. 2008. DNA分子標(biāo)記技術(shù)在水生動物遺傳多樣性研究中的應(yīng)用[J]. 科技通報，24（5）：623-630. [Wang Z H，Dong X Z，Qian G Y. 2008. Application in genetic diversity of aquatic of DNA molecular markers[J]. Bulletin of Science and Technology，24（5）：623-630.]

謝子強(qiáng)，劉雪媚， 郭昱嵩，劉楚吾. 2013. 北部灣皮氏叫姑魚的遺傳多樣性分析[J]. 南方農(nóng)業(yè)學(xué)報，44（1）：140-144. [Xie Z Q，Liu X M，Guo Y S，Liu C W. 2013. Genetic diversity analysis of Beibu Gulf Johnius belengerii[J]. Journal of Southern Agriculture，44（1）：140-144.]

熊翠玲，張璐，付中民，王鴻權(quán)，侯志賢，童新宇，李汶東，鄭燕珍，陳大福，郭睿. 2017. 基于RNA-seq數(shù)據(jù)大規(guī)模開發(fā)中華蜜蜂幼蟲的SSR分子標(biāo)記[J]. 環(huán)境昆蟲學(xué)報，39（1）：68-74. [Xiong C L，Zhang L，F(xiàn)u Z M，Wang H Q，Hou Z X，Tong X Y，Li W D，Zheng Y Z，Chen D F，Guo R. 2017. Large-scale development of SSR primers for Apis cerana cerana larvae based on its RNA-seq datasets[J]. Journal of Environmental Entomology，39（1）：68-74.]

熊良偉，王帥兵，岳麗佳，王建國，陶桂慶，徐亮，王權(quán). 2018. 寬體金線蛭基因組SSR序列特征分析及其分子標(biāo)記開發(fā)[J]. 南方農(nóng)業(yè)學(xué)報，49（11）：2298-2303. [Xiong L W，Wang S B，Yue L J，Wang J G，Tao G Q，Xu L，Wang Q. 2018. SSR sequence characters for genome of Whitmania pigra Whitman and development of molecular markers[J]. Journal of Southern Agriculture，49（11）：2298-2303.]

徐艷. 2014. 橄欖蟶蚌微衛(wèi)星標(biāo)記開發(fā)及群體遺傳多樣性分析[D]. 武漢：華中農(nóng)業(yè)大學(xué). [Xu Y. 2014. Isolation of microsatellite markers and population genetic diversity analysis in Solenaia oleivora[D]. Wuhan：Huazhong Agricultural University.]

楊兵，林琳，李純厚，徐姍楠，劉永，肖雅元，陳作志. 2015. 基于高通量測序的二長棘鯛微衛(wèi)星標(biāo)記開發(fā)與評價[J]. 南方水產(chǎn)科學(xué)，11（4）：116-120. [Yang B，Lin L，Li C H，Xu S N，Liu Y，Xiao Y Y，Chen Z Z. 2015. Development and evaluation of microsatellite markers in Parargyrops edita[J]. South China Fisheries Science，11（4）：116-120.]

袁文成. 2015. 基于轉(zhuǎn)錄組測序的翹嘴鱖微衛(wèi)星標(biāo)記的開發(fā)及MHC classI基因的克隆表達(dá)[D]. 蘇州：蘇州大學(xué). [Yuan W C. 2015. Development of SSR makers and MHC classI cloning and expressing of Siniperca chuatsi by RNA-Seq[D]. Suzhou：Soochow University.]

趙志英，梁麗運(yùn)，白麗蓉. 2018. 斑節(jié)對蝦3個野生群體遺傳多樣性的微衛(wèi)星標(biāo)記分析[J]. 熱帶海洋學(xué)報，37（3）：65-72. [Zhao Z Y，Liang L Y，Bai L R. 2018. Analysis of genetic diversity among three wild populations of Penaeus monodon using microsatellite marker[J]. Journal of Tro-pical Oceanography，37（3）：65-72.]

曾聰，高澤霞，羅偉，劉肖蓮，王衛(wèi)民. 2013. 基于454 GS FLX高通量測序的團(tuán)頭魴ESTs中微衛(wèi)星特征分析[J]. 水生生物學(xué)報，37（5）：982-988. [Zeng C，Gao Z X，Luo W，Liu X L，Wang W M. 2013. Characteristics of microsate-llites in blunt snout bream（Megalobrama amblycephala） EST sequence using 454 FLX[J]. Acta Hydrobiologica Sinica，37（5）：982-988.]

Clark T B，Ma L，Saillant E，Gold J R. 2004. Microsatellite DNA markers for population-genetic studies of Atlantic bluefin tuna（Thunnus thynnus thynnus） and other species of genus Thunnus[J]. Molecular Ecology Notes，4（1）：70-73.

Furukawa S，Takeshima H，Otaka T，Mitsuboshi T，Shirasu K，Ikeda D，Kaneko G，Nishida M，Watabe S. 2004. Isolation of microsatellite markers by in silico screening implicated for genetic linkage mapping in Japanese pufferfish Takifugu rubripes[J]. Fisheries Science，70（4）：620-628.

King T L，Eackles M S，Letcher B H. 2005. Microsatellite DNA markers for the study of Atlantic salmon（Salmo salar） kinship，population structure，and mixed fishery ana-lyses[J]. Molecular Ecology Notes，5（1）：130-132.

Pardo G B，Casas L，F(xiàn)ortes G G，Bouza C，Martínez P，Clark M S，Sánchez L. 2005. New microsatellite markers in turbot（Scophthalmus maximus） derived from an enriched genomic library and sequence databases[J]. Molecular Eco-logy Notes，5（1）：62-64.

Rivera M A J，Graham G C，Roderick G K. 2003. Isolation and characterization of nine microsatellite loci from the Hawaiian grouper Epinephelus quernus（Serranidae） for po-pulation genetic analyses[J]. Marine Biotechnology，5（2）：126-129.

Zane L，Bargelloni L，Patarnello T. 2002. Strategies for microsatellite isolation：A review[J]. Molecular Ecology，11（1）：1-16.

（責(zé)任編輯 蘭宗寶）