• <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ī)律研究
    国产激情偷乱视频一区二区| 黄色片一级片一级黄色片| 亚洲精品av麻豆狂野| 午夜福利成人在线免费观看| av片东京热男人的天堂| 亚洲成a人片在线一区二区| 亚洲欧美激情综合另类| 久久久久久亚洲精品国产蜜桃av| 国产精品亚洲av一区麻豆| 免费看日本二区| 精品欧美一区二区三区在线| x7x7x7水蜜桃| 精品午夜福利视频在线观看一区| 国产精品一区二区免费欧美| 免费在线观看视频国产中文字幕亚洲| 此物有八面人人有两片| 久久久国产精品麻豆| 欧美色视频一区免费| 正在播放国产对白刺激| 日韩中文字幕欧美一区二区| 亚洲成av片中文字幕在线观看| 欧美成人一区二区免费高清观看 | a在线观看视频网站| 亚洲欧美精品综合久久99| 亚洲国产日韩欧美精品在线观看 | 老司机福利观看| 韩国av一区二区三区四区| 日本免费一区二区三区高清不卡| 欧美日韩中文字幕国产精品一区二区三区| 午夜激情福利司机影院| 十八禁人妻一区二区| 黄片大片在线免费观看| 亚洲精华国产精华精| 亚洲成人免费电影在线观看| 手机成人av网站| 亚洲精品一卡2卡三卡4卡5卡| 不卡一级毛片| 久久久久久人人人人人| 亚洲精品在线观看二区| 欧美日韩黄片免| 久久久久久九九精品二区国产 | 色哟哟哟哟哟哟| 2021天堂中文幕一二区在线观 | 99热这里只有精品一区 | 午夜激情av网站| 午夜视频精品福利| 亚洲国产欧洲综合997久久, | 亚洲精品在线观看二区| 亚洲性夜色夜夜综合| 国产麻豆成人av免费视频| 国产99久久九九免费精品| 国产一级毛片七仙女欲春2 | 在线十欧美十亚洲十日本专区| 午夜福利18| 亚洲狠狠婷婷综合久久图片| 亚洲成国产人片在线观看| 欧美绝顶高潮抽搐喷水| 国产男靠女视频免费网站| a级毛片a级免费在线| 三级毛片av免费| 色播在线永久视频| 国产成人系列免费观看| 69av精品久久久久久| 特大巨黑吊av在线直播 | 亚洲成av人片免费观看| 中文字幕av电影在线播放| 久久精品人妻少妇| 少妇裸体淫交视频免费看高清 | 91大片在线观看| 亚洲五月色婷婷综合| 免费av毛片视频| 丝袜人妻中文字幕| 中文字幕最新亚洲高清| 老熟妇仑乱视频hdxx| 国产黄a三级三级三级人| 18禁黄网站禁片午夜丰满| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看| 亚洲av电影在线进入| 国内精品久久久久久久电影| www日本在线高清视频| 国产真人三级小视频在线观看| 久久精品国产亚洲av高清一级| 校园春色视频在线观看| 妹子高潮喷水视频| 免费高清视频大片| 亚洲欧美激情综合另类| 欧美激情久久久久久爽电影| 日本 av在线| 亚洲国产欧美一区二区综合| 久久人妻av系列| 精品国产亚洲在线| 女性生殖器流出的白浆| 国产成人av教育| 国产熟女xx| 91成人精品电影| 欧美性猛交╳xxx乱大交人| 国产亚洲欧美98| 国产精品亚洲美女久久久| 久久国产乱子伦精品免费另类| 亚洲精品美女久久av网站| 欧美黄色片欧美黄色片| 精品第一国产精品| 精品免费久久久久久久清纯| 欧美日韩中文字幕国产精品一区二区三区| 少妇裸体淫交视频免费看高清 | 久久亚洲精品不卡| 老熟妇仑乱视频hdxx| 亚洲av成人av| 成年免费大片在线观看| 欧美日韩精品网址| 国产精品,欧美在线| 香蕉国产在线看| 夜夜躁狠狠躁天天躁| 一区二区三区高清视频在线| 黄色视频,在线免费观看| 久久草成人影院| 激情在线观看视频在线高清| 99在线人妻在线中文字幕| 日韩一卡2卡3卡4卡2021年| 91在线观看av| 91麻豆精品激情在线观看国产| 亚洲五月婷婷丁香| 亚洲午夜理论影院| 美女大奶头视频| 免费人成视频x8x8入口观看| 国产片内射在线| 999久久久精品免费观看国产| 琪琪午夜伦伦电影理论片6080| 香蕉国产在线看| 99热6这里只有精品| 久久久国产欧美日韩av| av中文乱码字幕在线| 久久国产精品男人的天堂亚洲| 免费av毛片视频| 精品国产乱码久久久久久男人| 久久香蕉激情| e午夜精品久久久久久久| 国产又黄又爽又无遮挡在线| 久久久久国产精品人妻aⅴ院| 侵犯人妻中文字幕一二三四区| 天天躁夜夜躁狠狠躁躁| 亚洲国产精品sss在线观看| 欧美激情极品国产一区二区三区| 一边摸一边抽搐一进一小说| 桃红色精品国产亚洲av| 热99re8久久精品国产| 欧美另类亚洲清纯唯美| 国产亚洲精品综合一区在线观看 | 天天一区二区日本电影三级| 神马国产精品三级电影在线观看 | 国产一区二区三区在线臀色熟女| 久久性视频一级片| 精品久久久久久成人av| 欧美乱码精品一区二区三区| 亚洲成人久久性| 国产高清激情床上av| 真人一进一出gif抽搐免费| 久久 成人 亚洲| 亚洲国产毛片av蜜桃av| 午夜精品在线福利| 又黄又爽又免费观看的视频| 亚洲成人国产一区在线观看| 99久久精品国产亚洲精品| 亚洲av五月六月丁香网| 亚洲精品美女久久av网站| 欧美中文日本在线观看视频| 麻豆成人av在线观看| 长腿黑丝高跟| 欧美另类亚洲清纯唯美| 少妇粗大呻吟视频| 久久久久久国产a免费观看| www日本黄色视频网| 国产乱人伦免费视频| 欧美成人午夜精品| 好看av亚洲va欧美ⅴa在| 超碰成人久久| 国产精品久久久久久亚洲av鲁大| 欧美激情极品国产一区二区三区| 日本免费a在线| 日本免费一区二区三区高清不卡| 波多野结衣高清无吗| 岛国视频午夜一区免费看| 国产成人精品无人区| 两人在一起打扑克的视频| 真人做人爱边吃奶动态| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品一卡2卡三卡4卡5卡| 国产野战对白在线观看| 黄色视频不卡| 亚洲色图 男人天堂 中文字幕| 人妻久久中文字幕网| 超碰成人久久| 精品国产一区二区三区四区第35| 18禁黄网站禁片免费观看直播| 免费一级毛片在线播放高清视频| 50天的宝宝边吃奶边哭怎么回事| 熟女电影av网| 国产主播在线观看一区二区| 国产黄片美女视频| 精品人妻1区二区| 免费高清在线观看日韩| 香蕉国产在线看| 欧美一级a爱片免费观看看 | 成人永久免费在线观看视频| av中文乱码字幕在线| av在线天堂中文字幕| 免费在线观看黄色视频的| 夜夜躁狠狠躁天天躁| 亚洲精品中文字幕在线视频| 国产精品二区激情视频| 黄片小视频在线播放| ponron亚洲| 日韩欧美一区二区三区在线观看| 日韩欧美免费精品| 桃红色精品国产亚洲av| 国产激情久久老熟女| 一本大道久久a久久精品| av天堂在线播放| av在线播放免费不卡| 中文字幕另类日韩欧美亚洲嫩草| 精品熟女少妇八av免费久了| 精品国内亚洲2022精品成人| 老汉色∧v一级毛片| 在线免费观看的www视频| 午夜福利视频1000在线观看| 欧美久久黑人一区二区| 夜夜看夜夜爽夜夜摸| 国产av一区二区精品久久| 久久精品成人免费网站| 中文字幕最新亚洲高清| 熟女电影av网| 国产av一区二区精品久久| 一边摸一边抽搐一进一小说| 国产又黄又爽又无遮挡在线| 久久国产精品影院| 国产亚洲精品久久久久5区| 欧美久久黑人一区二区| 少妇被粗大的猛进出69影院| 人妻久久中文字幕网| 天天躁夜夜躁狠狠躁躁| 嫁个100分男人电影在线观看| 无限看片的www在线观看| 久久精品影院6| www.精华液| 欧美+亚洲+日韩+国产| 最近最新免费中文字幕在线| 在线观看舔阴道视频| 黄片小视频在线播放| 精品欧美一区二区三区在线| 午夜影院日韩av| 欧美色视频一区免费| 精品福利观看| 性色av乱码一区二区三区2| 啦啦啦韩国在线观看视频| av欧美777| 亚洲成a人片在线一区二区| 久久久久久免费高清国产稀缺| av视频在线观看入口| 国产成人精品久久二区二区免费| 十八禁网站免费在线| 日本在线视频免费播放| av在线播放免费不卡| aaaaa片日本免费| 在线观看www视频免费| 国产高清有码在线观看视频 | 91老司机精品| 精品国产美女av久久久久小说| 两人在一起打扑克的视频| 亚洲精品在线美女| 午夜免费激情av| 欧美性长视频在线观看| 日韩欧美免费精品| 热re99久久国产66热| 国产蜜桃级精品一区二区三区| 丝袜人妻中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 国产精品二区激情视频| 最新在线观看一区二区三区| 精品不卡国产一区二区三区| 亚洲黑人精品在线| 色老头精品视频在线观看| 久久久久久大精品| 久久热在线av| 18禁国产床啪视频网站| 曰老女人黄片| 国产伦人伦偷精品视频| 国产亚洲精品av在线| 99在线人妻在线中文字幕| 天堂影院成人在线观看| 国产精品香港三级国产av潘金莲| 欧美成人免费av一区二区三区| 精品人妻1区二区| 一级作爱视频免费观看| 久久婷婷成人综合色麻豆| 丝袜人妻中文字幕| 日本成人三级电影网站| 在线观看午夜福利视频| 国产黄片美女视频| 国产成人精品久久二区二区91| 色综合婷婷激情| 99在线人妻在线中文字幕| 欧美大码av| 久久精品aⅴ一区二区三区四区| 亚洲成a人片在线一区二区| 欧美一区二区精品小视频在线| 禁无遮挡网站| 亚洲熟女毛片儿| 人成视频在线观看免费观看| 少妇熟女aⅴ在线视频| 一个人免费在线观看的高清视频| 高潮久久久久久久久久久不卡| 久久99热这里只有精品18| 免费看日本二区| 国产男靠女视频免费网站| 国产麻豆成人av免费视频| 色哟哟哟哟哟哟| 十八禁人妻一区二区| 婷婷精品国产亚洲av| 国产黄色小视频在线观看| 国产成人精品久久二区二区免费| 淫秽高清视频在线观看| 一级毛片高清免费大全| 90打野战视频偷拍视频| 国产成人精品久久二区二区免费| 亚洲中文av在线| 亚洲在线自拍视频| 亚洲av美国av| 欧美久久黑人一区二区| 波多野结衣高清作品| 久久久久精品国产欧美久久久| 国产真人三级小视频在线观看| 亚洲一区高清亚洲精品| 国产三级在线视频| 久久精品国产99精品国产亚洲性色| 午夜福利在线在线| 亚洲人成伊人成综合网2020| 亚洲精品一区av在线观看| 日本一区二区免费在线视频| 日韩精品青青久久久久久| 免费在线观看黄色视频的| 国产成人欧美| 欧美人与性动交α欧美精品济南到| 妹子高潮喷水视频| 久久久国产成人免费| 欧美亚洲日本最大视频资源| 亚洲一卡2卡3卡4卡5卡精品中文| 悠悠久久av| 后天国语完整版免费观看| 国产一级毛片七仙女欲春2 | 一级a爱片免费观看的视频| 老司机午夜十八禁免费视频| 日本五十路高清| 午夜老司机福利片| 久久欧美精品欧美久久欧美| 很黄的视频免费| 91字幕亚洲| 18禁黄网站禁片午夜丰满| 国内少妇人妻偷人精品xxx网站 | 国产亚洲av嫩草精品影院| 欧美成人午夜精品| 欧美成狂野欧美在线观看| 动漫黄色视频在线观看| 久久香蕉激情| 自线自在国产av| 无人区码免费观看不卡| 欧美又色又爽又黄视频| 两人在一起打扑克的视频| 2021天堂中文幕一二区在线观 | 香蕉av资源在线| 中文字幕最新亚洲高清| 国产成年人精品一区二区| 精品久久久久久久久久免费视频| 日韩三级视频一区二区三区| 国产激情偷乱视频一区二区| 极品教师在线免费播放| 又黄又粗又硬又大视频| 国产高清有码在线观看视频 | 丝袜在线中文字幕| 欧美成狂野欧美在线观看| 色哟哟哟哟哟哟| 久久久国产成人免费| 少妇粗大呻吟视频| 亚洲人成伊人成综合网2020| 天天躁狠狠躁夜夜躁狠狠躁| 欧美成人免费av一区二区三区| 午夜福利欧美成人| 国产久久久一区二区三区| 国产主播在线观看一区二区| 精品第一国产精品| 亚洲精品粉嫩美女一区| 午夜日韩欧美国产| 一夜夜www| 亚洲va日本ⅴa欧美va伊人久久| 最近最新中文字幕大全免费视频| 成年版毛片免费区| 视频在线观看一区二区三区| 精品国产亚洲在线| 国产极品粉嫩免费观看在线| 婷婷精品国产亚洲av在线| 热99re8久久精品国产| 人人妻,人人澡人人爽秒播| 两人在一起打扑克的视频| 国产在线观看jvid| www.熟女人妻精品国产| www.999成人在线观看| 可以在线观看毛片的网站| 在线天堂中文资源库| 可以在线观看毛片的网站| a级毛片a级免费在线| 一本综合久久免费| 日韩av在线大香蕉| 夜夜看夜夜爽夜夜摸| 1024手机看黄色片| 久久久久亚洲av毛片大全| 99久久精品国产亚洲精品| 国产在线观看jvid| 精品久久久久久久久久久久久 | 国内久久婷婷六月综合欲色啪| 午夜免费观看网址| 99精品在免费线老司机午夜| 露出奶头的视频| 女人爽到高潮嗷嗷叫在线视频| 天天添夜夜摸| 曰老女人黄片| 青草久久国产| 欧美一区二区精品小视频在线| 亚洲人成伊人成综合网2020| 精品国产一区二区三区四区第35| 国产又黄又爽又无遮挡在线| 波多野结衣高清作品| 精品高清国产在线一区| а√天堂www在线а√下载| 国产乱人伦免费视频| 一本久久中文字幕| 亚洲全国av大片| 在线观看免费午夜福利视频| 欧美日本亚洲视频在线播放| 国产亚洲欧美精品永久| 国产一卡二卡三卡精品| 国产激情偷乱视频一区二区| 免费女性裸体啪啪无遮挡网站| 香蕉久久夜色| 搡老妇女老女人老熟妇| 亚洲专区国产一区二区| 在线观看日韩欧美| 国产亚洲精品第一综合不卡| 日韩视频一区二区在线观看| 人人妻人人看人人澡| 国产精品av久久久久免费| 国产精品乱码一区二三区的特点| 中国美女看黄片| 国产亚洲精品久久久久久毛片| 久久久国产成人免费| 亚洲真实伦在线观看| 欧美一级毛片孕妇| 午夜免费观看网址| 波多野结衣巨乳人妻| 国产一区二区在线av高清观看| 久久 成人 亚洲| 欧美性猛交╳xxx乱大交人| www日本在线高清视频| www.www免费av| 少妇粗大呻吟视频| 狂野欧美激情性xxxx| 日韩欧美一区二区三区在线观看| 国产激情欧美一区二区| 51午夜福利影视在线观看| 老司机在亚洲福利影院| 久9热在线精品视频| 人人妻,人人澡人人爽秒播| 欧美日韩亚洲国产一区二区在线观看| 88av欧美| 国产亚洲精品一区二区www| 国产成人精品久久二区二区91| 亚洲三区欧美一区| 亚洲人成伊人成综合网2020| 国产主播在线观看一区二区| 精品国内亚洲2022精品成人| 18禁美女被吸乳视频| 熟女电影av网| 在线看三级毛片| 一区二区日韩欧美中文字幕| 中文字幕av电影在线播放| 搡老妇女老女人老熟妇| 动漫黄色视频在线观看| 国产亚洲欧美98| 日韩欧美一区二区三区在线观看| 真人做人爱边吃奶动态| 色播亚洲综合网| 麻豆久久精品国产亚洲av| 日韩一卡2卡3卡4卡2021年| 欧美黑人精品巨大| 啦啦啦观看免费观看视频高清| 日韩视频一区二区在线观看| 少妇 在线观看| bbb黄色大片| av在线播放免费不卡| 日韩大尺度精品在线看网址| 一夜夜www| 欧美三级亚洲精品| 日韩一卡2卡3卡4卡2021年| 久久精品aⅴ一区二区三区四区| 国产99久久九九免费精品| 搡老妇女老女人老熟妇| 欧美日韩黄片免| 首页视频小说图片口味搜索| 亚洲av美国av| 国产成+人综合+亚洲专区| 视频在线观看一区二区三区| 久久久久亚洲av毛片大全| svipshipincom国产片| 国产三级黄色录像| 亚洲国产精品合色在线| 午夜精品久久久久久毛片777| 久久国产乱子伦精品免费另类| 国产又黄又爽又无遮挡在线| av在线天堂中文字幕| 欧美成人午夜精品| 桃色一区二区三区在线观看| 亚洲第一青青草原| 在线观看舔阴道视频| 99riav亚洲国产免费| av中文乱码字幕在线| 精品久久久久久久久久久久久 | 免费在线观看日本一区| 听说在线观看完整版免费高清| 午夜激情av网站| 国产极品粉嫩免费观看在线| 久热这里只有精品99| 成人精品一区二区免费| 精品久久久久久久末码| 人人妻,人人澡人人爽秒播| 国产片内射在线| 亚洲人成77777在线视频| 麻豆成人av在线观看| 在线看三级毛片| 国产私拍福利视频在线观看| 9191精品国产免费久久| 高清毛片免费观看视频网站| 午夜影院日韩av| 日本一本二区三区精品| 18禁美女被吸乳视频| 老汉色∧v一级毛片| 老司机福利观看| 黄色视频不卡| av免费在线观看网站| 一进一出抽搐动态| 日本 av在线| 亚洲中文日韩欧美视频| 午夜福利在线观看吧| av中文乱码字幕在线| 日本三级黄在线观看| 久久久久久久午夜电影| 国产精品99久久99久久久不卡| 非洲黑人性xxxx精品又粗又长| 日韩欧美一区视频在线观看| 亚洲av中文字字幕乱码综合 | 成人av一区二区三区在线看| 老司机福利观看| 国产欧美日韩精品亚洲av| 欧美国产精品va在线观看不卡| 别揉我奶头~嗯~啊~动态视频| 国产真人三级小视频在线观看| 一级毛片精品| 欧美zozozo另类| 伊人久久大香线蕉亚洲五| 黄色丝袜av网址大全| 久久中文字幕一级| 999精品在线视频| 美女扒开内裤让男人捅视频| 国产精品爽爽va在线观看网站 | 最近最新中文字幕大全免费视频| 午夜福利高清视频| 伦理电影免费视频| 亚洲性夜色夜夜综合| 久久久久久免费高清国产稀缺| 老司机福利观看| 亚洲无线在线观看| 夜夜看夜夜爽夜夜摸| 婷婷丁香在线五月| 国产精品日韩av在线免费观看| 麻豆成人av在线观看| 女人被狂操c到高潮| 国产黄a三级三级三级人| 国产精品免费视频内射| 欧美黄色淫秽网站| 亚洲久久久国产精品| 精品福利观看| 国产区一区二久久| 天堂动漫精品| 亚洲一区中文字幕在线| 国产成人av教育| 日韩欧美三级三区| 国产精品国产高清国产av| 日本一区二区免费在线视频| 国产私拍福利视频在线观看| 国产精品乱码一区二三区的特点| 女警被强在线播放| 亚洲第一青青草原| 亚洲一区中文字幕在线| 国产日本99.免费观看| 免费在线观看完整版高清| 亚洲午夜精品一区,二区,三区| 人人妻人人澡人人看| 久久精品国产清高在天天线| 日本熟妇午夜| 一进一出抽搐动态| 99热只有精品国产| 亚洲av五月六月丁香网| 两人在一起打扑克的视频| 亚洲精品美女久久久久99蜜臀| 白带黄色成豆腐渣| aaaaa片日本免费| 亚洲精品美女久久av网站|