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

    茶樹PPR基因家族全基因組鑒定及白化相關基因表達分析

    2021-04-20 02:28:56劉丁丁王君雅湯榕津陳亮馬春雷
    茶葉科學 2021年2期
    關鍵詞:白化葉綠體結構域

    劉丁丁,王君雅,湯榕津,陳亮,馬春雷

    茶樹基因家族全基因組鑒定及白化相關基因表達分析

    劉丁丁1,2,王君雅1,2,湯榕津1,2,陳亮1*,馬春雷1*

    1. 中國農(nóng)業(yè)科學院茶葉研究所/國家茶樹改良中心/農(nóng)業(yè)部茶樹生物學與資源利用重點實驗室,浙江 杭州 310008; 2. 中國農(nóng)業(yè)科學院研究生院,北京 100081

    三角狀五肽(PPR)蛋白作為一類靶向定位于半自主細胞器的序列特異性RNA結合蛋白,在植物的生長發(fā)育過程中具有重要的作用。本研究基于茶樹基因組數(shù)據(jù),利用生物信息學方法對茶樹基因家族進行系統(tǒng)鑒定,并對其蛋白質(zhì)理化性質(zhì)、結構域保守性、亞細胞定位、基因結構、染色體定位與分布等進行了分析。結果表明,在茶樹基因組數(shù)據(jù)中共鑒定到858個基因,分屬于P和PLS兩個亞家族;結構域分析表明各個結構域在茶樹中比較保守;亞細胞定位結果顯示超過一半的CsPPR蛋白定位于葉綠體;基因結構分析表明茶樹基因家族中共有31%的基因不具有內(nèi)含子結構,并且家族成員在進化過程中發(fā)生過多次基因復制事件。另外,為探究茶樹基因家族在調(diào)控茶樹白化相關基因表達過程中的作用,對正常葉色品種舒茶早和安吉黃茶等5個白化茶樹品種進行轉(zhuǎn)錄組測序,共篩選到24個差異共表達基因,并利用實時熒光定量PCR技術對其在不同品種和不同組織中的表達模式進行了分析,結果顯示大多數(shù)基因在新梢、成熟葉和莖中高表達,但在花和根中痕量表達。本研究結果可為茶樹家族成員的基因克隆和功能研究提供參考。

    茶樹;基因家族;轉(zhuǎn)錄組;生物信息學;實時熒光定量PCR

    三角狀五肽(Pentatricopeptide repeat,PPR)蛋白是一類在陸生植物中廣泛存在,數(shù)量龐大的蛋白家族。一般情況下,PPR蛋白由N端長度可變的信號序列、PPR重復基序單元和C端結構域組成;根據(jù)PPR基序的種類、數(shù)量和排列順序,可分為P亞家族和PLS亞家族,其中PLS亞家族蛋白根據(jù)C末端結構域的不同,又可分為PLS類、E類、E+類和DYW類4類蛋白質(zhì)[1-3]。研究表明,PPR蛋白大多定位于葉綠體或線粒體,通過參與調(diào)控RNA代謝過程,從而在植物生長發(fā)育、胚胎形成和調(diào)控細胞質(zhì)雄性不育等方面發(fā)揮作用[4-6]。除此之外,PPR蛋白也被鑒定到與植物的葉綠體發(fā)育密切相關[7-12]。如簡磊等[10]利用圖位克隆的方法在水稻第3號染色體InDel標記3G5和SG3之間定位到1個編碼PPR蛋白的基因,該基因在1?387位發(fā)生了1個堿基G的插入,使編碼的蛋白失去功能,產(chǎn)生白化表型。何鵬等[12]利用基因沉默技術篩選到1個與棉花葉色形成相關的基因,該基因沉默后,調(diào)控葉綠體基因轉(zhuǎn)錄的和基因的表達水平降低,最終導致棉花葉色變黃。目前,袁玲[13]在白葉1號中克隆了1個基因,并推測該基因可能負向調(diào)控葉綠體基因的表達,導致白化現(xiàn)象產(chǎn)生;其他關于茶樹中與調(diào)控葉綠體發(fā)育相關的基因的研究還鮮有報道。

    鑒于基因家族部分成員具有靶向調(diào)控葉綠體基因表達,進而影響葉綠體發(fā)育及光合作用的功能,本研究通過對正常綠色品種舒茶早和安吉黃茶等5個白化品種進行轉(zhuǎn)錄組測序,篩選到24個可能與茶樹白化表型相關的基因,并對其在不同品種和不同組織中的表達特性進行了檢測;同時結合舒茶早基因組數(shù)據(jù),利用生物信息學方法對茶樹基因家族成員進行系統(tǒng)鑒定,研究結果將為今后功能的深入挖掘提供幫助。

    1 材料與方法

    1.1 材料

    以保存于國家種質(zhì)杭州茶樹圃內(nèi)的白雞冠、舒茶早、安吉黃茶、景白1號、景白2號、中白1號和中黃3號為試驗材料,取其春季一芽二葉新梢,于液氮中速凍,后用于CsPPR基因序列分析和轉(zhuǎn)錄組測序分析。并以龍井43為試驗材料,分別取其一芽二葉、莖、成熟葉、花和根等組織器官,于液氮中速凍,用于組織表達模式分析。

    在茶樹基因組數(shù)據(jù)庫TPIA(http://tpia.teaplant.org/index.html)下載已測序完成的舒茶早2018版和2020版基因組序列、基因組注釋信息、預測CDS序列以及預測蛋白質(zhì)序列等數(shù)據(jù),用于茶樹基因家族的生物信息學分析。

    1.2 茶樹CsPPR基因家族篩選與鑒定

    在Pfam[14]數(shù)據(jù)庫下載PPR種子序列(PF01535),基于Linux系統(tǒng)利用隱馬爾可夫模型(profile hidden Markov models,profile HMMS)軟件包HMMER3.1b1[15-16]的搜索功能鑒定茶樹中的基因家族,定義閾值為1e-10;并利用NCBI的CDD數(shù)據(jù)庫分析基因的編碼蛋白特征,以確保家族成員僅具有PPR及其相關結構域。

    利用hmmpress程序格式化擬南芥基因家族的P(PPR)、S、L、L2、E、E+和DYW 7個保守結構域定義的HMMER矩陣,得到一個HMM數(shù)據(jù)庫;通過hmmscan程序檢索CsPPR蛋白家族候選序列的PPR結構域分布,定義閾值為1e-10,并人工逐條對候選序列進行結構分析和分類。

    1.3 茶樹CsPPR蛋白保守結構域分析

    基于CsPPR蛋白序列結構,提取CsPPR蛋白序列中7種PPR相關結構域,并分別利用ClustalW程序進行多序列比對,并利用WebLogo(http://weblogo.berkeley.edu/l ogo. cgi)在線軟件繪制出motif logo圖。

    1.4 茶樹CsPPR蛋白分析

    利用ProtParam在線軟件分析CsPPR的理化性質(zhì);利用SignalP-5.0 Server(http://www.cbs.dtu.dk/services/SignalP)在線軟件對CsPPR蛋白的信號肽進行分析;利用PSORT(https://wolfpsort.hgc.jp)在線軟件對CsPPR蛋白進行亞細胞定位預測;利用ClustalW程序包對CsPPR蛋白家族序列進行多序列比對,并通過MEGA 7.0軟件對多序列比對結果進行分析并構建序列進化樹,通過EvolView(http://120.202.110.254:8280/evolview)在線工具對進化樹進行可視化作圖。

    1.5 茶樹CsPPR家族的基因結構分析及染色體定位與分布

    基于舒茶早基因組數(shù)據(jù),利用GSDS在線軟件分析CsPPR家族成員的基因結構;并利用blast程序,分析茶樹基因家族成員的同源基因,之后利用Tbtools軟件基于基因注釋信息對同源關系進行可視化。

    1.6 茶樹部分CsPPR基因序列鑒定

    采用艾德萊植物RNA快速提取試劑盒和天根植物總DNA提取試劑盒提取白雞冠總RNA和總DNA,之后利用微量紫外檢測儀Nano Drop檢測RNA和DNA的濃度,并利用1.2%的瓊脂糖凝膠電泳檢測其質(zhì)量。采用天根Fastking gDNA Dispelling RT SuperMix kit將RNA反轉(zhuǎn)錄為cDNA。以白雞冠的DNA和cDNA為模板分別進行PCR擴增,PCR引物利用NCBI網(wǎng)站中的primer blast程序設計。PCR反應體系為:5?μL KOD buffer,5?μL dNTP,3?μL MgSO4,上下游PCR引物各1?μL,1?μL模板基因組DNA/cDNA,1?μL KOD酶,最后補加33?μL ddH2O至終體積50?μL。PCR反應程序為:94℃預變性2?min;98℃變性30?s;設置52~64℃退火溫度梯度,退火1?min,68℃延伸30?s,40個循環(huán);最后,68℃再延伸7?min,4℃保溫。利用1.2%的瓊脂糖凝膠電泳對PCR產(chǎn)物進行檢測,對檢測到含有單一目的條帶序列的PCR產(chǎn)物送浙江尚亞生物技術公司測序,驗證片段的準確性。

    1.7 不同茶樹品種中差異共表達的CsPPR分析

    本研究利用6個茶樹品種的轉(zhuǎn)錄組數(shù)據(jù),以正常綠色品種舒茶早為對照組,以白化品種安吉黃茶、景白1號、景白2號、中白1號和中黃3號為試驗組,篩選差異共表達的基因,差異基因的篩選條件為:表達差異倍數(shù)|log2Fold Change|>1(Fold change=試驗組/對照組)。根據(jù)獲得的表達數(shù)據(jù)利用R軟件包繪制韋恩圖。

    1.8 實時熒光定量PCR分析

    表1 熒光定量PCR引物序列

    2 結果與分析

    2.1 PPR基因家族的鑒定與分類

    近幾年,隨著測序技術的不斷進步,多個茶樹基因組草圖和染色體級別的茶樹參考基因組相繼完成[17-21],其中由安徽農(nóng)業(yè)大學于2018年和2020年分別發(fā)表的舒茶早基因組圖譜,因其較高的測序質(zhì)量,更具使用價值;兩個版本的基因組分別包含33?932個和50?525個蛋白編碼基因,為了進一步了解兩版注釋基因的差別,本研究對兩版注釋中的茶樹PPR家族進行了比較分析,分別鑒定得到676個和858個PPR蛋白序列(圖1)。隨后基于各序列中的PPR結構域及其相關結構域的分布情況,初步將PPR家族成員分成P和PLS兩個亞家族,其中PLS亞家族依據(jù)不同結構域的分布情況,可分為PLS類、E類、E+類和DYW類。初步分析表明,2020版基因組注釋得到的茶樹PPR亞家族均包含更多成員,其中P亞家族包含383個基因,PLS亞家族包含475個基因,將其分別命名為;而2018版注釋中P亞家和PLS亞家族分別包含331和345個基因家族成員。兩版注釋的差別主要來自P類、E+類和DYW類基因數(shù)量,2020版基因組鑒定到3類基因的數(shù)量較2018版基因組多50個左右,說明新版基因組對復雜結構域的注釋能力更強了。

    圖1 茶樹PPR基因家族分類及數(shù)目

    2.2 茶樹CsPPR家族的蛋白序列結構域保守性分析

    通過比較不同分類PPR蛋白的結構域序列特征,對茶樹CsPPR蛋白的保守性進行分析(圖2),結果顯示,茶樹CsPPR蛋白中的P、L、S、L2結構域序列與PF01535種子序列之間存在較高的序列相似性,特別是P結構域序列,與PF01535種子序列之間僅有5個堿基位點的差異;并且發(fā)現(xiàn)P、S結構域的保守性要高于L、L2結構域的保守性。進一步對每個結構域進行單獨分析,發(fā)現(xiàn)P和L結構域的第35位和L2結構域的第36位分別存在一個保守氨基酸位點,分別是Pro35、Ser35和Pro36;而E、E+和DYW結構域的一些堿基位點也分別具有較高的保守性,特別是DYW結構域序列末端具有3個保守氨基酸位點:Asp-Tyr-Trp,且這3種結構域序列長度大約為76、31、94個氨基酸殘基。

    2.3 茶樹CsPPR家族基因編碼蛋白的理化性質(zhì)分析

    對茶樹中的858個CsPPR蛋白進行理化性質(zhì)分析,結果表明,蛋白質(zhì)的序列長度為106~2?284個氨基酸;相對分子質(zhì)量最小的是CsPPR613,為11.68?kDa,最大的是CsPPR50,為261.97?kDa;理論等電點為4.48~9.93,其中有443個CsPPR蛋白的pI小于7,表明茶樹中有超過一半的CsPPR蛋白為酸性蛋白;親疏水性分析表明,CsPPR蛋白家族中有583條親水性蛋白,273條疏水性蛋白;不穩(wěn)定系數(shù)分析表明,在CsPPR蛋白家族中,有557個CsPPR蛋白是穩(wěn)定蛋白,301個屬于不穩(wěn)定蛋白。

    2.4 茶樹CsPPR家族基因編碼蛋白的信號肽分析

    信號肽是位于新合成肽鏈的N端,由15~30個氨基酸殘基組成的一段多肽,通常含有6~15個帶正電荷的非極性氨基酸,其作用是將新生肽鏈轉(zhuǎn)移到內(nèi)質(zhì)網(wǎng)上,具有可判定一個蛋白質(zhì)是否為分泌蛋白的功能,也與蛋白質(zhì)的定位有關[22]。858個CsPPR蛋白質(zhì)的信號肽預測結果表明,幾乎所有的CsPPR蛋白的Signal peptide預測分值都小于0.5,僅有4個CsPPR蛋白的預測分值大于0.5,分別是CsPPR22、CsPPR62、CsPPR261和CsPPR668。表明在CsPPR蛋白家族中僅有4個蛋白存在信號肽,屬于分泌蛋白,其他均屬于非分泌蛋白,合成后不進入內(nèi)質(zhì)網(wǎng)和高爾基體分泌途徑,而是直接釋放到細胞質(zhì)中。

    2.5 茶樹CsPPR家族基因編碼蛋白的亞細胞定位分析

    通過PSORT軟件對CsPPR基因家族的蛋白進行亞細胞定位預測,預測結果表明(表2),超過一半的CsPPR蛋白定位在了葉綠體上;其次是定位在細胞質(zhì)、線粒體和細胞核上;少數(shù)CsPPR蛋白定位在細胞質(zhì)膜、液泡和內(nèi)質(zhì)網(wǎng)等區(qū)域;個別CsPPR蛋白同時定位在葉綠體和線粒體,以及細胞質(zhì)與細胞核上。對各類型CsPPR蛋白的亞細胞定位預測結果單獨分析發(fā)現(xiàn),總體上各類型CsPPR蛋白在不同細胞器的分布情況與總的CsPPR蛋白具有相似的趨勢。總體來看,茶樹中的CsPPR蛋白以定位于葉綠體為主,少量定位于線粒體,且定位于半自主性細胞器的數(shù)量占CsPPR蛋白總數(shù)的一半以上,這與番茄、煙草、甘藍型油菜等植物的PPR蛋白定位結果存在一些差異,如譚暉等[23]利用Prot Comp 9.0軟件預測到甘藍型油菜中接近50%的PPR基因定位于線粒體。這說明PPR蛋白雖具有極高的分化程度,但是仍具有典型的定位于半自主性細胞器的特點。茶樹CsPPR蛋白大多定位于葉綠體,可能是因為茶樹屬于多年生常綠木本植物,且以葉用為主,其葉片進行光合作用的分子機制要復雜于其他植物。

    注:字母高度與氨基酸的保守性相關,字母高度越高,氨基酸越保守

    表2 茶樹CsPPR蛋白亞細胞定位預測

    2.6 茶樹CsPPR家族基因結構分析

    現(xiàn)有研究表明,陸生植物中的多數(shù)基因具有無內(nèi)含子的特點。為分析茶樹基因家族成員的序列結構,本研究利用茶樹基因內(nèi)含子、外顯子的注釋信息,以及家族序列比對結果,繪制基因的聚類樹(圖3)。結果表明,存在31%的基因不具有內(nèi)含子,而含有1~5個內(nèi)含子的基因數(shù)量最多,占整個基因家族的65%;少部分基因含有6~10個內(nèi)含子,僅14個基因含有內(nèi)含子的數(shù)量在11個以上。此外,對其非編碼區(qū)(Untranslation region,UTR)進行分析(圖3-C)表明,接近40%的基因不存在UTR區(qū)域;對基因家族成員的基因序列長度分析發(fā)現(xiàn),序列長度為0.4~50?kbp,存在很大的差異。

    2.7 CsPPR基因家族的染色體定位與分布

    將鑒定到的858條基因與2020年發(fā)表的染色體級舒茶早基因組進行比對分析,發(fā)現(xiàn)其中613條基因分布于15條染色體上(圖4-A),245條基因分布于211條contig序列上。15條染色體上包含的基因數(shù)量為20~74,其中12號染色體上分布數(shù)量最少,1號染色體上分布數(shù)量最多。隨后對分布于染色體上的基因家族成員進行共線性分析,發(fā)現(xiàn)不同染色體間具有高度同源性的基因有26對,而在染色體內(nèi)具有高度同源性的基因有128對(圖4-B),表明茶樹基因家族在進化過程中可能存在較多染色體內(nèi)部的基因復制事件。

    注:A:茶樹中有無內(nèi)含子的CsPPR基因分布(橘色,無內(nèi)含子;藍色,有內(nèi)含子);B:茶樹CsPPR基因家族中包含不同數(shù)量內(nèi)含子的CsPPR基因分布;C:茶樹中有無UTR區(qū)域的CsPPR基因分布(綠色,無UTR;紅色,有UTR)

    注:A. 茶樹各個染色體上PPR基因數(shù)量的分布;B. 茶樹各個染色體間PPR基因的共線性分析

    2.8 茶樹CsPPR基因家族成員的進化分析

    茶樹基因家族成員之間的系統(tǒng)進化結果表明(圖5),P和PLS兩個亞家族蛋白分別集中聚集成簇,但兩個亞家族內(nèi)均存在一定數(shù)量的離群序列,這在一定程度上反映了P亞家族和PLS亞家族成員在進化上的差異,但也表明了兩個亞家族成員之間存在基因交流現(xiàn)象。就PLS亞家族而言,僅DYW類蛋白分布相對集中,而E、E+和PLS類蛋白均分布較離散,這可能是與不同結構域之間的序列差異性和分化程度相關,也進一步反映了茶樹CsPPR家族蛋白序列結構的復雜性,可能正是這種序列結構的復雜性以及序列之間的差異性造就了基因家族功能的多樣性。

    2.9 部分CsPPR家族成員序列鑒定

    為進一步對生信分析結果的準確性進行驗證,我們隨機挑選了78個和55個基因,分別以DNA和cDNA為模板進行擴增。通過反復優(yōu)化PCR反應體系,最終成功擴增出76個基因的DNA序列和49個基因的cDNA序列(表3);測序結果表明,擴增序列與基因組數(shù)據(jù)庫中的序列一致,這說明目前發(fā)布的茶樹基因組拼接質(zhì)量較好,已能夠滿足常規(guī)的基因家族分析及基因克隆試驗。

    2.10 不同茶樹品種中差異CsPPR基因分析

    為進一步探究基因家族在調(diào)控茶樹白化相關基因表達過程中的作用,我們對正常葉色品種舒茶早和安吉黃茶等5個白化茶樹品種進行了轉(zhuǎn)錄組測序,測序質(zhì)量分析表明(表4),幾個品種的clean reads與參考基因組的比對率均在80%以上,其中比對到基因外顯子區(qū)域的reads數(shù)超過93%,Q30值在82.93%~89.51%,表明整體測序質(zhì)量較好、與參考基因組比對率高,符合后續(xù)分析要求。

    隨后,我們對5個組合的差異基因進行了初步分析(圖6),共篩選到379個差異表達的基因,其中24個基因在所有白化茶樹品種中差異共表達;在24個基因中下調(diào)表達的有20個,上調(diào)表達的有4個,推測這些基因可能與茶樹的白化表型相關。

    圖5 茶樹PPR基因家族的進化分析

    表3 部分CsPPR家族基因序列的PCR擴增結果

    表4 測序數(shù)據(jù)統(tǒng)計

    注:SCZ:舒茶早;AJHC:安吉黃茶;JB1:景白1號;JB2:景白2號;ZB1:中白1號;ZH3:中黃3號

    Note: SCZ: Shuchazao. AJHC: Anji Huangcha. JB1: Jingbai 1. JB2: Jingbai 2. ZB1: Zhongbai 1. ZH3: Zhonghuang 3

    2.11 實時熒光定量PCR分析

    為驗證轉(zhuǎn)錄組測序結果的可靠性,我們采用熒光定量PCR技術對篩選到的24個基因在6個茶樹品種的表達情況進行了檢測(圖7)。結果表明,120組熒光定量數(shù)據(jù)中,與轉(zhuǎn)錄組表達趨勢相同的有103組,占所有檢測數(shù)據(jù)的85.8%。其中17個基因在5個白化茶樹品種中的熒光定量結果與轉(zhuǎn)錄組表達趨勢相一致;和在4個白化茶樹品種中的熒光定量結果與轉(zhuǎn)錄組表達趨勢一致;而、和等5個基因僅在2個或3個品種的表達趨勢與轉(zhuǎn)錄組相一致,這可能是由于熒光定量PCR技術和轉(zhuǎn)錄組數(shù)據(jù)在檢測基因表達差異上的靈敏度不同造成的試驗誤差,因此上述差異在可接受范圍內(nèi),說明本文利用轉(zhuǎn)錄組數(shù)據(jù)篩選出的結果是可靠的。

    2.12 差異共表達的基因在茶樹不同組織中的表達模式分析

    為了進一步探究基因在茶樹生長發(fā)育過程中的潛在功能,本文分析了24個基因在5個茶樹組織中的表達模式。表達熱圖顯示(圖8),24個基因在茶樹各個組織中均有表達,但不同基因在不同組織中的表達量存在差異,且多數(shù)基因在茶樹新梢、成熟葉和莖中具有較高的表達量,但在花和根中痕量表達。其中基因在新梢中表達量最低;、、和在成熟葉中表達量較低;在莖中表達量較低;而、和在根中具有相對較高的表達量,表明這3個基因可能對茶樹根系具有重要調(diào)控作用。

    3 討論

    隨著多種植物基因組計劃的完成,從全基因組層面對基因家族進行鑒定和研究已成為現(xiàn)實。目前,除了在模式植物擬南芥、水稻、番茄、煙草和白楊中分別鑒定到441、491、471、487個和626個基因[1,24-26],在非模式植物甘藍型油菜、小麥和谷子中也已分別鑒定到1?079、1?351個和486個基因[23,27-28]。本研究基于已測序完成的2018版和2020版舒茶早茶樹基因組數(shù)據(jù)[18,20],利用生物信息學的方法分別篩選到676個和858個茶樹基因,均分屬于P和PLS兩個亞家族,其中PLS亞家族又包括PLS、E、E+和DYW 4類,該基因家族成員數(shù)量差異表明新版基因組對復雜結構域的注釋能力更強;亞細胞定位預測結果顯示,茶樹CsPPR蛋白主要定位于葉綠體上,少數(shù)定位于細胞質(zhì)、線粒體和細胞核中,表明茶樹CsPPR蛋白具有調(diào)控半自主細胞器中RNA代謝的功能;基因在染色體上的分布結果顯示,茶樹基因在所有的染色體上均有分布,且比較分散[27],這與小麥基因在染色體上的分布趨勢相一致,而與水稻基因在染色體上的分布有差異,水稻基因在染色體上的分布則比較集中[29]。系統(tǒng)進化分析發(fā)現(xiàn),家族成員中P類和DYW類分別聚集成簇,而PLS、E和E+類分布比較離散,這可能與基因的進化過程有關,為了更好的適應環(huán)境變化,基因之間分化程度較高,序列之間差異較大,說明基因家族成員的功能具有多樣性和復雜性;基因結構分析表明,茶樹中超過一半的基因不具有內(nèi)含子,這可能與基因在進化過程中發(fā)生了反轉(zhuǎn)錄轉(zhuǎn)座事件有關[3]。

    注:RNA-Seq數(shù)據(jù)是log2轉(zhuǎn)換后的fold change值;RT-qPCR數(shù)據(jù)是log2轉(zhuǎn)換后的基因相對表達量

    圖8 茶樹不同組織間基因的表達模式分析

    為了探究基因家族成員在調(diào)控茶樹白化相關基因表達過程中的作用,本研究從綠色品種舒茶早和安吉黃茶等5個白化品種的轉(zhuǎn)錄組數(shù)據(jù)中,篩選到24個差異共表達的基因;隨后對這24個基因在不同品種和不同組織中的表達情況進行分析,結果表明,多數(shù)基因在茶樹新梢、成熟葉和莖中具有較高的表達量,這與水稻中定位到的多個調(diào)控白化表型的基因表達模式一致,如[30]、[31]等基因均在葉片中具有較高的表達水平。此外,我們將這24個差異基因與擬南芥基因家族成員進行序列比對,發(fā)現(xiàn)、、、基因分別與擬南芥的、、、基因的同源性達到65%以上,表明它們可能具有相似的生物學功能?;蚨ㄎ挥谫|(zhì)體,通過與細胞器RNA編輯因子MORF2和MORF9相互作用,進而調(diào)控與葉綠體發(fā)育相關的質(zhì)體基因的表達,基因發(fā)生T-DNA插入導致擬南芥子葉色素缺失,進而出現(xiàn)子葉白化現(xiàn)象[32];基因定位于葉綠體,是擬南芥中第一個發(fā)現(xiàn)與早期胚胎致死相關的基因,其突變會導致細胞持續(xù)分裂和胚膨大,形成異常的組織結構,進而影響子葉發(fā)育[33];基因定位于葉綠體中,通過與葉綠體的操縱子相結合,從而調(diào)控葉綠體RNA代謝過程,該基因突變會影響細胞色素復合體的積累[34];是編輯葉綠體轉(zhuǎn)錄本和所需要的基因,該基因發(fā)生突變會導致葉綠體發(fā)育受損,幼苗呈現(xiàn)黃化表型,甚至死亡[35]。因此推測茶樹葉色白化變異也可能與某些基因的突變有關,但具體的分子機理,還需進一步的深入研究。

    [1] Claire L, Charles A, Sebastien A, et al. Genome-wide analysis ofpentatricopeptide repeat proteins reveals their essential role in organelle biogenesis [J]. The Plant Cell, 2004, 16(8): 2089-2103.

    [2] Small I D, Peeters N. The PPR motif-a TPR-related motif prevalent in plant organellar proteins [J]. Trends Biochem Sci, 2000, 25(2): 46-47.

    [3] Dipnarayan S, Prasad A M, Ramamurthy S. Pentatricopeptide repeat proteins and their emerging roles in plants [J]. Plant Physiology and Biochemistry, 2007, 45 (8): 521-534.

    [4] Barkan A, Waller M, Nolasco M. A nuclear mutation in maize blocks the processing and translation of several chloroplast mRNAs and provides evidence for the differential translation of alternative mRNA forms [J]. EMBO J, 1994, 13(13): 3170-3181.

    [5] Sam M. An overview of pentatricopeptide repeat proteins and their applications [J]. Bochimine, 2015, 113: 93-99.

    [6] Schmitz-Linneweber C, Small I. Pentatricopeptide repeat proteins: a socket for organelle gene expression [J]. Trends Plant Sci, 2008, 13(12): 663-670.

    [7] 陳龍. 水稻淡綠葉基因PGL12的克隆與功能分析[D]. 武漢: 華中農(nóng)業(yè)大學, 2019: 40-42. Chen L. Cloning and functional analysis of12 in rice [D]. Wuhan: Huazhong Agricultural University, 2019: 40-42.

    [8] 陳璐. 水稻低溫條件下葉綠體發(fā)育必需基因TCM8的克隆與功能分析[D]. 上海: 上海師范大學, 2019: 27-29. Chen L. Cloning and functional analysis of essential gene TCM8 for chloroplast development under low temperature conditions in rice () [D]. Shanghai: Shanghai Normal University, 2019: 27-29.

    [9] 葛生珍. 水稻黃化突變體xnt7的生理特性和基因精細定位[D]. 重慶: 西南大學, 2014: 28-33. Ge S Z. Physiological characteristics of a xanthic mutant and fine-mapping of its related genein rice [D]. Chongqing: Southwest University, 2014: 28-33.

    [10] 簡磊, 王仲康, 曾冬冬, 等. 水稻白化轉(zhuǎn)綠突變體albg的鑒定和基因精細定位[J]. 核農(nóng)學報, 2017, 31(12): 2289-2297.Jian L, Wang Z K, Zeng D D, et al. Characterization and fine mapping of a green-revertible albino (albg) mutant in rice [J]. Journal of Nuclear Agricultural Sciences, 2017, 31(12): 2289-2297.

    [11] 尚麗娜, 陳新龍, 米勝南, 等. 水稻溫敏型葉片白化轉(zhuǎn)綠突變體tsa2的表型鑒定與基因定位[J]. 作物學報, 2019, 45(5): 662-675. Shang L N, Chen X L, Mi S N, et al. Phenotypic identification and gene mapping of temperature-sensitive green revertible albino mutant tsa2 in rice (L.) [J]. Acta Agronomica Sinica, 2019, 45(5): 662-675.

    [12] 何鵬, 張麗華, 牟亞楠, 等. 棉花葉片黃化基因的篩選、鑒定及功能分析[C]//中國農(nóng)學會棉花分會2017年年會暨第九次會員代表大會論文匯編. 安陽: 棉花學報, 2017: 38. He P, Zhang L H, Mou Y N, et al. Dentification and functional analysis of yellow variegated leaf genein cotton [C]//2017 Annual Meeting of the Cotton Branch of Chinese Agricultural Association and the Ninth Member Congress. Anyang: Cotton Science, 2017: 38.

    [13] 袁玲. 白葉1號階段性返白過程中差異表達基因的分離及部分基因全長cDNA克隆[D]. 長沙: 湖南農(nóng)業(yè)大學, 2012: 29-39. Yuan L. Separation of differentially expressed genes and full length cDNA cloning of selective genes during periodic albinism in Anjibaicha() [D]. Changsha: Hunan Agricultural University, 2012: 29-39.

    [14] Punta M, Coggill P C, Eberhardt R Y. The pfam protein families database [J]. Nucleic Acids Research, 2012, 40(1): 290-301.

    [15] Simon C P, Aurelien L, Sean R. E, et al. HMMER web server: 2018 update [J]. Nucleic Acids Research, 2018, 46: 200-204.

    [16] Finn R D, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching [J]. Nucleic Acids Res, 2011, 39: 29-37.

    [17] Xia E H, Zhang H B, Sheng J, et al. The tea tree genome provides insights into tea flavor and independent evolution of caffeine biosynthesis [J]. Molecular Plant, 2017, 10(6): 866-877.

    [18] Wei C L, Yang H, Wang S B, et al. Draft genome sequence ofvarprovides insights into the evolution of the tea genome and tea quality [J]. PNAS, 2018, 115(18): E4151-E4158.

    [19] Chen J D, Zheng C, Ma J Q, et al. The chromosome-scale genome reveals the evolution and diversification after the recent tetraploidization ebent in tea plant [J]. Horticulture Research, 2020, 7: 63. doi: 10.1038/s41438-020-0288-2.

    [20] Xia E H, Tong W, Hou Y, et al. The reference genome of tea plant and resequencing of 18 diverse accessions provide insights into genome evolution and adaptation of tea plants [J]. Molecular Plant, 2020, 13(6/7): 1013-1026.

    [21] Zhang Q J, Li W, Li K, et al. The chromosome-level reference genome of tea tree unveils recent bursts of non-autonomous LTR retrotransposons to drive genome size evolution [J]. Molecular Plant, 2020, 13(7): 935-938.

    [22] 韋雪芳, 王冬梅, 劉思, 等. 信號肽及其在蛋白質(zhì)表達中的應用[J]. 生物技術通報, 2006(6): 38-42. Wei X F, Wang D M, Liu S, et al. Signal sequence and its application to protein expression [J]. Journal of Biotechnology, 2006(6): 38-42.

    [23] 譚暉, 官春云. 甘藍型油菜PPR家族生物信息學分析與新疆野生油菜候選育性基因克隆[J]. 作物研究, 2017, 31(3): 246-255. Tan H, Guan C Y. Bioinformatics analysis of PPR family inL. and cloning of candidate restorer gene of Xinjiang wild rapeseed [J]. Crop Research, 2006, 31(3): 246-255.

    [24] Nicholas O, Mitsuru H, Charles A, et al. On the expansion of the pentatricopeptide repeat gene family in plants [J]. Society of Molecular Biology and Evolution, 2008, 25: 1120-1128.

    [25] 丁安明. 番茄與煙草PPR基因家族分析及育性相關基因功能研究[D]. 北京: 中國農(nóng)業(yè)科學院, 2014. Ding A M. Identification of PPR gene family inand tomato and function analysis of Rf-related genes in tobacco [D]. Beijing: Chinese Academy of Agricultural Sciences, 2014.

    [26] Xing H T, Fu X K, Yang C. et al. Genome-wide investigation of pentatricopeptide repeat gene family in poplar and their expression analysis in response to biotic and abiotic stresses [J]. Scientific Reports, 2018, 8: 2817. doi: 10.1038/s41598-018-21269-1.

    [27] 郭彩娟, 公杰, 劉永杰, 等. 全基因組小麥PPR基因家族鑒定及表達分析[J]. 生物技術通報, 2019, 35(8): 1-8. Guo C J, Gong J, Liu Y J, et al. Genome-wide bioinformatics identification of PPR gene family and expression profiles analysis in wheat [J]. Journal of Biotechnology, 2019, 35(8): 1-8.

    [28] Liu J M, Xu Z S, Lu P P, et al. Genome-wide investigation and expression analyses of the pentatricopeptide repeat protein gene family in foxtail millet [J]. BMC Genomics. 2016, 17(1): 840. doi: 10.1186/s12864-016-3184-2.

    [29] Chen G, Zou Y, Hu J, et al. Genome-wide analysis of the rice PPR gene family and their expression profiles under different stress treatments [J]. BMC Genomics, 2018, 19(1): 720. doi: 10.1186/s12864-018-5088-9.

    [30] 倪晨子. 水稻PPR基因功能的初步分析[D]. 武漢: 武漢大學, 2017: 38-39. Ni C Z. An initial function analysis of the PPR genein riceL) [D]. Wuhan: Wuhan University, 2017: 38-39.

    [31] 高媛媛. 水稻PPR蛋白基因TCD34的克隆與功能分析[D]. 上海: 上海師范大學, 2020: 27-36. Gao Y Y. Cloning and functional analysis of PPR genein rice () [D]. Shanghai: Shanghai Normal University, 2020: 27-36.

    [32] Du L, Zhang J, Qu S F, et al. The pentratricopeptide repeat protein pigment-defective mutant2 is involved in the regulation of chloroplast development and chloroplast gene expression in[J]. Plant and Cell Physiology, 2017, 4(58): 747-759.

    [33] Cushing D A, Forsthoefel N R, Gestaut D R, et al. Arabidopsis emb175 and other ppr knockout mutants reveal essential roles for pentatricopeptide repeat (PPR) proteins in plant embryogenesis [J]. Planta, 2005, 221: 424-436.

    [34] Karin M, Susanne F, Takahiro N, et al. HCF152, an Arabidopsis RNA Binding Pentatricopeptide Repeat Protein Involved in the Processing of Chloroplast psbB-psbT-psbH-petB-petD RNAs [J]. Plant Cell, 2003, 15(6): 1480-1495.

    [35] Hein A, Knoop V. Expected and unexpected evolution of plant RNA editing factors CLB19, CRR28 and RARE1: retention of CLB19 despite a phylogenetically deep loss of its two known editing targets in Poaceae [J]. BMC Evol Biol, 2018, 18(1): 85. doi:10.1186/s12862-018-1203-4.

    Genome-wide Identification ofGene Family and Expression Analysis of Albino Related Genes in Tea Plants

    LIU Dingding1,2, WANG Junya1,2, TANG Rongjin1,2, CHEN Liang1*, MA Chunlei1*

    1. Tea Research Institute of the Chinese of Agricultural Sciences, National Center for Tea Improvement, Key Laboratory of Tea Biology and Resource Utilization, Ministry of Agriculture, Hangzhou 310008, China; 2. Graduate School of Chinese Academy of Agriculture Science, Beijing 100081, China

    Pentatricopeptide repeat (PPR) proteins are a kind of sequence-specific RNA binding proteins and targeted at semi-autonomous organelles, which play essential roles in play growth and development. In this study, thegenes were systematically identified by bioinformatics analysis based on the tea genome data. Then, their subcellular localization, physicochemical properties, gene structures, the chromosome locations and distribution were analyzed. The results show that a total of 858 putativemembers were obtained from the genome data, which belong to P and PLS subfamilies. Domain analysis shows that each domains were relatively conservative in tea plants. Subcellular localization prediction indicates more than half of CsPPR proteins were located in the chloroplasts. Gene structure analysis shows that 31% ofgenes lacked intron and the gene family had undergone extensive gene duplication events in the process of evolution. Subsequently, In order to investigate the role ofgene family in regulating the gene expressions of albino tea plants, transcriptome analysis was performed on the normal leaf color cultivar ‘Shuchazao’ and five albino tea cultivars such as ‘Anji Huangcha’. And 24 differential co-expressedgenes were identified from transcriptome data of five groups, and the real-time quantitative PCR technology was used to analysis the expression pattern of the 24genes in different cultivars and tissues of tea plants, and the results show that the majority of them were highly expressed in shoots, mature leaves and stems. The research results would provide a basis forgene cloning and functional research.

    tea plants,gene family, transcriptome, bioinformatics, real-time quantitative PCR

    S571.1;Q946.84+1;Q939.1

    A

    1000-369X(2021)02-159-14

    2020-07-13

    2020-09-15

    中國農(nóng)業(yè)科學院科技創(chuàng)新工程(CAAS-ASTIP-2017-TRICAAS)、現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系(CARS-19)、浙江省農(nóng)業(yè)(茶樹)新品種選育重大科技專項子課題(2016C02053-5)、中國農(nóng)業(yè)科學院茶葉研究所基本科研業(yè)務費(1610212017008、1610212019004)

    劉丁丁,女,碩士研究生,主要從事茶樹資源育種與遺傳改良研究。*通信作者:liangchen@tricaas.com;malei220@tricaas.com

    (責任編輯:趙鋒).

    猜你喜歡
    白化葉綠體結構域
    白化黃喉擬水龜人工培育研究①
    最嚴重白化
    蛋白質(zhì)結構域劃分方法及在線服務綜述
    南方紅豆杉葉綠體非編碼序列PCR體系優(yōu)化及引物篩選
    重組綠豆BBI(6-33)結構域的抗腫瘤作用分析
    組蛋白甲基化酶Set2片段調(diào)控SET結構域催化活性的探討
    白化茶種質(zhì)資源分類研究
    茶葉(2015年3期)2015-12-13 06:31:06
    白化和紫化茶種質(zhì)資源開發(fā)進展
    泛素結合結構域與泛素化信號的識別
    茶樹葉綠體DNA的PCR-RFLP反應體系優(yōu)化
    食品科學(2013年6期)2013-03-11 18:20:13
    亚洲成av片中文字幕在线观看| av片东京热男人的天堂| 91麻豆av在线| 老熟妇乱子伦视频在线观看| 欧美日本亚洲视频在线播放| 精品国产乱码久久久久久男人| 午夜免费激情av| av网站在线播放免费| 免费在线观看日本一区| 一边摸一边做爽爽视频免费| 国产精品久久久av美女十八| 嫩草影院精品99| 欧美日韩视频精品一区| 国产在线精品亚洲第一网站| 精品国产超薄肉色丝袜足j| av视频免费观看在线观看| 天堂中文最新版在线下载| 国产欧美日韩综合在线一区二区| 日日爽夜夜爽网站| 亚洲全国av大片| 日本精品一区二区三区蜜桃| 99在线视频只有这里精品首页| 久久伊人香网站| 久久精品亚洲熟妇少妇任你| 搡老岳熟女国产| 国产精品自产拍在线观看55亚洲| 精品久久久久久电影网| 美女福利国产在线| 亚洲自拍偷在线| 免费高清在线观看日韩| 女人被狂操c到高潮| 亚洲av片天天在线观看| 涩涩av久久男人的天堂| 亚洲人成77777在线视频| bbb黄色大片| 啦啦啦在线免费观看视频4| 精品国产超薄肉色丝袜足j| 日韩成人在线观看一区二区三区| 老司机亚洲免费影院| 最新美女视频免费是黄的| 熟女少妇亚洲综合色aaa.| 成人永久免费在线观看视频| 国产亚洲精品综合一区在线观看 | 免费在线观看日本一区| 精品卡一卡二卡四卡免费| 久久精品亚洲精品国产色婷小说| 国产精品乱码一区二三区的特点 | 日韩中文字幕欧美一区二区| 超碰97精品在线观看| 少妇粗大呻吟视频| 无遮挡黄片免费观看| 欧美激情 高清一区二区三区| 国产精品 国内视频| 久久伊人香网站| 亚洲专区国产一区二区| 日韩欧美三级三区| 亚洲人成电影观看| 窝窝影院91人妻| 亚洲成人久久性| 99在线人妻在线中文字幕| 欧美人与性动交α欧美精品济南到| 色婷婷av一区二区三区视频| 91九色精品人成在线观看| 神马国产精品三级电影在线观看 | 欧美激情久久久久久爽电影 | 久久久久久久久久久久大奶| 国产精品一区二区免费欧美| 91字幕亚洲| 十八禁网站免费在线| 国产精品av久久久久免费| 国产xxxxx性猛交| 国产成人欧美| 久久人人97超碰香蕉20202| 在线观看免费视频网站a站| 国产精品久久久久久人妻精品电影| 性欧美人与动物交配| 又黄又粗又硬又大视频| 在线永久观看黄色视频| 中文欧美无线码| 精品一区二区三区av网在线观看| 伊人久久大香线蕉亚洲五| 国产免费现黄频在线看| 日韩有码中文字幕| 精品国产国语对白av| 久久人人精品亚洲av| 精品国产乱码久久久久久男人| 少妇的丰满在线观看| 欧美日韩瑟瑟在线播放| 美女高潮喷水抽搐中文字幕| 色婷婷av一区二区三区视频| 亚洲欧美精品综合一区二区三区| 99在线人妻在线中文字幕| 国产精品美女特级片免费视频播放器 | 男女之事视频高清在线观看| 亚洲精华国产精华精| 国产野战对白在线观看| 999久久久国产精品视频| 在线观看一区二区三区激情| 国产精品免费一区二区三区在线| 天堂中文最新版在线下载| 两性夫妻黄色片| √禁漫天堂资源中文www| 黄片播放在线免费| 窝窝影院91人妻| 曰老女人黄片| 热re99久久精品国产66热6| 免费看a级黄色片| 欧美黄色片欧美黄色片| 国产欧美日韩一区二区精品| 成人永久免费在线观看视频| 80岁老熟妇乱子伦牲交| 欧美乱妇无乱码| 国产99久久九九免费精品| 纯流量卡能插随身wifi吗| 色综合站精品国产| 欧美丝袜亚洲另类 | 日日摸夜夜添夜夜添小说| 亚洲中文日韩欧美视频| 日韩 欧美 亚洲 中文字幕| 午夜福利在线观看吧| 最近最新中文字幕大全电影3 | 中文字幕高清在线视频| 一级a爱片免费观看的视频| 国产精品久久久久久人妻精品电影| 搡老乐熟女国产| 免费在线观看影片大全网站| 男女之事视频高清在线观看| 淫妇啪啪啪对白视频| 日韩免费av在线播放| 国产精品av久久久久免费| 熟女少妇亚洲综合色aaa.| 夜夜看夜夜爽夜夜摸 | 天堂影院成人在线观看| 国产av在哪里看| 久久99一区二区三区| 正在播放国产对白刺激| 午夜免费成人在线视频| 免费搜索国产男女视频| 怎么达到女性高潮| 一a级毛片在线观看| 又大又爽又粗| 久久久久久久久中文| 涩涩av久久男人的天堂| 色综合站精品国产| 成年人黄色毛片网站| 美女 人体艺术 gogo| 久久精品国产99精品国产亚洲性色 | 亚洲精品美女久久av网站| 国产午夜精品久久久久久| 视频区欧美日本亚洲| 国产精品偷伦视频观看了| 国产av一区在线观看免费| 午夜福利影视在线免费观看| 神马国产精品三级电影在线观看 | aaaaa片日本免费| 精品高清国产在线一区| 99久久久亚洲精品蜜臀av| 中文字幕人妻丝袜一区二区| 亚洲五月婷婷丁香| 成人国语在线视频| √禁漫天堂资源中文www| 亚洲国产精品合色在线| 天天躁夜夜躁狠狠躁躁| 男女床上黄色一级片免费看| 欧美成人性av电影在线观看| 欧美日韩福利视频一区二区| 亚洲成国产人片在线观看| 免费在线观看日本一区| 欧美成人午夜精品| 一级a爱片免费观看的视频| 欧美日韩中文字幕国产精品一区二区三区 | 女人被躁到高潮嗷嗷叫费观| 久热这里只有精品99| а√天堂www在线а√下载| aaaaa片日本免费| 国产成人精品久久二区二区免费| 每晚都被弄得嗷嗷叫到高潮| 母亲3免费完整高清在线观看| 超碰成人久久| 乱人伦中国视频| 精品久久久久久久毛片微露脸| 色哟哟哟哟哟哟| 精品一区二区三卡| 琪琪午夜伦伦电影理论片6080| av超薄肉色丝袜交足视频| 国产精品成人在线| 高清av免费在线| 一级毛片女人18水好多| 99国产精品99久久久久| 高清欧美精品videossex| 欧美国产精品va在线观看不卡| 99久久人妻综合| 亚洲精品粉嫩美女一区| 日韩大码丰满熟妇| 老司机福利观看| 国产免费男女视频| 久久久久亚洲av毛片大全| 成年人免费黄色播放视频| 在线av久久热| 亚洲国产欧美日韩在线播放| 交换朋友夫妻互换小说| 日日夜夜操网爽| 美女午夜性视频免费| 国产熟女午夜一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 五月开心婷婷网| 精品久久久久久成人av| 最好的美女福利视频网| 青草久久国产| 国产精品99久久99久久久不卡| 五月开心婷婷网| 日韩欧美在线二视频| 亚洲全国av大片| a级毛片在线看网站| 色综合婷婷激情| 久久精品国产清高在天天线| 欧美一级毛片孕妇| 午夜免费成人在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 精品少妇一区二区三区视频日本电影| 日本a在线网址| 久久青草综合色| 久久国产精品男人的天堂亚洲| 波多野结衣一区麻豆| 欧美激情久久久久久爽电影 | 99国产精品一区二区蜜桃av| 欧美乱码精品一区二区三区| 深夜精品福利| 免费av中文字幕在线| 伊人久久大香线蕉亚洲五| 18美女黄网站色大片免费观看| 精品卡一卡二卡四卡免费| 亚洲精品一二三| 欧美丝袜亚洲另类 | 不卡av一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲人成电影免费在线| 久久 成人 亚洲| 天堂俺去俺来也www色官网| 亚洲av成人av| 少妇的丰满在线观看| 亚洲人成网站在线播放欧美日韩| 黄色视频不卡| 老司机午夜福利在线观看视频| 在线观看66精品国产| 国产免费现黄频在线看| 色尼玛亚洲综合影院| 韩国精品一区二区三区| 免费人成视频x8x8入口观看| 国产精品一区二区免费欧美| 精品久久久久久,| 夜夜爽天天搞| 日日夜夜操网爽| 精品人妻1区二区| 亚洲激情在线av| 美国免费a级毛片| 91字幕亚洲| 国内毛片毛片毛片毛片毛片| 久久这里只有精品19| 久久精品亚洲精品国产色婷小说| 国产精品影院久久| 亚洲九九香蕉| 久久久久国产一级毛片高清牌| 亚洲国产精品合色在线| 在线观看一区二区三区| 中国美女看黄片| 两个人看的免费小视频| 久久婷婷成人综合色麻豆| 国产精品九九99| 国产精品秋霞免费鲁丝片| 99国产精品99久久久久| 久久国产乱子伦精品免费另类| 香蕉国产在线看| 水蜜桃什么品种好| 免费人成视频x8x8入口观看| 变态另类成人亚洲欧美熟女 | 两性夫妻黄色片| 国产精品一区二区精品视频观看| 国产伦人伦偷精品视频| 精品久久久精品久久久| 久久久精品欧美日韩精品| 国产aⅴ精品一区二区三区波| 精品久久久久久,| 午夜精品在线福利| 亚洲欧美激情在线| 在线观看午夜福利视频| 欧美精品亚洲一区二区| 精品福利永久在线观看| 免费日韩欧美在线观看| 国产亚洲精品久久久久5区| 啦啦啦在线免费观看视频4| 宅男免费午夜| 免费不卡黄色视频| 精品福利观看| 99久久99久久久精品蜜桃| www日本在线高清视频| 搡老岳熟女国产| 天堂中文最新版在线下载| 国产单亲对白刺激| 国产乱人伦免费视频| 国产免费av片在线观看野外av| 久久中文字幕人妻熟女| 性欧美人与动物交配| 欧美不卡视频在线免费观看 | 叶爱在线成人免费视频播放| 免费不卡黄色视频| 久久久久国产精品人妻aⅴ院| 午夜免费观看网址| 欧美日韩视频精品一区| 黑人操中国人逼视频| 天堂√8在线中文| 最新美女视频免费是黄的| 一区二区日韩欧美中文字幕| 黄色片一级片一级黄色片| 一边摸一边抽搐一进一出视频| 极品教师在线免费播放| 可以在线观看毛片的网站| 99热只有精品国产| 国产日韩一区二区三区精品不卡| 亚洲,欧美精品.| 夫妻午夜视频| 国产av一区二区精品久久| 精品久久蜜臀av无| 亚洲美女黄片视频| 可以免费在线观看a视频的电影网站| 久久久久亚洲av毛片大全| 精品熟女少妇八av免费久了| 亚洲全国av大片| 一级作爱视频免费观看| 天天躁夜夜躁狠狠躁躁| 久久这里只有精品19| 天天躁夜夜躁狠狠躁躁| 国产精品秋霞免费鲁丝片| 午夜精品国产一区二区电影| 国产成人一区二区三区免费视频网站| 亚洲精品在线观看二区| 国产区一区二久久| 久久国产乱子伦精品免费另类| 亚洲欧美激情综合另类| 色在线成人网| 色综合欧美亚洲国产小说| 欧美日本亚洲视频在线播放| 欧美黄色淫秽网站| 久久国产精品男人的天堂亚洲| 国产成人欧美在线观看| 99re在线观看精品视频| 欧美另类亚洲清纯唯美| 日韩成人在线观看一区二区三区| 大码成人一级视频| 午夜a级毛片| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 亚洲精品久久成人aⅴ小说| 欧美日韩亚洲国产一区二区在线观看| 一a级毛片在线观看| 成人三级做爰电影| 两个人免费观看高清视频| 欧美日韩亚洲国产一区二区在线观看| 男女下面进入的视频免费午夜 | 免费高清在线观看日韩| 精品国产国语对白av| 国产一区二区三区在线臀色熟女 | 男女之事视频高清在线观看| 热99国产精品久久久久久7| 日韩欧美免费精品| 不卡av一区二区三区| 18禁美女被吸乳视频| 日韩免费高清中文字幕av| 91av网站免费观看| 99国产精品一区二区三区| 日韩视频一区二区在线观看| 在线观看舔阴道视频| 极品教师在线免费播放| 免费日韩欧美在线观看| 精品福利永久在线观看| 欧美成人午夜精品| 久久九九热精品免费| 国产一区二区三区视频了| 超色免费av| 精品国产亚洲在线| 在线观看免费日韩欧美大片| 国产精品乱码一区二三区的特点 | 日韩成人在线观看一区二区三区| 国产精品久久久久久人妻精品电影| 丝袜美足系列| 黄色视频不卡| 国内久久婷婷六月综合欲色啪| 亚洲五月色婷婷综合| 国产精品av久久久久免费| 亚洲精品国产色婷婷电影| www.自偷自拍.com| 99国产精品一区二区蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 国产精品影院久久| 18美女黄网站色大片免费观看| 亚洲情色 制服丝袜| 婷婷精品国产亚洲av在线| 免费搜索国产男女视频| 成人18禁高潮啪啪吃奶动态图| 91字幕亚洲| 中亚洲国语对白在线视频| 可以免费在线观看a视频的电影网站| 色婷婷久久久亚洲欧美| 欧美黄色片欧美黄色片| 久久久久国产一级毛片高清牌| 精品一区二区三区视频在线观看免费 | 亚洲精品国产精品久久久不卡| e午夜精品久久久久久久| a级毛片黄视频| 午夜福利免费观看在线| 我的亚洲天堂| 国产高清videossex| 精品高清国产在线一区| 88av欧美| 黄色丝袜av网址大全| 啦啦啦免费观看视频1| 国产熟女午夜一区二区三区| 成人特级黄色片久久久久久久| 俄罗斯特黄特色一大片| 成人影院久久| 国产高清激情床上av| 人妻久久中文字幕网| 日韩欧美国产一区二区入口| 国产精品偷伦视频观看了| 又紧又爽又黄一区二区| 久久香蕉国产精品| 久久久国产成人精品二区 | 露出奶头的视频| a级毛片黄视频| 大陆偷拍与自拍| 50天的宝宝边吃奶边哭怎么回事| 日韩 欧美 亚洲 中文字幕| 午夜免费激情av| 国产乱人伦免费视频| 亚洲人成77777在线视频| 久久热在线av| 波多野结衣av一区二区av| 亚洲人成电影免费在线| 亚洲中文字幕日韩| 51午夜福利影视在线观看| 男女高潮啪啪啪动态图| 久久国产乱子伦精品免费另类| 99香蕉大伊视频| 水蜜桃什么品种好| 亚洲精品中文字幕一二三四区| 亚洲专区字幕在线| 男女下面插进去视频免费观看| 五月开心婷婷网| 怎么达到女性高潮| 一个人免费在线观看的高清视频| 又黄又爽又免费观看的视频| 高清毛片免费观看视频网站 | 十八禁人妻一区二区| 一边摸一边做爽爽视频免费| 在线观看日韩欧美| 99国产极品粉嫩在线观看| 中文字幕色久视频| 1024香蕉在线观看| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲| 亚洲av五月六月丁香网| 青草久久国产| 男女午夜视频在线观看| 精品久久久久久成人av| 长腿黑丝高跟| 精品国产国语对白av| 99精国产麻豆久久婷婷| 国产成人免费无遮挡视频| 亚洲国产看品久久| 成人三级黄色视频| 波多野结衣一区麻豆| 欧美日韩av久久| av有码第一页| 在线观看免费高清a一片| 亚洲成av片中文字幕在线观看| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣高清无吗| 亚洲五月婷婷丁香| 欧美日韩精品网址| 国产成人精品久久二区二区91| 在线观看舔阴道视频| 在线天堂中文资源库| 黄片播放在线免费| 久久精品国产清高在天天线| 夫妻午夜视频| 在线观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产男靠女视频免费网站| x7x7x7水蜜桃| 亚洲成人久久性| 91国产中文字幕| 亚洲精品在线美女| 99riav亚洲国产免费| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av熟女| 国产免费男女视频| 性色av乱码一区二区三区2| 国产精品久久久久成人av| 一级片免费观看大全| 久久久久久人人人人人| 久久久国产成人精品二区 | 亚洲精品国产一区二区精华液| 日韩欧美三级三区| 99久久99久久久精品蜜桃| 日韩人妻精品一区2区三区| 亚洲一区二区三区色噜噜 | 91字幕亚洲| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲第一青青草原| 午夜福利影视在线免费观看| 最近最新免费中文字幕在线| 嫩草影视91久久| 久久久国产成人精品二区 | 9191精品国产免费久久| 久久久久久久精品吃奶| 午夜精品国产一区二区电影| 人妻久久中文字幕网| 长腿黑丝高跟| 免费在线观看视频国产中文字幕亚洲| www.www免费av| 午夜免费成人在线视频| www.熟女人妻精品国产| 国产免费现黄频在线看| 国产精品99久久99久久久不卡| 操出白浆在线播放| av在线播放免费不卡| 高潮久久久久久久久久久不卡| 日本免费a在线| 亚洲欧美日韩另类电影网站| 18禁黄网站禁片午夜丰满| 久久狼人影院| 久久国产亚洲av麻豆专区| 亚洲一码二码三码区别大吗| 三上悠亚av全集在线观看| 国产亚洲欧美精品永久| 亚洲久久久国产精品| 欧美性长视频在线观看| 如日韩欧美国产精品一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品国产色婷婷电影| 久久久精品欧美日韩精品| 国产视频一区二区在线看| 午夜老司机福利片| 亚洲视频免费观看视频| 久久久久久久久中文| 久久精品成人免费网站| av中文乱码字幕在线| 正在播放国产对白刺激| 亚洲九九香蕉| 在线观看www视频免费| 成年女人毛片免费观看观看9| 十八禁网站免费在线| 两性午夜刺激爽爽歪歪视频在线观看 | 91九色精品人成在线观看| 又大又爽又粗| 日韩成人在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 日本免费a在线| 亚洲国产欧美日韩在线播放| 人妻久久中文字幕网| 久久久久久久午夜电影 | 免费在线观看黄色视频的| 午夜a级毛片| 久久 成人 亚洲| 中国美女看黄片| 午夜影院日韩av| xxx96com| 大型av网站在线播放| 欧美一区二区精品小视频在线| √禁漫天堂资源中文www| 亚洲人成77777在线视频| 国产精品亚洲一级av第二区| 久久精品国产99精品国产亚洲性色 | 国产一区在线观看成人免费| 免费在线观看日本一区| 狠狠狠狠99中文字幕| 18禁观看日本| 成人影院久久| 91av网站免费观看| 少妇粗大呻吟视频| 欧洲精品卡2卡3卡4卡5卡区| 成在线人永久免费视频| 成人特级黄色片久久久久久久| 女人被狂操c到高潮| 国产黄a三级三级三级人| 午夜激情av网站| 两人在一起打扑克的视频| 十八禁网站免费在线| 美女国产高潮福利片在线看| 无人区码免费观看不卡| 一进一出好大好爽视频| 老司机深夜福利视频在线观看| 欧美日韩福利视频一区二区| 两个人免费观看高清视频| 午夜福利影视在线免费观看| 青草久久国产| 在线观看免费视频网站a站| 欧美日韩国产mv在线观看视频| 久久性视频一级片| 久久狼人影院| 欧美乱码精品一区二区三区| 欧美成狂野欧美在线观看| 久久狼人影院| 在线国产一区二区在线| 欧美成狂野欧美在线观看| 一区二区日韩欧美中文字幕| 国产亚洲精品综合一区在线观看 | 波多野结衣av一区二区av| 国产99久久九九免费精品| av超薄肉色丝袜交足视频| 波多野结衣av一区二区av| 久久国产亚洲av麻豆专区| 国产成人影院久久av| 日韩 欧美 亚洲 中文字幕|