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

    凡納濱對蝦微衛(wèi)星分子標記的開發(fā)及不同養(yǎng)殖家系遺傳多態(tài)性分析

    2020-04-14 04:59:03李強勇李旻曾地剛朱威霖彭敏楊春玲劉青云趙永貞陳秀荔陳曉漢
    南方農業(yè)學報 2020年2期
    關鍵詞:測序

    李強勇 李旻 曾地剛 朱威霖 彭敏 楊春玲 劉青云 趙永貞 陳秀荔 陳曉漢

    摘要:【目的】基于第三代高通量測序技術——單分子實時(SMRT)測序開發(fā)凡納濱對蝦微衛(wèi)星分子標記,為凡納濱對蝦的遺傳育種研究提供基礎資料。【方法】采用SMRT測序對凡納濱對蝦轉錄組進行測序,經(jīng)Illumina測序糾錯后利用MISA對獲得的轉錄組數(shù)據(jù)進行分析;以Primer 3.0設計微衛(wèi)星引物,隨機挑選40對微衛(wèi)星引物進行PCR擴增驗證,并選用擴增成功且具有多態(tài)性的微衛(wèi)星引物用于不同養(yǎng)殖家系凡納濱對蝦遺傳多態(tài)性分析?!窘Y果】凡納濱對蝦轉錄組SMRT測序共獲得51367條非冗余全長轉錄本序列,鑒定出39674條微衛(wèi)星序列;微衛(wèi)星的分布密度為0.232 SSR/kb;以二核苷酸重復微衛(wèi)星序列分布最多,共有22223條(占56.01%)。隨機選取的40對微衛(wèi)星引物中有26對微衛(wèi)星引物能成功擴增出特異性條帶;微衛(wèi)星熒光分型結果顯示,有16個微衛(wèi)星位點具有多態(tài)性,共獲得67個等位基因,平均等位基因數(shù)(Na)為4.2個,計算獲得的平均觀測雜合度(Ho)為0.511,平均期望雜合度(He)為0.451,平均多態(tài)性信息含量(PIC)為0.489。在16個微衛(wèi)星位點中,有2個微衛(wèi)星位點為低度多態(tài)性(PIC<0.25),6個微衛(wèi)星位點為中度多態(tài)性(0.250.50);有4個微衛(wèi)星位點偏移Hardy-Weinberg平衡(P<0.05),其余12個微衛(wèi)星位點均未偏離Hardy-Weinberg平衡(P>0.05)。【結論】采用SMRT測序開發(fā)凡納濱對蝦微衛(wèi)星分子標記是一種簡單而高效的途徑,有助于開展凡納濱對蝦的遺傳育種研究工作。

    關鍵詞: 凡納濱對蝦;單分子實時(SMRT)測序;微衛(wèi)星分子標記;多態(tài)性分析

    中圖分類號: S917.4? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻標志碼: A 文章編號:2095-1191(2020)02-0429-08

    Development of microsatellite markers of Litopenaeus vannamei and genetic polymorphism analysis of different cultured families

    LI Qiang-yong, LI Min, ZENG Di-gang, ZHU Wei-lin, PENG Min, YANG Chun-ling, LIU Qing-yun, ZHAO Yong-zhen, CHEN Xiu-li, CHEN Xiao-han*

    (Guangxi Academy of Fishery Sciences/Guangxi Key Laboratory of Aquatic Genetic Breeding and Healthy

    Aquaculture, Nanning? 530021, China)

    Abstract:【Objective】Microsatellite markers of Litopenaeus vannamei based on the third-generation high-throughput sequencing technology, single-molecule real-time (SMRT) sequencing were developed in order to provide basic data for genetic and breeding research of L. vannamei. 【Method】The transcriptome of L. vannamei was sequenced using SMRT. The obtained transcriptome data were corrected by Illumina sequencing, then analyzed by MISA software. The microsatellite primers were designed with Primer 3.0 software, and 40 pairs of microsatellite primers were randomly selected for PCR amplification.? Microsatellite primers with successful amplification and polymorphism were selected for genetic polymorphism analysis of L. vannamei in different breeding families. 【Result】In total, 51367 non-redundant full length transcript sequences and 39674 SSR loci were obtained? by SMRT sequencing of the L. vannamei transcriptome. The distribution density of microsatellites was 0.232 SSR/kb; the dinucleotide repeat microsatellite sequences were the most distribu-ted, with a total of 22223(56.01%). Of the 40 randomly selected microsatellite primers, 26 pairs of microsatellite primers successfully amplified specific bands. The microsatellite fluorescence phenotyping showed that 16 microsatellite loci were polymorphic. A total of 67 alleles were obtained, and the average number of alleles(Na) was 4.2. The calculated average observed heterozygosity(Ho) was 0.511, the average expected heterozygosity(He) was 0.451,and the average polymorphism information content(PIC) was 0.489. Among the 16 microsatellite loci, 2 loci were low polymorphisms (PIC<0.25), 6 loci were moderate polymorphisms(0.25 0.50);4 microsatellite loci deviated from Hardy-Weinberg equilibrium(P<0.05), and the remaining 12 loci did not deviate from Hardy-Weinberg equilibrium(P>0.05). 【Conclusion】Using SMRT sequencing to develop microsatellite markers of L. vannamei is a simple and efficient way, which is helpful for the genetic and breeding research of L. vannamei.

    Key words: Litopenaeus vannamei; single molecule real-time sequencing; microsatellite markers; polymorphism analysis

    Foundation item: Guangxi Innovation Driven Development Project(Guike AA17204088-1); National Fishery Science Data Center Observation and Detection Project of Agricultural Basic Long-term Scientific and Technological Work(ZX08S211664); Science and Technology Plan Project of Guangxi Aquatic Animal Husbandry and Veterinary Bureau (Guiyumuke 201633015); Guangxi Key Laboratory of Genetic Breeding and Healthy Breeding Opening Project (17-A-04-04)

    0 引言

    【研究意義】凡納濱對蝦(Litopenaeus vannamei)又稱為南美白對蝦,是世界上養(yǎng)殖最廣泛的甲殼類動物(陳錨等,2008)。凡納濱對蝦具有耐高溫、抗病力強、出肉率高、營養(yǎng)豐富、生長速度快及鹽度適應范圍廣等特點,現(xiàn)已發(fā)展成為對蝦養(yǎng)殖的重要對象(王興強等,2004;黃薇等,2014)。1988年,中國科學院海洋研究所首次引進凡納濱對蝦,并于1992年突破人工繁殖難關,之后凡納濱對蝦在全國范圍內廣泛養(yǎng)殖,且各地不斷從國外引進親蝦用于蝦苗生產(chǎn)(王興強等,2004),如今凡納濱對蝦已成為我國對蝦養(yǎng)殖的主導品種。但近年來凡納濱對蝦病害發(fā)生頻繁,包括白斑綜合征、桃拉綜合征、傳染性皮下及造血組織壞死病和弧菌病等,給其養(yǎng)殖業(yè)造成巨大經(jīng)濟損失(Lightner,2011)。我國養(yǎng)殖的凡納濱對蝦主要是從國外引進的品種或品系,且多為雜交后代,遺傳背景較復雜,繁育多代后其性狀出現(xiàn)分離,表現(xiàn)為群體生長速度減慢、抗病和抗逆能力減弱等(熊建華等,2011)。因此,急需研發(fā)出有效的分子標記技術用于凡納濱對蝦種質資源鑒定及分子標記選育等研究工作?!厩叭搜芯窟M展】微衛(wèi)星(Microsatellite)又稱簡單重復序列(Simple sequence repeats,SSR),由2~6個核苷酸的串聯(lián)重復片段構成,是目前應用較廣泛的分子標記之一(盛巖等,2002;周大顏等,2019)。因其穩(wěn)定性和可重復性高,且具有特異性和共顯性的特點,現(xiàn)已廣泛應用于動植物的生物遺傳多樣性分析、種質資源鑒定及分子標記等研究領域(孫效文等,2008;張麗穎等,2009;凌士鵬等,2018)。凡納濱對蝦微衛(wèi)星分子標記在1996年已有報道(Garcia et al.,1996);而Meehan等(2003)最早實現(xiàn)批量開發(fā)微衛(wèi)星分子標記,通過構建基因組DNA文庫和使用探針篩選的方法共得包含573個微衛(wèi)星的173條序列。隨著第二代高通量測序技術的發(fā)展,已有學者采用高通量測序技術開發(fā)凡納濱對蝦微衛(wèi)星分子標記的研究報道。馬寧和曾地剛(2013)運用454高通量測序技術篩選出588條凡納濱對蝦微衛(wèi)星序列;Yu等(2014)利用Illumina HiSeq 2000測序技術查找凡納濱對蝦SNP位點;楊銘等(2017)采用Illumina測序技術進行凡納濱對蝦轉錄組測序,結果鑒定獲得14767條微衛(wèi)星序列;Lu等(2018)通過對凡納濱對蝦進行轉錄組測序,篩選出49個與氨相關的SNPs位點,且進一步鑒定發(fā)現(xiàn)有12個SNPs位點與氨耐受性相關;Santos等(2018)通過RNA測序技術(RNA-seq)分析鑒定了與凡納濱對蝦免疫應答及生長性能相關的SNP位點?!颈狙芯壳腥朦c】單分子實時(Single-molecule real-time,SMRT)測序是近年來由Pacific Biosciences公司開發(fā)的第三代高通量測序技術,為實現(xiàn)全長單分子DNA序列的直接測序提供了平臺(Eid et al.,2009)。SMRT測序長度可達20 kb,能對完整的cDNA分子進行測序,而無需片段化或測序后組裝(Roberts et al.,2013),但至今未見將SMRT測序應用于凡納濱對蝦微衛(wèi)星分子標記開發(fā)的研究報道?!緮M解決的關鍵問題】利用SMRT測序對凡納濱對蝦全長轉錄組進行測序分析,從獲得的全長RNA序列中篩選微衛(wèi)星序列并設計微衛(wèi)星引物,隨機挑選微衛(wèi)星引物對不同養(yǎng)殖家系凡納濱對蝦進行遺傳多態(tài)性分析,為凡納濱對蝦的遺傳育種研究提供基礎資料。

    1 材料與方法

    1. 1 試驗材料

    供試桂海1號凡納濱對蝦由廣西水產(chǎn)科學研究院南美白對蝦遺傳育種中心提供。取6尾凡納濱對蝦的肌肉、心臟和肝胰腺組織,混合后置于液氮中保存?zhèn)溆?。從無特定病原的桂海1號凡納濱對蝦種質庫中隨機取16個不同養(yǎng)殖家系,對蝦平均體長16.0 cm,平均體質量26.8 g/尾,取肌肉組織置于-80 ℃冰箱中保存?zhèn)溆谩?/p>

    1. 2 總RNA提取

    按TRIzol LS Reagent試劑盒(Invitrogen,美國)說明從混合的凡納濱對蝦組織中提取總RNA,使用DNase I(Invitrogen,美國)除去基因組DNA。采用NanoDrop 2000(Thermo Scientific,美國)測量總RNA的純度(OD260/280)、濃度和吸收峰,并以Agilent 2100 Bioanalyser(美國)測定總RNA質量,RIN(RNA integrity number)>7的總RNA用于構建SMRT測序的cDNA文庫。

    1. 3 SMART測序

    使用SMARTer PCR cDNA Synthesis Kit(TaKaRa,日本)將10 μg總RNA反轉錄合成cDNA,以BluePippinTM系統(tǒng)(Sage Science,美國)構建5個大小分別為<1 kb、1~2 kb、2~3 kb、3~6 kb和>6 kb的cDNA文庫。按照PacBio的步驟每個cDNA文庫在PacBio RSII平臺上使用C4試劑進行SMART測序;同時使用TruSeq RNA試劑盒(Illumina,美國)構建Illumina凡納濱對蝦cDNA文庫,在HiSeq2500高通量測序儀上進行測序,測序結果用于SMART測序數(shù)據(jù)的錯誤糾正。

    1. 4 微衛(wèi)星序列鑒定及引物設計

    對SMRT測序的原始數(shù)據(jù)進行質控及Illumina測序糾錯后獲得高質量的全長轉錄本序列,利用MISA(http://pgrc.ipk-gatersleben.de/misa/)搜尋SSR重復基元的區(qū)間。MISA參數(shù)設置如下:2-6、3-5、4-5、5-5、6-5,復合微衛(wèi)星分子標記的最大長度100 bp;然后采用Primer 3.0設計微衛(wèi)星序列的引物。

    1. 5 微衛(wèi)星引物驗證

    從批量設計的微衛(wèi)星引物中隨機挑選40對(表1),委托生工生物工程(上海)股份有限公司合成。不同養(yǎng)殖家系對蝦肌肉組織分別采用動物基因組提取試劑盒[天根生化科技(北京)有限公司]提取DNA,DNA經(jīng)1.0%瓊脂糖凝膠電泳檢查無降解后使用紫外分光光度計測其濃度,以TE(pH 8.0)稀釋至50 ng/μL,-20 ℃保存?zhèn)溆谩?/p>

    PCR反應體系50.0 μL:基因組DNA模板20 ng,2×Taq MasterMix 25.0 μL ,正、反向引物(10 μmol/L)各2.0 μL,超純水補足至50.0 μL。擴增程序:94 ℃預變性1 min;94 ℃ 45 s,58 ℃ 45 s,72 ℃ 45 s,進行35個循環(huán);72 ℃延伸5 min。以2.0%瓊脂糖凝膠在1×TAE緩沖液中150 V下電泳30 min,凝膠成像系統(tǒng)拍照。選取PCR擴增成功的微衛(wèi)星引物在其正向引物5'端標記FAM熒光基團,并對16尾凡納濱對蝦DNA進行PCR擴增。PCR產(chǎn)物經(jīng)2.0%瓊脂糖凝膠電泳檢測后送至上海邁浦生物科技有限公司,在ABI 3730測序儀上進行熒光分型。

    1. 6 微衛(wèi)星數(shù)據(jù)分析

    獲得不同養(yǎng)殖家系凡納濱對蝦的微衛(wèi)星分型數(shù)據(jù)后,使用PopGen32計算每個微衛(wèi)星位點的等位基因數(shù)(Na)、觀測雜合度(Ho)、期望雜合度(He)及多態(tài)信息含量(PIC)。

    2 結果與分析

    2. 1 凡納濱對蝦微衛(wèi)星數(shù)量及其分布密度

    凡納濱對蝦轉錄組SMRT測序共獲得51367條非冗余全長轉錄本序列(GenBank登錄號SRX3267788~ SRX3267801),利用MISA對這些轉錄本序列進行分析,共獲得39674條微衛(wèi)星序列,分布在29461條轉錄本序列中,包含微衛(wèi)星的序列占總序列的57.35%。序列拼接總長度為171322994 bp,微衛(wèi)星的分布密度為0.232 SSR/kb,平均每4.32 kb出現(xiàn)1條微衛(wèi)星序列。各核苷酸重復類型的凡納濱對蝦微衛(wèi)星數(shù)量差異明顯,其分布密度也存在明顯差異。在所鑒定的39674條凡納濱對蝦微衛(wèi)星序列(表2)中,二核苷酸重復微衛(wèi)星序列分布最多,共有22223條(占56.01%);三核苷酸重復微衛(wèi)星序列14781條(占37.26%),四核苷酸重復微衛(wèi)星序列1947條(占4.91%),六核苷酸重復微衛(wèi)星序列441條(占1.11%);分布最少的是五核苷酸重復微衛(wèi)星序列,僅有282條(占0.71%)。

    2. 2 凡納濱對蝦微衛(wèi)星不同重復類型中各核苷酸重復單元的數(shù)量與頻率

    凡納濱對蝦微衛(wèi)星二核苷酸重復類型中不同重復單元的數(shù)量分布情況如圖1所示,以含AT/TA重復單元的微衛(wèi)星序列數(shù)量最多,其次是CA/TG重復單元和GA/TC重復單元,AC/GT重復單元和AG/CT重復單元的數(shù)量差異不明顯。凡納濱對蝦微衛(wèi)星三核苷酸重復類型中不同重復單元的數(shù)量分布情況如圖2所示,以含CTC/GAG重復單元的微衛(wèi)星序列數(shù)量最多,其次是AAT、TTA/TAA、CTG/CAG和AGC/GCT等重復單元。

    2. 3 微衛(wèi)星分子標記驗證結果

    根據(jù)鑒定出的39674條微衛(wèi)星序列,利用Pri-mer 3.0設計39367對微衛(wèi)星引物。為驗證設計的微衛(wèi)星引物能否通過PCR擴增出微衛(wèi)星序列,隨機選取40對微衛(wèi)星引物,在桂海1號凡納濱對蝦DNA混池中進行PCR擴增,結果發(fā)現(xiàn)有26對微衛(wèi)星引物能成功擴增出特異性條帶(圖3),占所選PCR擴增引物數(shù)量的65.00%。為進一步驗證微衛(wèi)星的多態(tài)性,以篩選出的26對微衛(wèi)星引物分別對16個不同養(yǎng)殖家系凡納濱對蝦個體DNA進行PCR擴增,PCR擴增產(chǎn)物在ABI 3730測序儀上進行熒光分型(圖4),并采用GeneMarker 3.0.0分析分型結果。經(jīng)微衛(wèi)星分型后有16個微衛(wèi)星位點顯示出多態(tài)性(表3),多態(tài)性位點的等位基因數(shù)在2~8,共獲得67個等位基因,平均Na為4.2個,計算獲得的平均Ho為0.511,平均He為0.451,平均PIC為0.489。依據(jù)Botstein等(1980)的判斷標準:PIC<0.25代表該位點具有低度多態(tài)性,0.250.50代表該位點具有高度多態(tài)性。在16個微衛(wèi)星位點中,有2個微衛(wèi)星位點為低度多態(tài)性,6個微衛(wèi)星位點為中度多態(tài)性,其余8個微衛(wèi)星位點為高度多態(tài)性。此外,Hardy-Weinberg平衡檢驗結果顯示有4個微衛(wèi)星位點偏移Hardy-Weinberg平衡外(P<0.05),其余12個微衛(wèi)星位點均未偏離Hardy-Weinberg平衡(P>0.05)。

    3 討論

    微衛(wèi)星分子標記對開展水產(chǎn)動物遺傳學研究具有重要意義。為開發(fā)凡納濱對蝦微衛(wèi)星分子標記,本研究采用第三代高通量測序技術——SMRT測序對凡納濱對蝦轉錄組進行測序,結果在獲得51367條非冗余全長轉錄本序列中鑒定出39674條微衛(wèi)星序列,并設計39367對微衛(wèi)星引物。Meehan等(2003)通過傳統(tǒng)的構建文庫和探針篩選僅發(fā)現(xiàn)573個凡納濱對蝦微衛(wèi)星,楊銘等(2017)采用高通量測序技術則鑒定獲得14767條凡納濱對蝦微衛(wèi)星序列。可見,利用高通量測序技術開發(fā)微衛(wèi)星分子標記的優(yōu)勢在于能大規(guī)模發(fā)現(xiàn)微衛(wèi)星序列,極大豐富凡納濱對蝦微衛(wèi)星分子標記的數(shù)量。本研究在分析凡納濱對蝦轉錄組微衛(wèi)星數(shù)量及其分布密度時發(fā)現(xiàn),不同類型重復單元的微衛(wèi)星數(shù)量差異較明顯,其分布密度也存在明顯差異,其中,二核苷酸重復微衛(wèi)星序列分布最多,與馬寧和曾地剛(2013)、楊銘等(2017)采用高通量測序得到的凡納濱對蝦轉錄組微衛(wèi)星結果一致。微衛(wèi)星的分布密度為0.232 SSR/kb,高于楊銘等(2017)采用第二代高通量測序技術獲得的凡納濱對蝦轉錄組微衛(wèi)星分布密度(0.215 SSR/kb),究其原因可能是本研究測序獲得的RNA序列較長,雖然獲得的非冗余全長轉錄本序列數(shù)量(51367條)略少于楊銘等(2017)采用第二代高通量測序技術獲得的非冗余全長轉錄本序列數(shù)量(66815條),但相對較長的RNA序列可鑒定出更多微衛(wèi)星位點。說明采用第三代測序技術開發(fā)微衛(wèi)星分子標記更具優(yōu)勢。

    為驗證設計的微衛(wèi)星引物,隨機挑選40對微衛(wèi)星引物對桂海1號凡納濱對蝦DNA混池進行PCR擴增,結果發(fā)現(xiàn)有26對微衛(wèi)星引物擴增成功(占65.00%)。進一步的微衛(wèi)星分型結果顯示,在26對成功擴增的微衛(wèi)星引物中有16對在不同養(yǎng)殖家系凡納濱對蝦個體DNA的擴增產(chǎn)物中具有多態(tài)性,表明本研究開發(fā)獲得的這些微衛(wèi)星分子標記有相當高的比例可用于凡納濱對蝦遺傳學研究。Na、PIC及遺傳雜合度(H)是檢測種群遺傳多樣性的主要指標。本研究挑選的16個微衛(wèi)星位點的Na在2~9個,平均為4.2個,符合用于遺傳多樣性評估的微衛(wèi)星至少具有4個等位基因的標準(Wajid et al.,2014)。Botstein等(1980)認為:PIC<0.25代表該位點具有低度多態(tài)性,0.250.50代表該位點具有高度多態(tài)性。本研究發(fā)現(xiàn)在16個微衛(wèi)星位點中,有2個微衛(wèi)星位點為低度多態(tài)性,6個微衛(wèi)星位點為中度多態(tài)性,其余8個微衛(wèi)星位點為高度多態(tài)性,表明桂海1號凡納濱對蝦群體尚具有較豐富的遺傳變異。

    H是衡量群體遺傳變異的最適參數(shù),是指群體中某個位點上的雜合子頻率,主要反映群體的遺傳變異程度(Shikano and Taniguchi,2002);PIC與基因突變率有關(Gyapay et al.,1994),因此,通過H和PIC能反映群體內個體遺傳變異程度,其數(shù)值高低與遺傳變異程度呈正相關。本研究中,不同養(yǎng)殖家系凡納濱對蝦的平均Ho為0.511、平均He為0.451,群體的雜合度及遺傳多樣處于中等水平。雜合子缺失現(xiàn)象的產(chǎn)生可能是研究對象數(shù)量有限、近親雜交或人為干擾等因素導致稀有堿基缺失所致(Antoro et al.,2006)。Hardy-Weinberg平衡檢驗是評估樣本代表性及其基因分型質量的重要工具(Crane et al.,2004),樣本群體小、近親交配、等位基因突變等均可導致所檢測的種群基因型偏離Hardy-Weinberg平衡(宋煒等,2017)。本研究所分析的16個微衛(wèi)星多態(tài)性位點在16尾不同養(yǎng)殖家系凡納濱對蝦中,有12個微衛(wèi)星位點滿足Hardy-Weinberg平衡,有4個微衛(wèi)星點位偏離Hardy-Weinberg平衡,表明不同養(yǎng)殖家系凡納濱對蝦群體遺傳結構尚處于相對穩(wěn)定狀態(tài)。

    SMRT測序是近年來新出現(xiàn)的第三代測序技術,是通過觀測單個DNA聚合酶合成過程而實現(xiàn)單分子實時測序的技術。SMRT測序不僅具有測序長度優(yōu)勢,還不存在GC偏好(Roberts et al.,2013),因此開發(fā)的微衛(wèi)星分子標記在基因組中的分布相對較均勻,而有利于構建遺傳連鎖圖譜。本研究首次利用SMRT測序開發(fā)凡納濱對蝦微衛(wèi)星分子標記,有效增加了凡納濱對蝦的可用微衛(wèi)星分子標記數(shù)量,同時證實采用SMRT測序開發(fā)微衛(wèi)星分子標記是一種高效可行的方法。

    4 結論

    通過對凡納濱對蝦轉錄組進行SMRT測序、Illumina測序糾錯及MISA分析,共獲得39674條微衛(wèi)星序列,從設計的39367對微衛(wèi)星引物中隨機選取40對微衛(wèi)星引物進行驗證,發(fā)現(xiàn)有16對微衛(wèi)星引物顯示出多態(tài)性,說明采用SMRT測序開發(fā)獲得的微衛(wèi)星分子標記有相當高的比例可用于凡納濱對蝦遺傳學研究。即采用SMRT測序開發(fā)凡納濱對蝦微衛(wèi)星分子標記是一種簡單而高效的途徑,有助于開展凡納濱對蝦的遺傳育種研究工作。

    參考文獻:

    陳錨,吳長功,相建海,黃皓,劉小林,劉翠紅,何建國,孫成波. 2008. 凡納濱對蝦的選育與家系的建立[J]. 海洋科學,32(11):5-8. [Chen M,Wu C G,Xiang J H,Huang H,Liu X L,Liu C H,He J G,Sun C B. 2008. Selective breeding and pedigree foundation of Litopenaeus vannamei[J]. Marine Sciences,32(11):5-8.]

    黃薇,李劍鋒,張根玉,李永權,施永海,張海明,徐嘉波,駱志強. 2014. 池塘溫室南美白對蝦、羅氏沼蝦三茬輪養(yǎng)技術[J]. 中國水產(chǎn),(6):69-71. [Huang W,Li J F,Zhang G Y,Li Y Q,Shi Y H,Zhang H M,Xu J B,Luo Z Q. 2014. Three-crop rotation culture of Litopenaeus vannamei and Macrobrachium rosenbergii in pond greenhouse[J]. China Fisheries,(6):69-71.]

    凌士鵬,孫萍,林賢銳,沈建生. 2018. 基于SSR標記的桃種質資源遺傳多樣性研究[J]. 江西農業(yè)學報,30(11):14-18. [Ling S P,Sun P,Lin X R,Shen J S. 2018. A study on genetic diversity of peach germplasm resources based on SSR markers[J]. Acta Agriculturae Jiangxi,30(11):14-18.]

    馬寧,曾地剛. 2013. 凡納濱對蝦微衛(wèi)星序列的篩選[J]. 西南農業(yè)學報,26(6):2629-2633. [Ma N,Zeng D G. 2013. Isolation of microsatellite sequences from Litopenaeus vannamei[J]. Southwest China Journal of Agricultural Sciences,26(6):2629-2633.]

    盛巖,鄭蔚虹,裴克全,馬克平. 2002. 微衛(wèi)星標記在種群生物學研究中的應用[J]. 植物生態(tài)學報,26(S):119-126. [Sheng Y,Zheng W H,Pei K Q,Ma K P. 2002. Applications of micros atellites in population biology[J]. Acta Phytoecologica Sinica,26(S):119-126.]

    宋煒,孟永永,蔣科技,趙明,任桂靜,張鳳英,馬凌波. 2017. 棘頭梅童魚七個野生群體遺傳多樣性的微衛(wèi)星分析[J]. 水產(chǎn)學報,41(1):31-39. [Song W,Meng Y Y,Jiang K J,Zhao M,Ren G J,Zhang F Y,Ma L B. 2017. Analysis of genetic diversity among seven wild Collichthys lucidus populations by using microsatellite marker[J]. Journal of Fisheries of China,41(1):31-39.]

    孫效文,張曉鋒,趙瑩瑩,張研,賈智英,常玉梅,魯翠云,梁利群. 2008. 水產(chǎn)生物微衛(wèi)星標記技術研究進展及其應用[J]. 中國水產(chǎn)科學,15(4):689-703. [Sun X W,Zhang X F,Zhao Y Y,Zhang Y,Jia Z Y,Chang Y M,Lu C Y,Liang L Q. 2008. Development and application of microsatellite markers in aquatic species[J]. Journal of Fishery Sciences of China,15(4):689-703.]

    王興強,馬甡,董雙林. 2004. 凡納濱對蝦生物學及養(yǎng)殖生態(tài)學研究進展[J]. 海洋湖沼通報,(4):94-100. [Wang X Q,Ma S,Dong S L. 2004. Studies on the biology and cultural ecology of Litopenaeus vannamei:A review[J]. Transaction of Oceanology and Limnology,(4):94-100.]

    熊建華,趙永貞,高永華,謝達祥,張彬,陳曉漢. 2011. 凡納濱對蝦良種培育研究進展[J]. 南方農業(yè)學報,42(5):556-561. [Xiong J H,Zhao Y Z,Gao Y H,Xie D X,Zhang B,Chen X H. 2011. Advances in breeding improved varie-ties of Litopenaeus vannamei[J]. Journal of Southern Agriculture,42(5):556-561.]

    楊銘,于洋,張曉軍,王全超,劉敬文,李富花,相建海. 2017. 基于轉錄組數(shù)據(jù)的凡納濱對蝦微衛(wèi)星標記開發(fā)[J]. 海洋科學,41(2):96-102. [Yang M,Yu Y,Zhang X J,Wang Q C,Liu J W,Li F H,Xiang J H. 2017. Development of microsatellite markers from the transcriptome sequences of Pacific white shrimp (Litopenaus vannamei)[J]. Marine Sciences,41(2):96-102.]

    張麗穎,張忠旭,郝憲斌,姚繼攀. 2009. 微衛(wèi)星標記技術在水稻抗病育種中的應用[J]. 科技信息,(19):767. [Zhang L Y,Zhang Z X,Hao X B,Yao J P. 2009. Application of microsatellite marker technology in rice disease resistance breeding[J]. Science & Technology Information,(19):767.]

    周大顏,張志新,黃彩林,招志杰,莫飛龍. 2019. 3種白鯽雜交子代的轉錄組學分析[J]. 南方農業(yè)學報,50(6):1328-1338. [Zhou D Y,Zhang Z X,Huang C L,Zhao Z J,Mo F L. 2019. Transcriptomic analysis of three Carassius auratus cuvieri hybrids[J]. Journal of Southern Agriculture,50(6):1328-1338.

    Antoro S,Na-Nakorn U,Koedprang W. 2006. Study of gene-tic diversity of orange-spotted grouper,Epinephelus coioides,from Thailand and Indonesia using microsatellite markers[J]. Marine Biotechnology(New York),8(1):17-26.

    Botstein D,White R L,Skolnick M,Davis R W. 1980. Construction genetic linkage map in man using restriction fragment length polymorphisms[J]. American Journal of Human Genetics,32(3):314-331.

    Crane P A,Lewis C J,Kretschmer E J,Miller S J,Spearman W J,DeCicco A L,Lisac M J,Wenburg J K. 2004. Chara-cterization and inheritance of seven microsatellite loci from Dolly Varden,Salvelinus malma,and cross-species amplification in Arctic char,S. alpinus[J]. Conservation Genetics,5:737-741.

    Eid J,F(xiàn)ehr A,Gray J,Luong K,Lyle J,Otto G,Peluso P,Rank D,Baybayan P,Bettman B,Bibillo A,Bjornson K,Chaudhuri B,Christians F,Cicero R,Clark S,Dalal R,Dewinter A,Dixon J,F(xiàn)oquet M,Gaertner A,Hardenbol P,Heiner C,Hester K,Holden D,Kearns G,Kong X, Kuse R,Lacroix Y,Lin S,Lundquist P,Ma C,Marks P,Maxham M,Murphy D,Park I,Pham T,Phillips M,Roy J,Sebra R,Shen G,Sorenson J,Tomaney A,Travers K,Trulson M,Vieceli J,Wegener J,Wu D,Yang A,Zaccarin D,Zhao P,Zhong F,Korlach J,Turner S. 2009. Real-time DNA sequencing from single polymerase molecules[J]. Science,323(5910):133-138.

    Garcia D K,Dhar A K,Alcivar-Warren A. 1996. Molecular analysis of a RAPD marker (B20) reveals two microsa-tellites and differential mRNA expression in Penaeus vannamei[J]. Molecular Marine Biology and Biotechnology,5(1):71-83.

    Gyapay G,Morissette J,Vignal A,Dib C,F(xiàn)izames C,Millas-seau P,Marc S,Bernardi G,Lathrop M,Weissenbach J. 1994. The 1993-94 Généthon human genetic linkage map[J]. Nature Genetics,7(2):246-339.

    Lightner D V. 2011. Virus diseases of farmed shrimp in the Western Hemisphere (the Americas):A review[J]. Journal of Invertebrate Pathology,106(1):110-130.

    Lu X,Kong J,Meng X H,Cao B X,Luo K,Dai P,Luan S. 2018. Identification of SNP markers associated with to-lerance to ammonia toxicity by selective genotyping from de novo assembled transcriptome in Litopenaeus vannamei[J]. Fish & Shellfish Immunology,73:158-166.

    Meehan D,Xu Z,Zuniga G,Alcivar-Warren A. 2003. High frequency and large number of polymorphic microsate-llites in cultured shrimp,Penaeus (Litopenaeus) vannamei[Crustacea:Decapoda][J]. Marine Biotechnology(New York),5(4):311-330.

    Roberts R J,Carneiro M O,Schatz M C. 2013. The advanta-ges of SMRT sequencing[J]. Genome Biology,14(7):405.

    Shikano T,Taniguchi N. 2002. Relationships between genetic variation measured by microsatellite DNA markers and a fitness-related trait in the guppy(Poecilia reticulata)[J]. Aquaculture,209(1):77-90.

    Santos C A,Andrade S C S,F(xiàn)reitas P D. 2018. Identification of SNPs potentially related to immune responses and growth performance in Litopenaeus vannamei by RNA-seq analyses[J]. PeerJ,6:e5154. doi:10.7717/peerj.5154.

    Wajid A,Wasim M,Yaqub T,F(xiàn)iryal S,Tayyab M,Siddique S,Hussain T. 2014. Assessment of genetic diversity in Balochi and Rakhshani sheep breeds of Balochistan using microsatellite DNA markers[J]. The Journal of Animal & Plant Sciences,24(5):1348-1354.

    Yu Y,Wei J K,Zhang X J,Liu J W,Liu C Z,Li F H,Xiang J H. 2014. SNP discovery in the transcriptome of white Pacific shrimp Litopenaeus vannamei by next generation sequencing[J]. PLoS One,9(1):e87218.

    猜你喜歡
    測序
    杰 Sir 帶你認識宏基因二代測序(mNGS)
    新民周刊(2022年27期)2022-08-01 07:04:49
    二代測序協(xié)助診斷AIDS合并馬爾尼菲籃狀菌腦膜炎1例
    傳染病信息(2021年6期)2021-02-12 01:52:58
    外顯子組測序助力產(chǎn)前診斷胎兒骨骼發(fā)育不良
    中草藥DNA條形碼高通量基因測序一體機驗收會在京召開
    基因測序技術研究進展
    中美聯(lián)合成功研發(fā)全球最準三代基因測序儀
    綿羊骨骼肌轉錄組高通量測序從頭組裝和特征分析
    DNA微陣列芯片法與直接測序法檢測CYP2C19基因型的比較研究
    基因捕獲測序診斷血癌
    單細胞測序技術研究進展
    网址你懂的国产日韩在线| 特大巨黑吊av在线直播| 一边摸一边抽搐一进一小说| 九九热线精品视视频播放| 久久久久久久久大av| 亚洲真实伦在线观看| 女同久久另类99精品国产91| 国内精品一区二区在线观看| 亚洲av中文字字幕乱码综合| 亚洲人成网站高清观看| 99热这里只有精品一区| 国产精品伦人一区二区| 最近的中文字幕免费完整| 18禁在线无遮挡免费观看视频 | 69av精品久久久久久| 狂野欧美激情性xxxx在线观看| 丰满人妻一区二区三区视频av| 欧美日韩在线观看h| 久久韩国三级中文字幕| 日韩欧美一区二区三区在线观看| 亚洲成人精品中文字幕电影| ponron亚洲| 国产一级毛片七仙女欲春2| 高清午夜精品一区二区三区 | 搞女人的毛片| 天堂av国产一区二区熟女人妻| 午夜精品一区二区三区免费看| 极品教师在线视频| 久久99热6这里只有精品| 亚洲国产色片| 日本欧美国产在线视频| 18禁裸乳无遮挡免费网站照片| av卡一久久| 观看美女的网站| 亚州av有码| 高清日韩中文字幕在线| 一进一出好大好爽视频| 日韩制服骚丝袜av| 黄片wwwwww| 日韩精品有码人妻一区| 黄片wwwwww| 男人舔女人下体高潮全视频| 国产淫片久久久久久久久| 日韩精品有码人妻一区| 久久精品国产亚洲av香蕉五月| 国产 一区精品| 2021天堂中文幕一二区在线观| ponron亚洲| 男女啪啪激烈高潮av片| 精品人妻偷拍中文字幕| 亚洲一区二区三区色噜噜| 午夜福利在线观看免费完整高清在 | 插阴视频在线观看视频| 插阴视频在线观看视频| 国产精华一区二区三区| 男插女下体视频免费在线播放| 亚洲婷婷狠狠爱综合网| av在线播放精品| 国产成人a∨麻豆精品| 亚洲最大成人av| 丰满乱子伦码专区| 99在线人妻在线中文字幕| 国产精品人妻久久久影院| 亚州av有码| 久久精品影院6| 99在线人妻在线中文字幕| 男女边吃奶边做爰视频| 丰满乱子伦码专区| 黄色视频,在线免费观看| 丰满乱子伦码专区| 国产一区二区亚洲精品在线观看| 色哟哟·www| 亚洲性夜色夜夜综合| 舔av片在线| 日韩强制内射视频| 黄色视频,在线免费观看| 舔av片在线| 日韩欧美免费精品| 少妇裸体淫交视频免费看高清| 综合色丁香网| 日日干狠狠操夜夜爽| 夜夜爽天天搞| 一个人观看的视频www高清免费观看| 男女视频在线观看网站免费| 97热精品久久久久久| 国产亚洲av嫩草精品影院| 日本欧美国产在线视频| 熟女人妻精品中文字幕| 99热这里只有是精品在线观看| 老司机午夜福利在线观看视频| 校园春色视频在线观看| 大又大粗又爽又黄少妇毛片口| 国产精品人妻久久久久久| 看免费成人av毛片| .国产精品久久| 寂寞人妻少妇视频99o| 男女啪啪激烈高潮av片| 搡老岳熟女国产| 搡女人真爽免费视频火全软件 | 国产精品久久久久久av不卡| 成人无遮挡网站| 男人的好看免费观看在线视频| 一本一本综合久久| 欧美色视频一区免费| 国产精品一区www在线观看| 国产中年淑女户外野战色| 成人鲁丝片一二三区免费| 人妻久久中文字幕网| 欧美成人免费av一区二区三区| 国产91av在线免费观看| 国产av一区在线观看免费| 成人美女网站在线观看视频| 一进一出抽搐动态| 精品一区二区三区视频在线| 欧美3d第一页| 亚洲第一电影网av| 午夜精品一区二区三区免费看| 人人妻人人看人人澡| 久久久久久久久久成人| 国产精品国产高清国产av| 在线观看免费视频日本深夜| 99久久久亚洲精品蜜臀av| 大型黄色视频在线免费观看| 国产av在哪里看| www日本黄色视频网| 国产精品av视频在线免费观看| 香蕉av资源在线| 天堂√8在线中文| 99九九线精品视频在线观看视频| 精品免费久久久久久久清纯| 性插视频无遮挡在线免费观看| 国产亚洲精品久久久com| 1024手机看黄色片| 久久精品国产亚洲网站| 成年女人毛片免费观看观看9| 精品99又大又爽又粗少妇毛片| 一本精品99久久精品77| 免费人成在线观看视频色| 尾随美女入室| 人人妻人人澡人人爽人人夜夜 | 淫妇啪啪啪对白视频| 美女被艹到高潮喷水动态| 亚洲欧美成人精品一区二区| 日韩,欧美,国产一区二区三区 | 丝袜美腿在线中文| 亚洲精品国产av成人精品 | 欧美日韩综合久久久久久| 国内少妇人妻偷人精品xxx网站| 2021天堂中文幕一二区在线观| 成人国产麻豆网| 精品午夜福利在线看| 亚州av有码| 少妇的逼好多水| 久久久精品欧美日韩精品| 99九九线精品视频在线观看视频| 超碰av人人做人人爽久久| 亚洲av第一区精品v没综合| 久久久久国内视频| 精品一区二区三区视频在线| 蜜桃亚洲精品一区二区三区| 久久久精品大字幕| 亚洲性久久影院| 99视频精品全部免费 在线| 国产乱人视频| 搡老岳熟女国产| 97超视频在线观看视频| 激情 狠狠 欧美| 色播亚洲综合网| 老司机影院成人| 亚洲一区高清亚洲精品| 日韩欧美国产在线观看| 赤兔流量卡办理| 精品福利观看| 天天躁夜夜躁狠狠久久av| 亚洲精品成人久久久久久| 亚洲综合色惰| 国产一区二区三区在线臀色熟女| 欧美成人a在线观看| 三级毛片av免费| 99久久成人亚洲精品观看| 成人精品一区二区免费| 亚洲av五月六月丁香网| 老熟妇仑乱视频hdxx| 久久人人精品亚洲av| 亚洲综合色惰| 18禁在线播放成人免费| 99久久精品国产国产毛片| 午夜激情福利司机影院| 日韩欧美国产在线观看| 色5月婷婷丁香| 欧美三级亚洲精品| 亚洲精品日韩av片在线观看| 国产男人的电影天堂91| а√天堂www在线а√下载| 成人性生交大片免费视频hd| 久久精品国产99精品国产亚洲性色| 日本一本二区三区精品| 久久久久久久久久成人| 夜夜看夜夜爽夜夜摸| 成人三级黄色视频| 国产黄a三级三级三级人| 观看美女的网站| 午夜免费激情av| 啦啦啦啦在线视频资源| 淫妇啪啪啪对白视频| 亚洲av五月六月丁香网| 久久欧美精品欧美久久欧美| 国产精品免费一区二区三区在线| 亚洲性夜色夜夜综合| 男人的好看免费观看在线视频| 干丝袜人妻中文字幕| 久久久久久久亚洲中文字幕| 少妇熟女aⅴ在线视频| 12—13女人毛片做爰片一| 亚洲成人中文字幕在线播放| 亚洲国产欧美人成| videossex国产| 99热6这里只有精品| 岛国在线免费视频观看| 一边摸一边抽搐一进一小说| 日本撒尿小便嘘嘘汇集6| 亚洲欧美精品综合久久99| 精品久久久久久久人妻蜜臀av| 亚洲国产精品合色在线| 国产男靠女视频免费网站| 99热只有精品国产| 亚洲aⅴ乱码一区二区在线播放| 国产精品嫩草影院av在线观看| 免费黄网站久久成人精品| 18+在线观看网站| 国产高潮美女av| 男人舔女人下体高潮全视频| 91久久精品国产一区二区三区| 高清毛片免费观看视频网站| 成人二区视频| 天堂av国产一区二区熟女人妻| 国产视频一区二区在线看| 日本五十路高清| 亚洲18禁久久av| 我要搜黄色片| 亚洲国产精品成人久久小说 | 最近在线观看免费完整版| 在线播放国产精品三级| 伦精品一区二区三区| 亚洲国产高清在线一区二区三| 久久精品久久久久久噜噜老黄 | 国产欧美日韩精品一区二区| 国产黄片美女视频| 亚洲美女黄片视频| 久久精品影院6| 午夜福利视频1000在线观看| 久久草成人影院| 午夜视频国产福利| 永久网站在线| 国产不卡一卡二| 乱人视频在线观看| 久久国内精品自在自线图片| 晚上一个人看的免费电影| 国产精品野战在线观看| 在线a可以看的网站| 美女大奶头视频| 中文字幕熟女人妻在线| 1024手机看黄色片| 国产大屁股一区二区在线视频| 伦理电影大哥的女人| 亚洲精品亚洲一区二区| 国产一区二区激情短视频| 午夜福利视频1000在线观看| 别揉我奶头 嗯啊视频| 一个人看视频在线观看www免费| 国产男靠女视频免费网站| 我要看日韩黄色一级片| 美女cb高潮喷水在线观看| 尤物成人国产欧美一区二区三区| 少妇高潮的动态图| 欧美日韩精品成人综合77777| 我要看日韩黄色一级片| 午夜激情福利司机影院| 免费一级毛片在线播放高清视频| 免费在线观看影片大全网站| 亚洲一级一片aⅴ在线观看| av在线播放精品| 国产v大片淫在线免费观看| 一个人看视频在线观看www免费| 亚洲乱码一区二区免费版| 我要看日韩黄色一级片| 国产视频内射| 少妇的逼水好多| 日本免费一区二区三区高清不卡| 国产国拍精品亚洲av在线观看| 18禁裸乳无遮挡免费网站照片| 日韩欧美国产在线观看| 国产精品1区2区在线观看.| 亚洲成人中文字幕在线播放| 日本免费a在线| 久久精品国产亚洲网站| 成人精品一区二区免费| 亚洲激情五月婷婷啪啪| 亚洲精品粉嫩美女一区| 国产亚洲精品久久久com| 男女那种视频在线观看| 亚洲电影在线观看av| 亚洲自偷自拍三级| 禁无遮挡网站| 亚洲精品456在线播放app| 午夜久久久久精精品| 三级经典国产精品| 精品一区二区三区人妻视频| 久久久久久久久久黄片| 午夜老司机福利剧场| 精品人妻一区二区三区麻豆 | 亚洲精品一区av在线观看| 嫩草影院精品99| 最新在线观看一区二区三区| 久久久久久久久久成人| 亚洲成人精品中文字幕电影| 免费观看人在逋| 国产v大片淫在线免费观看| 精品少妇黑人巨大在线播放 | 91精品国产九色| 又黄又爽又免费观看的视频| 久久亚洲精品不卡| 美女内射精品一级片tv| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 婷婷色综合大香蕉| 免费人成视频x8x8入口观看| 久久久久久久久久久丰满| 亚洲成人av在线免费| 日日撸夜夜添| 午夜福利视频1000在线观看| 亚洲av成人精品一区久久| 三级国产精品欧美在线观看| 床上黄色一级片| 国产中年淑女户外野战色| 尾随美女入室| 精华霜和精华液先用哪个| 少妇猛男粗大的猛烈进出视频 | 亚洲在线观看片| 亚洲色图av天堂| 无遮挡黄片免费观看| 超碰av人人做人人爽久久| 国产精华一区二区三区| 午夜精品一区二区三区免费看| 国内少妇人妻偷人精品xxx网站| a级毛片a级免费在线| 日本黄色片子视频| 日本熟妇午夜| 国产私拍福利视频在线观看| 日韩亚洲欧美综合| 波多野结衣巨乳人妻| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 国产 一区 欧美 日韩| 九九在线视频观看精品| 免费不卡的大黄色大毛片视频在线观看 | 男女下面进入的视频免费午夜| 久久韩国三级中文字幕| 国产精品一区二区免费欧美| ponron亚洲| 俄罗斯特黄特色一大片| 一级av片app| 国产黄片美女视频| 97超视频在线观看视频| 亚洲国产色片| 久久久久免费精品人妻一区二区| 无遮挡黄片免费观看| 性插视频无遮挡在线免费观看| 久久久久久久亚洲中文字幕| 国产精品综合久久久久久久免费| 久久久久久国产a免费观看| 成人鲁丝片一二三区免费| 成年版毛片免费区| 黑人高潮一二区| 少妇人妻精品综合一区二区 | 51国产日韩欧美| 在线国产一区二区在线| 午夜影院日韩av| 久久韩国三级中文字幕| 观看免费一级毛片| av在线蜜桃| 长腿黑丝高跟| 久久99热6这里只有精品| 成人永久免费在线观看视频| 亚洲成人中文字幕在线播放| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 国产乱人视频| 岛国在线免费视频观看| 丝袜喷水一区| 91午夜精品亚洲一区二区三区| 精品免费久久久久久久清纯| 一本精品99久久精品77| 天天躁夜夜躁狠狠久久av| 我要搜黄色片| 中文字幕人妻熟人妻熟丝袜美| 1000部很黄的大片| 国产精品三级大全| 麻豆av噜噜一区二区三区| 欧美最新免费一区二区三区| 高清日韩中文字幕在线| 亚洲人与动物交配视频| 91在线观看av| 最新在线观看一区二区三区| 久久久欧美国产精品| 成人三级黄色视频| 性色avwww在线观看| 日韩一区二区视频免费看| 国产精品无大码| 亚洲五月天丁香| 成人欧美大片| 亚洲精品成人久久久久久| 日本一二三区视频观看| av在线亚洲专区| 久久99热这里只有精品18| 99在线视频只有这里精品首页| 91在线观看av| 久久久久久伊人网av| 国产亚洲欧美98| 国产亚洲精品久久久久久毛片| 亚洲无线在线观看| 亚洲不卡免费看| 免费在线观看影片大全网站| 欧美性感艳星| 别揉我奶头 嗯啊视频| 久久久久久九九精品二区国产| 岛国在线免费视频观看| 国产欧美日韩精品一区二区| 欧美bdsm另类| 午夜精品在线福利| 欧美高清成人免费视频www| 一进一出好大好爽视频| 在线免费观看不下载黄p国产| 亚洲欧美日韩东京热| 国产欧美日韩精品亚洲av| 老司机福利观看| 国产成人精品久久久久久| 我要搜黄色片| 欧美一级a爱片免费观看看| av在线观看视频网站免费| 国产成人影院久久av| 久久精品国产亚洲av涩爱 | 日本 av在线| 亚洲成a人片在线一区二区| 国产熟女欧美一区二区| 国国产精品蜜臀av免费| 人妻久久中文字幕网| 日本免费一区二区三区高清不卡| 狠狠狠狠99中文字幕| 啦啦啦韩国在线观看视频| avwww免费| 最近的中文字幕免费完整| 欧美日韩综合久久久久久| 亚洲精品色激情综合| 日本免费a在线| 成人一区二区视频在线观看| av国产免费在线观看| 亚洲成人久久爱视频| 婷婷精品国产亚洲av| 淫秽高清视频在线观看| 国产又黄又爽又无遮挡在线| 久久欧美精品欧美久久欧美| 久久久久国内视频| 欧美绝顶高潮抽搐喷水| 国产一区亚洲一区在线观看| 色av中文字幕| 嫩草影视91久久| 高清毛片免费观看视频网站| 亚洲最大成人av| 亚洲人成网站在线观看播放| 亚洲精品粉嫩美女一区| 欧美潮喷喷水| 99热只有精品国产| 亚洲一区高清亚洲精品| 国产精品一区二区性色av| 免费人成视频x8x8入口观看| 免费av观看视频| 99热只有精品国产| 亚洲av电影不卡..在线观看| 两个人的视频大全免费| 嫩草影院入口| 特级一级黄色大片| 国产精品日韩av在线免费观看| 成年av动漫网址| 日本五十路高清| 国产欧美日韩一区二区精品| 欧美区成人在线视频| 欧美日韩乱码在线| 真人做人爱边吃奶动态| 国产精品免费一区二区三区在线| 性欧美人与动物交配| 国产高清视频在线播放一区| 97热精品久久久久久| 日韩,欧美,国产一区二区三区 | 国产欧美日韩精品一区二区| www.色视频.com| 成人亚洲精品av一区二区| 久久欧美精品欧美久久欧美| 久久久久九九精品影院| 内射极品少妇av片p| 综合色av麻豆| 成人综合一区亚洲| 丰满的人妻完整版| 晚上一个人看的免费电影| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 成人亚洲欧美一区二区av| 一级毛片久久久久久久久女| 午夜视频国产福利| 十八禁网站免费在线| 国产综合懂色| 久久精品夜夜夜夜夜久久蜜豆| 99热全是精品| 亚洲精品亚洲一区二区| 最近在线观看免费完整版| av福利片在线观看| 岛国在线免费视频观看| 久久久久久久久久成人| a级毛片免费高清观看在线播放| 久久久a久久爽久久v久久| 热99在线观看视频| 国产视频内射| 91久久精品电影网| 国产精华一区二区三区| 老熟妇乱子伦视频在线观看| 日本爱情动作片www.在线观看 | 国产一区亚洲一区在线观看| 日韩欧美国产在线观看| 午夜a级毛片| 欧美丝袜亚洲另类| 久久精品人妻少妇| 春色校园在线视频观看| 国产白丝娇喘喷水9色精品| 日韩成人av中文字幕在线观看 | 99热只有精品国产| 国产在视频线在精品| 国产精品av视频在线免费观看| 人妻久久中文字幕网| 亚洲成人中文字幕在线播放| 啦啦啦啦在线视频资源| 亚洲精品日韩在线中文字幕 | 亚洲四区av| 国产高清不卡午夜福利| 久久6这里有精品| 特级一级黄色大片| 少妇裸体淫交视频免费看高清| 午夜日韩欧美国产| АⅤ资源中文在线天堂| 老女人水多毛片| 国产私拍福利视频在线观看| 精品人妻一区二区三区麻豆 | 岛国在线免费视频观看| 91狼人影院| 中文字幕熟女人妻在线| 看黄色毛片网站| 亚洲无线在线观看| 看黄色毛片网站| 国产精品爽爽va在线观看网站| 真人做人爱边吃奶动态| 久久精品综合一区二区三区| 国产伦精品一区二区三区四那| 99在线人妻在线中文字幕| 最近最新中文字幕大全电影3| 国产精品一及| 老女人水多毛片| 在线免费观看的www视频| 国内久久婷婷六月综合欲色啪| 黄片wwwwww| 一级av片app| av国产免费在线观看| 女的被弄到高潮叫床怎么办| 国产一区二区在线观看日韩| 女人被狂操c到高潮| 午夜影院日韩av| 麻豆精品久久久久久蜜桃| av中文乱码字幕在线| 99热这里只有是精品在线观看| av天堂在线播放| 国产69精品久久久久777片| 男人舔奶头视频| 亚洲成人久久爱视频| 嫩草影视91久久| 中文字幕精品亚洲无线码一区| 日韩精品青青久久久久久| 狂野欧美白嫩少妇大欣赏| 国产精品综合久久久久久久免费| 亚洲欧美中文字幕日韩二区| 日本免费一区二区三区高清不卡| 最近2019中文字幕mv第一页| 99久久久亚洲精品蜜臀av| 国产精品一区www在线观看| 国产高清视频在线播放一区| 欧美色视频一区免费| 欧美一区二区精品小视频在线| 日韩 亚洲 欧美在线| 日日啪夜夜撸| 欧美3d第一页| 3wmmmm亚洲av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产清高在天天线| 真实男女啪啪啪动态图| 国产男人的电影天堂91| 成人欧美大片| 久久久精品94久久精品| av在线亚洲专区| 一夜夜www| 波野结衣二区三区在线| 男插女下体视频免费在线播放| 精品一区二区三区人妻视频| 乱系列少妇在线播放| 免费看a级黄色片| 色av中文字幕| 日韩欧美一区二区三区在线观看| 精品久久久久久久久久免费视频| 九色成人免费人妻av|