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

    Using SSR Marker to Trace Chinese Shrimp Fenneropenaeus chinensis Released in Natural Sea — A Feasible Strategy for Assessment of Release Effect in Natural Resources Recovery Program

    2018-04-04 02:40:27WANGWeijiWANGMosangXIAOGuangxiaHUYulongSONGWenpingKONGJieJINXianshi
    漁業(yè)科學(xué)進(jìn)展 2018年1期
    關(guān)鍵詞:王偉微衛(wèi)星對(duì)蝦

    WANG Weiji, WANG Mosang,2, XIAO Guangxia, HU Yulong, SONG Wenping, KONG Jie, JIN Xianshi

    ?

    Using SSR Marker to Trace Chinese ShrimpReleased in Natural Sea — A Feasible Strategy for Assessment of Release Effect in Natural Resources Recovery Program

    WANG Weiji1, WANG Mosang1,2, XIAO Guangxia3, HU Yulong1, SONG Wenping3, KONG Jie1, JIN Xianshi1①

    (1. Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071; 2. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306; 3. Bohai Fisheries Research Institute of Tianjin, Tianjin 300457)

    Enhancement release has been proven effective in natural resources recovery of Chinese shrimpin the last several decades in China, however, to assess the effectiveness of enhancement release, we still need to develop a high-efficient and easy-operational method to replace those physical-tagging release method with labor intensive, size- and number-limited. In the present study, single (with maternal known) parentage identification using eight simple sequence repeat (SSR) markers genotype fingerprint was used to trace Chinese shrimp released in Bohai Bay in 2013. A total of 884 shrimp spawners were collected from two hatcheries in Tianjin City respectively after enhancement release of shrimp larvae in May 2013. A total of 844 shrimp samples were recaptured around the release location approximately 4 months after the shrimp larvae were released into the natural sea. Genotype data of 8 SSR loci of the 1,726 samples were used for maternal-offspring parentage identification using CEervus 3.0 software. The allele number in each locus ranged from 8 to 68 with an average value of 33.25, which produced the cumulative exclusion probability with one parent known of all these sight loci up to 99%. Among the 844 recaptured shrimp samples, 448 (♂:♀=212:232, gender information was lost for 4 samples) were successfully traced to their 337 maternal parents using a logarithm of odds (LOD) > 3.0 threshold. Among these 337 maternal parents, 253 had a single offspring, 62 had two offspring, 18 had three offspring, 3 had four offspring, and 1 had five offspring. For the first time, a large number of released shrimp were identified from recapture samples, and this study showed that it is possible to trace all released Chinese shrimp without using any type of physical tag in enhancement release activities. This not only means more precise recapture ratio assessment than ever expected, but also this method demonstrates an effective method for large-scale hatchery release as well as for organisms used in hatchery enhancement which are not suitable for physical tagging.

    ; Parentage identification; SSR; Enhancement release

    Introduction

    One of the most important native marine shrimp species in China, with smaller populations also found along the west coast of the Korean Peninsula, the Chinese fleshy shrimp () has played an indispensable role in Chinese fishery and marine culture (Deng, 1998). However, since the 1980s, the once plentiful captures of Chinese shrimp in the Yellow and Bohai Seas has dramatically declined. This decline has been attributed to an increase in fishing, pollution and natural habitat destruction. Additionally, the traditional spring shrimp fishing season disappeared approximately in 1986 (Deng, 2001). After several decades of continuous practice, hatchery release enhancement was proven to be an alternative strategy for natural population recovery, with the fishing industry also gradually benefiting from the strategy (Wang, 2006; Li, 2012). During this period, different methods were developed to assess the recapture rate. These methods mainly include marker- recapture methods with physical markers such as eye tags, visible implant elastomers (VIE, http://www.nmt. us), code wire tags (CWT, http://www.nmt.us) and radio frequency identification (RFID, http://www.biomark. com) (Klima, 1965; Luo, 2008; Gao, 2014). However, they do have deficiencies in common, such as labor-intensive tagging processes, damage to tagged individuals, and impracticality for large-scale application, all of which make it difficult to calculate precise recapture rates for assessing the effect of the release. Therefore, there is an urgency to develop a new marker for precise recapture rate assessments.

    Wang(2014) introduced a new method that was based on individual identification and pedigree tracing using simple sequence repeat (SSR) markers, and it was successfully applied in release effect assessments in Jiaozhou Bay and Bohai Bay in China in 2012. In that study, a large number of shrimp from several full-sib families were replaced and functioned as the previously physically tagged individuals, thereby skipping the labor-intensive tagging process and significantly increasing the tagged individuals’ number at the same time, consequently producing a more precise recapture rate. However, during this “molecular marker-recapture” process, several problems emerged, particularly regarding family development, which was tedious and prone to failure, and the number of shrimp carrying the “molecular marker” was still finite compared to the whole released stock (Wang, 2014). This prompted us to develop a more thorough and precise assessment method, where all released shrimp could be traced through parentage identification using an SSR marker. A previous study has verified the accuracy of SSR markers on single- parent parentage identification in Chinese shrimp (Zhang, 2014). In the current study, all of the 884 spawners from the two hatcheries responsible for larvae hatchery release in Tianjin City were collected and involved in a four-month-later single-parent parentage identification using 844 recaptured samples and genotype information from 8 SSR loci. For the first time, all offspring from those 884 spawners were identified from recapture samples; further, we could theoretically expect that all released shrimp could also be traced.

    Materials and methods

    Sample collection

    In early March 2013, approximately 2000 shrimp spawners (sexual matured female shrimp carrying spermatophore)were collected along their spawning migration route in the Yellow Sea. The spawners were cultured at the Bohai Fishery Enhancement Station of the Bohai Fisheries Research Institute of Tianjin and the Dashentang Shrimp Seedling Plant in Tianjin. These stations are the two exclusively designated hatcheries for Chinese shrimp larvae production used for enhancement release in Bohai Bay in Tianjin City. After one and half month indoor domestication and nutrient enrichment period, a total of 884 spawners in both hatcheries started a 15-day spawning period. During this time, any dead individuals were collected and stored at –20℃. After the spawning period, all 884 spawners (successful and unsuccessful) were numbered and transported at –20℃ to the laboratory in the Yellow Sea Fisheries Research Institute (YSFRI), Qingdao, P. R. China and stored at –75℃ until use.

    In early June, after a 45-day growth period, the fertilized eggs in the two hatcheries had developed into the post larval phase with body lengths of 1.2 cm, and an estimated 1.487 billion post-larvae were released into the Bohai Sea near Tianjin City. By the end of August, just before the natural Chinese shrimp population initiated their wintering migration, 844 shrimp were recaptured from the area around the original release location. The recaptured shrimp were transported alive to the YSFRI laboratory. The body length, weight, and gender were recorded for further analysis. All samples were maintained at –75℃ for DNA isolation.

    Methods

    The swimming apparatus of each shrimp sample was dissected for genomic DNA isolation using the protocol described by Wang (2008). Pure DNA was dissolved in double-distilled water to obtain a final concentration of 50 ng/μl. A total of 8 unlinked mutually SSR loci were used in genotyping all DNA samples, among which, 6 were the same as those used in previous reports by Zhang(2014) and Wang(2014). The first four of these six were EN0032, RS0622, FCKR002, and FCKR013, which formed a high-temperature group (HTG) of the quadruple SSR-PCR system to save PCR and genotyping processing time. The other two were RS1101 and FCKR005. In the present study, two new SSR loci, FC019 and RS0683, were introduced to replace the previous FCKR005 and FCKR007 loci for more ideal genotyping results. The PCR reaction components and procedures were adopted according to the description of Li(2012), and all other SSR primer details were displayed in Tab. 1 (Li, 2012; Dong, 2014). An ABI 3130 auto-sequencer (Applied Biosystems Inc.) was used for the PCR amplification product separation. Allele size estimation and recording were completed using the DNA marker GeneScan-500 Liz (Applied Biosystems Inc.) and Genemapper software (Applied Biosystems Inc.). Cervus 3.0 software (http://www. fieldgenetics.com), a computer program for the assignment of parents to their offspring using genetic markers, was used for single-parent (maternal only) parentage identification.

    Tab.1 Details of 4 SSR markers composing one set of quadruple-fold fluorescent-labeled PCR system and four single SSR-PCR systems

    Note: 6-FAM, ROX, HEX and TAMRA were the labeled fluorescent dyes

    Results

    The genotyping results of 8 SSR loci indicated that 1728 shrimp samples, including 884 spawners and 844 recapture samples, displayed high genetic variability. The allele number at each locus ranged from 8 to 68 with an average value of 33.25 (Tab.2). The polymorphic information content (PIC) values were 0.776~0.951 with an average value of 0.857. Further more, the cumulative exclusion probability with one parent known [combined non-exclusion probability (NE-1P): 8.913′10?5] of all 8 SSR loci was in excess of 0.99% (Tab. 2), with which high reliable parentage identification results were expected.

    Subsequently, 1728 samples were involved in single-parent parentage identification analysis between the 884 spawners and 844 recapture samples using the 8 SSR loci and Cervus 3.0 software. The results indicate that a total of 448 offspring from the 844 recapture samples were traced to 337 of the 884 spawners with a LOD value>3.0. The identified maternal-offspring relationship was comprised of the six following cases: among most of traced spawners (253), only one offspring was detected; in 62 traced spawners, two offspring were detected; in 18 traced spawners, three offspring were detected; in 3 traced spawners, four offspring were detected; and in one traced spawner, one offspring was detected.

    Discussion

    Parentage identification

    With both parents known, only several SSR loci (given there are abundant alleles in each locus) could produce a high enough cumulative exclusion probability to meet authentic parentage identification, even in high-capacity sample analysis, as shown in previous research (Li, 2012; Dong, 2014). However, in practice, it is not always the case that both parents are simultaneously available, especially in some crustacean species that mature asynchronously. In the instance of Chinese shrimp, the male first transfers the spermatophore to the female during winter migration and then dies in the following autumn. In the following spring, the females bearing spermatophores simultaneously produce eggs and release the matured sperm from spermatophore at the spawning sites. In this case, single parentage identification using more SSR markers than used when oth parents are known is an alternative for tracing released individuals or other related studies. The accuracy of parentage identification using SSR markers has been confirmed by an indoor sealed pond experiment using physically tagged individuals (VIE) as the reference (Li, 2012). The accuracy of maternal only parentage identification was indirectly confirmed by two enhancement release studies of Chinese shrimp in Bohai Bay and Jiaozhou Bay in 2012 (Wang, 2014). In the present study, two SSR loci from the previous study, FCKR005 and FCKR007, were replaced with two newly introduced loci, FC019 and RS0683 for better genotyping resolution (Wang, 2014; Dong, 2014). There were 43 and 8 alleles detected at FC019 and RS0683, respectively. The sum of alleles is slightly less than what was detected at loci FCKR005 and FCKR007 in the previous study, indicating that the cumulative exclusion probability produced by the present 8 SSR loci is slightly less than or equal to that produced by the previous 8 SSR loci. Moreover, in the present study, at each assignment where the 8 SSR loci matched between the candidate parent and offspring, only when the locus-by-locus logarithm (base 10) of odds (LOD) score was≥3.0. Could it be declared as a true parentage relationship?

    Tab.2 Summary of genetic variability of 8 SSR loci crossing the whole 1728 samples and exclusion probability they could offer

    Note:: Number of alleles at the locus;: Number of individuals typed at the locus;Obs: Observed heterozygosity;Exp: Expected heterozygosity; NE-1P: Average non-exclusion probability for one candidate parent; NE-2P: Average non-exclusion probability for one parent given the genotype of a known parent of the opposite sex; NE-PP: Average non-exclusion probability for a candidate parent pair

    Previous studies have shown that where the LOD score was not considered in the parentage assignment, even when alleles at each locus from offspring matched the parents or a single parent, the chance of a false positive would significantly increase. In regards to the precision of tracing released individuals using SSR markers, the following measures were suggested in the study by Bravington(2004): (i) Use as few brood- stock as possible. Generally, a wild Chinese shrimp spawner can produce 400,000 to 1,800,000 fertilized eggs, among which, 40%~60% could develop to the post larvae phase in plants and be released into the sea. In practice, plants fishing shrimp spawners along their spawning migration route at each spring for shrimp larvae production. Then, for a better profit, they usually produce as many post larvae as possible using as few spawners as possible, which in fact increases the chance to produce authentic parentage identification results, meanwhile increasing the risk of inbreeding depression. (ii) Obtain paternal as well as maternal genotypes. Chinese shrimp are a type of asynchronous-maturing crustacean, in other words, the male transfers the spermatophore to the female and subsequently dies. The female bearing the spermatophore simultaneously releases the egg and the sperm the next spring after they have finished their over-wintering phase. Removing the spermatophore directly would seriously damage the female. Therefore, in the present study, only the genotype information from the female was involved in the parentage identification analysis. (iii) Stock as many as possible, are similar to premise (i). In fact, the two shrimp seedling plants produced and released approximately 1.2×109post larvae shrimp using only approximately 1,200 captured spawners in the spring of 2012. (iv) Finally, capture and genotype are as many as prawns as possible. The number of recaptured samples in this study was 844, among which, a total of 448 released individuals were identified, representing two large sample sizes. On account of the above-mentioned measures, the accuracy of the parentage identification results in this study is reliable.

    Release effect assessment

    Using SSR markers to trace released individuals played an important part in assessing the accuracy of the effect of hatchery enhancement of Chinese shrimp. Due to the tiny size (the released shrimp post larvae were approximately 1.2 cm), and the large scale release number (approximately 0.5×109~1.5×109shrimp post larvae were released into Bohai Bay every year), there was no appropriate physical marker that could be used. This problem has persisted over several decades which was one of the most puzzling problems in the accuracy assessment of the effect of hatchery release since the Chinese shrimp hatchery release activity was initiated. Bravington(2004) evaluated the potential of using SSR markers for estimating the proportion of hatchery-reared offspring in an enhancement program, which prompted us to introduce the use of SSR markers in the assessment of the effect of Chinese shrimp hatchery enhancement. The first attempts were carried out in Jiaozhou Bay and Bohai Bay in 2012. In each bay, six full-sib families were released with a hatchery stock. During the late autumn fishing season, a total of 12 individuals from 7 families were identified from 5739 recaptured samples using SSR markers. This was the first time that molecular marker recapture assessment was applied to the Chinese shrimp hatchery enhancement (Wang, 2014). However, due to the small number of molecular-marked individuals released from the 12 full-sib families (but far more than previous studies that released physically tagged shrimp), the ratios between the molecular-marked individuals and the hatchery stock in Bohai Bay and Jiaozhou Bay were approximately 1:800 and 1:300, respectively. The reliability of the recapture ratio using molecular marker recapture in that study still needs be considered, although it produced a molecular marker recapture ratio very close to those of traditional methods(Wang, 2014).

    To improve the accuracy of molecular marker recapture assessment in Chinese shrimp, one strategy could be the release of all stock throughout full-family development. In that case, all released shrimp could be traced using SSR genotyping information under the condition that all parents were known. However, it was only theoretical, since large scale artificial fertilization is labor-extensive and easily fails. Additionally, mated female over-wintering culturesand specific family culture installations are only available at small plants. Since hatcheries release spawner shrimp every spring and the spawners that finished producing could be collected, the alternative is parentage identification with only maternal parents known. In the present study, the number of shrimp spawners used in the seedling development was close to 2000. The original plan was to sample all of these spawners after they finished spawning. However, for uncontrolled reasons, only 844 spawners were sampled for subsequent SSR genotyping, and only 337 individuals were identified as the previously released individuals from the 844 recaptured samples, with the remaining 505 individuals originating from both the wild population and the released stock. The genetic variation between the released stock and the wild population could not be further assessed in this study. This problem has been assessed in the Chinese shrimp hatchery enhancement field for the last thirty years.

    Before they initiate their over-wintering migration in the late autumn, all shrimp, including the wild population and the released stock concentrate and forage in a relatively fixed area. When all spawner shrimp used in the hatchery enhancement are sampled in an area (for example, Bohai Bay), then there is no difficulty in discriminating the wild population from the released stock in the recaptured samples. Additionally, It could also be clarified whether the released stock “dilute” the genetic diversity of the wild population (Kitada, 2009). Further, the number and proportion of the released stock in the following year’s spawner population could be accurately identified by investigation over consecutive years, which is also critical to evaluate the effect of the hatchery enhancement and has not yet been assessed. Chinese shrimp exhibit migration behavior over their entire life period, and their migration route spans northern latitudes from 33°(over-wintering) to 41°(spawning) (Deng, 1990). Deng(1990) investigated the distribution of the migration routes of the species, including spawning migration, foraging migration and over-wintering migration using physical tagging release assessment. Due to the limited number of physically tagged shrimp, however, the different geographic distribution of the population still needs to be evaluated. Using the parentage identification strategy in this study, this puzzle may be clarified in the future.

    Bravington MV, Ward RD. Microsatellite DNA markers: Evaluating their potential for estimating the proportion of hatchery-reared offspring in a stock enhancement programme. Molecular Ecology, 2004, 13(5): 1287–1297

    Deng JY, Ye CC, Liu YC. Penaeid shrimp—Its biology and management in the Yellow-Bohai Seas. Beijing: Ocean Press, 1990, 256–258 [鄧景耀, 葉昌臣, 劉永昌. 渤黃海的對(duì)蝦及其資源管理. 北京: 海洋出版社, 1990, 256–258]

    Deng JY, Zhuang ZM. The cause of recruitment variation ofin the Bohai Sea. Journal of Fishery Sciences of China, 2001, 7(4): 125–128 [鄧景耀, 莊志猛. 渤海對(duì)蝦補(bǔ)充量變動(dòng)原因的分析及對(duì)策研究. 中國(guó)水產(chǎn)科學(xué), 2001, 7(4): 125–128]

    Deng JY. The research status of shrimp fishery biology. Chinese Bulletin of Life Sciences, 1998, 10(4): 191–194, 1997 [鄧景耀. 對(duì)蝦漁業(yè)生物學(xué)研究現(xiàn)狀. 生命科學(xué), 1998, 10(4): 191–194, 1997]

    Dong SR, Luan S, Kong J. Screening of microsatellite DNA markers associated with growth traits in a singlefamily. Journal of Fishery Sciences of China, 2014, 21(5): 936–943 [董世瑞, 欒生, 孔杰. 中國(guó)明對(duì)蝦單個(gè)家系中與生長(zhǎng)性狀相關(guān)微衛(wèi)星標(biāo)記的初步篩選. 中國(guó)水產(chǎn)科學(xué), 2014, 21(5): 936–943]

    Gao H, Yan BL, Lai XF,. Progress on the identification technology used in stock enhancement for crustacean organisms. Transactions of Oceanology and Limnology. 2014(1): 94–100 [高煥, 閻斌倫, 賴曉芳, 等. 甲殼類生物增殖放流標(biāo)志技術(shù)研究進(jìn)展. 海洋湖沼通報(bào), 2014(1): 94–100]

    Kitada S, Shishidou H, Sugaya T,. Genetic effects of long-term stock enhancement programs. Aquaculture, 2009, 290(1–2): 69–79

    Klima EF. Evaluation of biological stains, inks, and fluorescent pigments as marks for shrimp. Washington: U.S. Department of the Interior, Bureau of Commercial Fisheries, 1965

    Li WY, Wang WJ, Kong J,. Establishment of the microsatellite quadruple PCR technology forand its application in the evaluation on the effect of releasing simulation. Acta Oceanologica Sinica, 2012, 34(5): 213–220 [李偉亞, 王偉繼, 孔杰, 等. 中國(guó)對(duì)蝦微衛(wèi)星四重PCR技術(shù)的建立及其在模擬放流效果評(píng)估方面的應(yīng)用. 海洋學(xué)報(bào), 2012, 34(5): 213–220]

    Li ZY, Wang J, Zhao ZL,. Resources enhancement ofin the Bohai Sea. Progress in Fishery Sciences, 2012, 33(3): 1–7 [李忠義, 王俊, 趙振良, 等. 渤海中國(guó)對(duì)蝦資源增殖調(diào)查. 漁業(yè)科學(xué)進(jìn)展, 2012, 33(3): 1–7]

    Luo K, Zhang TS, Kong J,. Study on techniques of internal tagging by visible implant fluorescent elastomer in. Marine Fisheries Research, 2008, 29(3): 48–52 [羅坤, 張?zhí)鞎r(shí), 孔杰, 等. 中國(guó)對(duì)蝦幼蝦熒光體內(nèi)標(biāo)記技術(shù)研究. 海洋水產(chǎn)研究, 2008, 29(3): 48–52]

    Wang QY, Zhuang ZM, Deng JY,. Stock enhancement and translocation of the shrimpin China. Fisheries Research, 2006, 80(1): 67–69

    Wang WJ, Zhang K, Luo K,. Assessment of recapture rates after hatchery release of Chinese shrimpin Jiaozhou Bay and Bohai Bay in 2012 using pedigree tracing based on SSR markers. Fisheries Science, 2014, 80(4): 749–755

    Wang WJ. Ⅰ Genetic mapping of the Chinese shrimpusing AFLP markers and commercial traits QTL mapping. Ⅱ Genetic linkage mapping using AFLP markers and primarily study on sex-determination of bluegill sunfish (). Doctoral Dissertation of Ocean University of China, 2008 [王偉繼. Ⅰ中國(guó)對(duì)蝦()AFLP分子標(biāo)記遺傳連鎖圖譜的構(gòu)建以及相關(guān)性狀的QTL定位分析Ⅱ藍(lán)鰓太陽魚()AFLP分子標(biāo)記遺傳連鎖圖譜的構(gòu)建及性別決定機(jī)制初探. 中國(guó)海洋大學(xué)博士研究生學(xué)位論文, 2008]

    Zhang K, Wang WJ, Kong J,. Accuracy of short sequence repeats on single-parent parentage identification in Chinese shrimp. Aquatic Biology, 2014, 20(1): 45–51

    (編輯 馬璀艷)

    10.11758/yykxjz.20170829001

    http://www.yykxjz.cn/

    * This work was supported by National Basic Research Program of China (973 Program)(2015CB453303) and International Science and Technology Cooperation Program of China (2013DFA31410). WANG Weiji, E-mail: wangwj@ysfri.ac.cn

    JIN Xianshi, E-mail: jin@ysfri.ac.cn

    2017-08-29,

    2017-09-15

    A

    Chinese Library Classification Number S931

    Article No. 2095-9869(2018)01-0021-06

    Wang WJ, Wang MS, Xiao GX, Hu YL, Song WP, Kong J, Jin XS.Using SSR marker to trace Chinese shrimpreleased in natural sea — A feasible strategy in release effect assessment in natural resources recovery program. Progress in Fishery Sciences, 2018, 39(1): 21–26

    猜你喜歡
    王偉微衛(wèi)星對(duì)蝦
    梅花鹿基因組微衛(wèi)星分布特征研究
    CHARACTERIZATION OF RESIDUATED LATTICES VIA MULTIPLIERS*
    對(duì)蝦養(yǎng)殖弱勢(shì)群體的管理
    對(duì)蝦吃料慢的原因分析和處理
    對(duì)蝦免疫增強(qiáng)劑研究進(jìn)展
    對(duì)蝦常見環(huán)境性疾病的防治
    小 蝌 蚪 的 尾 巴
    Convection: a neglected pathway for downward transfer of wind energy in the oceanic mixed layer*
    藝術(shù)百家 王偉
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    嘟嘟电影网在线观看| 日本猛色少妇xxxxx猛交久久| 国产av精品麻豆| 久久这里有精品视频免费| 久久久久久久久久久丰满| 中文字幕亚洲精品专区| 国产精品久久久久久久久免| 另类精品久久| 国产中年淑女户外野战色| 一级毛片久久久久久久久女| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 国产成人午夜福利电影在线观看| 2022亚洲国产成人精品| 黑人猛操日本美女一级片| 成人无遮挡网站| 日韩一区二区三区影片| 午夜精品国产一区二区电影| 亚洲成色77777| 人妻夜夜爽99麻豆av| 国产黄色免费在线视频| 青春草亚洲视频在线观看| videossex国产| 国产91av在线免费观看| 亚洲性久久影院| 亚洲无线观看免费| 美女中出高潮动态图| 多毛熟女@视频| 这个男人来自地球电影免费观看 | 国产一区亚洲一区在线观看| 国产成人精品婷婷| 亚洲av福利一区| 精品久久久久久电影网| 国产极品粉嫩免费观看在线 | 不卡视频在线观看欧美| 久久狼人影院| 日本91视频免费播放| 亚洲av不卡在线观看| 国产美女午夜福利| 欧美 亚洲 国产 日韩一| av国产久精品久网站免费入址| 国产高清有码在线观看视频| 卡戴珊不雅视频在线播放| 丰满迷人的少妇在线观看| 日韩欧美精品免费久久| 国产日韩欧美亚洲二区| 国产成人精品久久久久久| av福利片在线| 少妇裸体淫交视频免费看高清| 亚洲婷婷狠狠爱综合网| 成年av动漫网址| 国产欧美日韩一区二区三区在线 | 91在线精品国自产拍蜜月| 一本—道久久a久久精品蜜桃钙片| 国产一区二区三区av在线| 国产精品久久久久久久电影| 欧美高清成人免费视频www| 久久久a久久爽久久v久久| 免费黄网站久久成人精品| 十八禁高潮呻吟视频 | 22中文网久久字幕| 在线观看人妻少妇| 日本色播在线视频| 亚洲美女黄色视频免费看| 国产一区二区三区av在线| 国产精品无大码| 亚洲欧美精品专区久久| 十分钟在线观看高清视频www | 国产精品一区二区在线不卡| 午夜久久久在线观看| 黄色毛片三级朝国网站 | 精品国产露脸久久av麻豆| 男的添女的下面高潮视频| 一级毛片黄色毛片免费观看视频| 亚洲无线观看免费| 国产精品人妻久久久久久| 乱人伦中国视频| 国产精品蜜桃在线观看| 国产视频内射| 国产视频内射| 丝瓜视频免费看黄片| 国产一区有黄有色的免费视频| 在线亚洲精品国产二区图片欧美 | 性色av一级| 免费黄网站久久成人精品| 免费黄网站久久成人精品| 精品一区二区三卡| 日韩,欧美,国产一区二区三区| 99re6热这里在线精品视频| 国产高清三级在线| 久久久国产欧美日韩av| 中文字幕久久专区| 男人添女人高潮全过程视频| 国产欧美日韩精品一区二区| 丁香六月天网| 久久精品国产鲁丝片午夜精品| 青春草国产在线视频| 乱系列少妇在线播放| av不卡在线播放| 丰满少妇做爰视频| 免费观看的影片在线观看| 女性被躁到高潮视频| 免费黄色在线免费观看| 男人狂女人下面高潮的视频| 国产欧美日韩一区二区三区在线 | 亚洲av中文av极速乱| 免费av不卡在线播放| 黄色欧美视频在线观看| 国产探花极品一区二区| 亚洲第一区二区三区不卡| 街头女战士在线观看网站| 亚洲av男天堂| 欧美精品国产亚洲| 少妇的逼好多水| 成人漫画全彩无遮挡| 亚洲综合精品二区| 国产一级毛片在线| 99久久精品一区二区三区| 国产高清不卡午夜福利| 欧美精品一区二区大全| 3wmmmm亚洲av在线观看| 免费观看a级毛片全部| 黄片无遮挡物在线观看| 日韩不卡一区二区三区视频在线| 九九久久精品国产亚洲av麻豆| 精品久久久久久久久av| 在线观看国产h片| 高清欧美精品videossex| 熟女人妻精品中文字幕| 啦啦啦视频在线资源免费观看| 国产91av在线免费观看| 高清在线视频一区二区三区| 日本欧美视频一区| 亚洲,欧美,日韩| 久久女婷五月综合色啪小说| 新久久久久国产一级毛片| 亚洲第一av免费看| tube8黄色片| 欧美 日韩 精品 国产| 蜜桃久久精品国产亚洲av| 精品人妻一区二区三区麻豆| 精品一区二区三卡| 男人狂女人下面高潮的视频| 日本欧美视频一区| 人人妻人人爽人人添夜夜欢视频 | 人人妻人人澡人人爽人人夜夜| 永久网站在线| 亚洲欧美日韩东京热| 亚洲成人av在线免费| 婷婷色麻豆天堂久久| av天堂久久9| 日韩欧美精品免费久久| 午夜激情久久久久久久| 久久国产乱子免费精品| 国产黄频视频在线观看| 午夜av观看不卡| 国产黄色免费在线视频| 日本午夜av视频| 亚洲经典国产精华液单| 精品酒店卫生间| 人人妻人人澡人人爽人人夜夜| 夫妻性生交免费视频一级片| 美女大奶头黄色视频| 免费久久久久久久精品成人欧美视频 | 久热久热在线精品观看| 免费大片黄手机在线观看| 只有这里有精品99| 亚洲婷婷狠狠爱综合网| 久久久久久久国产电影| 高清黄色对白视频在线免费看 | 日日啪夜夜爽| 午夜免费男女啪啪视频观看| 一区二区三区精品91| 精品一区在线观看国产| 色视频www国产| 国产亚洲av片在线观看秒播厂| 一区二区三区免费毛片| 国产高清国产精品国产三级| 黑人高潮一二区| 男女啪啪激烈高潮av片| 你懂的网址亚洲精品在线观看| 久久久久久久久久成人| 男人添女人高潮全过程视频| 黑人巨大精品欧美一区二区蜜桃 | 777米奇影视久久| 精品少妇内射三级| 国产日韩欧美亚洲二区| 在线亚洲精品国产二区图片欧美 | 婷婷色综合www| 蜜桃在线观看..| 久久久久久久久久人人人人人人| 久久精品久久久久久噜噜老黄| 美女内射精品一级片tv| 高清视频免费观看一区二区| 人人妻人人澡人人看| 中国国产av一级| av黄色大香蕉| 久久免费观看电影| 男人舔奶头视频| 亚洲欧洲日产国产| 天堂俺去俺来也www色官网| 人人妻人人爽人人添夜夜欢视频 | 国内精品宾馆在线| 欧美高清成人免费视频www| 日本午夜av视频| 中文乱码字字幕精品一区二区三区| 热re99久久精品国产66热6| 久久久久国产网址| 丰满人妻一区二区三区视频av| 三级国产精品片| 两个人免费观看高清视频 | 一级av片app| 国产毛片在线视频| 免费观看a级毛片全部| 亚洲精品日本国产第一区| 国产av一区二区精品久久| 只有这里有精品99| 黄色毛片三级朝国网站 | 色吧在线观看| 亚洲精品中文字幕在线视频 | 人人妻人人澡人人看| 大片免费播放器 马上看| 国产午夜精品一二区理论片| 汤姆久久久久久久影院中文字幕| 美女福利国产在线| 亚洲国产精品999| 免费观看a级毛片全部| 中文精品一卡2卡3卡4更新| 性色avwww在线观看| 日本91视频免费播放| 欧美日韩综合久久久久久| 3wmmmm亚洲av在线观看| 精品国产一区二区久久| 丝袜在线中文字幕| 国产亚洲一区二区精品| 成人综合一区亚洲| 成人亚洲精品一区在线观看| 九草在线视频观看| .国产精品久久| 永久网站在线| 在线观看三级黄色| 亚洲综合精品二区| 丰满饥渴人妻一区二区三| 午夜免费男女啪啪视频观看| 久久国产精品男人的天堂亚洲 | 看免费成人av毛片| 亚洲av国产av综合av卡| 欧美bdsm另类| 亚洲内射少妇av| 麻豆乱淫一区二区| 国产成人一区二区在线| 午夜久久久在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲欧美日韩东京热| 桃花免费在线播放| 三上悠亚av全集在线观看 | 国产精品女同一区二区软件| 成人二区视频| 亚洲激情五月婷婷啪啪| 国产黄片视频在线免费观看| 日韩视频在线欧美| 免费观看a级毛片全部| 国产淫片久久久久久久久| 久久久久久久久久成人| 亚洲av福利一区| 黑人巨大精品欧美一区二区蜜桃 | 纵有疾风起免费观看全集完整版| 一本大道久久a久久精品| 欧美激情极品国产一区二区三区 | 一边亲一边摸免费视频| 久久久久久久久久久丰满| 亚洲美女黄色视频免费看| 国产精品久久久久久精品电影小说| 国产黄片美女视频| 日日爽夜夜爽网站| 99国产精品免费福利视频| 国产亚洲午夜精品一区二区久久| 一级片'在线观看视频| 蜜桃久久精品国产亚洲av| 亚洲欧美精品专区久久| 国产成人aa在线观看| 夜夜看夜夜爽夜夜摸| 免费观看无遮挡的男女| 熟女av电影| 九色成人免费人妻av| 亚洲国产精品国产精品| 久久国产乱子免费精品| 亚洲成人一二三区av| 男女边摸边吃奶| 日韩亚洲欧美综合| av.在线天堂| 午夜福利视频精品| 成年人免费黄色播放视频 | 在线观看免费高清a一片| av视频免费观看在线观看| 日韩大片免费观看网站| 亚洲精品456在线播放app| 久久久亚洲精品成人影院| 91精品伊人久久大香线蕉| 成人综合一区亚洲| 精品国产乱码久久久久久小说| 国产成人aa在线观看| 久久久久久久久久久免费av| av国产精品久久久久影院| 国产高清三级在线| 欧美日韩在线观看h| av女优亚洲男人天堂| 精品少妇黑人巨大在线播放| 日韩三级伦理在线观看| 欧美 日韩 精品 国产| 亚洲精品,欧美精品| 亚洲高清免费不卡视频| 两个人免费观看高清视频 | 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲综合精品二区| 水蜜桃什么品种好| 丁香六月天网| 人人妻人人澡人人爽人人夜夜| 免费观看a级毛片全部| 综合色丁香网| 久热久热在线精品观看| 国产在线免费精品| 爱豆传媒免费全集在线观看| 精品99又大又爽又粗少妇毛片| 国产精品一区二区性色av| 黄色怎么调成土黄色| 一个人免费看片子| 日日啪夜夜撸| 在线观看免费高清a一片| 欧美人与善性xxx| 在线观看三级黄色| 伊人久久国产一区二区| 久久久午夜欧美精品| 我要看日韩黄色一级片| 麻豆成人午夜福利视频| h视频一区二区三区| 亚洲精品自拍成人| 国产精品欧美亚洲77777| 我要看黄色一级片免费的| 国产高清三级在线| 免费观看在线日韩| 菩萨蛮人人尽说江南好唐韦庄| 狂野欧美激情性xxxx在线观看| 亚洲四区av| 男人添女人高潮全过程视频| 你懂的网址亚洲精品在线观看| 日本vs欧美在线观看视频 | 中文字幕人妻熟人妻熟丝袜美| 亚洲精品一二三| 成人黄色视频免费在线看| 久久久久网色| 黑丝袜美女国产一区| 国产精品欧美亚洲77777| 免费少妇av软件| av福利片在线| 亚洲精品成人av观看孕妇| 老司机影院毛片| 噜噜噜噜噜久久久久久91| 精品人妻一区二区三区麻豆| av免费观看日本| 国产 精品1| 亚洲成人av在线免费| 亚洲国产欧美在线一区| 国产精品国产三级专区第一集| 亚洲av成人精品一区久久| 国产男女超爽视频在线观看| 亚洲国产av新网站| 99久国产av精品国产电影| 免费人成在线观看视频色| 亚洲精品,欧美精品| 一级毛片电影观看| 精品久久久精品久久久| 国产精品久久久久成人av| 成人综合一区亚洲| 男人舔奶头视频| 高清不卡的av网站| 中文字幕人妻熟人妻熟丝袜美| 人妻夜夜爽99麻豆av| 成人漫画全彩无遮挡| 亚洲国产成人一精品久久久| 亚洲中文av在线| 又大又黄又爽视频免费| av不卡在线播放| 精品一品国产午夜福利视频| 日本午夜av视频| 中文字幕精品免费在线观看视频 | 人人妻人人添人人爽欧美一区卜| 国产精品麻豆人妻色哟哟久久| av.在线天堂| 亚洲av成人精品一二三区| av国产久精品久网站免费入址| 99热这里只有是精品50| av天堂中文字幕网| 夫妻性生交免费视频一级片| 日本vs欧美在线观看视频 | 国产av码专区亚洲av| 69精品国产乱码久久久| 一本色道久久久久久精品综合| 精品久久国产蜜桃| 午夜免费鲁丝| 国产日韩欧美视频二区| 欧美一级a爱片免费观看看| 久久久精品免费免费高清| 日韩成人伦理影院| 久久久久久久久久久久大奶| 一本—道久久a久久精品蜜桃钙片| 精品一区二区免费观看| 狂野欧美激情性bbbbbb| 精品一区二区三卡| 在线观看av片永久免费下载| 亚洲欧洲精品一区二区精品久久久 | 国产黄色视频一区二区在线观看| 成年av动漫网址| 九草在线视频观看| 国产伦精品一区二区三区视频9| 春色校园在线视频观看| 欧美变态另类bdsm刘玥| 免费观看a级毛片全部| av一本久久久久| 成人二区视频| 99视频精品全部免费 在线| 亚洲av福利一区| 久久精品夜色国产| 99国产精品免费福利视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av免费高清在线观看| 日韩欧美一区视频在线观看 | 国产极品粉嫩免费观看在线 | 十八禁高潮呻吟视频 | 观看美女的网站| 中文字幕免费在线视频6| 高清午夜精品一区二区三区| 高清黄色对白视频在线免费看 | 在线观看www视频免费| 午夜福利,免费看| 最新中文字幕久久久久| 99久国产av精品国产电影| 精品人妻一区二区三区麻豆| 看免费成人av毛片| 99久久中文字幕三级久久日本| 熟女人妻精品中文字幕| 色5月婷婷丁香| 亚洲怡红院男人天堂| 国语对白做爰xxxⅹ性视频网站| 午夜免费观看性视频| 成人黄色视频免费在线看| 性色avwww在线观看| 人妻 亚洲 视频| av免费观看日本| 中文字幕精品免费在线观看视频 | 夜夜爽夜夜爽视频| 黄色一级大片看看| av福利片在线| 女人精品久久久久毛片| 下体分泌物呈黄色| 国产精品不卡视频一区二区| 汤姆久久久久久久影院中文字幕| 亚洲综合精品二区| 国产精品久久久久久久电影| 久久午夜福利片| 少妇高潮的动态图| 黄色怎么调成土黄色| 午夜免费男女啪啪视频观看| 成人特级av手机在线观看| 欧美区成人在线视频| 欧美性感艳星| 另类亚洲欧美激情| 国产精品免费大片| 秋霞伦理黄片| 免费不卡的大黄色大毛片视频在线观看| 免费黄网站久久成人精品| 人人妻人人澡人人看| 色视频www国产| 国产精品国产三级专区第一集| 亚洲精品国产av蜜桃| 国产熟女午夜一区二区三区 | 我的女老师完整版在线观看| av女优亚洲男人天堂| 欧美三级亚洲精品| 国产在视频线精品| 日本黄色日本黄色录像| 91在线精品国自产拍蜜月| 亚洲在久久综合| av又黄又爽大尺度在线免费看| videos熟女内射| 七月丁香在线播放| 免费黄频网站在线观看国产| 18禁裸乳无遮挡动漫免费视频| 内地一区二区视频在线| 街头女战士在线观看网站| 中文字幕精品免费在线观看视频 | 欧美激情国产日韩精品一区| 亚州av有码| 国产精品久久久久久精品古装| 男女免费视频国产| 天堂8中文在线网| 国产在线免费精品| 亚洲国产最新在线播放| 久久女婷五月综合色啪小说| 乱系列少妇在线播放| 天天躁夜夜躁狠狠久久av| 久久精品夜色国产| 国产极品粉嫩免费观看在线 | 国产免费视频播放在线视频| www.av在线官网国产| 丰满少妇做爰视频| 黄色日韩在线| 精品少妇内射三级| .国产精品久久| 亚洲欧美清纯卡通| 免费观看av网站的网址| 性色avwww在线观看| 制服丝袜香蕉在线| 国产精品久久久久久精品电影小说| 最新的欧美精品一区二区| 建设人人有责人人尽责人人享有的| 国产黄频视频在线观看| 国产淫片久久久久久久久| 亚洲美女黄色视频免费看| 在线观看国产h片| 国产精品99久久久久久久久| 国产欧美日韩综合在线一区二区 | 99热这里只有精品一区| 女性生殖器流出的白浆| 国产欧美另类精品又又久久亚洲欧美| 熟女人妻精品中文字幕| 久久久精品免费免费高清| 亚洲精品自拍成人| 五月玫瑰六月丁香| 国产淫语在线视频| 午夜久久久在线观看| 国产淫语在线视频| 最近最新中文字幕免费大全7| 久久久久久人妻| 亚洲av.av天堂| 久久精品国产a三级三级三级| 黑人猛操日本美女一级片| 久久婷婷青草| 男女国产视频网站| 五月开心婷婷网| 午夜福利视频精品| 国产高清不卡午夜福利| 少妇猛男粗大的猛烈进出视频| 纵有疾风起免费观看全集完整版| 久久久久国产精品人妻一区二区| 亚洲精品乱码久久久久久按摩| 国产精品99久久99久久久不卡 | 大香蕉97超碰在线| 精品熟女少妇av免费看| 国产精品久久久久久精品电影小说| 男女边摸边吃奶| 日本黄色片子视频| 日韩中字成人| 国产乱人偷精品视频| 国产精品人妻久久久影院| 一级毛片aaaaaa免费看小| 噜噜噜噜噜久久久久久91| 99国产精品免费福利视频| 六月丁香七月| 天天操日日干夜夜撸| 老司机影院成人| 成人毛片a级毛片在线播放| 男人爽女人下面视频在线观看| 亚洲av日韩在线播放| 久久久久国产精品人妻一区二区| 国产精品偷伦视频观看了| 如何舔出高潮| 美女国产视频在线观看| 日日啪夜夜爽| 黑人高潮一二区| 日韩精品免费视频一区二区三区 | 91在线精品国自产拍蜜月| 亚洲丝袜综合中文字幕| 国产成人免费观看mmmm| 建设人人有责人人尽责人人享有的| 在线看a的网站| 六月丁香七月| 99久久精品一区二区三区| 少妇人妻久久综合中文| 久久精品久久久久久久性| 亚洲国产精品一区二区三区在线| 啦啦啦中文免费视频观看日本| 成年女人在线观看亚洲视频| 日韩制服骚丝袜av| 建设人人有责人人尽责人人享有的| 亚洲av日韩在线播放| 亚洲色图综合在线观看| a级片在线免费高清观看视频| av又黄又爽大尺度在线免费看| 欧美成人午夜免费资源| 这个男人来自地球电影免费观看 | 中文字幕免费在线视频6| 最新中文字幕久久久久| 国产成人精品一,二区| 一级毛片电影观看| 欧美人与善性xxx| 免费观看无遮挡的男女| 美女脱内裤让男人舔精品视频| 国产精品一区www在线观看| 三级国产精品欧美在线观看| av卡一久久| 黄色日韩在线| av女优亚洲男人天堂| 街头女战士在线观看网站| 成年av动漫网址| 亚洲精品日韩av片在线观看| 丁香六月天网| 欧美国产精品一级二级三级 | 一边亲一边摸免费视频| 中文资源天堂在线| 好男人视频免费观看在线| 色94色欧美一区二区| 卡戴珊不雅视频在线播放| 熟女av电影| 亚洲性久久影院|