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

    過表達MRP8提高釀酒酵母乙酸耐性及乙醇發(fā)酵效率

    2017-09-30 01:14:41萬青青趙心清白鳳武
    生物加工過程 2017年5期
    關鍵詞:耐性釀酒乙酸

    萬 春,萬青青,熊 亮,趙心清,白鳳武

    (1.上海交通大學 生命科學技術學院 微生物代謝國家重點實驗室,上海 200240;2.大連理工大學 生命科學與技術學院,遼寧 大連 116024)

    過表達MRP8提高釀酒酵母乙酸耐性及乙醇發(fā)酵效率

    萬 春1,2,萬青青1,熊 亮2,趙心清1,白鳳武1

    (1.上海交通大學 生命科學技術學院 微生物代謝國家重點實驗室,上海 200240;2.大連理工大學 生命科學與技術學院,遼寧 大連 116024)

    乙酸是木質纖維素水解液中含量較多的抑制物,因此提高釀酒酵母菌株對乙酸的耐受性有助于提高纖維素乙醇生產效率。本文中,筆者利用基于CRISPR/Cas9系統(tǒng)的基因組編輯技術過表達了釀酒酵母(Saccharomyces cerevisiae)S288c線粒體核糖體蛋白編碼基因MRP8,并比較了過表達MRP8的菌株與對照菌株的生長和發(fā)酵特性。平板耐性檢測發(fā)現(xiàn),MRP8過表達明顯提高了菌株的乙酸脅迫耐受性;乙醇發(fā)酵結果表明,在4.8 g/L乙酸脅迫條件下,過表達菌株MRP8-3在51 h消耗全部的葡萄糖,發(fā)酵時間縮短了25 h,顯著優(yōu)于相同時間的對照菌株。本研究結果為構建高效纖維素乙醇發(fā)酵的釀酒酵母菌株提供了新思路。

    釀酒酵母;MRP8;乙酸脅迫耐受性;乙醇發(fā)酵;CRISPR/Cas9;基因組編輯

    纖維素乙醇作為清潔可再生能源,是重要的傳統(tǒng)化石燃料替代品,具有良好的應用前景;然而纖維素原料預處理過程產生弱酸類、酚類和醛類等抑制物,顯著影響微生物的生長和發(fā)酵特性[1]。在纖維素預處理過程中產生的抑制物中,乙酸的量最高,可達15 g/L[2]。高濃度的乙酸能引起細胞內線粒體的損傷以及胞內活性氧(reactive oxygen species,ROS)積累,從而降低釀酒酵母生產纖維素乙醇的效率[3]。雖然可通過調節(jié)pH緩解乙酸對釀酒酵母細胞的抑制,但是堿液的使用會提高生產成本。如果釀酒酵母可以在高濃度乙酸抑制條件下正常生長和發(fā)酵,能更有利于工業(yè)化應用[4]。因此,提高酵母細胞對乙酸脅迫的耐受性成為亟需解決的問題。

    近年來,基于CRISPR/Cas9的基因組編輯技術在釀酒酵母中得到快速發(fā)展[5-6]。利用基因組編輯技術可實現(xiàn)快速多基因同時操作,可避免在工業(yè)釀酒酵母基因組中留下抗性標記,因此具有廣泛的應用前景。

    筆者所在課題組前期利用過表達ARO7、GRX5、PPR1和SET5,或者敲除ADY2,構建了乙酸耐性以及乙醇發(fā)酵效率提高的釀酒酵母菌株[7-10]。近期蛋白組分析發(fā)現(xiàn),在乙酸耐性提高的菌株中,Mrp8p表達量顯著上調(未發(fā)表結果)。Mrp8p屬于線粒體核糖體蛋白家族,其編碼基因MRP8位于釀酒酵母XI染色體上[11],當細胞壁出現(xiàn)損傷時,該基因上調1.2~1.6 倍[12],但是,MRP8 過表達對釀酒酵母乙酸耐性的影響目前還不清楚。因此,本文中,筆者利用基于CRISPR/Cas9的基因組編輯技術,在釀酒酵母S288c中過表達MRP8,研究它對乙酸耐性和乙醇發(fā)酵性能的影響,以期為進一步高效纖維素乙醇生產菌株的選育提供理論參考。

    1 材料與方法

    1.1 材料

    1.1.1 菌種與質粒

    釀酒酵母S.cerevisiae S288c和DNA克隆宿主Escherichia coli DH5α保存于微生物代謝國家重點實驗室;Cas9表達質粒為具有諾爾絲菌素抗性的Cas9-NAT[13],gRNA表達質粒為具有潮霉素 B抗性的pRS42H_gRNA_20以及含有MRP8基因的載體pRS42H-Ppgk1-MRP8-Tcyc1,均為筆者所在實驗室構建保存。其中pRS42H_gRNA_20靶向基因間隔區(qū)YPRCdelta15位點,其23堿基靶向序列為GGAATATTGGGTCAGATGAATGG。

    1.1.2 主要試劑和培養(yǎng)基

    諾爾絲菌素,北京中生瑞泰科技有限公司;基因組和質粒提取試劑盒、酵母粉蛋白胨等試劑及其他抗生素等,上海生工生物工程股份有限公司;實時定量PCR試劑盒,日本TaKaRa公司;PCR試劑,北京擎科新業(yè)生物技術有限公司;RT-qPCR儀器,美國Bio-Rad公司;引物委托北京擎科新業(yè)生物技術有限公司合成。

    YPD培養(yǎng)基(g/L):蛋白胨20,葡萄糖20,酵母粉10;

    乙醇發(fā)酵培養(yǎng)基(g/L):蛋白胨 3,葡萄糖 100,酵母粉4。

    1.2 MRP8過表達菌株構建

    用20_Ppgk-F(CTTAAGATGCTCTTCTTATTCTATTAAAAATAGAAAATGAAGGCATTTGCAAGAATTACTCGTG)和20_Tcyc1-R(GTTTGTTTGCGAAACCCTATGCTCTGTTGTTCGGATTTGATCATGTAATTAGTTATGTCACG)引物從pRS42H-Ppgk1-MRP8-Tcyc1載體上擴增轉化元件,其中包含釀酒酵母S288c染色體XVI上整合位點(YPRCdelta15)的上游同源臂、PGK1啟動子、MRP8基因、CYC1終止子以及整合位點的下游同源臂。

    通過醋酸鋰方法[14]將Cas9-NAT質粒轉入釀酒酵母S288c后,用100 mg/L諾爾絲菌素篩選得到轉化子。將pRS42H_gRNA_20質粒和含有MRP8基因純化后的PCR產物以相同方法同時轉化到含有Cas9-NAT質粒的轉化子中;然后用100 mg/L諾爾絲菌素和300 mg/L潮霉素B篩選轉化子。隨機挑選單克隆轉化子命名為MRP81-5,并提取基因組,用Ppgk-F和Tcyc1-R引物通過PCR擴增轉化元件,驗證獲得正確的轉化子;然后將MRP8-1~4在無抗性的YPD培養(yǎng)基中連續(xù)轉接5次(每24 h轉接1次),丟失Cas9-NAT和pRS42H_gRNA_20質粒;然后保存于-80℃冰箱中。

    1.3 菌種活化、耐性測試和發(fā)酵實驗

    1.3.1 菌種活化

    從-80℃冰箱中取出保存的轉化子MRP8-1~4在液體YPD培養(yǎng)基中,在30℃、150 r/min過夜培養(yǎng)至對數(shù)期進行下一步實驗。

    1.3.2 耐性試驗

    將對數(shù)期的轉化子培養(yǎng)液OD調為一致(OD600為1.0),按照10倍梯度稀釋,取 2 μL分別點樣于無脅迫條件的5 g/L乙酸和5 mmol/L H2O2的YPD平板,在菌株充分生長后照相記錄。

    1.3.3 發(fā)酵實驗

    將轉化子MRP8-3活化后轉接到50 mL YPD液體培養(yǎng)基中繼續(xù)培養(yǎng)過夜,然后收集菌體,接種于含100 mL發(fā)酵培養(yǎng)基的250 mL搖瓶中,初始OD600均調節(jié)為0.05左右;在30℃、150 r/min條件下培養(yǎng)。發(fā)酵培養(yǎng)基中不添加或者添加4.8 g/L乙酸。發(fā)酵過程中每3~10 h定時取樣測OD600、葡萄糖和乙醇含量;當葡萄糖低于1.0 g/L時認為發(fā)酵結束。實驗設置2個平行,重復2次。根據(jù)實驗室已建立的方法分析發(fā)酵液中的代謝物[9]。

    1.4 RNA提取及實時定量(RT-qPCR)分析

    利用RNA提取試劑盒提取對數(shù)期菌株的細胞總RNA,根據(jù)實驗室已有的方法進行RT-qPCR分析[9]。由于MSN2和SOD1均和乙酸脅迫響應相關,因此也做了分析;所用引物序列分別為rtMSN2-F(TATCACCATTTCCCACAGCA)和rtMSN2-R(CACCATTGCCAGTTGAGTTG)、rtSOD1-F(TGTCTCTGCTGGTCCTCACTTC)和rtSOD1-R(CACACCATTTTCGTCCGTCTT)、rtMRP8-F(CTGGAAGGTAAGGTGGAAGC)和

    rtMRP8-R(AGATAGCGGGCAACATCATT)。以ALG9為內參基因[15],所用引物為rtALG9-F(ATCGTGAAATTGCAGGCAGCTTGG)和rtALG9-R(CATGGCAACGGCAGAAGGCAATAA),以 2-ΔΔCt方法計算基因相對表達量[15]。

    2 結果與討論

    2.1 基因MRP8過表達對酵母細胞耐受性影響的檢測

    利用基于CRISPR/Cas9的基因組編輯技術將含有MRP8基因的轉化元件定點整合到釀酒酵母S288c染色體XVI的YPRCdelta15位點,并且通過Ppgk-F(GTCGTACGACTGTAATTGCTTTTAG)和Tcyc1-R(GCAAATTAAAGCCTTCGAGCGT)引物進行驗證,表明得到了正確的轉化子;然后通過點板實驗初步評價了隨機挑選的4個轉化子的生長性能。圖1為釀酒酵母S288c中過表達MRP8對生長性能的影響。由圖1可知:在無脅迫條件下,過表達菌株MRP8-1~4的生長性能和對照相差不大;但是當培養(yǎng)基中添加5 g/L乙酸時(圖1(b)),過表達菌株生長狀態(tài)明顯優(yōu)于對照,其中MRP8-2、3和4提高了約1 000倍,MRP8-1提高了約100倍。但是在含有5 mmol/L H2O2的平板上過表達菌株與對照相差不大。

    文獻[16]指出,在YPRCdelta15號位點整合基因具有較高的綠色熒光蛋白表達量,因此本實驗在此強整合位點實現(xiàn)了MRP8的過表達。隨機選擇了轉化子MRP8-3進行發(fā)酵實驗,研究對乙酸脅迫條件下乙醇發(fā)酵的影響。

    圖1 基因MRP8過表達對釀酒酵母細胞生長性能的影響Fig.1 Effects of overexpressing MRP8 on S.cerevisiae S288c growth ability

    2.2 基因MRP8過表達對釀酒酵母發(fā)酵性能的影響

    考察不同濃度乙酸對釀酒酵母菌株的影響,發(fā)現(xiàn)2.4和3.6 g/L乙酸添加條件下,不同菌株的發(fā)酵性能差別不大,推測在這兩個濃度下對釀酒酵母未產生明顯脅迫抑制。而添加5 g/L乙酸導致菌株生長非常緩慢,所以選擇4.8 g/L乙酸脅迫進行液體發(fā)酵實驗比較,先考察無乙酸脅迫條件下釀酒酵母菌株及過表達菌株的生長及發(fā)酵性能,結果見圖2。由圖2可知:在無乙酸脅迫條件下,過表達菌株MRP8-3在培養(yǎng)6~27 h時的OD稍微高于對照組,二者的最大OD600分別為2.98和2.76,但是葡萄糖消耗和乙醇產量差別不大。其中,選取對數(shù)期(9 h)的酵母細胞進行RT-qPCR分析,過表達菌株MRP8-3和對照菌株的 OD600分別為 1.38 和 1.00。

    考察在4.8 g/L乙酸脅迫條件下,過表達菌株與對照菌株的發(fā)酵性能,結果如圖3所示。由圖3可知:過表達菌株MRP8-3的延滯期比對照菌株縮短24 h;選取對數(shù)期的酵母細胞用于RT-qPCR分析,其中過表達菌株MRP8-3和對照菌株的取樣時間點分別為29 和 54 h,對應的 OD600值分別為1.11 和 1.01,盡量選擇OD值相近的酵母細胞進行RT-qPCR,以降低不同生長狀態(tài)對實驗結果的影響。由于發(fā)酵后期酵母細胞在乙酸脅迫條件下出現(xiàn)了聚集現(xiàn)象,導致OD值誤差較大。但是從葡萄糖和乙醇含量變化亦能觀察到過表達菌株MRP8-3的發(fā)酵性能明顯優(yōu)于對照菌株,并且在51 h消耗完所有葡萄糖。過表達菌株MRP8-3在54 h達到最高乙醇質量濃度43.59 g/L,略低于對照菌株在76 h的44.73 g/L,但對照菌株在76 h仍殘留中仍剩余約5 g/L葡萄糖。

    圖2 過表達菌株與對照菌株在無乙酸條件下的發(fā)酵性能比較Fig.2 Comparison of fermentation performances between overexpression strains and control strain without acetic acid

    圖3 過表達菌株與對照菌株在4.8 g/L乙酸條件下的發(fā)酵性能比較Fig.3 Comparison of fermentation performances between overexpression strains and control strain with 4.8 g/L acetic acid

    表1為過表達菌株和對照菌株在含有4.8 g/L乙酸的培養(yǎng)基中乙醇發(fā)酵的各指標。由表1可以看出:過表達MRP8基因可明顯縮短發(fā)酵時間,提高乙醇生產強度,且對乙醇得率影響不大。細胞通過合成分泌甘油來抵抗外界脅迫[17],因此在過表達菌株和對照菌株的發(fā)酵液中均檢測到甘油;此外,對照菌株中的甘油含量比過表達菌株高出了約86%,表明MRP8基因的過表達能提高細胞對乙酸的耐性,無需合成較多的甘油作為保護劑。由此也說明Mrp8p可能在酵母細胞抵抗乙酸脅迫中發(fā)揮重要作用。

    表1 在4.8 g/L乙酸脅迫條件下MRP8過表達菌株和對照菌株的發(fā)酵性能匯總Table 1 Summary of fermentation performance between MRP8 overexpressing strains and control strain with 4.8 g/L acetic acid

    2.3 RT-qPCR分析

    圖4為RT-qPCR分析過表達MRP8菌株中基因表達量結果。由圖4可知:在無乙酸脅迫條件下,過表達菌株MRP8-3中MRP8相對表達量稍微提高,在4.8 g/L乙酸脅迫條件下,MRP8相對表達量在過表達菌株MRP8-3中提高了約11倍,這揭示該基因轉錄受到乙酸脅迫的誘導。過表達菌株MRP8-3中SOD1表達量在無脅迫條件下沒有太大變化,而在乙酸脅迫條件下明顯下調。高濃度乙酸會引起細胞器損傷,脅迫DNA復制而加速細胞衰老[3,18],同時高濃度乙酸會導致胞內活性氧自由基(ROS)積累[19],而與胞內氧化脅迫有關的SOD1的相對表達量明顯下調,表明Mrp8p過表達有利于清除胞內ROS,因此不需要過表達SOD1。過表達菌株MRP8-3中MSN2基因在非脅迫條件下顯著上調,但在乙酸脅迫條件下 MSN2變化不明顯。Msn2p是脅迫響應相關的轉錄調控因子[20],推測MSN2在MRP8過表達菌株脅迫條件下的耐性提高沒有起到關鍵的調控作用。此外,MRP8表達量與Bur6p、Cst6p、Gcn5p 和 Sua7p 調節(jié)因子有關[21-22],這些調節(jié)因子也是以后進一步研究該基因調控機理的重要內容。

    圖4 RT-qPCR分析過表達MRP8菌株中基因表達量Fig.4 Analysis of genes expression level in the overexpressing MRP8 strains via RT-qPCR

    3 結論

    本研究首次報道了過表達MRP8基因能提高釀酒酵母的乙酸耐受性,同時提高了乙酸脅迫條件下的乙醇發(fā)酵水平。在4.8 g/L乙酸脅迫條件下,過表達菌株乙醇強度最高達到0.85 g/(L·h),比對照菌株提高了44%。本實驗室已發(fā)現(xiàn)了多種提高乙酸耐性的基因操作靶點,本研究結果進一步豐富了酵母細胞乙酸耐性機理,對MRP8的深入研究有望進一步完善釀酒酵母的乙酸耐性機理,構建穩(wěn)定高效的纖維素乙醇發(fā)酵菌株。

    [1] ZHAO X,XIONG L,ZHANG M,etal.Towardsefficient bioethanol production from agricultural and forestry residues:exploration of unique natural microorganisms in combination with advanced strain engineering[J].Bioresour Technol,2016,215:84-91.

    [2] WEI N,QUARTERMAN J,KIM S R,et al.Enhanced biofuel production through coupled acetic acid and xylose consumption by engineered yeast[J].Nat Commun,2013,4:2580.

    [3] 趙心清,張明明,徐桂紅,等.釀酒酵母乙酸耐性分子機制的功能基因組進展[J].生物工程學報,2014,30(3):368-380.

    [4] GRAVES T,NARENDRANATH N V,DAWSON K,et al.Effect of pH and lactic oracetic acid on ethanolproductivity by Saccharomyces cerevisiae in corn mash[J].J Ind Microbiol Biotechnol,2006,33(6):469-474.

    [5] APEL A R,D'ESPAUX L,WEHRS M,et al.A Cas9-based toolkit to program gene expression in Saccharomyces cerevisiae[J].Nucleic Acids Res,2017,45(1):496-508.

    [6] 謝澤雄,李炳志,沈玥,等.酵母基因組工程[J].中國科學:生命科學,2015,45(10):985-992.

    [7] ZHANG M M,ZHAO X Q,CHENG C,et al.Improved growth and ethanol fermentation of Saccharomyces cerevisiae in the presence of acetic acid by overexpression of SET5 and PPR1[J].Biotechnol J,2015,10(12):1903-1911.

    [8] ZHANG M,ZHANG K,MEHMOOD M A,et al.Deletion of acetate transporter gene ADY2 improved tolerance of Saccharomyces cerevisiae against multiple stresses and enhanced ethanol production in the presence of acetic acid[J].Bioresour Technol,2017,doi:10.1016/j.biortech.2017.05.191.

    [9] 張明明,萬青青,張克俞,等.過表達分支酸歧化酶編碼基因ARO7對釀酒酵母抑制物耐受性的影響[J].應用與環(huán)境生物學報,2016,22(2):201-205.

    [10] 方青,張明明,陳洪奇,等.過表達谷氧還蛋白基因GRX5提高釀酒酵母乙酸耐性[J].化工學報,2015,66(4):1434-1439.

    [11] ABRAHAM P R,MULDER A,VAN'T RIET J,et al.Molecular cloning and physicalanalysisofan 8.2 kb segmentof chromosome XI of Saccharomyces cerevisiae reveals five tightly linked genes[J].Yeast,1992,8(3):227-238.

    [12] BOORSMA A,DE NOBEL H,TER RIET B,etal.Characterization of the transcriptional response to cell wall stress in Saccharomyces cerevisiae[J].Yeast,2004,21(5):413-427.

    [13] ZHANG G C,KONG I I,KIM H,et al.Construction of a quadruple auxotrophic mutant of an industrial polyploid Saccharomycescerevisiaestrain byusingRNA-guided Cas9 nuclease[J].Appl Environ Microbiol,2014,80(24):7694-7701.

    [14] GIETZ R D,SCHIESTL R H.High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method[J].Nat Protoc,2007,2(1):31-34.

    [15] TESTE M A,DUQUENNE M,F(xiàn)RANCOIS J M,et al.Validation of reference genes for quantitative expression analysis by real-time RT-PCR in Saccharomyces cerevisiae[J].BMC Mol Biol,2009,10:99.

    [16] FLAGFELDT D,SIEWERS V,HUANG L,et al.Characterization of chromosomal integration sites for heterologous gene expression in Sacchromyces cerevisiae[J].Yeast,2009,26:545-551.

    [17] MEDINA V G,ALMERING M J,VAN MARIS A J,et al.Elimination of glycerol production in anaerobic cultures of a Saccharomyces cerevisiae strain engineered to use acetic acid as an electron acceptor[J].Appl Environ Microbiol,2010,76(1):190-195.

    [18] MURAKAMI C,DELANEY J R,CHOU,A et al.pH neutralization protectsagainstreduction in replicative lifespan following chronological aging in yeast[J].Cell Cycle,2012,11(16):3087-3096.

    [19] WAN C,ZHAO X,ZHANG M,et al.Impact of zinc sulfate addition on dynamic metabolic profiling ofSaccharomyces cerevisiae subjected to long term acetic acid stress treatment and identification of key metabolites involved in antioxidant effect of zinc[J].Metallomics,2015,7(2):322-332.

    [20] SCHMITT A P,MCENTEE K.Msn2p,a zinc finger DNA-binding protein,is the transcriptional activator of the multistress response in Saccharomyces cerevisiae[J].Proc Natl Acad Sci USA,1996,93(12):5777-5782.

    [21] HU Z,KILLION P J,IYER V R.Genetic reconstruction of a functional transcriptional regulatory network[J].Nat Genet,2007,39(5):683-687.

    [22] VENTERSBJ,WACHIS,MAVRICH TN,etal.A comprehensive genomic binding map of gene and chromatin regulatory proteins in Saccharomyces[J].Mol Cell,2011,41(4):480-492.

    (責任編輯 荀志金)

    Improvement of acetic acid tolerance of Saccharomyces cerevisiae by overexpression of mitochondrial ribosomal protein encoding gene MRP8 for efficient lignocellulosic ethanol production

    WAN Chun1,2,WAN Qingqing1,XIONG Liang2,ZHAO Xinqing1,BAI Fengwu1
    (1.State Key Laboratory of Microbial Metabolism,School of Life Sciences and Biotechnology,Shanghai Jiao Tong University,Shanghai 200240,China;2.School of Life Science and Biotechnology,Dalian University of Technology,Dalian 116024,China)

    Acetic acid is one of the major inhibitors derived from the hydrolysis of lignocellulosic biomass,and enhancing the tolerance of Saccharomyces cerevisiae to acetic acid is of great importance for improving the lignocellulosic ethanol production.Here,we overexpressed mitochondrial ribosomal protein encoding gene MRP8 via CRISPR/Cas9-based genome editing in S.cerevisiae S288c,and compared the growth and fermentation performance between the mutants with MRP8 and the wild type.The tolerance to acetic acid was increased by overexpression MRP8 on agar plates.Also,mutant MRP8-3 consumed all glucose within 51 h in the presence of 4.8 g/L acetic acid,and the fermentation time was shortened by 25 h,indicating the mutant was superior to the wild type.Our findings provide a novel strategy for constructing robust yeast strains for efficient lignocellulosic ethanol fermentation.

    Saccharomyces cerevisiae;MRP8;acetic acid tolerance;ethanol fermentation;CRISPR/Cas9;genome editing

    TQ920.1

    A

    1672-3678(2017)05-0080-06

    10.3969/j.issn.1672-3678.2017.05.010

    2017-07-14

    國家自然科學基金(21376043、31461143029);國家高技術研究發(fā)展計劃(863計劃)(2012AA021205、2012AA101805)

    萬 春(1989—),男,河南信陽人,博士,研究方向:生物化工;趙心清(聯(lián)系人),教授,E-mail:xqzhao@sjtu.edu.cn

    猜你喜歡
    耐性釀酒乙酸
    乙醇和乙酸常見考點例忻
    上半年釀酒產業(yè)產、銷、利均增長
    釀酒科技(2021年8期)2021-12-06 15:28:22
    為什么酵母菌既能做面包也能釀酒?
    沒耐性的小猴
    快樂語文(2020年26期)2020-10-16 03:37:00
    DMAC水溶液乙酸吸附分離過程
    釀酒忘米
    乙酸仲丁酯的催化合成及分析
    調虧灌溉在釀酒葡萄上的應用
    2-(N-甲氧基)亞氨基-2-苯基乙酸異松蒎酯的合成及表征
    水稻苗期不同階段與低氮耐性相關的QTL分析
    美女xxoo啪啪120秒动态图 | 国产av在哪里看| 亚洲色图av天堂| 婷婷精品国产亚洲av| 动漫黄色视频在线观看| 亚洲人成网站在线播放欧美日韩| 麻豆成人av在线观看| 在线观看免费视频日本深夜| 国产亚洲精品久久久久久毛片| 永久网站在线| 日本免费一区二区三区高清不卡| 久久精品人妻少妇| 亚洲精品亚洲一区二区| av天堂在线播放| 美女免费视频网站| 国产午夜精品久久久久久一区二区三区 | 国产色爽女视频免费观看| 可以在线观看毛片的网站| 99久久精品一区二区三区| 看免费av毛片| 18+在线观看网站| 一个人看视频在线观看www免费| 日韩精品中文字幕看吧| 免费av不卡在线播放| 动漫黄色视频在线观看| 男人舔奶头视频| 搡老岳熟女国产| 精品一区二区免费观看| 男女床上黄色一级片免费看| 亚洲精品粉嫩美女一区| 欧美日韩乱码在线| 69av精品久久久久久| 美女cb高潮喷水在线观看| 18禁黄网站禁片午夜丰满| 狂野欧美白嫩少妇大欣赏| 99久久精品国产亚洲精品| 高清在线国产一区| 亚洲三级黄色毛片| 一本综合久久免费| 亚洲色图av天堂| 国产极品精品免费视频能看的| 99国产极品粉嫩在线观看| www日本黄色视频网| 黄片小视频在线播放| 亚洲国产欧美人成| а√天堂www在线а√下载| 老鸭窝网址在线观看| 美女cb高潮喷水在线观看| 一级黄片播放器| 亚洲avbb在线观看| 国产高清三级在线| 欧美日本视频| 欧美激情久久久久久爽电影| 看十八女毛片水多多多| av天堂中文字幕网| 51午夜福利影视在线观看| 国产精品久久视频播放| 禁无遮挡网站| 精品熟女少妇八av免费久了| 首页视频小说图片口味搜索| 在线播放国产精品三级| 国产在线男女| 中文资源天堂在线| 免费无遮挡裸体视频| 一区二区三区高清视频在线| 欧美+亚洲+日韩+国产| 婷婷色综合大香蕉| avwww免费| 美女高潮喷水抽搐中文字幕| 我要看日韩黄色一级片| or卡值多少钱| 一边摸一边抽搐一进一小说| 可以在线观看的亚洲视频| 成人亚洲精品av一区二区| 精品一区二区三区av网在线观看| 日日摸夜夜添夜夜添av毛片 | 女同久久另类99精品国产91| 国产乱人视频| 婷婷六月久久综合丁香| 我的女老师完整版在线观看| 亚洲av日韩精品久久久久久密| 高潮久久久久久久久久久不卡| 波多野结衣巨乳人妻| av福利片在线观看| 神马国产精品三级电影在线观看| 香蕉av资源在线| 偷拍熟女少妇极品色| 日韩欧美精品免费久久 | 午夜影院日韩av| 亚州av有码| 日本a在线网址| 欧美成狂野欧美在线观看| 欧美中文日本在线观看视频| 精品久久久久久久久久免费视频| 少妇熟女aⅴ在线视频| 麻豆成人av在线观看| 精品久久久久久久久亚洲 | 亚洲黑人精品在线| 国产精品美女特级片免费视频播放器| 亚洲久久久久久中文字幕| 麻豆一二三区av精品| 又黄又爽又免费观看的视频| 99久久久亚洲精品蜜臀av| 最新中文字幕久久久久| 变态另类丝袜制服| 夜夜爽天天搞| 99精品久久久久人妻精品| 中文字幕人妻熟人妻熟丝袜美| 亚洲一区二区三区不卡视频| .国产精品久久| 2021天堂中文幕一二区在线观| 最新中文字幕久久久久| 国产精品爽爽va在线观看网站| 亚洲成av人片免费观看| 欧美中文日本在线观看视频| 少妇的逼水好多| 国产精品亚洲一级av第二区| 美女大奶头视频| 国产成人福利小说| 日韩高清综合在线| 校园春色视频在线观看| 亚洲,欧美,日韩| 又粗又爽又猛毛片免费看| 99久久精品热视频| 搡老妇女老女人老熟妇| bbb黄色大片| 久久久久久久亚洲中文字幕 | 欧美成人a在线观看| 色视频www国产| 亚洲午夜理论影院| 日本黄大片高清| 亚洲综合色惰| av视频在线观看入口| 99热这里只有是精品在线观看 | 精品久久国产蜜桃| 18美女黄网站色大片免费观看| 久久久国产成人精品二区| 午夜视频国产福利| 女同久久另类99精品国产91| 男人舔女人下体高潮全视频| 免费看美女性在线毛片视频| 欧美xxxx性猛交bbbb| 怎么达到女性高潮| 伦理电影大哥的女人| 国产三级在线视频| 色哟哟·www| 日韩大尺度精品在线看网址| 亚洲va日本ⅴa欧美va伊人久久| 欧美又色又爽又黄视频| 亚洲内射少妇av| 最近在线观看免费完整版| 精品久久久久久久久久久久久| 成人高潮视频无遮挡免费网站| 亚洲精品在线观看二区| 国产一区二区激情短视频| 成人av在线播放网站| 亚洲天堂国产精品一区在线| 欧美激情国产日韩精品一区| 国产极品精品免费视频能看的| 亚洲熟妇中文字幕五十中出| 一个人免费在线观看电影| 在线免费观看不下载黄p国产 | 天堂av国产一区二区熟女人妻| 精品国产亚洲在线| 精华霜和精华液先用哪个| 国产精品日韩av在线免费观看| 亚洲五月婷婷丁香| 亚洲欧美精品综合久久99| 亚洲午夜理论影院| 三级国产精品欧美在线观看| 亚州av有码| 午夜福利在线观看吧| 成人亚洲精品av一区二区| 好看av亚洲va欧美ⅴa在| av在线观看视频网站免费| or卡值多少钱| 国产亚洲精品久久久com| 国产在线精品亚洲第一网站| 色综合婷婷激情| 久久欧美精品欧美久久欧美| 精品久久国产蜜桃| 亚洲最大成人手机在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 无遮挡黄片免费观看| 一个人免费在线观看的高清视频| 国产淫片久久久久久久久 | 一区二区三区激情视频| 国产三级在线视频| 亚洲内射少妇av| 中国美女看黄片| 非洲黑人性xxxx精品又粗又长| 久久亚洲真实| 最新在线观看一区二区三区| 国产又黄又爽又无遮挡在线| 看十八女毛片水多多多| 亚洲美女黄片视频| 在线观看美女被高潮喷水网站 | 丁香六月欧美| 国产精品久久久久久久电影| 久久99热这里只有精品18| 嫩草影院入口| 老女人水多毛片| 中文资源天堂在线| 亚洲国产精品sss在线观看| 亚洲性夜色夜夜综合| 久久99热6这里只有精品| 久久精品国产亚洲av涩爱 | 女人被狂操c到高潮| 在线十欧美十亚洲十日本专区| 欧美3d第一页| 99国产精品一区二区三区| 少妇丰满av| 性插视频无遮挡在线免费观看| 亚洲欧美激情综合另类| 国产综合懂色| 国产老妇女一区| 精品乱码久久久久久99久播| 性色av乱码一区二区三区2| 99热这里只有精品一区| 一个人看的www免费观看视频| 国产精品亚洲一级av第二区| 国内少妇人妻偷人精品xxx网站| 久久精品国产清高在天天线| 在线观看美女被高潮喷水网站 | 久久久久国产精品人妻aⅴ院| 精品国产三级普通话版| 久久性视频一级片| 日日夜夜操网爽| 日韩精品青青久久久久久| 中文字幕熟女人妻在线| 日韩 亚洲 欧美在线| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 亚洲18禁久久av| 欧美色视频一区免费| 免费观看精品视频网站| 亚洲成人久久性| 床上黄色一级片| 国产中年淑女户外野战色| 麻豆一二三区av精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲精品粉嫩美女一区| 国内精品久久久久久久电影| 特级一级黄色大片| 一a级毛片在线观看| 国产精品久久视频播放| 十八禁国产超污无遮挡网站| 不卡一级毛片| 琪琪午夜伦伦电影理论片6080| x7x7x7水蜜桃| 久久精品国产亚洲av涩爱 | 欧美成人一区二区免费高清观看| 又爽又黄a免费视频| 夜夜看夜夜爽夜夜摸| 天堂√8在线中文| 国产成人影院久久av| 日本黄色片子视频| 白带黄色成豆腐渣| 88av欧美| 97超级碰碰碰精品色视频在线观看| 国内精品一区二区在线观看| www.色视频.com| 色尼玛亚洲综合影院| 深夜a级毛片| 国产精品人妻久久久久久| 国产毛片a区久久久久| 久久久久久国产a免费观看| 天美传媒精品一区二区| 在线观看午夜福利视频| 午夜福利在线观看吧| 99热这里只有是精品50| 999久久久精品免费观看国产| 免费高清视频大片| 能在线免费观看的黄片| 午夜视频国产福利| 2021天堂中文幕一二区在线观| 搡老妇女老女人老熟妇| 久久草成人影院| av福利片在线观看| 欧美精品国产亚洲| 别揉我奶头~嗯~啊~动态视频| 成人精品一区二区免费| 美女大奶头视频| 久久亚洲真实| 十八禁国产超污无遮挡网站| 黄色配什么色好看| 久久这里只有精品中国| 免费看美女性在线毛片视频| 特大巨黑吊av在线直播| 国产午夜精品论理片| 人妻夜夜爽99麻豆av| 久久中文看片网| www.999成人在线观看| 中文字幕精品亚洲无线码一区| 国产单亲对白刺激| 99视频精品全部免费 在线| 亚洲真实伦在线观看| 一个人观看的视频www高清免费观看| 午夜久久久久精精品| 少妇高潮的动态图| 国产乱人视频| 精品日产1卡2卡| 国产精品综合久久久久久久免费| 久久香蕉精品热| 亚洲av免费在线观看| 精品一区二区三区视频在线| 18禁黄网站禁片免费观看直播| 18禁裸乳无遮挡免费网站照片| 国产在视频线在精品| 在线观看免费视频日本深夜| 在线观看一区二区三区| 日韩欧美三级三区| 狠狠狠狠99中文字幕| 日韩欧美国产在线观看| 日韩免费av在线播放| 18禁在线播放成人免费| 日本黄色片子视频| 国产黄片美女视频| 国产精品av视频在线免费观看| 性插视频无遮挡在线免费观看| a级毛片免费高清观看在线播放| 精品99又大又爽又粗少妇毛片 | 亚洲天堂国产精品一区在线| 国产高清视频在线播放一区| 九九热线精品视视频播放| 日韩欧美在线二视频| 成人毛片a级毛片在线播放| 久久精品人妻少妇| 亚洲国产精品成人综合色| 国产大屁股一区二区在线视频| 黄色配什么色好看| 精品午夜福利在线看| 一级作爱视频免费观看| a级一级毛片免费在线观看| 麻豆久久精品国产亚洲av| 国产精品,欧美在线| 国内精品一区二区在线观看| eeuss影院久久| 亚洲内射少妇av| 性色av乱码一区二区三区2| 精品一区二区三区视频在线| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 韩国av一区二区三区四区| 成人性生交大片免费视频hd| 成年免费大片在线观看| 国产伦在线观看视频一区| 日本a在线网址| 国产精品一区二区免费欧美| 国产午夜精品久久久久久一区二区三区 | 日韩欧美精品免费久久 | 免费看美女性在线毛片视频| 国产高清三级在线| 丁香欧美五月| 99在线视频只有这里精品首页| 国产毛片a区久久久久| 九色成人免费人妻av| 中文字幕免费在线视频6| 国产精品综合久久久久久久免费| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 亚洲精品乱码久久久v下载方式| 亚洲va日本ⅴa欧美va伊人久久| 国内少妇人妻偷人精品xxx网站| 啪啪无遮挡十八禁网站| 国产视频内射| 一区二区三区激情视频| 99在线人妻在线中文字幕| 精品99又大又爽又粗少妇毛片 | 亚洲美女黄片视频| 午夜老司机福利剧场| 很黄的视频免费| 全区人妻精品视频| 噜噜噜噜噜久久久久久91| 久久精品国产99精品国产亚洲性色| 国产精品乱码一区二三区的特点| 波多野结衣高清作品| 色av中文字幕| 看免费av毛片| www.色视频.com| 亚洲av免费高清在线观看| av专区在线播放| 国产精品永久免费网站| 久久人人精品亚洲av| 午夜亚洲福利在线播放| 欧美xxxx性猛交bbbb| 人人妻人人看人人澡| 国产欧美日韩一区二区三| 日韩亚洲欧美综合| 免费无遮挡裸体视频| 色综合站精品国产| 99久久成人亚洲精品观看| 蜜桃亚洲精品一区二区三区| 成人毛片a级毛片在线播放| 色哟哟·www| 搡老妇女老女人老熟妇| 天堂影院成人在线观看| 国语自产精品视频在线第100页| 久久亚洲精品不卡| 性插视频无遮挡在线免费观看| 成人av在线播放网站| 亚洲专区国产一区二区| 午夜福利免费观看在线| 精品人妻一区二区三区麻豆 | 国产精品98久久久久久宅男小说| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美98| 国产精品一区二区三区四区免费观看 | 毛片女人毛片| 国产一区二区在线av高清观看| 日本三级黄在线观看| 网址你懂的国产日韩在线| 日韩欧美精品v在线| 精品一区二区三区视频在线| 99国产精品一区二区三区| 一个人看视频在线观看www免费| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成网站在线播| 国产乱人视频| 嫩草影院入口| 亚洲av电影不卡..在线观看| 国产男靠女视频免费网站| 99热只有精品国产| 别揉我奶头 嗯啊视频| 悠悠久久av| 成人无遮挡网站| 国产亚洲欧美98| 亚洲专区中文字幕在线| www.www免费av| 黄色一级大片看看| 日韩中字成人| 很黄的视频免费| 国产高潮美女av| 亚洲国产欧洲综合997久久,| 精品午夜福利视频在线观看一区| 国产精品久久久久久久电影| 国产伦人伦偷精品视频| 国产精品亚洲av一区麻豆| 精品人妻1区二区| 亚洲精品久久国产高清桃花| 人妻制服诱惑在线中文字幕| 永久网站在线| 亚洲精华国产精华精| 国产av一区在线观看免费| 午夜福利在线在线| 中亚洲国语对白在线视频| 99久久九九国产精品国产免费| 精品人妻视频免费看| 国产精品久久视频播放| 少妇高潮的动态图| 国产伦精品一区二区三区视频9| 免费观看的影片在线观看| 免费无遮挡裸体视频| 久久这里只有精品中国| 日韩欧美三级三区| 别揉我奶头 嗯啊视频| 欧美bdsm另类| 久久久久久久午夜电影| 99热精品在线国产| 亚洲av美国av| 少妇人妻精品综合一区二区 | 好男人在线观看高清免费视频| 亚洲精品一卡2卡三卡4卡5卡| 听说在线观看完整版免费高清| 久久99热6这里只有精品| 亚洲成人久久性| 禁无遮挡网站| 日韩国内少妇激情av| 三级国产精品欧美在线观看| 男插女下体视频免费在线播放| 网址你懂的国产日韩在线| 午夜影院日韩av| 中文字幕精品亚洲无线码一区| 老熟妇仑乱视频hdxx| 久久久久久国产a免费观看| 极品教师在线视频| 亚洲中文字幕一区二区三区有码在线看| 欧美成人一区二区免费高清观看| 欧美日本视频| 日韩欧美三级三区| 真实男女啪啪啪动态图| 一进一出好大好爽视频| 久久草成人影院| 亚洲自拍偷在线| 免费看日本二区| 少妇人妻一区二区三区视频| 亚洲精品影视一区二区三区av| 国产精品1区2区在线观看.| 观看美女的网站| 色5月婷婷丁香| 亚洲综合色惰| www.熟女人妻精品国产| 国产精品久久久久久精品电影| 91在线精品国自产拍蜜月| 亚洲最大成人中文| 亚洲熟妇中文字幕五十中出| 天天躁日日操中文字幕| 看片在线看免费视频| 琪琪午夜伦伦电影理论片6080| 精品人妻一区二区三区麻豆 | 非洲黑人性xxxx精品又粗又长| 国产又黄又爽又无遮挡在线| 我的女老师完整版在线观看| 精品久久久久久成人av| 欧美三级亚洲精品| 男女视频在线观看网站免费| 日本成人三级电影网站| 最近中文字幕高清免费大全6 | 亚洲aⅴ乱码一区二区在线播放| 国产成人啪精品午夜网站| 97人妻精品一区二区三区麻豆| 国产精品伦人一区二区| av在线蜜桃| 亚洲av免费在线观看| 精品午夜福利视频在线观看一区| www日本黄色视频网| 99热这里只有精品一区| 国产欧美日韩精品一区二区| 变态另类丝袜制服| 免费无遮挡裸体视频| a级毛片a级免费在线| 欧美最黄视频在线播放免费| 99热只有精品国产| 亚洲欧美日韩高清在线视频| 午夜福利18| 青草久久国产| 亚洲专区中文字幕在线| 午夜激情福利司机影院| 久久人人精品亚洲av| 成人特级av手机在线观看| 国产精品不卡视频一区二区 | 国产高清三级在线| 久久久久久久久大av| 九色国产91popny在线| 长腿黑丝高跟| 国产在视频线在精品| 精品一区二区免费观看| 日本一本二区三区精品| 婷婷六月久久综合丁香| 亚洲天堂国产精品一区在线| 国产爱豆传媒在线观看| 日本撒尿小便嘘嘘汇集6| 久久人人精品亚洲av| 欧美成人性av电影在线观看| 白带黄色成豆腐渣| 18美女黄网站色大片免费观看| 999久久久精品免费观看国产| 女生性感内裤真人,穿戴方法视频| 18+在线观看网站| 51午夜福利影视在线观看| 国产乱人伦免费视频| 亚洲成av人片免费观看| 国产乱人视频| 午夜福利高清视频| 国产日本99.免费观看| 国产真实伦视频高清在线观看 | 国产亚洲av嫩草精品影院| 中文字幕av成人在线电影| 欧美中文日本在线观看视频| 亚洲国产精品久久男人天堂| 国产成年人精品一区二区| 欧美潮喷喷水| 国产精品免费一区二区三区在线| 18美女黄网站色大片免费观看| 日本成人三级电影网站| 亚洲国产色片| 国产精品影院久久| 真人一进一出gif抽搐免费| 国产av在哪里看| 欧美日韩福利视频一区二区| 99国产极品粉嫩在线观看| 国产欧美日韩一区二区精品| 午夜两性在线视频| 亚洲成人久久爱视频| 欧美最新免费一区二区三区 | 久久久久久国产a免费观看| 国产午夜精品论理片| 69av精品久久久久久| 亚洲av成人精品一区久久| 90打野战视频偷拍视频| 夜夜躁狠狠躁天天躁| 嫩草影院精品99| av黄色大香蕉| 成年女人永久免费观看视频| 日本熟妇午夜| 久久久久九九精品影院| av在线蜜桃| 国产69精品久久久久777片| 日韩欧美 国产精品| 性欧美人与动物交配| 91麻豆精品激情在线观看国产| 久久草成人影院| 大型黄色视频在线免费观看| 极品教师在线视频| 国产主播在线观看一区二区| 男女那种视频在线观看| 亚洲最大成人av| 国产成人aa在线观看| 国产不卡一卡二| 国产色婷婷99| 国内毛片毛片毛片毛片毛片| 女同久久另类99精品国产91| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| 在现免费观看毛片| 国产野战对白在线观看| 国产精品永久免费网站| 日本与韩国留学比较| 亚洲美女搞黄在线观看 | 国产精品野战在线观看| 最近中文字幕高清免费大全6 | 一本综合久久免费| 九色成人免费人妻av| 天天一区二区日本电影三级| 国产亚洲精品久久久com|