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

    基于轉錄組測序的烏桕葉色相關基因的挖掘

    2024-09-24 00:00:00王曉曉鄭旭趙文靜繆美華劉興滿
    江蘇農業(yè)學報 2024年8期

    收稿日期:2023-10-11

    基金項目:中央財政林業(yè)科技推廣示范基金項目[蘇(2021)TG12];江蘇省蘇北科技專項(SZ-LYG202142);連云港市財政局專項基金項目(QNJJ2205)

    作者簡介:王曉曉(1996-),女,山東濟南人,碩士,研究實習員,主要從事觀賞樹種栽培及遺傳育種研究。(E-mail)2402287168@qq.com

    通訊作者:劉興滿,(E-mail)yyyy28@126.com

    摘要: 烏桕(Triadica sebifera)耐鹽堿能力強,葉色隨季節(jié)變化顯著,尤以秋季觀賞性最佳,是園林綠化中常用的彩葉植物。為揭示烏桕葉片呈色的分子機制,本研究對自主選育的3個烏桕品種連桕1號、云臺紅桕、云臺金桕轉色前后的葉片進行轉錄組測序,共獲得123.78 GB測序數據,70 815個高質量單一基因序列(Unigenes),其中48 367個Unigenes得到了功能注釋。將3個烏桕品種轉色前后的轉錄組數據進行兩兩對比,各比較組之間鑒定出共有的差異表達基因(DEG)2 674個。GO分析發(fā)現,DEG主要與細胞過程、細胞、結合等生物學功能相關。通過KEGG Pathway功能分析,篩選出53個與烏桕葉色相關的DEG,包括與花青素合成有關的基因14個,與卟啉和葉綠素代謝有關的基因22個,與類胡蘿卜素生物合成有關的基因17個。其中花青素合成途徑相關基因在烏桕轉色后相對表達量上升,葉綠素、類胡蘿卜素合成途徑相關基因相對表達量下降,而葉綠素、類胡蘿卜素降解途徑相關基因相對表達量上升,該結果與葉片中色素含量的變化趨勢一致。本研究結果可為探究烏桕葉片呈色機制提供理論依據。

    關鍵詞: 烏桕;葉色變化;轉錄組;基因挖掘;分子機制;呈色

    中圖分類號: S687.9 文獻標識碼: A 文章編號: 1000-4440(2024)08-1521-12

    Mining of genes related to Triadica sebifera leaf color based on transcriptome sequencing

    WANG Xiaoxiao, ZHENG Xu, ZHAO Wenjing, MIAO Meihua, LIU Xingman

    (Lianyungang Academy of Agricultural Sciences, Lianyungang 222000, China)

    Abstract: Triadica sebifera is a frequently-used color-leafed plant in landscaping because of its strong salt-tolerance and remarkable leaf color change with seasons, and it has the best ornamental properties in autumn. In order to reveal the molecular mechanism of leaf coloration of Triadica sebifera, the transcriptomes of leaves of three self-selected varieties of Triadica sebifera (Lianjiu No.1, Yuntaihongjiu and Yuntaijinjiu) were sequenced before and after leaf color changed. A total of 123.78 GB sequencing data were obtained, with 70 815 high-quality Unigenes, of which 48 367 were functionally annotated. A total of 2 674 shared differentially expressed genes (DEGs) of different comparison groups were identified from the transcriptome data of three Triadica sebifera varieties before and after color transformation, through pairwise comparison. Go analysis revealed that DEGs were mainly related to biological functions such as Cell process, Cell and Binding. Through functional analysis of KEGG Pathway, 53 DEGs related to leaf color were screened, including 14 genes related to anthocyanin synthesis, 22 genes related to porphyrin and chlorophyll metabolism and 17 genes related to biosynthesis of carotenoid. The relative expression levels of genes related to anthocyanin biosynthesis pathway were increased and those related to chlorophyll and carotenoid biosynthesis pathways were decreased, the relative expression levels of genes related to chlorophyll and carotenoid degradation pathways were increased, which results were consistent with the change trends of pigment content in leaves. The results of this study can provide a theoretical basis for exploring the mechanism of leaf coloration of Triadica sebifera.

    Key words: Triadica sebifera;leaf color change;transcriptome;gene mining;molecular mechanism;coloration

    烏桕(Triadica sebifera)隸屬大戟科(Euphorbiaceae)烏桕屬(Triadica),在中國栽培歷史悠久,耐鹽堿能力強,具有觀賞、油用、藥用等多種用途。烏桕樹體高大,在秋季呈現綠、紅、黃、紫等豐富的葉色,因此常用于道路綠化、公園造景,極具觀賞價值[1]。目前對烏桕的研究主要集中在繁殖栽培[2]、油用品質[3]、種子性狀[4]等方面。倪正[5]通過對比金黃4號、合紅6號和滁紫1號不同發(fā)育期葉片色素含量和N元素含量,闡明了兩者對葉色的影響。吳飛洋等[6]通過探究影響烏桕葉片轉色的環(huán)境因素,發(fā)現山地紅壤條件下的全光照處理更有利于秋季烏桕葉色的呈現。目前對烏桕葉片呈色機理的研究主要集中在生理特性和環(huán)境因子上,在分子水平上對烏桕呈色機制的研究還鮮有報道。

    色素是植物葉片呈色的物質基礎,葉綠素、類胡蘿卜素、花色素苷是決定葉色的主要色素,其中由花青素和糖組成的糖苷是紅葉植物的主要呈色物質[7]。當色素含量和比例發(fā)生變化時,葉色也會隨之改變[8]。近年來,轉錄組測序(RNA-Seq)技術不斷迭代升級,已成為研究彩葉植物呈色機制的重要工具。利用轉錄組測序技術,Dong等[9]發(fā)現,在熱脅迫下雞爪槭(Acer palmatum)中花青素和類胡蘿卜的生物合成相關基因表達量顯著上調,而葉綠素的生物合成相關基因表達量顯著下調,色素合成相關基因表達量的變化是葉色由紅變綠的直接分子機制。Liu等[10]對三色觀賞甘藍(Brassica oleracea L. var. acephala)轉錄組進行測序發(fā)現,甘藍葉片從綠色轉為白色是由于葉綠素生物合成過程受到抑制。

    連桕1號、云臺紅桕及云臺金桕均是從烏桕實生苗中選育出的新品種。其中,連桕1號葉片碩大,葉面積為普通烏桕的2倍,秋季葉片呈紫紅色,最佳觀葉期25 d。云臺紅桕早秋葉片呈亮紅色,最佳觀葉期比連桕1號提早1周。秋季云臺金桕葉片呈金黃色,葉片轉色早,觀賞期長。此外,3個烏桕品種耐鹽能力強,均可在含鹽量5‰的沿海地區(qū)種植。本研究以上述品種轉色前后的葉片(綠色葉、紅色葉、黃色葉)為研究對象,探究色素合成相關基因的差異表達情況,以揭示烏桕葉片呈色機制,同時為調控烏桕葉色、延長觀賞期提供理論依據。

    1 材料與方法

    1.1 試驗材料

    試驗材料為連云港市農業(yè)科學院東辛試驗基地自主選育的3個烏桕品種:連桕1號、云臺紅桕、云臺金桕。于2021年9月至11月分別采集3個烏桕品種的轉色前和轉色后的葉片,3個烏桕品種的葉片轉色前均為綠色,轉色后分別為紫紅色、紅色和金黃色,每種葉片設置3次生物學重復,依次進行編號(表1)。

    1.2 試驗方法

    1.2.1 色素含量的測定 花青素、葉綠素和類胡蘿卜素含量分別參照植物原花青素試劑盒、植物葉綠素含量試劑盒、類胡蘿卜素含量試劑盒(蘇州科銘生物技術有限公司產品)說明書,使用分光光度法進行測定,均為鮮重含量。

    1.2.2 總RNA的提取與質量檢測 使用MiniBEST Plant RNA Extraction試劑盒(日本TaKaRa公司產品)提取總RNA,隨后使用核酸蛋白質濃度測定儀(德國Eppendorf公司產品,產品型號為BioPhotometer Plus 6132)檢測RNA的純度、濃度和完整度,將RNA保存于-80 ℃冰箱中備用。

    1.2.3 cDNA文庫的構建及測序 用帶有Oligo dT的磁珠特異性結合帶有PolyA尾的mRNA,利用打斷試劑使mRNA片段化,以打斷后的mRNA為模板合成一鏈cDNA,隨后利用二鏈合成反應體系合成二鏈cDNA。使用試劑盒純化回收、末端修復并連接測序接頭,然后對片段大小進行選擇,最后進行PCR擴增。cDNA文庫的測序在Illumina平臺上進行。

    1.2.4 測序數據的過濾與組裝 對原始數據(Raw reads)進行過濾,得到高質量序列(Clean reads)。使用Trinity軟件進行de novo組裝,利用CD-HIT軟件聚類得到非冗余序列(Unigenes)。使用eXpress軟件計算Unigenes的表達水平,表達水平用FPKM值表示。

    1.2.5 Unigenes功能注釋 用7個功能數據庫對組裝得到的Unigenes進行注釋。其中KEGG、GO、NR、eggNOG、SwissProt、KOG注釋通過Diamond軟件進行,Pfam注釋通過HMMER軟件進行。

    1.2.6 差異表達基因(DEG)的篩選 首先篩選保留counts>2的基因,利用DESeq2軟件對各樣本基因的counts數目進行標準化處理,計算差異表達倍數(Fold change),將q<0.05且Fold change>2作為篩選條件,采用負二項分布檢驗的方式進行差異顯著性檢驗,最終根據檢驗結果確定DEG。

    1.2.7 DEGs的表達模式分析 使用TBtools軟件可視化展示候選DEG的相對表達量。

    2 結果與分析

    2.1 烏桕轉色前后色素含量的變化

    由圖1可以看出,轉色后連桕1號、云臺紅桕、云臺金桕葉片的花青素含量顯著高于轉色前,而葉綠素a、葉綠素b和類胡蘿卜素含量均顯著低于轉色前(P<0.05)。其中,連桕1號轉色后的花青素含量是轉色前的9.84倍,葉綠素(a+b)和類胡蘿卜素含量分別為轉色前的64.30%和51.70%;云臺紅桕轉色后的花青素含量是轉色前的13.20倍,葉綠素(a+b)和類胡蘿卜素含量分別為轉色前的23.94%和42.62%;云臺金桕轉色后的花青素含量是轉色前的5.54倍,葉綠素(a+b)和類胡蘿卜素含量分別為轉色前的4.56%和49.18%。

    2.2 轉錄組測序結果及質量評估

    以3個烏桕品種轉色前后的葉片為材料,構建18個cDNA文庫,使用Illumina平臺進行轉錄組測序,共獲得123.78 GB序列數據,拼接出70 815個Unigenes。將Raw reads過濾后,平均每個文庫的Clean reads為47.48 MB,有效數據量為6.69~7.14 GB,Q30堿基均不低于94.75%,平均GC含量為42.74%(表2)。以上研究結果表明,本次測序數據質量良好,可用于后續(xù)分析。

    2.3 Unigenes功能注釋

    將得到的Unigenes在7個數據庫(NR、SwissProt、KEGG、KOG、eggNOG、GO和Pfam)中進行功能注釋的結果見表3。有6 633個Unigenes(9.37%)被7個數據庫同時注釋,有48 367個Unigenes(68.30%)被任一數據庫注釋。注釋到Unigenes數量最多的為NR數據庫,共有48 103個Unigenes與其顯著相關,占總數的67.93%,其中,序列長度為300~<1 000 bp的Unigenes占24.53%,序列長度≥1 000 bp的占43.40%。

    NR注釋結果顯示,與烏桕具有同源性的物種有10個,其中巴西橡膠樹同源性較高,達28.72%,其次為蓖麻、麻風樹、木薯,同源性均在12.54%及以上(圖2)。

    隨后將烏桕Unigenes序列注釋到KOG數據庫,得到基因同源物的分類信息,結果如圖3所示。32 164個Unigenes根據基因功能歸納為25個組,其中注釋到功能預測(General function prediction only)類別的Unigenes數量最多,有6 635個,占注釋總數的20.63%;其次是翻譯后修飾(Posttranslational modification)、蛋白質轉換(Protein turnover)、分子伴侶(Chaperones)類別,有3 295個Unigenes,占總數的10.24%,注釋到信號轉導機制(Signal transduction mechanisms)類別的Unigenes有2 799個,占總數的8.70%。

    2.4 Unigenes的差異表達分析

    為獲得烏桕葉色變化的相關基因,將連桕1號轉色前(LGA)、連桕1號轉色后(LGB)、云臺紅桕轉色前(LGC)、云臺紅桕轉色后(LGD)、云臺金桕轉色前(LGE)、云臺金桕轉色后(LGF)的轉錄組數據進行兩兩對比,分為LGB和LGA、LGD和LGC、LGF和LGE、LGD和LGB和LGF和LGB 5個比較組,并將q<0.05且Fold change>2作為篩選條件來篩選DEG,共鑒定出52 491個DEG,DEG數量統(tǒng)計結果如圖4所示。LGB和LGA的比較組共有24 616個DEG,其中有11 765個DEG表達上調,12 851個表達下調;LGD和LGC的比較組共有28 686個DEG,其中有13 316個DEG表達上調,15 370個DEG表達下調;LGF和LGE的比較組共有37 669個DEG,其中有19 789個DEG表達上調,17 880個DEG表達下調;LGD和LGB的比較組共有11 648個DEG,其中有5 469個DEG表達上調,6 179個DEG表達下調;LGF和LGB的比較組共有24 492個DEG,其中有13 515個DEG表達上調,10 977個DEG表達下調。5個比較組共有的DEG有2 674個(圖4、圖5)。

    2.5 DEG的GO富集分析

    對3個烏桕品種葉片轉色前后的DEG進行GO富集分析發(fā)現,所有的DEG被富集在生物過程、細胞組分和分子功能三大類別。在生物過程類別中,DEG主要富集在細胞過程和代謝過程中,LGB 和LGA的比較組中的DEG數量分別為7 936個和6 746個,LGD和LGC的比較組中DEG數量分別為9 582個和8 224個,LGF和LGE的比較組中DEG數量分別為12 836個和10 923個,LGD和LGB的比較組中DEG數量分別為3 938個和3 450個,LGF和LGB的比較組中DEG數量分別為8 423個和7 296個。在細胞組分類別中,DEG主要富集在細胞和細胞成分中,LGB和LGA的比較組中的DEG數量分別為10 101個和10 090個,LGD和LGC的比較組中DEG數量分別為12 215個和12 197個,LGF和LGE的比較組中DEG數量分別為16 216個和16 191個,LGD和LGB的比較組中DEG數量分別為5 058個和5 051個,LGF和LGB的比較組中DEG數量分別為10 631個和10 613個;此外,在5個比較組中有大量DEG富集在分子功能類別中的結合和催化活性中(圖6)。

    2.6 DEG的KEGG Pathway分析

    KEGG Pathway注釋分類對明確DEG的生物學功能起到重要作用[11]。LGB和LGA、LGD和LGC、LGF和LGE、LGD和LGB、LGF和LGB 5個比較組分別得到了126條、127條、127條、125條、126條富集代謝通路,富集的DEG分別為6 515個、7 765個、10 825個、3 378個、7 428個。5個比較組中DEG富集程度排名前20條的代謝通路如圖7所示,可以看出,DEG在類黃酮生物合成、苯丙素生物合成、卟啉和葉綠素代謝、類胡蘿卜素生物合成等通路中的富集程度較高,可能與烏桕葉色變化過程相關。

    2.7 3種烏桕葉色相關基因的表達模式分析

    基于KEGG富集結果,對花青素生物合成途徑、卟啉和葉綠素代謝途徑以及類胡蘿卜素生物合成途徑相關的53個DEG的相對表達量進行可視化展示。

    本研究在6種不同葉色狀態(tài)的烏桕葉片中,篩選出14個與花青素合成相關的DEG,包括查爾酮合成酶基因(CHS)、黃烷酮-3-羥化酶基因(F3H)、二氫黃酮醇還原酶基因(DFR)、花青素合成酶基因(ANS)(圖8A)。隨著連桕1號、云臺紅桕和云臺金桕葉片在秋季轉色,這些基因的表達水平均上調,且在連桕1號轉色后查爾酮合成酶基因(CHS)、黃烷酮-3-羥化酶基因(F3H)、二氫黃酮醇還原酶基因(DFR)、花青素合成酶基因(ANS)表達水平的提升最為明顯。在卟啉和葉綠素代謝途徑中篩選出22個DEG,隨著葉片發(fā)育,葉色轉紅、轉黃,鎂原卟啉Ⅸ甲基轉移酶基因(CHLM)、葉綠素合成酶基因(CHLG)、二乙烯還原酶基因(DVR)、葉綠素羥甲基還原酶基因(HCAR)、δ-氨基酮戊酸脫水酶(ALAD)、尿卟啉原Ⅲ合成酶基因(UROS)、尿卟啉原Ⅲ脫羧酶基因(UROD)、原葉綠素酸酯氧化還原酶基因(POR)、原卟啉原氧化酶基因(PPOX)等與葉綠素合成相關的基因相對表達量下降;而葉綠素b還原酶基因(NOL)、脫鎂葉綠素甲酯酸a加氧酶基因(PAO)等與葉綠素降解相關的基因相對表達量上升(圖8B)。從類胡蘿卜素生物合成途徑中篩選出17個DEG。與轉色前相比,烏桕轉色后葉片中參與類胡蘿卜素生物合成的多種關鍵酶基因相對表達量下降,如八氫番茄紅素脫氫酶基因(PDS)、類胡蘿卜素異構酶基因(CRTISO)、番茄紅素環(huán)化酶基因(LCYE、LCYB)等,參與類胡蘿卜素降解途徑的部分9-順式環(huán)氧化類胡蘿卜素雙加氧酶(NCED)基因表達量上升,且各個基因之間的表達量相差較大,CYP707A基因家族相關基因的表達量上升(圖8C)。

    3 討論

    Illumina Solexa測序技術是利用邊合成邊測序和可逆終止子原理的雙端測序,具有較高的準確性及靈敏度,是研究植物生長發(fā)育的有力工具[12]。目前,人們已對許多彩葉植物進行了轉錄組測序,如文心蘭(Oncidium hybridum)[13]、羅漢松(Podocarpus costalis)[14]、杜梨(Pyrus betulifolia Bunge)[15]、紫薇(Lagerstroemia indica)[16]等,為研究彩葉植物呈色機理奠定了基礎。本研究對連桕1號、云臺紅桕、云臺金桕轉色前后的葉片進行轉錄組測序,共構建18個cDNA文庫,獲得123.78 GB數據,Q30堿基不低于94.75%,平均GC含量42.74%,表明本次測序質量較高,烏桕轉錄組數據滿足后續(xù)分析要求。

    本研究在NR、SwissProt、KEGG、KOG、eggNOG、GO和Pfam 7個數據庫中共成功注釋到48 367個Unigenes。隨后通過差異表達分析共鑒定到124 437個DEG,其中5個比較組共有的DEG為2 674個。GO功能分類結果表明,烏桕葉色變化與細胞過程和代謝過程(生物過程)、細胞和細胞成分(細胞組分)、結合和催化活性(分子功能)等過程密切相關;KEGG富集分析結果顯示,在類黃酮生物合成、苯丙素的生物合成、卟啉和葉綠素代謝、類胡蘿卜素生物合成等通路中的DEG顯著富集,這些代謝途徑可能與烏桕的葉色變化有關。

    花青素的生物合成受多種結構酶基因的控制,這些基因可分為早期生物合成基因(CHS、F3H)和晚期生物合成基因(DFR、ANS)[17-18]。其中CHS是啟動花青素合成的端口基因,負責調控花青素等類黃酮物質的合成和積累[19]。F3H位于CHS下游,是控制二氫黃酮醇合成的關鍵基因,二氫黃酮醇在DFR的作用下,進一步生成無色的花青素,然后由ANS催化成有顏色的花青素,從而影響植物組織著色[20-22]。本研究中,相較于綠色葉,紅色葉和黃色葉中CHS、F3H、DFR、ANS等基因相對表達量和花青素含量顯著上升,這與前人關于雞爪槭(Acer palmatum)[23]和彩葉桂(Osmanthus fragrans)[24]葉色的研究結果類似。由此推測,這些基因表達上升可能導致花青素的合成與積累量升高,從而促使烏桕紅色和黃色葉片的形成。

    葉綠素是植物光合色素的重要組成之一[25-32],呈綠色或藍綠色。Zhou等[33]分別對綠色和白色花椰菜進行研究發(fā)現,綠色植株中參與葉綠素生物合成的相關基因表達水平更高。POR和CHLG均為葉綠素合成途徑中的關鍵酶,前者能夠催化原葉綠素酸酯向葉綠素的轉化,其亞型PORB還具有平衡植物體內原葉綠素酸酯的作用;后者負責脫植基葉綠素酯化,完成葉綠素生物合成的最終步驟[34-35]。本研究發(fā)現烏桕轉色后POR、CHLG基因相對表達量下降,且在云臺金桕中最為明顯,這可能導致葉綠素生物合成過程受阻,從而使烏桕葉片褪綠變黃。已有研究結果表明,在紫花槭中NOL基因表達量上升是造成葉綠素降解、葉片呈現紅色的重要原因[36],本研究結果與之相似。同時,本研究測定的烏桕轉色前后的葉綠素含量,其變化趨勢與轉錄組分析結果一致。

    類胡蘿卜素的生物合成與葉綠素的生物合成相互協調,其含量會對葉色造成影響[37]。八氫番茄紅素是植物體內率先合成的類胡蘿卜素,經PDS和CRTISO等催化生成番茄紅素,然后在LCYB、LCYE的作用下轉化為α-胡蘿卜素和β-胡蘿卜素[38]。Chai等[39]研究水稻斑馬葉突變體時發(fā)現,編碼CRTISO的ZEBRA2基因突變,會導致保護性類胡蘿卜素含量降低。銀杏葉片轉色后期,LCYB顯著下調表達,同時伴隨著類胡蘿卜素含量下降[40]。9-順式環(huán)氧化類胡蘿卜素雙加氧酶(NCED)是參與類胡蘿卜素降解的重要酶之一,也是控制類胡蘿卜素向脫落酸(ABA)轉化的限速中樞[41]。CYP707A基因家族編碼的脫落酸8’-羥基化酶對ABA含量變化起重要作用。本研究發(fā)現,烏桕在類胡蘿卜素合成階段,PDS、CRTISO、LCYB、LCYE等基因相對表達量下降,而在類胡蘿卜素降解階段NCED、CYP707A等基因相對表達量上升,且類胡蘿卜素含量顯著下降,這可能會導致烏桕葉片中ABA含量增多,從而加速葉片衰老,降低活性氧(ROS)的清除效率,因此光氧化損傷可能是造成突變體葉色變化的原因之一。

    4 結論

    本研究對3個烏桕品種連桕1號、云臺紅桕、云臺金桕轉色前后的葉片進行轉錄組測序及分析,共獲得123.78 GB測序數據,70 815個高質量Unigenes;將全部Unigenes與7個數據庫進行比對,有48 367個Unigenes得到功能注釋。GO分析結果顯示,DEG主要與細胞過程和代謝過程(生物過程)、細胞和細胞成分(細胞組分)、結合和催化活性(分子功能)等生物學功能相關。通過KEGG分析,篩選出與花青素、葉綠素、類胡蘿卜素等相關的53個DEG,其中花青素合成途徑相關基因在烏桕轉色后相對表達量上升,葉綠素、類胡蘿卜素合成途徑相關基因相對表達量下降,而葉綠素、類胡蘿卜素降解途徑相關基因相對表達量上升。同時對3個烏桕品種轉色前后的主要色素含量進行測定發(fā)現,花青素含量在轉色后顯著上升,葉綠素和類胡蘿卜素含量在轉色后顯著下降。

    參考文獻:

    [1] 陳百強,曹受金,沈 鑫,等. 觀賞烏桕品種表型多樣性分析及綜合評價[J]. 分子植物育種,20(12):4108-4120.

    [2] 張廣濤,包厚天,黃衛(wèi)和,等. 不同施氮處理對烏桕容器苗養(yǎng)分狀況的影響[J]. 西北林學院學報,2023,38(1):80-87.

    [3] 王文靜,董泰瑋. 超聲波輔助提取烏桕籽油的工藝優(yōu)化研究[J]. 中國油脂,2019,44(4):6-9.

    [4] 向珊珊,杜洋文,程軍勇,等. 長江流域烏桕優(yōu)良種質資源種子品質性狀評價[J]. 經濟林研究,2023,41(2):30-37.

    [5] 倪 正. 觀葉烏桕優(yōu)株選擇及其葉片呈色生理生化分析[D]. 合肥:安徽農業(yè)大學,2023.

    [6] 吳飛洋,柳新紅,董峰平,等. 光照和土壤對烏桕秋季葉片色素及可溶性糖的影響[J]. 西南林業(yè)大學學報(自然科學),2019,39(6):41-48.

    [7] 張 敏,錢 猛,倪競德,等. 烏桕秋葉轉色前后生理特性及超微結構的變化[J]. 東北林業(yè)大學學報,2012,40(1):20-24.

    [8] 沈周高,謝應華,范曉娟,等. 葉色變異茶樹生理生化特征及其變異機制的研究進展[J]. 安徽農業(yè)大學學報,2023,50(4):582-594.

    [9] DONG H Y, XU B, JI K S. Comparative transcriptome analysis of genes involved in response to thermal stress and leaf colour change of Acer palmatum[J]. Scientia Horticulturae,2019,255:77-85.

    [10]LIU X P, ZHANG B, WU J, et al. Pigment variation and transcriptional response of the pigment synthesis pathway in the S2309 triple-color ornamental kale (Brassica oleracea L. var. acephala) line[J]. Genomics,2020,112(3):2658-2665.

    [11]王曉曉. 鳳丹對高溫脅迫的響應及緩解措施研究[D]. 揚州:揚州大學,2023.

    [12]LIANG Y, CHEN S Y, LIU G S. Application of next generation sequencing techniques in plant transcriptome[J]. Hereditas,2011,33(12):1317-1326.

    [13]SAKURABA Y, PARK S Y, PACK N C. The divergent roles of STAYGREEN (SGR) homologs in chlorophyll degradation[J]. Molecules and Cells,2015,38(5):390-395.

    [14]范 美. 海島羅漢松嫩葉葉色變化過程及分子機制研究[D]. 杭州:浙江農林大學,2023.

    [15]陳 燕,郭 聰,王 瑩,等. 杜梨葉片轉錄組分析[J]. 江蘇農業(yè)科學,2022,50(3):58-67.

    [16]王淑安,王 鵬,張振宇,等. 紫薇葉片轉錄組分析及葉綠素和類胡蘿卜素等途徑基因的鑒定[J]. 分子植物育種,2016,14(5):1059-1066.

    [17]LIU Y, TIKUNOV Y, SCHOUTEN R E, et al. Anthocyanin biosynthesis and degradation mechanisms in solanaceous vegetables:a review[J]. Frontiers in Chemistry,2018,6:52.

    [18]劉淑華,臧丹丹,孫 燕,等. 花青素生物合成途徑及關鍵酶研究進展[J]. 土壤與作物,2022,11(3):336-346.

    [19]LIM W S, and LI J R. Co-expression of onion chalcone isomerase in Del/Ros1-expressing tomato enhances anthocyanin and flavonol production[J]. Plant Cell Tissue Organ Culture,2017,128(1):113-124.

    [20]BUER C S, MUDAY G K, DJORDJEVIC M A. Flavonoids are differentially taken up and transported long distances in Arabidopsis[J]. Plant Physiology,2007,145:478-490.

    [21]李 瑩,高振蕊,張 馳,等. 花青素合成途徑中分子調控機制的研究進展[J]. 生態(tài)學雜志,2015,34(10):2937-2942.

    [22]杜婷婷,宋治華,董碧瑩,等. 木豆類黃酮代謝通路關鍵基因家族的鑒定與表達分析[J]. 農業(yè)生物技術學報,2021,29(12):2289-2303.

    [23]ZHU L, WEN J, MA Q Y, et al. Transcriptome profiling provides insights into leaf color changes in two Acer palmatum genotypes[J]. Bmc Plant Biology,2022,22(1):589.

    [24]黃均華,徐蘊晨,黃子洋,等. 彩葉桂(Osmanthus fragrans)品種‘虔南桂妃’葉片紅色消退過程中的轉錄組分析[J]. 基因組學與應用生物學,2022,41(增刊1):1964-1972.

    [25]李 碩,李美玄,劉英卉,等. 干旱-復水處理對吉林文冠果幼苗生理特性和光合參數的影響[J]. 江蘇農業(yè)科學,2023,51(18):134-143.

    [26]儂文蓮,李秀芳,吳世清,等. 咸水灌溉下生物炭對玉米光合特性及產量的影響[J]. 排灌機械工程學報,2023,41(6):605-613.

    [27]楊 洪,李旭毅,卿發(fā)紅,等. 不同產量水平水稻群體光合特性和產量構成差異[J]. 江蘇農業(yè)學報,2023,39(5):1089-1096.

    [28]李金威,原慧芳,田耀華,等. 5個珠芽魔芋品種(系)的生長發(fā)育及其光合特性分析[J]. 南方農業(yè)學報,2023,54(4):1166-1174.

    [29]王克曉,周 蕊,黃 祥,等. 油菜角果高光譜成熟指數構建[J]. 江蘇農業(yè)學報,2023, 39(3):716-723.

    [30]馬 玲,杜明華,孟 露,等. 基于高光譜成像技術的番茄葉片葉綠素含量檢測[J]. 江蘇農業(yè)科學,2023,51(11):167-174.

    [31]張雪芹,歐陽海波,謝志南,等. 不同品種油梨果實品質、光合特性比較及其相關分析[J]. 南方農業(yè)學報,2023,54(2):405-413.

    [32]王聲鋒,徐亞澤,彭致功,等. 基于PROSAIL模型的夏玉米葉綠素含量垂直分層模擬[J]. 排灌機械工程學報,2022,40(6):603-610.

    [33]ZHOU X J, FEI Z J, THANNHAUSER T W, et al. Transcriptome analysis of ectopic chloroplast development in green curd cauliflower (Brassica oleracea L. var. botrytis)[J]. BMC Plant Biology,2011,11 (1):169.

    [34]REINBOTHE C, BAKKOURI M E, BUHR F, et al. Chlorophyll biosynthesis:spotlight on protochlorophyllide reduction[J]. Trends Plant Science,2010,15:614-624.

    [35]田密霞,周福慧,姜愛麗,等. 蕓薹屬植物呈色機理研究進展[J]. 園藝學報,2023,50(9):1971-1986.

    [36]翁 卓,高玉福,李佳霖,等. 基于轉錄組的紫花槭光合色素相關基因篩選與分析[J]. 分子植物育種,2024,22(6):1881-1891.

    [37]郭 婷. 欒樹彩葉突變體葉色變異機理的初步探究[D]. 北京:北京林業(yè)大學,2022.

    [38]王紫璇,李佳佳,于旭東,等. 高等植物類胡蘿卜素生物合成研究進展[J]. 分子植物育種,2021,19(8):2627-2637.

    [39]CHAI C, FANG J, LIU Y, et al. ZEBRA2, encoding a carotenoid isomerase, is involved in photoprotection in rice[J]. Plant Molecular Biology,2011,75(3):211-221.

    [40]何智沖. 銀杏葉片衰老轉色進程及調控研究[D]. 揚州:揚州大學,2020.

    [41]JI K, KAI W B, ZHAO B, et al. SlNCED1 and SlCYP707A2:key genes involved in ABA metabolism during tomato fruit ripening[J]. Journal of Experimental Botany,2014,65:5243-5255.

    (責任編輯:成紓寒)

    色综合站精品国产| 色av中文字幕| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| 国产精品av视频在线免费观看| 99在线人妻在线中文字幕| 99国产精品一区二区蜜桃av| 97超级碰碰碰精品色视频在线观看| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 午夜免费观看网址| 男人的好看免费观看在线视频| 亚洲欧美精品综合久久99| 人人妻人人看人人澡| 老汉色∧v一级毛片| 日本 av在线| 搡老妇女老女人老熟妇| 99久国产av精品| 国产伦一二天堂av在线观看| 91av网一区二区| 久久久国产精品麻豆| 别揉我奶头~嗯~啊~动态视频| www.熟女人妻精品国产| 欧美最新免费一区二区三区 | 夜夜看夜夜爽夜夜摸| 国产主播在线观看一区二区| 久久精品91无色码中文字幕| 国产综合懂色| 欧美日韩黄片免| 亚洲,欧美精品.| 国产精品久久视频播放| 国产又黄又爽又无遮挡在线| 动漫黄色视频在线观看| 国产成年人精品一区二区| 免费看十八禁软件| 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 搡女人真爽免费视频火全软件 | 日日摸夜夜添夜夜添小说| 九九在线视频观看精品| 久久久久久久久大av| 一个人观看的视频www高清免费观看| 久久草成人影院| 91九色精品人成在线观看| 成人国产一区最新在线观看| 亚洲av一区综合| 99久久精品国产亚洲精品| 一区二区三区国产精品乱码| 亚洲自拍偷在线| 18禁在线播放成人免费| 欧美日韩亚洲国产一区二区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成人中文字幕在线播放| 夜夜爽天天搞| 国内久久婷婷六月综合欲色啪| 老司机在亚洲福利影院| 男女做爰动态图高潮gif福利片| 69av精品久久久久久| 色综合欧美亚洲国产小说| 亚洲最大成人中文| 国产伦精品一区二区三区视频9 | 婷婷亚洲欧美| 久久久久久久久大av| www.www免费av| 天堂影院成人在线观看| 人妻夜夜爽99麻豆av| 日本精品一区二区三区蜜桃| 久久久国产成人免费| 国产探花在线观看一区二区| av天堂在线播放| 亚洲性夜色夜夜综合| 非洲黑人性xxxx精品又粗又长| 亚洲国产欧美网| 久久久久亚洲av毛片大全| 天天添夜夜摸| 午夜福利成人在线免费观看| 亚洲成人免费电影在线观看| 97超视频在线观看视频| 女警被强在线播放| 一个人看视频在线观看www免费 | 午夜免费男女啪啪视频观看 | 亚洲人成网站在线播| 国产精品99久久99久久久不卡| 婷婷六月久久综合丁香| 精品一区二区三区视频在线 | а√天堂www在线а√下载| 99久国产av精品| 亚洲av五月六月丁香网| 亚洲,欧美精品.| tocl精华| 久久久久免费精品人妻一区二区| 久久婷婷人人爽人人干人人爱| 性色avwww在线观看| 久久人妻av系列| 国产99白浆流出| www.999成人在线观看| 亚洲欧美日韩高清专用| 亚洲中文日韩欧美视频| 亚洲美女黄片视频| 日韩欧美精品免费久久 | 色哟哟哟哟哟哟| 亚洲内射少妇av| 一进一出抽搐动态| 成人永久免费在线观看视频| 色综合亚洲欧美另类图片| www.www免费av| 人妻丰满熟妇av一区二区三区| 一进一出抽搐gif免费好疼| 久久久久久国产a免费观看| 国产麻豆成人av免费视频| 在线国产一区二区在线| 美女免费视频网站| 免费观看的影片在线观看| 欧美一级a爱片免费观看看| 悠悠久久av| 亚洲av电影在线进入| 国产97色在线日韩免费| 三级男女做爰猛烈吃奶摸视频| 宅男免费午夜| 97超级碰碰碰精品色视频在线观看| 国产探花极品一区二区| 18禁国产床啪视频网站| 在线免费观看不下载黄p国产 | 国产激情欧美一区二区| 国产av在哪里看| 一夜夜www| 99热这里只有精品一区| 亚洲成人久久性| 亚洲av日韩精品久久久久久密| 成人鲁丝片一二三区免费| www日本在线高清视频| 国产高清激情床上av| 天天躁日日操中文字幕| 国产激情欧美一区二区| 成人三级黄色视频| 男插女下体视频免费在线播放| 怎么达到女性高潮| 日韩欧美国产一区二区入口| 狠狠狠狠99中文字幕| 久久久久精品国产欧美久久久| 精品一区二区三区视频在线观看免费| 亚洲av一区综合| 亚洲国产精品久久男人天堂| 国产精品综合久久久久久久免费| 女人十人毛片免费观看3o分钟| 成人18禁在线播放| 亚洲 国产 在线| 91在线精品国自产拍蜜月 | 色综合站精品国产| 亚洲avbb在线观看| 国产高清视频在线播放一区| 日本免费a在线| 法律面前人人平等表现在哪些方面| 日韩中文字幕欧美一区二区| 99国产综合亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 精品不卡国产一区二区三区| 麻豆国产97在线/欧美| 欧美丝袜亚洲另类 | 日韩欧美在线乱码| 国产免费av片在线观看野外av| 欧美成狂野欧美在线观看| 一区二区三区国产精品乱码| 久久99热这里只有精品18| 亚洲欧美日韩高清专用| 激情在线观看视频在线高清| 国产一区二区激情短视频| 最新在线观看一区二区三区| 不卡一级毛片| 亚洲aⅴ乱码一区二区在线播放| 国产色爽女视频免费观看| 午夜福利高清视频| 嫩草影院入口| 色综合婷婷激情| 在线观看一区二区三区| 青草久久国产| 熟妇人妻久久中文字幕3abv| 亚洲av免费在线观看| 欧美成狂野欧美在线观看| 热99re8久久精品国产| 88av欧美| 亚洲精品国产精品久久久不卡| 十八禁网站免费在线| 欧美在线一区亚洲| 99在线视频只有这里精品首页| 国产av麻豆久久久久久久| 日本与韩国留学比较| 男女视频在线观看网站免费| 99热精品在线国产| 色综合亚洲欧美另类图片| 国产成人av教育| 国内精品一区二区在线观看| 欧美一区二区国产精品久久精品| 国产在视频线在精品| 亚洲av成人不卡在线观看播放网| 两个人视频免费观看高清| 非洲黑人性xxxx精品又粗又长| 757午夜福利合集在线观看| 少妇人妻一区二区三区视频| 免费人成视频x8x8入口观看| 搡老熟女国产l中国老女人| 色综合亚洲欧美另类图片| 国产成人av教育| 婷婷精品国产亚洲av| 国产精品日韩av在线免费观看| 99久久综合精品五月天人人| 在线观看av片永久免费下载| 欧美在线一区亚洲| 中文资源天堂在线| 亚洲精品美女久久久久99蜜臀| 久久精品国产99精品国产亚洲性色| 亚洲人成网站在线播放欧美日韩| 国产黄色小视频在线观看| 观看免费一级毛片| av欧美777| 麻豆国产av国片精品| 久久精品国产亚洲av涩爱 | 亚洲七黄色美女视频| 淫妇啪啪啪对白视频| 亚洲自拍偷在线| 欧美成狂野欧美在线观看| 老司机深夜福利视频在线观看| 国产精品一及| 少妇的丰满在线观看| 国产午夜福利久久久久久| 18禁美女被吸乳视频| 日本一本二区三区精品| 国产主播在线观看一区二区| 欧美中文日本在线观看视频| 日韩人妻高清精品专区| 亚洲av免费高清在线观看| 欧美绝顶高潮抽搐喷水| 国产精品日韩av在线免费观看| 欧美中文综合在线视频| 美女大奶头视频| 久久国产精品影院| 亚洲精品国产精品久久久不卡| 国产精品久久电影中文字幕| 最好的美女福利视频网| 麻豆成人午夜福利视频| 欧美日韩精品网址| 51午夜福利影视在线观看| 久久久精品大字幕| 欧美午夜高清在线| 亚洲国产中文字幕在线视频| 久久国产乱子伦精品免费另类| 欧美成狂野欧美在线观看| 欧美精品啪啪一区二区三区| 亚洲av日韩精品久久久久久密| 级片在线观看| 午夜福利欧美成人| 欧美3d第一页| 国产爱豆传媒在线观看| 久久性视频一级片| 黄色日韩在线| av视频在线观看入口| 一级毛片高清免费大全| 国产欧美日韩精品亚洲av| 熟女少妇亚洲综合色aaa.| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 欧美不卡视频在线免费观看| 一个人免费在线观看电影| 人妻丰满熟妇av一区二区三区| 亚洲片人在线观看| 狂野欧美激情性xxxx| 久久久久久人人人人人| 色视频www国产| 欧美黑人欧美精品刺激| 老司机福利观看| 久久精品人妻少妇| 国产精品嫩草影院av在线观看 | 久久久久国内视频| 日韩精品中文字幕看吧| 一二三四社区在线视频社区8| 国产精品久久电影中文字幕| 精华霜和精华液先用哪个| 亚洲成人久久爱视频| 国产精品久久久久久久电影 | 亚洲精品一卡2卡三卡4卡5卡| 亚洲av电影在线进入| 老司机福利观看| 黄色成人免费大全| 亚洲国产精品成人综合色| 91九色精品人成在线观看| 国产av不卡久久| 日本三级黄在线观看| 亚洲真实伦在线观看| 久久久色成人| 亚洲国产高清在线一区二区三| 丝袜美腿在线中文| 99国产精品一区二区三区| 丰满乱子伦码专区| 看片在线看免费视频| 深夜精品福利| 观看免费一级毛片| 美女黄网站色视频| 制服人妻中文乱码| 久久久成人免费电影| 婷婷丁香在线五月| 亚洲成人中文字幕在线播放| 99久久精品国产亚洲精品| 欧美一区二区精品小视频在线| 一本久久中文字幕| 亚洲中文字幕日韩| 黄片大片在线免费观看| 国产老妇女一区| 欧美zozozo另类| 国产精品99久久久久久久久| 极品教师在线免费播放| 一级毛片高清免费大全| 欧美日韩中文字幕国产精品一区二区三区| 免费人成在线观看视频色| 18禁国产床啪视频网站| 国产aⅴ精品一区二区三区波| 亚洲片人在线观看| 法律面前人人平等表现在哪些方面| 少妇人妻一区二区三区视频| 又爽又黄无遮挡网站| 午夜福利在线在线| av片东京热男人的天堂| 久久人人精品亚洲av| 99热6这里只有精品| 身体一侧抽搐| av女优亚洲男人天堂| 在线天堂最新版资源| 啪啪无遮挡十八禁网站| 亚洲av五月六月丁香网| 欧美成人a在线观看| 欧美激情久久久久久爽电影| 99国产精品一区二区蜜桃av| 老熟妇乱子伦视频在线观看| 老鸭窝网址在线观看| 日韩国内少妇激情av| 亚洲狠狠婷婷综合久久图片| 午夜日韩欧美国产| 午夜福利18| 久9热在线精品视频| 国产精品久久久久久久电影 | 在线观看av片永久免费下载| 国产麻豆成人av免费视频| 日韩欧美 国产精品| 免费观看的影片在线观看| 欧美区成人在线视频| 精品久久久久久久末码| 又紧又爽又黄一区二区| 桃红色精品国产亚洲av| 亚洲中文字幕日韩| 久久久色成人| 久久久久久九九精品二区国产| 男人和女人高潮做爰伦理| 中文亚洲av片在线观看爽| 中亚洲国语对白在线视频| 国产乱人伦免费视频| 日韩av在线大香蕉| 欧美日韩精品网址| 中亚洲国语对白在线视频| 国语自产精品视频在线第100页| av视频在线观看入口| 亚洲欧美日韩东京热| 超碰av人人做人人爽久久 | 国产av不卡久久| 人妻丰满熟妇av一区二区三区| 日本精品一区二区三区蜜桃| 成人欧美大片| 国产精品爽爽va在线观看网站| 一区二区三区国产精品乱码| 乱人视频在线观看| 伊人久久精品亚洲午夜| 亚洲aⅴ乱码一区二区在线播放| 成人高潮视频无遮挡免费网站| 69人妻影院| www日本在线高清视频| 观看美女的网站| 一个人看视频在线观看www免费 | 日韩欧美 国产精品| av女优亚洲男人天堂| 久久久久久人人人人人| 很黄的视频免费| 日本五十路高清| aaaaa片日本免费| 国产精品久久久久久亚洲av鲁大| 国产97色在线日韩免费| 日韩 欧美 亚洲 中文字幕| 91麻豆av在线| 国产免费av片在线观看野外av| 精品久久久久久久久久久久久| 国产欧美日韩一区二区三| 国产综合懂色| 国产成人福利小说| 色综合婷婷激情| 舔av片在线| netflix在线观看网站| 国产精品精品国产色婷婷| 国产伦一二天堂av在线观看| 69av精品久久久久久| 国产精品98久久久久久宅男小说| 日本一二三区视频观看| 免费在线观看影片大全网站| 18禁黄网站禁片午夜丰满| 亚洲精品在线观看二区| 婷婷丁香在线五月| 他把我摸到了高潮在线观看| 亚洲欧美日韩高清在线视频| 亚洲精品国产精品久久久不卡| 少妇的逼水好多| 精品久久久久久久毛片微露脸| 香蕉久久夜色| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品一卡2卡三卡4卡5卡| 性色avwww在线观看| 亚洲中文字幕一区二区三区有码在线看| 日韩亚洲欧美综合| 亚洲精品在线观看二区| 老司机福利观看| 久久精品国产亚洲av香蕉五月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一二三四社区在线视频社区8| 国产欧美日韩精品一区二区| 国内毛片毛片毛片毛片毛片| 精品福利观看| 一本一本综合久久| 亚洲精品在线观看二区| 色av中文字幕| 18禁国产床啪视频网站| 国产伦精品一区二区三区四那| 亚洲无线在线观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品成人综合色| 精品不卡国产一区二区三区| 成人性生交大片免费视频hd| 国产三级中文精品| 欧美成人免费av一区二区三区| 久久婷婷人人爽人人干人人爱| 久久国产精品人妻蜜桃| 琪琪午夜伦伦电影理论片6080| 久久久久免费精品人妻一区二区| 99久久无色码亚洲精品果冻| 色哟哟哟哟哟哟| 非洲黑人性xxxx精品又粗又长| 亚洲精品日韩av片在线观看 | 人人妻人人看人人澡| 好男人电影高清在线观看| 在线十欧美十亚洲十日本专区| 国产精品爽爽va在线观看网站| 国产探花在线观看一区二区| 久久中文看片网| 亚洲男人的天堂狠狠| 午夜免费成人在线视频| 狂野欧美白嫩少妇大欣赏| 我的老师免费观看完整版| 亚洲成av人片免费观看| 小说图片视频综合网站| 国产精品自产拍在线观看55亚洲| 国产精品电影一区二区三区| 国产精品一及| 国内久久婷婷六月综合欲色啪| 国产一区二区在线观看日韩 | 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看| 精品一区二区三区视频在线 | 99国产精品一区二区蜜桃av| 久久久久免费精品人妻一区二区| 国产成人福利小说| 国产色婷婷99| 色吧在线观看| 亚洲不卡免费看| av国产免费在线观看| 最新在线观看一区二区三区| 嫩草影视91久久| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 俺也久久电影网| 变态另类成人亚洲欧美熟女| 国产精品久久久久久人妻精品电影| 国产伦一二天堂av在线观看| 亚洲av电影在线进入| 天堂av国产一区二区熟女人妻| 变态另类丝袜制服| 午夜精品一区二区三区免费看| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 国产私拍福利视频在线观看| 国产欧美日韩精品一区二区| 男女之事视频高清在线观看| 日本成人三级电影网站| 三级毛片av免费| av国产免费在线观看| 两个人看的免费小视频| 黄色视频,在线免费观看| 国产黄片美女视频| 别揉我奶头~嗯~啊~动态视频| 日韩欧美精品v在线| 尤物成人国产欧美一区二区三区| 国产一区二区亚洲精品在线观看| 国产探花极品一区二区| 毛片女人毛片| 日韩欧美精品v在线| 一本精品99久久精品77| 成人特级av手机在线观看| 岛国在线观看网站| 亚洲成av人片免费观看| 精品久久久久久,| 国产 一区 欧美 日韩| 色综合站精品国产| 一个人看的www免费观看视频| 精品一区二区三区视频在线 | 最新美女视频免费是黄的| 综合色av麻豆| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩一级在线毛片| eeuss影院久久| 淫妇啪啪啪对白视频| 香蕉丝袜av| 亚洲中文字幕日韩| 两人在一起打扑克的视频| 国产精品久久久人人做人人爽| 真人做人爱边吃奶动态| 免费搜索国产男女视频| 国内揄拍国产精品人妻在线| 国产高清三级在线| 叶爱在线成人免费视频播放| 日韩中文字幕欧美一区二区| 很黄的视频免费| 色综合亚洲欧美另类图片| 精品久久久久久久人妻蜜臀av| 91久久精品国产一区二区成人 | 母亲3免费完整高清在线观看| 国产免费一级a男人的天堂| 久久精品国产99精品国产亚洲性色| 97超视频在线观看视频| 国产精品亚洲美女久久久| 免费观看的影片在线观看| 亚洲自拍偷在线| 国产 一区 欧美 日韩| 听说在线观看完整版免费高清| 欧美区成人在线视频| xxxwww97欧美| 91麻豆精品激情在线观看国产| 日韩欧美在线乱码| 性色avwww在线观看| 最好的美女福利视频网| 日本一本二区三区精品| 亚洲欧美一区二区三区黑人| 12—13女人毛片做爰片一| 成熟少妇高潮喷水视频| 国产一区二区亚洲精品在线观看| 黄色日韩在线| 久久久久性生活片| 18禁国产床啪视频网站| 内射极品少妇av片p| 亚洲成av人片免费观看| 国产一区在线观看成人免费| 日日摸夜夜添夜夜添小说| 久久6这里有精品| 亚洲专区中文字幕在线| 欧洲精品卡2卡3卡4卡5卡区| 1024手机看黄色片| 日本免费一区二区三区高清不卡| bbb黄色大片| 男女午夜视频在线观看| 熟女少妇亚洲综合色aaa.| 久久中文看片网| 99精品在免费线老司机午夜| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区| 久久香蕉国产精品| 欧美日本视频| 国产伦人伦偷精品视频| 国产乱人伦免费视频| 嫁个100分男人电影在线观看| 亚洲电影在线观看av| 欧美日韩福利视频一区二区| 人人妻,人人澡人人爽秒播| 中出人妻视频一区二区| 久久婷婷人人爽人人干人人爱| 男女那种视频在线观看| 国产欧美日韩精品亚洲av| 欧美高清成人免费视频www| 亚洲精品粉嫩美女一区| xxx96com| 亚洲精品色激情综合| 波多野结衣巨乳人妻| 国产成年人精品一区二区| 亚洲第一电影网av| 亚洲熟妇中文字幕五十中出| 99久久九九国产精品国产免费| 国产中年淑女户外野战色| 美女黄网站色视频| 久久精品国产亚洲av涩爱 | 男女视频在线观看网站免费| 午夜免费成人在线视频| 久久久精品大字幕| 母亲3免费完整高清在线观看| 亚洲精品一区av在线观看| 人妻丰满熟妇av一区二区三区| 老鸭窝网址在线观看| 国产三级在线视频| 国内精品久久久久精免费| 天天躁日日操中文字幕| 国产野战对白在线观看| 国内精品久久久久精免费| 嫩草影院入口| 亚洲欧美日韩无卡精品| 免费一级毛片在线播放高清视频| 99riav亚洲国产免费| 3wmmmm亚洲av在线观看| 欧美高清成人免费视频www| tocl精华| 久久久成人免费电影| 久久精品国产亚洲av涩爱 | 国产伦一二天堂av在线观看| 国产真实乱freesex| 日韩欧美一区二区三区在线观看| 国产探花极品一区二区|