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

    基于Cyt b基因的雅魯藏布江下游墨脫江段及察隅河墨脫裂腹魚的遺傳多樣性及種群歷史動(dòng)態(tài)分析

    2019-10-18 08:42:12林鵬程熊淑蓉唐芬蓮劉煥章
    水生生物學(xué)報(bào) 2019年5期
    關(guān)鍵詞:腹魚裂腹雅魯藏布江

    俞 丹 張 智, 張 健, 林鵬程 熊淑蓉 唐芬蓮 劉煥章

    (1. 中國科學(xué)院水生生物研究所水生生物多樣性與保護(hù)重點(diǎn)實(shí)驗(yàn)室, 武漢 430072;2. 中國科學(xué)院大學(xué), 北京 100049; 3. 西華師范大學(xué), 南充 637002)

    生物多樣性是地球上寶貴的生物資源[1], 包括遺傳多樣性(Genetic diversity)、功能多樣性(Functional diversity)和系統(tǒng)發(fā)育多樣性(Phylogenetic diversity)等不同的維度[2]。遺傳多樣性又稱為基因多樣性(Genes diversity), 廣義上是指地球上所有生物攜帶的遺傳信息的總和; 狹義上是指種內(nèi)不同種群或同一種群內(nèi)不同個(gè)體的遺傳變異的總和。通常所講的遺傳多樣性是指狹義的遺傳多樣性。生物的遺傳多樣性水平是生物長期進(jìn)化的產(chǎn)物, 是評估生物資源現(xiàn)狀的一個(gè)重要參數(shù), 遺傳多樣性越高,即遺傳變異越豐富, 表示對環(huán)境的適應(yīng)能力就越強(qiáng)。反之, 在一些瀕危物種或種群中, 由于生境惡化、過度捕撈以及生物入侵等影響, 會(huì)導(dǎo)致遺傳多樣性降低, 出現(xiàn)遺傳衰退, 最終可能滅絕[3]。因此,對遺傳多樣性的研究不僅可以揭示物種或種群的起源時(shí)間、方式等進(jìn)化歷史, 還可以進(jìn)一步分析其進(jìn)化潛力和發(fā)展趨勢。

    雅魯藏布江是青藏高原上的第一大河, 也是世界上海拔最高的大河之一, 它發(fā)源于喜馬拉雅北麓的杰馬央宗冰川, 由西向東橫貫西藏南部, 繞過喜馬拉雅山脈最東端的南迦巴瓦峰轉(zhuǎn)向南流, 從巴昔卡進(jìn)入印度境內(nèi)后稱布拉馬普特拉河[4]。雅魯藏布江下游從派鎮(zhèn)到巴昔卡, 呈高山峽谷, 同時(shí)匯集了金珠藏布、西莫河、西貢河等眾多支流, 孕育了豐富的魚類資源, 據(jù)調(diào)查雅魯藏布江下游分布有魚類近 20 種[5]。

    裂腹魚類是適應(yīng)高海拔冷水環(huán)境的一個(gè)重要類群, 是硬骨魚綱(Osteichthyes)鯉形目(Cyprini-formes)鯉科(Cyprinidae)裂腹魚亞科(Schizothoracinae)魚類的統(tǒng)稱, 其共同的外部態(tài)特征是具有臀鱗, 即在肛門和臀鰭基部的兩側(cè)各自排列著一行特化的大型鱗片, 形成了腹中線上的一條裂縫, 故稱“裂腹”。曹文宣等[6]依據(jù)裂腹魚類中體鱗退化程度、下咽齒行數(shù)以及觸須數(shù)目等, 將裂腹魚類分成適應(yīng)高原不同環(huán)境的 3 個(gè)等級類群, 即原始類群、特化類群和高度特化類群。每一個(gè)等級分別代表了青藏高原在隆起過程中的特定的歷史階段。墨脫裂腹魚(Schizothorax molesworthi)屬于裂腹魚3 種等級類群中的原始類群[7], 分布范圍較廣, 在我國境內(nèi)主要分布于雅魯藏布江下游及其支流和伊洛瓦底江上游支流, 國外也分布于印度、尼泊爾及緬甸。墨脫裂腹魚主要棲息于湍急河流之中, 主要以硅藻為食[5]。目前對墨脫裂腹魚的研究甚少, 僅見于少量的魚類資源調(diào)查, 尚未見種群遺傳多樣性、種群結(jié)構(gòu)方面的研究。

    線粒體細(xì)胞色素b(Cytb)基因不僅能有效評估物種和種群的遺傳多樣性, 同時(shí)也能敏銳地檢測出潛在的種群遺傳分化問題[8,9]。本研究對采自雅魯藏布江下游墨脫江段和察隅河的墨脫裂腹魚種群進(jìn)行Cytb基因測序, 分析墨脫裂腹魚不同地理種群的遺傳多樣性、種群分化及種群歷史動(dòng)態(tài)等, 為雅魯藏布江魚類資源的保護(hù)提供科學(xué)的理論依據(jù)。

    1 材料與方法

    1.1 研究材料

    本研究所用樣本于2015—2016年采自雅魯藏布江下游墨脫江段和察隅河, 共計(jì)167尾。對雅魯藏布江下游墨脫江段的干流江段、5條重要支流(白馬西路河、哈果河、金珠藏布、西貢河、西莫河)及察隅河進(jìn)行采樣, 根據(jù)其地理分布, 最終將其劃分為 7個(gè)地理種群。所有干流采樣點(diǎn)視為干流種群, 其他每個(gè)支流采樣點(diǎn)視為一個(gè)地理種群, 即用支流名稱對種群進(jìn)行命名, 察隅河為單獨(dú)一個(gè)地理種群。采集標(biāo)本均用 95% 酒精固定, 保存于中國科學(xué)院水生生物研究所水生生物博物館。樣本詳細(xì)信息及分布圖見表 1和圖 1。

    1.2 基因組 DNA 的提取、PCR 擴(kuò)增及測序

    基因組DNA的提取采用高鹽抽提法[10]。Cytb基因擴(kuò)增和測序引物為L14724 (5′-GACTTGAA AAACCACCGTTG-3′)和H15915 (5′-CTCCGATCT CCGGATTACAAGAC-3′)[11]。PCR反應(yīng)總體積為20 μL, 反應(yīng)體系包括: 10×Buffer 2 μL, dNTPs 0.25 μL(10 mmol/L), 引物各0.5 μL (10 μmol/L),TaqDNA聚合酶 0.5 U, 基因組DNA模板1—2 μL, 加滅菌雙蒸水至 20 μL。PCR反應(yīng)條件為: 94℃預(yù)變性 3min;然后循環(huán)包括: 94℃ 變性 45s, 52—58℃退火 45s,72℃ 延伸1min, 共 35 個(gè)循環(huán); 最后在 72℃ 延伸8min。PCR產(chǎn)物用 0.8% 瓊脂糖凝膠檢測目的片段, 成功后送上海生工生物工程股份有限公司完成純化與測序。

    1.3 數(shù)據(jù)分析

    序列比對采用CLUSTALX 2.1軟件[12]。采用MEGA 7.0 將堿基序列翻譯成相應(yīng)氨基酸序列進(jìn)行序列校正。采用PAUP* 4.0b10[13]對序列進(jìn)行初步分析。在MEGA中采用Kimura雙參數(shù)模型進(jìn)行遺傳距離計(jì)算[14]。堿基替代模型選擇采用Jmodeltest 2.0[15]。以弧唇裂腹魚(Schizothorax curilabiatus)為外類群, 采用鄰接法(Neighbor-joining, NJ)、貝葉斯法(Bayesian inferences, BI) 和最大似然法(Maximum likelihood, ML) 構(gòu)建分子系統(tǒng)發(fā)育樹。NJ分析在PAUP*軟件中進(jìn)行, BI分析在MrBayes 3.2 軟件[16]中進(jìn)行, ML分析在PhyML 3.0 軟件[17]中進(jìn)行。采用DnaSP v5.10 軟件統(tǒng)計(jì)種群遺傳多樣性參數(shù)[18]。利用Network5.0 軟件以中接法(Median-joining)構(gòu)建單倍型網(wǎng)絡(luò)關(guān)系圖, 分析單倍型之間的進(jìn)化關(guān)系[19]。使用Arlequin3.5軟件中的分子方差分析(Analysis of Molecular Variance, AMOVA)[20]估算種群遺傳結(jié)構(gòu)和地理群體遺傳變異的分布, 計(jì)算種群間遺傳分化指數(shù)FST。種群擴(kuò)張檢驗(yàn)包括: 在Arlequin3.5 計(jì)算錯(cuò)配分布值, Tajima’sD和Fu’sFs值, 通過公式τ=2ut估算種群擴(kuò)張發(fā)生時(shí)間 [t為種群擴(kuò)張以來的時(shí)間,u為序列的變異速率;u= 2μk,μ為基因序列的突變率,k分析序列的堿基數(shù)(bp), 擴(kuò)張時(shí)間T=t×世代時(shí)間], 進(jìn)化速率采用 1.69% /百萬年[21]。使用BEAST v1.8.2軟件[22]進(jìn)行種群歷史動(dòng)態(tài)的Bayesian skyline plot (BSP)分析, 最后在Tracer v 1.4軟件中構(gòu)圖[23]。

    表 1 雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚的樣本采集信息及遺傳多樣性信息Tab. 1 Genetic diversity for the Schizothorax molesworthi from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River

    2 結(jié)果

    2.1 Cyt b基因序列變異分析

    本研究共獲得167尾墨脫裂腹魚的線粒體Cytb基因序列。比對后得到序列長1140 bp, 變異位點(diǎn)18個(gè), 其中包含12個(gè)簡約信息位點(diǎn)。所有Cytb基因序列的轉(zhuǎn)換和顛換均未達(dá)飽和, 轉(zhuǎn)換數(shù)明顯大于顛換數(shù), 其平均Ti/Tv的值為2.94。所有個(gè)體平均堿基組成: A=27.30%、T=29.83%、C=26.68%、G=16.19%。A+T的含量(57.13%)高于G+C的含量(42.97%)。堿基組成存在強(qiáng)烈的偏倚, G的含量尤其在密碼子第三位的含量最低(9.90%)。

    2.2 單倍型多樣性和核苷酸多樣性

    在167個(gè)墨脫裂腹魚樣本中, 共檢測到21個(gè)單倍型, 單倍型多樣性為0.768 (0.030), 核苷酸多樣性為0.00167 (0.00011), 不同地理種群的單倍型多樣性和核苷酸多樣性見表 1。其中來自干流(GL)種群的單倍型多樣性和核苷酸多樣性最高, 分別為0.830(0.048)和0.00196 (0.00021); 西莫河(XM)種群的單倍型多樣性和核苷酸多樣性最低, 分別為0.462 (0.136)和0.00054 (0.00020)。從單倍型分布情況來看, 除了西莫河(XM)種群, 其他種群均含有特有單倍型,特有單倍型共有14個(gè), 其中察隅河種群含有4個(gè)特有單倍型。另外, 不同種群之間具有共享單倍型,例如除了金珠藏布(JZZB)種群, 其他所有種群都共享單倍型Hap 4。

    2.3 系統(tǒng)發(fā)育關(guān)系

    利用NJ、BI和ML構(gòu)建墨脫裂腹魚21個(gè)單倍型之間的系統(tǒng)發(fā)育樹。結(jié)果顯示, 3種方法分析得到的系統(tǒng)樹拓?fù)浣Y(jié)構(gòu)基本一致, 所有的單倍型可以聚為一個(gè)單系, 但在系統(tǒng)樹上不按照地理分布聚類,各結(jié)點(diǎn)支持率較低, 存在多歧分支(圖 2)。單倍型進(jìn)化網(wǎng)絡(luò)關(guān)系圖顯示, 單倍型的網(wǎng)絡(luò)進(jìn)化關(guān)系呈星狀分布, 除了金珠藏布(JZZB), 其他所有種群共享單倍型Hap4, 屬于原始單倍型和進(jìn)化中心, 部分單倍型之間通過缺失的中間單倍型相互連接, 不存在顯著的地理分布格局(圖 3)。

    2.4 種群間遺傳分化與種群結(jié)構(gòu)

    基于Cytb基因序列兩兩群體之間的遺傳分化指數(shù)(FST)范圍為-0.014—0.771 (表 2)。由表 2可知, 金珠藏布(JZZB)種群與其他種群之間的FST最大(0.372—0.771,P<0.01), 干流和察隅種群與其他種群之間的FST次之(0.093—0.221,P<0.01), 剩余種群之間的FST較小。在不分組的情況下對所有種群進(jìn)行AMOVA分析, 發(fā)現(xiàn)遺傳變異主要來自于種群內(nèi)的遺傳差異(種群內(nèi)的變異百分比占76.12%,P<0.01), 種群間的遺傳變異相對較小(23.88%)。將所有墨脫裂腹魚樣品按地理分布劃分為2—6個(gè)組,具有最大組間變異(FCT值)的劃分方案被假定為最符合自然的群體結(jié)構(gòu)模式。在所有分組方案中, 當(dāng)所有群體劃分為2個(gè)組時(shí), 即金珠藏布(JZZB)種群為一組, 剩余6個(gè)種群為一組, 得到的組間遺傳差異FCT最大, 組間變異占50.63% (P<0.01), 組內(nèi)種群間變異僅占 6.13%, 種群內(nèi)變異占43.25% (表 3)。這說明雅魯藏布江下游墨脫江段和察隅河的墨脫裂腹魚中存在較明顯的種群遺傳結(jié)構(gòu), 金珠藏布(JZZB)種群與其他種群具有顯著分化。

    圖 1 墨脫裂腹魚樣本采集地點(diǎn)Fig. 1 Sample locations for Schizothorax molesworthi

    2.5 種群歷史動(dòng)態(tài)分析

    核苷酸錯(cuò)配分析發(fā)現(xiàn), 所有群體均呈現(xiàn)單峰分布, 表明這些群體經(jīng)歷了近期擴(kuò)張; 而中性檢驗(yàn)顯示, 除了西莫河(XM)的Fu’sFs和Tajima’sD值呈現(xiàn)顯著的負(fù)值, 其他群體均符合中性進(jìn)化的假設(shè), 未檢測到種群擴(kuò)張,SSD檢驗(yàn)結(jié)果顯示4個(gè)種群(CY、BMXL、HG和JZZB)顯著偏離了種群擴(kuò)張模型, 但是Hri檢驗(yàn)結(jié)果顯示所有墨脫裂腹魚種群未顯著偏離種群擴(kuò)張模型, 根據(jù)1.69%/百萬年的突變速率,推算出墨脫裂腹魚不同地理種群發(fā)生擴(kuò)張的大致時(shí)間為0.027—0.153 Ma (表 4)。BSP分析了有效種群大小隨溯祖時(shí)間的動(dòng)態(tài)變化曲線, 所有個(gè)體組成的種群在0—0.01 Ma期間存在種群擴(kuò)張, 其他地理種群未檢測到種群有明顯的擴(kuò)張與收縮(圖 4)。

    3 討論

    3.1 雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚的遺傳多樣性和種群分化

    遺傳多樣性是生物多樣性的基礎(chǔ), 也是物種長期生存和進(jìn)化的前提。單倍型多樣性(h)和核苷酸多樣性(π)是衡量一個(gè)物種群體DNA變異程度的重要指標(biāo)[24]。本研究在167個(gè)墨脫裂腹魚樣本中共檢測到21個(gè)單倍型, 其整體單倍型多樣性和核苷酸多樣性分別為0.768和0.00167, 低于一些裂腹魚類的遺傳多樣性, 例如怒江的怒江裂腹魚(S.nukiangensis)的h和π分別為 0.965和 0.0024[25], 雅魯藏布江中上游異齒裂腹魚(S. oconnori)的h和π分別為 0.833—0.985 和 0.00311—0.01491[26], 但高于瑪柯河齊口裂腹魚(S. prenanti)的h和π (0.5606和0.0007)[27]。因此, 墨脫裂腹魚呈現(xiàn)單倍型多樣性高, 核苷酸多樣性低的模式, 這可能是群體受到瓶頸效應(yīng)后又發(fā)生了快速擴(kuò)張, 在擴(kuò)張過程中積累的遺傳變異會(huì)增加群體的單倍型多樣性, 但是核苷酸多樣性依然較低[28,29]; 亦可能是墨脫裂腹魚的種群進(jìn)化歷史較短的原因[30]。

    遺傳分化指數(shù)(F-statistic,FST)可在一定程度上指示種群間基因流和遺傳漂變的程度, 是反映群體進(jìn)化歷史的重要參數(shù)。Wright等[31]認(rèn)為, 若種群遺傳分化指數(shù)FST值在0—0.05, 說明不存在分化;FST值在0.05—0.15, 為中度分化;FST值在0.15—0.25, 則為高度分化。本研究中FST和AMOVA分析結(jié)果顯示不同地理種群間的遺傳分化指數(shù)在-0.014—0.771。其中, 哈果河(HG)和西貢河(XG)、白馬西路河(BMXL)和西莫河(XM)種群之間不存在種群分化; 察隅河(CY)和西貢河(XG)、干流(GL)和其他種群等存在具有中度分化(FST:0.099—0.221,P<0.01); 而金珠藏布(JZZB)與其他種群之間呈現(xiàn)高度分化(FST: 0.372—0.771,P<0.01),推測金珠藏布河口段的一個(gè)跌水造成了該種群與干流及其他種群的隔離, 長時(shí)間的隔離和生態(tài)適應(yīng)是造成金珠藏布(JZZB)種群與其他種群分化的主要原因。此外, 單倍型分析發(fā)現(xiàn)除了金珠藏布(JZZB)種群, 其他所有種群都共享單倍型Hap 4, 出現(xiàn)這個(gè)結(jié)果的原因也可能是因?yàn)樵摲N群取樣數(shù)量較少(11個(gè)個(gè)體), 后期需要增加樣本量確認(rèn)。通常來說, 裂腹魚類的種群分化程度相對較低[32,33], 但在怒江的怒江裂腹魚中同樣檢測到了中等程度的分化[25], 與本研究結(jié)果相似。另外, 單倍型分析發(fā)現(xiàn)66.67%的單倍型為特有單倍型, 這也暗示了不同地理種群之間的基因交流較少, 分化較高。本研究結(jié)果中值得一提的是, 雖然察隅河與雅魯藏布江下游地理距離較遠(yuǎn), 但是察隅河與墨脫江段其他種群之間(除了JZZB)的FST為0.093—0.169 (P<0.01), 屬于中等分化水平, 表明察隅河種群與雅魯藏布江下游種群尚有少量的基因交流。

    圖 2 基于線粒體Cyt b基因單倍型構(gòu)建的雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚的鄰接法、貝葉斯法和最大似然法樹Fig. 2 Neighbor-joining (NJ)/Bayesian inferences (BI)/maximum likelihood (ML) tree of Schizothorax molesworthi based on Cyt b gene from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River節(jié)點(diǎn)處的數(shù)值代表NJ/BI/ML分析獲得的支持率, “-”表示支持率低于10%, 表示枝長代表分歧度?!癏ap”表示單倍型, 其后的序號為單倍型的編號Values at the node represent support values in NJ/BI/ML analysis respectively. “-” represent support values are less than 10%. The length of branch indicates the percentage of divergence. “Hap”denotes Haplotype, subsequent number means haplotype code

    3.2 雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚的種群歷史動(dòng)態(tài)

    對于本研究中發(fā)現(xiàn)的墨脫裂腹魚部分地理種群核苷酸多樣性較低的現(xiàn)象, 我們認(rèn)為地質(zhì)、冰期事件影響下的種群歷史動(dòng)態(tài)變化(如近期的種群擴(kuò)張)可能是導(dǎo)致部分地理種群遺傳多樣性偏低的重要原因, 尤其是分布在青藏高原的裂腹魚類, 屬于冷水魚類, 冰期與間冰期的旋回會(huì)對物種的遺傳多樣性和種群分化產(chǎn)生重要影響。在種群動(dòng)態(tài)分析中, Tajima’sD和Fu’sFs值是兩種比較常用的通過中性檢驗(yàn)研究種群在歷史上是否發(fā)生過擴(kuò)張的參數(shù)。Tajima’sD值和Fu’sFs值接近零時(shí)表明種群處于穩(wěn)定狀態(tài), 負(fù)的D值和Fs值以及差異顯著的P值被認(rèn)為種群在歷史上有擴(kuò)張跡象[34]。但相對而言,Fu’sFs檢驗(yàn)對種群擴(kuò)張更加敏感[35]。在本研究中,所有個(gè)體構(gòu)成的種群Tajima’sD值為不顯著負(fù)值(P>0.05), 而Fu’sFs值則為差異極顯著負(fù)值(P<0.01),暗示雅魯藏布江下游墨脫江段和察隅河的墨脫裂腹魚可能在近期發(fā)生過種群擴(kuò)張。西莫河(XM)的Tajima’sD和Fu’sFs值均呈現(xiàn)顯著的負(fù)值, 表明西莫河種群經(jīng)歷了種群擴(kuò)張。此外, 偏差平方和(SSD)及粗糙指數(shù)(Hri)的統(tǒng)計(jì)檢驗(yàn)不顯著時(shí)(P>0.05), 表明不能拒絕群體擴(kuò)張的假說, 即符合群體擴(kuò)張假說[36]。SSD檢驗(yàn)結(jié)果顯示4個(gè)種群(CY、BMXL、HG和JZZB)顯著偏離了種群擴(kuò)張模型, 但是Hri檢驗(yàn)結(jié)果顯示所有墨脫裂腹魚群體未顯著偏離種群擴(kuò)張模型, 并估算獲得不同地理種群的擴(kuò)張時(shí)間約在0.004—0.153 Ma左右, 雅魯藏布江下游墨脫江段和察隅河的所有墨脫裂腹魚構(gòu)成的群體擴(kuò)張時(shí)間為0.137 Ma, 剛好與末次間冰期發(fā)生時(shí)間一致(MIS5: 78—140 Ka)[37,38]。最后, 墨脫裂腹魚種群動(dòng)態(tài)隨時(shí)間變化的BSP圖顯示, 所有墨脫裂腹魚組成的群體在0—0.01 Ma期間存在種群擴(kuò)張, 與最后一次冰期結(jié)束時(shí)間一致(0.075—0.01 Ma)[37,38], 但BSP分析在其他地理種群未檢測到種群有明顯的擴(kuò)張與收縮。上述結(jié)果進(jìn)一步驗(yàn)證了末次冰期事件對雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚種群歷史動(dòng)態(tài)的影響。

    圖 3 基于線粒體Cyt b基因的雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚的單倍型網(wǎng)絡(luò)圖Fig. 3 The haplotype network of Schizothorax molesworthi based on Cyt b gene from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River圓圈的大小表示單倍型包含的個(gè)體數(shù)多少; 中間白色圓圈代表缺失的單倍型The circle size represents the number of individual haplotypes. The white circle represents lacked haplotype

    表 2 雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚不同地理種群的遺傳分化系數(shù)(FST)Tab. 2 Pairwise FST values among different geographic populations of Schizothorax molesworthi from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River

    表 3 雅羅藏布江下游墨脫江段和察隅河墨脫裂腹魚的分子方差分析Tab. 3 Analysis of molecular variance analysis among populations of Schizothorax molesworthi from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River

    表 4 雅羅藏布江下游墨脫江段和察隅河墨脫裂腹魚的中性檢驗(yàn)和錯(cuò)配分析Tab. 4 Neutrality tests and mismatch distribution values for all populations of Schizothorax molesworthi from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River

    3.3 雅魯藏布江墨脫裂腹魚的保護(hù)

    墨脫裂腹魚作為雅魯藏布江下游支流水生生態(tài)系統(tǒng)的重要組成部分, 對維系整個(gè)生態(tài)系統(tǒng)結(jié)構(gòu)完整和功能的正常發(fā)揮具有重要作用。本研究結(jié)果顯示墨脫裂腹魚整體遺傳多樣性水平較低, 暗示墨脫裂腹魚容易受到環(huán)境變化等因素的影響。此外, 種群動(dòng)態(tài)歷史分析也發(fā)現(xiàn)墨脫裂腹魚在冰期受到了氣候變化的影響。因此, 在全球氣候變化的背景下, 墨脫裂腹魚的物種保護(hù)工作尤為重要。在本研究中除了西莫河(XM)種群外, 其他種群均含有特有單倍型, 尤其是察隅河種群含有4個(gè)特有單倍型。因此, 為了防止遺傳資源的喪失, 在遺傳資源保護(hù)管理中應(yīng)對察隅河資源加強(qiáng)保護(hù)。此外, 金珠藏布種群與其他6個(gè)地理種群出現(xiàn)了較高水平分化,這與金珠藏布河流沿線正好是波密到墨脫的公路有關(guān), 地質(zhì)條件復(fù)雜, 災(zāi)害頻繁, 河道生境也履遭破壞, 建議加強(qiáng)金珠藏布的生境修復(fù)工作, 提高河流的生境適宜性, 保護(hù)種群資源。

    圖 4 雅魯藏布江下游墨脫江段和察隅河墨脫裂腹魚種群動(dòng)態(tài)隨時(shí)間變化的BSP圖Fig. 4 Bayesian skyline plots (BSP) for population dynamics with time of Schizothorax molesworthi from the Motuo area of lower reaches of the Yarlung Zangbo River and Lohit River黑色線表示中值, 陰影部分為95%的置信區(qū)間(HPD)The black line indicates the median; the shadow part shows the 95% highest probability density (HPD)

    致謝:

    感謝中國科學(xué)院水生生物研究所劉飛、黎明政、鞏政、王健在采樣方面給予的幫助。

    猜你喜歡
    腹魚裂腹雅魯藏布江
    人工養(yǎng)殖2 種裂腹魚生物學(xué)指數(shù)和肌肉營養(yǎng)價(jià)值的比較分析
    齊口裂腹魚集群行為對流態(tài)的響應(yīng)
    南疆地區(qū)四種土著裂腹魚類幼魚對氯化鈉和碳酸氫鈉的急性耐受能力
    仿生態(tài)養(yǎng)殖短須裂腹魚母本繁殖特性研究
    淺析畢節(jié)市裂腹魚繁養(yǎng)現(xiàn)狀及發(fā)展對策探究
    異齒裂腹魚產(chǎn)后護(hù)理技術(shù)
    中國與孟加拉國在雅魯藏布江河流治理中的合作與問題探究
    商(2016年32期)2016-11-24 15:28:37
    西藏雅魯藏布江區(qū)橋梁水文計(jì)算
    西藏科技(2015年2期)2015-09-26 12:18:32
    雅魯藏布江特大橋的臨時(shí)鋼棧橋結(jié)構(gòu)計(jì)算
    西藏科技(2015年3期)2015-09-26 12:11:11
    大美尼洋河
    色5月婷婷丁香| av专区在线播放| 国产淫语在线视频| 少妇高潮的动态图| 蜜臀久久99精品久久宅男| 美女国产视频在线观看| 欧美性感艳星| 亚洲av免费高清在线观看| 国产av国产精品国产| 老司机影院成人| 亚洲av不卡在线观看| 欧美一级a爱片免费观看看| kizo精华| 欧美日韩综合久久久久久| 国产精品无大码| 久久久精品94久久精品| 韩国av在线不卡| 天堂网av新在线| 国产老妇女一区| 男人舔女人下体高潮全视频| 美女内射精品一级片tv| 丝瓜视频免费看黄片| 十八禁网站网址无遮挡 | 亚洲美女视频黄频| 亚洲aⅴ乱码一区二区在线播放| 美女高潮的动态| 18+在线观看网站| 久久热精品热| 国产视频首页在线观看| 韩国高清视频一区二区三区| 精品少妇黑人巨大在线播放| 少妇被粗大猛烈的视频| 毛片一级片免费看久久久久| 亚洲美女搞黄在线观看| 男女那种视频在线观看| 白带黄色成豆腐渣| 久久久久久久久大av| 性插视频无遮挡在线免费观看| 女的被弄到高潮叫床怎么办| 联通29元200g的流量卡| 日韩成人伦理影院| 在线a可以看的网站| 久久99蜜桃精品久久| 最新中文字幕久久久久| av免费在线看不卡| 日日撸夜夜添| 夫妻午夜视频| 啦啦啦中文免费视频观看日本| 好男人视频免费观看在线| 一本一本综合久久| 免费在线观看成人毛片| 午夜精品国产一区二区电影 | 亚洲欧美中文字幕日韩二区| 女人被狂操c到高潮| av线在线观看网站| 一级毛片aaaaaa免费看小| 尾随美女入室| 欧美激情久久久久久爽电影| 成人毛片60女人毛片免费| 丰满人妻一区二区三区视频av| 亚洲精品影视一区二区三区av| 午夜精品国产一区二区电影 | 寂寞人妻少妇视频99o| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 麻豆成人午夜福利视频| 久久久久久久久中文| 久久热精品热| 亚洲精品久久午夜乱码| 亚洲国产成人一精品久久久| 天天一区二区日本电影三级| 欧美日韩精品成人综合77777| 日本av手机在线免费观看| 国产单亲对白刺激| 国产毛片a区久久久久| 亚洲欧美日韩无卡精品| 欧美成人一区二区免费高清观看| 国产久久久一区二区三区| 看十八女毛片水多多多| 亚洲一级一片aⅴ在线观看| 性色avwww在线观看| 亚洲图色成人| 成人毛片60女人毛片免费| 有码 亚洲区| 熟女人妻精品中文字幕| 欧美变态另类bdsm刘玥| 97人妻精品一区二区三区麻豆| 国产黄片视频在线免费观看| 日日啪夜夜爽| 久久久亚洲精品成人影院| 欧美另类一区| 18禁动态无遮挡网站| 成年版毛片免费区| 一级毛片 在线播放| 2021少妇久久久久久久久久久| 中文字幕久久专区| 黄片无遮挡物在线观看| 国产免费视频播放在线视频 | 精品久久国产蜜桃| 最近2019中文字幕mv第一页| 日韩国内少妇激情av| 少妇丰满av| 性插视频无遮挡在线免费观看| 中文在线观看免费www的网站| 国产熟女欧美一区二区| 精品久久国产蜜桃| 又爽又黄无遮挡网站| 日本熟妇午夜| 爱豆传媒免费全集在线观看| 麻豆乱淫一区二区| 欧美最新免费一区二区三区| 日韩av免费高清视频| 97精品久久久久久久久久精品| 亚洲精华国产精华液的使用体验| 街头女战士在线观看网站| 日本爱情动作片www.在线观看| 午夜精品国产一区二区电影 | 国产乱人视频| 69人妻影院| 亚洲国产成人一精品久久久| 两个人的视频大全免费| 欧美性感艳星| 国产亚洲av片在线观看秒播厂 | 欧美xxxx黑人xx丫x性爽| 国产亚洲av片在线观看秒播厂 | 精品久久久精品久久久| 可以在线观看毛片的网站| 亚洲精品乱码久久久久久按摩| 国产精品爽爽va在线观看网站| 免费看日本二区| 成年版毛片免费区| 2018国产大陆天天弄谢| 婷婷色av中文字幕| 一级毛片黄色毛片免费观看视频| 综合色av麻豆| 久久精品国产鲁丝片午夜精品| 你懂的网址亚洲精品在线观看| 国产国拍精品亚洲av在线观看| 午夜激情福利司机影院| 国产黄片美女视频| 2018国产大陆天天弄谢| 极品教师在线视频| 黄色日韩在线| 高清av免费在线| 国产中年淑女户外野战色| 内射极品少妇av片p| 搡女人真爽免费视频火全软件| 一区二区三区免费毛片| av专区在线播放| 只有这里有精品99| 午夜老司机福利剧场| 国产淫语在线视频| 人人妻人人澡欧美一区二区| 亚洲国产精品成人综合色| 中文在线观看免费www的网站| 寂寞人妻少妇视频99o| 国产一区亚洲一区在线观看| 欧美xxxx黑人xx丫x性爽| 男插女下体视频免费在线播放| 自拍偷自拍亚洲精品老妇| 亚洲精品亚洲一区二区| 亚洲成人av在线免费| 永久免费av网站大全| 三级国产精品欧美在线观看| 九九在线视频观看精品| 亚洲激情五月婷婷啪啪| 永久免费av网站大全| 成人特级av手机在线观看| 亚洲国产精品成人久久小说| 日韩欧美一区视频在线观看 | 麻豆久久精品国产亚洲av| 亚洲精品成人久久久久久| 青春草国产在线视频| 卡戴珊不雅视频在线播放| 欧美高清成人免费视频www| 婷婷六月久久综合丁香| 亚洲欧美日韩无卡精品| 日日摸夜夜添夜夜添av毛片| 超碰97精品在线观看| 男插女下体视频免费在线播放| 18禁裸乳无遮挡免费网站照片| 91在线精品国自产拍蜜月| 一个人看的www免费观看视频| 天堂网av新在线| 国产综合懂色| 午夜免费观看性视频| 久久久久久久亚洲中文字幕| 国产综合懂色| 日韩成人伦理影院| av在线亚洲专区| 国产色婷婷99| 国产又色又爽无遮挡免| 麻豆成人av视频| 亚洲三级黄色毛片| 久久亚洲国产成人精品v| 欧美+日韩+精品| 人妻制服诱惑在线中文字幕| 亚洲欧美精品自产自拍| 日韩精品青青久久久久久| 国产精品无大码| 2021天堂中文幕一二区在线观| 欧美日韩在线观看h| 国产亚洲91精品色在线| 国国产精品蜜臀av免费| 亚洲精品影视一区二区三区av| 韩国av在线不卡| 女的被弄到高潮叫床怎么办| 国产色爽女视频免费观看| 男女那种视频在线观看| 日韩欧美一区视频在线观看 | 嫩草影院入口| 亚洲人与动物交配视频| 午夜激情欧美在线| 亚洲av电影不卡..在线观看| 国产高清三级在线| 熟妇人妻久久中文字幕3abv| 精品国产露脸久久av麻豆 | 美女主播在线视频| 亚洲自偷自拍三级| 夫妻午夜视频| 嫩草影院入口| 亚洲成人久久爱视频| 中文字幕制服av| 国产在线一区二区三区精| 欧美日本视频| 十八禁网站网址无遮挡 | 插阴视频在线观看视频| 亚洲高清免费不卡视频| 精品欧美国产一区二区三| 成年人午夜在线观看视频 | 精品一区二区三卡| 亚洲国产日韩欧美精品在线观看| 国产毛片a区久久久久| 欧美日韩综合久久久久久| 大话2 男鬼变身卡| 亚洲熟女精品中文字幕| 久久久久久伊人网av| 国国产精品蜜臀av免费| 国语对白做爰xxxⅹ性视频网站| 狂野欧美白嫩少妇大欣赏| 久久久久久久国产电影| 亚洲人成网站在线观看播放| 国产精品一区二区性色av| 国产成人免费观看mmmm| 卡戴珊不雅视频在线播放| 免费看av在线观看网站| 天堂av国产一区二区熟女人妻| 欧美不卡视频在线免费观看| 日本与韩国留学比较| 亚洲内射少妇av| 亚洲四区av| 亚洲久久久久久中文字幕| av在线亚洲专区| 国产精品国产三级专区第一集| 一级a做视频免费观看| 午夜精品在线福利| 亚洲国产精品国产精品| 久久久久国产网址| 国产一区二区在线观看日韩| a级毛片免费高清观看在线播放| 国产成人aa在线观看| 国产日韩欧美在线精品| 亚洲精品aⅴ在线观看| 看非洲黑人一级黄片| av又黄又爽大尺度在线免费看| 高清欧美精品videossex| 国产在线一区二区三区精| 免费无遮挡裸体视频| 亚洲国产精品专区欧美| 欧美xxxx黑人xx丫x性爽| 熟妇人妻不卡中文字幕| 久久久久久久久久久免费av| 欧美日韩在线观看h| 综合色丁香网| 成人av在线播放网站| 欧美一级a爱片免费观看看| 日本-黄色视频高清免费观看| 99久久中文字幕三级久久日本| 精品久久久久久久末码| 又爽又黄a免费视频| 看免费成人av毛片| 亚洲欧洲国产日韩| 国国产精品蜜臀av免费| 少妇人妻精品综合一区二区| 看黄色毛片网站| 成人毛片a级毛片在线播放| 亚洲av福利一区| 简卡轻食公司| 午夜免费观看性视频| av.在线天堂| 亚洲一级一片aⅴ在线观看| 久久国内精品自在自线图片| 老师上课跳d突然被开到最大视频| 波野结衣二区三区在线| 美女主播在线视频| 亚洲av一区综合| 国产一级毛片在线| 欧美高清成人免费视频www| ponron亚洲| 国产黄a三级三级三级人| 熟女人妻精品中文字幕| 又黄又爽又刺激的免费视频.| 人体艺术视频欧美日本| 日韩中字成人| 国产免费视频播放在线视频 | 免费观看a级毛片全部| 最近中文字幕高清免费大全6| 亚洲精品国产av蜜桃| 色网站视频免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日日啪夜夜撸| 搡老妇女老女人老熟妇| 国产精品伦人一区二区| 亚洲精华国产精华液的使用体验| 国产乱人视频| 淫秽高清视频在线观看| 99热全是精品| 亚洲在久久综合| 国内揄拍国产精品人妻在线| 欧美成人一区二区免费高清观看| 女的被弄到高潮叫床怎么办| 黄色日韩在线| 国产成年人精品一区二区| 日韩亚洲欧美综合| 赤兔流量卡办理| av在线蜜桃| 国产成人福利小说| 久久久色成人| 国产免费又黄又爽又色| 国内少妇人妻偷人精品xxx网站| 精品久久国产蜜桃| 国产精品国产三级国产专区5o| 亚洲欧美成人精品一区二区| 久久热精品热| 一级毛片电影观看| 久久99热6这里只有精品| 久久久久免费精品人妻一区二区| 人人妻人人澡欧美一区二区| 超碰av人人做人人爽久久| 国产单亲对白刺激| 久久久久久久久久成人| 日韩欧美精品v在线| 久久99热这里只频精品6学生| 又爽又黄a免费视频| 免费观看在线日韩| 国产精品无大码| 久久久久精品久久久久真实原创| 91精品伊人久久大香线蕉| 亚洲图色成人| 人体艺术视频欧美日本| 爱豆传媒免费全集在线观看| 精品欧美国产一区二区三| 久久久精品免费免费高清| 国内少妇人妻偷人精品xxx网站| 日韩强制内射视频| 99热全是精品| 久久久久国产网址| www.色视频.com| 在线 av 中文字幕| 精品国内亚洲2022精品成人| 成人毛片60女人毛片免费| 色哟哟·www| 亚洲成人久久爱视频| 国产伦理片在线播放av一区| 26uuu在线亚洲综合色| videossex国产| 中文字幕久久专区| 亚洲精品aⅴ在线观看| 欧美性猛交╳xxx乱大交人| 又爽又黄无遮挡网站| 久久鲁丝午夜福利片| 97超视频在线观看视频| 国产麻豆成人av免费视频| 久久久久久久久久成人| 亚洲av男天堂| 成人性生交大片免费视频hd| 亚洲最大成人中文| 国产精品一区二区三区四区久久| 简卡轻食公司| 看免费成人av毛片| 久久久久久国产a免费观看| 如何舔出高潮| av在线蜜桃| 69人妻影院| av在线老鸭窝| 亚洲真实伦在线观看| 亚洲高清免费不卡视频| 国产色婷婷99| 26uuu在线亚洲综合色| 中文乱码字字幕精品一区二区三区 | 亚洲美女搞黄在线观看| av网站免费在线观看视频 | 狠狠精品人妻久久久久久综合| 人妻夜夜爽99麻豆av| 日本av手机在线免费观看| 亚洲不卡免费看| 22中文网久久字幕| 亚洲高清免费不卡视频| 人人妻人人看人人澡| 26uuu在线亚洲综合色| 精品人妻熟女av久视频| 国产一区二区三区综合在线观看 | videossex国产| 午夜免费观看性视频| 最近最新中文字幕大全电影3| 纵有疾风起免费观看全集完整版 | 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 欧美性猛交╳xxx乱大交人| 精品久久久久久久久久久久久| 男的添女的下面高潮视频| 一个人免费在线观看电影| 免费少妇av软件| 五月伊人婷婷丁香| 中文乱码字字幕精品一区二区三区 | 大话2 男鬼变身卡| 亚洲国产精品成人综合色| 小蜜桃在线观看免费完整版高清| 黄片wwwwww| 三级毛片av免费| 国产老妇女一区| 免费大片18禁| 成人午夜精彩视频在线观看| 最近最新中文字幕免费大全7| 99久国产av精品| 国产乱人视频| 18禁在线无遮挡免费观看视频| 欧美人与善性xxx| 久久久精品94久久精品| 亚洲美女搞黄在线观看| 蜜臀久久99精品久久宅男| 精品欧美国产一区二区三| 色综合色国产| 淫秽高清视频在线观看| av在线老鸭窝| 亚洲欧美成人综合另类久久久| 最近2019中文字幕mv第一页| 国产老妇女一区| 激情 狠狠 欧美| 精品久久久噜噜| 亚洲精品成人久久久久久| 亚洲人与动物交配视频| .国产精品久久| 69人妻影院| 日日摸夜夜添夜夜添av毛片| 亚洲av成人av| 国产在线一区二区三区精| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 老司机影院毛片| 国产av国产精品国产| 国产精品人妻久久久久久| 国产精品久久久久久精品电影| 国产亚洲一区二区精品| av在线蜜桃| 成人无遮挡网站| 精品少妇黑人巨大在线播放| 亚洲最大成人av| 成人综合一区亚洲| 99久久精品热视频| 久久草成人影院| 国产黄色免费在线视频| 国产老妇伦熟女老妇高清| 日本wwww免费看| 国产在视频线精品| 久久精品国产亚洲av涩爱| 国国产精品蜜臀av免费| 午夜久久久久精精品| 亚洲自拍偷在线| 久久人人爽人人片av| 国内精品美女久久久久久| 一区二区三区免费毛片| 熟妇人妻久久中文字幕3abv| 久久精品综合一区二区三区| 中文在线观看免费www的网站| 黄片wwwwww| 国产亚洲一区二区精品| 欧美日韩视频高清一区二区三区二| 99热这里只有是精品50| 日韩不卡一区二区三区视频在线| 麻豆av噜噜一区二区三区| 22中文网久久字幕| 亚洲综合精品二区| 国内精品一区二区在线观看| 免费看美女性在线毛片视频| 久久这里有精品视频免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品三级大全| 在线观看人妻少妇| 黄片wwwwww| 一级av片app| 国产黄频视频在线观看| 最近2019中文字幕mv第一页| 天堂中文最新版在线下载 | 国产精品国产三级国产av玫瑰| 欧美成人午夜免费资源| 2021天堂中文幕一二区在线观| 久久久久精品性色| 成年女人看的毛片在线观看| 午夜福利视频精品| 天天一区二区日本电影三级| 欧美性感艳星| 国产成人午夜福利电影在线观看| 男女边吃奶边做爰视频| 欧美一级a爱片免费观看看| 国产精品一及| 亚洲av电影不卡..在线观看| 深夜a级毛片| 久久精品久久久久久久性| 22中文网久久字幕| 一本一本综合久久| 欧美区成人在线视频| 两个人视频免费观看高清| 中文资源天堂在线| 男女国产视频网站| 久久精品熟女亚洲av麻豆精品 | 天天躁日日操中文字幕| 可以在线观看毛片的网站| 一级毛片 在线播放| 中文字幕人妻熟人妻熟丝袜美| 精华霜和精华液先用哪个| 欧美3d第一页| 汤姆久久久久久久影院中文字幕 | 亚洲人成网站在线观看播放| 国产成人福利小说| 国产成人a区在线观看| 美女黄网站色视频| 亚洲天堂国产精品一区在线| 久久精品久久久久久噜噜老黄| 极品教师在线视频| 亚洲乱码一区二区免费版| 少妇熟女欧美另类| 精品人妻偷拍中文字幕| 成人性生交大片免费视频hd| 1000部很黄的大片| 中文乱码字字幕精品一区二区三区 | 26uuu在线亚洲综合色| 国产v大片淫在线免费观看| 97人妻精品一区二区三区麻豆| 男人狂女人下面高潮的视频| 91aial.com中文字幕在线观看| 麻豆av噜噜一区二区三区| 汤姆久久久久久久影院中文字幕 | 亚洲国产精品国产精品| 日日撸夜夜添| 精品亚洲乱码少妇综合久久| 国产一区二区三区av在线| 啦啦啦啦在线视频资源| 欧美日韩在线观看h| 亚洲精品中文字幕在线视频 | av福利片在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品一区二区三区四区免费观看| 成人午夜高清在线视频| 99re6热这里在线精品视频| 观看美女的网站| 亚洲欧美日韩东京热| 久久97久久精品| 九草在线视频观看| 国产黄片视频在线免费观看| 久久精品国产亚洲av涩爱| 免费观看性生交大片5| 一个人看的www免费观看视频| 国产成人a∨麻豆精品| h日本视频在线播放| 精品久久久久久久末码| 99视频精品全部免费 在线| 国产亚洲精品久久久com| 寂寞人妻少妇视频99o| 国产老妇伦熟女老妇高清| 日韩一本色道免费dvd| 美女高潮的动态| 男女下面进入的视频免费午夜| 日本三级黄在线观看| 国产 亚洲一区二区三区 | 欧美日韩亚洲高清精品| 黄片无遮挡物在线观看| 超碰97精品在线观看| 五月伊人婷婷丁香| 国产久久久一区二区三区| av免费在线看不卡| 免费少妇av软件| 久久人人爽人人爽人人片va| 亚洲人成网站在线观看播放| 日本av手机在线免费观看| 岛国毛片在线播放| 最近手机中文字幕大全| 欧美精品一区二区大全| 亚洲精品成人av观看孕妇| 2022亚洲国产成人精品| 如何舔出高潮| 七月丁香在线播放| 国产淫语在线视频| 久久久精品94久久精品| 人妻一区二区av| 丝袜美腿在线中文| 亚洲成人av在线免费| 国内精品宾馆在线| 国产爱豆传媒在线观看| 人妻系列 视频| 亚洲精品久久久久久婷婷小说| 午夜激情福利司机影院| 亚洲精品aⅴ在线观看| 亚洲精品成人av观看孕妇| 久久国产乱子免费精品| 亚洲美女搞黄在线观看| 亚洲av福利一区| 国产亚洲精品av在线| 国产精品久久视频播放| 国产一级毛片七仙女欲春2| 国产亚洲精品久久久com| 日韩av在线免费看完整版不卡| 美女高潮的动态| 国产伦精品一区二区三区视频9| 观看免费一级毛片| 激情五月婷婷亚洲|