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

    基于5993個核基因的被子植物系統(tǒng)發(fā)育關系研究

    2020-03-02 07:46:59金鑫程書楊拓余慷段肖霞倪雪梅李世明張耕耘
    廣西植物 2020年1期

    金鑫 程書 楊拓 余慷 段肖霞 倪雪梅 李世明 張耕耘

    摘 要:系統(tǒng)發(fā)育關系的構建對被子植物分類及進化研究非常重要。長期以來,被子植物系統(tǒng)發(fā)育的研究,大多使用質體基因、線粒體基因或少數保守的單拷貝核基因。該研究從已注釋基因組或轉錄組中搜集88種被子植物(包含58目)的核基因集;通過對其進行同源基因聚類及去旁系同源基因,獲得了5 993個一對一的直系同源基因家族(即對于每個基因家族,每種植物最多一條序列,最少包含50個物種);使用截取各種不同數目基因集的DNA或氨基酸序列,采用串聯(lián)法(concatenation)和溯祖法(coalescence),共構建了20棵進化樹。比較這些進化樹,雖然大部分結果支持APG IV中描述的被子植物主要支系之間的關系[(真雙子葉植物,單子葉植物),木蘭類植物],但真雙子葉植物內部各目分支的演化關系與APG IV有一個很大的不同,即認為檀香目和石竹目是薔薇類植物的姊妹群?;谶@些進化樹,估算了被子植物各目分支的分化時間,結果表明被子植物的起源時間為237.78百萬年前(95%置信區(qū)間為202.6~278.08),與主流觀點認為的225百萬年至240百萬年前一致。以上結果為構建進化樹提供了一種可行性策略,這種方法允許使用基因數目更多而計算速度更快。

    關鍵詞:系統(tǒng)發(fā)育關系,被子植物,核基因,同源基因聚類,串聯(lián)法,溯祖法,分化時間

    中圖分類號:Q949.4文獻標識碼:A

    文章編號:1000-3142(2020)01-0044-16

    Abstract:Construction of phylogeny is important for classification and research of angiosperms. For a long time,angiosperm phylogeny has been analysed using plastid genes,mitochondrial genes or a few conserved single-copy nuclear genes. We collected nuclear gene sets of 88 species of angiosperm (contains 58 orders) from annotated genomes or transcriptomes. By using a combined homology- and phylogeny tree-based approach,we obtained a total of 5 993 one-to-one ortholog groups (one sequence of each species for each ortholog group),each of which was represented by at least 50 species. Then,a total of 20 species trees were reconstructed using? methods with different combinations of reconstruction (concatenation-based and coalescence-based) and sequence type (nucleotide or amino acid) for gene data sets with different gene occupancy values. Most of the resulting topologies support the relationships of the major clades of angiosperm as described in APG IV,but present different deep relationships among major clades in eudicots phylogeny such as the placement of Santalales and Caryophyllales as sisters to Rosids. We estimate the divergence times of the major clades of angiosperm and conclude that the origin of angiosperm is about 237.78 million years ago(95% confidence interval is 202.6-278.08),which is in accordance with the previously accepted 225 million years to 240 million years ago. This study provides an efficient strategy for building phylogenetic trees using thousands of genes with ultrafast calculation.

    Key words:phylogeny,angiosperms,nuclear genes,ortholog inference,concatenation,coalescence,divergence time

    系統(tǒng)發(fā)育樹的正確構建對植物分類及進化研究非常重要。進化樹構建的準確度主要受以下因素的影響。其一,所使用的數據集的種類及大小。不僅使用形態(tài)性狀數據、質體基因、線粒體基因及核基因序列建立的進化樹不一樣(Endress & Doyle,2009; Soltis et al.,2011; Ruhfel et al.,2014; Zeng et al.,2014),使用全長核酸序列或僅使用基因密碼子某個位點的核酸序列及氨基酸序列所構建的進化樹也有所不同(Wickett et al.,2014);其二,構建樹的方法及模型。方法有串聯(lián)法(concatenation)和溯祖法(coalescence)。串聯(lián)法是將所有基因串聯(lián)作為一個整體,使用軟件RAxML(Stamatakis,2014)或iqtree(Nguyen et al.,2015)構建系統(tǒng)發(fā)育樹;溯祖法是先對每個基因建樹,再使用軟件ASTRAL(Zhang et al.,2017)建立所有基因樹的共有樹(Wickett et al.,2014)。而構建系統(tǒng)發(fā)育樹使用的模型更是多種多樣,如核酸模型GTR、HKY、JC、F81、K2P、K3P、K81uf等,蛋白質模型LG、Poisson、cpREV、mtREV、Dayhoff、mtMAM、JTT、WAG等(Nguyen et al.,2015)。

    被子植物是植物界最高等且種類最多的一類,它們在地球上占據著絕對優(yōu)勢?,F在已報道被子植物有352 000種(http://www.theplantlist.org/),屬于416科和64目,各目之間的演化關系一直是研究的熱點和爭論的焦點。被子植物除了最基部的三個目:無油樟目(Amborellales)、睡蓮目(Amborellales)和木蘭藤目(Austrobaileyales),又稱ANITA組,其余的(99.95%)可以分為五類:木蘭類植物(magnoliids)、單子葉植物(monocots)、真雙子葉植物(eudicots)、金栗蘭科(Chloranthaceae)和金魚藻科(Ceratophyllaceae)。這五類的系統(tǒng)演化拓撲關系一直存在爭論,Zeng et al.(2014)總結了已經發(fā)表的五種主要的拓撲關系(圖1:A-E),其中A是最主流的,也是APG IV(The Angiosperm Phylogeny Group,2016)的拓撲結構。Soltie et al.(2011)使用17個基因串聯(lián)(包括質體基因、線粒體基因和核基因)為640種植物構建的系統(tǒng)發(fā)育進化樹和Ruhfel et al.(2014)使用78個質體基因串聯(lián)為360種植物構建的進化樹,支持主流A拓撲結構。Wickett et al.(2014)使用674個核基因串聯(lián)為92種植物構建的進化樹和Zeng et al.(2014)使用59個核基因串聯(lián)為61種植物構建的進化樹,支持B拓撲結構。除此之外,Qiu et al.(2010)使用4個線粒體基因為380種植物構建的進化樹,支持C拓撲結構;Endress & Doyle(2009)使用形態(tài)性狀構建的進化樹,支持D拓撲結構;Zhang et al.(2012)使用5個核基因為91種植物構建的進化樹,支持E拓撲結構。

    去掉金栗蘭科和金魚藻科后,單子葉植物、木蘭類植物、真雙子葉植物之間的系統(tǒng)發(fā)育關系有三種:(真雙子葉植物,單子葉植物),木蘭類植物;(真雙子葉植物,木蘭類植物),單子葉植物;(單子葉植物,木蘭類植物),真雙子葉植物。Lu et al.(2018)使用4個質體基因和1個線粒體基因分析了5 864種中國被子植物(幾乎包括所有中國地區(qū)被子植物)的系統(tǒng)發(fā)育關系,其構建的進化樹支持拓撲結構[(真雙子葉植物,單子葉植物),木蘭類植物]。Chen et al.(2019)發(fā)布了木蘭類植物鵝掌楸(Liriodendron)基因組,使用其502個核基因及溯祖法為18種植物構建的進化樹,同樣支持拓撲結構[(真雙子葉植物,單子葉植物),木蘭類植物]。Chaw et al.(2019)發(fā)布了另一個木蘭類植物牛樟(stout camphor tree)基因組,使用其211個核基因為13種植物構建的進化樹,支持拓撲結構[(真雙子葉植物,木蘭類植物),單子葉植物]。Li et al.(2019)使用2 881種被子植物的質體基因組的80個基因重建了被子植物高分辨率的系統(tǒng)發(fā)育樹,支持拓撲結構[(真雙子葉植物,單子葉植物),木蘭類植物]。從上述已有的研究中,我們發(fā)現,使用核基因串聯(lián)法建立的進化樹基本都支持拓撲結構[(真雙子葉植物,木蘭類植物),單子葉植物],使用核基因溯祖法、質體和線粒體基因建立的進化樹基本都支持拓撲結構[(真雙子葉植物,單子葉植物),木蘭類植物]。

    真雙子葉植物內部各目的系統(tǒng)發(fā)育關系也存在爭論(圖1:F-K),真雙子葉植物除了最基部的毛茛目(Ranunculales)、山龍眼目(Proteales)、昆欄樹目(Trochodendrales)、黃楊目(Buxales)和洋二仙草目(Gunnerales),其余的可以分為兩類:薔薇類植物(Rosids)和菊類植物(Asterids)。這兩類植物的基部有6個目的系統(tǒng)發(fā)育關系比較混亂,即五椏果目(Dilleniales)、虎耳草目(Saxifragales)、葡萄目(Vitales)、檀香目(Santalales)、智利藤目(Berberidopsidales)及石竹目(Caryophyllales)。Zeng et al.(2017)總結了已經發(fā)表的六種主要的拓撲關系(圖1:F-K),其中K是APG IV中所認可的拓撲結構。Moore et al.(2010)使用83個質體基因為86種植物構建的進化樹,支持“五椏果目是薔薇類植物的姊妹群”;Soltie et al.(2011)使用17個基因串聯(lián)(包括質體基因、線粒體基因和核基因)為640種植物構建的進化樹和Moore et al.(2011)使用質體IR序列為87種植物構建的進化樹,支持“五椏果目是菊類植物的姊妹群”;Worberg et al.(2007)等使用五個基因組區(qū)域序列為56種植物構建的進化樹和Moore et al.(2011)使用質體IR序列為244種植物構建的進化樹,及APG IV都支持“五椏果目是薔薇類植物和菊類植物共同的姊妹群”。大部分研究都支持“葡萄目和虎耳草目是薔薇類植物的姊妹群,智利藤目、檀香目和石竹目是菊類植物的姊妹群”(Worberg et al.,2007; Moore et al.,2010,2011; Yang et al.,2015);Zeng et al.(2017)使用504個核基因和溯祖法為100種植物構建的進化樹,支持“檀香目和智利藤目是薔薇類植物的姊妹群”。

    被子植物的起源及進化一直是植物學界研究和爭論的熱點。在古生物學界,很長時期內,被子植物的最早化石記錄都是白堊紀125百萬年前,也是最早的真雙子葉植物化石記錄(Herendeen,1995)。Fu et al.(2018)發(fā)現了早侏羅紀地層(約175百萬年前)中的“南京花”,其具有花萼、花瓣、雌蕊,有明顯的杯托、下位子房上位花、樹狀的花柱,其種子/胚珠確實是被完全包裹著,子房壁將種子與外界完全隔絕,這都滿足了被子植物判斷標準。“南京花”的發(fā)現,將被子植物最早化石記錄向前推進了約5 000萬年,并填補了被子植物化石記錄(125百萬年前)與分子鐘推算時間(225百萬年至240百萬年前)之間的“侏羅紀空缺”(Jurassic gap)(Li et al.,2019)。目前,大多數基于系統(tǒng)進化樹的被子植物分化時間估計研究,都認為被子植物的起源為三疊紀225百萬年至240百萬年前(Magallon,2010; Smith et al.,2010; Zeng et al.,2014; Mandel,2019),這與起傳粉作用的核心植食性鱗翅目昆蟲的起源時間(約230百萬年前)一致(Zeng et al.,2014; Li et al.,2019)。

    本研究使用超過5 000個核基因的核酸及蛋白序列,用兩種進化樹構建方法分析了88種被子植物的系統(tǒng)發(fā)育關系(包括87科58目),并對各進化分支的分化時間進行了估計(總流程如圖2)。為了得到準確可靠的被子植物系統(tǒng)發(fā)育進化樹,我們對5 000多個核基因進行了拆分,得到了包含不同基因數量的多個數據集,并對各個數據集進行系統(tǒng)發(fā)育樹的構建,最后比較了所得到的20棵系統(tǒng)發(fā)育進化樹之間的一致性。

    1 材料與方法

    1.1 材料

    我們收集了1個裸子植物(Ginkgo biloba作為外類群)基因組、43個被子植物基因組(主要來自NCBI和PHYTOZOME數據庫)、43個被子植物已拼接轉錄組(http://www.onekp.com/public_data.html )及2個被子植物RNA-seq數據(其中無葉蓮Petrosavia sakurai是本研究測序的物種),A-E. 五類被子植物間[金栗蘭科(Chl)、金魚藻科(Cer)、木蘭類植物(Mag)、單子葉植物(Mon)及真雙子葉植物(Eud)]5種代表性拓撲結構; F-K. 真雙子葉植物內部各目間6種代表性拓撲結構。

    A-E. Five representative topologies among eudicots (Eud),monocots (Mon),magnoliids (Mag),Ceratophyllaceae (Cer) and Chloranthaceae (Chl); F-K. Six representative topologies among eudicots.其中被子植物共包含87科58目(表1)。

    1.2 基于基因組序列的直系同源基因鑒定

    我們使用Yang & Smith(2014)報道的方法,對43個植物基因組的基因集進行同源基因聚類分析。首先,使用軟件BLASTN v2.6.0+ 對43個基因集CDS序列進行all-by-all blast,每條序列取最佳的1 000條比對結果,去掉比對長度小于1/3總長的序列,修剪未比對上的末端序列。然后,使用MCL軟件(Van,2000)進行同源基因聚類(inflation value = 1.4),去除少于20個植物的基因家族,剩余基因家族使用MAFFT v7.310軟件(Katoh & Standley,2013)進行多序列比對(maximum iterative refinement cycles = 1 000),使用PHYUTILITY v2.2.6軟件(Smith & Dunn,2008)修剪缺失率大于90%的位點,使用軟件RAXML v8.2.11(Stamatakis,2014)對修剪后的多序列比對數據估算系統(tǒng)進化樹(model = GTRCAT)。最后修剪掉進化樹上的所有旁系同源基因枝,修剪枝長大于0.6的枝、比姐妹枝長十倍的末端枝,單源且全部同樣品的枝只保留一個,修剪枝長比預期堿基替換率大0.3倍的內部枝,再使用MO方法(Yang & Smith,2014)去除所有剩余的旁系同源枝,獲得one-to-one同源基因家族(即每個樣品最多一條序列),只保留大于20個樣品的基因家族。

    1.3 轉錄組及外類群數據處理

    我們對兩個來自兩個科(無葉蓮科Petrosavia sakurai和鴉跖草科Cyanotis arachnoidea)的RNA-seq數據從頭拼接。首先使用Trimmomatic v0.38軟件(Bolger et al.,2014)過濾原始reads數據(參數HEADCROP:15 LEADING:20 TRAILING:20 SLIDINGWINDOW:5:20 MINLEN:50 AVGQUAL:20),再使用Trinity v2.6.6軟件(Grabherr et al.,2011)拼接(min contig length=150 bp),最后使用TransDecoder v5.5.0(https://github.com/TransDecoder/TransDecoder/releases/tag/TransDecoder-v5.5.0 )進行CDS和蛋白質序列預測(參考數據庫為Swissprot和Pfam-A)。將得到的這兩個物種的基因集、從onekp數據庫下載得到的43種被子植物的基因集和1個裸子植物(Ginkgo biloba)的基因集,使用HaMStR v13.2.6軟件(Ebersberger et al.,2009)合并到利用基因組數據得到的同源基因家族中,最終只保留大于50個樣品的基因家族。

    1.4 系統(tǒng)發(fā)育進化樹構建

    我們采用兩種方法串聯(lián)法 (concatenation)和溯祖法(coalescence),并分別使用CDS序列和氨基酸序列構建進化樹。無論是CDS序列還是蛋白質序列,都使用PRANK v.170427軟件(http://wasabiapp.org/software/prank/)進行多序列比對,使用PHYUTILITY v2.2.6軟件(Smith & Dunn,2008)修剪缺失率大于70%的位點,其中CDS序列需去除長度小于300個堿基的序列,蛋白質序列需去除長度小于100個氨基酸的序列。

    溯祖法,先對每個基因使用RAxML v8.2.11軟件(默認參數)(Stamatakis,2014)畫樹,再使用ASTRAL v5.5.9軟件(Zhang et al.,2017)處理所有基因樹,得到共有樹,參數設置“-t 1--gene-only”以獲得bootstrap值和基因支持率,枝長使用iqtree v1.5.5軟件(Nguyen et al.,2015)獲得。

    串聯(lián)法,先使用PartitionFinder v2.1.1軟件(Lanfear et al.,2009)對串聯(lián)序列進行分區(qū)和進化模型檢測,從而設置較合理的分區(qū)和為每個分區(qū)選擇合理的進化模型。對CDS序列檢測下列的四個分區(qū)策略(表2):no partitioning,partitioning by each codon position (three partitions),partitioning by gene和partitioning by each codon position within each gene。對蛋白質序列檢測下列兩個分區(qū)策略:no partitioning和 partitioning by gene。參數設置如下:branch lengths = linked;model_selection = aicc;search = user;models = GTR,GTR+G,GTR+I+G(CDS序列)或者models = LG+G,LG+I+G,WAG+G,WAG+I+G(蛋白質序列)。再使用iqtree v1.5.5軟件畫樹[1000 ultrafast bootstrap replicates(Von Haeseler et al.,2013),-spp設置最優(yōu)分區(qū)策略],基因支持率使用ASTRAL v5.5.9軟件(-t 1)獲得。最后使用軟件Evolview v2(He et al.,2016)對獲得的所有進化樹進行美化。

    A. 每個同源基因家族含有的基因數目; B. 每個樣品含有的同源基因家族數目。

    A. Gene number of each orthologous gene family; B. Number of orthologous genes for each sample.

    1.5 分化時間估計

    我們使用PAML v4.9軟件包 (Yang,2007)的MCMCTREE程序進行分化時間估計,輸入拓撲結構為綜合20棵進化樹的最佳拓撲結構(即使用742個基因的CDS序列串聯(lián)法獲得的拓撲結構),輸入序列為742個基因的CDS序列。我們先對每個基因都分別估計分化時間,再綜合742個基因的分析結果(即每個節(jié)點取所有基因的平均值)獲得最終的分化時間樹。拓撲結構的枝長使用JONES+gamma堿基替換模型獲得;rgene gamma設定為G(1,4.5);sigma2 gamma設定為G(1,4.5);clock設定為3;Markov chain Monte Carlo(MCMC)設定為burnin=50 000,sampfreq=100,nsample=10 000。對每個基因,都是分別運行兩次獨立的MCMC(即不同的random seeds),使用Tracer v1.7軟件(https://github.com/beast-dev/tracer/releases/tag/v1.7.1 )觀察運行結果是否穩(wěn)定和收斂,所有節(jié)點及參數的effective sample size是否大于200。九個化石校準設定如下:銀杏分化時間為290百萬年至310百萬年前(Gao et al.,1989);單子葉植物和真雙子葉植物分化時間為130百萬年至200百萬年前(Kumar et al.,2017);真雙子葉植物共同祖先(即最早的雙子葉植物化石記錄)為125 百萬年前(Herendeen,1995; Zeng et al.,2014);山龍眼目(Proteales)的共同祖先為108.8百萬年前(Crane et al.,1996);葡萄目(Vitales)與其余薔薇類植物間分化時間為105百萬年至115百萬年前(Fawcett et al.,2009; Kumar et al.,2017);Arabidopsis thaliana與Populus trichocarpa間分化時間為97百萬年至109百萬年前(Kumar et al.,2017);豆目(Fabales)與殼斗目(Fagales)間分化時間為93.5百萬年前(Friis et al.,1996);山茱萸目(Cornales)共同祖先為85.8百萬年前(Takahashi et al.,2002);唇形目(Lamiales)共同祖先為44.3百萬年前(Call et al.,1992)。

    2 結果與分析

    2.1 直系同源基因鑒定

    我們對44個植物基因組基因集和45個已拼接轉錄組CDS序列進行同源基因聚類,并使用Yang & Smith(2014)報道的方法,去除所有旁系同源基因,最終獲得大于50個樣品的one-to-one 基因家族(即每個樣品最多一條序列)共5 993個(圖3:A),各種植物的基因覆蓋率從33.57%到97.85%,平均為80.40%(圖3:B)。

    2.2 系統(tǒng)發(fā)育進化樹構建

    我們采用串聯(lián)和溯祖法共構建了20棵進化樹,并比較它們之間的不同(圖4),以評估樹的穩(wěn)定性。CDS序列和蛋白質序列,都分別使用五個數據集,總共構建20棵樹(5棵CDS串聯(lián)法樹,5棵CDS溯祖法樹,5棵AA串聯(lián)法樹和5棵AA溯祖法樹)。這5個數據集分別包含5 928個orthologs(≥50 samples)、3 384個orthologs(≥70 samples)、1 791個orthologs(≥80 samples)、742個orthologs(≥85 samples)及42個orthologs(≥89 samples)。

    這20棵進化樹主要是為了進一步確定圖1中五類被子植物間演化關系和真雙子葉植物內部各目間系統(tǒng)發(fā)育關系。這些進化樹中的大多數,是與使用742個基因CDS序列(共4 069 848位點)串聯(lián)方法建立的進化樹高度一致的(圖5)(使用3 384個基因AA序列建立的進化樹,和使用1 791個基因AA序列建立的進化樹,也是相同的最佳拓撲結構)。

    2.2.1 木蘭類植物、單子葉植物及雙子葉植物間演化關系 無論核酸序列還是蛋白質序列,使用串聯(lián)法和溯祖法建立的進化樹基本都支持拓撲結構[(真雙子葉植物,單子葉植物),木蘭類植物](圖4)。

    2.2.2 金栗蘭科與金魚藻科 我們的研究表示,金魚藻科是真雙子葉植物的姊妹群,這與前人的研究結果一致(圖4)。但金栗蘭科是所有被子植物(除ANITA外)的基底旁系群,這與APG IV認為的“金栗蘭科是木蘭類植物的姊妹群”是不同的。

    2.2.3 雙子葉植物內部各目的系統(tǒng)發(fā)育關系 我們的研究認為,五椏果科是薔薇類植物和菊類植物共同的姊妹群,虎耳草目是薔薇類植物的姊妹群,這都與APG IV一致(圖4)。

    APG IV認為“檀香目和石竹目是菊類植物的姊妹群”,而我們的研究否定了這一結論:20棵進化樹中,所有結果都支持“石竹目是薔薇類植物的姊妹群”;大部分支持“檀香目是薔薇類植物的姊妹群”,這與Zeng et al.( 2017)的研究結果一致;少部分支持“檀香目是薔薇類植物和菊類植物共同的姊妹群”(圖4)。

    APG IV認為“智利藤目是菊類植物的姊妹群”,而我們的研究只有少部分支持這一結論。使用蛋白質序列建立的進化樹,無論串聯(lián)還是溯祖法,都支持“智利藤目是薔薇類植物和菊類植物共同的姊妹群”。使用核酸序列建立的進化樹,隨著基因數目的增多,逐漸轉變?yōu)橹С帧爸抢倌渴蔷疹愔参锏逆⒚萌骸?,與APG IV一致(圖4)。

    2.3 分化時間估計

    基于742個基因CDS序列串聯(lián)方法建立的進化樹,我們估計了被子植物的分化時間(圖6)。我們認為被子植物的起源時間為237.78百萬年前(95%置信區(qū)間為202.6~278.08),與主流觀點認為的225百萬年至240百萬年前一致 (Magallon,2010; Smith et al.,2010; Zeng et al.,2014)。木蘭類植物與單子葉植物和真雙子葉植物的分化時間約為166.11百萬年前;五椏果科與薔薇類和菊類植物的分化時間約為124.23百萬年前;薔薇類植物與菊類植物的分化時間約為116.98百萬年前;唇形類植物(Lmiids)與桔梗類植物(Campanulids)的分化時間約為102.37百萬年前。

    3 討論與結論

    長期以來,被子植物的系統(tǒng)發(fā)育關系重建,都是使用質體基因、線粒體基因或少數保守的單拷貝核基因。Yang & Smith(2014)報道了一種基于系統(tǒng)進化樹的同源基因聚類及去旁系同源基因的方法,我們使用此種方法對收集的88種植物核基因集進行聚類,共獲得了多達5 993個one-to-one基因家族,并從這個數據集里面截取各種大小的數據進行進化樹重建,以測定進化樹的穩(wěn)定性。

    獲得比以前更多的核基因家族后,制約系統(tǒng)演化關系構建的另一個因素就是大量的計算資源和計算時間。構建系統(tǒng)進化樹時,一般需要設置bootstrap值(100~1 000)迭代,此步驟非常耗費計算時間。Nguyen et al.(2015)發(fā)表的軟件iqtree,采用ultrafast bootstrap approximation(UFBoot)方法獲得bootstrap值(Von Haeseler et al.,2013),比RAxML軟件的傳統(tǒng)方法,計算速度快10~40倍,并且獲得的bootstrap值更精確。

    我們使用多達5 993個one-to-one基因家族構建的進化樹,與APG IV報道的主要差異為檀香目和石竹目在系統(tǒng)發(fā)育樹中的位置,本研究認為“檀香目和石竹目是薔薇類植物的姊妹群”,而APG IV認為“檀香目和石竹目是菊類植物的姊妹群”??赡茉蛴幸韵聝蓚€:一是基因數目的增多;二是本研究所選88個植物只有一半使用的基因組序列,另一半為轉錄組序列,而轉錄組序列一般存在大量的基因缺失(即未表達基因較多)。

    A. Bootstrap值; B. 基因支持率。紅色表示支持,藍色表示拒絕,紅色星星表示支持率最高的拓撲結構。

    A. Bootstrap value; B. Gene trees support ratio. Red represents support,blue represents rejection,the topology labeled with a red star represents the most support.枝上斜線左邊數字為bootstrap值,右邊數字為基因支持率。

    灰色條紋為分化時間的95%置信區(qū)間,九個化石校準時間為(1)銀杏分化時間為290百萬年至310百萬年前;(2)單子葉植物和真雙子葉植物分化時間為130百萬年至200百萬年前;(3)真雙子葉植物共同祖先(即最早的雙子葉植物化石記錄)為125 百萬年前;(4)山龍眼目(Proteales)的共同祖先為108.8百萬年前;(5)山茱萸目(Cornales)共同祖先為85.8百萬年前;(6)唇形目(Lamiales)共同祖先為44.3百萬年前;(7)葡萄目(Vitales)與其余薔薇類植物間分化時間為105百萬年至115百萬年前;(8)豆目(Fabales)與殼斗目(Fagales)間分化時間為93.5百萬年前;(9)Arabidopsis thaliana與Populus trichocarpa間分化時間為97百萬年至109百萬年前。

    Grey bars are 95% confidence intervals,nine fossil calibration points are as follows:(1) The divergence time of Ginkgo biloba is 290 million years to 310 million years ago; (2) The divergence time of eudicots and monocots is 130 million years to 200 million years ago; (3) The divergence time of eudicots is 125 million years ago; (4) The divergence time of Proteales is 108.8 million years ago; (5) The divergence time of Cornales is 85.8 million years ago; (6) The divergence time of Lamiales is 44.3 million years ago; (7) The divergence time of Vitales from Rosids is 105 million years to 115 million years ago; (8) The divergence time of Fabales and Fagales is 93.5 million years ago; (9) The divergence time of Arabidopsis thaliana and Populus trichocarpa is 97 million years to 109 million years ago.總的來說,本研究不僅進一步確定了被子植物各目間系統(tǒng)發(fā)育關系,而且為“使用更多的基因和計算速度更快的方法構建進化樹”探討了一種可行性策略:即使用Yang & Smith(2014)報道的同源基因聚類及去旁系同源基因方法,獲得大量的one-to-one基因家族,再使用IQ-TREE(串聯(lián)法)和ASTRAL(溯祖法)軟件,能快速精確的計算出進化樹。隨著更多植物基因組的測序和基因聚類及系統(tǒng)發(fā)育關系構建方法的進一步優(yōu)化,被子植物系統(tǒng)發(fā)育關系將越來越精確,例如進一步準確確定檀香目和石竹目在被子植物中與其他進化分支之間的關系。

    參考文獻:

    BOLGERAM,LOHSE M,USADEL B,2014. Trimmomatic:A flexible trimmer for Illumina sequence data[J]. Bioinforma-tics,30(15):2114-2120.

    CALL VB,DILCHER DL,1992. Investigations of angiosperms from the Eocene of southeastern North America:Samaras of Fraxinus wilcoxiana Berry[J]. Rev Palaeobot Palynol,74:249-266.

    CHAW SM,LIU YC,WU YW,et al.,2019. Stout camphor tree genome fills gaps in understanding of flowering plant genome evolution[J]. Natl Plants,5(1):63-73.

    CHEN JH,HAO ZD,GUANG XM,et al.,2019. Liriodendron genome sheds light on angiosperm phylogeny and species-pair differentiation[J]. Nat Plants,5(1):18-25.

    CRANE PR,HERENDEEN PS,1996. Cretaceous floras containing angiosperm flowers and fruits from eastern North America[J]. Rev Palaeobot Palynol,90:319-337.

    EBERSBERGER I,STRAUSS S,VON HAESELER A,2009. HaMStR:Profile hidden markov model based search for orthologs in ESTs[J]. Bmc Evol Biol,9(1):157-157.

    ENDRESS PK,DOYLE JA,2009. Reconstructing the ancestral angiosperm flower and its initial specializations[J]. Am J Bot,96(1):22-66.

    FAWCETT JA,MAERE S,VAN DE PEER Y,2009. Plants with double genomes might have had a better chance to survive the Cretaceous-Tertiary extinction event[J]. Proc Natl Acad Sci USA,106(14):5737-5742.

    FRIIS EM,PEDERSEN KR,SCHNENBERGER J,2006. Normapolles plants:A prominent component of the Cretaceous rosid diversification[J]. Plant Syst Evol,260:107-140.

    FU Q,DIEZ JB,POLE M,et al.,2018. An unexpected noncarpellate epigynous flower from the Jurassic of China[J]. Elife,7:e38827.

    GAO Z,BARRY AT,1989. A review of fossil cycad megasporophylls,with new evidence of Crossozamia pomel and its associated leaves from the lower permian of Taiyuan,China[J]. REV Palaeobot Palynol,60(3-4):205-223.

    GRABHERR MG,HAAS BJ,YASSOUR M,et al.,2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome[J]. Nat Biotechnol,29(7):644-652.

    HE ZL,ZHANG HK,GAO SH,et al.,2016. Evolview v2:An online visualization and management tool for customized and annotated phylogenetic trees[J]. Nucl Acid Res,44(W1):236-241.

    HERENDEEN PS,1995. The enigma of angiosperm origins[J]. Earth-Sci Rev,39(1):253-254.

    KATOH K,STANDLEY DM,2013. MAFFT multiple sequence alignment software version 7:Improvements in performance and usability[J]. Mol Biol Evol,30(4):772-780.

    KUMAR S,STECHER G,SULESKI M,et al.,2017. TimeTree:A Resource for 598 timelines,timetrees,and divergence times[J]. Mol Biol Evol,34:1812-1819.

    LANFEAR R,FRANDSEN PB,WRIGHT AM,et al.,2016. PartitionFinder 2:New methods for selecting partitioned models of evolution formolecular and morphological phylogenetic analyses[J]. Mol Biol Evol,34(3): 772-773.

    LU LM,MAO LF,YANG T,et al.,2018. Evolutionary history of the angiosperm flora of China[J]. Nature,554(1):234-238.

    LI HT,YI TS,GAO LM,et al.,2019. Origin of angiosperms and the puzzle of the Jurassic gap[J]. Nat Plants,5(1):461-470.

    MAGALLON S,2010. Using fossils to break long branches in molecular dating:A comparison of relaxed clocks applied to the origin of angiosperms[J]. Syst Biol,59(4):384-399.

    MOORE MJ,HASSAN N,GITZENDANNER MA,et al.,2011. Phylogenetic analysis of the plastid inverted repeat for 244 species:Insights into deeper-level angiosperm relationships from a long,slowly evolving sequence region[J]. Int J Plant Sci,172(4):541-558.

    MOORE MJ,SOLTIS PS,BELL CD,et al.,2010. Phylogenetic analysis of 83 plastid genes further resolves the early diversification of eudicots[J]. Proc Natl Acad Sci USA,107(10):4623-4628.

    NGUYEN LT,SCHMIDT HA,VON HAESELER A,et al.,2015. IQ-TREE:A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies[J]. Mol Biol Evol,32(1):268-274.

    QIU YL,LI LB,WANG B,et al.,2010. Angiosperm phylogeny inferred from sequences of four mitochondrial genes[J]. JSE,48(6):391-425.

    RUHFEL BR,GITZENDANNER MA,SOLTIS PS,et al.,2014. From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes[J]. Bmc Evol Biol,14(1):23.

    SMITH SA,BEAULIEU JM,DONOGHUE MJ,2010. An uncorrelated relaxed-clock analysis suggests an earlier origin for flowering plants[J]. Proc Natl Acad Sci USA,107(13): 5897-5902.

    SMITH SA,DUNN CW,2008. Phyutility:A phyloinformatics tool for trees,alignments and molecular data[J]. Bioinformatics,24(5):715-716.

    SOLTIS DE,SMITH SA,CELLINESE N,et al.,2011. Angiosperm phylogeny:17 genes,640taxa[J]. Am J Bot,98(4): 704-730.

    STAMATAKIS A,2014. RAxML Version 8:A tool for phylogenetic analysis and post—analysis of large phylogenies[J]. Bioinformatics,30(9):1312-1313.

    TAKAHASHI M,CRANE PR,MANCHESTER SR,2002. Hironoia fusiformis gen. et sp. nov.,a cornalean fruit from the Kamikitaba locality (Upper Cretaceous,Lower Coniacian) in northeastern Japan[J]. J Plant Res,115:463-473.

    THE ANGIOSPERM PHYLOGENY GROUP,2016. An update of the angiosperm phylogeny group classification for the orders and families of flowering plants:APGIV[J]. Bot J Linn Soc,181(1):1-20.

    VANDS,2000. Graph Clustering by Flow Simulation[M]. University of Utrecht.

    VONHAESELER A,MINH BQ,NGUYEN MAT,2013. Ultrafast approximation for phylogenetic bootstrap[J]. Mol Biol Evol,30(5):1188-1195.

    WICKETT NJ,MIRARAB S,NGUYEN N,et al.,2014. Phylotranscriptomic analysis of the origin and early diversification of land plants[J]. Proc Natl Acad Sci USA,111(45): 4859-4868.

    WORBERG A,QUANDT D,BARNISKE AM,et al.,2007. Phylogeny of basal eudicots:Insights from non-coding and rapidly evolving DNA[J]. Org Divers Evol,7(1):55-77.

    YANG Z,2007. PAML 4:Phylogenetic analysis by maximum likelihood[J]. Mol Biol Evol,24:1586-1591.

    YANG Y,MOORE MJ,BROCKINGTON SF,et al.,2015. Dissecting molecular evolution in the highly diverse plant clade caryophyllales using transcriptome sequencing[J]. Mol Biol Evol,32(8):2001-2014.

    YANG Y,SMITH SA,2014. Orthology inference in nonmodel organisms using transcriptomes and low-coverage genomes:Improving accuracy and matrix occupancy for phylogenomics[J]. Mol Biol Evol,31(11):3081-3092.

    ZENG LP,ZHANG N,ZHANG Q,et al.,2017. Resolution of deep eudicot phylogeny and their temporal diversification using nuclear genes from transcriptomic and genomic datasets[J]. New Phytol,214(3):1338-1354.

    ZENG LP,ZHANG Q,SUN RR,et al.,2014. Resolution of deep angiosperm phylogeny using conserved nuclear genes and estimates of early divergence times[J]. Nat Comm,5(1): 4956.

    ZHANG C,SAYYARI E,MIRARAB S,2017. ASTRAL-III:Increased scalability and impacts of contracting low support branches[J]. RECOMB-CG,Springer,Cham:53-75.

    ZHANG N,ZENG LP,SHAN HY,et al.,2012. Highly conserved low-copy nuclear genes as effective markers for phylogenetic analyses in angiosperms[J]. New Phytol,195(4):923-937.

    (責任編輯 周翠鳴)

    精品人妻偷拍中文字幕| 国产中年淑女户外野战色| 亚洲美女搞黄在线观看| 欧美不卡视频在线免费观看| 国产一区二区在线观看日韩| 搡老妇女老女人老熟妇| 日韩一区二区三区影片| 欧美日本视频| av网站免费在线观看视频 | 麻豆乱淫一区二区| 亚洲av中文字字幕乱码综合| 国产高潮美女av| 在线免费观看的www视频| 91久久精品国产一区二区成人| 99热这里只有是精品在线观看| 午夜精品一区二区三区免费看| 精品国产一区二区三区久久久樱花 | 视频中文字幕在线观看| 如何舔出高潮| 91狼人影院| 久久热精品热| 久久这里只有精品中国| 777米奇影视久久| 日本av手机在线免费观看| 最近2019中文字幕mv第一页| 80岁老熟妇乱子伦牲交| 汤姆久久久久久久影院中文字幕 | 在线观看人妻少妇| 国产一级毛片七仙女欲春2| 一级毛片 在线播放| 毛片一级片免费看久久久久| 又爽又黄a免费视频| 女人久久www免费人成看片| 91av网一区二区| 一级二级三级毛片免费看| 97精品久久久久久久久久精品| 成人综合一区亚洲| 纵有疾风起免费观看全集完整版 | 亚洲精品一区蜜桃| 永久网站在线| 日韩三级伦理在线观看| 国产精品美女特级片免费视频播放器| 国产伦理片在线播放av一区| 国内揄拍国产精品人妻在线| 亚洲精品国产av成人精品| av线在线观看网站| 肉色欧美久久久久久久蜜桃 | 久久国内精品自在自线图片| 99久国产av精品| 国产色爽女视频免费观看| av在线亚洲专区| www.av在线官网国产| av在线蜜桃| 91av网一区二区| 亚洲国产欧美在线一区| 搡老妇女老女人老熟妇| 午夜福利视频1000在线观看| 三级国产精品片| 久久久久久久久久久丰满| 少妇的逼水好多| 女人久久www免费人成看片| 在线观看av片永久免费下载| 日韩伦理黄色片| 综合色av麻豆| 最近手机中文字幕大全| 建设人人有责人人尽责人人享有的 | 欧美性猛交╳xxx乱大交人| 成人性生交大片免费视频hd| 久久精品国产亚洲av天美| 亚洲欧美成人精品一区二区| 不卡视频在线观看欧美| 欧美三级亚洲精品| 日韩欧美精品v在线| 国产精品国产三级国产av玫瑰| 久久综合国产亚洲精品| 一级毛片我不卡| 久久精品久久久久久久性| 国产 一区 欧美 日韩| 欧美激情国产日韩精品一区| 精品一区二区三区人妻视频| 激情 狠狠 欧美| 久久精品久久久久久久性| 一个人看视频在线观看www免费| 欧美日韩在线观看h| 一个人免费在线观看电影| 亚洲国产成人一精品久久久| 大陆偷拍与自拍| 久久久久久久久久黄片| 免费在线观看成人毛片| 亚洲最大成人中文| 免费av毛片视频| 国产亚洲精品久久久com| 免费少妇av软件| av网站免费在线观看视频 | 国产一区有黄有色的免费视频 | 好男人视频免费观看在线| videos熟女内射| 亚洲欧美一区二区三区黑人 | 白带黄色成豆腐渣| 日韩欧美三级三区| 久久人人爽人人片av| 欧美人与善性xxx| 国产不卡一卡二| 国产精品久久视频播放| 一级毛片aaaaaa免费看小| 国产精品综合久久久久久久免费| 国产综合精华液| 嘟嘟电影网在线观看| 免费无遮挡裸体视频| 精品一区二区免费观看| 国产人妻一区二区三区在| 免费看不卡的av| 国产精品福利在线免费观看| 卡戴珊不雅视频在线播放| 亚洲精品一区蜜桃| 亚洲精品亚洲一区二区| 一二三四中文在线观看免费高清| 亚洲在线观看片| 久久久亚洲精品成人影院| 听说在线观看完整版免费高清| 99久国产av精品| av福利片在线观看| xxx大片免费视频| 久久久久精品性色| 国产一区二区三区综合在线观看 | 亚洲熟女精品中文字幕| 91久久精品国产一区二区成人| 中文在线观看免费www的网站| 国内精品美女久久久久久| 久久精品综合一区二区三区| 亚洲久久久久久中文字幕| 夜夜看夜夜爽夜夜摸| 精品久久国产蜜桃| 欧美极品一区二区三区四区| 91久久精品电影网| 91在线精品国自产拍蜜月| 狠狠精品人妻久久久久久综合| 国产有黄有色有爽视频| 日韩av免费高清视频| 国产中年淑女户外野战色| 伊人久久精品亚洲午夜| 亚洲综合色惰| 久久人人爽人人片av| 久久鲁丝午夜福利片| 色综合亚洲欧美另类图片| 亚洲国产色片| 欧美不卡视频在线免费观看| 久久久久久国产a免费观看| 国产亚洲一区二区精品| 肉色欧美久久久久久久蜜桃 | 国产精品国产三级国产专区5o| 人妻制服诱惑在线中文字幕| 精华霜和精华液先用哪个| 成人漫画全彩无遮挡| 国产精品嫩草影院av在线观看| 又黄又爽又刺激的免费视频.| 国产精品综合久久久久久久免费| 亚洲欧美日韩卡通动漫| 午夜精品一区二区三区免费看| 国产亚洲5aaaaa淫片| 在线观看免费高清a一片| av在线老鸭窝| eeuss影院久久| 99久久精品国产国产毛片| 精品国内亚洲2022精品成人| av在线蜜桃| 又粗又硬又长又爽又黄的视频| 最近中文字幕高清免费大全6| 亚洲国产色片| 国产成人a区在线观看| 国产男女超爽视频在线观看| 免费播放大片免费观看视频在线观看| 午夜亚洲福利在线播放| 久久精品夜色国产| 美女被艹到高潮喷水动态| 永久网站在线| 亚州av有码| 一级毛片我不卡| 亚洲国产精品国产精品| 亚洲成色77777| 九九久久精品国产亚洲av麻豆| 免费观看无遮挡的男女| 晚上一个人看的免费电影| 久久久久免费精品人妻一区二区| 日产精品乱码卡一卡2卡三| 久久久久久国产a免费观看| 能在线免费观看的黄片| 成年版毛片免费区| 97在线视频观看| a级毛色黄片| 色5月婷婷丁香| 少妇人妻一区二区三区视频| 国产精品女同一区二区软件| 日韩一区二区三区影片| 国产亚洲精品av在线| 少妇丰满av| 好男人在线观看高清免费视频| 亚洲欧美成人精品一区二区| 亚州av有码| 国产69精品久久久久777片| 色网站视频免费| 最近的中文字幕免费完整| 看免费成人av毛片| 亚洲精品色激情综合| 亚洲一区高清亚洲精品| 国产免费一级a男人的天堂| 日本免费a在线| 丝袜喷水一区| 黄片wwwwww| 精品一区二区三区视频在线| 国产男女超爽视频在线观看| 乱人视频在线观看| 免费人成在线观看视频色| 五月伊人婷婷丁香| a级毛色黄片| 精品国产露脸久久av麻豆 | videos熟女内射| 亚洲精品,欧美精品| 亚洲精品国产av成人精品| 国产伦在线观看视频一区| 最近最新中文字幕免费大全7| 99热这里只有是精品在线观看| 国产有黄有色有爽视频| 在线观看免费高清a一片| 少妇的逼水好多| 九九久久精品国产亚洲av麻豆| 免费黄色在线免费观看| 欧美xxⅹ黑人| 亚洲高清免费不卡视频| 欧美激情国产日韩精品一区| av一本久久久久| 高清在线视频一区二区三区| 99视频精品全部免费 在线| 日韩欧美三级三区| 简卡轻食公司| 欧美区成人在线视频| 我要看日韩黄色一级片| 亚洲在线自拍视频| 国产精品av视频在线免费观看| 国产黄色视频一区二区在线观看| 午夜福利在线观看免费完整高清在| 亚洲精品久久午夜乱码| 中文资源天堂在线| 床上黄色一级片| 99热这里只有精品一区| 久久99蜜桃精品久久| 99久国产av精品国产电影| 一个人看视频在线观看www免费| 国产av不卡久久| 直男gayav资源| 亚洲真实伦在线观看| 禁无遮挡网站| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交╳xxx乱大交人| 久久久久久久亚洲中文字幕| 亚洲美女视频黄频| 欧美成人精品欧美一级黄| 国产爱豆传媒在线观看| 国产亚洲精品久久久com| 九草在线视频观看| 偷拍熟女少妇极品色| 国产探花在线观看一区二区| 久久久久久久久久黄片| 国产在视频线在精品| 亚洲自拍偷在线| 五月天丁香电影| 国产精品1区2区在线观看.| 国产综合精华液| 可以在线观看毛片的网站| 国产精品久久久久久精品电影| 欧美丝袜亚洲另类| 久久热精品热| 男插女下体视频免费在线播放| 看免费成人av毛片| 嫩草影院新地址| 午夜视频国产福利| 一级片'在线观看视频| 国产91av在线免费观看| 亚洲精品成人久久久久久| 在线免费观看不下载黄p国产| 人体艺术视频欧美日本| 午夜福利视频精品| 国产毛片a区久久久久| 亚洲最大成人av| 中文字幕久久专区| 精品午夜福利在线看| 丝瓜视频免费看黄片| 国产精品无大码| 亚洲人与动物交配视频| 在线观看av片永久免费下载| 简卡轻食公司| 边亲边吃奶的免费视频| 亚洲精品日本国产第一区| 成人亚洲欧美一区二区av| 午夜亚洲福利在线播放| 欧美人与善性xxx| 男女国产视频网站| 99久久精品一区二区三区| 美女黄网站色视频| 最近最新中文字幕大全电影3| 欧美日韩精品成人综合77777| 久久人人爽人人片av| 国产高潮美女av| 最近手机中文字幕大全| 国产有黄有色有爽视频| 亚洲va在线va天堂va国产| 少妇高潮的动态图| 亚洲aⅴ乱码一区二区在线播放| 亚洲性久久影院| 久久精品综合一区二区三区| 国产伦精品一区二区三区四那| 亚洲精品乱久久久久久| 久久久久久久亚洲中文字幕| 国产成人91sexporn| 免费黄色在线免费观看| 国产高清不卡午夜福利| 国产一区二区亚洲精品在线观看| 干丝袜人妻中文字幕| 亚洲最大成人中文| 亚洲综合色惰| 国产黄色视频一区二区在线观看| 国产成人freesex在线| 高清日韩中文字幕在线| 最近2019中文字幕mv第一页| 亚洲最大成人中文| 国产男人的电影天堂91| 国产成人aa在线观看| 免费大片黄手机在线观看| 我的女老师完整版在线观看| 免费黄色在线免费观看| h日本视频在线播放| 欧美另类一区| 蜜桃久久精品国产亚洲av| 边亲边吃奶的免费视频| 亚洲aⅴ乱码一区二区在线播放| 免费人成在线观看视频色| 少妇被粗大猛烈的视频| av国产久精品久网站免费入址| 国产精品久久视频播放| av在线老鸭窝| 精品国产三级普通话版| 成人无遮挡网站| 联通29元200g的流量卡| 人体艺术视频欧美日本| 乱人视频在线观看| 天堂影院成人在线观看| 男女啪啪激烈高潮av片| 国产精品一区二区三区四区久久| 国产男人的电影天堂91| 国产一级毛片在线| 国产永久视频网站| 日本一二三区视频观看| 精品酒店卫生间| 久久人人爽人人爽人人片va| 亚洲精品第二区| 精品酒店卫生间| 少妇的逼好多水| 国产91av在线免费观看| 亚洲精品自拍成人| 亚洲av一区综合| 天堂影院成人在线观看| 久久热精品热| 日日摸夜夜添夜夜爱| videos熟女内射| 亚洲精品,欧美精品| 国产在视频线精品| 久久久久久久久久成人| 高清欧美精品videossex| 日韩成人伦理影院| 日韩大片免费观看网站| 女人被狂操c到高潮| 国产免费一级a男人的天堂| 日日摸夜夜添夜夜爱| 天堂av国产一区二区熟女人妻| 日韩精品有码人妻一区| 欧美xxxx性猛交bbbb| 在线天堂最新版资源| 久久综合国产亚洲精品| 国产探花在线观看一区二区| 亚洲欧美日韩无卡精品| 天天躁日日操中文字幕| 99久久精品一区二区三区| 国产伦精品一区二区三区视频9| 久久人人爽人人片av| 日本免费在线观看一区| 校园人妻丝袜中文字幕| 国产日韩欧美在线精品| 亚洲成人中文字幕在线播放| av女优亚洲男人天堂| 欧美不卡视频在线免费观看| 久久久久久国产a免费观看| 天堂影院成人在线观看| 国产午夜精品久久久久久一区二区三区| 春色校园在线视频观看| 黄色配什么色好看| 色综合站精品国产| 老女人水多毛片| 日本午夜av视频| 午夜激情福利司机影院| 国产成人精品久久久久久| 成人午夜高清在线视频| .国产精品久久| 亚洲aⅴ乱码一区二区在线播放| 大香蕉久久网| av国产久精品久网站免费入址| 99视频精品全部免费 在线| 国产亚洲av嫩草精品影院| 我的老师免费观看完整版| 国产人妻一区二区三区在| 亚洲精品亚洲一区二区| 久久草成人影院| 少妇丰满av| 又爽又黄a免费视频| 黑人高潮一二区| 国产探花在线观看一区二区| 嫩草影院精品99| 伦精品一区二区三区| 欧美激情在线99| 国产精品1区2区在线观看.| 亚洲va在线va天堂va国产| 国内精品美女久久久久久| 国产成人aa在线观看| 免费高清在线观看视频在线观看| 欧美日韩亚洲高清精品| 精品一区二区免费观看| 国产探花在线观看一区二区| 亚洲av电影在线观看一区二区三区 | 麻豆精品久久久久久蜜桃| 亚洲av二区三区四区| 国产成人精品福利久久| 淫秽高清视频在线观看| 欧美97在线视频| 国产不卡一卡二| 干丝袜人妻中文字幕| av在线观看视频网站免费| 午夜精品在线福利| 亚洲经典国产精华液单| 国产精品三级大全| 久久精品国产鲁丝片午夜精品| 亚洲精品视频女| 在线观看免费高清a一片| 日韩欧美一区视频在线观看 | 亚洲综合色惰| 亚洲精品国产av蜜桃| 日产精品乱码卡一卡2卡三| 中文乱码字字幕精品一区二区三区 | 亚洲国产精品sss在线观看| 狂野欧美白嫩少妇大欣赏| 日本黄大片高清| 最后的刺客免费高清国语| 成人高潮视频无遮挡免费网站| 美女xxoo啪啪120秒动态图| 草草在线视频免费看| 成人二区视频| 1000部很黄的大片| 久久精品久久久久久噜噜老黄| 色吧在线观看| 成人美女网站在线观看视频| 亚洲av成人av| 亚洲精品中文字幕在线视频 | 身体一侧抽搐| 最后的刺客免费高清国语| 十八禁网站网址无遮挡 | 国产片特级美女逼逼视频| 最近中文字幕2019免费版| 丝袜美腿在线中文| 亚洲精华国产精华液的使用体验| 午夜激情久久久久久久| 丝袜喷水一区| 亚洲综合色惰| 99久久人妻综合| 日本色播在线视频| 国产精品一二三区在线看| 久久久久久九九精品二区国产| 高清欧美精品videossex| 美女脱内裤让男人舔精品视频| 91久久精品国产一区二区成人| 国产探花极品一区二区| 综合色丁香网| 国产av码专区亚洲av| 国产黄色视频一区二区在线观看| 少妇被粗大猛烈的视频| 成人综合一区亚洲| 搡老妇女老女人老熟妇| 成年女人看的毛片在线观看| 建设人人有责人人尽责人人享有的 | 观看免费一级毛片| 搡老乐熟女国产| 亚洲色图av天堂| 中文乱码字字幕精品一区二区三区 | 国产精品熟女久久久久浪| 蜜臀久久99精品久久宅男| 久久精品国产亚洲av天美| 九色成人免费人妻av| 国产成人一区二区在线| 午夜免费激情av| 婷婷色麻豆天堂久久| 蜜桃久久精品国产亚洲av| 久久久国产一区二区| 91狼人影院| 日韩欧美精品免费久久| 男女边摸边吃奶| 内地一区二区视频在线| 夫妻午夜视频| 三级经典国产精品| 熟女人妻精品中文字幕| 又大又黄又爽视频免费| 免费人成在线观看视频色| 日本爱情动作片www.在线观看| 免费观看av网站的网址| 成人鲁丝片一二三区免费| 国产 亚洲一区二区三区 | 久久这里有精品视频免费| 免费观看性生交大片5| 免费观看在线日韩| 亚洲一区高清亚洲精品| 免费不卡的大黄色大毛片视频在线观看 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费大片黄手机在线观看| 国内精品宾馆在线| 97人妻精品一区二区三区麻豆| 一级爰片在线观看| 国产综合精华液| 肉色欧美久久久久久久蜜桃 | 国产成人freesex在线| 亚洲最大成人手机在线| 黄片无遮挡物在线观看| 一级爰片在线观看| 日日干狠狠操夜夜爽| 日本-黄色视频高清免费观看| 啦啦啦啦在线视频资源| 日韩制服骚丝袜av| 欧美潮喷喷水| 久久韩国三级中文字幕| 91午夜精品亚洲一区二区三区| 高清毛片免费看| 国产一区二区亚洲精品在线观看| h日本视频在线播放| 国产人妻一区二区三区在| 男女下面进入的视频免费午夜| 国产色爽女视频免费观看| 在线观看美女被高潮喷水网站| 夜夜看夜夜爽夜夜摸| 亚洲成人一二三区av| 久久久久久久久大av| 尾随美女入室| 大香蕉97超碰在线| 五月伊人婷婷丁香| 午夜免费观看性视频| 啦啦啦中文免费视频观看日本| 国产有黄有色有爽视频| 国产单亲对白刺激| 少妇的逼好多水| 白带黄色成豆腐渣| 亚洲精品成人av观看孕妇| 七月丁香在线播放| 少妇的逼水好多| 日韩欧美精品免费久久| 中文欧美无线码| 午夜精品一区二区三区免费看| 一区二区三区四区激情视频| 日本-黄色视频高清免费观看| 成年免费大片在线观看| 亚洲av福利一区| 精品酒店卫生间| 久久精品夜色国产| 一级毛片 在线播放| 欧美激情国产日韩精品一区| 久久午夜福利片| 午夜福利高清视频| 尤物成人国产欧美一区二区三区| 亚洲av电影在线观看一区二区三区 | 久久精品人妻少妇| 又爽又黄a免费视频| 免费大片18禁| 国产高潮美女av| 色尼玛亚洲综合影院| 欧美性猛交╳xxx乱大交人| 国产午夜福利久久久久久| 成人亚洲精品av一区二区| 激情五月婷婷亚洲| 你懂的网址亚洲精品在线观看| 国产v大片淫在线免费观看| 精品久久久久久久久av| 国产高清不卡午夜福利| 国产单亲对白刺激| 亚洲国产欧美人成| 国产熟女欧美一区二区| 久久久久性生活片| 看非洲黑人一级黄片| 国产精品久久久久久久久免| 亚洲国产精品专区欧美| 精品人妻熟女av久视频| www.av在线官网国产| 91在线精品国自产拍蜜月| 精品不卡国产一区二区三区| 国产精品美女特级片免费视频播放器| 91av网一区二区| 观看免费一级毛片| 麻豆国产97在线/欧美| 伦精品一区二区三区| 国模一区二区三区四区视频| 两个人的视频大全免费| 麻豆精品久久久久久蜜桃| 久久精品夜夜夜夜夜久久蜜豆| 免费人成在线观看视频色| 久久97久久精品| 大话2 男鬼变身卡| 久久草成人影院| 国产男人的电影天堂91| 日日摸夜夜添夜夜添av毛片| 好男人视频免费观看在线| 亚洲精品视频女| 高清在线视频一区二区三区| 国产精品爽爽va在线观看网站|