• <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基因突變檢測
    亚洲五月天丁香| 国产黄片美女视频| 久久久久久久午夜电影| 人妻丰满熟妇av一区二区三区| 午夜福利18| 大又大粗又爽又黄少妇毛片口| 69人妻影院| 女同久久另类99精品国产91| 美女cb高潮喷水在线观看| 成年女人看的毛片在线观看| 国产精品一区www在线观看 | av福利片在线观看| 成年免费大片在线观看| 国产亚洲91精品色在线| 看十八女毛片水多多多| 久久精品国产亚洲av天美| 国产成人福利小说| 亚洲色图av天堂| 很黄的视频免费| 久久久成人免费电影| 日韩大尺度精品在线看网址| 91久久精品国产一区二区成人| 亚洲一区二区三区色噜噜| 欧美另类亚洲清纯唯美| 亚洲av日韩精品久久久久久密| 黄色一级大片看看| av女优亚洲男人天堂| 自拍偷自拍亚洲精品老妇| 免费av观看视频| 欧美另类亚洲清纯唯美| 久久人人爽人人爽人人片va| 少妇熟女aⅴ在线视频| 在线观看午夜福利视频| eeuss影院久久| 赤兔流量卡办理| 国产欧美日韩一区二区精品| ponron亚洲| 午夜久久久久精精品| 日日干狠狠操夜夜爽| 午夜福利在线观看免费完整高清在 | 能在线免费观看的黄片| 成人特级黄色片久久久久久久| 热99re8久久精品国产| 97碰自拍视频| 亚洲人成网站高清观看| 99riav亚洲国产免费| 久久精品人妻少妇| 中文字幕高清在线视频| 国产麻豆成人av免费视频| 舔av片在线| 深夜a级毛片| 在线观看舔阴道视频| 18禁在线播放成人免费| 99热这里只有是精品50| 此物有八面人人有两片| 一个人免费在线观看电影| 精品久久久久久久人妻蜜臀av| 色噜噜av男人的天堂激情| 床上黄色一级片| 熟女人妻精品中文字幕| 国产精品一区二区三区四区免费观看 | 国产一级毛片七仙女欲春2| 久久久色成人| 欧美性感艳星| 精华霜和精华液先用哪个| 欧美不卡视频在线免费观看| 啦啦啦观看免费观看视频高清| 亚洲成人久久爱视频| 综合色av麻豆| 搡老岳熟女国产| av天堂中文字幕网| 国产在视频线在精品| 美女 人体艺术 gogo| 亚洲av一区综合| 日韩精品有码人妻一区| 亚洲人成网站在线播放欧美日韩| 国产亚洲精品综合一区在线观看| 性欧美人与动物交配| 亚洲七黄色美女视频| 国产av不卡久久| 国产精品久久久久久亚洲av鲁大| 老司机深夜福利视频在线观看| 色尼玛亚洲综合影院| 国产精品久久视频播放| 色噜噜av男人的天堂激情| 久久精品91蜜桃| 中文资源天堂在线| 国产又黄又爽又无遮挡在线| 国产精品一区www在线观看 | 丰满乱子伦码专区| 51国产日韩欧美| 国内久久婷婷六月综合欲色啪| 无人区码免费观看不卡| 国产精品一及| 久久人人精品亚洲av| 亚洲成人久久性| 久久精品91蜜桃| 久久久久久九九精品二区国产| 亚洲经典国产精华液单| 免费av观看视频| 国产在线精品亚洲第一网站| 日本黄色片子视频| 久99久视频精品免费| 夜夜爽天天搞| а√天堂www在线а√下载| 久久久久久国产a免费观看| 国产伦精品一区二区三区视频9| 日日啪夜夜撸| 亚洲av五月六月丁香网| 国产精品亚洲美女久久久| 联通29元200g的流量卡| 亚洲国产精品合色在线| 成人美女网站在线观看视频| 国内久久婷婷六月综合欲色啪| 国产精品美女特级片免费视频播放器| 美女cb高潮喷水在线观看| 免费大片18禁| 精品免费久久久久久久清纯| 五月玫瑰六月丁香| 久久精品国产自在天天线| av在线观看视频网站免费| 色综合色国产| 真人做人爱边吃奶动态| 国产v大片淫在线免费观看| 久久天躁狠狠躁夜夜2o2o| 国产精品三级大全| 久久婷婷人人爽人人干人人爱| 日本色播在线视频| 亚洲性久久影院| x7x7x7水蜜桃| 久久久久九九精品影院| 国产精品爽爽va在线观看网站| 久久人人精品亚洲av| 国产一区二区三区视频了| 国产成人一区二区在线| 欧美国产日韩亚洲一区| 热99re8久久精品国产| 99riav亚洲国产免费| 免费不卡的大黄色大毛片视频在线观看 | 悠悠久久av| 亚洲在线观看片| 国产真实乱freesex| 男女之事视频高清在线观看| 国产精品一及| av中文乱码字幕在线| 久久精品91蜜桃| 国产伦在线观看视频一区| 国产精品野战在线观看| 极品教师在线视频| 国产人妻一区二区三区在| 欧美激情在线99| 中文字幕熟女人妻在线| 久久精品国产清高在天天线| 22中文网久久字幕| 亚洲成av人片在线播放无| 国产综合懂色| 桃色一区二区三区在线观看| 成人无遮挡网站| 变态另类丝袜制服| 久久中文看片网| 精华霜和精华液先用哪个| 老熟妇仑乱视频hdxx| 日韩精品中文字幕看吧| 在线播放无遮挡| 日日夜夜操网爽| 久久久精品大字幕| 一本精品99久久精品77| 91精品国产九色| 两个人视频免费观看高清| 欧美三级亚洲精品| 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看| 人妻久久中文字幕网| 3wmmmm亚洲av在线观看| 小蜜桃在线观看免费完整版高清| 噜噜噜噜噜久久久久久91| 22中文网久久字幕| 别揉我奶头 嗯啊视频| 国产精品国产三级国产av玫瑰| 免费人成在线观看视频色| 国产麻豆成人av免费视频| 精品久久久久久久久av| 三级国产精品欧美在线观看| 欧美成人免费av一区二区三区| 亚洲国产色片| 搡老妇女老女人老熟妇| 欧美日韩中文字幕国产精品一区二区三区| 久久精品国产亚洲av香蕉五月| 少妇人妻一区二区三区视频| 久久欧美精品欧美久久欧美| 亚洲自拍偷在线| 国产午夜精品久久久久久一区二区三区 | 99九九线精品视频在线观看视频| 成人三级黄色视频| 欧美性猛交╳xxx乱大交人| 18禁在线播放成人免费| 国产精品精品国产色婷婷| 观看美女的网站| 亚洲av熟女| 色播亚洲综合网| 免费观看的影片在线观看| 欧美黑人巨大hd| 午夜福利在线观看吧| 嫁个100分男人电影在线观看| 久久久久国内视频| 日韩亚洲欧美综合| 久久九九热精品免费| 亚洲国产色片| 国产美女午夜福利| 成人鲁丝片一二三区免费| 精品日产1卡2卡| 久久精品国产亚洲av天美| 高清在线国产一区| 国内精品久久久久精免费| 成人一区二区视频在线观看| 波多野结衣巨乳人妻| 日韩欧美在线乱码| 日日摸夜夜添夜夜添av毛片 | 国产午夜精品久久久久久一区二区三区 | 国产欧美日韩精品亚洲av| 在线观看av片永久免费下载| 少妇的逼水好多| 国产欧美日韩一区二区精品| 国产亚洲91精品色在线| 亚洲专区中文字幕在线| 欧美激情国产日韩精品一区| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲av嫩草精品影院| 久久久久久久久大av| 999久久久精品免费观看国产| 天堂av国产一区二区熟女人妻| 狂野欧美白嫩少妇大欣赏| 九九热线精品视视频播放| 无遮挡黄片免费观看| 国产精品1区2区在线观看.| 国产精品一区二区免费欧美| 热99re8久久精品国产| 蜜桃亚洲精品一区二区三区| 精品日产1卡2卡| 老女人水多毛片| 午夜精品一区二区三区免费看| 亚洲在线观看片| av在线亚洲专区| 国产精品98久久久久久宅男小说| 欧美三级亚洲精品| 国产免费一级a男人的天堂| 自拍偷自拍亚洲精品老妇| 人人妻,人人澡人人爽秒播| 最新中文字幕久久久久| 丰满人妻一区二区三区视频av| 久久精品影院6| 成年女人永久免费观看视频| 久久久久久久午夜电影| 一级毛片久久久久久久久女| 久久久久久久精品吃奶| 亚洲精品久久国产高清桃花| av女优亚洲男人天堂| 在线看三级毛片| 日本三级黄在线观看| 国产aⅴ精品一区二区三区波| 黄色欧美视频在线观看| 国产伦精品一区二区三区四那| 伦精品一区二区三区| 网址你懂的国产日韩在线| 亚洲中文字幕一区二区三区有码在线看| 精品一区二区三区人妻视频| 日本五十路高清| 久久婷婷人人爽人人干人人爱| 男女做爰动态图高潮gif福利片| 国产男人的电影天堂91| 在线观看美女被高潮喷水网站| 日本 欧美在线| 国产成人aa在线观看| 日日摸夜夜添夜夜添小说| 欧美激情久久久久久爽电影| 人妻制服诱惑在线中文字幕| 九九在线视频观看精品| 99国产极品粉嫩在线观看| eeuss影院久久| 99久久中文字幕三级久久日本| 女生性感内裤真人,穿戴方法视频| 两个人视频免费观看高清| 久久6这里有精品| 欧美成人一区二区免费高清观看| 亚洲av免费高清在线观看| 国产视频内射| 国产又黄又爽又无遮挡在线| 99热精品在线国产| 国产探花极品一区二区| 12—13女人毛片做爰片一| 又爽又黄a免费视频| 国产aⅴ精品一区二区三区波| 中文字幕久久专区| 国产精品野战在线观看| 看片在线看免费视频| 中国美白少妇内射xxxbb| 亚洲国产日韩欧美精品在线观看| 亚洲经典国产精华液单| 色综合婷婷激情| 色视频www国产| 日本欧美国产在线视频| 99热只有精品国产| 99国产极品粉嫩在线观看| 天天一区二区日本电影三级| 久久久久精品国产欧美久久久| 日韩亚洲欧美综合| 露出奶头的视频| 亚洲美女黄片视频| 欧美不卡视频在线免费观看| 日本一本二区三区精品| 麻豆国产97在线/欧美| 国产日本99.免费观看| 少妇高潮的动态图| 国产成人影院久久av| 又紧又爽又黄一区二区| 国产高清不卡午夜福利| 久久欧美精品欧美久久欧美| 精品一区二区免费观看| 亚洲av一区综合| 1024手机看黄色片| 老女人水多毛片| 内射极品少妇av片p| 人妻夜夜爽99麻豆av| 午夜a级毛片| 综合色av麻豆| 精品人妻熟女av久视频| 久久精品国产亚洲av天美| 国产亚洲精品av在线| 国产一级毛片七仙女欲春2| 少妇丰满av| 亚洲最大成人av| 最近在线观看免费完整版| 亚洲图色成人| 欧美绝顶高潮抽搐喷水| av在线蜜桃| 一区福利在线观看| 97超视频在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 无人区码免费观看不卡| 真实男女啪啪啪动态图| av福利片在线观看| 一a级毛片在线观看| 国产精品爽爽va在线观看网站| 一本久久中文字幕| 女的被弄到高潮叫床怎么办 | 草草在线视频免费看| 欧美性猛交╳xxx乱大交人| 亚洲人与动物交配视频| 中文在线观看免费www的网站| 国产毛片a区久久久久| 午夜福利高清视频| avwww免费| 欧美性猛交╳xxx乱大交人| 成人性生交大片免费视频hd| 欧美另类亚洲清纯唯美| 欧美日韩国产亚洲二区| 成年女人看的毛片在线观看| 搡女人真爽免费视频火全软件 | 国产精品野战在线观看| 婷婷精品国产亚洲av| 午夜免费成人在线视频| 国产国拍精品亚洲av在线观看| 性色avwww在线观看| 欧美成人免费av一区二区三区| 成人欧美大片| 欧美最新免费一区二区三区| 亚洲真实伦在线观看| 欧美黑人欧美精品刺激| 久久6这里有精品| 亚洲国产色片| 亚洲在线自拍视频| 久久香蕉精品热| 国产高清激情床上av| 免费观看人在逋| 久久精品国产亚洲网站| 可以在线观看毛片的网站| 国产亚洲精品久久久久久毛片| 99久久精品热视频| 国产精品精品国产色婷婷| 欧美3d第一页| 亚洲午夜理论影院| 免费观看在线日韩| 国产免费av片在线观看野外av| 亚洲欧美日韩卡通动漫| 精品午夜福利视频在线观看一区| av在线亚洲专区| 2021天堂中文幕一二区在线观| 成人特级av手机在线观看| 一本一本综合久久| 精品久久久久久久久久免费视频| 亚洲第一区二区三区不卡| 十八禁网站免费在线| 99久久精品国产国产毛片| 真实男女啪啪啪动态图| 日韩中文字幕欧美一区二区| 久久久久久久久久黄片| 亚洲四区av| 亚洲在线自拍视频| 人妻夜夜爽99麻豆av| 如何舔出高潮| 伊人久久精品亚洲午夜| 欧美bdsm另类| 性色avwww在线观看| av在线老鸭窝| 亚洲中文字幕一区二区三区有码在线看| 欧美精品国产亚洲| x7x7x7水蜜桃| 亚洲性久久影院| 97人妻精品一区二区三区麻豆| 99国产精品一区二区蜜桃av| 日本在线视频免费播放| 女同久久另类99精品国产91| 国产私拍福利视频在线观看| 成年人黄色毛片网站| 18+在线观看网站| 男人舔奶头视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲欧美日韩高清专用| 亚洲人成网站在线播| 在线a可以看的网站| 久久久久久伊人网av| 中国美女看黄片| 成熟少妇高潮喷水视频| 婷婷精品国产亚洲av| 亚洲综合色惰| 久久久久久久精品吃奶| 亚洲精品乱码久久久v下载方式| 尾随美女入室| 男女啪啪激烈高潮av片| 国产精品女同一区二区软件 | 日韩精品中文字幕看吧| 在线观看一区二区三区| 欧美在线一区亚洲| 干丝袜人妻中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久久久久久久亚洲 | 成年女人毛片免费观看观看9| 看片在线看免费视频| 日韩人妻高清精品专区| 亚洲国产欧洲综合997久久,| 99riav亚洲国产免费| 免费看日本二区| 老司机深夜福利视频在线观看| 亚洲美女黄片视频| 午夜亚洲福利在线播放| 少妇人妻精品综合一区二区 | 人人妻人人看人人澡| av.在线天堂| 深夜精品福利| 亚洲成人久久爱视频| 看黄色毛片网站| 日韩一本色道免费dvd| 99九九线精品视频在线观看视频| 免费看av在线观看网站| 国产精品国产高清国产av| 黄色女人牲交| 午夜福利在线观看吧| 中文字幕人妻熟人妻熟丝袜美| 国产精品98久久久久久宅男小说| 草草在线视频免费看| 1024手机看黄色片| 精品国内亚洲2022精品成人| 99久久久亚洲精品蜜臀av| 欧美性猛交黑人性爽| 麻豆av噜噜一区二区三区| 欧美日韩国产亚洲二区| 男女下面进入的视频免费午夜| 欧美成人一区二区免费高清观看| 长腿黑丝高跟| 日韩精品有码人妻一区| 观看美女的网站| 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 午夜亚洲福利在线播放| 久久草成人影院| 久久久久久久久中文| 欧美一区二区精品小视频在线| 国产精品三级大全| 国产精品久久视频播放| 亚洲人成网站高清观看| 神马国产精品三级电影在线观看| 制服丝袜大香蕉在线| 欧美成人免费av一区二区三区| 在线天堂最新版资源| 国产成人福利小说| 精品午夜福利视频在线观看一区| 全区人妻精品视频| 欧美日韩中文字幕国产精品一区二区三区| 国产v大片淫在线免费观看| 免费av不卡在线播放| 日韩强制内射视频| 精品国内亚洲2022精品成人| 九九在线视频观看精品| 极品教师在线免费播放| 精品日产1卡2卡| 国产老妇女一区| 色综合站精品国产| 最新在线观看一区二区三区| 成人国产综合亚洲| 久久人人爽人人爽人人片va| 99热这里只有是精品在线观看| 99视频精品全部免费 在线| 成年女人永久免费观看视频| 女的被弄到高潮叫床怎么办 | 国产精品女同一区二区软件 | 亚洲精品亚洲一区二区| 午夜亚洲福利在线播放| 又黄又爽又刺激的免费视频.| 午夜福利18| 午夜精品在线福利| 亚洲最大成人中文| 亚洲无线在线观看| 男女那种视频在线观看| 日本三级黄在线观看| 人人妻人人澡欧美一区二区| 毛片女人毛片| 伦理电影大哥的女人| 国产一区二区三区在线臀色熟女| 人妻久久中文字幕网| 不卡一级毛片| 久久亚洲真实| 午夜日韩欧美国产| 亚洲av成人av| 蜜桃久久精品国产亚洲av| 精品人妻偷拍中文字幕| 99久久无色码亚洲精品果冻| 熟妇人妻久久中文字幕3abv| 亚州av有码| 极品教师在线视频| 中文字幕av成人在线电影| 桃红色精品国产亚洲av| 日韩精品中文字幕看吧| 国产精品爽爽va在线观看网站| 精品久久久久久久久久久久久| 美女大奶头视频| 国产麻豆成人av免费视频| 日韩精品中文字幕看吧| 乱人视频在线观看| 成人一区二区视频在线观看| 99久久精品热视频| 婷婷丁香在线五月| av在线观看视频网站免费| 人妻少妇偷人精品九色| 亚洲人成网站在线播放欧美日韩| 校园人妻丝袜中文字幕| 看片在线看免费视频| 久久久久久久久久黄片| 内地一区二区视频在线| 亚洲成人中文字幕在线播放| 久久精品国产亚洲网站| 91久久精品电影网| 久久国内精品自在自线图片| 亚洲国产精品合色在线| 乱人视频在线观看| 波多野结衣巨乳人妻| 国产国拍精品亚洲av在线观看| 中文在线观看免费www的网站| avwww免费| 亚洲一区高清亚洲精品| 亚洲 国产 在线| 别揉我奶头~嗯~啊~动态视频| .国产精品久久| 噜噜噜噜噜久久久久久91| 国产又黄又爽又无遮挡在线| 久久这里只有精品中国| 国产69精品久久久久777片| 成年女人看的毛片在线观看| 啪啪无遮挡十八禁网站| 黄色配什么色好看| 超碰av人人做人人爽久久| 99精品在免费线老司机午夜| 国产伦精品一区二区三区视频9| 夜夜夜夜夜久久久久| 身体一侧抽搐| 免费观看的影片在线观看| 国产黄片美女视频| 欧美性猛交╳xxx乱大交人| 国产激情偷乱视频一区二区| 亚洲18禁久久av| 久久人妻av系列| 91麻豆精品激情在线观看国产| 亚洲人成网站在线播| 日本与韩国留学比较| 在线免费观看不下载黄p国产 | 日日干狠狠操夜夜爽| 亚洲无线观看免费| 我要看日韩黄色一级片| 日本一本二区三区精品| 黄色女人牲交| 亚洲avbb在线观看| 午夜a级毛片| 色av中文字幕| 18禁黄网站禁片免费观看直播| 国产单亲对白刺激| 村上凉子中文字幕在线| 久久久久性生活片| 欧美日本亚洲视频在线播放| aaaaa片日本免费| 日韩精品青青久久久久久| 午夜激情福利司机影院| ponron亚洲| 亚洲四区av| 1000部很黄的大片| 毛片女人毛片| 免费在线观看影片大全网站| 俺也久久电影网| 女同久久另类99精品国产91| 欧美xxxx黑人xx丫x性爽| 婷婷精品国产亚洲av| 免费搜索国产男女视频| 少妇高潮的动态图| 嫩草影院新地址| 国产免费男女视频| 嫩草影院精品99| 天堂av国产一区二区熟女人妻|