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

    升金湖和菜子湖越冬白額雁腸道寄生線蟲多樣性研究

    2022-09-17 03:21:16李清月龔治忠馮佳慧劉剛
    生態(tài)科學(xué) 2022年6期
    關(guān)鍵詞:菜子金湖線蟲

    李清月, 龔治忠, 馮佳慧, 劉剛

    升金湖和菜子湖越冬白額雁腸道寄生線蟲多樣性研究

    李清月, 龔治忠, 馮佳慧, 劉剛*

    安徽醫(yī)科大學(xué)生命科學(xué)學(xué)院, 合肥 230032

    腸道寄生蟲是鳥類中常見的病原體之一, 可以傳播給人類及其他動物, 引起嚴重的人禽共患病。通過采集菜子湖和升金湖60份越冬白額雁糞便樣本, 以腸道寄生線蟲特異引物作為標記基因進行高通量測序, 分析兩個湖泊越冬白額雁的腸道寄生線蟲感染種類、多樣性及其群落組成結(jié)構(gòu)。將3299430條reads歸屬為腸道寄生線蟲門, 共定義3249個腸道寄生線蟲的OTUs, 鑒定到4個目, 91屬和191種; 4個目分別為小桿目(Rhabditida)(62.12 %)、墊刃目(Tylenchida)(32.20%)、疏毛目(Araeolaimida)(3.40%)和單宮目(Monhysterida)(2.27%), 其中小桿目為優(yōu)勢物種。Shannon-wiener指數(shù)和Simpson指數(shù)顯示, 菜子湖越冬白額雁的腸道寄生線蟲alpha多樣性高于升金湖越冬白額雁, 而差異性不顯著, PCA聚類分析表明, 兩個湖泊中越冬白額雁腸道寄生線蟲群落的組成結(jié)構(gòu)相似。系統(tǒng)發(fā)育關(guān)系顯示, 樣品與湖泊匹配度不高, 部分樣品分布不均勻。升金湖和菜子湖越冬白額雁的腸道寄生線蟲群落群落結(jié)構(gòu)相似, 可能是由于食物來源相似、食性相似、棲息地環(huán)境類似等因素引起寄生蟲間的交叉?zhèn)鞑ヅc感染導(dǎo)致。通過本研究,初步掌握升金湖和菜子湖越冬白額雁的腸道寄生線蟲感染現(xiàn)狀, 為進一步探討越冬白額雁的寄生蟲流行病學(xué)研究積累基礎(chǔ)資料。

    腸道寄生線蟲; 白額雁; 多樣性

    0 前言

    鳥類是許多疾病的宿主和傳播媒介, 攜帶多種病毒、細菌和寄生蟲等病原體, 可以傳播給其他動物甚至人類, 引起嚴重的人禽共患病, 由于具有極強的擴散能力和廣泛的分布范圍, 鳥類為疾病的傳播帶來了潛在風(fēng)險[1-3]。腸道寄生蟲是鳥類中最常見、對其威脅最為明顯的病原體之一, 腸道寄生蟲能夠增加鳥類種群的死亡率, 降低出生率, 進而調(diào)節(jié)鳥類的種群數(shù)量[1-2,4-6]。野生鳥類在遷徙過程中接觸不同生境和中間媒介, 促使寄生蟲傳染給新的個體, 并傳播到新的地區(qū), 使寄生蟲的傳播途徑更加多元化和復(fù)雜化。野生鳥類作為寄生蟲的傳播媒介, 為寄生蟲的傳播提供了極大的便利, 為跨區(qū)域人禽共患病的傳播埋下隱患。

    研究表明, 雁形目等長距離遷徙的水鳥攜帶大量球蟲、線蟲、吸蟲和絳蟲等腸道寄生蟲感染, 其中線蟲是鳥類攜帶的腸道寄生蟲中最常見的病原體之一[5,7]。然而, 目前對鳥類腸道寄生線蟲的研究仍然落后于其他脊椎動物腸道寄生蟲和微生物[8-9]。近年來, 對鳥類腸道寄生蟲的研究主要集中在珍稀鳥類如鶴類, 目前有研究發(fā)現(xiàn)沙丘鶴()、灰鶴()和白頭鶴()易感染艾美爾球蟲()、禽蛔蟲()、前盲囊線蟲(Porrocaecum)、四棱線蟲()等腸道寄生蟲[5,10]; 黑鸛()感染呵欠鳥凱塞瑪吸蟲()等寄生蟲[11]; 而丹頂鶴()和東方白鸛()則普遍感染棘口吸蟲()等寄生蟲[12-13]; 也有學(xué)者在野生鴻雁()糞便中檢出隱孢子蟲(), 初步鑒定為貝氏隱孢子蟲()[14], 然而人們對腸道寄生線蟲的流行病學(xué)了解甚少[5,15]。一直以來, 腸道寄生線蟲與野生候鳥以及人禽共患病的等公共衛(wèi)生問題息息相關(guān), 對宿主的健康起著至關(guān)重要的作用, 在過去的二十年中, 腸道寄生線蟲引起的健康問題如導(dǎo)致人類和動物嚴重腹瀉現(xiàn)象, 已受到世界衛(wèi)生組織的高度重視[3,16]。

    菜子湖(30.75°—30.97°N, 117.00°—117.15°E)和升金湖(30.25°—30.50°N, 116.92°—117.25°E)及其密切相聯(lián)的江河共同構(gòu)成的江湖復(fù)合濕地生態(tài)系統(tǒng), 蘊藏著豐富的水生資源, 是長江中下游沿江重要的濕地, 也是東亞-澳大利亞遷徙水鳥的重要遷徙停歇地和越冬地[17-18]。季節(jié)性的水位變化為大量水鳥提供了充足的活動空間和食物, 其中包括雁形目(Anseriformes)中的小天鵝()、白額雁)、鴻雁()、小白額雁()、豆雁()和灰雁()等常見的越冬遷徙水鳥[17-19]。白額雁作為長距離遷徙的水鳥, 是重要的濕地指示物種(IUCN, 2019), 由于其強大的飛行能力和廣泛的地理分布, 在傳播寄生蟲病、皰疹病毒和禽流感等傳染病中扮演重要角色[20-22]。近年來, 隨著鳥類腸道寄生蟲研究的不斷深入, 白額雁腸道寄生線蟲已成為公共衛(wèi)生和鳥類學(xué)的重要研究課題[23]。本研究使用高通量技術(shù)平臺對白額雁腸道線蟲進行測序, 由此鑒定菜子湖和升金湖越冬白額雁腸道寄生線蟲的感染種類, 并結(jié)合生物信息學(xué)技術(shù)分析兩個湖泊之間的腸道線蟲的豐富度及其多樣性, 揭示白額雁腸道線蟲多樣性機制, 確定影響白額雁腸道線蟲種類組成的驅(qū)動因素, 探討越冬白額雁腸道線蟲與外界環(huán)境之間的關(guān)系, 為更深層次研究越冬白額雁的寄生蟲流行病學(xué)奠定了一定的基礎(chǔ)。

    1 材料和方法

    1.1 樣品的采集與保存

    2018年12月, 從菜子湖的梅花大圩村和于家嘴村以及升金湖的楊峨頭村和王壩村分別采集越冬白額雁的糞便樣本共60份。在采集糞便之前, 用雙筒望遠鏡觀察覓食雁群至少1小時, 確定周圍沒有其他鳥群活動, 或與其他集群有明顯的障礙物(如田埂、河溝等)間隔后方能采集, 等到白額雁取食結(jié)束后, 根據(jù)白額雁留下的足跡或覓食坑, 確定為白額雁覓食地點, 采樣時佩戴一次性PE手套, 仔細辨別并采集新鮮糞便樣本, 每采集一份更換一次手套, 采集后迅速將樣品放入無菌20 mL樣品收集管中, 并做好標簽, 放于冰上運到實驗室后置-80℃超低溫冰箱冷凍保存。

    1.2 鳥類物種鑒定、DNA提取和PCR 擴增

    采用德國 QIAamp Fast DNA Stool Mini Kit 糞便微生物基因組 DNA 提取試劑盒來提取糞便樣品中的總DNA, 提取方法嚴格按照試劑盒說明書的具體操作步驟進行, DNA 提取樣品保存在-20℃冰箱中以備后續(xù)實驗使用。利用線粒體16S rRNA特異性引物進行PCR擴增測序, 根據(jù)測序結(jié)果確認鳥類為白額雁。

    從每個樣本中提取50 ng的純化DNA作為擴增模板, 選擇線蟲特異引物區(qū)作為標記基因進行PCR擴增, 引物序列為A‐NF1(5’-GCC TCC CTC GCG CCA TCA GGG TGG TGC ATG GCC GTT CTT AGT T-3’)和 B-18S r2b(5’-GCC TTG CCA GCC CGC TCA GTA CAA AGG GCA GGG ACG TAA T-3’)。PCR反應(yīng)條件: 95 ℃預(yù)變性10 min, 95 ℃變性30 s, 50 ℃退火30 s, 72 ℃延伸30 s, 40個循環(huán)后, 72 ℃延伸10 min, 放于 4 ℃冰箱保存, 隨后利用2%(w/v)的Tris-硼酸-EDTA(TBE)瓊脂糖凝膠來評估擴增質(zhì)量。

    1.3 Illumina Miseq 測序和數(shù)據(jù)分析

    使用MiniElute PCR純化試劑盒(Axygen)純化PCR擴增產(chǎn)物, 并使用Applied Biosystems GeneAmp PCR System 9700進行定量。將PCR產(chǎn)物在Illumina Miseq測序平臺進行高通量測序。

    基于高通量測序得到Fastq的原始數(shù)據(jù), 使用Trimmomatic (version 0.33)軟件對這些原始數(shù)據(jù)進行數(shù)據(jù)質(zhì)量控制, 采用滑動窗口法, 過濾其中質(zhì)量值小于20并且讀長小于50 bp的低質(zhì)量數(shù)據(jù), 得到可用于后續(xù)分析使用的 clean reads; 通過Flash軟件對 clean reads進行拼接組裝, 組裝時設(shè)置overlap參數(shù)為200 bp, 以保證拼接序列的質(zhì)量, 使用QIIME(version 1.9.1)軟件對組裝結(jié)果進行質(zhì)量篩選, 保證后續(xù)分析中用到的基因序列同時滿足一下多個標準: 1)序列中不包含缺失, 2)單堿基重復(fù)小于8, 3)長度大于200 bp; 為避免嵌合體對后續(xù)分析結(jié)果的影響, 使用UCHIME (version 4.1)軟件去除 clean reads 中的嵌合體數(shù)據(jù), 再通過USEARCH (version 7.1 htp:// drive5.com/uparse)將高質(zhì)量序列以97%的序列相似度為閾值進行OTUs聚類。以上經(jīng)過多次質(zhì)量控制的高質(zhì)量OTU數(shù)據(jù)均已提交至NCBI (MN114085-MN114092)。

    寄生線蟲的多樣性通過Shannon-wiener指數(shù)進行計算, 使用Mothur分析軟件, 在97% OTU水平上計算了Rarefaction分析和alpha多樣性指數(shù)(基于豐度的覆蓋率估算Shannon-Wiener和Simpson指數(shù))。稀釋性曲線由R v.2.8.1中的vegan工具包完成, 以檢測數(shù)據(jù)的可靠性。運用R語言工具對97% 相似水平的OTU樣本進行統(tǒng)計構(gòu)建Venn圖。通過PCA聚類分析兩個湖泊的白額雁種群之間寄生蟲群落結(jié)構(gòu)的差異, 使用R v.2.8.1中的vegan工具包(版本2.0-2)(R Development Core Team, 2006)繪制熱圖。統(tǒng)計多個樣本中所共有以及獨有的OTU數(shù)目, 通過層次聚類中的加權(quán)組平均法UPGMA構(gòu)建進化樹。

    2 結(jié)果與分析

    2.1 高通量測序結(jié)果分析

    對60個糞便樣品進行特異性進行測序, 進行序列的過濾和修剪后, 共獲得3319530條reads, 對得到的序列進行不同分類水平上的鑒定, 其中3299430條reads鑒定為腸道線蟲。采用對序列進行隨機抽樣的方法, 以抽到的序列數(shù)與它們所能代表OTU的數(shù)目構(gòu)建稀釋性曲線, 稀釋性曲線在測序量較少時, 隨著測序量的增加OTU的數(shù)量會顯著增加, 而隨著測序量的逐漸增大新增OTU的數(shù)量逐漸變緩, 并趨于平坦, 說明測序數(shù)據(jù)量漸進合理, 本研究中的樣品測序數(shù)據(jù)量達到20000條時, 稀釋曲線趨于平坦, 說明本研究中的寄生蟲多樣性符合實際情況(圖1)。所有序列經(jīng)過聚類后共得到3249個腸道寄生線蟲OTU, 每個樣本的OTU數(shù)目從1 到1408不等(97 %相似), 其中3028個OTU(93.2 %)屬于兩個湖泊共有, 而其他19個OTU和16個OTU分別隸屬于菜子湖和升金湖(圖2)。

    2.2 腸道寄生線蟲的組成

    從60個糞便樣品中獲得的所有序列隸屬于4個目, 91屬和191種(圖3)。通過對兩個湖泊越冬白額雁腸道線蟲OTU含量構(gòu)建韋恩圖, 在屬的水平, 菜子湖和升金湖分別鑒定到84個屬和51個屬, 占總數(shù)的89.6%和63.7%, 其中51個屬為兩湖共有(56.0%); 而在種的水平, 菜子湖鑒定到171個線蟲物種(64.5%), 升金湖鑒定到93個物種(51.4%), 其中87個種為兩湖共有(48.1%)。鑒定出線蟲的主要目為小桿目(Rhabditida)(62.12 %)、墊刃目(Tylenchida)(32.20%)、疏毛目(Araeolaimida)(3.40%)和單宮目(Monhysterida)(2.27%), 其中小桿目為優(yōu)勢物種(圖3)。

    圖1 菜子湖(CZ)和升金湖(SJ)越冬白額雁腸道寄生線蟲高通量測序的稀釋曲線

    Figure 1 Intestinal parasitic nematode characteristics between Caizi Lake (CZ)and Shengjin Lake(SJ)

    圖2 菜子湖(CZ)和升金湖(SJ)白額雁腸道寄生線蟲重疊和特異的OTUs

    Figure 2 OTU richness for intestinal parasitic nematode from lesser white-fronted geese from CZ and SJ Lakes

    2.3 Alpha多樣性分析

    為進一步了解菜子湖和升金湖越冬白額雁腸道寄生蟲豐度和多樣性變化, 基于97% 相似度歸類的OTU計算60個樣品alpha多樣性指數(shù), Shannon- wiener指數(shù)和Simpson指數(shù)的結(jié)果表明, 菜子湖和升金湖越冬白額雁腸道寄生線蟲的多樣性無顯著差異(> 0.05; 圖5)。

    圖3 60個白額雁糞便樣本中的線蟲分布圖

    Figure 3 Distribution of the 4 main nematode orders found in 60 fecal samples from wintering greater white-fronted geese.

    2.4 Beta多樣性分析

    PCA聚類結(jié)果顯示, 兩個湖泊白額雁腸道寄生線蟲的群落組成相似性較高, 部分樣本在兩個湖泊之間呈現(xiàn)不規(guī)則分布狀態(tài)(圖6)?;谌郝浣M成在OTU水平上進行層次聚類分析顯示, 系統(tǒng)發(fā)育關(guān)系分為5大類群, 形成的拓撲結(jié)構(gòu)與菜子湖和升金湖的地理區(qū)域匹配度不高(圖7)。

    3 討論

    寄生蟲的感染往往導(dǎo)致宿主免疫力下降、體重減輕, 甚至死亡, 感染較嚴重的宿主在面臨疾病, 捕食者等不利因素影響時極易受到威脅[5,24-25]。白額雁作為一種重要的濕地指示物種, 隨著人為干擾和生境的變化, 在遷徙途中, 增加了寄生蟲感染的風(fēng)險, 然而, 相較于其他鳥類, 該物種的腸道寄生線蟲更具有代表性[26]。本研究對菜子湖和升金湖越冬白額雁的腸道寄生線蟲多樣性研究表明, 在本研究兩個湖泊越冬白額雁種群中, 最常見的腸道寄生線蟲為小桿目(64.62%,=60), 其中包括幾種重要的小桿目物種, 例如、、和等。小桿目通常存在于土壤環(huán)境、脊椎動物和人類體內(nèi), 其中濕地是其主要棲息環(huán)境, 生境中豐富的水生動物和水資源為該物種提供了極為豐富的生存資源。此外, 線蟲中較為常見的還有蛔蟲屬(), 主要寄生在家禽與遷徙水鳥體內(nèi), 可通過水生環(huán)境與周圍生境中的物種傳播, 能夠引起等嚴重的人禽共患病[27]。先前有研究在白頭鶴糞便樣品中檢測出5種寄生線蟲, 分別為毛細線蟲(spp.)、類圓線蟲(spp.)、禽蛔蟲(spp.)、毛圓線蟲(spp.)和鉤蟲(Ancylostomatidae), 其中感染率最高的為禽蛔蟲[28], 該屬物種也是線蟲中較為常見的人禽共患寄生蟲病的病原體之一, 在本次研究中也發(fā)現(xiàn)了該屬物種在菜子湖和升金湖的越冬白額雁種群中廣泛分布。經(jīng)過以上研究的初步統(tǒng)計, 已發(fā)現(xiàn)多個嚴重致病病原體的存在, 其中本研究結(jié)果中所發(fā)現(xiàn)的和等小桿目物種均極易感染人類, 可能造成人蛔蟲病的發(fā)生和流行, 人蛔蟲病是一種易被忽視的寄生蟲疾病, 它是由腸道寄生線蟲感染所引起, 全球約有7億人感染, 成人感染蛔蟲病時常出現(xiàn)食欲不振、腹痛, 腹瀉、便秘、蕁麻疹等癥狀, 嚴重時可導(dǎo)致腸道穿孔甚至死亡; 兒童感染時則會出現(xiàn)流涎、磨牙、煩躁不安等, 重者出現(xiàn)營養(yǎng)不良等癥狀[29-31]。近些年隨著人類活動對濕地生態(tài)的干擾, 白額雁與人類直接或間接接觸的幾率大大增高, 兩個湖泊濕地附近人類生活區(qū)域的公共衛(wèi)生風(fēng)險也隨之增加。

    圖4 韋恩圖顯示菜子湖(CZ)和升金湖(SJ)越冬白額雁糞便樣品中線蟲屬水平(A)和種水平(B)的異同

    Figure 4 Nematode genus (A) and species (B) compositions from the samples

    圖5 菜子湖(CZ)和升金湖(SJ)越冬白額雁腸道寄生線蟲的Shannon-Wiener指數(shù)(左)和Simpson指數(shù)(右)

    Figure 5 Shannon-Wiener index (left) and Simpson index (right) for the intestinal parasitic nematodes in the wintering greater white-fronted goose fecal samples from CZ lake and SJ lake

    圖6 兩個湖泊物種基于加權(quán)UniFrac的PCA分析(藍色代表來自升金湖的樣本, 紅色代表來自菜子湖的樣本)

    Figure 6 Principal coordinate analysis plot of weighted UniFrac distances for the species sampled from two lake s(Each point represents the gut parasitic nematode community of an individual goose. Blue=greater white-fronted goose from Shengjin Lake; Red=greater white-fronted goose from Caizi Lake)

    腸道寄生蟲群落具有一定的可塑性, 腸道寄生蟲的組成和結(jié)構(gòu)可能會受到遺傳因素(宿主本身)的影響, 有研究人員發(fā)現(xiàn), 棲息于鄱陽湖的幾種雁體內(nèi)腸道寄生蟲和微生物的差異性顯著, 表明了遺傳因素確實會影響腸道寄生蟲和微生物的組成和結(jié)構(gòu)[26]。本研究通過alpha多樣性分析發(fā)現(xiàn), 升金湖和菜子湖越冬白額雁的腸道寄生線蟲群落差異性不顯著, 說明在兩個湖泊之間的越冬白額雁腸道寄生蟲之間的群落結(jié)構(gòu)相似, 這與腸道菌群結(jié)構(gòu)有所不同[9], 說明宿主并非白額雁腸道寄生線蟲群落結(jié)構(gòu)形成的決定性因素, 可能是由于菜子湖和升金湖為越冬白額雁提供了相似的食物來源。腸道寄生蟲的分化受到了這些因素的約束, 因此在兩個研究區(qū)域間并未產(chǎn)生明顯的差異。

    過高的種群密度被認為能夠增加寄生蟲感染的風(fēng)險, 并且在寄生蟲傳播中起到關(guān)鍵性作用, 有研究報道稱, 日本越冬白頭鶴種群球蟲感染率高達90% 以上, 而中國長江中下游越冬白頭鶴的球蟲感染率明顯較低, 表明了種群密度在寄生蟲傳播中起到關(guān)鍵性作用[26]。近些年來, 隨著經(jīng)濟發(fā)展, 大量原始濕地被不斷開墾利用, 人為擾動不斷增強, 候鳥越冬湖泊濕地生態(tài)過程受到嚴重干擾, 迫使部分越冬白額雁種群集中在更小的棲息地中, 產(chǎn)生了較高的種群密度, 增加了個體間的接觸, 使得白額雁種群面臨寄生蟲疾病暴發(fā)的風(fēng)險增加。除此之外, 棲息地環(huán)境和食性也是影響動物腸道寄生蟲多樣性的主要因素, 能夠影響腸道寄生蟲群落的組成結(jié)構(gòu), 有研究表明, 在同一棲息地中, 宿主間的腸道寄生蟲的組成、種類和豐度具有高度相似性[8,26,31]。本研究發(fā)現(xiàn)的所有OTUs中, 有93.2%的OTUs是兩個湖的白額雁共有的, 菜子湖和升金湖分別僅發(fā)現(xiàn)了19個(0.58%)和16個(0.49%)不同的OTUs, 說明兩個湖泊的白額雁腸道寄生線蟲組成及豐度存在差異, 但在OTU水平的韋恩圖(圖 4)顯示, 兩個湖泊白額雁的腸道線蟲相似性更高。另外有研究表明, 越冬白額雁腸道中的線蟲群落結(jié)構(gòu)可能與從覓食生境(例如稻田和草地)獲得的食物有關(guān)[26]。隨著人類社會經(jīng)濟的發(fā)展, 鳥類覓食生境正面臨著不斷衰退、減少和碎片化的威脅, 菜子湖和升金湖作為淺水內(nèi)陸湖泊, 環(huán)境復(fù)雜, 附近不僅有大量的水稻田, 還有植物資源豐富的草灘, 以及冬季枯水期暴露出來的大量的泥灘, 本研究發(fā)現(xiàn)兩個湖泊中白額雁的越冬食物來源差異不大, 且食性相似, 均以菜子湖和升金湖的稻田和草地為覓食地, 喜食苔草、禾本科植物和雙子葉植物等植物。此外, 生態(tài)環(huán)境因素如宿主對棲息地的利用或者其他資源的可獲得性, 也能夠影響寄生蟲集群的組成[32]。有研究發(fā)現(xiàn), 寄生蟲群落組成的相似性往往隨著地理距離或環(huán)境因素呈負相關(guān)[33]。本研究中兩個地區(qū)寄生蟲的感染種類和感染率雖然存在變化, 但其豐富度指數(shù)、多樣性指數(shù)并不具有顯著的統(tǒng)計學(xué)差異, 相似的多樣性系數(shù)表明, 大多數(shù)物種均存在于兩個湖泊, 且一些種類處于較高的優(yōu)勢地位。以上結(jié)果表明, 棲息地的環(huán)境因素對白額雁腸道線蟲群落組成影響不明顯, 并不是決定線蟲系統(tǒng)發(fā)育和群落結(jié)構(gòu)的重要因素, 可能是由于菜子湖與升金湖提供的生態(tài)棲息環(huán)境相似。

    注: 圖中CZ與SJ后標數(shù)字均為菜子湖與升金湖的糞便樣品名稱。

    Figure 7 Weighted UniFrac tree analysis of the intestinal parasitic nematodes isolated from the wintering greater white-fronted goose fecal samples

    生境干擾可能為部分病原體創(chuàng)造適宜的生長環(huán)境和傳播條件, 從而導(dǎo)致疾病爆發(fā), 因此, 監(jiān)測寄生蟲多樣性對于完善濕地管理和公共衛(wèi)生具有非常重要的研究意義[5]。鑒于兩個湖泊中寄生線蟲群落結(jié)構(gòu)的相似性以及與白額雁遷徙的限制, 結(jié)果表明, 由于相似的濕地環(huán)境與食物來源, 寄生蟲可能無法被隔離在某一區(qū)域內(nèi)。另外本研究發(fā)現(xiàn)了許多潛在的致病性腸道寄生線蟲, 例如,和, 其中包括引起哺乳動物和鳥類疾病的幾種寄生蟲。本研究結(jié)果表明, 在菜子湖和升金湖的越冬白額雁種群中, 腸道寄生線蟲感染對人和動物更為敏感, 兩個湖泊的寄生蟲流行極有可能在越冬期發(fā)生, 從而引起人類疾病和潛在的公共健康風(fēng)險。

    4 結(jié)論

    本研究通過對越冬白額雁腸道寄生線蟲的特異性區(qū)域高通量測序, 分析其感染種類、多樣性及其群落組成結(jié)構(gòu)。將3299430條reads歸屬為腸道寄生線蟲門, 共定義3249個腸道寄生線蟲的OTUs, 鑒定為4個目, 91屬和191種, 分別為疏毛目、單宮目、小桿目和墊刃目, 最常見的是小桿目。結(jié)果顯示, 菜子湖越冬白額雁的腸道寄生線蟲alpha多樣性高于升金湖越冬白額雁, 且兩個湖泊中白額雁腸道寄生線蟲群落的組成結(jié)構(gòu)相似。升金湖和菜子湖越冬白額雁的腸道寄生線蟲群落相似性較高, 可能是由于食物來源相似、食性相似、棲息地環(huán)境類似等因素導(dǎo)致兩個湖泊的越冬白額雁群落結(jié)構(gòu)相似。

    [1] RAMEY A, PEARCE J, FLINT P, et alIntercontinental reassortment and genomic variation of low pathogenic avian influenza viruses isolated from northern pintails () in Alaska: examining the evidence through space and time[J]. Virology, 2010, 401: 179–189.

    [2] REED C, BRUDEN D, BYRD K K, et alCharacterizing wild bird contact and seropositivity to highly pathogenic avian influenza A (H5N1) virus in Alaskan residents[J]. Influenza and Other Respiratory Viruses, 2014, 8(5): 516–523.

    [3] CASALINSA L, MOLINAB M I, GULLOC M W, et al(Nematoda: Diplotriae-nidae) causing severe disease in a new host from Argentine Patagonia:(Aves: Picidae)[J]. Parasites and Wildlife, 2019, 8: 106–110.

    [4] CHAGAS C R, VALKIūNAS G, MONTEIRO E F, et alDiversity and distribution of avian malaria and related haemosporidian parasites in captive birdsfrom a Brazilian megalopolis[J]. Malaria journal, 2017, 1716(1): 83.

    [5] HUANG Wei, ZHOU Lizhi, ZHAO Niannian. Temporal- spatial patterns of intestinal parasites of the Hooded crane () wintering in lakes of the middle and lower Yangtze River floodplain[J]. Avian Research, 2014, 5(1): 40–48.

    [6] LIU Gang, LI Qingyue, WANG Chong, et al. The complete mitochondrial genome offrom the wintering greater white-fronted goose in Shengjin Lake, China, and phylogenetic relationships amongspecies [J]. Parasitology Research, 2019, 118: 1299–1306.

    [7] HONMA H, SUYAMA Y, WATANABE Y, et alAccurate analysis of prevalence of coccidiosis in individually identified wild cranes in inhabiting and migrating popula-tions in Japan[J]. Environmental microbiology, 2011, 13: 2876–2887.

    [8] ZHAO Guanghong, DONG Yuanqiu, ZHOU Lizhi, et alThe gut microbiome of hooded cranes () wintering at Shengjin Lake, China [J]. Microbiology Open. 2017, 6: e447.

    [9] XIANG Xingjia, ZHANG Fengling, FU Rong, et alSignificant differences in bacterial and potentially pathogenic communities between sympatric hooded crane and greater white-fronted goose[J]. Fronties in Micro-biology, 2019, 5(10): 163.

    [10] 王榮軍, 朱金鳳, 張龍現(xiàn), 等. 鄭州動物園鶴球蟲感染調(diào)查[J]. 河南農(nóng)業(yè)大學(xué)學(xué)報, 2007, 41(1): 61–63.

    [11] 趙世華, 林福德, 常建華, 等. 黑鸛的呵欠鳥凱撒瑪吸蟲病[J]. 野生動物, 1990, 11(2): 25.

    [12] 李清武, 姜建民, 陳明非, 等.丹頂鶴卷棘口吸蟲病[J]. 中國獸醫(yī)雜志, 2004, 40(6): 53.

    [13] 車桂翠, 張恩珠, 責(zé)亞華. 動物園動物體內(nèi)寄生蟲的調(diào)查[J]. 中國獸醫(yī)雜志, 2000, 26(4): 26.

    [14] 孫銘飛. 河南省禽類隱孢子蟲種類鑒定及分子種系發(fā)育關(guān)系研究[D]. 鄭州: 河南農(nóng)業(yè)大學(xué), 2005.

    [15] FANKE J, WIBBELT G, KRONE O. Mortality factors and diseases in free-ranging Eurasian cranes () in Germany[J]. Journal of Wildlife Diseases, 2011, 47: 627–637.

    [16] ATKINSON C T, DUSEK R J, WOODS K L, et alPathogenicity of avian malaria inexperimentally-infected[J]. Journal of Wildlife Diseases, 2000, 36(2): 197–204.

    [17] BARTER M, CAO L, CHEN L,et al. Results of a survey for waterbirds in the lower Yangtze floodplain, China, in January-February 2004[J]. Forktail, 2005, 21: 1–7.

    [18] CAO Lei, BARTER M, ZHAO Meijuan, et al. A systematic scheme for monitoring waterbird populations at Shengjin Lake, China: methodology and preliminary results[J]. Chinese Birds, 2011, 2: 1–17.

    [19] ZHANG Dongmei, ZHOU Lizhi, SONG Yunwei. E?ect of water level ?uctuations on temporal-spatial patterns of foraging activities by the wintering Hooded Crane ()[J]. Avian Research, 2015, 6:16.

    [20] LIU Jinzhen, XIAO Haixia, LEI Fumin, et alHighly pathogenic H5N1 influenza virus infection in migratory birds[J]. Science, 2005, 309: 1206.

    [21] OLSEN B, MUNSTER V J, WALLENSTEN A, et alGlobal patterns of influenza a virus in wild birds[J].Science, 2006, 312(5772): 384–388.

    [22] GAUTHIER-CLERC M, LEBARBENCHON C, THOMAS F. Recent expansion of highly pathogenic avian influenza H5N1: a critical review [J].Ibis, 2007, 149 (2): 202–214.

    [23] DAGLEISH M P, RYAN P G, GIRLING S, et al.Clinical pathology of the critically endangered gough bunting ()[J]. Journal of Comparative Patho-logy, 2017, 156(2-3): 264–274.

    [24] MOWLAVI G R, MASSOUD J, MOBEDI I, et alTetrameres (Tetrameres) grusi (Shumakovich, 1946) (Nematoda: Tetrameridae) in Eurasian cranes () in Central Iran[J]. Journal of Wildlife Diseases, 2006, 42: 397–401.

    [25] NOVILLA M N, CARPENTER J W. Pathology and pathogenesis of disseminated visceral coccidiosis in cranes[J]. Avian Pathology, 2004, 33: 275–280.

    [26] YANG Yuzhan, DENG Ye, CAO Lei. Characterising the interspecific variations and convergence of gut microbiota in Anseriformes herbivores at wintering areas[J]. Scientific Reports, 2016, 6: 32655.

    [27] PIETERSE A, MALAN A P, ROSS J L. Nematodes that associate with terrestrial molluscs as definitive hosts, including(Rhabditida: Rhabditidae) and its development as a biological molluscicide[J]. Journal of Helminthology, 2017, 91(5): 517–527.

    [28] 李博. 越冬白頭鶴()腸道寄生線蟲分子檢測方法的建立[D]. 合肥. 安徽大學(xué). 2015.

    [29] PULLAN R L, BROOKER S. Ascaris lumbricoides and ascariasis: estimating numbers infected and burden of disease//Holland C. Ascaris the neglected parasite[J]. London: Academic Press, 2013.

    [30] 許隆祺, 陳穎丹, 孫鳳華, 等. 全國人體重要寄生蟲病現(xiàn)狀調(diào)查報告[J]. 中國寄生蟲學(xué)與寄生蟲病雜志, 2005, 23(5): 14–22.

    [31] 尚煜, 湯林華. 兒童土源性線蟲感染疾病負擔(dān)的評價[J]. 中國婦幼保健, 2010, 25(34): 5149–5151.

    [32] ECKBURG P, BIK E M, BERNSTEIN C N, et alDiversity of the human intestinal microbial flora[J]. Science, 2005, 308: 1635–1638.

    [33] ESCH G, SHOSTAK A W, MARCOGLIESE D J, et alPatterns and processes in helminth parasite communities: an overview[J]. Springer Netherlands, 1990, 1–19.

    [34] POULIN R, BLANAR C A, THIELTGES D, et alThe biogeography of parasitism in sticklebacks: distance, habitat differences and the similarity in parasite occurrence and abundance[J]. Ecography, 2011, 34: 540–551.

    Intestinal parasitic nematode diversity in communities of Wintering Greater White-Fronted Geese () between Caizi Lake and Shengjin Lake in China

    LI Qingyue, Gong Zhizhong, Feng Jiahui, LIU Gang*

    School of Life Sciences, Anhui Medical University, Hefei 230032, China

    Intestinal parasites are commonly reported in birds, which can spread to other animals including humans, causing serious zoonotic diseases. In this study, 60 fecal samples were collected fromoverwintering at Caizi Lake (CZ) and Shengjin Lake (SJ) in China. High-throughput sequencing was performed from intestinal parasitic nematodes in these samples. We analyzed the species, diversity and community composition of intestinal parasitic nematode infections of white-fronted geese in two lakes. In total, 3299430 valid reads were classified as belonging to intestinal parasitic nematode phyla, with 3249 intestinal parasitic nematode operational taxonomic units. Nearly all reads from the 60 fecal samples were assignable to 4 orders, 91 genera, and 191 species. The four orders were Rhabditida (62.12%), Tylenchida (32.20%), Araeolaimida (3.40%) and Monhysterida (2.27%), with the predominant order being Tylenchida. The Shannon-Wiener and Simpson indices showed that intestinal parasitic Nematode alpha diversity in the wintering greater white-fronted geese from CZ was not significantly higher than that of the geese from SJ (>0.05). Principle component analysis (PCA) showed that the intestinal parasitic nematode community compositions differed nonsignificantly between the guts of the greater white-fronted geese from both lakes. The phylogenetic trees showed that the samples were poorly matched with their lakes, and some samples were irregularly distributed. The intestinal parasitic nematode community composition and diversity among thediffered nonsignificantly between CZ and SJ, suggesting that the hosts were not the crucial factor in shaping the intestinal parasitic nematode structure. This is because CZ and SJ provide the geese with abundant and diverse food sources, similar diets and habitat environment.

    intestinal parasitic nematodes;; diversity

    10.14108/j.cnki.1008-8873.2022.06.024

    Q341

    A

    1008-8873(2022)06-202-09

    2020-10-26;

    2020-11-20

    國家自然科學(xué)基金項目(31702030); 安徽省高校自然科學(xué)基金重點項目(KJ2015A002, KJ2021A0246)

    李清月(1994—), 女, 江蘇連云港人, 碩士, 主要從事進化與分子遺傳學(xué)研究, E-mail: 1051897431@qq.com

    通信作者:劉剛, 男, 博士, 副教授, 主要從事進化與分子遺傳學(xué)研究, E-mail: liugang8966@163.com

    李清月, 龔治忠, 馮佳慧, 等. 升金湖和菜子湖越冬白額雁腸道寄生線蟲多樣性研究[J]. 生態(tài)科學(xué), 2022, 41(6): 202–210.

    LI Qingyue, Gong Zhizhong, Feng Jiahui, et al. Intestinal parasitic nematode diversity in communities of Wintering Greater White-Fronted Geese () between Caizi Lake and Shengjin Lake in China[J]. Ecological Science, 2022, 41(6): 202–210.

    猜你喜歡
    菜子金湖線蟲
    衣食父母(外一首)
    夏季蔬菜換茬期線蟲防治要注意
    朱光武和他的金湖宏源國際大酒店
    華人時刊(2021年21期)2021-03-09 05:31:28
    菜子湖紀事(外一首)
    遼河(2020年11期)2020-11-28 07:19:34
    雨日寄友
    地黃花對秀麗線蟲壽命的影響
    中成藥(2018年2期)2018-05-09 07:20:04
    朝鮮孢囊線蟲——浙江省孢囊線蟲新記錄種
    線蟲共生菌Xenorhabdus budapestensis SN19次生代謝產(chǎn)物的分離純化與結(jié)構(gòu)鑒定
    荷都金湖明天更美好
    安徽園林(2015年3期)2015-10-20 02:38:07
    亚洲男人的天堂狠狠| 亚洲精品色激情综合| 亚洲最大成人中文| 国产av一区在线观看免费| www.www免费av| 国产av不卡久久| 亚洲国产高清在线一区二区三| 色尼玛亚洲综合影院| 一区二区三区四区激情视频 | 丰满人妻一区二区三区视频av| 国产白丝娇喘喷水9色精品| 十八禁国产超污无遮挡网站| aaaaa片日本免费| 免费大片18禁| 亚洲一级一片aⅴ在线观看| 少妇高潮的动态图| 国产真实乱freesex| 国产男人的电影天堂91| 久久国产精品人妻蜜桃| 国产淫片久久久久久久久| 一夜夜www| 一区福利在线观看| 成人综合一区亚洲| 日韩高清综合在线| 深夜精品福利| 日本五十路高清| 日本成人三级电影网站| 国产探花在线观看一区二区| 久9热在线精品视频| 亚洲美女视频黄频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲男人的天堂狠狠| 一个人观看的视频www高清免费观看| 午夜免费成人在线视频| 午夜免费激情av| a级毛片a级免费在线| 日日啪夜夜撸| 天堂影院成人在线观看| 国产淫片久久久久久久久| 午夜免费激情av| 色播亚洲综合网| 免费看美女性在线毛片视频| 久久久久久国产a免费观看| 亚洲国产欧洲综合997久久,| 白带黄色成豆腐渣| 欧美高清成人免费视频www| 日韩大尺度精品在线看网址| 十八禁国产超污无遮挡网站| 亚洲不卡免费看| 美女高潮的动态| 麻豆av噜噜一区二区三区| 在线播放国产精品三级| 蜜桃久久精品国产亚洲av| 国产精品自产拍在线观看55亚洲| 久久久久精品国产欧美久久久| 午夜福利欧美成人| 热99在线观看视频| 亚洲综合色惰| 精品国产三级普通话版| 亚洲成人免费电影在线观看| 国产白丝娇喘喷水9色精品| 久99久视频精品免费| 99久久精品热视频| 日韩欧美国产一区二区入口| 成人国产麻豆网| 三级国产精品欧美在线观看| 国产伦在线观看视频一区| 国产午夜精品论理片| 亚洲va日本ⅴa欧美va伊人久久| 一进一出抽搐动态| 久9热在线精品视频| 简卡轻食公司| 非洲黑人性xxxx精品又粗又长| 欧美3d第一页| 欧美成人一区二区免费高清观看| 性色avwww在线观看| 美女黄网站色视频| 草草在线视频免费看| 一区二区三区激情视频| 国产精品一区二区三区四区免费观看 | 色综合站精品国产| 国产精品日韩av在线免费观看| 91麻豆精品激情在线观看国产| 久久99热这里只有精品18| 久久婷婷人人爽人人干人人爱| 国产不卡一卡二| 日韩中字成人| 人人妻人人澡欧美一区二区| 免费在线观看影片大全网站| 国产精品亚洲一级av第二区| 国产毛片a区久久久久| 亚洲专区国产一区二区| 伊人久久精品亚洲午夜| 91狼人影院| 日韩欧美国产在线观看| 无人区码免费观看不卡| 无人区码免费观看不卡| 真人一进一出gif抽搐免费| 欧美xxxx性猛交bbbb| 一进一出抽搐动态| 国内精品久久久久久久电影| 日韩亚洲欧美综合| 国产不卡一卡二| 亚洲精品日韩av片在线观看| 1024手机看黄色片| 人妻少妇偷人精品九色| 欧美色欧美亚洲另类二区| 婷婷丁香在线五月| 很黄的视频免费| 男插女下体视频免费在线播放| 国产精品爽爽va在线观看网站| 亚洲精品粉嫩美女一区| 极品教师在线免费播放| 免费高清视频大片| 国产一区二区激情短视频| 亚洲专区中文字幕在线| 免费无遮挡裸体视频| 亚洲专区国产一区二区| 国产亚洲精品综合一区在线观看| 国产探花在线观看一区二区| 内地一区二区视频在线| 婷婷丁香在线五月| 别揉我奶头 嗯啊视频| 亚洲人成网站在线播| 精品日产1卡2卡| 在线免费观看的www视频| www日本黄色视频网| 国产一区二区三区av在线 | 在线免费观看的www视频| 亚洲一区二区三区色噜噜| 男女啪啪激烈高潮av片| 久久香蕉精品热| 国产精品亚洲一级av第二区| 欧美激情久久久久久爽电影| 精品日产1卡2卡| 免费在线观看成人毛片| 欧美xxxx性猛交bbbb| 一进一出抽搐动态| 国产精品伦人一区二区| 亚洲人与动物交配视频| 99热6这里只有精品| bbb黄色大片| 欧美色欧美亚洲另类二区| 日本黄色片子视频| 精品久久久久久久末码| 亚州av有码| 村上凉子中文字幕在线| 久久精品影院6| 久久婷婷人人爽人人干人人爱| 国产精品亚洲美女久久久| 久久久国产成人免费| 精品久久久久久,| 变态另类成人亚洲欧美熟女| 黄色欧美视频在线观看| 成人欧美大片| 免费观看精品视频网站| 久久久久久久久中文| 精品一区二区三区av网在线观看| 日日摸夜夜添夜夜添小说| 麻豆精品久久久久久蜜桃| 老司机午夜福利在线观看视频| 亚洲国产精品久久男人天堂| 联通29元200g的流量卡| 天堂av国产一区二区熟女人妻| 免费av观看视频| 18+在线观看网站| 欧美激情久久久久久爽电影| 不卡一级毛片| 精品福利观看| 日韩中文字幕欧美一区二区| 精品人妻1区二区| 搡女人真爽免费视频火全软件 | 最近最新免费中文字幕在线| 亚洲va在线va天堂va国产| 午夜精品在线福利| 国产视频内射| 极品教师在线视频| 淫妇啪啪啪对白视频| 亚洲精品456在线播放app | 日本-黄色视频高清免费观看| 赤兔流量卡办理| 特级一级黄色大片| 永久网站在线| 99久国产av精品| 亚洲人成网站高清观看| 国产伦一二天堂av在线观看| h日本视频在线播放| 看免费成人av毛片| 两个人的视频大全免费| 男女下面进入的视频免费午夜| 精品久久久久久久久亚洲 | 蜜桃亚洲精品一区二区三区| 精品午夜福利视频在线观看一区| 天堂影院成人在线观看| 亚洲性久久影院| 看十八女毛片水多多多| 欧美日韩综合久久久久久 | 网址你懂的国产日韩在线| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品在线观看二区| 色在线成人网| 少妇被粗大猛烈的视频| 亚洲国产欧洲综合997久久,| 国产精品爽爽va在线观看网站| 久久久久久九九精品二区国产| 欧美+亚洲+日韩+国产| 最新中文字幕久久久久| 天堂√8在线中文| 亚洲无线在线观看| 国产精品爽爽va在线观看网站| 一区二区三区四区激情视频 | 在线a可以看的网站| 联通29元200g的流量卡| 观看免费一级毛片| 国产精品不卡视频一区二区| 小蜜桃在线观看免费完整版高清| 91在线观看av| 免费av不卡在线播放| 亚洲经典国产精华液单| 一区二区三区高清视频在线| 深夜a级毛片| 91麻豆精品激情在线观看国产| 国产精品,欧美在线| 嫩草影院精品99| 日本 欧美在线| 制服丝袜大香蕉在线| 日韩欧美免费精品| 一边摸一边抽搐一进一小说| 无人区码免费观看不卡| 蜜桃亚洲精品一区二区三区| 免费大片18禁| 女人十人毛片免费观看3o分钟| 97碰自拍视频| 久久精品国产清高在天天线| 淫秽高清视频在线观看| 日本免费a在线| 一区二区三区免费毛片| 熟女人妻精品中文字幕| 日本在线视频免费播放| 91精品国产九色| 小蜜桃在线观看免费完整版高清| 午夜视频国产福利| 尾随美女入室| 成人特级黄色片久久久久久久| 久久欧美精品欧美久久欧美| av在线观看视频网站免费| 久久久久久伊人网av| 国产高清视频在线观看网站| 亚洲狠狠婷婷综合久久图片| 在线免费观看的www视频| a级毛片a级免费在线| 午夜福利欧美成人| 国内揄拍国产精品人妻在线| АⅤ资源中文在线天堂| www日本黄色视频网| 亚洲国产精品sss在线观看| 亚洲国产高清在线一区二区三| 欧美区成人在线视频| 国产精品久久久久久亚洲av鲁大| 亚洲综合色惰| 五月伊人婷婷丁香| 九色成人免费人妻av| 免费高清视频大片| 免费大片18禁| 亚洲av免费在线观看| 欧美在线一区亚洲| 狂野欧美白嫩少妇大欣赏| 日韩欧美在线二视频| 中文在线观看免费www的网站| 可以在线观看毛片的网站| 国产精品国产高清国产av| 成人三级黄色视频| 免费不卡的大黄色大毛片视频在线观看 | 日韩欧美在线乱码| 搞女人的毛片| 成人av一区二区三区在线看| 美女大奶头视频| 日日撸夜夜添| 我要看日韩黄色一级片| 女人十人毛片免费观看3o分钟| 亚洲精华国产精华精| 国产精品美女特级片免费视频播放器| 成人亚洲精品av一区二区| 国产探花在线观看一区二区| 级片在线观看| 成人综合一区亚洲| 欧美激情久久久久久爽电影| 永久网站在线| 精品人妻视频免费看| 国产国拍精品亚洲av在线观看| 天美传媒精品一区二区| 禁无遮挡网站| 免费看日本二区| 给我免费播放毛片高清在线观看| 国产欧美日韩精品亚洲av| 国产一区二区亚洲精品在线观看| 91久久精品国产一区二区成人| 能在线免费观看的黄片| 欧美日韩亚洲国产一区二区在线观看| av在线老鸭窝| 中文字幕av在线有码专区| 欧美最新免费一区二区三区| 波多野结衣高清无吗| 我的女老师完整版在线观看| 午夜亚洲福利在线播放| 美女cb高潮喷水在线观看| 中出人妻视频一区二区| 欧美潮喷喷水| 看十八女毛片水多多多| 亚洲专区中文字幕在线| 日本 av在线| 亚洲最大成人av| 亚洲av免费高清在线观看| 中文字幕免费在线视频6| 直男gayav资源| 伦精品一区二区三区| 国产一级毛片七仙女欲春2| 亚洲精品日韩av片在线观看| 国产男人的电影天堂91| 久久久精品欧美日韩精品| 丝袜美腿在线中文| 美女 人体艺术 gogo| 国产 一区 欧美 日韩| 日本熟妇午夜| 91久久精品电影网| 嫩草影视91久久| 久久精品久久久久久噜噜老黄 | 日韩精品有码人妻一区| 露出奶头的视频| 亚洲四区av| 女的被弄到高潮叫床怎么办 | 男女那种视频在线观看| 欧美+日韩+精品| 大型黄色视频在线免费观看| 午夜a级毛片| 制服丝袜大香蕉在线| 免费在线观看成人毛片| 大又大粗又爽又黄少妇毛片口| 小说图片视频综合网站| 色噜噜av男人的天堂激情| 免费观看人在逋| 春色校园在线视频观看| 国产视频内射| 又黄又爽又刺激的免费视频.| 免费无遮挡裸体视频| 久久欧美精品欧美久久欧美| 久久婷婷人人爽人人干人人爱| 国产精品久久久久久精品电影| 少妇裸体淫交视频免费看高清| 少妇人妻精品综合一区二区 | 午夜激情欧美在线| 内地一区二区视频在线| 看黄色毛片网站| 成年女人看的毛片在线观看| 色噜噜av男人的天堂激情| 超碰av人人做人人爽久久| aaaaa片日本免费| 国产色爽女视频免费观看| 黄色一级大片看看| 亚洲真实伦在线观看| 亚洲av成人av| www.www免费av| 啦啦啦啦在线视频资源| a级一级毛片免费在线观看| 国产一区二区三区视频了| 午夜福利欧美成人| 噜噜噜噜噜久久久久久91| 久久久久国内视频| 极品教师在线视频| 国产av麻豆久久久久久久| 熟妇人妻久久中文字幕3abv| 欧美绝顶高潮抽搐喷水| 久久精品人妻少妇| 亚洲18禁久久av| 性插视频无遮挡在线免费观看| 久9热在线精品视频| 最近中文字幕高清免费大全6 | 亚洲人成网站高清观看| 欧美又色又爽又黄视频| 在线观看av片永久免费下载| 极品教师在线免费播放| 日韩 亚洲 欧美在线| 九色国产91popny在线| 黄色欧美视频在线观看| 午夜免费男女啪啪视频观看 | .国产精品久久| 成人综合一区亚洲| 免费看日本二区| 直男gayav资源| 国产又黄又爽又无遮挡在线| 悠悠久久av| 国产精品三级大全| 最近在线观看免费完整版| 国产激情偷乱视频一区二区| 永久网站在线| www.色视频.com| 在线免费十八禁| 日韩大尺度精品在线看网址| 中文字幕久久专区| 日本成人三级电影网站| 五月伊人婷婷丁香| 日本一二三区视频观看| 久久国产精品人妻蜜桃| 色综合色国产| 亚洲欧美日韩无卡精品| 一区二区三区四区激情视频 | 欧美黑人巨大hd| 女人被狂操c到高潮| 在线观看66精品国产| 色吧在线观看| 国产成人av教育| 看黄色毛片网站| 日韩精品青青久久久久久| 最后的刺客免费高清国语| 免费人成在线观看视频色| 日韩中字成人| 国产精品三级大全| 日本欧美国产在线视频| 国产亚洲精品综合一区在线观看| 欧美一区二区精品小视频在线| 国内揄拍国产精品人妻在线| 毛片女人毛片| 久久久久久久亚洲中文字幕| 久久香蕉精品热| 国产精品久久久久久精品电影| 亚洲欧美日韩无卡精品| 一级a爱片免费观看的视频| 美女黄网站色视频| 亚洲成人精品中文字幕电影| 97超视频在线观看视频| av福利片在线观看| 国产 一区精品| 少妇的逼好多水| 一级黄片播放器| 国产亚洲欧美98| 悠悠久久av| 色综合亚洲欧美另类图片| 亚洲黑人精品在线| 成熟少妇高潮喷水视频| 精品一区二区三区av网在线观看| 人妻制服诱惑在线中文字幕| 日本五十路高清| 亚洲精品一区av在线观看| 亚洲av第一区精品v没综合| 免费看日本二区| 久久精品国产亚洲av涩爱 | 亚洲美女视频黄频| 久久久久久久久久黄片| 色综合亚洲欧美另类图片| 给我免费播放毛片高清在线观看| 久久草成人影院| 国产精品不卡视频一区二区| 美女黄网站色视频| 高清毛片免费观看视频网站| 日本熟妇午夜| 亚洲av中文av极速乱 | 国产黄a三级三级三级人| 俺也久久电影网| 老女人水多毛片| 啦啦啦观看免费观看视频高清| 国产 一区精品| 男女做爰动态图高潮gif福利片| 日韩av在线大香蕉| 嫩草影院入口| 国产精品亚洲美女久久久| 国产高清激情床上av| 欧美成人一区二区免费高清观看| 国产探花在线观看一区二区| 神马国产精品三级电影在线观看| 人妻夜夜爽99麻豆av| 最好的美女福利视频网| 欧美成人性av电影在线观看| 精品久久久久久久久久久久久| 人人妻人人澡欧美一区二区| 国产高清有码在线观看视频| 男人狂女人下面高潮的视频| 在线观看免费视频日本深夜| 我要看日韩黄色一级片| 午夜福利成人在线免费观看| 看十八女毛片水多多多| 午夜福利18| 国产高清激情床上av| а√天堂www在线а√下载| 亚洲精品久久国产高清桃花| 一卡2卡三卡四卡精品乱码亚洲| 成人一区二区视频在线观看| 伊人久久精品亚洲午夜| 精品一区二区免费观看| 欧美成人性av电影在线观看| 欧美中文日本在线观看视频| 免费人成在线观看视频色| bbb黄色大片| 久久精品国产鲁丝片午夜精品 | 欧美激情国产日韩精品一区| 久久婷婷人人爽人人干人人爱| 极品教师在线免费播放| 校园春色视频在线观看| 国产精品一区二区性色av| 日韩av在线大香蕉| 国产aⅴ精品一区二区三区波| 老熟妇仑乱视频hdxx| 国产 一区 欧美 日韩| 精品国产三级普通话版| 日本黄色片子视频| 麻豆久久精品国产亚洲av| 一级毛片久久久久久久久女| 日韩强制内射视频| 最后的刺客免费高清国语| 国产精品人妻久久久久久| 国产精品久久久久久亚洲av鲁大| 露出奶头的视频| 日本三级黄在线观看| 国产真实伦视频高清在线观看 | 黄色视频,在线免费观看| 很黄的视频免费| 国内久久婷婷六月综合欲色啪| 成年版毛片免费区| 97碰自拍视频| 亚洲18禁久久av| 亚洲av不卡在线观看| 成人毛片a级毛片在线播放| 制服丝袜大香蕉在线| 午夜激情福利司机影院| av在线亚洲专区| www日本黄色视频网| 国产精品精品国产色婷婷| 亚洲avbb在线观看| 国产成人av教育| 国产麻豆成人av免费视频| 女人被狂操c到高潮| 免费av不卡在线播放| 色哟哟·www| 干丝袜人妻中文字幕| 国产精品福利在线免费观看| 男女做爰动态图高潮gif福利片| 欧美另类亚洲清纯唯美| 成人欧美大片| 国产主播在线观看一区二区| 亚洲国产精品成人综合色| 天堂影院成人在线观看| 国产男靠女视频免费网站| 国产三级中文精品| 成年女人永久免费观看视频| 久久久久久久午夜电影| 人妻丰满熟妇av一区二区三区| 全区人妻精品视频| 国产精品爽爽va在线观看网站| 国产日本99.免费观看| 波野结衣二区三区在线| 亚洲国产欧美人成| 亚洲欧美精品综合久久99| 国产免费一级a男人的天堂| 最近中文字幕高清免费大全6 | 国产精品久久久久久精品电影| 99国产精品一区二区蜜桃av| 乱码一卡2卡4卡精品| 国产三级在线视频| 亚洲四区av| 成年女人看的毛片在线观看| 国产精品女同一区二区软件 | 又爽又黄a免费视频| 级片在线观看| 国产精品99久久久久久久久| 婷婷六月久久综合丁香| 俄罗斯特黄特色一大片| 九九爱精品视频在线观看| 婷婷精品国产亚洲av在线| 91久久精品国产一区二区成人| 久久久色成人| 天堂av国产一区二区熟女人妻| 最好的美女福利视频网| 亚洲欧美日韩高清专用| 国产亚洲精品av在线| 97超视频在线观看视频| 一夜夜www| 99久久精品一区二区三区| 日韩欧美国产一区二区入口| 听说在线观看完整版免费高清| 国产又黄又爽又无遮挡在线| 丰满人妻一区二区三区视频av| 成人鲁丝片一二三区免费| 欧美激情久久久久久爽电影| 久久婷婷人人爽人人干人人爱| 热99re8久久精品国产| 日韩精品中文字幕看吧| 欧美性猛交╳xxx乱大交人| 日本一本二区三区精品| 有码 亚洲区| 国产精品野战在线观看| 五月玫瑰六月丁香| 欧美潮喷喷水| 国产高清视频在线观看网站| 国产精品野战在线观看| 色在线成人网| 内地一区二区视频在线| 久久久久久伊人网av| 校园春色视频在线观看| 一进一出好大好爽视频| 女人被狂操c到高潮| 亚洲国产精品sss在线观看| 亚洲久久久久久中文字幕| 色在线成人网| 亚洲美女搞黄在线观看 | 国产精品嫩草影院av在线观看 | 久久久久国产精品人妻aⅴ院| 国产精品国产三级国产av玫瑰| 两人在一起打扑克的视频| 天堂av国产一区二区熟女人妻| 91麻豆精品激情在线观看国产| 亚洲av不卡在线观看| 日本五十路高清| 亚洲国产色片| 午夜福利在线观看免费完整高清在 | 欧美激情国产日韩精品一区| 亚洲18禁久久av|