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

    三七基因組SSR位點(diǎn)分析和多態(tài)性引物開發(fā)

    2018-05-14 14:44:48揭應(yīng)碧盧迎春宋婉玲韋坤華張廣輝楊生超
    熱帶作物學(xué)報(bào) 2018年7期
    關(guān)鍵詞:三七多態(tài)性基因組

    揭應(yīng)碧 盧迎春 宋婉玲 韋坤華 張廣輝 楊生超

    摘 要 本研究通過分析三七(Panax notoginseng)基因組中簡單重復(fù)序列(simple sequence repeats,SSRs)位點(diǎn)信息,對基因組中SSR位點(diǎn)的分布特征進(jìn)行描述,并開發(fā)多態(tài)性引物,為評價(jià)三七遺傳多樣性和群體結(jié)構(gòu)提供理論依據(jù)。利用MISA軟件對重新組裝的三七基因組進(jìn)行SSR位點(diǎn)搜索,primer3設(shè)計(jì)引物。用隨機(jī)合成的200對引物,在16個(gè)三七單株中進(jìn)行PCR擴(kuò)增,來篩選具有多態(tài)性的引物。在去除冗余序列后組裝得到2.39 Gb大小的三七基因組中總共識(shí)別到314 060個(gè)SSR位點(diǎn),SSR分布頻率為0.13/kb,其中占主導(dǎo)地位的是二核苷酸重復(fù)基序,占所有重復(fù)基序類型的68.55%,其次是三核苷酸56 250(17.91%)和單核苷酸22 475(7.16%)。在二核苷酸為主的重復(fù)類型中,出現(xiàn)頻率最高的是AT/AT重復(fù)基序占所有二核苷酸重復(fù)的63.66%,最低的是CG/CG重復(fù)基序占0.15%。同時(shí),在單核苷酸重復(fù)中出現(xiàn)頻率最高的也是A/T重復(fù)類型(79.01%)。經(jīng)過多態(tài)性篩選最終得到41對多態(tài)性較好的引物,PIC值范圍為0.11~0.78,平均值為0.46。本研究獲得的大量SSR位點(diǎn)信息為三七全基因組SSR標(biāo)記的開發(fā),三七分子標(biāo)記輔助育種以及遺傳多樣性分析提供方便。

    關(guān)鍵詞 三七;基因組;SSR;多態(tài)性

    中圖分類號(hào) S31 文獻(xiàn)標(biāo)識(shí)碼 A

    Abstract The analysis of simple sequence repeat (simple sequence repeats; SSRs) site information in the genome of P. notoginseng was used to describe the distribution characteristics of genome SSR loci, and to develop polymorphic primers to evaluate the genetic diversity and population structure of P notoginseng. MISA software was used to search for the SSR loci of the assembled genome of P. notoginseng. Primers were designed by Primer3, and 200 pairs of primers randomly synthesized were used to amplify PCR in 16 P. notoginseng plants to screen primers with polymorphism. Totally, 314, 060 SSR loci were identified, with an average distance of 0.13 kb between each loci. A large proportion of the SSRs included di-, tri-and mono- repeat motifs, which accounted for 68.55%, 17.91% and 7.16% of all repeats, respectively. Dinucleotide is the main type of repeated categories, the highest frequency of occurrence was AT/AT which accounted for 63.66% in total di- repeats, whereas the lowest frequency of occurrence was CG/CG (0.15%). A/T was the most common in mononucleotide repeat (79.01%). After polymorphic screening, we obtained 41 primers with preferable polymorphic primers. The PIC value was 0.11~0.78, with an average of 0.46. Plenty of the SSR locus information and SSR markers were provided in this research, not only provide valuable resources for the development of SSR markers in P. notoginseng, but also provide convenience for P. notoginseng marker-assisted breeding and genetic diversity analysis.

    Keywords Panax notoginseng; genome; SSR; polymorphism

    DOI 10.3969/j.issn.1000-2561.2018.07.014

    簡單重復(fù)序列(simple sequence repeats,SSRs)或微衛(wèi)星標(biāo)記,是由1~6個(gè)核苷酸多次串聯(lián)重復(fù)的DNA序列,是基因組中豐富存在的、均勻分布的共顯性和多等位基因的分子標(biāo)記,而且容易被檢測[1~4]。這些特性使得SSR標(biāo)記成為植物分子標(biāo)記輔助育種、遺傳資源DNA指紋圖譜分析、圖位克隆基因等研究的重要工具之一。最初,大多數(shù)植物的SSR標(biāo)記主要來自表達(dá)序列標(biāo)簽(expressed sequence tags,ESTs)和細(xì)菌人工染色體(bacterial artificial chromosome,BAC)序列。近年來,隨著越來越多的植物基因組測序的完成,也從基因組中鑒定了大量SSRs標(biāo)記[5~8]。但是,多集中于模式植物和主要的農(nóng)作物等,其他植物由于缺乏基因組序列,基于基因組的SSR標(biāo)記開發(fā)也鮮有報(bào)道。

    三七(Panax notoginseng (Burk.) F. H. Chen)為五加科人參屬多年生草本植物,以根及根莖入藥,具有散瘀止血、消腫定痛、防止腦缺血、止血、降脂和保護(hù)肝臟等功效[9-10]。據(jù)考證,三七至少已有400年以上的栽培歷史[11]。但是,幾百年來,一直沒有選育出用于推廣種植的品種,生產(chǎn)上長期使用混雜的種質(zhì),導(dǎo)致三七遺傳多樣性逐漸喪失[12],直接影響三七藥材的產(chǎn)量和品質(zhì),表現(xiàn)為三七植株個(gè)體間表型差異巨大[13],主要皂苷組分含量也相差3.6~10.2倍[14]。同時(shí),居群內(nèi)遺傳多樣性豐富,而居群間遺傳分化水平低,遺傳差異主要存在于居群內(nèi)[15-16]。因此,三七新品種選育是生產(chǎn)上亟待解決的關(guān)鍵科學(xué)問題之一。

    開發(fā)有效的分子標(biāo)記是植物育種的重要手段之一。前人曾利用fAFLP(fluorescent amplified fragment length polymorphism)和EST-SSR等對三七群體的遺傳多樣性進(jìn)行分析[15-16]。楊云[17]利用ISSR(inter-simple sequence repeat)和SRAP(sequence-related amplified polymorphism)對三七集團(tuán)選擇后代群體進(jìn)行了遺傳一致性分析,并選育了2個(gè)新品種“滇七1號(hào)”和“苗鄉(xiāng)三七1號(hào)”,其中“滇七1號(hào)”的主要特征為莖翠綠色、存苗率高;“苗鄉(xiāng)三七1號(hào)”莖則為紫色,抗性強(qiáng),皂苷含量高。近年來,利用轉(zhuǎn)錄組測序開發(fā)了人參屬植物大量EST-SSR標(biāo)記[18-20],但EST-SSR的覆蓋度低、不包含非基因非編碼區(qū),從而限制了其應(yīng)用。

    2017年,中藥材課題組率先發(fā)布了三七基因組(包括轉(zhuǎn)錄本),組裝的三七基因組草圖大小為2.39 Gb,contig N50和scaffold N50分別為16 kb和96 kb[21]?;蚪M測序的完成為三七育種和藥效成分生物合成關(guān)鍵酶基因的發(fā)現(xiàn)提供了基礎(chǔ)。本研究分析比較了三七基因組和轉(zhuǎn)錄組的SSR位點(diǎn)信息差異,利用三七基因組SSR位點(diǎn)開發(fā)了大量引物,并隨機(jī)選擇了200對引物進(jìn)行了多態(tài)性驗(yàn)證,從中篩選出了41對多態(tài)性引物。本研究為三七SSR標(biāo)記開發(fā)提供堅(jiān)實(shí)基礎(chǔ),也為三七種源遺傳多樣性鑒定等提供了有效的工具。

    1 材料與方法

    1.1 三七基因組中SSR位點(diǎn)信息的搜索

    利用軟件MISA(microsatellite identification tool, http://pgrc.ipk-gatersleben.de/misa/)搜索三七基因組中的SSR位點(diǎn)。搜索標(biāo)準(zhǔn)選擇重復(fù)單元為1、2、3、4、5、6個(gè)堿基的基序,最小重復(fù)數(shù)分別為12、6、5、5、4、4。2個(gè)相鄰SSR之間的堿基數(shù)≤100,則被定義為復(fù)合型SSR。為比較轉(zhuǎn)錄組中SSR位點(diǎn)和基因組SSR位點(diǎn)差異,本研究下載了三七不同部位轉(zhuǎn)錄組原始測序reads[19](NCBI Sequence Read Archive登錄號(hào)為PRJNA228978),混合組裝,并分析其SSR位點(diǎn)(不包含單核苷酸重復(fù))。

    1.2 三七SSR引物設(shè)計(jì)與合成

    根據(jù)SSR位點(diǎn)兩端的保守序列,利用primer3(http://primer3.sourceforge.net/)設(shè)計(jì)引物,引物長度20~23 bp,退火溫度60 °C左右;引物中堿基重復(fù)數(shù)小于等于4,引物兩端不能有2個(gè)連續(xù)的A/T堿基;PCR產(chǎn)物大小在150~300 bp。隨機(jī)選擇200對引物,由昆明碩陽科技有限公司合成,以篩選多態(tài)性引物。

    1.3 引物多態(tài)性驗(yàn)證

    在云南文山市苗鄉(xiāng)三七科技園(文山市硯山縣),隨機(jī)采集16個(gè)單株的新鮮葉片用于多態(tài)性引物篩選,其中自然群體6個(gè)單株、苗鄉(xiāng)三七1號(hào)6個(gè)單株和滇七1號(hào)4個(gè)單株。EasyPure Plant Genomic DNA Kit(北京全式金生物技術(shù)有限公司)提取基因組DNA。引物合成時(shí)在正向引物5′端加M13F通用引物序列(5′-TGTAAAACGA?CGG?CCAGT-3′),第一次PCR擴(kuò)增后產(chǎn)物帶上M13F序列。選擇第一次PCR中條帶理想的產(chǎn)物進(jìn)行二次擴(kuò)增,用帶有不同顏色且序列和M13F一致的熒光引物(A603012557Hex和A6030125?56?Fam)替換原來的M13F序列,目的是使得二次擴(kuò)增產(chǎn)物帶有熒光便于檢測。

    第一輪PCR擴(kuò)增的反應(yīng)體系為15 ?L,包含7.5 μL mix,正向和反向引物各1.0 μL,基因組DNA 1.0 μL和4.5 μL ddH2O;擴(kuò)增程序94 ℃預(yù)變性5 min,然后進(jìn)行30個(gè)循環(huán),每個(gè)循環(huán)包含94 ℃變性30 s,最佳退火溫度退火30 s,72 ℃延伸30 s;最后72 ℃延伸10 min。擴(kuò)增產(chǎn)物0.8%的瓊脂糖凝膠電泳檢測。以第一輪PCR的產(chǎn)物為模板進(jìn)行第二輪PCR擴(kuò)增,15 μL PCR反應(yīng)體系包含7.5 μL mix,1 μL反向引物、一次擴(kuò)增產(chǎn)物1.0 μL、熒光引物1.0 μL和4.5 μL ddH2O。擴(kuò)增程序和一次擴(kuò)增相同,將二次擴(kuò)增的產(chǎn)物進(jìn)行0.8%瓊脂糖凝膠電泳檢測。將PCR產(chǎn)物進(jìn)行相應(yīng)的稀釋定量,將稀釋后的兩種不同的熒光引物混合加入一個(gè)孔內(nèi),體積各為0.5 μL,加入10 μL的GS500LIZ組成11 μL反應(yīng)體系,離心,95 ℃變性5 min,–20 ℃冰箱保存。

    1.4 數(shù)據(jù)分析

    用軟件GeneMapper v4.1對測序結(jié)果分析位點(diǎn)多態(tài)性,根據(jù)產(chǎn)生的峰圖確定核心位點(diǎn)片段大小,選擇特異性高、多態(tài)性好的位點(diǎn)峰圖產(chǎn)生的Excel數(shù)據(jù),POPGENE32計(jì)算每個(gè)標(biāo)記的等位基因數(shù)、Shannon指數(shù),軟件PIC CALC計(jì)算每個(gè)標(biāo)記的多態(tài)性信息含量,統(tǒng)計(jì)分析時(shí)將SSR視為共顯性標(biāo)記。

    2 結(jié)果與分析

    2.1 三七基因組中SSR信息位點(diǎn)的分布特征

    三七基因組122 131條scaffold中共包含314 060個(gè)SSR位點(diǎn),分布于36 763條scaffold上,其中26 500條scaffold含有一個(gè)以上SSR位點(diǎn),復(fù)合型SSR有46 463個(gè)(表1)。相反,在三七轉(zhuǎn)錄組的45 678條unigenes中,識(shí)別到了4 756個(gè)SSR位點(diǎn),同樣地,在454測序的三七根部轉(zhuǎn)錄組的2 361條特異序列中識(shí)別到2 772個(gè)SSR基序。在三七轉(zhuǎn)錄組中每Mb有134.73個(gè)SSR,位點(diǎn)發(fā)生頻率為7.43 kb/SSR。三七基因組每Mb有131.00個(gè)SSR,SSR位點(diǎn)的發(fā)生頻率為7.62 kb/SSR。

    三七基因組中,二核苷酸重復(fù)單元占有絕對優(yōu)勢,占總數(shù)的68.55%;其次是三核苷酸重復(fù),占總數(shù)的17.91%;同樣,在三七轉(zhuǎn)錄組中,二核甘酸重復(fù)也最多(40.81%),三核苷酸重復(fù)所占比例也很高(36.67%)。在三七根部轉(zhuǎn)錄組中二核苷酸重復(fù)(67.53%)也是最多的,幾乎是三核苷酸重復(fù)(25%)的3倍(表2)。和許多植物基因組相似,三七基因組中較短重復(fù)基序(mono-、di-、tri-)所占比例明顯高于較長重復(fù)基序(tetra-、penta-、hexa-)所占比例,這一結(jié)果表明在大多數(shù)植物中基序長度的分布基本一致。從SSR重復(fù)次數(shù)來看,基因組和轉(zhuǎn)錄組中的SSR位點(diǎn)都是6次重復(fù)最多,分別是11.64%和27.12%;5次重復(fù)和7次重復(fù)的比例也很高。特別地,轉(zhuǎn)錄組中4次重復(fù)的SSR的比例也很高,達(dá)到16.97%(表2),但是不論在基因組還是轉(zhuǎn)錄組中SSR頻率的整體變化規(guī)律是隨著基序長度的增加而減少。

    在單核苷酸重復(fù)中和C/G相比出現(xiàn)數(shù)目最多的是A/T重復(fù)單元,占所有單核苷酸重復(fù)的79.01%。在二核苷酸重復(fù)中出現(xiàn)數(shù)目最多的是AT/AT重復(fù)基序,占所有二核苷酸重復(fù)基序的63.66%;在三七轉(zhuǎn)錄組中除了AT/AT還有AG/CT。三七基因組中二核苷酸重復(fù)類型中數(shù)目最少的是CG/CG,只有330個(gè),占0.15%。在三七基因組三核苷酸重復(fù)序列中,出現(xiàn)最普遍的重復(fù)類型是AAT/ATT(53.33%)和AAG/CTT(21.91%),出現(xiàn)數(shù)目最少的是ACG/CGT和CCG/CGG,分別只有145和319個(gè),占所有三核苷酸重復(fù)基序的0.82%(圖1)。然而在三七轉(zhuǎn)錄組的三核苷酸重復(fù)中,AAG/CTT所占比例最高,CCG/CGG和AAC/GTT最少。

    2.2 SSR多態(tài)性引物的篩選和驗(yàn)證

    本研究隨機(jī)選擇并合成200對三七引物對樣本位點(diǎn)進(jìn)行多態(tài)性檢測,將所有PCR擴(kuò)增產(chǎn)物進(jìn)行毛細(xì)管電泳之后用Genemapper v4.1產(chǎn)生的峰圖檢測樣品位點(diǎn)信息,SSR重復(fù)基序的重復(fù)數(shù)目差異以及短的插入缺失事件形成SSR多態(tài)性,峰圖所顯示的信息是帶有熒光引物的PCR產(chǎn)物在經(jīng)過毛細(xì)管電泳之后所反饋的產(chǎn)物片段信息。通過分析每一對引物在各個(gè)樣品中產(chǎn)生的峰圖顯示的各自特定位點(diǎn)片段大小信息之后,其中有27對引物的峰圖大小相同不存在差異,有132對引物沒有擴(kuò)增出片段,最終得到的有41對引物的產(chǎn)物符合預(yù)期大小,將這41對引物峰圖生成Excel數(shù)據(jù)表格,用于計(jì)算引物多態(tài)性信息含量(圖2)。

    通過軟件POPGENE32計(jì)算41對引物的等位基因觀察值、預(yù)期等位基因數(shù)、Shannon指數(shù)和多態(tài)性信息含量,各個(gè)位點(diǎn)等位基因數(shù)目范圍為2~10,平均數(shù)目為3.9。計(jì)算了41對引物的PIC值,結(jié)果見表3,結(jié)果發(fā)現(xiàn)16對引物具有高度多態(tài)性(PIC>0.50),23對引物具有中度多態(tài)性(0.25

    3 討論

    SSR具有共顯性、豐富性和高多態(tài)性的特點(diǎn),并且普遍存在于植物生物基因組中。基于不同植物基因組的調(diào)查發(fā)現(xiàn),在擬南芥(Arabidopsis thaliana (L.) Heynh.)中平均1.14 kb發(fā)現(xiàn)一個(gè)SSR[22],在水稻(Oryza sativa L.)中平均3.6 kb出現(xiàn)一個(gè)SSR[23],在甘藍(lán)(Brassica oleracea L.)中平均4 kb出現(xiàn)一個(gè)SSR[24],在大豆(Glycine max (Linn.) Merr.)中平均4.5 kb出現(xiàn)一個(gè)SSR[25],在高粱(Sorghum bicolor (L.) Moench)中平均220 kb出現(xiàn)一個(gè)SSR[26],在小麥(Triticum aestivum L.)中平均578 kb出現(xiàn)一個(gè)SSR[27]。本研究從三七基因組的122 131條序列、總長度為239 483 436 bp中識(shí)別到314 060個(gè)SSR信息位點(diǎn),分布在36 763條序列中,SSR的平均密度為每0.13 kb出現(xiàn)一個(gè)SSR。這些結(jié)果可能反映出了存在于植物基因組間的DNA水平的真實(shí)遺傳差異,或者各項(xiàng)研究使用了不同的測序方法和測序程序。

    三七的基因組中識(shí)別到的SSR基序以二核苷酸為主,占所有識(shí)別到的SSR基序的68.55%,四核苷酸、五核苷酸和六核苷酸所占比例的總和為6.39%。這一結(jié)果和許多植物中二核苷酸重復(fù)基序在植物基因組中所占比例較高是一致的[28-31]。和許多植物中單核苷酸重復(fù)占主導(dǎo)位置的結(jié)果不同[32],這一結(jié)果表明SSR在不同物種以及相同物種的不同重復(fù)類型之間是存在差異的。從另一方面講,擁有大量短重復(fù)類型的SSR位點(diǎn)的物種通常表現(xiàn)出較高的基因組突變率[33-35]。

    在單核苷酸重復(fù)中A/T基序明顯高于C/G重復(fù),同時(shí)本研究發(fā)現(xiàn)AT/AT重復(fù)基序是二核苷酸重復(fù)類型中出現(xiàn)頻率最高的,然而CG/GC重復(fù)基序是出現(xiàn)頻率最低的,這一結(jié)果和毛果楊(Populus trichocarpa Torr. & Gray)和落花生(Arachis hypogaea Linn.)中AT/AT重復(fù)基序在二核苷酸中出現(xiàn)頻率最高是一致的[36-37]。許多研究表明,在植物基因組中二核苷酸基序中,(AT)n是最普遍的[37-38]。這一結(jié)果和許多植物基因組中富含AT堿基的重復(fù)類型所占比例高于富含CG的重復(fù)類型是一致的[39-43]。

    綜上所述,本研究從三七的基因組中分析得到了豐富的SSR位點(diǎn)信息。對200對引物進(jìn)行多態(tài)性篩選后,計(jì)算PIC值范圍是0.06~0.78,平均值為0.46,最終獲得41對特異性高、多態(tài)性好的引物,多態(tài)性位點(diǎn)獲得率為20.5%。同時(shí)用3種不同居群的三七材料對篩選的41對引物進(jìn)行多態(tài)性驗(yàn)證,最終將材料大致分為3類:紫莖、綠莖和對照。其中部分對照分別和綠莖、紫莖聚在一起,這可能跟紫莖和綠莖三七是由對照三七經(jīng)過長期栽培選育所得有關(guān)。因此本研究獲得的引物具有較高的多態(tài)性潛能,可以用于區(qū)分不同居群性狀差異的三七,同時(shí)這些豐富的SSR位點(diǎn)將從基因?qū)用娼沂救呦嚓P(guān)性狀的本質(zhì),對進(jìn)一步研究三七的遺傳多樣性、加速三七功能基因資源的開發(fā)利用、豐富其分子標(biāo)記類型、遺傳資源評價(jià)、繪制遺傳圖譜、實(shí)現(xiàn)特定性狀的輔助選擇和進(jìn)行比較基因組學(xué)研究都具有重要的意義。

    致謝 感謝提供實(shí)驗(yàn)材料的云南文山苗鄉(xiāng)三七科技園,特此謝意!

    參考文獻(xiàn)

    [1] Powell W, Machray G C, Provan J. Polymorphism revealed by simple sequence repeats[J]. Trends in Plant Science, 1996, 1(7): 215-222.

    [2] Chambers G K, Macavoy E S. Microsatellites: consensus and controversy[J]. Comparative Biochemistry & Physiology Part B Biochemistry & Molecular Biology, 2000, 126(4): 455-476.

    [3] Ellegren H. Microsatellites: simple sequences with complex evolution[J]. Nature Reviews Genetics, 2004, 5(6): 435.

    [4] Buschiazzo E, Gemmell N J. The rise, fall and renaissance of microsatellites in eukaryotic genomes[J]. Bioessays News & Reviews in Molecular Cellular & Developmental Biology, 2006, 28(10): 1 040-1 050.

    [5] Cavagnaro P F, Senalik D A, Yang L M, et al. Genome-wide characterization of simple sequence repeats in cucumber (Cucumis sativus L.)[J]. BMC Genomics, 2010, 11(1): 569.

    [6] Sonah H, Deshmukh R K, Sharma A, et al. Genome-wide distribution and organization of microsatellites in plants: an insight into marker development in Brachypodium[J]. PLoS ONE, 2011, 6(6): e21298.

    [7] Han B, Wang C, Tang Z, et al. Genome-wide analysis of microsatellite markers based on sequenced database in Chinese spring wheat (Triticum aestivum L.)[J]. PLoS ONE, 2015, 10(11): e0141540.

    [8] Zhu H, Song P, Koo D H, et al. Genome wide characterization of simple sequence repeats in watermelon genome and their application in comparative mapping and genetic diversity analysis[J]. BM Genomics, 2016, 17(1): 557.

    [9] Ng T B. Pharmacological activity of sanchi ginseng (Panax notoginseng)[J]. Journal of Pharmacy & Pharmacology, 2006, 58(8): 1 007.

    [10] 國家藥典委員會(huì). 中國藥典:一部[M]. 北京: 中國醫(yī)藥科技出版社, 2010: 11-12.

    [11] Guo H B, Cui X M, An N, et al. Sanchi ginseng (Panax notoginseng (Burkill) F. H. Chen) in China: distribution, cultivation and variations[J]. Genetic Resources & Crop Evolution, 2010, 57(3): 453-460.

    [12] Zhou S L, Xiong G M, Li Z Y, et al. Loss of genetic diversity of domesticated panax notoginseng F H Chen as evidenced by ITS sequence and AFLP polymorphism: A comparative study with P. stipuleanatus H Tsai et K M Feng[J]. Journal of Integrative Plant Biology, 2005, 47(1): 107–115.

    [13] 張金渝, 肖 慧, 金 航, 等. 三七栽培居群間表型變異式樣研究[J]. 中國中藥雜志, 2009, 34(24): 3 295-3 298.

    [14] Hong D Y, Lau A J, Yeo C L, et al. Genetic diversity and variation of saponin contents in Panax notoginseng roots from a single farm.[J]. Journal of Agricultural & Food Chemistry, 2005, 53(22): 8 460-8 467.

    [15] Wang D, Hong D, Koh H L, et al. Biodiversity in cultivated Panax notoginseng populations[J]. 中國藥理學(xué)報(bào), 2008, 29(9): 1 137-1 140.

    [16] 張金渝, 楊維澤, 崔秀明. 三七栽培居群遺傳多樣性的EST-SSR分析[J]. 植物遺傳資源學(xué)報(bào), 2011, 12(2): 249-254.

    [17] 楊 云. 三七集團(tuán)選擇后代群體的鑒定與評價(jià)[D]. 昆明: 云南農(nóng)業(yè)大學(xué), 2013.

    [18] Luo H, Sun C, Sun Y, et al. Analysis of the transcriptome of Panax notoginseng root uncovers putative triterpene saponin-biosynthetic genes and genetic markers[J]. Bmc Genomics, 2011, 12(S5): 1-15.

    [19] Liu M H, Yang B R, Cheung W F, et al. Transcriptome analysis of leaves, roots and flowers of Panax notoginseng identifies genes involved in ginsenoside and alkaloid biosynthesis[J]. Bmc Genomics, 2015, 16(1): 265.

    [20] 鄒麗秋, 匡雪君, 李 瀅, 等. 人參屬藥用植物轉(zhuǎn)錄組研究進(jìn)展[J]. 中國中藥雜志, 2016, 62(22): 4 138-4 143.

    [21] Chen W, Kui L, Zhang G H, et al. Whole-genome sequencing and analysis of the Chinese herbal plant Panax notoginseng[J]. Molecular Plant, 2017, 10(6): 899-902.

    [22] Lawson M J, Zhang L. Distinct patterns of SSR distribution in the Arabidopsis thaliana and rice genomes[J]. Genome Biology, 2006, 7(2): R14.

    [23] Zhang Z, Deng Y, Tan J, et al. A genome-wide microsatellite polymorphism database for the Indica and Japonica rice[J]. DNA Research, 2007, 14(1): 37-45.

    [24] Iniguez-Luy F L, Voort A V, Osborn T C. Development of a set of public SSR markers derived from genomic sequence of a rapid cycling Brassica oleracea L. genotype[J]. Theoretical & Applied Genetics, 2008, 117(6): 977-985.

    [25] Song Q, Jia G, Zhu Y, et al. Abundance of SSR motifs and development of candidate polymorphic SSR markers (BARCSOYSSR_1.0) in soybean[J]. Crop Science, 2010, 50(5): 1 950-1 960.

    [26] Yonemaru J, Ando T, Mizubayashi T, et al. Development of genome-wide simple sequence repeat markers using whole-genome shotgun sequences of Sorghum (Sorghum bicolor (L.) Moench)[J]. DNA Research, 2009, 16(3): 187-193.

    [27] Morgante M, Hanafey M, Powell W. Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes[J]. Nature Genetics, 2002, 30(2): 194-200.

    [28] Dettori M T, Micali S, Giovinazzi J, et al. Mining microsatellites in the peach genome: development of new long-core SSR markers for genetic analyses in five Prunus, species[J]. Springerplus, 2015, 4(1): 337.

    [29] Liu S R, Li W Y, Long D, et al. Development and characterization of genomic and expressed SSRs in citrus by genome-wide analysis[J]. PLoS ONE, 2013, 8(10): e75149.

    [30] Zhu H, Song P, Koo D H, et al. Genome wide characterization of simple sequence repeats in watermelon genome and their application in comparative mapping and genetic diversity analysis[J]. BMC Genomics, 2016, 17(1): 557.

    [31] Zhang Q, Ma B, Li H, et al. Identification, characterization, and utilization of genome-wide simple sequence repeats to identify a QTL for acidity in apple[J]. BMC Genomics, 2012, 13(1): 537.

    [32] Sonah H, Deshmukh R K, Sharma A, et al. Genome-wide distribution and organization of microsatellites in plants: an insight into marker development in Brachypodium[J]. PLoS ONE, 2011, 6(6): e21298.

    [33] Tóth G, Gáspári Z, Jurka J. Microsatellites in different eukaryotic genomes: survey and analysis[J]. Genome Research, 2000, 10(7): 967.

    [34] Sia E A, Kokoska R J, Dominska M, et al. Microsatellite instability in yeast: dependence on repeat unit size and DNA mismatch repair genes[J]. Molecular & Cellular Biology, 1997, 17(5): 2 851.

    [35] Nagata T, Kaho K, Nishikawa S, et al. Long microsatellite alleles in drosophila melanogaster have a downward mutation bias and short[J]. Genetics, 2000, 155(3): 1 213.

    [36] Tuskan G A, Gunter L E, Yang Z K, et al. Characterization of microsatellites revealed by genomic sequencing of populus trichocarpa[J]. Canadian Journal of Forest Research, 2013, 34(1): 85-93.

    [37] Hui W, Penmetsa R V, Mei Y, et al. Development and characterization of BAC-end sequence derived SSRs, and their incorporation into a new higher density genetic map for cultivated peanut (Arachis hypogaea L.)[J]. BMC Plant Biology, 2012, 12(1): 10.

    [38] Varshney R K, Graner A, Sorrells M E. Genic microsatellite markers in plants: features and applications[J]. Trends in Biotechnology, 2005, 23(1): 48.

    [39] Xiao J, Zhao J, Liu M, et al. Genome-wide characterization of simple sequence repeat (SSR) loci in Chinese jujube and jujube SSR primer transferability[J]. PLoS ONE, 2015, 10(5): e0127812.

    [40] Biswas M K, Xu Q, Mayer C, et al. Genome wide characterization of short tandem repeat markers in sweet orange (Citrus sinensis)[J]. PLoS ONE, 2014, 9(8): e104182.

    [41] Dossa K, Yu J, Liao B, et al. Development of highly informative genome-wide single sequence repeat markers for breeding applications in sesame and construction of a Web resource: SisatBase[J]. Frontiers in Plant Science, 2017, 8(1470): 1 470.

    [42] Cavagnaro P F, Senalik D A, Yang L M, et al. Genome-wide characterization of simple sequence repeats in cucumber (Cucumis sativus L.)[J]. BMC Genomics, 2010, 11(1): 569.

    [43] Cui J, Cheng J, Nong D, et al. Genome-wide analysis of simple sequence repeats in bitter gourd (Momordica charantia)[J]. Frontiers in Plant Science, 2017, 8: 1 103.

    猜你喜歡
    三七多態(tài)性基因組
    單核苷酸多態(tài)性與中醫(yī)證候相關(guān)性研究進(jìn)展
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    人參屬藥用植物轉(zhuǎn)錄組研究進(jìn)展
    超高液相色譜法測定三七中三七總皂苷的方法研究
    我所知道的“三七”
    馬鈴薯cpDNA/mtDNA多態(tài)性的多重PCR檢測
    GlobalFiler~? PCR擴(kuò)增試劑盒驗(yàn)證及其STR遺傳多態(tài)性
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    CYP3A4*1G基因多態(tài)性及功能的初步探討
    一级毛片电影观看 | 午夜福利在线在线| 变态另类成人亚洲欧美熟女| 精品久久久久久久人妻蜜臀av| 国产v大片淫在线免费观看| 亚洲成人久久爱视频| 国产精品久久视频播放| 黄色视频,在线免费观看| 亚州av有码| 国产一区二区三区在线臀色熟女| 99久久中文字幕三级久久日本| 国产精品无大码| 午夜福利高清视频| 国产白丝娇喘喷水9色精品| 中文字幕av在线有码专区| 成人美女网站在线观看视频| 日本三级黄在线观看| 亚洲一区高清亚洲精品| 国产精品嫩草影院av在线观看| 中国美女看黄片| 校园人妻丝袜中文字幕| 淫秽高清视频在线观看| 一级毛片电影观看 | 午夜老司机福利剧场| 一个人免费在线观看电影| 97在线视频观看| 国产国拍精品亚洲av在线观看| 日日撸夜夜添| 久久久久久伊人网av| 国产一区二区在线观看日韩| 免费看a级黄色片| 国产成人精品久久久久久| 波野结衣二区三区在线| 又黄又爽又免费观看的视频| 国产精品99久久久久久久久| 看黄色毛片网站| 国产欧美日韩精品一区二区| 一级毛片aaaaaa免费看小| 最近最新中文字幕大全电影3| 久久精品国产亚洲av天美| 男人狂女人下面高潮的视频| 日本黄色片子视频| 最近2019中文字幕mv第一页| 久久午夜亚洲精品久久| 插阴视频在线观看视频| 欧美xxxx黑人xx丫x性爽| 91久久精品国产一区二区成人| 一区二区三区四区激情视频 | 一级毛片电影观看 | 可以在线观看的亚洲视频| 51国产日韩欧美| 国产淫片久久久久久久久| 日韩亚洲欧美综合| 国产精品嫩草影院av在线观看| 麻豆成人午夜福利视频| 少妇熟女欧美另类| 亚洲人与动物交配视频| av女优亚洲男人天堂| 国产片特级美女逼逼视频| 国产欧美日韩一区二区精品| 别揉我奶头~嗯~啊~动态视频| 可以在线观看的亚洲视频| 人妻少妇偷人精品九色| 亚洲欧美日韩高清专用| 久久午夜亚洲精品久久| 午夜爱爱视频在线播放| 99热6这里只有精品| 国产亚洲欧美98| a级毛色黄片| 久久精品国产99精品国产亚洲性色| 中出人妻视频一区二区| 欧美成人免费av一区二区三区| 特大巨黑吊av在线直播| 久久草成人影院| 国产v大片淫在线免费观看| 99热这里只有精品一区| 亚洲精品乱码久久久v下载方式| 欧美另类亚洲清纯唯美| 国内久久婷婷六月综合欲色啪| 国内久久婷婷六月综合欲色啪| 一本精品99久久精品77| 久久久a久久爽久久v久久| 十八禁国产超污无遮挡网站| 黄色欧美视频在线观看| 精品久久久久久久久久免费视频| 久久久精品大字幕| 亚洲欧美日韩高清在线视频| 国产色婷婷99| 真人做人爱边吃奶动态| 久久久久免费精品人妻一区二区| 国产精品乱码一区二三区的特点| 日本熟妇午夜| 久久久久国产精品人妻aⅴ院| 亚洲成人av在线免费| 欧美日韩一区二区视频在线观看视频在线 | 免费一级毛片在线播放高清视频| 美女被艹到高潮喷水动态| 美女被艹到高潮喷水动态| 午夜视频国产福利| 免费高清视频大片| 国产高清激情床上av| 国产高潮美女av| 成人亚洲欧美一区二区av| 婷婷精品国产亚洲av| 在线观看免费视频日本深夜| 免费观看在线日韩| 久久精品人妻少妇| 日本与韩国留学比较| 成人三级黄色视频| 国产 一区精品| 久久中文看片网| 国产成人福利小说| 俄罗斯特黄特色一大片| 少妇裸体淫交视频免费看高清| 免费观看人在逋| 卡戴珊不雅视频在线播放| 亚洲精品国产av成人精品 | .国产精品久久| 高清午夜精品一区二区三区 | 国产精品久久久久久久电影| 久久热精品热| 精品久久久久久久久亚洲| 亚洲av.av天堂| 露出奶头的视频| 狠狠狠狠99中文字幕| 男插女下体视频免费在线播放| 免费观看的影片在线观看| 老司机影院成人| 韩国av在线不卡| 精品日产1卡2卡| 老师上课跳d突然被开到最大视频| 国内精品久久久久精免费| 麻豆一二三区av精品| 亚洲无线观看免费| 黄色配什么色好看| 国产精品亚洲一级av第二区| 三级国产精品欧美在线观看| 午夜福利在线在线| 色噜噜av男人的天堂激情| 久久久久久伊人网av| 国产黄a三级三级三级人| 俄罗斯特黄特色一大片| 欧美成人一区二区免费高清观看| 日本欧美国产在线视频| 午夜爱爱视频在线播放| 精品日产1卡2卡| 亚洲欧美日韩高清在线视频| 亚洲无线观看免费| 国产精品1区2区在线观看.| 一区二区三区高清视频在线| 午夜福利在线在线| 亚洲美女视频黄频| 深夜精品福利| 91久久精品国产一区二区三区| 综合色丁香网| 国产成人影院久久av| 午夜爱爱视频在线播放| 国产精品久久久久久久久免| 亚洲第一电影网av| 天堂av国产一区二区熟女人妻| 欧美激情国产日韩精品一区| 亚洲av二区三区四区| 免费在线观看影片大全网站| 午夜福利18| 在线免费观看的www视频| 国产一区二区亚洲精品在线观看| 久久综合国产亚洲精品| 最近手机中文字幕大全| 成人鲁丝片一二三区免费| 日韩欧美国产在线观看| 中国美白少妇内射xxxbb| h日本视频在线播放| 午夜福利在线观看免费完整高清在 | 午夜福利成人在线免费观看| 真实男女啪啪啪动态图| 搡老熟女国产l中国老女人| 你懂的网址亚洲精品在线观看 | 国产黄色视频一区二区在线观看 | 波多野结衣高清无吗| 精品久久久久久久久亚洲| 九色成人免费人妻av| 日本免费a在线| 最后的刺客免费高清国语| 国产视频内射| 日本一本二区三区精品| 国产亚洲精品av在线| 成人亚洲欧美一区二区av| 日本与韩国留学比较| 高清日韩中文字幕在线| 黄色视频,在线免费观看| 国产精华一区二区三区| 99热6这里只有精品| 精品久久久久久久久久久久久| 看黄色毛片网站| 熟女电影av网| 人人妻,人人澡人人爽秒播| 久久久成人免费电影| 舔av片在线| videossex国产| 99热这里只有是精品在线观看| 国产大屁股一区二区在线视频| 婷婷精品国产亚洲av| 在线观看av片永久免费下载| 少妇被粗大猛烈的视频| 精品久久国产蜜桃| 18禁裸乳无遮挡免费网站照片| 啦啦啦啦在线视频资源| 天堂网av新在线| 人人妻人人澡人人爽人人夜夜 | 国产在线男女| 小说图片视频综合网站| 此物有八面人人有两片| 深夜精品福利| 久久天躁狠狠躁夜夜2o2o| 国产精品美女特级片免费视频播放器| 久久精品国产鲁丝片午夜精品| 中文资源天堂在线| 亚洲第一区二区三区不卡| 国产精品无大码| 国产成人福利小说| 国产精品久久视频播放| 精品午夜福利视频在线观看一区| 亚洲av中文字字幕乱码综合| 国产高清不卡午夜福利| 一级a爱片免费观看的视频| 九色成人免费人妻av| 精品久久久久久久末码| 97超级碰碰碰精品色视频在线观看| 99热精品在线国产| 日韩精品中文字幕看吧| 国产精品福利在线免费观看| 一a级毛片在线观看| 免费电影在线观看免费观看| 我要搜黄色片| 精品人妻熟女av久视频| 国产av一区在线观看免费| 国产成人一区二区在线| 久久热精品热| 亚洲一级一片aⅴ在线观看| 一区二区三区四区激情视频 | 欧美色视频一区免费| 级片在线观看| 国产精品亚洲美女久久久| 一进一出抽搐gif免费好疼| 丰满的人妻完整版| 亚洲av第一区精品v没综合| 午夜爱爱视频在线播放| 亚洲天堂国产精品一区在线| 国产探花极品一区二区| 日韩人妻高清精品专区| 国产欧美日韩精品亚洲av| 男女视频在线观看网站免费| 老司机午夜福利在线观看视频| 免费av观看视频| 中文字幕av成人在线电影| 日韩欧美三级三区| 18禁在线无遮挡免费观看视频 | 国产精品久久久久久亚洲av鲁大| 国产一区二区激情短视频| 国产精华一区二区三区| 日本在线视频免费播放| 日韩国内少妇激情av| 亚洲激情五月婷婷啪啪| 有码 亚洲区| 国产精品嫩草影院av在线观看| 日本三级黄在线观看| 色综合色国产| 色在线成人网| 久久久久久大精品| 久99久视频精品免费| 免费人成在线观看视频色| 国产精品亚洲一级av第二区| 久久精品夜夜夜夜夜久久蜜豆| 美女xxoo啪啪120秒动态图| 成人av在线播放网站| 黄色欧美视频在线观看| 成人特级黄色片久久久久久久| 九九久久精品国产亚洲av麻豆| а√天堂www在线а√下载| 免费在线观看成人毛片| 99精品在免费线老司机午夜| 免费在线观看影片大全网站| 国产精品一区二区三区四区久久| 91久久精品国产一区二区成人| 性欧美人与动物交配| 久久亚洲国产成人精品v| 午夜福利在线在线| 久久久久久国产a免费观看| 国产亚洲精品久久久com| 国产 一区 欧美 日韩| .国产精品久久| 草草在线视频免费看| 亚洲av第一区精品v没综合| 91精品国产九色| 内射极品少妇av片p| av国产免费在线观看| 日韩成人av中文字幕在线观看 | ponron亚洲| 床上黄色一级片| 麻豆国产av国片精品| 大又大粗又爽又黄少妇毛片口| 身体一侧抽搐| 国产男人的电影天堂91| 老司机影院成人| 免费一级毛片在线播放高清视频| 国产高清视频在线播放一区| 欧美精品国产亚洲| 简卡轻食公司| 精品一区二区三区av网在线观看| 在线观看66精品国产| 久久精品久久久久久噜噜老黄 | 不卡视频在线观看欧美| 国产伦一二天堂av在线观看| 此物有八面人人有两片| 激情 狠狠 欧美| 中文字幕精品亚洲无线码一区| 99久久精品一区二区三区| 三级国产精品欧美在线观看| 国产国拍精品亚洲av在线观看| 午夜爱爱视频在线播放| 日韩一区二区视频免费看| 国产午夜福利久久久久久| 国产高清视频在线观看网站| 99在线视频只有这里精品首页| 一级毛片aaaaaa免费看小| 男女啪啪激烈高潮av片| or卡值多少钱| 久久久久性生活片| 亚洲无线观看免费| 国产一区二区三区av在线 | 美女xxoo啪啪120秒动态图| 久久99热6这里只有精品| 成人永久免费在线观看视频| 国产一区二区三区av在线 | 嫩草影院精品99| 在线免费观看不下载黄p国产| 小蜜桃在线观看免费完整版高清| 人人妻,人人澡人人爽秒播| 一级毛片久久久久久久久女| 中出人妻视频一区二区| 亚洲av五月六月丁香网| 99久久精品热视频| 免费人妻精品一区二区三区视频| 夫妻午夜视频| 在线 av 中文字幕| 久久久精品免费免费高清| 男女啪啪激烈高潮av片| 亚洲精品国产色婷婷电影| a 毛片基地| 亚洲无线观看免费| 亚洲av国产av综合av卡| 18禁在线播放成人免费| 欧美成人精品欧美一级黄| 看十八女毛片水多多多| 亚洲综合色惰| 日韩亚洲欧美综合| 午夜免费观看性视频| 在线观看免费日韩欧美大片 | 日韩制服骚丝袜av| 人妻人人澡人人爽人人| 日韩精品免费视频一区二区三区 | 亚洲精品国产色婷婷电影| 丰满饥渴人妻一区二区三| 亚洲中文av在线| 一级毛片 在线播放| 一本一本综合久久| av视频免费观看在线观看| 亚洲天堂av无毛| av卡一久久| 日韩电影二区| 国产极品天堂在线| 观看美女的网站| 日日啪夜夜爽| 天天躁夜夜躁狠狠久久av| 中国国产av一级| 在线亚洲精品国产二区图片欧美 | 国产一区有黄有色的免费视频| 中文字幕人妻丝袜制服| 国产亚洲欧美精品永久| 韩国av在线不卡| 亚洲性久久影院| 国产极品粉嫩免费观看在线 | 国产极品粉嫩免费观看在线 | 国产一区二区在线观看av| 国产精品熟女久久久久浪| 亚洲人成网站在线观看播放| 亚洲美女黄色视频免费看| 嫩草影院入口| 青春草视频在线免费观看| 亚洲欧美精品专区久久| 男男h啪啪无遮挡| 丝袜喷水一区| 这个男人来自地球电影免费观看 | 精品酒店卫生间| 日本wwww免费看| 国产精品熟女久久久久浪| 人人妻人人看人人澡| 久久精品国产自在天天线| 在现免费观看毛片| 国产毛片在线视频| 久久久久久久久大av| 久久久久精品性色| 亚洲国产精品999| 亚洲国产精品专区欧美| 丝袜脚勾引网站| 看免费成人av毛片| 中国三级夫妇交换| 亚洲精品成人av观看孕妇| 三级经典国产精品| 多毛熟女@视频| 18禁动态无遮挡网站| 你懂的网址亚洲精品在线观看| 搡老乐熟女国产| 亚洲精品乱码久久久v下载方式| 午夜福利在线观看免费完整高清在| 精品久久久久久久久av| 91精品国产九色| 天堂中文最新版在线下载| 欧美性感艳星| 老熟女久久久| 精品熟女少妇av免费看| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品第二区| av天堂中文字幕网| 夫妻午夜视频| 亚洲国产成人一精品久久久| 天堂中文最新版在线下载| 天美传媒精品一区二区| 精品亚洲成国产av| 亚洲精品,欧美精品| 国语对白做爰xxxⅹ性视频网站| 一级毛片aaaaaa免费看小| 69精品国产乱码久久久| 曰老女人黄片| 内射极品少妇av片p| 各种免费的搞黄视频| 热re99久久精品国产66热6| 在线天堂最新版资源| 偷拍熟女少妇极品色| 一级爰片在线观看| 观看av在线不卡| 亚洲综合精品二区| a级片在线免费高清观看视频| 九九久久精品国产亚洲av麻豆| 97精品久久久久久久久久精品| 久久女婷五月综合色啪小说| 免费av中文字幕在线| 久久精品国产亚洲av天美| 边亲边吃奶的免费视频| 在线观看免费高清a一片| 80岁老熟妇乱子伦牲交| 一区二区av电影网| 一级黄片播放器| 国产黄片美女视频| av天堂久久9| 少妇精品久久久久久久| 亚洲精品日韩av片在线观看| 麻豆成人午夜福利视频| 日韩熟女老妇一区二区性免费视频| 国产精品无大码| 啦啦啦视频在线资源免费观看| 午夜免费观看性视频| 搡女人真爽免费视频火全软件| 少妇熟女欧美另类| 最近中文字幕2019免费版| 亚洲欧洲精品一区二区精品久久久 | 中文字幕亚洲精品专区| 男女免费视频国产| 亚洲色图综合在线观看| 成年女人在线观看亚洲视频| 亚洲欧美成人综合另类久久久| 国产视频内射| 天堂中文最新版在线下载| 中文字幕av电影在线播放| 国产色婷婷99| 九九久久精品国产亚洲av麻豆| 久久精品国产亚洲网站| 天堂中文最新版在线下载| 久久99热这里只频精品6学生| 两个人免费观看高清视频 | 亚洲精品国产av蜜桃| 国产永久视频网站| 一区二区三区乱码不卡18| 亚洲精华国产精华液的使用体验| 丰满人妻一区二区三区视频av| 久热久热在线精品观看| 亚洲精品视频女| 国产精品熟女久久久久浪| 欧美最新免费一区二区三区| 观看美女的网站| 国产在线男女| 自拍欧美九色日韩亚洲蝌蚪91 | 波野结衣二区三区在线| 国产日韩一区二区三区精品不卡 | 日本欧美国产在线视频| 亚洲中文av在线| 人人澡人人妻人| a级一级毛片免费在线观看| 国产成人免费观看mmmm| 欧美 日韩 精品 国产| av在线观看视频网站免费| 国产精品麻豆人妻色哟哟久久| 婷婷色综合大香蕉| 国产精品一二三区在线看| 欧美精品一区二区大全| 亚洲综合精品二区| www.av在线官网国产| 国产日韩一区二区三区精品不卡 | 中文字幕精品免费在线观看视频 | 欧美性感艳星| 两个人免费观看高清视频 | 亚洲精品国产av蜜桃| 国产精品成人在线| 国产精品熟女久久久久浪| 人体艺术视频欧美日本| 日日爽夜夜爽网站| 欧美老熟妇乱子伦牲交| 一级毛片电影观看| 免费少妇av软件| 国产成人免费观看mmmm| 亚洲人成网站在线播| 五月伊人婷婷丁香| 亚洲精品成人av观看孕妇| 在线看a的网站| 一本久久精品| 看十八女毛片水多多多| 欧美日韩视频高清一区二区三区二| 免费看不卡的av| 精品久久国产蜜桃| 中文字幕久久专区| 黑人高潮一二区| av在线app专区| 日韩制服骚丝袜av| 久久久久久久大尺度免费视频| 日本与韩国留学比较| 麻豆精品久久久久久蜜桃| 中文天堂在线官网| 日本欧美国产在线视频| av福利片在线观看| 久久6这里有精品| 在线精品无人区一区二区三| 成人影院久久| 国产av一区二区精品久久| av又黄又爽大尺度在线免费看| 中文在线观看免费www的网站| 日韩伦理黄色片| 免费看日本二区| 深夜a级毛片| 啦啦啦啦在线视频资源| 亚洲精品久久久久久婷婷小说| 简卡轻食公司| 美女大奶头黄色视频| 亚洲精品日本国产第一区| 欧美日韩一区二区视频在线观看视频在线| 热re99久久国产66热| 成人亚洲精品一区在线观看| 国产av国产精品国产| 亚洲国产精品999| 我要看日韩黄色一级片| 97超碰精品成人国产| 日韩强制内射视频| 亚洲人成网站在线观看播放| 精品久久久久久久久av| 久久国产精品男人的天堂亚洲 | 狠狠精品人妻久久久久久综合| 日本av免费视频播放| 国产精品一二三区在线看| 久久99精品国语久久久| 日韩一区二区三区影片| 日产精品乱码卡一卡2卡三| 亚洲精品成人av观看孕妇| 伊人亚洲综合成人网| 亚洲国产精品成人久久小说| 交换朋友夫妻互换小说| 国产精品久久久久久精品电影小说| 在线看a的网站| 在现免费观看毛片| 成年美女黄网站色视频大全免费 | 美女福利国产在线| av有码第一页| 国精品久久久久久国模美| 久久久久国产网址| 久久 成人 亚洲| h日本视频在线播放| videossex国产| 日本免费在线观看一区| 亚洲四区av| 两个人的视频大全免费| 丝袜脚勾引网站| 亚洲精品乱码久久久久久按摩| 久久人人爽人人片av| 免费高清在线观看视频在线观看| 国产精品一二三区在线看| 国产精品偷伦视频观看了| 国产日韩欧美视频二区| 91aial.com中文字幕在线观看| 色视频www国产| 亚洲天堂av无毛| 三上悠亚av全集在线观看 | 免费看av在线观看网站| 日本午夜av视频| 欧美日韩亚洲高清精品| 少妇被粗大猛烈的视频| 中文字幕亚洲精品专区| 日韩不卡一区二区三区视频在线| 内地一区二区视频在线| 精品久久久久久电影网| 在线观看人妻少妇| 亚洲av福利一区| 日韩一区二区视频免费看| a级毛片免费高清观看在线播放| 久久久欧美国产精品| 一级毛片久久久久久久久女| 伊人久久精品亚洲午夜| 美女福利国产在线|