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

    A phylogeny of the Tylototriton asperrimus group (Caudata: Salamandridae) based on a mitochondrial study: suggestions for a taxonomic revision

    2011-12-25 06:40:50YUANZhiYongJIANGKeShunQingYANGJunXiaoNGUYENQuangTruongNGUYENThienTaoJINJieQiongCHEJing
    Zoological Research 2011年6期
    關(guān)鍵詞:蠑螈支系智勇

    YUAN Zhi-Yong, JIANG Ke, Lü Shun-Qing, YANG Jun-Xiao, NGUYEN Quang Truong, NGUYEN Thien Tao, JIN Jie-Qiong, CHE Jing,

    (1. State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, the Chinese Academy of Sciences, Kunming Yunnan 650223, China; 2. College of Life and Environment Sciences, Huangshan University, Huangshan 245021, China; 3. Institute of Ecology and Biological Resources, 18 Hoang Quoc Viet St., Hanoi, Vietnam; 4. Vietnam National Museum of Nature, 18 Hoang Quoc Viet St., Hanoi, Vietnam)

    A phylogeny of theTylototriton asperrimusgroup (Caudata: Salamandridae) based on a mitochondrial study: suggestions for a taxonomic revision

    YUAN Zhi-Yong1, JIANG Ke1, Lü Shun-Qing2, YANG Jun-Xiao1, NGUYEN Quang Truong3, NGUYEN Thien Tao4, JIN Jie-Qiong1, CHE Jing1,*

    (1. State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, the Chinese Academy of Sciences, Kunming Yunnan 650223, China; 2. College of Life and Environment Sciences, Huangshan University, Huangshan 245021, China; 3. Institute of Ecology and Biological Resources, 18 Hoang Quoc Viet St., Hanoi, Vietnam; 4. Vietnam National Museum of Nature, 18 Hoang Quoc Viet St., Hanoi, Vietnam)

    A phylogenetic hypothesis for the Asian newts of theTylototriton asperrimusgroup was generated using data from two mitochondrial fragments including COI and the ND1-ND2 regions. Four distinct clades (A, B, C, D) were resolved with high nodal support within this monophyletic group. Clade A includedT. asperrimus,T.hainanensis,T. notialis, “T. vietnamensis”, and two unnamed salamander populations from Vietnam. Clade A, constituted the sister group of clades B + C. Newly identified clade C likely represents a new cryptic species. Clade C was the sister group ofT. wenxianensis. The trueT. vietnamensisexclusively constituted clade D. Our results bring into question some previous taxonomic decisions, and a revision is required. This study illustrates the necessity to include samples from type localities in taxonomic studies, and highlights the importance of fine-grained geographical sampling.

    Tylototriton notialis;Tylototriton hainanensis; Salamander; Southeast Asia; Tonkin; Cryptic diversity

    The Asian newt genusTylototritonAnderson, 1871 occurs from Nepal, through north-central India, Sikkim, Myanmar, northern Thailand, Laos, and Vietnam into southern China. Presently it contains nine species (AmphibiaWeb, 2011). The monophyly of the genus and two groups (or subgenera: Dubois & Raffa?lli, 2009), theT. verrucosusgroup andT. asperrimusgroup, have been consistently supported in recent studies (Stuart et al, 2010; Weisrock et al, 2006).

    With the addition of two newly described species,T. vietnamensisB?hme et al, 2005 andT. notialisStuart et al, 2010, theT. asperrimusgroup includes five species in total, along withT. asperrimusUnterstein, 1930,T. hainanensisFei et al, 1984, andT. wenxianensisFei et al, 1984. However, despite recent progress, uncovering the dimensions of amphibian diversity and inferring evolutionary history remain difficult, in part because of extreme morphological conservation and associated homoplasy (Bossuyt & Milinkovitch, 2000; Wake, 1991). Identifying species is challenging for theT. asperrimusgroup, and determining the exact distributions of some species has been confusing (Fei et al, 2005, 2009; Nguyen et al, 2009), for example,T. asperrimus,T. wenxianensis, andT. vietnamensis. Moreover, previous molecular phylogenetic studies on this group are chiefly based on limited samples and at the same time, unfortunately, materials from the type localities are often absent (Stuart et al, 2010; Weisrock et al, 2006).

    Here we report results of our use of mitochondrial data to investigate the phylogenetic relationships ofTylototriton,with an emphasis on understanding the evolutionary history of theT. asperrimusgroup.

    1 Material and Methods

    1.1 Sampling

    Tissue samples from adults, including toe tips, muscle, and liver, and whole juvenile specimen were collected in the field following Animal Use Protocols approved by the Kunming Institute of Zoology Animal Care and Ethics Committee. The tissues were preserved in 95% ethanol. A total of 36 individuals ofTylototritonfrom 17 localities, representing all 9 currently recognized species in this genus, were examined in this study. Nineteen sequences were downloaded from GenBank for certain comparisons and validations. Outgroup taxa wereEchinotriton andersoniandE. chinhaiensis,chosen mainly based on the study of Weisrock et al (2006). Fig. 1 provides an overview on the distribution of our samples. Details of taxonomic sampling and GenBank accession numbers are provided in Tab. 1. Species identity was based on morphological characters determined in the field. Some specimens could not be reliably identified to species and we refer to them as“sp.”.

    Fig. 1 Map showing the geographic locations of samples representing the Tytolotriton asperrimus group used in present analyses

    Tab. 1 Voucher specimens, localities, and GenBank accession numbers for sample used for the Tylototriton

    1.2 Extraction, amplification and sequencing

    Genomic DNA was extracted using the phenol-chloroform method (Sambrook et al, 1989). Two mitochondrial sequences were amplified. One fragment covers a 1 014 bp region that encodes part of the NADH dehydrogenase subunit 1, tRNA Ile, tRNA Gln, tRNA Met, and part of the NADH dehydrogenase subunit 2 genes (hereafter, ‘ND1 to ND2’). PCR primers KIZ1S (5'-TGACCAATAGCARCAATA-3') and KIZ1A (5'-GTGGGCRATAGATGARTA-3') were designed specifically for successful amplification. Internal primers KIZL4437 (5'-AAGCTTTCGGGCCCATACC-3') and KIZH4419 (5'-AAGCTTTTGGGCCCATACC-3') were applied to sequencing. Another gene is a 561 bp region from cytochrome oxidase subunit I (COI) which was amplified following the standard protocol (Che et al, in press). Amplification was performed in a 25 μL volume reaction with the following procedures: initial denaturation step for 5 min at 94 °C, 35 cycles of denaturation for 1 min at 94 °C, annealing for 1 min at 59 °C for ND1 to ND2 and 45 °C for COI, extension for 1 min at 72 °C. Final extension at 72 °C was conducted for 10 min.

    PCR products were purified with a Gel Extraction Mini Kit (Watson BioTechnologies, Shanghai), and sequenced in both directions using a BigDye Terminator Cycle Sequencing Kit (version 2.0, Applied Biosystems) and a ABI PRISM 3730 automated DNA sequencer (Applied Biosystems). The sequence data were translated into amino acids to confirm the absence of premature stop codons, i.e. to preclude the sequenced of nuclear DNA pseudogenes, and also submitted to a BLAST search in GenBank to confirm that the required sequences had been amplified. Nucleotide sequences were aligned using Clustal X 1.81 (Thompson et al, 1997) with default parameters, and then optimized by eye in MEGA 4.0. (Tamura et al, 2007)

    1.3 Molecular analyses

    Considering that all mtDNA gene sequences are virtually inherited as one linkage group, the two mtDNA gene segments were concatenated into a single partition. The ten samples for which COI sequences are lacking were treated as missing data when combined. Mixedmodel Bayesian phylogenetic inference was performed using Bayesian inference (BI) with MrBayes version 3.1.2 (Ronquist & Huelsenbeck, 2003). Data were partitioned by ND1 codon position, tRNA, ND2 codon position, and COI codon position. The model of sequence evolution that best described each of the data partitions was inferred using the Akaike Information Criterion as implemented in MrModeltest 2.3 (Posada & Crandall, 1998). These were TIM+G for ND1 first codon position; HKY+I for ND1 second codon position; GTR+I for ND1 third codon position and tRNA; TRN+G for ND2 first codon position; HKY+G for ND2 second codon position; GTR+G for ND2 third codon position and COI third codon position; and F81 for COI first and second codon positions. Four independent Bayesian analyses were performed on each dataset. In each analysis, four chains were run for 5 000 000 generations using the default priors, trees were sampled every 100 generations, log-likelihood scores were tracked for stabilization and the first 25% of the trees were discarded as burn-in.

    2 Results and Discussion

    2.1 Molecular analyses

    The aligned sequences yielded 1 014 base pairs of ND1 to ND2 (305 variable, 233 parsimony-informative sites) and 561 bp base pairs of COI (169 variable, 133 parsimony-informative sites). The expected bias against guanine nucleotides was observed (G=12.0%, A= 35.1%, T=25.6 %, and C=27.3 % for ND1 to ND2; and A=25.9%, T=27.4%, C=29.4%, G=17.3% for COI). These data indicated the mitochondrial gene was sequenced, and not a nuclear copy. Forty-seven new sequences were determined and deposited in GenBank (Tab. 1).

    The BI tree is shown in Fig. 2. Consistent with Weisrock et al (2006) and Stuart et al (2010), the monophly of the genus ofTylototritonis strongly supported. A major difference is that our data did not recover the two major clades withinTylototritonreported in previous studies. This is not surprising given our sequencing of a shorter fragment and our far denser sampling forT. asperrimusgroup. The monophyly of theT. asperimmusgroup as shown by Stuart et al (2010) collapsed into two subclades (A+B+C and D) in the present analyses. Clade A including seven strongly supported subclades (A1?A7) was recovered as the sister group of B+C. The position of clade D was not resolved in the present study.

    Fig. 2 A Bayesian tree resulting from mixed-model analysis of mitochondrial DNA

    Our increased sampling reveals previously unrecorded phylogenetic structure. On the one hand,T. asperrimusdoes not contain a single maternal lineage (Fig. 2). The nominal form ofT. asperrimus(type locality: Dayao shan, Guangxi; locality 2) is in lineage A6, which is the sister taxon to the two apparently unnamed samplesT.sp. from Hoa Binh, Vietnam (A5). Lineages A2, A3, and B2 have also been assigned toT. asperrimus. Salamanders from Baise, Guangxi (locality 1) form lineage A2, which is closely related toT. hainanensis(lineage A1). Lineage A3 is from Ha Giang, Vietnam (locality 3), and it is the sister group of “T. vietnamensis” from Cao Bang, Vietnam (locality 6). Lineage B2 is from Libo and Leigongshan, Guizhou, China (localities 4 and 5), is the sister group ofT. wenxianensis(lineage B1). Additionally, the trueT.vietnamensisfrom Bac Giang, Vietnam, exclusively constitutes clade D, which is distantly related to the population from Cao Bang (locality 6).

    On the other hand, compared with previous studies, our data unexpectedly identify a distinct evolutionary clade C from Mangshan, Hunan (locality 9), which is the sister group of clade B.

    2.2 Taxonomic implications

    Species identification within theT. asperrimusgroup is challenging because all members are superficially very similar and conservative in morphology (Stuart et al, 2010). Except forT. hainanensisandT. notialis, which occur in limited geographic regions, Hainan Island and Laos, respectively,T. wenxianensis,T. asperrimus, andT. vietnamensishave more dispersed distributions (Chen, 1991; Dai et al, 2011; Fei et al, 2005, 2009; Nguyen et al, 2009) and species borders are not clear.

    Tylototriton wenxianensiswas reported to occur in Gansu, Sichuan, Chongqing, southeastern Guizhou, Anhui, and nothern and central Hunan provinces in China (Fei et al, 2005, 2009). Our limited samples supported the monophyly ofT. wenxianensisas a separate species, which can be interpreted as support for the taxonomic revision by Fei et al (1984). Moreover, we consider two samples from southeastern and extremely southern Guizhou, originally identified as “T. asperimmus” by Weisrock et al (2006) and our team,to be conspecific withT. wenxianensis,based on the phylogeny and genetic distance estimation. Among these, one sample from Guizhou (locality 4) was unfortunately used to represent the nominalT. asperimmusin comparisons of genetic distances with other species, for example,T. notialis(Stuart et al, 2010). Based on our data, we refer clade B to beT. wenxianensis(Tab. 1, Fig. 3). The level of genetic divergence betweenT. wenxianensis(B) and other clades is relatively higher (0.073?0.103) than other comparisons withinTylototriton(Tab. 2). Denser sampling within this clade is required, including investigation of the recent report of a population ofT. wenxianensisrelatively far to the north in Henan province, China (Chen et al, 2010).

    Fig. 3 Unrooted phylogenetic tree of mitochondrial gene sequences reconstructed using Bayesian methods for the TylototritonPosterior probabilities are displayed only at the major nodes. Taxonomic change suggestions were indicated on this tree.

    Tab. 2 K2P pairwise distances of the ND1-ND2 sequences used in this study for the Tylototriton

    Currently,T. asperrimusis reported to be distributed in Guangxi, Guangdong, northern, western, and southern Guizhou, China (Fei et al, 2005, 2009) and northern Vietnam, for example, Lao Cai and Ha Giang (Nguyen et al, 2009). The distributions ofT. asperrimusandT. vietnamensisin Vietnam require further assessments. Presently,T. vietnamensisis only reported to occur in Cao Bang, Lang Song, and Bac Giang (Nguyen et al, 2009). Based on our new molecular analyses and calculated genetic distances (Tab. 2),T. vietnamensisfrom the type locality of Bac Giang was distinctly separated as our clade D, as hypothesized by Stuart et al (2010). However, the assumed “T. vietnamensis” (orT. cf. vietnamensis: Stuart et al, 2010; Weisrock et al, 2006) from Cao Bang (A4) was distantly separated from the trueT. vietnamensis. Together,T. hainanensis,T. notialis, “T. vietnamensis”, and three unnamed salamanders from Vietnam clustered with the trueT. asperrimusfrom Daoyaoshan, China, which constituted clade A. Based on the evolutionary relationships and genetic distances, we assign all members of our clade A to one species,T. asperrimus(Tab. 1, Fig. 3). Accordingly, we assignT. notialisStuart et al, 2010 andT. hainanensisFei et al, 1984 as junior objective synonyms ofT. asperrimusUnterstein 1930.

    The level of genetic divergence among A1-A7 ranged from 0.021 to 0.047 (not shown here), which is much lower than the hypothesized inter-specific divergence (0.052?0.120, Tab. 2). Samples from clade A are distributed around the Gulf of Tonkin (Fig. 1). The closure of the Tonkin Strait in the past probably allowed gene flow exchanges between populations of clade A.

    Samples from clade C could be treated as a cryptic species, and the level of genetic divergence between lineages C and other clades withinTylototritonranged from 0.073?0.109; however, it also is similar toT. asperrimusin morphology and we could not easily identify it when we are in the field. Further detailed studies including a combination of various datasets should be done to ascertain whether it represents a distinct species (Hou et al, unpublished data), or, like clade B2, it should be assigned toT. wenxianensis. A denser sampling covering the entire distribution ranges regarding clade B and C is desirable.

    Altogether, the description of amphibian diversity and evolutionary history remains a difficult task in some cases because of extremely conservative morphology and in other cases homoplasy (Bossuyt & Milinkovitch, 2000; Wake, 1991). In this case, molecular tools do help, however samples from type localities should be fully considered when formal taxonomic revisions and descriptions of new species are undertaken. Inclusion of additional genes, especially nuclear genes, may prove to be essential. MtDNA data alone cannot assess introgression and gene flow.

    Acknowledgements:We thank HOU M, SHI HT, ORLOV NL, MURPHY RW, YAN F, ZHOU WW, DING L, and LI JT for specimen collection, or providing tissues. We thank WAKE DB and PAPENFUSS TJ for reviewing and offering helpful comments on this manuscript. Technical support from the laboratory of Prof. ZHANG YP is gratefully acknowledged.

    AmphibiaWeb Information on amphibian biology and conservation. [web application]. 2011. Berkeley, California: AmphibiaWeb [DB/OL]. Available: http://amphibiaweb.org/. (Accessed: Nov 16, 2011).

    Anderson J. 1871. Description of a new genus of newts from western Yunan [J].Proc Zool SocLond, 1871: 423-425.

    Bossuyt F, Milinkovitch MC. 2000. Convergent adaptive radiations in Madagascan and Asian ranid frogs reveal covariation between larval and adult traits [J].Proc Natl Acad Sci USA, 97: 6585-6590. B?hme W, Schottler T, Nguyen TQ, K?hler J. 2005. A new species of salamander, genusTylototriton(Urodela: Salamandridae), from northern Vietnam [J].Salamandra, 41: 215-220.

    Che J, Cheng HM, Yang JX, Jing JQ, Jiang K, Yuan ZY, Murphy RW, Zhang YP. Universal COI primers for DNA barcoding amphibians [J].Mol Ecol Resour(in press).

    Chen BH. 1991. The Amphibians and Reptiles Fauna of Anhui [M]. Hefei: Publishing House of Science and Technology. (in Chinese)

    Chen XH, Wang XW, Tao J. 2010. A new subspecies of genusTylototritonfrom China (Caudata: Salamandridae) [J].Acta Zootax Sin, 35: 666-670. (in Chinese).

    Dai ZX, Zhang XM, Zhou Q, Liu ZS. 2011. Amphibian resources in Hubei province [J].Sichuan J Zool, 30: 144-147. (in Chinese).

    Dubois A, Raffa?lli J. 2009. A new ergotaxonomy of the family Salamandridae Goldfuss, 1820 (Amphibia, Urodela) [J].Alytes, 26: 1-85.

    Fei L, Ye CY, Yang R. 1984. A new species and a new subspecies of the genusTylototriton(Caudata: Salamandridae) [J].Acta Zool Sin,30: 85-91. (in Chinese).

    Fei L, Ye CY, Jiang JP, Xie F. 2005. An Illustrated Key to Chinese Amphibians [M]. Chengdu: Sichuan Publishing House of Science and Technology. (in Chinese).

    Fei L, Ye CY, Jiang JP. 2009. Colored Atlas of Chinese Amphibians [M]. Chengdu: Sichuan Publishing House of Science and Technology (in Chinese).

    Nguyen SV, Ho CT, Nguyen TQ. 2009. Herpetofauna of Vietnam [M]. Germany: Edition Chimaira, Frankfurtam Main.

    Posada D, Crandall KA. 1998. MODELTEST: testing the model of DNA substitution [J].Bioinformatics, 14: 817-818.

    Ronquist F, Huelsenbeck JP. 2003. MRBAYES 3: Bayesian phylogenetic inference under mixed models [J].Bioinformatics, 19: 1572-1574.

    Sambrook J, Fritsch E, Maniatis T. 1989. Molecular Cloning: A Laboratory Manual 2nd ed. [M] New York: Cold Spring Harbor Laboratory Pres.

    Stuart BL, Phimmachak S, Sivongxay N, Robichaud WG. 2010. A new species in theTylototriton asperrimusgroup (Caudata: Salamandridae) from central Laos [J].Zootaxa, 2060: 59-68.

    Tamura K, Dudley J, Nei M, Kumar S. 2007. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0 [J].Mol Biol Evol, 24: 1596-1599.

    Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools [J].Nucleic Acids Res, 25: 4876-4882.

    Wake DB. 1991. Homoplasy: the result of natural selection, or evidence of design limitations? [J].Am Nat, 138: 543-567.

    Weisrock DW, Papenfuss TJ, Macey JR, Litvinchuk SN, Polymeni R, Ugurtas IH, Zhao EM, Jowkar H, Larson A. 2006. A molecular assessment of phylogenetic relationships and lineage accumulation rates within the family Salamandridae (Amphibia, Caudata) [J].Mol Phylogenet Evol, 41: 368-383.

    Zhang P, Papenfuss TJ, Wake MH, Qu LH, Wake DB. 2008. Phylogeny and biogeography of the family Salamandridae (Amphibia: Caudata) inferred from complete mitochondrial genomes [J].Mol Phylogenet Evol, 49: 586-597.

    基于線粒體基因?qū)?xì)痣疣螈種組系統(tǒng)發(fā)育關(guān)系的探討及物種修訂建議

    袁智勇1, 蔣 柯1, 呂順清2, 楊軍校1, NGUYEN Quang Truong3, NGUYEN Thien Tao4, 金潔瓊1, 車 靜1,*

    (1.中國(guó)科學(xué)院昆明動(dòng)物研究所 遺傳資源與進(jìn)化國(guó)家重點(diǎn)實(shí)驗(yàn)室,云南 昆明650223; 2.黃山學(xué)院 生命與環(huán)境科學(xué)學(xué)院,安徽 黃山245021; 3. Institute of Ecology and Biological Resources, 18 Hoang Quoc Viet St., Hanoi, Vietnam; 4. Vietnam National Museum of Nature, 18 Hoang Quoc Viet St., Hanoi, Vietnam)

    該文基于兩段線粒體基因COI 和ND1-ND2部分序列片段構(gòu)建了細(xì)痣疣螈種組的系統(tǒng)發(fā)育關(guān)系, 結(jié)果顯示該單系群內(nèi)包含了四個(gè)明顯的進(jìn)化支系(A、B、C、D)。A支系與B、C支系構(gòu)成姐妹群關(guān)系; A支系包含了細(xì)痣疣螈、海南疣螈、老撾疣螈和“越南疣螈”以及另外兩個(gè)來(lái)自越南未被命名的種群; C支系代表了一個(gè)新的獨(dú)立進(jìn)化支, 可能是一個(gè)隱存種;C支系與B支系構(gòu)成姐妹群關(guān)系。此外,來(lái)自模式產(chǎn)地的越南疣螈單獨(dú)構(gòu)成支系D。該文結(jié)果不支持以往細(xì)痣疣螈種組部分物種的劃分, 建議對(duì)該種組進(jìn)行新的分類修訂。該文強(qiáng)調(diào)了在分類研究中應(yīng)該重視使用來(lái)自模式產(chǎn)地的樣品并進(jìn)行廣泛的地理居群采樣。

    老撾疣螈; 海南疣螈; 蠑螈; 東南亞; 東京灣; 隱存多樣性

    Q959.520.9; Q349.1; Q951.3

    A

    0254-5853-(2011)06-0577-08

    2011-11-18;接受日期:2011-11-28

    袁智勇,男,碩士研究生,主要從事東亞和東南亞蠑螈科物種遺傳多樣性及演化研究

    10.3724/SP.J.1141.2011.06577

    date: 2011-11-18; Accepted date: 2011-11-28

    s: This work was supported by grants from the Ministry of Science and Technology of China (MOST Grant 2011FY120200), the National Natural Science Foundation of China (31090100, 30870281), the Chinese Academy of Sciences (KSCX2-YW-Z-0807, KSCX2-EW-Z-2), and the Bureau of Science and Technology of Yunnan province.

    *通訊作者(

    ),女,副研究員。Tel: +86-871-5199030; Fax: +86-871-5133185; E-mail: chej@mail.kiz.ac.cn

    猜你喜歡
    蠑螈支系智勇
    揭秘蠑螈神奇的組織再生能力
    蠑螈的影子
    蠑螈的皮膚受傷后可快速愈合
    High-performance and fabrication friendly polarization demultiplexer
    楊智勇藝術(shù)作品欣賞
    論邊疆哈尼族支系語(yǔ)言
    ——以云南墨江自治縣為例
    廣西凌云地區(qū)背隴瑤的支系研究
    神州民俗(2018年9期)2018-11-21 11:10:02
    身家50億的智勇堅(jiān)守
    苗族支系服飾色彩探析
    流行色(2017年1期)2017-05-31 19:18:01
    中考題中的整式
    中文字幕最新亚洲高清| 视频区图区小说| 国产精品一国产av| 亚洲国产av新网站| 久久综合国产亚洲精品| 女性被躁到高潮视频| 在线 av 中文字幕| 久久久久久久国产电影| 久久精品国产鲁丝片午夜精品| 色婷婷久久久亚洲欧美| 丝瓜视频免费看黄片| 精品亚洲成a人片在线观看| 青春草视频在线免费观看| 亚洲国产毛片av蜜桃av| 黄频高清免费视频| 99香蕉大伊视频| 最近最新中文字幕免费大全7| 中文精品一卡2卡3卡4更新| 菩萨蛮人人尽说江南好唐韦庄| 人人妻人人添人人爽欧美一区卜| 免费日韩欧美在线观看| 狠狠精品人妻久久久久久综合| 亚洲国产毛片av蜜桃av| 天天躁狠狠躁夜夜躁狠狠躁| 伦精品一区二区三区| 丝袜美腿诱惑在线| 亚洲色图 男人天堂 中文字幕| 成人国语在线视频| 久久综合国产亚洲精品| 在现免费观看毛片| 日日啪夜夜爽| 亚洲人成网站在线观看播放| 欧美日韩一区二区视频在线观看视频在线| 一区二区三区精品91| 搡女人真爽免费视频火全软件| 亚洲综合色网址| 在线观看美女被高潮喷水网站| 精品少妇内射三级| 中文字幕最新亚洲高清| 嫩草影院入口| 国产亚洲av片在线观看秒播厂| 80岁老熟妇乱子伦牲交| 国产男女超爽视频在线观看| 欧美日韩视频精品一区| 国产精品亚洲av一区麻豆 | 在线观看免费日韩欧美大片| 久久这里只有精品19| 熟女少妇亚洲综合色aaa.| 国产精品女同一区二区软件| 欧美 亚洲 国产 日韩一| 国产乱来视频区| 日日啪夜夜爽| 精品人妻熟女毛片av久久网站| 叶爱在线成人免费视频播放| 国产一区二区三区av在线| 女人精品久久久久毛片| 日韩三级伦理在线观看| 人人妻人人添人人爽欧美一区卜| 久久精品夜色国产| 亚洲美女搞黄在线观看| 精品国产一区二区三区四区第35| 免费大片黄手机在线观看| 中文字幕人妻丝袜制服| 久久99精品国语久久久| 精品久久久久久电影网| 如日韩欧美国产精品一区二区三区| 色吧在线观看| 国产免费一区二区三区四区乱码| 蜜桃国产av成人99| 午夜福利一区二区在线看| 国产在线视频一区二区| 视频在线观看一区二区三区| 99热网站在线观看| 日韩成人av中文字幕在线观看| 美国免费a级毛片| 一级黄片播放器| 女人高潮潮喷娇喘18禁视频| 九草在线视频观看| 成人毛片60女人毛片免费| 如日韩欧美国产精品一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 成人18禁高潮啪啪吃奶动态图| 一二三四中文在线观看免费高清| 激情视频va一区二区三区| 国产亚洲av片在线观看秒播厂| 久久精品国产亚洲av涩爱| 大话2 男鬼变身卡| 国产一区二区在线观看av| 一级片'在线观看视频| 久久97久久精品| 国产 精品1| 韩国精品一区二区三区| 宅男免费午夜| 精品人妻偷拍中文字幕| 成人国产av品久久久| 99久久人妻综合| 少妇的逼水好多| 黑丝袜美女国产一区| 国产精品成人在线| 国产男女内射视频| 不卡av一区二区三区| 如何舔出高潮| 久久婷婷青草| 黄色一级大片看看| 免费高清在线观看视频在线观看| av网站在线播放免费| 午夜久久久在线观看| 一二三四中文在线观看免费高清| 在线观看www视频免费| 欧美精品高潮呻吟av久久| 亚洲在久久综合| 亚洲一码二码三码区别大吗| 亚洲精品自拍成人| 国产一区亚洲一区在线观看| 卡戴珊不雅视频在线播放| 国产精品99久久99久久久不卡 | 美国免费a级毛片| 亚洲精品国产av成人精品| 少妇猛男粗大的猛烈进出视频| 欧美精品人与动牲交sv欧美| 国产一区二区激情短视频 | 搡女人真爽免费视频火全软件| 中文字幕人妻熟女乱码| 黑人猛操日本美女一级片| 国产在视频线精品| 精品少妇一区二区三区视频日本电影 | 欧美人与性动交α欧美精品济南到 | 亚洲国产av新网站| 一区二区三区四区激情视频| 国产成人aa在线观看| 色播在线永久视频| 99久久人妻综合| 97在线人人人人妻| 久久久久久久久久久免费av| 最新中文字幕久久久久| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品电影小说| 日日啪夜夜爽| 国产男女超爽视频在线观看| 人人妻人人澡人人看| 亚洲欧洲精品一区二区精品久久久 | 欧美人与善性xxx| 巨乳人妻的诱惑在线观看| 99精国产麻豆久久婷婷| 免费在线观看黄色视频的| 精品午夜福利在线看| 嫩草影院入口| 中文字幕亚洲精品专区| 宅男免费午夜| 国产精品 国内视频| 国产成人免费观看mmmm| 综合色丁香网| 老汉色∧v一级毛片| 女人高潮潮喷娇喘18禁视频| videosex国产| 日韩不卡一区二区三区视频在线| 综合色丁香网| 国产精品免费视频内射| 久久久欧美国产精品| 国产免费福利视频在线观看| 亚洲一级一片aⅴ在线观看| 亚洲,一卡二卡三卡| 一边摸一边做爽爽视频免费| 国产av一区二区精品久久| 熟女av电影| 永久免费av网站大全| 成人国产av品久久久| 婷婷色麻豆天堂久久| 婷婷色麻豆天堂久久| 男女边吃奶边做爰视频| 最近2019中文字幕mv第一页| 纯流量卡能插随身wifi吗| 18禁国产床啪视频网站| 超色免费av| 黑人猛操日本美女一级片| 色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| 狠狠婷婷综合久久久久久88av| 亚洲欧洲国产日韩| 亚洲国产最新在线播放| 亚洲精品国产av蜜桃| 欧美激情 高清一区二区三区| 一级片免费观看大全| 十分钟在线观看高清视频www| 啦啦啦视频在线资源免费观看| 午夜日韩欧美国产| 91成人精品电影| 亚洲中文av在线| 久久精品国产自在天天线| 青春草亚洲视频在线观看| 老司机亚洲免费影院| 国产日韩欧美在线精品| 边亲边吃奶的免费视频| 免费av中文字幕在线| 9热在线视频观看99| 欧美在线黄色| 久久国产精品大桥未久av| 久久久精品国产亚洲av高清涩受| av天堂久久9| 久久精品久久精品一区二区三区| 777久久人妻少妇嫩草av网站| 少妇人妻精品综合一区二区| 国产淫语在线视频| 精品人妻在线不人妻| 欧美日韩一区二区视频在线观看视频在线| 黄片无遮挡物在线观看| 两个人看的免费小视频| 久久久久久人妻| 午夜影院在线不卡| 日韩视频在线欧美| 人体艺术视频欧美日本| 亚洲美女黄色视频免费看| 2018国产大陆天天弄谢| 国产在视频线精品| 国产一级毛片在线| 国产高清不卡午夜福利| 日韩三级伦理在线观看| 亚洲欧美精品综合一区二区三区 | 亚洲经典国产精华液单| 精品久久久久久电影网| 国产精品嫩草影院av在线观看| 国产精品二区激情视频| 中文字幕色久视频| 午夜91福利影院| 欧美成人午夜免费资源| 成人亚洲欧美一区二区av| 日韩不卡一区二区三区视频在线| 午夜免费鲁丝| 少妇人妻精品综合一区二区| 热99久久久久精品小说推荐| 黄片无遮挡物在线观看| 丰满迷人的少妇在线观看| 在线免费观看不下载黄p国产| 国产精品国产av在线观看| 视频在线观看一区二区三区| 黑人欧美特级aaaaaa片| 亚洲精品国产av蜜桃| 成人免费观看视频高清| 亚洲av男天堂| 又黄又粗又硬又大视频| 日本欧美视频一区| 久久久亚洲精品成人影院| 欧美中文综合在线视频| 日韩成人av中文字幕在线观看| 亚洲精品国产av成人精品| 在线免费观看不下载黄p国产| 高清在线视频一区二区三区| 久久这里有精品视频免费| 十八禁网站网址无遮挡| 欧美激情高清一区二区三区 | 久久人人97超碰香蕉20202| 丝瓜视频免费看黄片| 欧美老熟妇乱子伦牲交| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 十八禁高潮呻吟视频| 亚洲第一青青草原| 久久久久精品久久久久真实原创| av卡一久久| 国产亚洲精品第一综合不卡| 熟女电影av网| 亚洲综合色惰| 三级国产精品片| 午夜免费男女啪啪视频观看| 日韩在线高清观看一区二区三区| 亚洲精品视频女| 欧美日韩视频高清一区二区三区二| 亚洲精品国产av蜜桃| 亚洲人成77777在线视频| 中国三级夫妇交换| 男人爽女人下面视频在线观看| 美女中出高潮动态图| 欧美精品一区二区大全| 久久人人爽av亚洲精品天堂| 国产日韩欧美亚洲二区| 满18在线观看网站| 国产日韩一区二区三区精品不卡| 97在线人人人人妻| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 18禁动态无遮挡网站| 国产av码专区亚洲av| 高清av免费在线| 一区二区三区激情视频| av有码第一页| 成人国语在线视频| 国产 精品1| 欧美黄色片欧美黄色片| 一区二区av电影网| 人妻人人澡人人爽人人| 人妻 亚洲 视频| 高清黄色对白视频在线免费看| 激情五月婷婷亚洲| 在线免费观看不下载黄p国产| 亚洲成av片中文字幕在线观看 | 一区二区av电影网| a 毛片基地| 成人午夜精彩视频在线观看| 欧美精品一区二区免费开放| 国产成人免费无遮挡视频| 午夜福利视频在线观看免费| 一级爰片在线观看| 一级毛片黄色毛片免费观看视频| 99热网站在线观看| 亚洲欧美日韩另类电影网站| 在线观看一区二区三区激情| 汤姆久久久久久久影院中文字幕| av在线播放精品| 狠狠精品人妻久久久久久综合| 国产成人精品久久久久久| 国产一区亚洲一区在线观看| 久久精品夜色国产| 老汉色av国产亚洲站长工具| 亚洲精品自拍成人| 久久影院123| 黄色 视频免费看| 午夜福利,免费看| 日韩伦理黄色片| 亚洲成人手机| 国产日韩欧美视频二区| 两个人免费观看高清视频| 欧美亚洲 丝袜 人妻 在线| 日本爱情动作片www.在线观看| 精品亚洲乱码少妇综合久久| 午夜日本视频在线| 日日摸夜夜添夜夜爱| 巨乳人妻的诱惑在线观看| 日韩av在线免费看完整版不卡| 啦啦啦视频在线资源免费观看| 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| xxxhd国产人妻xxx| 熟妇人妻不卡中文字幕| 成年美女黄网站色视频大全免费| 亚洲国产精品一区三区| 日韩制服丝袜自拍偷拍| 国产女主播在线喷水免费视频网站| 黑人欧美特级aaaaaa片| av在线app专区| 欧美激情高清一区二区三区 | 国产福利在线免费观看视频| 91久久精品国产一区二区三区| 999精品在线视频| 国产精品一区二区在线观看99| 国产亚洲午夜精品一区二区久久| 18禁动态无遮挡网站| 18禁国产床啪视频网站| 亚洲av男天堂| 婷婷成人精品国产| 满18在线观看网站| 国产一区二区 视频在线| 自线自在国产av| 午夜福利在线观看免费完整高清在| 久久久国产一区二区| 建设人人有责人人尽责人人享有的| 国产一区二区激情短视频 | 久久精品aⅴ一区二区三区四区 | 日日爽夜夜爽网站| 巨乳人妻的诱惑在线观看| 国产熟女欧美一区二区| 丰满迷人的少妇在线观看| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| 超碰成人久久| 国产白丝娇喘喷水9色精品| 亚洲内射少妇av| 日本欧美国产在线视频| 亚洲,欧美,日韩| 精品国产乱码久久久久久小说| 久久久久久久大尺度免费视频| 香蕉国产在线看| 久久精品久久久久久久性| 亚洲激情五月婷婷啪啪| 一边摸一边做爽爽视频免费| 制服诱惑二区| 女人被躁到高潮嗷嗷叫费观| 99热国产这里只有精品6| 欧美成人精品欧美一级黄| 午夜av观看不卡| 美女视频免费永久观看网站| 永久免费av网站大全| 欧美日韩av久久| 色94色欧美一区二区| 国产亚洲欧美精品永久| 97精品久久久久久久久久精品| 亚洲中文av在线| 亚洲综合色惰| 热99国产精品久久久久久7| 赤兔流量卡办理| 亚洲国产av新网站| 男女午夜视频在线观看| 欧美日韩视频精品一区| 国产精品嫩草影院av在线观看| 亚洲欧美精品自产自拍| 国产免费视频播放在线视频| 宅男免费午夜| 亚洲av福利一区| 婷婷色综合www| h视频一区二区三区| 国产国语露脸激情在线看| 久久久久精品性色| 久久国产精品大桥未久av| 青青草视频在线视频观看| 乱人伦中国视频| 九草在线视频观看| 午夜激情久久久久久久| 韩国精品一区二区三区| 日本av免费视频播放| 在线 av 中文字幕| 最近中文字幕2019免费版| 久久久久国产精品人妻一区二区| 女性生殖器流出的白浆| 一级毛片我不卡| 久久这里有精品视频免费| 国产精品秋霞免费鲁丝片| 一边摸一边做爽爽视频免费| 日本-黄色视频高清免费观看| 国产一区二区 视频在线| 精品一区二区三卡| 亚洲欧美成人综合另类久久久| 建设人人有责人人尽责人人享有的| 在线看a的网站| 在线观看美女被高潮喷水网站| 人妻一区二区av| 日韩av免费高清视频| 亚洲精品国产av蜜桃| 国产精品国产三级专区第一集| 不卡av一区二区三区| 在线天堂最新版资源| 国产av国产精品国产| 久热这里只有精品99| 热re99久久国产66热| 日韩制服丝袜自拍偷拍| 黄网站色视频无遮挡免费观看| 另类亚洲欧美激情| 少妇被粗大猛烈的视频| 超碰成人久久| 欧美黄色片欧美黄色片| 亚洲欧洲国产日韩| 亚洲av免费高清在线观看| 久久av网站| 一级毛片 在线播放| 国产1区2区3区精品| 精品亚洲成a人片在线观看| 国产有黄有色有爽视频| 中文字幕人妻熟女乱码| 丝袜脚勾引网站| 国产精品二区激情视频| 看非洲黑人一级黄片| 中文字幕精品免费在线观看视频| 尾随美女入室| 欧美精品人与动牲交sv欧美| 久久 成人 亚洲| 久久久国产欧美日韩av| 中文天堂在线官网| 黄色毛片三级朝国网站| 精品人妻一区二区三区麻豆| 欧美精品人与动牲交sv欧美| 桃花免费在线播放| 老司机影院成人| 欧美精品高潮呻吟av久久| 国产高清不卡午夜福利| 亚洲精品国产色婷婷电影| 日韩中文字幕欧美一区二区 | av卡一久久| 丝袜喷水一区| 久久99一区二区三区| 亚洲综合色网址| 建设人人有责人人尽责人人享有的| 伦精品一区二区三区| 日韩中字成人| 人妻人人澡人人爽人人| 成年女人在线观看亚洲视频| 国产成人精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 久久久久久人妻| 伊人亚洲综合成人网| 黄色怎么调成土黄色| 丝袜喷水一区| 中文字幕人妻熟女乱码| 丝袜在线中文字幕| 在线观看免费日韩欧美大片| av不卡在线播放| 亚洲第一区二区三区不卡| 亚洲婷婷狠狠爱综合网| 久久女婷五月综合色啪小说| 日日摸夜夜添夜夜爱| 亚洲成人av在线免费| 97人妻天天添夜夜摸| 五月伊人婷婷丁香| 欧美av亚洲av综合av国产av | 国产成人精品一,二区| 亚洲久久久国产精品| 激情五月婷婷亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 国产熟女午夜一区二区三区| 激情视频va一区二区三区| 亚洲av国产av综合av卡| 亚洲情色 制服丝袜| 夫妻性生交免费视频一级片| 曰老女人黄片| 美女主播在线视频| 欧美精品高潮呻吟av久久| 久久午夜综合久久蜜桃| 国产高清不卡午夜福利| 久久久久国产精品人妻一区二区| 又大又黄又爽视频免费| 欧美97在线视频| 成人漫画全彩无遮挡| 国产精品.久久久| 亚洲av免费高清在线观看| 国产精品偷伦视频观看了| 精品亚洲乱码少妇综合久久| 亚洲精品国产色婷婷电影| 久久久国产精品麻豆| 五月天丁香电影| 91aial.com中文字幕在线观看| 亚洲精品一区蜜桃| 亚洲精品国产色婷婷电影| 国产成人精品福利久久| 精品久久久精品久久久| 国产成人精品久久二区二区91 | 精品国产超薄肉色丝袜足j| 久久久久久伊人网av| kizo精华| 亚洲欧美清纯卡通| 桃花免费在线播放| 一级爰片在线观看| 久久久久久久精品精品| 亚洲av日韩在线播放| 丝袜在线中文字幕| 精品午夜福利在线看| 精品人妻熟女毛片av久久网站| 热99国产精品久久久久久7| 如日韩欧美国产精品一区二区三区| 一个人免费看片子| 日本av手机在线免费观看| 久久久欧美国产精品| 高清不卡的av网站| 国产精品 国内视频| 国产老妇伦熟女老妇高清| 国产精品一国产av| 男人添女人高潮全过程视频| 大片电影免费在线观看免费| 曰老女人黄片| 亚洲激情五月婷婷啪啪| 亚洲图色成人| 日韩大片免费观看网站| a级片在线免费高清观看视频| 最新中文字幕久久久久| 欧美变态另类bdsm刘玥| 久久精品国产综合久久久| 超色免费av| av在线老鸭窝| 免费观看性生交大片5| 欧美精品一区二区大全| 国产精品久久久久久久久免| 国产国语露脸激情在线看| 如何舔出高潮| 免费av中文字幕在线| 宅男免费午夜| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 国产精品久久久久久精品电影小说| 久久精品久久久久久噜噜老黄| 乱人伦中国视频| 亚洲欧洲国产日韩| 午夜影院在线不卡| 中文字幕色久视频| 国产欧美亚洲国产| 午夜福利视频精品| 熟女电影av网| 如何舔出高潮| 免费不卡的大黄色大毛片视频在线观看| 天天影视国产精品| 久久热在线av| a级片在线免费高清观看视频| 国语对白做爰xxxⅹ性视频网站| 国产女主播在线喷水免费视频网站| 欧美日韩av久久| 日韩一卡2卡3卡4卡2021年| 久久久久久久亚洲中文字幕| 国产综合精华液| 波野结衣二区三区在线| 超碰成人久久| 欧美日韩综合久久久久久| 十分钟在线观看高清视频www| 一区福利在线观看| 丰满迷人的少妇在线观看| 免费观看在线日韩| 水蜜桃什么品种好| 两性夫妻黄色片| av在线老鸭窝| 国产精品亚洲av一区麻豆 | 伊人久久大香线蕉亚洲五| 午夜久久久在线观看| 黄片小视频在线播放| 欧美精品av麻豆av| 日本免费在线观看一区| 黄色视频在线播放观看不卡| 国产成人免费观看mmmm| 日韩精品免费视频一区二区三区| 热re99久久国产66热| 国产福利在线免费观看视频| 免费黄色在线免费观看| 少妇人妻久久综合中文| 国产精品久久久av美女十八| 一区二区三区四区激情视频| 叶爱在线成人免费视频播放| 亚洲国产色片| 美女主播在线视频| 亚洲av在线观看美女高潮| 久久影院123| 性高湖久久久久久久久免费观看| 亚洲精品成人av观看孕妇|