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

    基于BSA-seq技術(shù)定位調(diào)控西瓜紅色果肉形成基因

    2025-07-11 00:00:00張曼王瑩瑩劉金秋婁麗娜徐建羊杏平徐錦華
    江蘇農(nóng)業(yè)學(xué)報 2025年5期

    中圖分類號:S642.9 文獻標識碼:A 文章編號: 1000-4440(2025)05-0952-08

    Abstract:Flesh colorisanimportant trait thatdetermines thequalityof watermelon fruit.Inthisstudy,weused the pink-fleshed watermelon MW119 and the red-fleshed watermelon MW120 as parents to construct an F2 segregating population. Fromthis population,weselected individuals withredflesh and pink flesh phenotypes tocreate progenybulks.Wethen performed whole-genome resequencing with a coverage depth of 30 1× on both the parental and progeny bulks using BSA-seq technology.By using the Δ (SNP-index)method,we identified 12 intervals significantlyassociated with the red flesh trait on chromosomes2,4,5,6,7,9,nd1,withatotallengthof56.66Mb.Furtherfunctionalaotationandexpresionaalyis wereconductedon116genes withinthesignificantlyassociated intervalsthatcariednon-synonymousSNPandInDelframeshift mutations,ultimatelyleading totheselectionoffivekeycandidategenes:the WRKY transcriptionfactorgene Cla97C02G028430,the o -methyltransferase gene Cla97C02G028680,the F-box protein genes Cla97C02G048400 and Cla97C06Gl22530,and the pentatricopeptide repeat protein gene Cla97Co6G122120.These genes playa significant role in regulating the formation ofred fleshinwatermelon.The resultsof thisstudy furtherrefine themolecular mechanisms underlyingtheformation ofredfleshinwatermelonand providea theoretical basis for marker-assisted breeding.

    Key words: watermelon;flesh color; gene map-ping

    西瓜(Citrulluslanatus)是重要的園藝作物,因其色澤鮮艷,果肉甘甜多汁而深受消費者喜愛。隨著經(jīng)濟發(fā)展,人們對西瓜品質(zhì)的要求不斷提高。果肉顏色是衡量西瓜果實品質(zhì)的重要指標,選育果肉顏色鮮艷的品種已成為育種工作者的重點研究方向。研究結(jié)果表明,西瓜果肉顏色與類胡蘿卜素類型和含量有關(guān),紅色果肉是西瓜定向馴化選擇的結(jié)果[1],也是目前市場上最為常見的果肉顏色。紅瓢和粉瓢西瓜富含番茄紅素[;黃瓢西瓜的主要類胡蘿卜素為紫黃質(zhì)和新黃質(zhì)3;橙瓢西瓜的主要成分為前番茄紅素或 β -胡蘿卜素[45];白瓢西瓜則富含六氫番茄紅素和 ζ 胡蘿卜素[6-7]

    研究發(fā)現(xiàn),西瓜果肉顏色由多基因控制[8-9],其遺傳規(guī)律為:紅瓢(Y)對橙襄 (yo) 和淺黃瓢 (y) 為顯性;橙瓢 (yo) 對淺黃 (y) 為顯性;亮黃瓢 (C) 對除白瓢(uf)外的其他果肉顏色均為顯性[10-14]。在果肉顏色基因定位研究方面,研究人員通過構(gòu)建遺傳連鎖圖譜,分別在4號、3號和11號染色體上鑒定到調(diào)控番茄紅色合成的主效數(shù)量性狀位點(QTL)[15-17]。Branham等[4]在1號染色體上定位到調(diào)控 β -胡蘿卜素合成的主效數(shù)量性狀位點(QTL)。Nie等[18]通過重新組裝西瓜基因組,鑒定出15個與果實品質(zhì)相關(guān)的QTL,包含33個CDS區(qū)發(fā)生變異的候選基因,其中位于1號染色體9539129bp的八氫番茄紅素合酶基因Cla97C01G008760第一個外顯子處存在 單核苷酸變異,推測可能與果肉顏色形成有關(guān)。基于泛基因組分析,Sun等[]在400份西瓜種質(zhì)中鑒定到8個與果肉顏色關(guān)聯(lián)的缺失變異,其中4個與胡蘿卜素積累相關(guān)的無參考候選基因(Non-referencecandidategenes)在白色果肉中分布頻率較高。Yi等[14]通過群體遺傳分析將白色果肉基因定位于6號染色體 132.3kb 區(qū)間內(nèi)。Fang等[5]和Jin等[20]分別將調(diào)控橙色果肉形成基因Clorf定位于10號染色體 82.51kb 區(qū)間內(nèi)和1號染色體 39.08kb 區(qū)間內(nèi)。Liu等1將調(diào)控金黃色果肉形成基因精細定位于1號染色體 39.08kb 區(qū)間內(nèi)。Park等[13]將調(diào)控紅黃混合果肉(ICY)形成基因定位于2號染色體27.60\~27.88Mb,并開發(fā)了與ICY 性狀共分離的標記 M7 。Wang等[21]將調(diào)控紅色果肉形成基因定位于4號染色體 41.233kb 區(qū)間內(nèi) ?8887652~8928885bp), 。Li等[22]將調(diào)控紅色果肉形成基因定位于6號染色體 40.288kb 區(qū)間內(nèi)(21851079\~21891367bp)。果肉顏色是一個復(fù)雜的數(shù)量性狀,涉及類胡蘿卜素、番茄紅素等色素的合成和積累,仍需進一步開展基因定位研究以挖掘其更多調(diào)控位點,為多種顏色西瓜新品種的選育提供分子標記。本研究擬以紅色果肉西瓜MW120和粉色果肉西瓜MW119為親本構(gòu)建 F2 分離群體,通過極端混池全基因組重測序篩選調(diào)控紅色果肉形成基因,旨在為西瓜種質(zhì)改良提供理論依據(jù)。

    1材料與方法

    1.1 試驗材料

    本研究以紅色果肉西瓜MW120和粉色果肉西瓜MW119為親本,通過人工雜交構(gòu)建 F1 分離群體,包含130個單株。該群體于2022年在江蘇東臺種植,采用常規(guī)田間管理。果肉顏色由4名研究人員在相同光照條件下通過肉眼進行鑒定,篩選出果肉顏色分離的單株。

    1.2 親本池和子代池的構(gòu)建

    取5株紅色果肉西瓜MW120的基因組DNA,等量混合構(gòu)建MW120親本池;取5株粉色果肉西瓜MW119的基因組DNA,等量混合構(gòu)建MW119親本池。從 F2 群體中選取20株紅色果肉單株的基因組DNA,等量混合構(gòu)建紅色果肉西瓜子代池;選取20株粉色果肉單株的基因組DNA,等量混合構(gòu)建粉色果肉西瓜子代池。利用超微量紫外分光光度計(型號NanoDrop8000,ThermoScientific公司產(chǎn)品)測定DNA濃度。

    1.3 DNA文庫的BSA-Seq測序

    親本池和子代池的文庫構(gòu)建和重測序由廣州基迪奧生物科技有限公司完成。DNA文庫在IlluminaHiSeq×10 平臺上進行雙末端測序(PE150),基因組測序深度為 30× 。對測序原始數(shù)據(jù)(Rawreads)進行質(zhì)量控制和過濾,得到高質(zhì)量數(shù)據(jù)(Cleanreads)。使用BWA軟件將高質(zhì)量數(shù)據(jù)與97103西瓜參考基因組(http://cucurbitgenomics.org/)進行比對。

    1.4 SNP檢測與關(guān)聯(lián)分析

    利用GATKUnifiedGenotyper(v3.5)進行SNP檢測,并使用ANNOVAR進行功能注釋。采用 Δ (SNP-index)法進行關(guān)聯(lián)分析,以 99% 置信水平篩選顯著區(qū)間。通過Gene Ontology(GO)和 KEGG Pathway顯著性富集分析對候選基因進行功能注釋。

    1.5 熒光定量PCR分析

    以不同發(fā)育階段(授粉后 10d,18d,28d,34d)

    的西瓜果肉為材料,使用RNApurePlantKit(DNa-seI)試劑盒(康為世紀生物科技有限公司產(chǎn)品)提取總RNA,利用BU-SuperScriptRTKit試劑盒(Bi-ouniquer公司產(chǎn)品)合成cDNA。利用SuperRealPreMix[天根生化科技(北京)有限公司產(chǎn)品]進行熒光定量PCR分析。以 I8SrRNA 為內(nèi)參基因,采用 2-ΔΔCt 法計算相對表達量[23]。每個樣品設(shè)3次生物學(xué)重復(fù)。

    2 結(jié)果與分析

    2.1 測序數(shù)據(jù)質(zhì)量評估

    如表1所示,原始序列(RawReads)過濾后得到高質(zhì)量有效序列71.44G,測序序列中質(zhì)量值 ? 20的堿基占比( Q20) 為 95.35%~96.06% ,測序序列中質(zhì)量值 ?30 的堿基占比( Q30 )為 88.11%~ 89.84% , G+C 含量為 35.20%~35.61% ,比對率為99.01%~99.16% 。測序深度為 30× 時,參考基因組中至少被1個測序短序列(reads)覆蓋的位點占比為 98.21%~98.26% ,至少被4個測序短序列(reads)覆蓋的位點占比為 98.13%~98.20% ,至少被10個測序短序列(reads)覆蓋的位點占比為 97.94% 298.07% ,至少被20個測序短序列(reads)覆蓋的位點占比為 94.58%~97.63% 。綜上所述,所有樣本的測序數(shù)據(jù)量充足,測序序列質(zhì)量較高。

    表1測序數(shù)據(jù)質(zhì)量Table1 Sequencing dataquality
    :測序序列中質(zhì)量值 的堿基所占比例; :測序序列中質(zhì)量值 的堿基所占比例; :參考基因組中至少被1個測序短序列(reads)覆蓋的位點占比; 4× :參考基因組中至少被4個測序短序列(reads)覆蓋的位點占比; 10× :參考基因組中至少被10個測序短序列(reads)覆蓋的位點占比; 20× :參考基因組中至少被20個測序短序列(reads)覆蓋的位點占比。

    2.2 SNP位點分析

    結(jié)果表明,4個樣本池共檢測到193187個SNP位點,其中粉色果肉子代池中SNP位點數(shù)量最多,MW119親本池(粉色果肉)中SNP位點數(shù)量最少。親本間共有2958個非同義突變SNP位點,紅色果肉西瓜子代池和粉色果肉西瓜子代池之間共有4055個非同義突變SNP位點。SNP位點如表2所示,親本池之間和子代混池之間同一類型的SNP位點數(shù)量相當(dāng),變異位點主要集中于基因間區(qū)、內(nèi)含子、基因上游和基因下游。

    2.3 基于 Δ (SNP-index)法的關(guān)聯(lián)分析

    Δ (SNP-index)法是一種通過計算混池間基因型頻率的差異進行標記關(guān)聯(lián)分析的方法。如圖1所示,基于 Δ (SNP-index)方法( 99% 置信水平)鑒定到12個顯著關(guān)聯(lián)區(qū)間,總長度 56.66Mb ,分布于7條染色體。西瓜2號染色體包含5個區(qū)間( 1.04~4.34 Mb.8.24~10.36Mb.10.97~13.07Mb.14.51~16.98 Mb、27.41~37.92Mb);4號染色體包含1個區(qū)間(3.68~16.84Mb) ;5號染色體包含1個區(qū)間(20.77~23.57Mb);6號染色體包含2個區(qū)間(0.96~6.33Mb ,22.69~24.86Mb );7號染色體包含1個區(qū)間( 7.60~9.60Mb 0;9號染色體包含1個區(qū)間( 19.74~26.50Mb );10號染色體包含1個區(qū)間( 14.14~18.04Mb) 。候選區(qū)間內(nèi)共包含2822個基因,其中116個基因存在非同義突變/移碼突變。

    表2SNP位點信息Table2 SNPsiteinformation
    圖1基于 Δ (SNP-index)法的顯著關(guān)聯(lián)區(qū)間 Fig.1Significantlyassociated intervalsbased on the Δ (SNP-index)method

    黑色曲線為對 Δ (SNP-index)進行擬合后得到的曲線。

    2.4候選基因的GO和KEGG富集分析

    為了進一步解析候選基因功能,對其進行GO富集分析,對候選基因及其編碼產(chǎn)物的屬性按照細胞組成、分子功能和生物過程進行分類注釋。如圖2所示,候選基因主要富集在細胞組成(346個基因)和生物過程(342個基因)中,參與分子功能的候選基因(28個基因)最少。其中,脂質(zhì)生物合成過程(GO:0008610)、細胞間連接(G0:0005911)、細胞連接(G0:0030054)是候選基因富集程度最高的功能分類。進一步對候選基因進行KEGG通路富集分析,發(fā)現(xiàn)植物激素信號轉(zhuǎn)導(dǎo)、亞油酸代謝為候選基因主要富集通路。

    2.5 候選基因的功能注釋

    如表3所示,對關(guān)聯(lián)區(qū)間內(nèi)116個具有SNP非同義突變和InDel移碼突變的基因進行進一步的功能注釋和相關(guān)文獻查閱,從中篩選出17個可能與西瓜紅色果肉相關(guān)的候選基因,分別位于2號、4號和6號染色體上。其中Cla97C02G027340編碼紫黃質(zhì)脫環(huán)氧化酶,該酶可以催化類胡蘿卜素途徑中紫黃質(zhì)、單環(huán)氧花藥黃質(zhì)和玉米黃質(zhì)之間的相互轉(zhuǎn)化,從而調(diào)控果實成熟過程中 β 胡蘿素的積累。Cla97C02G027870、Cla97C02G028430 和 Cla97C04G070570分別編碼MYB轉(zhuǎn)錄因子、WRKY轉(zhuǎn)錄因子和MADS-box轉(zhuǎn)錄因子,這些轉(zhuǎn)錄因子參與果實顏色的調(diào)控。Cla97C02G028400編碼引導(dǎo)蛋白(DIR)。Cla97C02G028680編碼鄰甲基轉(zhuǎn)移酶(OMT),該酶與花青素的生物合成相關(guān)。Cla97C02G028640和Cla97C02G029020編碼尿苷二磷酸-糖基轉(zhuǎn)移酶。Cla97C02G048400、Cla97C06G122510、Cla97C06G122520 和Cla97C06G122530編碼F-box 蛋白。Cla97C02G049320、Cla97C04G070020、Cla97C04G070590 和 Cla97C06G122120編碼三角狀五肽重復(fù)蛋白。這些基因編碼區(qū)均發(fā)生了SNP非同義突變,可能是調(diào)控紅色果肉形成的關(guān)鍵基因。

    2.6候選基因在西瓜果實發(fā)育過程中的表達

    如圖3所示,17個候選基因在西瓜果實發(fā)育過程中呈現(xiàn)不同的表達模式。在整個果實發(fā)育期間,Cla97C02G028430、 Cla97C06G122520、 Cla97C02G048400、Cla97C06G122530和Cla97C06G122120的相對表達量均保持較高水平。在授粉后18d,Cla97C04G070590的相對表達量降低,隨后又逐漸升高。在授粉后18d,Cla97C02G028680的相對表達量達到最大值。隨著果實發(fā)育,Cla97C04G070570的相對表達量呈降低趨勢。在整個果實發(fā)育期間,Cla97C02G027870、Cla97C02G028400、Cla97C02G028640、Cla97C02G028650、 Cla97C06G122510、Cla97C02G049320的相對表達量均保持在較低水平。表明Cla97C02G028430、Cla97C06G122520、Cla97C02G048400、Cla97C06G122530和Cla97C06G122120為調(diào)控西瓜紅色果肉形成的關(guān)鍵基因。

    3 討論與結(jié)論

    紅色果肉因其廣受消費者青睞而成為市場主流。目前關(guān)于紅色果肉形成的遺傳機制研究已取得重要進展,但不同研究人員基于不同遺傳群體獲得的結(jié)果存在差異。Wang等[2i]將紅色果肉基因精細定位到4號染色體 區(qū)間內(nèi)(8886138~8926873bp ),番茄紅素 β -環(huán)化酶(LCYB)基因ClaO05011為關(guān)鍵候選基因。Li等[22]在6號染色體40kb 的區(qū)間內(nèi),發(fā)現(xiàn)番茄紅素生物合成相關(guān)基因GO:0045165:細胞命運決定;GO:0006511:泛素依賴性蛋白質(zhì)分解代謝過程;"GO;0035966":對拓撲錯誤蛋白質(zhì)的響應(yīng);"GO;0006633":脂肪酸生物合成過程;GO:0072330:單羧酸生物合成過程;GO:008610;脂質(zhì)生物合成過程;GO:0016701:氧化還原酶活性(作用于單供體并摻入分子氧);"G0:0005942":磷脂酰肌醇3-激酶復(fù)合體;GO:0031010:ISWI型復(fù)合體;GO:0035032:II類磷脂酰肌醇3-激酶復(fù)合體;GO:0070603:SWI/SNF超家族型復(fù)合體;GO:0019898:膜的外在成分;GO:0015629:肌動蛋白細胞骨架; G0:0098687 :染色體區(qū)域;GO:0030119:AP型膜被適配體復(fù)合體;GO:0031225:膜錨定成分;GO:0009527:質(zhì)體外膜;GO:005694;染色體;GO:005911;細胞間連接;GO:0030054;細胞連接。Q值:多重假設(shè)檢驗校正后的 P 值。

    表3候選基因的功能注釋Table3 Functional annotation of candidate genes
    圖3果實不同發(fā)育過程中候選基因相對表達量Fig.3Relative expression levels of candidate genes during different fruit development stages

    Cla018767、Cla018768、Cla018769、Cla018770和ClaOI877I 。Zhang等[24發(fā)現(xiàn),三角狀五肽重復(fù)蛋白基因 Cla97CO6GI22I20 可能通過調(diào)控類胡蘿卜素代謝影響果肉顏色,并將調(diào)控西瓜粉色果肉形成的基因定位于6號染色體 55.26kb 區(qū)間內(nèi),最終確定Cla97C06G122120為關(guān)鍵候選基因。本研究以粉色果肉西瓜MW119和紅色果肉西瓜MW120為親本,構(gòu)建 F2 分離群體,利用BSA-seq方法將調(diào)控紅色果肉形成的基因定位于2號、4號、5號、6號、7號、9號、10號染色體上。進一步對關(guān)聯(lián)區(qū)間內(nèi)116個攜帶非同義突變SNP和InDel移碼突變的基因進行功能注釋和表達分析,最終篩選出6個關(guān)鍵候選基因:WRKY轉(zhuǎn)錄因子基因Cla97C02G028430、鄰甲基轉(zhuǎn)移酶基因Cla97C02G028680、F-box蛋白基因

    Cla97C02G048400和Cla97C06G122530,以及三角狀五肽重復(fù)蛋白基因 Cla97CO6GI22I20

    目前已報道參與西瓜果肉顏色調(diào)控的基因主要包括:八氫番茄紅素合酶基因 (PSY)[25] 、番茄紅素β -環(huán)化酶基因(LCYB)[1]、9-順式環(huán)氧類胡蘿卜素雙加氧酶基因(NCED)[5,14,21]、類胡蘿 h 素異構(gòu)酶基因(ZISO)[26]等類胡蘿卜素通路相關(guān)基因以及MYB、WRKY和MADS-box 等轉(zhuǎn)錄因子基因[11,27-28]和質(zhì)體發(fā)育相關(guān)基因[29]。本研究首先篩選了編碼區(qū)存在非同義突變SNP的基因,進而通過熒光定量PCR分析其在果實發(fā)育過程中的表達模式,最終鑒定出5個可能參與紅色果肉形成的候選基因。WRKY轉(zhuǎn)錄因子可通過協(xié)同調(diào)控葉綠素降解和類胡蘿卜素的合成促進柑橘果實著色[30],或通過乙烯信號轉(zhuǎn)導(dǎo)途徑調(diào)控番茄果實顏色的形成[31」。本研究中,在果實發(fā)育中期,Cla97C02G028430的相對表達量保持較高水平。Liu等[27]發(fā)現(xiàn),WRYK轉(zhuǎn)錄因子參與果實成熟的調(diào)控。而其與果肉顏色轉(zhuǎn)變的相關(guān)性仍需進行進一步探究。Cla97C02G028680編碼鄰甲基轉(zhuǎn)移酶(OMT),在果實發(fā)育中期, Cla97CO2GO28680 的相對表達量保持較高水平。已有研究結(jié)果表明,鄰甲基轉(zhuǎn)移酶基因表達量與花色苷的積累量呈正相關(guān),而花色苷是決定葡萄[32]和辣椒[33果皮的呈色物質(zhì)。Cla97C02G048400和Cla97C06G122530編碼F-box家族蛋白。F-box在真核生物中普遍存在,參與多種生理過程的調(diào)控。研究結(jié)果表明,F(xiàn)-box參與調(diào)控有色小麥中的花青素積累[34]、辣椒果實轉(zhuǎn)色[35]以及甜瓜白色果皮的形成[36]。Cla97C06G122120 編碼三角狀五肽重復(fù)蛋白PPR,在最新發(fā)表的研究結(jié)果中,該基因被鑒定為調(diào)控顏色果肉的候選基因[12.24]。這些候選基因的發(fā)現(xiàn)為深入解析西瓜紅色果肉形成的分子機制提供了重要線索。后續(xù)研究將通過基因編輯等技術(shù)對這些基因進行功能驗證,以期完善西瓜果肉顏色形成的遺傳調(diào)控網(wǎng)絡(luò)。

    參考文獻:

    [1]GUO SG,ZHAO SJ,SUNHH,et al.Resequencing of 414 cultivatedandwildwatermelonaccessionsidentifiesselectionforfruit qualitytraits[J].NatureGenetics,2019,5:1616-1623.

    [2] 李艷杰,康保珊,趙文恩,等.西瓜中類胡蘿卜素的高效液相色 譜法分析[J].中國瓜菜,2005(4):12-15.

    [3] 王楠,劉識,朱子成,等.紅色和橙黃色果肉西瓜番茄紅素 (5):918-926.

    [4] BRANHAM S, VEXLER L, MEIR A, et al. Genetic mapping of a major codominant QTL associated with β -carotene accumulation in watermelon[J].Molecular Breeding,2017,37(12):146.

    [5]FANG XF,LI SL,ZHU ZC,et al.Clorf encodes carotenoid isomerase and regulates orange flesh color in watermelon(Citrullus lanatusL.)[J].Journal of Agricultural andFood Chemistry, 2023,71(42):15445-15455.

    [6]KANG B S, ZHAO W E,HOU Y B,et al. Expression of carotenogenic genes during the development and ripening of watermelon fruit[J]. Scientia Horticulturae,2010,124:368-375.

    [7]GRASSI S,PIRO G,LEE JM,et al. Comparative genomics reveals candidatecarotenoid pathway regulators of ripening watermelon fruit[J].BMC Genomics,2013,14:781.

    [8] 顧衛(wèi)紅,張燕,鄭洪建.西瓜的色遺傳及其環(huán)境調(diào)控作用 [J].上海農(nóng)業(yè)學(xué)報,2002,18(3):43-46.

    [9] MARAGAL S, RAL E S, REDDY D C L. Genetic mapping and meta-analysis identifies several candidate genes for watermelon (Citrullus lanatus)fruit quality traits[J]. Scientia Horticulturae, 2023,308;111545.

    [10]刁衛(wèi)楠,袁平麗,龔成勝,等.西瓜果肉檸檬黃色的遺傳分析和 基因定位[J].中國農(nóng)業(yè)科學(xué),2021,54(18):3945-3958.

    [11]LIU S,GAO ZQ,WANG X Z,et al.Nucleotide variation in the phytoene synthase (ClPsy1) gene contributes to golden flesh in watermelon (Citrullus lanatus L.)[J]. Theoretical and Applied Genetics,2022,135(1) :185-200.

    [12]LIU S,LIUMQ,CAO Y,et al.Identification of chromosome regionand candidate genes for canary-yellow flesh(Cyf) locus in watermelon(Citrullus lanatus)[J].Plant Science,2O23,329: 111594.

    [13]PARK G,SHAHWAR D,JANG G,et al. Identification of a novel locus C2 controlling canary yellow flesh color in watermelons[J]. Frontiers in Genetics,2023,14:1256627.

    [14]YI L C, ZHOU W, ZHANG Y,et al. Genetic mappng of a single nuclear locus determines the white flesh color in watermelon(Citrullus lanatusL.)[J].Frontiers in Plant Science,2023,14: 1090009.

    [15]LIU S,GAO P,WANG X Z,et al. Mapping of quantitative trait loci for lycopene content and fruit traits in Citrulus lanatus[J]. Euphytica,2015,202(3):411-426.

    [16]LIU S,GAO P,ZHUQL,et al.Development of cleaved amplified polymorphic sequence markers and a CAPS-based genetic linkage map in watermelon(Citrullus lanatus (Thunb.)Matsum. and Nakai) constructed using whole-genome re-sequencing data [J].Breeding Science,2016,66(2):244-259.

    [17]FALL LA,PERKINS-VEAZIEP,MA G Y. QTLs associated with flesh quality traits in an elite_elite watermelon population[J].Euphytica,2019,215:30.

    [18]NIE HL,KIM M,LEE S,et al. High-quality genome assembly and genetic mapping reveal a gene regulating flesh color in watermelon(Citrulluslanatus)[J].Frontiers inPlant Science,2023, 14:1142856.

    [19]SUNY,KOUDR,LIY,et al.Pan-genome of Citrullus genus highlights the extent of presence/absence variation during domesticationand selection[J].BMCGenomics,2023,24:332.

    [20]JINBK,JANGG,PARKG,etal.Developmentof a gene-based marker set for orange-colored watermelon flesh with a high β -carotene content[J]. International Journal of Molecular Sciences, 2024,25:210.

    [21]WANG CN,QIAO AH,F(xiàn)ANG XF,et al. Fine mapping of lycopene content and flesh color related gene and development of molecular marker-assisted selection for flesh color in watermelon(Citrullus lanatus)[J].Frontiers in Plant Science,2019,10:1240.

    [22]LIN,SHANGJL,WANGJM,et al.Discovery of the genomic region and candidate genes of the scarlet red flesh color(Yscr)locusinwatermelon(Citrullus LanatusL.)[J].Frontiers inPlant Science,2020,11:116.

    [23]LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔGt method [J].Methods,2001,25(4):402-408.

    [24]ZHANGWC,JIZQ,HUGQ,etal.Clpf encodes pentatricopeptide repeat protein(PPR5)and regulatespink flesh color in watermelon(Citrullus lanatusL.)[J].Theoretical and Applied Genetics,2024,137:126.

    [25]FANG XF,LIU S,GAO P,et al.Expression of CLPAP and CLPSYl in watermelon correlates with chromoplast differentiation, carotenoid accumulation,and flesh color formation[J].Scientia Horticulturae,2020,270:10943.

    [26]ZHANG J,SUNHH,GUO SG,et al.ClZISO mutation leads to photosensitive flesh in watermelon[J].Theoretical and Applied Genetics,2022,135(5):1565-1578.

    [27]LIU YD,XIAO W,LIAO L,et al. A PpEIL2/3-PpNAC1-PpWRKY14 module regulates fruit ripening by modulating ethylene production in peach[J]. Journal of Integrative Plant Biology, 2024,66(11) :2470-2489.

    [28]CAO Y,F(xiàn)ANGXF,LIU S,et al. Isolation and activity analysis of phytoene synthase(ClPsy1) gene promoter of canary-yellow and golden flesh-color watermelon[J].Agronomy,2023,13(4):1080.

    [29]ZHANGJ,GUO SG,RENY,et al.High-level expression of a novel chromoplast phosphate transporter CLPHT4;2 is required for flesh color development in watermelon[J].New Phytologist,2017, 213:1208-1221.

    [30]CHENHY,JIHY,HUANGWK,et al.Transcription factor CrWRKY42 coregulates chlorophyll degradation and carotenoid biosynthesis in citrus[J].PlantPhysiology,2024,195(1):728-744.

    [31]ZHAOWH,LIYH,F(xiàn)AN SZ,etal.The transcription factor WRKY32 affects tomato fruitcolour byregulating YELLOW FRUITED-TOMATO1,a core component of ethylene signal transduction[J]. Journal of Experimental Botany,2021,72(12):4269- 4282.

    [32]蘇玲,楊陽,王鵬飛,等.葡萄中調(diào)控果皮顏色相關(guān)基因的 研究進展[J].安徽農(nóng)業(yè)科學(xué),2018,46(2):19-22,26.

    [33]張寶璽,郭家珍,楊佳梅,等.辣椒綠熟期果色及主要色素含量 的遺傳[J].園藝學(xué)報,1996,23(1):94-96.

    [34]HONGMJ,KIMDY,CHOIH,et al.Isolation and characterization of kelch repeat-containing F-box proteins from colored wheat [J].Molecular Biology Reports,2020,47:1129-1141.

    [35]金威恒,MUMTAZMA,何成堯,等.基于加權(quán)關(guān)聯(lián)網(wǎng)絡(luò)分析的 辣椒(CapsicumannuumL.)果實顏色轉(zhuǎn)變相關(guān)基因的鑒定(英 文)[J].熱帶作物學(xué)報,2022,43(6):1129-1143.

    [36]FEDER A,BURGERJ,GAO S,et al.A kelch domain-containing F-box coding gene negatively regulates flavonoid accumulation in muskmelon[J].Plant Physiology,2015,169(3):1714-1726. (責(zé)仟編輯.成紓賓)

    日本色播在线视频| 综合色丁香网| 免费看不卡的av| av国产精品久久久久影院| 伦精品一区二区三区| 日日爽夜夜爽网站| 高清黄色对白视频在线免费看| 日韩欧美一区视频在线观看| 最近2019中文字幕mv第一页| 美女福利国产在线| 国产亚洲午夜精品一区二区久久| 国产黄频视频在线观看| 久久精品国产亚洲av高清一级| 欧美人与善性xxx| 久久亚洲国产成人精品v| 久久久久久人妻| 国产黄色视频一区二区在线观看| 亚洲一区二区三区欧美精品| 国产深夜福利视频在线观看| 久久影院123| 男女边摸边吃奶| 欧美人与善性xxx| www.自偷自拍.com| 国语对白做爰xxxⅹ性视频网站| 国产成人精品婷婷| 麻豆av在线久日| 久久久久精品人妻al黑| 新久久久久国产一级毛片| 91精品伊人久久大香线蕉| av免费在线看不卡| 欧美精品一区二区大全| av视频免费观看在线观看| 999久久久国产精品视频| 美女福利国产在线| 亚洲,一卡二卡三卡| 久久影院123| 国产综合精华液| 一级毛片电影观看| av又黄又爽大尺度在线免费看| 69精品国产乱码久久久| 午夜福利一区二区在线看| a 毛片基地| 丝袜美足系列| 国产成人免费无遮挡视频| 男女高潮啪啪啪动态图| 视频区图区小说| 亚洲精品自拍成人| 交换朋友夫妻互换小说| 伊人亚洲综合成人网| 丝袜美腿诱惑在线| 黄色一级大片看看| 18+在线观看网站| 久久婷婷青草| 美女中出高潮动态图| 亚洲美女搞黄在线观看| 如何舔出高潮| 999精品在线视频| 激情视频va一区二区三区| 免费在线观看黄色视频的| 人妻一区二区av| 国产乱来视频区| 国产精品 欧美亚洲| 免费高清在线观看日韩| 婷婷成人精品国产| 丰满乱子伦码专区| 成人毛片a级毛片在线播放| 亚洲,一卡二卡三卡| 久久这里只有精品19| 超色免费av| 只有这里有精品99| 免费在线观看黄色视频的| 日韩伦理黄色片| 如日韩欧美国产精品一区二区三区| 成人手机av| 大香蕉久久成人网| 91国产中文字幕| 丝袜人妻中文字幕| 欧美xxⅹ黑人| 午夜福利,免费看| 国产精品二区激情视频| 亚洲欧美色中文字幕在线| 99热全是精品| 国产免费视频播放在线视频| 性少妇av在线| 伦理电影大哥的女人| 成人午夜精彩视频在线观看| 中文字幕制服av| 丝袜脚勾引网站| 啦啦啦视频在线资源免费观看| 免费久久久久久久精品成人欧美视频| 精品国产乱码久久久久久小说| 成人手机av| 人成视频在线观看免费观看| www.av在线官网国产| 美女午夜性视频免费| 中文字幕亚洲精品专区| 午夜福利影视在线免费观看| 啦啦啦视频在线资源免费观看| 高清av免费在线| 中文精品一卡2卡3卡4更新| 精品久久蜜臀av无| av又黄又爽大尺度在线免费看| 一级毛片 在线播放| 免费不卡的大黄色大毛片视频在线观看| 国产麻豆69| 日韩人妻精品一区2区三区| 99国产综合亚洲精品| 性高湖久久久久久久久免费观看| 免费在线观看完整版高清| 国产精品不卡视频一区二区| 18+在线观看网站| 国产 精品1| 黄色毛片三级朝国网站| 亚洲 欧美一区二区三区| 男女国产视频网站| 波多野结衣一区麻豆| 男女啪啪激烈高潮av片| 看免费成人av毛片| 90打野战视频偷拍视频| 欧美最新免费一区二区三区| 99精国产麻豆久久婷婷| 中文字幕另类日韩欧美亚洲嫩草| 国产精品av久久久久免费| 黄频高清免费视频| 亚洲国产欧美日韩在线播放| 高清欧美精品videossex| 自拍欧美九色日韩亚洲蝌蚪91| 狠狠婷婷综合久久久久久88av| tube8黄色片| 女人高潮潮喷娇喘18禁视频| 色播在线永久视频| 色婷婷av一区二区三区视频| 寂寞人妻少妇视频99o| 亚洲成人av在线免费| 欧美人与善性xxx| 大香蕉久久网| 丝袜美腿诱惑在线| 午夜日本视频在线| 制服诱惑二区| 色婷婷av一区二区三区视频| 久久久久人妻精品一区果冻| 欧美日韩视频高清一区二区三区二| 久久综合国产亚洲精品| 国产成人免费无遮挡视频| 咕卡用的链子| 亚洲精品成人av观看孕妇| 久久人人爽人人片av| 国产成人a∨麻豆精品| 成年女人毛片免费观看观看9 | 亚洲激情五月婷婷啪啪| 最近中文字幕2019免费版| 美女高潮到喷水免费观看| 亚洲,欧美,日韩| 亚洲精品一二三| 午夜免费男女啪啪视频观看| 一级黄片播放器| 777久久人妻少妇嫩草av网站| 日本欧美国产在线视频| 免费黄色在线免费观看| 十八禁高潮呻吟视频| 人妻一区二区av| 日韩人妻精品一区2区三区| 高清av免费在线| 美女国产高潮福利片在线看| 日本欧美视频一区| av.在线天堂| 桃花免费在线播放| 亚洲国产最新在线播放| 久久人人97超碰香蕉20202| 深夜精品福利| 视频在线观看一区二区三区| 在线观看三级黄色| 黄色毛片三级朝国网站| 成人二区视频| 午夜福利乱码中文字幕| 国产精品.久久久| 男女啪啪激烈高潮av片| 久久这里只有精品19| 有码 亚洲区| 国产激情久久老熟女| 成人18禁高潮啪啪吃奶动态图| 一本—道久久a久久精品蜜桃钙片| 色婷婷av一区二区三区视频| 亚洲内射少妇av| 在线观看免费视频网站a站| 日韩 亚洲 欧美在线| 一级毛片我不卡| 国产成人a∨麻豆精品| 91成人精品电影| 久久人人爽人人片av| 在线亚洲精品国产二区图片欧美| 国产日韩欧美亚洲二区| 亚洲精品美女久久av网站| 免费高清在线观看视频在线观看| 女人久久www免费人成看片| 欧美人与善性xxx| 国产精品麻豆人妻色哟哟久久| 国产成人欧美| 午夜影院在线不卡| 精品一品国产午夜福利视频| 日本黄色日本黄色录像| av网站免费在线观看视频| 精品国产乱码久久久久久小说| 亚洲欧洲日产国产| 国产成人免费无遮挡视频| 少妇熟女欧美另类| 少妇精品久久久久久久| 亚洲成av片中文字幕在线观看 | 黄片小视频在线播放| 黑丝袜美女国产一区| 最新的欧美精品一区二区| 日本色播在线视频| 七月丁香在线播放| 26uuu在线亚洲综合色| 久久精品久久久久久久性| 国产一级毛片在线| 亚洲国产欧美日韩在线播放| 欧美亚洲日本最大视频资源| 999精品在线视频| 国产在线一区二区三区精| 欧美av亚洲av综合av国产av | 国产精品一二三区在线看| 国产在线一区二区三区精| 久久久国产欧美日韩av| 成年动漫av网址| 伦理电影大哥的女人| 最近最新中文字幕大全免费视频 | 久久热在线av| 日本色播在线视频| 母亲3免费完整高清在线观看 | 亚洲国产av影院在线观看| 国产精品三级大全| 99九九在线精品视频| 欧美成人午夜精品| 国产亚洲精品第一综合不卡| 中文字幕人妻丝袜一区二区 | 亚洲精品久久成人aⅴ小说| 1024香蕉在线观看| 精品酒店卫生间| 日韩av不卡免费在线播放| 日本wwww免费看| 多毛熟女@视频| 国语对白做爰xxxⅹ性视频网站| 黄网站色视频无遮挡免费观看| 蜜桃国产av成人99| 亚洲国产精品国产精品| 尾随美女入室| 高清黄色对白视频在线免费看| 亚洲国产av影院在线观看| 青春草亚洲视频在线观看| 精品久久久精品久久久| 999久久久国产精品视频| 亚洲精品久久成人aⅴ小说| 国产成人精品福利久久| av在线观看视频网站免费| 99热网站在线观看| 久久人人爽av亚洲精品天堂| 国产成人aa在线观看| 日日啪夜夜爽| 久久久久精品人妻al黑| 两个人免费观看高清视频| 大香蕉久久网| 十八禁网站网址无遮挡| 亚洲精品一二三| 国产不卡av网站在线观看| 精品一品国产午夜福利视频| 亚洲成av片中文字幕在线观看 | 人妻一区二区av| 日韩在线高清观看一区二区三区| 久久久久久久久久久免费av| 成年美女黄网站色视频大全免费| 国产熟女午夜一区二区三区| 国产亚洲最大av| 2021少妇久久久久久久久久久| 中文字幕制服av| 欧美黄色片欧美黄色片| 男的添女的下面高潮视频| 成人国语在线视频| 午夜福利视频精品| 久久久久久免费高清国产稀缺| 考比视频在线观看| 男女下面插进去视频免费观看| 亚洲综合色惰| 校园人妻丝袜中文字幕| 精品卡一卡二卡四卡免费| 在线 av 中文字幕| 精品卡一卡二卡四卡免费| 国产片特级美女逼逼视频| 婷婷成人精品国产| 国产成人精品一,二区| 亚洲视频免费观看视频| 精品国产露脸久久av麻豆| 满18在线观看网站| 日本色播在线视频| 精品午夜福利在线看| 观看美女的网站| 99久久精品国产国产毛片| 超碰97精品在线观看| 日本欧美视频一区| 女性生殖器流出的白浆| freevideosex欧美| 亚洲成国产人片在线观看| 亚洲精品久久午夜乱码| 狂野欧美激情性bbbbbb| 自线自在国产av| 男人添女人高潮全过程视频| 久久鲁丝午夜福利片| 亚洲精品aⅴ在线观看| 久久久久久伊人网av| 日韩成人av中文字幕在线观看| 99九九在线精品视频| 熟妇人妻不卡中文字幕| a级片在线免费高清观看视频| 男女免费视频国产| 18禁裸乳无遮挡动漫免费视频| 99久久综合免费| 欧美日韩亚洲高清精品| 免费看不卡的av| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 久久国产亚洲av麻豆专区| 国产乱来视频区| 精品国产一区二区三区四区第35| 欧美激情极品国产一区二区三区| 日本av免费视频播放| 人人澡人人妻人| 亚洲三区欧美一区| 亚洲第一青青草原| 最近最新中文字幕大全免费视频 | av在线播放精品| 亚洲精品中文字幕在线视频| 在线精品无人区一区二区三| 欧美精品一区二区免费开放| 综合色丁香网| 啦啦啦啦在线视频资源| 高清视频免费观看一区二区| 黄片无遮挡物在线观看| 亚洲情色 制服丝袜| 侵犯人妻中文字幕一二三四区| 国产欧美亚洲国产| 午夜日本视频在线| 黑人猛操日本美女一级片| 侵犯人妻中文字幕一二三四区| 国产精品久久久久久精品电影小说| 午夜91福利影院| 国产片特级美女逼逼视频| 国产精品一二三区在线看| 国产在视频线精品| 狂野欧美激情性bbbbbb| 我要看黄色一级片免费的| 韩国av在线不卡| 男人舔女人的私密视频| 欧美在线黄色| 亚洲第一青青草原| av在线观看视频网站免费| 久久综合国产亚洲精品| 99热全是精品| 日韩中文字幕欧美一区二区 | 亚洲精品自拍成人| 一区二区三区四区激情视频| 男女午夜视频在线观看| 亚洲人成电影观看| 在现免费观看毛片| kizo精华| 国产国语露脸激情在线看| 春色校园在线视频观看| 国产又色又爽无遮挡免| 久久99蜜桃精品久久| 日本爱情动作片www.在线观看| 免费看av在线观看网站| 日韩熟女老妇一区二区性免费视频| 亚洲色图 男人天堂 中文字幕| 91久久精品国产一区二区三区| 下体分泌物呈黄色| 国产xxxxx性猛交| 狠狠婷婷综合久久久久久88av| 日韩制服骚丝袜av| 日日爽夜夜爽网站| 国产欧美日韩一区二区三区在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av欧美aⅴ国产| 免费不卡的大黄色大毛片视频在线观看| 99久久综合免费| 久久久久久久国产电影| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区二区三区精品91| 国产视频首页在线观看| 久久精品夜色国产| 日韩一区二区视频免费看| 国产成人精品久久久久久| 一级毛片 在线播放| 亚洲综合精品二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | √禁漫天堂资源中文www| 999精品在线视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品乱久久久久久| 亚洲精品第二区| 久久久久国产网址| 日韩中文字幕视频在线看片| 侵犯人妻中文字幕一二三四区| 亚洲,欧美,日韩| 亚洲美女视频黄频| 精品久久久精品久久久| 考比视频在线观看| 欧美日韩视频精品一区| 9色porny在线观看| 麻豆乱淫一区二区| 久久狼人影院| 国产白丝娇喘喷水9色精品| 丝袜美腿诱惑在线| 亚洲av.av天堂| 99九九在线精品视频| 亚洲精品成人av观看孕妇| 新久久久久国产一级毛片| 日韩在线高清观看一区二区三区| 欧美激情 高清一区二区三区| 亚洲在久久综合| 国产有黄有色有爽视频| 91午夜精品亚洲一区二区三区| 亚洲美女黄色视频免费看| 最新中文字幕久久久久| 免费不卡的大黄色大毛片视频在线观看| 2018国产大陆天天弄谢| 久久热在线av| 亚洲成人手机| a 毛片基地| 欧美成人精品欧美一级黄| 涩涩av久久男人的天堂| 久热这里只有精品99| 十八禁网站网址无遮挡| 最近中文字幕2019免费版| 精品人妻偷拍中文字幕| 可以免费在线观看a视频的电影网站 | 少妇被粗大的猛进出69影院| 男女啪啪激烈高潮av片| 国产日韩一区二区三区精品不卡| 天堂俺去俺来也www色官网| 亚洲,欧美,日韩| 少妇被粗大猛烈的视频| 午夜福利视频精品| 涩涩av久久男人的天堂| 亚洲欧美一区二区三区黑人 | 久久人人97超碰香蕉20202| 香蕉丝袜av| 爱豆传媒免费全集在线观看| 亚洲精品国产av成人精品| 一区福利在线观看| 91aial.com中文字幕在线观看| 国产免费视频播放在线视频| 在线观看人妻少妇| 亚洲图色成人| 成人毛片a级毛片在线播放| 国产日韩一区二区三区精品不卡| 精品国产超薄肉色丝袜足j| 老汉色av国产亚洲站长工具| 国产国语露脸激情在线看| 欧美日韩国产mv在线观看视频| 波多野结衣av一区二区av| 亚洲激情五月婷婷啪啪| 高清欧美精品videossex| 久久婷婷青草| 亚洲精品视频女| 日本爱情动作片www.在线观看| av线在线观看网站| 亚洲欧美一区二区三区黑人 | 久久久精品区二区三区| 九草在线视频观看| 色吧在线观看| a 毛片基地| 国产精品二区激情视频| 又大又黄又爽视频免费| 日韩免费高清中文字幕av| 久热这里只有精品99| 成人手机av| av卡一久久| 免费黄频网站在线观看国产| 亚洲人成77777在线视频| 国产黄频视频在线观看| 亚洲国产欧美在线一区| 天堂8中文在线网| 国产在视频线精品| 美女视频免费永久观看网站| 人妻一区二区av| 2022亚洲国产成人精品| √禁漫天堂资源中文www| 日本欧美国产在线视频| 极品少妇高潮喷水抽搐| 国产精品一区二区在线观看99| 国产精品一国产av| 黑人猛操日本美女一级片| 精品国产国语对白av| 亚洲精品国产av蜜桃| 欧美人与善性xxx| 日本色播在线视频| 黄色毛片三级朝国网站| 999久久久国产精品视频| www.av在线官网国产| av福利片在线| 91精品国产国语对白视频| 国产欧美日韩一区二区三区在线| 亚洲第一青青草原| 成人国产麻豆网| 亚洲欧美一区二区三区国产| 亚洲成色77777| 久久热在线av| 欧美在线黄色| av免费观看日本| 亚洲欧美中文字幕日韩二区| 90打野战视频偷拍视频| 国产综合精华液| 亚洲精品久久成人aⅴ小说| 久久久久久伊人网av| 午夜精品国产一区二区电影| 久久免费观看电影| 婷婷色综合www| 亚洲男人天堂网一区| 日韩中文字幕视频在线看片| 免费看av在线观看网站| 男女免费视频国产| 高清不卡的av网站| 亚洲人成网站在线观看播放| 色婷婷av一区二区三区视频| 丰满饥渴人妻一区二区三| 高清黄色对白视频在线免费看| 看免费av毛片| 精品一区二区三卡| 免费播放大片免费观看视频在线观看| 久久精品夜色国产| 亚洲国产精品一区二区三区在线| 毛片一级片免费看久久久久| 深夜精品福利| 国产片内射在线| 久久国产精品男人的天堂亚洲| 又大又黄又爽视频免费| 欧美精品亚洲一区二区| 黑丝袜美女国产一区| 久久97久久精品| 欧美黄色片欧美黄色片| 2018国产大陆天天弄谢| 午夜精品国产一区二区电影| 亚洲精品日韩在线中文字幕| 一本大道久久a久久精品| 黄色一级大片看看| 亚洲欧洲精品一区二区精品久久久 | 亚洲精华国产精华液的使用体验| 各种免费的搞黄视频| 在线亚洲精品国产二区图片欧美| 日韩av在线免费看完整版不卡| 一区二区三区乱码不卡18| 18禁国产床啪视频网站| 一本色道久久久久久精品综合| 丝袜在线中文字幕| 亚洲一区二区三区欧美精品| 国产精品熟女久久久久浪| 日韩熟女老妇一区二区性免费视频| 制服人妻中文乱码| 日韩精品免费视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 国产成人欧美| 亚洲精品久久久久久婷婷小说| 国产精品嫩草影院av在线观看| 国产精品一区二区在线观看99| 日本av手机在线免费观看| 国产又色又爽无遮挡免| 午夜激情久久久久久久| 亚洲精品一区蜜桃| 成人黄色视频免费在线看| 午夜福利视频精品| 亚洲精品美女久久久久99蜜臀 | 国产野战对白在线观看| 久久精品久久久久久噜噜老黄| 黄网站色视频无遮挡免费观看| 久久久久久久久久久久大奶| 精品久久久精品久久久| 久久精品人人爽人人爽视色| 人人妻人人添人人爽欧美一区卜| av视频免费观看在线观看| 国产男人的电影天堂91| 在线观看人妻少妇| 亚洲av男天堂| 乱人伦中国视频| 精品国产一区二区久久| 一级,二级,三级黄色视频| 国产男女内射视频| 亚洲精品一区蜜桃| 精品一区二区三区四区五区乱码 | 丰满饥渴人妻一区二区三| 美女大奶头黄色视频| 99久国产av精品国产电影| 免费观看av网站的网址| av不卡在线播放| 亚洲欧美精品自产自拍| 日本欧美国产在线视频| 亚洲欧美一区二区三区黑人 | 久久久久国产网址| 啦啦啦中文免费视频观看日本| 久久亚洲国产成人精品v| 2018国产大陆天天弄谢| 久久国产精品男人的天堂亚洲| 亚洲一级一片aⅴ在线观看| 国产日韩欧美亚洲二区| 99热网站在线观看| av在线app专区| 国产成人91sexporn| 高清在线视频一区二区三区| 精品一区二区三区四区五区乱码 | 另类亚洲欧美激情| 日本猛色少妇xxxxx猛交久久| 男人添女人高潮全过程视频| 少妇的丰满在线观看| 在线观看一区二区三区激情| 可以免费在线观看a视频的电影网站 | 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利视频在线观看免费|