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

    脊尾白蝦血細胞ESTs的生物信息學與微衛(wèi)星序列特征分析

    2016-12-19 06:41:59段亞飛李吉濤
    水產(chǎn)科學 2016年5期
    關(guān)鍵詞:脊尾白蝦微衛(wèi)星

    段亞飛, 張 喆, 李吉濤, 李 健, 劉 萍

    ( 1.中國水產(chǎn)科學研究院 南海水產(chǎn)研究所,農(nóng)業(yè)部南海漁業(yè)資源開發(fā)利用重點實驗室,廣東 廣州 510300;2. 中國水產(chǎn)科學研究院 黃海水產(chǎn)研究所,農(nóng)業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點實驗室,山東 青島 266071;3. 青島海洋科學與技術(shù)國家實驗室,海洋漁業(yè)科學與食物產(chǎn)出過程功能實驗室,山東 青島 266071 )

    脊尾白蝦血細胞ESTs的生物信息學與微衛(wèi)星序列特征分析

    段亞飛1, 張 喆1, 李吉濤2, 3, 李 健2, 3, 劉 萍2, 3

    ( 1.中國水產(chǎn)科學研究院 南海水產(chǎn)研究所,農(nóng)業(yè)部南海漁業(yè)資源開發(fā)利用重點實驗室,廣東 廣州 510300;2. 中國水產(chǎn)科學研究院 黃海水產(chǎn)研究所,農(nóng)業(yè)部海洋漁業(yè)可持續(xù)發(fā)展重點實驗室,山東 青島 266071;3. 青島海洋科學與技術(shù)國家實驗室,海洋漁業(yè)科學與食物產(chǎn)出過程功能實驗室,山東 青島 266071 )

    對前期測序得到的脊尾白蝦血細胞2853條EST序列進行了生物信息學和微衛(wèi)星序列特征分析。EST序列拼接得到1053條Unigenes,包括329條Contigs和724條Singlets。BLAST分析表明,593 (56.3%)條Uingenes與數(shù)據(jù)庫中已知基因具有相似性。KEGG代謝途徑分析表明,181條Unigenes映射到120條代謝途徑。通過EST-SSR分析,共得到416條微衛(wèi)星序列,檢出率為14.58%。其中,兩堿基重復序列374條,占89.90%,AG重復類型最多;三堿基35條,占8.41%,AAT重復類型最多;四堿基7條,占1.68%,AAGT重復類型最多。本研究可為脊尾白蝦功能基因資源挖掘及分子標記篩選提供有效數(shù)據(jù)。

    脊尾白蝦;血細胞;EST;生物信息學;微衛(wèi)星

    脊尾白蝦(Exopalaemoncarinicauda)是我國黃、渤海地區(qū)重要的經(jīng)濟蝦類,具有繁殖能力強、生長速度快和環(huán)境適應性強等優(yōu)點[1-2]。近年來,脊尾白蝦養(yǎng)殖產(chǎn)業(yè)發(fā)展迅速,其遺傳育種研究日益受到重視。脊尾白蝦基因信息和分子遺傳標記的發(fā)掘,可以為其功能基因篩選、遺傳多樣性分析、遺傳連鎖圖譜構(gòu)建及標記輔助育種等奠定基礎(chǔ)。

    表達序列標簽(EST)是指從特定組織來源的cDNA文庫中隨機挑取單克隆進行單向測序而得到的部分cDNA序列,代表了生物體特定組織在某一時期的表達基因[3]。表達序列標簽分析可以用于新基因發(fā)現(xiàn)、功能基因篩選、分子標記開發(fā)和基因時空表達狀況分析等[3]。微衛(wèi)星標記(SSR)是水產(chǎn)動物遺傳育種重要的分子標記,具有多態(tài)性高、共顯性、通用性好、數(shù)量豐富和在基因組上分布均勻等特點[4-5]。EST-SSR可以避免傳統(tǒng)方法構(gòu)建、篩選文庫的繁瑣步驟和人力物力的耗費,廣泛應用于水生動物新基因的分離鑒定、比較基因組學、遺傳圖譜構(gòu)建和基因芯片制備等研究[6-7]。

    近年來,許多甲殼動物已獲得EST-SSRs標記,如中國明對蝦(Fenneropenaeuschinensis)[8]、凡納濱對蝦(Litopenaeusvannamei)[9]、斑節(jié)對蝦(Penaeusmonodon)[10]、三疣梭子蟹(Portunustrituberculatus)[11]和擬穴青蟹(Scyllaparamamosain)[12]等。而目前脊尾白蝦微衛(wèi)星標記僅見通過構(gòu)建基因組微衛(wèi)星富集文庫篩選獲得[13]。本研究利用前期構(gòu)建的脊尾白蝦血細胞全長cDNA文庫,對測序獲得的EST序列進行生物信息學和微衛(wèi)星序列特征分析,以期為脊尾白蝦功能基因資源挖掘及分子標記篩選提供有效數(shù)據(jù)。

    1 材料與方法

    1.1脊尾白蝦血細胞EST序列測定和分析

    利用前期構(gòu)建的脊尾白蝦血細胞全長cDNA文庫[3],對其進行批量測序,引物為5′pTriPLEx2 (5′-CTCCGAGATCTGGACGAGC-3′)。測得的EST序列,利用Cross-match軟件進行引物序列、載體序列及低質(zhì)量序列(<100 bp)去除;利用Phrap軟件對高質(zhì)量ESTs進行序列拼接,拼接結(jié)果為由Contigs和Singlets組成的Unigenes。

    1.2 ESTs的相似性比對、分類及代謝途徑分析

    利用BLAST軟件對拼接到的Unigenes與NCBI核苷酸和蛋白質(zhì)數(shù)據(jù)庫進行相似性比對(E-value < 1e-5)。將所有的Unigenes與日本京都大學的《京都基因和基因組百科全書》(KEGG)的PATHWAY子數(shù)據(jù)庫進行比較和代謝途徑分析,獲得脊尾白蝦血細胞的ESTs的部分代謝網(wǎng)絡(luò)信息。

    1.3 ESTs的微衛(wèi)星序列特征分析

    利用RepeatMasker軟件從拼接的Unigenes中查找微衛(wèi)星序列,最小重復參數(shù)設(shè)置為重復7次以上的兩堿基序列,重復5次以上的三堿基序列,重復5次以上的四堿基序列。單核苷酸重復序列未進行篩選,因為一般認為其作為多態(tài)性標記的意義不大[18]。

    2 結(jié) 果

    2.1 ESTs的生物信息學分析

    2.1.1 ESTs序列分析

    對脊尾白蝦血細胞全長cDNA文庫批量測序獲得的EST序列,去除載體序列、引物序列和低質(zhì)量載體序列(<100 bp)后,共得到2853條高質(zhì)量EST序列,平均長度為493 bp(表1)。利用Phrap軟件進行序列拼接,共獲得1053條Unigenes,包括329條(31.24%)Contigs和724條(68.76%)Singlets。對Unigenes中的ESTs數(shù)量分析可知,Contigs由2~284條ESTs拼接而成,其中大部分是由2~5條ESTs拼接而成,10條Contigs是由多于50條的ESTs拼接而成,只有1條Contigs是由284條ESTs拼接而成(圖1)。

    表1 脊尾白蝦血細胞ESTs基本信息

    BLAST分析表明,593(56.3%)條Uingenes與數(shù)據(jù)庫中已知基因具有相似性,22(2.1%)條Uingenes顯示為未知基因,438(41.6%)條Uingenes沒有任何相似性。ESTs相似比對的E-value分布表明,276(44.9%)條Uingenes具有高度相似性(E-value < 1e-100),256(41.6%)條Uingenes具有中度相似性(1e-99

    圖1 脊尾白蝦血細胞ESTs與其所拼接形成的Unigene之間的關(guān)系

    圖2 脊尾白蝦血細胞ESTs相似比對的E-value分布 圖3 脊尾白蝦血細胞ESTs匹配的相似物種分布

    2.1.2 ESTs的代謝途徑分析

    脊尾白蝦血細胞ESTs的代謝途徑分析表明,共有181條Unigenes得到注釋,相關(guān)酶類參與120條代謝途徑,主要包括代謝、細胞過程、環(huán)境信息的加工、遺傳信息的加工和人類疾病5類,其相關(guān)的基因序列分別占33.1%、21.6%、12.2%、16.0%和17.1%。其中,代謝途徑最多的類別是氨基酸代謝(11%),其次為傳染病(10%),第三為信號轉(zhuǎn)導(9%)(表2)。

    表2 脊尾白蝦血細胞ESTs的代謝途徑分析

    2.2 ESTs的微衛(wèi)星序列分析

    從血細胞2853條ESTs序列中,共篩選到微衛(wèi)星序列416個,占全部ESTs總數(shù)的14.58%。考慮到堿基互補配對和計數(shù)拷貝數(shù)起始堿基順序的差異,本研究將同類重復兼并為一種重復類型代表,如AAG代表六種重復類型,分別為AAG、AGA、GAA、TTC、TCT和CTT。統(tǒng)計所有重復類型,以兩堿基重復序列數(shù)量最多,為374條,占重復序列總數(shù)的89.90%;其次為三堿基,35個,占8.41%;再次是四堿基,7個,占1.68%(表3)。

    表3 不同重復類型的重復序列數(shù)目及其百分比

    同類型的SSR重復序列中,各重復拷貝類別所占的比例也各不相同。在兩堿基重復類型中,AG重復拷貝類別最多,為355個,占94.92%;其次為AC和AT,分別為12(3.21%)和7(1.87%);未發(fā)現(xiàn)GC重復拷貝類別。三堿基重復類型中,共有6種重復拷貝類別,分別為AAT、AAG、ATC、ACT、AGC和AGG。其中,以AAT最多,為21個,占60%;其次為ACT,為5個,占14.29%;再次為AAG和ATC,這兩種重復拷貝類別數(shù)量相同,為3個,占8.57%。四堿基重復類型中,共有2種重復拷貝類別,分別為AAGT和GCCG,其中以AAGT最多,為7個,占85.71%;其次為GCCG,為1個,占14.29%(表4)。

    表4脊尾白蝦血細胞ESTs篩選到的微衛(wèi)星序列

    重復類型重復單元數(shù)量比例/%二堿基AT71.87AG35594.92AC123.21三堿基AAT2160.00AAG38.57ATC38.57ACT514.29AGC12.86AGG25.71四堿基AAGT685.71GCCG114.29

    3 討 論

    3.1 ESTs生物信息學分析

    通過對前期測序得到的脊尾白蝦血細胞2853條ESTs進行序列拼接,共得到1053條Uingenes。BLAST分析表明,41.6%的Uingenes與美國國立生物技術(shù)信息中心數(shù)據(jù)庫中的已知基因序列無相似性,與三疣梭子蟹[14]、中華絨螯蟹(Eriocheirsinensis)[15]等的研究一致。此類無相似性的ESTs可能為新基因序列,表明甲殼類尚有大量基因未被發(fā)掘;又可能為已知基因,由于蝦類與其他脊椎動物的親緣關(guān)系太遠,而無法通過BLAST比對確定其相似關(guān)系。此外,Xu等[14]認為,ESTs是通過一次測序得到的序列,未必包含編碼區(qū),可能是ESTs未能注釋的原因之一。因此,此類未得到注釋的ESTs,有待進一步的深入研究。

    KEGG數(shù)據(jù)庫涵蓋了物種的基因與基因組、酶促反應及其途徑和各種生化物質(zhì)的信息,在基因與基因組、代謝途徑與代謝網(wǎng)絡(luò)的相關(guān)研究中具有重要作用[16]。KEGG PATHWAY數(shù)據(jù)庫不僅能夠提供生化物質(zhì)相互轉(zhuǎn)化所有可能的代謝途徑,還可以對催化各步反應的酶進行全面注解[16]。通過代謝途徑分析,本研究共有181條Unigenes(17.2%)獲得注釋,發(fā)現(xiàn)酶類參與120條代謝途徑。代謝途徑分析在多個水產(chǎn)物種中已有相關(guān)報道。Meng等[17]對蝦夷扇貝(Patinopectenyessoensis)肌肉EST序列進行KEGG代謝途徑分析,共有11.3%的Unigenes得到注釋,涉及103個代謝途徑;其中以代謝方面的Unigenes比例最大,占30.5%,與本研究結(jié)果一致。

    3.2 ESTs的微衛(wèi)星序列特征分析

    本研究通過對脊尾白蝦血細胞ESTs進行SSR篩選,共得到416條(14.58%)SSR序列,表明脊尾白蝦ESTs-SSR序列比例較高。研究表明,動物的ESTs-SSR以兩堿基重復為主,大部分植物則以三堿基重復為主[18]。脊尾白蝦血細胞ESTs-SSR共發(fā)現(xiàn)3種堿基重復類型,以兩堿基重復類型最多,占89.90%;其次為三堿基重復類型(8.41%)和四堿基重復類型(1.68%)。本研究未發(fā)現(xiàn)五堿基和六堿基重復類型,可能與所分析的ESTs序列的組織差異有關(guān)。

    兩堿基重復類型中,共發(fā)現(xiàn)3種重復拷貝類別,其中以AG最多,其次是AC,與斑節(jié)對蝦[19]、羅氏沼蝦(Macrobrachiumrosenbergii)[20]和三疣梭子蟹[21]等一致。賈舒雯等[13]通過構(gòu)建脊尾白蝦微衛(wèi)星富集文庫對多態(tài)性標記進行了篩選,發(fā)現(xiàn)脊尾白蝦基因組微衛(wèi)星核心序列以AG重復最多,占88.94%;其次是AC重復,占5.03%。這與本研究兩堿基重復類型結(jié)果(AG 94.92%,AC 3.21%)一致,由此說明AG重復拷貝類別在脊尾白蝦基因組中含量豐富。但本研究并未發(fā)現(xiàn)GC重復,而GC兩堿基重復在其他物種基因組中的含量均較少[22-23]。

    三堿基重復類型中,共發(fā)現(xiàn)6種重復拷貝類別,其中以AAT最多,占60%。研究表明,中國明對蝦[18]和紅鰭東方鲀(Takifugurubripes)[24]中,以AAT重復最多。高煥等[18]通過分析中國明對蝦基因組微衛(wèi)星特征,發(fā)現(xiàn)三堿基重復類型中以AAT重復最多,占50.34%。四堿基重復類型中,共發(fā)現(xiàn)7種重復拷貝類別,總體數(shù)量較少。而本研究中未發(fā)現(xiàn)五堿基和六堿基重復類型,可能與所分析的ESTs序列相對較少、覆蓋面較小有關(guān)。

    目前脊尾白蝦SSR的研究,主要是利用SSR核心序列進行探針雜交或以核心重復序列為引物進行PCR篩選,而從ESTs中篩選SSR序列的方法則更為簡單,不僅能夠節(jié)省人力和物力,而且由于ESTs-SSR標記位于基因編碼區(qū)中,因此具有高度的物種間轉(zhuǎn)移性[25]。本研究所篩選獲得的SSR序列,可以用于脊尾白蝦SSR標記篩選,對脊尾白蝦遺傳圖譜、基因定位和遺傳多樣性分析等奠定基礎(chǔ)。

    [1] 梁俊平, 李健, 劉萍, 等. 脊尾白蝦生物學特性與人工繁育的研究進展[J]. 中國農(nóng)學通報, 2012,28(17):109-116.

    [2] 梁俊平, 李健, 李吉濤, 等. 氨氮對脊尾白蝦幼蝦和成蝦的毒性試驗[J]. 水產(chǎn)科學, 2012,31(9):526-529.

    [3] 段亞飛, 劉萍, 李吉濤, 等. 脊尾白蝦血細胞全長cDNA文庫的構(gòu)建及EST序列分析[J]. 中國水產(chǎn)科學, 2013, 20(2):243-249.

    [4] 孫立元, 郭華陽, 朱彩艷, 等. 卵形鯧鲹育種群體遺傳多樣性分析[J]. 南方水產(chǎn)科學, 2014,10(2):67-71.

    [5] 楊兵, 林琳, 李純厚, 等. 基于高通量測序的二長棘鯛微衛(wèi)星標記開發(fā)與評價[J]. 南方水產(chǎn)科學, 2015,11(4):116-120.

    [6] Ji X S, Chen S L, Ma H Y, et al. Isolation and characterization of 19 EST-linked ploymorphic microsatellite loci for olive flounder (Paralichthysolivaceus)[J]. Aquaculture Research, 2009,40(8): 980-983.

    [7] Yu H, Li Q. Exploiting EST databases for the development and characterization of EST-SSRs in the Pacific oyster (Crassostreagigas)[J]. J Heredity,2008,99(2):208-214.

    [8] Wang H X, Li F H, Xiang J H. Polymorphic EST-SSR markers and their mode of inheritance inFenneropenaeuschinensis[J]. Aquaculture, 2005,249(1/4):107-114.

    [9] Perez F, Ortiz J,Zhinaula M, et al. Development of EST-SSR markers by data mining in three species of shrimp:Litopenaeusvannamei,Litopenaeusstylirostris, andTrachypenaeusbirdy[J]. Marine Biotechnol, 2005,7(5):554-569.

    [10] Maneeruttanarungroj C, Pongsomboon S, Wuthisuthimethavee S, et al. Development of polymorphic expressed sequence tag-derived microsatellites for the extension of the genetic linkage map of the black tiger shrimp (Penaeusmonodon) [J]. Animal Genetics, 2006,37(4):363-368.

    [11] 呂建建, 王渝, 高保全, 等. 三疣梭子蟹Ⅰ型微衛(wèi)星標記的發(fā)掘及多態(tài)性分析[J]. 水產(chǎn)學報, 2013,37(6):816-822.

    [12] Zhou Z H, Zhang Z P, Wang Y L, et al. EST analysis on the gonad development related organs and microarray screen for differentially expressed genes in mature ovary and testis ofScyllaparamamosain[J]. Comp Biochem and Physiol D, 2011,6(2):150-157.

    [13] 賈舒雯, 劉萍, 韓智科, 等.脊尾白蝦微衛(wèi)星富集文庫的構(gòu)建與多態(tài)性標記的篩選[J]. 水產(chǎn)學報, 2011,35(12):1787-1794.

    [14] Xu Q H, Liu Y, Liu R L. Expressed sequence tags from cDNA library prepared from gills of the swimming crab,Portunustrituberculatus[J]. J Experimental Marine Biol Ecol, 2010,394(1/2):105-115.

    [15] Gai Y C, Wang L L, Zhao J M, et al. The construction of a cDNA library enriched for immune genes and the analysis of 7535 ESTs from Chinese mitten crabEriocheirsinensis[J]. Fish and Shellfish Immunol, 2009,27(6):684-694.

    [16] Priebe S, Linde J, Albrecht D, et al. FungiFun:a web-based application for functional categorization of fungal genes and proteins [J]. Fungal Genetics Biol, 2011,48(4):353-358.

    [17] Meng X Y, Chang Y Q, Qiu X M, et al. Generation and analysis of expressed sequence tags from adductor muscle of Japanese scallopMizuhopectenyessoensis[J]. Comp Biochem Physiol D, 2010,5(4):288-294.

    [18] 高煥, 劉萍, 孟憲紅, 等. 中國對蝦基因組微衛(wèi)星特征分析[J]. 海洋與湖沼, 2004,35(5):424-431.

    [19] Xu Z,Dhar A K, Wyrzykowski J, et al. Identification of abundant and informative microsatellites from shrimp (Penaeusmonodon) genome [J]. Animal Genetics, 1999,30(2):150-156.

    [20] Charoentawee K, Poompuang S, Nakom U. Isolation and characterization of microsatellites in giant freshwater prawnMacrobrachiumrosenbergii[J]. Mol Ecol Resources, 2006,6(3):823-825.

    [21] 宋來鵬, 劉萍, 李健, 等. 三疣梭子蟹基因組微衛(wèi)星特征分析[J]. 中國水產(chǎn)科學, 2008,15(5):738-743.

    [22] Zhang L,Bao Z, Cheng J, et al. Fosmid library construction and initial analysis of end sequences in Zhikong scallop (Chlamysfarreri) [J]. Marine Biotechnol, 2007,9(5):606-612.

    [23] Wang Y, Ren R, Yu Z, et al.Bioinformatic mining of EST-SSR loci in the Pacific oyster,Crassostreagigas[J]. Animal Genetics, 2008,39(3):287-289.

    [24] 崔建洲, 申雪艷, 楊官品, 等. 紅鰭東方鲀基因組微衛(wèi)星特征分析[J]. 中國海洋大學學報:自然科學版, 2006,36(2):249-254.

    [25] 吳文婷, 張磊, 宋志民, 等. 5個紫菜屬物種絲狀體的微衛(wèi)星初步遺傳分析[J]. 南方水產(chǎn)科學, 2012,8(4):29-36.

    BioinformaticsandMicrosatelliteSequencesAnalysisofESTSequenceinRidgeTailShrimpExopalaemoncarinicauda

    DUAN Yafei1, ZHANG Zhe1, LI Jitao2, 3, LI Jian2, 3, LIU Ping2, 3

    ( 1. Key Laboratory of South China Sea Fishery Resources Exploitation & Utilization, Ministry of Agriculture, South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou 510300, China; 2. Key Laboratory of Sustainable Development of Marine Fisheries, Ministry of Agriculture, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, China; 3. Function Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China )

    The analysis of bioinformatics and microsatellite sequences characteristics of expressed sequence tag (EST) sequence were studied, basing on the hemocyte cDNA library of ridge tail shrimpExopalaemoncarinicaudaconstructed by our laboratory. The 1053 unique sequences were yielded from 2853 high quality ESTs, including 329 contigs and 724 singletons. BLAST analysis revealed 593 (56.3%) of the unique sequences as orthologs of genes from other organisms, 22 (2.1%) of which were hypothetical protein, and the remaining 438 (41.6%) had no homology to any sequences. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that 181 of the unique sequences mapped to 120 of metabolic pathways. A total of 416 (14.58%) of the microsatellite repeat sequences were found using EST-SSR analysis. In the 416 repeat sequences, most was dinucleotide with repeats of 374 (89.90%), followed by trinucleotide with repeats of 35 (8.41%), and the third tetranucleotide with repeats of 7 (1.68%). The most common repeat types were AG, AAT and AAGT repeat. The findings provide important information for screening the functional gene resource and molecular markers in ridge tail shrimp.

    Exopalaemoncarinicauda; hemocyte; EST; bioinformatics ; microsatellite

    10.16378/j.cnki.1003-1111.2016.05.018

    S917

    A

    1003-1111(2016)05-0562-06

    2015-11-03;

    2015-12-30.

    國家蝦產(chǎn)業(yè)技術(shù)體系項目(CARS-47);國家自然科學基金資助項目(31472275);山東省泰山產(chǎn)業(yè)領(lǐng)軍人才工程項目(LNJY2015002);青島海洋科學與技術(shù)國家實驗室鰲山科技創(chuàng)新計劃項目(2015ASKJ02).

    段亞飛(1989—),男,助理研究員,碩士;研究方向: 蝦類免疫學. E-mail: duanyafei89@163.com. 通訊作者:劉萍(1962—),女,研究員;研究方向: 海水養(yǎng)殖生物種質(zhì)資源與遺傳育種. E-mail: liuping@ysfri.ac.cn.

    猜你喜歡
    脊尾白蝦微衛(wèi)星
    梅花鹿基因組微衛(wèi)星分布特征研究
    疫情當下,上半年華東地區(qū)小棚、工廠化、土塘白蝦如何應對?
    脊尾白蝦“紅”起來了畝產(chǎn)可達五百斤!
    林麝全基因組微衛(wèi)星分布規(guī)律研究
    四川動物(2017年4期)2017-07-31 23:54:19
    六價鉻離子在脊尾白蝦和三疣梭子蟹體內(nèi)的富集動力學
    基于轉(zhuǎn)錄組測序的波紋巴非蛤微衛(wèi)星標記研究
    苯并[a]芘在脊尾白蝦(Exopalaemon carinicauda)體內(nèi)富集的動力學研究
    中國“一箭雙星”成功將“遙感衛(wèi)星二十一號”與“天拓二號視頻微衛(wèi)星”發(fā)射升空
    河北遙感(2014年3期)2014-07-10 13:16:48
    脊尾白蝦仿生態(tài)規(guī)模養(yǎng)殖試驗
    午夜两性在线视频| 宅男免费午夜| 国产在线免费精品| 日韩一本色道免费dvd| 日本vs欧美在线观看视频| 日韩一卡2卡3卡4卡2021年| 婷婷丁香在线五月| 新久久久久国产一级毛片| 日日爽夜夜爽网站| 国产老妇伦熟女老妇高清| 亚洲中文字幕日韩| 精品人妻1区二区| av天堂在线播放| 另类亚洲欧美激情| 黑丝袜美女国产一区| 人成视频在线观看免费观看| 制服诱惑二区| 亚洲国产中文字幕在线视频| 亚洲精品av麻豆狂野| 老司机影院成人| 久久久久久人人人人人| 国产淫语在线视频| 九草在线视频观看| 免费观看av网站的网址| 中国美女看黄片| av福利片在线| 一本—道久久a久久精品蜜桃钙片| xxx大片免费视频| 日韩大码丰满熟妇| 国产成人精品在线电影| 十八禁人妻一区二区| 国产精品秋霞免费鲁丝片| 亚洲欧洲精品一区二区精品久久久| 国产av国产精品国产| 婷婷色综合www| 麻豆乱淫一区二区| 波多野结衣av一区二区av| 久久久精品免费免费高清| 久久精品熟女亚洲av麻豆精品| 99热全是精品| 亚洲国产最新在线播放| 国产精品欧美亚洲77777| 男女高潮啪啪啪动态图| 亚洲 欧美一区二区三区| 亚洲精品国产av成人精品| 午夜老司机福利片| 国产精品香港三级国产av潘金莲 | 欧美国产精品一级二级三级| 亚洲av美国av| 亚洲av片天天在线观看| 黄频高清免费视频| 国产精品免费视频内射| 久9热在线精品视频| 新久久久久国产一级毛片| 欧美少妇被猛烈插入视频| 欧美黑人欧美精品刺激| 精品卡一卡二卡四卡免费| 亚洲一码二码三码区别大吗| 黄网站色视频无遮挡免费观看| 在线观看免费午夜福利视频| 国产一区有黄有色的免费视频| 韩国高清视频一区二区三区| 久久国产精品影院| 亚洲国产av新网站| 精品国产一区二区三区四区第35| 在线观看国产h片| 女人久久www免费人成看片| 欧美性长视频在线观看| 一区二区三区精品91| www.自偷自拍.com| 一级毛片电影观看| 久久久国产精品麻豆| 丝袜脚勾引网站| 中文字幕制服av| 考比视频在线观看| 黑丝袜美女国产一区| 亚洲色图 男人天堂 中文字幕| 亚洲欧美色中文字幕在线| 久久国产精品人妻蜜桃| 成在线人永久免费视频| 又紧又爽又黄一区二区| 一区二区av电影网| 国产免费又黄又爽又色| 校园人妻丝袜中文字幕| 国产亚洲精品久久久久5区| 又紧又爽又黄一区二区| 日韩一卡2卡3卡4卡2021年| 久久毛片免费看一区二区三区| 日本av手机在线免费观看| 最近中文字幕2019免费版| 亚洲精品在线美女| 日本vs欧美在线观看视频| 精品人妻在线不人妻| 日日夜夜操网爽| 国产福利在线免费观看视频| 美女视频免费永久观看网站| xxx大片免费视频| 国产一区二区激情短视频 | 男女国产视频网站| 色网站视频免费| 亚洲精品在线美女| 欧美日韩亚洲国产一区二区在线观看 | 日韩熟女老妇一区二区性免费视频| 精品福利观看| e午夜精品久久久久久久| 亚洲国产精品999| 男人爽女人下面视频在线观看| 操出白浆在线播放| 国产不卡av网站在线观看| 咕卡用的链子| 一区在线观看完整版| 精品一品国产午夜福利视频| 热re99久久国产66热| 国产色视频综合| 日韩欧美一区视频在线观看| 99九九在线精品视频| 日本色播在线视频| 久久精品aⅴ一区二区三区四区| 丝袜美足系列| 国产精品成人在线| 午夜福利乱码中文字幕| 亚洲av片天天在线观看| 免费观看a级毛片全部| 赤兔流量卡办理| 中文乱码字字幕精品一区二区三区| 99热国产这里只有精品6| 国产成人精品在线电影| 18在线观看网站| 18禁裸乳无遮挡动漫免费视频| 国产成人av教育| 久9热在线精品视频| 国产精品成人在线| 超色免费av| 免费久久久久久久精品成人欧美视频| 日韩av在线免费看完整版不卡| 亚洲免费av在线视频| 亚洲男人天堂网一区| 秋霞在线观看毛片| 精品一区在线观看国产| 交换朋友夫妻互换小说| 国产xxxxx性猛交| 中文字幕精品免费在线观看视频| 亚洲欧洲国产日韩| 大香蕉久久网| 国产欧美日韩一区二区三区在线| 搡老岳熟女国产| 99re6热这里在线精品视频| 蜜桃在线观看..| 人妻 亚洲 视频| h视频一区二区三区| 男人添女人高潮全过程视频| 国产在线观看jvid| 日韩中文字幕欧美一区二区 | 老汉色∧v一级毛片| 久久热在线av| 男人操女人黄网站| 男女床上黄色一级片免费看| 超碰97精品在线观看| 校园人妻丝袜中文字幕| 一区在线观看完整版| 黄色一级大片看看| 亚洲 国产 在线| 午夜福利视频精品| 亚洲欧美成人综合另类久久久| 波野结衣二区三区在线| 丝瓜视频免费看黄片| 亚洲中文日韩欧美视频| 菩萨蛮人人尽说江南好唐韦庄| 老司机影院毛片| 亚洲图色成人| 中文字幕高清在线视频| e午夜精品久久久久久久| 你懂的网址亚洲精品在线观看| 亚洲欧美色中文字幕在线| www日本在线高清视频| 2018国产大陆天天弄谢| 天天躁夜夜躁狠狠久久av| 日韩电影二区| 水蜜桃什么品种好| 精品第一国产精品| 麻豆国产av国片精品| 亚洲av日韩精品久久久久久密 | 搡老岳熟女国产| 水蜜桃什么品种好| 亚洲国产精品一区三区| 亚洲色图 男人天堂 中文字幕| 亚洲精品美女久久av网站| 好男人电影高清在线观看| 王馨瑶露胸无遮挡在线观看| 精品熟女少妇八av免费久了| 天天添夜夜摸| 精品国产一区二区三区四区第35| 91九色精品人成在线观看| xxx大片免费视频| 免费看十八禁软件| 欧美黑人欧美精品刺激| 国产爽快片一区二区三区| 王馨瑶露胸无遮挡在线观看| 狂野欧美激情性xxxx| 电影成人av| 精品少妇黑人巨大在线播放| 成年av动漫网址| av电影中文网址| 夫妻午夜视频| 免费观看a级毛片全部| 夫妻性生交免费视频一级片| 亚洲国产毛片av蜜桃av| 一二三四在线观看免费中文在| 丰满迷人的少妇在线观看| 免费在线观看黄色视频的| 精品国产一区二区三区久久久樱花| 青春草视频在线免费观看| 日韩一区二区三区影片| 亚洲 欧美一区二区三区| 18禁裸乳无遮挡动漫免费视频| 成年美女黄网站色视频大全免费| h视频一区二区三区| 免费在线观看影片大全网站 | 中国国产av一级| 男人操女人黄网站| videosex国产| 亚洲欧美色中文字幕在线| 国产日韩欧美视频二区| 中文字幕高清在线视频| 午夜激情久久久久久久| 欧美精品一区二区免费开放| 国产成人精品在线电影| 日韩视频在线欧美| 国产日韩欧美视频二区| 在线观看免费高清a一片| 国产成人欧美在线观看 | 国产成人系列免费观看| 国产黄色视频一区二区在线观看| 久久综合国产亚洲精品| 两人在一起打扑克的视频| 亚洲av欧美aⅴ国产| 亚洲男人天堂网一区| 免费黄频网站在线观看国产| 亚洲成人免费电影在线观看 | 国产麻豆69| 又粗又硬又长又爽又黄的视频| 制服人妻中文乱码| 丝袜美足系列| 亚洲专区国产一区二区| 伊人久久大香线蕉亚洲五| 日韩中文字幕欧美一区二区 | 亚洲成av片中文字幕在线观看| 久久久国产欧美日韩av| 悠悠久久av| 老司机在亚洲福利影院| 好男人视频免费观看在线| 两个人看的免费小视频| 中文字幕色久视频| 男人操女人黄网站| 久久久亚洲精品成人影院| 欧美日韩综合久久久久久| 人人妻人人爽人人添夜夜欢视频| 中文字幕另类日韩欧美亚洲嫩草| 免费看十八禁软件| 国产麻豆69| 国产高清国产精品国产三级| 黄色毛片三级朝国网站| 欧美 亚洲 国产 日韩一| 国产精品久久久久久精品电影小说| 99久久99久久久精品蜜桃| 在线观看国产h片| 国产精品人妻久久久影院| 国产国语露脸激情在线看| 亚洲精品久久久久久婷婷小说| 国产熟女欧美一区二区| 国产在线视频一区二区| 国产精品香港三级国产av潘金莲 | 精品免费久久久久久久清纯 | 老司机靠b影院| 一级a爱视频在线免费观看| 欧美 亚洲 国产 日韩一| 精品少妇内射三级| 午夜免费鲁丝| 欧美国产精品一级二级三级| 亚洲成国产人片在线观看| 18禁黄网站禁片午夜丰满| 一区二区三区精品91| 国产成人精品无人区| 国产高清不卡午夜福利| 在线天堂中文资源库| 亚洲视频免费观看视频| 在线观看一区二区三区激情| 天堂俺去俺来也www色官网| 色婷婷久久久亚洲欧美| 一级毛片我不卡| 亚洲成人手机| 中文字幕av电影在线播放| 中国国产av一级| 一级毛片电影观看| 美女扒开内裤让男人捅视频| 各种免费的搞黄视频| 亚洲一卡2卡3卡4卡5卡精品中文| 视频在线观看一区二区三区| 男的添女的下面高潮视频| 天堂中文最新版在线下载| 美女脱内裤让男人舔精品视频| av有码第一页| av有码第一页| 19禁男女啪啪无遮挡网站| 欧美黄色片欧美黄色片| 国产不卡av网站在线观看| 999久久久国产精品视频| 建设人人有责人人尽责人人享有的| 国产精品免费视频内射| a 毛片基地| 下体分泌物呈黄色| 亚洲国产精品成人久久小说| 亚洲精品日本国产第一区| 亚洲欧美清纯卡通| 99九九在线精品视频| 精品卡一卡二卡四卡免费| 丝袜美足系列| 王馨瑶露胸无遮挡在线观看| 亚洲av片天天在线观看| 王馨瑶露胸无遮挡在线观看| 国产成人精品久久二区二区免费| 少妇猛男粗大的猛烈进出视频| 丁香六月欧美| 中文字幕高清在线视频| 热re99久久国产66热| 大片免费播放器 马上看| 国产成人免费无遮挡视频| 久久青草综合色| 校园人妻丝袜中文字幕| 国产av国产精品国产| 亚洲,一卡二卡三卡| 日日夜夜操网爽| 一级毛片女人18水好多 | 精品亚洲成国产av| 波多野结衣av一区二区av| 巨乳人妻的诱惑在线观看| 91成人精品电影| 欧美在线黄色| 亚洲人成网站在线观看播放| 18禁国产床啪视频网站| 尾随美女入室| 看免费成人av毛片| 人人妻人人澡人人爽人人夜夜| 99国产综合亚洲精品| 天堂中文最新版在线下载| 99国产精品一区二区蜜桃av | 另类亚洲欧美激情| 免费在线观看完整版高清| 老汉色av国产亚洲站长工具| 性高湖久久久久久久久免费观看| 99精国产麻豆久久婷婷| 免费观看av网站的网址| 欧美精品一区二区大全| 久久久亚洲精品成人影院| 亚洲人成电影观看| 如日韩欧美国产精品一区二区三区| 99精品久久久久人妻精品| 99久久精品国产亚洲精品| 中文字幕色久视频| 免费av中文字幕在线| 精品第一国产精品| 国产成人a∨麻豆精品| 2021少妇久久久久久久久久久| 亚洲国产精品成人久久小说| 国产精品一二三区在线看| 久久久久久久国产电影| 国产在线免费精品| 操美女的视频在线观看| 亚洲成人手机| 亚洲av在线观看美女高潮| 成人亚洲精品一区在线观看| 69精品国产乱码久久久| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品古装| 人人妻,人人澡人人爽秒播 | 午夜免费男女啪啪视频观看| 在线精品无人区一区二区三| 在线天堂中文资源库| 亚洲欧美精品自产自拍| 亚洲精品自拍成人| 人人妻人人澡人人看| 色综合欧美亚洲国产小说| 精品一区二区三区av网在线观看 | 精品国产国语对白av| 老司机亚洲免费影院| 午夜福利视频在线观看免费| 中文字幕另类日韩欧美亚洲嫩草| 十分钟在线观看高清视频www| 丰满少妇做爰视频| 精品人妻一区二区三区麻豆| 国产av精品麻豆| 18禁裸乳无遮挡动漫免费视频| 热99国产精品久久久久久7| 国产无遮挡羞羞视频在线观看| 午夜影院在线不卡| 成年人免费黄色播放视频| 国产三级黄色录像| 久久免费观看电影| 亚洲欧美成人综合另类久久久| 欧美另类一区| 女人爽到高潮嗷嗷叫在线视频| 黄网站色视频无遮挡免费观看| 亚洲av综合色区一区| 亚洲欧美色中文字幕在线| 欧美日本中文国产一区发布| 天天操日日干夜夜撸| 少妇人妻 视频| 在线观看国产h片| 欧美 亚洲 国产 日韩一| 久久精品熟女亚洲av麻豆精品| 日韩,欧美,国产一区二区三区| cao死你这个sao货| 国产熟女午夜一区二区三区| 亚洲综合色网址| 丁香六月欧美| 国产日韩欧美视频二区| 国产亚洲精品第一综合不卡| 国产又色又爽无遮挡免| 久久国产精品大桥未久av| 精品福利永久在线观看| 亚洲av成人精品一二三区| 日本猛色少妇xxxxx猛交久久| 色播在线永久视频| 国产视频一区二区在线看| 亚洲国产毛片av蜜桃av| 国产三级黄色录像| 男女国产视频网站| 亚洲第一av免费看| 人人妻人人澡人人看| 国产黄色视频一区二区在线观看| 久久中文字幕一级| 色播在线永久视频| 男女无遮挡免费网站观看| 午夜av观看不卡| 欧美精品啪啪一区二区三区 | 亚洲免费av在线视频| 久久精品久久精品一区二区三区| 亚洲欧美精品综合一区二区三区| 免费黄频网站在线观看国产| 五月开心婷婷网| 后天国语完整版免费观看| 欧美少妇被猛烈插入视频| 欧美性长视频在线观看| 一本色道久久久久久精品综合| 日韩大片免费观看网站| 国产av国产精品国产| 亚洲三区欧美一区| 久热这里只有精品99| 亚洲欧美精品自产自拍| 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 老司机影院毛片| 欧美成狂野欧美在线观看| 黑人猛操日本美女一级片| www日本在线高清视频| 午夜免费成人在线视频| 亚洲中文av在线| 制服诱惑二区| 国产成人av激情在线播放| 秋霞在线观看毛片| 久久久精品免费免费高清| 国产成人精品久久二区二区91| 爱豆传媒免费全集在线观看| 中文字幕高清在线视频| 久久久久久久国产电影| 国产一区二区三区综合在线观看| 一级,二级,三级黄色视频| 国产男人的电影天堂91| 免费在线观看视频国产中文字幕亚洲 | 女人久久www免费人成看片| 波多野结衣一区麻豆| 国产黄色免费在线视频| 日本av手机在线免费观看| 国产真人三级小视频在线观看| 国产精品偷伦视频观看了| 91精品伊人久久大香线蕉| 久久国产精品男人的天堂亚洲| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| av在线老鸭窝| 满18在线观看网站| 大片电影免费在线观看免费| www日本在线高清视频| 无遮挡黄片免费观看| 黑丝袜美女国产一区| 一区二区av电影网| 久久精品国产综合久久久| 丝瓜视频免费看黄片| 成年动漫av网址| 大陆偷拍与自拍| 成人手机av| 熟女少妇亚洲综合色aaa.| 黄色视频在线播放观看不卡| 赤兔流量卡办理| 91国产中文字幕| 校园人妻丝袜中文字幕| 在线亚洲精品国产二区图片欧美| av有码第一页| 人人澡人人妻人| 精品一品国产午夜福利视频| 人人妻人人澡人人看| a 毛片基地| 99精国产麻豆久久婷婷| 久久综合国产亚洲精品| 天天影视国产精品| 免费在线观看黄色视频的| 少妇 在线观看| 黄频高清免费视频| xxxhd国产人妻xxx| 精品亚洲成a人片在线观看| 国产成人精品久久二区二区91| 欧美少妇被猛烈插入视频| 黄网站色视频无遮挡免费观看| 精品少妇黑人巨大在线播放| 久久久精品国产亚洲av高清涩受| 欧美日韩国产mv在线观看视频| 搡老岳熟女国产| 这个男人来自地球电影免费观看| av国产精品久久久久影院| 亚洲国产精品国产精品| 亚洲七黄色美女视频| 在线观看国产h片| 好男人电影高清在线观看| 大话2 男鬼变身卡| 中国国产av一级| 嫩草影视91久久| 久久精品国产亚洲av高清一级| 男女国产视频网站| 久久久久久久大尺度免费视频| 美女脱内裤让男人舔精品视频| 午夜久久久在线观看| 国产一级毛片在线| 国产熟女欧美一区二区| 欧美另类一区| 男人舔女人的私密视频| 国产成人精品无人区| 校园人妻丝袜中文字幕| 亚洲成人免费av在线播放| 丝袜人妻中文字幕| 亚洲中文日韩欧美视频| 日韩精品免费视频一区二区三区| 大香蕉久久成人网| 热re99久久国产66热| 在线观看免费视频网站a站| 如日韩欧美国产精品一区二区三区| 丰满少妇做爰视频| 午夜免费男女啪啪视频观看| 又大又黄又爽视频免费| 老司机影院成人| 国产在线视频一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 另类精品久久| 欧美xxⅹ黑人| 黄色怎么调成土黄色| 日本91视频免费播放| 欧美日韩成人在线一区二区| 91九色精品人成在线观看| kizo精华| 欧美黄色淫秽网站| 热re99久久精品国产66热6| 少妇裸体淫交视频免费看高清 | 黄色 视频免费看| 亚洲av日韩在线播放| 欧美日韩视频高清一区二区三区二| 新久久久久国产一级毛片| 免费高清在线观看视频在线观看| 99久久99久久久精品蜜桃| 国产精品国产av在线观看| 少妇人妻久久综合中文| 丰满饥渴人妻一区二区三| 人体艺术视频欧美日本| 国产精品一区二区免费欧美 | 少妇 在线观看| 欧美激情高清一区二区三区| 精品一区二区三区av网在线观看 | 久久精品熟女亚洲av麻豆精品| 99热国产这里只有精品6| 一级黄色大片毛片| 中文字幕另类日韩欧美亚洲嫩草| 校园人妻丝袜中文字幕| 两个人免费观看高清视频| 亚洲国产av影院在线观看| 黄网站色视频无遮挡免费观看| 男女之事视频高清在线观看 | 五月天丁香电影| 五月开心婷婷网| 欧美久久黑人一区二区| 青春草视频在线免费观看| 少妇人妻久久综合中文| 老汉色∧v一级毛片| 视频在线观看一区二区三区| 超碰成人久久| 波野结衣二区三区在线| 我的亚洲天堂| 80岁老熟妇乱子伦牲交| 可以免费在线观看a视频的电影网站| 欧美日韩亚洲高清精品| 亚洲人成电影观看| 亚洲午夜精品一区,二区,三区| 久久国产精品大桥未久av| 亚洲 欧美一区二区三区| 亚洲,欧美精品.| 夫妻午夜视频| 国产xxxxx性猛交| 成人国产av品久久久| 精品视频人人做人人爽| 两人在一起打扑克的视频| 欧美日韩亚洲高清精品| 大香蕉久久成人网| 丰满迷人的少妇在线观看| 精品一品国产午夜福利视频| 中文字幕色久视频| 99久久人妻综合| 日韩伦理黄色片| 国产一级毛片在线|