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

    草魚單養(yǎng)和混養(yǎng)池塘的水質(zhì)與生物組成特征

    2023-04-29 20:08:07肖述文劉興國陸詩敏趙宇曦顧兆俊周潤(rùn)鋒
    水生態(tài)學(xué)雜志 2023年6期
    關(guān)鍵詞:養(yǎng)殖模式

    肖述文 劉興國 陸詩敏 趙宇曦 顧兆俊 周潤(rùn)鋒

    摘要:為比較單養(yǎng)、混養(yǎng)草魚(Ctenopharyngodon idella)養(yǎng)殖池塘的水質(zhì)與生物組成特點(diǎn),采取水質(zhì)分析、環(huán)境DNA與傳統(tǒng)鑒別方法對(duì)草魚單養(yǎng)、混養(yǎng)(草魚:混養(yǎng)魚為80:20)兩種池塘的水質(zhì)變化、浮游生物、底棲生物、菌群結(jié)構(gòu)進(jìn)行了分析。結(jié)果顯示:混養(yǎng)池塘的水質(zhì)優(yōu)于單養(yǎng)池塘,混養(yǎng)池塘水體中總氮(TN)、硝態(tài)氮(NO[-3]-N)、氨氮(NH[+4]-N)、亞硝態(tài)氮(NO[-2]-N)的濃度比單養(yǎng)池塘分別低10.15%、3.78%、5.07%、80.18%,總磷(TP)和活性磷(SRP)的濃度分別低27.14%和56.26%;兩類池塘中浮游植物均以綠藻門(Chlorophyta)、藍(lán)藻門(Cyanophyta)、隱藻門(Cryptophyta)為優(yōu)勢(shì)種,但單養(yǎng)池塘中的藻類密度為3.00×107 個(gè)/L,低于混養(yǎng)池塘1.04×108 個(gè)/L;兩類池塘中的浮游動(dòng)物均以輪蟲和原生動(dòng)物為優(yōu)勢(shì)種,枝角類和橈足類生物數(shù)量較少,單養(yǎng)池塘中浮游動(dòng)物密度高于混養(yǎng)池塘;在底棲動(dòng)物方面,單養(yǎng)池塘存在螺類、水蚯蚓和搖蚊幼蟲,而混養(yǎng)池塘僅有螺類和搖蚊幼蟲。在菌群組成方面,單養(yǎng)池塘水體中以厚壁菌門(Firmicutes)為優(yōu)勢(shì)類群,混養(yǎng)池塘水體中以變形菌門(Proteobacteria)為優(yōu)勢(shì)類群;但在兩種池塘底泥中,均以變形菌門為優(yōu)勢(shì)類群。研究表明,草魚混養(yǎng)有利于改善養(yǎng)殖池塘水質(zhì),增加浮游植物豐富度,改變養(yǎng)殖水體菌群的結(jié)構(gòu)。研究結(jié)果為優(yōu)化草魚池塘養(yǎng)殖結(jié)構(gòu)、改善水質(zhì)、構(gòu)建高效池塘養(yǎng)殖模式提供了依據(jù)。

    關(guān)鍵詞:草魚池塘;養(yǎng)殖模式;水質(zhì)指標(biāo);生物組成;環(huán)境DNA

    中圖分類號(hào):S964.3,Q178.1? ? ? ? 文獻(xiàn)標(biāo)志碼:A? ? ? ? 文章編號(hào):1674-3075(2023)06-0079-09

    草魚(Ctenopharyngodon idella)是中國年產(chǎn)量最高的淡水魚類,有池塘養(yǎng)殖、網(wǎng)箱養(yǎng)殖、水庫養(yǎng)殖、稻田養(yǎng)殖等多種模式,其中池塘養(yǎng)殖占草魚養(yǎng)殖總產(chǎn)量的74%(國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系,2021);在華南、華東、華中等草魚主產(chǎn)區(qū),其池塘養(yǎng)殖模式主要有單養(yǎng)和混養(yǎng)兩種(常華,2009;楊帆,2016;陳小江和汪祈超,2016;梅順和任青松,2017)。草魚單養(yǎng)模式由于放養(yǎng)密度大,輪捕輪放能獲得較高的收益,但存在水質(zhì)管理難度大、易暴發(fā)疾病等問題(劉玲儀和金有坤,1994;高攀等,2009;劉朋, 2012);草魚混養(yǎng)模式主要為80:20方式(草魚:混養(yǎng)魚),混養(yǎng)品種多為鰱(Hypophthalmichthys molitrix)、鳙(Aristichthys nobilis)、鯽(Carassius auratus)等(占家智和羊茜,2015);其中鰱、鳙可攝食浮游動(dòng)植物和有機(jī)碎屑,抑制藍(lán)綠藻的過度生長(zhǎng)并改善養(yǎng)殖水質(zhì)(孫云飛,2013),鯽作為底層魚類能夠充分利用水體空間,攝食底層殘餌。近年來,隨著養(yǎng)殖尾水問題加重,這種多營(yíng)養(yǎng)層級(jí)的草魚混養(yǎng)模式被廣泛應(yīng)用,通過充分利用系統(tǒng)空間和資源,提高了系統(tǒng)內(nèi)營(yíng)養(yǎng)物質(zhì)的循環(huán)利用率,防止了病害發(fā)生并改善水質(zhì)(董雙林,2014;周波,2015)。

    環(huán)境DNA(Environmental DNA, eDNA)方法是從環(huán)境樣本中提取DNA,針對(duì)目標(biāo)類群設(shè)計(jì)通用引物,通過PCR擴(kuò)增并結(jié)合高通量測(cè)序,對(duì)環(huán)境中的目標(biāo)物種進(jìn)行識(shí)別和分析(Taberlet et al, 2012)。環(huán)境DNA方法能夠測(cè)定傳統(tǒng)方法中難以觀察到的物種,不會(huì)對(duì)物種和生態(tài)系統(tǒng)造成破壞,也能更加準(zhǔn)確、快速地判斷目標(biāo)物種(巴日斯,2018)。近年來,環(huán)境DNA分析技術(shù)已被廣泛應(yīng)用于海洋、湖泊、池塘等生態(tài)系統(tǒng)中的無脊椎動(dòng)物、浮游植物、細(xì)菌等生物檢測(cè)(李萌等,2019;舒璐等,2020;張麗娟等,2021;許鄭超,2021;Singer et al, 2021);而我國應(yīng)用環(huán)境DNA技術(shù)調(diào)查池塘生態(tài)系統(tǒng)中生物的研究并不常見。目前主要采用傳統(tǒng)方法來監(jiān)測(cè)池塘中浮游動(dòng)植物及底棲生物等(李瑞嬌,2014;葉曉彤等,2020;陳佳林等,2021),但傳統(tǒng)方法存在監(jiān)測(cè)方法復(fù)雜、物種難以鑒別等問題,故本文將兩種方法相結(jié)合以快速、有效地開展池塘浮游生物的監(jiān)測(cè)工作。

    目前,我國草魚池塘養(yǎng)殖依然存在飼料利用率低、水質(zhì)調(diào)控難等問題,亟需系統(tǒng)研究不同養(yǎng)殖模式的生態(tài)特征,優(yōu)化池塘生態(tài)系統(tǒng)(Li et al,2019;Niu et al,2023)。本文采用水質(zhì)分析、環(huán)境DNA技術(shù)及傳統(tǒng)鑒別方法,對(duì)單養(yǎng)和混養(yǎng)草魚養(yǎng)殖模式中的水質(zhì)、浮游生物、底棲生物、細(xì)菌等進(jìn)行研究,旨在為評(píng)價(jià)草魚養(yǎng)殖模式、構(gòu)建高效草魚池塘系統(tǒng)提供依據(jù),并對(duì)環(huán)境DNA方法監(jiān)測(cè)池塘生態(tài)系統(tǒng)中出現(xiàn)的問題進(jìn)行探究。

    1? ?材料與方法

    1.1? ?試驗(yàn)池塘

    在上海市松江區(qū)三泖水產(chǎn)養(yǎng)殖場(chǎng)選取兩類草魚養(yǎng)殖池塘進(jìn)行試驗(yàn),池塘規(guī)格為100 m × 50 m,平均水深1.6 m,水面積0.47 hm2。

    單養(yǎng)池塘放養(yǎng)0.57 kg/尾的草魚1 322尾,出于實(shí)際養(yǎng)殖需要,考慮到鰱并不攝食飼料,因此在單養(yǎng)塘投入0.15 kg/尾的鰱277尾;混養(yǎng)(80:20)池塘放養(yǎng)草魚897尾、0.50 kg/尾,鰱123尾、1.25 kg/尾,鳙529尾、0.26 kg/尾,鯽158尾、0.09 kg/尾,鳙魚苗2 000尾、0.5 g/尾,鯽魚苗1 000尾、0.25 g/尾。

    試驗(yàn)期間的增氧、投飼等管理方式一致,養(yǎng)殖期間不施肥、不噴灑殺蟲劑。

    1.2? ?樣本采集及測(cè)定

    選取池塘對(duì)角線上的3個(gè)點(diǎn),使用3點(diǎn)法對(duì)池塘水體及底泥樣品進(jìn)行采集,測(cè)定理化指標(biāo)并采集池塘中的生物樣本。水樣使用5 L有機(jī)玻璃采水器于池塘1 m深處進(jìn)行采集,底泥使用規(guī)格為1/40 m2彼得森采泥器抓取。

    水體理化指標(biāo)采取國標(biāo)法測(cè)定??偟═N)采用GB 11894-89(國家環(huán)境保護(hù)局規(guī)劃標(biāo)準(zhǔn)處,1989)過硫酸鉀氧化紫外分光光度法,硝態(tài)氮(NO3--N)采用HJ/T346-2007(國家環(huán)境保護(hù)總局科技標(biāo)準(zhǔn)司,2007)紫外分光光度法,氨氮(NH4+-N)采用HJ 535-2009(國家環(huán)境保護(hù)總局科技標(biāo)準(zhǔn)司,2009)納氏試劑光度法,亞硝態(tài)氮(NO2--N)采用GB 7493-87(國家環(huán)境保護(hù)局規(guī)劃標(biāo)準(zhǔn)處,1987)N-(1-萘基)-乙二胺光度法,總磷 (TP)、活性磷(SRP)采用GB 11893-89(國家環(huán)境保護(hù)局標(biāo)準(zhǔn)處,1989)鉬銻抗分光光度法。

    底泥總氮、總磷采取水浸法預(yù)處理,濾液采取國標(biāo)法測(cè)定。

    浮游生物與底棲生物的采樣定量參考《水和廢水監(jiān)測(cè)分析方法》(國家環(huán)境保護(hù)總局,2002)。

    1.3? ?環(huán)境DNA樣品

    3點(diǎn)法對(duì)池塘水樣和底泥取樣后,等體積均勻混合,取1.0 L水體及50.0 g底泥各為1個(gè)檢測(cè)樣本。每100 mL水樣使用0.22 μm的濾膜進(jìn)行抽濾處理,濾膜和底泥樣本置于-20℃冷凍保存,送至上海派森諾生物科技有限公司提取DNA并進(jìn)行后續(xù)測(cè)序分析。

    分別選取引物16S V3V4、18S V4、COI用于細(xì)菌、浮游植物、浮游動(dòng)物的PCR擴(kuò)增。細(xì)菌擴(kuò)增選取16S V3V4區(qū)域,上游引物:ACTCCTACGGGAGGCAGCA,下游引物:GGACTACNVGGGTWTCTAAT;浮游植物選取18S V4區(qū)域,上游引物:CCAGCASCYGCGGTAATTCC,下游引物:ACTTTCGTTCTTGATYRA;浮游動(dòng)物選取COI基因上的區(qū)域,上游引物:GGWACWGGWTGAACWGTWTAYCCYCC,下游引物:TANACYTCNGGRTGNCCRAARAAYCA。

    擴(kuò)增體系(25 μL):5 μL 5×reaction緩沖液,5 μL 5×GC緩沖液,2 μL dNTPs(2.5 mM),正反引物(10 μM)各1 μL,2 μL DNA模板(20 ng/μL),8.75 μL ddH2O,0.25 μL Q5 DNA聚合酶。反應(yīng)條件為98℃預(yù)變性2 min,98℃變性15 s、55℃退火30 s、72℃延伸30 s、25~30個(gè)循環(huán),72℃最后延伸5 min。

    1.2%瓊脂糖凝膠電泳檢測(cè)PCR擴(kuò)增結(jié)果,之后通過Illumina NovaSeq測(cè)序平臺(tái)進(jìn)行高通量測(cè)序和序列處理工作。

    1.4? ?數(shù)據(jù)處理與分析

    測(cè)序的原始數(shù)據(jù)利用Vsearch軟件進(jìn)行去引物、拼接、質(zhì)量過濾、去重、去嵌合體、聚類等步驟后,將得到的序列與數(shù)據(jù)庫進(jìn)行對(duì)比,進(jìn)行α-多樣性指數(shù)分析。細(xì)菌選取Silvar數(shù)據(jù)庫,浮游植物和浮游動(dòng)物選取NCBI數(shù)據(jù)庫。

    α-多樣性指數(shù)包括物種豐富度指數(shù)(Chao1指數(shù))、物種多樣性指數(shù)(Shannon指數(shù))、覆蓋率(Good's Coverage),均使用QIIME軟件計(jì)算。

    利用Excel對(duì)數(shù)據(jù)進(jìn)行初步統(tǒng)計(jì),通過SPSS 26.0軟件對(duì)池塘數(shù)據(jù)進(jìn)行獨(dú)立樣本T檢驗(yàn)分析,使用Origin 2016作圖。

    2? ?結(jié)果與分析

    2.1? ?池塘環(huán)境氮磷濃度

    對(duì)試驗(yàn)期間兩池塘水質(zhì)進(jìn)行檢測(cè),單養(yǎng)塘中的總氮、硝態(tài)氮、氨氮、亞硝態(tài)氮的濃度分別為(1.633± 0.018)、(0.164±0.015)、(0.207±0.016)、(0.031±0.004)mg/L,混養(yǎng)塘中的對(duì)應(yīng)濃度分別為(1.368±0.083)、(0.158±0.029)、(0.197±0.0240)、(0.004±0.002)mg/L,混養(yǎng)池塘水體中總氮、硝態(tài)氮、氨氮、亞硝態(tài)氮的濃度比單養(yǎng)池塘分別低10.15%、3.78%、5.07%、80.18%(圖1-a)。兩池塘中的總氮、亞硝態(tài)氮濃度差異顯著(P<0.05),硝態(tài)氮、氨氮濃度差異不顯著(P>0.05)。

    混養(yǎng)池塘中總磷和活性磷的濃度比單養(yǎng)池塘低27.14%和56.26%(圖1-b)。單養(yǎng)池塘總磷和活性磷的濃度分別為(0.788±0.204)、(0.373 ± 0.075)mg/L,混養(yǎng)池塘為(0.403±0.022)、(0.112 ± 0.018)mg/L,兩池塘中的總氮和活性磷濃度均存在顯著差異(P<0.05)。

    兩池塘底泥中的氮、磷含量檢測(cè)結(jié)果見圖2。單養(yǎng)池塘3個(gè)采樣點(diǎn)(點(diǎn)1、點(diǎn)2、點(diǎn)3)的總氮濃度分別為9.509、10.484、14.651 mg/kg,混養(yǎng)池塘分別為10.774、7.192、7.772 mg/kg;單養(yǎng)池塘3個(gè)采樣點(diǎn)的總磷濃度分別為1.100、0.670、4.606 mg/kg,混養(yǎng)池塘分別為0.273、0.094、3.869 mg/kg。兩池塘各區(qū)域總氮、總磷的含量均有所差異,總氮在各區(qū)域分布無明顯規(guī)律而總磷的含量差異較大,位于出水口附近區(qū)域的總磷含量尤其高,可能是由于飼料中未分解的磷元素沉積到底泥并被水流沖刷到出水口造成的。

    2.2? ?浮游植物相對(duì)豐度和密度

    采用環(huán)境DNA檢測(cè)兩池塘的浮游植物組成,其相對(duì)豐度結(jié)果見圖3。池塘中的主要藻類為綠藻門(Chlorophyta)、甲藻門(Dinophyta)、隱藻門(Cryptophyta)、硅藻門(Bacillariophyta)、金藻門(Chrysophyta)、黃藻門(Xanthophyta)、異鞭藻門(Heterokontophyta)。

    草魚單養(yǎng)池塘共計(jì)9門14綱17目,以綠藻門(31.06%)為優(yōu)勢(shì)類群,門水平檢測(cè)出來的藻類達(dá)41.73%。綱分類水平檢測(cè)出來的藻類占72.07%,以隱藻綱(Cryptophyceae)30.85%為優(yōu)勢(shì),其次是綠藻綱(Chlorophyceae)25.97%、圓篩藻綱(Coscinodiscophyceae)6.32%、共球藻綱(Trebouxiophyceae)4.57%、甲藻綱(Dinophyceae)2.22%。目分類水平檢出藻類占78.28%,以隱鞭藻目(Cryptomonadales)30.85%為優(yōu)勢(shì)。

    混養(yǎng)池塘共計(jì)10門17綱19目,以甲藻門(26.93%)為優(yōu)勢(shì)類群,門水平檢測(cè)出來的藻類可達(dá)36.30%。綱分類水平檢測(cè)出來的藻類占74.60%,隱藻綱38.68%、甲藻綱26.93%、綠藻綱4.81%、圓篩藻綱1.50%、共球藻綱0.95%。目分類水平檢出藻類占56.30%,以隱鞭藻目38.68%為優(yōu)勢(shì)。

    試驗(yàn)期間,使用傳統(tǒng)方法對(duì)池塘中藻類進(jìn)行鑒別計(jì)數(shù)(圖4)。結(jié)果表明,單養(yǎng)池塘的浮游植物為3.00×107個(gè)/L,混養(yǎng)池塘可達(dá)1.04×108個(gè)/L,混養(yǎng)池塘中的浮游植物數(shù)目高于單養(yǎng)池塘,兩池塘浮游植物密度存在顯著性差異(P<0.05)。

    2.3? ?浮游動(dòng)物和底棲動(dòng)物豐度

    兩類池塘浮游動(dòng)物和底棲動(dòng)物組成及相對(duì)豐度見圖5。單養(yǎng)池塘浮游動(dòng)物和底棲動(dòng)物共10門,混養(yǎng)池塘9門,兩類池塘均以袋形動(dòng)物門(Aschelminthes)、節(jié)肢動(dòng)物門(Arthropoda)、軟體動(dòng)物門(Mollusca)、脊索動(dòng)物門(Chordata)、刺胞動(dòng)物門(Cnidaria)為主,棘皮動(dòng)物門(Echinodermata)、紐形動(dòng)物門(Nemertea)、海綿動(dòng)物門(Porifera)等檢測(cè)量均不超過1%。目分類水平上,兩類池塘中浮游動(dòng)物均以游泳輪蟲目(Ploima)為優(yōu)勢(shì),豐度分別可以達(dá)到21.54%和20.83%,存在劍水蚤目(Cyclopoida)、雙甲目(Diplostraca)、哲水蚤目(Calanoida)。屬分類水平上,兩池塘以多肢輪屬(Polyarthra)、無柄輪屬(Ascomorpha)為優(yōu)勢(shì),存在溫劍水蚤屬(Thermocyclops)、華哲水蚤屬(Sinocalanus)、臂尾輪屬(Brachionus)、龜甲輪屬(Keratella),單養(yǎng)池塘發(fā)現(xiàn)少量偽鏢水蚤屬(Pseudodiaptomus),混養(yǎng)池塘未發(fā)現(xiàn)。

    采取傳統(tǒng)方法對(duì)兩池塘中原生動(dòng)物(Protozoa)、輪蟲(Rotifera)、枝角類(Cladocera)、橈足類(Copepoda)、底棲動(dòng)物的鑒別計(jì)數(shù)結(jié)果見圖6-a。兩池塘均以原生動(dòng)物和輪蟲為優(yōu)勢(shì),存在枝角類、橈足類等甲殼動(dòng)物,輪蟲均以多肢輪屬為優(yōu)勢(shì),其次是無柄輪屬、臂尾輪和少量異尾輪屬(Trichocercidae)和龜甲輪屬,橈足類多為劍水蚤和哲水蚤,單養(yǎng)池塘的浮游動(dòng)物數(shù)目高于混養(yǎng)池塘,兩池塘橈足類密度差異顯著(P<0.05)。

    使用傳統(tǒng)方法對(duì)底棲生物的計(jì)數(shù)結(jié)果見圖6-b。單養(yǎng)池塘發(fā)現(xiàn)了腹足動(dòng)物(Gastropoda)、寡毛類(Oligochaeta)和搖蚊幼蟲(Chironomidae),而混養(yǎng)池塘觀察到腹足動(dòng)物(螺類)和搖蚊幼蟲。混養(yǎng)池塘的螺類多于單養(yǎng)池塘。

    2.4? ?菌群組成

    2.4.1? ?菌群豐度? ?使用環(huán)境DNA技術(shù)檢測(cè)池塘中細(xì)菌的相對(duì)豐度,結(jié)果如圖7所示。

    單養(yǎng)池塘水體中微生物共計(jì)15門27綱38目,優(yōu)勢(shì)菌群為厚壁菌門(Firmicutes)81.86%,其次是變形菌門(Proteobacteria)6.40%、放線菌門(Actinobacteria)5.40%、擬桿菌門(Bacteroidetes)2.20%。綱分類水平以芽孢桿菌綱(Bacilli)81.76%占優(yōu),其次是γ-變形菌綱(Gammaproteobacteria)4.33%,放線菌綱(Actinobacteria)4.26%。屬分類水平,以微小桿菌屬(Exiguobacterium)77.40%占優(yōu)。

    單養(yǎng)池塘底泥中微生物共檢測(cè)出46門,以變形菌門44.61%為優(yōu)勢(shì)菌,其次為厚壁菌門23.42%、綠彎菌門(Chloroflexi)11.22%、擬桿菌門5.67%、酸桿菌門(Acidobacteria)24.97%、放線菌門2.92%。綱分類水平以γ-變形菌綱34.64%占優(yōu),其次是芽孢桿菌綱18.41%,厭氧繩菌綱(Anaerolineae)10.36%,δ-變形菌綱(Deltaproteobacteria)8.06%,擬桿菌綱(Bacteroidia)4.69%,梭狀芽胞桿菌(Clostridia)4.69%。屬分類水平以Paenisporosarcina(13.90%)占優(yōu)。

    混養(yǎng)池塘水體中微生物共16門34綱46目,優(yōu)勢(shì)菌群為變形菌門44.15%,其次是厚壁菌門35.08%、放線菌門8.66%、擬桿菌門2.68%、疣微菌門(Verrucomicrobia)2.20%。綱分類水平以γ-變形菌綱41.02%占優(yōu),其次是芽孢桿菌綱34.89%,放線菌綱6.96%,產(chǎn)氧光細(xì)菌綱(Oxyphotobacteria)5.75%。屬分類水平,以不動(dòng)桿菌屬(Acinetobacter)(35.52%)占優(yōu)。

    底泥中微生物共37門,變形菌門54.56%為優(yōu)勢(shì)菌,還存在厚壁菌門16.11%、放線菌門8.20%、綠彎菌門6.77%、擬桿菌門4.75%、酸桿菌門4.68%。綱分類水平,以γ-變形菌綱45.25%占優(yōu),其次是芽孢桿菌綱15.09%,放線菌綱6.80%,α-變形菌綱(Alphaproteobacteria)4.86%,δ-變形菌綱4.41%,擬桿菌綱4.24%。屬分類水平以假單胞菌屬(Pseudomonas)(16.51%)占優(yōu)。

    2.4.2? ?α-多樣性指數(shù)? ?以97%序列相似度為閾值生成操作性分類單元(Operational Taxonomic Units, OTUs)后,計(jì)算得出α-多樣性指數(shù)如表1所示。水體中Chao1指數(shù)和Shannon指數(shù)均為混養(yǎng)池塘>單養(yǎng)池塘,而底泥中Chao1指數(shù)和Shannon指數(shù)均為單養(yǎng)池塘>混養(yǎng)池塘。無論是Chao1指數(shù)還是Shannon指數(shù)均為底泥樣本>水體樣本,這表明兩池塘底泥中菌群的豐富度和多樣性均高于水體,混養(yǎng)池塘水體中菌群的豐富度和多樣性高于單養(yǎng)池塘,而單養(yǎng)池塘底泥中菌群的豐富度和多樣性高于混養(yǎng)池塘。

    3? ?討論

    本研究結(jié)合環(huán)境DNA和傳統(tǒng)鑒別方法對(duì)單養(yǎng)和混養(yǎng)兩種草魚池塘養(yǎng)殖模式的生物組成進(jìn)行了鑒別和分析,從生態(tài)角度對(duì)兩種養(yǎng)殖模式進(jìn)行了比較,在水質(zhì)、浮游生物、菌群結(jié)構(gòu)等方面具體分析了兩類草魚池塘的差異,為構(gòu)建物質(zhì)能量模型、進(jìn)一步分析養(yǎng)殖生態(tài)系統(tǒng)的結(jié)構(gòu)特征提供了依據(jù)。使用傳統(tǒng)方法對(duì)池塘物種鑒別較為準(zhǔn)確,但采樣操作復(fù)雜、觀測(cè)難度大,且會(huì)對(duì)生境造成破壞;環(huán)境DNA方法的效率和靈敏度高,能夠監(jiān)測(cè)環(huán)境中隱蔽和難以識(shí)別的物種,但易受環(huán)境中其他因素污染。環(huán)境DNA結(jié)合傳統(tǒng)方法進(jìn)行物種鑒別,不僅能比較兩種方法的差異,二者也相互補(bǔ)充,提高了對(duì)養(yǎng)殖系統(tǒng)內(nèi)生物分析的準(zhǔn)確性。

    3.1? ?環(huán)境DNA鑒別方法存在一定的局限性

    本研究結(jié)果中,環(huán)境DNA方法和傳統(tǒng)鑒別方法得出的結(jié)論有部分出入,如二者對(duì)藻類的結(jié)果有出入,環(huán)境DNA方法監(jiān)測(cè)到池塘中不存在的物種等。這種差異可能由于鳥類、昆蟲的??考敖邓畮淼腄NA等對(duì)池塘DNA造成了一定污染;底泥沉積中未被降解的DNA也可能會(huì)使不存在的物種被檢測(cè)到,從而導(dǎo)致環(huán)境DNA監(jiān)測(cè)結(jié)果出現(xiàn)一定偏差;同時(shí),池塘水體過小且流動(dòng)能力弱,使DNA無法均勻分散而在某一區(qū)域堆積,也可能造成調(diào)查的偏差;而環(huán)境DNA的靈敏性也使得其在擴(kuò)增測(cè)序及對(duì)比分析的過程中易產(chǎn)生偏差。因此,對(duì)于目標(biāo)物種引物以及數(shù)據(jù)庫的選擇需要慎重,使用傳統(tǒng)監(jiān)測(cè)手段對(duì)環(huán)境DNA技術(shù)進(jìn)行補(bǔ)充驗(yàn)證是必要的。

    3.2? ?混養(yǎng)池塘高密度藻類有助于吸收水體氮磷

    本試驗(yàn)結(jié)果表明,單養(yǎng)池塘水體和底泥中的氮磷含量均高于混養(yǎng)池塘,這與其他學(xué)者的研究結(jié)論相似(劉朋,2012;孫云飛,2013)。養(yǎng)殖水體中的氮磷輸入以飼料投放為主,過度投放的飼料未被養(yǎng)殖生物完全攝食,會(huì)增加水體中有機(jī)質(zhì)的含量。單養(yǎng)池塘本身生物組成簡(jiǎn)單、食物鏈短,也會(huì)造成營(yíng)養(yǎng)物質(zhì)的浪費(fèi),水質(zhì)容易受到影響;而混養(yǎng)池塘中的鰱、鳙等通過濾食作用過濾掉水中的殘餌,減少飼料中有機(jī)物的溶解釋放,對(duì)水質(zhì)凈化起到了一定作用。氮磷是藻類增殖所必需的營(yíng)養(yǎng)鹽(孫靈毅等,2003),養(yǎng)殖水體中的溶解性氮磷在藻類生長(zhǎng)過程中被大量消耗吸收,能夠降低池塘中氮磷含量?;祓B(yǎng)池塘的藻類密度遠(yuǎn)高于單養(yǎng)池塘,因此對(duì)水中氮磷的吸收利用更好(張雷鳴等,2020),對(duì)養(yǎng)殖水體能起到一定的凈化效果。

    3.3? ?池塘混養(yǎng)可減少藍(lán)綠藻并使浮游動(dòng)物小型化

    試驗(yàn)發(fā)現(xiàn),單養(yǎng)池塘和混養(yǎng)池塘的浮游植物主要為綠藻門、藍(lán)藻門和隱藻門,硅藻、甲藻和少量金藻也在兩類池塘觀察到,單養(yǎng)池塘以衣藻(Chlamydomonas)、隱藻屬(Cryptomonas)、海鏈藻(Thalassiosira)為優(yōu)勢(shì),混養(yǎng)池塘以隱藻屬為優(yōu)勢(shì)。這與李瑞嬌(2014)的研究結(jié)果相似,即草魚池塘養(yǎng)殖期間,浮游植物優(yōu)勢(shì)種集中在藍(lán)藻門、隱藻門和綠藻門;Turker等(2003)研究發(fā)現(xiàn),鰱、鳙能夠通過濾食作用顯著減少藍(lán)綠藻的數(shù)量,這可能是混養(yǎng)池塘以隱藻屬為優(yōu)勢(shì)種的原因之一。大量研究表明,適量放養(yǎng)鰱、鳙等濾食性魚類,可促進(jìn)水體養(yǎng)分循環(huán),有助于浮游植物密度和生物量的增加(趙玉寶,1993;谷孝鴻和劉桂英,1996;李瑞嬌,2014),這與本試驗(yàn)混養(yǎng)池塘的藻類密度高于單養(yǎng)池塘的結(jié)果相符。

    單養(yǎng)池塘和混養(yǎng)池塘中浮游動(dòng)物、底棲生物組成相似,均以多肢輪屬、無柄輪屬、溫劍水蚤屬、華哲水蚤屬、異尾輪科等為優(yōu)勢(shì)。輪蟲和原生動(dòng)物均具有孤雌生殖方式,生長(zhǎng)發(fā)育迅速,能在短期內(nèi)達(dá)到較高的豐度(陳佳林等,2021),因此在池塘生態(tài)系統(tǒng)的浮游動(dòng)物中占優(yōu)勢(shì)。兩類池塘中輪蟲較多且大都為廣溫性種類,水體均傾向于中營(yíng)養(yǎng)化狀態(tài)(葉曉彤等,2020)。單養(yǎng)池塘的浮游動(dòng)物數(shù)目明顯高于混養(yǎng)池塘,與陳雷等(2009)研究結(jié)果不符,可能由于混養(yǎng)池塘中的甲藻對(duì)浮游動(dòng)物產(chǎn)生一定的危害性;而混養(yǎng)池塘放養(yǎng)的濾食性鰱、鳙對(duì)浮游動(dòng)物的攝食壓力增大(王玉彬,2007),也會(huì)導(dǎo)致浮游動(dòng)物數(shù)量減少并趨向小型化(陳光榮等,2008)。

    3.4? ?池塘水體和底泥優(yōu)勢(shì)菌群及其生態(tài)環(huán)境功效

    兩類池塘菌群均以變形菌門和厚壁菌門為優(yōu)勢(shì)類群,但菌群的結(jié)構(gòu)和相對(duì)豐度有所差異。單養(yǎng)池塘水體以厚壁菌門(81.86%)、芽孢桿菌綱(81.76%)、微小桿菌屬(77.40%)為優(yōu)勢(shì)菌,混養(yǎng)池塘水體以變形菌門(44.15%)、γ-變形菌綱(41.02%)、不動(dòng)桿菌屬(35.52%)為優(yōu)勢(shì)菌。兩類池塘底泥均以變形菌門、γ-變形菌綱和厚壁菌門、芽孢桿菌綱為優(yōu)勢(shì)亞群,以Paenisporosarcina和假單胞菌屬為優(yōu)勢(shì)菌群。底泥中的菌群比水體豐富,存在大量的厭氧型菌群,與氮循環(huán)關(guān)系密切的Nitrospirae和Nitrospinae在底泥中的含量遠(yuǎn)高于水體。池塘中的微小桿菌屬在降解有機(jī)物、去除重金屬、脫蛋白以及處理有機(jī)高鹽廢水等方面發(fā)揮著重要的作用(徐蓮等,2012;張瑩等,2013)。γ-變形菌作為海洋生態(tài)系統(tǒng)中普遍存在的細(xì)菌,在有機(jī)質(zhì)含量較高的水域中多以硫氧化細(xì)菌及固氮細(xì)菌為優(yōu)勢(shì)菌群(Ivanova et al, 2003;Payne et al, 2006),能夠在厭氧條件下進(jìn)行光合作用,對(duì)環(huán)境中的碳、硫循環(huán)起著重要作用(白潔等,2009;Bakunina et al, 2020)。本試驗(yàn)還發(fā)現(xiàn),混養(yǎng)池塘水體的菌群多樣性和豐富度均高于單養(yǎng)池塘,而在底泥中呈現(xiàn)相反的趨勢(shì),這與田相利等(2012)對(duì)草魚養(yǎng)殖初期菌群多樣性的研究結(jié)論不一致,與養(yǎng)殖末期菌群多樣性的結(jié)論一致。可見在養(yǎng)殖過程中菌群的演替已經(jīng)開始,試驗(yàn)結(jié)果同時(shí)也受到溫度、飼料投放、養(yǎng)殖水體等外界條件的影響。

    綜上所述,合理混養(yǎng)可以充分利用飼料(劉朋,2012)、調(diào)節(jié)養(yǎng)殖系統(tǒng)中的菌落結(jié)構(gòu)(李秋芬等,2002)、改變浮游動(dòng)植物的豐富度(張雷鳴等,2020)。作為我國目前應(yīng)用較廣的兩種草魚池塘養(yǎng)殖模式,單養(yǎng)和混養(yǎng)草魚在浮游生物、菌群組成與豐富度上存在著一定的差異性。草魚養(yǎng)殖池塘中合理混養(yǎng)鰱、鳙等濾食性魚類,能夠有效地降低水體中氮磷的含量,增加浮游植物的數(shù)目,改變菌群的結(jié)構(gòu)和豐富度。草魚混養(yǎng)對(duì)改善池塘的水質(zhì)、浮游生物以及菌群的組成結(jié)構(gòu)等方面存在一定的正面作用,而兩種養(yǎng)殖方式在生態(tài)系統(tǒng)結(jié)構(gòu)、能量流動(dòng)以及營(yíng)養(yǎng)動(dòng)力學(xué)等方面的具體特征和差異,以及能否通過改善魚種的搭配提高生態(tài)效益、經(jīng)濟(jì)效益還有待進(jìn)一步探究。

    參考文獻(xiàn)

    巴日斯,2018. 環(huán)境DNA及其在水生生態(tài)系統(tǒng)保護(hù)中的應(yīng)用研究[J]. 資源節(jié)約與環(huán)保, (3):17,25.

    白潔,李海艷,趙陽國,2009. 黃海北部不同站位海洋細(xì)菌群落分布特征[J]. 微生物學(xué)報(bào), 49(3):343-350.

    常華,2009. 廣東草魚養(yǎng)殖概況及現(xiàn)狀分析[J]. 海洋與漁業(yè):水產(chǎn)前沿, (6):50-52.

    陳光榮,鐘萍,張修峰,等,2008. 惠州西湖浮游動(dòng)物及其與水質(zhì)的關(guān)系[J]. 湖泊科學(xué), 20(3):351-356.

    陳佳林,余海軍,王茜,2021. 大清河流域浮游動(dòng)物的群落特征研究[J]. 水生態(tài)學(xué)雜志, 42(3):72-78.

    陳雷,徐兆禮,姚煒民,等,2009. 甌江口春季營(yíng)養(yǎng)鹽、浮游植物和浮游動(dòng)物的分布[J]. 生態(tài)學(xué)報(bào), (3):1571-1577.

    陳小江,汪祈超,2016. 湖北通??诘貐^(qū)草魚養(yǎng)殖模式存在的問題及對(duì)策[J]. 漁業(yè)致富指南, 4(20):34-35.

    董雙林,2014. 中國綜合水產(chǎn)養(yǎng)殖的生態(tài)學(xué)基礎(chǔ)[M]. 北京:科學(xué)出版社.

    高攀,蔣明,趙宇江,等,2009. 主養(yǎng)草魚池塘水質(zhì)指標(biāo)的變化規(guī)律和氮磷收支[J]. 云南農(nóng)業(yè)大學(xué)學(xué)報(bào), 24(1):71-77.

    谷孝鴻,劉桂英,1996. 濾食性鰱鳙魚對(duì)池塘浮游生物的影響[J]. 農(nóng)村生態(tài)環(huán)境, 12(1):6-10,41.

    國家環(huán)境保護(hù)局標(biāo)準(zhǔn)處,1989. 水質(zhì) 總磷的測(cè)定: 鉬酸銨分光光度法: GB/T 11893-1989[S].

    國家環(huán)境保護(hù)局規(guī)劃標(biāo)準(zhǔn)處,1987. 水質(zhì) 亞硝酸鹽氮的測(cè)定: 分光光度法: GB 7493-87[S].

    國家環(huán)境保護(hù)局規(guī)劃標(biāo)準(zhǔn)處,1989. 水質(zhì) 總氮的測(cè)定: 堿性過硫酸鉀消解紫外分光光度法: GB 11894-89[S].

    國家環(huán)境保護(hù)總局,2002. 水和廢水監(jiān)測(cè)分析方法[M].4版. 北京:中國環(huán)境科學(xué)出版社.

    國家環(huán)境保護(hù)總局科技標(biāo)準(zhǔn)司,2007. 水質(zhì) 硝酸鹽氮的測(cè)定: 紫外分光光度法(試行): HJ/T346-2007[S]. 北京: 中國環(huán)境科學(xué)出版社.

    國家環(huán)境保護(hù)總局科技標(biāo)準(zhǔn)司,2009.水質(zhì) 氨氮的測(cè)定: 納氏試劑分光光度法: HJ 535-2009[S]. 北京: 中國環(huán)境科學(xué)出版社.

    國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系,2021. 草魚產(chǎn)業(yè)發(fā)展報(bào)告[J]. 中國水產(chǎn), (2):26-37.

    李萌,尉婷婷,史博洋,等,2019. 環(huán)境DNA技術(shù)在淡水底棲大型無脊椎動(dòng)物多樣性監(jiān)測(cè)中的應(yīng)用[J].生物多樣性, 27(5):480-490.

    李秋芬,陳碧鵑,曲克明,等,2002. 魚蝦混養(yǎng)生態(tài)系中細(xì)菌動(dòng)態(tài)變化的研究[J]. 應(yīng)用生態(tài)學(xué)報(bào), 13(6):731-734.

    李瑞嬌,2014. 草魚池塘浮游植物群落結(jié)構(gòu)及初級(jí)生產(chǎn)力、能量轉(zhuǎn)換的初步研究[D]. 武漢:華中農(nóng)業(yè)大學(xué).

    劉玲儀,金有坤,1994. 養(yǎng)殖魚類細(xì)菌性敗血癥暴發(fā)期間魚池水化學(xué)環(huán)境的初步調(diào)查[J]. 水產(chǎn)學(xué)報(bào), 18(2):153-156.

    劉朋,2012. 草魚不同混養(yǎng)系統(tǒng)有機(jī)碳、氮和磷收支及變動(dòng)的初步研究[D]. 青島:中國海洋大學(xué).

    梅順,任青松,2017. 南陵縣草魚實(shí)用養(yǎng)殖模式調(diào)查研究[J]. 現(xiàn)代農(nóng)業(yè)科技,(15):236,240.

    舒璐,林佳艷,徐源,等,2020. 基于環(huán)境DNA宏條形碼的洱海魚類多樣性研究[J]. 水生生物學(xué)報(bào), 44(5):1080-1086.

    孫靈毅,王力勇,徐惠章,2003. 單細(xì)胞藻類對(duì)營(yíng)養(yǎng)鹽類的吸收[J]. 水產(chǎn)科技情報(bào), (3):119-120.

    孫云飛,2013. 草魚(Ctenopharyngodon idellus)混養(yǎng)系統(tǒng)氮磷收支和池塘水質(zhì)與底質(zhì)的比較研究[D]. 青島:中國海洋大學(xué).

    田相利,鄭瑤瑤,柳炳俊,等,2012. 草魚混養(yǎng)系統(tǒng)細(xì)菌數(shù)量變動(dòng)和群落功能多樣性研究[J]. 中國海洋大學(xué)學(xué)報(bào), 42(11):19-27.

    王玉彬,2007. 茅蓮湖水產(chǎn)養(yǎng)殖池塘中浮游生物的研究[D]. 南昌:南昌大學(xué).

    徐蓮,孫紀(jì)全,吳曉磊,等,2012. 菌株Rhodococcus sp. Chr-9和Exiguobacterium sp. Chr-43的除鉻(VI)特性[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 18 (6):971-977.

    許鄭超,鄧燕飛,周軍,等,2021. 中華絨螯蟹養(yǎng)殖池塘菌群結(jié)構(gòu)和功能特征[J]. 生態(tài)學(xué)雜志, 40(7):2223-2233.

    楊帆,2016. 湖北草魚養(yǎng)殖模式總結(jié)之新洲篇[J]. 當(dāng)代水產(chǎn), 41(2):94-95.

    葉曉彤,梁迪文,王慶,等,2020. 洞庭湖流域常德柳葉湖及其鄰近水體輪蟲群落結(jié)構(gòu)變化及其對(duì)環(huán)境因子的響應(yīng)[J]. 湖泊科學(xué), 32(4):1126-1139.

    占家智,羊茜,2015. 高效養(yǎng)淡水魚[M]. 北京:機(jī)械工業(yè)出版社.

    張雷鳴,原居林,倪蒙,等,2020. 兩種池塘養(yǎng)殖模式水質(zhì)因子和浮游植物群落比較分析[J]. 浙江農(nóng)業(yè)學(xué)報(bào), 32(2):317-326.

    張麗娟,徐杉,趙崢,等,2021. 環(huán)境DNA宏條形碼監(jiān)測(cè)湖泊真核浮游植物的精準(zhǔn)性[J]. 環(huán)境科學(xué), 42(2):796-807.

    張瑩,石萍,馬炯,2013. 微小桿菌Exiguobacterium spp.及其環(huán)境應(yīng)用研究進(jìn)展[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào), 19(5):898-904.

    趙玉寶,1993. 鯉魚種和鰱鳙對(duì)池塘浮游生物的影響[J]. 生態(tài)學(xué)報(bào),13(4):348-355.

    周波,2015. 基于EwE模型的草魚綜合養(yǎng)殖池塘生態(tài)系統(tǒng)研究[D]. 青島:中國海洋大學(xué).

    Bakunina I I, Shevchenko L S, Nedashkovskaia O I, et al, 2000. Screening of marine bacteria for fucoidan hydrolases[J]. Mikrobiologiia, 69(3):370.

    Harper L R, Buxton A S, Rees H C, et al, 2019. Prospects and challenges of environmental DNA (eDNA) monitoring in freshwater ponds[J]. Hydrobiologia, 826(1):25-41.

    Ivanova E P, Sawabe T, Zhukova N V, et al, 2003. Occurrence and diversity of mesophilic Shewanella strains isolated from the North-West Pacific Ocean[J]. Systematic & Applied Microbiology, 26(2):293-301.

    Li Z, Wang G, Yu E, et al, 2019. Artificial substrata increase pond farming density of grass carp (Ctenopharyngodon idella) by increasing the bacteria that participate in nitrogen and phosphorus cycles in pond water. PeerJ, 7:e7906.

    Niu S, Zhang K, Li Z, et al, 2023. Analysis of the structure and function of microbial community in late-stage of grass carp (Ctenopharyngodon idella) farming ponds[J]. Aquaculture Reports, 30:101556.

    Payne M S, Hall M R, Bannister R, et al, 2006. Microbial diversity within the water column of a larval rearing system for the ornate rock lobster (Panulirus ornatus)[J]. Aquaculture, 258(1/4):80-90.

    Singer D, Seppey C, Lentendu G, et al, 2021. Protist taxonomic and functional diversity in soil, freshwater and marine ecosystems[J]. Environment International, 146:106262.

    Taberlet P, Coissac E, Hajibabaei M, et al, 2012. Environmental DNA[J]. Molecular Ecology, 21(8):1789-1793.

    Turker H, Eversole A G, Brune D E, 2003. Comparative Nile tilapia and silver carp filtration rates of Partitioned Aquaculture System phytoplankton[J]. Aquaculture, 220(1/4):449-457.

    (責(zé)任編輯? ?張俊友? ?熊美華)

    Water Quality and Bio-composition in a Monoculture

    and Polyculture Grass Carp Pond

    XIAO Shu‐wen1,2, LIU Xing‐guo2, LU Shi‐min2, ZHAO Yu‐xi1,2, GU Zhao‐jun2, ZHOU Run‐feng1,2

    (1. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai? ?201306, P.R. China;

    2. Fishery Machinery and Instrument Research Institute, Chinese Academy of Fishery Sciences,

    Shanghai? ?200092, P.R. China)

    Abstract:In this study, we analyzed the differences in water quality, plankton, benthos and bacterial community structure between monoculture and polyculture grass carp ponds (80:20). Water quality analysis, environmental DNA and traditional identification methods were used, and the water quality and biological composition of the two culturing methods were compared. We aimed to provide evidence for evaluating grass carp culture modes and construction of highly efficient grass carp culture systems. This study also explored the existing problems of using environmental DNA to monitor culture pond ecosystems. Results show that water quality in the polyculture pond was better than in the monoculture pond. Concentrations of TN, NO3--N, NH4+-N and NO2--N in the polyculture pond were lower, respectively, by 10.15%, 3.78%, 5.07% and 80.18%, and the concentrations of TP and SRP were 27.14% and 56.26% lower than in the monoculture pond. Chlorophyta, Cyanophyta and Cryptophyta were dominated in both ponds. The density of phytoplankton in the monoculture pond was 3.00×107 ind/L, lower than in polyculture pond (1.04×108 ind/L). Rotifers and protozoa were the dominant zooplankton groups in both ponds, with lower densities of Cladophora and Copepoda. The zooplankton density in the monoculture pond was higher than in the polyculture pond. In terms of benthos, Gastropoda, Oligochaeta and Chironomidae were observed in the monoculture pond, while only Gastropoda and Chironomidae were found in the polyculture pond. In terms of the microbial composition, Firmicutes was the dominant group in the monoculture pond, Proteobacteri was the dominant group in the polyculture pond and Proteobacteria was the dominant group in the sediments of both ponds. In conclusion, grass carp polyculture can improve the water quality of aquaculture ponds, increase phytoplankton richness and change the structure of aquaculture pond flora. This study provides a basis for optimizing grass carp culture, improving water quality and constructing an efficient culture system.

    Key words:grass carp pond; culture mode; water quality parameters; biological composition; environmental DNA

    收稿日期:2022-01-10? ? ? 修回日期:2022-10-21

    基金項(xiàng)目:國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)項(xiàng)目(CARS-46);國家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2019YFD0900300)。

    作者簡(jiǎn)介:肖述文,1998年生,女,碩士研究生,研究方向?yàn)樗蛏鷳B(tài)修復(fù)。E-mail:949960731@qq.com

    通信作者:劉興國。E-mail:liuxingguo@fmiri.ac.cn

    猜你喜歡
    養(yǎng)殖模式
    試析我國生豬養(yǎng)殖模式與疫病防控
    家禽規(guī)范養(yǎng)殖的幾點(diǎn)看法
    淺談不同養(yǎng)殖模式及存欄規(guī)模牛場(chǎng)奶牛的生產(chǎn)性能
    3種不同養(yǎng)殖模式對(duì)精養(yǎng)水體尼羅羅非魚生長(zhǎng)性能及養(yǎng)殖水質(zhì)的影響
    長(zhǎng)江江灘林間不同養(yǎng)殖模式草雞蛋品質(zhì)研究
    英山茶鄉(xiāng)雞及其養(yǎng)殖模式探討
    不同養(yǎng)殖模式及光照度對(duì)蛋雞產(chǎn)蛋率的影響
    絲尾鳠與大宗淡水魚網(wǎng)箱養(yǎng)殖模式比較研究
    吉富羅非魚單養(yǎng)與混養(yǎng)模式養(yǎng)殖效果的比較
    湖區(qū)波爾山羊生態(tài)型健康養(yǎng)殖模式
    一级av片app| 两个人的视频大全免费| 亚洲美女视频黄频| 午夜福利视频精品| 青春草国产在线视频| 伊人久久国产一区二区| 夫妻性生交免费视频一级片| 搡女人真爽免费视频火全软件| 免费大片18禁| 乱系列少妇在线播放| 色尼玛亚洲综合影院| 高清视频免费观看一区二区 | 色综合站精品国产| 亚洲最大成人av| 美女xxoo啪啪120秒动态图| 高清av免费在线| 亚洲av成人av| 天堂√8在线中文| 日韩强制内射视频| 综合色丁香网| 嫩草影院精品99| 久久精品久久精品一区二区三区| 国产精品久久久久久av不卡| 18禁在线无遮挡免费观看视频| 日日摸夜夜添夜夜爱| 婷婷色综合www| 亚洲精品第二区| 欧美日韩国产mv在线观看视频 | 午夜精品一区二区三区免费看| 久久久a久久爽久久v久久| 亚洲av二区三区四区| 日产精品乱码卡一卡2卡三| 少妇的逼好多水| 婷婷色综合www| 日韩欧美国产在线观看| 成人亚洲欧美一区二区av| 看非洲黑人一级黄片| 一级毛片久久久久久久久女| 蜜桃亚洲精品一区二区三区| 少妇高潮的动态图| 亚洲怡红院男人天堂| 视频中文字幕在线观看| 18禁动态无遮挡网站| 亚洲性久久影院| 免费高清在线观看视频在线观看| 男女国产视频网站| 国产亚洲一区二区精品| 欧美激情国产日韩精品一区| 精品久久久久久久久亚洲| 亚洲av一区综合| 国产亚洲一区二区精品| 国产精品一区二区三区四区免费观看| 国产伦精品一区二区三区视频9| 久久韩国三级中文字幕| 午夜福利成人在线免费观看| 老司机影院毛片| 久久这里只有精品中国| 毛片一级片免费看久久久久| av黄色大香蕉| 色网站视频免费| 少妇被粗大猛烈的视频| 丝袜美腿在线中文| 精品亚洲乱码少妇综合久久| 国产男女超爽视频在线观看| 一级毛片 在线播放| 国产精品久久视频播放| 99热这里只有是精品50| 日韩av不卡免费在线播放| 夫妻午夜视频| 欧美丝袜亚洲另类| 人妻制服诱惑在线中文字幕| 亚洲av男天堂| 97热精品久久久久久| 直男gayav资源| 波野结衣二区三区在线| 亚洲熟妇中文字幕五十中出| 2021少妇久久久久久久久久久| 亚洲在线观看片| 国产成人a∨麻豆精品| 乱系列少妇在线播放| 国产av不卡久久| 在现免费观看毛片| 国产高潮美女av| 欧美成人a在线观看| 人妻制服诱惑在线中文字幕| 日韩欧美三级三区| 国产精品麻豆人妻色哟哟久久 | 国内揄拍国产精品人妻在线| 最近的中文字幕免费完整| 三级经典国产精品| 一夜夜www| 日韩,欧美,国产一区二区三区| 久久久久久久午夜电影| 搞女人的毛片| 一区二区三区免费毛片| 亚洲欧美精品自产自拍| 日韩 亚洲 欧美在线| 国产爱豆传媒在线观看| 五月天丁香电影| 日本午夜av视频| 欧美+日韩+精品| 18禁在线播放成人免费| 日韩成人伦理影院| 久久精品国产自在天天线| 女人十人毛片免费观看3o分钟| 国产伦精品一区二区三区四那| 成人午夜高清在线视频| 少妇高潮的动态图| 最近中文字幕高清免费大全6| 久久久精品94久久精品| 91久久精品国产一区二区三区| 观看美女的网站| 国产精品久久视频播放| 麻豆成人午夜福利视频| 韩国高清视频一区二区三区| 51国产日韩欧美| 男女国产视频网站| 狂野欧美激情性xxxx在线观看| 国产一区亚洲一区在线观看| av在线蜜桃| 少妇的逼水好多| 久热久热在线精品观看| 午夜老司机福利剧场| 亚洲av一区综合| 色吧在线观看| 精华霜和精华液先用哪个| 插逼视频在线观看| 国产午夜精品一二区理论片| 欧美成人精品欧美一级黄| 亚洲va在线va天堂va国产| 一级片'在线观看视频| av国产免费在线观看| 成年版毛片免费区| 天堂av国产一区二区熟女人妻| 日韩欧美精品v在线| 午夜久久久久精精品| av线在线观看网站| 成人午夜高清在线视频| 免费不卡的大黄色大毛片视频在线观看 | 国产在视频线精品| 青青草视频在线视频观看| 亚洲av福利一区| 99久久中文字幕三级久久日本| 91aial.com中文字幕在线观看| 国产午夜精品论理片| 内地一区二区视频在线| 老司机影院成人| 国产精品久久久久久av不卡| 亚洲欧洲日产国产| 男人舔女人下体高潮全视频| av在线蜜桃| 国产毛片a区久久久久| 99久久精品一区二区三区| 国产久久久一区二区三区| 99热6这里只有精品| 欧美日韩国产mv在线观看视频 | 一级毛片我不卡| 国精品久久久久久国模美| 久久人人爽人人爽人人片va| 自拍偷自拍亚洲精品老妇| 午夜福利网站1000一区二区三区| 永久网站在线| 国产熟女欧美一区二区| 国产成人freesex在线| 欧美97在线视频| 中文字幕免费在线视频6| 中文乱码字字幕精品一区二区三区 | 一二三四中文在线观看免费高清| 校园人妻丝袜中文字幕| 夫妻午夜视频| 国产精品国产三级国产av玫瑰| 亚洲精华国产精华液的使用体验| 欧美激情在线99| 国产一区亚洲一区在线观看| h日本视频在线播放| 伊人久久国产一区二区| 丰满少妇做爰视频| 噜噜噜噜噜久久久久久91| 热99在线观看视频| 久久99蜜桃精品久久| 成人毛片60女人毛片免费| 午夜福利网站1000一区二区三区| 午夜福利高清视频| 热99在线观看视频| 免费观看无遮挡的男女| 免费在线观看成人毛片| 日本一本二区三区精品| 亚洲性久久影院| 国产熟女欧美一区二区| 免费播放大片免费观看视频在线观看| 美女xxoo啪啪120秒动态图| 毛片一级片免费看久久久久| 国产亚洲一区二区精品| 日本三级黄在线观看| 亚洲精品一二三| 三级国产精品欧美在线观看| 国产不卡一卡二| 欧美性猛交╳xxx乱大交人| 久久久久精品性色| av线在线观看网站| 成人高潮视频无遮挡免费网站| 日韩av不卡免费在线播放| 99久久中文字幕三级久久日本| 一区二区三区四区激情视频| 欧美日韩国产mv在线观看视频 | 国产精品.久久久| 日韩一本色道免费dvd| 国产黄a三级三级三级人| 亚洲精品影视一区二区三区av| 久久99热6这里只有精品| 美女高潮的动态| 大又大粗又爽又黄少妇毛片口| 看十八女毛片水多多多| 免费观看a级毛片全部| 日韩强制内射视频| 美女cb高潮喷水在线观看| 精品99又大又爽又粗少妇毛片| 日日摸夜夜添夜夜添av毛片| 在线观看av片永久免费下载| 精品国产三级普通话版| 一夜夜www| 超碰97精品在线观看| 国产一区有黄有色的免费视频 | 少妇丰满av| or卡值多少钱| 少妇人妻精品综合一区二区| 老师上课跳d突然被开到最大视频| 国产精品.久久久| 97人妻精品一区二区三区麻豆| 少妇的逼好多水| 国产片特级美女逼逼视频| 极品少妇高潮喷水抽搐| 九草在线视频观看| 免费观看av网站的网址| 婷婷六月久久综合丁香| 成人午夜高清在线视频| 亚洲av电影不卡..在线观看| 久久鲁丝午夜福利片| 久久99蜜桃精品久久| 亚洲欧美成人综合另类久久久| 日韩av在线大香蕉| 高清在线视频一区二区三区| 全区人妻精品视频| 中文字幕av成人在线电影| 国产成人免费观看mmmm| 国产成人免费观看mmmm| 久久久久久久午夜电影| 成人性生交大片免费视频hd| 中文在线观看免费www的网站| 伊人久久精品亚洲午夜| 神马国产精品三级电影在线观看| 国产成人精品婷婷| 禁无遮挡网站| 99热这里只有是精品50| 日本熟妇午夜| 亚洲欧美日韩无卡精品| 欧美性感艳星| 麻豆乱淫一区二区| 少妇的逼水好多| 我的老师免费观看完整版| 亚洲国产av新网站| 久久6这里有精品| 久久草成人影院| 九草在线视频观看| 国产一区二区在线观看日韩| 观看美女的网站| 国产精品一二三区在线看| 日韩在线高清观看一区二区三区| 中文天堂在线官网| av又黄又爽大尺度在线免费看| 美女cb高潮喷水在线观看| 欧美一级a爱片免费观看看| 国产熟女欧美一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品综合久久久久久久免费| 亚洲自偷自拍三级| 美女cb高潮喷水在线观看| 久久久久久久亚洲中文字幕| 久久久久久久久久人人人人人人| 少妇丰满av| 国产精品一及| or卡值多少钱| 人妻系列 视频| 色5月婷婷丁香| 国产亚洲5aaaaa淫片| 精品久久国产蜜桃| 一级黄片播放器| 国产亚洲av片在线观看秒播厂 | 国产有黄有色有爽视频| 美女国产视频在线观看| 高清欧美精品videossex| 老女人水多毛片| 色综合色国产| 亚洲美女视频黄频| 亚洲欧美清纯卡通| 午夜久久久久精精品| av.在线天堂| 国产一区二区亚洲精品在线观看| 亚洲精品成人av观看孕妇| 国产精品人妻久久久久久| 男女边摸边吃奶| 淫秽高清视频在线观看| 三级国产精品片| 熟妇人妻久久中文字幕3abv| 少妇人妻一区二区三区视频| 少妇猛男粗大的猛烈进出视频 | 日韩欧美国产在线观看| 高清午夜精品一区二区三区| 中文精品一卡2卡3卡4更新| 成年人午夜在线观看视频 | 身体一侧抽搐| 激情 狠狠 欧美| 中文字幕免费在线视频6| 久久鲁丝午夜福利片| 亚洲欧洲国产日韩| 亚洲婷婷狠狠爱综合网| 乱码一卡2卡4卡精品| 两个人的视频大全免费| a级毛片免费高清观看在线播放| 久久99热6这里只有精品| 久久精品国产亚洲av涩爱| 97人妻精品一区二区三区麻豆| 91精品伊人久久大香线蕉| 日本三级黄在线观看| 欧美xxⅹ黑人| 一边亲一边摸免费视频| 免费黄色在线免费观看| 欧美激情国产日韩精品一区| 热99在线观看视频| 亚洲欧美清纯卡通| 中文字幕av在线有码专区| 免费高清在线观看视频在线观看| 国产精品精品国产色婷婷| 午夜福利成人在线免费观看| 欧美日韩综合久久久久久| av又黄又爽大尺度在线免费看| 最近2019中文字幕mv第一页| 亚洲最大成人中文| 久久精品夜夜夜夜夜久久蜜豆| 性色avwww在线观看| 成人一区二区视频在线观看| av天堂中文字幕网| 激情五月婷婷亚洲| 国产一级毛片在线| 夫妻性生交免费视频一级片| 亚洲av免费在线观看| 亚洲av免费高清在线观看| 亚洲,欧美,日韩| av免费观看日本| 亚洲内射少妇av| 狂野欧美白嫩少妇大欣赏| 九草在线视频观看| 国产精品一及| 秋霞在线观看毛片| 免费看光身美女| 日本wwww免费看| 在线 av 中文字幕| 欧美日本视频| 26uuu在线亚洲综合色| 五月天丁香电影| 特级一级黄色大片| 亚洲欧洲国产日韩| 在现免费观看毛片| 成年av动漫网址| 免费观看av网站的网址| 亚洲国产欧美人成| 蜜桃亚洲精品一区二区三区| 美女脱内裤让男人舔精品视频| 麻豆国产97在线/欧美| 女人被狂操c到高潮| 国产有黄有色有爽视频| a级毛片免费高清观看在线播放| 免费黄网站久久成人精品| 国产在线一区二区三区精| 22中文网久久字幕| 国产一区二区亚洲精品在线观看| 蜜桃亚洲精品一区二区三区| 午夜福利在线观看免费完整高清在| 综合色丁香网| 极品少妇高潮喷水抽搐| 国产高潮美女av| 91久久精品国产一区二区成人| 又爽又黄a免费视频| 亚洲精品乱码久久久久久按摩| 欧美一区二区亚洲| 亚洲精品一区蜜桃| 午夜福利网站1000一区二区三区| 免费观看性生交大片5| 天堂影院成人在线观看| 欧美精品一区二区大全| 日本爱情动作片www.在线观看| 久久久a久久爽久久v久久| 久久久午夜欧美精品| 性插视频无遮挡在线免费观看| 午夜视频国产福利| 亚洲美女视频黄频| 在线免费十八禁| 伦理电影大哥的女人| 国产黄a三级三级三级人| 午夜福利成人在线免费观看| 大话2 男鬼变身卡| 亚洲精品国产成人久久av| 精品亚洲乱码少妇综合久久| 日本-黄色视频高清免费观看| or卡值多少钱| 尤物成人国产欧美一区二区三区| 国产精品一区二区性色av| 亚洲精品久久午夜乱码| 极品少妇高潮喷水抽搐| 2021天堂中文幕一二区在线观| 日本一二三区视频观看| 亚洲,欧美,日韩| 99久久中文字幕三级久久日本| 三级毛片av免费| 三级男女做爰猛烈吃奶摸视频| 亚洲精品,欧美精品| 观看美女的网站| 丝袜喷水一区| 日韩一本色道免费dvd| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站在线观看播放| 免费观看a级毛片全部| 一本一本综合久久| 日韩成人伦理影院| 亚洲欧美成人综合另类久久久| 最近的中文字幕免费完整| 国产精品久久久久久久电影| 日本免费a在线| 99久久精品一区二区三区| 在线免费观看的www视频| 啦啦啦啦在线视频资源| 黄色日韩在线| 亚洲自偷自拍三级| 欧美极品一区二区三区四区| 亚洲精品自拍成人| videossex国产| 麻豆成人av视频| 亚洲av男天堂| 青春草国产在线视频| 国产女主播在线喷水免费视频网站 | 亚洲av电影在线观看一区二区三区 | av卡一久久| 欧美日本视频| 久久久久久伊人网av| 最近的中文字幕免费完整| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久 | 水蜜桃什么品种好| 欧美日韩综合久久久久久| 一区二区三区四区激情视频| 日韩电影二区| 国内揄拍国产精品人妻在线| 禁无遮挡网站| 国产一区二区三区综合在线观看 | 有码 亚洲区| 久久精品国产亚洲av天美| 男人舔奶头视频| 性色avwww在线观看| 日韩国内少妇激情av| 精品久久久久久久久av| 中文字幕人妻熟人妻熟丝袜美| 中文欧美无线码| 亚洲最大成人手机在线| 白带黄色成豆腐渣| 青春草亚洲视频在线观看| 日韩视频在线欧美| 国产色爽女视频免费观看| 免费看光身美女| 丰满乱子伦码专区| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 好男人在线观看高清免费视频| 少妇高潮的动态图| 激情五月婷婷亚洲| 如何舔出高潮| 国产精品国产三级专区第一集| 国产精品久久视频播放| 亚洲精品一二三| 色吧在线观看| 亚洲精品成人av观看孕妇| 亚洲欧美精品自产自拍| 午夜福利在线观看吧| 国产白丝娇喘喷水9色精品| 搡老乐熟女国产| 日韩一区二区视频免费看| 国产精品人妻久久久影院| 最近最新中文字幕大全电影3| 99久久精品热视频| 欧美极品一区二区三区四区| 亚洲自拍偷在线| 熟女人妻精品中文字幕| 高清午夜精品一区二区三区| 欧美zozozo另类| 久久久久久久久中文| 七月丁香在线播放| 精品久久久久久成人av| freevideosex欧美| av女优亚洲男人天堂| 精品国产露脸久久av麻豆 | 国产人妻一区二区三区在| 国产白丝娇喘喷水9色精品| 国产精品女同一区二区软件| 老女人水多毛片| 精品国产三级普通话版| 欧美丝袜亚洲另类| 国产精品.久久久| 欧美一区二区亚洲| 一级爰片在线观看| 国产伦精品一区二区三区视频9| 亚洲精品第二区| 男的添女的下面高潮视频| 高清欧美精品videossex| 免费大片黄手机在线观看| 青春草亚洲视频在线观看| 国产激情偷乱视频一区二区| 26uuu在线亚洲综合色| 久久精品久久久久久久性| 婷婷色麻豆天堂久久| 日韩 亚洲 欧美在线| 肉色欧美久久久久久久蜜桃 | 亚洲aⅴ乱码一区二区在线播放| 国产探花极品一区二区| 欧美另类一区| 联通29元200g的流量卡| 久久久久久久久大av| 内地一区二区视频在线| 日本午夜av视频| 18禁动态无遮挡网站| 人人妻人人澡欧美一区二区| 午夜激情欧美在线| 中文字幕av成人在线电影| 精品人妻熟女av久视频| 国产爱豆传媒在线观看| 国产麻豆成人av免费视频| 男女啪啪激烈高潮av片| 亚洲欧美一区二区三区国产| 日韩大片免费观看网站| 人妻少妇偷人精品九色| 日本-黄色视频高清免费观看| 美女cb高潮喷水在线观看| 国产av码专区亚洲av| 在线观看免费高清a一片| 日日摸夜夜添夜夜添av毛片| 我要看日韩黄色一级片| 日韩欧美一区视频在线观看 | 久久久久久久久久黄片| 亚洲国产精品国产精品| 国产精品久久久久久精品电影| 久久久久久久久久久丰满| 男女下面进入的视频免费午夜| 成人亚洲精品av一区二区| 日韩精品有码人妻一区| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 国产精品蜜桃在线观看| 美女大奶头视频| 午夜福利成人在线免费观看| 亚洲欧洲国产日韩| 日日啪夜夜爽| 菩萨蛮人人尽说江南好唐韦庄| 国产乱人偷精品视频| 亚洲伊人久久精品综合| 乱码一卡2卡4卡精品| 国产精品人妻久久久影院| 午夜精品一区二区三区免费看| 亚洲精品色激情综合| 国产日韩欧美在线精品| 少妇的逼好多水| 小蜜桃在线观看免费完整版高清| 久久精品综合一区二区三区| 赤兔流量卡办理| 久久综合国产亚洲精品| 日韩,欧美,国产一区二区三区| 毛片一级片免费看久久久久| 欧美日韩综合久久久久久| 免费观看的影片在线观看| 麻豆久久精品国产亚洲av| 成年人午夜在线观看视频 | 中文精品一卡2卡3卡4更新| 联通29元200g的流量卡| 少妇的逼好多水| 看黄色毛片网站| 18禁在线播放成人免费| 国产午夜精品论理片| 校园人妻丝袜中文字幕| 女人久久www免费人成看片| 午夜老司机福利剧场| 日本黄大片高清| 精品一区二区三区视频在线| 欧美精品一区二区大全| 麻豆乱淫一区二区| 人人妻人人澡欧美一区二区| 成人鲁丝片一二三区免费| 黄片wwwwww| 老女人水多毛片| 麻豆国产97在线/欧美| 亚洲精品一二三| 久久精品久久精品一区二区三区| 久久久久久久久久久免费av| 能在线免费看毛片的网站| 精品一区二区三区人妻视频| 久99久视频精品免费| 婷婷色麻豆天堂久久| 女的被弄到高潮叫床怎么办| 国产精品久久久久久精品电影| 一级爰片在线观看| 成人特级av手机在线观看| 搞女人的毛片| 久久午夜福利片| 麻豆久久精品国产亚洲av| 青春草国产在线视频| 少妇被粗大猛烈的视频| 午夜久久久久精精品| 久久人人爽人人爽人人片va| 国产亚洲精品av在线|