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

    Complete mitochondrial genome of the leaf muntjac(Muntiacus putaoensis) and phylogenetics of the genus Muntiacus

    2017-12-14 03:35:02GuoGangLiMingXiaZhangKyawSwaKyawWinMaungRuiChangQuanSoutheastAsiaBiodiversityResearchInstituteChineseAcademyofSciencesYezinNayPyiTaw058MyanmarCenterforIntegrativeConservationXishuangbannaTropicalBotanicalGardenChinese
    Zoological Research 2017年5期

    Guo-Gang Li, Ming-Xia Zhang, Kyaw Swa, Kyaw-Win Maung, Rui-Chang Quan,*1 Southeast Asia Biodiversity Research Institute, Chinese Academy of Sciences, Yezin Nay Pyi Taw 058, Myanmar Center for Integrative Conservation, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla Yunnan 6660,China Hponkan Razi Wildlife Sanctuary Offices, Putao Kachin 01051, Myanmar Forest Research Institute, Forest Department Ministry of Environmental Conservation and Forestry, Yezin Nay Pyi Taw 058, Myanmar

    Complete mitochondrial genome of the leaf muntjac(Muntiacus putaoensis) and phylogenetics of the genus Muntiacus

    Guo-Gang Li1,2, Ming-Xia Zhang1,2, Kyaw Swa3, Kyaw-Win Maung4, Rui-Chang Quan1,2,*1Southeast Asia Biodiversity Research Institute, Chinese Academy of Sciences, Yezin Nay Pyi Taw 05282, Myanmar2Center for Integrative Conservation, Xishuangbanna Tropical Botanical Garden, Chinese Academy of Sciences, Mengla Yunnan 666303,China3Hponkan Razi Wildlife Sanctuary Offices, Putao Kachin 01051, Myanmar4Forest Research Institute, Forest Department Ministry of Environmental Conservation and Forestry, Yezin Nay Pyi Taw 05282, Myanmar

    The leaf muntjac (Muntiacus putaoensis) is an endemic deer species found in the east trans-Himalayan region. In recent years, population numbers have decreased due to heavy hunting and habitat loss, and little genetic data exists for this species, thus our knowledge of distribution rangs and population sizes likewise remain limited. We obtained mtDNA genes and the complete mitochondrial genome sequence of M. putaoensis using PCR, followed by direct sequencing. The complete mitogenome sequence was determined as a circular 16 349 bp mitochondrial genome,containing 13 protein-coding genes, two rRNA genes,22 tRNA genes, and one control region, the gene composition and order of which were similar to most other vertebrates so far reported. Most mitochondrial genes, except for ND6 and eight tRNAs, were encoded on the heavy strand. The overall base composition of the heavy strand was 33.1% A,29.3% T, 24.2% C, and 13.4% G, with a strong AT bias of 62.4%. There were seven regions of gene overlap totaling 95 bp and 11 intergenic spacer regions totaling 74 bp. Phylogenetic analyses (ML and BI) among the Muntiacus genus based on the sequenced of mitogenome and ND4L-ND4 supported M. putaoensis as a member of Muntiacus,most closely related to M. vuquangensis. However,when analyses based on cyt b included two more muntjacs, M. truongsonensis was most closely related to M. putaoensis rather than M.vuquangensis, and together with M. rooseveltorum,likely forming a M. rooseveltorum complex of the species. This study will help in the exploration of the evolutionary history and taxonomic status of the leafmuntjac, as well as its protection as a genetic resource.

    Muntiacus; Muntiacus putaoensis;Mitogenome; Phylogenetics

    INTRODUCTION

    Muntjacs (Muntiacus spp., Muntiacinae, Cervidae) are a group of small solitary deer occurring in the forests of Asia. They are of great interest to evolutionary biologists and mammalogists due to their considerable chromosome variations (from 2n=6(Muntiacus reevesi) to 2n=46 (M. muntjak)) (Fontana & Rubini,1990) and discovery of new species in the last decades of the twentieth century, including M. vuquangensis (Tuoc et al., 1994),M. truongsonensis (Giao et al., 1998), and M. putaoensis(Amato et al., 1999; Rabinowitz et al., 1999).1

    Muntiacus putaoensis is the most recently discovered species in this muntjac genus. It is named after the town of Putao, located in the most northern part of Myanmar, which is also its most recognizable reference point (Amato et al., 1999).Compared with other species, it is the smallest muntjac (mean adult body mass 12 kg), and is only half the size of the sympatric M. muntjak (22–29 kg) (Rabinowitz et al., 1999).Local people call it the ‘leaf deer’ because it is so small it can be wrapped in a single leaf of Phrynium capitatum (Rabinowitz &Khaing, 1998; Rabinowitz et al., 1999). Until recently,information on the distribution range of the species was limitedto northern Myanmar. However, it has since been found in the adjoining hill forests of southeastern Tibet, with the local name similar in meaning to its Myanmar name (Choudhury, 2007;Choudhury, 2009; Datta et al., 2003; James et al., 2008).Muntiacus putaoensis is endemic to the east trans-Himalayan region (including northern Myanmar and southeastern Tibet,China), the population of the leaf muntjac has decreased dramatically in recent years due to considerable degradation of its natural habitat and heavy hunting pressure (Rabinowitz &Khaing, 1998; Rao et al., 2010; Schaller & Rabinowitz, 2004).This species is listed as data deficient (DD) in the IUCN Red List of Threatened species as there is a lack of certainty about its taxonomy, distribution, population, natural history, and threats (Timmins & Duckworth, 2016). Previously the species was only known to occur along a 70 km east to west stretch in northern Myanmar, the new discovery means a two-fold increase in the total east-west range of this species (James et al., 2008). It also indicates that the extent of its geographic range is still being determined.

    The leaf muntjac was characterized and confirmed primarily by its diagnostic mitochondrial DNA (mtDNA) with several partial fragments (Amato et al., 1999). However, its taxonomic status remains controversial as it falls within a group of closely related ‘little’ muntjacs with strong morphological similarities,including M. rooseveltorum and M. truongsonensis (James, et al., 2008; Timmins & Duckworth, 2016). This species is characterized by short thin pedicles, small unbranched antlers,and relatively large preorbital fossa, with both males and females possessing canines (Rabinowitz et al., 1999). Because hunters usually cut off the lower jaw and part of the upper jaw, it is difficult to obtain complete skulls to compare with other muntjacs. In addition, the similarities in morphology also create difficulties when marking comparisons, thus requiring molecular genetic analyses. To date however, genetic data remain scarce,with only partial segments of mtDNA currently sequenced(Amato et al., 1999; James et al., 2008), which are insufficient for confirming the identity of this species or clarifying its phylogenetic relationship relative to other muntjacs.

    Here, we sequenced and analyzed the complete mitochondrial genome and two mtDNA segments of M. putaoensis, the aim of this study was to: 1) provide fundamental genetic data for further conservation genetic studies for this near cryptic mammal; 2) discuss the possible relationships among the genus of Muntiacus, especially that of M. putaoensis, M.rooseveltorum, and M. truongsonensis; and, 3) provide further discussion on the evolutionary history of this species.

    MATERIALS AND METHODS

    Sampling and laboratory procedures

    Seventeen M. putaoensis specimens were collected from Putao,Kachin state, northern Myanmar in three field trips (December 2015, December 2016 to January 2017, and May 2017) (Figure 1; Supplementary Table S1). Three skin samples were preserved in 95% ethanol, and 14 dry meat samples were preserved using silica gel for subsequent analyses. Voucher specimens were deposited in the Southeast Asia Biodiversity Research Institute, Chinese Academy of Sciences, in Nay Pyi Taw, Myanmar (Supplementary Table S1).

    Total genomic DNA was extracted from tissue using a DNeasy Blood & Tissue Kit (Qiagen, Shanghai, China). The master mixture for polymerase chain reaction (PCR) contained approximately 100 ng of template DNA, 1 μL (10 pmol) of each primer, 5 μL of 10× reaction buffer, 2 μL of dNTPs (2.5 mmol/L of each), and 2.0 U of Taq DNA polymerase, in a total volume of 50 μL. Reactions were carried out on a Veriti Thermal Cycler(Applied Biosystems, Carlsbad, CA, USA) and always included a negative control. The M. putaoensis mtDNA segments (ND4LND4 and cyt b) were amplified using PCR with universal primers from previous studies (Irwin et al., 1991; Wang & Lan,2000). After visualization of the fragments using 1% agarose gel,the PCR products were sequenced from both ends using an ABI PRISM 3700 sequencing system with using the same primers as for PCR (Beijing Tianyi Huiyuan Bioscience and Technology Incorporation, Beijing, China).

    To obtain the whole mitogenomic sequence of M. putaoensis,we designed 16 pairs of primers based on previous studies(Hassanin et al., 2012; Larkin et al., 2007; Martins et al., 2017;Shi et al., 2003; Wu & Fang, 2005; Zhang et al., 2004a, b). The primer sequences and approximate lengths of the amplified fragments from the whole genome are listed in Supplementary Table S2. All the primers were synthesized by Beijing Tianyi Huiyuan Bioscience and Technology Incorporation (Beijing, China).

    Sequence analyses

    Two mitochondrial DNA segments sequences were edited using the program DNASTAR 5.0 (DNASTAR Inc.), and were aligned using the CLUSTALW algorithm implemented in MEGA 5.05 with default parameters (Larkin et al., 2007; Tamura et al., 2011).Positions of ambiguous alignment and those from variable length intergenic regions were excluded from further analyses.Identical haplotypes were collapsed using DNASP 5.1 (Librado &Rozas, 2009). The base composition of the mitogenomic sequence was analyzed using MEGA 5.05 (Tamura et al., 2011).We annotated the genome sequence using DOGMA (Wyman et al., 2004).

    Phylogenetic analyses

    Phylogenies of the two mtDNA segments and mitogenome were constructed using maximum likelihood (ML) implemented in PHYML 3.0 (Guindon et al., 2010). Bayesian inference (BI)was implemented in MRBAYES 3.2.1 (Ronquist et al., 2012)using different parameter estimates for the two mtDNA segments and mitogenome. The most appropriate nucleotide substitution models for the two segments and genome were selected using the Akaike Information Criterion in JMODELTEST 2.1.4 (Darriba et al., 2012). To assess the statistical significance of the hypothesized lineages, bootstrap analysis with 200 replicates was used for the ML analyses, with other settings set to default. The posterior distributions were obtained by Markov Chain Monte Carlo (MCMC) analysis with one cold chain and three heated chains. Samples of the trees and parameters were drawn every 100 steps from a total of one million MCMC generations. Three additional runs were conducted, beginning with random trees. A 50% majority rule consensus of the postburn (using a burn-in of 25%) for all generations was computed for all four runs. Elaphodus cephalophus and Tragulus kanchil were chosen as outgroups to root the tree.

    Figure 1 Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees (GTR+G model) for Muntiacus based on 20 complete genomes from mitochondrial DNA

    RESULTS AND DISCUSSION

    Sequence data and genome content

    Cyt b sequences with a total length 1 140 bp were obtained from all 17 M. putaoensis individuals with four variable sites.We obtained 1 844 bp of unambiguous sequences (ND4LND4) from 13 individuals, which contained the complete sequences of NADH dehydrogenase subunit 4L (ND4L),NADH dehydrogenase subunit 4 (ND4), tRNASer, and tRNAHisgenes, and a partial sequence of the tRNALeugene, with 15 variable sites and six parsimony-informative sites. We detected three cyt b haplotypes and eight ND4L-ND4 haplotypes separately (GenBank accession Nos. MF737179–MF737190, Supplementary Table S1).

    The complete mitochondrial genome features of M. putaoensis were identical to those of most other muntjac deer (e.g.Hassanin et al., 2012; Shi et al., 2003; Wu & Fang, 2005;Zhang et al., 2004a, b). The complete mitogenome sequence was a circular molecule 16 349 bp in length, and included 13 protein-coding genes (ATP6, ATP8, COI-III, cyt b, ND1-6, and ND4L), two rRNA genes (12S and 16S rRNA), 22 tRNA genes,and one non-coding control region (D-loop), most of which were encoded on the heavy strand, except for ND6 and eight tRNAs(tRNAGln, tRNAAla, tRNAAsn, tRNACys, tRNATyr, tRNASer, tRNAGlu,and tRNAPro). Notably, gene overlap and separation, common characteristics found in other vertebrate mitochondrial genomes,were also observed in the M. putaoensis mitogenome. There were seven regions of gene overlap totaling 95 bp (varying from 1 to 40 bp) and 11 intergenic spacer regions totaling 74 bp(varying from 1 to 32 bp) (Table 1). The overall base composition of the heavy strand was 33.1% A, 29.3% T, 24.2%C, and 13.4% G with a strong AT bias of 62.4%. Nearly all 13 protein-coding genes started with the common initiation codon ATG, whereas ND2, ND3, and ND5 started with ATA. In addition,10 protein-coding genes shared the complete stop codons TAA or TAG (TAA for ND1, COI, COII, ATP6, ND4L, ND5, and ND6,and TAG for ND2, ATP8, and ND3), whereas the remaining three protein-coding genes (ND4, COIII and cyt b) possessed incomplete stop codons with a terminal T or TA. The 22 tRNA genes were interspersed along the genome, with lengths varying from 60 to 75 bp. The lengths of 12S rRNA and 16S rRNA were 958 and 1 568 bp, respectively, located between tRNAPheand tRNALeu(UUR)and separated by tRNAVal. The D-loop region was located between tRNAProand tRNAPhe, and was 914 bp in length. The base composition of the D-loop was A-29.9%,T-31.9%, C-23.1%, and G-15.1%, re fl ecting a strong feature in A and T (A+T=61.8%).

    Table 1 Characteristics of the mitochondrial genome of Muntiacus putaoensis

    Continued

    Phylogenetic analysis

    The phylogeny of the genus Muntiacus has long been debated(James et al., 2008; Timmins & Duckworth, 2016). Here, we conducted a phylogenetic analysis using complete mitochondrial genomes from five muntjac species (others have not yet been sequenced for mitochondrial genomes) (Figure 1). Additional sequences and species were included in phylogenetic analysis based on ND4L-ND4 (Figure 2A) and cyt b (Figure 2B). The ML and BI analyses recovered the same tree topologies (Figure 1 and 2).

    Figure 2 Maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees (GTR + G model) for Muntiacus based on ND4L-ND4 genes (A) and cyt b gene (B)

    The three different level sequences obtained in this study (i.e.,complete mitogenome, ND4L-ND4, and cyt b) provide insight into the taxonomic status of Muntiacus.

    Based on the complete genomes (Figure 1), the data showed two major clades within the genus of Muntiacus, one of which included M. putaoensis, M. vuquangensis and M. reevesi, and the other M. crinifrons and M. muntjak. In the latter clade, M.crinifrons was positioned on a distinct branch as sister to M.muntjak. Within M. muntjak, three phylogenetically and geographically distinct groups were previously identified by Martins et al. (2017). We determined the leaf muntjac (M.putaoensis) to be a member of genus Muntiacus, most closely related to M. vuquangensis, with M. reevesi as their sister species. However, as genomic data for other muntjacs is still unavailable, further genomic research on more species may illuminate different relationships within this taxon.

    Previous phylogenetic information has been provided for the genus Muntiacus from mtDNA genes (Amato et al., 1999;Wang & Lan, 2000). Based on eight haplotypes of M.putaoensis obtained in this study, combined with the ND4LND4 sequences of six other muntjacs downloaded from GenBank, our phylogenetic results showed M. putaoensis as a member of genus Muntiacus (Figure 2A), consistent with the above mitochondrial genome findings.

    However, when M. truongsonensis and M. rooseveltorm were included in analyses based on cyt b (Figure 2B), we found that M. truongsonensis, rather than M. vuquangensis, was most closely related to M. putaoensis. Furthermore, we observed only weak support (bootstrap supports for ML/posterior probability in BI=71/0.60; Figure 2B) among these three species,and that they all likely belonged to the M. rooseveltorum complex of species. This result is consistent with previous morphological studies (James et al., 2008; Timmins & Duckworth,2016).

    Complete mitogenome data of M. putaoensis would be fundamental to research phylogentic relationship and conservation genetics of muntjac genus. As the current conclusion is based on analysis of only a few mtDNA fragments, comparison of longer sequences or genomes will be more convincing. More studies based on genome are required to better illuminate wholesome view of phylogenetics and evolution in this extraordinary taxon.

    ACKNOWLEDGEMENTS

    The authors sincerely thank Prof. Shu-Qiang Li and Dr. Yu Song for their support with the laboratory work. We also thank Mr. Kyaw Kyaw Han and Shuyin Huang for sample collection, thank Ms. Nan Sun and Alice Hughes for proofreading the manuscript. This work was possible because of the support from the Forest Research Institute, Hponkan Razi Wildlife Sanctuary and Hkakaborazi National Park Offices, Myanmar.

    Amato G, Egan MG, Rabinowitz A. 1999. A new species of muntjac,Muntiacus putaoensis (Artiodactyla: Cervidae) from northern Myanmar.Animal Conservation, 2(1): 1-7.

    Choudhury A. 2007. Discovery of leaf deer Muntiacus putaoensis in Nagaland with a new northerly record from Arunachal Pradesh. Journal of the Bombay Natural History Society, 104: 205-208.

    Choudhury A. 2009. Records and distribution of Gongshan and leaf muntjacs in India. Deer Specialist Group News, 23: 2-7.

    Darriba D, Taboada GL, Doallo R, Posada D. 2012. jModelTest 2: more models, new heuristics and parallel computing. Nature Methods, 9(8): 772.Datta A, Pansa J, Madhusudan MD, Mishra C. 2003. Discovery of the leaf deer Muntiacus putaoensis in Arunachal Pradesh: An addition to the large mammals of India. Current Science, 84(3): 454-458.

    Fontana F, Rubini M. 1990. Chromosomal evolution in Cervidae.Biosystems, 24(2): 157-174.

    Giao PM, Tuoc D, Dung VV, Wikramanayake ED, Amato G, Arctander P,MacKinnon JR. 1998. Description of Muntiacus truongsonensis, a new species of muntjac (Artiodactyla: Muntiacidae) from central Vietnam, and implications for conservation. Animal Conservation, 1(1): 61-68.

    Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O.2010. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic Biology,59(3): 307-321.

    Hassanin A, Delsuc F, Ropiquet A, Hammer C, Van Vuuren BJ, Matthee C,Ruiz-Garcia M, Catzeflis F, Areskoug V, Nguyen TT. 2012. Pattern and timing of diversification of Cetartiodactyla (Mammalia, Laurasiatheria), as revealed by a comprehensive analysis of mitochondrial genomes. Comptes Rendus Biologies, 335(1): 32-50.

    Irwin DM, Kocher TD, Wilson AC. 1991. Evolution of the cytochrome b gene of mammals. Journal of Molecular Evolution, 32(2): 128-144.

    James J, Ramakrishnan U, Datta A. 2008. Molecular evidence for the occurrence of the leaf deer Muntiacus putaoensis in Arunachal Pradesh,north-east India. Conservation Genetics, 9(4): 927-931.

    Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA,McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD,Gibson TJ, Higgins DG. 2007. Clustal W and clustal X version 2.0.Bioinformatics, 23(21): 2947-2948.

    Librado P, Rozas J. 2009. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25(11): 1451-1452.

    Martins RF, Fickel J, Le M, Van Nguyen T, Nguyen HM, Timmins R, Gan HM, Rovie-Ryan JJ, Lenz D, Forster DW, Wilting A. 2017. Phylogeography of red muntjacs reveals three distinct mitochondrial lineages. BMC Evolutionary Biology, 17: 34.

    Rabinowitz A, Khaing ST. 1998. Status of selected mammal species in North Myanmar. Oryx, 32(3): 201-208.

    Rabinowitz A, Myint T, Khaing ST, Rabinowitz S. 1999. Description of the leaf deer (Muntiacus putaoensis), a new species of muntjac from northern Myanmar. Journal of Zoology, 249(4): 427-435.

    Rao M, Htun S, Zaw T, Myint T. 2010. Hunting, livelihoods and declining wildlife in the Hponkanrazi Wildlife Sanctuary, North Myanmar.Environmental Management, 46(2): 143-153.

    Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Hohna S,Larget B, Liu L, Suchard MA, Huelsenbeck JP. 2012. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Systematic Biology, 61(3): 539-542.

    Schaller GB, Rabinowitz A. 2004. Species of barking deer (genus Muntiacus) in the eastern Himalayan region. The Journal of the Bombay Natural History Society, 101(3): 442-444.

    Shi YF, Shan XN, Li J, Zhang XM, Zhang HJ. 2003. Sequence and organization of the complete mitochondrial genome of the Indian muntjac(Muntiacus muntjak). Acta Zoologica Sinica, 49(5): 629-636. (in Chinese)

    Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. 2011.MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods.Molecular Biology and Evolution, 28(10): 2731-2739.

    Timmins R J, Duckworth JW. 2016. Muntiacus putaoensis. The IUCN Red List of Threatened Species 2016: e.T136479A22159478.http://dx.doi.org/10.2305/IUCN.UK.2016-1.RLTS.T136479A22159478.en

    Tuoc D, Dung VV, Dawson S, Arctander P, MacKinnon J. 1994. Introduction of a new large mammal species in Vietnam. Science and Technology News,3: 4-12.

    Wang W, Lan H. 2000. Rapid and parallel chromosomal number reductions in muntjac deer inferred from mitochondrial DNA phylogeny. Molecular Biology and Evolution, 17(9): 1326-1333.

    Wu HL, Fang SG. 2005. Mitochondrial DNA genetic diversity of black muntjac (Muntiacus crinifrons), an endangered species endemic to China.Biochemical Genetics, 43(7-8): 407-416.

    Wyman SK, Jansen RK, Boore JL. 2004. Automatic annotation of organellar genomes with DOGMA. Bioinformatics, 20(17): 3252-3255.

    Zhang HJ, Li J, Shi YF, Zhang XM, Xu CH, Shan XN. 2004a. Complete sequence of black muntjac (Muntiacus crinifrons) mitochondrial genome.Chinese Journal of Biochemistry and Molecular Biology, 20(4): 513-518. (in Chinese)

    Zhang XM, Shan XN, Shi YF, Zhang HJ, Li J, Zheng AL. 2004b. Sequence and organization of Muntiacus reevesi mitochondrial genome. Hereditas(Beijing), 26(6): 849-853.

    10 August 2017; Accepted: 07 September 2017

    s: This work was supported by a grant from the Southeast Asia Biodiversity Research Institute, Chinese Academy of Sciences (Y4ZK111B01)

    *Corresponding author, E-mail: quanrc@xtbg.ac.cn

    10.24272/j.issn.2095-8137.2017.058

    两个人看的免费小视频| 日本熟妇午夜| 好男人电影高清在线观看| 日本 av在线| 天堂√8在线中文| 免费人成视频x8x8入口观看| 免费在线观看成人毛片| a级一级毛片免费在线观看| 在线观看66精品国产| 亚洲av成人不卡在线观看播放网| 91麻豆精品激情在线观看国产| 国内毛片毛片毛片毛片毛片| 午夜激情福利司机影院| a级一级毛片免费在线观看| 亚洲 欧美 日韩 在线 免费| 欧美黄色淫秽网站| 中出人妻视频一区二区| 国产免费男女视频| 此物有八面人人有两片| 国产国拍精品亚洲av在线观看 | 一本久久中文字幕| 亚洲av电影不卡..在线观看| 久久精品国产综合久久久| 日本黄大片高清| 亚洲成人中文字幕在线播放| 中文字幕熟女人妻在线| 一级黄片播放器| 手机成人av网站| 亚洲精品在线观看二区| 久久久国产精品麻豆| 亚洲最大成人手机在线| 国产av一区在线观看免费| 国产伦在线观看视频一区| 成人特级黄色片久久久久久久| 亚洲自拍偷在线| 久久香蕉国产精品| 欧美日韩福利视频一区二区| 国产亚洲欧美98| 国产精品久久久久久久久免 | 波多野结衣高清作品| 亚洲av免费在线观看| 哪里可以看免费的av片| 国产精品爽爽va在线观看网站| 免费大片18禁| а√天堂www在线а√下载| 中国美女看黄片| 亚洲av免费高清在线观看| aaaaa片日本免费| 免费看日本二区| 中文字幕人妻熟人妻熟丝袜美 | 偷拍熟女少妇极品色| 18美女黄网站色大片免费观看| 非洲黑人性xxxx精品又粗又长| 国产成人a区在线观看| 欧美最新免费一区二区三区 | 亚洲激情在线av| 色播亚洲综合网| 亚洲精品456在线播放app | 两个人看的免费小视频| 人妻夜夜爽99麻豆av| 久久精品国产亚洲av香蕉五月| 18禁美女被吸乳视频| 一级黄片播放器| 中文字幕高清在线视频| 欧美又色又爽又黄视频| 国内毛片毛片毛片毛片毛片| e午夜精品久久久久久久| 免费一级毛片在线播放高清视频| 99国产极品粉嫩在线观看| 国产亚洲欧美在线一区二区| 国内精品久久久久久久电影| 成人永久免费在线观看视频| 精品久久久久久,| 亚洲狠狠婷婷综合久久图片| 亚洲人与动物交配视频| netflix在线观看网站| 又粗又爽又猛毛片免费看| 国产精品一及| 国产精华一区二区三区| 久久人妻av系列| 日韩精品中文字幕看吧| 国产免费av片在线观看野外av| 欧美zozozo另类| 欧美一级毛片孕妇| 香蕉久久夜色| 亚洲成人免费电影在线观看| 久久久久久人人人人人| 亚洲av成人不卡在线观看播放网| 亚洲精品一区av在线观看| 亚洲乱码一区二区免费版| 欧美bdsm另类| 亚洲av成人av| 制服丝袜大香蕉在线| 最后的刺客免费高清国语| 女警被强在线播放| 欧洲精品卡2卡3卡4卡5卡区| 国产伦精品一区二区三区四那| 欧美日韩瑟瑟在线播放| 麻豆久久精品国产亚洲av| 精品久久久久久,| 欧美黑人欧美精品刺激| 免费观看人在逋| 两人在一起打扑克的视频| 91九色精品人成在线观看| 精品人妻1区二区| 免费看美女性在线毛片视频| 国产v大片淫在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲自拍偷在线| 亚洲欧美一区二区三区黑人| 精品99又大又爽又粗少妇毛片 | 成人欧美大片| 亚洲美女黄片视频| 天天躁日日操中文字幕| 天堂网av新在线| 欧美日韩综合久久久久久 | 国产黄片美女视频| 国产免费av片在线观看野外av| 国产成人影院久久av| 制服人妻中文乱码| 精品国产超薄肉色丝袜足j| 男女床上黄色一级片免费看| 免费在线观看成人毛片| 少妇熟女aⅴ在线视频| 国产一区二区在线观看日韩 | 在线观看舔阴道视频| 日本与韩国留学比较| 性色av乱码一区二区三区2| 国产精品久久久久久久久免 | 久久香蕉精品热| 亚洲国产欧洲综合997久久,| 精品无人区乱码1区二区| 91九色精品人成在线观看| 村上凉子中文字幕在线| а√天堂www在线а√下载| 国产成人啪精品午夜网站| 国产老妇女一区| 国产高清视频在线观看网站| 美女高潮的动态| 两性午夜刺激爽爽歪歪视频在线观看| 免费av毛片视频| 亚洲欧美日韩高清专用| 特级一级黄色大片| 757午夜福利合集在线观看| 在线视频色国产色| 俺也久久电影网| 一级毛片女人18水好多| 少妇人妻一区二区三区视频| 一个人免费在线观看电影| 亚洲av电影不卡..在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久99热这里只有精品18| 高清毛片免费观看视频网站| 中文字幕高清在线视频| 免费看a级黄色片| 日本精品一区二区三区蜜桃| 久久久久久久精品吃奶| 嫁个100分男人电影在线观看| 99精品久久久久人妻精品| 老司机福利观看| 一区二区三区国产精品乱码| 日本黄色视频三级网站网址| 五月伊人婷婷丁香| 综合色av麻豆| 免费观看的影片在线观看| 老司机午夜福利在线观看视频| 人妻丰满熟妇av一区二区三区| 成人高潮视频无遮挡免费网站| 全区人妻精品视频| 久久精品夜夜夜夜夜久久蜜豆| a级毛片a级免费在线| 精品欧美国产一区二区三| 高潮久久久久久久久久久不卡| 好男人电影高清在线观看| 制服丝袜大香蕉在线| 成人18禁在线播放| av中文乱码字幕在线| 91久久精品国产一区二区成人 | 1000部很黄的大片| 日韩欧美三级三区| 国产精品精品国产色婷婷| 日本 欧美在线| 亚洲男人的天堂狠狠| 国产伦一二天堂av在线观看| 国产高潮美女av| 搡女人真爽免费视频火全软件 | 午夜免费观看网址| 久久6这里有精品| 国产视频内射| 精品久久久久久久人妻蜜臀av| h日本视频在线播放| 无遮挡黄片免费观看| 天堂动漫精品| 在线观看午夜福利视频| 亚洲国产精品久久男人天堂| 欧美绝顶高潮抽搐喷水| 亚洲 国产 在线| 99在线人妻在线中文字幕| 午夜免费激情av| 国产69精品久久久久777片| 网址你懂的国产日韩在线| www.999成人在线观看| 夜夜夜夜夜久久久久| 欧美乱色亚洲激情| 长腿黑丝高跟| 精品99又大又爽又粗少妇毛片 | 欧美+日韩+精品| 久久久久久久久大av| 操出白浆在线播放| 久久香蕉精品热| 午夜福利免费观看在线| 制服人妻中文乱码| 精华霜和精华液先用哪个| 99riav亚洲国产免费| 麻豆一二三区av精品| 丁香欧美五月| 国产成人aa在线观看| 国产免费男女视频| 国产黄a三级三级三级人| 老司机午夜福利在线观看视频| 成年女人看的毛片在线观看| 亚洲欧美日韩卡通动漫| 亚洲午夜理论影院| 美女黄网站色视频| 色综合欧美亚洲国产小说| 亚洲精品在线美女| 亚洲aⅴ乱码一区二区在线播放| 99热这里只有是精品50| 亚洲精品日韩av片在线观看 | 国产精品乱码一区二三区的特点| 色吧在线观看| 精品99又大又爽又粗少妇毛片 | 国产精品99久久99久久久不卡| 亚洲精品粉嫩美女一区| 高清在线国产一区| 精华霜和精华液先用哪个| 久久久色成人| 麻豆成人av在线观看| 国产精品日韩av在线免费观看| 色av中文字幕| 国产三级中文精品| 欧美成狂野欧美在线观看| 国产熟女xx| 久久久久久久久大av| 免费在线观看亚洲国产| 精品一区二区三区视频在线观看免费| 狠狠狠狠99中文字幕| av黄色大香蕉| 在线观看一区二区三区| 国产成人av教育| 一个人看的www免费观看视频| 国产国拍精品亚洲av在线观看 | 婷婷精品国产亚洲av在线| 嫩草影院入口| 国产免费男女视频| 露出奶头的视频| 一本精品99久久精品77| 国产一区在线观看成人免费| av片东京热男人的天堂| 国产精品综合久久久久久久免费| 一区二区三区免费毛片| 一个人看视频在线观看www免费 | 亚洲成av人片在线播放无| 12—13女人毛片做爰片一| 九九热线精品视视频播放| 国产v大片淫在线免费观看| 国产探花极品一区二区| 男女床上黄色一级片免费看| 久久久久精品国产欧美久久久| 欧美+亚洲+日韩+国产| 精品久久久久久久末码| 国产免费男女视频| 嫩草影院精品99| 最新美女视频免费是黄的| 亚洲精品一区av在线观看| 免费观看的影片在线观看| 国产精品乱码一区二三区的特点| 国产精品日韩av在线免费观看| 在线播放国产精品三级| 色噜噜av男人的天堂激情| 香蕉av资源在线| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| 少妇裸体淫交视频免费看高清| 最近最新中文字幕大全电影3| av在线蜜桃| 中亚洲国语对白在线视频| 亚洲无线观看免费| 久久精品综合一区二区三区| 国产色爽女视频免费观看| 中文字幕av成人在线电影| 午夜老司机福利剧场| 国产精品亚洲美女久久久| 久久香蕉精品热| 久久精品国产亚洲av香蕉五月| 十八禁人妻一区二区| 日本在线视频免费播放| 精品不卡国产一区二区三区| 99精品欧美一区二区三区四区| 亚洲欧美日韩高清专用| 看片在线看免费视频| 国产野战对白在线观看| 69av精品久久久久久| 高清在线国产一区| 国产99白浆流出| 亚洲av日韩精品久久久久久密| 成年免费大片在线观看| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产亚洲av香蕉五月| 男插女下体视频免费在线播放| 国产精品电影一区二区三区| 村上凉子中文字幕在线| 亚洲欧美激情综合另类| 一区二区三区高清视频在线| 男人的好看免费观看在线视频| 国产精品一区二区三区四区免费观看 | 他把我摸到了高潮在线观看| 成人鲁丝片一二三区免费| 日本免费一区二区三区高清不卡| 高清毛片免费观看视频网站| 在线观看免费视频日本深夜| 两个人视频免费观看高清| 91在线观看av| 淫秽高清视频在线观看| 99精品在免费线老司机午夜| 少妇的逼好多水| 尤物成人国产欧美一区二区三区| 国产伦精品一区二区三区视频9 | 日本与韩国留学比较| 操出白浆在线播放| 亚洲精品在线美女| 午夜福利视频1000在线观看| 国产精品久久久久久久电影 | 日韩欧美 国产精品| 一本久久中文字幕| 亚洲中文字幕日韩| 99riav亚洲国产免费| 亚洲av五月六月丁香网| 女人十人毛片免费观看3o分钟| 免费av观看视频| 亚洲熟妇熟女久久| 国产探花在线观看一区二区| 天堂动漫精品| 精品一区二区三区av网在线观看| 国内揄拍国产精品人妻在线| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣高清无吗| 亚洲国产色片| 国产一区二区三区视频了| 成人无遮挡网站| 亚洲av五月六月丁香网| 母亲3免费完整高清在线观看| xxx96com| 男人舔奶头视频| 在线观看美女被高潮喷水网站 | 俄罗斯特黄特色一大片| 欧美另类亚洲清纯唯美| 色哟哟哟哟哟哟| 国产精品自产拍在线观看55亚洲| 亚洲精品456在线播放app | 成年免费大片在线观看| 黄色日韩在线| 久久久久久久午夜电影| 桃色一区二区三区在线观看| av天堂中文字幕网| 久久久国产精品麻豆| 欧洲精品卡2卡3卡4卡5卡区| 欧美三级亚洲精品| 欧美性猛交╳xxx乱大交人| 天堂影院成人在线观看| 国产精品影院久久| www.www免费av| 亚洲国产欧美网| 国产成年人精品一区二区| 亚洲天堂国产精品一区在线| 亚洲av美国av| 国产成人a区在线观看| 黄片小视频在线播放| 一本精品99久久精品77| 亚洲av免费在线观看| 最近最新中文字幕大全电影3| 丰满的人妻完整版| 国产精品 欧美亚洲| 一区福利在线观看| 色av中文字幕| 久久99热这里只有精品18| 国产精品 欧美亚洲| 男女做爰动态图高潮gif福利片| 国产精品久久电影中文字幕| 他把我摸到了高潮在线观看| 国产精品一区二区免费欧美| 亚洲激情在线av| 制服人妻中文乱码| 亚洲熟妇中文字幕五十中出| 国产伦精品一区二区三区四那| 亚洲精品国产精品久久久不卡| 高潮久久久久久久久久久不卡| 十八禁网站免费在线| 国产精品一区二区免费欧美| 亚洲内射少妇av| av专区在线播放| 久久久久久国产a免费观看| 窝窝影院91人妻| 女生性感内裤真人,穿戴方法视频| 窝窝影院91人妻| 久久精品综合一区二区三区| 91九色精品人成在线观看| 深爱激情五月婷婷| 熟妇人妻久久中文字幕3abv| 丝袜美腿在线中文| 国产视频一区二区在线看| 男插女下体视频免费在线播放| 欧美三级亚洲精品| 久99久视频精品免费| 1000部很黄的大片| 精品无人区乱码1区二区| av中文乱码字幕在线| 成人鲁丝片一二三区免费| 国产成人系列免费观看| 亚洲不卡免费看| 国产蜜桃级精品一区二区三区| 级片在线观看| 一二三四社区在线视频社区8| 窝窝影院91人妻| 亚洲无线在线观看| 老司机深夜福利视频在线观看| 69人妻影院| 成年版毛片免费区| 亚洲五月婷婷丁香| 两个人的视频大全免费| 亚洲欧美激情综合另类| 五月伊人婷婷丁香| 国产高清videossex| 久久久久久久久久黄片| 国产探花在线观看一区二区| 亚洲av电影在线进入| 国产伦一二天堂av在线观看| 99久久无色码亚洲精品果冻| 欧美区成人在线视频| 午夜免费激情av| 人妻夜夜爽99麻豆av| 女警被强在线播放| 日韩成人在线观看一区二区三区| 国产成人啪精品午夜网站| 久久精品国产清高在天天线| 99精品在免费线老司机午夜| 欧美性感艳星| 日本五十路高清| 亚洲精品久久国产高清桃花| 两性午夜刺激爽爽歪歪视频在线观看| 国产成人福利小说| 午夜福利成人在线免费观看| 精品一区二区三区人妻视频| 欧美性感艳星| 窝窝影院91人妻| 国产精品免费一区二区三区在线| 午夜福利成人在线免费观看| av中文乱码字幕在线| av女优亚洲男人天堂| 一级黄色大片毛片| 69人妻影院| 韩国av一区二区三区四区| 亚洲内射少妇av| 国产精品一区二区免费欧美| 色噜噜av男人的天堂激情| 国产精品亚洲av一区麻豆| 成人av一区二区三区在线看| 性色avwww在线观看| 手机成人av网站| 欧美性感艳星| 国产激情偷乱视频一区二区| 高清日韩中文字幕在线| 亚洲第一电影网av| 亚洲18禁久久av| 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯| 在线a可以看的网站| 午夜免费激情av| 免费在线观看影片大全网站| 97超级碰碰碰精品色视频在线观看| 3wmmmm亚洲av在线观看| 男女那种视频在线观看| 亚洲,欧美精品.| 中文字幕人妻熟人妻熟丝袜美 | 国产成人啪精品午夜网站| 亚洲,欧美精品.| av专区在线播放| 俄罗斯特黄特色一大片| 国产精品1区2区在线观看.| 午夜精品在线福利| 精品人妻1区二区| 校园春色视频在线观看| 国产高清视频在线观看网站| 亚洲无线观看免费| 天天躁日日操中文字幕| 国产高清视频在线播放一区| 成年女人毛片免费观看观看9| 国产伦精品一区二区三区视频9 | 99久久精品热视频| 国产精品一区二区免费欧美| 乱人视频在线观看| 手机成人av网站| 一本一本综合久久| 亚洲真实伦在线观看| 国产精品久久电影中文字幕| 男女床上黄色一级片免费看| 中文资源天堂在线| 一个人观看的视频www高清免费观看| 熟女电影av网| 欧美日韩一级在线毛片| 美女免费视频网站| 日本在线视频免费播放| 琪琪午夜伦伦电影理论片6080| 内射极品少妇av片p| 午夜福利在线观看吧| 日本黄大片高清| 国内精品久久久久久久电影| 久久99热这里只有精品18| 国内精品久久久久精免费| 国产精品精品国产色婷婷| 少妇高潮的动态图| 男人和女人高潮做爰伦理| 欧美又色又爽又黄视频| 国内久久婷婷六月综合欲色啪| 天堂动漫精品| 老司机深夜福利视频在线观看| 亚洲精品久久国产高清桃花| 国产精品久久久久久亚洲av鲁大| 精品午夜福利视频在线观看一区| 日韩欧美三级三区| 欧美乱码精品一区二区三区| 精品久久久久久久末码| 久久精品人妻少妇| 欧美极品一区二区三区四区| www.999成人在线观看| 久久久久免费精品人妻一区二区| 最近最新免费中文字幕在线| 人妻夜夜爽99麻豆av| 成年女人毛片免费观看观看9| 国产欧美日韩一区二区三| 麻豆一二三区av精品| 亚洲七黄色美女视频| 精品国产超薄肉色丝袜足j| x7x7x7水蜜桃| 一区福利在线观看| 脱女人内裤的视频| 99国产精品一区二区蜜桃av| 不卡一级毛片| 亚洲欧美精品综合久久99| 美女黄网站色视频| 国产成+人综合+亚洲专区| 国模一区二区三区四区视频| 久久久久亚洲av毛片大全| 欧美另类亚洲清纯唯美| 国产一区二区在线av高清观看| 在线十欧美十亚洲十日本专区| 久久午夜亚洲精品久久| 在线视频色国产色| 欧美一级a爱片免费观看看| av在线天堂中文字幕| 久久国产精品影院| 国产中年淑女户外野战色| 久久久成人免费电影| 国产亚洲精品综合一区在线观看| 中亚洲国语对白在线视频| 91字幕亚洲| 国内久久婷婷六月综合欲色啪| 久久精品综合一区二区三区| 波野结衣二区三区在线 | 在线观看免费视频日本深夜| 非洲黑人性xxxx精品又粗又长| 草草在线视频免费看| 日韩欧美在线乱码| 免费看十八禁软件| 男女午夜视频在线观看| 久久久久久久亚洲中文字幕 | 成人鲁丝片一二三区免费| 国产精品免费一区二区三区在线| 熟妇人妻久久中文字幕3abv| 一本一本综合久久| 久久精品国产自在天天线| 一级毛片女人18水好多| 亚洲国产精品成人综合色| 91久久精品电影网| 蜜桃久久精品国产亚洲av| 综合色av麻豆| 国产老妇女一区| 国产69精品久久久久777片| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩东京热| 亚洲人成电影免费在线| 久久亚洲精品不卡| 在线播放无遮挡| 日本一二三区视频观看| 男女午夜视频在线观看| 高潮久久久久久久久久久不卡| а√天堂www在线а√下载| 日韩免费av在线播放| 高清在线国产一区| 国产精品自产拍在线观看55亚洲| 精品国产超薄肉色丝袜足j| 亚洲第一欧美日韩一区二区三区| 床上黄色一级片| 亚洲av二区三区四区| 国内精品久久久久久久电影| 岛国视频午夜一区免费看| 国产亚洲精品一区二区www| 男人和女人高潮做爰伦理| 一区福利在线观看| 午夜免费男女啪啪视频观看 | 不卡一级毛片| 成熟少妇高潮喷水视频| a级毛片a级免费在线| 午夜两性在线视频| 91麻豆av在线|