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

    基于mtDNA Cyt b 序列分析齊口裂腹魚群體遺傳多樣性

    2017-05-16 09:07:42張爭世胡冰潔葉祥益劉桂嘉鄭曙明
    水生生物學報 2017年3期
    關(guān)鍵詞:腹魚巫溪組群

    張爭世 胡冰潔 葉祥益 劉桂嘉 鄭曙明 葉 華

    (西南大學魚類繁育與健康養(yǎng)殖研究中心, 重慶 402460)

    基于mtDNA Cyt b 序列分析齊口裂腹魚群體遺傳多樣性

    張爭世 胡冰潔 葉祥益 劉桂嘉 鄭曙明 葉 華

    (西南大學魚類繁育與健康養(yǎng)殖研究中心, 重慶 402460)

    研究測定了長江上游4個齊口裂腹魚(Schizothorax prenanti)野生群體(重慶巫溪、重慶城口、四川雅安、四川阿壩)共104個個體的線粒體Cyt b基因部分序列, 以探討齊口裂腹魚野生群體的遺傳多樣性和遺傳結(jié)構(gòu)。結(jié)果表明: 在104個個體Cyt b序列中共檢測到43個多態(tài)性位點, 25個單倍型。4個齊口裂腹魚群體的單倍型多樣性介于0.704—0.884, 核苷酸多樣性介于0.007—0.012。群體間Kimura雙參數(shù)遺傳距離介于0.008—0.017, 其中四川雅安群體與四川阿壩群體間遺傳距離最近, 基因交流頻繁。重慶城口群體與四川雅安群體間遺傳距離最遠, 基因交流受阻。AMOVA分析表明, 齊口裂腹魚的遺傳分化主要來自群體內(nèi)部, 且組群間、組群內(nèi)群體間和群體內(nèi)存在顯著的遺傳分化。中性檢驗得到Tajima’s D和 Fu’s Fs的值不顯著, 且歧點分布圖呈多峰, 表明長江上游4個齊口裂腹魚野生群體未經(jīng)歷過種群擴張。研究旨為齊口裂腹魚野生資源保護提供必要參考意見, 同時為齊口裂腹魚種質(zhì)資源合理開發(fā)和利用提供理論依據(jù)。

    齊口裂腹魚; 細胞色素b; 遺傳結(jié)構(gòu); 遺傳多樣性

    齊口裂腹魚(Schizothorax prenanti)俗稱雅魚、洋魚、齊口細鱗魚, 屬鯉形目、鯉科、裂腹魚亞科、裂腹魚屬, 主要分布于我國長江上游的金沙、岷江、大渡河、青衣江及烏江下游等水域, 在長江上游淡水生態(tài)系統(tǒng)的食物鏈中具有重要的地位[1,2]。近年來, 由于自然環(huán)境惡化及人為活動的不利影響,如過度捕撈、修建電站、水體污染等導致齊口裂腹魚野生資源遭到嚴重破壞, 其天然資源量逐漸減少, 種群數(shù)量急劇下降[3]。因此, 分析齊口裂腹魚不同群體的遺傳多樣性和遺傳結(jié)構(gòu), 對其資源的保護有著重要的理論價值和實際意義。

    迄今對齊口裂腹魚的研究多集中于生物學特性、營養(yǎng)與生長、繁殖和疾病防治等方面[4—7], 而對于其遺傳資源的調(diào)查、評估, 及種群結(jié)構(gòu)的研究較少, 研究者僅采用AFLP和mtDNA控制區(qū)標記分析了位于四川境內(nèi)的青衣江、大渡河、雅礱江上游的木里河及伍須湖齊口裂腹魚群體的遺傳多樣性[8—10], 而長江上游重慶段干流及支流的野生群體群體結(jié)構(gòu)情況等未見報道。魚類mtDNA具有分子較小、為共價閉合環(huán)狀雙鏈的DNA分子、嚴格的母系遺傳、快速的進化速度等特點[11,12], 因此被廣泛用于魚類親緣關(guān)系、種群遺傳分化、遺傳多樣性和種群歷史動態(tài)等研究[13—16]。在線粒體的眾多基因中, 細胞色素b (Cytochrome b, Cyt b)基因進化速度適中, 轉(zhuǎn)換或顛換比例較大, 被廣泛應用于系統(tǒng)發(fā)育、遺傳多樣性、物種鑒定研究中[17—20]。本研究以重慶和四川齊口裂腹魚野生群體為研究對象, 基于Cyt b分析齊口裂腹魚種群遺傳多樣性和種群遺傳結(jié)構(gòu), 以期為齊口裂腹魚種群資源的保護和合理開發(fā)利用提供理論依據(jù)。

    1 材料與方法

    1.1 實驗材料與DNA提取

    在2014年3月至2015年5月期間, 收集了重慶巫溪和城口、四川雅安和阿壩州4個齊口裂腹魚野生群體(表 1), 剪取鰭條用無水乙醇固定, 置于–20℃冰箱低溫保存?;蚪MDNA采用動物基因組DNA抽提試劑盒(生工生物工程股份有限公司, 上海)提取, 使用核酸蛋白定量儀(Eppendorf, BioPhotometer Plus, 德國)檢測所提DNA濃度與質(zhì)量, 并用1%瓊脂糖電泳檢測其完整性。DNA樣品稀釋成30 ng/mL的終濃度, 保存于–20℃?zhèn)溆谩?/p>

    表 1 齊口裂腹魚樣本采集信息Tab. 1 Sampling information of S. prenanti

    1.2 PCR擴增與測序

    根據(jù)齊口裂腹魚mtDNA全序列[21](GenBank序列號: NC_023829)采用Primer 5.0軟件設計Cyt b基因的引物, 并由英俊生物技術(shù)服務有限公司合成。上游引物和下游引物分別為: CYF: 5′-CGTGAACT ACGGCTGACT-3′; CYR: 5′-TAGGAACTGGGTG ATTGG-3′。

    PCR擴增體系為50 μL, 包括10×PCR buffer (with Mg2+) 5 μL、dNTPs (10 mmol/L) 4 μL、Taq酶(5 U/μL) 0.5 μL、上游及下游引物(10 μmol/L)各2.5 μL、模板DNA(30 ng/μL)2.5 μL, 補充無菌水至50 μL。所需試劑購自大連TaKaRa生物工程有限公司。PCR擴增程序為: 94℃預變性5min; 94℃變性30s, 55℃退火30s, 72℃延伸45s, 反應25個循環(huán); 72℃延伸10min, 15℃保存。擴增產(chǎn)物用1%瓊脂糖凝膠電泳檢測, 送上海英駿生物技術(shù)有限公司進行PCR擴增產(chǎn)物的純化、回收與單向測序, 測序引物為擴增引物。

    1.3 序列整理與數(shù)據(jù)分析

    采用BioEdit 7.0軟件[22]和Clustal X 1.81軟件[23]對測序結(jié)果進行編輯和同源比對。群體序列遺傳多樣性參數(shù), 包括多態(tài)性位點數(shù)(S)、單倍型數(shù)(H)、單倍型多樣性(Hd)、核苷酸多樣性(π)、平均核苷酸差異數(shù)(K)和Tajima’s D值等由DnaSP 5.0[24]軟件獲得。

    利用軟件Arlequin 3.5[25]進行核苷酸組成統(tǒng)計、分子方差分析(AMOVA)和遺傳分化指數(shù)(Φst)計算。使用Nm= (1–Φst)/(2Φst)計算兩兩群體間的基因流。利用Network 4.6.1.0軟件[26]構(gòu)建單倍型的簡約中介(Reduced-Median, MJ)網(wǎng)絡圖, 用以檢測單倍型之間的進化關(guān)系。利用MEGA 5.1軟件[27]計算群體間的Kimura雙參數(shù)模型(Kimura 2 Parameter, K2P)遺傳距離, 并構(gòu)建鄰接(Neighbor-Joinning, NJ)進化樹和最大簡約(Maximum Parsimony, MP)進化樹。使用Arlequin 3.5進行Tajima’s D檢驗[28]和Fu’s Fs檢驗[29], 同時結(jié)合核苷酸歧點分布分析(Mismatch distributions)來檢驗群體歷史動態(tài), 以確定是否存在瓶頸效應或群體擴張[30]。對顯著性檢驗的P值進行Bonferonni校正[31]。

    2 結(jié)果

    2.1 序列變異特征

    共測序獲得104尾齊口裂腹魚的Cyt b基因部分序列, 序列長750 bp。104條序列中共檢測到68個變異位點, 變異比例為9.86%。其中單一信息位點有33個, 簡約信息位點35個。變異均為轉(zhuǎn)換或顛換, 轉(zhuǎn)換與顛換的比值為54.33。4個種群所有個體序列中的A、T、C、G堿基組成比例為25.50%、29.70%、27.80%以及17.00%, A+T(55.20%)含量明顯高于G+C(44.80%)含量, 其中G含量最低(17.00%), 表現(xiàn)出明顯的反G偏倚特點。

    2.2 單倍型分布及遺傳關(guān)系

    圖5為標量脫靶量測量的參考坐標系,即脫靶量測量示意圖。測量系統(tǒng)位于坐標系原點,導彈以一定偏角和傾角與裝有脫靶量測量系統(tǒng)的靶標進行交匯。設交匯時間內(nèi),導彈相對于靶標作勻速直線運動,速度為v0,其運動軌跡與測量系統(tǒng)的最短距離,即脫靶量為ρ。記參考時間t =0時刻導彈與脫靶點的距離為D0,再記任意時刻導彈與測量系統(tǒng)的距離為R(t)。

    104個樣本中共檢測到25個單倍型, 其中Hap1、Hap8、Hap9為優(yōu)勢單倍型, 分別占個體數(shù)的20.18%、10.57%和14.42%; 單倍型Hap1、Hap2、Hap3、Hap7、Hap9和Hap15被2—3個群體共享, 其他19個單倍型為各群體特有單倍型, 其中, Hap2是分布最廣泛的單倍型, 除了城口群體, 其他3個群體都有分布(圖 1)。在4個群體中, 巫溪和雅安群體單倍型比例分別為43.75%和41.67%, 表現(xiàn)了較為豐富的單倍型數(shù)。而阿壩群體單倍型比例為24.14%, 表現(xiàn)了較為貧瘠的單倍型數(shù)。

    基于齊口裂腹魚單倍型構(gòu)建NJ分子系統(tǒng)樹(圖 1)可知, 25個單倍型被分成2個主支。其中一個主支包含14個單倍型, 有66個個體, 為優(yōu)勢分支; 另一個主支包含11個單倍型, 有48個個體。分子系統(tǒng)樹表明, 盡管不同群體的單倍型沒有嚴密的表現(xiàn)出與其地理位置分布的關(guān)系, 但是雅安和阿壩州群體之間,以及巫溪和城口群體之間表現(xiàn)出了較近的親緣關(guān)系。基于25個單倍型構(gòu)建的MP分子系統(tǒng)樹(圖 2)與MJ網(wǎng)絡圖(圖 3)顯示的結(jié)果相吻合。

    2.3 群體遺傳多樣性和遺傳分化

    圖 1 齊口裂腹魚25種Cyt b單倍型NJ系統(tǒng)樹及其在4個群體中的分布Fig. 1 NJ phylogenetic tree of 25 Cyt b haplotypes and its distribution of 4 S. prenanti populations

    圖 2 齊口裂腹魚25種Cyt b單倍型MP系統(tǒng)樹Fig. 2 MP phylogenetic tree of 25 Cyt b haplotypes of S. prenanti

    巫溪和雅安群體單倍型多樣性明顯高于城口和阿壩群體, 其中阿壩群體的值最低(0.704±0.079);從核苷酸多樣性和平均核苷酸差異數(shù)來看, 巫溪群體的值最高(0.012±0.001; 8.033), 阿壩群體的值最低(0.007±0.001; 4.621)。從以上結(jié)果可以看出(表 2),在4個群體中, 巫溪群體遺傳多樣性最高, 阿壩群體遺傳多樣性最低。

    從齊口裂腹魚4個地理群體的個體間遺傳距離來看(表 3), 群體內(nèi)遺傳距離大小次序為: 巫溪(0.012)>城口(0.010)>雅安(0.008)>阿壩州(0.007), 群體間的遺傳距離在0.008—0.017。其中, 雅安、阿壩群體與城口群體間遺傳距離較大, 雅安與阿壩群體之間遺傳距離最小。這說明, 雅安群體與城口群體親緣關(guān)系較遠, 雅安群體與阿壩群體親緣關(guān)系較近?;谌后w間K2P遺傳距離構(gòu)建UPGMA樹, 結(jié)果顯示,阿壩群體和雅安群體聚類為一支(四川組群), 巫溪群體和城口群體聚為另外一支(重慶組群)(圖 4)。

    圖 3 齊口裂腹魚群體Cyt b序列單倍型網(wǎng)絡結(jié)構(gòu)圖Fig. 3 Haplotypes network based on the different mutationsof Cyt b sequences of S. prenanti populations

    表 2 不同群體齊口裂腹魚遺傳多樣性參數(shù)Tab. 2 Genetic diversity index of four S. prenanti populations

    表 3 齊口裂腹魚群體內(nèi)遺傳距離(對角線, 粗體)及群體間遺傳距離(對角線下)和遺傳分化系數(shù)(Φst)(對角線上)Tab. 3 K2P genetic distances within population (diagonal, bold), and pairwise genetic distance (below diagonal) and fixation index (Φst) (above diagonal) among populations of S. prenanti populations

    齊口裂腹魚群體間的遺傳分化指數(shù)如表 3(對角線上)所示。4個群體遺傳分化指數(shù)在0.054—0.468。根據(jù)Wright關(guān)于群體遺傳分化程度的理論[32],雅安與阿壩群體之間的遺傳分化指數(shù)最小(Φst=0.054, P>0.05), 其值在0.05—0.15, 屬于低度分化, 表明雅安和阿壩群體間遺傳分化程度較低。巫溪與城口群體之間的遺傳分化指數(shù)(Φst=0.228, P=0)在0.15—0.25, 屬于中度分化。城口群體與四川組群(雅安、阿壩群體)之間, 以及巫溪與四川組群(雅安、阿壩群體)之間遺傳分化指數(shù)都大于0.25,屬于高度分化。其中, 城口和阿壩群體之間遺傳分化指數(shù)最高(Φst=0.468, P=0)。

    在AMOVA分析中(表 4), 將4個群體歸為一個組群分析表明, 長江上游齊口裂腹魚群體內(nèi)的變異組分(64.23%)大于群體間的變異組分(35.77%), 各群體間存在著顯著遺傳分化(Fst=0.357, P=0)?;谶z傳分化分析結(jié)果, 將4個群體分成四川和重慶2個組群, 結(jié)果表明, 31.50%的變異組分來自2個組群之間, 兩組群間存在顯著的遺傳分化(Fct=0.315, P=0), 10.69%的變異組分來自組群內(nèi)群體間(Fsc=0.156, P=0), 而各個群體內(nèi)個體間的變異組分達到57.81% (Fst=0.422, P=0), 由此說明, 長江上游齊口裂腹魚的變異主要來自群體內(nèi)個體間, 不同組群間及組群內(nèi)群體間存在顯著遺傳分化。

    圖 4 基于Kimura遺傳距離的齊口裂腹魚4個群體之間的UPGMA樹Fig. 4 UPMA tree of 4 S. prenanti populations based on K-2-P genetic distance

    表 4 齊口裂腹魚群體Cyt b序列變異的分子方差分析(AMOVA)Tab. 4 Analysis of molecular variances (AMOVA) of Cyt b sequences of S. prenanti populations

    對4個群體間基因流Nm分析可知(表 5), 巫溪群體與城口群體, 雅安群體與阿壩群體之間的Nm值均大于1, 說明重慶組群以及四川組群內(nèi)都有基因交流。其中, 雅安群體與阿壩群體間Nm值高達8.791,表明兩群體間存在較為頻繁的基因交流。而城口群體與四川組群(雅安、阿壩群體)之間的Nm值均小于1, 巫溪群體與四川組群(雅安、阿壩群體)之間的Nm值略大于1, 表明四川組群和重慶組群之間基因交流受阻。

    2.4 群體歷史動態(tài)分析

    對所有個體進行中性檢驗(表 6)和歧點分布分析, 結(jié)果顯示4個群體的Fu’s Fs檢驗及Tajima’s D檢驗結(jié)果都不顯著(P>0.05), 且其歧點分布圖均呈明顯的多峰形(圖 5), 表明4個齊口裂腹魚野生群體在歷史進化的過程中未經(jīng)歷過種群擴張。

    3 討論

    3.1 Cyt b序列變異特征

    本研究對長江上游4個野生群體齊口裂腹魚Cyt b序列分析表明, A+T達到55.20%, 明顯量高于G+C含量, 這與晁燕等[33]在齊口裂腹魚中發(fā)現(xiàn)的54.9% A+T含量相符。這可能與生物的進化地位有關(guān), 即線粒體DNA的進化地位隨著A+T含量的增加而增高[34]。同時, 本研究低G(16.90%)含量與Qi等[20]和謝佳燕等[35]對長江干流及青藏高原的齊口裂腹魚的研究結(jié)果一致, 表明Cyt b有很強的堿基組成偏向性。

    表 5 齊口裂腹魚群體間的基因流Nm(左下角)Tab. 5 Gene flow value (Nm, below diagonal) between every two populations of S. prenanti

    表 6 齊口裂腹魚群體的Fu’S Fs和Tajima’s D中性檢驗Tab. 6 Fu’S Fsand Tajima’s D of S. prenanti populations

    圖 5 齊口裂腹魚群體歧點分布圖Fig. 5 Mismatch distributions of S. prenanti populations實線為理論的期望分布, 虛線為觀測的分布The solid line represents theoretical expected distribution, and the dashed line represents observed distribution

    研究發(fā)現(xiàn)長江上游104個齊口裂腹魚Cyt b序列有68個變異位點, 共檢測到單倍型數(shù)為25個。單倍型數(shù)占全部個體數(shù)的比例為24.04%, 高于孟瑋等[36]對塔里木裂腹魚的研究結(jié)果(10.53%), 低于Qi等[20]對長江干流齊口裂腹魚的研究結(jié)果(53.85%)。與謝佳燕等[35]對青藏高原齊口裂腹魚的研究結(jié)果相似(29.63%)。25個單倍型中Hap2分布最廣泛, 存在于除重慶城口以外的其他3個群體中。長江上游4個野生齊口裂腹魚群體獨享19個單倍型, 遠高于晁燕等[33]的研究結(jié)果(2個), 說明長江上游四川、重慶江段干流及支流齊口裂腹魚群體單倍型較青海班瑪縣瑪柯河的豐富, 這些群體更能夠有代表性的說明齊口裂腹魚群體的遺傳多樣性。

    3.2 群體遺傳多樣性

    單倍型多樣性指數(shù)(Hd)和核苷酸多樣性指數(shù)(π)是評估遺傳多樣性的兩個重要參數(shù)。本研究中的齊口裂腹魚群體多樣性指數(shù)(Hd=0.704—0.884, π=0.007—0.012)高于青?,斂潞雍退拇ǖぐ偷凝R口裂腹魚群體(Hd=0.56—0.83, π=0.0007—0.0032)[36], 低于木里河、大渡河和青衣江群體(Hd= 0.93—0.99, π=0.014—0.02)[9]; 與其他裂腹魚亞科魚類相比, 高于塔里木裂腹魚(Hd=0.19—0.74, π= 0.0004—0.0029)[34], 和黃河和托索湖的黃河裸裂尻(Hd=0.67—0.70, π=0.0012—0.0026)[37], 略低于柴達木內(nèi)流水系和4個黃河支流的黃河裸裂尻(Hd= 0.79—0.97, π=0.0027—0.0099)[38]。綜上可知, 齊口裂腹魚遺傳多樣性處于中等水平。在4個群體中, 巫溪和雅安群體的Hd和π值較高, 巫溪群體群體內(nèi)遺傳距離最高, 遺傳變異最豐富, 遺傳多樣性最高, 而阿壩群體遺傳變異最小, 多樣性最低。重慶巫溪群體遺傳多樣性相對較高, 一方面可能與重慶巫溪有著較為適宜的水文條件和優(yōu)越的生存環(huán)境有關(guān); 另一方面可能是巫溪地勢險峻, 人為活動對齊口裂腹魚野生種質(zhì)資源的影響較小。四川阿壩群體遺傳多樣性較低, 這可能與近年來自然環(huán)境的惡化及以及棲息水域受到污染、繁育環(huán)境發(fā)生改變、過度捕撈等人為活動的影響有關(guān)。

    利用中性檢驗和核苷酸歧點分布研究物種的種群歷史, 若Fu’s Fs的檢測呈顯著負值且歧點分布分析圖譜呈泊松狀分布的單峰, 說明種群趨于發(fā)生了瓶頸效應或者種群擴張[39]。本研究Fu’s Fs和Tajima’s D均不顯著(P>0.05), 且歧點分布圖均呈明顯的多峰形, 表明長江上游4個野生齊口裂腹魚群體在歷史進化的過程中相對較為穩(wěn)定, 未經(jīng)歷過種群擴張。

    3.3 群體遺傳結(jié)構(gòu)

    研究群體遺傳結(jié)構(gòu)不僅可以用于評價物種群體的變異水平和不同地理群體之間的關(guān)系, 還可以確定群體中的進化顯著單元和管理單元, 以及制定資源的保護和管理策略[40]。在本研究中, 除雅安與阿壩群體其余群體間遺傳分化均達極顯著水平, 但群體間最大遺傳距離(0.017)仍小于亞種的界限(0.02), 其均未達到亞種分化[41,42]。根據(jù)Wright[32]關(guān)于對群體遺傳分化和基因交流程度的理論, 4個野生群體中, 城口和雅安群體(Φst=0.458, P=0, Nm= 0.592)、城口和阿壩群體(Φst=0.468, P=0, Nm= 0.567)之間基因交流受阻, 群體間遺傳分化達極顯著水平, 這可能與分布流域較遠和地理隔離等因素有關(guān)。阿壩和雅安群體(Φst=0.054, P>0.05 Nm= 8.791)基因交流頻繁, 群體間遺傳分化不顯著。同時, 阿壩和雅安群體之間共享多個單倍型, 有較小的遺傳距離, 可能是因為腳木足河是大渡河的支流,地理分布特點決定了兩群體間可頻繁的進行基因交流。

    本研究對齊口裂腹魚主要分布區(qū)域, 包括長江上游四川、重慶江段干流及支流齊口裂腹魚種群進行了多樣性分析, 4個群體遺傳多樣性整體不高,應該加強齊口裂腹魚各群體的資源保護力度。由于四川組群與重慶組群之間存在顯著的遺傳分化,建議將四川群體、重慶群體分別作為2個進化顯著單元(Evolutionarily Significant Unit, ESU)進行保護。本研究初步反映了齊口裂腹魚的遺傳背景、遺傳結(jié)構(gòu)及遺傳變異水平, 后續(xù)可以通過利用多種分子標記技術(shù)如SSR、SNP等聯(lián)合分析, 更全面地為齊口裂腹魚資源的保護和合理開發(fā)提供理論依據(jù)。

    [1]Wu Y F, Wu C Z. The Qinghai-Tibet Plateau Fish [M]. Chengdu: Sichuan Science and Technology Press. 1992, 361—365 [武云飛, 吳翠珍. 青藏高原魚類. 成都: 四川科學技術(shù)出版社. 1992, 361—365]

    [2]Chen Y F, Cao W X. Schizothoracinae. In: Le P Q, Shan X H, Lin R D, et al. Chinese Animal: Osteichthyes Carp Shape Mesh Volume [M]. Beijing: Science and Technology Press. 2000, 273—390 [陳毅峰, 曹文宣. 裂腹魚亞科. 樂佩琦, 單鄉(xiāng)紅, 林人端, 等. 中國動物志?硬骨魚綱?鯉形目(下). 北京: 科學出版社. 2000, 273—390]

    [3]Liu S P, Duan X B, Chen D Q, et al. Studies on status of fishery resources in the middle reach of the Yangtze River [J]. Acta Hydrobiologica Sinica, 2005, 29(6): 708—711 [劉紹平, 段辛斌, 陳大慶, 等. 長江中游漁業(yè)資源現(xiàn)狀研究. 水生生物學報, 2005, 29(6): 708—711]

    [4]Ruo M, W H T, Y Q Y, et al. Artificial breeding research of Schizothorax prenanti [J]. Freshwater Fisheries, 2001, 31(6): 3—5 [若木, 王鴻泰, 殷啟云, 等. 齊口裂腹魚人工繁殖的研究. 淡水漁業(yè), 2001, 31(6): 3—5]

    [5]Wang T, Zhou C W, Yuan D Y, et al. Schizothorax prenanti corticotropin-releasing hormone (CRH): molecular cloning, tissue expression, and the function of feeding regulation [J]. Fish Physiology & Biochemistry, 2014, 40(5): 1407—1415

    [6]Lin F J, Zhou C W, Chen H, et al. Molecular characterization, tissue distribution and feeding related changes of NUCB2A/nesfatin-1 in Ya-fish (Schizothorax prenanti) [J]. Gene, 2014, 536(2): 238—246

    [7]Zhang L, Wu Y L, Wang L, et al. Effects of Oxidized Konjac glucomannan (OKGM) on growth and immune function of Schizothorax prenanti [J]. Fish & Shellfish Immunology, 2013, 35(4): 1105—1110

    [8]Song J, Song Z B, Yue B S, et al. Assessing genetic diversity of wild populations of Prenant’s schizothoracin, Schizothorax prenanti, using AFLP markers [J]. Environmental Biology of Fishes, 2006, 77(1): 79—86

    [9]Song Z B, Song J, Yue B S. Population genetic diversity of Prenant’s schizothoracin Schizothorax prenanti, inferred from the mitochondrial DNA control region [J].Environmental Biology of Fishes, 2008, 81(3): 247—252

    [10]Liang J J, Liu Y, Zhang X F. An observation of the loss of genetic variability in Schizothorax prenanti, Schizothorax prenanti, inhabiting a plateau lake [J]. Biochemical Systematics and Ecology, 2011, 39(4—6): 361—370

    [11]Meland S, Johansen S, Johansen T, et al. Rapid disappearance of one parental mitochondrial genotype alfter isogamous mating in the myxomycete Physarum polycephalum [J]. Current Genetics, 1991, 19(1): 55—60

    [12]Zhao K. Mitochondrial DNA (mtDNA) of fish and application in molecular phylogenetic [J]. Journal of Qinghai University (Nature Science), 2006, 24(4): 49—53 [趙凱.魚類線粒體DNA 及其在分子系統(tǒng)學中的應用. 青海大學學報: 自然科學版, 2006, 24(4): 49—53]

    [13]Lu X, Wang H, Liu B, et al. An effective method for parentage determination of the clam (Meretrix meretrix) based on SSR and COI markers [J]. Aquaculture, 2011, 318(1—2): 223—228

    [14]Mabuchi K, Miya M, Senou H, et al. Complete mitochondrial DNA sequence of the Lake Biwa wild strain of common carp (Cyprinus carpio L.): further evidence for an ancient origin [J]. Aquaculture, 2006, 257(1—4): 67—77

    [15]Xiao W H, Zhang Y P. Genetics and evolution of mitochondrial DNA in fish [J]. Acta Hydrobiologica Sinica, 2000, 24(4): 384—391 [肖武漢, 張亞平. 魚類線粒體DNA的遺傳與進化. 水生生物學報, 2000, 24(4): 384—391]

    [16]Guo X H, Liu S J, Liu Q, et al. New progresses on mitochondrial DNA in fish [J]. Acta Genetica Sinica, 2004, 31(9): 983—1000 [郭新紅, 劉少軍, 劉巧, 等. 魚類線粒體DNA研究新進展. 遺傳學報, 2004, 31(9): 983—1000]

    [17]Hu J, Hou X Y, Yin S W, et al. Genetic diversity and divergence of Cheilinus undulates of different geographic populations of the south China sea revealed by COI and CYT B gene analyses [J]. Acta Hydrobiologica Sinica, 2014, 38(6): 1008—1016 [胡靜, 侯新遠, 尹紹武, 等. 基于mtDNA COI和Cyt b基因序列對中國南海不同海域波紋唇魚群體遺傳多樣性的研究. 水生生物學報, 2014, 38(6): 1008—1016]

    [18]Mcveigh H P, Davidson W S. A salmonid phylogeny inferred from mitochondrial cytochrome b gene sequences [J]. Journal of Fish Biology, 2006, 39(sA): 277—282

    [19]Durand J D, Tsigenopoulos C S, Unlu E, et al. Phylogeny and biogeography of the family cyprinidae in the middle east inferred from cytochrome b DNA- evolutionary significance of this region [J]. Molecular Phylogenetics & Evolution, 2002, 22(1): 91—100

    [20]Qi D L, Li T P, Zhao X Q, et al. Mitochondrial cytochrome b sequence variation and phylogenetics of the highly specialized Schizothoracine fishes (Teleostei: Cyprinidae) in the Qinghai-Tibet Plateau [J]. Biochemical Genetics, 2006, 744(5): 270—285

    [21]Chen C N, Huang Y Y, Li H, et al. The complete mitochondrial genome of Schizothorax prenanti (Tchang) (Teleostei, Cyprinidae, Schizothoracinae) [J]. Mitochondrial DNA, 2016, 27(1): 253—254

    [22]Hall T A. BioEdit: a user-friendly biological sequence alignment editor and analysis program for windows 95/98/NT [J]. Nucleic Acids Symposium Series, 1998, 41: 95—98

    [23]Thompson J D, Gibson T J, Plewniak F, et al. The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools [J]. Nucleic Acids Research, 1997, 25(24): 4876—4882

    [24]Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data [J]. Bioinformatics, 2009, 25(11): 1451—1452(2)

    [25]Excoffier L, Lischer H E L. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows [J]. Molecular Ecology Resources, 2010, 10(3): 564—567

    [26]Bandelt H J, Forster P, R?hl A A. Median-joining network for inferring intraspecific phylogenies [J]. Molecular Biology & Evolution, 1999, 16(1): 37—48

    [27]Tamura K, Peterson D, Peterson N, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods [J]. Molecular Biology and Evolution, 2011, 28(10): 2731—2739

    [28]Hickerson M J, Meyer C P. Testing comparative phylogeographic models of marine vicariance and dispersal using a hierarchical Bayesian approach [J]. BMC Evolutionary Biology, 2007, 8(3): 322—340

    [29]Fu Y X. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection [J]. Genetics, 1997, 147(2): 915—925

    [30]Aris-Brosou S, Excoffier L. The impact of population expansion and mutation rate heterogeneity on DNA sequence polymorphism [J]. Molecular Biology and Evolution, 1996, 13(3): 494—504

    [31]Rice W R. Analyzing tables of statistical tests [J]. Evolution, 1989, 43(1): 223—225

    [32]Wright S. Evolution in Mendelian populations [J]. Bulletin of Mathematical Biology, 1990, 52(1—2): 241—295

    [33]Chao Y, Yang C, Shen Z X, et al. Genetic Diversity of Schizothorax prenanti inferred from mitochondrial cytochrome b gene sequences [J]. Fisheries Science, 2010, 29(7): 421—424 [晁燕, 楊成, 申志新, 等. 齊口裂腹魚基于Cyt b基因序列多態(tài)性研究. 水產(chǎn)科學, 2010, 29(7): 421—424]

    [34]Folmer O, Black M, Hoeh W, et al. DNA primer for amplification of mitochondrial cytochrome coxidase subunit I from diverse metazoan invertebrates [J]. Molecular Marine Biotechnology, 1994, 3(5): 294—299

    [35]Xie J Y. Structure analysis of mtDNA control region in Schizothorax prenanti [J]. Chinese Journal of Zoology, 2011, 46(2): 97—101 [謝佳燕. 齊口裂腹魚線粒體DNA控制區(qū)結(jié)構(gòu)分析. 動物學雜志, 2011, 46(2): 97—101]

    [36]Meng W, Guo Y, Hai S, et al. Genetic structure and diversity of Schizothorax bidduphi populations [J]. Acta Hydrobiologica Sinica, 2012, 36(5): 851—857 [孟瑋, 郭焱, 海薩, 等. 塔里木裂腹魚群體遺傳結(jié)構(gòu)及遺傳多樣性分析. 水生生物學報, 2012, 36(5): 851—857]

    [37]Zhao K, Yang G S, Li J B, et al. Phylogenetic structure of Schizopygopsis pylzovi populations from mitochondrial B gene sequence variations [J]. Acta Hydrobiologica Sinica, 2006, 30(2): 129—133 [趙凱, 楊公社, 李俊冰, 等. 黃河裸裂尻魚群體遺傳結(jié)構(gòu)和Cyt b序列變異. 水生生物學報, 2006, 30(2): 129—133]

    [38]Qi D L, Chao Y, Guo S, et al. Genetic structure of five Huanghe schizothoracin Schizopygopsis pylzovi populations based on mtDNA control region sequences [J]. Acta Zoologica Sinica, 2008, 54(6): 972—980 [祁得林, 晁燕,郭松, 等. 黃河裸裂尻魚五種群mtDNA控制區(qū)的遺傳結(jié)構(gòu). 動物學報, 2008, 54(6): 972—980]

    [39]Barbosa A M, Real R, Mu?oz A R, et al. New measures for assessing model equilibrium and prediction mismatch in species distribution models [J]. Diversity and Distributions, (Diversity Distrib.) 2013, 19(10): 1333—1338

    [40]O’Brien S J. A role for molecular genetics in biological conservation [J]. Proceedings of the National Academy of Sciences of the United States of America, 1994, 91(13): 5748—5755

    [41]Shaklee J B, Tamaru C S, Waples R S. Speciation and evolution of marine fishes studied by the electrophoretic analysis of proteins [J]. Pacificence, 1982, 36(2): 141—157

    [42]Nei, Masatoshi. Molecular Population Genetics and Evolution [M]. American Elsevier Publishing Company. 1975, 128—132

    GENETIC DIVERSITY OF THE PRENANT’S SCHIZOTHORACIN (SCHIZOTHORAX PRENANTI) BASED ON PARTIAL MTDNA CYT B SEQUENCES

    ZHANG Zheng-Shi, HU Bing-Jie, YE Xiang-Yi, LIU Gui-Jia, ZHENG Shu-Ming and YE Hua
    (Department of Fisheries in Rongchang District Southwest University, Chongqing 402460, China)

    Prenant’s schizothoracin, Schizothorax prenanti (Tchang), a valuable and economic fish, distributes in the main streams and tributaries of the upper reaches of the Yangtze River in China. To assess genetic diversity and genetic structure of the wild populations of S. prenanti, we amplified partial Mitochondrial Cyt b gene sequences from 104 individuals of the four wild populations of S. pregnant (Wuxi, Chengkou, Ya’an and A’ba) from the upper reaches of the Yangtze River. The results showed that 43 polymorphic nucleotide sites and 25 haplotypes were detected among the Cyt b sequences of 104 individuals. The haplotype diversity was ranged from 0.704 to 0.884, and nucleotide diversity index was ranged from 0.007 to 0.012 among the four populations. Kimura 2-parameter distance was ranged from 0.008 and 0.017. The genetic distance between Ya’an and A’ba population was the smallest and the gene exchange between the two populations was frequent. The largest genetic distance was observed between Ya’an and Chengkou population, which might be explained by the little gene exchange between the two populations. AMOVA analysis showed that the genetic differentiation of S. prenanti was mainly occurred within populations variation. There were significant genetic differentiations between the groups, among populations within groups, and within populations. The results of neutrality test showed that the values of Tajima’s D and Fu’s Fswere not significant and the mismatch distribution was multimodal, indicating that the four populations of S. prenanti from the upper reaches of the Yangtze River had not experienced population expansion. This study provides valuable knowledge for protecting Schizothorax prenanti and theoretic information for a rational utilization of the germplasm resources.

    Schizothorax prenanti; Cyt b; Genetic structure; Genetic diversity

    Q347

    A

    1000-3207(2017)03-0609-08

    10.7541/2017.78

    2016-05-16;

    2016-11-01

    國家自然科學基金(31402302); 中央高?;究蒲袠I(yè)務費項目(XDJK2014C152)資助 [Supported by the National Natural Science Foundation of China (31402302); Fundamental Research Funds for the Central Universities (XDJK2014C152)]

    張爭世(1989—), 男, 山東臨沂人; 碩士研究生; 研究方向為水產(chǎn)動物遺傳育種。E-mail: zhzs0906@163.com

    鄭曙明(1957—), 女, 江蘇南京人; 教授, 碩導; 主要從事水產(chǎn)動物遺傳育種與病害防治技術(shù)研究。E-mail: zhsm199@163.com;葉華(1981—), 女, 四川廣元人; 副教授, 碩導; 主要從事水產(chǎn)動物遺傳育種與生物技術(shù)研究。E-mail: yhlh2000@126.com

    猜你喜歡
    腹魚巫溪組群
    人工養(yǎng)殖2 種裂腹魚生物學指數(shù)和肌肉營養(yǎng)價值的比較分析
    下一站,巫溪
    城市地理(2023年9期)2023-12-02 05:55:30
    齊口裂腹魚集群行為對流態(tài)的響應
    南疆地區(qū)四種土著裂腹魚類幼魚對氯化鈉和碳酸氫鈉的急性耐受能力
    在十八梯遇見巫溪
    今日重慶(2021年11期)2021-12-21 21:33:15
    73個傳統(tǒng)建筑組群組團出道!帶你活進從前的慢時光
    “組群”“妙比”“知人”:小學語文古詩群文閱讀的三個途徑
    “巫溪洋芋”拿到出口歐盟“通行證”
    重慶與世界(2020年9期)2020-10-20 06:27:03
    磁盤組群組及iSCSI Target設置
    齊口裂腹魚腸道氣單胞菌的分離鑒定
    国产精品秋霞免费鲁丝片| 丰满乱子伦码专区| 亚洲精品日本国产第一区| 精品99又大又爽又粗少妇毛片| 99久久中文字幕三级久久日本| 成人漫画全彩无遮挡| 中文乱码字字幕精品一区二区三区| 国产欧美日韩一区二区三区在线 | 国产高清有码在线观看视频| 久久99热6这里只有精品| 在线观看免费视频网站a站| 国产黄片美女视频| 国产精品av视频在线免费观看| 久久人人爽av亚洲精品天堂 | 纵有疾风起免费观看全集完整版| av免费在线看不卡| 日本一二三区视频观看| 国产高清三级在线| 成人亚洲欧美一区二区av| 婷婷色综合大香蕉| 女人十人毛片免费观看3o分钟| 麻豆乱淫一区二区| 99热这里只有是精品50| 夜夜看夜夜爽夜夜摸| 免费播放大片免费观看视频在线观看| 在线观看国产h片| 精品一区二区三区视频在线| 男人添女人高潮全过程视频| .国产精品久久| 全区人妻精品视频| 欧美一级a爱片免费观看看| 黄色日韩在线| 一本色道久久久久久精品综合| 在线观看免费高清a一片| 国产一级毛片在线| 伦理电影大哥的女人| 亚洲av欧美aⅴ国产| 婷婷色综合大香蕉| 午夜福利视频精品| 特大巨黑吊av在线直播| 中文字幕久久专区| 肉色欧美久久久久久久蜜桃| 一本色道久久久久久精品综合| 国精品久久久久久国模美| 建设人人有责人人尽责人人享有的 | 日韩一本色道免费dvd| 人妻系列 视频| 免费黄色在线免费观看| 最近2019中文字幕mv第一页| 久久久久精品性色| 黄色视频在线播放观看不卡| 蜜桃亚洲精品一区二区三区| 18禁动态无遮挡网站| 一级爰片在线观看| 99re6热这里在线精品视频| 在线天堂最新版资源| 亚洲av国产av综合av卡| 亚洲av成人精品一二三区| 少妇的逼好多水| 色视频在线一区二区三区| 丰满少妇做爰视频| 国产黄色视频一区二区在线观看| 内地一区二区视频在线| 精品酒店卫生间| 最黄视频免费看| 一级片'在线观看视频| 高清不卡的av网站| 欧美xxxx黑人xx丫x性爽| 日本欧美国产在线视频| 天天躁日日操中文字幕| 国产av国产精品国产| 老司机影院成人| 国产综合精华液| 国产亚洲精品久久久com| 国产中年淑女户外野战色| 少妇的逼水好多| 久久精品人妻少妇| 综合色丁香网| 女的被弄到高潮叫床怎么办| 成人亚洲精品一区在线观看 | 欧美xxxx黑人xx丫x性爽| av黄色大香蕉| 亚洲av欧美aⅴ国产| 欧美国产精品一级二级三级 | 国产欧美日韩精品一区二区| 久久热精品热| 久久青草综合色| 欧美极品一区二区三区四区| 国产熟女欧美一区二区| 欧美日韩精品成人综合77777| 久久精品久久久久久噜噜老黄| 欧美国产精品一级二级三级 | 国产免费又黄又爽又色| 午夜福利高清视频| 青春草亚洲视频在线观看| 国产视频内射| 亚洲av二区三区四区| 少妇人妻久久综合中文| 久久久久网色| 亚洲性久久影院| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 国产成人精品一,二区| 日韩欧美 国产精品| 多毛熟女@视频| 99热6这里只有精品| 欧美日本视频| 成人国产麻豆网| 国产免费一区二区三区四区乱码| 成年女人在线观看亚洲视频| 国产一区二区三区av在线| av国产精品久久久久影院| 九九爱精品视频在线观看| 夫妻性生交免费视频一级片| 色哟哟·www| 日韩一区二区三区影片| 黄色日韩在线| 亚洲熟女精品中文字幕| 成人无遮挡网站| 又黄又爽又刺激的免费视频.| 制服丝袜香蕉在线| 久久久精品94久久精品| 国产无遮挡羞羞视频在线观看| 国产欧美亚洲国产| 精品午夜福利在线看| 人妻少妇偷人精品九色| 成人高潮视频无遮挡免费网站| 26uuu在线亚洲综合色| 成人无遮挡网站| 国产免费又黄又爽又色| 少妇人妻一区二区三区视频| 毛片一级片免费看久久久久| 人妻 亚洲 视频| 日本与韩国留学比较| 亚洲,一卡二卡三卡| 精品亚洲乱码少妇综合久久| 身体一侧抽搐| 日本爱情动作片www.在线观看| 熟妇人妻不卡中文字幕| 亚洲av欧美aⅴ国产| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 精品少妇黑人巨大在线播放| 国产人妻一区二区三区在| 一级毛片aaaaaa免费看小| 久久99精品国语久久久| 高清在线视频一区二区三区| 美女视频免费永久观看网站| 最近中文字幕2019免费版| 啦啦啦在线观看免费高清www| 水蜜桃什么品种好| 国产探花极品一区二区| 国产一区亚洲一区在线观看| 国产69精品久久久久777片| 欧美激情极品国产一区二区三区 | 激情 狠狠 欧美| 午夜福利视频精品| 久久精品国产亚洲网站| av线在线观看网站| 赤兔流量卡办理| 国产精品一区二区性色av| 草草在线视频免费看| 亚洲人成网站在线观看播放| 日韩在线高清观看一区二区三区| 亚洲欧美一区二区三区黑人 | 免费大片黄手机在线观看| 久久精品人妻少妇| 一级a做视频免费观看| 国产成人aa在线观看| 国产成人freesex在线| av免费在线看不卡| 蜜桃亚洲精品一区二区三区| 亚洲最大成人中文| 这个男人来自地球电影免费观看 | 国产免费福利视频在线观看| 黄色怎么调成土黄色| 国产在线一区二区三区精| 国产女主播在线喷水免费视频网站| 99九九线精品视频在线观看视频| 亚洲av成人精品一区久久| 亚洲一区二区三区欧美精品| 能在线免费看毛片的网站| av专区在线播放| 少妇人妻一区二区三区视频| 丝袜脚勾引网站| 欧美激情国产日韩精品一区| 免费久久久久久久精品成人欧美视频 | 老师上课跳d突然被开到最大视频| av网站免费在线观看视频| 亚洲精品久久午夜乱码| 国产淫语在线视频| 亚洲国产av新网站| 性色av一级| 国产男女超爽视频在线观看| 一级毛片久久久久久久久女| 成人一区二区视频在线观看| 五月伊人婷婷丁香| 狂野欧美激情性bbbbbb| 在线天堂最新版资源| 久久久久人妻精品一区果冻| 国产成人a区在线观看| av播播在线观看一区| 午夜福利影视在线免费观看| 日韩,欧美,国产一区二区三区| 中国三级夫妇交换| 国产精品蜜桃在线观看| 久久精品人妻少妇| 99久国产av精品国产电影| 人人妻人人爽人人添夜夜欢视频 | 一边亲一边摸免费视频| 丝瓜视频免费看黄片| 久久久久久九九精品二区国产| 久久97久久精品| 美女中出高潮动态图| 欧美日韩精品成人综合77777| 国产老妇伦熟女老妇高清| 一个人看视频在线观看www免费| www.色视频.com| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费 | 国产精品人妻久久久影院| 啦啦啦啦在线视频资源| 成人毛片a级毛片在线播放| 啦啦啦中文免费视频观看日本| 国产日韩欧美在线精品| 99热国产这里只有精品6| 精品人妻熟女av久视频| 精品熟女少妇av免费看| 午夜福利视频精品| 国产探花极品一区二区| 噜噜噜噜噜久久久久久91| 新久久久久国产一级毛片| 日产精品乱码卡一卡2卡三| 最后的刺客免费高清国语| 欧美老熟妇乱子伦牲交| 国产亚洲91精品色在线| 欧美成人午夜免费资源| 国产精品无大码| 亚洲精品日本国产第一区| 国产在线视频一区二区| 国产精品一区二区在线不卡| 久久婷婷青草| 爱豆传媒免费全集在线观看| 国产精品一二三区在线看| 全区人妻精品视频| 高清不卡的av网站| 熟女av电影| 一本色道久久久久久精品综合| 美女视频免费永久观看网站| 免费高清在线观看视频在线观看| 人妻少妇偷人精品九色| 永久网站在线| 日日摸夜夜添夜夜爱| 在线观看一区二区三区激情| 久久久成人免费电影| 99热全是精品| 亚洲欧美一区二区三区黑人 | 好男人视频免费观看在线| 国产精品人妻久久久久久| 久久久久久伊人网av| 国产精品久久久久久精品古装| 高清黄色对白视频在线免费看 | 啦啦啦在线观看免费高清www| 99热这里只有精品一区| 精品人妻视频免费看| 亚洲国产精品成人久久小说| 亚洲欧美一区二区三区黑人 | 青春草国产在线视频| av福利片在线观看| 亚洲精品国产色婷婷电影| 在线观看人妻少妇| 亚洲天堂av无毛| 最近2019中文字幕mv第一页| 日本黄大片高清| 亚洲色图综合在线观看| 亚洲久久久国产精品| 精品久久久久久久久av| 亚洲三级黄色毛片| 国产精品国产三级国产av玫瑰| 国产v大片淫在线免费观看| av网站免费在线观看视频| 日韩欧美精品免费久久| 看十八女毛片水多多多| av在线播放精品| 成人一区二区视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | av福利片在线观看| 卡戴珊不雅视频在线播放| 蜜桃在线观看..| 日日撸夜夜添| 美女高潮的动态| 亚洲色图av天堂| 国产精品一及| 99久久精品国产国产毛片| 国产日韩欧美亚洲二区| 国产男女超爽视频在线观看| 精品午夜福利在线看| 亚洲欧美日韩另类电影网站 | 蜜臀久久99精品久久宅男| 免费观看性生交大片5| 久久久久久久久久久丰满| 人妻夜夜爽99麻豆av| 青春草亚洲视频在线观看| 内地一区二区视频在线| 天堂中文最新版在线下载| 久久久精品94久久精品| 香蕉精品网在线| 欧美日韩一区二区视频在线观看视频在线| 国产精品欧美亚洲77777| 干丝袜人妻中文字幕| 黑丝袜美女国产一区| 2018国产大陆天天弄谢| 精品人妻熟女av久视频| 日韩欧美 国产精品| 亚洲精品国产av成人精品| 亚洲中文av在线| 欧美性感艳星| 日韩成人伦理影院| 成人18禁高潮啪啪吃奶动态图 | 性色avwww在线观看| 亚洲最大成人中文| 黄色怎么调成土黄色| 欧美老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 国产在线一区二区三区精| 高清黄色对白视频在线免费看 | 亚洲成人手机| 国产精品一区二区三区四区免费观看| 嫩草影院新地址| 免费看不卡的av| 国产一区亚洲一区在线观看| 国产在线一区二区三区精| 亚洲精品,欧美精品| 久久精品夜色国产| 亚洲精品视频女| av女优亚洲男人天堂| 美女脱内裤让男人舔精品视频| 卡戴珊不雅视频在线播放| 不卡视频在线观看欧美| 久久久久人妻精品一区果冻| 小蜜桃在线观看免费完整版高清| 色哟哟·www| 精品一区二区三卡| 免费久久久久久久精品成人欧美视频 | 综合色丁香网| 亚洲不卡免费看| 亚洲av欧美aⅴ国产| 久久久成人免费电影| 九九久久精品国产亚洲av麻豆| 亚洲av成人精品一二三区| 国产美女午夜福利| 亚洲色图av天堂| 国产精品三级大全| 免费看日本二区| 免费观看av网站的网址| 国产亚洲午夜精品一区二区久久| 免费看av在线观看网站| 80岁老熟妇乱子伦牲交| 免费观看a级毛片全部| 啦啦啦视频在线资源免费观看| 日韩电影二区| 亚洲精品自拍成人| 91精品国产国语对白视频| 国产亚洲av片在线观看秒播厂| 精品人妻偷拍中文字幕| 精品一区在线观看国产| 蜜桃久久精品国产亚洲av| 国产永久视频网站| 亚洲精品久久午夜乱码| 国产成人精品一,二区| 中文字幕av成人在线电影| av卡一久久| 免费黄色在线免费观看| 18+在线观看网站| 日本wwww免费看| 久久久久国产精品人妻一区二区| 国产精品一区二区三区四区免费观看| 九九久久精品国产亚洲av麻豆| 乱系列少妇在线播放| 日本免费在线观看一区| 秋霞在线观看毛片| 亚洲精品,欧美精品| 韩国高清视频一区二区三区| 亚洲天堂av无毛| 亚洲精品中文字幕在线视频 | 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 女性生殖器流出的白浆| 午夜福利视频精品| 人人妻人人爽人人添夜夜欢视频 | 偷拍熟女少妇极品色| 亚洲性久久影院| 亚洲精品久久午夜乱码| 亚洲av二区三区四区| 亚洲婷婷狠狠爱综合网| 免费观看在线日韩| 中文乱码字字幕精品一区二区三区| 制服丝袜香蕉在线| 在线免费十八禁| 欧美精品人与动牲交sv欧美| 免费在线观看成人毛片| av免费观看日本| 欧美成人精品欧美一级黄| 久久韩国三级中文字幕| 内地一区二区视频在线| 最黄视频免费看| 18禁动态无遮挡网站| 日韩不卡一区二区三区视频在线| 中文在线观看免费www的网站| 爱豆传媒免费全集在线观看| av专区在线播放| 伦理电影免费视频| 一个人看视频在线观看www免费| 一区二区三区四区激情视频| 亚洲最大成人中文| 久久ye,这里只有精品| 成年女人在线观看亚洲视频| 青春草亚洲视频在线观看| 大又大粗又爽又黄少妇毛片口| 在线观看一区二区三区| 久久久久国产精品人妻一区二区| 国产高清不卡午夜福利| 91久久精品国产一区二区三区| 在线观看人妻少妇| 日韩三级伦理在线观看| 80岁老熟妇乱子伦牲交| 黄色怎么调成土黄色| 色综合色国产| 国产精品久久久久成人av| 大话2 男鬼变身卡| 免费在线观看成人毛片| 日韩中字成人| 国产乱人偷精品视频| 亚洲自偷自拍三级| 国产精品人妻久久久久久| 少妇裸体淫交视频免费看高清| 免费大片黄手机在线观看| 国产免费一区二区三区四区乱码| 一级毛片久久久久久久久女| 欧美精品一区二区大全| 18禁在线无遮挡免费观看视频| 亚洲无线观看免费| 国产成人精品婷婷| av在线app专区| 亚洲成人中文字幕在线播放| 亚洲精品成人av观看孕妇| 人妻夜夜爽99麻豆av| 欧美成人一区二区免费高清观看| 人妻一区二区av| 男女国产视频网站| 黄色欧美视频在线观看| 亚洲国产av新网站| 免费看av在线观看网站| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 成人亚洲精品一区在线观看 | 精品一区二区三区视频在线| 国产精品一区二区性色av| 久久国内精品自在自线图片| 亚洲经典国产精华液单| 久久国内精品自在自线图片| av网站免费在线观看视频| 亚洲欧美一区二区三区国产| 天天躁夜夜躁狠狠久久av| 久久精品久久久久久久性| 亚洲欧美成人综合另类久久久| 久久精品国产鲁丝片午夜精品| 亚洲精品乱久久久久久| 人人妻人人看人人澡| 久久精品久久精品一区二区三区| 国产av码专区亚洲av| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区| 内地一区二区视频在线| 2018国产大陆天天弄谢| 欧美+日韩+精品| 久久6这里有精品| 国产有黄有色有爽视频| 人妻少妇偷人精品九色| 美女主播在线视频| 国产成人a区在线观看| 97在线视频观看| 国产又色又爽无遮挡免| 少妇裸体淫交视频免费看高清| 亚洲色图av天堂| 日韩亚洲欧美综合| 欧美精品一区二区大全| 香蕉精品网在线| 中国三级夫妇交换| 欧美老熟妇乱子伦牲交| 国产免费一区二区三区四区乱码| 在线观看av片永久免费下载| 纯流量卡能插随身wifi吗| 亚洲精品乱久久久久久| 国产伦理片在线播放av一区| 能在线免费看毛片的网站| 日韩av不卡免费在线播放| 久久久久久久大尺度免费视频| 18禁裸乳无遮挡动漫免费视频| 国产成人午夜福利电影在线观看| 久久 成人 亚洲| 国产成人精品久久久久久| 青春草国产在线视频| 99久久精品热视频| 少妇的逼水好多| 少妇猛男粗大的猛烈进出视频| 国产免费一级a男人的天堂| 超碰97精品在线观看| 成人国产av品久久久| 毛片一级片免费看久久久久| 纵有疾风起免费观看全集完整版| 免费看日本二区| 亚洲精品乱久久久久久| 国产伦在线观看视频一区| 精品久久久精品久久久| 欧美亚洲 丝袜 人妻 在线| 国产深夜福利视频在线观看| 欧美另类一区| 国产毛片在线视频| 深爱激情五月婷婷| 91精品一卡2卡3卡4卡| 舔av片在线| 干丝袜人妻中文字幕| 亚洲国产色片| 午夜日本视频在线| 青春草国产在线视频| 黄片wwwwww| 亚洲美女搞黄在线观看| 欧美bdsm另类| 欧美成人午夜免费资源| 久久久成人免费电影| 亚洲人成网站高清观看| 少妇的逼好多水| 欧美老熟妇乱子伦牲交| 视频中文字幕在线观看| 久久婷婷青草| 亚洲国产欧美人成| av在线老鸭窝| 在线观看一区二区三区| 久久精品国产鲁丝片午夜精品| 18禁裸乳无遮挡动漫免费视频| 美女中出高潮动态图| 成人18禁高潮啪啪吃奶动态图 | av播播在线观看一区| 麻豆乱淫一区二区| 精华霜和精华液先用哪个| 亚洲av不卡在线观看| 亚洲伊人久久精品综合| 精品人妻一区二区三区麻豆| 日韩不卡一区二区三区视频在线| 亚洲美女搞黄在线观看| 26uuu在线亚洲综合色| 直男gayav资源| 99热这里只有是精品50| 日韩成人伦理影院| 一级毛片久久久久久久久女| 日韩国内少妇激情av| 久久6这里有精品| 国产欧美亚洲国产| 国产精品久久久久久av不卡| 国产av精品麻豆| 日本欧美国产在线视频| 午夜福利高清视频| 在线免费观看不下载黄p国产| 精品人妻偷拍中文字幕| a级毛色黄片| 成人毛片a级毛片在线播放| 欧美精品亚洲一区二区| 日本免费在线观看一区| 国产精品熟女久久久久浪| 在线看a的网站| 精品人妻熟女av久视频| freevideosex欧美| 亚洲成人一二三区av| 国产一区二区在线观看日韩| av在线蜜桃| 熟女人妻精品中文字幕| 欧美亚洲 丝袜 人妻 在线| 91精品国产九色| 99久久精品热视频| 国产成人a区在线观看| 国产精品欧美亚洲77777| 亚洲无线观看免费| 国产av国产精品国产| 精品国产一区二区三区久久久樱花 | 国产精品免费大片| 欧美日本视频| 久久久久久久精品精品| 黄色配什么色好看| 国产免费一区二区三区四区乱码| 欧美xxⅹ黑人| 午夜激情福利司机影院| 亚洲成人中文字幕在线播放| 亚洲国产精品国产精品| 国产男人的电影天堂91| 18+在线观看网站| 啦啦啦中文免费视频观看日本| 免费av不卡在线播放| 色婷婷久久久亚洲欧美| 91aial.com中文字幕在线观看| 美女中出高潮动态图| 18禁在线无遮挡免费观看视频| av播播在线观看一区| 久久精品熟女亚洲av麻豆精品| 好男人视频免费观看在线| 91aial.com中文字幕在线观看| 成年女人在线观看亚洲视频| av国产久精品久网站免费入址| av在线观看视频网站免费| 免费少妇av软件| 国产成人a∨麻豆精品| 日本与韩国留学比较| 亚洲丝袜综合中文字幕| 性色avwww在线观看| 黄片wwwwww|