陳燕園 ，林志華 ，劉圣，姚韓韓，董迎輝

(1.浙江萬里學院 生物與環(huán)境學院 浙江省水產(chǎn)種質資源高效利用技術研究重點實驗室，浙江 寧波 315100；2.上海海洋大學 水產(chǎn)科學國家級實驗教學示范中心，上海 201306；3.浙江萬里學院 寧海海洋生物種業(yè)研究院，浙江 寧海 315604)

1 引言

糖原是儲存于動物肝臟、肌肉中的支鏈多糖，是一種易于儲存、隨時可用的葡萄糖儲備[1]。糖原的存在既保證了機體的正常運行，又對維持能量動態(tài)平衡起著至關重要的作用。在海洋雙殼貝類中，有些種類的糖原儲存于特定囊泡細胞中，其含量的變化與配子發(fā)育、生長、繁殖、抗逆、風味等密切相關[2]。在滑項薄殼鳥蛤（Fulvia mutica）[3]、魁蚶（Scapharca broughtonii）[4]、縊蟶（Sinonovacula constricta）[5]等多種貝類的研究中發(fā)現(xiàn)，性腺成熟過程中糖原含量顯著下降；在夏季不利的環(huán)境中，糖原含量高的貝類抵抗能力強、死亡率低[6]；在長牡蠣（Crassostrea gigas）[7]、縊蟶[8]中，高糖原含量能顯著增加肉質的鮮味和口感。

糖原磷酸化酶（Glycogen Phosphorylase，GPH）是糖原分解過程中的關鍵酶，催化糖原的磷酸解反應[9]。GPH基因最早在家兔肌肉中發(fā)現(xiàn)，隨后獲得糖原磷酸化酶基因的cDNA 序列[10-11]。在海水貝類中，已克隆獲得長牡蠣和福建牡蠣（Crassostrea angulata）GPH基因的cDNA 全長序列，并發(fā)現(xiàn)GPH基因表達水平與糖原含量變化有關[12-13]。此外，在長牡蠣GPH基因中篩查到了1 個與糖原含量關聯(lián)的SNP 位點[14]。目前，尚未見對縊蟶糖原磷酸化酶基因的相關研究。

縊蟶廣泛分布于我國沿海灘涂，是浙江、福建等省份的主養(yǎng)貝類[15]。不同季節(jié)縊蟶的營養(yǎng)物質、鮮味變化很大，一般春季3-6 月最佳，與其糖原含量密切相關[16]，然而目前有關縊蟶糖原調控的分子機制尚不明晰。本研究克隆了縊蟶糖原磷酸化酶基因（Sc-GPH）cDNA 全長，分析其在縊蟶不同組織和月份的表達水平，并篩選與糖原含量相關的SNP 位點，為高糖原縊蟶的分子育種提供候選基因和標記。

2 材料與方法

2.1 實驗材料

以浙江宏野海產(chǎn)品有限公司養(yǎng)殖的縊蟶“甬樂1 號”為材料，自2019 年2 月開始采集縊蟶樣品（16 月齡，殼長為（61.63±3.69）mm，殼寬為（15.28±2.35）mm，殼高為（19.53±2.26）mm，體重為（14.08±2.11）g），此后每兩個月定時采樣，樣品均為同一批，直至12 月采樣結束。每次取樣，隨機選取12 顆縊蟶，活體解剖取8 個組織（足、外套膜、性腺、閉殼肌、肝胰腺、水管、鰓和唇瓣），經(jīng)液氮速凍后保存于-80℃冰箱，用于糖原含量測定和RNA 提取。

用于SNP 分析的縊蟶群體分別為“甬樂1 號”（YL1）和臺州市野生群體（TZ），各取100 顆個體的足肌，一部分用于RNA 提取，一部分經(jīng)真空冷凍干燥置于-20℃保存，用于糖原含量測定。

2.2 糖原含量分析

取縊蟶不同組織、不同月份及不同群體樣品，按照南京建成生物研究所生產(chǎn)的肝/肌糖原檢測試劑盒（A043）的說明書進行糖原含量的測定。凍干樣品中加入堿液，100°C 加熱20 min 進行皂化。冷卻后，加入蒽酮-硫酸溶液進行反應后，在Spark 多功能酶標儀（TECAN，瑞士）620 nm 波長處測定吸光度。

2.3 Sc-GPH cDNA 全長克隆

根據(jù)縊蟶轉錄組文庫中獲得的Sc-GPH轉錄本序列，設 計3′-RACE 和5′-RACE 特異性引物（表1）。以成貝足中總RNA 為模板，用SMART RACE 5′/3′試劑盒（Clontech，美國）合成cDNA 第一條鏈。以cDNA第一鏈為模板進行3′-RACE和5′-RACE 擴增，其產(chǎn)物用1%瓊脂糖電泳檢測，割膠回收純化后的PCR 產(chǎn)物與pEASY?-T1 克隆載體（全式金，北京）進行連接，接著轉化到DH5α 感受態(tài)細胞（全式金，北京）中，挑選陽性菌落送上海生工公司進行測序。

表1 所用的引物及信息Table 1 The information of primers used in this experiment

2.4 Sc-GPH 序列及系統(tǒng)進化分析

利用DNAMAN 6.6 軟件將所獲得Sc-GPH的5′-RACE 和3′-RACE 片段序列與轉錄本拼接，獲得cDNA 全長；通過ORF Finder 網(wǎng)站（https://www.ncbi.nlm.nih.gov/orffinder/）預測基因的開放閱讀框（OpeanReading Frame,ORF）及氨基酸序列；利用Expasy（https://web.expasy.org/compute_pi/）預測蛋白質的分子量和等電點、Signalp（http://www.cbs.dtu.dk/services/SignalP/）預測信號肽；用Smart（http://smart embl-heidelberg.de/）預測功能域；通過BLAST（http://www.ncbi.nlm.nih.gov/BLAST/）在線工具對Sc-GPH所編碼的氨基酸序列進行相似性分析，并用DNAMAN 6.6 軟件進行多重序列比對；用MEGA 6.0 軟件鄰接法構建系統(tǒng)進化樹。

2.5 實時熒光定量PCR 分析

將縊蟶8 個組織和不同月份縊蟶足、外套膜的cDNA 稀釋50 倍作為模板，用18S rDNA 基因作為內參基因，利用LightCycler?480（Roche，德國）熒光定量PCR儀進行實時定量PCR。反應體系（20μL）：2×SYBR GreenⅠMaster（Roche，德國）10 μL，上、下游引物各1 μL，cDNA 模板2 μL，ddH2O 6 μL；反應程序：95°C 預變性10 min，然后95°C 10 s、60°C 1 min，45 個循環(huán)。為了保證熒光定量PCR 的可靠性，每個樣品做3 個技術重復，每組做4 個生物學重復，Sc-GPH的相對表達量按照2-ΔΔct法計算，并用平均值±標準差表示。數(shù)據(jù)通過SPSS 22.0 軟件進行方差分析，p＜0.05為差異顯著，p＜0.01 為差異極顯著。

2.6 糖原含量相關SNP 的篩選

取縊蟶“甬樂1 號”（YL1）和臺州市野生群體（TZ）2 個群體各100 顆的足肌，提取RNA 并反轉錄合成cDNA。根據(jù)Sc-GPHcDNA 序列，用Primer Premier 5軟件設計引物H-S-F 和H-S-R（表1），PCR 擴增產(chǎn)物進行Sanger 測序。使用MEGA 6.0 軟件進行序列比對，運用SPSS 22.0 軟件將得到SNP 位點的基因型值與基因型值對應的糖原含量進行線性回歸分析，將其中一種等位基因純合子賦值為2，另一種等位基因純合子基因型值賦值為0，雜合子賦值為1[17]。

3 結果

3.1 縊蟶糖原含量的組織差異和月份變化分析

取浙江縊蟶軟體部糖原含量最高的4 月份樣品[17]的不同組織（鰓、閉殼肌、足、外套膜、肝胰腺、性腺、水管、唇瓣），用于糖原含量組織差異分析。結果顯示，足和外套膜糖原含量分別為（103.1±8.3）mg/g 和（111.0±10.5）mg/g，極顯著高于除閉殼肌外的其他組織（p＜0.01，圖1A）?？O蟶足和外套膜周年的糖原含量測定結果表明，上半年的糖原含量顯著高于下半年，從2 月份逐漸增加，4 月份達到峰值，足和外套膜糖原含量分別為（141.3±9.8）mg/g 和（148.3±4.7）mg/g（p＜0.01），隨后逐漸下降至10 月份的最低值，足和外套膜糖原含量分別為（13.0±1.2）mg/g 和（15.3±1.6）mg/g（圖1B）。

圖1 縊蟶不同組織（A）和不同月份足和外套膜（B）的糖原含量變化Fig.1 Changes of glycogen content in different tissues (A) and different months in foot and mantle (B) of Sinonovacula constricta

3.2 Sc-GPH 全長cDNA 序列分析

Sc-GPHcDNA 全長3 963 bp（GenBank 登錄號：MT125681），包括5′非編碼區(qū)（Untranslated Region，UTR）112 bp、3′UTR 1 310 bp、開放閱讀框2 541 bp，共編碼846 個氨基酸。預測其蛋白分子質量（Mw）為97.23 kDa，理論等電點（pI）為6.11；該蛋白質的氨基酸組成中，極性氨基酸所占比例較高，表現(xiàn)為親水性，沒有明顯的信號肽；預測Sc-GPH 蛋白在113～828 aa含有1 個功能域，為Phosphorylase。

選取16 種動物GPH 氨基酸序列與Sc-GPH 氨基酸序列作多重比對。結果發(fā)現(xiàn)，Sc-GPH 與蝦夷扇貝（Mizuhopecten yessoensis，OWF50424.1）、歐洲大扇貝（Pecten maximus，XP_033734012.1）、長牡蠣（Crassostrea gigas，CCN27372.1）GPH 的氨基酸序列相似性較高，分別為78.34%、78.11%和77.87%，與日本對蝦（Penaeus japonicus，BAJ23879.1）、家蠶（Bombyx mori，ACB41088.1）、智人（Homo sapiens，NP_059125.2）等物種的相似性為70.45%～75.86%。系統(tǒng)進化樹結果顯示，縊蟶與歐洲大扇貝、蝦夷扇貝、長牡蠣等貝類親緣關系較近，聚為一支；與昆蟲類、甲殼類、哺乳動物類親緣關系較遠（圖2）。

圖2 縊蟶與其他物種GPH 氨基酸序列系統(tǒng)進化分析Fig.2 Neighbor-joining phylogenetic tree of GPH among Sinonovacula constricta and other species

3.3 不同組織、不同月份Sc-GPH 的表達模式

取8 月份縊蟶的不同組織樣品，利用qRT-PCR 檢測Sc-GPH的組織表達特征，結果顯示，Sc-GPH在8 個組織均有表達，其中足、外套膜的表達量最高（p＜0.01，圖3A）。Sc-GPH在縊蟶不同月份外套膜和足的表達結果表明，Sc-GPH在8 月的表達水平極顯著高于其他月份（p＜0.01，圖3B）。

圖3 Sc-GPH 在縊蟶不同組織（A）和不同月份足和外套膜（B）的的表達特征Fig.3 The expression characteristics of Sc-GPH in different tissues (A) and different months in foot and mantle (B)

3.4 Sc-GPH SNP 位點與糖原含量的關聯(lián)分析

在縊蟶YL1 群體中，Sc-GPH編碼區(qū)域共發(fā)現(xiàn)33 個SNP 位點，其中4 個位點（c.874C＞T、c.930T＞C、c.1133T＞C、c.1284T＞A）與糖原含量的相關性較為顯著，c.1133T＞C、c.1284T＞A 位點差異顯著（p＜0.05，表2）；在TZ 群體同樣篩選到c.930T＞C 位點（表2），且該位點為同義突變。在縊蟶兩個群體中，c.930T＞C 位點與糖原含量相關性均為差異顯著（p＜0.05，表2），CC 基因型個體的糖原含量比TT 基因型的個體高（p＜0.05），說明該位點與縊蟶糖原含量相關。

表2 縊蟶“甬樂1 號”（YL1）和臺州野生群體（TZ）Sc-GPH SNPs 位點與糖原含量的相關性分析Table 2 Association of SNPs of Sc-GPH with glycogen content in YL1 and TZ of Sinonovacula constricta

4 討論

糖原是一種可以直接利用的儲能物質，包括儲存于肌肉中的肌糖原和肝臟中的肝糖原。在哺乳動物中，骨骼肌中的肌糖原占身體糖原總量的2/3，肝臟中的肝糖原占糖原總量的1/3。研究發(fā)現(xiàn)，不同海洋貝類貯存糖原的部位不盡相同，如長牡蠣的糖原主要貯存于性腺、唇瓣中[12]，紫貽貝（Mytilus edulis）的糖原則貯存在外套膜中[18]，而墨西哥灣扇貝（Argopecten irradians concentricus）的糖原貯存在閉殼肌中，并因繁殖周期、環(huán)境而發(fā)生變化[19]。本研究發(fā)現(xiàn)，縊蟶足和外套膜組織的糖原含量高于其他組織，與Yan 等[20]的研究結果基本一致，說明這兩個組織是縊蟶糖原的主要貯存部位，它們具有較高的肌肉含量，肌糖原分解可為其收縮供給充足能量。另外，在縊蟶肝胰腺中亦有較高的糖原含量，推測主要存在于肝細胞中的肝糖原，其分解可以維持血糖濃度的恒定，從而保證機體的正常生理功能。

糖原代謝的多個調節(jié)過程，如葡萄糖的增加和糖原降解的減少，均可能由糖原磷酸化酶引起[21]。研究表明，糖原磷酸化酶激活催化糖原水解為葡萄糖，導致GPH基因表達量增加，從而為機體提供能量[22-23]。此外，GPH基因在不同發(fā)育時期和組織中均有表達，但表達量存在明顯差異，如福建牡蠣GPH基因在性腺和閉殼肌中表達量較高[13]。本研究結果表明，在縊蟶8 個組織中均檢測到Sc-GPH的表達，其中足和外套膜呈現(xiàn)高表達，與其對應的高糖原含量及其存儲糖原能力有關。在許多物種中發(fā)現(xiàn)，肌肉是活性糖原的儲存庫，以便快速提供ATP[24]。2-5 月為縊蟶性腺發(fā)育的休止期，糖原含量較高，與此同時檢測到Sc-GPH在足和外套膜中低表達，說明縊蟶體內正在合成儲存糖原；隨著配子的發(fā)生，Sc-GPH的表達量顯著增加，而糖原含量降低，表明糖原儲存向糖原分解轉變。因此，Sc-GPHmRNA 表達水平的變化與縊蟶繁殖周期關系緊密，其在配子發(fā)生期間可利用足和外套膜中的糖原。

多數(shù)SNP 位于基因組的非編碼區(qū)，在群體遺傳和生物進化研究中起著重要作用[25-26]，而編碼區(qū)的SNP 位點則具有較高的遺傳穩(wěn)定性[27]。本實驗在縊蟶Sc-GPH的編碼區(qū)共檢測到個33 個SNP 位點，SNP 平均密度達到1/72 bp，這與縊蟶生長因子受體結合蛋白2 基因外顯子上 SNP 平均密度（1/65 bp）[28]、長牡蠣糖原磷酸化酶基因外顯子SNP 平均密度（1/35 bp）[29]等存在差異，說明同一物種不同基因或不同物種之間SNP 分布密度不同。在縊蟶Sc-GPH與糖原含量的關聯(lián)分析中發(fā)現(xiàn)，1 個SNP 位點與糖原含量相關，c.930T＞C 位點的CC 基因型個體的糖原含量比TT 基因型的個體高，該位點可作為縊蟶高糖原遺傳育種的重要候選分子標記。