• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    不同水交換量下皺紋盤鮑轉(zhuǎn)錄組測序與分析

    2022-03-01 07:26:06周佳崢郭永軍李永仁黃亞冬
    海洋科學 2022年12期
    關(guān)鍵詞:交換量換水皺紋

    周佳崢, 郭永軍, 全 峰, 梁 健, 李永仁, 黃亞冬, 梁 爽

    不同水交換量下皺紋盤鮑轉(zhuǎn)錄組測序與分析

    周佳崢1, 郭永軍1, 全 峰1, 梁 健1, 李永仁1, 黃亞冬2, 梁 爽1

    (1. 天津農(nóng)學院水產(chǎn)學院, 天津 300380; 2. 天津市海升水產(chǎn)養(yǎng)殖有限公司, 天津 300270)

    養(yǎng)殖過程的循環(huán)水交換量能夠顯著影響皺紋盤鮑()的生理狀態(tài), 進而影響其養(yǎng)殖存活率和生長率, 然而不同循環(huán)水量對于皺紋盤鮑的生理狀態(tài)影響的分子機制尚不明確。本實驗以皺紋盤鮑()為研究對象, 將鮑分別置于循環(huán)水交換頻率為1次/d、2次/d、4次/d的條件下(編號為C0、T2和T4組)暫養(yǎng)3 d, 然后利用Illumina Hiseq 2500 技術(shù)測定不同處理組鮑魚的肌肉組織中轉(zhuǎn)錄組表達情況。結(jié)果表明, 單次換水組(C0)和多次換水組(T2和T4)間差異表達的基因個數(shù)較多, 分別為4 095和3 617個, 而T2與T4組之間的差異表達基因個數(shù)則較少, 僅為232個。對差異表達基因進行GO和KEGG功能分析顯示, 單次換水組(C0)和多次換水組(T2和T4)間差異表達的基因主要集中在能量代謝和免疫反應相關(guān)通路或基因, 隨機選取10個差異基因進行熒光定量PCR驗證, 定量PCR結(jié)果與轉(zhuǎn)錄譜結(jié)果一致, 證實了轉(zhuǎn)錄組測序結(jié)果的可靠性。本研究說明了更高的水交換量下, 鮑的代謝和免疫相關(guān)通路被激活, 可能有助于鮑魚維持較好的生理和免疫狀態(tài), 從而表現(xiàn)出更好的生長、存活等表型。

    皺紋盤鮑; 水交換量; 轉(zhuǎn)錄組; 能量代謝; 免疫

    在水產(chǎn)生物養(yǎng)殖過程中, 水交換量對養(yǎng)殖生物的生理狀態(tài)有顯著影響。通常情況下, 水交換量的增加能夠改善生物的生長和存活等狀態(tài)[1-2]。例如, 在皺紋盤鮑幼鮑的養(yǎng)殖過程中, 當水交換量率由流水、24 h/次、48 h/次和72 h/次逐漸降低時, 幼鮑的存活率和增重率也隨之下降[3]。提高水交換量被認為能夠顯著減少養(yǎng)殖生物在攝食、排泄等過程中產(chǎn)生的氨氮、亞硝酸鹽等自污染因子的濃度, 同時提高溶解氧、透明度等指標, 從而使水生生物維持一個較好的生理狀態(tài), 最終改善其生長、存活情況等[4-5]。例如, 研究換水頻率對刺參養(yǎng)殖水質(zhì)的影響時發(fā)現(xiàn), 日換水率由0、10%、20%、30%、100%逐步提高時, 刺參養(yǎng)殖水體中的氨氮和亞硝酸鹽的濃度逐漸降低[2]; 當換水頻率為10 d/月時, 則能維持優(yōu)良的刺參池塘養(yǎng)殖環(huán)境[6]。

    轉(zhuǎn)錄組分析能夠了解特定條件下生物體內(nèi)不同基因的表達情況, 從而解析生物在生長發(fā)育過程中或外界環(huán)境因子影響下, 其體內(nèi)的生長發(fā)育、抗逆應激等相關(guān)基因的反應情況[7-9], 目前已被廣泛用于研究水產(chǎn)動物的生理功能(如水產(chǎn)動物毒理學、發(fā)育生物學、系統(tǒng)進化生物學等)及免疫功能(包括病毒病、細菌病、寄生蟲病免疫生物學)等研究[10-11]。轉(zhuǎn)錄組技術(shù)在研究水體自污染物質(zhì)(氨氮、亞硝酸鹽等)對水產(chǎn)動物的脅迫效應時也有應用[12]。例如, 氨氮脅迫下的日本沼蝦轉(zhuǎn)錄組結(jié)果顯示, 脅迫前后有887個基因呈差異性表達, 且主要富集在16條免疫相關(guān)通路中, 如補體和凝血級聯(lián)、T細胞受體信號通路、Toll樣受體信號通路與B細胞受體等通路[13]。

    本課題組過去的研究發(fā)現(xiàn), 水交換量的變化能夠顯著影響皺紋盤鮑的成活率、攝食率和耗氧率情況, 具體表現(xiàn)為: 隨著水交換量由1次/d提升到4次/d, 鮑的存活率和攝食量均逐漸提高; 且隨著水交換量的提高, 鮑的耗氧率逐漸降低, 1次換水組的鮑耗氧率(0.47±0.05) mg/(g·h)顯著高于4 次換水組(0.37±0.05) mg/(g·h); 相應地, 1次換水組的水體溶解氧量則小于4次換水組。此外, 水交換量提升降低了鮑養(yǎng)殖水體中的氨氮、亞硝酸氮等自污染因子的濃度[14]。然而, 水交換量變化影響鮑生理狀態(tài)的分子機制尚不明晰。本文擬采用轉(zhuǎn)錄組測序的手段, 研究皺紋盤鮑在不同水交換量的條件下, 其肌肉組織內(nèi)轉(zhuǎn)錄本表達變化, 以闡明不同水交換量對于皺紋盤鮑的生理狀態(tài)影響的分子機制, 為皺紋盤鮑養(yǎng)殖條件的優(yōu)化提供理論支撐, 為貝類健康生態(tài)養(yǎng)殖的理論奠定基礎(chǔ)。

    1 材料與方法

    1.1 材料

    實驗所用的皺紋盤鮑, 引自福建晉江某養(yǎng)殖公司。殼長(2.24±0.30) cm, 殼寬(1.44±0.21) cm, 殼高(1.32±0.19) cm。低溫空運至天津。實驗前進行7 d的暫養(yǎng), 暫養(yǎng)條件為鹽度33、pH 7.8、溫度(24±1) ℃。暫養(yǎng)期間按皺紋盤鮑體質(zhì)量的1%投喂鮑配合飼料。暫養(yǎng)期間每天早晚全量換水, 并及時清除死鮑以及殘餌糞便, 實驗前24 h暫停投喂餌料。

    1.2 方法

    1.2.1 不同水交換量處理

    試由上述馴養(yǎng)的群體中, 隨機取180只皺紋盤鮑, 平均分成3組, 平均分配于9個100 L的PVC養(yǎng)殖箱中, 按換水頻率分為1次每天組、2次每天組、4次每天組(每天全量換水1次、2次、4次, 每組3個平行), 每個換水組實驗持續(xù)3 d, 實驗期間無鮑魚死亡。

    1.2.2 樣品取樣

    換水實驗結(jié)束后, 從每個換水處理組中的PVC整理養(yǎng)殖箱中各取9只鮑(每個平行箱3只), 對其進行肌肉組織的取樣: 使用已滅菌并去除RNA酶的剪刀、鑷子剪取100 mg左右的肌肉組織, 放入干凈的凍存管中, 然后迅速放入液氮中速凍, 之后轉(zhuǎn)移至?80 ℃冰箱保存。

    1.2.3 RNA提取

    對保存的組織進行液氮研磨, 充分研磨后使用Trizol法[15]提取各換水組鮑魚肌肉組織的總RNA。對每個換水組的9個RNA樣品進行三三等量混合, 每組建立3個平行RNA混合樣品。各混樣總RNA利用Aligent 2000、Nanodrop、瓊脂糖凝膠電泳3種方法進行濃度及純度檢測。1次每天組、2次每天組、4次每天組的混合樣本分別標記為C0、T2和T4。

    1.2.4 轉(zhuǎn)錄組測序及分析

    轉(zhuǎn)錄組的建庫及測序工作委托北京百邁克生物科技有限公司完成, 每個換水組的鮑建立3個平行的轉(zhuǎn)錄組測序文庫(共建庫9個)并進行測序Illumina 2000雙端建庫測序。獲得原始測序序列(Sequenced reads)后, 進行以下生物信息學分析: Clean Data的分析; Unigene的拼接和功能注釋; Unigene功能注釋, 包括與NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam數(shù)據(jù)庫的比對進行; Unigene庫基因結(jié)構(gòu)分析主要包括SSR分析、CDS預測和SNP分析; 不同樣品中的基因進行基于FPKM的表達量分析; 基于基因在不同樣品中的表達量, DEseq識別鑒定差異表達基因; 對差異表達基因進行模式聚類、GO功能注釋以及KEGG富集分析。

    1.2.5 熒光定量PCR驗證

    從轉(zhuǎn)錄組結(jié)果里隨機選取10個差異基因進行qRT-PCR驗證, 根據(jù)所挑選的差異表達基因的序列, 遵循引物設(shè)計的原則, 使用引物設(shè)計軟件Beacon designer 7設(shè)計兩對定量引物。經(jīng)標準曲線法測定引物的擴增效率, 保證各定量引物的擴增效率類似, 且在95%~105%范圍內(nèi)。以上述反轉(zhuǎn)錄合成的cDNA為模板, β-actin基因為內(nèi)參[16], 在熒光定量PCR儀上進行定量PCR反應。反應體系包括: 1 μL的cDNA模板(10倍稀釋), 各1 μL的上下游引物(10pmol L-1), 10 μL的2×DNAmo ColorFlash Master Mix(Thermo, USA), 以及8 μL的滅菌超純水。反應程序設(shè)置為: 95 ℃ 7 min, (95 ℃ 20 s, 60 ℃ 60 s)×35, 72 ℃ 10 min, 每個反應后均檢查熔解曲線以確保產(chǎn)物單一性, 每個反應均設(shè)置3個重復。

    2 結(jié)果與分析

    對9個樣品進行Illumina 2000雙端建庫測序, 共得到207 435 516條原始測序序列, 測序正確率高達99.9%。對原始測序序列進行質(zhì)量篩選, 即截除Reads中的測序接頭、引物序列和低質(zhì)量值數(shù)據(jù)等, 共獲得61.81 Gb Clean data, 其中各樣品GC含量在46.60%~ 47.76%之間, Q30堿基百分比不小于94.39%, 表明測序數(shù)據(jù)完整性好, 質(zhì)量高(表1)。

    表1 樣品測序數(shù)據(jù)質(zhì)量評估

    將得到的所有Unigene序列分別與七個數(shù)據(jù)庫對比, 比對結(jié)果如表2所示, 可知在COG、GO、KEGG、KOG、Pfam、Swissprot、eggNOG七個數(shù)據(jù)庫中成功注釋的Unigene數(shù)量分別為6 518、7 750、9 116、12 139、15 772、11 014、16 863, 最終總共獲得20 480個至少在一個數(shù)據(jù)庫中有注釋信息的Unigene。其中長度大于等于1 000的Unigene數(shù)為16 025個, 長度介于300至1 000之間的Unigene數(shù)為4 455個。

    表2 Unigene注釋統(tǒng)計

    利用皮爾遜相關(guān)系數(shù)(Pearson’s correlation coefficient)對樣品間相關(guān)性的指標評估。如圖1所示, C0、T2、T4各組的生物學重復文庫之間的相關(guān)性均較高(2均大于0.98), 而C0和T2, C0和T4的樣品文庫之間的相關(guān)性則較低(2均小于0.25)。

    圖1 樣品的相關(guān)性熱圖

    如表3所示, C0和T2、C0和T4、T4和T2的差異表達基因數(shù)目分別為4 095、3 617、232, 且1次換水組與多次換水組間(C0和T2和C0和T4)的差異表達基因中, 上調(diào)基因的數(shù)目明顯小于下調(diào)基因的數(shù)目。

    表3 差異表達基因數(shù)目

    分別將單次換水組(C0)和多次換水組(T2和T4)的差異表達基因進行GO功能注釋, 按照GO三大類型再次進行分類, 共分為52個功能亞類, 如圖2所示, C0和T2和C0和T4的差異表達基因所注釋到的GO類群比較相似: 在生物學過程分類下, 主要富集在細胞學過程(cellular process)、定位(localization)、代謝過程(metabolic process)和單生物進程(single-organism process)等功能; 在細胞組分分類下, 主要富集在細胞(cell)、細胞部分(cell part)、細胞膜部分(membrane part)和細胞膜部分(membrane)等功能; 在分子功能中, 主要富集在催化活性(catalytic activi-ty)、附著(bin-ding)和運輸活性(transporter activity)等功能。

    分別將單次換水組(C0)和多次換水組(T2和T4)的差異表達基因進行KEGG功能注釋, 以此獲知不同換水組之間的差異基因可能參與的生理生化代謝通路。結(jié)果如圖3所示, C0和T2的差異表達基因主要顯著富集于淀粉和蔗糖代謝(starch and sucrose meta-bolism), 溶酶體(lysosome), 糖胺聚糖降解(glycosa-minoglycan degradation), 戊糖和葡萄糖醛酸酯的相互轉(zhuǎn)化(pentose and glucuronate interconversions), 氨基酸和核苷酸糖的代謝(amino sugar and nucleotide sugar metabolism)等通路, C0和T4的差異表達基因主要顯著富集于淀粉和蔗糖代謝(starch and sucrose metabolism), 溶酶體(lysosome), 糖胺聚糖降解(gly-cosa-mi-no-glycan degradation), 半乳糖代謝(galactose meta-bolism), 氨基酸和核苷酸糖代謝(amino sugar and nucleotide sugar metabolism)等通路。

    根據(jù)FPKM(表達量)值, 選取不同水交換量處理組之間差異表達最顯著的基因, 對其表達量及NR注釋信息進行分析、篩選和統(tǒng)計, 結(jié)果如表4和表5所示, 在C0和T2差異表達最顯著基因中, 與生物消化代謝(胰蛋白酶trypsin等)與組織構(gòu)成(透明質(zhì)蛋白perlucin, 黏液相關(guān)蛋白mucus-associated protein, 肌球蛋白myosin heavy chain等)等功能相關(guān)的基因顯著在多次換水組T2和T4中上調(diào)表達, 而在1次換水組C0中表達量較低; 而透明質(zhì)素樣蛋白(perlucin-like protein)海藻酸裂解酶(alginate lyase)內(nèi)切葡聚糖酶4樣亞型(endoglucanase 4-like isoform X1)等基因則在多次換水組中表達量較低(表4)。C0和T4差異表達最顯著基因中, 與組織構(gòu)成(睪蛋白testisin等), 生物消化代謝(胰蛋白酶trypsin等), 信號轉(zhuǎn)導(wnt信號通路受體frizzled 1, RyR受體ryanodine receptor-like等)在多次換水組T2和T4中上調(diào)表達, 而在1次換水組C0中表達量較低; 而透明質(zhì)蛋白(perlucin), 纖維素酶(cellulase), 黏蛋白(mucin-5AC-like isoform X3), 內(nèi)切葡聚糖酶(endo-glucanase 14-like)等基因則在多次換水組中表達量較低(表5)。

    從轉(zhuǎn)錄組數(shù)據(jù)庫隨機選取10個差異基因進行熒光定量PCR驗證。運用2–ΔΔCt方法對不同水交換量組中的差異表達基因進行相對表達水平定量分析。結(jié)果表示, 選取的10個差異表達基因(Unigene編號c100088, c100118, c100155, c100277, c100299, c100496, c100761, c100829, c100882)在三個換水組鮑中的轉(zhuǎn)錄本表達量(以FPKM的對數(shù)表示)和熒光定量PCR表達量(以倍數(shù)變化FC的對數(shù)表示)一致(圖4), 表明轉(zhuǎn)錄組測序結(jié)果是可靠的。

    3 討論

    隨著新一代高通分子測序技術(shù)的不斷完善和發(fā)展與成熟, 轉(zhuǎn)錄組研究在水產(chǎn)生物生理過程機制的研究中愈發(fā)常見, 已有大量研究將轉(zhuǎn)錄組學技術(shù)應用于水產(chǎn)生物疾病與免疫學研究[17-18], 繁殖和發(fā)育研究[19-20], 毒性和壓力過程研究[21-22], 生長和營養(yǎng)[23-24]等方面。本研究利用高通量測序平臺, 對在三種不同水交換量條件下養(yǎng)殖的皺紋盤鮑的肌肉進行轉(zhuǎn)錄組測序, 得到了質(zhì)量良好、結(jié)果可信的轉(zhuǎn)錄組測序結(jié)果數(shù)據(jù)。比較C0和T2、C0和T4與T2和T4數(shù)據(jù), 分別檢測到4 095、3 617、232個差異表達基因, 說明隨著水交換量的增加, 尤其是由1次/d的水交換量提升至2次/d或4次/d時, 顯著改變鮑中相關(guān)基因的表達水平, 說明換水量的增加對鮑的生理生化過程產(chǎn)生了比較明顯的影響, 也印證了使用轉(zhuǎn)錄組的手段研究水交換量對于鮑的生理狀態(tài)的影響機制是可行的。

    圖3 差異表達基因京都基因與基因組數(shù)據(jù)庫(KEGG)通路富集散點圖

    Fig. 3 Scatter diagram of the differentially expressed genes enriched in the KEGG pathways

    表4 1次換水組和2次換水組差異表達最顯著基因功能和表達量

    表5 1次換水組和4次換水組差異表達最顯著基因功能和表達量

    圖4 熒光定量聚合酶鏈式反應(PCR)驗證結(jié)果

    注: FC: 定量PCR差異倍數(shù); FPKM: 轉(zhuǎn)錄組表達量

    本研究中, 通過對C0組和T2、T4組的差異表達基因進行GO和KEGG富集分析, 發(fā)現(xiàn)差異表達基因主要富集到了與消化和代謝相關(guān)的基因通路上, 如淀粉和蔗糖代謝, 糖胺聚糖降解, 戊糖和葡萄糖醛酸酯的相互轉(zhuǎn)化, 氨基酸和核苷酸糖的代謝, 半乳糖代謝等通路。在對不同換水組之間的差異表達最顯著的基因進行功能分析時, 同樣發(fā)現(xiàn)了許多基因的功能屬于消化、代謝類, 如胰蛋白酶、內(nèi)切葡聚糖酶、纖維素酶等。結(jié)合我們過去的研究, 在投喂顆粒狀飼料的鮑養(yǎng)殖水體中, 1次換水組、2次換水組、3次換水組、4次換水組的水體中的氨氮含量分別為(1.46±0.42) mg/L、(0.93±0.59) mg/L、(0.46± 0.11) mg/L、(0.39±0.26) mg/L[14]。說明了換水量的提高很有可能通過改善水質(zhì), 提升了鮑魚機體對于營養(yǎng)物質(zhì)的消化和吸收效率。類似地, 日本沼蝦在應對不同濃度的氨氮脅迫時, 其糖酵解相關(guān)通路、蛋白質(zhì)消化和吸收相關(guān)通路也會發(fā)生相應的變化[13], 說明水生生物在水質(zhì)變化時, 其消化代謝相關(guān)生理過程的活性會發(fā)生變化, 從而從能量供給的角度幫助機體應對由水質(zhì)變化帶來的不利影響。

    此外, 我們發(fā)現(xiàn)溶酶體相關(guān)通路和基因在不同換水組鮑中有顯著的變化。溶酶體是無脊椎動物免疫系統(tǒng)中一個關(guān)鍵性的酶復合體, 由多種抗體、受體和酶成分組成, 包括溶酶體酶、凝集素、抗菌肽等體液免疫因子, 它以殺菌、促進吞噬等方式參與貝類的免疫防御[25-26]。有研究表明, 皺紋盤鮑體內(nèi)的溶酶體凝集素能夠參與對外源病原的識別和消滅過程[27]。雜色鮑的溶酶體酚氧化酶活性在注射副溶血弧菌菌液24 h后顯著增加, 且顯著受到其他外部環(huán)境因子如溫度和溶解氧的影響[28-29]。結(jié)合本實驗的結(jié)果, 說明了不同水交換量不僅能夠影響皺紋盤鮑的消化代謝功能, 而且可能通過影響其溶酶體相關(guān)的免疫系統(tǒng)功能, 從而影響皺紋盤鮑整體的免疫生理狀態(tài), 最終使得鮑表現(xiàn)出更好的生理狀態(tài)。

    4 結(jié)論

    本文研究了不同循環(huán)水交換量對皺紋盤鮑()轉(zhuǎn)錄組表達的影響, 研究結(jié)果證實了不同水交換頻率下, 鮑魚的肌肉轉(zhuǎn)錄本呈現(xiàn)顯著的差異性表達。這些差異表達的基因主要功能為能量代謝或免疫反應相關(guān), 預示更高的水交換量可能對皺紋盤鮑代謝和免疫相關(guān)通路的激活有利, 未來應進一步聚焦代表性差異表達基因的具體生物學功能。

    [1] 馬貴范, 郭文學, 王玲玲, 等. 養(yǎng)殖密度和換水量及頻率對凡納濱對蝦生長的影響[J]. 漁業(yè)現(xiàn)代化, 2018, 45(4): 21-26.

    MA Guifan, GUO Wenxue, WANG Lingling, et al. The effects of aquaculture density and water exchange rate and frequency on the growth of[J]. Fishery Modernization, 2018, 45(4): 21-26.

    [2] 王光玉, 姜佳惠, 祝楠, 等. 換水率和密度對刺參生長和水質(zhì)的影響[J]. 漁業(yè)現(xiàn)代化, 2017, 44(6): 55-61.

    WANG Guangyu, JIANG Jiahui, ZHU Nan, et al. Effects of water exchange rate and density on the growth and water quality of sea cucumbers[J]. Fishery Moder-nization, 2017, 44(6): 55-61.

    [3] 楊鳳, 雷衍之, 王仁波, 等. 皺紋盤鮑自污染及其對幼鮑生長及成活率的影響[J]. 大連水產(chǎn)學院學報, 2003, 18(1): 1-6.

    YANG Feng, LEI Yanzhi, WANG Renbo, et al. Self- contamination ofand its effect on growth and survival rate of juvenile abalone[J]. Journal of Dalian Fisheries University, 2003, 18(1): 1-6.

    [4] 唐黎標. 水質(zhì)對水產(chǎn)養(yǎng)殖的影響和對策[J]. 漁業(yè)致富指南, 2017, 18: 36-37.

    TANG Libiao. The impact of water quality on aquaculture and its countermeasures[J]. Guide to Getting Rich in Fisheries, 2017, 18: 36-37.

    [5] 張偉. 水產(chǎn)養(yǎng)殖系統(tǒng)的水質(zhì)預測預警的應用研究[D]. 大連: 大連海洋大學, 2019.

    ZHANG Wei. Application research on water quality prediction and early warning of aquaculture system[D]. Dalian: Dalian Ocean University, 2019.

    [6] 盧超超, 孫永軍, 高勤峰, 等. 不同換水頻率對刺參養(yǎng)殖池塘底泥環(huán)境的影響[J]. 海洋湖沼通報, 2014(4): 59-66.

    LU Chaochao, SUN Yongjun, GAO Qinfeng, et al. The impact of different water exchange frequencies on the sediment environment of sea cucumber culture ponds[J]. Transactions of Oceanology and Limnology, 2014(4): 59-66.

    [7] CINGOLANI P. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome ofstrain w1118; iso-2; iso-3[J]. Fly, 2012, 6(2): 1-13.

    [8] HUSSING C, KAMPMANN M L, MOGENSEN H S, et al. Comparison of techniques for quantification of next- generation sequencing libraries[J]. Forensic ence Inter-national Genetics Supplement Series. 2015, 5(4): e276- e278.

    [9] 王旭. 基于比較基因組學和轉(zhuǎn)錄組學技術(shù)揭示單增李斯特菌毒力因子的研究[D]. 上海: 上海海洋大學, 2016.

    WANG Xu. Research on the virulence factors ofmonocytogenes revealed based on comparative genomics and transcriptomics techniques[D]. Shanghai: Shanghai Ocean University, 2016.

    [10] 趙婷. 青蛤()轉(zhuǎn)錄組文庫測序及免疫相關(guān)序列分析[D]. 天津: 天津師范大學, 2014.

    ZHAO Ting. Sequencing and immune-related sequence analysis oftranscriptome library[D]. Tianjin: Tianjin Normal University, 2014.

    [11] 黃玉婷. 基于轉(zhuǎn)錄組數(shù)據(jù)的長牡蠣未折疊蛋白反應(UPR)熱脅迫響應機制初步研究[D]. 大連: 大連海洋大學, 2018.

    HUANG Yuting. Preliminary study on the heat stress response mechanism of the unfolded protein response (UPR) of Pacific oyster based on transcriptome data[D]. Dalian: Dalian Ocean University, 2018.

    [12] 章瓊, 蔣高中, 李冰. 水產(chǎn)動物對氨氮脅迫響應的轉(zhuǎn)錄組分析研究進展[J]. 江蘇農(nóng)業(yè)科學, 2015, 43(3): 227-230.

    ZHANG Qiong, JIANG Gaozhong, LI Bing. Transcri-ptome analysis research progress of aquatic animals’ response to ammonia nitrogen stress[J]. Jiangsu Agricultural Sciences, 2015, 43(3): 227-230.

    [13] 孫健. 日本沼蝦抗氨氮脅迫轉(zhuǎn)錄組與代謝組分析[D]. 泰安: 山東農(nóng)業(yè)大學, 2017.

    SUN Jian. Analysis of transcriptome and metabolome ofresistance to ammonia nitrogen stress[D]. Taian: Shandong Agricultural University, 2017.

    [14] 陳俊男, 郭永軍, 梁健, 等. 不同養(yǎng)殖方式對皺紋盤鮑養(yǎng)殖效果及養(yǎng)殖水質(zhì)的影響[J]. 海洋科學, 2020, 44(3): 104-112.

    CHEN Junnan, GUO Yongjun, LIANG Jian, et al. The effects of different culture methods on the aquaculture effect and the quality of aquaculture water of[J]. Marine Sciences, 2020, 44(3): 104-112.

    [15] HUMMON A B, LIM S R, DIFILIPPANTONIO M J, et al. Isolation and solubilization of proteins after TRIzol extraction of RNA and DNA from patient material following prolonged storage[J]. Biotechniques, 2007, 42(4): 472.

    [16] LI J, SUN H, LIU X. Acclimation-dependent expression of heat shock protein 70 in Pacific abalone (Ino) and its acute response to thermal exposure[J]. Chinese Journal of Oceanology & Limnology, 2012, 30(1): 146-151.

    [17] HUANG L X, LIU W J, JIANG Q L, et al. Integration of transcriptomic and proteomic approaches reveals the temperature-dependent virulence of[J]. Frontiers in Cellular and Infection Microbiology, 2018, 8: 207.

    [18] YE H, LIN Q, LUO H. Applications of transcriptomics and proteomics in understanding fish immunity[J]. Fish & Shellfish Immunology, 2018, 77: 319-327.

    [19] CHAPMAN R W, READING B J, SULLIVAN C V, et al. Ovary Transcriptome profiling via artificial intelligence reveals a transcriptomic fingerprint predicting egg quality in striped bass,[J]. PLoS One, 2014, 9(5): e96818.

    [20] AVARRE J C, DUGUé R, ALONSO P, et al. Analysis of the black-chinned tilapiareproducing under a wide range of salinities: from RNA-seq to candidate genes.[J]. Molecular Ecology Resources, 2013, 14(1): 139-149.

    [21] COLLI-DULA R C, FANG X, MORAGA-AMADOR D, et al. Transcriptome analysis reveals novel insights into the response of low-dose benzo(a)pyrene exposure in male tilapia[J]. Aquatic Toxicology, 2018, 201(3/4): 162-173.

    [22] CHUEH T C, HSU L S, KAO C M, et al. Transcriptome analysis of zebrafish embryos exposed to deltamethrin[J]. Environmental Toxicology, 2016, 32(5): 1548-1557.

    [23] SUN S, XUAN F, GE X, et al. Identification of differentially expressed genes in hepatopancreas of oriental river prawn,exposed to environmental hypoxia[J]. Gene, 2014, 534(2): 298- 306.

    [24] DIEGO R, CARLOS F, MIGUEL H, et al. Integrative transcriptome, genome and quantitative trait loci resources identify single nucleotide polymorphisms in candidate genes for growth traits in turbot[J]. International journal of molecular sciences, 2016, 17(2): 243.

    [25] CHANDHINI S, REJISH KUMAR V J. Transcripto-mics in aquaculture: current status and applications[J]. Reviews in Aquaculture, 2019, 11(4): 1379-1397.

    [26] 杜麗, 張巍, 陸逵, 等. 貝類免疫機制研究進展[J]. 動物醫(yī)學進展, 2008(3): 77-81.

    DU Li, ZHANG Wei, LU Kui, et al. Research progress in shellfish immune mechanism[J]. Advances in Veterinary Medicine, 2008(3): 77-81.

    [27] WANG N, WHANG I, LEE J. A novel C-type lectin from abalone,, agglutinates[J]. Developmental & Comparative Immunology, 2008, 32(9): 1034-1040.

    [28] CHENG W, HSIAO I S, HSU C H, et al. Change in water temperature on the immune response of Taiwan abalonesupertexta and its susceptibility to[J]. Fish & Shellfish Immunology, 2004, 77: 319-327.

    [29] CHENG W, LI C H, CHEN J C. Effect of dissolved oxygen on the immune response ofsupertexta and its susceptibility to[J]. Aquaculture, 2004, 232(14): 103-115.

    Sequencing and analysis of thetranscriptome under different water exchange rates

    ZHOU Jia-zheng1, GUO Yong-jun1, QUAN Feng1, LIANG Jian1, LI Yong-ren1, HUANG Ya-dong2, LIANG Shuang1

    (1. Fisheries College, Tianjin University of Agriculture, Tianjin 300380, China; 2. Tianjin Haisheng Aqua-cu-l-ture Limited Company, Tianjin 300270, China)

    The rate of circulating water exchange can remarkably affect the physiology of abalone, including survival and growth. However, the molecular mechanism of this effect is unclear. In this study, abalones were placed under a circulating water exchange frequency of 1, 2, or 4 times/day (called C0, T2, and T4, respectively) for 3 days. Illumina Hiseq 2500 technology was used to investigate the transcriptome expression in the muscle tissues of abalone from the different water exchange groups. The results showed that many differentially expressed genes were detected between the single water exchange group (C0) and the multiple water exchange groups (T2 and T4), such as 4, 095 and 3, 617, respectively, while the number of differentially expressed genes between the T2 and T4 groups was relatively small at 232. The GO and KEGG functional analyses showed that the differentially expressed genes between the single water exchange group (C0) and the multiple water exchange group (T2 and T4) were mainly classified into energy metabolism and immune response-related pathways or genes. Ten differentially expressed genes were selected for verification by fluorescence quantitative polymerase chain reaction, and the results were consistent with the transcription profile results, confirming the reliability of the transcriptome sequencing results. This study shows that metabolism and the immune-related pathways of abalone were activated under the higher water exchange rate, which may help abalone maintain better physiological and immune status, thus leading to better growth, survival, and other phenotypes.

    ; water exchange rate; transcriptome; energy metabolism; immunity

    Oct. 3, 2020

    S966.9

    A

    1000-3096(2022)12-0159-10

    10.11759/hykx20201003001

    2020-10-03;

    2020-11-09

    天津市種業(yè)科技重大專項(19ZXZYSN00080); 天津市自然科學基金(18JCQNJC84500); 廣東省水產(chǎn)健康安全養(yǎng)殖重點實驗室開放基金(GDKLHSA0806); 2020年天津市大學生創(chuàng)新創(chuàng)業(yè)項目(202010061010); 2020年現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項資金(CARS49)

    [Tianjin Seed Industry Science and Technology Major Project No. 19ZXZYSN00080; Tianjin Natural Science Foundation, No. 18JCQNJC84500; Guangdong Key Laboratory of Aquatic Health and Safety Aquaculture Open Fund, No. GDKLHSA0806; Tianjin University Student Innovation and Entrepreneurship Project in 2020, No. 202010061010; Modern Special Fund for Construction of Agricultural Industrial Technology System, No. CARS49]

    周佳崢(1999—), 男, 本科生, 水產(chǎn)養(yǎng)殖學專業(yè), E-mail: 3055912470@qq.com; 梁爽(1989—), 男,通信作者, 博士, 講師, 主要從事貝類養(yǎng)殖生理生態(tài)研究, E-mail: shuangliang@tjau.edu.cn

    (本文編輯: 楊 悅)

    猜你喜歡
    交換量換水皺紋
    非屏蔽門地鐵車站隧道與站臺空氣交換量模擬研究
    杭州亞運泳池不換水的“黑科技”
    軍事文摘(2023年24期)2023-12-19 06:49:50
    清涼峰自然保護區(qū)土壤陽離子交換量的剖面分布特征及其影響因素1)
    水培植物霎換水
    母親臉上的皺紋
    自動換水
    皺紋
    黑龍江西部地區(qū)土壤pH值、陽離子交換量和有機質(zhì)的分布特征
    防護林科技(2020年4期)2020-06-17 03:02:36
    淺談土壤陽離子交換量測試方法
    智富時代(2019年2期)2019-04-18 07:44:42
    身高線與皺紋線
    最后的刺客免费高清国语| 麻豆久久精品国产亚洲av| 国产亚洲精品久久久com| 亚州av有码| 日本午夜av视频| 久久鲁丝午夜福利片| 日韩一区二区三区影片| 成年女人永久免费观看视频| 欧美精品一区二区大全| 国产大屁股一区二区在线视频| 一二三四中文在线观看免费高清| 国产三级在线视频| 免费观看在线日韩| 黄色配什么色好看| 极品教师在线视频| 日日摸夜夜添夜夜爱| 午夜a级毛片| 精品不卡国产一区二区三区| 韩国高清视频一区二区三区| 国产成人精品婷婷| 亚洲电影在线观看av| 亚洲人成网站在线播| 精品免费久久久久久久清纯| 亚洲中文字幕一区二区三区有码在线看| 真实男女啪啪啪动态图| 日韩亚洲欧美综合| 有码 亚洲区| 国产精品国产三级专区第一集| 可以在线观看毛片的网站| 成人鲁丝片一二三区免费| 欧美成人a在线观看| 亚洲av男天堂| 日本黄色视频三级网站网址| 国产探花在线观看一区二区| 美女高潮的动态| 一个人免费在线观看电影| 午夜a级毛片| 亚洲av二区三区四区| 99久久成人亚洲精品观看| 99久国产av精品国产电影| 国产一区二区在线av高清观看| 久久99精品国语久久久| 少妇高潮的动态图| 亚洲人成网站高清观看| 亚洲av免费高清在线观看| 欧美性猛交黑人性爽| 亚洲成人精品中文字幕电影| 一个人看视频在线观看www免费| 又爽又黄a免费视频| 性色avwww在线观看| 亚洲18禁久久av| 久久久久久九九精品二区国产| 麻豆成人av视频| 晚上一个人看的免费电影| 国产黄片视频在线免费观看| 精华霜和精华液先用哪个| 97超碰精品成人国产| 麻豆av噜噜一区二区三区| 一夜夜www| 久久久久久久久大av| 免费看a级黄色片| 熟妇人妻久久中文字幕3abv| 中文字幕免费在线视频6| 精品欧美国产一区二区三| 综合色av麻豆| 热99在线观看视频| 日本免费一区二区三区高清不卡| 特级一级黄色大片| av黄色大香蕉| 国产中年淑女户外野战色| 久久精品久久精品一区二区三区| 99久久人妻综合| 狂野欧美激情性xxxx在线观看| 国产伦精品一区二区三区视频9| 久久热精品热| 最近中文字幕2019免费版| 国产成人freesex在线| 精品人妻一区二区三区麻豆| 国产真实乱freesex| 亚洲av一区综合| 丰满人妻一区二区三区视频av| 国产av在哪里看| 亚洲欧美日韩东京热| 99久久人妻综合| 韩国av在线不卡| 亚洲熟妇中文字幕五十中出| 尾随美女入室| 亚洲久久久久久中文字幕| 99久久精品国产国产毛片| 国产精品精品国产色婷婷| av又黄又爽大尺度在线免费看 | 欧美不卡视频在线免费观看| 日韩视频在线欧美| 嘟嘟电影网在线观看| 亚洲精品乱码久久久v下载方式| 免费黄色在线免费观看| 91精品伊人久久大香线蕉| 国产一区二区亚洲精品在线观看| 蜜臀久久99精品久久宅男| 国产伦精品一区二区三区四那| 一本一本综合久久| 免费观看人在逋| 99在线人妻在线中文字幕| 婷婷色综合大香蕉| 成人毛片a级毛片在线播放| 成人高潮视频无遮挡免费网站| 青春草视频在线免费观看| 干丝袜人妻中文字幕| 国内精品宾馆在线| 毛片一级片免费看久久久久| 国产精品av视频在线免费观看| 男女啪啪激烈高潮av片| 你懂的网址亚洲精品在线观看 | 一边亲一边摸免费视频| av免费观看日本| 91久久精品电影网| www.av在线官网国产| 97超碰精品成人国产| 亚洲天堂国产精品一区在线| 久久99热6这里只有精品| 国产成人a∨麻豆精品| 国产亚洲午夜精品一区二区久久 | 久久亚洲精品不卡| 黄色日韩在线| 99热这里只有是精品在线观看| 最近中文字幕2019免费版| 永久网站在线| eeuss影院久久| 国产精品美女特级片免费视频播放器| 日韩强制内射视频| 美女国产视频在线观看| 大话2 男鬼变身卡| 少妇人妻精品综合一区二区| 看非洲黑人一级黄片| 久久这里有精品视频免费| 亚洲欧美日韩卡通动漫| 男女啪啪激烈高潮av片| 成人二区视频| 久久亚洲精品不卡| 国产午夜精品一二区理论片| 国产精品久久视频播放| 亚洲精品日韩在线中文字幕| 又爽又黄无遮挡网站| 久久99蜜桃精品久久| 国产 一区 欧美 日韩| 精品午夜福利在线看| 免费一级毛片在线播放高清视频| 国产高清国产精品国产三级 | 久久欧美精品欧美久久欧美| 国产精品乱码一区二三区的特点| 精品国内亚洲2022精品成人| 别揉我奶头 嗯啊视频| 精品无人区乱码1区二区| 国产精品美女特级片免费视频播放器| 一本一本综合久久| 一本久久精品| 能在线免费看毛片的网站| 免费搜索国产男女视频| 日韩av在线大香蕉| 男女边吃奶边做爰视频| 高清日韩中文字幕在线| 久久精品久久久久久噜噜老黄 | 欧美zozozo另类| 乱人视频在线观看| 成年女人永久免费观看视频| 高清av免费在线| 日韩三级伦理在线观看| 七月丁香在线播放| 久久精品夜夜夜夜夜久久蜜豆| www日本黄色视频网| 久久久午夜欧美精品| 午夜福利在线在线| 精品一区二区免费观看| 丰满人妻一区二区三区视频av| 国产高清视频在线观看网站| 亚州av有码| 国产精品一及| www.色视频.com| 草草在线视频免费看| 村上凉子中文字幕在线| 亚洲欧美精品专区久久| 国产精品无大码| 日韩欧美精品v在线| 国产真实伦视频高清在线观看| 日本午夜av视频| 美女大奶头视频| 日韩中字成人| 久久久久精品久久久久真实原创| 久久国产乱子免费精品| 国产一区二区在线观看日韩| 精品一区二区免费观看| 国产私拍福利视频在线观看| 久久久久久伊人网av| 亚洲综合色惰| 中文精品一卡2卡3卡4更新| 国产精品一区二区三区四区久久| av又黄又爽大尺度在线免费看 | 欧美97在线视频| 又黄又爽又刺激的免费视频.| 日本三级黄在线观看| 国产精品一区www在线观看| 六月丁香七月| 两个人视频免费观看高清| 偷拍熟女少妇极品色| 国产精品日韩av在线免费观看| 爱豆传媒免费全集在线观看| 亚洲欧美成人综合另类久久久 | 久久国内精品自在自线图片| 春色校园在线视频观看| 插阴视频在线观看视频| 男人的好看免费观看在线视频| 一边摸一边抽搐一进一小说| 国产国拍精品亚洲av在线观看| 国产一区亚洲一区在线观看| 国产国拍精品亚洲av在线观看| av又黄又爽大尺度在线免费看 | 亚洲真实伦在线观看| 午夜福利网站1000一区二区三区| 夜夜看夜夜爽夜夜摸| 国产免费一级a男人的天堂| 午夜福利网站1000一区二区三区| 久久久亚洲精品成人影院| 日韩av在线免费看完整版不卡| 男人舔奶头视频| 最新中文字幕久久久久| 五月玫瑰六月丁香| 欧美性感艳星| 成人无遮挡网站| 熟女人妻精品中文字幕| 国产单亲对白刺激| 国产综合懂色| 赤兔流量卡办理| 白带黄色成豆腐渣| 国产高清不卡午夜福利| 欧美一区二区精品小视频在线| 国产精品国产三级专区第一集| 一本一本综合久久| 亚洲四区av| 热99re8久久精品国产| 99久国产av精品| 国产 一区 欧美 日韩| 91精品一卡2卡3卡4卡| 国语自产精品视频在线第100页| 精品一区二区三区人妻视频| 日本与韩国留学比较| 色综合亚洲欧美另类图片| 国产视频首页在线观看| av在线播放精品| 国产精品熟女久久久久浪| 久久精品国产亚洲网站| 超碰97精品在线观看| 老司机福利观看| 国产 一区 欧美 日韩| kizo精华| 亚洲av二区三区四区| 国内精品美女久久久久久| ponron亚洲| 国产三级在线视频| 免费av毛片视频| 欧美一区二区国产精品久久精品| 国产黄片视频在线免费观看| 午夜免费激情av| 成人鲁丝片一二三区免费| 欧美高清性xxxxhd video| av专区在线播放| 亚洲精华国产精华液的使用体验| 天堂av国产一区二区熟女人妻| 亚洲精品亚洲一区二区| 国产精品一及| 色综合站精品国产| a级毛片免费高清观看在线播放| 国产黄色小视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲国产精品成人综合色| 亚洲国产最新在线播放| 亚洲欧美日韩高清专用| 国产片特级美女逼逼视频| 亚洲国产色片| 欧美成人免费av一区二区三区| 观看免费一级毛片| 国产日韩欧美在线精品| 国产免费又黄又爽又色| 国产大屁股一区二区在线视频| 亚洲三级黄色毛片| 国内精品一区二区在线观看| 网址你懂的国产日韩在线| 91久久精品电影网| 免费观看在线日韩| 亚洲av二区三区四区| 在线天堂最新版资源| 国产极品天堂在线| 七月丁香在线播放| 两个人视频免费观看高清| 男人的好看免费观看在线视频| 麻豆精品久久久久久蜜桃| 十八禁国产超污无遮挡网站| 熟妇人妻久久中文字幕3abv| 国产免费一级a男人的天堂| 亚洲av.av天堂| 国产淫片久久久久久久久| 亚洲av中文字字幕乱码综合| 亚洲欧美日韩无卡精品| 一边摸一边抽搐一进一小说| 观看免费一级毛片| 亚洲精品亚洲一区二区| 嫩草影院新地址| 久久精品国产亚洲av天美| 国产女主播在线喷水免费视频网站 | 精品久久久久久久久久久久久| 赤兔流量卡办理| 欧美区成人在线视频| av免费观看日本| 噜噜噜噜噜久久久久久91| 国语对白做爰xxxⅹ性视频网站| 欧美3d第一页| 国产亚洲一区二区精品| 最后的刺客免费高清国语| 少妇人妻精品综合一区二区| 亚洲av成人精品一区久久| 国产精品久久电影中文字幕| 成人美女网站在线观看视频| 精品99又大又爽又粗少妇毛片| 久久人妻av系列| 国产午夜精品论理片| 三级国产精品欧美在线观看| av国产免费在线观看| 久久久久久国产a免费观看| 午夜老司机福利剧场| 日本黄色片子视频| 人人妻人人看人人澡| 熟女电影av网| 在线观看一区二区三区| 乱人视频在线观看| 老司机影院成人| 国产探花在线观看一区二区| 亚洲国产高清在线一区二区三| 色噜噜av男人的天堂激情| 国产精品爽爽va在线观看网站| 亚洲最大成人手机在线| 赤兔流量卡办理| 亚洲av成人精品一二三区| 午夜爱爱视频在线播放| 国产精品.久久久| 九九久久精品国产亚洲av麻豆| 国产视频首页在线观看| 精品一区二区三区视频在线| 亚洲国产精品久久男人天堂| 国产欧美日韩精品一区二区| 黄色一级大片看看| av播播在线观看一区| АⅤ资源中文在线天堂| 狠狠狠狠99中文字幕| 少妇熟女欧美另类| 99久久成人亚洲精品观看| 伊人久久精品亚洲午夜| 欧美又色又爽又黄视频| 最近最新中文字幕免费大全7| 91精品伊人久久大香线蕉| 极品教师在线视频| 日日摸夜夜添夜夜添av毛片| 精品一区二区三区人妻视频| 蜜臀久久99精品久久宅男| 久久久久久国产a免费观看| 男人和女人高潮做爰伦理| 亚洲天堂国产精品一区在线| 亚洲色图av天堂| 哪个播放器可以免费观看大片| 国产一区二区在线av高清观看| 国产在视频线精品| 免费黄网站久久成人精品| 亚洲欧美日韩东京热| 91久久精品电影网| 国产免费福利视频在线观看| 级片在线观看| 毛片一级片免费看久久久久| h日本视频在线播放| 嘟嘟电影网在线观看| 亚洲在线观看片| 在线播放无遮挡| 国产69精品久久久久777片| 成年免费大片在线观看| 大又大粗又爽又黄少妇毛片口| 欧美不卡视频在线免费观看| 天天躁夜夜躁狠狠久久av| 少妇的逼好多水| 亚洲欧美日韩卡通动漫| 国产人妻一区二区三区在| 成人亚洲精品av一区二区| 国产伦一二天堂av在线观看| 日韩一本色道免费dvd| www日本黄色视频网| 欧美一区二区亚洲| 青春草亚洲视频在线观看| 麻豆成人午夜福利视频| 天堂av国产一区二区熟女人妻| 成年女人看的毛片在线观看| 日本wwww免费看| 久久久精品欧美日韩精品| 午夜免费激情av| 久久草成人影院| 91狼人影院| 久久久亚洲精品成人影院| 亚洲精品色激情综合| 亚洲婷婷狠狠爱综合网| 美女高潮的动态| 欧美变态另类bdsm刘玥| 国产高清不卡午夜福利| 哪个播放器可以免费观看大片| 亚洲欧洲日产国产| 午夜免费激情av| 亚洲图色成人| 亚洲在久久综合| 91午夜精品亚洲一区二区三区| 长腿黑丝高跟| 成人午夜高清在线视频| 三级国产精品片| 夜夜爽夜夜爽视频| 亚洲成色77777| 神马国产精品三级电影在线观看| 国产精品久久久久久精品电影| 一夜夜www| 免费黄色在线免费观看| 高清在线视频一区二区三区 | 91精品国产九色| 亚洲国产精品国产精品| 99热这里只有是精品在线观看| 18禁裸乳无遮挡免费网站照片| 成年av动漫网址| 一级毛片aaaaaa免费看小| 两性午夜刺激爽爽歪歪视频在线观看| 看十八女毛片水多多多| 日韩精品青青久久久久久| 色综合亚洲欧美另类图片| 久久久精品94久久精品| 午夜福利视频1000在线观看| 国产精品一二三区在线看| 啦啦啦观看免费观看视频高清| 一个人看视频在线观看www免费| 久久热精品热| 精品午夜福利在线看| 日日摸夜夜添夜夜爱| 成人av在线播放网站| 级片在线观看| 看十八女毛片水多多多| 欧美日本亚洲视频在线播放| 国产成人精品婷婷| 日韩三级伦理在线观看| 草草在线视频免费看| 一夜夜www| 久久精品国产99精品国产亚洲性色| 91狼人影院| 国产精品无大码| 激情 狠狠 欧美| 亚洲成人久久爱视频| 国产在线男女| ponron亚洲| 久久久精品欧美日韩精品| 少妇熟女欧美另类| 99热这里只有精品一区| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 免费看美女性在线毛片视频| a级毛色黄片| 国产午夜精品论理片| 亚洲欧美日韩无卡精品| 女人十人毛片免费观看3o分钟| 免费看日本二区| 国语自产精品视频在线第100页| 亚洲精品亚洲一区二区| 九九热线精品视视频播放| 欧美xxxx性猛交bbbb| 久久久久免费精品人妻一区二区| 久久欧美精品欧美久久欧美| 少妇裸体淫交视频免费看高清| 国产69精品久久久久777片| 亚洲成av人片在线播放无| 亚洲欧美精品自产自拍| 少妇的逼好多水| 久久精品夜色国产| 特级一级黄色大片| 日韩精品有码人妻一区| 少妇猛男粗大的猛烈进出视频 | 两个人视频免费观看高清| 又爽又黄无遮挡网站| 日韩,欧美,国产一区二区三区 | 91在线精品国自产拍蜜月| 国产一级毛片在线| 高清在线视频一区二区三区 | 国产免费一级a男人的天堂| 爱豆传媒免费全集在线观看| 亚洲国产最新在线播放| 国产亚洲av片在线观看秒播厂 | 国产亚洲午夜精品一区二区久久 | 美女cb高潮喷水在线观看| 久久精品国产99精品国产亚洲性色| 联通29元200g的流量卡| 国产精品久久久久久久电影| 日韩大片免费观看网站 | 五月玫瑰六月丁香| 一边摸一边抽搐一进一小说| 国产大屁股一区二区在线视频| 国产熟女欧美一区二区| 欧美性猛交╳xxx乱大交人| 久久人人爽人人片av| 熟妇人妻久久中文字幕3abv| 日韩精品有码人妻一区| 日韩亚洲欧美综合| 九草在线视频观看| 国产一区亚洲一区在线观看| 久久久色成人| 人妻制服诱惑在线中文字幕| 国产私拍福利视频在线观看| 丝袜美腿在线中文| 少妇的逼水好多| 哪个播放器可以免费观看大片| 好男人视频免费观看在线| 麻豆一二三区av精品| 美女大奶头视频| 中文乱码字字幕精品一区二区三区 | 一个人看的www免费观看视频| 国产精品三级大全| 国产精品一二三区在线看| 免费av不卡在线播放| 激情 狠狠 欧美| 内地一区二区视频在线| 男的添女的下面高潮视频| 日韩欧美精品免费久久| 欧美潮喷喷水| 男人舔女人下体高潮全视频| 日本欧美国产在线视频| 一级毛片aaaaaa免费看小| 2021天堂中文幕一二区在线观| 能在线免费观看的黄片| 国产伦一二天堂av在线观看| 又粗又爽又猛毛片免费看| 中文字幕熟女人妻在线| 久久久国产成人免费| 成人二区视频| 免费av不卡在线播放| 免费播放大片免费观看视频在线观看 | 免费看av在线观看网站| 久久这里有精品视频免费| 三级国产精品欧美在线观看| 国产成人免费观看mmmm| 欧美色视频一区免费| 中文字幕久久专区| 人人妻人人澡欧美一区二区| 免费搜索国产男女视频| 国产大屁股一区二区在线视频| 在线a可以看的网站| 精品午夜福利在线看| 国产在视频线精品| 亚洲不卡免费看| av.在线天堂| 精品国内亚洲2022精品成人| 欧美人与善性xxx| 久久99热这里只有精品18| 午夜老司机福利剧场| 免费人成在线观看视频色| 日韩人妻高清精品专区| 日本免费a在线| 精品久久久久久成人av| 欧美日本亚洲视频在线播放| 男人和女人高潮做爰伦理| 一级二级三级毛片免费看| 亚洲天堂国产精品一区在线| 村上凉子中文字幕在线| 国产视频首页在线观看| 99久国产av精品| av在线蜜桃| 亚洲精品国产av成人精品| 国产伦在线观看视频一区| 中文字幕av成人在线电影| 久久久色成人| 欧美+日韩+精品| 日韩强制内射视频| 观看免费一级毛片| 美女国产视频在线观看| a级毛色黄片| 男人舔奶头视频| 精品久久久久久久末码| 尾随美女入室| 成人特级av手机在线观看| 久久精品国产鲁丝片午夜精品| 3wmmmm亚洲av在线观看| 精品一区二区三区视频在线| 成年版毛片免费区| 91在线精品国自产拍蜜月| 欧美一级a爱片免费观看看| 丰满人妻一区二区三区视频av| 国产精品国产三级国产专区5o | 久久久久久久午夜电影| 91精品国产九色| 日韩av在线免费看完整版不卡| 好男人在线观看高清免费视频| 国产精品一区二区三区四区免费观看| 国产精品三级大全| 免费人成在线观看视频色| 成人国产麻豆网| 成人亚洲精品av一区二区| 91午夜精品亚洲一区二区三区| 男插女下体视频免费在线播放| 在线观看66精品国产| 日本av手机在线免费观看| av黄色大香蕉| 亚洲欧美日韩无卡精品| 午夜精品在线福利| 国语自产精品视频在线第100页| 亚洲国产精品成人综合色| 日韩欧美在线乱码| 99久久中文字幕三级久久日本| 99热这里只有精品一区| 国产三级在线视频| 国产精品熟女久久久久浪| 国产精品电影一区二区三区|