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

    木薯及其同源四倍體基因組與DNA甲基化差異分析

    2018-05-14 14:44:47肖亮周慧文謝向譽(yù)尚小紅陳松筆嚴(yán)華兵
    熱帶作物學(xué)報(bào) 2018年5期
    關(guān)鍵詞:DNA甲基化二倍體木薯

    肖亮 周慧文 謝向譽(yù) 尚小紅 陳松筆 嚴(yán)華兵

    摘 要 為研究染色體加倍對木薯基因組的影響,以及多倍化與2種不同倍性木薯之間DNA堿基序列的變化的聯(lián)系。本研究以木薯品種“新選048”二倍體及其同源四倍體為研究材料,使用微衛(wèi)星(simple sequence repeat, SSR)、目標(biāo)起始密碼子多態(tài)性(stand codon targeted polymorphism, SCoT)和甲基化敏感擴(kuò)增多態(tài)性(methylation sensitive amplification polymorphism,MSAP)技術(shù)分別對木薯倍性間的遺傳差異、DNA甲基化水平和模式進(jìn)行分析。結(jié)果表明:利用SCoT和SSR均未檢測出多態(tài)性條帶產(chǎn)生,利用MSAP技術(shù)可檢測出在DNA甲基化水平和模式均發(fā)生較大改變。二倍體和同源四倍體的總甲基化率分別為60.21%、59.52%,全甲基化率分別為39.92%、41.42%,半甲基化率分別為20.29%、18.10%。木薯多倍化后,其中有27.30%的位點(diǎn)發(fā)生了過甲基化,四倍體有25.00%的位點(diǎn)表現(xiàn)出去甲基化。本研究結(jié)果為探究木薯倍性間遺傳差異和表觀遺傳變化規(guī)律進(jìn)行了前期探索。

    關(guān)鍵詞 木薯;二倍體;同源四倍體;DNA甲基化;遺傳差異

    中圖分類號(hào) S533 文獻(xiàn)標(biāo)識(shí)碼 A

    Abstract To study the effect of chromosome doubling in cassava genome and, explore the relationship between changes of DNA sequence in two different ploidy levels and polyploidization variation in cassava, the cassava variety ‘Xinxuan 048 and its autotetraploid were taken as the materials to analyze the genomic mutation and DNA methylation diversity using microsatellite (Simple Sequence Repeat, SSR) and Stand codon targeted polymorphism (SCoT) molecular marker and Methylation Sensitive Amplification Polymorphism (MSAP) technique. The results showed that no polymorphic band occurred using SCoT and SSR molecular marker. In the level of DNA methylation,both of the level and pattern of the autotetraploid cassava changed. The total methylation rate of the diploid and the autotetraploid was 60.21% and 59.52%, respectively;the full methylation rate was 39.92% and 41.42%, respectively;the hemi methylation rate was 20.29% and 18.10%, respectively. In comparison to the diploid, the pattern of the methylation showed that about 25.00% sites was demethylation, 30.63% sites has been methylated in autotetraploid cassava. This research is a foundation for the further exploration of genetic differences and epigenetic regulation after cassava polyploidization.

    Key words cassava; diploid; autotetraploid; DNA methylation; genetic differences

    doi 10.3969/j.issn.1000-2561.2018.05.001

    木薯(Manihot esculenta Crantz)(2n=36)是大戟科(Euphorbiaceae)木薯屬(Manihot)植物,是熱帶地區(qū)第三大糧食作物,全球第六大糧食作物,被譽(yù)為“淀粉之王”,是世界近6億人賴以生存的糧食。除此之外,木薯還可飼用和加工成各種工業(yè)產(chǎn)品,如淀粉、酒精等[1]。因此,高產(chǎn)穩(wěn)產(chǎn)和特色優(yōu)質(zhì)木薯品種選育顯得尤為必要[2]。

    染色體多倍化在植物進(jìn)化過程中頻繁發(fā)生,對于產(chǎn)生新物種貢獻(xiàn)很大。同源四倍體植株能在細(xì)胞學(xué)、生理生化及表型方面發(fā)生顯著變化[3]。木薯為無性繁殖作物,基因組高度雜合,人工誘導(dǎo)多倍體是木薯種質(zhì)創(chuàng)新的理想途徑[4]。伴隨著同源多倍化,植物在基因組和表觀遺傳等水平都發(fā)生了顯著變化,如染色體重組、DNA序列消除、DNA甲基化等[5]。

    分子標(biāo)記技術(shù)是目前研究上述領(lǐng)域的有效方法之一。梁武軍等[6]和鐘鳴等[7]利用簡單重復(fù)序列(simple sequence repeat, SSR)鑒定葡萄柚(Citrus paradisi Macf)后代和柳枝稷(Panicum virgatum)的遺傳多樣性。韓國輝等[8]和周慧文等[9]分別利用目標(biāo)起始密碼子多態(tài)性分子標(biāo)記(start condon targeted polymorphism,SCoT)分析枇杷(Eriobotrya japonica (Thunb.) Lindl.)和木薯的遺傳多樣性。然而,目前利用SSR和SCoT技術(shù)對不同倍性木薯間遺傳差異的相關(guān)研究卻鮮見報(bào)道。

    植物中甲基化位點(diǎn)包括3種:CG、CHH和CHG(H表示除G以外的核苷酸)[10]。李雅婷等[11]曾利用甲基化敏感擴(kuò)增多態(tài)性(methylation sensitive amplification polymorphism,MSAP)技術(shù)研究發(fā)現(xiàn)鐵皮石斛(Dendrobium officinale Kimura et Migo)同源四倍體有28.63%的位點(diǎn)發(fā)生了去甲基化,30.63%的位點(diǎn)發(fā)生了過甲基化。王紅娟等[12]以20對引物研究了二倍體與同源四倍體茅蒼術(shù)(Rhizoma Areactylodis Lanceae)的甲基化差異的擴(kuò)增總甲基化率為56.92%和55.03%,四倍體基因組DNA半甲基化率高于二倍體,而全甲基化率低于二倍體;四倍體與二倍體相比,DNA甲基化模式發(fā)生調(diào)整的基因位點(diǎn)占總檢測位點(diǎn)的52.53%。

    本研究利用SSR和SCoT技術(shù)分析木薯品種“新選048”二倍體及其同源四倍體間的遺傳差異,并利用MSAP技術(shù)探究木薯二倍體與同源四倍體的DNA甲基化水平與模式的變化,旨在探究木薯倍性間DNA堿基序列變化與木薯同源多倍化之間的聯(lián)系以及為揭示木薯不同倍性間甲基化變異趨勢奠定理論基礎(chǔ)。

    1 材料與方法

    1.1 材料

    本研究在廣西大學(xué)亞熱帶農(nóng)業(yè)生物資源保護(hù)與利用國家重點(diǎn)實(shí)驗(yàn)室進(jìn)行,供試的木薯品種“新選048”二倍體親本由廣西農(nóng)業(yè)科學(xué)院經(jīng)濟(jì)作物研究所提供。其同源四倍體由本課題組前期使用秋水仙素加倍該品種獲得[13]。

    19條SCoT引物序列參照嚴(yán)華兵等[14]文中所列(表1),36對SSR引物參照王明等[15]文中所列(表2)。MSAP接頭序列與引物序列見表3。

    1.2 方法

    1.2.1 材料獲得 將二倍體親本及其同源四倍體的6個(gè)不同株系的組培苗種植于溫室大棚內(nèi),取葉片用于分子標(biāo)記分析;在同源四倍體中挑選一個(gè)株系,待莖稈成熟后,將莖稈砍成20 cm一段的稈扦插至盆中,2個(gè)月后取二倍體和四倍體的頂端嫩葉,用于甲基化分析。

    1.2.2 木薯全基因組DNA提取 使用常規(guī)CTAB法提取其DNA。將提取好的DNA稀釋成60 ng/μL濃度,放于-20 ℃冰箱保存。

    1.2.3 SSR標(biāo)記擴(kuò)增與檢測 PCR擴(kuò)增體系與程序參照王明等[15]。SSR擴(kuò)增產(chǎn)物均用7%變性聚丙烯酰胺凝膠電泳進(jìn)行檢測,以恒定電壓600 V跑膠20 min。電泳結(jié)束后用1 g/L AgNO3染色8~10 min,顯影后拍照保存。

    1.2.4 SCoT標(biāo)記擴(kuò)增與檢測 PCR擴(kuò)增體系與程序參照相關(guān)文獻(xiàn)[14]。擴(kuò)增反應(yīng)結(jié)束后,取5 μL擴(kuò)增產(chǎn)物,經(jīng)1.5%瓊脂糖凝膠電泳、GoldView核酸染色劑染色后在紫外凝膠成像系統(tǒng)上拍照保存。

    1.2.5 MSAP分析 酶切體系完全參照李雅婷等[11]的方法。酶切完畢后,在酶切片段的兩端加上人工設(shè)計(jì)的與EcoR I和Hpa II/Msp I酶切位點(diǎn)互補(bǔ)的人工接頭(表3)。接頭連接體系為20 μL,其中包含:2 μL 10×T4 Buffer、0.4 μL 20 μmol/L EcoR I和Hpa II/Msp I 接頭,其余用水補(bǔ)齊,16 ℃過夜。隨后進(jìn)行的預(yù)擴(kuò)增和選擇性擴(kuò)增使用的體系和PCR程序同樣參照李雅婷等[11]的方法。擴(kuò)增產(chǎn)物經(jīng)94 ℃變性10 min后于6%變性聚丙烯酰胺凝膠進(jìn)行垂直電泳,銀染后觀察結(jié)果。每份材料擁有2條泳道,為同裂酶HpaⅡ、MspⅠ與EcoRⅠ組合進(jìn)行雙酶切所產(chǎn)生。

    參照王春國等[16]的方法,將木薯二倍體與其同源四倍體的DNA甲基化水平變化分為4大類共15個(gè)亞類。A類是指在同一位點(diǎn),二倍體與四倍體的擴(kuò)增模式未發(fā)生改變,即不存在擴(kuò)增多態(tài)性;B類表明在同一位點(diǎn)四倍體表現(xiàn)出過甲基化或超甲基化,即四倍體的甲基化水平較二倍體上調(diào);C類表明四倍體表現(xiàn)出去甲基化,即四倍體較二倍體的甲基化水平下調(diào);D類表明MSAP技術(shù)無法確定該位點(diǎn)甲基化狀態(tài)調(diào)整,但該位點(diǎn)的甲基化模式同樣產(chǎn)生多態(tài)性。

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

    在泳道的同一位點(diǎn)處,有條帶記為“1”,無條帶記為“0”,篩選條帶數(shù)量豐富,擴(kuò)增清晰的引物進(jìn)行條帶數(shù)據(jù)統(tǒng)計(jì)分析。MSAP條帶統(tǒng)計(jì)分析方法參考王春國等[16]和邊禹[17]的方法。

    2 結(jié)果與分析

    2.1 木薯二倍體與同源四倍體基因組的遺傳差異

    用SSR標(biāo)記的36對引物對木薯二倍體及其同源四倍體共7份材料進(jìn)行擴(kuò)增,部分?jǐn)U增結(jié)果如圖1所示。從圖2中可見,所使用的引物可以在這7份材料中擴(kuò)增出清晰的條帶。每對引物擴(kuò)增的條帶數(shù)在1~5之間,在7份材料中一共擴(kuò)增出672條條帶,平均每對引物在每份材料中擴(kuò)增出2.7條條帶,條帶大小主要集中于250~750 bp之間,但是沒有一條引物可以檢測到有新增位點(diǎn)或者缺失位點(diǎn)。

    用19條SCoT分子標(biāo)記引物對木薯二倍體及其同源四倍體共7份材料進(jìn)行擴(kuò)增,部分?jǐn)U增結(jié)果如圖2所示。從圖2中可見,所使用的引物均能在這7份材料中擴(kuò)增出清晰的條帶。每條引物擴(kuò)增條帶數(shù)在3~12條之間,大多數(shù)條帶大小集中在500~3 000 bp之間。19條SCoT引物共擴(kuò)增出1 008條條帶,平均每條引物在每份材料中擴(kuò)增出7.6條條帶。但結(jié)果同使用SSR分子標(biāo)記分析倍性間遺傳差異相同,均未發(fā)現(xiàn)有多態(tài)性條帶產(chǎn)生。雖然條帶有亮有暗,推測可能是由于基因的劑量效應(yīng)導(dǎo)致,而與堿基序列變化無關(guān)。所使用的6份同源四倍體材料與其二倍體親本沒有檢測出新增或缺失位點(diǎn),即利用SCoT分子標(biāo)記也未能檢測出木薯倍性間的遺傳差異。

    2.2 木薯二倍體與四倍體的MSAP結(jié)果

    從50個(gè)引物組合中篩選出了23個(gè)擴(kuò)增條帶豐富,清晰可辨的引物組合,以這23個(gè)引物組合所擴(kuò)增出的位點(diǎn)進(jìn)行分析,結(jié)果如表4所示。23個(gè)引物組合共擴(kuò)增出1 500個(gè)位點(diǎn),其中,二倍體共擴(kuò)增出754個(gè)位點(diǎn),同源四倍體共擴(kuò)增出746個(gè)位點(diǎn)。雖然倍性間的總甲基化率幾近持平,四倍體略低于二倍體0.69%,但是甲基化類型發(fā)生了變化。同源四倍體雙鏈內(nèi)甲基化的位點(diǎn)低于二倍體4.84%,單鏈外甲基化的位點(diǎn)低于二倍體2.19%,完全甲基化的位點(diǎn)高于二倍體6.34%。其中,四倍體半甲基化率低于二倍體2.19%,全甲基化率高于二倍體1.5%。由此可見,木薯二倍體與其同源四倍體的總甲基化水平相差較小,但是各個(gè)甲基化類型發(fā)生了一定規(guī)模的改變。從DNA甲基化水平類型看(表5),木薯二倍體及其同源四倍體擴(kuò)增類型為A類的位點(diǎn)占總擴(kuò)增位點(diǎn)的43.62%,B類擴(kuò)增位點(diǎn)數(shù)占總數(shù)的27.30%。C類的擴(kuò)增位點(diǎn)數(shù)的比例是25.00%,D類擴(kuò)增條帶比例最小,為4.08%。綜上所述,多態(tài)性擴(kuò)增位點(diǎn)數(shù)目占總擴(kuò)增數(shù)目的56.38%,即木薯多倍化后,共有56.38%的位點(diǎn)甲基化模式發(fā)生了調(diào)整。

    H2和M2分別代表木薯二倍體EcoR I/Hpa II和EcoR I/Msp I酶切擴(kuò)增結(jié)果。H4和M4代表木薯四倍體EcoR I/Hpa II和EcoR I/Msp I酶切擴(kuò)增結(jié)果。箭頭和字母表示不同甲基化擴(kuò)增模式(與表4對應(yīng))。E39Msp59-E44Msp39為不同引物組合。M為marker。

    3 討論

    木薯為無性繁殖作物,基因組高度雜合,這決定了多倍體育種可作為木薯種質(zhì)創(chuàng)新的重要途徑[4]。2017年,有報(bào)道稱為了滿足人們生活需要,木薯長期的無性繁殖固定了大量的有害突變,同源多倍化產(chǎn)生冗余基因?qū)е嘛@性基因掩蓋隱性等位基因的有害突變,大大減少了純合隱性突變,降低有害個(gè)體發(fā)生的頻率;另一方面,同源多倍化還可以通過改變生物體重要基因拷貝數(shù)的增加會(huì)導(dǎo)致基因功能的多樣化[18]。本研究中利用SCoT和SSR分子標(biāo)記分析木薯“新選048”二倍體及同源四倍體的DNA堿基序列進(jìn)行檢測,并未發(fā)現(xiàn)木薯二倍體與其同源四倍體間的遺傳差異,這與曾少華[19]對柑橘國慶1號(hào)的研究結(jié)論一致,原因可能為:(1)選取標(biāo)記數(shù)量太少,如果在全基因組范圍內(nèi)多挑選一些SSR或許二倍體和四倍體之間能表現(xiàn)出遺傳差異。(2)本研究使用的檢測手段——SSR分子標(biāo)記和SCoT分子標(biāo)記具有局限性,未能對木薯二倍體與同源四倍體進(jìn)行全面檢測。或許更換分子標(biāo)記種類,如使用AFLP或許可以檢測出遺傳差異。

    與通過高通量測序的重亞硫酸鹽處理DNA來檢測全基因組甲基化信息的方法相比,在AFLP分子標(biāo)記技術(shù)上改進(jìn)的甲基化敏感限制性酶切法(MSAP法)更為經(jīng)濟(jì)實(shí)惠,且該方法操作簡便,可在全基因組水平上進(jìn)行DNA甲基化模式分析[20],目前已經(jīng)被廣泛應(yīng)用于分析植物脅迫處理[21-22]、輻射誘變[23]以及化學(xué)方法誘變的DNA甲基化變異中。在本研究中應(yīng)用MSAP技術(shù)分析木薯二倍體與其同源四倍體的甲基化變異時(shí)發(fā)現(xiàn),兩個(gè)倍性的甲基化水平十分相近,僅相差0.69%,但是有超過一半的位點(diǎn)的甲基化狀態(tài)發(fā)生了改變,這與王春國等[16]結(jié)果相似。

    前人研究發(fā)現(xiàn),去甲基化可以引起轉(zhuǎn)基因擬南芥的表型產(chǎn)生巨大變異[24],本課題組觀察到經(jīng)秋水仙素誘導(dǎo)產(chǎn)生的同源四倍體木薯品種“新選048”有葉片變寬圓,葉色深綠,植株變矮和節(jié)間距變短現(xiàn)象[13]。根據(jù)MSAP甲基化分析結(jié)果推測,產(chǎn)生表型變異的原因可能是由于染色體加倍導(dǎo)致基因的“劑量效應(yīng)”或是DNA甲基化水平和模式改變而引起。

    本研究利用SSR與SCoT分子標(biāo)記技術(shù)未能檢測出木薯品種“新選048”二倍體與其同源四倍體之間的基因組差異。利用MSAP技術(shù)可檢測出木薯品種“新選048”二倍體與其同源四倍體之間的DNA甲基化水平和模式發(fā)生了變化;其中,二倍體的總甲基化水平為60.21%;四倍體的總甲基化水平為59.52%,低于二倍體0.69%;木薯基因組加倍后,有56.38%的位點(diǎn)的甲基化模式發(fā)生了變異,其中有27.30%的位點(diǎn)發(fā)生了過甲基化,四倍體有25.00%的位點(diǎn)表現(xiàn)出去甲基化,此結(jié)論與鐵皮石斛[11],茅蒼術(shù)[12],水稻[25],小麥[26]等研究結(jié)論基本一致。本研究為揭示木薯多倍化后表觀遺傳調(diào)控提供了理論依據(jù)。

    參考文獻(xiàn)

    [1] Balagopalan C. Cassava utilization in food, feed and industry [C]// Hillocks R J,Thresh J M,Bellotti A C (eds). Cassava: biology, production and utilization. CAB International,2002: 301-318.

    [2] 嚴(yán)華兵,葉劍秋,李開綿. 中國木薯育種研究進(jìn)展[J]. 中國農(nóng)學(xué)通報(bào),2015,31(15):63-70.

    [3] Abdoli M,Moieni A,Naghdi Badi H. Morphological,physiological,cytological and phytochemical studies in diploid and colchicine-induced tetraploid plants of Echinacea purpurea(L.)[J]. Acta Physiologiae Plantarum,2013,35(7):2 075-2 083.

    [4] Yan H ,Lu L,Hershey C,et al. Cassava mutation breeding:current status and trends[J]. Plant Mutation Reports,2013,33(11):643-645.

    [5] Doyle J J,F(xiàn)lagel L E,Paterson A H,et al. Evolutionary genetics of genome merger and doubling in plants[J]. Annual Review of Genetics,2008,42:443-461.

    [6] 梁武軍,解凱東,謝宗周,等. 三倍體葡萄柚實(shí)生后代多倍體的發(fā)掘與SSR遺傳鑒定[J]. 果樹學(xué)報(bào),2015,32(1):13-18.

    [7] 鐘 鳴,周思凡,張新全,等. 利用EST-SSR、IT-ISJ分子標(biāo)記研究四倍體與八倍體柳枝稷的遺傳多樣性[J]. 草業(yè)學(xué)報(bào),2016,25(10):113-123.

    [8] 韓國輝,汪衛(wèi)星,向素瓊,等. 不同倍性枇杷SCoT標(biāo)記技術(shù)體系的建立與應(yīng)用[C]//全國枇杷學(xué)術(shù)研討會(huì)論文(摘要)集. 2011.

    [9] 周慧文,單建偉,馮 斗,等. SCoT分子標(biāo)記技術(shù)在木薯遺傳多樣性分析中的應(yīng)用[J]. 熱帶作物學(xué)報(bào),2015,36(8):1 440-1 444.

    [10] Cokus S J,F(xiàn)eng S,Zhang X,et al. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning[J]. Nature,2008,452(7184):215-219.

    [11] 李雅婷,王紅娟,向增旭. 鐵皮石斛同源四倍體與二倍體基因組及甲基化差異分析[J]. 核農(nóng)學(xué)報(bào),2015,29(10):1 901-1 908.

    [12] 王紅娟,巢建國,李雅婷,等. 二倍體與同源四倍體茅蒼術(shù)基因組DNA甲基化水平與模式的MSAP分析[J]. 核農(nóng)學(xué)報(bào),2015,29(8):1 502-1 508.

    [13] Zhou H W,Zeng W D,Yan H B. In vitro induction of tetraploids in cassava variety ‘Xinxuan 048 using colchicine[J]. Plant Cell,Tissue and Organ Culture,2016,128(3):723-729.

    [14] 嚴(yán)華兵,周慧文,單建偉,等. 木薯SCoT-PCR反應(yīng)體系優(yōu)化及引物篩選[J]. 南方農(nóng)業(yè)學(xué)報(bào),2016,47(10):1 648-1 652.

    [15] 王 明,肖鑫輝,安飛飛,等. 利用SSR標(biāo)記分析木薯遺傳多樣性[J]. 熱帶農(nóng)業(yè)科學(xué),2015,35(11):38-44.

    [16] 王春國,古 瑜,陳成彬,等. 不同倍性西瓜基因組DNA甲基化水平與模式的MSAP分析[J]. 分子細(xì)胞生物學(xué)報(bào),2009(2):118-126.

    [17] 邊 禹. 不同倍性多倍體枇杷基因組變異及其DNA甲基化分析[D]. 重慶:西南大學(xué),2010.

    [18] Ramu P,Esuma W,Kawuki R,et al. Cassava haplotype map highlights fixation of deleterious mutations during clonal propagation[J]. Nature Genetics,2017,49(6):959-963.

    [19] 曾少華. 離體加倍創(chuàng)造柑橘多倍體種質(zhì)及其細(xì)胞學(xué)和分子鑒定[D]. 武漢:華中農(nóng)業(yè)大學(xué),2005.

    [20] 馬浪浪,江 舟,黃小波,等. 植物DNA甲基化調(diào)控研究進(jìn)展[J]. 中國生物工程雜志,2013,33(9):101-110.

    [21] 孫麗芳,胡凱鳳,高樹仁,等. 鹽脅迫誘導(dǎo)玉米苗期DNA甲基化變異的研究[J]. 核農(nóng)學(xué)報(bào),2017,31(4):627-634.

    [22] 朱紅菊,劉文革,趙勝杰,等. NaCl脅迫下二倍體和同源四倍體西瓜幼苗DNA甲基化差異分析[J]. 中國農(nóng)業(yè)科學(xué),2014,47(20):4 045-4 055.

    [23] 趙 倩,王 巍,張 萌,等. 重離子輻射誘導(dǎo)水稻DNA甲基化變化的分析[J]. 核農(nóng)學(xué)報(bào),2016,30(9):1 665-1 671.

    [24] Ronemus M J,Galbiati M,Ticknor C,et al. Demethylation-induced developmental pleiotropy in Arabidopsis[J]. Science,1996,273(5 275):654-657.

    [25] 楊 瓊,李 輝,左欽月,等. 同源四倍體與二倍體水稻基因組及甲基化差異分析[J]. 應(yīng)用與環(huán)境生物學(xué)報(bào),2013,19(4):630-636.

    [26] 聶利紅,王延召,孫其信,等. 人工異源六倍體小麥中甲基化差異的MSAP分析[J]. 華北農(nóng)學(xué)報(bào),2012,27(2):77-80.

    猜你喜歡
    DNA甲基化二倍體木薯
    刮木薯
    柬埔寨拜靈木薯喜獲大豐收,市場價(jià)格保持穩(wěn)定
    染色體加倍對黃毛草莓葉片形態(tài)的影響
    二倍體太子參叢生芽誘導(dǎo)研究
    挖木薯
    胃癌中DNA甲基化對microRNA調(diào)控作用的研究進(jìn)展
    阿爾茨海默病小鼠模型海馬組織AtP5a1基因甲基化改變
    DNA甲基化對植物生長發(fā)育的調(diào)控研究
    “寒富”蘋果與其同源四倍體耐鹽差異研究
    MSRE-qPCR法檢測少弱精子癥患者父源印記基因H19上游區(qū)域甲基化水平
    国产在线观看jvid| 啪啪无遮挡十八禁网站| 夜夜看夜夜爽夜夜摸 | 精品电影一区二区在线| 国产精品秋霞免费鲁丝片| 久久久国产精品麻豆| 超碰97精品在线观看| 久久 成人 亚洲| 亚洲av第一区精品v没综合| 美女高潮喷水抽搐中文字幕| 亚洲熟女毛片儿| 国产精品98久久久久久宅男小说| 99久久综合精品五月天人人| а√天堂www在线а√下载| 在线观看免费视频网站a站| 韩国av一区二区三区四区| 欧美丝袜亚洲另类 | 超碰成人久久| 99久久99久久久精品蜜桃| 99国产精品一区二区三区| 欧美日韩精品网址| 午夜两性在线视频| 在线观看66精品国产| 黄网站色视频无遮挡免费观看| 大型黄色视频在线免费观看| 亚洲少妇的诱惑av| 久久人人97超碰香蕉20202| 91字幕亚洲| 巨乳人妻的诱惑在线观看| 精品久久久精品久久久| 波多野结衣av一区二区av| 亚洲人成77777在线视频| 99香蕉大伊视频| 丰满迷人的少妇在线观看| 水蜜桃什么品种好| 久久亚洲精品不卡| 又大又爽又粗| www.精华液| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 天堂中文最新版在线下载| x7x7x7水蜜桃| 无限看片的www在线观看| 免费看a级黄色片| 国产精品爽爽va在线观看网站 | 国产蜜桃级精品一区二区三区| 亚洲中文av在线| 51午夜福利影视在线观看| 1024视频免费在线观看| 老汉色av国产亚洲站长工具| 久久久久久久午夜电影 | 美女福利国产在线| 亚洲一区二区三区不卡视频| 欧美成人免费av一区二区三区| 超色免费av| 美女高潮到喷水免费观看| 久久国产乱子伦精品免费另类| 变态另类成人亚洲欧美熟女 | 多毛熟女@视频| 美女 人体艺术 gogo| 一级毛片高清免费大全| av网站在线播放免费| 国产高清videossex| 人妻丰满熟妇av一区二区三区| 成人精品一区二区免费| 欧美日本亚洲视频在线播放| 俄罗斯特黄特色一大片| 国产av又大| 高清毛片免费观看视频网站 | 日韩视频一区二区在线观看| 久久亚洲真实| 99久久精品国产亚洲精品| 夜夜躁狠狠躁天天躁| 亚洲欧美精品综合一区二区三区| 乱人伦中国视频| 他把我摸到了高潮在线观看| 亚洲一码二码三码区别大吗| 精品熟女少妇八av免费久了| 欧美成人免费av一区二区三区| 亚洲午夜理论影院| 大型黄色视频在线免费观看| 国产精品二区激情视频| 亚洲黑人精品在线| 两性夫妻黄色片| 法律面前人人平等表现在哪些方面| 国产熟女xx| 黄色 视频免费看| 少妇被粗大的猛进出69影院| 亚洲中文av在线| av网站免费在线观看视频| 操出白浆在线播放| 久久久国产成人精品二区 | 免费看a级黄色片| 久久狼人影院| 国产一区二区在线av高清观看| 国产精华一区二区三区| 精品久久蜜臀av无| 美国免费a级毛片| 亚洲国产精品sss在线观看 | 免费日韩欧美在线观看| 亚洲精品一二三| 欧美乱妇无乱码| 一级毛片女人18水好多| 成人三级黄色视频| 在线观看免费日韩欧美大片| 国产亚洲精品综合一区在线观看 | 国产欧美日韩一区二区三区在线| 亚洲精品一区av在线观看| 欧美乱色亚洲激情| 18禁黄网站禁片午夜丰满| 69av精品久久久久久| 国产一区二区三区在线臀色熟女 | 女性生殖器流出的白浆| 中文字幕人妻丝袜一区二区| 一边摸一边做爽爽视频免费| 久9热在线精品视频| 真人一进一出gif抽搐免费| 欧美成狂野欧美在线观看| 高清毛片免费观看视频网站 | 免费高清在线观看日韩| 天堂动漫精品| 国产精品香港三级国产av潘金莲| 巨乳人妻的诱惑在线观看| 午夜福利一区二区在线看| 一二三四社区在线视频社区8| 美女高潮喷水抽搐中文字幕| 久久久国产精品麻豆| 欧美成狂野欧美在线观看| 操美女的视频在线观看| 午夜91福利影院| 一进一出抽搐gif免费好疼 | av天堂久久9| 欧美国产精品va在线观看不卡| 亚洲专区国产一区二区| 欧美一区二区精品小视频在线| 亚洲片人在线观看| 国产成年人精品一区二区 | 香蕉久久夜色| 每晚都被弄得嗷嗷叫到高潮| 久久久久久大精品| 日本三级黄在线观看| 欧美成人午夜精品| 日韩三级视频一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 水蜜桃什么品种好| 99国产精品免费福利视频| 精品无人区乱码1区二区| 国产精品久久久久成人av| 在线观看一区二区三区激情| 麻豆av在线久日| 黄色a级毛片大全视频| 日韩免费av在线播放| av在线天堂中文字幕 | 亚洲自拍偷在线| 精品福利永久在线观看| 91老司机精品| 涩涩av久久男人的天堂| 中文欧美无线码| 午夜日韩欧美国产| 久久精品国产清高在天天线| 亚洲一卡2卡3卡4卡5卡精品中文| 好看av亚洲va欧美ⅴa在| 男女下面进入的视频免费午夜 | 99国产精品免费福利视频| 成人手机av| 亚洲欧美日韩高清在线视频| 超碰97精品在线观看| 国产一区二区三区综合在线观看| 熟女少妇亚洲综合色aaa.| 亚洲人成网站在线播放欧美日韩| 天天添夜夜摸| 法律面前人人平等表现在哪些方面| 国产99白浆流出| 露出奶头的视频| 熟女少妇亚洲综合色aaa.| 国产99白浆流出| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 亚洲欧美精品综合一区二区三区| 国产精品自产拍在线观看55亚洲| 91大片在线观看| 91字幕亚洲| 大型av网站在线播放| 两个人看的免费小视频| 丁香欧美五月| 久99久视频精品免费| 老司机亚洲免费影院| 久久中文看片网| 如日韩欧美国产精品一区二区三区| 欧美乱码精品一区二区三区| 午夜福利,免费看| 性少妇av在线| 水蜜桃什么品种好| 99精国产麻豆久久婷婷| 最近最新中文字幕大全电影3 | 精品卡一卡二卡四卡免费| 欧美在线黄色| 三级毛片av免费| 国产乱人伦免费视频| 在线天堂中文资源库| 精品久久蜜臀av无| 嫩草影视91久久| 亚洲av熟女| 免费在线观看日本一区| 美女大奶头视频| 黄片大片在线免费观看| 亚洲一区中文字幕在线| 99在线人妻在线中文字幕| 欧美中文日本在线观看视频| 亚洲专区中文字幕在线| 高清毛片免费观看视频网站 | 亚洲精品粉嫩美女一区| 多毛熟女@视频| 亚洲 国产 在线| 免费搜索国产男女视频| 搡老熟女国产l中国老女人| 丰满的人妻完整版| 女人被狂操c到高潮| 亚洲中文字幕日韩| 俄罗斯特黄特色一大片| 国产精品野战在线观看 | 看免费av毛片| 一a级毛片在线观看| 亚洲人成电影免费在线| 久久香蕉精品热| 大陆偷拍与自拍| 正在播放国产对白刺激| 精品国产一区二区三区四区第35| 国产成人一区二区三区免费视频网站| 色老头精品视频在线观看| 欧美av亚洲av综合av国产av| 久久久久久人人人人人| 国产熟女xx| 欧美乱色亚洲激情| 老熟妇仑乱视频hdxx| 91成年电影在线观看| 老司机午夜十八禁免费视频| 久久影院123| 一边摸一边抽搐一进一出视频| 涩涩av久久男人的天堂| 日本黄色日本黄色录像| 欧美乱妇无乱码| 精品久久久久久久久久免费视频 | av中文乱码字幕在线| 久久久久久大精品| 男人操女人黄网站| 丰满迷人的少妇在线观看| 精品福利观看| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩中文字幕国产精品一区二区三区 | 久久久精品国产亚洲av高清涩受| 在线观看免费高清a一片| 午夜精品国产一区二区电影| 桃红色精品国产亚洲av| 777久久人妻少妇嫩草av网站| 一个人免费在线观看的高清视频| 脱女人内裤的视频| 黄片小视频在线播放| 首页视频小说图片口味搜索| 亚洲精品在线观看二区| 国产欧美日韩综合在线一区二区| 精品一区二区三区视频在线观看免费 | 女警被强在线播放| 欧美亚洲日本最大视频资源| www.www免费av| 超碰成人久久| 叶爱在线成人免费视频播放| 看黄色毛片网站| 1024香蕉在线观看| 成人18禁在线播放| 在线观看免费视频网站a站| 免费在线观看亚洲国产| 三上悠亚av全集在线观看| 亚洲精华国产精华精| 亚洲专区字幕在线| av在线天堂中文字幕 | 一进一出好大好爽视频| 美女扒开内裤让男人捅视频| 最好的美女福利视频网| 欧美人与性动交α欧美精品济南到| 一边摸一边抽搐一进一出视频| 久久久精品欧美日韩精品| 久久人妻av系列| 久久国产精品男人的天堂亚洲| 日本一区二区免费在线视频| 老司机靠b影院| a级毛片在线看网站| 亚洲一区二区三区不卡视频| 嫩草影院精品99| 精品久久久久久久久久免费视频 | 国产成人欧美在线观看| 亚洲全国av大片| 国产成人影院久久av| 99国产精品一区二区三区| 五月开心婷婷网| 丰满迷人的少妇在线观看| 99re在线观看精品视频| 在线av久久热| 国产精品影院久久| 久久午夜亚洲精品久久| 日韩免费高清中文字幕av| 国产一区在线观看成人免费| 国产精华一区二区三区| 91九色精品人成在线观看| 成人三级黄色视频| 国产色视频综合| 日日爽夜夜爽网站| 国产一区在线观看成人免费| 国产欧美日韩一区二区精品| 不卡av一区二区三区| 一边摸一边做爽爽视频免费| 亚洲精品美女久久av网站| av片东京热男人的天堂| 麻豆成人av在线观看| 亚洲精品粉嫩美女一区| 精品久久久久久成人av| 999久久久国产精品视频| av网站免费在线观看视频| 国产精品乱码一区二三区的特点 | 久久久国产成人精品二区 | 国产av一区在线观看免费| 巨乳人妻的诱惑在线观看| 久久久久久亚洲精品国产蜜桃av| 最好的美女福利视频网| 成人国语在线视频| 久久久国产欧美日韩av| 国产单亲对白刺激| 久久国产精品男人的天堂亚洲| 国产亚洲欧美98| 亚洲 国产 在线| 99香蕉大伊视频| 在线观看www视频免费| 中国美女看黄片| 成年版毛片免费区| 久久亚洲真实| 久久国产精品影院| 免费看十八禁软件| 亚洲欧美激情在线| 国产色视频综合| 夜夜看夜夜爽夜夜摸 | 欧美日韩瑟瑟在线播放| cao死你这个sao货| 亚洲色图av天堂| 欧美成人性av电影在线观看| 亚洲av五月六月丁香网| av网站免费在线观看视频| 亚洲国产精品合色在线| 母亲3免费完整高清在线观看| 日韩精品中文字幕看吧| 99国产精品一区二区蜜桃av| 两性午夜刺激爽爽歪歪视频在线观看 | 国产激情久久老熟女| 精品一区二区三卡| 精品久久久久久成人av| 国产精品久久久人人做人人爽| 国产成人啪精品午夜网站| av欧美777| www.精华液| av欧美777| 国产精品久久电影中文字幕| 精品国产美女av久久久久小说| 国产精品 欧美亚洲| 国产精品香港三级国产av潘金莲| svipshipincom国产片| 久久狼人影院| 国产片内射在线| 一进一出好大好爽视频| 精品一区二区三区av网在线观看| 免费高清视频大片| 国产激情欧美一区二区| videosex国产| 搡老熟女国产l中国老女人| 黄色视频,在线免费观看| 91在线观看av| а√天堂www在线а√下载| 欧美日韩瑟瑟在线播放| 午夜福利在线观看吧| 侵犯人妻中文字幕一二三四区| 老鸭窝网址在线观看| 欧美日韩av久久| 高清黄色对白视频在线免费看| 亚洲国产精品合色在线| 色精品久久人妻99蜜桃| 午夜两性在线视频| 如日韩欧美国产精品一区二区三区| 身体一侧抽搐| 夜夜爽天天搞| 男男h啪啪无遮挡| 精品人妻在线不人妻| 国产极品粉嫩免费观看在线| 97人妻天天添夜夜摸| 国产一区二区在线av高清观看| 欧美另类亚洲清纯唯美| 三上悠亚av全集在线观看| 久久精品91无色码中文字幕| netflix在线观看网站| 美女高潮喷水抽搐中文字幕| 成人三级做爰电影| 无遮挡黄片免费观看| 国产精品成人在线| 日韩精品中文字幕看吧| 91麻豆精品激情在线观看国产 | 两人在一起打扑克的视频| 国产精品乱码一区二三区的特点 | 高清黄色对白视频在线免费看| 国产精品野战在线观看 | 最近最新中文字幕大全免费视频| 狂野欧美激情性xxxx| 黄色 视频免费看| 国产精品 国内视频| 黄色a级毛片大全视频| 精品国产超薄肉色丝袜足j| 免费观看人在逋| 国产高清国产精品国产三级| 午夜两性在线视频| 亚洲精品美女久久久久99蜜臀| 亚洲熟妇中文字幕五十中出 | 一级黄色大片毛片| 国产精品久久久久久人妻精品电影| 午夜福利影视在线免费观看| 女生性感内裤真人,穿戴方法视频| 亚洲男人的天堂狠狠| 涩涩av久久男人的天堂| 国产精品美女特级片免费视频播放器 | 黄频高清免费视频| 欧美日韩国产mv在线观看视频| 中文字幕人妻丝袜一区二区| 午夜老司机福利片| 久久人人97超碰香蕉20202| 国产av又大| 精品国产亚洲在线| 99国产精品一区二区蜜桃av| 亚洲国产欧美网| 日本黄色日本黄色录像| 国产aⅴ精品一区二区三区波| 国产亚洲精品第一综合不卡| 久久精品成人免费网站| 久久精品国产综合久久久| 亚洲成av片中文字幕在线观看| 法律面前人人平等表现在哪些方面| 无遮挡黄片免费观看| 黄色片一级片一级黄色片| 精品国产亚洲在线| av超薄肉色丝袜交足视频| av电影中文网址| 亚洲成av片中文字幕在线观看| 黑人猛操日本美女一级片| 国产在线观看jvid| 97超级碰碰碰精品色视频在线观看| 少妇粗大呻吟视频| www.www免费av| 手机成人av网站| 亚洲欧美一区二区三区黑人| 咕卡用的链子| 啦啦啦 在线观看视频| 亚洲情色 制服丝袜| avwww免费| 精品午夜福利视频在线观看一区| 久久性视频一级片| 欧美大码av| 亚洲精品粉嫩美女一区| 国产免费男女视频| 午夜福利在线观看吧| 黄色 视频免费看| 亚洲欧美激情综合另类| 亚洲国产毛片av蜜桃av| 精品欧美一区二区三区在线| 亚洲,欧美精品.| 啦啦啦在线免费观看视频4| 国内毛片毛片毛片毛片毛片| 午夜激情av网站| 欧美成狂野欧美在线观看| 欧美+亚洲+日韩+国产| 精品久久蜜臀av无| 欧美黄色片欧美黄色片| 黑人操中国人逼视频| 亚洲aⅴ乱码一区二区在线播放 | 黄色视频,在线免费观看| a级毛片黄视频| 两性夫妻黄色片| 精品人妻1区二区| 成人三级做爰电影| 在线免费观看的www视频| 久久国产亚洲av麻豆专区| 99国产综合亚洲精品| 色老头精品视频在线观看| 久久精品亚洲av国产电影网| 12—13女人毛片做爰片一| 真人一进一出gif抽搐免费| 久久久久久免费高清国产稀缺| 中国美女看黄片| 黑人巨大精品欧美一区二区mp4| 国产精品 国内视频| 侵犯人妻中文字幕一二三四区| 69精品国产乱码久久久| 成人三级黄色视频| 亚洲情色 制服丝袜| 日韩视频一区二区在线观看| 老司机在亚洲福利影院| 波多野结衣一区麻豆| 亚洲成人精品中文字幕电影 | 啦啦啦免费观看视频1| 色综合欧美亚洲国产小说| 日韩成人在线观看一区二区三区| www国产在线视频色| 亚洲国产毛片av蜜桃av| 自拍欧美九色日韩亚洲蝌蚪91| 啦啦啦免费观看视频1| 精品午夜福利视频在线观看一区| 免费少妇av软件| 狂野欧美激情性xxxx| 一级a爱片免费观看的视频| 久久久国产成人精品二区 | 老司机午夜福利在线观看视频| 亚洲欧美日韩无卡精品| 精品日产1卡2卡| 一级a爱片免费观看的视频| 中出人妻视频一区二区| 在线观看一区二区三区| 男女高潮啪啪啪动态图| 欧美激情 高清一区二区三区| 亚洲情色 制服丝袜| 天堂动漫精品| 国产1区2区3区精品| 成人亚洲精品一区在线观看| 欧美日韩国产mv在线观看视频| 日韩欧美免费精品| a级毛片黄视频| 亚洲五月天丁香| 在线观看免费视频日本深夜| 久久午夜亚洲精品久久| 91字幕亚洲| 日本撒尿小便嘘嘘汇集6| 波多野结衣av一区二区av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品久久成人aⅴ小说| 日韩欧美在线二视频| 国产欧美日韩一区二区三| 欧美久久黑人一区二区| 国产欧美日韩一区二区三区在线| 老汉色av国产亚洲站长工具| 淫妇啪啪啪对白视频| 亚洲精品一二三| www.精华液| 亚洲少妇的诱惑av| 这个男人来自地球电影免费观看| 免费在线观看亚洲国产| 欧美一区二区精品小视频在线| 欧美日韩视频精品一区| 动漫黄色视频在线观看| 欧美乱妇无乱码| 亚洲国产欧美日韩在线播放| 99riav亚洲国产免费| 日韩中文字幕欧美一区二区| 国产精华一区二区三区| 女人被躁到高潮嗷嗷叫费观| 国产成人av教育| 免费看a级黄色片| 久久草成人影院| 久久亚洲精品不卡| 丰满迷人的少妇在线观看| 国产一区二区三区综合在线观看| 丝袜美足系列| 国产亚洲精品久久久久5区| 免费日韩欧美在线观看| 免费看十八禁软件| 丝袜在线中文字幕| 日本撒尿小便嘘嘘汇集6| 国产一区在线观看成人免费| 老司机福利观看| a级片在线免费高清观看视频| 怎么达到女性高潮| av有码第一页| 一本综合久久免费| a级毛片在线看网站| 国产一区二区三区视频了| 18禁黄网站禁片午夜丰满| 国产精品香港三级国产av潘金莲| 国产成人精品久久二区二区免费| 午夜免费激情av| 国产又爽黄色视频| 美女福利国产在线| 一区二区三区精品91| 国产真人三级小视频在线观看| 麻豆一二三区av精品| www.www免费av| 天天影视国产精品| 中文字幕av电影在线播放| 久久精品成人免费网站| 国产精品久久久av美女十八| 国产精品野战在线观看 | 在线国产一区二区在线| 久久精品亚洲av国产电影网| 精品国产一区二区三区四区第35| 色综合婷婷激情| 成年人黄色毛片网站| 国产高清激情床上av| 国产亚洲精品一区二区www| 久9热在线精品视频| 日韩欧美一区二区三区在线观看| 亚洲午夜理论影院| 成人黄色视频免费在线看| 久久精品aⅴ一区二区三区四区| 国产亚洲精品一区二区www| 99香蕉大伊视频| 国产一卡二卡三卡精品| 免费在线观看视频国产中文字幕亚洲| 在线观看免费午夜福利视频| 一进一出好大好爽视频| 在线十欧美十亚洲十日本专区| 久久久精品国产亚洲av高清涩受| 99精品欧美一区二区三区四区| 777久久人妻少妇嫩草av网站| 欧美丝袜亚洲另类 | 美女 人体艺术 gogo|