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

    單葉省藤和黃藤基因組組裝研究

    2018-03-17 05:37:00王思寧趙韓生高志民
    世界竹藤通訊 2018年6期
    關鍵詞:單葉棕櫚基因組

    王思寧 趙韓生 高志民

    (國際竹藤中心 北京 100102 )

    棕櫚藤是世界上最重要的非木材森林植物資源之一,屬于棕櫚科(Arecaceae),是自然進化中出現(xiàn)的攀爬植物的代表。近年來的研究表明,在棕櫚科省藤族省藤亞科中棕櫚藤植物可分為11屬,共有631種,其中省藤屬(NCBITaxonID:4711)和黃藤屬(NCBITaxonID:93268)種類最多,分別占藤種的約65%和20%。這2個屬也是最重要的棕櫚藤藤材來源,為工業(yè)生產提供了95%以上的藤條。全球有超過500萬人在經濟上依賴棕櫚藤,藤工業(yè)產值每年約70億美元,包括國內工業(yè)生產、國際貿易,以及藤編織的籃子、座椅和家具等。同時,對棕櫚藤遺傳育種的關注在持續(xù)增長,并且預計將來在幾年內棕櫚藤的種植面積將逐漸超過天然棕櫚藤面積。

    單葉省藤(Calamussimplicifolius)(NCBITaxonID:746888)是中國原產、種植面積最廣的藤種之一,莖長約50m,直徑約15mm。單葉省藤是海南島特有的棕櫚藤,可生產高品質、中等直徑的藤條,用于捆扎和編織。黃藤(Daemonoropsjenkinsiana)(NCBITaxonID:1510057),是常綠攀援棕櫚藤的代表,自然生長在海拔1000m以下的低地雨林,分布于孟加拉國、不丹、柬埔寨、印度、老撾、緬甸、尼泊爾、泰國、越南和中國的東南部。黃藤叢生生長,莖可達50m, 節(jié)間長40cm,直徑30mm。這2種產值最高的藤種現(xiàn)種植面積超過1000hm2,栽培在中國小于北緯23°30′的地區(qū),例如海南島、廣東、廣西、云南、福建等華南地區(qū)。

    單葉省藤和黃藤由于其藤條韌性強和耐久度好而具有廣泛的用途和巨大的開發(fā)潛力。需要應用分子育種技術開展棕櫚藤育種,來滿足對棕櫚藤質量和數(shù)量不斷增長的需求。然而,棕櫚藤遺傳背景不清,嚴重阻礙了對其分子生物學的研究和實際生產的應用,更影響了相關物種之間比較基因組分析的深入開展。因此,我們使用最新的高通量測序技術,完成了單葉省藤和黃藤的測序和染色體水平基因組組裝。隨著這2個棕櫚藤基因組的完成,使許多比較基因組分析和其他下游分析成為可能,包括分子標記的開發(fā)、功能基因的鑒定和分子輔助育種等。此外,高質量的基因組將有助于對棕櫚藤生物性狀進行基因組學、轉錄組學和代謝組學分析。本文選取并鑒定了棕櫚藤中改良藤條性能的一些重要候選基因——木質素生物合成基因家族,為今后利用木質素生物合成基因來改良藤條的品質以及研究物種多樣性奠定了基礎。

    1 材料方法與數(shù)據(jù)處理

    1.1 DNA分離、文庫構建和高通量測序

    單葉省藤和黃藤材料取自中國林業(yè)科學研究院熱帶林業(yè)研究所(北緯23°11′29″,東經113°22′40″,海拔87m),于2015年春季采集處于營養(yǎng)生長階段的幼嫩葉片組織用于總DNA提取。使用DNeasyPlantMiniKits(Qiagen)分離和提取總DNA。根據(jù)高分子量核DNA分離法純化基因組DNA獲得多個DNA文庫,并在IlluminaHiSeq4000和PacBioSequel平臺上進行測序(原文表1)。即:在基于Illumina的方法上構建了3個不同插入片段大小的文庫(270bp、500bp和800bp),用于PE測序,以及4個不同插入片段大小的文庫(2kb、5kb、10kb和20kb)用于MP測序。另外,按照PacBio方法,構建了5個插入片段大小為20kb的PacBioSequel文庫,進行測序。通過對測序數(shù)據(jù)清理和預處理,最終得到約494.08Gb的單葉省藤高質量數(shù)據(jù)(322.3GbPEreads,93.4GbMPreads,78.38GbPacBiodata),覆蓋率為單葉省藤基因組的252倍, 約426.17Gb的黃藤高質量數(shù)據(jù)(244.58GbPEreads,103.21GbMPreads,and78.38GbPacBiodata),覆蓋率為黃藤基因組的266倍。

    作為與Illumina和PacBio文庫構建測序,進行了另一平行分析,即使用相同組織材料為單葉省藤和黃藤建立了2個高通量測序文庫。使用MboI限制性內切酶在甲醛固定后消化DNA,然后用生物素標記5個殘基。在原位連接平末端片段后,對分離出的DNA進行反向交聯(lián)、純化和過濾,得到包含生物素標記的片段。隨后,按此順序進行DNA片段端修復、適配體連接和聚合酶鏈式反應。然后,根據(jù)BGISEQ-500的使用說明,對其進行100PE的測序。在分別獲得了大約148Gb和154Gb的原始數(shù)據(jù)后,經過低質量數(shù)據(jù)的過濾和其他預處理,最終獲得約6.7Gb和13.1Gb的高質量數(shù)據(jù),并使用HiC-Pro(version2.8.0devel)軟件分別對單葉省藤和黃藤進行了評價和分析(原文表1)。

    1.2 基因組調查

    對給定物種基因組大小、雜合性等基因組特征的了解解,將有助于開發(fā)一套適合的測序方案和組裝策略。因此,本研究采用4種獨立的方法來估算基因組大小,即KmerSpectrumPlot.plinALLPATHS-LG(版本r52488)、GCE(GenomeCharacteristicsEstimation,released7Jan.2015)、JELLYFISH(版本2.0)和流式細胞術(原文表S1、S2,圖S1、S2)。在調研中,我們測序獲得單葉省藤約為98Gb,黃藤約為60Gb的二代原始測序數(shù)據(jù)作為調研數(shù)據(jù)。經過數(shù)據(jù)預處理和組裝,結果顯示(原文表S1),最終預測單葉省藤和黃藤基因組大小分別約為1.98Gb和1.61Gb,相關雜合度分別為1.32%~1.52%和1.19%~1.31%。因此,基因組分析結果表明,這2種棕櫚藤的基因組比較適用于Illumina和PacBio測序數(shù)據(jù)的聯(lián)合組裝。

    1.3 利用Illumina和PacBio高通量測序數(shù)據(jù)聯(lián)合組裝基因組

    在Illumina數(shù)據(jù)的預處理過程中,過濾掉了低質量的讀長和接頭序列,獲得單葉省藤和黃藤的有效數(shù)據(jù)分別為416Gb和348Gb。對于PacBio數(shù)據(jù),使用了MECAT(發(fā)布于2017年6月27日)軟件,得到單葉省藤和黃藤分別為52Gb和32Gb的高質量數(shù)據(jù)。隨后,使用FALCON(版本0.3)軟件進行第1次contig組裝。如原文表S3所示,使用不同參數(shù)為單葉省藤基因組生成了2個組裝結果,一個ContigN50為67.2kb,基因組約為1.59Gb(約為預測基因組大小的80%),另一個ContigN50為66.7kb,基因組約為1.53Gb(約為預測基因組大小的77%)。此外,獲得一個黃藤組裝結果ContigN50為81.5kb,基因組約為1.27Gb(約為預測基因組大小的79%)。通過以上結果顯示,僅用三代數(shù)據(jù)和MECAT軟件不能獲得對2種高質量棕櫚藤基因組,可能是由于2種棕櫚藤的高雜度(單葉省藤1.32%~1.52%,黃藤1.19%~1.31%)、高比例重復序列(單葉省藤為54.15%,黃藤為70%;有關詳細信息請參閱原文中的后續(xù)分析)和使用了測序深度不足的三代數(shù)據(jù)(PacBio數(shù)據(jù)經錯誤校正后分別為26×和20×)造成的。因此結合以上結果,我們嘗試利用Illumina和PacBio測序數(shù)據(jù),對單葉省藤和黃藤進行了聯(lián)合基因組組裝。首先,利用Platanus(版本1.2.4,是一種用于高雜合數(shù)據(jù)的全新基因組匯編器,可以通過構建具有自動優(yōu)化的k-mer大小的DeBruijn圖)將PE片段整合到contig序列中。其次,用DBG2OLC(2015年7月11日發(fā)布,參照以下參數(shù):DBG2OLCContigscontig.faLD0K17KmerCovTh4MinOverlap25AdaptiveTh0.007RemoveChimera1fscaffold.fa.) 矯正PacBio測序數(shù)據(jù)并組裝scaffolds,以分別獲得單葉省藤和黃藤約1.92Gb和1.56Gb的初始組裝結果。再次,在進行SSPACE分析前,參照DBG2OLC(原文表S4)進行了平行分析,這一步驟有助于提高基因組組裝的質量,減少SSPACE過程中的錯誤。然后,使用SSPACE(版本3.0)和MP大片度測序數(shù)據(jù)延長scaffolds,并通過使用GapCloser(版本1.12)和PBJelly(2015年8月24日發(fā)布) 進行補洞。最終,獲得了1.96Gb的單葉省藤組裝數(shù)據(jù),包含5116個Contigs,ContigN50長度為107kb,ScaffoldN50長度為803kb,同時獲得了一個1.60Gb的黃藤組裝數(shù)據(jù),包含ContigN50和ScaffoldN50長度分別為108kb和784kb(原文表2)。

    隨后,通過3D-DNA分析流程(版本170123)對有效的高通量測序數(shù)據(jù)與上述裝配數(shù)據(jù)一起進行處理,得到了一個合理的連接模式,準確的染色體組裝結果。如原文圖2所示,使用Juicebox(版本1.5.2)顯示連接模式圖譜。原文表S5展示了單葉省藤12個最長的染色體水平Scaffold和黃藤13個最長的染色體水平Scaffold。預測染色體總長度分別占單葉省藤和黃藤基因組總長度的92.08%和92.01%,ScaffoldN50分別為169Mb和119Mb。

    1.4 基因組評估

    采用3種獨立的方法評估了單葉省藤和黃藤基因組組裝的準確性和完整性。首先,對基因組中不確定堿基的百分比(Ns)和鳥嘌呤、胞嘧啶含量(GC)進行了分析。結果顯示在每個基因組中低豐度的Ns(單葉省藤約為0.6%和黃藤約為0.7%)和總GC含量(單葉省藤為41.78%和黃藤為41.07%)與相關轉錄組數(shù)據(jù)相似(單葉省藤為41.68%和黃藤為41.89%)。此外,利用BasicLocalAlignmentSearchTool(BLAST)-like比對工具(版本1.0)(默認參數(shù))進行RNA測序(RNAseq)數(shù)據(jù)所組裝出來的功能基因序列與組裝得到的基因組序列比對,結果顯示,在一個大片段上90%以上的序列可以與基因組的組裝序列匹配(單葉省藤為92.89%和黃藤為81.81%)(原文表S6)。最后,利用BUSCO(版本3.0)對2種棕櫚藤基因組的完整性進行了評估,這是一種從進化的角度定量對單拷貝直系同源基因含量進行計算,以評估基因組完整性的方法。BUSCO結果顯示,在單葉省藤基因組組裝結果中檢測到保守性為96.4%的BUSCO蛋白(基于embryophytaodb9庫文件)(包括3.8%的BUSCO蛋白片段)。同時,在黃藤基因組組裝結果中檢測到保守性為87.3%的BUSCO蛋白(包括4.0%的BUSCO蛋白片段)(原文表S7)。

    1.5 重復注釋

    本研究在蛋白質編碼基因模型預測之前,完成了單葉省藤和黃藤基因組中轉座因子(TEs)和串聯(lián)重復序列的鑒定。采用2種獨立的方法來預測序列重復:同源注釋和從頭預測(denovo)法。同源注釋結果顯示,利用RepeatMasker(版本4.0.5)和RepeatProteinMasker(版本4.0.5)對Repbase文庫(發(fā)布于2017年12月01日)進行搜索,查找轉座因子序列。從頭預測結果顯示,利用Repeat-Modeler(版本1.0.8)和LTRFINDER方法在去除污染和多拷貝基因后,可以構建一個重復序列文庫。然后,使用RepeatMasker對基因組序列進行了分類,并與重復序列文庫進行了對比。此外,串聯(lián)重復序列由TandemRepeatFinder(版本4.09)識別,其參數(shù)為:Match=2,Mismatch=7,Delta=7,PM=80,PI=10,Minscore=50andMaxPeriod=2000。最終結果表明,長末端重復序列(longterminalrepeat,LTR)是最豐富的重復序列類型,2種非LTR逆轉座子(短穿插核元件和長穿插核元件)在2種棕櫚藤基因組中占比例最少(原文表S8)。TEs在單葉省藤和黃藤基因組中分別占54.15%和70%。TEs序列分歧度結果表明,從頭預測的重復序列比Repbase預測的重復序列更多(原文圖3)。

    1.6 RNA樣本收集、文庫構建和轉錄組組裝

    分別選取單葉省藤和黃藤3個不同發(fā)育階段的纖鞭樣本,每個樣本3個生物學重復(原文表S9)。為保證與2種棕櫚藤基因組研究的一致性,RNA取樣的位置和DNA取樣的位置一致。利用RNA提取試劑盒(TaKaRa, 日本)提取各樣本的總RNA,并使用NanoDrop2000分光光度計測定RNA的質量和濃度。使用RecombinantDNaseI于37℃孵育30min去除殘余DNA,并使用反轉錄試劑盒(Promega, 美國)合成cDNA第一條鏈,操作參照相關試劑盒說明書進行。然后,使用BGISEQ-500平臺對文庫池中的短100PE片段測序。在對轉錄組數(shù)據(jù)進行預處理時,使用SOAPnuke(版本1.5.6)對匹配序列和低質量的讀長序列進行過濾,過濾參數(shù)為:-n0.001-l20-q0.4-q2。所有樣本均使用Trinity(版本2.0.6)按以下參數(shù)組裝:(1)grouppairsdistance500, (2)mincontiglength200, (3)minkmercov2, (4)minglue2, (5)bflyopts-V5, (6)edge-thr=0.1, (7)stderr,and(8)SSlibtypeRF。然后,對Trinity的輸出結果進行聚類分析,使用TGI聚類工具(版本v2.0.6)生成一組非冗余集合,參數(shù)如下:(1)aminimumof95%identitybetweenthecontigs, (2)aminimumof35overlappingbases, (3)aminimumscoreof35,and(4)amaximumof20unmatchedoverhangingbasesatthesequenceends。最后,根據(jù)序列相似度將整合的轉錄本分為2類:成簇轉錄本和單一轉錄本。在成簇轉錄本中,轉錄本之間的序列相似區(qū)都超過了70%,由這些轉錄本被剪接成同源基因或旁系同源基因。此外,所有鑒定的基因將被用于后續(xù)的功能分析。

    1.7 基因建模與預測

    使用3種獨立的方法對完整的蛋白質編碼基因模型進行了綜合預測,包括從頭預測、同源蛋白預測和轉錄組數(shù)據(jù)預測。使用AUGUSTUS(版本3.3)的默認參數(shù)對重復區(qū)段進行第一次注釋,該程序是一個基于自我調試模型的從頭預測軟件。通過進行多次調試優(yōu)化數(shù)據(jù),在單葉省藤和黃藤中分別預測獲得了85246和87613個基因模型。在基于同源性的預測中,以油棕(Elaeisguineensis)、長葉刺葵(Phoenixdactylifera)、二穗短柄草(Brachypodiumdistachyon)、水稻(Oryzasativa)、粟(Setariaitalica)、高粱(Sorghumbicolor)和玉米(Zeamays)等7個物種作為參考基因組。從ENSEMBL數(shù)據(jù)庫下載它們的蛋白序列并應用TBLASTN(版本2.2.26)對單葉省藤和黃藤中E-value不小于1e-5的組裝序列進行比對。然后,由GeneWise(版本2.0)生成拼接結果。轉錄組數(shù)據(jù)預測結果顯示,使用HISAT2(版本2.0.2)來識別外顯子-內含子剪接位置,并進一步進行轉錄組序列與基因組序列的比對。使用Cufflinks(版本2.2.1)分別在單葉省藤和黃藤中分別鑒定了56024和58134個蛋白質編碼基因模型(原文表S10)。最后,使用MAKER(版本2)整合了上述3個獨立預測的數(shù)據(jù)結果,隨著相同基因序列在單葉省藤和黃藤基因組中的出現(xiàn),最終共分別鑒定出51235和53342個完整的蛋白編碼基因模型。

    1.8 基因注釋和基因功能預測

    使用2種獨立的方法對預測的基因集合進行評估:基因功能評估和BUSCO完整性評估。在基因功能評估中,比照基因在近緣物種中同源蛋白序列和預測的蛋白集合,評估了預測的注釋和蛋白質序列的一致性。參照NCBI非冗余蛋白質數(shù)據(jù)庫(2018年3月13日發(fā)布)、SWISS-PROT(2018年1月1日發(fā)布)、GeneOntology(GO)(2013年10月30日發(fā)布)、KyotoEncyclopediaofGenesandGenomes(KEGG)(數(shù)據(jù)集v81)和InterPro(數(shù)據(jù)集v.53)等5個權威蛋白數(shù)據(jù)庫進行注釋比對,結果表明(原文表S11)在單葉省藤和黃藤預測的基因模型中分別有5.34%和2.89%的基因為未注釋的基因。此外,BUSCO的評估顯示,在單葉省藤和黃藤中分別有88.7%和91.3%的保守BUSCO蛋白(embryophytaodb9)存在。在這些保守的BUSCO蛋白中,分別有76.2%和81.2%是完整的。此外,tRNA,rRNA,miRNA和snRNA4種非編碼RNA也被預測出來(原文表S12)。

    1.9 基因家族結構和棕櫚藤特異性基因家族分析

    本研究中,通過逐一序列比對的方法來預測基因組水平的直系同源基因。這種方法因速度快而常用于處理大數(shù)據(jù)量分析。采用BLAST-based的OrthoMCL(版本2.0.9)來識別單葉省藤和黃藤中E-value≤1e-5且匹配系數(shù)≥80的同源基因。Markovchain聚類(默認膨脹參數(shù))也共同用于與其他8個物種的全面BLASTP分析,包括無油樟(Amborellatrichopoda)、油棕、擬南芥、二穗短柄草、水稻、紫萍(Spirodelapolyrhiza)、長葉刺葵和高粱。在所有10個物種中鑒定得到30936個基因家族,其中單葉省藤和黃藤基因組中分別檢測到44700和44537個同源基因。10個物種中共有的基因家族近6132個(占19.8%),在單葉省藤和黃藤基因組中別檢測到的2366和2707個特有的基因家族(原文圖4b)。此外,家族分析結果顯示其中637個基因家族是棕櫚藤所特有的。這些特異的基因家族大多富集在膜組成成分、轉錄因子活性相關的基因簇上(原文表S13),且聚集在植物與病原體相互作用、植物激素信號轉導相關的KEGG通路中(原文表S14)。

    1.10 系統(tǒng)發(fā)育進化和分化時間分析

    利用獲得的962個種間高度保守的單拷貝直系同源基因,進行棕櫚藤和其他物種的進化關系分析。首先,通過MUSCLE(版本3.8.31)進行蛋白質序列的多序列比對;在此基礎上,基于蛋白質序列進行了DNA編碼序列(CDS)比對。隨后,將所有比對的CDS序列連接起來,使用Perl腳本為每個物種生成一個超級基因。然后,提取了每個密碼子第2個位置上(相位1)的核苷酸使用RAxML(版本8.2.3)構建系統(tǒng)發(fā)育樹,模型為“GTRGAMMA”。結果表明,4種棕櫚科植物分布在一支上,由2個獨立的子分支組成,其中一個分支上包含單葉省藤和黃藤,另一個分支上包含油棕和長葉刺葵(原文圖4a)。

    此外,使用了PAML(版本4.5)的MCMCTree程序并按照參數(shù)“-nsample200000-burnin40000”估計單葉省藤和黃藤與其他8個物種的分化時間。時間標準是根據(jù)參考物種發(fā)生分化的時間推算得出的。結果表明,單葉省藤和黃藤之間的分化時間約在1930萬年前,而棕櫚科的油棕和長葉刺葵的分化時間約在4080萬年前(原文圖4c)。

    1.11 木質素生物合成通路基因家族全基因組鑒定分析

    木質素是一種由木質素單體組成的復雜芳香族聚合物,與纖維素和半纖維素基相互作用共同構成次生細胞壁。木質素通常由3種木質素單體對-羥基苯基(p-hydroxyphenyl,H)、香草醛(vanillin,G)和丁香醛(syringaldehyde,S)組成。本研究利用擬南芥、二穗短柄草、水稻、高梁、毛竹(Phyllostachysedulis)、毛果楊(Populustrichocarpa)、單葉省藤和黃藤等8個物種的基因組進行木質素生物合成通路13個基因家族的全基因鑒定。多數(shù)基因組序列(擬南芥、二穗短柄草、水稻、高梁和毛果楊)從ENSEMBL數(shù)據(jù)庫下載獲得,從毛竹基因組數(shù)據(jù)庫(BambooGenomeDatabase)中下載了毛竹的基因組序列。在前人研究的基礎上,收集得到了經過實驗驗證的140個木質素生物合成通路基因(原文表S15),這些基因作為查詢序列進行進一步的基因功能鑒定。如上所述,在全基因組基因鑒定過程中使用了BLAST比對和保守結構域分析方法。簡單地說,使用已知基因的編碼序列,對所有基因組序列(包括兩種棕櫚藤的基因組序列)進行了目標蛋白的BLAST比對(版本2.2.26),選取E-value<1e-10,相似度>40%且序列的覆蓋率>95%的序列。隨后,利用Pfam-A.hmm數(shù)據(jù)庫(2017年3月31日發(fā)布)對已過濾的序列進行hmmsearch(版本3.1b2)分析,通過人工校正的方法去除結果中結構域不完整的序列。結果表明,單葉省藤和黃藤中木質素相關基因數(shù)量均有明顯增多(原文表3)。單葉省藤和黃藤中每個木質素生物合成通路基因家族平均成員數(shù)量分別約為15和13個,總數(shù)分別為193和172個。在2種棕櫚藤中過氧化物酶基因家族是最常見的基因,單葉省藤的苯丙氨酸解氨酶基因和黃藤中的香豆酸3-羥化酶基因與肉桂酸4-羥化酶基因是最不常見的。棕櫚藤中木質素生物合成基因的數(shù)目擴增可能是由于全基因組復制(WGD)事件造成的,因為全基因組復制可以提供更多的基因拷貝,從而促進具有新功能基因的進化。

    2 結論

    本研究采用多種類型的測序數(shù)據(jù),首次完成了單葉省藤和黃藤染色體水平基因組的組裝,2種棕櫚藤的基因組作為重要的參考基因組,將有助于促進其他藤種的從頭基因組測序組裝和重測序研究,同時通過與不同物種進行比較研究,為物種進化研究提供證據(jù)。借助兩種棕櫚藤高質量基因組數(shù)據(jù),使木質素生物合成通路關鍵基因的鑒定更加便捷,這些候選基因對棕櫚藤的生長發(fā)育非常重要。本研究為進一步對棕櫚藤及相關物種基因組的研究奠定了基礎。

    猜你喜歡
    單葉棕櫚基因組
    亞純函數(shù)關于單葉離散值的正規(guī)定理
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    算子作用下調和函數(shù)類的單葉半徑
    不同因素對單葉蔓荊無性繁殖育苗的影響
    它們可以用來書寫嗎
    棕櫚樹
    棕櫚
    散文詩(2017年17期)2018-01-31 02:34:05
    棕櫚設計有限公司
    風景園林(2017年3期)2017-06-22 13:39:43
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    世界科學(2014年8期)2014-02-28 14:58:31
    我的女老师完整版在线观看| 校园人妻丝袜中文字幕| 国产成人a∨麻豆精品| 99热6这里只有精品| 欧美另类亚洲清纯唯美| 亚洲国产日韩欧美精品在线观看| 你懂的网址亚洲精品在线观看 | 久久久精品欧美日韩精品| av视频在线观看入口| 亚洲欧美中文字幕日韩二区| 亚洲av成人精品一区久久| 白带黄色成豆腐渣| 国内揄拍国产精品人妻在线| 少妇猛男粗大的猛烈进出视频 | 好男人在线观看高清免费视频| 欧美成人午夜免费资源| 亚洲色图av天堂| 日日啪夜夜撸| 久久久久久久午夜电影| 韩国av在线不卡| 国产伦精品一区二区三区四那| 成年女人看的毛片在线观看| 国产成人精品婷婷| 欧美色视频一区免费| 亚洲天堂国产精品一区在线| 精品熟女少妇av免费看| 久久精品影院6| 99久久精品国产国产毛片| av线在线观看网站| 精品国产露脸久久av麻豆 | 亚洲国产欧洲综合997久久,| 欧美色视频一区免费| 亚洲人成网站高清观看| 身体一侧抽搐| 国产高潮美女av| 国产欧美日韩精品一区二区| 男人舔女人下体高潮全视频| 免费人成在线观看视频色| 中文字幕免费在线视频6| 久久久久久久久久黄片| 亚洲欧美成人综合另类久久久 | 欧美xxxx黑人xx丫x性爽| 亚洲精品影视一区二区三区av| 久久韩国三级中文字幕| 国产高清有码在线观看视频| av在线天堂中文字幕| a级毛片免费高清观看在线播放| 国产一区二区三区av在线| 精品人妻偷拍中文字幕| 中文欧美无线码| 网址你懂的国产日韩在线| 纵有疾风起免费观看全集完整版 | 乱系列少妇在线播放| av视频在线观看入口| 青春草亚洲视频在线观看| 日韩三级伦理在线观看| 97超视频在线观看视频| 又粗又硬又长又爽又黄的视频| 国产私拍福利视频在线观看| 国产一区亚洲一区在线观看| 国产av不卡久久| or卡值多少钱| 日韩成人av中文字幕在线观看| 久久久久久大精品| 国产成人aa在线观看| 天堂网av新在线| 国产午夜福利久久久久久| 婷婷六月久久综合丁香| 亚洲av免费高清在线观看| 一级爰片在线观看| 国产爱豆传媒在线观看| 免费观看人在逋| 成年免费大片在线观看| 日韩,欧美,国产一区二区三区 | 国产在视频线在精品| 免费无遮挡裸体视频| 久久精品夜夜夜夜夜久久蜜豆| 一二三四中文在线观看免费高清| 国产91av在线免费观看| 小说图片视频综合网站| 99久久中文字幕三级久久日本| 一级毛片我不卡| 国产亚洲av片在线观看秒播厂 | 91久久精品电影网| 嫩草影院精品99| 最近的中文字幕免费完整| 国产精品野战在线观看| 亚洲精品亚洲一区二区| 国产高清国产精品国产三级 | 六月丁香七月| 精品久久久久久久久久久久久| 亚洲久久久久久中文字幕| 亚洲欧美清纯卡通| 少妇高潮的动态图| 插逼视频在线观看| 欧美又色又爽又黄视频| 国产一级毛片七仙女欲春2| 久久久久久九九精品二区国产| 国产精品av视频在线免费观看| 中文字幕熟女人妻在线| 免费无遮挡裸体视频| 插阴视频在线观看视频| 免费看av在线观看网站| 亚洲综合精品二区| 男插女下体视频免费在线播放| 能在线免费观看的黄片| 亚洲综合精品二区| 三级国产精品片| 亚洲国产日韩欧美精品在线观看| 亚洲人成网站在线观看播放| 日韩一区二区视频免费看| 精品久久久久久久久av| 噜噜噜噜噜久久久久久91| 能在线免费观看的黄片| 亚洲最大成人中文| 亚洲婷婷狠狠爱综合网| 免费黄网站久久成人精品| 亚洲在久久综合| 97热精品久久久久久| 国产精品美女特级片免费视频播放器| 中文资源天堂在线| www日本黄色视频网| 欧美高清性xxxxhd video| 国产欧美日韩精品一区二区| 欧美又色又爽又黄视频| 国产免费视频播放在线视频 | 少妇高潮的动态图| 国产成人免费观看mmmm| 久久人人爽人人片av| 国产伦精品一区二区三区视频9| av在线播放精品| 人人澡人人妻人| 日韩大片免费观看网站| 性高湖久久久久久久久免费观看| 夜夜爽夜夜爽视频| 亚洲av综合色区一区| 亚洲欧美清纯卡通| 日本色播在线视频| 欧美最新免费一区二区三区| 极品少妇高潮喷水抽搐| 搡老乐熟女国产| 日韩精品有码人妻一区| 日韩视频在线欧美| 国产精品国产三级国产专区5o| 黑人欧美特级aaaaaa片| 国产极品粉嫩免费观看在线| 大香蕉久久网| 亚洲美女搞黄在线观看| 午夜老司机福利剧场| 又大又黄又爽视频免费| 精品少妇黑人巨大在线播放| 99热这里只有是精品在线观看| 美国免费a级毛片| 国产黄色视频一区二区在线观看| 国产1区2区3区精品| 观看av在线不卡| 久久久久久久久久人人人人人人| 亚洲国产精品999| 一区二区三区四区激情视频| 爱豆传媒免费全集在线观看| 国产精品成人在线| 少妇的逼水好多| 亚洲国产精品成人久久小说| 午夜福利网站1000一区二区三区| av线在线观看网站| 黄色毛片三级朝国网站| 男人操女人黄网站| 亚洲精华国产精华液的使用体验| 精品一品国产午夜福利视频| 在线观看免费视频网站a站| videos熟女内射| 亚洲av在线观看美女高潮| 另类精品久久| 亚洲高清免费不卡视频| 国产精品三级大全| 国产视频首页在线观看| 国产黄频视频在线观看| 亚洲精品自拍成人| 久久久久久久久久久久大奶| 丝瓜视频免费看黄片| 精品第一国产精品| 超色免费av| 人成视频在线观看免费观看| 伦精品一区二区三区| 亚洲美女搞黄在线观看| 久久久久精品性色| 久久97久久精品| 18禁裸乳无遮挡动漫免费视频| h视频一区二区三区| 免费观看性生交大片5| 国产国拍精品亚洲av在线观看| 亚洲欧美中文字幕日韩二区| 中文字幕av电影在线播放| 亚洲av.av天堂| 男女下面插进去视频免费观看 | 亚洲欧洲精品一区二区精品久久久 | 日韩欧美精品免费久久| 亚洲欧洲精品一区二区精品久久久 | 色网站视频免费| 免费观看a级毛片全部| 考比视频在线观看| a级毛片黄视频| 9色porny在线观看| 夫妻午夜视频| 久久久久视频综合| 国产亚洲欧美精品永久| 国产极品粉嫩免费观看在线| 国产成人精品无人区| 国产精品国产av在线观看| 亚洲精品日本国产第一区| 十分钟在线观看高清视频www| 日本wwww免费看| 黄片无遮挡物在线观看| 汤姆久久久久久久影院中文字幕| 久久这里有精品视频免费| av黄色大香蕉| 丝袜在线中文字幕| 天堂俺去俺来也www色官网| 女性被躁到高潮视频| 久久 成人 亚洲| 久久97久久精品| xxxhd国产人妻xxx| 狂野欧美激情性xxxx在线观看| 日本av免费视频播放| 全区人妻精品视频| 成人毛片a级毛片在线播放| av.在线天堂| 亚洲av电影在线进入| 另类亚洲欧美激情| h视频一区二区三区| 国产毛片在线视频| 成人免费观看视频高清| 免费观看av网站的网址| 秋霞伦理黄片| 国产乱人偷精品视频| 妹子高潮喷水视频| 一二三四在线观看免费中文在 | 国产精品久久久久久精品电影小说| 丝袜脚勾引网站| 亚洲国产成人一精品久久久| 在线观看一区二区三区激情| 日韩av不卡免费在线播放| 高清不卡的av网站| tube8黄色片| 中文欧美无线码| 我要看黄色一级片免费的| 久久久欧美国产精品| 久久久久精品久久久久真实原创| 国产 精品1| 在线天堂中文资源库| 啦啦啦在线观看免费高清www| 亚洲精品456在线播放app| 90打野战视频偷拍视频| 搡老乐熟女国产| 亚洲精品一区蜜桃| 一区在线观看完整版| 男人爽女人下面视频在线观看| 午夜日本视频在线| 成人无遮挡网站| 久久久精品免费免费高清| 成人国产麻豆网| 亚洲国产色片| 成人免费观看视频高清| 两个人看的免费小视频| 日韩视频在线欧美| www.色视频.com| 欧美xxxx性猛交bbbb| 久久久国产精品麻豆| 一级毛片电影观看| 丁香六月天网| 国产一区二区三区av在线| 欧美人与性动交α欧美软件 | 免费大片18禁| 赤兔流量卡办理| 国产欧美日韩综合在线一区二区| 丝袜美足系列| 人妻一区二区av| 我的女老师完整版在线观看| 欧美成人午夜精品| 男女下面插进去视频免费观看 | 国产av码专区亚洲av| 国产综合精华液| 午夜影院在线不卡| 人人妻人人澡人人看| 免费黄色在线免费观看| 亚洲国产精品一区二区三区在线| 综合色丁香网| 久久精品aⅴ一区二区三区四区 | 少妇的逼好多水| 亚洲精品国产av蜜桃| 国产永久视频网站| 最后的刺客免费高清国语| 亚洲成av片中文字幕在线观看 | 男女国产视频网站| 涩涩av久久男人的天堂| 精品国产乱码久久久久久小说| 国产精品99久久99久久久不卡 | 国产精品人妻久久久久久| 国产麻豆69| a级毛片黄视频| 国产成人精品福利久久| 人人澡人人妻人| 亚洲三级黄色毛片| av片东京热男人的天堂| 晚上一个人看的免费电影| 99视频精品全部免费 在线| 9色porny在线观看| 日韩一区二区视频免费看| 久久精品国产自在天天线| 国产 精品1| 精品酒店卫生间| 亚洲综合色惰| 亚洲婷婷狠狠爱综合网| 热re99久久国产66热| 午夜91福利影院| 熟女人妻精品中文字幕| a级毛色黄片| 精品视频人人做人人爽| 69精品国产乱码久久久| 99热6这里只有精品| 亚洲人成77777在线视频| 国产精品女同一区二区软件| 亚洲丝袜综合中文字幕| 五月玫瑰六月丁香| 十八禁高潮呻吟视频| 欧美日韩视频精品一区| 毛片一级片免费看久久久久| 国产精品久久久久久精品古装| av黄色大香蕉| 在线天堂中文资源库| 国产精品 国内视频| 亚洲av男天堂| 午夜福利影视在线免费观看| 精品酒店卫生间| 成人18禁高潮啪啪吃奶动态图| 午夜影院在线不卡| 日韩三级伦理在线观看| 老熟女久久久| 日韩精品有码人妻一区| 晚上一个人看的免费电影| 欧美日韩成人在线一区二区| 最近的中文字幕免费完整| 亚洲av日韩在线播放| 久久韩国三级中文字幕| 在线精品无人区一区二区三| 国国产精品蜜臀av免费| tube8黄色片| 国产国语露脸激情在线看| 欧美xxxx性猛交bbbb| 亚洲欧美精品自产自拍| 99九九在线精品视频| 亚洲在久久综合| 一级毛片电影观看| 成年动漫av网址| 黑人高潮一二区| 国产在线一区二区三区精| 国产一区二区三区av在线| 在线免费观看不下载黄p国产| 久久99一区二区三区| 99九九在线精品视频| 免费看不卡的av| 五月伊人婷婷丁香| 精品第一国产精品| 色视频在线一区二区三区| 亚洲av免费高清在线观看| 国产免费又黄又爽又色| 日韩成人伦理影院| 国产精品国产三级专区第一集| 午夜福利影视在线免费观看| 精品国产一区二区三区久久久樱花| 黄色毛片三级朝国网站| 免费大片黄手机在线观看| 两个人看的免费小视频| 国产伦理片在线播放av一区| 狂野欧美激情性xxxx在线观看| 久热这里只有精品99| 大片电影免费在线观看免费| 国产xxxxx性猛交| 99香蕉大伊视频| 欧美 亚洲 国产 日韩一| 日日啪夜夜爽| 久久精品aⅴ一区二区三区四区 | 欧美xxxx性猛交bbbb| 亚洲性久久影院| 一区二区av电影网| 亚洲欧美清纯卡通| 少妇的逼水好多| 国产熟女午夜一区二区三区| 久久久久国产精品人妻一区二区| 亚洲一码二码三码区别大吗| 在现免费观看毛片| 成人毛片60女人毛片免费| 春色校园在线视频观看| 欧美最新免费一区二区三区| 99九九在线精品视频| 久久精品人人爽人人爽视色| av卡一久久| www.av在线官网国产| 久久久久久久精品精品| 欧美日本中文国产一区发布| 亚洲 欧美一区二区三区| 亚洲国产精品999| 国产日韩一区二区三区精品不卡| 不卡视频在线观看欧美| 美女福利国产在线| 国产av精品麻豆| 99九九在线精品视频| 国产成人精品在线电影| 午夜91福利影院| 99久久综合免费| 韩国精品一区二区三区 | 亚洲伊人久久精品综合| 婷婷色综合www| 国产日韩欧美亚洲二区| 国产在线免费精品| 中文字幕人妻熟女乱码| 人妻人人澡人人爽人人| h视频一区二区三区| 草草在线视频免费看| 99热全是精品| 插逼视频在线观看| 高清不卡的av网站| 女人久久www免费人成看片| 久久韩国三级中文字幕| 精品熟女少妇av免费看| 亚洲国产精品成人久久小说| 99国产综合亚洲精品| 夜夜骑夜夜射夜夜干| 免费观看在线日韩| 黄色视频在线播放观看不卡| 免费高清在线观看视频在线观看| 妹子高潮喷水视频| 视频区图区小说| 九九爱精品视频在线观看| 亚洲国产av新网站| 亚洲国产毛片av蜜桃av| 午夜久久久在线观看| 精品一品国产午夜福利视频| 青春草亚洲视频在线观看| 高清视频免费观看一区二区| 国产av国产精品国产| 成人亚洲精品一区在线观看| 精品久久久精品久久久| 午夜福利网站1000一区二区三区| 亚洲精品成人av观看孕妇| 街头女战士在线观看网站| 蜜桃国产av成人99| 久久这里只有精品19| 国产亚洲av片在线观看秒播厂| 国产欧美另类精品又又久久亚洲欧美| 久热这里只有精品99| 青春草亚洲视频在线观看| 99re6热这里在线精品视频| 亚洲精品av麻豆狂野| 亚洲人成77777在线视频| 十八禁高潮呻吟视频| 日韩一本色道免费dvd| 少妇人妻精品综合一区二区| 日产精品乱码卡一卡2卡三| 女的被弄到高潮叫床怎么办| 极品人妻少妇av视频| 亚洲欧美色中文字幕在线| 亚洲综合色网址| 97超碰精品成人国产| 亚洲精品av麻豆狂野| 在线天堂最新版资源| 在现免费观看毛片| 国产免费福利视频在线观看| 99久久精品国产国产毛片| 久久 成人 亚洲| 美女中出高潮动态图| 国产一区二区在线观看日韩| a级毛色黄片| 制服丝袜香蕉在线| 在线观看免费高清a一片| 成人综合一区亚洲| 国产精品三级大全| 精品一区二区三卡| 十分钟在线观看高清视频www| 日本91视频免费播放| 亚洲国产精品专区欧美| 亚洲图色成人| 久久久久久伊人网av| 久久久国产欧美日韩av| 桃花免费在线播放| 99久久人妻综合| 一区二区av电影网| 欧美成人精品欧美一级黄| 在线观看免费日韩欧美大片| 一级黄片播放器| 2022亚洲国产成人精品| 免费黄色在线免费观看| 亚洲av中文av极速乱| 狂野欧美激情性bbbbbb| 久久人妻熟女aⅴ| 免费日韩欧美在线观看| 99视频精品全部免费 在线| av电影中文网址| 欧美日韩综合久久久久久| 亚洲综合精品二区| 亚洲经典国产精华液单| 午夜福利影视在线免费观看| 国产男女超爽视频在线观看| 99久久综合免费| 99国产综合亚洲精品| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 精品国产国语对白av| 色网站视频免费| av不卡在线播放| 日本黄大片高清| 丝袜脚勾引网站| 亚洲精品aⅴ在线观看| 欧美 亚洲 国产 日韩一| 亚洲av日韩在线播放| 九色亚洲精品在线播放| 久久久亚洲精品成人影院| 国产精品嫩草影院av在线观看| 九九在线视频观看精品| 性色avwww在线观看| 欧美xxⅹ黑人| 欧美国产精品va在线观看不卡| 黑丝袜美女国产一区| 美女内射精品一级片tv| 国产av精品麻豆| 最近手机中文字幕大全| 最近中文字幕2019免费版| 大片免费播放器 马上看| 免费av不卡在线播放| 国产免费视频播放在线视频| 久久国产精品男人的天堂亚洲 | 国产成人精品婷婷| 我要看黄色一级片免费的| 美女脱内裤让男人舔精品视频| 好男人视频免费观看在线| 看免费成人av毛片| 黑人高潮一二区| 国产成人精品一,二区| 国产成人精品无人区| 美女国产高潮福利片在线看| 久久精品国产亚洲av涩爱| 日日摸夜夜添夜夜爱| 精品第一国产精品| 纯流量卡能插随身wifi吗| 高清黄色对白视频在线免费看| 国产片特级美女逼逼视频| 国产毛片在线视频| 黄片播放在线免费| 国产色婷婷99| 国产成人精品婷婷| 国产69精品久久久久777片| 免费高清在线观看视频在线观看| 中文乱码字字幕精品一区二区三区| 久久久亚洲精品成人影院| 免费日韩欧美在线观看| 天天躁夜夜躁狠狠久久av| 国产成人免费无遮挡视频| 国产成人精品在线电影| 观看美女的网站| 五月开心婷婷网| 美女主播在线视频| 国产亚洲一区二区精品| 国产av精品麻豆| 日日爽夜夜爽网站| 久久久久久久久久人人人人人人| 91成人精品电影| 18禁动态无遮挡网站| 欧美日韩视频精品一区| 精品99又大又爽又粗少妇毛片| 少妇被粗大猛烈的视频| 天天躁夜夜躁狠狠躁躁| 日日撸夜夜添| 男女无遮挡免费网站观看| 在线观看国产h片| 在线天堂中文资源库| 18禁在线无遮挡免费观看视频| 久久人妻熟女aⅴ| 侵犯人妻中文字幕一二三四区| 街头女战士在线观看网站| 中文字幕免费在线视频6| 国产麻豆69| 国产日韩欧美在线精品| 欧美精品av麻豆av| 欧美+日韩+精品| 国产精品久久久久久av不卡| 中文字幕另类日韩欧美亚洲嫩草| 五月天丁香电影| 亚洲av国产av综合av卡| 香蕉精品网在线| www.色视频.com| h视频一区二区三区| 日韩成人av中文字幕在线观看| 夫妻性生交免费视频一级片| 久久久久网色| 男人操女人黄网站| 日本91视频免费播放| 夫妻午夜视频| 国产一区二区激情短视频 | 国产一区二区在线观看日韩| 国产日韩欧美在线精品| 久久ye,这里只有精品| 黄色配什么色好看| 春色校园在线视频观看| 最新的欧美精品一区二区| 欧美成人精品欧美一级黄| 我的女老师完整版在线观看| 80岁老熟妇乱子伦牲交| av国产精品久久久久影院| 99国产综合亚洲精品| 一二三四在线观看免费中文在 | 亚洲,一卡二卡三卡| 高清黄色对白视频在线免费看| 在线观看www视频免费| 亚洲伊人久久精品综合| 精品人妻在线不人妻|