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

    大青葉蟬內(nèi)共生菌的檢測及其16S rRNA基因序列分析

    2017-06-05 14:59:42練啟仙劉健鋒楊茂發(fā)周小蕓
    環(huán)境昆蟲學報 2017年2期
    關鍵詞:檢測

    練啟仙,劉健鋒,楊茂發(fā)*,羅 亭,吳 廣,周小蕓

    (1. 貴州大學昆蟲研究所,貴州省山地農(nóng)業(yè)病蟲害重點實驗室,貴陽 550025;2. 興義民族師范學院,貴州興義 562400)

    大青葉蟬內(nèi)共生菌的檢測及其16S rRNA基因序列分析

    練啟仙1,2,劉健鋒1,楊茂發(fā)1*,羅 亭1,吳 廣1,周小蕓1

    (1. 貴州大學昆蟲研究所,貴州省山地農(nóng)業(yè)病蟲害重點實驗室,貴陽 550025;2. 興義民族師范學院,貴州興義 562400)

    用細菌通用引物27F-1513R對大青葉蟬內(nèi)共生菌的16S rRNA基因進行PCR擴增,瓊脂糖凝膠電泳,克隆、測序、序列比對與分析,用貝葉斯法和極大似然法構建系統(tǒng)發(fā)育樹。結果表明:從大青葉蟬體內(nèi)鑒定出Sulciamuelleri和Sodalissp.2種內(nèi)共生菌,Sulciamuelleri隸屬于擬桿菌門Bacteroidetes黃桿菌綱Flavobacteria黃桿菌目Flavobacteriales黃桿菌科Flavobacteriaceae,該內(nèi)共生菌獲得13個新單倍型,14個變異位點,其序列間的一致率為和遺傳距離分別為99.40%-99.80%和0.001-0.004;系統(tǒng)發(fā)育樹顯示,Sulciamuelleri與大葉蟬亞科中檢測到的Sulciamuelleri聚在一起,其變異程度低,所有單倍型為同一個種。Sodalissp.屬于變形桿菌門Proteobacteria伽馬變形菌綱Gammaproteobacteria腸桿菌目Enterobacteriale腸桿菌科Enterobacteriaceae,本研究獲得5個新單倍型,20個變異位點,其序列間的一致率為96.10%-99.70%,遺傳距離范圍為0.012-0.040;系統(tǒng)發(fā)育樹顯示,Sodalissp.與采采蠅檢測到的Sodalisglossinidius聚在一起,其變異程度相對較大,屬于同一個屬。

    大青葉蟬;內(nèi)共生菌;16S rRNA基因;單倍型;系統(tǒng)發(fā)育

    內(nèi)共生菌(Endosymbiont)和昆蟲之間的共生關系是一個普遍存在的現(xiàn)象。Buchner(1965)用組織學的方法檢測到昆蟲具有兩種內(nèi)共生菌,即初生內(nèi)共生菌(Primary endosymbiont)和次生內(nèi)共生菌(Secondary endosymbiont)。初生內(nèi)共生菌生活在昆蟲的特定細胞中,與昆蟲幾乎同時發(fā)生,協(xié)同進化(Clarketal., 1992),為寄主提供必要的營養(yǎng)成分,如氨基酸等(Baumann, 2005);次生內(nèi)共生菌可使雄蟲雌性化(Werren and Windsor, 2000; Negrietal., 2008, 2009)、有殺雄作用(Werrenetal., 1994; Lawson,etal., 2001; von der Schulenburgetal., 2001)以及影響昆蟲的生長發(fā)育等(Kontsedalov, 2008; Nodaetal., 2012)。刺吸式口器的昆蟲主要是刺吸植物韌皮部和木質部的汁液為食,然而植物的汁液主要成分是蔗糖形式的碳水化合物,缺乏脂類和蛋白質,碳水化合物可合成脂類,但缺乏了必須氨基酸就不能合成蛋白質,而植物汁液中存在的氨基酸多數(shù)是不必要的氨基酸(Moranetal., 2008),為此,大多數(shù)刺吸植物汁液的昆蟲需要體內(nèi)的共生菌為其提供必需氨基酸等營養(yǎng)成分(Montlloretal., 2002; Daleetal., 2003; Douglas, 2009)。

    大青葉蟬Cicadellaviridis隸屬于半翅目Hemiptera葉蟬科Cicadellidae大葉蟬亞科Cicadellinae,可為害160多種植物,在為害過程中還傳帶一些病毒病,是重要的農(nóng)林害蟲(朱弘復和鄧國藩,1950;康芝仙等,1996)。該害蟲主要取食寄主植物嫩綠的枝梢莖葉汁液,和其他昆蟲一樣體內(nèi)具有共生菌,具有為其提供營養(yǎng)等作用。Takiya等(2006)在大青葉蟬中檢測到Sulciamuelleri內(nèi)共生菌;Michalik等(2014)檢測有Sulciamuelleri初生內(nèi)共生菌和Sodalis-likeendosymbiont次生內(nèi)共生菌,Lian等(2016)檢測有Rickettsiasp.次生內(nèi)共生菌。在中國,目前對大青葉蟬的生物學特性和化學防治研究較多,而對其內(nèi)共生菌的研究較少,本實驗利用通用引物對大青葉蟬內(nèi)共生菌的16S rRNA 基因進行PCR擴增、克隆及測序等,然后對檢測出的內(nèi)共生菌Sulciamuelleri和Sodalissp.序列的特征進行了詳細分析,以期為內(nèi)共生菌對寄主的生長發(fā)育的影響、不同地理種群大青葉蟬內(nèi)共生菌的差異以及大青葉蟬地理種群的遺傳分化與共生菌的協(xié)同進化等進一步深入研究奠定基礎。

    1 材料與方法

    1.1 材料

    1.1.1 供試蟲源

    于2013-2014年在貴州省貴陽市花溪區(qū)采集大青葉蟬成蟲,75%酒精清洗3次后用無水乙醇浸泡,于-20℃下保存,供檢測其內(nèi)共生菌用。

    1.2 方法

    1.2.1 大青葉蟬DNA的提取

    將-20℃下保存的大青葉蟬雌雄成蟲放在無菌濾紙上,待乙醇揮發(fā)后用無菌水清洗5次,然后置于1.5 mL的離心管中,用研磨棒將其磨碎,最后按照血液/組織、細胞基因組提取試劑盒說明書(天根生化科技(北京)有限公司)提取DNA,提取的DNA置于-20℃下保存?zhèn)溆谩?/p>

    1.2.2 PCR擴增

    用細菌通用引物27F(5′-AGAGTTTGATCM TGGCTCAG-3′)和1513R(5′-ACGGYTACCTTGTT ACGACTT-3′)(Weisburgetal., 1991)對提取的DNA進行PCR擴增。反應體系為25 μL,其中包括2.5 μL 10×buffer(含Mg2+),2 μL dNTP(Takara,寶生物工程(大連)有限公司),引物各0.5 μL(10 μM),0.5 μLTaqDNA聚合酶(Takara),1 μL DNA模板,最后加18 μL ddH2O。PCR反應程序為95℃預變性5 min,95℃變性30 s,63.5℃退火30 s,72℃延伸1 min,30個循環(huán),最后72℃延伸10 min。擴增完成后取4 μL用1.5%瓊脂糖凝膠電泳檢測。

    1.2.3 克隆與序列的測定

    PCR產(chǎn)物電泳回收純化后,連接到pMD?18-T載體上(Takara),連接產(chǎn)物轉化到Trans1-T2感受態(tài)細胞中,涂板并37℃培養(yǎng)過夜,然后挑取白色菌落至600 μL的LB培養(yǎng)液(含5 μL,50 μg/mL Amp)中,37℃下震蕩10 h左右。隨機挑選PCR鑒定的陽性克隆送上海生工公司測序,每個PCR產(chǎn)物挑取10-12個單菌落用于序列測定。

    1.2.4 序列分析

    測序后獲得的序列用DNAstar 5.03軟件包中的SeqMan軟件去除載體,結合測序峰圖,對序列進行人工校對,并拼接正反序列。然后把序列放在GenBank中進行比對,下載與之對應的大青葉蟬內(nèi)共生菌Sulciamuelleri和Sodalissp.的16S rRNA基因序列(AY676915,KJ660084,KJ660085)(Takiyaetal., 2006; Michaliketal., 2014),用Mega 6.05在默認設置下進行比對,并用DnaSP 5軟件包計算所獲序列的單倍型個數(shù)以及單倍型變異位點、多樣度(Hd)、多樣度分歧、核苷酸多樣度(Pi)和核苷酸差異均數(shù)(K),用Bioedit 7.0對大青葉蟬單倍型序列進行一致性檢驗,用Mega 6.05計算其遺傳距離。

    從GenBank中另外選取一致性較高的序列與之相對應的測序序列用Mega 6.05進行比對,用Jmodetest軟件計算出最佳模型,最后用MrBayes 3.2 和PAUP 4.0建系統(tǒng)發(fā)育樹,極大似然法(Maximum-likelihood)分支置信度采用Bootstrap重復1000次。

    2 結果與分析

    2.1 大青葉蟬內(nèi)共生菌16S rRNA基因序列在GenBank中Blast同源序列檢索

    對10頭大青葉蟬(雌雄各5頭)提取的DNA,用27F-1513R引物對16S rRNA基因進行PCR擴增后,用1.5%瓊脂糖凝膠電泳檢測,其片段大約為1500 bp(圖1),然后對PCR產(chǎn)物進行純化、克隆與測序,并在GenBank中Blast同源序列檢索,結果表明,克隆的所有序列中有50條克隆序列(31條來自于雌蟲,19條來自于雄蟲)與擬桿菌門Bacteroidetes黃桿菌綱Flavobacteria黃桿菌目Flavobacteriales黃桿菌科Flavobacteriaceae的Sulciamuelleri同源性很高,DNA序列片段大小為1522-1523 bp;28條克隆序列(18條來自于雌蟲,10條來自于雄蟲)與變形桿菌門Proteobacteria伽馬變形菌綱Gammaproteobacteria腸桿菌目Enterobacteriales腸桿菌科Enterobacteriaceae的Sodalis屬同源性很高,DNA序列片段大小為1499-1511 bp。同時,這兩種內(nèi)共生菌在10頭大青葉蟬中均被檢測到。

    圖1 大青葉蟬內(nèi)共生菌16S rRNA基因擴增產(chǎn)物Fig.1 PCR prouducts of endosymbiotic bacteria 16S rRNA gene from Cicadella viridis注:M,2000 bp分子質量標準;3/5/7/9/10,雌蟲;1/2/4/6/8,雄蟲;11,空白對照。Note: M, 2000 bp DNA marker; 3/5/7/9/10, Female;1/2/4/6/8, Male; 11, Blank control.

    2.2 單倍型序列特征分析

    2.2.1Sulciamuelleri16S rRNA基因序列特征分析

    獲得的50條Sulciamuelleri克隆序列與GenBank中原有大青葉蟬內(nèi)共生菌Sulciamuelleri的兩條序列(AY676915,KJ660084)進行分析,結果表明,從10個樣品中獲得14個單倍型,新單倍型有13個(Hap2-14)。來自10個樣品的35條克隆序列(雌22條,雄13條)與AY676915和KJ660084共享單倍型1(Hap1),該單倍型的比率為70%;Hap10有來自兩頭雌蟲的兩條序列共享,其余12個單倍型均為來自同一雌或雄個體的克隆序列(表1)。這些單倍型序列中檢測到14個變異位點(Variable sites),該單倍型多樣度(Haplotype diversity,Hd)為0.496,單倍型多樣度分歧(Variance of Haplotype diversity)為0.0075,核苷酸多樣度(Nucleotide diversity,Pi)為0.00107±0.00025,平均核苷酸差異(Average number of pairwise nucleotide differences,k)為0.64932。序列中堿基T、C、A、G的平均含量分別為25.3%、18.0%、30.3%、26.5%,A+T含量(55.5%)高于G+C含量(44.5%);其序列的一致率在99.40%-99.80%,一致性非常高(表2)。此外,表1顯示,所有單倍型無插入缺失,但有堿基變換,既有顛換,也有轉換。表3顯示,單倍型的遺傳距離范圍為0.001-0.004。

    2.2.2Sodalissp. 16S rRNA基因序列特征分析

    用Dnasp 5.0軟件對本實驗獲得的28條Sodalissp.序列與GenBank中的大青葉蟬Sodalis-like endosymbiont序列(KJ660085)進行單倍型分析。結果表明,共獲得6個單倍型;新單倍型5個(Hap2-6),有2個單倍型(Hap2,4)均有來自雌雄個體的克隆序列共享,KJ660085單獨為一個單倍型(表4)。這些單倍型序列中檢測到20個變異位點,該單倍型多樣度Hd為0.766,單倍型多樣度分歧為0.00221,核苷酸多樣度Pi為0.001016±0.00192,平均核苷酸差異k為4.42857。新單倍型序列中堿基T、C、A、G的平均含量分別為20.6%、22.5%、25.4%、31.5%,此基因序列中富含G+C(54.0%)。從表4中可看出,Sodalissp.單倍型堿基發(fā)生了變換,但沒有插入缺失。表5顯示,該單倍型序列的一致率為96.10%-99.70%。遺傳距離最小為0.002,最大為0.040(表6)。

    表 1 大青葉葉蟬內(nèi)共生菌Sulcia muelleri不同單倍型的變異位點

    表2 大青葉蟬內(nèi)共生菌Sulcia muelleri單倍型序列一致性檢驗

    表3 大青葉蟬內(nèi)共生菌Sulcia muelleri單倍型序列的遺傳距離

    表4 大青葉葉蟬內(nèi)共生菌Sodalis sp.不同單倍型變異位點

    表 5 大青葉蟬內(nèi)共生菌Sodalis sp.單倍型序列一致性檢驗

    表 6 大青葉蟬內(nèi)共生菌Sodalis sp.單倍型序列的遺傳距離

    2.3 單倍型序列GenBank登錄號

    大青葉蟬內(nèi)共生菌Sulciamuelleri和Sodalis sp. 16S rRNA基因序列提交到GenBank數(shù)據(jù)庫中,獲得了相應的GenBank登錄號,其中Sulciamuelleri16S rRNA基因序列的登錄號為KY765568-KY765581,Sodalis sp. 16S rRNA基因序列登錄號為KY765563-KY765567。

    2.4 單倍型序列系統(tǒng)發(fā)育樹的構建

    通過 DnaSP 5.0軟件計算獲得的單倍型序列與GenBank中的序列用Mega 6.05進行比對,然后用jModeltest計算,獲得貝葉斯(Bayesian,BA)和極大似然法(Maximum-likelihood,ML)建樹的最佳模型,用Mrbayes 3.2和PAUP 4.0構建Sulciamuelleri和Sodalis sp. 16S rRNA基因序列的BA樹和ML樹,結果表明:用兩種方法構建的系統(tǒng)發(fā)育樹是一致的。大青葉蟬Sulciamuelleri內(nèi)共生菌的14個單倍型序列與來自大葉蟬亞科的Sulciamuelleri(Bacteroidetes)序列形成同一個進化支,其支持率為1/98(BA/ML);同時可從圖中看出,半翅目昆蟲的Sulciamuelleri內(nèi)共生菌16S rRNA基因序列以相同亞科的昆蟲形成同一個進化分支(圖2)。大青葉蟬內(nèi)共生菌Sodalissp.與采采蠅中檢測到的Sodalisglossinidius聚在一起,其支持率為1/87(BA/ML),但6個單倍型序列在不同的進化支上;Hap1(KJ660084)、Hap4單獨為一支,Hap2和Hap3處在一個進化支上,Hap5和Hap6形成一個進化支(圖3)。

    圖2 大青葉蟬內(nèi)共生菌Sulcia muelleri 16S rRNA 基因系列系統(tǒng)發(fā)育分析Fig.2 Molecular phylogenetic analysis of the Sulcia muelleri endosymbiont from Cicadella viridis based on the 16S ribosomal RNA gene sequence

    圖3 大青葉蟬內(nèi)共生菌Sodalis sp. 16S rRNA 基因系列系統(tǒng)發(fā)育分析.Fig.3 Phylogenetic placement of Sodalis sp.endosymbiont of Cicadella viridis on the basis of the 16S ribosomal RNA gene sequence

    3 結論與討論

    Buchner(1965)研究認為,每種葉蟬中至少有兩種類型的內(nèi)共生菌。本研究用細菌通用引物27F-1513R檢測到大青葉蟬中有Sulciamuelleri、Sodalissp.和Rickettsiasp.(Lianetal. , 2016)3種內(nèi)共生菌,和前人的結論一致。本研究與國外的大青葉蟬中檢測到的內(nèi)共生菌Sulciamuelleri,其16S r RNA基因序列之間的一致率非常高,達99.0%以上,而Sodalissp. 16S rRNA基因序列的一致率在96.1%-99.7%;內(nèi)共生菌Sulciamuelleri16S r RNA基因序列之間遺傳距離范圍為0.001-0.004,Sodalissp. 16S rRNA基因序列的遺傳距離最大為0.04,最小為0.002。由此說明,雖然大青葉蟬的地理分布差異很大,但Sulciamuelleri16S rRNA基因的變異不大,相對保守,這與McCutcheon等(2010)的研究結果(Sulciamuelleri的基因有非常高的保守性)是一致的;而Sodalissp. 16S rRNA基因的變異程度相對較大。Devereux等(1990)認為,16SrDNA序列一致率低于98%,屬于不同的種,一致率低于95%,則屬于不同的屬; 王佳玉(1983)根據(jù)遺傳距離,提出劃分種的范圍:0.05-0.30,由此推測,Sulciasp.所有單倍型屬于同一個種,即Sulciamuelleri;本研究獲得的Sodalissp.屬于同一個屬,基因序列之間發(fā)生了一定的變異,導致其變異的原因和機制有待研究。

    本研究基于16S rRNA基因的進行了系統(tǒng)發(fā)育分析,發(fā)現(xiàn)大青葉蟬體內(nèi)的Sulciamuelleri共生菌與葉蟬科昆蟲體內(nèi)的Sulciamuelleri共生菌最為接近,而與菱蠟蟬科昆蟲體內(nèi)的Sulciamuelleri共生菌距離較遠,同時相同亞科昆蟲的Sulciamuelleri聚在一起,與昆蟲的分類系統(tǒng)一致,這一結果與Noda等(2012)分析構建的黑尾葉蟬Nephotettixcincticeps內(nèi)共生菌Sulciamuelleri的系統(tǒng)發(fā)育樹的結果相同。Sodalissp.為次生內(nèi)共生菌,最早在采采蠅Glossinamorsitansmorsitans中發(fā)現(xiàn)(Dale and Maudlin, 1999),后來在象鼻蟲(Tojuetal., 2013)等鞘翅目昆蟲中也檢測到了該內(nèi)共生菌,目前在已研究的葉蟬科昆蟲中,僅在大青葉蟬中檢測到。

    據(jù)研究,Sulciamuelleri是葉蟬、沫蟬、蠟蟬等半翅目刺吸式昆蟲中普遍存在的初生內(nèi)共生菌(Moranetal., 2005; Urban and Cryan, 2012; Kogaetal., 2013),是一種古老的內(nèi)共生菌;Sodalissp.在采采蠅中報道為次生內(nèi)共生菌,可分離純培養(yǎng)(Dale and Maudlin, 1999)。Thao等(2000)認為初生內(nèi)共生菌16S rRNA基因序列的(A+T)mol%含量高,而次生內(nèi)共生菌16S rRNA基因序列富含(G+C)mol%,Moran等(1998)認為低的(G+C)mol%是一個反應寄生物與寄主在遠古時代就結合的特征,而高的(G+C)mol%則反應其接近自由生活細菌的方面。大青葉蟬內(nèi)共生菌Sulciamuelleri富含(A+T)mol%(55.5%),Sodalissp. 16S rRNA 基因序列(G+C)mol%含量高(54.0%),由此可見,本實驗結果與前人的研究結果一致,Sulciamuelleri為大青葉蟬的初生內(nèi)共生菌,低的(G+C)mol%含量,說明與宿主起源的一致性;Sodalissp.為次生內(nèi)共生菌,但是Sodalissp.能否分離純培養(yǎng)仍需進一步研究。

    蟬亞目Cicadomorpha的昆蟲具有雙重初生內(nèi)共生菌,即Sulciamuelleri具有伴隨初生內(nèi)共生菌(Urban and Cryan, 2012),如:Sulciamuelleri+Baumanniacicadellinicola(Takiyaetal., 2006; Nodaetal., 2012)、Sulciamuelleri+Zinderiainsecticola(McCutcheon and Moran, 2010)、Sulciamuelleri+Vidaniafulgoroideae(Urban and Cryan, 2012)、Sulciamuelleri+Nasuiadeltocephalinicola(Ishiietal., 2013)、Sulciamuelleri+Hodgkiniacicadicola(McCutcheonetal. , 2009)。大青葉蟬Sulciamuelleri內(nèi)共生菌是否具有伴隨初生內(nèi)共生菌,本次對測序獲得的所有序列進行了Blast檢索與分析,未檢測到,有待進一步研究。

    本研究對Sulciamuelleri和Sodalissp. 16S rRNA基因序列做了詳細分析,其結果為內(nèi)共生菌對寄主的生長發(fā)育的影響以及大青葉蟬地理種群的遺傳分化與共生菌的協(xié)同進化等進一步深入研究提供了基礎資料。

    References)

    Baumann P. Biology of bacteriocyte-associated endosymbionts of plant sap-sucking insects [J].AnnualReviewofMicrobiology, 2005, 59: 155-189.

    Buchner P. Endosymbiosis of Animals with Plant Microorganisms [M]. New York: Wiley, 1965.

    Clark MA, Baumann L, Munson MA,etal. The eubacterial endosymbionts of whiteflies (Homoptera: Aleyrodoidea) constitute a lineage distinct from the endosymbionts of aphids and mealybugs [J].CurrentMicrobiology, 1992, 25 (2): 119-123.

    Dale C, Maudlin I.Sodalisgen. nov. andSodalisglossinidiussp. nov., a microaerophilic secondary endosymbiont of the tsetse flyGlossinamorsitansmorsitans[J].InternationalJournalofSystematicandEvolutionaryMicrobiology, 1999, 49 (1): 267-275.

    Dale C, Wang B, Moran N,etal. Loss of DNA recombinational repair enzymes in the initial stages of genome degeneration [J].MolecularBiologyandEvolution, 2003, 20 (8): 1188-1194.

    Devereux R, He SH, Doyle CEA,etal. Diversity and origin ofDesulfovibriospecies: Phylogenetic definition of a family [J].JournalofBacteriology, 1990, 172 (7): 3609-3619.

    Douglas AE. The microbial dimension in insect nutritional ecology [J].FunctionalEcology, 2009, 23 (1): 38-47.

    Ishii Y, Matsuura Y, Kakizawa S,etal. Diversity of bacterial endosymbionts associated withMacrostelesleafhoppers vectoring phytopathogenic phytoplasmas [J].AppliedandEnvironmentalMicrobiology, 2013, 79 (16): 5013-5022.

    Kang ZX, Lu H, Yi BR,etal. Study on the biological charactors ofTettigellaviridis[J].JournalofJilinAgriculturalUniversity, 1996, 18 (3): 19-26. [康芝仙, 路紅, 伊伯仁, 等. 大青葉蟬生物學特性的研究 [J]. 吉林農(nóng)業(yè)大學學報, 1996, 18 (3): 19-26]

    Koga R, Bennett GM, Cryan JR,etal. Evolutionary replacement of obligate symbionts in an ancient and diverse insect lineage [J].EnvironmentalMicrobiology, 2013, 15 (7): 2073-2081.

    Kontsedalov S, Zchori-Fein E, Chiel E,etal. The presence ofRickettsiais associated with increased susceptibility ofBemisiatabaci(Homoptera: Aleyrodidae) to insecticides [J].PestManagementScience, 2008, 64 (8): 789-792.

    Lawson ET, Mousseau TA, Klaper R,etal.Rickettsiaassociated with male-killing in a buprestid beetle [J].Heredity, 2001, 86 (4): 497-505.

    Lian QX, Liu JF, Yang MF,etal. Molecular identification and detection ofRickettsiaendosymbiont in the green leafhopper:Cicadellaviridis(Hemiptera: Cicadellidae) [J].JournalofAsia-PacificEntomology, 2016, 19: 659-664.

    McCutcheon JP, McDonald BR, Moran NA. Convergent evolution of metabolic roles in bacterial co-symbionts of insects [J].ProceedingsoftheNationalAcademyofSciences, 2009, 106 (36): 15394-15399.

    McCutcheon JP, Moran NA. Functional convergence in reduced genomes of bacterial symbionts spanning 200 My of evolution [J].GenomeBiologyandEvolution, 2010, 2: 708-718.

    Michalik A, Jankowska W, Kot M,etal. Symbiosis in the green leafhopper,Cicadellaviridis(Hemiptera, Cicadellidae). Association instatunascendi?[J]ArthropodStructure&Development, 2014, 43 (6): 579-587.

    Montllor CB, Maxmen A, Purcell AH. Facultative bacterial endosymbionts benefit pea aphidsAcyrthosiphonpisumunder heat stress [J].EcologicalEntomology, 2002, 27 (2): 189-195.

    Moran NA, McCutcheon JP, Nakabachi A. Genomics and evolution of heritable bacterial symbionts [J].AnnualReviewofGenetics, 2008, 42: 165-190.

    Moran NA, Telang A. Bacteriocyte-associated symbionts of insects [J].Bioscience, 1998, 48 (4): 295-304.

    Moran NA, Tran P, Gerardo NM. Symbiosis and insect diversification: An ancient symbiont of sap-feeding insects from the bacterial phylum Bacteroidetes [J].AppliedandEnvironmentalMicrobiology, 2005, 71 (12): 8802-8810.

    Negri I, Franchini A, Gonella E,etal. Unravelling theWolbachiaevolutionary role: The reprogramming of the host genomic imprinting [J].ProceedingsoftheRoyalSocietyofLondonB:BiologicalSciences, 2009, 276: 2485-2491.

    Negri I, Franchini A, Mandrioli M,etal. The gonads ofZyginidiapullulamales feminized byWolbachiapipientis[J].BulletinofInsectology, 2008, 61 (1): 213-214.

    Noda H, Watanabe K, Kawai S,etal. Bacteriome-associated endosymbionts of the green rice leafhopperNephotettixcincticeps(Hemiptera: Cicadellidae) [J].AppliedEntomologyandZoology, 2012, 47 (3): 217-225.

    Takiya DM, Tran PL, Dietrich CH,etal. Co-cladogenesis spanning three phyla: leafhoppers (Insecta: Hemiptera: Cicadellidae) and their dual bacterial symbionts [J].MolecularEcology, 2006, 15 (13): 4175-4191.

    Thao ML, Moran NA, Abbot P,etal. Cospceiation of psyllids and their primary Prokaryotic endosymbionts [J].AppliedandEnvironmentalMicrobiology, 2000, 66 (7): 2898-2905.

    Toju H, Tanabe AS, Notsu Y,etal. Diversification of endosymbiosis: Replacements, co-speciation and promiscuity of bacteriocyte symbionts in weevils [J].TheISMEJournal, 2013, 7 (7): 1378-1390.

    Urban JM, Cryan JR. Two ancient bacterial endosymbionts have coevolved with the planthoppers (Insecta: Hemiptera: Fulgoroidea) [J].BMCEvolutionaryBiology, 2012, 12 (1): 87-106.

    von der Schulenburg JHG, Habig M, Sloggett JJ,etal. Incidence of male killingRickettsiaspp. (α-proteobacteria) in the ten-spot ladybird beetleAdaliadecempunctataL. (Coleoptera: Coccinellidae) [J].AppliedandEnvironmentalMicrobiology, 2001, 67 (1): 270-277.

    Wang JY. Molecular Population Genetics and Evolution [M]. Beijing: China Agriculture Press, 1983:169-196. [王家玉(譯).分子群體遺傳學與進化論[M].北京: 中國農(nóng)業(yè)出版社, 1983: 169-196]

    Weisburg WG, Barns SM, Pelletier DA,etal. 16S ribosomal DNA amplification for phylogenetic study [J].JournalofBacteriology, 1991, 173 (2): 697-703.

    Werren JH, Hurst GDD, Zheng W,etal.Rickettsialrelative associated with male killing in the ladybird beetle (Adaliadipunctata) [J].JournalofBacteriology, 1994, 176 (2): 388-394.

    Werren, JH, Windsor D.Wolbachiainfection frequencies in insects: Evidence of a global equilibrium?[J].ProceedingsoftheRoyalSocietyofLondonB:BiologicalSciencesB, 2000, 267: 1277-1285.

    Zhu HF, Deng GF. Lifespan of the green leafhopper:Cicadellaviridis(Hemiptera: Cicadellidae) [J].ActaEntomologicaSinica, 1950, 1 (1): 14-40. [朱弘復, 鄧國藩. 青葉跳蟬的生活史[J]. 中國昆蟲學報, 1950,1 (1): 14-40]

    Detection of endosymbionts and the sequences analysis of 16S rRNA gene in the green leafhopper (Cicadellaviridis)

    LIAN Qi-Xian1,2, LIU Jian-Feng1, YANG Mao-Fa1*, LUO Ting1, WU Guang1, ZHOU Xiao-Yun1

    (1. Guizhou Provincial Key Laboratory for Agricultural Pest Management of the Mountainous Regions, Institute of Entomology,Guizhou University, Guiyang 550025, China; 2. Department of Chemical Biology, Xingyi Normal University for Nationalities, Xingyi 562400, Guizhou Province, China)

    PCR was performed using universal primers (27F-1513R) for 16S rRNA gene of endosymbiotic bacteria fromCicadellaviridis, followed by agarose gel electrophoresis, cloning, sequencing, sequence alignment and analysis. The phylogenetic trees were generated by the Bayesian and maximum-likelihood methods. The results showed thatSulciasp. was classified as Flavobacteriaceae, Flavobacteriales, Flavobacteria Bacteroidetes. The 16S rRNA gene sequences ofSulciasp. got 13 new haplotypes and 14 variable sites, and identity rate and genetic distance of which were 99.40%-99.80% and 0.001-0.004, respectively. The phylogenetic tree revealed thatSulciasp. was associated withSulciamuelleridetected from the leafhoppers of Cicadellinae.Sulcissp. had a low degree of variation, and all haplotypes belonged toSulciamuellerispecies.Sodalissp. belonged to Enterobacteriaceae, Enterobacteriales, Gammaproteobacteria, Proteobacteria. The 16S rRNA gene sequences ofSodalissp. obtained 5 new haplotypes and 20 variable sites in this study, and identity rate and genetic distance of which were 96.10%-99.70% and 0.002-0.040, respectively. The phylogenetic tree showed thatSodalissp. was most closely related toSodalisglossinidiusdetected from the Tsetse flies.Sodalissp. has a high degree of variation, and all haplotypes belong one genus.

    Cicadellaviridis; endosymbiont; 16S rRNA gene; haplotype; phylogeny

    貴州省高層次創(chuàng)新型人才(“百”層次)項目(黔科合人才20164022);貴州省研究生卓越人才計劃項目(黔教研合ZYRC字[2013]010號);中國煙草總公司貴州省公司科技項目(中煙黔科[2016]6號);貴州省農(nóng)科院專項(黔農(nóng)科院自主創(chuàng)新科研專項字[2014]009號);正安縣綠色產(chǎn)業(yè)發(fā)展辦公室委托項目。

    練啟仙,女,1978年生,貴州湄潭人,在讀博士,副教授,研究方向為昆蟲生態(tài)與害蟲防治,E-mail: llx0713@163.com

    *通訊作者Author for correspondence, E-mail: gdgdly@126.com

    Received: 2016-12-08;接受日期Accepted: 2017-03-24

    Q963;S433.39

    A

    1674-0858(2017)02-0342-09

    練啟仙,劉健鋒,楊茂發(fā),等.大青葉蟬內(nèi)共生菌的檢測及其16S rRNA基因序列分析[J].環(huán)境昆蟲學報,2017,39(2):342-350.

    猜你喜歡
    檢測
    QC 檢測
    “不等式”檢測題
    “一元一次不等式”檢測題
    “一元一次不等式組”檢測題
    “幾何圖形”檢測題
    “角”檢測題
    “有理數(shù)的乘除法”檢測題
    “有理數(shù)”檢測題
    “角”檢測題
    “幾何圖形”檢測題
    啦啦啦免费观看视频1| 免费观看精品视频网站| 免费搜索国产男女视频| 人人妻人人爽人人添夜夜欢视频| 亚洲一区二区三区色噜噜| 成人国产综合亚洲| 看黄色毛片网站| www.熟女人妻精品国产| 中出人妻视频一区二区| 欧美日韩一级在线毛片| 一二三四在线观看免费中文在| 亚洲成国产人片在线观看| 成熟少妇高潮喷水视频| 免费不卡黄色视频| 久久久久国内视频| 亚洲成av人片免费观看| 亚洲精华国产精华精| 亚洲国产欧美一区二区综合| 18禁美女被吸乳视频| 亚洲无线在线观看| 黑丝袜美女国产一区| 91麻豆精品激情在线观看国产| 一级a爱片免费观看的视频| 乱人伦中国视频| 精品人妻1区二区| 女人高潮潮喷娇喘18禁视频| 国产精品香港三级国产av潘金莲| 久久狼人影院| 欧美 亚洲 国产 日韩一| 久久精品国产综合久久久| 精品一区二区三区av网在线观看| 日本 欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 高清黄色对白视频在线免费看| 国内毛片毛片毛片毛片毛片| 淫秽高清视频在线观看| 露出奶头的视频| 欧美激情高清一区二区三区| 黄频高清免费视频| 淫秽高清视频在线观看| 午夜福利免费观看在线| 日韩视频一区二区在线观看| 99精品欧美一区二区三区四区| 久久中文看片网| 国产一卡二卡三卡精品| 精品人妻在线不人妻| 日韩欧美一区视频在线观看| 一区福利在线观看| av欧美777| 亚洲美女黄片视频| 午夜福利18| 久久人人精品亚洲av| 婷婷精品国产亚洲av在线| 波多野结衣高清无吗| 久久精品91蜜桃| 亚洲精品国产区一区二| 亚洲中文字幕一区二区三区有码在线看 | 国产精品国产高清国产av| 少妇熟女aⅴ在线视频| 天天一区二区日本电影三级 | 久久香蕉国产精品| 久久精品人人爽人人爽视色| 成人亚洲精品一区在线观看| 成人免费观看视频高清| 国产在线精品亚洲第一网站| www.熟女人妻精品国产| 久久久久久久久久久久大奶| 亚洲精品久久国产高清桃花| 亚洲自拍偷在线| 久久人人97超碰香蕉20202| 久久婷婷成人综合色麻豆| 亚洲全国av大片| av天堂在线播放| 日本欧美视频一区| 亚洲精品久久成人aⅴ小说| 日韩欧美三级三区| 天堂动漫精品| 亚洲中文av在线| 亚洲av第一区精品v没综合| 亚洲欧美一区二区三区黑人| 久久精品91无色码中文字幕| 亚洲色图 男人天堂 中文字幕| 日本一区二区免费在线视频| 久久久精品国产亚洲av高清涩受| 精品熟女少妇八av免费久了| 可以在线观看的亚洲视频| 精品无人区乱码1区二区| 性少妇av在线| 精品午夜福利视频在线观看一区| 精品久久蜜臀av无| 亚洲第一欧美日韩一区二区三区| 欧美中文综合在线视频| 免费在线观看视频国产中文字幕亚洲| 精品高清国产在线一区| 欧美国产精品va在线观看不卡| 欧美不卡视频在线免费观看 | 国产亚洲精品久久久久5区| cao死你这个sao货| 成年人黄色毛片网站| 99久久精品国产亚洲精品| 国产激情久久老熟女| 国产日韩一区二区三区精品不卡| 亚洲精品av麻豆狂野| 午夜福利影视在线免费观看| 国产精品99久久99久久久不卡| 国产亚洲av高清不卡| 亚洲专区中文字幕在线| 涩涩av久久男人的天堂| 国产不卡一卡二| 国产成年人精品一区二区| or卡值多少钱| 国产成人精品久久二区二区91| 成人手机av| 久久这里只有精品19| 熟女少妇亚洲综合色aaa.| 国产人伦9x9x在线观看| 91大片在线观看| 丁香六月欧美| 国产精品乱码一区二三区的特点 | 欧美日韩亚洲综合一区二区三区_| 久久人妻熟女aⅴ| a在线观看视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 99久久国产精品久久久| 可以免费在线观看a视频的电影网站| 91大片在线观看| 国产亚洲精品av在线| 亚洲午夜精品一区,二区,三区| 亚洲精品在线美女| 纯流量卡能插随身wifi吗| 亚洲av成人不卡在线观看播放网| 亚洲欧美日韩无卡精品| 色综合亚洲欧美另类图片| 国产成年人精品一区二区| 色综合站精品国产| 成在线人永久免费视频| 搡老妇女老女人老熟妇| 久久精品国产亚洲av高清一级| 夜夜爽天天搞| 一区在线观看完整版| 免费不卡黄色视频| 法律面前人人平等表现在哪些方面| av天堂久久9| 91精品国产国语对白视频| 激情在线观看视频在线高清| 12—13女人毛片做爰片一| 成人三级黄色视频| 精品乱码久久久久久99久播| 久久久久久国产a免费观看| cao死你这个sao货| 在线观看日韩欧美| 久久久国产成人精品二区| 国内久久婷婷六月综合欲色啪| 日日摸夜夜添夜夜添小说| 欧美乱色亚洲激情| 国产亚洲av高清不卡| 嫩草影视91久久| 午夜a级毛片| 色综合站精品国产| 免费搜索国产男女视频| 一级作爱视频免费观看| 大码成人一级视频| 夜夜躁狠狠躁天天躁| 亚洲国产日韩欧美精品在线观看 | 亚洲一码二码三码区别大吗| 国产成年人精品一区二区| 在线免费观看的www视频| 黄色视频,在线免费观看| 国产成人影院久久av| 欧美一级毛片孕妇| 精品国产乱子伦一区二区三区| 国产精品亚洲美女久久久| 亚洲精品一卡2卡三卡4卡5卡| 日韩精品青青久久久久久| 日韩精品免费视频一区二区三区| 长腿黑丝高跟| 欧美日韩亚洲综合一区二区三区_| 桃红色精品国产亚洲av| 在线av久久热| 精品久久久久久久人妻蜜臀av | 99国产极品粉嫩在线观看| 国产私拍福利视频在线观看| 午夜亚洲福利在线播放| 午夜精品久久久久久毛片777| 成熟少妇高潮喷水视频| 真人做人爱边吃奶动态| 亚洲人成网站在线播放欧美日韩| 美女免费视频网站| 91麻豆精品激情在线观看国产| 制服丝袜大香蕉在线| 国产熟女午夜一区二区三区| 久久久久久久久免费视频了| 成人国产综合亚洲| 神马国产精品三级电影在线观看 | 一区福利在线观看| 午夜福利一区二区在线看| 午夜亚洲福利在线播放| av片东京热男人的天堂| 欧美最黄视频在线播放免费| 侵犯人妻中文字幕一二三四区| 亚洲成人久久性| 久久中文看片网| 51午夜福利影视在线观看| 精品国产一区二区久久| 天堂影院成人在线观看| 手机成人av网站| 欧美日本亚洲视频在线播放| 亚洲中文av在线| 亚洲五月婷婷丁香| 狂野欧美激情性xxxx| 久久国产亚洲av麻豆专区| 曰老女人黄片| 乱人伦中国视频| 国产私拍福利视频在线观看| 亚洲欧美一区二区三区黑人| 日本撒尿小便嘘嘘汇集6| 国产在线精品亚洲第一网站| 精品久久久久久久久久免费视频| 久久人人爽av亚洲精品天堂| 麻豆一二三区av精品| 久久久久久人人人人人| 色哟哟哟哟哟哟| 精品人妻在线不人妻| 午夜a级毛片| 91精品三级在线观看| 国产精品 国内视频| 操出白浆在线播放| www.精华液| av电影中文网址| 亚洲精品在线美女| 国产精品久久久久久精品电影 | 色在线成人网| 在线免费观看的www视频| 大型av网站在线播放| x7x7x7水蜜桃| 国产aⅴ精品一区二区三区波| av电影中文网址| 少妇粗大呻吟视频| 午夜福利欧美成人| 国产精品精品国产色婷婷| 国产欧美日韩一区二区精品| 国产精品 欧美亚洲| 欧美 亚洲 国产 日韩一| 日本五十路高清| 在线天堂中文资源库| 久久久久九九精品影院| 精品欧美国产一区二区三| 日本三级黄在线观看| 久久精品国产99精品国产亚洲性色 | av福利片在线| 亚洲avbb在线观看| 老司机深夜福利视频在线观看| 久久精品91蜜桃| 欧美午夜高清在线| 夜夜爽天天搞| 中出人妻视频一区二区| 91在线观看av| 纯流量卡能插随身wifi吗| 在线永久观看黄色视频| 亚洲性夜色夜夜综合| 精品熟女少妇八av免费久了| 久久国产精品人妻蜜桃| 久久精品亚洲熟妇少妇任你| 精品电影一区二区在线| 久久亚洲精品不卡| 国产高清有码在线观看视频 | 久久国产精品人妻蜜桃| 精品不卡国产一区二区三区| 久久亚洲精品不卡| 欧美成人性av电影在线观看| av福利片在线| 亚洲 欧美一区二区三区| 亚洲一区二区三区色噜噜| 美女高潮喷水抽搐中文字幕| 一级片免费观看大全| 精品午夜福利视频在线观看一区| 久久精品国产综合久久久| 人人妻,人人澡人人爽秒播| 亚洲av电影在线进入| 99久久国产精品久久久| 窝窝影院91人妻| 在线观看www视频免费| 国产精品,欧美在线| 精品久久久久久久毛片微露脸| 欧美色欧美亚洲另类二区 | 啦啦啦 在线观看视频| 欧美人与性动交α欧美精品济南到| √禁漫天堂资源中文www| 香蕉久久夜色| 亚洲国产毛片av蜜桃av| 午夜亚洲福利在线播放| 很黄的视频免费| 真人做人爱边吃奶动态| 少妇粗大呻吟视频| 美女高潮到喷水免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 91老司机精品| 亚洲精品国产区一区二| 久久久水蜜桃国产精品网| 热99re8久久精品国产| 波多野结衣巨乳人妻| 激情视频va一区二区三区| 久久精品国产亚洲av香蕉五月| www.www免费av| 国产人伦9x9x在线观看| 国产午夜福利久久久久久| 欧美黄色片欧美黄色片| 制服诱惑二区| 亚洲 欧美一区二区三区| 丝袜在线中文字幕| 国产一区二区激情短视频| 怎么达到女性高潮| 久久国产精品影院| 免费av毛片视频| 青草久久国产| 中文字幕av电影在线播放| 久久久久久人人人人人| 黑丝袜美女国产一区| 午夜精品久久久久久毛片777| 免费在线观看日本一区| 久久久久久久久中文| 午夜视频精品福利| 制服丝袜大香蕉在线| 精品久久久久久成人av| 少妇 在线观看| 99热只有精品国产| 黄片大片在线免费观看| 国产aⅴ精品一区二区三区波| 亚洲国产高清在线一区二区三 | 精品少妇一区二区三区视频日本电影| 欧美午夜高清在线| 黄片大片在线免费观看| 怎么达到女性高潮| 国产一区二区三区视频了| 亚洲精品美女久久久久99蜜臀| 久99久视频精品免费| 亚洲午夜理论影院| 欧美av亚洲av综合av国产av| 亚洲黑人精品在线| 人妻丰满熟妇av一区二区三区| 国产欧美日韩一区二区三区在线| 国内毛片毛片毛片毛片毛片| 久久草成人影院| 9热在线视频观看99| 日韩欧美一区视频在线观看| 又大又爽又粗| 国内精品久久久久精免费| 97碰自拍视频| 色播亚洲综合网| 人人澡人人妻人| 免费看a级黄色片| 18美女黄网站色大片免费观看| 亚洲欧美日韩高清在线视频| 国产精品久久久久久亚洲av鲁大| 免费看美女性在线毛片视频| 欧美中文综合在线视频| 美女国产高潮福利片在线看| 国产91精品成人一区二区三区| 亚洲成av人片免费观看| 真人做人爱边吃奶动态| 涩涩av久久男人的天堂| 国产高清有码在线观看视频 | 中文亚洲av片在线观看爽| 老汉色av国产亚洲站长工具| 波多野结衣一区麻豆| 欧美日本视频| 久久狼人影院| 亚洲成av人片免费观看| 看黄色毛片网站| 国产不卡一卡二| 中国美女看黄片| 国产精品野战在线观看| 欧美日本亚洲视频在线播放| 久久亚洲真实| 啦啦啦韩国在线观看视频| 亚洲少妇的诱惑av| 高清黄色对白视频在线免费看| 制服人妻中文乱码| 老司机靠b影院| 丝袜美腿诱惑在线| 欧美国产日韩亚洲一区| 老熟妇仑乱视频hdxx| а√天堂www在线а√下载| 正在播放国产对白刺激| 黄色 视频免费看| 琪琪午夜伦伦电影理论片6080| 亚洲成国产人片在线观看| 国产亚洲av高清不卡| 91成年电影在线观看| 亚洲精品国产色婷婷电影| 日本免费一区二区三区高清不卡 | 黑人操中国人逼视频| 日韩欧美一区二区三区在线观看| 国产极品粉嫩免费观看在线| 欧美绝顶高潮抽搐喷水| 18禁国产床啪视频网站| 国产成人av教育| 欧美激情久久久久久爽电影 | 亚洲第一电影网av| 18禁观看日本| 51午夜福利影视在线观看| 免费在线观看日本一区| 97人妻精品一区二区三区麻豆 | 欧美黑人欧美精品刺激| 久久久精品国产亚洲av高清涩受| 在线十欧美十亚洲十日本专区| 亚洲久久久国产精品| 成年女人毛片免费观看观看9| 欧美成人午夜精品| 90打野战视频偷拍视频| 国产精品 欧美亚洲| 免费看a级黄色片| 国产成人影院久久av| 亚洲国产日韩欧美精品在线观看 | 国产伦一二天堂av在线观看| 久久天堂一区二区三区四区| 国产国语露脸激情在线看| 曰老女人黄片| 久久久久久大精品| √禁漫天堂资源中文www| 乱人伦中国视频| 成年人黄色毛片网站| 在线观看舔阴道视频| 欧美黑人精品巨大| 一级作爱视频免费观看| 岛国在线观看网站| 色综合站精品国产| 亚洲欧洲精品一区二区精品久久久| 精品人妻1区二区| 在线观看免费日韩欧美大片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲第一青青草原| 老司机午夜十八禁免费视频| 一级黄色大片毛片| 90打野战视频偷拍视频| 成人18禁在线播放| 亚洲精品在线观看二区| 免费人成视频x8x8入口观看| 免费观看精品视频网站| 99re在线观看精品视频| av福利片在线| 国产精品99久久99久久久不卡| 久久久久久久午夜电影| 国产精品 欧美亚洲| 涩涩av久久男人的天堂| 两性夫妻黄色片| 亚洲专区字幕在线| АⅤ资源中文在线天堂| 一本大道久久a久久精品| 韩国精品一区二区三区| 最近最新免费中文字幕在线| 成人永久免费在线观看视频| 国产精品,欧美在线| 国产激情久久老熟女| 午夜视频精品福利| 人成视频在线观看免费观看| 在线观看一区二区三区| 国产av在哪里看| 午夜精品在线福利| 麻豆成人av在线观看| 欧美日韩亚洲综合一区二区三区_| 高潮久久久久久久久久久不卡| 久久亚洲精品不卡| 非洲黑人性xxxx精品又粗又长| 久久久久久免费高清国产稀缺| 18禁美女被吸乳视频| 免费观看精品视频网站| 国产精品免费一区二区三区在线| 亚洲一区高清亚洲精品| 亚洲国产欧美日韩在线播放| 国产成人影院久久av| 两性午夜刺激爽爽歪歪视频在线观看 | 日韩精品免费视频一区二区三区| 免费在线观看黄色视频的| 一边摸一边做爽爽视频免费| 精品免费久久久久久久清纯| 国产成人精品久久二区二区免费| 国产高清激情床上av| 国产97色在线日韩免费| 自线自在国产av| 高清毛片免费观看视频网站| 不卡av一区二区三区| 久久精品人人爽人人爽视色| 亚洲情色 制服丝袜| 婷婷六月久久综合丁香| 日本一区二区免费在线视频| 亚洲自偷自拍图片 自拍| 一a级毛片在线观看| 免费女性裸体啪啪无遮挡网站| 欧美成人一区二区免费高清观看 | 国语自产精品视频在线第100页| 久久久久亚洲av毛片大全| 大型av网站在线播放| 久久 成人 亚洲| 午夜福利,免费看| 久久影院123| 搡老岳熟女国产| 在线观看日韩欧美| 久久精品91蜜桃| 这个男人来自地球电影免费观看| 久久精品亚洲精品国产色婷小说| 美女高潮喷水抽搐中文字幕| 国产精品 欧美亚洲| 亚洲欧美日韩无卡精品| 国产成人一区二区三区免费视频网站| 日韩精品中文字幕看吧| 亚洲狠狠婷婷综合久久图片| 欧美一级a爱片免费观看看 | 级片在线观看| av在线播放免费不卡| 欧美丝袜亚洲另类 | 自拍欧美九色日韩亚洲蝌蚪91| 亚洲全国av大片| 国产亚洲av嫩草精品影院| 免费无遮挡裸体视频| 日韩欧美免费精品| 亚洲中文日韩欧美视频| 99国产精品免费福利视频| 亚洲欧洲精品一区二区精品久久久| 一级a爱片免费观看的视频| 精品国产一区二区久久| 欧美亚洲日本最大视频资源| 日韩大码丰满熟妇| 岛国在线观看网站| 99国产综合亚洲精品| 日日爽夜夜爽网站| 99国产精品免费福利视频| 最新美女视频免费是黄的| 婷婷六月久久综合丁香| 国产亚洲欧美精品永久| 91成人精品电影| 久久久国产成人免费| 亚洲欧美日韩另类电影网站| 亚洲国产精品成人综合色| 欧美乱码精品一区二区三区| 国产一区二区激情短视频| 国产精品野战在线观看| 变态另类丝袜制服| 免费观看精品视频网站| 国产精品,欧美在线| 黄色丝袜av网址大全| 久久亚洲精品不卡| 亚洲欧美精品综合久久99| 国产三级在线视频| 精品一区二区三区四区五区乱码| 黄色 视频免费看| 高清毛片免费观看视频网站| 黑人巨大精品欧美一区二区mp4| 午夜福利在线观看吧| 亚洲av五月六月丁香网| 大型av网站在线播放| 美女国产高潮福利片在线看| 久久精品国产亚洲av高清一级| 91九色精品人成在线观看| 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 午夜福利免费观看在线| 国产成人精品久久二区二区免费| 女性生殖器流出的白浆| 成人精品一区二区免费| 久久久精品欧美日韩精品| 精品久久久精品久久久| 免费高清在线观看日韩| 咕卡用的链子| 国产黄a三级三级三级人| 亚洲欧美日韩高清在线视频| 老汉色av国产亚洲站长工具| 国产99久久九九免费精品| 19禁男女啪啪无遮挡网站| 午夜久久久久精精品| 国产一卡二卡三卡精品| 久久午夜综合久久蜜桃| 国产精品免费一区二区三区在线| 日本在线视频免费播放| 国产私拍福利视频在线观看| 欧美久久黑人一区二区| 亚洲精品美女久久av网站| 久久人人爽av亚洲精品天堂| 国产欧美日韩精品亚洲av| 免费一级毛片在线播放高清视频 | 人妻久久中文字幕网| 国产高清激情床上av| 性欧美人与动物交配| 老司机福利观看| 在线观看免费视频网站a站| 国产aⅴ精品一区二区三区波| 国产黄a三级三级三级人| 香蕉久久夜色| 亚洲aⅴ乱码一区二区在线播放 | 欧美激情极品国产一区二区三区| 在线国产一区二区在线| 91av网站免费观看| 亚洲最大成人中文| 一级,二级,三级黄色视频| 中文亚洲av片在线观看爽| 亚洲色图 男人天堂 中文字幕| 午夜精品国产一区二区电影| 国产97色在线日韩免费| 国产真人三级小视频在线观看| 久久草成人影院| 午夜福利一区二区在线看| 9色porny在线观看| 久久精品影院6| 久久久久久久久中文| 国内毛片毛片毛片毛片毛片| 久久久久国产精品人妻aⅴ院| 一级,二级,三级黄色视频| 黄网站色视频无遮挡免费观看| 久9热在线精品视频| 国产精品综合久久久久久久免费 | 色播在线永久视频| 久久人妻熟女aⅴ| 国产精品98久久久久久宅男小说| www.www免费av| 不卡一级毛片|