董忠典 黎學(xué)友 廖 健 張 寧 郭昱嵩 王中鐸

(廣東海洋大學(xué)水產(chǎn)學(xué)院 南海水產(chǎn)經(jīng)濟(jì)動(dòng)物增養(yǎng)殖重點(diǎn)實(shí)驗(yàn)室 湛江 524088)

肝臟是脊椎動(dòng)物最重要的器官之一， 在物質(zhì)能量代謝、激素合成及免疫反應(yīng)等生理過(guò)程中具有重要意義(van der Oostet al， 2003； Tomet al， 2004； Laiet al， 2015； Zhanget al， 2017)。在魚(yú)類中， 雌性肝臟還是卵黃蛋白原(vitellogenins，vtgs)、透明帶蛋白前體(choriogenins，chgs)、雌激素受體(estrogen receptor，er)和細(xì)胞色素P450(cytochrome P450s，cyp450s)等基因的主要表達(dá)器官， 這些基因?qū)Νh(huán)境中雌激素物質(zhì)較為敏感， 當(dāng)水體環(huán)境中存在雌激素類物質(zhì)時(shí)， 可以誘導(dǎo)雄魚(yú)肝臟表達(dá)這些基因(Leeet al， 2002； Chenet al， 2008； Laiet al， 2015)。

弓背青鳉(Oryzias curvinotus)是一種廣泛分布于我國(guó)南海沿岸水域的小型魚(yú)類， 為紅樹(shù)林水域的常見(jiàn)物種(Hayakawaet al， 2015； Wanget al， 2017)。弓背青鳉具有體型小、生長(zhǎng)快、性成熟時(shí)間短和胚胎發(fā)育透明的特點(diǎn)， 可以用于實(shí)驗(yàn)室飼養(yǎng)繁殖； 其對(duì)鹽度有較強(qiáng)的適應(yīng)能力， 具備開(kāi)發(fā)成海洋環(huán)境毒理及鹽度適應(yīng)研究新模式種的潛力(Hayakawaet al， 2015)。同時(shí)， 弓背青鳉的性染色體為XX/XY 型， 具有性別標(biāo)記， 是研究海洋魚(yú)類性別決定和分化機(jī)制的潛在模型(Matsudaet al， 2003； 董忠典等， 2018)。同屬的日本青鳉(Oryziaslatipes)和海水青鳉(Oryziasmelastigma)已被廣泛應(yīng)用于淡水和海水環(huán)境檢測(cè)研究(Metcalfeet al， 2000； Chenet al， 2008； Honget al， 2015； Conget al， 2017； Abdel-moneimet al， 2018)。然而， 缺乏基因組和轉(zhuǎn)錄組信息， 在很大程度上限制了弓背青鳉資源的開(kāi)發(fā)和利用。本研究采用RNA-Seq 技術(shù)分別對(duì)雌、雄弓背青鳉肝臟進(jìn)行了轉(zhuǎn)錄組測(cè)序、基因功能注釋和表達(dá)譜分析， 以期為弓背青鳉基因資源的開(kāi)發(fā)提供基礎(chǔ)數(shù)據(jù)。

1 材料與方法

1.1 弓背青鳉的采集及總肝臟RNA 提取

本研究所用性成熟弓背青鳉(平均體長(zhǎng)29.25mm， 平均體質(zhì)量0.27g)采集自湛江雷州半島近岸的高橋國(guó)家級(jí)紅樹(shù)林保護(hù)區(qū)(21°36′24"N， 109°47′8"E)， 實(shí)驗(yàn)室 暫 養(yǎng) 30d， 溫 度 25±1°C， 鹽 度 15， 光 暗 周 期14h：10h。解剖通過(guò)性腺類型確定弓背青鳉生理性別， 雌雄各 3 尾取肝臟分別進(jìn)行總 RNA 提取， 使用DNase I(寶生物， 大連)去除基因組污染。使用1%瓊脂糖凝膠電泳分析RNA 完整性， Nanodrop 檢測(cè)RNA的純度(OD260/OD280比值)， Qubit 檢測(cè)總RNA 濃度， RNA 滿足進(jìn)行下一步的文庫(kù)構(gòu)建的要求。對(duì)3 尾雌性肝臟總RNA， 分別取等量(1.0μg)混合為一個(gè)RNA樣品池(FL)， 用同樣的方法獲得雄性肝臟RNA 樣品池(ML)， 對(duì)FL 和ML 各取1.5μg 總RNA， 分別進(jìn)行cDNA 文庫(kù)構(gòu)建和轉(zhuǎn)錄組測(cè)序。

1.2 文庫(kù)構(gòu)建及測(cè)序

使用帶Oligo(dT)的磁珠富集真核生物mRNA， 加入fragmentation buffer 將mRNA 打斷成短片段， 以mRNA 為模板， 用六堿基隨機(jī)引物合成第一鏈cDNA， 然后加入緩沖液、dNTPs 和DNA 聚合酶和RNase H合成第二鏈cDNA， 再用AMPure XP beads 純化雙鏈cDNA。純化的雙鏈cDNA 先進(jìn)行末端修復(fù)、加A 尾并連接測(cè)序接頭， 再用 AMPure XP beads 選擇150—200bp 片段。隨后進(jìn)行PCR 擴(kuò)增， 并用AMPure XP beads 純化PCR 產(chǎn)物， 得到最終的文庫(kù)。文庫(kù)構(gòu)建完成后， 用Qubit2.0 進(jìn)行初步定量并稀釋至1.5ng/μL， 用 Agilent 2100 確認(rèn)文庫(kù)插入片段大小， 最后用Q-PCR 方法對(duì)文庫(kù)的有效濃度進(jìn)行準(zhǔn)確定量(文庫(kù)有效濃度＞2nmol/L)， 以保證文庫(kù)質(zhì)量。庫(kù)檢合格后， 進(jìn)行Illumina HiSeq 測(cè)序， cDNA 文庫(kù)的制備和測(cè)序由諾禾致源生物信息科技有限公司(天津)完成。

1.3 轉(zhuǎn)錄組測(cè)序質(zhì)量控制和組裝

為了保證轉(zhuǎn)錄組分析質(zhì)量， 對(duì)原始測(cè)序數(shù)據(jù)(raw reads)過(guò)濾， 去除帶接頭的reads、N(N 表示無(wú)法確定堿基信息)的比例大于0.1%的reads、低質(zhì)量(質(zhì)量值Qphred＜=20 的堿基數(shù)占整個(gè)reads 的50%以上)的reads， 得到高質(zhì)量測(cè)序數(shù)據(jù)(clean reads)。采用Trinity軟件對(duì)clean reads 進(jìn)行拼接(Grabherret al， 2011)。

1.4 基因功能注釋

為獲得全面的基因功能信息， 使用7 大數(shù)據(jù)庫(kù)對(duì)Trinity 拼接獲得的unigenes 進(jìn)行功能注釋。使用在線BLAST 程序?qū)nigenes 在Nr (http：//www.ncbi.nlm. nih.gov)期望值e-value ＜ e-5、Nt (http：//www.ncbi.nlm. nih.gov) e-value ＜ e-5、Swiss-Prot (http：//www.expasy. ch/sprot) e-value ＜ e-5、KOG (http：//www.ncbi.nlm.nih. gov) e-value ＜ e-3、KO (https：//www.kegg.jp/kegg/ko. html) e-value ＜ e-10和GO (http：//www.geneontology. org/) e-value ＜ e-6進(jìn)行比對(duì)注釋。通過(guò)HMMER3.0 (e-value ＜ 0.01)程序?qū)nigenes 在Pfam 數(shù)據(jù)庫(kù)中進(jìn)行蛋白家族注釋(Eddy， 2011)。

1.5 基因表達(dá)水平和差異分析

將Trinity 拼接得到的轉(zhuǎn)錄組作為參考序列， 通過(guò)RSEM 對(duì)每個(gè)樣本的基因表達(dá)水平進(jìn)行估計(jì)： 將Clean data 比對(duì)到組裝好的參考序列上， 根據(jù)比對(duì)結(jié)果得到每個(gè)基因的Readcount 數(shù)目， 采用TMM 對(duì)read count 數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理， 再用DEGseq 進(jìn)行差異分析， 篩選閾值為qvalue ＜0.005 且|log2(foldchange)|＞1 (Wanget al， 2010)。

1.6 差異基因GO 和KEGG 富集分析

使用R 軟件中的GOseq 對(duì)弓背青鳉雌雄肝臟轉(zhuǎn)錄組差異基因進(jìn)行顯著性富集分析(Younget al， 2010)。使用KOBAS 軟件檢測(cè)差異基因在KEGG 通路中富集情況(Maoet al， 2005)。

1.7 實(shí)時(shí)熒光定量PCR(RT-qPCR)驗(yàn)證

利用RT-qPCR 技術(shù)， 弓背青鳉雌雄肝臟轉(zhuǎn)錄組中選擇 12 個(gè)基因進(jìn)行定量驗(yàn)證(引物見(jiàn)表 1)。RT-qPCR 反應(yīng)體系15μL， 包含7.5μL 2×Power Green qPCR Mix (東盛生物， 廣州)， 0.6μL 上下游引物(10μmol/L)， 1.5μL cDNA， 4.8μL ddH2O。每個(gè)樣品技術(shù)重復(fù)3 次， 反應(yīng)置于Roche LightCycler 96(羅氏， 瑞士)上運(yùn)行。反應(yīng)程序如下： 95°C 3min； 95°C 15s， 60°C 15s， 72°C 30s (采集熒光)， 40 個(gè)循環(huán)； 熔解曲線分析檢測(cè)PCR 產(chǎn)物的特異性。gapdh和adp作為內(nèi)參基因， 使用SPSS17.0 采用2-ΔΔCt方法對(duì)目的基因表達(dá)水平進(jìn)行統(tǒng)計(jì)(Livaket al， 2001)。

表1 本研究所用引物及序列 Tab.1 The sequences of primers used in this study

續(xù)表

1.8 微衛(wèi)星標(biāo)記(SSR)檢測(cè)

采 用 MISA (http：//pgrc.ipk-gatersleben.de/misa/ misa.html)對(duì)弓背青鳉肝臟轉(zhuǎn)錄組拼接獲得的unigenes進(jìn)行微衛(wèi)星標(biāo)記(SSR)檢測(cè)， 設(shè)置各重復(fù)單元類型的最少重復(fù)次數(shù)為： 1—10、2—6、3—5、4—5、5—5、6—5。以1—10 為例， 該設(shè)置表示單核苷酸重復(fù)類型至少重復(fù)10 次才被算為微衛(wèi)星。

2 結(jié)果

2.1 測(cè)序結(jié)果與組裝

本研究分別構(gòu)建了雌、雄弓背青鳉肝臟cDNA 文庫(kù)， 并進(jìn)行了轉(zhuǎn)錄組測(cè)序。測(cè)序結(jié)果經(jīng)過(guò)質(zhì)控后， 分別獲得80095044 和87043984 條Clean reads。經(jīng)Trinity 組裝， 共得到49912 個(gè)unigenes， 其中N50 長(zhǎng)度為2394bp， 其中長(zhǎng)度大于1000bp 的unigenes 共計(jì)25423 個(gè)(圖1)。

圖1 unigenes 序列長(zhǎng)度分布圖 Fig.1 The distribution of sequence length of the unigenes

2.2 轉(zhuǎn)錄組注釋

通過(guò)BLAST 比對(duì)， 分別在NR、NT、SwissProt、PFAM、KOG 數(shù)據(jù)庫(kù)中注釋到32651、47988、19005、25587 和19005 個(gè)unigenes。與NR 數(shù)據(jù)庫(kù)比對(duì)發(fā)現(xiàn)有32649 個(gè)unigenes 與390 個(gè)物種基因序列高度同源， 其中77.6%(25349 個(gè))的unigenes 和日本青鳉基因序列同源， 其次有 3.6%的 unigenes 和深裂眶鋸雀鯛(Stegastes partitus)同源(圖2)。為進(jìn)一步了解弓背青鳉雌、雄肝臟表達(dá)基因集的功能， 將拼接的unigenes在GO 和KO 數(shù)據(jù)庫(kù)中進(jìn)行了比對(duì)分析。25630 個(gè)unigenes 被歸類到GO 三個(gè)大類(生物學(xué)過(guò)程、分子功能和細(xì)胞成分)55 個(gè)功能分類中(圖3)。生物學(xué)過(guò)程包 括25 個(gè)功能分類， 參與“細(xì)胞過(guò)程”(15094)和“代謝過(guò)程”(12821)的基因最多； 細(xì)胞成分包含20 個(gè)功能分類，“細(xì)胞”(8419)和“細(xì)胞部分”(8419)基因最多； 分子功能包含10 個(gè)功能分類， “結(jié)合”(16155)和“催化活性”(10114)基因最多。對(duì)基因做KO 注釋后， 可根據(jù)它們參與的KEGG 代謝通路進(jìn)行分類。共有18887個(gè)unigenes 在KO 數(shù)據(jù)庫(kù)中獲得注釋， 富集到232 個(gè)KEGG 信號(hào)通路中， 其中“PI3K-Akt”信號(hào)通路含有最多(657)的 unigenes， 其次是“內(nèi)吞”(615)、“斑黏連”(521)和“肌動(dòng)蛋白骨架”的調(diào)節(jié)通路(508)。

圖2 弓背青鳉unigenes 同源物種分布 Fig.2 Distribution of homologous species of O. curvinotus unigenes

2.3 差異基因表達(dá)及富集分析

通過(guò)DEGseq 對(duì)弓背青鳉雌、雄肝臟轉(zhuǎn)錄組基因進(jìn)行差異表達(dá)分析(差異基因篩選條件為： qvalue ＜ 0.005 & |log2(foldchange)| ＞ 1)， 結(jié)果如圖4 所示， 有207 個(gè)unigenes 在雌魚(yú)肝臟中上調(diào)， 364 個(gè)unigenes 在雄魚(yú)肝臟上調(diào)。對(duì)差異基因(DEGs)進(jìn)行GO 和KEGG 富集分析， GO 富集結(jié)果顯示雌性肝臟高表達(dá)基因主要 參與蛋白質(zhì)合成過(guò)程， 雄性高表達(dá)基因主要參與免疫及氧化還原反應(yīng)(圖5)。KEGG 分析結(jié)果顯示DEGs 富集到244 個(gè)KEGG 通路中， 涉及核糖體、PPAR 信號(hào)通路、卵巢類固醇激素生成、類固醇生物合成、不飽和脂肪酸合成和內(nèi)質(zhì)網(wǎng)蛋白質(zhì)加工等通路(圖6)。

2.4 RT-qPCR 驗(yàn)證

本研究選擇了在雌、雄弓背青鳉肝臟RNA-Seq分析中的12 個(gè)差異表達(dá)基因， 通過(guò)RT-qPCR 技術(shù)對(duì)RNA-Seq 結(jié)果進(jìn)行驗(yàn)證。所選的驗(yàn)證基因包括環(huán)腺苷三磷酸依賴的轉(zhuǎn)錄因子camp， 內(nèi)質(zhì)網(wǎng)基因endoplasmin， UDP-葡萄糖醛酸轉(zhuǎn)移酶基因utp， 磺基轉(zhuǎn)移酶家族2B 基因sult2b及部分檢測(cè)環(huán)境雌激素物質(zhì)的標(biāo)記基因如卵黃蛋白原基因vtg1、vtg2、vtglike， 透明帶蛋白前體基因chgl、chgh、chghm， 細(xì)胞色素P450 基因p4502k1、p45027c。結(jié)果表明所選基因的表達(dá)模式與RNA-Seq 分析相一致， 說(shuō)明RNA-Seq 分析結(jié)果可信(見(jiàn)圖7)。

圖6 弓背青鳉雌雄肝臟差異基因KEGG 富集 Fig.6 KEGG assignment of differentially expressed genes of O. curvinotus

圖7 RT-qPCR 檢測(cè)RNA-Seq 結(jié)果 Fig.7 Validation of RNA-Seq data using RT-qPCR

2.5 簡(jiǎn)單重復(fù)序列SSR 分析

利用MISA 軟件對(duì)弓背青鳉進(jìn)行SSR 標(biāo)記鑒定， 用于遺傳多樣性研究。有18326 條unigenes 序列含有SSR 標(biāo)記， 共計(jì)28195 個(gè)， 其中6373 條unigenes 包含一個(gè)以上的SSR。各種類型的SSR 標(biāo)記出現(xiàn)的頻率不同(表2)， 單核苷酸重復(fù)SSR 出現(xiàn)的頻率最高， 占總SSR 的73.17%， 其次是三核苷酸重復(fù)SSR， 占總SSR 的19%。另外， 不同重復(fù)序列的相對(duì)豐度差異很大。A/T 在單核苷酸SSR 中最常見(jiàn)， 在雙核苷酸SSR中， AG/CT 最為常見(jiàn)。AGG/CCT 和AAAC/GTTT 是三、四核苷酸SSR 中最常見(jiàn)的元件。

3 討論

3.1 弓背青轉(zhuǎn)錄組測(cè)序分析

為了闡明弓背青鳉肝臟表達(dá)基因的功能和涉及的生物學(xué)過(guò)程， 本研究分別構(gòu)建了性成熟雌、雄弓背青鳉肝臟 cDNA 文庫(kù)， 并在 Illumina 平臺(tái)進(jìn)行了RNA-Seq 測(cè)序， 分別獲得了80095044 和87043984 條Clean reads。N50 長(zhǎng)度是評(píng)價(jià)RNA-Seq 組裝質(zhì)量的重要參數(shù)， 本研究弓背青鳉肝臟轉(zhuǎn)錄組N50長(zhǎng)度為2394bp， 高于多個(gè)已發(fā)表的水生動(dòng)物肝臟轉(zhuǎn)錄組， 包括海水青鳉(2162bp) (Laiet al， 2015)， 七鰓鰻(lampetra japonica)(1447bp) (李 慶 偉 等， 2018)， 脂 鯉(Pterygoplichthys anisitsi) (1571bp) (Parenteet al， 2017)， 低 于 海 鱸(Dicentrarchus labrax) (3257bp) (Magnanouet al， 2014)。對(duì)NR 數(shù)據(jù)庫(kù)比對(duì)， 發(fā)現(xiàn)有77.64% (32649 條unigenes 中的25349 個(gè))的可比對(duì)序列與日本青鳉基因同源(圖2)， 可能是因?yàn)楣城圜毢腿毡厩圜毦鶎儆谇圜殞伲?親緣關(guān)系密切(Wanget al， 2017)。

表2 弓背青鳉肝臟轉(zhuǎn)錄組SSR 類型及數(shù)量 Tab.2 SSR summary from the O. curvinotus liver transcriptome

肝臟在魚(yú)類生長(zhǎng)和物質(zhì)能量代謝中起著重要的作用。在弓背青鳉肝臟轉(zhuǎn)錄組中， “細(xì)胞過(guò)程”、“細(xì)胞”和“結(jié)合”分別在GO 三大功能類別中富集(圖3)。這與其他水生動(dòng)物的肝臟RNA-Seq 結(jié)果類似， 說(shuō)明屬于這幾類功能的基因在物種間保守(Magnanouet al， 2014； Chenet al， 2016； 張燕萍等， 2018)。KEGG 的注釋結(jié)果表明， 弓背青鳉肝臟轉(zhuǎn)錄組中unigenes 主要參與了信號(hào)轉(zhuǎn)導(dǎo)(3130)、內(nèi)分泌系統(tǒng)(1536)及免疫系統(tǒng)(1491)的生物學(xué)過(guò)程。結(jié)果可能提示了肝臟在弓背青鳉生命過(guò)程中的主要功能， 為了解肝臟的功能演變奠定了基礎(chǔ)。

3.2 弓背青 免疫相關(guān)基因

肝臟是魚(yú)類的重要免疫器官， 通過(guò)KEGG 分析， 從弓背青鳉肝臟轉(zhuǎn)錄組中篩選出了多個(gè)可能參與免疫相關(guān)防御途徑的unigenes(表3)： 373 個(gè)unigenes 參與了血小板活化， 272 個(gè)unigenes 參與了白細(xì)胞跨內(nèi)皮遷移， 245 個(gè)unigenes 參與了T 細(xì)胞受體信號(hào)通路， 136 個(gè)unigenes 參與了補(bǔ)體和凝血級(jí)聯(lián)， 188 個(gè)unigenes 參與了自然殺傷細(xì)胞介導(dǎo)的細(xì)胞毒性， 163 個(gè)unigenes 參與了Toll-樣受體信號(hào)通路等過(guò)程。在七鰓鰻的肝臟轉(zhuǎn)錄組中也有相似的結(jié)果(李慶偉等， 2018)。

表3 弓背青鳉肝臟轉(zhuǎn)錄組中參與免疫反應(yīng)相關(guān)KEGG 通路的unigenes Tab.3 KEGG pathways with immune-related genes enrichment in the liver transcriptomes of O. curvinotus

血小板上具有toll 樣受體， 在機(jī)體的先天免疫應(yīng)答中起著關(guān)鍵作用。當(dāng)病原菌與血小板結(jié)合時(shí)， 血小板活化同時(shí)分泌抗菌肽， 引起機(jī)體免疫反應(yīng)(Coxet al， 2011)。白細(xì)胞跨內(nèi)皮遷移在先天性和適應(yīng)性免疫系統(tǒng)應(yīng)答中起重要作用(Muller， 2011)。當(dāng)機(jī)體受損時(shí)， 就會(huì)發(fā)生炎癥反應(yīng)， 它涉及血液中預(yù)先形成的可溶性元素迅速而短暫地運(yùn)送到受傷部位， 隨后是較長(zhǎng)時(shí)間的白細(xì)胞運(yùn)送， 經(jīng)內(nèi)皮細(xì)胞遷移可能是炎癥反應(yīng)的不可逆點(diǎn)(Jimenezet al， 2010； Muller， 2011)。補(bǔ)體和凝集通路在機(jī)體固有免疫防御中具有重要意義(Bajicet al， 2016)， 在 小 瓜 蟲(chóng) 感 染 的 大 黃 魚(yú)(Larimichthys crocea)早期免疫中發(fā)揮重要作用； 另外， 補(bǔ)體和凝集通路在鰻弧菌感染后半滑舌鰨的免疫反應(yīng)中也發(fā)揮了重要作用(Zhanget al， 2015； Yinet al， 2016)。綜上， 我們推測(cè)弓背青鳉的肝臟可能在抵御病原體入侵的免疫防御中發(fā)揮重要作用。

3.3 環(huán)境內(nèi)分泌干擾物的生物標(biāo)志物

魚(yú)類肝臟常被用作檢測(cè)環(huán)境中雌激素物質(zhì)的靶器官， 通常用于檢測(cè)環(huán)境雌激素物質(zhì)的標(biāo)記基因在雌魚(yú)肝臟中大量表達(dá)， 如chgs、vtg、cyp450s和er(Bucheliet al， 1995； Arukweet al， 2001； Leeet al， 2002； Chenet al， 2008)等。Chgs基因編碼放射帶蛋白(ZPs)的前體， ZPs 占據(jù)了魚(yú)卵殼的大部分， 由ZP-1、2和ZP-3 兩個(gè)主要亞基群組成(Hamazakiet al，1989； Murataet al， 1991)。在性成熟雌魚(yú)中， 雌激素誘導(dǎo)Chgs 在肝臟中合成， 隨血液流動(dòng)整合到放射帶中。當(dāng)雄魚(yú)被環(huán)境雌激素物質(zhì)刺激時(shí)， 肝臟也產(chǎn)生Chgs。在日本青鳉中， 雙酚A、壬基酚和乙炔雌二醇均可以誘導(dǎo)Chgl 和Chgh 的表達(dá)， 表達(dá)量的增加與雌激素劑量呈正相關(guān)(Leeet al， 2002)。Chen 等(2008)研究表明海水青鳉chgl和chgh基因?qū)Νh(huán)境雌激素物質(zhì)也有相似的響應(yīng)。卵黃蛋白原是硬骨魚(yú)類卵黃蛋白前體， 雖然Vtgs 在成熟的雌性體內(nèi)自然存在， 而在雄性體內(nèi)則不存在， 但雌性激素或環(huán)境雌激素分泌物可以誘導(dǎo)雄性產(chǎn)生Vtgs。因此， 雄性肝臟中的vtgs 轉(zhuǎn)錄產(chǎn)物， 可用于指示檢測(cè)對(duì)象過(guò)去或當(dāng)前暴露于雌激素或具有雌激素效應(yīng)的環(huán)境干擾物。本研究在弓背青鳉肝臟轉(zhuǎn)錄產(chǎn)物中發(fā)現(xiàn)了多種環(huán)境雌激素的生物標(biāo)志物， 這些基因在性成熟雌性肝臟中的表達(dá)水平明顯高于雄性肝臟， 可用于將弓背青鳉開(kāi)發(fā)成廣鹽性環(huán)境監(jiān)測(cè)模式物種。

3.4 SSR 標(biāo)記

SSRs 是一種重要的分子標(biāo)記， 可作為群體遺傳學(xué)研究的資源。本研究， 從18326 個(gè)unigenes 中共獲得28195 個(gè)SSRs。在已鑒定的SSRs 中， 單核苷酸重復(fù)基序最為豐富， 重復(fù)核苷酸較多的SSRs 較少(表2)， 這與斑石鯛(Oplegnathus punctatus)(Duet al， 2017)、花蟹(Charybdis feriatus)(Zhanget al， 2017)及多鱗鱚(Sillago sihama)(Tianet al， 2019)的SSRs 分析相一致。本研究獲得的弓背青鳉SSR 標(biāo)記， 為遺傳多樣性研究提供了數(shù)據(jù)， 也有助于進(jìn)一步繪制弓背青鳉遺傳連鎖圖譜和開(kāi)發(fā)弓背青鳉遺傳資源。

4 結(jié)論

本研究首次獲得了弓背青鳉雌雄肝臟轉(zhuǎn)錄組。共獲得49912 個(gè)unigenes， 48391 個(gè)unigenes 被定位到主要數(shù)據(jù)庫(kù)， 豐富了弓背青鳉的功能基因資源。通過(guò)功能分析推測(cè)肝臟可能在 弓背青 鳉抵御病原體入侵的免疫防御中發(fā)揮作用， 同時(shí)鑒定到了多個(gè)環(huán)境雌激素生物標(biāo)志基因(vtgs，chgs，er，cyp450等)， 為將其開(kāi)發(fā)成環(huán)境監(jiān)測(cè)物種提供了基礎(chǔ)。此外， 從18326 個(gè)含有SSR 序列的unigenes 中獲得了28195 個(gè)SSR， 為弓背青鳉遺傳多樣性的研究提供數(shù)據(jù)。