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

    高榕榕果內(nèi)Eupristina屬兩種榕小蜂的遺傳進化關(guān)系

    2013-12-09 06:35:36陳友鈴孫伶俐武蕾蕾吳文珊嚴菊媛
    生態(tài)學報 2013年19期
    關(guān)鍵詞:雌蜂雄蜂小蜂

    陳友鈴,孫伶俐,武蕾蕾,吳文珊,嚴菊媛,劉 燕

    (福建師范大學生命科學學院,福建省發(fā)育與神經(jīng)生物學重點實驗室,福州 350117)

    高榕榕果內(nèi)Eupristina屬兩種榕小蜂的遺傳進化關(guān)系

    陳友鈴,孫伶俐,武蕾蕾,吳文珊*,嚴菊媛,劉 燕

    (福建師范大學生命科學學院,福建省發(fā)育與神經(jīng)生物學重點實驗室,福州 350117)

    對生長在福州地區(qū)的高榕進行長期追蹤觀察,發(fā)現(xiàn)高榕榕果內(nèi)僅生活著Eupristinaaltissima和Eupristinasp.榕小蜂,前者為高榕的傳粉小蜂,后者無傳粉行為,兩者雌蜂之間在體色、觸角、花粉袋和花粉刷等部位存在細微的差異,而兩者雄蜂之間無形態(tài)差異。通過克隆福建地區(qū)5個樣地的高榕榕果內(nèi)收集到的E.altissima和Eupristinasp.榕小蜂,以及細葉榕的傳粉小蜂Eupristinaverticillata(外群)的Cytb及COI基因,并進行堿基組成及遺傳距離分析,用鄰接法構(gòu)建系統(tǒng)發(fā)育樹,分析兩榕小蜂群體之間的遺傳進化關(guān)系,結(jié)果顯示:(1)榕小蜂COI及Cytb序列堿基組成中A+T的含量(Cytb序列中A+T=75.3%,COI序列中A+T=75.5%)顯著高于G+C,符合膜翅目昆蟲線粒體基因堿基組成特征。(2)對兩群體小蜂進行遺傳距離分析顯示,Cytb序列中E.altissima和Eupristinasp. 群體內(nèi)各樣本之間的平均遺傳距離分別為0.0092和0.0030,而E.altissima與Eupristinasp. 群體間的平均距離為0.1588;COI序列中E.altissima與Eupristinasp. 群體內(nèi)各樣本之間的平均遺傳距離分別為0.0065和0.0205,而二者群體間的平均遺傳距離為0.1043,表明兩者群體間的遺傳距離明顯大于各自群體內(nèi)各樣本間的遺傳距離。統(tǒng)計GenBank中下載的6個屬34種榕小蜂Cytb序列的種間遺傳距離為0.0811—0.1723,6個屬28種榕小蜂COI序列的種間遺傳距離為0.0939—0.1986。由此認為E.altissima和Eupristinasp.之間的遺傳距離差異已經(jīng)達到了種間水平,即E.altissima與Eupristinasp.為兩個不同的種。(3)在形態(tài)上,兩種小蜂的雌蜂之間有微小差異,而二者雄蜂之間無差異,但Cytb與COI序列分析結(jié)果一致表明:E.altissima與Eupristinasp.雄蜂之間,以及二者雌蜂之間的遺傳距離均差異顯著,表明形態(tài)變異滯后于基因變異。雌蜂在表型上進化快于雄蜂,可能是由于雌蜂羽化后從榕果出飛,受到外界環(huán)境因素的影響較大,且兩種雌蜂在傳粉功能上存在差異,故二者之間的形態(tài)差異較大,而雄蜂壽命短,又終生生活在黑暗封閉、環(huán)境變化相對恒定的榕果內(nèi),兩種雄蜂在行為上不存在差異,故二者表型變異較為緩慢。E.altissima和Eupristinasp.小蜂對宿主的專一性不強,在榕-蜂協(xié)同進化過程中,可能發(fā)生過宿主轉(zhuǎn)移事件。

    高榕;高榕傳粉小蜂;DNA條形碼;遺傳結(jié)構(gòu);Cytb;COI

    高榕(FicusaltissimaBlume)隸屬于??崎艑伲菩弁?,常綠喬木,高可達30m。榕果成對腋生,橢圓狀或卵圓形,幼時包藏于早落的風帽狀苞片內(nèi),無柄,成熟時紅色或帶黃色。高榕原產(chǎn)于印度、尼泊爾、緬甸、泰國、馬來西亞、印度尼西亞、菲律賓等地區(qū),在我國主要分布于海南、廣西、四川、云南南部至中部等地[1- 2]。1981年Nair等人首次報道了印度地區(qū)的高榕榕果中的傳粉小蜂為Eupristinaaltissima[3],但未發(fā)現(xiàn)Eupristinasp. 的存在。引種到美國的弗羅里達州栽培的高榕,其傳粉小蜂亦為E.altissima,同樣在該地的高榕榕果內(nèi)未發(fā)現(xiàn)Eupristinasp.[4- 6]。高榕在中國引種栽培的歷史悠久,在長江以南地區(qū),高榕作為行道樹而廣為栽培。2003年谷海燕等發(fā)現(xiàn)云南西雙版納地區(qū)的高榕榕果內(nèi)有23種非傳粉小蜂,以及兩種傳粉小蜂:E.altissima和Eupristinasp.,兩者差異主要體現(xiàn)在體色和觸角索節(jié)數(shù)量的不同[7];2008年彭艷瓊等對這兩種小蜂的傳粉效率進行了比較,發(fā)現(xiàn)Eupristinasp.幾乎沒有傳粉能力,在形態(tài)上Eupristinasp. 無傳粉小蜂所特有的花粉刷,且花粉筐極度退化[8]。目前未發(fā)現(xiàn)兩種小蜂的雄蜂在形態(tài)上存在差異。本實驗室自2005年至2012年間對生長在福建地區(qū)的高榕進行長期的追蹤觀察,結(jié)果表明,福建高榕榕果中只有兩種小蜂:即Eupristinaaltissima和Eupristinasp.,而無其他種類的非傳粉小蜂,本實驗擬利用分子標記技術(shù),揭示高榕榕果內(nèi)這兩個形態(tài)高度相似,功能性行為迥異的小蜂種群之間的遺傳進化關(guān)系。

    2003年,加拿大的Hebert教授以鱗翅目昆蟲為研究對象,提出了用基因片段來區(qū)分物種的DNA條形碼技術(shù)(DNA barcoding)[9]。DNA條形碼技術(shù)是用短的DNA片段對物種進行識別和鑒定的分子生物學技術(shù)。線粒體細胞色素b(Cytb)和細胞色素c氧化酶亞基I(COI)是廣泛應(yīng)用于動物物種鑒定的DNA條形碼[10- 13]。線粒體Cytb與COI基因具有進化速度適中,序列沒有內(nèi)含子,很少存在插入和缺失的特點,且在細胞內(nèi)為多拷貝能很容易地利用引物擴增,是進行物種鑒定的優(yōu)秀標記基因[14- 15]。近年,Cytb與COI基因被廣泛運用于榕小蜂與榕樹的系統(tǒng)演化[16- 17]、榕小蜂新種和隱種的發(fā)現(xiàn)[18- 20]、榕小蜂宿主特異性分析[21- 22]及榕小蜂雄性多型鑒定[23]等的相關(guān)研究中。

    本文通過克隆高榕榕果內(nèi)的Eupristinaaltissima和Eupristinasp.榕小蜂,以及細葉榕(Ficusmicrocarpa)榕果的傳粉小蜂Eupristinaverticillata(外群)的COI和Cytb基因,并與通過GenBank下載的榕小蜂科(Agaonidae)榕小蜂的COI和Cytb序列進行系統(tǒng)發(fā)育樹的構(gòu)建,探討共存于高榕榕果內(nèi)的兩種榕小蜂之間的遺傳關(guān)系,并進一步推測它們可能的遺傳進化模式。

    1 材料和方法

    1.1 實驗材料

    實驗中所涉及的實驗材料及采集樣地見表1。

    表1 實驗材料、樣本代碼及采集樣地

    “-”表示無此樣本

    1.2 實驗方法

    1.2.1 實驗樣本的采集與保存

    將野外采集處于雄花期即將出蜂的高榕和細葉榕榕果帶回到實驗室,利用單果收蜂的方法[24]收集榕小蜂。在解剖鏡下,對出飛的雌蜂形態(tài)進行觀察和鑒定,收集單個榕果內(nèi)出飛的雌蜂為單一種類(E.altissima與Eupristinasp.兩種小蜂僅存在一種)的榕果內(nèi)的雌、雄蜂,分別存放在含95%乙醇溶液的1.5 mL的離心管中,使小蜂迅速死亡,貼上標簽,置于4 ℃保存。為了保證保存液始終為高濃度的乙醇,在保存期間要定期更換乙醇。

    1.2.2 基因組DNA的提取及目的基因的擴增

    取95%酒精下儲存的榕小蜂樣品用酚氯仿抽提法進行基因組DNA提取。對提取的基因組DNA進行目的片段的PCR擴增。擴增引物如下:

    COI-F 5′CAACATTTATTTTGATTTTTTGG3′

    COI-R 5′TCCAATGCACTAATCTGCCATATTA3′

    Cytb-F 5′TATGTACTACCATGAG GACAAATATC3′

    Cytb-R 5′ATTACACCTCCTAATTTATTAGGAAT3′

    COI基因擴增條件為:94 ℃預變性3 min;94 ℃變性30 s,54℃退火30 s,72℃延伸40 s(34個循環(huán));72 ℃延伸7 min。Cytb基因擴增條件為:94 ℃預變性4 min;94 ℃變性30 s,57 ℃退火60 s,72 ℃延伸60 s (35個循環(huán));72 ℃延伸5 min。

    對PCR擴增的產(chǎn)物直接進行測序,所有PCR產(chǎn)物均送至生工生物公司進行測序。將測序所得序列上傳GenBank上。

    1.2.3 序列比對和系統(tǒng)發(fā)育樹的建立

    測序結(jié)果利用CLUSTAL X 1.83進行序列比對,BioEdit 3.3校正比對結(jié)果,用MEGA4.0軟件進行實驗樣本序列間的堿基替換機堿基組成分析,并用Mega4.0軟件中的鄰接(Neighbor-Joining NJ)法對COI和Cytb序列建立系統(tǒng)發(fā)育樹并進行分析。

    2 結(jié)果

    2.1 Cytb與COI基因序列分析

    2.1.1 Cytb基因片段堿基組成分析

    采用Mega4.0軟件,對E.altissima和Eupristinasp.小蜂(20個樣本)校準后的371 bp序列進行堿基組成分析,結(jié)果表明:Cytb基因序列中,不變位點317個,變異位點52個,簡約信息位點47個,自裔位點5個,變異性為14.0%。

    T,C,A,G四種堿基的平均含量分別為41.5%,14.9%,33.8%,9.8%,A+T的含量為75.3%,而G+C的含量為24.7%,A+T的相對含量明顯高于G+C的含量,這是膜翅目昆蟲線粒體基因序列的共同特點。

    Cytb密碼子不同位點的核苷酸頻率及堿基替換情況見表2,由表2可知,在密碼子的各個位點上A+T含量均偏高,尤其在第二位點上A+T的含量高達85.4%,偏向最為明顯,其中T的使用最為頻繁達到46.7%,反映出Cyt b 基因在密碼子使用上具有偏向性。E.altissima和Eupristinasp.小蜂的Cytb基因序列上核苷酸的轉(zhuǎn)換數(shù)(si)大于顛換數(shù)(sv)。密碼子第二位點的轉(zhuǎn)換和顛換發(fā)生的頻率明顯高于第一、三位點。

    表2 Cytb密碼子不同位點的堿基頻率及堿基替換

    2.1.2 COI基因的堿基組成分析

    用Mega4.0軟件對E.altissima和Eupristinasp.小蜂(17個樣本)校準后的501bp序列進行堿基組成分析,結(jié)果表明:COI基因序列中,不變位點434個,變異位點60個,簡約信息位點55個,自裔位點4個,變異性為12.0%,與Cytb序列的變異性(14.0%)相似。

    A,T,G,C的平均含量分別為32.4%,43.0%,15.9%,8.7%,A+T的含量為75.5%,G+C含量為24.5%,A+T的含量明顯高于G+C的含量,符合膜翅目昆蟲線粒體基因序列的共同特點。

    COI密碼子不同位點的核苷酸頻率及堿基替換統(tǒng)計見表3,表中可知,與Cytb序列相似COI序列密碼子的各個位點上A+T含量均偏高,第一位點A+T的含量為85.3%,A+T偏向尤為明顯,其中T的使用最為頻繁為44.5%。說明COI基因在密碼子使用上也具有偏向性。E.altissima和Eupristinasp.小蜂COI基因序列上核苷酸顛換數(shù)(sv)大于轉(zhuǎn)換數(shù)(si)。密碼子第一位點的轉(zhuǎn)換和顛換發(fā)生的頻率明顯高于第二、三位點。

    表3 COI密碼子各堿基頻率及堿基替換

    2.2 Cytb及COI序列遺傳距離分析

    2.2.1 Cytb基因遺傳距離分析

    用Mega4.0軟件中對Eupristina屬小蜂Cytb序列計算K2P遺傳距離計算,結(jié)果表明:E.altissima10個樣本之間的遺傳距離為0—0.0287,平均遺傳距離為0.0092;Eupristinasp. 10個樣本之間的的遺傳距離為0—0.0055,平均遺傳距離為0.0030,表明同種個體之間遺傳距離差異小,而E.altissima與Eupristinasp.之間的遺傳距離為0.1527—0.1678,平均距離為0.1588,表明E.altissima與Eupristinasp.之間的遺傳距離明顯大于它們各自個體之間的遺傳距離。

    E.altissima雌蜂樣本間的遺傳距離為0,雄蜂樣本之間的遺傳距離為0.0000—0.0287,平均遺傳距離為0.0157;Eupristinasp.雌蜂樣本中除了旗山采集的樣本與別的樣本之間的遺傳距離為0.0055,其余各樣本之間的遺傳距離都為0;雄蜂樣本中除了倉山采集的樣本與別的樣本之間的遺傳距離為0.0055,其余樣本之間的距離都為0。兩種小蜂的雌蜂之間的遺傳距離在0.1582—0.1587之間,平均遺傳距離為0.1586,兩種雄蜂之間的遺傳距離為0.1527—0.1678,平均遺傳距離為0.1590,表明E.altissima與Eupristinasp.之間,雄蜂之間Cytb遺傳距離大于雌蜂。

    2.2.2 COI基因遺傳距離分析

    用Mega4.0軟件中對Eupristina屬小蜂COI序列計算K2P遺傳距離計算,結(jié)果表明:E.altissima7個樣本之間的遺傳距離為0—0.0165,平均遺傳距離為0.0065;Eupristinasp. 10個樣本之間的的遺傳距離為0—0.0379,平均遺傳距離為0.0205,表明種內(nèi)個體之間遺傳距離差異較小,而E.altissima與Eupristinasp.之間的遺傳距離為0.0940—0.1149,平均距離為0.1043,表明E.altissima與Eupristinasp.之間的遺傳距離遠遠大于它們各自個體之間遺傳距離。

    E.altissima雌蜂樣本間的遺傳距離為0,雄蜂樣本之間的遺傳距離為0.0020—0.0144,平均遺傳距離為0.0110;Eupristinasp.雌蜂樣本之間的遺傳距離為0—0.0102,平均遺傳距離為0.0057,雄蜂樣本之間的遺傳距離為0.0061—0.0252,平均遺傳距離為0.0139,表明種內(nèi)雄蜂個體之間的遺傳距離大于雌蜂。兩小蜂的雌蜂之間的遺傳距離在0.1015—0.1042之間,平均為0.1036,兩雄蜂之間的遺傳距離為0.0940—0.1149,平均遺傳距離為0.1052,表明E.altissima與Eupristinasp. 雄蜂之間的COI序列遺傳距離大于雌蜂。

    2.3 系統(tǒng)發(fā)育樹的構(gòu)建

    從GenBank上下載Eupristina屬和Blastophaga屬(作為外群)榕小蜂的Cytb和COI序列(表4),與實驗樣本一起構(gòu)建系統(tǒng)發(fā)育樹(圖1,圖2)。

    表4 GeneBank中下載的相關(guān)榕小蜂序列

    用Mega4.0軟件對實驗所得與GenBank中下載的榕小蜂Cytb與COI序列用NJ(Neighbor-joining)法,同時用重復抽樣1000次檢驗分子系統(tǒng)樹各分支的置信值,建立系統(tǒng)發(fā)育樹見圖1、圖2。圖中顯示,E.altissima的各樣本與Eupristinasp.各樣本各自匯聚成一支,顯示出較近的親緣關(guān)系。基于兩種分子標記序列建立的系統(tǒng)發(fā)育樹顯示出相似的拓撲結(jié)構(gòu),E.altissima與Eupristinasp.榕小蜂個體先分別匯聚成一支然后兩群體再匯聚到一起與該屬的別的種類的小蜂匯聚,最后與外群Blastophaga屬小蜂匯聚,這顯示出兩群體親緣關(guān)系較近的同時兩群體之間產(chǎn)生了變異。兩個屬及屬內(nèi)的各種小蜂都能夠明顯的區(qū)分開來,這說明Cytb與COI基因可以作為榕小蜂鑒別的分子標記。

    3 討論

    3.1 E. altissima和Eupristina sp.小蜂之間的分子遺傳關(guān)系

    本研究用線粒體Cytb及COI序列對福州地區(qū)高榕榕果內(nèi)存在的E.altissima與Eupristinasp.小蜂進行分析,Cytb序列分析顯示:E.altissima和Eupristinasp. 群體內(nèi)各樣本間的平均遺傳距離分別為0.0092和0.0030;而E.altissima與Eupristinasp. 群體間的平均距離為0.1588,表明兩者群體間的遺傳距離明顯大于各自群體內(nèi)各樣本間的遺傳距離。統(tǒng)計GenBank中下載的6個屬34種榕小蜂Cytb序列的種間遺傳距離為0.0811—0.1723,因此,Cytb序列分析表明E.altissima與Eupristinasp.為兩個不同的種。COI序列分析顯示:E.altissima7個樣本之間的平均遺傳距離為0.0065,Eupristinasp. 10個樣本之間的平均遺傳距離為0.0205;而二者之間的平均遺傳距離為0.1043,表明二者之間的遺傳距離明顯大于它們各自個體之間遺傳距離。統(tǒng)計GenBank中下載的6個屬28種榕小蜂COI序列的種間遺傳距離為0.0939—0.1986,表明E.altissima和Eupristinasp.之間的遺傳距離差異已經(jīng)達到了種間水平,因此,分子標記的實驗結(jié)果支持形態(tài)標記的鑒定結(jié)果,即E.altissima與Eupristinasp.為兩個不同的種。由于同物種的不同個體間差異小,不存在地理差異,而同屬不同種間遺傳距離的差異明顯,因此可以將本實驗序列作為分析及鑒別E.altissima與Eupristinasp.的依據(jù)。

    3.2 雌、雄蜂進化速度分析

    在形態(tài)上,E.altissima與Eupristinasp.的雌蜂之間有微小差異,而二者雄蜂之間無差異,但Cytb與COI序列分析結(jié)果一致表明:E.altissima與Eupristinasp.雄蜂之間的遺傳距離與二者雌蜂之間的遺傳距離接近,甚至略大于雌蜂。因此,僅比較Cytb與COI基因序列,雌、雄的進化速度是基本一致的,同時,也表明形態(tài)變異滯后于基因變異。雌蜂在表型上進化快于雄蜂,可能是由于雌蜂羽化后從榕果出飛,尋找雌花期榕果產(chǎn)卵,雌蜂在果外生活期間受到外界溫度,濕度,光照等環(huán)境因素的影響較大,尤其在傳粉行為上的差異,導致E.altissima與Eupristinasp.雌蜂之間的形態(tài)發(fā)生變異,而雄蜂終生生活在黑暗封閉、溫度和濕度都相對恒定的榕果內(nèi),且雄蜂羽化后存活壽命短。E.altissima與Eupristinasp.雄蜂之間生存周期、生活環(huán)境和活動行為的高度一致性,導致其表型變異極為緩慢。

    圖1 鄰接法構(gòu)建的Eupristina屬榕小蜂Cytb基因系統(tǒng)發(fā)育樹 Fig.1 NJ dendrograms based on the Cyt b gene sequence of Eupristina fig wasps

    圖2 鄰接法構(gòu)建的Eupristina屬榕小蜂COI基因系統(tǒng)發(fā)育樹 Fig.2 NJ dendrograms based on the COI gene sequence of Eupristina fig wasps

    3.3 兩種小蜂可能的進化模式

    段注標等對云南西雙版納地區(qū)高榕榕果中小蜂種群全年動態(tài)追蹤觀察結(jié)果表明,在4—7月間高榕的榕果只有Eupristinasp.,未見E.altissima,推測E.altissima小蜂可能發(fā)生了周期性的宿主轉(zhuǎn)移的現(xiàn)象[25]。本課題組從2005年到2012年對福州地區(qū)高榕榕果中小蜂種群全年動態(tài)追蹤觀察結(jié)果發(fā)現(xiàn),每年3—7月份生長在福州地區(qū)的高榕榕果內(nèi)沒有E.altissima和Eupristinasp.小蜂,而在8月份兩種小蜂又重新出現(xiàn),還觀察到在實驗室內(nèi)E.altissima小蜂快速進入聚果榕雌花期花序中(本實驗室未發(fā)表數(shù)據(jù)),上述觀察結(jié)果表明,兩種小蜂對宿主的專一性不強,可能易發(fā)生宿主轉(zhuǎn)移現(xiàn)象。由于E.altissima和Eupristinasp.之間的形態(tài)極為相似,推測這兩種小蜂早期可能是同一種蜂,在宿主轉(zhuǎn)移過程中發(fā)生種群分化。可能性一,Eupristinasp.原本是一種傳粉小蜂,其原宿主可能并不是高榕,由于原宿主滅絕,Eupristinasp.通過宿主轉(zhuǎn)移成功定居于高榕榕果中,高榕專一性傳粉小蜂E.altissima的存在保障了高榕的繁殖,Eupristinasp.逐漸喪失了傳粉功能,盡管Eupristinasp.不傳粉,也沒有因選擇作用而被排斥,兩種小蜂由于生活在同一個宿主內(nèi),在形態(tài)上發(fā)生了趨同進化,因而二者形態(tài)相似。另一種可能是,高榕的傳粉小蜂E.altissima通過宿主轉(zhuǎn)移定居于新宿主榕果中,由于新宿主原有的專一性傳粉小蜂行傳粉功能,保障了榕樹的繁殖,使E.altissima小蜂在新宿主榕果內(nèi)逐漸喪失了傳粉功能,而逐漸演化成為“Eupristinasp.”。Eupristinasp.可能既可寄生于新宿主榕樹,也可寄生于原宿主高榕,其雌蜂因未行傳粉作用,在形態(tài)上與傳粉相關(guān)的結(jié)構(gòu)退化或消失,而雄蜂的功能沒有改變,故形態(tài)尚未變異。

    致謝:感謝福建師范大學南方生物醫(yī)學研究中心陳騏教授對本文寫作的幫助。

    [1] Hill D S. Figs (Ficusspp.) of Hong Kong. Hong Kong: Hong Kong University Press, 1967.

    [2] Zhang X S, Wu Z Y, Cao Z Y. Flora of China. Vol. 23, No. 1. Beijing: Science Press, 1998: 95- 97.

    [3] Nair P B, Joseph M, Abdurahiman U. New fig wasps (Hymenoptera: Chalcidoidea) fromFicusaltissima. Entomology (Proceedings C), 1981, 84(2): 145- 153.

    [4] McKey D. Population biology of figs: applications for conservation. Experientia, 1989, 45(7): 661- 673.

    [5] Patel A, Hossaert-McKey M, McKey D. Ficus-pollinator research in India: past, present and future. Current Science, 1993, 65(3): 243- 243.

    [7] Gu H Y, Yang D R, Zhang G M, Peng Y Q, Song Q S. Species of fig wasps inFicusaltissimaand their ecological characters. Chinese Journal of Ecology, 2003, 22(2): 70- 73.

    [8] Peng Y Q, Duan Z B, Yang D R, Rasplus J Y. Co-occurrence of twoEupristinaspecies onFicusaltissimain Xishuangbanna, SW China. Symbiosis, 2008, 45(1/3): 9- 14.

    [9] Hebert P D N, Cywinska A, Ball S L, de Waard J R. Biological identifications through DNA barcodes. Proceedings of the Royal Society of London. Proceedings of the Royal Society B: Biological Sciences, 2003, 270(1512): 313- 321.

    [10] Pestano J, Brown R P, Suárez N M, Báez M. Diversification of sympatricSapromyza(Diptera: Lauxaniidae) from Madeira: six morphological species but only four mtDNA lineages. Molecular Phylogenetics and Evolution, 2003, 27(3): 422- 428.

    [11] Hebert P D, Penton E H, Burns J M, Janzen D H, Hallwachs W. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterflyAstraptesfulgerator. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(41): 14812- 14817.

    [12] Foster B T, Cognato A I, Gold R E. DNA-based identification of the eastern subterranean termite,Reticulitermesflavipes(Isoptera: Rhinotermitidae). Journal of Economic Entomology, 2004, 97(1): 95- 101.

    [13] Ward R D, Zemlak T S, Innes B H, Last P R, Hebert P D N. DNA barcoding Australia′s fish species. Philosophical Transactions of the Royal Society B: Biological Sciences, 2005, 360(1462): 1847- 1857.

    [14] Cao X R, Shu F Y, Zhang X R, Bi C M, Li C J, Hu J, Fang J Y. Phylogenetic relationships ofElaphodusCephalophusand threeMuntiacusspecies revealed by mitochondrial cytochrome B nucleotide sequence. Acta Zoologica Sinica, 2002, 48(1): 44- 49.

    [15] Xiao J H, Xiao H, Huang D W. DNA barcoding: new approach of biological taxonomy. Acta Zoologica Sinica, 2004, 50(5): 852- 855.

    [16] Herre E A, Machado C A, Bermingham E, Nason J D, Windsor D M, McCafferty S S, van Houten W, Bachmann K. Molecular phylogenies of figs and their pollinator wasps. Journal of Biogeography, 1996, 23(4): 521- 530.

    [17] Su Z H, Iino H, Nakamura K, Serrato A, Oyama K. Breakdown of the one-to-one rule in Mexican fig-wasp associations inferred by molecular phylogenetic analysis. Symbiosis, 2008, 45(1/3): 73- 81.

    [18] Molbo D, Machado C A, Sevenster J G, Keller L, Herre E A. Cryptic species of fig-pollinating wasps: Implications for the evolution of the fig-wasp mutualism, sex allocation, and precision of adaptation. Proceedings of the National Academy of Sciences of the United States of America, 2003, 100(10): 5867- 5872.

    [19] Xiao J H, Wang N X, Li Y W, Murphy R W, Wan D G, Niu L M, Hu H Y, Fu Y G, Huang D W. Molecular approaches to identify cryptic species and polymorphic species within a complex community of fig wasps. PLoS One, 2010, 5(11): 1- 7.

    [20] Zhou M J, Xiao J H, Bian S N, Li Y W, Niu L M, Hu H Y, Wu W S, Murphy R W, Huang D W. Molecular approaches identify known species, reveal cryptic species and verify host specificity of Chinese Philotrypesis (Hymenoptera: Pteromalidae). Molecular Ecology Resources, 2012, 12(4): 598- 606.

    [21] Li Y W, Zhou X, Feng G, Hu H Y, Niu L M, Hebert P D, Huang D W. COI and ITS2 sequences delimit species, reveal cryptic taxa and host specificity of fig-associatedSycophila(Hymenoptera, Eurytomidae). Molecular Ecology Resources, 2010, 10(1): 31- 40.

    [22] McLeish M J, Van Noort S. Codivergence and multiple host species use by fig wasp populations of theFicuspollination mutualism. BMC Evolutionary Biology, 2012, 12(1): 1- 12.

    [23] Jousselin E, van Noort S, Greeff J M. Labile male morphology and intraspecific male polymorphism in thePhilotrypesisfig wasps. Molecular Phylogenetics and Evolution, 2004, 33(3): 706- 718.

    [24] Wu W S, Chen Y L, Wu X G, Zhang Y J, Zheng C F. Chalcid wasps onFicusvirensAit. Var.sublanceolata(Miq.) Corner in Fuzhou. Journal of Environmental Entomology, 2012, 34(1): 57- 62.

    [25] Duan B Z, Peng Y Q, Yang D R, Xu L. Fig wasp community in the syconia ofFicusaltissima. Acta Ecologica Sinica, 2005, 25(10): 2589- 2594.

    參考文獻:

    [2] 張秀實, 吳征鎰, 曹子余. 中國植物志 (第二十三卷, 第一分冊). 北京: 科學出版社, 1998: 95- 97.

    [7] 谷海燕, 楊大榮, 張光明, 彭艷瓊, 宋啟示. 高榕隱頭果內(nèi)寄生蜂種類及生態(tài)學特征初步觀察. 生態(tài)學雜志, 2003, 22(2): 70- 73.

    [14] 曹祥榮, 束峰玨, 張錫然, 畢春明, 李朝軍, 胡均, 方箋陽. 毛冠鹿與3種麂屬動物的線粒體細胞色素b的系統(tǒng)進化分析. 動物學報, 2002, 48(1): 44- 49.

    [15] 肖金花, 肖暉, 黃大衛(wèi). 生物分類學的新動向——DNA條形編碼. 動物學報, 2004, 50(5): 852- 855.

    [24] 吳文珊, 陳友鈴, 伍小剛, 張彥杰, 鄭翠芳. 福州大葉榕隱頭果內(nèi)榕小蜂的分類. 環(huán)境昆蟲學報, 2012, 34(1): 57- 62.

    [25] 段柱標, 彭艷瓊, 楊大榮, 徐磊. 高榕隱頭果內(nèi)小蜂群落的動態(tài)變化. 生態(tài)學報, 2005, 25(10): 2589- 2594.

    ThegeneticevolutionaryrelationshipsoftwoEupristinaspeciesonFicusaltissima

    CHEN Youling, SUN Lingli, WU Leilei, WU Wenshan*, YAN Juyuan, LIU Yan

    ProvincialKeyLaboratoryforDevelopmentalBiologyandNeurosciences,CollegeofLifeScience,FujianNormalUniversity,Fuzhou350117,China

    We observed the habits of twoEupristinaspecies living inFicusaltissima. Only two fig wasp species,EupristinaaltissimaandEupristinasp., are found to be associated with the syconia ofFicusaltissimaaround the Fuzhou area, China. The former wasp is theFicusaltissimapollinator, while the latter wasp does not have pollination behavior. Some subtle differences are found in females of both wasp species regarding body color, antennae, pollen pockets, comb of setae on the fore coax. So far, no clear characters have been found to be used for distinguishing the males of these two wasp species. To better understand the genetic relationship betweenE.altissimaandEupristinasp., we cloned the mitochondrial Cytb and COI genes derived from bothE.altissimaandEupristinasp. collected from five different locations in Fuzhou.Eupristinaverticillata, the pollinator ofFicusmicrocarpa, was used as an outgroup control. We analyzed the base composition and genetic distance of the gene sequences and constructed a neighbor-joining phylogenetic tree to investigate the genetic evolutionary relationships. The data are summarized as follows: (1) The base composition of Cytb and COI show a high adenine (A) and thymine (T) bias (Cytb, A+T=75.3%; COI, A+T=75.5%), which pattern is consistent with the typical characteristics of mitochondrial gene sequences of hymenopteran. (2) Genetic distance analyses of two wasps groups show that the average distance of the Cytb sequence within each group ofE.altissimaandEupristinasp. are 0.0092 and 0.0030, respectively, far beyond the average genetic distance betweenE.altissimaandEupristinasp. that is 0.1588; COI sequence analyses show that the average genetic distances within each group ofE.altissimaandEupristinasp. are 0.0065 and 0.0205, respectively, although the average genetic distance betweenE.altissimaandEupristinasp. is 0.1043. Thus, the genetic distance between the two groups is significantly greater than that found within each respective group. Statistic analyses of the Cytb sequences derived from 6 genera and 34 species of wasps from GenBank show that the interspecies genetic distances of the Cytb sequence are between 0.0811 and 0.1723, and the interspecies genetic distances in the COI sequences derived from 6 genera and 28 species of wasps are between 0.0939 and 0.1986. These data suggest that the genetic distance betweenE.altissimaandEupristinasp. has reached to an interspecies difference level. Therefore, our molecular marker studies support the morphological characterization, namelyE.altissimaandEupristinasp. are indeed two different species. (3) Only a slight difference is found in females between the two species regarding the morphological characteristics, and no tangible difference is found in the males. However, the data of genetic distances derived from Cytb and COI sequence analyses consistently show the significant differences between the two species regardless their genders. It is likely that the morphological differentiation lags behind the genetic differentiation and the female wasps may evolve faster than the male wasps in the phenotypic features. The females fly out from figs after the emergence, which may be significantly affected by environmental factors resulting in the morphological differentiation between the female wasps of the two species. In contrast, the male wasps of the two species live in the closed darkness syconia with constantly environment lifelong, thus, no behavior differentiation may occur along with the extremely slow events of phenotypic mutating. The host specificity ofE.altissimaandEupristinasp. is not as stringent as previously assumed. The host switching events may have occurred during the fig-wasps co-evolution process.

    Ficusaltissima;Eupristinaaltissima; DNA barcoding; genetic structure; Cytb; COI

    國家自然科學基金資助項目(31270440);福建省科技廳重點項目(2011N0014);福建省高等學校學科帶頭人培養(yǎng)計劃資助項目;福建師范大學生命科學學院生物學拔尖生培養(yǎng)計劃項目

    2013- 05- 23;

    2013- 07- 25

    *通訊作者Corresponding author.E-mail: wuwenshan@126.com

    10.5846/stxb201305231153

    陳友鈴,孫伶俐,武蕾蕾,吳文珊,嚴菊媛,劉燕.高榕榕果內(nèi)Eupristina屬兩種榕小蜂的遺傳進化關(guān)系.生態(tài)學報,2013,33(19):6058- 6064.

    Chen Y L, Sun L L, Wu L L, Wu W S, Yan J Y, Liu Y.The genetic evolutionary relationships of twoEupristinaspecies onFicusaltissima.Acta Ecologica Sinica,2013,33(19):6058- 6064.

    猜你喜歡
    雌蜂雄蜂小蜂
    管氏腫腿蜂密度對其寄生及繁衍能力的影響
    孟氏胯姬小蜂種群動態(tài)及寄生特性
    楚雄腮扁葉蜂聚集產(chǎn)卵引誘物質(zhì)初探
    割雄蜂蛹經(jīng)驗談
    古橋槳角蚜小蜂與雙斑恩蚜小蜂的發(fā)育及形態(tài)比較
    桉樹枝癭姬小蜂風險分析及風險管理
    筆管榕榕小蜂的群落結(jié)構(gòu)與物種多樣性
    蠅蛹金小蜂的交配行為及雄蜂交配次數(shù)對雌蜂繁殖的影響
    割雄蜂房要做到事半功倍
    談雄蜂蛹的取與舍
    婷婷色av中文字幕| 国产欧美亚洲国产| 男女之事视频高清在线观看 | 亚洲天堂av无毛| 亚洲人成电影观看| 亚洲熟女精品中文字幕| 欧美日韩成人在线一区二区| 在线观看免费日韩欧美大片| 日韩av在线免费看完整版不卡| 欧美在线一区亚洲| 免费观看人在逋| 性色av乱码一区二区三区2| 天堂中文最新版在线下载| 亚洲国产欧美一区二区综合| 18在线观看网站| 99热国产这里只有精品6| 首页视频小说图片口味搜索 | 精品亚洲成a人片在线观看| 欧美黑人精品巨大| 黄色一级大片看看| 亚洲九九香蕉| 久久国产精品大桥未久av| 男女下面插进去视频免费观看| 免费高清在线观看视频在线观看| 欧美日韩福利视频一区二区| 狂野欧美激情性xxxx| 青春草亚洲视频在线观看| 亚洲专区中文字幕在线| 少妇猛男粗大的猛烈进出视频| 美女国产高潮福利片在线看| 久久精品国产a三级三级三级| 亚洲男人天堂网一区| 永久免费av网站大全| 极品人妻少妇av视频| 久久综合国产亚洲精品| 国产精品久久久久久精品电影小说| 亚洲国产精品国产精品| 捣出白浆h1v1| 国产成人精品在线电影| h视频一区二区三区| 1024香蕉在线观看| 亚洲欧美一区二区三区黑人| 久久精品人人爽人人爽视色| 男女下面插进去视频免费观看| 制服诱惑二区| 国产高清国产精品国产三级| 一本久久精品| 三上悠亚av全集在线观看| 久久国产精品男人的天堂亚洲| 亚洲精品一二三| 国产高清不卡午夜福利| 成人国产av品久久久| 一区在线观看完整版| 久久久久久久久免费视频了| 狠狠精品人妻久久久久久综合| 欧美日本中文国产一区发布| 日韩制服骚丝袜av| 曰老女人黄片| 久久av网站| 亚洲av国产av综合av卡| 精品少妇久久久久久888优播| 热re99久久精品国产66热6| 深夜精品福利| 国产麻豆69| 国产成人精品久久二区二区免费| 亚洲国产毛片av蜜桃av| 中文字幕人妻丝袜制服| 夫妻午夜视频| 性色av乱码一区二区三区2| 国产麻豆69| 天天躁日日躁夜夜躁夜夜| 国产亚洲av片在线观看秒播厂| 99精国产麻豆久久婷婷| 夫妻午夜视频| 亚洲国产欧美日韩在线播放| 成年人免费黄色播放视频| 欧美+亚洲+日韩+国产| 色综合欧美亚洲国产小说| 伦理电影免费视频| 国产精品熟女久久久久浪| 悠悠久久av| 日韩电影二区| 美女高潮到喷水免费观看| 黄频高清免费视频| 欧美日韩成人在线一区二区| 七月丁香在线播放| 欧美大码av| 欧美乱码精品一区二区三区| 午夜两性在线视频| 国产亚洲av片在线观看秒播厂| 免费在线观看完整版高清| 中文字幕人妻丝袜制服| 欧美人与性动交α欧美精品济南到| 中文精品一卡2卡3卡4更新| 麻豆国产av国片精品| 久久99热这里只频精品6学生| 欧美 日韩 精品 国产| 黄色一级大片看看| 亚洲国产欧美在线一区| 麻豆乱淫一区二区| 国产在线免费精品| 自拍欧美九色日韩亚洲蝌蚪91| 日韩 亚洲 欧美在线| 蜜桃在线观看..| 久久九九热精品免费| 免费在线观看日本一区| 18禁观看日本| 亚洲av日韩在线播放| 欧美中文综合在线视频| 欧美日韩亚洲高清精品| 成年美女黄网站色视频大全免费| 亚洲国产中文字幕在线视频| 久久久久久久大尺度免费视频| 九草在线视频观看| 一本色道久久久久久精品综合| 不卡av一区二区三区| 青春草视频在线免费观看| 亚洲一区二区三区欧美精品| 真人做人爱边吃奶动态| 男人添女人高潮全过程视频| 欧美日韩视频精品一区| 又黄又粗又硬又大视频| 亚洲成人国产一区在线观看 | 亚洲精品国产区一区二| 狠狠婷婷综合久久久久久88av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一卡2卡三卡4卡5卡 | 90打野战视频偷拍视频| 男女边吃奶边做爰视频| 丝袜在线中文字幕| 精品免费久久久久久久清纯 | 久久精品aⅴ一区二区三区四区| 婷婷色麻豆天堂久久| 黄色怎么调成土黄色| 欧美+亚洲+日韩+国产| 黄频高清免费视频| 亚洲精品久久午夜乱码| 免费黄频网站在线观看国产| 久久精品成人免费网站| 黄色毛片三级朝国网站| 精品人妻在线不人妻| 亚洲色图 男人天堂 中文字幕| 亚洲av片天天在线观看| 欧美在线黄色| 别揉我奶头~嗯~啊~动态视频 | 波多野结衣av一区二区av| 美女大奶头黄色视频| 两个人看的免费小视频| 日本91视频免费播放| 纵有疾风起免费观看全集完整版| 日日夜夜操网爽| 国产亚洲av片在线观看秒播厂| a级毛片在线看网站| 天堂8中文在线网| 男人舔女人的私密视频| 色精品久久人妻99蜜桃| xxxhd国产人妻xxx| avwww免费| 久久九九热精品免费| netflix在线观看网站| h视频一区二区三区| 午夜av观看不卡| 9191精品国产免费久久| 人人妻人人爽人人添夜夜欢视频| av又黄又爽大尺度在线免费看| 亚洲欧美精品综合一区二区三区| 国产成人一区二区在线| 91麻豆av在线| 欧美亚洲日本最大视频资源| 纵有疾风起免费观看全集完整版| 99久久综合免费| 亚洲国产看品久久| 午夜老司机福利片| 久久久久国产精品人妻一区二区| 老司机影院毛片| 国产男女内射视频| 好男人电影高清在线观看| 大型av网站在线播放| 黑人欧美特级aaaaaa片| 免费少妇av软件| 欧美精品一区二区免费开放| 男人爽女人下面视频在线观看| 亚洲成av片中文字幕在线观看| 少妇人妻久久综合中文| 啦啦啦啦在线视频资源| 人成视频在线观看免费观看| av国产久精品久网站免费入址| 久久亚洲精品不卡| 极品人妻少妇av视频| 97精品久久久久久久久久精品| 国产精品九九99| 高清av免费在线| 另类精品久久| 国产精品久久久久久人妻精品电影 | av片东京热男人的天堂| 曰老女人黄片| 99精国产麻豆久久婷婷| 亚洲精品国产色婷婷电影| 亚洲熟女毛片儿| 欧美人与善性xxx| 最近最新中文字幕大全免费视频 | 九草在线视频观看| 最新在线观看一区二区三区 | 欧美av亚洲av综合av国产av| 国产精品av久久久久免费| 欧美黄色淫秽网站| 中文字幕精品免费在线观看视频| 午夜福利视频在线观看免费| 久久女婷五月综合色啪小说| 亚洲天堂av无毛| 高潮久久久久久久久久久不卡| 一本大道久久a久久精品| 午夜免费成人在线视频| 五月开心婷婷网| 国产极品粉嫩免费观看在线| 脱女人内裤的视频| 欧美亚洲日本最大视频资源| 国产欧美日韩精品亚洲av| 日本av免费视频播放| 久久99热这里只频精品6学生| 丝袜在线中文字幕| 岛国毛片在线播放| 日韩伦理黄色片| 亚洲av欧美aⅴ国产| 国产成人啪精品午夜网站| 一边摸一边做爽爽视频免费| 性高湖久久久久久久久免费观看| 欧美久久黑人一区二区| 国产又爽黄色视频| 少妇人妻久久综合中文| 三上悠亚av全集在线观看| 母亲3免费完整高清在线观看| 黄色毛片三级朝国网站| 欧美人与善性xxx| 成人三级做爰电影| 日本91视频免费播放| 欧美日本中文国产一区发布| 欧美久久黑人一区二区| 国产午夜精品一二区理论片| av又黄又爽大尺度在线免费看| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站| av在线播放精品| 成人午夜精彩视频在线观看| 91麻豆av在线| 无遮挡黄片免费观看| 成人手机av| 色婷婷久久久亚洲欧美| 一级毛片黄色毛片免费观看视频| 丝袜美足系列| 日本黄色日本黄色录像| 久久ye,这里只有精品| 涩涩av久久男人的天堂| 日本wwww免费看| 精品少妇一区二区三区视频日本电影| 国产女主播在线喷水免费视频网站| 少妇猛男粗大的猛烈进出视频| 国产一卡二卡三卡精品| 亚洲国产中文字幕在线视频| 搡老岳熟女国产| 欧美日韩福利视频一区二区| 999精品在线视频| 亚洲国产欧美一区二区综合| 9191精品国产免费久久| 久久久国产精品麻豆| 久久影院123| 亚洲欧美色中文字幕在线| 婷婷色综合大香蕉| 亚洲精品在线美女| 亚洲人成电影观看| 亚洲专区国产一区二区| 男女下面插进去视频免费观看| 99久久99久久久精品蜜桃| 亚洲,欧美,日韩| 一边亲一边摸免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 桃花免费在线播放| 精品亚洲乱码少妇综合久久| bbb黄色大片| 另类精品久久| 久久99精品国语久久久| 久久精品亚洲熟妇少妇任你| 夫妻性生交免费视频一级片| 亚洲成人手机| 国产在线免费精品| 午夜免费鲁丝| 中文字幕人妻丝袜一区二区| svipshipincom国产片| 欧美成人午夜精品| 欧美xxⅹ黑人| 欧美日韩国产mv在线观看视频| 午夜福利视频在线观看免费| 99香蕉大伊视频| 日韩一区二区三区影片| 日韩av在线免费看完整版不卡| 日本欧美国产在线视频| 黄频高清免费视频| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 欧美日韩视频高清一区二区三区二| 久久久精品免费免费高清| 国产一区亚洲一区在线观看| 91成人精品电影| 国产精品一国产av| 韩国高清视频一区二区三区| 一区二区三区精品91| 国产片特级美女逼逼视频| 欧美日韩亚洲综合一区二区三区_| 国精品久久久久久国模美| 宅男免费午夜| 久久国产精品人妻蜜桃| 亚洲精品一二三| 国产精品欧美亚洲77777| 国产日韩欧美在线精品| 婷婷丁香在线五月| 午夜影院在线不卡| 各种免费的搞黄视频| 黑人巨大精品欧美一区二区蜜桃| 久久久精品国产亚洲av高清涩受| 国产深夜福利视频在线观看| 精品久久久久久电影网| 亚洲精品一区蜜桃| 亚洲伊人久久精品综合| 亚洲国产精品国产精品| 性少妇av在线| 亚洲少妇的诱惑av| 国产成人av教育| 久久久久久久大尺度免费视频| 久久精品亚洲av国产电影网| 建设人人有责人人尽责人人享有的| 美女大奶头黄色视频| www日本在线高清视频| 欧美日韩综合久久久久久| 精品高清国产在线一区| 午夜福利在线免费观看网站| 国产成人a∨麻豆精品| 捣出白浆h1v1| 一区二区三区四区激情视频| 成人免费观看视频高清| 丝袜美足系列| 美女脱内裤让男人舔精品视频| 真人做人爱边吃奶动态| 午夜福利在线免费观看网站| 女人爽到高潮嗷嗷叫在线视频| 成人亚洲精品一区在线观看| 亚洲国产欧美网| 国产亚洲欧美在线一区二区| 高清不卡的av网站| 国产三级黄色录像| 精品高清国产在线一区| 成人亚洲欧美一区二区av| 久久国产亚洲av麻豆专区| 一二三四社区在线视频社区8| 久久影院123| 国产一卡二卡三卡精品| 精品一区二区三区av网在线观看 | 一级黄片播放器| cao死你这个sao货| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 亚洲精品一卡2卡三卡4卡5卡 | 日本午夜av视频| 国产成人欧美| 亚洲人成电影免费在线| 亚洲av在线观看美女高潮| 国产国语露脸激情在线看| 黄色 视频免费看| 国产精品国产三级专区第一集| 国产精品一区二区精品视频观看| 国产成人影院久久av| 亚洲五月色婷婷综合| 国产熟女欧美一区二区| 欧美 亚洲 国产 日韩一| 叶爱在线成人免费视频播放| 久久精品亚洲熟妇少妇任你| 电影成人av| 黄色视频在线播放观看不卡| 亚洲人成电影观看| 丝袜在线中文字幕| 99热国产这里只有精品6| 纯流量卡能插随身wifi吗| 精品国产乱码久久久久久小说| 久久久久久久国产电影| 亚洲国产av影院在线观看| 久久精品国产a三级三级三级| 极品人妻少妇av视频| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 国产淫语在线视频| 亚洲精品国产色婷婷电影| 国产91精品成人一区二区三区 | 大片电影免费在线观看免费| 另类亚洲欧美激情| 夫妻午夜视频| 热99久久久久精品小说推荐| 亚洲欧美日韩高清在线视频 | 欧美精品亚洲一区二区| 亚洲,欧美,日韩| 男女国产视频网站| 日韩精品免费视频一区二区三区| 国产麻豆69| 国产高清视频在线播放一区 | 国产精品一国产av| 男人爽女人下面视频在线观看| 成人亚洲欧美一区二区av| 校园人妻丝袜中文字幕| 亚洲成人国产一区在线观看 | 丝袜喷水一区| 色婷婷av一区二区三区视频| 成人国产一区最新在线观看 | 亚洲图色成人| 男人操女人黄网站| 热99国产精品久久久久久7| 老汉色av国产亚洲站长工具| 欧美国产精品一级二级三级| 国产亚洲精品第一综合不卡| 菩萨蛮人人尽说江南好唐韦庄| 丰满饥渴人妻一区二区三| 亚洲欧美一区二区三区黑人| 男男h啪啪无遮挡| 女人久久www免费人成看片| 国产免费一区二区三区四区乱码| 男女免费视频国产| 一边摸一边做爽爽视频免费| 亚洲人成电影免费在线| 亚洲精品国产区一区二| 午夜视频精品福利| 国产成人精品久久二区二区91| 亚洲国产最新在线播放| 久久久国产精品麻豆| 校园人妻丝袜中文字幕| 久久亚洲国产成人精品v| 亚洲综合色网址| 男人舔女人的私密视频| 视频在线观看一区二区三区| 超碰成人久久| 波野结衣二区三区在线| 久久久久久久精品精品| 亚洲人成网站在线观看播放| 久久久久国产一级毛片高清牌| 欧美+亚洲+日韩+国产| 香蕉国产在线看| 精品少妇一区二区三区视频日本电影| 亚洲国产精品一区二区三区在线| 国产成人系列免费观看| 91麻豆av在线| 午夜福利,免费看| 久久国产亚洲av麻豆专区| 国产熟女欧美一区二区| 一级a爱视频在线免费观看| 一区二区三区精品91| 三上悠亚av全集在线观看| 最新在线观看一区二区三区 | 亚洲精品自拍成人| 国产高清视频在线播放一区 | 亚洲自偷自拍图片 自拍| 男女下面插进去视频免费观看| 欧美成人午夜精品| 国产成人av教育| 免费少妇av软件| av天堂久久9| 女人久久www免费人成看片| 高清欧美精品videossex| 男的添女的下面高潮视频| 天堂俺去俺来也www色官网| 精品国产乱码久久久久久小说| 午夜免费成人在线视频| 国产精品欧美亚洲77777| 十八禁网站网址无遮挡| 天天躁日日躁夜夜躁夜夜| 亚洲av欧美aⅴ国产| 蜜桃国产av成人99| 18禁裸乳无遮挡动漫免费视频| 日韩欧美一区视频在线观看| 黄色a级毛片大全视频| 国产高清不卡午夜福利| 一级a爱视频在线免费观看| 欧美人与性动交α欧美精品济南到| videos熟女内射| 日本猛色少妇xxxxx猛交久久| 欧美人与善性xxx| 亚洲国产最新在线播放| 天天躁夜夜躁狠狠久久av| 黄色 视频免费看| 成人国产av品久久久| 亚洲欧美激情在线| 久久久精品国产亚洲av高清涩受| 一区福利在线观看| 久久久久久久精品精品| 国产伦人伦偷精品视频| 久久精品久久久久久久性| 日韩一卡2卡3卡4卡2021年| 精品一品国产午夜福利视频| av在线播放精品| 精品免费久久久久久久清纯 | 人人妻人人澡人人看| av在线老鸭窝| 久久国产精品人妻蜜桃| 最近中文字幕2019免费版| 一级a爱视频在线免费观看| 69精品国产乱码久久久| 午夜影院在线不卡| 精品一区二区三区av网在线观看 | 乱人伦中国视频| 国产男女内射视频| 欧美激情极品国产一区二区三区| 日本欧美国产在线视频| 亚洲美女黄色视频免费看| 丁香六月天网| 欧美久久黑人一区二区| 你懂的网址亚洲精品在线观看| 国产日韩欧美视频二区| e午夜精品久久久久久久| 国产午夜精品一二区理论片| 亚洲成人手机| 视频在线观看一区二区三区| 久久精品国产a三级三级三级| √禁漫天堂资源中文www| 啦啦啦 在线观看视频| 五月天丁香电影| 一区在线观看完整版| 午夜影院在线不卡| 日本av手机在线免费观看| 男女之事视频高清在线观看 | 色精品久久人妻99蜜桃| 国产深夜福利视频在线观看| 激情五月婷婷亚洲| 国产高清不卡午夜福利| 亚洲三区欧美一区| 极品人妻少妇av视频| 国产欧美日韩一区二区三 | 亚洲欧美一区二区三区国产| 人人妻人人澡人人看| 人妻 亚洲 视频| 免费一级毛片在线播放高清视频 | 国产亚洲精品久久久久5区| 丰满饥渴人妻一区二区三| 首页视频小说图片口味搜索 | 亚洲欧美一区二区三区久久| 亚洲人成电影观看| av在线老鸭窝| 一级毛片女人18水好多 | 妹子高潮喷水视频| 各种免费的搞黄视频| 精品少妇黑人巨大在线播放| 成年人午夜在线观看视频| 一边亲一边摸免费视频| 在线天堂中文资源库| 丝袜脚勾引网站| bbb黄色大片| 国产精品久久久人人做人人爽| 午夜福利视频精品| 国产黄色免费在线视频| 国产日韩欧美亚洲二区| 久久影院123| 狠狠婷婷综合久久久久久88av| 婷婷成人精品国产| 又粗又硬又长又爽又黄的视频| 黄色一级大片看看| 久热爱精品视频在线9| 久久这里只有精品19| 91麻豆精品激情在线观看国产 | 午夜影院在线不卡| av在线播放精品| 天天躁夜夜躁狠狠久久av| 人人妻,人人澡人人爽秒播 | 亚洲av电影在线进入| 亚洲精品一卡2卡三卡4卡5卡 | 欧美精品av麻豆av| 黑丝袜美女国产一区| 一边摸一边做爽爽视频免费| 一本色道久久久久久精品综合| 国产欧美日韩精品亚洲av| 9热在线视频观看99| 国产精品秋霞免费鲁丝片| 欧美老熟妇乱子伦牲交| 99国产精品一区二区三区| 中文字幕人妻熟女乱码| 两个人免费观看高清视频| 国产成人精品久久二区二区免费| 亚洲人成电影观看| 免费少妇av软件| 91精品国产国语对白视频| 日韩人妻精品一区2区三区| 女人高潮潮喷娇喘18禁视频| 精品第一国产精品| 国产福利在线免费观看视频| 国产深夜福利视频在线观看| 欧美久久黑人一区二区| 永久免费av网站大全| 韩国精品一区二区三区| 精品第一国产精品| 国产福利在线免费观看视频| 黑人巨大精品欧美一区二区蜜桃| 两个人看的免费小视频| cao死你这个sao货| 十八禁高潮呻吟视频| 黄色怎么调成土黄色| 久久久亚洲精品成人影院| 日本五十路高清| 亚洲欧美日韩高清在线视频 | 亚洲国产精品一区二区三区在线| 亚洲 国产 在线| 国产人伦9x9x在线观看| 日本a在线网址| 丝瓜视频免费看黄片| 国产精品一区二区免费欧美 | www.熟女人妻精品国产| 少妇裸体淫交视频免费看高清 | h视频一区二区三区| 一级毛片我不卡| 日本午夜av视频| 久久久欧美国产精品| 免费日韩欧美在线观看| 亚洲国产av新网站|