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

    毛果楊Rubisco活化酶基因的克隆與功能分析*

    2017-06-05 09:07:55孫偉博諸葛強
    林業(yè)科學(xué) 2017年4期
    關(guān)鍵詞:飽和點光化學(xué)株系

    尹 吳 孫偉博 周 燕 諸葛強

    (南京林業(yè)大學(xué) 南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心 林木遺傳與生物技術(shù)省部共建教育部重點實驗室 南京 210037)

    ?

    毛果楊Rubisco活化酶基因的克隆與功能分析*

    尹 吳 孫偉博 周 燕 諸葛強

    (南京林業(yè)大學(xué) 南方現(xiàn)代林業(yè)協(xié)同創(chuàng)新中心 林木遺傳與生物技術(shù)省部共建教育部重點實驗室 南京 210037)

    【目的】 核酮糖-1, 5-二磷酸羧化酶(Rubisco)是參與植物光合作用的第1步碳同化的關(guān)鍵酶,而Rubisco活化酶(RCA)能夠使Rubisco處于穩(wěn)定的催化活性狀態(tài),從而提高光合效率。本研究從毛果楊中克隆RCA基因并通過遺傳轉(zhuǎn)化南林895楊,獲得PtRCA高表達(dá)的轉(zhuǎn)基因株系,并進(jìn)行分子檢測及相關(guān)功能分析,為培育楊樹新型高光效抗逆品種提供依據(jù)?!痉椒ā?基于毛果楊基因組數(shù)據(jù)庫信息克隆PtRCA基因序列,利用生物信息學(xué)對PtRCA基因進(jìn)行功能和結(jié)構(gòu)分析。采用Gateway技術(shù)將PtRCA構(gòu)建到植物表達(dá)載體pGWB406中,以南林895楊為受體材料,通過農(nóng)桿菌介導(dǎo)法進(jìn)行遺傳轉(zhuǎn)化。以轉(zhuǎn)PtRCA基因南林895楊和對照(南林895楊)為材料,測定分析在高溫脅迫下轉(zhuǎn)基因植株的基因表達(dá)、光合參數(shù)和葉綠素?zé)晒鈪?shù)的差異?!窘Y(jié)果】 從毛果楊中克隆獲得PtRCA的CDS序列長1 323 bp,編碼440個氨基酸殘基,其蛋白相對分子質(zhì)量為48 315.9 Da,等電點pI為5.57,為疏水性蛋白,無信號肽以及跨膜結(jié)構(gòu); 通過序列對比,發(fā)現(xiàn)PtRCA屬AAA+超級家族一員,且與大豆、擬南芥RCA蛋白同源性較高。轉(zhuǎn)PtRCA基因南林895楊RCA表達(dá)量均高于對照,且能夠利用強光充分進(jìn)行光合作用,其光飽和點也均高于對照,增加約12.5%~37.5%; 光補償點除3號株系外,其他轉(zhuǎn)基因株系均略低于對照; 轉(zhuǎn)基因植株在光飽和點的光合速率比對照增高24.6%~55.7%。轉(zhuǎn)基因株系對CO2的利用能力和羧化效率均強于對照組,CO2飽和點除1號和4號株系外,其他株系均比對照低12.5%~25.0%; CO2補償點比對照降低53.1%~80.4%; 光呼吸比對照低37.7%~79.3%,并且轉(zhuǎn)基因楊樹在CO2飽和點的光合速率比對照增高4.4%~26.4%。另外,轉(zhuǎn)基因楊樹表現(xiàn)出耐光氧化的能力,光氧化處理后,對照PSⅡ原初光化學(xué)效率下降61.7%,而轉(zhuǎn)基因楊樹下降45.0%~53.1%; 對照PSⅡ?qū)嶋H光化學(xué)效率下降54.1%,而轉(zhuǎn)基因楊樹下降38.7%~52.0%; 光化學(xué)猝滅系數(shù)對照下降68.3%,而轉(zhuǎn)基因楊樹下降51.0%~65.8%; 非光化學(xué)猝滅系數(shù)對照僅增加3.0%,而轉(zhuǎn)基因楊樹增加6.0%~26.5%?!窘Y(jié)論】 楊樹Rubisco活化酶(PtRCA)蛋白與大豆、擬南芥RCA蛋白同源性較高。通過實時定量及相關(guān)生理分析表明,轉(zhuǎn)PtRCA基因南林895楊具耐高溫特性,利用CO2和強光進(jìn)行光合作用的能力較強,催化RuBP進(jìn)行羧化反應(yīng)的效率較高。轉(zhuǎn)PtRCA基因楊樹能夠充分利用吸收的光子,PSⅡ反應(yīng)中心效率較高,過剩光能量得到較好的耗散,表現(xiàn)出耐光氧化的能力。研究結(jié)果表明,PtRCA基因的高效表達(dá)提高了轉(zhuǎn)基因植株的光合效率,并對高溫高光強具有調(diào)節(jié)能力。

    楊樹; Rubisco活化酶基因; 轉(zhuǎn)基因; 光合作用

    光合作用是植物學(xué)中最活躍的研究領(lǐng)域,也是地球上最重要的化學(xué)反應(yīng)之一。Calvin循環(huán)是高等植物進(jìn)行光合作用碳固定的最基本途徑,其反應(yīng)共有13步,由RuBP(核酮糖-1,5-二磷酸)開始到RuBP再生結(jié)束,所有反應(yīng)均在葉綠體基質(zhì)中進(jìn)行(Rainesetal., 2006; Smithetal., 2007)。Calvin循環(huán)的第1步反應(yīng)由核酮糖-1,5二磷酸羧化酶(EC 4.1.1.39, 簡稱Rubisco)催化。Rubisco是光合碳固定的關(guān)鍵酶,其催化RuBP羧化反應(yīng)和加氧反應(yīng),并進(jìn)一步調(diào)節(jié)2種反應(yīng)的關(guān)系(Jurczyketal., 2015)。在植物體內(nèi)環(huán)境中,Rubisco的活性受Rubisco活化酶(Rubisco activase, RCA)的調(diào)節(jié)與控制,即Rubisco必須經(jīng)過RCA的活化,才能表現(xiàn)出其羧化與加氧活性(Parryetal., 2012)。因此,RCA的發(fā)現(xiàn)為提高植物Rubisco活力,從而提高碳固定的效率提供了新的思路。因此,RCA的生理與分子生物學(xué)研究已經(jīng)成為光合作用領(lǐng)域的熱門課題(李海霞等, 2010)。

    隨著基因工程的發(fā)展,克隆、表達(dá)、轉(zhuǎn)化RCA基因等相關(guān)研究工作已經(jīng)逐步展開??茖W(xué)家們成功克隆擬南芥(Arabidopsisthaliana)、大麥(Hordeumvulgare)、玉米(Zeamays)、水稻(Oryzasativa)、棉花(Gossypiumspp.)等植物RCA基因的cDNA(Wernekeetal., 1989; Rundleetal., 1991; Toetal., 1999; Salvuccietal., 2003)。研究發(fā)現(xiàn),RCA基因?qū)Χ唐诟邷孛{迫反應(yīng)較為敏感,但其表達(dá)存在轉(zhuǎn)錄后調(diào)控,這種機制對于提高植物光合作用,從而應(yīng)對長期高溫脅迫具有重要意義(DeRidderetal., 2007)。光照強度對RCA活性影響的相關(guān)研究表明: RCA活性會隨著光量子通量密度的增加而隨之增加,而光強達(dá)到300 μmol·m-2s-1時迅速趨于飽和狀態(tài); 然而Rubisco活性要在1 000 μmol·m-2s-1光照強度時才趨于飽和(姜振升等, 2010)。對水稻進(jìn)行高溫脅迫的試驗表明,RCA對植物抗高溫有調(diào)節(jié)作用,主要由于面臨高溫脅迫時,植物體內(nèi)的RCA能夠與熱激蛋白等分子伴侶蛋白互作,從而保證高溫處理下RCA的活性(Yamorietal., 2012)。轉(zhuǎn)RCA基因植株在熱脅迫時仍具有較高的光合速率(Bayramovetal., 2014)。另有報道,RCA也是一種能夠與赤霉素相結(jié)合的磷蛋白,對GA信號的傳遞有積極作用(Stotzetal., 2011)。RCA與植物光合作用速率(Tsaietal., 2015)和植物產(chǎn)量(Yinetal., 2014)具有密切關(guān)系。這些研究為使用基因工程技術(shù)改良作物及經(jīng)濟林木,提高其光合速率和抗逆性提供了新的思路。

    楊樹(Populus)是一種適應(yīng)能力強、品種繁多、生長速度快的優(yōu)良樹種,可作為人造板、造紙原料等,基因工程技術(shù)對提高楊樹抗逆性提供了新途徑。本研究選擇楊樹為研究材料,利用楊樹基因組序列信息,克隆毛果楊(Populustrichocarpa)PtRCA基因,通過農(nóng)桿菌(Agrobacterium)介導(dǎo)遺傳轉(zhuǎn)化南林895楊(P.deltoides×P.euramericana‘Nanlin895’),獲得PtRCA高效表達(dá)植株,并進(jìn)行分子檢測及相關(guān)功能分析,為培育楊樹新型高光效抗逆品種提供依據(jù)。

    1 材料與方法

    1.1 植物材料

    毛果楊及南林895楊無菌組培苗均于溫度22 ℃、光強4 000 lx、日光照16 h的溫室內(nèi)培養(yǎng),分別用于目的基因克隆以及相關(guān)功能分析。

    1.2 毛果楊PtRCA基因的克隆及分析

    采用BIOMIGA公司的EZgeneTM Plant Easy Spin RNA Miniprep Kit (R6611)試劑盒提取毛果楊總RNA; 使用TAKARA公司提供的試劑反轉(zhuǎn)錄cDNA。根據(jù)已有的PtRCA基因序列的保守區(qū),設(shè)計5′端及3′端特異引物(表1),擴增PtRCA基因的編碼區(qū),并進(jìn)行D-TOPO PCR反應(yīng)。將所得PCR產(chǎn)物進(jìn)行瓊脂糖凝膠電泳,切膠回收,連接到PMD 19載體上,轉(zhuǎn)化大腸桿菌(Escherichiacoli)TOP 10,經(jīng)過藍(lán)白斑篩選得到陽性克隆,挑選單克隆,交金斯瑞公司測序,并進(jìn)行相關(guān)序列比對,利用生物信息學(xué)軟件對PtRCA基因進(jìn)行基因功能及其結(jié)構(gòu)分析(表2)。

    表1 克隆擴增引物及RT-PCR表達(dá)分析引物Tab. 1 Primers used in full length cloning and RT-PCR analysis

    表2 基因序列分析軟件以及相關(guān)網(wǎng)址Tab. 2 Gene sequence analysis softwares and related web sites

    1.3 溫度脅迫下RCA基因表達(dá)模式分析

    用于溫度脅迫試驗的植物材料為南林895楊溫室生長苗(22 ℃),6周正常生長后進(jìn)行脅迫處理。第1種處理: 脅迫溫度梯度設(shè)置為10, 20,30 ℃,分別在脅迫1,3,5,7,9天后取葉片,進(jìn)行相關(guān)定量分析; 第2種處理: 37 ℃高溫條件下,分別在0,2,4,8,12 h時取葉片; 再將溫度調(diào)回室溫22 ℃,分別在2,4,8,12 h時取葉片,對所有葉片樣本進(jìn)行相關(guān)定量分析。每次試驗相互獨立,分別重復(fù)3次。

    1.4 楊樹PtRCA基因表達(dá)載體構(gòu)建

    PMD 19 T質(zhì)粒; 用于入門載體擴增的TOP 10大腸桿菌菌株; 用于表達(dá)載體擴增的DB 3.1大腸桿菌菌株; 表達(dá)載體pGWB 406(日本RIKEN理化研究所贈予)。

    按照Gateway反應(yīng)程序,設(shè)計BP反應(yīng)引物,F(xiàn)端: GGGGACAAGTTTGTACAAAAAAGCAAGGCTCC ATGGCAGCCACCATCTCCAC; R端: GGGGACCAC TTTGTACAAGAAAGCTGGGTCTAACCATAGAAAGTT CCTC。PCR 產(chǎn)物回收純化后進(jìn)行BP 反應(yīng)和LR 反應(yīng)。將LR反應(yīng)產(chǎn)物轉(zhuǎn)化到TOP 10感受態(tài)細(xì)胞,進(jìn)行菌落PCR驗證并測序,提取質(zhì)粒,得到目的基因植物表達(dá)載體,-20 ℃保存。

    1.5 楊樹PtRCA基因的遺傳轉(zhuǎn)化與分子檢測

    將含有sGFP基因的表達(dá)載體轉(zhuǎn)入農(nóng)桿菌EHA 105,通過葉盤法進(jìn)行遺傳轉(zhuǎn)化。侵染濃度OD600值設(shè)置為0.6,預(yù)培養(yǎng)時間3天,重懸時間60 min。

    將預(yù)培養(yǎng)后的南林895楊外植體浸入重懸液中輕晃振蕩15 min,用無菌濾紙吸去葉片上多余的菌液,平鋪于分化培養(yǎng)基(MS + 6-BA 0.5 mg·L-1+TDZ 0.002 mg·L-1+蔗糖 30 g·L-1+瓊脂 6.5 g·L-1)上,待分化篩選培養(yǎng)基(分化培養(yǎng)基 + Kan 40 mg·L-1+ 頭孢 200 mg·L-1+ 特美汀 200 mg·L-1)上長出狀態(tài)較好的抗性芽后,轉(zhuǎn)至芽伸長篩選培養(yǎng)基(芽伸長培養(yǎng)基 + Kan 20 mg·L-1+ 頭孢 100 mg·L-1+ 特美汀 100 mg·L-1),最后轉(zhuǎn)移至生根篩選培養(yǎng)基(生根培養(yǎng)基 + Kan 10 mg·L-1+ 頭孢 100 mg·L-1+ 特美汀 100 mg·L-1)上,于22 ℃培養(yǎng)溫度、4 000 lx光強、每天16 h光照的條件下溫室培養(yǎng),最后獲得完整抗性植株。

    利用PCR對初步篩選出來的抗性植株進(jìn)行再次篩選,以陽性質(zhì)粒為正對照,未轉(zhuǎn)基因南林895楊為負(fù)對照進(jìn)行分子檢測。

    1.6 光合速率相關(guān)參數(shù)的測定

    在2016年6月采用GFS-3000(便攜式氣體交換熒光測定儀)測定轉(zhuǎn)PtRCA基因南林895楊和對照的光合速率相關(guān)參數(shù)。選擇與對照(3株)相同時間移栽、長勢良好且無病蟲害的盆栽轉(zhuǎn)基因植株30株(每個株系3株)作為樣本,測定時采用樹冠中部外圍著生于枝條頂端的3~5對生長健壯的葉片。每個株系共測3片,每個處理重復(fù)5次。單次測定在15 min內(nèi)完成。

    1.6.1 光合速率(Pn)日變化 晴朗無云條件下,在7: 00—19: 00之間,每隔1 h測定1次光合速率值,選用植株頂端的3對葉片,每株測定3片,重復(fù)記錄5次。

    1.6.2 光合作用對光響應(yīng) 晴朗無云條件下,光合有效輻射在0~1 400 μmol·m-2s-1范圍內(nèi)設(shè)定若干梯度,測定光合速率值。測定時,葉室溫度30 ℃,葉室CO2濃度350 μmol·mol-1(與大氣中CO2濃度相同),相對濕度60%; 相關(guān)參數(shù)通過回歸分析獲得。

    1.6.3 光合作用對CO2響應(yīng) 晴朗無云條件下,CO2濃度在0~1 000 μmol·mol-1范圍內(nèi)設(shè)定若干梯度,測定光合速率值。測定時,葉室溫度30 ℃,光合有效輻射1 200 μmol·m-2s-1,相對濕度60%; 相關(guān)參數(shù)通過回歸分析獲得。

    1.7 葉綠素?zé)晒鈪?shù)的測定

    在2016年8月采用GFS-3000(便攜式氣體交換熒光測定儀)測定葉綠素?zé)晒鈪?shù),選擇與對照(長勢最好的1株)相同時間移栽、長勢良好且無病蟲害的轉(zhuǎn)PtRCA基因南林895楊10株(每個株系中長勢最好的1株)作為樣本,測定時采用樹冠中部外圍著生于枝條頂端的第3對生長健壯的葉片,每個數(shù)值均重復(fù)測5次。

    在光照強度300 μmol·m-2s-1條件下, 葉片暗適應(yīng)20 min后再將光照強度增至1 800 μmol·m-2s-1,在2種光強條件下測定PSⅡ原初光化學(xué)效率(Fv/Fm)、PSⅡ的實際光化學(xué)效率(Fv′/Fm′)、光化學(xué)猝滅(qP)、非光化學(xué)猝滅(qN)等葉綠素?zé)晒鈪?shù)。1.8 數(shù)據(jù)處理

    采用Excel軟件進(jìn)行數(shù)據(jù)處理和圖表繪制,采用SAS統(tǒng)計軟件的Duncan法進(jìn)行處理間的差異顯著性分析。

    2 結(jié)果與分析

    2.1 楊樹PtRCA基因序列分析

    根據(jù)已有的PtRCA基因序列的保守區(qū),設(shè)計5′端及3′端特異引物,進(jìn)行D-TOPO PCR轉(zhuǎn)化,通過藍(lán)白斑篩選挑選白色陽性克隆,經(jīng)過菌落PCR,所有菌樣結(jié)果均出現(xiàn)目的條帶; 將測序結(jié)果通過NCBI中BLAST,與克隆的毛果楊PtRCA蛋白序列進(jìn)行對比,發(fā)現(xiàn)克隆片段與楊樹預(yù)測片段相一致,具完整閱讀框。將PtRCA基因序列在GenBank中登錄,獲得登錄號為KY307834。

    通過3次獨立的測序,PtRCA基因編碼序列(CDS)長1 323 bp,編碼區(qū)包含440個氨基酸,利用ProtScale對PtRCA基因編碼區(qū)CDS結(jié)構(gòu)域分析,發(fā)現(xiàn)楊樹RCA蛋白N-末端存在1個轉(zhuǎn)運肽結(jié)合位點,含有AAA+家族蛋白保守區(qū)P-loop NTPase 超家族結(jié)構(gòu)域,以及2個ATP結(jié)合域、1個ATPase響應(yīng)結(jié)合域(圖1)。

    圖1 PtRCA基因編碼區(qū)推測的核苷酸序列以及氨基酸序列Fig.1 Nucleic acid sequence and amino acid sequence of PtRCA粗實線區(qū)表示N-末端; 箭頭表示轉(zhuǎn)運肽結(jié)合位點; 陰影區(qū)表示AAA+家族蛋白保守區(qū); 雙實線且加框表示ATP結(jié)合域; 加框表示ATPase sensor 2 motif; 左邊數(shù)字表示核苷酸數(shù)字。Heavy line area refers to N-terminal; The arrow refers to transit peptide binding sites; The shadow area refers to the conservative district of AAA+; Double solid line and framed area refers to ATP domain; Framed area refers to ATPase sensor 2 motif; The number refers to the number of nucleotides.

    利用ExPASy ProtParam在線分析PtRCA的理化特性: PtRCA蛋白分子式為C2135H3367N579O656S21,相對分子質(zhì)量為48 315.9 Da。共包括440個氨基酸,主要氨基酸為Asp,Glu,Arg,Lys,正電荷殘基數(shù)為49,負(fù)電荷殘基數(shù)為53,等電點pI為5.57。 不穩(wěn)定指數(shù)為33.4,屬于穩(wěn)定性蛋白。脂肪系數(shù)為77.55,平均親水系數(shù)為-0.311,為疏水性蛋白。

    2.2 楊樹PtRCA基因蛋白同源性及系統(tǒng)進(jìn)化分析

    利用MAGA 5.0軟件結(jié)合NCBI同源蛋白序列搜索結(jié)果構(gòu)建進(jìn)化樹以分析楊樹PtRCA基因與其他物種RCA基因親緣關(guān)系,選取的物種分別為擬南芥、大豆(Glycinemax)、葡萄(Vitisvinifera)、菜豆(Phaseolusvulgaris)、小籽海紅豆(Adenantheramicrosperma)、桃(Amygdaluspersica)、薺菜(Capsellabursa-pastoris)、麻瘋樹(Jatrophacurcas)、南芥(Arabis)、水田碎米薺(Cardaminelyrata),楊樹PtRCA蛋白與擬南芥、麻瘋樹、大豆等具較高同源性,應(yīng)具備與其相似功能,這為楊樹PtRCA蛋白功能分析奠定了良好基礎(chǔ)。同時,將楊樹PtRCA基因編碼的蛋白與楊樹中其他蛋白質(zhì)序列進(jìn)行同源比對,并分析與PtRCA基因編碼的蛋白同源性較高的蛋白,發(fā)現(xiàn)其均屬于AAA+超級蛋白家族。這些都對PtRCA基因表達(dá)模式分析提供了參考依據(jù)。

    2.3 楊樹RCA基因溫度梯度脅迫誘導(dǎo)表達(dá)

    利用qRT-PCR檢測南林895楊(非轉(zhuǎn)基因)的RCA基因溫度梯度脅迫誘導(dǎo)表達(dá)情況,結(jié)果(圖2)顯示: 10 ℃低溫處理1天后,RCA基因表達(dá)量下降; 處理3天表達(dá)量極顯著提高,隨后顯著下降; 處理9天,表達(dá)量也顯著提高。20 ℃處理,RCA基因表達(dá)量基本呈現(xiàn)穩(wěn)定上升趨勢。30 ℃高溫處理1天后,表達(dá)量顯著升高,3天后顯著下降,5天后表達(dá)量逐步上升。具較高耐熱性的RCA轉(zhuǎn)基因植株,受到中度熱脅迫時,仍具有較高的光合速率以及生產(chǎn)量,故提高植物RCA熱穩(wěn)定性可作為改善植物抵御高溫脅迫的一種可行性方案(Kureketal., 2007)。

    圖2 溫度脅迫下RCA基因表達(dá)模式分析Fig.2 Real-time PCR results of RCA expression under diverse temperature stresses

    實時定量分析(圖3)顯示高溫(37 ℃)處理2 h后,RCA基因表達(dá)量顯著上調(diào),隨后逐步降低,在高溫處理12 h后,表達(dá)量達(dá)到最低,甚至于不表達(dá)。而在經(jīng)過22 ℃室溫恢復(fù)處理后,RCA基因表達(dá)量開始呈線性上升趨勢,且在恢復(fù)處理12 h時,RCA表達(dá)量超過了未經(jīng)高溫處理的初始表達(dá)量。這說明適度高溫處理可誘導(dǎo)RCA表達(dá)(Wangetal., 2010),這為推測RCA具耐高溫能力提供了依據(jù)。

    圖3 37 ℃高溫誘導(dǎo)RCA基因表達(dá)模式分析Fig.3 Real-time PCR results of the expression of RCA induced by 37 ℃*,P < 0.05; **,P<0.01.

    2.4 載體的構(gòu)建

    目的基因先經(jīng)過BP反應(yīng)整合到入門載體,再經(jīng)過LR反應(yīng),整合到表達(dá)載體pGWB 406上,通過PCR凝膠電泳,結(jié)果顯示目的基因成功構(gòu)入表達(dá)載體。

    2.5 轉(zhuǎn)PtRCA基因南林895楊的獲得

    待分化篩選培養(yǎng)基上長出狀態(tài)較好的抗性芽后,需及時轉(zhuǎn)移至芽伸長篩選培養(yǎng)基上,于22 ℃培養(yǎng)溫度、4 000 lx光強、每天16 h光照的條件下溫室培養(yǎng)。每隔10~15天,更換1次培養(yǎng)基,并實時觀察其生長狀況。

    待抗性芽生長至3 cm高時,需及時轉(zhuǎn)移至生根篩選培養(yǎng)基上,于22 ℃培養(yǎng)溫度、4 000 lx光強、每天16 h光照的條件下溫室培養(yǎng),最后獲得轉(zhuǎn)PtRCA基因南林895楊(圖4)。

    圖4 轉(zhuǎn)PtRCA基因南林895楊的獲得Fig.4 Obtaining of PtRCA transgenic ‘Nanlin895’ poplarA: 葉盤分化; B: 抗性芽伸長; C: 抗性芽的壯大; D: 抗性植株。A: Leaf differentiation; B: Resistant shoot elongation; C. Resistant shoot growing; D. Transgenic plant.

    利用PCR分子檢測對初步篩選出來的抗性植株進(jìn)行再次篩選,以陽性質(zhì)粒為正對照,未轉(zhuǎn)化南林895楊為負(fù)對照,顯示部分結(jié)果(圖5),PCR擴增出的特異條帶與陽性對照一致,本試驗共篩選出14個株系為陽性。

    圖5 轉(zhuǎn)PtRCA基因南林895楊PCR分子檢測電泳結(jié)果Fig.5 Electrophoresis results for PCR molecular detection of PtRCA transgenic ‘Nanlin895’ poplarM: DL 2000 marker. 1-4: 轉(zhuǎn)基因楊樹(sGFP表達(dá)載體片段); -CK: 野生型; +CK: 陽性質(zhì)粒。1-4: Transgenic lines; -CK: Negative control; +CK: Positive control.

    2.6 高溫脅迫下轉(zhuǎn)PtRCA基因南林895楊與對照RCA基因表達(dá)比較分析

    以室溫(22 ℃)正常生長且未轉(zhuǎn)化植株(南林895楊)組培苗為對照組,對生長45天的轉(zhuǎn)PtRCA基因南林895楊和未轉(zhuǎn)化植株(南林895楊)的組培苗進(jìn)行高溫脅迫處理,每升高5 ℃設(shè)置1個處理溫度,37 ℃封頂,設(shè)置的3個高溫脅迫溫度為27,32,37 ℃,分別共同處理1,3,5天后,取組培苗頂端第3片葉進(jìn)行采樣分析。

    實時定量分析結(jié)果(圖6)顯示: 高溫處理1天后,27 ℃轉(zhuǎn)基因植株RCA基因表達(dá)量顯著上升,32,37 ℃處理植株RCA基因表達(dá)量均有所下降,未轉(zhuǎn)化植株RCA基因表達(dá)量為極顯著下降; 高溫處理3天后,所有處理植株RCA基因表達(dá)量均顯著上升,27 ℃處理的植株RCA基因表達(dá)量為極顯著上升,且轉(zhuǎn)基因植株RCA基因表達(dá)量是未轉(zhuǎn)化植株RCA基因表達(dá)量的2倍左右; 高溫處理5天后,所有處理植株RCA基因表達(dá)量均有所下降,其中27 ℃處理的轉(zhuǎn)基因植株RCA基因表達(dá)量極顯著高于對照組,32 ℃處理的轉(zhuǎn)基因植株RCA基因表達(dá)量也高于對照組,而37 ℃處理的轉(zhuǎn)基因植株RCA基因表達(dá)量與對照組持平。此外,任一高溫和時間處理下的轉(zhuǎn)基因植株RCA基因表達(dá)量均顯著高于未轉(zhuǎn)化植株RCA基因表達(dá)量。

    圖6 高溫脅迫下南林895楊RCA基因表達(dá)模式分析Fig.6 Real-time PCR results of RCA expression in ‘Nanlin895’ poplar under high temperature stress*,P < 0.05; **,P<0.01.

    2.7 光合速率相關(guān)參數(shù)的測定

    2.7.1 轉(zhuǎn)PtRCA基因楊樹與對照的光合速率日變化規(guī)律比較 對光合速率日變化測定數(shù)據(jù)的分析表明(選取轉(zhuǎn)基因株系RCA1、RCA2、RCA6和對照進(jìn)行對比分析): 1)日變化曲線形態(tài)的不同: 對照植株的峰值出現(xiàn)在11: 00和14: 00,日變化呈雙峰曲線,最高凈光合速率(Pn)為12.7 μmol·m-2s-1左右,其低谷出現(xiàn)在12: 00—13: 00之間; 而3個轉(zhuǎn)PtRCA基因株系的平均峰值出現(xiàn)在12: 00,日變化基本呈單峰曲線,最高凈光合速率(Pn)為14.9 μmol·m-2s-1(RCA6),平均持續(xù)1 h左右開始逐漸下降(圖7)。2)與對照相比,轉(zhuǎn)PtRCA基因楊樹在7: 00的凈光合速率(Pn)并無明顯差異,為3.9~4.3 μmol·m-2s-1之間。

    根據(jù)上述分析,推測可能是由于PtRCA基因?qū)耄纳屏宿D(zhuǎn)PtRCA基因楊樹的C3光合循環(huán),并提高了轉(zhuǎn)基因植株利用強光的能力,使其充分利用強光進(jìn)行光合作用; 相比之下,對照則出現(xiàn)明顯的光抑制現(xiàn)象。

    圖7 轉(zhuǎn)PtRCA基因南林895楊和對照凈光合速率(Pn)的日變化Fig.7 Diurnal change of net photosynthetic rate(Pn) between transgenic poplars and controlRCA1,RCA2,RCA6: 轉(zhuǎn)基因株系Transgenic lines.下同。The same below.

    2.7.2 轉(zhuǎn)PtRCA基因楊樹與對照的光合作用對光響應(yīng)曲線的比較 光合作用對光的響應(yīng)測定數(shù)據(jù)(選取轉(zhuǎn)基因株系RCA1、RCA2、RCA6和對照進(jìn)行對比分析)表明: 光照強度0~400 μmol·m-2s-1時,轉(zhuǎn)PtRCA基因楊樹與對照的凈光合速率均呈線性增長,并隨著光照強度的不斷增加逐漸達(dá)到飽和狀態(tài)(圖8)。因此,依據(jù)凈光合速率與0~400 μmol·m-2s-1范圍內(nèi)光照強度呈線性關(guān)系,建立10個轉(zhuǎn)基因株系和對照的直線回歸方程,計算出光飽和點、光補償點、表觀量子效率、暗呼吸速率以及光飽和時的光合速率(表3)。結(jié)果表明: 1)光飽和點: 轉(zhuǎn)PtRCA基因楊樹所有株系的光飽和點均比對照高,增加約12.5%~37.5%,其中8號株系的光飽和點最高。這表現(xiàn)出轉(zhuǎn)PtRCA基因楊樹能夠更好地利用強光進(jìn)行光合作用。2)光補償點: 轉(zhuǎn)PtRCA基因楊樹(除3號株系外)的光補償點均略低于對照楊樹。這表明轉(zhuǎn)PtRCA基因楊樹具有較強的利用弱光能力。3)表觀量子效率: 轉(zhuǎn)PtRCA基因楊樹(除2號和6號株系外)的表觀量子效率均略高于對照楊樹。這表明轉(zhuǎn)PtRCA基因楊樹的光合作用更加活躍。4)暗呼吸速率: 轉(zhuǎn)PtRCA基因楊樹(除3號和9號株系外)的暗呼吸速率均略低于對照楊樹。這表明轉(zhuǎn)PtRCA基因楊樹在沒有光照時消耗的能量偏少。5)光飽和時的光合速率: 轉(zhuǎn)PtRCA基因楊樹在光飽和點時的光合速率比對照增高24.6%~55.7%(1號和2號株系的光合速率最高)。這進(jìn)一步表明高光強能夠提高轉(zhuǎn)基因植株的羧化反應(yīng)能力,使光合速率得到提高。

    2.7.3 轉(zhuǎn)PtRCA基因楊樹與對照的光合作用對CO2濃度響應(yīng)曲線的比較 光合作用對CO2濃度響應(yīng)測定數(shù)據(jù)(選取轉(zhuǎn)基因株系RCA1、RCA2、RCA6和對照進(jìn)行對比分析)表明: CO2濃度0~400 μmol·mol-1時,轉(zhuǎn)PtRCA基因楊樹和對照的凈光合速率呈線性增長,并隨著CO2濃度的不斷增加逐漸達(dá)到飽和狀態(tài)(圖9)。因此,依據(jù)凈光合速率與0~400 μmol·mol-1范圍內(nèi)CO2濃度呈線性關(guān)系,建立10個轉(zhuǎn)基因株系和對照的直線回歸方程,計算出CO2飽和點、CO2補償點、羧化效率、光呼吸速率、CO2飽和時的光合速率(表4)。結(jié)果表明: 1)CO2飽和點: 轉(zhuǎn)PtRCA基因楊樹CO2飽和點除1號和4號株系外,其他均比對照低12.5%~25.0%,在略高于大氣CO2濃度下就達(dá)到飽和。2)CO2補償點: 轉(zhuǎn)PtRCA基因楊樹的CO2補償點比對照大幅度降低,低53.1%~80.4%。這表明由于PtRCA基因的導(dǎo)入,轉(zhuǎn)基因植株對CO2親和力較高,這也是轉(zhuǎn)基因植株CO2飽和點低的原因。3)羧化效率: 轉(zhuǎn)PtRCA基因楊樹的羧化效率均略高于對照。這表明轉(zhuǎn)PtRCA基因楊樹的羧化效率較強,能更好地催化RuBP進(jìn)行羧化反應(yīng)。4)光呼吸速率: 轉(zhuǎn)PtRCA基因楊樹的光呼吸比對照低37.7%~79.3%。這表明轉(zhuǎn)PtRCA基因楊樹有利于光合作用有機物積累,光呼吸消耗較少,其生物量增加。5)CO2飽和時的光合速率: 轉(zhuǎn)PtRCA基因楊樹在CO2飽和點時的光合速率比對照增加4.4%~26.4%(2號株系的光合速率最高)。轉(zhuǎn)基因植株在CO2飽和點較低的情況下,光合速率卻比對照高,這也進(jìn)一步說明轉(zhuǎn)基因植株有較強的羧化反應(yīng)能力,光合速率能夠得到提高。

    圖8 轉(zhuǎn)PtRCA基因楊樹和對照的光合響應(yīng)變化Fig.8 Changes of photosynthetic response between transgenic poplars and the control設(shè)定條件: 葉室溫度30 ℃,葉室CO2濃度350 μmol·mol-1,相對濕度60%。Testing conditions: Hamuro temperature of 30 ℃, Hamuro CO2 concentration of 350 μmol·mol-1, relative humidity of 60%.

    表3 轉(zhuǎn)PtRCA基因楊樹和對照的光合響應(yīng)參數(shù)指標(biāo)Tab. 3 Index of photosynthetic response parameters of transgenic poplar and the control

    2.8 葉綠素?zé)晒鈪?shù)的比較

    2.8.1 轉(zhuǎn)PtRCA基因楊樹與對照的PSⅡ原初光化學(xué)效率(Fv/Fm)和實際光化學(xué)效率(Fv′/Fm′)比較 由表5的數(shù)據(jù)進(jìn)行對比分析可知: 1)PSⅡ原初光化學(xué)效率(Fv/Fm): 在低光照強度下,與對照相比,轉(zhuǎn)PtRCA基因楊樹的Fv/Fm沒有明顯差異; 而經(jīng)過光氧化處理后,F(xiàn)v/Fm均存在不同程度的降低; 相比較而言,對照下降得最多,為61.7%,轉(zhuǎn)基因楊樹下降45.0%~53.1%(4號株系下降得最少,為45%)。2)PSⅡ?qū)嶋H光化學(xué)效率(Fv′/Fm′): 在低光照強度下,轉(zhuǎn)PtRCA基因楊樹和對照的Fv′/Fm′沒有明顯差異; 然而經(jīng)過光氧化處理后,二者的Fv′/Fm′均有不同程度的降低,對照下降得最多,為54.1%,而轉(zhuǎn)PtRCA基因楊樹下降38.7%~52.0%(6號株系下降得最少,為38.7%)。這表明轉(zhuǎn)PtRCA基因植株利用吸收的光子的能力較強,吸收的光子供給PSⅡ反應(yīng)中心的效率較高,推測這可能是由于PtRCA基因的導(dǎo)入,提高植株P(guān)SⅡ反應(yīng)中心的效率。

    圖9 轉(zhuǎn)PtRCA基因楊樹和對照在不同CO2濃度下的光合響應(yīng)Fig.9 Photosynthetic response to CO2 concentration between transgenic poplars and the control設(shè)定條件: 葉室溫度30 ℃,光照強度1 200 μmol·m-2s-1,相對濕度60%。Testing conditions: Hamuro temperature of 30 ℃, light intensity of 1 200 μmol·m-2s-1, relative humidity of 60%.

    表4 轉(zhuǎn)PtRCA基因楊樹和對照的CO2光合響應(yīng)參數(shù)指標(biāo)Tab.4 Index of photosynthetic response to CO2 concentration between transgenic poplars and control

    2.8.2 轉(zhuǎn)PtRCA基因楊樹與對照的光化學(xué)猝滅(qP)和非光化學(xué)猝滅(qN)比較 由表6的數(shù)據(jù)進(jìn)行對比分析可知: 1)光化學(xué)猝滅系數(shù)(qP): 在低光照強度下,轉(zhuǎn)PtRCA基因楊樹與對照的qP沒有明顯差異; 經(jīng)過光氧化處理后,二者的qP均下降,對照下降了68.3%,而轉(zhuǎn)基因楊樹下降較少,為51.0%~65.8%(2號株系下降的最少,為51.0%)。這表明在強光下,與對照相比,轉(zhuǎn)PtRCA基因楊樹的光化學(xué)電子傳遞比較順暢,用于光化學(xué)電子傳遞的光能份額比較高,而對照出現(xiàn)了明顯的受阻現(xiàn)象。2)非光化學(xué)猝滅系數(shù)(qN): 在低光照強度下,轉(zhuǎn)PtRCA基因楊樹與對照的qN沒有明顯差異; 經(jīng)過光氧化處理后,對照增加僅為3.0%,而轉(zhuǎn)PtRCA基因楊樹增加較多,為6.0%~26.5%(2號株系增加的最多,為26.5%)。這表明PtRCA基因的導(dǎo)入,使得轉(zhuǎn)PtRCA基因楊樹的光合機構(gòu)沒有受到破壞,過剩光能得到較好的耗散,表現(xiàn)出較好的耐光氧化能力。

    表5 轉(zhuǎn)PtRCA基因楊樹和對照的PSⅡ原初光化學(xué)效率(Fv/Fm)和實際光化學(xué)效率(Fv′/Fm′)在不同光強下的比較①Tab.5 PS Ⅱ primary photochemical efficiency (Fv/Fm) and PSⅡ actual photochemical efficiency (Fv′/ Fm′) in different light intensity between transgenic poplars and control

    ① 平均值±標(biāo)準(zhǔn)偏差(5組數(shù)據(jù))。Mean ±SD (n=5).

    表6 轉(zhuǎn)PtRCA基因楊樹和對照的光化學(xué)猝滅(qP)和非光化學(xué)猝滅(qN)在不同光強下的比較①Tab.6 Photochemical quenching (qP) and non-photochemical quenching (qN) in different light intensity between transgenic poplars and the control

    ① 平均值±標(biāo)準(zhǔn)偏差(5組數(shù)據(jù))。Mean±SD (n=5).

    3 討論

    高光強一直是抑制植物生長、降低光合速率的主要因素,然而,近幾年一系列研究表明: Rubisco活化酶(RCA)具抵御高光強脅迫的作用,這為培育高產(chǎn)量抗逆植株提供了理論依據(jù)。

    3.1 轉(zhuǎn)PtRCA基因南林895楊與其他轉(zhuǎn)RCA基因植物的比較分析

    在研究擬南芥Rubisco活化問題時首先發(fā)現(xiàn)RCA,經(jīng)過對RCA進(jìn)行一系列抗體研究,發(fā)現(xiàn)它存在于高等植物及一些單細(xì)胞光合生物中,這說明在光合生物中,RCA對Rubisco活性調(diào)節(jié)機制普遍存在,RCA與光合作用有著密切關(guān)系(Hasseetal., 2015)。

    利用農(nóng)桿菌遺傳介導(dǎo)獲得了RCA基因正義擬南芥,結(jié)果顯示僅表達(dá)大亞基的正義植株各項指標(biāo)與野生型并無明顯差異,而只表達(dá)小亞基的正義轉(zhuǎn)化苗長勢明顯優(yōu)于野生型(Zhangetal., 2002)。同樣,相似正義遺傳轉(zhuǎn)化試驗結(jié)果顯示,轉(zhuǎn)小亞基正義植株長勢顯著優(yōu)于野生型植株,另外,該研究發(fā)現(xiàn)RCA與GA信號傳導(dǎo)有密切聯(lián)系(Sharmaetal., 2002; Wangetal., 2010)。雖然以上試驗均表明RCA基因主要是由小亞基表達(dá)功能,但對獲得的正義水稻研究表明,轉(zhuǎn)大亞基植株的Rubisco活性、光合效率指標(biāo)均明顯高于野生型,作物產(chǎn)量也有所提高(Wuetal., 2007)。這些研究為通過基因工程,增強植物的抗逆性,進(jìn)而提高植物光合速率和產(chǎn)量提供了新的思路。研究轉(zhuǎn)基因水稻RCA基因調(diào)控時,發(fā)現(xiàn)RCA對葉綠體的發(fā)育有調(diào)節(jié)作用,且在水稻發(fā)育生長后期,轉(zhuǎn)基因植株葉片中RCA含量一直偏高,隨著葉片的衰老,RCA含量逐步降低,光合速率下降,說明植物光合速率受RCA影響,RCA對凈光合速率起重要調(diào)控作用(蔣德安等, 2000; 翁曉燕等, 2001; 張國等, 2005)。

    本研究嘗試?yán)棉D(zhuǎn)基因技術(shù)使毛果楊PtRCA基因在楊樹中高效表達(dá),通過PtRCA基因的特性,使Rubisco處于穩(wěn)定的催化活性狀態(tài),讓Rubisco更加有效地參與到光合碳同化過程中。本研究結(jié)果也表明,克隆得到的PtRCA基因蛋白與擬南芥等RCA蛋白同源性較高,可作為研究楊樹PtRCA蛋白功能的參考依據(jù)。分析證實了RCA具耐高溫耐熱能力, 轉(zhuǎn)PtRCA基因楊樹能夠充分利用中午的強光進(jìn)行光合作用,利用CO2的能力較強,能更好地催化RuBP進(jìn)行羧化反應(yīng),羧化效率較強。在強光下,轉(zhuǎn)PtRCA基因楊樹的光合機構(gòu)沒有受到破壞,過剩光能得到較好的耗散,表現(xiàn)出一定的耐光氧化能力。

    上述結(jié)果表明具較高耐高溫耐熱性的RCA轉(zhuǎn)基因植株,受到中度熱脅迫時,仍具有較高的光合速率以及生產(chǎn)量。這與在研究擬南芥轉(zhuǎn)基因植株時,發(fā)現(xiàn)RCA對凈光合速率熱恢復(fù)具重要作用(Kureketal., 2007)等相關(guān)研究結(jié)果一致。這為利用轉(zhuǎn)基因手段進(jìn)行楊樹新型高光效抗逆新品種的研究提供了重要依據(jù)。

    3.2 轉(zhuǎn)PtRCA基因南林895楊與轉(zhuǎn)PEPC基因(C4植物)南林895楊的比較分析

    與C3植物只利用卡爾文循環(huán)中1,5-二磷酸核酮糖直接固定CO2相比,C4植物首先利用磷酸烯醇式丙酮酸(PEP),經(jīng)PEP羧化酶(PEPC)的作用與CO2結(jié)合,形成蘋果酸或天門冬氨酸。這些四碳雙羧酸轉(zhuǎn)移到鞘細(xì)胞里,通過脫羧酶的作用釋放CO2,后者在鞘細(xì)胞葉綠體內(nèi)經(jīng)核酮糖1,5-二磷酸(RuBP)羧化酶作用,進(jìn)入卡爾文循環(huán)。因為PEP與CO2的結(jié)合能力遠(yuǎn)遠(yuǎn)大于RuBP,這使得C4植物即使在CO2濃度比較低的時候,依然可以通過這一點源源不斷為光合作用提供充足的CO2,使得C4植物光合作用也較強。因此,通過分子技術(shù)讓C3植物與C4植物一樣,具有更加高效的光合作用已成為現(xiàn)今研究的熱門課題(Leegood, 2013)。

    據(jù)報導(dǎo),水生植物Egeriadensa和Hydrillaverticilata體內(nèi)也發(fā)現(xiàn)有C4光合途徑,這一發(fā)現(xiàn)很大程度上預(yù)示著,或許可以通過C4光合酶基因向C3植物的遺傳轉(zhuǎn)化從而提高C3植物光合作用的可操作性(Tengetal., 2001; Lietal., 2005)。作者利用農(nóng)桿菌介導(dǎo)法,將光合作用C4途徑關(guān)鍵酶基因PEPC導(dǎo)入楊樹并進(jìn)行相關(guān)特性的分析,其研究結(jié)果(尹吳等, 2012)與本研究相比: 轉(zhuǎn)PEPC基因楊樹具有C4植物的某些高光效特性,在光合速率、光能利用效率、羧化效率等方面的表現(xiàn)與其對照相比具有一定的優(yōu)勢,這些表現(xiàn)在一定程度上優(yōu)于本試驗中轉(zhuǎn)PtRCA基因株系對應(yīng)的數(shù)據(jù); 然而,在耐光氧化方面,轉(zhuǎn)PtRCA基因楊樹比轉(zhuǎn)PEPC基因楊樹具有更加明顯的優(yōu)勢。上述對比說明PEPC基因能夠更好地固定CO2,對CO2的親合力較強,從而促使磷酸烯醇式丙酮酸(PEP)固定大氣中濃度較低的CO2,供C3途徑利用,因此轉(zhuǎn)PEPC基因楊樹的羧化效率更高,利用CO2的能力更強,而RCA基因具有很強的耐高溫抗高光強的特性。

    4 結(jié)論

    轉(zhuǎn)PtRCA基因楊樹具有較高的凈光合速率,推測是由于PtRCA基因的導(dǎo)入,提高了轉(zhuǎn)PtRCA基因楊樹利用強光的能力,使其能夠充分利用中午的強光進(jìn)行光合作用。轉(zhuǎn)PtRCA基因楊樹所有株系的光飽和點較高,且在光飽和點時的光合速率較高。這說明轉(zhuǎn)PtRCA基因楊樹能夠較好地利用強光進(jìn)行光合作用。轉(zhuǎn)PtRCA基因楊樹CO2飽和點較低(除1號和4號株系外); 而在CO2飽和點時的光合速率,轉(zhuǎn)PtRCA基因楊樹較高; 此外,轉(zhuǎn)PtRCA基因楊樹的羧化效率較高。這些均表明轉(zhuǎn)PtRCA基因楊樹利用CO2的能力較強,能夠更好地催化RuBP進(jìn)行羧化反應(yīng),羧化效率較強。

    轉(zhuǎn)PtRCA基因楊樹利用吸收的光子的能力較強,吸收的光子供給PSⅡ反應(yīng)中心的效率較高,推測這可能是由于PtRCA基因的導(dǎo)入提高了PSⅡ反應(yīng)中心的效率。在高光照強度下,與未轉(zhuǎn)基因?qū)φ障啾?,轉(zhuǎn)PtRCA基因楊樹的光化學(xué)電子傳遞比較順暢,其PSⅡ天線色素吸收的光能用于光化學(xué)電子傳遞的份額較高,且轉(zhuǎn)PtRCA基因楊樹表現(xiàn)出耐光氧化的能力,其過剩光能量得到較好的耗散。

    蔣德安, 陸 慶, 翁曉燕, 等. 2000. 水稻劍葉衰老期Rubisco活化酶對Rubisco活力和光合速率的調(diào)節(jié). 浙江大學(xué)學(xué)報: 農(nóng)業(yè)與生命科學(xué)版, 26 (2): 119-124.

    (Jiang D A, Lu Q, Weng X Y,etal. 2000. Regulation of Rubisco carboxylation activity and photosynthetic rate by Rubisco activase during leaf senescence in rice. Journal of Zhejiang University: Agriculture and Life Sciences, 26(2): 119-124.[in Chinese])

    姜振升, 孫曉琦, 艾希珍, 等.2010.低溫弱光對黃瓜幼苗Rubisco與Rubisco 活化酶的影響.應(yīng)用生態(tài)學(xué)報,21(8): 2045-2050.

    (Jiang Z S, Sun X Q, Ai X Z,etal.2010.Responses of Rubisco and Rubisco activase in cucumber seedlings to low temperature and weak light. Chinese Journal of Applied Ecology,21(8): 2045-2050.[in Chinese])

    李海霞, 王真梅, 曾漢來. 2010. 植物Rubisco 活化酶的研究進(jìn)展.植物生理學(xué)報, (11): 1092-1100.

    (Li H X,Wang Z M,Zeng H L.2010. The research progresses in Rubisco activase in plant.Plant Physiology Journal,(11): 1092-1100.[in Chinese])

    翁曉燕, 陸 慶, 蔣德安. 2001. 水稻Rubisco活化酶在調(diào)節(jié)Rubisco活性和光合日變化中的作用. 中國水稻科學(xué), 15 (1): 35-40.

    (Weng X Y, Lu Q, Jiang D A. 2001. Rubisco activase and its regulation on diurnal changes of photosynthetic rate and the activity of ribulose 1,5-bisphosphate carboxyase/oxygenase (Rubisco). Chinese Journal of Rice Science, 15 (1): 35-40.[in Chinese])

    尹 吳, 李麗莎, 王立科, 等. 2012. 轉(zhuǎn)玉米PEPC基因楊樹的光合生理特性分析.林業(yè)科學(xué), 48(6): 63-71.

    (Yin W, Li L S, Wang L K,etal.2012. Analysis of photosynthetic characteristics of transgenic poplars with maizePEPCgene. Scientia Silvae Sinicae, 48(6): 63-71.[in Chinese])

    張 國, 李 濱, 鄒 琦. 2005. 小麥Rubisco活化酶基因的克隆和表達(dá)特性. 植物學(xué)通報, 22(3): 313-319.

    (Zhang G, Li B, Zou Q. 2005. Cloning and expression of Rubisco activase gene in wheat. Chinese Bulletin of Botany, 22 (3): 313-319.[in Chinese])

    Bayramov S, Guliyev N. 2014. Changes in Rubisco activase gene expression and polypeptide content inBrachypodiumdistachyon. Plant Physiol Biochem, 81: 61-66.

    DeRidder B P, Salvucci M E. 2007. Modulation of Rubisco activase gene expression during heat stress in cotton (GossypiumhirsutumL.) involves post-transcriptional mechanisms. Plant Science, 172: 246-254.

    Hasse D, Larsson A M, Andersson I. 2015. Structure ofArabidopsisthalianaRubisco activase. Acta Crystallographica, 71(4): 800-808.

    Jurczyk B, Hura K, Trzemecka A,etal. 2015. Evidence for alternative splicing mechanisms in meadow fescue(Festucapratensis) and perennial ryegrass (Loliumperenne) Rubisco activase gene. J Plant Physiol, 176: 61-64.

    Kurek I, Chang T K, Bertain S M,etal. 2007. Enhanced thermostability ofArabidopsisRubisco activase improves photosynthesis and growth rates under moderate heat stress. Plant Cell, 19(10): 3230-3241.

    Leegood R C. 2013. Strategies for engineering C4photosynthesis. J Plant Physiol, 170(4): 378-388.

    Li X, Jiao D. 2005. Transgenic rice overexpressing C4photosynthetic genes and their application in breeding. Mol Plant Breed, 3 (4): 550-556.

    Mahmoud S S, Wang S L, Moloney M M,etal.1998. Production of a biologically active novel goldfish growth hormone inEscherichiacoli. Comp Biochem Physiol B, 120:657-663.

    Motohashi K, Kondoh A, Stumpp M T,etal. 2001. Comprehensive survey of proteins targeted by chloroplast thioredoxin. Proc Natl Acad Sci USA, 98(20): 11224-11229.

    Parry M A, Andralojc P J, Scales J C,etal. 2012. Rubisco activity and regulation as targets for crop improvement. J Exp Bot, 64(3): 717-730.

    Raines C A, Paul M J. 2006. Products of leaf primary carbon metabolism modulate the developmental programme determining plant morphology. Journal of Experimental Botany, 57: 1857-1862.

    Rundle S J, Zielinski R E. 1991. Organization and expression of two tandemly oriented genes encoding ribulose bisphosphate carboxylase/oxygenase activase in barley. Journal of Biological Chemistry, 266: 4677-4685.

    Salvucci M E, van de Loo F J, Stecher D. 2003. Two isoforms of Rubisco activase in cotton, the products of separate genes not alternative splicing. Planta, 216: 736-744.

    Sharma A, Komatsu S. 2002. Involvement of a Ca2+-dependent protein kinase component downstream to the gibberellin-binding phosphoprotein, Rubisco activase, in rice. Biochemical and Biophysical Resarch Communications, 290: 690-695.

    Smith A M, Stitt M. 2007. Coordination of carbon supply and plant growth. Plant Cell Environment, 30: 1126-1149.

    Stotz M, Mueller O, Ciniawsky S,etal. 2011. Structure of green-type Rubisco activase from tobacco.Nat Strut Mol Biol, 18(12):1366-1370.

    Teng S, Qian Q, Huang D N. 2001. Advance in molecular biology and genetic engineering of C4photosynthesis pathway. J Agric Biotechnol, 9(2): 198-201.

    To K Y, Suen D F, Grace-Chen S C. 1999. Molecular characterization of ribulose-l,5-bisphosphate carboxylase/oxygenase activase in rice leaves. Planta, 209(1): 66-76.

    Tsai Y C C, Lapina M C, Bhushan S,etal. 2015. Identification and characterization of multiple rubisco activases in chemoautotrophic bacteria. Nat Commun, 6(3):8883.

    Wang D, Li X F, Zhou Z J,etal. 2010. Two Rubisco activase isoforms may play different roles in photosynthetic heat acclimation in the rice plant. Physiol Plantarum, 139(1): 55-67.

    Werneke J M, Chatfield J M, Ogren W L. 1989. Alternative mRNA splicing generates the two ribulose bisphosphate carboxylase/oxygenase activase polypeptides in spinach and arabidopsis. Plant Cell, 1(8): 815-825.

    Wu H R, Li L B, Jing Y X,etal. 2007. Over- and anti-sense expressions of the large isoform of ribulose-1, 5-bisphosphate carboxylase/oxygenase activase gene inOryzasativaaffect the photosynthetic capacity. Photosynthetica, 45 (2): 194-201.

    Yamori W, Masumoto C, Fukayama H,etal. 2012. Rubisco activase is a key regulator of non-steady-state photosynthesis at any leaf temperature and, to a lesser extent, of steady-state photosynthesis at high temperature. Plant J, 71(6): 871-880.

    Yin Z T, Zhang Z L, Deng D X,etal. 2014. Characterization of Rubisco activase genes in maize: an α-isoform gene functions alongside a β-isoform gene. Plant Physiol, 164(4): 2096-2106.

    Zhang N, Kallis R P, Ewy R G,etal. 2002. Light modulation of Rubisco inArabidopsisrequires a capacity for redox regulation of the larger Rubisco activase isoform. Proceeding of the National Academy of Sciences of the USA, 99(5): 3330-3334.

    (責(zé)任編輯 徐 紅)

    Cloning and Functional Analysis of Rubisco Activase Gene fromPopulustrichocarpa

    Yin Wu Sun Weibo Zhou Yan Zhuge Qiang

    (KeyLaboratoryofForestGenetics&BiotechnologyofMinistryofEducationCo-InnovationCenterforSustainableForestryinSouthernChinaNanjingForestryUniversityNanjing210037)

    【Objective】 Ribulose-1,5-bishosphate (Rubisco) is the key enzymes in the first step of carbon assimilation involved in plant photosynthesis. Rubisco was kept in a catalytic active state under the control of Rubisco activase(RCA), as a result, the efficiency of photosynthesis will be improved. In this research, thePtRCAgene of poplars was cloned. Transgenic poplars with high expression ofPtRCAwere obtained through this research. As a result, a new type of high efficiency photosynthesis with stress tolerance was proved to be feasible through the analysis of molecular detection and functional analysis. 【Method】According to the sequence ofPtRCAgene cloned fromPopulustrichocarpa, the analysis of function and structure ofPtRCAwere conducted through Bioinformatics software.PtRCAwas constructed into the expression vector pGWB406 by usingAgrobacteriummediated method. ‘Nanlin895’poplar plants (P.deltoides×P.euramericana‘Nanlin895’) were used in this research. Gene expression of both transgenic poplars and ‘Nanlin895’plants were compared under high temperature stress. The parameters of photosynthesis and chlorophyll fluorescence were compared between them as well. 【Result】The CDS sequence ofPtRCAwas 1 323 bp. There were 440 amino acid residues, and protein molecular weight is 48 315.9 Da. The isoelectric point is 5.57, which is hydrophobic protein without signal peptide and membrane structure.PtRCAis proved to be in the same AAA+ family according to the result of sequence comparison. It has same RCA homology protein with soybean and arabidopsis as well. The expression quantity of transgenic poplars were higher than ‘Nanlin895’ in average. Transgenic poplars showed advantage of photosynthesis during the noon, the light saturation point was 12.5% to 37.5% higher. Moreover, the efficiency of photosynthesis was 24.6% to 55.7% higher than the plants used as control. However, the light compensation point is lower than that of plants without transgenosis expect the No.3 strain. Besides, the usage of CO2and the carboxylation efficiency of the transgenic poplars were better. Expect No.1 and No.4 strains, CO2saturation point was around 12.5%-25.0% less than control plants. CO2compensation point and light respiration were proved to be 53.1%-80.4% and 37.7%-79.3% lower compared to the controls respectively. Moreover, the efficiency of photosynthesis at the CO2saturation point was around 4.4%-26.4% higher. The tolerance of transgenic poplars to light oxidation was increased according to this research. Under the light oxidation treatment, the primary photochemical efficiency of PSⅡ decreased 61.7% in control lines. And transgenic poplar decreased by 45.0%- 53.1%.The actual photochemical efficiency of PSⅡ in controls decreased 54.1%. And transgenic poplar decreased by 38.7%-52.0%.Photochemical quenching coefficient dropped 68.3%, and 51.0%-65.8% declined in transgenic poplars. Non photochemical quenching coefficient increased by 3.0% compared to the control, and 6.0%-26.5% increased in transgenic poplars.【Conclusion】 Poplar Rubisco activase (PtRCA) protein and soy protein,ArabidopsisthalianaRCA homology is higher. According to the real-time quantitative PCR and related physiological analysis, transgenic poplars (P.deltoides×P.euramericana‘Nanlin895’) showed high-temperature resistance, higher efficiency of photosynthesis, increased ability to use CO2and better catalyzing effect on carboxylation reaction. The efficiency of assimilating photons and supplying to PSⅡ was highly developed, and as a result, the light oxidation ability of transgenic poplars was proved to be enhanced. The efficiency of photosynthesis was improved by high expression ofPtRACgene in transgenic poplars, which also showed ability to adjust for high temperature and strong light.

    Populus; Rubisco activase gene; transgenic; photosynthesis

    10.11707/j.1001-7488.20170410

    2016-09-20;

    2016-12-10。

    國家科技部“863計劃”課題(2013AA102703); 國家國際科技合作專項(2014DFG32440); 國家自然科學(xué)基金項目(31570650); 江蘇省高校優(yōu)勢學(xué)科建設(shè)工程項目(PAPD)。

    S718.46; S718.43

    A

    1001-7488(2017)04-0083-13

    *諸葛強為通訊作者。

    猜你喜歡
    飽和點光化學(xué)株系
    過表達(dá)NtMYB4a基因增強煙草抗旱能力
    光化學(xué)蒸汽發(fā)生法在分析化學(xué)實驗教學(xué)中的應(yīng)用
    云南化工(2021年9期)2021-12-21 07:44:10
    安順山藥光合生理特性研究
    相似材料極限密度及抗壓強度穩(wěn)定性分析
    嫦娥5號返回式試驗衛(wèi)星小麥育種材料研究進(jìn)展情況
    衢州椪柑變異株系—黃皮椪柑相關(guān)特性研究
    浙江柑橘(2016年1期)2016-03-11 20:12:31
    對一道課后練習(xí)題的商榷
    中國化學(xué)會第29屆學(xué)術(shù)年會——來自光化學(xué)分會場的報道
    BiFeO3空間選擇性光化學(xué)還原Ag及光催化活性
    乳糜血對亞甲藍(lán)光化學(xué)法制備病毒滅活血漿的影響
    亚洲精品成人久久久久久| 国产成人91sexporn| av在线天堂中文字幕| 免费电影在线观看免费观看| 男女下面进入的视频免费午夜| 亚洲国产日韩欧美精品在线观看| 网址你懂的国产日韩在线| 亚洲成色77777| 国产av码专区亚洲av| 国产三级在线视频| 午夜日本视频在线| 欧美精品一区二区大全| 免费无遮挡裸体视频| 国产一区有黄有色的免费视频 | 99热网站在线观看| 2018国产大陆天天弄谢| 国产麻豆成人av免费视频| 久久午夜福利片| 欧美zozozo另类| 在线 av 中文字幕| 色尼玛亚洲综合影院| 国产午夜精品一二区理论片| 久久久久国产网址| 国产精品一区二区性色av| 成人亚洲欧美一区二区av| 国产一区二区三区综合在线观看 | 91午夜精品亚洲一区二区三区| 汤姆久久久久久久影院中文字幕 | 午夜福利在线在线| 成年版毛片免费区| 十八禁国产超污无遮挡网站| 一区二区三区乱码不卡18| 在线观看美女被高潮喷水网站| 免费av毛片视频| 又大又黄又爽视频免费| 国产 一区精品| 白带黄色成豆腐渣| 97超碰精品成人国产| 久久久久久久久中文| 精品酒店卫生间| 精品酒店卫生间| 亚洲av日韩在线播放| 身体一侧抽搐| 欧美最新免费一区二区三区| 成年女人看的毛片在线观看| av在线观看视频网站免费| 最近最新中文字幕免费大全7| 一级毛片电影观看| 免费观看无遮挡的男女| 欧美xxxx性猛交bbbb| 免费在线观看成人毛片| 国产麻豆成人av免费视频| 乱系列少妇在线播放| 亚洲精品第二区| 亚洲精品中文字幕在线视频 | 一区二区三区四区激情视频| 777米奇影视久久| 亚洲av.av天堂| 亚洲国产欧美人成| 一夜夜www| 国产精品嫩草影院av在线观看| 又爽又黄无遮挡网站| 国产精品一二三区在线看| 亚洲精品成人av观看孕妇| 欧美成人a在线观看| 日日摸夜夜添夜夜爱| 久久人人爽人人片av| 天堂网av新在线| 亚洲av日韩在线播放| 又粗又硬又长又爽又黄的视频| 日韩亚洲欧美综合| 久久精品熟女亚洲av麻豆精品 | 久久韩国三级中文字幕| 天堂影院成人在线观看| 全区人妻精品视频| 国产亚洲91精品色在线| 亚洲国产精品成人久久小说| 亚洲丝袜综合中文字幕| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 免费黄色在线免费观看| 久久久久久久久中文| 只有这里有精品99| 久久99热这里只频精品6学生| 久热久热在线精品观看| 亚洲经典国产精华液单| 美女黄网站色视频| 中文乱码字字幕精品一区二区三区 | 联通29元200g的流量卡| 日韩av在线免费看完整版不卡| 国产视频首页在线观看| 嘟嘟电影网在线观看| www.av在线官网国产| 精品久久久久久电影网| 国内精品宾馆在线| 午夜视频国产福利| 国产精品无大码| 91在线精品国自产拍蜜月| 成人亚洲精品av一区二区| 一个人看的www免费观看视频| 日韩制服骚丝袜av| 在现免费观看毛片| 小蜜桃在线观看免费完整版高清| 街头女战士在线观看网站| 一个人免费在线观看电影| 99久国产av精品| 看十八女毛片水多多多| 日韩大片免费观看网站| 久久久久精品性色| 中国美白少妇内射xxxbb| 91aial.com中文字幕在线观看| 日本av手机在线免费观看| 免费看日本二区| 看十八女毛片水多多多| 亚洲18禁久久av| 日韩 亚洲 欧美在线| 中文资源天堂在线| 久久精品久久久久久噜噜老黄| 好男人视频免费观看在线| 男女那种视频在线观看| 亚洲精品乱码久久久久久按摩| 日韩三级伦理在线观看| av免费在线看不卡| 夫妻午夜视频| 成人亚洲欧美一区二区av| 免费高清在线观看视频在线观看| 中国国产av一级| 亚洲国产精品成人综合色| 欧美一区二区亚洲| 少妇裸体淫交视频免费看高清| 午夜免费激情av| 国产69精品久久久久777片| 久久精品国产亚洲av涩爱| 国产一区二区在线观看日韩| 2021少妇久久久久久久久久久| 啦啦啦韩国在线观看视频| 午夜福利在线在线| 校园人妻丝袜中文字幕| 久久精品国产亚洲av天美| 在线天堂最新版资源| 床上黄色一级片| 成年免费大片在线观看| 亚洲av二区三区四区| 777米奇影视久久| 一级黄片播放器| 少妇熟女aⅴ在线视频| 色综合亚洲欧美另类图片| 国产淫片久久久久久久久| 在线播放无遮挡| 少妇人妻一区二区三区视频| 男插女下体视频免费在线播放| 亚洲丝袜综合中文字幕| av国产免费在线观看| 18禁在线无遮挡免费观看视频| 国产色婷婷99| 亚洲成人精品中文字幕电影| 一级毛片我不卡| 久久精品熟女亚洲av麻豆精品 | 久久久久精品性色| 日韩国内少妇激情av| 日韩制服骚丝袜av| 国产成人aa在线观看| 十八禁网站网址无遮挡 | 日韩电影二区| 97人妻精品一区二区三区麻豆| 精品久久久精品久久久| 国产欧美日韩精品一区二区| 少妇被粗大猛烈的视频| 在现免费观看毛片| 国产成人freesex在线| 亚洲av男天堂| 三级国产精品片| .国产精品久久| 麻豆国产97在线/欧美| 婷婷色av中文字幕| 在线天堂最新版资源| 国产精品爽爽va在线观看网站| 亚洲精品国产av蜜桃| 欧美xxⅹ黑人| 国产综合精华液| 极品少妇高潮喷水抽搐| 久久99精品国语久久久| 亚洲精品国产av成人精品| kizo精华| 欧美日韩精品成人综合77777| 一级毛片 在线播放| 欧美日韩视频高清一区二区三区二| 一本久久精品| 99久久人妻综合| 精品少妇黑人巨大在线播放| 伊人久久精品亚洲午夜| 精品不卡国产一区二区三区| 看非洲黑人一级黄片| 我要看日韩黄色一级片| 亚洲欧美清纯卡通| 免费黄网站久久成人精品| 免费看不卡的av| 国产精品1区2区在线观看.| 国产一区二区三区av在线| 丰满人妻一区二区三区视频av| 久99久视频精品免费| 久久精品国产亚洲av涩爱| 国产欧美另类精品又又久久亚洲欧美| 婷婷色av中文字幕| 成人毛片60女人毛片免费| 校园人妻丝袜中文字幕| 国产探花在线观看一区二区| 精品久久久久久久久亚洲| 久久久色成人| 女人久久www免费人成看片| 亚洲18禁久久av| 嘟嘟电影网在线观看| 寂寞人妻少妇视频99o| 免费看美女性在线毛片视频| 久久久久久久亚洲中文字幕| 亚洲精品成人久久久久久| 老司机影院成人| 国产成人freesex在线| 亚洲自偷自拍三级| 午夜福利成人在线免费观看| 国产永久视频网站| 一区二区三区免费毛片| 人妻一区二区av| 免费不卡的大黄色大毛片视频在线观看 | 在线播放无遮挡| 男女下面进入的视频免费午夜| 久久精品综合一区二区三区| 成人亚洲欧美一区二区av| 免费观看av网站的网址| 欧美日韩精品成人综合77777| 国产黄色小视频在线观看| 直男gayav资源| 插逼视频在线观看| 亚洲精品日韩av片在线观看| 午夜福利视频精品| 精品一区二区三卡| 国产av在哪里看| 少妇猛男粗大的猛烈进出视频 | 麻豆精品久久久久久蜜桃| 成人午夜精彩视频在线观看| 亚洲av中文av极速乱| 国产精品精品国产色婷婷| 九九爱精品视频在线观看| av在线蜜桃| 日韩不卡一区二区三区视频在线| 丰满少妇做爰视频| 又黄又爽又刺激的免费视频.| 男女啪啪激烈高潮av片| 69av精品久久久久久| 晚上一个人看的免费电影| 成年女人在线观看亚洲视频 | 亚洲色图av天堂| 国产不卡一卡二| 色播亚洲综合网| 男女那种视频在线观看| 2022亚洲国产成人精品| or卡值多少钱| 波多野结衣巨乳人妻| 国产一区二区三区综合在线观看 | 春色校园在线视频观看| 亚洲在线自拍视频| 日韩国内少妇激情av| 22中文网久久字幕| 国产成人免费观看mmmm| 一区二区三区四区激情视频| 女的被弄到高潮叫床怎么办| 久久久久久久久久成人| 久久99精品国语久久久| 校园人妻丝袜中文字幕| 中文在线观看免费www的网站| 身体一侧抽搐| 丰满人妻一区二区三区视频av| 高清欧美精品videossex| 国产精品麻豆人妻色哟哟久久 | 国产精品一区www在线观看| 欧美高清成人免费视频www| 深夜a级毛片| 国产精品久久视频播放| 天堂网av新在线| 三级国产精品欧美在线观看| 一级二级三级毛片免费看| 国产单亲对白刺激| 亚洲性久久影院| 亚洲成人精品中文字幕电影| 我要看日韩黄色一级片| 国产69精品久久久久777片| 赤兔流量卡办理| 婷婷色麻豆天堂久久| 九色成人免费人妻av| 欧美日韩国产mv在线观看视频 | 国产伦精品一区二区三区视频9| 男女国产视频网站| 国产亚洲av嫩草精品影院| 国产精品蜜桃在线观看| 国产又色又爽无遮挡免| 人妻夜夜爽99麻豆av| 国产精品爽爽va在线观看网站| 天堂俺去俺来也www色官网 | 三级男女做爰猛烈吃奶摸视频| 国产 一区精品| 日韩制服骚丝袜av| 99久久精品国产国产毛片| 亚洲国产色片| 超碰97精品在线观看| 免费看不卡的av| 国产亚洲5aaaaa淫片| 免费av毛片视频| 免费观看无遮挡的男女| 韩国高清视频一区二区三区| 国产淫语在线视频| 精品酒店卫生间| 久久精品人妻少妇| 黑人高潮一二区| 最后的刺客免费高清国语| 欧美丝袜亚洲另类| 亚洲精品一区蜜桃| 国产激情偷乱视频一区二区| 女人久久www免费人成看片| 亚洲熟女精品中文字幕| 欧美激情在线99| 欧美日韩亚洲高清精品| 三级经典国产精品| 免费av毛片视频| 免费看av在线观看网站| 欧美区成人在线视频| 亚洲真实伦在线观看| 听说在线观看完整版免费高清| 美女高潮的动态| 我的老师免费观看完整版| 视频中文字幕在线观看| 亚洲久久久久久中文字幕| 三级经典国产精品| 亚洲精品色激情综合| 免费黄色在线免费观看| 日韩,欧美,国产一区二区三区| 国产精品三级大全| 秋霞在线观看毛片| 亚洲色图av天堂| 亚洲av福利一区| 亚洲精品第二区| 99视频精品全部免费 在线| 超碰av人人做人人爽久久| 熟妇人妻不卡中文字幕| 亚洲国产高清在线一区二区三| 一级片'在线观看视频| 亚洲人与动物交配视频| 看免费成人av毛片| 26uuu在线亚洲综合色| 毛片女人毛片| 少妇的逼水好多| 97人妻精品一区二区三区麻豆| 亚洲综合精品二区| 看黄色毛片网站| 菩萨蛮人人尽说江南好唐韦庄| 97热精品久久久久久| 国产色爽女视频免费观看| 午夜福利在线观看免费完整高清在| 国产精品人妻久久久影院| 99久国产av精品| 在线a可以看的网站| 国产成人免费观看mmmm| 99久国产av精品国产电影| 99热这里只有是精品50| 精品人妻偷拍中文字幕| 免费大片18禁| 欧美高清性xxxxhd video| 肉色欧美久久久久久久蜜桃 | www.色视频.com| 老司机影院成人| 91精品一卡2卡3卡4卡| 久久国产乱子免费精品| 好男人视频免费观看在线| 国产高清三级在线| 精品人妻偷拍中文字幕| 综合色丁香网| 高清在线视频一区二区三区| 日本-黄色视频高清免费观看| 久久久久久国产a免费观看| 国产亚洲一区二区精品| 欧美激情久久久久久爽电影| 日韩大片免费观看网站| 国产精品久久久久久精品电影小说 | 久久国内精品自在自线图片| 日本一本二区三区精品| 国产永久视频网站| 美女xxoo啪啪120秒动态图| 成年av动漫网址| 欧美xxⅹ黑人| 舔av片在线| 丝瓜视频免费看黄片| 国产午夜福利久久久久久| 夫妻午夜视频| 如何舔出高潮| 成年女人看的毛片在线观看| 精品一区在线观看国产| 中文字幕制服av| 天堂中文最新版在线下载 | 街头女战士在线观看网站| 欧美日韩国产mv在线观看视频 | 天美传媒精品一区二区| 亚洲精品一区蜜桃| 久久久亚洲精品成人影院| 久久6这里有精品| 日韩视频在线欧美| 国产午夜精品久久久久久一区二区三区| 69av精品久久久久久| 欧美97在线视频| 欧美日韩在线观看h| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线观看播放| 少妇被粗大猛烈的视频| 成年人午夜在线观看视频 | 中文欧美无线码| 一二三四中文在线观看免费高清| 国产伦在线观看视频一区| 少妇的逼好多水| av在线老鸭窝| 99久久中文字幕三级久久日本| 中文字幕亚洲精品专区| 国产一区二区三区av在线| 五月天丁香电影| 欧美最新免费一区二区三区| 成人漫画全彩无遮挡| av专区在线播放| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av天美| 26uuu在线亚洲综合色| 一级av片app| 日日啪夜夜撸| 春色校园在线视频观看| 久久鲁丝午夜福利片| 久久久精品欧美日韩精品| 国产在视频线在精品| 91av网一区二区| 中文乱码字字幕精品一区二区三区 | 国产熟女欧美一区二区| 观看美女的网站| 黄色欧美视频在线观看| 国产精品国产三级专区第一集| 精品久久久精品久久久| 国内精品一区二区在线观看| 又粗又硬又长又爽又黄的视频| 少妇熟女欧美另类| 欧美不卡视频在线免费观看| 免费电影在线观看免费观看| 在线免费十八禁| 内射极品少妇av片p| 免费电影在线观看免费观看| 午夜激情福利司机影院| 国产精品.久久久| 在线免费观看的www视频| 麻豆乱淫一区二区| 波野结衣二区三区在线| 色综合站精品国产| 亚洲av不卡在线观看| 高清午夜精品一区二区三区| 日韩强制内射视频| 精品亚洲乱码少妇综合久久| 久久久久久九九精品二区国产| 嫩草影院新地址| 日韩欧美国产在线观看| 亚洲av中文av极速乱| 性色avwww在线观看| 老司机影院毛片| 天天躁日日操中文字幕| 欧美激情久久久久久爽电影| 久久久精品欧美日韩精品| 免费少妇av软件| 亚州av有码| 国产综合懂色| 日本wwww免费看| 亚洲人成网站在线播| av播播在线观看一区| a级毛色黄片| 亚洲激情五月婷婷啪啪| 久久久久精品性色| 99视频精品全部免费 在线| 亚洲高清免费不卡视频| 少妇高潮的动态图| 久久97久久精品| 精品少妇黑人巨大在线播放| 国产伦精品一区二区三区四那| 男人舔奶头视频| 狂野欧美白嫩少妇大欣赏| 久久久亚洲精品成人影院| 简卡轻食公司| 亚洲欧美精品自产自拍| 男插女下体视频免费在线播放| 插逼视频在线观看| 国产精品久久视频播放| 精品不卡国产一区二区三区| 亚洲综合精品二区| 天堂中文最新版在线下载 | 久久99热这里只频精品6学生| 熟妇人妻久久中文字幕3abv| 一夜夜www| 国产毛片a区久久久久| 久久精品久久精品一区二区三区| 麻豆成人av视频| 国产精品三级大全| 精品久久久久久久人妻蜜臀av| 国产国拍精品亚洲av在线观看| 天堂av国产一区二区熟女人妻| 久久久久久久午夜电影| 中文资源天堂在线| 久久97久久精品| 国产高清不卡午夜福利| 一本一本综合久久| 一区二区三区高清视频在线| 韩国av在线不卡| 日本色播在线视频| 日本三级黄在线观看| 九九久久精品国产亚洲av麻豆| 国产一区有黄有色的免费视频 | 亚洲国产精品成人久久小说| 汤姆久久久久久久影院中文字幕 | 国产老妇女一区| 中文乱码字字幕精品一区二区三区 | 色综合站精品国产| 蜜桃久久精品国产亚洲av| 国内揄拍国产精品人妻在线| 日韩欧美国产在线观看| 少妇的逼好多水| 久久久久免费精品人妻一区二区| 嫩草影院新地址| 大陆偷拍与自拍| 午夜免费激情av| 日韩欧美 国产精品| 国产美女午夜福利| 人人妻人人澡欧美一区二区| 欧美 日韩 精品 国产| 国产伦精品一区二区三区视频9| 老师上课跳d突然被开到最大视频| 亚洲自偷自拍三级| 国产伦精品一区二区三区四那| 99九九线精品视频在线观看视频| 久久鲁丝午夜福利片| 午夜福利视频精品| 九九在线视频观看精品| 99久久精品国产国产毛片| 久久热精品热| 又爽又黄无遮挡网站| 亚洲av成人精品一二三区| 国产 一区 欧美 日韩| 亚洲熟妇中文字幕五十中出| 国产午夜精品一二区理论片| 一个人免费在线观看电影| 亚洲av二区三区四区| 日本欧美国产在线视频| 亚洲精品视频女| 亚洲精品第二区| av国产免费在线观看| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 亚洲国产成人一精品久久久| 国产一区二区在线观看日韩| 美女高潮的动态| 久久久久精品久久久久真实原创| 国产亚洲av片在线观看秒播厂 | 久久精品夜色国产| 蜜臀久久99精品久久宅男| 男人狂女人下面高潮的视频| 日韩伦理黄色片| 天堂影院成人在线观看| 99久久中文字幕三级久久日本| 久久精品国产鲁丝片午夜精品| 国产伦精品一区二区三区视频9| 久久精品国产鲁丝片午夜精品| 亚洲怡红院男人天堂| 久久久久久久国产电影| 少妇丰满av| av黄色大香蕉| 欧美另类一区| 99re6热这里在线精品视频| 18禁在线无遮挡免费观看视频| 亚洲电影在线观看av| 水蜜桃什么品种好| 2021天堂中文幕一二区在线观| 熟妇人妻不卡中文字幕| 国产欧美另类精品又又久久亚洲欧美| 波多野结衣巨乳人妻| 美女xxoo啪啪120秒动态图| 看免费成人av毛片| 成人无遮挡网站| 国产精品久久久久久久电影| 人妻少妇偷人精品九色| 国产亚洲5aaaaa淫片| 欧美成人午夜免费资源| 中文精品一卡2卡3卡4更新| 男女那种视频在线观看| 国产亚洲91精品色在线| 午夜久久久久精精品| 久久久久久久久久黄片| 小蜜桃在线观看免费完整版高清| 我要看日韩黄色一级片| 国产午夜精品论理片| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区在线观看99 | 亚洲一区高清亚洲精品| 久久99精品国语久久久| 欧美 日韩 精品 国产| av免费在线看不卡| 国产精品不卡视频一区二区| 国产av国产精品国产| 插逼视频在线观看| 国产精品熟女久久久久浪| 男女啪啪激烈高潮av片| 99久国产av精品| 国产成人免费观看mmmm| 秋霞伦理黄片| 一级毛片久久久久久久久女| 亚洲av二区三区四区| 成人av在线播放网站| 人人妻人人澡欧美一区二区| 别揉我奶头 嗯啊视频|