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

    杉木磷轉(zhuǎn)運(yùn)蛋白基因ClPht1;1的克隆及表達(dá)分析*

    2017-06-23 12:08:44蘇爍爍吳鵬飛馬祥慶
    林業(yè)科學(xué) 2017年5期
    關(guān)鍵詞:家系杉木利用效率

    蘇爍爍 李 明,2 吳鵬飛,2 張 穎 馬祥慶,2

    (1.福建農(nóng)林大學(xué)林學(xué)院 福州 350002; 2.國家林業(yè)局杉木工程技術(shù)研究中心 福州 350002)

    ?

    杉木磷轉(zhuǎn)運(yùn)蛋白基因ClPht1;1的克隆及表達(dá)分析*

    蘇爍爍1李 明1,2吳鵬飛1,2張 穎1馬祥慶1,2

    (1.福建農(nóng)林大學(xué)林學(xué)院 福州 350002; 2.國家林業(yè)局杉木工程技術(shù)研究中心 福州 350002)

    【目的】PHT1基因家族是影響植物吸收磷營養(yǎng)的重要磷轉(zhuǎn)運(yùn)子之一。從杉木32號(hào)磷高效家系cDNA中克隆得到PHT1基因家族的1個(gè)杉木磷轉(zhuǎn)運(yùn)蛋白基因ClPht1;1,并對(duì)不同程度磷脅迫下ClPht1;1的時(shí)空表達(dá)進(jìn)行研究,為杉木PHT1基因序列特征和功能結(jié)構(gòu)的研究以及磷高效利用杉木基因型的選育奠定基礎(chǔ)。【方法】 根據(jù)PHT1基因家族序列保守性設(shè)計(jì)簡并引物,以32號(hào)磷高效杉木基因型根系cDNA為模板進(jìn)行擴(kuò)增獲得目的基因ClPht1;1的cDNA序列,使用RACE 技術(shù)對(duì)目的基因進(jìn)行全長克隆,并對(duì)其序列特征、同源性和編碼磷轉(zhuǎn)運(yùn)蛋白結(jié)構(gòu)進(jìn)行分析。實(shí)時(shí)熒光定量PCR檢測ClPht1;1在32號(hào)磷高效杉木家系根、莖、葉中的表達(dá),檢測中度缺磷脅迫下ClPht1;1在不同磷利用效率杉木4號(hào)、15號(hào)、25號(hào)、27號(hào)、28號(hào)、32號(hào)家系根系中的表達(dá)差異,以及在中度、重度缺磷脅迫下ClPht1;1在32號(hào)磷高效杉木家系根系中隨時(shí)間序列的表達(dá)量變化?!窘Y(jié)果】 克隆得到1個(gè)杉木磷轉(zhuǎn)運(yùn)蛋白PHT1基因,命名為ClPht1;1(GenBank登錄號(hào): KX302006),基因序列編碼區(qū)長1 638 bp,編碼545 aa的蛋白質(zhì)。ClPht1;1所編碼蛋白質(zhì)由12個(gè)疏水的跨膜區(qū)域組成,1個(gè)疑似跨膜域。每個(gè)跨膜結(jié)構(gòu)域基本由17~25個(gè)氨基酸殘基組成螺旋,同時(shí)跨膜蛋白的N 端和C端均位于細(xì)胞質(zhì)內(nèi),保守序列位于第4個(gè)跨膜域。構(gòu)成蛋白質(zhì)的主要骨架是α-螺旋,無信號(hào)肽序列。ClPht1;1基因編碼蛋白與日本柳杉PHT基因編碼蛋白的氨基酸序列相似性達(dá)到87.0%,與胡楊、油茶、馬尾松等PHT家族基因編碼蛋白的氨基酸序列相似性均在75%以上。ClPht1;1基因在杉木的根、莖、葉組織中均有表達(dá),其中在根中的表達(dá)量最高,在葉中的表達(dá)量最低。在中度缺磷脅迫下,ClPht1;1基因在杉木不同家系根部的表達(dá)量為25號(hào)>27號(hào)>4號(hào)>15號(hào)>32號(hào)>28號(hào)。在中度和重度缺磷脅迫下,ClPht1;1基因在32號(hào)杉木家系根部的表達(dá)量隨脅迫時(shí)間的延長而逐漸上升; 恢復(fù)供磷后,ClPht1;1基因表達(dá)量逐漸下降至正常水平; 重度缺磷脅迫下,ClPht1;1基因表達(dá)量要高于其在中度缺磷脅迫下的表達(dá)量?!窘Y(jié)論】ClPht1;1基因具有PHT1基因家族的典型結(jié)構(gòu),其編碼蛋白的氨基酸序列與日本柳杉等磷轉(zhuǎn)運(yùn)蛋白氨基酸序列具有高度相似性,為杉木高親和磷轉(zhuǎn)運(yùn)蛋白PHT1基因家族成員。ClPht1;1基因主要在杉木的根部表達(dá),在葉片中的表達(dá)量較低; 杉木磷利用效率越強(qiáng),ClPht1;1基因在其根部的表達(dá)量越高。在不同磷利用效率的杉木家系中ClPht1;1基因表達(dá)量存在較大差異。ClPht1;1 基因的表達(dá)受低磷脅迫的誘導(dǎo),缺磷脅迫下ClPht1;1基因表達(dá)量明顯升高,恢復(fù)供磷后ClPht1;1基因表達(dá)量明顯降低。

    杉木; 磷轉(zhuǎn)運(yùn)蛋白基因;Pht1;1; RACE; RT-qPCR; 磷利用效率

    植物PHT1家族成員基因表達(dá)水平的升高是植物適應(yīng)低磷脅迫的普遍機(jī)制。對(duì)擬南芥(Arabidopsisthaliana)、小麥(Triticumaestivum)、水稻(Oryzasativa)、大豆(Glycinemax)等植物PHT1基因的表達(dá)研究表明: 在介質(zhì)磷濃度滿足植物生長情況下,PHT1基因表達(dá)受到抑制,少量表達(dá)或不表達(dá); 在介質(zhì)磷濃度較低時(shí),植物高親和力轉(zhuǎn)運(yùn)系統(tǒng)啟動(dòng),從而誘導(dǎo)PHT1基因強(qiáng)烈表達(dá); 對(duì)PHT1 基因表達(dá)進(jìn)行干擾將會(huì)導(dǎo)致根系吸收磷能力的降低(Yang etal., 2010)。因此,PHT1基因的高表達(dá)是植物磷高效吸收利用的重要指標(biāo)之一,挖掘植物磷高效吸收利用的基因資源進(jìn)行遺傳改良,是緩解土壤有效磷缺乏和提高磷肥利用效率的有效途徑之一。研究表明,植物PHT1家族的成員組成和表達(dá)受植物種及品種的影響,水稻、玉米(Zeamays)等植物不同品種在低磷脅迫下的PHT1 基因表達(dá)具有明顯差異(Aietal., 2009; Calderon-Vazquezetal., 2008)。在水稻、大豆、歐洲油菜(Brassicanapus)、煙草(Nicotianatabacum)、油茶(Camelliaoleifera)、蒺藜苜蓿(Medicagotruncatula)等農(nóng)作物中都克隆得到了部分PHT1 家族的磷轉(zhuǎn)運(yùn)蛋白基因成員,并根據(jù)PHT1基因的表達(dá)量變化來篩選和培育了一些磷高效農(nóng)作物品種或無性系,為揭示植物的磷轉(zhuǎn)運(yùn)機(jī)制和提高磷素利用效率奠定了基礎(chǔ)。目前,對(duì)植物PHT1基因的研究主要集中在1年生農(nóng)作物上,在森林樹種上的研究報(bào)道較少,僅見毛果楊(Populustrichocarpa)、橡膠樹(Heveabrasiliensis)等喬木中。

    杉木(Cunninghamialanceolata)是我國南方最主要的用材林造林樹種之一,其面積和蓄積量均居我國人工林樹種的首位(俞新妥, 2006)。由于南方酸性土壤對(duì)磷具有強(qiáng)烈化學(xué)固定作用,土壤磷素多以Fe-P、Al-P和O-P等形態(tài)存在(賈興永等, 2011),這使得南方林區(qū)土壤固定態(tài)磷豐富,但是有效磷極度缺乏,杉木等南方林木表現(xiàn)出“遺傳學(xué)缺磷”現(xiàn)象(梁霞等, 2006)。長期以來,南方森林土壤有效磷不足一直是限制杉木人工林產(chǎn)量最重要的因素之一。因此篩選和培育磷高效杉木基因型,提高杉木對(duì)土壤磷素的吸收和利用效率,是培育杉木速生高效人工林和緩解多代連栽下杉木人工林地力衰退的重要方向。然而目前在杉木磷轉(zhuǎn)運(yùn)蛋白基因方面研究報(bào)道較少,據(jù)此課題組在長期的杉木耐低磷脅迫試驗(yàn)研究中,選育了一些磷高效吸收利用杉木基因型,這些杉木基因型或主要通過根系皮層細(xì)胞溶解、加快體內(nèi)磷素的轉(zhuǎn)運(yùn)來抵抗磷脅迫逆境(4號(hào)、32號(hào)家系等); 或主要通過在貧磷斑塊中的根系向四周拓展增生和分泌化學(xué)物質(zhì),增強(qiáng)對(duì)介質(zhì)中磷的吸收效率來適應(yīng)環(huán)境的磷脅迫(25號(hào)、32號(hào)家系等)(汪攀, 2015; 吳鵬飛等, 2012)。以這些磷高效利用杉木基因型為研究材料,對(duì)于發(fā)掘杉木高親和磷轉(zhuǎn)運(yùn)子PHT1基因起到重要作用并具有實(shí)踐意義。

    本試驗(yàn)以磷高效利用杉木4號(hào)、15號(hào)、25號(hào)、27號(hào)、28號(hào)、32號(hào)家系為試驗(yàn)材料,根據(jù)植物PHT1基因的保守性氨基酸序列設(shè)計(jì)簡并引物,利用RT-PCR和RACE技術(shù)對(duì)杉木PHT1基因全長進(jìn)行克隆,對(duì)其序列特征、同源性和編碼磷轉(zhuǎn)運(yùn)蛋白結(jié)構(gòu)進(jìn)行分析,并使用熒光定量PCR檢測杉木PHT1基因在磷高效杉木基因型中的時(shí)空表達(dá)情況,以期為揭示杉木吸收磷營養(yǎng)的分子機(jī)制和選育磷高效利用杉木基因型奠定基礎(chǔ)。

    1 材料與方法

    1.1 試驗(yàn)材料

    選取種植在福建農(nóng)林大學(xué)妙峰山試驗(yàn)苗圃的1年生不同利用效率的磷高效杉木4號(hào)(1.14±0.08)kg·g-1、15號(hào)(1.10±0.06)kg·g-1、25號(hào)(1.77±0.22)kg·g-1、27號(hào)(1.38±0.03)kg·g-1、28號(hào)(0.73±0.06)kg·g-1、32號(hào)(1.09±0.02)kg·g-1(汪攀, 2015)家系幼苗,根系經(jīng)純水清洗后固定于塑料小桶進(jìn)行水培供磷試驗(yàn)。幼苗經(jīng)過7天正常供磷營養(yǎng)液培養(yǎng)后,進(jìn)行不同程度缺磷脅迫處理。分別在缺磷脅迫0,12,24,48 h,以及恢復(fù)正常供磷12,24,48 h時(shí)采集幼苗的根、莖和葉,液氮速凍后-80 ℃保存待用。

    利用KH2PO4設(shè)置3個(gè)供磷水平: 正常供磷(1.0 mmol·L-1KH2PO4)、中度缺磷(0.5 mmol·L-1KH2PO4,0.5 mmol·L-1KCl)和重度缺磷(0 mmol·L-1KH2PO4,1.0 mmol·L-1KCl),脅迫中不足的K+用等量KCl代替。其他營養(yǎng)元素含量根據(jù)吳鵬飛(2009)改良的1/3Hoagland營養(yǎng)液配方進(jìn)行補(bǔ)充,即5.0 mmol·L-1KNO3、 2.0 mmol·L-1MgSO4·7H2O、5.0 mmol·L-1Ca(NO3)2·4H2O、l mL·L-1Fe-EDTA, Arnon微量元素(46.3 μmol·L-1H3BO3、0.3 μmol·L-1CuSO4·5H2O、0.8 μmol·L-1ZnSO4·7H2O、9.1 μmol·L-1MnC12·4H2O、0.4 μmol·L-1H2MoO4·4H2O),調(diào)節(jié)營養(yǎng)液pH為5.8。

    1.2ClPht1;1基因全長的克隆

    1.2.1ClPht1;1基因保守區(qū)域擴(kuò)增 采用RNA plant plus植物總RNA提取試劑盒(天根公司)提取重度缺磷處理24 h的32號(hào)磷高效杉木根系總RNA,并使用Fermentas公司RT-PCR試劑盒合成cDNA第1條鏈。

    根據(jù)油茶(GenBank: EU496869.1)、擬南芥(GenBank: AB005746)、玉米(GenBank: NM001111799.1)等植物PHT1基因氨基酸序列中的保守區(qū),利用Primer Premier 5.0軟件設(shè)計(jì)簡并引物ClPht1-R1∶5′-GTGCTGAATGCACTCGATGT-3′和ClPht1-F1∶5′-CCTAGCTGGGAAAATCTCAGC-3′對(duì)杉木根系cDNA進(jìn)行PCR擴(kuò)增,反應(yīng)程序?yàn)椋?95 ℃預(yù)變性3 min; 94 ℃ 變性40 s,54 ℃退火40 s,72 ℃ 延伸100 s,35個(gè)循環(huán); 72 ℃延伸 7 min。擴(kuò)增產(chǎn)物經(jīng)天根公司DNA凝膠回收試劑盒回收、純化后,連接到pGEM-T Easy(TaKaRa公司)載體上進(jìn)行測序,獲得目的基因的核心片段。

    1.2.2ClPht1;1基因全長cDNA的獲得 采用Invitrogen公司5′RACE(Version2.0)試劑盒,設(shè)計(jì)特異引物(5′-1:5′-TCGGAGGCAGAGAACC-3′; 5′-2:5′-TCCCAGGAGCCTCGTCAC-3′; 5′-3:5′-GCAGAAG AGGTCATACGCATC-3′)克隆目的基因的5′端序列。使用Superscript Ⅱ RT 酶和引物5′-1對(duì)總RNA進(jìn)行目的基因第1鏈cDNA的合成,使用引物5′-2和橋連鉚釘引物AAP(5′-GGCCACGCGTCGACTAGT ACGGGIIGGGIIGGGIIG-3′)對(duì)已經(jīng)加dC尾的cDNA進(jìn)行PCR第1輪擴(kuò)增,使用引物5′-3和橋連通用擴(kuò)增引物AUAP(5′-GGCCACGCGTC GACTAGTAC-3′)進(jìn)行巢式PCR第2輪擴(kuò)增,擴(kuò)增產(chǎn)物連接載體及測序過程同上。

    采用Clontech 公司的SMARTerTMRACE試劑盒,設(shè)計(jì)特異引物(3′-1:5′-ACCTGGTTTCTTCT GGACATCGCCTACT-3′; 3′-2: 5′-AACCTGTTTCAA AAGGACATCTTCACGG-3′)克隆目的基因的3′端序列。使用逆轉(zhuǎn)錄酶SmartScribeTMReverse Transcriptase 和引物3′CDS primer A(5′-AAGCAGTGGTATC AACGCAGACTAC(T)30 V N-3′)對(duì)總RNA進(jìn)行逆轉(zhuǎn)錄合成cDNA,并以此cDNA為模板,使用引物3′-1和UPM(5′-CTAATACGACTCACTATAGGGCA AGCAGTGGTATCAACGCAGAGT-3′)進(jìn)行第1輪PCR擴(kuò)增。擴(kuò)增產(chǎn)物稀釋50倍后,用引物3′-2和 UPM進(jìn)行第2輪PCR擴(kuò)增,擴(kuò)增產(chǎn)物連接載體及測序過程同上。

    根據(jù)獲得的目的基因核心片段、3′末端和5′末端3 個(gè)序列信息,利用Vector NTI10.3.0軟件進(jìn)行序列拼接,獲得目的全長基因的cDNA 序列。設(shè)計(jì)上游引物QC288F:5′-TTGAACCATGGCGGAGAA GC-3′和下游引物QC288R:5′-CTAAGCTGGT GCGGTTCTGG-3′,采用Takara Taq polymerase對(duì)目的基因全長進(jìn)行PCR擴(kuò)增,反應(yīng)程序?yàn)?4 ℃預(yù)變性5 min; 94 ℃變性60 s,50 ℃退火60 s,72 ℃ 延伸180 s,30個(gè)循環(huán); 72 ℃延伸 10 min。擴(kuò)增產(chǎn)物連接載體及測序過程同上。

    1.2.3ClPht1;1基因序列分析 通過DNAMAN軟件預(yù)測測序所得目的基因編碼的氨基酸序列,并使用ORF Finder預(yù)測杉木PHT1基因的最長開放閱讀框。使用NCBI在線Blast檢索其他植物PHT1基因編碼氨基酸序列,并采用Clustal X2.1軟件進(jìn)行多序列比對(duì)。使用ProtParam軟件和ProtScale軟件預(yù)測目的基因編碼蛋白質(zhì)的分子量、等電點(diǎn)、不穩(wěn)定系數(shù)、疏水性等理化性質(zhì); 使用TMpredServer軟件和MobylePortal軟件預(yù)測目的基因編碼蛋白質(zhì)的跨膜結(jié)構(gòu)和拓?fù)浣Y(jié)構(gòu); 使用SOPMA軟件和Swiss-model軟件預(yù)測目的基因編碼蛋白質(zhì)的二級(jí)和三級(jí)結(jié)構(gòu)。

    1.3ClPht1;1基因表達(dá)分析

    采用實(shí)時(shí)熒光定量PCR檢測目的基因在32號(hào)磷高效杉木家系根、莖、葉中的表達(dá)情況; 檢測中度缺磷脅迫下目的基因在不同磷利用效率杉木4號(hào)(1.14±0.08)kg·g-1、15號(hào)(1.10±0.06)kg·g-1、25號(hào)(1.77±0.22)kg·g-1、27號(hào)(1.38±0.03)kg·g-1、28號(hào)(0.73±0.06)kg·g-1、32號(hào)(1.09±0.02)kg·g-1家系根系中的表達(dá)差異; 檢測中度、重度供磷脅迫下0,12,24,48 h,以及恢復(fù)正常供磷12,24,48 h目的基因在32號(hào)磷高效杉木家系根系中的表達(dá)量變化,以上試驗(yàn)均設(shè)置3次生物學(xué)重復(fù)。實(shí)時(shí)熒光定量PCR在Roch480熒光定量PCR儀上進(jìn)行,熒光染料采用TakaRa公司SYBR Green PCR Master Mix,上游引物5′-GGCCATTCCGTACAACCACTG-3′,下游引物5′-AATTCTTGACCCCAATACCCT-3′。 使用β-Actin基因作為內(nèi)參基因,其上游引物為5′-TACATTGCCGGTGTATTGAACGTC-3′,下游引物為5′-AGTTGCTCCAGAAGAACATCC-3′。PCR反應(yīng)程序?yàn)?5 ℃預(yù)變性3 min; 95 ℃變性20 s,60 ℃退火30 s,72 ℃延伸30 s,45個(gè)循環(huán); 95 ℃ 30 s,60 ℃ 30 s,逐漸升溫到95 ℃,速度0.11 ℃·s-1。擴(kuò)增結(jié)果采用2-△△Ct法計(jì)算目的基因的相對(duì)表達(dá)水平。

    2 結(jié)果與分析

    2.1ClPht1;1基因全長克隆

    以重度缺磷處理24 h的32號(hào)磷高效杉木根系cDNA為模板,采用設(shè)計(jì)的簡并引物ClPht1-R1和ClPht1-F1進(jìn)行PCR擴(kuò)增測序,獲得1個(gè)長度為1 004 bp的杉木PHT1基因保守片段序列(圖1A)。根據(jù)獲得的杉木PHT1基因保守序列分別設(shè)計(jì)目的基因的5′和3′ RACE特異引物,5′RACE-PCR擴(kuò)增得到長度為250 bp的片段(圖1B),3′RACE-PCR擴(kuò)增得到長度為1 000 bp片段(圖1C)。將目的基因保守序列、5′RACE片段和3′RACE片段拼接后獲得全長為1 967 bp的基因序列。依據(jù)拼接獲得的基因序列設(shè)計(jì)全長引物H1和H2,擴(kuò)增目的基因全長并測序驗(yàn)證,在預(yù)測位置獲得單一清晰的目的條帶(圖1D),測序結(jié)果與拼接的序列一致,命名為杉木高親和磷轉(zhuǎn)運(yùn)蛋白基因ClPht1;1(GenBank登錄號(hào): KX302006)。

    圖1 杉木PHT1基因保守片段(A)、cDNA 5′RACE(B)、cDNA 3′RACE(C)和全長基因(D)Fig.1 Cunninghamia lanceolata PHT1 gene conservative fragment (A), cDNA 5′RACE (B), cDNA 3′RACE (C) and full length gene (D)

    2.2ClPht1;1基因生物信息學(xué)分析

    2.2.1ClPht1;1基因序列分析ClPht1;1基因(GenBank登錄號(hào): KX302006)全長序列分析表明,其長度為1 967 bp,開放閱讀框?yàn)? 638 bp,編碼1個(gè)含545個(gè)氨基酸的蛋白質(zhì),編碼蛋白分子量為59.95 kDa,等電點(diǎn)為9.04。ClPht1;1基因編碼多肽的原子組成為C2750H4198N698O751S28,不穩(wěn)定系數(shù)是33.44,表明其為穩(wěn)定蛋白質(zhì)分子。ClPht1;1基因的脂肪系數(shù)為82.57,總平均親水性為0.226,ProtScale 軟件分析表明ClPht1;1基因編碼蛋白含有較多的疏水區(qū)域,為疏水性蛋白(圖2)。

    圖3 ClPht1;1蛋白跨膜拓?fù)浣Y(jié)構(gòu)預(yù)測Fig.3 Prediction of ClPht1;1 protein transmembrane topology

    圖2 ClPht1;1疏水性曲線Fig.2 ClPht;1 hydrophobic curve

    2.2.2ClPht1;1基因編碼蛋白的結(jié)構(gòu)預(yù)測 利用在線TMpred Server 和Mobyle portal軟件對(duì)ClPht1;1基因編碼蛋白進(jìn)行跨膜結(jié)構(gòu)和跨膜區(qū)拓?fù)浣Y(jié)構(gòu)預(yù)測,結(jié)果表明ClPht1;1基因編碼蛋白具有12個(gè)跨膜域(TM),其中含1個(gè)疑似跨膜域,包括6 個(gè)N端的跨膜區(qū)以及6個(gè)C端的跨膜區(qū),中部(TM6和TM7之間)由1個(gè)中央親水環(huán)(hydrophilic loop)分隔(圖3)。每個(gè)跨膜結(jié)構(gòu)域由17~25個(gè)氨基酸殘基組成的螺旋構(gòu)成,跨膜蛋白的N 端、C端和中間的大環(huán)都分布在細(xì)胞質(zhì)中,保守序列位于第4個(gè)跨膜域,ClPht1;1基因編碼蛋白的跨膜區(qū)結(jié)構(gòu)與其他物種PHT1基因編碼蛋白的跨膜區(qū)結(jié)構(gòu)基本一致(Nussaumeetal., 2011)。使用SOPMA軟件和Swiss-model軟件對(duì)ClPht1;1基因編碼蛋白質(zhì)結(jié)構(gòu)進(jìn)行分析,結(jié)果表明其蛋白質(zhì)二級(jí)結(jié)構(gòu)中41.65%為α-螺旋,26.42%為與蛋白質(zhì)功能相關(guān)的無規(guī)則卷曲,21.47%為延伸鏈,10.46%為β-轉(zhuǎn)角; 三級(jí)結(jié)構(gòu)中構(gòu)成蛋白質(zhì)的主要骨架是α-螺旋(圖4); 保守域分析表明ClPht1;1基因編碼蛋白有典型的MFS結(jié)構(gòu)域。

    圖4 ClPht1;1蛋白三級(jí)結(jié)構(gòu)預(yù)測Fig.4 Prediction of ClPht1;1 protein tertiary structure

    2.2.3 ClPht1;1蛋白序列同源性分析 使用ClustalX2.1軟件對(duì)ClPht1;1基因和Blast獲得的其他植物PHT基因編碼氨基酸序列進(jìn)行多重比對(duì)(圖5)。結(jié)果表明,ClPht1;1基因編碼產(chǎn)物與日本柳杉(Cryptomeriajaponica)PHT基因編碼氨基酸序列的相似性為87%,與其他植物PHT基因編碼氨基酸序列的相似性均在70%以上。其中,與馬尾松(Pinusmassoniana)和可可樹(Theobromacacao)等南方喬木PHT基因編碼氨基酸序列相似性達(dá)到77%,與油茶、葡萄(Vitisvinifera)、胡楊(Populuseuphratica)等相似性達(dá)到76%。

    2.3ClPht1;1基因的表達(dá)分析

    采用實(shí)時(shí)熒光定量PCR檢測正常供磷情況下ClPht1;1基因在32號(hào)磷高效杉木家系根、莖、葉中的表達(dá)情況,結(jié)果表明ClPht1;1基因在杉木各組織部位均有表達(dá),但在根中的表達(dá)量要遠(yuǎn)高于莖和葉,在葉片中的表達(dá)量較低(圖6)。中度缺磷脅迫(0.5 mmol·L-1KH2PO4)下,ClPht1;1基因在不同磷利用效率杉木4號(hào)(1.14±0.08)kg·g-1、15號(hào)(1.10±0.06)kg·g-1、25號(hào)(1.77±0.22)kg·g-1、27號(hào)(1.38±0.03)kg·g-1、28號(hào)(0.73±0.06)kg·g-1、32號(hào)(1.09±0.02)kg·g-1家系根系中的表達(dá)量為25號(hào)>27號(hào)>4號(hào)>15號(hào)>32號(hào)>28號(hào)(圖7),表明ClPht1;1基因在杉木根部的表達(dá)量隨著不同家系杉木磷利用效率的增強(qiáng)而升高。中度、重度(0 mmol·L-1KH2PO4)缺磷脅迫下ClPht1;1基因在32號(hào)杉木家系根系中的表達(dá)量隨誘導(dǎo)時(shí)間(0,12,24,48 h)的延長均呈現(xiàn)明顯的上調(diào)趨勢; 恢復(fù)供磷12 h后ClPht1;1基因表達(dá)量明顯下降,并在恢復(fù)供磷24 h和48 h后呈現(xiàn)與正常供磷培養(yǎng)相近的表達(dá)量水平。重度缺磷脅迫下,ClPht1;1基因表達(dá)量要高于其在中度缺磷脅迫下的表達(dá)量,表明ClPht1;1基因的表達(dá)明顯受介質(zhì)低磷脅迫誘導(dǎo)(圖8)。

    圖5 ClPht1;1及部分植物PHT基因編碼蛋白的氨基酸序列多重比對(duì)Fig.5 Multiple alignments of the amino acid sequences of ClPht1;1 and some plant PHT genes encoding proteins羅馬數(shù)字Ⅰ-Ⅻ表示12 個(gè)跨膜區(qū),虛框內(nèi)表示疑似跨膜區(qū); ■表示N-糖基化位點(diǎn); ▲表示蛋白激酶C磷酸化位點(diǎn); ★表示酪蛋白激酶Ⅱ磷酸化位點(diǎn)。Roman numerals Ⅰ-Ⅻ represent 12 transmembrane regions, virtual boxes represent the suspected transmembrane areas; ■ represents N-glycosylation site; ▲ represents protein kinase C phosphorylation site; ★ represents casein kinase Ⅱ phosphorylation site.

    圖6 ClPht1;1在32號(hào)杉木不同組織中的表達(dá)分析Fig.6 Expression analysis of ClPht1;1 in different tissues of Cunninghamia lanceolata No. 32

    圖7 ClPht1;1在不同家系杉木根部的表達(dá)分析Fig.7 Expression analysis of ClPht1;1 in roots in different families of Cunninghamia lanceolata

    圖8 不同磷脅迫和處理時(shí)間下ClPht1;1在32號(hào)杉木根系的表達(dá)分析Fig.8 Expression analysis of ClPht1;1 in roots of Cunninghamia lanceolata No. 32 under different P stresses at different treatment time

    3 討論

    PHT1 基因家族成員大多主要在植物根部表達(dá),一些PHT1 基因不僅在植物根部表達(dá),還在植物的莖、葉、花、塊莖和子葉中表達(dá),這些PHT1基因可能涉及到植株體內(nèi)的磷轉(zhuǎn)運(yùn)功能(Liuetal., 2005)。在擬南芥中,PHT1;1和PHT1;9均受缺磷脅迫的誘導(dǎo)表達(dá),除了PHT1;6在花粉中表達(dá)外,其他成員均主要在根中表達(dá)(Bayleetal., 2011)。馬鈴薯(Solanumtuberosum)StPT1和番茄(Lycopersiconesculentum)LePT1(Nagyetal., 2005)在缺磷或不缺磷的根部和地上部均能大量表達(dá),但缺磷能誘導(dǎo)增強(qiáng)其表達(dá)量;StPT2和LePT2則僅在植株的根部受缺磷環(huán)境的誘導(dǎo)表達(dá)(Leggewieetal., 1997; Liuetal., 1998; Gordon-Weeksetal., 2003)。油茶BhPHT1;4基因在嫩根、老根、莖、新葉、老葉中都有表達(dá),但表達(dá)量依次遞減(趙彩芝, 2014)。本試驗(yàn)結(jié)果表明ClPht1;1基因與大部分植物PHT1基因類似,主要在杉木根部表達(dá),葉片中僅能檢測到少量表達(dá)信號(hào)。ClPht1;1基因在杉木根、莖和葉片中均有表達(dá),說明其不僅參與杉木對(duì)土壤磷的吸收,而且可能參與杉木體內(nèi)磷素的轉(zhuǎn)運(yùn)。植物對(duì)磷的轉(zhuǎn)運(yùn)效率受細(xì)胞膜上磷轉(zhuǎn)運(yùn)蛋白的數(shù)量及其對(duì)磷的親和性決定,在中度缺磷脅迫下,ClPht1;1基因在不同磷利用效率4號(hào)、15號(hào)、25號(hào)、27號(hào)、28號(hào)、32號(hào)家系杉木根系中的表達(dá)量存在差異,表達(dá)量隨著不同家系杉木磷利用效率的增強(qiáng)而升高,這可能與不同磷利用效率杉木根部細(xì)胞膜表面ClPht1;1磷轉(zhuǎn)運(yùn)蛋白數(shù)量多少有關(guān)。ClPht1;1基因的表達(dá)量隨根系缺磷脅迫時(shí)間的延長和脅迫程度的加重而升高,在恢復(fù)供磷培養(yǎng)后表達(dá)量下降逐漸至正常水平,表明ClPht1;1基因的表達(dá)受低磷脅迫的誘導(dǎo)。綜上研究,ClPht1;1基因的時(shí)空表達(dá)情況符合PHT1家族基因的表達(dá)特性,對(duì)杉木PHT1基因的表達(dá)量分析為磷高效轉(zhuǎn)運(yùn)基因型的篩選提供了參考。

    4 結(jié)論

    磷是我國杉木人工林產(chǎn)量和質(zhì)量的主要限制因子,由于南方紅壤區(qū)土壤有效磷長期不足,杉木長期處于低磷的生長環(huán)境中。本試驗(yàn)從32號(hào)磷高效利用杉木家系的根部克隆獲得一個(gè)高親和磷轉(zhuǎn)運(yùn)蛋白基因ClPht1;1,其具有植物PHT1基因家族的典型結(jié)構(gòu)特點(diǎn)和表達(dá)模式。對(duì)于長期處于低磷環(huán)境中的杉木,ClPht1;1在杉木高親和轉(zhuǎn)運(yùn)土壤低濃度有效磷中起到重要的作用。目前,杉木PHT1基因家族的成員組成和各自的磷轉(zhuǎn)運(yùn)功能仍有待研究,深入解析杉木高親和磷轉(zhuǎn)運(yùn)蛋白基因的結(jié)構(gòu)和功能,對(duì)于了解杉木磷吸收的代謝機(jī)制,培育磷高效利用杉木品種和增強(qiáng)林分磷素管理都具有重要意義。

    賈興永, 李菊梅. 2011. 土壤磷有效性及其與土壤性質(zhì)關(guān)系的研究. 中國土壤與肥料, (6): 76-82.

    (Jia X Y, Li J M. 2011. Study on the relationship between soil phosphorus availability and soil properties. Soil and Fertilizer in China, (6): 76-82. [in Chinese ])

    李 慧,叢 郁,常有宏,等. 2013. 豆梨磷轉(zhuǎn)運(yùn)蛋白質(zhì)基因(PcPht1)的克隆、表達(dá)及啟動(dòng)子分析. 江蘇農(nóng)業(yè)學(xué)報(bào),29 (4): 842-850.

    (Li H, Cong Y, Chang Y H,etal. 2013. Cloning, expression and promoter analysis of phosphorus transport protein gene (PcPht1) of soybean and pear. Jiangsu Journal of Agricultural Sciences,29 (4): 842-850. [in Chinese ])

    李立芹. 2011. 農(nóng)作物Pht1 家族磷轉(zhuǎn)運(yùn)體蛋白的生物信息學(xué)分析. 作物雜志,(3): 320-324.

    (Li L Q. 2011. Bioinformatics analysisi ofPht1 phosphate transporter protein family in crops. Crops, (3): 320-324. [in Chinese ])

    梁 霞,劉愛琴,馬祥慶,等. 2006. 不同杉木無性系磷素特性的比較. 植物生態(tài)學(xué)報(bào),30(6): 1005-1011.

    (Liang X, Liu A Q, Ma X Q,etal. 2006. Comparison of phosphorus characteristics of different clones of Chinese fir. Journal of Plant Ecology, 30(6): 1005-1011. [in Chinese ])

    湯銀輝,何 鵬. 2010. 巴西橡膠樹磷轉(zhuǎn)運(yùn)蛋白基因的克隆及生物信息學(xué)分析. 熱帶作物學(xué)報(bào),31(5): 758-766.

    (Tang Y H,He P. 2010. Molecular cloning and bioinformatic analysis of a gene encoding a phosphate transporter fromHeveabrasiliensis. Chinese Journal of Tropic Crops,31(5): 758-766. [in Chinese ])

    王 策,秦靜靜,甘紅豪,等. 2012. 毛果楊全基因組磷酸根轉(zhuǎn)運(yùn)蛋白家族成員序列分析. 浙江農(nóng)林大學(xué)學(xué)報(bào),29 (4): 516-526.

    (Wang C,Qin J J,Gan H H,etal. 2012. Genome wide analysis of phosphate transporter gene family inPopolustrichocarpa. Journal of Zhejiang Ariculure and Forestry University,29 (4): 516-526. [in Chinese ])

    王 萍, 陳愛群, 余 玲, 等. 2006. 植物磷轉(zhuǎn)運(yùn)蛋白基因及其表達(dá)調(diào)控的研究進(jìn)展. 植物營養(yǎng)與肥料學(xué)報(bào), 12 (4): 584-591.

    (Wang P,Chen A Q,Yu L,etal. 2006. Advance of plant phosphate transporter genes and their regulated expression. Plant Nutrition and Fertilizer Science,12 (4): 584-591. [in Chinese ])

    汪 攀. 2015. 低磷脅迫下杉木根系形成通氣組織與磷利用效率的關(guān)系研究. 福州: 福建農(nóng)林大學(xué)碩士學(xué)位論文.

    (Wang P. 2015. Study on the relationship between the formation of the root system of Chinese fir and the utilization efficiency of phosphorus under low phosphorus stress. Fuzhou: MS thesis of Fujian Agriculture and Forestry University. [in Chinese ])

    吳鵬飛,馬祥慶,陳友力,等. 2012. 杉木無性系測定林磷素利用效率的比較. 福建農(nóng)林大學(xué)學(xué)報(bào): 自然科學(xué)版, 41 (1): 40-45.

    (Wu P F, Ma X Q, Chen Y L,etal. 2012. Comparison of phosphorus use efficiency among clonal test plantations of Chinese fir. Journal of Fujian Agriculture and Forestry University: Natural Science Edition, 41 (1): 40-45. [in Chinese ])

    吳鵬飛. 2009. P高效利用杉木無性系適應(yīng)環(huán)境磷脅迫的機(jī)制研究. 福州: 福建農(nóng)林大學(xué)博士學(xué)位論文.

    (Wu P F. 2009. Adaptation mechanism of Chinese fir clones with high phosphorus-use-efficiency to environmental phosphorus stress. Fuzhou: PhD of Fujian Agriculture and Forestry University. [in Chinese ])

    楊存義, 劉 靈, 沈 宏, 等. 2006. 植物Pht1 家族磷轉(zhuǎn)運(yùn)子的分子生物學(xué)研究進(jìn)展. 分子植物育種, 4(2): 153-159.

    (Yang C Y, Liu L, Shen H,etal. 2006. Research progress on molecular biology of phosphorus transport in plantPht1 family. Molecular Plant Breeding, 4(2): 153-159.[in Chinese ])

    俞新妥. 2006. 中國杉木研究進(jìn)展 (2000—2005) I. 杉木生理生態(tài)研究綜述. 福建林學(xué)院學(xué)報(bào), 26 (2): 177-185.

    (Yu X T. Progress in research on Chinese fir (2000—2005) I. Chinese fir physiological and ecological research review. Journal of Fujian Forestry College, 26 (2): 177-185. [in Chinese ])

    趙彩芝. 2014. 油菜高親和磷轉(zhuǎn)運(yùn)蛋白基因BhPHT1;4的功能研究. 武漢: 華中師范大學(xué)碩士學(xué)位論文.

    (Zhao C Z. 2014. High affinity phosphate transporter geneBhPHT1;4 function of oilseed rape. Wuhan: MS thesis of Huazhong Normal University. [in Chinese ])

    周俊琴,譚曉風(fēng),袁 軍,等. 2013. 油茶Pht1;2的克隆及生物信息學(xué)分析. 生物技術(shù), (1): 37-42.

    (Zhou J Q, Tan X F, Yuan J,etal. 2013. Isolation and bioinformatics analysis ofCoPht1;2 from oiltea. Biotechnology, (1): 37-42. [in Chinese ])

    Ai P,Sun S,Zhao J,etal. 2009. Two rice phosphate transporters,OsPht1;2 and OsPht1;6,have different functions and kinetic properties in uptake and translocation. Plant Journal for Cell & Molecular Biology, 57 (5): 798-809.

    Bayle V, Arrighi J F, Creff A,etal. 2011.Arabidopsisthalianahigh affinity phosphate transporters exhibit multiple levels of posttranslational regulation. Plant Cell, 23(4): 1523-1535.

    Calderon-Vazquez C, Ibarra-Laclette E, Caballero-Perez J,etal. 2008. Transcript profiling ofZeamaysroots reveals gene responses to phosphate deficiency at the plant- and species-specific levels. Journal of Experimental Botany, 59(9): 2479-2497.

    Gordon-Weeks R, Tong Y, Davies T G,etal. 2003. Restricted spatial expression of a high-affinity phosphate transporter in potato roots. Journal of Cell Science, 116(15): 3135-3144.

    Leggewie G, Willmitzer L, Riesmeier J W. 1997. Two cDNAs from potato are able to complement a phosphate uptake-deficient yeast mutant: Identification of phosphate transporters from higher plants. Plant Cell, 9(3): 381- 392.

    Liu C M, Muchhal U S, Uthappa M,etal. 1998. Tomato phosphate transporter genes are differentially regulated in plant tissues by phosphorus. Plant Physiology, 116(1): 91-99.

    Liu J Q,Samac D A,Bucciarelli B,etal. 2005. Signaling of phosphorus deficiency induced gene expression in white lupin requires sugar and phloem transport. Plant Journal, 41 (2): 257 -268.

    Mitsukawa N,Okumura S,Shirano Y,etal. 1997. High-affinity phosphate transporter genes ofArabidopsisthaliana. Proceedings of the National Academy of Sciences, 94(13): 7098-7102.

    Nagy R, Karandashov V, Chague V,etal. 2005. The characterization of novel mycorrhiza-specific phosphate transporters fromLycopersiconesculentumandSolanumtuberosumuncovers functional redundancy in symbiotic phosphate transport in solanaceous species. Plant Journal, 42 (2): 236-250.

    Nussaume L, Kanno S, Javot H,etal. 2011. Phosphate import in plants: focus on thePHT1 transporters. Frontiers in Plant Science, 2(12): 1085-1091.

    Poirier Y, Jung J. 2015. Phosphate transporters// Plaxton W C, Lambers H. Annual plant reviews Volume 48: Phosphorus metabolism in plants. Wiley,125-158.

    Preuss C P,Huang C Y,Gilliham M,etal.2010. Channel-like characteristics of the low-affinity barley phosphate transporterPHT1;6 when expressed inXenopusoocytes.Plant Physiology,152(3): 1431-1441.

    Scheible W R, Rojas-Triana M. 2015. Sensing, signaling, and control of phosphate starvation in plants: molecular players and applications// Plaxton W C, Lambers H. Annual plant reviews Volume 48: Phosphorus metabolism in plants. Wiley, 25-64.Wu Z Y,Zhao J M,Gao R F,etal. 2011. Molecular cloning, characterization and expression analysis of two members of thePht1 family of phosphate transporters inGlycinemax. PLoS One, 6 (6): e19752.Yang X J,F(xiàn)innegan P M. 2010. Regulation of phosphate starvation responses in higher plants. Annals of Botany, 105(4): 513-526.

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

    Cloning and Expression Analysis of Phosphate Transporter GeneClPht1;1 inCunninghamialanceolata

    Su Shuoshuo1Li Ming1,2Wu Pengfei1,2Zhang Ying1Ma Xiangqing1,2

    (1.College of Forestry, Fujian Agriculture and Forestry University Fuzhou 350002; 2.State Forestry Administration Engineering Research Center of Chinese Fir Fuzhou 350002)

    【Objective】PHT1 gene family is one of the important phosphorus transporters which affect the uptake of phosphorus in plants. ACunninghamialanceolataphosphate transporter gene ofPHT1 gene family was cloned from No.32 phosphorus efficientC.lanceolatacDNA,and the temporal and spatial expression ofClPht1;1 under different levels of phosphorus stress was studied to lay the foundation for the research ofPHT1 gene sequence characteristics, functional structure and the selection of theC.lanceolatagenotype which is efficient in phosphorus utilization.【Method】According to the designed degenerate primer based on the conserved sequence ofPHT1 gene family, the target geneClPht1;1 cDNA sequence was obtained by amplification of the root cDNA of the P-efficient family No.32 ofC.lanceolataas the template, the full length of the target gene was cloned by RACE technique, and analyzed its sequence characteristics, homology and protein structure were then analyzed. The real-time fluorescence quantitative PCR was used to detect the expression ofClPht1;1 in roots, stems and leaves of No.32C.lanceolatawith high P-efficiency, differential expressions ofClPht1;1 in the roots ofC.lanceolatafamilies of No.4,No.15,No.25,No.27,No.28 and No.32 in different P-utilization efficiency under moderate P deficiency stress, and theClPht1;1 expression change with time in the root system of phosphorus efficientC.lanceolataNo.32 under moderate and severe phosphorus deficiency stress.【Result】AC.lanceolataPHT1 gene was cloned and namedClPht1;1 (GenBank accession number KJ493165). The coding region of the gene sequence was 1 638 bp, code 545 aa protein.ClPht1;1 encoded proteins were composed of 12 hydrophobic transmembrane domains, including 1 suspected transmembrane domains. Each transmembrane domain was composed of 17-25 amino acid residues, and the N and C ends of transmembrane proteins were located in the cytoplasm, and the conserved sequences were located in the fourth transmembrane domains. The main skeleton of the protein was α-helix, no signal peptide sequence. The similarity of amino acid sequences of encoding proteins ofClPht1;1 was 87% withCryptomeriajaponicaPHTgene, and above 75% withPopulus,Camellia,Pinusmassonianaetc.PHTfamily gene.ClPht1;1 genes were expressed in root, stem and leaf ofC.lanceolata, and the highest expression was in root,the lowest was in leaf; Under moderate P deficiency stress, the expression level ofClPht1;1 in the root ofC.lanceolatavaried among families in a order of No.25>No.27>No.4>No.15>No.32>No.28. Under moderate and severe P deficiency stresses, the expression ofClPht1;1 in the root of No.32 increased with the extension of stress time; The expression level ofClPht1;1 was gradually decreased to normal level after the recovery of phosphorus supply;ClPht1;1 gene expression in severe P deficiency stress was higher than that under moderate P deficiency stress.【Conclusion】ClPht1;1 gene has the typical structure ofPHT1 gene family, and its protein amino acid sequence is highly similar to the amino acid sequences ofCryptomeriajaponicaphosphate transporter protein, it is a member of the high affinity phosphate transporterPHT1 gene family ofC.lanceolata.ClPht1;1 gene was mainly expressed in root ofC.lanceolata, and its expression was lower in leaf. The stronger the phosphorus use efficiency, the higher theClPht1;1 gene expression in its roots.There were large differences in the expression level ofClPht1;1 gene among families with different phosphorus use efficiencies.ClPht1;1 gene expression was induced by low phosphorus stress; Under P deficiency stress, the expression level ofClPht1;1 was significantly increased, and the expression ofClPht1;1 gene was significantly decreased after the recovery of phosphorus supply.

    Cunninghamialanceolata; phosphate transporter;PHT1;1; RACE; RT-qPCR; phosphorus use efficiency

    10.11707/j.1001-7488.20170505

    2016-06-15;

    2016-08-12。

    國家自然科學(xué)基金海峽聯(lián)合基金項(xiàng)目(U1405211); 國家林業(yè)局杉木工程技術(shù)研究中心科技成果孵化基金(6213C0111)。

    S718.46

    A

    1001-7488(2017)05-0033-10

    *李明為通訊作者。

    猜你喜歡
    家系杉木利用效率
    杉木黃化病的防治技術(shù)措施研究
    避免肥料流失 提高利用效率
    杉木萌芽更新關(guān)鍵技術(shù)
    杉木育苗化學(xué)防除雜草技術(shù)
    體制改革前后塔里木河流域水資源利用效率對(duì)比分析
    杉木半同胞24年生優(yōu)良家系選擇
    馬氏珠母貝紅色閉殼肌F1代的家系選育及家系評(píng)定
    一個(gè)非綜合征型聾家系的分子病因?qū)W研究
    肝豆?fàn)詈俗冃?個(gè)家系的基因突變分析
    兩個(gè)雄激素不敏感綜合征家系中AR基因突變檢測
    一区在线观看完整版| 美国免费a级毛片| 国产乱人伦免费视频| 人人妻,人人澡人人爽秒播| 国产私拍福利视频在线观看| 亚洲五月天丁香| 欧美中文日本在线观看视频| 亚洲 欧美一区二区三区| 亚洲av五月六月丁香网| 男女下面进入的视频免费午夜 | 一a级毛片在线观看| 露出奶头的视频| 精品欧美国产一区二区三| 国产成人啪精品午夜网站| 好看av亚洲va欧美ⅴa在| 亚洲 国产 在线| 欧美黑人欧美精品刺激| 99riav亚洲国产免费| 中出人妻视频一区二区| 老司机深夜福利视频在线观看| 亚洲成人久久性| 日本 av在线| 搡老妇女老女人老熟妇| 成人亚洲精品av一区二区| 丝袜美足系列| 亚洲少妇的诱惑av| 日韩欧美免费精品| 欧美黄色淫秽网站| 亚洲国产精品合色在线| 久久国产乱子伦精品免费另类| 亚洲av日韩精品久久久久久密| 欧美成狂野欧美在线观看| 黄色丝袜av网址大全| 色老头精品视频在线观看| 日韩欧美国产一区二区入口| 亚洲九九香蕉| 97超级碰碰碰精品色视频在线观看| 窝窝影院91人妻| 妹子高潮喷水视频| 国产黄a三级三级三级人| 又大又爽又粗| 久久久久精品国产欧美久久久| 成人国产综合亚洲| 午夜福利一区二区在线看| 狠狠狠狠99中文字幕| 99久久精品国产亚洲精品| 国产精品野战在线观看| 国产精品综合久久久久久久免费 | 日本撒尿小便嘘嘘汇集6| 啦啦啦韩国在线观看视频| 国产亚洲av高清不卡| avwww免费| 国产精品自产拍在线观看55亚洲| 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 欧美+亚洲+日韩+国产| 88av欧美| 黄片小视频在线播放| 大型黄色视频在线免费观看| 国产激情欧美一区二区| 丝袜美足系列| 国产精品98久久久久久宅男小说| 亚洲人成网站在线播放欧美日韩| 纯流量卡能插随身wifi吗| 99国产极品粉嫩在线观看| 一区福利在线观看| 男女床上黄色一级片免费看| 国产精品久久久久久人妻精品电影| 麻豆国产av国片精品| 自拍欧美九色日韩亚洲蝌蚪91| www.熟女人妻精品国产| 亚洲成国产人片在线观看| 国产麻豆成人av免费视频| 亚洲国产中文字幕在线视频| 啦啦啦韩国在线观看视频| 国语自产精品视频在线第100页| 亚洲狠狠婷婷综合久久图片| 国产一区二区三区综合在线观看| 国产三级黄色录像| 精品国产乱码久久久久久男人| 免费看十八禁软件| 俄罗斯特黄特色一大片| 久久精品成人免费网站| 国产成人av教育| 成人亚洲精品av一区二区| 久久人人爽av亚洲精品天堂| 亚洲国产日韩欧美精品在线观看 | 18禁美女被吸乳视频| 国产成人精品久久二区二区91| 国产熟女午夜一区二区三区| 国产精品永久免费网站| 波多野结衣巨乳人妻| 欧美不卡视频在线免费观看 | 操出白浆在线播放| 久久影院123| 女性被躁到高潮视频| 久久人人97超碰香蕉20202| av在线播放免费不卡| 99精品欧美一区二区三区四区| 欧美激情高清一区二区三区| 欧美乱妇无乱码| 成熟少妇高潮喷水视频| 91麻豆av在线| 国产精品野战在线观看| 日韩国内少妇激情av| 国产亚洲精品久久久久久毛片| 女人被躁到高潮嗷嗷叫费观| 亚洲午夜精品一区,二区,三区| 日韩中文字幕欧美一区二区| 日本欧美视频一区| 一边摸一边抽搐一进一小说| 久久 成人 亚洲| 婷婷精品国产亚洲av在线| 极品教师在线免费播放| 亚洲第一电影网av| 夜夜看夜夜爽夜夜摸| 天堂动漫精品| 十分钟在线观看高清视频www| 搡老熟女国产l中国老女人| 日本黄色视频三级网站网址| 亚洲 欧美 日韩 在线 免费| 波多野结衣高清无吗| 一级a爱片免费观看的视频| 国产欧美日韩一区二区三| 国产免费av片在线观看野外av| 久久中文字幕人妻熟女| 免费观看人在逋| 桃红色精品国产亚洲av| 制服人妻中文乱码| netflix在线观看网站| 国产精品98久久久久久宅男小说| 此物有八面人人有两片| 婷婷六月久久综合丁香| 91av网站免费观看| 好男人在线观看高清免费视频 | 麻豆国产av国片精品| 色综合欧美亚洲国产小说| 9热在线视频观看99| 国产欧美日韩一区二区三| 一区在线观看完整版| 色播在线永久视频| 国产又爽黄色视频| 欧美色视频一区免费| 亚洲中文日韩欧美视频| 久久狼人影院| 国产熟女午夜一区二区三区| 好看av亚洲va欧美ⅴa在| 一区福利在线观看| 精品欧美国产一区二区三| 黑人操中国人逼视频| 一级a爱片免费观看的视频| 欧美成人一区二区免费高清观看 | 黑丝袜美女国产一区| 成年女人毛片免费观看观看9| 国产精品1区2区在线观看.| 日本精品一区二区三区蜜桃| 可以免费在线观看a视频的电影网站| 亚洲av五月六月丁香网| 视频区欧美日本亚洲| 欧美中文日本在线观看视频| 黑人巨大精品欧美一区二区mp4| 中文字幕人妻丝袜一区二区| 精品日产1卡2卡| 精品不卡国产一区二区三区| 亚洲片人在线观看| 狠狠狠狠99中文字幕| 欧美色视频一区免费| 亚洲情色 制服丝袜| 国产亚洲精品av在线| 亚洲中文字幕日韩| 久久精品国产清高在天天线| 日韩av在线大香蕉| 国产单亲对白刺激| 国产激情欧美一区二区| 女人被躁到高潮嗷嗷叫费观| 国产精品二区激情视频| 欧美成人性av电影在线观看| 99久久综合精品五月天人人| 免费观看精品视频网站| 亚洲五月色婷婷综合| 亚洲成人精品中文字幕电影| 亚洲中文av在线| a级毛片在线看网站| 久久香蕉激情| 亚洲欧美精品综合久久99| 精品卡一卡二卡四卡免费| 女同久久另类99精品国产91| avwww免费| 亚洲人成77777在线视频| 亚洲精华国产精华精| 亚洲第一欧美日韩一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 成人国产一区最新在线观看| 亚洲欧美一区二区三区黑人| 久久青草综合色| 村上凉子中文字幕在线| 午夜福利在线观看吧| 丁香六月欧美| 男女做爰动态图高潮gif福利片 | 色婷婷久久久亚洲欧美| 欧美日韩中文字幕国产精品一区二区三区 | 欧美大码av| 人人妻,人人澡人人爽秒播| a级毛片在线看网站| 麻豆一二三区av精品| 精品国产一区二区三区四区第35| 麻豆国产av国片精品| 少妇熟女aⅴ在线视频| 日本精品一区二区三区蜜桃| 亚洲avbb在线观看| 精品人妻1区二区| 又大又爽又粗| 757午夜福利合集在线观看| 欧美黑人欧美精品刺激| 国产精品 欧美亚洲| 国产精品综合久久久久久久免费 | 在线观看舔阴道视频| 国产精品久久久久久精品电影 | 99国产精品99久久久久| 88av欧美| 亚洲国产中文字幕在线视频| 亚洲自拍偷在线| 两性夫妻黄色片| 不卡av一区二区三区| 中文字幕人妻丝袜一区二区| www.熟女人妻精品国产| 免费在线观看影片大全网站| 国产午夜精品久久久久久| 十八禁人妻一区二区| 婷婷丁香在线五月| 中文字幕高清在线视频| 亚洲电影在线观看av| 亚洲激情在线av| 最近最新免费中文字幕在线| 国内久久婷婷六月综合欲色啪| 美国免费a级毛片| 国产一区二区三区视频了| 精品乱码久久久久久99久播| 香蕉国产在线看| 久久人妻av系列| 精品一区二区三区视频在线观看免费| 免费av毛片视频| 日韩高清综合在线| 校园春色视频在线观看| 久久伊人香网站| 757午夜福利合集在线观看| 精品电影一区二区在线| 91精品三级在线观看| 精品欧美国产一区二区三| 一级毛片精品| 韩国精品一区二区三区| 国产精品爽爽va在线观看网站 | 亚洲性夜色夜夜综合| 久久久精品国产亚洲av高清涩受| 成人欧美大片| 大型黄色视频在线免费观看| 90打野战视频偷拍视频| 免费在线观看亚洲国产| 中文字幕精品免费在线观看视频| 欧美黑人精品巨大| 老司机午夜福利在线观看视频| 亚洲电影在线观看av| 日日干狠狠操夜夜爽| 波多野结衣巨乳人妻| 久热这里只有精品99| 一进一出抽搐动态| 色哟哟哟哟哟哟| 国产精品 国内视频| 国产亚洲精品一区二区www| 一二三四社区在线视频社区8| 成人三级做爰电影| 波多野结衣av一区二区av| 男人操女人黄网站| 午夜两性在线视频| 色精品久久人妻99蜜桃| 手机成人av网站| 亚洲av电影在线进入| 久久久久久久久免费视频了| 18禁国产床啪视频网站| 日韩欧美在线二视频| 99国产极品粉嫩在线观看| 国产av又大| 国产精品免费视频内射| 亚洲国产精品成人综合色| 亚洲中文日韩欧美视频| av有码第一页| 免费观看精品视频网站| 99国产精品免费福利视频| 欧美国产日韩亚洲一区| avwww免费| 久久久久久亚洲精品国产蜜桃av| 嫩草影视91久久| 亚洲人成电影观看| 在线观看免费视频日本深夜| 国产精品久久视频播放| 九色亚洲精品在线播放| 久久久久久久久久久久大奶| 久久久久国产精品人妻aⅴ院| 久久久久久大精品| 一区二区三区国产精品乱码| 国产三级黄色录像| 国产成人免费无遮挡视频| 深夜精品福利| 久久国产精品人妻蜜桃| 欧美午夜高清在线| 日本在线视频免费播放| 日本 av在线| 日本vs欧美在线观看视频| 女人被躁到高潮嗷嗷叫费观| 69精品国产乱码久久久| 国产一级毛片七仙女欲春2 | 久久国产精品人妻蜜桃| 欧美老熟妇乱子伦牲交| 日本三级黄在线观看| 久久国产精品男人的天堂亚洲| 免费在线观看亚洲国产| 日本免费一区二区三区高清不卡 | 午夜影院日韩av| 国产欧美日韩综合在线一区二区| 日本免费一区二区三区高清不卡 | 亚洲美女黄片视频| aaaaa片日本免费| ponron亚洲| 在线观看免费日韩欧美大片| 亚洲三区欧美一区| 午夜福利视频1000在线观看 | 女警被强在线播放| 亚洲 国产 在线| 69精品国产乱码久久久| 午夜视频精品福利| 久久香蕉激情| 亚洲国产欧美日韩在线播放| 国产熟女午夜一区二区三区| 女人被躁到高潮嗷嗷叫费观| 在线观看日韩欧美| 日韩免费av在线播放| 男人操女人黄网站| 久久久久久人人人人人| 久久久国产成人免费| 999久久久精品免费观看国产| 午夜福利在线观看吧| 欧美激情久久久久久爽电影 | 高潮久久久久久久久久久不卡| 色综合亚洲欧美另类图片| √禁漫天堂资源中文www| 精品国产一区二区三区四区第35| 国产一区二区三区在线臀色熟女| 亚洲av片天天在线观看| 叶爱在线成人免费视频播放| 香蕉久久夜色| 香蕉丝袜av| 91在线观看av| 免费女性裸体啪啪无遮挡网站| 香蕉久久夜色| 欧美日韩乱码在线| 成人国产综合亚洲| 日本 欧美在线| 欧美 亚洲 国产 日韩一| 成在线人永久免费视频| 熟妇人妻久久中文字幕3abv| 黑丝袜美女国产一区| netflix在线观看网站| 亚洲狠狠婷婷综合久久图片| 成人永久免费在线观看视频| 亚洲七黄色美女视频| 亚洲最大成人中文| 亚洲国产欧美一区二区综合| 亚洲欧美日韩另类电影网站| 色精品久久人妻99蜜桃| 在线观看免费视频日本深夜| 欧美乱妇无乱码| 亚洲精品中文字幕在线视频| 免费在线观看黄色视频的| 老司机靠b影院| 一区在线观看完整版| 每晚都被弄得嗷嗷叫到高潮| 国产成人一区二区三区免费视频网站| 国产av在哪里看| 桃色一区二区三区在线观看| 国产又爽黄色视频| svipshipincom国产片| 精品国产美女av久久久久小说| 天堂影院成人在线观看| 国产精品二区激情视频| 看片在线看免费视频| 天堂√8在线中文| 中文字幕av电影在线播放| av有码第一页| av电影中文网址| 久久人妻福利社区极品人妻图片| 久久精品亚洲熟妇少妇任你| 中亚洲国语对白在线视频| 麻豆成人av在线观看| 日韩高清综合在线| 12—13女人毛片做爰片一| 69精品国产乱码久久久| av天堂在线播放| 在线av久久热| 国产xxxxx性猛交| 成熟少妇高潮喷水视频| 欧美午夜高清在线| 免费在线观看日本一区| 久久久国产精品麻豆| 韩国av一区二区三区四区| 精品国产乱码久久久久久男人| 两性午夜刺激爽爽歪歪视频在线观看 | 身体一侧抽搐| 国产极品粉嫩免费观看在线| 日本五十路高清| 国产97色在线日韩免费| 99re在线观看精品视频| av天堂久久9| 天堂动漫精品| 亚洲av成人一区二区三| 欧美日韩亚洲综合一区二区三区_| 一本大道久久a久久精品| 亚洲欧美一区二区三区黑人| 国语自产精品视频在线第100页| 一级片免费观看大全| 国产精品免费视频内射| 高清毛片免费观看视频网站| 悠悠久久av| 亚洲人成电影免费在线| 精品电影一区二区在线| 亚洲免费av在线视频| 亚洲国产看品久久| 天堂影院成人在线观看| 亚洲,欧美精品.| 国产亚洲欧美精品永久| 国产野战对白在线观看| 18禁美女被吸乳视频| 国产精品二区激情视频| 69精品国产乱码久久久| 亚洲国产精品合色在线| 国产欧美日韩精品亚洲av| 精品久久蜜臀av无| 国产一级毛片七仙女欲春2 | 俄罗斯特黄特色一大片| 午夜日韩欧美国产| 精品久久久久久久毛片微露脸| svipshipincom国产片| 少妇被粗大的猛进出69影院| 性欧美人与动物交配| 丝袜在线中文字幕| 午夜免费激情av| 女同久久另类99精品国产91| 这个男人来自地球电影免费观看| 在线观看66精品国产| 51午夜福利影视在线观看| 亚洲自偷自拍图片 自拍| 欧美成人性av电影在线观看| 午夜福利欧美成人| 大型av网站在线播放| www日本在线高清视频| 国产亚洲欧美在线一区二区| 亚洲欧美日韩高清在线视频| 亚洲,欧美精品.| 精品久久蜜臀av无| 搡老岳熟女国产| 亚洲av电影不卡..在线观看| 波多野结衣av一区二区av| 桃色一区二区三区在线观看| 国产精品久久久久久亚洲av鲁大| 日本撒尿小便嘘嘘汇集6| 午夜免费成人在线视频| 久久午夜亚洲精品久久| 国产精品亚洲一级av第二区| 一区二区三区精品91| 亚洲中文字幕日韩| 欧美乱妇无乱码| 日本黄色视频三级网站网址| 麻豆av在线久日| 国产精品1区2区在线观看.| 久久天躁狠狠躁夜夜2o2o| www日本在线高清视频| 亚洲av第一区精品v没综合| 精品卡一卡二卡四卡免费| 国产av精品麻豆| 两性午夜刺激爽爽歪歪视频在线观看 | 一进一出抽搐动态| 久久人人97超碰香蕉20202| 人人妻人人澡人人看| 亚洲av成人不卡在线观看播放网| 免费av毛片视频| 午夜久久久在线观看| 99国产极品粉嫩在线观看| x7x7x7水蜜桃| 欧美日本亚洲视频在线播放| 一进一出抽搐动态| 久久人人97超碰香蕉20202| 精品一区二区三区四区五区乱码| 波多野结衣高清无吗| 日本撒尿小便嘘嘘汇集6| 亚洲一区二区三区色噜噜| 色播亚洲综合网| 午夜日韩欧美国产| 免费高清在线观看日韩| 在线观看舔阴道视频| 香蕉丝袜av| 天天一区二区日本电影三级 | 天天一区二区日本电影三级 | 成年版毛片免费区| av网站免费在线观看视频| 亚洲色图av天堂| 99riav亚洲国产免费| 国产亚洲av高清不卡| 禁无遮挡网站| 悠悠久久av| 午夜老司机福利片| 村上凉子中文字幕在线| 99久久综合精品五月天人人| 亚洲国产看品久久| 丰满人妻熟妇乱又伦精品不卡| 亚洲国产看品久久| 国产成人啪精品午夜网站| 亚洲视频免费观看视频| 欧美大码av| 久久精品影院6| 波多野结衣巨乳人妻| 色综合站精品国产| 99国产精品99久久久久| videosex国产| 香蕉久久夜色| 精品免费久久久久久久清纯| 免费无遮挡裸体视频| 一区二区三区精品91| 久久久久久人人人人人| 亚洲成a人片在线一区二区| 国产真人三级小视频在线观看| 免费久久久久久久精品成人欧美视频| 婷婷六月久久综合丁香| 亚洲情色 制服丝袜| 久久精品影院6| 91精品国产国语对白视频| 天天添夜夜摸| 国产精品亚洲美女久久久| 男女之事视频高清在线观看| 成人国产综合亚洲| 麻豆av在线久日| 国产成人欧美| 日韩一卡2卡3卡4卡2021年| 亚洲国产精品999在线| 色综合欧美亚洲国产小说| 国产精品久久久久久精品电影 | 亚洲中文日韩欧美视频| 亚洲av日韩精品久久久久久密| 久久久久久久久中文| 久久人妻熟女aⅴ| 久久精品亚洲熟妇少妇任你| av超薄肉色丝袜交足视频| 亚洲人成网站在线播放欧美日韩| 久9热在线精品视频| 悠悠久久av| 亚洲免费av在线视频| 中文字幕人成人乱码亚洲影| 亚洲第一欧美日韩一区二区三区| 国产在线观看jvid| √禁漫天堂资源中文www| 国产主播在线观看一区二区| 首页视频小说图片口味搜索| 久热这里只有精品99| 亚洲情色 制服丝袜| 成人免费观看视频高清| 少妇熟女aⅴ在线视频| 69av精品久久久久久| 国产乱人伦免费视频| 亚洲色图 男人天堂 中文字幕| 神马国产精品三级电影在线观看 | 日韩精品免费视频一区二区三区| 多毛熟女@视频| 99国产精品免费福利视频| 麻豆一二三区av精品| 亚洲中文av在线| 午夜福利影视在线免费观看| 69av精品久久久久久| 午夜激情av网站| 禁无遮挡网站| 老熟妇仑乱视频hdxx| 18禁黄网站禁片午夜丰满| 国产欧美日韩一区二区三区在线| 99久久国产精品久久久| 亚洲美女黄片视频| 法律面前人人平等表现在哪些方面| 99精品欧美一区二区三区四区| 两个人看的免费小视频| 亚洲av成人av| 亚洲中文日韩欧美视频| 叶爱在线成人免费视频播放| 一级a爱视频在线免费观看| 后天国语完整版免费观看| 国产亚洲精品第一综合不卡| 国产精品久久电影中文字幕| 国产欧美日韩精品亚洲av| 中文字幕精品免费在线观看视频| 一本久久中文字幕| www.自偷自拍.com| 欧美日韩亚洲国产一区二区在线观看| 男女床上黄色一级片免费看| 又黄又粗又硬又大视频| 国产精华一区二区三区| 国产三级在线视频| 久久影院123| 国产精品av久久久久免费| 久久香蕉精品热| 夜夜夜夜夜久久久久| 宅男免费午夜| 99久久99久久久精品蜜桃| 1024视频免费在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲国产精品久久男人天堂| 女性生殖器流出的白浆| 国产精品一区二区精品视频观看| 可以在线观看的亚洲视频| 亚洲精品久久成人aⅴ小说| 免费一级毛片在线播放高清视频 |