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

    糜子分子遺傳研究進(jìn)展與展望

    2022-04-29 00:00:00楊天育
    甘肅農(nóng)業(yè)科技 2022年10期

    摘要:糜子是最古老的馴化栽培作物之一,具有抗旱、耐瘠、水分利用效率高、營養(yǎng)豐富等特點,在北方旱作農(nóng)業(yè)區(qū)特色糧食生產(chǎn)中占有重要地位。為深入推動糜子農(nóng)藝、產(chǎn)量及品質(zhì)相關(guān)性狀遺傳解析,在介紹國內(nèi)糜子遺傳研究現(xiàn)狀基礎(chǔ)上,通過文獻(xiàn)分析綜述了糜子分子遺傳研究進(jìn)展,并展望了未來發(fā)展趨勢。

    關(guān)鍵詞:糜子;基因組;遺傳圖譜;BSA;轉(zhuǎn)錄組

    中圖分類號:S516" " " " " " " 文獻(xiàn)標(biāo)志碼:A" " " " " " " 文章編號:2097-2172(2022)01-0032-05

    doi:10.3969/j.issn.2097-2172.2022.01.006

    Research Progress and Prospect on Molecular Genetics

    of Millet(Panicum miliaceum L.)

    YANG Tianyu

    (Gansu Academy of Agricultural Sciences, Lanzhou Gansu 730070, China)

    Abstract: Broom corn millet (Panicum miliaceum L.) is one of the most ancient crops domestically cultivated in the world, it possesses characters of drought resistance, barren-tolerance, high efficiency in water utilization and dense nutrients which plays an important role in the characteristic grain productionin northern dryland zone. To further promote the genetic analysis of correlated characters in agronomy traits, yields and nutrients of broom corn millet, based on thedomestic research status of broom corn millet genetics, advances in molecular genetics of broom corn milletwerereviewed and prospect was made in this paper.

    Key words: Panicum miliaceum L.; Genome; Genetic map; BSA; Transcriptome

    糜子(Panicum miliaceum L.)是禾本科黍?qū)俚囊粋€異源四倍體種[1 ],是最古老的馴化栽培作物之一,主要分布在歐亞大陸干旱半干旱地區(qū)[2 ]??脊抛C明,糜子在我國有超過10 000年的栽培歷史[3 ],廣泛地分布在我國各個地區(qū),從北緯19° 15′(海南)到49° 18′(新疆),東經(jīng)76° (新疆)到143° (黑龍江)及海拔200 m(山東)到3 000 m(西藏)都有種植[4 ]。據(jù)國家谷子高粱產(chǎn)業(yè)技術(shù)體系統(tǒng)計,我國糜子播種面積53.3萬多hm2,主要分布在山西、陜西、內(nèi)蒙古、甘肅、河北、寧夏等6個省區(qū),平均量產(chǎn)1 500~2 250 kg/hm2。由于長期的馴化栽培和廣泛的地理分布,形成了我國豐富的糜子種質(zhì)資源,目前我國庫存糜子種質(zhì)資源8 515份[5 ],保存量世界第1位。糜子具有生育期短、水分利用效率高、耐鹽堿、養(yǎng)分利用率高、C4光合屬性及營養(yǎng)豐富等特點,是輪作填閑、邊緣瘠薄地種植的先鋒作物[6 ],也是我國北方傳統(tǒng)的制米作物,在北方旱作農(nóng)業(yè)生產(chǎn)中占有重要地位。

    糜子基礎(chǔ)研究與產(chǎn)業(yè)發(fā)展水平與玉米、水稻、小麥等大宗作物相比差距很大,究其原因是糜子基礎(chǔ)研究相對落后,尤其是產(chǎn)量、品質(zhì)與抗性等遺傳基礎(chǔ)研究不透不深,無法在關(guān)鍵核心技術(shù)上有效支撐產(chǎn)業(yè)發(fā)展。糜子產(chǎn)量和品質(zhì)等遺傳研究大多集中品種資源遺傳差異與多樣性等方面[7 - 11 ],關(guān)于質(zhì)量和數(shù)量性狀位點/基因等分子遺傳的深入研究很少。因此,總結(jié)糜子分子遺傳學(xué)研究進(jìn)展并進(jìn)行科學(xué)研究問題展望,有助于深入解析糜子農(nóng)藝、產(chǎn)量、品質(zhì)及抗逆性等相關(guān)性狀的遺傳調(diào)控機理,推動糜子應(yīng)用基礎(chǔ)研究和產(chǎn)業(yè)技術(shù)的發(fā)展。

    1" "國內(nèi)糜子遺傳研究現(xiàn)狀

    據(jù)國家谷子高粱產(chǎn)業(yè)技術(shù)體系科技發(fā)展報告顯示,近年來在國家產(chǎn)業(yè)技術(shù)體系支持下,國內(nèi)從事糜子研究的單位開展了糜子種質(zhì)資源分子水平遺傳多樣性研究,全面評價了我國現(xiàn)有糜子地方資源、育成品種和野生材料之間的遺傳差異和群體遺傳結(jié)構(gòu);采用SSR、ITS、簡化基因組測序等手段對野生糜子分類地位、栽培糜子的起源、馴化與傳播進(jìn)行了深入研究;開展糜子全基因組測序工作和糜子核心種質(zhì)的關(guān)聯(lián)分析,基本明確了糜子群體遺傳結(jié)構(gòu),構(gòu)建了第一張糜子SSR標(biāo)記的遺傳圖譜,研究了野生糜子起源進(jìn)化的路徑和光周期轉(zhuǎn)錄組基因的差異,初步定位了控制糜子雄性不育和株高性狀的基因,基因組、轉(zhuǎn)錄組學(xué)的糜子基礎(chǔ)研究為促進(jìn)分子育種技術(shù)開發(fā)奠定了基礎(chǔ),也為認(rèn)識糜子從野生種到栽培種的進(jìn)化提供了科學(xué)依據(jù)。

    糜子是禾谷類作物中最耐旱的作物,具有C4植物光合途徑的3種亞類型[12 ],不僅具有高的光合效率,還具有水分和養(yǎng)分利用效率高等特點。糜子曾經(jīng)是我國北方重要栽培作物,但隨著農(nóng)業(yè)生產(chǎn)發(fā)展逐漸退出主栽作物地位成為區(qū)域特色作物,根本的原因是糜子產(chǎn)量沒有大的突破,品質(zhì)沒有根本提升,輕簡栽培技術(shù)無法快速推廣,制約了糜子產(chǎn)業(yè)的發(fā)展。這其中內(nèi)在的原因,是糜子基礎(chǔ)研究薄弱,對產(chǎn)量、品質(zhì)和抗性的研究不深,重要性狀功能基因挖掘不夠,遺傳調(diào)控機理不清,基礎(chǔ)研究創(chuàng)新成果的缺乏無法支撐技術(shù)的重大進(jìn)步。因此,了解糜子分子遺傳研究進(jìn)展對深入解析糜子農(nóng)藝、產(chǎn)量、品質(zhì)及抗逆性等相關(guān)性狀的遺傳調(diào)控機制具有重要意義。

    2" "糜子分子遺傳研究進(jìn)展

    2.1" "糜子基因組學(xué)的研究

    2019年,糜子2個參考基因組的釋放為產(chǎn)量、品質(zhì)和抗逆性等功能基因挖掘提供了重要研究平臺,對于解析糜子重要性狀分子機制及育種應(yīng)用也產(chǎn)生了重要影響[12 - 13 ]。Shi等[13 ]以“隴糜4號”為材料,聯(lián)合PacBio測序、BioNano和Hi-C作圖進(jìn)行了糜子基因組組裝,18個超級支架覆蓋了95.6%基因組(約887.8 Mb),注釋了63 671個蛋白質(zhì)編碼基因。進(jìn)化分析顯示,約86.2%的谷子共線基因在糜子中有2個同源拷貝,這表明谷子的近源種糜子在四倍體化后基因很少丟失。系統(tǒng)發(fā)育分析表明,谷子和糜子大約在1 310萬a前分化,而糜子的異源四倍體可能在591萬年內(nèi)發(fā)生。Zou等[12 ]以1個糜子地方品種為材料(國家種質(zhì)資源庫編號00000390),聯(lián)合short-read測序、單分子實時測序、Hi-C和高密度遺傳圖譜來組裝了另一個含有55 930個蛋白質(zhì)編碼基因和339個microRNA基因的高質(zhì)量、染色體水平糜子參考基因組。系統(tǒng)發(fā)育分析顯示,糜子異源四倍體的兩組同源染色體的合并約發(fā)生在560萬a,這兩組染色體都與玉米、水稻、谷子、高粱等其他禾本科種具有很強的同源性。他們還發(fā)現(xiàn),泛素E3連接酶BTB亞基的幾個亞家族中存在泛素類特異性擴張,表明蛋白質(zhì)動力學(xué)的強化調(diào)控可能有助于糜子的進(jìn)化。此外,還明確碳固定候選基因中均存在3個C4亞型相關(guān)基因,這為研究特殊非生物脅迫耐性和C4生物學(xué)奠定了基礎(chǔ)。

    2.2" "糜子遺傳圖譜的構(gòu)建

    遺傳圖譜(Genetic map)或稱為遺傳連鎖圖譜(Genetic linkage map),是通過遺傳距離來標(biāo)識已知的遺傳標(biāo)記或基因在基因組或染色體的相對位置。構(gòu)建連鎖圖譜的基礎(chǔ)是以同源染色體之間發(fā)生重組和交換為基礎(chǔ),染色體片段/基因間隔越大,發(fā)生重組和交換的概率就越高[14 ]。通常以重組率為遺傳距離,即1個厘摩(cM)大小相當(dāng)于百分之一的連鎖重組率[15 ] 。構(gòu)建遺傳圖譜包括4個步驟:選擇遺傳標(biāo)記、建立作圖群體、分析群體間不同品系或植株的標(biāo)記基因型、確定標(biāo)記位點與標(biāo)記間連鎖群[16 ]。隨著分子標(biāo)記的發(fā)展,連鎖圖譜構(gòu)建經(jīng)歷了以形態(tài)標(biāo)記、同工酶標(biāo)記、細(xì)胞學(xué)標(biāo)記、DNA標(biāo)記為主的研究過程,目前隨著新一代測序技術(shù)的發(fā)展,以SNP為主要DNA標(biāo)記的連鎖圖得到了普遍應(yīng)用。在遺傳研究中,連鎖圖譜的構(gòu)建在比較基因組作圖、基因定位與克隆、分子標(biāo)記輔助育種等方面有著非常重要的應(yīng)用價值[17 ]。2015年,王銀月等[18 - 19 ]采用高通量測序與磁珠富集法相結(jié)合開發(fā)出1210對SSR引物并在4個品種(會寧大黃糜、隴糜8號、太原55號和會寧野糜子)間篩選出116對多態(tài)性較好的引物用于糜子遺傳圖譜的構(gòu)建。翌年,從3506對引物中篩選出235對多態(tài)性引物,構(gòu)建了兩個F2群體(♀會寧大黃糜×♂隴糜8號、♀太原55號×♂會寧野糜子),分別繪制出覆蓋660.3 cM、251.3 cM,含13、3個連鎖群的遺傳圖譜。2016年,Rajput等[20 ]以Huntsman與Minsum 2個糜子品種雜交構(gòu)建的包含93個家系的群體為材料,共用833 個GBS-SNP標(biāo)記進(jìn)行了連鎖圖譜構(gòu)建,結(jié)果519個標(biāo)記未形成連鎖群,共有117個標(biāo)記分布在18個主連鎖群上,連鎖群基因組長度2 137 cM,標(biāo)記之間平均距離為18 cM。18個主連鎖群的標(biāo)記長度54.6~236 cM,數(shù)量為4~12。在14個連鎖群上共檢測到8個形態(tài)農(nóng)藝性狀的18個QTL,每個QTL都可以解釋13.2%~34.7%表型變異。2019年,Zou等[12 ]為驗證組裝基因組的準(zhǔn)確性,以一個亞洲品種和一個北美品種雜交構(gòu)建的含132個家系的RIL群體為材料,通過全基因組測序,獲得了221 787個高質(zhì)量SNP標(biāo)記并構(gòu)建了覆蓋18條染色體的連鎖圖譜,其中母本總遺傳長度為" 2 811 cM,父本總遺傳長度為3 092 cM。上述研究工作為糜子分子標(biāo)記開發(fā)奠定了良好基礎(chǔ),但目標(biāo)性狀標(biāo)記實用性不夠,且此后糜子遺傳連鎖圖譜構(gòu)建及其相關(guān)研究很少報道,說明尚有待加強糜子分子遺傳的研究。

    2.3" "糜子集群分離分析法的應(yīng)用

    集群分離分析法又稱為混合分組分析法(BSA),是一項快速挖掘QTL的簡易方法,主要用于快速識別到影響目的性狀的遺傳位點。Michelmore等[21 ]首次提出BSA法并成功應(yīng)用在萵苣中,篩選到3 個與目標(biāo)性狀相連鎖的標(biāo)記。與傳統(tǒng)的QTL定位方法相比較,BSA只需根據(jù)目標(biāo)性狀挑選出的極端差異個體來進(jìn)行測序就能夠得到目的性狀的關(guān)聯(lián)基因組候選區(qū)間,很大程度上降低了研究成本,因而廣泛應(yīng)用在水稻、玉米、西瓜、油菜、甜瓜等作物中[22 - 26 ],具有高效、快速、節(jié)省成本等特點[21 ]。

    近年來,隨著測序成本的降低和測序技術(shù)不斷發(fā)展,BSA-seq技術(shù)成為QTL定位更快捷和有效的工具[27 ],其在基因定位方面的作用越來越突出。Venuprasad等[28 ]在水稻重組自交系群體中使用BSA技術(shù)檢測出在干旱脅迫下與產(chǎn)量密切相關(guān)的2個SSR標(biāo)記,第3號染色體上與RM 416連鎖的QTL解釋了31%的產(chǎn)量性狀遺傳變異率,成為水稻中第一個在有氧環(huán)境和干旱中對產(chǎn)量影響較大的QTL。李麗[16 ]使用BSA 技術(shù)將與花生株型性狀相關(guān)的12 個QTLs定位到15號染色體上。李玉穎等[29 ]通過BSA技術(shù)在第1染色體上定位到1個與花生籽仁含油量相關(guān)的候選區(qū)域,其中的1個候選基因Arahy.55ECQ6,檢測到50個注釋基因和309個SNP。趙鈺涵等[30 ]通過SNP芯片與BSA 相結(jié)合,將控制花生紫色種皮的基因定位到第10號染色體上,并篩選出一個與花生黑種皮性狀緊密連鎖的SSR標(biāo)記。陳鵬云[31 ]將BSA-seq法和構(gòu)建連鎖圖譜的方法結(jié)合對棉花的早花性狀進(jìn)行QTL定位,最后將候選區(qū)間定位在第17染色體的40.33~41.40 Mb和15.96~25.68 Mb區(qū)段內(nèi),同源基因進(jìn)化分析和基因結(jié)構(gòu)分析確定了與開花相關(guān)的3個候選基因。在糜子中,目前僅有甘肅省農(nóng)業(yè)科學(xué)院小雜糧課題組利用高稈品種“隴糜12號”和EMS誘變矮稈系“張778”雜交構(gòu)建了的包含939個家系的F2群體,采用BSA-seq在第1染色體分析鑒定到調(diào)控糜子株高的關(guān)聯(lián)位點,進(jìn)一步開發(fā)InDel標(biāo)記構(gòu)建了連鎖圖譜,結(jié)合F2、F2∶3 939個家系株高表型數(shù)據(jù),將該QTL區(qū)間縮小到109 kb,QTL分析發(fā)現(xiàn)該位點是調(diào)控糜子株高的一個主效位點,加性效應(yīng)達(dá)15 cm,同時鑒定到2個CDS區(qū)發(fā)生單堿基非同義突變和1個啟動子區(qū)發(fā)生堿基插入/缺失的候選基因,并利用該區(qū)間內(nèi)的InDel標(biāo)記關(guān)聯(lián)分析了512份糜子種質(zhì)資源,驗證了該QTL是一個調(diào)控糜子株高的關(guān)鍵位點,為糜子株高育種及遺傳改良提供了分子遺傳基礎(chǔ)[32 ],同時,甘肅省農(nóng)業(yè)科學(xué)院小雜糧課題組基于該群體,利用BSA-seq解析了調(diào)控糜子黃、褐粒色的候選基因,進(jìn)一步驗證了利用BSA-seq結(jié)合QTL定位,能快速地解析個體間存在極端差異的性狀的調(diào)控基因(待發(fā)表)。

    2.4" "糜子轉(zhuǎn)錄組學(xué)的研究

    轉(zhuǎn)錄是指將遺傳信息從DNA傳遞到 RNA的過程,而轉(zhuǎn)錄組是指生物體在某一階段下,細(xì)胞內(nèi)所有進(jìn)行轉(zhuǎn)錄的總 RNA,包括非編碼RNA 和編碼 RNA。轉(zhuǎn)錄組由Vekulescu等[33 ]提出,不同于基因組,轉(zhuǎn)錄組并不是相對固定的,基因在表達(dá)時會受到生長環(huán)境、生育時期等條件的影響,從而具有明顯的空間和時間性[34 ]。轉(zhuǎn)錄組學(xué)是一個從整體水平上來研究細(xì)胞中所有基因的基因結(jié)構(gòu)、功能、轉(zhuǎn)錄以及調(diào)控規(guī)律的一門學(xué)科[35 ]。近年來,隨著測序技術(shù)的發(fā)展,轉(zhuǎn)錄組學(xué)被廣泛應(yīng)用到不同作物多種環(huán)境下各類組織基因的表達(dá)研究中。Zhang等[36 ]選用干旱敏感材料晉黍6號(JS6)和PEG誘導(dǎo)水分脅迫下的耐旱材料內(nèi)糜5號(NM5),轉(zhuǎn)錄組測序后發(fā)現(xiàn),在無PEG誘導(dǎo)水分脅迫條件下,JS6和NM5中觀察到1 695個差異表達(dá)基因(DEGs);2個品種分別使用20% PEG-6000處理6 h和24 h,在模擬干旱條件下,分別檢測到833、2 166個差異表達(dá)基因。在PEG-6000處理6 h和24 h下,晉黍6號比內(nèi)糜5號差異表達(dá)基因分別高0.298和0.754倍。此外,在NM5中觀察到ROS清除系統(tǒng)對模擬干旱處理的轉(zhuǎn)錄反應(yīng)延遲,而光合作用相關(guān)基因的表達(dá)相對容易恢復(fù);NM5的茉莉酸(JA)信號轉(zhuǎn)導(dǎo)途徑與JS6相比也有不同調(diào)控策略。Shan等[37 ]通過轉(zhuǎn)錄組分析,發(fā)現(xiàn)NAC基因家族在不同干旱脅迫時間內(nèi),表達(dá)量在不同組織中存在極大差異。甘肅省農(nóng)業(yè)科學(xué)院小雜糧課題研究發(fā)現(xiàn),在PEG-6000模擬干旱脅迫下,轉(zhuǎn)錄組測序結(jié)果YABBY及bZIP基因家族成員均在糜子幼苗期表達(dá)量存在顯著差異[38 - 39 ],說明了基因家族普遍參與了糜子的抗旱調(diào)控。糜子轉(zhuǎn)錄組測序技術(shù)已經(jīng)應(yīng)用到其抗旱性相關(guān)基因挖掘,而利用轉(zhuǎn)錄組挖掘糜子其他特性方面的相關(guān)基因或分子機制的研究未見報道,反映出糜子轉(zhuǎn)錄組學(xué)方面的研究仍相對缺乏。

    3" "糜子分子遺傳研究發(fā)展趨勢

    世界生物育種技術(shù)發(fā)展已經(jīng)歷了原始馴化選育、常規(guī)育種、分子育種時代,正向設(shè)計育種或智能化育種時代發(fā)展,而智能化育種將基因編輯、生物育種及人工智能相互融合,將實現(xiàn)性狀的精準(zhǔn)定向改良[40 ]。我國糜子育種目前仍處在雜交育種、誘變育種為主要育種方法的常規(guī)育種階段,糜子雜種優(yōu)勢利用剛剛起步,分子標(biāo)記輔助選擇育種僅是個例[38 ],轉(zhuǎn)基因育種相關(guān)技術(shù)的報道更少。未來作物育種呼吁更高的技術(shù), 因此糜子育種同樣需要緊跟時代要求,加強常規(guī)育種技術(shù)與現(xiàn)代生物育種技術(shù)的結(jié)合,綜合應(yīng)用功能基因組學(xué)、轉(zhuǎn)錄組學(xué)、全基因組關(guān)聯(lián)分析、基因工程、染色體工程等方法,推進(jìn)糜子育種技術(shù)不斷革新進(jìn)步,從而培育高產(chǎn)優(yōu)質(zhì)多抗新品種滿足市場的不同需求。

    糜子具有節(jié)水耐旱耐瘠薄的優(yōu)良特性,盡管相關(guān)研究已經(jīng)報道了涉及糜子水分利用、耐旱性及土壤養(yǎng)分利用效率的相關(guān)基因,但研究深度遠(yuǎn)遠(yuǎn)滿足不了分子育種的基本需求,仍是糜子中亟待研究的兩個重要方面。隨著我國農(nóng)業(yè)生產(chǎn)組織方式逐漸向種植大戶、專業(yè)合作社發(fā)展,開展適應(yīng)機械化管理和收獲的作物株型研究成為熱點[41 ],糜子因生育期水肥條件充足生產(chǎn)上極易出現(xiàn)倒伏現(xiàn)象,影響機械化收獲,因此要實現(xiàn)抗倒伏、宜機收的糜子品種選育的突破,抗倒伏相關(guān)基因及分子調(diào)控網(wǎng)絡(luò)的解析、株型研究及株型育種同樣是未來糜子分子遺傳研究的一個重要方向。糜子是C4作物,在C4反應(yīng)過程中同時兼具3種脫羧途徑,因此光合效率更高。近年來,糜子近緣種谷子已發(fā)展成為解析禾谷類作物C4光合作用的模式作物[42 ],因此,糜子C4光合作用研究也將是未來一個重要研究方向。糜子營養(yǎng)價值很高,微量元素Fe及必需氨基酸含量高于其他禾谷類作物[43 ],糜子和其他小宗作物一樣受市場歡迎,正是滿足了人們吃得好、吃得健康而對谷物營養(yǎng)品質(zhì)的需要。因此,解析糜子品質(zhì)性狀的分子遺傳機制,是現(xiàn)代育種必須關(guān)注的一個方向。

    參考文獻(xiàn):

    [1] HUNT HV, BADAKSHI F, ROMANOVA O, et al. Reticulate evolution in Panicum(Poaceae): the origin of tetraploid broomcorn millet, P. miliaceum[J]." Journal of Experimental Botany, 2014, 65: 3165-3175.

    [2] FULLER D Q. Agricultural origins and frontiersin South Asia: A working synthesis[J]." Journal of World Prehistory, 2006, 20: 1-86.

    [3] YANG X, WAN Z, PERRY L, et al. Early millet use in northern China[J]." PNAS, 2012, 109(10): 3726-3730.

    [4] WANG R, HUNT HV, QIAO Z, et al." Diversity and Cultivation of Broomcorn Millet(Panicum miliaceum L.) in China: A Review[J]." Economic Botany, 2016, 70(3): 332-342.

    [5] 王" "綸,王星玉,溫琪汾,等." 中國黍稷種質(zhì)資源研究與利用[J]." 植物遺傳資源學(xué)報,2005,6(4):474-477.

    [6] YUE H, WANG M, LIU S, et al." Transcriptome-wide identification and expression profiles of the WRKY transcription factor family in Broomcorn millet(Panicum miliaceum L.)[J]." BMC Genomics, 2016, 17: 343.

    [7] HABIYAREMYE C, MATANGUIHAN J B, GUEDES J D, et al. Proso Millet(Panicum miliaceum L.) and Its Potentialfor Cultivation in the Pacific North west, U. S.: A Review[J]." Frontiers in Plant Science, 2016, 7: 1961.

    [8] GHIMIRE B K, YU C Y, KIM S H, et al. Diversity in Accessions of Panicum miliaceum L. Based on Agro-Morphological, Antioxidative, and Genetic Traits[J]." Molecules, 2019, 24: 1012.

    [9] AMADOU I, GOUNGA M E, LE G W." Millets: Nutritional composition, some health benefits and processing-A Review[J]." Emirates Journal of Food and Agriculture, 2013, 25(7): 501-508.

    [10] VETRIVENTHAN M, UPADHYAYA HD." Diversity and trait-specific sources for productivity and nutritional traits in the global proso millet(Panicum miliaceum L.) germplasm collection[J]." The Crop Journal, 2018, 6: 451-463.

    [11] MANIMOZHI S V, NIRMALAKUMARI A, SENTHIL N. Genetic Diversity for Zinc, Calcium and Iron Content of Selected Little Millet Genotypes[J]." Nutrition and Food Sciences, 2015, 5:6.

    [12] ZOU C, LI L, MIKI D, et al. The genome of broomcorn millet[J]." Nature communications, 2019, 10: 436.

    [13] SHI J, MA X, ZHANG J, et al." Chromosome conformation capture resolved nearcomplete genome assembly of broomcorn millet[J]." Nature communications, 2019, 10: 464.

    [14] 張國棟." 不同磷環(huán)境下玉米行粒數(shù)的精細(xì)定位[D]. 泰安:山東農(nóng)業(yè)大學(xué),2012.

    [15] 趙文明." 玉米株型相關(guān)性狀QTL 定位與分析[D]." 鄭州:河南農(nóng)業(yè)大學(xué),2008.

    [16] 李" "麗." 利用連鎖和全基因組關(guān)聯(lián)分析鑒定花生株型相關(guān)性狀的QTLs[D]." 保定:河北農(nóng)業(yè)大學(xué),2019.

    [17] 王偉繼,岳志芹,孔" "杰,等." AFLP分子標(biāo)記技術(shù)的發(fā)展及其在海洋生物中的應(yīng)用[J]." 海洋水產(chǎn)研究,2005,26(1):80-85.

    [18] 王銀月,劉敏軒,陸" "平,等." 構(gòu)建黍稷分子遺傳圖譜SSR引物的篩選[J]." 作物雜志,2014(4):32-38.

    [19] 王銀月." 黍稷種質(zhì)遺傳資源多樣性及遺傳圖譜構(gòu)建研究[D]." 長春:吉林農(nóng)業(yè)大學(xué),2015.

    [20] RAJPUT S G, SANTRA D K, SCHNABLE J. Mapping QTLs for morpho-agronomic traits in proso millet(Panicum miliaceum L.)[J]." Molecular Breeding, 2016, 36:37.

    [21] MICHELMORE R W, KESSELI I P V." Identification of Markers Linked to Disease-Resistance Genes by Bulked Segregant Analysis: A Rapid Method to Detect Markers in Specific Genomic Regions by Using Segregating Populations[J]." Proceedings of the National Academy of Sciences of the United States of America, 1991, 88(21): 9828-9832.

    [22] 張云輝,張所兵,林" "靜,等." 利用BSA法檢測水稻條紋葉枯病高效應(yīng)抗性位點[J]." 華北農(nóng)學(xué)報,2014,29(2):85-88.

    [23] 呂香玲,鄭克志,李" "元,等." 利用BSA法發(fā)掘玉米抗灰斑病主效QTL[J]." 玉米科學(xué),2015,23(5):16-20.

    [24] DONG W, WU D F, LI G S, et al." Next-generation sequencing from bulked segregant analysis identifies a dwarfism gene in watermelon[J]." Scientific Reports, 2018, 8(1): 2908-2914.

    [25] 張堯鋒,張冬青,余華勝,等." 基于極端混合池(BSA)全基因組重測序的甘藍(lán)型油菜有限花序基因定位[J]." 中國農(nóng)業(yè)科學(xué),2018,51(16):3029-3039.

    [26] 王學(xué)征,邱" "果,陳克農(nóng),等. 基于BSA法開發(fā)CAPS標(biāo)記定位甜瓜果面溝相關(guān)基因研究[J]." 東北農(nóng)業(yè)大學(xué)學(xué)報,2018,49(5):20-26.

    [27] ZHENG Y, XU F, LI Q, et al." QTL mapping combined with bulked segregant analysis identify SNP markers linked to leaf shape traits in Pisum sativum using SLAF sequencing[J]." Frontiers in Genetics, 2018, 5(9): 615-625.

    [28] VENUPRASAD R, DALID C O, VALLE M D, et al." Identification and characterization of large-effect quantitative trait loci for grain yield under lowland drought stress in rice using bulk-segregant analysis[J]." Theoretical and Applied Genetics, 2009, 120(1): 177-190.

    [29] 李玉穎,于海洋,呂玉英,等." 基于BSA 重測序定位花生含油量相關(guān)基因位點[J]." 山東農(nóng)業(yè)科學(xué),2021,53(1):1-6.

    [30] 趙鈺涵,周希萌,趙慧玲,等. 利用SNP芯片結(jié)合 BSA 的方法定位花生黑色種皮顏色基因[J]." 分子植物育種,2021,19(9):2977-2984.

    [31] 陳鵬云." 棉花早熟性狀相關(guān)基因的定位與鑒定[D]. 長沙:湖南農(nóng)業(yè)大學(xué),2017.

    [32] LIU T, LIU X, HE J, et al. Identification and fine-mappingof a major QTL(PH1.1) conferringplant height in broomcorn millet(Panicum miliaceum)[J]." Frontiers in Plant Science, 2022, 13:1010057.

    [33] VELCULESCU V E, ZHANG L, ZHOU W, et al. Characterization of the yeast transcriptome[J]." Cell, 1997," 88: 243-251.

    [34] 周金龍. 玉米光周期敏感近等基因系的轉(zhuǎn)錄組分析[D]." 鄭州:河南農(nóng)業(yè)大學(xué),2016.

    [35] 高" "寧,劉" "博,楊德強,等." 轉(zhuǎn)錄組學(xué)在中藥研究中的應(yīng)用現(xiàn)狀[J]." 化學(xué)工程師,2017,31(6):50-53.

    [36] ZHANG Y, GAO X, LI J, et al. Comparative analysis of proso millet(Panicum miliaceum L.) leaf transcriptomes for insight into drought tolerance mechanisms[J]." BMC Plant Biology, 2019, 19: 397.

    [37] SHAN Z, JIANG Y, LI H, et al." Genome-wide analysis of the NA Ctranscription factor family in broomcornmillet (Panicum miliaceum L.) and expression analysis under drought stress[J]." BMC Genomics, 2020, 21: 96.

    [38] 盤婉向,劉天鵬,何繼紅,等." 糜子(Panicum miliaceum L.)全基因組 YABBY 基因家族鑒定與高滲溶液脅迫下表達(dá)特征[J]." 基因組學(xué)與應(yīng)用生物學(xué),2022,41(5):144-155.

    [39] 王" "媚,劉天鵬,何繼紅,等. 糜子bZIP基因家族鑒定及幼苗期聚乙二醇6000處理下的表達(dá)特征[J]." 應(yīng)用與環(huán)境生物學(xué)報,2022,28(4):1-12.

    [40] 景海春,田志喜,種" "康,等." 分子設(shè)計育種的科技問題及其展望概論[J]." 中國科學(xué):生命科學(xué),2021,51(10):1356-1365;1355.

    [41] 王文廣,王永紅." 作物株型與產(chǎn)量研究進(jìn)展與展望[J]." 中國科學(xué): 生命科學(xué),2021,51(10):1366-1375.

    [42] BRUTNELL T P, WANG L, SWARTWOOD K, et al. Setariaviridis: A Model for C4 Photosynthesis[J]." The Plant Cell, 2010, 22: 2537-2544.

    [43] DAS S, KHOUND R, SANTRA M, et al. Beyond Bird Feed: Proso Millet for Human Health and Environment[J]." Agriculture, 2019, 9: 64.

    收稿日期:2022 - 09 - 30

    基金項目:國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系項目(CARS06-14.5-A8);甘肅省農(nóng)業(yè)科學(xué)院現(xiàn)代生物育種項目(2021GAAS02);甘肅省拔尖人才項目(2021);甘肅省青年科技基金計劃項目(20JR10RA457)。

    作者簡介:楊天育(1968 — ),男,甘肅渭源人,研究員,主要從事小雜糧育種與種質(zhì)資源研究工作。Email:13519638111@163.com。

    99精品久久久久人妻精品| 黑人巨大精品欧美一区二区蜜桃| 日日爽夜夜爽网站| 夜夜夜夜夜久久久久| 麻豆av在线久日| 69av精品久久久久久| 国产欧美日韩精品亚洲av| 久久久国产欧美日韩av| 妹子高潮喷水视频| 老司机深夜福利视频在线观看| 精品国内亚洲2022精品成人 | bbb黄色大片| 高清毛片免费观看视频网站 | 免费在线观看影片大全网站| 欧美激情久久久久久爽电影 | 精品亚洲成a人片在线观看| 成人特级黄色片久久久久久久| 色综合婷婷激情| 女警被强在线播放| 一级片免费观看大全| 啦啦啦 在线观看视频| 国产成人免费无遮挡视频| 狠狠狠狠99中文字幕| 国产欧美日韩精品亚洲av| 老汉色av国产亚洲站长工具| 久久久精品免费免费高清| 日韩免费av在线播放| 麻豆乱淫一区二区| 高清欧美精品videossex| 久久青草综合色| 欧美国产精品一级二级三级| 性少妇av在线| 欧美日韩一级在线毛片| 这个男人来自地球电影免费观看| 午夜福利在线免费观看网站| 亚洲成av片中文字幕在线观看| 精品一品国产午夜福利视频| 精品久久久久久,| 国产精品影院久久| 国产野战对白在线观看| 亚洲精品在线观看二区| 中文字幕制服av| 老司机靠b影院| 国产99白浆流出| 欧美精品啪啪一区二区三区| 人人澡人人妻人| 久久久久精品国产欧美久久久| 91大片在线观看| 热99re8久久精品国产| 桃红色精品国产亚洲av| 女人精品久久久久毛片| 亚洲精品国产一区二区精华液| 欧美日韩福利视频一区二区| 中文字幕人妻熟女乱码| 亚洲五月婷婷丁香| 国产黄色免费在线视频| av网站免费在线观看视频| 国产av精品麻豆| 国产午夜精品久久久久久| 久久久水蜜桃国产精品网| 高清视频免费观看一区二区| 久久国产精品男人的天堂亚洲| 午夜福利视频在线观看免费| 亚洲美女黄片视频| 精品亚洲成a人片在线观看| 国产日韩欧美亚洲二区| 极品少妇高潮喷水抽搐| av中文乱码字幕在线| 999精品在线视频| 亚洲精品国产区一区二| 久久婷婷成人综合色麻豆| 纯流量卡能插随身wifi吗| 欧美 日韩 精品 国产| 中文字幕制服av| 国产精品综合久久久久久久免费 | 人人妻,人人澡人人爽秒播| 午夜91福利影院| 国产精品 欧美亚洲| 成人三级做爰电影| 啪啪无遮挡十八禁网站| 精品福利永久在线观看| 久久亚洲精品不卡| 午夜福利乱码中文字幕| 亚洲精品一二三| av欧美777| 欧美黑人精品巨大| 精品一区二区三卡| 天天影视国产精品| 天天躁夜夜躁狠狠躁躁| 999久久久国产精品视频| 丝袜在线中文字幕| 丝袜在线中文字幕| 999久久久国产精品视频| 天天躁夜夜躁狠狠躁躁| 国产男靠女视频免费网站| 国产野战对白在线观看| 巨乳人妻的诱惑在线观看| 老司机在亚洲福利影院| 亚洲熟女精品中文字幕| 免费看十八禁软件| 欧美亚洲 丝袜 人妻 在线| 国产精品秋霞免费鲁丝片| 极品人妻少妇av视频| 美女高潮到喷水免费观看| 国产亚洲精品久久久久久毛片 | 80岁老熟妇乱子伦牲交| 午夜免费成人在线视频| 亚洲精品在线观看二区| 欧美日韩亚洲高清精品| 国产精品亚洲一级av第二区| 精品福利观看| 免费日韩欧美在线观看| 成人永久免费在线观看视频| 亚洲av片天天在线观看| 女人爽到高潮嗷嗷叫在线视频| 无限看片的www在线观看| 久久精品亚洲熟妇少妇任你| 午夜激情av网站| 建设人人有责人人尽责人人享有的| 国产成人欧美| 天天躁夜夜躁狠狠躁躁| 午夜视频精品福利| 日本五十路高清| 波多野结衣av一区二区av| 天天躁狠狠躁夜夜躁狠狠躁| 国产免费av片在线观看野外av| 久久亚洲精品不卡| 亚洲男人天堂网一区| 国产精品.久久久| 少妇 在线观看| 最近最新中文字幕大全电影3 | 视频区欧美日本亚洲| 天堂动漫精品| 成人免费观看视频高清| 亚洲欧美一区二区三区黑人| 一级,二级,三级黄色视频| 大香蕉久久成人网| 黄色片一级片一级黄色片| 欧美丝袜亚洲另类 | 亚洲精品一卡2卡三卡4卡5卡| 丁香六月欧美| 建设人人有责人人尽责人人享有的| 精品国产国语对白av| 精品一区二区三区视频在线观看免费 | 亚洲色图av天堂| 免费在线观看日本一区| 我的亚洲天堂| xxxhd国产人妻xxx| 看黄色毛片网站| 咕卡用的链子| 黑人巨大精品欧美一区二区mp4| 亚洲精品国产区一区二| 久久精品国产a三级三级三级| 日韩欧美国产一区二区入口| 亚洲伊人色综图| 欧美在线黄色| 欧美黑人精品巨大| 亚洲精品一二三| 国产精品九九99| 国产成人精品在线电影| 国产男女超爽视频在线观看| 满18在线观看网站| 搡老岳熟女国产| 精品视频人人做人人爽| 黄片小视频在线播放| 99在线人妻在线中文字幕 | 人妻丰满熟妇av一区二区三区 | 天堂俺去俺来也www色官网| 一区二区三区激情视频| 韩国精品一区二区三区| 日本欧美视频一区| 日本精品一区二区三区蜜桃| 久久99一区二区三区| 国产欧美日韩一区二区精品| 国产欧美日韩一区二区精品| 国产欧美日韩综合在线一区二区| 国产精品秋霞免费鲁丝片| 自拍欧美九色日韩亚洲蝌蚪91| 50天的宝宝边吃奶边哭怎么回事| 一本综合久久免费| 国产精品98久久久久久宅男小说| 国产无遮挡羞羞视频在线观看| 精品人妻1区二区| 欧美黑人精品巨大| 五月开心婷婷网| 狠狠狠狠99中文字幕| 免费在线观看黄色视频的| 在线播放国产精品三级| 91精品国产国语对白视频| 看黄色毛片网站| 国产精品国产高清国产av | 久久精品成人免费网站| 一进一出抽搐动态| 大码成人一级视频| 中文字幕av电影在线播放| 亚洲一区中文字幕在线| 黑人操中国人逼视频| 另类亚洲欧美激情| 午夜福利免费观看在线| 男人的好看免费观看在线视频 | 亚洲第一青青草原| 亚洲成人国产一区在线观看| 欧美日本中文国产一区发布| 精品国产一区二区三区久久久樱花| 亚洲五月天丁香| 无人区码免费观看不卡| 黄色丝袜av网址大全| 国产亚洲精品一区二区www | 亚洲精华国产精华精| 久99久视频精品免费| 99久久综合精品五月天人人| av国产精品久久久久影院| 国产精品98久久久久久宅男小说| 大型黄色视频在线免费观看| 成年动漫av网址| 国产精品九九99| a级毛片在线看网站| 热99re8久久精品国产| 亚洲综合色网址| 亚洲欧美日韩高清在线视频| 99国产精品一区二区三区| 少妇粗大呻吟视频| 啦啦啦免费观看视频1| 一二三四社区在线视频社区8| а√天堂www在线а√下载 | 一本综合久久免费| 女同久久另类99精品国产91| 午夜视频精品福利| 亚洲精品久久成人aⅴ小说| 大陆偷拍与自拍| 十八禁人妻一区二区| e午夜精品久久久久久久| 久久热在线av| 人妻 亚洲 视频| 成人免费观看视频高清| 色精品久久人妻99蜜桃| 两人在一起打扑克的视频| 欧美中文综合在线视频| 一区二区三区国产精品乱码| 国产视频一区二区在线看| 极品少妇高潮喷水抽搐| 中文字幕人妻丝袜一区二区| 亚洲av欧美aⅴ国产| 欧美日韩一级在线毛片| 成人18禁在线播放| 亚洲色图综合在线观看| 手机成人av网站| 国产欧美日韩一区二区三区在线| 国产欧美日韩一区二区三| 国产精品av久久久久免费| 欧美最黄视频在线播放免费 | 嫁个100分男人电影在线观看| 亚洲色图 男人天堂 中文字幕| 婷婷丁香在线五月| 9191精品国产免费久久| 美女视频免费永久观看网站| 国产av一区二区精品久久| 在线看a的网站| 亚洲成国产人片在线观看| 亚洲精品av麻豆狂野| 中文亚洲av片在线观看爽 | 午夜精品久久久久久毛片777| 成熟少妇高潮喷水视频| 王馨瑶露胸无遮挡在线观看| 欧美日韩黄片免| 国产高清激情床上av| 狠狠婷婷综合久久久久久88av| 中文欧美无线码| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 久久精品国产a三级三级三级| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品二区激情视频| 久久人妻福利社区极品人妻图片| 香蕉国产在线看| av一本久久久久| 亚洲一区中文字幕在线| videos熟女内射| 国产三级黄色录像| 高清毛片免费观看视频网站 | 老司机影院毛片| 中文字幕人妻熟女乱码| 欧美精品av麻豆av| 亚洲中文日韩欧美视频| 午夜久久久在线观看| 美女福利国产在线| 精品电影一区二区在线| 99在线人妻在线中文字幕 | av电影中文网址| 成人特级黄色片久久久久久久| 国产真人三级小视频在线观看| 亚洲第一青青草原| 天天操日日干夜夜撸| 午夜老司机福利片| 久久人人爽av亚洲精品天堂| 欧美乱色亚洲激情| 飞空精品影院首页| 国产主播在线观看一区二区| 亚洲一码二码三码区别大吗| 成人国语在线视频| 天天躁日日躁夜夜躁夜夜| 在线观看免费视频网站a站| 中亚洲国语对白在线视频| 精品一区二区三区av网在线观看| 少妇 在线观看| 久久狼人影院| 亚洲av熟女| 老司机午夜福利在线观看视频| 美女福利国产在线| 精品一区二区三区视频在线观看免费 | 熟女少妇亚洲综合色aaa.| 99精品在免费线老司机午夜| 欧美国产精品va在线观看不卡| 自线自在国产av| 亚洲男人天堂网一区| 18禁裸乳无遮挡免费网站照片 | 午夜影院日韩av| 国产午夜精品久久久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 国产99久久九九免费精品| 中文亚洲av片在线观看爽 | 操美女的视频在线观看| 高潮久久久久久久久久久不卡| 国产精品98久久久久久宅男小说| 极品少妇高潮喷水抽搐| 女人久久www免费人成看片| 久久久国产成人免费| 国产成人免费无遮挡视频| 精品电影一区二区在线| 午夜福利免费观看在线| 亚洲欧美激情综合另类| 男女高潮啪啪啪动态图| 亚洲色图av天堂| 亚洲精品一二三| 视频在线观看一区二区三区| 久久国产精品人妻蜜桃| 麻豆乱淫一区二区| 天堂√8在线中文| 亚洲人成电影免费在线| 亚洲精品av麻豆狂野| 纯流量卡能插随身wifi吗| 人成视频在线观看免费观看| xxx96com| 久久久国产欧美日韩av| 嫁个100分男人电影在线观看| 99久久99久久久精品蜜桃| videos熟女内射| 亚洲国产欧美一区二区综合| 久久 成人 亚洲| 久久精品国产清高在天天线| 亚洲精品自拍成人| 黑人猛操日本美女一级片| 国产精品亚洲一级av第二区| 国产91精品成人一区二区三区| 午夜福利在线观看吧| 精品第一国产精品| 久久人人97超碰香蕉20202| 国内毛片毛片毛片毛片毛片| 精品欧美一区二区三区在线| av不卡在线播放| 一区福利在线观看| 日韩大码丰满熟妇| 黄色怎么调成土黄色| 亚洲av成人不卡在线观看播放网| 久久人妻福利社区极品人妻图片| 婷婷丁香在线五月| 亚洲专区国产一区二区| 日韩视频一区二区在线观看| www.熟女人妻精品国产| 亚洲精品国产区一区二| 久久久久久久午夜电影 | 国产精华一区二区三区| 黄色a级毛片大全视频| 女人被躁到高潮嗷嗷叫费观| videosex国产| 在线观看日韩欧美| 一级片'在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| a级片在线免费高清观看视频| 久久精品亚洲av国产电影网| 精品久久久精品久久久| 男女床上黄色一级片免费看| 国产高清videossex| 国产欧美日韩一区二区三区在线| 免费高清在线观看日韩| 日韩欧美在线二视频 | 国产精品偷伦视频观看了| 精品国产国语对白av| 精品人妻在线不人妻| 亚洲黑人精品在线| 成熟少妇高潮喷水视频| 九色亚洲精品在线播放| 亚洲精品美女久久av网站| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| 一进一出好大好爽视频| 18禁裸乳无遮挡免费网站照片 | 国产精品偷伦视频观看了| 免费一级毛片在线播放高清视频 | 青草久久国产| 精品少妇久久久久久888优播| 自线自在国产av| 满18在线观看网站| 精品一区二区三区av网在线观看| 三上悠亚av全集在线观看| 精品国产乱码久久久久久男人| 欧美成人免费av一区二区三区 | 三上悠亚av全集在线观看| 国产av精品麻豆| 国产单亲对白刺激| 亚洲欧美精品综合一区二区三区| 黄色女人牲交| 日本一区二区免费在线视频| 欧美国产精品va在线观看不卡| 美女高潮喷水抽搐中文字幕| 99久久国产精品久久久| 久久久久国产精品人妻aⅴ院 | 老司机午夜福利在线观看视频| 中文字幕制服av| 亚洲国产欧美一区二区综合| 日韩免费高清中文字幕av| 久久久久精品人妻al黑| 首页视频小说图片口味搜索| 久久国产精品影院| 欧美激情极品国产一区二区三区| 国产成人免费观看mmmm| 欧美日韩精品网址| 日韩欧美三级三区| 日本wwww免费看| 亚洲精品国产色婷婷电影| 一本大道久久a久久精品| 日韩制服丝袜自拍偷拍| 一进一出好大好爽视频| 最近最新中文字幕大全免费视频| 免费在线观看黄色视频的| 高潮久久久久久久久久久不卡| 亚洲色图av天堂| 免费黄频网站在线观看国产| 国产男女超爽视频在线观看| 丰满的人妻完整版| 一夜夜www| 国产区一区二久久| 精品卡一卡二卡四卡免费| 日韩 欧美 亚洲 中文字幕| 我的亚洲天堂| 久久国产乱子伦精品免费另类| 69av精品久久久久久| 看免费av毛片| 亚洲国产精品sss在线观看 | 在线永久观看黄色视频| 美女 人体艺术 gogo| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| 熟女少妇亚洲综合色aaa.| 少妇粗大呻吟视频| av超薄肉色丝袜交足视频| 美女福利国产在线| 国产精品香港三级国产av潘金莲| 69精品国产乱码久久久| 老司机亚洲免费影院| 日本精品一区二区三区蜜桃| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品久久久久久毛片 | xxxhd国产人妻xxx| 亚洲伊人色综图| 超色免费av| 国产亚洲精品一区二区www | www.精华液| 亚洲伊人色综图| 18禁裸乳无遮挡免费网站照片 | av天堂久久9| 日韩制服丝袜自拍偷拍| 天天操日日干夜夜撸| 校园春色视频在线观看| 亚洲精品一二三| 啪啪无遮挡十八禁网站| 99国产精品99久久久久| 91麻豆精品激情在线观看国产 | 极品人妻少妇av视频| 成人特级黄色片久久久久久久| 国产av精品麻豆| 欧美久久黑人一区二区| 啦啦啦免费观看视频1| 操出白浆在线播放| 国产视频一区二区在线看| 国产伦人伦偷精品视频| 日日摸夜夜添夜夜添小说| 1024香蕉在线观看| 色婷婷av一区二区三区视频| 欧美亚洲日本最大视频资源| cao死你这个sao货| 美女福利国产在线| 欧美午夜高清在线| 女人被狂操c到高潮| 日韩免费av在线播放| 18在线观看网站| 老司机在亚洲福利影院| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看 | 麻豆av在线久日| 无限看片的www在线观看| 成年人免费黄色播放视频| 亚洲avbb在线观看| 好男人电影高清在线观看| 一级片免费观看大全| 韩国av一区二区三区四区| 亚洲成av片中文字幕在线观看| 国产成人啪精品午夜网站| 免费少妇av软件| 国产精品秋霞免费鲁丝片| 亚洲熟女毛片儿| 亚洲一码二码三码区别大吗| 咕卡用的链子| 岛国毛片在线播放| xxx96com| av视频免费观看在线观看| 夜夜夜夜夜久久久久| 日韩免费av在线播放| 欧美乱妇无乱码| 一级片'在线观看视频| www.自偷自拍.com| 国产熟女午夜一区二区三区| 国产乱人伦免费视频| 免费观看人在逋| 午夜福利影视在线免费观看| 大陆偷拍与自拍| 美国免费a级毛片| 欧美日韩亚洲综合一区二区三区_| 欧美乱码精品一区二区三区| 国产欧美亚洲国产| 久久国产精品人妻蜜桃| 久久国产精品男人的天堂亚洲| 久久久久久亚洲精品国产蜜桃av| 身体一侧抽搐| 桃红色精品国产亚洲av| 国产精品一区二区精品视频观看| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 国产亚洲欧美在线一区二区| 9色porny在线观看| 91麻豆精品激情在线观看国产 | 一进一出抽搐gif免费好疼 | 亚洲精品自拍成人| 美女国产高潮福利片在线看| 免费观看人在逋| 成人永久免费在线观看视频| 色94色欧美一区二区| 黄色毛片三级朝国网站| 91精品三级在线观看| 高潮久久久久久久久久久不卡| 亚洲精品美女久久久久99蜜臀| xxxhd国产人妻xxx| 国产三级黄色录像| 久久 成人 亚洲| 在线天堂中文资源库| 一级毛片女人18水好多| 久久中文看片网| 亚洲精品国产色婷婷电影| 少妇被粗大的猛进出69影院| 成年女人毛片免费观看观看9 | 香蕉丝袜av| 国产在视频线精品| 国内毛片毛片毛片毛片毛片| 香蕉久久夜色| 久久中文字幕一级| 欧美乱色亚洲激情| 中文字幕另类日韩欧美亚洲嫩草| 搡老熟女国产l中国老女人| 一个人免费在线观看的高清视频| 国精品久久久久久国模美| 亚洲av成人av| 国产精品电影一区二区三区 | 黄片播放在线免费| 91精品三级在线观看| 久久青草综合色| 亚洲精品粉嫩美女一区| 国精品久久久久久国模美| 国产一区二区激情短视频| 丁香六月欧美| 黄色视频不卡| 午夜91福利影院| 国产三级黄色录像| 热re99久久国产66热| 一个人免费在线观看的高清视频| 亚洲午夜理论影院| 丝袜美腿诱惑在线| 欧美性长视频在线观看| 亚洲七黄色美女视频| 麻豆成人av在线观看| 久久九九热精品免费| 黄色 视频免费看| 中文字幕人妻熟女乱码| 999久久久国产精品视频| 精品无人区乱码1区二区| 日本五十路高清| 最新在线观看一区二区三区| 色在线成人网| 久久久国产一区二区| 久久亚洲精品不卡| 国内久久婷婷六月综合欲色啪| 国产欧美日韩一区二区三| 久久热在线av| 捣出白浆h1v1| 王馨瑶露胸无遮挡在线观看| 天堂√8在线中文| 久久人人97超碰香蕉20202| 日韩有码中文字幕| 嫩草影视91久久| 国产区一区二久久| 精品亚洲成a人片在线观看| 如日韩欧美国产精品一区二区三区| 国产精品久久久久成人av| 19禁男女啪啪无遮挡网站| 亚洲精品国产色婷婷电影| 国产xxxxx性猛交|