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

    克什米爾熊蜂轉(zhuǎn)錄組分析及SSR位點開發(fā)

    2024-01-01 00:00:00孫國閆京艷梁程博馬曉璇劉道鑫
    草地學報 2024年6期

    摘要:作為眾多農(nóng)作物和野生植物的重要傳粉者,熊蜂在維持農(nóng)業(yè)和生態(tài)系統(tǒng)平衡中發(fā)揮十分重要的作用。本研究以克什米爾熊蜂(Bombus kashmirensis)為研究對象,使用Illumina Novaseq平臺進行轉(zhuǎn)錄組測序,并基于轉(zhuǎn)錄組數(shù)據(jù)進行微衛(wèi)星位點的開發(fā)。結(jié)果顯示,克什米爾熊蜂轉(zhuǎn)錄組測序共獲得95 639條測序良好的unigene,GO注釋和KEGG通路富集主要集中在分子功能、細胞進程和碳水化合物代謝上?;谵D(zhuǎn)錄組數(shù)據(jù)共獲得27 414個微衛(wèi)星位點,位點發(fā)生頻率17.81%。簡單重復序列主要重復類型為單堿基,其中A/T為優(yōu)勢基元,占總微衛(wèi)星數(shù)量的55.89%,且微衛(wèi)星主要重復次數(shù)范圍在11~15次(48.61%)。用10頭不同克什米爾熊蜂樣品對30個微衛(wèi)星位點進行驗證,有17對引物可穩(wěn)定擴增。本研究為克什米爾熊蜂的適應與進化等種群遺傳學研究提供了重要參考。

    關鍵詞:克什米爾熊蜂;轉(zhuǎn)錄組測序;微衛(wèi)星分析

    中圖分類號:Q969.557.7""" 文獻標識碼:A"""" 文章編號:1007-0435(2024)06-1742-10

    Transcriptome Analysis and SSR Locus Development in Bombus kashmirensis

    SUN Guo, YAN Jing-yan, LIANG Cheng-bo, MA Xiao-xuan, LIU Dao-xin*

    (School of Agriculture and Animal Husbandry, Qinghai University, Xining, Qinghai Province 810016, China)

    Abstract:As the essential pollinators of many crops and wild plants,bumblebees play an important role in maintaining the balance of agriculture and ecosystems. In this study,the Bombus kashmirensis was used for transcriptome sequencing by using Illumina Novaseq platform and microsatellite locus were development based on the transcriptome data. The results indicated that a total of 95 639 well-sequenced unigenes were obtained from the transcriptome sequencing of Bombus kashmirensis,with GO annotations and KEGG pathway enrichment focusing on molecular functions,cellular processes and carbohydrate metabolism. A total of 27 414 microsatellite loci were obtained based on transcriptome data,with a locus frequency of 17.81%. The major repeat type of microsatellites was Mono-nucleotide repeats,of which A/T was the dominant motif,accounting for 55.89% of the total microsatellites number,and the major microsatellites base repeats were in the range of 11~15 repeats (48.61%). Thirty microsatellite loci were validated using samples from 10 different Bombus kashmirensis,of which 17 primer pairs were stably amplified. This study provides an important reference for population genetics studies such as adaptation and evolution of Bombus kashmirensis.

    Key words:Bombus kashmirensis;Transcriptome sequencing;Microsatellite analysis

    熊蜂(bombus spp.)隸屬于節(jié)肢動物門(Arthropoda)昆蟲綱(Insecta)膜翅目(Hymenoptera)蜜蜂科(Apidae)熊蜂屬(Bombus),是一類社會性傳粉昆蟲[1]。相比于蜜蜂,熊蜂體型較大,渾身長滿絨毛,耐寒性強,能夠分布在寒冷的區(qū)域為一些瀕危植物進行傳粉[2]。此外,由于具有耐濕性強、耐低光照以及聲震傳粉等特性,熊蜂成為了溫室作物的理想授粉昆蟲,尤其適宜為設施作物及深花冠植物授粉,能夠顯著提高農(nóng)林果蔬產(chǎn)量及品質(zhì)[3-4]。研究顯示,熊蜂授粉可使大棚試驗農(nóng)作物坐果率提高30%以上,單株產(chǎn)量提高17%以上[5]。目前全球已知熊蜂種類260余種,我國有120多種,約占全球熊蜂資源的50%[6-7]。然而,氣候變暖和頻繁的極端天氣導致全球熊蜂資源正在減少,包括我國熊蜂資源豐富的一些地區(qū)的熊蜂數(shù)量也在減少[8-9]。有研究表明,氣候變暖引起的植物單一化已使某些熊蜂的飲食寬度發(fā)生改變,導致它們與植物互作關系中傳粉功能的失配,已威脅到生態(tài)系統(tǒng)的平衡[10]。另外,農(nóng)藥的廣泛使用更是直接威脅到各類傳粉昆蟲的生存,不僅引起熊蜂訪花行為的異常,而且極大地造成其傳粉能力的下降[11]。

    為了維持熊蜂種群數(shù)量及多樣性,除了要采取各種措施保護熊蜂賴以生存的自然環(huán)境外,更多科學家使用分子手段為熊蜂的保護提供支持[12]。例如,Wang等[13]和Williams等[14]利用多種分子標記重建全球15個亞屬120多種熊蜂的系統(tǒng)發(fā)育關系,為熊蜂的分類提供了參考依據(jù);Lin等[15]和Sadd等[16]通過分析重要熊蜂物種的基因組序列,探討熊蜂基因組進化速率與海拔的關系以及熊蜂基因組進化過程中一些關鍵基因潛在的調(diào)控過程,為熊蜂的適應與進化奠定了遺傳學基礎。簡單重復序列(Single sequence reperts,SSR)又稱微衛(wèi)星序列(microsatellite),是研究動植物基因遺傳多樣性的一種重要分子標記[17]。作為一種多態(tài)性高、穩(wěn)定性好、共顯性表達且易檢測的分子遺傳標記[18],微衛(wèi)星被廣泛應用于遺傳多樣性分析[19]、基因定位[20]和親緣關系鑒定[21]等眾多領域。目前熊蜂有關微衛(wèi)星的研究主要集中在群體遺傳方面,如Takeuchi等[22]利用微衛(wèi)星標記研究北海道小峰熊蜂(Bombus hypocrita sapporoensis)的遺傳結(jié)構(gòu);Cejas等[23]利用微衛(wèi)星標記結(jié)合線粒體基因標記研究伊比利亞半島地熊蜂(Bombus terrestris)的時空格局等;而不同熊蜂轉(zhuǎn)錄組分析以及基于轉(zhuǎn)錄組數(shù)據(jù)開發(fā)SSR位點的報道很少。

    本研究對克什米爾熊蜂(Bombus kashmirensis)進行了轉(zhuǎn)錄組測序,分析了其基因表達特征,并基于轉(zhuǎn)錄組數(shù)據(jù)分析了SSR位點的分布特征,同時設計引物對部分SSR位點進行了有效性驗證。研究結(jié)果可以為克什米爾熊蜂的適應與進化等種群遺傳學研究提供參考。

    1 材料與方法

    1.1 供試昆蟲

    本試驗中克什米爾熊蜂樣品于2021年8—9月采集于三江源地區(qū)(31°39′~36°12′ N,89°45′~102°23′ E),將熊蜂樣品置于液氮速凍后,保存在-80℃條件下,以便于后續(xù)RNA提取。另取10頭工蜂用于微衛(wèi)星引物的驗證,使用Ezup 柱式動物基因組 DNA 抽提試劑盒(上海生工生物工程有限公司)提取DNA,經(jīng)過微量核酸分光光度計和1%的瓊脂糖凝膠電泳檢測,結(jié)果滿足進一步實驗要求,置于-20℃保存?zhèn)溆谩?/p>

    1.2 克什米爾熊蜂總RNA的提取、建庫及數(shù)據(jù)評估

    轉(zhuǎn)錄組測序委托上海美吉生物醫(yī)藥科技有限公司進行。取克什米爾熊蜂的頭、胸部肌肉且用MJZol total RNA extraction kit試劑盒(上海美吉生物醫(yī)藥科技有限公司)提取總RNA,使用Nanodrop 2000對所提RNA的濃度和純度進行檢測并用瓊脂糖凝膠電泳檢測RNA完整性;采用Illumina TruseqTM RNA sample prep Kit方法[24]進行文庫構(gòu)建;使用Illumina平臺完成轉(zhuǎn)錄組測序,并對獲得的測序數(shù)據(jù)進行質(zhì)量控制;獲得 RNA-seq高質(zhì)量測序數(shù)據(jù)后,使用Trinity軟件對轉(zhuǎn)錄組數(shù)據(jù)從頭組裝,并對初始組裝序列進行優(yōu)化過濾(TransRate軟件和CD-HIT軟件)并再次評估(BUSCO軟件)。

    1.3 克什米爾熊蜂轉(zhuǎn)錄組unigenes的功能注釋

    將轉(zhuǎn)錄組測序獲得的所有轉(zhuǎn)錄本與NR(https://www.ncbi.nlm.nih.gov/public/),Swiss-Prot(ftp://ftp.uniprot.org/pub/databases/uniprot/current_release/knowledgebase/complete/uniprot_sprot.fasta.gz),Pfam(http://pfam.xfam.org/),eggNOG(http://eggnogdb.embl.de/#/app/home),GO(http://www.geneontology.org/),KEGG(http://www.genome.jp/kegg/)等六大數(shù)據(jù)庫進行比對,獲得在各數(shù)據(jù)庫的注釋信息。

    1.4 SSR的篩選及引物的設計

    利用 MISA 軟件(http://pgrc.ipk-gatresleben.de/misa/)搜尋克什米爾熊蜂轉(zhuǎn)錄組組裝后 unigene 數(shù)據(jù)中潛在的SSR位點。軟件中參數(shù)設置為單堿基、二堿基、三堿基、四堿基、五堿基、六堿基的最短重復次數(shù)分別為10,6,5,5,5,5。將結(jié)果文件導入Primer3軟件,設計SSR位點的相關引物。

    1.5 SSR引物的驗證

    為了較全面的驗證引物的有效性,從批量設計的引物中隨機選出30對引物(表1)。(挑選原則為兩堿基重復的16個,三堿基重復的7個,四堿基重復的7個;位點片段長度在150 bp以下的位點10個、150~200 bp的位點10個、200~300 bp的位點10個;重復次數(shù)15~20次的5個、21~25次的5個、26~30次的5個、31~35次的5個、36次以上的5個;完美型微衛(wèi)星)交由生工生物工程(上海)股份有限公司合成。PCR反應體系為 25 μL:克什米爾熊蜂DNA模板 0.5 μL,上下游引物(0.2~1.0 μmol·L-1)各加 0.5 μL,2×Taq Master Mix* 12.5 μL,ddH2O 11 μL。反應程序:94℃預變性 1 min 30 s;94℃變性 20 s,59~60℃退火 20 s,72℃延伸 20 s,總共 35 個循環(huán);72℃延伸 5 min,4℃保溫。擴增產(chǎn)物用1.2%瓊脂糖凝膠電泳檢測(110V,35 min)。

    2 結(jié)果與分析

    2.1 克什米爾熊蜂轉(zhuǎn)錄組測序、組裝及注釋

    通過對克什米爾熊蜂樣品轉(zhuǎn)錄組測序,獲得119.45 Gb Clean Data,各樣品Clean Data均達到6.04 Gb以上,Q30堿基質(zhì)量比值在92.98%以上,GC含量在37.36%~47.39%之間,說明測序質(zhì)量較好。經(jīng)de novo組裝后,通用單拷貝同源基因基準值(Benchmarking universal single-copy orthologs,BUSCO)中有94.5%的測試基因被覆蓋,其中有78.9%的基因經(jīng)過比對發(fā)現(xiàn)是單拷貝,有15.6%的基因經(jīng)過比對發(fā)現(xiàn)包含多拷貝(測試基因覆蓋數(shù)=比對后的單拷貝數(shù)+比對后的多拷貝數(shù))[BUSCO值:C:94.5%(S:78.9%:D:15.6%),C=S+D],共獲得unigene 95 639條,平均長度為981.88 bp,N50長度為1 741 bp,unigene長度在200~500 bp之間,占總unigene的50%,隨著序列長度的增加,unigene總體呈逐步遞減趨勢(表2)。將克什米爾熊蜂轉(zhuǎn)錄組測序獲得的所有轉(zhuǎn)錄本分別與NR,Swiss-Prot,Pfam,eggNOG,GO 和 KEGG數(shù)據(jù)庫進行比對注釋,結(jié)果表明,在148 656條Transcript中,共有64.62%的Transcript被注釋到相應的數(shù)據(jù)庫,組裝評估后58.41%的基因被注釋到公共數(shù)據(jù)庫;同時有58.64%表達unigene 被注釋到公共數(shù)據(jù)庫中(表3)。轉(zhuǎn)錄組數(shù)據(jù)被存入NCBI SRA 數(shù)據(jù)庫中,克什米爾熊蜂的登錄號為SRR19913284 至SRR19913300。

    2.2 unigenes的GO功能注釋

    GO數(shù)據(jù)庫的28 443個unigenes的功能注釋結(jié)果表明,克什米爾熊蜂轉(zhuǎn)錄組功能分為分子功能(molecular function)、細胞組分(cellular component)和生物學進程(biological process)三大類和54個功能小類。注釋到的生物學進程有22個功能小類,注釋到的細胞組分有15個功能小類,注釋到的分子功能有17個功能小類。在生物學進程中注釋到的細胞進程(cellular process)功能的基因最多(13 091,36.46%),在細胞組分中注釋到的細胞部分(cell part)功能的基因最多(13 758,30.62%),在分子功能中注釋到的結(jié)合(binding)功能的基因最多(16 316,44.13%)(圖1)。

    2.3 unigenes的KEGG通路富集

    KEGG數(shù)據(jù)庫的31 212條unigenes 代謝通路富集分析表明,克什米爾熊蜂轉(zhuǎn)錄組代謝通路分別富集到環(huán)境信息處理(Envieonmental information processing)、新陳代謝(Metabolism)、遺傳信息處理(Genetic information processing)、細胞進程(Cellular processes)和有機系統(tǒng)(Organismal system)五大類別、33個小類別中。其中富集到新陳代謝的基因數(shù)最多(9 552),富集到環(huán)境信息處理中的基因最少(3 428)。在新陳代謝途徑中富集到氨基酸代謝(Amino acid metabolism)的基因最多(2 921,30.58%),在有機系統(tǒng)途徑中富集到內(nèi)分泌系統(tǒng)(Endocrine system)的基因最多(1 709,20.87%),在細胞進程途徑中富集到運輸與分解代謝(Transport and catabolism)的基因最多(2 223,44.65%),在遺傳信息處理途徑中富集到信號轉(zhuǎn)導(Signal transduction)的基因最多(3 107,90.64%)(圖2)。富集到的所有基因分布在342條代謝通路中,其中富集最多的10條代謝通路分別是核糖體(Ribosome)、戊糖和葡萄糖醛酸的相互轉(zhuǎn)化(Pentose and glucuronate interconversions)、內(nèi)質(zhì)網(wǎng)中的蛋白質(zhì)加工(Protein processing in endoplasmic reticulum)、氧化磷酸化(Oxidative phosphorylation)、產(chǎn)熱(Thermogenesis)、胞吞作用(Endocytosis)、剪接體(Spliceosome)、輔助因子的生物合成(Biosynthesis of cofactors)、泛素介導的蛋白水解(Ubiquitin mediated proteolysis)、AMPK 信號通路(AMPK signaling pathway)(表4)。

    2.4 克什米爾熊蜂轉(zhuǎn)錄組中SSR分布特征

    搜尋克什米爾熊蜂95 639條unigenes中的SSR位點,結(jié)果表明:共找到SSR點27 414個,分布在17 029個unigenes中,SSR位點發(fā)生的頻率為17.81%(SSR位點發(fā)生頻率=含有SSR位點的unigenes/全部序列數(shù)量),含有1個以上SSR位點的unigenes序列有4 109條,單堿基至六堿基的SSR位點依次減少,占比分別為56.26%,23.31%,17.08%,2.17%,0.61%,0.57%(表5)。

    單堿基至六堿基的堿基重復類型分布結(jié)果表明:單堿基、二堿基和三堿基重復基元總共16種。單堿基中共有2種重復基元,(A/T)n重復11~15次的最多(6 836);二堿基中共有4種重復基元,(AT/AT)n重復6~10次的最多(2 477);單堿基(A/T)n,(C/G)n與二堿基(AC/GT)n,(AG/CT)n,(AT/AT)n,(CG/CG)n重復均在5次以上。三堿基中共有10種重復基元,(AAG/CTT)n重復1~5次的最多(498)(圖3A)。四堿基重復基元共有29種,(ACAT/ATGT)n重復6~10次的最多(60)(圖3B)。五堿基重復基元共有51種,(AAAAG/CTTTT)n重復6~10次的最多(13)(圖3C),六堿基重復基元共有83種,(AAGATG/ATCTTC)n和(ACACGC/CGTGTG)n重復1~5次的最多,均為5次(圖3D)。三堿基至六堿基基元類型的重復次數(shù)在1~5,6~10,11~15和大于15次區(qū)間內(nèi)都有分布。

    基元重復次數(shù)在11~15次的最多,為13 327(48.61%)(百分數(shù)=該重復基元次數(shù)的SSR位點數(shù)/SSR位點總數(shù)),基元重復次數(shù)在1~5次的最少,為2 777(10.13%)。單堿基的優(yōu)勢重復基元為A/T,共出現(xiàn)15 321次(55.89%)(百分數(shù)=該堿基所有重復次數(shù)的SSR位點數(shù)/SSR位點總數(shù));二堿基的優(yōu)勢重復基元為AT/AT,共出現(xiàn)2 820次(10.29%);三堿基的優(yōu)勢重復基元為AAG/CTT,共出現(xiàn)1 033次(3.77%);四堿基的優(yōu)勢重復基元為AAAG/CTTT,共出現(xiàn)106次(0.61%);五堿基的重復基元為AAAAG/CTTTT,共出現(xiàn)24次(0.14%),六堿基由于各重復基元次數(shù)均未超過10次,優(yōu)勢重復基元不明顯。

    2.5 SSR分子標記開發(fā)

    按照1.5中的挑選原則,從篩選出的微衛(wèi)星引物中挑選30對引物進行驗證,每對引物分別用10個不同地點的克什米爾熊蜂樣品進行PCR擴增。結(jié)果顯示:有17對引物可以在10個不同樣品中穩(wěn)定擴增(表1前17對SSR引物),SSR引物的有效擴增率為56.67%,部分驗證結(jié)果如圖4所示。同時部分引物呈現(xiàn)一定的多態(tài)性,初步表明上述SSR引物可靠性較好。(由于是用10個不同采樣地的克什米爾熊蜂篩選SSR引物,部分引物出現(xiàn)1至2條非特異性或未擴增出來的條帶時可將其視為偶然情況,仍可用于后續(xù)分析。)

    3 討論

    3.1 克什米爾熊蜂轉(zhuǎn)錄組特點

    轉(zhuǎn)錄組測序技術通過對不同對象、不同發(fā)育時期的不同組織中的轉(zhuǎn)錄組進行分析,構(gòu)建基因調(diào)控網(wǎng)絡,篩選出關鍵基因,用于相關的研究[25]。目前,克什米爾熊蜂轉(zhuǎn)錄組信息處于研究空白。使用不同樣本材料測序,再將測序結(jié)果進行組裝,可以保證數(shù)據(jù)盡量覆蓋度廣、信息全面、準確度高[26]。本研究基于此方法對克什米爾熊蜂進行轉(zhuǎn)錄組測序,95 639條unigene被注釋到公共數(shù)據(jù)庫,占unigene 總數(shù)的58.41%。所有試驗樣本Q30堿基質(zhì)量比值均高于92.98%,N50長度為1 741 bp;組裝后BUSCO值C:94.5%[S:78.9%;D:15.6%],說明轉(zhuǎn)錄組組裝較完整;unigene平均長度為981.88 bp,長于中華蜜蜂(Apis cerana)(平均unigene 898 bp)[27]和澤蘭實蠅(Procecidochares utilis)(平均unigene 457 bp)[28],但是短于椰心葉甲嚙小蜂(Tetrastichus brontispae)(平均unigene 1 737.95 bp)[29]。有研究表明,序列長度是影響注釋的重要因素[30],本次測序獲得的N50、序列長度和堿基質(zhì)量比值等均顯示克什米爾熊蜂轉(zhuǎn)錄組的拼接注釋效果良好。

    本研究中,克什米爾熊蜂的GO功能注釋主要集中在結(jié)合、催化活性、細胞部分、細胞進程和代謝進程等功能中,與夜蛾黑卵蜂[31]、綠眼賽繭蜂[32]和椰心葉甲嚙小蜂[29]等膜翅目昆蟲轉(zhuǎn)錄組功能注釋結(jié)果相似。夜蛾黑卵蜂unigenes的GO功能注釋主要參與結(jié)合(4 925條unigenes)、催化活性(3 542)、細胞進程(5 124)、細胞解剖實體(4 314)和代謝進程(4 479)等功能[31];綠眼賽繭蜂unigenes的GO功能注釋主要參與細胞進程(5 230)、細胞部分(4 591)、代謝過程(4 068)、結(jié)合(6 652)和催化活性(5 302)等功能[32];椰心葉甲嚙小蜂unigenes的GO功能注釋主要參與細胞過程(1 527)、結(jié)合(1 733)和膜(1 324)等功能[29],這表明克什米爾熊蜂與這幾種蜂類的基因功能類似??耸裁谞栃芊涞腒EGG通路富集主要在信號轉(zhuǎn)導、碳水化合物代謝、氨基酸代謝、翻譯、運輸與分解代謝和內(nèi)分泌系統(tǒng)等過程中,富集到的主要通路有核糖體、戊糖和葡萄糖醛酸的相互轉(zhuǎn)化、內(nèi)質(zhì)網(wǎng)中的蛋白質(zhì)加工、氧化磷酸化和產(chǎn)熱等。其中富集最多的是與核糖體有關的基因,表明這些基因與克什米爾熊蜂蛋白的合成有關。

    3.2 克什米爾熊蜂SSR位點發(fā)生頻率及分布

    本研究中,從克什米爾熊蜂的轉(zhuǎn)錄組中共識別SSR位點27 414個,以單堿基、二堿基和三堿基核苷酸重復為主。SSR位點的發(fā)生頻率為17.81%,略高于意大利蜜蜂(Apis mellifera)(17.73%)[33],略低于中蜂(Apis cerana)(17.82%)[27],遠高于黑翅土白蟻 (Odontotermes formosaanus) (9.98%)[34]和煙粉虱 (Bemisia tabaci) (5.07%)[35]等昆蟲,表明其轉(zhuǎn)錄組中SSR相對比較豐富。不同昆蟲SSR位點出現(xiàn)頻率的差異與物種基因組特點有關,也與搜尋SSR位點時MISA軟件的參數(shù)設置有關[36]。如搜尋意大利蜜蜂[33]和中蜂[27]的SSR時,參數(shù)設置為二堿基至六堿基的最小重復為6,5,5,5,5,本次搜尋SSR位點與它們的參數(shù)一致,說明研究結(jié)果可靠。

    克什米爾熊蜂SSR單核苷酸重復中共有(A/T)n和(C/G)n兩種,A/T為優(yōu)勢基元,這與梨小食心蟲[37]和窄足真蚋[38]中優(yōu)勢堿基情況相同。且有研究表明,單堿基A/T既是單核苷酸重復中的優(yōu)勢堿基,也是所有重復類型的優(yōu)勢堿基[39]。克什米爾熊蜂SSR二堿基重復共有(AC/GT)n,(AG/CT)n,(AT/AT)n,(CG/CG)n四種,AT/AT為優(yōu)勢基元,與中華蜜蜂(Apis cerana cerana)[27]和意大利蜜蜂[33]相同,與青海草原毛蟲(Gynaephora qinghaiensis)(AC/GT)[40]不同??耸裁谞栃芊淙龎A基重復的優(yōu)勢基元為AAG/CTT,相比于單堿基和二堿基的重復基元出現(xiàn)的次數(shù),三堿基重復基元出現(xiàn)的次數(shù)更少,隨著重復基元堿基數(shù)的增加,重復基元的類型越來越多,但是重復基元出現(xiàn)的次數(shù)越來越少。另外有研究表明,當SSR重復次數(shù)高于12次時,多態(tài)性位點的比例較大[41]??耸裁谞栃芊滢D(zhuǎn)錄組SSR單堿基至六堿基的重復次數(shù)大多集中在11~15次(48.61%),說明其轉(zhuǎn)錄組SSR具有較高的多態(tài)性位點潛能。

    4 結(jié)論

    克什米爾熊蜂轉(zhuǎn)錄組的測序質(zhì)量良好,組裝注釋完整,其中GO功能注釋和KEGG通路富集主要集中于分子功能、細胞進程和碳水化合物代謝等一些途徑中?;谵D(zhuǎn)錄組數(shù)據(jù)的SSR位點出現(xiàn)頻率高、基元特征明顯。本研究通過對克什米爾熊蜂轉(zhuǎn)錄組測序數(shù)據(jù)的分析和微衛(wèi)星位點的篩選,豐富了克什米爾熊蜂的分子生物學研究,并設計開發(fā)出了17對能夠穩(wěn)定擴增的SSR引物,為進一步開展克什米爾熊蜂的遺傳學研究奠定了基礎。

    參考文獻

    [1] 梁程博,王久利,孫國,等. 三江源地區(qū)熊蜂物種多樣性研究[J]. 草地學報,2022,30(8):2126-2134

    [2] 周峰,姚麗媛,石涵,等. 傳粉熊蜂訪花行為的研究進展[J]. 昆蟲學報,2023,66(3):419-438

    [3] 黃家興,安建東. 中國熊蜂多樣性、人工利用與保護策略[J]. 生物多樣性,2018,26(5):486-497

    [4] 黃家興,安建東,吳杰,等. 熊蜂為溫室茄屬作物授粉的優(yōu)越性[J]. 中國農(nóng)學通報,2007(3):5-9

    [5] 黃訓兵,李輝,陳浩,等. 熊蜂授粉對北方冬季溫室不同品種番茄品質(zhì)的影響[J]. 北方園藝,2022(4):56-61

    [6] CAMERON S A,SADD B M. Global trends in bumble bee health[J]. Annual Review of Entomology,2020,65:209-232

    [7] SUN C,HUANG J,WANG Y,et al. Genus-wide characterization of bumblebee genomes provides insights into their evolution and variation in ecological and behavioral traits[J]. Molecular Biology and Evolution,2021,38(2):486-501

    [8] KWONG W K,ENGEL P,KOCH H,et al. Genomics and host specialization of honey bee and bumble bee gut symbionts[J]. Proceedings of the National Academy of Sciences,2014,111(31):11509-11514

    [9] SOROYE P,NEWBOLD T,KERR J. Climate change contributes to widespread declines among bumblebees across continents[J]. Science,2020,367(6478):685-688

    [10]MILLER-STRUTTMANN N E,GEIB J C,F(xiàn)RANKLIN J D,et al. Functional mismatch in a bumble bee pollination mutualism under climate change[J]. Science,2015,349(6255):1541-1544

    [11]SPRAYBERRY J D,RITTER K A,RIFFELL J A. The effect of olfactory exposure to non-insecticidal agrochemicals on bumblebee foraging behavior[J]. Plos One,2013,8(10):e76273

    [12]孫國,梁程博,趙芳,等. 15種熊蜂全基因組SSR分布規(guī)律研究[J]. 南方農(nóng)業(yè)學報,2023,54(11):3378-3387

    [13]WANG L,LIU S,TANG Y,et al. Using the combined gene approach and multiple analytical methods to improve the phylogeny and classification of Bombus (Hymenoptera,Apidae) in China[J]. Zookeys,2020,1007:1

    [14]WILLIAMS P H,BROWN M J,CAROLAN J C,et al. Unveiling cryptic species of the bumblebee subgenus Bombus s. str. worldwide with COI barcodes (Hymenoptera:Apidae)[J]. Systematics and Biodiversity,2012,10(1):21-56

    [15]LIN G,HUANG Z,WANG L,et al. Evolutionary rates of bumblebee genomes are faster at lower elevations[J]. Molecular Biology and Evolution,2019,36(6):1215-1219

    [16]SADD B M,BARRIBEAU S M,BLOCH G,et al. The genomes of two key bumblebee species with primitive eusocial organization[J]. Genome biology,2015,16(1):1-32

    [17]毛軒睿,劉玉萍,蘇旭,等. 沙鞭轉(zhuǎn)錄組簡單重復序列(SSR)位點特征分析[J]. 草地學報,2022,30(8):1990-2001

    [18]甘麗萍,田輝,唐恒,等. 6種鱗翅目昆蟲全基因組SSR分布規(guī)律[J]. 基因組學與應用生物學,2021,40(3):1022-1030

    [19]王日芳,李健,李吉濤,等. 脊尾白蝦(Exopalaemon carinicauda)近交系3個家系33個微衛(wèi)星座位的遺傳分析[J]. 漁業(yè)科學進展,2017,38(4):78-86

    [20]ZHOU L,LIU T,CHENG Y,et al. Molecular mapping of a stripe rust resistance gene in Chinese wheat landrace “Hejiangyizai” using SSR,RGAP,TRAP,and SRAP markers[J]. Crop Protection,2017,94:178-184

    [21]李東宇,孔杰,孟憲紅,等. 凡納濱對蝦(Litopenaeus vannamei)微衛(wèi)星多重PCR體系的建立及其在家系親權鑒定中的應用[J]. 漁業(yè)科學進展,2016,37(3):58-67

    [22]TAKEUCHI T,TAKAHASHI M,NISHIMOTO M,et al. Genetic structure of the bumblebee Bombus hypocrita sapporoensis,a potential domestic pollinator for crops in Japan[J]. Journal of Apicultural Research,2018,57(2):203-212

    [23]CEJAS D,DE LA RU′A P,ORNOSA C,et al. Spatial and temporal patterns of genetic diversity in Bombus terrestris populations of the Iberian Peninsula and their conservation implications[J]. Scientific Reports,2021,11(1):22471

    [24]JIANG L,SCHLESINGER F,DAVIS C A,et al. Synthetic spike-in standards for RNA-seq experiments[J]. Genome research,2011,21(9):1543-1551

    [25]郁偉杰,苗一凡,穆海婷,等. 轉(zhuǎn)錄組學解析紫花苜蓿雄性不育系花粉敗育機制[J]. 草地學報,2023,31(6):1702-1713

    [26]WANG Y,ZENG X,IYER N J,et al. Exploring the switchgrass transcriptome using second-generation sequencing technology[J]. Plos One,2012,7(3):e34225

    [27]熊翠玲,張璐,付中民,等. 基于RNA-seq數(shù)據(jù)大規(guī)模開發(fā)中華蜜蜂幼蟲的SSR分子標記[J]. 環(huán)境昆蟲學報,2017,39(1):68-74

    [28]高平,李麗芳,蘭明先,等. 澤蘭實蠅雌雄成蟲轉(zhuǎn)錄組及差異分析[J]. 應用昆蟲學報,2021,58(2):289-299

    [29]劉華偉,李朝緒,李芬,等. 椰心葉甲嚙小蜂轉(zhuǎn)錄組分析及基因功能注釋[J]. 中國生物防治學報,2021,37(3):412-419

    [30]ZHANG L,YANG X,QI X,et al. Characterizing the transcriptome and microsatellite markers for almond (Amygdalus communis L.) using the Illumina sequencing platform[J]. Hereditas,2018,155:1-9

    [31]周婉,徐琢,楊文俊,等. 夜蛾黑卵蜂轉(zhuǎn)錄組及化學感受相關基因的分析[J]. 植物保護,2022,48(6):264-277

    [32]王予彤,徐林波,段立清,等. 綠眼賽繭蜂轉(zhuǎn)錄組分析及嗅覺相關蛋白基因的鑒定[J]. 應用昆蟲學報,2021,58(4):846-855

    [33]郭睿,陳華枝,莊天藝,等. 利用轉(zhuǎn)錄組數(shù)據(jù)開發(fā)意大利蜜蜂的SSR分子標記[J]. 安徽農(nóng)業(yè)大學學報,2018,45(3):404-408

    [34]HUANG Q,SUN P,ZHOU X,et al. Characterization of head transcriptome and analysis of gene expression involved in caste differentiation and aggression in Odontotermes formosanus (Shiraki)[J]. Plos One,2012,7(11):e50383

    [35]XIE W,MENG Q,WU Q,et al. Pyrosequencing the Bemisia tabaci transcriptome reveals a highly diverse bacterial community and a robust system for insecticide resistance[J]. Plos One,2012,7(4):e35181

    [36]BISWAS M K,CHAI L,MAYER C,et al. Exploiting BAC-end sequences for the mining,characterization and utility of new short sequences repeat (SSR) markers in Citrus[J]. Molecular Biology Reports,2012,39:5373-5386

    [37]冷春蒙,李引,胡迪,等. 梨小食心蟲幼蟲中腸轉(zhuǎn)錄組及SSR分子標記分析[J]. 昆蟲學報,2018,61(11):1272-1283

    [38]郭歡,王剛,張樹田,等. 基于RNA-seq數(shù)據(jù)的窄足真蚋SSR分子標記開發(fā)[J]. 昆蟲學報,2018,61(7):815-824

    [39]王暉,王敬,高妍夏,等. 基于家蠶轉(zhuǎn)錄組測序的SSR序列分析[J]. 山東農(nóng)業(yè)科學,2022,54(1):14-20

    [40]南彥斌,許嘉誠,何啟玥,等. 青海草原毛蟲轉(zhuǎn)錄組分析及SSR位點開發(fā)[J]. 草地學報,2023,31(9):2653-2662

    [41]THAO D,YAMASHITA M,WATANABE A,et al. Development of tetranucleotide microsatellite markers in Pinus kesiya Royle ex Gordon[J]. Conservation Genetics Resources,2013,5:405-407

    (責任編輯 閔芝智)

    tube8黄色片| 亚洲av欧美aⅴ国产| 在线观看一区二区三区激情| 国产精品影院久久| 美女视频免费永久观看网站| 国产精品亚洲一级av第二区| 美国免费a级毛片| 99精品在免费线老司机午夜| 久久 成人 亚洲| 在线观看舔阴道视频| 久久精品国产综合久久久| 50天的宝宝边吃奶边哭怎么回事| 午夜成年电影在线免费观看| 亚洲人成电影观看| 成人18禁高潮啪啪吃奶动态图| 日韩欧美在线二视频 | 99精品在免费线老司机午夜| 黄色丝袜av网址大全| 黑人欧美特级aaaaaa片| av一本久久久久| 免费久久久久久久精品成人欧美视频| 欧美黄色淫秽网站| 婷婷丁香在线五月| 国产精品免费大片| 久久人妻av系列| 久久青草综合色| 黑人操中国人逼视频| 久久人妻av系列| 搡老熟女国产l中国老女人| 91九色精品人成在线观看| 欧美乱码精品一区二区三区| 久久国产亚洲av麻豆专区| 人人澡人人妻人| 精品少妇久久久久久888优播| 啦啦啦视频在线资源免费观看| 少妇猛男粗大的猛烈进出视频| 看黄色毛片网站| 亚洲欧洲精品一区二区精品久久久| 国产高清videossex| 正在播放国产对白刺激| 91大片在线观看| 国产精品永久免费网站| 色尼玛亚洲综合影院| 久久午夜综合久久蜜桃| 91大片在线观看| 精品国内亚洲2022精品成人 | 91麻豆av在线| 国产麻豆69| e午夜精品久久久久久久| bbb黄色大片| 热99re8久久精品国产| 国产欧美日韩综合在线一区二区| 国产黄色免费在线视频| 桃红色精品国产亚洲av| 男人的好看免费观看在线视频 | 国产乱人伦免费视频| 国产精品国产高清国产av | 国产精品二区激情视频| 岛国毛片在线播放| 亚洲国产精品一区二区三区在线| 精品福利永久在线观看| 电影成人av| 男人的好看免费观看在线视频 | 不卡av一区二区三区| 欧美精品啪啪一区二区三区| 精品福利观看| 黄色成人免费大全| 超碰97精品在线观看| 丁香欧美五月| avwww免费| 精品乱码久久久久久99久播| 99国产精品免费福利视频| 大香蕉久久网| 捣出白浆h1v1| 亚洲av日韩精品久久久久久密| 狠狠狠狠99中文字幕| 亚洲熟女精品中文字幕| e午夜精品久久久久久久| 夜夜夜夜夜久久久久| 法律面前人人平等表现在哪些方面| 精品少妇久久久久久888优播| 一本综合久久免费| 在线观看免费视频日本深夜| 午夜视频精品福利| www.自偷自拍.com| 自线自在国产av| 亚洲成人免费av在线播放| 久久国产精品男人的天堂亚洲| 天堂动漫精品| 一夜夜www| 99精品久久久久人妻精品| 亚洲综合色网址| 黑人欧美特级aaaaaa片| 亚洲午夜理论影院| 侵犯人妻中文字幕一二三四区| 老汉色∧v一级毛片| 日韩欧美一区视频在线观看| av视频免费观看在线观看| 国产免费现黄频在线看| 久久久久久久精品吃奶| 国产精品 国内视频| 国产欧美日韩综合在线一区二区| 老司机福利观看| 韩国av一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 精品一区二区三区视频在线观看免费 | 人人妻,人人澡人人爽秒播| 国产蜜桃级精品一区二区三区 | 成人18禁在线播放| 欧美激情 高清一区二区三区| 视频区欧美日本亚洲| 成人免费观看视频高清| 中国美女看黄片| 成人国语在线视频| 人人澡人人妻人| 久久精品亚洲熟妇少妇任你| 妹子高潮喷水视频| 亚洲国产毛片av蜜桃av| 亚洲成av片中文字幕在线观看| 国内毛片毛片毛片毛片毛片| 久久天躁狠狠躁夜夜2o2o| 一进一出抽搐动态| 18禁黄网站禁片午夜丰满| 夜夜夜夜夜久久久久| 午夜久久久在线观看| 脱女人内裤的视频| 欧美精品啪啪一区二区三区| 不卡一级毛片| 日本五十路高清| 精品乱码久久久久久99久播| 亚洲性夜色夜夜综合| 久久久久久久精品吃奶| 99香蕉大伊视频| 国产男女内射视频| 日日夜夜操网爽| 侵犯人妻中文字幕一二三四区| 老司机靠b影院| 亚洲国产看品久久| 99re6热这里在线精品视频| 久久人妻熟女aⅴ| 这个男人来自地球电影免费观看| 免费在线观看亚洲国产| 欧美一级毛片孕妇| av一本久久久久| 看片在线看免费视频| av不卡在线播放| 亚洲欧美激情在线| 久久国产乱子伦精品免费另类| 日本精品一区二区三区蜜桃| 黄色视频,在线免费观看| 高清毛片免费观看视频网站 | 久久久国产欧美日韩av| 精品国产美女av久久久久小说| 亚洲av日韩精品久久久久久密| 久久天堂一区二区三区四区| 成人亚洲精品一区在线观看| 亚洲国产精品一区二区三区在线| 国产成人av教育| 伦理电影免费视频| 午夜福利影视在线免费观看| 丰满饥渴人妻一区二区三| 亚洲国产毛片av蜜桃av| 老汉色av国产亚洲站长工具| 国产成人精品无人区| 欧美日韩一级在线毛片| 欧美国产精品一级二级三级| 国产在视频线精品| 午夜精品国产一区二区电影| 亚洲久久久国产精品| 一边摸一边抽搐一进一小说 | www.自偷自拍.com| 亚洲片人在线观看| 国产三级黄色录像| 久久性视频一级片| 美女扒开内裤让男人捅视频| 99国产精品99久久久久| 悠悠久久av| 国产精品久久电影中文字幕 | 亚洲欧美精品综合一区二区三区| 99国产极品粉嫩在线观看| 十八禁网站免费在线| 日韩人妻精品一区2区三区| 99久久精品国产亚洲精品| 久久久久久久久免费视频了| 免费在线观看日本一区| 99国产精品一区二区三区| www.自偷自拍.com| 欧美激情极品国产一区二区三区| a级毛片在线看网站| 女人精品久久久久毛片| 色播在线永久视频| 国产精品一区二区精品视频观看| 国产深夜福利视频在线观看| 久久中文看片网| 狠狠狠狠99中文字幕| 久99久视频精品免费| 国产精品一区二区免费欧美| 777久久人妻少妇嫩草av网站| 亚洲精品在线观看二区| 欧美激情 高清一区二区三区| 国产伦人伦偷精品视频| 国产欧美日韩综合在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 亚洲专区字幕在线| 成人国产一区最新在线观看| 欧美黑人精品巨大| 性色av乱码一区二区三区2| 制服诱惑二区| √禁漫天堂资源中文www| 亚洲av片天天在线观看| 国产亚洲精品久久久久久毛片 | 天天操日日干夜夜撸| 91麻豆av在线| 一本综合久久免费| 老熟女久久久| 国产精品九九99| 久9热在线精品视频| 亚洲在线自拍视频| 少妇 在线观看| 80岁老熟妇乱子伦牲交| 国产精品久久视频播放| 老司机午夜福利在线观看视频| 99国产精品一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲熟妇中文字幕五十中出 | 青草久久国产| 91大片在线观看| 婷婷成人精品国产| 大香蕉久久网| 极品人妻少妇av视频| 91九色精品人成在线观看| 99国产极品粉嫩在线观看| 女性生殖器流出的白浆| 美女 人体艺术 gogo| 国产精品久久久久成人av| 中文字幕制服av| 在线天堂中文资源库| 国产欧美日韩一区二区精品| 丁香六月欧美| 日韩制服丝袜自拍偷拍| 欧美日韩中文字幕国产精品一区二区三区 | 日本一区二区免费在线视频| 久久国产亚洲av麻豆专区| 飞空精品影院首页| 亚洲精品国产精品久久久不卡| 极品人妻少妇av视频| 精品国产国语对白av| 婷婷精品国产亚洲av在线 | 91在线观看av| 日韩欧美国产一区二区入口| 国产精品自产拍在线观看55亚洲 | 下体分泌物呈黄色| 一a级毛片在线观看| 日韩视频一区二区在线观看| 免费人成视频x8x8入口观看| 国产精品永久免费网站| 曰老女人黄片| 精品人妻在线不人妻| 99riav亚洲国产免费| 欧洲精品卡2卡3卡4卡5卡区| 97人妻天天添夜夜摸| 国产91精品成人一区二区三区| 最新的欧美精品一区二区| 日日爽夜夜爽网站| av中文乱码字幕在线| 18禁美女被吸乳视频| 免费久久久久久久精品成人欧美视频| 欧美日韩精品网址| 国产不卡av网站在线观看| 69av精品久久久久久| 亚洲情色 制服丝袜| 丝袜在线中文字幕| 国产一区二区三区在线臀色熟女 | 欧美另类亚洲清纯唯美| 国产精品99久久99久久久不卡| 热99久久久久精品小说推荐| 麻豆国产av国片精品| 国产激情欧美一区二区| 色在线成人网| 亚洲国产精品一区二区三区在线| 久久久久精品国产欧美久久久| 黄色女人牲交| 亚洲一区二区三区欧美精品| 国产精品永久免费网站| 999精品在线视频| 久久久国产精品麻豆| 色综合欧美亚洲国产小说| 亚洲五月色婷婷综合| 亚洲精品美女久久av网站| 性少妇av在线| 9热在线视频观看99| 日韩欧美三级三区| 精品一区二区三区四区五区乱码| 欧美性长视频在线观看| 他把我摸到了高潮在线观看| 日韩有码中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 热99久久久久精品小说推荐| 18禁黄网站禁片午夜丰满| 国产精品久久久久成人av| 亚洲第一青青草原| 99热网站在线观看| 亚洲欧美激情在线| 成人亚洲精品一区在线观看| 久久国产精品人妻蜜桃| 久久中文看片网| 国产精品 欧美亚洲| 国产一区有黄有色的免费视频| 精品国产一区二区久久| 日韩制服丝袜自拍偷拍| 精品无人区乱码1区二区| 不卡一级毛片| 成人黄色视频免费在线看| 国产精品亚洲一级av第二区| 亚洲精华国产精华精| 国产精品乱码一区二三区的特点 | 人妻 亚洲 视频| 欧美黄色淫秽网站| 亚洲国产毛片av蜜桃av| 精品欧美一区二区三区在线| 亚洲第一av免费看| 国产亚洲欧美精品永久| 亚洲片人在线观看| 免费av中文字幕在线| 国产不卡一卡二| 久久久久视频综合| 午夜视频精品福利| 国产精品成人在线| 精品电影一区二区在线| 一进一出抽搐动态| av线在线观看网站| 操美女的视频在线观看| 新久久久久国产一级毛片| 热99re8久久精品国产| 欧美日韩亚洲高清精品| 久9热在线精品视频| 精品一区二区三区视频在线观看免费 | 91在线观看av| 久久中文字幕人妻熟女| 又黄又粗又硬又大视频| 午夜亚洲福利在线播放| 免费在线观看视频国产中文字幕亚洲| 99国产综合亚洲精品| 国产1区2区3区精品| 啦啦啦 在线观看视频| 精品少妇久久久久久888优播| 99热只有精品国产| 精品少妇久久久久久888优播| 黄片小视频在线播放| 国内久久婷婷六月综合欲色啪| avwww免费| 大片电影免费在线观看免费| 一区二区日韩欧美中文字幕| 成人影院久久| 国产欧美日韩一区二区三| 国产精品永久免费网站| 80岁老熟妇乱子伦牲交| 91大片在线观看| 18禁国产床啪视频网站| 夜夜躁狠狠躁天天躁| 18禁国产床啪视频网站| 我的亚洲天堂| 五月开心婷婷网| av国产精品久久久久影院| 天堂√8在线中文| xxxhd国产人妻xxx| 女性生殖器流出的白浆| 久久精品国产综合久久久| 亚洲人成77777在线视频| 99久久99久久久精品蜜桃| 免费高清在线观看日韩| 国产免费男女视频| 大型av网站在线播放| 一边摸一边抽搐一进一出视频| 亚洲aⅴ乱码一区二区在线播放 | 黄片小视频在线播放| 成人18禁高潮啪啪吃奶动态图| 欧美老熟妇乱子伦牲交| 免费观看精品视频网站| 欧美日韩中文字幕国产精品一区二区三区 | 国产一区二区三区在线臀色熟女 | 日韩成人在线观看一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 99re在线观看精品视频| 国产片内射在线| 欧美老熟妇乱子伦牲交| 国产99久久九九免费精品| 黄色视频不卡| 久久精品国产综合久久久| 亚洲九九香蕉| a级毛片黄视频| 老鸭窝网址在线观看| 久热爱精品视频在线9| 亚洲人成电影观看| 欧美日韩亚洲高清精品| 日本撒尿小便嘘嘘汇集6| 国产野战对白在线观看| 91麻豆av在线| 一进一出好大好爽视频| xxx96com| 日韩熟女老妇一区二区性免费视频| 看片在线看免费视频| 人人妻人人爽人人添夜夜欢视频| 亚洲avbb在线观看| 精品一品国产午夜福利视频| 91大片在线观看| 巨乳人妻的诱惑在线观看| 女警被强在线播放| 精品国产亚洲在线| 18禁黄网站禁片午夜丰满| av天堂在线播放| 国产单亲对白刺激| 两个人免费观看高清视频| 搡老乐熟女国产| 天堂动漫精品| 精品国产国语对白av| 香蕉丝袜av| 色老头精品视频在线观看| 久久天躁狠狠躁夜夜2o2o| 成人特级黄色片久久久久久久| 亚洲色图 男人天堂 中文字幕| 久久久久久久国产电影| 久久久国产成人精品二区 | a级毛片在线看网站| 国产亚洲精品久久久久久毛片 | 在线视频色国产色| 国产精品永久免费网站| 亚洲av成人av| 天天躁日日躁夜夜躁夜夜| 亚洲精品中文字幕一二三四区| 亚洲一区中文字幕在线| 91国产中文字幕| 纯流量卡能插随身wifi吗| 日本撒尿小便嘘嘘汇集6| 亚洲第一av免费看| 制服人妻中文乱码| 久久久久精品人妻al黑| 色播在线永久视频| 男女下面插进去视频免费观看| 亚洲五月色婷婷综合| 国产欧美日韩一区二区精品| 极品少妇高潮喷水抽搐| 日韩欧美在线二视频 | 韩国av一区二区三区四区| 男女免费视频国产| 天堂俺去俺来也www色官网| 亚洲人成电影免费在线| 国产一区有黄有色的免费视频| 多毛熟女@视频| 99国产精品一区二区蜜桃av | 视频区图区小说| 无遮挡黄片免费观看| 夫妻午夜视频| 亚洲精品美女久久av网站| 久久精品国产综合久久久| 99国产精品一区二区蜜桃av | 性少妇av在线| 国产高清激情床上av| 我的亚洲天堂| 欧美日韩亚洲综合一区二区三区_| 欧美不卡视频在线免费观看 | 国产一区有黄有色的免费视频| 欧美日韩亚洲国产一区二区在线观看 | 国产精品久久视频播放| 老熟妇仑乱视频hdxx| 很黄的视频免费| 国产又色又爽无遮挡免费看| 我的亚洲天堂| 久久中文字幕一级| 亚洲五月色婷婷综合| 中文字幕人妻熟女乱码| 一级a爱片免费观看的视频| 国产人伦9x9x在线观看| 90打野战视频偷拍视频| 午夜激情av网站| 一区在线观看完整版| 午夜影院日韩av| 丰满人妻熟妇乱又伦精品不卡| 亚洲熟妇熟女久久| 热99国产精品久久久久久7| 日韩大码丰满熟妇| 伊人久久大香线蕉亚洲五| 老司机福利观看| 一本一本久久a久久精品综合妖精| 中文字幕高清在线视频| av网站免费在线观看视频| 欧美亚洲日本最大视频资源| 欧美丝袜亚洲另类 | 搡老岳熟女国产| 91大片在线观看| 露出奶头的视频| 女性被躁到高潮视频| 成人手机av| 极品少妇高潮喷水抽搐| 成人影院久久| 久久国产精品大桥未久av| 777久久人妻少妇嫩草av网站| 久久久国产成人免费| 国产成人系列免费观看| 91老司机精品| 久久人妻熟女aⅴ| 天堂俺去俺来也www色官网| 极品人妻少妇av视频| 久热这里只有精品99| 日日夜夜操网爽| 男女免费视频国产| 免费在线观看影片大全网站| 精品午夜福利视频在线观看一区| 日韩欧美三级三区| 欧美乱色亚洲激情| 两性夫妻黄色片| 色在线成人网| 精品久久久久久久久久免费视频 | 一边摸一边抽搐一进一小说 | 国产精品 欧美亚洲| 岛国毛片在线播放| 欧美不卡视频在线免费观看 | 色综合婷婷激情| av超薄肉色丝袜交足视频| 999久久久国产精品视频| 高清av免费在线| 日韩三级视频一区二区三区| 国产精品乱码一区二三区的特点 | 亚洲一区二区三区欧美精品| 99久久人妻综合| 黄色丝袜av网址大全| 色94色欧美一区二区| 国产精品1区2区在线观看. | 老汉色∧v一级毛片| 老司机午夜福利在线观看视频| 欧美黑人欧美精品刺激| 99国产精品一区二区三区| 欧美成狂野欧美在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 夜夜躁狠狠躁天天躁| 天堂中文最新版在线下载| 国产一区二区三区视频了| 成人国语在线视频| 亚洲成av片中文字幕在线观看| 中文字幕精品免费在线观看视频| 欧美精品一区二区免费开放| 亚洲欧美色中文字幕在线| 精品国产乱子伦一区二区三区| а√天堂www在线а√下载 | 亚洲综合色网址| 天堂动漫精品| 又黄又爽又免费观看的视频| 国产欧美日韩综合在线一区二区| 人妻 亚洲 视频| 亚洲午夜理论影院| 在线观看免费日韩欧美大片| 国产真人三级小视频在线观看| 成人一区二区视频在线观看| 夜夜爽天天搞| 99国产极品粉嫩在线观看| 国产综合懂色| 性色av乱码一区二区三区2| 麻豆久久精品国产亚洲av| 亚洲av不卡在线观看| 天堂√8在线中文| 亚洲av五月六月丁香网| 中文字幕精品亚洲无线码一区| 久久久久久久久大av| 成人永久免费在线观看视频| 国产亚洲av嫩草精品影院| 啦啦啦免费观看视频1| av专区在线播放| 精品久久久久久久末码| 日本与韩国留学比较| 夜夜躁狠狠躁天天躁| 欧美在线一区亚洲| 91字幕亚洲| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲av熟女| 久久精品人妻少妇| 最新美女视频免费是黄的| 搡老妇女老女人老熟妇| 午夜福利高清视频| 日韩欧美国产在线观看| 男女之事视频高清在线观看| ponron亚洲| 欧美国产日韩亚洲一区| 免费大片18禁| 一区二区三区高清视频在线| 国产av一区在线观看免费| 亚洲av一区综合| 久久久久久人人人人人| 高清在线国产一区| 观看美女的网站| 日本免费一区二区三区高清不卡| av黄色大香蕉| 日本a在线网址| 伊人久久大香线蕉亚洲五| 久久人妻av系列| 欧美激情在线99| 午夜福利在线在线| 人妻丰满熟妇av一区二区三区| 夜夜爽天天搞| 午夜影院日韩av| 成人欧美大片| 俺也久久电影网| 国产不卡一卡二| 91在线观看av| a在线观看视频网站| 精华霜和精华液先用哪个| 久久久久久久久中文| 亚洲最大成人手机在线| eeuss影院久久| 久久精品国产综合久久久| 好男人在线观看高清免费视频| 男女做爰动态图高潮gif福利片| 欧美高清成人免费视频www| 午夜精品在线福利| 九九久久精品国产亚洲av麻豆| 欧美绝顶高潮抽搐喷水|