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

    融水小型豬線粒體DNA結(jié)構(gòu)和系統(tǒng)進(jìn)化的分析

    2015-11-27 05:27:02徐嘉悅施赫赫王緒敏余細(xì)勇唐小江
    關(guān)鍵詞:融水豬種密碼子

    徐嘉悅,施赫赫,王緒敏,高 翔,余細(xì)勇,于 軍,唐小江

    (1.中國(guó)科學(xué)院北京基因研究所,北京 100101;2.廣東省醫(yī)學(xué)實(shí)驗(yàn)動(dòng)物中心,廣東佛山 528248;3.廣東貝格生物科技有限公司,廣東佛山 528100;4.南京大學(xué)模式動(dòng)物研究所,南京 210061;5.廣東省醫(yī)學(xué)科學(xué)院,廣州 510060)

    融水小型豬是源自廣西苗寨的小黑豬,具有體型微小、性情溫馴、產(chǎn)仔力強(qiáng)、母性較好等特點(diǎn),具備培育成標(biāo)準(zhǔn)化實(shí)驗(yàn)用小型豬的有利特點(diǎn)。我們于2012年將該豬種源引到廣東三水繁育,建立了基礎(chǔ)種群和的系譜檔案,已繁殖出F1代和F2代,參照北京市地方標(biāo)準(zhǔn)[1]初步建立了融水小型豬的質(zhì)量控制標(biāo)準(zhǔn)。為了進(jìn)一步掌握其遺傳背景,本文對(duì)線粒體DNA進(jìn)行了分析。

    1 材料和方法

    1.1 樣本的收集和DNA的提取

    融水小型豬由廣東貝格生物科技有限公司(廣東省醫(yī)學(xué)實(shí)驗(yàn)動(dòng)物中心小型豬基地)提供。用無(wú)菌肝素鈉抗凝采血管取融水小型豬外周靜脈血10 mL,用Qiagen公司的QlAamp DNA試劑盒按操作說(shuō)明書(shū)提取DNA。

    1.2 序列的測(cè)序和組裝

    采用全基因組鳥(niǎo)槍法的測(cè)序策略,使用Illumina基因組分析測(cè)序技術(shù)進(jìn)行測(cè)序研究。為了提高組裝的質(zhì)量,對(duì)融水小型豬選取了插入片段大小為300 bp和500 bp的paired-end文庫(kù)進(jìn)行末端測(cè)序,用基于 Brujin圖算法[2]的 SOAPdenovo軟件(http://soap.genomeics.org.cn)對(duì)測(cè)序得到的短片段進(jìn)行拼接組裝[3-4],得到融水小型豬的完整線粒體DNA序列。

    1.3 線粒體DNA的注釋和分析

    為了確定融水小型豬的線粒體結(jié)構(gòu),首先以Susscrofa的線粒體序列為參考序列,采用BLAST軟件[4]完成了Susscrofa與融水小型豬的比對(duì),確定線粒體的編碼蛋白基因、rRNA和 tRNA,用軟件CodonW對(duì)編碼蛋白基因的密碼子使用情況進(jìn)行統(tǒng)計(jì)。另外用 tRNAscan-SE v.1.21軟件[5]對(duì)線粒體tRNA的二級(jí)結(jié)構(gòu)及各個(gè)tRNA的反義密碼子進(jìn)行預(yù)測(cè)。用clustal W[6]軟件將融水小型豬的編碼蛋白區(qū)序列與GenBank中的15頭國(guó)內(nèi)外的豬線粒體編碼蛋白區(qū)序列(表1)進(jìn)行比對(duì)。

    用Dnasp5.10[7]軟件計(jì)算實(shí)驗(yàn)材料的單倍型、遺傳多樣性和遺傳分化的相關(guān)數(shù)值[8]?;谝陨系膯伪缎途垲?lèi)的分類(lèi),采用 Tajima’s D[9]和 Fu’s[10]兩個(gè)中性檢驗(yàn)?zāi)P蛯?duì)16個(gè)豬的單倍型分布進(jìn)行中性檢驗(yàn),其中種群遺傳格局及變異的檢測(cè)在Arkequin3.11[11]中 實(shí) 施。 使 用 分 子 變 異 分 析(AMOVA)的方法以10000次重復(fù)隨機(jī)抽樣單倍型重排后進(jìn)行顯著性檢驗(yàn),以此來(lái)評(píng)價(jià)種群遺傳變異水平與種群地理位置的相關(guān)性,使用統(tǒng)計(jì)量將遺傳差異分為3個(gè)不同層次的地理等級(jí)。用Network 4.6(http://www.fluxus-engineering.com/sharenet.htm)對(duì)16個(gè)豬種進(jìn)行聚類(lèi),并用MEGA 5.0[12]的最大似然法和鄰接法對(duì)16個(gè)豬種的編碼蛋白區(qū)進(jìn)行系統(tǒng)進(jìn)化分析。

    2 結(jié)果

    2.1 融水小型豬線粒體結(jié)構(gòu)

    融水小型豬與其他豬的結(jié)構(gòu)框架相同,都是由13個(gè)編碼蛋白,22個(gè)tRNA和2個(gè)rRNA組成(表2)。線粒體全長(zhǎng)為16888 bp。線粒體DNA的核苷酸組成為A>C>T>G。線粒體的GC含量總是低于AT含量。

    2.2 編碼蛋白基因

    融水小型豬由13個(gè)基因構(gòu)成,除基因ND6位于線粒體的輕鏈上之外,其余12個(gè)基因包括ND1、ND2、ND3、ND4、ND4L、ND5、COX1、COX2、COX3、CYTB、ATP6、ATP8都在線粒體的重鏈(表2)。這13個(gè)基因的總長(zhǎng)為11255 bp,其中基因ATP6和基因ATP8有42 bp的重疊區(qū)域,基因ND4L和基因ND4有7 bp的重疊區(qū)域,基因ND5和基因ND6有17 bp的重疊區(qū)域。位于融水小型豬重鏈的12個(gè)編碼蛋白基因的堿基組分為A>C >T>G,其中ND1、ATP6、COX3、ND3、ND4、ND5、ND2,CYTB8 個(gè)基因的堿基組份是一致的:A >C >T >G,而基因 COX1、COX2、ATP8、ND4L的堿基組份為A>T>C>G;另外,位于融水小型豬輕鏈的ND6基因的堿基組份為C>A>T>G。

    2.3 非編碼蛋白區(qū)域:D-loop控制區(qū)和基因間區(qū)

    在拼接的融水小型豬線粒體中,D-loop區(qū)的長(zhǎng)度為1254 bp,tRNA-Pro和tRNA-Phe排布于D-loop區(qū)的兩端。在D-loop區(qū)堿基組份為A>C>T>G,與線粒體DNA一致,AT的含量高于GC含量。此外,融水小型豬線粒體DNA中非編碼的基因間區(qū)長(zhǎng)度為654 bp,這些基因間區(qū)的主要原因是線粒體中tRNA簇的重排。其中最長(zhǎng)的一段基因間區(qū)是位于ND1和tRNA之間的區(qū)域,長(zhǎng)為139 bp。基因間區(qū)為高級(jí)別的系統(tǒng)進(jìn)化研究提供了可能。

    表1 15個(gè)豬種線粒體DNA的基本信息及其在NCBI數(shù)據(jù)庫(kù)中的檢索號(hào)Tab.1 List of mitochondrial genome sequences of Swine breeds and their accession numbers as retrieved from NCBI database

    表2 融水小型豬線粒體DNA的結(jié)構(gòu)特征Tab.2 Location of feature in the mitochondrial DNA of Rongshui miniature pig

    2.4 核糖體RNA

    線粒體中包括兩個(gè)核糖體 RNA,分別為12S rRNA and 16S rRNA(表1),12S rRNA的長(zhǎng)度為959 bp,16S rRNA的長(zhǎng)度為1512 bp。16S rRNA編碼于纈氨酸和亮氨酸之間,其堿基組份為A>T>C>G;12S rRNA編碼于丙氨酸和纈氨酸之間,其堿基組份為A>C>T>G。12S rRNA和16S rRNA的GC含量均在40%左右,同樣保持了GC含量低于AT含量這一特點(diǎn)。

    2.5 轉(zhuǎn)運(yùn)RNA

    對(duì)融水小型豬線粒體,我們預(yù)測(cè)了其22個(gè)轉(zhuǎn)運(yùn)RNA的二級(jí)結(jié)構(gòu),并確定了各轉(zhuǎn)運(yùn)RNA的反義密碼子(如表2所示)。這22個(gè)轉(zhuǎn)運(yùn)RNA的GC含量都低于AT含量,其中精氨酸的GC含量最低,約21%;甲硫氨酸的GC含量最高,約47%。

    融水小型豬線粒體DNA的密碼子使用情況(表3)。包括終止密碼子在內(nèi),共有3772個(gè)密碼子,其中使用密碼子CUA的數(shù)目最多,為306個(gè),而AGG和CGG兩個(gè)密碼子沒(méi)有被使用。在密碼子的統(tǒng)計(jì)中,密碼子第一位核苷酸使用的頻率為A>C>T>G;而在密碼子第二位置,堿基使用頻率為T(mén)>C>A>G;在密碼子的第三個(gè)位置,堿基使用頻率為A>C>U>G。眾所周知,密碼子在不同位置使用頻率的差異是由于密碼子使用的偏好性造成的,上述的密碼子使用情況也為研究線粒體DNA的進(jìn)化提供了有用的信息[13]。

    2.6 同義突變和非同義突變

    同義突變和非同義突變之間的比例(Ka/Ks值)可用于分析不同物種之間的的系統(tǒng)進(jìn)化關(guān)系[14]。通過(guò)對(duì)16個(gè)豬13個(gè)編碼蛋白區(qū)的所有序列進(jìn)行比對(duì),構(gòu)建了一個(gè)參考序列(Ref),分別計(jì)算每頭豬的各個(gè)編碼蛋白基因序列與Ref的Ka/Ks比值。結(jié)果如封2圖1所示,16個(gè)豬種在編碼蛋白區(qū)的Ka/KS值均小于1,說(shuō)明所有13個(gè)編碼蛋白基因都在承受著正選擇,其中基因COX1的選擇壓力最大,這意味著該基因正在承受著很強(qiáng)的純化選擇;而ND5的選擇壓力相對(duì)較小,說(shuō)明該基因是在進(jìn)化中舍棄一些較為不利的突變。

    2.7 系統(tǒng)進(jìn)化分析

    使用 MEGA軟件中的 Kimura’s two-parameter模型(K2P)計(jì)算16個(gè)豬種的編碼蛋白區(qū)的遺傳距離。融水小型豬與其他豬的遺傳距離最遠(yuǎn),堿基的平均遺傳距離為1.141。而對(duì)總體的遺傳距離統(tǒng)計(jì)中,物種間以恒定速率進(jìn)行核苷酸替換的標(biāo)準(zhǔn)差僅為0.0015?;谏鲜鲞z傳距離矩陣,通過(guò)MEGA軟件中的最大似然法和鄰接法對(duì)16頭豬的編碼蛋白區(qū)構(gòu)建系統(tǒng)進(jìn)化樹(shù)(圖2),從圖2可以看出,16個(gè)豬種可以分成四個(gè)類(lèi)群,即歐洲群(E群)、亞洲1群(A1群)、亞洲2群(A2群)、亞洲3群(A3群)。其中,E群主要是由歐洲豬種組成,其中意大利的野豬是相對(duì)古老的豬種;A1群是由中國(guó)、日本、韓國(guó)等亞洲豬種及Berkshire豬共同組成,在這一分支中,也可以證明豬的進(jìn)化是由亞洲起源后分化到歐洲的[15];A2群是由滇南7號(hào)、海南野豬和蘭嶼豬構(gòu)成,這一結(jié)論也驗(yàn)證了之前的研究,即云南豬種和臺(tái)灣地區(qū)蘭嶼豬是古老的豬種[15]。在系統(tǒng)進(jìn)化分析中最值得關(guān)注的是融水小型豬,可以看出他比蘭嶼豬更為古老。

    2.8 線粒體序列的單倍型分析及中性檢驗(yàn)

    系統(tǒng)進(jìn)化分析是對(duì)過(guò)去已經(jīng)發(fā)生的進(jìn)化時(shí)間的模擬分析,為了驗(yàn)證系統(tǒng)進(jìn)化分析結(jié)果的準(zhǔn)確性,用單倍型的方法,對(duì)16個(gè)豬種的進(jìn)化關(guān)系做了進(jìn)一步分析。用Dnasp5.10[7]分析16個(gè)豬種的線粒體DNA的編碼區(qū)序列,共得到14個(gè)單倍型,但單倍型的多樣性指數(shù)為0.596。將分析后的單倍型數(shù)據(jù)利用Network軟件作圖(封2圖3)。14個(gè)單倍型可以分成3大類(lèi):漢普郡豬、蘇塞克斯家養(yǎng)豬、長(zhǎng)白豬、蘇塞克斯豬、杜洛克豬、意大利野豬、歐洲伯克希爾豬、聚為一類(lèi)(類(lèi)群1);濟(jì)州8#豬、濟(jì)州10#豬、韓國(guó)野豬、梅山豬和臺(tái)灣地區(qū)野豬聚為一類(lèi)(類(lèi)群2);蘭嶼豬、滇南7號(hào)豬、海南野豬、融水小型豬聚為一類(lèi)(類(lèi)群3;其中類(lèi)群1主要是歐洲豬種,類(lèi)群2主要以亞洲的韓國(guó)和臺(tái)灣地區(qū)豬種為主,類(lèi)群3主要是中國(guó)4頭古老的豬種??梢钥闯鯪etwork作圖結(jié)果與系統(tǒng)進(jìn)化分析的結(jié)果大致是一致的,這也進(jìn)一步驗(yàn)證了融水小型豬是一個(gè)非常古老的豬種。

    基于單倍型的分析結(jié)果,分別對(duì)3個(gè)類(lèi)群進(jìn)行了中性檢驗(yàn)(如表4)。類(lèi)群3不但具有較高的單倍型多樣性(h=0.899),也具有較高的核苷酸多樣性(Pi=0.03417),表明類(lèi)群3中的豬種仍具有豐富的遺傳多樣性,存在著較高的進(jìn)化潛力。對(duì)于類(lèi)群1和類(lèi)群2,單倍型多樣性較高,但是核苷酸多樣性很低。說(shuō)明這兩個(gè)類(lèi)群的豬種經(jīng)歷過(guò)近期擴(kuò)張或者進(jìn)化歷史很短[17]。本研究中類(lèi)群1和類(lèi)群3的Tajin a’s D與Fs值為負(fù)值且差異顯著,這表明類(lèi)群1和類(lèi)群3中的豬種在歷史上都發(fā)生了種群擴(kuò)張。

    表3 融水小型豬密碼子Tab.3 Codon usage for Rongshui miniature pig

    圖2 16個(gè)豬種的系統(tǒng)進(jìn)化分析(最大似然法和鄰接法)Fig.2 The phylogenetic analysis of 16 pig breeds(NJ,ML analyses)

    表4 對(duì)16個(gè)豬種編碼蛋白基因的中性檢驗(yàn)和遺傳多樣性分析Tab.4 Neutrality test and genetic diversity for main haplogroups in coding region for 16 swine mitochondrial genomes

    3 討論

    線粒體DNA存在于細(xì)胞質(zhì)線粒體機(jī)質(zhì)中的閉合環(huán)狀雙鏈超螺旋DNA分子,是獨(dú)立于細(xì)胞核染色體外的基因組,具有自我復(fù)制、轉(zhuǎn)錄和編碼功能,是母系遺傳,但同時(shí)受核 DNA的控制[2]。與核酸DNA相比,線粒體DNA具有分子結(jié)構(gòu)簡(jiǎn)單,以母性方式遺傳,核苷酸歧義大,進(jìn)化速度快等特點(diǎn)。因此,使用線粒體研究哺乳動(dòng)物群體遺傳多樣性和種內(nèi)、種間親緣關(guān)系等,一直受到很多學(xué)者的關(guān)注[15,18-20]。本研究采用新一代測(cè)序技術(shù),對(duì)融水小型豬進(jìn)行測(cè)序分析,拼接成完整的線粒體DNA序列,從分子遺傳學(xué)角度評(píng)估了融水小型豬的結(jié)構(gòu)質(zhì)特性、系統(tǒng)進(jìn)化關(guān)系及其遺傳多樣性。為融水小型豬的實(shí)驗(yàn)動(dòng)物化提供科學(xué)依據(jù)。

    隨著新一代測(cè)序技術(shù)的發(fā)展,越來(lái)越多的學(xué)者對(duì)哺乳動(dòng)物遺傳機(jī)制的研究不再僅僅停留在核DNA水平上,而是更多的綜合細(xì)胞質(zhì)DNA進(jìn)行分析。線粒體DNA結(jié)構(gòu)簡(jiǎn)單、穩(wěn)定,是母系遺傳方式,在遺傳過(guò)程中未曾發(fā)生重組。

    本文使用新一代測(cè)序技術(shù),用Illumina測(cè)序平臺(tái)對(duì)融水小型豬進(jìn)行從頭測(cè)序,使用SOAPdenovo拼接成完整線粒體序列并進(jìn)行注釋。融水小型豬的線粒體DNA成環(huán)形,總長(zhǎng)為16888 bp,由13個(gè)基因,22個(gè)tRNA和2個(gè)rRNA組成。其中GC的平均含量約為44%。線粒體的GC含量總是低于AT含量,這一結(jié)論與其他物種一致,如人的GC含量為44.4%[21],小鼠的 GC 含量為 26.6%[22]。本文結(jié)合公共數(shù)據(jù)庫(kù)發(fā)布的15個(gè)豬種線粒體序列,以線粒體的編碼蛋白區(qū)為研究對(duì)象,對(duì)豬種的遺傳多樣性進(jìn)行了研究,進(jìn)而對(duì)豬種的起源和進(jìn)化進(jìn)行探討。從系統(tǒng)進(jìn)化分析中我們可以得出,融水小型豬在系統(tǒng)進(jìn)化樹(shù)中成為一個(gè)新的獨(dú)立外群,我們認(rèn)為融水小型豬是一個(gè)古老于蘭嶼豬的傳統(tǒng)中國(guó)豬種。為了進(jìn)一步驗(yàn)證這個(gè)結(jié)論,我們對(duì)系統(tǒng)進(jìn)化樹(shù)中的三個(gè)類(lèi)群進(jìn)行了遺傳多樣性分析,Network的結(jié)果與系統(tǒng)進(jìn)化分析的結(jié)果是吻合的。由滇南7號(hào)、海南野豬、蘭嶼豬和融水小型豬構(gòu)成的類(lèi)群3中有4個(gè)單倍型,其具有較高的單倍型多樣性和核苷酸多樣性顯示出類(lèi)群3具有豐富的遺傳多樣性。

    Tajin a‘s D與Fu’s Fs中性檢驗(yàn)常用以推測(cè)種群擴(kuò)張的歷史事件。如果Tajin a’s D與Fs值呈負(fù)值,且在統(tǒng)計(jì)學(xué)上達(dá)到顯著標(biāo)準(zhǔn),則說(shuō)明序列中含有比中性進(jìn)化模型中更多的核苷位點(diǎn)變化,可能預(yù)示著該種群曾經(jīng)有過(guò)擴(kuò)張的歷史。融水小型豬的Tajin a’s D與Fs值為負(fù)值且差異顯著,說(shuō)明了類(lèi)群3的古老豬種在長(zhǎng)期進(jìn)化過(guò)程中還在與其他的豬種還在發(fā)生基因交流和種群擴(kuò)張。

    另外,在系統(tǒng)進(jìn)化分析時(shí),蘭嶼豬、滇南7號(hào)豬、海南野豬與融水小型豬分成了不同的分支,而在Network分析中成為一個(gè)類(lèi)群,出現(xiàn)這種現(xiàn)象的原因可能是融水小型豬在長(zhǎng)期的進(jìn)化過(guò)程中,發(fā)生了基因交流。中性檢驗(yàn)的結(jié)果不但從統(tǒng)計(jì)水平上驗(yàn)證了融水小型豬在豬線粒體發(fā)展中的祖先地位,而且還解釋了融水小型豬作為祖先豬種是如何與其他豬種發(fā)生基因流和種群擴(kuò)張的。遺傳多樣性分析進(jìn)一步提示融水小型豬是中國(guó)傳統(tǒng)豬種中一種古老的小型豬。

    本文對(duì)融水小型豬DNA進(jìn)行了初步的分析,限于篇幅和GenBank的數(shù)據(jù),未與哥廷根小型豬、西藏小型豬、巴馬小型豬等其他實(shí)驗(yàn)用小型豬的線粒體DNA進(jìn)行比較,有待進(jìn)一步研究。

    [1] 北京市質(zhì)量技術(shù)監(jiān)督局.實(shí)驗(yàn)用小型豬:微生物學(xué)等級(jí)及監(jiān)測(cè)DB11T 828.1-2011[S];寄生蟲(chóng)學(xué)等級(jí)及監(jiān)測(cè) DB11T 828.2-2011[S];遺傳質(zhì)量控制 DB11T 828.3 -2011[S];病理學(xué)診斷規(guī)范 DB11T 828.4-2011[S];配合飼料 DB11T 828.5-2011[S];環(huán)境及設(shè)施 DB11T828.6 -2011[S],2011-11-10發(fā)布,2012-03-01實(shí)施.

    [2] Pevzner,P.A.,H.Tang,M.S.Waterman,An Eulerian path approach to DNA fragment assembly[J].Proc Natl Acad Sci U S A,2001,98(17):74-85.

    [3] Angiero and Natalie.,Pigs Prove to Be Smart,if Not Vain.The New York Times,2009.

    [4] Altschul,S.F.,Gish,W.,Miller,W.,et al.,Basic local alignment search tool[J].J Mol Biol,1990,215(3):403-410.

    [5] Lowe,T.M.S.R.Eddy,tRNAscan-SE:a program for improved detection of transfer RNA genes in genomic sequence[J].Nucleic Acids Res,1997,25(5):95 -98.

    [6] Thompson,J.D.,D.G.Higgins,T.J.Gibson,CLUSTAL W:improving the sensitivity ofprogressive multiple sequence alignment through sequence weighting,position-specific gap penalties and weight matrix choice[J].Nucleic Acids Res,1994,22(22):73-80.

    [7] Rozas,J.,Sanchez-DelBarrio,J.C.,Messeguer,X.,et al.,DnaSP,DNA polymorphism analyses by the coalescent and other methods[J].Bioinformatics,2003.19(18):49 -67.

    [8] Kimura,M.,Molecular evolutionary clock and the neutral theory[J].J Mol Evol,1987.26(1-2):24-33.

    [9] Tajima,F(xiàn).,Statistical method for testing the neutral mutation hypothesis by DNA polymorphism[J].Genetics,1989,123(3):58-63.

    [10] Fu,Y.X.,Statistical tests of neutrality of mutations against population growth,hitchhiking and background selection[J].Genetics,1997,147(2):91 -105.

    [11] Excoffier L,L.L., Schneider S, Arlequin ver.3.0:an integrated software package for population genetics data analysis[J].Evolutionary Bioinformatics Online,2005,1:47 -50.

    [12] Tamura,K.,Peterson,D.,Peterson,N.,et al.,MEGA5:molecular evolutionary genetics analysis using maximum likelihood, evolutionarydistance, and maximum parsimony methods[J].Mol Biol Evol.28(10):73 -91.

    [13] Dubey,B., MeganathanP.R., HaqueI., DNA minibarcoding:an approach for forensic identification of some endangered Indian snake species[J].Forensic Sci Int Genet,5(3):18-24.

    [14] Yang, Z. R. Nielsen, Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models[J].Mol Biol Evol,2000,17(1):32 -43.

    [15] Chen,Huang.,H.L,Yang H.Y.,et al.Mitochondrial genome of Taiwan pig (Susscrofa) [J]. African Journalof Biotechnology,2011,10(13):2556-2561.

    [16] Wu,G.S.,Yao,Y.G.,Qu,K.X.,et al.,Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia[J].Genome Biol,2007,8(11):245 -253.

    [17] Femandez,A.I.,Alves,E.,Ovilo,C.,et al.,Divergence time estimates of East Asian and European pigs based on multiple near complete mitochondrial DNA sequences[J].Anim Genet,42(1):86-88.

    [18] Kijas,J.M.Andersson L.,A phylogenetic study of the origin of the domestic pig estimated from the near-complete mtDNA genome[J].J Mol Evol,2001,52(3):302 -8.

    [19] Femandez,A.I.,Alves,E.,Ovilo,C.,et al.,Divergence time estimates of East Asian and European pigs based on multiple near complete mitochondrial DNA sequences[J].Anim Genet,42(1):86-89.

    [20] Yang,S.,Zhang,H.,Mao,H.,et al.,The local origin of the Tibetan pig and additional insights into the origin of Asian pigs[J].PLoS One,6(12):218 -215.

    [21] Anderson,S.,Bankier,A.T.,Barrell,B.G.,et al.,Sequence and organization of the human mitochondrial genome[J].Nature,1981,290(5806):45 -65.

    [22] Bibb,M.J.,Van Etten,R.A.,Wright,C.T.,et al.,Sequence and gene organization of mouse mitochondrial DNA[J].Cell,1981,26(2 Pt 2):67 -80.

    猜你喜歡
    融水豬種密碼子
    貴州地方特色豬種培育
    可視化中國(guó)地方豬種地理分布圖
    我在融水唱貝江
    河南精旺豬種改良有限公司
    國(guó)外豬種大量引進(jìn)導(dǎo)致本土豬種瀕臨滅絕的原因分析及應(yīng)對(duì)建議(以玉山黑豬為例)
    融水糯米柚種植栽培管理與技術(shù)探究
    密碼子與反密碼子的本質(zhì)與拓展
    10種藏藥材ccmFN基因片段密碼子偏好性分析
    中成藥(2018年7期)2018-08-04 06:04:10
    例舉高中數(shù)學(xué)解題切入點(diǎn)的找尋
    廣西苗鄉(xiāng)“蘆笙鬥馬節(jié)”驚險(xiǎn)刺激
    老司机福利观看| 91精品三级在线观看| 亚洲人成77777在线视频| 制服诱惑二区| 亚洲欧洲日产国产| 色婷婷久久久亚洲欧美| 啦啦啦免费观看视频1| 黄色视频不卡| 国产极品粉嫩免费观看在线| 久久人人爽av亚洲精品天堂| 另类精品久久| 少妇被粗大的猛进出69影院| 这个男人来自地球电影免费观看| 久久久久久久久免费视频了| 国产一级毛片在线| 深夜精品福利| 久久久国产成人免费| 伦理电影免费视频| 一区二区三区激情视频| 人人妻人人爽人人添夜夜欢视频| 亚洲熟女精品中文字幕| 中文字幕av电影在线播放| 热re99久久精品国产66热6| 国产免费av片在线观看野外av| 丝袜美腿诱惑在线| 可以免费在线观看a视频的电影网站| 欧美日韩国产mv在线观看视频| 亚洲精品国产av蜜桃| 一本一本久久a久久精品综合妖精| 麻豆av在线久日| 777米奇影视久久| 亚洲专区字幕在线| 午夜福利视频在线观看免费| 99精国产麻豆久久婷婷| 精品欧美一区二区三区在线| 丝袜脚勾引网站| 欧美97在线视频| 天堂中文最新版在线下载| 国产成人免费观看mmmm| 欧美性长视频在线观看| 一区福利在线观看| 亚洲人成电影观看| 亚洲精品久久午夜乱码| 黑人巨大精品欧美一区二区蜜桃| 久久中文字幕一级| 可以免费在线观看a视频的电影网站| 国产成+人综合+亚洲专区| 99久久人妻综合| 亚洲av日韩精品久久久久久密| 亚洲熟女精品中文字幕| 王馨瑶露胸无遮挡在线观看| 亚洲 国产 在线| 国产在线免费精品| 免费不卡黄色视频| 久久精品国产综合久久久| 嫩草影视91久久| 色94色欧美一区二区| 国产成人av教育| 美女福利国产在线| 91成人精品电影| 免费高清在线观看日韩| 欧美一级毛片孕妇| 久久国产精品影院| 国产在线视频一区二区| 99久久精品国产亚洲精品| 午夜91福利影院| 女人高潮潮喷娇喘18禁视频| 精品福利永久在线观看| 丰满迷人的少妇在线观看| 亚洲中文字幕日韩| 国产亚洲欧美在线一区二区| 18禁黄网站禁片午夜丰满| 久久 成人 亚洲| 亚洲全国av大片| 久久精品人人爽人人爽视色| 国产高清国产精品国产三级| 飞空精品影院首页| 午夜福利免费观看在线| 夜夜夜夜夜久久久久| 精品国产国语对白av| 精品人妻一区二区三区麻豆| 欧美精品一区二区大全| 亚洲情色 制服丝袜| 动漫黄色视频在线观看| 99精品欧美一区二区三区四区| 久久精品亚洲av国产电影网| 啦啦啦视频在线资源免费观看| 狠狠精品人妻久久久久久综合| 国产在线免费精品| 日韩大片免费观看网站| 一级黄色大片毛片| 国产欧美日韩精品亚洲av| 亚洲一卡2卡3卡4卡5卡精品中文| netflix在线观看网站| 精品亚洲成a人片在线观看| 日韩欧美免费精品| 欧美少妇被猛烈插入视频| 精品少妇一区二区三区视频日本电影| 国产欧美日韩精品亚洲av| 操出白浆在线播放| netflix在线观看网站| 波多野结衣av一区二区av| 亚洲欧美一区二区三区黑人| 久久午夜综合久久蜜桃| 飞空精品影院首页| 男女之事视频高清在线观看| 国产精品一区二区在线不卡| 欧美+亚洲+日韩+国产| 午夜福利免费观看在线| 日本猛色少妇xxxxx猛交久久| 黑人巨大精品欧美一区二区mp4| 欧美激情 高清一区二区三区| 久久九九热精品免费| 一级毛片女人18水好多| 欧美av亚洲av综合av国产av| 这个男人来自地球电影免费观看| 男男h啪啪无遮挡| av天堂在线播放| 肉色欧美久久久久久久蜜桃| 黑人猛操日本美女一级片| 性色av一级| 在线天堂中文资源库| 国产精品秋霞免费鲁丝片| 咕卡用的链子| 午夜福利乱码中文字幕| 青草久久国产| 咕卡用的链子| 亚洲成人国产一区在线观看| 亚洲av男天堂| 又紧又爽又黄一区二区| av有码第一页| 一级片'在线观看视频| 少妇 在线观看| 中文字幕人妻丝袜一区二区| 各种免费的搞黄视频| a 毛片基地| 久久av网站| 免费在线观看黄色视频的| 亚洲欧美一区二区三区久久| 久久久欧美国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲午夜精品一区,二区,三区| 亚洲精品一二三| 69精品国产乱码久久久| 在线 av 中文字幕| 日韩,欧美,国产一区二区三区| 热99re8久久精品国产| 一本久久精品| kizo精华| 日韩大片免费观看网站| 婷婷色av中文字幕| 日韩视频一区二区在线观看| 国产一区二区在线观看av| 日日摸夜夜添夜夜添小说| 亚洲成人免费电影在线观看| 久久天堂一区二区三区四区| 欧美xxⅹ黑人| 久久99一区二区三区| 最新在线观看一区二区三区| 少妇猛男粗大的猛烈进出视频| 国产成人欧美| 亚洲第一青青草原| 视频在线观看一区二区三区| 亚洲国产中文字幕在线视频| 精品久久久精品久久久| 飞空精品影院首页| 亚洲精品一二三| 男女免费视频国产| 中文欧美无线码| 亚洲国产av新网站| 久久 成人 亚洲| 建设人人有责人人尽责人人享有的| tube8黄色片| 性高湖久久久久久久久免费观看| 欧美成狂野欧美在线观看| 中文字幕色久视频| 国产免费视频播放在线视频| 亚洲国产看品久久| 亚洲一卡2卡3卡4卡5卡精品中文| 老司机在亚洲福利影院| 亚洲av电影在线进入| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品一卡2卡三卡4卡5卡 | 亚洲欧洲日产国产| 乱人伦中国视频| 亚洲av欧美aⅴ国产| 黑人欧美特级aaaaaa片| 国产日韩欧美亚洲二区| 桃红色精品国产亚洲av| av电影中文网址| kizo精华| 桃花免费在线播放| 久久人人爽人人片av| 日韩熟女老妇一区二区性免费视频| 2018国产大陆天天弄谢| 久久人人97超碰香蕉20202| 最近最新中文字幕大全免费视频| 热re99久久国产66热| 精品国产一区二区三区四区第35| 三级毛片av免费| 丝袜脚勾引网站| 美女高潮到喷水免费观看| 亚洲第一欧美日韩一区二区三区 | 乱人伦中国视频| 久久天躁狠狠躁夜夜2o2o| 国精品久久久久久国模美| 欧美在线黄色| 国产一区二区三区av在线| 日韩免费高清中文字幕av| 久久中文看片网| 两人在一起打扑克的视频| 国产精品免费视频内射| 91成年电影在线观看| 国产成人啪精品午夜网站| 12—13女人毛片做爰片一| 久久久久久久大尺度免费视频| 大陆偷拍与自拍| 高清在线国产一区| 午夜免费成人在线视频| 成人黄色视频免费在线看| 国产在视频线精品| av有码第一页| 国产一区二区 视频在线| 91成人精品电影| 亚洲美女黄色视频免费看| 欧美日韩国产mv在线观看视频| 91麻豆精品激情在线观看国产 | 国产主播在线观看一区二区| 亚洲欧美成人综合另类久久久| 久热爱精品视频在线9| 亚洲人成77777在线视频| 搡老乐熟女国产| 在线观看人妻少妇| av福利片在线| 性色av一级| 99精国产麻豆久久婷婷| 天堂8中文在线网| 国产日韩欧美在线精品| 久久久精品国产亚洲av高清涩受| 老汉色∧v一级毛片| 亚洲一码二码三码区别大吗| 搡老熟女国产l中国老女人| 国产伦理片在线播放av一区| 国产成人啪精品午夜网站| 国产在线一区二区三区精| 国产精品偷伦视频观看了| 欧美国产精品va在线观看不卡| 老汉色av国产亚洲站长工具| 制服人妻中文乱码| av欧美777| 精品一区二区三区四区五区乱码| 成人av一区二区三区在线看 | 亚洲精品国产区一区二| 免费看十八禁软件| 老司机在亚洲福利影院| 亚洲七黄色美女视频| 亚洲 欧美一区二区三区| 亚洲精品久久成人aⅴ小说| 日本wwww免费看| 麻豆乱淫一区二区| 久久人妻熟女aⅴ| 亚洲精品av麻豆狂野| 久久精品国产亚洲av香蕉五月 | 国产熟女午夜一区二区三区| 99热全是精品| 国产免费福利视频在线观看| 一本久久精品| 国产精品自产拍在线观看55亚洲 | 国产成人啪精品午夜网站| 99国产极品粉嫩在线观看| 国产精品成人在线| 国产亚洲av高清不卡| 肉色欧美久久久久久久蜜桃| 免费一级毛片在线播放高清视频 | 成年人黄色毛片网站| 国产精品久久久久久精品电影小说| 90打野战视频偷拍视频| 亚洲欧美日韩高清在线视频 | 日韩大码丰满熟妇| 水蜜桃什么品种好| 69av精品久久久久久 | 久久久久久人人人人人| 国内毛片毛片毛片毛片毛片| 精品欧美一区二区三区在线| 曰老女人黄片| 国产又色又爽无遮挡免| 欧美成狂野欧美在线观看| 大片免费播放器 马上看| av在线app专区| 久久中文字幕一级| 91麻豆av在线| 丝袜脚勾引网站| 欧美变态另类bdsm刘玥| 国产伦人伦偷精品视频| 午夜成年电影在线免费观看| 亚洲精品一卡2卡三卡4卡5卡 | 成人手机av| 黄色毛片三级朝国网站| 亚洲少妇的诱惑av| 国产福利在线免费观看视频| av天堂久久9| 在线亚洲精品国产二区图片欧美| 亚洲性夜色夜夜综合| 日本vs欧美在线观看视频| 成年人午夜在线观看视频| 欧美变态另类bdsm刘玥| 操出白浆在线播放| 免费日韩欧美在线观看| 国产av一区二区精品久久| 69精品国产乱码久久久| av电影中文网址| 在线永久观看黄色视频| 别揉我奶头~嗯~啊~动态视频 | 最近中文字幕2019免费版| 18禁裸乳无遮挡动漫免费视频| 又黄又粗又硬又大视频| 日韩电影二区| 超色免费av| 咕卡用的链子| 国产成人系列免费观看| 一级片免费观看大全| xxxhd国产人妻xxx| 99国产极品粉嫩在线观看| 伊人亚洲综合成人网| 91老司机精品| 精品人妻1区二区| 国产精品久久久久久人妻精品电影 | 操出白浆在线播放| 夜夜骑夜夜射夜夜干| 日韩欧美一区视频在线观看| 欧美人与性动交α欧美精品济南到| 欧美激情久久久久久爽电影 | 可以免费在线观看a视频的电影网站| 国产精品影院久久| 日韩一卡2卡3卡4卡2021年| 天天躁夜夜躁狠狠躁躁| 国产av又大| 国产日韩欧美视频二区| 久久精品亚洲av国产电影网| 亚洲精品国产区一区二| 国产熟女午夜一区二区三区| 狠狠精品人妻久久久久久综合| 69精品国产乱码久久久| 999久久久国产精品视频| 国产一区二区三区av在线| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 韩国高清视频一区二区三区| 国产一区二区三区av在线| 天堂俺去俺来也www色官网| 亚洲人成电影观看| 999久久久精品免费观看国产| 日韩大码丰满熟妇| 亚洲精品在线美女| 美女主播在线视频| 成人国产一区最新在线观看| 亚洲av男天堂| 搡老岳熟女国产| 欧美黄色淫秽网站| 国产成人系列免费观看| 亚洲国产欧美在线一区| 国产主播在线观看一区二区| 亚洲欧洲精品一区二区精品久久久| 欧美变态另类bdsm刘玥| 高清av免费在线| 在线观看一区二区三区激情| 色综合欧美亚洲国产小说| 亚洲国产毛片av蜜桃av| 久久久久国产一级毛片高清牌| 亚洲成人免费av在线播放| 极品人妻少妇av视频| 乱人伦中国视频| 日本撒尿小便嘘嘘汇集6| 国产在线一区二区三区精| 国产一区二区在线观看av| 99国产极品粉嫩在线观看| 日本撒尿小便嘘嘘汇集6| 极品少妇高潮喷水抽搐| 精品一区二区三区四区五区乱码| 国产野战对白在线观看| 人妻久久中文字幕网| 日韩视频在线欧美| 国产精品 欧美亚洲| 午夜福利乱码中文字幕| 中文字幕最新亚洲高清| www.999成人在线观看| 淫妇啪啪啪对白视频 | 国产精品一区二区在线不卡| 99re6热这里在线精品视频| xxxhd国产人妻xxx| 人人妻人人澡人人爽人人夜夜| 亚洲中文日韩欧美视频| 91av网站免费观看| 精品福利观看| 色精品久久人妻99蜜桃| 午夜两性在线视频| 久久香蕉激情| 国产精品影院久久| av又黄又爽大尺度在线免费看| 国产精品久久久av美女十八| 午夜久久久在线观看| 少妇被粗大的猛进出69影院| 精品一区二区三区av网在线观看 | 欧美xxⅹ黑人| 狂野欧美激情性bbbbbb| 日本av免费视频播放| 韩国高清视频一区二区三区| 久久精品久久久久久噜噜老黄| 自线自在国产av| 两人在一起打扑克的视频| 欧美av亚洲av综合av国产av| 国产精品 欧美亚洲| 蜜桃在线观看..| 制服人妻中文乱码| www.av在线官网国产| 午夜福利视频精品| 少妇裸体淫交视频免费看高清 | 99国产综合亚洲精品| 亚洲精品国产av成人精品| 男女免费视频国产| 黑人巨大精品欧美一区二区mp4| 中文欧美无线码| 中文字幕制服av| 国产麻豆69| 国产一区有黄有色的免费视频| 精品欧美一区二区三区在线| 久久天堂一区二区三区四区| 一区二区三区四区激情视频| 国产精品自产拍在线观看55亚洲 | 午夜日韩欧美国产| 黄频高清免费视频| 又紧又爽又黄一区二区| 最黄视频免费看| 国产在线一区二区三区精| 久久久久视频综合| h视频一区二区三区| 搡老乐熟女国产| 欧美日本中文国产一区发布| 美国免费a级毛片| 亚洲一码二码三码区别大吗| 一区二区三区激情视频| 午夜老司机福利片| 精品国产一区二区三区四区第35| 亚洲av电影在线观看一区二区三区| 日本精品一区二区三区蜜桃| netflix在线观看网站| 久久精品国产亚洲av高清一级| 亚洲专区国产一区二区| 麻豆国产av国片精品| 精品亚洲乱码少妇综合久久| 欧美精品啪啪一区二区三区 | 国产高清国产精品国产三级| 亚洲精品国产精品久久久不卡| 精品久久久精品久久久| 满18在线观看网站| 激情视频va一区二区三区| 欧美日韩一级在线毛片| 国产成人系列免费观看| 91大片在线观看| 久久久久精品人妻al黑| 欧美乱码精品一区二区三区| 老鸭窝网址在线观看| 建设人人有责人人尽责人人享有的| 日韩有码中文字幕| 国产精品久久久av美女十八| 亚洲国产精品一区三区| 色94色欧美一区二区| 在线十欧美十亚洲十日本专区| 久久久久国内视频| 丝瓜视频免费看黄片| 日韩一卡2卡3卡4卡2021年| 丝瓜视频免费看黄片| 操出白浆在线播放| 在线亚洲精品国产二区图片欧美| 男人操女人黄网站| 美女视频免费永久观看网站| 大香蕉久久成人网| 亚洲熟女毛片儿| 在线观看免费日韩欧美大片| 啦啦啦啦在线视频资源| 如日韩欧美国产精品一区二区三区| 不卡一级毛片| 如日韩欧美国产精品一区二区三区| 97精品久久久久久久久久精品| 国产福利在线免费观看视频| 99国产精品一区二区三区| 99国产综合亚洲精品| 丁香六月欧美| 99精品欧美一区二区三区四区| 国产精品 欧美亚洲| 天天影视国产精品| 99久久精品国产亚洲精品| 久久久久久久国产电影| 中文欧美无线码| 手机成人av网站| 免费高清在线观看日韩| 91精品伊人久久大香线蕉| 亚洲精品美女久久久久99蜜臀| 丰满人妻熟妇乱又伦精品不卡| 精品一区二区三卡| 久久久久视频综合| 在线观看人妻少妇| 久久久久久亚洲精品国产蜜桃av| 男女午夜视频在线观看| 免费不卡黄色视频| 男女午夜视频在线观看| 精品一区二区三区四区五区乱码| 亚洲精品一区蜜桃| 欧美黄色淫秽网站| videos熟女内射| 中文字幕另类日韩欧美亚洲嫩草| 日本wwww免费看| 另类精品久久| 国产成人av教育| 欧美av亚洲av综合av国产av| 亚洲欧美精品自产自拍| 不卡av一区二区三区| 纵有疾风起免费观看全集完整版| 国产成人a∨麻豆精品| 国产av精品麻豆| 欧美 亚洲 国产 日韩一| 国产福利在线免费观看视频| 国产亚洲精品一区二区www | 国产成人精品久久二区二区免费| 90打野战视频偷拍视频| 99九九在线精品视频| 欧美日韩亚洲综合一区二区三区_| 窝窝影院91人妻| 精品人妻熟女毛片av久久网站| 50天的宝宝边吃奶边哭怎么回事| 国产主播在线观看一区二区| 久久久久精品国产欧美久久久 | 嫩草影视91久久| 少妇人妻久久综合中文| 国产亚洲av片在线观看秒播厂| 首页视频小说图片口味搜索| 国产真人三级小视频在线观看| 午夜久久久在线观看| 国产欧美日韩一区二区三 | 悠悠久久av| 一区二区av电影网| 国产野战对白在线观看| 精品久久蜜臀av无| 久久精品aⅴ一区二区三区四区| 欧美精品av麻豆av| 亚洲,欧美精品.| 久久ye,这里只有精品| 欧美老熟妇乱子伦牲交| 高潮久久久久久久久久久不卡| 最近中文字幕2019免费版| 午夜两性在线视频| 国产成人欧美在线观看 | 国产又色又爽无遮挡免| 欧美乱码精品一区二区三区| 国产色视频综合| 国产又爽黄色视频| 丁香六月欧美| 亚洲一区二区三区欧美精品| 国产野战对白在线观看| 我的亚洲天堂| 亚洲精品中文字幕一二三四区 | 欧美黑人欧美精品刺激| 欧美精品啪啪一区二区三区 | 性色av一级| av又黄又爽大尺度在线免费看| 久久久水蜜桃国产精品网| 老司机在亚洲福利影院| 亚洲伊人久久精品综合| 国产91精品成人一区二区三区 | 男人添女人高潮全过程视频| 国产免费视频播放在线视频| 久久久久国内视频| 丝袜喷水一区| 国精品久久久久久国模美| 亚洲精华国产精华精| 乱人伦中国视频| 亚洲国产精品一区三区| 亚洲精品中文字幕一二三四区 | 国产男女超爽视频在线观看| 国产亚洲一区二区精品| 可以免费在线观看a视频的电影网站| 亚洲精品美女久久av网站| 777久久人妻少妇嫩草av网站| 久久久精品国产亚洲av高清涩受| a级片在线免费高清观看视频| 高清欧美精品videossex| 麻豆国产av国片精品| 人妻 亚洲 视频| 一本久久精品| 欧美人与性动交α欧美精品济南到| 日韩欧美一区视频在线观看| 精品久久久精品久久久| 69精品国产乱码久久久| 国产精品久久久人人做人人爽| 欧美日韩国产mv在线观看视频| 日本av免费视频播放| 欧美黑人精品巨大| 女人爽到高潮嗷嗷叫在线视频| 成人黄色视频免费在线看| 人妻 亚洲 视频| 一本大道久久a久久精品| 69精品国产乱码久久久| 热99久久久久精品小说推荐| 亚洲专区字幕在线| 久久香蕉激情| 1024香蕉在线观看| 老司机影院成人| 日本wwww免费看| 亚洲成av片中文字幕在线观看| 精品福利观看| 日韩大片免费观看网站| 国产三级黄色录像| 黄网站色视频无遮挡免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 免费高清在线观看日韩|