• <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
    身高線與皺紋線
    最后的刺客免费高清国语| 国产欧美日韩一区二区精品| 在线观看一区二区三区| 欧美激情久久久久久爽电影| 婷婷丁香在线五月| а√天堂www在线а√下载| 深夜精品福利| 亚洲精品一区av在线观看| 亚洲欧美日韩东京热| 制服人妻中文乱码| 午夜激情福利司机影院| 一本久久中文字幕| 久久久精品欧美日韩精品| 日韩 欧美 亚洲 中文字幕| 国产精品久久久久久久电影 | 久久久久国产精品人妻aⅴ院| 久久久久精品国产欧美久久久| 黄色片一级片一级黄色片| 国产伦人伦偷精品视频| 99热只有精品国产| www日本黄色视频网| 一进一出抽搐gif免费好疼| 老熟妇乱子伦视频在线观看| 一级毛片高清免费大全| 亚洲男人的天堂狠狠| 黄色片一级片一级黄色片| 一边摸一边抽搐一进一小说| 听说在线观看完整版免费高清| 国产精品久久久久久精品电影| 波多野结衣高清无吗| 国产乱人伦免费视频| 色综合婷婷激情| 欧美极品一区二区三区四区| 毛片女人毛片| 国产99白浆流出| 桃红色精品国产亚洲av| 97超视频在线观看视频| 国产爱豆传媒在线观看| 亚洲人成网站高清观看| 国产亚洲欧美98| 麻豆国产av国片精品| 最近最新中文字幕大全免费视频| 在线观看一区二区三区| 男人舔奶头视频| 亚洲精品一卡2卡三卡4卡5卡| 18+在线观看网站| 最新美女视频免费是黄的| 一本一本综合久久| 国产黄色小视频在线观看| 男女下面进入的视频免费午夜| 亚洲av五月六月丁香网| 国产av一区在线观看免费| 欧美+日韩+精品| 欧美成人性av电影在线观看| 亚洲精品美女久久久久99蜜臀| 18禁黄网站禁片免费观看直播| 女人高潮潮喷娇喘18禁视频| 一个人看的www免费观看视频| av福利片在线观看| 国产激情偷乱视频一区二区| 成年免费大片在线观看| 国产真实乱freesex| 亚洲精华国产精华精| 天天添夜夜摸| 一二三四社区在线视频社区8| 麻豆久久精品国产亚洲av| 国内精品一区二区在线观看| 看免费av毛片| 欧美大码av| 99热这里只有精品一区| 久99久视频精品免费| 久久亚洲真实| 超碰av人人做人人爽久久 | 88av欧美| 亚洲欧美日韩无卡精品| 91久久精品电影网| 国产毛片a区久久久久| 亚洲午夜理论影院| 最新美女视频免费是黄的| 人妻久久中文字幕网| aaaaa片日本免费| 国模一区二区三区四区视频| 亚洲精品影视一区二区三区av| 在线观看美女被高潮喷水网站 | 色综合婷婷激情| 高潮久久久久久久久久久不卡| 精品日产1卡2卡| 国产中年淑女户外野战色| 狠狠狠狠99中文字幕| 亚洲av电影在线进入| 十八禁网站免费在线| 免费电影在线观看免费观看| 热99在线观看视频| 好男人电影高清在线观看| 亚洲色图av天堂| 韩国av一区二区三区四区| 蜜桃久久精品国产亚洲av| 欧美日韩一级在线毛片| 最近在线观看免费完整版| 免费观看的影片在线观看| 99久久无色码亚洲精品果冻| 久久久久久久精品吃奶| 舔av片在线| 九九热线精品视视频播放| 成年女人永久免费观看视频| 亚洲精品色激情综合| 男女做爰动态图高潮gif福利片| 亚洲自拍偷在线| 成人特级黄色片久久久久久久| 又黄又爽又免费观看的视频| 成人欧美大片| 成人18禁在线播放| 亚洲国产色片| 99在线视频只有这里精品首页| 久久亚洲真实| 国产伦人伦偷精品视频| 舔av片在线| 床上黄色一级片| 亚洲国产欧洲综合997久久,| 国产精品一及| 18禁黄网站禁片免费观看直播| 大型黄色视频在线免费观看| 麻豆久久精品国产亚洲av| 精品国产超薄肉色丝袜足j| 欧美区成人在线视频| 内地一区二区视频在线| 欧美黄色淫秽网站| 欧美中文日本在线观看视频| 国产亚洲欧美98| 成年女人毛片免费观看观看9| 桃色一区二区三区在线观看| 亚洲成人久久爱视频| 欧美在线一区亚洲| 国产免费av片在线观看野外av| 国产一区二区在线观看日韩 | 欧美黑人欧美精品刺激| 日本与韩国留学比较| 亚洲一区二区三区色噜噜| 精品福利观看| 成年版毛片免费区| av天堂中文字幕网| 老司机在亚洲福利影院| 男女之事视频高清在线观看| 床上黄色一级片| 亚洲精品国产精品久久久不卡| 中文字幕熟女人妻在线| 久久久久久久精品吃奶| 91久久精品国产一区二区成人 | 热99re8久久精品国产| 精品一区二区三区视频在线 | 婷婷亚洲欧美| 欧美日韩瑟瑟在线播放| 亚洲av五月六月丁香网| av福利片在线观看| 成人国产综合亚洲| 免费无遮挡裸体视频| av福利片在线观看| av在线蜜桃| 在线观看66精品国产| 亚洲精品一卡2卡三卡4卡5卡| 免费无遮挡裸体视频| 一进一出抽搐gif免费好疼| 无遮挡黄片免费观看| 男人舔女人下体高潮全视频| 一进一出好大好爽视频| 久久这里只有精品中国| 日本与韩国留学比较| 亚洲专区中文字幕在线| 亚洲国产日韩欧美精品在线观看 | 看片在线看免费视频| 国产高清视频在线播放一区| 久久精品综合一区二区三区| 动漫黄色视频在线观看| 脱女人内裤的视频| 美女 人体艺术 gogo| 国产黄色小视频在线观看| 亚洲一区高清亚洲精品| a在线观看视频网站| 国产精品爽爽va在线观看网站| 五月伊人婷婷丁香| 悠悠久久av| 日韩精品中文字幕看吧| 国产在线精品亚洲第一网站| 在线观看一区二区三区| 少妇高潮的动态图| 国产淫片久久久久久久久 | 老汉色av国产亚洲站长工具| 黄色女人牲交| 最新中文字幕久久久久| 亚洲一区二区三区色噜噜| 白带黄色成豆腐渣| 一级a爱片免费观看的视频| 亚洲性夜色夜夜综合| 亚洲人成电影免费在线| 中文字幕高清在线视频| 日本免费一区二区三区高清不卡| 九色成人免费人妻av| 国产真实乱freesex| 国产 一区 欧美 日韩| 岛国在线观看网站| 丁香欧美五月| 日本免费a在线| 欧美丝袜亚洲另类 | 欧美一区二区亚洲| 可以在线观看的亚洲视频| 亚洲在线自拍视频| 91久久精品电影网| 国产成人福利小说| 日韩亚洲欧美综合| 久久久久国产精品人妻aⅴ院| 男女午夜视频在线观看| 亚洲狠狠婷婷综合久久图片| aaaaa片日本免费| 一级毛片高清免费大全| 国产乱人视频| 日韩av在线大香蕉| 欧美国产日韩亚洲一区| 亚洲电影在线观看av| 日本精品一区二区三区蜜桃| 亚洲成av人片免费观看| 欧美绝顶高潮抽搐喷水| 亚洲成人精品中文字幕电影| 天美传媒精品一区二区| 国产精品久久电影中文字幕| 国产精品1区2区在线观看.| 欧美xxxx黑人xx丫x性爽| 国产精品一区二区三区四区久久| 1024手机看黄色片| 国产在线精品亚洲第一网站| 麻豆成人av在线观看| 丰满的人妻完整版| 一本一本综合久久| 午夜视频国产福利| 日韩精品青青久久久久久| 床上黄色一级片| 久久精品国产99精品国产亚洲性色| 午夜免费男女啪啪视频观看 | 深爱激情五月婷婷| 国产日本99.免费观看| 99久久精品一区二区三区| 一级毛片高清免费大全| 亚洲精品乱码久久久v下载方式 | 亚洲av中文字字幕乱码综合| 中国美女看黄片| 色综合站精品国产| 深夜精品福利| 午夜福利高清视频| 免费av不卡在线播放| 日日夜夜操网爽| 最新在线观看一区二区三区| 欧美日韩一级在线毛片| 一区二区三区免费毛片| 免费人成视频x8x8入口观看| 嫩草影院入口| 亚洲成a人片在线一区二区| 在线播放无遮挡| 婷婷精品国产亚洲av在线| 叶爱在线成人免费视频播放| 日本黄色片子视频| 又黄又爽又免费观看的视频| 宅男免费午夜| 国产欧美日韩一区二区精品| 丁香六月欧美| 亚洲av成人不卡在线观看播放网| 俺也久久电影网| 国产aⅴ精品一区二区三区波| 又紧又爽又黄一区二区| 亚洲av免费在线观看| 美女cb高潮喷水在线观看| 夜夜躁狠狠躁天天躁| 老熟妇仑乱视频hdxx| 九色成人免费人妻av| 欧美色欧美亚洲另类二区| 人人妻人人看人人澡| av中文乱码字幕在线| 18禁美女被吸乳视频| av片东京热男人的天堂| www日本在线高清视频| 国产国拍精品亚洲av在线观看 | 欧美中文综合在线视频| 午夜福利在线观看免费完整高清在 | 久9热在线精品视频| 国产精品一区二区免费欧美| 午夜免费成人在线视频| 国产伦人伦偷精品视频| 成熟少妇高潮喷水视频| 国产三级在线视频| 两个人视频免费观看高清| 人妻久久中文字幕网| 国内精品一区二区在线观看| 午夜免费观看网址| 午夜激情欧美在线| 天天一区二区日本电影三级| 日韩欧美一区二区三区在线观看| 亚洲最大成人手机在线| 麻豆成人av在线观看| АⅤ资源中文在线天堂| 久久久国产成人精品二区| 757午夜福利合集在线观看| 国产精品av视频在线免费观看| 国模一区二区三区四区视频| 啦啦啦观看免费观看视频高清| 欧美日韩综合久久久久久 | 色av中文字幕| 国产视频一区二区在线看| 长腿黑丝高跟| 亚洲一区高清亚洲精品| 亚洲成人免费电影在线观看| 欧美激情在线99| 中文字幕人妻熟人妻熟丝袜美 | 午夜免费激情av| 丰满人妻熟妇乱又伦精品不卡| av天堂中文字幕网| 18禁在线播放成人免费| 精品99又大又爽又粗少妇毛片 | 成年版毛片免费区| av在线天堂中文字幕| 亚洲人成网站高清观看| 老熟妇乱子伦视频在线观看| 51午夜福利影视在线观看| 麻豆一二三区av精品| 欧美日本亚洲视频在线播放| 88av欧美| 精品日产1卡2卡| 在线a可以看的网站| 91麻豆精品激情在线观看国产| 亚洲av免费高清在线观看| 看片在线看免费视频| 亚洲 国产 在线| 丰满乱子伦码专区| 亚洲欧美激情综合另类| 嫩草影视91久久| 69av精品久久久久久| 真实男女啪啪啪动态图| 国产午夜精品久久久久久一区二区三区 | 亚洲精华国产精华精| 欧美一级a爱片免费观看看| 亚洲av成人不卡在线观看播放网| 精品免费久久久久久久清纯| 国产在线精品亚洲第一网站| 精品久久久久久,| www.www免费av| 精品午夜福利视频在线观看一区| 日日摸夜夜添夜夜添小说| 久久久久久久久中文| 欧美乱码精品一区二区三区| 精品电影一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 久久国产精品人妻蜜桃| 国产三级在线视频| 国产成人影院久久av| 一级毛片女人18水好多| 亚洲av成人av| 婷婷丁香在线五月| 日韩欧美 国产精品| 舔av片在线| 久久久国产成人精品二区| 精品电影一区二区在线| 免费观看人在逋| 天堂影院成人在线观看| 午夜免费成人在线视频| 色播亚洲综合网| 欧美乱色亚洲激情| 99久久精品热视频| 美女 人体艺术 gogo| 99久久九九国产精品国产免费| 99在线人妻在线中文字幕| 麻豆国产av国片精品| 黄色成人免费大全| 成年女人永久免费观看视频| 亚洲久久久久久中文字幕| 1000部很黄的大片| 免费人成视频x8x8入口观看| 日韩精品中文字幕看吧| 少妇的逼好多水| 色av中文字幕| av天堂在线播放| 国产精品久久久久久人妻精品电影| 久久精品国产亚洲av涩爱 | 一区二区三区高清视频在线| 99久久无色码亚洲精品果冻| 亚洲国产日韩欧美精品在线观看 | 国产色爽女视频免费观看| 又粗又爽又猛毛片免费看| 嫩草影院入口| 久久久久免费精品人妻一区二区| 欧美午夜高清在线| 亚洲人成网站在线播| 亚洲专区中文字幕在线| 一本久久中文字幕| 激情在线观看视频在线高清| 欧美xxxx黑人xx丫x性爽| 韩国av一区二区三区四区| 88av欧美| 国产成人啪精品午夜网站| 国产亚洲精品久久久久久毛片| 观看美女的网站| 一级毛片高清免费大全| av国产免费在线观看| 美女 人体艺术 gogo| 久久香蕉精品热| 琪琪午夜伦伦电影理论片6080| 午夜福利成人在线免费观看| 一进一出抽搐gif免费好疼| 啪啪无遮挡十八禁网站| 欧美乱色亚洲激情| 五月玫瑰六月丁香| 亚洲天堂国产精品一区在线| 在线视频色国产色| 最近在线观看免费完整版| 淫妇啪啪啪对白视频| 五月伊人婷婷丁香| 99久久精品国产亚洲精品| 日本成人三级电影网站| 亚洲av一区综合| 九色国产91popny在线| 久久中文看片网| 婷婷丁香在线五月| 人人妻,人人澡人人爽秒播| 一本久久中文字幕| 在线国产一区二区在线| 一二三四社区在线视频社区8| 精品99又大又爽又粗少妇毛片 | 久久久久久久亚洲中文字幕 | 搡女人真爽免费视频火全软件 | 精品久久久久久久末码| 免费一级毛片在线播放高清视频| 国产不卡一卡二| 国产精品香港三级国产av潘金莲| 性色avwww在线观看| 欧美一级a爱片免费观看看| 无限看片的www在线观看| 3wmmmm亚洲av在线观看| 男人的好看免费观看在线视频| 成人永久免费在线观看视频| 丰满的人妻完整版| e午夜精品久久久久久久| 欧美日韩精品网址| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添小说| 欧美性猛交╳xxx乱大交人| av专区在线播放| 19禁男女啪啪无遮挡网站| 婷婷丁香在线五月| 国产欧美日韩一区二区三| 尤物成人国产欧美一区二区三区| 国产一区在线观看成人免费| 亚洲av电影不卡..在线观看| 中文字幕av成人在线电影| 国产精品嫩草影院av在线观看 | 国产一级毛片七仙女欲春2| 日本黄色视频三级网站网址| 美女黄网站色视频| 国产精品 欧美亚洲| 丁香欧美五月| 色尼玛亚洲综合影院| 精品久久久久久,| 好男人在线观看高清免费视频| 好男人电影高清在线观看| 国产精品 国内视频| 小说图片视频综合网站| 亚洲av第一区精品v没综合| 国产成+人综合+亚洲专区| 国产精品免费一区二区三区在线| 成人18禁在线播放| 亚洲精品成人久久久久久| 欧美在线一区亚洲| 日韩有码中文字幕| 十八禁人妻一区二区| 夜夜夜夜夜久久久久| 脱女人内裤的视频| 国产精品久久久久久久电影 | 法律面前人人平等表现在哪些方面| 国产老妇女一区| 88av欧美| 国产精品亚洲av一区麻豆| 在线观看av片永久免费下载| 欧美性感艳星| 欧美成人免费av一区二区三区| 激情在线观看视频在线高清| 男人舔女人下体高潮全视频| 99国产精品一区二区蜜桃av| 亚洲国产欧美网| 久久精品夜夜夜夜夜久久蜜豆| 一个人观看的视频www高清免费观看| av黄色大香蕉| 午夜福利高清视频| 国产精品自产拍在线观看55亚洲| av专区在线播放| 久久亚洲精品不卡| 特大巨黑吊av在线直播| 网址你懂的国产日韩在线| 免费搜索国产男女视频| 国产精品国产高清国产av| 亚洲精华国产精华精| 色综合婷婷激情| 欧美性猛交黑人性爽| 好看av亚洲va欧美ⅴa在| 亚洲av日韩精品久久久久久密| 国产精品永久免费网站| 色噜噜av男人的天堂激情| 一级黄片播放器| 亚洲最大成人中文| 日韩欧美免费精品| 一级毛片女人18水好多| 国产精品一区二区三区四区久久| 国产亚洲精品av在线| 亚洲av电影在线进入| 国产熟女xx| 日韩欧美三级三区| 国产v大片淫在线免费观看| 久久婷婷人人爽人人干人人爱| av女优亚洲男人天堂| 亚洲人成网站高清观看| 极品教师在线免费播放| 免费av观看视频| 亚洲精品日韩av片在线观看 | 不卡一级毛片| 欧美乱妇无乱码| 免费看a级黄色片| 欧美乱妇无乱码| 最好的美女福利视频网| 欧美色欧美亚洲另类二区| 国产精品久久久久久久电影 | 久久久久久久久中文| 女人十人毛片免费观看3o分钟| 国产探花极品一区二区| 国产精华一区二区三区| 久久精品91蜜桃| 午夜福利在线观看吧| 国内少妇人妻偷人精品xxx网站| tocl精华| 日韩欧美精品v在线| 在线免费观看不下载黄p国产 | 女人被狂操c到高潮| 伊人久久大香线蕉亚洲五| 午夜老司机福利剧场| 亚洲成人免费电影在线观看| 99riav亚洲国产免费| 日韩欧美国产一区二区入口| 99热6这里只有精品| 亚洲av五月六月丁香网| 欧美日韩一级在线毛片| 成人特级黄色片久久久久久久| 一级作爱视频免费观看| 老熟妇仑乱视频hdxx| 午夜福利成人在线免费观看| 人人妻人人澡欧美一区二区| 中文字幕人妻熟人妻熟丝袜美 | 欧美另类亚洲清纯唯美| 亚洲熟妇熟女久久| 久久久久久久亚洲中文字幕 | 国产黄a三级三级三级人| 久久精品国产综合久久久| 亚洲国产欧美人成| 国产黄片美女视频| 噜噜噜噜噜久久久久久91| 在线观看66精品国产| 一边摸一边抽搐一进一小说| 精品人妻偷拍中文字幕| x7x7x7水蜜桃| 国产乱人视频| 国产精品 欧美亚洲| 黄色女人牲交| 亚洲午夜理论影院| 性色av乱码一区二区三区2| 欧美三级亚洲精品| 日韩亚洲欧美综合| 中文字幕人成人乱码亚洲影| 美女 人体艺术 gogo| 久久久成人免费电影| 国产av一区在线观看免费| 久久久国产精品麻豆| 国产亚洲精品一区二区www| 久久久久国内视频| 国产一级毛片七仙女欲春2| 狂野欧美激情性xxxx| 亚洲精品亚洲一区二区| 日本在线视频免费播放| 亚洲美女视频黄频| 在线播放国产精品三级| 午夜免费观看网址| 国内久久婷婷六月综合欲色啪| 久久久久久大精品| 一个人看视频在线观看www免费 | 人妻久久中文字幕网| 美女免费视频网站| tocl精华| 91字幕亚洲| 欧美在线一区亚洲| 在线播放无遮挡| 搡老妇女老女人老熟妇| 国产伦精品一区二区三区四那| 人人妻人人澡欧美一区二区| 亚洲精品影视一区二区三区av| 午夜免费男女啪啪视频观看 | 国产欧美日韩一区二区精品| 亚洲无线观看免费| 悠悠久久av| 国内精品久久久久久久电影| av天堂中文字幕网| 久久6这里有精品| 国产色婷婷99| 波多野结衣高清作品| 国产精品永久免费网站| 日韩欧美国产一区二区入口| 久久久久性生活片| 国产熟女xx| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 精品久久久久久久毛片微露脸| 国产又黄又爽又无遮挡在线| 久久精品人妻少妇| 桃红色精品国产亚洲av| 变态另类成人亚洲欧美熟女|