• <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| 亚洲人成电影观看| 又大又爽又粗| 亚洲熟女精品中文字幕| 男人的好看免费观看在线视频 | 又黄又粗又硬又大视频| 欧美日韩福利视频一区二区| 美女 人体艺术 gogo| 精品久久蜜臀av无| 9色porny在线观看| 91在线观看av| 一二三四在线观看免费中文在| 欧美 日韩 精品 国产| 天天躁夜夜躁狠狠躁躁| 精品国产超薄肉色丝袜足j| 国产在线观看jvid| 90打野战视频偷拍视频| 超色免费av| e午夜精品久久久久久久| 大片电影免费在线观看免费| 成人18禁在线播放| 91精品国产国语对白视频| 久久精品国产综合久久久| 欧美精品亚洲一区二区| 国产精品综合久久久久久久免费 | av中文乱码字幕在线| 日本五十路高清| 18禁国产床啪视频网站| 19禁男女啪啪无遮挡网站| 免费高清在线观看日韩| 亚洲av熟女| 欧美人与性动交α欧美精品济南到| 国产国语露脸激情在线看| 中文字幕最新亚洲高清| 日韩免费高清中文字幕av| 下体分泌物呈黄色| 一级毛片女人18水好多| 亚洲自偷自拍图片 自拍| 欧美+亚洲+日韩+国产| 国产免费av片在线观看野外av| 一区二区日韩欧美中文字幕| 欧美av亚洲av综合av国产av| 婷婷精品国产亚洲av在线 | 在线观看免费日韩欧美大片| 国产精品美女特级片免费视频播放器 | 国产真人三级小视频在线观看| 精品国产美女av久久久久小说| 老司机影院毛片| 亚洲av欧美aⅴ国产| 久9热在线精品视频| 国产黄色免费在线视频| 久久性视频一级片| 国产欧美日韩一区二区精品| 国产成人av教育| 精品国产乱子伦一区二区三区| 不卡一级毛片| 五月开心婷婷网| 18禁黄网站禁片午夜丰满| 日日爽夜夜爽网站| 欧美日韩亚洲高清精品| 日本黄色日本黄色录像| 国产国语露脸激情在线看| 窝窝影院91人妻| 中文字幕制服av| 日韩三级视频一区二区三区| 久99久视频精品免费| 免费在线观看日本一区| 日韩制服丝袜自拍偷拍| 久久人妻熟女aⅴ| 成年女人毛片免费观看观看9 | 啦啦啦 在线观看视频| 中文字幕人妻熟女乱码| 咕卡用的链子| 日韩精品免费视频一区二区三区| av福利片在线| 人人妻人人添人人爽欧美一区卜| 国产精品久久视频播放| 精品熟女少妇八av免费久了| 一级片'在线观看视频| 啦啦啦在线免费观看视频4| 久久久久国产精品人妻aⅴ院 | 日韩成人在线观看一区二区三区| 亚洲精品国产精品久久久不卡| 在线观看午夜福利视频| 亚洲五月婷婷丁香| 久久久久国产一级毛片高清牌| 国产成人免费无遮挡视频| 777久久人妻少妇嫩草av网站| 欧美日韩精品网址| 免费黄频网站在线观看国产| 99在线人妻在线中文字幕 | 午夜日韩欧美国产| 久久精品亚洲精品国产色婷小说| 日本五十路高清| 日韩 欧美 亚洲 中文字幕| 美女视频免费永久观看网站| 亚洲国产欧美日韩在线播放| 久热这里只有精品99| 天堂√8在线中文| 在线观看午夜福利视频| 日韩欧美在线二视频 | 黄色成人免费大全| 精品一区二区三卡| 高清黄色对白视频在线免费看| 人人妻人人澡人人看| 国产97色在线日韩免费| 脱女人内裤的视频| 不卡一级毛片| 捣出白浆h1v1| 男女床上黄色一级片免费看| 国产伦人伦偷精品视频| 青草久久国产| 免费在线观看黄色视频的| 国内久久婷婷六月综合欲色啪| 亚洲欧美日韩另类电影网站| 岛国毛片在线播放| 人妻久久中文字幕网| 国产乱人伦免费视频| 久热这里只有精品99| 亚洲五月婷婷丁香| 日日爽夜夜爽网站| 99re在线观看精品视频| 成年人黄色毛片网站| 亚洲 欧美一区二区三区| 国产在线一区二区三区精| 香蕉丝袜av| 欧美另类亚洲清纯唯美| 亚洲精品中文字幕在线视频| 叶爱在线成人免费视频播放| 99久久综合精品五月天人人| 国产视频一区二区在线看| 日本精品一区二区三区蜜桃| 国产成人欧美在线观看 | 手机成人av网站| 国产淫语在线视频| 亚洲人成电影观看| 亚洲国产欧美网| 亚洲欧美色中文字幕在线| 久久亚洲真实| 人妻久久中文字幕网| 久久久久国产一级毛片高清牌| 日韩成人在线观看一区二区三区| 精品国产美女av久久久久小说| 国产精品一区二区在线观看99| 国产av又大| 国产精品一区二区在线不卡| 美女国产高潮福利片在线看| 99久久人妻综合| 啦啦啦免费观看视频1| 久久精品熟女亚洲av麻豆精品| av网站免费在线观看视频| 黑人猛操日本美女一级片| 美女视频免费永久观看网站| 91av网站免费观看| 高清在线国产一区| 天天躁夜夜躁狠狠躁躁| а√天堂www在线а√下载 | 国产蜜桃级精品一区二区三区 | 精品人妻1区二区| 天天操日日干夜夜撸| 91大片在线观看| 国产精品综合久久久久久久免费 | 亚洲伊人色综图| 狠狠狠狠99中文字幕| 久久国产精品人妻蜜桃| 日韩制服丝袜自拍偷拍| 交换朋友夫妻互换小说| 看免费av毛片| 国产高清激情床上av| 久久午夜亚洲精品久久| a级毛片黄视频| 天堂动漫精品| 欧美黄色淫秽网站| 亚洲精品国产区一区二| 99国产精品一区二区蜜桃av | 人人妻人人澡人人爽人人夜夜| 亚洲精品乱久久久久久| 欧美日韩黄片免| 黑丝袜美女国产一区| √禁漫天堂资源中文www| 亚洲精品中文字幕在线视频| 欧美精品人与动牲交sv欧美| 在线国产一区二区在线| 99热只有精品国产| 中出人妻视频一区二区| 999久久久国产精品视频| 国产伦人伦偷精品视频| 亚洲在线自拍视频| 午夜福利乱码中文字幕| 免费不卡黄色视频| 丰满饥渴人妻一区二区三| 精品国产一区二区三区久久久樱花| 午夜久久久在线观看| 美女视频免费永久观看网站| 中亚洲国语对白在线视频| 精品国内亚洲2022精品成人 | 免费在线观看黄色视频的| 成年人午夜在线观看视频| 捣出白浆h1v1| 国产人伦9x9x在线观看| 一级a爱片免费观看的视频| 欧美乱妇无乱码| 午夜福利影视在线免费观看| 精品人妻熟女毛片av久久网站| 夜夜夜夜夜久久久久| 少妇 在线观看| 人人妻人人澡人人爽人人夜夜| 国产精品 欧美亚洲| 亚洲 国产 在线| 国产精品 国内视频| 亚洲精品av麻豆狂野| 国产一区在线观看成人免费| 久久久精品国产亚洲av高清涩受| 亚洲专区国产一区二区| 久久精品国产亚洲av高清一级| 午夜免费鲁丝| 国产xxxxx性猛交| 99精品久久久久人妻精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲少妇的诱惑av| 午夜免费观看网址| 国产视频一区二区在线看| 国产主播在线观看一区二区| 成年人黄色毛片网站| 午夜免费鲁丝| 日韩制服丝袜自拍偷拍| 免费女性裸体啪啪无遮挡网站| www.精华液| 国产精品久久久久久精品古装| 久久精品人人爽人人爽视色| 亚洲一卡2卡3卡4卡5卡精品中文| 脱女人内裤的视频| 亚洲国产欧美日韩在线播放| 十八禁网站免费在线| 高清毛片免费观看视频网站 | 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜制服| 人妻丰满熟妇av一区二区三区 | 精品无人区乱码1区二区| 下体分泌物呈黄色| 亚洲人成77777在线视频| 最新美女视频免费是黄的| 女人高潮潮喷娇喘18禁视频| 欧美色视频一区免费| 老司机深夜福利视频在线观看| tube8黄色片| av超薄肉色丝袜交足视频| 人人妻人人添人人爽欧美一区卜| 亚洲av成人不卡在线观看播放网| 亚洲av欧美aⅴ国产| 精品卡一卡二卡四卡免费| 91麻豆精品激情在线观看国产 | 久久精品亚洲精品国产色婷小说| 19禁男女啪啪无遮挡网站| 国内久久婷婷六月综合欲色啪| 新久久久久国产一级毛片| 欧美色视频一区免费| 97人妻天天添夜夜摸| 日韩视频一区二区在线观看| 一级毛片精品| 国产欧美日韩精品亚洲av| 国产男女内射视频| 啦啦啦免费观看视频1| 热99re8久久精品国产| 丰满饥渴人妻一区二区三| 国产成人精品久久二区二区91| 亚洲中文日韩欧美视频| 一二三四社区在线视频社区8| 亚洲综合色网址| 伊人久久大香线蕉亚洲五| 香蕉丝袜av| 亚洲精品乱久久久久久| 国产欧美亚洲国产| 久久久国产欧美日韩av| 69av精品久久久久久| 最新美女视频免费是黄的| 999精品在线视频| 91成年电影在线观看| aaaaa片日本免费| 99国产精品一区二区三区| av视频免费观看在线观看| 99riav亚洲国产免费| 丁香六月欧美| 天天躁狠狠躁夜夜躁狠狠躁| 黄色视频,在线免费观看| 精品国产一区二区久久| 侵犯人妻中文字幕一二三四区| 一夜夜www| 欧美久久黑人一区二区| 久久婷婷成人综合色麻豆| 91成年电影在线观看| 色婷婷av一区二区三区视频| 色综合婷婷激情| av中文乱码字幕在线| 18禁美女被吸乳视频| 免费久久久久久久精品成人欧美视频| 女同久久另类99精品国产91| 色94色欧美一区二区| 丰满迷人的少妇在线观看| 亚洲精品久久成人aⅴ小说| 久久久国产成人精品二区 | 村上凉子中文字幕在线| 99香蕉大伊视频| 午夜两性在线视频| 亚洲一区二区三区欧美精品| 黄色成人免费大全| 久久精品91无色码中文字幕| 国产在视频线精品| 99精国产麻豆久久婷婷| 一区二区三区国产精品乱码| 亚洲av日韩精品久久久久久密| 精品国产乱子伦一区二区三区| 国产一卡二卡三卡精品| 国内久久婷婷六月综合欲色啪| 啦啦啦免费观看视频1| 国产亚洲欧美98| 午夜视频精品福利| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区三区在线| 在线观看免费午夜福利视频| 丰满饥渴人妻一区二区三| 久久精品91无色码中文字幕| 一二三四在线观看免费中文在| 最近最新中文字幕大全免费视频| 亚洲自偷自拍图片 自拍| 大陆偷拍与自拍| 窝窝影院91人妻| 高清视频免费观看一区二区| 久久久久久久久免费视频了| 91精品三级在线观看| 丰满迷人的少妇在线观看| 国产精品国产av在线观看| av网站在线播放免费| 国产成人系列免费观看| 色综合欧美亚洲国产小说| 精品亚洲成a人片在线观看| 国产成人av教育| 三上悠亚av全集在线观看| 男人的好看免费观看在线视频 | 午夜福利乱码中文字幕| 国产精品 欧美亚洲| 久久精品熟女亚洲av麻豆精品| 可以免费在线观看a视频的电影网站| 日韩欧美国产一区二区入口| 极品人妻少妇av视频| 午夜精品国产一区二区电影| 日韩有码中文字幕| 黑人巨大精品欧美一区二区蜜桃| 日韩一卡2卡3卡4卡2021年| 午夜精品国产一区二区电影| 精品无人区乱码1区二区| 好男人电影高清在线观看| 在线观看免费午夜福利视频| 丝瓜视频免费看黄片| 国产精品美女特级片免费视频播放器 | 老司机午夜十八禁免费视频| 极品少妇高潮喷水抽搐| 免费在线观看黄色视频的| 中文字幕制服av| 在线观看日韩欧美| 嫩草影视91久久| 三级毛片av免费| 亚洲 国产 在线| 美女 人体艺术 gogo| 久久婷婷成人综合色麻豆| 亚洲一码二码三码区别大吗| 99久久精品国产亚洲精品| 国产精品美女特级片免费视频播放器 | 国产蜜桃级精品一区二区三区 | 少妇的丰满在线观看| 美女视频免费永久观看网站| 一级片免费观看大全| 韩国av一区二区三区四区| 亚洲第一青青草原| 中文字幕色久视频| 久久久久久人人人人人| 欧美精品人与动牲交sv欧美| 中文字幕av电影在线播放| 国产精品免费一区二区三区在线 | 老熟妇乱子伦视频在线观看| 日韩视频一区二区在线观看| 99精品在免费线老司机午夜| 欧美激情 高清一区二区三区| 欧美中文综合在线视频| 免费看a级黄色片| 在线观看午夜福利视频| 在线十欧美十亚洲十日本专区| 99在线人妻在线中文字幕 | 在线视频色国产色| 欧美色视频一区免费| 亚洲三区欧美一区| 18禁美女被吸乳视频| 亚洲在线自拍视频| 国产精品久久久av美女十八| 69精品国产乱码久久久| 色老头精品视频在线观看| 人妻一区二区av| tocl精华| 亚洲自偷自拍图片 自拍| 亚洲,欧美精品.| aaaaa片日本免费| 侵犯人妻中文字幕一二三四区| 亚洲人成77777在线视频| 欧美精品av麻豆av| 日韩欧美国产一区二区入口| 在线观看午夜福利视频| 咕卡用的链子| 99精品欧美一区二区三区四区| 精品午夜福利视频在线观看一区| 丝袜在线中文字幕| 69av精品久久久久久| 精品高清国产在线一区| 99热只有精品国产| 成熟少妇高潮喷水视频| 日韩 欧美 亚洲 中文字幕| svipshipincom国产片| 国产精品国产av在线观看| 两个人看的免费小视频| 久久久久久久午夜电影 | 丝瓜视频免费看黄片| 国产区一区二久久| 久久亚洲精品不卡| 我的亚洲天堂| 国产精品二区激情视频| 超碰成人久久| 亚洲,欧美精品.| 欧美日韩亚洲国产一区二区在线观看 | 80岁老熟妇乱子伦牲交| 久久午夜亚洲精品久久| 老熟妇乱子伦视频在线观看| videosex国产| 久久国产精品人妻蜜桃| 精品人妻在线不人妻| 国产精品久久电影中文字幕 | 久久久久久久午夜电影 | 亚洲色图综合在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美黄色淫秽网站| 人人妻,人人澡人人爽秒播| 成人黄色视频免费在线看| 黑人巨大精品欧美一区二区mp4| 精品久久久久久久毛片微露脸| 精品久久久久久久久久免费视频 | 色在线成人网| 亚洲国产中文字幕在线视频| 美女高潮喷水抽搐中文字幕| 色综合婷婷激情| 亚洲七黄色美女视频| 一级黄色大片毛片| 亚洲人成电影免费在线| 欧美精品高潮呻吟av久久| 法律面前人人平等表现在哪些方面| 黄片小视频在线播放| 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| 曰老女人黄片| 亚洲熟女精品中文字幕| 涩涩av久久男人的天堂| 黑人猛操日本美女一级片| 99国产精品免费福利视频| 在线观看免费高清a一片| 国产精品av久久久久免费| 亚洲色图av天堂| 成年人黄色毛片网站| av线在线观看网站| 18禁裸乳无遮挡动漫免费视频| av在线播放免费不卡| 9色porny在线观看| 久久国产精品大桥未久av| 国产精品国产高清国产av | www.自偷自拍.com| 国产亚洲欧美精品永久| 中文字幕人妻丝袜一区二区| 王馨瑶露胸无遮挡在线观看| 亚洲国产欧美一区二区综合| 成人av一区二区三区在线看| 欧美黄色淫秽网站| 亚洲精品一二三| 在线观看免费高清a一片| 亚洲精品中文字幕一二三四区| 欧美成人午夜精品| x7x7x7水蜜桃| 久久精品熟女亚洲av麻豆精品| 久久久久久久久免费视频了| 久久中文字幕人妻熟女| 亚洲欧美日韩另类电影网站| 三级毛片av免费| 国产淫语在线视频| 免费观看精品视频网站| 免费女性裸体啪啪无遮挡网站| 在线观看免费视频日本深夜| 国产极品粉嫩免费观看在线| 精品视频人人做人人爽| 我的亚洲天堂| 黄频高清免费视频| 久久久久久免费高清国产稀缺| 免费在线观看亚洲国产| 午夜福利一区二区在线看| 亚洲欧美日韩高清在线视频| 免费女性裸体啪啪无遮挡网站| 久久香蕉精品热| 国产麻豆69| 国产淫语在线视频| 窝窝影院91人妻| 日本a在线网址| 色婷婷久久久亚洲欧美| 中文字幕另类日韩欧美亚洲嫩草| 久久精品91无色码中文字幕| 亚洲欧美一区二区三区久久| 亚洲精品美女久久av网站| 精品一区二区三区av网在线观看| 少妇被粗大的猛进出69影院| 法律面前人人平等表现在哪些方面| 日本欧美视频一区| 成年女人毛片免费观看观看9 | tube8黄色片| 久久草成人影院| 少妇被粗大的猛进出69影院| 免费久久久久久久精品成人欧美视频| 精品国产一区二区三区四区第35| 色综合欧美亚洲国产小说| 国产成人系列免费观看| 久久久久国产一级毛片高清牌| 国产淫语在线视频| 欧洲精品卡2卡3卡4卡5卡区| 国产精品香港三级国产av潘金莲| 成熟少妇高潮喷水视频| 精品国产一区二区三区四区第35| 色综合欧美亚洲国产小说| 一区福利在线观看| 午夜久久久在线观看| tube8黄色片| 欧美黄色淫秽网站| 精品亚洲成国产av| 日本vs欧美在线观看视频| 免费黄频网站在线观看国产| 亚洲人成电影免费在线| 亚洲 国产 在线| 三级毛片av免费| 嫩草影视91久久| 很黄的视频免费| 人人妻,人人澡人人爽秒播| 国产区一区二久久| 亚洲国产精品sss在线观看 | 久久久国产成人免费| 日韩 欧美 亚洲 中文字幕| 窝窝影院91人妻| 久久亚洲真实| 乱人伦中国视频| 国产高清国产精品国产三级| 久久中文字幕人妻熟女| 精品一区二区三区视频在线观看免费 | 亚洲一区二区三区欧美精品| 99国产极品粉嫩在线观看| 国产片内射在线| 人人妻人人澡人人爽人人夜夜| 麻豆成人av在线观看| 很黄的视频免费| 欧美激情 高清一区二区三区| 久久精品亚洲精品国产色婷小说| 午夜精品在线福利| 久久狼人影院| 热re99久久国产66热| 久久久久国内视频| 日韩成人在线观看一区二区三区| 亚洲欧美日韩高清在线视频| cao死你这个sao货| 少妇粗大呻吟视频| 国产精品欧美亚洲77777| 午夜福利免费观看在线| 午夜两性在线视频| 久久久久国产一级毛片高清牌| 日韩三级视频一区二区三区| 水蜜桃什么品种好| 久久久精品国产亚洲av高清涩受| 91麻豆精品激情在线观看国产 | 日韩成人在线观看一区二区三区| 人妻丰满熟妇av一区二区三区 | 精品久久久精品久久久| 香蕉国产在线看| 天天影视国产精品| 久久国产亚洲av麻豆专区| 久久精品成人免费网站| 日韩成人在线观看一区二区三区| 精品国内亚洲2022精品成人 | 亚洲av欧美aⅴ国产| 久久草成人影院| 巨乳人妻的诱惑在线观看| 夫妻午夜视频| 亚洲精品一二三| 看免费av毛片| 欧美日韩乱码在线| 国产乱人伦免费视频| 丝袜人妻中文字幕| 悠悠久久av| 国产激情久久老熟女| 国产男女超爽视频在线观看| 超碰97精品在线观看| 午夜久久久在线观看| 99国产综合亚洲精品| 飞空精品影院首页| 亚洲精品国产色婷婷电影| 天堂中文最新版在线下载| 999久久久国产精品视频| 久久影院123| а√天堂www在线а√下载 | 久久香蕉激情| 一进一出好大好爽视频| 欧美国产精品va在线观看不卡|